Perspective on Theoretical and Experimental Advances in Atmospheric Photochemistry

Research that explores the chemistry of Earth's atmosphere is central to the current understanding of global challenges such as climate change, stratospheric ozone depletion, and poor air quality in urban areas. This research is a synergistic combination of three established domains: earth observation, for example, using satellites, and in situ field measurements; computer modeling of the atmosphere and its chemistry; and laboratory measurements of the properties and reactivity of gas-phase molecules and aerosol particles. The complexity of the interconnected chemical and photochemical reactions which determine the composition of the atmosphere challenges the capacity of laboratory studies to provide the spectroscopic, photochemical, and kinetic data required for computer models. Here, we consider whether predictions from computational chemistry using modern electronic structure theory and nonadiabatic dynamics simulations are becoming sufficiently accurate to supplement quantitative laboratory data for wavelength-dependent absorption cross-sections, photochemical quantum yields, and reaction rate coefficients. Drawing on presentations and discussions from the CECAM workshop on Theoretical and Experimental Advances in Atmospheric Photochemistry held in March 2024, we describe key concepts in the theory of photochemistry, survey the state-of-the-art in computational photochemistry methods, and compare their capabilities with modern experimental laboratory techniques. From such considerations, we offer a perspective on the scope of computational (photo)chemistry methods based on rigorous electronic structure theory to become a fourth core domain of research in atmospheric chemistry.

Ozonolysis can produce long-lived greenhouse gases from commercial refrigerants

Hydrofluoroolefins are being adopted as sustainable alternatives to long-lived fluorine- and chlorine-containing gases and are finding current or potential mass-market applications as refrigerants, among a myriad of other uses. Their olefinic bond affords relatively rapid reaction with hydroxyl radicals present in the atmosphere, leading to short lifetimes and proportionally small global warming potentials. However, this type of functionality also allows reaction with ozone, and whilst these reactions are slow, we show that the products of these reactions can be extremely long-lived. Our chamber measurements show that several industrially important hydrofluoroolefins produce CHF3 (fluoroform, HFC-23), a potent, long-lived greenhouse gas. When this process is accounted for in atmospheric chemical and transport modeling simulations, we find that the total radiative effect of certain compounds can be several times that of the direct radiative effect currently recommended by the World Meteorological Organization. Our supporting quantum chemical calculations indicate that a large range of exothermicity is exhibited in the initial stages of ozonolysis, which has a powerful influence on the CHF3 yield. Furthermore, we identify certain molecular configurations that preclude the formation of long-lived greenhouse gases. This demonstrates the importance of product quantification and ozonolysis kinetics in determining the overall environmental impact of hydrofluoroolefin emissions.

Gas Phase and Gas-Solid Interface Ozonolysis of Nitrogen Containing Alkenes: Nitroalkenes, Enamines, and Nitroenamines

Emerging contaminants are of concern due to their rapidly increasing numbers and potential ecological and human health effects. In this study, the synergistic effects of the presence of multifunctional nitro, amino and carbon-carbon double bond (C═C) groups on the gas phase ozonolysis in O₂ or at the air/solid interface were investigated using five simple model compounds. The gas phase ozonolysis rate constants at 296 K were (3.5 ± 0.9) × 10^-20 cm³ molecule^-1 s^-1 for 2-methyl-1-nitroprop-1-ene and (6.8 ± 0.8) × 10^-19 cm³ molecule^-1 s^-1 for 4-methyl-4-nitro-1-pentene, with lifetimes of 134 and 7 days in the presence of 100 ppb ozone in the atmosphere, respectively. The rate constants for gas phase E-N,N-dimethyl-1-propenylamine and N,N-dimethylallylamine reactions with ozone were too fast (>10^-18 cm³ molecule^-1 s^-1) to be measured, implying lifetimes of less than 5 days. A multiphase kinetics model (KM-GAP) was used to probe the gas-solid kinetics of 1-dimethylamino-2-nitroethylene, yielding a rate constant for the surface reaction of 1.8 × 10^-9 cm² molecule^-1 s^-1 and in the bulk 1× 10^-16 cm³ molecule^-1 s^-1. These results show that a nitro group attached to the C═C lowers the gas phase rate constant by 2-3 orders of magnitude compared to the simple alkenes, while amino groups have the opposite effect. The presence of both groups provides counterbalancing effects. Products with deleterious health effects including dimethylformamide and formaldehyde were identified by FTIR. The identified products differentiate whether the initial site of ozone attack is C═C and/or the amino group. This study provides a basis for predicting the environmental fates of emerging contaminants and shows that both the toxicity of both the parent compounds and the products should be taken into account in assessing their environmental impacts.

Gas-Phase Ozone Reaction Kinetics of C5-C8 Unsaturated Alcohols of Biogenic Interest

Unsaturated alcohols are volatile organic compounds (VOCs) that characterize the emissions of plants. Changes in climate together with related increases of biotic and abiotic stresses are expected to increase these emissions in the future. Ozonolysis is one of the oxidation pathways that control the fate of unsaturated alcohols in the atmosphere. The rate coefficients of the gas-phase O₃ reaction with seven C₅-C₈ unsaturated alcohols were determined at 296 K using both absolute and relative kinetic methods. The following rate coefficients (cm³ molecule^-1 s^-1) were obtained using the absolute method: (1.1 ± 0.2) × 10^-16 for cis-2-penten-1-ol, (1.2 ± 0.2) × 10^-16 for trans-2-hexen-1-ol, (6.4 ± 1.0) × 10^-17 for trans-3-hexen-1-ol, (5.8 ± 0.9) × 10^-17 for cis-3-hexen-1-ol, (2.0 ± 0.3) × 10^-17 for 1-octen-3-ol, and (8.4 ± 1.3) × 10^-17 for trans-2-octen-1-ol. The following rate coefficients (cm³ molecule^-1 s^-1) were obtained using the relative method: (1.27 ± 0.11) × 10^-16 for trans-2-hexen-1-ol, (5.01 ± 0.30) × 10^-17 for trans-3-hexen-1-ol, (4.13 ± 0.34) × 10^-17 for cis-3-hexen-1-ol, and (1.40 ± 0.12) × 10^-16 for trans-4-hexen-1-ol. Alkenols display high reactivities with ozone with lifetimes in the hour range. Rate coefficients show a strong and complex dependence on the structure of the alkenol, particularly the relative position of the OH group toward the C═C double bond. The results are discussed and compared to both the available literature data and four structure-activity relationship (SAR) methods.

Predicting the rate constants of volatile organic compounds (VOCs) with ozone reaction at different temperatures

Based on the bond order, fukui indices and other related descriptors, as well as temperature, a new QSAR model was established to predict the rate constant (kO3) of VOCs degradation by O3. 302 logkO3 values (178-409 K) of 149 VOCs were divided into training set (242 logkO3) and test set (60 logkO3), respectively, which were used to construct and verify the QSAR model. The optimal model (R2 = 0.83, q2 = 0.82, Qext2 = 0.72) shows that EHOMO, BOx and q(C-)n have a greater influence on the value of logkO3. In addition, fukui indices and logkO3 are well correlated. The applicability domains of the current models can be used to predict kO3 of a wide range of VOCs at different temperatures.

Kinetic study of the gas-phase reaction of O3 with three unsaturated alcohols

Rate constants for the reactions of ozone with 1-octen-3-ol, 1-nonen-3-ol and 1-nonen-4-ol have been determined at 298±1K and atmospheric pressure for the first time. The experiments were performed in a 100-L FEP Teflon film bag using absolute rate method; the rate constants were (1.91±0.19)×10^-17, (1.89±0.20)×10^-17, and (0.83±0.08)×10^-17cm³/(molecule·sec) for 1-octen-3-ol, 1-nonen-3-ol, and 1-nonen-4-ol, respectively. The rate constants have been compared with those of unsaturated alcohols structural homologs, and used to estimate the reaction reactivity. The electronegativity of carbon-carbon double bond was calculated by atomic charges analysis. The calculated results show that the electronic effect of the lone pair electrons of hydroxyl oxygen is the main cause of the difference in rate coefficient. According to the obtained rate constants, the atmospheric lifetimes of studied unsaturated alcohols were also estimated, which indicates that the reaction with ozone is an important loss pathway in the atmosphere, especially in polluted areas.

Direct observation of vinyl hydroperoxide

Many alkyl-substituted Criegee intermediates are predicted to undergo an intramolecular 1,4-hydrogen transfer to form isomeric vinyl hydroperoxide species (C[double bond, length as m-dash]COOH moiety), which break apart to release OH and vinoxy radicals. We report direct detection of stabilized vinyl hydroperoxides formed via carboxylic acid-catalyzed tautomerization of Criegee intermediates. A doubly hydrogen-bonded interaction between the Criegee intermediate and carboxylic acid facilitates efficient hydrogen transfer through a double hydrogen shift. Deuteration of formic or acetic acid permits migration of a D atom to yield partially deuterated vinyl hydroperoxides, which are distinguished from the CH3CHOO, (CH3)2COO, and CH3CH2CHOO Criegee intermediates by mass. Using 10.5 eV photoionization, three prototypical vinyl hydroperoxides, CH2[double bond, length as m-dash]CHOOD, CH2[double bond, length as m-dash]C(CH3)OOD, and CH3CH[double bond, length as m-dash]CHOOD, are detected directly. Complementary electronic structure calculations reveal several reaction pathways, including the barrierless acid-catalyzed tautomerization reaction predicted previously and a barrierless addition reaction that yields hydroperoxy alkyl formate.

Gas-phase rate coefficients of the reaction of ozone with four sesquiterpenes at 295 ± 2 K

Room temperature rate coefficients of the reaction of ozone with four sesquiterpenes were measured using absolute and relative rate techniques.

Development of a model for predicting reaction rate constants of organic chemicals with ozone at different temperatures

To assess the persistence and fate of volatile organic compounds in the troposphere, the rate constants for the reaction with ozone (kO3) are needed. As kO3 values are only available for hundreds of compounds, and experimental determination of kO3 is costly and time-consuming, it is of importance to develop predictive models on kO3. In this study, a total of 379 logkO3 values at different temperatures were used to develop and validate a model for the prediction of kO3, based on quantum chemical descriptors, Dragon descriptors and structural fragments. Molecular descriptors were screened by stepwise multiple linear regression, and the model was constructed by partial least-squares regression. The cross validation coefficient QCUM(2) of the model is 0.836, and the external validation coefficient Qext(2) is 0.811, indicating that the model has high robustness and good predictive performance. The most significant descriptor explaining logkO3 is the BELm2 descriptor with connectivity information weighted atomic masses. kO3 increases with increasing BELm2, and decreases with increasing ionization potential. The applicability domain of the proposed model was visualized by the Williams plot. The developed model can be used to predict kO3 at different temperatures for a wide range of organic chemicals, including alkenes, cycloalkenes, haloalkenes, alkynes, oxygen-containing compounds, nitrogen-containing compounds (except primary amines) and aromatic compounds.

Rate acceleration of the heterogeneous reaction of ozone with a model alkene at the air-ice interface at low temperatures

The kinetics of the ozonation reaction of 1,1-diphenylethylene (DPE) on the surface of ice grains (also called "artificial snow"), produced by shock-freezing of DPE aqueous solutions or DPE vapor-deposition on pure ice grains, was studied in the temperature range of 268 to 188 K. A remarkable and unexpected increase in the apparent ozonation rates with decreasing temperature was evaluated using the Langmuir-Hinshelwood and Eley-Rideal kinetic models, and by estimating the apparent specific surface area of the ice grains. We suggest that an increase of the number of surface reactive sites, and possibly higher ozone uptake coefficients are responsible for the apparent rate acceleration of DPE ozonation at the air-ice interface at lower temperatures. The increasing number of reactive sites is probably related to the fact that organic molecules are displaced more to the top of a disordered interface (or quasi-liquid) layer on the ice surface, which makes them more accessible to the gas-phase reactants. The effect of NaCl as a cocontaminant on ozonation rates was also investigated. The environmental implications of this phenomenon for natural ice/snow are discussed. DPE was selected as an example of environmentally relevant species which can react with ozone. For typical atmospheric ozone concentrations in polar areas (20 ppbv), we estimated that its half-life on the ice surface would decrease from ∼5 days at 258 K to ∼13 h at 188 K at submonolayer DPE loadings.

Development of a new flow reactor for kinetic studies. Application to the ozonolysis of a series of alkenes

A new flow reactor has been developed to study ozonolysis reactions at ambient pressure and room temperature (297 ± 2 K). The reaction kinetics of O(3) with 4-methyl-1-pentene (4M1P), 2-methyl-2-pentene (2M2P), 2,4,4-trimethyl-1-pentene (tM1P), 2,4,4-trimethyl-2-pentene (tM2P) and α-pinene have been investigated under pseudo-first-order conditions. Absolute measurements of the rate coefficients have been carried out by recording O(3) consumption in excess of organic compound. Alkene concentrations have been determined by sampling adsorbent cartridges that were thermodesorbed and analyzed by gas-chromatography coupled to flame ionization detection. Complementary experimental data have been obtained using a 250 L Teflon smog chamber. The following ozonolysis rate coefficients can be proposed (in cm(3) molecule(-1) s(-1)): k(4M1P) = (8.23 ± 0.50) × 10(-18), k(2M2P) = (4.54 ± 0.96) × 10(-16), k(tM1P) = (1.48 ± 0.11) × 10(-17), k(tM2P) = (1.25 ± 0.10) × 10(-16), and k(α-pinene) = (1.29 ± 0.16) × 10(-16), in very good agreement with literature values. The products of tM2P ozonolysis have been investigated, and branching ratios of (21.4 ± 2.8)% and (73.9 ± 7.3)% have been determined for acetone and 2,2-dimethyl-propanal, respectively. Additionally, a new nonoxidized intermediate, 2-methyl-1-propene, has been identified and quantified. A topological SAR analysis was also performed to strengthen the consistency of the kinetic data obtained with this new flow reactor.

Determination of gas-phase ozonolysis rate coefficients of a number of sesquiterpenes at elevated temperatures using the relative rate method

The rates of ozonolysis of four sesquiterpenes, β-caryophyllene, α-humulene, isolongifolene and α-cedrene, are determined in the gas phase at an elevated temperature of 366 ± 3 K and a pressure of ~780 Torr using the EXTreme RAnge chamber (EXTRA). The experimentally obtained rate coefficients agree with extrapolated room temperature rate coefficients for isolongifolene and α-cedrene but not for β-caryophyllene and α-humulene, which were found to be three orders of magnitude slower than this in the literature. These new measurements support the hypothesis that operating under ambient conditions, kinetic measurements of condensable species can be influenced adversely by heterogeneous processes and should therefore be treated with caution.

Structure-activity relationship (SAR) for the prediction of gas-phase ozonolysis rate coefficients: an extension towards heteroatomic unsaturated species

Heteroatomic unsaturated volatile organic compounds (HUVOCs) are common trace components of the atmosphere, yet their diverse chemical behaviour presents difficulties for predicting their oxidation kinetics using structure-activity relationships (SARs). An existing SAR is adapted to help meet this challenge, enabling the prediction of ozonolysis rates with unprecedented accuracy. The new SAR index, x(H), correlates strongly with available literature measurements of ozonolysis rate coefficients (R(2) = 0.87), a database representing 110 species. It was found that capturing the inductive effect rather than the steric effect is of primary importance in predicting the reactivity of these species, which is to be anticipated since HUVOCs can possess a variety of functional groups with a range of electron-withdrawing and donating tendencies. New experimental measurements of ozonolysis rate coefficients were conducted for 1-penten-3-ol, 3-methyl; ethene, 1,1-dimethoxy; E-2-pentenoic acid; E-1,2-dichloroethene; Z-1,2-dichloroethene; trichloroethene; tetrachloroethene; 1-butene, 3-chloro and 2-chloropropene, and were determined to be 5.15 × 10(-18), 4.82 × 10(-16), 3.07 × 10(-18), 8.05 × 10(-20), 4.88 × 10(-21), 6.04 × 10(-22), 1.56 × 10(-24), 2.26 × 10(-18) and 1.13 × 10(-19) cm(3) molecule(-1) s(-1), respectively. The index of the inductive effect, i(H), is compared with other indices of the electron-withdrawing capacity of a substitution, notably the Taft σ* constants and the rate of reaction of a given species with the hydroxyl radical, both of which are expected to be unaffected by steric factors. i(H) correlates strongly in both cases and suggests a universal response by olefinic species towards electrophilic addition.