Platinum-group free Fe-N-C catalysts have been shown to be a promising class of electrocatalysts for the oxygen reduction reaction (ORR) in acidic media. However, typical synthesis methods can result in a variety of FeN4 sites embedded in graphene, including bulk-hosted and edge-hosted FeN4. The planar structure of the FeN4 active site allows two additional ligands to adsorb to Fe on either side of the graphene, which can be present from ORR reaction intermediates and other solution phase species. Here, we studied the impact of axially adsorbed ligands on bulk-hosted and edge-hosted FeN4 sites by screening over all pairs of H, O, OH, OOH, O2, H2O, NO, SO4, HSO4, and ClO4 axial ligands using solvated grand canonical density functional theory (GCDFT) to self-consistently account for applied potential. The first axial Fe adsorbate can strongly influence the binding of a second axial adsorbate and increase its adsorption energy by up to 2 eV, with more destabilization occurring for a more positive applied potential and the O* ORR intermediate being affected most. However, GCDFT calculations on these same sites without axial spectator ligands show that O* and OH* are both thermodynamically persistent, suggesting that axial spectator ligands might play a role in explaining this catalyst’s high experimental activity. Additionally, we show how the combined effects of potential, nature of the FeN4 site, and axial ligand impact O2 adsorption favorability and availability of sites in the presence of competing molecules. This work highlights the need for the continued development of advanced electrocatalyst models that account for a wide range of active site structures to facilitate comparison with experiment.
Read full abstract