Radiatively limited dephasing and exciton dynamics in MoSe{}_{2} monolayers

Radiatively limited dephasing and exciton dynamics in MoSe monolayers

Tomasz Jakubczyk Université Grenoble Alpes, Institut Néel, F-38000 Grenoble, France    Valentin Delmonte Université Grenoble Alpes, Institut Néel, F-38000 Grenoble, France    Maciej Koperski Laboratoire National des Champs Magnétiques Intenses, CNRS-UGA-UPS-INSA-EMFL, 25 Av. des Martyrs, 38042 Grenoble, France    Karol Nogajewski Laboratoire National des Champs Magnétiques Intenses, CNRS-UGA-UPS-INSA-EMFL, 25 Av. des Martyrs, 38042 Grenoble, France    Clément Faugeras Laboratoire National des Champs Magnétiques Intenses, CNRS-UGA-UPS-INSA-EMFL, 25 Av. des Martyrs, 38042 Grenoble, France    Wolfgang Langbein Cardiff University School of Physics and Astronomy, The Parade, Cardiff CF24 3AA, UK    Marek Potemski Laboratoire National des Champs Magnétiques Intenses, CNRS-UGA-UPS-INSA-EMFL, 25 Av. des Martyrs, 38042 Grenoble, France    Jacek Kasprzak Université Grenoble Alpes, Institut Néel, F-38000 Grenoble, France jacek.kasprzak@neel.cnrs.fr

By implementing four-wave mixing (FWM) micro-spectroscopy we measure coherence and population dynamics of the exciton transitions in monolayers of MoSe. We reveal their dephasing times T and radiative lifetime T in a sub-picosecond (ps) range, approaching T=2T, and thus indicating radiatively limited dephasing at a temperature of 6 K. We elucidate the dephasing mechanisms by varying the temperature and by probing various locations on the flake exhibiting a different local disorder. At a nanosecond range, we observe the residual FWM produced by the incoherent excitons, which initially disperse towards the dark states, but then relax back to the optically active states within the light cone. By introducing polarization-resolved excitation, we infer inter-valley exciton dynamics, showing an initial polarization degree of around 30%, constant during the initial sub-picosecond decay, followed by the depolarization on a picosecond timescale. The FWM hyperspectral imaging reveals the doped and undoped areas of the sample, allowing to investigate the neutral exciton, the charged one or both transitions at the same time. In the latter, we observe the exciton-trion beating in the coherence evolution indicating their coherent coupling.


CNRS, Institut Néel, "Nanophysique et semiconducteurs" group, F-38000 Grenoble, France \alsoaffiliationCNRS, Institut Néel, "Nanophysique et semiconducteurs" group, F-38000 Grenoble, France \alsoaffiliationInstitute of Experimental Physics, Faculty of Physics, University of Warsaw, ul. Pasteura 5, 02-093 Warsaw, Poland \alsoaffiliationCNRS, Institut Néel, "Nanophysique et semiconducteurs" group, F-38000 Grenoble, France


excitons, coherent nonlinear spectroscopy, transition-metal dichalcogenides, MoSe, four-wave mixing, coherent nonlinear spectroscopy


The identification of atomically-thin solids 1 resulted in the development of the intriguing physics of graphene, followed by the emerging technological applications 2, 3. Also, it stimulated a rapid progress in fundamental studies of thin films extracted from other than graphite layered materials, such as, semiconducting transition metal dichalcogenides (S-TMDs). The bandgap of S-TMDs, converts from indirect to direct, when reducing the material thickness to a single monolayer 4, enabling exceptionally strong excitonic transitions. Owing to the breakdown of out of plane translational symmetry for two dimensional systems, the coupling of excitons with light is boosted 5, resulting in their short radiative lifetime and thus increased oscillator strength . Excitons in S-TMD monolayers display large binding energies E of several hundreds of meV 6, 7, 8, 9 - an increase by one to two orders of magnitude with respect to a typical semiconductor quantum well 10, 11. The observed excitonic absorption in S-TMDs reaches the values as high as 10 %12, illustrating an exceptionally strong in these systems and implying a radiative lifetime in the sub-picosecond (sub-ps) range 13, as recently observed in monolayers of WSe14, 15 and other emerging two-dimensional systems, namely nano-platelets 27. The robust coupling with light is attractive in prospective applications, especially for photodetectors 3. It is also appealing in optical fundamental studies, for instance in polaritonics 16 and nonlinear spectroscopy 15, 17, 19.

The large in S-TMDs gives rise to a giant nonlinear optical response, which is investigated in this work via three-pulse four-wave mixing (FWM) micro-spectroscopy. This approach offers a direct access to the exciton coherence and population dynamics with a time resolution limited only by the duration of the laser pulses, which resonantly excite targeted optical transitions. Indeed, time-resolved non-resonant photoluminescence measurements 18, offer insufficient time-resolution to infer sub-ps evolution and involve complex relaxation pathways. Instead, resonant experiments 15, 19, also investigating internal transitions 14 occurring in the THz domain, have recently revealed ultrafast radiative recombination of the exciton ground state in WSe.

In this work, we show that the optical dephasing time T of the neutral exciton transition (EX) in a monolayer of MoSe is intrinsically limited by the EX ultrafast radiative recombination T. In contrast to previous works inferring optical coherence in S-TMDs 15, 19, we take advantage of the micro-spectroscopy approach, employing the laser beams focused down to the diffraction limited size of 0.7 m (full width at half-maximum, FWHM). Such implementation is used to perform the FWM hyperspectral imaging 20, which discriminates the exciton charge state across the sample and helps to reveal the striking features of coherent coupling between the neutral and charged excitons. We demonstrate that the EX linear response inferred via micro-reflectance, is affected by the inhomogeneous broadening also on a sub -m scale, evidenced by the photon echo formation in the FWM transients 21. The analysis of the FWM delay dependence yields T2T, occurring on a sub-picosecond timescale. A spatially-resolved study shows that T depends on the local disorder, generating spatially varying localization potentials for excitons and thus varying , directly influencing T. From the temperature dependent homogenous linewidth, retrieved from the measured coherence dynamics, we determine the dephasing due to phonon interaction, which is described by a linear part and a thermally activated part, as is quantum wells. We also monitor the density dynamics of excitons, which is governed by the interplay between bright states within the light cone and various available dark states.

Sample & experiment

An image of the studied MoSe monolayer is shown in Fig. 2. It was exfoliated from bulk material and transferred onto a Si/SiO substrate, with a nominal SiO thickness of 86 nm. Relatively large monolayer flakes of up to m size have been fabricated and placed in an optical He-flow cryostat. To infer both the coherence and the population dynamics of excitons, we retrieved their FWM by implementing a three-beam configuration 20 of the heterodyne spectral interferometry 21. This technique has been proven as an efficient detection scheme of optical nonlinearities in a solid, until now only employed to retrieve wave mixing signals generated by individual transitions in semiconductor quantum dots 21, 25, 20. FWM is an optical polarization created with short, resonant driving pulses; , and , as depicted in Fig. 1 a. In the third-order regime its amplitude is proportional to . Thus, owing to a large , a dramatic enhancement of FWM is expected in MoSe monolayers. Note that the stands for the complex conjugate, which is the origin of the FWM rephasing in inhomogeneously broadened systems 21, generating the photon echo (as also sketched in Fig. 1 a), which acts as a probe of the microscopic dephasing. Micro-FWM spectroscopy requires the co-linear arrangement of the driving fields , which is enabled by phase-selecting the signal through optical heterodyning 21. By employing acousto-optic modulation, are frequency up-shifted by radio-frequencies, introducing controlled phase-drifts in their respective pulse trains generated by a Ti:Sapphire femto-second laser. As we intend to measure sub-ps dynamics, pulseshaping 20 has been applied to correct the temporal chirp. After acquiring the delays and , introduced by a pair of mechanical delay lines, are recombined into a common spatial mode and are focused on the sample with the microscope objective. The reference beam , used in the heterodyne mixing and interferometric detection, is also focused at the sample, yet is displaced with respect to , as depicted in Fig. 2 a. The time-ordering of the pulses is presented on Fig. 1 a: measuring time-integrated FWM detected at the heterodyne frequency, as a function of (), yields the coherence (population) dynamics of an optical transition. The FWM signal is measured in reflectance, attaining a shot-noise detection limit and rejecting the resonant driving fields with () selectivity in field (intensity). The interference between the heterodyned signal and is spectrally resolved with an imaging spectrometer. Further details regarding the current experimental implementation are given in Ref. [20].

In Fig. 1 b we present a typical spectral interference between and the FWM at a temperature T=5 K for . The FWM intensity is retrieved by spectral interferometry and shown in Fig. 1 c. The signal consists of two transitions identified as the ground state - exciton (EX) and the single electron - trion transitions (TR) 22. In Fig. 1 d we present the FWM intensity of EX as a function of intensity demonstrating a linear dependence in the regime, as expected. Note that the FWM can be driven with a intensity as low as a few tens of nW, corresponding to a few hundreds of photons per pulse and generating a low carrier density, less than  cm. Such density is far below the saturation density 23 estimated at around  cm. In our study, we therefore consider the creation of excitons with K center of mass momentum, which appear at the, K+ and K-, points of the Brillouin zone of the S-TMD crystal. Those excitons can either decay radiatively or disperse out of the light cone (K) via phonon scattering. The latter process occupies dark exciton states, which relax back to K on a pico-second time scale and eventually recombine. Further relaxation pathways are scattering electrons and holes between the K-points.

Hyperspectral imaging

The micro-spectroscopy approach enables to perform FWM hyperspectral imaging 25, 20, as shown in Fig. 2. It allows to identify regions of the flake dominated by the FWM of EX or TR (see Fig. 2 b and c). Clearly, the two images are complementary, which permits to distinguish regions of different resident carrier concentration. In Fig. 2 b we note that the FWM amplitude remains virtually constant over the areas of more than m, indicating weak disorder and thus enabling an extremely fast radiative decay of excitons. In the following experiments, the performed imaging allowed us to selectively address EX or TR, or to drive simultaneously both transitions. Performing micro-spectroscopy also permits to locally address sub -m regions of smaller spectral inhomogeneous broadening with respect to the total area of the flake. To exemplify this, we have performed a statistical analysis of the micro-photoluminescence (PL) hyperspectral imaging. Similarly as in Fig. 1 c, confocally detected PL spectra yield EX and TR transitions, as displayed at the top of Fig. 2 d. We observe particularly bright emission, with spectrally integrated count rate of typically 350 kHz from each transition. For every spatial position, we have determined the integrated intensity for both transitions and their center energies. For the latter, we observe the spread over meV, as displayed at the bottom of Fig. 2 d. Interestingly, despite this large, macroscopic - which we attribute to the strain distribution across the flake - the trion binding energy remains well defined, meV. We note that the lower the TR transition energy, the higher its intensity is measured, which is attributed to the the distribution of residual electrons in the sample. Interestingly, such correlation is not observed for EX, i.e., the EX intensity is not sensitive to the apparent, in our sample changes in the charge density. In Fig. 2 e and f we present the PL imaging spectrally integrated over EX and TR transition, respectively. By comparing it with Fig. 2 b and c, we point out two advantages of the FWM imaging with respect to the PL one. Firstly, owing to the third-power scaling of the FWM intensity with the excitation power, combined with heterodyning with , the spatial resolution in the FWM imaging is enhanced 26 up to  nm, surmounting the standard diffraction limit by a factor of 2. Secondly, the FWM yields a significantly improved imaging selectivity of EX and TR across the flake. In fact, different properties are inferred in both experiments. The FWM directly probes of the resonantly generated excitons at the K-points of the valleys. Instead, the non-resonant PL reflects more complex carrier relaxation along the valleys towards their K-points, prior to the exciton formation, followed by their radiative recombination.

Coherence dynamics

The strength of the FWM spectroscopy in assessing the coherence in solids lies in its capability to separate homogenous () and inhomogeneous () contributions of the transition’s spectral width. In particular, in a presence of , the time-resolved FWM amplitude exhibits a photon echo 15, 27, which decays as . Hence, to investigate the exciton coherence dynamics, we measured degenerate FWM , as a function of . Time-resolved FWM amplitude of the EX transition, displayed in Fig. 3 a, clearly demonstrates formation of the photon echo. From its temporal width we estimate the local to be in the meV range (around 3 meV for the case shown in Fig. 3 a). Fig. 3 b shows the time-integrated FWM as a function of . The data are modeled by a convolution of a Gaussian profile with an exponential decay. The former exhibits the FWHM width of 0.16 ps, reflecting duration of impinging the sample. Instead, from the latter we retrieve at T=6 K the dephasing time T=fs, and thus meV (FWHM). For simplicity, the dynamics owing to the echo formation process close to zero delay, has been here disregarded, yet it explains a minor deviation between the measurement and prediction. The homogenous broadening is around twice smaller than the transition linewidth directly measured via micro-reflectivity, as shown in the inset of Fig. 3 b and also in FWM (see Fig. 1 c). We thus conclude that the line-shape remains affected by the inhomogeneous broadening , even though a sub-m area is probed. With increasing temperature, T is expected to decrease 15, owing to phonon scattering. This is highlighted in Fig. 3 b. At T=45 K the dephasing accelerates and we measure Tfs. At ambient temperature the dephasing occurs at a timescale faster than 100 fs and is not resolved by our setup, although a pronounced FWM is still measured. We note that for the TR transition at T=6 K we find Tfs and similar inhomogeneous broadening as for EX, also generating a photon echo as in Fig. 3 c (not shown). This shorter dephasing of the charged exciton, is tentatively interpreted in terms of the final state damping, due to the energy distribution of final state energies of the leftover electron. In Fig. 3 c we present the coherence dynamics measured at the boundary of the doped and undoped regions of the flake, marked with a cross in Fig. 2 d, such that the FWM of both EX and TR is driven in tandem. We observe a beating, withstanding on both transitions during initial positive delays , up to around 700 fs when the TR coherence virtually vanishes. The beating period of fs (marked with a pair of vertical lines), well corresponds to the EX - TR binding energy meV and thus indicates their coherent coupling 28. To support this statement - and in particular, to distinguish between polarization interference and coherent coupling 25 - the data have been Fourier-transformed along the delay to yield the two-dimensional FWM spectrum, indeed revealing off-diagonal coupling terms between EX and TR (not shown).

To gain a deeper understanding of the exciton ultrafast dynamics, we have combined the enhanced spatial and temporal resolution of our experiment and we have performed spatially-resolved dephasing study at T=6 K. Within an area of 8 m by 8 m, displaying uniquely the EX transition, we have scanned the coherence dynamics with a spatial step of m. The analysis of the obtained statistics of dephasing (169 traces), reveals variations of T within the probed area from around 0.5 ps up to 1.5 ps. Interestingly, the locations on the flake yielding the shortest T also display the broadest photon echo, and thus the smallest . Conversely, the longest T is measured on the areas characterized by a larger , and thus showing the narrowest photon echo, here limited by the temporal duration of the laser pulses. An example of such two representative cases is displayed in Fig. 4 a and b. The spatially-resolved dephasing experiment indicates that T is governed by a local disorder, realizing various localization potentials. The resulting spatially-dependent coherence volume of EX has a direct impact on its radiative lifetime T: the fastest recombination is expected at the areas of the smallest , as indeed measured. Full spatial correlations between T, T and will be reported in a forthcoming publication. In Fig. 4 c we present the coherence dynamics on the chosen area exhibiting an increased T, measured from 6 K to 150 K. The data clearly show a gradual decrease of T with temperature from 1.4 ps to 0.2 ps respectively. The retrieved is plotted in the inset. It can be modeled with a linear dependence and an additional bosonic term 41: . The linear term (meV, meV/K) is due to low energy acoustic phonons. The latter term, with the energy meV and meV, could be attributed to thermal activation of higher energy optical phonons 41.

Population dynamics

In the following, the FWM is employed to infer the EX population dynamics after their resonant and selective excitation in a given valley ( and , denoted as , are co-circularly polarized). They arrive at the flake with virtually no delay, fs, generating exciton population. Owing to the large , the excitons exhibit fast radiative decay, yet also are subject to scattering on lattice defects and phonons. The matrix element for the Coulomb induced, parametric exciton-exciton scattering is proportional to the square of an exciton Bohr radius 29. We thus propose that the Coulomb scattering is reduced for these spatially compact excitons. In particular, this channel is negligible for the small exciton densities employed here. Nevertheless, excitons still experience scattering 24 for instance induced by defects, disorder and phonons, which can efficiently redistribute them out of the light cone (marked with dashed lines in the inset of Fig. 3 c), toward K. This type of scattering is particularly efficient in S-TMDs: due to heavy masses (and thus flat bands at the bottom of the branches) and large , the excitons’ center of mass scatter out of the light cone at practically no cost in energy. As a result, one part of the created population instantly decays with its radiative lifetime T, while the remaining part spreads along the dispersion branch populating dark states (i.e. outside the light cone) or is scattered into one of the counter-polarized valleys involving spin-flip processes. Also a part of these excited excitons eventually relax toward the bottom of the valley with a characteristic time T and recombine radiatively, as depicted in Fig. 3 c. In the dynamics probed with FWM, we disregard the influence of strongly localized, quantum-dot like states, due to their small density (not detected on the region shown in Fig. 2). For simplicity, we also disregard non-radiative Auger processes 14, as we operate at low exciton densities.

The FWM, triggered by from the density grating generated by and , therefore probes all the above mentioned processes via dependence. This is shown in Fig. 3 c, for the EX (obtained the same spatial position as in Fig. 3 b) for driving upon co-circular polarization of . The data are modeled with a double exponential decay convoluted with the Gaussian laser pulse, showing temporal width of 0.16 ps (FWHM). From the initial decay we retrieve the exciton lifetime of Tfs, while the FWM for further delays yields the intra-valley relaxation time of Tps. At T=45 K we measure instead (T, T)=ps. An increase of T with temperature is attributed to Boltzmann distribution of excitons, allowing for reaching higher energies and K-vectors out of the light cone, and also to access different dark states offered by the complex structure of the valley-excitons. A physical picture arising from the FWM experiment at low temperature is that the initial exciton decay with the time T, covering an order of magnitude in amplitude (see Fig. 3 c), is due to the radiative recombination, while the non-radiative processes are of minor impact. Importantly, comparing the T with T times, we conclude that the dephasing is principally due to the radiative decay, nearing to the radiative limit - the data can also be well modeled by fixing T=2T and using T as the only fitting parameter. Observation of the radiatively limited dephasing is a prerequisite for implementing more advanced optical coherent control schemes in S-TMDs.

We point out that the dynamics of the secondary, incoherent excitons - here probed by dependence of FWM - is particularly complex and might be influenced by a set of unconventional features present in S-TMDs. One should bear in mind coupling with and scattering toward counter-polarized K-valleys and spin-split bands, relevant in formation of exciton complexes. In fact, triple-degeneracy of the K-valleys, enables various configurations for bright and dark excitons states, as previously considered for other systems hosting multi-excitons 31, 30. We also point out exotic dispersion relations for the center-of-mass momentum 32. Thus, we refrain from firm interpretation and first-principle modeling of the secondary exciton dynamics, displayed in Fig. 5 on a 100 ps time scale.

Inter-valley dynamics

Below, we present the FWM results obtained upon polarization-resolved driving, employed to investigate exciton scattering between counter-polarized valleys, testing the robustness of the pseudo-spin degree of freedom. The reciprocal space of TMDs monolayers displays nonequivalent bands with the extrema at -points 33, 34, 35, labeled as and . The dipole-allowed transitions in these two valleys can be selectively addressed by circularly polarized light, and , respectively. The light helicity is therefore considered as an asset to drive, manipulate and read the state of the valley subspaces. The valley polarization 34, although protected by the strong spin-orbit splitting in the valence band, decays mainly due to the to electron-hole exchange interaction 36, 37. Since the stability of this degree of freedom is required in prospective applications of TMDs, intense efforts are currently devoted to study its dynamics and to reveal the mechanisms that govern it 38, in particular in a presence of spin-forbidden transitions 12, 39, 42, 43, 44.

To infer the exciton inter-valley dynamics, we implemented polarization-resolved excitation of the FWM signal. The polarization state of the beams is adjusted by a set of and plates, by monitoring spectral interferences of with . Specifically, and are set as co-circular and thus yield a maximal contrast of the heterodyne spectral fringes. Instead, the polarization state of is rotated until the interference contrast is minimized, yielding opposite circular polarizations for (creating exciton population in the ) and (converting the population, redistributed from into , toward the FWM). To probe the exciton dynamics in the same valley (as in Fig. 3 c of the manuscript), we set the co-circular polarization of and , by maximizing the spectral interference of and . We estimate the polarization cross-talk between and around . Thus, polarization selective driving permits to probe the dynamics of the exciton transfer towards a counter-polarized valley, as depicted in Fig. 5, and thus to test the robustness of the valley pseudo-spin 35.

The measurement is presented in Fig. 5 (bottom). As in Fig. 3 c, we observe an initial fast decay of the radiatively recombining excitons within the first ps and subsequent recombination of the secondary excitons. The early dynamics for both polarization configuration is not the same during initial several ps of . Instead, for longer delays, ps, the measured FWM (generated by the secondary excitons) is the same for both driving configurations. We calculate the FWM circular polarization degree, which we define as , where denotes the spectral FWM amplitude. The resulting is plotted at the top of Fig. 5. We thus do observe a significant pseudo-spin polarisation around zero delay, , which however vanishes extremely rapidly, within around 5 ps, such that for ps. This result brings a new input (demonstration of a possibility to create valley-polarized excitons, but their rapid depolarization) towards a firm understanding of an intriguingly weak efficiency of optical pumping in MoSe monolayers (in contrast to other S-TMDs) 45. The polarization degree, relevant for valleytronics applications, could be stabilized by using magnetic fields 38 or by fabricating more involved hetero-structures 46, 47 based on S-TMDs. When analyzing Fig. 5 it is worth to note a large fraction of the FWM response within initial several ps for both polarization configurations, as compared to the subsequent dynamics. This shows that the exciton density is mainly removed radiatively within the sub-ps radiative lifetime T, as discussed before.


By performing FWM spectroscopy, we have demonstrated a giant, optical, coherent, nonlinear response of exciton transitions in MoSe monolayers. We find an increase in the FWM amplitude by two orders of magnitude, as compared to a high quality CdTe semiconductor quantum well (not shown), emitting at the similar wavelength and driven with a comparable excitation power. Taking advantage of the microscopy approach, we have performed the hyperspectral imaging, which allowed us to accurately determine the areas yielding FWM responses of either neutral, or charged excitons, or both. Electrical gating of such structures could enable spectral control of coherent responses from TMDs, providing insights into the mutual influence of TR and EX onto their dynamics. Using two-beam FWM micro-spectroscopy, we have measured the excitons’ coherence evolution, accessing dephasing times T, which turn out to be intrinsically limited by the radiative lifetime. The inhomogeneous broadening on a sub-m range is reduced by an order of magnitude with respect to measured on the entire size of the flake, although still gives rise to the pronounced photon echo. We have also shown that T spatially varies across the flake, as the T is governed by the local disorder, determining the exciton localization. Prospective experiments, exploiting two-dimensional FWM spectroscopy, will ascertain coupling mechanisms between exotic valley-exciton species (excitons, trions, biexcitons and multiexcitons) offered by S-TMDs. Our approach could be used to retrieve coherent responses of the localized exciton states in S-TMDs 48, which are expected to exhibit ultra-long dephasing 49. Finally, by exploiting polarization-resolved, three-beam FWM we assessed the exciton population dynamics, revealing the interplay between the sub-ps radiative decay T and the intra-valley relaxation T of the excitons scattered towards the dark states, revisiting the light-matter coupling in two-dimensional systems 5.

We acknowledge the support by the ERC Starting Grant PICSEN contract no. 306387 and the ERC Advanced Grant MOMB contract no. 320590. We thank G. Nogues and M. Richard for helpful remarks on the manuscript. We also acknowledge the support from Nanofab facility of the Institute Néel, CNRS UGA. Monolayers of MoSe were obtained by means of polydimethylsiloxane-based exfoliation of bulk crystals purchased from HQ Graphene.


  • Novoselov et al. 2004 Novoselov, K. S.; Geim, A. K.; Morozov, S. V.; Jiang, D.; Zhang, Y.; Dubonos, S. V.; Grigorieva, I. V.; Firsov, A. A. Science 2004, 306, 666–669.
  • Fiori et al. 2014 Fiori, G.; Bonaccorso, F.; Iannaccone, G.; T.Palacios,; Neumaier, D.; Seabaugh, A.; Banerjee, S. K.; Colombo, L. Nat. Nanotech. 2014, 10, 768–779.
  • Koppens et al. 2014 Koppens, F. H. L.; Mueller, T.; Avouris, P.; Ferrari, A. C.; Vitiello, M. S.; Polini, M. Nat. Nanotech. 2014, 10, 780–793.
  • Mak et al. 2010 Mak, K. F.; Lee, C.; Hone, J.; Shan, J.; Heinz, T. F. Phys. Rev. Lett. 2010, 105, 136805.
  • Deveaud et al. 1991 Deveaud, B.; Clérot, F.; Roy, N.; Satzke, K.; Sermage, B.; Katzer, D. S. Phys. Rev. Lett. 1991, 67, 2355.
  • Ye et al. 2014 Ye, Z.; Cao, T.; O’Brien, K.; Zhu, H.; Yin, X.; Wang, Y.; Louie, S.G.; Zhang, X.; Nature 2014, 513, 214–218.
  • Chernikov et al. 2014 Chernikov, A.; Berkelbach, T.C.; Hill, H.M.; Rigosi, A.; Li, Y.; Aslan, O.B.; Reichman, D.R.; Hybertsen, M.S.; Heinz, T.F. Phys. Rev. Lett. 2014, 113, 076802.
  • Wang et al. 2015 Wang, G.; Marie, X.; Gerber, I.; Amand, T.; Lagarde, D.; Bouet, L.; Vidal, M.; Balocchi, A.; Urbaszek, B. Phys. Rev. Lett. 2015, 114, 097403.
  • Olsen et al. 2016 Olsen, T.; Latini, S.; Rasmussen, F.; Thygesen, K. S. Phys. Rev. Lett. 2016, 116, 056401.
  • Leavitt and Little 1990 Leavitt, R. P.; Little, J. W. Phys. Rev. B 1990, 42, 11774.
  • Andreani 1995 Andreani, L. C. In Confined Electrons and Photons: New Physics and Applications; Burstein, E., Weisbuch, C., Eds.; Nato Science Series: B Physics; Plenum Press: New York, 1995; Vol. 340; pp 57–112.
  • Arora et al. 2015 Arora, A.; Koperski, M.; Nogajewski, K.; Marcus, J.; Faugeras, C.; Potemski, M. Nanoscale 2015, 7, 10421–10429.
  • Palummo et al. 2015 Palummo, M.; Bernardi, M.; Grossman, J. C. Nano Lett. 2015, 15, 2794.
  • Poellmann et al. 2015 Poellmann, C.; Steinleitner, P.; Leierseder, U.; Nagler, P.; Plechinger, G.; Porer, M.; Bratschitsch, R.; Schüller, C.; Korn, T.; Huber, R. Nat. Mater. 2015, 14, 889–893.
  • Moody et al. 2015 Moody, G.; Dass, C. K.; Hao, K.; Chen, C.-H.; Li, L.-J.; Singh, A.; Tran, K.; Clark, G.; Xu, X.; Berghaüser, G.; Malic, E.; Knorr, A.; Li, X. Nat. Comm. 2015, 6, 8315.
  • Liu et al. 2014 Liu, X.; Galfsky, T.; Sun, Z.; Xia, F.; chen Lin, E.; Lee, Y.-H.; Kéna-Cohen, S.; Menon, V. M. Nat. Phot. 2014, 9, 30–34.
  • Malard et al. 2013 Malard, L. M.; Alencar, T. V.; Barboza, A. P. M.; Mak, K. F.; de Paula, A. M. Phs. Rev. B 2013, 87, 201401(R).
  • Lagarde et al. 2014 Lagarde, D.; Bouet, L.; Marie, X.; Zhu, C. R.; Liu, B. L.; Amand, T.; Tan, P. H.; Urbaszek, B. Phys. Rev. Lett. 2014, 112, 047401.
  • Hao et al. 2015 Hao, K.; Moody, G.; Wu, F.; Dass, C. K.; Xu, L.; Chen, C.-H.; Li, M.-Y.; Li, L.-J.; MacDonald, A. H.; Li, X. Nat. Phys. doi:10.1038/nphys3674 2016,
  • Fras et al. 2015 Fras, F.; Mermillod, Q.; Nogues, G.; Hoarau, C.; Schneider, C.; Kamp, M.; Höfing, S.; Langbein, W.; Kasprzak, J. Nat. Phot. 2016, 10, 155.
  • Langbein and Patton 2005 Langbein, W.; Patton, B. Phys. Rev. Lett. 2005, 95, 017403.
  • Wang et al. 2015 Wang, G.; Palleau, E.; Amand, T.; Tongay, S.; Marie, X.; Urbaszek, B. Appl. Phys. Lett. 2015, 106, 112101.
  • Jones et al. 2016 Jones, A. M.; Yu, H.; Ghimire, N. J.; Wu, S.; Aivazian, G.; Ross, J. S.; Zhao, B.; Yan, J.; Mandrus, D. G.; Xiao, D.; Yao, W.; ; Xu, X. Nat. Nanotech. 2016, 8, 634.
  • Langbein et al. 2002 Langbein, W.; Runge, E.; Savona, V.; Zimmermann, R. Phys. Rev. Lett. 2002, 89, 157401.
  • Kasprzak et al. 2011 Kasprzak, J.; Patton, B.; Savona, V.; Langbein, W. Nat. Phot. 2011, 5, 57-63.
  • Langbein 2010 Langbein, W. Rivista del nuovo cimento 2010, 33, 255–312.
  • Naeem et al. 2015 Naeem, A.; Masia, F.; Christodoulou, S.; Moreels, I.; Borri, P.; Langbein, W. Phys. Rev. B 2015, 91, 121302(R).
  • Singh et al. 2014 Singh, A.; Moody, G.; Wu, S.; Wu, Y.; Ghimire, N. J.; Yan, J.; Mandrus, D. G.; Xu, X.; Li, X. Phys. Rev. Lett. 2014, 112, 216804.
  • Porras et al. 2002 Porras, D.; Ciuti, C.; Baumberg, J. J.; Tejedor, C. Phys. Rev. B 2002, 66, 085304.
  • Masia et al. 2011 Masia, F.; Langbein, W.; Moreels, I.; Hens, Z.; and Borri, P. Phys. Rev. B 2011, 83, 201305(R).
  • An et al. 2011 An, J. M.; Franceschetti, A.; and Zunger, A. Nano. Lett. 2007, 7, 2129.
  • Qiu et al. 2015 Qiu, D. Y.; Cao, T.; Louie, S. G. Phys. Rev. Lett. 2015, 115, 176801.
  • Xiao et al. 2012 Xiao, D.; Liu, G.-B.; Feng, W.; Xu, X.; Yao, W. Phys. Rev. Lett. 2012, 108, 196802.
  • Jones et al. 2013 Jones, A. M.; Yu, H.; Ghimire, N. J.; SanfengWu,; Aivazian, G.; Ross, J. S.; Zhao, B.; Yan, J.; Mandrus, D. G.; Xiao, D.; Yao, W.; Xu, X. Nat. Nanotech. 2013, 8, 634–638.
  • Xu et al. 2014 Xu, X.; Yao, W.; Xiao, D.; Heinz, T. F. Nat. Phys. 2014, 10, 343–350.
  • Yu and Wu 2014 Yu, T.; Wu, M. W. Phys. Rev. B 2014, 89, 205303.
  • Glazov et al. 2014 Glazov, M. M.; Amand, T.; Marie, X.; Lagarde, D.; Bouet, L.; Urbaszek, B. Phys. Rev. B 2014, 89, 201302(R).
  • Smoleński et al. 2016 Smoleński, T.; Goryca, M.; Koperski, M.; Faugeras, C.; Kazimierczuk, T.; Nogajewski, K.; Kossacki, P.; Potemski, M. Phys. Rev. X 2016, 6, 021024
  • Zhang et al. 2015 Zhang, X.-X.; You, Y.; Zhao, S. Y. F.; Heinz, T. F. Phys. Rev. Lett. 2015, 115, 257403.
  • Wang et al. 2015 Wang, G.; Robert, C.; Suslu, A.; Chen, B.; Yang, S.; Alamdari, S.; Gerber, I. C.; Amand, T.; Marie, X.; Tongay, S.; Urbaszek, B. Nat. Comm. 2015, 6, 10110.
  • Horzum et al. 2013 Horzum, S.; Sahin, H.; Cahangirov, S.; Cudazzo, P.; Rubio, A.; Serin, T.; Peeters. F.M.; Phys. Rev. B. 2013, 87, 125415.
  • Wang et al. 2015 Wang, G.; Robert, C.; Suslu, A.; Chen, B.; Yang, S.; Alamdari, S.; Gerber, I. C.; Amand, T.; Marie, X.; Tongay, S.; Urbaszek, B. Nat. Comm. 2015, 6, 10110.
  • Echeverry et al. 2016 Echeverry, J. P.; Urbaszek, B.; Amand, T.; Marie, X.; Gerber, I. C. arXiv:1601.07351 2016,
  • Withers et al. 2015 Withers, F. et al. Nano Lett. 2015, 15, 8223–8228.
  • Dery and Song 2015 Dery, H.; and Song, Y. Phys. Rev. B 2015, 92, 125431.
  • Rivera et al. 2015 Rivera, P.; Schaibley, J. R.; Jones, A. M.; Ross, J. S.; Wu, S.; Aivazian, G.; Klement, P.; Seyler, K.; Clark, G.; Ghimire, N. J.; Yan, J.; Mandrus, D. G.; Yao, W.; Xu, X. Nat. Comm. 2015, 6, 6242.
  • Rivera et al. 2016 Rivera, P.; Seyler, K. L.; Yu, H.; Schaibley, J. R.; Yan, J.; Mandrus, D. G.; Yao, W.; Xu, X. arXiv:1601.02641 2016,
  • Koperski et al. 2015 Koperski, M.; Nogajewski, K.; Arora, A.; Cherkez, V.; Mallet, P.; Veuillen, J.-Y.; Marcus, J.; Kossacki, P.; Potemski, M. Nat. Nanotech. 2015, 10, 503–506.
  • Schaibley et al. 2015 Schaibley, J. R.; Karin, T.; Yu, H.; Ross, J. S.; Rivera, P.; Jones, A. M.; Scott, M. E.; Yan, J.; Mandrus, D. G.; Yao, W.; Fu, K.-M.; Xu, X. Phys. Rev. Lett. 2015, 114.

Figures and captions:

Figure 1: Four-wave mixing spectroscopy of the MoSe monolayer. a) Pulse sequence employed in FWM experiments and related observables. b) FWM spectral interferogram obtained on the flake position displaying the exciton (EX) and trion (TR) transitions. c) FWM intensity (black) retrieved from b) via spectral interferometry. is shown with a blue line. d) Spectrally integrated FWM intensity of EX as a function of intensity, showing the driving range yielding the regime of the FWM.
Figure 2: Hyperspectral mapping of the MoSe monolayer. a) Optical contrast of the sample observed in reflectance. b)-c) Spatial imaging of the FWM amplitude, spectrally averaged over the exciton (b) response, centered at meV and the trion (c) at meV. Heterodyning at the FWM frequency with fs, TK. d) Top: A typical PL spectrum, non-resonantly excited at at meV (nm) with W average power arriving at the sample surface, displaying bright emission from EX and TR, with integrated count rate of 350kHz for each transition. Bottom: Correlation between the PL intensity of the TR (black) and its binding energy (green), as a function of its transition energy. e)-f) PL imaging of TR and EX, respectively. Excitation conditions as in d). Linear color scale, as shown by horizontal bars.
Figure 3: Exciton dynamics in the MoSe monolayer inferred with FWM microscopy. a) Time-resolved FWM amplitude as a function of , showing formation of a photon echo: a Gaussian signal with a maximum for . Logarithmic color scale given by the color bar. b) FWM amplitude as a function of retrieved from the EX at TK given by (black circles, purple squares and green triangles), respectively. The simulations yielding Tfs are given by red traces. Inset: reflectance measured at T=6 K yielding meV FWHM. c) Cartoon of the considered radiative recombination and intra-valley relaxation processes is presented in the inset. The initial dynamics of the EX population measured at (5, 45, 300) K, color coding as in b). The results yield radiative lifetime T and relaxation time (see main text). d)  dependence of the FWM amplitude when simultaneously driving the EX (blue) and TR (brown) transitions, revealing a beating with 140 fs period, and thus indicating EX-TR coherent coupling. The location of the excitation is marked with a cross in Fig. 2 b, while the corresponding micro-reflectance spectrum is given in the inset.
Figure 4: Impact of a local disorder and temperature on the exciton dephasing in the MoSe monolayer. a) Time-resolved FWM amplitude, measured on the location displaying a larger , showing a temporarily narrowed photon echo. The observed increased dephasing time with respect to Fig. 3 attributed to a localization induced increase of the radiative lifetime T. b)  as a) but measured on the area showing a broader echo and thus smaller . A weaker localization yields a shorter T=T/2 than in a). Logarithmic color scale over two orders of magnitude. c) FWM amplitude dynamics as a function of temperature measured on the area as in a). Above T=125 K the dephasing is faster than the temporal resolution of the experiment. The temperature dependent homogenous broadening (inset) indicates the phonon-induced dephasing mechanism.
Figure 5: Inter-valley population dynamics in the MoSe monolayers investigated with the polarization-resolved driving of the FWM. Bottom: FWM amplitude versus of the EX transition obtained at T=K upon co- (green circles) and cross-circular (blue squares) setting of . The modeling, as in Fig. 3 c, is displayed as red lines. Equivalent dynamics have been measured (not shown) for a counter-polarized driving . Inset: Rationale of the inter-valley dynamics probed with FWM using polarization-resolved driving: induces the FWM signal in a given valley, which originates from the population generated in the counter- or co-polarized valley, by setting counter- or co-circular , respectively. Top:  delay dependence of the measured circular polarization degree, . We measure of several tens of for initial delays , followed by a total depolarization of the FWM, generated by the secondary (see the main text) exciton population.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description