A pyrolysis-free path toward superiorly catalytic nitrogen-coordinated single atom

See allHide authors and affiliations

Science Advances  02 Aug 2019:
Vol. 5, no. 8, eaaw2322
DOI: 10.1126/sciadv.aaw2322


Nitrogen-coordinated single-atom catalysts (SACs) have emerged as a frontier for electrocatalysis (such as oxygen reduction) with maximized atom utilization and highly catalytic activity. The precise design and operable synthesis of SACs are vital for practical applications but remain challenging because the commonly used high-temperature treatments always result in unpredictable structural changes and randomly created single atoms. Here, we develop a pyrolysis-free synthetic approach to prepare SACs with a high electrocatalytic activity using a fully π-conjugated iron phthalocyanine (FePc)–rich covalent organic framework (COF). Instead of randomly creating Fe-nitrogen moieties on a carbon matrix (Fe-N-C) through pyrolysis, we rivet the atomically well-designed Fe-N-C centers via intermolecular interactions between the COF network and the graphene matrix. The as-synthesized catalysts demonstrate exceptional kinetic current density in oxygen reduction catalysis (four times higher than the benchmark Pt/C) and superior power density and cycling stability in Zn-air batteries compared with Pt/C as air electrodes.


Carbon-based catalysts with nitrogen-coordinated single transition metal atoms as active sites have emerged as promising systems for catalytic processes because of their maximal atom utilization and highly catalytic activity (18). As revealed both experimentally and theoretically, the nitrogen-coordinated single transition metal moieties supported in carbon provide versatile active sites, while the carbon matrix ensures a large specific area and efficient mass transport, provides a stable matrix for the metal atoms, and affects the electronic density owing to strong intermolecular interactions (916). Currently, reported strategies to prepare the N-coordinated single-atom catalysts (SACs) dominantly depend on high-temperature pyrolysis, which highly demands the rigorous manipulation of the carbonization process and precise elemental ratio of precursors (1721). For instance, during the pyrolysis of zeolitic imidazolate frameworks to prepare catalysts with single Fe-N-C sites, slight changes in pyrolysis temperature will lead to notable differences in catalytic performance (2225). Meanwhile, because the isolated metal atoms are thermodynamically unstable, for many situations, the metal atoms tend to agglomerate as their ratio in precursors increases but hardly form single-atom centers when the ratio decreases (14, 26). By this means, randomly creating Fe-N-C sites on the carbon matrix through pyrolysis is unpredictable and difficult to repeat in an atomic manner. There are many researches working on precisely tuning the nitrogen-coordinated SAC (2528). However, developing precisely controllable synthesis of nitrogen-coordinated SACs in the long term for practical applications by means of high-temperature pyrolysis remains a challenge. To this end, we developed a pyrolysis-free synthetic approach SACs with a high electrocatalytic activity using a fully closed π-conjugated iron phthalocyanine (FePc)–rich covalent organic framework (COF).


Instead of randomly creating Fe-nitrogen moieties on the carbon matrix (Fe-N-C) through pyrolysis, we rivet the atomically well-designed Fe-N-C centers via intermolecular interactions between the Fe-N-C–rich pyrolysis-free COF network and the graphene matrix. Typically, as shown in Fig. 1A, we assembled benzene-1,2,4,5-tetracarbonitrile with Fe centers into a fully closed π-conjugated COF (termed COFBTC) according to our previous method (29) to fabricate atomically dispersed Fe-N-C moieties (for additional details, see Materials and Methods), which were subsequently assembled with the graphene matrix via intermolecular interactions [see the extended x-ray absorption fine structure (EXAFS) and resistivity measurements discussed below]. Different from unpredictable pyrolysis, our single Fe-N-C centers were pre-assembled and directly riveted onto the graphene matrix; thus, the structures of the pyrolysis-free synthesized SACs (marked as pfSAC-Fe-X, where X represents the added mass ratio of COFBTC) were well atomically arranged in the intact COF network. As expected, the graphene provided a stable matrix for the pre-assembled single Fe-N-C centers of COFBTC. On the basis of the transmission electron microscopy (TEM) images (associated with elemental mapping), the as-obtained pfSAC-Fe demonstrated homogeneous elemental distribution (figs. S1 and S2), which was also confirmed by x-ray photoelectron spectroscopy (XPS) elemental analyses (fig. S3). We further performed atomic-resolution high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) and identified that the coordinated Fe atoms, as represented by bright dots, were exclusively in the single-atom format and uniformly anchored throughout the graphene matrix (Fig. 1B and fig. S4).

Fig. 1 Synthesis and structural characterization of pfSAC-Fe.

(A) Synthesis route of the pfSAC-Fe catalyst. The intermolecular interactions are constructed during the process. (B) HAADF-STEM of the pfSAC-Fe-0.2. The graphene matrix is observed, on which plentiful Fe atoms (bright dots) are anchored. (C and D) Fe K-edge XANES spectra and Fourier transform (FT) of pfSAC-Fe-0.2, FePc, and Fe foil. a.u., arbitrary units. (E) Fe K-edge EXAFS analysis result in R spaces. The χ signals are the two-body backscattering paths of Fe-N, Fe-C, and Fe-Fe. (F) Simulated structures according to the Fe K-edge EXAFS analysis result.

The riveting of single Fe-C-N centers was successfully realized via the intermolecular interactions between the COFBTC and the graphene matrix. As investigated by x-ray absorption fine structure (XAFS) spectroscopy, the Fe K-edge of x-ray absorption near-edge structure (XANES) of these specimens suggested that the pfSAC-Fe contained near-edge structures similar to those of the original FePc but totally different from those of the Fe foil (Fig. 1C), indicating a N-coordinated chemical state of single Fe atoms. In the Fourier transform analyses, the curve of pfSAC-Fe demonstrated a main peak belonging to the Fe-N scattering paths, while the Fe-Fe peak at about 2.2 Å was not detected (Fig. 1D). We further simulated the Fe K-edge EXAFS, finding the signals of the two-body backscattering path of Fe-N in both pfSAC-Fe and COFBTC while Fe-C was uniquely formed in the pfSAC-Fe (Fig. 1E and figs. S5 to S7). Accordingly, only van der Waals interaction was observed during the electron localization function (ELF) analysis (Fig. 2A) with the density functional theory (DFT) calculations. As confirmed by the comparison of charge density differences, the electrons of graphene were attracted to the N-coordinated Fe sites, forming an Fe-C electron pathway (Fig. 2B). Moreover, once oxygen was introduced, we observed that it was only absorbed onto the single Fe atoms (Fig. 2C), which offered superiorly catalytic sites for the oxygen reduction reaction (ORR) (vide infra).

Fig. 2 Atomic structure analysis of pfSAC-Fe.

(A) ELF of pfSAC-Fe. (B) Differential charge density distribution on pfSAC-Fe. (C) Differential charge density distribution on pfSAC-Fe with absorption of oxygen.

Similar to randomly creating Fe-N-C centers on the carbon matrix through pyrolysis, our approach to rivet the pre-assembled Fe-N-C onto graphene also generated a facile electron transfer environment. In the powder electrical resistivity measurement, we observed that once the intermolecular interactions were constructed, the resistivity of as-prepared pfSAC-Fe-0.2 decreased markedly, becoming very close to that of the raw graphene (Fig. 3A). By constructing the equivalent circuit based on the experimental data, we noted that the resistivity of the pfSAC-Fe-0.2 was in between the parallel connection and the series-parallel connection at low pressure (<12 MPa) while climbing higher than either the parallel connection or the series-parallel connection as the pressure increases (Fig. 3B). Because of the high pressing force leading to a much tighter connection and better conductivity, this phenomenon suggested that the intermolecular interactions generated an electrically conducting bridge between the Fe-N-C moieties and the graphene matrix, which was not as efficient as the conjugated systems indwelling the graphene, but still greatly benefitted electron transport. As observed from the density of state (DOS), the COFBTC was a semiconductor with a small direct bandgap. With the intermolecular interactions with graphene, the pfSAC-Fe exhibited excellent electrical conductivity similar to metal, which prominently enhanced the charge transfer (fig. S12). Furthermore, the work function of pfSAC-Fe-0.2 obtained from the ultraviolet photoelectron spectroscopy (UPS) spectra in the valence band emission region was as small as 4.60 eV (Fig. 3C), suggesting a high driving force for donating electrons from the catalytic surface and facilitating the catalytic performance during ORR (30, 31).

Fig. 3 Electronic conductivity studies.

(A) The resistivity of COFBTC, pfSAC-Fe-0.2, and pure graphene at different pressures. (B) The calculated equivalent circuit as series, parallel, and series&parallel connections compared with the pfSAC-Fe-0.2. Below: Schematic structures for series, parallel, and series&parallel connections. (C) The UPS spectra in the valence band emission region for pfSAC-Fe-0.2. The work function was calculated as 4.60 eV.

The electrocatalytic performance of pfSAC-Fe catalysts with different added ratios of COFBTC was evaluated with a traditional three-electrode system. In cyclic voltammetry (CV) and linear sweep voltammetry (LSV) measurements with a rotating disk electrode, all of the pfSAC-Fe catalysts exhibited excellent catalytic performance for ORR compared to the benchmark Pt/C (Fig. 4A). In addition, the performance of pfSAC-Fe was enhanced with the increase in the added ratio of COFBTC, demonstrating considerable high half-wave (E1/2) potentials from 0.88 to 0.910 V versus reversible hydrogen electrode (RHE) (Fig. 4B). Meanwhile, the E1/2 of the benchmark Pt/C was 0.86 V versus RHE, demonstrating the same activity as that in state-of-the-art literature (3234). Consistent with the DFT calculations, the single Fe atoms offered superior catalytic sites for the absorption and electron transfer of oxygen molecules during ORR. The optimized catalyst, pfSAC-Fe-0.2, exhibited the best ORR activity with a high kinetic current density (Jk) of 25.86 mA cm−2 at 0.85 V (versus RHE), which was about four times higher than that of Pt/C (6.49 mA cm−2), and with a low Tafel slope of 31.7 mV decade−1 (Fig. 4, B and C). The efficient kinetics were also confirmed by the Koutecky-Levich (K-L) plots and rotating ring disk electrode (RRDE) test, revealing a highly efficient four-electron pathway (3.85 to 4.0 per O2) and a low H2O2 yield (below 7%) over 0.5 to 0.9 V (Fig. 4, D and E) (35). We also noted that the electrocatalytic performance, such as half-wave potential and kinetic current density, was enhanced with increased catalytic loading (fig. S13), which was consistent with pioneering work (36, 37). Superior to commercial Pt/C, the pfSAC-Fe catalysts exhibited better methanol tolerance during the methanol crossover test (Fig. 4F) and much longer durability during the continuous cycles (fig. S14A). The XAFS measurement of the pfSAC-Fe-0.2 test demonstrated similar curves after ORR tests, indicating that the structure of the catalysts was quite stable during the reaction (fig. S15). In addition, TEM and HAADF-STEM analyses confirmed that the single Fe atoms were still uniformly anchored throughout the graphene matrix after the ORR test (figs. S16 and S17).

Fig. 4 Electrochemical characterization.

(A) LSV curves of pfSAC-Fe-X and 20% Pt/C in O2-saturated 0.1 M KOH solution at a scan rate of 5 mV s−1 and a rotation speed of 1600 rpm. (B) Kinetic current density (at 0.85 V versus RHE) and half-wave potentials of pfSAC-Fe-X and 20% Pt/C. (C) Tafel plots of graphene, pfSAC-Fe-0.2, and Pt/C. (D) LSV curves of pfSAC-Fe-0.2 at different rotation speeds. Inset: The corresponding K-L plots and electron transfer number. (E) H2O2 yields and electrode transfer number of pfSAC-Fe-0.2 and Pt/C in O2-saturated 0.1 M KOH. (F) CV curves of pfSAC-Fe-0.2 in O2-saturated 0.1 M KOH without and with 1.0 M methanol solutions.

Furthermore, in the homemade Zn-air battery with pfSAC-Fe-0.2 as the electrocatalysts for air cathode (Fig. 5A), we observed an open circuit potential as high as 1.41 V (Fig. 5B). Compared with the battery using conventional Pt/C as the electrocatalysts, the Zn-air battery driven by pfSAC-Fe-0.2 delivered similar discharge behaviors, higher power density (123.43 mW cm−2, compared to 113.81 mW cm−2, the power density of the battery driven by Pt/C), and larger specific capacity (732 mAh gZn−1 at 100 mA cm−2 with a zinc utilization of 89.3%) at the same loading (0.2 mg cm−2) (Fig. 5, C and D). Besides, the Zn-air battery based on pfSAC-Fe-0.2 exhibited a notably long life cycle over 300 hours with a less than 0.1% decrease, suggesting the considerable stability of the pfSAC-Fe and the practical potential of the as-designed batteries, while the discharge voltage of the Pt/C-based system markedly decreased within 20 hours under the same conditions (Fig. 5E). We further periodically changed the current density from 5 to 40 mA cm−2 to evaluate the stability of cycle discharge performance and only observed a small voltage difference over 30 hours (Fig. 5F). Unexpectedly, even when the current density was as high as 100 mA cm−2, the battery was still able to operate for more than 300 min with little potential drop (Fig. 5G). By increasing the loading of pfSAC-Fe catalysts to 0.5 mg cm−2, the peak power density of the battery was raised up to 126.83 mW cm−2 (fig. S23), demonstrating a promising growth potential to replace the commercial Pt/C catalysts for metal-air batteries with high power density.

Fig. 5 Zn-air battery performance of pfSAC-Fe.

(A) A schematic configuration of the homemade Zn-air battery. (B) Open circuit plot of the Zn-air battery using pfSAC-Fe-0.2 as catalyst. (C) Comparison of polarization and power density curves using Pt/C and pfSAC-Fe-0.2 as catalysts. (D) Comparison of specific capacities of the Zn-air batteries. (E) Long-term discharge/charge cycling performance of Zn-air batteries with Pt/C and pfSAC-Fe-0.2 as the ORR catalyst and IrO2 as the oxygen evolution reaction catalyst at a current density of 5 mA cm−2. Inset: Voltage efficiency (VE) of Zn-air batteries based on pfSAC-Fe-0.2 during the cycling test. (F) Cycle discharge curves of pfSAC-Fe-0.2–driven Zn-air batteries at periodically changed current densities of 5, 10, 20, and 40 mA cm−2. (G) Discharge curves of Zn-air batteries using Pt/C and pfSAC-Fe-0.2 at current densities of 100 mA cm−2. For each battery, a Zn plate was used as the anode, and 8 M KOH + 0.5 M ZnO was used as the electrolyte. All of the catalytic mass loading is at 0.2 mg cm−2.


Rather than randomly creating single atoms on carbon through pyrolysis, we developed a pyrolysis-free synthetic approach to prepare SACs by riveting the pre-assembled COFBTC on a conductive matrix (e.g., graphene) via intermolecular interactions. The as-synthesized pfSAC-Fe demonstrated a small work function and highly efficient catalytic performance for ORR. The Zn-air battery driven by the pfSAC-Fe exhibited superior power density and cycling stability for practical applications. The invented pyrolysis-free approach toward SACs not only paves a promising way to atomically optimize the SAC structures but also simplifies the accurate theoretical prediction calculations based on well-defined structures to develop more excellent catalysts for various energy conversions and storage devices (e.g., CO2 reduction, N2 fixation, and photoelectronic catalysis) besides oxygen reduction.


Synthesis of pfSAC-Fe-X

To synthesize the pfSAC-Fe-X electrocatalysts, different mass ratios of as-synthesized COFBTC (1, 5, 10, and 20%) were dispersed in glycol solution, followed by grinding with corresponding amounts of graphene. After grinding the mixture for 1 hour, the as-obtained composite material was kept in argon at 250°C to remove the glycol solvent. The as-prepared composite electrocatalysts were marked as pfSAC-Fe-X (X = 0.01, 0.05, 0.1, and 0.2, which represents the added mass ratios of the COFBTC).

DFT calculation

We used the Vienna ab initio simulation package herein for all DFT calculations. To describe the interactions of the core electrons within different atoms, projected augmented wave was used. For exchange-correlation interactions, we adapted the generalized gradient approximation by Perdew-Burke-Ernzerhof. The plane wave basis sets had a cutoff energy of 500 eV to describe all atoms’ valence electrons. The total energy of all the structures converged to 10−4 eV, and the Hellmann-Feynman force converged to 0.02 eV/Å. The Brillouin zone was set as the γ point. The graphene and COFBTC were optimized separately before building the composite structure model, which was built by periodical single-layer graphene-based systems with a 15-Å vacuum to avoid the interactions between the layers. The lattice constants of these models were 19.7 Å by 19.7 Å by 15 Å. Van der Waals correction was applied in all composite structure calculations. ISPIN values were considered when calculating the models with magnetism.

Assembly of the Zn-air battery

Zn-air battery was assembled with homemade frameworks. On the air electrode, the catalysts were coated along with a gas diffusion layer at an effective size of 1.0 cm by 1.0 cm. A polished zinc foil was assembled as the negative electrode with the same effective size. The electrolyte was composed of 8 M KOH and 0.5 M ZnO. The battery performance was tested using a constant current charging-discharging process by Land Instruments at room temperature.


Supplementary material for this article is available at

Supplementary Materials and Methods

Fig. S1. TEM images (associated with elemental mapping) of pfSAC, showing the existence and homogeneous distribution of C, Fe and N.

Fig. S2. TEM images (associated with elemental mapping) of pure graphene, in which no Fe and N was observed.

Fig. S3. XPS spectra and high-resolution XPS spectra of pfSAC-Fe-0.2.

Fig. S4. HAADF-STEM images of the pfSAC-Fe-0.2 sample.

Fig. S5. Fe K-edge EXAFS analysis of the Fe foil.

Fig. S6. Fe K-edge EXAFS analysis of pure COFBTC.

Fig. S7. Fe K-edge EXAFS analysis of pfSAC-Fe-0.2.

Fig. S8. Relaxed structures for top-vertical absorption sites.

Fig. S9. Relaxed structures for top-horizontal absorption sites.

Fig. S10. Relaxed structures for bridge-horizontal absorption sites.

Fig. S11. Relaxed structures for bridge-vertical absorption sites.

Fig. S12. DOS for pfSAC-Fe and COFBTC.

Fig. S13. The RRDE test using different pfSAC-Fe-0.2 loadings (from 40 to 304 μg cm−2) in 0.1 M KOH.

Fig. S14. Stability tests for pfSAC-Fe-0.2.

Fig. S15. The Fe K-edge XANES spectra of pfSAC-Fe-0.2 before and after ORR.

Fig. S16. TEM images (associated with elemental mapping) of the pfSAC-0.2 sample after 6000 cycles of CV test, showing the homogeneous distribution of C, Fe, and N.

Fig. S17. HAADF-STEM images of the pfSAC-Fe-0.2 sample after 6000 cycles of CV test.

Fig. S18. CVs, H2O2 yields, and electron transfer number of pfSAC-Fe-0.01, pfSAC-Fe-0.05, pfSAC-Fe-0.1, and pfSAC-Fe-0.2.

Fig. S19. LSV curves of graphene, pfSAC-Fe-0.2, and 20% Pt/C in O2-saturated 0.1 M KOH solution at a scan rate of 5 mV s−1 and a rotation speed of 1600 rpm.

Fig. S20. Kinetic current density of pfSAC-Fe-0.01, pfSAC-Fe-0.05, pfSAC-Fe-0.1, pfSAC-Fe-0.2, Pt/C, and graphene over 0.82–1.0 V versus RHE.

Fig. S21. LSV curves and the corresponding K-L plots of pfSAC-Fe-0.01, pfSAC-Fe-0.05, and pfSAC-Fe-0.1.

Fig. S22. Tafel plot and methanol crossover tests of the samples.

Fig. S23. Polarization and power density curves of assembled Zn-air batteries using Pt/C and pfSAC-Fe-0.2 as ORR catalysts with different loading masses.

Table S1. Fe K-edge EXAFS curve fitting parameters.

Table S2. Comparison of the ORR activity between pfSAC-Fe-0.2 and other nonprecious catalysts under basic conditions (0.1 M KOH) in literature.

Table S3. Inductively coupled plasma atomic emission spectroscopy results for Fe contents in COFBTC and the as-obtained pfSAC-Fe-0.2.

Table S4. Comparison of the performance of the zinc-air batteries based on pfSAC-Fe-0.2 and other electrocatalysts.

References (3853)

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


Acknowledgments: We thank L. Dai from Case Western Reserve University and J.-F. Chen and D. Cao at Beijing University of Chemical Technology for helpful discussion and suggestions; Z. Wen from Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences for the fitting and analysis of the XAFS results; and the 1W1B station for XAFS measurements in Beijing Synchrotron Radiation Facility (BSRF). Funding: This work was supported by the National Key Research and Development Program of China (2017YFA0206500), the Natural Science Foundation of China (21676020 and 21606232), the Beijing Natural Science Foundation (17L20060), the Fundamental Research Funds for the Central Universities (buctrc201420 and buctrc201714), the BUCT Fund for Disciplines Construction and Development (XK1502), the Young Elite Scientists Sponsorship Program by CAST (2017QNRC001), the Talent cultivation and open project (OIC-201801007) of State Key Laboratory of Organic-Inorganic Composites, the Distinguished Scientist Program at BUCT (buctylkxj02), and the “111” Project of China (B14004). Author contributions: P.P., L.S., and C.M. performed synthesis, structural characterizations, electrochemical tests, and Zn-Air battery tests. P.P. and Z.X. wrote this paper. F.H., X.W., and S.Z. performed DFT calculations. Z.X. supervised and led this project. All authors provided critical feedback, helped shape the research and manuscript, and commented on the manuscript. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.

Stay Connected to Science Advances

Navigate This Article