Statistical behavior of adaptive multilevel splitting algorithms in simple models
For researchers using adaptive multilevel splitting in rare event simulation, this work highlights the critical importance of reaction coordinate choice and provides theoretical insights into algorithm behavior.
The paper analyzes adaptive multilevel splitting algorithms for estimating tail distributions and reactive trajectories in simple stochastic models. It shows that with non-optimal reaction coordinates, the algorithm's variance peaks at phase transitions, convergence slows, and reactive trajectories deviate from the most probable ones, while optimal coordinates (committor) resolve these issues.
Adaptive multilevel splitting algorithms have been introduced rather recently for estimating tail distributions in a fast and efficient way. In particular, they can be used for computing the so-called reactive trajectories corresponding to direct transitions from one metastable state to another. The algorithm is based on successive selection-mutation steps performed on the system in a controlled way. It has two intrinsic parameters, the number of particles/trajectories and the reaction coordinate used for discriminating good or bad trajectories. We investigate first the convergence in law of the algorithm as a function of the timestep for several simple stochastic models. Second, we consider the average duration of reactive trajectories for which no theoretical predictions exist. The most important aspect of this work concerns some systems with two degrees of freedom. They are studied in details as a function of the reaction coordinate in the asymptotic regime where the number of trajectories goes to infinity. We show that during phase transitions, the statistics of the algorithm deviate significatively from known theoretical results when using non-optimal reaction coordinates. In this case, the variance of the algorithm is peaking at the transition and the convergence of the algorithm can be much slower than the usual expected central limit behavior. The duration of trajectories is affected as well. Moreover, reactive trajectories do not correspond to the most probable ones. Such behavior disappears when using the optimal reaction coordinate called committor as predicted by the theory. We finally investigate a three-state Markov chain which reproduces this phenomenon and show logarithmic convergence of the trajectory durations.