Free energy reconstruction from steered dynamics without post-processing
arXiv:1001.2039 · doi:10.1016/j.jcp.2010.06.003
Abstract
Various methods achieving importance sampling in ensembles of nonequilibrium trajectories enable to estimate free energy differences and, by maximum-likelihood post-processing, to reconstruct free energy landscapes. Here, based on Bayes theorem, we propose a more direct method in which a posterior likelihood function is used both to construct the steered dynamics and to infer the contribution to equilibrium of all the sampled states. The method is implemented with two steering schedules. First, using non-autonomous steering, we calculate the migration barrier of the vacancy in Fe-alpha. Second, using an autonomous scheduling related to metadynamics and equivalent to temperature-accelerated molecular dynamics, we accurately reconstruct the two-dimensional free energy landscape of the 38-atom Lennard-Jones cluster as a function of an orientational bond-order parameter and energy, down to the solid-solid structural transition temperature of the cluster and without maximum-likelihood post-processing.
Accepted manuscript in Journal of Computational Physics, 7 figures
References in corpus (12)
- Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method
- Statistically optimal analysis of samples from multiple equilibrium states
- Isothermal-isobaric molecular dynamics using stochastic velocity rescaling
- Rare events and the convergence of exponentially averaged work values
- Path-integral analysis of fluctuation theorems for general Langevin processes
- Cumulants and large deviations of the current through non-equilibrium steady states
- Computation of free energy profiles with parallel adaptive dynamics
- Optimized free energies from bidirectional single-molecule force spectroscopy
- Comparison of work fluctuation relations
- On the efficiency of path sampling methods for the calculation of free energies from non-equilibrium simulations
- Optimum bias for fast-switching free energy calculations
- Path sampling with Stochastic Dynamics: Some new Algorithms
Cited by in corpus (9)
- Nonequilibrium candidate Monte Carlo: A new tool for efficient equilibrium simulation
- Simulating structural transitions by direct transition current sampling: the example of LJ38
- Free energy calculations from adaptive molecular dynamics simulations with adiabatic reweighting
- 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
- Waste-recycling Monte Carlo with optimal estimates: application to free energy calculations in alloys
- Using Bayes formula to estimate rates of rare events in transition path sampling simulations
- Estimating time-correlation functions by sampling and unbiasing dynamically activated events
- Using Markov transition matrices to generate trial configurations in Markov chain Monte Carlo simulations
- Sampling metastable systems using collective variables and Jarzynski-Crooks paths