377 articles – 175 references  [version française]
HAL: hal-00506047, version 1

Detailed view  Export this paper
Journal of Chemical Physics 134, 5 (2011) 054108 (16 pages)
A multiple replica approach to simulate reactive trajectories
Frédéric Cérou 1, Arnaud Guyader 1, 2, Tony Lelièvre 3, 4, David Pommier 3, 4
(2011)

A method to generate reactive trajectories, namely equilibrium trajectories leaving a metastable state and ending in another one is proposed. The algorithm is based on simulating in parallel many copies of the system, and selecting the replicas which have reached the highest values along a chosen one-dimensional reaction coordinate. This reaction coordinate does not need to precisely describe all the metastabilities of the system for the method to give reliable results. An extension of the algorithm to compute transition times from one metastable state to another one is also presented. We demonstrate the interest of the method on two simple cases: a one-dimensional two-well potential and a two-dimensional potential exhibiting two channels to pass from one metastable state to another one.
1:  ASPI (INRIA - IRISA)
CNRS : UMR6074 – INRIA – Université de Rennes 1
2:  Département Mathématiques appliquées et sciences sociales - Rennes 2 (MASS)
Université de Rennes II - Haute Bretagne
3:  Centre d'Enseignement et de Recherche en Mathématiques et Calcul Scientifique (CERMICS)
Ecole des Ponts ParisTech
4:  MICMAC (INRIA Paris - Rocquencourt)
Ecole des Ponts ParisTech – INRIA
Physics/Condensed Matter/Statistical Mechanics
Reactive trajectories – multilevel splitting algorithm
Fulltext link: 
http://fr.arXiv.org/abs/1007.4481