A Generalized Parallel Replica Dynamics
arXiv:1404.6191 · doi:10.1016/j.jcp.2015.01.002
Abstract
Metastability is a common obstacle to performing long molecular dynamics simulations. Many numerical methods have been proposed to overcome it. One method is parallel replica dynamics, which relies on the rapid convergence of the underlying stochastic process to a quasi-stationary distribution. Two requirements for applying parallel replica dynamics are knowledge of the time scale on which the process converges to the quasi-stationary distribution and a mechanism for generating samples from this distribution. By combining a Fleming-Viot particle system with convergence diagnostics to simultaneously identify when the process converges while also generating samples, we can address both points. This variation on the algorithm is illustrated with various numerical examples, including those with entropic barriers and the 2D Lennard-Jones cluster of seven atoms.
24 pages, 15 figures, 6 tables
References in corpus (1)
Cited by in corpus (16)
- The Seven Deadly Sins: when computing crystal nucleation rates, the devil is in the details
- Jump Markov models and transition state theory: the Quasi-Stationary Distribution approach
- Local and Global Perspectives on Diffusion Maps in the Analysis of Molecular Systems
- Accelerated dynamics: Mathematical foundations and algorithmic improvements
- A framework of nonequilibrium statistical mechanics. II. Coarse-graining
- The parallel replica method for simulating long trajectories of Markov chains
- gen.parRep: a first implementation of the Generalized Parallel Replica dynamics for the long time simulation of metastable biochemical systems
- Spin-Diffusions and Diffusive Molecular Dynamics
- Analysis of estimators for adaptive Kinetic Monte Carlo
- Parallel replica dynamics method for bistable stochastic reaction networks: simulation and sensitivity analysis
- Computing Long Timescale Biomolecular Dynamics using Quasi-Stationary Distribution Kinetic Monte Carlo (QSD-KMC)
- Coarse-graining via the fluctuation-dissipation theorem and large-deviation theory
- Sharp asymptotics of the first exit point density
- Stationary averaging for multi-scale continuous time Markov chains using parallel replica dynamics
- Generalizing Parallel Replica Dynamics: Trajectory Fragments, Asynchronous Computing, and PDMPs
- Simulating the stochastic dynamics and cascade failure of power networks