Abstract:We describe the generalized weighted thermodynamic perturbation (gwTP) method for estimating the free energy surface of an expensive "high-level" potential energy function from the umbrella sampling performed with multiple inexpensive "low-level" reference potentials. The gwTP method is a generalization of the weighted thermodynamic perturbation (wTP) method developed by Li and co-workers [J. Chem. Theory Comput. 2018, 14, 5583-5596] that uses a single "low-level" reference potential. The gwTP method offers new possibilities in model design whereby the sampling generated from several low-level potentials may be combined (e.g., specific reaction parameter models that might have variable accuracy at different stages of a multistep reaction). The gwTP method is especially well suited for use with machine learning potentials (MLPs) that are trained against computationally expensive ab initio quantum mechanical/molecular mechanical (QM/MM) energies and forces using active learning procedures that naturally produce multiple distinct neural network potentials. Simulations can be performed with greater sampling using the fast MLPs and then corrected to the ab initio level using gwTP. The capabilities of the gwTP method are demonstrated by creating reference potentials based on the MNDO/d and DFTB2/MIO semiempirical models supplemented with the "range-corrected deep potential" (DPRc). The DPRc parameters are trained to ab initio QM/MM data, and the potentials are used to calculate the free energy surface of stepwise mechanisms for nonenzymatic RNA 2'--transesterification model reactions. The extended sampling made possible by the reference potentials allows one to identify unequilibrated portions of the simulations that are not always evident from the short time scale commonly used with ab initio QM/MM potentials. We show that the reference potential approach can yield more accurate ab initio free energy predictions than the wTP method or what can be reasonably afforded from explicit ab initio QM/MM sampling.

An Efficient Umbrella Potential for the Accurate Calculation of Free Energies by Molecular Simulation.

Combine Umbrella Sampling with Integrated Tempering Method for Efficient and Accurate Calculation of Free Energy Changes of Complex Energy Surface

Buckets Instead of Umbrellas for Enhanced Sampling and Free Energy Calculations

Accurate Calculations of Free-Energy Differences by the Distribution Method

Understanding Free-Energy Perturbation Calculations Through A Model of Harmonic Oscillators: Theory and Implications to Improve the Sampling Efficiency by Molecular Simulation

Tutorial on Umbrella Sampling Simulation with a Combined QM/MM Potential: The Potential of Mean Force for an SN2 Reaction in Water

Chapter 10. Free Energy Analysis Algorithms along Transition Paths and Transmembrane Ion Permeation

Weighted Least Square Analysis Method For Free Energy Calculations

Tutorial on Umbrella Sampling Simulation with a Combined QM/MM Potential: The Potential of Mean Force for an S N 2 Reaction in Water

Extension to the weighted histogram analysis method: combining umbrella sampling with free energy calculations

Do all roads lead to Rome? Convergence issues in umbrella sampling simulations

Efficient Sampling over Rough Energy Landscapes with High Barriers: A Combination of Metadynamics with Integrated Tempering Sampling.

Multireference Generalization of the Weighted Thermodynamic Perturbation Method.

Amber Free Energy Tools: Interoperable Software for Free Energy Simulations Using Generalized Quantum Mechanical/molecular Mechanical and Machine Learning Potentials.

A Transferrable Range-Separated Force Field for Water: Combining the Power of Both Physically-Motivated Models and Machine Learning Techniques

Sampling free energy surfaces as slices by combining umbrella sampling and metadynamics

Accelerated weight histogram method for exploring free energy landscapes

Do all paths lead to Rome? How reliable is umbrella sampling along a single path?

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

Eigenvector method for umbrella sampling enables error analysis

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