Search NASA⌕ Search

SEARCH · Search NASA

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Calculations of the pH-Dependent Onset Potential for CO Electrooxidation on Au(111)

CO electrooxidation on the Au(111) surface exhibits an onset potential that depends strongly on the pH of the electrolyte. In alkaline solution, the onset potential for CO electrooxidation is found, experimentally, to be 0.5 V lower than in acidic media on the reversible hydrogen electrode scale. This phenomenon is explained here with density functional theory which is used to calculate adsorbate binding energies including the electric double layer of the Au(111)/aqueous interface. Our model consists of a charged Au(111) slab and implicit solvation for the electrolyte. The double reference method is used to determine the potential-dependent CO electrooxidation reaction energetics. A microkinetic model, based upon the calculated reaction energetics, confi rms the lower onset potential in alkaline media. Our results show that there are three factors contributing to the lower onset potential in base: (1) stronger CO adsorption, (2) attraction between adsorbed CO and OH - , and (3) the high concentration of OH - in base.

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

On the nature of active sites for formic acid decomposition on gold catalysts

Formic acid decomposition has been studied experimentally on supported gold nanoparticles with strong evidence showing the critical role of gold clusters in the subnanometer range in catalyzing the reaction. However, there is a lack of theoretical studies capable of explaining these experimental observations. In this work, without accounting for support effects, vapor phase formic acid decomposition was studied systematically on sub-nanometric gold clusters from Au 4 to Au 25 , among which several candidate Au clusters were identified as a promising active site model for the Au/SiC catalysts. Combining theoretical and experimental results suggested that the active site on Au/SiC catalysts could be represented by an Au 18 cluster, on which the reaction rates calculated from the microkinetic model closely match the experimentally measured rates. On Au 18 , formic acid decomposition proceeds through the same formate mediated pathway as that on extended Au surfaces (HCOOH → HCOO + H → CO 2 + 2H → CO 2 + H 2 ), with the reaction taking place on a triangular Au 3 site where the reactive Au atoms have a coordination number of 5. Despite the fact that other Au clusters among those studied, including Au17 and Au19 which also expose the same triangular Au 3 site, they were not found to be an accurate representation of the active sites, suggesting an atom-specific activity of gold clusters for formic acid decomposition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2025

Multi-Scale Modeling for Plasma-Enhanced Ammonia Decomposition over Carbides and Nitrides

Ammonia is a carbon-free hydrogen carrier, but its decomposition typically requires high temperatures over costly Ru-based catalysts due to the large barrier for N≡N bond formation. We develop a multiscale framework combining density functional theory, zero-dimensional plasma kinetics, and microkinetic modeling to elucidate how non-thermal plasma (NTP) enables low-temperature NH 3 decomposition over Co-based carbides and nitrides, benchmarked against Ru and Co. Under thermal conditions, all catalysts are limited by N≡N bond formation, with Co 3 C(001) most active owing to its negatively charged surface, strong N* binding, and low activation barriers of N≡N bond formation. Plasma-induced vibrational excitation of NH 3 and its reactive radicals promotes a radical-driven •NH 2 –N* coupling pathway that dominates on Co 3 C(001) and Co 3 N(001), shifting the rate-limiting step to NH 3 (v1) dissociation, increasing turnover frequencies by up to 6 orders of magnitude, and reducing the temperature needed to reach a turnover frequency of 5 s –1 from >680 °C (Ru and Co under thermal condition) to 267 °C (Co 3 C) and 415 °C (Co 3 N). These results identify Co-based carbides and nitrides as promising plasma-active catalysts for energy-efficient hydrogen production from ammonia.

ammonia decomposition↗
DOE OSTI2020

Oxidative dehydrogenation of propane on the oxygen adsorbed edges of boron nitride nanoribbons

Metal-free boron nitride has recently been reported to exhibit high catalytic activity and selectivity towards oxidative dehydrogenation of propane (ODHP). Despite several experimental studies, the exact nature and function of the active sites in this emerging catalyst are unknown. In the present work, density functional theory calculations are combined with microkinetic modeling to systematically explore the ODHP on oxygen passivated boron nitride nanoribbons (BNNRs). Here, the relative stabilities of different edge structures of BNNR were examined using ab initio thermodynamic analysis, and the most stable oxygen passivated edge structure was used for a mechanistic study. Microkinetic analysis revealed that a N 2 O or NO x -type active site is active and selective for the conversion of propane to propene. In addition to a heterogeneous catalytic cycle, the proposed N 2 O/NO x -type active sites are expected to be able to generate gas phase C 3 H 7 ˙, which in turn can trigger gas phase reactions of ODHP as experimentally speculated.

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

Modeling Hydrogen Evolution Reaction Kinetics through Explicit Water–Metal Interfaces

Despite the apparent simplicity of the hydrogen evolution reaction (HER) and the decades of research into it, controversy remains in the literature regarding the identity of the active site and the competition between the Heyrovsky and Tafel steps. In this work, we use charge-extrapolated ab initio simulations with explicit water in conjunction with mean-field microkinetic modeling to explore the mechanism for HER on both close-packed (111) and stepped (211) transition metals. First, we show that atop H*, beyond a monolayer of hollow H*, is unlikely to play a role in the HER mechanism, given its very positive adsorption energies. The energetics suggests the Volmer–Heyrovsky mechanism to predominate on fcc transition metals under typical operating conditions. Here, we evaluate our theoretical results vs several experimental observations. We show that the Volmer–Heyrovsky mechanism predicts an activity volcano with its peak at a H* binding ΔG H* ≈ 0 eV, consistent with experiment. In contrast, the Volmer–Tafel volcano shows a broad rate plateau between ΔG H* ≈ 0 eV and ΔG H* ≈ – 0.4 eV. We find our theoretical Tafel slopes to be consistent with experimental ones on a range of transition metals. We show that, in line with experimental observations, the introduction of a CO(g) atmosphere shifts the strong binding metals toward the weak binding leg. Our study suggests that the simple thermodynamic approach to HER activity still holds, even when a detailed kinetic picture is considered.

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

Multiscale modeling of hydrogenolysis of ethane and propane on Ru(0001): Implications for plastics recycling

Plastic waste presents an environmental threat. Chemical recycling via hydrogenolysis can convert plastic waste into waxes, lubricants, and fuels. Among catalysts, Ru stands out for its superior activity and selectivity. The chemistry of light alkane hydrogenolysis can help understanding plastics deconstruction. Here, we perform first-principles calculations, develop descriptor-based relations, and conduct microkinetic modeling and analysis on ethane and propane. Predictions are in excellent agreement with experimental data. We identify a similar cracking pattern for both hydrocarbons entailing a deeply dehydrogenated species with the removal of four hydrogen atoms: CHCH*+* → 2CH* for ethane and CH 3 CCH*+2* → CH 3 C* + CH* for propane. We find that the rate-determining step is the C-C cracking for ethane and the first dehydrogenation from the terminal carbon (CH 3 CH 2 CH 3 *+* → CH 3 CH 2 CH 2 *+H*) for propane. The vinyl species CH 2 CH* produced from propane cracking is responsible for whether a single or multiple cracking events occur and effectively controls the selectivity. Specifically, ethane formation in propane hydrogenolysis is suppressed at elevated temperatures due to over-cracking via multiple (two here) cracking events being preferred over hydrogenation and desorption of ethane from the catalyst. Our workflow and models provide a baseline for future studies on heavier hydrocarbons. Insights into recent experimental studies of polyethylene over Ru-based catalysts are discussed.

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

Institute for Catalysis in Energy Processes (ICEP) (Final Report)

The Institute for Catalysis in Energy Processes (ICEP) was a multi-PI program located at the Northwestern University Center for Catalysis and Surface Science from 2006-2024. Over several renewal cycles and several organizing themes, ICEP addressed fundamental questions in catalysis science. In turn, the scientific questions addressed were directly relevant to the efficient use of the nation’s resources to produce fuels and commodity chemicals, harness alternate energy sources for chemical reactions such as light or electricity, reduce emissions and waste, and minimize the impact of our use of plastics. Selective oxidation, (oxidative) dehydrogenation, deNOx, CO 2 reduction, hydrogen release, polymer decomposition / depolymerization, and many other reactions of intense interest to the US Department of Energy were studied in the center. Major ICEP strengths were in catalyst synthesis, in measurements that elucidated catalyst properties, in reaction mechanisms and kinetics, and in predictive theory and modeling. ICEP researchers were often the inventors or developers of materials and methods that were brought to bear on the Institute’s catalytic systems. Key characterization tools advanced through this project included resonance Raman and related techniques, sum frequency generation, and X-ray standing wave spectroscopy. Center members relied extensively on computational tools such as density functional theory and microkinetic modeling to fully investigate catalyst structures and catalytic mechanisms. ICEP researchers were early developers of atomic layer deposition (ALD) for catalyst synthesis and were early pioneers of metal organic frameworks for catalysis. ICEP innovations are found in areas of intense research such as single-atom catalysts, catalytic deconstruction of plastics, and catalytic MOFs, to name a few. Research areas initiated by ICEP indirectly led to DOE EFRCs, startup companies, and other achievements. Over its 18 years of existence, ICEP involved two different PIs (Stair and Notestein) and 24 senior investigators across several academic disciplines at Northwestern University. There were tight collaborations with Argonne National Laboratory, including specialized equipment and experiments permanently located there. The center supported ~400 person-years of effort by graduate students and postdocs, either directly through DOE support or indirectly by leveraging independent support provided to the trainees, e.g. through the NSF GRFP or internal fellowships. The project resulted in 328 manuscripts and numerous conference presentations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Hydrogen Dissociation Controls 1-Hexyne Selective Hydrogenation on Dilute Pd-in-Au Catalysts

Increasing the selectivity of the catalytic hydrogenation of alkynes to alkenes is of major importance for the processing of petrochemicals and the production of fine chemicals. Achieving high selectivity for alkene formation at high conversions, however, remains a long-standing challenge in heterogeneous catalysis. Here, the mechanism and origin of the high selectivity of dilute Pd-in-Au catalysts has been studied by a combination of first-principles calculations, microkinetic simulations, and isotopic exchange hydrogenation experiments using Pd 0.04 Au 0.96 nanoparticles embedded in raspberry colloid-templated silica. The Pd is predominantly in the form of isolated atoms, only surrounded by Au atoms, based on prior studies. The simulations indicate that the rate-limiting process for 1-hexyne hydrogenation on Pd monomers in Au(111) is H 2 dissociation, which has a large free energy barrier of 0.86 eV at 363 K and 0.2 bar of H 2 . The C–H bond formation steps, on the other hand, proceed with lower barriers, which contrasts with previous studies of extended Pd catalysts. The microkinetic simulations identify the sizable H 2 dissociation barrier and the small barrier for the hydrogenation of 1-hexyne as key factors that lead to a high selectivity for the production of 1-hexene from 1-hexyne, even at high conversion. The unconventional H 2 dissociation limiting process in combination with the low coverage of weakly bound hydrocarbon intermediates explains the near-zero order of 1-hexyne found experimentally. Furthermore, the partial hydrogenation of 1-hexyne to form 1-hexene is shown to be an irreversible process from our isotopic exchange hydrogenation experiments and is explained by the strongly exothermic nature of the reaction. Diluting active species, like Pd, in a less active host metal, like Au, hence appears promising as a means of tuning the binding energy of reactants and altering reaction profiles, leading to distinct kinetic behavior for an optimal catalytic activity and selectivity. Furthermore, the combination of microkinetic modeling, density functional theory calculations, and isotopic exchange experiments is thus demonstrated to be an effective approach to modeling important catalytic phenomena.

1-hexyne↗
DOE OSTI2023

Operational and scaling-up barriers of SOEC and mitigation strategies to boost H 2 production- a comprehensive review

Solid oxide electrolysis cells (SOECs) proffer a flexible and environmentally friendly solution to reduce CO 2 emissions and/or produce H 2 fuel for various applications. Here, the integration of this technology with renewable energy systems makes them even more promising in realizing a near-perpetual green energy infrastructure in the near future. SOEC technology has achieved stupendous improvement in recent years, and some of the limiting factors to the commercialization of this technology have been explored. It is, therefore, necessary to have a succinct outlook of the different factors that contribute to the recent developments in SOEC technology. The current work presents a comprehensive review of the state-of-the-art SOEC component materials and their degradation mechanisms, SOEC modeling essentials and efficiencies, and the scalability of SOECs. Although SOEC is considered a highly promising technology and sufficient multifaceted research has been done, yet not been widely commercialized. Despite the improvement to mitigate degradation of the SOEC materials, a large share of the degradation concerns sprang from charge transfer reaction at the triple-phase boundaries of the H 2 -electrode. Part of which can be improved by microkinetic modeling of these materials supported by their characterization to give insight into their development. Some other directions for the future in advancing this technology toward the commercialization stage are thoroughly addressed.

08 HYDROGEN↗
DOE OSTI2022

Modified Energy Span Analysis of Catalytic Parallel Pathways and Selectivity

Mechanistic modeling provides vital insights into catalytic reactions. To analyze complex reaction networks with parallel pathways, we leverage the graph theory approach of the Energy Span Model (ESM) to develop a modified energy span analysis (MESA). A new method of cycle plots is proposed to perform reaction pathways analysis visually. We demonstrate this method on two published models: one describing carbon monoxide oxidation and the other simulating ethylene conversion to propanal via hydroformylation or ethane via hydrogenation. Fundamental insights explain kinetic observables, such as a reactant’s negative reaction order. General principles are revealed, such as rate-determining surface species being outside the primary reaction flux cycle and pathway selectivity being a purely kinetic property when reaction conditions are not near equilibrium. Lastly, we demonstrate MESA’s consistency with published microkinetic modeling results, highlighting this technique’s extension of the ESM to heterogeneous catalysts using collision theory to describe adsorption steps and concentration effects.

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

Unifying principles for catalytic hydrotreating processes (Final Technical Report)

This project builds on the hypothesis that the hydrotreating processes for the removal of oxygen and sulfur are fundamentally similar at the atomic-scale and existing knowledge from the treatment of petroleum derived feedstock can be leveraged for the design of novel catalysts for the upgrade of bio-oil. We tested this hypothesis by comparing computed potential energy diagrams for hydrodesulfurization (HDS) of thiophene over MoS 2 with hydrodeoxygenation (HDO) of furan over MoO 3 and concluded that certain aspects, such as catalyst promotion with transition metals, are valid strategies for both reactions. On the other hand, we also noticed significant differences in the mechanism for hydrogen (H 2 ) activation, which requires sites with metallic character. While MoS 2 is known to have metallic edge states that can catalyze H 2 dissociation, this elementary step is prohibitively slow on defect-free oxides. Only in the presence of vacancies or by creating metal/oxide interfaces can efficient H 2 activation sites during HDO be formed. The need for bifunctional catalyst when it comes to efficient and selective HDO or dehydrogenation reactions was further corroborated in joint experimental and theoretical studies of the Guerbet reaction for the coupling of biomass derived oxygenates over PdCu alloys, nitrate reduction over In-promoted Pd nanoparticles, and ethylene dehydroaromatization over Ga-exchanged ZSM-5 zeolites. All of these catalytic systems have in common that catalytic sites with distinct functional requirements are needed to create a working catalyst. Detailed computational studies were carried out for HDO of m -cresol and phenol on Ru-modified TiO 2 surfaces, which allowed us to attribute catalytic activity to the metal/oxide interface. A surprising finding was that heterolytic cleavage of the H-H bond across the Ru/TiO 2 interface was critically important, despite lower barriers for homolytic H 2 activation on Ru metal. The explanation lies in the high barriers for hydrogen spillover from Ru to TiO 2 , which becomes unnecessary in the heterolytic activation pathway. Moreover, we also reported that proton transfer steps between metal and oxide sites are mediated by weakly adsorbed surface water. During attempts to develop and validate a kinetic Monte Carlo (kMC) model for HDO reactions at the Ru/TiO 2 interface, it became clear that lateral interactions are paramount to describe realistic surface chemistry and without these interactions, the reduction and hydroxylation behavior from our simulations was inconsistent with reported experiments. To assess the importance of lateral interactions in popular computational catalyst design strategies relying on the identification of reactivity descriptors, which can be used along with Brønsted–Evans–Polanyi (BEP) and scaling relations as input to a microkinetic model (MKM) to make predictions for activity or selectivity trends, we compared predicted trends with those obtained from descriptor-based kMC models. We critically evaluated the benefits of kMC over MKM in terms of trend predictions and computational cost when using only a small set of input parameters. After confirming that in the absence of lateral interactions the kMC and MKM approaches yield identical trends and mechanistic information, we observed substantial differences between the two kinetic models when lateral interactions were introduced. The mean-field implementation applies coverage corrections directly to the descriptors, causing an artificial overprediction of the activity of strongly binding metals. In contrast, the cluster expansion in kMC implementation can differentiate among the highly active metals but it is very sensitive to the set of included interaction parameters. Considering that computational screening relies on a minimal set of descriptors, for which MKM makes reasonable trend predictions at a ca. three orders of magnitude lower computational cost than kMC, we concluded that the MKM approach does provide an overall better entry point for computational catalyst design. Overall, this project has led to 11 peer-reviewed publications, and their scientific is impact is well illustrated by their combined 702 citations.

08 HYDROGEN↗
DOE OSTI2019

Barium-Promoted Ruthenium Catalysts on Yittria-Stabilized Zirconia Supports for Ammonia Synthesis

The cost-effective, small-scale, distributed synthesis of ammonia depends on effective catalysts and processes that operate under modest elevated pressure (i.e., p < 20 bar) conditions. The present paper considers Ru as the active catalyst supported on yittria-stabilized zirconia (YSZ). The addition of alkali and alkaline-earth metal promoters is found to increase synthesis rates by an order of magnitude. The rate enhancement is largely insensitive to the promoter concentration, with Cs outperforming Ba and K by a factor of 2. However, Ba is found to be stable whereas Cs degrades more rapidly, which is attributed to the low melting point of its oxide. At 400 °C and 1.0 MPa, the specific synthesis rate over Ba–Ru is measured to be approximately 1410 mmol gRu –1 h –1 , higher than the most active oxide-supported Ru catalysts reported in the literature. The rate becomes inhibited by H 2 absorption at low temperature (below 350 °C), but lower H 2 /N 2 ratios enable the rate to remain comparable to what is observed in stoichiometric mixtures at temperatures below 400 °C. The paper reports a new detailed microkinetic model that accurately captures the observed behavior, revealing that adsorption is coverage dependent. Furthermore, these results provide insight and direction into developing alternatives to Haber–Bosch for distributed synthesis of green ammonia.

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

Selective Benzyl Alcohol Oxidation over Pd Catalysts

In the last decades, the selective liquid phase oxidation of alcohols to the corresponding carbonyl compounds has been a subject of growing interest. Research has focused on green methods that use “clean” oxidants such as O 2 in combination with supported metal nanoparticles as the catalyst. Among the alcohols, benzyl alcohol is one of the most studied substrates. Indeed, benzyl alcohol can be converted to benzaldehyde, largely for use in the pharmaceutical and agricultural industries. This conversion serves as model reaction in testing new potential catalysts, that can then be applied to other systems. Pd based catalysts have been extensively studied as active catalytic metals for alcohol oxidation for their high activity and selectivity to the corresponding aldehyde. Several catalytic materials obtained by careful control of the morphology of Pd nanoparticles, (including bimetallic systems) and by tuning the support properties have been developed. Moreover, reaction conditions, including solvent, temperature, pressure and alcohol concentration have been investigated to tune the selectivity to the desired products. Different reaction mechanisms and microkinetic models have been proposed. The aim of this review is to provide a critical description of the recent advances on Pd catalyzed benzyl alcohol oxidation.

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

Theoretical Modeling of Electrochemical Proton-Coupled Electron Transfer

Proton-coupled electron transfer (PCET) plays an essential role in a wide range of electrocatalytic processes. Here, a vast array of theoretical and computational methods have been developed to study electrochemical PCET. These methods can be used to calculate redox potentials and pKa values for molecular electrocatalysts, proton-coupled redox potentials and bond dissociation free energies for PCET at metal and semiconductor interfaces, and reorganization energies associated with electrochemical PCET. Periodic density functional theory can also be used to compute PCET activation energies and perform molecular dynamics simulations of electrochemical interfaces. Various approaches for maintaining a constant electrode potential in electronic structure calculations and modeling complex interactions in the electric double layer (EDL) have been developed. Theoretical formulations for both homogeneous and heterogeneous electrochemical PCET spanning the adiabatic, nonadiabatic, and solvent-controlled regimes have been developed and provide analytical expressions for the rate constants and current densities as functions of applied potential. The quantum mechanical treatment of the proton and inclusion of excited vibronic states have been shown to be critical for describing experimental data, such as Tafel slopes and potential dependent kinetic isotope effects. The calculated rate constants can be used as input to microkinetic models and voltammogram simulations to elucidate complex electrocatalytic processes.

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

Toward Carbon Monoxide Methanation at Mild Conditions on Dual-Site Catalysts

The catalytic carbon monoxide (CO) methanation is an ideal model reaction for the fundamental understanding of catalysis on the gas–solid interface and is crucial for various industrial processes. However, the harsh operating conditions make the reaction unsustainable, and the limitations set by the scaling relations between the dissociation energy barrier and dissociative binding energy of CO further increase the difficulty in designing high-performance methanation catalysts operating under milder conditions. In this work, we proposed a theoretical strategy to circumvent the limitations elegantly and achieve both facile CO dissociation and C/O hydrogenation on the catalyst containing a confined dual site. The DFT-based microkinetic modeling (MKM) reveals that the designed Co-Cr 2 /G dual-site catalyst could provide 4–6 orders of magnitude higher turnover frequency for CH 4 production than the cobalt step sites. We believe that the proposed strategy in the current work will provide essential guidance for designing state-of-the-art methanation catalysts under mild conditions.

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

How the Bioinspired Fe 2 Mo 6 S 8 Chevrel Breaks Electrocatalytic Nitrogen Reduction Scaling Relations

The nitrogen reduction reaction (NRR) is a renewable alternative to the energy- and CO 2 -intensive Haber–Bosch NH 3 synthesis process but is severely limited by the low activity and selectivity of studied electrocatalysts. The Chevrel phase Fe 2 Mo 6 S 8 has a surface Fe–S–Mo coordination environment that mimics the nitrogenase FeMo-cofactor and was recently shown to provide state-of-the-art activity and selectivity for NRR. Here, we elucidate the previously unknown NRR mechanism on Fe 2 Mo 6 S 8 via grand-canonical density functional theory (GC-DFT) that realistically models solvated and biased surfaces. Fe sites of Fe 2 Mo 6 S 8 selectively stabilize the key *NNH intermediate via a narrow band of free-atom-like surface d-states that selectively hybridize with p-states of *NNH, which results in Fe sites breaking NRR scaling relationships. These sharp d-states arise from an Fe–S bond dissociation during N 2 adsorption that mimics the mechanism of the nitrogenase FeMo-cofactor. Furthermore, we developed a new GC-DFT-based approach for calculating transition states as a function of bias (GC-NEB) and applied it to produce a microkinetic model for NRR at Fe 2 Mo 6 S 8 that predicts high activity and selectivity, in close agreement with experiments. Furthermore, our results suggest new design principles that may identify effective NRR electrocatalysts that minimize the barriers for *N 2 protonation and *NH 3 desorption and that may be broadly applied to the rational discovery of stable, multinary electrocatalysts for other reactions where narrow bands of surface d-states can be tuned to selectively stabilize key reaction intermediates and guide selectivity toward a target product. Furthermore, our results highlight the importance of using GC-DFT and GC-NEB to accurately model electrocatalytic reactions.

10 SYNTHETIC FUELS↗