1 Introduction
Low-dimensional materials with extraordinary peculiarities, such as zero-dimensional buckyballs and one-dimensional carbon nanotubes, have empowered their potentiality in many important fields [
1,
2]. Their deficiency in distinguishing metal and semiconductor phases has led to the boomingly development of two-dimensional (2D) nanomaterials, including graphene, boron nitride, and transitional metal dichalcogenides (TMDCs), etc. [
3-
9]. In comparison with zero-bandgap graphene, layered TMDCs (MoS
, WS
, MoSe
, WSe
, etc.) exhibit a tunable bandgap of about 1−2 eV, showing an interesting capability for miniaturizing electronics. Among them, MoS
is the most optimistic TMDCs because of its surplus and innoxious and easily feasible synthesis [
10,
11]. Structurally, MoS
is a typical hexagonal layered compound with strong in-plane covalent and weak out-of-plane van der Waals (vdW) interaction, resulting in the absence of dangling bonds on the surface [
12]. When cleaved from bulk to monolayer (1L), energy band structure of MoS
transforms from indirect to direct due to the quantum confinement effect [
6]. Field-effect transistors (FETs) based on monolayer MoS
have been demonstrated to exhibit a high on/off current ratio 10
8 and a high mobility of 200 cm
2·V
−1·s
−1 at room temperature [
13]. These intriguing properties give it a wide range of potential applications in transistors [
13-
15], memories [
16,
17], flexible electronics [
18,
19], photodetectors [
20-
23], energy storage and conversion and hydrogen evolution reaction [
24].
In the past few decades, wide bandgap semiconductors (WBSs), such as GaN and SiC, have drawn considerable attention because of their outstanding properties, and widely used in the emerging electronics and optoelectronics. GaN possesses wide bandgap (3.4 eV), high room-temperature electron mobility (1300 cm
2·V
−1·s
−1), and large short-wavelength absorption coefficient, making it an important semiconductor for applications in light emitting diodes (LED), FETs, solar cells, lasers, and photodetectors [
25-
29]. SiC exists as over 170 polytypes with the 4H polytype is most commonly used in high-temperature and high-frequency electronic devices due to its high carrier mobility (900 cm
2·V
−1·s
−1), large band gap (3.26 eV) and low dopant ionization energy [
30].
Beyond separate studies of 3D WBSs or 2D TMDCs, it is favourable to integrate 2D/3D TMDCs/WBSs heterostructures, which benefits for improving the performance of devices, such as LEDs, photodetectors, sensors and energy storage devices, and further providing the material platform of novel physical phenomena [
25]. Due to the great light absorption of multilayer MoS
, Goel et al. have successfully fabricated the MoS
/GaN UV photodetector with the detectivity of 10
11 Jones and the external spectral responsivity of 10
3 A·W
−1 [
31]. The MoS
/GaN heterostructure device fabricated by Moun et al. demonstrated a detectivity of 10
14 Jones and an ultra-high optical responsivity of 10
5 A·W
−1 under a 405 nm laser radiation [
32]. Because of its suitable band gap, SiC has always been an important material for UV photodetectors [
33]. Then a MoS
/4H-SiC photodetector operating in both the UV and visible regions were fabricated, which exhibits a responsivity of 5.7 A·W
−1 under UV irradiation and a noise equivalent power (NEP) of 10
−13−10
−15 W·Hz
−1/2 [
20]. Based on the dual-photogating effect, MoS
/SiC photodetectors were demonstrated to have a ultra-high responsivity of 1.6 × 10
3 A·W
−1 and 1 × 10
4 A·W
−1 in VIS and UV irradiation, respectively [
21]. The MoSe
/WSe
heterojunction photodetector fabricated on n-doped 4H-SiC also shows the excellent performance with the maximum responsivity of 7.17 A·W
−1, the corresponding maximum EQE of 1.67 × 10
3 %, and detectivity of 5.51 × 10
11 Jones with a gain of 10
3 [
34].
In the aspects of photocatalyst, gas sensor and other applications, the MoSSe/SiC and WSSe/SiC heterostructures were certified to be a high-efficient photocatalyst due to their specific electronic and optical properties based on the density functional theory (DFT) calculations [
35]. And the MoS
/GaN photoanode could enhance the photocurrent density to 5.2 mA·cm
−2, which is 2.6 times higher than that of GaN photoanode [
36]. As compared to the MoS
or GaN gas sensors, the sensitivity can be significantly enhanced (20 times) when integrating the MoS
/GaN heterostructure [
37,
38]. The MoS
/GaN heterojunction nanowires could be used in memory device, which possesses the good retention characteristics of 3.4 × 10
3 s under a low switching voltage [
39]. Although above experimental and theoretical works have confirmed the enormous advantage of the integrated TMDCs/WBSs, the detailed study and analysis of the interfacial properties, especially multiple interfacial interactions modulated energy band alignment, of these heterostructures are still quite limited.
In this study, monolayer (1L), bilayer (2L) and quadlayer (4L) MoS were deposited on 4H-SiC (0001) substrates by reactive sputtering deposition. The interfacial properties of MoS/SiC heterostructures were studied by combining the first-principles calculations and X-ray photoelectron spectroscopy. Experimental (theoretical) VBO increases from 1.49 (1.46) to 2.19 (2.36) eV, with increasing MoS layers up to 4L. Theoretical calculations further reveal that a strong interlayer interaction was demonstrated at 1L MoS/SiC interface and Fermi level pinning and totally surface passivation were realized for the 4H-SiC (0001) Si surface. 1L MoS/SiC interface exhibits type I band alignment with the asymmetric CBO and VBO. For 2L and 4L MoS/SiC, Fermi level was just pinning at the lower MoS monolayer. They exhibit the type II band alignments and the enlarged CBOs and VBOs. High efficiency of charge separation will emerge due to the asymmetric band alignment and built-in electric field for all the MoS/SiC interfaces, and enhance the device performance significantly.
2 Experimental and theoretical methods
The 4H-SiC substrates were prepared using a four-step cleaning process. It was ultra-sonicated in acetone, ethanol and DI water for 10 mins each at 45 °C to remove particle and organic contaminants. After which, dilute HF (1:1) was used to etch the surface to remove any SiO oxidation. MoS layers were grown via reactive magnetron sputtering whereby a molybdenum target was sputtered using argon in a chamber with ambient sulfur. The partial pressure of sulfur vapour was controlled at 1.2 × 10−5 mbar through the heating of sulfur powder in a separate compartment. One of the substrates was left bare while the rest were sputtered with a DC power source at 10 W, 3.2 × 10−4 mbar and 700 °C for 5 mins, 9 mins, 12 mins and 2 hours (1L, 2L, 4L and bulk MoS).
Thermo Scientific DXR microscope with a 514.5 nm laser were used to take the Raman and photoluminescence (PL) spectra of the as-grown heterojunction samples. Atomic force microscopy (AFM) images were obtained by using the Bruker Dimension Icon system and the tapping mode. High-resolution transmission electron microscopy (TEM) images were taken from samples grown directly on 8 nm SiO support membranes. High resolution X-ray photoelectron spectroscopy (XPS) spectra were taken on a VG ESCALAB 220i-XL XPS system. The binding energy scale was calibrated with pure Au, Ag, and Cu standard samples by setting the Au 4f, Ag 3d, and Cu 2p at binding energies of 83.98 eV, 368.26 eV, and 932.67 eV, respectively. C 1s signal from the sample surface was used to correct the core-level binding energy.
The first-principles calculations were performed by using the Vienna ab initio simulation package (VASP) [
40]. The exchange-correlation energy is described by the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional [
41]. The semi-empirical correction scheme of Grimme (DFT-D3) were utilized to correct the effect of vdW interactions in multilayer TMDCs system [
42]. Furthermore, Heyd−Scuseria−Ernzerhof (HSE06) was used to simulate all the density of states (DOSs), because of its accuracy for obtaining the bandgap of the calculated system [
43]. We also considered the dipole correction in the calculation to eliminate possible artificial dipole interactions [
44]. The charge transfer was calculated by using the Bader analysis [
45].
3 Results and discussion
Based on the above sputtering parameters, we could obtain the extremely low growth rate of MoS
layers. The successful deposition of monolayer and few-layer MoS
on different substrates have been achieved by precisely controlling the growing time, which has also been demonstrated by our previous works [
46-
50]. Raman spectroscopy was further used to characterize the number of layers in our samples. It has been shown in previous studies that the difference in Raman shift between the in-plane (
) and out-of-plane (
) characteristic vibration modes increases with the number of layers of MoS
and hence can be used to accurately determine the number of layers within 1 to 5 MoS
layers [
51–
53]. Fig.1(a) shows the
and
vibration modes of 1L, 2L, 4L and bulk MoS
/SiC samples. A Raman shift difference of 19.5, 21.7, 23.9 and 25.5 cm
−1 was observed from the samples, which corresponds to 1L, 2L, 4L and bulk MoS
respectively.
AFM images of SiC substrate before etching, after etching, and 1L MoS
grown on SiC substrate were taken to determine the uniformity of the samples, as shown in Fig.1(d)−(f). The step-terrace structure on the surface is due to different etching rates of the atomic layers in the 4H-SiC crystal structure during the material processing [
54]. The dilute HF wet etching process only resulted in a slight increase in RMS value from 0.09 nm to 0.11 nm and the deposition of MoS
did not result in an increase in surface roughness showing good monolayer uniformity.
Fig.1(b) exhibits annular dark field STEM (ADF-STEM) image of monolayer MoS with alternating Mo and S columns arranging into hexagonal structure, which can be identified to be the 2H-MoS polytype. Meanwhile, moir patterns (overlap on the lower left) are also observed. Therefore, TEM results confirm the high quality 2D MoS. Fig.1(c) shows the PL spectra of 1L MoS grown on SiC and SiO as a comparison. The PL intensity of the 1L MoS on SiC was about 500× weaker in comparison to 1L MoS grown on a 100 nm SiO/Si substrate. Though a tiny PL peaks were observed at 649 nm from the 1L MoS/SiC sample, the quenching of PL intensity forebodes a strong interaction between the as-grown MoS and SiC substrate, which agrees well with the following theoretical observations. And no PL peaks were observed from 2L, 4L and bulk samples.
Surface chemistry was probed with XPS measurements. The fitting rule adopted in this work is that the full width at half maximum (FWHM) of both components are comparable and the intensity ratios follow the expect quantum mechanically predicted ratio when XPS doublet have narrow separations. Fig.2 shows the Si 2p scan for SiC, 1L, 2L and 4L MoS
/SiC. The SiC Si 2p peak was resolvable into 2p
and 2p
peaks due to the asymmetry of the 2p peak and was done to enhance the accuracy of the subsequent band alignment calculations. The Si 2p
peak for all samples except bulk MoS
/SiC was aligned to the literature value of 100.4 eV to account for charging effects [
55]. The spectra for the bare SiC showed no signs of atmospheric oxidation which supports that the etching process removed any surface oxidation on the SiC surface. A convoluted Si
1+ peak can be observed at a binding energy of 102.07,102.16 and 102.21 eV for the 1L, 2L and 4L MoS
grown samples. This is because the XPS characterization was done ex-situ and the MoS
grown samples were left in air for a longer time for deposition whereas the bare SiC substrate was characterized immediately after etching. Nonetheless, there was only a small amount of Si
1+ present.
Fig.3 shows the Mo 3d and S 2p scan for 1L, 2L, 4L and bulk MoS
. The Mo
4+ 3d
and 3d
peaks were located at 228.69, 228.59, 228.60 and 231.83, 231.78, 231.76 eV for 1L, 2L and 4L MoS
/SiC. Whereas the S
2− 2p
and 2p
peaks were located at 161.50, 161.43, 161.42 and 162.69, 162.62, 162.62 eV respectively. This agrees well with the current literature values [
56-
59]. Due to the absence of the SiC peak in bulk MoS
, the spectrum was aligned to the Mo 3d
of 4L MoS
for presentation purposes. Due to incomplete sulfurization, small amounts of metallic Mo was detected on 2L (3.5 at %), 4L (3.0 at %) and bulk (2.0 at %) MoS
.
To derive the band alignment of the MoS
/SiC heterostructures, we adopted a direct method suggested by Santoni
et al. [
60], which requires to align the the XPS spectra of SiC, 1L, 2L and 4L MoS
to a common reference. The alignment has already been done to the Si 2p
peak at 100.4 eV. The Si 2p
peak in SiC was chosen to be the reference as it can be detected in all the required samples in sufficient quantity. The valence band maximum (VBM) was then derived by the intersection of the linear regressions of the leading edge of valence band and the baseline of its spectra, as shown in Fig.4. The VBM of bare SiC was identified to be 1.26 eV. Interestingly, the valence band spectra of 1L, 2L and 4L MoS
/SiC exhibit two leading edges with different slope. According the following theoretical density of states (DOS), the red linear regression of the leading edge comes from the Mo 4d
orbital, while the blue one originates from the Mo 4d
orbital. It is the Mo 4d
orbital that contributes the VBM of 1L MoS
. Then, the VBM of the 1L, 2L and 4L MoS
/SiC could be identified to be −0.23, −0.78, and −0.93 eV, respectively. Further considering the VBM of SiC substrate of 1.26 eV, the valence band offsets (VBOs) could be obtained to be 1.49, 2.04, and 2.19 eV for the 1L, 2L and 4L MoS
/SiC, respectively. Therefore, the VBOs increase from 1.49 eV up to 2.19 eV with increasing MoS
layers up to 4L.
To understand above experimental observations and illuminate the interfacial interactions and electronic properties in the MoS
/4H-SiC heterostructures, the first-principles calculations were performed by using the Vienna ab initio simulation package (VASP). Fig.5 shows the optimized atomic structure of 1L MoS
and the HSE06 total DOSs of 1L MoS
and partial DOSs of Mo and S atom. The lattice constants of 3.19 Å and the HSE06 bandgap of 2.10 eV are obtained for 1L MoS
. These basic parameters are in good agreement with the previous theoretical and experimental results [
12,
61]. The conduction band minimum (CBM) and VBM mainly come from the Mo atom. Careful observation shows that the CBM mainly comes from the Mo 4d
orbital. The valence band mainly originates from the 4d
and 4d
orbitals, while it’s the 4d
orbital that contributes the VBM of 1L MoS
. The different orbital origins of the valence band lead to the formation of two different slope distributions in the DOS of MoS
at the valence band. This is consistent well with the above XPS valence band spectra of MoS
, as exhibited in Fig.4, and is beneficial to the accurate identification of the VBM of 1L, 2L and 4L MoS
in XPS valence bands.
The lattice constant of bulk 4H-SiC was obtained to be 3.09 Å, consistent with the reported value [
62,
63]. The HSE06 bandgap is 3.49 eV, as shown in Fig.6(b), which is only 7 % larger than the reported bandgap of 3.26 eV. The lattice constant of the substrate 4H-SiC was slightly expanded in the model constructed to make it consistent with the lattice constants of MoS
, due to the weak vdW interlayer interaction in multilayer MoS
, which will be constructed in our calculations. The bandgap of the strained bulk 4H-SiC is 3.10 eV, as shown in Fig.6(b), which is only about 5 % smaller than the reported bandgap of 3.26 eV due to the slight tensile strain. Therefore, the model based on the strained 4H-SiC can effectively ensure the accurate results of the band alignment.
The 4H-SiC (0001) surface is usually formed by a Si or C atomic layer due to the polarity. Previous researches have demonstrated that the Si-face is energy stabilized surface as compared to the C-face [
64]. Then the geometrical and electronic structures of 4H-SiC (0001) Si surface under slight tensile strain were adopted and calculated in this study. Fig.6(a) shows the strained atomic structure of 4H-SiC (0001) Si surface. Due to the Si dangling bond formed on the surface, the HSE06 total DOSs, as shown in Fig.6(b), exhibit the strong surface states in the middle of the bandgap. And the Fermi level just passes through the surface states, indicating the Fermi level pinning, which is in agreement with the case of 2H-SiC (0001) Si surface [
59]. Fig.6(c) shows the partial DOSs of the outmost and inner Si and C atoms as indicated by the red and green circles in Fig.6(a), respectively. It is found that the mid-gap surface states mainly come from the outmost unsaturated Si atoms due to the broken Si-C bonds at the surface. It seems like that the outmost Si and C atomic layers contribute greatly to these surface states, while almost no surface states were observed in the bandgap for the inner Si and C atomic layers. Therefore, the inner Si and C atoms keep the electronic properties of bulk 4H-SiC.
Furthermore, 1L MoS/4H-SiC configuration was first constructed by a MoS monolayer (containing 1 Mo, 2 S atomic layers) and 8 atomic layers of the 4H-SiC(0001) Si surface (containing 4 Si and 4 C atomic layers). Considering the lattice matching and hexagonal symmetry, six representational high-symmetric 1L MoS/4H-SiC configurations were constructed as shown in Fig.7(a)−(f). The binding energies of these configurations were calculated by using the following equation:
where
,
and
are the total energies of MoS
/4H-SiC system, isolated 4H-SiC and MoS
, respectively. Among above six possible configurations, the preferred atomic structure of 1L MoS
/4H-SiC, as shown in Fig.7(c) and (f), has the lowest binding energy of −0.96 eV. In this energy favorable configuration, S atom locates on the top of the Si atom of the 4H-SiC (0001) Si surface, while Mo atom is aligned with the C atom below. As shown in Fig.7(f), a distinct Si-S chemical bond is formed at the 1L MoS
/4H-SiC interface. Moreover, the Si-S bond length is calculated to be 2.20 Å. This value is shorter than the interlayer distance of 2.5−3.5 Å in the vdW systems and similar with the Si-S bond length of 2.0 Å in other covalent bonding system [
65,
66]. Above large binding energy and formation of a chemical Si-S bond indicate a significantly strong interlayer interaction at the 1L MoS
/4H-SiC interface, agreement well with the above experimental observation of the quenching of PL.
To identify the interface properties, especially the band alignment, atomic DOSs of 1L MoS
/SiC interface were calculated and aligned from bottom to up according to the atomic order, as shown in Fig.8(a). As compared to that of bare 4H-SiC (0001) Si surface, the original mid-gap surface states vanish when 1L MoS
layer adsorbed on the surface. That is, the surface states originated from the Si dangling bond was passivated by the S atom by forming the chemical Si-S bond. Though the absence of the dangling bond of the MoS
surface, the existing of lone pair at the S surface accounts for the surface passivation of the SiC substrate, according to our previous calculations [
67]. Then the Fermi level shifts up and passes through the bottom of conductance band, indicating the Fermi level pinning at the interface states of the 1L MoS
/4H-SiC interface.
According to the real space atomic DOSs alignments, we could directly identify the VBM and CBM of each atomic layer and further obtained the band alignment at the 1L MoS
/SiC interface, as indicated by the red dashed lines in Fig.8. Atomic DOSs of bottom C1, Si1 and C2 exhibit the inner SiC electronic properties with the bandgap of about 3.23 eV. The CBM and VBM lie at the −0.13 and −3.36 eV, respectively. Looking up from Si2 to S2, the CBM shifts a little down to −0.33 eV. However, the VBM shifts largely up to −1.90 eV from C4 to S2. Then the bandgap shrinks from 3.23 eV to 3.03 eV, further to 1.57 eV from bottom to up in the 1L MoS
/SiC heterostructure. And the CBM and VBM exhibits special asymmetric evolution, which would play important role in the carrier distribution. The new bandgap of 1L MoS
(1.57 eV) in 1L MoS
/SiC heterostructure is much smaller than that of bare MoS
(2.10 eV), further indicating a strong interaction between MoS
and SiC. The interaction at the MoS
/SiC interface is totally different from the previous weak interaction among the WS
/HfO
, MoS
/HfO
and WS
/SiO
interfaces [
68–
70]. The CBM and VBM distributions demonstrate the formation of type I band alignment at the 1L MoS
/SiC interface. However, due to the special asymmetric evolution, the CBM locates at the SiC and MoS
region, while the VBM only lies at the MoS
region, different from the traditional type I band alignment.
Moreover, the charge transfer between 1L MoS
and SiC substrate was identified by using the Bader analysis [
45]. The calculation results show that about 0.96e transfer from the SiC substrate to the MoS
layer. Since that, an electric field forms and it directs from SiC to MoS
. The built-in electric filed is significantly important for the high performance of the electronic and optoelectronic devices based on the MoS
/SiC heterostructure [
20]. Combined with the built-in electric field and the special band alignment, the photogenerated holes tend to accumulate to the MoS
layer, while electrons tend to accumulate to the SiC layer. The efficiency of separation of photogenerated carriers is benefitting for enhancing the performance of the photodetectors and photocatalytic devices.
In order to understand the evolution of VBOs as the MoS
thickness increasing, atomic DOSs of 2L and 4L MoS
/SiC interfaces were further calculated and aligned from bottom to up according to the atomic order, as shown in Fig.8(b) and (c). According to our group’s works by using the X-ray diffraction (XRD) and TEM [
46,
47], the as-grown multilayer MoS
on different substrates tend to exhibits the 2H-phase. And our and others previous first-principles calculations [
12,
71] also demonstrated that the 2H-phase MoS
generates the lowest formation energy in the multilayer stacking systems. Therefore, 2H-phase MoS
were mainly constructed to clarify the issue of band alignment evolution in this work. For 2L MoS
/SiC interface, the Fermi level also passes through the interface states locating at the region between the lower MoS
monolayer and SiC substrate, similar with the case at 1L MoS
/SiC interface. However for the upper MoS
layer, almost no interface states were observed in the bandgap and the Fermi level just passes through the center of the band gap. Therefore, Fermi level was just pinning at the interface states in the lower 1L MoS
/SiC heterostructure. The interaction between the upper and lower MoS
layer keeps the weak vdW interaction. That is, the lower MoS
monolayer can effectively passivate the dangling bonds of the SiC surface and obtain a new weak vdW surface.
We could also identify the VBM and CBM of each atomic layer and obtain the band alignment at the 2L MoS
/SiC interface, as indicated by the red dashed lines in Fig.8(b). Similar with 1L MoS
/SiC, atomic DOSs of bottom C1, Si1 and C2 also exhibit the inner SiC electronic properties with the bandgap of about 3.23 eV. The CBM and VBM also lie at the −0.13 and −3.36 eV, respectively. Looking up from Si2 to S2, the CBM firstly shifts a little down to −0.33 eV. Then, the CBM shifts largely up to 0.67 eV from S3 to S4, further indicating the weak vdW interaction between the upper and lower MoS
layer. For the VBM, it firstly shifts largely up to −1.95 eV from C4 to S1. Then, the VBM continues to shift up to −1.14 eV from Mo1 to S4 due to the interlayer orbital coupling interaction in bilayer MoS
. The interlayer orbital coupling tends to split and shift the VBM to the higher energy region, which has been demonstrated in the previous theoretical calculations of multilayer MoS
[
12]. Then the bandgap changes from 3.23 eV to 3.03 eV, 1.62 eV, 0.81 eV and 1.81 eV from bottom to up in the 2L MoS
/SiC heterostructure. Therefore, the CBM and VBM exhibit more complicated asymmetric evolution, as compared with that of 1L MoS
/SiC. The CBM locates at the interface from Si2 to S2, while the VBM locates at the MoS
layer from Mo1 to S4. On the whole, the complicated CBM and VBM evolution demonstrates the type II band alignment of the 2L MoS
/SiC heterostructure. Compared with 1L MoS
/SiC, 2L MoS
/SiC possess thicker layers, smaller bandgap distribution, and type II band alignment, which would more effectively improve the light absorption, separate the photogenerated carriers and enhance the performance of optoelectronic devices.
For 4L MoS/SiC interface, the interface states also locate at the region from Si2 to S2, and Fermi level is pinned at them, similar with the cases in the 1L and 2L MoS/SiC interfaces. Except the lowest MoS monolayer, the upper MoS layers show almost no interface states in the bandgap and the Fermi level just passes through the center of the bandgap. Therefore, Fermi level was just pinning at the interface states in the lowest 1L MoS/SiC heterostructure. The interaction between the upper MoS layers and lowest MoS monolayer keeps the weak vdW interaction, similar with the case of 2L MoS/SiC interface.
The evolution of CBM and VBM of 4L MoS/SiC interface, as indicated by the red dashed lines in Fig.8(c), exhibits the similar trend as that of 2L MoS/SiC interface. Atomic DOSs of bottom C1, Si1 and C2 also exhibit the inner SiC electronic properties with the bandgap of about 3.23 eV. The CBM and VBM also lie at the −0.13 and −3.36 eV, respectively. Looking up from Si2 to S2, the CBM firstly shifts a little down to −0.37 eV. Then, the CBM shifts largely up to 0.48 eV from S3 to S6, indicating the weak vdW interaction between the upper and lowest MoS layer. And from S7 to S8, the CBM shifts a little down to 0.19 eV. For the VBM, it firstly shifts largely up to −1.90 eV from C4 to S1. Then, the VBM continues to shift up to −1.00 eV from Mo1 to S6 due to the interlayer orbital coupling interaction in quadlayer MoS. Finally, the VBM shifts a little down to −1.29 eV. Then the bandgap changes from 3.23 eV to 2.99 eV, 1.53 eV, 0.63 eV, 1.48 eV, and 1.48 eV from bottom to up in the 4L MoS/SiC heterostructure. The CBM locates at the interface of lowest MoS and SiC substrate from Si2 to S2, while the VBM locates at the MoS layer from Mo1 to S6. The 4L MoS/SiC heterostructure also exhibits the type II band alignment. It will have similar advantages to 2L MoS/SiC in enhancing electronic and optoelectronic device performance.
Fig.9 (a)−(c) exhibits the theoretical energy band alignments at 1L, 2L, and 4L MoS/SiC interfaces. The 1L MoS/SiC interface exhibits the type I band alignment with the CBO of 0.20 eV and VBO of 1.46 eV. The large VBO and small CBO are attributed to the strong interaction between 1L MoS and SiC substrate and the Fermi level pinning at the interface states. Due to the asymmetric CBM and VBM evolution, it is different from the traditional type I band alignment. Further due to the built-in electric field, the photogenerated electrons accumulate at SiC region, while the photogenerated holes accumulate at MoS region.
For the 2L (4L) MoS/SiC interface, they exhibit the type II band alignments and the CBO and VBO are 0.80 (0.61) and 2.22 (2.36) eV, respectively. The enlarged CBOs and VBOs, as compared to that of 1L MoS/SiC, are attributed to the weak vdW interaction between the lower MoS monolayer and upper MoS layers and the strong interlayer orbital coupling. Then the photogenerated electrons and holes tend to accumulate to the SiC and MoS layer, respectively, due to type II band alignment and built-in electric field. The high efficiency of charge separation in these heterostructures will enhance the device performance significantly. Finally, both the experimental and theoretical (2H-phase stacking) VBOs of the 1L, 2L and 4L MoS/SiC interfaces were ploted in Fig.9(d). As the MoS thickness increasing from 1L to 4L, the experimental VBO increases from 1.49 eV to 2.04 eV and 2.19 eV. Theoretically, the VBO (2H-phase stacking) increases from 1.46 eV to 2.22 eV and 2.36 eV. Therefore, the experimental and theoretical (2H-phase stacking) VBOs are in good agreement with each other, confirming the validity of the interface properties demonstrated in this study.
Since the metastable 3R-phase stacking MoS
were also observed in experiments, we further calculated the atomic DOSs of the 3R-phase stacking 2L and 4L MoS
/SiC heterostructures and the theoretical VBOs (3R-phase stacking) of the 2L and 4L MoS
/SiC interfaces were also exhibited in Fig.9(d) for comparison. The VBOs of 3R-phase stacking structure increase from 1.46 eV to 2.25 eV and 2.44 eV. As compared to that of 2H-phase, the increased theoretical VBO of 3R-phase stacking is attributed to the small shrunken bandgap and enlarged VBM of the 3R-phase MoS
, which is agreement with the previous study [
71]. Therefore, both the theoretical 2H- and 3R-phase VBOs are in good agreement with each other, and the 2H- and 3R-phase stacking structures do not influence significantly the interfacial properties, especially the band alignment, of multilayer MoS
/SiC heterostructures.
4 Conclusions
In summary, the interfacial properties of 2D/3D MoS/4H-SiC heterostructures were studied by XPS and first-principles calculations. It is found that the VBO increases from 1.49 eV to 2.19 eV with increasing MoS layers up to 4L, consistent with the first-principles calculations. For 1L MoS/SiC, a strong interlayer interaction between the 1L MoS and SiC sheet was demonstrated by the formation of a chemical Si-S bond (2.20 Å) and a large binding energy (−0.96 eV/unit cell). Then Fermi level pinning at the interface states and a passivation of 4H-SiC (0001) Si surface were realized. About 0.96e per unit cell transferring from SiC to MoS monolayer indicates an electric field forms, and it directs from SiC to MoS. Then, the 1L MoS/SiC interface exhibits the type I band alignment with the CBO of 0.20 eV and VBO of 1.46 eV. The large VBO and small CBO are attributed to the strong interaction between 1L MoS and SiC substrate and the Fermi level pinning at the interface states. For 2L and 4L MoS/SiC, Fermi level was just pinning at the lower MoS monolayer. The interaction between the upper MoS layers and lower MoS monolayer keeps the weak vdW interaction. Then, they exhibit the type II band alignments and the CBO and VBO are 0.80 (0.61) and 2.22 (2.36) eV, respectively. The enlarged CBOs and VBOs are attributed to the weak vdW interaction between the lower MoS monolayer and upper MoS layers and the strong interlayer orbital coupling. Despite of the different band alignment type, the photogenerated electrons and holes tend to accumulate to the SiC and MoS layer, respectively, due to the asymmetric band alignment and built-in electric field in all the 1L, 2L and 4L MoS/SiC interfaces. The high efficiency of charge separation in these heterostructures will enhance the device performance significantly. The multiple interfacial interactions discussed in this study provide a new modulated perspective for the next-generation electronics and optoelectronics based on the 2D/3D semiconductors heterojunctions.