Robust and efficient configurational molecular sampling via Langevin Dynamics

Abstract

A wide variety of numerical methods are evaluated and compared for solving the stochastic differential equations encountered in molecular dynamics. The methods are based on the application of deterministic impulses, drifts, and Brownian motions in some combination. The Baker-Campbell-Hausdorff expansion is used to study sampling accuracy following recent work by the authors, which allows determination of the stepsize-dependent bias in configurational averaging. For harmonic oscillators, configurational averaging is exact for certain schemes, which may result in improved performance in the modelling of biomolecules where bond stretches play a prominent role. For general systems, an optimal method can be identified that has very low bias compared to alternatives. In simulations of the alanine dipeptide reported here (both solvated and unsolvated), higher accuracy is obtained without loss of computational efficiency, while allowing large timestep, and with no impairment of the conformational exploration rate (the effective diffusion rate observed in simulation). The optimal scheme is a uniformly better performing algorithm for molecular sampling, with overall efficiency improvements of 25% or more in practical timestep size achievable in vacuum, and with reductions in the error of configurational averages of a factor of ten or more attainable in solvated simulations at large timestep.

The article has been submitted to the Journal of Chemical Physics

References in corpus (5)

• Accurate sampling using Langevin dynamics
• Robust and efficient configurational molecular sampling via Langevin Dynamics
• Design of quasi-symplectic propagators for Langevin dynamics
• Trotter Derivation of Algorithms for Brownian and Dissipative Particle Dynamics
• Discretization errors in molecular dynamics simulations with deterministic and stochastic thermostats

Cited by in corpus (70)

• GPUMD: A package for constructing accurate machine-learned potentials and performing highly efficient atomistic simulations
• Robust and efficient configurational molecular sampling via Langevin Dynamics
• A simple and accurate algorithm for path integral molecular dynamics with the Langevin thermostat
• On the Random Batch Method for second order interacting particle systems
• The computation of averages from equilibrium and nonequilibrium Langevin molecular dynamics
• A versatile framework for simulating the dynamic mechanical structure of cytoskeletal networks
• A unified thermostat scheme for efficient configurational sampling for classical/quantum canonical ensembles via molecular dynamics
• On the numerical treatment of dissipative particle dynamics and related systems
• The Generalized Langevin Equation: An efficient approach to non-equilibrium molecular dynamics of open systems
• i-PI 3.0: a flexible and efficient framework for advanced atomistic simulations
• Graph Neural Networks Accelerated Molecular Dynamics
• Cayley modification for strongly stable path-integral and ring-polymer molecular dynamics
• Numerical integrators for the Hybrid Monte Carlo method
• Scalable Hybrid Deep Neural Networks/Polarizable Potentials Biomolecular Simulations including long-range effects
• Design principles for selective self-assembly of active networks
• Applications of the Generalised Langevin Equation: towards a realistic description of the baths
• Dynamic density functional theories for inhomogeneous polymer systems compared to Brownian dynamics simulations
• Probing quantum coherence in ultrafast molecular processes: an ab initio approach to open quantum systems
• Modeling the structural and thermal properties of loaded metal-organic frameworks. An interplay of quantum and anharmonic fluctuations
• Finite-temperature materials modeling from the quantum nuclei to the hot electrons regime
• Local and Global Perspectives on Diffusion Maps in the Analysis of Molecular Systems
• Adaptive Thermostats for Noisy Gradient Systems
• Stationary state distribution and efficiency analysis of the Langevin equation via real or virtual dynamics
• Assessing mixed quantum-classical molecular dynamics methods for nonadiabatic dynamics of molecules on metal surfaces
• Path integral molecular dynamics with surface hopping for thermal equilibrium sampling of nonadiabatic systems
• A random batch Ewald method for charged particles in the isothermal-isobaric ensemble
• Sampling the isothermal-isobaric ensemble by Langevin dynamics
• Structure-preserving integrators for dissipative systems based on reversible-irreversible splitting
• Pairwise adaptive thermostats for improved accuracy and stability in dissipative particle dynamics
• Molecular dynamics simulations in hybrid particle-continuum schemes: Pitfalls and caveats
• Path probability ratios for Langevin dynamics -- exact and approximate
• Plasmonic enhancement of molecular hydrogen dissociation on metallic magnesium nanoclusters
• Assessing numerical methods for molecular and particle simulation
• Nonequilibrium processes from Generalised Langevin Equations: realistic nanoscale systems connected to two thermal baths
• Systematic derivation of hybrid coarse-grained models
• Understanding Molecular Dynamics with Stochastic Processes via Real or Virtual Dynamics
• Comparison of Modern Langevin Integrators for Simulations of Coarse-Grained Polymer Melts
• Fast-Forward Langevin Dynamics with Momentum Flips
• Accurate and robust splitting methods for the generalized Langevin equation with a positive Prony series memory kernel
• Improved torque estimator for condensed-phase quasicentroid molecular dynamics
• Contraction and Convergence Rates for Discretized Kinetic Langevin Dynamics
• MIST: A Simple and Efficient Molecular Dynamics Abstraction Library for Integrator Development
• Optimal friction matrix for underdamped Langevin sampling
• Time correlation functions of equilibrium and nonequilibrium Langevin dynamics: Derivations and numerics using random numbers
• Accelerated sampling by infinite swapping of path integral molecular dynamics with surface hopping
• Tunable Spatiotemporal Orders in Driven Insulators
• The strong convergence of operator-splitting methods for the Langevin dynamics model
• Dynamic principle for ensemble control tools
• Efficient Sampling of Thermal Averages of Interacting Quantum Particle Systems with Random Batches
• Data clustering based on Langevin annealing with a self-consistent potential
• Continuum limit and preconditioned Langevin sampling of the path integral molecular dynamics
• Accurate and efficient splitting methods for dissipative particle dynamics
• Harmonic Oscillator Staging Coordinates for Efficient Path Integral Simulations of Quantum Oscillators and Crystals
• Velocity jump processes : an alternative to multi-timestep methods for faster and accurate molecular dynamics simulations
• The Bayesian Inversion Problem for Thermal Average Sampling of Quantum Systems
• Comparison of effective and stable Langevin dynamics integrators
• On the relation between the velocity- and position-Verlet integrators
• Velocity Jumps for Molecular Dynamics
• Non-unique Hamiltonians for Discrete Symplectic Dynamics
• Robust Kirkwood-Buff inversion in complex mixtures via reciprocal-space methods
• Understanding the Sources of Error in MBAR through Asymptotic Analysis
• Modelling a Bistable System Strongly Coupled to a Debye Bath: A Quasiclassical Approach Based on the Generalised Langevin Equation
• On the Numerical Stationary Distribution of Overdamped Langevin Equation in Harmonic System
• Dynamics with Simultaneous Dissipations to Fermionic and Bosonic Reservoirs
• Variational Langevin Hamiltonian Monte Carlo for Distant Multi-modal Sampling
• Random-batch list algorithm for short-range molecular dynamics simulations
• GROMACS Stochastic Dynamics and BAOAB are equivalent configurational sampling algorithms
• Faster Molecular Dynamics with Neural Network Potentials via Distilled Multiple Time-Stepping and Non-Conservative Forces
• Multi-level Monte Carlo path integral molecular dynamics for thermal average calculation in the nonadiabatic regime
• NQCDynamics.jl: A Julia Package for Nonadiabatic Quantum Classical Molecular Dynamics in the Condensed Phase