Mapping the Binary Covalent Alloy Space to Pursue Superior Nitrogen Reduction Reaction Catalysts

Abstract The electrochemical nitrogen reduction reaction (NRR) has the potential to decarbonize industrial ammonia production. However, NRR has poor activity and selectivity versus the competing hydrogen evolution reaction for catalysts that adhere to scaling relations. Overcoming the limitations imposed by scaling relations requires more complex catalyst materials, however, evaluating materials beyond simple metal systems is a large combinatorial problem that requires an improved understanding of the electrocatalyst surface to rationally guide the discovery of superior catalysts. The study uses grand canonical density functional theory to uncover NRR trends on a large and disparate set of binary covalent alloys (BCA) with variable compositions and active‐site geometries. The studied BCAs generally follow scaling relations, albeit with larger variance and several systems that significantly break scaling. BCAs with early‐ to mid‐transition metals tend to lie near the volcano peak and activate the N 2 triple bond via a side‐on binding configuration. Trends in the BCA space cannot be readily predicted using simple electronic descriptors, which is ascribed to the large geometric variability of the BCA surfaces. It is anticipated that these findings will provide a foundation for the rational design of superior NRR electrocatalysts with increasing material complexity.