Free energy reconstruction from steered dynamics without post-processing

Various methods achieving importance sampling in ensembles of nonequilibrium trajectories enable one 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 [email protected] 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.

[1]  C. Jarzynski.  Nonequilibrium Equality for Free Energy Differences , 1996, cond-mat/9610209.

[2]  C. Dellago,et al.  Transition Path Sampling , 2005 .

[3]  Graeme Ackland,et al.  Development of an interatomic potential for phosphorus impurities in α-iron , 2004 .

[4]  Gabriel Stoltz,et al.  Computation of free energy profiles with parallel adaptive dynamics. , 2007, The Journal of chemical physics.

[5]  M. Parrinello,et al.  Isothermal-isobaric molecular dynamics using stochastic velocity rescaling. , 2009, The Journal of chemical physics.

[6]  D. Landau,et al.  Efficient, multiple-range random walk algorithm to calculate the density of states. , 2000, Physical review letters.

[7]  G. Ciccotti,et al.  Mapping the network of pathways of CO diffusion in myoglobin. , 2010, Journal of the American Chemical Society.

[8]  Pascal Bellon,et al.  Identification of novel diffusion cycles in B2 ordered phases by Monte Carlo simulation , 1997 .

[9]  D. Ceperley,et al.  Monte Carlo simulation of a many-fermion study , 1977 .

[10]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[11]  E. Vanden-Eijnden,et al.  A temperature accelerated method for sampling free energy and determining reaction pathways in rare events simulations , 2006 .

[12]  Mihai-Cosmin Marinica,et al.  Orientation of Interstitials in Clusters in α-Fe: A Comparison between Empirical Potentials , 2007 .

[13]  A. Laio,et al.  Escaping free-energy minima , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[14]  M. Kalos,et al.  First-passage Monte Carlo algorithm: diffusion without all the hops. , 2006, Physical review letters.

[15]  B. Jourdain,et al.  DOES WASTE-RECYCLING REALLY IMPROVE THE MULTI-PROPOSAL METROPOLIS-HASTINGS MONTE CARLO ALGORITHM ? , 2009 .

[16]  D. Frenkel.  Speed-up of Monte Carlo simulations by sampling of rejected states. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[17]  C. Jarzynski.  Nonequilibrium work relations : foundations and applications , 2008 .

[18]  M. Athènes.  A path-sampling scheme for computing thermodynamic properties of a many-body system in a generalized ensemble , 2004 .

[19]  C. Jarzynski,et al.  Path-integral analysis of fluctuation theorems for general Langevin processes , 2006, cond-mat/0605471.

[20]  Christopher Jarzynski,et al.  Comparison of work fluctuation relations , 2007, 0710.2261.

[21]  D. L. Freeman,et al.  Phase changes in 38-atom Lennard-Jones clusters. I. A parallel tempering study in the canonical ensemble , 2000, physics/0003068.

[22]  Daan Frenkel,et al.  Virtual-move parallel tempering. , 2005, Chemphyschem : a European journal of chemical physics and physical chemistry.

[23]  Manuel Athènes,et al.  Thermodynamic modelling of glasses at atomistic scale , 2008 .

[24]  David J Wales,et al.  Potential energy and free energy landscapes. , 2006, The journal of physical chemistry. B.

[25]  Michael R. Shirts,et al.  Equilibrium free energies from nonequilibrium measurements using maximum-likelihood methods. , 2003, Physical review letters.

[26]  F. Calvo,et al.  Sampling along reaction coordinates with the Wang-Landau method , 2002, cond-mat/0205428.

[27]  M. Athènes,et al.  Web ensemble averages for retrieving relevant information from rejected Monte Carlo moves , 2007 .

[28]  Gabriel Stoltz.  Path sampling with stochastic dynamics: Some new algorithms , 2007, J. Comput. Phys..

[29]  Christoph Dellago,et al.  Optimum bias for fast-switching free energy calculations , 2007, Comput. Phys. Commun..

[30]  D. Frenkel.  Waste-Recycling Monte Carlo , 2006 .

[31]  G. Torrie,et al.  Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling , 1977 .

[32]  Michael R. Shirts,et al.  Statistically optimal analysis of samples from multiple equilibrium states. , 2008, The Journal of chemical physics.

[33]  G. Adjanor,et al.  Free energy landscape from path-sampling: application to the structural transition in LJ38 , 2006 .

[34]  Berend Smit,et al.  Understanding Molecular Simulation , 2001 .

[35]  Y. Kuzovlev,et al.  Nonlinear fluctuation-dissipation relations and stochastic models in nonequilibrium thermodynamics , 1981 .

[36]  G. Crooks.  Nonequilibrium Measurements of Free Energy Differences for Microscopically Reversible Markovian Systems , 1998 .

[37]  G. Crooks.  Path-ensemble averages in systems driven far from equilibrium , 1999, cond-mat/9908420.

[38]  Charles H. Bennett,et al.  Efficient estimation of free energy differences from Monte Carlo data , 1976 .

[39]  C. Jarzynskia.  Nonequilibrium work relations: foundations and applications , 2008 .

[40]  M Athènes.  Parallel Monte Carlo simulations using a residence weight algorithm. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[41]  M. Parrinello,et al.  Well-tempered metadynamics: a smoothly converging and tunable free-energy method. , 2008, Physical review letters.

[42]  P. Steinhardt,et al.  Bond-orientational order in liquids and glasses , 1983 .

[43]  Flynn,et al.  Jump dynamics and the isotope effect in solid-state diffusion. , 1986, Physical review letters.

[44]  Manuel Athènes,et al.  Multiple-replica exchange with information retrieval. , 2008, Chemphyschem : a European journal of chemical physics and physical chemistry.

[45]  Manuel Athènes,et al.  Measurement of nonequilibrium entropy from space-time thermodynamic integration. , 2008, The Journal of chemical physics.

[46]  G. Hummer,et al.  Free energy reconstruction from nonequilibrium single-molecule pulling experiments , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[47]  M. Athènes.  Computation of a chemical potential using a residence weight algorithm. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[48]  D. R. Mason,et al.  Stochastic kinetic Monte Carlo algorithms for long-range Hamiltonians , 2004, Comput. Phys. Commun..

[49]  D. L. Freeman,et al.  Phase changes in 38-atom Lennard-Jones clusters. II. A parallel tempering study of equilibrium and dynamic properties in the molecular dynamics and microcanonical ensembles , 2000, physics/0003072.

[50]  David D L Minh,et al.  Optimized free energies from bidirectional single-molecule force spectroscopy. , 2008, Physical review letters.

[51]  M. Tuckerman,et al.  On the use of the adiabatic molecular dynamics technique in the calculation of free energy profiles , 2002 .

[52]  Benjamin Jourdain,et al.  Does Waste Recycling Really Improve the Multi-Proposal Metropolis–Hastings algorithm? an Analysis Based on Control Variates , 2009, Journal of Applied Probability.

[53]  C. Dellago,et al.  Biased sampling of nonequilibrium trajectories: can fast switching simulations outperform conventional free energy calculation methods? , 2005, The journal of physical chemistry. B.

[54]  C. Jarzynski.  Rare events and the convergence of exponentially averaged work values. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[55]  Daan Frenkel,et al.  Monte Carlo Sampling of a Markov Web. , 2005, Journal of chemical theory and computation.

[56]  Alan M. Ferrenberg,et al.  Optimized Monte Carlo data analysis. , 1989, Physical Review Letters.

[57]  J. M. Lanore.  Simulation of the point defect behaviour in a lattice by the Monte-Carlo method , 1974 .

[58]  Christopher Jarzynski,et al.  Nonequilibrium work theorem for a system strongly coupled to a thermal environment , 2004 .

[59]  Michele Parrinello,et al.  Molecular dynamics study of the solvation of calcium carbonate in water. , 2007, The journal of physical chemistry. B.

[60]  Jeanne-Marie Lanore,et al.  Simulation de l'evolution des defauts dans un reseau par le methode de monte-carlo , 1974 .

[61]  T. Lelièvre,et al.  Free Energy Computations: A Mathematical Perspective , 2010 .

[62]  J. VandeVondele,et al.  Canonical Adiabatic Free Energy Sampling (CAFES): A Novel Method for the Exploration of Free Energy Surfaces , 2002 .

[63]  E. Vanden-Eijnden,et al.  Large-scale conformational sampling of proteins using temperature-accelerated molecular dynamics , 2010, Proceedings of the National Academy of Sciences.

[64]  B. Derrida,et al.  Cumulants and large deviations of the current through non-equilibrium steady states , 2007, 0704.2726.

[65]  Eric Vanden-Eijnden,et al.  Some recent techniques for free energy calculations , 2009, J. Comput. Chem..

[66]  J. Doye,et al.  THE DOUBLE-FUNNEL ENERGY LANDSCAPE OF THE 38-ATOM LENNARD-JONES CLUSTER , 1998, cond-mat/9808265.

[67]  David J Wales,et al.  Calculating rate constants and committor probabilities for transition networks by graph transformation. , 2009, The Journal of chemical physics.

[68]  M. Athènes,et al.  Gibbs free-energy estimates from direct path-sampling computations. , 2005, The Journal of chemical physics.

[69]  A. B. Bortz,et al.  A new algorithm for Monte Carlo simulation of Ising spin systems , 1975 .

[70]  Jorge Kurchan,et al.  Fluctuation theorem for stochastic dynamics , 1998 .

[71]  Esselink,et al.  Parallel Monte Carlo simulations. , 1995, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[72]  Christoph Dellago,et al.  On the efficiency of path sampling methods for the calculation of free energies from non-equilibrium simulations , 2007 .