We introduce a new approach to evaluate transition rates for rare events in complex many-particle systems. Building on a path-integral representation of transition probabilities for Markov processes, the rate is first expressed in terms of a free energy in the transition-path ensemble. We then define an auxiliary process where a suitably defined reaction variable is dynamically decoupled from all the others, whose dynamics is left unchanged. For this system the transition rates coincide with those of a unidimensional process whose only coordinate is the reaction variable. The free-energy difference between the auxiliary and the physical transition-path ensembles is finally evaluated using standard techniques. The efficiency of this method is deemed to be optimal because the physical and auxiliary dynamics differ by one degree of freedom only at any system size. Our method is demonstrated numerically on a simple model of Lennard-Jones particles ruled by the overdamped Langevin equation.

Absolute Transition Rates for Rare Events from Dynamical Decoupling of Reaction Variables / Gobbo, Gianpaolo; Laio, A.; Maleki, Armin; Baroni, S.. - In: PHYSICAL REVIEW LETTERS. - ISSN 0031-9007. - 109:15(2012), pp. 150601.1-150601.5. [10.1103/PhysRevLett.109.150601]

Absolute Transition Rates for Rare Events from Dynamical Decoupling of Reaction Variables

Gobbo, Gianpaolo;Laio, A.;Maleki, Armin;Baroni, S.
2012-01-01

Abstract

We introduce a new approach to evaluate transition rates for rare events in complex many-particle systems. Building on a path-integral representation of transition probabilities for Markov processes, the rate is first expressed in terms of a free energy in the transition-path ensemble. We then define an auxiliary process where a suitably defined reaction variable is dynamically decoupled from all the others, whose dynamics is left unchanged. For this system the transition rates coincide with those of a unidimensional process whose only coordinate is the reaction variable. The free-energy difference between the auxiliary and the physical transition-path ensembles is finally evaluated using standard techniques. The efficiency of this method is deemed to be optimal because the physical and auxiliary dynamics differ by one degree of freedom only at any system size. Our method is demonstrated numerically on a simple model of Lennard-Jones particles ruled by the overdamped Langevin equation.
2012
109
15
1
5
150601
https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.109.150601
Gobbo, Gianpaolo; Laio, A.; Maleki, Armin; Baroni, S.
File in questo prodotto:
File Dimensione Formato  
PhysRevLett.109.150601.pdf

non disponibili

Tipologia: Versione Editoriale (PDF)
Licenza: Non specificato
Dimensione 532.7 kB
Formato Adobe PDF
532.7 kB Adobe PDF   Visualizza/Apri   Richiedi una copia

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/20.500.11767/12414
Citazioni
  • ???jsp.display-item.citation.pmc??? 0
  • Scopus 9
  • ???jsp.display-item.citation.isi??? 9
social impact