Skip to main content
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian
  • Log in
  • My Cart

Main menu

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home

Advanced Search

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses

New Research In

Physical Sciences

Featured Portals

  • Physics
  • Chemistry
  • Sustainability Science

Articles by Topic

  • Applied Mathematics
  • Applied Physical Sciences
  • Astronomy
  • Computer Sciences
  • Earth, Atmospheric, and Planetary Sciences
  • Engineering
  • Environmental Sciences
  • Mathematics
  • Statistics

Social Sciences

Featured Portals

  • Anthropology
  • Sustainability Science

Articles by Topic

  • Economic Sciences
  • Environmental Sciences
  • Political Sciences
  • Psychological and Cognitive Sciences
  • Social Sciences

Biological Sciences

Featured Portals

  • Sustainability Science

Articles by Topic

  • Agricultural Sciences
  • Anthropology
  • Applied Biological Sciences
  • Biochemistry
  • Biophysics and Computational Biology
  • Cell Biology
  • Developmental Biology
  • Ecology
  • Environmental Sciences
  • Evolution
  • Genetics
  • Immunology and Inflammation
  • Medical Sciences
  • Microbiology
  • Neuroscience
  • Pharmacology
  • Physiology
  • Plant Biology
  • Population Biology
  • Psychological and Cognitive Sciences
  • Sustainability Science
  • Systems Biology
Research Article

Power-law decay of the spatial correlation function in exciton-polariton condensates

Georgios Roumpos, Michael Lohse, Wolfgang H. Nitsche, Jonathan Keeling, Marzena Hanna Szymańska, Peter B. Littlewood, Andreas Löffler, Sven Höfling, Lukas Worschech, Alfred Forchel, and Yoshihisa Yamamoto
PNAS April 24, 2012 109 (17) 6467-6472; https://doi.org/10.1073/pnas.1107970109
Georgios Roumpos
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: georgios.roumpos@jila.colorado.edu
Michael Lohse
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Wolfgang H. Nitsche
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jonathan Keeling
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Marzena Hanna Szymańska
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Peter B. Littlewood
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Andreas Löffler
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sven Höfling
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Lukas Worschech
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Alfred Forchel
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Yoshihisa Yamamoto
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  1. Edited by Peter Zoller, University of Innsbruck, Austria, and approved February 27, 2012 (received for review May 20, 2011)

  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Abstract

We create a large exciton-polariton condensate and employ a Michelson interferometer setup to characterize the short- and long-distance behavior of the first order spatial correlation function. Our experimental results show distinct features of both the two-dimensional and nonequilibrium characters of the condensate. We find that the gaussian short-distance decay is followed by a power-law decay at longer distances, as expected for a two-dimensional condensate. The exponent of the power law is measured in the range 0.9–1.2, larger than is possible in equilibrium. We compare the experimental results to a theoretical model to understand the features required to observe a power law and to clarify the influence of external noise on spatial coherence in nonequilibrium phase transitions. Our results indicate that Berezinskii–Kosterlitz–Thouless-like phase order survives in open-dissipative systems.

  • quantum well excitons
  • semiconductor microcavities

The spatial correlation function quantifies the coherence properties of a system (1). In a 3D Bose-condensed gas, long range order is observed, and the correlation function decays toward a plateau at large distances (2, 3). In the homogeneous 2D Bose gas (4), however, no long range order can be established (5). Instead, Berezinskii–Kosterlitz–Thouless (BKT) theory of the equilibrium interacting gas predicts a transition to a low-temperature superfluid phase, which shows a power-law decay of the correlation function (6, 7). Unfortunately, it is frequently hard to directly measure this, and only very recently (8) was indication of the power-law decay of coherence seen in a 2D atomic gas. It has been theoretically predicted (9, 10) that power-law decay of coherence survives in the nonequilibrium problem, and it is this prediction that the current experiment sets out to test.

Exciton polaritons are short-lived quasiparticles formed in a semiconductor quantum well strongly coupled to a planar microcavity (11). Each one is a superposition of a quantum well exciton and a microcavity photon, and they behave as 2D bosons below the Mott density. Above a threshold particle density, condensation is observed (12). Due to the nonequilibrium nature of polariton condensation, understanding its coherence properties is quite revealing regarding the different roles of fluctuations in the equilibrium and nonequilibrium problems.

Previous measurements on polariton condensates have demonstrated coherence at large distances but were limited by large experimental uncertainties (13) or highly disordered samples (14, 15), and the long-distance behavior could not be fully extracted. Recently, the correlation function at large distances was studied in 1D condensates confined in a quantum wire (16) and in a valley of the disorder potential (17). In ref. 16, the data was energy-resolved so that excited states were filtered out, while in ref. 17 a rare area on the sample was chosen in which a single mode condensate is seen. The purpose of both those experiments was to investigate how long the coherence length of a spectrally isolated 1D condensate state can be. We, on the other hand, are interested in the functional form of the correlation function in a 2D condensate and how the excitations populated by the pumping and decay processes can modify it.

With our setup, we can measure values of g(1)(r) as low as 0.02, so we can reliably extract the long-distance behavior. We find that, although true thermal equilibrium is not established, an effective thermal de Broglie wavelength can still be defined from the short-distance gaussian decay of g(1)(r). Furthermore, g(1)(r) at long distances r decays according to a power law, in analogy to the equilibrium BKT superfluid phase. The exponent of the power-law decay is, however, higher than can be possible within the BKT theory. We apply a nonequilibrium theory (9, 10) to identify the source of the large exponent. We argue that, although the spectrum is modified due to dissipation, the exponent would still have the equilibrium value if the spectrum was thermally populated. If, on the other hand, a white noise source acts on the system and induces a flat occupation of the excited states, the exponent can have a large value, proportional to the noise strength. We therefore conclude that the pumping and decay processes, which introduce a nonthermal occupation of the excited states, can be responsible for the large value of the exponent.

Results

In our study, we use a weak-disorder GaAs-based sample, the same one as in our recent experiments (18). The condensate is generated nonresonantly by the multimode laser, which creates free electon-hole pairs at an excitation energy approximately 100 meV above the lower polariton (LP) energy. Carriers suffer multiple scatterings before reaching the LP energy, so coherence is established spontaneously in the condensate and cannot be inherited from the laser pump. The laser is continuously on and replaces LPs that leak out of the microcavity at a ps rate. We are interested in the limit of the homogeneous 2D polariton gas. For this purpose, we employ a setup based on a refractive beam shaper that forms a large laser excitation spot with uniform intensity. There is no confining potential on our sample. Because of the short lifetime, however, the condensate density follows the photon density of the excitation spot, so we can create circular condensates with almost flat density and diameters ranging from 14 μm to 44 μm (see ref. 18 and SI Appendix). LP luminescence in the steady state is observed through a combination of a long-pass and a band-pass interference filter, which reject scattered laser light without distorting the LP spectrum.

We confirmed that the sample disorder potential is weak in two ways (see SI Appendix). First, the lineshape of the luminescence at low excitation power is Lorentzian, which is characteristic of a homogeneously broadened line. Second, we measured a 2D map of the disorder potential with resolution approximately 1 μm and found that its spatial fluctuations are indeed weaker than the homogeneous broadening and also much weaker than the energy shift due to polariton–polariton interactions. Therefore, we can ignore the sample disorder in our experiment. The condensate is still localized in space, though, following the shape of the laser excitation spot.

The first order spatial correlation function is defined as Embedded Image[1]where Graphic and ψi are the creation and annihilation field operators at space-time point (ri,ti). To measure this function, we built a Michelson interferometer setup. A schematic is shown in Fig. 1A. It includes a mirror in one arm, and a right angle prism in the other. We overlap the condensate real-space image with its reflected version, so that fringes similar to that of Fig. 1B are observed on the camera. By changing the length of one interferometer arm, as shown in Fig. 1A, the relative phase of the two beams is shifted. As a result, the intensity measured at one pixel point shows a sinusoidal modulation (Fig. 1C). From the data of Fig. 1C, we extract the phase difference of the two images at a particular pixel point, as well as the fringe visibility. The latter is proportional to the first order correlation function, which is the physical quantity we are interested in in this experiment.

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

Michelson interferometer. (A) Schematic of the setup for measurement of the correlation function. The laser is linearly polarized, and we record luminescence of the orthogonal linear polarization through a polarizing beamsplitter (PBS). We then employ a 50-50 nonpolarizing beamsplitter (NPBS), a mirror (M1) and a right-angle prism (M2). The latter creates the reflection of the original image along one axis, depending on the prism orientation. A two-lens microscope setup overlaps the two real space images of the polariton condensate on the camera. (B) Typical interference pattern observed above the polariton condensation threshold along with a schematic showing the orientation of the two overlapping images. (C) Blue circles: measured intensity on one pixel of the camera as a function of the prism (M2) position in normalized units. Red line: fitting to a sine function.

The prism M2 in Fig. 1A forms the reflection of the condensate image along the prism axis. Therefore, point (x,y) overlaps with either (-x,y), or (x,-y) on the camera, depending on the orientation of the prism. This allows us to measure Embedded Image[2]or Embedded Image[3]where Graphic denotes time average. In this experiment, we are mainly interested in interference at τ = 0, so when the time argument is not mentioned explicitly, we imply τ = 0.

We repeat the procedure explained in Fig. 1 for every pixel, so that we measure the phase difference between the two interfering images in addition to the correlation function across the whole spot. Representative data are shown in Fig. 2. Recording both these quantities allows us to identify useful signal from systematic or random noise. Because the prism displaces the beam that is incident on it, the images from the mirror and the prism are focused on the camera from different angles, so the two phase fronts are tilted with respect to each other. As a consequence, we expect to measure a constant phase tilt. This is the case in Fig. 2B, in which the laser power is above threshold and a condensate has formed. We conclude that our measurement of the correlation function in Fig. 2D is reliable over this whole area. On the other hand, at a pump rate below threshold, only short-range correlations exist. Fig. 2A shows that in this case the phase difference is measured correctly only over a small area around the center, (|x| ≤ 1 μm). So the measured values of g(1)(x,-x) outside this area are not reliable and give an estimate of our measurement uncertainties. As is clear from Fig. 2C, the experimental error can be suppressed down to 0.01.

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

Phase map measured for laser power (A) below and (B) above the threshold power Pth. The prism in the Michelson interferometer is oriented horizontally. The schematics on the top right of A and B show the orientation of the two interfering images. (C and D) Measured g(1)(x,-x) corresponding to A and B, respectively, averaged over the y axis inside the excitation spot area of 19-μm radius. Blue circles are experimental data. The continuous red and dashed yellow fitting lines are explained in Figs. 3 and 6, respectively.

Fig. 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 6.

(A) g(1)(Δx) vs. Δx for increasing laser power. The laser pumping spot radius is R0 = 19 μm and the threshold power Pth = 55 mW. (B) g(1)(Δx) vs. Δx for one particular laser power and for x both positive (blue circles) and negative (red squares). Dashed line is a power-law fit. (C) Exponent ap of the power-law decay as a function of laser power.

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

(A) Decay of g(1)(x,-x) at short distances. Blue dots are experimental data, the red line is a gaussian fit. Data at |x| > 1 μm is noise. (B) Effective de Broglie wavelength λeff as a function of laser pumping power. λeff is extracted from the width of the gaussian fit as shown in A. Blue circles and red squares correspond to orthogonal orientations of the prism in the Michelson interferometer (see text). The condensation threshold is at approximately 55 mW.

Phase maps such as those in Fig. 2 have been used to identify localized phase defects—namely, quantum vortices (19). The data of Fig. 2B show that such localized defects are not present in our sample. At points with large fringe visibility (near x = 0 μm), fringes are perfectly parallel, whereas defects that appear for large |x| could be due to a numerical uncertainty in the measurement of the local phase due to the small fringe visibility. In any case, localized stationary phase defects cannot influence g(1)(r), because it is their motion that destroys spatial correlations and not their mere presence. It has been found that vortices appear in large disorder samples (19), when a direct external perturbation is introduced (20), before the condensate reaches its steady state (21), or when the condensate moves against an obstacle (22, 23). None of these conditions is satisfied in our experiment. On the other hand, we have found that, under the same conditions as the current experiment, mobile bound vortex pairs appear spontaneously due to the special form of the pumping spot and the pumping and decay noise (18). In ref. 18), we found that a single mobile bound vortex–antivortex pair is visible in a small condensate. In the current experiment, we probe larger condensate sizes, so it is likely that several vortex pairs are present at the same time. Mobile bound vortex pairs are in general invisible in time-integrated phase maps, like the one in Fig. 2B, and they are consistent with a power-law decay of g(1)(r).

Fig. 3A shows the short-distance dependence of g(1)(x,-x) for the same pumping power as in Fig. 2 A and C. Every dot in Fig. 3A corresponds to one pixel on the camera, and the x axis is its distance from the axis of reflection (slightly tilted with respect to the columns of the charge-coupled device array). Data at distances |x| > 1 μm is noise, because the measured phase in this area is random (Fig. 2A). At shorter distances, we can measure g(1)(x,-x) reliably, and we find that the correlation function has a gaussian form. This is the same functional dependence as for a thermalized Bose gas when the temperature is sufficiently high or the density sufficiently small (2, 4). In that equilibrium case, the width of the gaussian decay is proportional to the thermal de Broglie wavelength. Although our nonequilibrium system is quite different than the thermalized Bose gas, we will use this analogy to define a thermal de Broglie wavelength and therefore also a temperature. We note that the temperature measured from the short-distance behavior of g(1)(x,-x) is a measure of the occupation of the higher energy part of the spectrum (i.e., the particle-like part of the spectrum). For an insufficiently thermalized system, it is quite possible that excitations in different energy ranges have different effective temperatures. Therefore, the temperature measured this way will not necessarily agree with other measures of temperature.

In Fig. 3B we plot the effective wavelength λeff as a function of pumping power. If σ is the standard deviation of the gaussian fit for g(1)(x,-x), Graphic in analogy to the thermal de Broglie wavelength. λeff shows a smooth increase for increasing pumping power with no obvious threshold, analogous to the theory of equilibrium noninteracting 2D Bose gas as the particle density is increased (4). We performed the same experiment for two orthogonal prism orientations as shown in the legend of Fig. 3B. In one case we measured g(1)(x,-x), whereas in the other case we measured g(1)(y,-y). We found that λeff is shorter along the y axis and attribute this difference to a small asymmetry of the laser pumping spot. The occupation of excited states (which determines λeff) depends on their spatial overlap with the laser pumping spot, so states of equal energy are not always equally populated. This asymmetry shows that λeff is not simply related to the cryostat temperature and depends on the spatial and energy profiles of the high-energy states involved in producing this correlation length. We also note that the resolution limit of our imaging setup is approximately 1 μm, hence the measurement of λeff at small pumping power is resolution-limited.

It is known that an ideal autocorrelation measurement with a Michelson interferometer provides the same information as an ideal measurement of the spectrum. In particular, g(1)(x,-x; t) is the Fourier transform of the power spectrum in momentum space S(k,ω) (24). However, systematic noise in measurement of S(k,ω) currently makes the direct measurement of g(1)(x,-x; t) the only way to reliably extract λeff of Fig. 3B as well as the power-law decay at long distances to be explained later. The Fourier-transform relationship between g(1)(x,-x; t) and S(k,ω) is illustrated in Fig. 4. The measured g(1)(x,-x; t) at very low pumping power is shown in Fig. 4A. At time delay t = 0, it has a gaussian form as a function of x, but for increasing t it broadens and acquires a multipeak structure. This unusual space-time dependence is reproduced by the numerical Fourier transform (Fig. 4C) of measured S(k,ω) (Fig. 4B). As explained in SI Appendix, measurement of the time dependence of g(1)(x,-x; t) is limited by inhomogeneous broadening due to time-integrated data, so it cannot provide an estimate of the homogeneous dephasing time.

Fig. 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 4.

(A) Measured g(1)(x,-x; t) for very low pumping power. (B) Measured momentum-space spectrum S(kx,ℏω) for very low pumping power. As explained in the text, g(1)(x,-x; t) is the Fourier transform of S(kx,ℏω). (C) Fourier transform of the experimental data shown in B. The result indeed reproduces accurately A. In B and C, the data is plotted in linear color scale in arbitrary units.

At long distances, the behavior of the correlation function at zero time delay t = 0 is no longer gaussian. We found that it is influenced by the edge of the condensate. In Fig. 5, we plot the measured g(1)(Δx) = g(1)(|2x|) ≡ g(1)(x,-x) at pumping power P ∼ 3 × Pth for increasing pumping spot radius. The measured g(1)(Δx) at long distances decreases as the spot size is increased and eventually converges to a power-law decay for large condensates.

Fig. 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 5.

Measured g(1)(Δx) vs. Δx for various pumping spot radii R0. All data is taken above threshold and is chosen such that λeff ∼ 4.1 μm. As the condensate size increases, g(1)(Δx) converges to a power-law decay.

We note that the condensate size is slightly smaller than the pump laser spot radius (3–4 μm smaller from each side for a large spot). Because of the repulsive interaction between polaritons, and between polaritons and reservoir excitons, the large density of the condensate and reservoir creates an antitrapping potential that pushes LPs away from the center. This effect is stronger for a gaussian or a very small pumping spot and in long-lifetime samples (16, 25), whereas in the present experiment it only influences LPs that are close to the edge.

In the case of a large condensate, we should recover the limit of (infinitely large) homogeneous polariton gas. Therefore, we consider a pumping spot radius R0 = 19 μm. In Fig. 6A, we plot the correlation function g(1)(Δx) versus Δx as the pumping power is increased. Only short-range correlations exist for small pumping power, whereas above the condensation threshold of approximately 55 mW (4.8 kW/cm2), substantial phase coherence appears across the whole spot. The functional form of the long-distance decay is measured to be a power law over about one decade, as can be seen in Fig. 6B, in which we plot the data at one specific laser power. We fit the data to a function Graphic and plot the exponent ap as a function of pumping power in Fig. 6C. It is found to be in the range 0.9–1.2. λp is a parameter with units of length and is not related to λeff, which is plotted in Fig. 3B*.

It has been claimed that a criterion for polariton condensation is the appearance of a second threshold as the pumping power is increased (26, 27). The state after the first threshold was called a “polariton BEC,” “polariton laser,” or “polariton condensate,” whereas the state after the second threshold has not been fully understood yet. It might be a Bardeen–Cooper–Schrieffer (BCS) crossover (28, 29), photon BEC (30), or photon laser (26). This double threshold behavior has been observed in micropillar structures (26), and using a stress trap (27). In the supplementary information of ref. 18, we also reported the observation of double threshold using the same sample and excitation conditions as in the present experiment. We found that in our sample the window of intensities between the two thresholds is not very wide, and can only be witnessed using a flat laser excitation spot. This spot creates a uniform polariton density over a large area, as opposed to the more common gaussian spot, where the density changes a lot across the pumping spot.

Finally, we repeated the same measurement of g(1)(x,-x) using an identical sample at a temperature of 200 K. Because of the small binding energy, the GaAs excitonic effect is weak at this temperature. Also, the lasing energy was well above the bandgap. Therefore, only standard photon lasing was possible. In this case, we only found exponential decay of the correlation function and no power law. The details of this measurement are reported in SI Appendix. This suggests that the interactions of the strongly coupled exciton-polaritons are essential in the observation of the reported phenomena

Discussion

In ref. 31, it was found that excitation with a low-noise single mode laser revealed the formation of multimode condensation, and the different condensate modes could be spectrally separated. The authors of ref. 31 argued that, if one wants to measure the intrinsic linewidth of polariton condensates, single mode laser excitation and energy-resolved data are required. Indeed, it has been shown (32) that the temporal coherence properties can be understood based on the idea that laser intensity noise introduces population fluctuations, which modulate the interaction energy accordingly, leading to decoherence. However, it is not clear how pump and decay noise influences spatial coherence. Under single mode laser excitations, and when the lowest-energy state is spectrally isolated, long coherence lengths can be observed (15⇓–17, 33). In the current experiment, we are interested in how robust spatial correlations are when excitations are included. We study the “worst case” scenario of multimode laser excitation, which gives broader spectra compared to single mode laser. However, as shown in SI Appendix, single mode laser excitation gives similar results in energy-integrated data. We note that laser phase noise cannot be an issue in our experiment, because the laser energy is approximately 100 meV above the LP energy, so the generated quasiparticles suffer multiple scatterings before forming the condensate.

Focusing on the multimode laser excitation case, let us now explore the interpretation of the power law that we observe and consider what it means for the properties of the nonequilibrium polariton condensate. In particular, we discuss under what conditions a power-law decay should be seen and what may control the value and pump power dependence of the observed exponent. As has been discussed previously (9, 10, 34), power-law decay of spatial correlations are not an artifact of equilibrium condensates but survive more generally in a nonequilibrium condensate. Because the power-law decay at long distances arises from the long wavelength collective modes, this statement is not trivial, because dissipation can modify the spectrum at long wavelengths (9, 10, 34).

Let us first recall the results that would apply if one were to consider an equilibrium interacting 2D Bose gas. In this case, the exponent is given by ap = 1/nsλ2 ≤ 1/4, where ns is the superfluid density and λ is the thermal de Broglie wavelength Graphic. The restriction ap < 1/4 occurs because increasing temperature has two effects: It excites long wavelength phase fluctuations, which are responsible for the power-law decay, and it can also excite vortex pairs. The maximum value of ap occurs at the transition when vortex–antivortex pairs unbind, so that vortices would proliferate, and cause the BKT transition to a phase with short-range correlations. The observation here of a power law ap > 1/4 implies that effects beyond thermal equilibrium are required to explain the data; i.e. there is noise that excites phase fluctuations without leading to vortex proliferation. In addition, because the equilibrium exponent ap ∝ 1/ns, one would expect the exponent to decrease with pump power, as the condensate density increases; the absence of such a decrease again implies effects beyond thermal equilibrium are relevant and suggests that pumping noise is indeed affecting the observed exponent.

Whereas the existence of power-law decay in a nonequilibrium condensate was discussed previously, the value of the exponent and its pump power dependence were not given in those previous works. Using the formalism described in refs. 9 and 10, the exponent can be found by calculating Embedded Image[4]where Graphic, and Graphic are the Keldysh, retarded, and advanced Green’s functions for phase fluctuations. The advantage of writing the correlation function in this formal way is that it allows one to disentangle the effects of changes to the spectrum of long wavelength excitations from the effects of how this spectrum is populated.

The retarded and advanced Green’s functions are independent of how the spectrum is occupied, and following quite general arguments (9, 10, 34) we can prove that they have poles describing the low energy spectrum Graphic where ϵk = ℏ2k2/2m is the long wavelength polariton dispersion, μ the chemical potential (or blueshift) and γ is the linewidth.† Despite the modification of the equilibrium spectrum introduced by γ, it nevertheless remains the case that if this spectrum is occupied thermally (i.e., if the Keldysh Green’s function is chosen to obey the equilibrium fluctuation dissipation theorem), one finds † Embedded Image[5]which is independent of γ and matches the equilibrium form of f(k). Thus, despite the modifications to the long wavelength spectrum, a sufficiently thermalized polariton condensate has the equilibrium exponent.

In order to explain the larger exponent observed, and the flat dependence on pump power, we consider a crude model of a system with excess pumping noise as an opposite extreme to the thermalized case. We thus consider a case where the occupation of excitations is set by a Markovian noise source of strength ζ. Namely, we take the inverse Keldysh Green’s function to be energy independent.† This differs significantly from thermal noise correlations, which are frequency dependent, and diverge at the chemical potential. The measured spectra shown in SI Appendix are broad and their linewidth increases as the pumping power is increased. This occupation of excited states could be induced by an energy-independent noise source whose strength increases with the pumping power. In this case, the function f(k) is given by Embedded Image[6]which, despite the changed occupation spectrum, still yields a power-law decay. The exponent becomes ap = (mζ/2πℏ2ns)[(μ2 + γ2)/2μγ]. Because this has the form ap ∝ ζ/ns then if the noise strength and polariton density both increase with the pump power, then this would explain the absence of a 1/ns decrease of the exponent, as seen on Fig. 6C.

One point not addressed so far regards the process of vortex proliferation in a noisy nonequilibrium condensate. As pump noise increases, it is likely eventually to lead to proliferation of vortices, and a transition to a state with only short-range correlations, just as occurs at high temperatures in equilibrium.

In conclusion, the measured power-law decay of the correlation function suggests that some form of the BKT superfluid phase survives in nonequilibrium condensates; namely, phase flucutations are excited but no vortices. The large value of the exponent implies that, in the current experiment, this ordered phase is more robust against external noise than would be expected in equilibrium, in which equipartition holds. We conjecture that the main noise source is pump and decay noise, which create a nonthermal occupation of excited states, and apply a nonequilibrium theory to show that a power-law decay with a large exponent is possible in an open system with excess noise. One may anticipate that sufficient noise could induce vortex proliferation and a transition to short-range coherence. This fascinating possibility remains an open question for future studies.

Materials and Methods

Our GaAs-based sample shows a Rabi splitting of 2ℏΩRabi = 14 meV and LP lifetime of τLP ∼ 2–4 ps near photon-exciton detuning δ = 0, where the data presented in this paper is taken. From the curvature of the measured energy versus momentum dispersion at low pumping power, the LP effective mass was found to be m∗ = 9.5 × 10-5me at this detuning, where me is the electron rest mass. The sample is the same as in our recent experiments (18, 33), and the experimental setup is very similar to ref. 18. We pump the system with a multimode Ti-Sapphire laser operated in the continuous wave mode, combined with a chopper that creates 0.5-ms pulses at 100-Hz repetition rate. All powers quoted in the text and SI Appendix refer to the unchopped laser beam. We employ a commercial refractive beam shaper to generate a flat-top pumping profile of varying size. The Michelson interferometer consists of a 50-50 nonpolarizing cube beamsplitter, a dielectric mirror in the first arm, and an uncoated glass right angle prism in the second one. The position of the prism is controlled by a combination of a translation stage and a piezoelectric actuator.

Acknowledgments

We acknowledge support from National Science Foundation Grant ECCS-09 25549, Navy/SPAWAR Grant N66001-09-1-2024, FIRST, MEXT, Special Coordination Funds for Promoting Science and Technology, EPSRC, Department of Energy support under FWP 70069, and the State of Bavaria. G.R. thanks M. D. Fraser, F. M. Marchetti, I. Carusotto, and F. Amet for helpful discussions.

Footnotes

  • ↵1Present address: JILA, University of Colorado, Boulder, CO 80309

  • ↵2To whom correspondence should be addressed. E-mail: georgios.roumpos{at}jila.colorado.edu.
  • Author contributions: P.B.L. and Y.Y. designed research; G.R., M.L., W.H.N., J.K., and M.H.S. performed research; A.L., S.H., L.W., and A.F. contributed new reagents/analytic tools; G.R. and M.L. analyzed data; and G.R. and J.K. wrote the paper.

  • The authors declare no conflict of interest.

  • This article is a PNAS Direct Submission.

  • This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1107970109/-/DCSupplemental.

  • ↵*See SI Appendix for a discussion of λp, data at different detunings, the orthogonal prism orientation, as well as for time-resolved data.

  • ↵†See SI Appendix for further details.

References

  1. ↵
    1. Chaikin PM,
    2. Lubensky TC
    (2000) Principles of Condensed Matter Physics (Cambridge Univ Press).
  2. ↵
    1. Pitaevskii L,
    2. Stringari S
    (2003) Bose-Einstein Condensation (Oxford Univ Press).
  3. ↵
    1. Bloch I,
    2. Hänsch TW,
    3. Esslinger T
    (2000) Measurement of the spatial coherence of a trapped Bose gas at the phase transition. Nature 403:166–170.
    OpenUrlCrossRefPubMed
  4. ↵
    1. Hadzibabic Z,
    2. Dalibard J
    (2011) Two-dimensional bose fluids: An atomic physics perspective. Rivista del Nuovo Cimento 34:389.
    OpenUrl
  5. ↵
    1. Hohenberg PC
    (1967) Existence of long-range order in one and two dimensions. Phys Rev 158:383–386.
    OpenUrlCrossRef
  6. ↵
    1. Berezinskii VL
    (1972) Destruction of long-range order in one-dimensional and two-dimensional systems possessing a continuous symmetry group. ii. quantum systems. Sov Phys JETP-USSR 34:1144–1156.
    OpenUrl
  7. ↵
    1. Kosterlitz JM,
    2. Thouless DJ
    (1973) Ordering, metastability and phase transitions in two-dimensional systems. J Phys C Solid State 6:1181–1203.
    OpenUrlCrossRef
  8. ↵
    1. Hadzibabic Z,
    2. Krüger ,
    3. Cheneau M,
    4. Battelier B,
    5. Dalibard J
    (2006) Berezinskii-Kosterlitz-Thouless crossover in a trapped atomic gas. Nature 441:1118–1121.
    OpenUrlCrossRefPubMed
  9. ↵
    1. Szymańska MH,
    2. Keeling J,
    3. Littlewood PB
    (2006) Nonequilibrium quantum condensation in an incoherently pumped dissipative system. Phys Rev Lett 96:230602.
    OpenUrlCrossRefPubMed
  10. ↵
    1. Szymańska MH,
    2. Keeling J,
    3. Littlewood PB
    (2007) Mean-field theory and fluctuation spectrum of a pumped decaying Bose-Fermi system across the quantum condensation transition. Phys Rev B Condens Matter Mater Phys 75:195331.
    OpenUrlCrossRef
  11. ↵
    1. Weisbuch C,
    2. Nishioka M,
    3. Ishikawa A,
    4. Arakawa Y
    (1992) Observation of the coupled exciton-photon mode splitting in a semiconductor quantum microcavity. Phys Rev Lett 69:3314–3317.
    OpenUrlCrossRefPubMed
  12. ↵
    1. Deng H,
    2. Haug H,
    3. Yamamoto Y
    (2010) Exciton-polariton Bose-Einstein condensation. Rev Mod Phys 82:1489–1537.
    OpenUrlCrossRef
  13. ↵
    1. Deng H,
    2. Solomon GS,
    3. Hey R,
    4. Ploog KH,
    5. Yamamoto Y
    (2007) Spatial coherence of a polariton condensate. Phys Rev Lett 99:126403.
    OpenUrlCrossRefPubMed
  14. ↵
    1. Kasprzak J,
    2. et al.
    (2006) Bose-Einstein condensation of exciton-polaritons. Nature 443:409–414.
    OpenUrlCrossRefPubMed
  15. ↵
    1. Krizhanovskii DN,
    2. et al.
    (2009) Coexisting nonequilibrium condensates with long-range spatial coherence in semiconductor microcavities. Phys Rev B Condens Matter Mater Phys 80:045317.
    OpenUrlCrossRef
  16. ↵
    1. Wertz E,
    2. et al.
    (2010) Spontaneous formation and optical manipulation of extended polariton condensates. Nat Phys 6:860–864.
    OpenUrlCrossRef
  17. ↵
    1. Manni F,
    2. et al.
    (2011) Polariton condensation in a onedimensional disordered potential. Phys Rev Lett 106:176401.
    OpenUrlCrossRefPubMed
  18. ↵
    1. Roumpos G,
    2. et al.
    (2011) Single vortex-antivortex pair in an exciton polariton condensate. Nat Phys 7:129–133.
    OpenUrlCrossRef
  19. ↵
    1. Lagoudakis KG,
    2. et al.
    (2008) Quantized vortices in an exciton-polariton condensate. Nat Phys 4:706–710.
    OpenUrlCrossRef
  20. ↵
    1. Tosi G,
    2. et al.
    (2011) Onset and dynamics of vortex-antivortex pairs in polariton optical parametric oscillator superfluids. Phys Rev Lett 107:036401.
    OpenUrlCrossRefPubMed
  21. ↵
    1. Lagoudakis KG,
    2. et al.
    (2011) Probing the dynamics of spontaneous quantum vortices in polariton superfluids. Phys Rev Lett 106:115301.
    OpenUrlCrossRefPubMed
  22. ↵
    1. Amo A,
    2. et al.
    (2011) Polariton superfluids reveal quantum hydrodynamic solitons. Science 332:1167–1170.
    OpenUrlAbstract/FREE Full Text
  23. ↵
    1. Nardin G,
    2. et al.
    (2011) Hydrodynamic nucleation of quantized vortex pairs in a polariton quantum fluid. Nat Phys 7:635–641.
    OpenUrlCrossRef
  24. ↵
    1. Weiner AM
    (2009) Ultrafast Optics (John Wiley & Sons, Hoboken, NJ).
  25. ↵
    1. Wouters M,
    2. Carusotto I,
    3. Ciuti C
    (2008) Spatial and spectral shape of inhomogeneous nonequilibrium exciton-polariton condensates. Phys Rev B Condens Matter Mater Phys 77:115340.
    OpenUrlCrossRef
  26. ↵
    1. Bajoni D,
    2. et al.
    (2008) Polariton laser using single micropillar GaAs-GaAlAs semiconductor cavities. Phys Rev Lett 100:047401.
    OpenUrlCrossRefPubMed
  27. ↵
    1. Balili R,
    2. Nelsen B,
    3. Snoke DW,
    4. Snoke L,
    5. West K
    (2009) Role of the stress trap in the polariton quasiequilibrium condensation in gaas microcavities. Phys Rev B Condens Matter Mater Phys 79:075319.
    OpenUrlCrossRef
  28. ↵
    1. Keeling J,
    2. Eastham PR,
    3. Szymanska MH,
    4. Littlewood PB
    (2005) BCS-BEC crossover in a system of microcavity polaritons. Phys Rev B Condens Matter Mater Phys 72:115320.
    OpenUrlCrossRef
  29. ↵
    1. Byrnes T,
    2. Horikiri T,
    3. Ishida N,
    4. Yamamoto Y
    (2010) BCS wave-function approach to the BEC-BCS crossover of exciton-polariton condensates. Phys Rev Lett 105:186402.
    OpenUrlCrossRefPubMed
  30. ↵
    1. Klaers J,
    2. Schmitt J,
    3. Vewinger F,
    4. Weitz M
    (2010) Bose-Einstein condensation of photons in an optical microcavity. Nature 468:545–548.
    OpenUrlCrossRefPubMed
  31. ↵
    1. Love APD,
    2. et al.
    (2008) Intrinsic decoherence mechanisms in the microcavity polariton condensate. Phys Rev Lett 101:067404.
    OpenUrlCrossRefPubMed
  32. ↵
    1. Whittaker DM,
    2. Eastham PR
    (2009) Coherence properties of the microcavity polariton condensate. EPL 87:27002.
    OpenUrlCrossRef
  33. ↵
    1. Roumpos G,
    2. Nitsche WH,
    3. Höfling S,
    4. Forchel A,
    5. Yamamoto Y
    (2010) Gaininduced trapping of microcavity exciton polariton condensates. Phys Rev Lett 104:126403.
    OpenUrlCrossRefPubMed
  34. ↵
    1. Wouters M,
    2. Carusotto I
    (2006) Absence of long-range coherence in the parametric emission of photonic wires. Phys Rev B Condens Matter Mater Phys 74:245316.
    OpenUrlCrossRef
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Power-law decay of the spatial correlation function in exciton-polariton condensates
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Power-law decay of g1(r) in polariton condensates
Georgios Roumpos, Michael Lohse, Wolfgang H. Nitsche, Jonathan Keeling, Marzena Hanna Szymańska, Peter B. Littlewood, Andreas Löffler, Sven Höfling, Lukas Worschech, Alfred Forchel, Yoshihisa Yamamoto
Proceedings of the National Academy of Sciences Apr 2012, 109 (17) 6467-6472; DOI: 10.1073/pnas.1107970109

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
Power-law decay of g1(r) in polariton condensates
Georgios Roumpos, Michael Lohse, Wolfgang H. Nitsche, Jonathan Keeling, Marzena Hanna Szymańska, Peter B. Littlewood, Andreas Löffler, Sven Höfling, Lukas Worschech, Alfred Forchel, Yoshihisa Yamamoto
Proceedings of the National Academy of Sciences Apr 2012, 109 (17) 6467-6472; DOI: 10.1073/pnas.1107970109
Digg logo Reddit logo Twitter logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley
Proceedings of the National Academy of Sciences: 109 (17)
Table of Contents

Submit

Sign up for Article Alerts

Article Classifications

  • Physical Sciences
  • Physics

Jump to section

  • Article
    • Abstract
    • Results
    • Discussion
    • Materials and Methods
    • Acknowledgments
    • Footnotes
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

Surgeons hands during surgery
Inner Workings: Advances in infectious disease treatment promise to expand the pool of donor organs
Despite myriad challenges, clinicians see room for progress.
Image credit: Shutterstock/David Tadevosian.
Setting sun over a sun-baked dirt landscape
Core Concept: Popular integrated assessment climate policy models have key caveats
Better explicating the strengths and shortcomings of these models will help refine projections and improve transparency in the years ahead.
Image credit: Witsawat.S.
Double helix
Journal Club: Noncoding DNA shown to underlie function, cause limb malformations
Using CRISPR, researchers showed that a region some used to label “junk DNA” has a major role in a rare genetic disorder.
Image credit: Nathan Devery.
Steamboat Geyser eruption.
Eruption of Steamboat Geyser
Mara Reed and Michael Manga explore why Yellowstone's Steamboat Geyser resumed erupting in 2018.
Listen
Past PodcastsSubscribe
Birds nestling on tree branches
Parent–offspring conflict in songbird fledging
Some songbird parents might improve their own fitness by manipulating their offspring into leaving the nest early, at the cost of fledgling survival, a study finds.
Image credit: Gil Eckrich (photographer).

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Special Feature Articles – Most Recent
  • List of Issues

PNAS Portals

  • Anthropology
  • Chemistry
  • Classics
  • Front Matter
  • Physics
  • Sustainability Science
  • Teaching Resources

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Librarians
  • Press
  • Site Map
  • PNAS Updates

Feedback    Privacy/Legal

Copyright © 2021 National Academy of Sciences. Online ISSN 1091-6490