A novel path sampling method for the calculation of rate constants

We derive a novel efficient scheme to measure the rate constant of transitions between stable states separated by high free energy barriers in a complex environment within the framework of transition path sampling. The method is based on directly and simultaneously measuring the fluxes through many phase space interfaces and increases the efficiency with at least a factor of 2 with respect to existing transition path sampling rate constant algorithms. The new algorithm is illustrated on the isomerization of a diatomic molecule immersed in a simple fluid.

[1]  W. Hase,et al.  Trajectory studies of SN2 nucleophilic substitution. II. Nonstatistical central barrier recrossing in the Cl−+CH3Cl system , 1992 .

[2]  James C. Keck,et al.  Statistical investigation of dissociation cross-sections for diatoms , 1962 .

[3]  J. Rivail,et al.  Molecular dynamics simulations of elementary chemical processes in liquid water using combined density functional and molecular mechanics potentials. II. Charge separation processes , 1997 .

[4]  James B. Anderson,et al.  Statistical theories of chemical reactions. Distributions in the transition region , 1973 .

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

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

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

[8]  Kihyung Song,et al.  Trajectory Studies of SN2 Nucleophilic Substitution. 8. Central Barrier Dynamics for Gas Phase Cl- + CH3Cl , 2001 .

[9]  D. Chandler,et al.  Introduction To Modern Statistical Mechanics , 1987 .

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

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

[12]  Y. Okuno.  Microscopic description of nonadiabatic, nonequilibrium, and equilibrium solvations for solvated cluster reactions: (H2O)nCl−+CH3Cl→ClCH3+Cl−(H2O)n , 1996 .

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