Sampling diffusive transition paths

Abstract

We address the problem of sampling double-ended diffusive paths. The ensemble of paths is expressed using a symmetric version of the Onsager-Machlup formula, which only requires evaluation of the force field and which, upon direct time discretization, gives rise to a symmetric integrator that is accurate to second order. Efficiently sampling this ensemble requires avoiding the well-known stiffness problem associated with sampling infinitesimal Brownian increments of the path, as well as a different type of stiffness associated with sampling the coarse features of long paths. The fine-feature sampling stiffness is eliminated with the use of the fast sampling algorithm (FSA), and the coarse-feature sampling stiffness is avoided by introducing the sliding and sampling (S&S) algorithm. A key feature of the S&S algorithm is that it enables massively parallel computers to sample diffusive trajectories that are long in time. We use the algorithm to sample the transition path ensemble for the structural interconversion of the 38-atom Lennard-Jones cluster at low temperature.

13 pages 5 figures

References in corpus (6)

• String Method for the Study of Rare Events
• A Doubly Nudged Elastic Band Method for Finding Transition States
• Optimal series representations for numerical path integral simulations
• Upon the existence of short-time approximations of any polynomial order for the computation of density matrices by path integral methods
• The fast sampling algorithm for Lie-Trotter products
• On the efficient Monte Carlo implementation of path integrals

Cited by in corpus (14)

• Cayley modification for strongly stable path-integral and ring-polymer molecular dynamics
• Pathways to self-organization: crystallization via nucleation and growth
• Stochastic actions for diffusive dynamics: Reweighting, sampling, and minimization
• Flows in Complex Networks: Theory, Algorithms, and Application to Lennard-Jones Cluster Rearrangement
• Minima Hopping Guided Path Search: An Efficient Method for Finding Complex Chemical Reaction Pathways
• Simulating structural transitions by direct transition current sampling: the example of LJ38
• Onsager-Machlup action-based path sampling and its combination with replica exchange for diffusive and multiple pathways
• Global perspectives on the energy landscapes of liquids, supercooled liquids, and glassy systems: Geodesic pathways through the potential energy landscape
• Multiscale enhanced path sampling based on the Onsager-Machlup action: Application to a model polymer
• Two-scale coupling for preconditioned Hamiltonian Monte Carlo in infinite dimensions
• Path-accelerated molecular dynamics: Parallel-in-time integration using path integrals
• Estimating Reciprocal Partition Functions to Enable Design Space Sampling
• Nonequilibrium dynamics and thermodynamics provide the underlying physical mechanism of the perceptual rivalry
• Entropic effects in large-scale Monte Carlo simulations