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.

On the GREM approximation of TAP free energies

On the GREM approximation of TAP free energies

On the TAP equations via the cavity approach in the generic mixed $p$-spin models

On the TAP Equations via the Cavity Approach in the Generic Mixed p-Spin Models

Contributions to Random Energy Models

Emergence of near-TAP free energy functional in the SK model at high temperature

Perturbative Interpretation of Adaptive Thouless-Anderson-Palmer Free Energy

On concavity of TAP free energy in the SK model

Dynamical Approach to the TAP Equations for the Sherrington-Kirkpatrick Model

TAP free energy, spin glasses, and variational inference

TAP variational principle for the constrained overlap multiple spherical Sherrington-Kirkpatrick model

A Constructive Proof of the Spherical Parisi Formula

Multireference Generalization of the Weighted Thermodynamic Perturbation Method.

Efficient simulation of Grassmann Tensor Product States

A simplified Parisi Ansatz II: REM universality

Free energies of Boltzmann Machines: self-averaging, annealed and replica symmetric approximations in the thermodynamic limit

Total Energy beyond GW: Exact Results and Guidelines for Approximations

On the Convexity and Reliability of the Bethe Free Energy Approximation

Targeted free energy perturbation revisited: Accurate free energies from mapped reference potentials

Assessing the accuracy and efficiency of free energy differences obtained from reweighted flow-based probabilistic generative models

Airy-averaged gradient corrections for two-dimensional fermion gases