English

Accelerating equilibrium isotope effect calculations: II. Stochastic implementation of direct estimators

Chemical Physics 2024-09-26 v2 Computational Physics Quantum Physics

Abstract

Path integral calculations of equilibrium isotope effects and isotopic fractionation are expensive due to the presence of path integral discretization errors, statistical errors, and thermodynamic integration errors. Whereas the discretization errors can be reduced by high-order factorization of the path integral and statistical errors by using centroid virial estimators, two recent papers proposed alternative ways to completely remove the thermodynamic integration errors: Cheng and Ceriotti [J. Chem. Phys. 141, 244112 (2015)] employed a variant of free-energy perturbation called "direct estimators," while Karandashev and Van\'{\i}\v{c}ek [J. Chem. Phys. 143, 194104 (2017)] combined the thermodynamic integration with a stochastic change of mass and piecewise-linear umbrella biasing potential. Here we combine the former approach with the stochastic change of mass in order to decrease its statistical errors when applied to larger isotope effects, and perform a thorough comparison of different methods by computing isotope effects first on a harmonic model, and then on methane and methanium, where we evaluate all isotope effects of the form CH4xDx/CH4\mathrm{CH}_{\mathrm{4-x}}\mathrm{D}_{\mathrm{x}}/\mathrm{CH}_{4} and CH5xDx+/CH5+\mathrm{CH}_{\mathrm{5-x}}\mathrm{D}^{+}_{\mathrm{x}}/\mathrm{CH}^{+}_{5}, respectively. We discuss thoroughly the reasons for a surprising behavior of the original method of direct estimators, which performed well for a much larger range of isotope effects than what had been expected previously.

Keywords

Cite

@article{arxiv.1909.02910,
  title  = {Accelerating equilibrium isotope effect calculations: II. Stochastic implementation of direct estimators},
  author = {Konstantin Karandashev and Jiri Vanicek},
  journal= {arXiv preprint arXiv:1909.02910},
  year   = {2024}
}

Comments

Please find Supplementary material for the article attached as an ancillary file available either by extracting the archive provided by "Other format" download option or by following the link under "Ancillary files."

R2 v1 2026-06-23T11:07:47.415Z