Search NASA⌕ Search

SEARCH · Search NASA

Results for “rate constants”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Pressure-dependent kinetics of cyclopentene consumption

Recent experiments have uncovered deficiencies in the ability of current kinetic models for cyclopentane combustion to predict species profiles as a function of O2 concentration. One likely source of these discrepancies flows from the limited availability of relevant kinetic data for key functionalized intermediates, requiring the use of distant analogies to estimate their rates. This work fills this gap for cyclopentene, the most prominent C5 intermediate in low-temperature cyclopentane combustion. Ab initio transition-state-theory based master equation methods are used to calculate pressure-dependent rate constants across four potential energy surfaces (C5H8 + OH, C5H8 + HO2, C5H7, and C5H7 + O2). Potential energy surfaces are characterized at the CCSD(T)-F12/cc-pVDZ-F12//revDSD-PBEP86-D3BJ/def2-TZVPP level of theory. Pressure-dependent rate constants and branching fractions are determined via 1D master equation calculations. The results are used to analyze the competition for radicals between cyclopentene and cyclopentane in the intermediate temperature region, and to identify the dominant pathways flowing from these initiation channels. On the OH initiation surface, we find good agreement with recently published experimental rate constants at high temperature, and we elucidate the complex temperature dependence below 800 K, where 𝜋𝜋 addition becomes prominent. On the HO2 initiation surface, we find an interesting alternative route to bicyclic ether formation via well-skipping over the HO2 adduct well. On the cyclopentenyl radical surface, the major product is cyclopentadiene plus H. On the cyclopentenyl + O2 surface, we find that resonance stabilization largely prevents the allylic radical from undergoing oxidation, whereas the alkylic radical has a somewhat elevated propensity for oxidation pathways relative to cyclopentyl.

ab initio calculations↗
DOE OSTI2022

An SECM-Based Spot Analysis for Redoxmer-Electrode Kinetics: Identifying Redox Asymmetries on Model Graphitic Carbon Interfaces

The fundamental process in non-aqueous redox flow battery (NRFB) operation revolves around electron transfer (ET) between a current collector electrode and redox-active organic molecules (redoxmers) in solution. Here, we present an approach utilizing scanning electrochemical microscopy (SECM) to evaluate interfacial ET kinetics between redoxmers and various electrode materials of interest at desired locations. This spot-analysis method relies on the measurement of heterogeneous electron transfer rate constants (k f or k b ) as a function of applied potential (E-E 0 '). As demonstrated by COMSOL simulations, this method enables the quantification of Butler-Volmer kinetic parameters, the standard heterogeneous rate constant, k 0 , and the transfer coefficient, α. Our method enabled the identification of inherent asymmetries in the ET kinetics arising during the reduction of ferrocene-based redoxmers, compared to their oxidation which displayed faster rate constants. Similar behavior was observed on a wide variety of carbon electrodes such as multi-layer graphene, highly ordered pyrolytic graphite, glassy carbon, and chemical vapor deposition-grown graphite films. However, aqueous systems and Pt do not exhibit such kinetic effects. Our analysis suggests that differential adsorption of the redoxmers is insufficient to account for our observations. Displaying a greater versatility than conventional electroanalytical methods, we demonstrate the operation of our spot analysis at concentrations up to 100 mM of redoxmer over graphite films. Looking forward, our method can be used to assess non-idealities in a variety of redoxmer/electrode/solvent systems with quantitative evaluation of kinetics for applications in redox-flow battery research.

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

Implementation of new mixture rules has a substantial impact on combustion predictions for H 2 and NH 3

Complex-forming reactions comprise a substantial fraction of all important combustion reactions and are central to combustion behavior. Despite being often called “pressure-dependent” reactions, their rate constants depend on not only the pressure but also the composition. While modern combustion codes allow arbitrarily high accuracy in treating pressure dependence, recent work has consistently demonstrated dramatic failures of essentially all available treatments of mixture dependence. In situations where mixture dependence is treated at all, it is inevitably treated through specification of pressure-dependent rate constants for a set of pure bath gases, which are then combined to estimate the rate constant in a mixture via a “mixture rule.” While there had been a generally unquestioning confidence in these mixture rules, they had, in reality, been scarcely tested until the last decade, when comparisons against master equation calculations revealed order-of-magnitude errors for important pressure-dependent reactions. New mixture rules, based on the reduced pressure, have recently been proposed and shown to reproduce master equation calculations for broad classes of complex-forming reactions very accurately. Here, in this work, we present an implementation of one such new mixture rule (“LMR-R”) in Cantera and then use it to enable simulations that use new high-accuracy ab initio data for individual bath gases (for the first time, since codes previously could not accommodate the complex bath gas dependence). Demonstrations focus on combustion of H 2 and NH 3 , where (1) high-accuracy ab initio data are available and (2) the impact is expected to be large due to the high fractions of efficient colliders (e.g., H 2 O and NH 3 ) in the burned and unburned gases. Indeed, we find the impact of this treatment to be substantial and may explain previous modeling difficulties for these important carbon-free fuels, particularly for NH 3 , whose extraordinarily high third-body efficiency (~20) is often omitted from kinetic models.

Ammonia↗
DOE OSTI2021

Role of a 3D Structure in Energy Transfer in Mixed-Ligand Metal–Organic Frameworks

We present a detailed investigation of the photophysical properties of mixed-ligand pyrene- and porphyrin-based metal–organic frameworks (MOFs) as a function of their 3D structure. Solvothermal reactions between metal salts (InCl 3 , Zr(acac) 4 , and ZrCl 4 ) and suitable ratios of 1,3,6,8-tetrakis(p-benzoic acid)pyrene (TBAPy) and meso-tetrakis- (4-carboxyphenyl)porphyrin (TCPP) were performed to prepare a series of mixed-ligand ROD-7, NU-901, and NU-1000 MOFs. Time-resolved and steady-state fluorescence measurements were conducted on the mixed-ligand MOFs to study their photophysics. Based on the results, we concluded that upon excitation of TBAPy linkers in the MOFs, singlet excitation energy migrates across TBAPy linkers until it finds a TCPP unit. TCPP acts as an energy trap and quenches the excitation. The efficiency of TBAPy-to-TBAPy energy transfer was found to be sensitive to the structural parameters of MOFs. Analysis of steady-state and time-resolved fluorescence data revealed that excitation energy transfer (EET) is most efficient in ROD-7, followed by NU-901 and NU-1000. We propose that topology that invokes the shorter interchromophoric distances between TBAPy linkers in ROD-7 is responsible for its higher EET efficiency. The distance dependence of the EET rate constant (kEET) was investigated to gain insight into the mechanistic aspects of energy transfer in MOFs. This study revealed that (a) energy transfer in MOFs deviates from the classical Förster model and (b) the geometrical arrangement of linkers influences the mechanism of EET in MOFs. A theoretical investigation was also performed to determine energy-transfer rate constants along different directions and assess the directionality of energy transfer in these MOFs. The magnitude of rate constants indicated that energy transfer in ROD-7 should be highly anisotropic along the stacking direction. Furthermore, these findings suggest that ROD-7 is a promising candidate to play the role of the light-harvesting and energy-transfer component in solar energy conversion devices, where directional energy transfer is required.

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 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 OSTI2024

Homomolecular Triplet–Triplet Annihilation in Metalloporphyrin Photosensitizers

Metalloporphyrins are ubiquitous in their applications as triplet photosensitizers, particularly for promoting sensitized photochemical upconversion processes. In this study, bimolecular excited state triplet–triplet quenching kinetics, termed homomolecular triplet–triplet annihilation (HTTA), exhibited by the traditional triplet photosensitizers–zinc(II) tetraphenylporphyrin (ZnTPP), palladium(II) octaethylporphyrin (PdOEP), platinum(II) octaethylporphyrin (PtOEP), and platinum(II) tetraphenyltetrabenzoporphyrin (PtTPBP)–were revealed using conventional transient absorption spectroscopy. Nickel(II) tetraphenylporphyrin was used as a control sample as it is known to be rapidly quenched intramolecularly through ligand-field state deactivation and, therefore, cannot result in triplet–triplet annihilation (TTA). The single wavelength transients associated with the metalloporphyrin triplet excited state decay–measured as a function of incident laser pulse energy in toluene–were well modeled using parallel first- and second-order kinetics, consistent with HTTA being operable. The combined transient kinetic data enabled the determination of the first-order rate constants (k T ) for excited triplet decay in ZnTPP (4.0 × 10 3 s –1 ), PdOEP (3.6 × 10 3 s –1 ), PtOEP (1.2 × 10 4 s –1 ), and PtTPBP (2.1 × 10 4 s –1 ) as well as the second-order rate constant (k TT ) for HTTA in ZnTPP (5.5 × 10 9 M –1 s –1 ), PdOEP (1.1 × 10 10 M –1 s –1 ), PtOEP (7.1 × 10 9 M –1 s –1 ), and PtTPBP (1.6 × 10 10 M –1 s –1 ). In most instances, triplet excited state extinction coefficients are either reported for the first time or have been revised using ultrafast transient absorption spectroscopy and singlet depletion: ZnTPP (78,000 M –1 cm –1 ) at 470 nm, PdOEP (67,000 M –1 cm –1 ) at 430 nm, PtOEP (51,000 M –1 cm –1 ) at 418 nm, and PtTPBP (100,000 M –1 cm –1 ) at 460 nm. Furthermore, the combined experimental results establish competitive time scales for homo- and heteromolecular TTA rate constants, implying the significance of considering HTTA processes in future research endeavors harnessing TTA photochemistry using common metalloporphyrin photosensitizers.

14 SOLAR ENERGY↗
DOE OSTI2021

Transport Modeling of Kinetically Limited Microscale Extraction Systems: Droplet and Supported Liquid Membrane Separations

Microfluidic droplet extraction devices can be used to efficiently measure solute phase transfer kinetics in liquid–liquid extraction systems. In this work, a droplet extraction apparatus was constructed using off-the-shelf and custom 3D printed parts for less than 110 USD. The interfacial kinetics of the transfer of Pu(IV) between an aqueous nitric acid phase and a 5 v/v% tridodecylamine (TDA) in toluene phase were characterized by determining the extent of Pu(IV) transfer between phases at different flow rates. Forward interfacial rate constants calculated from this data using an analytical model are 5.6 × 10 –3 mm s –1 under extraction (aqueous-to-organic solute transfer) conditions and 3.3 × 10 –3 mm s –1 under stripping (organic-to-aqueous solute transfer) conditions. The analytical model is valid when the rate of solute transfer across the interface between the two liquid phases has pseudo-first-order reversible reaction kinetics and is slow relative to the rate of solute transport to the interface. The change in the outlet concentration with residence time determined using the analytical model was compared with the results of a numerical model accounting for advection/diffusion to the interface in addition to kinetically limited interfacial solute transfer. These calculations revealed that for systems with forward interfacial rate constants greater than approximately 10 –2 mm s –1 , the analytical model must be supplemented with a numerical transport model accounting for an additional advection/diffusion mass transfer resistance. This threshold corresponds to a Damköhler number of one. The measured interfacial kinetic rate constants for Pu(IV)/TDA were used in a supported liquid membrane (SLM) numerical model, which showed that the Pu(IV) phase transfer kinetics are fast enough for use in a practical separation.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

A theoretical kinetic study of ĊH 3 + ṄH 2 : From electronic structure to NH 3 /CH 4 combustion modelling implications

Carbon–nitrogen interaction reactions play an important role in governing the reactivity of ammonia blended fuels. However, there remains uncertainties regarding their detailed reaction pathways and rate constants, hampering the development of high-fidelity chemical kinetic models. In this study, the kinetics of ĊH 3 + ṄH 2 , a key C–N interaction reaction in ammonia/methane blend combustion have been investigated. The potential energy surface has been explored using the high-level ANL0F method, yielding highly accurate stationary point energies that agree with ATcT values within 0.1 kcal mol –1 . Variable reaction coordinate transition state theory is used to treat the barrierless association and decomposition reaction channels, based on directly sampled radical-radical interaction energies at the CASPT2-F12(2e,2o)/cc-pVTZ-F12 level of theory. The minimum transitional mode numbers of states obtained are then coupled with the RRKM/master equation to calculate temperature- and pressure-dependent rate constants. Our a priori calculations capture available experimental measurements from the literature very well. The calculated rate constants have been incorporated into an NH 3 /CH 4 chemical kinetic model currently under development at the University of Galway. The effect of the updated kinetic data for ĊH 3 + ṄH 2 on model predicted NH 3 /CH 4 fuel reactivity is elucidated.

ab initio↗
DOE OSTI2024

Probabilistic and maximum entropy modeling of chemical reaction systems: Characteristics and comparisons to mass action kinetic models

We demonstrate and characterize a first-principles approach to modeling the mass action dynamics of metabolism. Starting from a basic definition of entropy expressed as a multinomial probability density using Boltzmann probabilities with standard chemical potentials, we derive and compare the free energy dissipation and the entropy production rates. We express the relation between entropy production and the chemical master equation for modeling metabolism, which unifies chemical kinetics and chemical thermodynamics. Because prediction uncertainty with respect to parameter variability is frequently a concern with mass action models utilizing rate constants, we compare and contrast the maximum entropy model, which has its own set of rate parameters, to a population of standard mass action models in which the rate constants are randomly chosen. We show that a maximum entropy model is characterized by a high probability of free energy dissipation rate and likewise entropy production rate, relative to other models. We then characterize the variability of the maximum entropy model predictions with respect to uncertainties in parameters (standard free energies of formation) and with respect to ionic strengths typically found in a cell.

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

Predicting third-body collision efficiencies for water and other polyatomic baths

Low-pressure-limit microcanonical (collisional activation) and thermal rate constants are predicted using a combination of automated ab initio potential energy surface construction, classical trajectories, transition state theory, and a detailed energy- and angular-momentum-resolved collision kernel. Several systems are considered, including CH 4 (+M) and HO 2 (+M), with an emphasis on systems where experimental information is available for comparison. The a priori approach involves no adjustable parameters, and we show that the predicted thermal rate constants are in excellent agreement with experiments, with average deviations of less than 25%. Notably, the a priori approach is shown to perform equally well for atomic, diatomic, and polyatomic baths, including M = H 2 O, CO 2 , and “fuel” baths like M = CH 4 and NH 3 . Additionally, the utility of microcanonical rate constants for interpreting trends and inferring mechanistic details in the thermal kinetics is demonstrated.

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

Effect of Strength of Soluble and Solid Acid Catalysts on Selectivity in Aldol Reactions

Aldol reactions can follow the classical condensation path and produce an α,β-unsaturated carbonyl compound and water, however in some cases, aldol chemistry may be steered towards a fission that gives an olefin and a carboxylic acid. This investigation examines how the strength of Brønsted acid sites—both in homogeneous and heterogeneous catalysis—controls these two pathways and their kinetics. The cross-aldol reaction between benzaldehyde and 3-pentanone served as test case. Batch reactions, conducted in toluene as solvent at a temperature of 140 °C under autogenous pressure, were analyzed by GC and in situ ATR-FTIR spectroscopy to determine product distributions and rate constants for condensation and fission pathways. A series of soluble acids including a family of sulfonic acids mostly favored the condensation pathway, with formic acid as the weakest in the series by its pKa being inactive for aldol chemistry. Significant amounts of fission products were rare except for the known selectivity of phosphoric acid, and the circumstances were difficult to delineate. For the sulfonic acid family, the logarithm of the first-order condensation rate constants scaled only roughly with pKa (water) values, whereas a good correlation was obtained with calculated deprotonation Gibbs energies in toluene. A series of H-forms of isomorphously substituted beta zeolites, HESiBEA with E = Al, Ga, Fe, or B, favored the fission pathway. Site density and site strength were characterized by calorimetric measurements of the heats of adsorption of isopropylamine, which decreased in the order Al > Ga, Fe > B. The logarithm of the site-normalized first-order fission rate constants scaled roughly with the heats of adsorption and correlated excellently with reported deprotonation energies. In conclusion, acid strength mainly affects activity and can be seen as a pre-requisite for either aldol condensation or fission chemistry, whereas additional, yet to be full clarified catalyst properties and reaction conditions are required to steer aldol chemistry towards fission selectivity.

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

Quantifying the Decomposition Kinetics of Linear C 2 –C 4 Perfluoroalkanes

The pyrolysis of C 2 –C 4 linear perfluoroalkanes was investigated using laser schlieren densitometry (LS) in a diaphragmless shock tube. Mixtures of 1, 2, and 4% perfluoroethane (PFE), perfluoropropane (PFP), and perfluorobutane (PFB) in excess krypton were shock-heated to 1400–2500 K at pressures of 60, 120, and 240 Torr. The initiation step for all perfluoroalkanes involved carbon–carbon bond fission. The rate constants for unimolecular decomposition were computed using microcanonical rate theory, and these rate constants were combined into a detailed chemical kinetic mechanism. For PFE and PFP, there was one dissociation channel, while PFB dissociated via two competing bond-fission reactions. The perfluoroalkyl radicals subsequently recombined or underwent further unimolecular decomposition. In conclusion, the kinetic model accurately simulated the density gradient profiles of all perfluoroalkanes.

Guzman, Eduardo H. [Brown Univ., Providence, RI (U↗
DOE OSTI2020

Photophysics of graphene quantum dot assemblies with axially coordinated cobaloxime catalysts

We report a study of chromophore-catalyst assemblies composed of light harvesting hexabenzocoronene (HBC) chromophores axially coordinated to two cobaloxime complexes. The chromophore-catalyst assemblies were prepared using bottom-up synthetic methodology and characterized using solid-state NMR, IR, and x-ray absorption spectroscopy. Detailed steady-state and time-resolved laser spectroscopy was utilized to identify the photophysical properties of the assemblies, coupled with time-dependent DFT calculations to characterize the relevant excited states. The HBC chromophores tend to assemble into aggregates that exhibit high exciton diffusion length (D = 18.5 molecule 2 /ps), indicating that over 50 chromophores can be sampled within their excited state lifetime. We find that the axial coordination of cobaloximes leads to a significant reduction in the excited state lifetime of the HBC moiety, and this finding was discussed in terms of possible electron and energy transfer pathways. By comparing the experimental quenching rate constant (1.0 × 10 9 s -1 ) with the rate constant estimates for Marcus electron transfer (5.7 × 10 8 s -1 ) and Förster/Dexter energy transfers (8.1 × 10 6 s -1 and 1.0 × 10 10 s -1 ), we conclude that both Dexter energy and Marcus electron transfer process are possible deactivation pathways in CoQD-A. No charge transfer or energy transfer intermediate was detected in transient absorption spectroscopy, indicating fast, subpicosecond return to the ground state. Finally, these results provide important insights into the factors that control the photophysical properties of photocatalytic chromophore-catalyst assemblies.

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

Selective chemical looping combustion of acetylene in ethylene-rich streams

Here, the requirement for C 2 H 2 concentrations below 2 parts per million (ppm) in gas streams for C 2 H 4 polymerization necessitates its semihydrogenation to C 2 H 4 . Here, we demonstrate selective chemical looping combustion of C 2 H 2 in C 2 H 4 -rich streams by Bi 2 O 3 as an alternative catalytic pathway to reduce C 2 H 2 concentration below 2 ppm. Bi 2 O 3 combusts C 2 H 2 with a first-order rate constant that is 3000 times greater than the rate constant for C 2 H 4 combustion. In successive redox cycles, the lattice O of Bi 2 O 3 can be fully replenished without discernible changes in local Bi coordination or C 2 H 2 combustion selectivity. Heterolytic activation of C–H bonds across Bi–O sites and the higher acidity of C 2 H 2 results in lower barriers for C 2 H 2 activation than C 2 H 4 , enabling selective catalytic hydrocarbon combustion leveraging differences in molecular deprotonation energies.

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

Stoichiometrically-informed symbolic regression for extracting chemical reaction mechanisms from data

A data-driven computational method is introduced to extract chemical reaction mechanisms from time series chemical concentration data. It is realized through the use of dynamic symbolic regression in which a sparse analytical form for a dynamical system is discoverable from the underlying data. We specifically develop the stoichiometrically-informed symbolic regression (SISR) method to address a standing challenge in complex chemical reaction networks: given a time-series dataset of concentrations of several components, what is the mechanism and the associated rate constants? SISR finds the optimal mechanism, kinetic equations and rate constants by combining differential optimization with a genetic optimization approach that searches a symbolic space of possible reaction mechanisms. Use of SISR in several paradigmatic examples spanning linear and nonlinear reaction schemes results in excellent agreement between true and predicted mechanisms, including when the method is applied to noisy data. The advantages of a stoichiometrically-informed approach such as SISR to address reaction discovery is illustrated through comparison with the use of generic state-of-the-art data-driven approaches.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Association of Cl with C 2 H 2 by unified variable-reaction-coordinate and reaction-path variational transition-state theory

Significance A key issue in chemical kinetics is advancing the theoretical framework to handle reactions beyond the domain of textbook transition-state theory by including––for example––anharmonicity, barrierless transition states, transition states in series, and the effect of conformational flexibility on equilibrium constants. Other key issues are validating affordable electronic structure methods for direct dynamics and the direct calculation of high-pressure limiting rate constants, which are often obtainable experimentally only by extrapolation. The present article addresses all these issues for a prototype radical–molecule association reaction and thereby demonstrates how to combine improved theoretical methods to provide rate constants in cases where experimental data are uncertain or missing.

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

From electronic structure to model application of key reactions for gasoline/alcohol combustion: Hydrogen-atom abstractions by $CH_3\dot{O}$ radical

H-atom abstraction by methoxy radical ($CH_3\dot{O}$) plays an important role in capturing the kinetics of reactions between gasoline components and alcohols. This study focuses on determining the reaction rates and thermodynamic properties of methoxy radical reactions with five gasoline fuel components: n-heptane, iso-octane, 1-hexene, cyclopentane and toluene. Electronic structure calculations were per-formed for all the stationary points with M06-2X/6 -311 ++ g(d,p) method. G3 composite method with atomization method is used for determining Δ f H 0 of all the closed shell and radical species, using which the necessary thermodynamic data of all the species was determined. Coupled cluster theory QCISD(T)/cc-pVXZ (where X = D and T) and Møller-Plesset perturbation theory MP2/cc-pVXZ (where X = D, T and Q) were used to calculate single point energies. Subsequently, rate constants for all hydrogen atom ab-straction channels have been performed using conventional transition state theory with unsymmetric tunneling corrections. A systematic comparison of rates for abstraction from different sites within the same species and same site from different species is done in order to get insights into this reaction class. Here, the computed thermodynamic properties and rate constants were incorporated into a recent gasoline mechanism to investigate the impact of the calculations performed in this work. A shift in predicted NTC (negative temperature coefficient) behavior and significant reduction in model reactivity is observed upon incorporating the rates calculated herein.

33 ADVANCED PROPULSION SYSTEMS↗