Stochastic, resonance-free multiple time-step algorithm for molecular dynamics with very large time steps

Molecular dynamics is one of the most commonly used approaches for studying the dynamics and statistical distributions of physical, chemical, and biological systems using atomistic or coarse-grained models. It is often the case, however, that the interparticle forces drive motion on many time scales, and the efficiency of a calculation is limited by the choice of time step, which must be sufficiently small that the fastest force components are accurately integrated. Multiple time-stepping algorithms partially alleviate this inefficiency by assigning to each time scale an appropriately chosen step-size. As the fast forces are often computationally cheaper to evaluate than the slow forces, this results in a significant gain in efficiency. However, such approaches are limited by resonance phenomena, wherein motion on the fastest time scales limits the step sizes associated with slower time scales. In atomistic models of biomolecular systems, for example, resonances limit the largest time step to around 5-6 fs. Stochastic processes promote mixing and ergodicity in dynamical systems and reduce the impact of resonant modes. In this paper, we introduce a set of stochastic isokinetic equations of motion that are shown to be rigorously ergodic, largely free of resonances, and can be integrated using a multiple time-stepping algorithm which is easily implemented in existing molecular dynamics codes. The technique is applied to a simple, illustrative problem and then to a more realistic system, namely, a flexible water model. Using this approach outer time steps as large as 100 fs are shown to be possible.

[1]  H. Yoshida.  Construction of higher order symplectic integrators , 1990 .

[2]  M. Suzuki,et al.  General theory of fractal path integrals with applications to many‐body theories and statistical physics , 1991 .

[3]  M. Klein,et al.  Nosé-Hoover chains : the canonical ensemble via continuous dynamics , 1992 .

[4]  Masuo Suzuki,et al.  General Nonsymmetric Higher-Order Decomposition of Exponential Operators and Symplectic Integrators , 1992 .

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

[6]  Richard L. Tweedie,et al.  Markov Chains and Stochastic Stability , 1993, Communications and Control Engineering Series.

[7]  Andrew E. Torda,et al.  Local elevation: A method for improving the searching properties of molecular dynamics simulation , 1994, J. Comput. Aided Mol. Des..

[8]  Grubmüller,et al.  Predicting slow structural transitions in macromolecular systems: Conformational flooding. , 1995, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[9]  B. Berne,et al.  A new molecular dynamics method combining the reference system propagator algorithm with a fast multipole method for simulating proteins and other complex systems , 1995 .

[10]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[11]  Mark E. Tuckerman,et al.  Explicit reversible integrators for extended systems dynamics , 1996 .

[12]  Ruhong Zhou,et al.  Molecular dynamics with multiple time scales: The selection of efficient reference system propagators , 1996 .

[13]  Fei Zhang.  Operator-splitting integrators for constant-temperature molecular dynamics , 1997 .

[14]  B. Berne,et al.  Large scale simulation of macromolecules in solution: Combining the periodic fast multipole method with multiple time step integrators , 1997 .

[15]  R. Skeel,et al.  Nonlinear Resonance Artifacts in Molecular Dynamics Simulations , 1998 .

[16]  G. Martyna,et al.  Electrostatic calculations and multiple time scales in molecular dynamics simulation of flexible molecular systems , 1998 .

[17]  M. Tuckerman,et al.  On the classical statistical mechanics of non-Hamiltonian systems , 1999 .

[18]  Mark E. Tuckerman,et al.  Exploiting multiple levels of parallelism in Molecular Dynamics based calculations via modern techniques and software paradigms on distributed memory computers , 2000 .

[19]  Eric F Darve,et al.  Calculating free energies using average force , 2001 .

[20]  Andrew G. Glen,et al.  APPL , 2001 .

[21]  G. Ciccotti,et al.  Non-Hamiltonian molecular dynamics: Generalizing Hamiltonian phase space principles to non-Hamiltonian systems , 2001 .

[22]  A. Laio,et al.  Escaping free-energy minima , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[23]  M. Tuckerman,et al.  On the use of the adiabatic molecular dynamics technique in the calculation of free energy profiles , 2002 .

[24]  M. Tuckerman,et al.  An Adiabatic Molecular Dynamics Method for the Calculation of Free Energy Profiles , 2002 .

[25]  Jonathan C. Mattingly,et al.  Ergodicity for SDEs and approximations: locally Lipschitz vector fields and degenerate noise , 2002 .

[26]  Jesús A. Izaguirre,et al.  Verlet-I/R-RESPA/Impulse is Limited by Nonlinear Instabilities , 2003, SIAM J. Sci. Comput..

[27]  Mark E. Tuckerman,et al.  Algorithms and novel applications based on the isokinetic ensemble. I. Biophysical and path integral molecular dynamics , 2003 .

[28]  C. Chipot,et al.  Overcoming free energy barriers using unconstrained molecular dynamics simulations. , 2004, The Journal of chemical physics.

[29]  M. Tuckerman,et al.  Long time molecular dynamics for enhanced conformational sampling in biomolecular systems. , 2004, Physical review letters.

[30]  Lula Rosso,et al.  Mapping the backbone dihedral free-energy surfaces in small peptides in solution using adiabatic free-energy dynamics. , 2005, The journal of physical chemistry. B.

[31]  S. Attal,et al.  Open Quantum Systems II , 2006 .

[32]  Wei Zhang,et al.  An accurate and simple quantum model for liquid water. , 2006, The Journal of chemical physics.

[33]  E. Vanden-Eijnden,et al.  A temperature accelerated method for sampling free energy and determining reaction pathways in rare events simulations , 2006 .

[34]  Luc Rey Bellet.  Ergodic Properties of Markov Processes , 2006 .

[35]  Luc Rey-Bellet,et al.  Ergodic properties of Markov processes , 2006 .

[36]  M. Chaplain,et al.  Thermostats for “Slow” Configurational Modes , 2004, physics/0412163.

[37]  Peter Minary,et al.  Dynamical Spatial Warping: A Novel Method for the Conformational Sampling of Biophysical Structure , 2008, SIAM J. Sci. Comput..

[38]  Eric Darve,et al.  Adaptive biasing force method for scalar and vector free energy calculations. , 2008, The Journal of chemical physics.

[39]  M. Tuckerman,et al.  Efficient and direct generation of multidimensional free energy surfaces via adiabatic dynamics without coordinate transformations. , 2008, The journal of physical chemistry. B.

[40]  Gary P. Morriss,et al.  Statistical Mechanics of Nonequilibrium Liquids , 2008 .

[41]  I. P. Omelyan,et al.  Advanced multiple time scale molecular dynamics , 2009 .

[42]  B. Leimkuhler,et al.  A Gentle Stochastic Thermostat for Molecular Dynamics , 2009 .

[43]  M. Tuckerman.  Statistical Mechanics: Theory and Molecular Simulation , 2010 .

[44]  Ruhong Zhou,et al.  Molecular Dynamics with Multiple Time Scales: How to Avoid Pitfalls. , 2010, Journal of chemical theory and computation.

[45]  Frédéric Legoll,et al.  Free energy calculations: an efficient adaptive biasing potential method. , 2009, The journal of physical chemistry. B.

[46]  Andriy Kovalenko,et al.  Multiple time scale molecular dynamics for fluids with orientational degrees of freedom. II. Canonical and isokinetic ensembles. , 2011, The Journal of chemical physics.

[47]  Mark E Tuckerman,et al.  Temperature-accelerated method for exploring polymorphism in molecular crystals based on free energy. , 2011, Physical review letters.

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

[49]  M. Tuckerman,et al.  Heating and flooding: a unified approach for rapid generation of free energy surfaces. , 2012, The Journal of chemical physics.

[50]  B. Leimkuhler,et al.  Rational Construction of Stochastic Numerical Methods for Molecular Sampling , 2012, 1203.5428.