Efficient methods and practical guidelines for simulating isotope effects

TL;DR

This paper introduces efficient path integral quantum mechanics estimators for simulating isotope effects, outperforming traditional methods in accuracy and computational cost, with practical guidelines for their application in liquid systems.

Contribution

It presents novel estimators based on free energy perturbation that enable isotope effect calculations from a single trajectory, improving efficiency and accuracy over existing approaches.

Findings

Single-trajectory estimators reduce computational cost

Enhanced accuracy over quasi-harmonic approximations

Guidelines for method selection in various scenarios

Abstract

The shift in chemical equilibria due to isotope substitution is often exploited to gain insight into a wide variety of chemical and physical processes. It is a purely quantum mechanical effect, which can be computed exactly using simulations based on the path integral formalism. Here we discuss how these techniques can be made dramatically more efficient, and how they ultimately outperform quasi-harmonic approximations to treat quantum liquids not only in terms of accuracy, but also in terms of computational efficiency. To achieve this goal we introduce path integral quantum mechanics estimators based on free energy perturbation, which enable the evaluation of isotope effects using only a single path integral molecular dynamics trajectory of the naturally abundant isotope. We use as an example the calculation of the free energy change associated with H/D and 16O/18O substitutions in liquid water, and of the fractionation of those isotopes between the liquid and the vapor phase. In doing so, we demonstrate and discuss quantitatively the relative benefits of each approach, thereby providing a set of guidelines that should facilitate the choice of the most appropriate method in different, commonly encountered scenarios. The efficiency of the estimators we introduce and the analysis that we perform should in particular facilitate accurate ab initio calculation of isotope effects in condensed phase systems.