Theoretical kinetics of HO2 + C5H5: A missing piece in cyclopentadienyl radical oxidation reactions

The resonantly-stabilized cyclopentadienyl radical (C₅H₅) is a key species in the combustion and molecular growth kinetics of mono and poly-aromatic hydrocarbons (M/PAHs). At intermediate-to-low temperatures, the C₅H₅ reaction with the hydroperoxyl radical (HO₂) strongly impacts the competition between oxidation to smaller products and growth to PAHs, precursors of soot. However, literature estimates for the HO₂ + C₅H₅ reaction rate are inaccurate and inconsistent with recent theoretical calculations, thus generating discrepancies in global combustion kinetic models. In this work, we perform state-of-the-art theoretical calculations for the HO₂ + C₅H₅ reaction including variable reaction coordinate transition state theory for barrierless channels, accurate thermochemistry, and multi-well master equation (ME) simulations. Contrary to previous studies, we predict that OH + 1,3-C₅H₅O is the main reaction channel. The new rate constants are introduced in two literature kinetic models exploiting our recently developed ME based lumping methodology and used to perform kinetic simulations of experimental data of MAHs oxidation. It is found that the resonantly-stabilized 1,3-C₅H₅O radical is the main C₅H₅O isomer, accumulating in relevant concentration in the system, and that the adopted lumping procedure is fully consistent with results obtained with detailed kinetics. The reactivity of C₅H₅O with OH and O₂ radicals is included in the kinetic mechanisms based on analogy rules. As a result, C₅H₅O mostly reacts with O₂ producing smaller C₃/C₄ species and large amounts of C₅H₄O, suggesting that further investigations of the reactivity of both C₅H₅O and C₅H₄O with oxygenated radicals is necessary. Overall, this work presents new reliable rate constants for the HO₂ + C₅H₅ reaction and provides indications for future investigations of relevant reactions in the sub-mechanisms of cyclopentadiene and MAH oxidation.     

Theoretical rate coefficients for the reaction of methyl radical with hydroperoxyl radical and for methylhydroperoxide decomposition

The kinetics of the CH₃ + HO₂ bimolecular reaction and the thermal decomposition of CH₃OOH are studied theoretically. Direct variable reaction coordinate transition state theory (VRC-TST), coupled with high level multireference electronic structure calculations, is used to compute capture rates for the CH₃ + HO₂ reaction and to characterize the transition state of the barrierless CH₃O + OH product channel. The CH₂O + H₂O product channel and the CH₃ + HO₂ → CH₄ + O₂ reaction are treated using variational transition state theory and the harmonic oscillator and rigid rotor approximations. Pressure dependence and product branching in the bimolecular and decomposition reactions are modeled using master equation simulations. The predicted rate coefficients for the major products channels of the bimolecular reaction, CH₃O + OH and CH₄ + O₂, are found to be in excellent agreement with values obtained in two recent modeling studies. The present calculations are also used to obtain rate coefficients for the CH₃O + OH association/decomposition reaction.     

A quantitative explanation for the apparent anomalous temperature dependence of OH + HO2 = H2O + O2 through multi-scale modeling

Kinetic models for complex chemical mechanisms are comprised of tens to thousands of reactions with rate constants informed by data from a wide variety of sources – rate constant measurements, global combustion experiments, and theoretical kinetics calculations. In order to integrate information from distinct data types in a self-consistent manner, a framework for combustion model development is presented that encapsulates behavior across a wide range of chemically relevant scales from fundamental molecular interactions to global combustion phenomena. The resulting kinetic model consists of a set of theoretical kinetics parameters (with constrained uncertainties), which are related through kinetics calculations to temperature/pressure/bath-gas-dependent rate constants (with propagated uncertainties), which in turn are related through physical models to combustion behavior (with propagated uncertainties). Direct incorporation of theory in combustion model development is expected to yield more reliable extrapolation of limited data to conditions outside the validation set, which is particularly useful for extrapolating to engine-relevant conditions where relatively limited data are available. Several key features of the approach are demonstrated for the H₂O₂ decomposition mechanism, where a number of its constituent reactions continue to have large uncertainties in their temperature and pressure dependence despite their relevance to high-pressure, low-temperature combustion of a variety of fuels. Here, we use the approach to provide a quantitative explanation for the apparent anomalous temperature dependence of OH + HO₂ = H₂O + O₂ – in a manner consistent with experimental data from the entire temperature range and ab initio transition-state theory within their associated uncertainties. Interestingly, we do find a rate minimum near 1200 K, although the temperature dependence is substantially less pronounced than previously suggested.     

Methanol oxidation up to 100 atm in a supercritical pressure jet-stirred reactor

Methanol (CH₃OH) has attracted considerable attention as a renewable fuel or fuel additive with low greenhouse gas emissions. Methanol oxidation was studied using a recently developed supercritical pressure jet-stirred reactor (SP-JSR) at pressures of 10 and 100 atm, at temperatures from 550 to 950 K, and at equivalence ratios of 0.1, 1.0, and 9.0 in experiments and simulations. The experimental results show that the onset temperature of CH₃OH oxidation at 100 atm is around 700 K, which is more than 100 K lower than the onset at 10 atm and this trend cannot be predicted by the existing kinetics models. Furthermore, a negative temperature coefficient (NTC) behavior was clearly observed at 100 atm at fuel rich conditions for methanol for the first time. To understand the observed temperature shift in the reactivity and the NTC effect, we updated some key elementary reaction rates of relevance to high pressure CH₃OH oxidation from the literature and added some new low-temperature reaction pathways such as CH₂O + HO₂ = HOCH₂O₂ (RO₂), RO₂ + RO₂ = HOCH₂O (RO) + HOCH₂O (RO) + O₂, and CH₃OH + RO₂ = CH₂OH + HOCH₂O₂H (ROOH). Although the model with these updates improves the prediction somewhat for the experimental data at 100 atm and reproduces well high-temperature ignition delay times and laminar flame speed data in the literature, discrepancies still exist for some aspects of the 100 atm low-temperature oxidation data. In addition, it was found that the pressure-dependent HO₂ chemistry shifts to lower temperature as the pressure increases such that the NTC effect at fuel-lean conditions is suppressed. Therefore, as shown in the experiments, the NTC phenomenon was only observed at the fuel-rich condition where fuel radicals are abundant and the HO₂ chemistry at high pressure is weakened by the lack of oxygen resulting in comparatively little HO₂ formation.     

Adjoint sensitivity analysis of kinetic,thermochemical, and transport data of nitrogen and ammonia chemistry

Many studies apply sensitivity analysis to explore the impact of reaction kinetic parameters on model predictions. The importance of thermochemical and transport data is often assumed to be relatively low. While this is true for specific combustion properties of hydrocarbons, the role of thermochemical and transport data in combustion processes of nitrogen-containing molecules remains to be investigated. Thus, this work applies adjoint sensitivity analysis to the complete set of parameters in combustion models, i.e., kinetics, thermodynamics, and transport data. This integral approach increases the number of parameters considered in the sensitivity analysis drastically. Compared to forward sensitivity analysis, the adjoint approach is very efficient for a large number of parameters, and analysis with several thousand parameters can be performed in seconds. Nitrogen oxide formation in methane/air flames and laminar burning velocities of ammonia/air flames are considered as prediction targets. Sensitivity analysis results for kinetic, thermochemical, and transport data are compared by jointly considering all appearing parameter uncertainties. The comparison reveals that, due to their importance for the equilibrium constants of elementary reactions, the optimization potential of thermodynamic properties is often similarly high as that of the kinetics parameters. Transport parameters are found to be of the lowest priority for the model development due to their low uncertainties, even though high sensitivities are determined for several of them. More specifically, the analysis for the laminar burning velocities of ammonia/air flames reveals a high optimization potential for parameters in the N₂-amine chemistry, including the molar heat capacities of N₂H₂, N₂H₃, and NH. Interestingly, analyses with different mechanisms reveal strongly diverging results, especially regarding the importance of reactions with OH, which is uncommon when considering the combustion of hydrocarbons.     

Kinetic modeling of the benzyl + HO2 reaction

Benzyl is a resonantly stabilized radical that commonly occurs as an intermediate in the combustion of aromatic compounds. The bimolecular reaction of benzyl with HO₂ is important in the oxidation of toluene, especially at low to moderate temperatures, where unimolecular decomposition of the benzyl radical is slow. We show that the addition of HO₂ to the methylene site in benzyl produces a vibrationally excited benzylhydroperoxide adduct, with over 60 kcal mol⁻¹ (251 kJ mol⁻¹) of excess energy above the ground state. RRKM simulations are performed on the benzyl + HO₂ reaction, using thermochemical and kinetic parameters obtained from ab initio calculations, with variational transition state theory (VTST) for treatment of barrierless radical + radical reaction kinetics. Our results reveal that the benzyl + HO₂ reaction proceeds predominantly to the benzoxyl radical + OH at temperatures of around 800 K and above, with the production of stabilized benzylhydroperoxide molecules dominating at lower temperatures. The heat of formation of the benzyl radical is calculated as 52.5 kcal mol⁻¹ (219.7 kJ mol⁻¹) at the G3B3 level of theory, in relative agreement with other recent determinations of this value.     

From quantum chemistry to dissociation kinetics: what we need to know †

The relationship between rate constants for dissociation and the reverse association reactions and their potential energy surfaces is illustrated. The reaction systems e^? + SF₆ ? SF₆ ^? →SF₅ ^? + F, H + CH₃ ?CH₄, 2 CF₂ ? C₂F₄, H + O₂ →HO₂, HO + O ?HO₂ ? H + O₂, and C + HO →CHO are chosen as representative examples. The necessity to know precise thermochemical data is emphasised. The interplay between attractive and anisotropic components of the potentials influences the rate constants. Spin–orbit and electronic–rotational coupling in reactions between electronic open-shell radicals so far generally has been neglected, but is shown to have a marked influence on low temperature rate constants.     

Prompt NO formation in flames: The influence of NCN thermochemistry

The influence of the route via the NCN radical on NO formation in flames was examined from a thermochemistry and reaction kinetics perspective. A detailed analysis of available experimental and theoretical thermochemical data combined with an Active Thermochemical Tables analysis suggests a heat of formation of 457.8 ± 2.0 kJ/mol for NCN, consistent with carefully executed theoretical work of Harding et al. (2008) [5]. This value is significantly different from other previously reported experimental and theoretical values. A combination of an extensively validated comprehensive hydrocarbon oxidation model extended by the GDFkin3.0_NCN-NO_x sub-mechanism reproduced NCN and NO mole fraction profiles in a recently characterized fuel-rich methane flame only when heat of formation values in the range of 445–453 kJ/mol are applied. Sensitivity analysis revealed that the sensitivities of contributing steps to NO and NCN formation are strongly dependent on the absolute value of the heat of formation of NCN being used. In all flames under study the applied NCN thermochemistry highly influences simulated NO and NCN mole fractions. The results of this work illustrate the thermochemistry constraints in the context of NCN chemistry which have to be taken into account for improving model predictions of NO concentrations in flames.     

Formation pathways of HO2 and OH changing as a function of temperature in photolytically initiated oxidation of dimethyl ether

The time resolved product formation in oxidation of dimethyl ether (DME) has been studied between 298-625 K and 20-90 torr total pressure. Near-infrared frequency modulation spectroscopy (FMS) with Herriott type multi pass optics and UV absorption spectroscopy (UV) were conducted in the same cell. The reaction was initiated by pulsed photolysis in a mixture of Cl₂, O₂, and DME via CH₃OCH₂ radical formation. The reaction process was investigated through FMS measurement of HO₂ and OH, and UV measurement of CH₃OCH₂O₂. The yields of HO₂ and OH are obtained by comparison with reference mixtures, Cl₂, O₂, and CH₃OH for HO₂, and Cl₂, O₂, CH₃OH, and NO for OH, which convert 100% of initial Cl to HO₂ and OH. The CH₃OCH₂O₂ yield is also obtained. It was found that the HO₂ yield increases sharply over 500 K mainly with a longer time constant than that of R + O₂ reaction, while a prompt component exists throughout the temperature range at a few percent yield. OH was found to be produced promptly at a yield considerably larger than that known for the simplest alkanes. The CH₃OCH₂O₂ profile has a prompt rise followed by a gradual decay whose rate is consistent with the slow HO₂ formation. The species profiles were successfully predicted with a model constructed by modifying the existing one to suit the reduced pressure condition. After modification, it was inferred that the HO₂ formation over 500 K is secondary from HCHO + OH and HCO + O₂ and a part of HCO is formed directly from the O₂ adduct, whereas the HO₂ formation below 500 K is governed by CH₃OCH₂O₂ chemistry. The HCO forming pathway via isomerization-decomposition of the O₂ adduct, which was not included in the former models, was supported by our quantum-chemical calculations.     

Hydrogen atom abstraction from Piperazine by hydroxyl radical: a theoretical investigation

Dual level of quantum mechanical calculations have been carried out for hydrogen abstraction from Piperazine [HN(CH₂CH₂)₂NH] initiated by OH radical. Geometry optimisation and frequency calculations of all species involved in the titled reaction have been performed at M06-2X/6-31+G(d,p) level of theory. For the accuracy in the thermochemistry and kinetics data, single-point energy calculations have been further carried out at coupled cluster CCSD(T) method along with 6-311G(d,p) basis set. An energy profile diagram for the reaction has been plotted along with pre-reactive and post-reactive complexes at entrance and exit channels. Intrinsic reaction coordinates (IRCs) calculations have been performed for identification of real transition states that connect it via reactant to product. Our result shows that the H-atom abstraction takes place from the C–H position of Piperazine. The rate constant is calculated using canonical transition state theory (CTST) is found to be 2.86 × 10^?10 cm³ molecule^?1 s^?1 which is in good agreement with the reported experimental rate constant (2.38 ± 0.28) × 10^?10 cm³ molecule^?1 s^?1 at 298 K. We have also reported rate constant for the temperature range 300–500 K. Using group-balance isodesmic reaction, the standard enthalpies of formation for Piperazine and product radicals generated by hydrogen abstraction are reported. The branching ratios for both reaction channel (i.e. H-abstraction from –CH₂ and –NH position of Piperazine) are found to be 93% and 7%, respectively. The calculated atmospheric life time of Piperazine is found to be 0.97 hour.     

On the formation and decomposition of C7H8

The kinetics of reactions on the C₇H₈ surface were studied with state-of-the-art ab initio transition state theory (TST) and master equation methodologies. A priori predictions of the capture rate for C₆H₅ + CH₃ and for C₇H₇ + H are obtained from direct variable reaction coordinate TST simulations. These simulations employ small basis set CASPT2 interaction energies coupled with one-dimensional reaction path corrections based on higher level simulations for related reactions. For the C₇H₇ + H reaction, predictions are obtained for both the total rate and for the branching between toluene, o-isotoluene and p-isotoluene. A mapping of the low energy pathways for isomerization from these three C₇H₈ isomers identifies a number of processes with barriers at or below the dissociation threshold. Nevertheless, at combustion temperatures the dissociation rates are predicted to exceed the isomerization rates, and it is reasonable to treat the kinetics of each isomer as a simple single well association/dissociation equilibrium. Master equation simulations yield predictions for the temperature and pressure dependence of each of the recombination and dissociation processes, as well as for the C₇H₇ + H → C₆H₅ + CH₃ bimolecular reaction. These simulations implement collisional energy transfer probabilities based on the work of Luther and co-workers. The theoretical predictions are found to be in satisfactory agreement with the available experimental data for the photodissociation of toluene, the temperature and pressure dependent dissociation of toluene, and the reaction of benzyl radical with H. For the C₆H₅ + CH₃ recombination, the theoretical predictions exceed the experimental measurements of Lin and coworkers by a factor of 2 or more for all temperatures.     

Theoretical investigation into the low-temperature oxidation of ethylbenzene

The classical topic on the oxidation of alkylbenzene has been revisited via performing accurate theoretical calculations to address the salient features for the initial oxidation of ethylbenzene. Potential energy surfaces are mapped out for all possible reactions in the systems of (1-phenylethyl + O₂ and 2-phenylethyl + O₂). Reaction rate constants at the high-pressure limit are calculated for all possible reactions in these two systems. Direct H abstraction from 1-phenylethyl radical by oxygen molecule appears to be an important route for the formation of styrene from the oxidation of ethylbenzene. Concerted elimination of HO₂ is predicted to contribute significantly the production of styrene from system of 2-phenylethyl + O₂; especially at the atmospheric pressure and intermediate temperatures. Formation of the other major experimental product, benzaldehyde, is attributed to the unimolecular decomposition of C₆H₅CH₂(O)CH₃ rather than to unimolecular isomerisation of the two initial peroxy adducts. Kinetic and mechanistic data presented herein are instrumental for better understanding of the oxidative decomposition of ethylbenzene, i.e., major constituents of commonly formulated fuel surrogates.     

Acetaldehyde oxidation at elevated pressure

A detailed chemical kinetic model for oxidation of CH₃CHO at intermediate to high temperature and elevated pressure has been developed and evaluated by comparing predictions to novel high-pressure flow reactor experiments as well as shock tube ignition delay measurements and jet-stirred reactor data from literature. The flow reactor experiments were conducted with a slightly lean CH₃CHO/O₂ mixture highly diluted in N₂ at 600–900 K and pressures of 25 and 100 bar. At the highest pressure, the oxidation of CH₃CHO was in the NTC regime, controlled to a large extent by the thermal stability and reactions of peroxide species such as HO₂, CH₃OO, and CH₃C(O)OO. Model predictions were generally in good agreement with the experimental data, even though the predicted temperature for onset of reaction was overpredicted at 100 bar. This discrepancy was attributed mainly to uncertainties in the CH₃C(O)OO reaction subset. Predictions of ignition delays in shock tubes and species profiles in JSR experiments were also satisfactory. At temperatures above the NTC regime, acetaldehyde ignition and oxidation is affected mainly by the competition between dissociation of CH₃CHO and reaction with the radical pool, and by reactions in the methane subset.     

H-abstractions by O2, NO2, NH2, and HO2 from H2NO: Theoretical study and implications for ammonia low-temperature kinetics

Unraveling the low-temperature chemistry of ammonia is still an open challenge in combustion kinetics, yet of primary importance because of the novel combustion concepts operating in these conditions, as well as of the rising interest on ammonia as an energy carrier. In this work, a fundamental investigation of the H-abstraction reactions from H₂NO by O₂, NO₂, NH₂, and HO₂ was performed. These reactions, which belong to the radical-radical abstraction class, associate a high sensitivity to the key low temperature ammonia combustion parameters, to a high uncertainty in rate constant values. Theoretically, the investigation of reactions belonging to this class is complicated by their intrinsic multireference nature. To address this issue, a structured theoretical methodology that relies heavily on the use of CASPT2 calculations was devised. The predicted rate constants highlighted significant deviations from the rates commonly adopted in the state-of-the-art mechanisms, most often based on analogies and estimations. In order to understand their impact on ammonia low-temperature kinetics, the obtained rates were integrated into a kinetic model, which was used to investigate ammonia oxidation and ignition at low-temperature and oxygen-rich conditions. It was found that O₂ and NH₂ play the major role, as abstractors, in regulating ammonia oxidation and ignition. In particular, ignition delay time predictions proved extremely sensitive to the adopted rates: modifying each of them within their theoretical uncertainty caused deviations by even an order of magnitude, and totally changed the predictive features of the mechanism. The kinetic analysis highlighted then the need of a targeted optimization of the critical rates, downstream of the present work and within their uncertainty boundaries, to further refine the mechanism capability over a wide range of operating conditions.     

Vibrational non-equilibrium in the hydrogen–oxygen reaction. Comparison with experiment

A theoretical model is proposed for the chemical and vibrational kinetics of hydrogen oxidation based on consistent accounting of the vibrational non-equilibrium of the HO₂ radical that forms as a result of the bimolecular recombination H+O₂ → HO₂. In the proposed model, the chain branching H+O₂ = O+OH and inhibiting H+O₂+M = HO₂+M formal reactions are treated (in the terms of elementary processes) as a single multi-channel process of forming, intramolecular energy redistribution between modes, relaxation, and unimolecular decay of the comparatively long-lived vibrationally excited HO₂ radical, which is able to react and exchange energy with the other components of the mixture. The model takes into account the vibrational non-equilibrium of the starting (primary) H₂ and O₂ molecules, as well as the most important molecular intermediates HO₂, OH, O₂(¹Δ), and the main reaction product H₂O. It is shown that the hydrogen–oxygen reaction proceeds in the absence of vibrational equilibrium, and the vibrationally excited HO₂(v) radical acts as a key intermediate in a fundamentally important chain branching process and in the generation of electronically excited species O₂(¹Δ), O(¹D), and OH(²Σ⁺). The calculated results are compared with the shock tube experimental data for strongly diluted H₂–O₂ mixtures at 1000 < T < 2500 K, 0.5 < p < 4 atm. It is demonstrated that this approach is promising from the standpoint of reconciling the predictions of the theoretical model with experimental data obtained by different authors for various compositions and conditions using different methods. For T < 1500 K, the nature of the hydrogen–oxygen reaction is especially non-equilibrium, and the vibrational non-equilibrium of the HO₂ radical is the essence of this process. The quantitative estimation of the vibrational relaxation characteristic time of the HO₂ radical in its collisions with H₂ molecules has been obtained as a result of the comparison of different experimental data on induction time measurements with the relevant calculations.     

Computational studies on the thermal decomposition of the CH2FOCHFO radical

Theoretical studies have been carried out on the kinetics and thermochemistry of the thermal decomposition of the CH₂FOCHFO radical formed during the photo-oxidation of CH₂FOCH₂F (HFE-152E) using the dual-level method of obtaining the optimised structure at DFT(M06-2X)/6-311++G(d,p) followed by a single-point energy calculation at the G3 level of theory. The rate constant for different reaction channels involved during the decomposition processes of CH₂FOCHFO is evaluated at 298 K and 1 atm using canonical transition-state theory. The results point out that the C–H bond scission is the dominant path involving an energy barrier of 9.5 kcal mol^?1 determined at the G3 level of theory. A potential energy diagram is constructed and the results are compared with the data available from the literature for a structurally similar molecule.     

On the mechanism of inhibition of a hydrogen flame by propane

Numerical simulations demonstrated that small additives of propane to rich hydrogen-air mixtures suppress the formation of HO₂ and OH in the low-temperature region of the flame zone, thereby causing a substantial decrease in the laminar flame speed. In the low-and high-temperature regions, propane interacts predominantly with OH and H, respectively. In the flame zone, propane is completely converted to CO, CO₂, CH₄, C₂H₂, H₂, and H₂O, being oxidized concurrently with hydrogen at that.     

Identification of homogeneous chemical kinetic regimes of methane-air ignition

Sensitivity analysis results for ignition delay time (IDT) may be very different depending on the initial temperature, pressure and equivalence ratio φ, but similar in some regions of these variables. This phenomenon was investigated systematically by carrying out ignition simulations and local sensitivity calculations of methane−air mixtures using the Aramco-II-2016 mechanism at 14,417 combinations of initial temperature (changed between 500 and 3000 K), initial pressure (0.05−500 atm) and φ (0.05−8.0) values. The cluster analysis of the sensitivity vectors identified five large kinetically homogeneous regions. Each region has well defined borders in the (T, p φ) space and can be characterized by different sets of important reactions. The related kinetic scheme is very different in each region. Regions 1 and 2 are dominated by catalytic cycles based on species CH₃O₂/CH₃O₂H and HO₂/H₂O₂/CH₃O, respectively. In regions 3, 4, and 5 the H atoms are converted to CH₃ in an identical chain branching sequence, but the back conversion is via three different routes. Literature experimental data on the IDTs of methane−air mixtures were sorted according to these five regions. Regions 1 to 5 contain 214, 328, 3, 0, and 237 experimental data points, respectively. In regions 1, 2 and 5 the data points are well reproduced by the Aramco-II-2016 mechanism, but little or no experimental information is available about kinetic regions 3 and 4. Further experimental exploration of the ignition of methane−air mixtures may aim the study of these regions. A similar approach can be used for the characterization of other combustion systems and sorting the related experimental data.     

Experimental and detailed kinetic modeling study of the high pressure oxidation of methanol sensitized by nitric oxide and nitrogen dioxide

The perturbation of the combustion by NO_x is important in several practical systems (recent NO_x-reduction strategies, combustion with exhaust-gas recirculation in diesel and HCCI engines and for mild combustion). New experimental results were obtained for the oxidation of methanol in absence and in presence of NO or NO₂ in a fused silica jet-stirred reactor operating at 10 atm, over the temperature range 700-1100 K. Probe sampling followed by on-line FTIR analyses and off-line GC-TCD/FID analyses permitted to measure the concentration profiles of the reactants, stable intermediates and the final products. A detailed chemical kinetic modeling of the present experiments was performed. An overall good agreement between the present data and this modeling was obtained. The oxidation of methanol is significantly sensitized by NO₂, whereas the effect of NO is more limited. According to the proposed model, the mutual sensitization of the oxidation of methanol and NO proceeds through the NO to NO₂ conversion by HO₂. The increased production of OH resulting from the oxidation of NO by HO₂ promotes the oxidation of the fuel. A simplified reaction scheme can be proposed for the NO-seeded oxidation of methanol: NO + HO₂ ⇒ NO₂ + OH followed by OH + CH₃OH ⇒ H₂O + CH₂OH and CH₃O. The enhanced oxidation of methanol by addition of NO₂ is also due to additional OH production through: NO₂ + HO₂ ⇒ HONO + O₂, NO₂ + H ⇒ NO + OH and HONO ⇒ NO + OH followed by OH + CH₃OH ⇒ CH₂OH and CH₃O. The further reactions CH₂OH + O₂ ⇒ CH₂O + HO₂; CH₃O ⇒ CH₂O + H; CH₂O + OH ⇒ HCO; HCO + O₂ ⇒ HO₂ and H + O₂ ⇒ HO₂ complete the sequence whether NO or NO₂ is added.     

An experimental and theoretical kinetic study of t-butyl radical reaction with molecular oxygen

In this work, a unifying picture of the kinetics of the t-C₄H₉ + O₂ reaction is presented by combining the current and previous experimental results with theory. Direct, time-resolved experiments were performed over a wide temperature range (200–500 K) at low pressures (0.3–6 Torr) using a photoionization mass spectrometry method. The kinetic measurements of the t-butyl + O₂ reaction were initiated by laser photolysis of pinacolone at 193 nm or t-butyl bromide at 248 nm to produce t-C₄H₉ radicals. Energies calculated by quantum chemistry at the CCSD(T)/CBS and CASPT2/CBS levels of theory were used in master equation simulations of the kinetics of the t-C₄H₉ + O₂ reaction. The calculations successfully reproduce the pressure and temperature dependencies of both the current low-pressure experiments and literature kinetic data at about atmospheric pressure as well as the literature kinetic data for the overwhelmingly most important bimolecular reaction channel, t-C₄H₉ + O₂ → i-C₄H₈ + HO₂ in the intermediate temperature range. The experimentally constrained master equation model was utilized to simulate the t-C₄H₉ + O₂ reaction kinetics over wide range of conditions. The results of these simulations are provided in ChemKin compatible PLOG format for later use.