Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transition state”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Reveal, A General Reverse Engineering Algorithm for Inference of Genetic Network Architectures

Given the immanent gene expression mapping covering whole genomes during development, health and disease, we seek computational methods to maximize functional inference from such large data sets. Is it possible, in principle, to completely infer a complex regulatory network architecture from input/output patterns of its variables? We investigated this possibility using binary models of genetic networks. Trajectories, or state transition tables of Boolean nets, resemble time series of gene expression. By systematically analyzing the mutual information between input states and output states, one is able to infer the sets of input elements controlling each element or gene in the network. This process is unequivocal and exact for complete state transition tables. We implemented this REVerse Engineering ALgorithm (REVEAL) in a C program, and found the problem to be tractable within the conditions tested so far. For n = 50 (elements) and k = 3 (inputs per element), the analysis of incomplete state transition tables (100 state transition pairs out of a possible 10(exp 15)) reliably produced the original rule and wiring sets. While this study is limited to synchronous Boolean networks, the algorithm is generalizable to include multi-state models, essentially allowing direct application to realistic biological data sets. The ability to adequately solve the inverse problem may enable in-depth analysis of complex dynamic systems in biology and other fields.

Liang, Shoudan↗
DOE OSTI2020

Reactivity descriptors in acid catalysis: acid strength, proton affinity and host–guest interactions

Brønsted acids mediate chemical transformations via proton transfer to bound species and interactions between the conjugate anion and bound cationic intermediates and transition states that are also stabilized by van der Waals forces within voids of molecular dimensions in inorganic hosts. This Feature Article describes the relevant descriptors of reactivity in terms of the properties of acids and molecules that determine their ability to donate and accept protons and to reorganize their respective charges to optimize their interactions at bound states. The deprotonation energy (DPE) of the acids and the protonation energy (Eprot) of the gaseous analogs of bound intermediates and transition states reflect their respective properties as species present at non-interacting distances. These properties accurately describe the reactivity of acids of a given type, such as polyoxometalates (POM) with a given type of addenda atom but different central atoms and heterosilicates, for different families of reactions. They do not fully capture, however, differences among acid types (e.g., Mo and W POM, heterosilicates, mineral acids) for diverse types of chemical transformations (e.g., elimination, isomerization, dimerization, condensation). The incompleteness of such descriptors reflects their inability to describe how protonated molecular species and conjugate anions restructure their respective charges when present as a binding pair at interacting distances. Such interaction energies represent electrostatic forces that depend on charge distributions in the cations and anions and the ability to reorganize the distributions to maximize the interactions. In the case of deprotonation, the electrostatic and charge reorganization components of DPE for various acids solely reflect the ability of the conjugate anion to accept and distribute the negative charge, a characteristic unique of each type of solid acid and specifically of the composition of its extended conjugate anion framework. The energy required to accept and rearrange the positive charge in bound intermediates and transition states reflects, in turn, their respective ability to recover the ionic and covalent components of DPE, the energy required to detach proton from conjugate anions. The DPE components and the recovery fractions together lead to a modified DPE, which captures only the part of DPE that remains unrecovered by the ion-pair interactions at bound intermediates and transition states, as the unifying descriptor for broad families of acids and reactions. The electrostatic and charge reorganization energies involved in these general descriptors are placed in historical context by assessing their connections to the heuristics of hard–soft acid–base displacements. Further development of these concepts requires benchmarking and extension of electrostatic and reorganization components of energies for a more diverse set of reaction types and acid families and advancement of methods for more efficient calculations of electrostatic interactions. Reactivity descriptors must also account for dispersive interactions between host cavities and guest molecules, requiring a framework analogous to the one described here for ion-pair interactions; these dispersive interactions depend on the fit between their shapes and sizes as well as their ‘‘structural stiffness’’ that determines the ability to modify the shapes of molecules and voids to minimize free energy. Entropy considerations and estimates of their dependence on properties of catalysts and molecules are also required for accurately determining Gibbs free energies that ultimately determine reaction rates.

Deshlahra, Prashant↗
DOE OSTI2025

Elucidation of Marcus Relationships for Hydride Transfer Reactions Involving Transition Metal Hydrides

The rate of hydride transfer from three Ir hydride complexes of the type Cp*Ir( R bpy)H + (Cp* = C 5 Me 5 ; R bpy = 4,4′-R-2,2′-bipyridine, R = OMe, H, CO 2 Me) to six N-methylacridinium ( R Acr + ) acceptors with electronically different substituents in the 2- or 2,7-positions were measured. Using the thermodynamic hydricity of the donors and the hydride affinity of the acceptors the thermodynamic driving forces for hydride transfer were determined. Brønsted plots, which correlate kinetic and thermodynamic hydricity, demonstrate distinct linear free energy relationships for each complex, with different Brønsted α values. Thus, at the same driving force hydride transfer from Cp*Ir( OMe bpy)H + is faster than for Cp*Ir(bpy)H + or Cp*Ir( CO2Me bpy)H + . Experimental and computational analyses are consistent with a concerted hydride transfer mechanism for all Ir complexes. As the thermodynamic driving force increases an earlier transition state is observed and all transition states also include π-stacking interactions between the donor and acceptor, which likely contribute to the different α values. The experimental data fits well to the Marcus model, enabling the determination of reorganization energies (λ) that range from 58 to 69 kcal mol -1 . These are lower than λ values for hydride transfer reactions involving organic donors and acceptors. This work provides a rare example of the correlation of kinetic and thermodynamic hydricity using only experimental data and shows that hydride transfer reactions involving metal hydrides can follow Marcus theory. Furthermore, the findings offer insight into controlling metal-catalyzed hydride transfer reactions, which is valuable for designing improved systems for a range of transformations.

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

Machine Learning Analysis of Direct Dynamics Trajectory Outcomes for Thermal Deazetization of 2,3-Diazabicyclo[2.2.1]hept-2-ene

Experimentally, the thermal gas-phase deazetization of 2,3-diazabicyclo[2.2.1]hept-2-ene (1) results in the loss of N2 and the formation of bicyclo products 3 (exo) and 4 (endo) in a nonstatistical ratio, with preference for the exo product. Here, we report unrestricted M06-2X quasiclassical trajectories initialized from the concerted N2 ejection transition state that were able to replicate the experimental preference to form 3. We found that the 3:4 ratio results from the relative amounts of very fast (ballistic) exotype trajectories versus trajectories that lead to the 1,3-diradical intermediate 2. These quasiclassical trajectories provided a set of transition-state vibrational, velocity, momenta, and geometric features for the machine learning analysis. Additionally, a selection of popular supervised classification algorithms (e.g., random forest) provided poor prediction of trajectory outcomes based on only transition-state vibrational quanta and energy features. However, these machine learning models provided more accurate predictions using atomic velocities and atomic positions, attaining ~70% accuracy using initial conditions and between 85 and 95% accuracy at later reaction time steps. This increased accuracy allowed the feature importance analysis to reveal that, at the later-time analysis, the methylene bridge out-of-plane bending is correlated with trajectory outcomes for the formation of either the exo product or toward the diradical intermediate. Possible reasons for the struggle of machine learning algorithms to classify trajectories based on transition-state features is the heavily overlapping feature values, the finite but very large possible vibrational mode combinations, and the possibility of chaos as trajectories propagate. We examined this chaos by comparing a set of nearly identical trajectories that differed by only a very small scaling of the kinetic energies resulting from the transition-state reaction coordinate.

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

Product specific thermal degradation kinetics of bisphenol F epoxy in inert and oxidative atmospheres using evolved gas analysis–mass spectrometry

Knowledge of the degradation kinetics for polymer materials is important for understanding thermal stability. In this study, evolved gas analysis–mass spectrometry and pyrolysis gas-chromatography-mass spectrometry were evaluated for the potential to deliver additional insight into thermal degradation kinetics of diglycidal ether of bisphenol F (DGEBF) epoxy thermoset under inert and oxidative atmospheres. Degradation products of selected precursor ions were evaluated for their uniqueness to the specific precursor using extracted ion thermographs. Unique mass peaks, solely attributed to a single reaction pathway of a specific product, were determined from extracted ion thermographs and used to determine both activation energy (E a ) and pre-exponential factors for the specific primary reaction pathways. These primary reaction pathways for DGEBF epoxy degradation were then evaluated in the context of transition state theory (TST) and related transition state enthalpies (ΔH ‡ ) and entropies (ΔS ‡ ) of activation to further elucidate the degradation process. It was determined under pyrolysis conditions, as suggested by the E a , the formation of bisphenol F monomer was the rate-limiting step toward the formation of xanthene and phenol. In contrast, under thermo-oxidative conditions, reactions involving oxygen containing species were identified as the rate-limiting step for all observed products based on the large negative ΔS ‡ calculated from TST. This work demonstrates a powerful combination of technique and theory that can provide new insight into the degradation of polymer materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2007

State-Chart Autocoder

A computer program translates Unified Modeling Language (UML) representations of state charts into source code in the C, C++, and Python computing languages. ( State charts signifies graphical descriptions of states and state transitions of a spacecraft or other complex system.) The UML representations constituting the input to this program are generated by using a UML-compliant graphical design program to draw the state charts. The generated source code is consistent with the "quantum programming" approach, which is so named because it involves discrete states and state transitions that have features in common with states and state transitions in quantum mechanics. Quantum programming enables efficient implementation of state charts, suitable for real-time embedded flight software. In addition to source code, the autocoder program generates a graphical-user-interface (GUI) program that, in turn, generates a display of state transitions in response to events triggered by the user. The GUI program is wrapped around, and can be used to exercise the state-chart behavior of, the generated source code. Once the expected state-chart behavior is confirmed, the generated source code can be augmented with a software interface to the rest of the software with which the source code is required to interact.

Clark, Kenneth↗
DOE OSTI2022

Nucleation Rate Theory for Coordination Number: Elucidating Water-Mediated Formation of a Zigzag Na 2 SO 4 Morphology

Predicting and controlling nanostructure formation during nucleation can pave the way to synthesizing novel energy materials via crystallization. However, such control over nucleation and crystallization remains challenging due to an inadequate understanding of critical factors that govern evolving atomistic structures and dynamics. Here, in this paper, we utilize coordination number as a reaction coordinate and rate theory to investigate how sodium sulfate, commonly known as a phase-change energy material, nucleates in a supersaturated aqueous solution. In conjunction with ab initio and force field-based molecular dynamics simulation, the rate theoretical analysis reveals that sodium sulfate from an initially dissolved metastable state transits to a heterogeneous mixture of prenucleated clusters and finally to a large cylindrical zigzag morphology. Measurements of Raman spectra and their ab initio modeling confirm that this nucleated morphology contains a few waters for every sulfate. Rate processes such as solvent exchange and desolvation exhibit high sensitivity to the evolving prenucleation/nucleation structures, providing a means to distinguish between critical nucleation precursors. Desolvation and forming the first-shell interionic coordination structure via monomer-by-monomer addition around sulfates are found to explain the formation of large nuclei. Thus, a detailed understanding of the step-by-step structure formation across scales has been achieved. This can be leveraged to predict nucleation-related structures and dynamics and potentially control the synthesis of novel phase-change materials for energy applications.

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

Structural trends in the dehydrogenation selectivity of palladium alloys

Alloying is well-known to improve the dehydrogenation selectivity of pure metals, but there remains considerable debate about the structural and electronic features of alloy surfaces that give rise to this behavior. To provide molecular-level insights into these effects, a series of Pd intermetallic alloy catalysts with Zn, Ga, In, Fe and Mn promoter elements was synthesized, and the structures were determined using in situ X-ray absorption spectroscopy (XAS) and synchrotron X-ray diffraction (XRD). The alloys all showed propane dehydrogenation turnover rates 5–8 times higher than monometallic Pd and selectivity to propylene of over 90%. Moreover, among the synthesized alloys, Pd 3 M alloy structures were less olefin selective than PdM alloys which were, in turn, almost 100% selective to propylene. This selectivity improvement was interpreted by changes in the DFT-calculated binding energies and activation energies for C–C and C–H bond activation, which are ultimately influenced by perturbation of the most stable adsorption site and changes to the d-band density of states. Furthermore, transition state analysis showed that the C–C bond breaking reactions require 4-fold ensemble sites, which are suggested to be required for non-selective, alkane hydrogenolysis reactions. These sites, which are not present on alloys with PdM structures, could be formed in the Pd 3 M alloy through substitution of one M atom with Pd, and this effect is suggested to be partially responsible for their slightly lower selectivity.

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

Assessment of enzyme active site positioning and tests of catalytic mechanisms through X-ray–derived conformational ensembles

How enzymes achieve their enormous rate enhancements remains a central question in biology, and our understanding to date has impacted drug development, influenced enzyme design, and deepened our appreciation of evolutionary processes. While enzymes position catalytic and reactant groups in active sites, physics requires that atoms undergo constant motion. Numerous proposals have invoked positioning or motions as central for enzyme function, but a scarcity of experimental data has limited our understanding of positioning and motion, their relative importance, and their changes through the enzyme’s reaction cycle. To examine positioning and motions and test catalytic proposals, we collected “room temperature” X-ray crystallography data for Pseudomonas putida ketosteroid isomerase (KSI), and we obtained conformational ensembles for this and a homologous KSI from multiple PDB crystal structures. Here, ensemble analyses indicated limited change through KSI’s reaction cycle. Active site positioning was on the 1- to 1.5-Å scale, and was not exceptional compared to noncatalytic groups. The KSI ensembles provided evidence against catalytic proposals invoking oxyanion hole geometric discrimination between the ground state and transition state or highly precise general base positioning. Instead, increasing or decreasing positioning of KSI’s general base reduced catalysis, suggesting optimized Ångstrom-scale conformational heterogeneity that allows KSI to efficiently catalyze multiple reaction steps. Ensemble analyses of surrounding groups for WT and mutant KSIs provided insights into the forces and interactions that allow and limit active-site motions. Most generally, this ensemble perspective extends traditional structure–function relationships, providing the basis for a new era of “ensemble–function” interrogation of enzymes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Fullerene Negative Ions: Formation and Catalysis

We first explore negative-ion formation in fullerenes C 44 to C 136 through low-energy electron elastic scattering total cross sections calculations using our Regge-pole methodology. Then, the formed negative ions C 44 – to C 136 – are used to investigate the catalysis of water oxidation to peroxide and water synthesis from H 2 and O 2 . The exploited fundamental mechanism underlying negative-ion catalysis involves hydrogen bond strength-weakening/breaking in the transition state. Density Functional Theory transition state calculations found C 60 – optimal for both water and peroxide synthesis, C 100 – increases the energy barrier the most, and C 136 – the most effective catalyst in both water synthesis and oxidation to H 2 O 2 .

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

Magnetic Flux Plays an Important Role during a Black Hole X-Ray Binary Outburst in Radiative Two-temperature General Relativistic Magnetohydrodynamic Simulations

Abstract Black hole (Bh) X-ray binaries cycle through different spectral states of accretion over the course of months to years. Although persistent changes in the Bh mass accretion rate are generally recognized as the most important component of state transitions, it is becoming increasingly evident that magnetic fields play a similarly important role. In this article, we present the first radiative two-temperature general relativistic magnetohydrodynamics simulations in which an accretion disk transitions from a quiescent state at an accretion rate of M ̇ ∼ 10 − 10 M ̇ Edd to a hard-intermediate state at an accretion rate of M ̇ ∼ 10 − 2 M ̇ Edd . This huge parameter space in mass accretion rate is bridged by artificially rescaling the gas density scale of the simulations. We present two jetted BH models with varying degrees of magnetic flux saturation. We demonstrate that in “standard and normal evolution” models, which are unsaturated with magnetic flux, the hot torus collapses into a thin and cold accretion disk when M ̇ ≳ 5 × 10 − 3 M ̇ Edd . On the other hand, in “magnetically arrested disk” models, which are fully saturated with vertical magnetic flux, the plasma remains mostly hot with substructures that condense into cold clumps of gas when M ̇ ≳ 1 × 10 − 2 M ̇ Edd . This suggests that the spectral signatures observed during state transitions are closely tied to the level of magnetic flux saturation.

Astronomy & Astrophysics↗
NASA NTRSDate not supplied

High Reynolds Number Testing of the NATO AVT-298 SWiFT Configuration at the National Transonic Facility

A high Reynolds number wind tunnel test of the NATO AVT-298 Swept Wing Flow Test (SWiFT) configuration was conducted in the National Transonic Facility at the NASA Langley Research Center during the summer of 2023. The SWiFT research model geometry has relevance to both blended/hybrid wing body (BWB/HWB) and unmanned combat aerial vehicle (UCAV) configurations, and the test campaign was the culmination of an international collaboration under the NATO AVT-298 research task group. The main objectives of the test were to investigate Reynolds number scaling effects on low-speed stability & control characteristics and to examine the onset and progression of flow separation on the wings particularly near the wing crank. Force & moment and surface pressure data were acquired at Mach numbers from 0.2 to 0.8, Reynolds numbers from 2.5 to 34 million, angles of attack from -3 to 20 degrees, and sideslip angles from -10 to 10 degrees. Boundary layer transition detection techniques utilizing static pressure taps, unsteady pressure transducers, and sublimating chemicals were used on the model in a free/natural transition state or a forced transition state using trip dots. Pressure sensitive paint was used to obtain a global surface pressure profile on the model and an advanced laser velocimetry technique was used to obtain velocity measurements in the wake downstream of the wing crank. The results from the test showed clear Reynolds number scaling effects on the maximum lift coefficient and pitching moment coefficient at Mach 0.2, while also capturing significant hysteresis effects. The data acquired from the test will ultimately help improve computational aerodynamic analysis and design tools for application to future BWB-type vehicle configurations.

NATO AVT-298↗
NASA NTRSDate not supplied

High Reynolds Number Testing of the NATO AVT-298 SWiFT Configuration at the National Transonic Facility

A high Reynolds number wind tunnel test of the NATO AVT-298 Swept Wing Flow Test (SWiFT) configuration was conducted in the National Transonic Facility at the NASA Langley Research Center during the summer of 2023. The SWiFT research model geometry has relevance to both blended/hybrid wing body (BWB/HWB) and unmanned combat aerial vehicle (UCAV) configurations, and the test campaign was the culmination of an international collaboration under the NATO AVT-298 research task group. The main objectives of the test were to investigate Reynolds number scaling effects on low-speed stability & control characteristics and to examine the onset and progression of flow separation on the wings particularly near the wing crank. Force & moment and surface pressure data were acquired at Mach numbers from 0.2 to 0.8, Reynolds numbers from 2.5 to 34 million, angles of attack from -3 to 20 degrees, and sideslip angles from -10 to 10 degrees. Boundary layer transition detection techniques utilizing static pressure taps, unsteady pressure transducers, and sublimating chemicals were used on the model in a free/natural transition state or a forced transition state using trip dots. Pressure sensitive paint was used to obtain a global surface pressure profile on the model and an advanced laser velocimetry technique was used to obtain velocity measurements in the wake downstream of the wing crank. The results from the test showed clear Reynolds number scaling effects on the maximum lift coefficient and pitching moment coefficient at Mach 0.2, while also capturing significant hysteresis effects. The data acquired from the test will ultimately help improve computational aerodynamic analysis and design tools for application to future BWB-type vehicle configurations.

BWB/HWB↗
DOE OSTI2023

The multichannel i -propyl + O2 reaction system: A model of secondary alkyl radical oxidation

The i-propyl + O2 reaction mechanism has been investigated by definitive quantum chemical methods to establish this system as a benchmark for the combustion of secondary alkyl radicals. Focal point analyses extrapolating to the ab initio limit were performed based on explicit computations with electron correlation treatments through coupled cluster single, double, triple, and quadruple excitations and basis sets up to cc-pV5Z. The rigorous coupled cluster single, double, and triple excitations/cc-pVTZ level of theory was used to fully optimize all reaction species and transition states, thus, removing some substantial flaws in reference geometries existing in the literature. The vital i-propylperoxy radical (MIN1) and its concerted elimination transition state (TS1) were found 34.8 and 4.4 kcal mol−1 below the reactants, respectively. Two β-hydrogen transfer transition states (TS2, TS2′) lie above the reactants by (1.4, 2.5) kcal mol−1 and display large Born–Oppenheimer diagonal corrections indicative of nearby surface crossings. An α-hydrogen transfer transition state (TS5) is discovered 5.7 kcal mol−1 above the reactants that bifurcates into equivalent α-peroxy radical hanging wells (MIN3) prior to a highly exothermic dissociation into acetone + OH. The reverse TS5 → MIN1 intrinsic reaction path also displays fascinating features, including another bifurcation and a conical intersection of potential energy surfaces. An exhaustive conformational search of two hydroperoxypropyl (QOOH) intermediates (MIN2 and MIN3) of the i-propyl + O2 system located nine rotamers within 0.9 kcal mol−1 of the corresponding lowest-energy minima.

Chemistry↗
NASA NTRS2005

X-Ray and Radio Studies of Black Hole X-Ray Transients During Outburst Decay

Black hole (BH) and black hole candidate (BHC) transients are X-ray binary systems that typically undergo bright outbursts that last a couple months with recurrence times of years to decades. For this ADP project, we are studying BH/BHC systems during the decaying phases of their outbursts using the Rossi X-ray Taming Explorer (RXTE), the Chandra X-ray Observatory, and multi-wavelength facilities. These systems usually undergo state transitions as they decay, and our observations are designed to catch the state transitions. The specific goals of this proposal include: 1. To determine the evolution of the characteristic frequencies present in the power spectrum (such as quasi-periodic oscillations, QPOs) during state transitions in order to place constraints on the accretion geometry; 2. To contemporaneously measure X-ray spectral and timing properties along with flux measurements in the radio band to determine the relationship between the accretion disk and radio jets; 3. To extend our studies of X-ray properties of BHCs to very low accretion rates using RXTE and Chandra. The work performed under this proposal has been highly successful, allowing the PI to lead, direct, or assist in the preparation of 7 related publications in refereed journals and 6 other conference presentations or reports. These items are listed below, and the abstracts for the refereed publications have also been included. Especially notable results include our detailed measurements of the characteristic frequencies and spectral parameters of BH/BHCs after the transition to the hard state (see All A3, and A5) and at low flux levels (see A4). Our measurements provide one of the strongest lines of evidence to date that the inner edge of the optically thick accretion disk gradually recedes from the black hole at low flux levels. In addition, we have succeeded in obtaining excellent multi-wavelength coverage of a BH system as its compact jet turned on (see Al). Our results show, somewhat unexpectedly, that the radio jet does not turn on until the hard X-ray emission is well past its peak hard state level, strongly constraining theoretical models for hard X-ray production and the spectrum emitted by the jet. Finally, the X-ray/radio results in A2 led us to propose a general picture about the relationship between jet production and X-ray spectral states .

Tomsick, John A.↗
DOE OSTI2022

A Combined DFT, Energy Decomposition, and Data Analysis Approach to Investigate the Relationship Between Noncovalent Interactions and Selectivity in a Flexible DABCOnium/Chiral Anion Catalyst System

Developing strategies to study reactivity and selectivity in flexible catalyst systems has become an important topic of research. Herein, we report a combined experimental and computational study aimed at understanding the mechanistic role of an achiral DABCOnium cofactor in a regio- and enantiodivergent bromocyclization reaction. In this study, it was found that electron-deficient aryl substituents enable rigidified transition states via an anion–π interaction with the catalyst, which drives the selectivity of the reaction. In contrast, electron-rich aryl groups on the DABCOnium result in significantly more flexible transition states, where interactions between the catalyst and substrate are more important. An analysis of not only the lowest-energy transition state structures but also an ensemble of low-energy transition state conformers via energy decomposition analysis and machine learning was crucial to revealing the dominant noncovalent interactions responsible for observed changes in selectivity in this flexible system.

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

The reaction NH 2 + CH 2 O: Kinetic measurements and computational studies

The reaction of amidogen with formaldehyde is relevant to astrochemistry and couples nitrogen and carbon chemistry in flames. The first measurements of the temperature dependence of the rate constant k were made, over 409–643 K, and are summarized as k = 7.1 × 10 −12 exp(−18.2 kJ mol −1 /RT) cm 3 molecule −1 s −1 with a confidence interval of ±20%. These results address a controversy over the presence of a significant barrier in the addition/elimination pathway leading to the production of formamide between theoretical models that include the zero-point energy of all modes at the transition state and a model that excludes transitional modes. The latter would lead to a negligible barrier, which is inconsistent with the experiments. Contrary to earlier claims, CCSDT(Q)-corrected energies combined with transition state theory yield quantitative accord with the measurements.

Astrochemistry↗
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↗