Publications

Methyl formate (MF; CH₃OCHO) 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—CH₃OCO and CH₂OCHO—through a combination of modeling, experiment, and theory. For the experimental effort, excimer laser photolysis of Cl₂ was used as a source of Cl atoms to initiate reactions with MF in the gas phase. Time-resolved species profiles of MF, Cl₂, HCl, CO₂, CH₃, CH₃Cl, CH₂O, and CH₂ClOCHO 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. Here, we calculated the rate coefficients for the H-abstraction reactions Cl + MF → HCl + CH₃OCO (R1a) and Cl + MF → HCl + CH₂OCHO (R1b): k_1a,theory = 6.71 × 10^–15·T^1.14·exp(—606/T) cm³/molecule·s; k_1b,theory = 4.67 × 10^–18·T^2.21·exp(—245/T) cm³/molecule·s over T = 200–2000 K. Electronic structure calculations indicate that the barriers to CH₃OCO and CH₂OCHO dissociation are 13.7 and 31.6 kcal/mol and lead to CH₃ + CO₂ (R3) and CH₂O + HCO (R5), respectively. The master equation-based theoretical rate coefficients are k_3,theory (P = ∞) = 2.94 × 10⁹·T^1.21·exp(—6209/T) s^–1 and k_5,theory (P = ∞) = 8.45 × 10⁸·T^1.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 CH₂OCHO + Cl, CH₂OCHO + Cl₂, and CH₃ + Cl₂ 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.

Methyl-ethyl-substituted Criegee intermediate (MECI) is a four-carbon carbonyl oxide that is formed in the ozonolysis of some asymmetric alkenes. MECI is structurally similar to the isoprene-derived methyl vinyl ketone oxide (MVK-oxide) but lacks resonance stabilization, making it a promising candidate to help us unravel the effects of size, structure, and resonance stabilization that influence the reactivity of atmospherically important, highly functionalized Criegee intermediates. We present experimental and theoretical results from the first bimolecular study of MECI in its reaction with SO2, a reaction that shows significant sensitivity to the Criegee intermediate structure. Using multiplexed photoionization mass spectrometry, we obtain a rate coefficient of (1.3 ± 0.3) × 10-10 cm3 s-1 (95% confidence limits, 298 K, 10 Torr) and demonstrate the formation of SO3 under our experimental conditions. Through high-level theory, we explore the effect of Criegee intermediate structure on the minimum energy pathways for their reactions with SO2 and obtain modified Arrhenius fits to our predictions for the reaction of both syn and anti conformers of MECI with SO2 (ksyn = 4.42 × 1011 T-7.80exp(−1401/T) cm3 s-1 and kanti = 1.26 × 1011 T-7.55exp(−1397/T) cm3 s-1). Our experimental and theoretical rate coefficients (which are in reasonable agreement at 298 K) show that the reaction of MECI with SO2 is significantly faster than MVK-oxide + SO2, demonstrating the substantial effect of resonance stabilization on Criegee intermediate reactivity.

Accurate fuel oxidation mechanisms can enable predictive capabilities that aid in advancing combustion technologies. High-level computational kinetics can yield reasonable rate coefficients with uncertainties, in some cases, below a factor of 2. Computed rate coefficients can be constrained further by optimizing against experimental data. Here, we explore the application of genetic algorithm (GA) optimization to constrain computed rate coefficients in complex fuel oxidation mechanisms in conjunction with temperature-dependent species mole fractions from jet-stirred reactor (JSR) measurements. Cyclohexane is a model candidate for understanding the reactivity of cyclic fuels. In this work, we optimize the rate coefficients of the most recent literature cyclohexane mechanism, which incorporates theoretically computed rate coefficients for the reaction networks stemming from the first and second O2 addition pathways, against the experimental results of two separate literature JSR studies. Optimization consistency is evaluated by carrying out three GA optimizations: fitting to the temperature-dependent species mole fractions in each JSR experiment separately and simultaneously fitting the species mole fractions in both experiments. Local sensitivity analyses are used to identify five influential low-temperature oxidation reactions for optimization. Although the three optimizations do not yield identical rate coefficients, the direction of change in all five rate coefficients is consistent among the three optimizations. Performance of the models from the three optimizations is assessed against literature ignition delay times with differences in the level of agreement observed among the different optimizations. Comparisons are made with our recent optimization work of a cyclopentane oxidation master-equation model against time-resolved species concentrations, and insights and improvements of the strategy for constraining rate coefficients using GA optimization are discussed.

Alkene ozonolysis generates short-lived Criegee intermediates that are a significant source of hydroxyl (OH) radicals. This study demonstrates that roaming of the separating OH radicals can yield alternate hydroxycarbonyl products, thereby reducing the OH yield. Specifically, hydroxybutanone has been detected as a stable product arising from roaming in the unimolecular decay of the methyl-ethyl-substituted Criegee intermediate (MECI) under thermal flow cell conditions. The dynamical features of this novel multistage dissociation plus a roaming unimolecular decay process have also been examined with ab initio kinetics calculations. Experimentally, hydroxybutanone isomers are distinguished from the isomeric MECI by their higher ionization threshold and distinctive photoionization spectra. Moreover, the exponential rise of the hydroxybutanone kinetic time profile matches that for the unimolecular decay of MECI. A weaker methyl vinyl ketone (MVK) photoionization signal is also attributed to OH roaming. Complementary multireference electronic structure calculations have been utilized to map the unimolecular decay pathways for MECI, starting with 1,4 H atom transfer from a methyl or methylene group to the terminal oxygen, followed by roaming of the separating OH and butanonyl radicals in the long-range region of the potential. Roaming via reorientation and the addition of OH to the vinyl group of butanonyl is shown to yield hydroxybutanone, and subsequent C-O elongation and H-transfer can lead to MVK. A comprehensive theoretical kinetic analysis has been conducted to evaluate rate constants and branching yields (ca. 10-11%) for thermal unimolecular decay of MECI to conventional and roaming products under laboratory and atmospheric conditions, consistent with the estimated experimental yield (ca. 7%).

The efficient utilization of lignin, the direct source of renewable aromatics, into value-added renewable chemicals is an important step towards sustainable biorefinery practices. Nevertheless, owing to the random heterogeneous structure and limited solubility, lignin utilization has been primarily limited to burning for energy. The catalytic depolymerization of lignin has been proposed and demonstrated as a viable route to sustainable biorefinery, however, low yields and poor selectivity of products, high char formation, and limited to no recycling of transition-metal-based catalyst involved in lignin depolymerization demands attention to enable practical-scale lignocellulosic biorefineries. In this study, we demonstrate the catalytic depolymerization of ionic liquid-based biorefinery poplar lignin into guaiacols over a reusable zirconium phosphate supported palladium catalyst. The essence of the study lies in the high conversion (>80 %), minimum char formation (7–16 %), high yields of guaiacols (up to 200 mg / g of lignin), and catalyst reusability. Both solid residue, liquid stream, and gaseous products were thoroughly characterized using ICP-OES, PXRD, CHN analysis, GC-MS, GPC, and 2D NMR to understand the hydrogenolysis pathway.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Using chemical kinetic modeling and statistical analysis, we investigate the possibility of correlating key chemical "markers"-typically small molecules-formed during very lean (φ ∼0.001) oxidation experiments with near-stoichiometric (φ ∼1) fuel ignition properties. One goal of this work is to evaluate the feasibility of designing a fuel-screening platform, based on small laboratory reactors that operate at low temperatures and use minimal fuel volume. Buras et al. [Combust. Flame 2020, 216, 472-484] have shown that convolutional neural net (CNN) fitting can be used to correlate first-stage ignition delay times (IDTs) with OH/HO2measurements during very lean oxidation in low-T flow reactors with better than factor-of-2 accuracy. In this work, we test the limits of applying this correlation-based approach to predict the low-temperature heat release (LTHR) and total IDT, including the sensitivity of total IDT to the equivalence ratio, φ. We demonstrate that first-stage IDT can be reliably correlated with very lean oxidation measurements using compressed sensing (CS), which is simpler to implement than CNN fitting. LTHR can also be predicted via CS analysis, although the correlation quality is somewhat lower than for first-stage IDT. In contrast, the accuracy of total IDT prediction at φ = 1 is significantly lower (within a factor of 4 or worse). These results can be rationalized by the fact that the first-stage IDT and LTHR are primarily determined by low-temperature chemistry, whereas total IDT depends on low-, intermediate-, and high-temperature chemistry. Oxidation reactions are most important at low temperatures, and therefore, measurements of universal molecular markers of oxidation do not capture the full chemical complexity required to accurately predict the total IDT even at a single equivalence ratio. As a result, we find that φ-sensitivity of ignition delay cannot be predicted at all using solely correlation with lean low-T chemical speciation measurements.

Time-resolved spectroscopies using high-energy photons in the vacuum ultraviolet (VUV) to the X-ray region of the electromagnetic spectrum, have proven to be powerful probes of chemical dynamics. These high-energy photons can access valence and core orbitals of molecules and materials, providing key information on molecular and electronic structure and their time evolution. This report details the development of table-top sources of extreme ultraviolet (XUV) and VUV pulses at Sandia National Laboratories for use in studies of gas phase chemical dynamics. Femtosecond duration XUV pulses are produced using laser-driven high harmonic generation and their detected range span ~40-140 eV photon energies. These pulses are used in conjunction with ultraviolet pulses in a pump-probe scheme to study excited state dynamics of gas phase molecules. VUV pulses at 7.75 eV are generated using a four-wave-mixing scheme driven by 800 nm and 266 nm pulses in an argon-filled hollow-core fiber.

We present a new experimental methodology for detailed experimental investigations of depolymerization reactions over solid catalysts. This project aims to address a critical need in fundamental research on chemical upcycling of polymers – the lack of rapid, sensitive, isomerselective probing techniques for the detection of reaction intermediates and products. Our method combines a heterogeneous catalysis reactor for the study of multiphase (gas/polymer melt/solid) systems, coupled to a vacuum UV photoionization time-of-flight mass spectrometer. This apparatus draws on our expertise in probing complex gas-phase chemistry and enables highthroughput, detailed chemical speciation measurements of the gas phase above the catalyst, providing valuable information on the heterogeneous catalytic reactions. Using this approach, we investigated the depolymerization of high-density polyethylene (HDPE) over Ir-doped zeolite catalysts. We showed that the product distribution was dominated by low-molecular weight alkenes with terminal C=C double bonds and revealed the presence of many methyl-substituted alkenes and alkanes, suggesting extensive methyl radical chemistry. In addition, we investigated the fundamental reactivity of model oligomer molecules n-butane and isobutane over ZSM-5 zeolites. We demonstrated the first direct detection of methyl radical intermediates, confirming the key role of methyl in zeolite-catalyzed activation of alkanes. Our results show the potential of this experimental method to achieve deep insight into the complex depolymerization reactions and pave the way for detailed mechanistic studies, leading to increased fundamental understanding of key processes in chemical upcycling of polymers.

We present a combined experimental and theoretical investigation of the autoignition chemistry of a prototypical cyclic hydrocarbon, cyclopentane. Experiments using a high-pressure photolysis reactor coupled to time-resolved synchrotron VUV photoionization mass spectrometry directly probe the short-lived radical intermediates and products in cyclopentane oxidation reactions. We detect key peroxy radical intermediates ROO and OOQOOH, as well as several hydroperoxides, formed by second O2 addition. Automated quantum chemical calculations map out the R + O2 + O2 reaction channels and demonstrate that the detected intermediates belong to the dominant radical chain-branching pathway: ROO (+ O2) → γ-QOOH + O2 → γ-OOQOOH → products. ROO, OOQOOH, and hydroperoxide products of second-O2 addition undergo extensive dissociative ionization, making their experimental assignment challenging. We use photoionization dynamics calculations to aid in their characterization and report the absolute photoionization spectra of isomerically pure ROO and γ-OOQOOH. A global statistical fit of the observed kinetics enables reliable quantification of the time-resolved concentrations of these elusive, yet critical species, paving the way for detailed comparisons with theoretical predictions from master-equation-based models.

Abstract not provided.

Oxiranes are a class of cyclic ethers formed in abundance during low-temperature combustion of hydrocarbons and biofuels, either via chain-propagating steps that occur from unimolecular decomposition of β-hydroperoxyalkyl radicals (β-̇QOOH) or from reactions of HOȮ with alkenes. Ethyloxirane is one of four alkyl-substituted cyclic ether isomers produced as an intermediate from n-butane oxidation. While rate coefficients for β-̇QOOH → ethyloxirane + ȮH are reported extensively, subsequent reaction mechanisms of the cyclic ether are not. As a result, chemical kinetics mechanisms commonly adopt simplified chemistry to describe ethyloxirane consumption by convoluting several elementary reactions into a single step, which may introduce mechanism truncation error—uncertainty derived from missing or incomplete chemistry. The present work provides fundamental insight on reaction mechanisms of ethyloxirane in support of ongoing efforts to minimize mechanism truncation error. Reaction mechanisms are inferred from the detection of products during chlorine atom-initiated oxidation experiments using multiplexed photoionization mass spectrometry conducted at 10 Torr and temperatures of 650 K and 800 K. To complement the experiments, calculations of stationary point energies were conducted using the ccCA-PS3 composite method on ̇R + O2 potential energy surfaces for the four ethyloxiranyl radical isomers, which produced barrier heights for 24 reaction pathways. In addition to products from ̇QOOH → cyclic ether + ȮH and ̇R + O2 → conjugate alkene + HOȮ, both of which were significant pathways and are prototypical to alkane oxidation, other species were identified from ring-opening of both ethyloxiranyl and ̇QOOH radicals. The latter occurs when the unpaired electron is localized on the ether group, causing the initial ̇QOOH structure to ring-open and form a resonance-stabilized ketohydroperoxide-type radical. The present work provides the first analysis of ethyloxirane oxidation chemistry, which reveals that consumption pathways are complex and may require an expansion of submechanisms to increase the fidelity of chemical kinetics mechanisms.

Abstract not provided.

Abstract not provided.

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 (HOOPO, 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, O2 concentration, and temperature on the competition among radical decomposition and second O2 addition pathways.