Accurate multiple time step in biased molecular simulations

Abstract

Many recently introduced enhanced sampling techniques are based on biasing coarse descriptors (collective variables) of a molecular system on the fly. Sometimes the calculation of such collective variables is expensive and becomes a bottleneck in molecular dynamics simulations. An algorithm to treat smooth biasing forces within a multiple time step framework is here discussed. The implementation is simple and allows a speed up when expensive collective variables are employed. The gain can be substantial when using massively parallel or GPU-based molecular dynamics software. Moreover, a theoretical framework to assess the sampling accuracy is introduced, which can be used to assess the choice of the integration time step in both single and multiple time step biased simulations.

Accepted for publication on J. Chem. Theory Comput

References in corpus (6)

• Canonical sampling through velocity-rescaling
• Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method
• PLUMED: a portable plugin for free-energy calculations with molecular dynamics
• Accurate sampling using Langevin dynamics
• Efficient multiple time scale molecular dynamics: using colored noise thermostats to stabilize resonances
• RNA unwinding from reweighted pulling simulations

Cited by in corpus (8)

• Combining simulations and solution experiments as a paradigm for RNA force field refinement
• Enhanced Conformational Sampling using Replica Exchange with Collective-Variable Tempering
• Analyzing and biasing simulations with PLUMED
• Accurate molecular dynamics and nuclear quantum effects at low cost by multiple steps in real and imaginary time: using density functional theory to accelerate wavefunction methods
• Reweighting of molecular simulations with explicit-solvent SAXS restraints elucidates ion-dependent RNA ensembles
• ATP dependent NS3 helicase interaction with RNA: insights from molecular simulations
• Reaction Coordinates for Conformational Transitions using Linear Discriminant Analysis on Positions
• Extending and validating bubble nucleation rate predictions in a Lennard-Jones fluid with enhanced sampling methods and transition state theory