Research ArticleChemistry

A single iron site confined in a graphene matrix for the catalytic oxidation of benzene at room temperature

See allHide authors and affiliations

Science Advances  04 Dec 2015:
Vol. 1, no. 11, e1500462
DOI: 10.1126/sciadv.1500462


Coordinatively unsaturated (CUS) iron sites are highly active in catalytic oxidation reactions; however, maintaining the CUS structure of iron during heterogeneous catalytic reactions is a great challenge. Here, we report a strategy to stabilize single-atom CUS iron sites by embedding highly dispersed FeN4 centers in the graphene matrix. The atomic structure of FeN4 centers in graphene was revealed for the first time by combining high-resolution transmission electron microscopy/high-angle annular dark-field scanning transmission electron microscopy with low-temperature scanning tunneling microscopy. These confined single-atom iron sites exhibit high performance in the direct catalytic oxidation of benzene to phenol at room temperature, with a conversion of 23.4% and a yield of 18.7%, and can even proceed efficiently at 0°C with a phenol yield of 8.3% after 24 hours. Both experimental measurements and density functional theory calculations indicate that the formation of the Fe═O intermediate structure is a key step to promoting the conversion of benzene to phenol. These findings could pave the way toward highly efficient nonprecious catalysts for low-temperature oxidation reactions in heterogeneous catalysis and electrocatalysis.

  • single atom site
  • coordinatively unsaturated iron
  • benzene oxidation
  • Graphene
  • non-precious catalyst


Earth-abundant transition metal centers, such as coordinatively unsaturated (CUS) iron sites, can exhibit higher catalytic activity for reactions than precious metals. Yet, because of the instability of CUS sites, it is difficult to maintain the active structure of transition metal centers during a heterogeneous catalytic reaction. On the other hand, many successful examples can be found in enzymes such as cytochrome P-450 (1, 2), nitrogenase (3), and methane monooxygenase (4), as well as some homogeneous catalysts, where the organic ligands and proteins confine these CUS iron sites, making them highly active and stable (1, 57). In heterogeneous catalysis, however, preparation of the analogous CUS iron sites in supported catalysts with robust structures and high activity remains an attractive challenge (811). Our previous work demonstrated that the CUS ferrous sites, confined at the interface of precious metal Pt, are highly active and stable in activating oxygen at low temperatures (12, 13). However, the high cost of Pt prevents the commercialization of these catalysts. A major research thrust has been made to replace Pt with earth-abundant materials while maintaining the CUS ferrous structure. Graphene with a well-defined two-dimensional (2D) structure and high specific surface area shows high mechanical strength and thermal stability under realistic catalytic conditions (14, 15). Its unique structural and electronic properties render it a promising host to confine the CUS metal atoms in the matrix. Several recent works have demonstrated that the single metal atom can be successfully embedded in a graphene matrix through in situ electron beam irradiation in a transmission electron microscopy (TEM) system (1618). However, the pure metal atoms in graphene are mobile under irradiation (16, 17), implying their instability under realistic catalytic conditions. Moreover, it is difficult to obtain a sufficient quantity for catalytic applications using the irradiation method. FeN4 centers with CUS Fe sites in organic macrocycles have been proven to be stable structures, whereas the supported FeN4 macrocycles on substrates tend to aggregate during catalytic reactions because of the weak interactions between these macrocycles and substrates (19). Therefore, one possible route to stabilizing the CUS Fe sites in the graphene matrix is via the introduction of N atoms as an “anchor,” because the C–N bond has been proven to be highly stable in N-doped graphene (20, 21).

Here, we report one strategy to achieve a highly dispersed single FeN4 center with CUS Fe sites confined in a graphene matrix at a large quantity via high-energy ball milling of iron phthalocyanine (FePc) and graphene nanosheets (GNs) under controllable conditions. High-energy ball milling has been demonstrated as a powerful method to cut and reconstruct the chemical bonds of materials or molecules with necessary energy input (20, 2224). We prepared a series of graphene-embedded FeN4 (FeN4/GN) catalysts with different Fe content, that is, FeN4/GN-1.5 (1.5% Fe, see table S1), FeN4/GN-2.7, and FeN4/GN-4.0, by ball milling the composites of FePc and GN with appropriately chosen energies (see Materials and Methods for more details).


The typical morphology of FeN4/GN is presented in Fig. 1 (A to D) and figs. S1 and S2, and was obtained using low-voltage (80 kV) spherical aberration–corrected HRTEM. One can see homogeneously dispersed small black dots in the graphene matrix. Some are tagged by the red arrows and circles in Fig. 1 (A to D) and fig. S2, which could be assigned as single Fe atoms. The structure of single Fe centers was further evidenced by sub-angstrom resolution HAADF-STEM images (Fig. 1, G and H, and fig. S3), which show the atomic size and homogeneous distribution of the bright dots within the graphene matrix. Through EELS atomic spectra of the bright dots (Fig. 1, H and I), one can clearly see the presence of both Fe and N elements in one bright dot, suggesting the formation of Fe–Nx bonding. This indicates that Fe atoms observed in Fig. 1 (A to D, G, and H) and figs. S1 to S3 should be bonded with N atoms in the surroundings and further contacted with the graphene matrix, as shown in the atomic models (Fig. 1E) for the experimental structures (Fig. 1D), which is also highly consistent with the density functional theory (DFT)–simulated HRTEM image (Fig. 1F). Note that some disordered structures can also be observed around Fe atoms in some areas (fig. S4), implying the introduction of defects in the graphene network around some iron atoms during the high-energy ball milling. X-ray diffraction (XRD) (fig. S5) and Raman spectra (fig. S6) further indicate that there is no characteristic structural information of FePc, Fe, or FeOx observed in FeN4/GN samples, implying a well-dispersed feature of these Fe sites in FeN4/GN samples, which is highly consistent with the TEM and HAADF-STEM analysis. To obtain more atomic and electronic structure information of FeN4 centers in the graphene matrix, we performed low-temperature scanning tunneling microscopy (LT-STM, 4 K). Figure 1J shows a typical atomic-resolution STM image of a single FeN4 center embedded in the graphene matrix. The iron center is resolved as a bright spot, whereas neighboring atoms exhibit a higher apparent height than other carbon atoms in the graphene matrix. STM simulation (Fig. 1K) of an FeN4 center embedded in the graphene lattice is in agreement with the measured STM image (Fig. 1J), suggesting that the iron center significantly modifies the density of states of adjacent atoms. The bright dot in Fig. 1J is attributed to the iron center, whose neighboring C and N atoms are also electronically rich and appear brighter than carbon atoms located further away. Accordingly, STM contours of the bright spot, the corresponding conductance spectra (Fig. 1L), and the stability of the bright spot during scanning tunneling spectroscopy (STS) measurements all suggest that the FeN4 center is in the plane of graphene and forms stable bonds with neighboring carbon atoms. FeN4 in macrocycles often exhibits sharp electronic states near the Fermi level, corresponding to their highest occupied molecular orbital (HOMO)–lowest unoccupied molecular orbital (LUMO) levels (2527). On graphene, because of the large gap between the HOMO and LUMO levels of the FeN4 center, such states are often 1.5 to 2 eV away from the Fermi level (28). In Fig. 1L, STS measurements across the FeN4 center also show a sharp resonance state at −0.63 eV below the Fermi level, suggesting that the iron center strongly interacts with the graphene lattice and thus introduces a new electronic state near the Fermi level, not seen in isolated FePc molecules or N-doped graphene (29). To the best of our knowledge, this is the first time the well-defined FeN4 atomic structure in graphene has been observed by combining HRTEM/STEM with LT-STM/STS.

Fig. 1 Structural analysis of graphene-embedded FeN4 (FeN4/GN) catalysts.

(A to D) High-resolution transmission electron microscopy (HRTEM) images of FeN4/GN-2.7. The area with arrows and the dashed circles shows some typical single Fe atoms in the nanosheets. (E and F) Atomic models (E) and the corresponding simulated HRTEM images (F) for the structures in (D), where the FeN4/GN structures have been optimized. (G and H) High-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) images of FeN4/GN-2.7. (I) The electron energy loss spectroscopy (EELS) atomic spectra of Fe and N elements from the bright dots as shown by the red arrow in (H). The red circles show Fe and N signals, respectively. a.u., arbitrary units. (J) Low-temperature scanning tunneling microscopy (LS-STM) image of FeN4/GN-2.7, measured at a bias of 1.0 V and a current (I) of 0.3 nA (2 nm × 2 nm). (K) Simulated STM image for (J). The inserted schematic structures represent the structure of the graphene-embedded FeN4. The gray, blue, and light blue balls in (E), (J), and (K) represent C, N, and Fe atoms, respectively. (L) dI/dV spectra acquired along the white line in the inset image. U, 1.0 V; I, 0.3 nA; modulation frequency, 500 Hz; amplitude, 20 millivolts peak to peak; RC, 7 Hz.

X-ray absorption fine structure (XAFS) spectroscopy was used to further probe the chemical state and coordination structure of these confined Fe centers. As shown in Fig. 2A, the Fe K-edge of XANES in FeN4/GN samples exhibits a near-edge structure similar to that of the original FePc but is very different from those of Fe foil and Fe2O3, indicating that the valence state of Fe remains the same with FePc, which can be further confirmed by Fe 2p XPS analysis (fig. S7B). EXAFS of the Fe K-edge (Fig. 2B) shows that the magnitude of the FT spectra of the FeN4/GN samples also closely resembles the original FePc reference curve (30, 31). From the shape and amplitude of the first strong peak (with phase shift correction) in the FT plot, one can see that the bonding environment in the first shell of FeN4/GN samples is the same as that of FePc, suggesting that one Fe site connects four N atoms as the FeN4 structure in its precursor FePc. Furthermore, the N 1s XPS (Fig. 2E) reveals that the intensity of pyrrolic Nα (400.4 eV) (bonding with Fe) is almost unchanged whereas that of pyridinic Nβ (398.6 eV) (bonding with carbon on the outside macrocycle) is significantly reduced compared with FePc; this indicates that part of the pyridinic Nβ species has been destroyed during ball milling whereas pyrrolic Nα species are well retained in FeN4/GN samples. We further investigated the C K-edge XAS spectra of FeN4/GN samples to study the macrocyclic structure change during ball milling. As shown in Fig. 2C, the FeN4/GN samples show a strong π* and σ* band structure, indicating that the graphene matrix is still graphitized. It can be seen that the intensities of B and C features, considered as contributions predominantly from carbon atoms of the pyrrole rings (3234), have been obviously reduced in FeN4/GN samples compared with the FePc sample, indicating that some parts of the carbon atoms in the outside macrocyclic structure have also been destroyed. Meanwhile, the N K-edge XAS spectra of FeN4/GN samples (Fig. 2D) show that the intensity of the π* band (at ca. 398 eV) significantly decreases relative to that of the σ* band (at ca. 406 eV) for N, suggesting that the number of C═N bonds was significantly reduced and the FeN4 structure remains almost unchanged.

Fig. 2 Chemical state and coordination information of FeN4/GN catalysts.

(A and B) Fe K-edge x-ray absorption near-edge structure (XANES) (A) and Fourier transform (FT) extended x-ray absorption fine structure (EXAFS) (B) signals of FeN4/GN samples with various Fe content in comparison to FePc, Fe foil, and Fe2O3. (C and D) C K-edge (C) and N K-edge (D) x-ray absorption spectroscopy (XAS) spectra of FeN4/GN samples with various Fe content in comparison to that of FePc. (E) N 1s x-ray photoelectron spectroscopy (XPS) spectra of FeN4/GN samples with various Fe content in comparison to FePc. The inserted schematic structures represent the FePc molecule, where the pyrrolic N with Fe bonding is denoted as Nα and the pyridinic N with carbon bonding on the outside macrocycle is denoted as Nβ.

The above results demonstrated that the FeN4 centers have been successfully embedded into the matrix of GNs via high-energy ball milling of FePc and GN. In one proposed mechanism, described in fig. S8, the outside macrocyclic structure of FePc can be destroyed during the ball milling, the residual isolated FeN4 centers will interact with the graphene at the defected site, and the adjacent carbon atoms of FeN4 can further reconstruct with the high energy of ball milling, finally leading to the formation of the FeN4 centers embedded into the graphene matrix. Our previous work indicated that isolated Fe atoms embedded within a silicide matrix showed high activity and long-term stability toward direct conversion of methane to ethylene and hydrocarbons (35). Thus, this graphene-confined single CUS iron site is expected to have high performance for catalytic reactions.

The direct catalytic conversion of benzene to phenol is one of the most active topics in fundamental and applied research (3639). Different catalysts including Ti-containing zeolites, palladium membranes, and transition metal (such as Fe, Cu, and V)–based oxides or chelates have been widely investigated for the direct conversion of benzene to phenol. This reaction is usually carried out at 50° to 140°C, because it is very difficult to directly proceed at room temperature owing to the highly stable C–H bond of benzene (36, 4042). Here, we found that the FeN4/GN samples showed a high activity and selectivity for phenol at room temperature. The oxidation of benzene was conducted at 25°C with hydrogen peroxide as the oxidant. With the increase of the Fe content in graphene, the activity and the yield of phenol first increased quickly and then decreased (Fig. 3A and table S2). This trend in performance with the Fe content in graphene can be attributed to the observation that a moderate amount of FeN4 can promote both the dispersion of FeN4 centers in graphene and their bonding with graphene, whereas a higher content of FeN4 will lead to the agglomeration of FeN4. The optimized FeN4/GN-2.7 catalyst has a turnover frequency of 84.7 hour−1 for benzene conversion within the initial 5 min (Fig. 3B) and can achieve a benzene conversion of 23.4% and phenol yield of 18.7% in 24 hours (Fig. 3B and table S3). Some residual FePc dissolved in reaction solution may have contributed to converting benzene to phenol, but the contribution should be minor because the FeN4/GN-2.7 sample shows significantly better activity compared with the FePc monomer, despite the latter having more Fe sites. For comparison, in a blank experiment without catalyst, no obvious activity was observed over a 24-hour run (table S2). When the graphite flake (GF) and GN were used as the catalyst, only low conversions of benzene were observed (Fig. 3A and table S2), that is, 0.6% for GF and 5.4% for GN under the same conditions. Considering that the edges and defects of graphene may contribute to the activity as reported previously by Deng et al. (20), it is reasonable that the GN exhibited a higher conversion of benzene than GF because the GNs have more edges and defects. In addition, the low-temperature O2 temperature-programmed desorption (TPD) measurement shows that the FeN4/GN-2.7 has a significantly higher adsorption capacity of O2 compared with GN and GF (fig. S9). O2 can be easily adsorbed on the FeN4 structure in metal porphyrin or phthalocyanine according to previous studies (43, 44), which indicates that the FeN4/GN-2.7 had more active sites. Furthermore, we found that the FeN4/GN-2.7 catalyst can even proceed efficiently at 0°C with a phenol yield of 8.3% under 24 hours (table S4) and can remain stable after six cycles (fig. S10), further supporting its excellent catalytic performance.

Fig. 3 The performance and reaction process of the catalytic oxidation of benzene to phenol over FeN4/GN catalysts.

(A) The performance of the direct oxidation of benzene to phenol by FeN4/GN samples compared with GF, GN, and FePc. Reaction conditions: 50 mg of catalyst, 0.4 ml of benzene, 6 ml of H2O2 (30%), and 3 ml of CH3CN in a pressure vessel at 25°C for 24 hours. (B) The phenol yield of FeN4/GN-2.7 for the direct oxidation of benzene to phenol with different reaction times. (C and D) Fe K-edge XANES (C) and FT EXAFS (D) signals of FeN4/GN samples with H2O2 treatment in comparison to their corresponding original samples.

To gain further insights into the activity of FeN4/GN toward benzene oxidation, we carried out DFT calculations to investigate the reaction mechanism of the process (see Materials and Methods and figs. S11 and S12 for more details on calculation). A model of the FeN4 structure embedded in graphene was adopted according to experimental characterization. Figure 4A shows that the formation energy of the FeN4 center in graphene matrix (FeN4/GN) is significantly lower than that of the single-atom Fe in a pure graphene matrix (Fe/GN), suggesting that the N atoms can be used as an anchor to enhance the stability of Fe atoms in graphene, which supports the experimental results. Furthermore, the free energy profile and reaction pathway of benzene oxidation on the confined iron site are depicted in Fig. 4 (B and C). A H2O2 molecule can be easily dissociated on the confined iron site by forming an Fe═O intermediate and releasing one H2O molecule, followed by the dissociation of another H2O2 on the other side of the iron atom with an energy barrier of 0.55 eV by forming an O═Fe═O center. The O species of the O═Fe═O is active for the adsorption of the benzene molecule via the formation of a C–O bond with an energy barrier of 0.59 eV. In comparison, the direct adsorption of benzene on the O site of the Fe═O site is energetically unfavorable and needs an additional free energy of 0.86 eV (fig. S12). The benzene adsorbed on the O═Fe═O site can transform to phenol via the transfer of one adjacent H atom from C to O with a barrier of 0.35 eV. The Fe═O site can be regenerated in a reaction cycle after desorption of the phenol from the iron. The highest energy barrier in the reaction pathway occurs at the adsorption of benzene on the O═Fe═O site, which is only 0.59 eV and moderate for low-temperature reactions. The formation of Fe═O/O═Fe═O intermediates on the FeN4 center is also evidenced by XAFS analysis of the FeN4/GN samples after the H2O2 treatment. As shown in Fig. 3C and fig. S13, after the H2O2 treatment, the XANES of the Fe K-edge shows almost no energy shift, whereas the pre-edge peak, that is, the Fe 1s-to-3d transition in all FeN4/GN samples, would increase and broaden, probably because the formation of Fe═O leads to Fe 3d mixing with O 2p and thus destroys the D4h symmetry of FeN4 according to the literature (30, 45). The Fe═O interaction (relative to the Fe–N bonding) will increase the unoccupied state of Fe because of the electronegativity of O; thus, the pre-edge is more intense. EXAFS of Fe K-edge further confirmed the hypothesis. One can see that the amplitude of the first strong peak in the FT plot of these samples was significantly enhanced after the H2O2 treatment (Fig. 3D), suggesting that the coordination number of the Fe center sharply increases, which likely originates from the formation of Fe═O/O═Fe═O bonds during the reaction. EXAFS fitting of these results shows that these original FeN4/GN samples have an average coordination number of about 4 (fig. S14 and table S5), which is almost the same as that of FePc but less than the maximal coordination number of 6. Thus, the iron sites in these samples are also CUS. With further H2O2 treatment of the FeN4/GN sample, the coordination number of the iron site increases (table S5). These results supported the DFT calculation that the CUS iron site can effectively activate H2O2 and form Fe═O bonds. 57Fe Mössbauer spectra (fig. S15 and table S6) also indicate that the symmetrical O═Fe═O structure will significantly increase in FeN4/GN when treated with H2O2, whereas the O═Fe═O structure will decrease again when further treated with benzene during the catalytic reaction. The above experimental results indicate that the FeN4 centers play an important role in the adsorption and activation of oxygen, which are in agreement with reaction cycles resulting from DFT calculations.

Fig. 4 Theoretical analysis of the FeN4/GN structure and the catalytic reaction process by DFT calculations.

(A) The formation energies of FeN4/GN and Fe/GN structures. The formation energy is calculated as follows: EFe-embeddedEFe-bulkE(N)GN, where EFe-embedded and EFe-bulk are the total energies of FeN4/GN and the Fe/GN structure and an Fe atom in Fe bulk, respectively, and E(N)GN is the total energy of the optimized structure of FeN4/GN or Fe/GN with the Fe atom removed from the system. (B) Free energy diagram of the oxidation of benzene to phenol on FeN4/NG. The gray, blue, light blue, red, and white balls represent C, N, Fe, O, and H atoms, respectively. (C) Scheme for the reaction mechanism of the oxidation of benzene to phenol on FeN4/NG.

The reaction pathway of benzene oxidation on the iron site of FePc was also calculated for comparison. Figure S12 shows that the reaction mechanism and free energy diagram of benzene oxidation on FePc is similar to that on FeN4/GN, but the dissociation energy of the first H2O2 forming Fe═O on FeN4/GN is lower than that on the FePc monomer. Bader charge analysis shows that the iron of FeN4/GN has an additional 0.14 electrons compared with that of FePc, leading to the O absorbed on FeN4/GN obtaining an additional 0.10 electrons compared with that on FePc, which then energetically favors the formation of the Fe═O bond on FeN4/GN. Therefore, the enhanced benzene oxidation activity of FeN4/GN can be attributed to both the intrinsic activity improvement of the active sites and the high dispersion of these CUS sites in FeN4/GN compared with bulk FePc.


In summary, the CUS single-iron site has been confined in the GN matrix through a one-step ball milling synthesis. The unique 2D structure of this catalyst provides a well-defined model for understanding the nature of the catalytic oxidation reaction on FeN4/GN catalysts by means of experiments and DFT calculations. In this system, the FeN4 center is highly dispersed and well stabilized by the graphene matrix, which subsequently enhances the activity and stability for the oxidation of benzene to phenol. This reaction can proceed efficiently at room temperature and even at temperatures as low as 0°C. DFT calculations indicate that the catalytic activity arises from the confined iron sites and the activation barriers are quite moderate for reactions to proceed at room temperature, in agreement with experimental results. These findings pave the way toward the design of highly efficient nonprecious catalysts for catalytic oxidation reactions at low temperatures.


Raw materials

GFs (99.8%, metals basis) were purchased from Alfa Aesar. FePc (96%) was purchased from Acros Organics.

Synthesis of GNs

GNs were prepared following the same procedure used in our previous report (20). In a typical experiment, 2.0 g of GF and 60 g of steel balls (1 to 1.3 cm in diameter) were put into a hardened steel vial inside a glove box and purged with high-purity Ar (99.999%) for 20 min before the vials were sealed. Ball milling was carried out at 450 rpm for 20 hours.

Synthesis of FeN4/GN

A combined mass of 2.0 g of FePc and GN composites with a desired ratio and 60 g of steel balls (1 to 1.3 cm in diameter) was ball-milled following the same procedure as the GN synthesis. A series of FeN4/GN samples with different Fe content were prepared, that is, FeN4/GN-1.5 [1.5% Fe, see inductively coupled plasma (ICP) data in table S1, the same below], FeN4/GN-2.7, and FeN4/GN-4.0 from the precursor FePc and GN with a ratio of 15, 30, and 45%, respectively. The utility ratio of FePc in the final FeN4/GN catalysts is around 90% according to the Fe content analysis using inductively coupled plasma atomic emission spectroscopy (ICP-AES).


HRTEM was carried out using an image spherical aberration–corrected TEM system (FEI Titan 80-300). An acceleration voltage of 80 kV was chosen to achieve enough resolution while maintaining the structure of the graphene. STEM and EELS were performed on a JEOL ARM200F equipped with double aberration correctors and a cold field emission gun operated at 80 kV. STEM images were recorded using a HAADF detector with a convergence angle of 30 mrad and a collection angle between 90 and 370 mrad. Under these conditions, the spatial resolution is ca. 0.08 nm. STM and STS were acquired using a commercial Createc LT-STM system with base pressures below 7.0 × 10−11 mbar. The sample was dispersed in petroleum ether and further dripped on the surface of HOPG (highly oriented pyrolitic graphite). The sample was then transferred to the Createc LT-STM system. Before imaging, the sample was degassed at ~450 K to remove impurities absorbed on the surface. STM experiments were performed at liquid He temperatures at a constant current mode using an electrochemically etched W tip. ICP-AES was conducted in Shimadzu ICPS-8100. The samples for ICP-AES analysis were first heated at 600°C for 12 hours in air, then treated with hydrochloric acid in Teflon-lined autoclaves at 120°C for 12 hours, and finally transferred to volumetric flasks. During XPS measurements, Mg Kα radiation (1253.6 eV) with a power of 200 W and a pass energy of 50.0 eV was used. XAS measurements were conducted at the SGM (11ID-1) beamline of the Canadian Light Source (CLS). Fe K-edge XAFS spectra of the catalysts were recorded at the SXRMB (06B1-1) beamline of the CLS and the BL14W1 beamline of the Shanghai Synchrotron Radiation Facility (SSRF). 57Fe Mössbauer spectroscopy analysis was conducted on a Topologic 500A spectrometer with a proportional counter. 57Co(Rh) was used as the radioactive source, and the Doppler velocity of the spectrometer was calibrated with α-Fe foil. The spectra were fitted with appropriate superpositions of Lorentzian lines using the MossWinn 3.0i program. XRD was performed on a Rigaku D/MAX 2500 diffractometer with Cu Kα radiation (λ = 1.5418 Å) at 40 kV and 200 mA. Raman spectroscopy was performed on a Jobin Yvon LabRAM HR 800 instrument with a 532-nm excitation laser at a power of 0.7 mW. O2 TPD measurements were carried out using the AutoChem II 2920 with a flowing 5% O2/He stream (50 ml min−1) at −50°C. The samples were pretreated with He at 250°C for 1 hour to remove the adsorbed gaseous impurities before the TPD test.

Catalytic benzene oxidation evaluation

Benzene oxidation reaction was carried out in a 50-ml Teflon-lined stainless steel reactor with 0.4 ml of benzene, 6 ml of H2O2 (30%), and 3 ml of CH3CN at 25° or 0°C. After the reaction, an additional 20 ml of CH3CN was added to transfer the products and 0.2 ml of toluene was also added as an internal standard. The products were analyzed with Agilent 1260 Infinity HPLC using a Unitary C-18 column. Before the analysis, the products were filtered by a syringe with a filter head.

DFT calculations

DFT calculations were performed using the Vienna Ab-initio Simulation Package (4648). The projector augmented-wave pseudopotentials and a cutoff energy of 400 eV for the plane-wave basis set were adopted (49, 50). The generalized gradient approximation method with Perdew-Burke-Ernzerhof functionals for the exchange-correlation term was used (51, 52). The Monkhorst-Pack scheme was used for sampling the Brillouin zone (53). The FeN4/GN model was set in a 6 × 6 supercell of graphene (fig. S11A). The vacuum thickness between graphene layers was set as 15 Å to avoid interlayer interactions. The FePc monomer model was set in a 25 × 25 × 16 Å rectangular box (fig. S11B). Spin polarization was considered throughout the calculations. The transition states were searched using the constrained minimization approach (5456). The free energies (G) of the reactants, surface intermediates, and products were obtained using the equation G = Etotal + ZPE − TS, where Etotal is the total energy of the species, ZPE is the zero point energy, and S is the entropy.

EXAFS analysis

FT EXAFS spectra of the FeN4/GN samples, as well the Fe foil and FePc reference materials, were generated and fitted using WinXAS (57). Scattering paths used in the fitting process were calculated ab initio using FEFF 8.2 (58) and models of the FePc and FeN4/GN structures (fig. S11). The Fe foil spectrum was used to determine an empirical S02 value (0.86), which was then fixed for the fitting of all subsequent samples.


Supplementary material for this article is available at

Fig. S1. HRTEM images of FeN4/GN-2.7.

Fig. S2. HRTEM image of FeN4/GN-2.7 with the red circles showing some typical single Fe atom positions in the graphene network.

Fig. S3. HAADF-STEM image of FeN4/GN-2.7.

Fig. S4. HRTEM image of FeN4/GN-2.7 with the red circles showing some Fe atoms with different defects in the surroundings.

Fig. S5. XRD patterns of graphite, GN, FeN4/GN-1.5, FeN4/GN-2.7, FeN4/GN-4.0, and FePc.

Fig. S6. Raman spectra of FeN4/GN samples in comparison to their parent materials FePc, GN, and graphite.

Fig. S7. XPS spectra of FePc, FeN4/GN-4.0, FeN4/GN-2.7, and FeN4/GN-1.5.

Fig. S8. Scheme of a proposed mechanism for synthesis of FeN4/GN via a facile ball milling method.

Fig. S9. Low-temperature O2 TPD profiles of FeN4/GN-2.7, GN, and GF.

Fig. S10. The recycling experiments of FeN4/GN-2.7.

Fig. S11. Models of FeN4/GN and the FePc monomer in the DFT calculations.

Fig. S12. Free energy profile of the benzene oxidation reaction intermediates on the iron site of the FePc monomer and FeN4/GN.

Fig. S13. Fe K-edge XANES signal of FeN4/GN samples with H2O2 treatment in comparison to their corresponding original samples.

Fig. S14. The Fe K-edge EXAFS analysis of FeN4/GN samples before and after H2O2 treatment.

Fig. S15. Room-temperature 57Fe Mössbauer spectra of FeN4/GN-2.7, FeN4/GN-2.7-H2O2, and FeN4/GN-2.7-H2O2-Ben.

Table S1. The elemental compositions of FePc, FeN4/GN-4.0, FeN4/GN-2.7, and FeN4/GN-1.5 estimated from XPS and ICP measurements.

Table S2. Catalytic performance of different samples for the direct oxidation of benzene to phenol.

Table S3. Catalytic performance of FeN4/GN-2.7 for the direct oxidation of benzene to phenol with different reaction times.

Table S4. Catalytic performance of different samples for the direct oxidation of benzene to phenol at 0°C.

Table S5. Fitting parameters for the analysis of the EXAFS spectra of FeN4/GN samples with H2O2 treatment in comparison to their corresponding original samples.

Table S6. Fitting parameters for the 57Fe Mössbauer spectra in fig. S15.

References (5961)

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 T. Regier at the CLS for his assistance on XAS measurements and BL14W1 beamline of the SSRF for assistance on XAFS measurements. Funding: This work was supported by the National Natural Science Foundation of China (grant nos. 21321002, 21303191, and 51420105003) and the Strategic Priority Research Program of the Chinese Academy of Sciences (grant no. XDA09030100). Author contributions: X.B. and D.D. supervised the work and designed the experiments. D.D. and X. Chen prepared the samples and performed most of the experiments. L.Y. performed the DFT calculations for the catalytic reactions. X.W., T.X., and L.S. performed HRTEM microscopy and simulation. H.Y., H.T., and J.L. performed HAADF-STEM microscopy. Q.L., Y.L., and F.Y. performed STM microscopy. H.L. and J.X. simulated the STM images. Y.H., R.S., P.D., and J.Z. performed XAFS characterization. J.W. performed Mössbauer spectroscopy analysis. P.N.D. and P.Z. performed EXAFS analysis and fitting. X. Cui, J.D., and X.P. helped with the sample preparation and evaluation of catalytic reactions. D.D., X. Chen, and X.B. interpreted the data and wrote the paper. All authors discussed the results 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 at dhdeng{at}

Stay Connected to Science Advances

Navigate This Article