4.7 Article

Enhanced path sampling using subtrajectory Monte Carlo moves

Journal

JOURNAL OF CHEMICAL PHYSICS
Volume 158, Issue 2, Pages -

Publisher

AIP Publishing
DOI: 10.1063/5.0127249

Keywords

-

Ask authors/readers for more resources

Path sampling is a method that allows the study of rare events in chemistry and biology by sampling short molecular dynamics trajectories. The efficiency of this method depends on the type of Monte Carlo moves used. Recent developments have proposed new moves, such as stone-skipping and web-throwing, which are more efficient but more challenging to implement. This article presents strategies to address this challenge, including a new move called wire fencing that is less efficient but more versatile, making it easier to integrate path sampling with external software and handle expensive order parameter calculations.
Path sampling allows the study of rare events, such as 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, 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.

Authors

I am an author on this paper
Click your name to claim this paper and add it to your profile.

Reviews

Primary Rating

4.7
Not enough ratings

Secondary Ratings

Novelty
-
Significance
-
Scientific rigor
-
Rate this paper

Recommended

No Data Available
No Data Available