A “fast growth” method of computing free energy differences

Let ΔF be the free energy difference between two equilibrium states of a system. An established method of numerically computing ΔF involves a single, long “switching simulation,” during which the system is driven reversibly from one state to the other (slow growth, or adiabatic switching). Here we study a method of obtaining the same result from numerous independent, irreversible simulations of a much shorter duration (fast growth). We illustrate the fast growth method, computing the excess chemical potential of a Lennard-Jones fluid as a test case, and we examine the performance of fast growth as a practical computational tool.

[1]  J. Kirkwood.  Statistical Mechanics of Fluid Mixtures , 1935 .

[2]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[3]  R. Zwanzig.  High‐Temperature Equation of State by a Perturbation Method. I. Nonpolar Gases , 1954 .

[4]  B. Widom,et al.  Some Topics in the Theory of Fluids , 1963 .

[5]  R. Reid,et al.  The Properties of Gases and Liquids , 1977 .

[6]  H. C. Andersen.  Molecular dynamics simulations at constant pressure and/or temperature , 1980 .

[7]  S. Nosé.  A unified formulation of the constant temperature molecular dynamics methods , 1984 .

[8]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

[9]  Griffiths,et al.  Chemical potential by gradual insertion of a particle in Monte Carlo simulation. , 1985, Physical review. A, General physics.

[10]  B. Efron.  The jackknife, the bootstrap, and other resampling plans , 1987 .

[11]  Jan Hermans,et al.  Simple analysis of noise and hysteresis in (slow-growth) free energy simulations , 1991 .

[12]  Robert H. Wood,et al.  Systematic errors in free energy perturbation calculations due to a finite sample of configuration space: sample-size hysteresis , 1991 .

[13]  J. Andrew McCammon,et al.  Free energy difference calculations by thermodynamic integration: Difficulties in obtaining a precise value , 1991 .

[14]  William P. Reinhardt,et al.  A finite‐time variational method for determining optimal paths and obtaining bounds on free energy changes from computer simulations , 1993 .

[15]  W. V. Gunsteren,et al.  Molecular dynamics free energy calculation in four dimensions , 1994 .

[16]  C. Mou,et al.  Absolute entropy of simple point charge model water by adiabatic switching processes , 1994 .

[17]  Berend Smit,et al.  Understanding molecular simulation: from algorithms to applications , 1996 .

[18]  C. Jarzynski.  Equilibrium free-energy differences from nonequilibrium measurements: A master-equation approach , 1997, cond-mat/9707325.

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

[20]  Ken Sekimoto,et al.  Langevin Equation and Thermodynamics , 1998 .

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

[22]  Gavin E. Crooks,et al.  Excursions in Statistical Dynamics , 1999 .

[23]  The fluctuation theorem and Green–Kubo relations , 1999, cond-mat/9902021.