Reweighted Jarzynski sampling: Acceleration of rare events and free energy calculation with a bias potential learned from nonequilibrium work

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 (7)

• Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method
• Accurate sampling using Langevin dynamics
• 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
• On practical applicability of the Jarzynski relation in statistical mechanics: a pedagogical example
• Making the best of a bad situation: a multiscale approach to free energy calculation