Temperature-Dependent, Site-Specific Rate Coefficients for the Reaction of OH (OD) with Methyl Formate Isotopologues via Experimental and Theoretical Studies.

Methyl esters are an important component of combustion and atmospheric systems. Reaction with the OH radical plays an important role in the removal of the simplest methyl ester, methyl formate (MF, CH3OCHO). In this paper, the overall rate coefficients for the reactions of OH and OD with MF isotopologues, studied under pseudo-first-order conditions, are reported using two different laser flash photolysis systems with the decay of OH monitored by laser-induced fluorescence. The room-temperature rate coefficient for OH + MF, (1.95 ± 0.34) × 10-13 cm3 molecule-1 s-1, is in good agreement with the literature. The rate coefficient exhibits curved Arrhenius behavior, and our results bridge the gap between previous low-temperature and shock tube studies. In combination with the literature, the rate coefficient for the reaction of OH with MF between 230 and 1400 K can be parametrized as kOH+MF = (3.2 × 10-13) × (T/300 K)2.3 × exp(-141.4 K/T) cm3 molecule-1 s-1 with an overall estimated uncertainty of ∼30%. The reactions of OD with MF isotopologues show a small enhancement (inverse secondary isotope effect) compared to the respective OH reactions. The reaction of OH/OD with MF shows a normal primary isotope effect, a decrease in the rate coefficient when MF is partially or fully deuterated. Experimental studies have been supported by ab initio calculations at the CCSD(T)-F12/aug-cc-pVTZ//M06-2X/6-31+G** level of theory. The calculated, zero-point-corrected, barrier heights for abstraction at the methyl and formate sites are 1.3 and 6.0 kJ mol-1, respectively, and the ab initio predictions of kinetic isotope effects are in agreement with experiment. Fitting the experimental isotopologue data refines these barriers to 0.9 ± 0.6 and 4.1 ± 0.9 kJ mol-1. The branching ratio is approximately 50:50 at 300 K. Between 300 and 500 K, abstraction via the higher-energy, higher-entropy formate transition state becomes more favored (60:40). However, experiment and calculations suggest that as the temperature increases further, with higher energy, less constrained conformers of the methyl transition state become more significant. The implications of the experimental and theoretical results for the mechanisms of MF atmospheric oxidation and low-temperature combustion are discussed.

Kinetics of the Gas-Phase Reactions of syn- and anti-CH3CHOO Criegee Intermediate Conformers with SO2 as a Function of Temperature and Pressure.

Kinetics of reactions between SO2 and CH3CHOO Criegee intermediate conformers have been measured at temperatures between 242 and 353 K and pressures between 10 and 600 Torr using laser flash photolysis of CH3CHI2/O2/N2/SO2 gas mixtures coupled with time-resolved broadband UV absorption spectroscopy. The kinetics of syn-CH3CHOO + SO2 are pressure-dependent and exhibit a negative temperature dependence, with the observed pressure dependence reconciling apparent discrepancies between previous measurements performed at ∼298 K. Results indicate a rate coefficient of (4.80 ± 0.46) × 10-11 cm3 s-1 for the reaction of syn-CH3CHOO with SO2 at 298 K and 760 Torr. In contrast to the behavior of the syn-conformer, the kinetics of anti-CH3CHOO + SO2 display no significant dependence on temperature or pressure over the ranges investigated, with a mean rate coefficient of (1.18 ± 0.21) × 10-10 cm3 s-1 over all conditions studied in this work. Results indicate that the reaction of syn-CH3CHOO with SO2 competes with unimolecular decomposition and reaction with water vapor in areas with high SO2 concentration and low humidity, particularly at lower temperatures.

Solving the OH + Glyoxal Problem: A Complete Theoretical Description of Post-Transition-State Energy Deposition in Activated Systems.

Activated chemistry in coupled reaction systems has broadened our understanding of the chemical kinetics. In the case of intermediates formed in gas phase abstraction reactions (e.g., OH + HC(O)C(O)H (glyoxal) →HC(O)CO + H2O), it is particularly crucial to understand how the reaction energy is partitioned between product species as this determines the propensity for a given product to undergo "prompt" dissociation (e.g., HC(O)CO → HCO + CO) before the excess reaction energy is removed. An example of such an activated system is the OH + glyoxal + O2 coupled reaction system. In this work, we develop a molecular dynamics pipeline, which, combined with a master equation analysis, accurately models previous experimental measurements. This new work resolves previous complexities and discrepancies from earlier master equation modeling for this reaction system. The detailed molecular dynamics approach employed here is a powerful new tool for modeling challenging activated reaction systems.

New Measurements and Calculations on the Kinetics of an Old Reaction: OH + HO2 → H2O + O2.

Literature rate coefficients for the prototypical radical-radical reaction at 298 K vary by close to an order of magnitude; such variations challenge our understanding of fundamental reaction kinetics. We have studied the title reaction at room temperature via the use of laser flash photolysis to generate OH and HO2 radicals, monitoring OH by laser-induced fluorescence using two different approaches, looking at the direct reaction and also the perturbation of the slow OH + H2O2 reaction with radical concentration, and over a wide range of pressures. Both approaches give a consistent measurement of k1,298K ∼1 × 10-11 cm3 molecule-1 s-1, at the lowest limit of previous determinations. We observe, experimentally, for the first time, a significant enhancement in the rate coefficient in the presence of water, k1,H2O, 298K = (2.17 ± 0.09) × 10-28 cm6 molecule-2 s-1, where the error is statistical at the 1σ level. This result is consistent with previous theoretical calculations, and the effect goes some way to explaining some, but not all, of the variation in previous determinations of k1,298K. Supporting master equation calculations, using calculated potential energy surfaces at the RCCSD(T)-F12b/CBS//RCCSD/aug-cc-pVTZ and UCCSD(T)/CBS//UCCSD/aug-cc-pVTZ levels, are in agreement with our experimental observations. However, realistic variations in barrier heights and transition state frequencies give a wide range of calculated rate coefficients showing that the current precision and accuracy of calculations are insufficient to resolve the experimental discrepancies. The lower value of k1,298K is consistent with experimental observations of the rate coefficient of the related reaction, Cl + HO2 → HCl + O2. The implications of these results in atmospheric models are discussed.

Studies on the Kinetics of the CH + H2 Reaction and Implications for the Reverse Reaction, 3CH2 + H.

The reaction of CH radicals with H2 has been studied by the use of laser flash photolysis, probing CH decays under pseudo-first-order conditions using laser-induced fluorescence (LIF) over the temperature range 298-748 K at pressures of ∼5-100 Torr. Careful data analysis was required to separate the CH LIF signal at ∼428 nm from broad background fluorescence, and this interference increased with temperature. We believe that this interference may have been the source of anomalous pressure behavior reported previously in the literature (Brownsword, R. A.; J. Chem. Phys. 1997, 106, 7662-7677). The rate coefficient k1 shows complex behavior: at low pressures, the main route for the CH3* formed from the insertion of CH into H2 is the formation of 3CH2 + H, and as the pressure is increased, CH3* is increasingly stabilized to CH3. The kinetic data on CH + H2 have been combined with experimental shock tube data on methyl decomposition and literature thermochemistry within a master equation program to precisely determine the rate coefficient of the reverse reaction, 3CH2 + H → CH + H2. The resulting parametrization is kCH2+H(T) = (1.69 ± 0.11) × 10-10 × (T/298 K)(0.05±0.010) cm3 molecule-1 s-1, where the errors are 1σ.

Extensive field evidence for the release of HONO from the photolysis of nitrate aerosols.

Particulate nitrate ([Formula: see text]) has long been considered a permanent sink for NOx (NO and NO2), removing a gaseous pollutant that is central to air quality and that influences the global self-cleansing capacity of the atmosphere. Evidence is emerging that photolysis of [Formula: see text] can recycle HONO and NOx back to the gas phase with potentially important implications for tropospheric ozone and OH budgets; however, there are substantial discrepancies in "renoxification" photolysis rate constants. Using aircraft and ground-based HONO observations in the remote Atlantic troposphere, we show evidence for renoxification occurring on mixed marine aerosols with an efficiency that increases with relative humidity and decreases with the concentration of [Formula: see text], thus largely reconciling the very large discrepancies in renoxification photolysis rate constants found across multiple laboratory and field studies. Active release of HONO from aerosol has important implications for atmospheric oxidants such as OH and O3 in both polluted and clean environments.

Kinetics and Product Branching Ratio Study of the CH3O2 Self-Reaction in the Highly Instrumented Reactor for Atmospheric Chemistry.

The fluorescence assay by gas expansion (FAGE) method for the measurement of the methyl peroxy radical (CH3O2) using the conversion of CH3O2 into methoxy radicals (CH3O) by excess NO, followed by the detection of CH3O, has been used to study the kinetics of the self-reaction of CH3O2. Fourier transform infrared (FTIR) spectroscopy has been employed to determine the products methanol and formaldehyde of the self-reaction. The kinetics and product studies were performed in the Highly Instrumented Reactor for Atmospheric Chemistry (HIRAC) in the temperature range 268-344 K at 1000 mbar of air. The product measurements were used to determine the branching ratio of the reaction channel forming methoxy radicals, rCH3O. A value of 0.34 ± 0.05 (errors at 2σ level) was determined for rCH3O at 295 K. The temperature dependence of rCH3O can be parametrized as rCH3O = 1/{1 + [exp(600 ± 85)/T]/(3.9 ± 1.1)}. An overall rate coefficient of the self-reaction of (2.0 ± 0.9) × 10-13 cm3 molecule-1 s-1 at 295 K was obtained by the kinetic analysis of the observed second-order decays of CH3O2. The temperature dependence of the overall rate coefficient can be characterized by koverall = (9.1 ± 5.3) × 10-14 × exp((252 ± 174)/T) cm3 molecule-1 s-1. The found values of koverall in the range 268-344 K are ∼40% lower than the values calculated using the recommendations of the Jet Propulsion Laboratory and IUPAC, which are based on the previous studies, all of them utilizing time-resolved UV-absorption spectroscopy to monitor CH3O2. A modeling study using a complex chemical mechanism to describe the reaction system showed that unaccounted secondary chemistry involving Cl species increased the values of koverall in the previous studies using flash photolysis to initiate the chemistry. The overestimation of the koverall values by the kinetic studies using molecular modulation to generate CH3O2 can be rationalized by a combination of underestimated optical absorbance of CH3O2 and unaccounted CH3O2 losses to the walls of the reaction cells employed.

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.

New Approach to the Detection of Short-Lived Radical Intermediates.

We report a new general method for trapping short-lived radicals, based on a homolytic substitution reaction SH2'. This departure from conventional radical trapping by addition or radical-radical cross-coupling results in high sensitivity, detailed structural information, and general applicability of the new approach. The radical traps in this method are terminal alkenes possessing a nitroxide leaving group (e.g., allyl-TEMPO derivatives). The trapping process thus yields stable products which can be stored and subsequently analyzed by mass spectrometry (MS) supported by well-established techniques such as isotope exchange, tandem MS, and high-performance liquid chromatography-MS. The new method was applied to a range of model radical reactions in both liquid and gas phases including a photoredox-catalyzed thiol-ene reaction and alkene ozonolysis. An unprecedented range of radical intermediates was observed in complex reaction mixtures, offering new mechanistic insights. Gas-phase radicals can be detected at concentrations relevant to atmospheric chemistry.

Direct Measurements of Isoprene Autoxidation: Pinpointing Atmospheric Oxidation in Tropical Forests.

2-Methyl-1,3-butadiene (isoprene), released from biogenic sources, accounts for approximately a third of hydrocarbon emissions and is mainly removed by hydroxyl radicals, OH, the primary initiator of atmospheric oxidation. In situ measurements in clean tropical forests (high isoprene and low NO x ) have measured OH concentrations up to an order of magnitude higher than model predictions, which impacts our understanding of global oxidation. In this study, direct, laser flash photolysis, laser-induced fluorescence measurements at elevated temperatures have observed OH recycling in the presence of isoprene and oxygen under conditions where interference from secondary or heterogeneous chemistry is minimal. Our results provide the first direct, time-resolved, experimental validation of the theory-based Leuven Isoprene Mechanism (LIM1), based on isomerization of isoprene-RO2 radicals and OH regeneration, that partially accounts for model:measurement divergence in OH. While our data can be fit with only minor alterations in important LIM1 parameters, and the overall rate of product formation is similar to LIM1, there are differences with the recent experimental study by Teng et al. J. Am. Chem. Soc. 2017, 139, 5367-5377. In addition, our study indicates that the dihydroperoxide products are significantly enhanced over previous estimates. Dihydroperoxides are chemical and photochemical sources of OH, and the implications of enhanced hydroperoxide formation on the agreement between models and observations in tropical forests are examined.

OH Kinetics with a Range of Nitrogen-Containing Compounds: N-Methylformamide, t-Butylamine, and N-Methyl-propane Diamine.

Emissions of amines and amides to the atmosphere are significant from both anthropogenic and natural sources, and amides can be formed as secondary pollutants. Relatively little kinetic data exist on overall rate coefficients with OH, the most important tropospheric oxidant, and even less on site-specific data which control the product distribution. Structure-activity relationships (SARs) can be used to estimate both quantities. Rate coefficients for the reaction of OH with t-butylamine (k1), N-methyl-1,3-propanediamine (k2), and N-methylformamide (k3) have been measured using laser flash photolysis coupled with laser-induced fluorescence. Proton-transfer-reaction mass spectrometry (PTR-MS) has been used to ensure the reliable introduction of these low-vapor pressure N-containing compounds and to give qualitative information on products. Supporting ab initio calculations are presented for the t-butylamine system. The following rate coefficients have been determined: k1,298K= (1.66 ± 0.20) × 10-11 cm3 molecule-1 s-1, k(T)1 = 1.65 × 10-11 (T/300)-0.69 cm3 molecule-1 s-1, k2,293K = (7.09 ± 0.22) × 10-11 cm3 molecule-1 s-1, and k3,298K = (1.03 ± 0.23) × 10-11 cm3 molecule-1 s-1. For OH + t-butylamine, ab initio calculations predict that the fraction of N-H abstraction is 0.87. The dominance of this channel was qualitatively confirmed using end-product analysis. The reaction of OH with N-methyl-1,3-propanediamine also had a negative temperature dependence, but the reduction in the rate coefficient was complicated by reagent loss. The measured rate coefficient for reaction 3 is in good agreement with a recent relative rate study. The results of this work and the literature data are compared with the recent SAR estimates for the reaction of OH with reduced nitrogen compounds. Although the SARs reproduce the overall rate coefficients for reactions, site-specific agreement with this work and other literature studies is less strong.

Global Master Equation Analysis of Rate Data for the Reaction C2H4 + H ⇄ C2H5: ΔfH0⊖C2H5.

While forward and reverse rate constants are frequently used to determine enthalpies of reaction and formation, this process is more difficult for pressure-dependent association/dissociation reactions, especially since the forward and reverse reactions are usually studied at very different temperatures. The problems can be overcome by using a data-fitting procedure based on a master equation model. This approach has been applied to existing experimental pressure-dependent forward and reverse rate coefficients for the reaction C2H4 + H ⇄ C2H5 (k1, k-1) using the MESMER code to determine ΔfH0⊖C2H5 from the enthalpy of the reaction. New measurements of k1, k-1 were included in analysis. They are based on laser flash photolysis with direct observation of H atom time profiles by vacuum ultraviolet laser-induced fluorescence under conditions where the approach to equilibrium could be observed. Measurements were made over the temperature range 798-828 K and with [He] from 2.33 to 7.21 × 1018 molecule cm-3. These data were then combined with a wide range of existing experimental data with helium as the bath gas (112 measurements of k1 and k-1, covering the temperature range 285-1094 K, and [He] = 7.1 × 1015-1.9 × 1019 molecule cm-3) and fitted using the master equation solver MESMER. The required vibrational frequencies and rotational constants of the system were obtained from ab initio calculations, and the activation threshold for association (ΔEthresh), enthalpy of reaction (ΔrH0⊖), imaginary frequency (υimag), and helium energy-transfer parameters (⟨ΔE⟩d,298(T/298)n) were optimized. The resulting parameters (errors are 2σ) are ΔEthresh = 11.43 ± 0.34 kJ mol-1, ΔrH0⊖ = -145.34 ± 0.60 kJ mol-1, υimag = 730 ± 130 cm-1, ⟨ΔE⟩d,298 = 54.2 ± 7.6 cm-1, and n = 1.17 ± 0.12. A value of ΔfH298.15⊖(C2H5) = 120.49 ± 0.57 kJ mol-1 is obtained by combining ΔrH0⊖ with standard enthalpies of formation for H and C2H4 and making the appropriate temperature corrections. The dependence of these parameters on how the internal rotor and CH2 inversion modes are treated has been explored. The experimental data for other bath gases have been analyzed, and data sets compatible with the potential energy surface parameters determined above have been identified. The parameters are virtually identical but with slightly smaller error limits. Parameterization of k1, k-1 using the Troe formalization has been used to investigate competition between ethyl decomposition and reaction with oxygen under combustion conditions.

Kinetics of the gas phase reaction of the Criegee intermediate CH2OO with SO2 as a function of temperature.

The kinetics of the gas phase reaction of the Criegee intermediate CH2OO with SO2 have been studied as a function of temperature in the range 223-344 K at 85 Torr using flash photolysis of CH2I2/O2/SO2/N2 mixtures at 248 nm coupled to time-resolved broadband UV absorption spectroscopy. Measurements were performed under pseudo-first-order conditions with respect to SO2, revealing a negative temperature dependence. Analysis of experimental results using the Master Equation Solver for Multi-Energy well Reactions (MESMER) indicates that the observed temperature dependence, combined with the reported lack of a pressure dependence in the range 1.5-760 Torr, can be described by a reaction mechanism consisting of the formation of a pre-reaction complex leading to a cyclic secondary ozonide which subsequently decomposes to produce HCHO + SO3. The temperature dependence can be characterised by kCH2OO+SO2 = (3.72 ± 0.13) × 10-11 (T/298)(-2.05±0.38) cm3 molecule-1 s-1. The observed negative temperature dependence for the title reaction in conjunction with the decrease in water dimer (the main competitor for the Criegee intermediate) concentration at lower temperatures means that Criegee intermediate chemistry can play an enhanced role in SO2 oxidation in the atmosphere at lower temperatures.

Kinetics of the Reactions of Hydroxyl Radicals with Furan and Its Alkylated Derivatives 2-Methyl Furan and 2,5-Dimethyl Furan.

Furans are promising second generation biofuels with comparable energy densities to conventional fossil fuels. Combustion of furans is initiated and controlled to a large part by reactions with OH radicals, the kinetics of which are critical to understand the processes occurring under conditions relevant to low-temperature combustion. The reactions of OH radicals with furan (OH + F, R1), 2-methyl furan (OH + 2-MF, R2), and 2,5-dimethyl furan (OH + 2,5-DMF, R3) have been studied in this work over the temperature range 294-668 K at pressures between 5 mbar and 10 bar using laser flash photolysis coupled with laser-induced fluorescence (LIF) spectroscopy to generate and monitor OH radicals under pseudo-first-order conditions. Measurements at p ≤ 200 mbar were made in N2, using H2O2 or (CH3)3COOH radical precursors, while those at p ≥ 2 bar were made in He, using HNO3 as the radical precursor. The kinetics of reactions R1-R3 were observed to display a negative dependence on temperature over the range investigated, indicating the dominance of addition reactions under such conditions, with no significant dependence on pressure observed. Master equation calculations are in good agreement with the observed kinetics, and a combined parametrization of addition channels and abstraction channels for R1-R3 is provided on the basis of this work and previous shock tube measurements at higher temperatures. This work significantly extends the temperature range previously investigated for R1 and represents the first temperature-dependent measurements of R2 and R3 at temperatures relevant for atmospheric chemistry and low-temperature combustion.

CH2OO Criegee intermediate UV absorption cross-sections and kinetics of CH2OO + CH2OO and CH2OO + I as a function of pressure.

The UV absorption cross-sections of the Criegee intermediate CH2OO, and kinetics of the CH2OO self-reaction and the reaction of CH2OO with I are reported as a function of pressure at 298 K. Measurements were made using pulsed laser flash photolysis of CH2I2/O2/N2 gas mixtures coupled with time-resolved broadband UV absorption spectroscopy at pressures between 6 and 300 Torr. Results give a peak absorption cross-section of (1.37 ± 0.29) × 10-17 cm2 at ∼340 nm and a rate coefficient for the CH2OO self-reaction of (8.0 ± 1.1) × 10-11 cm3 s-1, with no significant pressure dependence of the absorption cross-sections or the self-reaction kinetics over the range investigated. The rate coefficient for the reaction between CH2OO and I demonstrates pressure dependence over the range investigated, with a Lindemann fit giving k0 = (4.4 ± 1.0) × 10-29 cm6 s-1 and k∞ = (6.7 ± 0.6) × 10-11 cm3 s-1. The origins of IO in the system have been investigated, the implications of which are discussed.

Temperature and Pressure Dependent Kinetics of QOOH Decomposition and Reaction with O2: Experimental and Theoretical Investigations of QOOH Radicals Derived from Cl + (CH3)3COOH.

QOOH radicals are key species in autoignition, produced by internal isomerizations of RO2 radicals, and are central to chain branching reactions in low-temperature combustion. The kinetics of QOOH radical decomposition and reaction with O2 have been determined as a function of temperature and pressure, using observations of OH radical production and decay following H-atom abstraction from t-butyl hydroperoxide ((CH3)3COOH) by Cl atoms to produce QOOH (CH2(CH3)2COOH) radicals. The kinetics of QOOH decomposition have been investigated as a function of temperature (251-298 K) and pressure (10-350 Torr) in helium and nitrogen bath gases, and those of the reaction between QOOH and O2 have been investigated as a function of temperature (251-304 K) and pressure (10-100 Torr) in He and N2. Decomposition of the QOOH radical was observed to display temperature and pressure dependence, with a barrier height for the decomposition of (44.7 ± 4.0) kJ mol-1 determined by master equation fitting to the experimental data. The rate coefficient for the reaction between QOOH and O2 was determined to be (5.6 ± 1.7) × 10-13 cm3 s-1 at 298 K, with no significant dependence on pressure, and can be described by the Arrhenius parameters A = (7.3 ± 6.8) × 10-14 cm3 s-1 and Ea = -(5.4 ± 2.1) kJ mol-1 in the temperature range of 251-304 K. This work represents the first measurements of any QOOH radical kinetics as a function of temperature and pressure.