The electrocatalytic carbon dioxide reduction reaction (eCO2RR) driven by renewable energy has great potential for the sustainable production of chemicals and fuels at the gigaton scale that can be used any time, any place1–4. It also offers a promising way to store energy in chemical bonds due to having nearly two orders of magnitude higher energy density compared to the most advanced battery technologies5. However, reducing CO2 to value-added chemicals is both costly and slow based on intrinsic thermodynamics and kinetics, making the goal of an effective and feasible process a real challenge6–9.
Conventional pure metal catalysts such as gold (Au), palladium (Pd), silver (Ag) and newly developed transition metal dichalcogenides (TMDCs)8,10–18 are known to exhibit high activities for the CO2RR in different electrolyte solutions19–26. However, these catalysts are mainly selective for carbon monoxide (CO), known as an intermediate product8,27. Other catalysts such as copper (Cu) and Cu-based catalysts have the ability to reduce CO2 to various chemicals such as methane (CH4), ethylene (C2H4), formic acid (HCOOH), methanol (CH3OH), and ethanol (C2H5OH)28–36. Despite their good selectivity, these catalysts require high potentials –excess energy– to achieve suitable current densities –reaction rates– impeding their use for effective production of chemicals and fuels37,38. Therefore, a novel catalyst with outstanding electronic properties needs to be developed to selectively produce hydrocarbons at high rates at relatively low potentials.
Heteroatomic transition metal carbide (TMC) catalysts, also known as MXenes, have recently received great attention for various electrocatalytic reactions due to their unique structural and electronic properties39–42. In particular, M2C (M denotes transition metals) stoichiometry of this class of two-dimensional materials forms layered structures of M-C-M where a plane of carbon atoms is sandwiched between two hexagonal planes of metal atoms. This structure provides a high density of active metal atoms at the surface breaking conventional scaling relationships that limit electrocatalytic performance of their counterparts such as TMDCs and pure metals43. However, there is limited knowledge of their performance and characteristics as eCO2RR catalysts under actual experimental conditions.
Here, we have studied the performance of di-tungsten carbide (W2C), di-molybdenum carbide (Mo2C) di-niobium carbide (Nb2C) and di-vanadium carbide (V2C) nanoflakes (NFs) as inexpensive, non-precious members of TMCs for eCO2RR.
The TMC NFs i.e., W2C, Mo2C, Nb2C and V2C were synthesized using a carburization process followed by the liquid exfoliation technique (Supplementary section S1)27,44,45. The electrocatalytic performance of TMC NFs with the similar crystallite sizes (25.4±5 nm) were then studied in a three-electrode cell and compared with Au and Cu nanoparticles (NPs), conventional catalysts for this reaction,46 under identical experimental conditions (Supplementary section S2). To improve the CO2RR performance in competing with hydrogen evolution reaction (HER), we have employed a mixture of 3 M potassium hydroxide (KOH) and 2 M choline-chloride (CC) solution (KOH:CC (3M:2M)) as the electrolyte in this study47.
The linear sweep voltammetry (LSV) experiments and a real time product stream analysis show that CO2RR on the W2C surface starts at a potential of -122.7 mV vs reversible hydrogen electrode (RHE) by producing CO and H2 and reach maximum CO2RR current density (jCO2RR) of -548.9 mA/cm2 at -1.05 V vs RHE (Supplementary Figs. S2-4 and Fig. 1a). As shown in Fig. 1a, jCO2RR of -419.9, -381.9 and -350.8 mA/cm2 were observed for Mo2C, Nb2C, and V2C NFs, respectively, at this potential (Supplementary section S3). However, Au and Cu NPs exhibit a jCO2RR of -208.11 and -89.53 mA/cm2 at -1.05 V vs RHE (Fig. 1a). The selectivity analysis also indicates that TMC NFs produce hydrocarbons (i.e., CH4, C2H4, CH3OH and C2H5OH) at a potential range of -0.45 to -1.05 V vs RHE for W2C, Mo2C and Nb2C NFs and a potential range of -0.55 to -1.05 V vs RHE for V2C NFs where CH4 is identified as the main product (Supplementary section S3).
Figure 1b illustrates CH4 formation current densities (jCH4, mA/cm2) of the TMC NFs compared to Cu NPs, a conventional catalyst for hydrocarbon production. The partial current densities of different products (i.e., H2, CO, CH4, C2H4, CH3OH and C2H5OH) were calculated by multiplying FEs and total current densities at different potentials (Supplementary section S3 and Fig. S3). As shown in Fig. 1b, a maximum jCH4 of -421.63 mA/cm2 is obtained for W2C NFs at a potential of -1.05 V vs RHE where Nb2C NFs, Mo2C NFs, and V2C NFs show values of -219.16, -211.33, and -147.56 mA/cm2, respectively, at this potential. We also compared the CH4 formation activity of TMCs i.e., W2C, Nb2C, Mo2C, and V2C NFs with state-of-the-art catalysts in the literature (Supplementary Table S2)46,48–54. Supplementary Table S2b indicates that the jCH4 of W2C NFs is 3.6 and 4.2 times higher than recently studied La2CuO4 (-117 mA/cm2 at -1.4 V vs RHE)51 and Cu-N (-100 mA/cm2 at -1.0 V vs RHE)48, respectively. The partial current densities of other hydrocarbon products i.e., C2H4, CH3OH and C2H5OH are also shown in Fig. S3 (Supplementary section S3).
To evaluate the intrinsic activity of W2C NFs, we measured CH4 formation turnover frequency (TOFCH4) by normalizing its activity to the number of active atoms at the surface using roughness factor method and compared it with the other catalysts in this study (Supplementary section S5). Our calculations indicate a TOFCH4 of 10.42 s-1 at a potential of -1.05 V vs RHE for W2C NFs; by comparison, TOFCH4 of 4.54, 3.74 and 2.79 s-1 were calculated for Mo2C NFs, Nb2C NFs, and V2C NFs, respectively. The calculated TOFCH4 of W2C NFs at the potential of -1.05 V vs RHE is about 2 orders of magnitude higher than that of Cu NPs (0.0736 s-1) under identical experimental conditions (Supplementary Fig. S7a). Moreover, total CO2RR turnover frequencies (TOFCO2RR) of 19.09, 19.36, 17.82, 17.55 s-1 were calculated for W2C NFs, Mo2C NFs, Nb2C NFs, and V2C NFs, respectively, where Au NPs and Cu NPs exhibit TOFCO2RR of 4.35 and 0.1956 s-1, respectively (Supplementary Fig. S7f). These results suggest the superior CH4 selectivity of TMC catalysts compared to state-of-the-art catalysts48–51,54–57.
Furthermore, we performed a comparative mechanistic study by calculating Tafel slopes for different products to gain insight about the eCO2RR mechanism of the TMCs i.e., W2C, Mo2C, Nb2C, and V2C NFs in the two-compartment three-electrode electrochemical cell (Supplementary section S6, Fig. S8) 58. Our Tafel plot analyses show that the TMC NFs possess steeper Tafel slopes, and therefore a weaker potential dependence compared with Cu NPs for the formed products (i.e., CO, CH4 and C2H4) (Supplementary Fig. S8)58. The Tafel plot analyses suggest a different CO2RR mechanism for TMC NFs than that of Cu catalysts where C-O bond scission is the rate determining step58.
To gain more insight to the remarkable performance of these catalysts for electrocatalytic CO2RR, the structural and physicochemical properties of TMC NFs were characterized at molecular and atomic scales by performing X-ray diffraction (XRD), X-ray photoelectron spectroscopy (XPS), and scanning transmission electron microscopy (STEM) (Supplementary sections S7-9). At first, we have performed XPS experiments to analyze the surface chemistry of TMC NFs. XPS analysis (Supplementary Fig. S9) indicates that our NF samples contain metallic TMCs, with little or no evident surface oxidation. The results show that the chemical composition of the surface, the empirical formula of M2C (M: transition metal, C: Carbide) and the oxidation state of +2 for the transition metals i.e., W, Mo, Nb, and V are similar in all synthesized TMCs (Supplementary section S7). The lattice structure and crystallite size of the TMC NFs were then studied by performing XRD experiments. The XRD pattern of W2C NFs shows a sharp peak at 39.91º along with three pronounced peaks at 34.84º, 38.54º, and 52.65º corresponding to (101), (100), (002) and (102) crystal surfaces of W2C, respectively. The XRD spectra of the TMCs show all Bragg peaks of W2C, Mo2C, Nb2C, and V2C NFs; verifying their homogenous and pure structures. The XRD results indicate a constant dominant lattice plane of (101) and a similar average crystallite size of 25.4±5 nm for all synthesized TMCs (Supplementary Fig. S10) 59–61.
Furthermore, we performed atomic scale STEM experiments to study surface atom coordination, crystallite sizes and dominant plane structures of TMC NFs (Supplementary Fig. S11-18). Figures 2a-d show STEM results of W2C NFs. Figures 2a and b indicate high-angle annular dark-field (HAADF) image and corresponding fast Fourier transforms (FFT) of W2C NFs in the <101> zone axis. The atomic models of the <101> zone axis and bright field (BF) image of W2C NFs are represented in Figs. 2c and d, respectively. Figure 2d indicates the carbon atomic columns in the red box and the intensity profile across the red box region showing that the distance between two carbon atoms is 2.55 Å. The STEM results of other TMCs i.e., Mo2C, Nb2C, and V2C NFs are explained in Supplementary section S9. The STEM and XRD results of synthesized TMC NFs confirm that the structure of these materials is a perfect match with standard 1T structure, suggesting a tetragonal symmetry and octahedral coordination of the atoms (Fig. 2e) 43,62. Figure 2e indicates the schematic of 1T structure TMC NFs showing tetragonal symmetry, one layer per repeat unit with octahedral coordination. The lattice constant a is in the range of 3.07 to 3.15 Å for synthesized TMC NFs. The stacking index b indicates the interlayer spacing which is in the range of 4.53 to 5 Å for studied TMCs. As shown in Fig. 2e, 1T atomic coordination provides metal-terminated surface atoms that are known to be favorable binding sites of adsorbed intermediates in eCO2RR 43,62. Our atomic and molecular scale structural analyses indicate that the synthesized TMC NFs have fairly similar structural properties e.g., (1T) crystalline structure with a dominant plane of (101), crystallite sizes, and atomic coordination.
To further discern the difference between the observed electrocatalytic performance of the TMCs, we have studied their electronic properties by performing electrochemical impedance spectroscopy (EIS) (Supplementary section S10) and, work function measurements using ultraviolet photoelectron spectroscopy (UPS) (Supplementary section S11)8. At first, we have employed the EIS experiments to compare the overall electron-transfer properties of the TMC catalysts in the double layer region (Supplementary section S10). To do this, TMC NFs with similar structural and physical properties e.g., sizes, shapes, and mass loadings (0.1 mg/cm2) coated on glass carbon were used as the working electrodes. This results in similar roughness, morphology, intrinsic capacitance, and exposed surface area of the studied samples confirmed by our characterization results (Supplementary sections S5-9). The EIS experiments have been performed at a potential of -310 mV vs RHE for all TMCs under identical experimental conditions (Supplementary section S10). Figure 2f shows the fitted EIS spectra of each TMC catalyst using Randles circuit model, indicating a smaller charge transfer resistance (Rct) for W2C NFs (~17 ohms) compared to the other TMCs, i.e., Mo2C NFs (~25 ohm), Nb2C NFs (~33 ohm), and V2C NFs (~38 ohm)63. The UPS method also was used to compare the surface work function of TMCs (Fig. 2g). The results indicate a lower work function for W2C NFs (0.2 to 0.84 eV) compared to Mo2C NFs (3.92 eV), Nb2C NFs (4.44 eV), and V2C NFs (4.55 eV). The charge transfer resistance obtained by EIS experiments and the surface work function value measured by UPS experiments suggest the superior activity of W2C NFs compared to other TMCs in this study i.e., Mo2C, Nb2C, and V2C NFs.
In addition to our experimental observations, we have performed density functional theory (DFT) calculations to gain more insight to the electronic and catalytic properties of M2C compounds. The aim is to address the enhanced activity and selectivity of these TMCs and to explore both electrochemical (i.e., driven) and chemical (i.e., favorable or spontaneous) processes that distinguish them from other catalysts, such as Au and Cu.
With respect to activity, the electronic density of states (DOS) indicate that transition metal d states dominate at the Fermi level of these TMCs, much more so than elemental Au, another high activity catalyst. Bader charge calculations indicate that metal atoms at TMC surface are significantly more reduced compared to the bulk atoms (Fig. 2h, Supplementary Fig. S24). These results indicate the increased availability of electrons at metal-rich TMC surfaces, which may increase the catalytic activity of TMC NFs.
With respect to the increased selectivity of TMC NFs, especially for CH4 production, we have explored the CO2RR pathway on the W2C (101) surface in detail by using DFT calculations. Focusing initially on electrochemical processes, we employed the computational hydrogen electrode (CHE) model64–66 (Supplementary Tables S7 and S8) to explore the stepwise electronic reduction and protonation of adsorbed species in the low molecular coverage limit. The lowest free energy pathway to produce CH4 with only electrochemical steps is shown in Fig. 3 (see a more detailed path in Supplementary Fig. S25), indicating limiting steps (at zero potential with respect to RHE) that involve protonation of adsorbed CO2 and O, and, most importantly, the desorption of H2O following protonation of adsorbed OH. The difficulty of this final step is not surprising, as W2C (101) strongly adsorbs and spontaneously dissociates water ( , ) without electrochemical assistance. Similarly, our calculations indicate that W2C (101) strongly chemisorbs CO2 ( , bond length ) in contrast to normally weak physisorption on Cu ( ) and other catalyst surfaces28,67–69. Furthermore, the (101) surface of W2C enables favorable and unassisted dissociation of adsorbed CO2 ( , , Supplementary Table S9) suggesting that C–O bond scission may take place in the early stages of CO2 reduction, skipping the uphill production of adsorbed carboxyl. Based on these findings, we propose that W2C (101) distinguishes itself as a catalyst due to an interplay between surface-assisted chemical steps, whose energetics will depend on the local chemical equilibrium at the surface and electrochemical steps that reduce pre-existing surface reagents and open up new pathways for the overall reaction to proceed. More detailed studies of such cooperative catalytic processes and their limiting steps may be encouraged based on the promise of W2C as a high-performance CO2 reduction catalyst. Here, we highlight the plausible cooperative effects of these steps, which set apart W2C from conventional noble metal catalysts and the other TMCs, specially for CH4 production. The immediate benefit of the favorable chemical processes mentioned above should be a higher surface coverage of CO2 (and consequently CO) and an excess of surface protons. This may explain the high Faradaic efficiencies for production of both H2 and CO at low potentials (see Supplementary section S3 and Table S1). However, once a limiting potential (-0.74 V estimate) is reached, the readily protonated products of adsorbed CO that produce CH4 are no longer hindered by a build-up of adsorbed byproducts (O* then OH*), which can now be protonated and released from the surface.
We can divide the complex, multistep reaction into two key parts: initial conversion of adsorbed CO2 to adsorbed CO, followed by conversion of adsorbed CO to CH4 with the release of H2O (see Fig. 3). As indicated in Fig. 3, the first part, generation of adsorbed CO, can be achieved by chemical or electrochemical means. We have direct and favorable chemical conversion of adsorbed CO2 to adsorbed CO and O on W2C (101), but also two electrochemical pathways: production of HO-CO* in a single step (+0.66 eV, Fig. 3) or an alternative initially favorable protonation to OCHO* followed by two uphill electrochemical steps producing first OCH2O* followed by release of H2 and the final product of HO-CO* with a similar free energy cost (+0.68 eV, Supplementary Fig. S25 and Table S9). A final electrochemically driven protonation of HO-CO* favorably releases H2O and leaves CO*.
With chemically or electrochemically generated adsorbed CO, we can proceed to the second part of the overall reaction to produce CH4 from CO2, which involves multiple favorable protonation steps. The W2C catalyst distinguishes itself here. The electrochemical activation of remains thermodynamically favorable (ΔG = -0.26 eV) on W2C (101), whereas on other catalysts, such as Cu, this process is usually uphill with the potential ranging from -0.74 V to -0.97 V vs RHE66,70. Moreover, due to the spontaneous water dissociation, the direct H* transfer step on W2C could be even more favorable with a resultant ΔG= -0.433 eV (Supplementary Table S9). The next two electrochemical steps are downhill (ΔG = -0.04 and -0.58 eV): the first forming the unstable methoxy radical with oxygen attached to a surface W atom; the second leading to spontaneous dissociation into the methyl radical and a surface oxygen atom . The electrochemical conversion of the surface into CH4 is favorable (ΔG = -0.43 eV) and the protonation of the surface oxygen is only slightly uphill (ΔG = +0.03 eV). As we already stated, for the overall reaction on W2C (101) it is the final protonation of to release H2O that is the limiting step (ΔG = +0.74 eV).
We also compared W2C with the other TMCs studied by calculating the energies of adsorption of water and CO2 as well as the potentials of the rate determining step (i.e., protonation of ) for Nb2C, Mo2C, and V2C (Supplementary Table S10). Our calculations indicate that these TMCs also strongly chemisorb CO2 with adsorption energies of -1.32, -1.62, and -0.96 eV, respectively. Moreover, Nb2C also shows favorable C–O bond scission of adsorbed CO2. Additionally, Nb2C, Mo2C, and V2C strongly adsorb water with the energies of -1.87, -1.23 and -0.59 eV, respectively, where Nb2C is the only other catalyst that dissociates water. In contrast to W2C, the energies required for the protonation of are higher: +1.17, +1.25, +0.85 eV for Nb2C, Mo2C, and V2C, respectively (Supplementary Table S10). Therefore, we can conclude that, within this set of four TMCs, W2C possesses the optimal characteristics for efficient completion of CO2RR: (1) sufficiently strong adsorption of CO2, (2) spontaneous dissociation of water, and (3) the lowest limiting potential for OH* protonation. We conclude that the performance of Nb2C is reduced due to its stronger water adsorption, resulting in the protonation of requiring more energy. We would expect Mo2C to have a lower surface coverage of protons and higher costs for the protonation of . The weakest CO2 adsorption on V2C decreases its surface coverage, making it the worst TMC catalyst here, despite its relatively small limiting reaction potential of protonation of OH*.
As we mentioned before, for W2C the realistic network of pathways towards CH4 consists of a potential dependent combination of competing chemical and electrochemical steps with the actual limiting potential being in the range from -0.483 V to -0.744 V vs RHE (see full path. Supplementary Fig. S25), which is consistent with our three-electrode electrochemical experimental results (Supplementary section S3 and Table S1). A steeper Tafel slope for CH4 formation on W2C than other TMCs and Cu, (Supplementary Fig. S8) also indicates the competition between reactions for the active sites on the catalyst surface. Specifically, the spontaneous water dissociation on W2C (101) explains the ease of the HER in our non-acidic electrolyte where the source of protons is normally water. A weak potential dependence of the partial CO current and its small overpotential also originate from the interplay between chemical and electrochemical steps (see Supplementary Material for details).
Experimentally, we have studied the effect of choline chloride (CC) on the activity and selectivity of the TMC catalysts. To do this, we have performed electrochemical CO2RR in different choline chloride concentrations of i.e., 0.01, 0.1, 1, and 2 M mixed with 3 M KOH (Supplementary section S13). Figure 4 shows CO2RR overall current density and different products (i.e., CH4, C2H4, CO, alcohols-CH3OH and C2H5OH- and H2) partial current densities for W2C NFs in different CC concentration electrolytes. Figure 4a indicates that by increasing the concentration of CC in the electrolyte the CO2RR current density (jCO2RR) increases and reaches a maximum value of -548.89 mA/cm2 at a potential of -1.05 V vs RHE for 2 M of CC. The obtained value is about 32, 24, 17, 9% higher than that of 0, 0.01, 0.1, and 1 M of CC, respectively. Moreover, a maximum CH4 formation current density (jCH4) of -421.63 mA/cm2 is obtained for 2M CC at a potential of -1.05 V vs RHE that is about 1.41, 1.29, 1.19 and 1.1 times higher than that of 0, 0.01, 0.1, and 1 M, respectively (Fig. 4b).
The results also indicate using W2C NFs, maximum partial current densities of other products i.e., C2H4 (jC2H4 of -35.84 mA/cm2), CO (jCO of -78.48 mA/cm2), and alcohols (jAlcohols of -12.81 mA/cm2; -6.84 mA/cm2 for CH3OH and -5.97 mA/cm2 for C2H5OH) were obtained at the potential of -1.05 V vs RHE in 2M CC (Figs. 4b-d). In contrast, the measured H2 partial current densities indicate that by adding higher concentration of CC to the electrolyte solution the rate of H2 production decreases significantly where a minimum H2 formation current density of -4.48 mA/cm2 was obtained for 2 M CC at a potential of -0.85 V vs RHE that is 12.31, 8.97, 6.76, 3.23 times lower than that of 0, 0.01, 0.1 and 1 M CC, respectively.
These results suggest that adding choline chloride to the 3 M KOH electrolyte suppresses the competing hydrogen evolution reaction (HER) and increases the formation of CO2RR products more specifically CH4 27.
The stability of the choline chloride electrolytes was studied by conducting nuclear magnetic resonance (NMR) and 13CO2 isotope experiments (Supplementary section S14 and S15) 27,46,71. The 1H and 13C NMR spectra reveal similar peak areas and chemical shifts for fresh and used electrolytes indicating no generation of new diamagnetic species or change in the choline chloride structure under applied potential of -1.05 V vs RHE (Supplementary Figs. S36 and S37). The 13CO2 isotope experiments also show that the CO2 gas present inside the electrolyte is the only source of the formed products in the electrochemical CO2RR (Supplementary Fig. S39). These results confirm that choline chloride with different concentrations i.e., 0, 01, 0.1, 1 and 2 M remains stable at the range of applied potentials in the electrochemical CO2RR experiments.
Next, we studied the performance of W2C NFs in our developed solid polymer electrolyte flow electrolyzer for continuous electrochemical CO2RR using this catalyst as the cathode (Supplementary section S16). The flow electrolyzer used in this study consists of a two-compartment electrochemical setup with active area of 5 cm2 coated with W2C NFs at the cathode and iridium oxide nanoparticles (IrO2 NPs) as the anode and were then fed with humidified CO2 and KOH:CC (3M:2M) electrolyte, respectively (Supplementary section S16).
To study the CO2RR performance of W2C NFs in the flow electrolyzer, we performed chronoamperometry (CA) experiments at different cell potentials ranging from -1.5 to -2.3 V for W2C NFs (Supplementary section S16). As shown in Fig. 5a, the results show that at a cell potential of -1.5 V, hydrogen (H2, FE of 54.9%±1.4) and CO (FE of 40.1%±1.8) are the dominant products. However, our measurements indicate that by increasing the cell potential a system becomes more selective for CH4 formation with the maximum FE of 82.7%±2 at a cell potential of -2.1 V. At this potential, W2C NFs slightly produce other products such as C2H4, CH3OH, C2H5OH, CO, and H2 with FEs of 5.6%, 1.4%, 1.2%, 6.1%, and 1.4%, respectively. Figure 5b shows the maximum CH4, C2H4, CH3OH, and C2H5OH current densities of -421.28, -27.31, -5.95, and -5.19 mA/cm2 at the cell potential of -2.3 V, respectively, confirming high selectivity of W2C NFs towards CH4 as the main product.
Next, we coupled the electrolyzer to a triple junction photovoltaic (TJ-PV) cell with a maximum efficiency of 34.3% to determine the CO2RR performance and energy efficiency of W2C NFs in a solar-driven device (Supplementary section S17). The j-V characteristic curve of the TJ-PV cell under one sun illumination (100 mW/cm2) using a sun simulator light source is shown in Supplementary Fig. S46. The operating point is chosen to provide a photo current density of -394.3 mA/cm2 at a potential of -2.1 V which has the maximum FE of CH4 (82.7%±2) calculated in the flow electrolyzer (Supplementary Fig. S47).
Figure 5c shows the current density of the solar-driven electrolyzer for a 700-hour continuous process at a potential of -2.1 V. The results shown in Fig. 5c indicate a negligible decrease (approximately 2%) in the photo-current density of W2C NFs over the 700-hour experiment while the corresponding photo-potential fluctuates between -2.08 to -2.12 V, confirming the high stability of W2C NFs for CO2RR.
The measured sun to CO2RR products (CO, CH4, C2H4, CH3OH, and C2H5OH) as well as total solar-to-fuel efficiency (SFE) of W2C NFs over a 700-hour process are shown in Fig. 5d (Supplementary section S17). As shown in this figure, an average sun to CH4 production efficiency of 17.3% with negligible variation (2%) is achieved during the 700-hour continuous process. Considering other products, W2C NFs show a SFE of 20.7%.
We also calculated the energy efficiency of CO2RR in our developed flow electrolyzer and compared it with state-of-the-art catalytic systems in the literature (Supplementary section S16)29,30,37,57,72. As shown in this figure (Supplementary Fig. S45), a maximum energy efficiency of 62.3% was obtained for our developed flow electrolyzer using W2C catalyst that is about 67% and 73% more efficient than Cuoh (37.4%)50 and recently developed Cu-CIPH (36.1%)72 catalytic systems, respectively.
In summary, we have synthesized four members of TMCs with a formula of M2C, i.e., W2C, Mo2C, Nb2C, and V2C NFs using the carburization method followed by the liquid exfoliation technique and tested their catalytic performance for eCO2RR in KOH:CC (3M:2M) electrolyte. The electrocatalytic performance studies of TMCs shows these materials are mainly selective for CH4 formation with W2C NFs having the best CO2RR activity compared to the studied catalysts. For instance, a CO2RR current density of -548.89 mA/cm2 and a maximum CH4 current density of -421.63 mA/cm2 at the potential of -1.05 V vs RHE were observed for W2C NFs. Our electrochemical results also indicate that adding choline chloride to the electrolyte enhances the formation of CO2RR products by suppressing the HER for all studied TMCs. Moreover, the NMR and 13CO2 isotope experiments confirm that the choline chloride remains stable during the electrochemical experiments. Atomic and molecular scale characterizations such as XPS, XRD, and STEM indicate that all synthesized TMCs have similar lattice structure of 1T with a dominant plane of (101) and almost the same average crystallite size of 25.4 nm. Furthermore, the electronic property analyses of TMCs reveal superior electronic properties of W2C NFs: low work function; small charge transfer resistance in the electrochemical double layer region; and heavily reduced tungsten atoms at the surface, which may lead to the observed outstanding activity. Computational results also indicate that the studied TMCs spontaneously chemisorb CO2 and water as compared to Cu. However, among the TMCs studied, W2C exhibits the optimal combination for CH4 production, with favorable adsorption energies of water and CO2 coupled with spontaneous dissociation, and less costly protonation of OH*, which is the limiting step, with a low limiting potential in the range of -0.483 V to -0.744 V vs RHE. Using W2C NFs, we have demonstrated a solar-driven flow electrolyzer that can work up to 700 hours with a solar to CH4 efficiency and a total SFE of 17.3 and 20.7%, respectively, under one sun illumination. The demonstrated solar-driven flow electrolyzer using a non-precious metal catalyst (W2C NFs) in this study achieves a maximum efficiency of 62.3% making it a good candidate to approach the commercially relevant electrocatalytic CO2RR. This opens a new direction toward a low cost, sustainable large-scale production of fuels from CO2 that can be used any time any place.