Insights into the Reaction Mechanism of Criegee Intermediate CH2OO with Methane and Implications for the Formation of Methanol

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.

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.

Mechanistic and Kinetic Approach on the Propargyl Radical (C3H3) with the Criegee Intermediate (CH2OO)

The detailed reaction mechanism and kinetics of the C₃H₃ + CH₂OO system have been thoroughly investigated. The CBS-QB3 method in conjunction with the ME/vRRKM theory has been applied to figure out the potential energy surface and rate constants for the C₃H₃ + CH₂OO system. The C₃H₃ + CH₂OO reaction leading to the CH₂-[cyc-CCHCHOO] + H product dominates compared to the others. Rate constants of the reaction are dependent on temperatures (300-2000 K) and pressures (1-76,000 Torr), for which the rate constant of the channel C₃H₃ + CH₂OO → CH₂-[cyc-CCHCHOO] + H decreases at low pressures (1-76 Torr), but it increases with rising temperature if the pressure P ≥ 760 Torr. Rate constants of the three reaction channels C₃H₃ + CH₂OO → CHCCH₂CHO + OH, C₃H₃ + CH₂OO → OCHCHCHCHO + H, and C₃H₃ + CH₂OO → CHCHCHO + CH₂O fluctuate with temperatures. The branching ratio of the C₃H₃ + CH₂OO → CH₂-[cyc-CCHCHOO] + H channel is the highest, accounting for 51-98.7% in the temperature range of 300-2000 K and 760 Torr pressure, while those of the channels forming the products PR10 (OCHCHCHCHO + H) and PR11 (CHCHCHO + CH₂O) are the lowest, less than 0.1%, indicating that the contribution of these two reaction paths to the title reaction is insignificant. The proposed temperature- and pressure-dependent rate constants, together with the thermodynamic data of the species involved, can be confidently used for modeling CH₂OO-related systems under atmospheric and combustion conditions.

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.

Energetics and electronics of polar Diels-Alder reactions at the atomic level: QTAIM and IQA analyses of complete IRC paths

The mechanism of Diels-Alder reactions between cyclopentadiene and several cyanoethylenes was studied by means of Density Functional Theory calculations using QTAIM and IQA (Interacting Quantum Atoms) analyses along complete IRC paths. Each geometry from the IRC had its wavefunction computed and the topology of the electronic density for it was then evaluated. By means of IQA, the global energetic profile was partitioned among the various atoms in the molecule, providing insight into what atoms are the main ones responsible for the magnitude of the energy barriers. The (a)synchronicity of the reaction mechanisms featuring non-symmetrically substituted dienophiles was characterized, from QTAIM, by the electron densities and Laplacians over the LCP's as well as by the different atomic energy barriers obtained from IQA. The magnitude of the atomic barrier nicely explains the (a)synchronicity of the reaction mechanisms, and the degree of (a)synchronicity is nicely revealed by the difference between the earlier and later bond breaking and bond formations. The main conclusion is that important energetic and electronic changes are occurring before and after the position of the transition state structure, mainly for those asynchronous mechanisms, and although these provide essential insight into the reaction mechanism, most studies cannot assess this kind of information because they are focusing solely on reactants, transition states, and products. We advocate that the additional computational effort required for such analyses is more than compensated by the great amount of useful information it provides.

Theoretical study on the cycloaddition reaction of phosphenium cation and formaldehyde

The mechanism of cycloaddition reaction between phosphenium cation and phosphindene with formaldehyde has been systematically investigated at the B3LYP/6-311++G(d,p) level of theory in order to better understand the reactivity for the valence isoelectronic of carbene. Phosphenium cation acts as an electrophilic reagent and accepts σ electrons of formaldehyde to form a complex in the first addition step. The greater the positive charge on phosphorus atom in phosphenium cation, the more stable the formed complex is. Introduction of substituents will decrease positive charge on phosphorus atom in phosphenium cation. The order of positive charge on phosphorus atom is HP+–F >  HP+–OH >  HP+–NH2, which is consistent with their Lewis acidities. The complex transforms to a three-membered ring product via a transition state in the second cyclization step. The product is less stable than the complex due to its tension of small ring.

Computational mechanistic study of the unimolecular dissociation of ethyl hydroperoxide and its bimolecular reactions with atmospheric species

A detailed computational study of the atmospheric reaction of the simplest Criegee intermediate CH₂OO with methane has been performed using the density functional theory (DFT) method and high-level calculations. Solvation models were utilized to address the effect of water molecules on prominent reaction steps and their associated energies. The structures of all proposed mechanisms were optimized using B3LYP functional with several basis sets: 6-31G(d), 6-31G (2df,p), 6-311++G(3df,3pd) and at M06-2X/6-31G(d) and APFD/6-31G(d) levels of theory. Furthermore, all structures were optimized at the B3LYP/6-311++G(3df,3pd) level of theory. The intrinsic reaction coordinate (IRC) analysis was performed for characterizing the transition states on the potential energy surfaces. Fifteen different mechanistic pathways were studied for the reaction of Criegee intermediate with methane. Both thermodynamic functions (ΔH and ΔG), and activation parameters (activation energies E_a, enthalpies of activation ΔH^ǂ, and Gibbs energies of activation ΔG^ǂ) were calculated for all pathways investigated. The individual mechanisms for pathways A1, A2, B1, and B2, comprise two key steps: (i) the formation of ethyl hydroperoxide (EHP) accompanying with the hydrogen transfer from the alkanes to the terminal oxygen atom of CIs, and (ii) a following unimolecular dissociation of EHP. Pathways from C1 → H1 involve the bimolecular reaction of EHP with different atmospheric species. The photochemical reaction of methane with EHP (pathway E1) was found to be the most plausible reaction mechanism, exhibiting an overall activation energy of 7 kJ mol⁻¹, which was estimated in vacuum at the B3LYP/6-311++G(3df,3pd) level of theory. All of the reactions were found to be strongly exothermic, expect the case of the sulfur dioxide-involved pathway that is predicted to be endothermic. The solvent effect plays an important role in the reaction of EHP with ammonia (pathway F1). Compared with the gas phase reaction, the overall activation energy for the solution phase reaction is decreased by 162 and 140 kJ mol⁻¹ according to calculations done with the SMD and PCM solvation models, respectively.

Hydrogen bond, ring tension and π-conjugation effects: methyl and vinyl substitutions dramatically change the photodynamics of Criegee intermediates

The substitution effect in chemistry is a concept that is probably too common to mention, while for a molecule with an elusive electronic structure, substitution can introduce an unusual effect that dramatically tunes the chemical process. To reveal the substitution effects on the photodynamics of Criegee Intermediates (CIs), we carried out the multireference CASSCF trajectory surface-hopping (TSH) molecular dynamics and CASPT2 electronic-structure calculations for a methyl-substituted CI (MCI) and a vinyl-substituted CI (VCI). The results show that for different substituents, the hydrogen bond, ring tension and π-conjugation not only alter the relative stabilities of the conformers/configurations, but also dramatically change the photo-induced channel of CIs. For an anti-MCI, the dominant channel starting from the S₁ state is the ring-closure process leading to dioxirane, while in the syn configuration, the intramolecular CHO hydrogen bond hinders the rotation around the C-O bond and thus leads to a high yield of in-plane O-O dissociation towards acetaldehyde (X¹A') and the O(¹D) atom. In a VCI with an unsaturated substituent, the π-conjugation greatly strengthens the O-O bond and therefore no O-O dissociation is observed in all configurations. In addition, the CHO hydrogen bond in the syn_(CO)-VCI further stabilizes the S₁-state intermediates and makes them less reactive; in contrast, isomerization to dioxirane becomes the globally dominant channel in the anti_(CO)-VCI. The dramatic substitution effects by saturated and unsaturated substituents on CIs found here will deepen the understanding of Criegee-intermediate chemistry.

Theoretical study on the reaction between silacyclopropenylidene and three-membered heterocyclic compounds (azirane and oxirane): An alternative approach to the formation of heterocyclic silylene

The reaction mechanism between silacyclopropenylidene and three-membered heterocyclic compounds (azirane and oxirane) has been systematically investigated at the B3LYP/6-311+G* level of theory in order to better understand the reactivity of unsaturated cyclic silylene. Geometry optimizations and vibrational analyses have been conducted for the stationary points on the potential energy surface of the system. Calculations show that the Si-spiroheterocyclic intermediate and four-membered heterocyclic silylene compound could be produced through the insertion process and subsequent dissociation process between silacyclopropenylidene and three-membered heterocyclic compounds. For the insertion process, it is easier for silacyclopropenylidene to insert into C-N bond of azirane than into C-O bond of oxirane. This study is helpful to understand the reactivity of silacyclopropenylidene, the evolution of silicon-bearing molecules in space, and to offer an alternative approach to the formation of enlarged heterocyclic silylene compound.

The atmospheric importance of methylamine additions to Criegee intermediates

The methylamine addition to Criegee intermediates is investigated using high level ab initio methods.

Theoretical Insights into the Reaction Mechanism between 2,3,7,8-Tetrachlorodibenzofuran and Hydrogen Peroxide: A DFT Study

A detailed knowledge of the reactivity of 2,3,7,8-tetrachlorodibenzofuran (TCDF) at the molecular level is important to better understand the transformation of dioxins analogous to TCDF in the environment. To clarify the reactivity of the organic hydroperoxides toward TCDF, the reaction of the TCDF with hydrogen peroxide (H2O2) and its anion has been investigated theoretically. For the reaction of the neutral H2O2, a molecular complex can be formed between TCDF and H2O2 first. Then, the nucleophilic aromatic substitution of TCDF by H2O2 occurs in the presence of the water molecules to form an intermediate containing an O-O bond. Finally, the O-O bond cleavages homolytically for the above intermediate. On the other hand, as for the reaction of the anion of H2O2 (HO2 -), the nucleophilic addition of HO2 - to TCDF can also occur besides the nucleophilic aromatic substitution reaction mentioned above, resulting in the dissociation of the C-O bond of TCDF. Unlike the reaction involving neutral H2O2, no water molecules are required. In addition, the selected substitution effects, such as F-, Br-, and CH3-substituents, on the reactivity of the above reaction have also been explored. Hopefully, the present results can enable us to gain insights into the reactivity of the organic hydroperoxides with TCDF-like environmental pollutants.

Theoretical investigation on the reaction mechanism and kinetics of a Criegee intermediate with ethylene and acetylene

The competition among the possible pathways, the branching ratios of the adduct and the decomposition products at different temperatures and pressures have been evaluated.

An Extended Computational Study of Criegee Intermediate-Alcohol Reactions

High-level ab initio calculations (DF-LCCSD(T)-F12a//B3LYP/aug-cc-pVTZ) are performed on a range of stabilized Criegee intermediate (sCI)-alcohol reactions, computing reaction coordinate energies, leading to the formation of α-alkoxyalkyl hydroperoxides (AAAHs). These potential energy surfaces are used to model bimolecular reaction kinetics over a range of temperatures. The calculations performed in this work reproduce the complicated temperature-dependent reaction rates of CH₂OO and (CH₃)₂COO with methanol, which have previously been experimentally determined. This methodology is then extended to compute reaction rates of 22 different Criegee intermediates with methanol, including several intermediates derived from isoprene ozonolysis. In some cases, sCI-alcohol reaction rates approach those of sCI-(H₂O)₂. This suggests that in regions with elevated alcohol concentrations, such as urban Brazil, these reactions may generate significant quantities of AAAHs and may begin to compete with sCI reactions with other trace tropospheric pollutants such as SO₂. This work also demonstrates the ability of alcohols to catalyze the 1,4-H transfer unimolecular decomposition of α-methyl substituted sCIs.

Theoretical Investigations on the Reactivity of Hydrogen Peroxide toward 2,3,7,8-Tetrachlorodibenzo-p-dioxin

Acquiring full knowledge of the reactivity of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) is crucial for the better understanding of the transformation and degradation of TCDD-like dioxins in the environment. To clarify the reactivity of the organic hydroperoxides toward TCDD, in this study, the reactions between the neutral/anion of the hydrogen peroxide (H₂O₂) and TCDD have been systematically investigated theoretically. It was found that the neutral H₂O₂ is relatively difficult to react with TCDD compared with its anion, exhibiting the pH dependence of the title reaction. As for the anion of H₂O₂, it reacts with TCDD through two reaction mechanisms, i.e., nucleophilic substitution and nucleophilic addition. For the former, the terminal O atom of HO₂- nucleophilically attacks the C atom of the C-Cl bond in TCDD to form an intermediate containing an O-O bond, accompanying the dissociation of the chlorine atom. For the latter, the terminal O atom of HO₂- can be easily attached to the C atom of the C-O bond in TCDD, resulting in the decomposition of C-O bond and the formation of an intermediate containing an O-O bond. For these formed intermediates in both reaction mechanisms, their O-O bonds can be homolytically cleaved to produce different radicals. In addition, the selected substitution effects including F-, Br-, and CH₃- substituents on the above reactions have also been studied. Hopefully, the present results can provide new insights into the reactivity of the organic hydroperoxides toward TCDD-like environmental pollutants.

Theoretical Investigations on the Reactivity of Methylidyne Radical toward 2,3,7,8-Tetrachlorodibenzo-p-Dioxin: A DFT and Molecular Dynamics Study

To explore the potential reactivity of the methylidyne radical (CH) toward 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), the reaction mechanism between them has been systematically investigated employing the density functional theory (DFT) and ab initio molecular dynamics simulations. The relevant thermodynamic and kinetic parameters in the possible reaction pathways have been discussed as well as the IR spectra and hyperfine coupling constants (hfcc's) of the major products. Different from the reaction of the CH radical with 2,3,7,8-tetrachlorodibenzofuran, CH radical can attack all the C-C bonds of TCDD to form an initial intermediate barrierlessly via the cycloaddition mechanism. After then, the introduced C-H bond can be further inserted into the C-C bond of TCDD, resulting in the formation of a seven-membered ring structure. The whole reactions are favorable thermodynamically and kinetically. Moreover, the major products have been verified by ab initio molecular dynamics simulations. The distinct IR spectra and hyperfine coupling constants of the major products can provide some help for their experimental detection and identification. In addition, the reactivity of the CH radical toward the F- and Br-substituted TCDDs has also been investigated. Hopefully, the present findings can provide new insights into the reactivity of the CH radical in the transformation of TCDD-like dioxins.

Theoretical Insights into the Interaction Mechanisms between Nitric Acid and Nitrous Oxide Initiated by an Excess Electron

Nitric acid (HNO₃) and nitrous oxide (N₂O) play an important role in the atmospheric chemistry in regulating the global environment and climate changes. In this study, the interaction mechanisms between them have been systematically investigated before and after the electron capture employing the density functional theory in combination with the AIM, NBO, and ab initio molecular dynamics calculations. It was found that HNO₃ and N₂O can form transient complexes through intermolecular H-bonds. HNNO, OH, and NO₂ free radicals can be produced after the electron capture of the formed complexes, providing an alternative source of these radicals in the atmosphere. The present results not only can provide new insights into the transformation of the HNO₃ and N₂O atmospheric species but also can enable us to better understand the potential role of the free electron in the atmosphere.

Theoretical Insights into the Electron Capture Behavior of H₂SO₄···N₂O Complex: A DFT and Molecular Dynamics Study

Both sulfuric acid (H₂SO₄) and nitrous oxide (N₂O) play a central role in the atmospheric chemistry in regulating the global environment and climate changes. In this study, the interaction behavior between H₂SO₄ and N₂O before and after electron capture has been explored using the density functional theory (DFT) method as well as molecular dynamics simulation. The intermolecular interactions have been characterized by atoms in molecules (AIM), natural bond orbital (NBO), and reduced density gradient (RDG) analyses, respectively. It was found that H₂SO₄ and N₂O can form two transient molecular complexes via intermolecular H-bonds within a certain timescale. However, two molecular complexes can be transformed into OH radical, N₂, and HSO₄⁻ species upon electron capture, providing an alternative formation source of OH radical in the atmosphere. Expectedly, the present findings not only can provide new insights into the transformation behavior of H₂SO₄ and N₂O, but also can enable us to better understand the potential role of the free electron in driving the proceeding of the relevant reactions in the atmosphere.

Theoretical insights into the reaction mechanisms between 2,3,7,8-tetrachlorodibenzofuran and the methylidyne radical

To explore the potential role of the methylidyne radical (CH) in the transformation of 2,3,7,8-tetrachlorodibenzofuran (TCDF), in this study, the detailed reaction mechanisms between TCDF and CH radical have been systematically investigated employing the B3LYP method of density functional theory (DFT) in combination with the atoms in molecules (AIM) theory and ab initio molecular dynamics. It was found that the title reaction is a multi-channel reaction, i.e., the CH radical can attack the C-X (X = C, Cl, H, O) bonds of TCDF via the insertion modes, resulting in the formation of 13 products. Thermodynamically, the whole reaction processes are exothermic and spontaneous since all the enthalpy and Gibbs free energy changes are negative values in the formation processes. Moreover, the thermodynamic stability of the products is controlled by the distribution of the single unpaired electron. Kinetically, the most favorable reaction channel is the insertion of the CH radical into the C-C bond except for the C atoms attached to the chlorine atom. Moreover, the dominant products have been further confirmed by the molecular dynamics. Meanwhile, the IR spectra and hyperfine coupling constants of the dominant products have been investigated to provide helpful information for their identification experimentally. In addition, the reactivity of the CH radical toward the F- and Br-substituted TCDFs has also been investigated. Expectedly, the present findings can enable us to better understand the reactivity of the CH radical toward organic pollutants analogous to TCDF in the atmosphere.

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.