Fluctuations in the ensemble of reaction pathways.

The dominant reaction pathway is a rigorous framework to microscopically compute the most probable trajectories, in nonequilibrium transitions. In the low-temperature regime, such dominant pathways encode the information about the reaction mechanism and can be used to estimate nonequilibrium averages of arbitrary observables. On the other hand, at sufficiently high temperatures, the stochastic fluctuations around the dominant paths become important and have to be taken into account. In this work, we develop a technique to systematically include the effects of such stochastic fluctuations, to order k(B)T. This method is used to compute the probability for a transition to take place through a specific reaction channel and to evaluate the reaction rate.

[1]  H. Orland,et al.  Simulating stochastic dynamics using large time steps. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[2]  Erkang Wang,et al.  Quantifying kinetic paths of protein folding. , 2005, Biophysical journal.

[3]  H. Scheraga,et al.  An atomically detailed study of the folding pathways of protein A with the stochastic difference equation , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[4]  David Chandler,et al.  Statistical mechanics of isomerization dynamics in liquids and the transition state approximation , 1978 .

[5]  Communications: Ab initio dynamics of rare thermally activated reactions. , 2009, The Journal of chemical physics.

[6]  P. Faccioli.  Characterization of protein folding by dominant reaction pathways. , 2008, The journal of physical chemistry. B.

[7]  M Sega,et al.  Quantitative protein dynamics from dominant folding pathways. , 2007, Physical review letters.

[8]  Sebastian Doniach,et al.  Simulation of protein folding by reaction path annealing , 2001 .

[9]  T. Ala‐Nissila,et al.  Polymer escape from a metastable Kramers potential: path integral hyperdynamics study. , 2010, The Journal of chemical physics.

[10]  H. Orland,et al.  Dominant folding pathways of a peptide chain from ab initio quantum-mechanical simulations. , 2010, The Journal of chemical physics.

[11]  Pietro Faccioli,et al.  Dominant reaction pathways in protein folding: A direct validation against molecular dynamics simulations. , 2010, The Journal of chemical physics.

[12]  David Chandler,et al.  Transition path sampling: throwing ropes over rough mountain passes, in the dark. , 2002, Annual review of physical chemistry.

[13]  Thomas B. Woolf,et al.  Multiple pathways in conformational transitions of the alanine dipeptide: An application of dynamic importance sampling , 2006, J. Comput. Chem..

[14]  H. Orland,et al.  Dominant reaction pathways in high-dimensional systems. , 2009, The Journal of chemical physics.

[15]  B. C. Garrett,et al.  Current status of transition-state theory , 1983 .

[16]  R. Elber,et al.  Computing time scales from reaction coordinates by milestoning. , 2004, The Journal of chemical physics.

[17]  P. Hänggi,et al.  Reaction-rate theory: fifty years after Kramers , 1990 .

[18]  C. Dellago,et al.  Transition path sampling and the calculation of rate constants , 1998 .

[19]  Ron Elber,et al.  Temperature dependent reaction coordinates , 2000 .

[20]  D M Zuckerman,et al.  Efficient dynamic importance sampling of rare events in one dimension. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[21]  Ericka Stricklin-Parker,et al.  Ann , 2005 .

[22]  P. Faccioli.  Molecular dynamics at low time resolution. , 2010, The Journal of chemical physics.

[23]  H. Orland,et al.  Dominant pathways in protein folding. , 2005, Physical review letters.

[24]  P. Zanardi.  Quantum entanglement in fermionic lattices , 2002 .

[25]  Liao Y. Chen,et al.  Solution of the langevin equation for rare event rates using a path-integral formalism , 1999 .

[26]  L Y Chen,et al.  Efficient transition path sampling for systems with multiple reaction pathways. , 2005, The Journal of chemical physics.