A finite‐time variational method for determining optimal paths and obtaining bounds on free energy changes from computer simulations

This paper elaborates on the method introduced by Reinhardt and Hunter [J. Chem. Phys. 97, 1599 (1992)] for determining upper and lower bounds on free energy changes and systematically improving these bounds by variational optimization of the Hamiltonian switching path connecting the two states of interest. The advantages, conditions of applicability, and potential pitfalls of this finite‐time variational method are illustrated in detail in a Monte Carlo calculation of the free energy change accompanying the reversal of magnetic field in a two‐dimensional Ising model. The power of its flexibility in an inhomogeneous system is demonstrated by computing the boundary tension between domains of positive and negative magnetization in the Ising model. Its success in a problem involving molecular dynamics and continuous potentials is illustrated by computing the free energy difference between a Lennard‐Jones fluid and a soft‐sphere fluid. Finally, its connection with the method of dynamically modified slow growt...

[1]  Frank H. Stillinger,et al.  Cluster optimization simplified by interaction modification , 1990 .

[2]  Robert H. Wood,et al.  Estimation of errors in free energy calculations due to the lag between the hamiltonian and the system configuration , 1991 .

[3]  L. Onsager.  Crystal statistics. I. A two-dimensional model with an order-disorder transition , 1944 .

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

[5]  Wolfhard Janke,et al.  Monte Carlo calculation of the surface free energy for the two-dimensional 7-state Potts model, and an estimate for four-dimensional SU(3) gauge theory , 1992, hep-lat/9202001.

[6]  Kurt Binder,et al.  Monte Carlo calculation of the surface tension for two- and three-dimensional lattice-gas models , 1982 .

[7]  Alan E. Mark,et al.  Calculation of Relative Free-Energy Via Indirect Pathways , 1991 .

[8]  D. Chandler,et al.  Introduction To Modern Statistical Mechanics , 1987 .

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

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

[11]  Potvin,et al.  Simple and efficient procedure for the Monte Carlo calculation of the surface free energy in coexisting phases. , 1989, Physical review letters.

[12]  D. Beveridge,et al.  Free energy via molecular simulation: applications to chemical and biomolecular systems. , 1989, Annual review of biophysics and biophysical chemistry.

[13]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.

[14]  Peter A. Kollman,et al.  The lag between the Hamiltonian and the system configuration in free energy perturbation calculations , 1989 .

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

[16]  Watanabe,et al.  Direct dynamical calculation of entropy and free energy by adiabatic switching. , 1990, Physical review letters.

[17]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[18]  Peter A. Kollman,et al.  A new method for carrying out free energy perturbation calculations: Dynamically modified windows , 1989 .

[19]  C. Rebbi,et al.  Surface tension and universality in the Ising model , 1993 .

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

[21]  W. Reinhardt,et al.  Variational path optimization and upper and lower bounds to free energy changes via finite time minimization of external work , 1992 .

[22]  G. Parisi,et al.  Statistical Field Theory , 1988 .

[23]  T. Straatsma,et al.  Free energy of hydrophobic hydration: A molecular dynamics study of noble gases in water , 1986 .

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

[25]  Daan Frenkel,et al.  Simulation of liquids and solids , 1987 .

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

[27]  H. Scheraga,et al.  On the multiple-minima problem in the conformational analysis of molecules: deformation of the potential energy hypersurface by the diffusion equation method , 1989 .

[28]  Stanislaw Sieniutycz,et al.  Finite-time thermodynamics and thermoeconomics , 1990 .