Reweighted Jarzynski sampling: Acceleration of rare events and free energy calculation with a bias potential learned from nonequilibrium work
arXiv:2105.03483 · doi:10.1021/acs.jctc.1c00574
Abstract
We introduce a simple enhanced sampling approach for the calculation of free energy differences and barriers along a one-dimensional reaction coordinate. First, a small number of short nonequilibrium simulations are carried out along the reaction coordinate, and the Jarzynski equality is used to learn an approximate free energy surface from the nonequilibrium work distribution. This free energy estimate is represented in a compact form as an artificial neural network and used as an external bias potential to accelerate rare events in a subsequent molecular dynamics simulation. The final free energy estimate is then obtained by reweighting the equilibrium probability distribution of the reaction coordinate sampled under the influence of the external bias. We apply our reweighted Jarzynski sampling recipe to four processes of varying scales and complexities-spanning chemical reaction in the gas phase, pair association in solution, and droplet nucleation in supersaturated vapor. In all cases, we find reweighted Jarzynski sampling to be a very efficient strategy, resulting in rapid convergence of the free energy to high precision.
References in corpus (20)
- Escaping free-energy minima
- A nonequilibrium equality for free energy differences
- CP2K: An Electronic Structure and Molecular Dynamics Software Package -- Quickstep: Efficient and Accurate Electronic Structure Calculations
- PLUMED 2: New feathers for an old bird
- Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method
- The Entropy Production Fluctuation Theorem and the Nonequilibrium Work Relation for Free Energy Differences
- Accurate sampling using Langevin dynamics
- Rethinking Metadynamics: from bias potentials to probability distributions
- A Variational Approach to Enhanced Sampling and Free Energy Calculations
- Rare events and the convergence of exponentially averaged work values
- Colored-noise thermostats à la carte
- Langevin equation with colored noise for constant-temperature molecular dynamics simulations
- Neural networks-based variationally enhanced sampling
- On practical applicability of the Jarzynski relation in statistical mechanics: a pedagogical example
- Learning Free Energy Landscapes Using Artificial Neural Networks
- Gaussian Mixture Based Enhanced Sampling
- Overcoming timescale and finite-size limitations to compute nucleation rates from small scale Well Tempered Metadynamics simulations
- Making the best of a bad situation: a multiscale approach to free energy calculation
- Time-independent free energies from metadynamics via Mean Force Integration
- Nucleation rates from small scale atomistic simulations and transition state theory