Axial Coordination Effect on the Oxygen Reduction Reaction of FeN4 Electrocatalysts Based on Grand Canonical Density Functional Theory

The Fe–N–C materials are promising noble-free alternatives to Pt-based oxygen reduction reaction (ORR) electrocatalysts. Identifying the actual active structure of FeN4-based moieties is helpful for the regulation and design of high-performance electrocatalysts. However, the current theoretical researches were mostly based on the charge-neutral model (CNM), which cannot accurately describe the electrochemical interface. Herein, we employed the constant potential based grand canonical density functional theory (GC-DFT) computations to study the ORR mechanism of axially decorated FeN4 electrocatalysts. We studied around ten types of axial ligands, and the constant potential energetics and microkinetic modeling demonstrated that the Fe center can exhibit excellent activity for boosting the four-electron ORR via covalently linked −NH2 ligand. The lowering of the antibonding dz2–pz orbital is responsible for weakening the adsorption of oxygen in FeN4–Ls to promote ORR activity, and the −NH2 decoration led to the lowest antibonding orbital energy. In particular, the adsorption free energy (ΔG*O2) and O–O bond length (LO–O) of the adsorbed O2 reactant can be used as the effective energetic and geometric descriptors to describe the ORR activity of FeN4–Ls electrocatalysts. Our results not only elucidate the axial coordination effect on the ORR performance of FeN4 SACs but also demonstrate the importance of electrode potential in computational electrochemistry.