Yongheng
Xiong†
,
Yang
Yu†
,
Huangang
Shi
,
Jifa
Qu
and
Wenyi
Tan
*
Department of Environmental Engineering, Nanjing Institute of Technology, Nanjing 211167, PR China. E-mail: twy1102@gmail.com; Tel: +86 025-86118922
First published on 21st November 2025
Single atom Fe sites, doped in CN materials, exhibit outstanding electrochemical activity for CO2-to-CO conversion. The pyrolysis of ZIF8 is a controllable method for fabricating isolated single atom metal sites. In this study, we propose a new strategy to increase the ratio of Fe in ZIF8 precursors by synergistically replacing Zn2+ with Bi3+ and Fe3+. After precursor pyrolysis, the obtained Bi/Fe–N–C catalysts, consisting of Bi sites and pyrrole-type Fe–Nx sites, serve as efficient electrocatalysts for the CO2RR. The results show that the optimized catalyst loaded with 94.8 mg per kgCat Fe exhibits a high FECO of >90.1% over a wide potential range of −0.4 to −0.7 VRHE (98.2% at −0.5 VRHE). Insights into the electrochemical reaction mechanism show that this successful design of Bi/Fe–N–C catalysts can provide a stable catalytic site to form *COOH, thus achieving energy-efficient electrochemical CO2 reduction to CO.
Carbon-based materials have received extensive attention for the CO2RR due to their porous structure and controllable doping modification.15,16 These materials have distinct intrinsic properties, such as abundant natural resources, environmental friendliness, customized structures, high surface area, tailored structures, high-temperature stability, superior electric conductivity, and chemical inertness to acids and bases.17 Despite the low CO2RR activity of pure carbon-based catalysts, heteroatom doping can improve the current efficiency and conductivity of carbon-based materials owing to the strong electron connection between the heteroatom and carbon.18 The N-doped carbon-based materials outperform other doping atoms (S, B, F, and P) in terms of overall performance.19,20 Based on the reported results,21 the graphitic and pyridinic N of N-doped carbon-based materials play a key role in increasing the selectivity and decreasing the overpotential towards CO2-to-CO conversion. Meanwhile, calculations by density functional theory (DFT) confirmed that pyridinic N atoms with a lone pair of electrons can bind to CO2.22,23 Despite multiple promising results, the fatal flaw of metal-free catalysts is their poor catalytic activity, which limits further research on the CO2RR process.
Porous carbon materials doped with N-coordinated metal sites (M–N–C), usually derived from pyrolyzing the precursors composed of metal-(M-), N- and C-containing complexes, are derivative electrocatalysts of N-doped carbon for the CO2 reduction reaction (CO2RR).24–26 With remarkable catalytic activity and excellent selectivity, M/N-doped carbon has gained a lot of attention, especially for Fe–N–C catalysts.27 Many studies have confirmed that the effective method for improving catalytic performance towards the CO2RR is to achieve high exposure of Fe–Nx sites on Fe–N–C catalysts.21,28–30 For this purpose, Fe-doped ZIF8 is widely selected as a precursor to design Fe–N–C catalysts.31 Much effort has been put into precisely designing ZIF8 precursor with high surface areas and rich micropores, such as surface functionalization with ammonium ferric citrate,30 post-synthetic strategy,32,33 and carrier modification.34 However, the framework structure of ZIF8 determines that Fe3+ can only partially replace the Zn2+ center site, limiting the molar ratio of Fe3+/Zn2+.35 Few studies have been reported to dramatically increase the Fe ratio in the ZIF8 precursor towards Fe–N–C catalysts. In pursuit of better CO2RR performance, it is important to develop a method to increase the upper limit of the iron content in ZIF8 precursors.
In this paper, we report a chemical assisting approach to synthesize Bi/Fe-doped ZIF8 precursors. The Bi element is explored to increase the Fe ratio in ZIF8. One-step thermal activation at 950 °C is applied to successfully prepare the Fe–N–C catalysts with isolated Fe–Nx sites anchored on N-doped carbon. The Fe–N–C catalyst with the highest CO2RR activity was screened by regulating the Bi content in the mother liquor of the ZIF8 precursor. The Bi4/Fe–N–C catalyst with the loaded 94.8 mg per kgCat Fe exhibits a high FECO above 90% over a wide potential range from −0.4 to −0.7 VRHE (98% at −0.5 VRHE). A series of experimental measurements and in situ characterization reveal that a high Fe ratio in Bi/Fe–N–C is conducive to promoting the dissociation of H2O, accelerating the proton transfer process and forming *COOH intermediates.
![]() | ||
| Fig. 1 Illustration of the Bi-doping chemistry strategy to fabricate Bi/Fe-doped ZIF8 precursors for catalyst preparation. | ||
According to the preparation process of Bi/Fe-ZIF8 precursors, the mass ratios of Bi(NO3)3·5H2O and Fe(NO3)3·9H2O were configured to 1/8 (10 mg/80 mg), 2/8 (20 mg/80 mg), 3/8 (30 mg/80 mg), 4/8 (40 mg/80 mg) and 5/8 (50 mg/80 mg). After thermal activation, the obtained catalysts correspond to Bi1/Fe–N–C, Bi2/Fe–N–C, Bi3/Fe–N–C, Bi4/Fe–N–C and Bi5/Fe–N–C, respectively.
During the CO2 reduction experiments, the electrolyte in the cathodic compartment was stirred at 900 rpm. Linear sweep voltammetry (LSV) was performed with a scan rate of 10 mV s−1 from −0.6 V to −1.8 V (vs. Ag/AgCl) in Ar-saturated 0.5 M KHCO3 (pH = 8.8) and CO2-saturated 0.5 M KHCO3 (pH = 7.2) as the supporting electrolyte. The potentials in the study were reported versus RHE with the conversion E (vs. RHE) = E (vs. Ag/AgCl) + 0.1989 V + 0.0591 × pH.
The current density was obtained by normalizing it with the carbon fiber paper geometric surface area. CO2 gas was delivered at an average rate of 20 mL min−1 (at room temperature and ambient pressure) and routed directly into the gas sampling loop (1 mL) of a gas chromatograph. The gas phase composition was analyzed by GC every 30 min. The separated gas products were analyzed using a thermal conductivity detector (TCD) and a flame ionization detector (FID). The quantification of the products was performed using the conversion factor derived from the standard calibration gases. Liquid products were analyzed afterwards using quantitative NMR with dimethyl sulphoxide (DMSO) as an internal standard. The solvent pre-saturation technique was implemented to suppress the water peak.
485C mol−1; and FECO: faradaic efficiency for CO production.
485C mol−1; mcat: catalyst mass in the electrode, g; w: Fe loading in the catalyst from ICP; and Fe: atomic mass of iron, 56 g mol−1.
![]() | ||
| Fig. 2 (a–e) SEM images of Bi1/Fe-ZIF8, Bi2/Fe-ZIF8, Bi3/Fe-ZIF8, Bi4/Fe-ZIF8 and Bi5/Fe-ZIF8. (f) Ratios of Zn, Bi, and Fe to N elements based on the EDS results. | ||
After 2 hours of annealing at 950 °C under Ar protection, the N-coordinated Fe sites of Bi/Fe–N–C were successfully anchored on porous carbon, serving for the CO2RR process as catalytically active sites. Meanwhile, part of the Zn was reduced by high-temperature carbonization and evaporated away (b.p. 907 °C).26 For comparison, Bi–N–C and Fe–N–C were also treated with the same annealing process. As revealed by scanning electron microscopy (SEM) and transmission electron microscopy (TEM) (Fig. S8), the as-prepared Fe–N–C presents a uniform size (450–500 nm) and retains the polyhedral morphology of the reported ZIF-8.31,35 As shown by EDS mapping images, all key elements (C, N, and Fe) are evenly distributed in the Fe–N–C samples. Importantly, no obvious metal nanoparticles are observed through the above characterization. Based on previous studies,28–30 it could be concluded that the Fe element is anchored on porous carbon in the form of Fe–Nx.
However, the collapsed edge and dispersive nanoparticles of Bi4/Fe–N–C are clearly presented in the SEM and TEM (Fig. 3a and b). The morphology of the nanoparticles displays a core–shell nanostructure covered with graphite carbon (Fig. S23). The other major part, without agglomeration, remains the structure of porous carbon. This is also consistent with the results of XRD, as depicted in Fig. S9. The Bi4/Fe–N–C pattern has no characteristic diffraction peaks, which are similar to Fe–N–C and Bi–N–C. Two broad diffraction peaks are located at 25.3° and 43.7°, which correspond to graphitic carbon.37,38 There is no other obvious iron-containing phase in Bi4/Fe–N–C catalysts. The analysis of EDS shows the distributions of N, C, Bi and Fe, which directly proves that the aggregation of Fe results in the formation of the Fe nanoparticles (NPs) (Fig. S23(b–f)). Other bulk parts of the Bi4/Fe–N–C catalyst without Fe NPs also show a uniform distribution of Fe. The atomically dispersed Fe sites were evidenced by the aberration-corrected high-angle annular dark-field scanning TEM (HAADF-STEM), showing numerous bright dots on the low-atomic number hosts (Fig. 3c and d). C and N, as the two dominating elements, are uniformly distributed on the carbon matrix (Fig. 3f and h), indicating the successful synthesis of porous N-doped carbon. Compared with the EDS mapping results of Fe–N–C (Fig. S8c–f), the Fe region of Bi4/Fe–N–C is similar to that of C and N (Fig. 3e), indicating the different content of Fe in different catalysts.
![]() | ||
| Fig. 3 (a) SEM, (b) TEM and (c and d) HAADF-STEM image of Bi4/Fe–N–C and the corresponding EDS mappings of Fe (e), C (f), Bi (g), N (h) and Zn (i). | ||
The elemental information of Bi–N–C, Fe–N–C and Bi4/Fe–N–C surfaces is obtained from XPS (Fig. 4a). It is evident that the main difference is the appearance of the Fe signal and the disappearance of the Zn signal for the Bi4/Fe–N–C pattern. The details of the Fe 2p pattern (Fe 2p3/2, 711.9 eV; Fe 2p3/2, 725.2 eV) are displayed clearly in Fig. 4b. The Bi4/Fe–N–C catalyst shows a peak of Fe 2p, as opposed to Bi–N–C and Fe–N–C. This demonstrates that the Fe content of the Bi4/Fe–N–C surface is significantly higher than those of Bi–N–C and Fe–N–C. Other support from the ICP results also confirms the increase in Fe content in the Bi4/Fe–N–C catalyst (Fig. 4f). Thus, it also proves that the Bi-doping method can effectively increase the content of Fe in the ZIF8 precursor. In contrast, the high-resolution XPS Zn 2p spectra show an opposite result to the above (Fig. 4c). At 1021.9 eV, the Zn 2p signal of Bi4/Fe–N–C is smooth without any peak. Combined with the ICP results of the Zn element, the Zn content decreases remarkably in the Bi4/Fe–N–C catalyst, which may result from the formation of Fe NPs.39 The agglomeration of Fe destroys the N–C structure, resulting in more Zn being reduced and then volatilized. As Fe metal evaporates at a higher temperature than Zn, Fe–Nx sites are more easily anchored on carbon than Zn–Nx sites. Thus, the residual amount of Zn (Zn–Nx) is affected by the Fe-doping content.
As shown in Fig. 4d, there are no characteristic peaks in the Bi 4f XPS spectrum. The details of the chemical composition are confirmed by ICP (Fig. 4f). Therefore, the Bi element can be doped into ZIF8 but only in extremely small amounts. In spite of this, it is key for Fe-doping into ZIF8. Moreover, the high-resolution XPS N 1s spectra of all the samples were deconvoluted into pyridinic (398.5 eV), pyrrolic (400.9 eV), graphitic (402.0 eV), and oxidized (404.0 eV) N species (Fig. 4e).26,28,40 The porphyrin-like moieties at 399.7 eV correspond to metal-nitrogen (M–N) coordination. Visibly, the M–N peak area of Bi4/Fe–N–C is higher than those of Fe–N–C and Bi–N–C. Since the major metal increase is the Fe element after Bi-doping in Fe-ZIF8, the main part of the M–N site in the Fe–N–C catalyst should be the Fe–Nx site. It is generally evidenced to be an active site for the CO2RR. This reflects the superiority of the Bi-doping strategy for the synthesis of abundant Fe–Nx sites.
The CO2 electroreduction activities of Bi4/Fe–N–C, Fe–N–C and Bi–N–C were tested using a three-electrode H-cell with a 0.5 M KHCO3 solution as the electrolyte. As revealed by linear sweep voltammetry (LSV) shown in Fig. 5a, Bi4/Fe–N–C exhibits the best CO2RR activity with a current density of 8.55 mA cm−2 at −0.60 VRHE in CO2-saturated electrolyte, which is roughly 1.3 and 6.6 times as high as Fe–N–C and Bi–N–C. Besides, in an Ar-saturated KHCO3 electrolyte, Bi4/Fe–N–C still presents a very similar trend of CO2RR activity with a relatively small current density, indicating that the remarkable activity of Bi4/Fe–N–C was saturated with CO2. The gas products of carbon dioxide electroreduction analyzed by on-line gas chromatography (GC) and off-line nuclear magnetic resonance (NMR) were identified as only CO and H2. The CO faradaic efficiencies (FECO) of Bi4/Fe–N–C were compared in potentials ranging from −0.4 to −0.7 VRHE (Fig. 5c). The maximum FECO was achieved at approximately 98% at −0.5 and −0.55 VRHE, while the FECO slowly decreased at a higher potential matched with the current–time response of Bi4/Fe–N–C (Fig. 5b).
![]() | ||
| Fig. 5 (a) LSV curves of Bi–N–C, Fe–N–C and Bi4/Fe–N–C in 0.5 M KHCO3 solution. (b) Stability test of Bi4/Fe–N–C. (c) vco and FECO. (d) TOF of Bi4/Fe–N–C. | ||
As shown in Fig. 5d, the calculated TOFs for Bi4/Fe–N–C and Fe–N–C present a clear distinction at different potentials. At −0.5 VRHE, the TOF of Bi4/Fe–N–C is calculated to be 468.6 h−1, which is 8.1 times higher than that of Fe–N–C. Moreover, the disparity of TOF between Bi4/Fe–N–C and Fe–N–C gradually increases with the increase in potential. It also proves that an increase in Fe–Nx density promotes the activity of CO2 electroreduction.12–14
To accurately understand the effect of Bi on the Bi/Fe–N–C catalyst, different ratios of Bi/Fe were designed to prepare precursors of Bi/Fe–N–C. After annealing activation, the pyrolyzed M–N–C materials with various Fe contents were characterized by XRD, as illustrated in Fig. S15, and showed no significant differences. The ICP results (Fig. S16) reflect the opposite two-sided effect of Bi on Fe content, suppressing at a low Bi ratio and improving under a high Bi ratio. The CO2 electroreduction activities present similar results (Fig. 6a). CO yield (activity) and FECO (selectivity) improve with increasing Fe content. However, the destruction of ZIF8 morphology and the agglomeration of Fe are inevitable and increase during the pyrolysis of precursors with a higher Bi ratio, which can greatly reduce the formation of active sites (Fe–Nx) for CO2-to-CO conversion. Thus, the activity of the electrochemical CO2RR is significantly reduced at a high Bi ratio.
![]() | ||
| Fig. 6 (a) vco and FECO of Bi/Fe–N–C at −0.5 VRHE. (b) LSV difference of the samples under different saturated atmospheres (LSV(CO2)–LSV(N2)). | ||
To compare LSV in Ar- and CO2-saturated electrolyte, the difference of value is calculated using the following formula: δ = LSV(CO2) − LSV(N2). As shown in Fig. 6b, similar results tend to coincide except for Bi–N–C. The results indicate that the current density in Ar-saturated electrolyte increases with iron content.41 Herein, the hydrogen evolution reaction (HER) is the main competitive reaction of the CO2RR.24,42 However, in a CO2-saturated electrolyte, the Fe–Nx sites can provide an active site for the second CO2 activation and decrease the barrier of CO2(g) → CO(g).43 According to the reported DFT calculations, CO2 absorbed in the electrolyte can be adsorbed and reduced on the active site (Fe–Nx) for reaction.26,44 This means that the enhancement of CO2 activation results in high activity and FECO for electrochemical CO2RR on Fe–N sites. Thus, a similar value of δ proves that a high iron content in Bi4/Fe–N–C can markedly improve electrochemical CO2RR and FECO in a CO2-saturated electrolyte.
Furthermore, in situ ATR-FTIR spectroscopy was performed under CO2RR conditions to identify and monitor the reaction intermediates and provide mechanistic insights into CO production. It was conducted using a Thermo-Fisher Nicolet iS20 system equipped with a liquid nitrogen-cooled HgCdTe (MCT) detector and a PIKE electrochemical three-electrode cell (Fig. S17). According to the activity of electrochemical CO2RR, −0.5 VRHE was selected as the operating voltage. During the reaction process, the ATR-FTIR spectra for Bi4/Fe–N–C catalysts show several peaks (Fig. 7). The peaks at 1490 and 1060 cm−1 correspond to the atop and bridge configurations of *COOH adsorption, respectively, which increase with time.45 This consequence indicates that *COOH can be stably adsorbed on the catalyst surface during the electrochemical process of CO2(g) → CO(g). The increased peaks of CO and CO2 further support the above speculation. For the CO2RR, the pathway of CO2 protonation first adsorbs CO2 on the surface of the Fe–N–C catalyst, and the adsorbed CO2 converts to *COOH. Next, *COOH is broken to form CO.
Supplementary information is available. See DOI: https://doi.org/10.1039/d5se01074h.
Footnote |
| † These authors contributed equally to this work. |
| This journal is © The Royal Society of Chemistry 2026 |