Monte Carlo evaluation of the equilibrium isotope effects using the Takahashi–Imada factorization of the Feynman path integral

Abstract The Feynman path integral approach for computing equilibrium isotope effects and isotope fractionation corrects the approximations made in standard methods, although at significantly increased computational cost. We describe an accelerated path integral approach based on three ingredients: the fourth-order Takahashi–Imada factorization of the path integral, thermodynamic integration with respect to mass, and centroid virial estimators for relevant free energy derivatives. While the first ingredient speeds up convergence to the quantum limit, the second and third improve statistical convergence. The combined method is applied to compute the equilibrium constants for isotope exchange reactions H 2 + D ⇌ H + HD and H 2 + D 2 ⇌ 2 HD .

[1]  M. R. Peterson,et al.  A refined H3 potential energy surface , 1996 .

[2]  R. Feynman,et al.  Quantum Mechanics and Path Integrals , 1965 .

[3]  Michele Ceotto,et al.  Quantum instanton approximation for thermal rate constants of chemical reactions , 2003 .

[4]  Masuo Suzuki,et al.  Hybrid exponential product formulas for unbounded operators with possible applications to Monte Carlo simulations , 1995 .

[5]  L. Rebelo,et al.  Isotope Effects: in the Chemical, Geological, and Bio Sciences , 2009 .

[6]  M. Parrinello,et al.  Study of an F center in molten KCl , 1984 .

[7]  D. Major,et al.  Path-Integral Calculations of Nuclear Quantum Effects in Model Systems, Small Molecules, and Enzymes via Gradient-Based Forward Corrector Algorithms. , 2011, Journal of chemical theory and computation.

[8]  Berend Smit,et al.  Understanding Molecular Simulation , 2001 .

[9]  A 'path-by-path' monotone extrapolation sequence for Feynman path integral calculations of quantum mechanical free energies , 2003 .

[10]  William H Miller,et al.  Path integral evaluation of the quantum instanton rate constant for proton transfer in a polar solvent. , 2005, The Journal of chemical physics.

[11]  Takeshi Yamamoto,et al.  Path-integral virial estimator based on the scaling of fluctuation coordinates: application to quantum clusters with fourth-order propagators. , 2005, The Journal of chemical physics.

[12]  V. Vesovic,et al.  Ab initio intermolecular potential energy surface and thermophysical properties of hydrogen sulfide. , 2011, Physical chemistry chemical physics : PCCP.

[13]  M. Wolfsberg.  Theoretical evaluation of experimentally observed isotope effects , 1972 .

[14]  Klein,et al.  Staging: A sampling technique for the Monte Carlo evaluation of path integrals. , 1985, Physical review. B, Condensed matter.

[15]  D. Truhlar,et al.  High-precision quantum thermochemistry on nonquasiharmonic potentials: converged path-integral free energies and a systematically convergent family of generalized Pitzer-Gwinn approximations. , 2005, The journal of physical chemistry. A.

[16]  Hagen Kleinert,et al.  Path Integrals in Quantum Mechanics , 1995 .

[17]  Jirí Vanícek,et al.  Efficient estimators for quantum instanton evaluation of the kinetic isotope effects: application to the intramolecular hydrogen transfer in pentadiene. , 2007, The Journal of chemical physics.

[18]  B. Berne,et al.  On path integral Monte Carlo simulations , 1982 .

[19]  R. Coalson.  On the connection between Fourier coefficient and Discretized Cartesian path integration , 1986 .

[20]  O. A. von Lilienfeld,et al.  Path Integral Computation of Quantum Free Energy Differences Due to Alchemical Transformations Involving Mass and Potential. , 2011, Journal of chemical theory and computation.

[21]  Soonmin Jang,et al.  Applications of higher order composite factorization schemes in imaginary time path integral simulations , 2001 .

[22]  J. Vaníček,et al.  Three applications of path integrals: equilibrium and kinetic isotope effects, and the temperature dependence of the rate constant of the [1,5] sigmatropic hydrogen shift in (Z)-1,3-pentadiene , 2010, Journal of molecular modeling.

[23]  Cristian Predescu.  Reconstruction of thermally symmetrized quantum autocorrelation functions from imaginary-time data. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[24]  Michele Ceriotti,et al.  Efficient first-principles calculation of the quantum kinetic energy and momentum distribution of nuclei. , 2012, Physical review letters.

[25]  Michele Ceriotti,et al.  Efficient methods and practical guidelines for simulating isotope effects. , 2012, The Journal of chemical physics.

[26]  J. Vaníček,et al.  Path integral evaluation of equilibrium isotope effects. , 2009, The Journal of chemical physics.

[27]  E. M.,et al.  Statistical Mechanics , 2021, Manual for Theoretical Chemistry.

[28]  M. Buchowiecki.  Quantum calculations of the temperature dependence of the rate constant and the equilibrium constant for the NH3 + H ⇌ NH2 + H2 reaction , 2012 .

[29]  H. G. Petersen,et al.  Error estimates on averages of correlated data , 1989 .

[30]  D. Ceperley.  Path integrals in the theory of condensed helium , 1995 .

[31]  H. Kleinert.  Path Integrals in Quantum Mechanics, Statistics, Polymer Physics, and Financial Markets , 2006 .

[32]  H. Trotter.  On the product of semi-groups of operators , 1959 .

[33]  Alejandro Pérez,et al.  Improving the convergence of closed and open path integral molecular dynamics via higher order Trotter factorization schemes. , 2011, The Journal of chemical physics.

[34]  Masatoshi Imada,et al.  Monte Carlo Calculation of Quantum Systems. II. Higher Order Correction , 1984 .

[35]  Siu A. Chin,et al.  Symplectic integrators from composite operator factorizations , 1997 .

[36]  Quantum statistical calculations and symplectic corrector algorithms. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[37]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[38]  Oliver Riordan,et al.  The inefficiency of re-weighted sampling and the curse of system size in high-order path integration , 2011, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[39]  H. Urey,et al.  The Equilibrium Between the Three Hydrogens , 1934 .

[40]  Higher order and infinite Trotter-number extrapolations in path integral Monte Carlo. , 2004, The Journal of chemical physics.