Premium
TheRate: Program for ab initio direct dynamics calculations of thermal and vibrational‐state‐selected rate constants
Author(s) -
Duncan Wendell T.,
Bell Robert L.,
Truong Thanh N.
Publication year - 1998
Publication title -
journal of computational chemistry
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.907
H-Index - 188
eISSN - 1096-987X
pISSN - 0192-8651
DOI - 10.1002/(sici)1096-987x(19980715)19:9<1039::aid-jcc5>3.0.co;2-r
Subject(s) - hessian matrix , semiclassical physics , ab initio , density functional theory , interpolation (computer graphics) , electronic structure , potential energy surface , chemistry , computational chemistry , quantum mechanics , physics , mathematics , quantum , classical mechanics , motion (physics)
We introduce TheRate (THEoretical RATEs), a complete application program with a graphical user interface (GUI) for calculating rate constants from first principles. It is based on canonical variational transition‐state theory (CVT) augmented by multidimensional semiclassical zero and small curvature tunneling approximations. Conventional transition‐state theory (TST) with one‐dimensional Wigner or Eckart tunneling corrections is also available. Potential energy information needed for the rate calculations are obtained from ab initio molecular orbital and/or density functional electronic structure theory. Vibrational‐state‐selected rate constants may be calculated using a diabetic model. TheRate also introduces several technical advancements, namely the focusing technique and energy interpolation procedure. The focusing technique minimizes the number of Hessian calculations required by distributing more Hessian grid points in regions that are critical to the CVT and tunneling calculations and fewer Hessian grid points elsewhere. The energy interpolation procedure allows the use of a computationally less demanding electronic structure theory such as DFT to calculate the Hessians and geometries, while the energetics can be improved by performing a small number of single‐point energy calculations along the MEP at a more accurate level of theory. The CH 4 +H↔CH 3 +H 2 reaction is used as a model to demonstrate usage of the program, and the convergence of the rate constants with respect to the number of electronic structure calculations. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1039–1052, 1998