2,4,dimethyloxetane is an important cyclic ether intermediate that is produced from hydroperoxyalkyl (QOOH) radicals in the low-temperature combustion of n -pentane. However, the reaction mechanisms and rates of consumption pathways remain unclear. In the present work, the pressure- and temperature-dependent kinetics of seven cyclic ether peroxy radicals, which stem from 2,4,dimethyloxetane via H-abstraction and O 2 addition, were determined. The automated kinetic workflow code, KinBot, was used to model the complexity of the chemistry in a stereochemically resolved manner and solve the resulting master equations from 300–1000 K and from 0.01–100 atm. The main conclusions from the calculations include (i) diastereomeric cyclic ether peroxy radicals show significantly different reactivities, (ii) the stereochemistry of the peroxy radical determines which QOOH isomerization steps are possible, (iii) conventional QOOH decomposition pathways, such as cyclic ether formation and HO 2 elimination, compete with ring-opening reactions, which primarily produce OH radicals, the outcome of which is sensitive to stereochemistry. Ring-opening reactions lead to unique products, such as unsaturated, acyclic peroxy radicals, that form direct connections with species present in other chemical kinetics mechanisms through "cross-over" reactions that may complicate the interpretation of experimental results from combustion of n-pentane and, by extension, other alkanes. For example, one cross-over reaction involving 1-hydroperoxy-4-pentanone-2-yl produces 2-(hydroperoxymethyl)-3-butanone-1-yl, which is an iso-pentane-derived ketohydroperoxide (KHP). At atmospheric pressure, the rate of chemical reactions of all seven peroxy radicals compete with that of collisional stabilization, resulting in well-skipping reactions. However, at 100 atm, only one out of seven peroxy radicals undergoes significant well-skipping reactions. The rates produced from the master equation calculations provide the first foundation for the development of detailed sub-mechanisms for cyclic ether intermediates. In addition, analysis of the complex reaction mechanisms of 2,4-dimethyloxetane-derived peroxy radicals provides insights into the effects of stereoisomers on reaction pathways and product yields.
more »
« less
Pressure-dependent kinetics of peroxy radicals formed in isobutanol combustion
Bio-derived isobutanol has been approved as a gasoline additive in the US, but our understanding of its combustion chemistry still has significant uncertainties. Detailed quantum calculations could improve model accuracy leading to better estimation of isobutanol's combustion properties and its environmental impacts. This work examines 47 molecules and 38 reactions involved in the first oxygen addition to isobutanol's three alkyl radicals located α, β, and γ to the hydroxide. Quantum calculations are mostly done at CCSD(T)-F12/cc-pVTZ-F12//B3LYP/CBSB7, with 1-D hindered rotor corrections obtained at B3LYP/6-31G(d). The resulting potential energy surfaces are the most comprehensive isobutanol peroxy networks published to date. Canonical transition state theory and a 1-D microcanonical master equation are used to derive high-pressure-limit and pressure-dependent rate coefficients, respectively. At all conditions studied, the recombination of γ-isobutanol radical with O 2 forms HO 2 + isobutanal. The recombination of β-isobutanol radical with O 2 forms a stabilized hydroperoxy alkyl radical below 400 K, water + an alkoxy radical at higher temperatures, and HO 2 + an alkene above 1200 K. The recombination of β-isobutanol radical with O 2 results in a mixture of products between 700–1100 K, forming acetone + formaldehyde + OH at lower temperatures and forming HO 2 + alkenes at higher temperatures. The barrier heights, high-pressure-limit rates, and pressure-dependent kinetics generally agree with the results from previous quantum chemistry calculations. Six reaction rates in this work deviate by over three orders of magnitude from kinetics in detailed models of isobutanol combustion, suggesting the rates calculated here can help improve modeling of isobutanol combustion and its environmental fate.
more »
« less
- Award ID(s):
- 1709993
- NSF-PAR ID:
- 10296867
- Date Published:
- Journal Name:
- Physical Chemistry Chemical Physics
- Volume:
- 22
- Issue:
- 35
- ISSN:
- 1463-9076
- Page Range / eLocation ID:
- 19802 to 19815
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The reactions of thioformaldehyde (H 2 CS) with OH radicals and assisted by a single water molecule have been investigated using high level ab initio quantum chemistry calculations. The H 2 CS + ˙OH reaction can in principle proceed through: (1) abstraction, and (2) addition pathways. The barrier height for the addition reaction in the absence of a catalyst was found to be −0.8 kcal mol −1 , relative to the separated reactants, which has a ∼1.0 kcal mol −1 lower barrier than the abstraction channel. The H 2 CS + ˙OH reaction assisted by a single water molecule reduces the barrier heights significantly for both the addition and abstraction channels, to −5.5 and −6.7 kcal mol −1 respectively, compared to the un-catalyzed H 2 CS + ˙OH reaction. These values suggest that water lowers the barriers by ∼6.0 kcal mol −1 for both reaction paths. The rate constants for the H 2 CS⋯H 2 O + ˙OH and OH⋯H 2 O + H 2 CS bimolecular reaction channels were calculated using Canonical Variational Transition state theory (CVT) in conjunction with the Small Curvature Tunneling (SCT) method over the atmospherically relevant temperatures between 200 and 400 K. Rate constants for the H 2 CS + ˙OH reaction paths for comparison with the H 2 CS + ˙OH + H 2 O reaction in the same temperature range were also computed. The results suggest that the rate of the H 2 CS + ˙OH + H 2 O reaction is slower than that of the H 2 CS + ˙OH reaction by ∼1–4 orders of magnitude in the temperatures between 200 and 400 K. For example, at 300 K, the rates of the H 2 CS + ˙OH + H 2 O and H 2 CS + ˙OH reactions were found to be 2.2 × 10 −8 s −1 and 6.4 × 10 −6 s −1 , respectively, calculated using [OH] = 1.0 × 10 6 molecules cm −3 , and [H 2 O] = 8.2 × 10 17 molecules cm −3 (300 K, RH 100%) atmospheric conditions. Electronic structure calculations on the H 2 C(OH)S˙ product in the presence of 3 O 2 were also performed. The results show that H 2 CS is removed from the atmosphere primarily by reacting with ˙OH and O 2 to form thioformic acid, HO 2 , formaldehyde, and SO 2 as the main end products.more » « less
-
For more than half a century, pericyclic reactions have played an important role in advancing our fundamental understanding of cycloadditions, sigmatropic shifts, group transfer reactions, and electrocyclization reactions. However, the fundamental mechanisms of photochemically activated cheletropic reactions have remained contentious. Here we report on the simplest cheletropic reaction: the [2+1] addition of ground state 18 O-carbon monoxide (C 18 O, X 1 Σ + ) to D2-acetylene (C 2 D 2 ) photochemically excited to the first excited triplet (T1), second excited triplet (T2), and first excited singlet state (S1) at 5 K, leading to the formation of D2- 18 O-cyclopropenone (c-C 3 D 2 18 O). Supported by quantum-chemical calculations, our investigation provides persuasive testimony on stepwise cheletropic reaction pathways to cyclopropenone via excited state dynamics involving the T2 (non-adiabatic) and S1 state (adiabatic) of acetylene at 5 K, while the T1 state energetically favors an intermediate structure that directly dissociates after relaxing to the ground state. The agreement between experiments in low temperature ices and the excited state calculations signifies how photolysis experiments coupled with theoretical calculations can untangle polyatomic reactions with relevance to fundamental physical organic chemistry at the molecular level, thus affording a versatile strategy to unravel exotic non-equilibrium chemistries in cyclic, aromatic organics. Distinct from traditional radical–radical pathways leading to organic molecules on ice-coated interstellar nanoparticles (interstellar grains) in cold molecular clouds and star-forming regions, the photolytic formation of cyclopropenone as presented changes the perception of how we explain the formation of complex organics in the interstellar medium eventually leading to the molecular precursors of biorelevant molecules.more » « less
-
more » « less
Abstract The development and refinement of NOx chemical kinetic mechanisms have been instrumental in understanding and reducing NOx formation. However, relatively little work has been performed with NOx species as the oxidizer, and such experiments can provide unique insights into NOx kinetics. Furthermore, speciation data can often provide useful information that complements global measurements such as ignition delay times in facilitating mechanism refinement. To provide such speciation data in the H2‐NO2system, H2O measurements were performed using a fixed‐wavelength, direct absorption laser diagnostic near 1.39 µm behind reflected shock waves in fuel‐lean, near‐stoichiometric, and fuel‐rich mixtures of H2and NO2highly diluted in argon. Experiments were performed between 917 and 1782 K near atmospheric pressure. The H2O profiles obtained herein are markedly different from those using O2as the oxidizer obtained in a previous study. The GRI 3.0 mechanism was found to greatly underestimate the H2O formation, whereas two modern mechanisms were found to predict the H2O formation quite accurately except at colder temperatures for fuel‐rich conditions. Explanations for the differences between these mechanisms are given and discussed, with the conclusion that older mechanisms such as GRI 3.0 should not be used to model hydrocarbon/NOx combustion chemistry as they are lacking several key reactions and species, namely NO3and HONO. The discrepancy between models and data at lower temperatures could not be reconciled even when modifying two of the most‐sensitive reaction rates. To the best of the authors’ knowledge, this study presents the first shock‐tube speciation study in the H2‐NO2system.
-
more » « less
Abstract The mineral apatite, Ca10(PO4)6(F,OH,Cl)2, incorporates sulfur (S) during crystallization from S-bearing hydrothermal fluids and silicate melts. Our previous studies of natural and experimental apatite demonstrate that the oxidation state of S in apatite varies systematically as a function of oxygen fugacity (fO2). The S oxidation states –1 and –2 were quantitatively identified in apatite crystallized from reduced, S-bearing hydrothermal fluids and silicate melts by using sulfur K-edge X-ray absorption near-edge structure spectroscopy (S-XANES) where S 6+/ΣS in apatite increases from ~0 at FMQ-1 to ~1 at FMQ+2, where FMQ refers to the fayalite-magnetite-quartz fO2 buffer. In this study, we employ quantum-mechanical calculations to investigate the atomistic structure and energetics of S(-I) and S(-II) incorporated into apatite and elucidate incorporation mechanisms.
One S(-I) species (disulfide, S22−) and two S(-II) species (bisulfide, HS−, and sulfide, S2−) are investigated as possible forms of reduced S species in apatite. In configuration models for the simulation, these reduced S species are positioned along the c-axis channel, originally occupied by the column anions F, Cl, and OH in the end-member apatites. In the lowest-energy configurations of S-incorporated apatite, disulfide prefers to be positioned halfway between the mirror planes at z = 1/4 and 3/4. In contrast, the energy-optimized bisulfide is located slightly away from the mirror planes by ~0.04 fractional units in the c direction. The energetic stability of these reduced S species as a function of position along the c-axis can be explained by the geometric and electrostatic constraints of the Ca and O planes that constitute the c-axis channel.
The thermodynamics of incorporation of disulfide and bisulfide into apatite is evaluated by using solid-state reaction equations where the apatite host and a solid S-bearing source phase (pyrite and Na2S2(s) for disulfide; troilite and Na2S(s) for sulfide) are the reactants, and the S-incorporated apatite and an anion sink phase are the products. The Gibbs free energy (ΔG) is lower for incorporation with Na-bearing phases than with Fe-bearing phases, which is attributed to the higher energetic stability of the iron sulfide minerals as a source phase for S than the sodium sulfide phases. The thermodynamics of incorporation of reduced S is also evaluated by using reaction equations involving dissolved disulfide and sulfide species [HnS(aq)(2−n) and HnS(aq)(2−n); n = 0, 1, and 2] as a source phase. The ΔG of S-incorporation increases for fluorapatite and chlorapatite, and decreases for hydroxylapatite, as these species are protonated (i.e., as n changes from 0 to 2). These thermodynamic results demonstrate that the presence of reduced S in apatite is primarily controlled by the chemistry of magmatic and hydrothermal systems where apatite forms (e.g., an abundance of Fe; solution pH). Ultimately, our methodology developed for evaluating the thermodynamics of S incorporation in apatite as a function of temperature, pH, and composition is highly applicable to predicting the trace and volatile element incorporation in minerals in a variety of geological systems. In addition to solid-solid and solid-liquid equilibria treated here at different temperatures and pH, the methodology can be easily extended to different pressure conditions by just performing the quantum-mechanical calculations at elevated pressures.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1039/d0cp02872j
