Extrapolation versus impulse in multiple-timestepping schemes. II. Linear analysis and applications to Newtonian and Langevin dynamics
- 1 August 1998
- journal article
- research article
- Published by AIP Publishing in The Journal of Chemical Physics
- Vol. 109 (5), 1633-1642
- https://doi.org/10.1063/1.476737
Abstract
Force splitting or multiple timestep (MTS) methods are effective techniques that accelerate biomolecular dynamics simulations by updating the fast and slow forces at different frequencies. Since simple extrapolation formulas for incorporating the slow forces into the, discretization produced notable energy drifts, symplectic MTS variants based on periodic impulses became more popular. However, the efficiency gain possible with these impulse approaches, is limited by a timestep barrier due to resonance-a numerical artifact occurring when the timestep is related to the period of the fastest motion present in the dynamics. This limitation is lifted substantially for MTS methods based on extrapolation in combination with stochastic dynamics, as demonstrated for the LN method in the companion paper for protein dynamics. To explain our observations on those complex nonlinear systems, we examine here the stability of extrapolation and impulses to force-splitting in Newtonian and Langevin dynamics. We analyze for a simple linear test system the energy drift of the former and the resonance-related artifacts of the latter technique. We show that two-class impulse methods are generally stable except at integer multiples of half the period of the fastest motion, with the severity of the instability worse at larger timesteps. Extrapolation methods are generally unstable for the Newtonian model problem, but the instability is bounded for increasing timesteps. This boundedness ensures good long-timestep behavior of extrapolation methods for Langevin dynamics with moderate values of the collision parameter. We thus advocate extrapolation methods for efficient integration of the stochastic Langevin equations of motion, as in the LN method described in paper I. (C) 1998 American Institute of Physics.This publication has 28 references indexed in Scilit:
- BIOMOLECULAR DYNAMICS AT LONG TIMESTEPS:Bridging the Timescale Gap Between Simulation and ExperimentationAnnual Review of Biophysics, 1997
- A Multiple-Time-Step Molecular Dynamics Algorithm for MacromoleculesThe Journal of Physical Chemistry, 1994
- Finite representation of an infinite bulk system: Solvent boundary potential for computer simulationsThe Journal of Chemical Physics, 1994
- LIN: A new algorithm to simulate the dynamics of biomolecules by combining implicit‐integration and normal mode techniquesJournal of Computational Chemistry, 1993
- Unfolding and refolding of the native structure of bovine pancreatic trypsin inhibitor studied by computer simulationsBiochemistry, 1993
- Multiple Time Steps: Limits on the Speedup of Molecular Dynamics Simulations of Aqueous SystemsMolecular Simulation, 1993
- Reversible multiple time scale molecular dynamicsThe Journal of Chemical Physics, 1992
- Langevin dynamics of peptides: The frictional dependence of isomerization rates of N‐acetylalanyl‐N′‐methylamideBiopolymers, 1992
- Higher-order implicit strong numerical schemes for stochastic differential equationsJournal of Statistical Physics, 1992
- Generalized Verlet Algorithm for Efficient Molecular Dynamics Simulations with Long-range InteractionsMolecular Simulation, 1991