Karim
Khan†
*abcde,
Ayesha Khan
Tareen†
e,
Muhammad
Aslam†
e,
Yupeng
Zhang
cd,
Renheng
Wang
cd,
Zhengbiao
Ouyang
*bc,
Zhongyi
Gou
*a and
Han
Zhang
*cd
aAdvanced electromagnetic function laboratory, Dongguan University of Technology (DGUT), Dongguan, Guangdong Province, P.R. China. E-mail: karim_khan_niazi@yahoo.com; guozhongyi@hfut.edu.cn
bCollege of Electronic Science and Technology, and THz Technical Research Center of Shenzhen University, Shenzhen, 518060, P.R. China
cKey Laboratory of Optoelectronics Devices and Systems of Ministry of Education and Guangdong Province Shenzhen University, Shenzhen, 518060, P.R. China. E-mail: wangrh@szu.edu.cn; zbouyang@szu.edu.cn; hzhang@szu.edu.cn
dShenzhen Engineering Laboratory of Phosphorene and Optoelectronics, and SZU-NUS Collaborative Innovation Center for Optoelectronic Science and Technology, Shenzhen University, Shenzhen, 518060, P.R. China
eGovernment Degree college PaharPur, Gomel University, Dera Ismail Khan, K.P.K., Islamic Republic of Pakistan. E-mail: chemistayesha@yahoo.com
First published on 8th November 2019
Two-dimensional (2D) materials have a wide platform in research and expanding nano- and atomic-level applications. This study is motivated by the well-established 2D catalysts, which demonstrate high efficiency, selectivity and sustainability exceeding that of classical noble metal catalysts for the oxygen reduction reaction (ORR), oxygen evolution reaction (OER), and/or hydrogen evolution reaction (HER). Nowadays, the hydrogen evolution reaction (HER) in water electrolysis is crucial for the cost-efficient production of a pure hydrogen fuel. We will also discuss another important point related to electrochemical carbon dioxide and nitrogen reduction (ECR and N2RR) in detail. In this review, we mainly focused on the recent progress in the fuel cell technology based on 2D materials, including graphene, transition metal dichalcogenides, black phosphorus, MXenes, metal–organic frameworks, and metal oxide nanosheets. First, the basic attributes of the 2D materials were described, and their fuel cell mechanisms were also summarized. Finally, some effective methods for enhancing the performance of the fuel cells based on 2D materials were also discussed, and the opportunities and challenges of 2D material-based fuel cells at the commercial level were also provided. This review can provide new avenues for 2D materials with properties suitable for fuel cell technology development and related fields.
In fuel cells, fuels are oxidized on the anode, and released electrons transfer through an external circuit to the cathode where O2 is reduced. There are two major energy producing cycles, which rely on carbon and water (H2O) and have been extensively researched. The hydrolysis of H2O is caused hydrogen and oxygen evolution at the cathode (HER) and anode (OER), respectively, in the electrolytic cell. In the opposite process, the H2 (acts as fuel) produced at the anode is oxidized (HOR), while oxygen at the cathode is reduced (ORR) in H2/air fuel cell to produce electricity, with H2O as the major product. The HER and HOR proceed at kinetically high rates; whereas the ORR is slow because it entails multiple electron–proton coupled transfer steps. The combustion of fossil fuels releases CO2, which when emitted by carbonaceous materials can be recycled through the carbon cycle for utilization in the production of fuels (e.g., methanol, methane, and ethanol). Although well-known catalysts such as nanoscale Pt and costly metal oxides (e.g., RuO2 and IrO2) show high catalytic ability, their scarcity and high cost restrict their industrial applications. Thus, new NPMCs with confined dimensions (2D) have been introduced, but their atomic level understanding is still an ongoing debate.1,6–9
To fulfill future electrochemical energy demands, it is crucial to explore advanced electrocatalysts with well-controlled nanostructures for dramatic enhancements in energy and power density, catalytic activity, efficiency, and durability.1,6–10 One renewable energy generation/distribution route called the hydrogen economy concept is receiving significant attention from researchers. The hydrogen economy concept utilizes H2 gas as the main energy source for the efficient running of buildings, homes and vehicles. Unfortunately, due to the many inefficiencies associated with the energy conversion device (electrolyzer) needed to produce pure H2, compatible fuel cell devices to convert H2 into electricity and store the generated electricity are necessary (batteries and supercapacitors (SC)). Furthermore, the development of cheap and active catalysts that rival the state-of-the-art materials for these devices will make the hydrogen economy concept closer to reality. Therefore, the need for alternative, abundant, and financially feasible electrode materials is urgent. 2D materials compared to their bulk counterparts show improved activity as electrocatalysts.11 The energy crisis and global challenges provide excellent opportunities for designing and exploring new electrocatalysts.1,6–10 Electrocatalysis is influenced by the electrode materials, kinetics, and adsorption/desorption at the electrode/electrolyte interface, and it assists in the conversion of chemical and electrical energy in energy conversion and storage devices. Three seemingly simple reactions, the ORR, OER and HER, are significant for green and sustainable energy-producing devices, e.g., fuel cells, water-electrolysis devices, and batteries.12 Accordingly, cost-efficient, alternative, 2D nanosheets materials with large active sites are beneficial for prompt interfacial charge transfer with feasible catalysis. Among the 2D materials, graphene and its derivatives, carbides, nitrides, chalcogenides, and non-noble metal oxides introduce a new generation of electrocatalysts. These 2D materials can act as cathode materials in electrochemical reactions, and reduce the largest portion of electrochemical kinetic loss and system costs. Hence, this has promoted studies on 2D material-based electrocatalysts, and, to some extent, given rise to the recent research boom on the OER and HER for the CO2/N2 electrochemical reduction reactions (ECR/ENR). Most publications focused on discussing different approaches and strategies for the preparation of a specific class of ORR catalysts (e.g., Pt-based core–shell catalysts, nano-structured Pt-alloy electrocatalysts and doped-carbon materials), but few reviews are available regarding the electrocatalytic applications of 2D materials, only one or a few 2D families. In view of the rapid development of 2D electrocatalysts, there is a great need to provide timely updates in this field with emphasis on new design concepts and significant breakthroughs in electrocatalysts. Strategic applications of different families of 2D nanomaterials in electrocatalysis should be analyzed, focusing on the role of edges and surface characteristics. Since the last decade, attention on nanostructured materials (NSMs) has demonstrated that new characteristics are gained at the nanoscale, which change with a change in size or morphology because the size of NSMs determines novel and exciting aspects.13,14 In the middle of the 20th century, for the first time, Gleiter15 gave the classification of NSMs, which was promoted by Skorokhod.16 At that time, this classification was not complete since low-dimensional structures, i.e. fullerenes, nanoflakes and nanotubes, were not included because their effects determined by size are part of the evolution of the physicochemical features of these restricted systems with a change in size.17–20 With the development of new technologies, a new class of low-dimensional materials has resulted in the growth of nanoscale sciences, which are known as 2D materials. The fundamental material parameter is dimensionality, which decides their properties to a significant extent. 2D materials or compounds with the same chemical composition can exhibit different properties.21
Different dimension nanomaterials have different density of states (DOS) due to the number of different electrons present. For example, zero-dimensional (0D) quantum dots,23,24 one-dimensional (1D) nanoribbons, nanotubes, and nanowires, two-dimensional (2D) single-atom-thick materials, and three-dimensional (3D) nano-balls and nanocones (Fig. 1).21,25–42 For 0D NSMs, excited electrons are completely confined in all directions and they are equivalent to atomic/molecular clusters with discreet states, which are well-separated in dissimilar energy states. Michael Faraday in the middle of 19th century presented the first experiments with the nanomaterial colloidal gold, and determined its size-dependent optical properties.43 The 0D NSMs are fabricated via physical and chemical methods due to their well-ordered dimensions. Recently, 0D NSM heterogeneous particle arrays including uniform quantum dots, nanoplatelets, nano-onions, hollow spheres, and core–shell quantum dots have been developed with some potential application from the laboratory to the industrial scale.1,6–10,44–46 Furthermore, regarding the application of 0D materials, i.e. 0D quantum dots, they have been widely applied in solar cells,47 light emitting diodes (LEDs),48 lasers,49,50 and single-electron transistors.51 They are also used in routine applications such as pregnancy test-kits (gold NSMs), antibacterial bandages (silver NSMs), sun-screen lotions formulations (SiO2 NSMs or other oxide ceramic nanoparticles) as few examples of consumer products that contain 0D NSMs.
Fig. 1 Dimensionality spectrum of sp2 carbon materials. (a) Structure model of a fullerene, where the red circle shows it trapped inside the cavity of an MWCNT. (b) Structure model of CNT and STM image of an SWNT exposed at the surface of a rope. (c) Structure model of graphene and a snapshot (ordinary white light) of a comparatively large graphene multilayer flake of ∼3 nm thickness on an Si wafer substrate. (d) Structure model of graphite and a piece of natural graphite. This figure has been adapted/reproduced from ref. 22 with permission from the Royal Society of Chemistry. |
Similarly, the second category of NSMs is 1D nanomaterials, which have nanoscale sizes along two dimensions with a rod- or wire-like appearance. In these nanomaterials, quantum confinement and SSA-related nanoscale effects, which are absent in 0D nanomaterials, are pronounced, allowing them to be integrated into device architectures due to their single bulk-like dimensions.21 1D nanostructured materials, such as nanowires and nanotubes, have attracted substantial attention after groundbreaking studies, especially the discovery of the good transparency with prominent conductivity for transparent and flexible devices in 2D sheet materials. In this review, we mostly focus on the applications of 2D materials; therefore, from here we will discuss in detail the different kinds of 2D materials and their potential industrial applications, especially the establishment of metrics that can facilitate definite targets to be to achieved for technology challenges in energy production and storage devices.
The potential of future advanced technologies ensures that materials can be manipulated at the nanometer level and outcomes from new tools and technologies not visualized previously, ranging from nanoparticle-based cancer therapy61,62 to efficient composites, and thus, nanotechnology launches new boundaries for ideas in medicine, electronic devices, etc.23,33,38,62–78 Advanced research and development activity on 2D materials have already shown that these materials are appropriate for large-scale manufacturing. The fabrication of 2D materials is one conspicuous example of the fast progress toward their commercialization, with the random discovery of micro- to large-scale flakes in the laboratory and roll-to-roll manufacturing of nanosheets with sizes approaching the meter-scale. Additionally, the mechanical strength of 2D materials is very high, 200 times higher than that of steel, which gives these materials extraordinary breaking strengths and unpredictable flexible. Their mechanical properties can be further tuned with defects or functional groups, which makes 2D materials strong as lightweight polymer. Furthermore, as ultrathin stretchable membranes, single layer materials are suitable for use in tunable electrochemical systems. They can withstand more than 20% elastic distortion and also present extraordinarily high flexibility.79 At larger strain, these materials are considered to be brittle since they will break like glass because they are very thin. Chemical modification via chemical doping with various metal ions may impart more elastic distortion strength. Besides elasticity in atomic-level-thin layers, they also have confined quantum levels, and hence some remarkable electronic and optical properties.80
Graphene (2004) is a single-atom layer derived from graphite, where carbon atom are connected in a hexagonal honeycomb lattice, and it is a very prominent and mature 2D-layered material, which has opened a new world of 2D materials and owing to its exceptional properties not necessarily found in the familiar 3D world.25,28,30,75,81,82 Besides graphene, its other carbon-based 2D material siblings/family include, graphene oxide (GO), reduced graphene oxide (rGO), graphane, fluorographene, graphyne, and graphdiyne. Their thicknesses range from one to several atoms and their hallmark feature is their unique bonding interactions. These ultrathin sheets typically demonstrate properties that are dissimilar to the parent bulk material due to anisotropy. Notably, the novel structure of graphene endows it with various superior characteristics for example elevated electrical and thermal conductivities, good transparency, high mechanical strength, intrinsic flexibility and large specific surface area (SSA). The rapid development of graphene covers chemistry, physics, biology, materials science and related interdisciplinary fields; however, the large-scale synthesis of big-sized single-layered or few-layered 2D carbon nanosheets is requires for industrial revolution. Graphene is the first and best known 2D material1,8 with a zero band gap semimetal nature, where electrons are highly mobile and exhibit exceptional conductivity.74,83 Inspired by the attractive properties of graphene due to its 2D structure, researchers have been stimulated to hunt for new 2D material families. The 2D class of materials consists of a wide variety of materials with varying chemical compositions and tunable properties.
The unique characteristics of 2D structures such as high surface-volume ratio, shape, surface charge, anisotropic nature and tunable functionalities open up their potential application scope further. Generally 2D materials are categorized by their structure, including graphene, transition metal dichalcogenides (TMDCs), layered double hydroxides (LDHs), gC3N4, LAPONITE® clay, hexagonal boron nitride (h-BN), and black phosphorous (BP). 2D materials consistently have exposed crystal lattices compared to other materials, for example, nanowires, nanoporous materials, and nanotubes, and thus have superior exposed electrocatalytically active sites at a specific loaded catalyst concentration. The ordered and simple molecular framework of 2D materials also endow them with an active SSA according to both theoretical and experimental approaches. An additional reason to investigate and expand 2D electrocatalysts is they are usually cost effective, which make them novel substitutes for precious metal-based catalysts (PMCs). Fig. 2 shows the classification and structures of 2D nanomaterials. Therefore, the rapidly increasing attention on graphene is the last few years has led to the investigation of other 2D atomic crystals, including single- or few-layer TMDs (e.g. molybdenum disulfide (MoS2), tungsten disulfide (WS2), molybdenum diselenide (MoSe2), and tungsten diselenide (WSe2)),85 h-BN, borophene (2D boron), silicene (2D silicon), germanene (2D germanium), metal oxides, LDHs, graphitic carbon nitride (g-C3N4), a family of mono-elemental compounds, BP (or phosphorene), arsenene, antimonene and bismuthene (Xenes), and metal nitrides and carbides (collectively known as MXenes).6,27,64,86 2D nanomaterials have expanded to semiconductor III–VI compounds (InSe, GaSe, etc.),87 transition metal halides (TMHs) (e.g., PbI2 and MgBr2),88 metal oxides (e.g., MnO2 and MoO3),89 perovskite-type oxides (e.g., K2Ln2Ti3O10 and RbLnTa2O7 (Ln: lanthanide ion)),90 hexagonal boron nitride (h-BN or white graphene),91 graphitic carbon nitride (g-C3N4),92 2D polymers and frameworks,93 and MXenes, which present versatile opportunities.94 A TMDC single layer with MX2 stoichiometry consists of a hexagonally stacked transition metal (TM) (e.g. M = Ti, W, Nb, Mo, and Ta) inserted among two planes of chalcogen atoms (X = S, Se and Te).95 In contrast to the semimetal graphene with a zero band gap, 2D TMDs may be semiconductors (e.g., MoS2, MoSe, WS2, and WSe2), semimetals (e.g., TiSe2 and WTe2), or metallic conductors (e.g., VSe2 and NbS2), relying on the coordination setting of the individual metal atom, the particular metal d-electron configuration, or doping the relevant crystal structure.96 In this review, we will discuss in detail the recent progress made to date in the application of 2D materials as electrocatalysts for the ORR, OER, HER, and electrochemical CO2/N2 reduction (ECR/ENR), and also the limitations of traditional catalysts. In this section, we present each family in detail and their corresponding electrocatalytic applications based on important properties.
Fig. 2 Structures of 2D nanomaterials, including carbon-based nanomaterials (graphene, graphene oxide (GO and rGO), silicate clays, LDHs, TMDCs and TMOs. This figure has been adapted/reproduced from ref. 84 with permission from Wiley. |
Fig. 3 2D materials ranging from insulators and metals to semiconductors with diverse bandgaps, and electron/hole mobility summary scale diagram. This figure has been adapted/reproduced from ref. 104 with permission from Nature Publishing Group. |
Dai et al. in 2009 reported N-doped carbon nanotubes (N-CNT), and showed that they can replace Pt for the ORR in fuel cells. After this discovery, significant research progress was made to perfect the potential properties of 2D materials as metal-free electrocatalysts in various renewable energy formation processes. Among the these materials, the graphene family, graphdiyne, graphitic carbon nitride (g-C3N4), h-BN, and BP, will be discussed in detail in the following section. To date, among the 2D-based materials, graphene and graphene-based nanocomposite materials are promising electrode materials with high performances in the field of energy devices (sensors, electronics and optoelectronics, efficient catalysts, etc.) due to their high electrical conductivity and large SSA. The enormous achievement of graphene has been followed by remarkable search for further 2D materials at the atomic level for amazing feature/characteristics.1,6–9 The application of graphene in photovoltaic devices, supercapacitors (SCs), batteries, and fuel cells for hydrogen generation is an opportunity to deal with the increasing challenges globally in the energy field. The 2D graphene with a simulated ∼2600 m2 g−1 surface to mass ratio and high flexibility and electrical conductance show it can store electric ions, charge, or hydrogen. The benefit of using 2D crystals is the opportunity to generate and control layered synthetic structures with “customized” characteristics through their spins (spintronic). This results in the development of accessible materials with insulating, semiconducting, and metallic properties.
Fig. 4 Electrocatalytic applications of 2D materials in energy production (fuel cells) and energy saving devices (batteries and superconductors). This figure has been adapted/reproduced from ref. 129 with permission from Elsevier Inc. |
Fuel cell-based energy production involves multi-electron charge transfer processes, for instance, the ORR is classified as a 2e− or 4e− system. These electrons though external circuit are channeled via the anode to cathode, which produces a current. Along with produced electricity, water is a by-product from its constituent atoms, i.e. H2 and O2. In the opposite case, H2 (fuel) and O2 are generated via the application of comparatively low energy. Thus, the combination of forward and reverse reactions, consisting of consecutive hydrogen- and oxygen-associated catalytic processes, makes a cycle, which is called the water cycle. Progress in this technology originated from the increasing hydrogen market, which has caused the electrocatalytic water cycle to be an important research topic. The perfect technical water cycle is composed of water electrolysis via the OER and HER for fuel production, followed by power production via the ORR and HOR. For electrolytic cells, the HER and OER occur at the cathode and anode, creating H2 gas and O2 molecules, respectively (Fig. 5(a)). The ORR and HOR occur at the cathode and anode for H2 and O2, respectively (Fig. 5(b)). In the last two decades, considerable investigations have addressed the progress in achieving high electrocatalytic properties for both the water and carbon cycles. Generally, the kinetics of the 2e− transfer method in half-cell HER and HOR is straightforward. However, the slow rates for the multi-electron transfer ORR and OER decrease the efficiency of the energy equipment. Conventionally, noble metals (Pt, Ir, and Ru), are used as a electrocatalysts for these reactions, but these expensive electrocatalysts suffer from poor stability, which significantly limits their large-scale, long-term industrial application. Thus, the development of low cost, high performance electrocatalysts has become important for clean and sustainable energy production.128 The water cycle was a novel achievement summed up in Grove's letter, On a gaseous voltaic battery, which was published in 1842.136 His working principle relied on the explanation by Christian Friedrich Schönbein (1838).137 Subsequently, two systems, discharging mode (fuel cell) and charging mode (capacitors and batteries), have been examined independently in the last few decades, in alkaline (anion conducting membranes (AEM)) and acid media (PEM) (Fig. 5(c and d), respectively). Moreover, interest in the development of H2/O2 energy storage systems led to the combination of both systems, the so-called URC, which was reported by Mitlitsk et al. (Fig. 5(e)).138–140
Fig. 5 (a) Schematic illustration of an electrolytic/fuel cell. (b) Polarization curves of the H2- (red color curves) and O2-involved reactions (blue color curves), without a proper scale. This figure has been adapted/reproduced from ref. 128 with permission from the American Chemical Society. Innermost polymer electrolyte membrane (PEM) system (c) fuel cell (FC), (d) electrolyzer cell (EC), and (e) unitized regenerative cell (URC). Discharging mode: FC and recharging mode: EC. PEMFCs, where the fuel (e.g. H2) channel starting from one side, is split by the catalyst (e.g. graphene electrode) into H+ ions and e−. (f) More detail explanation of the current generation and combination of hydrogen ions and oxidant at the cathode, generating water and heat. This figure has been adapted/reproduced from ref. 135 with permission from the Royal Society of Chemistry. |
Electrochemical reactions play an essential role in energy conversion and energy storage systems, e.g., the ORR and HOR occurring on the cathode and anode of the H2–O2 fuel cell, whereas the OER and HER are the anodic and cathodic reactions in the electrolyzer, respectively.141 However, deep understanding of the basic principles and chemistry is required for understanding these future energy producing/storing devices. In the past few years, the advantages of fuel cell technology and its industrial application have drawn considerable attention especially in electric vehicles, and portable and residential power sources due to its high competence (specifically, 40–60% or up to 85% when dissipated heat is utilized), almost pollution free nature, minimal corrosion problems, and possible significant applications in many fields.142 The main advantages of fuel cells is that they are functionalized to produce electricity once a constant source of oxygen and fuel is applied. There are several types of fuel cells, including PEMFCs,143 solid oxide fuel cells (SOFCs),144 molten carbonate,145 and phosphoric acid,146 but they all share a cathode, anode, and electrolyte, which permit charges to travel in the fuel cell between two sides. A further detail explanation of the structure of the fuel cell, e.g., PEMFC, is shown in Fig. 5(f).
In fuel cells, rationally designed porous structures are the most advantageous electrocatalyst structures, especially hierarchical porous structures, where macro/meso-pores with a large number of electrocatalytically active sites make it the transportation of reactants and products (O2 and H2O) easy. However, major challenges exist regarding the synthesis of efficient electrocatalysts, and thus different synthetic methods, e.g., pyrolysis, are used. Another important factor regarding the efficiency of fuel cells is their durability, besides the activity of cathode catalysts, which can be enhanced either by synthesizing more stable electrocatalysts or anion exchange membranes (AEMs). Basically, AEMs consist of a polymer backbone with cationic ion-exchange groups that permit only the passage of OH−/H+ anions. The diffusion of OH− is considerably low compared to H+, which restricts the efficiency of fuel cells. An additional significant problem in fuel cells is their low durability, e.g., in AEMFCs due to the chemical degradation of OH− conducting ionomers. Therefore, suitable methods for membrane-electrode assembly (MEA) is of utmost importance for high performance. Despite the above-mentioned challenges, fuel cell technology holds great promise to reach the aim of developing affordable and sustainable power sources, and accordingly several breakthroughs have been made in fuel cell research in the last five years. Further increasing the stability, reducing cost, and elevating durability are hot topics in the development of a fuel cells, especially AEMFCs, which can work at a low-operating temperature.147 Similarly to batteries, fuel cells change fuel chemical energy in combination with oxidants into electric energy (Fig. 6).
Fig. 6 (a) Energy storage/conversion in the H2 and O2 cycle, where the energy storage two half-cell reactions in water splitting are the OER and HER. While for energy conversion, the half-cell reactions include the ORR and HOR. (b) Hydrogen electrocatalysis (HER, and HOR oxygen electrocatalysis (OER and ORR)) overpotential scheme. OER and ORR seriously limit the rate determination and energy efficiency in electrochemical energy devices because of the large η required to drive electrochemical reactions, even when utilizing benchmark catalysts in alkaline electrolytes. Anion-exchange membrane schematic representation (c) electrolyte and (d) alkaline electrolyte-based fuel cell. This figure has been adapted/reproduced from ref. 148 with permission from the Royal Society of Chemistry. |
The dissimilarity from batteries is that fuel cells do not require recharging if an oxidant and fuel are constantly provided. After the application of the fuel (hydrogen), the fuel cell produces mainly electricity, and some water and heat as a byproduct. In contrast to thermal engines, fuel cells have the advantages of enhanced effectiveness, no ecological contamination, and infinite availability of reactants. Thus, fuel cells are predicted to be applied in the future extensively for commercial uses in stationary, transportation and portable power production, and may facilitate overcoming the worldwide energy crisis with a clean atmosphere. Among the fuel cells, PEMFCs are actively utilized for moveable electronics, vehicles and combined heat and power (CHP) systems because of their elevated power density, low working temperature, and swift start-up. PEMFCs also have the advantages of low ecological impact and high energy density, and thus have gained worldwide attention. PEMFCs are particularly suitable as central power sources in automobiles and buses. Fuel cell vehicles (FCVs) are considered an absolute solution in the automotive trade and contain large scale benefits above battery-based electric vehicles (EVs). The first commercially produced FCVs, the Toyota Mirai, were commercially sold/leased in 2014, but with a high price due to the excessive Pt loading in the fuel cell.149 PEMFCs have already been widely recognized, and the commercialization of these devices is gradually progressing, but AEMFCs are also being intensively developed. As one of the subcategories of AEMFCs, alkaline direct methanol fuel cells (ADMFCs) are capable alternative power sources in portable electronics due to their easier storage and handling of alcohols compared to hydrogen.147 Fuel cells have many uses, such as in automobiles, backup power systems, smart phones, and smart textiles (implanted digital computing and electronics) since they present a long-lasting electricity supply. Hence, the worldwide increasing demand for clean energy has directed research efforts towards alternative and efficient energy conversion technologies. It is estimated that fuel cells will play an increasing role in future energy security. At peak times, it will be desirable to accumulate extra energy in the form of H2 gas (for example via water electrolysis), which allows one to store it for a longer time and then to convert the chemical energy into electricity in the fuel cell at any time.
Some automobile companies such as Toyota, Honda, and Hyundai have introduced fuel cell-based electrical vehicles (FCEVs). However, is a need to gain progress in this field for sophisticated energy exchange systems such as water electrolysis, metal–air batteries, and fuel cells. Considering the continuous decrease in fossil fuels and deterioration of the environment, it is significant and necessary to discover plentiful, environmentally friendly and renewable energy sources. The standard Nernst potential in the O2 half-cell reaction observed is 1.23 V versus RHE. This is determined by the standard H2 potential at a specific pH. At 0 pH, this potential is considered to be the standard hydrogen electrode (SHE) potential. The kinetically favourable ORR is preferentially below the half-cell potential since the OER lies above this potential. The ORR/OER kinetics is very poor due to immense deviation from the half-cell potential, which is called “η”, requiring considerable electricity. Thus, considering the physical response η is very important toward the further synthesis of proficient catalysts for the OER. The oxygen electrocatalysis benchmark catalysts in acidic solutions is Pt for the ORR (∼$50 per g) and ruthenium oxides for the OER (∼$2.50 per g). At practical i = 1.5 A cm−2geo, these benchmark electrocatalysts still have ∼0.4 V η. In contrast, in basic solution, the ORR is catalyzed using metal oxides that have earth abundant elements such as NPMCs (e.g. Fe: 0.0001 g−1; Co: 0.03 g−1; and Ni: ∼$0.02 g−1) and early rare-earth metals (e.g. La: ∼$0.02 g−1). Some good oxide catalysts contain MnOx for the ORR and Ni-Fe-Co oxides for the OER, needing an η of approximately 0.4 V at 10 mA cm−2geo, and recently our group introduced an rGO-coated abundant nature mayenite electride, 12CaO·7Al2O3 (C12A7:e−) for the ORR.1,148
Electrocatalysts have a major position in clean energy production, allowing a number of durable methods for potential pollution-free energy devices, such as fuel cells, electrolyzers and metal–air batteries. The basic reactions occurring in these devices are the ORR, OER and/or HER, which are multielectron-based electrochemical reactions with high η values and slow kinetics. Accordingly, many electrocatalysts are commercially available for these reactions, mostly based on PGMs; however, they are costly and unstable, and thus not feasible for common application. The requirements for an ideal type of electrocatalyst include easy preparation, very stable, inexpensive, and preferably have a bifunctional nature, which can simplify design protocols in rechargeable metal–air batteries and reproducible fuel cells/electrolyzers, and lead to the realization of more efficient but practical devices. The universal feedstock of abundant elements such as H2O, CO2, and N2 in the Earth's atmosphere provides a chance to convert the abovementioned products into energy by means of electrochemical methods coupled with inexhaustible energy if electrocatalysts with the required properties are developed. In the case of water splitting, hydrogen and oxygen evolution half-reactions occur, by which sustainable pure hydrogen can produced. Hydrogen is an energy carrier by which chemical energy can be transformed into clean electrical energy in fuel cells via the HOR and ORR reactions. Similarly, hydrogen peroxide can produced by the ORR, which is an important chemical used in the water treatment, pulp and paper bleaching industries.150 Similarly, atmospheric or open starting point sources of CO2 can also be used as fuel, commodity chemicals, final product chemicals, and starting materials for polymers through groundwork electric-reduction.151 Moreover, the electric reduction of N2 into NH3 can be utilized for the production of fertilizers and reduce distribution costs, which can be synthesized by the Haber–Bosch process and using the newly introduced C12A7:e− catalyst-based method.6 Based on these crucial ideas, the development of more active and stable electrocatalysts for further boosting the efficiency and selectivity in chemical conversions is important. Regarding improving activity and the reaction rate, increasing the number of active sites at a particular electrode or increasing the built-in activity of each active site are important to investigate. There are some physical limitations for catalysts loaded on the working electrode, which do not influence other significant factors, for example charge and mass transportation. In contrast, increasing inherent activity leads to a direct enhancement in electrode reactivity in a way that diminishes charge transport issues due to massive catalyst loadings, and with enhanced inherent reactivity, the catalyst components may be reduced, which can be cost-effective.2
The basic method for the manufacture of hydrogen gas is the HER. Water electrolysis is the electrochemical process traditionally used for the production of high-purity hydrogen. On the other hand, recently the HER acts as a connection between basic surface electrocatalysis and the behavior of recently reported materials. Thus, designing highly efficient HER electrocatalysts is necessary for clean energy production at the fundamental and utilization levels. Basically, the HER cathodic electrochemical reaction for water splitting is based on the 2e− transfer mechanism with one catalytic transitional species, H* (* = electrode active site), and it may occur via the Volmer–Heyrovsky or Volmer–Tafel mechanism.152,153 The overall HER rate is determined by the hydrogen adsorption free energy (ΔGH), which is strongly based on the electronic configuration of the electrocatalyst. The most favorable assessment of ΔGH = 0 eV, shows that the H2 bound on the catalyst surface may not be too strong nor too feeble. Noble metals, for instance, the benchmark Pt and Ru, have a ΔGH approaching 0 eV. In the case of earth-abundant metals, this value is either excessively high or excessively very low, resulting in insufficient catalytic activity, according to the prominent hydrogen volcano.154 Therefore, developing NPMCs that can replace Pt for the HER in acidic and basic electrolyte is an essential challenge. After catalyst activity, another important factor is lasting stability. The useful approaches to evaluate electrocatalyst durability include cyclic voltammetry (CV), stability tests that calculate the catalyst amount from drain to electrolyte, and utilizing thin-film-based catalyst configurations.153,155 Basically, in the HER, protons (acidic medium) or water molecules (alkaline medium) get reduced, and gaseous hydrogen is evaluated due to water splitting. Therefore, generally the HER process can be written as follows (catalyst active site “*”):
Acid medium
2H+ + 2e− → H2(g) | (1) |
Alkaline medium
2H2O + 2e− → H2(g) + 2OH− | (2) |
The HER process can also be explained further. The standard potentials (E°) depend on the nature of the active ions. Also, the HER consists of many electrochemical reactions, and the activation energy obstacle to advance the reaction is known as “η”. Consequently, materials having a lower η will have an enhanced reaction rate and high efficiency. The HER mechanisms are the Volmer–Heyrovsky and Volmer–Tafel mechanisms. The HER in an acid environment consists of the following steps:
1. Hydrogen atom adsorption (Had) due to the combination of electrons and protons on the electrode surface (proton release) via the Volmer reaction
H+ + * + e− ⇔ Had* | (3) |
2. Had atom interaction with electrons and protons resulting in electrochemical desorption (Heyrovsky reaction)
Had* + H+ + e− ⇔ H2 + * | (4) |
3. Two Had atoms coupling to result in the dissociative desorption of H2 vie the Tafel reaction
2*Had ⇔ H2 + 2* | (5) |
On the other hand, in alkaline electrolyte, owing to the OH− concentration, HER advances due to following steps:
1. Water molecule coupling to electrons, leading to the formation of Had atom at the electrode interface (Volmer reaction)
H2O + e− ⇔ Had* + OH− | (6) |
2. Had atom combines with molecular water and electron, allowing electrochemical H2 desorption (Heyrovsky reaction):
Had* + H2O + e− ⇔ H2 + * + OH− | (7) |
3. The Tafel reaction is similar to that previously discussed in acidic environment.
Hence, in alkaline and acid environment, for the HER, the formation of Had starts through the Volmer reaction (eqn (3) or (6)) while, the successive formation Had can continue through the Heyrovsky reaction (eqn (4) or (7)) or dissociative desorption through the Tafel reaction (eqn (5)). In contrast, the HER can be investigated via the Tafel slope obtained from the HER polarization curves. The Tafel slope analysis shows the native electrocatalyst properties, and the experiential magnitudes allow the mechanisms to be differentiated. In the Tafel rate-determining step (RDS), the Tafel slope is ∼30 mV dec−1, but for the Heyrovsky and Volmer RDS, the Tafel slopes are 40 and 120 mV dec−1, respectively. Co-based electrocatalysts follow the Volmer–Heyrovsky mechanism for the HER, where generally the Heyrovsky step is the RDS in an acid environment, whereas the Volmer step is the measured RDS in an alkaline environment.141,156,157
Acid medium
2H2O(aq) → 4H+(aq) + 4e− + O2(g) | (8) |
Alkaline medium
4OH−(aq) → 2H2O(aq) + 4e− + O2(g) | (9) |
In more detail, in the literature, various proposed OER mechanisms occur in acid (eqn (10)–(14)) and alkaline media (eqn (15)–(19)), involving MOH and MO intermediates. Fig. 7 shows two dissimilar pathways for the OER, where the black line shows that its progression involves the development of peroxide (M–OOH) intermediate and the green line shows the direct response for two nearby oxo (M–O) intermediates to produce an O2 molecule.141 In the green route, 2MO directly combines to make O2(g) (eqn (12)), including the generation of MOOH intermediate (eqn (13) and (18)), which later disintegrates via the black route to O2(g) (eqn (14) and (19)). In the heterogeneous OER process, all M–O bonding interactions in intermediates (MOH, MO and MOOH) are critical to the overall electrocatalytic activity.
Fig. 7 In loop blue line shows acidic and red line alkaline OER mechanism. This figure has been adapted/reproduced from ref. 141 with permission from the Royal Society of Chemistry. |
Acid medium
H2O + M → MOH + H+ + e+ | (10) |
MOH + OH− → MO + H2O + e− | (11) |
MO → 2M + O2(g) | (12) |
MO + HO → MOOH + H+ + e− | (13) |
MOOH + H2O(l) → M + O2 + H+ + e− | (14) |
Alkaline medium
OH + M− → MOH | (15) |
MOH + OH− → MO + H2O | (16) |
MO → 2M + O2(g) | (17) |
MO + OH− → MOOH + e− | (18) |
MOOH + OH− → M + O2(g) + H2O | (19) |
Thus far, the benchmark materials examined as potential NPMCs for the OER are Co- and Ni-based materials (both without support and with support on, e.g., carbon). Among the NPMCs for the OER, Co-based materials are only active in alkaline medium, e.g., CoOx and CoOOH.141 The limitation of the traditional OER process based on the four-step mechanism is that it is frequently restricted by its lethargic electron transfer kinetics. Fascinatingly, the progress of water oxidation catalysts with excellent performance and substantial stability is big challenge. Therefore, in the OER, noble-metal oxides, e.g., RuO2 and IrO2, are taken as the active catalysts, regardless of their tremendously high cost and shortage. Thus, the discovery of earth-abundant 2D materials as alternatives to replace the expensive Ir and Ru is required. Accordingly, substrate-free, transition metal-based catalysts and graphene show novel benefits due to their high charge transport, although the poor OER activity of graphene considerably prevents its general usage. Graphene-like TMDCs are also significantly active in the field of OER because of their novel structural-based electronic properties. However, previous studies showed that the most important shortcoming of layered TMDCs (WS2, MoS2, and MoSe2) is their low conductivity. Fortunately, BP is a metal-free layered semiconductor with comparatively elevated carrier transport (≥200 cm2 V−1 s−1) and substantial catalytic performances.
Direct 4e− pathway
(28) |
(29) |
Two-step 2e− pathway
(30) |
And
(31) |
O2 + 2H+ + 2e− → H2O2 (E° = +0.867 V vs. RHE, pH = 0) | (32) |
And
H2O2 + 2H+ + 2e− → H2O \cr \tab \quad (E° = +1.76 V vs. RHE, pH = 0), (acidic medium) | (33) |
In fuel cells, the ORR with the 4e− pathway is an advantageous cathodic method due to its positive potential pass on a large open circuit potential. To gain maximum energy from the reaction, the production of peroxide should be avoided. Peroxide formation not only decreases the operating potential and current efficiencies, but also causes the degradation of fuel cell components, which results in an unstable performance. Thus, in more detail, the ORR mechanism can be explained graphically. The typical ORR polarization curves can be divided into three regions, i.e., kinetically controlled region, diffusion-controlled region, and mixed kinetic/diffusion-controlled region (Fig. 8). JL denotes the diffusion limited i. In the kinetically restricted region the O2 reduction rate is very slow with a small increase in i with a decrease in potential. Conversely, a considerable boost in i is observed in the mixed kinetic and diffusion interface area. In this region, acceleration of the reaction occurs with a marked drop in potential value. In the diffusion interface region, i is determined by the rate at which the diffusion of the reactants occur. Quantitative analysis of a catalyst in terms of its activity can be done considering two parameters, i.e. the onset potential (Eonset) and half-wave potential (E1/2). The more positive the potential, more active the ORR electrocatalyst.163 However, fast and more suitable screening techniques to describe the electrochemical nature of recently synthesized materials are necessary. Among them, the rotating disk electrode (RDE) is the widely applied method to exemplify the properties of catalysts/catalysts maintained in liquid electrolyte, which was reported by Stonehart and Ross in 1976.164
Fig. 8 Characteristic ORR curve of a single catalyst. This figure has been adapted/reproduced from ref. 163 with permission from Elsevier. |
Generally, in RDE-based measurements, a thin film is fabricated, which was proposed by Gloaguen et al. (1994).165 Subsequently, the expansion of this method was reviewed by Schmidt and Gasteiger.166 Consequently, electrocatalysts are normally put in a water/alcohol/Nafion mixture to obtain a homogeneous ink. Actually, Nafion acts as a glue material to adhere the catalyst on the RDE during electrode rotation sets. There are two ways to add Nafion, either dilute Nafion solution is added directly to powder/solvent mixture at some stage in making the ink or just directly dropping Nafion on the catalyst film to make a Nafion coating. In either option, the quantity of Nafion must be very low to reduce extra dispersion, which limits the amount of O2 and decrease in IR content. Generally, the Nafion film thickness must be less the 0.2 μm when it is spread on the catalyst film or the quantity of solid Nafion must be less than 20 wt% of the catalyst weight when it is put in the ink. The obtained ink is deposited on glassy carbon electrode (GCE), which is known as the working electrode, to form catalyst films. To mollify the mass transfer in ORR activity measurements, the GCE is rotated at different rotation speeds to increase the O2 mass transfer rates at the electrode. The inherent ORR catalyst activity (kinetic current without mass transfer effect) can be calculated using the Koutecky–Levich (K–L) relation:
1/j = 1/jk + 1/jl,c = 1/jk + 1/0.62nFAC0 × D02/3v−1/6ω1/2 | (34) |
where, j, jk, and jl,c are the calculated i, kinetic i, and diffusion-limited i, respectively. The dispersion-limited current is calculated using the electron transfer (n) number, and Faraday constant (F), electrode geometric area (A), concentration of dissolved O2 in solution (C0*), diffusion coefficient of O2 (D0), kinetic viscosity of solution (v), and rotation speed of electrode (ω). jk is normally calculated using the K–L plot (j−1vs. ω−1/2) at different ration speeds. On the other hand, it can also be calculated using the calculated diffusion-limited current at fix rotation speed (normally, 1600 rpm). The K–L equation is ideally applied when the surface of the fabricated thin film on a GCE is smooth, under laminar flow hydrodynamics.
The catalyst thin film quality strongly affects RDE-based kinetic current measurements. Therefore, the fabricated catalyst film should be very thin, uniform, and very smooth. In case of thick films, mass transport is restricted throughout the film, which causes partial use of the catalyst. For PGM-based electrocatalysts, especially Pt, the thickness of the catalyst thin film on the GCE is generally calculated considering the quantity of carbon black. In the case of loading Pt electrocatalysts with different quantities, it must be considered that the weight of carbon black must remain the same on the working electrode because it will probably retain same the catalyst thin film thicknesses. This is noticeable in the case of a non-uniform, rough thin film surfaces for RDE calculations since the K–L relation is generally applied for a smooth film surface on the GCE. To get good and repeatable results by RDE calculations, the first thing is to make a well dispersed catalyst ink. However, due to the different surface natures of support materials used for catalytic measurements (such as carbon black or noncarbon) with electrocatalysts, it is not easy to develop a general method for all types of electrocatalysts. Besides, water and ethanol/isopropanol are generally use to make the ink because they help to keep the carbon surfaces wet. The water/alcohol ratios strongly depend on the type of carbon and electrocatalysts used, and also the pre-treatment conditions applied, such as annealing and acid treatment. Sometimes, even with a fine ink obtained using well-dispersed electrocatalyst and support material particles, a uniform catalyst film is not certain. The synthetic method, applied drying conditions, and environmental circumstances strongly affect the film quality. Therefore, Garsany et al.167,168 introduced new way to reproduce fine value electrocatalyst films by rotating the GCE at 700 rpm and adding the ink dropwise on top of the GCE surface. The normal variation was a lower electrochemical active area (ECA) and ORR value in the case when electrocatalyst films on GCE were prepared via the rotational drying process, and their activities were 70% greater than previously applied methods for the preparation of films. Similarly, Ke et al.169 also introduced the “intermittently micro-contact coating fine droplets” technique, in which the GCE can be uniformly covered with about 3 × 103 tiny droplets (3 nL per droplet). A solution for the reproducibility problem was also suggested by Shinozaki et al.,170 which when applied, they found that drying the electrocatalyst thin film in isopropanol gas atmosphere could help to make an excellent film with good ORR activity. The Nafion ionomer binder used for catalyst films has a negatively effect on the calculated catalytic activity. Therefore, due to the block effect using Nafion ionomer, the specific activity for Pt/C at 0.9 V was only 0.15 mA cm−2 Pt.149 However, most electrocatalysts have been evaluated using an RDE rather than real fuel-cell devices, i.e., MEA. It is well-known that good ORR activity from RDE tests cannot always transfer good MEA performance since the latter is in a much more complicated environment. Specifically, the RDE setup is more simplified and ideally designed to evaluate the intrinsic electrochemical property of electrocatalysts. Its limitation is also obvious, where it is not able to evaluate ion transfer, water management, architecture stability, and particularly mass transfer, which becomes even more severe for carbon-based metal-free cathodes due to their significantly thick catalyst layer. More importantly, there is very limited fundamental understanding of this gap because of much less effort placed on this.171
ηtherm = (Qreact − Qlat)/Qreact = ΔG/ΔH = Wel/Qreact | (35) |
ηv = Ui/E° = UiEMF | (36) |
Hence, the η of fuel cells can be written as, η = ηtherm × ηv. For H2/O2 fuel cells, the rate of EMF = E° is 1.229 V. In an actual system, it is important to know the total efficiency of the fuel cell, “ηtotal”, not only depends ηtherm, and ηv, but also on two other features, e.g., efficiency because of reactant consumption or coulombic efficiency, ηcoul; and the efficiency of design, ηdesign (ηtotal = η × ηcoul × ηdesign). All these relations show that to improve the fuel cell efficiency, firstly, we should increase the following three factors, first ηv (via electrocatalysis), second ηcoul (through selectivity/tolerance of materials), and third ηdesign (through system and layer by layer structure design).172 Here, we discuss in detail about the electrochemical properties of 2D materials, the ORR, OER, HER, etc. The integration of fuel cells in the electronics community faces numerous challenges, and a few of them are as follows:
1. Synthesis of suitable electrodes materials, especially also applicable to flexible electronics.
2. Replacement of high cost noble metals and their alloys as electrocatalysts, to make the fuel cell-based technology accessible to common people.
3. The requirement of decreasing metal electrode poisoning and enhancing their long term stability.
To solve these issues, a new category of electrocatalytic materials that is inexpensive, highly efficient and high sustainable need to be synthesized considering fuel cell-based technology is the future of humanity, where common people can benefit from it. Regarding the solution to the aforementioned issues, 2D materials are perfect to tackle these problems and make electron transfer easier in both fuel oxidation and ORR. Also, 2D materials should have potential because of proton membranes and their high hydrogen ion conductivity. This, together with water impermeability, methanol, and H2 may resolve the setback in future fuel cells and electrode poisoning. In next section, we explain in detail the electrochemical applications of 2D materials for energy production and storage.
Additionally, the tunable electronic properties of 2D materials can tailor their catalytic performance.175 In 2D nanomaterials with moderate electrocatalytic applications, the most important parameters are the adsorption energy (microscopic property), number of active sites, reaction kinetics (macroscopic property) and moderate conductivity. The obvious rates of catalytic response are mainly reliant on extrinsic properties (crystalline morphology, etc.) and inherent structure (work function, d-band position, etc.) of the catalyst, followed by the calculated transitional adsorption energies, kinetics, and reaction energy barriers. Considerable progress has been achieved in experimental nanotechnologies based on DFT simulation, and particularly their combination, resulting in new approaches for catalytic methods at the atomic/molecular level. Therefore, different synthetic protocols based on the modification of electrocatalysts such as heteroatom doping, tuning band gaps via defects, and interfacial band gap tailoring, have been applied to further boost the electrocatalytic activity of 2D nanomaterials. Fig. 9 shows how the electrocatalytic properties of 2D materials can be enhanced, which can be separated into two groups, growing active surfaces (top semicircle) and increasing intrinsic catalyst activity (base semicircle). These methods will help with detail explanations of proposed electrocatalyst principles for specific catalytic process. For example, for OER electrocatalysis, most metal-free nanosheets and 2D MXenes have been found to be promising alternatives to the benchmark catalysts.95
Fig. 9 Electrocatalysis scheme showing the importance of 2D materials as electrocatalysts. This figure has been adapted/reproduced from ref. 173 with permission from the Royal Society of Chemistry. |
After the introduction of graphene in 2004, other 2D materials have been investigated and their potential applications studied in various research fields.80 Therefore, Jiao et al. carried out DFT simulations to establish an accurate structural description, and electrochemical experiments to determine the origin of the ORR and HER activity in heteroatom-doped graphene. The formed doped graphene with B, N, O, S, and P showed different physicochemical characteristics and was applied to study the comparison between morphology, SSA, and defects of the different heteroatom-doped samples (Fig. 10(a and b)). Also a theoretical-based study was conducted based on different molecular models to compare computational outcomes with experimental ones (Fig. 10(c)). The calculated free energy illustration for the HER, ORR and experimentally measured visible catalytic activity are shown in Fig. 10(d–g). By comparing the proper theoretical and experimental results (e.g., ΔG and i0, respectively), a volcano plot can be built, which shows the intrinsic activity trend and gives predictions about catalysts with superior performances (Fig. 10(h and i)). Basically, the volcano plot is associated with the Sabatier principle, according to which the highest catalytic activity can obtained when the catalyst surface has suitable binding energies for reactive intermediates. Therefore, in the case of loosely bonded intermediates, surface activation and reaction are very difficult. On the other hand, in the case of the too strongly bonded intermediates, they may engage all existing active sites and poison them. Fig. 10(h and i) show that all the samples are situated on the right side of the volcano, which represents that intermediate adsorption is comparatively weak. Thus, obviously, if a catalyst has a ΔG value in the center of volcano (i.e., small ΔG value, stronger adsorption in doped graphene), it will own a larger i0 value from both computational and experimental sides. By relating the molecular orbital theory (MOT), a catalyst design method to modify graphene-based materials can be proposed (nearer to the volcano tip). Besides graphene, several new types of 2D materials have also been studied for electrocatalytic applications.176 In comparison to other electrocatalytic materials, the two very significant structures of the 2D family in catalysis are their variable and consistently exposed lattice planes and exclusive electronic levels. 2D materials such as MoS2 exhibit great activity for the HER since their H2 adsorption Gibbs free energy (ΔGH*) is merely 0.08 eV, which is lower than that of the benchmark Pt (ΔGH* ≈ 0.09 eV).160 Here, we discuss in detail the electrocatalytic applications of 2D materials.
Fig. 10 (a) TEM image of N2-graphene; with inset showing ID/IG and SSA values. (b) N2 K-edge NEXAFS spectrum. (c) Schematic heteroatom doping structure. (d) ORR and (e) HER free energy diagrams at equilibrium potential for diverse models. (f) ORR and (g) HER Tafel plots. (h) ORR and (i) HER volcano plots. This figure has been adapted/reproduced from ref. 128 with permission from the American Chemical Society. |
Fig. 11 (a) Schematic representation of graphene with defects. (b) HAADF figure of defective graphene with hexagons, pentagons, heptagons, and octagons, shown by orange, green, blue, and red colors, respectively. (c) AFM image of graphene with defects. (d) LSVs for pristine, N-doped and defective graphene for the ORR. This figure has been adapted/reproduced from ref. 181 with permission from Wiley-Blackwell. |
Regarding ORR applications, the electron-deficient B atoms assist the chemisorption of negatively polarized O atoms and encourage O–O splitting. Recently, Qiao et al. doped B in graphene and showed that it exhibited the lowest ORR activation barrier free energy, resulting in its higher catalytic activity compared to P- and S-doped graphene.178 The O2 adsorption on B-doped graphene was calculated to continue via the end-on dioxygen intermediate, form elevated energy intermediate for dissociative system (I) and hydroperoxo species for associative method (II),192 which can largely activate the ORR process. On the other hand, the B content also had a limit, if the B content >140 ppm, then it causes a decrease in ORR activity in alkaline electrolyte due to an increase in the electrical resistivity.193
On the other hand, N2 doping provokes polarization in graphene sheets because of the higher electronegativity of N2 (3.04) than that of C. The N doping causes opening of the band gap, which makes it semiconducting.194 There are three types of primary N-bonding configurations observed, basal plane pyrrolic N, quaternary N (central graphitic), and edge pyridinic N. The sp2-coupled quaternary N and pyridinic N have a small influence on graphene, while sp3-hybridized pyrrolic N causes a distortion in the planar geometry of graphene. N-Doped graphene semiconductor is related to the N-bonding configuration. For quaternary N, the three (sp2) valence electrons of N2 form three bonds with adjacent carbon atoms, one (2pz) orbital occupied in p hybridization, and the fifth one is delocalized in the CB. The n-type doping in quaternary N leads to the donation of 0.5 electrons to the graphene lattice.195 On the other hand, pyridinic and pyrrolic N-atoms located at defective vacancies cause p-doping by withdrawing graphene electrons.196 The carrier concentration increased to 2.6 × 1013 cm2 (4 times that of pure graphene) using a very small doping concentration, i.e. 0.6% atom of quaternary N-doping.197 Dai et al. synthesized N-doped graphene and studied its ORR catalytic activity (2010),198,199 and showed that the quaternary N and pyridinic N fraction can make the ORR activity possible. For further explanations about these doping, Ruoff et al. further studied doped graphene in detail. They suggested that the quaternary N can only settle the limiting i, but on the other side, the pyridinic N can mainly develop the onset potential in the case of electrochemical ORR activity.200 Quaternary N may boost the density of electrons by p-orbital delocalization (n-type doping), causing improved nucleophile strength in the adjacent carbon, and hence will promote the O2 adsorption, which will further increase the ORR activation in alkaline electrolytes.200 In a more detailed study, it was suggested that the pyridinic N can further enhance the active surface for boosting the ORR activity in N-graphene.201 The importance of N-doping was also studied by Zhang et al., where they applied g-C3N4 as a template and source of N2 to synthesize highly porous N-doped graphene, with ORR catalytic activity analogous to the benchmark Pt/C in alkaline solution.182
Similarly to the previous studies on doping in graphene, the effect of P doping in graphene was also studied, which causes sp3 hybridization and tetrahedral analogue structure along three available carbon atoms. Due to the longer bond length (1.77 Å) of P–C than C–C, the P atoms stick out from the graphene plane by around 0.09 Å, where the P electronegativity (2.19) decreases compared to that of C (2.55), and the C–P polarity is opposite to the C–N bond. When contrasting to 0.5 electron, quaternary N relocates to the p orbitals, and P just donates about 0.2 electrons (Hirshfeld population analysis).202 Thus, due to the weaker electron contribution than N-doped graphene, P-doped bilayers show n-type behavior, and more stable in an O2 environment than n-type doping resultant from N-doping, and with higher mobility. P-Doped graphene has a high ORR performance with a more positive η and higher i.203 To explain this theoretically, Wu et al. demonstrated that the ORR occurs through the 4e− pathway.204 Similarly, recent DFT simulations by Yang et al. suggested that the ORR on P-doped graphene first occurs through the 2e− pathway, forming OOH transitional species, and then through the 4e− pathway via the O–O bond in OOH. The second step concerning OH reduction to H2O, with an energy barrier of 0.88 eV (very large), is suggested to be the rate-limiting step. The P atoms and their nearby C atoms are recommended to be more suitable as adsorption and active sites in the ORR. The P atoms may also function as a bridge, promoting electron transfer from C to adsorbates.205
Besides the aforementioned nonmetallic dopants in graphene, it is also very significant to study and discuss the doping of metals, e.g., Fe, Co, Ni, and Cu,206,207 and metal oxide (MnO2)208 in graphene, which may further boost the catalytic activity of doped graphene. Therefore, Pumera et al. doped a small quantity of Mn (18 ppm, 0.0018 wt%), which led to a potential shift by 80 mV in the ORR.208 Similarly, substitution of multiple hetero-atoms (B/N, S/N, P/N, I/N, B/N/P, N/P/S, N/Fe, and N/S/Fe) may modify the electronic structure of graphene.95 This can create a synergistic effect and engineer new properties. Fig. 12 shows Co-doped B in N-doped graphene, and Qiao et al. showed it will be very helpful in the ORR in alkaline solution. Fig. 12(a) shows that the η of B, and N-doped graphene is about 110 mV and E1/2 about 150 mV, and that Pt/C has a negative value.176 Co-Doping with S and N further affected the electron donor and acceptors,209 and enhanced the active sites.210 Zhou et al. demonstrated a codoping effect using aminothiazole-based Fe, S, O, and N co-doped graphene, AT-N/Fe-G,183 with a doping quantity of 6.5%, 6.8%, 1.5%, and 2.6% for N, O, S, and Fe, respectively. These co-doped graphene showed good ORR activity in both alkaline and acidic conditions compared to that of the benchmark 20% Pt/C. Regarding the stability in acidic solution, the potential of AT-N/Fe-G shifted positively by 12 mV after 10000 cycles, but up to date, the potential of Pt/C shifted by 58 mV.
Fig. 12 (a) LSVs curves at rotating rate of 1500 rpm, in O2-saturated alkaline solution (0.1 M KOH scan rate: 10 mV s−1). (b) HO2 adsorbed on B, N-graphene optimized arrangement structure. This figure has been adapted/reproduced from ref. 176 with permission from John Wiley and Sons Ltd. |
Theoretical studies based on DFT calculations show that dopant and doping vacancies strongly affect the ORR activity. A theoretical method based on free energy adsorption intermediates (O*, OH*, and OOH*) was introduced to calculate the ORR activity on graphene and its derivatives.178 Generally, DFT simulation-based experimental investigations for ORR activity of mono-element doping in graphene display the following trend: N-graphene ≥ B-graphene > P-graphene > S-graphene. In the case of di-atom-doped graphene, the activity changes as (P,N)-graphene > (B,N)-graphene > (S,N)-graphene. The best metal-free 2D catalyst demonstrated to date is N-doped porous graphene with a higher η of 0.076 V versus Ag/AgCl and limiting i of 5.79 mA cm−2 than that of the benchmark 20 wt% Pt/C (0.105 V vs. Ag/AgCl; 5.72 mA cm2). However, the ternary, Fe, S, and N-doped graphene still shows better ORR performances than N-doped catalysts. The E1/2 = 0.926 V vs. RHE for this ternary doped graphene is higher than that of Pt/C (0.907 V vs. RHE).95
Moreover, the addition of two or multiple components endows enhanced properties possibly due to the synergistic effect, and thereby improves the performances of MoS2-based composite materials. In addition, doping of MoS2 with other materials or atoms can alter the basal plane or enlarge the MoS2 interlayer spacing in nanosheets, and thus modify their d-band electronic properties. Here, first we discuss the most promising applications of the graphene family in composites with mayenite electride, C12A7:e−, for the ORR in fuel cells, which was introduced for the first time by our group.1 In this work, Karim Khan et al., for the first time, initiated a facile synthetic method,1,6–9 and the synthesized nanosized C12A7−xSnx:e−, where the doping level, x = 0.20 to 1 (thereafter, Sn-doped C12A7:e−) composite exhibited a high SSA of 244 m2 g−1 and increasing conductivity of 24 S cm−1, 68 S cm−1, 190 S cm−1 and 290 S cm−1, with an Sn doping level of x = 0.2, 0.4, 0.8, and 1, respectively. Fig. 13(a) show the CV curves, where the cathodic reduction currents revealed obvious diffusion-controlled ORR peaks at 0.78 V, and 0.8 V (vs. RHE) in O2-saturated 0.1 M KOH solution for Sn-doped C12A7:e− with and without rGO coating, respectively. The diffusion-limited i increased with an increase in the rotation rate, as usually observed (Fig. 13(b and c)). These results demonstrate the adequately high intrinsic activity without/with rGO coating for the Sn-doped C12A7:e− composite in the ORR.
Fig. 13 (a) Sn-doped C12A7:e− CVs curves (with and without rGO coated, with doping level, x = 1) at 50 mV s−1 scan rate in O2-saturated 0.1 M KOH electrolyte. Sn-doped C12A7e− (x = 1) (b) without rGO coating and (c) with rGO coating in O2-saturated 0.1 M KOH with various rotation rates at a scan rate of 10 mV s−1. This figure has been adapted/reproduced from ref. 1 with permission from the Royal Society of Chemistry. |
Fig. 14(a) shows that the i of Sn-doped C12A7:e− is approximately similar to that of the benchmark 20% Pt/C, but that of the rGO-coated sample is about 5.9 mA cm−2 higher. This may be because of the work-function difference between rGO and Sn-doped C12A7:e−, which is the reason for the transfer of electrons from the Sn-doped C12A7:e− to rGO.211 Consequently, the Sn-doped C12A7:e− composite with moderate conductivity, led to a low resistance and higher binding energy (BE), minimizing the unfavorable adsorption of intermediates, resulting in enhanced catalytic properties.212 These results demonstrate satisfactory high intrinsic activity with rGO-coated Sn-doped C12A7:e− in the ORR. The suggested number of exchanged electrons was 4, representing that the rGO-coated Sn-doped C12A7:e− composite followed the four-electron pathway, with water as the byproduct (Fig. 14(b)).213–215Fig. 14(c) shows that the relative current without and with rGO-coated Sn-doped C12A7:e− did not change after 11 h, and about a 59% decrease was observed for the 20% Pt/C electrode.213 Furthermore, the rGO-coated Sn-doped C12A7:e− exhibited no change in its LSV curve upon the addition of CH3OH.213,216
Fig. 14 (a) ORR polarization curves of without/with rGO coated Sn-doped C12A7:e−, and Pt/C, in O2-saturated 0.1 M KOH, @1600 rmp, and sweep rate of 10 mV s−1. (b) rGO coated Sn-doped C12A7:e− K–L plots, and (c) Sn-doped C12A7:e− (with and without rGO coating, doping level, x = 1) chronoamperometry and Pt/C curves. This figure has been adapted/reproduced from ref. 1 with permission from the Royal Society of Chemistry. |
Although, a higher temperature generally supports the removal of free O2 species from C12A7 and improve the electrical conductivity of materials, in anticipation of certain edges.218 Zhou et al. theoretically calculated ORR measurements for N-doped graphene supported by MXenes (Ti2C, V2C, Nb2C, and Mo2C monolayers) (Fig. 15(A) and (B)). Graphitic sheets hybridized with V2C (G/V2C) and Mo2C (G/Mo2C) MXene monolayers had remarkable ORR catalytic activities. The most active site was the hollow-site C atoms located close to the pyridinic N dopants, giving rise to the low overpotentials of 0.36 V (G/V2C) and 0.39 V (G/Mo2C) (Fig. 15(C)), which are competitive to the benchmark Pt (0.45 V). The strong electronic coupling between MXene and graphitic materials as a potential multifunctional hybrid is promising for flexible energy devices. Additionally, Xie et al.219 experimentally calculated the ORR activity for Pt and Ti3C2Tx composite, which were stable even after 10000 cycles, and the half-wave potentials were constant (Fig. 15(D)), but Pt/C lost 21 mV (Fig. 15(E)). By compositing alkalization-intercalated MXene (Ti3C2(OH/ONa)2) with Ag composites, the urchin-structured MXene-Ag0.9Ti0.1 bimetallic nanowire featured pronounced electrocatalytic activity for the ORR (Fig. 15(F)),220 with onset- and half wave-potentials higher than that of the benchmark 20 wt% Ag/C and Ag-based catalysts.217Table 1 describes the details about the ORR properties of graphene-based 2D materials in alkaline medium. rGO-coated Sn-doped C12A7:e− exhibited increased catalytic activity for the ORR mostly due to the balanced relationship between the number of active sites and the electron transportation and stability. With increasing heat-treatment, the graphitization degree increased in the rGO-coated on Sn-doped C12A7:e− electride, which is beneficial for the transportation of electrons.
Fig. 15 MXenes as catalysts for the ORR. (A–C) Four different models of N-doped graphene on (A) V2C and (B) Mo2C MXene monolayer. (C) Volcano plot of ORR overpotential and ηORRvs. OH* binding energy. The black box indicates the dimensions of the supercells. The red numbers label the C sites. C, N, O, V, and Mo atoms are denoted by gray, blue, red, green, and turquoise, respectively. ORR current on (D) Pt/Ti3C2X2 and (E) Pt/C catalysts before and after the accelerated durability tests. (F) ORR polarization curves of MXene-Ag0.9Ti0.1 nanowire and other control samples in 0.1 M KOH. This figure has been adapted/reproduced from ref. 217 with permission from Elsevier Inc. |
Catalyst | Synthetic method (precursor) | Electrochemical performance | Ref. |
---|---|---|---|
Defect graphene | Ar plasma etching (GO) | Onset/half-wave potential: 0.912 V, 0.737 V, electrons transferred: 3.85 (0.4 V vs. RHE) | 221 |
Annealing at 1150 °C in N2 (N-doped graphene) | Onset/Half-wave potential: 0.91 V, 0.76 V electrons transferred: 3.87 (0.50–0.65 V, vs. RHE) | 181 | |
B-Doped graphene | Annealing at 900 °C in Ar (GO and B2O3) | Onset potential: 0.029 V i @ 0.3 V: 2.38 mA cm−2 (vs. NHE) | 178 |
Annealing at 1200–1600 °C in N2 (B4C) | Overpotential difference relative to 20% Pt/C: 59 mV n: 3.6–4.2 (0.3 to 0.7 V vs. Ag/AgCl) Limiting current: 16.32 mA cm−2 (0.6 V vs. Ag/AgCl) | 222 | |
N-Doped graphene | Fluorination and then annealing at 800 °C (reduced GO, XeF2, and NH3) | Onset potential: 0.05 V (vs. Ag/AgCl) i @ 1600 rpm: 4.99 mA cm−2 | 180 |
Annealing at 700 °C (reduced GO and NH3) | Onset potential: 0.91 V, i @ 300 rpm: 1.4 mA cm−2 at 0.5 V (vs. RHE) | 223 | |
Hydrothermal method (180 °C) (GO and NH3) | Onset potential: 0.08 V, i @ 0.3 V: 1.67 mA cm−2, electrons transferred: 3 (0.32 to 0.7 V) vs. Ag/AgCl | 224 | |
Porous N-doped graphene | Annealing at 800–1000 °C in the mixture of H2 and Ar (porous Ni sheets and pyridine) | Onset potential: 0.08 V, i @ 0.4 V: 8.2 mA cm−2n: 3.9 (0.4 V) vs. Ag/AgCl | 225 |
N-Doped porous carbon nanosheets | Hydrothermal treatment and then annealing at 900 °C in N2 (porous g-C3N4 and glucose) | Competitive with Pt/C catalysis Onset/Half-wave potential: 0.076 V, 0.21 V, electrons transferred: 3.8–3.92 (0.35 to 0.6 V vs. Ag/AgCl) Tafel slope: 143 mV dec−1 | 182 |
P-Doped G | Pyrolysis at 1000 C in Ar (GO and triphenylphosphine) | Onset potential: 0.92 V, electrons transferred: 3.0–3.8 (0.5–0 V) vs. RHE | 203 |
S-Doped G | Cycling of S-G composite-based Li–S batteries (S-G composite) | Onset/Half-wave potential: 0.15 V, 0.37, ORR reduction peak: 0.34 V, V (vs. SCE), electrons transferred: 3.13 (0.8 V) vs. SCE | 226 |
Mn/reduced GO | Staudenmaier oxidation then hydrazine reduction | Onset potential: 183 mV (vs. Ag/AgCl) | 208 |
B/N co-doped graphene | Ionic liquids in a Wurtz-type reductive coupling (CCl4, K, and C1C4ImBF4) | ORR reduction peak: 0.314 V, half-wave potential @ 1600 rpm: 86 mV relative to 40% Pt/C, i @ 1600 rpm: 5.2 mA cm−2, electrons transferred: 3.84–3.88 (0.8 to 0.4 V vs. SCE) | 227 |
Two-step pyrolysis (dicyandiamide, GO, and boric acid) | Onset/Half-wave potential (@ 900 rpm) 0.86 V, 0.61 V, ORR reduction peak: 0.61 V, (vs. RHE), mass activity @ 0.75 V: 0.53 mA mg−1 | 209 | |
S/N co-doped graphene | Pyrolysis at 800 °C in Ar (dopamine-cysteine-reduced GO) | Onset potential: 0.002 V, ORR reduction peak potential: 0.213 V, electrons transferred: 2.98 to 3.36 (0.30 to 0.60 V) vs. Hg/HgO | 228 |
P/N co-doped graphene | Two-step pyrolysis (dicyandiamide, GO, and phosphoric acid) | Onset/Half-wave (@ 900 rpm) potential: 0.87 V, 0.64 V, ORR reduction peak: 0.64 V, mass activity @ 0.75 V: 0.80 mA mg−1 (vs. RHE) | 209 |
(P/N co-doped graphene)/carbon nanosheets | Pyrolysis (GO, polyaniline and phytic acid) | Onset/Half-wave potential: 1.01 V, 20 mV, ORR reduction peak: 0.85 V, i @ 0.6 V: 5.56 mAcm−2; Tafel slope: 51 mV dec−1n: 3.96 | 229 |
I/N co-doped graphene | Annealing at 900 °C in Ar (iodine, poly(aniline), and activated G) | Onset potential: 0.08 V, i @ 0.7 V: 11.76 mA cm−2, electrons transferred: 3.93 (0.7 V vs. Ag/AgCl) | 230 |
I/P co-doped graphene | Reflux heating (GO and PI3) | Onset potential: 0.202 V (vs. Ag/AgCl) | 231 |
N/P/S tri-doped graphene | Pyrolysis at 800 °C in N2 (GO, triphenylphosphine, and thiourea) | Onset potential: 0.03 V, ORR peak: 0.20 V (vs. Ag/AgCl) i: 6.41 mA cm−2 | 232 |
N/P/F tri-doped graphene | Pyrolysis at 950 °C (polyaniline, GO, and NH4PF6) | Half-wave potential: 40 mV, electrons transferred: 3.85 (0–0.6 V vs. RHE) | 233 |
N/Fe co-doped graphene | Pyrolysis and activation (Ferrocene, rGO, and NH3) | Onset/Half-wave potential: 0.09 V, 0.29 V, i: 15.7 mA cm−2, n: 3.83 (0.25 V) vs. Ag/AgCl | 234 |
N/S/Fe tri-doped graphene | Polymerization and then annealing at 900 °C in Ar (carbon black, 2-aminothiazole, and FeCl3) | Competitive with Pt/C under acidic and alkaline conditions Half-wave potential: 0.926 V, mass activity @ 1.0 V: 0.56 A g−1n: 3.97(vs. RHE) | 183 |
P-MoS2 | Thermolysis process (280 °C) (C3N4, MoO3, thiourea, and red phosphorus) | Onset/Half-wave potential: 0.96 V, 0.80 V, electrons transferred: 3.4–3.6 (0.6 to 0.8 V vs. RHE) | 235 |
O-MoS2 | Pyrolysis synthesis followed by O incorporation (ammonium molybdate, thiourea, melamine, and H2O2) | Onset/Half-wave potential: 0.94 V, 0.80 V (vs. RHE), i: 3.49 mA cm−2 | 236 |
MoS2/N-doped porous carbon nanosheets | Functionalization and then pyrolysis (4-iodophenylfunctionalized MoS2 and conjugated microporous polymer) | Half-wave potential: 0.14 V, diffusion-limited i @ 1600 rpm: 5.4 mA cm−2, i: 26.3 mA cm−2, electrons transferred: 2 4.0 (0.4 to 0.9 V) vs. Ag/AgCl | 237 |
MoSe2/G | In situ solvothermal route (280 °C) (rGO, MoCl5, and Se) | Onset potential: 0.02 V, i @ 0.8 V: 5.35 mA cm−2, electrons transferred: 3.40 (0.2 V) vs. Ag/AgCl | 238 |
WSe2/graphene | Thermolysis process (280 °C) (reduced GO, WCl6, and Se) | Onset potential: 0.02, i @ 0.8 V: 5.26 mA cm−2n: 3.57 (0.2 V) vs. Ag/AgCl | 235 |
Co3O4/N-doped graphene | Hydrolysis and oxidation and then hydrothermal crystallization (GO, cobalt acetate, and NH4OH) | Half-wave potential: 0.83 V, electrons transferred: 4.0 (0.60 to 0.75 V vs. RHE) Tafel slope: 42 mV dec−1 | 239 |
Carbon nanotubes/N-doped graphene | In situ doping during CVD growth (CH4 and NH3) | Onset potential: 0.88 V, electrons transferred: 3.22 (0.47 vs. RHE), limiting i: 132.8 mAcm−2 | 240 |
Catalyst | Electrochemical performance | Ref. |
---|---|---|
Defective graphene | Potential: 1.57 V (@ 10 mA cm−2vs. RHE) Tafel slope: 97 mV dec−1 | 181 |
N-Doped graphene | i (@ 0.7 V vs. Ag/AgCl): 0.27 mA cm−2 | 224 |
P-Doped graphene | Onset potential: 1.48 V, potential @ 10 mA cm−2: 1.56 V (vs. RHE) | 242 |
(P/N co-doped graphene)/carbon nanosheets jORR, 3 mA cm−2: 0.71 V | Onset potential: 1.32 V, i: 70.75 mA cm−2(@ 0.9 V), potential: 1.57 V (@ 10 mA cm−2vs. RHE) Tafel slope: 70 mV dec−1, Difference in potential between jOER, 10 mA cm−2 and jORR, 3 mA cm−2: 0.71 V | 229 |
N/P/F tri-doped graphene | Onset potential: 1.62 V (vs. RHE) Tafel slope: 136 mV dec−1 | 233 |
La+ Ni-Fe LDH nanosheets | Over-potential: 0.26 V (@ 10 mAcm−2) | 246 |
Ni–Fe LDH nanosheets | Overpotential: 0.305 V (@ 10 mA cm−2), i: 9.5 mA cm−2 (@ h = 300 mV), Tafel slope: 40 mV dec−1 | 247 |
Ni–Co LDH nanosheets | Overpotential: 0.335 V (@ 10 mA cm−2), i: 1.5 mA cm−2(@ h = 300 mV), Tafel slope: 41 mV dec−1 | 247 |
Co–Co LDH nanosheets | Overpotential: 0.355 V (@ 10 mA cm−2), i: 0.85 mA cm−2 (@ h = 300 mV), Tafel slope: 45 mV dec−1 | 247 |
Ni–V LDH nanosheets | i: 27.0 1.6 mA cm−2 (@ h = 350 mV), Tafel slope: 50 mVdec−1, TOF: 0.054–0.003 s−1 | 248 |
TiO2/NiO nanosheets | Overpotential (h): 0.32 V, potential: 1.55 V (@ 10 mA cm−2, vs. RHE), Tafel slope: 52 mV dec−1 | 249 |
Porous b-Ni(OH)2 nanosheets | Overpotential: 0.415 V (@ 10 mA cm−2) Tafel slope: 60 mV dec−1 | 250 |
Ni3N nanosheets | Mass activity: 572 A g−1, overpotential: 0.35 V (@ 52.3 mA cm−2), Tafel slope: 45 mV dec−1 | 251 |
NiPS3@NiOOH nanosheets | Onset potential: 1.48 V (vs. RHE) Overpotential: 0.35 V (@ 10 mA cm−2), Tafel slope: 80 mV dec−1 | 95 |
Porous Co3O4 nanosheets | Onset potential: 1.51 V (vs. RHE), i: 10 mA cm−2 at 1.56 V (vs. RHE) Tafel slope: 69 mV dec−1 | 240 |
NiCo2O4 nano-sheets (O vacancies) | Overpotential: 0.32 V (@ 10 mA cm−2), i: 285 mA cm−2(@ 0.8 V, vs. Ag/AgCl), Tafel slope: 30 mV dec−1 | 252 |
g-CoOOH nano-sheets | Onset potential: 1.47 V (vs. RHE), overpotential: 0.30 V(@ 10 mA cm−2), mass activity: 66.6 A g−1, (@ 0.3 V, vs. RHE) Tafel slope: 38 mV dec−1 | 253 |
Ni/N-doped graphene | Overpotential: 0.32 V(@ 10 mA cm−2), i: 16.3 mA cm−2 (at an overpotential of 0.4 V) | 83 |
Co3O4/N-doped graphene | Overpotential: 0.31 V (@ 10 mA cm−2), Tafel slope: 67 mV dec−1 | 239 |
Carbon nano-tubes/N-doped graphene | Potential: 1.63 V (@ 10 mA cm−2, vs. RHE), Tafel slope: 83 mV dec−1 | 240 |
g-C3N4/graphene | Onset potential: 0.58 V (vs. Ag/AgCl), over-potential: 0.539 V (@ 10 mA cm−2) | 254 |
g-C3N4/Ti3C2 nanosheets | Onset potential: 1.44 V (vs. RHE), potential: 1.65 V (@ 10 mA cm−2, vs. RHE), i: 13 mA cm−2(@ 1.70 V, vs. RHE), Tafel slope: 74.6 mV dec−1 | 255 |
Few-layer BP | Onset potential: 1.45 V (vs. RHE), Tafel slope: 88 mV dec−1 | 14 |
BP bulk@Ti foil | Onset potential: 1.48 V (vs. RHE), Tafel slope: 91.5 mV dec−1 | 29 |
IrO2 | Onset potential: 1.5 V (vs. RHE), Tafel slope: 49 mV dec−1 | 256 |
N-Graphene/NiCo | Onset potential: 0.7 V (vs. Ag/AgCl), Tafel slope: 614 mV dec−1 | 257 |
RuO2 | Onset potential: 1.42 V vs. RHE, Tafel slope: 76 mV dec−1 | 258 |
Co3O4/N-graphene | Onset potential: 1.5 V (vs. RHE), Tafel slope: 67 mV dec−1 | 239 |
N,O-Graphene/CNT | Onset potential: 0.315 (vs. Ag/AgCl), Tafel slope: 141 mV dec−1 | 259 |
Co3O4 | Onset potential: 1.5 V (vs. RHE), Tafel slope: 234 mV dec−1 | 260 |
Pt/C | Onset potential: 1.5 V (vs. RHE), Tafel slope: 317 mV dec−1 | 261 |
Similarly, γ-CoOOH sheets with a thickness of 1.4 nm were used as an electrocatalyst for water oxidation, showing high mass activity (66.6 A g−1), which was 20-fold greater than that of bulk γ-CoOOH and 2.4-fold superior to that of the state-of-the-art IrO2,264 with the best OER activity among the NNMCs sheets published in the literature. This elevated activity is mainly due to their semi-metallic properties, resulting in improved H2O electrophilicity and quicker interfacial transfer of electrons between Co ions and adsorbed –OOH intermediate species to make O2. The NiOOH species resulting from the oxidation of Ni, and Ni–O–C or Ni–N–C (between Ni and N-doped graphene bonding) are active centers, which promote the oxidation of OH to O2. The porous g-C3N4 and MXene (Ti3C2) composite consists of double-sheets, which enhance the O2 electrochemistry, giving a sharp onset potential of 1.44 V (vs. RHE) and OER currents higher than that of the benchmark 20 wt% IrO2/C.255 This electrode offers a significant synergistic effect in comparison to the separate g-C3N4 and Ti3C2 electrodes. BP has gained significant attention due to its novel properties, which also result in enhanced electrocatalytic performances. Recently, Wang et al. synthesized bulk BP, which had an onset potential of 1.49 V (vs. RHE) and a low Tafel slope of 72.88 mV dec−1, showing comparative electrocatalytic characteristics to the RuO2 benchmark catalysts for the OER, but its bulk crystal structure still suffers from low active sites. Recently it was reported that a decrease in the number of layers in 2D materials will result in further exposed electrocatalytic active surface and amplified SSA than that of their bulk counterparts. These favorable results were first observed with the native ultrathin lamellar BP structure, which can open up alternative paths in using 2D layered materials with high electrocatalytic activity. Xiaohui Ren et al. used BP for the OER as a new rediscovered layered material (Fig. 16).14
Fig. 16 (a) Scheme of synthesis of BP nanosheets; (b) low- and (c) high-resolution SEM images of as-prepared BP nanosheets on SiO2/Si substrate; (d) TEM image (inset in (d): SAED pattern); (e) HR-TEM image; (f) Raman spectrum; (g) tapping-mode AFM image and comparative cross-sectional analysis (red, green, and blue solid line area). This figure has been adapted/reproduced from ref. 14 with permission from Wiley-VCH Verlag. |
Significantly, a decrease in the thickness of BP nanosheets will create extra active sites and boost their electrocatalytic activities.14Fig. 17 demonstrates the measured polarization curves of 2D BP under different C[OH−] electrolytes (1, 0.2, 0.1, and 0.05 M) in 1 M KOH, which showed improved i and lower η = 1.45 eV than the previously obtained 1.48 eV in bulk BP.14 Interestingly, the long-term robustness of the BP nanosheets was demonstrated in 1, 0.2, 0.1, and 0.05 M KOH. The i values of BP in 1 M KOH is shown in Fig. 17(a), which fluctuate near 6.2 mA cm−2, and becomes much larger than that in 0.2, 0.1, and 0.05 M KOH. After 10000 s, no noticeable decay in the i was observed in the 2D BP in the presence of 1, 0.2, and 0.1 M KOH, while a slight decline was detected in 0.05 m KOH electrolyte. In addition, linear sweep voltammetry (LSV) of the BP nanosheets in 1, 0.2, 0.1, and 0.05 M KOH was carried out. The stability before LSV and after LSV was studied for comparison, as shown in (Fig. 17(b)). The LSV curves demonstrate the high stability of the BP nanosheets in KOH solution, and comparatively, i in 1 M KOH is rather steady at 8 mA cm−2 (under a potential of 1.8 V) presenting no noticeable decay.14Fig. 17(c) shows the long-term stability of the BP nanosheets (Fig. 18).
Fig. 17 (a) Polarization curves (LSV) and inset corresponds to the Tafel plots of the BP nanosheets in KOH electrolyte with different concentrations (1, 0.2, 0.1, and 0.05 M). (b) Tafel slope at 10 mV cm−2 as a function of the concentration of OH− (C[OH−]).14 (c) Long-term stability of BP nanosheets in different basic electrolyte concentrations of 1, 0.2, 0.1, and 0.05 M KOH during the OER process. (d) LSVs of BP nanosheets in 1, 0.2, 0.1, and 0.05 M KOH before and after 10000 s measurement. This figure has been adapted/reproduced from ref. 14 with permission from Wiley-VCH Verlag. |
Fig. 18 (a) Synthesis of N and P co-doped porous foam, (i) aniline, (ii) phytic acid, and (iii) complex formation, via oxidative polymerization to form a 3D PANi hydrogel cross-linked with phytic acid. SEM images of (b) PANi aerogel and (c) NPMC-1000. Inset in c: PANi aerogel digital photograph before (left) and after (right) pyrolysis at 1000 °C. (d) HR-TEM image and (e) TEM image (e, left), with equivalent element NPMC-1000 mapping images. The element mapping images in the TEM analysis for C, N and P show a uniform element distribution. This figure has been adapted/reproduced from ref. 261 with permission from Nature Publishing Group. |
Jingfang Zhang et al. studied the OER activity of transition metal-based electrocatalysts, especially for single precious metal atoms supported on LDHs. The single-atom Au supported on NiFe LDH (sAu/NiFe LDH) had overpotential of 0.21 V in contrast to the calculated result (0.18 V).265 We ascribe the excellent OER activity of sAu/NiFe LDH to the charge redistribution of active Fe and its surrounding atoms caused by the neighboring sAu on NiFe oxyhydroxide stabilized by interfacial CO32− and H2O interfacing with LDH (Fig. 19).265
Fig. 19 (a) CV curves, (b) Overpotential (η) at 10 mA cm−2 (left) and Tafel slope (right) for sAu/NiFe LDH and pure NiFe LDH. (c) Polarization curves of sAu/NiFe LDH before and after 2000 cycles. (d) Raman spectra of sAu/NiFe LDH at different potentials during the anodic and cathodic sweep in a CV cycle. This figure has been adapted/reproduced from ref. 265 with permission from the American Chemical Society. |
Fig. 20 ORR and OER activity. (a) LSV curves for NPMC-900, 1000, 1100, 1000, NPC-1000 and benchmark Pt/C catalyst (1600 rpm.) in 0.1 M KOH O2-saturated solution at a scan rate of 5 mV s−1. (b) LSV curves of NPMC-1000 in 0.1 M KOH at different rotations. (c) NPMC-1000 and Pt/C K–L plots at different potentials. (d) Kinetic current of different samples at 0.65 V for O2 reduction. (e) Different catalyst loading RRDE measurements (1600 rpm.) for the ORR at the NPMC-1000 electrode. (f) LSV curves of NPMC1000, NPMC-1100, RuO2 and Pt/C catalyst (1600 rpm.) in 0.1 M KOH (5 mV s−1 scan rate). This figure has been adapted/reproduced from ref. 261 with permission from Nature Publishing Group. |
Fig. 21 (a) CVs for doped and un-doped graphene nanoporous structure using diverse CVD heat treatments and unusual dopants. (b) HER polarization curves of the different materials in 0.5 M H2SO4, at a scan rate of 5 mV s−1. This figure has been adapted/reproduced from ref. 95 with permission from the Royal Society of Chemistry. |
Here, we first focus on 2D materials, especially the typical 2D graphene independent from metal and metal-based TMDs and explain their properties and applications based on combined DFT computations and experimental observations, especially for the HER.128 2D crystals have been further explored to boost their HER nature, by improving the electron transport efficiency.97 However, the HER mechanism strongly depends on the type of material; thus, it is significant to understand the fundamental properties of electrocatalysts to optimize their use in the HER. To advance the activity, we need to apply strategies to further adjust the electronic structure of 2D materials to afford abundant active sites because loss in resistance is one of the main setbacks in non-metal electrodes. Fig. 22 shows that there are four main way to enhance the HER performance of 2D electrocatalyst materials.
Fig. 22 (a–d) Scheme for obtaining high efficiency 2D materials, (e) enhancing the catalytic efficiencies of MoS2 with S doping and strain, (f) O2 incorporation, (g) 1T0-MoS2, (h) Formation of heterostructures, and (i) summary of the overpotentials @10 mA cm2 with Tafel slopes of 2D HER electrocatalysts. Data points are broken down into two parameters, i.e., Volmer-Tafel route (purple) and Volmer–Heyrovsky route (blue). This figure has been adapted/reproduced from ref. 173 with permission from the Royal Society of Chemistry. |
1. Firstly, create surface defects (e.g., vacancies and edges) (Fig. 22(a)),274 which will cause high i with small coordinated active sites.
2. Secondly, create/introduced strain on the material (Fig. 22(b)), where strain and edge-site substitution with TMs (e.g., Fe, Ni, Cu, and Co) will further adjust the H bonding at the edges.275
3. The third is phase engineering (Fig. 22(c)).
4. Fourthly, interface modulation, by which the electronic structures and the DOS around the Fermi energy of 2D materials can be adjusted (Fig. 22(d)).
5. The fifth is the use of a highly conductive substrate or support coupling with 2D-materials.
In first case of surface defects, MoS2 is used as an example, which was initially known as a weak HER catalyst. Regarding improving the HER activity of MoS2 in acidic media, Jaramillo et al. explained the apparent contradictions160 that the HER normally occurs at the MoS2 edge sites with a less active basal plane. Thus, to get more potential use of the active edges, perpendicularly arranged MoS2 structures were fabricated, which exposed a huge quantity of edge sites, but the basal plane surface sites were not exploited and advance electron transportation to conductive substrates occurred. Using DFT simulation, Li et al.276 suggested that the fundamental MoS2 plane may become active by generating sulfur vacancies in 2D sheets (Fig. 22(e)).
Regarding exerting strain, Li et al.276 showed that the introduction of surface defects under applied strain will further boost the HER activity of MoS2 sheets (Fig. 22(e)). For example, in case of sulfur vacancies under applied strain, the fundamental plane ΔGH reached up to 0 eV and the DOS around the Fermi level increased, causing an enhancement in HER activity. Also, Xie et al.274 explained that the HER activity of 2D MoS2 was noticeably enhanced by simultaneously applying disorder manufacturing and oxygen doping (Fig. 22(f)), where, the disordered structure presented a large quantity of active unsaturated sulfur atoms, and O2 doping successfully modified the electronic structure, improving conductivity. The third strategy is phase engineering, for example, metallic 1T-phase MoS2 exhibited an enhanced HER performance compared to the 2H semiconducting phase. The main reason for this maybe that the electron transfer in metallic 1T phase MoS2 is greater at the fundamental planes than that in the 2H phase equivalent (Fig. 22(g)). The final plan related to interface modulation is to tune the electronic structures of 2D materials and their DOS near the Fermi energy. Li et al.277 observed that the HER performance in 2D MoS2 could further enhanced by synthesizing MoS2/rGO hetero-structures (Fig. 22(h)). Luo et al.278 also proved that the HER activity of MoS2 in alkaline medium can be enhanced by attaching OH nanoclusters in 2D MoS2. Moreover, besides MoS2, other 2D materials and their hetero-structures illustrate good HER activity, such as graphene, b-Mo2C, TaS2, and NbS2. The HER activities of different types of graphene, carbides, TMDCs, and 2D hetero-structures in acid solution are presented in Fig. 22(i).
Besides graphene, nanosheets of (N, S dual-doped) Mo2C,280 GaS,281 and TMDs e.g., MoS2, WS2, WSe2, ReS2,MoTe2, MoSe2, VS2, VSe2, and WTe2, also show HER activity.160,282,283 The layer-dependent electrocatalytic activity of MoS2 for the HER has been observed.284 The decrease in exchange i with a 4.47 factor was observed with an increase in the number of layers. Similarly, TM and chalcogen edges (MoSe2 and WS2) also show potential for the HER, while only the Se and Mo edges are active sites for MoS2 and WSe2, respectively. The W-edge sites of WSe2 and S-edge sites of MoS2, where the H binding is very weak and too strong, respectively, are unsuitable for the HER to occur.282 The experimental285 and theoretical studies based on the comprehensive performance of six TMDCs by Gholamvand et al.286 clearly showed the activity follows the trend selenides > sulphides > tellurides, with MoSe2 outperforming other materials.282,287 Monolayer VS2 is probably more important than all the above 6 TMDCs in catalyzing the HER with a performance analogous to the benchmark Pt,288 @ 10 mA cm−2i with an almost 30 mV more negative potential, and a Tafel slope of 34 mV dec−1 compared to 30 mV dec−1 for Pt.289 The TMDC edge planes are expected to support faster mixed electron transport than basal planes, and hence show higher HER activity. For example, Mo-edges with sulfide-end in 2H-MoS2 (trigonal prismatic phase) are HER active sites, but the basal planes are not.290 However, defects in MoS2 may significantly develop the HER activity of its “inert” in-plane domains.291 Furthermore electron hopping can affect the conductivity of TMDCs, which will also improve their HER properties. The meta-stable octahedral 1 T phase, relative to the trigonal prismatic 2H phase is thermodynamically stable, has intrinsically lower charge transfer resistance with elevated active sites due to strained lattice distortion, and causes a quicker HER.292 In 1T-TMDCs compared to 2H-TMDCs, their catalytic activity is independent of the edge sites.293 Especially, elevated fundamental HER activity is obtained on porous 1T-MoS2 nanosheets, with an overpotential of 153 mV@ 10 mA cm−2, 43 mV dec−1 Tafel slope, and a large electrochemical active SSA (Fig. 23).294
Fig. 23 (a) iR-corrected J–V curves of porous 1T MoS2 (P-1T MoS2), 1T MoS2, porous 2H MoS2 (P-2H MoS2), 2H MoS2, porous 2H MoS2 after S compensation (P-2H MoS2 + S), bulk MoS2 and a Pt wire. (b) Extraction of the double-layer capacitance in MoS2. This figure has been adapted/reproduced from ref. 95 with permission from the Royal Society of Chemistry. |
Metal dopants such as Co, Mn, Cu, Ta, Ru, Fe, Rh and Ni can alter the potency of sulfur binding on edges and modify the hydrogen adsorption free energy to move toward the most favorable ΔGH*.295 However, in HER applications, Pumera et al. illustrated that bulk doped TMDCs (MoS2 and WS2) with p-doped TMDCs (Nb and Ta) showed relatively lower activity than their un-doped counterparts.296 Theoretical computation showed that strain notably altered the free energy of H2 ion adsorption at the interface of defected 1T WS2.297 The ΔGH* value was about 0.28 eV without strain and 0 eV at a strain of 2.7%. Disorder engineering permits the number of low-coordinated atoms to be enhanced, which act as active sites in the HER. The introduction of a few disorders can hinder electron transfer, resulting in a decrease in electrical conductance and catalytic activity, but some suitable disorders are favorites for an improvement in catalytic activity. Oxygen plasma treatment and hydrogen heat treatment are effective and simple experiential strategies to improve the HER activity of MoS2 single sheets.298 In this case, when the MoS2 nanosheets contain 34–40% disorder, then they can replace Pt catalysis, with an η = 160 mV more negative than that of the benchmark 20% Pt/C.274 TMDCs and graphene composites, and also doping in graphene can further boost their conductivity, which will enhance the HER kinetics.299 For example, MoS2/N-GO needed only η = 56 mV @ 10 mA cm−2, which shows good HER activity as an NPM nanosheet catalyst.299 The calculated Tafel slope for 2D/2D hybrid materials was 41.3 mV dec−1, almost similar to benchmark 20% Pt/C (30 mV dec−1). Very thin non-lamellar TMCs (CoSe2,300 NiSe2,301 and Fe–Ni sulfide302) nanosheets, having a metallic nature, are potential NPMCs for the HER. Additionally, doping pyrite CoSe2 with metal atoms, i.e., Mn, will generate more opportunities for active sites for the HER. Furthermore, Mn doping decreases the kinetic energy (K.E) barrier, promoting H–H bond formation among two neighboring adsorbed H atoms and enhancing H2 gas evolution, resulting in comparable catalytic properties for HER with benchmark NPMCs. The reported best TMDC HER catalyst is porous NiSe2 nanosheets, which are obtained from the transformation of LDH precursors,303 with a lower onset by η ∼ 90 mV. The η value required to obtain current densities of 10, 100, and 200 mA cm−2 are 135, 183, and 202 mV, respectively. In addition, well-known 2D HER electrocatalysts contain inexpensive NPMNs, e.g., NiMoNx,304 and NPMHOs, e.g., Ni(OH)2.95
Yanmei Shi et al. demonstrated the synthesis of 1T′-MoS2 crystals and studied different electrochemical microcells with individual MoS2 nanosheets, which revealed that the basal plane of 1T′-MoS2 is highly catalytically active.306 Recently, Yifu Yu et al. studied the phase arrangement, which plays a significant role in the specific formation of inorganic materials (Fig. 24).305 They reported the large-scale formation of micrometre-sized metallic-phase 1T′-MoX2 (X = S and Se)-layered massive crystals with elevated purity. The electrochemical study revealed that the basal plane of 1T′-MoS2 is much more active than that of 2H-MoS2 for the HER in an acid environment, with an onset overpotential of 65 mV and current density of 607 mA cm−2 at 400 mV, which is superior among the MoS2-based catalysts. This outstanding HER activity of the synthesized 2D materials is due to more active nature of their basal plane and improved charge transport ability of 1T′-MoS2 compared to 2H-MoS2.305 Chao Zhang et al.307 experimentally and theoretically investigated the HER performance of Co3S4 porous nanosheets with sulfur vacancies (Co3S4 PNSvac) in alkaline medium, with an onset overpotential of 18 mV and an overpotential of 63 mV for 10 mA cm−2, and large mass activity of 1056.6 A g−1 at an overpotential of 200 mV.307
Fig. 24 (a) Scheme showing the fabrication of 3 types of electrochemical microcells. (b) HER measurements and (c) optical view of the electrochemical set-up. (d and e) Polarization curves obtained and Tafel-slops. This figure has been adapted/reproduced from ref. 305 with permission from Nature Publishing Group. |
The main key to designing materials with improved HER activity is based on understanding how to organize their reactive intermediate binding energies on surfaces.2Fig. 25(a) shows the volcano plot, where each element has a different electrochemical reaction representation and show a wider structure, by which catalysts cross a broad selection of chemical reactions. Besides this volcano plot, further aspects are not present in the straightforward description model, which are necessary to quantitatively confirm the complete reaction rates. For example, deviations in kinetic barrier dimensions from one material class to the next, and therefore MoS2 has lower exchange current densities than PMCs, although it has ΔGH close to a moderate value. The kinetic barrier variations are also a function of pH for a given potential vs. RHE, which shows the pH dependence of i.159 Although feasible alterations in this process are challenging (e.g., regarding kinetic barriers or pH), it is noticeable that the activity volcano does not shift toward the right or left, and somehow only up and down, showing that the descriptor still provides information on the bonding nature for the most selective HER catalysts.308 However, for complete explanations and understanding, additional detailed and capable techniques for measuring electrochemical barriers are necessary for hydrogen ion transfer reactions containing OH and hydronium ions.162,309,310 Since Pt sits near the top of the hydrogen volcano, with an approximately thermo-neutral ΔGH, and is the most active catalyst for the HER, it needs small η values for elevated response times in acidic solutions (Fig. 25(b) and (c)). On the other hand, scarcity and high cost of Pt may decrease its extensive technological utilizations. Thus, the boosted earth-abundant NPMCs can probably replace Pt.153 For decades, MoS2 was considered an inactive catalyst for the HER,311 but motivated by hydrogen-producing enzymes (hydrogenases and nitrogenases), DFT-based simulations were performed to study the Mo(1010) edge of MoS2, showing about 50% hydrogen coverage, and a ΔGH = 0.08 eV, near 0 eV (optimum value) (Fig. 25(a))312 but the basal plane shows poor activity (ΔGH = 1.92 eV) in bulk MoS2 crystals.161 Therefore, due to these calculations, MoS2 was synthesized on a carbon black support to increase its edge sites and it was successfully tested in a membrane electrode setup. The statistical area-normalized @ 10 mA cm−2geoi was attained at η = ∼175 mV, which at that time was the most active reported NPMC for the HER in acidic electrolyte.2 Shortly after, it was experimentally shown that the MoS2 edges are certainly catalytic active sites in the HER.160 An MoS2 monolayer was deposited on the surface of Au(111) nanoparticles, which showed that the HER activity linearly increased with MoS2 perimeter length, but independent of its SSA (Fig. 25(b) and (d)). To explain this, simulations and experimental observations showed that the periphery is more active, thus motivating the development of MoS2 catalysts with considerable bare edge sites. Therefore, the 3D mesoporous MoS2 nanostructure with a double-gyroid morphology was discovered (Fig. 25(d)), which minimized the formation of basal planes, and exposed the maximum active edge sites. Thus, average turnover frequency (TOFavg) increased from 2 to 4 from that of MoO3–MoS2 nanowires synthesized via the same sulfide method (Fig. 25(d)).290,313 A limitation of the double-gyroid structure is the elongated electron transport distance from the active site to conductive substrate, which enhances resistive loss since the electron mobility perpendicular to the MoS2 basal planes is lower than the in-plane charge mobility by almost thrice in magnitude. Thus, synthesized vertically aligned MoS2 nanostructures bare numerous edge sites at the surface and also enable facile electron transfer to the conducting substrate (Fig. 25(d)).314,315 Another smart method in catalyst growth is to disperse high SSA nanoparticles on supports.2 For example, MoS2 nanoparticles were prepared on rGO,277 which resulted in better dispersion and decreased MoS2 aggregation, leadin to superior activity because of the improved edge sites and superior ion transport (Fig. 25(b)). Lithium ion intercalation into the vdW gap of MoS2 has also been examined to increase the HER activity by changing the electronic properties of MoS2.293,316 This showed that MoS2 chemical exfoliation and phase change from 2H semiconducting phase into metallic phase 1T are other ways to tailor the catalytic activity Fig. 25(d).
Fig. 25 (a) HER volcano plot, (b) TOFavg plots with LSV, (c) order development in η of MoS2-based and phosphide HER catalysts, and (d) microscopy images of HER catalysts. This figure has been adapted/reproduced from ref. 2 with permission from the American Association for the Advancement of Science. |
It is suggested that the improved activity of 1T-MoS2 over similarly prepared 2H-MoS2 is due to the increase in active edge sites, and decrease in charge transfer resistance.316 Based on another study, the edge 1T-MoS2 is not the only major active sites, but the basal plane can also be possible catalytically active sites.293 Thus, recently, it was exposed that the vacancies in the MoS2 basal plane also show activity, which can also be modified by strain.276 Similarly, in the case of molybdenum sulfide, its high HER activity was studied, which was shown to be mainly because of its high SSA (Fig. 25(d)). The calculated as-synthesized amorphous MoS2 composition was close to MoS3, but considering its reduction potentials in an electrochemical cell, it was MoS2, which was also supported by in situ studies.317,318 Further doping with TMs (Ni, Fe, and Co), the catalyst performance is considerably improved in acidic medium319 because of the increase in SSA, while the TOFavg is enhanced in neutral environment. Different methods to tailor MoS2-based catalysts actually increase the active sites numbers; however, their overall activity is still limited, which may be due to the small portion of edge sites contributing to the rate of the reaction (Fig. 25(b–d)). This requires us to further modify our strategy and design molecular clusters with under coordinated S at the surface analogous to MoS2 edges, such as [Mo3S4]4+ cubanes, Mo-IV-disulfide, and thiomolybdate [Mo3S13]2− complexes (Fig. 25(d)). Another strategy to increase the number of active sites is proposed based on increasing the catalyst loading, but it results in the restriction of charge and mass transport. This will encourage advancement of further catalysts with higher intrinsic activity through experimental and simulation work on descriptors based on ΔGH ≈ 0. Hence, some TM phosphides, selenides, borides, carbides, and nitrides show HER activities closer to Pt in the form of η @ 10 mA cm−2geo. Although, due to high catalyst loadings and large SSAs, NPMCs lag behind the benchmark Pt regarding the TOFavg in acidic medium (Fig. 25(b–d)). In alkaline media, to date, NPMCs have also developed, e.g., Ni–Mo systems showed low η values @ 10 mA cm−2geo. However, they also exhibit lower TOFavg values than the benchmark Pt, similarly to NPMCs in acidic media. Homogeneous catalysts with high TOFavg have also been developed, even though they classically need a large η value to reach substantial current densities, thus leaving much room for improvement (Fig. 25(b)).2,320Table 3 shows recent studies on ultrathin 2D electrocatalysts for the HER.
2D material electrocatalyst | Electrochemical performance | Ref. |
---|---|---|
2H-phase MoS2 | Onset over-potential: ca. 100 mV, Tafel slope: ca. 50 mV dec−1 | 321 |
Defect- and S-rich MoS2 | Overpotential: 135 mV(@ 10 mA cm−2), Tafel slope: 48 mV dec−1 | 322 |
1T metallic MoS2 | Overpotential: 187 mV(@ 10 mA cm−2), Tafel slope: 43 mV dec−1 | 316 |
Defect-rich MoS2 | Onset overpotential: ca. 120 mV, overpotential: 190 mV (@ 10 mA cm−2), Tafel slope: 50 mV dec−1 | 323 |
MoS2 | Onset over-potential: ca. 100 mV, Tafel slope: 73 mV dec−1 | 324 |
O-Doped MoS2 | Onset over-potential: ca. 120 mV, Tafel slope: 55 mV dec−1 | 274 |
1T-VS2 | Overpotential: 68 mV (@ 10 mA cm−2), Tafel slope: 34 mV dec−1 | 325 |
WS2 | Onset over-potential: ca. 80 mV, overpotential: 240 mV (@ 10 mA cm−2), Tafel slope: 55 mV dec−1 | 297 |
WS2 | Onset over-potential: ca. 100 mV, Tafel slope: 48 mV dec−1 | 326 |
1T-WS2 nanosheets | Overpotential: 142 mV (@ 10 mA cm−2), Tafel slope: 70 mV dec−1 | 327 |
VSe2 | Onset over-potential: ca. 108 mV, overpotential: 206 mV (@ 10 mA cm−2), Tafel slope: 88 mV dec−1 | 328 |
1T-MoSe2 | Onset over-potential: ca. 60 mV, overpotential: 179 mV (@ 10 mA cm−2), Tafel slope: 78 mV dec−1 | 329 |
2H–1T MoSe2 mixture | Overpotential: 152 mV (@ 10 mA cm−2), Tafel slope: 52 mV dec−1 | 330 |
MoS2(1−x)Se2x | Onset over-potential: ca. 80 mV, overpotential: 164 mV (@ 10 mA cm−2), Tafel slope: 48 mV dec−1 | 331 |
Supramolecular polymer | Onset over-potential: ca. 110 mV, overpotential @ 10 mA cm−2: 333 mV, Tafel slope: 80.5 mV dec−1 | 332 |
Pt-MoS2 | Overpotential: 55 mV (@ 10 mA cm−2), Tafel slope: 40 mV dec−1 | 333 |
Single-atom Pt-doped MoS2 | Overpotential: 145 mV (@ 10 mA cm−2), Tafel slope: 96 mV dec−1 | 334 |
Ru@C2N | Onset over-potential: ca. 95 mV, overpotential: 22 mV (@ 10 mA cm−2), Tafel slope: 30 mV dec−1 | 335 |
NiCo2S4 | Onset over-potential: ca. 17 mV, overpotential: 65 mV (@ 10 mA cm−2), Tafel slope: 84.5 mV dec−1 | 336 |
Mn-doped CoSe2 | Onset over-potential: ca. 174 mV, overpotential: 195 mV (@ 10 mA cm−2), Tafel slope: 36 mV dec−1 | 337 |
Ultra-thin Fe–Ni–S | Overpotential @ 10 mA cm−2: 105–117 mV, Tafel slope: 40–48 mV dec−1 | 302 |
NiSe2 | Overpotential: 135 mV (@ 10 mA cm−2), Tafel slope: 37 mV dec−1 | 250 |
CoP | Onset over-potential: ca. 40 mV, overpotential: 90 mV (@ 10 mA cm−2), Tafel slope: 43 mV dec−1 | 338 |
Ni5P4–Ni2P | Onset over-potential: ca. 54 mV, overpotential: 120 mV (@ 10 mA cm−2), Tafel slope: 79.1 mV dec−1 | 339 |
Mo–W–P | Overpotential: 93 mV (@ 20 mA cm−2), Tafel slope: 52 mV dec−1 | 340 |
MoP | Onset over-potential: ca. 50 mV, overpotential: 124 mV (@ 10 mA cm−2), Tafel slope: 58 mV dec−1 | 341 |
CoP | Overpotential: 56 mV (@ 10 mA cm−2), Tafel slope: 44 mV dec−1 | 342 |
C@Ni8P3 | Overpotential: 110 mV (@ 10 mA cm−2), Tafel slope: 46 mV dec−1 | 343 |
Ni0.9Fe0.1PS3 | Overpotential: 72 mV (@ 10 mA cm−2), Tafel slope: 73 mV dec−1 | 344 |
NiMoNx | Onset overpotential: ca. 78 mV, Tafel slope: 35.9 mV dec−1 | 304 |
MoN | Onset over-potential: ca. 100 mV, overpotential: 220 mV (@ 10 mA cm−2), Tafel slope: 90 mV dec−1 | 221 |
N/S co-doped G | Onset potential: 0.13 V (vs. RHE) Operating potential: 0.28 V (@ 10 mA cm−2, vs. RHE) Tafel slope: 80.5 mV dec−1 | 270 |
(Co/N co-doped mesoporous carbon spheres)/N-doped carbon nanosheets | Over-potential: 0.22 V (@ 10 mA cm−2) Tafel slope: 81 mV dec−1 | 345 |
CoS2/CNT/G | Overpotential: 0.142 V (@ 10 mA cm−2) Tafel slope: 51 mV dec−1 Exchange i: 6.26 105 mA cm2 | 346 |
CoSe2 nanosheets | Overpotential: 0.27 V (@ 10 mA cm−2), TOF: 745 h−1, Tafel slope: 64 mV dec−1 | 300 |
CoxW(1−x)S2 nanosheets | Over-potential: 0.121 V (@ 10 mA cm−2) Tafel slope: 67 mV dec−1 | 275 |
Fe–Ni sulfide nanosheets | Overpotential: 0.105 V (@ 10 mA cm−2) Tafel slope: 40.0 mV dec1 | 302 |
g-C3N4/N-doped G | Overpotential: 0.24 V (@ 10 mA cm−2) Exchange i: 3.5 107 A cm−2 | 272 |
Co/N-Doped G | Overpotential: 0.03 V (@ 0.3 mA cm−2) Overpotential: 0.147 V (@ 10 mA cm−2), Tafel slope: 82 mV dec−1 | 273 |
GaS nanosheets | Onset potential: 0.48 V (@ 1 mA cm−2), i @ 0.6 V (vs. RHE): 22 mA cm−2, Tafel slope: 85 mV dec−1 | 281 |
Mn-CoSe2 nanosheets | Overpotential: 0.174 V (@ 10 mA cm−2), Tafel slope: 36 mV dec−1 Exchange i: 0.0683 mA cm−2 | 337 |
N/S co-doped Mo2C nanosheets | Onset potential: 0.046 V (vs. RHE) Operating potential: 0.086 V (@ 10 mA cm−2, vs. RHE), exchange i: 3.8 102 mA cm−2, Tafel slope: 47 mV dec−1 | 347 |
Strained MoS2 nanosheets (S vacancies) | Formal potential: 0.53 V (vs. Ag/AgCl), electron-transfer coefficient: 0.4, electron-transfer rate constant: 2.3 104 cm s−1 | 283 |
Hierarchical MoS2 nanosheets | Over-potential: 0.167 V (@ 10 mAcm−2), Tafel slope: 70 mV dec−1, turnover frequency (TOF):0.41 S1 at an over-potential of 0.125 V | 291 |
Metallic-phase MoS2 nanosheets | i:10 mA cm−2 (@ 0.175 V) Tafel slope: 41 mV dec−1 | 292 |
MoS2/mesoporous G | Onset overpotential: 0.14 V, Tafel slope: 42 mV dec−1 | 348 |
MoS2/N-doped G | Onset potential: 0.005 V (vs. RHE), over-potential: 0.056 V (@ 10 mA cm−2), Tafel slope: 41.3 mV dec−1, exchange i: 0.74 mA cm−2 | 299 |
MoS2/G | Onset overpotential: 0.03 V, i @ over-potential of 0.11 V: 10 mA cm−2 | 349 |
C/MoS2@G | Onset potential: 0.16 V (vs. RHE) Onset overpotential: 0. 165 V, Tafel slope: 46 mV dec1 | 350 |
N-CNT/MoS2+x nanosheets | Onset potential: 0.135 V (vs. RHE), overpotential: 0.08 V (@ 10 mA cm−2) | 351 |
MoS2/CoSe2 nanosheets | Onset potential: 0.011 V (vs. RHE), overpotential: 0.068 V(@ 10 mA cm−2), exchange i: 0.073 mA cm−2 | 352 |
MoSe2 nanosheets | Exchange i: 1.9 × 10−3 mA cm−2, Tafel slope: 76 mV dec−1 | 286 |
Ni nanosheets | Onset potential: 0.034 V (vs. RHE), Tafel slope: 114 mV dec−1 | 267 |
NiMoNx/C nanosheets | Onset potential: 0.078 V (vs. RHE), Tafel slope: 35.9 mV dec−1 Exchange i: 0.24 mA cm−2 | 304 |
NiSe2 nanosheets | Over-potential: 0.117 V (@ 10 mA cm−2), Tafel slope: 32 mV dec−1 | 353 |
Vertically oriented ReS2 nanosheets/Au | Over-potential: 0.2 V (@ 10 mA cm−2), Tafel slope: 85 mV dec−1 | 354 |
Ru nanosheets | Onset potential: 0 V (vs. RHE), overpotential: 0.02 V (@ 10 mA cm−2), Tafel slope: 46 mV dec−1 | 95 |
SnS/N-doped G | Overpotential: 0.125 V (@ 10 mA cm−2), TOF: 0.23 S−1 at an over-potential of 0.125 V, Tafel slope: 38 mV dec−1 | 355 |
VS2 nanosheets | Overpotential: 0.068 V (@ 10 mA cm−2), overpotential: 0.34 V (@ 100 mA cm−2), Tafel slope: 34 mV dec−1 | 289 |
VSe2 monolayers | Overpotential: 0.206 V (@ 10 mA cm2), Tafel slope: 88 mV dec−1 | 328 |
Hierarchical WS2 nanosheets | Overpotential: 0.16 V (@ 10 mA cm−2), Tafel slope: 60 mV dec−1,TOF: 1.24 S−1 at and over-potential of 0.15 V | |
Vertically aligned WS2 nanosheets | Onset overpotential: 0.03 V IR-corrected kinetic i: 48 mA cm−2 at and overpotential of 0.221 V, Tafel slope: 61 mV dec−1 | 356 |
WS2/G | Overpotential: 0.17 V (@ 10 mA cm−2), Tafel slope: 52 mV dec−1 | 357 |
2H-phase MoS2 | Onset overpotential: ca. 100 mV, Tafel slope: ca. 50 mV dec−1 | 321 |
Mo2CTx (rich-O termination group) | Overpotential: 189 V (@ 10 mA cm−2), Tafel slope: 75 mV dec−1 Stability: >1000 CV | 358 |
0.5W-MoxC | Overpotential: 148 V (@ 20 mA cm−2), Tafel slope: 56 mV dec−1, stability: >24 h | 359 |
Mo/Mo2C-HNS-750 | Onset overpotential: ca. 16 mV, overpotential: 89 V (@ 10 mA cm−2), Tafel slope: 70.72 mV dec−1, stability: >20 h | 360 |
N-Doped Mo2C | Onset overpotential: ca. 48.3 mV, overpotential: 99 V (@ 10 mA cm−2), Tafel slope: 44.5 mV dec−1, stability: >12 h | 361 |
Ti2CTx (rich-O termination group) | Onset overpotential: ca. 75 mV, overpotential: 170 V (@ 10 mA cm−2), Tafel slope: 100 mV dec−1, stability: >5 h | 362 |
MoC@2D-NPCs | Onset overpotential: ca. 0 mV, overpotential: 45 V (@ 10 mA cm−2), Tafel slope: 46 mV dec−1, stability: >20 h | 363 |
Mo2CTx (MXene with-O termination) | Overpotential: 189 V (@ 10 mA cm−2), Tafel slope: 70 mV dec−1, stability: >2 h | 364 |
MoS2/Ti3C2-MXene@C | Onset overpotential: ca. 20 mV, overpotential: 135 V (@ 10 mA cm−2), Tafel slope: 45 mV dec−1, stability: >20 h | 365 |
a-Mo2C | Overpotential: 198 V (@ 10 mA cm−2), Tafel slope: 56 mV dec−1, stability: >1000 CV | 366 |
NS-Doped Mo2C | Onset overpotential: ca. 46 mV, overpotential: 86 V (@ 10 mA cm−2), Tafel slope: 47 mV dec−1, stability: >25 h | 367 |
Fig. 26 (a) LSV curves of Ni2P at a scan rate of 5 mV s−1 in 40 mL of 1.0 M KOH with and without 0.5 mmol 1a. (b) Time-dependent evolution of 1a and its semi-dehydrogenation product 2a. (c) FEs and selectivity of 2a at different potentials. (d) Cycle-dependent yields and FEs of 2a. (e) LSV curves over an Ni2P∥Ni2P electrolyzer. (f) Cycle-dependent FEs of an Ni2P∥Ni2P electrolyzer for both H2 and 2a production. This figure has been adapted/reproduced from ref. 369 with permission from John Wiley and Sons Ltd. |
1. The chemical adsorption of CO2 molecules on the cathode catalyst active sites.
2. Transfer of electron and/or migration of protons to split C–O bonds and convert them into C–H bonds.
3. Products rearrangement followed by desorption from the electrode surface and discharge to the electrolyte.
However, the ECR is a complex process, and it is better to continue through proton-coupled multi-electron shift. Nevertheless, the definite full reaction mechanism mainly relies on numerous factors (Table 4).
Electron transfer | Reaction | E° (V vs. SHE) |
---|---|---|
e− | CO2 + e− → CO2− | −1.9 |
2e− | CO2 + 2H+ + 2e− → CO + H2O | −0.53 |
2CO2 + 2H+ + 2e− → HCOOH | −0.61 | |
2CO2 + 2H+ + 2e− → H2C2O4 | −0.913 | |
4e− | CO2 + 4H+ + 4e− → HCHO + H2O | −0.48 |
6e− | CO2 + 6H+ + 6e− → CH3OH + H2O | −0.38 |
8e− | CO2 + 8H+ + 8e− → CH4 + 2H2O | −0.24 |
12e− | 2CO2 + 12H+ + 12e− → C2H4 + 4H2O | −0.349 |
2CO2 + 12H+ + 12e− → C2H5OH + 3H2O | −0.329 | |
14e− | 2CO2 + 14H+ + 14e− → C2H6 + 4H2O | −0.27 |
18e− | 3CO2 + 18H+ + 18e− → C3H7OH + H2O | −0.31 |
As shown in Fig. 27, the chemical features manage the selectivity, for example redox properties, binding affinity, and acid/base character. Structural and material properties can affect the activity, for example, crystallinity, particle size, and surface area. Lastly, the efficiency is identified by optical and electrical properties, i.e., quantum FE, photonic, and morphology. In this perspective, nature gives very active and selective catalysts, but less competence. In contrast, inorganic catalysts, e.g., metals and metal oxides, exhibit elevated efficiency and activity. However, alteration of selectivity at the molecular level is still very complicated. Homogeneous molecular catalysts are capable of achieving good selectivity and efficiency, but are intrinsically inadequate regarding activity. The view of understanding all three features in a single system may only be attained by reticular chemistry. The Nernst equation shows that the ECR mechanism is more thermodynamically favorable for processes having very positive standard redox potentials (E°). However, the ECR rate depends on several factors, e.g., pressure, electrolyte, electrocatalyst, and temperature, despite the position of thermodynamic equilibrium. Catalysts particularly possess a vital character in the ECR kinetic process, e.g., bulk TMs are categorized in three divisions based on their capability of producing BE-definite response transition species. Also, previous studies proposed that electrocatalyst size/structure, morphology, and chemical state have a considerable impact on reaction kinetics. Additionally, the electrolyte has a significant effect on the ECR kinetics. Generally, aqueous solutions having salts of inorganic metals (e.g., KHCO3 and NaHCO3) are used as the electrolyte, and thus the HER (E° = −0.42 V vs. SHE) is not the preferred challenging reaction. Although non-aqueous organic electrolytes can efficiently restrain the competing HER and advance the dissolubility of CO2, their comparatively low viscosity are responsible for poor electrolyte diffusion, and their elevated price and toxic nature restrict their large-scale use.
Fig. 27 Strategy for improving the emerging ECR. This figure has been adapted/reproduced from ref. 370 with permission from the Nature Publishing Group. |
The ECR kinetics is also influenced by temperature and pressure because of the insolubility of CO2 in electrolyte. Usually, a high CO2 partial pressure and low temperature are helpful in increasing the CO2 dissolution. Generally, a high CO2 concentration causes high kinetics, and thus more energy is necessary to generate this circumstance. Thus, ECR studies are mostly performed in CO2-saturated aqueous electrolyte at room temperature. Also, faradaic efficiency (FE), η, and energy efficiency, j, are two additional parameters that are significant for assessing ECR efficiency because of the range of products created. The FE can calculated as follows:
FE = αnF/Q | (37) |
EE = (E°/E° + η) × FE | (38) |
The ECR has potential as a main approach for the renewable manufacture of fuels and chemicals through the use of renewable CO2 free energy sources. However, the development of efficient catalysts is critical since currently, there are no industrial scale operations that use this technology due to its low energetic efficiency. For ECR on a large scale, electrocatalysts must have high activity and selectivity to some product.
In the 1980s, carbon ECR was calculated on many heterogeneous elemental surfaces.371,372 The electrocatalyst results were mainly divided based on their selectivity toward their major reaction product, CO (e.g., Au and Ag), formate (e.g., Pb and Sn), hydrocarbons (e.g., Cu), and hydrogen (e.g., Pt and Ni). Here, we show some non-nanosheet state-of-the-art values and compare them with that of some 2D materials. During the past few decades, many TMs (e.g., first transition series, coinage/precious metals, and few normal metals such as Bi and In), TMOs (TiO2, Cu2O, and SnO2), and TMDCs as well as metal-free carbon materials have been discovered as electrodes for the ECR. It was planned that the starting step in ECR involved CO2 anion radical formation after one-electron transfer in a CO2 molecule, which needs satisfactory energy to rearrange a linear CO2 molecule into a non-linear anion at 1.9 V (vs. SHE), and in most cases this is considered rate-determining. Electrochemical conversion obstacles such as large η (difference between applied potential and equilibrium), sluggish kinetics, unsatisfactory selectivity and low efficiency must be tackled in proposed catalysts, where the recent research on inexpensive electrocatalysts with high activity/selectivity is shown in Table 5.
Catalyst | Synthetic method (precursor) | Electrolyte | Electrochemical performance | Ref. |
---|---|---|---|---|
Bismuth nanosheets | Electrochemical reduction (BiOCl nano-sheets) | 0.5 M KHCO3 | Onset potential: 1.3 V Max FE of formate: 92% at 1.5 V (vs. SCE) | 373 |
N-Doped graphene like material | Carbonation at 1000 °C in Ar (3-pyridine carbonitrile) | 97 wt% [Bmim] BF4, 3 wt% water | Max CH4 FE: 93.5% at 1.4 V (vs. SHE) CH4 partial i (@ 1.4 V, vs. SHE): 3.26 mA cm−2 | 374 |
4-Atomthick Co3O4/Co sheet | Ligand-confined growth (Co(iii) acetyl-acetonate) | 0.1 M Na2SO4 | i (@ 0.85 V): 10.59 mA cm−2, FE for formate:90.1% (@ 0.85 V, vs. SCE) Overpotential @ 10.59 mA cm−2: 0.24 V | 303 |
Mo-Bi bimetallic chalcogenide nanosheets | Solvothermal method (230 °C) (Bi(NO3)3·5H2O and (NH4)2MoS4) dec−1 | [Bmim]BF4/MeCN | Max FE of methanol: 71.2% at 0.7 V (vs. SHE), i at maximum FE of methanol: 12.1 mA cm−2 Tafel slope: 124.4 mV (@ overpotential range 0.16–0.26 V) | 375 |
Vertically aligned MoS2 nanoflakes | Mechanical exfoliation (bulk MoS2) | 96 mol% water, 4 mol% EMIM-BF4 | Onset potential: 0.164 V (vs. RHE), i: 30 mA cm−2 (@ 0.764 V, vs. RHE) | 376 |
WSe2 nanoflakes | CVD | 50 vol% EMIMBF4, 50 vol% deionized water | Onset potential: 0.164 V (overpotential of 54 mV), i: 18.95 mA cm−2(@ 0.164 V), i: 330 mA cm−2 (@ 0.764 V, vs. RHE), CO formation FE: 24% CO formation TOF: 0.28 s−1 | 377 |
SnO2/G | Hydrothermal method (SnCl2, ethylene glycol, and water) | 0.1 M NaHCO3 | i: 13.1 mA cm−2 (@ 1.8 V), specific i: 266 A g−1 FE for formate: 93.6% (@ 1.8 V, vs. SCE) | 378 |
Sn quantum sheets in G | Spatially confined reduction then the hydrothermal process (180 °C), calcinations at 500 °C and then 1000 °C (SnO2 layer and glucose) | 0.1 M NaHCO3 | i: 21.1 mA cm−2 (@ 1.8 V), onset potential: 0.85 V, Max FE for formate @ 1.8 V (vs. SCE): 89%, Tafel slope: 83 mV dec−1 | 379 |
Fig. 28 ECR to formate on four-atom thick partially oxidized Co (red) layers, four-atom thick pristine Co (blue), partially oxidized bulk Co (violet) and intact bulk Co (black). (a) LSVs in CO2-saturated (solid line) and N2-saturated (dashed line) 0.1 M Na2SO4 aqueous solution. (b) FE of formate for 4 h. (c) Chronoamperometry results at the equivalent potentials (in b) with the maximum FE. This figure has been adapted/reproduced from ref. 95 with permission from the Royal Society of Chemistry. |
Fig. 29 N2-Doped graphene catalysis method for CO2 (1) reduction to methane (8). This figure has been adapted/reproduced from ref. 95 with permission from the Royal Society of Chemistry. |
This established procedure involves maintaining the selectivity after the addition of a small quantity of H2O (3%) to the IL, with an increase in i from 1.4 to 3.3 mA cm2, which is about six-fold superior to that of Cu foil. The N2-doped electrocatalyst was significant since only CO and H2 were formed in its absence. This metal-free 2D material shows considerable success in improving metal-based electrocatalysis involving complex multi-electron methods. Many researchers use ILs, which exhibit electrolyte properties, in addition to coordinating and solvating CO2 (provided that pre-concentration occurs).
H+ + * + O2 + e− → OOH* | (39) |
OOH* + H+ + e− → H2O2 + * | (40) |
Thus, it is necessary to design electrocatalysts theoretically with almost zero η, which may have an excellent ΔGOOH and moderate OOH* intermediate BE.150 Although numerous catalysts, e.g., Pt, Ag, Au, Pd–Au alloys, N2-doped C, and hierarchically porous carbon, have been discovered, they were created to show merely reasonable efficiency in manufacturing H2O2. Appropriate catalysts are required to have high selectivity toward the 2e− in contrast to the 4e− pathway. DFT-based simulations were used to draw a volcano framework, which showed the theoretical η value for ΔGOOH for the 2e− reduction in O2 conversion into H2O2 and experimental η value at 1 mA cm−2 are overlaid in this plot (Fig. 30(a)).382 For metals that bind OOH* strongly, the 4e− ORR will follow the 2e− pathway. Alternatively, in the case of weak OOH* binding, the 2e− and 4e− volcano plots overlay on each another, which shows a decrease in activity for H2O2 selectivity with weaker OOH* binding. Thus, the very promising electrocatalysts having elevated activity/selectivity to H2O2 will lie on the top of the 2e− volcano plot. Simulation-based calculated Pt–Hg, Pd–Hg, and Ag–Hg alloys showed not only notable mass activity but also high selectivity (>95%).150,382 The normal design approach has exposed significant ideas to begin showing and recognizing very good electrocatalytic materials for H2O2 formation, especially to avoid the toxicity of Hg. Thus, the thermodynamic framework and expanding it to know the kinetic barriers and interfacial procedure along a wider array of materials and applied circumstances, will assist to give more information for the development of electrocatalysts that can operate at low overpotentials for the selective production of H2O2.2
Fig. 30 (a) Volcano plot for H2O2 production. (b) Volcano plot for ECR. (c) Volcano plot for N2 reduction. This figure has been adapted/reproduced from ref. 2 with permission from the American Association for the Advancement of Science. |
Being anisotropic, TMDCs are similar to graphene in that discernible electron transfer characteristic occurs from their edge and basal planes. Gerischer et al. showed the method of electron relocation at these two orthogonal planes of MoS2 using [Fe(CN)6]3−/4−, Fe3+/2+ and Cu2+/+ redox probes and recognized a connection to the electronic structure of MoS2.407 As one may anticipate, the HET rate of the MoS2 edge plane surpasses that of its basal surface, as also seen in the case of graphene. This fact is recognized due to the considerable overlap between the dxy and dxy22 orbitals of Mo CB and orbitals of the redox probes. In agreement with this, it was determined that the edge plane of macroscopic MoS2 crystals (Fig. 31) demonstrates a fast HET rate and k0obs = 4.96 × 10−5 cm s−1 for [Fe(CN)6]3−/4− and 1.1 × 10−3 cm s−1 for [Ru(NH3)6]3+/2+ redox probes.406 On the other hand, the pristine basal plane showed slow HET rates close to zero for both redox probes. Exfoliation of bulk TMDCs may either enlarge or deteriorate the HET rates towards [Fe(CN)6]3−/4−. In comparison to their individual bulk counterparts, faster HET rates were observed for exfoliated MoSe2 and WS2 across conventional organolithium reagents and aromatic intercalants, but exfoliated MoS2 and WSe2 showed diverse properties, where the rate was quicker in some cases and slower in others relying on the intercalant.408 Sometimes, exfoliation may initiate oxides that obstruct electron movement at the electrode–electrolyte interface. Electrochemical treatment results in structural and electronic alteration, which modify the HET rates of TMDs. Zhang et al. formed exfoliated 2D MoS2 from molybdenite crystals via electrochemical Li-intercalation. Electroreduction of exfoliated MoS2 resulted in enhanced conductivity, with faster HET rates for redox probes than before treatment. On the other hand, oxidation weakens the HET rates of both bulk and exfoliated MoS2. The DFT-based computational study showed that the electron doping during electroreduction is the primary factor stabilizing the 1T phase, which enhances electron transfer, and hence catalytic properties for the HER.396 It has become a general characteristic for TMDCs that electro-reduction enhances their HET rates for [Fe(CN)6]3−/4−, but electrooxidation may impede electron transfer (Fig. 32). Dopants and impurities are known to modify their electron transfer characteristics. N-G, with electron-donating nitrogen dopants set in by thermal exfoliation of graphite oxide in NH3-saturated atmosphere, shows faster HET rates for the [Fe(CN)6]3−/4− redox couple than undoped graphene.409 Also, N-G prepared through N2 plasma treatment of graphene shows superior catalytic activity in the reduction of H2O2 to the untreated graphene.410 It is also reported that metal-based contaminants of Ni, Fe and Co in parts per billion concentrations in graphene also improve its electrocatalytic activity for analytes for example hydrazine and NaHS.411 Also, the occurrence of TM dopants in TMDCs impacts their final HET performance for the [Fe(CN)6]3−/4− redox probe. Nb- and Ta-doped MoS2 showed slightly a lower HET rate than undoped MoS2, whereas improved HET rate upon doping WS2 with Nb or Ta.296 The anisotropic nature of BP also adds to the diverse electron transfer rates on the edge and basal planes. For both the [Fe(CN)6]3−/4− and [Ru(NH3)6]3+/2+ redox probes, quicker electron transfer was observed for the edge but slow electron transfer rates on the basal plane, with poorly defined redox signals for the basal plane. Similarly, the edge plane of BP is very responsive to the oxidation of ascorbic acid analyte, as evident by its larger current compared to the mild current of the basal plane.412 Other members of the pnictogen family show that shear exfoliation enhance their electron transfer nature. There is a noticeable boost in the HET rate upon shear exfoliation of bulk pnictogens, wherein the most accentuated increase is shown in bismuthene contrast with its bulk state. All shear-exfoliated pnictogens exhibit improved catalytic nature towards the oxidation of ascorbic acid. Mainly, antimonene showed a dramatic decrease in potential for the oxidation of ascorbic acid by 0.1 V in contrast with bulk Sb.413 The effect of surface nature on the electron transfer characteristics of h-BN is evident by the different current signals of h-BN immobilized on various carbon-based substrates. Using the [Ru(NH3)2]3+/2+ redox probe, the HET rate decreases with an increase in h-BN mass loading, which denotes slower electrode kinetics occurring on the h-BN surface relative to the underlying carbon-based substrates.44 When h-BN is tailored on smooth substrates, such as a glassy carbon electrode, the cathodic current signal exhibited only a minor increase in intensity. However, there is a considerable rise in cathodic current upon the use of h-BN tailored on a screen-printed electrode, which offered a rough and ridged surface. Recently, the HET rates of MXenes have also been reviewed. Ti3C2, the archetypal MXene, was modified with fluorine and oxygen, resulting in slower electron transfer kinetics for [Fe(CN)6]3−/4– (ref. 414) and sluggish electron transfer on a halogen-terminated diamond electrode because of the feeble interface between the anionic [Fe(CN)6]3−/4− and electronegative fluorine and hydrophobic surface.415 Upon alkalization of Ti3C2, its electronegative fluorine functional groups were replaced with hydroxyl groups, which are less electronegative. Thus, alk-Ti3C2 showed quicker electron movement than before, together with a decline in Rct.414
Fig. 31 Anisotropic effects that influence electron transfer in 2D materials. (a) Schematic illustration of the edge and basal planes of MoS2. Inset: macroscopic molybdenite crystal (approximately 190 mm × 280 mm). Surface morphologies of the (b) edge and (c) basal planes of MoS2 by optical microscopy. This figure has been adapted/reproduced from ref. 406 with permission from John Wiley & Sons, Inc. |
Fig. 32 Surface characteristics that influence electron transfer in 2D materials. (a) Scheme of the reduction of functional groups on GO by thiourea dioxide (TUDO). (b), XPS spectra of graphite oxide after reduction with NaOH in the absence of TUDO for 2 h at 90 °C (NaOH control), and after reduction with NaOH and TUDO for 1 h at 90 °C (chemically modified rGO, T1-CRGO), and their (c) CV recorded in [Fe(CN)6]3−/4− probe. Voltammograms depicting the changes in electron transfer in the treated (d) PdTe2 and (e) PtTe2 towards the [Fe(CN)6]3−/4− probe. This figure has been adapted/reproduced from ref. 175 with permission from John Wiley and Sons Ltd. |
Fig. 33 Mass transport affects 2D electrocatalytic materials. (a) Schematic illustration of engineered channels in MoS2. (b) Catalytic HER presentation of MoS2 nanosheets and MoS2 with engineered channels. This figure has been adapted/reproduced from ref. 175 with permission from John Wiley and Sons Ltd. |
It was also investigated that the basal planes of the metallic 1T-phase of TMDCs, similar to the edges, also have HER active nature.161 Exfoliated 1T-phase group VIB TMDC nanosheets showed better catalytic HER activity to that of the 2H-phase, with the low Tafel slope of 1T-MoS2 of 43 mV dec−1 because the increase in active sites led to higher conductivity in the 1T-phase of TMDCs.316 Besides TMDCs as HER electrocatalysts, MXenes equipped with O* or OH* termination are also promising candidates, although pure MXenes are seldomly used as electrocatalysts because of their low activity. It is exciting to show that the O* basal planes of the majority of MXenes are active for the HER. Especially, the delamination of Mo2CTx enhanced its HER activity despite of having a greater part of exposed basal planes.364 Unexpectedly, the edges of TMDCs are also active for the ECR in the presence of ionic liquids as co-catalysts, regardless of being active sites for the HER, an opposing reaction to the ECR in acidic conditions. In the ECR, CO2 reduction to CO occurs at the Mo-terminated edges of MoS2 due to their metallic nature and high d-electron density.376 The ionic liquid, 1-ethyl-3-methylimidazolium tetrafluoroborate (EMIM-BF4), stabilizes CO2via complex arrangement and prevents the HER. A series of TMDCs nano-flakes including WS2 and WSe2 showing a superior density of edges than their bulk counterparts showed improved ECR catalysis.377 This was the most prominent for the WSe2 nano-flakes terminated with W atoms. Among the 2D materials, layered metal oxides and LDHs are the most exploited OER catalysts recognized due to their wonderful OER nature, where edges account for their OER activity. Liquid-phase exfoliation of bulk LDHs (NiFe, NiCo, and CoCo) into monolayers resulted in up to 4.5 times higher OER catalysis than before, similar to IrO2 catalysts.247 Furthermore, the OER activity was further enhanced upon fracturing single-layered NiFe LDHs into ultrafine nanosheets with a thickness of <3 nm and maximum exposed edges.422 When BP is thinned into nanosheets by liquid exfoliation followed by centrifugation, additional OER active sites are formed. These produced sites are considered to be edges.14 Compared with bulk phosphorus, nanosheets exhibit improved OER performance with an onset potential of 1.45 V and Tafel slope of 88 mV dec−1.14,175
An exciting set of synergistic effects were exposed by graphene co-doped with elements close to carbon in the periodic table, e.g., P and S, besides B and N. N-Doped graphene supplied paramount non-metallic catalysts in the ORR to date, but the impacts of N on the quaternary basal plane against pyridinic N at the edges are still not established and demand additional critical analysis. The mixed metal/graphene composite approach has resulted in electrocatalysts comparable with IrO2, the benchmark OER electrocatalyst. Carbon/nitride MXene (e.g., TiCx carbide) has also been obtained as a better alternative than carbon-supported IrO2. Metal independent P-G, N-G and B/N-G have demonstrated encouraging OER activity. The codoping of graphene with N/S and N/P materials has shown very promising metal independent HER catalysts, approaching the Pt catalyst. The carbon nitrides supported on N-G emerged as comparable with the bulk MoS2 catalyst. In HER catalysis, TMDC nanosheets are being extensively explored, and to date, selenides have been proven to be the most efficient HER electrocatalysts. At the edge sites, the insertion of different metals actually offers bimetallic synergistic effects and enhanced HER rates. TMDCs/rGO combinations also show promise. The ECR is useful for obtaining useful chemicals products (methanol, ethanol, methane etc.). Moreover its reduction is also a key factor to minimize the carbon footprint from the combustion energy route because it is continually dependent on carbon oxidation. Converting CO2 back into fuel gas using 2D-Mxene-based electrocatalysts is in the early stages, and initial subsidiary achievements have been achieved with ultrathin metal and metal sulfide nanosheets. Metal-free carbon/graphene electrodes are also utilized in ionic liquid environment for CO2 reduction. The polymerized ionic liquids show outstanding CO2 adsorption capability, and thus it is probability that their synergism matching with graphene will create a number of thrilling outputs in CO2 reduction. The scalable utilization of 2D materials as electrocatalysts needs some consideration.
1. Understanding the nature of 2D material electrocatalysts through theoretical analysis and in situ characterizations will help in the design and development of customized diverse new nanosheets and atomic level electrocatalysts.
2. The mass production of 2D materials is key factor for their application. There is a need for the controllable production of a number of layers and flakes with dimension confinement for 2D-Mxene materials.
3. The modification of 2D materials in a controllable way, such as doping to regulate their electronic structures and surface engineering to bring in functional groups, is highly recommended for sustainable applications.
4. Manufacturing and functionalization of nanosheet heterostructures are required for tuning electronic properties.
Footnote |
† Equal contribution. |
This journal is © The Royal Society of Chemistry 2019 |