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.

Modeling Many-Body Interactions in Water with Gaussian Process Regression

Gaussian Mixture Models for Diatomic Gas−surface Interactions under Thermal Non-Equilibrium Conditions

FFLUX molecular simulations driven by atomic Gaussian process regression models

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

Modeling Molecular Interactions in Water: From Pairwise to Many-Body Potential Energy Functions

Toward Gaussian Process Regression Modeling of a Urea Force Field

Comparison of permutationally invariant polynomials, neural networks, and Gaussian approximation potentials in representing water interactions through many-body expansions

Completely Multipolar Model as a General Framework for Many-Body Interactions as Illustrated for Water

Electronic coarse graining enhances the predictive power of molecular simulation allowing challenges in water physics to be addressed

Multireference Generalization of the Weighted Thermodynamic Perturbation Method.

Simple and Accurate One-Body Energy and Dipole Moment Surfaces for Water and Beyond

Seamless integration of GEM, a density based-force field, for QM/MM simulations via LICHEM, Psi4 and Tinker-HP

Combining Force Fields and Neural Networks for an Accurate Representation of Chemically Diverse Molecular Interactions

Interplay of local hydrogen-bonding and long-ranged dipolar forces in simulations of confined water

Development of an Advanced Force Field for Water Using Variational Energy Decomposition Analysis

Modeling Liquid Water by Climbing Up Jacob's Ladder in Density Functional Theory Facilitated by Using Deep Neural Network Potentials

Interpolating many-body wave functions for accelerated molecular dynamics on the near-exact electronic surface

Gaussian Processes for Finite Size Extrapolation of Many-Body Simulations

A deep potential model with long-range electrostatic interactions

Modeling of thermophysical properties and vapor-liquid equilibrium using Gaussian process regression

Bottom-up Derived Flexible Water Model with Dipole and Quadrupole Moments for Coarse-Grained Molecular Simulations.