Abstract:Rare events such as nucleation processes are of ubiquitous importance in real systems. The most popular method for nonequilibrium systems, forward flux sampling (FFS), samples rare events by using interfaces to partition the whole transition process into sequence of steps along an order parameter connecting the initial and final states. FFS usually suffers from two main difficulties: low computational efficiency due to bad interface locations and even being not applicable when trapping into unknown intermediate metastable states. In the present work, we propose an approach to overcome these difficulties, by self-adaptively locating the interfaces on the fly in an optimized manner. Contrary to the conventional FFS which set the interfaces with equal distance of the order parameter, our approach determines the interfaces with equal transition probability which is shown to satisfy the optimization condition. This is done by firstly running long local trajectories starting from the current interface λ i to get the conditional probability distribution P c (λ > λ i |λ i ), and then determining λ i+1 by equaling P c ( λ i+1|λ i ) to a give value p 0. With these optimized interfaces, FFS can be run in a much more efficient way. In addition, our approach can conveniently find the intermediate metastable states by monitoring some special long trajectories that neither end at the initial state nor reach the next interface, the number of which will increase sharply from zero if such metastable states are encountered. We apply our approach to a two-state model system and a two-dimensional lattice gas Ising model. Our approach is shown to be much more efficient than the conventional FFS method without losing accuracy, and it can also well reproduce the two-step nucleation scenario of the Ising model with easy identification of the intermediate metastable state.

Contour forward flux sampling: Sampling rare events along multiple collective variables

Nonadiabatic Forward Flux Sampling for Excited-State Rare Events

Studying rare events using forward-flux sampling: Recent breakthroughs and future outlook

Phase space and collective variable based simulation methods for studies of rare events

Enhanced Sampling with Sub-optimal Collective Variables: Reconciling Accuracy and Convergence Speed

Everything everywhere all at once, a probability-based enhanced sampling approach to rare events

Unbiasing Enhanced Sampling on a High-Dimensional Free Energy Surface with a Deep Generative Model

Consensus-based construction of high-dimensional free energy surface

Exploring high dimensional free energy landscapes: Temperature accelerated sliced sampling

Descriptors-free Collective Variables From Geometric Graph Neural Networks

Efficient Rare Event Sampling with Unsupervised Normalising Flows

Acceleration with Interpretability: Surrogate Model Based Collective Variable for Enhanced Sampling

An Efficient Self-Optimized Sampling Method for Rare Events in Nonequilibrium Systems

Enhanced Sampling of Biomolecular Slow Conformational Transitions Using Adaptive Sampling and Machine Learning

Unbiasing Enhanced Sampling on a High-dimensional Free Energy Surface with Deep Generative Model

Automatic, optimized interface placement in forward flux sampling simulations

Fast Sampling of Protein Conformational Dynamics

Enhanced Sampling of Transition States

Free energy surface reconstruction from umbrella samples using Gaussian process regression II: Multiple collective variables