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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Kinetic Isotope Effect in the Radical Step-Growth Termination of On-Surface Synthesized Graphene Nanoribbons

Bottom-up on-surface synthesis has emerged as a versatile tool to access and finely tune the electronic structure of nanographenes. The controlled generation of reactive intermediates catalyzed and stabilized by a supporting substrate has enabled the design, assembly, and characterization of a wide range of exotic tailor-made quantum materials. Even under these tightly controlled conditions, the growth of extended structures remains limited by termination processes and undesired side reactions. Here, we identify an H atom transfer as one principal contributor to the radical step growth termination that limits the on-surface growth of N = 7 armchair graphene nanoribbons (7-AGNRs) on Au(111) surfaces. Analysis of 7-AGNR lengths grown from protiated and deuterated molecular precursors reveals a primary kinetic isotope effect of KIE ∼ 1.4. First-principles density functional theory calculations suggest that a concerted H atom transfer mechanism that involves the breaking of a C–H/D bond in the transition state is associated with radical chain termination.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Synergic Effects of Surface Chemistry and Applied Potentials on the Kinetics of the Electrocatalytic Oxygen Evolution Reaction in IrO 2

Improving the efficacy of the oxygen evolution reaction (OER) through water oxidation is critical for advancing photoelectrochemical water splitting. Among the catalysts, IrO 2 exhibits a high OER catalytic activity and stability under acidic conditions. The OER mechanism in this system has been a topic of intense research; however, many mechanistic understandings are lacking. In this work, we carried out first-principles calculations with an implicit solvation model at a constant potential to study the OER behavior on the IrO 2 (110) surface. We find that the surface hydrogen coverage has a significant effect on OER kinetics and transition states. We then develop a microkinetics model that accounts for the continuous evolution of both OER activation energy and hydrogen coverage as a function of an applied potential. We show that this inclusion leads to significant improvement in the simulated Tafel plot compared to available experiments. In conclusion, our results point to a complex interplay between surface chemistry and the applied potential on OER kinetics.

08 HYDROGEN↗
DOE OSTI2021

Reactive Species and Reaction Pathways for the Oxidative Cleavage of 4-Octene and Oleic Acid with H 2 O 2 over Tungsten Oxide Catalysts

The oxidative cleavage of carbon-carbon double bonds (C=C) in alkenes and fatty acids produces aldehydes and acids valued as chemical intermediates. Solid tungsten oxide catalysts are low cost, non-toxic, and selective for the oxidative cleavage of C=C bonds with hydrogen peroxide (H 2 O 2 ) and are, therefore, a promising option for continuous processes. Despite the relevance of these materials, the elementary steps involved and their sensitivity to the form of W sites present on surfaces have not been described. In this work, we combine in situ spectroscopy and rate measurements to identify significant steps in the reaction, the reactive species present on catalysts, and examine differences between the kinetics of this reaction on isolated W atoms grafted to alumina and on those exposed on crystalline WO 3 nanoparticles. Raman spectroscopy shows that W-peroxo complexes (W-(η 2 -O 2 ) formed from H 2 O 2 ) react in a kinetically relevant step to produce epoxides, which undergo hydrolysis at protic surface sites. Subsequently, the CH 3 CN solvent deprotonates diols to form alpha hydroxy ketones that react to form aldehydes and water following nucleophilic attack of H 2 O 2 . Turnover rates for oxidative cleavage, determined by in situ site titrations, on WO x -Al 2 O 3 are 75% greater than those on WO 3 at standard conditions. These differences reflect activation enthalpies (ΔH ‡ ) for oxidative cleavage of 4-octene that are much lower for the isolated WO x sites (36±3 and 60±6 kJ∙mol -1 for WO x -Al 2 O 3 and WO 3 , respectively) and correlate strongly with the difference between enthalpies of adsorption for epoxyoctane (ΔH ads,epox ), which resembles the transition state for epoxidation. The WO x -Al 2 O 3 catalysts mediate oxidative cleavage of oleic acid with H 2 O 2 following a mechanism comparable to that for 4-octene. The WO 3 materials, however, form only the epoxide and do not cleave the C-C bond or produce aldehydes and acids. These differences reflect the distinct site requirements for these reaction pathways and indicate that acid sites required for diol formation are strongly inhibited by oleic acids and epoxides on WO 3 whereas the Al 2 O 3 support provides sites competent for this reaction and increase the yield of the oxidative cleavage products.

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

Optimizing the Solvent Reorganization Free Energy by Metal Substitution for Nanocage Catalysis

The supramolecular capsule Ga 4 L 6 12– has been found to catalyze a number of important chemical reactions, but the close proximity of a poorly organized external solvent environment can serve as an impediment to the reaction chemistry. Here, using ab initio molecular dynamic calculations and electric field analyses, we find that metal substitution of indium for gallium lowers the total and electrostatic activation barriers by ~3–4 kcal/mol in water solvent. Using an energy decomposition analysis of the interaction of water solvent with the metal vertices, we find that Pauli repulsion between the metal and water decreases, allowing the water coordination with the metal to increase, upon substitution of In for Ga. We therefore determine that the stabilization of the transition state is due to a better arrangement of water molecules around the metal vertices of In 4 L 6 12– and provides a proof of concept of how to redesign the nanocage scaffold in order to reduce the solvent reorganization energy.

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

Modeling the Effect of Surface Platinum–Tin Alloys on Propane Dehydrogenation on Platinum–Tin Catalysts

Uncertainty analysis, reported experimental literature data, and density functional theory were synthesized to model the effect of surface tin coverage on platinum-based catalysts for nonoxidative propane dehydrogenation to propylene. Here, this study tests four different platinum–tin skin surface models as potential catalytic sites, Pt 3 Sn/Pt(100), PtSn/Pt(100), Pt 3 Sn/Pt(111), and Pt 2 Sn/Pt(211), and compares them to the corresponding pure Pt surface sites using an uncertainty analysis methodology that uses BEEF-vdW with its ensembles (BMwE) to generate the uncertainty for the energies of the intermediates and transition states. One experimental data set with two experimental observations, selectivity to propylene and turnover frequency of propylene, was used as a calibration data set to evaluate the impact of the experimental data on informing the models. This study finds that the prior model for Pt 3 Sn/Pt(100) is the most active and Pt 2 Sn/Pt(211) is the most selective toward propylene. Active sites on the (100) facet have the highest probability of being responsible for C 1 and C 2 product formations (C–C bond cleavage). Increasing the Sn coverage on the (100) surface facet to a PtSn/Pt(100) active site leads to a significantly reduced rate and might explain the experimentally observed higher selectivity of Sn-doped catalysts relative to pure Pt catalysts. Next, this study finds that for all surfaces, except PtSn/Pt(100), the rate-controlling steps are the initial dehydrogenation steps alongside some partially rate-controlling second dehydrogenation steps. For PtSn/Pt(100), only the initial terminal dehydrogenation step to CH 3 CH 2 CH 2 * and second dehydrogenation steps are rate-controlling. Next, the calibrated models for all surfaces were found to be selective toward propylene production and model the reported turnover frequency successfully. Nevertheless, Pt 2 Sn/Pt(211) emerges as the active site with some (minor) evidence as the main active site based on Jeffreys’ scale interpretation of Bayes factors. This observation agrees with prior studies that also found step sites to be most likely the most relevant active sites for pure Pt catalysts. Overall, the results indicate that tin, in addition to affecting the binding strength of the adsorbed species, prevents deeper dehydrogenation (reducing coking) and cracking reactions through increasing activation barriers for unwanted side reactions.

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

Formation of (Rh–Fe)–FeO x Complex Sites Enables Methanol Synthesis from CO 2

Here, we addressed the challenges of designing catalysts for selective CO 2 hydrogenation by incorporating oxide Fe species onto Rh nanoparticles. Nanoscopic FeO x domains created a “reverse catalyst” structure (i.e., a metal oxide supported on a metal) that increased the density of interfacial sites compared to traditional supported catalysts. The contact between the metal nanoparticle and the oxide overlayer induced the formation of a surface Rh-Fe alloy that stabilize methoxy groups while suppressing hydrogenolysis to methane. Sites at FeO x -metal interfaces interact with CO 2 sevenfold stronger than sites on metal surfaces, show larger energy barriers to cleave the C-O bonds, and offer a barrierless pathway for hydrogenation of methoxy species to methanol. Consequently, the multifunctional sites over FeO x /Rh-Fe catalysts highlight and meet the requirements of a selective methanol catalyst: strong interaction with CO 2 to ensure high density of transition states; metal sites to activate and make hydrogen available to surface intermediates; and high energy barriers for C-O bond cleavage to form carbides. These synthesis and catalytic chemistries, demonstrated for Rh-Fe-FeO x interfaces, enable us to overcome the limitations to the design of methanol production catalysts.

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

Structural Snapshots of Proteus vulgaris Tryptophan Indole-Lyase Reveal Insights into the Catalytic Mechanism

Tryptophan indole lyase (TIL; [E.C. 4.1.99.1]) is a bacterial pyridoxal-5′-phosphate (PLP)-dependent enzyme that catalyzes reversible β-elimination of indole from L-tryptophan. The mechanism of elimination of indole from L-tryptophan starts with the formation of an external aldimine of the substrate and PLP, followed by deprotonation of the α-CH of the substrate, forming a resonance-stabilized quinonoid intermediate. Proton transfer to C3 of the indole ring and carbon–carbon bond cleavage of the quinonoid intermediate provide indole and aminoacrylate bound to PLP, which then releases indole, followed by iminopyruvate. We have now determined the X-ray crystal structures of TIL complexes with (3S)-dioxindolyl-l-alanine, an inhibitor, and with substrates L-tryptophan, 7-aza-L-tryptophan, and S-ethyl-l-cysteine (SEC) in the presence of benzimidazole (BZI), an isostere of the product indole. These structures show a mixture of gem-diamine, external aldimine, quinonoid, and aminoacrylate intermediates, in both open and closed active site conformations. In the closed conformations of L-tryptophan, (3S)-dioxindolyl-l-alanine, and 7-aza-L-tryptophan complexes, hydrogen bonds form between Asp-133 with N1 of the ligand heterocyclic ring and NE2 of His-458 in the small domain of TIL. This hydrogen bond also forms in the BZI complex with the aminoacrylate intermediates formed from both L-tryptophan and SEC. The closed quinonoid complex of 7-aza-L-tryptophan shows that the azaindole ring in the closed conformation is bent out of plane of the Cβ–C3 bond by about 40°, putting it in a geometry that leads toward the transition-state geometry. Thus, both conformational dynamics and substrate activation play critical roles in the reaction mechanism of the TIL.

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

Catalytic Resonance Theory: Turnover Efficiency and the Resonance Frequency

Programmable catalysts exhibiting forced oscillation in the free energy of reacting surface species were simulated to understand the general mechanisms leading to efficient use of the input energy. Catalytic ratchets with either positive or negative adsorbate scaling exhibited oscillation conditions of both high and low turnover efficiency, yielding catalytic turnover frequencies either close to or significantly lower than the applied catalyst oscillation frequency, respectively. The “effective rate”, defined as the product of the catalytic turnover frequency and the turnover efficiency (ηTOE), was limited via two catalytic mechanisms: a leaky catalytic ratchet existed when molecules repeatedly traversed backward through the catalytic transition state upon catalyst oscillation, while a catalytic ratchet with low surface participation exhibited reduced formation of a gas-phase final product due to low surface product coverage. Furthermore, a single applied frequency yielding a maximum effective catalytic rate defined as the “resonance frequency” provided maximum combined benefit for catalytic rate and efficiency.

25 ENERGY STORAGE↗
DOE OSTI2025

Epoxide Alcoholysis over M-BEA Zeolites: Effects of Alcohol Chain Length on Rates and Regioselectivities

The structures of nucleophilic reactants affect their coordination behavior among solvent molecules and kinetics of reactions with surface intermediates within the confines of fluidfilled pores of zeolites and other microporous materials. Consequently, rates and regioselectivities of diverse chemistries may depend sensitively on nucleophile identity in manners not observed for classic fluid phase reactions. Here, we examine the impact of varying the primary alcohol (ROH) chain length on the kinetics of 1,2-epoxybutane (C 4 H 8 O) ring-opening within Brønsted (Al-BEA) and Lewis acid (Zr-BEA) zeolites. Turnover rates increase by factors of ∼6 (Al-BEA) and 4-fold (Zr-BEA) between methanol and 1-hexanol, yet the reaction mechanisms remain comparable. Despite modest rate differences, apparent activation enthalpies calculated from rates and activities of solvated reactants decrease linearly by 12 (Al-BEA) to 33 kJ mol −1 (Zr- BEA) with increased proton affinity, which suggests bond formation energies for the nucleophile strongly influence rate increases. The molecular interpretation of these trends demonstrates, however, that the solvation of ring-opening transition states by zeolite pore structures and solvent molecules also governs rates. The impact of local solvating interactions appears most directly as changes in regioselectivities, which tend to enhance terminal alcohol formation with increasing ROH chain length. Regioselectivities largely do not vary with differences in fluid composition for a given ROH. The addition of H 2 O increases the number of hydrogen bonds among reactive species, and trends in regioselectivities imply that the decreased hydrogen bonding ability of longer chain ROH, and not the nucleophile strength or steric bulk, determines the regioselectivities of the resulting products. This work provides direct experimental evidence that nucleophilicity and hydrogen bonding influence reaction barriers and regioselectivities in zeolitecatalyzed epoxide ring-opening, offering pathways to better control reaction kinetics.

acidic zeolites↗
DOE OSTI2026

Solvent Reorganization and Surface–Species Interactions Control Transfer Hydrogenation Rates in Sn-Beta Zeolites

Reorganization of confined alcohol networks reshapes adsorption and activation thermodynamics in Lewis acid zeolites, yet the roles of adsorbate sterics and non-H-bonding cosolvents remain unclear. Here, we investigate transfer hydrogenation of methyl- (mCH) and tert-butylcyclohexanone (tbCH) with 2-butanol over hydrophobic Sn-Beta. First-order rates per defect-open Sn site are 2–5× higher for mCH than tbCH, but converge in the zero-order regime (333–393 K), indicating that solvent reorganization entropically drives adsorption, while bond activation is insensitive to ketone structure. Reactions in toluene or acetonitrile (with 1 M 2-butanol) show that intraporous solvent organization governs transition-state stabilization: on hydrophobic Sn-Beta, apparent enthalpies and entropies rise with decreasing solvent polarity, whereas hydrophilic zeolites display nearly constant barriers across solvents. Liquid-phase IR spectroscopy corroborates these trends, showing that hydrophilic pores enforce 2-butanol-like solvation environments even in toluene. Together, these results reveal how zeolite pores regulate the intraporous solvent structure, whose reorganization upon adsorption reshapes the thermodynamic landscape of Lewis acid catalysis

Adsorption↗
DOE OSTI2026

Grand Challenges and Opportunities in Stimulated Dynamic and Resonant Catalysis

Traditional heterogeneous catalysis is constrained by kinetic and thermodynamic limits, such as the Sabatier principle and reaction equilibrium. Dynamic and resonant catalysts hold promise to overcome these limitations by actively oscillating a catalyst’s physical or electronic structure at the time scale of the catalytic cycle, allowing programmable control over reaction pathways, and leading to improved rate and selectivity. External stimuli such as temperature swing, mechanical strain, electric charge, and light can perturb catalyst surfaces in different ways, altering adsorbate coverage, binding energies, and transition states beyond what steady-state catalysis allows. This work surveys the current state of dynamic catalysis, introduces the concept of “ stimulando ” characterization for observing transient dynamics, and outlines key modeling, mechanistic, and benchmarking strategies to advance the field toward improved chemical transformation.

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

Confinement of Lewis Acid–Base Sites by Microporous Silica Layers on Titania for Enhanced Alkanol Dehydration Reactivity

Alkanol dehydration offers a pathway to upgrade biomass-derived short-chain oxygenates into alkenes, essential chemical building blocks widely used in industrial applications. Transition metal oxides with Lewis acid-base site pairs are attractive catalysts due to their high reactivity and cost-effectiveness. This work demonstrates a synthetic pathway to manipulate local environments around active Lewis acid-base pairs in anatase TiO 2 to enhance their reactivity in alkanol dehydration. Microporous SiO 2 layers with an average pore diameter of ~0.6 nm and a controlled thickness of 0.8-33 nm are deposited on anatase TiO 2 powders by using a molecular templated SiO 2 deposition method. The Lewis acid-base strength of accessible Ti-O pairs remains unchanged, as shown by temperature-programmed surface reactions of surface-bound formic acid-derived species and temperature-programmed desorption of pyridine. However, measured alkanol dehydration rates on confined Ti-O pairs are much higher (by up to 7-fold) than those on TiO 2 . The extent of rate enhancements depends on the reactant size and functional group positioning, suggesting that the rate enhancements reflect the interactions between the guest molecules (reactants and transition states) and the surrounding SiO 2 micropore environments. By providing a detailed synthetic procedure to tailor the local environments around active sites in bulk oxides, this approach offers an additional avenue for enhancing catalytic performance.

09 BIOMASS FUELS↗
DOE OSTI2026

Mechanistic Studies of an Iron-Catalyzed Intermolecular C–H Amination Reaction under Catalytic Conditions and Having a Large KIE

The conversion of C–H bonds into amines by nitrene insertion is an attractive transformation since it is both atom- and step-economical, and provides a direct route to functionalizing hydrocarbons. Using an iron catalyst [{( tBu pyrr) 2 pyr}Fe(OEt 2 )] (1-OEt 2 ) (( tBu pyrr) 2 pyr 2– = 3,5- t Bu 2 -bis(pyrrolyl)pyridine), we recently demonstrated the catalytic conversion of weak C–H bonds into secondary amines using aryl azides as the nitrene source [Zars, E.; Angew. Chem., Int. Ed. 2023, 62, e202311749]. Here, we describe detailed mechanistic studies of this intermolecular C–H amination reaction under catalytic conditions. We find by Variable Time Normalization Analysis (VTNA) that the conversion of xanthene (2-H 2 ) and 2,4,6-trimethyl-phenyl azide ( Me 3) catalyzed by 1-OEt 2 is an overall 3/2 order process, being 1 st order in 2-H 2 and half order in Me 3. A kinetic isotope effect study (KIE) using 2-d 2 results in a significant decrease in the rate (KIE = 61(15)), which clearly implicates the C–H insertion step as rate-determining. Furthermore, treatment of 1-OEt 2 with one equivalent of N 3 -2,6- i Pr 2 –C 6 H 3 yields the mixed-valence C–N coupled product [( tBu pyrr) 2 pyrFe-N═C(2,6 i Pr 2 –Ph)═N-(2,6 i Pr 2 –Ph))Fe tBu pyrrpyr(2-H-pyrr)] (5 iPr ). Quantum chemical calculations confirm the electronic structure of the mixed-valence dimer in 5 iPr and rationalize the Hammett correlation by a delicate balance in the dinuclearization of the catalytically active monomers. Calculations further indicate significant tunneling for the pivotal H atom abstraction by the iron-imidyl complex. Combining all these results allows us to propose a mechanism consisting of imido formation in equilibrium with a radical-coupled diiron system, followed by stepwise C–H insertion via a linear H atom abstraction transition state and subsequent radical rebound.

azides↗
DOE OSTI2026

Regio- and Stereoselective Lactone Polymerization: Divergent Effect of Catalyst Modification and Monomer Structure

Selective ring-opening polymerization (ROP) of chiral lactones enables access to biodegradable polyesters with precisely controlled microstructures. Here, combined DFT modeling and experimental validation elucidate how fine-tuning of enantiopure SalBinam aluminum catalysts (through introduction of bromine atoms and tert-butyl groups, respectively, in ortho- and para-positions) modulates regio- and stereocontrol in the ROP of methyl glycolide (MeG) and lactide (LA). Computations reveal that regioselectivity in MeG polymerization arises mainly from steric repulsion, with a small contribution from weak stabilizing C-H···Br interactions that favor ring-opening at the glycolic site, consistent with the experimentally enhanced regioselectivity for (R)-MeG. In contrast, the same steric congestion at the ligand’s ortho positions destabilizes the key transition states in rac-LA polymerization, reducing the calculated stereoselectivity. Experiments confirm the predicted loss of stereocontrol, yielding nearly atactic PLA under standard conditions. Extension of the computational framework to rac-MeG polymerization promoted by racemic catalyst identified a low-barrier, stepwise polymer chain exchange pathway that rationalizes the experimentally observed syndiotacticity of poly(lactic-co-glycolic acid).

DFT calculations↗
DOE OSTI2020

Demonstrating the Critical Role of Solvation in Supported Ti and Nb Epoxidation Catalysts via Vapor-Phase Kinetics

Catalytic oxidation of hydrocarbons with hydrogen peroxide (H 2 O 2 ) has been of the utmost importance for several decades. The vast majority of studies have been performed in the condensed phase, even though condensed phases introduce complex solvent effects and can promote the leaching of active sites. In response, we have built a custom reactor system to understand H 2 O 2 activation and selective oxidation in the vapor-phase. In this report, we study the epoxidation of cyclohexene with H 2 O 2 over four Lewis-acidic metal oxide catalysts: Ti and Nb grafted on SiO 2 and on the Zr based metal–organic framework, NU-1000. The M–SiO2 materials are highly selective to the formation of epoxides and diols, as they can be in the condensed phase, while the NU-1000 based materials are far more prone to overoxidation to CO 2 , which appears to be connected to their strong reactant adsorption. Apparent activation energies are calculated for all materials when operating in the same kinetic regime, and the heats of cyclohexene adsorption into their pores are then used to directly compare intrinsic enthalpies of activation in the vapor vs condensed phase for the M–SiO 2 catalysts. Nb–SiO 2 catalysts exhibit similar intrinsic enthalpies of activation in the vapor and condensed phases, whereas the condensed phase transition state in Ti–SiO 2 is 24 kJ/mol lower in energy than that of the same material in the vapor phase. These experiments establish another methodology for understanding the various roles of solvent in selective oxidation reactions and studying these reactions under conditions that differ significantly from the thousands of prior studies in the condensed phase.

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

Interfacial Reactivity and Speciation Emerging from Na-Montmorillonite Interactions with Water and Formic Acid at 200 °C: Insights from Reactive Molecular Dynamics Simulations, Infrared Spectroscopy, and X-ray Scattering Measurements

Reactive organic fluid - mineral interactions at elevated temperatures contribute to the evolution of planetary matter. One of the less studied but important transformations in this regard involves the reactions of formic acid with naturally occurring clays such as sodium montmorillonite. To advance a mechanistic understanding of these interactions, we use ReaxFF reactive molecular dynamics simulations in conjunction with infrared (IR) spectroscopy and X-ray scattering experiments to investigate the speciation behavior of water-formic acid mixtures on sodium montmorillonite interfaces at 473 K and 1 atm. Using a newly developed reactive forcefield, we show that the experimental IR spectra of unreacted and reacted mixture can be accurately reproduced by ReaxFF/MD. We further benchmark the simulation predictions of sodium carbonate and bicarbonate formation in the clay interlayers using Small and Wide-Angle X-ray Scattering measurements. Subsequently, leveraging the benchmarked forcefield, we interrogate the pathway of speciation reactions with emphasis on carbonate, formate, and hydroxide groups elucidating the energetics, transition states, intermediates, and preferred products. Further, we also delineate the differences in reactivities and catalytic effects of clay edges, facets, and interlayers owing to their local chemical environments, which have far reaching consequences in their speciation behavior. The experimental and simulation approaches described in this study and the transferable forcefields can be applied translationally to advance the science of clay-fluid interactions for several applications including subsurface fluid storage and recovery and clay-pollutant dynamics.

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

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

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

ab initio calculations↗