Accelerating equilibrium isotope effect calculations: I. Stochastic thermodynamic integration with respect to mass

Abstract

Accurate path integral Monte Carlo or molecular dynamics calculations of isotope effects have until recently been expensive because of the necessity to reduce three types of errors present in such calculations: statistical errors due to sampling, path integral discretization errors, and thermodynamic integration errors. While the statistical errors can be reduced with virial estimators and path integral discretization errors with high-order factorization of the Boltzmann operator, here we propose a method for accelerating isotope effect calculations by eliminating the integration error. We show that the integration error can be removed entirely by changing particle masses stochastically during the calculation and by using a piecewise linear umbrella biasing potential. Moreover, we demonstrate numerically that this approach does not increase the statistical error. The resulting acceleration of isotope effect calculations is demonstrated on a model harmonic system and on deuterated species of methane.

References in corpus (11)

• Reconstructing the Density of States by History-Dependent Metadynamics
• Efficient methods and practical guidelines for simulating isotope effects
• Path integral virial estimator based on the scaling of fluctuation coordinates: Application to quantum clusters with fourth-order propagators
• Optimal series representations for numerical path integral simulations
• Direct path integral estimators for isotope fractionation ratios
• Path integral evaluation of equilibrium isotope effects
• Monte Carlo evaluation of the equilibrium isotope effects using the Takahashi-Imada factorization of the Feynman path integral
• Accelerating quantum instanton calculations of the kinetic isotope effects
• 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
• Isotopic fractionation in proteins as a measure of hydrogen bond length
• Accelerating equilibrium isotope effect calculations: II. Stochastic implementation of direct estimators

Cited by in corpus (3)

• Nuclear quantum effects enter the mainstream
• Semiclassical analysis of the quantum instanton approximation
• Accelerating equilibrium isotope effect calculations: II. Stochastic implementation of direct estimators