Experimental and Master Equation Study of the Kinetics of OH + C2H2: Temperature Dependence of the Limiting High Pressure and Pressure Dependent Rate Coefficients

Dynamics studies for the multi-well and multi-channel reaction of OH with C2H2 on a full-dimensional global potential energy surface

The C2H2 + OH reaction is an important acetylene oxidation pathway in the combustion process, as well as a typical multi-well and multi-channel reaction. Here, we report an accurate full-dimensional machine learning-based potential energy surface (PES) for the C2H2 + OH reaction at the UCCSD(T)-F12b/cc-pVTZ-F12 level, based on about 475 000 ab initio points. Extensive quasi-classical trajectory (QCT) calculations were performed on the newly developed PES to obtain detailed dynamic data and analyze reaction mechanisms. Below 1000 K, the C2H2 + OH reaction produces H + OCCH2 and CO + CH3. With increasing temperature, the product channels H2O + C2H and H + HCCOH are accessible and the former dominates above 1900 K. It is found that the formation of H2O + C2H is dominated by a direct reaction process, while other channels belong to the indirect mechanism involving long-lived intermediates along the reaction pathways. At low temperatures, the C2H2 + OH reaction behaves like an unimolecular reaction due to the unique PES topographic features, of which the dynamic features are similar to the decomposition of energy-rich complexes formed by C2H2 + OH collision. The classification of trajectories that undergo different reaction pathways to generate each product and their product energy distributions were also reported in this work. This dynamic information may provide a deep understanding of the C2H2 + OH reaction.

Theoretical Investigations on the Hydroxyl-Initialized Oxidation of Hexafluoro-2-butyne in the Presence of Oxygen

Hexafluoro-2-butyne (C₄F₆) is a potential eco-friendly alternative gas in plasma, refrigerants, and electrical insulation applications. Mechanisms for the reactions of C₄F₆ with OH/O₂ have been revealed in detail using various theoretical methods including ROCBS-Q, RCCSD(T), multireference RS2C, and extrapolations to the complete basis-set limit with Aug-cc-pVnZ (n = T, Q, 5) basis sets. Rate coefficients and product branching ratios were predicted for a wide range of temperatures and pressures using the solution of master equations. The vibrationally adiabatic ground-state barrier for the initial C₄F₆ + OH association was best estimated to be 1.53-2.26 kcal/mol. Energetically preferable decomposition paths for the conformation-dependent C₄F₆OH adducts include six-center HF elimination, four-center proton migration, and C-C bond cleavage, but the collisional deactivation is dominant under ambient conditions. The subsequent oxidation of C₄F₆OH by O₂ bifurcates in two orientations and proceeds without any well-defined barrier followed by the successive isomerization/elimination steps, forming perfluorobiacetyl to regenerate OH radicals or trifluoroacetic acid with trifluoroacetyl radicals. The OH-recycling path accounts for a branching ratio of 70% under ambient conditions. Theoretical rate coefficients are in good agreement with the available experimental results. The effect of fluorination on the reactivity of alkynes toward OH/O₂ is discussed.

Atmospheric Ring-Closure and Dehydration Reactions of 1,4-Hydroxycarbonyls in the Gas Phase: The Impact of Catalysts

1,4-Hydroxycarbonyls can potentially undergo sequential reactions involving cyclization followed by dehydration to form dihydrofurans. As dihydrofurans contain a double bond, they are highly reactive toward atmospheric oxidants such as OH, O₃, and NO₃. In the present study, we use ab initio calculations to examine the impact of various atmospheric catalysts on the energetics and kinetics of the gas-phase cyclization and dehydration reaction steps associated with 4-hydroxybutanal, a prototypical 1,4-hydroxycarbonyl molecule. The cyclization step transforms 4-hydroxybutanal into 2-hydroxytetrahydrofuran, which can subsequently undergo dehydration to form 2,3-dihydrofuran. As the barriers associated with the cyclization and dehydration steps for 4-hydroxybutanal are, respectively, 34.8 and 63.0 kcal/mol in the absence of a catalyst, both reaction steps are inaccessible under atmospheric conditions in the gas phase. However, the presence of a suitable catalyst can significantly reduce the reaction barriers, and we have examined the impact of a single molecule of H₂O, HO₂ radical, HC(O)OH, HNO₃, and H₂SO₄ on these reactions. We find that H₂SO₄ reduces the reaction barriers the greatest, with the barrier for the cyclization step being reduced to -13.1 kcal/mol and that for the dehydration step going down to 9.2 kcal/mol, measured relative to their respective separated starting reactants. Interestingly, our kinetic study shows that HNO₃ gives the fastest rate due to the combined effects of a larger atmospheric concentration and a reduced barrier. Thus, our study suggests that, with acid catalysis, the cyclization reaction step can readily occur for 1,4-hydroxycarbonyls in the gas phase. Because the dehydration step exhibits a significant barrier even with acid catalysis, the 2-hydroxytetrahydrofuran products, once formed, are likely lost through their reaction with OH radicals in the atmosphere. We have investigated the reaction pathways and the rate constant for this bimolecular reaction in the presence of excess molecular oxygen (³O₂), as it would occur under tropospheric conditions, using computational chemistry over the 200-300 K temperature range. We find that the main products from these OH-initiated oxidation reactions are succinaldehyde + HO₂ and 2,3-dihydro-2-furanol + HO₂.

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.

Unimolecular decomposition kinetics of the stabilised Criegee intermediates CH2OO and CD2OO

Decomposition kinetics of stabilised CH2OO and CD2OO Criegee intermediates have been investigated as a function of temperature (450-650 K) and pressure (2-350 Torr) using flash photolysis coupled with time-resolved cavity-enhanced broadband UV absorption spectroscopy. Decomposition of CD2OO was observed to be faster than CH2OO under equivalent conditions. Production of OH radicals following CH2OO decomposition was also monitored using flash photolysis with laser-induced fluorescence (LIF), with results indicating direct production of OH in the v = 0 and v = 1 states in low yields. Master equation calculations performed using the Master Equation Solver for Multi-Energy well Reactions (MESMER) enabled fitting of the barriers for the decomposition of CH2OO and CD2OO to the experimental data. Parameterisations of the decomposition rate coefficients, calculated by MESMER, are provided for use in atmospheric models and implications of the results are discussed. For CH2OO, the MESMER fits require an increase in the calculated barrier height from 78.2 kJ mol-1 to 81.8 kJ mol-1 using a temperature-dependent exponential down model for collisional energy transfer with ΔEdown = 32.6(T/298 K)1.7 cm-1 in He. The low- and high-pressure limit rate coefficients are k1,0 = 3.2 × 10-4(T/298)-5.81exp(-12 770/T) cm3 s-1 and k1,∞ = 1.4 × 1013(T/298)0.06exp(-10 010/T) s-1, with median uncertainty of ∼12% over the range of experimental conditions used here. Extrapolation to atmospheric conditions yields k1(298 K, 760 Torr) = 1.1+1.5-1.1 × 10-3 s-1. For CD2OO, MESMER calculations result in ΔEdown = 39.6(T/298 K)1.3 cm-1 in He and a small decrease in the calculated barrier to decomposition from 81.0 kJ mol-1 to 80.1 kJ mol-1. The fitted rate coefficients for CD2OO are k2,0 = 5.2 × 10-5(T/298)-5.28exp(-11 610/T) cm3 s-1 and k2,∞ = 1.2 × 1013(T/298)0.06exp(-9800/T) s-1, with overall error of ∼6% over the present range of temperature and pressure. The extrapolated k2(298 K, 760 Torr) = 5.5+9.2-5.5 × 10-3 s-1. The master equation calculations for CH2OO indicate decomposition yields of 63.7% for H2 + CO2, 36.0% for H2O + CO and 0.3% for OH + HCO with no significant dependence on temperature between 400 and 1200 K or pressure between 1 and 3000 Torr.

Dissociative Photoionization of the Elusive Vinoxy Radical

These experiments report the dissociative photoionization of vinoxy radicals to m/z = 15 and 29. In a crossed laser-molecular beam scattering apparatus, we induce C-Cl bond fission in 2-chloroacetaldehyde by photoexcitation at 157 nm. Our velocity measurements, combined with conservation of angular momentum, show that 21% of the C-Cl photofission events form vinoxy radicals that are stable to subsequent dissociation to CH₃ + CO or H + ketene. Photoionization of these stable vinoxy radicals, identified by their velocities, which are momentum-matched with the higher-kinetic-energy Cl atom photofragments, shows that the vinoxy radicals dissociatively photoionize to give signal at m/z = 15 and 29. We calibrated the partial photoionization cross section of vinoxy to CH₃⁺ relative to the bandwidth-averaged photoionization cross section of the Cl atom at 13.68 eV to put the partial photoionization cross sections on an absolute scale. The resulting bandwidth-averaged partial cross sections are 0.63 and 1.3 Mb at 10.5 and 11.44 eV, respectively. These values are consistent with the upper limit to the cross section estimated from a study by Savee et al. on the O(³P) + propene bimolecular reaction. We note that the uncertainty in these values is primarily dependent on the signal attributed to C-Cl primary photofission in the m/z = 35 (Cl⁺) time-of-flight data. While the value is a rough estimate, the bandwidth-averaged partial photoionization cross section of vinoxy to HCO⁺ calculated from the signal at m/z = 29 at 11.53 eV is approximately half that of vinoxy to CH₃⁺. We also present critical points on the potential energy surface of the vinoxy cation calculated at the G4//B3LYP/6-311++G(3df,2p) level of theory to support the observation of dissociative ionization of vinoxy to both CH₃⁺ and HCO⁺.

An Experimental and Master Equation Study of the Kinetics of OH/OD + SO2: The Limiting High-Pressure Rate Coefficients

The kinetics of the reaction OH/OD + SO₂ were studied using a laser flash photolysis/laser-induced fluorescence technique. Evidence for two-photon photolysis of SO₂ at 248 nm is presented and quantified, and which appears to have been evident to some extent in most previous photolysis studies, potentially leading to values for the rate coefficient k₁ that are too large. The kinetics of the reaction OH(v = 0) + SO₂ (T = 295 K, p = 25-300 torr) were measured under conditions where SO₂ photolysis was taken into account. These results, together with literature data, were modeled using a master equation analysis. This analysis highlighted problems with the literature data: the rate coefficients derived from flash photolysis data were generally too high and from the flow tube data too low. Our best estimate of the high-pressure limiting rate coefficient k₁^∞ was obtained from selected data and gives a value of (7.8 ± 2.2) × 10^-13 cm³ molecule^-1 s^-1, which is lower than that recommended in the literature. A parametrized form of k₁([N₂],T) is provided. The OD(v = 0) + SO₂ (T = 295 K, p = 25-300 torr) data are reported for the first time, and master equation analysis reinforces our assignment of k₁^∞.

An Experimental Study of the Kinetics of OH/OD(v = 1,2,3) + SO2: The Limiting High-Pressure Rate Coefficients as a Function of Temperature

The kinetics of the reaction OH/OD(v = 1,2,3) + SO₂ were studied using a photolysis/laser-induced fluorescence technique. The rate coefficients OH/OD(v = 1,2,3) + SO₂, k₁, over the temperature range of 295-810 K were used to determine the limiting high-pressure limit k₁^∞. This method is usually applicable if the reaction samples the potential well of the adduct HOSO₂ and if intramolecular vibrational relaxation is fast. In the present case, however, the rate coefficients showed an additional fast removal contribution as evidenced by the increase in k₁ with vibrational level; this behavior together with its temperature dependence is consistent with the existence of a weakly bound complex on the potential energy surface prior to adduct formation. The data were analyzed using a composite mechanism that incoporates energy-transfer mechanisms via both the adduct and the complex, and yielded a value of k₁^∞(295 K) equal to (7.2 ± 3.3) × 10^-13 cm³ molecule^-1 s^-1 (errors at 1σ), a factor of between 2 and 3 smaller than the current recommended IUPAC and JPL values of (2.0_-1.0^+2.0) and (1.6 ± 0.4) × 10^-12 cm³ molecule^-1 s^-1 at 298 K, respectively, although the error bars do overlap. k₁^∞ was observed to only depend weakly on temperature. Further evidence for a smaller k₁^∞ is presented in the companion paper.

Obtaining effective rate coefficients to describe the decomposition kinetics of the corannulene oxyradical at high temperatures

This work analyzes the effect of overlapping eigenvalues on the high-temperature kinetics of a large oxyradical based on master equation solutions.

The Onset of H + Ketene Products from Vinoxy Radicals Prepared by Photodissociation of Chloroacetaldehyde at 157 nm

We investigate the unimolecular dissociation of the vinoxy radical (CH2CHO) prepared with high internal energy imparted from the photodissociation of chloroacetaldehyde (CH2ClCHO) at 157 nm. Using a velocity map imaging apparatus, we measured the speed distribution of the recoiling chlorine atoms, Cl((2)P3/2) and Cl((2)P1/2), and derived from this the resulting distribution of kinetic energy, P(ET), imparted to the Cl + vinoxy fragments upon dissociation. Using conservation of energy, the distribution of kinetic energy was used to determine the total internal energy distribution in the radical. The P(ET) derived for the C-Cl bond fission presented in this work suggests the vinoxy radicals are mostly formed in the à state. We also took ion images at m/z = 42 and m/z = 15 to characterize the branching between the unimolecular dissociation channels of the vinoxy radical to H + ketene and methyl + CO products. Our results show a marked change in the branching ratio between the two channels from the previous study on the photodissociation of chloroacetaldehyde at 193 nm by Miller et al. (J. Chem. Phys., 2004, 121, 1830) in that the production of ketene is now favored over the production of methyl. To help analyze the data, we developed a model for the branching between the two channels that takes into account how the change in rotational energy en route to the products affects the vibrational energy available to surmount the barriers to the channels. The model predicts the portion of the C-Cl bond fission P(ET) that produces dissociative vinoxy radicals, then predicts the branching ratio between the H + ketene and CH3 + CO product channels at each ET. The model uses Rice-Ramsperger-Kassel-Marcus rate constants at the correct sums and densities of vibrational states while accounting for angular momentum conservation. We find that the predicted portion of the P(ET) that produces H + ketene products best fits the experimental portion (that we derive by taking advantage of conservation of momentum) if we use a barrier height for the H + ketene channel that is 4.0 ± 0.5 kcal/mol higher than the isomerization barrier en route to CH3 + CO products. Using the G4 computed isomerization barrier of 40.6 kcal/mol, this gives an experimentally determined barrier to the H + ketene channel of 44.6 kcal/mol. From these calculations, we also predict the branching ratio between the H + ketene and methyl + CO channels to be ∼2.1:1.

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.

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.

Interactions between theory and experiment in the investigation of elementary reactions of importance in combustion

Elementary reactions are a central component of models of combustion processes. Rate constants and channel yields are needed for those models. Both experimental and theoretical methods used to determine such rate data are discussed in this tutorial review, which is of interest to reaction kinetics and combustion engineering communities. Applications to combustion present particular problems because the conditions required can be well outside the ranges of temperature and pressure accessible to experiment, and the rate data can show a complex dependence on conditions. Under these conditions, the interplay between theory and experiment becomes important.

H Atom Yields from the Reactions of CN Radicals with C2H2, C2H4, C3H6, trans-2-C4H8, and iso-C4H8

The kinetics and H atom channel yield at both 298 and 195 K have been determined for reactions of CN radicals with C2H2 (1.00+/-0.21, 0.97+/-0.20), C2H4 (0.96+/-0.032, 1.04+/-0.042), C3H6 (pressure dependent), iso-C4H8 (pressure dependent), and trans-2-C4H8 (0.039+/-0.019, 0.029+/-0.047) where the first figure in each bracket is the H atom yield at 298 K and the second is that at 195 K. The kinetics of all reactions were studied by monitoring both CN decay and H atom growth by laser-induced fluorescence at 357.7 and 121.6 nm, respectively. The results are in good agreement with previous studies where available. The rate coefficients for the reaction of CN with trans-2-butene and iso-butene have been measured at 298 and 195 K for the first time, and the rate coefficients are as follows: k298K=(2.93+/-0.23)x10(-10) cm3 molecule(-1) s(-1), k195K=(3.58+/-0.43)x10(-10) cm3 molecule(-1) s(-1) and k298K=(3.17+/-0.10)x10(-10) cm3 molecule(-1) s(-1), k195K=(4.32+/-0.35)x10(-10) cm3 molecule(-1) s(-1), respectively, where the errors represent a combination of statistical uncertainty (2sigma) and an estimate of possible systematic errors. A potential energy surface for the CN+C3H6 reaction has been constructed using G3X//UB3LYP electronic structure calculations identifying a number of reaction channels leading to either H, CH3, or HCN elimination following the formation of initial addition complexes. Results from the potential energy surface calculations have been used to run master equation calculations with the ratio of primary:secondary addition, the average amount of downward energy transferred in a collision DeltaEd, and the difference in barrier heights between H atom elimination and an H atom 1, 2 migration as variable parameters. Excellent agreement is obtained with the experimental 298 K H atom yields with the following parameter values: secondary addition complex formation equal to 80%, DeltaEd=145 cm(-1), and the barrier height for H atom elimination set 5 kJ mol(-1) lower than the barrier for migration. Finally, very low temperature master equation simulations using the best fit parameters have been carried out in an increased precision environment utilizing quad-double and double-double arithmetic to predict H and CH3 yields for the CN+C3H6 reaction at temperatures and pressures relevant to Titan. The H and CH3 yields predicted by the master equation have been parametrized in a simple equation for use in modeling.

A combined experimental and theoretical study of the reaction between methylglyoxal and OH/OD radical: OH regeneration

Experimental studies have been conducted to determine the rate coefficient and mechanism of the reaction between methylglyoxal (CH(3)COCHO, MGLY) and the OH radical over a wide range of temperatures (233-500 K) and pressures (5-300 Torr). The rate coefficient is pressure independent with the following temperature dependence: k(3)(T) = (1.83 +/- 0.48) x 10(-12) exp((560 +/- 70)/T) cm(3) molecule(-1) s(-1) (95% uncertainties). Addition of O(2) to the system leads to recycling of OH. The mechanism was investigated by varying the experimental conditions ([O(2)], [MGLY], temperature and pressure), and by modelling based on a G3X potential energy surface, rovibrational prior distribution calculations and master equation RRKM calculations. The mechanism can be described as follows: Addition of oxygen to the system shows that process (4) is fast and that CH(3)COCO completely dissociates. The acetyl radical formed from reaction (4) reacts with oxygen to regenerate OH radicals (5a). However, a significant fraction of acetyl radical formed by reaction (R4) is sufficiently energised to dissociate further to CH(3) + CO (R4b). Little or no pressure quenching of reaction (R4b) was observed. The rate coefficient for OD + MGLY was measured as k(9)(T) = (9.4 +/- 2.4) x 10(-13) exp((780 +/- 70)/T) cm(3) molecule(-1) s(-1) over the temperature range 233-500 K. The reaction shows a noticeable inverse (k(H)/k(D) < 1) kinetic isotope effect below room temperature and a slight normal kinetic isotope effect (k(H)/k(D) > 1) at high temperature. The potential atmospheric implications of this work are discussed.