Simulating rare events in equilibrium or non-equilibrium stochastic systems

Abstract

We present three algorithms for calculating rate constants and sampling transition paths for rare events in simulations with stochastic dynamics. The methods do not require a priori knowledge of the phase space density and are suitable for equilibrium or non-equilibrium systems in stationary state. All the methods use a series of interfaces in phase space, between the initial and final states, to generate transition paths as chains of connected partial paths, in a ratchet-like manner. No assumptions are made about the distribution of paths at the interfaces. The three methods differ in the way that the transition path ensemble is generated. We apply the algorithms to kinetic Monte Carlo simulations of a genetic switch and to Langevin Dynamics simulations of intermittently driven polymer translocation through a pore. We find that the three methods are all of comparable efficiency, and that all the methods are much more efficient than brute force simulation.

18 pages, 11 figures

References in corpus (7)

• String Method for the Study of Rare Events
• A Novel Path Sampling Method for the Calculation of Rate Constants
• Sampling rare switching events in biochemical networks
• Elaborating Transition Interface Sampling Methods
• Rate constants for diffusive processes by partial path sampling
• Rate of Homogeneous Crystal Nucleation in molten NaCl
• Enhancement of the stability of genetic switches by overlapping upstream regulatory domains

Cited by in corpus (96)

• Crystal Nucleation in Liquids: Open Questions and Future Challenges in Molecular Dynamics Simulations
• Forward Flux Sampling for rare event simulations
• The Many Faces of Heterogeneous Ice Nucleation: Interplay Between Surface Morphology and Hydrophobicity
• Forward Flux Sampling-type schemes for simulating rare events: Efficiency analysis
• The Statistical Mechanics of Dynamic Pathways to Self-assembly
• Crystal nucleation of hard spheres using molecular dynamics, umbrella sampling and forward flux sampling: A comparison of simulation techniques
• Evaporation Rate of Water in Hydrophobic Confinement
• Ice Nucleation on Carbon Surface Supports the Classical Theory for Heterogeneous Nucleation
• Steady-state simulations using weighted ensemble path sampling
• Local order parameters for use in driving homogeneous ice nucleation with all-atom models of water
• Thermodynamically Driven Assemblies and Liquid-liquid Phase Separations in Biology
• Computing stationary distributions in equilibrium and non-equilibrium systems with Forward Flux Sampling
• Effective interactions and large deviations in stochastic processes
• Evidence for out-of-equilibrium crystal nucleation in suspensions of oppositely charged colloids
• Free energy landscapes for homogeneous nucleation of ice for a monatomic water model
• Surface Induced Crystallization in Supercooled Tetrahedral Liquids
• Colloidal Hard Spheres: Triumphs, Challenges and Mysteries
• Homogeneous Bubble Nucleation driven by local hot spots: a Molecular Dynamics Study
• Galerkin Approximation of Dynamical Quantities using Trajectory Data
• Studying Rare Events using Forward-Flux Sampling: Recent Breakthroughs and Future Outlook
• Rare Event Kinetics from Adaptive Bias Enhanced Sampling
• A General Algorithm for Sampling Rare Events in Non-Equilibrium and Non-Stationary Systems
• Special Topic: Markov Models of Molecular Kinetics
• Path Similarity Analysis: a Method for Quantifying Macromolecular Pathways
• Homogeneous nucleation under shear in a two-dimensional Ising model: cluster growth, coalescence and breakup
• Limit of validity of Ostwald's rule of stages in a statistical mechanical model of crystallization
• State-of-the-art models for the phase diagram of carbon and diamond nucleation
• Pathways to self-organization: crystallization via nucleation and growth
• Automatic, optimized interface placement in forward flux sampling simulations
• Efficient Stochastic Simulation of Chemical Kinetics Networks using a Weighted Ensemble of Trajectories
• Forward-Flux Sampling with Jumpy Order Parameters
• DNA looping provides stability and robustness to the bacteriophage lambda switch
• Patterning a surface so as to speed nucleation from solution
• Reaction coordinates for the flipping of genetic switches
• Enhanced Sampling in the Age of Machine Learning: Algorithms and Applications
• Computation of Large Deviation Statistics via Iterative Measurement-and-Feedback Procedure
• Efficiency analysis of reaction rate calculation methods using analytical models I: The 2D sharp barrier
• Entropy-reduced retention times in magnetic memory elements: A case of the Meyer-Neldel Compensation Rule
• Generalized convex hull construction for materials discovery
• Eliminating fast reactions in stochastic simulations of biochemical networks: a bistable genetic switch
• Non-Stationary Forward Flux Sampling
• Kinetic Accessibility of Buried DNA Sites in Nucleosomes
• Localized orientational order chaperons the nucleation of Rotator phases in hard polyhedral particles
• Specificity and Completion Time Distributions of Biochemical Processes
• Adaptively accelerating reactive molecular dynamics using boxed molecular dynamics in energy space
• Induced Charge Anisotropy: a Hidden Variable Affecting Ion Transport through Membranes
• Path sampling for lifetimes of metastable magnetic skyrmions and direct comparison with Kramers' method
• How to Quantify and Avoid Finite Size Effects in Computational Studies of Crystal Nucleation: The Case of Heterogeneous Ice Nucleation
• Practical guide to replica exchange transition interface sampling and forward flux sampling
• Reaction coordinates and rate constants for liquid droplet nucleation: quantifying the interplay between driving force and memory
• Contour forward flux sampling: Sampling rare events along multiple collective variables
• Direct evaluation of rare events in active matter from variational path sampling
• Rare Event Sampling Improves Mercury Instability Statistics
• Selective nucleation in porous media
• Competitive nucleation and the Ostwald rule in a generalized Potts model with multiple metastable phases
• PySAGES: flexible, advanced sampling methods accelerated with GPUs
• How to quantify and avoid finite size effects in computational studies of crystal nucleation: The case of homogeneous crystal nucleation
• Path Weight Sampling: Exact Monte Carlo Computation of the Mutual Information between Stochastic Trajectories
• Optimal modularity for nucleation in network-organized Ising model
• Multiscale enhanced path sampling based on the Onsager-Machlup action: Application to a model polymer
• Path sampling with Stochastic Dynamics: Some new Algorithms
• Calculating thermal stability and attempt frequency of advanced recording structures without free parameters
• Computation of nucleation of a non-equilibrium first-order phase transition using a rare-event algorithm
• The Barrier Method: A Technique for Calculating Very Long Transition Times
• Quantifying Rare Events in Stochastic Reaction-Diffusion Dynamics Using Tensor Networks
• Numerical calculation of free-energy barriers for entangled polymer nucleation
• Homogeneous Nucleation of Sheared Liquids: Advances and Insights from Simulations and Theory
• Effect of the Composition on the Free Energy of Crystal Nucleation for CuPd Nanoalloys
• The Harmonic Measure of Diffusion-Limited Aggregates including Rare Events
• Prospects of Transition Interface Sampling simulations for the theoretical study of zeolite synthesis
• Rare switching events in non-stationary systems
• Using Bayes formula to estimate rates of rare events in transition path sampling simulations
• Nucleation in scale-free networks
• Ideal Conductor Model: An analytical finite-size correction for non-equilibrium molecular dynamics simulations of ion transport through nanoporous membranes
• A simple Boxed Molecular Kinetics approach to accelerate rare events in the stochastic kinetic master equation
• Quantifying the thermal stability in perpendicularly magnetized ferromagnetic nanodisks with forward flux sampling
• Efficient dynamical correction of the transition state theory rate estimate for a flat energy barrier
• Role of Pore Dilation in Molecular Transport through the Nuclear Pore Complex: Insights from Polymer Scaling Theory
• Long-wavelength density fluctuations as nucleation precursors
• Heterogeneous path ensembles for conformational transitions in semi-atomistic models of adenylate kinase
• Nucleation pathways on complex networks
• An Efficient Self-optimized Sampling Method for Rare Events in Nonequilibrium Systems
• LeaPP: Learning Pathways to Polymorphs through machine learning analysis of atomic trajectories
• Estimating Reciprocal Partition Functions to Enable Design Space Sampling
• Heterogeneous nucleation in the random field Ising model
• Heterogeneous nucleation on complex networks with mobile impurities
• The Impact of Hydration Shell Inclusion and Chain Exclusion in the Efficacy of Reaction Coordinates for Homogeneous and Heterogeneous Ice Nucleation
• Chemical master equation parameter exploration using DMRG
• Kinetic reconstruction of free energies as a function of multiple order parameters
• Noise induces rare events in granular media
• From high-dimensional committors to reactive insights
• NATPS: Nonadiabatic Transition Path Sampling Using Time-Reversible MASH Dynamics
• Robustness of classical nucleation theory to chemical heterogeneity of crystal nucleating substrates
• The pathway elaboration method for mean first passage time estimation in large continuous-time Markov chains with applications to nucleic acid kinetics
• Folding kinetics of a polymer [corrigendum]
• Non-equilibrium pathways between cluster morphologies in active phase separation: necking, rupture and cavitation