Langevin stabilization of molecular dynamics
Top Cited Papers
- 24 January 2001
- journal article
- research article
- Published by AIP Publishing in The Journal of Chemical Physics
- Vol. 114 (5), 2090-2098
- https://doi.org/10.1063/1.1332996
Abstract
In this paper we show the possibility of using very mild stochastic damping to stabilize long time step integrators for Newtonian molecular dynamics. More specifically, stable and accurate integrations are obtained for damping coefficients that are only a few percent of the natural decay rate of processes of interest, such as the velocity autocorrelation function. Two new multiple time stepping integrators, Langevin Molly (LM) and Brünger–Brooks–Karplus–Molly (BBK–M), are introduced in this paper. Both use the mollified impulse method for the Newtonian term. LM uses a discretization of the Langevin equation that is exact for the constant force, and BBK–M uses the popular Brünger–Brooks–Karplus integrator (BBK). These integrators, along with an extrapolative method called LN, are evaluated across a wide range of damping coefficient values. When large damping coefficients are used, as one would for the implicit modeling of solvent molecules, the method LN is superior, with LM closely following. However, with mild damping of 0.2 ps−1, LM produces the best results, allowing long time steps of 14 fs in simulations containing explicitly modeled flexible water. With BBK–M and the same damping coefficient, time steps of 12 fs are possible for the same system. Similar results are obtained for a solvated protein–DNA simulation of estrogen receptor ER with estrogen response element ERE. A parallel version of BBK–M runs nearly three times faster than the Verlet-I/r-RESPA (reversible reference system propagator algorithm) when using the largest stable time step on each one, and it also parallelizes well. The computation of diffusion coefficients for flexible water and ER/ERE shows that when mild damping of up to 0.2 ps−1 is used the dynamics are not significantly distorted.Keywords
This publication has 23 references indexed in Scilit:
- Longer time steps for molecular dynamicsThe Journal of Chemical Physics, 1999
- Masking Resonance Artifacts in Force-Splitting Methods for Biomolecular Simulations by Extrapolative Langevin DynamicsJournal of Computational Physics, 1999
- Extrapolation versus impulse in multiple-timestepping schemes. II. Linear analysis and applications to Newtonian and Langevin dynamicsThe Journal of Chemical Physics, 1998
- Overcoming stability limitations in biomolecular dynamics. I. Combining force splitting via extrapolation with Langevin dynamics in LNThe Journal of Chemical Physics, 1998
- Long-Time-Step Methods for Oscillatory Differential EquationsSIAM Journal on Scientific Computing, 1998
- Difficulties with multiple time stepping and fast multipole algorithm in molecular dynamicsJournal of Computational Chemistry, 1997
- Dangers of Multiple Time Step MethodsJournal of Computational Physics, 1993
- Reversible multiple time scale molecular dynamicsThe Journal of Chemical Physics, 1992
- Generalized Verlet Algorithm for Efficient Molecular Dynamics Simulations with Long-range InteractionsMolecular Simulation, 1991
- An analysis of the accuracy of Langevin and molecular dynamics algorithmsMolecular Physics, 1988