Skip to Main content Skip to Navigation
Journal articles

A generalized parallel replica dynamics

Andrew J. Binder 1 Tony Lelièvre 2, 3 Gideon Simpson 4
ENPC - École des Ponts ParisTech, CERMICS - Centre d'Enseignement et de Recherche en Mathématiques et Calcul Scientifique, Inria Paris-Rocquencourt
Abstract : Metastability is a common obstacle to performing long molecular dynamics simulations. Many numerical methods have been proposed to overcome it. One method is parallel replica dynamics, which relies on the rapid convergence of the underlying stochastic process to a quasi-stationary distribution. Two requirements for applying parallel replica dynamics are knowledge of the time scale on which the process converges to the quasi-stationary distribution and a mechanism for generating samples from this distribution. By combining a Fleming-Viot particle system with convergence diagnostics to simultaneously identify when the process converges while also generating samples, we can address both points. This variation on the algorithm is illustrated on various numerical examples, including those with entropic barriers and the 2D Lennard-Jones cluster of seven atoms.
Complete list of metadata
Contributor : Tony Lelièvre Connect in order to contact the contributor
Submitted on : Friday, April 25, 2014 - 9:31:55 AM
Last modification on : Thursday, March 25, 2021 - 8:58:24 AM

Links full text





Andrew J. Binder, Tony Lelièvre, Gideon Simpson. A generalized parallel replica dynamics. Journal of Computational Physics, Elsevier, 2015, 284, ⟨10.1016/⟩. ⟨hal-00983282⟩



Record views