Computational Studies of Intramolecular Hydrogen Atom Transfers in the β-Hydroxyethylperoxy and β-Hydroxyethoxy Radicals

Theoretical kinetic studies on intramolecular H-migration reactions of peroxy radicals of diethoxymethane

Diethoxymethane (DEM), a promising carbon-neutral fuel, has high reactivity at low temperatures. The intramolecular hydrogen migration reaction of the DEM peroxy radicals can be viewed as a critical step in the low temperature oxidation mechanism of DEM. In this work, multistructural transition state theory (MS-TST) was utilized to calculate the high-pressure limit rate constants of 1,5, 1,6 and 1,7 H-migration reactions for DEM peroxy radicals. In addition to the tunneling effects and anharmonic effects, the intramolecular effects, including steric hindrance, intramolecular hydrogen bonding and conformational changes in reactants and transition states, are also considered in the rate constant calculations. The calculated energy barriers and rate constants demonstrated the substantial impact of intramolecular effects on the kinetics of H-migration reactions in DEM peroxy radicals. Specifically, the distinct configurations of transition states could potentially influence the reaction kinetics. The pressure-dependent rate constants are computed using system-specific quantum RRK theory. The calculated results show that the falloff effect of 1,5 and 1,6 H-migration reactions is more pronounced than that of the 1,7 H-migration reaction. The thermodynamics and kinetics presented in this study could be instrumental in understanding the low-temperature oxidation mechanism of DEM and might prove crucial for future research on comprehensively analyzing the autoignition behavior.

Computational Study on the Reaction of β-Hydroxyethylperoxy Radical with HO2 and Effects of Water Vapor

The reaction of β-hydroxyethylperoxy radical (β-HEP) and HO₂ with and without water was studied using quantum chemistry and kinetic calculations. The main products are HOCH₂CH₂OOH and ³O₂ for the reaction with and without water, while all other reaction channels can be neglected. The rate coefficients of the reaction follow negative temperature dependence. The pseudo-second-order rate coefficients are 2-4 orders of magnitude smaller for the reaction with saturated water vapor, indicating the negligible contribution of water in this reaction. This is probably also true for other peroxy radicals (except for HO₂), indicating that a large part of previous results on the water enhancement of reaction rate coefficients might have overestimated the influence of water.

Theoretical Study on Reactions of α-Site Hydroxyethyl and Hydroxypropyl Radicals with O2

α-Site alcohol radicals are the most important products of H-abstract reactions from alcohols since the hydroxyl moiety weakens the α-site C-H bond. Reactions between α-site alcohol radicals and O2 play an important role in combustion of alcohols, especially at relatively low temperatures. However, reliable reaction pathways and rate constants for these reactions are still lacking. Theoretical studies on reactions in α-hydroxyethyl radical (CH3C•HOH) + O2 and α-hydroxypropyl radical (C2H5C•HOH and CH3C•OHCH3) + O2 reaction systems are performed in this work. Pressure-dependent rate constants for the involved reactions in a wide range of temperatures are determined using the Rice-Ramsperger-Kassel-Marcus/master equation (RRKM/ME) method. Our results show that rate constants for reactions in the α-hydroxypropyl radical + O2 system are quite different from those in the CH3C•HOH + O2 system. Detailed reaction pathways for these reaction systems are clarified, although combustion characteristics of ethanol and propanol do not change much with the obtained rate constants for these reactions. Important reaction channels in producing enols, especially in the combustion of propanol, are also provided. The obtained rate constants for these reactions over a wide range of temperatures and pressures are helpful in developing combustion mechanisms for ethanol and propanol.

Criegee Intermediates in Autoxidation Reactions: Mechanistic Considerations

Products of Criegee intermediate (CI) chemistry were recently detected in radical chain autoxidation reactions involving β-hydroxyperoxy radicals. Here, we demonstrate by means of accurate G4 computations that direct scission of the latter to CIs and radical byproducts is thermodynamically highly unfavorable. Instead, the reaction becomes possible through a hydrogen abstraction reaction that could proceed by reversible formation of a dimeric tetroxide and a subsequent [1,6] hydrogen shift of the hydroxy hydrogen.

A Systematic Theoretical Kinetics Analysis for the Waddington Mechanism in the Low-Temperature Oxidation of Butene and Butanol Isomers

The Waddington mechanism, or the Waddington-type reaction pathway, is crucial for low-temperature oxidation of both alkenes and alcohols. In this study, the Waddington mechanism in the oxidation chemistry of butene and butanol isomers was systematically investigated. Fundamental quantum chemical calculations were conducted for the rate constants and thermodynamic properties of the reactions and species in this mechanism. Calculations were performed using two different ab initio solvers: Gaussian 09 and Orca 4.0.0, and two different kinetic solvers: PAPR and MultiWell, comprehensively. Temperature- and pressure-dependent rate constants were performed based on the transition state theory, associated with the Rice Ramsperger Kassel Marcus and master equation theories. Temperature-dependent thermochemistry (enthalpies of formation, entropy, and heat capacity) of all major species was also conducted, based on the statistical thermodynamics. Of the two types of reaction, dissociation reactions were significantly faster than isomerization reactions, while the rate constants of both reactions converged toward higher temperatures. In comparison, between two ab initio solvers, the barrier height difference among all isomerization and dissociation reactions was about 2 and 0.5 kcal/mol, respectively, resulting in less than 50%, and a factor of 2–10 differences for the predicted rate coefficients of the two reaction types, respectively. Comparing the two kinetic solvers, the rate constants of the isomerization reactions showed less than a 32% difference, while the rate of one dissociation reaction (P1 ↔ WDT12) exhibited 1–2 orders of magnitude discrepancy. Compared with results from the literature, both reaction rate coefficients (R4 and R5 reaction systems) and species’ thermochemistry (all closed shell molecules and open shell radicals R4 and R5) showed good agreement with the corresponding values obtained from the literature. All calculated results can be directly used for the chemical kinetic model development of butene and butanol isomer oxidation.

Intramolecular effects on the kinetics of unimolecular reactions of β-HOROO˙ and HOQ˙OOH radicals

A theoretical study describing the influence of intramolecular effects on the energy barriers and rate constants of unimolecular reactions involving β-HOROO˙ and HOQ˙OOH radicals is proposed. The reactions considered are HO2˙ elimination, the Waddington mechanism, H-shift, cyclic ether formation and β-scission. All the calculations are performed at the CBS-QB3 level of theory along with canonical transition state theory and statistical thermodynamics, including a specific treatment of hindered rotors. Several structural parameters are investigated, such as the location of the hydroxyl function in the cyclic transition states or the substitution of H atoms by alkyl groups on carbon atoms involved in the reaction coordinate. It is shown that these molecular systems involve numerous transition states, especially for reactions such as 1,5 or 1,6 H-shift, and that, a priori simplification is not possible. It is also shown that the position of the -OH group in the transition state can largely modify both the barrier heights and the rate constants. However, opposite trends can be observed depending on the competition between energetic and entropic effects. Similar observations are made when H atoms are replaced by methyl or alkyl groups. These results can largely be explained by intramolecular effects such as hydrogen bonds, stabilization effects (from -OH or -CH3 groups), steric influences and by the coupling between them. The last point renders the classic establishment of the structure-reactivity relationship challenging.

Variational transition state theory: theoretical framework and recent developments

This article reviews the fundamentals of variational transition state theory (VTST), its recent theoretical development, and some modern applications.

A large source of low-volatility secondary organic aerosol

Forests emit large quantities of volatile organic compounds (VOCs) to the atmosphere. Their condensable oxidation products can form secondary organic aerosol, a significant and ubiquitous component of atmospheric aerosol, which is known to affect the Earth's radiation balance by scattering solar radiation and by acting as cloud condensation nuclei. The quantitative assessment of such climate effects remains hampered by a number of factors, including an incomplete understanding of how biogenic VOCs contribute to the formation of atmospheric secondary organic aerosol. The growth of newly formed particles from sizes of less than three nanometres up to the sizes of cloud condensation nuclei (about one hundred nanometres) in many continental ecosystems requires abundant, essentially non-volatile organic vapours, but the sources and compositions of such vapours remain unknown. Here we investigate the oxidation of VOCs, in particular the terpene α-pinene, under atmospherically relevant conditions in chamber experiments. We find that a direct pathway leads from several biogenic VOCs, such as monoterpenes, to the formation of large amounts of extremely low-volatility vapours. These vapours form at significant mass yield in the gas phase and condense irreversibly onto aerosol surfaces to produce secondary organic aerosol, helping to explain the discrepancy between the observed atmospheric burden of secondary organic aerosol and that reported by many model studies. We further demonstrate how these low-volatility vapours can enhance, or even dominate, the formation and growth of aerosol particles over forested regions, providing a missing link between biogenic VOCs and their conversion to aerosol particles. Our findings could help to improve assessments of biosphere-aerosol-climate feedback mechanisms, and the air quality and climate effects of biogenic emissions generally.

Site-specific rate coefficients for reaction of OH with ethanol from 298 to 900 K

The rate coefficients for reactions of OH with ethanol and partially deuterated ethanols have been measured by laser flash photolysis/laser-induced fluorescence over the temperature range 298-523 K and 5-100 Torr of helium bath gas. The rate coefficient, k(1.1), for reaction of OH with C(2)H(5)OH is given by the expression k(1.1) = 1.06 × 10(-22)T(3.58) exp(1126/T) cm(3) molecule(-1) s(-1), and the values are in good agreement with previous literature. Site-specific rate coefficients were determined from the measured kinetic isotope effects. Over the temperature region 298-523 K abstraction from the hydroxyl site is a minor channel. The reaction is dominated by abstraction of the α hydrogens (92 ± 8)% at 298 K decreasing to (76 ± 9)% with the balance being abstraction at the β position where the errors are 2σ. At higher temperatures decomposition of the CH(2)CH(2)OH product from β abstraction complicates the kinetics. From 575 to 650 K, biexponential decays were observed, allowing estimates to be made for k(1.1) and the fractional production of CH(2)CH(2)OH. Above 650 K, decomposition of the CH(2)CH(2)OH product was fast on the time scale of the measured kinetics and removal of OH corresponds to reaction at the α and OH sites. The kinetics agree (within ±20%) with previous measurements. Evidence suggests that reaction at the OH site is significant at our higher temperatures: 47-53% at 865 K.

Computational studies of the isomerization and hydration reactions of acetaldehyde oxide and methyl vinyl carbonyl oxide

Alkene ozonolysis is a major source of hydroxyl radical (*OH), the most important oxidant in the troposphere. Previous experimental and computational work suggests that for many alkenes the measured *OH yields should be attributed to the combined impact of both chemically activated and thermalized syn-alkyl Criegee intermediates (CIs), even though the thermalized CI should be susceptible to trapping by molecules such as water. We have used RRKM/master equation and variational transition state theory calculations to quantify the competition between unimolecular isomerization and bimolecular hydration reactions for the syn and anti acetaldehyde oxide formed in trans-2-butene ozonolysis and for the CIs formed in isoprene ozonolysis possessing syn-methyl groups. Statistical rate theory calculations were based on quantum chemical data provided by the B3LYP, QCISD, and multicoefficient G3 methods, and thermal rate constants were corrected for tunneling effects using the Eckart method. At tropospheric temperatures and pressures, all thermalized CIs with syn-methyl groups are predicted to undergo 1,4-hydrogen shifts from 2 to 8 orders of magnitude faster than they react with water monomer at its saturation number density. For thermalized anti acetaldehyde oxide, the rates of dioxirane formation and hydration should be comparable.

Unimolecular beta-hydroxyperoxy radical decomposition with OH recycling in the photochemical oxidation of isoprene

A novel process in the photochemical oxidation of isoprene that recycles hydroxyl (OH) radicals has been identified using first-principles computational chemistry. Isoprene is the dominant biogenic volatile organic compound (VOC), and its oxidation controls chemistry in the forest boundary layer and is also thought to contribute to cloud formation in marine environments. The mechanism described here involves rapid unimolecular decomposition of the two major peroxy radicals (beta-hydroxyperoxy radicals) produced by OH-initiated isoprene oxidation. Peroxy radicals are well-known as key intermediates in VOC oxidation, but up to now were only thought to be destroyed in bimolecular reactions. The process described here leads to OH recycling with up to around 60% efficiency in environments with low levels of peroxy radicals and NO(x). In forested environments reaction of the beta-hydroxyperoxy radicals with HO2 is expected to dominate, with a small contribution from the mechanism described here. Peroxy radical decomposition will be more important in the unpolluted marine boundary layer, where lower levels of NO and HO2 are encountered.

HOx radical regeneration in the oxidation of isoprene

We propose, and quantify from first principles, two novel HO(x)-regenerating unimolecular reactions in isoprene oxidation, which are estimated to yield in pristine tropical forest conditions about 0.7 HO(2) and 0.03 OH radicals per isoprene oxidized; it is further argued that the photolabile coproduct of HO(2) can be a major source of OH, with a yield of the order of 1. The newly proposed chemistry could provide a rationalization for the unexpectedly high OH concentrations often observed in forested environments, such as over the Amazon forest in the recent Gabriel campaign.