Open Access Article
This Open Access Article is licensed under a Creative Commons Attribution-Non Commercial 3.0 Unported Licence

Manipulating the spin state to activate the atomically dispersed Fe–N–C catalyst for oxygen reduction

Fan Liu abc, Chengxiang Shi abc, Lun Pan abc, Zhen-Feng Huang *abc, Xiangwen Zhang abc and Ji-Jun Zou *abcd
aKey Laboratory for Green Chemical Technology of the Ministry of Education, School of Chemical Engineering and Technology, Tianjin University, Tianjin 300072, China. E-mail: jj_zou@tju.edu.cn; zfhuang@tju.edu.cn
bCollaborative Innovative Center of Chemical Science and Engineering (Tianjin), Tianjin 300072, China
cZhejiang Institute of Tianjin University, Ningbo, Zhejiang 315201, China
dQinghai Quality Certification Consulting and Inspection Center Co., Ltd, Qinghai, China

Received 27th March 2023 , Accepted 21st April 2023

First published on 25th April 2023


Abstract

Atomically dispersed metal–nitrogen carbon (M–N–C) catalysts with e1g are believed to be very active for the oxygen reduction reaction; however, the activity origin of spin state manipulation is unclear. Here, using Fe–N–C as a model catalyst, we have developed a facile approach to manipulate its spin state from the low-spin to medium-spin state (FeSAC–NC, e1g) by utilizing the secondary coordination sphere effect of adjacent N[double bond, length as m-dash]C–N moieties. The presence of N[double bond, length as m-dash]C–N configurations with high electronegativity modifies the degree of hybridization between Fe 3dz2, 3dxz (3dyz), and oxygen-containing intermediate π* orbitals, leading to an ideal balance between O2 activation and *OH desorption. More impressively, we found that the catalysts with eg = 1 in their metal centres may show very different activity, and the magnetic moment of 3dxz + 3dyz as a descriptor can accurately predict the ORR activity. Benefiting from the regulated coordination and electronic structures, the designed FeSAC–NC catalyst exhibited 3-fold mass activity and 7-fold specific activity than Pt/C.



Broader context

Fuel cells have been identified as prospective clean power contenders due to their high power density, environmental friendliness, and efficiency. Metal–nitrogen carbon catalysts with atomically dispersed Fe–Nx active sites have garnered a lot of attention as a possible non-noble metal-based catalyst class for the oxygen reduction process. A highly active Fe–N–C catalyst may be obtained via using an atomically dispersed metal framework that accurately controls the spin state of the active metal. However, given the different catalytic activities of catalysts with the identical spin state, it is highly difficult to adequately establish the structure–activity relationship between the catalytic activity and active site structure. In this study, we unveil the activity origin of atomically dispersed Fe–N–C catalyst by spin state manipulation and construct the targeted Fe–N–C catalysts with medium-spin state through utilizing the secondary coordination sphere effect of adjacent N[double bond, length as m-dash]C–N moieties. Actually, the catalyst with eg = 1 may show very different activity while magnetic moment of 3dxz + 3dyz is a more accurate descriptor to predict activity. This work fills in the gaps of eg theory and unveils the ORR activity origin of M–N–C catalysts, providing new insights into precise customization of electronic structures of catalysts.

1. Introduction

The increasing energy crisis and environmental pollution caused by burning fossil fuels have motivated scientists to explore renewable energy conversions technologies such as fuel cells and metal–air batteries.1–5 Oxygen reduction reactions (ORR), as the key reaction, involves multi-step proton and electron transfer and necessitate significant overpotential due to the sluggish reaction kinetics, which has become a major bottleneck.6,7 Generally, noble metals such as Pt are the most active to accelerate ORR kinetics. Considering the scarcity, price, and low methanol crossover tolerance of noble metals, the development of stable and efficient non-precious metal catalysts is particularly important.8,9

M–N–C catalysts derived from low-cost and earth-abundant elements are very promising for ORR,10,11 where the filling degree of d orbital in the metal center dominated by the oxidation state and spin state plays a central role in ORR. However, the active sites usually exhibit unsuitable spin states; for example, FeN4 typically exhibits a low spin state, which has a relatively strong adsorption strength of oxygen-containing intermediates, limiting the O2 activation and *OH desorption.12 Actually, the reduction of O2 preferentially occurs on the Fe(III) site with a medium spin state, where one eg electron (t42ge1g) can readily penetrate the antibonding π-orbital of oxygen, resulting in relatively weak adsorption of oxygen-containing intermediates. Many efforts have been devoted to optimizing the spin-state of the Fe center, such as introducing a second metal,13,14 doping electronegative non-metal elements (e.g. N, P, S, O),15,16 and applying an external magnetic field.17 These methods bring many complicated factors to the activity, however the investigation is solely focused on promoting O2 activation or *OH desorption; thus, many important factors may be ignored. More significantly, the catalytic activity of medium spin (e1g) catalysts is not identical,18,19 suggesting that the eg theory does not apply in cases where there are several putative active sites with e1g and thus cannot adequately describe the ORR activity origin. Thus far, it is challenging to thoroughly establish the structure–activity relationship between the catalytic activity and active site structure, particularly its spin state, let alone accurately fabricating specific spin state catalysts. Therefore, a better knowledge of the active origin of Fe–N–C catalysts is required to design better ORR catalysts.

Theoretically, the local charge distribution and spin state of the Fe active center can be modulated and its ORR intrinsic activity can be thus improved by manipulating the near-range coordinated atom species (primary coordination sphere) or long-range heteroatom interactions (secondary coordination sphere). Herein, we report an approach to manipulate the spin state of the Fe center by tuning the secondary coordination spheres in atomically dispersed Fe–N–C catalysts. 57Mössbauer spectroscopy, X-ray absorption spectroscopy and density functional theory (DFT) calculations revealed that the insertion of neighbouring N[double bond, length as m-dash]C–N configurations with high electronegativity can drive the spin state of Fe(III) transition from a thermodynamically stable low spin to a medium spin and modify the degree of hybridization between Fe 3dxz (3dyz) and oxygen-containing intermediates π* orbitals, which can optimize the adsorption energy for oxygen-containing intermediates, thereby leading to an ideal balance between O2 activation and *OH desorption. More importantly, the catalysts with e1g in their metal centers are found to exhibit significantly diverse activity, while the magnetic moment of 3dxz + 3dyz as a descriptor can precisely predict the ORR activity. Benefiting from the regulated coordination and electronic structures, FeSAC–NC exhibits excellent ORR catalytic performance and stability in alkaline electrolytes, surpassing conventional FeN4 and commercial Pt/C catalysts. In a homemade zinc–air battery, the FeSAC–NC electrode delivered a higher maximum power density and discharge specific capacity than commercial Pt/C. This study fills in the gaps of the eg theory and unveils the ORR activity origin of M–N–C catalysts, providing new insight into the precise customization of the electronic structures of catalysts.

2. Experimental

Synthesis of electro-catalyst

For the preparation of nitrogen-doped carbon, HBCT (20 mmol) and melamine (20 mmol) were dispersed separately in 50 mL of distilled water at 90 °C. Then, the two solutions were mixed and stirred for 20 min, followed by natural cooling to room temperature under stirring. The polymer was collected by filtration and dried, subsequently annealed at 400 °C for 2 h with a heating rate of 5 °C min−1 under a flowing N2 atmosphere. The obtained nitrogen-doped carbon was denoted as NC. The as-synthesized NC material (500 mg) was dispersed in 200 mL of FeCl2·4H2O (0.5 mM) aqueous solution and stirred for 1 h. Then, the impregnated NC was separated by centrifugation and freeze-dried. After drying overnight, the dried material was annealed at 700 °C for 2 h with a heating rate of 5 °C min−1 under a flowing N2 gas. Afterward, the sample was naturally cooled to room temperature and labelled as FeSAC–NC. At the same time, a typical FeN4 catalyst was synthesized according to the literature for comparison.20

Characterizations

Scanning electron microscopy (SEM) observations were carried out on a Hitachi S-4800 microscope with an acceleration voltage of 530 KV. Transmission electron microscopy (TEM), high-resolution transmission electron microscopy (HR-TEM), and electron energy loss spectrum (EELS) observations were carried out on the Tecnai G2 F-20 microscope with an acceleration voltage of 200 KV. The high-angle annular dark field scanning transmission electron microscopy (HAADF-STEM) was carried out on a JEOL JEM-2100F with an acceleration voltage of 200 kV. X-ray photoelectron spectroscopy (XPS) analysis was conducted on a PHI-1600 X-ray photoelectron spectrometer using Al Kα radiation, and C 1s peak (284.8 eV) of contamination carbon was adopted to calibrate binding energy. The X-ray absorption fine structure spectroscopy (XAFS) measurement was carried out at Beijing Synchrotron Radiation Facility with a 1W2B beamline, using an incident beam monochromatized by Si (111) double crystal monochromators. X-ray diffraction (XRD) patterns were collected on a D8-Focus X-ray diffractometer system with Cu Kα radiation (λ = 1.5419 Å) with a scanning rate of 5° min−1. Fourier transform infrared (FT-IR) spectrum was recorded on a Bruker FTIR spectrophotometer. Nitrogen physisorption experiments were conducted on a Micomeritics ASAP 2420 volumetric absorption analyzer at −196 °C.

Electrochemical measurements

All electrochemical measurements were carried out on the IVIUMSTAT workstation (Ivium Technologies BV, Netherlands) under the three-electrode system at room temperature. The catalyst was supported on a glassy carbon electrode (GCE) with a diameter of 5 mm (Pine Research Instrumentation Inc.) was used as a working electrode, whereas graphite rod and Ag/AgCl electrode were used as counter and reference electrodes, respectively. The cyclic voltammetry (CV) was carried out in O2 or Ar-saturated 0.1 M KOH with a sweep rate of 50 mV s−1. Linear sweep voltammetry (LSV) polarization curves were measured at 1600 rpm in O2 saturated 0.1 M KOH electrolyte with a sweep rate of 10 mV s−1. The rotating ring disk electrode (RRDE) analysis was also run in Ar or O2 saturated 0.1 M KOH electrolyte with a rotation speed of 1600 rpm and the applied ring voltage was 1.125 V vs. RHE. LSVs rotated between 400 and 2025[thin space (1/6-em)]rpm were obtained under an uninterrupted O2 flow, in which the sweep window and rate were 0.2 to 1.2[thin space (1/6-em)]V vs. RHE and 10[thin space (1/6-em)]mV s−1, respectively. Catalyst inks were prepared by dispersing 2.5 mg of the catalyst and 0.5 mg carbon in 980 μL of a mixed solvent (water and isopropanol in a 3[thin space (1/6-em)]:[thin space (1/6-em)]1 volume ratio) and 20 μL of a 5 wt% Nafion solution, with 30 min of sonication to obtain a uniform suspension. Then, 10 μL of homogeneous ink was dropped on a glassy carbon electrode (area, 0.196 cm2) and fully dried in air at room temperature. The mass loading of the electrocatalyst on glassy carbon was 0.128 mg cm−2.

DFT calculations

All DFT calculations were performed using the Vienna ab initio simulation package (VASP).21 The generalized gradient approximation (GGA) and Perdew–Burker–Ernzerhof (PBE) functional were used to address the nonlocal exchange-correlation energy. The DFT+U method with U = 4.0 eV was applied to eliminate Fe 3d orbital self-interaction error. The cutoff energy and Brillouin zone integration were set to 520 eV and 2 × 2 × 1, respectively. The force and energy convergence criteria of structure optimization were 0.03 eV Å−1 and 1 × 10−5 eV. All electronic structure calculations were performed using a 6 × 6 × 1 k-point grid. Spin polarization was considered in our calculations. To eliminate spurious periodic interactions, a 15 Å vacuum slab was added in the z-direction.

3. Results and discussion

3.1 Synthesis and characterization of FeSAC–NC

FeSAC–NC was synthesized via a combined hydrothermal and pyrolysis strategy, as illustrated in Fig. 1(a). First, 1,3,5-benzenetricarboxylic acid (HBCT) and melamine were mixed at 90 °C to enable self-assembly. As shown in FT-IR spectroscopy (Fig. 1(b)), the decrease in the intensity of the peaks associated with amine and carboxy groups of the obtained polymer revealed the amidation reaction between melamine and HBCT. Subsequently, porous nitrogen-doped carbon (NC) was obtained by pyrolyzing the polymer at 400 °C. Next, NC and FeCl2·4H2O were dispersed in a solvent and freeze-dried, ensuring the impregnation of ferric salt on NC. Finally, FeSAC–NC was obtained by pyrolyzing impregnated NC at 700 °C under an N2 atmosphere. The nitrile group (–CN) at 2235 cm−1 of NC was not detected in FeSAC–NC, suggesting that the –CN group underwent a transformation and became involved in the creation of FeNx active sites. The FeN4 catalyst was synthesized as a reference using the same synthetic strategy.20
image file: d3ey00066d-f1.tif
Fig. 1 (a) Schematic representation of preparation for FeSAC–NC. (b) FT-IR spectra of HTBC, melamine, prepolymer, NC, FeN4, and FeSAC–NC. (c) and (d) TEM and HRTEM images of FeSAC–NC. (e) Aberration-corrected HAADF-STEM image of FeSAC–NC (left). Middle: The enlarged image of FeSAC–NC as marked by the red box. Right: Intensity distributions between two Fe atoms. (f) EDS element mapping of FeSAC–NC. (g) XRD pattern of FeSAC–NC and FeN4. (h) Raman spectra of FeSAC–NC and FeN4.

From scanning electron microscopy (SEM) images (Fig. S1, ESI), the NC presents a bark-like structure on the surface. After adsorbing Fe and pyrolysis in the N2 atmosphere, FeSAC–NC displays a flower-like nanostructure with a large number of mesopores. Bright-field transmission electron microscopy (TEM) images in Fig. 1(c) and (d) display amorphous carbon substrate without visible Fe nanoparticles in FeSAC–NC. Electron energy loss spectrometer (EELS) elemental analyses (Fig. S2, ESI) confirmed the loading of Fe element in the porous carbon skeleton. The aberration-corrected high-angle annular dark field-scanning transmission electron microscopy (HAADF-STEM) was performed to observe their structures in detail, and the results obtained are shown in Fig. 1(e). Due to the contrast difference of Fe with N and C, single Fe atoms are directly observed in the form of isolated brighter spots and homogeneously dispersed on the entire carbonaceous support. Furthermore, the space between two Fe atoms was 0.8 nm (Fig. 1(e), right) far exceeding the Fe–Fe bond length (0.248 nm), which consolidated the isolated atom feature of Fe in FeSAC–NC. EDS elemental mapping images (Fig. 1(f)), further manifest the uniform dispersion of C, N, and Fe elements in FeSAC–NC. X-ray diffraction (XRD) was employed to investigate the crystalline structure. As illustrated in Fig. 1(g), two broad peaks at approximately 25.0° and 44.0° correspond to (002) and (101) phases of carbon, respectively (JCPDS no. 41-1487).22 Both FeSAC–NC and FeN4 showed negligible Fe or FeOx signals, indicating that the Fe atoms are in atomic separation. As viewed in Fig. 1(h), the ID/IG value of FeSAC–NC is 1.02, with a graphitization degree similar to FeN4 (ID/IG = 1.05). The BET surface area of FeSAC–NC is 725 m2 g−1, much higher than that of NC (255 m2 g−1), suggesting increased surface area by creating mesopores. However, compared with FeN4, which features a high surface area (752 m2 g−1), FeSAC–NC possessed a slightly lower surface area (Table S1, ESI). The isotherm curves of FeSAC–NC present a typical IV profile with a classical H4-type hysteresis loop,23,24 revealing the well-developed micropores and the mesoporous structure (Fig. S3 and S4, ESI). The formation of a mesopore-dominated porous structure in FeSAC–NC is attributed to the thermal decomposition of NC in the second heat treatment. The mesopore-dominated structure can boost the ORR performance by providing more nanochannels for fast reactant transportation and more active sites, which can readily adsorb and react with the reactants. The metal content in FeSAC–NC was analyzed using inductively coupled plasma atomic emission spectrometry (ICP), as shown in Table S2 (ESI). The ICP combined with the elemental analyzer indicates the coexistence of C (81.38 wt%), N (8.97 wt%), O (3.21 wt%), Fe (5.79 wt%), and H (0.96 wt%) in FeSAC–NC. The Fe content measured from ICP is similar to that obtained from TGA analysis (Fig. S5, ESI). Notably, the N/Fe atomic ratio of FeSAC–NC is approximately 6[thin space (1/6-em)]:[thin space (1/6-em)]1, significantly higher than that of FeN4 (4[thin space (1/6-em)]:[thin space (1/6-em)]1), which means abundant N species are embedded in the structure.

X-ray absorption spectroscopy (XAS) measurements including X-ray absorption near-edge structure (XANES) and extended X-ray absorption fine structure (EXAFS) spectra were conducted to detect the electronic and coordination structure, with Fe foil, Fe2O3, iron(II) phthalocyanine (FePc), and FeN4 as references. The formation of single atom Fe species was confirmed from Fe K-edge EXAFS spectra, where both FeSAC–NC and FeN4 present a primary peak corresponding to the Fe–N scattering path (Fig. 2(a)). However, the major Fe–N peak in FeSAC–NC situated at 1.47 Å is larger than that in FeN4 (1.44 Å). The quantitative EXAFS fitting result shows that the Fe atom in FeSAC–NC is coordinated with four-N atoms, forming a Fe–N4 moiety in the first shell (Fig. S6 and Table S3, ESI). In the pre-edge region, a peak specified as a 1s–3d transition at 7113 eV is detected in FePc and FeN4, which is regarded as the fingerprint of the square-planner Fe–N4 configuration. A decreased pre-edge peak intensity of FeSAC–NC implies the emergence of a broken square-planar structure with D4h symmetry owing to the Fe–C scattering path.25 The normalized XANES spectra show that the near-edge absorption of the Fe K-edge of FeSAC–NC is between FePc and Fe2O3, indicating that the average oxidation state of Fe is between +2 and +3 (Fig. 2(b) and Table S4, ESI). Meanwhile, the absorption edge position for FeN4 is similar to that of FePc, elucidating that the oxidation state of Fe is approximately +2. This suggests that the oxidation state of Fe in FeSAC–NC is higher than that in FeN4.


image file: d3ey00066d-f2.tif
Fig. 2 (a) The k2-weighted Fourier transform EXAFS spectra of FeSAC–NC, Fe2O3, Fe foil, FePc, and FeN4. (b) The normalized XANES spectra at the Fe K-edge. (c) The normalized XANES spectra at the N K-edge. (d) The model of FeSAC–NC (gold: iron, silvery: nitrogen, brown: carbon). (e) XPS spectra of Fe 2p. (f) XPS spectra of C 1s. (g) Wavelet transform plots of Fe foil, FeN4, and FeSAC–NC. (h) XPS spectra of N 1s. (i) Room-temperature 57Fe Mössbauer spectrum of FeSAC–NC. (j) Schematic representation of the spin transition of Fe(III).

The N K-edge XANES spectra were collected to probe the N coordination structure. As shown in Fig. 2(c), three peaks were observed for FeN4, denoted as N1 to N3. The peak N1 at 399.8 eV is attributed to the N 1s → π* transition in aromatic N atoms of heterocyclic rings, π*(C[double bond, length as m-dash]N–C). The peak N2 at 401.7 eV is assigned to the graphitic threefold N atoms π*(N–3C). The peak N3 could not be assigned to any N species, derived from the general transition from the N 1s core level to C–N σ* states.26 It is worth noting that the peak N2 splits into double peaks for FeSAC–NC, confirming that a portion of graphitic N is bonded to extra N atoms and forms N[double bond, length as m-dash]C–N coordination-like DNA at 401.4 eV,27,28 rather than direct coordination with metal Fe in the first shell, which is substantially different from FeN4. Fig. 2(d) shows a detailed schematic of the N[double bond, length as m-dash]C–N configuration in FeSAC–NC. The bonding between C[double bond, length as m-dash]N and N atoms, thereafter initiates a slight asymmetrical expansion of the original Fe–N bonds (FeN4: 1.89 Å and FeSAC–NC: 1.90 Å), which is consistent with the EXAFS spectra (Fig. 2(a) and (d)). The normalized C K-edge XANES spectra of FeN4 and FeSAC–NC are shown in Fig. S7 (ESI). The slight positive shift of the binding energy and the decrease in the peak intensity on π* (C–N) in FeSAC–NC confirm the existence of the N[double bond, length as m-dash]C–N configuration because C exhibits lower electronegativity compared to N. As shown in Fig. 2(g), FeSAC–NC, and FeN4 show WT-maximum with the k values of 4.1 Å−1 and 4.0 Å−1 for the Fe–N path, respectively, which is clearly distinguished from the Fe–Fe path (k = 7.5 Å−1). Note that the deviation of the Fe–N path results presumably from the N[double bond, length as m-dash]C–N configuration for FeSAC–NC.

X-ray photoelectron spectroscopy (XPS) characterization was carried out to probe the composition and chemical structure. As shown in Fig. 2(e), FeN4 shows only two Fe 2p XPS peaks at 711.0 and 724.2 eV derived from Fe(II)–N4.29 In contrast, FeSAC–NC possesses a larger proportion of Fe(III)–N4 (714.5 and 732.2 eV), consistent with the higher oxidation state.30 The Fe(III)–N4/Fe(II)–N4 area ratio is 2[thin space (1/6-em)]:[thin space (1/6-em)]1, implying that Fe is mainly in the form of Fe3+ for FeSAC–NC. Fig. 2(f) displays the C 1s peaks at 284.8 eV, 286.1, and 288.0 eV for FeSAC–NC, which can be assigned to C–C[double bond, length as m-dash]C, C–N, and N–C[double bond, length as m-dash]N species,31,32 respectively, consistent with the C and N K-edge XANES spectra. Fig. 2(h) displays the N 1s peaks at 398.1, 399.2, and 400.5 eV for FeSAC–NC, which can be assigned to pyridine-N, Fe–N4, N–C[double bond, length as m-dash]N (graphite-N), respectively.20,33 Compared with FeN4, the position of pyridine-N and Fe–N4 of FeSAC–NC are shifted towards lower energy, implying that the electrons transfer from Fe to N (Fig. 2(h)). The considerable increase in the graphite-N content in FeSAC–NC over FeN4 further demonstrates that the presence of N–C[double bond, length as m-dash]N configuration.

57Mössbauer spectroscopy was further used to quantify the spin state of atomic Fe–Nx sites. The first two doublets were the characteristics of Fe–N–C catalysts (FeN4) ascribed to a ferrous low-spin Fe(II)–N4 (D1) and a ferrous mediate-spin Fe(II)–N4 (D2).34,35 Previous research identified the D1-related component as the ORR active site in Fe–N–C catalysts.36 Notably, FeSAC–NC contains three components: D1, D3, and sextet1, which are assigned to Fe(II)–N4, Fe(III)–N4, and α-Fe (Fig. 2(i) and Table S5, ESI).37,38 According to the values of quadrupole splitting (QS), the doublets D1 and D3 have been assigned to, mainly, low-spin Fe(II)–N4 species (S = 0) and medium-spin Fe(III)–N4 species (S = 3/2), respectively. Compared to the case of FeN4, it can be seen that the incorporation of the N[double bond, length as m-dash]C–N configuration induces a new medium spin Fe(III). The medium spin Fe3+ has an electron configuration of image file: d3ey00066d-t1.tif and shows lower eg filling with favorable adsorption strength of oxygen-containing intermediates (Fig. 2(j)), and is thus more active for ORR.37 Furthermore, the area of D1 and D3 is 30.2% and 61.7%, respectively, indicating that FeSAC–NC is mainly composed of medium spin Fe(III)–N4 species.

On the basis of the above analysis, we can conclude that the coordination structure of FeSAC–NC involves a Fe central atom surrounded by four nitrogen atoms in the first shell, with an N[double bond, length as m-dash]C–N configuration in the second shell. Such a coordination structure can readily trigger the electron transfer from Fe atoms to N and then to N[double bond, length as m-dash]C due to the relatively high electron affinity of N.

3.2 ORR electrocatalytic performance of FeSAC–NC

Fig. S8 (ESI) shows the cyclic voltammetry (CV) curves of FeSAC–NC in O2 and Ar-saturated solution. In contrast to the virtually featureless CV curves measured in the Ar-saturated solution, a notable cathodic peak appears in the O2-saturated solution, which can be attributed to the reduction of O2. The ORR peak of Pt/C appears at 0.84 V, while the ORR peak of FeSAC–NC reaches 0.89 V, approximately 50 mV higher than that of Pt/C, suggesting that FeSAC–NC shows higher ORR activity than commercial Pt/C.

The linear scanning voltammetry (LSV) was adopted to further investigate the ORR activity, as shown in Fig. 3(a). NC delivers an extremely low activity (Fig. S9, ESI). In Fig. 3(a) and (b), FeSAC–NC displays excellent activity with an onset potential of 1.06 V and a half-wave potential (E1/2) of 0.90 V and even exceeds commercial Pt/C. Remarkably, compared with FeN4, FeSAC–NC exhibits enhanced electrocatalytic performance toward ORR, indicating that the formation of medium spin Fe3+ further improves the ORR activity. As for the reaction kinetics, FeSAC–NC possesses an ultrahigh kinetic current density (Jk) of 44.4 mA cm−2 at 0.80 V, nearly an order of magnitude higher than that of commercial Pt/C (Jk, 7.89 mA cm−2). The superior ORR kinetics of FeSAC–NC was further confirmed by the Tafel slope. The Tafel slopes of FeSAC–NC, 20% Pt/C, and FeN4 were 66.0, 63.4, and 89.5 mV dec−1, respectively, implying that FeSAC–NC possessed more favorable ORR kinetics compared with FeN4 (Fig. 3(c)). The corresponding turnover frequency (TOF) of FeSAC–NC was 44.06 s−1, outperforming that of Pt/C (4.81 s−1),39 suggesting that the FeSAC–NC possessed outstanding intrinsic ORR activity. The mass activity and specific activity of FeSAC–NC are 0.64 mA μgFe−1 and 2.10 mA cmFe−1, respectively, three times and seven times higher than that of Pt/C (0.2 mA μgPt−1 and 0.31 mA cmPt−1). Additionally, the ORR activity of FeSAC–NC in this work surpasses that reported for most near-term reported single atoms-based electrocatalysts (Table S6, ESI).


image file: d3ey00066d-f3.tif
Fig. 3 (a) The LSV curves of FeSAC–NC, FeN4, and 20% Pt/C in 0.1 M KOH. (b) The onset potential and half-wave potential of FeSAC–NC, FeN4, and 20% Pt/C. (c) Tafel plots of FeSAC–NC, FeN4, and 20% Pt/C. (d) Electron transfer number (n) and H2O2 yield calculated from RRDE measurement of FeSAC–NC and 20% Pt/C. (e) Stability evaluation of FeSAC–NC and 20% Pt/C. (f) LSV curves of FeSAC–NC before and after 100k cyclic voltammetry cycles.

Fig. S10a (ESI) shows the LSV curves of FeSAC–NC measured at different rotating speeds. The current density goes up rapidly with the increase in the rotating speed. Fig. S10b (ESI) exhibits Koutecky–Levich (K–L) plots of FeSAC–NC at different polarizing potentials, which show excellent linearity under different potentials, suggesting that the transferred electron numbers per O2 molecule in the ORR process are almost the same. The average electron transfer number (n) of ORR for FeSAC–NC is 3.97 approaching that of Pt/C (n = 3.98), which manifests an efficient four-electron transfer process on FeSAC–NC. As shown in Fig. S11 and S12 (ESI), the calculated average electron transfer number for FeN4 is 3.89 because the Fe centers catalyze the formation of undesirable H2O2 products.30,40 Therefore, medium spin Fe3+ in FeSAC–NC can avoid the possible Fenton reaction. To quantitatively evaluate the intermediate peroxide product, the RRDE measurement was performed. As shown in Fig. 3(d), the measurement conducted on the RRDE also revealed that n for FeSAC–NC was 3.97 and the H2O2 yield of FeSAC–NC was measured to be below 4%, confirming high selectivity for the 4-electron ORR. However, the yield of H2O2 increases significantly in the high overpotential region of FeN4 (Fig. S12, ESI). The electrochemical active surface area of FeSAC–NC measured by the double-layer capacitance (Cdl) was 30.5 mF cm−2, lower than that of FeN4 (34.4 mF cm−2) (Fig. S13 and S14, ESI). All of these verify the key role of medium spin Fe3+ in improving the ORR intrinsic activity.

Chronoamperometric response curves measured for FeSAC–NC and commercial Pt/C are shown in Fig. 3(e). After 80[thin space (1/6-em)]000 seconds at a potential 0.6 V, the commercial Pt/C electrode displayed a dramatical reduction in the voltammetric current by approximately 31.8% with respect to the original value. In sharp contrast, FeSAC–NC retains up to 92.9% of its initial current under identical conditions, which manifests much-improved stability for ORR. The TEM, XRD, and Raman for 80[thin space (1/6-em)]000 seconds of electrocatalysis are provided, as shown in Fig. S15 (ESI). No obvious changes were observed in the above characterizations for post-catalysis FeSAC–NC, which suggested that both structures and properties of the catalyst were well retained. Fig. 3(f) shows the continuous cycling for 100k cycles between 0.6 and 1.1 V. After 100k cycles, the half-wave potential of FeSAC–NC has a slight shift of 10 mV, indicating a robust ORR performance.

The tolerance to methanol crossover is important for practical applications such as direct methanol fuel cells. The methanol tolerance was evaluated by adding a 3 M methanol solution during the stability test (Fig. S16, ESI). FeSAC–NC still maintained a high catalytic activity in the presence of methanol, demonstrating good tolerance to methanol. On the contrary, the commercial Pt/C showed a sharp jump in the current density due to the oxidation reaction of methanol.

We further assembled a zinc–air battery to demonstrate the potential of FeSAC–NC in practical applications. FeSAC–NC-based battery delivers a higher maximum power density and discharge voltage at the same current density than the Pt/C-based battery (Fig. 4(a)), confirming the high ORR activity of FeSAC–NC in the device. As shown in Fig. 4(b), the zinc–air battery with FeSAC–NC cathode shows a discharge-specific capacity of 843 mA h g−1 at a current density of 10 mA cm−2, higher than that of Pt/C (718 mA h g−1).


image file: d3ey00066d-f4.tif
Fig. 4 (a) Discharging polarization curves and the corresponding power density plots of zinc–air battery with FeSAC–NC and Pt/C. (b) Specific capacity of zinc–air battery with FeSAC–NC and Pt/C.

3.3 Activity origin for Fe–N–C with e1g

To understand the origin of the medium spin Fe(III) (e1g) activity in Fe–N–C catalysts, DFT calculations were carried out. Based on ICP, XPS, XAFS, and 57Mössbauer spectroscopy analyses, FeSAC–NC involves a Fe central atom surrounded by four N atoms in the first shell, with N[double bond, length as m-dash]C–N configuration in the second shell. Seventeen possible structures were constructed to assess their ORR activity, as shown in Fig. S17 (ESI), and marked as FeSAC–NXC (X = 1–17), where x stands for different configurations of the N[double bond, length as m-dash]C–N. Both the formation energy of FeSAC–NXC and the binding energy of the single-atom Fe sites indicate that FeSAC–NXC is thermodynamically stable (Fig. S18, ESI). As shown in Fig. S19 (ESI), the electrons are transferred from Fe to N and further migrated to the N[double bond, length as m-dash]C group, resulting in a distinct depletion of charge density on the Fe site, which is consistent with the results of XAS and XPS.

To probe the effect of the N[double bond, length as m-dash]C–N configuration on the electronic structure of metal Fe atoms, the density of states (DOS) of FeSAC–NC and FeN4 near the Fermi level, mainly originating from the 3d state, were investigated. Remarkably, the 3dx2y2 and 3dxy orbitals of all Fe–N–C catalysts are either fully unoccupied or occupied, indicating that these orbitals are not engaged in bonding with oxygen-containing intermediates (Fig. S20, ESI). The magnetic moment, however, comes from the spin splitting of partially occupied 3dz2, 3dxz, and 3dyz orbitals, implying the magnetic moments of these orbitals govern the ORR activity. The magnetic moments of 3dx2y2 + 3dxy, 3dz2, 3dxz + 3dyz orbitals of FeSAC–NXC were calculated, and it is concluded that the presence of the N[double bond, length as m-dash]C–N configurations obviously changes the magnetic moments of 3dxz + 3dyz orbitals and the magnetic moments of the 3dxz + 3dyz orbitals is nearly linear in relationship with that of the magnetic moments of metal Fe (Fig. 5(a) and (b)). Simultaneously, we also found a nearly linear relationship between the magnetic moment of 3dxz + 3dyz orbitals, ΔG*OH (an important ORR descriptor),41 and overpotential, revealing that the magnetic moment of 3dxz + 3dyz orbitals can reflect the origin of ORR activity (Fig. 5(c)). The change in the free energy of *OH@FeSAC–NXC ranged from 0.75 to 0.87 eV, which is much larger than the *OH@Fe–N4 site (0.58 eV). The weak adsorption of *OH on FeSAC–NXC is the key to improving the ORR activity. Compared with FeN4, the magnetic moment of 3dxz + 3dyz of FeSAC–NXC is significantly reduced, which can weaken the interaction between *OH and Fe sites and promote the desorption of *OH. However, the theoretical overpotential of the FeSAC–NXC ranges from 0.37 eV to 0.48 eV, revealing that the activities of the active sites with the e1g electron configuration are not identical. Meanwhile, the Fe–O bond of *OH@FeSAC–NXC decreased with the increase of 3dxz + 3dyz magnetic moment (Fig. 5(d)). This is due to the different hybridization modes brought about by the orbital heights of 3dz2 and 3dxz (3dyz) orbitals.42 Specifically, when the Fe–O bond length is short, 3dxz + 3dyz and 3dz2 orbitals participate in the hybridization. On the contrary, the contribution of 3dxz + 3dyz orbital is greatly reduced when the Fe–O bond is longer (Fig. S21, ESI).


image file: d3ey00066d-f5.tif
Fig. 5 (a) The magnetic moment of 3dz2, 3dxz + 3dyz, and 3dx2y2 + 3dxy orbital at the Fe site on FeSAC–NxC (X = 1–17). (b) Correlation between the magnetic moment of 3dxz + 3dyz and magnetic moment of Fe site. (c) Correlation between the magnetic moment of 3dxz + 3dyz, ΔG*OH (up) and overpotential (down) on FeSAC–NXC (X = 1–17). (d) Correlation between the magnetic moment of 3dxz + 3dyz and the Fe–O bond length on *OH@FeSAC–NXC (X = 1–17). (e) The free energy diagrams of ORR pathways on FeSAC–N1C, FeSAC–N3C, FeSAC–N5C, and FeSAC–N16C at U = 1.23 V vs. RHE. Note: 1–17 here represent FeSAC–N1C to FeSAC–N17C, respectively.

All of these verify that the insertion of neighbouring N[double bond, length as m-dash]C–N configurations with high electronegativity modifies the degree of hybridization between Fe 3dz2, 3dxz, 3dyz, and oxygen-containing intermediates π* orbitals, and weakens the interaction between the oxygen-containing intermediates and Fe centers, thus improving the ORR activity. However, the seventeen possible structures show significantly different activity although they all have e1g, indicating that a more accurate and intrinsic indicator is needed to predict the activity. Actually, our results show a linear relationship between the intrinsic ORR activity and magnetic moment of 3dxz + 3dyz, suggesting magnetic moment as an ideal activity descriptor. Based on the descriptor, we screened four favorable structures with high activity, namely FeSAC–N1C, FeSAC–N3C, FeSAC–N4C, and FeSAC–N16C. It should be noted that the calculated average Fe–N bond length of them is about 1.90 ± 0.01 Å, which is consistent with the experimental results of EXAFS from XAS (Fig. S21 and Table S7, ESI). We then used the four most possible structures to fully investigate the overall ORR pathway on FeSAC–NC. Fig. 5(e) shows the Gibbs free energy variation diagram of the reaction under U = 1.23 V vs. RHE (Tables S8 and S9, ESI). As seen, there is both high O2 activation barrier (0.70 eV) and *OH desorption barrier (0.63 eV) on FeN4, limiting the ORR activity. Contrarily, barriers for both the *OOH formation and *OH desorption are significantly decreased on FeSAC–NXC, indicating that the spin state manipulation promotes the effective balance between O2 activation and *OH desorption. The overpotentials of FeSAC–N1C, FeSAC–N3C, FeSAC–N4C, and FeSAC–N16C are determined to be 0.39, 0.39, 0.37, and 0.38 V, respectively. Overall, the N[double bond, length as m-dash]C–N configuration-induced medium spin Fe3+ altered Fe 3dz2, 3dxz, 3dyz hybridization with the oxygen-containing intermediates π* orbitals, and better-balanced O2 activation and *OH desorption, thus improving the catalytic performance of ORR. Furthermore, the magnetic moment of 3dxz + 3dyz can be considered as a universal descriptor to predict the ORR activity for Fe–N–C catalysts even with eg = 1 in their metal centers.

4. Conclusions

In summary, we developed a facile strategy for the synthesis of a single-atom Fe(III) N4-dominated N-rich carbon environment, which shows great promise as an oxygen electrode in renewable energy conversion technologies. The formed medium spin Fe(III)–N4 sites induced by innovative N[double bond, length as m-dash]C–N configuration endows Fe(III) with favorable 3dz2, 3dxz (3dyz) hybridization with the oxygen-containing intermediates π* orbitals, which is advantageous to balance O2 activation and *OH desorption. Furthermore, the ORR activity of Fe–N–C catalysts is highly correlated with the magnetic moment of 3dxz + 3dyz, and a smaller magnetic moment of 3dxz + 3dyz yields a higher ORR activity of the Fe–N–C even if their metal centers are all with e1g. The obtained FeSAC–NC exhibits extraordinary ORR catalytic performance in alkaline aqueous electrolytes, even surpassing the commercial Pt/C catalyst. In a homemade zinc–air battery, the FeSAC–NC electrode delivered a slightly higher maximum power density and discharge-specific capacity. Our finding offers new insight into the M–N–C ORR activity origin and provides a guideline to develop highly efficient noble metal-free catalysts.

Author contributions

J.-J. Z., and Z.-F. H. designed the studies. F. L. and C. S. synthesized the catalysts, performed the catalytic tests, and conducted STEM, XPS, Mössbauer spectrum, and XAFS measurements. F. L. performed the density functional theory calculations. F. L. wrote the paper with the help of L. P. and X. Z. All authors discussed the results and commented on the manuscript.

Conflicts of interest

The authors declare no conflicts of interest.

Acknowledgements

The authors appreciate the support from the National Key R&D Program of China (2020YFA0710000), the National Natural Science Foundation of China (22161142002, 22008170, and 22121004), and the Applied Basic Research Program of Qinghai Province (2023-ZJ-701). The authors acknowledge the assistance of XAS characterization and analyses from 1W1B XAFS at Beijing Synchrotron Radiation, Institute of High Energy Physics, Chinese Academy of Sciences, Beijing, China.

References

  1. R. Gao, J. Wang, Z.-F. Huang, R. Zhang, W. Wang, L. Pan, J. Zhang, W. Zhu, X. Zhang, C. Shi, J. Lim and J.-J. Zou, Nat. Energy, 2021, 6, 614–623 CrossRef CAS .
  2. K. Li, R. Zhang, R. Gao, G.-Q. Shen, L. Pan, Y. Yao, K. Yu, X. Zhang and J.-J. Zou, Appl. Catal., B, 2019, 244, 536–545 CrossRef CAS .
  3. M. Luo, Z. Zhao, Y. Zhang, Y. Sun, Y. Xing, F. Lv, Y. Yang, X. Zhang, S. Hwang, Y. Qin, J.-Y. Ma, F. Lin, D. Su, G. Lu and S. Guo, Nature, 2019, 574, 81–85 CrossRef CAS PubMed .
  4. L. Chong, J. Wen, J. Kubal, F. G. Sen, J. Zou, J. Greeley, M. Chan, H. Barkholtz, W. Ding and D.-J. Liu, Science, 2018, 362, 1276–1281 CrossRef CAS PubMed .
  5. H. T. Chung, D. A. Cullen, D. Higgins, B. T. Sneed, E. F. Holby, K. L. More and P. Zelenay, Science, 2017, 357, 479–483 CrossRef CAS PubMed .
  6. X. Ge, A. Sumboja, D. Wuu, T. An, B. Li, F. W. T. Goh, T. S. A. Hor, Y. Zong and Z. Liu, ACS Catal., 2015, 5, 4643–4667 CrossRef CAS .
  7. J. Woo, J. S. Lim, T. Lim, D. S. Baek, J. H. Kim, J. H. Lee, H. Y. Jeong, C. H. Choi and S. H. Joo, EES. Catal., 2023, 1, 62–73 RSC .
  8. N. Zhang, T. Zhou, M. Chen, H. Feng, R. Yuan, C. A. Zhong, W. Yan, Y. Tian, X. Wu, W. Chu, C. Wu and Y. Xie, Energy Environ. Sci., 2020, 13, 111–118 RSC .
  9. C. Tang, L. Chen, H. Li, L. Li, Y. Jiao, Y. Zheng, H. Xu, K. Davey and S. Z. Qiao, J. Am. Chem. Soc., 2021, 143, 7819–7827 CrossRef CAS PubMed .
  10. X. Wan, X. Liu, Y. Li, R. Yu, L. Zheng, W. Yan, H. Wang, M. Xu and J. Shui, Nat. Catal., 2019, 2, 259–268 CrossRef CAS .
  11. M. Zhu, C. Zhao, X. Liu, X. Wang, F. Zhou, J. Wang, Y. Hu, Y. Zhao, T. Yao, L.-M. Yang and Y. Wu, ACS Catal., 2021, 11, 3923–3929 CrossRef CAS .
  12. Y. Wang, Y.-J. Tang and K. Zhou, J. Am. Chem. Soc., 2019, 141, 14115–14119 CrossRef CAS PubMed .
  13. X. T. Wang, T. Ouyang, L. Wang, J. H. Zhong, T. Ma and Z. Q. Liu, Angew. Chem., Int. Ed., 2019, 58, 13291–13296 CrossRef CAS PubMed .
  14. W. Zhong, Y. Qiu, H. Shen, X. Wang, J. Yuan, C. Jia, S. Bi and J. Jiang, J. Am. Chem. Soc., 2021, 143, 4405–4413 CrossRef CAS PubMed .
  15. W. Cheng, P. Yuan, Z. Lv, Y. Guo, Y. Qiao, X. Xue, X. Liu, W. Bai, K. Wang, Q. Xu and J. Zhang, Appl. Catal., B, 2020, 260, 118198 CrossRef CAS .
  16. Z. Chen, H. Niu, J. Ding, H. Liu, P.-H. Chen, Y.-H. Lu, Y.-R. Lu, W. Zuo, L. Han, Y. Guo, S.-F. Hung and Y. Zhai, Angew. Chem., Int. Ed., 2021, 60, 25404–25410 CrossRef CAS PubMed .
  17. J. Yan, Y. Wang, Y. Zhang, S. Xia, J. Yu and B. Ding, Adv. Mater., 2021, 33, e2007525 CrossRef PubMed .
  18. G. Yang, J. Zhu, P. Yuan, Y. Hu, G. Qu, B.-A. Lu, X. Xue, H. Yin, W. Cheng, J. Cheng, W. Xu, J. Li, J. Hu, S. Mu and J.-N. Zhang, Nat. Commun., 2021, 12, 1734 CrossRef CAS PubMed .
  19. J. Suntivich, K. J. May, H. A. Gasteiger, J. B. Goodenough and Y. Shao-Horn, Science, 2011, 334, 1383–1385 CrossRef CAS PubMed .
  20. W. Liu, L. Zhang, X. Liu, X. Liu, X. Yang, S. Miao, W. Wang, A. Wang and T. Zhang, J. Am. Chem. Soc., 2017, 139, 10790–10798 CrossRef CAS PubMed .
  21. R. Nityananda, J. Sci. Educ., 2017, 22, 809 Search PubMed .
  22. X. Li, X. Yang, L. Liu, H. Zhao, Y. Li, H. Zhu, Y. Chen, S. Guo, Y. Liu, Q. Tan and G. Wu, ACS Catal., 2021, 11, 7450–7459 CrossRef CAS .
  23. M. Xiao, Z. Xing, Z. Jin, C. Liu, J. Ge, J. Zhu, Y. Wang, X. Zhao and Z. Chen, Adv. Mater., 2020, 32, 2004900 CrossRef CAS PubMed .
  24. J. Zhao, J. Chen, S. Xu, M. Shao, Q. Zhang, F. Wei, J. Ma, M. Wei, D. G. Evans and X. Duan, Adv. Funct. Mater., 2014, 24, 2938–2946 CrossRef CAS .
  25. H. Zhang, W. Cheng, D. Luan and X. W. Lou, Angew. Chem., Int. Ed., 2021, 60, 13177–13196 CrossRef CAS PubMed .
  26. Y. Zheng, Y. Jiao, Y. Zhu, L. H. Li, Y. Han, Y. Chen, A. Du, M. Jaroniec and S. Z. Qiao, Nat. Commun., 2014, 5, 3783 CrossRef PubMed .
  27. S. M. Kirtley, O. C. Mullins, J. Chen, J. Vanelp, S. J. George, C. T. Chen, T. Ohalloran and S. P. Cramer, Biochim. Biophys. Acta, 1992, 1132, 249–254 CrossRef CAS PubMed .
  28. K. Fujii, K. Akamatsu and A. Yokoya, J. Phys. Chem. B, 2004, 108, 8031–8035 CrossRef CAS .
  29. Y. Luo, Y. Chen, Y. Xue, J. Chen, G. Wang, R. Wang, M. Yu and J. Zhang, Small, 2022, 18, 2105594 CrossRef CAS PubMed .
  30. L. Lin, Z. K. Yang, Y.-F. Jiang and A.-W. Xu, ACS Catal., 2016, 6, 4449–4454 CrossRef CAS .
  31. J. Wan, Z. Zhao, H. Shang, B. Peng, W. Chen, J. Pei, L. Zheng, J. Dong, R. Cao, R. Sarangi, Z. Jiang, D. Zhou, Z. Zhuang, J. Zhang, D. Wang and Y. Li, J. Am. Chem. Soc., 2020, 142, 8431–8439 CrossRef CAS PubMed .
  32. J. Qin, H. Liu, P. Zou, R. Zhang, C. Wang and H. L. Xin, J. Am. Chem. Soc., 2022, 144, 2197–2207 CrossRef CAS PubMed .
  33. J. Feng, H. Gao, L. Zheng, Z. Chen, S. Zeng, C. Jiang, H. Dong, L. Liu, S. Zhang and X. Zhang, Nat. Commun., 2020, 11, 4341 CrossRef PubMed .
  34. U. I. Koslowski, I. Abs-Wurmbach, S. Fiechter and P. Bogdanoff, J. Phys. Chem. C, 2011, 115, 23417–23427 CrossRef .
  35. U. I. Koslowski, I. Abs-Wurmbach, S. Fiechter and P. Bogdanoff, J. Phys. Chem. C, 2008, 112, 15356–15366 CrossRef CAS .
  36. U. I. Kramm, M. Lefevre, N. Larouche, D. Schmeisser and J. P. Dodelet, J. Am. Chem. Soc., 2014, 136, 978–985 CrossRef CAS PubMed .
  37. A. Zitolo, V. Goellner, V. Armel, M. T. Sougrati, T. Mineva, L. Stievano, E. Fonda and F. Jaouen, Nat. Mater., 2015, 14, 937–942 CrossRef CAS PubMed .
  38. Y. Liu, X. Liu, Z. Lv, R. Liu, L. Li, J. Wang, W. Yang, X. Jiang, X. Feng and B. Wang, Angew. Chem., Int. Ed., 2022, 61, e202117617 CAS .
  39. Q. Cheng, L. Yang, L. Zou, Z. Zou, C. Chen, Z. Hu and H. Yang, ACS Catal., 2017, 7, 6864–6871 CrossRef CAS .
  40. E. Luo, H. Zhang, X. Wang, L. Gao, L. Gong, T. Zhao, Z. Jin, J. Ge, Z. Jiang, C. Liu and W. Xing, Angew. Chem., Int. Ed., 2019, 58, 12469–12475 CrossRef CAS PubMed .
  41. H.-x Zhong, J. Wang, Y.-w Zhang, W.-l Xu, W. Xing, D. Xu, Y.-f Zhang and X.-b Zhang, Angew. Chem., Int. Ed., 2014, 53, 14235–14239 CrossRef CAS PubMed .
  42. Z. Fu, B. Yang and R. Wu, Phys. Rev. Lett., 2020, 125, 156001 CrossRef CAS PubMed .

Footnote

Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d3ey00066d

This journal is © The Royal Society of Chemistry 2023
Click here to see how this site uses Cookies. View our privacy policy here.