A Stochastic, Resonance-Free Multiple Time-Step Algorithm for Polarizable Models That Permits Very Large Time Steps.

Molecular dynamics remains one of the most widely used computational tools in the theoretical molecular sciences to sample an equilibrium ensemble distribution and/or to study the dynamical properties of a system. The efficiency of a molecular dynamics calculation is limited by the size of the time step that can be employed, which is dictated by the highest frequencies in the system. However, many properties of interest are connected to low-frequency, long time-scale phenomena, requiring many small time steps to capture. This ubiquitous problem can be ameliorated by employing multiple time-step algorithms, which assign different time steps to forces acting on different time scales. In such a scheme, fast forces are evaluated more frequently than slow forces, and as the former are often computationally much cheaper to evaluate, the savings can be significant. Standard multiple time-step approaches are limited, however, by resonance phenomena, wherein motion on the fastest time scales limits the step sizes that can be chosen for the slower time scales. In atomistic models of biomolecular systems, for example, the largest time step is typically limited to around 5 fs. Previously, we introduced an isokinetic extended phase-space algorithm (Minary et al. Phys. Rev. Lett. 2004, 93, 150201) and its stochastic analog (Leimkuhler et al. Mol. Phys. 2013, 111, 3579) that eliminate resonance phenomena through a set of kinetic energy constraints. In simulations of a fixed-charge flexible model of liquid water, for example, the time step that could be assigned to the slow forces approached 100 fs. In this paper, we develop a stochastic isokinetic algorithm for multiple time-step molecular dynamics calculations using a polarizable model based on fluctuating dipoles. The scheme developed here employs two sets of induced dipole moments, specifically, those associated with short-range interactions and those associated with a full set of interactions. The scheme is demonstrated on the polarizable AMOEBA water model. As was seen with fixed-charge models, we are able to obtain large time steps exceeding 100 fs, allowing calculations to be performed 10 to 20 times faster than standard thermostated molecular dynamics.

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

[2]  Mark E. Tuckerman,et al.  Molecular dynamics algorithm for condensed systems with multiple time scales , 1990 .

[3]  B. Leimkuhler,et al.  Stochastic, resonance-free multiple time-step algorithm for molecular dynamics with very large time steps , 2013, 1307.1167.

[4]  Andriy Kovalenko,et al.  Overcoming the Barrier on Time Step Size in Multiscale Molecular Dynamics Simulation of Molecular Liquids. , 2012, Journal of chemical theory and computation.

[5]  R. Elber,et al.  Computing time scales from reaction coordinates by milestoning. , 2004, The Journal of chemical physics.

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

[7]  Christophe Chipot,et al.  Exploring Multidimensional Free Energy Landscapes Using Time-Dependent Biases on Collective Variables. , 2010, Journal of chemical theory and computation.

[8]  Mark E. Tuckerman,et al.  Molecular dynamics algorithm for multiple time scales: Systems with disparate masses , 1991 .

[9]  Jason Crain,et al.  Signature properties of water: Their molecular electronic origins , 2015, Proceedings of the National Academy of Sciences.

[10]  E. Vanden-Eijnden,et al.  Single-sweep methods for free energy calculations. , 2007, The Journal of chemical physics.

[11]  T. Schlick,et al.  Masking Resonance Artifacts in Force-Splitting Methods for Biomolecular Simulations by Extrapolative Langevin Dynamics , 1999 .

[12]  Pengyu Y. Ren,et al.  Temperature and Pressure Dependence of the AMOEBA Water Model , 2004 .

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

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

[15]  Jonathan Smith.  Constructing a Small Place , 1998 .

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

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

[18]  J. Izaguirre.  Longer Time Steps for Molecular Dynamics , 1999 .

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

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

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

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

[23]  Michele Parrinello,et al.  Integrating the Car–Parrinello equations. II. Multiple time scale techniques , 1994 .

[24]  Thomas E Markland,et al.  Multiple time step integrators in ab initio molecular dynamics. , 2014, The Journal of chemical physics.

[25]  Jií Kolafa,et al.  Time‐reversible always stable predictor–corrector method for molecular dynamics of polarizable molecules , 2004, J. Comput. Chem..

[26]  Pengyu Y. Ren,et al.  Polarizable Atomic Multipole Water Model for Molecular Mechanics Simulation , 2003 .

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

[28]  R. Skeel,et al.  Langevin stabilization of molecular dynamics , 2001 .

[29]  Celeste Sagui,et al.  Towards an accurate representation of electrostatics in classical force fields: efficient implementation of multipolar interactions in biomolecular simulations. , 2004, The Journal of chemical physics.

[30]  Edward D Harder,et al.  Efficient multiple time step method for use with Ewald and particle mesh Ewald for large biomolecular systems , 2001 .

[31]  Y. Sugita,et al.  Replica-exchange molecular dynamics method for protein folding , 1999 .

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

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

[34]  Benoît Roux,et al.  Modeling induced polarization with classical Drude oscillators: Theory and molecular dynamics simulation algorithm , 2003 .

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

[36]  Glenn J. Martyna,et al.  A unified formalism for many-body polarization and dispersion: The quantum Drude model applied to fluid xenon , 2006 .

[37]  Jason Crain,et al.  Electronically coarse-grained molecular dynamics using quantum Drude oscillators , 2013 .

[38]  Alexander D. MacKerell,et al.  Implementation of extended Lagrangian dynamics in GROMACS for polarizable simulations using the classical Drude oscillator model , 2015, J. Comput. Chem..

[39]  Ryan P Steele,et al.  Multiple-Timestep ab Initio Molecular Dynamics Using an Atomic Basis Set Partitioning. , 2015, The journal of physical chemistry. A.

[40]  Igor Omelyan,et al.  Interpretation of atomic motion in flexible molecules: accelerating molecular dynamics simulations. , 2012, Physical review. E, Statistical, nonlinear, and soft matter physics.

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

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

[43]  Teresa Head-Gordon,et al.  An efficient and stable hybrid extended Lagrangian/self-consistent field scheme for solving classical mutual induction. , 2015, The Journal of chemical physics.

[44]  Wei Wang,et al.  Fast evaluation of polarizable forces. , 2005, The Journal of chemical physics.

[45]  Michel Masella,et al.  The multiple time step r-RESPA procedure and polarizable potentials based on induced dipole moments , 2006 .

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

[47]  Ron Elber,et al.  Extending molecular dynamics time scales with milestoning: example of complex kinetics in a solvated peptide. , 2007, The Journal of chemical physics.

[48]  Ryan P Steele,et al.  Multiple-Time Step Ab Initio Molecular Dynamics Based on Two-Electron Integral Screening. , 2015, Journal of chemical theory and computation.

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

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

[51]  Joan-Emma Shea,et al.  The Thermodynamics of Folding of a β Hairpin Peptide Probed Through Replica Exchange Molecular Dynamics Simulations , 2006 .

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

[53]  Ryan P Steele,et al.  Communication: Multiple-timestep ab initio molecular dynamics with electron correlation. , 2013, The Journal of chemical physics.

[54]  Pengyu Y. Ren,et al.  Consistent treatment of inter‐ and intramolecular polarization in molecular mechanics calculations , 2002, J. Comput. Chem..

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

[56]  Denis J. Evans,et al.  Non-Newtonian molecular dynamics , 1984 .

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

[58]  Alexander D. MacKerell,et al.  A simple polarizable model of water based on classical Drude oscillators , 2003 .

[59]  Anders M N Niklasson,et al.  Extended Born-Oppenheimer molecular dynamics. , 2008, Physical review letters.

[60]  Mark E. Tuckerman,et al.  Molecular dynamics algorithm for multiple time scales: Systems with long range forces , 1991 .

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

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