Kinetics and Thermochemistry of 2,5-Dimethyltetrahydrofuran and Related Oxolanes: Next Next-Generation Biofuels

Time-resolved quantification of key species and mechanistic insights in low-temperature tetrahydrofuran oxidation

We investigate the kinetics and report the time-resolved concentrations of key chemical species in the oxidation of tetrahydrofuran (THF) at 7500 torr and 450-675 K. Experiments are carried out using high-pressure multiplexed photoionization mass spectrometry (MPIMS) combined with tunable vacuum ultraviolet radiation from the Berkely Lab Advanced Light Source. Intermediates and products are quantified using reference photoionization (PI) cross sections, when available, and constrained by a global carbon balance tracking approach at all experimental temperatures simultaneously for the species without reference cross sections. From carbon balancing, we determine time-resolved concentrations for the ROO˙ and ˙OOQOOH radical intermediates, butanedial, and the combined concentration of ketohydroperoxide (KHP) and unsaturated hydroperoxide (UHP) products stemming from the ˙QOOH + O2 reaction. Furthermore, we quantify a product that we tentatively assign as fumaraldehyde, which arises from UHP decomposition via H2O or ˙OH + H loss. The experimentally derived species concentrations are compared with model predictions using the most recent literature THF oxidation mechanism of Fenard et al., (Combust. Flame, 2018, 191, 252-269). Our results indicate that the literature mechanism significantly overestimates THF consumption and the UHP + KHP concentration at our conditions. The model predictions are sensitive to the rate coefficient for the ROO˙ isomerization to ˙QOOH, which is the gateway for radical chain propagating and branching pathways. Comparisons with our recent results for cyclopentane (Demireva et al., Combust. Flame, 2023, 257, 112506) provide insights into the effect of the ether group on reactivity and highlight the need to determine accurate rate coefficients of ROO˙ isomerization and subsequent reactions.

Unimolecular Reactions of 2,4-Dimethyloxetanyl Radicals

Alkyl-substituted oxetanes are cyclic ethers formed via unimolecular reactions of QOOH radicals produced via a six-membered transition state in the preceding isomerization step of organic peroxy radicals, ROO. Owing to radical isomer-specific formation pathways, cyclic ethers are unambiguous proxies for inferring QOOH reaction rates. Therefore, accounting for subsequent oxidation of cyclic ethers is important in order to accurately determine rates for QOOH → products. Cyclic ethers can react via unimolecular reaction (ring-opening) or via bimolecular reaction with O₂ to form cyclic ether-peroxy adducts. The computations herein provide reaction mechanisms and theoretical rate coefficients for the former type in order to determine competing pathways for the cyclic ether radicals. Rate coefficients of unimolecular reactions of 2,4-dimethyloxetanyl radicals were computed using master equation modeling from 0.01 to 100 atm and from 300 to 1000 K. Coupled-cluster methods were utilized for stationary-point energy calculations, and uncertainties in the computed rate coefficients were accounted for using variation in barrier heights and in well depths. The potential energy surfaces reveal accessible channels to several species via crossover reactions, such as 2-methyltetrahydrofuran-5-yl and pentanonyl isomers. For the range of temperature over which 2,4-dimethyloxetane forms during n-pentane oxidation, the following are the major channels: 2,4-dimethyloxetan-1-yl → acetaldehyde + allyl, 2,4-dimethyloxetan-2-yl → propene + acetyl, and 2,4-dimethyloxetan-3-yl → 3-butenal + methyl, or, 1-penten-3-yl-4-ol. Well-skipping reactions were significant in a number of channels and also exhibited a markedly different pressure dependence. The calculations show that rate coefficients for ring-opening are approximately an order of magnitude lower for the tertiary 2,4-dimethyloxetanyl radicals than for the primary and secondary 2,4-dimethyloxetanyl radicals. Unlike for reactions of the corresponding ROO radicals, however, unimolecular rate coefficients are independent of the stereochemistry. Moreover, rate coefficients of cyclic ether radical ring-opening are of the same order of magnitude as O₂ addition, underscoring the point that a competing network of reactions is necessary to include for accurate chemical kinetics modeling of species profiles for cyclic ethers.

Correcting the Experimental Enthalpies of Formation of Some Members of the Biologically Significant Sulfenic Acids Family

Sulfenic acids are important intermediates in the oxidation of cysteine thiol groups in proteins by reactive oxygen species. The mechanism is influenced heavily by the presence of polar groups, other thiol groups, and solvent, all of which determines the need to compute precisely the energies involved in the process. Surprisingly, very scarce experimental information exists about a very basic property of sulfenic acids, the enthalpies of formation. In this Article, we use high level quantum chemical methods to derive the enthalpy of formation at 298.15 K of methane-, ethene-, ethyne-, and benzenesulfenic acids, the only ones for which some experimental information exists. The methods employed were tested against well-known experimental data of related species and extensive CCSD(T) calculations. Our best results consistently point out to a much lower enthalpy of formation of methanesulfenic acid, CH₃SOH (Δ_fH⁰(298.15K) = -35.1 ± 0.4 kcal mol^-1), than the one reported in the NIST thermochemical data tables. The enthalpies of formation derived for ethynesulfenic acid, HC≡CSOH, +32.9 ± 1.0 kcal/mol, and benzenesulfenic acid, C₆H₅SOH, -2.6 ± 0.6 kcal mol^-1, also differ markedly from the experimental values, while the enthalpy of formation of ethenesulfenic acid CH₂CHSOH, not available experimentally, was calculated as -11.2 ± 0.7 kcal mol^-1.

Thermal Decomposition of 2-Methyltetrahydrofuran behind Reflected Shock Waves over the Temperature Range of 1179-1361 K

The thermal unimolecular decomposition of 2-methyltetrahydrofuran (2-MTHF) was studied behind reflected shock waves in a single-pulse shock tube over the temperature range of 1179-1361 K and pressure range of 9-17 atm. Methane, ethylene, ethane, 1,3-butadiene, propylene, acetaldehyde, and acetylene were identified as products in the decomposition of 2-MTHF. A reaction scheme was proposed to explain the mechanism for the observed products. The experimentally determined rate coefficients were best fit to an Arrhenius expression for the overall decomposition and is represented as k_total^exp(1179-1361 K) = (3.23 ± 0.59) × 10¹¹ s^-1 exp(-51.3 ± 1.4 kcal mol^-1/RT). Quantum chemistry methods were used to calculate the energetics and kinetics of various possible unimolecular dissociation pathways involved in the thermal decomposition of 2-MTHF. The initial decomposition of 2-MTHF occurs predominantly via ring-methyl (C-CH₃) single bond fission, leading to the formation of tetrahydrofuran (C₄H₇O) radical, and methyl radical was found to be the major reaction compared to all the possible initial bond fission, ring opening, and molecular elimination channels. The temperature-dependent rate coefficients for the unimolecular dissociation of 2-MTHF were calculated using the RRKM (Rice-Ramsperger-Kassel-Marcus) theory in combination with the CCSD(T)/cc-pVTZ//B3LYP/cc-pVTZ level of electronic structure calculations over the temperature range of 800-1500 K. The computed high-pressure limiting rate coefficients for the initial decomposition of 2-MTHF through C-CH₃ single bond fission channel were found to be ∼2 times higher in the temperatures between 800 and 900 K, and above this temperature, they agree well with the values reported in the literature.

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.

Reactions of Three Lactones with Cl, OD, and O3: Atmospheric Impact and Trends in Furan Reactivity

Lactones, cyclic esters of hydroxycarboxylic acids, are interesting biofuel candidates as they can be made from cellulosic biomass and have favorable physical and chemical properties for distribution and use. The reactions of γ-valerolactone (GVL), γ-crotonolactone (2(5H)-F), and α-methyl-γ-crotonolactone (3M-2(5H)-F) with Cl, OD, and O₃ were investigated in a static chamber at 700 Torr and 298 ± 2 K. The relative rate method was used to determine k_GVL+Cl = (4.56 ± 0.51) × 10^-11, k_GVL+OD = (2.94 ± 0.41) × 10^-11, k_2(5H)-F+Cl = (2.94 ± 0.41) × 10^-11, k_2(5H)-F+OD = (4.06 ± 0.073) × 10^-12, k_{3M-2(5H)-F+Cl} = (16.1 ± 1.8) × 10^-11, and k_{3M-2(5H)-F+OD} = (12.6 ± 0.52) × 10^-12, all rate coefficients in units of cm³ molecule^-1 s^-1. An absolute rate method was used to determine k_2(5H)-F+O3 = (6.73 ± 0.18) × 10^-20 and k_{3M-2(5H)-F+O3} = (5.42 ± 1.23) × 10^-19 in units of cm³ molecule^-1 s^-1. Products were identified for reactions of the lactones with Cl. In the presence of O₂ the products are formic acid (HCOOH), formyl chloride (CHClO), and phosgene (CCl₂O), and also maleic anhydride (C₂H₂(CO)₂O) for 2(5H)-F. In addition both reactions produced a number of unidentified products that likely belong to molecules with the ring-structure intact. A review of literature data for reactions of other furans show that the reactivity of the lactones are generally lower compared to that of corresponding compounds without the carbonyl group.

Atmospheric Chemistry of Tetrahydrofuran, 2-Methyltetrahydrofuran, and 2,5-Dimethyltetrahydrofuran: Kinetics of Reactions with Chlorine Atoms, OD Radicals, and Ozone

FTIR smog chamber techniques were used to study the kinetics of the gas-phase reactions of Cl atoms, OD radicals, and O3 with the five-membered ring-structured compounds tetrahydrofuran (C4H8O, THF), 2-methyltetrahydrofuran (CH3C4H7O, 2-MTHF), 2,5-dimethyltetrahydrofuran ((CH3)2C4H5O, 2,5-DMTHF), and furan (C4H4O). The rate coefficients determined using relative rate methods were kTHF+Cl = (1.96 ± 0.24) × 10(-10), kTHF+OD = (1.81 ± 0.27) × 10(-11), kTHF+O3 = (6.41 ± 2.90) × 10(-21), k2-MTHF+Cl = (2.65 ± 0.43) × 10(-10), k2-MTHF+OD = (2.41 ± 0.51) × 10(-11), k2-MTHF+O3 = (1.87 ± 0.82) × 10(-20), k2,5-DMTHF+OD = (4.56 ± 0.68) × 10(-11), k2,5-DMTHF+Cl = (2.84 ± 0.34) × 10(-10), k2,5-DMTHF+O3 = (4.58 ± 2.18), kfuran+Cl = (2.39 ± 0.27) × 10(-10), and kfuran+O3 = (2.60 ± 0.31) × 10(-18) molecules cm(-3) s(-1). Rate coefficients of the reactions with ozone were also determined using the absolute rate method under pseudo-first-order conditions. OD radicals, in place of OH radicals, were produced from CD3ONO to avoid spectral overlap of isopropyl and methyl nitrite with the reactants. The kinetics of OD radical reactions are expected to resemble the kinetics of OH radical reactions, and the rate coefficients of the reactions with OD radicals were used to calculate the atmospheric lifetimes with respect to reactions with OH radicals. The lifetimes of THF, 2-MTHF, and 2,5-DMTHF are approximately 15, 12, and 6 h, respectively.

Toward the Development of a Fundamentally Based Chemical Model for Cyclopentanone: High-Pressure-Limit Rate Constants for H Atom Abstraction and Fuel Radical Decomposition

Theoretical aspects of the development of a chemical kinetic model for the pyrolysis and combustion of a cyclic ketone, cyclopentanone, are considered. Calculated thermodynamic and kinetic data are presented for the first time for the principal species including 2- and 3-oxo-cyclopentyl radicals, which are in reasonable agreement with the literature. These radicals can be formed via H atom abstraction reactions by Ḣ and Ö atoms and ȮH, HȮ2, and ĊH3 radicals, the rate constants of which have been calculated. Abstraction from the β-hydrogen atom is the dominant process when ȮH is involved, but the reverse holds true for HȮ2 radicals. The subsequent β-scission of the radicals formed is also determined, and it is shown that recent tunable VUV photoionization mass spectrometry experiments can be interpreted in this light. The bulk of the calculations used the composite model chemistry G4, which was benchmarked in the simplest case with a coupled cluster treatment, CCSD(T), in the complete basis set limit.

A mechanistic and kinetic study on the decomposition of morpholine

The combustion chemistry of morpholine (C(4)H(8)ONH) has been experimentally investigated recently as a representative model compound for O- and N-containing structural entities in biomass. Detailed profiles of species indicate the self-breakdown reactions prevailing over oxidative decomposition reactions. In this study, we derive thermodynamic and kinetic properties pertinent to all plausible reactions involved in the self-decomposition of morpholine and its derived morphyl radicals as a crucial task in the development of comprehensive combustion mechanism. Potential energy surfaces have been mapped out for the decomposition of morpholine and the three morphyl radicals. RRKM-based calculations predict the self-decomposition of morpholine to be dominated by 1,3-intramolecular hydrogen shift into the NH group at all temperatures and pressures. Self-decomposition of morpholine is shown to provide pathways for the formation of the experimentally detected products such as ethenol and ethenamine. Energetic requirements of all self-decomposition of morphyl radicals are predicted to be of modest values (i.e., 20-40 kcal/mol) which in turn support the occurrence of breaking-down reactions into two-heavy-atom species and the generation of doubly unsaturated four-heavy-atom segments. Calculated thermochemical parameters (in terms of standard enthalpies of formation, standard entropies, and heat capacities) and kinetic parameters (in terms of reaction rate constants at a high pressure limit) should be instrumental in building a robust kinetic model for the oxidation of morpholine.