The computation of averages from equilibrium and nonequilibrium Langevin molecular dynamics

Abstract

We consider numerical methods for thermodynamic sampling, i.e. computing sequences of points distributed according to the Gibbs-Boltzmann distribution, using Langevin dynamics and overdamped Langevin dynamics (Brownian dynamics). A wide variety of numerical methods for Langevin dynamics may be constructed based on splitting the stochastic differential equations into various component parts, each of which may be propagated exactly in the sense of distributions. Each such method may be viewed as generating samples according to an associated invariant measure that differs from the exact canonical invariant measure by a stepsize-dependent perturbation. We provide error estimates a la Talay-Tubaro on the invariant distribution for small stepsize, and compare the sampling bias obtained for various choices of splitting method. We further investigate the overdamped limit and apply the methods in the context of driven systems where the goal is sampling with respect to a nonequilibrium steady state. Our analyses are illustrated by numerical experiments.

Accepted to IMA Journal of Numerical Analysis

References in corpus (7)

• Accurate sampling using Langevin dynamics
• Robust and efficient configurational molecular sampling via Langevin Dynamics
• Design of quasi-symplectic propagators for Langevin dynamics
• From ballistic to diffusive behavior in periodic potentials
• Trotter Derivation of Algorithms for Brownian and Dissipative Particle Dynamics
• Langevin dynamics with space-time periodic nonequilibrium forcing
• Weak backward error analysis for Langevin process

Cited by in corpus (55)

• A simple and accurate algorithm for path integral molecular dynamics with the Langevin thermostat
• On the numerical treatment of dissipative particle dynamics and related systems
• Global Convergence of Stochastic Gradient Hamiltonian Monte Carlo for Non-Convex Stochastic Optimization: Non-Asymptotic Performance Bounds and Momentum-Based Acceleration
• Local and Global Perspectives on Diffusion Maps in the Analysis of Molecular Systems
• Stationary state distribution and efficiency analysis of the Langevin equation via real or virtual dynamics
• A fast spectral method for electrostatics in doubly-periodic slit channels
• Structure-preserving integrators for dissipative systems based on reversible-irreversible splitting
• Decentralized Stochastic Gradient Langevin Dynamics and Hamiltonian Monte Carlo
• Pairwise adaptive thermostats for improved accuracy and stability in dissipative particle dynamics
• Convergence of unadjusted Hamiltonian Monte Carlo for mean-field models
• On the long-time integration of stochastic gradient systems
• Order conditions for sampling the invariant measure of ergodic stochastic differential equations on manifolds
• Langevin dynamics with space-time periodic nonequilibrium forcing
• Understanding Molecular Dynamics with Stochastic Processes via Real or Virtual Dynamics
• Nonreversible Langevin Samplers: Splitting Schemes, Analysis and Implementation
• Exotic aromatic B-series for the study of long time integrators for a class of ergodic SDEs
• Covariance-Controlled Adaptive Langevin Thermostat for Large-Scale Bayesian Sampling
• Stable schemes for dissipative particle dynamics with conserved energy
• Accurate and robust splitting methods for the generalized Langevin equation with a positive Prony series memory kernel
• New parallelizable schemes for integrating the Dissipative Particle Dynamics with Energy Conservation
• Local density dependent potential for compressible mesoparticles
• Constraint-Based Regularization of Neural Networks
• Estimating the speed-up of Adaptively Restrained Langevin Dynamics
• Size consistency in Smoothed Dissipative Particle Dynamics
• The strong convergence of operator-splitting methods for the Langevin dynamics model
• Langevin Markov Chain Monte Carlo with stochastic gradients
• Uniform minorization condition and convergence bounds for discretizations of kinetic Langevin dynamics
• Velocity Jumps for Molecular Dynamics
• Ergodic SDEs on submanifolds and related numerical sampling schemes
• Constructing sampling schemes via coupling: Markov semigroups and optimal transport
• Theoretical analysis and simulation methods for Hawkes processes and their diffusion approximation
• Velocity jump processes : an alternative to multi-timestep methods for faster and accurate molecular dynamics simulations
• Convergence of Stochastic-extended Lagrangian molecular dynamics method for polarizable force field simulation
• A uniformly accurate scheme for the numerical integration of penalized Langevin dynamics
• The shifted ODE method for underdamped Langevin MCMC
• On the Numerical Stationary Distribution of Overdamped Langevin Equation in Harmonic System
• A Unifying and Canonical Description of Measure-Preserving Diffusions
• Cayley Splitting for Second-Order Langevin Stochastic Partial Differential Equations
• Information Criteria for quantifying loss of reversibility in parallelized KMC
• A technique for studying strong and weak local errors of splitting stochastic integrators
• Ensemble preconditioning for Markov chain Monte Carlo simulation
• Information metrics for long-time errors in splitting schemes for stochastic dynamics and parallel KMC
• Some new results on relative entropy production, time reversal, and optimal control of time-inhomogeneous diffusion processes
• Estimating linear response statistics using orthogonal polynomials: An RKHS formulation
• Piecewise deterministic sampling with splitting schemes
• Convergence rates for nonequilibrium Langevin dynamics
• High order conformal symplectic and ergodic schemes for stochastic Langevin equation via generating functions
• Error estimates of the backward Euler-Maruyama method for multi-valued stochastic differential equations
• Multiscale Gentlest Ascent Dynamics for Saddle Point in Effective Dynamics of Slow-Fast System
• TATi-Thermodynamic Analytics ToolkIt: TensorFlow-based software for posterior sampling in machine learning applications
• Stochastic Norton dynamics: An alternative approach for the computation of transport coefficients in dissipative particle dynamics
• Quadrature Points via Heat Kernel Repulsion
• Convergence of the likelihood ratio method for linear response of non-equilibrium stationary states
• GROMACS Stochastic Dynamics and BAOAB are equivalent configurational sampling algorithms
• Ergodic properties of quasi-Markovian generalized Langevin equations with configuration dependent noise and non-conservative force