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

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

SN - 0027-8424

IS - 19

ER -