Effects of Olefin Group and Its Position on the Kinetics for Intramolecular H-Shift and HO2 Elimination of Alkenyl Peroxy Radicals

Comprehensive Theoretical Study on Four Typical Intramolecular Hydrogen Shift Reactions of Peroxy Radicals: Multireference Character, Recommended Model Chemistry, and Kinetics

Intramolecular hydrogen shift reactions in peroxy radicals (RO₂^• → ^•QOOH) play key roles in the low-temperature combustion and in the atmospheric chemistry. In the present study, we found that a mild-to-moderate multireference character of a potential energy surface (PES) is widely present in four typical hydrogen shift reactions of peroxy radicals (RO₂^•, R = ethyl, vinyl, formyl methyl, and acetyl) by a systematic assessment based on the T₁ diagnostic, %TAE diagnostic, M diagnostic, and contribution of the dominant configuration of the reference CASSCF wavefunction (C₀²). To assess the effects of these inherent multireference characters on electronic structure calculations, we compared the PESs of the four reactions calculated by the multireference method CASPT2 in the complete basis set (CBS) limit, single-reference method CCSD(T)-F12, and single-reference-based composite method WMS. The results showed that ignoring the multireference character will introduce a mean unsigned deviation (MUD) of 0.46-1.72 kcal/mol from CASPT2/CBS results by using the CCSD(T)-F12 method or a MUD of 0.49-1.37 kcal/mol by WMS for three RO₂^• reactions (R = vinyl, formyl methyl, and acetyl) with a stronger multireference character. Further tests by single-reference Kohn-Sham (KS) density functional theory methods showed even larger deviations. Therefore, we specifically developed a new hybrid meta-generalized gradient approximation (GGA) functional M06-HS for the four typical H-shift reactions of peroxy radicals based on the WMS results for the ethyl peroxy radical reaction and on the CASPT2/CBS results for the others. The M06-HS method has an averaged MUD of 0.34 kcal/mol over five tested basis sets against the benchmark PESs, performing best in the tested 38 KS functionals. Last, in a temperature range of 200-3000 K, with the new functional, we calculated the high-pressure-limit rate coefficients of these H-shift reactions by the multi-structural variational transition-state theory with the small-curvature tunneling approximation (MS-CVT/SCT) and the thermochemical properties of all of the involved key radicals by the multi-structural torsional (MS-T) anharmonicity approximation method.

Cl-Initiated oxidation of methacrolein under NOx-free conditions studied by VUV photoionization mass spectrometry

The Cl-initiated oxidation of methacrolein (MACR, C₄H₆O) under NO_x-free conditions has been investigated in a fast flow tube by using a home-made vacuum ultraviolet (VUV) photoionization mass spectrometer complemented by high-level theoretical calculations. The key species such as intermediates and radicals together with products involved in the oxidation are observed online and confirmed in photoionization mass spectra. The reaction potential energy surfaces of the transient C₄H₅O and C₄H₆OCl radicals, formed from the hydrogen-abstraction reaction and the addition reaction of MACR with Cl atoms, with oxygen have been theoretically calculated to illuminate the formation of the peroxy radicals of C₄H₅OO₂ and C₄H₆OClO₂. The photoionization processes of these peroxy radicals, whose cations are not stable, and their individual self-reactions as well as bimolecular reactions with HO₂ radical are studied and discussed. In addition, kinetic experiments are also performed to get the time evolution of specific products and compared with theoretical models, providing a detailed insight into the reaction mechanism of the Cl-initiated oxidation of MACR.

High-Pressure-Limit Rate Coefficients for HO2 Elimination Reactions of Hydroperoxyalkenylperoxy Radicals based on the Reaction Class Transition State Theory

Thermokinetic parameters and transport parameters are of great importance to the combustion model and the reaction rate rules are of great importance to construct the combustion reaction mechanism for hydrocarbon fuels. The HO₂ elimination reaction class for hydroperoxyalkenylperoxy radicals is one of the key reaction classes for olefin, for which the rate coefficients are lacking. Therefore, the rate coefficients and rate rules of the HO₂ elimination reaction class for hydroperoxyalkenylperoxy radicals are studied in this work. The reaction class transition state theory (RC-TST) is used to calculate the rate coefficients. In addition, the HO₂ elimination reaction class of hydroperoxyalkenylperoxy radicals is divided into four subclasses depending upon the type of H-C_β bond that is broken in the reactant molecules, and the rate rules are calculated by taking the average of rate coefficients from a representative set of reactions in a subclass. The calculated kinetics data would be valuable for the construction of the combustion reaction mechanism for olefin.

A structure activity relationship for ring closure reactions in unsaturated alkylperoxy radicals

Terpenoids are an important class of multi-unsaturated volatile organic compounds emitted to the atmosphere. During their oxidation in the troposphere, unsaturated peroxy radicals are formed, which may undergo ring closure reactions by an addition of the radical oxygen atom on either of the carbons in the C[double bond, length as m-dash]C double bond. This study describes a quantum chemical and theoretical kinetic study of the rate of ring closure, finding that the reactions are comparatively fast with rates often exceeding 1 s-1 at room temperature, making these reactions competitive in low-NOx environments and allowing for continued autoxidation by ring closure. A structure-activity relationship (SAR) is presented for 5- to 8-membered ring closure in unsaturated RO2 radicals with aliphatic substituents, with some analysis of the impact of oxygenated substituents. H-migration in the cycloperoxide peroxy radicals formed after the ring closure was found to be comparatively slow for unsubstituted RO2 radicals. In the related cycloperoxide alkoxy radicals, migration of H-atoms implanted on the ring was similarly found to be slower than for non-cyclic alkoxy radicals and is typically not competitive against decomposition reactions that lead to cycloperoxide ring breaking. Ring closure reactions may constitute an important reaction channel in the atmospheric oxidation of terpenoids and could promote continued autoxidation, though the impact is likely to be strongly dependent on the specific molecular backbone.

Atmospheric Chemistry of Allylic Radicals from Isoprene: A Successive Cyclization-Driven Autoxidation Mechanism

The atmospheric chemistry of isoprene has broad implications for regional air quality and the global climate. Allylic radicals, taking 13-17% yield in the isoprene oxidation by ^•Cl, can contribute as much as 3.6-4.9% to all possible formed intermediates in local regions at daytime. Considering the large quantity of isoprene emission, the chemistry of the allylic radicals is therefore highly desirable. Here, we investigated the atmospheric oxidation mechanism of the allylic radicals using quantum chemical calculations and kinetics modeling. The results indicate that the allylic radicals can barrierlessly combine with O₂ to form peroxy radicals (RO₂^•). Under ≤100 ppt NO and ≤50 ppt HO₂^• conditions, the formed RO₂^• mainly undergo two times "successive cyclization and O₂ addition" to finally form the product fragments 2-alkoxy-acetaldehyde (C₂H₃O₂^•) and 3-hydroperoxy-2-oxopropanal (C₃H₄O₄). The presented reaction illustrates a novel successive cyclization-driven autoxidation mechanism. The formed 3-hydroperoxy-2-oxopropanal product is a new isomer of the atmospheric C₃H₄O₄ family and a potential aqueous-phase secondary organic aerosol precursor. Under >100 ppt NO condition, NO can mediate the cyclization-driven autoxidation process to form C₅H₇NO₃, C₅H₇NO₇, and alkoxy radical-related products. The proposed novel autoxidation mechanism advances our current understanding of the atmospheric chemistry of both isoprene and RO₂^•.

Pressure-dependent rate rules for cycloaddition, intramolecular H-shift, and concerted elimination reactions of alkenyl peroxy radicals at low temperature

The reactions of cycloaddition, intramolecular H-shift and concerted elimination of alkenyl peroxy radicals are three kinds of important reactions in the low temperature combustion of alkenes. In this study, the cycloaddition reactions are divided into classes considering endo-cycloaddition, exo-cycloaddition and the size of the transition states; the intramolecular H-shift reactions are divided into classes depending upon the ring size of the transition states and the type of C-H bonds from which the hydrogen atom is transferred; the concerted elimination reactions are divided into classes according to the type of H-C_βC_αOO bond that is broken. All geometry optimizations are performed at the B3LYP/6-31G(2df,p) level. With the electronic structure calculations being performed using the composite Gaussian-4 (G4) method, high pressure limit rate constants and pressure-dependent rate constants at pressures varying from 0.01 to 100 atm are calculated by using canonical transition state theory and the Rice-Ramsberger-Kassel-Marcus/master equation method, respectively. All rate constants are given in the form of the modified Arrhenius expression. The high pressure limit rate rules and the pressure-dependent rate rules are derived by averaging the rate constants of a representative set of reactions in each class. The results show that the rate rules for these three classes of reactions have a large uncertainty and the impact of the pressure on the rate constants increases as temperature increases.

Assessing the Molecular Basis of the Fuel Octane Scale: A Detailed Investigation on the Rate Controlling Steps of the Autoignition of Heptane and Isooctane

N-Heptane and 2,2,4-trimethylpentane (isooctane) are the key species in the modeling of ignition of hydrocarbon-based fuel formulations. Isooctane is knock-resistant whereas n-heptane is a very knock-prone hydrocarbon. It has been suggested that interconversion of their associated alkylperoxy and hydroperoxyalkyl species via hydrogen-transfer isomerization reaction is the key step to understand their different knocking behavior. In this work, the kinetics of unimolecular hydrogen-transfer reactions of n-heptylperoxy and isooctylperoxy are determined using canonical variational transition-state theory and multidimensional small curvature tunneling. Internal rotation of involved molecules is taken explicitly into account in the molecular partition function. The rate coefficients are calculated in the temperature range 300-900 K, relevant to low-temperature autoignition. The concerted HO₂ elimination is an important reaction that competes with some H-transfer and is associated with chain termination. Thus, the branching ratio between these reaction channels is analyzed. We show that variational and multidimensional tunneling effects cannot be neglected for the H-transfer reaction. In particular, the pre-exponential Arrhenius fitting parameter derived from our rate constants shows a strong dependence on the temperature, because tunneling increases quickly at temperatures below 500 K. On the basis of our results, the existing qualitative model for the reasons for different knock behavior observed for n-heptane and isooctane is quantitatively validated at the molecular level.

First-Principles Chemical Kinetic Modeling of Methyl trans-3-Hexenoate Epoxidation by HO2

The design of innovative combustion processes relies on a comprehensive understanding of biodiesel oxidation kinetics. The present study aims at unraveling the reaction mechanism involved in the epoxidation of a realistic biodiesel surrogate, methyl trans-3-hexenoate, by hydroperoxy radicals using a bottom-up theoretical kinetics methodology. The obtained rate constants are in good agreement with experimental data for alkene epoxidation by HO₂. The impact of temperature and pressure on epoxidation pathways involving H-bonded and non-H-bonded conformers was assessed. The obtained rate constant was finally implemented into a state-of-the-art detailed combustion mechanism, resulting in fairly good agreement with engine experiments.

Computational investigation into the gas-phase ozonolysis of the conjugated monoterpene α-phellandrene

Reaction with ozone is a major atmospheric sink for α-phellandrene, a monoterpene found in both indoor and outdoor environments, however experimental literature concerning the reaction is scarce. In this study, high-level G4(MP2) quantum chemical calculations are used to theoretically characterise the reaction of ozone with both double bonds in α-phellandrene for the first time. Results show that addition of ozone to the least substituted double bond in the conjugated system is preferred. Following addition, thermal and chemically activated unimolecular reactions, including the so-called hydroperoxide and ester or 'hot' acid channels, and internal cyclisation reactions, are characterised to major first generation products. Conjugation present in α-phellandrene allows two favourable Criegee intermediate reaction pathways to proceed that have not previously been considered in the literature; namely a 1,6-allyl resonance stabilised hydrogen shift and intramolecular dioxirane isomerisation to an epoxide. These channels are expected to play an important role alongside conventional routes in the ozonolysis of a-phellandrene. Computational characterisation of the potential energy surface thus provides insight into this previously unstudied system, and will aid future mechanism development and experimental interpretation involving α-phellandrene and structurally similar species, to which the results are expected to extend.

Theoretical Analysis of the Effect of C═C Double Bonds on the Low-Temperature Reactivity of Alkenylperoxy Radicals

Biodiesel contains a large proportion of unsaturated fatty acid methyl esters. Its combustion characteristics, especially its ignition behavior at low temperatures, have been greatly affected by these C═C double bonds. In this work, we performed a theoretical analysis of the effect of C═C double bonds on the low-temperature reactivity of alkenylperoxy radicals, the key intermediates from the low-temperature combustion of biodiesel. To understand how double bonds affect the fate of peroxy radicals, we selected three representative peroxy radicals from heptane, heptene, and heptadiene having zero, one, and two double C═C bonds, respectively, for study. The potential energy surfaces were explored at the CBS-QB3 level, and the reaction rate constants were computed using canonical/variational transition state theories. We have found that the double bond is responsible for the very different bond dissociation energies of the various types of C-H bonds, which in turn affect significantly the reaction kinetics of alkenylperoxy radicals.

Direct observation and kinetics of a hydroperoxyalkyl radical (QOOH)

Oxidation of organic compounds in combustion and in Earth’s troposphere is mediated by reactive species formed by the addition of molecular oxygen (O2) to organic radicals. Among the most crucial and elusive of these intermediates are hydroperoxyalkyl radicals, often denoted “QOOH.” These species and their reactions with O2 are responsible for the radical chain branching that sustains autoignition and are implicated in tropospheric autoxidation that can form low-volatility, highly oxygenated organic aerosol precursors. We report direct observation and kinetics measurements of a QOOH intermediate in the oxidation of 1,3-cycloheptadiene, a molecule that offers insight into both resonance-stabilized and nonstabilized radical intermediates. The results establish that resonance stabilization dramatically changes QOOH reactivity and, hence, that oxidation of unsaturated organics can produce exceptionally long-lived QOOH intermediates.

Kinetics of homoallylic/homobenzylic rearrangement reactions under combustion conditions

Homoallylic/homobenzylic radicals refer to typical radicals with the radical site located at the β position from the vinyl/phenyl group. These radicals are largely involved in combustion systems, such as the pyrolysis or oxidation of alkenes, cycloalkanes, and aromatics. The 1,2-vinyl/phenyl migration via two steps (cyclization/fission) is a peculiar reaction type for the homoallylic/homobenzylic radicals, entitled homoallylic/homobenzylic rearrangement, which has been studied by theoretical calculations including the Hirshfeld atomic charge analysis in the present work. With the help of rate constant calculations, the competition between this reaction channel and other possible pathways under combustion temperatures (500-2000 K) were evaluated. Analogous 1,3- and 1,4-vinyl/phenyl migration reactions for similar radicals with the radical sites located at the γ and δ positions from the vinyl/phenyl group were also computed. The results indicate that the 1,2-vinyl/phenyl migration is particularly important for the kinetics of unimolecular reactions of homoallylic radicals under 1500 K; nevertheless, it still has noticeable contribution at higher temperature. For those radicals with the radical site at the γ or δ positions, the respective 1,3- or 1,4-vinyl/phenyl migration channel plays an insignificant role under combustion conditions.

The gas-phase reaction of methane sulfonic acid with the hydroxyl radical without and with water vapor

The gas phase reaction between methane sulfonic acid (CH3SO3H; MSA) and the hydroxyl radical (HO), without and with a water molecule, was investigated with DFT-B3LYP and CCSD(T)-F12 methods. For the bare reaction we have found two reaction mechanisms, involving proton coupled electron transfer and hydrogen atom transfer processes that produce CH3SO3 and H2O. We also found a third reaction mechanism involving the double proton transfer process, where the products and reactants are identical. The computed rate constant for the oxidation process is 8.3 × 10(-15) cm(3) s(-1) molecule(-1). CH3SO3H forms two very stable complexes with water with computed binding energies of about 10 kcal mol(-1). The presence of a single water molecule makes the reaction between CH3SO3H and HO much more complex, introducing a new reaction that consists in the interchange of H2O between HO and CH3SO3H. Our kinetic calculations show that 99.5% of the reaction involves this interchange of the water molecule and, consequently, water vapor does not play any role in the oxidation reaction of methane sulfonic acid by the hydroxyl radical.

Biofuels combustion

This review describes major features of current research in renewable fuels derived from plants and from fatty acids. Recent and ongoing fundamental studies of biofuel molecular structure, oxidation reactions, and biofuel chemical properties are reviewed, in addition to combustion applications of biofuels in the major types of engines in which biofuels are used. Biofuels and their combustion are compared with combustion features of conventional petroleum-based fuels. Two main classes of biofuels are described, those consisting of small, primarily alcohol, fuels (particularly ethanol, n-butanol, and iso-pentanol) that are used primarily to replace or supplement gasoline and those derived from fatty acids and used primarily to replace or supplement conventional diesel fuels. Research efforts on so-called second- and third-generation biofuels are discussed briefly.

Unusually fast 1,6-h shifts of enolic hydrogens in peroxy radicals: formation of the first-generation C2 and C3 carbonyls in the oxidation of isoprene

In a theoretical investigation using the CBS-QB3//UB3LYP/6-31+G** method supported by higher-level computations such as CBS-QB3//UQCISD/6-31+G**, the 1,6-H shifts of the enolic hydrogen in peroxy radicals of the type Z-HO-CH═CH-CH(2)-OO(•) were found to face exceptionally low energy barriers of only about 11 kcal mol(-1)--i.e., 6-9 kcal mol(-1) lower than the barriers for similar shifts of alkane hydrogens--such that they can proceed at unequaled rates of order 10(5) to 10(6) s(-1) at ambient temperatures. The unusually low barriers for enolic 1,6-H shifts in peroxy radicals, characterized here for the first time to our knowledge, are rationalized. As cases in point, the secondary peroxy radicals Z-HO-CH═C(CH(3))-CH(OO(•))-CH(2)OH (case A) and Z-HO-CH═CH-C(CH(3))(OO(•))-CH(2)OH (case B) derived from the primary Z-δ-hydroxy-peroxy radicals in the oxidation of isoprene, are predicted to undergo 1,6-H shifts of their enolic hydrogens at TST-calculated rates in the range 270-320 K of k(T)(A) = 5.4 × 10(-4) × T(5.04) × exp(-1990/T) s(-1) and k(T)(B) = 109 × T(3.13) × exp(-3420/T) s(-1), respectively, i.e., 2.0 × 10(6) and 6.2 × 10(4) s(-1), respectively, at 298 K, far outrunning in all relevant atmospheric and laboratory conditions their reactions with NO proposed earlier as their dominant pathways (Dibble J. Phys. Chem. A 2004, 108, 2199). These fast enolic-H shifts are shown to provide the explanation for the first-generation formation of methylglyoxal + glycolaldehyde, and glyoxal + hydroxyacetone in the oxidation of isoprene under high-NO conditions, recently determined by several groups. However, under moderate- and low-NO atmospheric conditions, the fast interconversion and equilibration of the various thermally labile, initial peroxy conformers/isomers from isoprene and the isomerization of the initial Z-δ-hydroxy-peroxy radicals, both recently proposed by us (Peeters et al. Phys. Chem. Chem. Phys. 2009, 11, 5935), are expected to substantially reduce the yields of the small carbonyls at issue.