The nearby atomic environment effect on an Fe–N–C catalyst for the oxygen reduction reaction: a density functional theory-based study

Abstract

Fe–N–C materials have emerged as highly promising non-noble metal catalysts for oxygen reduction reactions (ORRs) in polymer electrolyte membrane fuel cells. However, they still encounter several challenges that need to be addressed. One of these challenges is establishing an atomic environment near the Fe–N₄ site, which can significantly affect catalyst activity. To investigate this, herein, we employed density functional theory (DFT). According to our computational analysis of the Gibbs free energy of the reaction based on the computational hydrogen electrode (CHE) model, we successfully determined two C–O–C structures near the Fe–N₄ site (referred to as str-11) with the highest limiting potential (0.813 V). Specifically, in the case of O-doped structures, the neighboring eight carbon (C) atoms around nitrogen (N) can be categorized into two distinct types: four C atoms (type A) exhibiting high sensitivity to the limiting potential and the remaining four C atoms (type B) displaying inert behavior. Electronic structure analysis further elucidated that a structure will have strong activity if the valence band maximum (VBM) around its gamma point is mainly contributed by d_xz, d_yz or d_z² orbitals of Fe atoms. Constant-potential calculations showed that str-11 is suitable for the ORR under both acidic and alkaline conditions with a limiting potential of 0.695 V at pH = 1 and 0.926 V at pH = 14, respectively. Additionally, microkinetic simulations indicated the possibility of str-11 as the active site for the ORR under working potential at pH = 14.