Comparison of free energy methods for molecular systems

Abstract

We present a detailed comparison of computational efficiency and precision for several free energy difference ($\Delta F$) methods. The analysis includes both equilibrium and non-equilibrium approaches, and distinguishes between uni-directional and bi-directional methodologies. We are primarily interested in comparing two recently proposed approaches, adaptive integration and single-ensemble path sampling, to more established methodologies. As test cases, we study relative solvation free energies, of large changes to the size or charge of a Lennard-Jones particle in explicit water. The results show that, for the systems used in this study, both adaptive integration and path sampling offer unique advantages over the more traditional approaches. Specifically, adaptive integration is found to provide very precise long-simulation $\Delta F$ estimates as compared to other methods used in this report, while also offering rapid estimation of $\Delta F$. The results demonstrate that the adaptive integration approach is the best overall method for the systems studied here. The single-ensemble path sampling approach is found to be superior to ordinary Jarzynski averaging for the uni-directional, ``fast-growth'' non-equilibrium case. Closer examination of the path sampling approach on a two-dimensional system suggests it may be the overall method of choice when conformational sampling barriers are high. However, it appears that the free energy landscapes for the systems used in this study have rather modest configurational sampling barriers.

13 pages, 3 figures

References in corpus (7)

• A nonequilibrium equality for free energy differences
• An efficient, multiple range random walk algorithm to calculate the density of states
• Simulated Tempering: A New Monte Carlo Scheme
• Equilibrium free energy differences from nonequilibrium measurements: a master equation approach
• Path ensembles averages in systems driven far-from-equilibrium
• Theory of Systematic Computational Error in Free Energy Differences
• Single-ensemble nonequilibrium path-sampling estimates of free energy differences

Cited by in corpus (22)

• Optimal finite-time processes in stochastic thermodynamics
• Optimal protocols for minimal work processes in underdamped stochastic thermodynamics
• Generalized detailed Fluctuation Theorem under Nonequilibrium Feedback control
• A High-Throughput Steered Molecular Dynamics Study on the Free Energy Profile of Ion Permeation through Gramicidin A
• Geodesic path for the minimal energy cost in shortcuts to isothermality
• Stem-loop formation drives RNA folding in mechanical unzipping experiments
• Steady-State Homogeneous Nucleation and Growth of Water Droplets: Extended Numerical Treatment
• Comparison of free energy estimators and their dependence on dissipated work
• Modeling molecular and ionic absolute solvation free energies with quasi-chemical theory bounds
• A theoretical interpretation of variance-based convergence criteria in perturbation-based theories
• Stochastic thermodynamics of fluctuating density fields: Non-equilibrium free energy differences under coarse-graining
• Waste-recycling Monte Carlo with optimal estimates: application to free energy calculations in alloys
• Quantum Work Fluctuations in connection with Jarzynski Equality
• Predicting Solvation Free Energies in Non-Polar Solvents using Classical Density Functional Theory based on the PC-SAFT equation of state
• Jarzynski's equality, fluctuation theorems, and variance reduction: Mathematical analysis and numerical algorithms
• Satisfying the fluctuation theorem in free energy calculations with Hamiltonian replica exchange
• Reducing the Bias and Uncertainty of Free Energy Estimates by Using Regression to Fit Thermodynamic Integration Data
• Multidimensional integration through Markovian sampling under steered function morphing: a physical guise from statistical mechanics
• Free energy computations by minimization of Kullback-Leibler divergence: an efficient adaptive biasing potential method for sparse representations
• Accurate Estimation of Solvation Free Energy Using Polynomial Fitting Techniques
• Calculation of semiclassical free energy differences along non-equilibrium classical trajectories
• Equilibrium free energy differences at different temperatures from a single set of nonequilibrium transitions