Quantum Mechanical Pressure-Dependent Reaction and Recombination Rates for O + OH → H + O₂, HO₂

We extend recent flux−flux autocorrelation function methods for the direct computation of thermal reaction rate constants and unimolecular recombination rates to the case where both reaction and recombination are possible. Rather than a single transition state dividing surface, dividing surfaces are placed on both the reactant (r) and product (p) sides of the intermediate collision complex region. The thermal recombination rate expression then involves a flux cross-correlation function Crp(t) in addition to the usual autocorrelation function Crr(t), both of which are computed during a single quantum time propagation. This method is applied to the three-dimensional O + OH ⇌ H + O2 (J = 0) reactions, employing parallel computation because of the necessary large basis (218 grid points) and long propagation times (2−3 ps). Thermal rate constants (in the absence of recombination effects) are presented for T = 500−2000 K, using the J-shifting approximation to account for nonzero total angular momentum; good agr...