Abstract
The recently discovered monolayer transition metal dichalcogenides (TMDCs) provide a fertile playground to explore new coupled spin–valley physics1,2,3. Although robust spin and valley degrees of freedom are inferred from polarized photoluminescence (PL) experiments4,5,6,7,8, PL timescales are necessarily constrained by short-lived (3–100 ps) electron–hole recombination9,10. Direct probes of spin/valley polarization dynamics of resident carriers in electron (or hole)-doped TMDCs, which may persist long after recombination ceases, are at an early stage11,12,13. Here we directly measure the coupled spin–valley dynamics in electron-doped MoS2 and WS2 monolayers using optical Kerr spectroscopy, and reveal very long electron spin lifetimes, exceeding 3 ns at 5 K (two to three orders of magnitude longer than typical exciton recombination times). In contrast with conventional III–V or II–VI semiconductors, spin relaxation accelerates rapidly in small transverse magnetic fields. Supported by a model of coupled spin–valley dynamics, these results indicate a novel mechanism of itinerant electron spin dephasing in the rapidly fluctuating internal spin–orbit field in TMDCs, driven by fast inter-valley scattering. Additionally, a long-lived spin coherence is observed at lower energies, commensurate with localized states. These studies provide insight into the physics underpinning spin and valley dynamics of resident electrons in atomically thin TMDCs.
Similar content being viewed by others
Main
Studies of optical spin orientation and spin relaxation using polarized light have a long and exciting history in conventional III–V and II–VI semiconductors14,15. Early seminal works focused on magneto-optical studies of polarized PL from recombining excitons14, from which spin lifetimes could be indirectly inferred. However, it was the direct observation of very long-lived spin coherence of resident electrons in materials such as GaAs and ZnSe (refs 15, 16)—revealed unambiguously by time-resolved Faraday and Kerr rotation studies—that captured widespread interest and helped to launch the burgeoning field of ‘semiconductor spintronics’ in the late 1990s (ref. 15). With a view towards exploring coupled spin/valley physics of resident electrons in the new atomically thin and direct-bandgap TMDC semiconductors, here we apply related experimental methods and directly reveal surprisingly long-lived and coherent spin dynamics in monolayer MoS2 and WS2.
Figure 1a depicts the experimental set-up. High-quality monolayer crystals of n-type MoS2 and WS2, grown by chemical vapour deposition on SiO2/Si substrates17, were selected on the basis of low-temperature reflectance and PL studies (see Methods). Transverse magnetic fields (By) were applied using external coils. A weak pump laser illuminates individual crystals with right- or left-circularly polarized light (RCP or LCP) using wavelengths near the lowest-energy A exciton transition, which primarily photoexcites spin-polarized electrons and holes into the K or K′ valley, respectively1,2,3,4,5,6,7,8,9,10. Any induced spin and valley polarization is then detected by means of the optical Kerr rotation (KR) or Kerr ellipticity (KE) that is imparted to a linearly polarized and wavelength-tunable probe laser that is incident normally and focused on to the crystal. Either continuous-wave (c.w.) or pulsed pump/probe lasers can be used; both types of experiments will be discussed.
The diagrams in Fig. 1a depict the conduction and valence bands in the K and K′ valleys in a typical monolayer TMDC, along with the relevant spin/valley optical selection rules and scattering processes. In n-type material, the resident spin-up and spin-down electrons in the conduction band of the K (K′) valley have densities n↑ and n↓ (n↑′ and n↓′), which all share a common chemical potential in thermal equilibrium. However, following pulsed photoexcitation and fast (∼10 ps) recombination with photogenerated holes9,10, these resident electron densities may be unequal and out of equilibrium, as depicted. This can arise, for example, from many-body correlations while holes are present18,19, or from unequal nonradiative recombination rates of holes with the optically forbidden electron states (nonradiative processes account for the vast majority of recombination in TMDCs (ref. 5)). Thus, photoexcitation can impart a net spin polarization (Sz = n↑ − n↓ + n↑′ − n↓′) and/or valley polarization (Nv = n↑ + n↓ − n↑′ − n↓′) onto the resident electrons that may remain even after all holes have recombined. (Analogously, weak steady-state photoexcitation can establish a non-equilibrium steady-state polarization of resident electrons.) The intrinsic relaxation dynamics of this polarization, which proceeds without the perturbing influence of the holes, is the principal focus of these studies.
Crucially, any long-lived polarization of the resident electrons can be directly monitored using Kerr spectroscopy. This stands in marked contrast to polarized PL studies, which explicitly require the participation of (and recombination with) a photo-excited hole. Kerr effects depend only on the difference between a material’s RCP and LCP absorption and indices of refraction. According to the selection rules in monolayer MoS2 and related TMDCs (refs 1, 2), RCP light near the lowest-energy A exciton in TMDCs couples to the resident electron density n↑ in the K valley. Similarly, LCP light couples to n↓′ in the K′ valley. Perturbations to the densities n↑ and n↓′ shift their chemical potentials, which change the absorption and refraction of RCP and LCP probe light, particularly at wavelengths near optical transitions. Thus, to leading order, Kerr signals using light near the A exciton are proportional to n↑ − n↓′, to which both spin and valley polarization can contribute, namely (Sz + Nv)/2.
Figure 1b demonstrates optically induced polarization in monolayer MoS2 using c.w. pump and probe lasers. The wavelengths λpump (632.8 nm) and λprobe (675 nm) were chosen to address only the A exciton. At zero applied field, a steady-state KR of ∼15 μrad is induced by the pump onto the probe laser. Surprisingly, this polarization signal is sharply reduced by very small transverse fields By. The narrow and nearly Lorentzian-shaped dependence of the measured KR on By (only ∼7 mT half-width) strongly suggests a long-lived spin polarization of resident electrons. (In contrast, the PL polarization, which probes primarily short-lived excitons, is unchanged over this field range, in agreement with earlier studies7.) These Kerr data are very reminiscent of traditional Hanle-effect studies in conventional semiconductors such as n-type GaAs (refs 14, 20, 21), wherein By dephases an optically injected steady-state electron spin polarization (Sz) owing to spin precession about By. This leads to a Lorentzian dependence of Sz on By, from which the spin lifetime τs can be inferred using the half-width By0 of the Hanle peak if the g-factor ge is known; namely, τs−1 = geμBBy0/ℏ.
It is therefore tempting to associate the data in Fig. 1b with spin dephasing in MoS2 due solely to precession of resident electrons about By, and to infer a spin lifetime. However, in contrast to conventional III–V or II–VI semiconductors, electrons in the high-momentum K and K′ valleys of monolayer MoS2 are theoretically predicted to experience strong spin–orbit coupling, due to the different curvature of the spin-up and spin-down conduction bands, and also to any additional intrinsic splitting Δc of these bands1,2,22,23,24. This coupling can be viewed as a large out-of-plane effective magnetic field Bso (of the order of 10 T, depending on electron density). Importantly, Bso is oriented parallel or antiparallel to , depending on whether the electron resides in the K or K′ valley. Spin-conserving inter-valley electron scattering, which is not forbidden in the conduction band and which is expected to be fast in these materials (γv−1 ∼ 0.1–1 ps) given the low electron mobility and the dominant role of impurity scattering, therefore leads to a rapidly fluctuating effective magnetic field ‘seen’ by electrons.
This fluctuating field alone will not affect (dephase) electron spins that are also oriented along . However, in the additional presence of By, electron spins will precess about the total fluctuating field , which is no longer oriented along . This leads to a valley-dependent spin precession and associated dephasing that is analogous to the momentum-dependent spin precession and dephasing common in conventional semiconductors (for example, the Dyakonov–Perel mechanism14) or to electron spin depolarization in germanium, which is also driven by inter-valley scattering25. A direct and testable consequence is that electron spin relaxation in MoS2 is expected to depend strongly on By—even for small By—in marked contrast to ordinary III–V and II–VI bulk semiconductors.
The narrow Hanle–Kerr data in Fig. 1b are consistent with this scenario. However, to directly measure electron spin relaxation and test this hypothesis, we turn to time-resolved KR studies using ultrafast lasers. Figure 1c shows the measured decay of the optically induced electron polarization in the same MoS2 crystal, using photons with energy just below the A exciton resonance. A key finding of this work is that the polarization decay time is extremely long (∼3 ns) at zero field. This timescale exceeds typical PL recombination times by two to three orders of magnitude9,10, further implicating resident electrons as the source of long-lived polarization. Recent studies of monolayer MoS2 and WSe2 using related techniques did not reveal any long-lived polarization imparted to resident electrons, perhaps owing to elevated temperatures13 or to the use of probe photons with higher energy11,12. The nanosecond electron spin relaxation observed here at zero field may ultimately be limited by intra-valley Dyakonov–Perel or Elliot–Yafet processes26,27 (potentially due to long-wavelength flexural phonons28), or by spin-flip scattering with magnetic impurities.
Crucially, the polarization decay time decreases rapidly with small increasing By, and no prominent spin precession is observed. These two observations directly support the scenario described above, of depolarization due to a rapidly fluctuating Bso driven by fast inter-valley electron scattering, for the following two reasons: first, if Bso were small or did not exist, then pronounced spin precession about By would be observed (as is the case for resident electrons in GaAs, ZnSe, GaN, or CdTe (refs 15, 16, 29, 30)). Second, if Bso were static and not fluctuating, then spins would be effectively pinned along Bso and the small By (≪|Bso|) would have little influence on the long spin decay. Neither of these phenomena were observed. Rather, these data are consistent with a novel spin relaxation mechanism in monolayer TMDCs that is driven by fast inter-valley scattering (rapidly fluctuating Bso), and activated by small By.
A model of coupled spin–valley dynamics captures the underlying physics and reproduces essential features of the data. The phenomenological equation of motion for the electron spin polarizations SK and SK′ in the two valleys reads:
where τ = τ(t) = ±1 is the index for K and K′ valleys respectively (inter-valley scattering events flip τ). The first two terms describe spin precession about an applied field and the internal spin–orbit field . The third term describes intrinsic spin relaxation within a given valley with rate γs, and the fourth term is spin-conserving inter-valley scattering with rate γv, which is expected to be fast and which leads to fast relaxation of any valley polarization. We compute the total spin polarization S = SK + SK′. For large Bso (Ωso ≫ ΩL) and fast inter-valley scattering (γv > Ωso), and assuming initial polarization , the solution is Sz(t) = S0 ∑j=1,2. Here, and the two eigenmodes , where Γv ≡ Ωso2/4γv. There is a critical applied field ΩLc = Γv, below which the modes are purely decaying and above which the modes oscillate (see also Supplementary Fig. 4 and accompanying discussion of the complete model). The Hanle curve is calculated by integrating Sz(t) for each ΩL; namely ∫ 0∞Sz(t)dt = S0(γs + 2Γv)/[ΩL2 + (γs2 + 2γsΓv)]. Thus, Hanle–Kerr data are predicted to be Lorentzian (as observed) with half-width .
Figure 1d, e shows the calculated Hanle curves and electron spin dynamics. The model reproduces the rapid drop in electron spin polarization with increasing (small) By, and captures the shallow dip and subsequent recovery of the electron spin polarization at short timescales. The model does not, of course, capture the offset of the measured Hanle–Kerr data (see Fig. 1b), the origin of which suggests an additional long-lived and field-independent polarization, perhaps from localized states in MoS2, which have been studied recently. This offset is also manifested in the time-resolved data as the very slowly decaying and largely field-independent signal that persists at long delays.
A surprising additional observation is the appearance, at larger magnetic fields, of a small but long-lived oscillatory signal visible between ∼200–800 ps (red trace in Fig. 1c). Figure 1f shows an expanded view of this signal at By = 236, 270 and 304 mT. The oscillation frequency scales linearly with By, indicating that this signal arises from coherently precessing electrons with g-factor |ge| ≃ 1.8. Such long-lived coherence signals are not expected from itinerant resident electrons, for reasons described above. However, they may arise owing to contributions from an additional population of localized states that do not undergo rapid inter-valley scattering and which precess about the bare applied field By. Time-resolved measurements at lower photon energies below the A exciton are consistent with this scenario (Supplementary Fig. 1).
We confirm that these Kerr signals originate from MoS2 by showing their spectral dependence. Figure 2 shows the MoS2 reflectivity, below which are the peak Hanle–Kerr signals at By = 0. Both KR and KE show a strong resonance at the A exciton. Data using 632.8 nm and 543.5 nm c.w. pump lasers are shown; the latter allows the induced signals to be tracked out to the higher-energy B exciton, which elicits a smaller response. This behaviour was confirmed on many different MoS2 crystals. The resonances are redshifted ∼25 meV from the exciton peak, suggesting that the resident electrons’ polarization may be revealed preferentially through trion-related (rather than neutral exciton-related) optical transitions2, analogous to studies in III–V and II–VI semiconductor quantum wells29,30. The dependence of Hanle–Kerr data on probe wavelength is shown in Supplementary Fig. 2.
An essential consideration for spin-coherent phenomena in semiconductors is its dependence on temperature. Systematic studies on MoS2 (see Fig. 3) reveal very long spin lifetimes at 5 K (∼3 ns), which decrease rapidly to <200 ps by 40 K. This behaviour is in qualitative agreement with recent predictions28 of Elliot–Yafet spin relaxation processes (within a given valley) due to electron–phonon scattering with long-wavelength flexural phonons. A comparison of supported versus unsupported monolayer crystals could confirm this mechanism.
Finally, we show that these long-lived electron spin polarizations are not unique to MoS2, but also appear in other TMDCs. Figure 4 shows Hanle–Kerr and time-resolved KR measurements in monolayer WS2, along with their spectral dependence. WS2 also exhibits very narrow Hanle–Kerr signals and nanosecond polarization relaxation times. The peak Hanle–Kerr signal also occurs on the low-energy side of the A exciton, again suggesting that coupling to resident electrons may proceed preferentially through trion transitions. The field-independent background on the Hanle data is significantly larger (about 80% of the peak, rather than 30–40% for MoS2), which may result from lesser material quality or, possibly, the larger spin–orbit splitting in WS2. Together, these measurements open the door for systematic studies of the detailed interplay between spin–orbit coupling, chemical potential, and spin/valley dynamics of resident carriers in new doped two-dimensional TMDC semiconductors.
Methods
High-quality triangular crystals of n-type (electron doped) monolayer MoS2 and WS2 were grown by chemical vapour deposition on SiO2/Si substrates, following ref. 17. Typical lateral dimensions ranged from 10–20 μm. The background (residual) electron doping level in these crystals is estimated to be ne ∼ 5 × 1012 cm−2 on the basis of transport studies of field-effect transistor devices fabricated from similarly grown single-monolayer crystals. The samples were mounted on the vacuum cold finger of a small variable-temperature liquid helium optical cryostat (3–300 K). Applied magnetic fields up to 300 mT were generated using an external electromagnet. Individual crystals were screened for a high degree of circularly polarized exciton PL (>80%) and sharp reflectivity spectra.
All Kerr-effect measurements were performed in a reflection geometry, as depicted in Fig. 1a. For Hanle–Kerr measurements of steady-state polarization, continuous-wave (c.w.) pump and probe lasers were used. The 632.8 nm and 543.5 nm lines of HeNe lasers were used for the pump, whereas the probe laser was typically a tunable narrowband dye laser or fixed-wavelength diode lasers. Time-resolved measurements used ultrafast 250 fs pulses from a wavelength-tunable 76 MHz optical parametric oscillator (OPO). Time-resolved studies typically used wavelength-degenerate pump and probe pulses (an exception is the temperature-dependent data of Fig. 3, where a 635 nm pulsed diode laser with ∼100 ps pulse duration was synchronized to the OPO and used as a non-degenerate pump laser). In all experiments, the pump laser was weakly focused so as to illuminate the entire TMDC crystal, thereby mitigating any possible influence of density gradients or carrier diffusion, whereas the linearly polarized and normally incident probe laser was more tightly focused to a ∼4 μm spot in the centre of the crystal. It was verified that the small off-axis angle of the pump laser (about 10 degrees from sample normal) did not influence the results—similar data were obtained with co-propagating normally incident pump and probe lasers.
Except where otherwise noted, the pump laser polarization was modulated between RCP and LCP by a photoelastic modulator to facilitate lock-in detection. Detection of the induced optical polarization rotation (KR) and ellipticity (KE) imparted on to the probe beam was achieved with a standard optical bridge arrangement using balanced photodiodes. Low-power pump and probe beams were used, as it was observed that the Hanle widths and polarization decay rates increased with increasing pump and/or probe power. Typical average probe power was in the tens of μW, whereas the pump was in the range from 100–1,000 μW, depending on wavelength and temperature.
References
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).
Xu, X., Yao, W., Xiao, D. & Heinz, T. F. Spin and pseudospins in layered transition metal dichalcogenides. Nature Phys. 10, 343–350 (2014).
Mak, K. F., McGill, K. L., Park, J. & McEuen, P. L. The valley Hall effect in MoS2 transistors. Science 344, 1489–1492 (2014).
Splendiani, A. et al. Emerging photoluminescence in monolayer MoS2 . Nano Lett. 10, 1271–1275 (2010).
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).
Zeng, H., Dai, J., Yao, W., Xiao, D. & Cui, X. Valley polarization in MoS2 monolayers by optical pumping. Nature Nanotech. 7, 490–493 (2012).
Sallen, G. et al. Robust optical emission polarization in MoS2 monolayers through selective valley excitation. Phys. Rev. B 86, 081301(R) (2012).
Jones, A. M. et al. Optical generation of excitonic valley coherence in monolayer WSe2 . Nature Nanotech. 8, 634–638 (2013).
Lagarde, D. et al. Carrier and polarization dynamics in monolayer MoS2 . Phys. Rev. Lett. 112, 047401 (2014).
Wang, G. et al. Valley dynamics probed through charged and neutral exciton emission in monolayer WSe2 . Phys. Rev. B 90, 075413 (2014).
Zhu, C. R. et al. Exciton valley dynamics probed by Kerr rotation in WSe2 monolayers. Phys. Rev. B 90, 161302(R) (2014).
Plechinger, G., Nagler, P., Schuller, C. & Korn, T. Time-resolved Kerr rotation spectroscopy of valley dynamics in single-layer MoS2. Preprint at http://arxiv.org/abs/1404.7674 (2014)
Dal Conte, S. M. et al. Valley and spin dynamics in monolayer MoS2. Preprint at http://xxx.tau.ac.il/abs/1502.06817 (2015)
Meier, F. & Zakharchenya, B. P. (eds) Optical Orientation: Modern Problems in Condensed Matter Sciences Vol. 8 (North-Holland, 1984).
Awschalom, D. D., Loss, D. & Samarth, N. (eds) in Semiconductor Spintronics and Quantum Computation (Nanoscience and Technology Series, Springer, 2002).
Kikkawa, J. M., Smorchkova, I. P., Samarth, N. & Awschalom, D. D. Room-temperature spin memory in two-dimensional electron gases. Science 277, 1284–1287 (1997).
Najmaei, S. et al. Vapor phase growth and grain boundary structure of molybdenum disulfide atomic layers. Nature Mater. 12, 754–759 (2013).
Mai, C. et al. Many-body effects in valleytronics: Direct measurement of valley lifetimes in single-layer MoS2 . Nano Lett. 14, 202–206 (2014).
Singh, A. et al. Coherent electronic coupling in atomically thin MoSe2 . Phys. Rev. Lett. 112, 216804 (2014).
Dzhioev, R. I. et al. Low-temperature spin relaxation in n-type GaAs. Phys. Rev. B 66, 245204 (2002).
Furis, M. et al. Local Hanle-effect studies of spin drift and diffusion in n:GaAs epilayers and spin-transport devices. New J. Phys. 9, 347 (2009).
Kormányos, A. et al. Monolayer MoS2: Trigonal warping, the Γ valley, and spin–orbit coupling effects. Phys. Rev. B 88, 045416 (2013).
Kośmider, K., González, J. W. & Fernández-Rossier, J. Large spin splitting in the conduction band of transition metal dichalcogenide monolayers. Phys. Rev. B 88, 245436 (2013).
Tse, W.-K., Saxena, A., Smith, D. L. & Sinitsyn, N. A. Spin and valley noise in two-dimensional Dirac materials. Phys. Rev. Lett. 113, 046602 (2014).
Li, P. K., Li, J., Qing, L., Dery, H. & Appelbaum, I. Anisotropy-driven spin relaxation in germanium. Phys. Rev. Lett. 111, 257204 (2013).
Ochoa, H. & Roldan, R. Spin-orbit-mediated spin relaxation in monolayer MoS2 . Phys. Rev. B 87, 245421 (2013).
Wang, L. & Wu, M. W. Electron spin relaxation due to D’yakonov-Perel’ and Elliot-Yafet mechanisms in monolayer MoS2: Role of intravalley and intervalley processes. Phys. Rev. B 89, 115302 (2014).
Song, Y. & Dery, H. Transport theory of monolayer transition-metal dichalcogenides through symmetry. Phys. Rev. Lett. 111, 026601 (2013).
Zhukov, E. A. et al. Spin coherence of a two-dimensional electron gas induced by resonant excitation of trions and excitons in CdTe/(Cd, Mg)Te quantum wells. Phys. Rev. B 76, 205310 (2007).
Chen, Z., Carter, S. G., Bratschitsch, R. & Cundiff, S. T. Optical excitation and control of electron spins in semiconductor quantum wells. Physica E 42, 1803–1819 (2010).
Acknowledgements
We gratefully acknowledge D. L. Smith and H. Dery for helpful discussions, and W. D. Rice for laser expertise. This work was supported by the Los Alamos LDRD programme. These optical studies were performed at the National High Magnetic Field Laboratory, which is supported by NSF DMR-1157490 and the State of Florida. We also acknowledge the support from AFOSR (grant FA9550-14-1-0268) and the Welch Foundation (grant C1716).
Author information
Authors and Affiliations
Contributions
L.Y. and S.A.C. conceived and built the experiments. W.C., J.Y., J.Z. and J.L. grew the samples. L.Y. performed the optical measurements. N.A.S. provided theoretical insight. L.Y., N.A.S. and S.A.C. wrote the paper in consultation with all authors.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Supplementary Information
Supplementary Information (PDF 752 kb)
Rights and permissions
About this article
Cite this article
Yang, L., Sinitsyn, N., Chen, W. et al. Long-lived nanosecond spin relaxation and spin coherence of electrons in monolayer MoS2 and WS2. Nature Phys 11, 830–834 (2015). https://doi.org/10.1038/nphys3419
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/nphys3419
This article is cited by
-
Critical assessment of G0W0 calculations for 2D materials: the example of monolayer MoS2
npj Computational Materials (2024)
-
Hydrogenic spin-valley states of the bromine donor in 2H-MoTe2
Communications Physics (2023)
-
Twist angle dependent interlayer transfer of valley polarization from excitons to free charge carriers in WSe2/MoSe2 heterobilayers
npj 2D Materials and Applications (2023)
-
Spin filtering effect in all-van der Waals heterostructures with WSe2 barriers
npj 2D Materials and Applications (2022)
-
Ultrafast pseudospin quantum beats in multilayer WSe2 and MoSe2
Nature Communications (2022)