Efficient Dynamic Importance Sampling of Rare Events in One Dimension

Abstract

Exploiting stochastic path integral theory, we obtain \emph{by simulation} substantial gains in efficiency for the computation of reaction rates in one-dimensional, bistable, overdamped stochastic systems. Using a well-defined measure of efficiency, we compare implementations of ``Dynamic Importance Sampling'' (DIMS) methods to unbiased simulation. The best DIMS algorithms are shown to increase efficiency by factors of approximately 20 for a $5k_{B}T$ barrier height and 300 for $9k_{B}T$, compared to unbiased simulation. The gains result from close emulation of natural (unbiased), instanton-like crossing events with artificially decreased waiting times between events that are corrected for in rate calculations. The artificial crossing events are generated using the closed-form solution to the most probable crossing event described by the Onsager-Machlup action. While the best biasing methods require the second derivative of the potential (resulting from the ``Jacobian'' term in the action, which is discussed at length), algorithms employing solely the first derivative do nearly as well. We discuss the importance of one-dimensional models to larger systems, and suggest extensions to higher-dimensional systems.

version to be published in Phys. Rev. E

References in corpus (1)

• Computing probabilities of very rare events for Langevin processes: a new method based on importance sampling

Cited by in corpus (18)

• Transition-Event Durations in One Dimensional Activated Processes
• The mean shape of transition and first-passage paths
• Stochastic actions for diffusive dynamics: Reweighting, sampling, and minimization
• Transition Events in Butane Simulations: Similarities Across Models
• Weighted Ensemble Milestoning (WEM): A Combined Approach for Rare Event Simulations
• Dynamical reweighting methods for Markov models
• Steered Transition Path Sampling
• Summing over trajectories of stochastic dynamics with multiplicative noise
• A review of Girsanov Reweighting and of Square Root Approximation for building molecular Markov State Models
• Fluctuations in the Ensemble of Reaction Pathways
• Path probability ratios for Langevin dynamics -- exact and approximate
• Rapid Determination of Multiple Reaction Pathways in Molecular Systems: The Soft-Ratcheting Algorithm
• Estimating the Most Probable Transition Time for Stochastic Dynamical Systems
• Path-accelerated molecular dynamics: Parallel-in-time integration using path integrals
• Statistical Mechanics of Histories: A Cluster Monte Carlo Algorithm
• Iterative trajectory reweighting for estimation of equilibrium and non-equilibrium observables
• Characterization of the Most Probable Transition Paths of Stochastic Dynamical Systems with Stable Lévy Noise
• Unbiased simulation of structural transitions in calmodulin