Sakshi
Bhardwaj
a,
Samadhan
Kapse
b,
Soirik
Dan
c,
Ranjit
Thapa
*b and
Ramendra Sundar
Dey
*a
aInstitute of Nano Science and Technology (INST), Sector-81, Mohali-140306, Punjab, India. E-mail: rsdey@inst.ac.in
bDepartment of Physics, SRM University, Andhra Pradesh 522240, India
cIndian Association for the Cultivation of Science, Kolkata 700032, India
First published on 7th April 2023
The active sites of electrocatalysts play a crucial role in the material design and mechanistic exploration of an electrocatalytic reaction. Defect-tailored heteroatom-doped carbon-based electrocatalysts for oxygen reduction reaction (ORR) have been much explored, but there is ambiguity in the prediction of active sites responsible for the performance of the material. To find the origin of the activity of this class of catalysts towards ORR, in this work, we use the quantum mechanics/machine learning (QM/ML) approach to derive energy-optimized models with both N-atoms and 5-8-5 defect sites which manifest exemplary ORR activity. Following this approach, we synthesized defect-engineered graphene (DG) using the zinc template method at 1050 °C to achieve optimum N-dopants and intrinsic (5-8-5) defects. The obtained electrocatalyst exhibits hierarchical porosity, high surface area, low nitrogen content, good stability and a satisfying ORR performance with a half-wave potential (E1/2) of 0.82 V, comparable to that of commercial Pt/C (E1/2 = 0.82 V). Further, the full energy profile was deduced using density functional theory and the charge redistribution in the material cross-verified a reduced overpotential for ORR. This work provides a strategy for the synthesis of noble-metal-free high-performance electrocatalysts for energy conversion.
Pure graphene is non-polar, lacks active sites and shows poor electrochemical activity.11 Therefore, numerous strategies have been adopted to improve its electrocatalytic properties.11–13 Defect engineering in carbon nanomaterials is considered one of the most effective methods.14–17 According to the second law of thermodynamics,18,19 crystalline materials will always have lattice imperfections or disorder arising from the production process or impurities. The presence of atomic disorder or defects is the origin of the electronic, optical and mechanical properties of carbonaceous matrices. It is well known that the presence of heteroatoms including N, P, S, B and O (extrinsic defects) and topological defects (intrinsic defects) in nanocarbon materials is responsible for the disturbance in the local electronic structure that influences the charge distribution in the framework and results in improved electrocatalytic performance of the carbonaceous materials towards ORR.20,21 Moreover, the presence of N and defects in the sp2 framework cause alteration in the charge and spin distribution of the carbon framework, which facilitates the adsorption of oxygen molecules, reduces the energy barrier, subsequently leads to the breaking of the O–O bond and improves the overall reaction rate.22,23 Nonetheless, there is still controversy in the literature questioning which of graphitic-N,23,24 pyridinic-N25,26 or topological defects is the true active site for ORR when they all are present in carbonaceous matrices. Nitrogen-doped nanocarbon materials have gained massive attention, with various studies showing good electrocatalysis,27,28 after the report on nitrogen-doped carbon nanotube arrays for ORR from Dai's group.29 Intrinsic defects have been ignored for a long time during research progress into the ORR mechanism, but have now become a point of focus due to their significant contribution towards advanced ORR activity.30,31 Most literature has some solution for choosing active sites, but the design of the active sites is based on proof-of-concept research. Therefore, an efficient principle is required for the screening of carbon catalysts. For this, descriptors in combination with machine learning (ML) algorithms is a good approach to narrow down the ideal models for catalysts.32–35 Then, in the selected models, a predictive model can be used to find those with a higher number of active sites.31 This approach helps to narrow down the search for a catalyst and suggest the best possible active material to synthesize in the laboratory.
In this work, we have considered distinct types of models with N-content and defects. In a previous study, the 5-8-5 defect proved to be more active for ORR. Our main concern is to find the active sites within models having both N-content and the 5-8-5 defect. Using the QM/ML approach, we screened the models and narrowed the choices down to three models (I, II and III) with optimum N-content and 5-8-5 defects (Fig. 1a–c) and their active sites were determined for efficient ORR. Following this theoretical investigation, we have synthesized a N-doped and 5-8-5 defect-enriched graphene-like material (DG-T) by pyrolysis of C and N sources in the presence of a porogen at three different temperatures. The catalyst prepared at 1050 °C (DG-1050) possesses an optimal amount of defects and N-content. The presence of N-dopants and defects was evidenced by XPS and STM studies. The synergistic impact of both these defects is shown in the E1/2 potential for ORR of 0.82 V, very close to that of commercial Pt/C (0.829 V); it also achieved excellent stability. In both the theoretical and experimental results, DG-1050 is found to be most active towards ORR with the maximum number of active sites. Next, a full energy pathway was done for this material. We found that N present in the vicinity of an intrinsic defect (5-8-5) synergistically plays a critical role in the activity origin of the nanocarbon catalyst. This work sheds light on the method of developing metal-free carbon-based electrocatalysts by tuning the content of defects and N-dopants for their efficient usage in electrocatalytic applications.
This QM/ML approach is faster, more efficient and can help analyze any type of graphene-based configuration for ORR without performing adsorption calculations. Previous reports40,41 suggest that machine learning in combination with DFT can accurately find the sites which are responsible for the activity. However, until now there have been no reports, to the best of our knowledge, wherein machine learning and DFT were performed on defect-enriched N-doped graphene to find the true active sites for ORR. Nitrogen doped in carbon systems has already been reported as an active dopant for ORR. Through the QM/ML approach, herein, we discover the specific roles of N-doping and defects towards the formation of active sites for ORR (Fig. S1†). Among them, the 5-8-5 defect more effectively decreases the charge neutrality to increase active sites on the graphene surface.42 Based on the previous reports and our study on N-doping and the 5-8-5 defect, we modeled three realistic, random complex model structures termed systems I, II and III, as seen in Fig. 1. Model structure I is made with two pyrrolic-N, two pyridinic-N, two oxidized-N and two graphitic-N sites; model structure II is made with one pyrrolic-N, two pyridinic-N, one oxidized-N, one graphitic-N and a 5-8-5 defect; model structure III is made with one pyrrolic-N, one pyridinic-N, one oxidized-N and a 5-8-5 defect. For each system, the possible carbon active sites around the N-dopants and defects were selected for ORR study. Different types of energy, electronic and structural descriptors have been reported to define the ORR activity in carbon materials.43–49 Most of these descriptors have limitations, such as lower efficiency, inability to predict the site-specific activity, etc. In this context, the change in Gibbs free energy of OH adsorption (ΔGOH) is a well-known energy descriptor that effectively predicts site-specific ORR activity.4,50 From the volcano correlation between ΔGOH and thermodynamic overpotential (η), a value of ΔGOH near zero indicates a smaller η and higher ORR activity.75 Recently, the role of the π orbital of an active site in a carbon system has also been identified as an origin with electronic descriptors to predict site-specific catalytic activity for ORR.51 Further, these π-orbital based descriptors (Dπ(EF), R-Oπ) have been used as features in machine learning algorithms to develop a highly cost effective (less computational power) predictive model (ΔGO − ΔGOH) OER application.52 As for ORR, these electronic descriptors have been utilized in various machine learning algorithms to develop a predictive model of ΔGOH.39 Here, the predictive model of ΔGOH is also calculated by using the simple linear fit method (Fig. S2†). The calculations for the simple linear fit method are given in Section S1a.† Among all the models, the SVR predictive model showed optimal performance and a higher test R2 score (0.87) in 5-fold cross-validation (Fig. S3, Sections S1b and c†). This SVR predictive model of ΔGOH is faster and more cost effective, requiring only a single DFT calculation for a carbon system to predict the ORR activity of all the sites (the Python code script is given in Section S1d†). Without the SVR model, we would need to perform three DFT calculations for O, OH and OOH intermediate adsorption on each active site to compute the ORR performance of each site. Therefore, we used this reported SVR predictive model to analyze the active sites in systems I, II and III for ORR.
Fig. 1a represents the 22 different possible carbon active sites of model I. To find the theoretical overpotentials on 22 sites, we would need to do 66 DFT calculations (three intermediates for each site), which require resources and time. But, by applying our SVR model, we can predict the ΔGOH values and overpotential for these active sites by doing single DFT calculations of the host surface only. The story is the same for the other two models. So, the predictive model developed using the QM/ML approach and fewer DFT calculations helped us to find the best active site in a very efficient way. We found that in model I, the carbon site denoted as number 9 in Fig. 1a adjacent to the oxidized nitrogen shows the lowest ΔGOH (0.53) and is the most active site for ORR in this model structure53 (Table S1†), whereas the ΔGOH values for seven other sites are in the range of 0.6 eV to 0.8 eV and fourteen sites show higher ΔGOH values (>0.8 eV). As shown in Fig. 1b, model II has a 5-8-5 defect with a lower N-doping concentration compared to model I. Here, we observed that the combined effect of the 5-8-5 defect and N doping decreases the charge neutrality on graphene and creates more active sites with lower ΔGOH values (<0.6 eV) for ORR (Table S2†). In model II, we found seven ideal active sites for ORR near 5-8-5 defects. In the case of the model III defective graphene (Fig. 1c), it is noted that the further decrease of N-doping percentage reduces the number of ideal active sites near 5-8-5 defects (Table S3†). In Fig. 1d, for easy understanding, we display a histogram using the percentage of active sites with different ΔGOH ranges in models I, II and III. Therefore, we predict that the 5-8-5 defect and N doping play a combined role to improve the ORR performance in graphene systems.
Following the theoretical study, we synthesised the samples DGx:y-T (T = pyrolysis temperature, x:y = the ratio of melamine and zinc catalysts), as shown schematically in Fig. 2a. The experimental details are mentioned in the Experimental section. Briefly, DG was prepared by the pyrolysis of melamine and zinc dust in a particular ratio. Due to its alternating carbon and nitrogen atoms and powerful self-condensation capacity, melamine was selected as the carbon and nitrogen source to prepare DG. In this synthesis procedure, Zn not only serves as a template, but also acts as a porogen to construct a highly cross-linked 3D porous network structure, since Zn can be evaporated >907 °C during the carbonization process.22,54 During the carbonization step at 800 °C, metallic Zn embeds into the carbon matrix and leads to the expansion of carbon layer distance. Further heating was performed at different temperatures (950, 1050, 1150 °C) to remove zinc and nitrogen to create defects.22,55 Subsequently, the carbon atoms rearrange and self-organize themselves to form a carbon framework. Using only the carbon source dextrose, a defective carbon sample (C1:1) was also prepared to check the effect of N on the ORR activity.
Fig. 2 (a) Schematic illustration of synthesis of DG1:1-1050; (b) SEM image; (c and d) HRTEM images; (e) AFM image, inset: height profile graph; (f and g) STM images of DG1:1-1050. |
Field emission scanning electron microscopy (FESEM) and transmission electron microscopy (TEM) were used to analyse the morphology of the prepared samples. FESEM images (Fig. 2b) show a highly porous interconnected 3D network of graphene sheets that offers high specific surface area with a greater number of active sites to help in easy mass transport. The formation of flake-like wrinkled defective graphene sheets after carbonization at high temperature was confirmed by TEM (Fig. S4†). The edge defects produced by the removal of Zn56,57 can be clearly seen in the high-resolution TEM (HRTEM) images (Fig. 2c and d). No usual graphitic lattices were found, suggesting the disordered and defect-rich character of DG1:1-1050. The associated selected area electron diffraction (SAED) pattern (inset Fig. S4†) displayed the two typical rings, (002) and (101), that correspond to the amorphous nature of DG1:1-1050. Atomic force microscopy (AFM) evidenced the formation of ultrathin graphene sheets (Fig. 2e inset; height profile: ∼3–4 nm) with localised destruction due to the removal of N and Zn atoms from the framework (Fig. 2e). The created intrinsic defects can be pentagons, heptagons, even octagons and holes etc., formed by the rearrangement of atoms due to multiple single-atom or diatom vacancies created after the removal of N atoms.58 Scanning tunnelling microscopy (STM) was used to visualize the defects formed in the synthesised material. Fig. 2f shows the low magnification image presenting the DG1:1-1050 sheet; its highly resolved image (Fig. 2g) obtained at 1 nm clearly presents the various structural defects such as pentagons and 5-8-5 defects. These defects cause the local modulation of the electronic environment responsible for enhancing the electrocatalytic activity.
The X-ray diffraction (XRD) patterns for the different DG-Ts (Fig. 3a) exhibited one broad characteristic peak at 24.8° and another small hump around 44° indexed to the (002) and (101) planes of nanocarbon, respectively, suggesting an amorphous structure with a low degree of graphitization that is expected to be beneficial for the electrocatalytic activity of the catalyst. Raman analysis is the most widely used tool to identify the degree of graphitization. The introduction of defects in carbon nanomaterials in the form of heteroatom doping or vacancies gives rise to an increase in the D-band (1337 cm−1) intensity due to disruption of the hexagonal rings.59,60 The G-band at 1569 cm−1 is attributed to the C–C and C–N bond vibrations of E2g mode in sp2-bonded graphitic carbon. These two characteristic peaks were obtained in all three DG1:1-T samples (Fig. 3b). The intensity ratio of the D-band to the G-band (ID/IG) estimates the defect level present in carbon materials.15,61 The ID/IG ratios were 1.05, 1.11, and 1.18 for DG1:1-950, DG1:1-1050 and DG1:1-1150, respectively, showing that all the materials contained a high degree of carbon disorder. There was substantial decrease in N-content with increasing temperature and, thus, the surge in ID/IG ratio distinctly proves the extensive formation of defective domains. The abundant edge and topological defects along with the optimal N-content attuned the surface and electronic properties of DG1:1-1050 that regulate the adsorption energies of electrocatalysis steps. The defect density present in the sp2 domains of the graphene rose from 3.14 × 1011 to 3.54 × 1011 with increasing temperature, as calculated using the formula
(1) |
Although Raman spectroscopy helps define the relationship between defects and ORR performance, it does not provide much information about the C-speciation related to the D-band feature, such as the hybridization of the carbon atoms in the ring (sp3 and sp2). Therefore, near edge X-ray absorption fine structure (NEXAFS) measurements were carried out to investigate the local chemical bonding and electronic structure of the carbon-defect active sites. The C-K edge spectra obtained for all DG1:1-T samples were found to be similar (Fig. 3c), indicating the similar features of their carbon skeletons. The π* region is defined as up to ∼290 eV and the photon energy range above 290 eV is the σ* region.63 The peak at 285.5 eV was assigned to the presence of low coordination carbon atoms at the edges of DG1:1-T that led to the rehybridization and breakage of integrity of the π-conjugation.64,65 Therefore, the intensity of the π* feature decreased as temperature rose and electron transfer to the O2 molecule is facilitated by the active unpaired π-electrons located at each edge carbon atom. The excitation at 286.2 eV is attributed to N-substituted aromatic carbon and that at 287–288 eV is attributed to C–O/C–OH/CO.66
The specific surface area and porosity of DG1:1-T were determined by N2 adsorption and desorption isotherm measurements (Fig. 3d). The Brunauer–Emmett–Teller (BET) specific surface areas for the DG1:1-950, -1050 and -1150 samples were found to be 580, 822 and 831 cm2 g−1, respectively, which shows that the specific surface area increased with an increase in carbonization temperature. The obtained type I isotherm followed by a type IV isotherm with hysteresis suggested the formation of an inter-connected hierarchical porous graphene-like structure spanning both micro- and macropores, further confirmed by the pore size distribution graphs (inset Fig. 3d).67 The high BET surface area gives a high probability for active site exposure and the hierarchical porous structure can participate more efficiently in the formation of a triple-phase boundary that provides less-resistant diffusion channels for the facile movement of ORR reaction species with the electrolyte and enhances their interactions with catalytically active sites for ORR.67,68
That the nitrogen loss by heat treatment is responsible for the defects in the graphene sheets was confirmed by X-ray photoelectron spectroscopy (XPS), as it gives the chemical composition and bonding configuration of the catalyst. The XPS survey spectra (Fig. 4a) confirmed the presence of C, N and O elements in the DG1:1-T. The fitted high-resolution C 1s spectra indicated residual N and O atoms bonded to the C atoms after annealing, with peaks at about 284.7, 286.07 and 287.8 eV corresponding to C–C, C–N and CO, respectively, as presented in Fig. S5.†69 The successful removal of N atoms from the carbon skeleton by annealing was verified by the comparison of deconvoluted N 1s spectra of DG1:1-T that depict the characteristic peaks for pyridinic N (398.0 eV), pyrrolic N (399.9 eV), graphitic N (401.0 eV) and oxidized N (403.5 eV), as shown in Fig. 4c.70 However, on annealing at higher temperatures, the signals for N became weak in the XPS pattern.71 It is clearly represented in the bar diagram (Fig. 4b) that, with the increase in temperature, the N-content decreased from 12.63% to 2.91% and then to 1.36%, resulting in the rise in the formation of defects.72,73 The comparison bar graph of N-content indicates that the defects are mainly created by the removal of pyridinic and graphitic nitrogen, giving the 5-8-5 or pentagon defects (clearly seen in the STM image) that are highly favorable for ORR activity.74 Three models have been made based on the N-content obtained from the XPS analysis, as shown in Fig. 4d, and they are similar to theoretical models I, II and III.
It is considered that N dopants, particularly graphitic24,75 and pyridinic nitrogen,76,77 can act as active centers for ORR. However, Yao et al. presented in their studies that defects are more functional than N in electrocatalytic activity.22,32,78 Here, we suggest that N and a nearby defect can work synergistically to promote the active sites and the number of active sites also increases due to existence of both extrinsic and intrinsic defects. Optimal concentrations of defects and nitrogen content are required for excellent ORR activity. The presence of N species in the carbon skeleton helps to improve its charge mobility, due to the introduction of electron-donor characteristics, and raises the catalytic activity of carbon in electron-transfer reactions and the abundance of intrinsic defects in the graphene sheet attunes the π electron density in sp2 carbon systems, providing active sites for ORR.
ϕ = hν − (Ecut-off energy − EF) | (2) |
From the SVR models, we also found maximum ideal sites for model II. The reaction kinetics were studied using a rotating ring disk electrode (RRDE) under similar conditions by changing the rotation speed of the electrode from 625 to 4900 rpm and using the Koutecky–Levich (K–L) equation (Fig. 5d). The diffusion limiting current density increased with an increase in the speed of rotation due to more diffusion of oxygen and a shortening of the diffusion layer at the electrode surface.80 Therefore, the obtained K–L plot (inset Fig. 5d) at different potentials (0.43 to 0.63 V) showed linearity and near parallelism that accounted for the first order kinetics of the concentration of dissolved oxygen and the similar number of transferred electrons (n; eqn (4)†) over the whole potential range during the electroreduction of O2 (Fig. 5e), representing the superior 4 electron transfer in the ORR process. It is noticeable that the catalyst yielded 4 OH−, which means, following the 4e− transfer kinetics, a negligible number of unwanted intermediates must have formed. The production of H2O2 (2e− process) is an important concern as its decomposition can lead to a chain oxidation reaction due to the formation of highly reactive intermediates.81 The H2O2 yield can be calculated by measuring the ring current using a RRDE system (eqn (5)†). The reduction of peroxide species occurring in the ORR process can be seen from the ring current profile shown in Fig. S11.† The H2O2 produced over the potential range is very low (Fig. 5e). The superiority of the activity of DG1:1-1050 compared to other control samples and comparability to commercial Pt/C was reinforced by the mass transport corrected Tafel plot drawn from the low overpotential region of LSV at 1600 rpm at a 2 mV s−1 scan rate. The Tafel slope of DG1:1-1050 was 52 mV dec−1, better than that of Pt/C (64 mV dec−1), as shown in Fig. 5f, presenting fast ORR kinetics with a high electron transfer coefficient and fast first electron transfer from the catalytic site. The low overpotential, high catalytic activity and four-electron pathway were supported by the theoretical calculations. Also, for validation of the SVR model prediction of ΔGOH, we considered the most ideal carbon site number 21 (with predicted ΔGOH = 0.36 eV) in DG1:1-1050 (Fig. 6a) and performed DFT calculations for intermediate O, OH, and OOH adsorption. By plotting the free energy profile of ORR (Section S2†) for carbon site 21, we confirmed its lower overpotential and higher catalytic performance (Fig. 6b). The first step of OOH* formation is found to be the potential determining step with an overpotential of 0.45 V at equilibrium potential (U = 0.4 V). We also studied the possibility of a two-electron pathway for ORR on site 21 (Section S1e†). In this context, we show that the four-electron pathway is more favorable over the two electron pathway of ORR due to the higher free energy of free OOH/OHOH ion formation than of the O* intermediate (Fig. S12†).
Further, we computed the density of states to investigate the adsorption of ORR intermediates on ideal active site 21 (Fig. S13†). From this, we observed a reduction in the π orbital density of states near the Fermi level after adsorption of O, OH, or OOH, whereas oxygen adsorption occurs through an enolate configuration on site 21 with a large decrease in the density of states near the Fermi level. Later, we confirmed the charge transfer between the active site and the ORR intermediate ions through a Lowdin charge analysis (Fig. S14†). In summary, we unveiled the combined effect of the 5-8-5 defect with N doping and the origin of the higher ORR activity of DG1:1-1050 by using π orbital descriptors and the SVR model. The overall mechanism of ORR on active site 21 in DG1:1-1050 is shown in Fig. S15.†
A useful catalyst must have the potential to prevent the possible methanol crossover effect that is a major concern for fuel cell catalysts. To investigate the effect of methanol on our catalyst and Pt/C, the chronoamperometric response was checked after the addition of 1 M CH3OH in 0.1 M KOH. As seen in Fig. 7b, there occurred almost no change in current density for 850 s for our catalyst, while in the case of Pt/C, an instantaneous increase in current density from negative to positive current was observed after the addition of methanol. The LSVs of DG1:1-1050 for ORR taken before and after methanol addition exactly overlaid with each other (inset Fig. 7b). Therefore, the stability of DG1:1-1050 in contrast to that of Pt/C against methanol crossover presented the feasibility of its usage as a cathode material in both alkaline and direct methanol fuel cells.
For comparison, all the defective N-doped carbon based electrocatalysts for ORR are summarized in Table S4† and a bar graph in Fig. 7c. DG1:1-1050 is found to show similar ORR activity amongst the reported non-metal carbonaceous materials.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d3ta00871a |
This journal is © The Royal Society of Chemistry 2023 |