Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Microkinetic Modeling”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Impact of microkinetic modeling assumptions on predicted kinetics and mechanisms over undercoordinated sites

Accurate modeling of catalytic reactions on undercoordinated sites requires accounting for the structural and ensemble-specific nature of the active sites. This study examines how common microkinetic modeling (MKM) assumptions affect predicted kinetics and mechanisms on the stepped Pt(211) facet for the ethane dehydrogenation (EDH) and the ethane hydrogenolysis (EH). Six (211) MKMs were developed, differing in (i) the number of active sites represented, (ii) adsorbate site occupancy treatment, and (iii) inclusion of cross-facet interactions. These models are benchmarked against a particle-based microkinetic model (PB-MKM), which best represents step-edge behavior. MKM assumptions caused deviations in turnover frequencies exceeding ten orders of magnitude and led to contrasting mechanistic and selectivity predictions. Multi-site MKMs overestimate activity by inflating free site availability, single-site models underestimate activity, and uniform occupancy models overpredict coverage of multi-dentate intermediates, leading to reaction-specific artifacts. Overall, the Combined Site Edge Model (CSEM), a single-site MKM accounting for site occupancy and cross-facet interactions, most closely approximates PB-MKM predictions. All models predict similar kinetics when surfaces are clean or primarily occupied by monodentate species. This work provides practical guidance for selecting MKM frameworks for undercoordinated catalytic surfaces and highlights the critical role of modeling assumptions in catalytic predictions.

(211) facet↗
DOE OSTI2026

A Unified Workflow for Sensitivity-Based Kinetic Analysis in Microkinetic Models

Degrees of rate control (DRC), apparent activation energies, and apparent reaction orders are established local sensitivity diagnostics for interpreting microkinetic models, but applying them routinely to large mechanisms often requires substantial reaction-specific bookkeeping, perturbation design, and postprocessing. Here, in this study, we present a unified derivative-based workflow that evaluates these quantities from a single compiled reaction-network model and target-rate definition. For any user-provided microkinetic model, the workflow compiles the mechanism into stoichiometrically consistent mass-action rate equations, solves the surface dynamics, and uses automatic differentiation to compute sensitivities with respect to rate constants, temperature, and gas partial pressures. By combining their calculations in the same framework, the workflow clearly demonstrates the relationships between different DRCs and the apparent activation energy. Using existing examples of propylene partial oxidation and methane oxidation on Pd(100), we verify expected transient redistribution of rate control, distinguish net Campbell DRCs from one-sided directional sensitivities, and show how apparent activation energy can be reconstructed either from one-sided DRCs or from state-based DRCs while critical mechanistic insights are obtained consistently. In the methane oxidation case, a pathway-subset test further illustrates how a simplified mechanism preserves key kinetic signatures of a full model, showing the potential of our user-friendly tool for model construction beyond kinetic analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Pore-Scale Transport Effects in Electrochemical CO 2 Reduction on Gold via Coupled Microkinetic-Transport Modeling

A pore-resolved modeling framework is developed to quantify how pore-scale transport affects the intrinsic microkinetics of CO 2 -to-CO on Au. A DFT-informed microkinetic model is coupled self-consistently to a Generalized-Modified Poisson–Nernst–Planck (GMPNP) transport description in a single, electrolyte-filled cylindrical pore, allowing local concentrations and potential to feed back into site-specific reaction rates. FIB-SEM is used to determine pore sizes within realistic electrode materials. Across pore diameters, d p = 10–6000 nm, the surface-averaged CO 2 reduction rate is systematically reduced relative to the ideal microkinetic baseline where mass transport is not accounted for; the effectiveness factor 𝜂 𝑠,CO 2 , which quantifies this ratio, decreases rapidly at more negative potentials and is about 1% near −1.0 V vs SHE due to reactant depletion. Spatial maps reveal pore-bulk alkalization that emerges at higher cathodic bias, with a small, near-wall pH dip due to electrostatic repulsion of hydroxide at the cathode interface. For a fixed aspect ratio L p /d p , narrower pores exhibit larger 𝜂 𝑠,CO 2 by shortening diffusion paths, whereas variations in the aspect ratio L p /d p play a secondary role. A dimensionless analysis (surface/bulk Damköhler numbers) delineates operating regimes. In conclusion, this work offers a concept for incorporating microkinetic models into homogenized porous-electrode models through effectiveness factors and pore-size distribution.

Au-catalyst↗
DOE OSTI2025

Enhanced Selectivity for C 2 H 4 Production from C 2 H 6 on Partially Chlorinated IrO 2 (110) Surfaces

Modifying metal oxide surfaces to limit their oxidizing activity can provide a means of improving catalytic selectivity toward the partial oxidation of light alkanes. Here, in this study, we investigated the oxidation of C 2 H 6 on Cl-modified IrO 2 (110) surfaces using temperature-programmed reaction spectroscopy (TPRS) and first-principles microkinetic modeling. We find that substituting Cl for O in the IrO 2 (110) surface enhances the selectivity for C 2 H 6 conversion to C 2 H 4 during TPRS by suppressing extensive oxidation to CO x products, while also either enhancing C 2 H 4 production or altering it to a lesser extent, depending on the initial C 2 H 6 coverage. The C 2 H 4 selectivity increased with increasing C 2 H 6 and Cl coverage, but reached a limiting value below 50%. The Cl coverage changed negligibly during C 2 H 6 oxidation, and the surface reactivity decreased only marginally for Cl coverages up to 0.5 ML (monolayer). TPRS simulations using a microkinetic model predict C 2 H 4 and CO x product yields as a function of the Cl coverage that agree closely with the experimental results. According to the simulations, C 2 H 6 conversion to C 2 H 4 occurs on Cl-IrO 2 (110) by the hydrogenation of C 2 H 3 * species adsorbed in blocked states, in which neighboring sites are occupied only by unreactive HO and Cl species. The microkinetic modeling shows that H-hopping away from surface HO groups provides a relatively efficient route for C 2 H 3 * to escape blocked configurations and dehydrogenate, and that this process can limit the C 2 H 4 selectivity on Cl-IrO 2 (110) under the conditions studied. Overall, our results demonstrate that Cl-substitution into IrO 2 (110) enhances the selectivity for C 2 H 4 production from C 2 H 6 and provides insights into the reaction mechanism that can guide strategies to further improve the C 2 H 4 selectivity.

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

Decoupling plasma, catalyst, and gaseous mechanisms for non-oxidative methane conversion

Direct non-oxidative methane (CH 4 ) conversion to value-added hydrogen (H 2 ) and C 2 products remains hindered by fundamental catalytic scaling constraints and rapid surface deactivation at elevated temperatures. Plasma-enabled catalysis offers a promising route to overcome the thermodynamic and kinetic barriers of direct non-oxidative methane upgrading at mild conditions, yet control over C–C product selectivity and catalyst stability remains elusive. Here, we establish a unified mechanistic framework including Langmuir–Hinshelwood (L–H) and Langmuir–Rideal (L–R) mechanisms that disentangles the roles of plasma excitation (including vibrationally activated species and radicals), surface temperature (T sur ), and catalyst binding energy in steering CH 4 conversion to H 2 and C 2 hydrocarbons. Through a combination of density functional theory (DFT) informed microkinetic modeling, in situ and ex situ surface characterization, and product quantification under dielectric barrier discharge conditions, we show that vibrationally excited CH 4 lowers activation barriers selectively for dissociative chemisorption, enabling surface activation across a wide range of transition metal catalysts at low thermal energy input. We find that once CH 4 is dissociatively chemisorbed, the branching between C 2 H 2 , C 2 H 4 , and C 2 H 6 is governed by surface properties (carbon binding energy, T sur , etc), regardless of plasma excitation. The DFT informed microkinetic model decouples the effects of molecular activation from surface properties and indentifies operating windows that maximize target yields while suppressing carbon accumulation and subsequent catalytic inactivation. Experiments on polycrystalline Cu/Al 2 O 3 , Ni/Al 2 O 3 , and Pt/Al 2 O 3 validate these predictions, revealing catalyst-dependent branching toward ethane or ethylene and distinct deactivation profiles. We unify these trends into a generalized three-dimensional plasma-thermal-catalytic design space, from which reduced descriptors such as T vib /T sur in the limit of vibrationally excited L–H pathways emerge as predictive metrics. These results enable rational tuning of methane conversion pathways and unlock selective C 2 formation using earth-abundant metals under mild plasma conditions.

catalyst inactivation↗
DOE OSTI2025

Nitrate Reduction Modeling under Acidic Conditions with Late Transition Metals

The electrochemical reduction of nitrate (NO 3 R) to ammonia is a bold yet conceivable way of producing ammonia using renewable electricity. However, serious challenges remain in finding optimal electrocatalysts for the process. An atomistic understanding of the surface energetics behind the NO 3 R is needed in order to design an efficient catalyst. Herein, we combine energetics from density functional theory and microkinetic modeling to demonstrate how surface descriptors can help simplify the search for efficient NO 3 R electrocatalysts. We illustrate the strong correlations between transition-state energetics and O* binding energies for adsorbed nitrate and nitrite on transition metals. For intermediates from NO* and beyond, we compare the benefits of using either the N* or H* binding energies to predict reduction onset potentials. These insights enable us to develop a simple microkinetic model that elucidates the surface coverages of intermediates and the product selectivity of NO 3 R across a range of potentials and transition metals. As a result, we show that the model adequately corroborates with quasi-steady-state rates observed experimentally.

ammonia↗
DOE OSTI2025

Theoretical and Experimental Insights into CO 2 Capture and Methanation over Amine-Grafted Ru-Based Catalysts

Carbon capture and storage (CCS) technologies, along with CO 2 capture and conversion methods, have emerged as crucial research areas to address rising CO 2 emissions. In this study, we seek to understand the mechanistic role of amines in enabling lower-energy pathways for CO 2 conversion. Our research focuses on the development and analysis of dual-functional materials (DFMs) engineered for the reactive capture and conversion (RCC) of CO 2 into methane, utilizing Ru catalysts grafted with amine groups. We employ Density Functional Theory (DFT) calculations using methylamine as a model amine to investigate the impact of amine groups on CO 2 methanation on a Ru(0001) surface, both in the presence and absence of amine groups. The amine ligand alters the carbon coordination environment, promoting direct C–O dissociation and potentially destabilizing the CO* adsorbate, thereby reducing the risk of CO poisoning. Additionally, we observe a preference for hydrogenation, although it becomes more energetically uphill in the amine-bound scenario. Our experiments, however, report similar CO 2 conversion and CH 4 production rates over the synthesized catalysts “Ru/TiO 2 ” and the amine (N-(2-aminoethyl)-3-aminoproplytrimethoxysilane (“diaminosilane”)) deposited catalyst “Diamine−Ru/TiO 2 ”. By constructing comparative reaction-free energy diagrams and performing microkinetic modeling (MKM) simulations, we link our theoretical findings with experimentally observed CO 2 uptake, conversion, and methane production rates. A microkinetic model was employed to investigate the anomaly, showing reduced amine–carbon complex coverage and increased CO 2 coverage at all temperatures. The MKM simulations consistently confirmed these trends. In conclusion, this comprehensive approach offers key insights into the role of the amine-CO 2 bond in methanation, highlighting a pathway toward lower-energy, more efficient CO 2 capture and conversion processes.

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

Multi‐Scale Model‐Informed Deep Learning for Plasma‐Nanoparticle Interaction

The Overarching Goal of this proposed research is to understand and quantitively determine the interactions between non-thermal plasma (hot electrons, reactive radicals, vibrationally excited species) and surface reactions on influencing the activity and selectivity of the desired reactions via developing multi-scale model informed deep learning algorithm. Investigating non-thermal plasma-surface interaction is feasible due to the low bulk temperature in the discharge region. To investigate the role of plasma-nanoparticle interaction on enhancing the reaction kinetics, we will focus on ammonia cracking to generate clean hydrogen over earth-abundant, non-critical metallic nanoparticles, which is of great significance for decarbonization. We hypothesize that (1) reactive radicals interacting with surface reaction species via Eley–Rideal mechanism will significantly lower the energetics of the potential rate-limiting step of nitrogen formation; (2) the surface will be charged heterogeneously under non-thermal plasma conditions and the charged site will lower the energetics of ammonia cracking through Langmuir– Hinshelwood mechanism; (3) vibrationally excited ammonia will further promote the initial N-H bond cleavage. To access the hypothesis, we will (1) reveal the surface charge effects on tunning the reaction energetics via interpretable, physics-informed deep learning accelerated density functional theory (DFT) calculations; (2) determine the reactive radicals interacting with surface reaction species on tuning the reaction energetics via DFT; (3) reveal the surface charge effects on tunning the reaction energetics via DFT and deep learning models, (4) quantify how vibrationally excited species, reactive radicals, and surface charging effects on enhancing the catalysis via developing DFT-based microkinetic modeling (MKM) and active learning. Deep and active learning of plasma-nanoparticle interactions effects on enhancing ammonia cracking to generate hydrogen represents a new paradigm for designing high performance plasma materials. The fundamental science of how plasma-nanoparticle interactions will change the plasma kinetics and will improve the energy efficiency for decarbonization and sustainability. The interpretable and physics-informed machine learning model will accelerate low temperature plasma chemistry and material discovery with physics rules and model interpretation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Group-Additivity–Embedded Multiscale Modeling for Electric Field-Enhanced Nanocatalysis

Elucidating structure-performance relationships remains a central challenge in field-enhanced catalysis, where nanoparticles exhibit nonuniform surface sites with site-dependent responses to electric fields. Low-coordination sites (edges, corners, and tips) are particularly electric field-sensitive (EF), leading to nonuniform charge distribution, adsorption energies, and catalytic activity. Here, using ammonia decomposition on a ruthenium cluster as a model system, we develop a transferable multiscale framework integrating density functional theory, group additivity (GA), Brønsted-Evans-Polanyi scaling, and microkinetic modeling to predict EF-dependent activity across nonuniform cluster sites. Across sites and fields, the nitrogen adsorption energy (E N ) emerges as the governing descriptor, yielding robust volcano relationships whose optimum shifts systematically with field: negative fields strengthen N binding via electron accumulation, while positive fields weaken N binding via charge depletion, moving the optimal E N toward weaker binding. Microkinetic analysis shows that N≡N bond formation remains the key kinetic bottleneck over most conditions; positive fields lower the effective barrier and, critically, increase the fraction of near-optimal active sites, leading to a net enhancement in overall activity relative to zero-field and negative-field cases. By capturing EF- and site-dependent energetics with high accuracy and low computational cost, this GA-embedded multi-scale simulation workflow provides a physically interpretable route to predict and design field-enhanced nanocatalysis.

ammonia decomposition↗
DOE OSTI2026

Interpretable Deep Learning for Advancing Field-Enhanced Catalysis

This DOE Early Career project developed a physics-informed, interpretable AI-and-modeling framework to understand and exploit electric-field effects in heterogeneous catalysis, with ammonia cracking and synthesis as a representative pathway. The team built and validated methods to map local electric fields on metal surfaces and nanoparticles, showing that low-coordination features (tips/edges/corners) can concentrate fields by several-fold relative to flat facets. Using DFT-generated datasets, the project created physics-guided machine learning models that rapidly predict local electric fields and field-dependent adsorption energetics with near-DFT accuracy while reducing computational cost by orders of magnitude. These predictions were integrated with microkinetic modeling to quantify how field-dipole interactions reshape reaction energetics and mechanisms, enabling large increases in predicted catalytic rates and substantial reductions in operating temperature under favorable field conditions. To accelerate discovery of earth-abundant catalysts, the project combined interpretable ML screening (with electronic-structure descriptors identified as key drivers) with a generative inverse-design workflow based on diffusion models and physics constraints. The resulting closed-loop approach, linking simulation, mechanistic modeling, and AI, provides reusable tools and datasets for designing catalysts and operating conditions in field-enhanced catalysis, with broad relevance to electrostatic catalysis, plasma catalysis, electrocatalysis, and other energy-related chemical transformations.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2025

Understanding of Ag Nanocatalysts for Electrocatalytic CO2 Conversion: Effects of Particle Size and Carbon Support

In this talk, we combined ultrahigh vacuum (UHV) surface science techniques, electrochemical measurements, and computational modeling to investigate Ag based electrocatalysts for CO2 reduction reaction (CO2RR). Our goal is to understand the critical characteristics governing the activity and selectivity of Ag electrocatalysts. Ag electrocatalysts were grown on highly oriented pyrolytic graphite (HOPG) in the UHV chamber, characterized with X-ray photoelectron spectroscopy (XPS) and scanning tunneling microscopy (STM), and then tested in a custom-built gastight H-cell. Supported by computational modeling based on density functional theory (DFT) calculations and microkinetic modeling (MKM), our studies revealed a strong size-dependent electrocatalytic CO2-to-CO conversion of the Ag nanoparticle electrocatalysts with average particle diameter between 2 to 6 nm. Smaller diameter (< 3 nm) particles favored H2 evolution reaction (HER) due to a high population of Ag edge sites, whereas larger diameter particles favored CO2RR as the population of Ag(100) surface sites grew. We further discovered that electronic interactions between small diameter Ag particles and highly defective carbon supports could break the size-dependent CO2RR reactivity, resulting in highly selective (CO Faradaic Efficiency > 90%) and active Ag nanoparticle electrocatalysts with sizes < 2 nm diameter. This knowledge is key to understand electrocatalysts performance and to ultimately guide electrocatalyst design

Ag nanoparticles↗
DOE OSTI2025

Mechanistic insights into nitrogen activation on atomic Ru clusters in self-pillared pentasil using operando atomistic models and experimental kinetics

Alternative catalysts to the industrial Haber Bosch process have been of significant interest in the field of heterogeneous catalysis, yet realizing ammonia synthesis under mild conditions (e.g., 300 °C and 10 bar) is challenging due to the low per-pass conversion. One strategy is to promote the associative ammonia synthesis mechanism which eschews direct N-N bond cleavage. Here, in this work, we use self-pillared pentasil, a self-pillared hierarchical zeolite built by thin MFI zeolite nanosheets, as a support for subnanometric Ru clusters to synthesize ammonia. We show that Ru remains well-dispersed during reaction and further demonstrate that ammonia synthesis rates are higher than Cs-Ru/MgO. Reaction kinetics show a positive order in H 2 providing evidence for the associative mechanism, which then becomes negative in H 2 if Ru is allowed to aggregate into nanoparticles. Operando Density Functional Theory models for Ru speciation in SPP, free energy diagrams, and microkinetic modeling were then applied to develop a reaction mechanism that involves sequential hydrogenation of N 2 from metallic Ru clusters. For this hydrogenation to occur, there are site requirements for N 2 to adopt a bridge-bound configuration that facilitates sequential hydrogenation on single sites and metal clusters. These site requirements in turn inform the design of improved zeolite-supported ammonia synthesis catalysts.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Cross-Scale Catalyst Modeling Applied to H 2 Storage and Release via Formic Acid

Here, we propose the Systems-to-Atoms (S2A) modeling framework that integrates the kinetics of reaction chemistry and structural configurations across various length scales with the aim of establishing a versatile template for multiscale modeling of reactive flow problems and to predict the operando activity of catalyst materials. The approach encompasses a microkinetic model to analyze surface reactions on individual facets of catalyst nanoparticles coupled with the computation of average surface reaction rates for catalyst nanoparticles of specific size distributions. Macro-homogeneous surface reaction kinetics are derived as a function of catalyst loading and used as input parameters for the continuum-scale reactor model. The cross-scale framework enables the optimization of catalyst utilization through reactor design and operating strategy. To demonstrate the framework, we studied the storage and release of hydrogen from formic acid, a promising liquid organic hydrogen carrier (LOHC), over Pd, Pt, and Cu catalysts. The framework predicts observed trends in formic acid dehydrogenation activity for catalysts with comparable weight loadings and metal particle diameters, demonstrating satisfactory quantitative alignment. Finally, the seamless transmission of parameter uncertainties between scales is also discussed.

08 HYDROGEN↗
DOE OSTI2025

Catalytic resonance theory for parametric uncertainty of programmable catalysis

Microkinetic models are useful tools for screening catalytic materials; however, errors in their input parameters can lead to significant uncertainty in model predictions of catalyst performance. Here, in this work, we investigate the impact of linear scaling and Brønsted-Evans-Polanyi relation parametric uncertainty on microkinetic predictions of programmable-catalyst performance. Two case studies are considered: a generic A-to-B prototype reaction and the oxygen evolution reaction (OER). The results show that error-unaware models can accurately predict trends and, for the prototype reaction, values of optimal waveform parameters. The specific model parameters driving output uncertainty are identified via variance-based global sensitivity analysis. However, predictions of dynamic rate enhancement can decrease when uncertainty is propagated into the models. In both cases, we identify operating conditions where the programmable catalyst achieves a rate enhancement of at least one order of magnitude despite parametric uncertainty in the model, supporting programmable catalysis as a viable strategy for exceeding the Sabatier limit.

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

Understanding the Unique Reactivity of Cu for Electrochemical CO 2 Reduction with a 3-Site Model

Cu-based catalysts for the electrochemical reduction of CO 2 and CO exhibit a perplexingly unique reactivity toward multicarbon based products compared to other studied electrocatalysts. Here, in this study, we use insights gained from a recent phenomenological 3-site microkinetic model and grand-canonical density functional theory calculations to clarify the importance of an underemphasized aspect critical to Cu’s unique reactivity: a population of so-called “reservoir” sites. Using model Cu surface motifs, we discuss how these types can be represented by undercoordinated structural defects like step edges and grain boundaries which form a network of highly anisotropic migration channels. These pathways are found to be amenable for feeding *CO over time to reactive sites like Cu adatoms more active toward C–C coupling. These results highlight an often overlooked aspect of catalyst optimization: reservoir site engineering, which exploits surface mobility and presents an equally important avenue for electrocatalyst engineering oversimply maximizing active site densities.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

First-Principles Insights into the Thermocatalytic Cracking of Ammonia-Hydrogen Blends on Fe(110). 2. Kinetics

Ammonia (NH 3 ) is an energy-rich molecule that is routinely synthesized from nitrogen (N 2 ) and hydrogen (H 2 ). NH 3 ’s more favorable physical properties compared to H 2 suggests it may offer a way to more conveniently store, transport, and, when needed, extract H 2 via thermal decomposition. However, the high kinetic barrier and endoergicity to decompose to H 2 and N 2 require high temperatures. The standard reaction free energy indicates nearly 100% thermodynamic conversion to the diatomic molecules only at ~673 K and higher. However, even at these temperatures, a catalyst, e.g., iron (Fe), is needed for favorable kinetic conversion. Here, in this study, we explore via density functional theory the kinetics of NH 3 decomposition on the most stable facet of body-centered cubic Fe, namely, (110), under typical high-temperature and finite-pressure operando conditions. We predict coverage-dependent energetics of elementary surface reactions, often neglected in atomic-scale modeling. From these models, we find the recombinative desorption of adsorbed N as N 2 is rate-determining at 573.15–773.15 K and even at an extreme case of 1173.15 K. From microkinetic modeling, we find that the steady-state turnover frequencies (TOFs) for N 2 and H 2 generation rates (r$_{H_2}$) depend exponentially on temperature. The catalyst achieves a steady-state TOF of 36.4 s –1 and an r$_{H_2}$ of 0.107 μmol cm –2 s –1 for a feed of 1.8 bar NH 3 with 0.2 bar H 2 at 1173.15 K. However, at 773.15 K, with the same feed composition and velocity, the steady-state TOF and r$_{H_2}$ decrease to 0.14 s –1 and 4.10 × 10 –4 μmol cm –2 s –1 , respectively, as the process is significantly hindered by slow N 2 desorption. Although at first glance counterintuitive, our simulations suggest that surface modifications that reduce Fe’s reactivity toward NH x species should enhance its overall NH 3 decomposition activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗