A novel and facile decay path of Criegee intermediates by intramolecular insertion reactions via roaming transition states

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.

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.

OH-Formation following vibrationally induced reaction dynamics of H2COO

The reaction dynamics of H2COO to form HCOOH and dioxirane as first steps for OH-elimination is quantitatively investigated. Using a machine learned potential energy surface (PES) at the CASPT2/aug-cc-pVTZ level of theory vibrational excitation along the CH-normal mode νCH with energies up to 40.0 kcal mol-1 (∼5νCH) leads almost exclusively to HCOOH which further decomposes into OH + HCO. Although the barrier to form dioxirane is only 21.4 kcal mol-1 the reaction probability to form dioxirane is two orders of magnitude lower if the CH-stretch mode is excited. Following the dioxirane-formation pathway is facile, however, if the COO-bend vibration is excited together with energies equivalent to ∼2νCH or ∼3νCOO. For OH-formation in the atmosphere the pathway through HCOOH is probably most relevant because the alternative pathways (through dioxirane or formic acid) involve several intermediates that can de-excite through collisions, relax via internal vibrational relaxation (IVR), or pass through loose and vulnerable transition states (formic acid). This work demonstrates how, by selectively exciting particular vibrational modes, it is possible to dial into desired reaction channels with a high degree of specificity.

Unimolecular Kinetics of Stabilized CH3CHOO Criegee Intermediates: syn-CH3CHOO Decomposition and anti-CH3CHOO Isomerization

The kinetics of the unimolecular decomposition of the stabilized Criegee intermediate syn-CH3CHOO has been investigated at temperatures between 297 and 331 K and pressures between 12 and 300 Torr using laser flash photolysis of CH3CHI2/O2/N2 gas mixtures coupled with time-resolved broadband UV absorption spectroscopy. Fits to experimental results using the Master Equation Solver for Multi-Energy well Reactions (MESMER) indicate that the barrier height to decomposition is 67.2 ± 1.3 kJ mol-1 and that there is a strong tunneling component to the decomposition reaction under atmospheric conditions. At 298 K and 760 Torr, MESMER simulations indicate a rate coefficient of 150-81+176 s-1 when tunneling effects are included but only 5-2+3 s-1 when tunneling is not considered in the model. MESMER simulations were also performed for the unimolecular isomerization of the stabilized Criegee intermediate anti-CH3CHOO to methyldioxirane, indicating a rate coefficient of 54-21+34 s-1 at 298 K and 760 Torr, which is not impacted by tunneling effects. Expressions to describe the unimolecular kinetics of syn- and anti-CH3CHOO are provided for use in atmospheric models, and atmospheric implications are discussed.

Formation of extremely low-volatility organic compounds from styrene ozonolysis: Implication for nucleation

Styrene is a highly reactive compound with the dual nature of aromatics and olefins. This work presents evidence for formation of extremely low-volatility organic compounds (ELVOCs) from styrene ozonolysis. The molecules of ELVOCs were analyzed using a high-resolution orbitrap mass spectrometer. The results show that ELVOCs were oligomers characterized by stabilized Criegee radicals (SCIs) as chain units. The addition of C₆H₅CHOO (SCI1) or CH₂OO (SCI2) can dramatically decrease the oligomers' volatility. At low relative humidity (RH), ELVOCs are mainly formed from the reaction of RO₂ radical, C₆H₅OO·, with SCI1 and SCI2; however, ELVOCs are primarily produced by the reaction between benzoic acid and SCI1 at high RH. Ambient particles were also collected to propose the probable oligomers from styrene-SCI. Our results suggest that styrene-SCI derived ELVOCs may act as nucleating agents, potentially providing an experimental basis for nucleation events that frequently occur in urban areas.

Dynamics of Reaction CH3CHI + O2 Investigated via Infrared Emission of Products CO, CO2, and OH

The reaction CH₃CHI + O₂ has been commonly employed in laboratories to produce a methyl-substituted Criegee intermediate CH₃CHOO, but the detailed dynamics of this reaction remain unexplored. We carried out this reaction by irradiating a flowing mixture of CH₃CHI₂ (∼70 mTorr) and O₂ (∼4 and 8 Torr) at 308 or 248 nm and observed infrared emission of the products with a step-scan Fourier-transform spectrometer. Upon irradiation at 248 nm with O₂ ∼4 Torr, a Boltzmann distribution of CO (v ≤ 4, J ≤ 25) with average vibrational energy (12 ± 2) kJ mol^-1 and of OH (v = 1, J ≤ 5.5) were observed and assigned to be produced from the decomposition of CH₃C(O)OH* to form CO + CH₃OH and OH + CH₃CO, respectively. The observed broadband emission of CO₂ was simulated with two vibrational distributions of average energies (42 ± 3) and (114 ± 6) kJ mol^-1 and assigned to be produced from the decomposition of CH₃C(O)OH* and (methyl dioxirane)*, respectively. The results upon irradiation of the sample at 308 nm are similar, likely indicating a small fraction of energy partition into these products and rapid thermalization of CH₃CHI*. Compared with reaction CH₂I + O₂, the title reaction yielded products with much less internal excitation, consistent with the expectation that these observed products receive much less fraction of available energy upon fragmentation when an additional methyl moiety was present in the parent. The large-v component of CO observed in experiments of CH₂I + O₂ at 248 nm, produced from secondary reaction HCO + O₂, was absent in this work because the corresponding secondary reaction CH₃CO + O₂ in decomposition of CH₃CHOO* produces α-lactone + OH or H₂CO + CO + OH.

The favorable routes for the hydrolysis of CH2OO with (H2O)n (n = 1-4) investigated by global minimum searching combined with quantum chemical methods

The hydrolysis reaction of CH₂OO with water and water clusters is believed to be a dominant sink for the CH₂OO intermediate in the atmosphere. However, the favorable route for the hydrolysis of CH₂OO with water clusters is still unclear. Here global minimum searching using the Tsinghua Global Minimum program has been introduced to find the most stable geometry of the CH₂OO(H₂O)_n (n = 1-4) complex firstly. Then, based on these stable complexes, favorable hydrolysis of CH₂OO with (H₂O)_n (n = 1-4) has been investigated using the quantum chemical method of CCSD(T)-F12a/cc-pVDZ-F12//B3LYP/6-311+G(2d,2p) and canonical variational transition state theory with small curvature tunneling. The calculated results have revealed that, although the contribution of CH₂OO + (H₂O)₂ is the most obvious in the hydrolysis of CH₂OO with (H₂O)_n (n = 1-4), the hydrolysis of CH₂OO with (H₂O)₃ is not negligible in atmospheric gas-phase chemistry as its rate is close to the rate of the CH₂OO + H₂O reaction. The calculated results also show that, in a clean atmosphere, the CH₂OO + (H₂O)_n (n = 1-2) reaction competes well with the CH₂OO + SO₂ reaction at 298 K when the concentrations of (H₂O)_n (n = 1-2) range from 20% relative humidity (RH) to 100% RH, and SO₂ is 2.46 × 10¹¹ molecules per cm³. Meanwhile, when the RH is higher than 40%, it is a new prediction that the CH₂OO + (H₂O)₃ reaction can also compete well with the CH₂OO + SO₂ reaction at 298 K. Besides, Born-Oppenheimer molecular dynamics simulation results show that all the favorable channels of the CH₂OO + (H₂O)_n (n = 1-3) reaction cannot react on a time scale of 100 ps in the NVT simulation. However, the NVE simulation results show that the CH₂OO + (H₂O)₃ reaction can be finished well at 8.5 ps, indicating that the gas phase reaction of CH₂OO + (H₂O)₃ is not negligible in the atmosphere. Overall, the present results have provided a definitive example of how the favorable hydrolysis of important atmospheric species with (H₂O)_n (n = 1-4) takes place, which will stimulate one to consider the favorable hydrolysis of water and water clusters with other Criegee intermediates and other important atmospheric species.

Atmospheric Chemistry of CH2OO: The Hydrolysis of CH2OO in Small Clusters of Sulfuric Acid

The hydrolysis of CH₂OO is not only a dominant sink for the CH₂OO intermediate in the atmosphere but also a key process in the formation of aerosols. Herein, the reaction mechanism and kinetics for the hydrolysis of CH₂OO catalyzed by the precursors of atmospheric aerosols, including H₂SO₄, H₂SO₄···H₂O, and (H₂SO₄)₂, have been studied theoretically at the CCSD(T)-F12a/cc-pVDZ-F12//B3LYP/6-311+G(2df,2pd) level. The calculated results show that the three catalysts decrease the energy barrier by over 10.3 kcal·mol^-1; at the same time, the product formation of HOCH₂OOH is more strongly bonded to the three catalysts than to the reactants CH₂OO and H₂O, revealing that small clusters of sulfuric acid promote the hydrolysis of CH₂OO both kinetically and thermodynamically. Kinetic simulations show that the H₂SO₄-assisted reaction is more favorable than the H₂SO₄···H₂O- (the pseudo-first-order rate constant being 27.9-11.5 times larger) and (H₂SO₄)₂- (between 2.8 × 10⁴ and 3.4 × 10⁵ times larger) catalyzed reactions. Additionally, due to relatively lower concentration of H₂SO₄, the hydrolysis of CH₂OO with H₂SO₄ cannot compete with the CH₂OO + H₂O or (H₂O)₂ reaction within the temperature range of 280-320 K, since its pseudo-first-order rate ratio is smaller by 4-7 or 6-8 orders of magnitude, respectively. However, the present results provide a good example of how small clusters of sulfuric acid catalyze the hydrolysis of an important atmospheric species.

Photodissociation Dynamics of CO-Forming Channels on the Ground-State Surface of Methyl Formate at 248 nm: Direct Dynamics Study and Assessment of Generalized Multicenter Impulsive Models

The photodissociation dynamics of methyl formate in the electronic ground state S₀, initiated by a 248 nm-wavelength laser, is studied by direct dynamics simulations. We analyze five channels, where four of them have as products CH₃OH + CO, one leading to the formation of three fragments, H₂CO + H₂ + CO, and a channel characterized by a roaming transition state. The analysis of energy distribution among the degrees of freedom of the product and the comparison with experimental results previously published by other groups provide the ingredients to distinguish the examined dissociation pathways. The interpretation of the results proves that the characterization of dissociation mechanisms must rely on a dynamics approach involving multiple electronic states, including considerations on the features of the S₁/S₀ conical intersection. Here, we also assess the generalized multicenter impulsive model, GMCIM, that has been designed for dissociation processes with exit barriers, and the energy distribution in the products is predicted on the basis of information from the saddle points and the intrinsic reaction coordinates. Main features, advantages, limits, and future perspectives of the method are reported and discussed.

Dynamics of the reaction CH2I + O2 probed via infrared emission of CO, CO2, OH and H2CO

The reaction CH2I + O2 has been widely employed recently for the production of the simplest Criegee intermediate CH2OO in laboratories, but the detailed dynamics of this reaction have been little explored. Infrared emission of several products of this reaction, initiated on irradiation of CH2I2 and O2 (∼8 Torr) in a flowing mixture at 308 or 248 nm, was recorded with a step-scan Fourier-transform spectrometer; possible routes of formation were identified according to the observed vibrational distribution of products and published theoretical potential-energy schemes. Upon irradiation at 308 nm, Boltzmann distributions of CO (v ≤ 5, J ≤ 19) with an average vibrational energy of 32 ± 3 kJ mol-1 and OH (v ≤ 3, J ≤ 5.5) with an average vibrational energy of 29 ± 4 kJ mol-1 were observed and assigned to the decomposition of HCOOH* to form CO + H2O and OH + HCO, respectively. The broadband emission of CO2 was simulated with two vibrational distributions of average energies (91 ± 4) and (147 ± 8) kJ mol-1 and assigned to be produced from the decomposition of HCOOH* and methylene bis(oxy), respectively. Upon irradiation of samples at 248 nm, the emission of OH and CO2 showed similar distributions with slightly greater energies, but the distribution of CO (v ≤ 11, J ≤ 19) became bimodal with average vibrational energies of (23 ± 4) and (107 ± 29) kJ mol-1, and branching (56 ± 5) : (44 ± 5). The additional large-v component is assigned to be produced from a secondary reaction HCO + O2 to form CO + HO2; HCO is a coproduct of OH. The branching between CO and OH is (50 ± 5) : (50 ± 5) at 308 nm and (64 ± 5) : (36 ± 4) at 248 nm, consistent with the mechanism according to which an additional channel to produce CO opens at 248 nm. Highly internally excited H2CO was also observed. With O2 at 16 Torr, the extrapolated nascent internal distributions are similar to those with O2 at 8 Torr except for a slight quenching effect.

Time-resolved, broadband UV-absorption spectrometry measurements of Criegee intermediate kinetics using a new photolytic precursor: unimolecular decomposition of CH2OO and its reaction with formic acid

We present a time-resolved broadband cavity-enhanced UV-absorption spectrometer apparatus that we have constructed and utilized for temperature- and pressure-dependent kinetic measurements of formaldehyde oxide (CH₂OO) reactions. We also introduce and utilize a new photolytic precursor, bromoiodomethane (CH₂IBr), which photolysis at 213 nm in presence of O₂ produces CH₂OO. Importantly, this precursor appears to be free from secondary reactions that may regenerate CH₂OO in kinetic experiments. The unimolecular decomposition rate coefficient of CH₂OO has been measured over wide pressure (5-400 Torr) and temperature (296-600 K) ranges and master equation simulations of the decomposition kinetics have been performed using MESMER program. The MESMER simulations of the experimental data with the calculated zero-point energy corrected transition state energy 85.9 kJ mol^-1 for decomposition required no adjustment and returned 〈ΔE〉_down = 123.2 × (T/298 K)^0.74 cm^-1 for temperature-dependent exponential-down model of the collisional energy transfer in He. A very good agreement between results of simulations and experiments is obtained. The results are compared with the previously reported unimolecular decomposition study by Stone et al. (Phys. Chem. Chem. Phys., 2018, 20, 24940-24954). Current master equation simulations suggest about 61% decomposition yield for the predominant H₂ + CO₂ channel, whereas the yields of two other channels, H₂O + CO, and HCO + OH, are sensitive on the parameters involved in the simulations. The kinetics of CH₂OO reaction with formic acid has also been investigated as function of pressure (5-150 Torr) and temperature (296-458 K). The bimolecular rate coefficient for CH₂OO + HCOOH reaction shows a negative temperature dependency, decreasing from (1.0 ± 0.03) × 10^-10 cm³ molecule^-1 s^-1 at 296 K to (0.47 ± 0.05) × 10^-10 cm³ molecule^-1 s^-1 at 458 K with an Arrhenius activation energy of -4.9 ± 1.6 kJ mol^-1, where statistical uncertainties shown are 2σ. Estimated overall uncertainty in the measured rate coefficients is about ±20%. Current bimolecular rate coefficient at room temperature agrees with the previously reported rate coefficients from the direct kinetic experiments. The reaction is found to be pressure independent over the range between 5 and 150 Torr at 296 K in He.

Four-Carbon Criegee Intermediate from Isoprene Ozonolysis: Methyl Vinyl Ketone Oxide Synthesis, Infrared Spectrum, and OH Production

The reaction of ozone with isoprene, one of the most abundant volatile organic compounds in the atmosphere, produces three distinct carbonyl oxide species (RR'COO) known as Criegee intermediates: formaldehyde oxide (CH₂OO), methyl vinyl ketone oxide (MVK-OO), and methacrolein oxide (MACR-OO). The nature of the substituents (R,R' = H, CH₃, CH═CH₂) and conformations of the Criegee intermediates control their subsequent chemistry in the atmosphere. In particular, unimolecular decay of MVK-OO is predicted to be the major source of hydroxyl radicals (OH) in isoprene ozonolysis. This study reports the initial laboratory synthesis and direct detection of MVK-OO through reaction of a photolytically generated, resonance-stabilized monoiodoalkene radical with O₂. MVK-OO is characterized utilizing infrared (IR) action spectroscopy, in which IR activation of MVK-OO with two quanta of CH stretch at ca. 6000 cm^-1 is coupled with ultraviolet detection of the resultant OH products. MVK-OO is identified by comparison of the experimentally observed IR spectral features with theoretically predicted IR absorption spectra. For syn-MVK-OO, the rate of appearance of OH products agrees with the unimolecular decay rate predicted using statistical theory with tunneling. This validates the hydrogen atom transfer mechanism and computed transition-state barrier (18.0 kcal mol^-1) leading to OH products. Theoretical calculations reveal an additional roaming pathway between the separating radical fragments, which results in other products. Master equation modeling yields a thermal unimolecular decay rate for syn-MVK-OO of 33 s^-1 (298 K, 1 atm). For anti-MVK-OO, theoretical exploration of several unimolecular decay pathways predicts that isomerization to dioxole is the most likely initial step to products.

Theories and simulations of roaming

The phenomenon of roaming in chemical reactions has now become both commonly observed in experiment and extensively supported by theory and simulations. Roaming occurs in highly-excited molecules when the trajectories of atomic motion often bypass the minimum energy pathway and produce reaction in unexpected ways from unlikely geometries. The prototypical example is the unimolecular dissociation of formaldehyde (H₂CO), in which the "normal" reaction proceeds through a tight transition state to yield H₂ + CO but for which a high fraction of dissociations take place via a "roaming" mechanism in which one H atom moves far from the HCO, almost to dissociation, and then returns to abstract the second H atom. We review below the theories and simulations that have recently been developed to address and understand this new reaction phenomenon.

Photochemistry of the Simplest Criegee Intermediate, CH2OO: Photoisomerization Channel toward Dioxirane Revealed by CASPT2 Calculations and Trajectory Surface-Hopping Dynamics

The photochemistry of Criegee intermediates plays a significant role in atmospheric chemistry, but it is relatively less explored compared with their thermal reactions. Using multireference CASPT2 electronic structure calculations and CASSCF trajectory surface-hopping molecular dynamics, we have revealed a dark-state-involved A¹A → X¹A photoisomerization channel of the simple Criegee intermediate (CH₂OO) that leads to a cyclic dioxirane. The excited molecules on the A¹A state, which can have either originated from the B¹A state via B¹A → A¹A internal conversion or formed by state-selective electronic excitation, is driven by the out-of-plane motion toward a perpendicular A/X¹A minimal-energy crossing point (MECI) then radiationless decay to the ground state with an average time constant of ∼138 fs, finally forming dioxirane at ∼254 fs. The dynamics starting from the A¹A state show that the quantum yield of photoisomerization from the simple Criegee intermediate to dioxirane is 38%. The finding of the A¹A → X¹A photoisomerization channel is expected to broaden the reactivity profile and deepen the understanding of the photochemistry of Criegee intermediates.

Selective deuteration illuminates the importance of tunneling in the unimolecular decay of Criegee intermediates to hydroxyl radical products

Ozonolysis of alkenes, an important nonphotolytic source of hydroxyl (OH) radicals in the atmosphere, proceeds through unimolecular decay of Criegee intermediates. Here, we report a large kinetic isotope effect associated with the rate-limiting hydrogen-transfer step that releases OH radicals for a prototypical Criegee intermediate, CH3CHOO. IR excitation of selectively deuterated syn-CD3CHOO is shown to result in deuterium atom transfer and release OD radical products. Vibrational activation of syn-CD3CHOO is coupled with direct time-resolved detection of OD products to measure a 10-fold slower rate of unimolecular decay upon deuteration in the vicinity of the transition state barrier, which is confirmed by microcanonical statistical theory that incorporates quantum mechanical tunneling. The corresponding kinetic isotope effect of ∼10 is attributed primarily to the decreased probability of D-atom vs. H-atom transfer arising from tunneling. Master equation modeling is utilized to compute the thermal unimolecular decay rates for selectively and fully deuterated syn methyl-substituted Criegee intermediates under atmospheric conditions. At 298 K (1 atm), tunneling is predicted to enhance the thermal decay rate of syn-CH3CHOO compared with the deuterated species, giving rise to a significant kinetic isotope effect of ∼50.

How does substitution affect the unimolecular reaction rates of Criegee intermediates?

Unimolecular reaction rates of Criegee intermediates show substitution effect.