Nonlinear reaction coordinate analysis in the reweighted path ensemble.

We present a flexible nonlinear reaction coordinate analysis method for the transition path ensemble based on the likelihood maximization approach developed by Peters and Trout [J. Chem. Phys. 125, 054108 (2006)]. By parametrizing the reaction coordinate by a string of images in a collective variable space, we can optimize the likelihood that the string correctly models the committor data obtained from a path sampling simulation. The collective variable space with the maximum likelihood is considered to contain the best description of the reaction. The use of the reweighted path ensemble [J. Rogal et al., J. Chem. Phys. 133, 174109 (2010)] allows a complete reaction coordinate description from the initial to the final state. We illustrate the method on a z-shaped two-dimensional potential. While developed for use with path sampling, this analysis method can also be applied to regular molecular dynamics trajectories.

[1]  P. Bolhuis,et al.  Rate constant and reaction coordinate of Trp-cage folding in explicit water. , 2008, Biophysical journal.

[2]  P. Bolhuis.  Rare events via multiple reaction channels sampled by path replica exchange. , 2008, The Journal of chemical physics.

[3]  W. E,et al.  Towards a Theory of Transition Paths , 2006 .

[4]  G. Ciccotti,et al.  String method in collective variables: minimum free energy paths and isocommittor surfaces. , 2006, The Journal of chemical physics.

[5]  Eric Vanden-Eijnden,et al.  Revisiting the finite temperature string method for the calculation of reaction tubes and free energies. , 2009, The Journal of chemical physics.

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

[7]  Cindy Starbuck,et al.  Surface-mediated nucleation in the solid-state polymorph transformation of terephthalic acid. , 2007, Journal of the American Chemical Society.

[8]  G. Hummer,et al.  Reaction coordinates and rates from transition paths. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[9]  A. Laio,et al.  Escaping free-energy minima , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[10]  G. Ciccotti,et al.  Constrained reaction coordinate dynamics for the simulation of rare events , 1989 .

[11]  Ernesto E. Borrero,et al.  Reaction coordinates and transition pathways of rare events via forward flux sampling. , 2007, The Journal of chemical physics.

[12]  G. Hummer.  From transition paths to transition states and rate coefficients. , 2004, The Journal of chemical physics.

[13]  B. Trout,et al.  Obtaining reaction coordinates by likelihood maximization. , 2006, The Journal of chemical physics.

[14]  Philipp Metzner,et al.  Illustration of transition path theory on a collection of simple examples. , 2006, The Journal of chemical physics.

[15]  Aaron R Dinner,et al.  Automatic method for identifying reaction coordinates in complex systems. , 2005, The journal of physical chemistry. B.

[16]  Peter G. Bolhuis,et al.  Transition path sampling on diffusive barriers , 2003 .

[17]  A. Berezhkovskii,et al.  One-dimensional reaction coordinates for diffusive activated rate processes in many dimensions. , 2005, The Journal of chemical physics.

[18]  T. S. V. Erp,et al.  Reaction rate calculation by parallel path swapping. , 2007, Physical review letters.

[19]  C. Dellago,et al.  Reaction coordinates of biomolecular isomerization. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[20]  Eric Vanden-Eijnden,et al.  Transition-path theory and path-finding algorithms for the study of rare events. , 2010, Annual review of physical chemistry.

[21]  Peter G. Bolhuis,et al.  Predicting the reaction coordinates of millisecond light-induced conformational changes in photoactive yellow protein , 2010, Proceedings of the National Academy of Sciences.

[22]  P. Bolhuis,et al.  Sampling the multiple folding mechanisms of Trp-cage in explicit solvent , 2006, Proceedings of the National Academy of Sciences.

[23]  Peter G. Bolhuis,et al.  A novel path sampling method for the calculation of rate constants , 2003 .

[24]  J. Langer.  Statistical theory of the decay of metastable states , 1969 .

[25]  Bernhardt L Trout,et al.  Extensions to the likelihood maximization approach for finding reaction coordinates. , 2007, The Journal of chemical physics.

[26]  Alan M. Ferrenberg,et al.  Optimized Monte Carlo data analysis. , 1989, Physical Review Letters.

[27]  G. Torrie,et al.  Monte Carlo free energy estimates using non-Boltzmann sampling: Application to the sub-critical Lennard-Jones fluid , 1974 .

[28]  G. Schwarz.  Estimating the Dimension of a Model , 1978 .

[29]  E. Vanden-Eijnden,et al.  String method for the study of rare events , 2002, cond-mat/0205527.