Printer Friendly

Dopant Introduced Valley Polarization, Spin, and Valley Hall Conductivity in Doped Monolayer Mo[S.sub.2].

1. Introduction

The charge and spin degrees of freedom (DOF) of electrons are at heart of modern electronics as they form the basis of a wide range of applications, such as transistors, photodetectors, and magnetic memory devices. In addition to this, electrons in two-dimensional (2D) crystals that have a honeycomb lattice structure possess an extra valley DOF in addition to charge and spin [1]. Valley is a quantum number defined in an electronic system whose band structure contains energetically degenerate but nonequivalent structures due to a certain crystal structure. Among various material candidates for valleytronics, spatial inversion symmetry broken two-dimensional (2D) honeycomb lattice systems such as graphene and monolayer Mo[S.sub.2] are predicted to be the most useful. These systems have two valleys called K and K'. The use of valley indexes for a potential information carrier was first suggested in the studies of conventional semiconductors such as AlAs and Si [2].On the other hand, the presence of a valley-dependent orbital magnetic moment in graphene with gap suggests that currents flow perpendicular to applied electric field even in the absence of a magnetic field, named the 'valley Hall effect.' The valley Hall effect was first observed in graphene with substrate. However, due to the small value of the energy band gap this effect has not been yet observed experimentally. Recently, photogenerated valley Hall conductivity was observed in monolayer Mo[S.sub.2] [3]. In contrast, to graphene in monolayer Mo[S.sub.2] the inversion symmetry is explicitly broken [2], which can give rise to the valley Hall effect where carriers in different valleys flow to opposite transverse edges when an in-plane electric field is applied [4, 5]. In addition to this, Mo[S.sub.2] has a strong spin-orbit coupling (SOC) originating from the d orbitals of the heavy metal atoms(Mo)[4] and can be an interesting platform to explore spin physics and spintronics applications which are absent in graphene due to its vanishing SOC [6, 7]. However, as far as both valleys are concerned due to time reversal symmetry analogy to paramagnetic and paraelectric materials, pure monolayers Mo[S.sub.2] are paravlley material. Therefore, the major challenge in valleytronics is to break the degeneracy between the two prominent [K.sup.+] and [K.sup.-] valleys so as to achieve the valley polarization. Recent experimental investigation demonstrated optical means of breaking valley degeneracy [8-11]. In addition to this, other approaches like applied electric field [12], by electron-electron interaction [13], by means of external magnetic field [14-17], are also recently proposed mechanism to break valley degeneracy.

In this paper, we theoretically study dynamics of valley polarization, spin, and valley Hall conductivity in the vicinity of [K.sup.+] and [K.sup.-] valleys of 2D hexagonal Brillouin zone (BZ) of monolayer Mo[S.sub.2] considering spin-orbital coupling and dopant introduced exchange field. In second section using model Hamiltonian together with Berry phase approach we obtained explicit expression of Berry curvature, spin and valley split form of Hall conductivity and total valley and spin Hall conductivity. In section three, we present the numerical results showing effect of the spin-orbit coupling and dopant introduced exchange field on energy band structure of two prominent valleys ([K.sup.+] and [K.sup.-]) and total valley and spin Hall conductivity. Finally the conclusion is presented in last section.

2. Theoretical Model

We consider an effective three-band tight binding model in the vicinity of [K.sup.+] and [K.sup.-] valleys of hexagonal Mo[S.sub.2] Brillouin zone including spin-orbital coupling originating from transition element Mo 4d states and the dopant introduced Hamiltonian. In symmetry based three-band tight binding model, the conduction bands are dominated by 4d[z.sup.2] and the valence bands are dominated by [4d.sub.xy] and 4d[x.sup.2] - [y.sup.2]. Since a three-band TB model is enough to capture the band-edge properties in the vicinity of [K.sup.[+ or -]] valleys [15-17] including energy dispersions, Berry curvature, and valley Hall and spin Hall conductivity. Therefore, it is convenient to use this model to study desired properties in the vicinity of those valleys.

The Hamiltonian of the system will have the form

[mathematical expression not reproducible]. (1)

Following [2, 18] the low energy effective Hamiltonian in K.P theory of the system reads

[[??].sub.k.p] = at (v[[tau].sub.x][k.sub.x] + [[tau].sub.y][k.sub.y]) + [[DELTA].sub.m][[tau].sub.3], (2)

and spin-orbital term Hamiltonian can be written as

[[??]] = v[[lambda]] (1 - [[tau].sub.3])[s.sub.z], (3)

whereas exchange energy resulting from dopants is assumed to be

[[??].sup.dop.sub.ex] = -v[s.sub.z][h.sub.ex][[sigma].sub.z]. (4)

Here, [[??].sup.dop.sub.exc] is magnetic dopant (Mn, V) induced magnetic exchange energy and is the energy splitting resulting from exchange interaction between dopants, t and a are first-order valence (conduction) band effective hopping and lattice constant for Mo[S.sub.2], respectively, and [[DELTA].sub.m] and [[lambda]] are mass correction for Dirac Hamiltonian and spin-orbit coupling parameter, respectively. v = [+ or -]1 designates the value index for [K.sup.[+ or -]], [[tau].sub.x,y,z] and [[sigma].sub.x,y,z] are Pauli matrices in lattice space and spin space, respectively. Using the value of Dirac matrices the matrix form of Hamiltonian reads

[mathematical expression not reproducible], (5)

where we have introduced short notation [k.sup.v.sub.-] = [k.sub.x]v - i[k.sub.y] and its complex conjugates [k.sup.v.sub.+] = [k.sub.x]v + i[k.sub.y]. The eigenfunction of tight binding Hamiltonian can be expressed in terms of product of periodic part and plane wave as

[[psi].sub.p](k) = [1/[square root of A] exp (ik.r) [u.sup.v.sub.p] (k) (6)

where A is area of the two-dimensional system and (1[square root of A]) exp(ik.r) is the plane wave solution corresponding to kinetic (hopping term) of Hamiltonian in (5) and [u.sup.v.sub.s](k) is periodic part of Bloch state which can be written as a two-component way considering spin-up and spin-down components. After solving characteristic equation, the energy eigen value and energy eigenvector, respectively, is obtained as

[mathematical expression not reproducible]. (7)

[mathematical expression not reproducible], (8)


[mathematical expression not reproducible]. (9)

[mathematical expression not reproducible]. (10)

On the other hand, from our previous work [19], the general expression for z component Berry curvature which is analogical with real space magnetic field is given by expression

[mathematical expression not reproducible]. (11)

Now, fascinatingly, this expression of Berry curvature for particular system is under consideration together with using (7) and (8). Making transformation, [mathematical expression not reproducible], and [[nabla].sub.y] [right arrow] [partial derivative]/[partial derivative]ky using those variable changes in (11), we shall rewrite Berry curvature of conduction band minimum (CBM) and valance band maximum (VBM) in the vicinity of [K.sup.-] and [K.sup.+] valleys as

[mathematical expression not reproducible]. (12)

where x and y component of velocity in (12) can be obtained using the relation between wave vector dependent Hamiltonian in (5) and momentum operator as

[mathematical expression not reproducible], (13)

[mathematical expression not reproducible]. (14)

Substituting (13) and (14) into (12) together with simple manipulation yields the Berry curvature for upper and lower band ([+ or -]) as

[mathematical expression not reproducible], (15)

In (15), [e.sub.z] is unit vector along Z axis. To obtain the expression [e.sub.z] explicitly we introduce a three-dimensional unit vector along x, y, and z in a two-dimensional plane in k space in terms of eigenvalues as

[mathematical expression not reproducible]. (16)

Using (16) into (15), we find

[mathematical expression not reproducible], (17)

where p=1 for conduction band curvature and p=-1 for valence band curvature. On the other hand, using linear response theory the quantum Kubbo formula for Hall conductivity in current-Current correlated system is found to be

[mathematical expression not reproducible]. (18)

In dilute limit or in absence of disorder effect, after some mathematical algebra, z components of conductivity in (18) are simplified to

[mathematical expression not reproducible]. (19)

On account of (12), the right side of (19) inside the bracket gives (-[[OMEGA].sup.p](z)/2) which reproduces conductivity in terms of Berry curvature as

[mathematical expression not reproducible] (20)

Following [20, 21] z component of total spin and valley Hall conductivity in the vicinity of both Valleys, namely, [K.sup.+] and [K.sup.-], respectively, can be obtained using

[mathematical expression not reproducible]. (21)

[mathematical expression not reproducible]. (22)

where [+ or -] designates short notation for [K.sup.+] and [K.sup.-] valleys, respectively.

First we evaluate both spin and valley split form of anomalous Hall conductivity [[sigma].sup.[+ or -].sub.[up arrow][down arrow]] using (20) but considering two limiting cases.

Case I. At absolute zero and at low impurity concentration or in pure monolayer Mo[S.sub.2], in this scenario all electron states below the Fermi level will contribute (i.e., it is a Fermi sea property). We use the approximation f([[zeta].sub.k+]) = 0 and f([[zeta].sub.k-]) = 1. Here it should be known that the word 'anomalous' reflects Hall conductivity in absence of external magnetic field. Therefore, in this limit after transforming summation into integration, (20) becomes

[mathematical expression not reproducible]. (23)

After solving (23) the anomalous Hall conductivity in Valley [K.sup.+] in occupied valance band, for spin-up and spin-down electrons (holes)

[right arrow] [[sigma].sup.+.sub.[up arrow]](z) = [[sigma].sup.+.sub.[down arrow]] (z) = -[[q.sup.2]/h]. (24)

Similar to anomalous Hall conductivity in the vicinity of [K.sup.-] valleys in occupied valance band, for spin-up and spin-down electrons (holes),

[[sigma].sup.-.sub.[up arrow]](z) = [[sigma].sup.-.sub.[down arrow]](z) = [[q.sup.2]/h]. (25)

From these results we can understand that, in pure limit, anomalous Hall conductivity in a single spin state (spin-up or spin-down) and in single valley ([+ or -]) is quantized (depends only on fundamental unit of conductivity). On the other hand, as seen from (24) and (25) the signs of calculated anomalous Hall conductivity for [K.sup.+] and [K.sup.-] are opposite to each other which indicates that those valleys are related to time reversal symmetry.

The total spin Hall conductivity and valley Hall conductivity in the vicinity of both valleys ([K.sup.+] and [K.sup.-]) is obtained by substituting (24) and (25) into (21) and (22)

[[sigma].sup.spin] (z) = 0. (26)

[[sigma].sup.valley] (z) = -4[[q.sup.2]/h]. (27)

Case II. At finite doping concentration, in this limit both spin hall and valley Hall conductivity should be calculated by taking into consideration both occupied bands from -[infinity] to 0 plus the conduction band until ([[member of].sub.p]) or ([k.sub.F]. Hence in this limit (20) is rewritten as

[mathematical expression not reproducible]. (28)

After transforming summation into integration

[mathematical expression not reproducible] (29)

Substitute (17) into (29) and evaluating for [[sigma].sup.+.sub.[up arrow]](z), [[sigma].sup.+.sub.[down arrow]](z), [[sigma].sup.-.sub.[up arrow]](z), and[[sigma].sup.+.sub.[down arrow]](z) by varying the values of v and [s.sub.z] in (20). To evaluate [[sigma].sup.+.sub.[up arrow]](z) from (20), we use v = 1, [s.sub.z] = 1, p =- 1 and after series of steps

[mathematical expression not reproducible]. (30)

Similarly substitution of v = 1, [s.sub.z] = -1, results in the expression of spin-down electrons in the vicinity of [K.sup.+] valleys as

[mathematical expression not reproducible]. (31)

Following the same procedure, [[sigma].sup.-.sub.[up arrow]](z) and [[sigma].sup.+.sub.[down arrow]](z) are found to be

[mathematical expression not reproducible]. (32)

[mathematical expression not reproducible]. (33)

Total Spin Hall conductivity at finite doping concentration at zero temperature can be obtained by substituting ((30)-(33)) into (21) as

[mathematical expression not reproducible]. (34)

Finally, the total valley Hall conductivity can be calculated by making use of those ((30)-(33)) into (22)

[mathematical expression not reproducible]. (35)

Equation (34) reveals that unlike the case of pure, limit the total spin Hall conductivity is finite due to spin polarization of [K.sup.-] and [K.sup.+] valleys associated with dopant introduced magnetic exchange field ([h.sub.ex]).

3. Result and Discussion

The energy plotted as a function of dimensionless wave number (ka/[pi]) is displayed in Figure 1. In Figures 1(a) and 1(b), we find that for zero spin-orbital coupling (lo=0) and zero dopant introduced exchange field ([h.sub.ex]=0), both bands ([K.sup.+] and [K.sup.-]) valleys are degenerate. On the other hand, for nonzero spin-orbit coupling and zero exchange field as shown in Figures 1(c) and 1(d), conduction band [K.sup.+] point eigenvalues are identical to those of the [K.sup.-] point if spin-up and down are exchanged but the degeneracy in valence band hole state in [K.sup.+] and [K.sup.-] valleys are lifted due to spin- orbital coupling resulting from transition metal (Mo) orbitals. However, for finite spin-orbit coupling (lo) and finite dopant introduced exchange field degeneracy in both valence band holes and conduction band electrons in valley [K.sup.+] and [K.sup.-] is lifted Figures 1(e) and 1(f). Moreover, to understand the role of dopant introduced magnetic exchange field ([h.sub.ex]) on momentum space Berry curvature, we plot the Berry curvature as function of wave number (ka/[pi]) in the vicinity of [K.sup.+] and [K.sup.-] valleys for spin-up holes in valence band Figure 2, enhancement of pick of curvature with rise in dopant introduced exchange field (0 to 0.16 eV) for state in the vicinity [K.sup.+] valleys are observed Figure 2.

On the other hand, Berry curvature is related to both spin Hall and valley Hall through relation (20). Therefore, the spin and valley Hall conductivity is expected to increase with dopant introduced exchange field for state in the vicinity of [K.sup.+] valleys. In contrast with Berry curvature of state in the vicinity of [K.sup.+] valley, the Berry curvature of state in the vicinity of [K.sup.-] valley suppresses with rise in dopant introduced exchange field (0 to 0.16 eV). Therefore, the degeneracy of curvature of [K.sup.+] and [K.sup.-] is lifted under dopant introduced magnetic exchange field or these valleys become spin polarized (unequal number of spin-up states are accumulated in [K.sup.+] and [K.sup.-] curvature). In addition to this, at absolute zero (T=0K) and insulating regime where the chemical potential is in gap, (24) and (25) show that the universal value of conductivity [q.sup.2]/h is obtained for single spin and single valleys. Hence, we can say that Hall conductivity is quantized in this limit and the signs of conductivity in [K.sup.+] and [K.sup.-] are opposite indicating that both of valleys are related to time reversal symmetry. On the other hand, from (26) and (27) and the total spin Hall conductivity vanishes identically. Whereas total valley Hall conductivity is integral multiple of quantum conductivity (4([q.sup.2]/h)) the term 4 appears from twofold valley degrees of freedom (DOF) corresponding to [K.sup.-] and [K.sup.+] valleys of the BZ and twofold to spin degree of freedom of those carriers. We have also calculated both spin and Valley Hall conductivities at finite doping concentration. The calculated result reveals that in contrast to conductivity in insulating regime Hall conductivity is not quantized but depends on parameters ([[DELTA].sub.m], [[lambda]], [h.sub.ex], [k.sub.F], a, and t); see ((32)-(35)). The calculated total spin Hall conductivity and total valley Hall conductivity in (34) and (35) can be written in more compact form as

[mathematical expression not reproducible]. (36)

[mathematical expression not reproducible]. (37)

From (36) and (37), total spin Hall conductivity and valley Hall conductivity depend on parameters: [[DELTA].sub.m], [[lambda]], [h.sub.ex], [k.sub.F], a, and t. For particular limits [h.sub.ex] = [[lambda]] = 0 and [h.sub.ex] = 2[[lambda]] total spin Hall conductivity in (36) vanishes whereas the total valley Hall conductivity is simplified to -(4[q.sup.2]/h) [(([[DELTA].sup.2.sub.m] + [a.sup.2] [t.sup.2] [k.sup.2.sub.F])/at[k.sub.F]).sup.1/2], which indicates that the interplay between spin-orbit coupling and dopant introduced magnetic exchange field is crucial for existence of spin and valley Hall conductivity.

4. Conclusion

In conclusion, based on a three-band tight binding model we have studied valley polarization, spin, and valley Hall conductivity in the presence of dopants introduced magnetic exchange field. The result reveals that in pure limit Hall conductivity in single valley is quantized. In finite doping concentration the existence of dopant introduced exchange field breaks the time inversion symmetry and decouples the energetically degenerated valley into nondegenerate and elucidating the occurrence of valley polarization. We have shown that the Berry curvature of occupied state resulting from transverse motion unequal number of carriers in both valleys of 2D hexagonal lattice is responsible for existence of valley and spin Hall conductivity. Moreover, we have found that, at finite dopant concentration, the strength of the spin-orbit coupling together with the exchange energy determines the valley polarization, which in turn controls total valley and spin Hall conductivity in doped monolayer Mo[S.sub.2] system. The spin Hall and valley Hall conductivity are dissipationless in the absence of any external magnetic field. This effect could enable a new generation of low power electronic devices.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper. Submitting authors are responsible for coauthors declaring their interest.


This work was supported by the School of Graduate Studies of Addis Ababa University, Addis Ababa, Ethiopia, and Arbaminch University, Arbaminch, Ethiopia.


[1] A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, "The electronic properties of graphene," Reviews of Modern Physics, vol. 81, no. 1, pp. 109-162, 2009.

[2] D. Xiao, G. Liu, W. Feng, X. Xu, and W. Yao," Coupled Spin and Valley Physics in Monolayers of," Physical Review Letters, vol. 108, no. 19, 2012.

[3] K. F. Mak, K. L. McGill, J. Park, and P. L. McEuen, "The valley hall effect in MoS2 transistors," Science, vol. 344, no. 6191, pp. 1489-1492, 2014.

[4] W. Yao, A. H. MacDonald, and Q. Niu, "Optical Control of Topological Quantum Transport in Semiconductors," Physical Review Letters, vol. 99, no. 4, 2007.

[5] W. Yao, D. Xiao, and Q. Niu, "Valley-dependent optoelectronics from inversion symmetry breaking," Physical Review B: Condensed Matter and Materials Physics, vol. 77, no. 23, 2008.

[6] H. Min, J. E. Hill, N. A. Sinitsyn, B. R. Sahu, L. Kleinman, and A. H. MacDonald, "Intrinsic and Rashba spin-orbit interactions in graphene sheets," Physical Review B: Condensed Matter and Materials Physics, vol. 74, no. 16, 2006.

[7] Y. Yao, F. Ye, X. Qi, S. Zhang, and Z. Fang, "Spin-orbit gap of graphene: First-principles calculations," Physical Review B: Condensed Matter and Materials Physics, vol. 75, no. 4, 2007.

[8] A. Umer, S. Naveed, N. Ramzan, and M. S. Rafique, "Selection of a suitable method for the synthesis of copper nanoparticles," Nanotech, vol. 7, no. 5, Article ID 1230005, 2012.

[9] H. Zeng, J. Dai, W. Yao, D. Xiao, andX. Cui, "Valley polarization in MoS2 monolayers by optical pumping," Nat Nanotechnol, vol. 7, no. 8, pp. 490-493, 2012.

[10] T. Cao, G. Wang, W. Han, and H. Ye, "Valley-selective circular dichroism of monolayer molybdenum disulphide," Nature Communications, vol. 3, article no. 887, 2012.

[11] K. F. Mak, K. He, J. Shan, and T. F. Heinz, "Control of valley polarization in monolayer Mo[S.sub.2] by optical helicity," Nature Nanotechnology, 2012.

[12] G. Kioseoglou and A. T. Hanbicki, "Anomalous temperature dependent spin-valley polarization in monolayer WS2," Scientific Reports, vol. 6, article no. 18885, 2016.

[13] J. Shi, G. Vignale, D. Xiao, and Q. Niu, "Quantum Theory of Orbital Magnetization and Its Generalization to Interacting Systems," Physical Review Letters, vol. 99, no. 19, 2007.

[14] A. Srivastava, M. Sidler, A. V. Allain, and D. S. Lembke, "Valley Zeeman effect in elementary optical excitations of monolayer W[Se.sub.2]," Nature Physics, pp. 141-147, 2015.

[15] A. V. Chumak, V. I. Vasyuchka, A. A. Serga, and B. Hillebrands, "Magnetic control of valley pseudospin in monolayer W[Se.sub.2]," Nature Physics, vol. 11, no. 6, pp. 148-152, 2015.

[16] D. MacNeill, C. Heikes, K. F. Mak, and Z. Anderson, "Breaking of Valley Degeneracy by Magnetic Field in Monolayer Mo[Se.sub.2]," Physical Review Letters, vol. 114, no. 19, Article ID 037401, 2015.

[17] Y. Li, J. Ludwig, T. Low, A. Chernikov, and X. Cui, "Valley Splitting and Polarization by the Zeeman Effect in Monolayer Mo[Se.sub.2]," Physical Review Letters, vol. 113, Article ID 266804, 2014.

[18] A. Kormanyos, G. Burkard, M. Gmitra, J. Fabian, and V. Zolyomi, "k.p theory for two-dimensional transition metal dichalcogenide semiconductors," 2D Materials, vol. 2, Article ID 022001, 2015,

[19] S. Mekonnen and P. Singh, "Berry approach to intrinsic anomalous Hall conductivity in dilute magnetic semiconductors (Ga1xMnxAs)," World Journal of Condensed Matter Physics, vol. 5, pp. 179-186, 2015.

[20] Y. Yao, L. Kleinman, A. H. MacDonald, J. Sinova, and T. Jungwirth, "First Principles Calculation of Anomalous Hall Conductivity in Ferromagnetic bcc Fe," Physics Letters B, vol. 92, 2004.

[21] G. Liu, D. Xiao, Y. Yao, X. Xu, and W. Yao, "Electronic structures and theoretical modelling of two-dimensional group-VIB transition metal dichalcogenides," Chemical Society Reviews, vol. 44, no. 9, pp. 2643-2663, 2015.

Sintayehu Mekonnen (iD) and Pooran Singh

Addis Ababa University, Department of Physics, Addis Ababa, Ethiopia

Correspondence should be addressed to Sintayehu Mekonnen;

Received 7 February 2018; Revised 27 April 2018; Accepted 2 May 2018; Published 1 August 2018

Academic Editor: Bakhtiar U. Haq

Caption: Figure 1: Electronic band structure near valley [K.sup.[+ or -]]: (a) and (b) for zero both spin-orbit coupling and exchange field (lo=hex=0), (c) and (d) for spin-orbit coupling(lo=0.08 eV,) and zero exchange field, and (e) and (f) for spin-orbit coupling (lo=0.08 eV) and exchange field (hex=0.2 eV). Red: spin-up electron (hole) states, blue: spi- down electron (hole) states, right panel for [K.sup.+], and left panel for [K.sup.-] valley.

Caption: Figure 2: Berry curvature as function of (ka/[pi]) for different values of exchange fields ([h.sub.ex]), right panel for [K.sup.+], and left panel for [K.sup.-] valley.
COPYRIGHT 2018 Hindawi Limited
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2018 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Title Annotation:Research Article
Author:Mekonnen, Sintayehu; Singh, Pooran
Publication:Advances in Condensed Matter Physics
Article Type:Report
Geographic Code:1USA
Date:Jan 1, 2018
Previous Article:Self-Assembled and Artificial Surfaces/Interfaces: From Soft Matter to Metamaterials.
Next Article:High Sensitivity Refractive Index Sensor by D-Shaped Fibers and Titanium Dioxide Nanofilm.

Terms of use | Privacy policy | Copyright © 2021 Farlex, Inc. | Feedback | For webmasters |