Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetic parameters”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

On the Rate Constant for NH 2 +HO 2 and Third-Body Collision Efficiencies for NH 2 +H(+M) and NH 2 +NH 2 (+M)

In low-temperature flash photolysis of NH 3 /O 2 /N 2 mixtures, the NH 2 consumption rate and the product distribution is controlled by the reactions NH 2 + HO 2 → products (R1), NH 2 + H (+M) → NH 3 (+M) (R2), and NH 2 + NH 2 (+M) → N 2 H 4 (+M) (R3). In the present work, published flash photolysis experiments by, among others, Cheskis and co-workers, are re-interpreted using recent direct measurements of NH 2 + H (+N 2 ) and NH 2 + NH 2 (+N 2 ) from Altinay and Macdonald. To facilitate analysis of the FP data, relative third-body collision efficiencies compared to N 2 for R2 and R3 were calculated for O 2 and NH 3 as well as for other selected molecules. We report results were in good agreement with the limited experimental data. Based on reported NH 2 decay rates in flash photolysis of NH 3 /O 2 /N 2 , a rate constant for NH 2 + HO 2 → NH 3 + O 2 (R1a) of $k_{1\text{a}}$ = 1.5(±0.5) × 10 14 cm 3 mol –1 s –1 at 295 K was derived. This value is higher than earlier determinations based on the FP results but in good agreement with recent theoretical work. Kinetic modeling of reported N 2 O yields indicates that NH 2 + HO 2 → H 2 NO + O (R1c) is competing with R1a, but perturbation experiments with addition of CH4 indicate that it is not a dominating channel. Measured HNO profiles indicate that this component is formed directly by NH 2 + HO 2 → HNO + H 2 O (R1b), but theoretical work indicates that R1b is only a minor channel. Based on this analysis, we estimate $k_{1\text{c}}$ = 2.5 × 10 13 cm 3 mol –1 s –1 and $k_{1\text{b}}$ = 2.5 × 10 12 cm 3 mol –1 s –1 at 295 K, with significant uncertainty margins.

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

Theoretical Study of the Phenoxy Radical Recombination with the O( 3 P) Atom, Phenyl plus Molecular Oxygen Revisited

Quantum chemical calculations of the C 6 H 5 O 2 potential energy surface (PES) were carried out to study the mechanism of the phenoxy + O( 3 P) and phenyl + O 2 reactions. CASPT2(15e,13o)/CBS//CASSCF(15e,13o)/DZP multi-reference calculations were utilized to map out the minimum energy path for the entrance channels of the phenoxy + O( 3 P) reaction. Stationary points on the C 6 H 5 O 2 PES were explored at the CCSD(T)-F12/cc-pVTZ-f12//B3LYP/6-311++G** level for the species with a single-reference character of the wave function and at the CASPT2(15e,13o)/CBS//B3LYP/6-311++G** level of theory for the species with a multi-reference character of the wave function. Conventional, variational, and variable reaction coordinate transition state theories were employed in Rice-Ramsperger-Kassel-Marcus Master Equation calculations to assess temperature- and pressure-dependent phenomenological rate constants and product branching ratios. Here, the main bimolecular product channels of the phenoxy + O( 3 P) reaction are concluded to be para/ortho-benzoquinone + H, 2,4-cyclopentadienone + HCO, and, at high temperatures, also phenyl + O 2 . The main bimolecular product channels of the phenyl + O 2 reaction include 2,4-cyclopentadienone + HCO at lower temperatures, and phenoxy + O( 3 P) at higher temperatures. Both for the phenoxy + O( 3 P) and phenyl + O 2 reactions, the collisional stabilization of peroxybenzene at low temperatures and high pressures competes with the bimolecular product channels.

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

On the Mechanism of Soot Nucleation. III. The Fate and Facility of the E-Bridge

Rotationally-excited dimerization of aromatic moieties is a mechanism proposed recently to explain the initial steps of soot particle inception in combustion and pyrolysis of hydrocarbons. The product of such dimerization, termed E-bridge, is an angled molecular structure composed of two aromatic rings sharing a common bond. The present study explores the immediate fate of the E-bridge. The performed theoretical analysis indicates that abstraction of a bridge H atom by a gaseous H leads to a rapid transformation of the angled to planar structure. The implications of this result is that the collisionally-activated E-bridge formation followed by its flattening effectively increases the size of “planar” aromatic precursors by combining two aromatic moieties with essentially collisional rates, instead of a slower “atom-by-atom” buildup. The faster growth speeds up PAH reaching a size when physical dimerization takes over. Furthermore, the dimerization can be further assisted by the biradicaloid valence structure of the flattened E-bridge.

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

Kinetics, Products, and Brown Carbon Formation by Aqueous-Phase Reactions of Glycolaldehyde with Atmospheric Amines and Ammonium Sulfate

Glycolaldehyde (GAld) is a C 2 water-soluble aldehyde produced during the atmospheric oxidation of isoprene and many other species and is commonly found in cloudwater. Previous work has established that glycolaldehyde evaporates more readily from drying aerosol droplets containing ammonium sulfate (AS) than does glyoxal, methylglyoxal, or hydroxyacetone, which implies that it does not oligomerize as quickly as these other species. Here, we report NMR measurements of glycolaldehyde’s aqueous-phase reactions with AS, methylamine, and glycine. Reaction rate constants are smaller than those of respective glyoxal and methylglyoxal reactions in the pH range of 3–6. In follow-up cloud chamber experiments, deliquesced glycine and AS seed particles were found to take up glycolaldehyde and methylamine and form brown carbon. At very high relative humidity, these changes were more than 2 orders of magnitude faster than predicted by our bulk liquid NMR kinetics measurements, suggesting that reactions involving surface-active species at crowded air–water interfaces may play an important role. The high-resolution liquid chromatography–electrospray ionization–mass spectrometric analysis of filter extracts of unprocessed AS + GAld seed particles identified sugar-like C 6 and C 12 GAld oligomers, including proposed product 3-deoxyglucosone, with and without modification by reactions with ammonia to diimine and imidazole forms. Chamber exposure to methylamine gas, cloud processing, and simulated sunlight increased the incorporation of both ammonia and methylamine into oligomers. Many C 4 –C 16 imidazole derivatives were detected in an extract of chamber-exposed aerosol along with a predominance of N-derivatized C 6 and C 12 glycolaldehyde oligomers, suggesting that GAld is capable of forming brown carbon SOA.

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

Resolving Discrepancies between State-of-the-Art Theory and Experiment for HO 2 + HO 2 via Multiscale Informatics

Recent high-level theoretical calculations predict a mild temperature dependence for HO 2 + HO 2 inconsistent with state-of-the-art experimental determinations that upheld the stronger temperature dependence observed in early experiments. Via MultiScale Informatics analysis of the theoretical and experimental data, we identified an alternative interpretation of the raw experimental data that uses HO 2 + HO 2 rate constants nearly identical to theoretical predictions---implying that the theoretical and experimental data are actually consistent, at least when considering the raw data from experimental studies. Here, similar analyses of typical signals from low-temperature experiments indicate that an HOOOOH intermediate---identified by recent theory but absent from earlier interpretations---yields modest effects that are smaller than, but may have contributed to, the scatter in data among different experiments. More generally, the findings demonstrate that modern chemical theories and experiments have progressed to a point where meaningful comparison requires joint consideration of their data simultaneously.

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

Probing High-Temperature Amine Chemistry: Is the Reaction NH 3 + NH 2 ⇌ N 2 H 3 + H 2 Important?

The reaction NH 3 + NH 2 ⇌ N 2 H 3 + H 2 (R1) has been identified as a key step to explain experimental results for pyrolysis and oxidation of ammonia. However, no direct experimental or theoretical evidence for the reaction has been reported. In the present work, the reaction was studied by ab initio theory and by reinterpretation of experimental data. We could not locate a transition state for R1 occurring as a direct process, but alternative mechanisms yielded an upper bound to k 1 of 1.5 × 10 13 exp(–58.9 kcal mol –1 /RT) cm 3 mol –1 s –1 over 1000–2500 K, several orders of magnitude lower than values applied in modeling. Consistent with the theoretical work, re-evaluation of NH 3 pyrolysis data supported a very low value of k 1 . However, this finding opens up a novel unresolved issue. Current kinetic models cannot capture the NH 3 oxidation behavior in a number of laminar flow reactor and jet-stirred reactor experiments without adopting an improbably high value for k 1 . Furthermore, important oxidation steps might be underestimated or missing from mechanisms.

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

Re-Examination of the N 2 O + O Reaction

The reaction of N 2 O with O is a key step in consumption of nitrous oxide in thermal processes. It has two product channels, NO + NO (R2) and N 2 + O 2 (R3). The rate constant for R2 has been measured both in the forward and the reverse direction at elevated temperature and is well established. However, the rate constant for the N 2 + O 2 channel (R3) has been difficult to quantify and has significant error limits. Here, the direct reaction on the triplet surface has a barrier of around 40 kcal mol –1 , and it is too slow for the N 2 + O 2 channel to have any practical significance. Recently, Pham and Lin (2022) suggested an alternative low activation energy reaction path that involves intersystem crossing and reaction on the singlet surface. In the present work, we re-examined a wide range of experiments relevant for the N 2 O + O reaction through kinetic modeling, paying attention to the impact of artifacts such as impurities and surface reactions. Experimental results from shock tubes and batch reactors on the final NO yield in N 2 O decomposition, covering temperatures of 973–2200 K and pressures of 0.013–11.5 atm, support k 3 ~ 0, consistent with the high activation energy for reaction on the triplet surface and a low probability of ISC.

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

High-Pressure Rate Rules for Ether Alkylperoxy Radical Isomerization

The first isomerization reaction of an alkylperoxy (RO 2 ) radical holds significant importance in low-temperature oxidation, as it governs the branching ratios of the hydroperoxyalkyl (QOOH) radicals, which influence the competition between the chain-propagation and chain-branching reactions. In this study, we systematically calculated high-pressure rate rules for the RO 2 isomerization reaction of monoethers, exploring 5-, 6-, 7-, and 8-membered ring transition states. Primary, secondary, and tertiary carbon sites, where both the abstracting peroxy group and the abstracted hydrogen are located, were considered, with particular emphasis on distinguishing between secondary carbons adjacent (alpha) and nonadjacent to the ether functional group. Using the G4//B3LYP/6-311++G(2df,2pd) level of theory and the transition state theory, we estimated the rate constants and the Arrhenius coefficient for over 120 possible isomerization reactions. We examined the effect of ring size and ring atoms, revealing that 6- and 7-membered ring isomerizations were generally the fastest. The impact of the ether functional group on transition states was investigated by comparing reactions with identical ring size, peroxy, and radical positions, but with the ether functional group positioned either outside (i.e., out) or inside (i.e., in) the transition state ring, leading to differences in the rate constants. When comparing to analogous alkane rate constants, differences of up to an order of magnitude were observed, underscoring the need for caution when assigning rate rules by analogy. We applied our rate constants in the di-iso-butyl ether kinetic model and evaluated their influence on low-temperature chemistry finding that they altered the branching ratios by up to a factor of 9, highlighting the significance of site-specific rate constants for more accurate low-temperature modeling.

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 OSTI2025

The Reaction Rates of Amidogen and Ammonia with Nitrous Oxide: Implications for Combustion Mechanisms

Pulsed laser photolysis experiments with laser-induced fluorescence detection of NH 2 set an upper limit to the rate constant for reaction with N 2 O of k < 1 × 10 −15 cm 3 molecule −1 s −1 at 513 K. Computations were based on geometries and anharmonic frequency analysis (B2PLYP-D3/cc-pVTZ) followed by coupled cluster calculations extrapolated to the infinite basis set limit, with corrections for core−valence electron correlation, scalar relativistic effects, and correlation up to CCSDT(Q). Species that showed multireference character were quantified with MRCI(7,7)+Q/cc-pVTZ theory. Rate constants were obtained for the dominant product channel H 2 NN + NO, along with HN 3 + OH, H 2 NO + N 2 , NNH + HNO, ON(NH) 2 and HNNH + NO. The last channel is slow even at 2500 K, contrary to an early empirical estimate and confirming recent suggestions. Modeling of literature experiments on oxidation of NH 3 by N 2 O shows that all channels are too slow to make a significant impact on the loss of N 2 O in ammonia flames. Similarly, the direct NH 3 + N 2 O reaction is found to be negligibly slow.

Ammonia↗
DOE OSTI2026

Restricted Partition Function Semiclassical Transition State Theory (RPF-SCTST): Applications to Reactions Involving Hydrogen Cyanide, Hydrogen Peroxide, and Formaldehyde Oxide

We demonstrate the efficiency of the numerical calculation of thermal semiclassical transition state theory (SCTST) rates across several representative chemical reactions using the restricted partition function (RPF)─viz a function that depends only on the imaginary action associated with a reaction. Here, we treat the potential energy surfaces (PESs) through fourth order expansions around the saddle point and integrate the Hamiltonian up to second order in employing vibrational perturbation theory. We apply this formalism to uni- and bimolecular reactions with different─though relatively high─barrier heights to uncover the influence of the barrier properties, minimal energies, and separability of the rotational motion on rate constants and tunneling corrections. Although all modes are coupled within the RPF, we found in our numerical examples that the rotational component in the absence of significant rotational distortions can be separated from the vibrational contributions. Moreover, a classical treatment of the rotational contribution is adequate over the temperature range from approximately 100 to 1000 K. We also found that the choice of DFT basis set can lead to variations in rate constants of up to an order of magnitude. Lastly, different schemes for counting eligible energy levels result in rate constants that differ by no more than a factor of 2, with the remaining discrepancies attributed to the treatment of near-convergent levels in perturbation theory.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Parsimonious Potential Energy Surface Expansions Using Dictionary Learning with Multipass Greedy Selection

Potential energy surfaces fit with basis set expansions have been shown to provide accurate representations of electronic energies and have enabled a variety of high-accuracy dynamics, kinetics, and spectroscopy applications. The number of terms in these expansions scales poorly with system size, a drawback that challenges their use for systems with more than similar to 10 atoms. A solution is presented here using dictionary learning. Subsets of the full set of conventional basis functions are optimized using a newly developed multipass greedy regression method inspired by forward and backward selection methods from the statistics, signal processing, and machine learning literatures. Here, the optimized representations have accuracies comparable to the full set but are 1 or more orders of magnitude smaller, and notably, the number of terms in the optimized multipass greedy expansions scales approximately linearly with the number of atoms.

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

New Dimension in Ab Initio Electronic Structure Theory: Temperature, Pressure, and Chemical Potential

Ab initio electronic structure theory has transformed gas-phase molecular science with its predictive ability. In the attempt to bring such predictive ability to macroscopic systems and condensed matter, the theory must integrate quantum mechanics with statistical thermodynamics, so that thermodynamic functions such as free energy, internal energy, entropy, and chemical potentials are computed as functions of temperature in a systematically converging series of approximations. Here, a general, versatile strategy of elevating ab initio electronic structure theory to nonzero temperatures is introduced and discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Intertwined Topological and Magnetic Orders in Atomically Thin Chern Insulator MnBi 2 Te 4

The interplay between band topology and magnetic order plays a key role in quantum states of matter. MnBi 2 Te 4 , a van der Waals magnet, has recently emerged as an exciting platform for exploring Chern insulator physics. Its layered antiferromagnetic order was predicted to enable even-odd layer-number dependent topological states, supported by promising edge transport measurements. Furthermore, it becomes a Chern insulator when all spins are aligned by an applied magnetic field. However, the evolution of the bulk electronic structure as the magnetic state is continuously tuned and its dependence on layer number remains unexplored. Here, employing multimodal probes, we establish oneto-one correspondence between bulk electronic structure, magnetic state, topological order, and layer thickness in atomically thin MnBi 2 Te 4 devices. As the magnetic state is tuned through the canted magnetic phase, we observe a band crossing, i.e., the closing and reopening of the bulk bandgap, corresponding to the concurrent topological phase transition. Surprisingly, we find that the even- and odd-layer number devices exhibit a similar topological phase transition coupled to magnetic states, distinct from recent theoretical and experimental reports. Our findings shed new light on the interplay between band topology and magnetic order in this newly discovered topological magnet and validate the band crossing with concurrent measurements of topological invariant in a continuously tuned topological phase transition.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Relationship between the Nature of Monovalent Cations and Charge Recombination in Metal Halide Perovskites

We report charge recombination rates in lead iodide perovskites with alloys of formamidinium (FA) and methylammonium (MA) cations, FA x MA 1–x PbI 3 , 0 ≤ x ≤ 0.9, prepared through ion exchange with minimal differences in morphology and charge trap density. Our results show that the trap-mediated recombination rate increases over 2 orders of magnitude with decreasing MA content. These results are consistent with a proposed mechanism that MA molecules are more effective in screening charged defects. In contrast, band-to-band charge recombination rates are minimized at an intermediate alloy composition. Here, these findings reveal that the monovalent cations impact recombination processes through different mechanisms.

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

Thermal Hydroquinone Oxidation on Co/N-doped Carbon Proceeds by a Band-Mediated Electrochemical Mechanism

Molecules catalyze aerobic oxidation reactions via redox cycling at the metal center to effect sequential activation of O 2 and the substrate. Metal surfaces can catalyze the same transformations by coupling independent half-reactions for oxygen reduction and substrate oxidation mediated via the exchange of band-electrons. Metal/nitrogen-doped carbons (MNCs) are promising catalysts for aerobic oxidation that consist of molecule-like active sites embedded in a conductive carbon host. Owing to their combining molecular and metallic features, it remains unclear whether they proceed via the sequential redox cycling pathways of molecules or the band-mediated pathways of metals. Herein, we simultaneously track the potential of the catalyst and the rate of the turnover of aerobic hydroquinone oxidation on a cobalt-based MNC catalyst contacted to a catalytically inactive carbon electrode. By comparing operando measurements of rate and potential with the current-voltage behavior of each constituent half-reaction under identical conditions, we show that these molecular materials can display the band-mediated reaction mechanisms of extended metallic solids. Further, we show that the action of these band-mediated mechanisms explains fractional reaction orders in both oxygen and hydroquinone, the time evolution of catalyst potential and rate, and the dependence of rate on overall reaction free energy. Selective poisoning experiments suggest that oxygen reduction proceeds at cobalt sites, whereas quinone oxidation proceeds at native oxidic defects on the carbon. These findings highlight that molecule-like active sites can take advantage of band-mediated mechanisms when coupled to conductive hosts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗