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.

Simplified algorithm for the Worldvolume HMC and the Generalized-thimble HMC

Simplified algorithm for the Worldvolume HMC and the Generalized-thimble HMC

Generalized HMC using Nambu mechanics for lattice QCD

Multireference Generalization of the Weighted Thermodynamic Perturbation Method.

Lagrangian Dynamical Monte Carlo

HMC and underdamped Langevin united in the unadjusted convex smooth case

MG-MLMC++ as a Variance Reduction Method for Estimating the Trace of a Matrix Inverse

Applying the tempered Lefschetz thimble method to the Hubbard model away from half-filling

Effective and Efficient Algorithm for the Wigner Rotation Matrix at High Angular Momenta

Weight-adjusted discontinuous Galerkin methods: curvilinear meshes

Stand-alone Multigrid for Helmholtz Revisited: Towards Convergence Using Standard Components

Overcoming Ergodicity Problems of the Hybrid Monte Carlo Method using Radial Updates

Towards the efficient calculation of quantity of interest from steady Euler equations I: a dual-consistent DWR-based h-adaptive Newton-GMG solver

Tempered Lefschetz thimble method and its application to the Hubbard model away from half filling

The No-U-Turn Sampler: Adaptively Setting Path Lengths in Hamiltonian Monte Carlo

A multiscale hybrid Maxwellian-Monte-Carlo Coulomb collision algorithm for particle simulations

Geometric integrators and the Hamiltonian Monte Carlo method

Markov Chain Monte Carlo from Lagrangian Dynamics

Efficient energy minimization in finite-difference micromagnetics: Speeding up hysteresis computations

Riemannian Manifold HMC with fermions

Exponential speed up in Monte Carlo sampling through Radial Updates