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

David A. Sivak, John D. Chodera, and Gavin E. Crooks.
Phys. Rev. X 3:011007, 2013. [DOI] [PDF]

The finite-timestep errors in molecular dynamics simulations can be interpreted as a form of nonequilibrium work.  We show how this leads to straightforward schemes for correcting for these errors or assessing their impact.

Keywords: velocity verlet with Velocity randomization; VVVR; nonequilibrium free energy; integrator error; nonequilibrium integration

Nonequilibrium candidate Monte Carlo is an efficient tool for equilibrium simulation

Jerome P. Nilmeier, Gavin E. Crooks, David D. L. Minh, and John D. Chodera. 
Proc. Natl. Acad. Sci. USA 108:E1009, 2011. [DOI] [PDF]

We present a significant generalization of Monte Carlo methods that provide an enormously useful tool for enhancing the efficiency of molecular simulations and enabling molecular design.

Keywords: NCMC; Monte Carlo; Metropolis-Hastings; acceptance rates; molecular dynamics

Estimating equilibrium ensemble averages using multiple time slices from driven nonequilibrium processes

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
David D. L. Minh and John D. Chodera
J. Chem. Phys. 134:024111, 2011. [DOI] [PDF]

We derive a new estimator for estimating equilibrium expectations from nonequilibrium experiments, and show how it can be used to estimate a variety of useful quantities in simulated single-molecule force spectroscopy experiments.

Optimal estimators and asymptotic variances for nonequilibrium path-ensemble averages

David D. D. L. Minh and John D. Chodera.
J. Chem. Phys. 131:134110, 2009. [DOI] [PDF]

We derive an optimal estimator and corresponding statistical uncertainties for inferring expectations of bidirectional nonequilibrium processes.  These estimators have widespread applicability in single-molecule biophysical force-spectroscopy experiments and nonequilibrium molecular simulations.