Computing transition path theory quantities with trajectory stratification

Abstract

Transition path theory computes statistics from ensembles of reactive trajectories. A common strategy for sampling reactive trajectories is to control the branching and pruning of trajectories so as to enhance the sampling of low probability segments. However, it can be challenging to apply transition path theory to data from such methods because determining whether configurations and trajectory segments are part of reactive trajectories requires looking backward and forward in time. Here, we show how this issue can be overcome efficiently by introducing simple data structures. We illustrate the approach in the context of nonequilibrium umbrella sampling (NEUS), but the strategy is general and can be used to obtain transition path theory statistics from other methods that sample segments of unbiased trajectories.

arXiv comment12 pages, 11 figures

References in corpus (12)

• PLUMED 2: New feathers for an old bird
• PLUMED: a portable plugin for free-energy calculations with molecular dynamics
• A Novel Path Sampling Method for the Calculation of Rate Constants
• Sampling rare switching events in biochemical networks
• Forward Flux Sampling for rare event simulations
• Variational Koopman models: slow collective variables and molecular kinetics from short off-equilibrium simulations
• Steady-state simulations using weighted ensemble path sampling
• Galerkin Approximation of Dynamical Quantities using Trajectory Data
• Steered Transition Path Sampling
• Eigenvector method for umbrella sampling enables error analysis
• Augmented Transition Path Theory for Sequences of Events
• Unbiased estimation of equilibrium, rates, and committors from Markov state model analysis

Cited by in corpus (5)

• Computing committors via Mahalanobis diffusion maps with enhanced sampling data
• Inexact iterative numerical linear algebra for neural network-based spectral estimation and rare-event prediction
• Variational deep learning of equilibrium transition path ensembles
• Augmented Transition Path Theory for Sequences of Events
• From high-dimensional committors to reactive insights