Onsager-Machlup action-based path sampling and its combination with replica exchange for diffusive and multiple pathways

Abstract

For sampling multiple pathways in a rugged energy landscape, we propose a novel action-based path sampling method using the Onsager-Machlup action functional. Inspired by the Fourier-path integral simulation of a quantum mechanical system, a path in Cartesian space is transformed into that in Fourier space, and an overdamped Langevin equation is derived for the Fourier components to achieve a canonical ensemble of the path at a finite temperature. To avoid "path trapping" around an initially guessed path, the path sampling method is further combined with a powerful sampling technique, the replica exchange method. The principle and algorithm of our method is numerically demonstrated for a model two-dimensional system with a bifurcated potential landscape. The results are compared with those of conventional transition path sampling and the equilibrium theory, and the error due to path discretization is also discussed.

20 pages, 5 figures, submitted to J. Chem. Phys

References in corpus (9)

• Escaping free-energy minima
• A nonequilibrium equality for free energy differences
• Parallel Tempering Algorithm for Conformational Studies of Biological Molecules
• Dominant Pathways in Protein Folding
• Quantitative Protein Dynamics from Dominant Folding Pathways
• Stochastic actions for diffusive dynamics: Reweighting, sampling, and minimization
• The Construction of Double-Ended Classical Trajectories
• Dominant Reaction Pathways in High Dimensional Systems
• Sampling diffusive transition paths

Cited by in corpus (6)

• Reinforcement learning of rare diffusive dynamics
• Metadynamics of paths
• On the hierarchical parallelization of ab initio simulations
• Multiscale enhanced path sampling based on the Onsager-Machlup action: Application to a model polymer
• Path probability distribution of stochastic motion of non dissipative systems: a classical analog of Feynman factor of path integral
• A modified nudged elastic band algorithm with adaptive spring lengths