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.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.