Efficient calculation of free energy differences associated with isotopic substitution using path-integral molecular dynamics - Archive ouverte HAL Accéder directement au contenu
Article Dans Une Revue Journal of Chemical Theory and Computation Année : 2014

Efficient calculation of free energy differences associated with isotopic substitution using path-integral molecular dynamics

Résumé

The problem of computing free energy differences due to isotopic substitution in chemical systems is discussed. The shift in the equilibrium properties of a system upon isotopic substitution is a purely quantum mechanical effect that can be quantified using the Feynman path integral approach. In this paper, we explore two developments that lead to a highly efficient path integral scheme. First, we employ a mass switching function inspired by the work of Ceriotti and Markland [ J. Chem. Phys. 2013, 138, 014112] that is based on the inverse square root of the mass and which leads to a perfectly constant free energy derivative with respect to the switching parameter in the harmonic limit. We show that even for anharmonic systems, this scheme allows a single-point thermodynamic integration approach to be used in the construction of free energy differences. In order to improve the efficiency of the calculations even further, however, we derive a set of free energy derivative estimators based on the fourth-order scheme of Takahashi and Imada [ J. Phys. Soc. Jpn. 1984, 53, 3765]. The Takahashi-Imada procedure generates a primitive fourth-order estimator that allows the number of imaginary time slices in the path-integral approach to be reduced substantially. However, as with all primitive estimators, its convergence is plagued by numerical noise. In order to alleviate this problem, we derive a fourth-order virial estimator based on a transferring of the difference between second- and fourth-order primitive estimators, which remains relatively constant as a function of the number of configuration samples, to the second-order virial estimator. We show that this new estimator converges as smoothly as the second-order virial estimator but requires significantly fewer imaginary time points.
Fichier non déposé

Dates et versions

hal-01765162 , version 1 (12-04-2018)

Identifiants

Citer

Ondrej Marsalek, Pei-Yang Chen, Romain Dupuis, Magali Benoit, Merlin Méheut, et al.. Efficient calculation of free energy differences associated with isotopic substitution using path-integral molecular dynamics. Journal of Chemical Theory and Computation, 2014, 10 (4), pp.1440-1453. ⟨10.1021/ct400911m⟩. ⟨hal-01765162⟩
112 Consultations
0 Téléchargements

Altmetric

Partager

Gmail Facebook X LinkedIn More