Engineering topics

Jasper, Ahren W.

Publications and source records attributed to Jasper, Ahren W..

At least 19 records

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 OSTI2024

Dynamics of HCN, HNC, and HNCO Formation in the 193 nm Photodissociation of Formamide

Formamide (NH 2 CHO) is the simplest molecule containing a peptide linkage [-NH-C($=$O)-], and it plays an essential role in the study of prebiotic chemistry. Exposure to UV irradiation allows formamide to decompose and act as a prebiotic feedstock in the formation of nucleobases and other necessary starting materials. The photodissociation mechanism of gaseous formamide at 193 nm is studied using (a) chirped-pulse Fourier transform millimeter-wave spectroscopy in the 260-290 GHz spectral region in a room-temperature flow-tube reactor at 1 mu bar pressure, (b) a combination of electronic structure theory, transition state theory, and quasiclassical trajectories, and (c) the Active Thermochemical Tables. The HCN and HNC photoproducts of hydrogenated (NH 2 CHO) and deuterated (NH 2 CDO and ND 2 CHO) formamide precursors are examined to gain insight into the photodissociation mechanism. Here, the theoretical investigation has characterized the main pathway leading to each of the HCN/HNC isomers from the precursor isotopologues. The theoretical branching ratio [HNC]/[HCN] = 2.1 for nascent photofragments agrees with the experiment. The effect of the postphotolysis HNC ↔ HCN isomerization on the [HNC]/[HCN] ratio is predicted. We report the experimental branching ratio [HNCO]: ([HNC] + [HCN]) = 12 ± 3 and propose that most of HNCO originates from dissociation on the S-1 electronic state of formamide.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Tracking the reaction networks of acetaldehyde oxide and glyoxal oxide Criegee intermediates in the ozone-assisted oxidation reaction of crotonaldehyde

The reaction of unsaturated compounds with ozone (O 3 ) is recognized to lead to the formation of Criegee intermediates (CIs), which play a key role in controlling the atmospheric budget of hydroxyl radicals and secondary organic aerosols. The reaction network of two CIs with different functionality, i.e. acetaldehyde oxide (CH 3 CHOO) and glyoxal oxide (CHOCHOO) formed in the ozone-assisted oxidation reaction of crotanaldehyde (CA), is investigated over a temperature range between 390 K and 840 K in an atmospheric pressure jet-stirred reactor (JSR) at a residence time of 1.3 s, stoichiometry of 0.5 with a mixture of 1% crotonaldehyde, 10% O 2 , at an fixed ozone concentration of 1000 ppm and 89% Ar dilution. Molecular-beam mass spectrometry in conjunction with single photon tunable synchrotron vacuum-ultraviolet (VUV) radiation is used to identify elusive intermediates by means of experimental photoionization energy scans and ab initio threshold energy calculations for isomer identification. Addition of ozone (1000 ppm) is observed to trigger the oxidation of CA already at 390 K, which is below the temperature where the oxidation reaction of CA was observed in the absence of ozone. The observed CA + O 3 product, C 4 H 6 O 4 , is found to be linked to a ketohydroperoxide (2-hydroperoxy-3-oxobutanal) resulting from the isomerization of the primary ozonide. Products corresponding to the CIs uni- and bi-molecular reactions were observed and identified. A network of CI reactions is identified in the temperature region below 600 K, characterized by CIs bimolecular reactions with species like aldehydes, i.e., formaldehyde, acetaldehyde, and crotonaldehyde and alkenes, i.e., ethene and propene. The region below 600 K is also characterized by the formation of important amounts of typical low-temperature oxidation products, such as hydrogen peroxide (H 2 O 2 ), methyl hydroperoxide (CH 3 OOH), and ethyl hydroperoxide (C 2 H 5 OOH). Detection of additional oxygenated species such as alcohols, ketene, and aldehydes are indicative of multiple active oxidation routes. This study provides important information about the initial step involved in the CIs assisted oligomerization reactions in complex reactive environments where CIs with different functionalities are reacting simultaneously. It provides new mechanistic insights into ozone-assisted oxidation reactions of unsaturated aldehydes, which is critical for the development of improved atmospheric and combustion kinetics models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Initiation and Carbene Induced Radical Chain Reactions in CH 2 F 2 Pyrolysis

High temperature dissociations of organic molecules typically involve a competition between radical and molecular processes. In this work, we use a modeling, experiment, theory (MET) framework to characterize the high temperature thermal dissociation of CH 2 F 2 , a flammable hydrofluorocarbon (HFC) that finds widespread use as a refrigerant. Initiation in CH 2 F 2 proceeds via a molecular elimination channel; CH 2 F 2 →CHF+HF. Here we show that the subsequent self-reactions of the singlet carbene, CHF, are fast multichannel processes and a facile source of radicals that initiate rapid chain propagation reactions. These have a marked influence on the decomposition kinetics of CH 2 F 2 . The inclusion of these reactions brings the simulations into better agreement with the present and literature experiments. Additionally, flame simulations indicate that inclusion of the CHF+CHF multichannel reaction leads to a noticeable enhancement in predictions of laminar flame speeds, a key parameter that is used to determine the flammability of a refrigerant.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Quasi-Classical Trajectory Calculation of Rate Constants Using an Ab Initio Trained Machine Learning Model (aML-MD) with Multifidelity Data

Machine learning (ML) provides a great opportunity for the construction of models with improved accuracy in classical molecular dynamics (MD). However, the accuracy of a ML trained model is limited by the quality and quantity of the training data. Generating large sets of accurate ab initio training data can require significant computational resources. Furthermore, inconsistent or incompatible data with different accuracies obtained using different methods may lead to biased or unreliable ML models that do not accurately represent the underlying physics. Recently, transfer learning showed its potential for avoiding these problems as well as for improving the accuracy, efficiency, and generalization of ML models using multifidelity data. In this work, ab initio trained ML-based MD (aML-MD) models are developed through transfer learning using DFT and multireference data from multiple sources with varying accuracy within the Deep Potential MD framework. Further, the accuracy of the force field is demonstrated by calculating rate constants for the H + HO 2 → H 2 + 3 O 2 reaction using quasi-classical trajectories. We show that the aML-MD model with transfer learning can accurately predict the rate constants while reducing the computational cost by more than five times compared to the use of more expensive quantum chemistry training data sets. Hence, the aML-MD model with transfer learning shows great potential in using multifidelity data to reduce the computational cost involved in generating the training set for these potentials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Isomer-resolved unimolecular dynamics of the hydroperoxyalkyl intermediate (•QOOH) in cyclohexane oxidation

The oxidation of cycloalkanes is important in the combustion of transportation fuels and in atmospheric secondary organic aerosol formation. A transient carbon-centered radical intermediate (•QOOH) in the oxidation of cyclohexane is identified through its infrared fingerprint and time- and energy-resolved unimolecular dissociation dynamics to hydroxyl (OH) radical and bicyclic ether products. Although the cyclohexyl ring structure leads to three nearly degenerate •QOOH isomers (β-, γ-, and δ-QOOH), their transition state (TS) barriers to OH products are predicted to differ considerably. Selective characterization of the β-QOOH isomer is achieved at excitation energies associated with the lowest TS barrier, resulting in rapid unimolecular decay to OH products that are detected. A benchmarking approach is employed for the calculation of high-accuracy stationary point energies, in particular TS barriers, for cyclohexane oxidation (C 6 H 11 O 2 ), building on higher-level reference calculations for the smaller ethane oxidation (C 2 H 5 O 2 ) system. The isomer-specific characterization of β-QOOH is validated by comparison of experimental OH product appearance rates with computed statistical microcanonical rates, including significant heavy-atom tunneling, at energies in the vicinity of the TS barrier. Master-equation modeling is utilized to extend the results to thermal unimolecular decay rate constants at temperatures and pressures relevant to cyclohexane combustion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Phenalenyl growth reactions and implications for prenucleation chemistry of aromatics in flames

The energetics and kinetics of phenalene and phenalenyl growth reactions were studied theoretically. Rate constants of phenalene and phenalenyl H-abstraction and C 2 H 2 addition to the formed radicals were evaluated through quantum-chemical and rate-theory calculations. The obtained values, assigned to all π radicals, were tested in deterministic and kinetic Monte Carlo simulations of aromatics growth under conditions of laminar premixed flames. Kekulé and non-Kekulé structures of the polycyclic aromatic hydrocarbons (PAHs) evolving in the stochastic simulations were identified by on-the-fly constrained optimization. The numerical results demonstrated an increased PAH growth and qualitatively reproduced experimental observations of Homann and co-workers of non-decaying PAH concentrations with nearly equal abundances of even and odd carbon-atom PAHs. The analysis revealed that the PAH growth proceeds via alternating and sterically diverse acetylene and methyl HACA additions. Furthermore, the rapid and diverse spreading in the PAH population supports a nucleation model as PAH dimerization, assisted by the non-equilibrium phenomena, forming planar aromatics first and then transitioning to the PAH–PAH stacking with size.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Theoretical examination of nuclear spin diffusion in light-induced spin coherences in photosystem I

Microsecond coherence times are predicted in spin dynamics calculations of a donor–acceptor electron spin pair P + A 1A - created upon light excitation of photosystem I (PSI). The effects of nuclear spin diffusion (NSD) due to the anisotropic protein environment on the predicted coherence times T M are studied. Closely positioned pairs and triples of protons located 5–8 Å from the electron spin are shown to largely control T M . Additionally, knowledge of the crystal structure of PSI allows for spin dynamics calculations where specific cofactors and amino acid residues are removed or replaced and the identification of anisotropic environmental features controlling electronic decoherence. Finally, we show that NSD alone cannot explain the >3× shorter experimentally observed coherence times and suggest that methyl groups at key protein sites may explain this discrepancy.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2024

An experimental, theoretical, and kinetic modeling study of post-flame oxidation of ammonia

The post -flame oxidation rate of ammonia was investigated in a novel atmospheric pressure flow reactor at temperatures of 1280 ± 16 K and as a function of residence time and mixture composition (1-10% O 2 , dry and moist). The experimental results, as well as selected data from literature, were analyzed using an updated detailed chemical kinetic model. The medium temperature, very lean conditions enhance the importance of reactions of the nitroxyl (HNO) intermediate. High-level theory was used to calculate the rate constant for HNO + NH 2 , indicating that this step is significantly faster than values used in literature. Furthermore, a trajectory based approach was used to determine collision efficiencies for selected bath gases for HNO + M. The experimental results show that the NH 3 oxidation rate increases with temperature and O 2 concentration, while the presence of water vapor slightly inhibits reaction. Formation of NO and N 2 O was strongly promoted at higher levels of O 2 . Modeling results agreed well with the measurements, except at the lowest level of O 2 . The predicted oxidation rate of NH 3 was shown to result from a delicate balance between chain branching and terminating steps involving NH 2 , H 2 NO, and HNO. Recent theoretical work on reactions of these species by Klippenstein and coworkers and Stagni et al. was instrumental in improving modeling predictions. After initiation, NO reached a pseudo -steady-state level, where the pathways to NO were largely balanced by the NH 2 + NO reaction. Nitric oxide was partly oxidized to NO 2 , with the NH 2 + NO 2 reaction responsible for most of the N 2 O formation. Novelty and significance statement: This study provides the first detailed kinetic analysis of the lean postflame oxidation of ammonia, based on time -resolved flow reactor data in a novel reactor. In addition to the post -flame oxidation rate of ammonia, data for formation of NO and N 2 O were compared with modeling predictions. The medium temperature, very lean conditions enhance the importance of reactions of the HNO and H 2 NO intermediates. Inclusion in the model of results from recent high-level theoretical work, including present calculations for HNO + NH 2 and HNO + M, was crucial for capturing the observed behavior. It is argued that the post -initiation steady-state NH 3 oxidation rates constitute important data for model validation, along with ignition delays and laminar flame speeds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Can Third-Body Stabilization of Bimolecular Collision Complexes in Cold Molecular Clouds Happen?

For more than half a century, networks of radiative association, dissociative recombination, and bimolecular reactions have been postulated to drive the low-temperature chemistry of cold molecular clouds. Third-body stabilizations of collision complexes have been assumed to be 'irrelevant' due to short lifetime of such complexes. Here, we conduct crossed molecular beam studies of ground state atomic silicon with diacetylene in combination with electronic structure calculations and microcanonical kinetics models operating under cold molecular cloud conditions. Our combined experimental, electronic structure, and microcanonical kinetics modelling investigations provide compelling evidence that three-body collisions of molecular hydrogen with long-lived reaction intermediates accessed through intersystem crossing are prevalent deep inside molecular clouds. This concept might be exportable to reactions involving polycyclic aromatic hydrocarbons thus affording a versatile machinery to complex organics via third-body stabilizations of bimolecular collision complexes deep inside cold molecular clouds.

astrochemistry↗

DOE OSTI2023

Coherences of Photoinduced Electron Spin Qubit Pair States in Photosystem I

Here, this publication presents the first comprehensive experimental study of electron spin coherences in photosynthetic reaction center proteins, specifically focusing on photosystem I (PSI). The ultrafast electron transfer in PSI generates spin-correlated radical pairs (SCRPs), which are entangled spin pairs formed in well-defined spin states (Bell states). Since their discovery in our group in the 1980s, SCRPs have been extensively used to enhance our understanding of structure–function relationships in photosynthetic proteins. More recently, SCRPs have been utilized as tools for quantum sensing. Electron spin decoherence poses a significant challenge in realizing practical applications of electron spin qubits, particularly the creation of quantum entanglement between multiple electron spins. This work is focused on the systematic characterization of decoherence in SCRPs of PSI. These decoherence times were measured as electron spin echo decay times, termed phase memory times (T M ), at various temperatures. Decoherence was recorded on both transient SCRP states P 700 + A 1 – and thermalized states. Our study reveals that T M exhibits minimal dependence on the biological species, biochemical treatment, and paramagnetic species. The analysis indicates that nuclear spin diffusion and instantaneous diffusion mechanisms alone cannot explain the observed decoherence. As a plausible explanation we discuss the assumption that the low-temperature dynamics of methyl groups in the protein surrounding the unpaired electron spin centers is the main factor governing the loss of the spin coherence in PSI.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Direct observation of coherence transfer and rotational-to-vibrational energy exchange in optically centrifuged CO2 super-rotors

Abstract Optical centrifuges are laser-based molecular traps that can rotationally accelerate molecules to energies rivalling or exceeding molecular bond energies. Here we report time and frequency-resolved ultrafast coherent Raman measurements of optically centrifuged CO 2 at 380 Torr spun to energies beyond its bond dissociation energy of 5.5 eV ( J max = 364, E rot = 6.14 eV, E rot / k B = 71, 200 K). The entire rotational ladder from J = 24 to J = 364 was resolved simultaneously which enabled a more accurate measurement of the centrifugal distortion constants for CO 2 . Remarkably, coherence transfer was directly observed, and time-resolved, during the field-free relaxation of the trap as rotational energy flowed into bending-mode vibrational excitation. Vibrationally excited CO 2 ( ν 2 > 3) was observed in the time-resolved spectra to populate after 3 mean collision times as a result of rotational-to-vibrational (R-V) energy transfer. Trajectory simulations show an optimal range of J for R-V energy transfer. Dephasing rates for molecules rotating up to 5.5 times during one collision were quantified. Very slow decays of the vibrational hot band rotational coherences suggest that they are sustained by coherence transfer and line mixing.

74 ATOMIC AND MOLECULAR PHYSICS↗

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

The P( 4 S) + NH( 3 Σ – ) and N( 4 S) + PH( 3 Σ – )reactions as sources of interstellar phosphorus nitride

Phosphorus nitride (PN) is believed to be one of the major reservoirs of phosphorus in the interstellar medium (ISM). For this reason, understanding which reactions produce PN in space and predicting their rate coefficients is important for modelling the relative abundances of P-bearing species and clarifying the role of phosphorus in astrochemistry. Here, in this work, we explore the potential energy surfaces of the $\textrm{P}(^4\textrm{S}) + \textrm{NH}(^3\Sigma^-)$ and $\textrm{N}(^4\textrm{S}) + \textrm{PH}(^3\Sigma^-)$ reactions and the formation of $\textrm{H}(^2\textrm{S}) + \textrm{PN}(^1\Sigma^+)$ through high accuracy ab initio calculations and the variable reaction coordinate transition state theory (VRC-TST). We found that both reactions proceed without an activation barrier and with similar rate coefficients that can be described by a modified Arrhenius equation ( $k(T)=\alpha\!\left( T/300 \right)^{\beta} \exp\!{(\!-\!\gamma/T)})$ with $\alpha=0.93\times 10^{-10}\rm cm^3\,s^{-1}$ , $\beta=-0.18$ and $\gamma=0.24\, \rm K$ for the $\textrm{P} + \textrm{NH} \longrightarrow \textrm{H} + \textrm{PN}$ reaction and $\alpha=0.88\times 10^{-10}\rm cm^3\,s^{-1}$ , $\beta=-0.18$ and $\gamma=1.01\, \rm K$ for the $\textrm{N} + \textrm{PH} \longrightarrow \textrm{H} + \textrm{PN}$ one. Both reactions are expected to be relevant for modelling PN abundances even in the cold environments of the ISM. Given the abundance of hydrogen in space, we have also predicted rate coefficients for the destruction of PN via H + PN collisions.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

The role of energy transfer and competing bimolecular reactions in characterizing the unimolecular dissociations of allylic radicals

Resonantly stabilized radicals (RSR's) such as allylic radicals typically have bond energies exceeding 45 kcal/mol and are therefore reasonably long-lived up to quite high temperatures. Consequently, kinetics of the unimolecular dissociations and bimolecular reactions of these long-lived RSR's assume significance in high temperature chemically reacting systems such as in flames. Our recent analyses (J. Cho et al., Proc. Comb. Inst., 2022) of the uncertainty propagated by the energy transfer parameters for 1-methylallyl (1MA) dissociation to the laminar flame speed, indicates an intricate coupling between the kinetics of 1MA dissociation (chain propagation) and its reaction with H-atoms (chain-termination). Here, this work extends such analyses to other C 3 – C 5 allylic radicals: allyl, 2-methylallyl, 1,1-dimethylallyl, 1,2-dimethylallyl, and 1,3-dimethylallyl. Theoretical kinetics for the C 3 – C 4 allylic radicals were taken from prior literature studies (J.A. Miller et al., J. Phys. Chem., 2008; J. Cho et al., Proc. Comb. Inst., 2022; R.S. Tranter et al., Proc. Comb. Inst., 2017). Potential energy surfaces were newly computed for the dissociations of the larger C5 allylic radicals and the reactions of these radicals with H-atoms. These PES's were then used to calculate rate constants for the relevant unimolecular and bimolecular reactions. Energy transfer has a crucial role to play in these pressure-dependent reactions and therefore the present kinetics calculations relied on a-priori predictions of collisional energy transfer parameters to best characterize the branching between chain propagating unimolecular dissociations of these allylic radicals and chain terminating bimolecular reactions of these RSR's with H-atoms. Prompt dissociations are now established to be a universal feature of free radical dissociations. Here we have also characterized the propensity for prompt dissociations in these more-stable allylic RSR's and their role in combustion simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Molecular weight growth by the phenyl + cyclopentadienyl reaction: Well-skipping, ring-opening, and dissociation

Radical-radical reactions contribute to molecular weight growth of polycyclic aromatic hydrocarbons (PAHs), but their role is not well understood. The phenyl (Ph, C 6 H 5 ) + cyclopentadienyl (C 5 H 5 ) reaction has not been studied but can serve as a surrogate reaction for many reactions of larger aryl-σ radicals with resonance-stabilized π-radicals. The adducts of these reactions easily lose an H atom to produce a new π-radical. Here, through a combination of dilute flash pyrolysis experiments, theory, and simulation, we find that for the title reaction, substantial Ph—C 5 H 4 + H is formed, dominantly by a well-skipping pathway, at conditions around 30 Torr and 1400 K. This radical is thermally stable at moderate temperatures but decomposes at higher temperatures following ring opening of the 5-member ring. Products corresponding to the loss of H, C 2 H 2 , and C 3 H 3 are observed experimentally. Calculation of the C 11 H 9 potential energy surface affirms that these are expected decomposition products and reveals multiple possible pathways.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

The role of radical-radical chain-propagating pathways in the phenyl + propargyl reaction

Well-skipping radical-radical reactions can provide a chain-propagating pathway for formation of polycyclic radicals implicated in soot inception. Here we use controlled pyrolysis in a microreactor to isolate and examine the role of well-skipping channels in the phenyl (C 6 H 5 ) + propargyl (C 3 H 3 ) radical-radical reaction at temperatures of 800-1600 K and pressures near 25 Torr. The temperature and concentration dependence of the closed-shell (C 9 H 8 ) and radical (C 9 H 7 ) products are observed using electron-ionization mass spectrometry. The flow in the reactor is simulated using a boundary layer model employing a chemical mechanism based on recent rate coefficient calculations. Further, comparison between simulation and experiment shows reasonable agreement, within a factor of 3, while suggesting possible improvements to the model. In contrast, eliminating the well-skipping reactions from the chemistry mechanism causes a much larger discrepancy between simulation and experiment in the temperature dependence of the radical concentration, revealing that the well-skipping pathways, especially to form indenyl radical, are significant at temperatures of 1200 K and higher. While most C 9 H 7 forms by well-skipping at 25 Torr, an additional simulation indicates that the wellskipping channels only contribute around 3% of the C 9 H x yield at atmospheric pressure, thus indicating a negligible role of the well-skipping pathways at atmospheric and higher pressures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Inefficient intramolecular vibrational energy redistribution for the H + HO2 reaction and negative internal energy dependence for its rate constant

Quasiclassical trajectories (QCT) and newly constructed global potential energy surfaces are used to compute thermal and nonthermal rate constants for the H + HO 2 reaction. The thermal QCTs rate constants are up to 50% smaller than transition state theory (TST) rate constants based on the same level of electronic structure theory. This reduction is demonstrated to result from inefficient intramolecular vibrational energy redistribution (IVR) in the transient H 2 O 2 well, with a significant fraction of trajectories that reach the H 2 O 2 well promptly dissociating back to reactants instead of via the heavily statistically favored 2OH channel. The nonstatistical reduction factor, κ IVR , that quantifies this effect is shown to increase in importance with temperature, with κ IVR = 0.81 at 300 K and 0.47 at 2500 K. Finally, we show that inefficient IVR causes H + HO 2 rate constants mediated by H 2 O 2 to depend inversely on the initial vibrational excitation of HO 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗