Measurements, Theory, and Modeling of OH Formation in Ethyl + O2 and Propyl + O2 Reactions

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.

Stereo-electronic factors influencing the stability of hydroperoxyalkyl radicals: transferability of chemical trends across hydrocarbons and ab initio methods

The hydroperoxyalkyl radicals (˙QOOH) are known to play a significant role in combustion and tropospheric processes, yet their direct spectroscopic detection remains challenging. In this study, we investigate molecular stereo-electronic effects influencing the kinetic and thermodynamic stability of a ˙QOOH along its formation path from the precursor, alkylperoxyl radical (ROO˙), and the depletion path resulting in the formation of cyclic ether + ˙OH. We focus on reactive intermediates encountered in the oxidation of acyclic hydrocarbon radicals: ethyl, isopropyl, isobutyl, tert-butyl, neopentyl, and their alicyclic counterparts: cyclohexyl, cyclohexenyl, and cyclohexadienyl. We report reaction energies and barriers calculated with the highly accurate method Weizmann-1 (W1) for the channels: ROO˙ ⇌ ˙QOOH, ROO˙ ⇌ alkene + ˙OOH, ˙QOOH ⇌ alkene + ˙OOH, and ˙QOOH ⇌ cyclic ether + ˙OH. Using W1 results as a reference, we have systematically benchmarked the accuracy of popular density functional theory (DFT), composite thermochemistry methods, and an explicitly correlated coupled-cluster method. We ascertain inductive, resonance, and steric effects on the overall stability of ˙QOOH and computationally investigate the possibility of forming more stable species. With new reactions as test cases, we probe the capacity of various ab initio methods to yield quantitative insights on the elementary steps of combustion.

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.

Rate constant and branching ratio of the reaction of ethyl peroxy radicals with methyl peroxy radicals

The cross-reaction of ethyl peroxy radicals (C₂H₅O₂) with methyl peroxy radicals (CH₃O₂) (R1) has been studied using laser photolysis coupled to time resolved detection of the two different peroxy radicals by continuous wave cavity ring down spectroscopy (cw-CRDS) in their AÃ-X̃ electronic transition in the near-infrared region, C₂H₅O₂ at 7602.25 cm^-1, and CH₃O₂ at 7488.13 cm^-1. This detection scheme is not completely selective for both radicals, but it is demonstrated that it has great advantages compared to the widely used, but unselective UV absorption spectroscopy. Peroxy radicals were generated from the reaction of Cl-atoms with the appropriate hydrocarbon (CH₄ and C₂H₆) in the presence of O₂, whereby Cl-atoms were generated by 351 nm photolysis of Cl₂. For different reasons detailed in the manuscript, all experiments were carried out under excess of C₂H₅O₂ over CH₃O₂. The experimental results were best reproduced by an appropriate chemical model with a rate constant for the cross-reaction of k = (3.8 ± 1.0) × 10^-13 cm³ s^-1 and a yield for the radical channel, leading to CH₃O and C₂H₅O, of (ϕ_1a = 0.40 ± 0.20).

Investigating the Association Reactions of HOCH2CO and HOCHCHO with O2: A Quantum Computational and Master Equation Study

Glycolaldehyde, HOCH₂CHO, is an important multifunctional atmospheric trace gas formed in the oxidation of ethylene and isoprene and emitted directly from burning biomass. The initial step in the atmospheric photooxidation of HOCH₂CHO yields HOCH₂CO and HOCHCHO radicals; both of these radicals react rapidly with O₂ in the troposphere. This study presents a comprehensive theoretical investigation of the HOCH₂CO + O₂ and HOCHCHO + O₂ reactions using high-level quantum chemical calculations and energy-grained master equation simulations. The HOCH₂CO + O₂ reaction results in the formation of a HOCH₂C(O)O₂ radical, while the HOCHCHO + O₂ reaction yields (HCO)₂ + HO₂. Density functional theory calculations have identified two open unimolecular pathways associated with the HOCH₂C(O)O₂ radical that yield HCOCOOH + OH or HCHO + CO₂ + OH products; the former novel bimolecular product pathway has not been previously reported in the literature. Master equation simulations based on the potential energy surface calculated here for the HOCH₂CO + O₂ recombination reaction support experimental product yield data from the literature and indicate that, even at total pressures of 1 atm, the HOCH₂CO + O₂ reaction yields ∼11% OH at 298 K.

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.

Spiers Memorial Lecture: theory of unimolecular reactions

One hundred years ago, at an earlier Faraday Discussion meeting, Lindemann presented a mechanism that provides the foundation for contemplating the pressure dependence of unimolecular reactions. Since that time, our...

Vacuum ultraviolet photochemistry of the conformers of the ethyl peroxy radical

We study the conformers of the ethyl peroxy radical (C₂H₅O₂), the simplest peroxy radical having more than one conformer, by combining synchrotron radiation vacuum ultraviolet (VUV) photoionization mass spectrometry with theoretical calculations. The ethyl peroxy radical is formed in a microwave discharge flow tube through the reaction of the ethyl radical (C₂H₅) with oxygen molecules, where C₂H₅ is generated via the hydrogen-abstraction reaction of ethane with fluorine atoms. Two kinds of C₂H₅⁺, originating from photoionization of C₂H₅ and from dissociative photoionization of C₂H₅O₂, whose cation is not stable, have been identified and separated in photoionization mass spectra. The photoionization spectrum corresponding to C₂H₅O₂ is obtained and assigned with Franck-Condon calculations. The present findings show that the gauche conformer (G-C₂H₅O₂) of C₂H₅O₂ has favorable Franck-Condon factors in the ionization transitions, whereas the contribution of the trans conformer (T-C₂H₅O₂) to the photoionization spectrum is minor or negligible due to its large geometric changes in the photoionization process. Moreover, the reason for the instability of C₂H₅O₂⁺ and its detailed dissociation mechanisms have been unraveled with the aid of the calculated potential energy curves.

Accidental Combustion Phenomena at Cryogenic Conditions

The presented state of the art can be intended as an overview of the current understandings and the remaining challenges on the phenomenological aspects involving systems operating at ultra-low temperature, which typically characterize the cryogenic fuels, i.e., liquefied natural gas and liquefied hydrogen. To this aim, thermodynamic, kinetic, and technological aspects were included and integrated. Either experimental or numerical techniques currently available for the evaluation of safety parameters and the overall reactivity of systems at cryogenic temperatures were discussed. The main advantages and disadvantages of different alternatives were compared. Theoretical background and suitable models were reported given possible implementation to the analyzed conditions. Attention was paid to models describing peculiar phenomena mainly relevant at cryogenic temperatures (e.g., para-to-ortho transformation and thermal stratification in case of accidental release) as well as critical aspects involving standard phenomena (e.g., ultra-low temperature combustion and evaporation rate).

Ab initio rate coefficients for reactions of 2,5-dimethylhexyl isomers with O2: temperature- and pressure-dependent branching ratios

Chemical kinetics of O2-addition to alkyl radicals (R), termed first O2-addition in the oxidation mechanism of alkanes, are of central importance to next-generation combustion strategies designed for operations in the low- to intermediate-temperature region (<1000 K). In the present work, stationary points on potential energy surfaces (PES), temperature- and pressure-dependent rate coefficients, and branching fractions of product formation from R + O2 reactions initiated by the addition of molecular oxygen (3O2) to the three alkyl radicals of a branched alkane, 2,5-dimethylhexane, are reported. The stationary points were determined utilizing ab initio/DFT methods and the reaction energies were computed using the composite CBS-QB3 method. Rice-Ramsperger-Kassel-Marcus (RRKM)/master equation (ME) calculations were employed to compute rate coefficients, from which branching fractions were determined over the pressure range of 10-3-20 atm and the temperature range of 400-900 K on three different surfaces. The quantum chemistry results reveal several distinct features. For the addition of O2 to the tertiary alkyl radical 2,5-dimethylhex-2-yl, the most energetically favorable channel leads to the formation of 2,2,5,5,-tetramethyl-tetrahydrofuran, a cyclic ether intermediate formed coincident with OH in a chain-propagating step from the decomposition of tertiary-tertiary hydroperoxyalkyl (QOOH). On the R + O2 surface of the secondary radical, 2,5-dimethylhex-3-yl, the pathways for the formation of methyl-propanal + iso-butene + OH via concerted C-C and O-O bond scission of tertiary QOOH and that of cyclic ether + OH are the most energetically favorable pathways. The R + O2 surface for the reaction of the primary radical, 2,5-dimethylhex-1-yl, reveals two competitive chain-propagation channels, leading to 2-iso-propyl-4-methyl-tetrahydrofuran + OH and 2,2,5-trimethyltetrahydropyran + OH. Below 100 Torr, the formation of the aforementioned species dominates the respective total R + O2 rate coefficient, while at pressures above 1 atm collisionally stabilized alkylperoxy (ROO) dominates at the temperatures considered here. The results of this study are in very good agreement with the experimentally measured intermediates and products of the 2,5-dimethylhexyl radical + O2 reaction.

Theoretical Study of Radical-Molecule Reactions with Negative Activation Energies in Combustion: Hydroxyl Radical Addition to Alkenes

Many of the radical-molecule reactions are nonelementary reactions with negative activation energies, which usually proceed through two steps. They exist extensively in the atmospheric chemistry and hydrocarbon fuel combustion, so they are extensively studied both theoretically and experimentally. At the same time, various models, such as a two transition state model, a steady-state model, an equilibrium-state model, and a direct elementary dynamics model are proposed to get the kinetic parameters for the overall reaction. In this paper, a conversion temperature T _C1 is defined as the temperature at which the standard molar Gibbs free energy change of the formation of the reaction complex is equal to zero, and it is found that when T ≫ T _C1, the direct elementary dynamics model with an inclusion of the tunneling correction of the second step reaction is applicable to calculate the overall reaction rate constants for this kind of reaction system. The reaction class of hydroxyl radical addition to alkenes is chosen as the objects of this study, five reactions are chosen as the representative for the reaction class, and their single-point energies are calculated using the method of CCSD(T)/CBS, and it is shown that the highest conversion temperature for the five reactions is 139.89 K, far below the usual initial low-temperature (550 K) oxidation chemistry of hydrocarbon fuels; therefore, the steady-state approximation method is applicable. All geometry optimizations are performed at the BH&HLYP/6-311+G(d,p) level, and the result shows that the geometric parameters in the reaction centers are conserved; hence, the isodesmic reaction method is applicable to this reaction class. To validate the accuracy of this scheme, a comparison of electronic energy difference at the BH&HLYP/6-311+G(d,p) level and the corrected electronic energy difference with the electronic energy difference at the CCSD(T)/CBS level is performed for the five representative reactions, and it is shown that the maximum absolute deviation of electronic energy difference can be reduced from 2.54 kcal·mol^-1 before correction to 0.58 kcal·mol^-1 after correction, indicating that the isodesmic reaction method is applicable for the accurate calculation of the kinetic parameters for large-size molecular systems with a negative activation energy reaction. The overall rate constants for 44 reactions of the reaction class of hydroxyl radical addition to alkenes are calculated using the transition-state theory in combination with the isodesmic correction scheme, and high-pressure limit rate rules for the reaction class are developed. In addition, the thermodynamic parameter is calculated and the results indicate that our dynamics model is applicable for our studied reaction class. A chemical kinetic modeling and sensitivity analysis using the calculated kinetic data is performed for the combustion of ethene, and the results indicate the studied reaction is important for the low-to-medium temperature combustion modeling of ethene.

Direct Kinetics and Product Measurement of Phenyl Radical + Ethylene

The phenyl + ethylene (C₆H₅ + C₂H₄) reaction network was explored experimentally and theoretically to understand the temperature dependence of the reaction kinetics and product distribution under various temperature and pressure conditions. The flash photolysis apparatus combining laser absorbance spectroscopy (LAS) and time-resolved molecular beam mass spectrometry (MBMS) was used to study reactions on the C₈H₉ potential energy surface (PES). In LAS experiments, 505.3 nm laser light selectively probed C₆H₅ decay, and we measured the total C₆H₅ consumption rate coefficients in the intermediate temperature region (400–800 K), which connects previous experiments performed in high-temperature (pyrolysis) and low-temperature (cavity-ring-down methods) regions. From the quantum chemistry calculations by Tokmakov and Lin using the G2M(RCC5)//B3LYP method, we constructed a kinetic model and estimated phenomenological pressure-dependent rate coefficients, k(T, P), with the Arkane package in the reaction mechanism generator. The MBMS experiments, performed at 600–800 K and 10–50 Torr, revealed three major product peaks: m/z = 105 (adducts, mostly 2-phenylethyl radical, but also 1-phenylethyl radical, ortho-ethyl phenyl radical, and a spiro-fused ring radical), 104 (styrene, co-product with a H atom), and 78 (benzene, co-product with C₂H₃ radical). Product branching ratios were predicted by the model and validated by experiments for the first time. At 600 K and 10 Torr, the yield ratio of the H-abstraction reaction (forming benzene + C₂H₃) is measured to be 1.1% and the H-loss channel (styrene + H) has a 2.5% yield ratio. The model predicts 1.0% for H-abstraction and 2.3% for H-loss, which is within the experimental error bars. The branching ratio and formation of styrene increase at high temperature due to the favored formally direct channel (1.0% at 600 K and 10 Torr, 5.8% at 800 K and 10 Torr in the model prediction) and the faster β-scission reactions of C₈H₉ isomers. The importance of pressure dependence in kinetics is verified by the increase in the yield of the stabilized adduct from radical addition from 80.2% (800 K, 10 Torr) to 88.9% (800 K, 50 Torr), at the expense of styrene + H. The pressure-dependent model developed in this work is well validated by the LAS and MBMS measurements and gives a complete picture of the C₆H₅ + C₂H₄ reaction.

Dual Fuel Reaction Mechanism 2.0 including NOx Formation and Laminar Flame Speed Calculations Using Methane/Propane/n-Heptane Fuel Blends

This study presents the further development of the TU Wien dual fuel mechanism, which was optimized for simulating ignition and combustion in a rapid compression expansion machine (RCEM) in dual fuel mode using diesel and natural gas at pressures higher than 60 bar at the start of injection. The mechanism is based on the Complete San Diego mechanism with n-heptane extension and was attuned to the RCEM measurements to achieve high agreement between experiments and simulation. This resulted in a specific application area. To obtain a mechanism for a wider parameter range, the Arrhenius parameter changes performed were analyzed and updated. Furthermore, the San Diego nitrogen sub-mechanism was added to consider NOx formation. The ignition delay time-reducing effect of propane addition to methane was closely examined and improved. To investigate the propagation of the flame front, the laminar flame speed of methane–air mixtures was simulated and compared with measured values from literature. Deviations at stoichiometric and fuel-rich conditions were found and by further mechanism optimization reduced significantly. To be able to justify the parameter changes performed, the resulting reaction rate coefficients were compared with data from the National Institute of Standards and Technology chemical kinetics database.

Oxidation Kinetics and Thermodynamics of Resonance-Stabilized Radicals: The Pent-1-en-3-yl + O2 Reaction

The kinetics and thermochemistry of the pent-1-en-3-yl radical reaction with molecular oxygen (CH2CHCHCH2CH3 + O2) has been studied by both experimental and computational methods. The bimolecular rate coefficient of the reaction was measured as a function of temperature (198-370 K) and pressure (0.2-4.5 Torr) using laser photolysis-photoionization mass-spectrometry. Quantum chemical calculations were used to explore the potential energy surface of the reaction, after which Rice-Ramsperger-Kassel-Marcus theory/master equation simulations were performed to investigate the reaction. The experimental data were used to adjust key parameters, such as well depths, in the master equation model within methodological uncertainties. The master equation simulations suggest that the formation rates of the two potential RO2 adducts are equal and that the reaction to QOOH is slower than for saturated hydrocarbons. The initial addition reaction, CH2CHCHCH2CH3 + O2, is found to be barrierless when accounting for multireference effects. This is in agreement with the current experimental data, as well as with past experimental data for the allyl + O2 reaction. Finally, we conducted numerical simulations of the pent-1-en-3-yl + O2 reaction system and observed significant amounts of penta-1,3-diene being formed under engine-relevant conditions.

From benzene to naphthalene: direct measurement of reactions and intermediates of phenyl radicals and acetylene

First-time measurement of time evolution of the main products and critical intermediates on phenyl HACA pathways with a validated pressure-dependent model.

Photodissociation dynamics of the simplest alkyl peroxy radicals, CH3OO and C2H5OO, at 248 nm

The photodissociation dynamics of the simplest alkyl peroxy radicals, methyl peroxy (CH₃OO) and ethyl peroxy (C₂H₅OO), are investigated using fast beam photofragment translational spectroscopy. A fast beam of CH₃OO^- or C₂H₅OO^- anions is photodetached to generate neutral radicals that are subsequently dissociated using 248 nm photons. The coincident detection of the photofragment positions and arrival times allows for the determination of mass, translational energy, and angular distributions for both two-body and three-body dissociation events. CH₃OO exhibits repulsive O loss resulting in the formation of O(¹D) + CH₃O with high translational energy release. Minor two-body channels leading to OH + CH₂O and CH₃O + O(³P) formation are also detected. In addition, small amounts of H + O(³P) + CH₂O are observed and attributed to O loss followed by CH₃O dissociation. C₂H₅OO exhibits more complex dissociation dynamics, in which O loss and OH loss occur in roughly equivalent amounts with O(¹D) formed as the dominant O atom electronic state via dissociation on a repulsive surface. Minor two-body channels leading to the formation of O₂ + C₂H₅ and HO₂ + C₂H₄ are also observed and attributed to a ground state dissociation pathway following internal conversion. Additionally, C₂H₅OO dissociation yields a three-body product channel, CH₃ + O(³P) + CH₂O, for which the proposed mechanism is repulsive O loss followed by the dissociation of C₂H₅O over a barrier. These results are compared to a recent study of tert-butyl peroxy (t-BuOO) in which 248 nm excitation results in three-body dissociation and ground state two-body dissociation but no O(¹D) production.

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.

Flow reactor studies of non-equilibrium plasma-assisted oxidation of n-alkanes

The oxidation of n-alkanes (C1-C7) has been studied with and without the effects of a nanosecond, non-equilibrium plasma discharge at 1 atm pressure from 420 to 1250 K. Experiments have been performed under nearly isothermal conditions in a flow reactor, where reactive mixtures are diluted in Ar to minimize temperature changes from chemical reactions. Sample extraction performed at the exit of the reactor captures product and intermediate species and stores them in a multi-position valve for subsequent identification and quantification using gas chromatography. By fixing the flow rate in the reactor and varying the temperature, reactivity maps for the oxidation of fuels are achieved. Considering all the fuels studied, fuel consumption under the effects of the plasma is shown to have been enhanced significantly, particularly for the low-temperature regime (T<800 K). In fact, multiple transitions in the rates of fuel consumption are observed depending on fuel with the emergence of a negative-temperature-coefficient regime. For all fuels, the temperature for the transition into the high-temperature chemistry is lowered as a consequence of the plasma being able to increase the rate of fuel consumption. Using a phenomenological interpretation of the intermediate species formed, it can be shown that the active particles produced from the plasma enhance alkyl radical formation at all temperatures and enable low-temperature chain branching for fuels C3 and greater. The significance of this result demonstrates that the plasma provides an opportunity for low-temperature chain branching to occur at reduced pressures, which is typically observed at elevated pressures in thermal induced systems.

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.

Study of the low temperature oxidation of propane

The low-temperature oxidation of propane was investigated using a jet-stirred reactor at atmospheric pressure and two methods of analysis: gas chromatography and synchrotron vacuum ultraviolet photoionization mass spectrometry (SVUV-PIMS) with direct sampling through a molecular jet. The second method allowed the identification of products, such as molecules with hydroperoxy functions, which are not stable enough to be detected by gas chromatography. Mole fractions of the reactants and reaction products were measured as a function of the temperature (530-730 K), with a particular attention to reaction products involved in the low temperature oxidation, such as cyclic ethers, aldehydes, alcohols, ketones, and hydroperoxides. A new model has been obtained from an automatically generated one, which was used as a starting point, with a large number of re-estimated thermochemical and kinetic data. The kinetic data of the most sensitive reactions, i.e., isomerizations of alkylperoxy radicals and the subsequent decompositions, have been calculated at the CBS-QB3 level of theory. The model allows a satisfactory prediction of the experimental data. A flow rate analysis has allowed highlighting the important reaction channels.

Improvement of the modeling of the low-temperature oxidation of n-butane: study of the primary reactions

This paper revisits the primary reactions involved in the oxidation of n-butane from low to intermediate temperatures (550-800 K) including the negative temperature coefficient (NTC) zone. A model that was automatically generated is used as a starting point and a large number of thermochemical and kinetic data are then re-estimated. The kinetic data of the isomerization of alkylperoxy radicals giving (•)QOOH radicals and the subsequent decomposition to give cyclic ethers has been calculated at the CBS-QB3 level of theory. The newly obtained model allows a satisfactory prediction of experimental data recently obtained in a jet-stirred reactor and in rapid compression machines. A considerable improvement of the prediction of the selectivity of cyclic ethers is especially obtained compared to previous models. Linear and global sensitivity analyses have been performed to better understand which reactions are of influence in the NTC zone.

Low-temperature combustion chemistry of biofuels: pathways in the initial low-temperature (550 K-750 K) oxidation chemistry of isopentanol

The branched C(5) alcohol isopentanol (3-methylbutan-1-ol) has shown promise as a potential biofuel both because of new advanced biochemical routes for its production and because of its combustion characteristics, in particular as a fuel for homogeneous-charge compression ignition (HCCI) or related strategies. In the present work, the fundamental autoignition chemistry of isopentanol is investigated by using the technique of pulsed-photolytic Cl-initiated oxidation and by analyzing the reacting mixture by time-resolved tunable synchrotron photoionization mass spectrometry in low-pressure (8 Torr) experiments in the 550-750 K temperature range. The mass-spectrometric experiments reveal a rich chemistry for the initial steps of isopentanol oxidation and give new insight into the low-temperature oxidation mechanism of medium-chain alcohols. Formation of isopentanal (3-methylbutanal) and unsaturated alcohols (including enols) associated with HO(2) production was observed. Cyclic ether channels are not observed, although such channels dominate OH formation in alkane oxidation. Rather, products are observed that correspond to formation of OH viaβ-C-C bond fission pathways of QOOH species derived from β- and γ-hydroxyisopentylperoxy (RO(2)) radicals. In these pathways, internal hydrogen abstraction in the RO(2)⇄ QOOH isomerization reaction takes place from either the -OH group or the C-H bond in α-position to the -OH group. These pathways should be broadly characteristic for longer-chain alcohol oxidation. Isomer-resolved branching ratios are deduced, showing evolution of the main products from 550 to 750 K, which can be qualitatively explained by the dominance of RO(2) chemistry at lower temperature and hydroxyisopentyl decomposition at higher temperature.

Towards cleaner combustion engines through groundbreaking detailed chemical kinetic models

In the context of limiting the environmental impact of transportation, this critical review discusses new directions which are being followed in the development of more predictive and more accurate detailed chemical kinetic models for the combustion of fuels. In the first part, the performance of current models, especially in terms of the prediction of pollutant formation, is evaluated. In the next parts, recent methods and ways to improve these models are described. An emphasis is given on the development of detailed models based on elementary reactions, on the production of the related thermochemical and kinetic parameters, and on the experimental techniques available to produce the data necessary to evaluate model predictions under well defined conditions (212 references).