Abstract:We introduce an efficient and robust method to compute alchemical free energy differences, resulting from the application of multiple walker Adaptive Biasing Force (ABF) in conjunction with strongly damped Langevin $\lambda$-dynamics. Unbiased alchemical free energy surfaces are naturally recovered by Thermodynamic Integration (TI). No manual optimization of the $\lambda$ schedule is required as the sampling of the $\lambda$ variable is continuous and converges towards a uniform distribution. Free diffusion of $\lambda$ improves orthogonal relaxation compared to fixed $\lambda$ methods such as standard TI or Free Energy Perturbation (FEP). Furthermore, the multiple walker strategy provides coverage of orthogonal space in a generic way with minimal user input and negligible computational overhead. Of practical importance, no adiabatic decoupling between the alchemical and Cartesian degrees of freedom is assumed, ensuring unbiased estimates for a wide envelope of numerical parameters. We present two high-performance implementations of the method in production molecular dynamics engines, namely NAMD and Tinker-HP, through coupling with the Colvars open source library. These interfaces enable the combination of the rich feature sets of those packages. We demonstrate the correctness and efficiency of the approach on several real-world cases: from solvation free energies up to ligand-receptor binding (using a recently proposed binding restraint scheme) with both fixed-charge and polarizable models. We find that, for a chosen accuracy, the computational cost is strongly reduced compared to state-of-the-art fixed-lambda methods and that results within 1~kcal/mol of experimental value are recovered for the most complex system. The implementation is publicly available and readily usable by practitioners of current alchemical methods.

Reducing the bias and uncertainty of free energy estimates by using regression to fit thermodynamic integration data

Effect of the Integration Method on the Accuracy and Computational Efficiency of Free Energy Calculations Using Thermodynamic Integration

Skewed Thermodynamic Geometry and Optimal Free Energy Estimation

Convergence Criteria for Single-step Free-energy Calculations: The Relation between the Π Bias Measure and the Sample Variance

Lambda-ABF: Simplified, Portable, Accurate and Cost-effective Alchemical Free Energy Computations

Uncertainty quantification of tabulated supercritical thermodynamics for compressible Navier-Stokes solvers

A new high order mathematical model for calculating the energy conservation equation

Thermodynamic Integration for Dynamically Unstable Systems Using Interatomic Force Constants without Molecular Dynamics

Achieving Rapid Free Energy Estimates from λ-Dynamics with Bias Updated Gibbs Sampling

Neural Thermodynamic Integration: Free Energies from Energy-based Diffusion Models

Quantified uncertainty in thermodynamic modeling for materials design

Acceleration of Solvation Free Energy Calculation via Thermodynamic Integration Coupled with Gaussian Process Regression and Improved Gelman-Rubin Convergence Diagnostics

Nonequilibrium free-energy calculation of solids using LAMMPS

Extended Regression Analysis for Debye–Einstein Models Describing Low Temperature Heat Capacity Data of Solids

Automated free energy calculation from atomistic simulations

Free-energy estimates from nonequilibrium trajectories under varying-temperature protocols

Model for Small-Sample Bias of Free-Energy Calculations Applied to Gaussian-distributed Nonequilibrium Work Measurements.

Equilibrium Equality for Free Energy Difference

ESPEI for efficient thermodynamic database development, modification, and uncertainty quantification: application to Cu-Mg

Estimating Free Energy Surfaces and their Convergence from multiple, independent static and history-dependent biased molecular-dynamics simulations with Mean Force Integration.

Computational thermodynamics and its applications