Direct Trace Fitting of Experimental Data Using the Master Equation: Testing Theory and Experiments on the OH + C2H4 Reaction

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

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

OH(2Π) + C2H4 Reaction: A Combined Crossed Molecular Beam and Theoretical Study

The reaction between the ground-state hydroxyl radical, OH(²Π), and ethylene, C₂H₄, has been investigated under single-collision conditions by the crossed molecular beam scattering technique with mass-spectrometric detection and time-of-flight analysis at the collision energy of 50.4 kJ/mol. Electronic structure calculations of the underlying potential energy surface (PES) and statistical Rice-Ramsperger-Kassel-Marcus (RRKM) calculations of product branching fractions on the derived PES for the addition pathway have been performed. The theoretical results indicate a temperature-dependent competition between the anti-/syn-CH₂CHOH (vinyl alcohol) + H, CH₃CHO (acetaldehyde) + H, and H₂CO (formaldehyde) + CH₃ product channels. The yield of the H-abstraction channel could not be quantified with the employed methods. The RRKM results predict that under our experimental conditions, the anti- and syn-CH₂CHOH + H product channels account for 38% (in similar amounts) of the addition mechanism yield, the H₂CO + CH₃ channel for ∼58%, while the CH₃CHO + H channel is formed in negligible amount (<4%). The implications for combustion and astrochemical environments 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 H₂ 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 k₁ shows complex behavior: at low pressures, the main route for the CH₃* formed from the insertion of CH into H₂ is the formation of ³CH₂ + H, and as the pressure is increased, CH₃* is increasingly stabilized to CH₃. The kinetic data on CH + H₂ 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, ³CH₂ + H → CH + H₂. The resulting parametrization is k_CH2+H(T) = (1.69 ± 0.11) × 10^-10 × (T/298 K)^{(0.05±0.010)} cm³ molecule^-1 s^-1, where the errors are 1σ.

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-RO₂ 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.

Nitrene formation is the first step of the thermal and photochemical decomposition reactions of organic azides

In this work, the decomposition of a prototypical azide, isopropyl azide, both in the ground and excited states, has been investigated through the use of multiconfigurational CASSCF and MS-CASPT2 electronic structure approaches. Particular emphasis has been placed on the thermal reaction starting at the S₀ ground state surface. It has been found that the azide thermally decomposes via a stepwise mechanism, whose rate-determining step is the formation of isopropyl nitrene, which is, in turn, the first step of the global mechanism. After that, the nitrene isomerizes to the corresponding imine derivative. Two routes are possible for such a decomposition: (i) a spin-allowed path involving a transition state; and (ii) a spin-forbidden one via a S₀/T₀ intersystem crossing. Both intermediates have been determined and characterised. Their associated relative energies have been found to be quite similar, 45.75 and 45.52 kcal mol^-1, respectively. To complete this study, the kinetics of the singlet and triplet channels are modeled with the MESMER (Master Equation Solver for Multi-Energy Well Reactions) code by applying the RRKM and Landau-Zener (with WKB tunnelling correction) theories, respectively. It is found that the canonical rate-coefficients of the singlet path are 2-orders of magnitude higher than the rate-coefficients of the forbidden reaction. In addition, the concerted mechanism has been investigated that would lead to the formation of the imine derivative and nitrogen extrusion in the first step of the decomposition. After a careful analysis of CASSCF calculations with different active spaces and their comparison with single electronic configuration methods (MP2 and B3LYP), the concerted mechanism is discarded.

An experimental and computational study of the reaction between 2-methylallyl radicals and oxygen molecules: optimizing master equation parameters with trace fitting

We have investigated the reaction between 2-methylallyl radicals and oxygen molecules with experimental and computational methods. Kinetic experiments were conducted in a tubular laminar flow reactor using laser photolysis for radical production and photoionization mass spectrometry for detection. The reaction was investigated as a function of temperature (203-730 K) and pressure (0.2-9 torr) in helium and nitrogen bath gases. At low temperatures (T < 410 K), the reaction proceeds by a barrierless reaction to form 2-methylallylperoxyl. Equilibration of the peroxyl adduct and the reactants was observed between 350-410 K. Measurements were extended to even higher temperatures, up to 730 K, but no reaction could be observed. Master equation simulations of the reaction system were performed with the MESMER program. Kinetic parameters in the master equation model were optimized by direct fitting to time-resolved experimental 2-methylallyl traces. Trace fitting is a recently implemented novel feature in MESMER. The trace approach was compared with the more traditional approach where one uses experimental rate coefficients for parameter optimization. The optimized parameters yielded by the two approaches are very similar and do an excellent job at reproducing the experimental data. The optimized master equation model was then used to simulate the reaction under study over a wide temperature and pressure range, from 200 K and 0.01 bar to 1500 K and 100 bar. The simulations predict a small phenomenological rate coefficient under autoignition conditions; about 1 × 10^-18 cm³ s^-1 at 400 K and 5 × 10^-16 cm³ s^-1 at 1000 K. We provide modified Arrhenius expressions in PLOG format for the most important product channels to facilitate the use of our results in combustion models.

An Experimental and Computational Study of the Reaction between Pent-3-en-2-yl Radicals and Oxygen Molecules: Switching from Pure Stabilisation to Pure Decomposition with Increasing Temperature

We have used laser-photolysis – photoionization mass spectrometry, quantum chemical calculations, and master equation simulations to investigate the kinetics of the reaction between (E/Z)-pent-3-en-2-yl (CH3CHCHCHCH3), a resonance-stabilised hydrocarbon radical, and...

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 C₂H₄ + H ⇄ C₂H₅ (k₁, k_-1) using the MESMER code to determine Δ_fH₀^⊖C₂H₅ from the enthalpy of the reaction. New measurements of k₁, 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 × 10¹⁸ 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 k₁ and k_-1, covering the temperature range 285-1094 K, and [He] = 7.1 × 10¹⁵-1.9 × 10¹⁹ 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 (ΔE_thresh), enthalpy of reaction (Δ_rH₀^⊖), imaginary frequency (υ_imag), and helium energy-transfer parameters (⟨ΔE⟩_d,298(T/298)ⁿ) were optimized. The resulting parameters (errors are 2σ) are ΔE_thresh = 11.43 ± 0.34 kJ mol^-1, Δ_rH₀^⊖ = -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 Δ_fH_298.15^⊖(C₂H₅) = 120.49 ± 0.57 kJ mol^-1 is obtained by combining Δ_rH₀^⊖ with standard enthalpies of formation for H and C₂H₄ and making the appropriate temperature corrections. The dependence of these parameters on how the internal rotor and CH₂ 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 k₁, k_-1 using the Troe formalization has been used to investigate competition between ethyl decomposition and reaction with oxygen under combustion conditions.

Dissociative Photodetachment Dynamics of the OH-(C2H4) Anion Complex

Photoelectron-photofragment coincidence (PPC) measurements on OH^-(C₂H₄) anions at a photon energy of 3.20 eV revealed stable and dissociative photodetachment product channels, OH-C₂H₄ + e^- and OH + C₂H₄ + e^-, respectively. The main product channel observed was dissociation to the reactants (>67%), OH + C₂H₄ (v = 0, 1, 2) + e^-, where vibrational excitation in the C-H stretching modes of the C₂H₄ photofragments corresponds to a minor channel. The low kinetic energy release (KER) of the dissociating fragments is consistent with weak repulsion between the OH + C₂H₄ reactants near the transition state as well as the partitioning of energy into rotation of the dissociation products. An impulsive model was used to account for rotational energy partitioning in the dissociative photodetachment (DPD) process and showed good agreement with the experimental results. The low KER of the dissociating fragments and the similarities in the photoelectron spectra between stable and dissociative events support a mechanism involving the van der Waals complex formed upon photodetachment of OH^-(C₂H₄) as an intermediate in the dominant OH + C₂H₄ + e^- dissociative channel.

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 N₂, using H₂O₂ or (CH₃)₃COOH radical precursors, while those at p ≥ 2 bar were made in He, using HNO₃ 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.