Abstract
Phase transitions, where observable properties of a manybody system change discontinuously, can occur in both open and closed systems. By placing cold atoms in optical cavities and inducing strong coupling between light and excitations of the atoms, one can experimentally study phase transitions of open quantum systems. Here we observe and study a nonequilibrium phase transition, the condensation of supermodedensitywave polaritons. These polaritons are formed from a superposition of cavity photon eigenmodes (a supermode), coupled to atomic density waves of a quantum gas. As the cavity supports multiple photon spatial modes and because the light–matter coupling can be comparable to the energy splitting of these modes, the composition of the supermode polariton is changed by the light–matter coupling on condensation. By demonstrating the ability to observe and understand densitywavepolariton condensation in the fewmodedegenerate cavity regime, our results show the potential to study similar questions in fully multimode cavities.
Introduction
Ultracold atoms have provided an exemplary model system to demonstrate the physics of closedsystem, equilibrium phase transitions, confirming many theoretical models and results^{1}. A striking manifestation of the role quantum mechanics can play in the physics of equilibrium phase transitions is the wellknown Bose–Einstein condensation (BEC) of bosonic particles at low temperatures into a single, macroscopically populated quantum wave. Our understanding of dissipative phase transitions in quantum systems is less developed and experiments that probe this physics even less so. For example, condensation in quantum systems out of thermal equilibrium is far less understood than in equilibrium^{2,3}, yet becoming experimentally relevant, especially via the study of polariton condensates: when matter couples strongly to light, new collective modes called polaritons arise. Condensation of these quasiparticles has been actively studied in the form of exciton polaritons^{4,5,6,7,8,9,10,11}.
With only a single mode of light coupled to the BEC, selforganization has been observed, but is uniquely defined by cavity geometry alone and is best described as a densitywave (DW) polariton condensate^{12,13,14}. In contrast, for exciton polaritons, complex pattern formation^{8,9,10} has been seen, arising from the existence of many nearly degenerate modes in planar microcavities. In addition, related physics has also been observed in systems in which no DWpolariton condensation occurs: for example, superradiant emission using hyperfine states^{15}, supermode emission without selforganization^{16} and selforganization without a cavity^{17}. The Bose–Hubbard model with infiniterange interactions has been studied^{18,19}. We note that condensation of DW polaritons is quite distinct from lasing: standard lasing requires inversion of the gain medium, is driven by an incoherent pump and is generally incompatible with strong light–matter coupling, whereas the transition we study differs on all of these points.
We report a step towards the observation of an unusual form of nonequilibrium condensation, a ‘supermodeDWpolariton’ condensate. Here, supermode refers to the eigenmode built as a superposition of the bare cavity modes^{20}. The supermodepolariton dressed state is dependent on the pumpcavity detuning Δ_{c} and on the overlap of the barecavity modes with the BEC position and shape. The matter component is an atomic DW excitation rather than the electronic excitation of excitonpolariton condensates. It is important to distinguish the condensation of supermode DW polaritons from condensation of the atoms; for example, selfconsistent formation of atomic DWs and cavity light can be studied with thermal atoms^{21,22} and with BECs^{12}. However, some of the signatures we demonstrate, such as the atomic structure factor, are observable only with BECs of atoms in multimode cavities (see Methods for further discussion). By embedding the exquisite control available for ultracold atoms within multimode quantumoptical systems, our experiment opens avenues for experimentally studying quantum fluctuationdriven transitions and quantum criticality^{2,23,24,25}. In such a system, one can expect physics beyond meanfield theories such as a quantum Brazovskii transition or a shortrange spin glass^{23,26}. Thus, these systems will provide experimental access to nontrivial phase transitions in driven dissipative quantum systems^{2,3,23,24,27} and enabling the studies of exotic nonequilibrium spin glasses and neuromorphic computation^{26,28}.
Results
Supermode DW polaritons
In the work reported here, we see how with a few modes coupled to the BEC (see Fig. 1b,c) the supermodeDWpolariton condensate arises in a regime intermediate between BECs coupled to a singlemode cavity and those coupled to a confocal or concentric cavity supporting many hundreds of degenerate modes. When the atoms are pumped with a laser orthogonal to the cavity axis (see Fig. 1a), the normal modes of the atomcavity system evolve to become supermode DW polaritons, new superpositions of supermodes mixed by the DW fluctuations of the atoms. Above a critical pump threshold, we see defining characteristics of supermodeDWpolariton condensates, heralded by three observables as follows: first, superradiant emission of light from the cavity with the spatial pattern of one of these new supermodes; second, Z_{2} symmetry breaking of the phase of the cavity field, locking to either or +π with respect to the pump phase; and third, organization of the BEC wavefunction into one of two checkerboard lattice configurations—each corresponding to a specific phase of the cavity field. Such observables have been seen in singlemode cavities^{12}. In addition, we see how nontrivial transverse spatial structure of the supermode can result in lattice defects (matter–wave phase slips): these are seen in the structure factor in an atomic timeofflight measurement. Observations of all these defining characteristics are presented.
The nearconfocal optical cavity employed here supports families of optical modes that each lie within a small frequency bandwidth, as shown in Fig. 1b,c (see Methods). One can observe the superradiant emission of various supermodeDWpolariton condensates by pumping at different Δ_{c} tuned near or within a mode family, as can be seen in Fig. 1e,j. This is in contrast to the DWpolariton condensate of a singlemode cavity such as shown in Fig. 1d: This is not a supermode no matter the detuning Δ_{c} near this isolated Gaussian mode. The supermodes in Fig. 1e–j, differ from ideal Hermite–Gaussian modes due to three factors as follows: first, the bare cavity modes are themselves mixtures of ideal Hermite–Gaussian modes (due to mirror aberrations mixing modes when there is spectral overlap of modes near degeneracy)^{29}; second, these bare cavity modes are mixed by the dielectric atomic medium to form supermode polaritons; and third, these dressed states are remixed by the emergent DW to form new supermodes above the polariton condensation threshold.
Decomposition of photonic components
Figure 2 illustrates how remixed photonic components of the supermodeDWpolariton condensate can differ from the supermodepolariton dressed states. The three bare cavity modes of the l+m=2 family are mixed by a BEC at the cavity centre to produce the three supermode peaks. The photonic component, shown in Fig. 2h–l, of the supermodeDWpolariton condensates can differ from that of the supermode polaritons in Fig. 2c–g due to new supermode mixing by the macroscopically populated atomic DW above threshold. This is most pronounced away from resonance: see, for example, how Fig. 2c,h differ; the supermodeDWpolariton condensate at 46 MHz (Fig. 2h) is ∼81% TEM_{02} and 19% TEM_{20}, whereas the associated belowthreshold supermode polariton (Fig. 2c) is ∼52% TEM_{02}, 9% TEM_{20} and 39% TEM_{11}. The suppression of the TEM_{11} component above threshold can be understood as resulting from its poor overlap with the BEC and thus weaker mixing with the DW mode. A similar remixing of supermodes occurs on the bluedetuned side at 56 MHz.
The components closely follow the theory prediction based on a linear stability analysis for mode content at condensation threshold, except near the mode. We believe this discrepancy is due to dynamical effects not captured by this static stability analysis. See Fig. 8 in Methods for discussion.
Position dependence
For the limit we consider, where the atomic cloud is smaller than the beam waist, the position of the BEC with respect to the barecavity modes strongly affects the photonic mode composition of the supermodeDWpolariton condensate. That is, the neardegeneracy of the bare cavity modes means that the particular supermode selected is determined by overlap with the cloud. This is easily observed by moving the BEC within the transverse plane of the cavity with the pump tuned near the l+m=1 family (see Fig. 3). We pump the system with the BEC trapped by the optical dipole trap (ODT) at each of the four intracavity positions illustrated in Fig. 3a. With the BEC trapped near either the antinode of the TEM_{10} or the TEM_{01} mode, we observe superradiant emission with a spatial pattern nearly identical to these barecavity modes, as shown in Fig. 3c,d. However, a BEC at the intersection between the two modes’ antinodal lobes (Fig. 3b) yields an emitted spatial pattern at 45° to the bare cavity mode axes. The threshold for organization is the same for a BEC at a barecavity antinode or between the antinodes, demonstrating that the BEC has mixed these bare modes equally and has created a basis for supermode DW polaritons that are rotated 45° from the original l+m=1 family eigenbasis; see Methods.
Structure factor
To complete the description of this nonequilibrium condensate, we report that observations of the momentum distribution in timeofflight reveal the influence of the supermode structure on the matter wave component of the polariton condensate. For singlemode cavities pumped near the TEM_{00} mode, the atoms organize in one of two possible checkerboard patterns. However, organization in cavities supporting higherorder modes is more complicated. As illustrated in Fig. 4a versus Fig. 4b when the BEC overlaps with a node of the cavity mode, the effective DW—cavity mode coupling changes sign across the node, because the DW couples to the interference between pump field and cavity mode. This results in an organized state with a plane defect in the checkerboard lattice: that is, a π phase slip, and we confirm this configuration to be the optimal organized state via opensystem simulations of the pumped BECcavity system. See Fig. 4d and Methods.
Timeofflight expansion of the atoms yields the atomic momentum distribution and the lattice defect appears as a node in the (k_{y}, k_{z})=(±1, ±1)k_{r} Bragg peaks of the expanding BEC’s interference pattern; see simulation in Fig. 4g. The node may be understood as a structure factor resulting from the lowmomentum modulation of the organized atomic wavefunction caused by coupling to the transverse nodal structure of the supermode. The observability of a structure factor in the momentum distribution of BECs organized in orthogonally oriented and higherorder modes is shown in Fig. 5.
Stronger pumping modifies the matter wavefunction by increasing the light–matter coupling nonlinearity, as may be seen by the emerging node at the centre of the zerothorder and (±2, 0) Bragg peaks in Fig. 4k. This distortion of the condensate wavefunction is similar to that which happens in dilute gas BECs on increasing interaction energy. Lastly, we note that these supermodeDWpolariton condensates are observed to break the same Z_{2} symmetry observed in singlemode DWpolariton condensates^{21,22,30}. Figure 6 presents measurements of pumpcavity field phase locking.
Discussion
Our results demonstrate it is possible to study the DW polariton condensation phase transition when there are a few degenerate (or nearly degenerate) cavity modes. We have shown this has notable effects on properties such as the structure factor in the timeofflight image of atomic density. The demonstration of this phase in the fewmodedegenerate system paves the way for measurements of critical behaviour in nonequilibrium quantum systems employing fully multimode cavity quantum electrodynamics. The multimode regime can be easily realized in our existing apparatus by tuning the cavity mirror spacing to confocality in situ. Indeed, this capability has already been achieved in our system in a robust manner, as discussed in ref. 31. To explore critical behaviour in phase transitions of open quantum systems, the crucial requirement is the existence of a continuum of modes, so that the response of the system at different length and timescales is not all dependent on the same few modes. The number of modes should be large enough that the beyondmeanfield corrections to the threshold power are greater than the ability to measure and control this threshold power.
Multimode cavity quantum electrodynamics, in the limit of multimode collective ultrastrong coupling wherein the collective coupling is larger than the bandwidth of the degenerate modes, strongly mixing them, should provide access to a much more exotic condensation transition. This fluctuationinduced firstorder Brazovskii transition is predicted to yield a superfluid smecticlike quantum liquid crystalline order of the intracavity BEC^{23}. This opens avenues to study the interplay of quantum liquid crystallinity and unconventional superfluidity under controlled dimensionality and disorder, as well as the study of superfluid glasses and spin glasses^{23,32}, longstanding problems in statistical mechanics.
Methods
Apparatus
We prepare a nearly pure ^{87}Rb BEC of 3 × 10^{5} atoms at the centre of our cavity, confined in a crossed ODT with trap frequencies [ω_{x}, ω_{y}, ω_{z}]=2π × [59.1(6), 88.0(9), 89.4(5)] Hz. The atoms are prepared in the state and a 1.4 G magnetic field is oriented along the z axis. The Thomas–Fermi radii of the BEC [R_{x}, R_{y}, R_{z}]=[9.8(3), 8.3(2), 8.3(2)] μm are significantly smaller than the 35 μm waist (1/e radius of the cavity field) of the TEM_{00} cavity mode. The crossed ODT is formed by a pair of 1,064 nm laser beams with waists 39 and 20 μm intersecting at 45° in the xy plane. Acoustooptic modulators are used to stabilize the intensity of each ODT beam and control its position, allowing us to translate the BEC inside the cavity to control its overlap with the cavity modes.
The cavity is operated in a nearconfocal regime in which the length L is set to differ from the radiusofcurvature by 50 μm (ref. 31). The L=1 cmlong cavity has a freespectral range of 15 GHz and a singleatom TEM_{00} cooperativity of 2.5: g_{0}=2π × 1.04 MHz and κ=2π × 132 kHz. A weak 1,560 nm laser is used to stabilize the cavity length using the Pound–Drever–Hall method. In addition, light from this laser is amplified and doubled to generate 780 nm light for the transverse pump and longitudinal probe beams. The wavelength of the locking laser is chosen to achieve a large atomic detuning of Δ_{a}=−102 GHz between the pump and the 6 MHzwide D2 line of ^{87}Rb. An electrooptic modulator placed in the path of the locking laser beam allows us to tune the detuning Δ_{c} between the cavity modes and the pump or probe beams. The pump beam is polarized along x and is focused down to a waist of 80 μm at the BEC and retroreflected to create an optical lattice oriented along y. At the cavity mirror, the TEM_{00} beam has a waist of 50 μm, which is much smaller than the mirror size. We have checked elsewhere^{31} that this remains true up to modes of order l, m≃50, and that these modes also exhibit high finesse. This provides an ultimate limit of around 1,250 on the number of transverse modes that can be made resonant exactly at the confocal point.
The cavity modes can be probed using a large, nearly flat longitudinal probe beam propagating along the axis of the cavity. This probe couples to all transverse modes of the loworder families. Reference ^{31} describes the design and vibration isolation of the lengthadjustable cavity.
Measurements
The cavity output can be directed to three different detection channels. A singlephoton counting module can record photon numbers via a multimode fibre coupled to the multimode cavity output. We measure a detection efficiency—from cavity output mirror to detector, including quantum efficiency and losses—of 10% for the loworder modes discussed here. The dispersive shift data in Fig. 2a is taken for an intracavity photon number much less than one and the same Δ_{a} as for the rest of the data, −102 GHz.
Superradiant emission from the cavity is observed by monitoring the cavity output on the singlephoton counting module. A sharp rise in intracavity photon number as the pump power is increased heralds the condensation transition. See Figs 7 and 8.
Alternatively, we can image the emission using an electronmultiplying CCD (chargecoupled device) camera, to spatially resolve the transverse mode content of the cavity emission, although with no temporal resolution. All images of cavity emission are taken in a single experimental run (no averaging of shots with different BEC realizations) and with a camera integration time between 1 and 3.3 ms. The lower signaltonoise in the images of Fig. 2c–g versus Fig. 2h–l is due to lower intracavity photon number. The images in Fig. 2b with no atoms present are taken with intracavity photon number well above unity. The pump power in Fig. 4k is 70% larger than in Fig. 4j. See Fig. 9 for images of higherorder modes.
The phase difference between the cavity output and pump beam can be determined by performing a heterodyne measurement with a local oscillator beam^{22}. The data in Fig. 6 has a frequency offset between pump and signal beams of 11 MHz.
We calibrate the Rabi frequency of the transverse pump by measuring the depth of the pump lattice through Kapitza–Dirac diffraction of the BEC. With the cavity modes detuned faroff resonance, we pulse the pump lattice onto the BEC for a time Δt. The BEC is then released from the trap and the population P_{m}(Δt) of the diffracted orders at momenta is measured after a timeofflight expansion. By fitting the measured P_{m}(Δt) to theory^{33}, we extract a lattice depth V_{0}=Ω^{2}/Δ_{a}, which allows us to determine the Rabi frequency Ω of the pump.
The momentum distribution of the ^{87}Rb cloud is measured by releasing the cloud from the trap and performing resonant absorption imaging after an expansion time t_{TOF}=17 ms. The appearance of Bragg peaks at k=k_{r} in the atomic momentum distribution coincides with the onset of superradiant emission. All atomic timeofflight images are taken in a single experimental run (no averaging of shots with different BEC realizations). The paired atomic absorption and cavity output images in Fig. 4i,j and in Fig. 5a–f are each taken for the same experimental run. The cavity output image in Fig. 4l was taken for the same experimental run as the atomic absorption image in Fig. 4k.
We use the position dependence of the l+m=1 supermodeDWpolariton condensates (see Fig. 3) to place our BEC at the centre of the cavity modes. The BEC is translated in the xy plane using the acoustooptic modulators of the two ODT beams. By monitoring the orientation of the superradiant emission as a function of BEC position, we are able to infer the displacement between the BEC and cavity centre. Exploiting these effects allows us to position the BEC at the cavity centre in both x and y directions to within 4 μm.
Mode decomposition
We analyse the transverse mode content of the polaritons in the l+m=s family by decomposing the cavity field into a superposition of unitnormalized, Hermite–Gaussian modes Φ_{lm}(x, y; w_{0}, x_{0}, y_{0}). This is achieved by fitting the EMCCD image I(x, y) to the function
with the mode magnitudes A_{lm} and phases φ_{lm} determined as fit parameters. The waist w_{0} and centre positions (x_{0}, y_{0}) of Hermite–Gaussians are determined from an image of the TEM_{00} cavity mode and are held fixed during the fit. Fixing the phase of the Φ_{s0} modes to φ_{s0}=0 allows the fitting algorithm to converge to a local optimum. Using the fitted values of A_{lm}, we extract admixture fractions
of the Φ_{lm} mode in the cavity output.
Model Hamiltonian
We model the coupled dynamics of the atomic wavefunction and cavity modes (where μ=(l, m)) with the Hamiltonian
The first term represents the evolution of the cavity modes and second is the familiar Gross–Pitaevskii Hamiltonian for a weakly interacting BEC trapped in a harmonic potential V(r). The atomic contact interaction is accounted for in the term proportional to U. In the regime of large atomic detuning Δ_{a}, we can neglect the excited electronic state of the atom, and the atomlight interactions may be described solely through the dispersive light shifts. The cavity–atom interaction becomes
where g_{ν}(r)=g_{0}Φ_{ν}(r)/Φ_{00}(0) is the spatially dependent, singlephoton (vacuum) Rabi frequency for the cavity mode ν. Similarly, for a pump field with a Rabi frequency Ω(r), the atom–pump interaction is
The last term of equation (3) represents the light shift arising from the interference between the cavity and pump fields and is written as
Simulation of supermode composition at threshold
To predict the location of threshold, and the nature of the supermode at that point, one may make use of a linear stability analysis, assuming a small occupation of the cavity modes and atomic DW excitation^{34,35}. For the cavity mode, this is straightforward. For the atoms, this corresponds to assuming a condensate wavefunction
where μ_{n}(θ) are 2π periodic eigenfunctions of the Mathieu equation, with eigenvalue a_{n}, that is,
where q=−E_{Ω}/ω_{r} and E_{Ω}=Ω^{2}/Δ_{a}. These Mathieu functions describe the effects of the pump beam in the x direction and do not assume a weak pump lattice. In the cavity direction, z, the lattice is assumed weak and so we only consider the first two modes, that is, 1 and cos(kx). As the above expression encapsulates all effects of the longitudinal coordinate, we will suppress the label ⊥ on the transverse coordinates.
We must then solve coupled equations for the atomic transverse envelope functions and cavity mode amplitudes α_{μ}. To leading order in perturbation theory, the ground state envelope does not change and so corresponds to the solution of the Gross–Pitaevskii equation:
where μ is the chemical potential and N the number of atoms.
Meanfield equations of motion for the and are derived from the Hamiltonian in equation (3),
where we have taken Ω(r)=Ω cos(kx). The spatial dependence of the pump enters through the overlap , of the first two Mathieu functions due to the cross pump–cavity light field potential. The energy scale ω_{0}(q)=ω_{r}(1+a_{1}(q)−a_{0}(q)) corresponds to the effective recoil in pump and cavity directions, allowing for the possibility of a deep pump lattice. For a shallow lattice these functions become 1 and 2ω_{r}, respectively.
From these linearized equations, we may then determine when supermodeDWpolariton condensation occurs, by identifying the point at which the linearized fluctuations become unstable. There is some subtlety to this point, discussed further below. Calculating the growth/decay rates of linearized fluctuations is straightforward, corresponding to an eigenvalue equation. As there are anomalous coupling terms (that is, as α_{μ} depends on both and , and vice versa) one must use the Bogoliubov–de Gennes parametrization, that is, write and similarly for . It is convenient to resolve the function onto some set of basis states. We use the harmonic oscillator basis states, giving a particularly simple result in the limit U→0.
With the basis noted above, the eigenvalue problem is given by Det[A−λ1]=0 where the matrix A can be written in the block form in terms of blocks:
In this expression, the various block matrices are as follows: the matrix Δ_{c} is a diagonal matrix consisting of the detuning between the pump laser and each cavity mode. Δ_{dw}(q) is similarly a diagonal matrix describing the energy difference between a given atomic transverse mode function and the ground state mode function. This is a function of q as it also includes the energy ω_{0}(q) associated with the different scattering states. The matrices W denote the effect of atomatom interactions, corresponding to the overlap between two atomic modes and the atomic ground state density. That is, they describe scattering off the condensate causing transitions between modes. The matrices M denote the dielectric shift due to the atoms, corresponding to the overlap between two cavity modes and the atomic ground state density. (The matrices M and W differ in general because the cavity beam waist does not match the harmonic oscillator length of the atoms.) The matrix Q denotes atomcavity scattering and involves the overlap of the atomic groundstate mode function with a given excited mode and a given cavity mode. The energy scale is for cavity beam waist w.
The condition for the matrix A having unstable eigenvalues can be directly related to the idea of nonequilibrium condensation of polaritons^{36}. As discussed there, one may directly relate the real frequencies at which the real and imaginary parts of the inverse retarded Green’s function vanish to the polariton energies and the effective chemical potential. At the point of condensation, these frequencies meet; that is, there exists a real frequency at which the inverse Green’s function vanishes. This point is equivalent to the boundary between the Green’s function having unstable and stable poles. For the singlemode cavity, this connection to condensation has been discussed extensively^{24,34}, including the identification of the low energy effective temperature emerging at the transition. The Green’s function which defines all these quantities is directly related to equation (11), although conventional definitions of Green’s functions introduce minus signs in alternate rows and columns.
On solving equation (11), one finds that although there is a threshold for instability when Δ_{c}<0 of nearby modes, there is always an unstable eigenvalue as soon as Δ_{c}>0 for any mode (that is, the pump is bluedetuned of any mode). However, the growth rate of this instability varies widely with parameters. At low pump powers, the timescale for growth is very long (that is, seconds)^{35}. As pump strength increases, there is a sharp threshold where two eigenvalues of equation (11) cross, demarcating a transition to a state that rapidly orders (timescale of microseconds). The curves in Fig. 2 correspond to finding the eigenvector (i.e., mode composition) of this mode which becomes rapidly unstable.
Structure factor simulation
The calculated density and momentum distributions presented in Fig. 4 are evaluated by numerically integrating the meanfield equations of motion in three dimensions for the atomic wavefunction and cavity mode. Using the Hamiltonian in equation (3), we derive equations of motion for the atomic wavefunction and cavity field under a meanfield approximation where and . This gives us the coupled differential equations,
Furthermore, we adiabatically eliminate the cavity field α_{μ} under the assumption that it equilibrates on a timescale much faster than the atomic motion. To simulate the behaviour presented in Fig. 4, we restrict α_{μ} to a single mode, either TEM_{00} or TEM_{10}, and numerically integrate the equations of motion. The initial atomic wavefunction is set to the Thomas–Fermi distribution associated with our ODT parameters and the initial cavity field is set to α_{μ}(0)=0. The strength of the pump field is increased linearly in time from 0 to Ω(r, t) to simulate the transverse pumping of our cavity. The in situ density distributions shown in Fig. 4c–e are . The momentum distributions are obtained by a Fourier transform of the in situ atomic wavefunction.
Condensation of supermode DW polaritons
The model described above is specific to the case of condensation of supermode DW polaritons occurring on top of an existing BEC of atoms, as was studied experimentally. The phase transition to macroscopic occupation of cavity modes can however happen for both Bose condensed and thermal atoms^{22}. Indeed, for a singlemode cavity, the distinction of the phase transitions for BEC of atoms and superradiance was discussed by Piazza et al.^{37}, leading to a phase diagram where either phase can occur independently of the other. However, the signatures we study, such as the appearance of structure in the Bragg peaks from the timeofflight image, would no longer be visible in the absence of atomic coherence^{18,19}, that is, only for a coherent atomic state is the momentum distribution direction related by Fourier transformation to the density profile in real space.
The relation of the DWpolariton condensation transition to the coherence of atoms is complicated by the interpretation of singlemode DWpolariton condensation^{12} in terms of the Dicke model. This interpretation is based on the similarity to the proposed realisation of the Dicke model^{38} using twophoton transitions between hyperfine states of the atoms, later realized by Baden et al.^{15} In a singlemode cavity, DWpolariton states can be approximately mapped to this model^{12,13,34} as long as the atoms are Bosecondensed—in this case, two macroscopically occupied momentum states of the atoms play the rôle of the internal states in the proposal^{38}. For a noncondensed cloud^{37}, or for multimode cavities, such an approximation does not hold. Indeed, before the experiments on a condensate in a singlemode cavity, supermodeDWpolariton condensation in a multimode cavity had been discussed by Gopalakrishnan et al.^{23,39}
Data availability
Parts of the research data supporting this publication can be accessed from ( http://dx.doi.org/10.17630/eddbeb149bab40588033cc4672770e0c). The remaining data are available from the corresponding author on reasonable request.
Additional information
How to cite this article: Kollár, A. J. et al. Supermodedensitywavepolariton condensation with a Bose–Einstein condensate in a multimode cavity. Nat. Commun. 8, 14386 doi: 10.1038/ncomms14386 (2017).
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
References
 1
Bloch, I., Dalibard, J. & Zwerger, W. Manybody physics with ultracold gases. Rev. Mod. Phys. 80, 885–964 (2008).
 2
Diehl, S., Tomadin, A., Micheli, A., Fazio, R. & Zoller, P. Dynamical phase transitions and instabilities in open atomic manybody systems. Phys. Rev. Lett. 105, 015702 (2010).
 3
Sieberer, L. M., Huber, S. D., Altman, E. & Diehl, S. Dynamical critical phenomena in drivendissipative systems. Phys. Rev. Lett. 110, 195301 (2013).
 4
Deng, H., Weihs, G., Santori, C., Bloch, J. & Yamamoto, Y. Condensation of semiconductor microcavity exciton polaritons. Science 298, 199–202 (2002).
 5
Kasprzak, J. et al. BoseEinstein condensation of exciton polaritons. Nature 443, 409–414 (2006).
 6
Balili, R., Hartwell, V., Snoke, D., Pfeiffer, L. & West, K. Boseeinstein condensation of microcavity polaritons in a trap. Science 316, 1007–1010 (2007).
 7
Amo, A. et al. Collective fluid dynamics of a polariton condensate in a semiconductor microcavity. Nature 457, 291–295 (2009).
 8
Manni, F., Lagoudakis, K. G., Liew, T. C. H., André, R. & DeveaudPiédran, B. Spontaneous pattern formation in a polariton condensate. Phys. Rev. Lett. 107, 106401 (2011).
 9
Tosi, G. et al. Sculpting oscillators with light within a nonlinear quantum fluid. Nat. Phys. 8, 190–194 (2012).
 10
Carusotto, I. & Ciuti, C. Quantum fluids of light. Rev. Mod. Phys. 85, 299–366 (2013).
 11
Gao, T. et al. Observation of nonHermitian degeneracies in a chaotic excitonpolariton billiard. Nature 526, 554–558 (2015).
 12
Baumann, K., Guerlin, C., Brennecke, F. & Esslinger, T. Dicke quantum phase transition with a superfluid gas in an optical cavity. Nature 464, 1301–1306 (2010).
 13
Nagy, D., Konya, G., Szirmai, G. & Domokos, P. Dickemodel phase transition in the quantum motion of a BoseEinstein condensate in an optical cavity. Phys. Rev. Lett. 104, 130401 (2010).
 14
Ritsch, H., Domokos, P., Brennecke, F. & Esslinger, T. Cold atoms in cavitygenerated dynamical optical potentials. Rev. Mod. Phys. 85, 553–601 (2013).
 15
Baden, M. P., Arnold, K. J., Grimsmo, A. L., Parkins, S. & Barrett, M. D. Realization of the Dicke model using cavityassisted Raman transitions. Phys. Rev. Lett. 113, 020408 (2014).
 16
Wickenbrock, A., Hemmerling, M., Robb, G. R. M., Emary, C. & Renzoni, F. Collective strong coupling in multimode cavity QED. Phys. Rev. A 87, 043817 (2013).
 17
Labeyrie, G. et al. ‘Optomechanical selfstructuring in a cold atomic gas. Nat. Photonics 8, 321–325 (2014).
 18
Klinder, J., Keßler, H., Bakhtiari, M. R., Thorwart, M. & Hemmerich, A. Observation of a superradiant Mott715 insulator in the DickeHubbard model. Phys. Rev. Lett. 115, 230403 (2015).
 19
Landig, R. et al. Quantum phases from competing shortand longrange interactions in an optical lattice. Nature 532, 476–479 (2016).
 20
Kapon, E., Katz, J. & Yariv, A. Supermode analysis of phaselocked arrays of semiconductor lasers. Opt. Lett. 10, 125–127 (1984).
 21
Domokos, P. & Ritsch, H. Collective cooling and selforganization of atoms in a cavity. Phys. Rev. Lett. 89, 253003 (2002).
 22
Black, A. T., Chan, H. W. & Vuletić, V. Observation of collective friction forces due to spatial selforganization of atoms: from Rayleigh to Bragg scattering. Phys. Rev. Lett. 91, 203001 (2003).
 23
Gopalakrishnan, S., Lev, B. L. & Goldbart, P. M. Emergent crystallinity and frustration with BoseEinstein condensates in multimode cavities. Nat. Phys. 5, 845–850 (2009).
 24
Dalla Torre, E. G., Diehl, S., Lukin, M. D., Sachdev, S. & Strack, P. Keldysh approach for nonequilibrium phase transitions in quantum optics: beyond the Dicke model in optical cavities. Phys. Rev. A 87, 023831 (2013).
 25
Marino, J. & Diehl, S. Driven Markovian quantum criticality. Phys. Rev. Lett. 116, 070407 (2016).
 26
Gopalakrishnan, S., Lev, B. L. & Goldbart, P. M. Frustration and glassiness in spin models with cavitymediated interactions. Phys. Rev. Lett. 107, 277201 (2011).
 27
Chang, D. E., Cirac, J. I. & Kimble, H. J. Selforganization of atoms along a nanophotonic waveguide. Phys. Rev. Lett. 110, 113606 (2013).
 28
Strack, P. & Sachdev, S. Dicke quantum spin glass of atoms and photons. Phys. Rev. Lett. 107, 277202 (2011).
 29
Siegman, A. E. Lasers Univ. Science Books (1986).
 30
Asbóth, J. K., Domokos, P., Ritsch, H. & Vukics, A. Selforganization of atoms in a cavity field: Threshold, bistability, and scaling laws. Phys. Rev. A 72, 053417 (2005).
 31
Kollár, A. J., Papageorge, A. T., Baumann, K., Armen, M. A. & Lev, B. L. An adjustablelength cavity and BoseEinstein condensate apparatus for multimode cavity QED. New J. Phys. 17, 043012 (2015).
 32
Habibian, H., Winter, A., Paganelli, S., Rieger, H. & Morigi, G. Boseglass phases of ultracold atoms due to cavity backaction. Phys. Rev. Lett. 110, 075304 (2013).
 33
Denschlag, J. H. et al. A BoseEinstein condensate in an optical lattice. J. Phys. B 35, 3095–3110 (2002).
 34
Keeling, J., Bhaseen, M. J. & Simons, B. D. Collective dynamics of BoseEinstein condensates in optical cavities. Phys. Rev. Lett. 105, 043001 (2010).
 35
Bhaseen, M. J., Mayoh, J., Simons, B. D. & Keeling, J. Dynamics of nonequilibrium Dicke models. Phys. Rev. A 85, 013817 (2012).
 36
Szymanska, M. H., Keeling, J. & Littlewood, P. B. in Quantum Gases: Finite Temperature and NonEquilibrium Dynamics 1st edn eds Proukakis N. P., Gardiner S., Davis M. J., Szymanska M. H. Vol. 1, Ch. 30, Imperial College Press (2013).
 37
Piazza, F., Strack, P. & Zwerger, W. BoseEinstein condensation versus DickeHeppLieb transition in an optical cavity. Ann. Phys. 339, 135–159 (2013).
 38
Dimer, F., Estienne, B., Parkins, A. S. & Carmichael, H. J. Proposed realization of the Dickemodel quantum phase transition in an optical cavity QED system. Phys. Rev. A 75, 013804 (2007).
 39
Gopalakrishnan, S., Lev, B. L. & Goldbart, P. M. Atomlight crystallization of BoseEinstein condensates in multimode cavities: nonequilibrium classical and quantum phase transitions, emergent lattices, supersolidity, and frustration. Phys. Rev. A 82, 043612 (2010).
Acknowledgements
We thank S. Gopalakrishnan for stimulating discussions, and J. Witmer, B. Pichler and N. Sarpa for early experimental assistance. Funding was provided by the Army Research Office. B.L.L. acknowledges support from the David and Lucille Packard Foundation, and A.J.K. and A.T.P. acknowledge support the NDSEG fellowship program. J.K. acknowledges support from EPSRC programme ‘TOPNES’ (EP/I031014/1) and from the Leverhulme Trust (IAF2014025).
Author information
Affiliations
Contributions
B.L. and J.K. conceived of project. B.L., A.K. and A.P. designed and built experiment. A.K., V.V. and Y.G. took data. J.K. and A.P. performed simulations. B.L., J.K., V.V. and A.K. wrote the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Kollár, A., Papageorge, A., Vaidya, V. et al. Supermodedensitywavepolariton condensation with a Bose–Einstein condensate in a multimode cavity. Nat Commun 8, 14386 (2017). https://doi.org/10.1038/ncomms14386
Received:
Accepted:
Published:
Further reading

An optical lattice with sound
Nature (2021)

Strongly correlated Fermions strongly coupled to light
Nature Communications (2020)

Cavityless selforganization of ultracold atoms due to the feedbackinduced phase transition
Scientific Reports (2020)

Spontaneous lightmediated magnetism in cold atoms
Communications Physics (2018)

A strongly interacting polaritonic quantum dot
Nature Physics (2018)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.