Abstract:We present methods that introduce concepts from Rosenbluth sampling [M. N. Rosenbluth and A. W. Rosenbluth, J. Chem. Phys. 23, 356 (1955)] into the Jarzynski nonequilibrium work (NEW) free-energy calculation technique [C. Jarzynski, Phys. Rev. Lett. 78, 2690 (1997)]. The proposed hybrid modifies the way steps are taken in the NEW process. With it, each step is selected from a range of alternatives, with bias given to steps that contribute the least work. The definition of the work average is modified to account for the bias. We introduce two variants of this method, l-bias sampling and configuration-bias sampling, respectively; a combined l-and configuration-bias method is also considered. By reducing the likelihood that large nonequilibrated work values enter the ensemble average, the Rosenbluth sampling aids in remedying problems of inaccuracy of the calculation. We demonstrate the performance of the proposed methods through a model system of N independent harmonic oscillators. This model captures the difficulties involved in calculating free energies in real systems while retaining many tractable features that are helpful to the study. We examine four variants of this model that differ qualitatively in the nature of their phase-space overlap. Results indicate that the l-bias sampling method is most useful for systems with entropic sampling barriers, while the configuration-bias methods are best for systems with energetic sampling barriers. The Rosenbluth-sampling schemes yield much more accurate results than the unbiased nonequilibrium work method. Typically the accuracy can be improved by about an order of magnitude for a given amount of sampling; this improvement translates into two or more orders of magnitude less sampling required to obtain a given level of accuracy, owing to the generally slow convergence of the NEW calculation when the inaccuracy is large. 2005 American Institute of Physics.

On the Calculation of Absolute Free Energies from Molecular-Dynamics or Monte Carlo Data.

A method to compute absolute free energies or enthalpies of fluids

A novel method for calculating relative free energy of similar molecules in two environments

Accurate Calculations of Free-Energy Differences by the Distribution Method

Efficient Calculation of Free Energy Differences Associated with Isotopic Substitution Using Path-Integral Molecular Dynamics

Computing absolute free energies of disordered structures by molecular simulation

An Efficient Method for Computing Excess Free Energy of Liquid

Optimal Estimation of Free Energies and Stationary Densities from Multiple Biased Simulations

Nonequilibrium free-energy calculation of solids using LAMMPS

Configurational Space Discretization and Free Energy Calculation in Complex Molecular Systems

A New Monte Carlo Algorithm for Free Energy Calculation

Direct free energy calculation from ab initio path integral Monte Carlo simulations of warm dense matter

How to measure the free energy and partition function from atom-atom correlations

η-ensemble path integral Monte Carlo approach to the free energy of the warm dense electron gas and the uniform electron liquid

Free Energy Methods for the Description of Molecular Processes

Calculations Of Free Energies Based On The Combination Of Molecular Dynamics Simulations And Continuum Solvation Model

Rosenbluth-Sampled Nonequilibrium Work Method For Calculation Of Free Energies In Molecular Simulation

Accelerated weight histogram method for exploring free energy landscapes

Efficient Approaches To Solutions Of Partition Function For Condensed Matters

Comparison of Two Efficient Methods for Calculating Partition Functions

Using Overlap and Funnel Sampling to Obtain Accurate Free Energies from Nonequilibrium Work Measurements.