The Oxygen Reduction Reaction on Nitrogen-Doped Graphene

Non-PGM Electrocatalysts for PEM Fuel Cells: A DFT Study on the Effects of Fluorination of FeNx-Doped and N-Doped Carbon Catalysts

Fluorination is considered as a means of reducing the degradation of Fe/N/C, a highly active FeN_x-doped disorganized carbon catalyst for the oxygen reduction reaction (ORR) in PEM fuel cells. Our recent experiments have, however, revealed that fluorination poisons the FeN_x moiety of the Fe/N/C catalytic site, considerably reducing the activity of the resulting catalyst to that of carbon only doped with nitrogen. Using the density functional theory (DFT), we clarify in this work the mechanisms by which fluorine interacts with the catalyst. We studied 10 possible FeN_x site configurations as well as 2 metal-free sites in the absence or presence of fluorine molecules and atoms. When the FeN_x moiety is located on a single graphene layer accessible on both sides, we found that fluorine binds strongly to Fe but that two F atoms, one on each side of the FeN_x plane, are necessary to completely inhibit the catalytic activity of the FeN_x sites. When considering the more realistic model of a stack of graphene layers, only one F atom is needed to poison the FeN_x moiety on the top layer since ORR hardly takes place between carbon layers. We also found that metal-free catalytic N-sites are immune to poisoning by fluorination, in accordance with our experiments. Finally, we explain how most of the catalytic activity can be recovered by heating to 900 °C after fluorination. This research helps to clarify the role of metallic sites compared to non-metallic ones upon the fluorination of FeN_x-doped disorganized carbon catalysts.

Boron Coordination Effect in Ni-Nx Doped Graphene Catalysts on the ORR Performance Based on DFT Calculations

The coordination atoms of metal active site in transition metal N-doped carbon single atom electrocatalysts play a vital role in dominating the catalytic performance of oxygen reduction reaction (ORR) at the cathode of fuel cells or metal-air cells. In view of weak adsorption ability of Ni active site in NiN₄ -C catalysts to oxygen intermediate states, herein we introduce boron atoms with smaller electronegativity than N and C atoms to modulate the local coordination environment and electronic structures of Ni site. First-principles density functional calculations reveal that both B substitution for N atoms (NiN₂ B₂ -C) and B coordinating with N and C (NiN₄ B₈ -C) can effectively optimize the Gibbs free energy of oxygen intermediate states and hence improve the catalytic activity of the materials. In addition, we propose that the trend change in catalytic activity is mainly governed by the filling of antibonding orbitals between Ni-3d and O-2p states near the Fermi level.

Assessment of the Accuracy of Density Functionals for Calculating Oxygen Reduction Reaction on Nitrogen-Doped Graphene

Experimental studies of the oxygen reduction reaction (ORR) at nitrogen-doped graphene electrodes have reported a remarkably low overpotential, on the order of 0.5 V, similar to Pt-based electrodes. Theoretical calculations using density functional theory have lent support to this claim. However, other measurements have indicated that transition metal impurities are actually responsible for the ORR activity, thereby raising questions about the reliability of both the experiments and the calculations. To assess the accuracy of the theoretical calculations, various generalized gradient approximation (GGA), meta-GGA, and hybrid functionals are employed here and calibrated against high-level wave-function-based coupled-cluster calculations (CCSD(T)) of the overpotential as well as self-interaction corrected density functional calculations and published quantum Monte Carlo calculations of O adatom binding to graphene. The PBE0 and HSE06 hybrid functionals are found to give more accurate results than the GGA and meta-GGA functionals, as would be expected, and for a low dopant concentration, 3.1%, the overpotential is calculated to be 1.0 V. The GGA and meta-GGA functionals give a lower estimate by as much as 0.4 V. When the dopant concentration is doubled, the overpotential calculated with hybrid functionals decreases, while it increases in GGA functional calculations. The opposite trends result from different potential-determining steps, the *OOH species being of central importance in the hybrid functional calculations, while the reduction of *O determines the overpotential obtained in GGA and meta-GGA calculations. The results presented here are mainly based on calculations of periodic representations of the system, but a comparison is also made with molecular flake models that are found to give erratic results due to finite size effects and geometric distortions during energy minimization. The presence of the electrolyte has not been taken into account explicitly in the calculations presented here but is estimated to be important for definitive calculations of the overpotential.

Catalytic CO oxidation reaction over N-substituted graphene nanoribbon with edge defects

Density functional theory calculations, including dispersion effects, are used to demonstrate how substitutional nitrogen atoms can improve the catalytic reactivity of graphene nanoribbons (GNR) with edge defects in the CO oxidation process. It is demonstrated that the addition of nitrogen impurities significantly enhances O₂ adsorption on GNR. Carbon atoms near the edges of defects are the most active sites for capturing O₂ molecules. The lower adsorption energy of CO relative to O₂ implies that the N-modified GNR is resistant to CO poisoning. The Eley-Rideal (E-R) mechanism has activation energies as low as 0.38 eV, making it the most energetically relevant pathway for the CO + O₂ reaction. The findings of this study might help in the design of catalysts for metal-free catalysis of CO oxidation.

Synergic effects between boron and nitrogen atoms in BN-codoped C59-n BN n fullerenes (n = 1-3) for metal-free reduction of greenhouse N2O gas

The geometries, electronic structures, and catalytic properties of BN-codoped fullerenes C59-n BN n (n = 1-3) are studied using first-principles computations. The results showed that BN-codoping can significantly modify the properties of C60 fullerene by breaking local charge neutrality and creating active sites. The codoping of B and N enhances the formation energy of fullerenes, indicating that the synergistic effects of these atoms helps to stabilize the C59-n BN n structures. The stepwise addition of N atoms around the B atom improves catalytic activities of C59-n BN n in N2O reduction. The reduction of N2O over C58BN and C57BN2 begins with its chemisorption on the B-C bond of the fullerene, followed by the concerted interaction of CO with N2O and the release of N2. The resulting OCO intermediate is subsequently transformed into a CO2 molecule, which is weakly adsorbed on the B atom of the fullerene. On the contrary, nitrogen-rich C56BN3 fullerene is found to decompose N2O into N2 and O* species without the requirement for activation energy. The CO molecule then removes the O* species with a low activation barrier. The activation barrier of the N2O reduction on C56BN3 fullerene is just 0.28 eV, which is lower than that of noble metals.

Design of Local Atomic Environments in Single-Atom Electrocatalysts for Renewable Energy Conversions

Single-atom electrocatalysts (SAECs) have recently attracted tremendous research interest due to their often remarkable catalytic responses, unmatched by conventional catalysts. The electrocatalytic performance of SAECs is closely related to the specific metal species and their local atomic environments, including their coordination number, the determined structure of the coordination sites, and the chemical identity of nearest and second nearest neighboring atoms. The wide range of distinct chemical bonding configurations of a single-metal atom with its surrounding host atoms creates virtually limitless opportunities for the rational design and synthesis of SAECs with tunable local atomic environment for high-performance electrocatalysis. In this review, the authors first identify fundamental hurdles in electrochemical conversions and highlight the relevance of SAECs. They then critically examine the role of the local atomic structures, encompassing the first and second coordination spheres of the isolated metal atoms, on the design of high-performance SAECs. The relevance of single-atom dopants for host activation is also discussed. Insights into the correlation between local structures of SAECs and their catalytic response are analyzed and discussed. Finally, the authors summarize major challenges to be addressed in the field of SAECs and provide some perspectives in the rational construction of superior SAECs for a wide range of electrochemical conversions.

Catalytic role of graphitic nitrogen atoms in the CO oxidation reaction over N-containing graphene: a first-principles mechanistic evaluation

The catalytic role of graphitic nitrogen atoms of a series of nitrogen-doped graphene surfaces is explored for low-temperature oxidation of CO using periodic DFT calculations.

Novel joint catalytic properties of Fe and N co-doped graphene for CO oxidation

Using density functional theory, we have performed detailed calculations of the joint catalytic activity of graphene co-doped with Fe and N atoms. The Fe atom can be located on single vacancy graphene and acts as the active site. Due to the strong attraction of the Fe ion, the O-O bond length of the O2 molecule is elongated, which decreases the bonding energy between the O atoms. The energy barrier of CO oxidization is 0.84 eV. When N atoms are doped into the graphene, the interactions between the Fe ions and O2 molecules are stronger, and the O-O bond lengths are elongated further, which makes the desorption of the quasi-CO2 molecule easier. The energy barriers are reduced to 0.62 eV, 0.49 eV, and 0.33 eV for graphene doped with one, two and three N atoms, respectively. The O atom remaining on the Fe ion can form a CO2 molecule with an additional CO molecule. The produced CO2 molecule can be released with a small or even zero energy barrier by adsorbing an O2 molecule. The adsorbed O2 molecule is then involved in the next reaction process, and the material can be used as a recycled catalyst.

First-principles description of electrocatalytic characteristics of graphene-like materials

Graphene-like materials (GLMs) have received much attention as a potential alternative to precious metal-based electrocatalysts. However, the description of their electrocatalytic characteristics may still need to be improved, especially under constant chemical potential. Unlike the case of conventional metal electrodes, the potential drop across the electrical double layer (ϕ_D) at the electrode-electrolyte interface can deviate substantially from the applied voltage (ϕ_app) due to a shift of the Dirac point (eϕ_G) with charging. This may in turn significantly alter the interfacial capacitance (C_T) and the relationship between ϕ_app and free-energy change (ΔF). Hence, accurate evaluation of the electrode contribution is necessary to better understand and optimize the electrocatalytic properties of GLMs. In this work, we revisit and compare first-principles methods available to describe the ϕ_app-∆F relation. Grand-canonical density functional theory is used to determine ΔF as a function of ϕ_app or electrode potential (ϕ_q), from which the relative contribution of eϕ_G is estimated. In parallel, eϕ_G is directly extracted from a density functional theory analysis of the electronic structure of uncharged GLMs. The results of both methods are found to be in close agreement for pristine graphene, but their predictions deviate noticeably in the presence of adsorbates; the origin of the discrepancy is analyzed and explained. We then evaluate the application of the first-principle methods to prediction of the electrocatalytic processes, taking the reduction (hydrogenation) and oxidation (hydroxylation) reactions on pristine graphene as examples. Our work highlights the vital role of the modification of the electrode electronic structure in determining the electrocatalytic performance of GLMs.

Theoretical Study on a Nitrogen-Doped Graphene Nanoribbon with Edge Defects as the Electrocatalyst for Oxygen Reduction Reaction

Both theory and experiment show that sp2 carbon nanomaterials doped with N have great potential as high-efficiency catalysts for oxygen reduction reactions (ORR). At present, there are theoretical studies that believe that C-sites with positive charge or high-spin density values have higher adsorption capacity, but there are always some counter examples, such as the N-doped graphene nanoribbons with edge defects (ND-GNR) of this paper. In this study, the ORR mechanism of ND-GNR was studied by density functional theory (DFT) calculation, and then the carbon ring resonance energy was analyzed from the perspective of chemical graph theory to elucidate the cause and distribution of active sites in ND-GNR. Finally, it was found that the overpotential of the model can be adjusted by changing the width of the model or dopant atoms while still ensuring proper adsorption energy (between 0.5 and 2.0 eV). The minimum overpotential for these models is approximately 0.36 V. These findings could serve as guidelines for the construction of efficient ORR carbon nanomaterial catalysts.

N-Doped Graphene Supported on Metal-Iron Carbide as a Catalyst for the Oxygen Reduction Reaction: Density Functional Theory Study

The development of an efficient electrocatalyst for the oxygen reduction reaction (ORR) is essential for the commercialization of fuel-cell technologies. Iron carbide encapsulated in N-doped graphene (NG/Fe₃ C) has been recognized recently as a promising ORR catalyst. In this study, the stability and catalytic activity of N-doped graphene supported on metal-iron carbide (NG/M_Fe₃ C) toward the ORR are investigated by using DFT calculations. The NG/M_Fe₃ C heterostructure is modeled by substituting Fe atoms in the Fe₃ C substrate near the NG/Fe₃ C interface by metal atoms M (M=Cr-Mn, Co-Zn, Nb-Mo, Ta-W). The calculations show that the introduction of the metal atoms M alters the work function of the overlayer N-doped graphene, which is found to correlate with the binding strength of the ORR intermediates. The introduction of Ni or Co atoms at the interface improves the ORR activity of the NG/Fe₃ C and stabilizes the heterostructure. The ORR activity increases as the concentration of Ni or Co atoms near the interface increases, and the stable heterostructure is available in a wide range of substituted concentrations. These results suggest approaches to improve the ORR activity of NG/Fe₃ C catalysts.

Catalytic trends of nitrogen doped carbon nanotubes for oxygen reduction reaction

Replacing the state-of-the-art fuel cell catalyst platinum for a cheaper and abundant alternative would make the hydrogen economy viable. Both nitrogen-doped graphene and nitrogen-doped carbon nanotubes (N-CNT) have been shown to be capable of acting as a metal-free catalyst for the oxygen reduction reaction (ORR). Until now, most of the research has been focused on the nitrogen doping and less on the structure of the nanotubes. Here, density functional theory calculations are used to calculate trends in ORR catalytic activity of graphitic-N-doped CNTs of different sizes and chirality of selected tubes between (4,0) and (20,10). This includes 13 armchair tubes, 17 zig-zag tubes and 42 chiral tubes, or 72 N-CNTs in total. 22 tubes are predicted to have a lower overpotential than the platinum catalyst and 46 tubes have lower overpotential than nitrogen doped graphene. The most active tubes are (14,7), (12,6), and (8,8), and display an overpotential of around 0.35 V, or 0.1 V lower overpotential than predicted on Pt(111) with the same level of theory.

Catalytic performance of graphene quantum dot supported manganese phthalocyanine for efficient oxygen reduction: density functional theory approach

The thermodynamic free-energy diagrams predict that MnPc/GQD is more active toward ORR than the isolated MnPc, clearly highlighting the effect of the GQD matrix on ORR activity from a thermodynamic perspective.

Oxygen reduction reaction mechanism on a phosporus-doped pyrolyzed graphitic Fe/N/C catalyst

The interaction between neighbouring FeN₄ and P_subs sites of a P-doped pyrolyzed Fe/N/C catalyst promotes four-electron reduction through associative and dissociative mechanisms.

The chemical identity, state and structure of catalytically active centers during the electrochemical CO2 reduction on porous Fe-nitrogen-carbon (Fe-N-C) materials

We report novel structure–activity relationships and explore the chemical state and structure of catalytically active sites under operando conditions during the electrochemical CO₂ reduction reaction (CO₂RR) catalyzed by a series of porous iron–nitrogen–carbon (FeNC) catalysts.

We report novel structure–activity relationships and explore the chemical state and structure of catalytically active sites under operando conditions during the electrochemical CO₂ reduction reaction (CO₂RR) catalyzed by a series of porous iron–nitrogen–carbon (FeNC) catalysts. The FeNC catalysts were synthesized from different nitrogen precursors and, as a result of this, exhibited quite distinct physical properties, such as BET surface areas and distinct chemical N-functionalities in varying ratios. The chemical diversity of the FeNC catalysts was harnessed to set up correlations between the catalytic CO₂RR activity and their chemical nitrogen-functionalities, which provided a deeper understanding between catalyst chemistry and function. XPS measurements revealed a dominant role of porphyrin-like Fe–N_x motifs and pyridinic nitrogen species in catalyzing the overall reaction process. Operando EXAFS measurements revealed an unexpected change in the Fe oxidation state and associated coordination from Fe²⁺ to Fe¹⁺. This redox change coincides with the onset of catalytic CH₄ production around –0.9 V_RHE. The ability of the solid state coordinative Fe¹⁺–N_x moiety to form hydrocarbons from CO₂ is remarkable, as it represents the solid-state analogue to molecular Fe¹⁺ coordination compounds with the same catalytic capability under homogeneous catalytic environments. This finding highlights a conceptual bridge between heterogeneous and homogenous catalysis and contributes significantly to our fundamental understanding of the FeNC catalyst function in the CO₂RR.

Is iron nitride or carbide highly active for oxygen reduction reaction in acidic medium?

The high ORR activities in acidic medium for the so-called Fe₂N- and Fe₃C-based catalysts actually originate from the trace of Fe–N_x/C moieties (x ≥ 4) rather than the Fe₂N or Fe₃C phases themselves.

Theoretical Studies of Oxygen Reactivity of Free-Standing and Supported Boron-Doped Graphene

Graphene inertness towards chemical reactivity can be considered as an accepted postulate by the research community. This limit has been recently overcome by chemically and physically modifying graphene through non-metal doping or interfacing with acceptor/donor materials (metals or semiconductors). As a result, outstanding performances as catalytic, electrocatalytic, and photocatalytic material have been observed. In this critical Review we report computational work performed, by our group, on the reactivity of free-standing, metal- and semiconductor-supported B-doped graphene towards oxygen, which is at the basis of extremely important energy-related chemical processes, such as the oxygen reduction reaction. It appears that a combination of doping and interfacing approaches for the activation of graphene can open unconventional and unprecedented reaction paths, thus boosting the potential of modified graphene in many chemical applications.

DFT calculation of oxygen adsorption on a core-single shell ZnNb catalyst

Oxygen adsorption onto Pt, Nb, and core-single shell ZnNb particles is studied using ab initio molecular orbital (MO) calculations to realize their ability for the oxygen reduction reaction (ORR).

Experimental Observation of Redox-Induced Fe-N Switching Behavior as a Determinant Role for Oxygen Reduction Activity

The commercialization of electrochemical energy conversion and storage devices relies largely upon the development of highly active catalysts based on abundant and inexpensive materials. Despite recent achievements in this respect, further progress is hindered by the poor understanding of the nature of active sites and reaction mechanisms. Herein, by characterizing representative iron-based catalysts under reactive conditions, we identify three Fe-N4-like catalytic centers with distinctly different Fe-N switching behaviors (Fe moving toward or away from the N4-plane) during the oxygen reduction reaction (ORR), and show that their ORR activities are essentially governed by the dynamic structure associated with the Fe(2+/3+) redox transition, rather than the static structure of the bare sites. Our findings reveal the structural origin of the enhanced catalytic activity of pyrolyzed Fe-based catalysts compared to nonpyrolyzed Fe-macrocycle compounds. More generally, the fundamental insights into the dynamic nature of transition-metal compounds during electron-transfer reactions will potentially guide rational design of these materials for broad applications.

A systematic study of metal-supported boron nitride materials for the oxygen reduction reaction

Surfaces that efficiently catalyse the oxygen reduction reaction (ORR) are highly desirable for applications in energy utilization. Here, we computationally investigate the ORR on hexagonal boron nitride (h-BN) supported on Ni, Cu, and Co. We find a significant influence of the metal on the reaction energetics. In particular, h-BN/Cu is predicted to catalyse the ORR with a low overpotential, while on the other substrates the reaction is impeded by the formation of too stable surface hydroxyl species. Our results highlight trends in the reactivity of these heterostructures and may guide further rational design of O2-activating catalysts based on supported h-BN.

Activation of Graphenic Carbon Due to Substitutional Doping by Nitrogen: Mechanistic Understanding from First Principles

Nitrogen-doped graphene and carbon nanotubes are popularly in focus as metal-free electrocatalysts for oxygen reduction reactions (ORR) central to fuel cells. N-doped CNTs have been also reported to chemisorb mutually, promising a route to their robust predetermined assembly into devices and mechanical reinforcements. We propose from first principles a common mechanistic understanding of these two aspects pointing further to a generic chemical activation of carbon atoms due to substitution by nitrogen in experimentally observed configurations. Wannier-function based orbital resolved study of mechanisms suggests increase in C-N bond-orders in attempt to retain π-conjugation among carbon atoms, causing mechanical stress and loss of charge neutrality of nitrogen and carbon atoms, which remedially facilitate chemical activation of N-coordinated C atoms, enhancing sharply with increasing coordination to N and proximity to zigzag edges. Activated C atoms facilitate covalent adsorption of radicals in general, diradicals like O2 relevant to ORR, and also other similarly activated C atoms, leading to self-assembly of graphenic nanostructures while remaining inert to ordinary graphenic C atoms.

Activity of N-coordinated multi-metal-atom active site structures for Pt-free oxygen reduction reaction catalysis: role of *OH ligands

We report calculated oxygen reduction reaction energy pathways on multi-metal-atom structures that have previously been shown to be thermodynamically favorable. We predict that such sites have the ability to spontaneously cleave the O2 bond and then will proceed to over-bind reaction intermediates. In particular, the *OH bound state has lower energy than the final 2 H2O state at positive potentials. Contrary to traditional surface catalysts, this *OH binding does not poison the multi-metal-atom site but acts as a modifying ligand that will spontaneously form in aqueous environments leading to new active sites that have higher catalytic activities. These *OH bound structures have the highest calculated activity to date.

Oxygen reduction reaction on neighboring Fe–N4 and quaternary-N sites of pyrolized Fe/N/C catalyst

Interactions of Fe–N₄ and quaternary-N sites at the zigzag edge of graphene enhance the ORR free-energy profile of the quaternary-N site.

Nitrogen-doped graphene nanoribbons as efficient metal-free electrocatalysts for oxygen reduction

Nitrogen-doped graphene nanoribbon (N-GNR) nanomaterials with different nitrogen contents have been facilely prepared via high temperature pyrolysis of graphene nanoribbons (GNR)/polyaniline (PANI) composites. Here, the GNRs with excellent surface integration were prepared by longitudinally unzipping the multiwalled carbon nanotubes. With a high length-to-width ratio, the GNR sheets are prone to form a conductive network by connecting end-to-end to facilitate the transfer of electrons. Different amounts of PANI acting as a N source were deposited on the surface of GNRs via a layer-by-layer approach, resulting in the formation of N-GNR nanomaterials with different N contents after being pyrolyzed. Electrochemical characterizations reveal that the obtained N8.3-GNR nanomaterial has excellent catalytic activity toward an oxygen reduction reaction (ORR) in an alkaline electrolyte, including large kinetic-limiting current density and long-term stability as well as a desirable four-electron pathway for the formation of water. These superior properties make the N-GNR nanomaterials a promising kind of cathode catalyst for alkaline fuel cell applications.

A density functional theory study on oxygen reduction reaction on nitrogen-doped graphene

Nitrogen (N)-doped carbons reportedly exhibit good electrocatalytic activity for the oxygen reduction reaction (ORR) of fuel cells. This work provides theoretical insights into the ORR mechanism of N-doped graphene by using density functional theory calculations. All possible reaction pathways were investigated, and the transition state of each elementary step was identified. The results showed that OOH reduction was easier than O-OH breaking. OOH reduction followed a direct Eley-Rideal mechanism (the OOH species was in gas phase, but H was chemisorbed on the surface) with a significantly low reaction barrier of 0.09 eV. Pathways for both four-electron and two-electron reductions were possible. The rate-determining step of the two-electron pathway was the reduction of O₂ (formation of OOH), whereas that of the four-electron pathway was the reduction of OH into H₂O. After comparing the barriers of the rate-determining steps of the two pathways, we found that the two-electron pathway was more energetically favored than the four-electron pathway.