<<

arXiv:1602.04342v1 [cond-mat.mes-hall] 13 Feb 2016 cinsrnt teRb splitting inter- Rabi -matter (the large strength Extremely action [20]. h ial bandgap sizable with semiconductors the direct-zone are genides dichalco- molybdenum metal transition single-layer tundsen-based contrast and In graphene, ma- to terial, two-dimensional scale. near-perfect achievements another temperature to these room transfer the to oppor- unprecedented tunity provide crystals layer proposed. re- been have cently [19] and light [14–18], novel sources [13], po- interferometry [8– nonlinear computing of 12], optical applications in condensates Numerous lariton atomic [3–7]. ultra-cold for gases that six than about higher temperatures orders at superfluidity and (BEC) condensation Bose-Einstein demon- the polaritons strate mass, effective to small Thanks in their applications [2]. potential to and due system as this well as in [1] col- phenomena macroscopic quantum to lective interest fundamental to due attention significant attracts cav- , and ity excitons of superposition quantum a o ad omtmeauesprudt fectnpolarit exciton of superfluidity room-temperature Towards rniinmtldcacgnd mono- dichalcogenide Transition-metal polaritons, of physics decade, past the In notclmcoaiywt nebde MoS embedded an with microcavity optical an emnV Kolmakov, V. German > k iaefrteplrtngsdniy hc nbe n ode superfluidit to polariton one temperature enables room which of density, observation gas for polariton the for timate occu power transitions pumping superfluid the at that gas shown ton is It room-tempe which observed. at be parameters of monolayer, range relevant disulfide imentally molybdenum embedded an with crocavity ycnieigdie iuiednmc fectnpolarito exciton of dynamics diffusive driven considering By .INTRODUCTION I. 0eV/cm 50 h iyUiest fNwYr,NwYr,NwYr 10016 York New York, New York, New of University City The h iyUiest fNwYr,Boky,NwYr 11201 York New Brooklyn, York, New of University City The 1 hsc eatet e okCt olg fTechnology, of College City York New Department, Physics 2 eas rps ipeaayi oe htpoie usef a provides that model analytic simple a propose also We . 2 E h rdaeSho n nvriyCenter, University and School Graduate The g 1 ∼ endM Pomirchi, M. Leonid 1 − > P ∼ eV 2 0 Wadadtapn oeta strength potential trapping and and mW 600 50 uehsrcnl enrpre nRf [21]. Ref. in reported been recently has ture (MoS mi- fide disul- a molybdenum in embedded with polaritons crocavity of phe- Observation collective nomena. quantum of temperature observation room for materi- candidates these exceptional make als 23] [22, dichalcogenides e)[1 n ag xio idn energy binding exciton large ∼ and [21] meV) ia ircvt iha meddMoS embedded an op- with an in microcavity temperatures tical room can at gas achieved polariton a be of superfluidity which parameters, for relevant experimentally of range and fluctuations [24]. quantum interactions are polariton-polariton factors lead- ing the where conventional microcavities with semiconductor for experiments essential polariton temperature low-temperature not that po- the is note of We drop spreading a in polariton density. thus, a and lariton trap of a spreading in in cloud results sig- sys- that the is tem in factors fluctuations important temperature the nificant of One temper- room atures. transition-metal at po- with monolayers a cavities dichalcogenide for can in experiments gas transition actual lariton in superfluid satisfied a be for conditions 0 oee,i is it However, h olo hswr st eemn a determine to is work this of goal The . 3 1 − n oa a Kezerashvili Ya. Roman and 0 . 2 Vi ooae rniinmetal transition monolayer in eV 9 ooae ttero tempera- room the at monolayer ) y. auesprudt can superfuidity rature emn h conditions the termine si rpe polari- trapped a in rs edtrieexper- determine we si notclmi- optical an in ns priory a 2 monolayer nla fthe if unclear les- ul n in ons ,2 1, 2 2

gas of interacting particles, the condensate can be significantly depleted due to mutual particle scattering: the condensate fraction for a system of N 1 particles falls as N /N log(N)/N 1/2≫[27]. In effect, despite c ∝ T0 for an ideal gas of polaritons can be esti- mated as a few hundred K because of a small mass m, the actual BEC transition temper- ature of a trapped polariton gas lies in the helium temperature range [24]. Nevertheless, a Bose system with interac- FIG. 1. Temperature of the Kosterlitz-Thouless tions can transit to a superfluid state even (KT) phase transition, Tc, as a function of the if the condensate fraction is small or absent 4 polariton gas density n estimated from Eq. (1) et all. The famous example is superfluid He (He-II), in which the occupation of the lowest where the cutoff length r0 is set equal to the exci- energetic level is only 10%, and 90% of ton Bohr radius in molybdenum disulfide aB 1 − ≈ ∼ ∼ nm and the polariton mass is m = 3.14 10 5m , particles are “pushed” out of the condensate × 0 see text. It is seen that the KT transition at the to higher energetic levels due to strong in- teratomic interactions [28]. The relation be- room temperature T = 300 K occurs at the den- − tween the Bose condensation and superfluid- sity n 4.2 1012 m 2. ≈ × ity is delicate and is discussed in a number of publications (e.g., [26, 29, 30]). The su- monolayer. Specifically, by considering the perfluid transition in a two-dimensional sys- quasiclassical diffusive dynamics of a trapped tem with interactions is a Kosterlitz-Thouless polariton gas we determined the values of the (KT) phase transition, which is related to trapping potential strength and laser pump pairing of quantized vortices [31]. The KT power, for which the room temperature po- transition temperature for a gas of interact- lariton superfluidity can be observed. ing polaritons is given by the following equa- tion [32]

II. SUPERFLUID PHASE 2πh¯2n TRANSITION IN A POLARITON GAS Tc = 2 (1) m log log(1/nr0)

A homogenous two-dimensional Bose gas where n is the gas density and r0 is the small- can only transits to a Bose condensed state, scale cutoff length. The dependence (1) for in which the occupation of the lowest ener- polaritons in a cavity with embedded molyb- getic level is macroscopically large, at zero denum disulfide is shown in Fig. 1. It is seen temperature since the temperature fluctua- in Fig. 1 that a polariton gas can undergo tions destroy the long-range order in the sys- a superfluid transition at the room temper- tem. In contrast, an ideal trapped quantum ature T = 300 K if the gas density is high 12 −2 Bose gas undergoes the BEC transition at a enough, n > nc 4.2 10 m . At −1 ≈ × finite temperature T0 =h ¯(πkB) q6kN/m n nc, the thermal wavelength of polaritons ∼ 1/2 1/2 where k is the effective trapping potential λT = 2 πh/¯ (mkBT ) 1.3 µm is about strength, N is the total number of particles, five times larger than the≈ inter-particle dis- −1/2 and m is the particle mass [25]. At low tem- tance rs = (πn) 0.28 µm and thus, peratures T T virtually all particles oc- quantum collective effects≈ are important at ≪ 0 cupy the lowest energetic level that is, a con- densities above nc. However, the gas param- densate is formed [26]. However, for a Bose eter is small, n a2 4 10−6, that justifies c B ∼ × 3

the application of the weakly interacting gas g = 6E a2 X 4 is the effective polariton- B B| k| model (2), see the next sections. polariton interaction strength [1], with EB In what follows we identify the range of and aB to be the respective binding energy experimental parameters at which the KT and Bohr radius of the excitons, while Xk superfluid transition can be observed in a is the Hopfield coefficient that describes the trapped polariton gas in microcavities with excitonic fraction in the polariton wave func- an embedded transition-metal dichalcogenide tion. monolayer. In our consideration, we focus In our approach, the driven-diffusive dy- on microcavities with an embedded molyb- namics of the polariton gas is captured via denum disulfide monolayer that is, on a sys- a quasiclassic stochastic differential equation tem where room-temperature polaritons have for the center-of-mass coordinate r(t) of the recently been observed [21]. polariton wave packets,

III. SIMULATION METHOD dr(t)= µF (r(t), t)dt + √2DdW (t). (3) The dynamics of a normal polariton gas at the scales much larger than the polariton mean free path can be described by the driven The first term in Eq. (3) describes the ef- diffusion equation fect of the force F (r, t) on the polariton dif- fusion and the second term is the Brown- ∂n ian contribution with dW (t) to be the dif- = µ div(F n)+ D∆n Γn + Pn (2) ∂t − − ferential of a Wiener process with unit vari- where n n(r, t) is the two-dimensional den- ance. The ensemble-averaged distribution sity of the≡ polariton gas, which depends on of “tracer” particles (or quasiparticles) that the spatial coordinate r = (x, y) in the mi- obey Eq. (3) reproduces the driven-diffusive crocavity plane and time t, µ is the polariton dynamics (2) [34]. The advantage of the mobility, F F (r, t) is the external force model (3) is that, in contrast to Eq. (2), acting on polaritons,≡ D is the diffusivity, Γ is the stochastic equation describes the fluctu- ations of the particle density. The model (3) the width of the polariton resonance and Pn is the source (see Ref. [1] for extensive review). is off lattice and is computationally effective In what follows, we consider the dynamics for systems with a large number of particles. of the polariton gas in an in-plane trapping The stochastic model (3) was applied for the r 1 2 description of diffusive properties of quantum potential Utrap( ) = 2 kr , where k is the effective strength of the trap and r = r . gases (see Ref. [1]) as well as for nanoparticle The main mechanisms of the diffusive collec-| | diffusion in fluids [34]. tive motion of polaritons include their mu- We numerically integrated Eq. (3) by us- tual interactions, scattering on phonons and ing the Euler method. The force acting on scettring on excitonic and photonic defects polaritons was set equal to F = U where −∇ in the microcavity [33]. (We consider an un- U = Utrap + Upp. The finite-size polariton doped cavity.) In our simulations, we take source at the pumping laser light spot was into account the density-dependent polari- modeled by adding particles at each numer- ton interactions in the mean-field approxi- ical time step ∆t with the Gaussian spatial mation where the effective repulsive interac- probability distribution centered at the spot tion potential is Upp(r, t) = gn(r, t), which center. The finite lifetime of polaritons was describes exactly the same interactions that taken into account by randomly eliminating lead to the blue shift of the polariton excita- particles from the ensemble at each time step tion spectrum in a superfluid state [1]. Here, with the probability w = Γ∆t. 4

IV. MODEL PARAMETERS TABLE I. Simulations parameters for a cavity with embedded MoS layer In our studies, we use the parameters rel- 2 evant for the A-exciton polaritons in a cav- Exciton mass mex 0.78m0 −5 ity with an embedded molybdenum disulfide mass mph 1.57 10 m0 × −5 monolayer [21]. In our simulations, we focus Polariton mass m 3.14 10 m0 × −13 on the polariton dynamics at the room tem- Polariton momentum τpol 9.60 10 s perature T = 300 K. relaxation time × The effective polariton mass was calcu- Polariton diffusion D 139 m2/s lated as coefficient Polariton interaction g 1.44 10−3 −1 −1 × 2 m =2/(mex + mph ). (4) strength meVµm Polariton mobility µ 3.36 1022 × The exciton mass in MoS2 was taken equal to s/kg mex = me+mh =0.78m0 where me =0.35m0 Dielectric constant ǫ 4.3 (pol) and mh = 0.43m0 are the and hole Polariton life time τlife 84 ps mass respectively (with m0 to be the free Trapping potential k 15 300 electron mass) [35, 36]. The effective photon strength eV/cm− 2 mass is mph = π√ǫh/cL¯ C , where the cavity Efficiency of polariton q 10−3 length LC was found from the condition of generation the resonance of the photonic and excitonic Numerical unit of ∆x 1 µm states, Eph = πhc/¯ √ǫLC , where Eph is the length energy of the cavity photon resonance. The Numerical time step ∆t = 272 fs − latter was taken equal to Eph =1.87 eV [21]. m∆x2¯h 1 The dielectric constant of MoS2 is ǫ = 4.3 [35]. The exciton diffusion coefficient was taken The respective momentum relaxation times 2 for excitons and polaritons was estimated as equal to Dex = 14 cm /s [37]. This value is τ = D m /k T and τ = X −4τ . comparable with the exciton diffusion coef- ex ex ex B pol | | ex ficient in semiconductor quantum wells [1]. The polariton resonance width was esti- The polariton diffusivity was calculated as mated as [33, 38] 2 2 Γ= X Γex + C Γph (7) −4 mex | | | | D = X Dex. (5) | | m where C 2 =1 X 2 is the photon Hopfield | | −| | In what follows, we consider zero detun- coefficient. According to Ref. [39], the ra- (ex) −1 ing between the photonic and excitonic res- diative exciton lifetime τlife Γex in MoS2 onances thus, we set X = 1/√2. However, (ex≡) monolayer increases from τlife = 0.23 ps at the detuning in the experiments with polari- T 0 to τ (ex) 4 ps at T = 10 K and then, tons in a cavity with MoS monolayer can be life 2 to→ 270 ps at the≈ room temperature. Even varied in a wide range up to 50 meV [21]. larger radiative exciton life time of 850 ps was Thus, the diffusive properties∼ of polaritons reported for MoS monolayer crystals at the can be controlled by adjusting the detuning 2 room temperature [40]. On the other hand, value. large broadening of the excitonic resonance The polariton mobility was calculated as h¯Γ = 30 meV (with the respective time follows [33, 38], ex scale of 22 fs) was reported in Ref. [21]. The τ respective short-time processes are probably µ = pol . (6) m related to non-radiative trapping of excitons 5

at impurities or hot carrier thermalization, in a cavity with MoS2 is not known; we suggest agreement with Ref. [39]. In our simulation, that the value is close to that obtained for (ex) the quantum well structures. Thus, in our we set τlife = 270 ps. The photon resonance −1 width was taken equal to Γph = (50ps) model, the polariton injection rate, Pn, was [41]. The polariton life time was estimated calculated as −1 (pol) from Eq. (7) as Γ = τlife = 84 ps, which qP is in agreement with the observations for Pn = , (9) E semiconductor-based microcavities [41]. ph In the simulations, we varied the trap- where P is the pumping laser power. In ping potential strength from k = 15 eV/cm2 the simulations, we consider a Gaussian laser to 300 eV/cm2 that lies in the experimen- spot profile with full width at half maximum tally relevant parameter range [24, 42]. The (FWHM) of 30 µm. trapping potential can be formed by apply- The simulations parameters for T = 300 ing inhomogeneous stress [42], static elec- K are summarized in Table I. tric or magnetic field, or laser radiation (see Ref. [43] and references therein). The po- lariton interaction constant was calculated as g = 6E a2 X 4 = 1.44 10−3 meV µm2 V. RESULTS AND DISCUSSION B B| | × given the exciton binding energy EB = 0.96 First, we explore the case where the exci- eV and the exciton Bohr radius aB 1 nm [21, 40]. It is worth noting that,≈ despite tation laser spot is positioned at the center of the trap. The relaxation of the total number the exciton binding energy EB in transition metal dichalcogenides is significantly larger of polaritons in the system to its steady-state than that in quantum wells, the polariton in- value is shown in Fig. 2. It is seen that the system reaches a steady state 400 ps after teraction strength is of the same order in the ∼ both cases. the excitation is switched on. The polariton density at a point r at the moment t was calculated from the positions

5 of polariton quasiparticles ri(t) (where the 1.5x10 index i labels the quasiparticles) as follows

5 r r r 1.0x10 n( , t)= X f( i(t) ) (8) P = 550 mW | − | 2 i k = 50 eV/cm

4

5.0x10 with the weight function f(r) = 2 −1 2 2 (2πrav) exp( r /2rav) and the averag- Totalnumber of particles, N − 0.0 ing length rav = 10 µm. 0 1000 2000 3000 Qualitative analysis of experimental data t, ps in Refs. [15, 18, 42, 44, 45] for quantum well FIG. 2. Relaxation to a steady state of a trapped structures shows that the efficiency of the ex- polariton cloud after the laser pump is turned on citon generation by a continuous wave (cw) at the moment t = 0. The data points show the light source is q 0.1%; in other words, only total number of polaritons as a function of time 1 of 103 incident∼ photons creates a bound ex- t citon, which is, in its turn, forms a polariton. . The excitation power of the laser pump is P = 550 mW, the trapping potential strength is A single-layer transition metal dichalcogenide 2 monolayer can capture up to 10% of incident k = 50 eV/cm . The fluctuations of the total number of polaritons in the steady state t 400 photons [46, 47]. However, the efficiency of ≥ the polariton generation for cw pumping in ps are seen in the figure. 6

The steady-state spatial distribution of polaritons in the trap for the trapping po- tential strength k = 50 eV/cm2 is shown in Fig. 3. It is seen that the characteristic po- lariton cloud size reaches 200 µm that is, about an order larger than∼ the excitation spot size 30 µm. It is also much larger than the Thomas-Fermi size of a cloud at zero temper- ature aTF = q2nmaxg/k 1 10 µm for the maximum density at the≈ center− n n max ≈ c

[24]. The main reason for cloud spreading (b) -2

4 is thermal diffusion of polaritons towards the m edges of the cloud from the excitation area. 12

P=550 mW

3 Thus, at room temperatures, thermal spread- 2 k = 50 eV/cm ing is more pronounced than the effect of the interparticle interactions. 2

Figure 4 shows the maximum density of 1 polaritons, nmax, reached at the center of the Polaritondensity n, 10 cloud as a function of the excitation power 0 P . The data are shown for the trapping po- -400 -200 0 200 400 tential strength k = 50 eV/cm2 and tempera- x, m ture T = 300 K. It is seen that the maximum FIG. 3. (Color online) Spatial distribution of density n nearly linearly grows with the max trapped polaritons in a cavity with MoS2 mono- increase of the power P . For the excitation layer at T = 300 K. The power of the laser pump power P = P 600 mW, the maximum c ≈ and the trapping potential strength are marked density of polaritons reaches the critical den- in the figure. (a) The two-dimensional polariton sity n 4.2 1012 m−2. At the pumping c density profile n(r) averaged over the time do- power P≈ >P ×the central part of a polariton c main t = 400 2700 ps, in which the system has system at room temperatures should demon- − reached the steady state. The color bar shows strate superfluid properties previously known the polariton density in the units of 1012 m−2. for a low-temperature polariton superfluid in a semiconductor microcavity (see review [1]). (b) A one-dimensional cut of the graph (a) made The wings of the spatial distribution always along the line y = 0. Points show the results of the simulations, vertical bars show the standard have a density lower than nc thus, the super- fluid “lake” at the cloud center is surrounded deviation about the mean for the density fluctu- by a low-density gas of polaritons in a normal ations. In the simulations, the full width at half state. maximum (FWHM) of the laser spot is 30 µm; the center of the laser spot is positioned at the We emphasize that, after a part of the po- center of the parabolic trap r = (0, 0). lariton system transits to a superfluid state, quantum correlations in the polariton dy- namics should be taken into account. How- of a small superfluid fraction will not mod- ever, if the superfluid density is small com- ify the value of the maximum density in the pared to the total density of the system, cloud and thus, the criterion of the super- nmax nc < nc, the superfluid effects on the fluidity formation will not be changed. (We over-all− large-scale density distribution in a note that a few-percent density jump can oc- steady state can be disregarded to the first cur during the superfluid transition, in a sim- approximation. In particular, the presence ilarity with the normal-to-superfluid transi- 7

tion in He-II [48]; in our case however, this jump is less than or comparable with a typical 8 750 mW

value of the thermal fluctuations of the den- 600 mW sity 4% estimated from Fig. 3 and there- 300 mW 2 fore,∼ it can be disregarded.) We also disre- 6 cm

gard the dissipative dynamics at a very close 12 vicinity of the superfluid transition where the , 10

4

superfluid component can appear as a fluc- max tuation mode [49]. Our goal is to pinpoint n a domain in the parameter space where the room-temperature superfluidity of polaritons 2

0 100 200 300

can be experimentally achieved. 2

k, eV/cm

FIG. 5. Dependence of the maximum density 8 of the polariton cloud on the trapping potential

2 strength k. Points show the results of the sim- k=50 eV/cm -2 6 ulations, lines show the theoretical dependence m

12 Eq. (11) for the respective parameters. The val- ues for the laser pump power are listed in the 4 , 10 ,

max figure. n

2

P

c sults of the numerical simulations are sum-

0 marized in Fig. 5. It is seen that the density 0 250 500 750 1000 grows with the rise of the trapping potential P, mW strength. From the obtained results it fol- lows that at experimentally realistic pump- FIG. 4. Increase in the maximum density of po- ing power P 600 mW, the cloud can tran- laritons, n , at the center of the cloud r = max sit to a superfluid∼ state at T = 300 K for (0, 0) with the rise of the power of the laser pump the trapping strength k 50 eV/cm2. For a P . Points show the results of the simulations for ≥ 2 larger power P = 750 mW, the polariton den- k = 50 eV/cm , the line segments connecting the sity required for the room-temperature su- points are shown to guide the eye. The dashed perfluidity is reached at even smaller trap- horizontal line corresponds to the critical density ping strengths k 15 eV/cm2. On the other n = 4.2 1012 m−2, at which a gas of interacting ≥ c × hand, for P = 300 mW, the conditions for the polaritons in a cavity with MoS2 undergoes the room temperature polariton superfluidity are Kosterlitz-Thouless superfluid transition at the not satisfied even for large trapping strengths room temperature. A polaritons gas is super- k 300 eV/cm2. ∼ fluid at the central part of the cloud above the To further investigate the trapped polari- line (P > P 600 mW) and is in the normal c ≈ ton gas dynamics, we also studied the effect state below the line (P < Pc). of the spatial offset in the excitation spot po- sition compared to the center of the parabolic To better understand the optimal condi- trap. For this purpose, we numerically stud- tions for observations of room-temperature ied the polariton spatial distribution when polariton superfluidity, we studied the de- the center of the laser spot is shifted along pendence of the maximum density nmax on the x axis on the distance x0 = 200 µm the trapping potential strength. The re- that is, larger than the excitation spot size 8

of 30 µm. The obtained numerical results strengths, the limiting factor for the cloud (pol) are shown in Fig. 6. Fig. 6a shows the two- size is a finite polariton life time τlife . In this dimensional polariton density distribution for limit, the cloud size can be estimated as a = k = 150 eV/cm2, and Fig. 6b shows its cross- (pol) 1/2 (Dτlife ) . In the opposite case of strong section made along the x axis at y = 0 trapping potentials, the cloud size is limited 2 1 2 for k = 150 eV/cm (filled squares) and for by thermal fluctuations, kBT ka , and 2 2 k = 50 eV/cm (unfilled triangles). It is seen ∼ 1/2 thus, it can be estimated as a = (2kBT/k) . that the polaritons diffuse towards the cen- In the intermediate case, the cloud size can r ter of the trap = 0. It is also seen that be obtained from the interpolation equation the polariton density at the center of the trap increases with the rise of the trapping (pol) 1/2 (Dτlife ) potential strength k. However, for moder- a = (10) √1+ γ ate strengths k 150 eV/cm2, the polariton r ≤ density at = 0 remains small compared to where γ = kDτ (pol)/2k T is the dimension- the maximum density at the excitation spot life B r less trapping potential strength. We note area. Specifically, the density at = 0 is that for the parameters corresponding to an 19% of the maximum density in the cloud embedded MoS monolayer at T = 300 K for≈ k = 150 eV/cm2 and it is only 10% for 2 2 (Table I), one obtains γ = 1 for k 440 k = 50 eV/cm . eV/cm2. By approximating the steady-state≈ The geometry with a shifted excitation cloud shape as an axially-symmetric pro- 6 spot might provide an additional opportu- file n(r) = nmax/(1 + (r/a) ) and the to- nity for the experimental detection of room- tal number of polaritons in the system as (pol) temperature polariton superfluidity. One can N Pnτlife , one has the following expres- expect a larger occupation of the lowest en- sion≈ for the maximum polariton density, ergetic state near the trap center r =0 for a (pol) (pol) superfluid gas compared to the case of a nor- 3√3 Pnτ Pnτ n = life 0.26 life , (11) mal gas (these simulations). In other words, max 2π2 a2 ≈ a2 the superfluid polariton gas would “easier” flow towards the bottom of the trap that, in in which the parameter a is given in Eq. (10) its turn, should lead to higher densities n(0). and the source Pn is estimated in Eq. (9). Thus, observation of “anomalously” high po- The model Eq. (11) is in good agreement with lariton population n(0) with the increase of the results of the simulations, as is seen in the trapping strength k or the laser power P Fig. 5. in experiments with a shifted laser spot can be a simple experimental test of the normal- VI. CONCLUSIONS to-superfluid transition in the polariton gas. Finally, we provide an estimate of the In this paper, we studied the polariton gas cloud densities that facilitates the search of dynamics in a parabolic trapping potential parameters suitable for the polariton super- by using the quasiclassical stochastic equa- fluidity observation. As it follows from the tion for the polarion wave packets diffusion. above consideration, a key parameter that We showed that in a cavity with an embed- determines the feasibility of polariton super- ded molybdenum disulfide monolayer, room- fluidity at room temperatures is the maxi- temperature superfluidity in a trapped po- mum polariton density at the center of the lariton cloud can be observed for experimen- cloud nmax. To estimate the maximum den- tally realistic trapping strengths k 50 300 2 ∼ − sity nmax the following model can be uti- eV/cm and laser excitation power P 600 lized. In the case of small trapping potential mW. The experiments with a shifted≥ laser 9

shows that a trapped polariton gas is super- fluid at the central part of the polariton cloud if the laser pumping power exceeds the criti- cal value 3.8Dn E P = c ph . (12) c 1+ γ q In our studies, we are inspired by the re- cent observation of room temperature polari-

5 tons [21] as well as by the advances of low- (b) temperature polariton physics [1, 2]. Su- 4 perfluidity of excitons in transition metal -2 dichalcogenide-based systems has recently

3 m

12 been predicted for temperatures beyond the

2 helium temperature range [50]. In addi- n,10 tion to the fundamental importance of high-

1 temperature quantum effects, polariton su- perfluidity might be a key to new opto-

0

-200 0 200 400 electronic technologies including those men-

x, m tioned in Introduction. Room-temperature polariton superfluidity would open a new FIG. 6. (Color online) Polariton diffusion to- route towards the design of photonic inte- wards the center of the trap in the case where the grated circuits suitable for the use in opti- excitation spot is shifted from the trap center. cal information processing systems with low The pumping power is P = 600 mW. The offset energy consumption. Our hope is that the for the excitation spot center is x0 = 200 µm in experimental search for room-temperature the positive direction of the x axis. (a) Two- quantum fundamental effects in transition dimensional polariton density distribution for metal dichalcogenide-based heterostructures k = 150 eV/cm2. (b) The density distribution would facilitate the progress in this direction. calculated along the line y = 0 for k = 50 (un- filled triangles) and 150 eV/cm2 (filled squares). For the reference, the density distribution for VII. FUNDING INFORMATION k = 50 eV/cm2 with zero excitation spot offset, National Science Foundation (NSF) (Sup- x0 = 0, is shown by the dashed curve. Fluctu- ations of the polariton density ( 4%) are not plement to 1345219); Professional Staff ∼ shown. Congress – City University of New York (PSC-CUNY) (68090-00 46). The authors thank the Center for Theoret- ical Physics of the New York City College of spot would provide a simple qualitative test Technology for providing computational re- for superfluidity. The theoretical estimate sources.

[1] I. Carusotto and C. Ciuti, “Quantum flu- N. Peyghambarian, eds., Nonlinear Pho- ids of light,” Rev. Mod. Phys. 85, 299–366 tonics (Springer, London, 2011). (2013). [3] H. Deng, G. Weihs, C. Santori, J. Bloch, [2] H. M. Gibbs, G. Khitrova, and and Y. Yamamoto, “Condensation of semi- 10

conductor microcavity exciton polaritons,” Nat. Commun. 4, 1778 (2013). Science 298, 199–202 (2002). [12] O. L. Berman, R. Ya. Kezerashvili, and [4] J. Kasprzak, M. Richard, S. Kundermann, G. V. Kolmakov, “Harnessing the polari- A. Baas, P. Jeambrun, J. M. J. Keeling, ton drag effect to design an electrically con- F. M. Marchetti, M. H. Szymanska, R. An- trolled optical switch,” ACS Nano 8, 10437– dre, J. L. Staehli, V. Savona, P. B. Lit- 10447 (2014). tlewood, B. Deveaud, and L. S. Dang, [13] C. Sturm, D. Tanese, H. S. Nguyen, “Bose-Einstein condensation of exciton po- H. Flayac, E. Galopin, A. Lemaitre, laritons,” Nature 443, 409–414 (2006). I. Sagnes, D. Solnyshkov, A. Amo, [5] R. Balili, V. Hartwell, D. Snoke, L. Pfeiffer, G. Malpuech, and J. Bloch, “All-optical and K. West, “Bose-Einstein condensation phase modulation in a cavity-polariton of microcavity polaritons in a trap,” Science Mach–Zehnder interferometer,” Nat. Com- 316, 1007–1010 (2007). mun. 5, 3278 (2014). [6] S. Utsunomiya, L. Tian, G. Roumpos, [14] H. Deng, G. Weihs, D. Snoke, J. Bloch, and C. W. Lai, N. Kumada, T. Fujisawa, Y. Yamamoto, “Polariton lasing vs. pho- M. Kuwata-Gonokami, A. L¨offler, ton lasing in a semiconductor microcavity,” S. H¨ofling, A. Forchel, and Y. Yamamoto, Proc. Natl. Acad. Sci. U.S.A. 100, 15318– “Observation of Bogoliubov excitations 15323 (2003). in exciton-polariton condensates,” Nat. [15] G. Malpuech, A. Kavokin, A. Di Carlo, Physics 4, 700–705 (2008). and J. J. Baumberg, “Polariton lasing by [7] A. Amo, D. Sanvitto, F. P. Laussy, exciton-electron scattering in semiconduc- D. Ballarini, E. del Valle, M. D. Martin, tor microcavities,” Phys. Rev. B 65, 153310 A. Lemaitre, J. Bloch, D. N. Krizhanovskii, (2002). M. S. Skolnick, C. Tejedor, and L. Vina, [16] M. H. Szymanska and P. B. Littlewood, “Collective fluid dynamics of a polariton “The crossover between lasing and polari- condensate in a semiconductor microcav- ton condensation in optical microcavities,” ity,” Nature 457, 291–293 (2009). Solid State Communications 124, 103–107 [8] D. Bajoni, E. Semenova, A. Lemaitre, (2002). S. Bouchoule, E. Wertz, P. Senellart, [17] D. Bajoni, P. Senellart, A. Lemaˆıtre, and S. Barbay, R. Kuszelewicz, and J. Bloch, J. Bloch, “Photon lasing in GaAs microcav- “Optical bistability in a GaAs-based polari- ity: Similarities with a polariton conden- ton ,” Phys. Rev. Lett. 101, 266402 sate,” Phys. Rev. B 76, 201305 (2007). (2008). [18] S. Christopoulos, G. B. H. von H¨ogersthal, [9] V. M. Menon, L. I. Deych, and A. A. A. J. D. Grundy, P. G. Lagoudakis, A. V. Lisyansky, “Nonlinear : Towards po- Kavokin, J. J. Baumberg, G. Christmann, laritonic logic circuits,” Nat. Photonics 4, R. Butt´e, E. Feltin, J.-F. Carlin, and 345–346 (2010). N. Grandjean, “Room-temperature polari- [10] T. Gao, P. S. Eldridge, T. C. H. Liew, ton lasing in semiconductor microcavities,” S. I. Tsintzos, G. Stavrinidis, G. Deligeor- Phys. Rev. Lett. 98, 126405 (2007). gis, Z. Hatzopoulos, and P. G. Savvidis, [19] J. Restrepo, C. Ciuti, and I. Favero, “Polariton condensate switch,” “Single-polariton optomechanics,” Phys. Phys. Rev. B 85, 235102 (2012). Rev. Lett. 112, 013601 (2014). [11] D. Ballarini, M. De Giorgi, E. Cancel- [20] S. Das, A. Prakash, R. Salazar, and J. Ap- lieri, R. Houdr´e, E. Giacobino, R. Cin- penzeller, “Toward low-power electronics: golani, A. Bramati, G. Gigli, and D. San- tunneling phenomena in transition metal vitto, “All-optical polariton transistor,” dichalcogenides,” ACS Nano 8, 1681–1689 11

(2014). two-dimensional systems,” J. Phys. C: Solid [21] X. Liu, T. Galfsky, Z. Sun, F. Xia, E.- St. Phys. 6, 1181–1203 (1973). C. Lin, Y.-H. Lee, S. K´ena-Cohen, and [32] Yu. E. Lozovik, A. G. Semenov, and V. M. Menon, “Strong light–matter cou- M. Willander, “Kosterlitz-Thouless tran- pling in two-dimensional atomic crystals,” sition in microcavity polariton system,” Nat. Photon. 9, 30–34 (2015). JETP Lett. 84, 146–150 (2006). [22] D. Y. Qiu, F. H. da Jornada, and S. G. [33] G. Bley, U. Neukirch, and J. Gutowski, Louie, “Optical spectrum of MoS2: Many- “Exciton polariton diffusion in ZnSe lay- body effects and diversity of exciton states,” ers,” J. Crystal Growth 184, 1095–1099 Phys. Rev. Lett. 111, 216805 (2013). (1998). [23] A. Chernikov, T. C. Berkelbach, H. M. Hill, [34] P. Szymczak and A. J. C. Ladd, “Bound- A. Rigosi, Y. Li, O. B. Aslan, D. R. Reich- ary conditions for stochastic solutions of the man, M. S. Hybertsen, and T. F. Heinz, convection-diffusion equation,” Phys. Rev. “Exciton binding energy and nonhydro- E 68, 036704 (2003). genic Rydberg series in monolayer WS2,” [35] A. Ramasubramaniam, “Large excitonic ef- Phys. Rev. Lett. 113, 076802 (2014). fects in monolayers of molybdenum and [24] O. L. Berman, Yu. E. Lozovik, and D. W. tungsten dichalcogenides,” Phys. Rev. B Snoke, “Theory of Bose-Einstein condensa- 86, 115409 (2012). tion and superfluidity of two-dimensional [36] T. Cheiwchanchamnangij and W. R. L. polaritons in an in-plane harmonic poten- Lambrecht, “Quasiparticle band structure tial,” Phys. Rev. B 77, 155317 (2008). calculation of monolayer, bilayer, and bulk [25] V. Bagnato and D. Kleppner, “Bose- MoS2,” Phys. Rev. B 85, 205302 (2012). Einstein condensation in low-dimensional [37] F. Ceballos, M. Z. Bellus, H.-Y. Chiu, and traps,” Phys. Rev. A 44, 7439–7441 (1991). H. Zhao, “Probing charge transfer exci- [26] F. Dalfovo, S. Giorgini, L. P. Pitaevskii, tons in a MoSe2–WS2 van der Waals het- and S. Stringari, “Theory of Bose-Einstein erostructure,” Nanoscale 7, 17523–17528 condensation in trapped gases,” Rev. Mod. (2015). Phys. 71, 463–512 (1999). [38] O. L. Berman, R. Ya. Kezerashvili, and [27] M. Holzmann, G. Baym, J.-P. Blaizot, and Yu. E. Lozovik, “Drag effects in a system F. Lalo¨e, “Superfluid transition of homo- of and microcavity polaritons,” geneous and trapped two-dimensional Bose Phys. Rev. B 82, 125307 (2010). gases,” Proc. Natl. Acad. Sci. U.S.A. 104, [39] M. Palummo, M. Bernardi, and J. C. 1476–1481 (2007). Grossman, “Exciton radiative lifetimes in [28] K. R. Atkins, Liquid Helium (Cambridge two-dimensional transition metal dichalco- University Press, Cambridge, 1959). genides,” Nano Letters 15, 2794–2800 [29] VI Yukalov, “Basics of Bose-Einstein con- (2015). densation,” Physics of Particles and Nuclei [40] H. Shi, R. Yan, S. Bertolazzi, J. Brivio, 42, 460–513 (2011). B. Gao, A. Kis, D. Jena, H. G. Xing, and [30] O. L. Berman, R. Ya. Kezerashvili, and L. Huang, “Exciton dynamics in suspended G. V. Kolmakov, “Introduction to nonlinear monolayer and few-layer MoS2 2D crys- phenomena in superfluid liquids and Bose- tals,” ACS Nano 7, 1072–1080 (2013). Einstein condensates: Helium, semiconduc- [41] B. Nelsen, G. Liu, M. Steger, D. W. Snoke, tors and graphene,” Contemporary Phys. R. Balili, K. West, and L. Pfeiffer, “Dis- 52, 319–340 (2011). sipationless flow and sharp threshold of a [31] J. M. Kosterlitz and D. J. Thouless, “Order- polariton condensate with long lifetime,” ing, metastability and phase transitions in Phys. Rev. X 3, 041015 (2013). 12

[42] V. Negoita, D. W. Snoke, and K. Eberl, Lett. 105, 136805 (2010). “Harmonic-potential traps for indirect exci- [47] A. Splendiani, L. Sun, Y. Zhang, T. Li, tons in coupled quantum wells,” Phys. Rev. J. Kim, C.-Y. Chim, G. Galli, and B 60, 2661–2669 (1999). F. Wang, “Emerging in [43] A. A. High, J. R. Leonard, M. Remeika, monolayer MoS2,” Nano Letters 10, 1271– L. V. Butov, M. Hanson, and A. C. Gos- 1275 (2010). sard, “Condensation of excitons in a trap,” [48] I. M. Khalatnikov, An Introduction to the Nano Lett. 12, 2605–2609 (2012). Theory of Superfluidity (Benjamin, New [44] J. Feldmann, G. Peter, E. O. G¨obel, York, 1965). P. Dawson, K. Moore, C. Foxon, and [49] I. M. Khalatnikov, V. L. Pokrovskii, and R. J. Elliott, “Linewidth dependence of ra- D. M. Semiz, “High-frequency hydrody- diative exciton lifetimes in quantum wells,” namics near the λ point,” J. Low Temp. Phys. Rev. Lett. 59, 2337–2340 (1987). Phys. 6, 305–321 (1972). [45] L. V. Butov, A. C. Gossard, and D. S. [50] O. L. Berman and R. Ya. Kezerashvili, Chemla, “Macroscopically ordered state in “High-temperature superfuidity of the two- an exciton system,” Nature 418, 751–754 component Bose gas in a TMDC bi- (2002). layer,” arXiv:1601.07981 [cond-mat.mes- [46] K. F. Mak, C. Lee, J. Hone, J. Shan, and hall] (2016). T. F. Heinz, “Atomically thin MoS2: A new direct-gap semiconductor,” Phys. Rev.