Weakly bound molecules in the atmosphere: a case study of HOOO

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

Rapid Acceleration of Hydrogen Atom Abstraction Reactions of OH at Very Low Temperatures through Weakly Bound Complexes and Tunneling

A generally accepted principle of chemical kinetics is that a reaction will be very slow at low temperatures if there is an activation barrier on the potential energy surface to form products. However, this Account shows that the reverse is true for gas-phase hydrogen abstraction reactions of the hydroxyl radical, OH, with organic molecules with which it can form a weakly bound (5-30 kJ mol^-1) hydrogen-bonded complex. For hydrogen atom abstraction reactions of OH with volatile organic compounds (VOCs) containing alcohol, ether, carbonyl, and ester functional groups, the reaction accelerates rapidly at very low temperatures, with rate coefficients, k, that can be up to a 1000 times faster than those at room temperature, despite the barrier to products. The OH radical is a crucial intermediate in Earth's atmosphere, combustion processes, and the chemistry of the interstellar medium, where temperatures can reach as low as 10 K, so this behavior has very important implications for gas-phase chemistry in space. The key point is that at low temperatures the lifetime of the OH-VOC complex against re-dissociation back to reactants becomes much longer, and hence the probability of quantum mechanical tunneling under the reaction barrier to form products becomes much higher. These observations were made possible by using Laval nozzles to generate uniform supersonic flows at extremely low temperatures so that condensation of the reagents at reactor walls is avoided. In this Account, the use of laser flash-photolysis combined with laser-induced fluorescence spectroscopy within Laval flows is described to study the unusual kinetics of this type of reaction at temperatures down to 21 K and demonstrate the rapid upturn in the rate coefficient. For the reaction of OH with CH₃OH, further evidence for the precomplex and tunneling mechanism comes from observation of the CH₃O reaction product at very low temperatures, despite it being formed over the higher barrier to reaction. The experimental observations are supported by theoretical calculations using the MESMER master equation package to calculate k( T) and product yields as a function of temperature and which make use of potential energy surfaces determined using ab initio methods. The CH₃O product is formed over a narrower barrier with a larger imaginary frequency and is calculated to be the sole product at very low temperatures. The kinetics of the OH reaction with CH₃OH were measured to be independent of pressure, consistent with a tunneling mechanism rather than any collisional stabilization of the prereactive complex. In this Account, we collate available kinetic data and show that this newly discovered mechanism for H atom transfer reactions appears to be generally applicable for reactions of OH with organic molecules containing oxygenated functional groups, which have been observed in space by radio-astronomy. Rather than being ignored for a range of interstellar environments, these OH reactions are now being included in chemical networks in space and have been shown to significantly influence the abundance of OH, the organic molecules themselves, and reaction products and provide novel routes to forming even more complex functional groups, for example, precursors to prebiotic molecules.

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.

Silica Nanoparticle-Generated ROS as a Predictor of Cellular Toxicity: Mechanistic Insights and Safety by Design

Evaluating toxicological responses of engineered nanomaterials such as silica nanoparticles is critical in assessing health risks and exposure limits. Biological assays can be used to evaluate cytotoxicity of individual materials, but specific nano-bio interactions-which govern its physiological response-cannot currently be predicted from materials characterization and physicochemical properties. Understanding the role of free radical generation from nanomaterial surfaces facilitates understanding of a potential toxicity mechanism and provides insight into how toxic effects can be assessed. Size-matched mesoporous and nonporous silica nanoparticles in aminopropyl-functionalized and native forms were investigated to analyze the effects of porosity and surface functionalization on the observed cytotoxicity. In vitro cell viability data in a murine macrophage cell line (RAW 264.7) provides a model for what might be observed in terms of cellular toxicity upon an environmental or industrial exposure to silica nanoparticles. Electron paramagnetic resonance spectroscopy was implemented to study free radical species generated from the surface of these nanomaterials and the signal intensity was correlated with cellular toxicity. In addition, in vitro assay of intracellular reactive oxygen species (ROS) matched well with both the EPR and cell viability data. Overall, spectroscopic and in vitro studies correlate well and implicate production of ROS from a surface-catalyzed reaction as a predictor of cellular toxicity. The data demonstrate that mesoporous materials are intrinsically less toxic than nonporous materials, and that surface functionalization can mitigate toxicity in nonporous materials by reducing free radical production. The broader implications are in terms of safety by design of nanomaterials, which can only be extracted by mechanistic studies such as the ones reported here.

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.

Theoretical calculation of the OH vibrational overtone spectra of 1,5-pentanediol and 1,6-hexanediol

It is well-known that intramolecular hydrogen bonding affects the relative energetics of conformers, as well as the OH stretching peak positions, intensities, and width. In this study we simulated the Δv(OH) = 3, 4 overtone spectra of 1,5-pentanediol (PeD) and 1,6-hexanediol (HD) using the peak positions, intensities, and width calculated from the B3LYP/6-31+G(d,p) method. Furthermore, room temperature free energy calculations were performed using B3LYP/6-31+G(d,p) MP2/6-31+G(d,p), and MP2/6-311++G(3df,3pd) to obtain the relative population of the conformers. From the calculation of 109 and 381 distinct conformers for PeD and HD, respectively, we find that for these long chain diols the intramolecular hydrogen bonded conformers are not the most dominant conformation at room temperature. This is in stark contrast with shorter chain diols such as ethylene glycol for which the hydrogen bonded conformer dominates the population at room temperature. On the other hand, we found that the correlation between the hydrogen bonded OH red shift versus the homogeneous width, Γ = 0.0155(Δω)(1.36), which was derived for shorter chain diols, is valid even for these longer chain diols. We also showed that the intramolecular hydrogen bonded OH initially decays through the CCOH torsion and COH bending mode no matter how long the alkanediol chain length is for 1,n-alkanediols for n up to 6.