Complete basis set and density functional determination of the enthalpy of formation of the controversial HO3 radical: a discrepancy between theory and experiment

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.

Weakly Bound Clusters in Astrochemistry? Millimeter and Submillimeter Spectroscopy of trans-HO3 and Comparison to Astronomical Observations

The emergence of chemical complexity during star and planet formation is largely guided by the chemistry of unstable molecules that are reaction intermediates in terrestrial chemistry. Our knowledge of these intermediates is limited by both the lack of laboratory studies and the difficulty in their astronomical detection. In this work, we focus on the weakly bound cluster HO3 as an example of the connection between laboratory spectroscopic study and astronomical observations. Here, we present a fast-sweep spectroscopic technique in the millimeter and submillimeter range to facilitate the laboratory search for trans-HO3 and DO3 transitions in a discharge supersonic jet and report their rotational spectra from 70 to 450 GHz. These new measurements enable full determination of the molecular constants of HO3 and DO3. We also present a preliminary search for trans-HO3 in 32 star-forming regions using this new spectroscopic information. HO3 is not detected, and column density upper limits are reported. This work provides additional benchmark information for computational studies of this intriguing radical, as well as a reliable set of molecular constants for extrapolation of the transition frequencies of HO3 for future astronomical observations.

Is HO3− multiple-minimum and floppy? Covalent to van der Waals isomerization and bond rupture of a peculiar anion

The HO₃⁻ anion is multiple-minimum and floppy: the two main isomers and isomerization barrier all lie quite below dissociation.

The structures of ozone and HOx radicals in aqueous solution from combined quantum/classical molecular dynamics simulations

Ozone in aqueous solution decomposes through a complex mechanism that involves initial reaction with a hydroxide ion followed by formation of a variety of oxidizing species such as HO, HO(2), and HO(3) radicals. Though a number of hydrogen-bonded complexes have been described in the gas phase, both theoretically and experimentally, the structures of ozone and HO(x) in liquid water remain uncertain. In this work, combined quantum/classical computer simulations of aqueous solutions of these species have been reported. The results show that ozone undergoes noticeable electron polarization but it does not participate in hydrogen bonds with liquid water. The main contribution of the solvation energy comes from dispersion forces. In contrast, HO(x) radicals form strong hydrogen bonds. They are better proton donors but weaker proton acceptors than water. Their electronic and geometrical structures are significantly modified by the solvent, especially in the case of HO(3). In all cases, fluctuations in amplitudes of electronic properties are considerable, suggesting that solvent effects might play a crucial role on oxidation mechanisms initiated by ozone in liquid water. These mechanisms are important in a broad range of domains, such as atmospheric processes, plant response to ambient ozone, and medical and industrial applications.

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.

Thermochemistry of the Hypobromous and Hypochlorous Acids, HOBr and HOCl

The enthalpies of formation of HOBr and HOCl have been estimated by employing coupled cluster theory in conjunction with the correlation consistent basis sets and corrections for core-valence, relativistic, and anharmonic effects. We have employed three different reactions to estimate the DeltaH(o)(f,298)(HOBr), namely, the atomization reaction and two homodesmic reactions. Our best estimation is DeltaH(o)(f,298) (HOBr) = -15.3 +/- 0.6 kcal/mol and is very likely to lie toward the more negative values. The present value is 1.4 kcal/mol lower than the widely used experimental determination of Ruscic and Berkowitz (J. Chem. Phys. 1994, 101, 7795), DeltaH(o)(f,298)(HOBr) > -13.93 +/- 0.42 kcal/mol. However, it is closer to the more recent measurement of Lock et al. (J. Phys. Chem. 1996, 100, 7972), DeltaH(o)(f,298)(HOBr) = -14.8 +/- 1 kcal/mol. In the case of HOCl we have determined DeltaH(o)(f,298)(HOCl) = -18.1 +/- 0.3 kcal/mol, just in the middle of the two experimental values proposed, -17.8 +/- 0.5 kcal/mol (JANAF), obtained from equilibrium constant measurements, and -18.36 +/- 0.03 kcal/mol (Joens, J. A. J. Phys. Chem. A 2001, 105, 11041), determined from the measurements of the Cl-OH bond energy. If our conclusions are correct, several enthalpies of formation that have been determined by experimental chemists, Orlando and Burholder (J. Phys. Chem. 1995, 99, 1143), and theoretical chemists, Lee (J. Phys. Chem. 1995, 99, 15074), need to be revised, since a larger value was used for DeltaH(o)(f,298)(HOBr). Employing the results obtained by Orlando and Burkholder for Br(2)O we propose DeltaH(o)(f,298)(Br(2)O) = 24.9 +/- 0.6 kcal/mol, and employing Lee's enthalpies of reaction we propose the following DeltaH(o)(f,298): for BrBrO, HBrO, ClOBr, ClBrO, BrClO, BrCN, BrNC, BrNO, BrON, FOBr, and FBrO, 39.5 +/- 1, 41.0 +/- 1, 22.7 +/- 1.5, 34.2 +/- 1.5, 40.9 +/- 1.5, 43.7 +/- 1.5, 80.1 +/- 1.5, 22.3 +/- 1, 46.2 +/- 1, 17.3 +/- 1.5, and 6.3 +/- 1.5 kcal/mol, respectively. We expect that this work will stimulate new experimental measurements of the thermodynamic properties of HOBr and HOCl.

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.

Kinetics of the Radical−Radical Reaction, O(3PJ) + OH(X2ΠΩ) → O2 + H, at Temperatures down to 39 K

The kinetics of the reaction between O atoms and OH radicals, both in their electronic ground state, have been investigated at temperatures down to ca. 39 K. The experiments employed a CRESU (Cinétique deRéaction en Ecoulement Supersonique Uniforme) apparatus to attain low temperatures. Both reagents were created using pulsed laser photolysis at 157.6 nm of mixtures containing H2O and O2 diluted in N2 carrier gas. OH radicals were formed by both direct photolysis of H2O and the reaction between O(1D) atoms and H2O. O(3P) atoms were formed both as a direct product of O2 photolysis and by the rapid quenching of O(1D) atoms formed in that photolysis by N2 and O2. The rates of removal of OH radicals were observed by laser-induced fluorescence, and concentrations of O atoms were estimated from a knowledge of the absorption cross-section for O2 at 157.6 nm and of the measured fluence from the F2 laser at this wavelength. To obtain a best estimate of the rate constants for the O + OH reaction, we had to correct the raw experimental data for the following: (a) the decrease in the laser fluence along the jet due to the absorption by O2 in the gas mixture, (b) the increase in temperature, and consequent decrease in gas density, as a result of energy released in the photochemical and chemical processes that occurred, and (c) the formation of OH(v = 0) as a result of relaxation, particularly by O2, of OH radicals formed in levels v > 0. Once these corrections were made, the rate constant for reaction between OH and O(3P) atoms showed little variation in the temperature range of 142 to 39 K and had a value of (3.5 +/- 1.0) x 10(-11) cm3 molecule(-1) s(-1). It is recommended that this value is used in future chemical models of dense interstellar clouds.