Quasi 2D electronic states with high spin-polarization in centrosymmetric MoS2 bulk crystals.

Time reversal dictates that nonmagnetic, centrosymmetric crystals cannot be spin-polarized as a whole. However, it has been recently shown that the electronic structure in these crystals can in fact show regions of high spin-polarization, as long as it is probed locally in real and in reciprocal space. In this article we present the first observation of this type of compensated polarization in MoS2 bulk crystals. Using spin- and angle-resolved photoemission spectroscopy (ARPES), we directly observed a spin-polarization of more than 65% for distinct valleys in the electronic band structure. By additionally evaluating the probing depth of our method, we find that these valence band states at the point in the Brillouin zone are close to fully polarized for the individual atomic trilayers of MoS2, which is confirmed by our density functional theory calculations. Furthermore, we show that this spin-layer locking leads to the observation of highly spin-polarized bands in ARPES since these states are almost completely confined within two dimensions. Our findings prove that these highly desired properties of MoS2 can be accessed without thinning it down to the monolayer limit.

Two-dimensional materials are considered excellent candidates for next-generation electronic devices that could overcome the restrictions of classical, Si-based electronics. Apart from the possibility of a reduced size compared to state of the art transistor devices, thinning down materials like molybdenum disulphide (MoS2) gives rise to new phenomena that might allow entirely different approaches for applications such as spintronics, valleytronics, solar cells, or optical sensors1.

Time reversal symmetry dictates that electronic states in nonmagnetic, centrosymmetric crystals cannot be spin-polarized. However, a recent publication by Zhang and Liu et al.2 introduces a “hidden spin-polarization” resulting from specific site asymmetries instead of the crystal space group. MoS2 is a model material which exhibits such effects. It is a layered crystal that consists of two-dimensional S-Mo-S atomic trilayers with broken inversion symmetry within the plane. In A-B stacking these trilayers recover inversion symmetry [Fig. 1b] and therefore the bulk shows a centrosymmetric space group34. Spin-polarization in MoS2 is zero when averaged in real or reciprocal space. Nevertheless, as it was shown experimentally for other materials, electronic states resolved in space and in momentum can be spin-polarized56. Angle resolved photoemission spectroscopy (ARPES) meets both of these conditions due to its surface sensitivity and momentum-resolution capabilities and was therefore chosen a tool for the experimental study of the hidden spin-polarization in MoS2.

Figure 1

(a) Top view of a MoS2 monolayer, (b) side view of the 2H crystal structure of bulk MoS2, (c) optical reflectivity spectrum showing exciton transitions characteristic of 2H-MoS2, (d) AFM height profile along the red line in AFM image, (e) AFM amplitude image with single monolayer step edge (in the bottom left) separating two large, atomically flat areas. The smaller speckles are adsorbants from ambient air exposure, (f,g) LEED images taken at two positions on the sample ≈1 mm apart showing a reversing of the threefold symmetry.

The majority of electronic structure studies of MoS2 focuses on quasi freestanding trilayers, which are usually referred to as one monolayer of MoS2, and its integration into novel two-dimensional devices789. One of the advantageous properties of a freestanding MoS2 monolayer is its direct band gap of about 1.8 eV at the point, while bulk MoS2 has a smaller, indirect band gap with the valence band maximum (VBM) at the Γ point910. Furthermore, at the point the topmost valence band forms a distinct valley, that makes MoS2 a candidate in the novel field of valleytronics11. In MoS2, as in other transition metal dichalcogenides (TMDCs), the valley is split into two subbands. These split states are usually considered to be highly spin-polarized perpendicular to the surface in a freestanding monolayer, but unpolarized in bulk MoS231213. However, the local nature of spin-orbit coupling (SOC) leads to a locking of the spin to layers. An indication for this spin-layer locking is the high circular polarization dependence of the photoluminescence for inversion-symmetric TMDC films, e.g. consisting of an even number of layers or bulk-like crystals41114. Recently Liu et al.15 derived from first principle calculations that the observed circular polarization is a consequence of this hidden spin-polarization. Consistently with this report, our findings for bulk MoS2 and a recent publication6 concerning the similar layered crystal WSe2 show that a spin-polarization can be observed, as long as the probing tool focuses on a specific layer in the bulk. By the same mechanism, a significant spin-polarization has been measured in the bulk states of Bi5.

In this article we present a combined spin-resolved ARPES (spin-ARPES) and density functional theory (DFT) study of a cleaved MoS2 single crystal surface. We focus on the valence band around the point and in both experiment and theory. Specifically, we can show that these electronic states in the bulk material are almost completely confined within the plane of the layers and highly spin-polarized.

Results and Discussion

The spin-polarization of the states of MoS2 is hidden by its alternating sign from one layer to the next. Although low energy photoemission is sufficiently surface sensitive to probe it, an atomically flat surface (A or B layer terminated) is required for our experiment. The amplitude image of an atomic force microscope (AFM) in Fig. 1e shows the typical surface topography of a cleaved MoS2 single crystal and that the scotch tape method produces atomically flat areas at the MoS2 surface that largely exceed the maximum field of view of the AFM of 30 μm. The height profile cut in Fig. 1d reveals a step of ≈0.6 nm height, which fits a single monolayer step edge very well. The numerous small speckles that are visible within the terraces have a much lower height of ≈0.2 nm. Very likely these are adsorbants from air, since the sample has been exposed to ambient pressure during these measurements. However, we have no direct feedback for the microscopic position of the beam spot in the ARPES experiments and have to rely on the quality of the spectra as an indication for a flat area. Mahatha et al.16 reported an inhomogeneous band bending in ARPES from areas with high step densities, which we also observed in other areas of our samples.

Another consideration is concerned with the possible existence of a non-centrosymmetric 3R-MoS2 phase in our samples. To rule out the 3R-MoS2 phase as the origin of the spin-polarization that is the central focus of this article, we performed optical reflectivity measurements which are known to show distinguishable exciton transitions for 2H-MoS2 and 3R-MoS2121718. Figure 1c shows a reflectivity spectrum from our sample. The observed separation of the two exciton transitions perfectly matches the expected value of 0.18 eV for 2H-MoS2 compared to 0.14 eV for the 3R-phase. This observation is consistent with the reversal of the symmetry in the LEED images in Fig. 1f,g which corresponds to two terraces of the 2H-phase, as they were observed with the AFM, separated in height by an odd number of monolayers.

From these results, we conclude that in the mineral crystals that were used for our experiments, the centrosymmetric 2H-phase is clearly dominant. We discuss the possibility of a minor 3R-contamination in detail in the supplementary section and show that one can safely assume an upper limit of 10% 3R-MoS2 and that the effect on our concluded spin-polarization would only be 5%.

Figure 2 depicts a direct comparison of the measured ARPES data with the calculated bulk projected band structure along the direction for both 2H-MoS2 and 3R-MoS2. The same comparison is shown for the direction in Fig. 2 of the supplementary section. Although our experiment is extremely surface sensitive we find a very good agreement of the ARPES spectra (in particular at the point, where the differences between the two phases are most pronounced) with the calculated bulk band structure of 2H-MoS2, which supports the other evidences that we measured the centrosymmetric crystal phase. Finding that the bulk band structure matches ARPES spectra has been reported in several studies12161920 and is not surprising for this class of materials since the layers do not give rise to dangling bonds in the (0001) direction and show almost no surface reconstruction.

Figure 2

(a) ARPES spectrum along direction measured with hν = 21.2 eV using He-VUV source and the corresponding calculation of bulk projected band structure for 2H-MoS2 in (b) and 3R-MoS2 in (c). In the calculated maps bright areas indicate sharp bands with little out-of-plane dispersion, darker bands a broader projection. In (b,c) the Fermi energy is shifted for comparison with (a).

At the point, which is the line KH of the 3D bulk band structure projected into the surface plane, the bulk projection reveals two very sharp subbands of the topmost valence band. Although our bulk calculations are in good overall agreement with earlier studies19212223 only few studies show the out-of-plane dispersion and some show a much larger dispersion for the split valence band along the KH direction with degenerate bands at the H point1024. In Fig. 3d,e we show that this is an artifact caused by the omission of SOC effects in the calculations. The out-of-plane dispersion at the point drops drastically when SOC is included (Fig. 3d). This means that the strong SOC in MoS2 not only gives rise to a splitting of the valence band along the entire KH direction, but it also increases the confinement of these states within the monolayer plane. The resulting bands are quasi two-dimensional and therefore very similar to those reported for freestanding monolayers of MoS2 and similar TMDCs925.

Figure 3

(a–c) ARPES spectra of the VBM at the point taken at different photon energies to map the k⊥ dependence, (d,e) DFT calculation of the bulk band structure of MoS2 with and without SOC, color coded for different values of k⊥, (f,g) spin-ARPES spectra of the valence band at the / point. These spectra are corrected for the spin-detector efficiency and the non-zero emission angle. The filled area represents the statistical error, (h) 10-layer slab calculation of MoS2 band structure. The size of the red and blue circles represents the out-of-plane spin-polarization in the topmost monolayer.

In the ARPES spectrum in Fig. 2a the splitting of the valence band at the point seems obvious. Since the ARPES spectrum only cuts a 2D plane out of the 3D Brillouin Zone (BZ) the observation of the two sharp bands alone does not prove that they remain split along the entire out-of-plane KH direction. However, in photoemission the out-of-plane momentum of the photoelectron can be scanned by varying the photon energy. Figure 3a–c show the ARPES spectra along the high symmetry direction for three different photon energies. All three spectra represent a part of the 3D-BZ with a different out-of-plane photoelectron momentum k⊥. The values for k⊥ were estimated using the free electron final state model and assuming an inner potential of 10 eV. There is almost no change in the size of the splitting measured with different photon energies, which is a strong indication for a low out-of-plane dispersion and therefore an in-plane confinement as we find it in the DFT calculations. To our knowledge there are no other reports regarding an experimental observation of the out-of-plane dispersion of these states in bulk MoS2. However, our observation fits very well to recent reports that show the same behavior for WSe26 and WS220, which have a very similar band structure.

Figure 3f,g show the spin-ARPES spectra taken at the point and the point on the opposite side of the BZ. One can immediately see that these states are highly spin-polarized and that the polarization reverses at the point. The two observed peaks are separated by 180 meV, which is expected for 2H-MoS21220. This rules out the non-centrosymmetric 3R-MoS2 with an expected splitting of 140 meV as origin of the observed polarization12.

For the polarization we find ≈65% at and ≈50% at . The difference in the size of the polarization is very likely caused by a sightly different position of the beam spot on the sample. To reach the point, the sample had to be rotated in between the measurements and this could have caused a larger contribution of adjacent terraces.

An additional mechanism that effectively reduces the observed spin-polarization is photoemission from more than one monolayer, since the polarization reverses from one layer to the next. Although the low kinetic energy of the photoelectrons of ≈15–40 eV makes our method extremely surface sensitive, we can estimate the contribution of the individual layers. Using the universal curve26, we assume 4 Å as inelastic mean free path. Considering the thickness of one monolayer of MoS2 of ≈6 Å and the photoelectron emission angle of 36° for the spin-ARPES measurements, one can estimate that ≈83% of the photoemission signal comes from the first and the third monolayer and only ≈17% from the second (see supplementary section). In the case of a fully polarized state we would effectively observe only 66% spin-polarization. This simple model suggests that the 65% spin-polarization that we observed experimentally could in fact result from a polarization close to 100% in the individual layers. Our approximation for the inelastic mean free path is very rough, but a more reliable calculation is not trivial, especially for a highly anisotropic material like MoS2.

In Fig. 3h we show our results of a 10 monolayer slab calculation including the spin-polarization of the first monolayer. For the VBM at the point we calculated a spin-polarization larger than 99%. This result supports our approximation of the probing depth and the conclusion that the polarization of the individual layers is almost 100%. We want to point out that the polarization is not caused by the symmetry breaking at the surface. Also in our bulk calculations shown in Fig. 3d we calculated a spin-polarization larger than 60% in an individual layer.

So far there are only few publications to which we can compare our results for the spin-polarization. The only study of bulk MoS2 so far was published recently by Suzuki and Sakano et al.12. In this study they present a comparison of centrosymmetric 2H-MoS2 and its non-centrosymmetric polytype 3R-MoS2, which is rare in nature but Suzuki and Sakano et al. were able to grow bulk crystals of both phases. Their results seemingly contradict our findings since their spin-ARPES data shows no spin-polarization for 2H-MoS2 and close to 100% for the 3R-type. They do not directly show the microscopic topography of the 2H-MoS2 samples but an optical micrograph image suggests a much smaller domain size than we find for the exfoliated mineral crystals. Therefore it is very likely that the vanishing spin-polarization that Suzuki and Sakano et al. found is a result of the compensation by averaging over a large number of terraces or crystal grains with different orientations.

Our findings are in better agreement with the results for WSe2 that were published by Riley et al.6. They find close to 100% spin-polarization for centrosymmetric WSe2, which has the same crystal structure and therefore already proves the existence of the hidden spin-polarization in this type of material. In the context of this article the essential difference between MoS2 and WSe2 is the larger SOC mediated by tungsten, which is significantly heavier than molybdenum. The strength of the SOC directly determines the size of the valence band splitting at the point. We find a splitting of ≈180 meV for MoS2 compared with ≈500 meV that were reported by Riley et al. for WSe26. As another result of the larger spin-orbit splitting in WSe2, the topmost valence band at the point is almost degenerate with the global VBM at the point. However, a larger SOC does not necessarily result in higher spin-polarization: E.g. in topological insulators it was found that SOC-induced spin-mixing can effectively reduce the spin-polarization2728. In TMDCs the measured polarization is influenced by the localization of the states in single layers, which is larger in the W-compounds215. This can explain the higher spin-polarization observed in WSe26. Riley et al. show a dependence of the measured polarization on the photon energy. They perform a model calculation to consider the contribution of the second monolayer to the photoemission spectrum, that is in good agreement with the observed photon energy dependence. We were only able to measure the polarization at one photon energy, which prevents us from performing such analysis but it confirms our assumption that the real spin-polarization in MoS2 is at least ≈65% but possibly higher.

Conclusions

Using a combination of ARPES and DFT calculations we were able to show that the distinct valleys in the electronic band structure of bulk MoS2 are quasi two-dimensional and highly spin-polarized despite its three-dimensional, centrosymmetric space group. We have experimentally shown the confinement of these electronic states within the plane of the monolayers with the absence of out-of-plane dispersion in the photoemission spectra and confirmed these findings with our DFT calculations. The spin-polarization that we observed directly in the photoemission experiment was only 65%. However, by taking into account that the probing depth of our method is slightly higher than the thickness of one monolayer, we could estimate that in the individual layers the valleys at the K points are close to being fully polarized. This result is also in very good agreement with our DFT calculations, which show a spin-polarization larger than 99%.

Our findings represent the first observation of this recently introduced type of hidden spin-polarization in MoS2. This class of materials could broaden the playground for novel spintronic applications as long as the compensated nature of the spin-polarization, both in real and reciprocal space, is taken into account.

Methods

Experimental procedures

The sample was prepared by cleaving a MoS2 single crystal (SPI Supplies, USA) in air using the scotch tape method creating a surface with shiny and flat areas. After the cleaving, the samples were immediately put into ultra high vacuum (<3 ⋅ 10−10 mbar).

The ARPES spectra were measured at two different experimental setups. Some spectra were taken using a He discharge source and its spectral lines He I with a photon energy of hν = 21.22 eV and He II with hν = 40.81 eV. For this setup we used an MBS A-1 hemispherical analyzer and the sample temperature was kept at ≈15 K. The overall energy resolution was 20 meV.

Further ARPES spectra, including the spin-ARPES spectra, were taken at beamline 5 of DELTA at TU Dortmund. At DELTA we used a VG Scienta SES-2002 hemispherical analyzer and a Focus SPLEED detector for spin-ARPES measurements29. The spot size of the synchrotron light on the sample was ≈800 μm × 300 μm. We subtracted a Shirley background from the spin resolved spectra and the asymmetry has been corrected for the detector efficiency using a Sherman function of 0.25 and for the non-zero emission angle. These experiments were performed at room temperature with an overall energy resolution of 100 meV.

The reflectivity measurements were performed using an optical microscope with a tungsten krypton lamp light coupled to optic fibre as an illumination source and with an objective lens with a magnification 10, corresponding to the light spot size of 1 mm diameter. To compensate for the higher lateral resolution of this method compared to spin-ARPES we averaged reflectivity spectra taken at 20 positions from the area that was illuminated by the synchrotron light. The reflected light spectrum was coupled with an optic fibre to a spectrometer with a spectral resolution of 0.5 nm.

Band structure calculations

All calculations were carried out within the all-electron full-potential linearized augmented-plane-wave (FLAPW) formalism as implemented in the DFT code fleur30. The electron density was determined self-consistently employing the Perdew, Burke and Ernzerhof (PBE) parametrization of the GGA exchange-correlation functional31. The core electrons were treated fully relativistically by solving the Dirac equation with the spherically averaged effective potential around each nucleus. For the valence electrons, space is partitioned into muffin-tin (MT) spheres and an interstitial region. In the former we used an angular momentum cutoff lmax = 10 and 8 for the Mo and S atoms, respectively. In the latter, a plane-wave cutoff of 4.2 bohr−1 was used for the bulk calculations. In the MT spheres, relativistic effects are included in the scalar-relativistic approximation32 while the spin-orbit coupling (SOC) is incorporated self-consistently employing the “second variation” technique33. An 8 × 8 × 4 k-point grid was used to sample the bulk BZ. For the film calculations somewhat smaller cutoff (3.8 bohr−1) and a 7 × 7 k-point grid were used. We used the experimental lattice structures from a publication by Dickinson et al.34.

Additional Information

How to cite this article: Gehlmann, M. et al. Quasi 2D electronic states with high spin-polarization in centrosymmetric MoS2 bulk crystals. Sci. Rep. 6, 26197; doi: 10.1038/srep26197 (2016).