Observation of monolayer valence band spin-orbit

Full text

Turn on search term navigation

ARTICLE

Received 30 Dec 2013 | Accepted 11 Jul 2014 | Published 22 August 2014

Nasser Alidoust1, Guang Bian1, Su-Yang Xu1, Raman Sankar2, Madhab Neupane1, Chang Liu1, Ilya Belopolski1, Dong-Xia Qu3, Jonathan D. Denlinger4, Fang-Cheng Chou2 & M. Zahid Hasan1

MoS2 grown on highly ordered pyrolytic graphite substrate. Our experimental results resolve the Fermi surface trigonal warping of bulk MoSe2, and provide evidence for the critically important spin-orbit split valence bands of monolayer MoS2. Moreover, we systematically image the formation of quantum well states on the surfaces of these materials, and present a theoretical model to account for these experimental observations. Our ndings provide important insights into future applications of transition metal dichalcogenides in nanoelectronics, spintronics and photonics devices as they critically depend on the spin-orbit physics of these materials.

DOI: 10.1038/ncomms5673

Observation of monolayer valence band spin-orbit effect and induced quantum well states in MoX2

1 Joseph Henry Laboratory, Department of Physics, Princeton University, Princeton, New Jersey 08544, USA. 2 Center for Condensed Matter Sciences, National Taiwan University, Taipei 10617, Taiwan. 3 Condensed Matter and Materials Division, Lawrence Livermore National Laboratory, Livermore, California 94550, USA. 4 Advanced Light Source, Lawrence Berkeley National Laboratory, Berkeley, California 94305, USA. Correspondence and requests for materials should be addressed to M.Z.H. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 5:4673 | DOI: 10.1038/ncomms5673 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms5673

In the recent years, with the discovery of graphene1,2 and topological insulators38, there has been an outburst of research activities to understand the physical properties of

two-dimensional (2D) materials and their applications to optoelectronics and spintronics. One family of 2D materials that has recently attracted much attention is transition metal dichalcogenides (TMDCs) such as monolayers and few-layers of MX2 (M Mo, W and X S, Se)9,10. These materials are easy to

exfoliate and fabricate11,12, and have band gaps comparable in size to that of silicon (1.11.4 eV in bulk and 1.51.9 eV in monolayers) with the predicted property of transitioning from an indirect to a direct one upon reduction in size from bulk to monolayer9,1319, making them remarkably suitable for potential use in electronics applications such as thin lm transistors20. More excitingly, it has been predicted that the exotic spin-valley coupling phenomenon occurs in monolayers of TMDCs, thus making them appropriate for spintronics and valleytronics devices2124. Nevertheless, even though many recent studies have examined the possibility of using TMDCs in single-layer transistors2528, and have investigated their optical properties and potential use in optoelectronics18,19,23,2932, surprisingly there has only been extremely limited experimental work on uncovering the key aspects of the electronic band structures of bulk, few-layers and monolayers of these materials and their modication under external effects such as doping, surface deposition or heterostructuring. The few existing momentum-resolved spectroscopic studies on these materials also lack the resolution and depth to highlight the important features of their band structures3337. For instance, only recently, using micro-angle-resolved photoemission spectroscopy (ARPES), the evolution of the valence band at the corner of the Brillouin zone (BZ) from higher to lower binding energies upon reducing the number of layers to monolayer has been demonstrated37.

The interest in monolayers of the MX2 series lies in the emergence of the direct band gap as a consequence of the disappearance of van der Waals interlayer interactions in the monolayer limit9,19, as well as the valence band strong spin-orbit splitting because of the loss of inversion symmetry9,21,22,37. However, as lms of TMDCs on typical substrates such as Si/SiO2 grow in very small akes of B10 10 mm2 and are usually

composed of islands of variable thicknesses next to each other, tracking these two properties in a momentum-resolved manner has proven to be challenging. Furthermore, the splitting of the valence band top into two distinct degenerate bands in the bulk limit due to both strong spin-orbit coupling and interlayer hopping38 could be of interest and has not been resolved to date. Moreover, until very recently MoS2 has experimentally gotten the most attention among various members of the TMDCs family, such as MoSe2, WS2, and WSe2. This is despite the fact that the stronger spin-orbit coupling of the latter compounds make them more attractive for spintronics applications.

Thus, in this work, we utilize ARPES to directly map the electronic band structure of bulk MoSe2 across the entire BZ. We then study the evolution of its band structure upon in situ surface deposition with an alkali metal (potassium). Our results demonstrate the formation of a nearly free 2D electron gas (2DEG) within the potassium overlayers conned to the surface of MoSe2 in the form of robust quantum well states (QWSs). This quantum connement of potassium layers could potentially enhance the optoelectronics performance of bulk and thin lms of the TMDCs and be utilized in optoelectronics applications of these materials. We also resolve the band structure of a monolayer of MoS2 grown on highly ordered pyrolytic graphite (HOPG) substrate using the chemical vapour deposition (CVD) growth method and nd evidence suggesting the existence of the split valence band induced by strong spin-orbit coupling and

broken inversion symmetry. We compare our ARPES results with our rst-principles theoretical calculations to gain microscopic insights on the electron behaviour in these materials, which is critically important for optimizing their potential performance in device settings.

ResultsCrystal structure of MX2. The crystal structure of MX2 (M Mo,

W and X S, Se) is a layered structure with layers of M atoms

sandwiched between two layers of X atoms. These layers, which are held together through van der Waals forces, are shifted relative to each other in such a way that M atoms in one layer are placed directly above the X atoms in the neighbouring layers9,10,34. This crystal structure and the corresponding BZ are shown in Fig. 1b,c, respectively. The centre of the BZ is the G point, with M points the midpoints of the hexagonal BZ edges, and K points the corners of the hexagon. G, M and K are the projection of these points onto the 2D surface BZ as illustrated in Fig. 1c. X-ray photoemission spectroscopy measurements of the single crystals of MoSe2 used in our studies reveal sharp peaks corresponding to Mo 3s, 3p and 3d core levels, as well as those of Se 3p and 3d as shown in Fig. 1a, thus indicating the high quality of the studied samples.

Electronic structure of bulk MoSe2. We start by investigating the band structure of MoSe2. The kE maps of the electronic band structure measured by ARPES (left panels in Fig. 2a,b) along the two high symmetry directions of the BZ, G M and G K, show

extraordinary resemblance to the calculated bands from the rst-principles calculations (right panels in Fig. 2a,b). The ARPES measurements here were performed at a photon energy of hn 90 eV. The kinetic energy value of the Fermi level is deter

mined by measuring that of the polycrystalline gold separately at the same photon energies and under the same experimental conditions. The valence band maximum (VBM) is found to be located at the G point of the BZ at a binding energy of 1.2 eV. The valence band just below the topmost band at binding energies of B1.71.9 eV, which appears clearly in our rst-principle calculations, has a somewhat weaker intensity in our ARPES measurements because of matrix element effects.

One important feature of our measurements is the unambiguous observation of the two distinct bands near the top of the valence band along the G K direction of the BZ. These bands

are highlighted in Fig. 2c, a zoomed-in version of the top part of the valence band along the G K direction, and are well in

agreement with band structure calculations conducted here as well as the ones reported in earlier studies9,14,15 (see Supplementary Note 1 and Supplementary Figs 13 for more details). Constant binding energy contours at representative binding energies of 1.4, 1.7 and 2.0 eV are presented in Fig. 2df. At EB 1.7 eV, there exists only one electron-like pocket around

each of the K points of the BZ (see Fig. 2e), whereas two concentric pockets are observed at EB 2.0 eV (see Fig. 2f),

clearly resolving the two distinct bands at the BZ corners. Moreover, the hole pockets near K exhibit trigonal warping effect as seen in Fig. 2e,f. This effect has been theoretically predicted in monolayer TMDCs39,40, and here we have shown that it exists in the bulk materials as well.

Spin-orbit-induced splitting of the valence band in monolayer TMDCs. Figure 1d shows our rst-principles band structure calculations for bulk (top row) and monolayer (bottom row) MoS2, with (right column) and without (left column) considering spin-orbit coupling. These demonstrate the indirect to direct band gap transition on going to the monolayer limit and the

2 NATURE COMMUNICATIONS | 5:4673 | DOI: 10.1038/ncomms5673 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms5673 ARTICLE

M M

K M M

Intensity (a.u.)

25,000

20,000

15,000

10,000

5,000

Mo 3s

Mo 3p

C 1s

Mo 3d

Se 3p Se 3d

No SOC SOC

Bulk

E B(eV)

1,000 800 600 400 200 0

Binding Energy (eV)

K K

Mo, W

S, Se

C1 C2

No SOC SOC

1 ML

M K

Figure 1 | Electronic structure of transition metal dichalcogenides (TMDCs). (a) X-ray photoemission spectroscopy measurement of the studied MoSe2 crystals, showing both Mo and Se peaks. The inset shows a photograph of these crystals. (b) Three dimensional crystal structure of layered TMDCs, with the chalcogen atoms (X) in orange and the metal atoms (M) in purple. (c) The corresponding bulk and surface Brillouin zone (BZ) with thehigh symmetry points marked. (d) First-principles band structure calculations for bulk (top row) and monolayer (bottom row) MoS2 with (right column)

and without (left column) considering spin-orbit coupling, demonstrating the indirect to direct band gap transition on going to the monolayer limit and the spin-orbit splitting of the monolayer valence band. The Fermi level is indicated by the dashed lines.

spin-orbit splitting of the monolayer valence band in TMDCs. The Fermi level in these kE maps is indicated by the dashed lines. These results are in good agreement with the previously reported rst-principles calculations on this material21,41.

In the next step, we turn our attention to monolayer MoS2 on

HOPG substrate. This conducting substrate was deliberately chosen to avoid the electronic charging of the lm during the photoemission process. The kE map of the electronic band structure of this sample along the high symmetry direction G K

presented in Fig. 3a, reveals the electronic bands of monolayer MoS2 as well as the linear p bands of graphite from the HOPG substrate at the K point of the BZ. We have conducted rst-principles calculations of the hybrid structure of monolayer MoS2 and graphite (Fig. 3b), which indicates the general agreement with our ARPES ndings. Despite the existence of the graphites p bands at the K points, and their merging with the MoS2 valence band around these points of the BZ, we can still resolve a large section of the MoS2 valence band top near K along the high symmetry direction G K, as highlighted in Fig. 3c, a zoomed-in

version of Fig. 3a near the top of the valence band.

Moreover, we note the broadening of the valence band top along the G K high symmetry direction (see Fig. 3c). Perform

ing 2D curvature analysis42 of the obtained bands, we resolve two split branches of the valence band of MoS2 along this high symmetry direction (close to the K point) as shown in Fig. 3d. The green dotted lines are guides to the eyes for these two spin-orbit split branches of the valence band. This is consistent with the spin splitting scenario predicted to exist in this material, as shown by our rst-principles calculations in Fig. 3e. Monolayers of TMDCs lack inversion symmetry and are noncentrosymmetric. The in-plane connement of electron motion and high mass of the elements in these materials result in strong spin-orbit splitting of the valence band. This property makes

these materials promising candidates for applications in spintronics devices9,21,22,24. Our observation here is the rst momentum-resolved experimental evidence for such spin-orbit splitting. Further, spin-resolved ARPES studies are needed to identify the spin conguration of these split bands.

To further demonstrate the orbital property and spin texture of monolayer MoS2, we calculate the charge densities and spin expectation values of the lowest two conduction states and the topmost two valence states at K marked in Fig. 1d as C1, C2, V1 and V2, respectively. The charge density plots shown in Fig. 3f indicate that the dominant orbital components are Mo 4d3z2 r

and Mo 4dx2 y

2 4dxy for C1(C2) and V1(V2), respectively. The calculated spin expectation values of states near K show nearly out-of-plane spin polarization, which can be attributed to an interplay of the pure in-plane electron motion and the inversion symmetry breaking crystal potential. The out-of-plane character of Mo 4d3z2 r

2 explains the negligible spin splitting between the C1 and C2 states.

On the other hand, the calculated energy splitting between V1 and V2 states is 138 meV at K, as a consequence of the prominent in-plane charge density variation of Mo 4dx2 y

2 4dxy orbital, which is in good agreement with previous results from density functional theory (DFT) calculations21,41. According to the calculated spin expectation values, we deploy a straightforward method to determine the spin orientation of the V1 and V2 states, which are important when considering the selection rule of polarized photon excitations. Viewing monolayer MoS2 from top, as sulphur atoms sit on the hexagonal lattices it is alway possible to nd a primitive vector connecting to adjacent sulphur atoms in the same direction as G K. Calling this vector a and drawing

another vector b from the sulphur atom at the origin of a to the nearest Mo atom, the spin orientation of V2 is in the direction of a b and that of V1 points to the opposite direction (Fig. 3g).

NATURE COMMUNICATIONS | 5:4673 | DOI: 10.1038/ncomms5673 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms5673

K K

Experiment Calculation Experiment Calculation

High

Low

E B(eV)

1.0 0.5 0.0 0.5 1.0 1.0 0.5 0.0 0.5 1.0 1.0 0.5 0.0 0.5 1.0 1.0 0.5 0.0 0.5 1.0

k|| (1)

k|| (1) k|| (1) k|| (1)

c d e f

EB=1.4 eV EB=1.7 eV EB=2.0 eV

Cut 1 Cut 2 Cut 3

1.21.41.61.82.02.22.4

Cut 1

Cut 2

Cut 3

k ||(1)

1.51.00.50.0 0.5 1.0 1.5

0.5 1.0 0.5 0.0 0.5 1.0 1.0 0.5 0.0 0.5 1.0 1.0 0.5 0.0 0.5 1.0

0.5 1.5

k|| (1) k|| (1)

Figure 2 | Electronic structure and spin-degenerate bands of bulk MoSe2. (a) Angle-resolved photoemission spectroscopy (ARPES) electronic structure measurements (left panel) and rst-principles calculated electronic structure (right panel) of MoSe2 crystals along the high symmetry direction

G M. (b) Same as a for measurements and calculations along the high symmetry direction G K. (c) Zoomed-in version of the ARPES electronic

structure spectra along the G K direction near the top of the valence band, highlighting the two distinct spin-degenerate bands at the K point of the BZ.

(df) Constant binding energy contours at EB 1.4, 1.7 and 2.0 eV, which are indicated in panel c. The green hexagon in d represents the rst BZ.

The two distinct bands can also be distinguished in the contour at EB 2.0 eV in f from the electron-like pockets around the K points. The ARPES

measurements were performed at photon energy of hn 90 eV and temperature of T 20 K.

Observation of QWSs of a nearly free 2DEG by surface deposition. To study the effect of alkali metal deposition on the surface of bulk MoSe2, we deposit potassium (K) in situ on cleaved surfaces of crystals of this material, and track the evolution of the band structure spectra as the deposition occurs. Deposition was performed at T 60 K with an estimated

deposition rate of B0.080.09 min 1. The total deposition time is 1 h and the nominal thickness of potassium is about4.85.4 . Figure 4a shows the evolution of the band structure k E map along the G M direction of the BZ during the deposition

process. We observe that the VBM at G moves closer to the Fermi level by about 200 meV. We associate this shift with the formation of a conducting layer on the surface of our sample after potassium deposition. More specically, the VBM of the MoSe2 sample is at

EB 1.0 eV, but before deposition is shifted down by about

200 meV because of the electronic charging effect (a photon-induced effect, not intrinsic to the material under study, but rather caused by the photoemission process43). Next, we observe the formation of new electronic states in the band gap of MoSe2,

extending from the Fermi level to the valence band (see middle panels in Fig. 4a). Upon the completion of the deposition process these in-gap states are completely formed as shown in the last panel of Fig. 4a and b.

We note the formation of electron-like pockets on the Fermi surface after the completion of the deposition process. These pockets appear right at the K points of the BZ. As the energy spacing between these pockets and the VBM (0.8 eV) is smaller than the indirect gap of MoSe2 (B1.1 eV), we rule out the

possibility of these pockets being a part of the conduction band of MoSe2. The resulting Fermi surface after potassium deposition shows two concentric circular-shaped pockets around the centre of the BZ, as well as the aforementioned pockets at the K points (Fig. 4d). The circular pockets around G are identied as the rst and second QWSs, whereas the pockets at K are possibly resulting from the surface states of potassium. We furthermore conduct rst-principles calculations of potassium layers formed on the surface of bulk MoSe2, as shown in Fig. 4c. In the simulation, the relaxed interlayer spacing of potassium is determined to be3.51 , and so the nominal thickness corresponds to 1.351.53 layers of freestanding K lms. As an independent check of this value, we compare the ARPES spectra with DFT simulations, which gives a thickness between 1 and 2 layers, in agreement with the previous estimate (for more details see Supplementary Note 1 and Supplementary Fig. 4). The possibility of directly accessing the conduction band of bulk and monolayer TMDCs through intercalation and adsorption of alkali metals has been theoretically predicted44. Our results suggest that adsorption of alkali metals on the surface of bulk MoSe2 results in the formation of

QWSs and a nearly free 2DEG rather than electron doping and providing access to the conduction band.

DiscussionThe splitting of the valence band along the G K high symmetry

direction of the monolayer MoS2 and other TMDCs, which is due to strong spin-orbit coupling (because of the high mass of the

4 NATURE COMMUNICATIONS | 5:4673 | DOI: 10.1038/ncomms5673 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms5673 ARTICLE

ML MoS2 & graphite

calculation

V1,V2 : 4dMo (4dMo)

x2y2 xy

Experiment

E B(eV)

1.0 0.5 0.5 1.0 1.0 0.5 0.0 0.5 1.0

1.5

0.0

k|| (1) k|| (1)

k|| (1)

1.5

2.0

2.5

1.5

2.0

2.5

1.5

2.0

0.5 0.0 0.5 1.0 1.5

C1,C2 : 4dMo

3z 2r 2

2D Curvature

S Mo

sv1 || a

2.5

Figure 3 | Spin-orbit-induced splitting of the valence band in monolayer MoS2. (a) ARPES electronic structure measurements of monolayer MoS2 on highly ordered pyrolytic graphite (HOPG) substrate along the high symmetry direction G K, showing the observation of MoS2 monolayer electronic bands

as well as the p bands of graphite at the K point of the BZ. (b) First-principles calculated electronic structure measurements of monolayer MoS2 accompanied by graphite along the high symmetry direction G K. (c) Zoomed-in version of the ARPES electronic structure spectra along the G K

direction. (d) 2D curvature plot of the bands near the valence band top, showing the experimental observation of the predicted spin-orbit split valence band along G K. The green dotted lines are guides to the eyes for these two spin-orbit split branches of the valence band. Here the ARPES measurements were

performed at photon energy of hn 110 eV and temperature of T 200 K. (e) First-principles calculated bands near the valence band top, conrming the

experimental observation of the predicted spin-orbit split valence band. (f) Charge densities of conduction and valence states at K. (g) Top view of monolayer MoS2. The spin polarization of the V1 band is in the direction of a b, provided that a is aligned along the G K direction.

elements) and lack of inversion symmetry, has been suggested in previous theoretical studies as well as our rst-principles calculations here9,21. For the bulk, this splitting is caused by the combination of interlayer interaction and spin-orbit coupling. In this case, there is no spin splitting between these two bands and they are indeed spin-degenerate (Kramers degeneracy) because of the presence of inversion symmetry along with time reversal symmetry. But, the strong spin-orbit coupling in TMDCs results in enhanced splitting of these spin-degenerate bands at the K point of the BZ (see Fig. 1d)38. Our results resolve these two distinct bands at the corners of the BZ through ARPES measurements and theoretical rst-principles calculations. Moreover, we have resolved the trigonal warping of the valence band near the K point of the BZ in the bulk materials. This effect has been theoretically predicted in monolayer TMDCs and has been attributed to an anisotropic correction of order q3 to the energy. This warping effect reects the threefold rotational symmetry of the underlying crystal lattice39,40. From symmetry considerations, the same is expected to happen for bulk TMDCs as well, as we have shown in our measurements.

Furthermore, we have provided the rst momentum-resolved experimental evidence of the spin-split valence band of mono-layer MoS2 grown on HOPG substrate. However, further spin-resolved ARPES measurements are needed to conrm the spin conguration of these two split bands. According to our calculations, the maximal splitting occurs at K, and is about 148 meV. The magnitude of the splitting depends on the momentum. The splitting curve rises sharply near k 0.6 1,

and approaches to its maximum with a smaller slope at larger k. This momentum dependence of spin splitting can be attributed to the changes in the orbital components of the states in the valence and conduction bands (see Supplementary Note 1 and Supplementary Figs 5 and 6 for details).

Recently, experimental evidence and theoretical prediction of valley connement has been suggested in monolayer MoS2.

Specically, it has been suggested that the split valence band at neighbouring K points of the BZ (K and K0) posses opposite out-of-plane spin congurations2224,32. The spin-orbit splitting observed in our study of monolayers of MoS2 paves the way to more detailed studies of spin and valley physics in TMDCs.

NATURE COMMUNICATIONS | 5:4673 | DOI: 10.1038/ncomms5673 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms5673

Time evolution of the low-energy spectra during K deposition

0.0

0.5

1.0

1.5

E B(eV)

0.8 0.4 0.8 0.4 0.0 0.4 0.80.8 0.4 0.0 0.4 0.80.8 0.4 0.0 0.4 0.8

0.0 0.4 0.8

k|| (1)

k|| (1) k|| (1) k|| (1)

Experiment

Calculation

1 ML of K 2 ML of K

0.0

0.5

1.0

1.5

2.0

2.5

3.0

E B(eV)

0.0

0.5

1.0

1.5

2.0

2.5

3.0

1.0 0.5 0.0 0.5 1.0

k|| (1)

1.0 0.5 0.0 0.5 1.0

k|| (1)

1.0 0.5 0.0 0.5 1.0

d e

EB = 0 eV

1.0

0.5

0.0

0.5

k ||(1)

1.0 0.5 0.0 0.5 1.0

MoSe2

k|| (1)

Figure 4 | Quantum well states (QWSs) on TMDCs through surface potassium deposition. (a) Evolution of the low-energy electronic structure of MoSe2 during in situ surface potassium deposition as measured by ARPES, showing the process of formation of in-gap quantum well states. (b) ARPESband structure spectra upon the completion of the potassium deposition process. (c) First-principles calculated band structure of one and two atomic layers of potassium on the surface of bulk MoSe2, indicating the formation of two QWSs in the ARPES measurements. (d) The Fermi surface obtained after the potassium deposition process is completed, featuring the observation of two QWSs in the form of two electron-like circular pockets around G,as well as the surface states of potassium at K points. The green lines represent the G M and G K directions of the BZ. Here the ARPES measurements

were performed at photon energy of hn 70 eV and temperature of T 60 K. (e) Illustration of the in situ surface potassium deposition.

It also serves as a guide for spintronics and valleytronics applications of TMDCs.

Our results also track the evolution of the band structure of MoSe2 throughout the surface deposition process with an alkali metal (potassium). We note that short after the deposition process has started potassium atoms sit on the surface of MoSe2 in a manner that results to the partial formation of the rst QWS (see second panel in Fig. 4a). Consequently, the rst QWS appears in full, as shown in the third panel of Fig. 4a, and nally the second quantum well band is formed partially close to the Fermi level. The QWSs in the epitaxial potassium layers can be considered as being conned by vacuum barrier on one side and the bulk band gap of MoSe2 on the other side45. When the energy of the QWS is within the bulk band gap of MoSe2, the state is well

trapped in the epitaxial layer (for more details see Supplementary Note 1). In the process of surface deposition, 2D conduction channels (QWSs) emerge on the surface of the semiconducting MoSe2 (with an indirect insulating energy gap of B1.8 eV) and, as a consequence, the electronic behaviour of the system changes dramatically. This tunable feature can be utilized to quantitatively control the electronic properties of this quasi-2D material. Our investigation of the dynamical process in which a nearly free 2DEG is formed on TMDCs could have important implications for future optoelectronics studies and applications of TMDCs using surface manipulations by elements such as alkali metals on bulk, few-layers and monolayers of these materials.

In conclusion, we have conducted detailed high-resolution ARPES and rst-principles calculations studies of the electronic

6 NATURE COMMUNICATIONS | 5:4673 | DOI: 10.1038/ncomms5673 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms5673 ARTICLE

structure of both bulk and monolayer TMDCs. We have resolved

the two distinct spin-degenerate valence bands of the bulk MoSe2,

and have provided evidence for the spin-orbit split bands of the monolayer MoS2, at the corners of the BZ. We have furthermore investigated the process of the formation of QWSs on the surface of bulk MoSe2, and demonstrated the formation of a 2DEG on the cleaved surfaces of TMDCs. These results will be extremely valuable in designing, fabricating and optimizing these materials for future device applications.

Finally, we note that after the submission of this manuscript, another group reported the splitting in the valence band of MoSe2 monolayer and few-layers lms, but neither trigonal warping effect nor QWSs on the bulk MoSe2 were reported there46.

Methods

Electronic structure measurements. ARPES measurements were performed with incident photon energies of 35100 eV at beamline 4.0.3 of the Advanced Light Source in the Lawrence Berkeley National Laboratory. Samples were cleaved in situ between 10 and 20 K at chamber pressure better than 5 10 11 Torr resulting

in shiny surfaces. Energy resolution was better than 15 meV and momentum resolution was better than 1% of the surface BZ.

Sample growth and characterization. To prepare MoS2 lms, molybdenum chloride (MoCl5) and sulphur powder were used as precursors and the growth of monolayer lms was carried out in a CVD chamber. The ultrasonicated cleaned HOPG substrates were loaded into the CVD chamber 6 cm away from the central hot zone. Argon (Ar) gas was used as carrier gas in a pressure of 3 Torr (ref. 47). Single crystals of MoSe2 were grown by the chemical vapour transport method, using I2 as the transporting agent. They were formed in silver-coloured, graphite-like, single-crystalline platelets up to 10 10 mm2 in surface area and 2 mm in

thickness (inset in Fig. 1a). X-ray photoelectron spectroscopy measurements were performed by means of Physical Electronics 5500 photoelectron spectrometer. All spectra were collected using Mg Ka X-ray operating at 300 W. The binding energy scale was calibrated by setting the main component of the C 1s peakat 284.6 eV.

First-principles calculation methods. First-principles calculations of the electronic structures were performed using HGH-type pseudo-potentials48 and a plane-wave basis set. The main programme employed was developed by the ABINIT group49,50. Spin-orbit coupling was included, where appropriate, using the relativistic local-density approximation (LDA). The experimental resolution was taken into account in the spectral simulations of Figs 24, but the matrix element effect was not considered (see Supplementary Note 2 and Supplementary Figs 7 and 8 for more details).

References

1. Geim, A. K. & Novoselov, K. S. The rise of graphene. Nat. Mater. 6, 183191 (2007).

2. Geim, A. K. & Kim, P. Carbon wonderland. Sci. Am. 298, 9097 (2008).3. Hasan, M. Z. & Kane, C. L. Colloquium: topological insulators. Rev. Mod. Phys. 82, 30453067 (2010).

4. Qi, X.-L. & Zhang, S.-C. Topological insulators and superconductors. Rev. Mod. Phys. 83, 10571110 (2011).

5. Hsieh, D. et al. A topological Dirac insulator in a quantum spin Hall phase. Nature 452, 970974 (2008).

6. Xia, Y. et al. Observation of a large-gap topological-insulator class with a single Dirac cone on the surface. Nat. Phys. 5, 398402 (2009).

7. Xu, S.-Y. et al. Topological phase transition and texture inversion in a tunable topological insulator. Science 332, 560564 (2011).

8. Xu, S.-Y. et al. Observation of a topological crystalline insulator phase and topological phase transition in Pb1 xSnxTe. Nat. Commun. 3, 1192 (2012).

9. Wang, Q. H., Kalantar-Zadeh, K., Kis, A., Coleman, J. N. & Strano, M. S. Electronics and optoelectronics of two-dimensional transition metal dichalcogenides. Nat. Nanotechnol. 7, 699712 (2012).

10. Balendhran, S. et al. Two-dimensional molybdenum trioxide and dichalcogenides. Adv. Funct. Mater. 23, 39523970 (2013).

11. Coleman, J. N. et al. Two-dimensional nanosheets produced by liquid exfoliation of layered materials. Science 331, 568571 (2011).

12. Castellanos-Gomez, A. et al. Laser-thinning of MoS2: on demand generation of a single-layer semiconductor. Nano Lett. 12, 31873192 (2012).

13. Kobayashi, K. & Yamauchi, J. Electronic structure and scanning-tunneling-microscopy image of molybdenum dichalcogenide surfaces. Phys. Rev. B 51, 1708517095 (1995).

14. Lebgue, S. & Eriksson, O. Electronic structure of two-dimensional crystals from ab initio theory. Phys. Rev. B 79, 115409 (2009).

15. Ding, Y. et al. First principles study of structural, vibrational and electronic properties of graphene-like MM2 (M Mo, Nb, W, Ta; X S, Se, Te)

monolayers. Physica B Cond. Matt. 406, 22542260 (2011).16. Coehoorn, R. et al. Electronic structure of MoSe2, MoS2 and WSe2. I.

Band structure calculations and photoelectron spectroscopy. Phys. Rev. B 35, 61956202 (1987).

17. Kuc, A., Zibouche, N. & Heine, T. Inuence of quantum connement on the electronic structure of the transition metal sulde TS2. Phys. Rev. B 83, 245213 (2011).

18. 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).

19. Splendiani, A. et al. Emerging photoluminescence in monolayer MoS2. Nano Lett. 10, 12711275 (2010).

20. Podzorov, V., Gershenson, M. E., Kloc, C., Zeis, R. & Bucher, E. High-mobility eld-effect transistors based on transition metal dichalcogenides. Appl. Phys. Lett. 84, 33013303 (2004).

21. Zhu, Z. Y., Cheng, Y. C. & Schwingenschlgl, U. Giant spin-orbit-induced spin splitting in two-dimensional transition-metal dichalcogenide semiconductors. Phys. Rev. B 84, 153402 (2011).

22. Xiao, D., Liu, G.-B., Feng, W., Xu, X. & Yao, W. Coupled spin and valley physics in monolayers of MoS2 and other group-VI dichalcogenides. Phys. Rev.

Lett. 108, 196802 (2012).23. Zeng, H., Dai, J., Yao, W., Xiao, D. & Cui, X. Valley polarization in MoS2 monolayers by optical pumping. Nat. Nanotechnol. 7, 490493 (2012).

24. Zeng, H. et al. Optical signature of symmetry variations and spin-valley coupling in atomically thin tungsten dichalcogenides. Sci. Rep. 3, 1608 (2013).

25. Radisavljevic, B., Radenovic, A., Brivio, J., Giacometti, V. & Kis, A. Single-layer MoS2 transistors. Nat. Nanotechnol. 6, 147150 (2011).

26. Yoon, Y., Ganapathi, K. & Salahuddin, S. How good can monolayer MoS2 transistors be? Nano Lett. 11, 37683773 (2011).

27. Yin, Z. et al. Single-layer MoS2 phototransistors. ACS Nano 6, 7480 (2012).28. Zhang, Y., Ye, J., Matsuhashi, Y. & Iwasa, Y. Ambipolar MoS2 thin ake transistors. Nano Lett. 12, 11361140 (2012).

29. Frindt, R. F. The optical properties of single crystals of WSe2 and MoTe2.J. Phys. Chem. Solids 24, 11071108 (1963).30. Frindt, R. F. & Yoffe, A. D. Physical properties of layer structures: optical properties and photoconductivity of thin crystals of molybdenum disulphide. Proc. R. Soc. Lond Ser. A 273, 6983 (1963).

31. Kam, K. K. & Parkinson, B. A. Detailed photocurrent spectroscopy of the semiconducting group VIB transition metal dichalcogenides. J. Phys. Chem. A 86, 463467 (1982).

32. Mak, K. F., He, K., Shan, J. & Heinz, T. F. Control of valley polarization in monolayer MoS2 by optical helicity. Nat. Nanotechnol. 7, 494498 (2012).

33. Bker, Th. et al. Band structure of MoS2, MoSe2 and a-MoTe2: angle-resolved photoelectron spectroscopy and ab initio calculations. Phys. Rev. B 64, 235305 (2001).

34. Mahatha, S. K., Patel, K. D. & Menon, K. S. R. Electronic structure investigation of MoS2 and MoSe2 using angle-resolved photoemission spectroscopy andab initio band structure studies. J. Phys. Condens. Matter 24, 475504 (2012).

35. Mahatha, S. K. & Menon, K. S. R. Inhomogeneous band bending on MoS2(0001) arising from surface steps and dislocations. J. Phys. Condens.

Matter 24, 305502 (2012).36. Mahatha, S. K. & Menon, K. S. R. Quantum well states in Ag thin lms on MoS2(0001) surfaces. J. Phys.: Condens. Matter 25, 115501 (2013).

37. Jin, W. et al. Direct measurement of the thickness-dependent electronic band structure of MoS2 using angle-resolved photoemission spectroscopy. Phys. Rev.

Lett. 111, 106801 (2013).38. Zahid, F., Liu, L., Zhu, Y., Wang, J. & Guo, H. A generic tight-binding model for monolayer, bilayer and bulk MoS2. AIP Adv. 3, 052111 (2013).

39. Kormnyos, A. et al. Monolayer MoS2: trigonal warping, the l valley, and spin-orbit coupling effects. Phys. Rev. B 88, 045416 (2013).

40. Rostami, H., Moghaddam, A. G. & Asgari, R. Effective lattice Hamiltonian for monolayer MoS2: tailoring electronic structure with perpendicular electric and magnetic elds. Phys. Rev. B 88, 085440 (2013).

41. Chang, C.-H., Fan, X., Lin, S.-H. & Kuo, J.-L. Orbital analysis of electronic structure and phonon dispersion in MoS2, MoSe2, WS2, and WSe2 monolayers under strain. Phys. Rev. B 88, 195420 (2013).

42. Zhang, P. et al. A precise method for visualizing dispersive features in image plots. Rev. Sci. Instrum 82, 043712043712 (2011).

43. Hfner, S. Photoelectron Spectroscopy: Principles and Applications (Springer, Heidelberg, 2003).

44. Dolui, K., Rungger, I., Das Pemmaraju, C. & Sanvito, S. Possible doping strategies for MoS2 monolayers: an ab initio study. Phys. Rev. B 88, 075420 (2013).

45. Chiang, T.-C. Photoemission studies of quantum well states in thin lms. Surf. Sci. Rep. 39, 181235 (2000).

46. Zhang, Y. et al. Direct observation of the transition from indirect to direct bandgap in atomically thin epitaxial MoSe2. Nat. Nanotechnol. 9, 111115 (2014).

NATURE COMMUNICATIONS | 5:4673 | DOI: 10.1038/ncomms5673 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms5673

47. Yu, Y. et al. Controlled scalable synthesis of uniform, high-quality monolayer and few-layer MoS2 lms. Sci. Rep. 3, 1866 (2013).

48. Hartwigsen, C., Goedecker, S. & Hutter, J. Relativistic separable dual-space Gaussian pseudopotentials from H to Rn. Phys. Rev. B 58, 36413662 (1998).

49. Gonze, X. et al. First-principles computation of material properties: the ABINIT software project. Comput. Mater. Sci. 25, 478492 (2002).

50. Gonze, X., Rignanese, G.-M. & Caracas, R. First-principle studies of the lattice dynamics of crystals, and related properties. Z. Kristallogr. 220, 458472 (2005).

Acknowledgements

Work at Princeton University is supported by the US National Science Foundation

Grant, NSF-DMR-1006492. M.Z.H. acknowledges visiting-scientist support from

Lawrence Berkeley National Laboratory and additional partial support from the A.P.

Sloan Foundation and NSF-DMR-0819860. The photoemission measurements using

synchrotron X-ray facilities are supported by the Basic Energy Sciences of the US

Department of Energy. Theoretical computations are supported by the US Department

of Energy (supported by DE-FG02-05ER46200 and Feynman cluster at Princeton

University, Department of Physics). Sample growth and characterization are partially

supported by DMR-0819860, DMR-1006492 and by the National Science Council (NSC),

Taiwan, under project number NSC100-2119-M- 002-021.

Author contributions

N.A. and M.Z.H. conceived and designed the experiments. N.A. performed the experi

ments with assistance from S.-Y.X., M.N., C.L., I.B. and J.D.D.; G.B. performed DFT

calculations and spectral simulations; N.A., G.B. and M.Z.H. performed data analysis,

gure planning and draft preparation; D.-X.Q., R.S. and F.-C.C. provided samples and

performed sample characterization; M.Z.H. was responsible for the overall direction,

planning and integration among different research units.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/

http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsand

Web End =http://npg.nature.com/

http://npg.nature.com/reprintsand

Web End =reprintsand permissions/

How to cite this article: Alidoust, N. et al. Observation of monolayer valence band

spin-orbit effect and induced quantum well states in MoX2. Nat. Commun. 5:4673

doi: 10.1038/ncomms5673 (2014).

8 NATURE COMMUNICATIONS | 5:4673 | DOI: 10.1038/ncomms5673 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

Word count: 6485

Show less

Abstract

Translate

Transition metal dichalcogenides transition metal dichalcogenides have attracted much attention recently due to their potential applications in spintronics and photonics because of the indirect to direct band gap transition and the emergence of the spin-valley coupling phenomenon upon moving from the bulk to monolayer limit. Here, we report high-resolution angle-resolved photoemission spectroscopy on MoSe₂ single crystals and monolayer films of MoS₂ grown on highly ordered pyrolytic graphite substrate. Our experimental results resolve the Fermi surface trigonal warping of bulk MoSe₂ , and provide evidence for the critically important spin-orbit split valence bands of monolayer MoS₂ . Moreover, we systematically image the formation of quantum well states on the surfaces of these materials, and present a theoretical model to account for these experimental observations. Our findings provide important insights into future applications of transition metal dichalcogenides in nanoelectronics, spintronics and photonics devices as they critically depend on the spin-orbit physics of these materials.

Details

Title

Observation of monolayer valence band spin-orbit effect and induced quantum well states in MoX2

Author

Alidoust, Nasser; Bian, Guang; Xu, Su-yang; Sankar, Raman; Neupane, Madhab; Liu, Chang; Belopolski, Ilya; Qu, Dong-xia; Denlinger, Jonathan D; Chou, Fang-cheng; Hasan, M Zahid

Pages

4673

Publication year

2014

Publication date

Aug 2014

Publisher

Nature Publishing Group

e-ISSN

20411723

Source type

Scholarly Journal

Language of publication

English

DOI

https://doi.org/10.1038/ncomms5673

ProQuest document ID

1555275144

Observation of monolayer valence band spin-orbit effect and induced quantum well states in MoX2

Jump to:

Full text

Abstract

Details

Suggested sources