Comp Chem Research Developments

November 8, 2004

Quantum mechanical reaction rate constants by vibrational configuration interaction

Time-dependent quantum mechanics has a number of advantages for calculations of the reaction dynamics of polyatomic systems. It is especially efficient for generating results at a series of total energies, as required for thermally averaged rate constants, because a single wave packet carries information about a wide range of energies. Recent advances in the theory of vibrational configuration interaction method [A. Chakraborty, D. G. Truhlar, J. M. Bowman, and S. Carter, Journal of Chemical Physics 121, 2071 (2004), and references therein] have prompted a re-examination of the power of time-independent quantum mechanics, and now, for the first time, a converged quantum mechanical rate constant has been calculated for a polyyat00omic reaction by time-independent quantum mechanics. This work was carried out by graduate student Arindam Chakraborty and Professor Don Truhlar using the OH + H2 ® H2O + H reaction as a prototype reaction to demonstrate the method. This application is a four-body system, but the method used here is not special to four-body reactions, and the results demonstrate that the time-independent vibrational configuration method can be a very convenient way to calculate converged quantum mechanical rate constants. Implementation of hierarchical expansion of the potential energy function in the present method will make the calculations affordable for larger systems. This work lays the groundwork for calculating converged rate constants for much larger reactions than have been treated up to now.

The rate constant calculations were performed using flux autocorrelation functions, with a time-independent square-integrable basis set. Two modes that actively participate in bond forming and bond breaking were treated using two-dimensional distributed Gaussian functions, and the remaining (nonreactive) modes were treated 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 was represented in the basis formed by the eigenvectors of the Hamiltonian. The rate constant was obtained by integrating the flux autocorrelation function. An advantage of this procedure is that the vibrational configuration interaction method provides an intuitive picture of the mode coupling at the transion state while it also provides converged results. Furthermore, it suggests a number of well calibrated sublevels of theory that systematically approach the exact result; use of such sublevels for systematic convergence has been instrumental in the success of electronic structure theory but has so far played a lesser role in reaction dynamics.