B-MoSe2 3D nanoshells were obtained by a hard template method. Au nanospheres were first synthesized by reducing gold precursors in oil bath system with ethylene glycol (EG) as both reductant and solvent (Supplementary Fig. 1). A given amount of Au nanospheres/graphene (G) mixture underwent sequentially annealing treatment by Mo(CO)6 and Se powder, forming Au@MoSe2/G composite material (Supplementary Figs. 2–4). After removing the Au core, B-MoSe2 were synthesized in high yield of approaching 100%. Transmission electron microscopy (TEM) images emphasize that the product is characterized by monodispersed hollow nanoshells with an average diameter of 56.8 ± 2.5 nm (Fig. 2a and Supplementary Figs. 5–7). High-resolution transmission electron microscopy (HRTEM) image reveals the uniform few-atomic thick layers and the atomically curved MoSe2 from the cross-sectional view (Fig. 2b, c). Aberration-corrected high-angle annular dark-field scanning transmission electron microscopy (AC-HAADF-STEM) confirms the existence of 2H-MoSe2, which resolves a characteristic hexagonal lattice arrangement of Mo (brighter dots) and Se (darker dots) atoms (Fig. 2d and Supplementary Fig. 8). The Fast Fourier transform (FFT) image indicates the corresponding (100) and (110) planes of 2H phase. Moreover, the biaxial strain endows MoSe2 with abundant Se vacancies evidenced by the yellow circle. A closer inspection of the cross-sectional view of the nanoshell suggests that the typical atomic arrangement of Mo and Se in (002) plane of 2H-MoSe2 with expanded lattice spacing of 0.69 nm. To investigate their unique atomic configuration, the geometric phase analysis (GPA) is adopted to detect the lattice strain changes (Fig. 2f). The obvious evenly color can be observed in the strain mapping, indicating the atomic curvature generates a uniform biaxial strain of ~1.2%. HAADF-STEM and the corresponding EDS mapping images show the uniform distribution of Mo and Se elements (Fig. 2g). For detailed comparison, we further fabricated uniaxially strained (U-MoSe2) and flat (F-MoSe2) MoSe2 without strain using similar strategies (Fig. 2h–k and Supplementary Figs. 9–12). As shown in AC-HAADF-STEM image and corresponding scheme, the U-MoSe2 with an average diameter of 57.1 ± 2.8 nm is composed of curved (002) planes with a lattice spacing of 0.69 nm (Fig. 2i and Supplementary Fig. 7). The uniaxial strain was calculated to be ~1.2%, which is similar with the biaxial one. This enables us to investigate the effect of strain dimension on the HER behaviors. The F-MoSe2 present a flat (002) planes with decreased lattice spacing (0.65 nm). Overall, all these results conclude that MoSe2 with controlled strain dimensions has been successfully prepared.Fig. 2: Morphology characterization of B-MoSe2, U-MoSe2, and F-MoSe2.a–c TEM images under different magnifications. The inset in c is the scheme of B-MoSe2. d Aberration-corrected HAADF-STEM image of B-MoSe2. The inset is the FFT image. e Aberration-corrected HAADF-STEM image from cross-sectional view. f GPA analysis for (d). The change from blue to red in color scale indicates the gradually change from compressive to tensile strain. g STEM and corresponding EDS elemental mapping images. TEM, aberration-corrected HAADF-STEM image, and scheme of (h, i) U-MoSe2 and (j, k) F-MoSe2. The bule and yellow spheres in (d, e, i, k) indicate Mo and Se atoms, respectively. The yellow dash circle is Se vacancy.The electronic structure and chemical environment of MoSe2 with controlled strain dimensions were investigated by Raman spectrum, X-ray photoelectron spectroscopy (XPS), X-ray absorption spectroscopy (XAS), and electron paramagnetic resonance (EPR). The Raman spectrum confirms the 2H phase of MoSe2 of samples, indicating the strain does not change the crystal phase (Fig. 3a). It is known that the A1g vibration modes are sensitive to the structural information of MoSe2. Compared to the flat MoSe2, the A1g mode shows an obvious blueshift with the increase of strain dimension, which is usually generated from the increased vibration energy27. The high-resolution XPS Mo 3d signals show the double peaks centered at ~228.8 and ~232.2 eV can be assigned to 3d5/2 (Mo4+) and 3d3/2 (Mo4+), respectively (Fig. 3b and Supplementary Fig. 13)28. With increasing the strain dimension from F-MoSe2, U-MoSe2 to B-MoSe2, the binding energy of 3d5/2 gradually blueshifts from 228.8, 228.9 to 229.1 eV, indicating a more delocalized local electric structure. Another peak at about 235.8 eV for F-MoSe2 is attributed to Mo6+ species from the surface oxidation. Moreover, the Se 3d5/2 peaks of B-MoSe2 shift to lower binding energy (54.33 eV) when compared to those of F-MoSe2 (54.54 eV) and U-MoSe2 (54.47 eV), as shown in Fig. 3c. This means Se sites are pumped into abundant electronics by the effect of biaxial strain, which is favorable for the adsorption of proton and thus accelerating HER performance29,30. Moreover, the absence of oxidation peaks in U-MoSe2 and B-MoSe2 suggests their superior anti-oxidation ability and structure stability compared to the F-MoSe2 (Fig. 3b). Possible reason is, an enhanced electron interaction is induced in the strained MoSe2 compared to the flat counterpart, leading to a stronger Mo-Se bonding in the former.Fig. 3: Chemical state and electronic structure analysis.a Raman spectra, (b) Mo 3d and (c) Se 3d XPS spectra. d Mo K-edge XANES spectra. e The fitted average oxidation states and d-band hole counts of Mo from XANES spectra. f The corresponding FT-EXAFS k3χ(R) spectra. g EXAFS wavelet transform analysis. The color scale changed from blue to red indicate the increasing intensity of Mo-Se bonds. h Comparison of the mass normalized EPR spectra.X-ray absorption near-edge spectroscopy (XANES) was conducted for the MoSe2 composites at the Mo K-edge to study intrinsic electronic structures. As shown in Fig. 3d, all the synthesized MoSe2 samples show the absorption edges between the MoSe2 and MoO3 references, implying that their average valence of Mo species between +4 and +6 coincide with the XPS result. The peaks at ~20,010 eV are mainly due to the electron transition between 1s and 4d orbital in Mo-anion configuration ([MoO6] in MoO3 and [MoSe6] in MoSe2)31,32. With increasing the strain dimension from F-MoSe2, U-MoSe2 to B-MoSe2, this peak gradually grows, indicating the enhanced electron interaction between Mo and Se. Moreover, the adsorption edge of B-MoSe2 shifts to higher energy compared to F-MoSe2 and U-MoSe2, which is due to the strong electronic delocalization by biaxial strain effect. The white-line intensity of B-MoSe2 is higher than that of F-MoSe2 and U-MoSe2, and the oxidation states can be estimated to be +4.17 for F-MoSe2, +4.34 for U-MoSe2, and +4.45 for B-MoSe2, using MoSe2 (4d2) and MoO3 (4d0) as references (Fig. 3e). Moreover, the d-band hole number for F-MoSe2, U-MoSe2, and B-MoSe2 was estimated to be 8.07, 8.34, and 8.45, respectively, indicating the variation of the unoccupied Mo 4d-orbitals.In addition, the extended X-ray adsorption fine structure (EXAFS) (Fig. 3f) and the corresponding wavelet transform (Fig. 3g) both show gradual decrease in intensity. As listed in Supplementary Table 1, the Mo-Se coordination numbers of F-MoSe2, U-MoSe2, and B-MoSe2 are 5.9, 4.4, and 4.1, respectively, which are lower than that (6) of MoSe2 foil (Supplementary Fig. 14). Such a trend is mainly due to the increased strain dimension of MoSe2, resulting in the breaking of Mo-Se bonds and thus creating more Se vacancies. Moreover, compared with F-MoSe2, there appears a significant shift of Mo-Mo bonds to longer location in R-space for U-MoSe2 and B-MoSe2, which can be attributed to the tensile strain on the surface and the enlarged Mo-Mo bond length. Se vacancies are important factors for the HER activity of strained MoSe2. Hence, electron paramagnetic resonance (EPR) was further performed to investigate the condition of Se vacancies. The EPR signal intensity at g = 2.009 of F-MoSe2 approaches 0, while those of U-MoSe2 and B-MoSe2 are drastically enhanced. Moreover, even the tensile strain degree is similar, B-MoSe2 exhibits the higher concentration of Se vacancies than U-MoSe2 (Fig. 3h). Based on the above fine structure analysis, the tensile strain effect could vary the oxidation state, d band holes, and vacancies of atoms as well as the Mo-Mo bond distance. At the same strain degree, this effect is more immense in case of biaxial strain.To investigate the stabilized OH engineer for accelerating alkaline hydrogen evolution, three types of MoSe2 were initially evaluated to explore the underlying influence of tensile strain dimension on the catalytic activity by a three-electrode cell in 1 M KOH aqueous electrolyte. All the polarization curves were calibrated against reversible hydrogen electrode (RHE) as shown in Supplementary Fig. 16. Linear sweep voltammograms (LSV) obtained for F-MoSe2, U-MoSe2, B-MoSe2, and commercial Pt/C (20 wt%) are presented in Fig. 4a. B-MoSe2 shows the highest HER catalytic activity, as evident from the early onset potential and lowest overpotentials of only 58.2 and 190.5 mV to deliver current densities of 10 and 100 mA cm−2, respectively, superior to F-MoSe2 (144.3 and 330.1 mV), and U-MoSe2 (101.1 and 231.2 mV) (Fig. 4b). The overpotential of B-MoSe2 at 10 mA cm−2 is much lower than that of the reported work (most of them are higher than 100 mV), indicating superior activity (Supplementary Table 2). The reaction kinetics were further investigated by Tafel analysis (Fig. 4c). B-MoSe2 demonstrates the fastest HER reaction kinetics confirmed by the smallest Tafel slope value of 97.1 mV dec−1, much lower than that of F-MoSe2 (131.8 mV dec−1) and U-MoSe2 (122.4 mV dec−1). Even three catalysts all show the Volmer-Heyrovsky mechanism according to their Tafel slope values, the HER kinetics exhibited by F-MoSe2 and U-MoSe2 are more sluggish compared to that of B-MoSe2, revealing favorable surface electronic modulation effects via the biaxial stain. The lower Tafel slope value of B-MoSe2 indicates a fast water dissociation (Volmer step) behavior (H2O + * + e− → H* + OH−). The decreased energy barrier of the RDS may have substantially accelerated the overall hydrogen evolution rates. The double-layer capacitance (Cdl) measurement was carried out to evaluate the electrochemically active surface areas (ECSAs) of catalysts (Fig. 4d and Supplementary Fig. 17). The B-MoSe2 catalyst displays the highest Cdl value of 34.3 mF cm−2 among the other two samples, implying the highest ECSA and more accessible active sites. Moreover, the catalytic performance was normalized by ECSA, which shows the best intrinsic activity of B-MoSe2 (Supplementary Fig. 18). To further evaluate the intrinsic activity, the turnover frequency (TOF) value for B-MoSe2 at the overpotential of 200 mV was calculated to be 8.7 s−1, which is an order of magnitude higher than that of the U-MoSe2 catalysts (0.79 s−1) and ~2 order of magnitude higher than that of F-MoSe2 (0.09 s−1), showing the most favorable HER kinetics. The electrochemical impedance spectroscopy (EIS) measurements were then performed to track the charge transfer kinetics between H2O and catalyst surface (Fig. 4e). As expected, all MoSe2 samples exhibit the similar uncompensated solution resistance (Rs) value of ~3.0 Ω. Furthermore, the charge-transfer resistance (Rct) of the B-MoSe2 electrode was measured to be 11.6 Ω in Nyquist plot, significantly lower than those of F-MoSe2 (45.2 Ω) and U-MoSe2 (30.5 Ω), suggesting the faster charge transfer kinetics between biaxial strained surface and H2O molecules. The adsorption behavior of OH− on the catalyst surface could govern the thermodynamics and kinetics of alkaline HER33. The CO-stripping measurement was thus applied to assess the electrocatalyst’s capability for adsorbing OH, owing to the CO-oxidation process ignites by the reactive OH* 6,34. It means that the lower the onset potential, the stronger the OH− adsorption capability. By measuring the LSV curves in CO-saturated electrolyte (Fig. 4f), one can find that compared with those on U-MoSe2 and F-MoSe2, the onset potential of CO oxidization on B-MoSe2 has negatively shifted ~40 and ~190 mV, respectively. These significant shifts caused by the strain effect indicate the easier adsorption of OH− on the curved surface of MoSe2 than flat one. More importantly, biaxial strain could generate much lower OH− adsorption energy compared to uniaxial strain. Based on the above comparison results, we have successfully verified that the biaxial strain from B-MoSe2 plays a more significant role in accelerating the kinetics of Volmer step and decreasing OH− adsorption energy in alkaline HER activity. Moreover, we investigated the effect of the magnitude of biaxial strain on HER activity by obtaining MoSe2 nanoshells with other diameters (35.7 ± 3.6 and 115.0 ± 6.0 nm) (see Method and Supplementary Figs. 19, 20). LSV curves show that decreasing the B-MoSe2 diameters results in higher HER activity and reaction kinetics (Supplementary Fig. 21). The B-MoSe2 with diameter of 35.7 ± 3.6 nm displays the highest Cdl compared to the other two samples (Fig. 4d and Supplementary Fig. 21). Hence, a higher magnitude of biaxial strain is beneficial to the HER activity.Fig. 4: Alkaline HER performance.a HER polarization curves of F-MoSe2, U-MoSe2, B-MoSe2, and 20 wt% Pt/C. The potential is corrected by an automatic 90% of iR compensation. b Comparison of overpotentials at 10 and 100 mA cm−2, c Tafel plots, (d) half of current density differences (Ja-Jc) plotted against scan rates, and (e) Nyquist plot of electrochemical impedance spectra. charge-transfer resistance (Rct) of the B-MoSe2, U-MoSe2, and F-MoSe2electrode was measured to be 11.6, 30.5, and 45.2 Ω. f LSV curves in CO-saturated 1 M KOH. g Chronopotentiometric test at two current densities. h Durability test of the MEA electrolyser at the current density of 1 A cm−2. The insets in (h) are scheme of MEA electrolyser (left) and LSV curves operated at 25 and 60 °C (right).The B-MoSe2 catalyst exhibits advanced stability in the alkaline HER. We evaluated the durability through the chronopotentiometric method up to 200 h, during which only a slight decrease in potential occurred at a set current density of 200 mA cm−2 (Fig. 4g). Even at a higher current density of 500 mA cm−2, it still maintained 94.3% after more than 70 h continuous test. After the durability test, the hollow morphology, composition, and electronic structure of B-MoSe2 were evaluated in detail. There are no apparent shifts for the Mo 3d and Se 3d XPS spectra compared to the initial ones (Supplementary Fig. 22). Inductively coupled plasma-mass spectrometry (ICP-MS) indicates the molar ratio of Mo/Se was still approaching 1/2, verifying the stable structures (Supplementary Table 3). In addition, the unchanged hollow structures were confirmed by TEM and HRTEM images (Supplementary Fig. 23). Upon closer examination, it can be observed that the partial surface of B-MoSe2 has become thinner by stripping away certain layers, which leads us to conclude the reason of superior durability. We further conducted the aberration-corrected HAADF-STEM image of a single B-MoSe2 after durability test (Supplementary Fig. 24). The peeling off area is clearly observed. The lattice strain was determined to be ~1.2% by GPA, which is comparable to that of the pristine sample (Supplementary Fig. 24d). In addition to the stronger Mo-Se bonds, the onion-like shell of B-MoSe2 ensures the remarkable stability during catalysis. Catalysis primarily occurs on the external surface of the catalyst. For layered B-MoSe2, it possesses a nearly perfect spherical hollow structure, and each layer can function as an active layer due to the uniformity of biaxial strain. The outermost layer of B-MoSe2 can be stripped away due to structural damage during catalysis. The exposed surface from the inner layers becomes new catalytic active sites. This self-protection process will boost the stability of active structures. To test B-MoSe2 under more practical conditions, we constructed a membrane-electrode-assembly (MEA) based electrolyser, using anion exchange membrane (AEM) as the solid electrolyte and commercial IrO2/Ti as the anodic catalyst (in set in Fig. 4h and Supplementary Fig. 25). A higher temperature could enhance the MEA performance than that at 25 °C. Specifically, cell voltages of only 2.13 and 2.37 V are needed to reach current densities of 1 and 1.5 A cm−2 at an operation temperature of 60 °C, respectively. When the B-MoSe2-based AEM electrolyser runs 100 h at 1 A cm−2, only slight activity loss was observed, better than Pt/C60 and other advanced catalysts (Fig. 4h and Supplementary Table 4)35.The behavior of interfacial water on electrocatalysts is extremely important for water adsorption/dissociation during alkaline HER. The Operando Raman spectra under different applied potentials was performed to reveal the effects of strain dimensions on the reaction of interfacial water (Fig. 5 and Supplementary Fig. 26). The broad band in the region of 3000–3750 cm−1 can be ascribed to the O-H stretching mode (νO–H) of H2O36. It can be further deconvoluted into three different peaks: four hydrogen bonded water (4HB-H2O), two hydrogen bonded water (2HB-H2O), and hydrated K+ ion water (K+-H2O) (Fig. 5a, b and Supplementary Fig. 27)37. As the applied potentials negatively increases from 0 to −0.5 V, the content of K+-H2O on the surface of B-MoSe2 decreases from 8.9 % to 2.1%, while that of 4HB-H2O increases from 45.7 % to 56.0 %. For the U-MoSe2 and F-MoSe2, the contents of 4HB-H2O and K+-H2O shows negligible changes (Fig. 5c). It can be deduced that the biaxial strain can effectively induce the transition of K+-H2O to 4HB-H2O through desolvation, thereby confirming the accelerated dissociation of H2O and the occurrence of water reduction38. The relationship of the vibrational frequency of adsorbate as a function of applied potentials has been considered to be vibrational Stark effect, which can describe the sensitivity between intermediates and catalyst surfaces12. As shown in Fig. 5d, the Stark tuning rate of 4HB-H2O for B-MoSe2 was evaluated to be 23.2 cm−1/V, which is much higher than those of U-MoSe2 (6.4 cm−1/V) and F-MoSe2 (5.9 cm−1/V). Moreover, the vibrational frequency at −0.3 V of 4HB-H2O on B-MoSe2 (3252.0 cm−1) is lower than that of U-MoSe2 (3258.3 cm−1) and F-MoSe2 (3257.7 cm−1). This result can be explained by the stronger electrostatic interaction and closer distance between the 4HB-H2O and the B-MoSe2 surfaces at negative potentials, compared to the other two surfaces. Such an strong interaction is conducive to the breaking of the O-H bond of surficial water and protons transfer36. Therefore, we may conclude that the biaxial strain can promote the transformation of K+-H2O to 4HB-H2O via desolvation at negative potentials. This mechanism and its effect to alkaline HER performance is schematically illustrated in Fig. 5e, f.Fig. 5: Behavior of interfacial water on different MoSe2.Operando Raman spectra of interfacial water at (a) B-MoSe2 and (b) F-MoSe2. The 4HB-H2O, 2HB-H2O, and K+-H2O present three kinds of O-H stretching modes for H2O. c Stark tuning rates of the νO-H in Raman spectra for 4HB-H2O. d Comparison of the contents for 4HB-H2O and K+-H2O during alkaline HER. Scheme of interfacial water on (e) F-MoSe2 and (f) B-MoSe2. 4HB-H2O, 2HB-H2O, and K+-H2O indicate four hydrogen bonded water, two hydrogen bonded water, and hydrated K+ ion water, respectively. The blue (e) and red (f) shadings present sluggish and enhanced reaction between water and catalysts.Density functional theory (DFT) calculations were conducted to deep insight to the reaction mechanism. Three strain dimensionalities of 2H-MoSe2 (without strain, with uniaxial strain and biaxial strain) were constructed as the theoretical models, in which the strain was applied to be 1.2% (Supplementary Fig. 28). Considering the practical condition, the Mo coordination numbers of 2H-MoSe2 were specified as 6 (6Mo), 5 (5Mo), 4 (4Mo), and 3 (3Mo) (Supplementary Fig. 29). The OH- species prefer to stabilize at the Se top site of MoSe2 (Supplementary Fig. 30). As the ignition step for HER, the H2O adsorption (\({\Delta {{{\rm{G}}}}}_{{{{{\rm{H}}}}}_{2}{{{\rm{O}}}}}\)), was first systematically investigated (Fig. 6a and Supplementary Fig. 31). The optimization of MoSe2 surface with HO- modification shows the enhanced water adsorption, which is consistent with the experimental results. Moreover, in case of the biaxially strained MoSe2 with the modification of OH- adsorption, the adsorption of H2O strengthens as the coordination number decreases. To be more specific, the \({\Delta {{{\rm{G}}}}}_{{{{{\rm{H}}}}}_{2}{{{\rm{O}}}}}\) of H2O on B-3Mo-OH is as low as −0.10 eV, indicating a spontaneous exothermic reaction (Fig. 6b). Notably, H2O prefers to adsorb onto OH* via the formation of hydrogen bond between OH* and H2O (Fig. 6c). The change of adsorption configuration of H2O from O-down on Mo (B-3Mo) to O-horizontal mode on OH* (B-3Mo-OH) may be responsible for the negative \({\Delta {{{\rm{G}}}}}_{{{{{\rm{H}}}}}_{2}{{{\rm{O}}}}}\).Fig. 6: Computation of the biaxial strain effect on water adsorption and hydrogen evolution.a Calculated \({\Delta {{{\rm{G}}}}}_{{{{{\rm{H}}}}}_{2}{{{\rm{O}}}}}\) values of various Mo coordination numbers. b Comparation of \({\Delta {{{\rm{G}}}}}_{{{{{\rm{H}}}}}_{2}{{{\rm{O}}}}}\) values between different Mo configurations and strain dimension with (with OH) and without OH− (without OH) adsorption. c Model of H2O adsorption on B-3Mo-OH. d Calculated \({\Delta {{{\rm{G}}}}}_{{{{{\rm{H}}}}}^{*}}\) for U-5Mo-OH, B-5Mo-OH, B-4Mo-OH, and B-3Mo-OH. e, f The H* adsorbed optimized structures of B-4Mo and B-4Mo-OH, respectively. g The charge density difference plot of B-4Mo-OH. h Calculated pDOS of Mo-3d band and H-1s band. i Schematic of the tandem process (coexistence of 3Mo and 4Mo) for alkaline HER on the B-MoSe2. For the name of sample, the number before Mo is the coordination number of Mo; U and B presents uniaxial and biaxial strain, respectively; the OH after Mo indicates the OH adsorption on MoSe2.As one of the crucial descriptors for HER activity, the \({\Delta {{{\rm{G}}}}}_{{{{{\rm{H}}}}}^{*}}\) values for series of MoSe2 configurations are further screened (Supplementary Fig. 33 and Table 5). One can find that strain, Se vacancies, and OH− modification will influence the hydrogen binding. The \({\Delta {{{\rm{G}}}}}_{{{{{\rm{H}}}}}^{*}}\) of B-4Mo-OH was calculated to be −0.004 eV, which is better than compared samples (Fig. 6d)39. As shown in Fig. 6e, the H species prefer to adsorb and stabilize on the hollow sites bonded with three Mo atoms in B-4Mo. When the OH- ions stabilize on Se sites in B-4Mo-OH, H will preferentially adsorb onto the unsaturated 4Mo sites closest to OH* with a reduced Mo-H bond length of 1.72 Å (Fig. 6f), indicating an enhanced binding energy of H. Additionally, the electronic structures of the Mo species connected to the OH* adsorbed Se sites all undergo alterations (Fig. 6g). The Se site transforms into electron donor to stabilize OH- adsorbates. More electrons are thus transferred from the bridging Mo sites to maintain the charge balance of B-4Mo-OH, leading to substantial changes for electronic structure (Supplementary Table 6). From the projected density of states (pDOS, Fig. 6h), the B-4Mo-OH structure exhibits an stronger hybridization of Mo-3d and H-1s orbitals, thus leading to the strengthened Mo-H bond.By integrating the results of water and hydrogen adsorption, the B-3Mo-OH structure exhibits remarkable water adsorption capability via O-horizontal configuration, while the B-4Mo-OH structure displays the highest intrinsic hydrogen evolution activity. Considering that both 4Mo and 3Mo structures should coexist in our B-MoSe2, a collaborative promotion mechanism may occur during HER (Fig. 6i). 3Mo sites may accelerate the water adsorption, and the generated H* then preferentially spillover and adsorb onto the 4Mo sites, expediting the hydrogen adsorption/desorption. This tandem process within the B-MoSe2 structure is thus expected to facilitate the catalytic kinetics for alkaline HER.
