Reaction of ethylene with hydroxyl radicals: a theoretical study

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.

Theoretical investigation on isomerization and decomposition reactions of pentanol radicals-part II: linear pentanol isomers

High-level ab initio calculations are conducted for studying the kinetics of three linear pentanol radicals generated through H-atom abstraction reactions. The species involved are optimized using the M06-2X/6-311++G(d,p) level of theory, while a relaxed scan at the M06-2X/6-31g level of theory with 10° increments is used for the hindrance potential for low-frequency torsional modes. Single-point energies for all stationary points are obtained through the QCISD(T) and MP2 methods in combination with cc-pVDZ, cc-pVTZ, and cc-pVQZ basis sets, which can be extrapolated to the complete basis set (CBS) limit. The rate constants and branching ratios for isomerization and decomposition reactions are computed over a temperature range of 250-2000 K and a pressure range of 0.01-100 atm. Isomerization reactions are dominant at low temperatures, while decomposition reactions are more dominant at high temperatures. The branching ratio of the isomerization reaction exhibits a slight decrease with increasing pressure, while the trend for decomposition reactions depends on the type of the breaking bond. Based on the calculations for five branched pentanol radicals in part I, kinetics of linear and branched pentanol radicals are compared in this work and the results reveal that, for the same kind of β-scission reaction at similar positions of linear and branched pentanol radicals, the rate constants of branched ones are faster than those of linear ones at low temperatures. The hydroxyl group adjacent to the breaking bond can increase the β-scission reaction rate constants, while the effect can be ignored when the hydroxyl group is not adjacent to the breaking bond. Moreover, compared to when the hydroxyl group is located in the middle of the carbon chain, its positioning at the chain's end yields a more noticeable impact on the products and rate constants of C-O bond and O-H bond β-scission reactions. Besides, when incorporating calculated rate constants into the CRECK model, the updated mechanism shows a better performance for ignition delay times of 1-pentanol in the NTC range but exhibits lower reactivity at higher temperatures. The simulation of speciation profiles also shows better agreement with the experimental data obtained using a flow reactor.

Quantum chemistry and kinetics of hydrogen sulphide oxidation

A fundamental understanding of the acid gas (H2S and CO2) chemistry is key to efficiently implement the desulphurisation process and even the production of clean fuels such as hydrogen or syngas. In this work, we developed a new kinetic model for the pyrolysis and oxidation of hydrogen sulphide by merging two previously reported models with the goal of covering a wider range of conditions and including the effect of carbon dioxide. The resulting model, which consists of 75 species and 514 reactions, was used to conduct rate of production and sensitivity analysis in plug flow reactor simulations, and the results were used to determine the most prominent reactions in which hydrogen sulphide, molecular hydrogen, and sulphur monoxide are involved. The resulting list of important reactions was screened and the kinetics of three of them, i.e., SO2 + S2 → S2O + SO, S2O + S2 → S3 + SO, and SO + SH → S2 + OH, was found to warrant further investigation. With the goal of improving the accurancy of our new kinetic model, we carried out a robust quantum chemistry and Rice-Ramsperger-Kassel-Marcus master equation study to obtain, for the first time, the forward and reverse rate constants for those three reactions at temperatures and pressures of interest for combustion and atmospheric chemistry. This work is the first step of a kinetic study that is aimed at improving the understanding of the chemistry of the pyrolysis and oxidation of H2S, highlighting the importance of sulphur-sulphur interactions and providing a fundamental basis for future kinetic models of H2S not only in the field of combustion, but also in atmospheric chemistry.

Infrared spectroscopy of the α-hydroxyethyl radical isolated in cryogenic solid media

The α-hydroxyethyl radical (CH3·CHOH, 2A) is a key intermediate in ethanol biochemistry, combustion, atmospheric chemistry, radiation chemistry, and astrochemistry. Experimental data on the vibrational spectrum of this radical are crucially important for reliable detection and understanding of the chemical dynamics of this species. This study represents the first detailed experimental report on the infrared absorption bands of the α-hydroxyethyl radical complemented by ab initio computations. The radical was generated in solid para-H2 and Xe matrices via the reactions of hydrogen atoms with matrix-isolated ethanol molecules and radiolysis of isolated ethanol molecules with x rays. The absorption bands with maxima at 3654.6, 3052.1, 1425.7, 1247.9, 1195.6 (1177.4), and 1048.4 cm-1, observed in para-H2 matrices appearing upon the H· atom reaction, were attributed to the OHstr, α-CHstr, CCstr, COstr + CCObend, COstr, and CCstr + CCObend vibrational modes of the CH3·CHOH radical, respectively. The absorption bands with the positions slightly red-shifted from those observed in para-H2 were detected in both the irradiated and post-irradiation annealed Xe matrices containing C2H5OH. The results of the experiments with the isotopically substituted ethanol molecules (CH3CD2OH and CD3CD2OH) and the quantum-chemical computations at the UCCSD(T)/L2a_3 level support the assignment. The photolysis with ultraviolet light (240-300 nm) results in the decay of the α-hydroxyethyl radical, yielding acetaldehyde and its isomer, vinyl alcohol. A comparison of the experimental and theoretical results suggests that the radical adopts the thermodynamically more stable anti-conformation in both matrices.

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.

Mai T, Assali M, Nguyen TTD, Nguyen H, Szőri M, Huynh LK, Fittschen C, Farooq A. Reaction Kinetics of 1,4-Cyclohexadiene with OH radicals : An Experimental and Theoretical Study

This work presents OH-initiated oxidation kinetics of 1,4-cyclochexadiene (1,4-CHD). Temperature dependence of the reaction was investigated by utilizing laser flash photolysis flow reactor and laser-induced fluorescence (LPFR/LIF) technique over the...

The Effect of Hydrogen Bonding on the Reactivity of OH Radicals with Prenol and Isoprenol: A Shock Tube and Multi-Structural Torsional Variational Transition State Theory Study

The presence of two functional groups (OH and double bond) in C5 methyl-substituted enols (i.e., isopentenols), such as 3-methyl-2-buten-1-ol (prenol) and 3-methyl-3-buten-1-ol (isoprenol), makes them excellent biofuel candidates as fuel...

Effect of a single water molecule on ˙CH2OH + 3O2 reaction under atmospheric and combustion conditions

The hydroxymethyl (˙CH₂OH) radical is an important intermediate species in both atmosphere and combustion reaction systems. The rate coefficients for ˙CH₂OH + ³O₂ and (˙CH₂OH + ³O₂ (+H₂O)) reactions were calculated using the Rice-Ramsperger-Kassel-Marcus (RRKM)/master equation (ME) simulation and canonical variational transition state theory (CVT) between the temperature range of 200 to 1500 K based on the potential energy surface constructed using CCSD(T)//ωB97XD/6-311++G(3df,3pd). The results show that ˙CH₂OH + ³O₂ leads to the formation of CH₂O and HO₂ at temperatures below 800 K, and goes back to reactants at high temperature (>1000 K). When a water molecule is added to the reaction, the formation of CH₂O and HO₂ is favored at all temperatures. The calculated rate coefficient for the ˙CH₂OH + ³O₂ (2.8 × 10^-11 cm³ molecule^-1 s^-1 at 298 K) is in good agreement with the previous experimental values (∼1 × 10^-11 cm³ molecule^-1 s^-1 at 298 K). The rate coefficients for the water-assisted reaction (2.4 × 10^-16 cm³ molecule^-1 s^-1 at 1000 K) is at least 3-4 orders of magnitude smaller than the water-free reaction (6.2 × 10^-12 cm³ molecule^-1 s^-1 at 1000 K). This result is consistent with the similar types of reaction system. Our calculations also predict that the effect of a single water molecule favors the formation of CH₂O in the combustion condition. However, the water-free reaction favors the formation of CH₂O in the atmospheric condition. The current study helps to understand how a single water molecule changes the reaction mechanism and chemical kinetic behaviour under atmospheric and combustion conditions.

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.

Relative energetics of CH3CH2O, CH3CHOH, and CH2CH2OH radical products from ethanol dehydrogenation

This study has examined the relative energetics of nine stationary points associated with the three different radical isomers generated by removing a H atom from ethanol at the O atom (ethoxy, CH₃CH₂O), the α C atom (CH₃CHOH), and the β C atom (CH₂CH₂OH). For the first time, CCSD(T) geometry optimizations and harmonic vibrational frequency computations with the cc-pVTZ and aug-cc-pVTZ basis sets have been carried out to characterize two unique minima for each isomer along with three transition state structures with C_s symmetry. Explicitly correlated CCSD(T) computations were also performed to estimate the relative energetics of these nine stationary points near the complete basis set limit. These benchmark results were used to assess the performance of various density functional theory (DFT) and wave function theory methods, and they will help guide method selection for future studies of alcohols and their radicals. The structures generated by abstracting H from the α C atom have significantly lower electronic energies (by at least 7 kcal mol^-1) than the CH₃CH₂O and CH₂CH₂OH radicals. Although previously reported as a minimum on the ground-state surface, the ²A″ C_s structure of the ethoxy radical was found to be a transition state in this study with MP2, CCSD(T), and a number of DFT methods. An implicit solvation model used in conjunction with DFT and MP2 methods did not qualitatively change the relative energies of the isomers, but the results suggest that the local minima for the CH₃CHOH and CH₂CH₂OH radicals could become more energetically competitive in condensed phase environments, such as liquid water and ethanol.

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.

Collision Efficiency Parameter Influence on Pressure-Dependent Rate Constant Calculations Using the SS-QRRK Theory

The system-specific quantum Rice-Ramsperger-Kassel (SS-QRRK) theory (J. Am. Chem. Soc. 2016, 138, 2690) is suitable to determine rate constants below the high-pressure limit. Its current implementation allows incorporating variational effects, multidimensional tunneling, and multistructural torsional anharmonicity in rate constant calculations. Master equation solvers offer a more rigorous approach to compute pressure-dependent rate constants, but several implementations available in the literature do not incorporate the aforementioned effects. However, the SS-QRRK theory coupled with a formulation of the modified strong collision model underestimates the value of unimolecular pressure-dependent rate constants in the high-temperature regime for reactions involving large molecules. This underestimation is a consequence of the definition for collision efficiency, which is part of the energy transfer model. Selection of the energy transfer model and its parameters constitutes a common issue in pressure-dependent calculations. To overcome this underestimation problem, we evaluated and implemented in a bespoke Python code two alternative definitions for the collision efficiency using the SS-QRRK theory and tested their performance by comparing the pressure-dependent rate constants with the Rice-Ramsperger-Kassel-Marcus/Master Equation (RRKM/ME) results. The modeled systems were the tautomerization of propen-2-ol and the decomposition of 1-propyl, 1-butyl, and 1-pentyl radicals. One of the tested definitions, which Dean et al. explicitly derived (Z. Phys. Chem. 2000, 214, 1533), corrected the underestimation of the pressure-dependent rate constants and, in addition, qualitatively reproduced the trend of RRKM/ME data. Therefore, the used SS-QRRK theory with accurate definitions for the collision efficiency can yield results that are in agreement with those from more sophisticated methodologies such as RRKM/ME.

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

Laser flash photolysis coupled with laser-induced fluorescence observation of OH has been used to observe the equilibration of OH + C₂H₄ ↔ HOC₂H₄ over the temperature range 563-723 K and pressures of bath gas (N₂) from 58 to 250 Torr. The time-resolved OH traces have been directly and globally fitted with a master equation in order to extract Δ_RH₀⁰, the binding energy of the HOC₂H₄ adduct, with respect to reagents. The global approach allows the role that OH abstraction plays at higher temperatures to be identified. The resultant value ofΔ_RH₀⁰, 111.8 kJ mol^-1, is determined to be better than 2 kJ mol^-1 and is in agreement with our ab initio calculations (carried out at the CCSD(T)/CBS//M06-2X/aug-cc-pVTZ level), 111.4 kJ mol^-1, and other state of the art calculations. Parameters for the abstraction channel are also in good agreement with previous experimental studies. To effect this analysis, the MESMER master equation code was extended to directly incorporate secondary chemistry: diffusional loss from the observation region and reaction with the photolytic precursor. These extensions, which, among other things, resolve issues related to the merging of chemically significant and internal energy relaxation eigenvalues, are presented.

Kinetics of the benzyl + HO2 and benzoxyl + OH barrierless association reactions: fate of the benzyl hydroperoxide adduct under combustion and atmospheric conditions

Radical-radical association reactions are challenging to address theoretically due to difficulties finding the bottleneck that variationally minimizes the reactive flux. For this purpose, the variable reaction coordinate (VRC) formulation of the variational transition state theory (VTST) represents an appropriate tool. In this work, we revisited the kinetics of two radical-radical association reactions of importance in combustion modelling and poly-aromatic hydrocarbon (PAH) chemistry by performing VRC calculations: benzyl + HO2 and benzoxyl + OH, both forming the adduct benzyl hydroperoxide. Our calculated rate constants are significantly lower than those previously reported based on VTST calculations, which results from a more efficient minimization of the reactive flux through the bottleneck achieved by the VRC formulation. Both reactions show different trends in the variation of their rate constants with temperature. We observed that if the pair of single occupied molecular orbitals (SOMOs) of the associating radicals show a similar nature, i.e. similar character, and thereby a small energy gap, a highly stabilized transition state structure is formed as the result of a very efficient SOMO-SOMO overlap, which may cancel out the free energy bottleneck for the formation of the adduct and result in large rate constants with a negative temperature dependence. This is the case of the benzoxyl and OH radical pair, whose SOMOs show O2p nature with an energy gap of 20.2 kcal mol-1. On the other hand, the benzyl and HO2 radical pair shows lower rate constants with a positive temperature dependence due to the larger difference between both SOMOs (a 28.9 kcal mol-1 energy gap) as a consequence of the contribution of the multiple resonance structures of the benzyl radical. The reverse dissociation rate constants were also calculated using multi-structural torsional anharmonicity partition functions, which were not included in previous work, and the results show a much slower dissociation of benzyl hydroperoxide. Our work may help to improve kinetic models of interest in combustion and PAH formation, as well as to gain further understanding of radical-radical association reactions, which are ubiquitous in different environments.

Prediction of Rate Coefficients for the H2CO + OH → HCO + H2O Reaction at Combustion, Atmospheric and Interstellar Medium Conditions

Despite the relevance of the H₂CO + OH → HCO + H₂O reaction for combustion, atmospheric, and interstellar medium conditions, a large discrepancy on energetic and kinetic data for this reaction is still observed in the previous literature. In this work, this hydrogen abstraction reaction has been investigated at the CCSD(T)/CBS level of theory, suggesting that both the prebarrier complex and saddle point are stabilized in relation to the reactants by 3.31 and 1.35 kcal mol^-1, respectively. Moreover, from the Gibbs free energy profile of the reaction coordinate, it has been verified that the formation of the prebarrier complex is endergonic, for temperatures above 550 K. Hence, for temperatures lower than 550 K, the reaction is described by a mechanism consisting of three elementary steps, while for higher temperatures it can be assumed to be an elementary reaction. Finally, the prediction of rate coefficients suggests that unified statistical rate theory best applies to the low temperature regime, while canonical variational rate coefficients better fit experimental data at the high temperature regime.

Ab initio kinetics of the C2H2 + NH2 reaction: a revisited study

This work provides a rigorous detailed kinetic study on the C₂H₂ + NH₂ reaction in a wide range of conditions (T = 250-2000 K & P = 1-76000 Torr). In particular, the composite method W1U was used to construct the potential energy surface on which the kinetic behaviors were characterized within the state-of-the-art master equation/Rice-Ramsperger-Kassel-Marcus (ME/RRKM) framework. Corrections of the hindered internal rotation (HIR) treatment and quantum tunneling effect were included. A clear reaction mechanism shift with respect to both temperature and pressure was revealed via detailed kinetic and species analyses. In particular, bimolecular products (i.e., CH₂[double bond, length as m-dash]C[double bond, length as m-dash]NH + H, CH[triple bond, length as m-dash]CNH₂ + H, CH₃CN + H, CH[triple bond, length as m-dash]C· + NH₃ in the decreasing mole fraction order) can be formed directly from the reactants at high temperature and/or low pressure while they can be produced indirectly via intermediates (e.g., ·CH[double bond, length as m-dash]CHNH₂(cis), ·CH[double bond, length as m-dash]CHNH₂(trans), CH₂[double bond, length as m-dash]C·NH₂,…) at low temperature and/or high pressure. The calculated rate constants are in good agreement with the literature data from ab initio calculations without any adjustment; thus, the proposed temperature- and pressure-dependent rate constants, together with the thermodynamic data of the species involved, can be confidently used for modeling NH₂-related systems under atmospheric and combustion conditions.

High-Temperature Unimolecular Decomposition of Diethyl Ether: Shock-Tube and Theory Studies

The unimolecular decomposition of diethyl ether (DEE; C₂H₅OC₂H₅) is considered to be initiated via a molecular elimination and a C-O and a C-C bond fission step: C₂H₅OC₂H₅ → C₂H₄ + C₂H₅OH (1), C₂H₅OC₂H₅ → C₂H₅ + C₂H₅O (2), and C₂H₅OC₂H₅ → CH₃ + C₂H₅OCH₂ (3). In this work, two shock-tube facilities were used to investigate these reactions via (a) time-resolved H-atom concentration measurements by H-ARAS (atomic resonance absorption spectrometry), (b) time-resolved DEE-concentration measurements by high repetition-rate time-of-flight mass spectrometry (HRR-TOF-MS), and (c) product-composition measurements via gas chromatography/MS (GC/MS) after quenching the test gas. The experiments were conducted at temperatures ranging from 1054 to 1505 K and at pressures between 1.2 and 2.5 bar. Initial DEE mole fractions between 0.4 and 9300 ppm were used to perform the kinetics experiments by H-ARAS (0.4 ppm), GC/MS (200-500 ppm), and HRR-TOF-MS (7850-9300 ppm). The rate constants, k_total (k_total = k₁ + k₂ + k₃) derived from the GC/MS and HRR-TOF-MS experiments agree well with each other and can be described by the Arrhenius expression, k_total(1054-1467 K; 1.3-2.5 bar) = 10^12.81±0.22 exp(-240.27 ± 5.11 kJ mol^-1/RT) s^-1. From the H-ARAS experiments, overall rate constants for the bond fission channels, k₂₊₃ = k₂ + k₃ have been extracted. The k₂₊₃ data can be well described by the Arrhenius equation, k₂₊₃(1299-1505 K; 1.3-2.5 bar) = 10^14.43±0.33 exp(-283.27 ± 8.78 kJ mol^-1/RT) s^-1. A master-equation analysis was performed using CCSD(T)/aug-cc-pvtz//B3LYP/aug-cc-pvtz and CASPT2/aug-cc-pvtz//B3LYP/aug-cc-pvtz molecular properties and energies for the three primary thermal decomposition processes in DEE. The derived experimental data is very well reproduced by the simulations with the mechanism of this work. With regard to the branching ratios between bond fissions and elimination channels, uncertainties remain.

Temperature and Pressure Dependent Rate Coefficients for the Reaction of Ketene with Hydroxyl Radical

The reaction of ketene with hydroxyl radical is drawing growing attention, for it is found to constitute an important step during the combustion of hydrocarbon and oxygenated hydrocarbon fuels, e.g., acetylene, propyne, allene, acetone, gasoline, diesel, jet fuels, and biofuels. We studied the potential energy surface (PES) of this reaction using B2PLYP-D3/cc-PVTZ for geometry optimization and composite methods based on CCSD(T)-F12/cc-PVTZ-F12 for energy calculations. From this PES, temperature- and pressure-dependent rate coefficients and branching ratios at 200-3000 K and 0.01-100 atm were derived using the RRKM/ME approach. The reaction is dominated by four product channels: (i) OH addition on the olefinic carbon of ketene to form CH₂OH + CO, which is the most dominant under all conditions; (ii) H abstraction producing HCCO + H₂O, which is favored at high temperatures; (iii) OH addition on the carbonyl carbon to form CH₃ + CO₂, which is favored at low pressures and high temperatures; and (iv) collisional stabilization of CH₂COOH, which is favored at high pressures and low temperatures. With increasing temperatures, the overall rate constant k_overall exhibit first negative but then positive temperature dependency, with its switching point (also the minimum point) at ∼400 K. Both product channel CH₂OH + CO and HCCO + H₂O are independent of pressure, whereas formation of CH₃ + CO₂ and collisional stabilization of CH₂COOH are highly pressure dependent. Fitted modified Arrhenius expressions of the calculated rate constants are provided for the purpose of combustion modeling.

Effects of water vapor on the reaction of CH2OO with NH3

A strong synergic effect of water and ammonia molecules may enhance the formation of H₂NCH₂OOH.

A quantum theory investigation on atmospheric oxidation mechanisms of acrylic acid by OH radical and its implication for atmospheric chemistry

The hydroxyl radical, as the most important oxidant, controls the removal of some volatile organic compounds (VOCs) in the atmosphere. In this work, the atmospheric oxidation processes of acrylic acid by OH radical have been investigated by density functional theory (DFT). The energetic routes of the reaction of CH₂CHCOOH with OH radical have been calculated accurately at the CCSD(T)/cc-pVTZ//M06-2X/6-311++G(d,p) level. It is implicated that the oxidation has five elementary reaction pathways mostly hinging on how hydroxyl radical approaches to the carbon skeleton of acrylic acid. The atmospheric degradation mechanisms of the CH₂CHCOOH by OH radical are the formation of reactive intermediates IM1 and IM2. Meanwhile, the further oxidation mechanisms of IM1 and IM2 by O₃ and NO are also investigated. The rate coefficients have been computed using tight transition state theory of the variflex code. The calculated rate coefficient is 2.3 × 10^-11 cm³ molecule^-1 s^-1 at standard pressure and 298 K, which is very close to the laboratory data (1.75 ± 0.47 × 10^-11 cm³ molecule^-1 s^-1). Moreover, the atmospheric lifetime of acrylic acid is about 6 h at 298 K and 1 atm, implying that the fast sinks of acrylic acid by hydroxyl radical.

First kinetic study of the atmospherically important reactions BrHg˙ + NO2 and BrHg˙ + HOO

We use computational chemistry to determine the rate constants and product yields for the reactions of BrHg˙ with the atmospherically abundant radicals NO₂ and HOO. The reactants, products, and well-defined transition states are characterized using CCSD(T) with large basis sets. The potential energy profiles for the barrierless addition of HOO and NO₂ to BrHg˙ are characterized using CASPT2 and RHF-CCSDT, and the rate constants are computed as a function of temperature and pressure using variational transition state theory and master equation simulations. The calculated rate constant for the addition of NO₂ to BrHg˙ is larger than that for the addition of HOO by a factor of up to two under atmospheric conditions. For the reaction of HOO with BrHg˙ the addition reaction entirely dominates competing HOO + BrHg˙ reaction channels. The addition of NO₂ to BrHg˙ initially produces both BrHgNO₂ and BrHgONO, but after a few seconds under atmospheric conditions the sole product is syn-BrHgONO. A previously unsuspected reaction channel for BrHg˙ + NO₂ competes with the addition to yield Hg + BrNO₂. This reaction reduces the mercury oxidation state in BrHg˙ from Hg(i) to Hg(0) and slows the atmospheric oxidation of Hg(0). While the rate constant for this reduction channel is not well-constrained by the present calculations, it may be as much as 18% as large as the oxidation channel under some atmospheric conditions. As no experimental kinetic or product yield data are available for the reactions studied here, this work will provide guidance for atmospheric modelers and experimental kineticists.

Atmospheric oxidation of halogenated aromatics: comparative analysis of reaction mechanisms and reaction kinetics

Atmospheric transport is the major route for global distribution of semi-volatile compounds such as halogenated aromatics as well as their major exposure route for humans. Their major atmospheric removal process is oxidation by hydroxyl radicals. There is very little information on the reaction mechanism or reaction-path dynamics of atmospheric degradation of halogenated benzenes. Furthermore, the measured reaction rate constants are missing for the range of environmentally relevant temperatures, i.e. 230-330 K. A series of recent theoretical studies have provided those valuable missing information for fluorobenzene, chlorobenzene, hexafluorobenzene and hexachlorobenzene. Their comparative analysis has provided additional and more general insight into the mechanism of those important tropospheric degradation processes as well as into the mobility, transport and atmospheric fate of halogenated aromatic systems. It was demonstrated for the first time that the addition of hydroxyl radicals to monohalogenated as well as to perhalogenated benzenes proceeds indirectly, via a prereaction complex and its formation and dynamics have been characterized including the respective transition-state. However, in fluorobenzene and chlorobenzene reactions hydroxyl radical hydrogen is pointing approximately to the center of the aromatic ring while in the case of hexafluorobenzene and hexachlorobenzene, unexpectedly, the oxygen is directed towards the center of the aromatic ring. The reliable rate constants are now available for all environmentally relevant temperatures for the tropospheric oxidation of fluorobenzene, chlorobenzene, hexafluorobenzene and hexachlorobenzene while pentachlorophenol, a well-known organic micropollutant, seems to be a major stable product of tropospheric oxidation of hexachlorobenzene. Their calculated tropospheric lifetimes show that fluorobenzene and chlorobenzene are easily removed from the atmosphere and do not have long-range transport potential while hexafluorobenzene seems to be a potential POP chemical and hexachlorobenzene is clearly a typical persistent organic pollutant.

Comparison of Thiyl, Alkoxyl, and Alkyl Radical Addition to Double Bonds: The Unusual Contrasting Behavior of Sulfur and Oxygen Radical Chemistry

High-level ab initio calculations have been used to compare prototypical thiyl, alkoxyl, and alkyl radical addition reactions. Thiyl radical addition to the sulfur center of thioketones is exothermic and rapid, occurring with negative enthalpic barriers and only weakly positive Gibbs free energy barriers. In stark contrast, alkoxyl radical addition to the oxygen center of ketones is highly endothermic and occurs with very high reaction barriers, though these are also suppressed. On the basis of analysis of the corresponding alkyl radical additions to these substrates and the corresponding reactions of these heteroatom radicals with alkenes, it suggested that addition reactions involving thiyl radicals have low intrinsic barriers because their unpaired electrons are better able to undergo stabilizing resonance interactions with the π* orbitals of the substrate in the transition state.

Temperature and pressure dependent rate coefficients for the reaction of C2H4 + HO2 on the C2H4O2H potential energy surface

The potential energy surface (PES) for reaction C2H4 + HO2 was examined by using the quantum chemical methods. All rates were determined computationally using the CBS-QB3 composite method combined with conventional transition state theory(TST), variational transition-state theory (VTST) and Rice-Ramsberger-Kassel-Marcus/master-equation (RRKM/ME) theory. The geometries optimization and the vibrational frequency analysis of reactants, transition states, and products were performed at the B3LYP/CBSB7 level. The composite CBS-QB3 method was applied for energy calculations. The major product channel of reaction C2H4 + HO2 is the formation C2H4O2H via an OH(···)π complex with 3.7 kcal/mol binding energy which exhibits negative-temperature dependence. We further investigated the reactions related to this complex, which were ignored in previous studies. Thermochemical properties of the species involved in the reactions were determined using the CBS-QB3 method, and enthalpies of formation of species were compared with literature values. The calculated rate constants are in good agreement with those available from literature and given in modified Arrhenius equation form, which are serviceable in combustion modeling of hydrocarbons. Finally, in order to illustrate the effect for low-temperature ignition of our new rate constants, we have implemented them into the existing mechanisms, which can predict ethylene ignition in a shock tube with better performance.

Predicting the effect of angular momentum on the dissociation dynamics of highly rotationally excited radical intermediates

We present a model which accurately predicts the net speed distributions of products resulting from the unimolecular decomposition of rotationally excited radicals. The radicals are produced photolytically from a halogenated precursor under collision-free conditions so they are not in a thermal distribution of rotational states. The accuracy relies on the radical dissociating with negligible energetic barrier beyond the endoergicity. We test the model predictions using previous velocity map imaging and crossed laser-molecular beam scattering experiments that photolytically generated rotationally excited CD2CD2OH and C3H6OH radicals from brominated precursors; some of those radicals then undergo further dissociation to CD2CD2 + OH and C3H6 + OH, respectively. We model the rotational trajectories of these radicals, with high vibrational and rotational energy, first near their equilibrium geometry, and then by projecting each point during the rotation to the transition state (continuing the rotational dynamics at that geometry). This allows us to accurately predict the recoil velocity imparted in the subsequent dissociation of the radical by calculating the tangential velocities of the CD2CD2/C3H6 and OH fragments at the transition state. The model also gives a prediction for the distribution of angles between the dissociation fragments' velocity vectors and the initial radical's velocity vector. These results are used to generate fits to the previously measured time-of-flight distributions of the dissociation fragments; the fits are excellent. The results demonstrate the importance of considering the precession of the angular velocity vector for a rotating radical. We also show that if the initial angular momentum of the rotating radical lies nearly parallel to a principal axis, the very narrow range of tangential velocities predicted by this model must be convoluted with a J = 0 recoil velocity distribution to achieve a good result. The model relies on measuring the kinetic energy release when the halogenated precursor is photodissociated via a repulsive excited state but does not include any adjustable parameters. Even when different conformers of the photolytic precursor are populated, weighting the prediction by a thermal conformer population gives an accurate prediction for the relative velocity vectors of the fragments from the highly rotationally excited radical intermediates.

New insights into thermal decomposition of polycyclic aromatic hydrocarbon oxyradicals

Thermal decompositions of polycyclic aromatic hydrocarbon (PAH) oxyradicals on various surface sites including five-membered ring, free-edge, zigzag, and armchair have been systematically investigated by using ab initio density functional theory B3LYP/6-311+G(d,p) basis set. The calculation based on Hückel theory indicates that PAHs (3H-cydopenta[a]anthracene oxyradical) with oxyradicals on a five-membered ring site have high chemical reactivity. The rate coefficients of PAH oxyradical decomposition were evaluated by using Rice-Ramsperger-Kassel-Marcus theory and solving the master equations in the temperature range of 1500-2500 K and the pressure range of 0.1-10 atm. The kinetic calculations revealed that the rate coefficients of PAH oxyradical decomposition are temperature-, pressure-, and surface site-dependent, and the oxyradical on a five-membered ring is easier to decompose than that on a six-membered ring. Four-membered rings were found in decomposition of the five-membered ring, and a new reaction channel of PAH evolution involving four-membered rings is recommended.

Theoretical study on the gas phase reaction of propargyl alcohol with hydroxyl radical

The reaction of propargyl alcohol with hydroxyl radical has been studied extensively at CCSD(T)/aug-cc-pVTZ//MP2/cc-pVTZ level. This is the first time to gain a conclusive insight into the reaction mechanism and kinetics for this important reaction in detail. Two reaction mechanisms were revealed, namely addition/elimination and hydrogen abstraction mechanism. The reaction mechanism confirms that OH addition to C≡C triple bond forms the chemically activated adducts, IM1 (·CHCOHCH2OH) and IM2 (CHOH·CCH2OH), and the hydrogen abstraction pathways (-CH2OH bonded to the carbon atom and alcohol hydrogen) may occur via low barriers. Harmonic model of Rice-Ramsperger-Kassel-Marcus theory and variational transition state theory are used to calculate the overall and individual rate constants over a wide range of temperatures and pressures. The calculated rate constants are in good agreement with the experimental data. At atmospheric pressure with Ar as bath gas, IM1 (·CHCOHCH2OH) and IM2 (CHOH·CCH2OH) formed by collisional stabilization are dominant in the low temperature range. The production of CHCCHOH + H2O via hydrogen abstraction becomes dominate at higher temperature. The fraction of IM3 (CH2COHCH2·O) is very significant over the moderate temperature range.

Theoretical study on the gas phase reaction of allyl chloride with hydroxyl radical

The reaction of allyl chloride with the hydroxyl radical has been investigated on a sound theoretical basis. This is the first time to gain a conclusive insight into the reaction mechanism and kinetics for important pathways in detail. The reaction mechanism confirms that OH addition to the C=C double bond forms the chemically activated adducts, IM1 (CH2CHOHCH2Cl) and IM2 (CH2OHCHCH2Cl) via low barriers, and direct H-abstraction paths may also occur. Variational transition state model and multichannel RRKM theory are employed to calculate the temperature-, pressure-dependent rate constants. The calculated rate constants are in good agreement with the experimental data. At 100 Torr with He as bath gas, IM6 formed by collisional stabilization is the major products in the temperature range 200-600 K; the production of CH2CHCHCl via hydrogen abstractions becomes dominant at high temperatures (600-3000 K).

Mechanism and kinetic study of 3-fluoropropene with hydroxyl radical reaction

Potential energy surface for the reaction of hydroxyl radical (OH) with 3-fluoropropene (CH₂CHCH₂F) has been studied to evaluate the reaction mechanisms, possible products and rate constants. It has been shown that the CH₂CHCH₂F with OH reaction takes place via a barrierless addition/elimination and hydrogen abstraction mechanism. It is revealed for the first time that the initial step for the barrierless additional process involves a pre-reactive loosely bound complex (CR1) that is 1.60 kcal/mol below the energy of the reactants. Subsequently, the reaction bifurcates into two different pathways to form IM1 (CH₂CHOHCH₂F) and IM2 (CH₂OHCHCH₂F), which can decompose or isomerize to various products via complicated mechanisms. Variational transition state model and multichannel RRKM theory are employed to calculate the temperature-, pressure-dependent rate constants and branching ratios. At atmospheric pressure with He as bath gas, IM1 formed by collisional stabilization is dominated at T≤600 K; whereas the direct hydrogen abstraction leading to CH₂CHCHF and H₂O are the major products at temperatures between 600 and 3000 K, with estimated contribution of 72.9% at 1000 K. Furthermore, the predicted rate constants are in good agreement with the available experimental values.