Periodic Defect Boundary-Mediated Activity of Electrocatalytic Oxygen Reduction Reactions of Fe-N-C Catalysts

Ziwei Wang,Huiting Niu,Tiantian Wu,Shujiang Ding,Bao Yu Xia,Yaqiong Su

DOI: https://doi.org/10.31635/renewables.024.202400058

2024-01-01

Abstract:Open AccessRenewablesRESEARCH ARTICLES22 Jun 2024Periodic Defect Boundary-Mediated Activity of Electrocatalytic Oxygen Reduction Reactions of Fe-N-C Catalysts Ziwei Wang†, Huiting Niu†, Tiantian Wu, Shujiang Ding, Bao Yu Xia and Yaqiong Su Ziwei Wang† School of Chemistry, Engineering Research Center of Energy Storage Materials and Devices of Ministry of Education, National Innovation Platform (Center) for Industry-Education Integration of Energy Storage Technology, Xi'an Jiaotong University, Xi'an 710049 , Huiting Niu† School of Chemistry and Chemical Engineering, Key Laboratory of Material Chemistry for Energy Conversion and Storage (Ministry of Education), Hubei Key Laboratory of Material Chemistry and Service Failure, Hubei Engineering Research Center for Biomaterials and Medical Protective Materials, Wuhan National Laboratory for Optoelectronics, Huazhong University of Science and Technology (HUST), Wuhan 430074 , Tiantian Wu *Corresponding authors: E-mail Address: [email protected] E-mail Address: [email protected] School of Chemistry, Engineering Research Center of Energy Storage Materials and Devices of Ministry of Education, National Innovation Platform (Center) for Industry-Education Integration of Energy Storage Technology, Xi'an Jiaotong University, Xi'an 710049 , Shujiang Ding School of Chemistry, Engineering Research Center of Energy Storage Materials and Devices of Ministry of Education, National Innovation Platform (Center) for Industry-Education Integration of Energy Storage Technology, Xi'an Jiaotong University, Xi'an 710049 , Bao Yu Xia *Corresponding authors: E-mail Address: [email protected] E-mail Address: [email protected] School of Chemistry and Chemical Engineering, Key Laboratory of Material Chemistry for Energy Conversion and Storage (Ministry of Education), Hubei Key Laboratory of Material Chemistry and Service Failure, Hubei Engineering Research Center for Biomaterials and Medical Protective Materials, Wuhan National Laboratory for Optoelectronics, Huazhong University of Science and Technology (HUST), Wuhan 430074 and Yaqiong Su *Corresponding authors: E-mail Address: [email protected] E-mail Address: [email protected] School of Chemistry, Engineering Research Center of Energy Storage Materials and Devices of Ministry of Education, National Innovation Platform (Center) for Industry-Education Integration of Energy Storage Technology, Xi'an Jiaotong University, Xi'an 710049 https://doi.org/10.31635/rnwb.024.202400058 SectionsSupplemental MaterialAboutAbstractPDF ToolsAdd to favoritesDownload CitationsTrack Citations ShareFacebookTwitterLinked InEmail The introduction of defects can adjust the activity of graphene-based single-atom catalysts for oxygen reduction reactions (ORR). Herein, we for the first time investigate the ORR catalytic activity of FeN4 sites embedded on graphene with four types of line-defective boundary via density functional theory calculations. Our results show that periodic line defects consisting of pentagon-pentagon-octagon (C585) or quad-octagon chains (C484) can significantly enhance ORR activity, owing to the optimized electronic structures of FeN4 sites. The spin magnetic moment and the valence state of the Fe atom are both well correlated with the ORR overpotential. Experimental investigations further corroborate that FeN4 with a high degree of defects exhibits better ORR activity and stability compared to FeN4 sites of pristine graphene and commercial Pt/C. This work unravels the influence of the periodic defect boundary on the ORR performance of Fe-N-C catalysts and paves the way towards the rational design of highly effective single-atom electrocatalysts. Download figure Download PowerPoint Introduction Transition metal nitrogen-graphene materials (M-N-C) have emerged as viable candidates among nonprecious metal electrocatalysts for the oxygen reduction reaction (ORR) due to their ability to easily accommodate foreign electrons and impair the bonding strength of oxygen-containing intermediates, particularly with MNxCy catalysts possessing 3d unoccupied orbits.1–5 Among them, single-atom Fe-N-C catalysts are considered one of the most promising materials due to their maximum atomic utilization and impressive ORR activity.6–10 Both theoretical and experimental studies have indicated that the FeN4 site is the source of activity in Fe-N-C catalysts.11–14 However, the ORR activity of FeN4 sites greatly depends on the carbon supports.15 Previous studies have demonstrated that the ORR activity on the FeN4 site can be enhanced by manipulating the electric structure of the carbon supports.16,17 It is widely recognized that introducing heteroatoms can regulate the charge distribution of FeN4 and significantly improve its ORR activity.16–21 Alternatively, topological defects in carbon materials can also modulate their electronic structures to improve the ORR performance.22 Hou et al.23 investigated the electrical characteristics of native defects in graphene and their impact on N doping, revealing that preintroducing the defect before introducing the N dopant enhances the N doping of graphene. Liu et al.15 found that divacancy defects on graphene modulated the ORR activity of Fe-N4 sites, demonstrating that the hybridization between Fe 3dz2,3dyz (3dxz) and O2 π* orbitals is the origin of Fe-N4 ORR activity. Zhang et al.24 studied the ORR mechanisms on the nitrogen-doped graphene in acidic environments and proposed that the Stone–Wales defect facilitates the ORR on the doped graphene. However, these studies all concentrate on the point defects or their irregular clusters of graphene. The influence of periodic defect boundaries on the ORR activity of the FeN4 site still needs to be investigated. Herein, we used density functional theory (DFT) calculations to systematically explore the ORR of the FeN4 sites embedded into graphene with four types of periodic line-defect boundary. The results demonstrate that a periodic defect boundary consisting of pentagon-pentagon-octagon (C585) or quad-octagon chains (C484) evidently enhance ORR activity of FeN4 sites, and C585(3h) shows its highest efficiency in catalyzing the ORR through the 4e− pathway, with the overpotential at ηORR = 0.61 V, much lower than that of FeN4 sites (0.86 V) on perfect graphene. The correlation between the ORR catalytic activity and the spin magnetic moment of Fe, as well as its valence state, was further identified. Additionally, the density of states (DOS), highest occupied molecular orbital (HOMO)–lowest unoccupied molecular orbital (LUMO) energy gap, and charge difference were calculated to elucidate the change origin of the spin magnetic moment (MFe) and Fe valence state considered as descriptors for ORR activity. Further experimental studies verify that the FeN4 catalysts with high-density defects can significantly improve the ORR performance compared to pristine FeN4 catalysts and commercial Pt/C catalysts. Results and Discussion Based on previous research, four types of graphene clusters with periodic defect boundaries were built (as illustrated in Figure 1a–d).25–28 Referring to Supporting Information Figure S1, it is clear that, compared to pristine graphene, high-density C585, low-density C585, C57, and C484 exhibit heightened DOS values near the Fermi level, indicative of their boosted metallic conductivity. In order to add a single Fe atom to the structure, two carbon atoms next to it were taken out, and four nearby carbon atoms were substituted with nitrogen atoms. Figure 1e shows the structure of different carbon substrates embedded in FeN4. Given the carbon ring's structure, symmetry, and positioning, Supporting Information Figure S2a illustrates four methods for eliminating two neighboring carbon atoms from C585. Supporting Information Figure S2b,c demonstrates the removal of two adjacent carbon atoms from C57 employing three methods. A similar process is depicted for C484 using two methods. As a result, a total of 13 FeN4 structures were constructed, comprising four types on the high-density C585 ( Supporting Information Figure S3a), four types on the low-density C585 ( Supporting Information Figure S3b), three types on the C57 ( Supporting Information Figure S4), and two types on the C484 ( Supporting Information Figure S5). We conducted a comprehensive analysis of the formation energy of FeN4 systems (Figure 1f). The benchmark for our study was set at the formation energy of the pristine graphene FeN4 (Fe-N@C) configuration (2.31 eV), represented by the red dashed line. Our findings indicate that the formation energy of FeN4 structures range from −0.19 to 2.82 eV, with most periodic defected structures exhibiting more negative formation energies than the Fe-N@C configuration. And the binding energies (Eb) of all the systems were negative, ranging from −0.25 to −2.72 eV ( Supporting Information Figure S6). This suggests that periodic defected graphene, particularly high-density Fe-N@C585 and Fe-N@C484 structures, facilitated the incorporation of Fe and N atoms into the carbon matrix compared to Fe-N@C. Furthermore, we conducted ab initio molecular dynamics simulations to verify the dynamic stability of C585(3h), as illustrated in Supporting Information Figure S7. Figure 1 | Periodic defective graphene clusters (a–d) and (e) the corresponding FeN4 systems. Brown sphere: carbon; yellow sphere: iron; gray sphere: nitrogen. (f) The formation energies Ef of Fe-N4 systems. (g) The calculated overpotentials for ORR on Fe-N4 site supported by defected graphene. The red dashed line indicates the overpotential of FeN4@G (0.86 V). Download figure Download PowerPoint The catalytic performance of FeN4 structures for ORR is discussed in this section. Figure 1g illustrates the ORR overpotentials of Fe-N4 sites. We observe that Fe-N@C585 structures and Fe-N@C484 structures exhibit lower overpotentials compared to Fe-N@C (0.86 eV), with C585(3h) exhibiting its highest efficiency in catalyzing the ORR, achieving optimal activity with an overpotential of 0.61 V. However, Fe-N@C57 displays higher overpotentials, indicating inferior ORR activity. The negligible difference in overpotential between high and low defect densities of Fe-N@C585 structures suggests that the density of periodic defect boundaries on the carbon support has minimal impact on the ORR activity of FeN4. The ORR free energy diagrams of 13 models at U = 1.23 V are shown in Figure 2a,b and Supporting Information Figure S8a,b, in which the step with the maximum ΔG (ΔGmax) represents the potential-determining step (PDS). We found that the last step (OH* → H2O) is the PDS at all FeN4 configurations, and the periodic defect regulates the free energy difference of this step on the Fe-N4 site ( Supporting Information Table S2). As shown in the Supporting Information Figure S8, the Fe-N4 sites located on the C57 defect exhibit more negative values of ΔGOH* and the corresponding higher ηORR values, indicating the robust adsorption strength of OH* on these models, which corresponds to a higher overpotential. Figure 2 | The Gibbs free energy diagrams of ORR. The Gibbs free energy diagrams for 4e−ORR on FeN4 sites embedded the (a) high density C585 defect, (b) C484 defect. (c) Relationship curve of ORR (ηORR) as a function of MFe. (d) Activity volcano of ORR (ηORR) as a function of the valence state of Fe site. Download figure Download PowerPoint Recent research indicates that modulation of the spin magnetic moment of the active center can regulate ORR activity on SACS.29–34 Meanwhile, it has been proved that the Perdew–Burke–Ernzerhof (PBE) method can reliably generate the spin-related conclusions.35 Accordingly, we investigated the relationship between the spin magnetic moment and the overpotential, as presented in Figure 2c. The overpotential of Fe-N@DC systems exhibits an increasing trend with the spin magnetic moment until it reaches a maximum value of 1.04 V at 1.89 μB, followed by a decreasing trend to its minimum value of 0.61 V at 2.10 μB. The magnetic moment of iron (Fe) in its most stable state for each structure is listed in Supporting Information Table S3, and C585 (3h) was selected to test the reasonableness of the magnetic moment selection in Supporting Information Table S4. Then, the overpotential as a function of the valence state of the Fe atom is also depicted in Figure 2d, revealing a distinct volcano curve. In Figure 2d and Supporting Information Table S2, at the Fe valence state of 1.2, C585 (3h) occupies the top spot on the volcano curve, boasting the most minimal overpotential of ηORR = 0.61 V among all supports, indicating its superior ORR catalytic activity. Overall, periodic defects can modify the electronic structure of the active central atom of FeN4, which further affects the activity of the Fe-N4 sites. As mentioned in the previous section, periodic defect boundaries in graphene can modulate the activity of ORR catalytic reaction by adjusting the magnetic moment and valence state of the active center Fe atom. To further elucidate its underlying mechanism, we selected four structures—C585(3h), C57(1), C57(2), and C484(1)—for a detailed analysis of their electronic structure, as indicated by the red circles in Figure 2c,d. At first, we chose to examine how the magnetic moment affects the electrocatalytic performance of Fe-N-C catalysts by focusing on the C585(3h), C57(2), and C484(1) structures. Figure 2c illustrates that C57(2) and C484(1) exhibit higher overpotentials, corresponding to lower spin magnetic moments. Our calculated DOS in Figure 3a reveals an increased electron density in the spin-up channel below the Fermi energy (EF) of C585(3h), compared to C57(2) and C484(1), resulting in an increase in magnetic moment. Figure 3 | The origin of activity trends associated with the spin magnetic moment. (a) The reaction-free energy diagrams of ORR for C585(3h), C57(2), and C484(1). (b) Fe 3d DOS of C57(2) and C484(1). (c) The PDOS of Fe 3d orbitals for C585(3h), C57(2), and C484(1). The black dash line represents the Fermi energy level. (d) DOS of the five Fe 3d orbitals in C585(3h). The CDD as well as Bader charge analysis of (e) C57(1), (f) C57(2), and (g) C585(3h) for OH. The yellow (cyan) represents charge accumulation (deletion). The isosurface value is set to be 0.005 e/Å3. Download figure Download PowerPoint As depicted in Figure 3a,b, the analysis of the origin of the high catalytic activity of C585(3h) suggests that Fe sites on C585(3h) exhibits minimal energy difference between the HOMO and LUMO, indicating superior electron transport performance and an accelerated ORR reaction. Moreover, a positive quasi-linear correlation exists between ηORR and the energy gap values of these Fe sites (Figure 3c). The reduced band gaps of the Fe site contributed to enhancing the catalytic activity by lowering the activation energy for catalyzing the ORR.36,37 Additionally, as illustrated in Figure 3a, the DOS reveals that the state labeled with a red circle of C585(3h) is closer to the Fermi energy level than that of C484(1), indicating superior ORR activity in C585(3h) due to more active electrons available for reduction reactions. Figure 3d provides additional evidence supporting the notion that this distinct state correlates with the dz2 electrons of the Fe atom, aligning with findings from prior research.15,32 Therefore, the projected density of states plot provides evidence that periodic defects affect the spin magnetic moment of Fe by altering the distribution of spin electrons within the Fe 3d orbitals and subsequently influences ORR activity. To deeply understand the superior catalytic activity of C585(3h) shown in Figure 2d, we conducted charge density difference (CDD) and Bader charge analyses on OH adsorption for C585(3h), C57(1), and C57(2) carbon structures,38,39 as depicted in Figure 3e–g. The results from the CDD analysis indicate that the adsorbed intermediate (OH*) undergoes charge transfer and redistribution with the FeN4 surface. The charge accumulation of O on C585(3h)-OH* was observed to be the largest (as indicated by the yellow region in Figure 3g). The O atoms attached to C57(1)-OH*, C57(2)-OH*, and C585(3h)-OH* gain electrons from the directly bonded Fe and H atoms, with gains of 0.91 e, 0.88 e, and 1.01 e, respectively, as illustrated in Figure 3e–g. Based on the analysis presented in Figure 2b, we inferred that a larger electron accumulated by O in the adsorption intermediate OH* is advantageous for the subsequent desorption of OH, thereby further enhancing the catalytic activity at the Fe-N4 site. Furthermore, Supporting Information Figure S9 depicts the valence states of Fe atoms, N atoms, and surrounding C atoms in the three distinct structures—C585(3h), C57(1), and C57(2). Notably, the valence states of N atoms adjacent to the Fe site are more negatively charged in the C585(3h) structure compared to those in C57(1) and C57(2), exhibiting a consistent trend with their overpotential change as shown in Figure 2b. Considering the results obtained above from theoretical studies, Fe-N@DC rich in defects was obtained by pyrolyzing Fe-doped zeolitic imidazolate framework-8 with NaCl and LiCl as the template. In contrast, Fe-N@C with few defects was synthesized via the same method, except that there was no NaCl and LiCl in the synthesis process. The scanning electron microscopy (SEM) image shows that Fe-N@DC is in the form of a thin nanosheet due to the addition of NaCl and LiCl template, which is different from polyhedral particles of precursor (Figure 4a). The X-ray diffraction (XRD) pattern of Fe-N@DC has two broad peaks at 24.83° and 43.44°, which correspond to the (101) and (002) crystal planes of graphite, respectively (Figure 4b).40 In addition, there is no Fe-related peak, indicating that there is no crystalline Fe species in the synthesized catalysts. The corresponding Raman spectroscopy reveals the defective carbon (D) at 1337.9 cm−1 and the graphite carbon (G) at 1592 cm−1 (Figure 4c).41 The higher ID/IG value confirms the higher defect density of Fe-N@DC. The electron paramagnetic resonance (EPR) pattern of Fe-N@DC shows an obvious signal at the magnetic field value of 2.005 G, which proves that the defect degree of Fe-N@DC is higher than that of Fe-N@C (Figure 4d).42 Furthermore, the full X-ray photoelectron spectroscopy (XPS) spectrometer reveals the existence of Fe, N, O, and C elements in Fe-N@DC and Fe-N@C ( Supporting Information Figure S8). Moreover, the higher defect degree of Fe-N@DC is also confirmed by the higher sp3/sp2 of C 1s XPS profile (0.93 for Fe-N@DC, 0.76 for Fe-N@C), which corresponds to the Raman and EPR analysis (Figure 4e).43 The peaks of the Fe 2p XPS spectrum at 709.9 eV can be assigned to Fe 2p3/2 ( Supporting Information Figure S9). The N 1s spectrum is deconvoluted into pyridinic N (398.6 eV), Fe-N (399.3 eV), pyrrolic N (400.3 eV), and graphitic N (401.4 eV) (Figure 4f).44 The existence of Fe-N peaks indicates the coordination of Fe and N in the carbon support. Figure 4 | (a) SEM image of Fe-N@DC. (b) XRD pattern, (c) Raman spectrum, (d) EPR spectra, (e) C 1s spectra, and (f) N 1s spectra of Fe-N@DC and Fe-N@C. (g) LSV polarization curves and (h) Tafel slope for Fe-N@DC, Fe-N@C, and 20% Pt/C. (i) LSV polarization curves before and after 10,000 cycles over Fe-N@DC. Download figure Download PowerPoint ORR performance was evaluated by rotating disk electrode (RDE) in the three-electrode system and in the O2-saturated 0.1 m HClO4 solution. The cyclic voltammetry (CV) profile shows the obvious reduction peak of Fe-N@DC at 0.72 V versus reversible hydrogen electrode (RHE), higher than that of Fe-N@C (0.66 V) ( Supporting Information Figure S10). Linear sweep voltammetry (LSV) polarization curves reveal that the half-wave potential (E1/2) of Fe-N@DC is 0.82 V versus RHE, which is higher than that of commercial Pt/C (0.81 V) and Fe-N@C (0.79 V) (Figure 4g and Supporting Information Figure S11), indicating the excellent ORR activity of Fe-N@DC. Compared with Pt/C (107.2 mV s−1) and Fe-N@C (102.4 mV s−1), Fe-N@DC exhibits the smallest Tafel slope (84.26 mV s−1), suggesting the faster ORR kinetics (Figure 4h). The electrochemical stability of the catalyst is further evaluated by an accelerated durability test (ADT). The E1/2 of Fe-N@DC only decreases by 13 mV after 10,000 cycles, much lower than that of Fe-N@C (49 mV) (Figure 4i and Supporting Information Figures S12 and S13). Fe-N@DC with a high defect degree shows better ORR activity and stability than Fe-N@C with pristine graphene structure, corresponding to DFT calculation results. Conclusion In summary, our DFT results demonstrate that periodic C585 and C484 defects can significantly enhance the ORR activity of supported FeN4 sites compared to the pristine graphene structure. C585(3h) exhibits the lowest ORR overpotential (0.61 V), with the magnetic moment of the Fe atom of 2.10 μB. Additionally, we identified a volcano curve that successfully depicts the strong correlation of the spin and valence states of Fe atoms with the corresponding ORR overpotential. These findings suggest that the spin and charge states of the Fe atoms of Fe-N-C catalysts serve as effective catalytic descriptors of ORR. Electronic structure analysis further discloses that the periodic defects contribute to enhanced catalytic activity by reducing the energy gap between HOMO and LUMO, facilitating increased electron accumulation of oxygen atoms in the adsorbed intermediate OH*. Experimentally, the FeN4 sites with a high degree of defects further manifest the better ORR activity and stability compared to FeN4 sites of pristine graphene and commercial Pt/C. Our work provides insights into the ORR activity of Fe-N4 sites with a periodic defect boundary and theoretical guidance for the rational design of efficient graphene-based single-atom electrocatalysts by defect engineering. Methods Computational methods The calculations were carried out using the DFT method with plane waves, as implemented in the Vienna ab initio simulation packages (VASP) software.45,46 The interaction between ions and electrons was modeled using the projector augmented wave approach.47,48 A plane-wave basis set with a cutoff energy of 400 eV was used to expand the Kohn-Sham valence states. The generalized gradient approximation as implemented in the PBE functional is used for the exchange–correlation energy functional.49 In the optimization of the structure, we employed a Monkhorst-Pack 3 × 3 × 1 mesh for the Brillouin zone integration.50 Electronic structural computations utilized a 5 × 5 × 1 k-points, with additional convergence tests employing a larger mesh revealing minimal energy deviations (refer to Supporting Information Table S1). Van der Waals interactions were characterized using the DFT-D3 method with Becke-Johnson damping to correct for dispersion effects.51 A 15 Å vacuum along the z-direction was utilized to prevent periodic interaction. During the optimization of geometry, all atoms were relaxed until both the minimum energy and minimum force reached convergence criteria of 10−4 eV and 0.05 eV/Å, respectively. In an implicit solvent model, VASPsol, was employed because the reaction typically occurs in an aqueous environment.52 The transfer of charge between the substrate and adsorbates was quantitatively evaluated using the Bader charge analysis.38,39 More details concerning the Gibbs free energy calculation are shown in the Supporting Information. The formation energy (Ef) was calculated using a formula of: E f = E sub + 6 μ c − ( E Gra + 4 μ N + E Fe )where Esub and EGra are the total energies of the FeN4 systems and three types of graphene with periodic defect boundaries, respectively. The chemical potential of a single C atom in graphene is represented by μ c , μ N is defined as half the energy of an N2 molecule. The energy of an isolated Fe atom is represented by E Fe . The bind energy (Eb) was calculated using a formula of: E bind = E Fe + sub + E sub − E Fe where EFe+sub and Esub denote the overall energies of catalysts with or without a metal single atom, and EFe corresponds to the isolated Fe atom within the Fe bulk. Experimental methods Electrochemical measurements Details of the synthesis of FeN4 rich in defects (Fe-N@DC) and characterization are shown in the Supporting Information. The electrochemical workstation (Autolab PGSTAT302N, Metrohm, Switzerland) with a RDE was used to evaluate the electrocatalytic activity in a three-electrode electrochemical test system at room temperature. The glassy carbon electrode (0.196 cm2) covered with catalyst film is used as the working electrode, and the Pt plate and RHE are used as the counter electrode and reference electrode respectively. For the preparation of catalyst ink, 5 mg of catalyst was dispersed in ethanol solution containing Nafion (DuPont, United States) followed by ultrasonication for 30 min. Then, 10 μL of catalyst ink was piped out onto the clean glassy carbon surface and dried at room temperature to obtain the working electrode. All electrochemical tests were carried out in O2-saturated 0.1 M HClO4 electrolyte and at room temperature. CV curves were collected from 0.05 to 1.1 V at a scan rate of 50 mV s−1. The LSV curves were collected at 1600 rpm with a sweep rate of 10 mV s−1. The ADT was conducted by performing cyclic potential sweeps between 0.6 and 1.1 V at a sweep rate of 100 mV s−1 for 10,000 cycles. Supporting Information Supporting Information is available free of charge at CCS Chemistry and includes optimized structure, additional convergence details, additional dataset details, DOS, and AIMD result presentation. Conflict of Interests The authors declare no competing interests. Acknowledgments Y. Su acknowledges the "Young Talent Support Plan" of Xi'an Jiaotong University. Supercomputing facilities were provided by Hefei Advanced Computing Center. References 1. Wang H.; Liu R.; Li Y.; Lü X.; Wang Q.; Zhao S.; Yuan K.; Cui Z.; Li X.; Xin S.; Zhang R.; Lei M.; Lin Z.Durable and Efficient Hollow Porous Oxide Spinel Microspheres for Oxygen Reduction.Joule2018, 2, 337–348. Google Scholar 2. Yang G.; Zhu J.; Yuan P.; Hu Y.; Qu G.; Lu B.-A.; Xue X.; Yin H.; Cheng W.; Cheng J.; Xu W.; Li J.; Hu J.; Mu S.; Zhang J.-N.Regulating Fe-Spin State by Atomically Dispersed Mn-N in Fe-N-C Catalysts with High Oxygen Reduction Activity.Nat. Commun.2021, 12, 1734. Google Scholar 3. Li L.; Huang S.; Cao R.; Yuan K.; Lu C.; Huang B.; Tang X.; Hu T.; Zhuang X.; Chen Y.Optimizing Microenvironment of Asymmetric N,S-Coordinated Single-Atom Fe via Axial Fifth Coordination toward Efficient Oxygen Electroreduction.Small2022, 18, 2105387. Google Scholar 4. Jing H.; Zhu P.; Zheng X.; Zhang Z.; Wang D.; Li Y.Theory-Oriented Screening and Discovery of Advanced Energy Transformation Materials in Electrocatalysis.Adv. Powder Mater.2022, 1, 100013. Google Scholar 5. Chen M.; Li X.; Yang F.; Li B.; Stracensky T.; Karakalos S.; Mukerjee S.; Jia Q.; Su D.; Wang G.; Wu G.; Xu H.Atomically Dispersed MnN4 Catalysts via Environmentally Benign Aqueous Synthesis for Oxygen Reduction: Mechanistic Understanding of Activity and Stability Improvements.ACS Catal.2020, 10, 10523–10534. Google Scholar 6. Yang L.; Cheng D.; Xu H.; Zeng X.; Wan X.; Shui J.; Xiang Z.; Cao D.Unveiling the High-Activity Origin of Single-Atom Iron Catalysts for Oxygen Reduction Reaction.Proc. Natl. Acad. Sci. U.S,A,2018, 115, 6626–6631. Google Scholar 7. Peng H.; Liu F.; Liu X.; Liao S.; You C.; Tian X.; Nan H.; Luo F.; Song H.; Fu Z.; Huang P.Effect of Transition Metals on the Structure and Performance of the Doped Carbon Catalysts Derived From Polyaniline and Melamine for ORR Application.ACS Catal.2014, 4, 3797–3805. Google Scholar 8. Han J.; Meng X.; Lu L.; Bian J.; Li Z.; Sun C.Single-Atom Fe-Nx-C as an Efficient Electrocatalyst for Zinc–Air Batteries.Adv. Funct. Mater.2019, 29, 1808872. Google Scholar 9. Jiang R.; Li L.; Sheng T.; Hu G.; Chen Y.; Wang L.Edge-Site Engineering of Atomically Dispersed Fe–N4 by Selective C–N Bond Cleavage for Enhanced Oxygen Reduction Reaction Activities.J. Am. Chem. Soc.2018, 140, 11594–11598. Google Scholar 10. Li J.; Sougrati M. T.; Zitolo A.; Ablett J. M.; Oğuz I. C.; Mineva T.; Matanovic I.; Atanassov P.; Huang Y.; Zenyuk I.; Di Cicco A.; Kumar K.; Dubau L.; Maillard F.; Dražić G.; Jaouen F.Identification of Durable and Non-Durable FeNx Sites in Fe–N–C Materials for Proton Exchange Membrane Fuel Cells.Nat. Catal.2021, 4, 10–19. Google Scholar 11. Lin Y.; Liu