ON THE DETERMINATION OF RATE CONSTANTS FROM EQUILIBRIUM MOLECULAR-DYNAMICS SIMULATIONS

被引:21
作者
BROWN, D
CLARKE, JHR
机构
[1] Chemistry Department, U.M.I.S.T., Manchester M60 1QD, Sackville Street
关键词
D O I
10.1063/1.458744
中图分类号
O64 [物理化学(理论化学)、化学物理学];
学科分类号
070304 ; 081704 ;
摘要
In this article we consider the calculation of rate constants from equilibrium molecular dynamics simulations using the conformational isomerization in liquid butane as an example. In particular we discuss the method by which transition state theory estmates of the rate constant may be calculated from simulation data. Care has to be taken when calculating the average velocity along the reaction coordinate at the barrier. If this is computed just from actual barrier crossings we show that there is a systematic overestimation by a factor of exactly π/2. Further analysis has allowed the connection between the barrier crossing rate and the transition state theory value for the rate constant to be established. We compare rate constants determined by the "reactive flux" method with values obtained using alternative approaches applied to liquid n-butane [D. Brown and J. H. R. Clarke, J. Chem. Phys. 92, 3062 (1990)]. The results of this comparison have implications regarding the identification of the plateau value when using the reactive flux method. A straightforward procedure is described which can establish confidence limits for rate constants calculated in this way. © 1990 American Institute of Physics.
引用
收藏
页码:4117 / 4122
页数:6
相关论文
共 10 条