Using nonequilibrium fluctuation theorems to understand and correct errors in equilibrium and nonequilibrium discrete Langevin dynamics simulations

Abstract

Common algorithms for computationally simulating Langevin dynamics must discretize the stochastic differential equations of motion. These resulting finite time step integrators necessarily have several practical issues in common: Microscopic reversibility is violated, the sampled stationary distribution differs from the desired equilibrium distribution, and the work accumulated in nonequilibrium simulations is not directly usable in estimators based on nonequilibrium work theorems. Here, we show that even with a time-independent Hamiltonian, finite time step Langevin integrators can be thought of as a driven, nonequilibrium physical process. Once an appropriate work-like quantity is defined -- here called the shadow work -- recently developed nonequilibrium fluctuation theorems can be used to measure or correct for the errors introduced by the use of finite time steps. In particular, we demonstrate that amending estimators based on nonequilibrium work theorems to include this shadow work removes the time step dependent error from estimates of free energies. We also quantify, for the first time, the magnitude of deviations between the sampled stationary distribution and the desired equilibrium distribution for equilibrium Langevin simulations of solvated systems of varying size. While these deviations can be large, they can be eliminated altogether by Metropolization or greatly diminished by small reductions in the time step. Through this connection with driven processes, further developments in nonequilibrium fluctuation theorems can provide additional analytical tools for dealing with errors in finite time step integrators.

11 pages, 4 figures

References in corpus (10)

• Dissipation: The phase-space perspective
• Accurate sampling using Langevin dynamics
• Isothermal-isobaric molecular dynamics using stochastic velocity rescaling
• Thermodynamic metrics and optimal paths
• Rare events and the convergence of exponentially averaged work values
• Nonequilibrium candidate Monte Carlo: A new tool for efficient equilibrium simulation
• The length of time's arrow
• Unifying approach for fluctuation theorems from joint probability distributions
• Joint probability distributions and fluctuation theorems
• Estimating equilibrium ensemble averages using multiple time slices from driven nonequilibrium processes: theory and application to free energies, moments, and thermodynamic length in single-molecule pulling experiments

Cited by in corpus (18)

• Pressure control using stochastic cell rescaling
• Robust and efficient configurational molecular sampling via Langevin Dynamics
• Theory of nonequilibrium free energy transduction by molecular machines
• Optimal driving of isothermal processes close to equilibrium
• Active Brownian and inertial particles in disordered environments: short-time expansion of the mean-square displacement
• Nonlinear transport coefficients from large deviation functions
• Preserving Correlations Between Trajectories for Efficient Path Sampling
• Exactly solvable nonequilibrium Langevin relaxation of a trapped nanoparticle
• Thermodynamic inference based on coarse-grained data or noisy measurements
• Escape Time Characterization of Pendular Fabry-Perot
• Absence of dissipation in trajectory ensembles biased by currents
• Generalized equipartition for nonlinear multiplicative Langevin dynamics: application to laser-cooled atoms
• Limits on the Precision of Catenane Molecular Motors: Insights from Thermodynamics and Molecular Dynamics Simulations
• State dependent diffusion in a bistable potential: conditional probabilities and escape rates
• Time-Asymmetric Fluctuation Theorem and Efficient Free Energy Estimation
• Switching Times in Fabry-Perot Measurements
• Noise Estimate of Pendular Fabry-Perot through Reflectivity Change
• GROMACS Stochastic Dynamics and BAOAB are equivalent configurational sampling algorithms