Rate constants for diffusive processes by partial path sampling.

We introduce a path sampling method for the computation of rate constants for complex systems with a highly diffusive character. Based on the recently developed transition interface sampling (TIS) algorithm this procedure increases the efficiency by sampling only parts of complete transition trajectories. The algorithm assumes the loss of memory for diffusive progression along the reaction coordinate. We compare the new partial path technique to the TIS method for a simple diatomic system and show that the computational effort of the new method scales linearly, instead of quadratically, with the width of the diffusive barrier. The validity of the memory loss assumption is also discussed.

[1]  Christoph Dellago,et al.  Sampling ensembles of deterministic transition pathways , 1998 .

[2]  P. Bolhuis.  Transition-path sampling of β-hairpin folding , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[3]  H. C. Andersen,et al.  Role of Repulsive Forces in Determining the Equilibrium Structure of Simple Liquids , 1971 .

[4]  Christoph Dellago,et al.  On the calculation of reaction rate constants in the transition path ensemble , 1999 .

[5]  I. Kevrekidis,et al.  Coarse molecular dynamics of a peptide fragment: Free energy, kinetics, and long-time dynamics computations , 2002, physics/0212108.

[6]  M. Sprik.  Computation of the pK of liquid water using coordination constraints , 2000 .

[7]  W. H. Weinberg,et al.  Theoretical foundations of dynamical Monte Carlo simulations , 1991 .

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

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

[10]  A. B. Bortz,et al.  A new algorithm for Monte Carlo simulation of Ising spin systems , 1975 .

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

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

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

[14]  J. Hynes,et al.  Nonequilibrium free energy surfaces for hydrogen-bonded proton-transfer complexes in solution , 1991 .

[15]  M. J. Ruiz-Montero,et al.  Efficient schemes to compute diffusive barrier crossing rates , 1997 .

[16]  D. Gillespie.  A General Method for Numerically Simulating the Stochastic Time Evolution of Coupled Chemical Reactions , 1976 .