Search NASA⌕ Search

SEARCH · Search NASA

Results for “Master equation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Comprehensive Kinetics on the C 7 H 7 Potential Energy Surface under Combustion Conditions

The automated kinetics workflow code, KinBot, was used to explore and characterize the regions of the C 7 H 7 potential energy surface that are relevant to combustion environments and especially soot inception. We first explored the lowest-energy region, which includes the benzyl, fulvenallene + H, and cyclopentadienyl + acetylene entry points. We then expanded the model to include two higher-energy entry points, vinylpropargyl + acetylene and vinylacetylene + propargyl. The automated search was able to uncover the pathways from the literature. In addition, three important new routes were discovered: a lower-energy pathway connecting benzyl with vinylcyclopentadienyl, a decomposition mechanism from benzyl that results in side-chain hydrogen atom loss to produce fulvenallene + H, and shorter and lower energy routes to the dimethylene-cyclopentenyl intermediates. We systematically reduced the extended model to a chemically relevant domain composed of 63 wells, 10 bimolecular products, 87 barriers, and 1 barrierless channel and constructed a master equation using the CCSD(T)-F12a/cc-pVTZ//ωB97X-D/6-311++G(d,p) level of theory to provide rate coefficients for chemical modeling. Our calculated rate coefficients show excellent agreement with measured ones. We also simulated concentration profiles and calculated branching fractions from the important entry points to provide an interpretation of this important chemical landscape.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Radical–Radical Reactions in Molecular Weight Growth: The Phenyl + Propargyl Reaction

The mechanism for hydrocarbon ring growth in sooting environments is still the subject of considerable debate. The reaction of phenyl radical (C 6 H 5 ) with propargyl radical (H 2 CCCH) provides an important prototype for radical–radical ring-growth pathways. We studied this reaction experimentally over the temperature range of 300–1000 K and pressure range of 4–10 Torr using time-resolved multiplexed photoionization mass spectrometry. We detect both the C 9 H 8 and C 9 H 7 + H product channels and report experimental isomer-resolved product branching fractions for the C 9 H 8 product. We compare these experiments to theoretical kinetics predictions from a recently published study augmented by new calculations. Here, these ab initio transition state theory-based master equation calculations employ high-quality potential energy surfaces, conventional transition state theory for the tight transition states, and direct CASPT2-based variable reaction coordinate transition state theory (VRC-TST) for the barrierless channels. At 300 K only the direct adducts from radical–radical addition are observed, with good agreement between experimental and theoretical branching fractions, supporting the VRC-TST calculations of the barrierless entrance channel. As the temperature is increased to 1000 K we observe two additional isomers, including indene, a two-ring polycyclic aromatic hydrocarbon, and a small amount of bimolecular products C 9 H 7 + H. Our calculated branching fractions for the phenyl + propargyl reaction predict significantly less indene than observed experimentally. We present further calculations and experimental evidence that the most likely cause of this discrepancy is the contribution of H atom reactions, both H + indenyl (C 9 H 7 ) recombination to indene and H-assisted isomerization that converts less stable C 9 H 8 isomers into indene. Especially at low pressures typical of laboratory investigations, H-atom-assisted isomerization needs to be considered. Regardless, the experimental observation of indene demonstrates that the title reaction leads, either directly or indirectly, to the formation of the second ring in polycyclic aromatic hydrocarbons.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

OH Roaming during the Ozonolysis of α-Pinene: A New Route to Highly Oxygenated Molecules?

The formation of low-volatility organic compounds in the ozonolysis of α-pinene, the dominant atmospheric monoterpene, provides an important route to aerosol formation. Here, in this work, we consider a previously unexplored set of pathways for the formation of highly oxygenated molecules in α-pinene ozonolysis. Pioneering, direct experimental observations of Lester and co-workers have demonstrated a significant production of hydroxycarbonyl products in the dissociation of Criegee intermediates. Theoretical analyses indicate that this production arises from OH roaming-induced pathways during the OO fission of the vinylhydroperoxides (VHPs), which in turn come from internal H transfers in the Criegee intermediates. Ab initio kinetics computations are used here to explore the OH roaming-induced channels that arise from the ozonolysis of α-pinene. For computational reasons, the calculations consider a surrogate for α-pinene, where two spectator methyl groups are replaced with H atoms. Multireference electronic structure calculations are used to illustrate a variety of energetically accessible OH roaming pathways for the four VHPs arising from the ozonolysis of this α-pinene surrogate. Ab initio transition-state theory-based master equation calculations indicate that for the dissociation of stabilized VHPs, these OH roaming pathways are kinetically significant with a branching that generally increases from ~20% at room temperature up to ~70% at lower temperatures representative of the troposphere. For one of the VHPs, this branching already exceeds 60% at room temperature. For the overall ozonolysis process, these branching ratios would be greatly reduced by a limited branching to the stabilized VHP, although there would also be some modest roaming fraction for the nonthermal VHP dissociation process. The strong exothermicities of the roaming-induced isomerizations/additions and abstractions suggest new routes to fission of the cyclobutane rings. Such ring fissions would facilitate further autoxidation reactions, thereby providing a new route for producing highly oxygenated nonvolatile precursors to aerosol formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Theoretical Kinetics Predictions for Reactions on the NH 2 O Potential Energy Surface

Recent modeling studies of ammonia oxidation, which are motivated by the prospective role of ammonia as a zero-carbon fuel, have indicated significant discrepancies among the existing literature mechanisms. In this study, high-level theoretical kinetics predictions have been obtained for reactions on the NH 2 O potential energy surface, including the NH 2 + O, HNO + H, and NH + OH reactions. These reactions have previously been highlighted as important reactions in NH 3 oxidation with high sensitivity and high uncertainty. The potential energy surface is explored with coupled cluster calculations, including large basis sets and high-level corrections to yield high-accuracy (~0.2 kcal/mol 2σ uncertainty) estimates of the stationary point energies. Variational transition state theory is used to predict the microcanonical rate constants, which are then incorporated in master equation treatments of the temperature- and pressure-dependent kinetics. For radical–radical channels, the microcanonical rates are obtained from variable reaction coordinate transition state theory implementing directly evaluated multireference electronic energies. The analysis yields predictions for the total rate constants as well as the branching ratios. We find that the NO + H 2 channel contributes 10% of the total NH 2 + O flux at combustion temperatures, although this channel is not included in modern NH 3 oxidation mechanisms. Modeling is used to illustrate the ramifications of these rate predictions on the kinetics of NH 3 oxidation and NO x formation. Finally, the present results for NH 2 + O are important for predicting the chain branching and formation of NO in the oxidation of NH 3 and thermal DeNO x .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Infrared Fingerprint and Unimolecular Decay Dynamics of the Hydroperoxyalkyl Intermediate (•QOOH) in Cyclopentane Oxidation

A transient carbon-centered hydroperoxyalkyl intermediate (•QOOH) in the oxidation of cyclopentane is identified by IR action spectroscopy with time-resolved unimolecular decay to hydroxyl (OH) radical products that are detected by UV laser-induced fluorescence. Two nearly degenerate •QOOH isomers, β- and γ-QOOH, are generated by H atom abstraction of the cyclopentyl hydroperoxide precursor. Fundamental and first overtone OH stretch transitions and combination bands of •QOOH are observed and compared with anharmonic frequencies computed by second-order vibrational perturbation theory. An OH stretch transition is also observed for a conformer arising from torsion about a low-energy CCOO barrier. Definitive identification of the β-QOOH isomer relies on its significantly lower transition state (TS) barrier to OH products, which results in rapid unimolecular decay and near unity branching to OH products. A benchmarking approach is utilized to compute high-accuracy stationary point energies, most importantly TS barriers, for cyclopentane oxidation (C 5 H 9 O 2 ), building on higher level reference calculations for ethane oxidation (C 2 H 5 O 2 ). The experimental OH product appearance rates are compared with computed statistical microcanonical rates using RRKM theory, including heavy-atom tunneling, thereby validating the computed TS barrier. The results are extended to thermal unimolecular decay rate constants at temperatures and pressures relevant to cyclopentane combustion via master-equation modeling. Furthermore, the various torsional and ring puckering states of the wells and transition states are explicitly considered in these calculations.

Chemical calculations↗
DOE OSTI2026

Rate constants and third-body collision efficiencies for recombination of Na with OH and O2: Implications for flame inhibition by alkali salts

Flame inhibition by alkali metals has implications for solid fuel combustion, fire safety, and a number of industrial processes. While the mechanism of inhibition is fairly well understood, details related to thermodynamic properties and rate constants are still in question. In the present work, the recombination of Na with OH (R5) and O2 (R7), respectively, was characterized theoretically, and the implications for modeling laminar premixed hydrogen flames doped with sodium species were examined. Third-body collision efficiencies and low-pressure-limit rate constants were obtained using newly fitted ab initio-based potential energy surfaces, classical trajectories, and one-dimensional master equation calculations. The results are consistent with available experimental results and aid in the present modeling study by providing rate information for bath gases and conditions that remain unexplored experimentally. Chemical kinetic modeling of relative Na and absolute H and OH profiles in H2-fueled laminar, premixed flames doped with a sodium salt shows that the most important radical removal cycle is the sequence Na + OH (+M) → NaOH (+M) (R5), NaOH + H → Na + H2O (R12), even under oxidizing conditions. A secondary cycle, Na + O2 (+M) ⇄ NaO2 (+M) (R7), NaO2 + OH → NaOH + O2 (R14), is less important due to the low thermal stability of NaO2. In most flames, reaction R7 is partially equilibrated, and reaction R14 becomes rate-limiting for the second cycle. The flame analysis supports a lower value of k14 than indicated by recent work on KO2 + OH, but more work is required to confirm this.

Jasper, Ahren W.↗
DOE OSTI2023

On the Mechanism of Polaritonic Rate Suppression from Quantum Transition Paths

Polariton chemistry holds promise for facilitating mode-selective chemical reactions, but the underlying mechanism behind the rate modifications observed under strong vibrational coupling is not well understood. Using the recently developed quantum transition path theory, we have uncovered a mechanism of resonant suppression of a thermal reaction rate in a simple model polaritonic system consisting of a reactive mode in a bath confined to a lossless microcavity with a single photon mode. Here we observed the formation of a polariton during rate-limiting transitions on reactive pathways and identified the concomitant rate suppression as being due to hybridization between the reactive mode and the cavity mode, which inhibits bath-mediated tunneling. The transition probabilities that define the quantum master equation can be directly translated into a visualization of the corresponding polariton energy landscape. This landscape exhibits a double funnel structure with a large barrier between the initial and final states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Kinetics and Mechanism of the Singlet Oxygen Atom Reaction with Dimethyl Ether

Here, we combine in situ laser spectroscopy, quantum chemistry, and kinetic calculations to study the reaction of a singlet oxygen atom with dimethyl ether. Infrared laser absorption spectroscopy and Faraday rotation spectroscopy are used for the detection and quantification of the reaction products OH, H 2 O, HO 2 , and CH 2 O on submillisecond time scales. Fitting temporal profiles of products with simulations using an in-house reaction mechanism allows product branching to be quantified at 30, 60, and 150 Torr. The experimentally determined product branching agrees well with master equation calculations based on electronic structure data and transition state theory. The calculations demonstrate that the dimethyl peroxide (CH 3 OOCH 3 ) generated via O-insertion into the C–O bond undergoes subsequent dissociation to CH 3 O + CH 3 O through energetically favored reactions without an intrinsic barrier. This O-insertion mechanism can be important for understanding the fate of biofuels leaking into the atmosphere and for plasma-based biofuel processing technologies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Interplay between Ultrafast Electronic and Librational Dynamics in Liquid Nitrobenzene Probed with Two-Color Four-Wave Mixing

Here, we present an experimental and theoretical study of the interplay between ultrafast electron dynamics and librational dynamics in liquid nitrobenzene. A femtosecond ultraviolet pulse and two femtosecond near-infrared pulses interact with nitrobenzene molecules, generating a four-wave mixing nonlinear signal measured in the Optical Kerr Effect geometry. The signal is measured to be nonzero only at negative time delays, corresponding to the near-infrared pulses arriving before the ultraviolet pulse. We perform time-dependent Quantum Master Equation calculations with classical libration to simulate the experiment. The simulations support the conclusion that the near-infrared pulses launch librational motion while creating electronic coherences resulting in a libration-modulated electronic nonlinear response. The analysis of the phase-matched four-wave mixing signals suggests a nonparametric process leaving the molecules in an excited electronic state, providing new insight into ultrafast nonlinear optical interactions in liquids and advancing toward probing ultrafast electronic coherences in complex molecular liquids.

Shivaram, Niranjan [Lawrence Berkeley National Lab↗
DOE OSTI2021

Origin of Chirality Induced Spin Selectivity in Photoinduced Electron Transfer

Here we propose a mechanism by which spin-polarization can be generated dynamically in chiral molecular systems undergoing photoinduced electron transfer. The proposed mechanism explains how spin-polarization emerges in systems where charge transport is dominated by incoherent hopping, mediated by spin-orbit and electronic exchange couplings through an intermediate charge transfer state. We derive a simple expression for the spin-polarization that predicts a nonmonotonic temperature dependence, consistent with recent experiments, and a maximum spin-polarization that is independent of the magnitude of the spin-orbit coupling. We validate this theory using approximate quantum master equations and the numerically exact hierarchical equations of motion. The proposed mechanism of chirality induced spin selectivity should apply to many chiral systems, and the ideas presented here have implications for the study of spin transport at temperatures relevant to biology and provide simple principles for the molecular control of spins in fluctuating environments.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Predictive Chemical Kinetic Modeling: Where We Succeed, Where We Struggle, and What Comes Next

Chemical kinetic modeling plays a foundational role in fields ranging from energy to environmental science, pharmaceuticals, and advanced materials. The past two decades have seen remarkable progress, particularly in modeling gas-phase reactions for thermochemical processes, leading to impactful industrial applications such as steam cracking and air quality management. However, new challenges are emerging. The successful development of systematic methodologies for the description of gas-phase kinetics opens the possibility to apply the same approach to the study of more challenging systems. Here, we review recent advances, including ab initio transition state theory-based master equation estimation of elementary rates, automated mechanism generation, machine-learning-assisted kinetics, and uncertainty quantification, and discuss the advances needed to apply the same methodological approach in areas such as heterogeneous catalysis, electrochemistry, liquid-phase and solid-state reactivity, and multiscale model integration. We advocate for the development of targeted tools, especially methods that go beyond empirical tuning toward first-principles-based predictions. We highlight the need for accessible software and AIaugmented workflows to democratize modeling for industry and academia alike. In this perspective, we call attention to not only what has worked but also what remains unsolved, advocating to avoid overemphasizing successes in scientific works at the expense of realism. The next decade should focus on predictive capability, physical accuracy, and community infrastructure (e.g., databases and services) to enable innovation across diverse fields. We argue that kinetic modeling, properly equipped, can accelerate discovery far beyond its traditional domains.

ab initio calculations↗
DOE OSTI2021

Substitution Reactions in the Pyrolysis of Acetone Revealed through a Modeling, Experiment, Theory Paradigm

The development of high-fidelity mechanisms for chemically reactive systems is a challenging process that requires the compilation of rate descriptions for a large and somewhat ill-defined set of reactions. The present unified combination of modeling, experiment, and theory provides a paradigm for improving such mechanism development efforts. Here we combine broadband rotational spectroscopy with detailed chemical modeling based on rate constants obtained from automated ab initio transition state theory-based master equation calculations and high-level thermochemical parametrizations. Broadband rotational spectroscopy offers quantitative and isomer-specific detection by which branching ratios of polar reaction products may be obtained. Using this technique, we observe and characterize products arising from H atom substitution reactions in the flash pyrolysis of acetone (CH 3 C(O)CH 3 ) at a nominal temperature of 1800 K. The major product observed is ketene (CH 2 CO). Minor products identified include acetaldehyde (CH 3 CHO), propyne (CH 3 CCH), propene (CH 2 CHCH 3 ), and water (HDO). Literature mechanisms for the pyrolysis of acetone do not adequately describe the minor products. The inclusion of a variety of substitution reactions, with rate constants and thermochemistry obtained from automated ab initio kinetics predictions and Active Thermochemical Tables analyses, demonstrates an important role for such processes. Furthermore, the pathway to acetaldehyde is shown to be a direct result of substitution of acetone's methyl group by a free H atom, while propene formation arises from OH substitution in the enol form of acetone by a free H atom.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Rapid Allylic 1,6 H-Atom Transfer in an Unsaturated Criegee Intermediate

Here, a novel allylic 1,6 hydrogen-atom-transfer mechanism is established through infrared activation of the 2-butenal oxide Criegee intermediate, resulting in very rapid unimolecular decay to hydroxyl (OH) radical products. A new precursor, Z/E -1,3-diiodobut-1-ene, is synthesized and photolyzed in the presence of oxygen to generate anew four-carbon Criegee intermediate with extended conjugation across the vinyl and carbonyl oxide groups that facilitates rapid allylic 1,6 H-atom transfer. A low-energy reaction pathway involving isomerization of 2-butenal oxide from a lower-energy ( tZZ ) conformer to a higher-energy ( cZZ ) conformer followed by 1,6 hydrogen transfer via a seven-membered ring transition state is predicted theoretically and shown experimentally to yield OH products. The low-lying ( tZZ ) conformer of 2-butenaloxide is identified based on computed an harmonic frequencies and intensities of its conformers. Experimental IR action spectra recorded in the fundamental CH stretch region with OH product detection by UV laser-induced fluorescence reveal a distinctive IR transition of the low-lying ( tZZ ) conformer at 2996 cm -1 that results in rapid unimolecular decay to OH products. Statistical RRKM calculations involving a combination of conformational isomerization and unimolecular decay via 1,6 H-transfer yield an effective decay rate k eff ( E ) on the order of 10 8 s -1 at ca. 3000 cm -1 in good accord with the experiment. Unimolecular decay proceeds with significant enhancement due to quantum mechanical tunneling. A rapid thermal decay rate of ca. 10 6 s -1 is predicted by master-equation modeling of 2-butenal oxide at 298 K, 1 bar. This novel unimolecular decay pathway is expected to increase the nonphotolytic production of OH radicals upon alkene ozonolysis in the troposphere.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Chemistry of Sugar Formation in the Gas Phase: Following the Activated Aldehyde

Sugars are produced by living organisms, and are required building blocks for life as we know it, which raises the foundational question of how sugars formed in a prebiotic environment. The abiotic formose reaction produces sugars from formaldehyde, but our understanding of its initiation step remains murky, with chemists invoking the concept of an “activated aldehyde” to seed this reaction. Singlet hydroxycarbenes, high-energy isomers of aldehydes, were recently reported to facilitate sugar formation under cold, nonaqueous conditions relevant to interstellar environments. Here, we generate singlet methylhydroxycarbene ( 1 CH 3 –C̈–OH) from the photodissociation of pyruvic acid and experimentally measure its gas-phase reaction with d 4 -acetaldehyde using multiplexed photoionization mass spectrometry. The C 4 H 4 D 4 O 2 isomer d 4 -acetoin is the sole product, which we kinetically link to the reactant CH 3 –C̈–OH, and attribute to a carbonyl-ene formation mechanism. We see no evidence of 3-hydroxybutanal, the C–H insertion product expected in carbene chemistry. Using automated exploration we calculate stationary points on the potential energy surface and report master equation rate coefficients from T = 20–600 K, providing quantitative kinetics of this fast reaction for use in chemical models. The prereactive complex in this reaction is stabilized by both hydrogen bonding and electrophilic carbene C═O interactions. These effects create a short-range dynamical bottleneck for the reaction besides the long- and midrange barrierless bottlenecks. Combined with recent reports of 1 HC̈OH production from methanol photodissociation and pyruvic acid production in cold irradiated ices, this work provides evidence that singlet hydroxycarbene + aldehyde chemistry is a feasible path to prebiotic sugar formation.

aldehydes↗
DOE OSTI2022

Reconfigurable quantum phononic circuits via piezo-acoustomechanical interactions

Abstract We show that piezoelectric strain actuation of acoustomechanical interactions can produce large phase velocity changes in an existing quantum phononic platform: aluminum nitride on suspended silicon. Using finite element analysis, we demonstrate a piezo-acoustomechanical phase shifter waveguide capable of producing ± π phase shifts for GHz frequency phonons in 10s of μm with 10s of volts applied. Then, using the phase shifter as a building block, we demonstrate several phononic integrated circuit elements useful for quantum information processing. In particular, we show how to construct programmable multi-mode interferometers for linear phononic processing and a dynamically reconfigurable phononic memory that can switch between an ultra-long-lifetime state and a state strongly coupled to its bus waveguide. From the master equation for the full open quantum system of the reconfigurable phononic memory, we show that it is possible to perform read and write operations with over 90% quantum state transfer fidelity for an exponentially decaying pulse.

Taylor, Jeffrey C. (ORCID:0000000215982436)↗
DOE OSTI2020

Theoretical study of the reaction mechanism and kinetics of the phenyl + propargyl association

Potential energy surface for the phenyl + propargyl radical recombination reaction has been studied at the CCSD(T)-F12/cc-pVTZ-f12//B3LYP/6-311G** level of theory for the closed-shell singlet species and at the triplet–singlet gap CASPT2/cc-pVTZ-CCSD(T)-F12/cc-pVTZ-f12//CASSCF/cc-pVTZ level of theory for the diradical species. High-pressure limit rate constants for the barrierless channels were evaluated with variable reaction coordinate transition state theory (VRC-TST). Rice–Ramsperger–Kassel–Marcus Master Equation (RRKM-ME) calculations have been performed to assess temperature- and pressure-dependent phenomenological rate constants and product branching ratios. Here, the entrance channels of the radical association reaction produce 3-phenyl-1-propyne and phenylallene which can further dissociate/isomerize into a variety of unimolecular and bimolecular products. Theoretical evidence is presented that, at combustion relevant conditions, the phenyl + propargyl recombination provides a feasible mechanism for the addition of a second five-member ring to the first six-member aromatic ring producing the prototype two-ring species indene and indenyl. Rate expressions for all important reaction channels in a broad range of temperatures and pressures have been generated for kinetic modeling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Mechanism and kinetics of the oxidation of 1,3-butadien-1-yl (n-C 4 H 5 ): a theoretical study

Here, ab initio CCSD(T)-F12/cc-pVTZ-f12//B3LYP/6-311G(d,p) calculations of the C 4 H 5 O 2 potential energy surface have been combined with Rice–Ramsperger–Kassel–Marcus Master Equation (RRKM-ME) calculations of temperature- and pressure-dependent rate constants and product branching ratios to unravel the mechanism and kinetics of the n-C 4 H 5 + O 2 reaction. The results indicate that the reaction is fast, with the total rate constant being in the range of 3.4–5.6 × 10 –11 cm 3 molecule –1 s –1 . The main products include 1-oxo-n-butadienyl + O and acrolein + HCO, with their cumulative yield exceeding 90% at temperatures above 1500 K. Two conformers of 1-oxo-n-butadienyl + O are formed via a simple mechanism of O 2 addition to the radical site of n-C 4 H 5 followed by the cleavage of the O–O bond proceeding via a van der Waals C 4 H 5 O∙∙∙O complex. Alternatively, the pathways leading to acrolein + HCO involve significant reorganization of the heavy-atom skeleton either via formal migration of one O atom to the opposite end of the molecule or its insertion into the C1–C2 bond. Not counting thermal stabilization of the initial peroxy adducts, which prevails at low temperatures and high pressures, all other products share a minor yield of under 5%. Rate constants for the significant reaction channels have been fitted to modified Arrhenius expressions and are proposed for kinetic modeling of the oxidation of aromatic molecules and 1,3-butadiene. As a secondary reaction, n-C 4 H 5 + O 2 can be a source for the formation of acrolein observed experimentally in oxidation of the phenyl radical at low combustion temperatures, whereas another significant (secondary) product of the C 6 H 5 + O 2 reaction, furan, could be formed through unimolecular decomposition of 1-oxo-n-butadienyl. Both the n-C 4 H 5 + O 2 reaction and unimolecular decomposition of its 1-oxo-n-butadienyl primary product are shown not to be a substantial source of ketene.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

The CH(X 2 Π) + H 2 O reaction: two transition state kinetics

he reaction of ground state methylidyne (CH) with water vapor (H 2 O) is theoretically re-investigated using high-level coupled cluster computations in combination with semi-classical transition state theory (SCTST) and two-dimensional master equation simulations. Insertion of CH into a H–O bond of H 2 O over a submerged barrier via a well-skipping mechanism yielding solely H and CH 2 O is characterized. The reaction kinetics is effectively determined by the formation of a pre-reaction van der Waals complex (PRC, HC—OH 2 ) and its subsequent isomerization to activated CH 2 OH in competition with PRC re-dissociation. The tunneling effects are found to be minor, while variational effects in the PRC → CH 2 OH step are negligible. The calculated rate coefficient k ( T ) is nearly pressure-independent, but strongly depends on temperature with pronounced down-up behavior: a high value of 2 x 10 -10 cm 3 s -1 at 50 K, followed by a fairly steep decrease down to 8 × 10 -12 cm 3 s -1 at 900 K, but increasing again to 5 × 10 -11 cm 3 s -1 at 3500 K. Over the T-range of this work, k ( T ) can be expressed as: k ( T , P = 0) = 2.31 × 10 -11 ( T /300 K) -1.615 exp(-38.45/ T ) cm 3 s -1 for T = 50-400 K k ( T , P = 0) = 1.15 × 10 -12 ( T /300 K) 0.8637 exp(892.6/ T ) cm 3 s -1 for T = 400-1000 K k ( T , P = 0) = 4.57 × 10 -15 ( T /300 K) 3.375 exp(3477.4/ T ) cm 3 s -1 for T = 1000–3500 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗