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.

Simulating rare events using a Weighted Ensemble-based string method

Rare-Event Sampling using a Reinforcement Learning-Based Weighted Ensemble Method

Weighted ensemble: Recent mathematical developments

Weighted ensemble milestoning (WEM): A combined approach for rare event simulations

Structuring and sampling complex conformation space: Weighted ensemble dynamics simulations

Unsupervised learning of progress coordinates during weighted ensemble simulations: Application to millisecond protein folding

Augmenting Human Expertise in Weighted Ensemble Simulations through Deep Learning based Information Bottleneck

Finite Temperature String Method for the Study of Rare Events

Equilibrium sampling by re-weighting non-equilibrium simulation trajectories

Lattice Models for Magnetic Fluids: Correlations Between Order Parameters

Enhanced Sampling of Configuration and Path Space in a Generalized Ensemble by Shooting Point Exchange

BAD-NEUS: Rapidly converging trajectory stratification

Integrating Path Sampling with Enhanced Sampling for Rare-event Kinetics

RiteWeight: Randomized Iterative Trajectory Reweighting for Steady-State Distributions Without Discretization Error

Markov State Models with Weighted Ensemble Simulation: How to Eliminate the Trajectory Merging Bias

Path Integral Coarse-Graining Replica Exchange Method for Enhanced Sampling

Multireference Generalization of the Weighted Thermodynamic Perturbation Method.

Enhanced Sampling in the Well-Tempered Ensemble

Phase space and collective variable based simulation methods for studies of rare events

Discriminant Analysis Optimizes Progress Coordinate in Weighted Ensemble Simulations of Rare Event Kinetics