Identification of catalytic sites in cobalt-nitrogen-carbon materials for the oxygen reduction reaction.

Single-atom catalysts with full utilization of metal centers can bridge the gap between molecular and solid-state catalysis. Metal-nitrogen-carbon materials prepared via pyrolysis are promising single-atom catalysts but often also comprise metallic particles. Here, we pyrolytically synthesize a Co-N-C material only comprising atomically dispersed cobalt ions and identify with X-ray absorption spectroscopy, magnetic susceptibility measurements and density functional theory the structure and electronic state of three porphyrinic moieties, CoN4C12, CoN3C10,porp and CoN2C5. The O2 electro-reduction and operando X-ray absorption response are measured in acidic medium on Co-N-C and compared to those of a Fe-N-C catalyst prepared similarly. We show that cobalt moieties are unmodified from 0.0 to 1.0 V versus a reversible hydrogen electrode, while Fe-based moieties experience structural and electronic-state changes. On the basis of density functional theory analysis and established relationships between redox potential and O2-adsorption strength, we conclude that cobalt-based moieties bind O2 too weakly for efficient O2 reduction.Nitrogen-doped carbon materials with atomically dispersed iron or cobalt are promising for catalytic use. Here, the authors show that cobalt moieties have a higher redox potential, bind oxygen more weakly and are less active toward oxygen reduction than their iron counterpart, despite similar coordination.

Introduction

The transition from fossil to renewable energies is necessary to meet the rising energy demand while minimizing anthropogenic climate change and urban pollution[1-3]. Electrochemical energy conversion will play an increasing role for the storage of renewable electricity, production of fuels and their conversion into electricity. Hydrogen is an interesting energy vector since it can be produced via water electrolysis (hydrogen evolution reaction, HER) and later oxidized in H2/air fuel cells to reform water and electricity on demand (hydrogen oxidation reaction, HOR). However, the electrochemical oxygen reactions (oxygen reduction reaction (ORR) and oxygen evolution reaction (OER)) are slow, limiting the roundtrip efficiency[4]. While acidic electrolytes are more restrictive than alkaline ones regarding the breadth of catalysts that may be stable under ORR/OER conditions, the advent of highly conductive and stable proton exchange membranes (PEM) has hitherto favoured the development of acidic fuel cells and electrolyzers. Their major drawback is, however, the need for platinum-group metals to catalyze the ORR and OER.

Following the pioneering report of Jasinski on the ORR activity of cobalt phthalocyanine[5], advanced Metal–N–C materials have since 1989 been prepared by pyrolyzing separate metal, nitrogen and carbon precursors[6]. Although the activity and durability of such catalysts have been improved[7-11], the identification of the active-site structure has lagged behind due to the non-crystallographic order of metal atoms in the most active sites and the simultaneous presence of crystalline metal phases. Metal ions coordinated with pyridinic nitrogen atoms embedded in a graphene matrix with hexagonal atomic arrangement, such as the MeN4C10 moiety, have for a long time been viewed as the most probable site structure[12-18]. Other studies hypothesized moieties integrated in disordered carbon sheets involving non-hexagonal rings, such as the MeN4C12 moiety[19, 20]. Calle-Vallejo et al.[21] investigated with density functional theory (DFT) the adsorption energy of oxygen intermediates on MeN4C12 and MeN4C10 moieties, revealing differences up to 0.7 eV. This highlights the importance of precisely determining the local site structures in Me–N–C materials for deciphering their reactivity.

Using X-ray absorption near-edge structure (XANES) spectroscopy we recently identified the active-site structure in pyrolyzed Fe–N–C catalysts as being a porphyrin-like FeN4C12 moiety, in contrast with FeN4C10 or FeN2+2C4+4 moieties previously assumed[22]. The formation of FeN4C12 moieties requires a strongly-disordered host material, in line with experimental observations[9, 23]. Although moieties present in Fe–N–C materials are now better described[22, 24–26], active sites in Co–N–C materials are still poorly identified. Co–N–C catalysts are ORR-active and advantageous vs. Fe–N–C since they produce less radical oxygen species, leading to improved fuel cell durability[27]. They also have recently been reported to catalyze the HER[28, 29]. The existence of atomically dispersed cobalt atoms on N-doped graphene layers has been unambiguously demonstrated by mass spectroscopy and scanning transmission electron microscopy (STEM)[29, 30]. DFT investigations also concluded the energetically favoured formation of CoNC moieties[20, 21, 31–34]. Extended X-ray absorption fine structure (EXAFS), previously used to investigate CoNC moieties[35-37], is insufficiently sensitive to the spatial arrangement of coordinating light elements surrounding the absorbing cobalt nucleus.

The present study identifies with XANES the detailed structure of CoNC moieties in a Co–N–C catalyst free of cobalt particles. Similar to Fe–N–C, the porphyrin-like CoN4C12 moiety is found to be the likeliest structure. In contrast to Fe–N–C, defective cobalt porphyrinic moieties are also viable candidates. Also metal-discriminating, operando XANES spectroscopy shows no change with electrochemical potential in the entire ORR range for Co–N–C, but a major change for Fe–N–C. The spectral change for Fe–N–C is correlated with a change of oxidation state in that potential region, while the lack of spectral change for cobalt is correlated with an unmodified oxidation state in that potential region.

Results

Absence of Co–Co bonds in Co0.5 evidenced with EXAFS and TEM

We synthesized a Co–N–C catalyst with the same approach developed for preparing Fe–N–C catalysts comprising only FeNC moieties, and labeled it Co0.5, with 0.5% Co in the catalyst precursor before pyrolysis (see Methods).[22] The EXAFS spectrum of Co0.5 was analyzed assuming a variable number of nitrogen ligands in the equatorial plane and additional oxygen atoms in the axial position. The best-fit analysis is shown in Fig. 1, revealing 4- and/or 5-fold coordination structures. No Co–Co backscattering signal was needed to obtain an excellent fit. The dominant contribution to the EXAFS signal is given by the Co–N bonds of the first coordination shell (top of Fig. 1). The lower panels show the corresponding fits in the Fourier-transformed space. The peak at 1.45 Å is associated with the Co–N and/or Co–O first-shell contributions while the peak at 2.44 Å is due to Co–C. The structural parameters obtained from the fitting are reported in Supplementary Table 1. The Co–N distance, 1.95–1.96 Å, agrees with values determined by DFT for cobaltous-porphyrins.[38] However, the EXAFS experimental data can be reproduced with the same accuracy by a CoN4 or an O-CoN4 structure, as reported for Fe0.5 (ref. [22]). Moreover, the uncertainty on the nitrogen coordination number (Supplementary Table 1) leaves room for hypothesizing moieties with a lower coordination, as investigated later in this study. The absence of metallic Co–Co bonds in Co0.5 was independently confirmed by TEM (Supplementary Fig. 1). The selected area electron diffraction (SAED) pattern shows only broad rings, as expected for amorphous carbon. No sign of reflections originating from metallic particles were observed, and no evidence of particles was found in the images. Multiple areas of the sample were analyzed, and the images and SAED pattern in Supplementary Fig. 1 are representative.

Fig. 1

Co K-edge EXAFS analysis. Spectra of Co0.5 with a CoN4 moiety having zero (left) or one (right) oxygen atom in the axial direction. Cobalt, nitrogen and oxygen atoms are represented in purple, blue and red, respectively (carbon atoms in the second coordination sphere are not represented). a, b Curves from top to bottom: Co–N, Co–O and Co–C γ(2) two-body signals and the N-Co-N γ(3) three-body signal included in the fit, the total signal (red line) superimposed on the experimental one (black dots). c, d the fit in the Fourier-transformed space. No phase-shift correction was applied to the Fourier transforms

Electrocatalytic properties of Co0.5 towards the ORR

The ORR activity and selectivity of Co0.5 was then investigated with a rotating ring disk electrode (RRDE) in acid medium (Fig. 2a). The ORR current density at 0.8 V vs. a reversible hydrogen electrode (RHE) is 0.9 mA cm−2 for a loading of 800 µg cm−2 (12 µgCo cm−2), which compares well with state of art Co–N–C catalysts in acid medium, in spite of low cobalt content in Co0.5 (refs [39-41]). The Pt/C catalyst reaches 1.0 mA cm−2 at 0.91 V vs. RHE, and its higher activity is due to its higher atom-specific ORR activity (the higher Pt loading, 40 µgPt cm−2, resulting in identical number of metal atoms per cm2 as for cobalt, due to the 3.3 times higher molar mass of Pt vs. Co).

Fig. 2

ORR electrochemical characterization. Measurements of Co0.5, Fe0.5, N-C and a Pt/C benchmark catalyst. Upper panels are RRDE measurements at various rotation rates and loadings performed in O2-saturated 0.1 M H2SO4 aqueous solution. a ORR current density at the disk. The total catalyst loading for Co0.5, N-C and Pt/C (5 wt% Pt/C) was 0.8 mg cm−2. b % H2O2 measured during ORR at the ring. The total catalyst loading for Co0.5 was 0.8 or 0.2 mg cm−2. Lower panels are iR-corrected PEMFC polarization curves with Co0.5, Fe0.5 or Pt/C at the cathode, presented in linear scale c and semi-logarithmic scale d. The fuel cell temperature was 80 °C, pure H2 and O2 gases were fully humidified at the cell temperature, the gauge pressure was 1 bar. The cathode catalyst loading for Co0.5 and Fe0.5 was 4 mg cm−2 (60 µgmetal per cm2). For Pt/C at the cathode, 5% Pt/C was diluted to 1% Pt with uncatalysed Vulcan carbon black, and 4 mg cm−2 of 1% Pt/C (40 µgPt per cm2) deposited

Regarding selectivity, the % peroxide released by Co0.5 during ORR at 0.5–0.8 V vs. RHE is 3–5 % in 0.1 M H2SO4 at a loading of 800 µg cm−2 (Fig. 2b). Co0.5 is thus less selective to water formation than Fe0.5, for which we measured <0.5 % H2O2 in acid, at equal catalyst loading.[22] The peroxide production was also measured at a lower loading, since it has been reported that, for some Fe–N–C catalysts, this can result in greatly increased detection of H2O2 (ref. [41]). The % H2O2 released by Co0.5, however, only increased up to 11% at 200 µg cm−2 (Fig. 2b). By analogy with what is known for unpyrolyzed Metal-N4 macrocycles, this result may be explained by the generally higher formal potential of MetalIII/MetalII for cobalt vs. iron macrocycles, leading to weaker interaction with oxygen intermediates, including H2O2 (ref. [42]). Square-wave-voltammetry in acidic medium supports this possible explanation, identifying redox-peak positions at 1.25 V and 0.75 V vs. RHE for Co0.5 and Fe0.5, respectively (Supplementary Fig. 2). Although Co–N–C catalysts typically release more peroxide during the ORR, subsequent reactions between H2O2 and MeNC moieties produce less radical oxygen species on cobalt- than iron-centers, which is beneficial for fuel cell durability[27]. Also, the insignificant ORR activity in acidic medium of N–C (Fig. 2a) means that the cobalt moieties are the most active-site for ORR in Co0.5. This catalyst was then investigated in PEM fuel cell and compared to Fe0.5 and Pt/C (Fig. 2c, d). For comparison, we prepared a Pt-based cathode with similar metal-site density (1% Pt/C) and similar thickness as Co0.5 or Fe0.5 cathodes. For carbon-rich catalysts, the thickness is set by the catalyst loading, here 4 mg cm−2, typically resulting in 80–100-µm-thick layers[10]. Although the ORR activity at 0.8 V of Co0.5 is lower than that of Fe0.5 (Fig. 2d), the lower apparent Tafel slope (TS) observed with Co0.5 at high potential leads to a higher power performance at E < 0.55 V. The apparent TS of 83 mV/decade for Fe0.5 probably results from the outcome of an intrinsic TS with value > 83 mV per decade combined with the potential-dependence of the Fe(III) and Fe(II) coverages[43], with Fe(III) blocking the ORR due to overly strong O2 adsorption. A similar phenomenon occurs for Pt, switching from surface-oxidized state to a surface free of oxygen adsorbates at ~ 0.8 V vs. RHE[44, 45]. In contrast, Co0.5 does not change oxidation state throughout the ORR region (Supplementary Fig. 2) and the TS-value of 65 mV per decade observed at high potential is, therefore, intrinsic to the ORR mechanism on CoNC moieties. Overall, and for similar cathode thickness, the gap between the curves for 1% Pt/C and Co0.5 or Fe0.5 is only a factor two at 0.6 V. A key issue is that, while Pt can be deposited as nanoparticles up to 50 wt% on carbon, the ORR activity of Metal–N–C materials levels off at 2–3 wt% metal[46]. Further work is thus needed to (i) improve reactant transport in thick Me–N–C electrodes, (ii) increase the number of active sites or (iii) their turnover frequency. The latter two tasks require identifying the structure of active moieties, object of the present study.

Structural identification of Co-moieties in Co0.5 by XANES

To overcome the limitations of EXAFS, we resorted to XANES, which is more sensitive to the geometrical arrangement of atoms around the photo-absorber. Our XANES analysis was first validated with cobalt phthalocyanine (Co(II)Pc). The agreement between experimental and calculated spectra of Co(II)Pc is excellent in the whole energy range (Supplementary Fig. 3 ). For this fit, the square residual function, Rsq (Methods section) is 1.06. The optimized structure resulted in a Co–N distance of 1.90(2) Å, which agrees with an X-ray diffraction determination[47]. We then calculated the Co K-edge XANES spectra of various CoNC candidate sites for Co0.5, including structures previously investigated with DFT[20, 31–34]. Some of these clusters are depicted in Supplementary Fig. 4 together with the results of the fitting procedure. Supplementary Fig. 4a, c report the XANES fitting performed assuming a CoN4C10 moiety enclosed in a graphene plane and a CoN2+2C4+4 moiety bridging two graphene planes, respectively. For the calculations, only the in-plane distances between cobalt and nitrogen atoms were allowed to vary. For these hypothetical sites, the fit between experimental and calculated XANES spectra is unsatisfactory (Rsq = 2.97 and 2.81, respectively) and the fit quality was little improved by adding an axial oxygen (Supplementary Fig. 4b, d). The disordered arrangement of nitrogen and carbon atoms surrounding cobalt ions has been recently revealed by STEM for another Co–N–C catalyst[29]. This suggests that the arrangement of carbon atoms around cobalt ions is different from that in graphene.

We then considered two defective sites derived from CoN4C10, namely pyridinic CoN3C10,pyr and CoN2C8 (Supplementary Fig. 5a, c) and two defective sites derived from CoN2+2C4+4, namely CoN2+1C4+3 and CoN2C4 (Supplementary Fig. 5e, g). None of these four pyridinic-defective candidate sites correctly reproduced the experimental XANES spectrum of Co0.5. Again, the addition of an O2 molecule in end-on mode did not improve much the fit quality (Supplementary Fig. 5b, d, f, h). The structural parameters obtained from these analyses are summarized in Supplementary Table 2.

We then considered a cobalt moiety based on a porphyrinic architecture, analogous to that found for the Fe-based active-site geometry in Fe0.5 (ref. [22]). Fig. 3a, b show the XANES analyses performed on the square-planar CoN4C12 moiety, with and without an axial O2 adsorbed end-on. These fits show an excellent agreement between experimental and theoretical spectra and correspond to a moiety where the cobalt ion is coordinated in-plane by four nitrogen atoms at 1.96–1.97 Å, and possibly by one oxygen molecule at 2.23 Å (Table 1, rows 1 and 2). Two defective sites, derived from the porphyrinic CoN4C12 moiety by subtracting one or two nitrogen atoms (porphyrinic CoN3C10,porp and CoN2C5 model sites), also correctly reproduce the experimental spectrum (Fig. 3c, e, respectively). For the CoN3C10,porp motif, the XANES analysis resulted in a Co–N bond length of 1.96 Å when O2-free, and 1.99 Å when O2 is adsorbed end-on without out-of-plane displacement of cobalt (Table 1, rows 3 and 4). The XANES analysis also revealed the CoN2C5 edge-defect, representing a halved porphyrinic moiety, binding O2 end-on with a Co–O bond length of 1.90 Å with O2 in the same plane as the CoN2C5 motif (Fig. 3e, f and Table 1, rows 5 and 6). For the CoN2C5 motif only, the bound O2 molecule was necessary to reach a good fit (Rsq-value of 1.12–1.16, Table 1). For the CoN4C12 and CoN3C10,porp moieties (with or without adsorbed O2), the Co–N bond lengths determined by XANES and EXAFS are in good agreement, and so is the case for the CoN2C5 moiety free of adsorbed O2 (Table 1 and Supplementary Table 1).

Fig. 3

Experiment versus theory. a-f Comparison between the K-edge XANES experimental spectrum of Co0.5 (black hollow circles) and the theoretical spectrum calculated with the depicted structures (solid red lines). Cobalt, nitrogen, oxygen and carbon atoms are represented in purple, blue, red and gray, respectively

Table 1

XANES structural parameters for porphyrinic cobalt moieties

Row nbr.	Moiety	Co-N/Å	Co–O/Å	Bending/°	Rsq
	CoN 4 C 12
1	4-fold	1.96 (2)	—	—	1.12
2	5-fold	1.97 (3)	2.23 (5)	33 (6)	1.08
	CoN 3 C 10,porp
3	3-fold	1.96 (3)	—	—	1.16
4	4-fold	1.99 (3)	1.86 (4)	32 (7)	1.05
	CoN 2 C 5
5	2-fold	1.96 (3)	—	—	1.98
6	3-fold	2.01 (4)	1.90 (4)	50 (5)	1.13

Abbreviation: XANES, X-ray absorption near-edge structure

Best-fit structural parameters obtained from the analysis of the XANES spectrum of Co0.5 performed on the structures proposed in this work and depicted in Fig. 3. Bending is the angle between the Co–O vector and the O-O bond and Rsq is the residual function. Errors are given in parentheses

The energy of formation and O2 adsorption energy were then calculated by DFT-D approach for the possible sites identified by XANES analysis, namely CoN4C12, CoN3C10,por and CoN2C5, and the results were compared to those obtained for the CoN4C10 pyridinic moiety (Supplementary Fig. 6 ). The optimized bond distances for the above-mentioned porphyrinic sites, with and without a dioxygen molecule adsorbed in end-on mode, follow the trends of bond distances determined by EXAFS and XANES (Supplementary Table 3). The stability of these porphyrinic defective motifs is assessed by the Co(II) binding energy, −6.8 to −7.5 eV for the ground states, comparable to −8.2 eV obtained for CoN4C10 (Supplementary Table 4). Moreover, our DFT-D results reveal that the adsorption of an O2 molecule results in the displacement of the cobalt ion by 0.4 and 0.5 Å above the plane for CoN3C10,porp and CoN2C5, respectively. The O2 adsorption energies in the ground-state structures are −0.80, −1.23 and −1.26 eV for CoN4C12, CoN3C10,porp and CoN2C5, respectively. This is comparable to the O2 adsorption energies for the ground state of the CoN4C10 site, −0.97 eV. Although O2 adsorbs on all cobalt moieties, the obtained adsorption energies are, however, significantly smaller than those of the FeN4C10 and FeN4C12 moieties in ground state, −1.83 and −1.84 eV, respectively (Supporting Information of ref. [22]). This indicates that cobalt binds O2 more weakly compared to iron-based moieties. Due to the lower experimental ORR activity of Co0.5 vs. Fe0.5 (Fig. 2d) one can unambiguously conclude that cobalt moieties are situated on the weak-binding side of a volcano plot. While no definitive conclusion can be made on the volcano-plot position of iron moieties in Fe0.5 only from the present work (the weak and strong binding branches being both possible), yet unpublished work on Fe0.5 post-treated with H2O2 suggests that Fe0.5 is also positioned on the weak binding side of a volcano plot. The weaker binding of O2 (and thus of H2O2, via scaling relationships) by CoNC moieties also explains the higher amount of peroxide detected by RRDE for Co–N–C vs. Fe–N–C, as reported above. The O2 adsorption energy of molecular cobalt catalysts is often weaker than desired for an ideal metal-centered moiety,[42] and the same conclusion is made from the present combined experimental and theoretical study for the pyrolyzed Co0.5 catalyst. Hence, cobalt-moieties with the highest calculated O2 adsorption energy should be the most active ones, which gives the following expected ORR turnover frequency ranking CoN2C5 > CoN3C10,porp > CoN4C12. To better assess which moieties may be present in Co0.5, we measured the magnetic susceptibility of Co0.5 and N–C with superconducting quantum interference device (Supplementary Fig. 7). The plot of the inverse of the magnetic molar susceptibility χm of Co0.5 is linear with temperature (inset of Supplementary Fig. 7), as expected for Curie–Weiss paramagnetism. The positive intercept at 0 K leads to a negative Curie–Weiss temperature of −0.36 K, signifying weak antiferromagnetic interaction between magnetic moments[48]. The value of the slope 1/χ  = f(T) yields a magnetic moment µ eff = 3.52 Bohr magnetons (µ B), averaged on all cobalt moieties in Co0.5. For the first-row transition metals, to the first approximation, only the spin contribution to the effective magnetic moment can be considered (µ eff  = 2 µ B √(s2 + s)) leading to saverage = 1.33 in Co0.5, where saverage is the spin density averaged on all cobalt moieties present in Co0.5. According to the DFT spin density analysis and selecting the three best candidates on the basis of EXAFS and XANES analysis (Supplementary Table 4, CoN4C12, CoN3C10,porp and O2-CoN2C5 in their ground state, spin density of cobalt of 0.88, 0.58 and 1.83, respectively), one can propose that Co0.5 contains for example CoN4C12 and O2-CoN2C5 in a 53%/47% distribution, or CoN3C10,porp and O2-CoN2C5 in a 40%/60% distribution. Both distributions result in a theoretical saverage value of 1.33, similar to the experimental value (Supplementary Table 5). Considering only these three motifs, the fraction of O2-CoN2C5 sites must be in the range of 47–60%, the remainder being split between CoN4C12 and CoN3C10,porp sites.

We recently demonstrated that the porphyrinic FeN4C12 moiety with one or two axial ligands matched the XANES experimental spectrum of Fe0.5 and other Fe-based catalysts[22]. Since defective Fe-porphyrinic moieties were not investigated in that study, we extended here the XANES analysis of Fe0.5 in order to investigate whether such Fe-based moieties could also match the experimental XANES spectrum of Fe0.5. The analyses are shown in Supplementary Fig. 8 for FeN3C10,porp and FeN2C5. While O2-FeN3C10,porp reproduces the general features of the spectrum, the fit quality is poorer than that obtained with O2-FeN4C12 (Fig. 5e, f in ref. [22]). Although this does not preclude the existence of FeN3C10,porp moieties in Fe0.5, they may only represent a minor iron fraction. The present study therefore suggests that the structural disparity of MeNC moieties is larger in Co–N–C than Fe–N–C materials. Support for this is also provided by square-wave voltammetry (SWV) revealing a broader peak for Co0.5 (Supplementary Fig. 2).

Operando XANES signatures of Co- and Fe-moieties during ORR

We then investigated the behavior of Co- and Fe-based moieties in operando conditions. Figure 4a shows the normalized operando XANES spectra recorded at ORR potential range for Co0.5 in acid medium and Fig. 4b shows the same type of data for Fe0.5. Although the latter show a large change with electrochemical potential, similar to a previous report by Mukerjee’s group on other Fe catalysts[43], the magnitude of the change is insignificant for Co0.5 in the ORR region, even on the differential Δµ spectra (inset of Fig. 4a). The potential dependence of the XANES spectra observed on Fe0.5 in N2-saturated electrolyte supports the fact that the XANES variability is primarily controlled by the electrochemical potential, while the adsorption of O2 or oxygen intermediates has only a minor effect (Fig. 4d). The spectral changes observed in Fig. 4b probably result not only from a change in oxidation state, which would a priori only shift the edge position of the XANES spectra, but also from deeper structural changes and reorganization of the N (or C) ligands, as proposed in ref. [26] and/or an Fe(II) low-to-high spin-crossover[49]. SWV (Supplementary Fig. 2) shows that, while Fe0.5 displays a clear peak at 0.75 V vs. RHE in acidic medium, Co0.5 did not show any signal in the 0–1 V region, but instead shows a redox transition at 1.25 V vs. RHE. This largely explains the different potential-dependences of the XANES spectra recorded in the ORR region for Fe0.5 and Co0.5. In addition, the redox potential observed for Fe0.5 corresponds to its onset potential for the ORR (see Fig. 3b in ref. [22]). It is however surprising that the XANES spectra of Fe0.5 are still changing below 0.6 V vs. RHE, while most Fe ions should already be in + II oxidation state. Spin-state or conformation changes may still be occurring below 0.6 V.

Fig. 4

Operando XANES spectra. Taken in 0.5 M H2SO4 for Co0.5 (left) and Fe0.5 (right). The spectra were measured as a function of the electrochemical potential in N2-saturated electrolyte for a Co0.5 and b Fe0.5 (legend for potentials is in V vs RHE) and measured as a function of the saturating gas (O2 or N2) at either 0.2 V vs. RHE or 0.8 V vs. RHE for c Co0.5 and d Fe0.5. Insets of a and b are differential Δµ XANES spectra obtained by subtracting the normalized spectrum at every potential to the spectrum recorded at 0.2 V vs. RHE. Insets of c and d are Δµ spectra obtained by subtracting the normalized XANES spectra recorded in O2- to that recorded in N2-saturated electrolyte, at a fixed potential

The very different position of the redox peak for Fe0.5 and Co0.5 experimentally supports the fact that O2 adsorption on cobalt moieties is much weaker than on iron moieties. Linear relationships between the Me(III)/Me(II) redox potential and the adsorption energy for O2 have clearly been established for various macrocycles comprising first-row transition metals (Fig. 3 in ref. [42]). On this basis, the experimentally-determined relative positions of the redox peak for Fe0.5 and Co0.5 (0.75 and 1.25 V vs. RHE) are in line with our DFT calculations that predict O2 adsorption energies of ~ −1.8 eV for FeN4C12 and −0.8 to −1.2 eV for the porphyrinic Co-moieties.

Comparative operando XANES spectra collected under O2 and N2 provide additional insights in the reactivity of the active sites. Fig. 4c, d show the XANES spectra measured in O2- and N2-saturated acid electrolyte at two representative potentials (0.2 and 0.8 V vs. RHE) for Co0.5 and Fe0.5, respectively. Although for Fe0.5 we do not observe significant spectral changes in O2-free or O2-saturated electrolyte, for Co0.5 a clear variation is observed between 7720 and 7735 eV. This effect is magnified in the experimental ΔµO2-N2 spectra obtained by subtracting from the XANES spectrum of Co0.5 measured under O2 that measured under N2 (insets of Fig. 4c, d). A fit of the ΔµO2-N2 spectrum of Co0.5 at 0.8 V is shown in Supplementary Fig. 9, corresponding to an active-site structure with four nitrogen atoms at 1.95 Å with or without an oxygen molecule adsorbed end-on at 2.22 Å. In N2-saturated electrolyte, oxygen comes exclusively from water activation, whereas in O2-saturated electrolyte the oxygen may come from water activation or O2 adsorption. These Δµ analyses strongly suggest that, for Fe0.5, an Fe-O bond probably exists in both O2-saturated and N2-saturated electrolytes, with oxygen originating from O2 and H2O, respectively. In contrast, for Co0.5 the active sites are less oxophylic and a Co–O bond is formed only in O2-saturated electrolyte. To the best of our knowledge, this is the first direct operando evidence of molecular O2 adsorption on the cobalt centers in Co–N–C materials. The general picture that emerges from the present work supports the view that the active moieties formed at high temperature have a porphyrinic or defective porphyrinic architecture and also obey to the same general trends that have been reported on well-defined macrocycles comprising 3d transition-metal ions from the first row[22, 26, 42, 43]. However, while cobalt moieties in pyrolyzed materials do not change structure when catalyzing the ORR in acid medium, the iron moieties do change structure due to the change of oxidation and/or the spin state. The exact structural change induced by the redox switch will require extensive modeling work, which is beyond the scope of the present paper. Last, while general reactivity trends seem to apply for both pyrolyzed Metal–N–C materials and well-defined metal macrocycles, this does not mean however that their absolute ORR activity and O2 adsorption energy are the same. For example, the O2 adsorption energy on the three porphyrinic-derived CoNxCy moieties in ground state is −0.80 to −1.26 eV (Supplementary Table 4), higher than that reported for cobalt phthalocyanines and porphyrins (−0.60 to −0.35 eV).[42] The covalent integration of the Metal-N moieties in a conductive carbon matrix generally modifies the electron density at the metal relative to a metal macrocycle adsorbed on carbon.[50]

Electroactivity of cobalt moieties for other reactions

After having identified the porphyrinic structure of three possible cobalt moieties in Co0.5, we investigated their electrochemical activity for the OER, HER and HOR in acid medium. The cobalt moieties are the most active-site for OER in Co0.5, as indicated by the significantly lower OER activity of N–C (purple vs. Gray curve in Supplementary Fig. 10). Due to the low density of cobalt sites in Co0.5, the OER potential at a current density of 2 mA cm−2 is however positively shifted by 200 mV relative to that of unsupported IrO2 nanoparticles. While the OER activity of such cobalt moieties is interesting from a fundamental point of view, the low stability of carbon at such high electrochemical potential limits their practical application.[51] The low current observed at 1.3–1.6 V vs. RHE on Co0.5 might in fact entirely or partially originate from carbon corrosion, enhanced in the presence of cobalt relative to N-C. OER activity has been predicted from DFT theory for CoNC and FeNC moieties.[52] In contrast to the case for ORR, the oxygen-binding energy of the porphyrinic CoN4C12 moiety is close to the ideal value that may be expected for a Metal-N4 moiety.[52] OER activity has also been reported experimentally for heterogeneous Co–N[52]C materials in alkaline medium.[22] These heterogeneous catalysts mostly comprised metallic cobalt or cobalt oxide while the presence of CoNC moieties could not be certified.[53-55] The onset of OER in acidic medium on Co0.5 is not well defined but clearly situated above the redox peak associated with Co(II)/Co(III) at 1.25 V vs. RHE (Supplementary Fig. 2), suggesting Co(III)NxCy moieties could be the active-site for OER in Co0.5. The operando XANES spectra were also recorded at OER potentials, and these spectra showed only a small difference compared to the XANES spectrum recorded at 1.0 V vs. RHE (Supplementary Fig. 11), supporting the fact that CoNxCy moieties were still present at such high potential and are responsible for the OER activity of Co0.5. As for the ORR, it must be noted that the cobalt content in Co0.5 was not optimized to reach the highest activity.

In contrast, Co–N–C materials should be stable for long duration at low potentials needed to catalyze the HER. Two recent reports demonstrated the high activity of CoNC moieties toward HER in acid and alkaline medium.[28, 29] The HER polarization curve of Co0.5 is however negatively shifted by 200 mV relative to Pt/C (Supplementary Fig. 12). In addition, Co0.5 did not show any activity for the HOR (Supplementary Fig. 13), highlighting an irreversible behavior between HER and HOR on cobalt moieties. This intriguing fact seems to also apply to non-precious metal solid-state materials such as Mo and W nitrides, sulfides, and carbides. Such structures lead to efficient HER catalysts but were never reported to be efficient HOR catalysts.[56, 57]

Discussion

We identified with XANES, EXAFS, magnetic susceptibility measurements and DFT analysis the detailed structures and electronic state of single-atom catalytic sites in pyrolyzed Co–N–C materials. Three (defective) porphyrinic moieties were identified and showed some activity toward the ORR, OER and HER in acidic medium. DFT-D calculations showed that these cobalt moieties likely bind oxygen intermediates too weakly relative to an optimum Metal–N–C ORR catalyst, while iron moieties bind oxygen intermediates more strongly than cobalt ones. It however cannot be determined from the present study if iron moieties bind oxygen intermediates more or less strongly than optimally desired for the ORR. These insights support that the general activity descriptors for the active sites in pyrolyzed Metal–N–C catalysts are similar to those identified earlier for molecular Metal-N4 macrocycles. For Co–N–C catalysts, modifications of the carbon matrix with nitrogen or other light elements resulting in a stronger O2 binding energy at cobalt centers should improve their ORR activity. Generally, rational modification of the carbon matrix hosting such metal-centered moieties could be used for optimizing their activity toward a large variety of (electro)chemical reactions. These materials can bridge the gap between molecular and solid-state catalytic materials, offering perspectives for catalysts based on earth-abundant elements.

Methods

Synthesis

The catalyst precursor of Co0.5 was prepared from a Zn(II) zeolitic imidazolate framework (Basolite Z1200 from BASF, labeled ZIF-8), Co(II) acetate and 1,10-phenanthroline. Weighed amounts of the powders of Co(II)Ac (15.97 mg), phen (200 mg) and ZIF-8 (800 mg) were ball-milled (Pulverisette 7 premium, Fritsch) at 400 rpm for 2 h (4 × 30 min with 5 min rest in-between) in a ZrO2 crucible (45 cm3) filled with 100 ZrO2 balls (5 mm diameter). A split-hinge oven was equilibrated at 1,050 °C under Ar flow, and the catalyst precursor was quickly introduced and pyrolyzed at 1,050 °C in Ar for 1 h (see details in ref. [22]). The N–C material was prepared similarly but without cobalt. Due to a mass loss of 65–70 wt% during pyrolysis in Ar (unmodified by the presence of Co or Fe at 0.5 wt%) caused by volatile products formed from ZIF-8 and phen while Co does not form volatile compounds, the cobalt content in Co0.5 is 1.5 wt%.

Rotating ring disk electrode measurements

Electrochemical activity and selectivity toward the ORR was determined using the RRDE technique. For Co0.5, N–C and Pt/C (5 wt% Pt/C, provided by Johnson Matthey) catalysts, an ink including 10 mg of catalyst, 109 µl of a 5C wt% Nafion® solution containing 15–20% water, 300 µl of ethanol and 36 µl of de-ionized water was sonicated and mixed using a vortex. An aliquot of 9 µl was deposited on the glassy-carbon disk (0.247 cm2) with a Pt ring (Pine Instruments, nominal collection efficiency 0.37), resulting in a loading of 818 µg·cm−2. For lower Co0.5 loadings (200 µg cm−2), the catalyst mass and Nafion solution aliquot in the ink were proportionally decreased, and the decrease in Nafion solution balanced by increased ethanol aliquot. 9 µl of ink was again deposited on glassy carbon. The working electrode with the deposited catalyst layer was used in a four-electrode cell setup connected to a bipotentiostat (Biologic SP 300) and rotator (MSR, Pine Instruments). The counter and reference electrodes were a graphite rod and a reversible hydrogen electrode (RHE), respectively. The acidic electrolyte was an O2-saturated aqueous solution of 0.1 M H2SO4, except for ORR on Pt/C (0.1 M HClO4). For Pt/C only, the electrode was first cleaned by performing 300 cycles at 500 mV s−1 between 0 and 1 V vs. RHE. The RDE polarization curves were recorded with a scan rate of 10 (Fe–N–C) or 50 mV s−1 (Pt/C, to avoid contamination) at 1,600 r.p.m. and corrected for the background current measured in N2-saturated electrolyte. The second cycle was used for correction. For measuring the % H2O2 released during ORR, the potential of the Pt ring was held at 1.2 V vs. RHE.

Square-wave voltammetry

Catalysts inks were prepared by dispersing 3.1 mg catalyst in a solution prepared by mixing 150.2 µL of millipore water, 465 µL of isopropyl alcohol and 6.2 µl of 5 wt% Nafion® (Nafion to catalyst mass ratio of 10 wt%). The ink solution was then sonicated 60 min. 20 μl of the ink was pipetted on the glassy carbon disk (0.247 cm2) to reach a loading of 400 μg cm−2. This lower loading than used for most RRDE experiments was necessary to obtain a proper balance between the signal coming from the redox peak and that due to the double layer. The electrodes were dried through a rotational drying method at 800 r.p.m. SWV experiments were carried out in 0.1 M HClO4 aqueous electrolyte using an Autolab bipotentiostat (PGSTAT302N) with a step potential of 5 mV, potential amplitude of 20 mV and scan frequency of 10 Hz in a standard electrochemical cell (Chemglass).

Fuel cell measurements

For the membrane electrode assembly, cathode inks were prepared using the following formulation: 20 mg of catalyst, 652 μl of a 5 wt.% Nafion® solution containing 15–20% water, 326 μl of ethanol and 272 μl of de-ionized water. The inks were alternatively sonicated and agitated with a vortex mixer every 15 min, for a total of 1 h. Then, three aliquots of 405 μl of the catalyst ink were successively deposited on the microporous layer of an uncatalysed 4.84 cm2 gas diffusion layer (Sigracet S10-BC) to reach a catalyst loading of 4 mg cm−2. The cathode was then placed in a vacuum oven at 80 °C to dry for 2 h. The anode used for all PEM fuel cell tests performed in this work was 0.5 mgPt·cm−2 on Sigracet S10-BC. Assemblies were prepared by hot-pressing 4.84 cm2 anode and cathode against either side of a Nafion NRE-211 membrane at 135 °C for 2 min. PEMFC tests were performed with a single-cell fuel cell with serpentine flow field (Fuel Cell Technologies Inc.) using an in-house fuel cell bench and a Biologic Potentiostat with a 50 A load and EC-Lab software. For the tests, the fuel cell temperature was 80 °C, the humidifiers were set at 85 °C, and the inlet pressures were set to 1 bar gauge for both anode and cathode sides. The flow rates for humidified H2 and O2 were controlled downstream of the polymer electrolyte membrane fuel cell (PEMFC) and fixed at 50 s.c.c.m. Polarization curves were recorded by scanning the cell voltage at 0.5 mV s−1.

Spectroscopic experimental characterization

Co and Fe K-edge X-ray absorption spectra were collected at room temperature at SAMBA beamline (Synchrotron SOLEIL) equipped with a double crystal Si 220 monochromator. Catalyst inks were prepared by mixing 10 mg catalyst with 50 µl de-ionized water and 100 µL of 5 wt% Nafion® solution with ultrasounds. A 50 µl aliquot was then pipetted on ~3 cm2 circular area of a 100-µm-thick graphite foil (Goodfellow cat. C 000200/2), resulting in a catalyst loading of ~1 mg cm−2. The graphite foil then served as a working electrode, and was installed in an electrochemical cell (PECC2, from Zahner), see Supplementary Fig. 14. The cell also includes a reference Ag/AgCl electrode, and a Pt counter electrode. The cell was filled with a given electrolyte and saturated with either O2 or N2 by continuously bubbling gas in the electrolyte. Operando measurements were performed by recording the Kα X-ray fluorescence of the element under study (Co or Fe) with a Canberra 35-elements monolithic planar Ge pixel array detector, while the ex situ spectra of the reference compound Co(II)Pc and Co0.5 were recorded in transmission geometry on pelletized disks of 10 mm diameter with 1 mm thickness, using Teflon powder (1 µm particle size) as a binder.

Spectroscopic analysis and modeling

The EXAFS data analysis was performed with the GNXAS code[58, 59]. In that approach the interpretation of the experimental data is based on the decomposition of the EXAFS χ(k) signal (defined as the oscillation with respect to the atomic background cross-section normalized to the corresponding K-edge channel cross-section) into a summation over n-body distribution functions γ (n) calculated by means of the multiple-scattering (MS) theory. The cobalt first coordination shells have been modeled with Γ-like distribution functions which depend on four parameters, namely, the coordination number N, the average distance R, the mean-square variation σ 2 and the skewness β. Note that β is related to the third cumulant C3 through the relation C3 = σ 3 β. The standard deviations given for the refined parameters in Supplementary Table 1 are obtained from k 2-weighted least-squares refinements of the EXAFS function χ(k), and do not include systematic errors of the measurements. Least-square fits of the EXAFS raw experimental data have been performed by minimizing a residual function of the type:where N is the number of experimental points, E {λ} = (λ , λ ,…, λ ) are the p parameters to be refined and σ 2 is the variance associated with each experimental point α exp(E ). Additional non-structural parameters were minimized, namely E 0 (core ionization threshold energy) and S 0 2 (amplitude reduction factor taking into account intrinsic losses).

The XANES data analysis of Co0.5 was carried out with the MXAN code in the framework of the full MS scheme, following the same modeling approach applied to identify the structure of the catalytic active sites in Fe–N–C materials[22]. The MXAN method is based on the MT approximation for the shape of the potential and uses a complex optical potential, based on the local density approximation of the self-energy of the excited photoelectron[60]. The minimization of the Co(II)Pc XANES spectrum has been carried out starting from the X-ray structure of this compound[47], while in the case of Co0.5 we tested different coordination geometries. The cluster size used in the calculations was chosen on the basis of a convergence criterion. The fit includes a minimal number of selected parameters: the Co–N and Co–O distances, the Co-displacement, the Co-5th ligand distance and the bending angle between the Co–O vector and the O–O bond. During the fit, the outer carbon atoms of the ring rigidly followed the motion of the nitrogen atoms. Least-square fits of the experimental data in the space of the structural parameters were achieved by minimizing the residual function defined aswhere n is the number of independent parameters, m is the number of data points, y th and y exp are the theoretical and experimental values of absorption, respectively, ε is the individual error in the experimental data set, and wi is a statistical weight. Here, we assumed a constant experimental error, ε = 1.2%, for the whole experimental data set. Five non-structural parameters have been optimized, namely the Fermi energy level E F, the experimental resolution Γexp, the threshold energy E 0 and energy and amplitude of the plasmon, E s and A s. The fit of the Δµ XANES spectrum of Co0.5 in O2- or N2-saturated electrolyte was obtained by the MXAN code with the following procedure: (1) Theoretical XANES spectrum calculation of the CoN4C12 moiety as reference, using the spectrum under nitrogen at 0.8 V. (2) At each step of the fitting procedure the program calculated the theoretical spectrum of a CoN4C12-O2 moiety under oxygen. (3) The difference spectrum is calculated.

DFT-D computation

DFT-D computations were carried out on cluster models shown in Supplementary Fig. 6. The dangling bonds were saturated with hydrogen. The computational method is based on the DFT cluster approach augmented with an empirical dispersion term as implemented in deMon2k computer program used for the present computations[61]. For the exchange-correlation functional, the approximations of Perdew–Burke–Ernzerhof’s 1996 was used[62]. Triple-ζ basis sets were used for the C and H atoms and double-ζ bases plus polarization for Co(II) cations[63]. Automatically generated auxiliary functions up to l = 2 (for the metal atom) and 3 (for H, C and O atoms) were used for fitting the density. A quasi-Newton method in internal redundant coordinates with analytical energy gradients was used for structure optimization. The convergence was based on the Cartesian gradient and displacement vectors with thresholds of 10−3 a.u. The energy convergence was set to 10−7 a.u. The binding energy of cobalt cation per Co–N bond (BE) and adsorption energy of end-on adsorbed O2 molecule (E ads) are obtained as

In the above formula, the total energies for the ground state triplet dioxygen with total spin S = 1 (two unpaired electrons) and Co2+ in its quartet ground state (total spin S = 3/2) were used. The relative energies ΔE with spin number are obtained as a difference between the total energies of the higher-energy minima states and the ground state. The atomic spin density was computed for all the minima energy structures using Mulliken[64] and Hirshfeld[65] population analysis. Both population schemes led to nearly identical results. Cobalt divalent cations were used. The relative energies ΔE with spin number are obtained as a difference between the total energies of the higher-energy minima states and the ground state.

SQUID measurements

The molar magnetic susceptibility χ m was measured from 2 to 300 K for Co0.5 and N–C with a Superconducting QUantum Interference Device (SQUID) (MPMS XL-7T, Quantum Design) at a magnetic field of 5,000 Oe. A mass of 16 mg of Co0.5 or N–C was weighed and introduced in a polymer straw. The average effective magnetic moment of cobalt atoms (µ eff) was then obtained by fitting the plot of 1/χ (in mol of cobalt atoms per emu) vs. 1/T with a linear law in the region 15–77 K. The value for the slope obtained is defined as a 1/C . From the C -value, the average effective magnetic moment of cobalt atoms (µ eff) is directly calculated via the relation µ eff  = 2.82·C 1/2, in units of Bohr magneton (µ B). The average spin of cobalt ions is then obtained via µ S = g·µ B·sqrt(s(s + 1)), where g is equal to 2.002 and s is the average spin density for all cobalt moieties present in Co0.5.

TEM measurements

For transmission electron microscopy (TEM), the sample powders were dispersed in ethanol using an ultrasonic bath and deposited onto holey-carbon grids. TEM imaging and SAED were carried out on a JEM-2100HCKM microscope (JEOL Ltd., Japan) operating at 120 keV.

Data availability

The data that support the findings of this study are available from the corresponding authors upon reasonable request. The X-ray absorption Spectroscopy raw data associated with this work is permanently stored at SOLEIL and available upon request.