Atmospheric Chemistry of Criegee Intermediates: Unimolecular Reactions and Reactions with Water

Ozonolysis of ketoprofen in polluted water: Reaction pathways, kinetics, removal efficiency, and health effects

Ketoprofen (KET), as a non-steroidal anti-inflammatory drug frequently detected in aqueous environments, is a threat to human health due to its accumulation and low biodegradability, which requires the transformation and degradation of KET in aqueous environments. In this paper, the reaction process of ozone-initiated KET degradation in water was investigated using density functional theory (DFT) method at the M06-2X/6-311++g(3df,2p)//M06-2X/6-31+g(d,p) level. The detailed reaction path of KET ozonation is proposed. The thermodynamic results show that ozone-initiated KET degradation is feasible. Under ultraviolet irradiation, the reaction of ozone with water can also produce OH radicals (HO·) that can react with KET. The degradation reaction of KET caused by HO· was further studied. The kinetic calculation illustrates that the reaction rate (1.99 × 10-1 (mol/L)-1 sec-1) of KET ozonation is relatively slow, but the reaction rate of HO· reaction is relatively high, which can further improve the degradation efficiency. On this basis, the effects of pollutant concentration, ozone concentration, natural organic matter, and pH value on degradation efficiency under UV/O3 process were analyzed. The ozonolysis reaction of KET is not sensitive to pH and is basically unaffected. Finally, the toxicity prediction of oxidation compounds produced by degradation reaction indicates that most of the degradation products are harmless, and a few products containing benzene rings are still toxic and have to be concerned. This study serves as a theoretical basis for analyzing the migration and transformation process of anti-inflammatory compounds in the water environment.

Gas-phase and air-solid interface behavior of phthalate plasticizer and ozone: The influence of indoor mineral dust

Indoor environments serve as reservoirs for a variety of emerging pollutants (EPs), such as phthalates (PAE), with intricate interactions occurring between these compounds and indoor oxidants alongside dust particles. However, the precise mechanisms governing these interactions and their resulting environmental implications remain unclear. By theoretical simulations, this work uncovers multi-functional compounds and high oxygen molecules as important products arising from the interaction between DEP/DEHP and O3, which are closely linked to SOA formation. Further analysis reveals a strong affinity of DEP/DEHP for mineral dust surfaces, with an adsorption energy of 22.11/30.91 kcal mol-1, consistent with a higher concentration of DEHP on the dust surface. Importantly, mineral particles are found to inhibit every step of the reaction process, albeit resulting in lower product toxicity compared to the parent compounds. Thus, timely removal of dust in an indoor environment may reduce the accumulation and residue of PAEs indoors, and further reduce the combined exposure risk produced by PAEs-dust. This study aims to enhance our understanding of the interaction between PAEs and SOA formation, and to develop a fundamental reaction model at the air-solid surface, thereby shedding light on the microscopic behaviors and pollution mechanisms of phthalates on indoor dust surfaces.

Theoretical investigation of the OH-initiated atmospheric degradation mechanism of CX2CHX (X = H, F, Cl) by advanced quantum chemical and transition state theory methods

Halogenated olefins are anthropogenic compounds with many industrial applications but at the same time raising many environmental and health concerns. Gas-phase electrophilic addition of the OH radical to the olefinic CC bond represents the primary sink for these chemicals in the atmosphere, with the degree and type of halogenation playing a significant role in their overall reactivity. In this work, we present a theoretical investigation of the reaction mechanisms and kinetics for the reactions between the OH radical and CH2CH2 (ethylene, ETH), CF2CHF (trifluoroethylene, TFE) and CCl2CHCl (trichloroethylene, TCE), simulated by state-of-the-art protocols and methods, with the aim of providing a detailed interpretation of the available experimental results, as well as new data of relevance to tropospheric chemistry. Specifically, potential energy surfaces (PESs) are obtained using the jun-Cheap (jChS) composite scheme, whereas temperature and pressure dependent rate coefficients and product distributions in the 100-600 K temperature range are calculated within the Rice-Ramsperger-Kassel-Marcus/master equation (RRKM/ME) framework. The rates for barrierless channels are obtained from variable reaction coordinate-variational transition state theory (VRC-VTST) combined with the two transition state model. While the reactions with ETH and TFE proceed mainly via the formation of addition adducts at P = 1 atm and T = 298 K, the dominant channel for TCE is the Cl-elimination reaction. Global rate constants for the two halogenated olefins, TFE and TCE, are found to be pressure-independent, contrary to the case of ETH. The computed rate constants, as well as their temperature and pressure dependence, are in remarkable agreement with the available experimental data, and they are used to derive atmospheric lifetimes (τ) for both TFE and TCE as a function of altitude (h) in the atmosphere, by taking into account variations in the rate coefficients (k (T, P)) and [OH] concentration.

Insights into the catalytic effect of atmospheric organic trace species on the hydration of Criegee intermediates

The bimolecular reactions between Criegee intermediates (CIs) and atmospheric trace species have been extensively investigated, with a particular focus on the reaction with water, while the catalytic role of atmospheric organic compounds in hydration reactions was often neglected. In this study, we employed quantum chemical calculations and Born-Oppenheimer molecular dynamics (BOMD) simulations to investigate the catalytic effects of atmospheric organic amines, organic acids, and alcohols on the hydration reactions of CIs in the gas phase and at the gas-liquid interface. The catalytic reactions were found to follow a cyclic catalytic structure and a stepwise reaction mechanism. Gas-phase studies revealed that organic acids exhibited stronger catalytic effects compared to amines and alcohols, and the catalytic efficiency of amines and alcohols was similar to those of single water molecule. In addition, the catalytic reaction barriers of organic acids and alcohols were positively correlated with their gas-phase acidity (R2 = 0.94 to 0.97). A negative correlation was observed between the catalytic reaction barrier of amines and their gas-phase basicity (R2 = 0.84 to 0.90) and proton affinity (R2 = 0.84 to 0.92). At the gas-liquid interface, organic acids promoted the formation of hydroxyethyl hydroperoxide (HEHP, CH3CH(OH)(OOH)), organic acid ions, and H3O+, whereas the catalytic hydration of CIs by organic amines resulted in the formation of CH3CH(OH)OO and amine ions. Both HEHP and CH3CH(OH)OO can be further decomposed to form OH and HO2, or participate in new particles formation as precursors. This study complements the research gap on the reaction of CIs with water, providing valuable insights into the atmospheric sources of HEHP and HOx as well as the formation of secondary organic aerosols (SOAs).

How Do the Position and Number of Methyl Substituents Affect the Photochemical Process of Criegee Intermediate? Trajectory Surface-Hopping Dynamics of Four-Carbon CIs

Electronic-structure calculations combined with nonadiabatic trajectory surface-hopping (TSH) dynamic simulations were carried out on two alkenyl-substituted Criegee intermediates (CIs), i.e., propenyl-substituted CI (PCI) and 1-methyl-propenyl substituted CI (MPCI), in order to investigate the influence of the position and number of substituents on the photochemical process of CI in S1 states. It is found that they play critical roles in the reactivity, dominant product channel, and mechanism of the CIs. More specifically, introducing a methyl group on either C1 (α-C) or C3 (γ-C) position of a vinyl-substituted CI (VCI) skeleton facilitates the rotation of the C1═O1 bond and leads to the formation of a three-membered dioxirane ring; meanwhile, it evidently enhances the reactively of the S1-state molecule. Meanwhile, methyl substitution on the vinyl moiety [i.e., C2 (β-C) and C3 (γ-C) positions] is beneficial for the rotation of the C2═C3 bond and thus facilitates the formation of the five-membered 1,2-dioxole ring, and the substitution on C2 site decreases the reactivity. The cosubstitution of C2 and C3 atoms by methyl groups well balances the features of VCI in the sense of high reactivity, consistently predominant channel, and possible dioxole side-product. The findings here not only deepen the knowledge on the photochemical processes of the CI but also inspire the rethinking of the "old" concept of substitution effect.

Ab initio MO study on direct production of H2O, N2O and CO3 from the respective CH2OO "Bee-sting-like" attack at H2, N2 and CO2

CONTEXT

We have computationally elucidated the mechanism for formation of H2O, N2O and CO3 from the reactions of CH2OO with H2, N2 and CO2, respectively, by the direct attack of the terminal O atom of CH2OO. This unique mechanism, which is characteristically "bee-sting-like" in nature, was found to be closely parallel to their reactions with the O(1D) atom. Reactions with H2 and CO2 take place by side-on attack, while the N2 reaction occurs by end-on attack with predicted barriers, 19.4, 13.1 and 25.3 kcal.mol-1, respectively. The CO2 reaction with CH2OO was found to occur by producing the C2v CO3, O = C < (O)O, instead of its D3h conformer, essentially similar to the O(1D) + CO2 reaction. The rate constants for the three reactions have been computed by the transition state theory (TST) based on the predicted potential energy profiles. We have also utilized the isodesmic nature of the dative bond exchange in the N2 reaction, CH2O → O + N2 = CH2O + N2 → O, to estimate the heat of the formation of CH2OO. Based on the heat of reaction computed at the highest level of theory employed, we obtained ΔfHo0 (CH2OO) = 27.5 kcal.mol-1; the value agrees with the recent results within ± 1 kcal.mol-1.

METHODS

All calculations were performed using Gaussian 16 software. Geometry, frequency, and IRC analysis calculations were conducted at the M06-2X/aug-cc-pVTZ level of theory. The heats of reaction have been evaluated at the highest level, CCSD(T)/CBS(T,Q,5)//M06-2x/aug-cc-pvTz.

Revealing new pathways for the reaction of Criegee intermediate CH2OO with SO2

Criegee intermediates play an important role in the tropospheric oxidation models through their reactions with atmospheric trace chemicals. We develop a global full-dimensional potential energy surface for the CH2OO + SO2 system and reveal how the reaction happens step by step by quasi-classical trajectory simulations. A new pathway forming the main products (CH2O + SO3) and a new product channel (CO2 + H2 + SO2) are predicted in our simulations. The new pathway appears at collision energies greater than 10 kcal/mol whose behavior demonstrates a typical barrier-controlled reaction. This threshold is also consistent with the ab initio transition state barrier height. For the minor products, a loose complex OCH2O ∙ ∙ ∙ SO2 is formed first, and then in most cases it soon turns into HCOOH + SO2, in a few cases it decomposes into CO2 + H2 + SO2 which is a new product channel, and rarely it remains as ∙OCH2O ∙ + SO2.

Tracking the reaction networks of acetaldehyde oxide and glyoxal oxide Criegee intermediates in the ozone-assisted oxidation reaction of crotonaldehyde

The reaction of unsaturated compounds with ozone (O3) is recognized to lead to the formation of Criegee intermediates (CIs), which play a key role in controlling the atmospheric budget of hydroxyl radicals and secondary organic aerosols. The reaction network of two CIs with different functionality, i.e. acetaldehyde oxide (CH3CHOO) and glyoxal oxide (CHOCHOO) formed in the ozone-assisted oxidation reaction of crotanaldehyde (CA), is investigated over a temperature range between 390 K and 840 K in an atmospheric pressure jet-stirred reactor (JSR) at a residence time of 1.3 s, stoichiometry of 0.5 with a mixture of 1% crotonaldehyde, 10% O2, at an fixed ozone concentration of 1000 ppm and 89% Ar dilution. Molecular-beam mass spectrometry in conjunction with single photon tunable synchrotron vacuum-ultraviolet (VUV) radiation is used to identify elusive intermediates by means of experimental photoionization energy scans and ab initio threshold energy calculations for isomer identification. Addition of ozone (1000 ppm) is observed to trigger the oxidation of CA already at 390 K, which is below the temperature where the oxidation reaction of CA was observed in the absence of ozone. The observed CA + O3 product, C4H6O4, is found to be linked to a ketohydroperoxide (2-hydroperoxy-3-oxobutanal) resulting from the isomerization of the primary ozonide. Products corresponding to the CIs uni- and bi-molecular reactions were observed and identified. A network of CI reactions is identified in the temperature region below 600 K, characterized by CIs bimolecular reactions with species like aldehydes, i.e., formaldehyde, acetaldehyde, and crotonaldehyde and alkenes, i.e., ethene and propene. The region below 600 K is also characterized by the formation of important amounts of typical low-temperature oxidation products, such as hydrogen peroxide (H2O2), methyl hydroperoxide (CH3OOH), and ethyl hydroperoxide (C2H5OOH). Detection of additional oxygenated species such as alcohols, ketene, and aldehydes are indicative of multiple active oxidation routes. This study provides important information about the initial step involved in the CIs assisted oligomerization reactions in complex reactive environments where CIs with different functionalities are reacting simultaneously. It provides new mechanistic insights into ozone-assisted oxidation reactions of unsaturated aldehydes, which is critical for the development of improved atmospheric and combustion kinetics models.

Evaluating Possible Formation Mechanisms of Criegee Intermediates during the Heterogeneous Autoxidation of Squalene

Organic molecules in the environment oxidatively degrade by a variety of free radical, microbial, and biogeochemical pathways. A significant pathway is heterogeneous autoxidation, in which degradation occurs via a network of carbon and oxygen centered free radicals. Recently, we found evidence for a new heterogeneous autoxidation mechanism of squalene that is initiated by hydroxyl (OH) radical addition to a carbon-carbon double bond and apparently propagated through pathways involving Criegee Intermediates (CI) produced from β-hydroxy peroxy radicals (β-OH-RO2•). It remains unclear, however, exactly how CI are formed from β-OH-RO2•, which could occur by a unimolecular or bimolecular pathway. Combining kinetic models and multiphase OH oxidation measurements of squalene, we evaluate the kinetic viability of three mechanistic scenarios. Scenario 1 assumes that CI are formed by the unimolecular bond scission of β-OH-RO2•, whereas Scenarios 2 and 3 test bimolecular pathways of β-OH-RO2• to yield CI. Scenario 1 best replicates the entire experimental data set, which includes effective uptake coefficients vs [OH] as well as the formation kinetics of the major products (i.e., aldehydes and secondary ozonides). Although the unimolecular pathway appears to be kinetically viable, future high-level theory is needed to fully explain the mechanistic relationship between CI and β-OH-RO2• in the condensed phase.

Reaction Kinetics of CH2OO and syn-CH3CHOO Criegee Intermediates with Acetaldehyde

Criegee intermediates exert a crucial influence on atmospheric chemistry, functioning as powerful oxidants that facilitate the degradation of pollutants, and understanding their reaction kinetics is essential for accurate atmospheric modeling. In this study, the kinetics of CH2OO and syn-CH3CHOO reactions with acetaldehyde (CH3CHO) were investigated using a flash photolysis reaction tube coupled with the OH laser-induced fluorescence (LIF) method. The experimental results indicate that the reaction of syn-CH3CHOO with CH3CHO is independent of pressure in the range of 5-50 Torr when using Ar as the bath gas. However, the rate coefficient for the reaction between CH2OO and CH3CHO at 5.5 Torr was found to be lower compared to the near-constant values observed between 10 and 100 Torr. Furthermore, the reaction of syn-CH3CHOO with CH3CHO demonstrated positive temperature dependence from 283 to 330 K, with a rate coefficient of (2.11 ± 0.45) × 10-13 cm3 molecule-1 s-1 at 298 K. The activation energy and pre-exponential factor derived from the Arrhenius plot for this reaction were determined to be 2.32 ± 0.49 kcal mol-1 and (1.66 ± 0.61) × 10-11 cm3 molecule-1 s-1, respectively. In comparison, the reaction of CH2OO with CH3CHO exhibited negative temperature dependence, with a rate coefficient of (2.16 ± 0.39) × 10-12 cm3 molecule-1 s-1 at 100 Torr and 298 K and an activation energy and a pre-exponential factor of -1.73 ± 0.31 kcal mol-1 and (1.15 ± 0.21) × 10-13 cm3 molecule-1 s-1, respectively, over the temperature range of 280-333 K.

Quantitative kinetics reveal that reactions of HO2 are a significant sink for aldehydes in the atmosphere and may initiate the formation of highly oxygenated molecules via autoxidation

Large aldehydes are widespread in the atmosphere and their oxidation leads to secondary organic aerosols. The current understanding of their chemical transformation processes is limited to hydroxyl radical (OH) oxidation during daytime and nitrate radical (NO3) oxidation during nighttime. Here, we report quantitative kinetics calculations of the reactions of hexanal (C5H11CHO), pentanal (C4H9CHO), and butanal (C3H7CHO) with hydroperoxyl radical (HO2) at atmospheric temperatures and pressures. We find that neither tunneling nor multistructural torsion anharmonicity should be neglected in computing these rate constants; strong anharmonicity at the transition states is also important. We find rate constants for the three reactions in the range 3.2-7.7 × 10-14 cm3 molecule-1 s-1 at 298 K and 1 atm, showing that the HO2 reactions can be competitive with OH and NO3 oxidation under some conditions relevant to the atmosphere. Our findings reveal that HO2-initiated oxidation of large aldehydes may be responsible for the formation of highly oxygenated molecules via autoxidation. We augment the theoretic studies with laboratory flow-tube experiments using an iodide-adduct time-of-flight chemical ionization mass spectrometer to confirm the theoretical predictions of peroxy radicals and the autoxidation pathway. We find that the adduct from HO2 + C5H11CHO undergoes a fast unimolecular 1,7-hydrogen shift with a rate constant of 0.45 s-1. We suggest that the HO2 reactions make significant contributions to the sink of aldehydes.

Infrared spectroscopy of the syn-methyl-substituted Criegee intermediate: A combined experimental and theoretical study

An IR-vacuum ultraviolet (VUV) ion-dip spectroscopy method is utilized to examine the IR spectrum of acetaldehyde oxide (CH3CHOO) in the overtone CH stretch (2νCH) spectral region. IR activation creates a depletion of the ground state population that reduces the VUV photoionization signal on the parent mass channel. IR activation of the more stable and populated syn-CH3CHOO conformer results in rapid unimolecular decay to OH + vinoxy products and makes the most significant contribution to the observed spectrum. The resultant IR-VUV ion-dip spectrum of CH3CHOO is similar to that obtained previously for syn-CH3CHOO using IR action spectroscopy with UV laser-induced fluorescence detection of OH products. The prominent IR features at 5984 and 6081 cm-1 are also observed using UV + VUV photoionization of OH products. Complementary theoretical calculations utilizing a general implementation of second-order vibrational perturbation theory provide new insights on the vibrational transitions that give rise to the experimental spectrum in the overtone CH stretch region. The introduction of physically motivated small shifts of the harmonic frequencies yields remarkably improved agreement between experiment and theory in the overtone CH stretch region. The prominent features are assigned as highly mixed states with contributions from two quanta of CH stretch and/or a combination of CH stretch with an overtone in mode 4. The generality of this approach is demonstrated by applying it to three different levels of electronic structure theory/basis sets, all of which provide spectra that are virtually indistinguishable despite showing large deviations prior to introducing the shifts to the harmonic frequencies.

Pressure and temperature dependent kinetics and the reaction mechanism of Criegee intermediates with vinyl alcohol: a theoretical study

Criegee intermediates (CIs), the key intermediates in the ozonolysis of olefins in atmosphere, have received much attention due to their high activity. The reaction mechanism of the most simple Criegee intermediate CH2OO with vinyl alcohol (VA) was investigated by using the HL//M06-2X/def2TZVP method. The temperature and pressure dependent rate constant and product branching ratio were calculated using the master equation method. For CH2OO + syn-VA, 1,4-insertion is the main reaction channel while for the CH2OO + anti-VA, cycloaddition and 1,2-insertion into the O-H bond are more favorable than the 1,4-insertion reaction. The 1,4-insertion or cycloaddition intermediates are stabilized collisionally at 300 K and 760 torr, and the dissociation products involving OH are formed at higher temperature and lower pressure. The rate constants of the CH2OO reaction with syn-VA and anti-VA both show negative temperature effects, and they are 2.95 × 10-11 and 2.07 × 10-13 cm3 molecule-1 s-1 at 300 K, respectively, and the former is agreement with the prediction in the literature.

Experimental and theoretical study of Criegee intermediate (CH2OO) reactions with n-butyraldehyde and isobutyraldehyde: kinetics, implications and atmospheric fate

The reactions of the simplest Criegee intermediate (CH2OO) with n-butyraldehyde (nBD) and isobutyraldehyde (iBD) were studied at 253-318 K and (50 ± 2) torr, using Cavity Ring-down spectroscopy (CRDS). The rate coefficients obtained at room temperature were (2.63 ± 0.14) × 10-12 and (2.20 ± 0.21) × 10-12 cm3 molecule-1 s-1 for nBD and iBD, respectively. Both the reactions show negative temperature-dependency, following equations, knBD(T = 253-318 K) = (11.51 ± 4.33) × 10-14 × exp{(918.1 ± 107.2)/T} and kiBD(T = 253-318 K) = (6.23 ± 2.29) × 10-14 × exp{(1051.4 ± 105.2)/T} cm3 molecule-1 s-1. High-pressure limit rate coefficients were determined from theoretical calculations at the CCSD(T)-F12/cc-pVTZ-F12//B3LYP/6-311+G(2df, 2p) level of theory, with <40% deviation from the experimental results at room temperature and above. The kinetic simulations were performed using a master equation solver to predict the temperature-dependency of the rate coefficients at the experimental pressure, as well as to predict the contribution of individual pathways. The major products predicted from the theoretical calculations were formaldehyde and formic acid, along with butyric acid from nBD and isobutyric acid from iBD reactions.

Rapid Atmospheric Reactions between Criegee Intermediates and Hypochlorous Acid

Hypochlorous acid (HOCl) is a paramount compound in the atmosphere due to its significant contribution to both tropospheric oxidation capacity and ozone depletion. The main removal routes for HOCl are photolysis and the reaction with OH during the daytime, while these processes are unimportant during the nighttime. Here, we report the rapid reactions of Criegee intermediates (CH2OO and anti/syn-CH3CHOO) with HOCl by using high-level quantum chemical methods as the benchmark; their accuracy is close to coupled cluster theory with single, double, and triple excitations and quasiperturbative connected quadruple excitations with a complete basis limit by extrapolation [CCSDT(Q)/CBS]. Their rate constants have been calculated by using a dual-level strategy; this combines conventional transition state theory calculated at the benchmark level with variational transition state theory with small-curvature tunneling by a validated density functional method. We find that the introduction of the methyl group into Criegee intermediates not only affects their reactivities but also exerts a remarkable influence on anharmonicity. The calculated results uncover that anharmonicity increases the rate constants of CH2OO + HOCl by a factor of 18-5, while it is of minor importance in the anti/syn-CH3CHOO + HOCl reaction at 190-350 K. The present findings reveal that the loose transition state for anti-CH3CHOO and HOCl is a rate-determining step at 190-350 K. We also find that the reaction of Criegee intermediates with HOCl contributes significantly to the sink of HOCl during the nighttime in the atmosphere.

Ammonolysis of Glyoxal at the Air-Water Nanodroplet Interface

The reactions of glyoxal (CHO)2 ) with amines in cloud processes contribute to the formation of brown carbon and oligomer particles in the atmosphere. However, their molecular mechanisms remain unknown. Herein, we investigate the ammonolysis mechanisms of glyoxal with amines at the air-water nanodroplet interface. We identified three and two distinct pathways for the ammonolysis of glyoxal with dimethylamine and methylamine by using metadynamics simulations at the air-water nanodroplet interface, respectively. Notably, the stepwise pathways mediated by the water dimer for the reactions of glyoxal with dimethylamine and methylamine display the lowest free energy barriers of 3.6 and 4.9 kcal ⋅ mol-1 , respectively. These results showed that the air-water nanodroplet ammonolysis reactions of glyoxal with dimethylamine and methylamine were more feasible and occurred at faster rates than the corresponding gas phase ammonolysis, the OH+(CHO)2 reaction, and the aqueous phase reaction of glyoxal, leading to the dominant removal of glyoxal. Our results provide new and important insight into the reactions between carbonyl compounds and amines, which are crucial in forming nitrogen-containing aerosol particles.

Quantitative kinetics of the atmospheric reaction between isocyanic acid and hydroxyl radicals: post-CCSD(T) contribution, anharmonicity, recrossing effects, torsional anharmonicity, and tunneling

Hydroxyl radicals (OH) are the most important atmospheric oxidant, initiating atmospheric reactions for the chemical transformation of volatile organic compounds. Here, we choose the HNCO + OH reaction as a prototype reaction because it contains the fundamental reaction processes for OH radicals: H-abstraction reaction by OH and OH addition reaction. However, its kinetics are unknown under atmospheric conditions. We investigate the reaction of HNCO with OH by using the GMM(P).L method close to the accuracy of single, double, triple, and quadruple excitations and noniterative quintuple excitations with a complete basis set (CCSDTQ(P)/CBS) as benchmark results and a dual-level strategy for kinetics calculations. The calculated rate constant of HNCO + OH is in good agreement with the experimental data available at the temperatures between 620 and 2500 K. We find that the rate constant cannot be correctly obtained by using experimental data to extrapolate the atmospheric temperature ranges. We find that the post-CCSD(T) contribution is very large for the barrier height with the value of -0.85 kcal mol-1 for the H-abstraction reaction, while the previous investigations were done up to the CCSD(T) level. Moreover, we also find that recrossing effects, tunneling, torsional anharmonicity, and anharmonicity are important for obtaining quantitative kinetics in the OH + HNCO reaction.

Mechanistic insight into the N2O + O(1D,3P) reaction: role of post-CCSD(T) corrections and non-adiabatic effects

In the present work, we have studied the N2O + O(1D,3P) reaction using high level quantum chemical calculations along with non-adiabatic kinetics. For quantum chemical calculations, we used the post-CCSD(T) method, which includes corrections from full triple excitations and partial quadratic excitations at the coupled-cluster level. For both the paths (N2 + O2 and 2NO), we have computed the rate constants over a wide range of temperatures (100-500 K for singlet paths and 700-4000 K for triplet paths). To assess the accuracy of our computations, we have compared our results with various experimentally measured quantities (absolute rate constant, branching fraction, and crossover temperature) and found a good match with all of them. We recommend the Arrhenius expressions for singlet paths, which turn out to be 4.46 × 10-11 exp(0.022/RT) cm3 molecule-1 s-1 and 7.12 × 10-11 exp(0.024/RT) cm3 molecule-1 s-1 for N2 + O2 and NO paths, respectively. For triplet paths, our recommended Arrhenius expressions are 5.15 × 10-12 exp(-15.35/RT) cm3 molecule-1 s-1 and 1.59 × 10-10 exp(-27.76/RT) cm3 molecule-1 s-1 for N2 + O2 and NO paths, respectively.

Accurate fundamental invariant-neural network representation of ab initio potential energy surfaces

Highly accurate potential energy surfaces are critically important for chemical reaction dynamics. The large number of degrees of freedom and the intricate symmetry adaption pose a big challenge to accurately representing potential energy surfaces (PESs) for polyatomic reactions. Recently, our group has made substantial progress in this direction by developing the fundamental invariant-neural network (FI-NN) approach. Here, we review these advances, demonstrating that the FI-NN approach can represent highly accurate, global, full-dimensional PESs for reactive systems with even more than 10 atoms. These multi-channel reactions typically involve many intermediates, transition states, and products. The complexity and ruggedness of this potential energy landscape present even greater challenges for full-dimensional PES representation. These PESs exhibit a high level of complexity, molecular size, and accuracy of fit. Dynamics simulations based on these PESs have unveiled intriguing and novel reaction mechanisms, providing deep insights into the intricate dynamics involved in combustion, atmospheric, and organic chemistry.

Theoretical Kinetics studies of isoprene peroxy radical chemistry: The fate of Z-δ-(4-OH, 1-OO)-ISOPOO radical

The OH radical recycling mechanism in isoprene oxidation is one of the most exciting topics in atmospheric chemistry, and the corresponding studies expand our understanding of oxidation mechanisms of volatile organic compounds in the troposphere and provide reliable evidence to improve and develop conventional atmospheric models. In this work, we performed a detailed theoretical kinetics study on the Z-δ-(4-OH, 1-OO)-ISOPOO radical chemistry, which is proposed as the heart of OH recycling in isoprene oxidation. With the full consideration of its accumulation and consumption channels, we studied and discussed the fate of Z-δ-(4-OH, 1-OO)-ISOPOO radical by solving the energy-resolved master equation over a broad range of conditions, including not only room temperatures but also high temperatures of a forest fire or low temperatures and pressures of the upper troposphere. We found non-negligible pressure dependence of its fate at combustion temperatures (up to two orders of magnitude) and demonstrated the significance of both the multi-structural torsional anharmonicity and tunneling for accurately calculating kinetics of the studied system. More interestingly, the tunneling effect on the phenomenological rate constants of the H-shift reaction channel is also found to be pressure-dependent due to the competition with the O2 loss reaction. In addition, our time evolution calculations revealed a two-stage behavior of critical species in this reaction system and estimated the shortest half-lives for the Z-δ-(4-OH, 1-OO)-ISOPOO radical at various temperatures, pressures and altitudes. This detailed kinetics study of Z-δ-(4-OH, 1-OO)-ISOPOO radical chemistry offers a typical example to deeply understand the core mechanism of OH recycling pathways in isoprene oxidation, and provides valuable insights for promoting the development of relevant atmospheric models.

Significant influence of water molecules on the SO3 + HCl reaction in the gas phase and at the air-water interface

The products resulting from the reactions between atmospheric acids and SO3 have a catalytic effect on the formation of new particles in aerosols. However, the SO3 + HCl reaction in the gas-phase and at the air-water interface has not been considered. Herein, this reaction was explored exhaustively by using high-level quantum chemical calculations and Born Oppenheimer molecular dynamics (BOMD) simulations. The quantum calculations show that the gas-phase reaction of SO3 + HCl is highly unlikely to occur under atmospheric conditions with a high energy barrier of 22.6 kcal mol-1. H2O and (H2O)2 play obvious catalytic roles in reducing the energy barrier of the SO3 + HCl reaction by over 18.2 kcal mol-1. The atmospheric lifetimes of SO3 show that the (H2O)2-assisted reaction dominates over the H2O-assisted reaction within the altitude range of 0-5 km, whereas the H2O-assisted reaction is more favorable within an altitude range of 10-50 km. BOMD simulations show that H2O-induced formation of the ClSO3-⋯H3O+ ion pair and HCl-assisted formation of the HSO4-⋯H3O+ ion pair were identified at the air-water interface. These routes followed a stepwise reaction mechanism and proceeded at a picosecond time scale. Interestingly, the formed ClSO3H in the gas phase has a tendency to aggregate with sulfuric acids, ammonias, and water molecules to form stable clusters within 40 ns simulation time, while the interfacial ClSO3- and H3O+ can attract H2SO4, NH3, and HNO3 for particle formation from the gas phase to the water surface. Thus, this work will not only help in understanding the SO3 + HCl reaction driven by water molecules in the gas-phase and at the air-water interface, but it will also provide some potential routes of aerosol formation from the reaction between SO3 and inorganic acids.

Atmospheric Chemistry of NH2SO3H in Polluted Areas: An Unexpected Isomerization of NH2SO3H in Acid-Polluted Regions

NH2SO3H is an effective nucleation agent for the formation of atmospheric aerosols and cloud particles. So, the ammonolysis of SO3 to form NH2SO3H without and with neutral (H2O) and basic (NH3) trace gases has been extensively investigated. However, the acidic trace gas X (X = H2SO4 and CH3SO3H)-assisted ammonolysis of SO3 is still up for debate. In this work, a comprehensive theoretical investigation of X-assisted ammonolysis of SO3 and its reverse reaction (the isomerization of NH2SO3H to form SO3-···NH3+) was carried out in the gas phase and at the air-water interface. The gas-phase results show that X-assisted isomerization of NH2SO3H to form SO3-···NH3+ is more energetically and kinetically favorable than its reverse reaction and the isomerization of NH2SO3H in the presence of H2O and NH3. Such unexpected findings revealed that gas-phase NH2SO3H is highly reactive in the presence of acidic trace gas in contrast to the high stability of NH2SO3H in neutral and basic conditions. At the air-water interface, the X-assisted isomerization reaction of NH2SO3H involves multiple water molecules. The loop structure of the reaction center (X···NH2SO3H···3H2O) promotes the transfer of protons in the water molecules to form the SO3-···NH3+ ion pair, which can then interact with several interfacial water molecules to form ammonium bisulfate. These interfacial reaction channels follow a stepwise mechanism and proceed at the picosecond time-scale. The findings of this study will contribute to a better understanding of the atmospheric behavior of NH2SO3H in polluted acidic trace gases.

Sources, sinks, and chemistry of Stabilized Criegee Intermediates in the Indo-Gangetic Plain

Night-time oxidation significantly affects the atmospheric concentration of primary and secondary air pollutants but is poorly constrained over South Asia. Here, using a comprehensively measured and unprecedented set of precursors and sinks of Stabilized Criegee Intermediates (SCI), in the summertime air of the Indo-Gangetic Plain (IGP), we investigate the chemistry, and abundance in detail. This study reports the first summertime levels from the IGP of ethene, propene, 1-butene, cis-2-butene, trans-2-butene, 1-pentene, cis-2-pentene, trans-2-pentene, and 1-hexene and their possible roles in SCI chemistry. Ethene, propene, and 1-butene were the highest ambient alkenes in both the summer and winter seasons. Applying chemical steady-state to the measured precursors, the average calculated SCI concentrations were 4.4 (±3.6) × 103 molecules cm-3, with Z-CH3CHOO (55 %) as the major SCI. Z-RCHOO (35 %) and α-pinene derived PINOO (34 %) were identified as the largest contributors to SCI with a 7.8 × 105 molecules cm-3 s-1 production rate. The peak SCI occurred during the evenings. For all SCI species, the loss was dominated (>50 %) by unimolecular decomposition or reactions with water vapor or water vapor dimer. Pollution events influenced by crop burning resulted in significantly elevated SCI production (2.1 times higher relative to non-polluted periods) reaching as high as (7.4 ± 2.5) × 105 molecules cm-3 s-1. Among individual SCI species, Z-CH3CHOO was highest in all the plume events measured accounting for at least ~41 %. Among alkenes, trans-2-butene was the highest contributor to P(SCI) in plume events ranging from 22 to 32 %. SCIs dominated the night-time oxidation of sulfur dioxide with rates as high as 1.4 (±1.1) × 104 molecules cm-3 s-1 at midnight, suggesting that this oxidation pathway could be a significant source of fine mode sulfate aerosols over the Indo-Gangetic Plain, especially during summertime biomass burning pollution episodes.

Kinetics of Sulfur Trioxide Reaction with Water Vapor to Form Atmospheric Sulfuric Acid

Although experimental methods can be used to obtain the quantitative kinetics of atmospheric reactions, experimental data are often limited to a narrow temperature range. The reaction of SO3 with water vapor is important for elucidating the formation of sulfuric acid in the atmosphere; however, the kinetics is uncertain at low temperatures. Here, we calculate rate constants for reactions of sulfur trioxide with two water molecules. We consider two mechanisms: the SO3···H2O + H2O reaction and the SO3 + (H2O)2 reaction. We find that beyond-CCSD(T) contributions to the barrier heights are very large, and multidimensional tunneling, unusually large anharmonicity of high-frequency modes, and torsional anharmonicity are important for obtaining quantitative kinetics. We find that at lower temperatures, the formation of the termolecular precursor complexes, which is often neglected, is rate-limiting compared to passage through the tight transition states. Our calculations show that the SO3···H2O + H2O mechanism is more important than the SO3 + (H2O)2 mechanism at 5-50 km altitudes. We find that the rate ratio between SO3···H2O + H2O and SO3 + (H2O)2 is greater than 20 at altitudes between 10 and 35 km, where the concentration of SO3 is very high.

The simplest Criegee intermediate CH2OO reaction with dimethylamine and trimethylamine: kinetics and atmospheric implications

We have used the OH laser-induced fluorescence (LIF) method to measure the kinetics of the simplest Criegee intermediate (CH2OO) reacting with two abundant amines in the atmosphere: dimethylamine ((CH3)2NH) and trimethylamine ((CH3)3N). Our experiments were conducted under pseudo-first-order approximation conditions. The rate coefficients we report are (2.15 ± 0.28) × 10-11 cm3 molecule-1 s-1 for (CH3)2NH at 298 K and 10 Torr, and (1.56 ± 0.23) × 10-12 cm3 molecule-1 s-1 for (CH3)3N at 298 K and 25 Torr with Ar as the bath gas. Both reactions exhibit a negative temperature dependence. The activation energy and pre-exponential factors derived from the Arrhenius equation were (-2.03 ± 0.26) kcal mol-1 and (6.89 ± 0.90) × 10-13 cm3 molecule-1 s-1 for (CH3)2NH, and (-1.60 ± 0.24) kcal mol-1 and (1.06 ± 0.16) × 10-13 cm3 molecule-1 s-1 for (CH3)3N. We propose that the electronegativity of the atom in the co-reactant attached to the C atom of CH2OO, in addition to the dissociation energy of the fragile covalent bonds with H atoms (H-X bond), plays an important role in the 1,2-insertion reactions. Under certain circumstances, the title reactions can contribute to the sink of amines and Criegee intermediates and to the formation of secondary organic aerosol (SOA).

Unimolecular Reactions of E-Glycolaldehyde Oxide and Its Reactions with One and Two Water Molecules

The kinetics of Criegee intermediates are important for atmospheric modeling. However, the quantitative kinetics of Criegee intermediates are still very limited, especially for those with hydroxy groups. Here, we calculate rate constants for the unimolecular reaction of E-glycolaldehyde oxide [E-hydroxyethanal oxide, E-(CH₂OH)CHOO], for its reactions with H₂O and (H₂O)₂, and for the reaction of the E-(CH₂OH)CHOO…H₂O complex with H₂O. For the highest level of electronic structure, we use W3X-L//CCSD(T)-F12a/cc-pVDZ-F12 for the unimolecular reaction and the reaction with water and W3X-L//DF-CCSD(T)-F12b/jun-cc-pVDZ for the reaction with 2 water molecules. For the dynamics, we use a dual-level strategy that combines conventional transition state theory with the highest level of electronic structure and multistructural canonical variational transition state theory with small-curvature tunneling with a validated density functional for the electronic structure. This dynamical treatment includes high-frequency anharmonicity, torsional anharmonicity, recrossing effects, and tunneling. We find that the unimolecular reaction of E-(CH₂OH)CHOO depends on both temperature and pressure. The calculated results show that E-(CH₂OH)CHOO…H₂O + H₂O is the dominant entrance channel, while previous investigations only considered Criegee intermediates + (H₂O)₂. In addition, we find that the atmospheric lifetime of E-(CH₂OH)CHOO with respect to 2 water molecules is particularly short with a value of 1.71 × 10^-6 s at 0 km, which is about 2 orders of magnitude shorter than those usually assumed for Criegee intermediate reactions with water dimer. We also find that the OH group in E-(CH₂OH)CHOO enhances its reactivity.

The atmospheric relevance of primary alcohols and imidogen reactions

Organic alcohols as very volatile compounds play a crucial role in the air quality of the atmosphere. So, the removal processes of such compounds are an important atmospheric challenge. The main goal of this research is to discover the atmospheric relevance of degradation paths of linear alcohols by imidogen with the aid of simulation by quantum mechanical (QM) methods. To this end, we combine broad mechanistic and kinetic results to get more accurate information and to have a deeper insight into the behavior of the designed reactions. Thus, the main and necessary reaction pathways are explored by well-behaved QM methods for complete elucidation of the studying gaseous reactions. Moreover, the potential energy surfaces as  a main factor are computed for easier judging of the most probable pathways in the simulated reactions. Our attempt to find the occurrence of the considered reactions in the atmospheric conditions is completed by precisely evaluating the rate constants of all elementary reactions. All of the computed bimolecular rate constants have a positive dependency on both temperature and pressure. The kinetic results show that H-abstraction from the α carbon is dominant relative to the other sites. Finally, by the results of this study, we conclude that at moderate temperatures and pressures primary alcohols can degrade with imidogen, so they can get atmospheric relevance.

Master Equation Studies of the Unimolecular Decay of Thermalized Methacrolein Oxide: The Impact of Atmospheric Conditions

Master equation simulations of the unimolecular reaction dynamics of the Criegee intermediate methacrolein oxide (MACR oxide) have been performed under a variety of temperature and pressure conditions. These simulations provide insight into how the unimolecular kinetics vary across temperatures spanning the range 288-320 K. This work has incorporated a new potential energy surface and includes the anti-to-syn and cis-to-trans conformational dynamics of MACR oxide, as well as the unimolecular reactions to form dioxirane and dioxole species. The competition between the unimolecular reactivity of MACR oxide and previously documented bimolecular reactivity of MACR oxide with water vapor is explored, focusing on how this competition is affected by changes in atmospheric conditions. The impact on the role of MACR oxide as an atmospheric oxidant of SO₂ is noted.

Full-dimensional neural network potential energy surface and dynamics of the CH2OO + H2O reaction

An accurate global full-dimensional machine learning-based potential energy surface (PES) of the simplest Criegee intermediate (CH₂OO) reaction with water monomer was developed based on the high level of extensive CCSD(T)-F12a/aug-cc-pVTZ calculations. This analytical global PES not only covers the regions of reactants to hydroxymethyl hydroperoxide (HMHP) intermediates, but also different end product channels, which facilities both the reliable and efficient kinetics and dynamics calculations. The rate coefficients calculated by the transition state theory with the interface to the full-dimensional PES agree well with the experimental results, indicating the accuracy of the current PES. Extensive quasi-classical trajectory (QCT) calculations were performed both from the bimolecular reaction CH₂OO + H₂O and from HMHP intermediate on the new PES. The product branching ratios of hydroxymethoxy radical (HOCH₂O, HMO) + OH radical, formaldehyde (CH₂O) + H₂O₂ and formic acid (HCOOH) + H₂O were calculated. The reaction yields dominantly HMO + OH, because of the barrierless pathway from HMHP to this channel. The computed dynamical results for this product channel show the total available energy was deposited into the internal rovibrational excitation of HMO, and the energy release in OH and translational energy is limited. The large amount of OH radical found in the current study implies that the CH₂OO + H₂O reaction can provide crucially OH yield in Earth's atmosphere.

Water Plays Multifunctional Roles in the Intervening Formation of Secondary Organic Aerosols in Ozonolysis of Limonene: A Valence Photoelectron Spectroscopy and Density Functional Theory Study

Although water may affect aqueous aerosol chemistry, how it intervenes in the formation of secondary organic aerosols (SOAs) at the molecular level remains elusive. Ozonolysis of limonene is one of the most important sources of indoor SOAs. Here, we report the valence electronic properties of limonene aerosols and SOAs derived from limonene ozonolysis (Lim-SOAs) via aerosol vacuum ultraviolet photoelectron spectroscopy, with a focus on the effects of water on Lim-SOAs. The first vertical ionization energy of limonene aerosols is measured to be 8.79 ± 0.07 eV. While water significantly increases the total photoelectron yield of Lim-SOAs, three photoelectron features attributable to Lim-SOAs each exhibit distinct dependence on the fraction of water in aerosols, implying that different formation pathways and molecular origins are involved in the formation of Lim-SOAs. Combined with density functional theory calculation and mass spectrometry measurements, this study reveals that water, particularly the water dimer, enhances the formation of Lim-SOAs by altering the ozonolysis energetics and pathways by intervening in its Criegee chemistry, acting as both a catalyst and a reactant. The atmospheric implication is discussed.

Water Charge Transfer Accelerates Criegee Intermediate Reaction with H2O- Radical Anion at the Aqueous Interface

Criegee intermediates (CIs) are important carbonyl oxides that may react with atmospheric trace chemicals and impact the global climate. The CI reaction with water has been widely studied and is a main channel for trapping CIs in the troposphere. Previous experimental and computational reports have largely focused on reaction kinetic processes in various CI-water reactions. The molecular-level origin of CI's interfacial reactivity at the water microdroplet surface (e.g., as found in aerosols and clouds) is unclear. In this study, by employing the quantum mechanical/molecular mechanical (QM/MM) Born-Oppenheimer molecular dynamics with the local second-order Møller-Plesset perturbation theory, our computational results reveal a substantial water charge transfer up to ∼20% per water, which creates the surface H₂O⁺/H₂O^- radical pairs to enhance the CH₂OO and anti-CH₃CHOO reactivity with water: the resulting strong CI-H₂O^- electrostatic attraction at the microdroplet surface facilitates the nucleophilic attack to the CI carbonyl by water, which may counteract the apolar hindrance of the substituent to accelerate the CI-water reaction. Our statistical analysis of the molecular dynamics trajectories further resolves a relatively long-lived bound CI(H₂O^-) intermediate state at the air/water interface, which has not been observed in gaseous CI reactions. This work provides insights into what may alter the oxidizing power of the troposphere by the next larger CIs than simple CH₂OO and implicates a new perspective on the role of interfacial water charge transfer in accelerating molecular reactions at aqueous interfaces.

Kinetics of the Simplest Criegee Intermediate CH2OO Reaction with tert-Butylamine

The kinetics of the simplest Criegee intermediate (CH₂OO) reaction with tert-butylamine ((CH₃)₃CNH₂) was studied under pseudo-first-order conditions with the OH laser-induced fluorescence (LIF) method at the temperature range of 283-318 K and the pressure range of 5-75 Torr. Our pressure-dependent measurement showed that at 5 Torr─the lowest pressure measured in the current experiment─this reaction was under the high-pressure limit condition. At 298 K, the reaction rate coefficient was measured to be (4.95 ± 0.64) × 10^-12 cm³ molecule^-1 s^-1. The title reaction was observed to be negative temperature-dependent; the activation energy of (-2.82 ± 0.37) kcal mol^-1 and the pre-exponential factor of (4.21 ± 0.55) × 10^-14 cm³ molecule^-1 s^-1 were derived from the Arrhenius equation. The rate coefficient of the title reaction is slightly larger than (4.3 ± 0.5) × 10^-12 cm³ molecule^-1 s^-1 of the CH₂OO reaction with methylamine; the electron inductive effect and the steric hindrance effect might play a role in contributing to such difference.

Interfacial Extraction to Trap and Characterize the Criegee Intermediates from Phospholipid Ozonolysis

Criegee intermediates (CIs) play a significant role in cell membrane peroxidation, but their identification remains elusive at the molecular level. Herein, we combined interfacial extraction and sonic spray ionization mass spectrometry to study the oxidation reaction of 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-(1'-rac-glycerol) (POPG) mediated by ozone (O₃) at/near the surface of a hung water droplet. On-line interfacial extraction and ionization provided a snapshot of the short-lived CIs. Experiments in which the content of water was varied provided evidence for the formation of CIs, which has not been previously observed. Capture experiments using 5,5-dimethyl-pyrroline N-oxide (DMPO) indicated that CIs could be selectively characterized, and the extracted ion current (EICs) of CIs vs DMPO-CI adducts further confirmed the successful observation of CIs. Theoretical calculation suggested that surface ozonolysis of POPG was mainly mediated by anti-CI. These results open a new route for aqueous surface reactive species identification, and benefit toward the understanding of disease development associated with cell oxidative stress mediated by CIs.

Semiclassical Dynamics on Machine-Learned Coupled Multireference Potential Energy Surfaces: Application to the Photodissociation of the Simplest Criegee Intermediate

Determination of high-dimensional potential energy surfaces (PESs) and nonadiabatic couplings have always been quite challenging. To this end, machine learning (ML) models, trained with a finite set of ab initio data, allow accurate prediction of such properties. To express the PESs in terms of atomic contributions is the cornerstone of any ML based technique because it can be easily scaled to large systems. In this work, we have constructed high fidelity PESs and nonadiabatic coupling terms at the CASSCF level of ab initio data using a machine learning technique, namely, kernel-ridge regression. Additional MRCI-level calculations were carried out to assess the quality of the PESs. We use these machine-learned PESs and nonadiabatic couplings to simulate excited-state molecular dynamics based on Tully's fewest-switches surface hopping method (FSSH). FSSH is a semiclassical method in which nuclei move on the PESs due to the electrons according to the laws of classical mechanics. Nonadiabatic effects are taken into account in terms of transitions between PESs. We apply this scheme to study the O-O photodissociation of the simplest Criegee intermediate (CH₂OO). The FSSH trajectories were initiated on the lowest optically bright singlet excited state (S₂) and propagated along the three most important internal coordinates, namely, O-O and C-O bond distances and the COO bond angle. Some of the trajectories end up on energetically lower PESs as a result of radiationless transfer through conical intersections. All of the trajectories lead to the dissociation of the O-O bond due to the dissociative nature of the excited PESs through one of the two dissociative channels. The simulation reveals that there is about 88.4% probability of dissociation through the lower channel leading to the H₂CO (X¹A₁) and O (¹D) products, whereas there is only 11.6% probability of dissociation through the upper channel leading to H₂CO (a³A″) and O (³P) products.

Accurate determination of reaction energetics and kinetics of the HO2˙ + O3 → OH˙ + 2O2 reaction

In the present work, we have studied the HO₂˙ + O₃ → HO˙ + 2O₂ reaction using chemical kinetics and quantum chemical calculations. We have employed the post-CCSD(T) method to estimate the barrier height and reaction energy for the title reaction. In the post-CCSD(T) method, we have included zero point energy corrections, contributions from full triple excitations and partial quadratic excitations at the coupled-cluster level, and core corrections. We have also computed the reaction rate in the temperature range of 197-450 K and found good agreement with all the available experimental results. In addition, we have also fitted the computed rate constants with the Arrhenius expression and obtained an activation energy of 1.0 ± 0.1 kcal mol^-1, almost identical to the value recommended by IUPAC and JPL.

Predicting the environmental fates of emerging contaminants: Synergistic effects in ozone reactions of nitrogen-containing alkenes

While nitro and amino alkenes are common in pharmaceuticals, pesticides, and munitions, their environmental fates are not well known. Ozone is a ubiquitous atmospheric oxidant for alkenes, but the synergistic effects of nitrogen-containing groups on the reactions have not been measured. The kinetics and products of ozonolysis of a series of model compounds with different combinations of these functional groups have been measured in the condensed phase using stopped-flow and mass spectrometry methods. Rate constants span about six orders of magnitude with activation energies ranging from 4.3 to 28.2 kJ mol^-1. Vinyl nitro groups substantially decrease the reactivity, while amino groups have the opposite effect. The site of the initial ozone attack is highly structure dependent, consistent with local ionization energy calculations. The reaction of the neonicotinoid pesticide nitenpyram, which forms toxic N-nitroso compounds, was consistent with model compounds, confirming the utility of model compounds for assessing environmental fates of these emerging contaminants.

A spur to molecular geometry optimization: Gradient-enhanced universal kriging with on-the-fly adaptive ab initio prior mean functions in curvilinear coordinates

We present a molecular geometry optimization algorithm based on the gradient-enhanced universal kriging (GEUK) formalism with ab initio prior mean functions, which incorporates prior physical knowledge to surrogate-based optimization. In this formalism, we have demonstrated the advantage of allowing the prior mean functions to be adaptive during geometry optimization over a pre-fixed choice of prior functions. Our implementation is general and flexible in two senses. First, the optimizations on the surrogate surface can be in both Cartesian coordinates and curvilinear coordinates. We explore four representative curvilinear coordinates in this work, including the redundant Coulombic coordinates, the redundant internal coordinates, the non-redundant delocalized internal coordinates, and the non-redundant hybrid delocalized internal Z-matrix coordinates. We show that our GEUK optimizer accelerates geometry optimization as compared to conventional non-surrogate-based optimizers in internal coordinates. We further showcase the power of the GEUK with on-the-fly adaptive priors for efficient optimizations of challenging molecules (Criegee intermediates) with a high-accuracy electronic structure method (the coupled-cluster method). Second, we present the usage of internal coordinates under the complete curvilinear scheme. A complete curvilinear scheme performs both surrogate potential-energy surface (PES) fitting and structure optimization entirely in the curvilinear coordinates. Our benchmark indicates that the complete curvilinear scheme significantly reduces the cost of structure minimization on the surrogate compared to the incomplete curvilinear scheme, which fits the surrogate PES in curvilinear coordinates partially and optimizes a structure in Cartesian coordinates through curvilinear coordinates via the chain rule.

Organic Peroxides in Aerosol: Key Reactive Intermediates for Multiphase Processes in the Atmosphere

Organic peroxides (POs) are organic molecules with one or more peroxide (-O-O-) functional groups. POs are commonly regarded as chemically labile termination products from gas-phase radical chemistry and therefore serve as temporary reservoirs for oxidative radicals (HOx and ROx) in the atmosphere. Owing to their ubiquity, active gas-particle partitioning behavior, and reactivity, POs are key reactive intermediates in atmospheric multiphase processes determining the life cycle (formation, growth, and aging), climate, and health impacts of aerosol. However, there remain substantial gaps in the origin, molecular diversity, and fate of POs due to their complex nature and dynamic behavior. Here, we summarize the current understanding on atmospheric POs, with a focus on their identification and quantification, state-of-the-art analytical developments, molecular-level formation mechanisms, multiphase chemical transformation pathways, as well as environmental and health impacts. We find that interactions with SO2 and transition metal ions are generally the fast PO transformation pathways in atmospheric liquid water, with lifetimes estimated to be minutes to hours, while hydrolysis is particularly important for α-substituted hydroperoxides. Meanwhile, photolysis and thermolysis are likely minor sinks for POs. These multiphase PO transformation pathways are distinctly different from their gas-phase fates, such as photolysis and reaction with OH radicals, which highlights the need to understand the multiphase partitioning of POs. By summarizing the current advances and remaining challenges for the investigation of POs, we propose future research priorities regarding their origin, fate, and impacts in the atmosphere.

Electronic Absorption Spectroscopy and Photochemistry of Criegee Intermediates

Interest in Criegee intermediates (CIs), often termed carbonyl oxides, and their role in tropospheric chemistry has grown massively since the demonstration of laboratory-based routes to their formation and characterization in the gas phase. This article reviews current knowledge regarding the electronic spectroscopy of atmospherically relevant CIs like CH₂ OO, CH₃ CHOO, (CH₃ )₂ COO and larger CIs like methyl vinyl ketone oxide and methacrolein oxide that are formed in the ozonolysis of isoprene, and of selected conjugated carbene-derived CIs of interest in the synthetic chemistry community. Of the aforementioned atmospherically relevant CIs, all except CH₂ OO and (CH₃ )₂ COO exist in different conformers which, under tropospheric conditions, can display strikingly different thermal loss rates via unimolecular and bimolecular processes. Calculated photolysis rates based on their absorption properties suggest that solar photolysis will rarely be a significant contributor to the total loss rate for any CI under tropospheric conditions. Nonetheless, there is ever-growing interest in the absorption cross sections and primary photochemistry of CIs following excitation to the strongly absorbing ¹ ππ* state, and how this varies with CI, with conformer and with excitation wavelength. The later part of this review surveys the photochemical data reported to date, including a range of studies that demonstrate prompt photo-induced fission of the terminal O-O bond, and speculates about possible alternate decay processes that could occur following non-adiabatic coupling to, and dissociation from, highly internally excited levels of the electronic ground state of a CI.

Barrier heights, reaction energies and bond dissociation energies for RH + HO2 reactions with coupled-cluster theory, density functional theory and diffusion quantum Monte Carlo methods

Hydrogen abstraction reactions by the HO₂ radical from hydrocarbon molecules are an important class of reactions in the autoignition of hydrocarbon fuels. Performance of DLPNO-CC and DFT methods using three hybrids and four double hybrids as well as FN-DMC with the single-Slater-Jastrow trial wavefunction on barrier heights and reaction energies of RH + HO₂ reactions as well as bond dissociation energies of the involved X-H molecules is evaluated by comparison with the highly accurate CCSD(T)-F12b/CBS results in this study. Our results show that the DLPNO-CCSD(T)-F12 method can achieve highly accurate barrier heights, reaction energies and X-H bond energies for RH + HO₂ reactions at a relatively low computational cost, and it is applicable to the H-abstraction reactions of larger molecules. Among all DFAs, MN15 and the employed double hybrids can achieve accurate barrier heights and reaction energies with MADs of less than or around 2 kJ mol^-1, but their error on X-H bond energies is more pronounced. Only DSD-BLYP and DSD-PBEB95 can provide X-H bond energies with MADs less than 4 kJ mol^-1. Considering dispersion correction in DFT calculations does not improve these barrier heights and reaction energies. The error of FN-DMC on barrier heights and reaction energies is slightly larger than that of MN15 and those of double hybrids, but it can achieve results within chemical accuracy for these reactions and the X-H bond energies.

Peroxy Radical Autoxidation and Sequential Oxidation in Organic Nitrate Formation during Limonene Nighttime Oxidation

Limonene is an abundant monoterpene released into the atmosphere via biogenic emissions and biomass burning. However, the atmospheric oxidation and secondary organic aerosol (SOA) formation mechanisms of limonene, especially during nighttime, remain largely understudied. In this work, limonene was oxidized synergistically by ozone (O₃) and nitrate radicals (NO₃) in a flow tube reactor and a continuous flow stirred tank reactor. Upon oxidation, many highly oxidized organic nitrates and nitrooxy peroxy radicals (RO₂) were observed in the gas phase within 1 min. Combining quantum chemical calculations with kinetic simulations, we found that the primary nitrooxy RO₂ (C₁₀H₁₆NO₅) through NO₃ addition at the more substituted endocyclic double bond and at the exocyclic double bond (previously considered as minor pathways) can undergo autoxidation with rate constants of around 0.02 and 20 s^-1 at 298 K, respectively. These pathways could explain a major portion of the observed highly oxidized organic nitrates. In the SOA, highly oxidized mono- and dinitrates (e.g., C₁₀H₁₇NO_7-8 and C₁₀H₁₆,₁₈N₂O_8-10) make up a significant contribution, highlighting nitrooxy RO₂ autoxidation and sequential NO₃ oxidation of limonene. The same organic nitrates are also observed in ambient aerosol during biomass burning and nighttime in the southeastern United States. Therefore, the present work provides new insights into the nighttime oxidation of limonene and SOA formation in the atmosphere.

Quantitative Kinetics of HO2 Reactions with Aldehydes in the Atmosphere: High-Order Dynamic Correlation, Anharmonicity, and Falloff Effects Are All Important

Kinetics provides the fundamental parameters for elucidating sources and sinks of key atmospheric species and for atmospheric modeling more generally. Obtaining quantitative kinetics in the laboratory for the full range of atmospheric temperatures and pressures is quite difficult. Here, we use computational chemistry to obtain quantitative rate constants for the reactions of HO₂ with HCHO, CH₃CHO, and CF₃CHO. First, we calculate the high-pressure-limit rate constants by using a dual-level strategy that combines conventional transition state theory using a high level of electronic structure wave function theory with canonical variational transition state theory including small-curvature tunneling using density functional theory. The wave-function level is beyond-CCSD(T) for HCHO and CCSD(T)-F12a (Level-A) for XCHO (X = CH₃, CF₃), and the density functional (Level-B) is specifically validated for these reactions. Then, we calculate the pressure-dependent rate constants by using system-specific quantum RRK theory (SS-QRRK) and also by an energy-grained master equation. The two treatments of the pressure dependence agree well. We find that the Level-A//Level-B method gives good agreement with CCSDTQ(P)/CBS. We also find that anharmonicity is an important factor that increases the rate constants of all three reactions. We find that the HO₂ + HCHO reaction has a significant dependence on pressure, but the HO₂ + CF₃CHO reaction is almost independent of pressure. Our findings show that the HO₂ + HCHO reaction makes important contribution to the sink for HCHO, and the HO₂ + CF₃CHO reaction is the dominant sink for CF₃CHO in the atmosphere.

Antioxidant activity of caffeic acid: thermodynamic and kinetic aspects on the oxidative degradation pathway

Caffeic acid is a phenolic secondary metabolite from plants, which is known for its antioxidant properties. The effective mitigation of methanol-induced oxidative stress by caffeic acid depends on the direct radical scavenging as well as the formation of new metabolites via oxidative degradation. Herein, thermodynamic and kinetic aspects of the oxidative degradation pathway of caffeic acid in the presence of radical CH₃O^• and its isomer, ^•CH₂OH are discussed for the first time, employing density functional theory (DFT). The direct radical scavenging activity of caffeic acid against these radicals is verified via hydrogen atom transfer (HAT) and radical adduct formation (RAF) mechanisms. HAT is predicted to be more feasible than RAF mechanism as per the computed data. Additionally, energetic details of the proposed oxidative degradation pathway of radical adduct intermediates toward the formation of a cyclic metabolite is analyzed. Kinetic studies indicated a significant tunneling contribution to the H abstraction pathways having high activation barriers. Further, our results imply that the newly formed metabolites exhibit comparable antioxidant activity with that of caffeic acid.

Theoretical study on the formation of Criegee intermediates from ozonolysis of pentenal: An example of trans-2-pentenal

In this study, we investigated the reaction mechanism and kinetics of ozone with trans-2-pentenal using density functional theory (DFT) and conventional transition state theory (CTST). At 298 K and 1 atm, the gas-phase reaction mechanisms and kinetic parameters were calculated at the level of CCSD(T)/6-311+G(d,p)//M06-2X/6-311+G(d,p). Both CC and CO bond cycloaddition as well as hydrogen abstraction were found. The calculations indicated that the main reaction path is 1,3-dipole cycloaddition reactions of ozone with CC bond with the relatively lower syn-energy-barrier of 3.35 kcal mol^-1 to form primary ozonide which decomposed to produce a carbonyl oxide called a Criegee intermediate (CI) and an aldehyde. The subsequent reactions of CIs were analysed in detail. It is found that the reaction pathways of the novelty CIs containing an aldehyde group are extremely similar with general CIs when they react with NO, NO₂, SO₂, H₂O, CH₂O and O₂. The condensed Fukui function were calculated to identify the active site of the chosen molecules. At 298 K and 1 atm, the reaction rate coefficient was 9.13 × 10^-18 cm³ molecule^-1 s^-1 with atmospheric lifetime of 1.3 days. The calculated rate constant is in general agreement with the available experimental data. The branching ratios indicated that syn-addition pathways are prior to anti-addition. The atmospheric ratios for CIs formation and the bimolecular reaction rate constants for the Criegee intermediates with the variety of partners were calculated. Our theoretical results are of importance in atmospheric chemistry of unsaturated aldehyde oxidation by ozone.

Mechanistic Insights into Criegee Intermediate-Hydroperoxyl Radical Chemistry

The reaction between a Criegee intermediate and the hydroperoxyl radical (HO₂) is believed to play a role in the formation of new particles in the troposphere. Although the reaction has been previously studied in the gas phase, there are several knowledge gaps that still need to be filled. We simulated the reaction of anti-CH₃CHOO with HO₂ and HO₂-H₂O radical complexes in the gas phase at 0 K, which exhibited a low-barrier profile for water-containing systems and a barrierless profile for water-free systems. Moreover, the reaction was found to follow a proton-transfer mechanism, which challenges previous assumptions that the gas-phase reaction involves a hydrogen atom transfer. The HO₂ radical was found to mediate the Criegee hydration reaction in the gas phase. Metadynamics simulations further confirmed that the expected radical adduct formation between anti-CH₃CHOO and the HO₂ radical, as well as the HO₂- and H₂O-mediated reactions in the gas phase, followed a concerted mechanism. By combining constrained ab initio molecular dynamics simulations with thermodynamic integration, we quantitively evaluated the free-energy barriers at high temperatures. The barriers obtained for all three Criegee-HO₂ reaction systems were found to be temperature-dependent. We also compared the free-energy barriers of water-free and water-containing systems; the results revealed that water could hinder the reaction between the Criegee and HO₂ radical. These results suggest that HO₂ radicals may be involved in the formation of tropospheric radical adducts, and water molecules may also play important roles in the reactions of Criegee intermediates.

Oxidation of HOSO˙ by O2 (3Σg-): a key reaction deciding the fate of HOSO˙ in the atmosphere

In the present work, we have studied the oxidation of HOSO˙ by O₂ (³Σ_g^-) employing quantum chemical and kinetic calculations. The present work reveals that HOSO˙ + O₂ (³Σ_g^-) is a barrierless reaction which proceeds through a stable hydrogen-bonded complex. The estimated atmospheric lifetime of HOSO˙ in the presence of O₂ (³Σ_g^-) is found to be several orders of magnitude less compared to the other oxidation paths of HOSO˙, suggesting that the oxidation of HOSO˙ by O₂ (³Σ_g^-) might be the most dominant oxidation path of HOSO˙ in the atmosphere.