Error and timing analysis of multiple time-step integration methods for molecular dynamics

Abstract Molecular dynamics simulations of biomolecules performed using multiple time-step integration methods are hampered by resonance instabilities. We analyze the properties of a simple 1D linear system integrated with the symplectic reference system propagator MTS (r-RESPA) technique following earlier work by others. A closed form expression for the time step dependent Hamiltonian which corresponds to r-RESPA integration of the model is derived. This permits us to present an analytic formula for the dependence of the integration accuracy on short-range force cutoff range. A detailed analysis of the force decomposition for the standard Ewald summation method is then given as the Ewald method is a good candidate to achieve high scaling on modern massively parallel machines. We test the new analysis on a realistic system, a protein in water. Under Langevin dynamics with a weak friction coefficient ( ζ = 1 ps −1 ) to maintain temperature control and using the SHAKE algorithm to freeze out high frequency vibrations, we show that the 5 fs resonance barrier present when all degrees of freedom are unconstrained is postponed to ≈ 12 fs . An iso-error boundary with respect to the short-range cutoff range and multiple time step size agrees well with the analytical results which are valid due to dominance of the high frequency modes in determining integrator accuracy. Using r-RESPA to treat the long range interactions results in a 6× increase in efficiency for the decomposition described in the text.

[1]  Masakatsu Watanabe,et al.  Dynamics of molecules with internal degrees of freedom by multiple time-step methods , 1993 .

[2]  M. Karplus,et al.  SIMULATIONS OF MACROMOLECULES BY MULTIPLE TIME-STEP METHODS , 1995 .

[3]  Robert S. Germain,et al.  Blue Matter, an application framework for molecular simulation on Blue Gene , 2003, J. Parallel Distributed Comput..

[4]  Nikos L. Doltsinis,et al.  Quantum Simulations of complex many-body systems:from theory to algorithms , 2002 .

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

[6]  P. Procacci,et al.  Taming the Ewald sum in molecular dynamics simulations of solvated proteins via a multiple time step algorithm , 1996 .

[7]  Godehard Sutmann,et al.  Long-Range Interactions in Many-Particle Simulation , 2002 .

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

[9]  J. M. Sanz-Serna,et al.  Numerical Hamiltonian Problems , 1994 .

[10]  K. Schulten,et al.  Difficulties with multiple time stepping and fast multipole algorithm in molecular dynamics , 1997 .

[11]  P. Kollman,et al.  Settle: An analytical version of the SHAKE and RATTLE algorithm for rigid water models , 1992 .

[12]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

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

[14]  Helmut Grubmüller,et al.  Multiple time step algorithms for molecular dynamics simulations of proteins: How good are they? , 1998 .

[15]  T. Darden,et al.  The effect of long‐range electrostatic interactions in simulations of macromolecular crystals: A comparison of the Ewald and truncated list methods , 1993 .

[16]  M. Dixon,et al.  Interionic potentials in alkali halides and their use in simulations of the molten salts , 1976 .

[17]  Klaus Schulten,et al.  Generalized Verlet Algorithm for Efficient Molecular Dynamics Simulations with Long-range Interactions , 1991 .

[18]  Tamar Schlick,et al.  Inherent speedup limitations in multiple time step/particle mesh Ewald algorithms , 2003, J. Comput. Chem..

[19]  L. Greengard.  The Rapid Evaluation of Potential Fields in Particle Systems , 1988 .

[20]  T. Schlick,et al.  Extrapolation versus impulse in multiple-timestepping schemes. II. Linear analysis and applications to Newtonian and Langevin dynamics , 1998 .

[21]  T Schlick,et al.  Biomolecular dynamics at long timesteps: bridging the timescale gap between simulation and experimentation. , 1997, Annual review of biophysics and biomolecular structure.

[22]  Yuefan Deng,et al.  Performance Models on QCDOC for Molecular Dynamics with Coulomb Potentials , 2004, Int. J. High Perform. Comput. Appl..

[23]  Robert D. Skeel,et al.  Dangers of multiple time step methods , 1993 .

[24]  H. C. Andersen.  Rattle: A “velocity” version of the shake algorithm for molecular dynamics calculations , 1983 .

[25]  B. Berne,et al.  A Multiple-Time-Step Molecular Dynamics Algorithm for Macromolecules , 1994 .

[26]  Jim Glosli,et al.  Comments on P3M, FMM, and the Ewald method for large periodic Coulombic systems , 1996 .

[27]  T. Schlick,et al.  Efficient multiple-time-step integrators with distance-based force splitting for particle-mesh-Ewald molecular dynamics simulations , 2002 .

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

[29]  J. Perram,et al.  Cutoff Errors in the Ewald Summation Formulae for Point Charge Systems , 1992 .

[30]  J. Board,et al.  Ewald summation techniques in perspective: a survey , 1996 .

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

[32]  T. Schlick,et al.  Overcoming stability limitations in biomolecular dynamics. I. Combining force splitting via extrapolation with Langevin dynamics in LN , 1998 .

[33]  H. G. Petersen,et al.  An algorithm for the simulation of condensed matter which grows as the 3/2 power of the number of particles , 1988 .

[34]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[35]  Laxmikant V. Kale,et al.  Algorithmic Challenges in Computational Molecular Biophysics , 1999 .

[36]  T. Darden,et al.  Molecular dynamics simulations of biomolecules: long-range electrostatic effects. , 1999, Annual review of biophysics and biomolecular structure.

[37]  O. Steinhauser,et al.  Cutoff size does strongly influence molecular dynamics results on solvated polypeptides. , 1992, Biochemistry.

[38]  T. Darden,et al.  Toward the Accurate Modeling of DNA: The Importance of Long-Range Electrostatics , 1995 .

[39]  D. Tildesley,et al.  Multiple time-step methods in molecular dynamics , 1978 .

[40]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[41]  Siu A Chin.  Dynamical multiple-time stepping methods for overcoming resonance instabilities. , 2004, The Journal of chemical physics.

[42]  Avram Hayli.  Le problème des N corps dans un champ extérieur. Application à l'évolution dynamique des amas ouverts. I , 1967 .

[43]  P. P. Ewald.  Die Berechnung optischer und elektrostatischer Gitterpotentiale , 1921 .

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

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

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

[47]  T Schlick,et al.  Time-trimming tricks for dynamic simulations: splitting force updates to reduce computational work. , 2001, Structure.

[48]  Timothy R. Forester,et al.  On Multiple Time-step Algorithms and the Ewald Sum , 1994 .