Guangfu Qian ,Jinli Chen ,Tianqi Yu ,Jiacheng Liu ,Lin Luo ,Shibin Yin

ABSTRACT Constructing heterojunction is an effective strategy to develop high-performance non-preciousmetal-based catalysts for electrochemical water splitting(WS).Herein,we design and prepare an N-doped-carbon-encapsulated Ni/MoO2 nano-needle with three-phase heterojunction (Ni/MoO2@CN) for accelerating the WS under industrial alkaline condition.Density functional theory calculations reveal that the electrons are redistributed at the three-phase heterojunction interface,which optimizes the adsorption energy of H-and O-containing intermediates to obtain the best ΔGH* for hydrogen evolution reaction (HER) and decrease the ΔG value of ratedetermining step for oxygen evolution reaction (OER),thus enhancing the HER/OER catalytic activity.Electrochemical results confirm that Ni/MoO2@CN exhibits good activity for HER (?-10=33 mV,?-1000=267 mV)and OER (?10=250 mV,?1000=420 mV).It shows a low potential of 1.86 V at 1000 mA cm-2 for WS in 6.0 M KOH solution at 60 °C and can steadily operate for 330 h.This good HER/OER performance can be attributed to the three-phase heterojunction with high intrinsic activity and the self-supporting nano-needle with more active sites,faster mass diffusion,and bubbles release.This work provides a unique idea for designing high efficiency catalytic materials for WS.

KEYWORDS Three-phase heterojunction;Interfacial electronic redistribution;Catalyst;Large current density;Water splitting

1 Introduction

Electrochemical water splitting (WS) is considered as a green and convenient method for producing hydrogen (H2)[1-4].However,the practical water electrolysis is limited by the sluggish kinetics at cathode (hydrogen evolution reaction,HER) and anode (oxygen evolution reaction,OER)[5-11].For industrial WS catalyst,these requirements must be met:(1) stable and high active sites at large current density [12];(2) continuous intensive gas evolution [13,14];(3)fast electron transfer [15,16];(4) low-cost and accessibility[17,18].Although some progress has been made,it is still a challenge to explore high-performance WS catalysts,especially in industrial environment.

The strategies,such as constructing heterojunction [19],fabricating carbon-encapsulating structure [20,21],and manufacturing defects [22],are demonstrated to improve the WS performance of 3D transition-metals-based (TMB) catalysts.Notably,heterojunction is widely studied for enhancing the WS intrinsic activity through adjusting the electronic structure of catalysts [23,24].Mu et al.synthesized Mo-doped Ni3S2/NixPyheterojunction,which optimizes the absorption energy of H-and O-containing intermediates to promote the WS catalytic process [25].Jiao et al.[26] prepared the heterojunction of Co and nitrogen-rich nitride,which induces its interfacial electronic redistribution to obtain the optimal H* absorption and decrease the dissociation energy barrier.

Furthermore,carbon-encapsulated 3D TMB materials as a kind of heterojunction can not only modulate the electronic structure to boost HER/OER catalytic activity,but also avoid metal corrosion in harsh electrolyte to improve the stability [27].Deng et al.[28] proved that the surface electronic structure of FeNi alloy with ultrathin graphene layer is adjusted to obtain the good OER performance.Despite the significant progress,most of them are still far from the practical current density (≥ 500 mA cm-2).In our previous work,the graphene-coated NiCo alloy coupled with NiCoMoO and self-supported on 3D nickel foam (NiCo@C-NiCoMoO/NF)was synthesized to improve the intrinsic activity and stability at large current density (≥ 500 mA cm-2),which is different from other carbon-encapsulated 3D TMB materials due to the unique metal oxide anchoring alloy structure.Nevertheless,the understanding of its active sites is still not deep enough.Thus,developing 3D TMB materials with excellent WS performance at large current density and investigating the mechanism of its high intrinsic activity are significant but challenging.

Besides,in order to improve WS stability at large current density,constructing a reasonable structure (nanowires,nanosheets,etc.) on 3D porous substrate (e.g.,NF) is significant,which can accelerate the mass transfer and bubbles escape during the electrochemical reaction process,and provide large specific surface area to expose abundant active sites [15,29].MoS2/Ni3S2heterostructure nanowires selfsupported on NF were constructed by Cheng et al.[30].It exhibits a small overpotential of 182 mV at 500 mA cm-2for HER,which is comparable to Pt/C.

Inspired by the above strategies,the N-doped-carboncoated Ni/MoO2nano-needle with three-phase heterojunction (Ni/MoO2@CN) self-supported on 3D porous NF is rationally engineered and synthesized for the first time,and the HER/OER intrinsic activity is further investigated by density functional theory (DFT) calculations.This novel structure has the following advantages:(1) The electronic structure of N-doped-carbon (CN),Ni,and MoO2can be optimized at three-phase heterojunction interface for adjusting the adsorption energy of H-and O-containing intermediates,which can be beneficial to obtain the best H* Gibbs free energy (ΔGH*) for HER and decrease the Gibbs free energy (ΔG) value of rate-determining step (RDS) for OER,thus enhancing the HER/OER intrinsic activity and electron transfer ability under a large current density [31].(2) CN can prevent metal corrosion in electrolyte to boost catalytic stability [27,32].(3) Nano-needles self-supported on 3D NF with meso/macroporous structure can provide large specific surface area to expose more active sites for reducing the energy barriers of WS [33].Furthermore,it is beneficial to mass diffusion and bubbles release at large current density,thus improving the WS performance [34].

Benefitting from the synergistic effect among the CN,Ni,and MoO2at the three-phase heterojunction interface,as well as the unique self-supporting nano-needle structure,the prepared Ni/MoO2@CN exhibits low overpotentials at ± 10 and± 1000 mA cm-2for HER (?-10=33 mV,?-1000=267 mV)and OER (?10=250 mV,?1000=420 mV).As for WS,it shows a low potential (1.86 V) at 1000 mA cm-2in 6.0 M KOH solution at 60 °C with an ignorable depletion of activity after 330 h testing,showing excellent stability.Therefore,this work provides a promising catalyst for WS,and it sheds light on developing highly efficient non-precious-metal materials for industrial water electrolysis.

2 Result and Discussion

2.1 DFT Calculations for HER and OER

Firstly,the computational models of Ni/MoO2@CN,Ni/MoO2,Ni@CN,MoO2@CN,and CN are designed (the details are displayed in Fig.S1).Then,the three-phase heterojunction interface synergy effect and HER activity are studied via DFT computations.Particularly,the ΔGH*and electron distribution for the active atoms are investigated by DFT and density of states (DOS) [35].Hence,we consider all possible active sites of Ni/MoO2@CN for HER and calculate their ΔGH*(Figs.1a and S2).As shown in Fig.1a,the ΔGH*of Ni/MoO2@CN (0.056 eV) on the Mo-ortho-C between Ni and MoO2interface (named as Mo-C1) is close to zero,which is superior to other active sites in Ni/MoO2@CN model (Figs.1a and S3).Furthermore,the ΔGH*value of Mo-C1 in Ni/MoO2@CN model is better than the ortho-C of MoO2@CN,Ni@CN,and CN models (Figs.S4-S9),indicating that the three-phase heterojunction interface can effectively adjust the H adsorption energy for ortho-C to obtain the highest HER intrinsic activity.

Afterward,the DOS of CN,MoO2@CN,and Ni/MoO2@CN models are used to further confirm the electronic structure of C at the three-phase heterojunction interface.As displayed in Fig.1b,the strength of this ordinate for CN,MoO2@CN,and Ni/MoO2@CN models indicates that the introduction of Ni species can construct the three-phase heterojunction interface to adjust the electronic structure of CN,which optimizes the ΔGH*of HER.Besides,for revealing the bond capabilities,the crystal orbital Hamilton population (COHP) derived from partial DOS (PDOS) calculation (Fig.1c,d) are implemented [36].The energy integral of COHP (ICOHP) of C-H for Ni/MoO2@CN is -3.38(Fig.1c),which is larger than those in other models (Figs.S5,S7,and S9),suggesting that it has stronger hydrogen adsorption to obtain higher HER catalytic activity [37].

Fig.1 a ΔGH* calculated at different adsorb sites for Ni/MoO2@CN model,insert:schematic illustration of H adsorption;b PDOS analysis of C for all samples;c COHP;and d PDOS analysis for the Ni/MoO2@CN model with the H atom adsorbed on the sites

Subsequently,Fig.2a shows the four-step mechanism processes for OER in alkaline media [38].The surface of Ni/MoO2@CN (Fig.2a),Ni/MoO2(Fig.S10),Ni@CN (Fig.S11),MoO2@CN (Fig.S12),and CN (Fig.S13) models displays the adsorption of OER intermediates.During the four OER reaction processes,the ΔGis employed to investigate the intrinsic activity,and the largest ΔGvalue indicates the rate-determining step (RDS) [39].As shown in Fig.2b,step 4 is RDS for Ni/MoO2,which exhibits ΔGvalue of 2.39 eV.For Ni/MoO2@CN,we combine CN with Ni/MoO2to form the three-phase heterogeneous junction interface,the RDS is changed from step 4 to step 2,and its ΔGvalue is reduced to 1.65 eV.The smaller ΔGvalue of Ni/MoO2@CN implies the faster reaction kinetics to obtain lower potential during the OER process,and it is better than other catalysts in this work (Figs.S14-S16).Besides,it is found that the Nid-band center of Ni/MoO2@CN exhibits a negative shift compared with Ni/MoO2(Fig.2c),which is also good for adjusting the adsorption of O-intermediates.Meanwhile,the PDOS results (Fig.2c) demonstrate that the electrons are transferred from Ni to CN,indicating that the three-phase heterojunction formed by CN and Ni/MoO2exists electron interaction,which could enhance the intrinsic activity for OER.

Fig.2 a Schematic illustration of the *OH,*O,and *OOH adsorption for Ni/MoO2@CN model;b OER reaction pathway for Ni/MoO2@CN model;c PDOS analysis of Ni for all samples;d COHP;and e PDOS analysis for the Ni/MoO2@CN model with the O atom adsorbed on the sites

The COHP (Fig.2d) is used to uncover the bond capabilities,which is originated from PDOS calculation (Fig.2e)[40].When CN combines with Ni/MoO2to form the threephase heterojunction,the strength of interaction between Ni and O can be changed.As shown in Figs.2d and S17a,the ICOHP (-0.69,-0.50,-0.56) values of Ni/MoO2@CN for Ni suggest that it possesses better interaction between Ni and O,in relative to ICOHP values (-0.54,-0.54,-0.54)of Ni/MoO2.Furthermore,the Mo atom has good oxygen adsorption capacity (-2.58) at the three-phase heterojunction interface to provide more O-intermediates to boost the OER intrinsic activity.Thus,the COHP results indicate that the three-phase heterojunction interface composed of Ni,MoO2,and CN can effectively optimize the desorption of O-containing intermediates to enhance the OER intrinsic activity.

2.2 Physicochemical Characterization

Guided by the DFT predictions,the facile two-step methods (Fig.3,the detailed steps are shown in Supplementary Information) are used to synthesize the CN layer coated Ni/MoO2nano-needle with three-phase heterojunction (Ni/MoO2@CN) self-supported on 3D porous NF.The powder used for X-ray diffraction (XRD) characterization is obtained from NF that has been treated by ultrasonication rather than directly using NF samples.As shown in Supplementary Figure S18,firstly,the NiMoO4nano-needle is grown on NF via the solvothermal method,which can be determined by XRD (Supplementary Figure S19).Then,the NiMoO4nano-needle is reduced by the high temperature(450 °C) calcination in a mixed atmosphere (5% H2+95%Ar) to form the Ni nanoparticle and MoO2nano-needle,and this is confirmed by XRD and scanning electron microscopy(SEM,Fig.S20-S21).The reason can be attributed to the different enthalpy of Ni and Mo [41].Besides,Raman proves the existing of CN,which can be ascribed to the reduction of organic carbon to CN [the ratio of area D and G (AreaDand AreaG) is 1.20],this is because the catalysis of Ni atoms(Fig.S22c) [22,42].Thus,these results suggest that the catalyst is composed of CN,Ni,and MoO2.

Fig.3 Synthesis diagram of Ni/MoO2@CN

Subsequently,the XRD,Raman,and SEM images obtained by annealing the NiMoO4nano-needle at 350 and 550 °C are shown in Figs.S22-S23,respectively.When it calcined at 350 °C,the characteristic peak intensities of Ni and MoO2are too weak.In this case,the morphology of nano-needle is smooth,and the ratio of AreaD/AreaGis 1.10 with less defects of CN.When it calcined at 550 °C,the characteristic peaks of Ni and MoO2are replaced by the MoNi alloy,and the morphology of nano-needle is destroyed,under the circumstances,the ratio of AreaD/AreaGis 1.15.Therefore,these results indicate that the temperature plays an important role on the formation of this catalyst.

Then,summarizing the results of XRD,Raman,and SEM(Figs.S22-S23) obtained at different temperatures,the following reason for the formation of catalyst can be proposed:The different enthalpy of Ni and Mo can lead to the segregation of Ni atoms from NiMoO4to form Ni as well as the reduction of MoO42-to MoO2at 450 °C,and the segregated Ni can catalyze the organic carbon to form CN.However,when the NiMoO4nano-needle annealed at 350 °C,due to the low temperature,there is no enough energy to reduce NiMoO4forming Ni and MoO2,and the surface organic carbon cannot form CN.And when it calcined at 550 °C,due to the high temperature,the Ni and Mo atoms segregate from the NiMoO4to form the NiMo alloy,but organic carbon cannot form the CN on its surface,mainly attributed to the rapid reduction at such high temperature.Therefore,the NiMoO4calcined at 350 and 550 °C exhibits lower HER/OER intrinsic activity (Figs.S29 and S41).

Subsequently,the transmission electron microscopy(TEM) and high-resolution TEM (HRTEM) are employed to reveal the three-phase heterojunction and lattice fringe of Ni/MoO2@CN.As shown in Fig.4a-c,the 0.208 nm (111)of Ni nanoparticle couples with the 0.242 nm (200) of MoO2nano-needle (consistent with XRD results).Obviously,the carbon layers are coated on the surface of Ni (≈ 3 layers,Fig.4b) and MoO2(≈ 2 layers,Fig.4d),which can efficiently optimize the electron distribution to enhance the HER/OER intrinsic activity,and this is also confirmed by DFT [43,44].Furthermore,carbon layers can avoid metal corrosion in electrolyte to improve the electrochemical stability [27].More importantly,a clear heterojunction exists between the Ni nanoparticle and MoO2nano-needle (Fig.4c,d),which can realize the redistribution of electrons at the interface for further optimizing the adsorption energy of H-and O-containing intermediates to boost the HER/OER intrinsic activity [25].Additionally,heterojunction can enhance the electron transfer ability to improve the WS performance at large current density [45].The high-angle annular darkfield scanning TEM (HAADF-STEM) and elemental mapping images sugggest the uniform distribution of C,N,O,Ni,and Mo on the Ni/MoO2@CN (Fig.4e-j).

X-ray photoelectron spectroscopy (XPS) is used to evaluate the elemental compositions and valance states.XPS survey spectrum of Ni/MoO2@CN is shown in Fig.S24a,which confirms the existence of Ni,C,Mo,N,and O elements.In Fig.S24b,c,HRXPS of O 1sand C 1sis located at Mo-O (530.8 eV),absorbed-O (532.1 eV),C=O(288.8 eV),C-N (286.3 eV),and C-C (284.8 eV).For N 1sspectrum,the peaks at 396.1 and 398.6 eV belong to Mo 3pand pyridinic-N (Fig.S24d),as known,the doping of N atoms into carbon layer can enhance the HER/OER intrinsic activity [46,47].The Ni 2pand Mo 3dof Ni/MoO2@CN show a positive shift (≈ 0.3 eV,Fig.4k and Table S1) and negative shift (≈ 0.3 eV,Fig.4l and Table S2) relative to Ni/MoO2,respectively,displaying that the CN affects the electronic structure of Ni and MoO2,and this is in agreement with the DFT predictions.

Furthermore,Ni 2pspectrum of Ni/MoO2@CN presents a positive shift (≈ 0.5 eV,Fig.4m and Table S1) relative to Ni@CN,showing that the introduction of MoO2can adjust the electronic structure of Ni@CN.In Fig.4n and Table S2,Mo 3dspectrum of Ni/MoO2@CN exhibits a negative shift (≈ 0.5 eV) as against to MoO2@CN,which suggests that the electronic structure of Mo 3dfor MoO2@CN is also affected,because it is coupled with Ni@CN to form the three-phase heterojunction interface.In contrast,as displayed in Fig.S25a,b,Ni 2p,and Mo 3dpeaks of Ni@CN+MoO2@CN hybrids (mixing the Ni@CN and MoO2@CN together) are the same as Ni@CN and MoO2@CN,respectively,indicating that the mechanical mixing cannot affect the binding energy.However,compared with the Ni@CN+MoO2@CN hybrids,the Ni 2pand Mo 3dpeaks of Ni/MoO2@CN present positive and negative shift,respectively.These results demonstrate that coupling Ni@CN with MoO2@CN to form the three-phase heterojunction interface can lead to the strong electronic interaction,thus adjusting the absorption/desorption of H/O-intermediates to enhance the intrinsic activity [48],which is consistent with DFT calculation results.

Fig.4 a-j TEM,HRTEM,HAADF-STEM,and elemental mapping images of Ni/MoO2@CN;k-n HRXPS of Ni 2p and Mo 3d spectra for Ni/MoO2@CN,Ni/MoO2,Ni@CN,and MoO2@CN

It can be seen from the above results that there is indeed a three-phase heterojunction between Ni,MoO2,and CN.This unique three-phase heterojunction can make the interfacial electronic redistribution of Ni,MoO2,and CN,thus providing lots of defects to form abundant active sites for reducing the WS reaction energy barrier.Furthermore,the nano-needle self-supporting on 3D NF can expose more active sites to further increase WS catalytic activity,and it is also beneficial for mass diffusion at large current density to improve the WS performance.

2.3 HER Catalytic Performance

Based on the guidance of DFT results and the unique threephase heterojunction of Ni/MoO2@CN,HER mechanism can be illustrated in Fig.5a,and the HER activity is evaluated by the linear sweep voltammetry (LSV) curves with/withoutiR-correction (Fig.S26).As displayed in Fig.5b,c,the Ni/MoO2@CN possesses good HER activity with low overpotentials (η-10=33 mV,η-1000=267 mV),which is better than that of MoO2@CN (η-10=134 mV,η-1000=430 mV).This result suggests that the introduction of Ni into MoO2@CN to form the three-phase heterojunction of Ni/MoO2@CN can greatly enhance the HER activity,and it is consistent with the DFT results.It is also better than that of precursor (η-10=210 mV,η-1000=570 mV),NF (η-10=256 mV,η-1000=590 mV),and close to Pt/C(η-10=25 mV,η-1000=296 mV),and superior to most of the recently reported catalysts (Fig.S27 and Table S3).Furthermore,Ni/MoO2@CN shows smaller overpotentials compared with Ni@CN (η-10=173 mV,η-1000=474 mV)in Figs.5c and S28a.This proves that the three-phase heterojunction of Ni/MoO2@CN synthesized by combining MoO2with Ni@CN leads to further optimize the electronic structure of CN for obtaining an outstanding HER intrinsic activity,and it agrees with DFT results.Meanwhile,Ni/MoO2needs higher overpotentials (η-10=100 mV,η-1000=358 mV) compared with Ni/MoO2@CN in Figs.5c and S28a,implying that the CN decorated Ni/MoO2can reduce the HER reaction barrier by optimizing the ΔGH*.

Fig.5 a Schematic diagram of HER mechanism:Ni/MoO2@CN exhibits better HER kinetics compared to MoO2@CN;b LSV curves;c comparison of overpotentials;d Tafel plots;e TOF curves;f ECSA-normalized LSV curves of as-prepared samples;g CP curve of Ni/MoO2@CN

Through above results comparing of Ni/MoO2@CN,Ni/MoO2,Ni@CN,and Ni/MoO2,we can conclude that the interfacial charge redistribution at three-phase heterojunction is good for HER,and the Ni,MoO2,and CN are all indispensable.Moreover,the sample with a Ni/Mo molar ratio of 1:7 (Fig.S29) and the NiMoO4annealed at 450 °C(Supplementary Figure S30) exhibits higher HER intrinsic activity,which indicates that the Ni/Mo molar ratio as well as the calcination temperature has an important impact on the HER intrinsic activity.

Figures 5d and S28b show Tafel plots originated from the LSV curves.The Ni/MoO2@CN exhibits a lower Tafel slope of 45 mV dec-1than those of MoO2@CN (92 mV dec-1),Ni@CN (85 mV dec-1),Ni/MoO2(81 mV dec-1),precursor (108 mV dec-1),NF (165 mV dec-1),and like Pt/C(40 mV dec-1),indicating that Ni/MoO2@CN obtains faster HER catalytic kinetics.Then,the charge transfer kinetics is analyzed by electrochemical impedance spectroscopy(Fig.S31),and Ni/MoO2@CN has the smallest charge transfer resistance (Rct) among all samples,implying the fastest charge transfer rate of HER.The superiority of HER reaction kinetics and charge transfer can be ascribed to interface electronic transmission at the three-phase heterojunction among CN,Ni,and MoO2.

Besides,the turnover frequency (TOF) and electrochemical surface area (ECSA) are used to explore the HER intrinsic activity and the number of active sites,respectively [49].Cyclic voltammetry (CV) method is used to study the TOF of HER (Fig.S32) [50].In Fig.5e and Table S4,Ni/MoO2@CN possesses the largest TOF value (1.45 s-1) at the overpotential of 100 mV,as compared with MoO2@CN (0.38 s-1),Ni@CN (0.19 s-1),and Ni/MoO2(0.53 s-1),suggesting that Ni/MoO2@CN has the highest catalytic activity.Besides,it also has the highest TOF values at other overpotentials(Table S5),and it is better than most of the previously reported results in literature (Fig.S33 and Table S6).

ECSA can be evaluated by electrochemical double-layer capacitances (Cdl),due to the positive correlation between ECSA andCdl(Fig.S34a-d) [42].Obviously,Ni/MoO2@CN exhibits the largestCdl(97.94 mF) in Fig.S34e,which implies that it has the highest density of active sites to accelerate the HER process,and it is also reflected in the ECSAnormalized LSV curves (Fig.5f).The highest intrinsic activity can be due to the interfacial charge redistribution at the three-phase heterojunction interface,which obtains suitable adsorption energy of H and O-containing intermediates.Besides,the self-supporting nano-needle can provide a large specific surface area to expose more active sites.

Stability is another vital index to evaluate catalyst performance [51],especially at large current density.Ni/MoO2@CN displays outstanding durability at -1000 mA cm-2for 200 h with stable release of H2bubbles,which is tested by the chronopotentiometry (CP) method (Fig.5g).Furthermore,the activity before and after the CP test has negligible changes (Fig.S35),implying its good stability.Subsequently,the SEM images (Fig.S36) and HRXPS spectra (Fig.S37)of Ni/MoO2@CN after the CP test are also obtained,and they almost maintain the original state.The better stability can be attributed to the following reasons:(1) The interfacial charge redistribution at three-phase heterojunction for Ni/MoO2@CN can accelerate the electron transfer at large current density;(2) Carbon layers could avoid the dissolution of metal in harsh solution;(3) Self-supporting nano-needle enhances the wettability of catalyst to effectively boost the mass diffusion and bubbles release.

2.4 OER Catalytic Performance

The Ni/NiMoO2@CN not only has outstanding HER performance,but also exhibits good OER performance under the same condition;the schematic diagram of the OER mechanism is displayed in Fig.6a.The OER LSV curves of Ni/MoO2@CN with/withoutiR-correction are presented in Fig.S38.As shown in Fig.6b,c,Ni/MoO2@CN shows higher activity (η10=250 mV,η1000=420 mV) than Ni/MoO2(η10=300 mV,η1000=550 mV),which indicates that the decoration of CN to Ni/NiMoO2can effectively enhance the OER intrinsic activity.Meanwhile,the overpotentials of Ni/MoO2@CN are also smaller than that of MoO2@CN (η10=310 mV,η1000=590 mV) and Ni@CN(η10=323 mV,η1000=660 mV) in Figs.6c and S39a.These results indicate that the electronic structures of Ni,Mo,and CN are adjusted by interfacial electronic redistribution at the three-phase heterojunction to further optimize the adsorption of O-intermediates;this is also consistent with the DFT results.Besides,as shown in Fig.6b,c,the Ni/MoO2@CN exhibits smaller overpotentials than NF (η10=345 mV,η1000=700 mV),precursor (η10=327 mV,η1000=650 mV),IrO2/C (η10=280 mV,η1000=540 mV),and even superior to most of the recently reported catalysts (Fig.S40 and Table S7).Moreover,the sample with a Ni/Mo molar ratio of 1:7 (Fig.S41) and the NiMoO4annealed at 450 °C (Fig.S42) displays the lowest OER overpotential,suggesting that the Ni/Mo molar ratio and calcination temperature play an important role on the OER intrinsic activity.

In addition,the Tafel slope and EIS are used to assess the kinetics of OER (Figs.6d,S39b and S43).Obviously,Ni/MoO2@CN displays the smallest Tafel slope and smallestRctvalue among all catalysts,which implies its high electron transfer rate and fast kinetics for OER at the threephase heterojunction interface.Subsequently,the good OER intrinsic activity of Ni/MoO2@CN is also verified by TOF.The active surface redox sites method is used to study the TOF for OER by calculating the redox surface sites of Ni2+/Ni3+without the capacitive current (Fig.S44)[50].In Fig.6e and Table S8,the TOF value of Ni/MoO2@CN (1.23 s-1) is higher than Ni/MoO2(0.28 s-1) and Ni@CN (0.14 s-1) at the overpotential of 300 mV,suggesting its higher intrinsic activity.Besides,it also has higher TOF values at other overpotentials (Table S9),which is better than most of the reported works (Fig.S45 and Table S10).The current density is normalized by ECSA to further reflect the highest intrinsic activity of Ni/MoO2@CN (Fig.6f).The remarkable OER catalytic activity can be ascribed to the following reasons:(1) Interfacial electronic redistribution at three-phase heterojunction interface can be beneficial to optimize the adsorption energy of O-containing intermediates;(2) Interfacial electronic redistribution also can accelerate the charge transfer at large current density to enhance OER activity;(3) Self-supporting nano-needle has a large specific surface area to provide more active sites for boosting the OER performance.

Fig.6 a Schematic diagram of OER mechanism:Ni/MoO2@CN exhibits better OER kinetics compared to Ni/MoO2;b LSV curves;c comparison of overpotentials;d Tafel plots;e TOF curves;f ECSA-normalized LSV curves of as-prepared samples;g CP curve of Ni/MoO2@CN

Subsequently,its stability is studied by CP method.Impressively,the Ni/MoO2@CN maintains good stability for 200 h at 1000 mA cm-2(Fig.6g);the change of potential is about 70 mV.To investigate its outstanding durability,LSV curves and Nyquist plots before and after stability test are displayed in Fig.S46;the potential at 1000 mA cm-2andRctvalue at 1.5 V have minor changes.The reason for the good OER stability should be the same as HER.Then,the change of morphology and elements states after the CP test are studied by SEM and HRXPS.In Fig.S47,SEM images of Ni/MoO2@CN can keep the primary appearance,indicating its good stability.The HRXPS spectrum of Ni0disappeared (Fig.S48),implying that the Ni0is oxidized to NiOOH [52],and the Mo4+and Mo5+are also oxidized to Mo6+,which is a common phenomenon during OER process [42].As reported in literatures [52-54],the generated NiOOH can combine with the Ni to form the heterojunction interface between the NiOOH and Ni during OER process,thus further enhancing the OER intrinsic activity.Meanwhile,according to literatures,the generated Mo6+can be beneficial to the conversion of Ni2+to Ni3+and produce the NiOOH for enhancing the OER intrinsic activity [55,56].

2.5 WS Catalytic Performance

In considering of the outstanding HER and OER performance of Ni/MoO2@CN,the WS performance of Ni/MoO2@CN is also evaluated by two-electrode system to simulate the actual application,and the schematic diagram is shown in Fig.7a.For comparison,the WS intrinsic activity of Pt/C//IrO2/C and NF//NF are also measured by the same method with Ni/MoO2@CN.Under 1.0 M KOH+30 °C condition,the LSV curve of Ni/MoO2@CN (1.83 V) as cathode and anode electrodes display a higher activity than Pt/C//IrO2/C (1.89 V) and NF//NF (2.26 V) at 200 mA cm-2(Fig.7b).Interestingly,the Ni/MoO2@CN as bifunctional electrodes can drive the large current density of 500 and 1000 mA cm-2while only requiring 1.92 and 2.02 V,and it is better than most of the reported literatures at 500 mA cm-2(Fig.7c) [57-65].Furthermore,it can maintain for 300 h at 1000 mA cm-2(Fig.7d),exhibiting its excellent stability.

In addition,in order to further evaluate the industrial application potential of Ni/MoO2@CN,it is applied to WS under the simulated industrial condition (6.0 M KOH+60 °C).As shown in Fig.7b,it only requires 1.86 V for Ni/MoO2@CN to reach 1000 mA cm-2when used as cathode and anode electrodes.It can work stably for 330 h at 1000 mA cm-2with negligible potential change (Fig.7e),implying that it is a potential material for large-scale hydrogen production.The good WS performance of Ni/MoO2@CN under the large current density of industrial condition can be attributed to the three-phase heterojunction with high HER/OER intrinsic activity and the self-supporting nanoneedle structure with more active sites,faster bubbles desorption,and the diffusion of electrolyte.

Subsequently,the Faradic efficiency and contact angles(CA) of Ni/MoO2@CN are surveyed (Figs.7f,g and S49).The volume-time curve for H2and O2is obtained by water drainage method at 0,10,15,20,25,30,35,40,45,and 50 min measured at 30.0 mA.The volume-time curve for H2and O2is about 2:1,and they are close to the theoretical values,indicating the nearly 100% Faradic efficiency.Besides,the wettability of catalysts is very important for WS.The CA of Ni/MoO2@CN and NF with 1.0 M KOH solution are shown in Fig.7g.Obviously,the CA of Ni/MoO2@CN (9°)is smaller than NF (45°),suggesting the good gas release and mass transfer.Then,as shown in Fig.S50,the size of H2bubbles on NF electrode is much larger than it on Ni/MoO2@CN electrode,indicating that the H2bubbles are easier to desorb from the surface of Ni/MoO2@CN.Thus,the above results suggest that Ni/MoO2@CN has good hydrophilicity to enhance the contact with electrolyte and facilitate bubbles desorption for improving the performance of WS.

Fig.7 a Schematic diagram of WS for Ni/MoO2@CN;b LSV curves of the samples in 1.0 M KOH+30 °C,and LSV curve of Ni/MoO2@CN in 6.0 M KOH+60 °C;c comparison of potentials at 500 mA cm-2 of Ni/MoO2@CN and reported catalysts;CP curves of Ni/MoO2@CN d at 1.0 M KOH+30 °C and e 6.0 M KOH+60 °C;f Volume of H2 and O2 measured at 30.0 mA as a function of time for Ni/MoO2@CN;g CA of NF and Ni/MoO2@CN under 1.0 M KOH.Insert:optical photographs

3 Conclusions

In short,the CN layers encapsulated Ni/MoO2nano-needle with three-phase heterojunction is successfully engineered and synthesized.Interestingly,as predicted by DFT results,the electrons are redistributed among the CN,Ni,and MoO2for optimizing the adsorption energy of H-and O-containing intermediates in order to obtain the best ΔGH*for HER and decrease the ΔGvalue of RDS for OER,thus enhancing the HER/OER intrinsic activity.The physical and electrochemical results are in good agreement with the theoretical predictions.Compared with the XPS results of Ni/MoO2,Ni@CN,and MoO2@CN,the peaks of Ni 2pand Mo 3dfor Ni/MoO2@CN show negative/positive shift,which further demonstrate that there is a strong interaction between CN,Ni,and MoO2.Therefore,the Ni/MoO2@CN exhibits good activity for HER (?-10=33 mV) and OER (?10=250 mV) compared with Ni/MoO2(HER:?-10=100 mV,OER:?10=300 mV),MoO2@CN (HER:?-10=134 mV,OER:?10=310 mV),and Ni@CN (HER:?-10=173 mV,OER:?10=323 mV).Meanwhile,the three-phase heterojunction can facilitate electronic transfer at large current density and the carbon layers can avoid the metal dissolution in harsh solution.Self-supporting nanoneedle can enhance the wettability of catalyst to effectively boost the mass diffusion and bubbles release.Therefore,it can steadily operate for 330 h at 1000 mA cm-2for WS (6.0 M KOH+60 °C),suggesting an excellent stability.This work offers a unique idea of designing catalytic materials for industrial electrochemical water splitting.

AcknowledgementsThis work is supported by the National Natural Science Foundation of China (21872040,22162004),the Hundred Talents Program of Guangxi Universities,the Excellent Scholars and Innovation Team of Guangxi Universities,the Innovation Project of Guangxi Graduate Education (YCBZ2021011),and the High-performance Computing Platform of Guangxi University.

FundingOpen access funding provided by Shanghai Jiao Tong University.

Open AccessThis 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 licence,and indicate if changes were made.The images or other third party material in this article are included in the article’s Creative Commons licence,unless indicated otherwise in a credit line to the material.If material is not included in the article’s Creative Commons licence 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 licence,visit http:// creat iveco mmons.org/ licen ses/ by/4.0/.

Supplementary InformationThe online version contains supplementary material available at https:// doi.org/ 10.1007/s40820-021-00744-x.