Fine tuning classical and quantum molecular dynamics using a generalized Langevin equation.

Generalized Langevin Equation (GLE) thermostats have been used very effectively as a tool to manipulate and optimize the sampling of thermodynamic ensembles and the associated static properties. Here we show that a similar, exquisite level of control can be achieved for the dynamical properties computed from thermostatted trajectories. We develop quantitative measures of the disturbance induced by the GLE to the Hamiltonian dynamics of a harmonic oscillator, and show that these analytical results accurately predict the behavior of strongly anharmonic systems. We also show that it is possible to correct, to a significant extent, the effects of the GLE term onto the corresponding microcanonical dynamics, which puts on more solid grounds the use of non-equilibrium Langevin dynamics to approximate quantum nuclear effects and could help improve the prediction of dynamical quantities from techniques that use a Langevin term to stabilize dynamics. Finally we address the use of thermostats in the context of approximate path-integral-based models of quantum nuclear dynamics. We demonstrate that a custom-tailored GLE can alleviate some of the artifacts associated with these techniques, improving the quality of results for the modeling of vibrational dynamics of molecules, liquids, and solids.

[1]  R. Baer,et al.  Self-averaging stochastic Kohn-Sham density-functional theory. , 2013, Physical review letters.

[2]  Jianshu Cao,et al.  The formulation of quantum statistical mechanics based on the Feynman path centroid density. IV. Algorithms for centroid molecular dynamics , 1994 .

[3]  Michele Ceriotti,et al.  Efficient multiple time scale molecular dynamics: Using colored noise thermostats to stabilize resonances. , 2010, The Journal of chemical physics.

[4]  Luc Rey-Bellet,et al.  Uncertainty quantification for generalized Langevin dynamics. , 2016, The Journal of chemical physics.

[5]  J. Tennyson,et al.  A PROGRAM SUITE FOR THE CALCULATION OF RO-VIBRATIONAL SPECTRA OF TRIATOMIC-MOLECULES , 1989 .

[6]  Harry Partridge,et al.  The determination of an accurate isotope dependent potential energy surface for water from extensive ab initio calculations and experimental data , 1997 .

[7]  Timothy John Harvey Hele.  How to obtain thermostatted ring polymer molecular dynamics from exact quantum dynamics and when to use it , 2015, 1508.05820.

[8]  Ali Hassanali,et al.  Nuclear quantum effects in a HIV/cancer inhibitor: The case of ellipticine. , 2016, The Journal of chemical physics.

[9]  G. Mazzola,et al.  Unexpectedly high pressure for molecular dissociation in liquid hydrogen by electronic simulation , 2014, Nature Communications.

[10]  Michele Parrinello,et al.  Generalized neural-network representation of high-dimensional potential-energy surfaces. , 2007, Physical review letters.

[11]  Timothy J. H. Hele,et al.  On the relation between thermostatted ring-polymer molecular dynamics and exact quantum dynamics , 2016 .

[12]  M. Ceriotti.  A novel framework for enhanced molecular dynamics based on the generalized Langevin equation , 2010 .

[13]  Guohui Li,et al.  The Vibrational Level Spectrum of H2O(X̃1A′) on a Partridge–Schwenke Potential up to the Dissociation Limit , 2001 .

[14]  Chem. , 2020, Catalysis from A to Z.

[15]  Michele Parrinello,et al.  The δ-thermostat: selective normal-modes excitation by colored-noise Langevin dynamics , 2010 .

[16]  T. Morawietz,et al.  How van der Waals interactions determine the unique properties of water , 2016, Proceedings of the National Academy of Sciences.

[17]  Michele Ceriotti,et al.  i-PI: A Python interface for ab initio path integral molecular dynamics simulations , 2014, Comput. Phys. Commun..

[18]  Matthias Krack,et al.  Efficient and accurate Car-Parrinello-like approach to Born-Oppenheimer molecular dynamics. , 2007, Physical review letters.

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

[20]  L. Stella,et al.  Generalized Langevin equation: An efficient approach to nonequilibrium molecular dynamics of open systems , 2013, 1312.4903.

[21]  Robert Zwanzig,et al.  Memory Effects in Irreversible Thermodynamics , 1961 .

[22]  Michele Ceriotti,et al.  Communication: On the consistency of approximate quantum dynamics simulation methods for vibrational spectra in the condensed phase. , 2014, The Journal of chemical physics.

[23]  Michael J. Willatt,et al.  Communication: Relation of centroid molecular dynamics and ring-polymer molecular dynamics to exact quantum dynamics. , 2015, The Journal of chemical physics.

[24]  M. Daube-Witherspoon,et al.  An Iterative Image Space Reconstruction Algorthm Suitable for Volume ECT , 1986, IEEE Transactions on Medical Imaging.

[25]  Jörg Behler,et al.  Nuclear Quantum Effects in Water at the Triple Point: Using Theory as a Link Between Experiments. , 2016, The journal of physical chemistry letters.

[26]  L. Colombo,et al.  Simulating Energy Relaxation in Pump-Probe Vibrational Spectroscopy of Hydrogen-Bonded Liquids. , 2017, Journal of chemical theory and computation.

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

[28]  D. M. Titterington,et al.  On some Bayesian/regularization methods for image restoration , 1995, IEEE Trans. Image Process..

[29]  George C. Schatz,et al.  The journal of physical chemistry letters , 2009 .

[30]  Michele Parrinello,et al.  Linear scaling electronic structure calculations and accurate statistical mechanics sampling with noisy forces , 2006 .

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

[32]  Ondrej Marsalek,et al.  Quantum Dynamics and Spectroscopy of Ab Initio Liquid Water: The Interplay of Nuclear and Electronic Quantum Effects. , 2017, The journal of physical chemistry letters.

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

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

[35]  G. A. Pavliotis,et al.  Exponential return to equilibrium for hypoelliptic quadratic systems , 2011, 1106.2326.

[36]  David E Manolopoulos,et al.  Comparison of path integral molecular dynamics methods for the infrared absorption spectrum of liquid water. , 2008, The Journal of chemical physics.

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

[38]  Jean-Jacques Greffet,et al.  Quantum thermal bath for molecular dynamics simulation. , 2009, Physical review letters.

[39]  Jörg Behler,et al.  High order path integrals made easy. , 2016, The Journal of chemical physics.

[40]  H. Haubeck.  COMP , 2019, Springer Reference Medizin.

[41]  Alejandro Pérez,et al.  A comparative study of the centroid and ring-polymer molecular dynamics methods for approximating quantum time correlation functions from path integrals. , 2009, The Journal of chemical physics.

[42]  Oriol Vendrell,et al.  Full dimensional (15-dimensional) quantum-dynamical simulation of the protonated water dimer. II. Infrared spectrum and vibrational dynamics. , 2007, The Journal of chemical physics.

[43]  Joel M Bowman,et al.  Ab initio potential energy and dipole moment surfaces for H5O2 +. , 2005, The Journal of chemical physics.

[44]  Michele Ceriotti,et al.  Nuclear Quantum Effects in H(+) and OH(-) Diffusion along Confined Water Wires. , 2016, The journal of physical chemistry letters.

[45]  Michael F. Insana.  IEEE Transactions on Medical Imaging information for authors , 2016 .

[46]  D. Manolopoulos,et al.  How to remove the spurious resonances from ring polymer molecular dynamics. , 2014, The Journal of chemical physics.

[47]  Oriol Vendrell,et al.  Full dimensional (15-dimensional) quantum-dynamical simulation of the protonated water dimer. II. Infrared spectrum and vibrational dynamics. , 2007, The Journal of chemical physics.

[48]  S. Grimme,et al.  A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. , 2010, The Journal of chemical physics.

[49]  B. Berne,et al.  Comparison of rate theories for generalized Langevin dynamics , 1991 .

[50]  Dominik Marx,et al.  On the applicability of centroid and ring polymer path integral molecular dynamics for vibrational spectroscopy. , 2009, The Journal of chemical physics.

[51]  M. Wilkinson,et al.  Spherical Ornstein-Uhlenbeck Processes , 2011, 1105.4903.