Path probability ratios for Langevin dynamics -- exact and approximate

Abstract

Path reweighting is a principally exact method to estimate dynamic properties from biased simulations - provided that the path probability ratio matches the stochastic integrator used in the simulation. Previously reported path probability ratios match the Euler-Maruyama scheme for overdamped Langevin dynamics. Since MD simulations use Langevin dynamics rather than overdamped Langevin dynamics, this severely impedes the application of path reweighting methods. Here, we derive the path probability ratio $M_L$ for Langevin dynamics propagated by a variant of the Langevin Leapfrog integrator. This new path probability ratio allows for exact reweighting of Langevin dynamics propagated by this integrator. We also show that a previously derived approximate path probability ratio $M_{\mathrm{approx}}$ differs from the exact $M_L$ only by $\mathcal{O}({\xi }^4\Delta t^4)$, and thus yields highly accurate dynamic reweighting results. ($\Delta t$ is the integration time step, $\xi$ is the collision rate.) The results are tested and the efficiency of path-reweighting is explored using butane as an example.

References in corpus (17)

• Identification of slow molecular order parameters for Markov model construction
• From Metadynamics to Dynamics
• Accurate sampling using Langevin dynamics
• Best Practices for Alchemical Free Energy Calculations
• Robust and efficient configurational molecular sampling via Langevin Dynamics
• Langevin equation with colored noise for constant-temperature molecular dynamics simulations
• Maximum Caliber: a general variational principle for dynamical systems
• Statistically optimal analysis of state-discretized trajectory data from multiple thermodynamic states
• Efficient Dynamic Importance Sampling of Rare Events in One Dimension
• Time step rescaling recovers continuous-time dynamical properties for discrete-time Langevin integration of nonequilibrium systems
• Stochastic actions for diffusive dynamics: Reweighting, sampling, and minimization
• Girsanov reweighting for path ensembles and Markov state models
• Dynamical reweighting methods for Markov models
• xTRAM: Estimating equilibrium expectations from time-correlated simulation data at multiple thermodynamic states
• Estimation of the infinitesimal generator by square-root approximation
• Path sampling with Stochastic Dynamics: Some new Algorithms
• Microscopic reweighting for non-equilibrium steady states dynamics

Cited by in corpus (5)

• A review of Girsanov Reweighting and of Square Root Approximation for building molecular Markov State Models
• Assessing transition rates as functions of environmental variables
• Inferring nonlinear fractional diffusion processes from single trajectories
• Predicting the unobserved: a statistical mechanics framework for non-equilibrium material response with quantified uncertainty
• GROMACS Stochastic Dynamics and BAOAB are equivalent configurational sampling algorithms