Quantum mechanical reaction rate constants by vibrational configuration interaction: The OH + H₂→H₂O + H reaction as a function of temperature

The thermal rate constant of the 3D OH + H(2)-->H(2)O + 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 + H(2) 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.