Ultrafast self | Guangzhou Titanium Co.,Ltd

Nature Communications volume 13, Article number: 3338 (2022) Cite this article

7589 Accesses

30 Citations

13 Altmetric

Metrics details

Designing cost-effective and high-efficiency catalysts to electrolyze water is an effective way of producing hydrogen. Practical applications require highly active and stable hydrogen evolution reaction catalysts working at high current densities (≥1000 mA cm−2). However, it is challenging to simultaneously enhance the catalytic activity and interface stability of these catalysts. Herein, we report a rapid, energy-saving, and self-heating method to synthesize high-efficiency Mo2C/MoC/carbon nanotube hydrogen evolution reaction catalysts by ultrafast heating and cooling. The experiments and density functional theory calculations reveal that numerous Mo2C/MoC hetero-interfaces offer abundant active sites with a moderate hydrogen adsorption free energy ΔGH* (0.02 eV), and strong chemical bonding between the Mo2C/MoC catalysts and carbon nanotube heater/electrode significantly enhances the mechanical stability owing to instantaneous high temperature. As a result, the Mo2C/MoC/carbon nanotube catalyst achieves low overpotentials of 233 and 255 mV at 1000 and 1500 mA cm−2 in 1 M KOH, respectively, and the overpotential shows only a slight change after working at 1000 mA cm−2 for 14 days, suggesting the excellent activity and stability of the high-current-density hydrogen evolution reaction catalyst. The promising activity, excellent stability, and high productivity of our catalyst can fulfil the demands of hydrogen production in various applications.

Hydrogen evolution reaction (HER), which involves overall water splitting, has been considered as a promising approach for hydrogen production1. Catalytic activity and stability are important criteria for high-efficiency HER catalysts. In industrial applications, an HER catalyst must be highly active and stable on an electrode at high current densities (e.g., ≥ 500 or 1000 mA cm−2) over a long period of time (≥300 h)2. Pt is usually used for HER but limited for the high cost and scarcity, and thus Pt-group metal-free catalysts with the HER activity comparable to Pt have been extensively studied. However, the development of high-efficiency and Pt-group metal-free HER catalytic electrodes for high-current-density HER is challenging because it requires simultaneous high chemical activity, high chemical stability, and high mechanical stability of the electrodes. In recent years, high activity of HER catalysts at high current densities has been reported via heterogeneous atom doping3,4,5, composite synergy6,7,8,9,10, morphology engineering11,12, and strain engineering13. Nevertheless, many HER catalysts with high activity usually exhibit low chemical stability, and large mechanical forces suffered by active sites during the release of large quantities of H2 bubbles would continually exfoliate the catalyst from the electrodes, making their mechanical stability difficult to maintain at high current densities14,15,16.

It is difficult to simultaneously improve the chemical and mechanical stabilities without affecting the activity of an HER catalyst. To improve the chemical stability, chemically stable Pt-group metal-free catalysts, such as Mo2C17, MoS218,19, MoS2/Ni3S29, and Cr1–xMoxB220 were explored. To enhance the mechanical stability, robust catalytic electrodes were directly employed21, or binders were used to strengthen the adhesion between the catalyst and electrode22; moreover, specific channels or superaerophobic structures of electrodes were designed to reduce the mechanical forces generated in catalysts by facilitating the release of bubbles10,12. However, these approaches to enhancing the stability may weaken the chemical activity of the HER catalysts. For example, additional binders obstruct the exposure of active sites and reduce the overall activity14,23,24. Therefore, development of new methods for intrinsically enhancing the stability of HER catalysts and simultaneously maintaining their high activity is crucial for applications at high current densities.

Instantaneously creating chemical bonding between the active catalysts and electrodes may be an effective way to intrinsically improve the mechanical stability of HER catalysts and not influence their chemical activity and stability. In this regard, fast self-heating (Joule-heating) of a conductive matrix could be used to in situ synthesize chemically bonded catalysts on the matrix and avoid the decay of their activity caused by agglomeration, which is superior to traditional methods25,26,27. Herein, we develop a low-energy-consumption method using a carbon nanotube (CNT) film as a heat source and matrix, which rapidly changes its temperature in hundreds of milliseconds to in situ synthesize a robust Mo2C/MoC/CNT composite film in the presence of Mo and C precursors. The as-prepared uniformly dispersed Mo2C/MoC heterogeneous nanoparticles are tens of nanometres in size and form strong chemical bonds with the CNT film. Consequently, massive Mo2C/MoC interfaces offer abundant active sites for HER, resulting in the Mo2C/MoC/CNT film with a low overpotential of 255 mV at a high current density of 1500 mA cm−2 in 1 M KOH. The strong chemical bonds between Mo2C/MoC and CNTs significantly weaken the dissolution and shedding of the Mo2C/MoC nanoparticles during the HER at high current densities. As a result, the overpotential of the Mo2C/MoC/CNT film changes by only ~32 and ~47 mV after working at 500 and 1000 mA cm−2 for 14 days, respectively. Density functional theory (DFT) calculations demonstrate the moderate free energy (ΔGH*) of 0.02 eV for hydrogen adsorption at sites around Mo2C/MoC interfaces and a strong coupling between the MoxC and CNT matrix, which ensures the high activity and stability of the heterogeneous Mo2C/MoC/CNT film.

The synthesis of Mo2C/MoC/CNT catalysts by self-heating is shown in Fig. 1a. First, precursors including ammonium molybdate as metal source and glucose and urea as carbon source were loaded onto a laser-drilled CNT film10 by dip coating (step 1). Second, the CNT film loaded with precursors was dried at 60 °C for 10 min in air (step 2). Finally, self-heating synthesis was performed in a mixed atmosphere of 10% H2 and 90% Ar with a total flow rate of 200 sccm, in which the precursors in situ reacted on the rapidly Joule-heated CNT film (step 3), resulting in a Mo2C/MoC/CNT composite film as an HER electrode (Fig. 1b). It took approximately ~256 ms to heat up from room temperature to ~1770 K and ~330 ms to cool down from ~1770 K to ~600 K (Fig. 1c and Supplementary Fig. 1a). The CNT films emit visible light after the application of a voltage of ~ 0.5 V/mm and show a uniform distribution of temperature (Fig. 1d and inset). The as-prepared Mo2C/MoC/CNT film remains flexible, akin to the pure CNT, MoC/CNT, and Mo2C/CNT films before and after the rapid treatment at ~1770 K (Supplementary Fig. 2). We simulated the in-plane temperature distribution of the CNT film at ~1700 K and the central temperature under the same power density with different CNT film sizes (by fixing the aspect ratio at 15:8) using the COMSOL Multi-Physics software. The constant temperature zone increases with the size of the CNT film and is adjustable to the application requirements (Fig. 1e and Supplementary Fig. 1b, d), although the temperature rapidly changes within ~1 mm at the boundary in contact with the graphite clips. This method is also suitable for large-scale fabrication (inset of Fig. 1e).

a Schematic illustration of the synthesis of Mo2C/MoC/CNT films by the self-heating method. b Mo2C/MoC/CNT film as an HER electrode. c Temperature-Time curve of heating (left) and cooling (right) process during self-heating. d Optical pictures of the CNT film at high temperature in a large size (80 mm × 40 mm) and a small size (15 mm × 8 mm, inset). e Simulation of temperature distribution (colour contour) at ~1700 K via COMSOL Multi-Physics software for the large-size film. The top-layer curve shows the temperature distribution along the centerline.

Raman spectroscopy and X-ray diffraction (XRD) spectroscopy were performed to analyse the chemical composition of the Mo2C/MoC/CNT film. The Mo2C/MoC/CNT film exhibits four apparent Raman peaks, of which two peaks at 822 and 995 cm−1 belong to β-Mo2C, and two peaks at ~1350 and ~1580 cm−1 belong to CNT (Fig. 2a)28,29,30. As shown in Fig. 2b and Supplementary Fig. 3a, the strong XRD peaks originate from β-Mo2C, α-MoC, and CNT. The fitted weight percentages of α-MoC and β-Mo2C are 59.8% and 40.2%, respectively, corresponding to an α-MoC-to-β-Mo2C molar ratio of about 2.8:1. Besides these peaks, there is a weak XRD peak at about 2θ = 11.5° (Supplementary Fig. 3b), corresponding to a reduced c lattice parameter of a transition phase from 2D Mo2C (MXene) to β-Mo2C. This may result from the defunctionalization of the surface terminations, the removal of intercalated water, and the reestablishment of long-range order in 2D Mo2C31,32. These results indicate that β-Mo2C and α-MoC are major phases in the composite and the transition phase from 2D Mo2C to β-Mo2C is minor. Furthermore, we qualitatively studied the trend of component changes through the peak intensity in the XRD patterns. The ratio of ammonium molybdate to glucose in the precursor solution adjusts the final loading content of Mo2C and MoC as well as the ratio of MoC to Mo2C. As shown in Fig. 2c and Supplementary Fig. 4, the heating time at 135 W also regulates MoC:Mo2C. The XRD peak intensity ratio of MoC (111) to Mo2C (002) decreases with increased heating time, which indicates an increased content of Mo2C or a decreased content of MoC and verifies the conversion of MoC to Mo2C at high temperatures. When the heating time is long enough, the Mo2C phase completely replaces the hybrid Mo2C/MoC, suggesting that a short reaction time is crucial for synthesizing the Mo2C/MoC heterostructure (Supplementary Fig. 5 and Supplementary Fig. 6).

a Raman spectra of Mo2C/MoC/CNT film, Mo2C powder, MoC/CNT film, and CNT film. The inset shows D and G Raman peaks of CNTs. b XRD pattern of Mo2C/MoC/CNT film. c XRD spectra of Mo2C/MoC/CNT films synthesized with different heating time at 135 W (left panel) and the corresponding peak intensity ratios of MoC (111) to Mo2C (002) (right panel). The error bars represent standard deviations. d XPS spectra of Mo 3d in MoC powder (top), Mo2C/MoC/CNT film (middle), and Mo2C powder (bottom). e XPS spectra of C 1 s in the Mo2C/MoC/CNT film. f TG analysis of Mo2C/MoC/CNT film in air.

X-ray photoelectron spectroscopy (XPS) was performed to analyse the electronic states. As shown in Supplementary Fig. 7a, clear Mo and C peaks appeared in the full spectra. The peaks of Mo 3d are divided into three groups, including two peaks of Mo2+ at 228.5 and 231.7 eV, two peaks of Mo3+ at 228.9 and 232.1 V, and two peaks of Mo4+ at 229.9 and 233.1 eV, respectively (Fig. 2d)33,34,35. It is still debating for the exact valence states of Mo in Mo2C or MoC. Most studies suggest the dominance of Mo2+ in Mo2C and Mo3+ in MoC, and higher valence states of Mo4+, Mo5+, and Mo6+ result from partial oxidation36,37. Under this assumption, compared with the Mo2+ peak in Mo2C and the Mo3+ peak in MoC, in the Mo2C/MoC/CNT film, the Mo2+ peak apparently blueshifts (~0.6 eV) while the Mo3+ peak redshifts (~0.5 eV), suggesting the existence of electron transfer from Mo2C to MoC in the heterogeneous composite. Note that the strong Mo4+, Mo5+, and Mo6+ peaks in the MoC sample should result from the surface MoOx because MoC is very prone to oxidation. Furthermore, we investigated the charge distribution in Mo2C and MoC by DFT calculations. As shown in Supplementary Fig. 8, each Mo atom loses about 0.80 and 0.50 electron, respectively, in the body of the MoC (111) and Mo2C (100) regions, while at the Mo2C/MoC interface, each Mo atom loses about 0.75 and 0.70 electron because of the electron transfer from Mo2C to MoC regions. This result clearly proves that the charge transfer exists solely at the Mo2C/MoC interface. In addition, evident peaks in the C 1 s spectra of C–C, C–O, O = C–O, and Mo–C are derived from CNT and MoxC (Fig. 2e)38. However, the peaks of N 1 s are not observed in the XPS spectra though urea is a precursor. This may be because N cannot replace C as a doping atom at a high temperature of ~1770 K (Supplementary Fig. 7b) during rapid heating. We also verified the absence of N using a material element analyser, as shown in Supplementary Table. 1. Thermogravimetric (TG) analysis was performed in air at 30–700 °C to estimate the content of each component in the Mo2C/MoC/CNT film (Fig. 2f). A significant increase in weight from ~250 to ~450 °C is mainly due to the oxidation of MoC and Mo2C to MoO3, followed by a weight loss caused by the oxidation of CNTs to CO2. Combining XRD and TGA data, the wt% of Mo2C, MoC, and CNTs in the Mo2C/MoC/CNT film is estimated to be ~15.4%, 22.8%, and 61.8%, respectively.

The surface morphology of the Mo2C/MoC/CNT film was observed using scanning electron microscopy (SEM). To improve the electrode stability, the CNT film was drilled by a laser to construct periodic micropores as the channels for H2 release (Supplementary Fig. 9)10. As shown in Fig. 3a, b the catalyst particles are uniformly dispersed on the surface of the CNT film without visible agglomeration. Mo2C/MoC particles, functioning as solder joints, bind several CNTs together to form a bundle, which tightly connect adjacent CNTs and prevent CNTs from sliding, greatly enhancing the binding force between adjacent CNTs and strengthening the CNT matrix. The small size (~20 nm) and distribution states of Mo2C/MoC particles facilitate the effective strengthening (Fig. 3c). We also changed the heating ramp time for both heating steps from room temperature to ~1100 K (30 W) and from ~1100 to ~1770 K (135 W). As the heating ramp rate decreases, the Mo2C/MoC particles gradually agglomerate and increase in size because they are more likely to diffuse and aggregate during the heating process, especially at high temperatures (Supplementary Fig. 10).

a, b SEM images of Mo2C/MoC/CNT film. c TEM images of Mo2C/MoC/CNT film. The inset is the statistics of Mo2C/MoC particle size distribution. d HRTEM images of the interface between Mo2C (100) and MoC (111) in Mo2C/MoC/CNT film. e HAADF-STEM image of Mo2C/MoC. The inset is HAADF-STEM image of MoC phase and Mo2C phase, as well as the FT patterns. Green dots represent Mo atoms, and blue dots represent C atoms. f SAED pattern of Mo2C/MoC/CNT film. g EDS mapping shows the distribution of Mo and C elements.

Distinct interfaces and two Mo-C phases in the particles were observed by high-resolution transmission electron microscopy (HRTEM, Fig. 3d). The (100) crystal plane of Mo2C forms an interface with the (111) crystal plane of MoC, whose interplanar distances of ~0.260 and ~0.246 nm are relatively similar. To further detect the phase structure and interface, high-angle annular dark field scanning transmission electron microscopy (HAADF-STEM) was performed, displaying Mo atoms because of a large Z value. Composite phases of Mo2C/MoC are formed in most catalytic particles and one typical particle is shown in Fig. 3e, where Mo2C, MoC, and an interfacial transition zone simultaneously exist. The inset (top right) of Fig. 3e and Supplementary Fig. 11a show the crystal structure of MoC and the atomic sites of Mo, where the angle of (111) and (220) is 90°, which is consistent with the measurement results in HAADF-STEM images and Fourier transform (FT) patterns. The inset (bottom right) of Fig. 3e and Supplementary Fig. 11b show the crystal structure of Mo2C, where the angles of (100) and (002) are 90° from the measurements, which is consistent with the theoretical prediction. Figure 3f shows the results of selected electron area diffraction (SAED), in which the diffraction spots/circles of Mo2C (101), MoC (111), and CNT (002) suggest the composition of Mo2C, MoC, and CNT. Energy-dispersive spectroscopy (EDS) mapping images show the distribution of Mo, C, and N in Fig. 3g and Supplementary Fig. 12. Mo atoms are concentrated on the catalyst particles, and C atoms are concentrated on catalyst particles and CNTs. However, no clear signal of N atoms was detected, which is consistent with the XPS spectra.

The electrochemical performance of the catalysts was measured in 1 M KOH (pH = 14) in a three-electrode configuration. Different Mo2C/MoC/CNT samples were prepared by changing the type of carbon sources, e.g. glucose, sodium phenolate, and phenylalanine, the content of the precursor, the weight percent of CNT, and the heating time (Supplementary Figs. 13–18). Glucose is a suitable carbon source owing to its high activity, easily breaking C-C bonds and forming Mo-C bonds within a short time (Supplementary Fig. 13). Considering the ratio of Mo to C in the precursor, an insufficient C source cannot react with the Mo source on time before the Mo source is evaporated out, leading to a small Mo2C/MoC loading content and insufficient Mo2C content. However, excess C source may cover Mo2C/MoC, preventing the active sites at Mo2C/MoC from thorough exposure (Supplementary Fig. 14). In our method, the light content of urea fine-tunes the atomic ratio of Mo to C, and the fast self-heating process adjusts the conversion from MoC to Mo2C, which further optimizes the ratio of MoC to Mo2C (Supplementary Figs. 15 and 16). A suitable proportion of MoC and Mo2C offers richer interfaces and provides more HER active sites. The optimal ratio of Mo to C source and glucose to urea are 4:3 and 20:3, respectively. Besides the ratios of Mo to C source and glucose to urea, the CNT wt% also plays a great role in the performance of the Mo2C/MoC/CNT films. As shown in Supplementary Fig. 17, the sample with the moderate CNT wt% (~61.8 wt% CNT and ~38.2 wt% Mo2C/MoC) possesses the best HER activity because it exhibits no apparent agglomeration of nanoparticles and thus has abundant Mo2C/MoC interfacial area. More or less content of Mo2C/MoC will lead to poorer HER activity. When the content of CNT is high, the content of MoC or Mo2C is very few, and thereby the total Mo2C/MoC interfacial area becomes very limited. And for the samples with a low content of CNT, Mo2C/MoC will agglomerate during the self-heating process, also reducing the total Mo2C/MoC interfacial area (Supplementary Fig. 18). This result suggests that an appropriate CNT wt% is needed to increase the Mo2C/MoC interfacial area for a better HER.

As shown in Fig. 4a, the as-prepared Mo2C/MoC/CNT film requires the overpotentials of only 82, 201, 233, and 255 mV to achieve HER current densities of 10, 500, 1000, and 1500 mA cm−2, respectively. Compared to precious metal electrodes, Mo2C/MoC/CNT films perform better at current densities over ~450 mA cm−2 than Pt/C, which is considered as a benchmark for HER catalysts. The overpotential of the Mo2C/MoC/CNT film is notably smaller than that of the Mo2C/CNT, MoC/CNT, and CNT films throughout the entire HER process, for which the overpotentials are ~318, >350, and >500 mV at 1000 mA cm−2, respectively. Compared to the Mo2C/CNT film, the MoC/CNT film, and especially the physically mixed Mo2C/MoC/CNT film (m-Mo2C/MoC/CNT film, Supplementary Fig. 19), the superiority of the Mo2C/MoC/CNT film made by self-heating suggests the importance of the Mo2C/MoC interface, where Mo is more conducive to the adsorption and desorption of H* via interfacial charge transfer33. Otherwise, either strong adsorption or strong desorption of H* weaken the activity of HER, leading to the lower electrochemical activity of the MoC/CNT and Mo2C/CNT films, respectively. To investigate the kinetics, Tafel slopes were calculated from the polarization curves (Supplementary Fig. 20). The Tafel slope of the Mo2C/MoC/CNT film is 42 mV dec−1, which is smaller than that of the Mo2C/CNT film (50 mV dec−1), MoC/CNT film (58 mV dec−1), and is close to that of Pt/C (31 mV dec−1), indicating that HER is based on the Volmer-Heyrovsky mechanism. As the kinetic activity is affected by the electrochemical surface area (ECSA) and charge transfer resistance, we measured the electrochemical double-layer capacitances (Cdl) of various electrodes, which were proportional to the ECSAs, and charge transfer resistances. As shown in Supplementary Figs. 21 and 22, the Mo2C/MoC/CNT film exhibits a high Cdl of 119.9 mF cm−2 and a small charge transfer resistance of ~3.2 Ω, suggesting an evident advantage compared to the Mo2C/CNT (83.5 mF cm−2 and ~4.5 Ω) and MoC/CNT (66.0 mF cm−2 and ~6.9 Ω) films. The polarization curves of the three kinds of films normalized to the ECSA and the mass of active materials are shown in Supplementary Fig. 23 (Cs is set to be 40 μF cm−2), which demonstrate that the activity of Mo2C/MoC/CNT film does increase intrinsically compared with those of Mo2C/CNT film and MoC/CNT film. For Cdl, the measurement error includes systematic errors (0.01% for voltage application and 0.2% for current detection) and random error (~5%), which is negligible. For Cs, as the actual value cannot be determined precisely in this work. Cs for a flat surface is generally found to be in a range of 20–60 μF cm−2, and the value of 40 μF cm−2 is used in this work to calculate the ECSA. The actual value of Cs could be very different from the used value and thus an error may be introduced in the absolute value of ECSA. However, the relative values of the ECSA for the materials reported in this work are not affected by the absolute errors because the Mo2C/MoC/CNT film, Mo2C/CNT film, and MoC/CNT film are prepared on the same supporting material, the CNT films. The Cs values for all these materials should be nearly identical in principle, and thus the calculated ECSAs can be compared relatively. Moreover, we also evaluated the turnover frequency (TOF) of each catalyst film (Supplementary Fig. 24). At an overpotential of 250 mV, the TOF of Mo2C/MoC/CNT film, Mo2C/CNT film, and MoC/CNT film is 0.65, 0.30, and 0.22 s−1, respectively, which validates that Mo2C/MoC/CNT catalyst has higher intrinsic activity besides the larger ECSA. The increase of the intrinsic activity should be attributed to the increase in the number of Mo2C/MoC interfaces.

a Polarization curves of Mo2C/MoC/CNT film, Mo2C/CNT film, MoC/CNT film, pure CNT film and Pt/C (20 wt% of Pt) in 1 M KOH. b Comparison of overpotentials at 1000 mA cm−2 versus productivity (synthesis cycles in a day) among HER catalysts maintaining activity at high current density. c XRD patterns of self-heating samples only use CNT as carbon sources. d Polarization curves of a pristine Mo2C/MoC/CNT film and working after 10000 CV cycles. e Long-term test of self-heating and furnace-heating samples at 500 or 1000 mA cm−2 without iR compensation (left). A mass loss of self-heating and furnace-heating sample after 5 min of ultrasonic treatment in 50 mL alcohol (right). The error bars represent standard deviations. f1 The overpotential of self-heating and furnace-heating samples at 500, 1000, and 1500 mA cm−2 in 1 M KOH, respectively. f2 The ratios of Δη/Δlog|j| for self-heating and furnace-heating samples in different current density ranges. f3 The mechanical tensile curves of self-heating sample, furnace-heating sample, and pure CNT film. g Comparison of long-term stability of various HER catalysts at both small and high current densities.

Besides the excellent HER activity of the Mo2C/MoC/CNT film at high current densities, our self-heating method has notable advantages in the rapid synthesis process and high productivity, which takes a comparable pre-/post-synthesis processing time and a much shorter synthesis time than traditional methods (Fig. 4b and Supplementary Table. 2). The η1000 of the Mo2C/MoC/CNT film is close to MoS2/Mo2C2, MoS2/Ni3S29, and is much lower than Ni2P/NF12, Ni2(1-x)Mo2xP-NF39, α-MoB240, NiCo LDH/NF7, Ni2P-Fe2P41, 2H-Nb1.35S242, HC-MoS2/Mo2C18, Sn-Ni3S25, LiCoBPO/NF43, and IrFe/NC44. The total growth time of one-cycle self-heating synthesis is only 45 s, which means that 1920 synthesis cycles can be theoretically completed within a day, more efficient than traditional furnace heating or solvothermal methods that require several hours for one cycle. Therefore, our self-heating synthesis not only delivers Mo2C/MoC/CNT films that have excellent HER activity and tolerate high current densities but also achieves high productivity beyond traditional methods.

To study the formation of chemical bonds between the CNTs and MoxC, we used only CNTs as the carbon source to conduct an experiment under the same reaction conditions without glucose and urea. The XRD patterns show the successful synthesis of MoxC (Fig. 4c), indicating that CNTs are involved in chemical reactions and Mo-C bonds are formed. Thus, the interaction between CNTs and MoxC is strengthened because of the Mo-C bonds, which significantly improve the electron transfer and stability. As a result, the polarization curve of the Mo2C/MoC/CNT film after 10,000 CV cycles remains almost unchanged compared with the original polarization curve, with a deviation of only 10 mV from 0 to 1500 mA cm−2 (Fig. 4d), which suggests that the Mo2C/MoC/CNT film is highly stable under CV cycle conditions. As shown in Supplementary Table. 3, the CV cycle stability of the Mo2C/MoC/CNT film is higher than all of the previous catalysts. We also performed 14-day chronopotentiometry (CP, at fixed current) at large current densities of 500 and 1000 mA cm−2 to further verify the stability (Fig. 4e, left). The Mo2C/MoC/CNT film synthesized by self-heating shows an increase of overpotential by only ~32 and ~47 mV after 14 days (336 h) at 500 and 1000 mA cm−2, respectively. In contrast, the overpotential of Mo2C/MoC/CNT film synthesized by tube furnace heating (denoted as f-Mo2C/MoC/CNT film, Supplementary Fig. 25) increases by ~344 mV after working at 1000 mA cm−2 for 14 days, which is much inferior to that of the Mo2C/MoC/CNT film by self-heating. The long-term stability of the Mo2C/MoC/CNT film at 3000 mA cm−2 is also much better than that of the f-Mo2C/MoC/CNT film (Supplementary Fig. 26).

It is found that the HER performance of the Mo2C/MoC/CNT catalyst was improved at the initial several hours (Supplementary Fig. 27), which may result from the fact that the surface of the as-prepared samples was reconstructed owing to the adsorption of oxygen-containing groups during the HER2,45. After a long-term test at a high current density, the content of Pt in the electrolyte is below the detection limit of inductively coupled plasma (ICP) mass spectrometry (Supplementary Table. 4), which excludes the influence of Pt dissolution from counter electrodes during the electrochemical measurements. A similar long-term test at the high current density using a graphite counter electrode also exhibits high stability until the graphite electrode is dissolved (Supplementary Fig. 28). We also used a graphite rod as the counter electrode to measure the CV curves of Mo2C/MoC/CNT films for 50 cycles, and then changed to use a Pt counter electrode for the other 50-cycles CV measurement. The CV curves obtained by graphite and Pt counter electrodes are almost identical (Supplementary Fig. 29). To further investigate the higher stability of the Mo2C/MoC/CNT catalyst beyond the f-Mo2C/MoC/CNT catalyst, we used an ultrasonic machine to treat the catalysts. As shown in Fig. 4e (right), the mass loss of the Mo2C/MoC/CNT film (3.2 ± 2.9%) after the ultrasonic treatment is lower than that of the f-Mo2C/MoC/CNT film (8.7 ± 1.1%), suggesting a stronger binding between Mo2C/MoC and CNTs in the Mo2C/MoC/CNT film. On the microscopic level, unlike bulk Mo2C/MoC with poor crystallinity covering CNTs in the f-Mo2C/MoC/CNT film, the well-crystalline Mo2C/MoC nanoparticles disperse in the Mo2C/MoC/CNT film (Supplementary Figs. 30 and 31), leading to huge effective contact areas and sufficient binding strength at the interfaces due to chemical Mo-C bonding. Such strengthened interfaces are sufficient to prevent the Mo2C/MoC from being peeled off by large quantities of H2 bubbles when working at high current densities. After working for ~6 days at 1000 mA cm−2, the Mo2C/MoC/CNT film changes slightly in the microscopic morphology of Mo2C/MoC particles, while the f-Mo2C/MoC/CNT film only has a few agglomerated Mo2C/MoC remained surrounding the CNTs (Supplementary Fig. 32). The XRD spectrum of the Mo2C/MoC/CNT film after working at 1000 mA cm−2 for ~6 days still shows clear MoC and Mo2C peaks as well as unnoticeable MoO2 or MoO3 peaks, although the ratio of the peak intensity of MoC to Mo2C changes (Supplementary Fig. 33) owing to the different etching and shedding rates of the two phases. Although the material is corroded, the remained excellent HER performance indicates that the Mo2C/MoC interface, rather than single MoC or Mo2C, is crucial in HER. In contrast, XRD peaks belonging to MoO2 and MoO3 appear in the f-Mo2C/MoC/CNT film after the test in the same period. This indicates that the f-Mo2C/MoC/CNT film is more prone to oxidation because of its poorer crystallinity, because in the furnace-heating, it is difficult to obtain high crystallinity of the f-Mo2C/MoC/CNT film and simultaneously maintain the Mo2C/MoC composite phase due to the much longer heating process at a lower temperature.

As shown in Fig. 4f1, the η500, η1000, and η1500 of the Mo2C/MoC/CNT film are lower than those of the f-Mo2C/MoC/CNT film by 14, 28, and 43 mV, respectively. Throughout HER, the Tafel slope of the Mo2C/MoC/CNT film increases slightly and is below 120 mV dec−1, indicating a rapid kinetic process even at a high current density. In contrast, the Tafel slope of the f-Mo2C/MoC/CNT film increases visibly and exceeds that of the Mo2C/MoC/CNT film, suggesting a higher reaction resistance (Fig. 4f2). Owing to the good dispersion of the catalyst particles and their tight combination with the CNT matrix, the self-heating samples exhibit a larger Cdl and smaller charge transfer resistance than that of the furnace heating ones (94.3 mF cm−2 and ~10.4 Ω). Therefore, the self-heating samples are preferable for HER (Supplementary Figs. 34 and 35). Direct mechanical tensile experiments visually illustrate that the Mo2C/MoC/CNT film has much higher Young's modulus (E = 231.60 MPa), breaking strength (σb = 6.87 MPa), and toughness (26.16 N m−3) than pure CNT film and f-Mo2C/MoC/CNT film. After the self-heating, the breaking strain of the pure CNT film is notably reduced while the breaking strength only slightly increases to ~0.32 MPa (Supplementary Fig. 36). In contrast, the Mo2C/MoC/CNT film synthesized by the self-heating process has a much higher breaking strength. These results suggest that the enhancement of the mechanical strength of the Mo2C/MoC/CNT film should be mostly attributed to the strong interaction between Mo2C/MoC and CNTs, which improves the load transfer efficiency inside the film.

We have enumerated the long-term stability of the Mo2C/MoC/CNT film and other HER catalysts. For a quantitative comparison, we introduce a parameter, Csta, to evaluate the stability, which represents the discharge of HER per unit area resulting from per millivolt of overpotential change and can be described by a function,

where j is current density, t is the whole working time, and ΔV is the variation of overpotential. A larger Csta indicates that more H2 bubbles are produced per unit area when the same degradation occurs in overpotential, and thus the corresponding catalyst has better stability. As shown in Fig. 4g and Supplementary Table. 5, the Csta of the Mo2C/MoC/CNT film is as large as 2.57 × 107 C cm−2 V−1, which is times or even orders of magnitude higher than the values of other high-performance HER catalysts including 1T-MoS2, MoNi4/MoO2, Ni2P-Fe2P, Co-NiS2, CoP/NiCoP/NC, S-MoS2@C, single atom NiI, etc. that were measured in alkali, as well as IrFe/NC and Mo2CTx/2H-MoS2 that were measured in acid4,8,9,10,41,44,46,47,48,49,50,51,52,53,54,55,56,57,58.

DFT calculations were performed to further investigate the mechanisms of activity and stability. In general, ΔGH* is considered as an important descriptor for analysing the HER activity. If hydrogen binds to the surface weakly, the adsorption step limits the catalytic activity, whereas, if the binding is too strong, the desorption step limits the activity. The ΔGH* should be close to zero for an optimal catalyst1. In the calculations, oxygen-terminated MoC (111) surface, Mo2C (100) surface, and Mo2C/MoC heterostructure were considered for HER, because the existence of high-valence Mo states (Mo4+, Mo5+, or Mo6+) in the XPS spectra (Fig. 2d) suggests partial surface oxidation of Mo2C, MoC, or Mo2C/MoC. In addition, the calculated Pourbaix diagram also shows that, under the real reaction condition, the surface of Mo2C is terminated by oxygen, which plays a key role in the catalytic performance of HER in alkaline media2. As shown in Fig. 5a, b, the values of ΔGH* for the Mo2C (100), MoC (111), and Mo2C (100)/MoC (111) heterostructures were calculated, where β-Mo2C and α-MoC phases were used for modelling as these two phases are the main phases as identified by XRD patterns. The values of ΔGH* for H and T sites on Mo2C (100) surface are 0.25 and 0.70 eV, respectively, which shows weak hydrogen adsorption. For the MoC (111) surface, a ΔGH* of −0.56 eV indicates strong hydrogen adsorption, which is not favourable for HER. For the Mo2C/MoC heterostructure, the values of ΔGH* for B1, B2, T1, and T2 sites on the Mo2C/MoC heterostructure are −0.14, −0.16, 0.02, and 0.08 eV, which is much closer to the ideal ΔGH* of 0 eV, offering more favourable hydrogen adsorption kinetics toward HER than that of Mo2C (100) and MoC (111). To reveal why the Mo2C/MoC interfaces are more active, we performed further DFT calculations to investigate the electronic structures of MoC (111), Mo2C (100), and Mo2C/MoC interface. The projected density of states (PDOS) of pz orbital of O atoms adsorbed on MoC (111), Mo2C (100), and Mo2C/MoC interface before and after H adsorption are shown in Supplementary Fig. 37. The pz-band energy center \(({\varepsilon }_{{p}_{z}})\) of the O atoms is investigated to describe the bonding strength of H* on the catalyst surface. The \({\varepsilon }_{{p}_{z}}\) of O atoms on the MoC (111), Mo2C (100), and Mo2C/MoC interface are −2.49, −3.86, and −3.09 eV, respectively, which shows that the binding of H* on the oxidized MoC (111) surface is the strongest, while for the oxidized Mo2C (100) it is the weakest. The binding of H* at the Mo2C/MoC interface is moderate, which results in excellent thermodynamic activity in HER. These results exhibit a similar trend with the HER activity in the experiment and thus verify that the Mo atoms at the interface are more HER active than Mo3+ in MoC or Mo2+ in Mo2C. It is worth noting that the Mo2C/MoC/CNT film synthesized by the self-heating method has numerous defects, dislocations, and twin boundaries, as shown in the HAADF-STEM image (Supplementary Fig. 38). DFT calculations show that the absolute values of ΔGH* for Mo2C (100) surface with carbon vacancy is almost unchanged compared to those of Mo2C (100) surface without carbon defects, showing the poor HER activity. However, the absolute value of ΔGH* for MoC (111) surface with carbon vacancy decreases and is closer to 0 eV than that of MoC (111) surface without carbon defects, showing that the carbon vacancy defects in the MoC (111) surface can improve HER activity. Although carbon vacancy may promote HER activity57, the ΔGH* for Mo2C/MoC interface is as low as 0.02 eV, which is closer to 0 eV than those of Mo2C (−0.35 or 0.65 eV for (100)) and MoC (−0.04 eV for (111)) surfaces with the carbon vacancy (Supplementary Table. 6 and Fig. 39). In experiments, the HER performance of the Mo2C/MoC composite is much better than the single Mo2C or MoC phase. Therefore, in the Mo2C/MoC/CNT film, we believe that the Mo2C/MoC interface is the main contribution to the HER activity, and the defects may only additionally promote the HER. The suitable ΔGH* of Mo2C/MoC and its unique structure indicate that the interface interaction between Mo2C and MoC as well as the synthesis by self-heating is favourable for the effective adsorption and activation of the reactant, enhancing the HER performance.

a Adsorption structures of hydrogen on the Mo2C (100), MoC (111) surfaces, and Mo2C (100)/MoC (111) heterostructure, in which the adsorption sites of hydrogen are shown in the yellow dotted circle. The cyan, brown, red and white spheres present the Mo, C, O, and H atoms, respectively (B: bridge site; T: top site; H: hollow site). b Free energy diagram of the HER for different adsorption sites on Mo2C (100), MoC (111), and Mo2C (100)/MoC (111) heterostructure. c Charge density difference of graphene/Mo2C (001) and graphene/MoC (111) heterostructure. The cyan and yellow regions represent charge depletion and accumulation, respectively, and the isosurface value is 0.005 e/Bohr3.

To theoretically study the stability, we establish the heterostructures of Mo2C (001) or MoC (111) slabs with graphene, instead of CNTs, as a simplified model, as graphene has similar properties to CNTs. The binding energy and charge density differences were calculated, as shown in Fig. 5c. The charges transfer from the MoxC slab to graphene and the binding energies for graphene/Mo2C and graphene/MoC are −62.11 and −52.26 meV/Å2, respectively, which is approximately three times larger than that of bilayer graphene (−23.01 meV/Å2). These results indicate strong coupling between the MoxC slab and graphene. Moreover, the shortest distance between Mo2C (MoC) and graphene is 2.28 (2.25) Å, which is close to the Mo–C bond length of Mo2C (2.11 Å) (2.17 Å for MoC), indicating that they are bonded to each other.

Moreover, self-heating synthesis is a universal synthesis method suitable for single-phase and composite carbides. As shown in Supplementary Fig. 40, we successfully synthesized well-dispersed single-phase Nb4C3 and W2C/WC composite nano-carbides by self-heating within several seconds and regulated the production by adjusting parameters, such as metal/non-metal sources, self-heating time, and output power with high productivity.

In summary, we developed a fast, self-heating (Joule heating) method using a CNT film as heat source and matrix to in-situ synthesize a highly HER active and robust Mo2C/MoC/CNT film catalyst in the presence of Mo and C precursors. Although Pt-group metal-free and high-current-density HER catalysts are much demanded in practical applications, their fast and scalable synthesis methods have still been very limited. In our method, the heating and cooling processes rapidly occur in hundreds of milliseconds, and the entire synthesis process only lasts for tens of seconds, very advantageous for large-scale and low-cost production. The productivity of the self-heating synthesis reaches about 2000 growth cycles per day, much higher than those of traditional methods such as furnace heating and hydrothermal synthesis.

The as-prepared Mo2C/MoC/CNT catalyst possesses ultra-low overpotentials of 201, 233, and 255 mV at 500, 1000, and 1500 mA cm−2, respectively, in 1 M KOH. The overpotential increases by only 47 mV after working for 14 days at 1000 mA cm−2, and changes by less than 10 mV after 10,000 CV cycles. The HER activity of the catalyst even keeps stable at 3000 mA cm−2 for days. Combining the results of experiments and DFT calculations, we reveal the mechanism of the composite catalyst for maintaining both high activity and high stability as follows: 1) The ultrafast heating and cooling and the short growth time facilitate the uniform dispersion of Mo2C/MoC nanoparticles and the formation of abundant Mo2C/MoC hetero-interfaces. 2) The charge transfer at the Mo2C/MoC hetero-interfaces results in the formation of a moderate H* adsorption of Mo between Mo3+ (strong adsorption of H*) and Mo2+ (strong desorption of H*) states. Therefore, Mo2C/MoC hetero-interfaces serve as the main HER active sites with a moderate ΔGH*, in which the adsorption and desorption proceed smoothly even at high current densities. 3) Strong chemical bonds form between Mo2C/MoC particles and CNT matrix due to the instantaneous high temperature, which prevents the catalyst particles from being peeled off, offering the Mo2C/MoC/CNT film with long-term stability at high current densities. All of these benefits are difficult to achieve by conventional methods and catalysts. Our study introduces a scalable method to develop electrodes with both high chemical activity and high mechanical stability, which can be applied to the production of industry-scale catalytic electrodes, to ultrafast material synthesis and strengthening, and to many other energy storage applications.

CNT film was made from super-aligned CNTs (SACNTs), which were synthesized by our previous work through low-pressure chemical vapour deposition59. 20 mg SACNTs were dispersed into 250 mL ethanol by high-power probe ultrasonication (SCIENTZ-950E) for 5 min at 600 W. The suspension solution was immediately transferred to the suction filter to form a CNT film. After sufficiently dried in the air, the CNT film was cut into rectangular shapes and drilled with many microscale holes by a direct laser writing machine (1064 nm in wavelength). The drilled holes have a diameter of ~40 mm and a pitch of 800 μm. This kind of holey CNT film effectively releases H2 bubbles during HER, as revealed in our previous work10.

To synthesize Mo2C/MoC/CNT film, ammonium molybdate ((NH4)2MoO4·4H2O) and glucose (C6H12O6) with varied atomic ratios of Mo:C (1:1, 4:3, 2:1, and 4:1) were dissolved in a mixed solution of deionized water and ethanol. Urea (CH4N2O) was then dissolved into the mixed solution based on an optimal atomic ratio of Mo:C (4:3 in the experiment) and varied mole ratios of glucose to urea (0, 10:1, 20:3, 5:1, 5:2). To promote the dissolution of (NH4)2MoO4·4H2O and prevent the solution from precipitation, we added ammonia water to adjust the pH value of the solution to about 11.5. The precursor was loaded on a CNT film through three times of dip coating and then dried at 60 °C for 10 min in air. The CNT film loaded with precursor was then clamped at both ends by graphite clips and put into a quartz tube. The graphite clips were connected to a stabilized power supply through copper foil. The self-heating process of the CNT film took place in a reductive, mixed atmosphere of 10% H2 and 90% Ar with a certain voltage applied to both ends of the CNT film by a stabilized power supply outputting 30 W or 135 W. The heating time was 30 s at 30 W and ranged from 10 to 120 s at 135 W. The power supply kept an output of 30 W for 30 s at first for carbonization of glucose/urea and decomposition of ammonium molybdate, which prevented H2O molecules generated by the decomposition of glucose and ammonium molybdate from etching and disintegrating CNTs at high temperatures. Then the output was increased to 135 W and maintained for a short period (10–120 s) to form Mo-C bonds. The as-synthesized Mo2C/MoC/CNT film is washed with deionized water and ethanol orderly several times and dried before the following electrochemical tests. Besides the self-heating approach to Mo2C/MoC/CNT film, we also synthesize f-Mo2C/MoC/CNT films by furnace, with the same pretreatment process but a tube-furnace growth as shown in Supplementary Fig. 25. The loading of the Mo2C/MoC on the Mo2C/MoC/CNT film is ~1.8 mg cm−2, which indicated ~38 wt % of active material. The thickness is ~0.05 mm and the size is 15 mm × 8 mm.

To synthesize Mo2C/CNT film, the same pretreatment process was carried out as did in the synthesis of Mo2C/MoC/CNT film, except that the self-heating growth process was set at 30 W for 20 min. To synthesize MoC/CNT film, the self-heating growth process was set at 30 W for 45 s. The loading of the active materials, i.e., Mo2C on the Mo2C/CNT film, and MoC on the MoC/CNT film, are ~2.0 and ~1.7 mg cm−2, respectively.

Morphology of samples was observed by scanning electron microscope (SEM, ZEISS, Merlin Compact), transmission electron microscopy (TEM, JEOL, JEM-2010F, 200 kV), and spherical correction transmission electron microscope (sc-TEM, JEOL, JEM ARM 200 F, 200 kV). Surface states of the samples were analysed by X-Ray diffraction (XRD, D/max-2500/PC, Rigaku) with Cu Kα radiation (λ = 0.15406 nm) operated at 40 kV and 150 mA. Raman spectra were collected in the range of 200–800 cm−1 with an excitation wavelength of 532 nm (Horiba-iHR550) and X-ray photoelectron spectroscopy (XPS, Thermo Fisher, Escalab 250Xi, Al Kα) was carried out in the range of 0 to 1350 eV at a step of 1 eV. Thermogravimetry (TG, TA INSTRUMENTS, Q5000IR) was performed from 25 to 700 °C in air with a heating rate of 10 °C min−1. The temperature of CNT matrix was measured by an infrared thermometer (Optrics, PI640).

All electrochemical measures were performed on CHI 660e and CHI 760e electrochemical workstation by a standard three-electrode test. The CNT composite films directly served as the working electrodes. Hg/HgO served as the reference electrode and Pt as a counter electrode. All the potentials were converted to RHE. The polarization curves were measured at 5 mV s−1 with an 85% iR compensation. The cyclic voltammetry curves for fitting double-layer capacitance (Cdl) were measured at 10–50 mV s−1 from 0–0.1 V vs. RHE. The CV cycles for stability test ranged from 0–250 mV with a scan rate of 50 mV s−1. Chronopotentiometry (CP, fixed current) was performed at 500, 1000, or 3000 mA cm−2 for several days without iR compensation. In order to maintain the pH and liquid level of the electrolyte in the CP test, we introduced a microinjector to continuously replenish deionized water. The noble metal electrodes of Pt/C (20 wt%) were prepared with loading of ~1.9 mg/cm−2 on a CNT film via dropping catalyst ink containing Pt/C powder, 50 μL Nafion, 500 μL ethanol, and 450 μL distilled water. Electrochemical impedance spectroscopy (EIS) measurements were tested at an overpotential of 150 mV with a frequency ranging from 105 to 10−2 Hz with an AC amplitude of 5 mV.

The simulation via COMSOL Multi-Physics software included a geometric model, electric current module, heat transfer module, and electromagnetic heat module. The geometric model consisted of a CNT film, two copper foils, and four graphite clips with the same dimensions as used in experiments. The conductivity, thermal conductivity and density of CNT were set as 3136 S m−1, 28 × (300/T) W m−1 K−1, and 90 kg m−3, respectively. These physical parameters were set according to the experimental parameters, COMSOL database, and literature60. When investigating the central temperature of different-length CNT films, the power density was fixed at ~5.125 W mm−3.

All calculations based on density functional theory (DFT) were performed using the Vienna ab initio simulation package (VASP)61. The projector augmented wave (PAW) potentials62 and generalized gradient approximation (GGA) of the Perdew-Burke-Ernzerhof (PBE) functional63 were used for the electron-ion interaction and exchange-correlation energy, respectively. A 3 × 2 × 1 Mo2C (100) slab consisting of six layers of Mo and three layers of C, and a 2 × 2 × 1 MoC (111) slab composed of three layers of Mo and three layers of C were constructed to investigate the HER catalytic performance. Moreover, a 3 × 3 × 1 Mo2C (001) slab consisting of three layers Mo2C and a 2 × 2 × 1 MoC (111) slab were taken into account to establish the heterostructures with a 5 × 5 × 1 graphene supercell. During structural relaxation, the bottom two layers of Mo and C for the Mo2C (100) and MoC (111) and the bottom Mo2C layer for Mo2C (001) slab were fixed. The cutoff energy of the plane wave basis was set to 400 eV. The convergence criteria for the total energy and force were set to 10−5 eV and 0.01 eV Å−1, respectively. The vacuum layer of at least 15 Å in the z direction was selected to eliminate the interactions between the periodic images. A dipole correction along the z direction of the slab was used in all calculations. The DFT-D3 method64 was used to describe the Van der Waals interaction. The Gibbs free energy of adsorption hydrogen (ΔGH*) is calculated using ΔGH* = ΔEH* + ΔEZPE − TΔS, where ΔEH*, ΔEZPE and ΔS are the adsorption energy, zero-point energy change and entropy change of hydrogen adsorption, respectively. T is the temperature (T = 298.15 K). Moreover, according to our calculated results, the values of ΔGH* without the implicit solvent environment show the same conclusion as the implicit solvent environment as implemented in VASPsol65.

The data supporting this study are available within the paper and the Supplementary Information. All other relevant source data are available from the corresponding authors upon reasonable request. Source data are provided with this paper.

Seh, Z. W. et al. Combining theory and experiment in electrocatalysis: Insights into materials design. Science 355, eaad4998 (2017).

Article PubMed Google Scholar

Luo, Y. et al. Morphology and surface chemistry engineering toward pH-universal catalysts for hydrogen evolution at high current density. Nat. Commun. 10, 269 (2019).

Article ADS PubMed PubMed Central CAS Google Scholar

Yang, H., Chen, Z., Guo, P., Fei, B. & Wu, R. B-doping-induced amorphization of LDH for large-current-density hydrogen evolution reaction. Appl Catal. B. 261, 118240 (2020).

Article CAS Google Scholar

Zhao, Y. et al. Non-metal Single-Iodine-Atom Electrocatalysts for the Hydrogen Evolution Reaction. Angew. Chem. Int Ed. 58, 12252–12257 (2019).

Article CAS Google Scholar

Jian, J. et al. Sn-Ni3S2 Ultrathin Nanosheets as Efficient Bifunctional Water-Splitting Catalysts with a Large Current Density and Low Overpotential. ACS Appl Mater. Interfaces. 10, 40568–40576 (2018).

Article CAS PubMed Google Scholar

Staffell, I. et al. The role of hydrogen and fuel cells in the global energy system. Energy Environ. Sci. 12, 463–491 (2019).

Article CAS Google Scholar

Zhu, W. et al. NiCo/NiCo–OH and NiFe/NiFe–OH core shell nanostructures for water splitting electrocatalysis at large currents. Appl Catal. B. 278, 119326 (2020).

Article CAS Google Scholar

Jia, Z. et al. A Novel Multinary Intermetallic as an Active Electrocatalyst for Hydrogen Evolution. Adv. Mater. 32, e2000385 (2020).

Article PubMed CAS Google Scholar

Xue, S., Liu, Z., Ma, C., Cheng, H.-M. & Ren, W. A highly active and durable electrocatalyst for large current density hydrogen evolution reaction. Sci. Bull. 65, 123–130 (2020).

Article CAS Google Scholar

Li, C. et al. A lightly Fe-doped (NiS2/MoS2)/carbon nanotube hybrid electrocatalyst film with laser-drilled micropores for stabilized overall water splitting and pH-universal hydrogen evolution reaction. J. Mater. Chem. A. 8, 17527–17536 (2020).

Article CAS Google Scholar

Tang, T. et al. Electronic and Morphological Dual Modulation of Cobalt Carbonate Hydroxides by Mn Doping toward Highly Efficient and Stable Bifunctional Electrocatalysts for Overall Water Splitting. J. Am. Chem. Soc. 139, 8320–8328 (2017).

Article CAS PubMed Google Scholar

Yu, X. et al. "Superaerophobic" Nickel Phosphide Nanoarray Catalyst for Efficient Hydrogen Evolution at Ultrahigh Current Densities. J. Am. Chem. Soc. 141, 7537–7543 (2019).

Article CAS PubMed Google Scholar

Hu, C. et al. Synergism of Geometric Construction and Electronic Regulation: 3D Se-(NiCo)Sx/(OH)x Nanosheets for Highly Efficient Overall Water Splitting. Adv. Mater. 30, e1705538 (2018).

Article PubMed CAS Google Scholar

Yu, Q. et al. A Ta-TaS2 monolith catalyst with robust and metallic interface for superior hydrogen evolution. Nat. Commun. 12, 6051 (2021).

Article ADS CAS PubMed PubMed Central Google Scholar

Liu, Y. et al. Self-optimizing, highly surface-active layered metal dichalcogenide catalysts for hydrogen evolution. Nat. Energy. 2, 17172 (2017).

Article Google Scholar

Lu, Z. et al. Ultrahigh hydrogen evolution performance of under-water "superaerophobic" MoS2 nanostructured electrodes. Adv. Mater. 26, 2683–2687 (2014).

Article ADS CAS PubMed Google Scholar

Jia, J. et al. Ultrathin N-Doped Mo2C nanosheets with exposed active sites as efficient electrocatalyst for hydrogen evolution reactions. ACS Nano. 11, 12509–12518 (2017).

Article CAS PubMed Google Scholar

Zhang, C. et al. High-throughput production of cheap mineral-based two-dimensional electrocatalysts for high-current-density hydrogen evolution. Nat. Commun. 11, 3724 (2020).

Article ADS CAS PubMed PubMed Central Google Scholar

Li, C. et al. Flexible and free-standing hetero-electrocatalyst of high-valence-cation doped MoS2/MoO2/CNT foam with synergistically enhanced hydrogen evolution reaction catalytic activity. J. Mater. Chem. A 8, 14944–14954 (2020).

Article CAS Google Scholar

Park, H. et al. Canonic-Like HER activity of Cr1-xMoxB2 solid solution: overpowering Pt/C at high current density. Adv. Mater. 32, e2000855 (2020).

Article PubMed CAS Google Scholar

Peng, M. et al. 3D printed mechanically robust graphene/CNT electrodes for highly efficient overall water splitting. Adv. Mater. 32, e1908201 (2020).

Article PubMed CAS Google Scholar

Yang, W. et al. Conversion of intercalated MoO3 to multi-heteroatoms-doped MoS2 with high hydrogen evolution activity. Adv. Mater. 32, e2001167 (2020).

Article PubMed CAS Google Scholar

Zhuo, S. et al. Dual-template engineering of triple-layered nanoarray electrode of metal chalcogenides sandwiched with hydrogen-substituted graphdiyne. Nat. Commun. 9, 3132 (2018).

Article ADS PubMed PubMed Central CAS Google Scholar

Zhang, J., Zhang, Q. & Feng, X. Support and interface effects in water-splitting electrocatalysts. Adv. Mater. 31, e1808167 (2019).

Article PubMed CAS Google Scholar

Xie, H. et al. High-temperature pulse method for nanoparticle redispersion. J. Am. Chem. Soc. 142, 17364–17371 (2020).

Article CAS PubMed Google Scholar

Wang, C. et al. A general method to synthesize and sinter bulk ceramics in seconds. Science 368, 521–526 (2020).

Article ADS CAS PubMed Google Scholar

Yao, Y. et al. High temperature shockwave stabilized single atoms. Nat. Nanotechnol. 14, 851–857 (2019).

Article ADS CAS PubMed Google Scholar

Huang, Y. et al. Mo2C nanoparticles dispersed on hierarchical carbon microflowers for efficient electrocatalytic hydrogen evolution. ACS Nano. 10, 11337–11343 (2016).

Article CAS PubMed Google Scholar

Kim, H., Anasori, B., Gogotsi, Y. & Alshareef, H. N. Thermoelectric properties of two-dimensional molybdenum-based MXenes. Chem. Mater. 29, 6472–6479 (2017).

Article CAS Google Scholar

Geng, D. et al. Direct Synthesis Of Large-area 2D Mo2C on in situ grown graphene. Adv. Mater. 29, 1700072 (2017).

Article CAS Google Scholar

Zhou, H. et al. Two-dimensional molybdenum carbide 2D-Mo2C as a superior catalyst for CO2 hydrogenation. Nat. Commun. 12, 5510 (2021).

Article ADS CAS PubMed PubMed Central Google Scholar

Halim, J. et al. Synthesis and characterization of 2D molybdenum carbide (MXene). Adv. Funct. Mater. 26, 3118–3127 (2016).

Article CAS Google Scholar

Zhang, X. et al. Structure and phase regulation in MoxC (α-MoC1-x/β-Mo2C) to enhance hydrogen evolution. Appl Catal. B. 247, 78–85 (2019).

Article CAS Google Scholar

Ge, Y. et al. Maximizing the synergistic effect of CoNi catalyst on α-MoC for robust hydrogen production. J. Am. Chem. Soc. 143, 628–633 (2021).

Article CAS PubMed Google Scholar

Yao, S. et al. Atomic-layered Au clusters on α-MoC as catalysts for the low-temperature water-gas shift reaction. Science 357, 389–393 (2017).

Article ADS CAS PubMed Google Scholar

Murugappan, K. et al. Operando NAP-XPS unveils differences in MoO3 and Mo2C during hydrodeoxygenation. Nat. Catal. 1, 960–967 (2018).

Article CAS Google Scholar

Liu, W. et al. A durable and pH-universal self-standing MoC-Mo2C heterojunction electrode for efficient hydrogen evolution reaction. Nat. Commun. 12, 6776 (2021).

Article ADS CAS PubMed PubMed Central Google Scholar

Yu, B. et al. Mo2C nanodots anchored on N-doped porous CNT microspheres as electrode for efficient Li-Ion storage. Small Methods. 3, 1800287 (2019).

Article CAS Google Scholar

Yu, L. et al. Ternary Ni2(1-x)Mo2xP nanowire arrays toward efficient and stable hydrogen evolution electrocatalysis under large-current-density. Nano Energy. 53, 492–500 (2018).

Article CAS Google Scholar

Chen, Y. et al. Highly active, nonprecious electrocatalyst comprising borophene subunits for the hydrogen evolution reaction. J. Am. Chem. Soc. 139, 12370–12373 (2017).

Article CAS PubMed Google Scholar

Wu, L. et al. Heterogeneous bimetallic phosphide Ni2P‐Fe2P as an efficient bifunctional catalyst for water/seawater splitting. Adv. Funct. Mater. 31, 2006484 (2021).

Article CAS Google Scholar

Yang, J. et al. Ultrahigh-current-density niobium disulfide catalysts for hydrogen evolution. Nat. Mater. 18, 1309–1314 (2019).

Article ADS CAS PubMed Google Scholar

Lei, C. et al. Efficient alkaline hydrogen evolution on atomically dispersed Ni–Nx species anchored porous carbon with embedded Ni nanoparticles by accelerating water dissociation kinetics. Energy Environ. Sci. 12, 149–156 (2019).

Article CAS Google Scholar

Jiang, P. et al. Improving electrocatalytic activity of iridium for hydrogen evolution at high current densities above 1000 mA cm−2. Appl Catal. B 258, 117965 (2019).

Article CAS Google Scholar

He, L. et al. Molybdenum carbide-oxide heterostructures: in situ surface reconfiguration toward efficient electrocatalytic hydrogen evolution. Angew. Chem. Int Ed. 59, 3544–3548 (2020).

Article CAS Google Scholar

Chen, W. et al. Achieving rich and active alkaline hydrogen evolution heterostructures via interface engineering on 2D 1T‐MoS2 quantum sheets. Adv. Funct. Mater. 30, 2000551 (2020).

Article CAS Google Scholar

Senthil Raja, D., Chuah, X.-F. & Lu, S.-Y. In situ grown bimetallic MOF-based composite as highly efficient bifunctional electrocatalyst for overall water splitting with ultrastability at high current densities. Adv. Energy Mater. 8, 1801065 (2018).

Article CAS Google Scholar

Li, S. et al. Bifunctional CoNi/CoFe2O4/Ni foam electrodes for efficient overall water splitting at a high current density. J. Mater. Chem. A 6, 19221–19230 (2018).

Article CAS Google Scholar

Zhang, H. et al. Continuous oxygen vacancy engineering of the Co3O4 layer for an enhanced alkaline electrocatalytic hydrogen evolution reaction. J. Mater. Chem. A 7, 13506–13510 (2019).

Article CAS Google Scholar

Shen, F. et al. Bimetallic iron-iridium alloy nanoparticles supported on nickel foam as highly efficient and stable catalyst for overall water splitting at large current density. Appl Catal. B 278, 119327 (2020).

Article CAS Google Scholar

Zhang, J. et al. Efficient hydrogen production on MoNi4 electrocatalysts with fast water dissociation kinetics. Nat. Commun. 8, 1–8 (2017).

ADS CAS Google Scholar

Lim, K. R. G. et al. 2H-MoS2 on Mo2CTx MXene nanohybrid for efficient and durable electrocatalytic hydrogen evolution. ACS Nano. 14, 16140–16155 (2020).

Article PubMed CAS Google Scholar

Tang, C., Cheng, N., Pu, Z., Xing, W. & Sun, X. NiSe nanowire film supported on nickel foam: an efficient and stable 3d bifunctional electrode for full water splitting. Angew. Chem. Int Ed. 127, 9483–9487 (2015).

Article ADS Google Scholar

Yin, J. et al. Atomic arrangement in metal-doped NiS2 boosts the hydrogen evolution reaction in alkaline media. Angew. Chem. Int Ed. 58, 18676–18682 (2019).

Article CAS Google Scholar

Boppella, R., Tan, J., Yang, W. & Moon, J. Homologous CoP/NiCoP heterostructure on N‐doped carbon for highly efficient and pH‐universal hydrogen evolution electrocatalysis. Adv. Funct. Mater. 29, 1807976 (2019).

Article CAS Google Scholar

Li, R.-Q. et al. Monolithic electrode integrated of ultrathin NiFeP on 3D strutted graphene for bifunctionally efficient overall water splitting. Nano Energy 58, 870–876 (2019).

Article CAS Google Scholar

Xu, Q. et al. Unsaturated sulfur edge engineering of strongly coupled mos2 nanosheet-carbon macroporous hybrid catalyst for enhanced hydrogen generation. Adv. Energy Mater. 9, 1802553 (2019).

Article CAS Google Scholar

Ibupoto, Z. H. et al. MoSx@NiO composite nanostructures: an advanced nonprecious catalyst for hydrogen evolution reaction in alkaline media. Adv. Funct. Mater. 29, 1807562 (2019).

Article CAS Google Scholar

Zhang, X. et al. Spinning and processing continuous yarns from 4-inch wafer scale super-aligned carbon nanotube arrays. Adv. Mater. 18, 1505–1510 (2006).

Article CAS Google Scholar

Yu, W., Liu, C. & Fan, S. Advances of CNT-based systems in thermal management. Nano Res. 14, 2471–2490 (2021).

Article ADS Google Scholar

Kresse, G. & Furthmuller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B. 54, 11169 (1996).

Article ADS CAS Google Scholar

Blochl, P. E. Projector augmented-wave method. Phys. Rev. B. 50, 17953–17979 (1994).

Article ADS CAS Google Scholar

Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865 (1996).

Article ADS CAS PubMed Google Scholar

Grimme, S., Antony, J., Ehrlich, S. & Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 132, 154104 (2010).

Article ADS PubMed CAS Google Scholar

Mathew, K., Sundararaman, R., Letchworth-Weaver, K., Arias, T. A. & Hennig, R. G. Implicit solvation model for density-functional study of nanocrystal surfaces and reaction pathways. J. Chem. Phys. 140, 084106 (2014).

Article ADS PubMed CAS Google Scholar

Download references

K.L. acknowledges the financial support from the National Key R&D Program of China (2018YFA0208400), Basic Science Center Project of NSFC under grant No. 51788104, and National Natural Science Foundation of China (51972193). J.L. acknowledges the support from the National Natural Science Foundation of China (11874036), Local Innovative and Research Teams Project of Guangdong Pearl River Talents Program (2017BT01N111), and Basic Research Project of Shenzhen, China (JCYJ20200109142816479, WDZC20200819115243002). Y.S. acknowledges the support from the National Natural Science Foundation of China (11974041).

These authors contributed equally: Chenyu Li, Zhijie Wang, Mingda Liu.

State Key Laboratory of New Ceramics and Fine Processing, School of Materials Science and Engineering, Tsinghua University, Beijing, 100084, China

Chenyu Li, Mingda Liu, Enze Wang, Bolun Wang, Longlong Xu & Kai Liu

Shenzhen Geim Graphene Center and Institute of Materials Research, Tsinghua Shenzhen International Graduate School, Tsinghua University, Shenzhen, 518055, China

Zhijie Wang & Jia Li

Department of Physics and Tsinghua-Foxconn Nanotechnology Research Center, Tsinghua University, Beijing, 100084, China

Kaili Jiang & Shoushan Fan

Beijing Key Laboratory for Magneto-Photoelectrical Composite and Interface Science, School of Mathematics and Physics, University of Science and Technology Beijing, Beijing, 100083, China

Yinghui Sun

You can also search for this author in PubMed Google Scholar

K.L. and C.L. conceived the project and designed experiments. C.L., M.L., L.X., and Y.S. fabricated and tested samples. M.L. completed the COMSOL simulation. J.L. and Z.W. performed the DFT calculations. C.L., E.W., and B.W. discussed the working mechanism. C.L., K.L., J.L., Y.S., K.J., and S.F. analysed the data. C.L., K.L., Z.W., J.L., and Y.S. wrote and revised the manuscript. All authors discussed the results and contributed to the final version of the manuscript.

Correspondence to Yinghui Sun, Jia Li or Kai Liu.

The authors declare no competing interests.

Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. Peer reviewer reports are available.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Li, C., Wang, Z., Liu, M. et al. Ultrafast self-heating synthesis of robust heterogeneous nanocarbides for high current density hydrogen evolution reaction. Nat Commun 13, 3338 (2022). https://doi.org/10.1038/s41467-022-31077-x

Download citation

Received: 13 August 2021

Accepted: 31 May 2022

Published: 09 June 2022

DOI: https://doi.org/10.1038/s41467-022-31077-x

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

Nano Research (2023)

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.