A Theoretical Study of the Potential Energy Surface and Rate Constant for an O(3P) + HO2 Reaction

Reaction Pathway Control via Reactant Vibrational Excitation and Impact on Product Vibrational Distributions: The O + HO2 → OH + O2 Atmospheric Reaction

Chemical reactions often have multiple pathways, the control of which is of fundamental and practical importance. In this Letter, we examine the dynamics of the O + HO₂ → OH + O₂ reaction, which plays an important role in atmospheric chemistry, using quasi-classical trajectories on a recently developed full-dimensional potential energy surface (PES). This reaction has two pathways leading to the same products: the H abstraction pathway (O_a + HO_bO_c → O_aH + O_bO_c) and the O abstraction pathway (O_a + HO_bO_c → O_bH + O_aO_c). Under thermal conditions, the reaction is dominated by the latter channel, which is barrierless, leading to vibrational excitation of the O₂ product. However, we demonstrate that excitation of the HO₂ reactant in its O-H (v₁) vibrational mode results in dramatic switching of the reaction pathway to the activated H abstraction channel, which leads to a highly excited OH product vibrational state distribution. The implications of such dynamical effects in the atmospheric chemistry are discussed.

Theoretical H + O3 rate coefficients from ring polymer molecular dynamics on an accurate global potential energy surface: assessing experimental uncertainties

Thermal rate coefficients and kinetic isotope effects have been calculated for an important atmospheric reaction H/D + O3 → OH/OD + O2 based on an accurate permutation invariant polynomial-neural network potential energy surface, using ring polymer molecular dynamics (RPMD), quasi-classical trajectory (QCT) and variational transition-state theory (VTST) with multidimensional tunneling. The RPMD approach yielded results that are generally in better agreement with experimental rate coefficients than the VTST and QCT ones, especially at low temperatures, attributable to its capacity to capture quantum effects such as tunneling and zero-point energy. The theoretical results support one group of existing experiments over the other. In addition, rate coefficients for the D + O3 → OD + O2 reaction are also reported using the same methods, which will allow a stringent assessment of future experimental measurements, thus helping to reduce the uncertainty in the recommended rate coefficients of this reaction.

Riddles of the structure and vibrational dynamics of HO3 resolved near the ab initio limit

The hydridotrioxygen (HO₃) radical has been investigated in many previous theoretical and experimental studies over several decades, originally because of its possible relevance to the tropospheric HO_x cycle but more recently because of its fascinating chemical bonding, geometric structure, and vibrational dynamics. We have executed new, comprehensive research on this vexing molecule via focal point analyses (FPA) to approach the ab initio limit of optimized geometric structures, relative energies, complete quartic force fields, and the entire reaction path for cis-trans isomerization. High-order coupled cluster theory was applied through the CCSDT(Q) and even CCSDTQ(P) levels, and CBS extrapolations were performed using cc-pVXZ (X = 2-6) basis sets. The cis isomer proves to be higher than trans by 0.52 kcal mol^-1, but this energetic ordering is achieved only after the CCSDT(Q) milestone is reached; the barrier for cis → trans isomerization is a minute 0.27 kcal mol^-1. The FPA central r_e(O-O) bond length of trans-HO₃ is astonishingly long (1.670 Å), consistent with the semiexperimental r_e distance we extracted from microwave rotational constants of 10 isotopologues using FPA vibration-rotation interaction constants (α_i). The D₀(HO-O₂) dissociation energy converges to a mere 2.80 ± 0.25 kcal mol^-1. Contrary to expectation for such a weakly bound system, vibrational perturbation theory performs remarkably well with the FPA anharmonic force fields, even for the torsional fundamental near 130 cm^-1. Exact numerical procedures are applied to the potential energy function for the torsional reaction path to obtain energy levels, tunneling rates, and radiative lifetimes. The cis → trans isomerization occurs via tunneling with an inherent half-life of 1.4 × 10^-11 s and 8.6 × 10^-10 s for HO₃ and DO₃, respectively, thus resolving the mystery of why the cis species has not been observed in previous experiments executed in dissipative environments that allow collisional cooling of the trans-HO₃ product. In contrast, the pure ground eigenstate of the cis species in a vacuum is predicted to have a spontaneous radiative lifetime of about 1 h and 5 days for HO₃ and DO₃, respectively.

Anab initiobased full-dimensional potential energy surface for OH + O2⇄ HO3and low-lying vibrational levels of HO3

A full-dimensional potential energy surface for HO₃, including the HO + O₂dissociation asymptote, is developed and rigorous quantum dynamics calculations based on this PES have been carried out to compute the vibrational energy levels of HO₃.

An Evaluation of Gas Phase Enthalpies of Formation for Hydrogen-Oxygen (HxOy) Species

We have compiled gas phase enthalpies of formation for nine hydrogen-oxygen species (HxOy) and selected recommended values for H, O, OH, H2O, HO2, H2O2, O3, HO3, and H2O3. The compilation consists of values derived from experimental measurements, quantum chemical calculations, and prior evaluations. This work updates the recommended values in the NIST-JANAF (1985) and Gurvich et al. (1989) thermochemical tables for seven species. For two species, HO3 and H2O3 (important in atmospheric chemistry) and not found in prior thermochemical evaluations, we also provide supplementary data consisting of molecular geometries, vibrational frequencies, and torsional potentials which can be used to compute thermochemical functions. For all species, we also provide supplementary data consisting of zero point energies, vibrational frequencies, and ion reaction energetics.

Methane activation by metal-free radical cations: experimental insight into the reaction intermediate

A precise jab to methane: The SO(2)(*+) radical cation (see figure) effectively activates CH(4) at room temperature through a [H(3)C(*)...HOSO(+)] methyl intermediate isolated in the gas phase by mass spectrometry. Methanol and ionized methyl hydrogen sulfoxylate, CH(3)OSOH(*+), are formed by selective, direct attack of the incipient methyl radical at the O atom of the intermediate. The reaction shows radical and charge effects in the activation of methane by metal-free radical cations.

The Gas-Phase Reaction between O3 and HO Radical: A Theoretical Study

We report a theoretical study on the reaction of ozone with hydroxyl radical, which is important in the chemistry of the atmosphere and in particular participates in stratospheric ozone destruction. The reaction is a complex process that involves, in the first stage, a pre-reactive hydrogen-bonded complex (C1), which is formed previous to two transition states (TS1 and TS2) involving the addition of the hydroxyl radical to ozone, and leads to the formation of HO4 polyoxide radical before the release of the products HO2 and O2. The reaction is computed to be exothermic by 42.72 kcal mol(-1), which compares quite well with the experimental estimate, and the energy barriers of TS1 and TS2 with respect to C1 are computed to be 1.80 and 2.26 kcal mol(-1) at 0 K. A kinetic study based on the variational transition state theory (VTST) predicts a rate constant, at 298 K, of 7.37 x 10(-14) cm3 molecule(-1) s(-1), compared to the experimentally recommended value of 7.25 x 10(-14) cm3 molecule(-1) s(-1).

Stability of the Hydrogen Trioxy Radical via Infrared Action Spectroscopy

The hydrogen trioxy radical (HO3) has been proposed as an intermediate in several important chemical reactions and relaxation processes involving OH in the atmosphere. In this work, the gas-phase infrared action spectrum of HO3 is obtained in the OH overtone region, along with the product state distribution of the OH fragment following dissociation. The highest observed OH product channel sets an upper limit for the HO-O2 binding energy of 6.12 kcal mol(-1). The experimental stability of HO3 and derived equilibrium constant imply that up to 66% of atmospheric OH may be converted into HO3 in the tropopause region.

Ab initio study on the decomposition of first excited state HOOO radicals

The HOOO radical plays a crucial role in atmospheric processes involving the OH radical and O(2) molecule. We present an ab initio molecular orbital theory study on the decomposition reaction of the first excited state HOOO((2)A') with respect to OH and O(2). The geometries and harmonic vibrational frequencies of all stationary points are calculated at the CASSCF and MRCI levels of theory in conjunction with the 6-31+G(d,p) basis set. The potential energy profile of the decomposition reaction is studied at the CASSCF/6-31+G(d,p) level of theory, in which the complete valence orbitals and electrons are included in the active space. The energies of the potential energy profile are further refined at the CASPT2 and MRCI levels of the theory. Additionally, we have determined the interesting reaction process: the HOOO((2)A') radical with C(s) symmetry does not dissociate to OH((2)Pi) and O(2)((3)Sigma(-)(g)) directly as this is forbidden by orbital symmetry, but dissociates to OH((2)Pi) and O(2)((3)Sigma(-)(g)) via the change in symmetry from C(s) to C(infinity v) symmetry with a low barrier.

Structure of the HOOO Radical in Liquid Water: A Theoretical Study

The HOOO radical is supposed to play a role in ozone chemistry, both in the gas phase and aqueous media. We discuss the influence of the solvent on the electronic and geometrical structure of this radical using density functional and high-level ab initio calculations together with continuum, discrete, and discrete-continuum solvent models. Solute-solvent electrostatic interactions are shown to be fundamental, and lead to a noticeable stabilization of the radical, which should adopt a trans conformation in aqueous media. In fact, no energy minimum for the cis conformation is predicted in these conditions.