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.

Efficient treatment of local meta-generalized gradient density functionals via auxiliary density expansion: the density fitting (DF) J+X approximation

Meta-generalized gradient approximation: explanation of a realistic nonempirical density functional

Partially deorbitalized meta-GGA

Revisiting Gauge-Independent Kinetic Energy Densities in Meta-GGAs and Local Hybrid Calculations of Magnetizabilities

Towards accurate orbital-free simulations: A generalized gradient approximation for the noninteracting free energy density functional

A simple generalized gradient approximation for the noninteracting kinetic energy density functional

Atoms, molecules, solids, and surfaces: Applications of the generalized gradient approximation for exchange and correlation.

Framework for Laplacian-Level Noninteracting Free-Energy Density Functionals

Generalized Gradient Approximation Made Simple.

A More Accurate Generalized Gradient Approximation for Solids

Construction of a generalized gradient approximation by restoring the density-gradient expansion and enforcing a tight Lieb-Oxford bound

Regularized Gradient Expansion for Atoms, Molecules, and Solids

Prescription for the design and selection of density functional approximations: more constraint satisfaction with fewer fits

Meta-Generalized-Gradient Approximation made Magnetic

A Hybrid Meta On-Top Functional for Multiconfiguration Pair-Density Functional Theory

Multireference Generalization of the Weighted Thermodynamic Perturbation Method.

Doubly Hybrid Density Functional for Accurate Descriptions of Nonbond Interactions, Thermochemistry, and Thermochemical Kinetics

Analytic Gradients for the Long-Range-Corrected XYG3 Type of Doubly Hybrid Density Functionals: Theory, Implementation, and Assessment.

A General and Transferable Local Hybrid Functional for Electronic Structure Theory and Many-Fermion Approaches

On the challenge to improve the density response with unusual gradient approximations