Accelerating the convergence of path integral dynamics with a generalized Langevin equation.

The quantum nature of nuclei plays an important role in the accurate modelling of light atoms such as hydrogen, but it is often neglected in simulations due to the high computational overhead involved. It has recently been shown that zero-point energy effects can be included comparatively cheaply in simulations of harmonic and quasiharmonic systems by augmenting classical molecular dynamics with a generalized Langevin equation (GLE). Here we describe how a similar approach can be used to accelerate the convergence of path integral (PI) molecular dynamics to the exact quantum mechanical result in more strongly anharmonic systems exhibiting both zero point energy and tunnelling effects. The resulting PI-GLE method is illustrated with applications to a double-well tunnelling problem and to liquid water.

[1]  David E Manolopoulos,et al.  On the short-time limit of ring polymer molecular dynamics. , 2006, The Journal of chemical physics.

[2]  Jianshu Cao,et al.  The formulation of quantum statistical mechanics based on the Feynman path centroid density. II. Dynamical properties , 1994 .

[3]  Ian R. Craig,et al.  Quantum statistics and classical mechanics: real time correlation functions from ring polymer molecular dynamics. , 2004, The Journal of chemical physics.

[4]  Dominik Endres,et al.  A new metric for probability distributions , 2003, IEEE Transactions on Information Theory.

[5]  Michele Parrinello,et al.  Nuclear quantum effects in ab initio dynamics: Theory and experiments for lithium imide , 2010, 1009.2862.

[6]  B. Berne,et al.  Efficient molecular dynamics and hybrid Monte Carlo algorithms for path integrals , 1993 .

[7]  Mark E. Tuckerman,et al.  Reversible multiple time scale molecular dynamics , 1992 .

[8]  David E Manolopoulos,et al.  Zero point energy leakage in condensed phase dynamics: an assessment of quantum simulation methods for liquid water. , 2009, The Journal of chemical physics.

[9]  M. Kalos.  Monte Carlo Methods in Quantum Problems , 1984 .

[10]  R. Zwanzig.  Nonequilibrium statistical mechanics , 2001, Physics Subject Headings (PhySH).

[11]  M. Parrinello,et al.  Canonical sampling through velocity rescaling. , 2007, The Journal of chemical physics.

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

[13]  Peter G. Wolynes,et al.  Exploiting the isomorphism between quantum theory and classical statistical mechanics of polyatomic fluids , 1981 .

[14]  J. A. Barker.  A quantum‐statistical Monte Carlo method; path integrals with boundary conditions , 1979 .

[15]  Giovanni Bussi,et al.  Colored-Noise Thermostats à la Carte , 2010, 1204.0822.

[16]  D. Manolopoulos,et al.  An efficient ring polymer contraction scheme for imaginary time path integral simulations. , 2008, The Journal of chemical physics.

[17]  C. DeWitt-Morette,et al.  Techniques and Applications of Path Integration , 1981 .

[18]  Roberto Car,et al.  Nuclear quantum effects in water. , 2008, Physical review letters.

[19]  Voth.  Calculation of equilibrium averages with Feynman-Hibbs effective classical potentials and similar variational approximations. , 1991, Physical review. A, Atomic, molecular, and optical physics.

[20]  Michele Parrinello,et al.  Displaced path integral formulation for the momentum distribution of quantum particles. , 2010, Physical review letters.

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

[22]  Thomas E. Markland,et al.  Competing quantum effects in the dynamics of a flexible water model. , 2009, The Journal of chemical physics.

[23]  C. Gardiner.  Handbook of Stochastic Methods , 1983 .

[24]  Michele Parrinello,et al.  Efficient stochastic thermostatting of path integral molecular dynamics. , 2010, The Journal of chemical physics.

[25]  M. Shiga,et al.  Calculation of heat capacities of light and heavy water by path-integral molecular dynamics. , 2005, The Journal of chemical physics.

[26]  Thomas E Markland,et al.  A fast path integral method for polarizable force fields. , 2009, The Journal of chemical physics.

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

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

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

[30]  David E. Manolopoulos,et al.  A refined ring polymer contraction scheme for systems with electrostatic interactions , 2008 .

[31]  Giovanni Bussi,et al.  Langevin equation with colored noise for constant-temperature molecular dynamics simulations. , 2008, Physical review letters.

[32]  Gregory A. Voth,et al.  A derivation of centroid molecular dynamics and other approximate time evolution methods for path integral centroid variables , 1999 .

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

[34]  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.

[35]  Rustam Z. Khaliullin,et al.  Graphite-diamond phase coexistence study employing a neural-network mapping of the ab initio potential energy surface , 2010 .

[36]  Giovanni Bussi,et al.  Nuclear quantum effects in solids using a colored-noise thermostat. , 2009, Physical review letters.