TY - JOUR

T1 - Numerical analysis of parallel replica dynamics

AU - Simpson, Gideon

AU - Luskin, Mitchell

PY - 2013/9

Y1 - 2013/9

N2 - Parallel replica dynamics is a method for accelerating the computation of processes characterized by a sequence of infrequent events. In this work, the processes are governed by the overdamped Langevin equation. Such processes spend much of their time about the minima of the underlying potential, occasionally transitioning into different basins of attraction. The essential idea of parallel replica dynamics is that the exit distribution from a given well for a single process can be approximated by the distribution of the first exit of N independent identical processes, each run for only 1 / N-th the amount of time. While promising, this leads to a series of numerical analysis questions about the accuracy of the exit distributions. Building upon the recent work in [C. Le Bris, T. Lelièvre, M. Luskin and D. Perez, Monte Carlo Methods Appl. 18 (2012) 119-146], we prove a unified error estimate on the exit distributions of the algorithm against an unaccelerated process. Furthermore, we study a dephasing mechanism, and prove that it will successfully complete.

AB - Parallel replica dynamics is a method for accelerating the computation of processes characterized by a sequence of infrequent events. In this work, the processes are governed by the overdamped Langevin equation. Such processes spend much of their time about the minima of the underlying potential, occasionally transitioning into different basins of attraction. The essential idea of parallel replica dynamics is that the exit distribution from a given well for a single process can be approximated by the distribution of the first exit of N independent identical processes, each run for only 1 / N-th the amount of time. While promising, this leads to a series of numerical analysis questions about the accuracy of the exit distributions. Building upon the recent work in [C. Le Bris, T. Lelièvre, M. Luskin and D. Perez, Monte Carlo Methods Appl. 18 (2012) 119-146], we prove a unified error estimate on the exit distributions of the algorithm against an unaccelerated process. Furthermore, we study a dephasing mechanism, and prove that it will successfully complete.

KW - Accelerated dynamics

KW - Parallel replica

KW - Rare events

UR - http://www.scopus.com/inward/record.url?scp=84880194242&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84880194242&partnerID=8YFLogxK

U2 - 10.1051/m2an/2013068

DO - 10.1051/m2an/2013068

M3 - Article

AN - SCOPUS:84880194242

VL - 47

SP - 1287

EP - 1314

JO - Mathematical Modelling and Numerical Analysis

JF - Mathematical Modelling and Numerical Analysis

SN - 0764-583X

IS - 5

ER -