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

Abstract

Molecular dynamics is one of the most commonly used approaches for studying the dynamics and statistical distributions of many 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. 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. In this paper, we introduce a set of stochastic isokinetic equations of motion that are shown to be rigorously ergodic and that can be integrated using a multiple time-stepping algorithm that can be 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.

References in corpus (2)

• Thermostats for "slow" configurational modes
• Efficient multiple time scale molecular dynamics: using colored noise thermostats to stabilize resonances

Cited by in corpus (17)

• A unified thermostat scheme for efficient configurational sampling for classical/quantum canonical ensembles via molecular dynamics
• Accurate multiple time step in biased molecular simulations
• Adaptive Thermostats for Noisy Gradient Systems
• Adaptive multi-stage integrators for optimal energy conservation in molecular simulations
• Speeding-up Ab Initio Molecular Dynamics with Hybrid Functionals using Adaptively Compressed Exchange Operator based Multiple Timestepping
• From Classical to Quantum and Back: Hamiltonian Adaptive Resolution Path Integral, Ring Polymer, and Centroid Molecular Dynamics
• Wavefunction-based electrostatic-embedding QM/MM using CFOUR through MiMiC
• Enhanced solvation force extrapolation for speeding up molecular dynamics simulations of complex biochemical liquids
• Speeding-up Hybrid Functional based Ab Initio Molecular Dynamics using Multiple Time-stepping and Resonance Free Thermostat
• Compressing and forecasting atomic material simulations with descriptors
• Velocity jump processes : an alternative to multi-timestep methods for faster and accurate molecular dynamics simulations
• The Effects of Intrinsic Dynamical Ghost Modes in Discrete-Time Langevin Simulations
• Velocity Jumps for Molecular Dynamics
• On ergodicity for multi-dimensional harmonic oscillator systems with Nose-Hoover type thermostat
• Piecewise deterministic sampling with splitting schemes
• Faster Molecular Dynamics with Neural Network Potentials via Distilled Multiple Time-Stepping and Non-Conservative Forces
• GROMACS Stochastic Dynamics and BAOAB are equivalent configurational sampling algorithms