Experimental and Modeling Studies of the Pressure and Temperature Dependences of the Kinetics and the OH Yields in the Acetyl + O2 Reaction

A three-channel thermal dissociation cavity ring-down spectrometer for simultaneous measurement of ambient total peroxy nitrates, total alkyl nitrates, and NO2

A newly constructed thermal dissociation cavity ring-down spectrometer (TD-CRDS) for the simultaneous measurement of ambient total peroxy nitrates (ΣPNs, RO2NO2), total alkyl nitrates (ΣANs, RONO2), and NO2 was presented in this work. ΣPNs and ΣANs were detected as NO2 with the CRDS instrument after thermal dissociation. PNs and ANs completely dissociated at 180 °C and 360 °C, with conversion efficiencies of 96 % and 99 %, respectively. The effects of NO2 and NO on measurement in different temperatures and two types of thermal dissociation inlet (TDI) were further explored. The influence of ambient NO2 and NO on PNs and ANs in the improved TDI (TDI-2) was significantly improved. To further enhance the measurement accuracy, the consistency of the observed NO2 in the three channels was tested, which achieved good agreement. The detection limits of the TD-CRDS instrument for NO2, ΣPNs, and ΣANs were determined as 6.5, 6.8, and 8.6 pptv (10 s, 1σ), respectively. Observations of PNs and ANs were conducted in a suburban site in Hefei, China, from September 2-30, 2021, using the TD-CRDS instrument, and the consecutive time series of PNs and ANs were derived, verifying the capability of the TD-CRDS instrument for continuous field observations of ΣPNs and ΣANs.

Comprehensive Theoretical Study on Four Typical Intramolecular Hydrogen Shift Reactions of Peroxy Radicals: Multireference Character, Recommended Model Chemistry, and Kinetics

Intramolecular hydrogen shift reactions in peroxy radicals (RO₂^• → ^•QOOH) play key roles in the low-temperature combustion and in the atmospheric chemistry. In the present study, we found that a mild-to-moderate multireference character of a potential energy surface (PES) is widely present in four typical hydrogen shift reactions of peroxy radicals (RO₂^•, R = ethyl, vinyl, formyl methyl, and acetyl) by a systematic assessment based on the T₁ diagnostic, %TAE diagnostic, M diagnostic, and contribution of the dominant configuration of the reference CASSCF wavefunction (C₀²). To assess the effects of these inherent multireference characters on electronic structure calculations, we compared the PESs of the four reactions calculated by the multireference method CASPT2 in the complete basis set (CBS) limit, single-reference method CCSD(T)-F12, and single-reference-based composite method WMS. The results showed that ignoring the multireference character will introduce a mean unsigned deviation (MUD) of 0.46-1.72 kcal/mol from CASPT2/CBS results by using the CCSD(T)-F12 method or a MUD of 0.49-1.37 kcal/mol by WMS for three RO₂^• reactions (R = vinyl, formyl methyl, and acetyl) with a stronger multireference character. Further tests by single-reference Kohn-Sham (KS) density functional theory methods showed even larger deviations. Therefore, we specifically developed a new hybrid meta-generalized gradient approximation (GGA) functional M06-HS for the four typical H-shift reactions of peroxy radicals based on the WMS results for the ethyl peroxy radical reaction and on the CASPT2/CBS results for the others. The M06-HS method has an averaged MUD of 0.34 kcal/mol over five tested basis sets against the benchmark PESs, performing best in the tested 38 KS functionals. Last, in a temperature range of 200-3000 K, with the new functional, we calculated the high-pressure-limit rate coefficients of these H-shift reactions by the multi-structural variational transition-state theory with the small-curvature tunneling approximation (MS-CVT/SCT) and the thermochemical properties of all of the involved key radicals by the multi-structural torsional (MS-T) anharmonicity approximation method.

Investigating the Association Reactions of HOCH2CO and HOCHCHO with O2: A Quantum Computational and Master Equation Study

Glycolaldehyde, HOCH₂CHO, is an important multifunctional atmospheric trace gas formed in the oxidation of ethylene and isoprene and emitted directly from burning biomass. The initial step in the atmospheric photooxidation of HOCH₂CHO yields HOCH₂CO and HOCHCHO radicals; both of these radicals react rapidly with O₂ in the troposphere. This study presents a comprehensive theoretical investigation of the HOCH₂CO + O₂ and HOCHCHO + O₂ reactions using high-level quantum chemical calculations and energy-grained master equation simulations. The HOCH₂CO + O₂ reaction results in the formation of a HOCH₂C(O)O₂ radical, while the HOCHCHO + O₂ reaction yields (HCO)₂ + HO₂. Density functional theory calculations have identified two open unimolecular pathways associated with the HOCH₂C(O)O₂ radical that yield HCOCOOH + OH or HCHO + CO₂ + OH products; the former novel bimolecular product pathway has not been previously reported in the literature. Master equation simulations based on the potential energy surface calculated here for the HOCH₂CO + O₂ recombination reaction support experimental product yield data from the literature and indicate that, even at total pressures of 1 atm, the HOCH₂CO + O₂ reaction yields ∼11% OH at 298 K.

Infrared spectroscopic signature of a hydroperoxyalkyl radical (•QOOH)

Infrared (IR) action spectroscopy is utilized to characterize a prototypical carbon-centered hydroperoxyalkyl radical (•QOOH) transiently formed in the oxidation of volatile organic compounds. The •QOOH radical formed in isobutane oxidation, 2-hydroperoxy-2-methylprop-1-yl, •CH₂(CH₃)₂COOH, is generated in the laboratory by H-atom abstraction from tert-butyl hydroperoxide (TBHP). IR spectral features of jet-cooled and stabilized •QOOH radicals are observed from 2950 to 7050 cm^-1 at energies that lie below and above the transition state barrier leading to OH radical and cyclic ether products. The observed •QOOH features include overtone OH and CH stretch transitions, combination bands involving OH or CH stretch and a lower frequency mode, and fundamental OH and CH stretch transitions. Most features arise from a single vibrational transition with band contours well simulated at a rotational temperature of 10 K. In each case, the OH products resulting from unimolecular decay of vibrationally activated •QOOH are detected by UV laser-induced fluorescence. Assignments of observed •QOOH IR transitions are guided by anharmonic frequencies computed using second order vibrational perturbation theory, a 2 + 1 model that focuses on the coupling of the OH stretch with two low-frequency torsions, as well as recently predicted statistical •QOOH unimolecular decay rates that include heavy-atom tunneling. Most of the observed vibrational transitions of •QOOH are readily distinguished from those of the TBHP precursor. The distinctive IR transitions of •QOOH, including the strong fundamental OH stretch, provide a general means for detection of •QOOH under controlled laboratory and real-world conditions.

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.

Atmospheric Chemistry of Allylic Radicals from Isoprene: A Successive Cyclization-Driven Autoxidation Mechanism

The atmospheric chemistry of isoprene has broad implications for regional air quality and the global climate. Allylic radicals, taking 13-17% yield in the isoprene oxidation by ^•Cl, can contribute as much as 3.6-4.9% to all possible formed intermediates in local regions at daytime. Considering the large quantity of isoprene emission, the chemistry of the allylic radicals is therefore highly desirable. Here, we investigated the atmospheric oxidation mechanism of the allylic radicals using quantum chemical calculations and kinetics modeling. The results indicate that the allylic radicals can barrierlessly combine with O₂ to form peroxy radicals (RO₂^•). Under ≤100 ppt NO and ≤50 ppt HO₂^• conditions, the formed RO₂^• mainly undergo two times "successive cyclization and O₂ addition" to finally form the product fragments 2-alkoxy-acetaldehyde (C₂H₃O₂^•) and 3-hydroperoxy-2-oxopropanal (C₃H₄O₄). The presented reaction illustrates a novel successive cyclization-driven autoxidation mechanism. The formed 3-hydroperoxy-2-oxopropanal product is a new isomer of the atmospheric C₃H₄O₄ family and a potential aqueous-phase secondary organic aerosol precursor. Under >100 ppt NO condition, NO can mediate the cyclization-driven autoxidation process to form C₅H₇NO₃, C₅H₇NO₇, and alkoxy radical-related products. The proposed novel autoxidation mechanism advances our current understanding of the atmospheric chemistry of both isoprene and RO₂^•.

Experimental Evidence of Dioxole Unimolecular Decay Pathway for Isoprene-Derived Criegee Intermediates

Ozonolysis of isoprene, one of the most abundant volatile organic compounds emitted into the Earth's atmosphere, generates two four-carbon unsaturated Criegee intermediates, methyl vinyl ketone oxide (MVK-oxide) and methacrolein oxide (MACR-oxide). The extended conjugation between the vinyl substituent and carbonyl oxide groups of these Criegee intermediates facilitates rapid electrocyclic ring closures that form five-membered cyclic peroxides, known as dioxoles. This study reports the first experimental evidence of this novel decay pathway, which is predicted to be the dominant atmospheric sink for specific conformational forms of MVK-oxide (anti) and MACR-oxide (syn) with the vinyl substituent adjacent to the terminal O atom. The resulting dioxoles are predicted to undergo rapid unimolecular decay to oxygenated hydrocarbon radical products, including acetyl, vinoxy, formyl, and 2-methylvinoxy radicals. In the presence of O₂, these radicals rapidly react to form peroxy radicals (ROO), which quickly decay via carbon-centered radical intermediates (QOOH) to stable carbonyl products that were identified in this work. The carbonyl products were detected under thermal conditions (298 K, 10 Torr He) using multiplexed photoionization mass spectrometry (MPIMS). The main products (and associated relative abundances) originating from unimolecular decay of anti-MVK-oxide and subsequent reaction with O₂ are formaldehyde (88 ± 5%), ketene (9 ± 1%), and glyoxal (3 ± 1%). Those identified from the unimolecular decay of syn-MACR-oxide and subsequent reaction with O₂ are acetaldehyde (37 ± 7%), vinyl alcohol (9 ± 1%), methylketene (2 ± 1%), and acrolein (52 ± 5%). In addition to the stable carbonyl products, the secondary peroxy chemistry also generates OH or HO₂ radical coproducts.

Reaction between Peroxy and Alkoxy Radicals Can Form Stable Adducts

Peroxy (RO2) and alkoxy (RO) radicals are prototypical intermediates in any hydrocarbon oxidation. In this work, we use computational methods to (1) study the mechanism and kinetics of the RO2 + OH reaction for previously unexplored "R" structures (R = CH(O)CH2 and R = CH3C(O)) and (2) investigate a hitherto unaccounted channel of molecular growth, R'O2 + RO. On the singlet surface, these reactions rapidly form ROOOH and R'OOOR adducts, respectively. The former decomposes to RO + HO2 and R(O)OH + O2 products, while the main decomposition channel for the latter is back to the reactant radicals. Decomposition rates of R'OOOR adducts varied between 103 and 0.015 s-1 at 298 K and 1 atm. The most long-lived R'OOOR adducts likely account for some fraction of the elemental compositions detected in the atmosphere that are commonly assigned to stable covalently bound dimers.

Observation of a new channel, the production of CH3, in the abstraction reaction of OH radicals with acetaldehyde

Using laser flash photolysis coupled to photo-ionization time-of-flight mass spectrometry (PIMS), methyl radicals (CH₃) have been detected as primary products from the reaction of OH radicals with acetaldehyde (ethanal, CH₃CHO) with a yield of ∼15% at 1-2 Torr of helium bath gas. Supporting measurements based on laser induced fluorescence studies of OH recycling in the OH/CH₃CHO/O₂ system are consistent with the PIMS study. Master equation calculations suggest that the origin of the methyl radicals is from prompt dissociation of chemically activated acetyl products and hence is consistent with previous studies which have shown that abstraction, rather than addition/elimination, is the sole route for the OH + acetaldehyde reaction. However, the observation of a significant methyl product yield suggests that energy partitioning in the reaction is different from the typical early barrier mechanism where reaction exothermicity is channeled preferentially into the newly formed bond. The master equation calculations predict atmospheric yields of methyl radicals of ∼9%. The implications of the observations in atmospheric and combustion chemistry are briefly discussed.

CH3CO + O2 + M (M = He, N2) Reaction Rate Coefficient Measurements and Implications for the OH Radical Product Yield

The gas-phase CH3CO + O2 reaction is known to proceed via a chemical activation mechanism leading to the formation of OH and CH3C(O)OO radicals via bimolecular and termolecular reactive channels, respectively. In this work, rate coefficients, k, for the CH3CO + O2 reaction were measured over a range of temperature (241-373 K) and pressure (0.009-600 Torr) with He and N2 as the bath gas and used to characterize the bi- and ter-molecular reaction channels. Three independent experimental methods (pulsed laser photolysis-laser-induced fluorescence (PLP-LIF), pulsed laser photolysis-cavity ring-down spectroscopy (PLP-CRDS), and a very low-pressure reactor (VLPR)) were used to characterize k(T,M). PLP-LIF was the primary method used to measure k(T,M) in the high-pressure regime under pseudo-first-order conditions. CH3CO was produced by PLP, and LIF was used to monitor the OH radical bimolecular channel reaction product. CRDS, a complementary high-pressure method, measured k(295 K,M) over the pressure range 25-600 Torr (He) by monitoring the temporal CH3CO radical absorption following its production via PLP in the presence of excess O2. The VLPR technique was used in a relative rate mode to measure k(296 K,M) in the low-pressure regime (9-32 mTorr) with CH3CO + Cl2 used as the reference reaction. A kinetic mechanism analysis of the combined kinetic data set yielded a zero pressure limit rate coefficient, kint(T), of (6.4 ± 4) × 10(-14) exp((820 ± 150)/T) cm(3) molecule(-1) s(-1) (with kint(296 K) measured to be (9.94 ± 1.3) × 10(-13) cm(3) molecule(-1) s(-1)), k0(T) = (7.39 ± 0.3) × 10(-30) (T/300)(-2.2±0.3) cm(6) molecule(-2) s(-1), and k∞(T) = (4.88 ± 0.05) × 10(-12) (T/300)(-0.85±0.07) cm(3) molecule(-1) s(-1) with Fc = 0.8 and M = N2. A He/N2 collision efficiency ratio of 0.60 ± 0.05 was determined. The phenomenological kinetic results were used to define the pressure and temperature dependence of the OH radical yield in the CH3CO + O2 reaction. The present results are compared with results from previous studies and the discrepancies are discussed.

Global Uncertainty Propagation and Sensitivity Analysis in the CH3OCH2 + O2 System: Combining Experiment and Theory To Constrain Key Rate Coefficients in DME Combustion

Statistical rate theory calculations, in particular formulations of the chemical master equation, are widely used to calculate rate coefficients of interest in combustion environments as a function of temperature and pressure. However, despite the increasing accuracy of electronic structure calculations, small uncertainties in the input parameters for these master equation models can lead to relatively large uncertainties in the calculated rate coefficients. Master equation input parameters may be constrained further by using experimental data and the relationship between experiment and theory warrants further investigation. In this work, the CH3OCH2 + O2 system, of relevance to the combustion of dimethyl ether (DME), is used as an example and the input parameters for master equation calculations on this system are refined through fitting to experimental data. Complementing these fitting calculations, global sensitivity analysis is used to explore which input parameters are constrained by which experimental conditions, and which parameters need to be further constrained to accurately predict key elementary rate coefficients. Finally, uncertainties in the calculated rate coefficients are obtained using both correlated and uncorrelated distributions of input parameters.

Branching ratios for the reactions of OH with ethanol amines used in carbon capture and the potential impact on carcinogen formation in the emission plume from a carbon capture plant

Branching ratios for the OH reaction with ethanol amines and potential risk of carcinogenic formation in the carbon capture plume.

Branching ratios in reactions of OH radicals with methylamine, dimethylamine, and ethylamine

The branching ratios for the reaction of the OH radical with the primary and secondary alkylamines: methylamine (MA), dimethylamine (DMA), and ethylamine (EA), have been determined using the technique of pulsed laser photolysis-laser-induced fluorescence. Titration of the carbon-centered radical, formed following the initial OH abstraction, with oxygen to give HO2 and an imine, followed by conversion of HO2 to OH by reaction with NO, resulted in biexponential OH decay traces on a millisecond time scale. Analysis of the biexponential curves gave the HO2 yield, which equaled the branching ratio for abstraction at αC-H position, r(αC-H). The technique was validated by reproducing known branching ratios for OH abstraction for methanol and ethanol. For the amines studied in this work (all at 298 K): r(αC-H,MA) = 0.76 ± 0.08, r(αC-H,DMA) = 0.59 ± 0.07, and r(αC-H,EA) = 0.49 ± 0.06 where the errors are a combination in quadrature of statistical errors at the 2σ level and an estimated 10% systematic error. The branching ratios r(αC-H) for OH reacting with (CH3)2NH and CH3CH2NH2 are in agreement with those obtained for the OD reaction with (CH3)2ND (d-DMA) and CH3CH2ND2 (d-EA): r(αC-H,d-DMA) = 0.71 ± 0.12 and r(αC-H,d-EA) = 0.54 ± 0.07. A master equation analysis (using the MESMER package) based on potential energy surfaces from G4 theory was used to demonstrate that the experimental determinations are unaffected by formation of stabilized peroxy radicals and to estimate atmospheric pressure yields. The branching ratio for imine formation through the reaction of O2 with α carbon-centered radicals at 1 atm of N2 are estimated as r(CH2NH2) = 0.79 ± 0.15, r(CH2NHCH3) = 0.72 ± 0.19, and r(CH3CHNH2) = 0.50 ± 0.18. The implications of this work on the potential formation of nitrosamines and nitramines are briefly discussed.

Direct kinetic study of OH and O3 formation in the reaction of CH3C(O)O2 with HO2

The reaction between HO2 and CH3C(O)O2 has three exothermic product channels, forming OH (R3a), peracetic acid (R3b), and acetic acid plus O3 (R3c). The branching ratios of the OH- and ozone-forming reaction channels were determined using a combination of laser-induced fluorescence (LIF, for time-resolved OH concentration measurement) and transient absorption spectroscopy (TAS, for time-resolved O3 concentration measurement) following pulsed laser generation of HO2 and CH3C(O)O2 from suitable precursors. TAS was also used to determine the initial concentration of the reactant peroxy radicals. The data were evaluated by numerical simulation using kinetic models of the measured concentration profiles; a Monte Carlo approach was used to estimate the uncertainties of the rate constants (k3) and branching ratios (α) thus obtained. The reaction channel forming OH (R3a) was found to be the most important with α3a = 0.61 ± 0.09 and α3c = 0.16 ± 0.08. The overall rate coefficient of the title reaction was found to be k3 = (2.1 ± 0.4) × 10(-11) cm(3) molecule(-1) s(-1) for both HO2 and DO2. Use of DO2 resulted in an increase in α3a to 0.80 ± 0.14. Comparison with former studies shows that OH formation via (R3) has been underestimated significantly to date. Possible reasons for these discrepancies and atmospheric implications are discussed.

The fast C(3P) + CH3OH reaction as an efficient loss process for gas-phase interstellar methanol

The reaction between ground state atomic carbon and methanol is shown to be an efficient destruction mechanism for interstellar methanol.

Interception of excited vibrational quantum states by O2 in atmospheric association reactions

Bimolecular reactions in Earth's atmosphere are generally assumed to proceed between reactants whose internal quantum states are fully thermally relaxed. Here, we highlight a dramatic role for vibrationally excited bimolecular reactants in the oxidation of acetylene. The reaction proceeds by preliminary adduct formation between the alkyne and OH radical, with subsequent O(2) addition. Using a detailed theoretical model, we show that the product-branching ratio is determined by the excited vibrational quantum-state distribution of the adduct at the moment it reacts with O(2). Experimentally, we found that under the simulated atmospheric conditions O(2) intercepts ~25% of the excited adducts before their vibrational quantum states have fully relaxed. Analogous interception of excited-state radicals by O(2) is likely common to a range of atmospheric reactions that proceed through peroxy complexes.