Quantum-instanton evaluation of the kinetic isotope effects

A general quantum-mechanical method for computing kinetic isotope effects is presented. The method is based on the quantum-instanton approximation for the rate constant and on the path-integral Metropolis-Monte Carlo evaluation of the Boltzmann operator matrix elements. It computes the kinetic isotope effect directly, using a thermodynamic integration with respect to the mass of the isotope, thus avoiding the more computationally expensive process of computing the individual rate constants. The method should be more accurate than variational transition-state theories or the semiclassical instanton method since it does not assume a single tunneling path and does not use a semiclassical approximation of the Boltzmann operator. While the general Monte Carlo implementation makes the method accessible to systems with a large number of atoms, we present numerical results for the Eckart barrier and for the collinear and full three-dimensional isotope variants of the hydrogen exchange reaction H + H2 --> H2 + H. In all seven test cases, for temperatures between 250 and 600 K, the error of the quantum instanton approximation for the kinetic isotope effects is less than approximately 10%.     

Two-Term and multi-term approximation of the nonstationary electron velocity distribution in an electric field in a gas

A very efficient, strict nonstationary multi-term approach has been developed as a generalization of the conventional and the strict nonstationary two-term approximations used for solving the nonstationary , electron Boltzmann equation. As a first application the temporal relaxation of electrons in a model plasma acted upon by a de electric field has been investigated. The results are compared with corresponding ones obtained by the conventional and the strict nonstationary two-term approach as well as with very accurate Monte Carlo simulations. Perfect agreement between nonstationary eight-term Boltzmann and Monte Carlo calculations is found.     

A ground state potential energy surface for H2 using Monte Carlo methods

Using variational Monte Carlo and a simple explicitly correlated wave function we have computed the Born-Oppenheimer energy of the H2 ground state (X 1Sigmag+) at 24 internuclear distances. We have also calculated the diagonal correction to the Born-Oppenheimer approximation and the lowest-order relativistic corrections at each distance using variational Monte Carlo techniques. The nonadiabatic values are evaluated from numerical derivatives of the wave function with respect to the nuclear coordinates. With this potential energy surface we have computed several of the lowest vibrational-rotational energies for this system. Our results are in good agreement with the best values found in the literature.     

A simple approach to treat anisotropic elastic collisions in Monte Carlo calculations of the electron energy distribution function in cold plasmas

A simplified method allowing one to treat anisotropic electron heavy species elastic scattering in MonteCarlo models of gas discharges with the proper value for collision frequency is proposed The method is applied to an electric discharge in a Ne · Xe/HCl mixture, and the results are compared with the solution of the two-term expansion of the Boltzmann equation under the same conditions. Methods for reduction of computational time in Monte Carlo codes and the use of the Monte Carlo flux method are also discussed.     

Free energies and phase equilibria of chain molecules

In this Article, a review is presented of recent developments in Monte Carlo simulations of chain molecules. The Rosenbluth chain insertion technique is used to calculate the free energy of the chain molecules. Furthermore, this insertion method is used to generate biased Monte Carlo moves. It is shown that this bias can be removed by adjusting the acceptance rules such that configurations are generated with their correct Boltzmann weight. This configurational-bias Monte Carlo method can be combined with the Gibbs-ensemble technique which results in an efficient method to simulate phase equilibria of chain molecules.     

The generalized van der Waals partition function. III. Local composition models for a mixture of equal size square-well molecules

New local composition models for mixtures of equal size molecules with differing attractive forces are presented, and compared with our results of Monte Carlo computer simulations for mixtures of square-well molecules which are also reported here. Unlike most previous local composition models, these new models predict random mixing in the high density limit and Boltzmann factor nonrandomness at low density; both limiting behaviors are in agreement with statistical mechanical theory. The predictions of these new models as a function of composition, density and temperature are in good agreement with the Monte Carlo computer simulation results.     

On a class of extended mean spherical approximations

The extended mean spherical approximation of Martina B, and del Rio is solved variationally for the restricted primitive model of 1—1 and 2—2 electrolyte solutions. Thermodynamic properties are calculated and compared to Monte Carlo results for the same model. A modification of the extended mean spherical approximation is presented that is more internally self-consistent and is in better agreement with Monte Carlo data.     

O—酰基—α—酮肟光分解反应的Monte Carlo处理

本文首次对光化学反应体系用Monte Carlo方法进行模拟处理。通过5个O-酰基-α-酮肟光分解反应的Monte Carlo模拟,可避免解析解中由于对吸收光强须采用一级近似求解动力学微分方程组,而造成拟合反应在后期产生与实验结果的偏差。     

Monte Carlo simulation of liquid Br2: Test of an effective hamiltoniain

Monte Carlo simulations for the structure of a fluid of molecules interacting with a pair potential having one Lennard- Jones center plus a quadrupole interaction is reported. The numerical experiments are compared with an approximation method previously suggested and with experimental results for the liquid Br₂ structure. This effective pair potential turns out to be unable to reproduce the Br₂ structure while the approximation method is only useful for quadrupole moments less than 1.2 × 10^?40 C m².     

Quantum Monte Carlo calculations of the potential energy curve of the helium dimer

We report results of two quantum Monte Carlo methods -- variational Monte Carlo and diffusion Monte Carlo -- on the potential energy curve of the helium dimer. In contrast to previous quantum Monte Carlo calculations on this system, we have employed trial wave functions of the Slater-Jastrow form and used the fixed node approximation for the fermion nodal surface. We find both methods to be in excellent agreement with the best theoretical results at short range. In addition, the diffusion Monte Carlo results give very good agreement across the whole potential energy curve, while the Slater-Jastrow wave function fails to bind the dimer at all.     

自由基聚合反应的Monte Carlo模拟方法

本文将化学反应动力学的MonteCarlo模拟方法运用到引发剂引发的自由基聚合反应的非稳态动力学,针对自由基聚合反应动力学数值模拟所特有的"无伸缩问题",采用"偏倚抽样法"解决了MonteCarlo模拟中的"无伸缩问题",模拟结果与非稳态动力学解的结果完全一致,此算法易推广到研究更复杂的自由基聚合反应体系。     

DFT-LDA pseudopotentials in quantum Monte Carlo

We investigate the portability of standard norm-conserving pseudopotentials outside the density functional theory-local density approximation (DFT-LDA) framework, i.e., their use and interpretation as electron-ion effective potentials in valence-only diffusion Monte Carlo simulations. While first-principles many-body pseudopotentials are not available in the literature yet, the use of approximate pseudopotentials in quantum Monte Carlo simulations is becoming widespread. Here we attempt a systematic analysis of the portability of norm-conserving pseudopotentials generated within DFT-LDA, focusing on a model many-body system, the two-electron valence-only ion. Our results indicate that the portability is good in most cases, hence the use of pseudopotentials in quantum Monte Carlo simulations is in general a reasonable approximation but suggest that in some cases this approximation may be relevant. © 1997 John Wiley & Sons, Inc.     

Flexible polyelectrolyte simulations at the Poisson-Boltzmann level: a comparison of the kink-jump and multigrid configurational-bias Monte Carlo methods

We present a new approach for simulating the motions of flexible polyelectrolyte chains based on the continuous kink-jump Monte Carlo technique coupled to a lattice field theory based calculation of the Poisson-Boltzmann (PB) electrostatic free energy "on the fly." This approach is compared to the configurational-bias Monte Carlo technique, in which the chains are grown on a lattice and the PB equation is solved for each configuration with a linear scaling multigrid method to obtain the many-body free energy. The two approaches are used to calculate end-to-end distances of charged polymer chains in solutions with varying ionic strengths and give similar numerical results. The configurational-bias Monte Carlo/multigrid PB method is found to be more efficient, while the kink-jump Monte Carlo method shows potential utility for simulating nonequilibrium polyelectrolyte dynamics.     

A version of diffusion Monte Carlo method based on random grids of coherent states. II. Six-dimensional simulation of electronic states of H2

We report a new version of the diffusion Monte Carlo (DMC) method, based on coherent-state quantum mechanics. Randomly selected grids of coherent states in phase space are used to obtain numerical imaginary time solutions of the Schrodinger equation, with an iterative refinement technique to improve the quality of the Monte Carlo grid. Accurate results were obtained, for the appropriately symmetrized two lowest states of the hydrogen molecule, by Monte Carlo sampling and six-dimensional propagation in the full phase space.     

Monte Carlo global sensitivity analysis of Boltzmann equation for electron kinetics in H2 discharges

Global sensitivity analysis with the Monte Carlo method is applied to the Boltzmann equation for the electron energy distribution function (eedf). The results show the sensitivity of eedf and related quantities on the global variation of cross sections set. A new indicator of global sensitivity is used, which appears to be of easy applicability.     

Quantum Monte Carlo study of helium clusters doped with nitrous oxide: quantum solvation and rotational dynamics

Dynamical and structural properties of small (4)He(N)-N(2)O complexes have been analyzed using ground-state and finite-temperature Monte Carlo simulations. The effective rotational constants resulting from the ground-state calculations are in excellent agreement with the results of a recent spectroscopic study [Y. Xu et al., Phys. Rev. Lett. 91, 163401 (2003)]. After an initial decrease for cluster sizes up to N=8, the rotational constant increases, signaling a transition from a molecular complex to a quantum-solvated system. Such a turnaround is not present in the rotational constants extracted from the finite-temperature Monte Carlo calculations, performed for Boltzmann statistics, thus highlighting the importance of exchange effects to explain the decoupling between a solvated dopant and the helium motion.     

Comparison of ISO-GUM and Monte Carlo methods for the evaluation of measurement uncertainty: application to direct cadmium measurement in water by GFAAS

The propagation stage of uncertainty evaluation, known as the propagation of distributions, is in most cases approached by the GUM (Guide to the Expression of Uncertainty in Measurement) uncertainty framework which is based on the law of propagation of uncertainty assigned to various input quantities and the characterization of the measurand (output quantity) by a Gaussian or a t-distribution. Recently, a Supplement to the ISO-GUM was prepared by the JCGM (Joint Committee for Guides in Metrology). This Guide gives guidance on propagating probability distributions assigned to various input quantities through a numerical simulation (Monte Carlo Method) and determining a probability distribution for the measurand.In the present work the two approaches were used to estimate the uncertainty of the direct determination of cadmium in water by graphite furnace atomic absorption spectrometry (GFAAS). The expanded uncertainty results (at 95% confidence levels) obtained with the GUM Uncertainty Framework and the Monte Carlo Method at the concentration level of 3.01 μg/L were ±0.20 μg/L and ±0.18 μg/L, respectively. Thus, the GUM Uncertainty Framework slightly overestimates the overall uncertainty by 10%. Even after taking into account additional sources of uncertainty that the GUM Uncertainty Framework considers as negligible, the Monte Carlo gives again the same uncertainty result (±0.18 μg/L). The main source of this difference is the approximation used by the GUM Uncertainty Framework in estimating the standard uncertainty of the calibration curve produced by least squares regression. Although the GUM Uncertainty Framework proves to be adequate in this particular case, generally the Monte Carlo Method has features that avoid the assumptions and the limitations of the GUM Uncertainty Framework.     

Depletion potentials in colloidal mixtures of hard spheres and rods

The depletion potential between a hard sphere and a planar hard wall, or two hard spheres, imposed by suspended rigid spherocylindrical rods is computed by the acceptance ratio method through the application of Monte Carlo simulation. The accurate results and ideal-gas approximation results of the depletion potential are determined with the acceptance ratio method in our simulations. For comparison, the depletion potentials are also studied by using both the density functional theory and Derjaguin approximations. The density profile as a function of positions and orientations of rods, used in the density functional theory, is calculated by Monte Carlo simulation. The potential obtained by the acceptance ratio method is in good agreement with that of density functional theory under the ideal-gas approximation. The comparison between our results and those of other theories suggests that the acceptance ratio method is the only efficient method used to compute the depletion potential induced by nonspherical colloids with the volume fraction beyond the ideal-gas approximation.     

Fixed—node Quantum Monte Carlo：A Novel Approach

The fixed-node quantum Monte Carlo (FNQMC)1,2 method has made it possible to calculate the electronic structure of relatively large molecular systems. These large systems range from positron complexes [NH2, Ps] with ～10 electrons to C20 isomers with 120 electrons, silicon crystal structures of 250 atoms and 1000 valence electrons. In the practical calculation for FNQMC method, in general, a minimal basis set of Slater-type atomic orbital (STO) and Jastrow functions are taken to cons…     

Higher order and infinite Trotter-number extrapolations in path integral Monte Carlo

Improvements beyond the primitive approximation in the path integral Monte Carlo method are explored both in a model problem and in real systems. Two different strategies are studied: The Richardson extrapolation on top of the path integral Monte Carlo data and the Takahashi-Imada action. The Richardson extrapolation, mainly combined with the primitive action, always reduces the number-of-beads dependence, helps in determining the approach to the dominant power law behavior, and all without additional computational cost. The Takahashi-Imada action has been tested in two hard-core interacting quantum liquids at low temperature. The results obtained show that the fourth-order behavior near the asymptote is conserved, and that the use of this improved action reduces the computing time with respect to the primitive approximation.