TY - JOUR
T1 - Quantum mechanical reaction rate constants by vibrational configuration interaction
T2 - The OH + H2→H2O + H reaction as a function of temperature
AU - Chakraborty, Arindam
AU - Truhlar, Donald G.
PY - 2005/5/10
Y1 - 2005/5/10
N2 - The thermal rate constant of the 3D OH + H2→H2O + H reaction was computed by using the flux autocorrelation function, with a time-independent square-integrable basis set. Two modes that actively participate in bond making and bond breaking were treated by using 2D distributed Gaussian functions, and the remaining (nonreactive) modes were treated by using harmonic oscillator functions. The finite-basis eigenvalues and eigenvectors of the Hamiltonian were obtained by solving the resulting generalized eigenvalue equation, and the flux autocorrelation function for a dividing surface optimized in reduced-dimensionality calculations was represented in the basis formed by the eigenvectors of the Hamiltonian. The rate constant was obtained by integrating the flux autocorrelation function. The choice of the final time to which the integration is carried was determined by a plateau criterion. The potential energy surface was from Wu, Schatz, Lendvay, Fang, and Harding (WSLFH). We also studied the collinear H + Hz reaction by using the Liu-Siegbahn-Truhlar-Horowitz (LSTH) potential energy surface. The calculated thermal rate constant results were compared with reported values on the same surfaces. The success of these calculations demonstrates that time-independent vibrational configuration interaction can be a very convenient way to calculate converged quantum mechanical rate constants, and it opens the possibility of calculating converged rate constants for much larger reactions than have been treated until now.
AB - The thermal rate constant of the 3D OH + H2→H2O + H reaction was computed by using the flux autocorrelation function, with a time-independent square-integrable basis set. Two modes that actively participate in bond making and bond breaking were treated by using 2D distributed Gaussian functions, and the remaining (nonreactive) modes were treated by using harmonic oscillator functions. The finite-basis eigenvalues and eigenvectors of the Hamiltonian were obtained by solving the resulting generalized eigenvalue equation, and the flux autocorrelation function for a dividing surface optimized in reduced-dimensionality calculations was represented in the basis formed by the eigenvectors of the Hamiltonian. The rate constant was obtained by integrating the flux autocorrelation function. The choice of the final time to which the integration is carried was determined by a plateau criterion. The potential energy surface was from Wu, Schatz, Lendvay, Fang, and Harding (WSLFH). We also studied the collinear H + Hz reaction by using the Liu-Siegbahn-Truhlar-Horowitz (LSTH) potential energy surface. The calculated thermal rate constant results were compared with reported values on the same surfaces. The success of these calculations demonstrates that time-independent vibrational configuration interaction can be a very convenient way to calculate converged quantum mechanical rate constants, and it opens the possibility of calculating converged rate constants for much larger reactions than have been treated until now.
KW - Chemical kinetics
KW - Dynamics
KW - Flux autocorrelation
KW - Localized basis functions
KW - Many-body problem
UR - http://www.scopus.com/inward/record.url?scp=18744363042&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=18744363042&partnerID=8YFLogxK
U2 - 10.1073/pnas.0408048102
DO - 10.1073/pnas.0408048102
M3 - Article
C2 - 15774583
AN - SCOPUS:18744363042
SN - 0027-8424
VL - 102
SP - 6744
EP - 6749
JO - Proceedings of the National Academy of Sciences of the United States of America
JF - Proceedings of the National Academy of Sciences of the United States of America
IS - 19
ER -