Open access peer-reviewed chapter

Valley Polarized Single Photon Source Based on Transition Metal Dichalcogenides Quantum Dots

By Fanyao Qu, Alexandre Cavalheiro Dias, Antonio Luciano de Almeida Fonseca, Marco Cezar Barbosa Fernandes and Xiangmu Kong

Submitted: November 22nd 2016Reviewed: July 4th 2017Published: October 25th 2017

DOI: 10.5772/intechopen.70300

Downloaded: 1503


Photonic quantum computer, quantum communication, quantum metrology, and optical quantum information processing require a development of efficient solid‐state single photon sources. However, it still remains a challenge. We report theoretical framework and experimental development on a novel kind of valley‐polarized single‐photon emitter (SPE) based on two‐dimensional transition metal dichalcogenides (TMDCs) quantum dots. In order to reveal the principle of the SPE, we make a brief review on the electronic structure of the TMDCs and excitonic behavior in photoluminescence (PL) and in magneto‐PL of these materials. We also discuss coupled spin and valley physics, valley‐polarized optical absorption, and magneto‐optical absorption in TMDC quantum dots. We demonstrate that the valley‐polarization is robust against dot size and magnetic field, but optical transition energies show sizable size‐effect. Three versatile models, including density functional theory, tight‐binding and effective k⋅p method, have been adopted in our calculations and the corresponding results have been presented.


  • single-photon source
  • quantum dots
  • transition metal dichalcogenides

1. Introduction

Traditional semiconductors have been used for decades for making all sorts of devices like diodes, transistors, light emitting diodes, and lasers [1]. Due to the advances of technology in fabrication, it is possible not only to make ever pure semiconductor crystals, but also to study heterostructures, in which carriers (electrons or holes) are confined in thin sheets, narrow lines, or even a point [1, 2]. Quantum dots (QDs) are zero‐dimensional objects where all the three spatial dimensions are quantized with sizes smaller than some specific characteristic lengths, e.g., the exciton Bohr radius [1, 2]. Because of confinement, electrons in the QDs occupy discrete energy levels, in a similar way as they do in atoms [2, 3]. For these reasons, QDs are also referred to as artificial atoms [1, 2]. In spite of some similarities between the QD and the real atom, the former demonstrates several special characteristics. For instance, its size can vary from a few to hundreds of nanometers, and it can trap from a very small number of electrons (Ne<10) to 50100electrons or more [1]. In addition, the shape of the QD that can be tuned at will determine its spatial symmetry. In turn, the change in spatial symmetry modifies physical properties of the system. As known, the three‐dimensional spherically symmetric QDs possess degenerate electron shells, 1s, 2s, 2p, 3s, 3p,… [3]. When the number of electrons is equal to 2, 10, 18, 36,…, the electron shells are completely filled [1], yielding a particularly stable configuration. In contrast, a two‐dimensional cylindrically symmetric parabolic potential leads to formation of a two‐dimensional shell structure with the magic numbers 2, 6, 12, 20,…[2, 4, 5]. Hence, the lower degree of symmetry in two‐dimensional QDs leads to a smaller magic number sequence. Since the shape and size of the QDs can be precisely manufactured, the energy structure of the carriers in the QDs as well as their optical, transport, magnetic, and thermal properties can be engineered in a large scale [19].

Various techniques have been developed to produce the QDs such as etching, regrowth from quantum well structures, beam epitaxy, lithography, holograph patterning, chemical synthesis, etc [1, 2]. Consequently, many kinds of QDs emerge. According to the electrical property of their parent material, they can be classified into metal, semiconductor, or super‐conducting dots. From geometry point of view, the QDs form two groups: two‐dimensional [2, 4, 5, 8] or three‐dimensional (3D) [1] dots. The former can be further divided into conventional 2D semiconductor QDs, such as self‐assembled‐ and gated‐QDs based on traditional semiconductor quantum wells [48, 10], and the novel QDs made from two‐dimensional‐layered materials (2DLMs) [1115].

Atomically thin 2DLMs have revolutionized nanoscale materials science [16]. The interatomic interaction within layers is covalent in nature, while the layers are held together by weak van der Waals (vdW) forces. The family of 2D materials, which started with graphene [16], has expanded rapidly over the past few years and now includes insulators, semiconductors, semimetals, metals, and superconductors [1719]. The most well studied 2D systems beyond graphene, are the silicene, germanene, stanine, and borophene, organic‐inorganic hybrid perovskites, insulator hexagonal boron‐nitride [17, 18], the anisotropic semiconductor phosphorene, transition metal‐carbides, ‐nitrides, ‐oxides, and ‐halides, as well as the transition metal dichalcogenides (TMDCs) [2028]. Compared with traditional semiconducting materials, the 2DLMs take advantage of inherent flexibility and an atomically‐thin geometry. Moreover, because of their free dangling bonds at interfaces [25, 29, 30], two‐dimensional‐layered materials can easily be integrated with various substrates [17]. They can also be fabricated in complex‐sandwiched structures or even suspended to avoid the influence of the substrate [31]. The monolayer TMDCs with infinite geometry exhibit strong carrier confinement in one dimension but preserve the bulk‐like dispersion in the 2D plane. In contrast, electrons in a TMDC QD are restricted in three dimensions, which present size tunable electronic and optical properties in addition to the remarkable characteristics related to spin‐valley degree of freedom inherited from its 2D bulk materials. Very recently, graphene QDs (GQDs) have attracted intensive research interest due to their high transparency and high surface area. Many remarkable applications ranging from energy conversion to display to biomedicine are prospected [11]. Nevertheless, from quantum nano‐devices point of view, the TMDCs have advantages over graphene. For instance, the semiconducting TMDCs have a band gap large enough to form a QD using the electric field, as shown in Figure 1, unlike etched GQDs made on semi‐metallic graphene.

Figure 1.

(a) Scanning electron microscope image of theWS2quantum dot studied. TheWS2flake is highlighted by the white dotted line, and the four top gates are labeled as MG, LB, PG, RB. The scale bar represents 5μm. (b) Three‐dimensional schematic view of the device [20]. (Copyright 2015 by the Royal Society of Chemistry. Reprinted with permission).

The applications of quantum dots are still mostly restricted to research laboratories, but they are remarkable due to the fact that QDs provide access to the quantum mechanical degrees of freedom of few carriers. Single electron transistors [16], the manipulation of one [47] or two [13] electron spins, manipulation of a single spin in a single magnetic ion‐doped QD [47] are only some examples. Optically active quantum dots can also be used in both quantum communication and quantum computation [47, 1215]. The emerging field of quantum information technology, as unconditionally quantum cryptography, quantum‐photonic communication and computation, needs the development of individual photon sources [1215, 32, 33]. Recently, individual photon emitters based on defects in TMDC monolayers with different sample types (WSe2and MoSe2) have been reported, but only operate at cryogenic temperatures [12, 13, 3235]. In addition, it should be kept in mind that the presence of defects is not always beneficial for the PL signal. For instance, defect‐mediated nonradiative recombination might result in an internal quantum yield droop in the defective TMDCs [2628]. In this context, single quantum emitter based on QD is desirable [14, 15].

In this chapter, we show the optical and magneto‐optical properties of the TMDC QD's. We choose MoS2, which has been widely studied in the literature as our example. Three versatile models including density functional theory, tight‐binding, and effective kpapproach have been adopted in our calculations and the corresponding results have been presented. We show that the valley‐polarization is robust against dot size and magnetic field, but the optical transition energies show sizable size effect. Based on the computed optical absorption spectra, a novel kind of valley‐polarized single‐photon source based on TMDC quantum dot is proposed [15].


2. Physical properties of transition metal dichalcogenides

2.1. Electronic band structure of transition metal dichalcogenides

Layered TMDCs have the generic formula MX2, where M stands for a metal and X represents a chalcogen. The monolayer crystal structure of the MoS2, which is one of the most studied TMDCs in the literature is shown in Figure 2. Notice that the monolayer MoS2has a trigonal prism crystal structure. An inversion asymmetry results in a large direct gap semiconductor with the gap lying at the two inequivalent K‐points of the hexagonal Brillouin zone.

Figure 2.

(a) Schematic diagram of crystal structure ofMoS2. (b) Coordination environment of Mo (blue sphere) in the structure (the left panel); the middle and the right panels correspond to side‐ and top‐ views of the monolayerMoS2lattice. Sulfur is shown as golden spheres.

The major orbital contribution at the edge of the conduction band (CB) is from d3z2r2orbital of the metal, plus minor contributions from px, and pyorbitals of chalcogens. On the other hand, at the edge of the valence band (VB) at the K‐point, the most important orbital contribution is due to a combination of dxyand dx2y2of the metal, which hybridize to pxand pyorbitals of the chalcogen atoms [24]. In addition, there is a strong spin‐orbit interaction (SOI), especially in the valence band. The spatial inversion asymmetry along with strong SOI leads to a spin‐valley coupling, which results in the spin of electron being locked to the valley [20, 21, 36], as demonstrated in Figure 3(a). Furthermore, the band structure in the K‐valley presents time reversal symmetry with that in K‐valley. And, it changes dramatically as the number of the layers of the TMDC increases. When the thickness increases, the band gap in MoS2and other group VI TMDCs decreases, and more importantly, the material becomes an indirect gap semiconductor [37], as shown in Figure 3(b).

Figure 3.

(a) Band structure of monolayerMoS2. The solid curves were obtained using the QUANTUM ESPRESSO package with fully relativistic pseudopotentials under the Perdew‐Burke‐Ernzerhof generalized‐gradient approximation, and a16×16×1kgrid. The dashed curves were calculated from the tight‐binding model, with cyan (red) representing states that are even (odd) under mirror operation with respect to the Mo plane.v1,2andc1,2label the bands close to the valence and conduction band edges near theK‐ andK'‐points. The inset shows the hexagonal Brillouin zone (pink) associated with the triangular Bravais lattice ofMoS2and an alternate rhombohedral primitive zone (black), and labels the principle high‐symmetry points in reciprocal space [36] (Copyright 2015 by American Physical Society, Reprinted with permission). (b) Calculated band structures of (b1) bulkMoS2, (b2) quadrilayerMoS2, (b3) bilayerMoS2, and (b4) monolayerMoS2. The solid arrows indicate the lowest energy transitions. BulkMoS2is characterized by an indirect bandgap. The direct excitonic transition occurs atKpoint with a higher transition energy than that of indirect one [30]. (Copyright 2010 by the American Chemical Society. Reprinted with permission).

In order to get insight into the physical origins of the band gap variations with the number of layers, Figure 4 shows evolutions of the band gaps (a) and band edges (b) of MoS2as a function of the number of layers. Notice that with increasing layer thickness, the indirect band gap (ΓK, ΓΛ) becomes smaller, while the direct excitonic transition (KK)only slightly changes, as illustrated in Figure 4(a). Note also that the monolayer (n=1)MoS2is a direct band gap semiconductor, but it becomes an indirect band gap semiconductor when the number of layers is larger than one. This phase transition is also clearly demonstrated in Figure 4(b). In addition, as the number of layers increases, both the VB and CB edges at K‐valley exhibit only slight variation while the degeneracy at Γis lifted and a splitting of the bands occurs, which pushes the VB maximum to higher energy. The variation of the band gap is largely driven by the variation of the VB at Γpoint. Going from a monolayer to a bilayer significantly raises the VB at Γ, resulting in a transition from the direct KKgap to an indirect ΓKgap. This dramatic change of electronic structure in monolayer MoS2results in the jump in monolayer photoluminescence efficiency.

Figure 4.

(a) Evolution of the band gaps as a function of the number of layers (n). The black circles(KK), red squares(ΓK), and green diamonds(ΓΛ)indicate the magnitude of the different band gaps. Hollow symbols indicate the bulk band gaps. (b) Position of the band edge with respect to the vacuum level for the VB at theK‐point (orange crosses), VB atΓ(red squares), CB at theK‐point (black circles), and CB atΛ(green diamonds) [29]. (Copyright 2014 by the American Physical Society. Reprinted with permission).

2.1.1. Massive Dirac fermions

To gain insight of physics around the K‐ and K‐points, one can reduce multi‐band tight‐binding model to a two band k.pmodel, using Löwdin partitioning method [38]. For the monolayer TMDCs, one gets the Hamiltonian in the first order of kapproximation,


where Δ=1.6eVis the band gap, t=1.1eVis the effective hopping parameter, λ=0.075eVis the SOI parameter, a=3.19Åis the lattice parameter [39], τ=±1is the valley index. The eigenvalues and eigenvectors can be derived straightforwardly as follows:



|c,k τ,sz=|sz(cos(ϑn2)τsin(ϑn2)eiτφk)|v,k τ,sz=|sz(τsin(ϑn2)eiτφkcos(ϑn2))E3

The energy dispersion around the K‐ and K‐ points, described by Eq. 2, is shown in Figure 5. In order to see the reliability of the kpapproach, Figure 6 plots the energy spectrum of monolayer MoS2calculated by the first principle and kpmodel. It can be found that in the vicinity of the K–(K)–point, they have very good agreement.

Figure 5.

Energy dispersion in (a)K‐valley and (b)K'‐valley obtained by thekpmodel.

Figure 6.

(a) Quasiparticle band structure of monolayer MoS2 calculated by the density functional theory (DFT). (b) A blowup of the rectangular red area in (a). The blue (solid), red (dashed), green (dotted), and purple (dot‐dashed) curves correspond to the results obtained by the DFT, andkptheory of the first order, second order, and third order, respectively. Notice that around theK‐point, all methods give almost identical results.

2.1.2. Landau levels of monolayer MoS2

For a perpendicular magnetic field applied to the MoS2sheet, we use Peierls substitution KiΠi=Ki+(e/)Ai, where eis the elementary charge. In the Landau Gauge A=(0,Bx), we define the operators Π±=τΠx±iΠy, which have the following properties:


where the magnetic length is lB=/eB(25.6/B). Using these operators, the destruction and creation operators are introduced in the following way:


in the K‐(τ=1) valley, and


in the K(τ=1) valley. Then, the Hamiltonian in the presence of a perpendicular magnetic field can be well described by


It is worth to mention that since the Zeeman effect is vanishingly small (<5meV), it is neglected. After some algebra calculations, the Landau levels are obtained as follows:


where nis an integer and n1, ωc=2/lB. The corresponding Landau fan diagrams of the monolayer MoS2are shown in Figure 7. The corresponding eigenfunctions are given by

Figure 7.

Conduction band Landau levels of monolayerMoS2in the vicinity ofK‐ (a) andK‐ (b) valleys with spin‐orbit interaction (SOI). (c) and (d) are the same as (a) and (b), but for the valence band. The blue and red lines correspond to the spin‐up and spin‐down Landau levels, respectively.






The eigenfunctions can be written in a compact form as,


For the special case in which n=0, the eigenvalues become


and corresponding eigenfunctions turn out to be


2.2. Optical selection rules

In monolayer TMDCs, both the top of valence bands and the bottom of conduction bands are constructed primarily by the d‐orbits of the transition metal atoms. The giant spin‐orbit coupling splits the valence bands around the K(K') valley by 0.5eV, for MoS2while the conduction band splitting is neglectable. In addition, time reversal symmetry (TRS) leads to the opposite spin splitting at the K‐ and K'‐valleys. Namely the Kramers doublet (K,) and (K',) are separated from the other doublet (K',) and (K, ) by the spin–orbit interaction (SOI) splitting, as shown in Figure 8.

Figure 8.

Schematic illustration of optical transition rules of the valley and spin in theK‐ (the left) andK‐ (the right) valleys. The red (blue) color represents spin‐up (down) states.

We assume that the monolayer TMDCs are exposed to light fields with the energy ωand wave vector kl, which is orthogonal to the monolayer plane and much smaller than 1/a. Up to the first order approximation, the light‐matter interaction Hamiltonian is described by,


with the light field Ain=A0α^cos(klrωt)and v=(1/)kHbeing propagation velocity [15]. Here, A0and α^stand for the amplitude and orientation of the polarization field, respectively. For an electron being excited by an incident photon from its initial state |ito a final state |f, the transition probability is given by the Fermi's golden rule as, Wfi=2π/|f|HLM|i|2n(E), where n(E)is the density of states available for the final state. Since the absorption intensity Iis proportional to the transition rate, it can be evaluated by,


where Ψc(Ψv)is the conduction (valence) band wave function, Λ=γ/π{[ω(Ec(mc,nc)Ev(mv,nv))]2+γ2}, and ϒ=fcfv, with fiFermi‐Dirac distribution function, nithe principal quantum number, mithe quantum number associating with orbital angular momentum, i=c, and vreferring to conduction and valence band, respectively, mc=mand mv=m, γa parameter determined by the Lorentzian distribution.

For a circularly polarized (CP) light, α^=(1,cos(ωtσπ/2),0)T, with σ=±1denoting the corresponding positive and negative helicities and T stands for the transpose of a matrix, then the perturbed Hamiltonian becomes,


with W^σ=etaA0(τσx+iσσy)/2. For example, for τ=1,sz=±1, the transition rate for the monolayer TMDCs is determined by,




The optical transition rate for τ=1and sz=±1can be obtained by replacing + with −in the HLM. From Eqs. (29) and (30), notice that under CP light excitation, a valley and spin polarized emission or absorption light is expected in monolayer MoS2, as shown in Figure 8. In contrast, linearly polarized light does not present valley‐selected emission and absorption spectra because both K‐ and K‐ valleys absorb light simultaneously.

Figure 9.

Reflection and photoluminescence spectra of ultrathinMoS2layers. (a) Reflection difference due to an ultrathinMoS2layer on a quartz substrate, which is proportional to theMoS2absorption constant. The observed absorption peaks at 1.85 eV (670 nm) and 1.98 eV (627 nm) correspond to the A and B direct excitonic transitions with the energy split from valence band spin‐orbital coupling. The inset shows the bulkMoS2band structure neglecting the relatively weak spin‐orbital coupling, which has an indirect bandgap around 1 eV and a single higher energy direct excitonic transition at theKpoint denoted by an arrow. (b) A strong photoluminescence is observed at the direct excitonic transitions energies in a monolayer, respectively [30]. (Copyright 2010 by American Chemical Society, reprinted with permission).

2.3. Valley polarized photoluminescence and excitonic effects of the monolayer TMDCs

In monolayer TMDCs, strong Coulomb interactions due to reduced screening and strong 2D confinement lead to exceptionally high binding energies for excitons [23, 24, 36], which allow them be able to survive even at room temperature. Hence, the typical absorption spectra are usually characterized by strong excitonic peaks marked by A and B, located at 670 and 627 nm, respectively. The strong spin‐orbit interaction in the valence band gives rise to a separation between them, as shown in Figure 9. In addition, an injection of electrons into the conduction band of MoS2, which can be realized by gate‐doping [26], photoionization of impurities [28], substrates [25] or functionalization layers [22, 27], leads to the formation of negatively charged excitons (X). The peak of the Xis positioned at a lower energy side of neutral exciton with a binding energy about 36meV for MoS2, see the peak indicated by Xin Figure 10. In addition, the emergence of the charged exciton is accompanied by a transfer of spectral weight from the exciton. Therefore, the intensity ratio between a neutral and charged exciton can be tuned externally. Besides, with increasing the nonequilibrium excess electron density, a red‐shift of the excitonic ground‐state absorption due to Coulomb‐induced band gap shrinkage occurs. It is also worth to point out that on the one hand the trion can provide a novel channel for exciton relaxation, and on the other hand, it can also be excited by an optical phonon into an excitonic state to realize an upconversion process in monolayer WSe2.

Figure 10.

Photoluminescence spectra (PL) of monolayerWSe2at 50 K for pulsed excitation under applied pump fluences of0.8μJcm2(black curve) and12μJcm2(red curve), respectively. The spectra are normalized to yield the same emission strength for the neutral exciton [40]. (Copyright 2015 by the Nature Publishing Group. Reprinted with permission.).

In the regime of high exciton density, the exciton‐exciton collision leads to exciton annihilation through Auger process or formation of biexciton in the monolayer TMDCs. The biexciton is identified as a sharply defined state in the PL, see P0in Figure 10 and also XX‐peak in Figure 11. The nature of the biexcitonic state is supported by the dependence of its PL intensity on the excitation laser power. At low excitation laser intensity, the peaks P0and Xgrow superlinearly and linearly with incident laser power, whereas they increase sub‐quadratically and sublinearly with the laser power at sufficiently high laser fluence. The large circular polarization of P0emission provides a further support for this assignment.

Figure 11.

Circularly polarized photoluminescence (PL) spectra of monolayerWSe2excited by near‐resonant circularly polarized radiation at 15 K. (a) PL for low exciton density with continuous wave excitation at a photon energy of1.92eV. (b) PL for high exciton density with pulsed excitation at a photon energy of1.82eV. Blue and red curves correspond to the same and opposite circularly polarized states. The emission energies for neutral (X) and charged (X) excitons and the biexciton (XX) state are indicated by dashed lines. Inset shows schematic diagram of energy dispersion inK‐ andK‐valleys and valley‐polarized emission. The vertical arrows indicate the electron spin direction. The circles represent conduction and valence band electrons [40]. (Copyright 2015 by the Nature Publishing Group. Reprinted with permission).

The polarization of the photoluminescence from the TMDCs, which is defined by η=(Iσ+Iσ)/(Iσ++Iσ), inherits that of the excitation source, where Iσ+(Iσ) is PL intensity of right (left) hand circularly polarized light. Figure 11 illustrates photoluminescence spectra of monolayer WSe2excited by near‐resonant circularly polarized radiation at 15K. Notice that the peaks for X, X, and XXall exhibit significant circular polarization. In addition to biexciton emission, the Xand Xemission bands also exhibit strong valley polarization.

2.4. Defect induced photoluminescence and single photon source

As known, vacancy defects, impurities, potential wells created by structural defects or local strain or other disorders might be introduced in the growth process of the TMDC materials [12, 13, 32, 33]. They can produce localized states to participate the optical emission and absorption as manifested by P1to P3in Figure 10, and the emission bands on the lower energy side of the peak XX in Figure 11. Since the point defects can induce intervalley coupling, the defect‐related emission peaks show no measurable circular polarization character. Besides, the excitons, trions, and even biexcitons can be trapped by theses crystal structure imperfections to form corresponding bound quasiparticles. Therefore, delocalized excition, charged exciton and biexciton emissions, and localized ones can coexist in the TMDCs. Interestingly, these carrier trapping centers can act as single‐photon emitters to emit stable and sharp emission line [12, 13, 32, 33]. For this kind of single quantum emitter, since the maximum number of emitted single photons is limited by the lifetime of the excited state, a saturation of the PL intensity at high excitation laser power is expected.

2.5. Magneto‐optical properties of the monolayer TMDCs

The presence of a magnetic field induces a quantization of the energy levels. At high magnetic field, the Landau levels (LLs) form. The transition rate between the conduction–and valence–band Landau levels can be calculated using the eigenfunctions in Eq. 21.




In the presence of magnetic field, both σ+and σabsorptions take place in each valley. However, there is a great difference in their intensity. For instance, the absorption spectrum intensity of σ+‐light is 104times larger than σ‐light in the K–valley. And, the lowest transition energy absorption spectrum demonstrates a valley polarization, i.e., σ+in the K–valley and σin the K–valley, as showed in Figure 12 [41]. In addition, the transition occurs from n=0to n=1LLs in the K–valley, whereas n=1to n=0LLs in the K–valley. Therefore, the valley polarization remains in the magneto‐optical absorption, as showed in Figure 13. It is worth to argue that the higher‐order terms in the effective k.pmodel only induces about 0.1% correction to absorption spectrum intensity [41], which allows us neglect them safely.

Figure 12.

Selection rules for the interband transitions between Landau levels inK‐valley (the left panel) and inK‐valley–valley (the right panel) of monolayerWS2subjected to a magnetic field along thez^‐direction, excited by a circularly polarized light. The blue and red arrows correspond toσ+andσabsorptions, respectively [41]. (Copyright 2013 by the American Physical Society. Reprinted with permission).

Figure 13.

σ+absorption spectrum of monolayerMoS2for spin–up states in theK‐valley (the left panels) andσabsorption spectrum for spin–down states in theK‐valley (the right panels), for a magnetic fieldB=2((a) and (d)),10((b) and (e)),15T((c) and (f)), respectively.

2.6. TMDC quantum dots and valley polarized single‐photon source

The Hamiltonian of the TMDC QDs in polar coordinates is given by [15]


As a matter of convenience, we get rid of the angular part by using the following ansatzfor the eigenfunctions


where the quantum number m=jτ/2, jis the quantum number related to the effective angular momenta Jeffτ=Lz+τσz/2with Lzbeing the orbital angular momenta in the z^direction. Then the Schrödinger equation for TMDC QDs, showed in Figure 14, becomes

Figure 14.

(a) Schematic of a monolayerMoS2QD with radius R. (b) Top view ofMoS2crystal structure.


After some algebra calculations, we get two decoupled equations. They are




where χ=(2EΔ)(Δ+2E2λsτ)4t2a2. These two second order differential equations can be straightforwardly resolved. Finally we obtain,




where N¯is the normalization constant and Jn(x)is Bessel Function of the first kind. Applying the infinite mass boundary condition ψ2(R,θ)ψ1(R,θ)=iτeiτθ, we obtains the secular equation


where Ris the QD radius, see Figure 14.

From Figures 15 and 16, we see that the bound states formed in a single valley, and Eτ(j)Eτ(j)for both conduction– and valence– bands. In addition, the electron‐hole symmetry is broken. Due to the confinement potential, the effective time reversal symmetry (TRS) is broken within a single valley, even without the magnetic field, similarly to graphene QDs [42]. On the other hand, the inverse asymmetry of the crystalline structure, the terms of spin‐orbit interaction and the confinement potential, do not commute with the effective inversion operator, defined in a single valley Pe=Iτσx. Consequently, the QDs do not preserve the electron‐hole symmetry in the same valley. However, comparing Figure 15(a, b), we can found that Eτ(j)=Eτ(j)was still true. This is attributed to the TRS, where THT1=H, the TRS operator is defined as T=iτxsyC, with Cbeing the conjugate complex operator.

Figure 15.

(a) Wave function profileΨnand its two components (an,bn) of the lowest two (n=1, 2) conduction band states in theK–valley with spin–up atj=1.5. The pink circle indicates that the wave functionΨnis nonzero (even though small) atr=R. Energy spectrum of the lowest four (n = 1, 2, 3, 4) conduction bands (solid curves) and the highest four valence bands (dashed curves) as a function of j, in theK‐valley with spin‐up (b) andK‐valley with spin­down (c) of aMoS2dot. The dot size isR=40nmand the magnetic field isB=0[14]. (Copyright 2016 by the IOP publishing. Reprinted with permission).

Figure 16.

Energy spectrum of monolayerMoS2QD with a radius of 40 nm for the conduction–band in theK‐ (a) andK‐ (b) valleys. (c) and (d) are corresponding figures for the valence–band. The blue and red curves correspond to spin–up and spin–down energy levels, respectively.

2.7. Landau levels in monolayer MoS2quantum dots

Similarly to what we did in the case of the monolayer TMDCs, for quantum dot subjected to a perpendicular magnetic field, we do the Peierls substitution and use now the symmetric gauge, i.e., A=(By/2,Bx/2). The part of magnetic field Hamiltonian, in polar coordinates is given by


Then, the total Hamiltonian becomes H+HB, which leads to the following two‐coupled differential equations:


In order to solve this eigenvalue problem, let us first decouple these two equations into






Solving these equations, we obtain the following two components of the eigenfunctions




S(x)is the sign function, 1F1(a,b,x)is the confluent hypergeometric function of the first kind, εis an arbitrary constant (0<ε<1) used to avoid a singularity in the function S(x).

With the eigenfunctions at hand, we can derive the secular equation for the eigenvalues by applying infinite mass boundary condition, i.e.,


Figure 17(a) illustrates the energy spectrum of the lowest four spin­up conduction bands in the K‐valley (τ=1, s=1) as a function of magnetic field (B), for the orbital angular momentum m=±0,1,±2, in the 70 nm dot. At zero magnetic field, the electronic shells of an artificial atom such as s, p, and d shells emerge. In a certain valley, say the K‐valley, the atomic states possess both spin‐ and orbital–degeneracies such as E(m)=E(m). In addition, we should emphasize that the energy spectrum in different valleys have time reversal symmetry at B=0. However, the presence of magnetic field breaks down this symmetry and leads to splittings of the atomic orbitals of the dot. Moreover, in the regime of weak magnetic fields, unlike the monolayer MoS2in which the energy shows linear B response, valley dependent energy levels with nonlinear B response are observed. Such effect is attributed to the competition of the QD confinement with the magnetic field effect [14].

Figure 17.

Energy spectrum of the lowest four conduction band states with spin–up in theK‐valley (a) and spin­down in theK‐valley (b), and the highest four valence band states with spin–up in theK‐valley (e) and spin‐­down in theK‐valley (f) of a 70 nmMoS2QD, as a function of the magnetic field, for angular momentumm=0(red curves),1pink curves),1(green curves),2(black curves), and2(blue curves). The corresponding analogues for a 40 nm dot are shown in (c), (d), (g), and (h), respectively [14]. (Copyright 2016 by the IOP publishing. Reprinted with permission).

As B increases, an effective confinement induced by the magnetic field gradually becomes comparable to that of the dot. Hence, their contributions to the electronic energy are balanced. With a further increasing of B, magnetic field effect starts to dominate the features of the energy spectrum. Accordingly, the LLs which show a linear dependence on B, became of the heavily massive Dirac character, are formed just like in the pristine monolayer MoS2. The higher the energy level is, the stronger the magnetic field needed to form the corresponding LL. For instance, the lowest LL is formed around a critical value B=Bc=2T. Interestingly, in the Hall regime, an energy locked (energy independent on B) mode referring to the lowest spin‐down conduction band in the K‐valley emerges, as shown in Figure 17(b). It is expected to be an analog to the zero energy mode in gapless graphene, associating with certain topological properties. These novel features of the QD energy spectrum are tunable by QD size [14].

Figure 17(c, d) are the corresponding analogs of Figure 17(a, b), but for a dot with R=40nm. Comparing with the 70nm dot (Bc=2T), here the energy locked mode takes place until Bc=4T, which turns out to be larger than the one (2 T) for the 70 nm dot, indicating the locked energy modes arise from the competition between the dot confinement and the applied magnetic field [14].

Let us turn to the energy spectrum of the valence band in the dot of R=70nm. Figure 17(e, f) show the B field dependence of the energy spectrum of the highest four valence bands for several values of angular momentum ms with m=0,±1,±2in the K‐valley with spin‐up (τ=1,s=1) and K‐valley with spin­down (τ=1,s=1), respectively. In contrast to the conduction band, here we find E(m)E(m)even at B=0due to the spin–orbit coupling, see Figure 17(a, b, e, f). Besides that, we also observe the emergent locked energy mode around Bc=2T, but referring to the highest valence band. We should emphasize that the locked energy modes for the conduction and valence bands appear in distinct valleys with opposite spin, see Figure 17(b, e). The corresponding valence band energy spectra for the 40 nm dot are shown in Figure 17(g–h), where the locked energy mode associated with a larger Bc=4T, also arises from the combined effect of the dot confinement and the applied magnetic field, similar to the conduction band. Therefore, the flat band or energy locked modes appear only in the valence band of the K‐valley and the conduction band of the K‐valley [14]. The effect of changing the magnetic field direction is showed in Figure 19.

Figure 18.

Energy spectrum of the conduction band states with spin–up in theK‐valley (a) and spin‐down in theK‐valley (b) of a monolayerMoS2(dotted curves) and of aMoS2QD with radiusR=70nm, as a function of the magnetic field.

Figure 19.

Energy spectrum of the states with spin‐up in theK‐valley (a) and spin­down in theK‐valley (b) of a monolayerMoS2QD with radiusR=70nm, as a function of the magnetic field along both positive and negativez^–direction, for angular momentumm=0(red curves),1(pink curves),1(green curves),2(black curves), and2(blue curves).

A comparison of the energy spectrum of the 70 nm dot with that of the bulk TMDC (i.e., infinite geometry) is shown in Figure 18. Because of the large effective mass at the band edges, the LLs of the bulk TMDC scale as E±(ωc,n)=λτs2±(Δλτs)24+t2a2ωc2nin the low energy region, which resembles conventional 2D semiconductors more than Dirac fermions. The n=0LL appears only in the conduction band of K‐valley (and the valence band of K‐valley, not shown), implying the lifting of valley degeneracy for the ground state. As magnetic field increases, there is an evolution of the energy spectrum from atomic energies to Landau levels. More specifically, in the regime of weak B field, the atomic structure emerges, where the energy is distinct from that of the bulk case. On the other hand, in the strong field regime, the energy in the bulk TMDC and quantum dot becomes identical because the effective confinement due to the magnetic field dominates physical behaviors of the dot. Note that although the two ­band model we used here is widely adopted in the literature [39], the model itself still has limitations, e.g., it cannot properly describe the spin splitting of the conduction band, the trigonal warping of the spectrum and the degeneracy breakdown by applied magnetic field for Landau levels with the same quantum number. However, for usual experimental setups, these effects in the vicinity of the K or Kvalley in which we are interested play a minor role. Hence it can be safely neglected.

2.8. Optical selection rules in monolayer MoS2quantum dots

In the QDs as demonstrated in Figure 20, the optical transition matrix elements in the QDs are computed by,

Figure 20.

Schematic of a monolayerMoS2circular quantum dot with radius R indicated by a red circle, excited by a light field. A magnetic field B is applied perpendicularly to theMoS2sheet [15]. (Copyright 2017 by the Nature Publishing Group. Reprinted with permission).


where szv(szc) denote the valence (conduction) band spin state, Rσ=0Rbc*avrdrand Rσ=0Rac*bvrdrwith ac/vand bc/vthe radial components of the conduction/valence band spinor. The selection rule for optical transitions in TMDC QDs is defined by mvmc=±τand szv=szc, i.e., the angular momentum of the initial and final states differs by ±1, but the spin of these two states is the same. The magnitude of transition rates is determined by the integral Rσand Rσfor the transitions taking place in the valley τ=σand σ, respectively. Since |ac(r)|>|bc(r)|and |av(r)|<|bv(r)|, the integral Rσis much smaller than the Rσ. As a consequence, the absorption in the valley τ=σis stronger than that in the τ=σ, which leads to a valley selected absorption. In a special case in which one component of the wave function spinor is equal to zero such as nll=0LL, only photons with the helicity σ=τis absorbed. Then, one obtains a dichroism η=1. For the linear polarized light, however, there is no valley polarization in the absorption spectrum, in contrast to the case of CPL [15].

In the 2D bulk MoS2, the bottom of the conduction band at the two valleys is characterized by the orbital angular momentum m=0. In contrast, at the top of the valence band, the orbitals with m=2in the K‐valley, whereas mis equal to −2 in the K‐valley. The valley dependent angular momentum in the valence band allows one to address different valleys by controlling the photon angular momentum, i.e., the helicity of the CP light. Such valley‐specific circular dichroism of interband transitions in the 2D bulk has been confirmed [19, 43, 44]. The question is whether its counterpart QD also possesses this exotic optical property. Figure 21(a, b) depict zero‐field band‐edge optical absorption spectrum of a 70‐nm dot pumped by the CP light field. Interestingly, one notices that (i) the polarization of the absorption spectrum is locked with the valley degree of freedom, manifested by the intensity of absorption spectrum with σ=τbeing about 106 times stronger than that with σ=τ, and (ii) the spectrum is spin‐polarized. Thus, the QDs indeed inherit the valley and spin dependent optical selection rule from their counterpart of 2D bulk MoS2. In spite of the distinction in the spin‐ and valley‐polarization of absorption spectra in the distinct valleys, their patterns are the same required by the time reversal symmetry. In Figure 21(c, f) we show the zero‐field optical absorption spectrum as a function of excitation energy for several values of dot‐radius within R=20–80 nm. The involved transitions in Figure 21(c, f) lagged by the numbers have been schematically illustrated in Figure 21(g). Alike conventional semiconductor QDs, several peaks stemmed from discrete excitations of the MoS2QD are observed. As the dot size decreases, the peaks of absorption spectrum undergo a blue shift. In other words, a reduction of the dot size pushes the electron excitations to take place between higher energy states, as a result of the enhancement of the confinement on carriers induced by a shrink of the dot. Therefore, the spin‐coupled valley selective absorption with a tunable transition frequency can be achieved in QDs by varying dot geometry, in contrast to the 2D bulk where a fixed transition frequency is uniquely determined by the bulk band structure. In addition to the blue shift of the transition frequency, the absorption intensity can also be controlled by dot geometry due to size dependence of the peak interval. In fact, an increasing of dot size results in a reduction of the energy separation among confined states. Thus, as the dot size increases, the absorption peaks get closer and closer, see Figure 21(c–f). Eventually, several individual absorption peaks merge together to yield a single composite‐peak with an enhanced intensity. For instance, for the 20‐nm dot (Figure 21c), the lowest energy peak is generated by only one transition, labeled by (1) (see also Figure 21g). However, for the dot with R=80nm, we observe a highly enhanced absorption intensity labeled by (1+2+3), which in the 20‐nm dot refers to three separated peaks with weak absorption intensities indicated by (1), (2), (3), respectively.

Figure 21.

Zero‐field optical absorption spectrum associated with interband transitions between conduction and valence‐band ground states of a 70‐nmMoS2dot for the spin‐up state in theK‐valley (a) and for the spin‐down state in theK‐valley (b), pumped by both clockwise (σ+, blue curve) and anticlockwise (σ, red curve) circularly polarized light fields. (c–f) Absorption intensity for the spin‐down state in theK‐valley under the excitation ofσ, for QDs ofR=20, 35, 50, and 80 nm, respectively. (g) Schematic diagram of the involved interband transitions in (c‐f) tagged by the numbers. The color of the curves is used to highlight the principle quantum number (n) of transition involved states [15]. (Copyright 2015 by the Nature Publishing Group. Reprinted with permission).

2.9. Magneto‐optical properties of monolayer MoS2quantum dots

With the knowledge of the energy spectrum and eigenfunctions of the QD, we are ready to study its magneto‐optical properties. The optical transition matrix in Eq. 27 is applicable to the current case provided that we use newly obtained wavefunctions presented in Eq. 47 and Eq. 48.

Figure 22 shows the magneto‐optical absorption spectra for the spin‐down states in the K‐valley of a monolayer MoS2QD with R=40nm excited by left‐hand circularly polarized light σ, for several values of magnetic field ranging from 0to 15T. A few interesting features are observed. Firstly, the magneto‐optical absorption is also spin‐ and valley‐dependent, as it does in optical absorption spectrum at zero field. In particular the lowest transition energy absorption peak related to the interband transition involving nll=0LL is totally valley polarized with dichroism equal to 1. Thus, the polarization of magneto‐optical absorption locks with the valley. Secondly, for a fixed value of dot size, increasing (decreasing) the strength of the magnetic field results in a blue (red) shift in the absorption spectrum. This arises from the fact that the magnetic field induces an effective confinement characterized by the magnetic length lB=/eB, which is more pronounced for a stronger magnetic field. In addition, the magnetic quantization induced by the magnetic field favors the QD to absorb photons with higher energies. The stronger the magnetic field, the greater the capacity of the QDs to absorb the photons of higher energy. Thirdly, in the high magnetic‐field regime the absorption intensity can be highly enhanced by increasing the strength of magnetic field due to increased degeneracy of the LLs.

Figure 22.

Absorption spectrum for the spin‐down states in theK‐valley of a monolayerMoS2QD withR=40nm excited by left‐hand circularly polarized lightσ, at the magnetic fieldB=0(a), 2T (b), 10T (c), 15T (d), respectively. The corresponding enumerated optical transitions are schematically shown inFigure 21(g)[15]. (Copyright 2015 by the Nature Publishing Group. Reprinted with permission).

2.10. Excitonic effect in monolayer MoS2quantum dots

The optical and magneto‐optical absorptions that we have discussed so far are based on the independent electron‐hole picture. In reality, there is a strong Coulomb interaction between the electron and hole in MoS2. The full treatment of the Coulomb interaction using many‐body theory is beyond the scope of this chapter. However, the excitonic effects in absorption spectrum of the monolayer MoS2QDs can be properly addressed by an exact diagonalization, which is adopted in this work. The Hamiltonian used to describe the exciton is, H(re,rh)=He(re)+Hh(rh)+Veh(rerh), where He(h)is the single electron (hole) Hamiltonian in QDs (see Eq. 1) and Vehis the electron‐hole Coulomb interaction given by Veh(|rerh|)=(1/4πϵrϵ0)(e2/|rerh|). Here ϵ0is the permittivity, ϵris the dielectric constant, and reand rhrespectively stand for the position of electron and hole. An exciton can be understood as a coherent combination of electron‐hole pairs. Thus, the wavefunctions of an exciton can be constructed based on a direct product of single‐particle wave functions for the electron and hole. Since the electron‐hole Coulomb interaction is larger than or at least in the same magnitude as the quantum confinement in the MoS2QDs, the excitonic effect play an very important role. Hence the single‐particle wavefunctions should be modified due to the presence of the Coulomb interaction. Therefore, to get quick convergence of numerical calculation, instead of using simple single‐particle wavefunction, we use the modified one to construct the exciton states, i.e., χj(re,h)=Nexp(re,h/rb)Ψj(re,h), where the exponential factor is a hydrogen‐like s‐wave state, Nis the normalization constant, Ψjis the wave function of the Hamiltonian He,h, and rbis the exciton bohr radius, which has the value of rb1nm in MoS2[24]. Then, the exciton wave function Ψexccan be straightforwardly written as,


with the superscript νreferring to the ν‐th exciton state. Within the Hilbert space made of states of electron‐hole pairs {χi(re)χj(rh)}, the matrix element of Coulomb integral reads,


which involves the wave function χiof the electron and hole in the absence of Coulomb interaction, with j=p,r,s,qindicating the state index. We should emphasize that Eq. 55 contains both the direct interaction between the electron and hole, i.e., Eq. 55 at p=qand r=s,


and the exchange interaction, i.e., Eq. 55 at s=pand r=q,


To calculate Coulomb interaction defined in Eq. 55, we expand 1/|rerh|in terms of half‐integer Legendre function of the second kind Qm1/2, i.e.,


which is widely used in many‐body calculations [45]. Here θe(h)is the polar angle of the position vector re(h)in the 2D plane of MoS2, ξ=(re2+rh2)/2rerh, ϵm=1for m=0and ϵm=2for m0. The exciton energy and the corresponding wavefunction can be obtained by diagonalization of the matrix of many‐particle Hamiltonian H(re,rh). With the exciton state at hand (Eq. 54), we are ready to the determine the excitonic absorption involving a transition from the ground state |0to exciton state |f=|Ψexc,


where Eexcνand ωare the exciton‐ and photon‐energy, respectively, and P=i,jδsz,szχi|HLM|χjai,szhj,szis the polarization operator, with ai,szand hj,szbeing the electron and hole annihilation operators. After straightforward calculation, we finally obtain


In our numerical calculations, we have used five modified single‐particle basis functions with angular momentum ranging from −2.5 to 1.5. Figure 23 shows that there is an exciton absorption peak located at around 550 meV (i.e., exciton binding energy) below the band‐edge absorption. And, the excitonic absorption peak shifts monotonically to higher absorption energy as the dot size is increased. Above the band gap, however, the spectrum is similar to what we found in previous sections in the band‐to‐band transitions using the independent electron‐hole model. Since the exciton absorption peak is far away from the band‐edge absorption, one can in principle study them separately. And, the Coulomb interaction between electron‐hole pair does not change the valley selectivity and our general conclusion. Finally, it is worth to remarking that one can shift the excitonic absorption peak to a higher absorption energy, by varying the band gap parameter (Δ) in our model Hamiltonian.

Figure 23.

Excitonic effects on zero‐field optical absorption spectra for the spin‐down states in theK‐valley of monolayerMoS2QD withR=20nm (a), 35 nm (b), 50 nm (c), and 80 nm (d) under the excitation of the left‐hand circularly polarized lightσ. The Fermi energy is chosen as 0. The vertical dashed line separates two distinct regions of the absorption spectrum: left‐hand side for excitonic absorption and right‐hand side for single‐particle like absorption. To be consistent with the experimental report of excitonic absorption energy, a right shift of photon energy of 700 meV is made [15]. (Copyright 2015 by the Nature Publishing Group. Reprinted with permission.).

© 2017 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Fanyao Qu, Alexandre Cavalheiro Dias, Antonio Luciano de Almeida Fonseca, Marco Cezar Barbosa Fernandes and Xiangmu Kong (October 25th 2017). Valley Polarized Single Photon Source Based on Transition Metal Dichalcogenides Quantum Dots, Quantum-dot Based Light-emitting Diodes, Morteza Sasani Ghamsari, IntechOpen, DOI: 10.5772/intechopen.70300. Available from:

chapter statistics

1503total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Widely Tunable Quantum-Dot Source Around 3 μm

By Alice Bernard, Marco Ravaro, Jean-Michel Gerard, Michel Krakowski, Olivier Parillaud, Bruno Gérard, Ivan Favero and Giuseppe Leo

Related Book

First chapter

Prologue: Nanorods – Recent Advances and Future Perspective

By Soumen Dhara

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us