The multichannel n-propyl + O2 reaction surface: Definitive theory on a model hydrocarbon oxidation mechanism

Dimeric product formation in the self-reaction of small peroxy radicals using synchrotron radiation vacuum ultraviolet photoionization mass spectrometry

Peroxy radicals (RO2) are key reactive intermediates in atmospheric oxidation processes and yet their chemistry is not fully unraveled. Little is known about their structures and the structures of the dimeric products (ROOR) in the self-reaction of small RO2, which are among the most abundant RO2 in the atmosphere. The product branching ratios of ROOR and their atmospheric roles are still in controversy. Here, the self-reaction of propyl peroxy radicals (C3H7O2), a typical small RO2 radical in the atmosphere, has been studied using synchrotron radiation vacuum ultraviolet photoionization mass spectrometry. Both radical (C3H7O) and closed-shell molecular (C3H6O, C3H7OH, C3H7OOC3H7) products in the self-reaction are observed in photoionization mass spectra and their elusive isomers are definitely identified in mass-selected photoionization spectra. Three isomers of the C3H7OOC3H7 dimeric products, R1OOR1, R1OOR2, and R2OOR2 (R1 and R2 represent 1-C3H7 and 2-C3H7, respectively), as well as their complex structures have been determined for the first time. Kinetic experiments are performed and compared with chemical simulations to reveal the sources of specific products. The branching ratio of the C3H7OOC3H7 dimeric channel is measured at 10 ± 5%. This work demonstrates that the dimeric product formation in the self-reaction of small RO2 radicals is non-negligible and should provide valuable new insight into atmospheric modelling.

Ethynyl Radical Hydrogen Abstraction Energetics and Kinetics Utilizing High-Level Theory

The ethynyl radical, C2H, is found in a variety of different environments ranging from interstellar space and planetary atmospheres to playing an important role in the combustion of various alkynes under fuel-rich conditions. Hydrogen-atom abstraction reactions are common for the ethynyl radical in these contrasting environments. In this study, the C2H + HX → C2H2 + X, where HX = HNCO, trans-HONO, cis-HONO, C2H4, and CH3OH, reactions have been investigated at rigorously high levels of theory, including CCSD(T)-F12a/cc-pVTZ-F12. For the stationary points thus located, much higher levels of theory have been used, with basis sets as large as aug-cc-pV5Z and methods up to CCSDT(Q), and core correlation was also included. These molecules were chosen because they can be found in either interstellar or combustion environments. Various additive energy corrections have been included to converge the relative enthalpies of the stationary points to subchemical accuracy (≤0.5 kcal mol-1). Barriers predicted here (2.19 kcal mol-1 for the HNCO reaction and 0.47 kcal mol-1 for C2H4) are significantly lower than previous predictions. Reliable kinetics were acquired over a wide range of temperatures (50-5000 K), which may be useful for future experimental studies of these reactions.

The multichannel i-propyl + O2 reaction system: A model of secondary alkyl radical oxidation

The i-propyl + O2 reaction mechanism has been investigated by definitive quantum chemical methods to establish this system as a benchmark for the combustion of secondary alkyl radicals. Focal point analyses extrapolating to the ab initio limit were performed based on explicit computations with electron correlation treatments through coupled cluster single, double, triple, and quadruple excitations and basis sets up to cc-pV5Z. The rigorous coupled cluster single, double, and triple excitations/cc-pVTZ level of theory was used to fully optimize all reaction species and transition states, thus, removing some substantial flaws in reference geometries existing in the literature. The vital i-propylperoxy radical (MIN1) and its concerted elimination transition state (TS1) were found 34.8 and 4.4 kcal mol-1 below the reactants, respectively. Two β-hydrogen transfer transition states (TS2, TS2') lie above the reactants by (1.4, 2.5) kcal mol-1 and display large Born-Oppenheimer diagonal corrections indicative of nearby surface crossings. An α-hydrogen transfer transition state (TS5) is discovered 5.7 kcal mol-1 above the reactants that bifurcates into equivalent α-peroxy radical hanging wells (MIN3) prior to a highly exothermic dissociation into acetone + OH. The reverse TS5 → MIN1 intrinsic reaction path also displays fascinating features, including another bifurcation and a conical intersection of potential energy surfaces. An exhaustive conformational search of two hydroperoxypropyl (QOOH) intermediates (MIN2 and MIN3) of the i-propyl + O2 system located nine rotamers within 0.9 kcal mol-1 of the corresponding lowest-energy minima.

Concordant Mode Approach for Molecular Vibrations

The Concordant Mode Approach (CMA) is advanced as a novel hierarchy for increasing the system size and level of theory feasible for quantum chemical computations of harmonic vibrational frequencies. The key concept behind CMA is that transferrable, internal-coordinate normal modes computed at an appropriate lower level of theory (B) comprise a superb basis for converging to vibrational frequencies given by a higher level of theory (A). Accordingly, high-level harmonic frequencies can be evaluated via CMA from a collection of single-point energies that essentially scales linearly in the number of atoms, providing nearly order-of-magnitude CPU time speedups. The accuracy of CMA methods was established by comprehensive tests on over 120 molecules for target Level A = CCSD(T)/cc-pVTZ with auxiliary Level B choices of both CCSD(T)/cc-pVDZ and B3LYP/6-31G(2df,p). Remarkably, the frequency residuals given by the diagonal CMA-0A(nc) scheme exhibit mean absolute deviations (MADs) of only 0.2 cm^-1 and standard deviations less than 0.5 cm^-1; the corresponding zero-point vibrational energies (ZPVEs) have negligible errors in the vicinity of 0.3 cm^-1.

Diastereomers and Low-Temperature Oxidation

Diastereomers have historically been ignored when building kinetic mechanisms for combustion. Low-temperature oxidation kinetics, which continues to gain interest in both combustion and atmospheric communities, may be affected by the inclusion of diastereomers in radical chain-branching pathways. In this work, key intermediates and transition states lacking stereochemical specification in an existing diethyl ether low-temperature oxidation mechanism were replaced with their diastereomeric counterparts. Rate coefficients for reactions involving diastereomers were computed with ab initio transition state theory master equation calculations. The presence of diastereomers increased rate coefficients by factors of 1.2-1.6 across various temperatures and pressures. Ignition delay simulations incorporating these revised rate coefficients indicate that the diastereomers enhanced the overall reactivity of the mechanism by almost 15% and increased the peak ketohydroperoxide concentration by 30% in the negative temperature coefficient region at combustion-relevant pressures. These results provide an illustrative indication of the important role of stereomeric effects in oxidation kinetics.

Identification and Reactivity of s-cis,s-cis-Dihydroxycarbene, a New [CH2O2] Intermediate

We report the first preparation of the s-cis,s-cis conformer of dihydroxycarbene (1cc) by means of pyrolysis of oxalic acid, isolation of the lower-energy s-trans,s-trans (1tt) and s-cis,s-trans (1ct) product conformers at cryogenic temperatures in a N₂ matrix, and subsequent narrow-band near-infrared (NIR) laser excitation to give 1cc. Carbene 1cc converts quickly to 1ct via quantum-mechanical tunneling with an effective half-life of 22 min at 3 K. The potential energy surface features around 1 were pinpointed by convergent focal point analysis targeting the AE-CCSDT(Q)/CBS level of electronic structure theory. Computations of the tunneling kinetics confirm the time scale of the 1cc → 1ct rotamerization and suggest that direct 1cc → H₂ + CO₂ decomposition may also be a minor pathway. The intriguing latter possibility cannot be confirmed spectroscopically, but hints of it may be present in the measured kinetic profiles.

An Intramolecular Hydrogen-Shift in a Peroxy Radical at Cryogenic Temperatures: The Reaction of 2-Hydroxyphenyl Radical with O2

Peroxy radical hydrogen-shifts are pivotal elementary reaction steps in the oxidation of small hydrocarbons in autoignition and the lower atmosphere. Although these reactions are typically associated with a substantial barrier, we demonstrate that the [1,5]H-shift in the peroxy species derived from the 2-hydroxyphenyl radical 1 is so facile that it even proceeds rapidly in an argon matrix at 35 K through a proton-coupled electron transfer mechanism. Hydrogen-bound complexes of o-benzoquinone are identified as the main reaction products by infrared spectroscopy although their formation through O-O bond scission is hampered by a barrier of 11.9 kcal mol^-1 at the ROCCSD(T)/cc-pVTZ/UB3LYP/6-311G(d,p) level of theory.

Energetics and mechanisms for the acetonyl radical + O2 reaction: An important system for atmospheric and combustion chemistry

The acetonyl radical (^•CH₂COCH₃) is relevant to atmospheric and combustion chemistry due to its prevalence in many important reaction mechanisms. One such reaction mechanism is the decomposition of Criegee intermediates in the atmosphere that can produce acetonyl radical and OH. In order to understand the fate of the acetonyl radical in these environments and to create more accurate kinetics models, we have examined the reaction system of the acetonyl radical with O₂ using highly reliable theoretical methods. Structures were optimized using coupled cluster theory with singles, doubles, and perturbative triples [CCSD(T)] with an atomic natural orbital (ANO0) basis set. Energetics were computed to chemical accuracy using the focal point approach involving perturbative treatment of quadruple excitations [CCSDT(Q)] and basis sets as large as cc-pV5Z. The addition of O₂ to the acetonyl radical produces the acetonylperoxy radical, and multireference computations on this reaction suggest it to be barrierless. No submerged pathways were found for the unimolecular isomerization of the acetonylperoxy radical. Besides dissociation to reactants, the lowest energy pathway available for the acetonylperoxy radical is a 1-5 H shift from the methyl group to the peroxy group through a transition state that is 3.3 kcal mol^-1 higher in energy than acetonyl radical + O₂. The ultimate products from this pathway are the enol tautomer of the acetonyl radical along with O₂. Multiple pathways that lead to OH formation are considered; however, all of these pathways are predicted to be energetically inaccessible, except at high temperatures.

Riddles of the structure and vibrational dynamics of HO3 resolved near the ab initio limit

The hydridotrioxygen (HO₃) radical has been investigated in many previous theoretical and experimental studies over several decades, originally because of its possible relevance to the tropospheric HO_x cycle but more recently because of its fascinating chemical bonding, geometric structure, and vibrational dynamics. We have executed new, comprehensive research on this vexing molecule via focal point analyses (FPA) to approach the ab initio limit of optimized geometric structures, relative energies, complete quartic force fields, and the entire reaction path for cis-trans isomerization. High-order coupled cluster theory was applied through the CCSDT(Q) and even CCSDTQ(P) levels, and CBS extrapolations were performed using cc-pVXZ (X = 2-6) basis sets. The cis isomer proves to be higher than trans by 0.52 kcal mol^-1, but this energetic ordering is achieved only after the CCSDT(Q) milestone is reached; the barrier for cis → trans isomerization is a minute 0.27 kcal mol^-1. The FPA central r_e(O-O) bond length of trans-HO₃ is astonishingly long (1.670 Å), consistent with the semiexperimental r_e distance we extracted from microwave rotational constants of 10 isotopologues using FPA vibration-rotation interaction constants (α_i). The D₀(HO-O₂) dissociation energy converges to a mere 2.80 ± 0.25 kcal mol^-1. Contrary to expectation for such a weakly bound system, vibrational perturbation theory performs remarkably well with the FPA anharmonic force fields, even for the torsional fundamental near 130 cm^-1. Exact numerical procedures are applied to the potential energy function for the torsional reaction path to obtain energy levels, tunneling rates, and radiative lifetimes. The cis → trans isomerization occurs via tunneling with an inherent half-life of 1.4 × 10^-11 s and 8.6 × 10^-10 s for HO₃ and DO₃, respectively, thus resolving the mystery of why the cis species has not been observed in previous experiments executed in dissipative environments that allow collisional cooling of the trans-HO₃ product. In contrast, the pure ground eigenstate of the cis species in a vacuum is predicted to have a spontaneous radiative lifetime of about 1 h and 5 days for HO₃ and DO₃, respectively.