New Insights into Low-Temperature Oxidation of Propane from Synchrotron Photoionization Mass Spectrometry and Multiscale Informatics Modeling

Modeling-Experiment-Theory Analysis of Reactions Initiated from Cl + Methyl Formate

Methyl formate (MF; CH3OCHO) is the smallest representative of esters, which are common components of biodiesel. The present study characterizes the thermal dissociation kinetics of the radicals formed by H atom abstraction from MF─CH3OCO and CH2OCHO─through a combination of modeling, experiment, and theory. For the experimental effort, excimer laser photolysis of Cl2 was used as a source of Cl atoms to initiate reactions with MF in the gas phase. Time-resolved species profiles of MF, Cl2, HCl, CO2, CH3, CH3Cl, CH2O, and CH2ClOCHO were measured and quantified using photoionization mass spectrometry at temperatures of 400-750 K and 10 Torr. The experimental data were simulated using a kinetic model, which was informed by ab initio-based theoretical kinetics calculations and included chlorine chemistry and secondary reactions of radical decomposition products. We calculated the rate coefficients for the H-abstraction reactions Cl + MF → HCl + CH3OCO (R1a) and Cl + MF → HCl + CH2OCHO (R1b): k1a,theory = 6.71 × 10-15·T1.14·exp(-606/T) cm3/molecule·s; k1b,theory = 4.67 × 10-18·T2.21·exp(-245/T) cm3/molecule·s over T = 200-2000 K. Electronic structure calculations indicate that the barriers to CH3OCO and CH2OCHO dissociation are 13.7 and 31.6 kcal/mol and lead to CH3 + CO2 (R3) and CH2O + HCO (R5), respectively. The master equation-based theoretical rate coefficients are k3,theory (P = ∞) = 2.94 × 109·T1.21·exp(-6209/T) s-1 and k5,theory (P = ∞) = 8.45 × 108·T1.39·exp(-15132/T) s-1 over T = 300-1500 K. The calculated branching fractions into R1a and R1b and the rate coefficient for R5 were validated by modeling of the experimental species time profiles and found to be in excellent agreement with theory. Additionally, we found that the bimolecular reactions CH2OCHO + Cl, CH2OCHO + Cl2, and CH3 + Cl2 were critical to accurately model the experimental data and constrain the kinetics of MF-radicals. Inclusion of the kinetic parameters determined in this study showed a significant impact on combustion simulations of larger methyl esters, which are considered as biodiesel surrogates.

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.

Direct time-resolved detection and quantification of key reactive intermediates in diethyl ether oxidation at T = 450-600 K

High-pressure multiplexed photoionization mass spectrometry (MPIMS) with tunable vacuum ultraviolet (VUV) ionization radiation from the Lawrence Berkeley Labs Advanced Light Source is used to investigate the oxidation of diethyl ether (DEE). Kinetics and photoionization (PI) spectra are simultaneously measured for the species formed. Several stable products from DEE oxidation are identified and quantified using reference PI cross-sections. In addition, we directly detect and quantify three key chemical intermediates: peroxy (ROO˙), hydroperoxyalkyl peroxy (˙OOQOOH), and ketohydroperoxide (HOOP[double bond, length as m-dash]O, KHP). These intermediates undergo dissociative ionization (DI) into smaller fragments, making their identification by mass spectrometry challenging. With the aid of quantum chemical calculations, we identify the DI channels of these key chemical species and quantify their time-resolved concentrations from the overall carbon atom balance at T = 450 K and P = 7500 torr. This allows the determination of the absolute PI cross-sections of ROO˙, ˙OOQOOH, and KHP into each DI channel directly from experiment. The PI cross-sections in turn enable the quantification of ROO˙, ˙OOQOOH, and KHP from DEE oxidation over a range of experimental conditions that reveal the effects of pressure, O₂ concentration, and temperature on the competition among radical decomposition and second O₂ addition pathways.

Identifying isomers of peroxy radicals in the gas phase: 1-C3H7O2vs. 2-C3H7O2

The two isomers of propylperoxy radical 1-C₃H₇O₂ and 2-C₃H₇O₂, together with their rotamers, are individually identified and assigned.

Shock Tube Laser Schlieren Study of the Pyrolysis of Isopropyl Nitrate

The decomposition of isopropyl nitrate was measured behind incident shock waves using laser schlieren densitometry in a diaphragmless shock tube. Experiments were conducted over the temperature range of 700-1000 K and at pressures of 71, 126, and 240 Torr. Electronic structure theory and RRKM Master Equation methods were used to predict the decomposition kinetics. RRKM/ME parameters were optimized against the experimental data to provide an accurate prediction over a broader range of conditions. The initial decomposition i-C₃H₇ONO₂ ⇌ i-C₃H₇O + NO₂ has a high-pressure limit rate coefficient of 5.70 × 10²²T^-1.80 exp[-21287.5/T] s^-1. A new chemical kinetic mechanism was developed to model the chemistry after the initial dissociation. A new shock tube module was developed for Cantera, which allows for arbitrarily large mechanisms in the simulation of laser schlieren experiments. The present work is in good agreement with previous experimental studies.

Decomposition kinetics for HONO and HNO2

This work presents a detailed investigation into the isomerization and decomposition of HONO and HNO₂.

Kinetics of the a-C3H5 + O2 reaction, investigated by photoionization using synchrotron radiation

The kinetics of the combustion-relevant reaction of the allyl radical, a-C3H5, with molecular oxygen has been studied in a flow tube reactor at the vacuum ultraviolet (VUV) beamline of the Swiss Light Source storage ring, using the CRF-PEPICO (Combustion Reactions Followed by Photoelectron Photoion Coincidence Spectroscopy) setup. The ability to measure threshold photoelectron spectra enables a background-free detection of reactive species as well as an isomer-specific analysis of reaction products. Allyl was generated by direct photodissociation of allyl iodide at 266 nm and 213 nm and indirectly by the reaction of propene with Cl atoms, which were generated by photolysis from oxalyl chloride at 266 nm. Experiments were conducted at room temperature at low pressures between 0.8 and 3 mbar using Ar as the buffer gas and with excess O2 to maintain nearly pseudo-first-order reaction conditions. Whereas allyl was detected by photoionisation using synchrotron radiation, the main reaction product allyl peroxy was not observed due to dissociative ionisation of this weakly bound species. From the concentration-time profiles of the allyl signal, second-order rate constants between 1.35 × 1011 cm3 mol-1 s-1 at 0.8 mbar and 1.75 × 1011 cm3 mol-1 s-1 at 3 mbar were determined. The rates obtained for the different allyl radical generation schemes agree well with each other, but are about a factor of 2 higher than the ones reported previously using He as a buffer gas. The discrepancy is partly attributed to the higher collision efficiency of Ar causing a varying fall-off behavior. When allyl is produced by the reaction of propene with Cl atom, an unexpected product is observed at m/z = 68, which was identified as 1,3-butadienal in the threshold photoelectron spectrum. It is formed in a secondary reaction of allyl with the OCCl radical, which is generated in the 266 nm photolysis of oxalyl chloride.

The reaction of Criegee intermediate CH2OO with water dimer: primary products and atmospheric impact

The rapid reaction of the smallest Criegee intermediate, CH₂OO, with water dimers is the dominant removal mechanism for CH₂OO in the Earth's atmosphere, but its products are not well understood. This reaction was recently suggested as a significant source of the most abundant tropospheric organic acid, formic acid (HCOOH), which is consistently underpredicted by atmospheric models. However, using time-resolved measurements of reaction kinetics by UV absorption and product analysis by photoionization mass spectrometry, we show that the primary products of this reaction are formaldehyde and hydroxymethyl hydroperoxide (HMHP), with direct HCOOH yields of less than 10%. Incorporating our results into a global chemistry-transport model further reduces HCOOH levels by 10-90%, relative to previous modeling assumptions, which indicates that the reaction CH₂OO + water dimer by itself cannot resolve the discrepancy between the measured and predicted HCOOH levels.