Polariton condensation in photonic molecules
Marta Galbiati, Lydie Ferrier, Dmitry D. Solnyshkov, Dimitrii Tanese, Esther Wertz, Alberto Amo, Pascale Senellart, Isabelle Sagnes, Aristide Lemaitre, Elisabeth Galopin, Guillaume Malpuech, Jacqueline Bloch
aa r X i v : . [ c ond - m a t . o t h e r] O c t Polariton Condensation in photonic molecules
Marta Galbiati, Lydie Ferrier, Dmitry D. Solnyshkov, Dimitrii Tanese, Esther Wertz, Alberto Amo, PascaleSenellart, Isabelle Sagnes, Aristide Lemaˆıtre, Elisabeth Galopin, Guillaume Malpuech, and Jacqueline Bloch CNRS-Laboratoire de Photonique et Nanostructures, Route de Nozay, 91460 Marcoussis, France LASMEA, CNRS, Clermont University, University Blaise Pascal,24 avenue des Landais, 63177 Aubi`ere cedex, France (Dated: October 10, 2018)We report on polariton condensation in photonic molecules formed by two coupled micropillars.We show that the condensation process is strongly affected by the interaction with the cloud ofuncondensed excitons. Depending on the spatial position of these excitons within the molecule,condensation can be triggered on both binding and anti-binding polariton states of the molecule, ona metastable state or a total transfer of the condensate into one of the micropillars can be obtained.Our results highlight the crucial role played by relaxation kinetics in the condensation process.
PACS numbers: 71.36.+c, 67.85.Hj, 78.67.Pt, 78.55.Cr
Most of the experimental studies in atomic Bose con-densates have explored conditions of thermodynamicequilibrium since typical condensate lifetimes are muchlonger than interaction times. Recent theoretical pro-posals have shown that out of equilibrium bosonic sys-tems present qualitatively new behaviors [1]. One pro-posed way to reach this regime is the use of photonic sys-tems with effective photon-photon interactions and dissi-pation provided by inherent optical losses [2]. Localizedto delocalized phase transitions [3, 4], highly entangledstates [5], or fermionisation effects in a ring of coupledsites [6] are predicted in such systems.Microcavity polaritons are a model system for the in-vestigation of the physics of driven-dissipative boson con-densates [7–13]. They are the quasi-particules arisingfrom the strong coupling between excitons confined inquantum wells and the optical mode of a microcavity.Because of their light-matter nature, polaritons presentpeculiar properties: they interact efficiently with theirenvironment through their excitonic part [14, 15] whiletheir photonic part enables efficient coupling with the freespace optical modes. Polariton condensates can be gen-erated in zero dimensional micropillars [11] or in arrays ofpillars with fully controlled coupling [16, 17]. In this con-figuration, the non-equilibrium nature of polariton con-densates should allow the realization of metastable col-lective states, such as the self-trapped states in a bosonicJosephson junction [18–20].In the present paper we investigate polariton conden-sation in photonic molecules obtained by coupling twomicropillars. We demonstrate that polariton interactionsstrongly affect the way condensation occurs in such cou-pled system, not only modifying the wavefunction of thepolariton condensate, but also the relaxation dynamics.This effect, specific to an out-of-equilibrium bosonic sys-tem, is illustrated by considering different positions of thenon resonant excitation within the molecule. When theexcitation spot is placed at the center of the molecule,polariton condensation is observed on both binding and anti-binding states. Interactions induce strong changesin the condensate wavefunction, the most important onebeing the change in its spatial anisotropy.When the excitation spot is positioned on one of thetwo coupled micropillars, condensation occurs in a verydifferent way. As the excitation power is increased, thepolariton condensate is first created in a metastable statelocalized in the excited pillar by an effect similar to selftrapping[18]. Further increasing the pumping power,the condensation tends to progressively evolve from akinetic regime to a regime closer to thermodynamicalequilibrium: massive occupation is observed on the low-est energy state of the system, mainly localized in thenon-excited pillar. These results are simulated usinga relaxation model including semi-classical Boltzmannequations and nonlinear Schr¨odinger equation in a self-consistent way. The key features of our experiments canonly be reproduced if one properly includes the changesof the relaxation rates due to interaction-induced modi-fications of the overlaps between the different polaritonstates.The microcavity sample, described in Ref. [11], consistsin a λ/ d of the micropillars ranges from 3 to 4 µm and their cen-ter to center distance d CC is varied from 2.3 to 3 . µm ,corresponding to an expected coupling constant between0 . meV . Note that d CC is always kept smallerthan twice the radius of the pillars, ensuring the directcoupling of the polariton modes of the two micropillars.A scanning electron micrography of an array of such pho-tonic molecules is shown in Fig.1.a. Microphotolumines-cence experiments are performed on single molecules us-ing a single mode cw Ti:Saph laser focused onto a 2 µ mdiameter spot with a microscope objective. The sampleis maintained at 10 K and the excitation laser energy istuned typically 100 meV above the lower polariton reso- FIG. 1: (Color online) (a) Emission spectrum measured ona single molecule at low excitation power ( d = 4 µm and d CC = 3 . µm ); A and B indicate the binding and anti-binding states. Inset : scanning electron micrograph of anarray of pillars and molecules. b) Emission spectra measuredon a 4 µ m round micropillar (black line) and on photonicmolecules with d = 4 µm and various values of d CC . c) left(resp. right) column : measured (resp. calculated) emissionpattern of the four lowest energy modes of a photonic moleculewith d = 4 µm and d CC = 3 . µm ; d) Schematic of thehybridization of the individual pillar modes within a photonicmolecule. nance, thus providing non-resonant optical excitation ofthe system. The emission is collected through the sameobjective and imaged on the entrance slit of a monochro-mator. The spectrally dispersed emission is detected witha nitrogen-cooled CCD camera. We define the detuning δ = E c − E x as the energy difference between the lowestenergy photonic mode and the exciton resonance.The polariton modes in these photonic molecules areinvestigated by photoluminescence measurements at lowexcitation power. An example of an emission spec-trum measured on a single molecule is presented inFig. 1.a. Discrete emission peaks are observed corre-sponding to polariton quantum states fully confined inthe microstructure. The broader line at higher energyis due to emission of the excitonic reservoir. The twolowest energy modes (labeled B and A) are attributed tothe binding and anti-binding states arising from the hy-bridization of the lowest energy mode of each micropillar.The splitting between these states is proportional to thecoupling between the two micropillars. It can be contin- FIG. 2: (Color online)
Excitation at the center ofthe molecule : a)-c)Spectrally resolved emission distributionalong the molecule axis measured for several excitation pow-ers; d)-f) Calculated emission patterns for several excitationpowers. d = 4 µm and d CC = 3 . µm . δ = − meV uously tuned by changing the center to center distance d CC , as illustrated in Fig. 1.b. The left column of Fig. 1.cshows the spatial distribution (near field) of the four low-est energy modes in a molecule made of two 4 µm pillarswith d CC = 3 . µm (the first two panels correspond tothe B and A modes). The right column shows the calcu-lated polariton wavefunctions considering a confinementpotential taken as infinite outside the photonic moleculeand equal to zero inside. As schematically illustrated inFig. 1.d, all these states result from the hybridization ofthe optical modes of each individual micropillar, showingthe strong analogy between our system and the orbitalsof a diatomic molecule.We now discuss polariton condensation in thesemolecules under non resonant optical excitation. As wehave shown in our recent work on single micropillars [21],such an excitation scheme not only populates the con-fined polariton states, but also creates a population ofuncondensed excitons called the excitonic reservoir. Be-cause of the limited diffusion length of excitons, thisreservoir remains localized in the excitation area. Re-pulsive interactions between polaritons and the excitonicreservoir strongly influence the precise quantum state inwhich polaritons accumulate. In the following we willconsider two different locations of the excitonic reservoirwithin the molecule, which can be selected via the posi-tion of the excitation laser spot. First we will considerexcitation conditions where the reservoir of uncondensedexcitons is at the center of the molecule. Then we willaddress the case of asymmetric excitation, in which theexcitonic reservoir is injected only in one of the two mi-cropillars. Specific spatial behavior of the polariton con-densates is observed in each case, driven by the polaritoninteraction with the excitonic reservoir.Figure 2.a-c present the measured emission distribu-tion along the axis of a photonic molecule excited at its FIG. 3: (Color online)
Excitation at the center of themolecule : a) Integrated intensity measured as a function ofthe excitation power; b) Splitting between the two lowest en-ergy emission lines measured as a function of the excitationpower; c) Lowest energy emission lines measured with a polar-ization parallel (solid line) and perpendicular (squares) to themolecule axis for P = 0 . P th ; d) same as c) for P = 7 P th .The inset shows the calculated spatial shape of the bindingstate center. As summarized in Fig. 3.a, above a well definedexcitation threshold, a strong non-linear increase of theemission intensity is observed, one of the signatures ofpolariton condensation [11]. Under this excitation con-dition, a massive accumulation of polaritons occurs bothin the binding and anti-binding lowest energy polaritonstates. Interestingly, as the excitation power is increased,a progressive change of the spatial shape of the polari-ton wavefunction is observed, with vanishing probabilitydensity at the center of the molecule. This behavior isdue to repulsive interactions with the excitonic reservoirinjected in that region.To describe our experiments, we simulate polariton re-laxation using coupled semi-classical Boltzman equationsand nonlinear Schr¨odinger equations in a self-consistentway [22]. As described in the supplementary informa-tion, on one hand, a nonlinear Schr¨odinger equation issolved at each time iteration, in the presence of the po-tential induced by the excitonic reservoir and by the oc-cupation of the polariton states, to find the eigenstatesof the system. On the other hand, the time evolution ofthese occupation numbers are calculated using the semi-classical Boltzmann equations in which the scatteringrates are obtained using the exact shape of the polaritonstates obtained from the non-linear Schr¨odinger equa-tion. Both exciton-phonon and exciton-exciton scatter-ing mechanisms are considered. The reservoir is assumedto be thermalized at the lattice temperature. The key in-gredient to describe our experiments is to consider scat- tering rates proportional to the spatial overlap betweeninitial and final states. For instance, the scattering rateof two excitons from the reservoir resulting in one po-lariton in state i and one exciton in the reservoir is pro-portional to R | ψ i ( x, y ) | ( ρ R ( x, y )) dxdy , where ψ i ( x, y )is the wavefunction of the polariton state i and ρ R ( x, y )is the spatial distribution of the excitonic reservoir, de-scribed by a gaussian of width w R . To fit the data weuse an exciton-exciton scattering rate W XX = 2 10 s − and an exciton-phonon scattering rate W XP = 10 s − .The lifetime of the reservoir is equal to 400 ps , and thatof the polariton states to 30 ps , a value extracted fromcoherence measurements in ref.[12].The results of the model for central pumping are shownin Fig. 2d-f. The potential of the reservoir increases withthe optical pumping and has a maximum at the centerof the molecule, inducing the spatial separation of thestates of each micropillar. As a result, both the cou-pling and the associated splitting between the two lowestenergy states decrease. Indeed, this feature is shown inFig. 3.b: simultaneous to the spatial separation of the po-lariton distributions, the measured splitting goes down toa minimum at about P = 4 P th . At higher powers, themeasured splitting tends to increase again. In the follow-ing, we will show that at such high density the observedsplitting is no longer related to the binding anti-bindingsplitting but is an anisotropic splitting induced by thelateral spatial shrinking of the condensate wavefunctionin each micropillar. This feature is not reproduced by ourmodel which does not include the polarization degree offreedom.Polarization resolved measurements were performedfor different excitation densities. Figure 3.c and 3.d dis-play emission spectra linearly polarized parallel and per-pendicular to the molecule axis, below and above thresh-old. In the low density regime (Fig. 3.c), the presenceof uncondensed excitons has a negligible effect, and weobserve a polarization splitting in both the binding andanti-binding states. As previously reported theoreticallyand experimentally [23], this splitting is larger for thebinding state and amounts to 70 µeV . The inset ofFig.3.c shows that the overall shape of the polariton bind-ing state is elongated along the molecule. This is whythe lowest energy line of the binding state doublet is po-larized parallel to the molecule axis. The situation isdifferent at high density ( P > P th ), where a strongerpotential barrier is induced by the excitonic reservoir atthe center of the molecule. In this case, the coupling be-tween the states of each micropillar becomes too small tobe resolved in our experiment. The remaining splitting isthen purely a polarization splitting. In each micropillar,the polariton wavefunction is shrunk by the excitonic re-pulsive potential and thus it is more elongated perpendic-ularly to the molecule axis. This change in anisotropy isevidenced in the polarization of the emission.The lowestenergy emission line is now strongly polarized perpendic- FIG. 4: (Color online)
Asymmetric Excitation : a)-c)Spectrally resolved emission distribution along the moleculeaxis measured for several excitation powers; d)-f) Calculatedemission patterns for several excitation powers. d = 3 . µm and d CC = 3 . µm . δ = − meV ular to the molecule axis.Different dynamics is observed when the laser spot ispositioned on top of one of the micropillars forming themolecule. Emission distributions measured under theseconditions are reported in Fig. 4.a-c. Close to thresh-old, condensation occurs in the pillar that is opticallypumped. Interparticle interactions blue shift the energyof the ground state in the pumped pillar, decoupling itfrom the states of the other pillar , and severely limitingthe Josephson transfer of particles from one pillar to theother. This localisation is very analogous to the origi-nal self-trapping effect, except that here the blue shiftis mainly induced by the interaction between the exci-ton reservoir and the condensate and not by interactionswithin the condensate itself. As the excitation poweris increased, the better relaxation kinetic and the pres-ence of excited states in the un-pumped pillar destroythe metastable state[24] and allow the system reachingits ground state. Condensation then occurs in the non-pumped micropillar (Fig. 4.c).We can reproduce this overall behavior using the self-consistent polariton relaxation model presented above.As shown in Fig. 4.d, at low power the reservoir poten-tial creates an asymmetry in the two bottom states of themolecule, whereas the upper states are not affected. Thepopulations of the excitonic states in the reservoir aresmall and their relative distribution is given by the latticetemperature. When pumping increases (Fig. 4.e), con-densation starts at the eigenstate possessing the largestoverlap integral with the reservoir. This corresponds tothe third polariton state (population n ). The reser-voir potential also significantly perturbs the two bottomstates, which become both confined in the right pillar.This blocks the direct relaxation from the reservoir to- wards these states and further enforce condensation inthe third metastable state localized in the pumped pil-lar. At high pumping (Fig. 4.f) a new effect is observed:condensation is triggered onto the ground state (of popu-lation n ). This does not arise from direct scattering fromthe reservoir, since the overlap between the ground stateand the reservoir remains very small. This condensationis the result of efficient scattering from the intermediatethird state toward the ground state via acoustic phonons.This regime occurs when the phonon assisted scatteringrate W XP n (1 + n ) from the intermediate state into theground state overcomes the ground state radiative losses n τ . This happens when n ≈
100 and further increaseof the pumping does not change n , but only increases n .The second quantized state remains weakly populated:it has a smaller overlap with the reservoir than the thirdstate and a higher energy than the ground state. Over-all the complete condensation dynamics can only be un-derstood considering the modifications of the relaxationrates induced by interactions (as shown in the supple-mentary [22] when constant scattering rates are consid-ered, condensation in the metastable states can not bereproduced).To conclude, we have demonstrated polariton conden-sation in photonic molecules, a fully controlled systemin which the coupling constant between two condensatescan be adjusted. Depending on the precise location ofthe excitonic reservoir, strong renormalization of the po-lariton states is induced, resulting in strong modificationof the condensation dynamics. In particular condensa-tion can occur in a metastable localized state , or inthe ground state of the system. These results open theway towards a detailed investigation of Josephson oscil-lations [25] or, more generally, the pumped-dissipativephysics of arrays of coupled condensates in an engineeredenvironment.This work was supported by the C’Nano Ile de France(Sophiie2), by the ANR (PNANO- 07-005 GEMINI), bythe ”Triangle de la Physique” (Picorre), by the FP7ITN ”Clermont4” (235114) and by the FP7 ITN ”Spin-Optronics”(237252). [1] P. Werner et al., Phys. Rev. Lett. , 047201 (2005).[2] D. Gerace et al., Nature Physics , 281 (2009)[3] A. D. Greentree et al., Nature Phys. , 856 (2006).[4] S. Schmidt et al., Phys. Rev. B , 100507(R) (2010).[5] M. J. Hartmann et al., Nature Phys. , 849 (2006).[6] I. Carusotto et al., Phys. Rev. Lett. , 033601 (2009).[7] J. Kasprzak et al., Nature , 409 (2006).[8] R. Balili et al., Science , 1007 (2007).[9] C.W. Lai et al., Nature , 529-532 (2007).[10] S. Christopoulos et al., Phys. Rev. Lett. , 126405(2007).[11] D. Bajoni et al., Phys. Rev. Lett. , 047401 (2008).[12] E. Wertz et al., Nature Phys. ,860 (2010). [13] E. A. Cerda-M´endez et al., Phys. Rev. Lett. , 116402(2010)[14] F. Tassone and Y. Yamamoto Phys. Rev. B , 10830(1999)[15] C. Ciuti et al., Phys. Rev. B , R4825 (2000)[16] M. Bayer et al., Phys. Rev. Lett. , 2582 (1998)[17] G. Guttroff et al., Phys. Rev. E , 036611 (2001).[18] A. Smerzi et al., Phys. Rev. Lett. , 4950 (1997).[19] D. Sarchi et al., Phys. Rev. B , 125324 (2008)[20] I. A. Shelykh et al., Phys. Rev. B , 041302 (2008) [21] L. Ferrier et al., Phys. Rev. Lett. , 126401 (2011).[22] Description of the model used to calculate polariton re-laxation is found in the supplementary information.[23] S. Michaelis de Vasconcellos et al., Appl. Phys. Lett. ,101103 (2011)[24] M.T. Trujillo-Martinez et al., Phys. Rev. Lett. ,105302 (2009).[25] K. G. Lagoudakis et al., Phys. Rev. Lett. , 120403(2010). r X i v : . [ c ond - m a t . o t h e r] O c t Supplementary Information on ”Polariton Condensation inphotonic molecules”
Marta Galbiati, Lydie Ferrier, Dmitry D. Solnyshkov, Dimitrii Tanese, Esther Wertz, Alberto Amo, Pascale Senellart, Isabelle Sagnes, AristideLemaˆıtre, Elisabeth Galopin, Guillaume Malpuech, and Jacqueline Bloch CNRS-Laboratoire de Photonique et Nanostructures,Route de Nozay, 91460 Marcoussis, France LASMEA, CNRS, Clermont University, University Blaise Pascal,24 avenue des Landais, 63177 Aubi`ere cedex, France (Dated: October 10, 2018)
Abstract
We describe in this supplementary information the model we have developed to simulate polaritonrelaxation in photonic molecules under the excitation conditions described in the paper. The resultsof these calculations are shown on figure 2 and 4 of the paper. Finally we show that the overallcondensation dynamics is not reproduced if one does not take into account the interaction inducedchanges in the relaxation rates.
PACS numbers:
1e simulate polariton relaxation in a self-consistent way, using coupled semi-classicalBoltzmann equations and nonlinear Schr¨odinger equations. At t = 0 the system is emptyand the pumping is switched on. We calculate the time evolution of the system until theestablishment of a steady state. The nonlinear Schr¨odinger equation written below is usedto find the spatial shapes and energies E i,t j of the polariton wave functions ψ i,t j ( x, y ). E i,t j ψ i,t j = − ~ m ∆ ψ i,t j + αρ R ( x, y ) N R ( t j ) + α (cid:12)(cid:12)(cid:12)(cid:12)(cid:12)X i ψ i,t j − e − iE i,tj − t j / ~ (cid:12)(cid:12)(cid:12)(cid:12)(cid:12) ψ i,t j (1)m is the polariton effective mass. This equation is solved at time step t j to findthe eigenstates in presence of the potential created by the exciton reservoir (written as αρ R ( x, y ) N R ( t j )) and by the occupation of each polariton quantized states (last term of theequation). α = 6 xE b ( a B ) is the exciton-exciton interaction constant, where x the excitonfraction of the considered quantum state (1 in the reservoir, 0.5 in the polariton states), E b the exciton binding energy, a B the exciton Bohr radius. ρ R ( x, y ) is the normalized reservoirdensity ( R ρ R ( x, y ) dxdy = 1), N R ( t j ) the number of excitons in the reservoir.The reservoir is assumed to be thermalized at the lattice temperature: n R ( E, t j ) = n R (0 , t j ) exp ( − E/k B T ). The total number of particles in the reservoir is N R ( t j ) =2 P E n R ( E, t j ), where 2 stands for the spin degeneracy of bright heavy-hole excitons. Sinceexcitons in the reservoir are confined in a quantum well, N R ( t j ) = n R (0 , t j ) Smk B T / π ~ ,where S is the surface covered by the reservoir. ρ R ( x, y ) is assumed to be a Gaussian functionwith a width w R close to the one of the excitation laser spot.The time evolution of n i , the occupation numbers of the quantized polariton states and of N R , the reservoir population are calculated using semi-classical Boltzmann equations, takinginto account both exciton-exciton and exciton-phonon interactions, as well as the presenceof the thermalized excitonic reservoir.For the quantized polariton states i the equation reads:2 n i dt = − n i τ i − n i P k (cid:0) W XX I ikRR N R + W ikXP (cid:1) ( n k + 1) , k < ie − E ik /k B T , k > i +( n i + 1) P k (cid:0) W XX I ikRR N R + W ikXP I ik (cid:1) n k e − E ik /k B T , k < i , k > i + (cid:0) W XX I iRRR N R + W iRXP I iR (cid:1) (cid:0) N R ( n i + 1) − n i Smk B T π ~ e − E iR /k B T (cid:1) and for the reservoir state R the equation is: dN R dt = P − N R τ R − X i (cid:0) W XX I iRRR N R + W iRXP I iR (cid:1) (cid:18) N R ( n i + 1) − n i Smk B T π ~ e − E iR /k B T (cid:19) (2) P is the non-resonant pumping populating the reservoir, τ R = 400 ps (resp. τ i = 30 ps ) is the excitonic reservoir (resp. polariton) lifetime. W XX describes the exciton-excitoninteraction rates. The exciton-phonon scattering rate reads: W ik,RXP = W XP e − ( Ei − Ek,R ) ~ c πaB The wave vector conservation rule is replaced by an overlap integral between the spatialextensions of the initial and final states. The overlap integrals change with time, togetherwith the shapes and the energies of the states according to equation (1). The overlapintegrals involved are defined as follows: I ikRR = S Z | ψ ∗ i ψ k | ρ R dxdy,I iRRR = S Z | ψ i | ρ R dxdy (3)The overlap integrals obtained by summation of scattering rates with different phononwavevectors read: I ik = S Z | ψ ∗ i ψ k | dxdy,I iR = S Z | ψ i | ρ R dxdy (4)These overlap integrals are dimensionless.Finally, the scattering rates W XX and W XP describing exciton-exciton and exciton-phonon interaction were used as the fitting parameters of the model, together with thewidth of the gaussian spatial distribution of the reservoir which we authorize to slightly3 IG. 1: Schematic representation of the exciton-exciton interaction in the reservoir leading to thescattering towards the quantized polariton states. The reservoir has a 2D (constant) density ofstates with a thermal distribution function. vary above the width given by the pumping laser spot size. We used W XX = 2 10 s − , W XP = 10 s − , w r = 3 µm .In order to understand the terms in the Boltzmann equation, let us consider for examplethe case of the scattering from the reservoir towards a quantized state i and backward assistedby exciton-phonon and exciton-exciton interaction. We neglect the final state stimulationwithin the reservoir which means that we assume a weak occupation ( n R ( E ) ≪