Analysis of the HOOO torsional potential

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.

Spectroscopic characterization of the first excited state and photochemistry of the HO3 radical

We report the one-dimensional cuts of the six-dimensional potential energy surfaces (PESs) of the ground and lowest doublet and quartet electronic states of trans-HO₃ at the MRCI-F12/aug-cc-pVTZ level of theory. Theoretical calculations predict that the first excited state (A²A) presents a real minimum on its PES and possesses a nonplanar structure. The adiabatic excitation energy at the MRCI+Q and MRCI-F12 levels shows that the A²A state lies in the near-infrared region. Both the transition dipole moment and the oscillator strength were predicted to be weak, which suggests that photodissociation of HO₃ to produce OH and O₂ after UV-Vis absorption is not a plausible mechanism. The harmonic vibrational frequencies and rotational constants of the weakly bound complex OH-O₂ in the two electronic states were predicted to help in its detection. Our PES shows that the reactions of H + O₃ or HO₂ + O in their ground states do not lead to trans-HO₃ in its ground electronic state if one of the component fragments, i.e., HO₂(A²A') + O(³P) or H(²S) + O₃(³B₂), is excited.

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₃.

Torsional Potentials of Glyoxal, Oxalyl Halides, and Their Thiocarbonyl Derivatives: Challenges for Popular Density Functional Approximations

The reliability of popular density functionals was studied for the description of torsional profiles of 36 molecules: glyoxal, oxalyl halides, and their thiocarbonyl derivatives. HF and 18 functionals of varying complexity, from local density to range-separated hybrid approximations and double-hybrid, have been considered and benchmarked against CCSD(T)-level rotational profiles. For molecules containing heavy halogens, most functionals fail to reproduce barrier heights accurately and a number of functionals introduce spurious minima. Dispersion corrections show no improvement. Calibrated torsion-corrected atom-centered potentials rectify the shortcomings of PBE and also improve on σ-hole based intermolecular binding in dimers and crystals.

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.