Abstract:The Haber–Bosch (HB) process is the primary chemical synthesis technique for industrial production of ammonia (NH3) for manufacturing nitrate-based fertilizer and as a potential hydrogen carrier. The HB process alone is responsible for over 2% of all global energy usage to produce more than 160 million tons of NH3 annually. Iron catalysts are utilized to accelerate the reaction, but high temperatures and pressures of atmospheric nitrogen gas (N2) and hydrogen gas (H2) are required. A great deal of research has aimed at increased performance over the last century, but the rate of progress has been slow. This Account focuses on determining the atomic-level reaction mechanism for HB synthesis of NH3 on the Fe catalysts used in industry and how to use this knowledge to suggest greatly improved catalysts via a novel paradigm of catalyst rational design.We determined the full reaction mechanism on the two most active surfaces for the HB process, Fe(111) and Fe(211)R. We used density functional theory (DFT) to predict the free-energy barriers for all 12 important reactions and the 34 most important 2 × 2 surface configurations. Then we incorporated the mechanism into kinetic Monte Carlo (kMC) simulations run for several hours of real time to predict turnover frequencies (TOFs). The predicted TOFs are within experimental error, indicating that the predicted barriers are within 0.04 eV of experiment.With this level of accuracy, we are poised to use DFT to improve the catalyst. Rather than forming bulk alloys with uniform concentration, we aimed at finding additives that strongly prefer near-surface sites so that minor amounts of the additive might lead to dramatic improvements. However, even for a single additive, the combinations of surface species and reactions multiplies significantly, with ∼48 reaction steps to examine and nearly 100 surface configurations per 2 × 2 site. To make it practical to examine tens of dopant candidates, we developed the hierarchical high-throughput catalysis screening (HHTCS) approach, which we applied to both the Fe(111) and Fe(211) surfaces. For HHTCS, we identified the most important 4 reaction steps out of 12 for the two surfaces to examine >50 dopant cases, where we required performance at each step no worse than for pure Fe. With HHTCS, the computational cost is about 1% of that for doing the full reaction mechanism, allowing us to do ≈50 cases in about 1/2 the time it took to do pure Fe(111). The new leads identified with HHTCS are then validated with full mechanistic studies.For Fe(111), we predict three high-performance dopants that strongly prefer the second layer: Co with a rate 8 times higher, Ni with a rate 16 times higher, and Si with a rate 43 times higher, at 400 °C and 20 atm. We also found four dopants that strongly prefer the top layer and improve performance: Pt or Rh 3 times faster and Pd or Cu 2 times faster. For Fe(211), the best dopant was found to be second-layer Co with a rate 3 times faster than that for the undoped surface.The DFT/kMC data were used to predict reshaping of the catalyst particles under reaction conditions and how to tune dopant content so as to maximize catalytic area and thus activity. Finally, we show how to validate our mechanistic modeling via a comparison between theoretical and experimental operando spectroscopic signatures.The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.accounts.1c00789.Schematic depiction of the bcc surface structures here investigated also showing 4-fold binding sites (PDF)This article has not yet been cited by other publications.

Reaction Mechanism and Kinetics for Ammonia Synthesis on the Fe(211) Reconstructed Surface.

Reaction Mechanism and Kinetics for Ammonia Synthesis on the Fe(111) Surface

Reaction Mechanisms, Kinetics, and Improved Catalysts for Ammonia Synthesis from Hierarchical High Throughput Catalyst Design

Highly Efficient Ni-Doped Iron Catalyst for Ammonia Synthesis from Quantum-Mechanics-Based Hierarchical High-Throughput Catalyst Screening

Operando probing of the surface chemistry during the Haber–Bosch process

QM-Mechanism-Based Hierarchical High-Throughput in Silico Screening Catalyst Design for Ammonia Synthesis.

Discovery of Dramatically Improved Ammonia Synthesis Catalysts through Hierarchical High-Throughput Catalyst Screening of the Fe(211) Surface

Correlation Between the Structural Features and Intrinsic Activity Trend of Fe Surfaces for Ammonia Synthesis

Si-Doped Fe Catalyst for Ammonia Synthesis at Dramatically Decreased Pressures and Temperatures

Ammonia Decomposition on Fe(1 1 0), Co(1 1 1) and Ni(1 1 1) Surfaces: A Density Functional Theory Study

Ammonia decomposition on Fe(1 1 0), Co(1 1 1) and Ni(1 1 1) surfaces: A density functional theory study

Theoretical Investigation of the Ammonia Adsorption Process on Si(100)-2×1 Surface

Adsorption Of Ammonia On Si(100), (111) Surfaces - A Theoretical Approach

Reaction Kinetics for Ammonia Synthesis using Ruthenium and Iron based Catalysts under Low Temperature and Pressure Conditions

Atmospheric Pressure Ammonia Synthesis on AuRu Catalysts Enabled by Plasmon-Controlled Hydrogenation and Nitrogen-species Desorption

Plasma‐catalytic Ammonia Synthesis via Eley‐Rideal Reactions: a Kinetic Analysis

Theoretical Kinetics Predictions for Reactions on the NH2O Potential Energy Surface

The synthesis mechanism of nitrogen to ammonia on the Fe, Co, Ni-doped Cu(1 0 0) Surface: A DFT study

Formation and decomposition of chemically activated and stabilized hydrazine

Dissociative Chemisorption Dynamics of Ammonia on Ni(111)

Mechanocatalytic Synthesis of Ammonia: State of the Catalyst During Reaction and Deactivation Pathway