Pitfalls of choosing an order parameter for rare event calculations.

The mechanism of rare events in complex systems can be found by sampling dynamical paths that connect stable states. To calculate a rate using transition paths, an order parameter is required to describe the progress of the reaction and to distinguish the initial and final states. In this work, we compare two implementations of transition path sampling for Langevin paths, one for which paths are sampled in configuration space and the other in the space of the random variables that describe the thermostat. These two approaches are found to give different rates for the rearrangement of a seven-particle cluster despite the fact that both are formally exact. The difference is understood in terms of how efficiently the methods sample states along the order parameter. The more efficient approach takes the system to unexpected states that are allowed by a poor choice of order parameter. While transition path sampling is formally correct, we show how mistakes can be made when the system escapes to unknown states along an order parameter represented in terms of collective variables.

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

[2]  Alessandro Laio,et al.  A recipe for the computation of the free energy barrier and the lowest free energy path of concerted reactions. , 2005, The journal of physical chemistry. B.

[3]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[4]  G. Crooks,et al.  Efficient transition path sampling for nonequilibrium stochastic dynamics. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

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

[6]  Christoph Dellago,et al.  Efficient transition path sampling: Application to Lennard-Jones cluster rearrangements , 1998 .

[7]  G. Henkelman,et al.  Improved tangent estimate in the nudged elastic band method for finding minimum energy paths and saddle points , 2000 .

[8]  D. Wales Discrete path sampling , 2002 .

[9]  Eugene P. Wigner,et al.  The transition state method , 1938 .

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

[11]  S. Schwartz,et al.  How enzyme dynamics helps catalyze a reaction in atomic detail: a transition path sampling study. , 2005, Journal of the American Chemical Society.

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

[13]  Roberto Car,et al.  Free energy profile along a discretized reaction path via the hyperplane constraint force and torque. , 2005, The Journal of chemical physics.

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

[15]  Eric Vanden-Eijnden,et al.  Transition state theory: variational formulation, dynamical corrections, and error estimates. , 2005, The Journal of chemical physics.

[16]  B M Dickson,et al.  Rapid sampling of reactive Langevin trajectories via noise-space Monte Carlo. , 2007, The Journal of chemical physics.

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

[18]  Dmitrii E Makarov,et al.  The rate constant of polymer reversal inside a pore. , 2008, The Journal of chemical physics.

[19]  L Y Chen,et al.  Transition rate prefactors for systems of many degrees of freedom. , 2006, The Journal of chemical physics.

[20]  G. Henkelman,et al.  A dimer method for finding saddle points on high dimensional potential surfaces using only first derivatives , 1999 .

[21]  H. Eyring The Activated Complex in Chemical Reactions , 1935 .

[22]  Herbert Levine,et al.  Effective stochastic dynamics on a protein folding energy landscape. , 2006, The Journal of chemical physics.

[23]  Richard J. Dimelow,et al.  Exploring reaction pathways with transition path and umbrella sampling: application to methyl maltoside. , 2006, The Journal of chemical physics.