# Electrocatalytic hydrogen evolution on the noble metal-free MoS 2 /carbon nanotube heterostructure: a theoretical study

Farhad Keivanimehr, Sajjad Habibzadeh, Alireza Baghban, Amin Esmaeili, Ahmad Mohaddespour, Amin Hamed Mashhadzadeh, Mohammad Reza Ganjali, Mohammad Reza Saeb, Vanessa Fierro, Alain Celzard

► To cite this version:

Farhad Keivanimehr, Sajjad Habibzadeh, Alireza Baghban, Amin Esmaeili, Ahmad Mohaddespour, et al.. Electrocatalytic hydrogen evolution on the noble metal-free MoS 2 /carbon nanotube heterostructure: a theoretical study. Scientific Reports, 2021, 11, pp.3958. 10.1038/s41598-021-83562-w. hal-03156955

HAL Id: hal-03156955

https://hal.univ-lorraine.fr/hal-03156955

Submitted on 2 Mar 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

# Electrocatalytic hydrogen evolution on the noble metal-free MoS₂/carbon nanotube heterostructure: a theoretical study

Farhad Keivanimehr<fnref n="1" />, Sajjad Habibzadeh<fnref n="1" />,²✉, Alireza Baghban<fnref n="1" />, Amin Esmaeili³, Ahmad Mohaddespour⁴, Amin Hamed Mashhadzadeh⁵, Mohammad Reza Ganjali⁵, Mohammad Reza Saeb⁵, Vanessa Fierro⁶ & Alain Celzard⁶

Molybdenum disulfide (MoS₂) is considered as a promising noble-metal-free electrocatalyst for the Hydrogen Evolution Reaction (HER). However, to effectively employ such material in the HER process, the corresponding electrocatalytic activity should be comparable or even higher than that of Pt-based materials. Thus, efforts in structural design of MoS₂ electrocatalyst should be taken to enhance the respective physico-chemical properties, particularly, the electronic properties. Indeed, no report has yet appeared about the possibility of an HER electrocatalytic association between the MoS₂ and carbon nanotubes (CNT). Hence, this paper investigates the synergistic electrocatalytic activity of MoS₂/CNT heterostructure for HER by Density Functional Theory simulations. The characteristics of the heterostructure, including density of states, binding energies, charge transfer, bandgap structure and minimum-energy path for the HER process were discussed. It was found that regardless of its configuration, CNT is bound to MoS₂ with an atomic interlayer gap of 3.37 Å and binding energy of 0.467 eV per carbon atom, suggesting a weak interaction between CNT and MoS₂. In addition, the energy barrier of HER process was calculated lower in MoS₂/CNT, 0.024 eV, than in the MoS₂ monolayer, 0.067 eV. Thus, the study elaborately predicts that the proposed heterostructure improves the intrinsic electrocatalytic activity of MoS₂.

Hydrogen production from the water-splitting process has attracted increasing attention to meet the global energy demand and provide a viable solution to environmental issues<fnref n="1" />. An economical process for hydrogen production is based on a high-performance surface Hydrogen Evolution Reaction (HER) on an appropriate electrocatalyst². Molybdenum and tungsten sulfides have been identified as promising noble-metal-free electro-catalysts, particularly for the HER process³⁻⁶. The basic concept of HER mechanisms has also been understood through relationships between computational approaches and corresponding experiments⁷⁻¹⁵. Hinnemann et al. proposed an approach based on the Density Functional Theory (DFT) in which they showed that the exposed edges of MoS₂ sheets are the active sites for the adsorption of hydrogen if the binding free energy of atomic hydrogen to the electrocatalyst is close to zero⁹. Nevertheless, a critical issue for the application of MoS₂ as an electrocatalyst in electrochemical reactions is attributed to its low electronic conductivity between two neighbor-ing S-Mo-S sheets, bonded by van der Waals (vdW) forces¹⁶. The resistivity through the basal planes was indeed determined to be 2200 times larger than that parallel to the planes¹⁶.

There have been basically two proposed ways to improve the MoS₂ electrocatalyst towards the HER: (1) increasing the density of active sites at the surface of the electrocatalyst; and (2) enhancing the electrical contact at these sites by reducing the number of layers and by placing MoS₂ on highly conductive substrates, such as carbon-based materials¹⁶. Consequently, the stacking of MoS₂ nanosheets with only a few layers perpendicular

<footnote n="1">¹Surface Reaction and Advanced Energy Materials Laboratory, Chemical Engineering Department, Amirkabir University of Technology (Tehran Polytechnic), Tehran, Iran. ²Department of Chemical Engineering, McGill University, 3610 University Street, Montreal, QC H3A 0C5, Canada. ³Department of Chemical Engineering, School of Engineering Technology and Industrial Trades, College of the North Atlantic - Qatar, Doha, Qatar. ⁴College of Engineering and Technology, American University of Middle East, Egaila, Kuwait. ⁵Center of Excellence in Electrochemistry, School of Chemistry, College of Science, University of Tehran, P.O. Box: 14155-6455, Tehran, Iran. ⁶CNRS, IJL, Université de Lorraine, 88000 Épinal, France. ✉email: Sajjad.habibzadeh@mail.mcgill.ca</footnote>
to a conductive substrate is expected to be an effective electrocatalyst. This can facilitate charge transfer along the edge of the electrode substrate to the active sites with minimal resistance¹⁶ while suppressing MoS₂ aggregation at the same time¹⁷,¹⁸. Recent progress in heterostructures based on vdW forces at the atomic level led to new categories of vertical quantum heterostructures with sharp atomic interfaces between materials having different physicochemical properties¹⁹. Such heterostructures, including two-dimensional (2D) crystalline layers, can provide interfaces with new physical and chemical characteristics that can be potentially employed in certain applications²⁰,²¹. Moreover, these 2D structures, thanks to the presence of strong covalent bonds, can already provide adequate in-plane stability. However, maintaining the stacking of these heterostructures together requires relatively weak vdW interactions²².

The design of atomic layers based on vdW heterostructures can be quite challenging when it comes to adjusting their electrocatalytic activity to the HER process, in order to render them more efficient than electrocatalysts based on metals and metal oxides²³,²⁴. It was found that no electrocatalytic activity was observed for the defect-free basal planes in most 2D layers²⁴. However, the intrinsic electrocatalytic properties of the individual layers for a specific reaction can be significantly changed depending on the design of the various heterostructures. Such changes can be due to the electric field created between the different layers of the active electrocatalysts and the respective conductive substrates¹⁹,²⁵,²⁶. Namely, the vdW stacking of hexagonal boron nitride and graphene makes the corresponding heterostructure a quite active electrocatalyst towards HER despite the inactive sites of each layer²,²³. In addition, to evaluate the electrochemical catalytic performance of vdW solids, the selection of the layers in the heterostructure and their sequence must be carefully considered. The importance of the stack layer sequence was addressed experimentally when the electrocatalytic activity of the vdW heterostructure of carbon nanotube (CNT) on MoS₂ showed a higher HER efficiency than the layered structure of MoS₂ on CNT²⁷,²⁸. In addition, a recent study has indicated that, when graphene is placed on top of MoS₂, a higher electrocatalytic efficiency can be obtained for the HER process in acidic solution compared to the inverse configuration²⁶.

Despite the excellent electronic properties of MoS₂/carbon-derived electrocatalysts for HER, only a few studies have been carried out on the electrochemical reaction pathways of these advanced materials²⁴,²⁶. Specifically, a DFT calculation on MoS₂/graphene complexes indicated that the presence of graphene as an underlayer of MoS₂ significantly affects the charge density distribution of MoS₂²⁴. In addition, the induced electric field of the MoS₂/graphene hybrid provides an excess negative charge density to the system, thus improving its HER activity²⁴. Furthermore, this sandwich configuration can also make the MoS₂ basal plane near the thermo-neutral Gibbs free energy change (ΔG<sub>H</sub> ~ 0), facilitating the activation of the MoS₂ basal plane towards the HER process. Although a few studies have focused on MoS₂/graphene, no report was found on the use of MoS₂/CNT as a HER electrocatalyst. In the present study, we used DFT calculations to understand the molecular mechanism of the stacking sequence and subsequent layers of MoS₂/CNT in the HER process. Namely, the electrocatalytic activity of a new MoS₂/CNT heterostructure and the effect of MoS₂ on the structural and electronic characteristics of the CNT substrate for HER were theoretically explored.

## Methodology

The relaxation of the geometry and the calculation of the electronic structure were performed by the DFT approach. The stacked heterostructure of MoS₂/CNT was investigated as our main system. In addition, an 8 × 8 CNT (a: 14.01 b: 13.91 c: 17.3 Å) containing 224 carbon atoms was applied to match a 4 × 4 MoS₂ monolayer (a: 12.66 b: 12.66 c: 18.4 Å) comprising 16 molybdenum and 32 sulfur atoms. The lattice mismatch of the MoS₂ and CNT layers was approximately 5%. In addition, a void space of 15 Å was considered on the Z-axis to ignore possible interactions between the periodic structures. The relaxation process was carried out for both the MoS₂ and CNT layers, atoms, and cells, and then the binding energy of the heterostructure was computed using the following equation:

$$E_b = E_{total} - E_{CNT} - E_{MoS_2} \quad (1)$$

where E<sub>total</sub>, E<sub>CNT</sub>, and E<sub>MoS₂</sub> refer to the total energy of the MoS₂/CNT heterostructure, the energy of a single CNT, and the energy of a single-layer of MoS₂, respectively. It is worth mentioning that the stable heterostructure was such that the total energy of the MoS₂/CNT heterostructure is lower than the energy of individual CNT and MoS₂ (see Eq. (1)). Moreover, the DFT method was used based on the Dmol³ code with the Generalized Gradient Approximation (GGA) approach in the Materials Studio version 7.0 package²⁹,³⁰. Furthermore, the exchange-correlation functional used in this study was based on the work of Perdew, Burke and Ernzerhof (PBE)³¹ with the Gaussian double zeta plus polarization numerical base (DNP) set. To treat the core electrons, DFT semi-core pseudopotentials (DSPPs) were selected. The geometry relaxation and energy computations were chosen with 1 × 10⁻⁵ Ha, 0.002 Ha/Å, and 0.005 Å for energy, force, and displacement tolerances, respectively. Thanks to Grimme's semi-empirical dispersion-corrected density functional theory (DFT-D2)³²,³³ for considering weak interactions with high accuracy, we used DFT-D2 instead of the standard PBE functional. Besides, to analyze the characteristics of the electron density difference, the CASTEP code³⁴ of plane wave and ultra-soft pseudopotentials³⁵ was applied with a plane-wave cutoff energy of 400 eV.

## Results and discussion

**Optimized structures.** The optimized configurations of the CNT, MoS₂, and MoS₂ deposited on the CNT substrate are shown in Fig. 1. There is a good agreement between the Mo-S bond length of our optimized structure (2.42 Å) and the experimental one (2.41 Å)³⁶,³⁷. In addition, a distance of 3.37 Å was optimally achieved between the CNT and the nearest sulfur layer. Since both of the above distances are greater than 1.81 Å (the sum of the covalent radii of carbon and sulfur atoms), vdW forces might be established between the MoS₂ and the CNT substrate as the main interactions. Furthermore, the average C-C optimized bond length in the
<figure><img src="image_3.png" /><figcaption>**Figure 1.** The most stable configurations of: (a) CNT; (b) MoS₂; (c) top view of MoS₂/CNT; and (d) side view of MoS₂/CNT. Gray, green and yellow spheres stand for C, Mo and S atoms, respectively.</figcaption></figure>

CNT nanostructure studied in our work was calculated to be about 1.42 Å, which is in line with the previously reported DFT studies (1.43 Å)³⁸,³⁹. Moreover, these results agree with previously published works on graphene/ZnO⁴⁰ and phosphorene/graphene⁴¹,⁴² equilibrium distances.

The stability of the MoS₂/CNT interface was evaluated by computing the binding energy per C atom between MoS₂ and the CNT substrate. The binding energy obtained was approximately 0.467 eV/atom, similar to the values reported for MoS₂/graphene⁴³ and MoS₂/N-graphene⁴⁴. Therefore, the value of binding energy can also be considered as a further evidence of the weak vdW interaction forces existing in the MoS₂/CNT system.

**Density of states (DOS).** The electronic properties of the system can be determined through the DOS while considering the HOMO–LUMO bandgap energy ($E_g$). This term is the minimum energy required to excite electrons from the valence band to the conduction band. A lower $E_g$ can result in greater system conductivity and hydrogen adsorption capacity. During the formation of the vdW heterostructure, the electronic structure of the CNT on the one hand and of the MoS₂ monolayer on the other hand were altered near the Fermi energy level. As seen in Fig. 2, the DOS of the MoS₂/CNT structures was compared to the DOS of the isolated constituent monolayers.

Furthermore, the projected density of states (PDOS) for the C-p orbitals in the CNT and the S-p and Mo-d orbitals in MoS₂ are shown in Fig. 2a,b. It should be mentioned that CNTs are two dimensional Dirac materials with a linear dispersion near the Fermi energy level⁴⁵,⁴⁶, as observed in Fig. 2a. In addition, the MoS₂ monolayer as a semiconductor possesses a bandgap of about 1.8 eV⁴⁷. However, the bandgap obtained for MoS₂ was calculated to be 2.08 eV in the current simulation (see Fig. 2b). This difference of $E_g$ can be attributed to the unfilled d-orbital in the Mo atoms. Thus, semi-local DFT functionals cannot be used to calculate the corresponding $E_g$ correctly, which requires the application of many-body corrections by Green's function (GW)⁴⁸ or hybrid Heyd, Scuseria, and Ernzerhof (HSE) functionals⁴⁹ to compensate for the bandgap differences of 0.28 eV in the semi-local DFT functionals. However, to reduce computational costs, such bandgap corrections were ignored because the semi-empirical DFT-D2 method can achieve sufficient accuracy for the calculation of structural and electronic properties of the MoS₂/CNT interface⁵⁰.

Furthermore, the bottom of the conduction band and the top of the valence band originate mainly from the Mo-d orbitals and both the Mo-d and S-p orbitals, respectively (see Fig. 2b). In addition, it can be seen that the Mo-d and S-p orbitals were hybridized together at the top of the valence band (see Fig. 2b). It can be concluded from Fig. 2a that the half-filled p orbitals perpendicular to the planar structure create the π and π* bands in the electronic configuration of CNT. Moreover, at the corner of the Brillouin zone of the CNT, both bonding and antibonding bands touch at a single point near the Fermi energy level. Furthermore, the PDOS of MoS₂/CNT (see Fig. 2c) is displayed relative to the partial d-DOS (blue) of Mo in MoS₂ and the p-DOS (red) of C in the CNT. It can be seen that the Fermi energy level of the MoS₂/CNT heterostructure is characterized by the Dirac-cone-like characteristic from CNT and a gap-like characteristic from MoS₂ (see Fig. 3c). Furthermore, the location of PDOSs for CNT is similar to that of MoS₂/CNT where there is also a remarkable change in the intensity and profile of the PDOSs of MoS₂.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** PDOS diagrams for: (a) CNT; (b) MoS₂; and (c) MoS₂/CNT. The Fermi energy level is set at zero.</figcaption></figure>

**Band structure analysis.** To evaluate the impact of stacking configurations on the bandgap and electronic properties of the MoS₂/CNT heterostructure, structural and electronic calculations were performed separately on the freestanding CNT and the MoS₂ monolayer. As seen in Fig. 3a, the bandgap obtained for an isolated CNT is determined to be 0.017 eV. Besides, a bandgap structure based on a linear Dirac dispersion can be observed for the isolated CNT near the Fermi energy level. The band structure of MoS₂ is illustrated in Fig. 3b, which shows a direct bandgap of 2.01 eV with the conduction and valence bands positioned at the K point. The resultant bandgap is greater than the reported experimental value (about 1.80 eV)⁴⁷. However, this discrepancy can be corrected by applying the GW approximation technique, which is not the subject of our study.

The band structure of the MoS₂/CNT heterostructure can be determined by the energy bands of the CNT and MoS₂. Figure 3c shows the linear dispersion bands of the CNT, which are located in the large energy gap of MoS₂ while the electronic energy band of the pristine CNT can be found without any major change. However, a significant change near the Fermi energy level can be observed. The results obtained from the band structure can evidently show the efficient interactions between CNT and MoS₂ that improve the electronic properties of the CNT, the band gap of the heterostructure having dropped to 0.008 eV.

**Determination of charge density and charge transfer.** Figure S1 illustrates the variations in the average atomic charge (Δq) on the sulfur and molybdenum atoms specified by Hirshfeld's charge analysis. S1–S4 and M1–M4 represent different positions of sulfur and molybdenum atoms, respectively, in the MoS₂/CNT heterostructure. The average atomic charge obtained on the MoS₂ monolayer for S and Mo atoms was calculated to be –0.114 and 0.229 a.u., respectively. Table S1 summarizes the charges on Mo and S atoms in the MoS₂/CNT heterostructure. It should be noted that the charge on Mo atoms remains almost unchanged at the value of ~0.229 a.u, whereas it changes negatively and significantly on sulfur atoms.

Figure 4 illustrates the difference in charge density of the current MoS₂/CNT system. The charge depletion is observed at the two middle neighboring S and C planes. Moreover, no orbital overlap between the MoS₂ layer and the CNT can be observed due to the weak vdW interactions between MoS₂ and the CNT. Besides, the charge transfer, involving the total sum of the Hirshfeld charge populations, was analyzed and computed for the MoS₂ layer. A negative charge value means that the charge is transferred from the CNT to MoS₂, while the charge
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Band structures of: (a) CNT; (b) MoS₂; and (c) MoS₂/CNT.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Spatial mapping of charge density differences for: (a) a MoS₂ monolayer; and (b) the MoS₂/CNT heterostructure. Regions of electron accumulation and depletion are denoted by blue and red lobes, respectively. BIOVIA, Dassault Systèmes, Materials Studio, version 7. https://bit.ly/38lRRQR.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>Figure 5. Minimum-energy pathway of the Volmer reaction on: (a) MoS₂; and (b) the MoS₂/CNT heterostructure. IS, TS and FS stand for initial, transition and final states, respectively. BIOVIA, Dassault Systèmes, Materials Studio, version 7. https://bit.ly/381RRQR.</figcaption></figure>

transfer from MoS₂ to the CNT reflects its positive charge value. Hence, a negative charge value of -0.051 a.u. obtained for MoS₂/CNT implies a charge transfer from the CNT to the MoS₂ monolayer.

**Mechanism of the Hydrogen Evolution Reaction (HER).** The superior electronic structure of MoS₂/CNT relative to pristine MoS₂ means that this heterostructure can be used effectively as an electrocatalyst for the HER process. A well-known mechanism of HER, known as Volmer-Heyrovsky or Volmer-Tafel mechanism, can be expressed as follows⁵¹,⁵²:

$$2\text{H}^+ + 2\text{e}^- \rightarrow 2\text{H}^* \quad (\text{Volmer}) \qquad (2)$$

$$\text{H}^* + \text{H}^+ \rightarrow \text{H}_2(\text{g}) \quad (\text{Heyrovsky}) \qquad (3)$$

$$\text{H}^* + \text{H}^* \rightarrow \text{H}_2(\text{g}) \quad (\text{Tafel}) \qquad (4)$$

According to this mechanism, a hydronium ion first adsorbs on the surface and forms a hydrogen radical. Then, two hydrogen radicals combine and form a hydrogen molecule. Basically, HER activity can be investigated by the adsorption of hydronium ions and the activation energy barrier. The adsorption energy of the hydronium ion adsorbed on MoS₂ and MoS₂/CNT monolayer are calculated to be 0.0057 eV and 0.0039 eV, respectively, indicating an easier adsorption of hydronium ion on the MoS₂/CNT heterostructure than on the MoS₂ monolayer.

To characterize the HER activity of the proposed heterostructure, we considered the initial step of the HER process (Volmer reaction) as the rate-determining step (RDS). As illustrated in Fig. 5, the path of minimum energy for the transfer of one of the solvated protons to the MoS₂ surface in the 4×4 supercell consists of three steps: initial, transition and final states (IS, TS and FS, respectively). It can be seen that the adsorption of H atom at each step occurs on the edge of S atom, signifying the electrocatalytic activity of the edge S atoms towards HER. In addition, the energy barrier of MoS₂/CNT and MoS₂ monolayer is calculated to be 0.024 eV and 0.067 eV, respectively. This observation might be ascribed to the electron redistribution of the edge S atom after the adsorption of the CNT⁵³. From the Hirshfeld charge analysis, the edge S atom of MoS₂/CNT gains -0.190 e (see Table S1) while that of pristine MoS₂ acquires -0.114 e. This suggests that due to electrostatic attraction, the edge S atom in MoS₂/CNT has a more negative charge and thus a greater interaction with the H atom. As a result, it can be seen that the CNT reduces the energy barrier (or the onset potential of HER), thus improving the intrinsic activity of MoS₂.

Moreover, the other two primary steps which can be plausible for H₂ evolution in the second step of the HER process encompass Heyrovsky and Tafel reactions. In the case of Heyrovsky reaction (as the initial step), H₂ molecule is formed through the reaction of the proton (in the water layer) with an adsorbed hydrogen (see Eq. (3)). Figure 6 displays the estimated minimum-energy paths at two dissimilar structures. It is evident that an adsorbed H atom on a sulfur one reaches an H atom of a hydronium ion in the water layer. This is followed by breaking the adsorbed H from the surface, forming H₂ molecule within the water layer. In the TS (Fig. 6), the interfacial adsorbed H is separated from the surface and the S-H distance increased from 1.301 Å in the IS to 2.391 Å in the MoS₂ monolayer. At the same time, the proton from H₃O⁺ travels toward the separated H atom, forming a molecule with H-H bond lengths of 0.760 Å and 0.756 Å in the MoS₂ monolayer and MoS₂/CNT heterostructure, respectively. The evolved H₂ molecule is detached from the surface in the FS. An activation energy of 0.68 eV was determined for MoS₂ while such energy barrier drops to 0.41 eV for the MoS₂/CNT heterostructure (Fig. 6). Thus, the energy barrier is far greater for the Heyrovsky reaction than that for the Volmer reaction, suggesting that the H desorption procedure would be the RDS of the Volmer-Heyrovsky pathway.
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Minimum-energy pathway of the Heyrovsky reaction on: (a) MoS₂; and (b) MoS₂/CNT heterostructure. BIOVIA, Dassault Systèmes, Materials Studio, version 7. https://bit.ly/38lRRQR.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 7.** Minimum-energy pathway of the Tafel reaction on: (a) MoS₂; and (b) the MoS₂/CNT heterostructure. BIOVIA, Dassault Systèmes, Materials Studio, version 7. https://bit.ly/38lRRQR.</figcaption></figure>

Furthermore, adjoining H atoms coupling to surface S atoms perform the Tafel reaction on MoS₂ (see Eq. (4)). Figure 7 depicts the minimum-energy path of direct recombination of two adjoining protons adsorbed on S active sites. Due to the absence of charge transfer over the interface, the entire energies along the path are not amended for possible deviance. In the IS, the distances between the two H atoms are 3.922 Å and 3.817 Å in the MoS₂ monolayer and MoS₂/CNT heterostructure, respectively. The two S–H bond lengths in the MoS₂ monolayer and MoS₂/CNT heterostructure are 1.309 Å and 1.324 Å, respectively. However, the formed H₂ molecule undergoes desorption from the surface with a H–H bond length of 0.750 Å for both of the developed structures in the FS. In the MoS₂ monolayer and MoS₂/CNT heterostructure, the estimated energy barrier of direct recombination values are 1.27 eV and 0.56 eV, respectively, which are markedly greater than that for the Volmer reaction. Therefore, the RDS in the Volmer–Tafel reaction is the H desorption (Tafel step).

Overall, a comparison of the energy barriers for the Heyrovsky and Tafel reactions on MoS₂ monolayer and MoS₂/CNT heterostructure reveals that the Tafel reaction should surmount elevated barriers of 1.27 and 0.56 eV, respectively, while it is possible for the Heyrovsky reaction to proceed more straightforwardly. This indicates more efficiency of the Heyrovsky reaction where the Volmer–Heyrovsky mechanism is the major pathway of HER. In addition, finding from the recent reports shows that HER potentially occurs through the Volmer–Heyrovsky process on the basal plane of 1 T-MoS₂, 2H-MoS₂ and Ni-MoS₂/RGO⁵²,⁵⁴,⁵⁵.

## Conclusion

The present study theoretically investigates the synergistic electrocatalytic activity of the MoS₂/CNT heterostructure towards the HER process by applying DFT simulations. The results indicated a weak van der Waals interaction between the CNT and the MoS₂ monolayer. Moreover, a distance of 3.37 Å was determined between them, and the binding energy per C atom in this system was found to be approximately 0.467 eV. The bandgap structure indicated that the linear Dirac-like dispersion of CNT near the Fermi energy level remains unchanged in the MoS₂/CNT interface as well. However, it was found that a bandgap around 8 meV was calculated at the Dirac K-point of the CNT in the MoS₂/CNT interface. Finally, it was confirmed that the presence of CNT can
improve the electronic conductivity while reducing the energy barrier in the MoS₂/CNT heterostructure for the HER process.

Received: 24 October 2020; Accepted: 18 January 2021

Published online: 17 February 2021

References

1. Dresselhaus, M. S. & Thomas, I. L. Alternative energy technologies. *Nature* **414**, 332 (2001).

2. Seh, Z. W. et al. Combining theory and experiment in electrocatalysis: Insights into materials design. *Science* (80-. ) 355, eAAD4998 (2017).

3. Xue, X. et al. Rational inert-basal-plane activating design of ultrathin 1T' phase MoS2 with a MoO₃ heterostructure for enhancing hydrogen evolution performances. *Nanoscale* **10**, 16531–16538 (2018).

4. Cao, J., Zhou, J., Zhang, Y., Wang, Y. & Liu, X. Dominating role of aligned MoS2/Ni3S2 nanoarrays supported on three-dimensional Ni foam with hydrophilic interface for highly enhanced hydrogen evolution reaction. *ACS Appl. Mater. Interfaces* **10**, 1752–1760 (2018).

5. Gnanasekar, P., Periyanagounder, D. & Kulandaivel, J. Vertically aligned MoS2 nanosheets on graphene for highly stable electro-catalytic hydrogen evolution reactions. *Nanoscale* **11**, 2439–2446 (2019).

6. Bae, C. et al. Bulk layered heterojunction as an efficient electrocatalyst for hydrogen evolution. *Sci. Adv.* **3**, e1602215 (2017).

7. Kibsgaard, J., Chen, Z., Reinecke, B. N. & Jaramillo, T. F. Engineering the surface structure of MoS2 to preferentially expose active edge sites for electrocatalysis. *Nat. Mater.* **11**, 963 (2012).

8. Xie, J. et al. Defect-rich MoS2 ultrathin nanosheets with additional active edge sites for enhanced electrocatalytic hydrogen evolution. *Adv. Mater.* **25**, 5807–5813 (2013).

9. Hinnemann, B. et al. Biomimetic hydrogen evolution: MoS2 nanoparticles as catalyst for hydrogen evolution. *J. Am. Chem. Soc.* **127**, 5308–5309 (2005).

10. Laursen, A. B., Kegnæs, S., Dahl, S. & Chorkendorff, I. Molybdenum sulfides—efficient and viable materials for electro-and photo-electrocatalytic hydrogen evolution. *Energy Environ. Sci.* **5**, 5577–5591 (2012).

11. Merki, D. & Hu, X. Recent developments of molybdenum and tungsten sulfides as hydrogen evolution catalysts. *Energy Environ. Sci.* **4**, 3878–3888 (2011).

12. Benck, J. D., Chen, Z., Kuritzky, L. Y., Forman, A. J. & Jaramillo, T. F. Amorphous molybdenum sulfide catalysts for electrochemical hydrogen production: insights into the origin of their catalytic activity. *ACS Catal.* **2**, 1916–1923 (2012).

13. Merki, D., Fierro, S., Vrubel, H. & Hu, X. Amorphous molybdenum sulfide films as catalysts for electrochemical hydrogen production in water. *Chem. Sci.* **2**, 1262–1267 (2011).

14. Li, Y. et al. MoS2 nanoparticles grown on graphene: an advanced catalyst for the hydrogen evolution reaction. *J. Am. Chem. Soc.* **133**, 7296–7299 (2011).

15. Jaramillo, T. F. et al. Identification of active edge sites for electrochemical H2 evolution from MoS2 nanocatalysts. *Science* (80-. ) **317**, 100–102 (2007).

16. Tributsch, H. & Bennett, J. C. Electrochemistry and photochemistry of MoS2 layer crystals. I. *J. Electroanal. Chem. Interfacial Electrochem.* **81**, 97–111 (1977).

17. Yan, Y. et al. Facile synthesis of low crystalline MoS2 nanosheet-coated CNTs for enhanced hydrogen evolution reaction. *Nanoscale* **5**, 7768–7771 (2013).

18. Yuan, H., Li, J., Yuan, C. & He, Z. Facile synthesis of MoS2@ CNT as an effective catalyst for hydrogen production in microbial electrolysis cells. *ChemElectroChem* **1**, 1828–1833 (2014).

19. Shi, J. et al. Temperature-mediated selective growth of MoS2/WS2 and WS2/MoS2 vertical stacks on Au foils for direct photocatalytic applications. *Adv. Mater.* **28**, 10664–10672 (2016).

20. Geim, A. K. & Van der Grigorieva, I. V. Waals heterostructures. *Nature* **499**, 419–425 (2013).

21. Shi, Y., Zhang, H., Chang, W.-H., Shin, H. S. & Li, L.-J. Synthesis and structure of two-dimensional transition-metal dichalcogenides. *MRS Bull.* **40**, 566–576 (2015).

22. Ajayan, P., Kim, P. & van der Banerjee, K. Waals materials. *Phys. Today* **69**, 9–38 (2016).

23. Bawari, S. et al. On the hydrogen evolution reaction activity of graphene–hBN van der Waals heterostructures. *Phys. Chem. Chem. Phys.* **20**, 15007–15014 (2018).

24. Li, H. et al. Charge-transfer induced high efficient hydrogen evolution of MoS2/graphene cocatalyst. *Sci. Rep.* **5**, 18730 (2015).

25. Araujo, P. T., Terrones, M. & Dresselhaus, M. S. Defects and impurities in graphene-like materials. *Mater. Today* **15**, 98–109 (2012).

26. Biroju, R. K., Pal, S., Sharma, R., Giri, P. K. & Narayanan, T. N. Stacking sequence dependent photo-electrocatalytic performance of CVD grown MoS2/graphene van der Waals solids. *Nanotechnology* **28**, 85101 (2017).

27. Zhang, J. et al. SWCNT-MoS2-SWCNT vertical point heterostructures. *Adv. Mater.* **29**, 1604469 (2017).

28. Biroju, R. K. et al. Hydrogen evolution reaction activity of graphene–MoS2 van der Waals heterostructures. *ACS Energy Lett.* **2**, 1355–1361 (2017).

29. Delley, B. An all-electron numerical method for solving the local density functional for polyatomic molecules. *J. Chem. Phys.* **92**, 508–517 (1990).

30. Delley, B. From molecules to solids with the DMol3 approach. *J. Chem. Phys.* **113**, 7756–7764 (2000).

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

32. Grimme, S. Accurate description of van der Waals complexes by density functional theory including empirical corrections. *J. Comput. Chem.* **25**, 1463–1473 (2004).

33. Grimme, S. Semiempirical GGA-type density functional constructed with a long-range dispersion correction. *J. Comput. Chem.* **27**, 1787–1799 (2006).

34. Segall, M. D. et al. First-principles simulation: ideas, illustrations and the CASTEP code. *J. Phys. Condens. Matter* **14**, 2717 (2002).

35. Vanderbilt, D. Soft self-consistent pseudopotentials in a generalized eigenvalue formalism. *Phys. Rev. B* **41**, 7892 (1990).

36. Wilson, J. A. & Yoffe, A. D. The transition metal dichalcogenides discussion and interpretation of the observed optical, electrical and structural properties. *Adv. Phys.* **18**, 193–335 (1969).

37. Lévy, F. A. *Intercalated Layered Materials* Vol. 6 (Springer, Berlin, 2012).

38. Mashhadzadeh, A. H., Fereidoon, A. & Ahangari, M. G. Surface modification of carbon nanotubes using 3-aminopropyltriethoxysilane to improve mechanical properties of nanocomposite based polymer matrix: experimental and density functional theory study. *Appl. Surf. Sci.* **420**, 167–179 (2017).

39. Ghalkhani, M., Beheshtian, J. & Salehi, M. Electrochemical and DFT study of an anticancer and active anthelmintic drug at carbon nanostructured modified electrode. *Mater. Sci. Eng. C* **69**, 1345–1353 (2016).

40. Xu, P., Tang, Q. & Zhou, Z. Structural and electronic properties of graphene–ZnO interfaces: dispersion-corrected density functional theory investigations. *Nanotechnology* **24**, 305401 (2013).
41. Cai, Y., Zhang, G. & Zhang, Y.-W. Electronic properties of phosphorene/graphene and phosphorene/hexagonal boron nitride heterostructures. *J. Phys. Chem. C* **119**, 13929–13936 (2015).

42. Guo, G.-C. et al. First-principles study of phosphorene and graphene heterostructure as anode materials for rechargeable Li batteries. *J. Phys. Chem. Lett.* **6**, 5002–5008 (2015).

43. Ahmed, T., Modine, N. A. & Zhu, J.-X. Graphene/MoS2 van der Waals Bilayer as the Anode Material for Next Generation Li-ion Battery: A First-Principles Investigation. *arXiv Prepr. arXiv1502.07398* (2015).

44. Zan, W., Geng, W., Liu, H. & Yao, X. Influence of interface structures on the properties of molybdenum disulfide/graphene composites: a density functional theory study. *J. Alloys Compd.* **649**, 961–967 (2015).

45. Wehling, T. O., Black-Schaffer, A. M. & Balatsky, A. V. Dirac materials. *Adv. Phys.* **63**, 1–76 (2014).

46. Abergel, D. S. L., Apalkov, V., Berashevich, J., Ziegler, K. & Chakraborty, T. Properties of graphene: a theoretical perspective. *Adv. Phys.* **59**, 261–482 (2010).

47. Mak, K. F., Lee, C., Hone, J., Shan, J. & Heinz, T. F. Atomically thin MoS2: a new direct-gap semiconductor. *Phys. Rev. Lett.* **105**, 136805 (2010).

48. Shi, H., Pan, H., Zhang, Y.-W. & Yakobson, B. I. Quasiparticle band structures and optical properties of strained monolayer MoS2 and WS2. *Phys. Rev. B* **87**, 155304 (2013).

49. Ellis, J. K., Lucero, M. J. & Scuseria, G. E. The indirect to direct band gap transition in multilayered MoS2 as predicted by screened hybrid density functional theory. *Appl. Phys. Lett.* **99**, 261908 (2011).

50. Hieu, N. N. et al. First-principles study of the structural and electronic properties of graphene/MoS2 interfaces. *J. Appl. Phys.* **122**, 104301 (2017).

51. Lei, X., Yu, K., Li, H., Tang, Z. & Zhu, Z. First-principle and experiment framework for charge distribution at the Interface of the molybdenum dichalcogenide hybrid for enhanced electrochemical hydrogen generation. *J. Phys. Chem. C* **120**, 15096–15104 (2016).

52. Tang, Q. & Jiang, D. Mechanism of hydrogen evolution reaction on 1T-MoS2 from first principles. *Acs Catal.* **6**, 4953–4961 (2016).

53. Ma, X.-C., Dai, Y., Yu, L. & Huang, B.-B. Energy transfer in plasmonic photocatalytic composites. *Light Sci. Appl.* **5**, e16017–e16017 (2016).

54. Chen, L. X., Chen, Z. W., Wang, Y., Yang, C. C. & Jiang, Q. Design of dual-modified MoS2 with nanoporous Ni and graphene as efficient catalysts for the hydrogen evolution reaction. *ACS Catal.* **8**, 8107–8114 (2018).

55. Li, W. et al. Hydrogen evolution reaction mechanism on 2H-MoS2 electrocatalyst. *Appl. Surf. Sci.* **498**, 143869 (2019).

## Author contributions

F.K., and A.B. contributed to the conceptualization, model development, results and discussion. S.H. conceived the theory and supervised the project. A.E., A.M., A.H.M., M.R.G., M.R.S. worked on the first draft and co-wrote the manuscript. V.F., A.C. contributed to the data analysis.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-021-83562-w.

**Correspondence** and requests for materials should be addressed to S.H.

**Reprints and permissions information** is available at www.nature.com/reprints.

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

<figure><img src="image_10.png" /></figure>

**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 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://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021