Enhanced path sampling using subtrajectory Monte Carlo moves

Abstract

Path sampling allows the study of rare events like chemical reactions, nucleation and protein folding via a Monte Carlo (MC) exploration in path space. Instead of configuration points, this method samples short molecular dynamics (MD) trajectories with specific start- and end-conditions. As in configuration MC, its efficiency highly depends on the types of MC moves. Since the last two decades, the central MC move for path sampling has been the so-called shooting move in which a perturbed phase point of the old path is propagated backward and forward in time to generate a new path. Recently, we proposed the subtrajectory moves, stone-skipping (SS) and web-throwing (WT), that are demonstrably more efficient. However, the one-step crossing requirement makes them somewhat more difficult to implement in combination with external MD programs or when the order parameter determination is expensive. In this article, we present strategies to address the issue. The most generic solution is a new member of subtrajectory moves, wire fencing (WF), that is less thrifty than the SS, but more versatile. This makes it easier to link path sampling codes with external MD packages and provides a practical solution for cases where the calculation of the order parameter is expensive or not a simple function of geometry. We demonstrate the WF move in a double well Langevin model, a thin film breaking transition based on classical force fields, and a smaller ruthenium redox reaction at the ab initio level in which the order parameter explicitly depends on the electron density.

arXiv comment20 pages, 7 figures

References in corpus (11)

• Canonical sampling through velocity-rescaling
• Maximally localized Wannier functions: Theory and applications
• Structure identification methods for atomistic simulations of crystalline materials
• A Novel Path Sampling Method for the Calculation of Rate Constants
• Forward Flux Sampling for rare event simulations
• Elaborating Transition Interface Sampling Methods
• Rate constants for diffusive processes by partial path sampling
• Reaction rate calculation by parallel path swapping
• Forward-Flux Sampling with Jumpy Order Parameters
• Efficiency analysis of reaction rate calculation methods using analytical models I: The 2D sharp barrier
• Exchanging replicas with unequal cost, infinitely and permanently

Cited by in corpus (3)

• Exchanging replicas with unequal cost, infinitely and permanently
• Path sampling challenges in large biomolecular systems: RETIS and REPPTIS for ABL-imatinib kinetics
• Generalized Path Reweighting and History-Dependent Free Energies