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

Automated theoretical chemical kinetics: Predicting the kinetics for the initial stages of pyrolysis

Large scale implementation of high level computational theoretical chemical kinetics offers the prospect for dramatically improving the fidelity of combustion chemical modeling. To facilitate such efforts, we have developed a suite of codes, collectively referred to as AutoMech, that allow for the automatic prediction of the kinetics for large sets of reactions via ab initio transition-state-theory based master-equation calculations. The primary input is simply the mechanism, a dictionary relating chemically identifiable species descriptors (e.g., SMILES or InChIs) to species labels in the mechanism, and a specification of the electronic structure and transition state theory models to be implemented. Here we illustrate the current utility of AutoMech through a study of the initial stages of pyrolysis for 3 sets of fuels: sequences of alkanes, alcohols, and aldehydes. For simplicity, the analysis focuses on abstractions from the fuel by H, CH 3 , and OH, and the decomposition of the resulting radicals. Altogether, there are a total of 166 input channels in these sets (more than 363 forward reactions when expanded to the full set of elementary reactions). The code successfully produces high quality rate estimates (with apparent uncertainties less than a factor of two in limited comparisons with experiment) for > 95% of these. For the radical decomposition reactions, the analysis includes predictions for the pressure dependence of the kinetics. This wide-ranging exploration illustrates (i) the effect of different levels of prediction on the expected accuracy, (ii) the branching between abstractions at different sites for different abstractors, (iii) the dependence of the rates on the chemical structure, and (iv) the variation in radical stabilities across chemical families. These results, as well as the demonstrated feasibility of the methodology, should find further utility in the development of accurate rate expressions for arbitrary fuels. (c) 2020 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Elliott, Sarah N.↗
DOE OSTI2021

Theoretical kinetics predictions for NH 2 + HO 2

Recent modeling studies of NH 3 oxidation, which are motivated by the prospective role of ammonia as a zero-carbon fuel, have indicated significant discrepancies between existing literature mechanisms. In this study high level theoretical kinetics predictions have been obtained for the reaction of NH 2 with HO 2 , which has previously been highlighted as an important reaction with high sensitivity and high uncertainty. The potential energy surface is explored with coupled cluster calculations including large basis sets and high-level corrections to yield high accuracy (~0.2 kcal/mol) estimates of the stationary point energies. Variational transition state theory is used to predict the microcanonical rate constants, which are then incorporated in master equation treatments of the temperature and pressure dependent kinetics. For the radical-radical channels, the microcanonical rates are obtained from variable reaction coordinate transition state theory implementing directly evaluated multireference electronic energies. The analysis yields predictions for the total rate constant as well as the branching to the NH 3 + O 2 , H 2 NO + OH, and HNO + H 2 O channels. Rate constants are also reported for the H 2 NO + OH reaction as they arise naturally from the analysis. The rate constant and branching fraction determined in this work for the NH 2 + HO 2 reaction deviate significantly from values used in most previous modeling studies. The fact that the main product channel is chain terminating, rather than propagating, has strong implications for modeling NH 3 ignition and oxidation, in particular at intermediate temperatures and elevated pressure.

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

Transition Path Sampling Study of the Feruloyl Esterase Mechanism

Serine hydrolases cleave peptide and ester bonds and are ubiquitous in nature, with applications in biotechnology, in materials, and as drug targets. The serine hydrolase two-step mechanism employs a serine–histidine–aspartate/glutamate catalytic triad, where the histidine residue acts as a base to activate poor nucleophiles (a serine residue or a water molecule) and as an acid to allow the dissociation of poor leaving groups. This mechanism has been the subject of debate regarding how histidine shuttles the proton from the nucleophile to the leaving group. To elucidate the reaction mechanism of serine hydrolases, we employ quantum mechanics/molecular mechanics-based transition path sampling to obtain the reaction coordinate using the Aspergillus niger feruloyl esterase A (AnFaeA) as a model enzyme. The optimal reaction coordinates include terms involving nucleophilic attack on the carbonyl carbon and proton transfer to, and dissociation of, the leaving group. During the reaction, the histidine residue undergoes a reorientation on the time scale of hundreds of femtoseconds that supports the “moving histidine” mechanism, thus calling into question the “ring flip” mechanism. We find a concerted mechanism, where the transition state coincides with the tetrahedral intermediate with the histidine residue pointed between the nucleophile and the leaving group. Moreover, motions of the catalytic aspartate toward the histidine occur concertedly with proton abstraction by the catalytic histidine and help stabilize the transition state, thus partially explaining how serine hydrolases enable poor nucleophiles to attack the substrate carbonyl carbon. Rate calculations indicate that the second step (deacylation) is rate-determining, with a calculated rate constant of 66 s –1 . Altogether, these results reveal the pivotal role of active-site dynamics in the catalytic mechanism of AnFaeA, which is likely similar in other serine hydrolases.

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

A Conical Intersection Influences the Ground State Rearrangement of Fulvene to Benzene

The rearrangement of fulvene to benzene is believed to play an important role in the formation of soot during hydrocarbon combustion. Previous work has identified two possible mechanisms for the rearrangement-a unimolecular path and a hydrogen-atom-assisted, bimolecular path. Computational results to date have suggested that the unimolecular mechanism faces a barrier of about 74 kcal/mol, which makes it unable to compete with the bimolecular mechanism under typical combustion conditions. This computed barrier is about 10 kcal/mol higher than the experimental value, which is an unusually large discrepancy for modern electronic structure theory. In the present work, we have reinvestigated the unimolecular mechanism computationally, and we have found a second transition state that is approximately 10 kcal/mol lower in energy than the previously identified one and, therefore, in excellent agreement with the experimental value. The existence of two transition states for the same rearrangement arises because there is a conical intersection between the two lowest singlet states which occurs in the vicinity of the reaction coordinates. The two possible paths around the cone on the lower adiabatic surface give rise to the two distinct saddle points. The lower barrier for the unimolecular mechanism now makes it competitive with the bimolecular one, according to our calculations. In support of this conclusion, we have reanalyzed some previous experimental results on anisole pyrolysis, which leads to benzene as a significant product and have shown that the unimolecular and bimolecular mechanisms for fulvene ? benzene must be occurring competitively in that system. Finally, we have identified that similar conical intersections arise during the isomerizations of benzofulvene and isobenzofulvene to naphthalene.

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

Kinetic analyses for solid-state phase transition of metastable amorphous-AlO x (2.5 < x ≤ 3.0) nanostructures into crystalline alumina polymorphs

Solid-solid phase change materials (SS-PCMs) hold promise for energy storage/dissipation in batteries and energetic materials. Yet, phase change kinetics for SS-PCMs undergoing metastable to semi-stable/stable phase transformations remain relatively ill-studied because trapping metastable phases remain challenging. Recently, we demonstrated the kinetic entrapment and stabilization of a highly disordered and amorphous Al-oxide phase m-AlO x @C (x~2.5-3.0) via laser ablation synthesis in solution (LASiS). We report here, to our knowledge, the first chemical kinetics analysis for S-S phase transition of the m-AlO 3 @C nanocomposites (< 5–8 nm sizes) into semi-stable equilibrium alumina phases (θ/γ-Al 2 O 3 ) via disproportionation reaction, while releasing excess trapped gases. Our results indicate the atomic density of the AlO 3 structures to be ~5–10 times less than that of the final Al 2 O 3 phases, which led to the hypothesis of a volume shrinkage process during their phase transition. Temperature-dependent X-ray diffraction studies reveal the high-temperature phase transition for m-AlO 3 → θ/γ-Al 2 O 3 to follow contracting volume kinetics model, thereby validating our earlier hypothesis. Using the geometric volume contraction model, reaction kinetics analyses from Arrhenius plots reveal the activation energy barrier for the phase transition to be ~270±11 kJ/mol. This makes the activation energy barrier nearly identical to the oxidation of micron-sized Al particles.

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

Facet-Dependent Water Inhibition of Alkanol Dehydration on TiO2 via Distinct Water–Alkanol Complexes

Water is ubiquitous in biomass-derived feeds, yet its molecular impact on oxygen-elimination reactions remains poorly understood, particularly for catalysts exposing different facets. Here, we utilize well-defined TiO2 nanocrystals with dominant (101) and (001) facets to reveal a pronounced facet-dependent effect of water, where inhibition for dehydration of isopropanol (IPA) on the TiO2(001) surface is about four times more severe than TiO2(101). Through a combination of in situ solid state NMR, in situ infrared spectroscopy, kinetics studies, and theoretical calculations, we demonstrate that this disparity arises from the formation of distinct alkanol-water complex intermediates. On TiO2(001), IPA undergoes dissociative adsorption to form an isopropoxide-H2O complex that readily drives the surface into a complex-dominated regime. This pathway increases the activation barrier for C–H cleavage by 40 kJ mol-1 by inducing a disordered transition state. In contrast, TiO2(101) favors molecular IPA adsorption with weak hydrogen bonding to water, resulting in a smaller complex formation constant and a much smaller activation barrier increase (25 kJ mol-1). By quantitatively linking facet-dependent complex coverage to transition-state destabilization, this work moves beyond simple site-blocking models and provides a conceptual framework for designing catalysts that remain active in water-containing environments

Hu, Wenda↗
DOE OSTI2024

Demystifying group-4 polyolefin hydrogenolysis catalysis: Gaseous propane hydrogenolysis mechanism over the same catalysts

A kinetic/mechanistic investigation of gaseous propane hydrogenolysis over the single-site heterogeneous polyolefin depolymerization catalysts AlS/ZrNp 2 and AlS/HfNp 2 (AlS = sulfated alumina, Np = neopentyl), is use to probe intrinsic catalyst properties without the complexities introduced by time- and viscosity-dependent polymer medium effects. In a polymer-free automated plug-flow catalytic reactor, propane hydrogenolysis turnover frequencies approach 3,000 h −1 at 150 °C. Both catalysts exhibit approximately linear relationships between rate and [H 2 ] at substoichiometric [H 2 ] with rate law orders of 0.66 ± 0.09 and 0.48 ± 0.07 for Hf and Zr, respectively; at higher [H 2 ], the rates approach zero-order in [H 2 ]. Reaction orders in [C 3 H 8 ] and [catalyst] are essentially zero-order under all conditions, with the former implying rapid, irreversible alkane binding/activation. This rate law, activation parameter, and DFT energy span analysis support a scenario in which [H 2 ] is pivotal in one of two plausible and competing rate-determining transition states—bimolecular metal-alkyl bond hydrogenolysis vs. unimolecular β-alkyl elimination. The Zr and Hf catalyst activation parameters, ΔH ‡ = 16.8 ± 0.2 kcal mol −1 and 18.2 ± 0.6 kcal mol −1 , respectively, track the relative turnover frequencies, while ΔS ‡ = −19.1 ± 0.8 and −16.7 ± 1.4 cal mol −1 K −1 , respectively, imply highly organized transition states. These catalysts maintain activity up to 200 °C, while time-on-stream data indicate multiday activities with an extrapolated turnover number ~92,000 at 150 °C for the Zr catalyst. This methodology is attractive for depolymerization catalyst discovery and process optimization.

03 NATURAL GAS↗
DOE OSTI2026

Transient hydroperoxyalkyl intermediates (•QOOH) in isopentane oxidation. II. Isomer-resolved unimolecular dynamics

Transient carbon-centered hydroperoxyalkyl intermediates (•QOOH) in isopentane oxidation are characterized by their time- and energy-resolved unimolecular dissociation dynamics to hydroxyl (OH) and cyclic ether products. Two distinct •QOOH isomers are examined with radical sites at a primary carbon of one of the methyl groups (β-Me) or a secondary carbon (β-Et) of the ethyl group. Energy-dependent unimolecular rates are obtained from the time-dependent appearance of OH products for the two isomers and compared with statistical microcanonical rates computed using RRKM theory, including heavy-atom tunneling, based on high-level theoretical calculations. A benchmark-corrected approach is utilized to compute high-accuracy stationary-point energies, most importantly, transition-state barriers, for the •QOOH_Me and •QOOH_Et isomers in isopentane oxidation, building on higher-level reference calculations for the oxidation of ethane (C 2 H 5 O 2 ) and propane (C 3 H 7 O 2 ), respectively. Here, the measured rates are compared with RRKM calculations incorporating the benchmark-corrected transition-state parameters, a vibrationally adiabatic multidimensional hindered-rotor treatment of key torsions, and quantum tunneling. Agreement between experiment and theory validates the statistical description and shows faster decay for •QOOH_Et due to its lower barrier. Both β-QOOH isomers decay almost exclusively to OH + cyclic ether products under the conditions studied.

Oxidation processes↗
DOE OSTI2026

Nuclear–Electronic Orbital General Rate Theory: Predicting Hydrogen Kinetic Isotope Effects in the Deep Tunneling Regime

Hydrogen transfer is a critical component of many chemical and biological processes. The ratio of rate constants for hydrogen and deuterium transfer defines the H/D kinetic isotope effect (KIE), which is a powerful tool for elucidating hydrogen transfer mechanisms. Interpretation of experimental H/D KIEs relies on accurate and affordable computational methods. However, due to their light mass, hydrogen and deuterium can undergo tunneling, which is challenging to describe in multidimensional molecular systems. Herein, we introduce the nuclear–electronic orbital general rate theory (NEO-GRT), which enables the efficient prediction of H/D KIEs based on full-dimensional molecular quantum chemistry calculations. The NEO-GRT approach describes the hydrogen transfer rate constant with a general expression that spans the vibrationally adiabatic and nonadiabatic hydrogen tunneling regimes. The input quantities are computed using NEO density functional theory, which treats the transferring hydrogen or deuterium nucleus quantum mechanically on the same level as the electrons. We investigate two intramolecular proton transfer reactions in organic molecules at temperatures down to 50 K to evaluate the performance of NEO-GRT by comparison to transition state theory and ring-polymer instanton theory. The KIEs computed with NEO-GRT agree with those calculated using ring-polymer instanton theory for the full-dimensional molecular systems at the same level of electronic structure theory. This agreement indicates that NEO-GRT captures the deep hydrogen tunneling effects, in contrast to transition state theory, which neglects such effects. Given its relatively low computational cost, NEO-GRT is a promising approach for predicting H/D KIEs in large organic and organometallic systems.

Hydrogen↗
DOE OSTI2021

Role of Hydrogen-bonded Bimolecular Formic Acid–Formate Complexes for Formic Acid Decomposition on Copper: A Combined First-Principles and Microkinetic Modeling Study

Hydrogen bonding interactions alter the nanoscale reaction mechanisms of many chemistries. Yet, it remains unclear how they affect heterogeneously catalyzed decomposition of formic acid (FA), a reaction of intense interest since FA is a promising hydrogen carrier. In this work, we elucidate how hydrogen bonding affects the reaction mechanisms for FA decomposition on Cu(111) by combining first-principles density functional theory calculations to calculate reaction energetics, Latin-hypercube sampling to elucidate stable high-coverage adsorbate configurations, and coverage self-consistent mean-field microkinetic models to predict reaction kinetics. We demonstrate that hydrogen-bonded complexes of FA with formate (bimolecular FA–HCOO complexes) can play a dominant role in FA decomposition. Specifically, our first-principles calculations show that hydrogen bonding of FA with HCOO may stabilize the crucial monodentate HCOO intermediate and the transition states for HCOO decomposition, especially at low coverages. We predict that, depending on the reaction conditions, 40–80% of the reaction flux goes through pathways involving the bimolecular FA–HCOO complexes. Additionally, the active site for FA decomposition on Cu(111) involves a high coverage (~0.4 monolayers (ML)) of these complexes, which unexpectedly stabilize intermediates and transition states via van der Waals interactions. Our work provides molecular insights consistent with previous experimental observations on supported Cu/Al 2 O 3 catalysts. This paves the way toward the development of novel catalysts for FA decomposition as well as for other industrially important chemistries with intermediates capable of hydrogen bonding, such as ammonia electrooxidation and CO 2 hydrogenation.

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

Machine Learning Accelerated First-Principles Study of the Hydrodeoxygenation of Propanoic Acid

The complex reaction network of catalytic biomass conversions often involves hundreds of surface intermediates and thousands of reaction steps, greatly hindering the rational design of metal catalysts for these conversions. Here, we present a framework of machine learning (ML)-accelerated first-principles studies for the hydrodeoxygenation (HDO) of propanoic acid over transition metal surfaces. The microkinetic model (MKM) is initially parametrized by ML-predicted energies and iteratively improved by identifying the rate-determining species and steps (RDS), computing their energies by density functional theory (DFT), and reparameterizing the MKM until all the RDS are computed by DFT. The Gaussian process (GP) model performs significantly better than the linear ridge regression model for predicting both the adsorption free energies and transition state free energies. Parameterized with energies from the GP model, only 5–20% of the full reaction network has to be computed by DFT for the MKM to possess DFT-level accuracy for the TOF and dominant reaction pathway. While the linear ridge regression model performs worse than the GP model, its performance is greatly improved when only transition states are predicted by the regression model and adsorption energies are computed by DFT. Overall, we find that a high accuracy in adsorption free energies is more important for a reliable MKM than a high accuracy in TS free energies. Lastly, based on the GP model with GOH and GCHCHCO as catalyst descriptors, we build two-dimensional volcano plots in activity and selectivity that can help design promising alloy catalysts for HDO reactions of organic acids.

adsorption↗
DOE OSTI2026

Selective Chemical Looping Combustion of Terminal Alkynes in Mixtures with Alkenes

The selective combustion of terminal alkynes in mixtures with alkenes is demonstrated during anaerobic reduction half-cycles on bulk bismuth oxide (Bi 2 O 3 ) as an approach to remove alkynes, which act as inhibitors in olefin polymerization. Bi 2 O 3 combusts phenylacetylene in styrene, 3-methylphenylacetylene in 3-methylstyrene, propyne in propylene, 1-hexyne in 1-hexene, and 1-octyne in 1-octene, with alkyne combustion selectivities exceeding 96%. Near unity reaction orders for hydrocarbon consumption during reduction half-cycles are consistent with combustion pathways initiated by rate-determining initial C–H activation, which drive selective alkyne combustion through intrinsic differences in the first-order rate constants for alkyne and alkene combustion rather than preferential adsorption of alkynes on Bi 2 O 3 surfaces. Computational assessments of initial C–H activation pathways for alkynes and alkenes on (010) α-Bi 2 O 3 surfaces using density functional theory illustrate that heterolytic transition states which form proton-carbanion pairs on Bi–O sites kinetically favor the activation of alkynes rather than alkenes due to differences in C–H bond acidity, and the barrier for heterolytic C–H activation is dictated in part by the sum of the molecular deprotonation energy and the energy to bind an R – carbanion to a Bi site in its transition-state geometry. Finally, these heterolytic reactivity channels during selective chemical looping combustion present novel routes for purifying olefin gas streams containing alkyne impurities.

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

Probing the dynamics and bottleneck of the key atmospheric SO 2 oxidation reaction by the hydroxyl radical

SO 2 (Sulfur dioxide) is the major precursor to the production of sulfuric acid (H 2 SO 4 ), contributing to acid rain and atmospheric aerosols. Sulfuric acid formed from SO 2 generates light-reflecting sulfate aerosol particles in the atmosphere. This property has prompted recent geoengineering proposals to inject sulfuric acid or its precursors into the Earth’s atmosphere to increase the planetary albedo to counteract global warming. SO 2 oxidation in the atmosphere by the hydroxyl radical HO to form HOSO 2 is a key rate-limiting step in the mechanism for forming acid rain. However, the dynamics of the HO + SO 2 → HOSO 2 reaction and its slow rate in the atmosphere are poorly understood to date. Herein, we use photoelectron spectroscopy of cryogenically cooled HOSO 2 – anion to access the neutral HOSO 2 radical near the transition state of the HO + SO 2 reaction. Spectroscopic and dynamic calculations are conducted on the first ab initio-based full-dimensional potential energy surface to interpret the photoelectron spectra of HOSO 2 – and to probe the dynamics of the HO + SO 2 reaction. In addition to the finding of a unique pre-reaction complex (HO⋯SO 2 ) directly connected to the transition state, dynamic calculations reveal that the accessible phase space for the HO + SO 2 → HOSO 2 reaction is extremely narrow, forming a key reaction bottleneck and slowing the reaction rate in the atmosphere, despite the low reaction barrier. This study underlines the importance of understanding the full multidimensional potential energy surface to elucidate the dynamics of complex bimolecular reactions involving polyatomic reactants.

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

Combined Experimental and Computational Efforts to Establish Ion Mobility, Solubility and Stability of Functional Liquids for Electrochemical Energy Storage

This work provides a computation-driven investigation of the stability of organic electrolytes for lithium-air batteries. Electrolyte instability is currently a key challenge that limits practical use of aprotic Li-air batteries, and the chemical processes that cause this instability are often kinetically-driven. Computational screening for kinetic stability involves the determination of reaction barriers for the numerous potential reaction mechanisms, barriers that are challenging to calculate due to the difficulty of locating transition state structures. Here we screen a broad set of substituted electrolytes for susceptibility to nucleophilic attack by superoxide. We find that carbonates are not typically expected to be stable and that sulfones are generally stable, validating literature trends. We study the effects of chemical functionalization with electron-donating and withdrawing groups and their interplay with steric factors, identifying functional groups and other chemical modifications that increase stability in these groups. User-input driven transition state identification is used for these initial calculations, and an automated computational pipeline is subsequently presented and validated as a means to perform further high-throughput searches across mechanisms and chemistries. The pipeline integrates cheminformatics-based reaction encoding, relaxed potential energy scans, and nudged elastic band calculations for an end-to-end approach to barrier calculations. We review this automated search approach and its current limitations, and discuss challenges and further work.

25 ENERGY STORAGE↗
DOE OSTI2022

High accuracy barrier heights, enthalpies, and rate coefficients for chemical reactions

Abstract Quantitative chemical reaction data, including activation energies and reaction rates, are crucial for developing detailed kinetic mechanisms and accurately predicting reaction outcomes. However, such data are often difficult to find, and high-quality datasets are especially rare. Here, we use CCSD(T)-F12a/cc-pVDZ-F12// ω B97X-D3/def2-TZVP to obtain high-quality single point calculations for nearly 22,000 unique stable species and transition states. We report the results from these quantum chemistry calculations and extract the barrier heights and reaction enthalpies to create a kinetics dataset of nearly 12,000 gas-phase reactions. These reactions involve H, C, N, and O, contain up to seven heavy atoms, and have cleaned atom-mapped SMILES. Our higher-accuracy coupled-cluster barrier heights differ significantly (RMSE of ∼5 kcal mol −1 ) relative to those calculated at ω B97X-D3/def2-TZVP. We also report accurate transition state theory rate coefficients $${k}_{\infty }(T)$$ k ∞ ( T ) between 300 K and 2000 K and the corresponding Arrhenius parameters for a subset of rigid reactions. We believe this data will accelerate development of automated and reliable methods for quantitative reaction prediction.

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

Resolving the Coverage Dependence of Surface Reaction Kinetics with Machine Learning and Automated Quantum Chemistry Workflows

Microkinetic models for catalytic systems require estimation of many thermodynamic and kinetic parameters that can be calculated for isolated species and transition states using ab initio methods. However, the presence of nearby coadsorbates on the surface can dramatically alter these thermodynamic and kinetic parameters causing them to be dependent on species coverage fractions. As there are combinatorially many coadsorbed configurations on the surface, computing the coverage dependence of these parameters is far less straightforward. We present a framework for generating and applying machine learning models to predict coverage-dependent parameters for microkinetic models. Our toolkit enables automatic calculation and evaluation of coadsorbed configurations allowing us to sample 2,000 coadsorbed adsorbates and transition states (TSs) for a diverse set of 9 reactions on Cu(111), a challenging surface, with four possible coadsorbates. This dataset was then used to train subgraph isomorphic decision trees (SIDTs) to predict the stability and association energy of configurations. We were able to achieve mean absolute errors (MAEs) of 0.106 eV on adsorbates, 0.172 eV on TSs, and due to natural error cancellation in SIDTs for relative properties, 0.130 eV on reaction energies and 0.180 eV on activation barriers. In conclusion, we describe how to use these models to predict coverage-dependent corrections for adsorbates and TSs and demonstrate on H*, HO*, and O* comparing the generated SIDT model with an iteratively refined version.

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

Controlled Placement of Trivalent Heteroatoms in MFI Zeolite Frameworks Using Structure-Directing Agents

MFI zeolites contain distinct confining pore environments, either smaller (∼0.55 nm) channels or larger (∼0.70 nm) intersections. The substitution of trivalent heteroatoms generates Brønsted acid sites of varying acid strength, which, together with their location within different pore environments, dictates the stability of transition states via electrostatic (ion-pair) and van der Waals interactions. Here, we report that structure-directing agents (SDAs) that control Al 3+ substitutional patterns during MFI crystallization analogously position other trivalent heteroatoms (Ga 3+ , Fe 3+ , B 3+ ), altering the distribution of acid sites between channels and intersections, as probed by low-temperature (403 K) toluene methylation kinetics. Heteroatom-substituted MFI crystallized with tetrapropylammonium result in low selectivity toward p-xylene (<30%), while those crystallized using (co-)SDAs containing peripheral hydrogen-bonding groups (e.g., ethylenediamine) show high p-xylene selectivity (∼80%). These selectivities are consistent with DFT-calculated Gibbs free energy barriers for intersection-dominant and channel-dominant active site distributions, respectively. Transition states to form each xylene isomer are similar in size and charge; thus, their rate constants decrease with acid strength similarly, causing isomer selectivity to depend strongly on confinement but not on acid strength. These generalizable synthetic strategies enable independently controlling acid site strength and location in MFI zeolites and, in turn, catalytic rates and selectivities.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Role of the ionic environment in enhancing the activity of reacting molecules in zeolite pores

Tailoring the molecular environment around catalytically active site allows the enhancement of catalytic reactivity via a hitherto unexplored pathway. In zeolites, the presence of water creates an ionic environment via the formation of hydrated hydronium ions and the negatively charged framework Al tetrahedra. The high density of cation-anion pairs determined by the aluminum concentration of a zeolite induces a high local ionic strength that increases the excess chemical potential of sorbed and uncharged organic reactants. Charged transition states (as for example, the carbenium ions in the discussed alcohol dehydration) are stabilized, reducing the energy barrier and leading to a higher reaction rate. Using the intramolecular dehydration of cyclohexanol on H-MFI in water, we show quantitatively the enhancement of the reaction rate by the presence of high ionic strength as well as potential limitations of this strategy. The approach opens a new pathway to systematically enhance catalytic reactivity rates and has wide-ranging implications for understanding catalysis in condensed phase.

Pfriem, Niklas↗