Using Bayes formula to estimate rates of rare events in transition path sampling simulations
arXiv:1509.08371 · doi:10.1063/1.4932389
Abstract
Transition path sampling is a method for estimating the rates of rare events in molecular systems based on the gradual transformation of a path distribution containing a small fraction of reactive trajectories into a biased distribution in which these rare trajectories have become frequent. Then, a multistate reweighting scheme is implemented to postprocess data collected from the staged simulations. Herein, we show how Bayes formula allows to directly construct a biased sample containing an enhanced fraction of reactive trajectories and to concomitantly estimate the transition rate from this sample. The approach can remediate the convergence issues encountered in free energy perturbation or umbrella sampling simulations when the transformed distribution insufficiently overlaps with the reference distribution.
arXiv comment11 pages, 8 figures
References in corpus (15)
- Statistically optimal analysis of samples from multiple equilibrium states
- Development of an interatomic potential for phosphorus impurities in iron
- A Novel Path Sampling Method for the Calculation of Rate Constants
- Simulating rare events in equilibrium or non-equilibrium stochastic systems
- Forward Flux Sampling-type schemes for simulating rare events: Efficiency analysis
- Extended Ensemble Monte Carlo
- Some improvements of the ART method for finding transition pathways on potential energy surfaces
- Finite-temperature critical point of a glass transition
- Path Factorization Approach to Stochastic Simulations
- Free energy reconstruction from steered dynamics without post-processing
- Free energy calculations from adaptive molecular dynamics simulations with adiabatic reweighting
- Identifying rare chaotic and regular trajectories in dynamical systems with Lyapunov weighted path sampling
- Waste-recycling Monte Carlo with optimal estimates: application to free energy calculations in alloys
- Path sampling with Stochastic Dynamics: Some new Algorithms
- Estimating time-correlation functions by sampling and unbiasing dynamically activated events