Quasiclassical trajectory calculations of the OH+NO2 association reaction

PHYS 576

Chao Chen, cchen7@emory.edu1, Benjamin C. Shepler, bcshepl@emory.edu1, Bastiaan J. Braams2, and Joel M. Bowman, jmbowma@emory.edu3. (1) Department of Chemistry, Emory University, atlanta, GA 30322, (2) Cherry L. Emerson Center of Scientific Computation and Department of Chemistry, Emory University, Atlanta, GA 30322, (3) Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, GA 30322
We report a full-dimensional potential energy surface (PES) for the OH+NO2 reaction based on fitting more than 55,000 energies obtained with Density Functional Theory-B3LYP/6-311G(d,p)-calculations. The PES is invariant with respect to permutation of like nuclei and describes all isomers of HOONO, HONO2, and the fragments OH+NO2 and HO2+NO. Detailed comparison of the structures, energies and harmonic frequencies of various stationary points on the PES are made with previous and present high-level ab initio calculations. Two hydrogen-bond complexes are found on the PES and confirmed by new ab initio CASPT2 calculations. Quasiclassical trajectory calculations of the cross sections for ground ro-vibrational OH+NO2 association reactions to form HOONO and HONO2 are done using this PES. The cross section to form HOONO is larger than the one to form HONO2 at low collision energies but the reverse is found at higher energies. The enhancement of the HOONO complex at low collision energies is shown to be due, in large part, to the transient formation an H-bond complex, which decays preferentially to HOONO. The association cross sections are used to obtain rate constants for formation of HOONO and HONO2 for the ground ro-vibrational states in the high-pressure limit.