since the prediction of keldysh and kopaev ,  @xcite there has been a long and intense pursuit to realise a condensed phase in solid state excitonic , and related , systems .
in particular , polaritons in semiconductor microcavities , the coupled eigenstates of an exciton with a cavity photon ,  @xcite represent ideal candidates for observing condensation phenomena .
the very light mass of these composite bosonic particles promises relatively high transition temperatures . in the past decade , improvements in the growth technology of semiconductor heterostructures have made the study of high - quality strongly - coupled planar microcavities almost routine for iii - v and ii - vi semiconductors .
the high degree of external control of these systems , and the possibility of direct detection for them has opened the route towards a new generation of fast optical matter - wave lasers and amplifiers .
@xcite more recently , concerted experimental efforts have been devoted to the realisation of a bose - einstein condensate of microcavity polaritons .
@xcite    on the experimental side , a challenge to the realisation of a condensed polariton phase might be represented by the finite quality of the cavity mirrors and the resultant short polariton lifetime , of the order of picoseconds .
in addition , due to the ` bottleneck effect ' ,  @xcite the relaxation of polaritons to the zero momentum state can be delayed , hindering the creation of a thermal population in the lowest energy state .
it has however been recently shown  @xcite that thermalization processes due to particle - particle scattering can be dramatically magnified by increasing the value of the ( non - resonant ) pump power , and by positively detuning the cavity energy above the excitonic energy . under these conditions ,
the progress towards a zero momentum quasi - equilibrium condensate has been significant ,  @xcite including a non - linear threshold behaviour in the emission intensity at zero momentum ,  @xcite the investigation of the second order coherence function ,  @xcite a characteristic change in the momentum space distribution above threshold ,  @xcite and evidence that the equilibration time is shorter than the polariton lifetime has been seen .
@xcite finally , very recently , a clear demonstration of condensation of cavity polaritons has been demonstrated in cdte .
@xcite kasprzak and collaborators have shown that condensation of equilibrated polaritons can be achieved for effective temperatures around @xmath0k and evidence for condensation has been seen in the occupation function , the first order coherence ( both in time and in space ) and in the spontaneous appearance of linear polarisation of the condensate emission .    alongside the experimental effort ,
a significant theoretical effort has been invested in analysing properties and predicting signatures of polariton condensation .
@xcite much of this work focuses on modelling the conditions under which condensation can occur , both in equilibrium ,  @xcite and considering the effects of pumping and decay .
@xcite possible signatures include the nonlinear relation of emission at zero momentum to pumping power ,  @xcite changes to the linewidth ,  @xcite the pl spectrum and the angular distribution of radiation ,  @xcite and spontaneous polarisation of emitted radiation .
@xcite    in this paper , we discuss the optical properties of condensed polaritons , focusing our interests on absorption , photoluminescence ( pl ) and resonant rayleigh scattering ( rrs ) . from our study
we conclude that rrs , the coherent scattering by disorder of polaritons into directions other than that of the original probe , represents a powerful tool for investigating the condensed phase .
we will show that signatures of condensation are visible in the rrs spectrum , allowing a direct probe of the collective excitation properties of the polariton condensate .
in particular , we will show that , above the threshold for condensation , strong emission from the collective sound mode both above and below the chemical potential can be seen in the rrs spectra ( see fig .
[ fig : rrsca ] ) .
in contrast such features are expected to be much harder to observe in usual pl emission spectra , where the spectrum is dominated by the very strong condensate emission at the chemical potential , which is likely to mask these more subtle features .
in addition , we will show that the rrs spectra directly reflects the disorder averaged density of states of excitonic particle - hole like excitations , i.e. bound excitons coupled to the coherent photon field .     for @xmath1 as a function of the dimensionless momentum @xmath2 and rescaled energy @xmath3 , for zero detuning , rabi splitting @xmath4mev , temperature @xmath5k , and a disorder strength characterize by an inverse scattering time @xmath6mev .
( a ) non - condensed regime ( dimensionless density @xmath7 ) ; ( b ) condensed regime ( @xmath8 ) ; ( c ) condensed regime ( @xmath9 ) . ( the parameters chosen for these plots are the same as those used later for spectral weight and photoluminescence . ) the value of the chemical potential is explicitly marked [ horizontal green ( gray ) line ] . while in the non - condensed regime ,
rrs emission is always above the chemical potential , in the condensed phase , emission from the collective sound mode is seen both above and below the chemical potential . ]
resonant rayleigh scattering depends on disorder to scatter polaritons between momentum states . in order to carry on our analysis
, we introduce a realistic description of disorder and analyse its effects on rrs and more generally on other optical probes , like absorption and pl .
to do this , we make use of a quantitatively accurate model for exciton disorder ,  @xcite and numerically evaluate the distribution of excitonic energies and oscillator strengths associated with a given disorder potential .
an accurate treatment of exciton disorder for quantum wells in microcavities is important because the large ratio of exciton mass to photon mass means that those exciton states contributing most to the thermally populated polaritons are strongly localised , i.e. can not be treated within the first born approximation .
such a treatment shows that , for a given exciton energy and momentum , there is still a distribution of possible oscillator strengths . for the exciton states relevant
to polariton formation , the distribution of oscillator strengths varies from a narrow gaussian at low energies ( below the band edge ) to a porter - thomas distribution at high energies , as is known from previous works .
@xcite here , we study the effect of this distribution on the many - body physics , and in particular , on rrs and pl spectra .    in this paper
, we consider the case relevant to the existing cdte or gaas microcavities of high - quality quantum wells , where the typical excitonic disorder amplitude is smaller than the rabi splitting . in this regime
, we can make use of the _ coupled oscillator model _ or wave vector conserving approximation , as explained in ref .  .
as will be explained in more detail later on , this corresponds to approximating the full exciton green s function by its momentum - diagonal part ( or equivalently considering its disorder average ) , and treating perturbatively the off - diagonal terms . in this way
, the translational invariance broken by the presence of the disorder is restored , polaritons have a defined wave vector and , at high enough densities , condense in the lowest momentum state .
such a treatment implies that excitonic disorder , being on short length scales , does not lead to spatial inhomogeneity on the length scale associated with the polariton .
extended polaritons are formed from a superposition of many localised excitons , therefore recovering translational invariance at the level of polaritons . at low densities , in the non - condensed state ,
many observable properties can be adequately described by the coupled oscillator model and so can be found from the mean squared oscillator strength at a given energy . in this
limit our method recovers the well known results for the inhomogeneous broadening of the polariton pl ,  @xcite and  by considering in addition the mean fourth power of the oscillator strength  the averaged rrs response .
@xcite however , when condensed , there are observable effects associated with the full distribution of oscillator strengths . in particular , the non - vanishing probability of excitons to have arbitrarily small oscillator strengths
has direct consequences for optical probes , including both rrs and pl .    in ref .  , we considered specific aspects of resonant rayleigh scattering arising from the model discussed in this paper . here ,
we provide and compare further experimental probes of condensation , and discuss the underlying physical mechanisms involved .
all these optical probes , including rrs , are in addition to a non - resonant pumping of the microcavity polaritons .
a related problem is studied in ref .  . there ,
the rayleigh scattering of a strong resonant pump is considered ( treating disorder perturbatively ) , thus the probe is the pump . in this paper , in contrast , one has to distinguish two types of coherence .
the first type is the internal coherence of the condensed system , which arises spontaneously following non - resonant pumping .
the second type is the coherent scattering of an external laser probe on disorder , which can be strongly modified by the presence of the condensate .
the form of the density of states , optical density , and distribution of oscillator strengths found from the numerical calculations are also important in other thermodynamic properties and probes of polaritons .
one such example is the calculation of the linear blue - shift vs. density in the low density regime .
the inclusion of exciton states that couple weakly to photons extends the validity of the current model to densities beyond those in similar treatments ,  @xcite and provides a stronger basis for the use of those models in the regimes where they are valid .
the paper is arranged in the following sections : we describe the effect of quantum well disorder on the exciton oscillator strength in sec .
[ sec : disor ] .
there , we also give details of the numerical analysis ( sec .
[ sec : numer ] ) , which will be used in the many - body model introduced in sec .
[ sec : model ] to describe the thermodynamics of polariton condensation ( sec .
[ sec : phase ] ) .
section  [ sec : photo ] describes the optical probes such as spectral weight and photoluminescence ( sec .
[ sec : spect ] ) and resonant rayleigh scattering ( sec .
[ sec : rrsca ] ) .
conclusions are collected in sec .
[ sec : concl ] , while appendix  [ sec : phase - sens - detect ] explains how one can detect rrs using phase sensitive measurements .
the problem of an exciton in a disordered quantum well has been studied at some length in the last two decades .
@xcite quantum well disorder can arise due to interface and alloy fluctuations and affects the properties of the excitonic linewidth and absorption spectrum .
similarly to refs .  , we will assume the external disorder potential to be correlated on a length scale larger than the exciton bohr radius .
accordingly , we factorise the excitonic in - plane relative and centre of mass coordinates @xmath10 and assume that the disorder affects only the excitonic centre of mass motion @xmath11 , while the internal degrees of freedom can be restricted to the ground state hydrogenic state @xmath12 . here , @xmath13 is the exciton bohr radius and @xmath14 the reduced mass ( henceforth we will set @xmath15 ) . neglecting the transverse degrees of freedom related to the confinement of the excitons in the qw
, the energy associated to the wavefunction @xmath16 is given by the sum of the relative motion energy @xmath17 [ i.e. , the band gap minus the exciton binding energy , @xmath18 and the energy related to the centre of mass motion @xmath19 : @xmath20 \phi_{\alpha } ( { { \mathbf r } } ) = \varepsilon_{\alpha }    \phi_{\alpha } ( { { \mathbf r}})\ ; . \label{eq : centr}\ ] ] here , the effective disorder potential @xmath21 represents the microscopic structural disorder averaged over the electron - hole motion .
@xcite this can be approximated , e.g. , with a gaussian noise correlated on a length scale @xmath22 with variance equal to @xmath23 : @xmath24 where @xmath25 is the quantisation area .
it will be convenient , later on , to introduce the scattering time @xmath26 where @xmath27 is the two - dimensional density of states in the clean limit .    in two dimensional non - interacting systems
, all states are localised by the disorder potential .
however , the localisation length and the character of the excitonic wavefunction change significantly from below to above the band edge @xmath17 .
well below the band edge , low - energy lifshitz tail states ,  @xcite are well localised in deep potential minima , with a nodeless ( roughly gaussian ) shape .
these states are rare , because they occur due to large , rare , fluctuations of the disorder potential .
in contrast , higher energy states above the band edge have a fractal - like shape with many nodes and can be approximated by a random superposition of plane waves with the same momentum @xmath28 . here
, the localisation mechanism is closely related to quantum mechanical interference effects .
accordingly , the change of the shape of the centre of mass wavefunctions across the band edge is accompanied by a dramatic change in the excitonic oscillator strength .
@xcite the oscillator strength describing coupling of a quantum well exciton to light , @xmath29 , is given by the probability amplitude of finding an electron and a hole at the same position and with centre of mass momentum equal to the photon momentum @xmath30 and is therefore proportional to the fourier transform of the centre of mass wave function , @xmath31 : @xmath32 where , @xmath33 is the dipole matrix element .
the dispersion for photons in a microcavity of width @xmath34 is given by @xmath35 , where @xmath36 , and can be approximated , for small momenta , by a parabolic dispersion @xmath37 , with the photon mass given by @xmath38 .
because there may be many different exciton wavefunctions corresponding to similar exciton energies , the oscillator strength @xmath39 is a random quantity , which varies both in phase and magnitude .
considering many disorder realisations , we find a distribution of squared oscillator strengths , @xmath40 as a function of the energy @xmath19 and momentum @xmath30 .
this distribution reflects the statistical properties of the centre of mass excitonic wave functions .    to compare to experimentally relevant observables
, we introduce the density of states ( dos ) , @xmath41 and the mean squared oscillator strength : @xmath42 where @xmath43 is the average over different disorder realisations .
these quantities are related to the excitonic optical density by the relation : @xmath44 experimentally , the excitonic optical density can be measured by dividing the pl emission by the excitonic occupation .
@xcite    focusing for the moment on average properties , rather than on the entire distribution , a simple expression exists for the mean squared oscillator strength in the high energy limit . in this limit ,
where the dos is flat and energy independent , @xmath45 , making use of the born approximation , the squared average oscillator strength @xmath46 does not depend separately on the disorder potential correlation length @xmath47 and variance @xmath48 , but instead only on the scattering time @xmath49 ( see fig .
[ fig : gepsp ] ) : @xmath50 where @xmath51 is the free particle dispersion . from this form , and the effectively constant density of states
, one can see that both the mean squared oscillator strength and the excitonic optical density are symmetrical . the comparison of this approximate form to the numerical simulation
is shown in fig .  [
fig : gepsp ] ,    in contrast , for energies much below the band edge , the specific asymptotic expression of @xmath46 depends whether the correlation length @xmath47 is smaller or larger than the localisation length @xmath52 . in the white noise
limit,@xcite @xmath53 , one can show that the centre of mass wavefunction @xmath11 can be approximated by a gaussian centred at a randomly distributed site @xmath54 , @xmath55 and therefore giving a squared oscillator strength proportional to : @xmath56 thus , here the distribution of squared oscillator strengths is very narrow , with a mean square value given by the above form .
similarly , much theoretical ( and numerical ) work has been done to establish the energy dependence of the density of states . in the low energy tail , as before , the specific asymptotic form of the dos depends on the value of the correlation length @xmath47 . in the white noise limit , one can show that @xmath57 , while , in the opposite ( classical ) limit , @xmath58 , instead one has @xmath59 ( see , e.g. , ref .  ) .
in general , for the finite values of the disorder correlation length corresponding to typical experiments , the regions in energy where one of these two analytical regimes apply are very restricted , and therefore a numerical analysis is required .
numerical analysis is also essential in order to account for the distribution of squared oscillator strengths @xmath60 near the band edge .
this distribution changes substantially from low energy states to high energy ones .
as we discuss more in detail in the next section , for high energy states the oscillator strength distribution is governed by a porter - thomas law , while for lifshitz tail states the distribution follows a narrow gaussian - like distribution .
neither of these distributions apply for energies around the band edge .
however , it is states near the band edge that have the largest optical density , and so the distribution of oscillator strengths in this region has a significant impact on derived quantities .
it thus becomes essential to use numerical analysis to find the entire distribution of oscillator strengths .
to photons of momentum @xmath61 ( a ) ( from ref .  , for comparison ) and @xmath62 ( b ) , where all exciton states found numerically .
results are taken from 160 different realizations of disorder potential , and for ( b ) , coupling strengths from eight different photon momenta @xmath63 with the same value of @xmath64 are combined .
the mean squared averaged oscillator strength @xmath65 for the same value of momentum [ lower red ( gray ) points ] and the density of states @xmath66 [ upper green ( gray ) points ] are also explicitly plotted .
inset : fit of @xmath67 to the expression   with a renormalized energy @xmath68 . ]     to photons of momentum @xmath61 ( a ) ( from ref .
, for comparison ) and @xmath62 ( b ) , where all exciton states found numerically .
results are taken from 160 different realizations of disorder potential , and for ( b ) , coupling strengths from eight different photon momenta @xmath63 with the same value of @xmath64 are combined .
the mean squared averaged oscillator strength @xmath65 for the same value of momentum [ lower red ( gray ) points ] and the density of states @xmath66 [ upper green ( gray ) points ] are also explicitly plotted .
inset : fit of @xmath67 to the expression   with a renormalized energy @xmath68 . ]     versus energy and momentum [ or equivalently angle @xmath69 .
note that the scale in angle is not linear . the free particle dispersion @xmath70 [ solid green ( gray ) line ] , and the a trace of the mean squared oscillator strength for two representative values of momenta , @xmath71 ( @xmath72 ) and @xmath73 ( @xmath74 ) [ red ( gray ) plus symbols ] are explicitly plotted ( cf fig .
[ fig : gepsp ] ) . the figure is adapted from ref .  .
]      to solve eq .
, we exactly diagonalise this problem within a finite basis set , using conjugate - gradient minimisation with pre - conditioning of the steepest descent wavevector ( for a detailed discussion of this method , see ref .  ) .
we find the exact eigenvalues and eigenfunctions in a finite basis on a grid of @xmath75 points for a system of size @xmath76 m , @xmath77mev , @xmath78  and excitonic mass @xmath79 .
these parameters give an inverse scattering time of @xmath80mev . for this choice of the grid
, one can show that convergence is reached . from the evaluated eigenvalues @xmath81 and eigenstates @xmath82 over 160 realisations of the disorder potential , we can derive the excitonic density of states , the oscillator coupling strength and its squared average , eq .  , which we plot in figs .
[ fig : gepsp ] , and  [ fig : geps2 ] , while the corresponding optical density is plotted in fig .
[ fig : optic ] .    the lower panel of fig .
[ fig : gepsp ] shows the squared coupling strength @xmath60 versus energy for a fixed value of momentum , @xmath83 , corresponding , for a cavity of @xmath84ev , to an angle of @xmath85 .
note that , because of the presence of disorder in the quantum well , one photon with a given momentum couples with many exciton states with different energies .
we will see later on therefore that a polariton with a given momentum is formed by the superposition of one photon state @xmath86 and many exciton states @xmath87 .
these states are more or less strongly coupled depending on the distribution of oscillator strength for that given momentum . as fig .  [
fig : gepsp ] shows , by probing a quantum well at a large angle , i.e. with high momenta photons , the excited excitons with larger oscillator strength are the ones which are almost - delocalised in nature and with a many - node fractal - like shape . by taking the average over many ( 160 ) disorder realisations , the squared average oscillator strength
@xmath65 shown in the inset in fig .
[ fig : gepsp ] is well described by the lorentzian shape predicted by the born approximation of eq .  , with a fitted width of @xmath88mev in good agreement with the theoretical value @xmath80mev .
however , the peak of the lorentzian does not coincide with the energy of the clean limit , @xmath68 , but is renormalised down in value , as can be shown by employing a self - consistent born approximation .    in contrast , as shown in the upper panel of fig .
[ fig : gepsp ] , for photons with zero momentum , the maximum value of the oscillator strength characterises excitonic states below the band edge which are more localised in nature .
this can be easily understood by the following qualitative argument : at very low energies , the excitonic state is strongly localised in a deep potential minimum and has no nodes . increasing the energy , at first the localisation length increases [ e.g. , @xmath89 in the white noise limit ] and thus increases the oscillator strength . however , eventually , the wavefunction starts developing nodes and consequently the squared average oscillator strength decreases .
when @xmath90 , only the high energy side of @xmath91 can be described by the born limit  . for the chosen values of the disorder potential correlation length @xmath47 and variance @xmath48 , an analytical expression for the low energy ( lifshitz ) tail is not known .
the crossover from localised to more plane - wave - like excitonic states obtained by increasing the value of the photon momentum is plotted in fig .
[ fig : geps2 ] . at high momenta ,
the maximum of the squared average oscillator strength @xmath65 shown in fig .
[ fig : geps2 ] follows the free particle dispersion , plus a renormalisation down in energy , which decreases for higher photon momenta , i.e. these states are described well by free particles , including disorder in the first born approximation . at low momenta , the states which have the stronger oscillator strength
are effectively localised .
the crossover , as seen in the contour plot , happens at relatively large angle , @xmath92 , because of the large ratio of exciton to photon mass : the crossover momentum is set by exciton mass , but its conversion to an angle depends on effective photon mass . for a connected reason , the thermally populated polariton states are formed out of strongly localised ( i.e. beyond first born approximation ) excitonic states : for values of temperature and photon mass relevant for experiments , thermal population of polaritons extends up to around @xmath93 in momentum , which as seen in fig .
[ fig : geps2 ] corresponds to exciton states not accurately described by the born approximation .
for this reason , in the following we will concentrate on the oscillator strength corresponding to @xmath94 , @xmath95 .
mev and the fwhm @xmath96mev . ]    finally we plot in fig .
[ fig : optic ] the optical density . for the chosen values of @xmath47 and @xmath48
, the optical density shows a maximum around @xmath97mev below the band edge , a full width at half maximum ( fwhm ) of approximatively @xmath96mev and a clear asymmetry of the line - shape .
it is well known  @xcite that , by indicating with @xmath98 the confinement energy of the lowest state in a typical potential minimum , the excitonic line shape is determined by the ratio @xmath99 and , for a finite value of this , the optical density develops an asymmetry towards higher energies . for our choices of parameters , @xmath100 .
asymmetry of the optical density of quantum well excitons has been also measured experimentally ( see , e.g. , ref .  ) , by dividing the measured pl by the excitonic occupation .    as mentioned earlier ,
the full distribution of oscillator strengths , and not just its mean squared value will be important .
it is useful to discuss here some technical details of how this is extracted .
the numerical analysis provides the excitonic eigenvalues and eigenstates only within a finite interval in energy .
( the lowest energy states arise from rare potential fluctuations , which would require a larger region of space to be sampled ; the highest energy states have spatial variation on lengthscales finer than our grid . ) within the interval of energies found , averages are performed making use of the raw data coming from 160 realisations of the disorder potential ; outside this interval , averages are taken by extrapolating the numerics . in particular , in the low energy lifshitz tail , as the oscillator strengths have a narrow gaussian distribution around its squared averaged value , we approximate the distribution of @xmath95 with a delta function at its extrapolated value @xmath101 . in the very high energy region , instead , we make use of the porter - thomas distribution , @xmath102}{\sqrt{2 \pi x \bar{x } } } } \ ; , \ ] ] where @xmath103 , again extrapolating the fitted value for the squared averaged oscillator strength .
in addition , we fix the overall scale of @xmath95 to match integrated optical intensity with half the experimentally measured rabi splitting @xmath104 squared : @xmath105 this normalisation accounts for factors other than the wavefunction in eq .
( [ eq : oscil ] ) .    having described the effect of quantum well disorder on the energies of excitons and on their dipole coupling to the cavity photons , we will next introduce the many - body hamiltonian that will allow us to describe a system of thermalized interacting polaritons which , for a given density and temperature , can condense to a superfluid phase .
the hamiltonian describing a model of excitons with random energies @xmath19 , corresponding to the centre of mass eigenstate from eq .
, dipole coupled to a cavity mode can be written as @xmath106 the description of exciton - exciton interaction in this model is somewhat simplified : it imposes the constraint that two excitons can not occupy the same energy level @xmath19 , because of their composite fermionic nature .
this constraint is implemented by representing each energy level [ where energy levels are distributed according to the density of states , @xmath66 ] by two possible states , respectively @xmath108 and @xmath109 .
these two states are then represented by occupation of two fermionic states , so that the ground state @xmath110 , and an excitonic state @xmath111 , eliminates the unphysical states @xmath112 and @xmath113 . in this way
saturation effects are introduced to all orders .
note that these fermionic operators @xmath114 can not be written as a linear combination of free electron and hole operators  this is a result of having already included the intra - exciton coulomb interaction .
the transition matrix element between a photon state @xmath115 and an exciton state , @xmath116 , is given by @xmath29 from eq .   and evaluated numerically , as explained in the previous section .
it is convenient to rescale this coupling according to @xmath117 where @xmath118 is the inverse level spacing @xmath119 , measured in units of the excitonic rydberg energy .
this corresponds to measuring the density of particles per bohr radius squared . using these units , we may write the total number of excitations @xmath120 n a dimensionless form by introducing a dimensionless density of particles @xmath121 , or equivalently @xmath122 , where @xmath123 is the physical areal density of particles .     and thermal occupation factor [ red ( gray ) line starting at top left ] for @xmath124mev and @xmath125k .
]    in the limit of vanishing density , our occupancy constraint has no effect , and so this description is equivalent to one treating excitons as bosons . as will be shown later , in sec .
[ sec : blues ] , for low densities , our occupancy constraint is also equivalent to a model of excitons as bosons , but with some effective saturation interaction .
thus , as in the bosonic case , saturation effects are linear  @xcite and the blue - shift increases linearly in density .
there is however an important difference to those previous works .
in contrast to the clean systems studied there , our effective saturation interaction has a coefficient related to disorder and temperature rather than to the rabi splitting and bohr radius only .
however , increasing the density , saturation effects , and thus the underlying fermionic degrees of freedom become more and more important , and the conventional linear treatment becomes inadequate . in our treatment saturation effects are introduced exactly and taken into account to all orders .
our model is accurate as far as the occupancy is restricted to the strongly localised lifshitz states in the tail below the band edge ( see fig .  [
fig : occup ] ) .
all the results described here respect this limit .
one can moreover show that it is this density regime which is relevant for on - going experiments in cdte .
@xcite higher energy states , beyond the lifshitz tail , contribute to the optical response , however they are only weakly occupied , and so the above limit remains valid .    what the model eq .
does not include is the coulomb interaction between excitons with different energies @xmath19 , i.e. , at low densities , on different localisation sites @xmath54 .
this contribution is expected to be small in the low density regime .
this low density regime will be observed all the way through this paper . at the same time
we are not including double - occupancy of a single exciton energy level , which could be important at higher densities .    in the following ,
we will make extensive use of the numerical analysis of oscillator strength distributions explained in the previous section .
let us note here again that in the non - condensed regime , the full distribution is unnecessary  the excitonic optical density is sufficient to derive the polariton dispersion and its inhomogeneous broadening .
however , in the condensed phase , and in order to correctly describe resonant rayleigh scattering mechanisms , we will see that it is of fundamental importance to consider the entire distribution of oscillator strength .
we consider a thermal equilibrium system of polaritons , as describing the situation with non - resonant pumping and strong thermalization processes .
this limit has been recently demonstrated to be accessible in on - going experiments .
@xcite the total density of polaritons is fixed by introducing a chemical potential @xmath126 . making use of standard path integral techniques and of the grand canonical ensemble , @xmath127 , we integrate over the fermionic fields , and thus write the partition function @xmath128 $ ] , with @xmath129 , in terms of the following imaginary time action : @xmath130 = \int_0^\beta d\tau \sum_{{{\mathbf p } } } \psi_{{{\mathbf p}}}^ *    \left(\partial_\tau + \tilde{\omega}_{{{\mathbf p}}}\right )    \psi_{{{\mathbf p } } } - \operatorname{tr}\ln g^{-1 } \ ; .
\label{eq : actio}\ ] ] here , @xmath131 is the , energy level diagonal , inverse single - particle green s function : @xmath132 where @xmath133 and @xmath134 are respectively the excitonic and photonic energies measured with respect to the chemical potential . within this path integral formulation , it is possible to show  @xcite that the constraint on the site occupation can be taken into account elegantly by shifting the fermionic matsubara frequencies according to @xmath135 which we will assume from here on in .
we will see later in section  [ sec : photo ] how the formation of polaritons , the result of strong coupling between photons and excitons , influences the optical response and in particular how the optical response is affected when the system enters a condensed phase .
this section will discuss the mean - field equations , and resultant phase diagram ; however first we wish to stress an important point .
as observed previously , the relevant polariton states form from the mixing of a photon with a given momentum to many localised excitonic states with different energies .
the strength of this mixing is determined by the distribution of oscillator strengths .
although the quantum well excitons may be disorder - localised , for a weak disorder potential , with @xmath136 , and in the absence of strong photon disorder ( which in some cases can also be relevant  @xcite ) the resultant polaritons are delocalised . as such
, polaritons will be described by a momentum quantum number , and condense , as in the usual picture , in the lowest momentum state .
the static and uniform minimum @xmath137 of the action   can be found from the equation of motion , or saddle - point equation , which has to be solved together with the mean - field equation for the total number of excitations  : @xmath138\ ; .
\label{eq : homo2}\end{gathered}\ ] ] here , the spectrum of particle - hole excitations , i.e. the quasi - particle spectrum discussed later , is controlled by the energy of an exciton in presence of a coherent photon field @xmath139 , and is given by : @xmath140    in the following , we will solve simultaneously these coupled equations , both above ( @xmath141 ) and below ( @xmath142 ) the critical temperature for condensation .
it is interesting to compare these with similar equations which have been already considered , e.g. , in the bcs theory of superconductivity ,  @xcite or in models describing gases of fermionic atoms close to a feshbach resonance.@xcite considering particularly the case of feshbach resonance , the density of particles is typically kept fixed , while the interaction strength can be varied externally to span across a bec condensate of tightly bound fermionic molecules and a bcs state of loosely bound fermionic pairs . in our case ,
two main important differences arise .
firstly , density is the parameter which is changed .
secondly , the coupling strength of excitons with cavity light depends on the excitonic energy .
there is a wide variety of behaviour that can come from these coupled equations   and  , controlled by the distribution of oscillator strengths and by changing the density of states .    by solving the ` gap equation ' eq .   in the limit @xmath143 ,
let us concentrate for the moment on the critical temperature @xmath144 as a function of the chemical potential @xmath126 , forgetting about the corresponding excitation density .
at the phase boundary , where @xmath145 , in contrast to inside the condensed phase , the squared oscillator strength appears only linearly in eq .  , this expression therefore involving only the optical density @xmath146 .
it is clear that the standard bcs result can be thought of as a limiting case where the effective optical density is much wider than the temperature .
@xcite in the opposite limit , of a narrow delta - like optical density , one can solve the gap equation analytically .
@xcite    by taking into account the density equation  , the critical temperature can be expressed as a function of the excitation density rather than as a function of the chemical potential . in the standard bcs theory
there is no modulation of the interaction strength and therefore no difference between optical density @xmath146 and density of states @xmath66 . in our system instead one has to introduce separately both quantities , and further a realistic description of the dos is required in order to measure the density in physical units . in previous treatments ,
@xcite where the dos has been attributed an arbitrary gaussian - like shape ( and moreover where the coupling strength was kept fixed ) , the density was measured in arbitrary units .
note however that , the fact that the density of states asymptotes to a constant at large densities , and the presence of those exciton states that couple only weakly to light , is responsible only for the asymptotic form of the phase boundary at large densities ( well beyond the range shown in fig .  [
fig : criti ] and beyond the validity of our model ) . at smaller densities ,
such states have a little effect on the phase boundary . in order to calculate the phase boundary using our numerically calculated distributions
, we will solve eq .
( [ eq : homo1 ] ) and eq .
( [ eq : homo2 ] ) numerically .
how the disorder averaged sums in these equations are to be evaluated has been discussed sec .  [
sec : numer ] .     versus the dimensionless density @xmath147 for effective zero detuning @xmath148mev ( @xmath149 is the energy at which the excitonic optical density has its maximum , and so as seen in fig .
[ fig : optic ] , @xmath150 implies @xmath151mev ) ( solid black ) and for positive detuning @xmath152mev ( @xmath153mev ) ( dotted - dashed black ) and @xmath4mev .
the mean - field boundaries for the two different detunings are cut - off by the expected linear dependence of the critical temperature , as indicated at extremely low densities [ orange ( gray ) solid and dotted - dashed lines ] .
the horizontal dashed line marks the temperature of @xmath125k which will be used in later figures . ]
the critical temperature @xmath144 at which the system , for a given density , goes from a non - condensed phase to a condensed phase , is plotted in fig .
[ fig : criti ] for two different values of the detuning @xmath154 , where @xmath149 is the energy at which the excitonic optical density has its maximum ( see fig .  [
fig : optic ] ) .
it is clear that increasing the detuning of the cavity mode above the excitonic one decreases the overall photonic fraction and therefore decreases the value of the critical temperature at a given density of excitations .
conversely , if one fixes temperature , the critical density for condensation is higher for larger detuning , as fig .
[ fig : criti ] clearly shows .
this dependence is easy to understand in the extremely low density limit , where the critical temperature varies linearly with density , and is inversely proportional to the polariton mass : by increasing the detuning , the polariton mass also increases .
having discussed the thermodynamics via the mean - field equations , we now turn in the rest of this paper to the analysis of the optical responses and probes of microcavities , both below and above threshold for condensation , with particular emphasis on the resonant rayleigh scattering emission .
such optical responses are described by considering fluctuations about the mean - field theory outlined in this section .
the spectrum of these fluctuations , and the various optical response functions that these fluctuations describe is the subject of sec .
[ sec : photo ] .    as well as describing the optical response , fluctuations about the mean field theory can also be important in describing , in the extremely low density limit , corrections to the mean field thermodynamics .
since the subject of this paper is optical probes , we will not discuss the subject of fluctuation corrections to the critical temperature here , as they have been discussed elsewhere  @xcite .    for low densities
, a bosonic form can be recovered from our hamiltonian by using the holstein - primakoff transformation to express the two - level systems in terms of bosonic operators , and so this low density limit can of course be treated within our model .
a smooth crossover between the low density ( fluctuation dominated ) , and the higher density ( long - range interaction dominated ) regimes can be found by considering fluctuation corrections to the density , as was studied in ref  .
the optical absorption and emission spectra of a microcavity can be derived by first finding the green s function describing photon propagation .
the non - condensate part of these response functions ( i.e. away from zero momentum ) may be found in practice by considering the green s function for fluctuations about the mean - field solution , @xmath155 , where @xmath156 is a bosonic matsubara frequencies , and expanding the action   up to quadratic terms : @xmath157 this gives the matsubara inverse photonic green s function : @xmath158 physical response functions can be found by analytic continuation of the imaginary time matsubara green s function to real times .
@xcite the matrix elements of @xmath159 can be expressed in terms of the bare photon energy @xmath160 and the excitonic quasi - particle energy @xmath161 as follows ,    @xmath162    where @xmath163    it is useful to decompose the photonic green s function into its momentum diagonal and off - diagonal contributions : @xmath164 by treating the off - diagonal terms perturbatively , the translational invariance can be recovered and polariton eigenstates can be labelled by momentum vectors . in the next section , we will focus on the diagonal terms , which characterise both the spectral weight and the photoluminescence emission .
section  [ sec : rrsca ] is then dedicated to rrs , for which we will see that the off - diagonal terms are necessary in order to describe the scattering of an incident photon ( via the excitonic component of the microcavity polariton mode ) into directions other than its original direction .
@xcite      secondary emission from a semiconductor microcavity after optical excitation is the source of both incoherent pl and coherent rrs . at short times , this emission is dominated by rrs , the coherent scattering from disorder , and so is at the energy of , and coherent with , the incident radiation . at longer times ,
phonon and particle - particle scattering destroy coherence , and redistribute the energy , leading to a quasi - equilibrium distribution of energies , and thus the incoherent pl emission intensity , @xmath165 is given by the bose occupation factor @xmath166 times the spectral weight : @xmath167 the spectral weight can be interpreted as an absorption coefficient  @xcite ( the probability to absorb a photon minus the probability to emit a photon ) , where negative values of @xmath168 represent gain .
in contrast , the pl @xmath169 is always positive .    in calculating the pl , it is convenient to make an approximation by neglecting multiple polariton scattering , while still including the effects of exciton - disorder scattering to all orders .
this is discussed in ref .  , where comparison between this approximation and exact numerical calculations show this approximation to be remarkably good .
physically , this is a good limit to consider because the typical exciton - disorder scattering times are very short compared to the inverse frequencies considered in pl .
as pl depends on the momentum diagonal part of the photon green s function , neglect of multiple scattering means in practice averaging over disorder realisations at the level of the inverse photon green s function , @xmath159 .
since off - diagonal terms in the inverse green s function , eq .  , break translational invariance , they average to zero , and so scattering between different photon momentum states can thus be neglected .
the off - diagonal terms neglected here will however play a crucial role in the case of rrs response , as discussed below .
the spectral weight calculated from this formula , along with the quasi - particle density of states is shown in fig .
[ fig : spwei ] .     as a function of the dimensionless momentum @xmath2 and
rescaled energy @xmath3 , for zero detuning @xmath170 = 0 ( @xmath171mev ) , @xmath4mev , @xmath5k and photon mass @xmath172 .
right column : plot of the quasi - particle dos for the same choice of parameters as the left column , respectively : ( a ) non - condensed ( @xmath7 , @xmath173mev ) ( the bare exciton and photon dispersions [ blue ( dark gray ) dotted line ] and the upper and lower polariton dispersions [ red ( gray ) solid line ] obtained from the effective coupled oscillator model are shown for comparison ) ; ( b ) condensed ( @xmath174 , @xmath175mev ) ; ( c ) condensed ( @xmath9 , @xmath176mev ) . ]
we first discuss this in the non - condensed case , where @xmath177 , and so @xmath178 , and the photon s green function becomes diagonal in particle - hole space , and simplifies to : @xmath179 this expression describes the coupling of one harmonic oscillator ( the photon mode ) to many harmonic oscillators ( the exciton modes ) . in the limit of small density ,
the chemical potential is far below all exciton modes , and so @xmath180 . in this limit ,
the result is identical to a bosonic description of excitons sometimes also called the _ linear dispersion model_. the underlying fermionic structure appears as a reduction of the effective exciton - photon coupling , due to saturation effects described by the @xmath181 term , and discussed in more detail below .    from this _ linear dispersion model
_ , describing one oscillator coupled to many , there are in general two broadened modes at high and low energies  here these are the lower polariton ( lp ) and upper polariton ( up )  and a continuum of modes associated with the exciton optical density . however , when considering the corresponding spectral weight , modes are weighted by their photonic component , so one primarily sees the lp and up modes , and only weak emission near the bare excitonic states  the excitonic ` dark ' states .
these three features are clearly visible in fig .
[ fig : spwei ] , and have been previously predicted by houdr _
et al . _
@xcite making use of a simplified and exactly solvable model
. note that in the limit of an infinitely narrow spectral width , these dark exciton states become entirely dark , having vanishing photon component .
the presence of ` dark ' exciton states , coexisting with strong coupling polaritons can explain the simultaneous observation of large rabi splittings and long decay times seen in some experiments .
@xcite    when the rabi splitting is substantially larger than the exciton inhomogeneous broadening , ( in our case we have @xmath182mev and the fwhm of the optical density @xmath183mev ) , there is a substantial difference between broadening of lower and upper polaritons .
this is because the high energy tail of the optical density decays as a power law , while the low energy lifshitz tail decays exponentially .
thus , the optical density has a larger value at the up mode than at the lp mode , giving a larger broadening compared to the almost vanishing width of the lp .
this description of the polariton linewidth due to the excitonic inhomogeneous broadening coincides with that of whittaker  @xcite , and it has been well tested experimentally .
@xcite    in addition , the location of the lp and up can be found by making use of an effective two - oscillator model , i.e. assuming a narrow delta - like optical density , @xmath184 where @xmath149 is the location of the maximum optical density and effectively , the exciton energy . in this case
the system reduces to two coupled oscillators , giving unbroadened lp and up poles at @xmath185 where , @xmath186 is the reduced rabi splitting , due to saturation effects at higher densities .
the reduction of @xmath187 splitting thus translates directly into a reduction of the lp  up splitting , and thus a blue - shift of the lp , which at small densities can be shown to be linear . for comparison ,
the results of this formula are shown by the solid ( red ) lines in panel ( a ) of fig .
[ fig : spwei ] . in sec .
[ sec : blues ] we will discuss the calculation of this linear blue - shift in the low density regime , by evaluating @xmath187 as a function of density , making use of the full density of states .    let us now turn to the signatures of condensation as seen in the spectral weight , and thus in the pl emission ( but which are most probably masked in the pl emission by strong emission from the condensate mode ) .
when condensed , the polariton modes are replaced by new collective modes :  @xcite the lower polariton becomes the linear goldstone mode , and two branches appear below the chemical potential .
the appearance of new excitation branches below the chemical potential is generic to condensation ,  @xcite however the experiments required to probe these modes are not easy in other bose condensed systems , such as atomic gases . for this reason , let us briefly discuss the physical origin of these new branches , and the reason they may be observed in optical response of polariton systems .
the bogoliubov spectrum arises because of the possibility of processes that spontaneously either create or destroy two non - condensed particles .
( such processes arise due to scattering from or to the condensate . ) as a result of these processes , there is mixing between the propagation of an extra particle , or propagation of a missing particle ( i.e. a hole ) .
( such a language of particle and hole refers to the normal state quasi - particles , in the current case the normal state polaritons . ) in the normal state , one can separately calculate the spectral weight of particle excitations , which have weight only above the chemical potential , and hole excitations which have weight only below the chemical potential . when condensed , as the `` particle '' and `` hole '' spectral weights become mixed , this mixing can lead to spectral weight below the chemical potential that is associated with particle propagation . to observe this weight
it is however necessary to have the ability to inject a particle which is not a quasi - particle of the condensed system , i.e. not a bogoliubov quasiparticle . in atomic experiments
this is hard to achieve , but for polaritons can be naturally achieved by injecting a photon .
these new branches below the chemical potential are seen as optical gain in the spectral weight ( see fig .  [
fig : spwei ] ) . note that the presence of pumping and decay will modify the linear dispersion of the goldstone mode at low momentum , making it diffusive ,  @xcite and that quantisation by disorder may also have some effect .
@xcite    the spectral weight also contains information about the excitonic quasi - particle dos of the system .
as discussed above for the non - condensed case , this is visible via the appearance of ` dark ' exciton states [ see figs .  [
fig : vsdet ] and  [ fig : spwei](a ) ] .
when condensed , there is a coherent field that modifies the energies of these excitonic quasi - particle states , as given in eq .
( [ eq : quasi ] ) , and thus modifies their density of states .
the density of states given by taking the energy of such modes is shown by the right column of fig .
[ fig : spwei ] , and can be compared to the corresponding faint features seen in the photon spectral weight in the left column .     for @xmath90 as a function of the detuning @xmath170 and energy @xmath188 , for a fixed and very low value of density ( @xmath189 , @xmath190mev ) and @xmath191k ( @xmath192 ) .
the bare exciton and photon dispersions [ blue ( dark gray ) dotted line ] and the upper and lower polariton dispersions [ red ( gray ) solid line ] obtained from the effective coupled oscillator model are explicitly shown . ]
the change to the quasi - particle density of states that occurs on condensation requires some explanation . as mentioned earlier , the change of this spectrum is similar to that seen in the bcs theory of superconductivity , as such it is surprising that there is no gap in the density of states of fig .  [
fig : spwei ] .
in fact , there would be a gap if the distribution of oscillator strengths were replaced with the mean square oscillator strength . in that case
, there would be a branch cut , and a @xmath193 singularity in the density of states .
however , with the full distribution of oscillator strengths , one finds that there is always a non - vanishing probability of arbitrarily small oscillator strengths .
( note that if the oscillator strength were zero , the corresponding state would not contribute to any photon response , but for an arbitrarily small coupling , it has some contribution . ) thus , the contribution of these weakly coupled exciton states in effect smoothes out the gap .
thus , this system can be vaguely described as ` gapless fermion condensation ' , analogous to ` gapless superconductivity ' , but through a mechanism very different to the standard abrikosov - gorkov mechanism considered in superconductors .
@xcite    as one moves away from the chemical potential , the mixing of ` particle ' and ` hole ' modes described by the bogoliubov spectrum reduces . as a result , far above the chemical potential , the modes are much as the uncondensed case , and the far below the chemical potential , the new modes disappear . to explain this quantitatively , it is clearer to discuss the case of weakly interacting bose gas .  @xcite writing the bogoliubov mode energy as @xmath194 , where @xmath195 is the mean - field value , the spectral weight is given by : @xmath196 it is thus clear that at large momenta , where @xmath197 , the coefficient of the first delta function ( modes above the chemical potential ) will go to one , and the the coefficient of the second ( modes below ) will be suppressed to zero ( roughly quadratically in energy ) .
in contrast to the power law suppression of the spectral weight of modes far below the chemical potential , the pl signal from modes far above zero is suppressed exponentially by the thermal occupation of these modes [ see eq .  ] , while below the chemical potential , there is no such decay .
however , this discussion neglects emission from the condensate mode , which should included at zero momentum , as defined in eq .  .
to see a noticeable change of the normal modes requires a relatively large condensate density .
due to instrumental broadening , the presence of emission from this large condensate density might obscure emission from the normal modes , so they may only be weakly visible , as shown in fig .
[ fig : plumi ] . as discussed below
, rrs may provide a means of escape from this problem .     as a function of the dimensionless momentum @xmath2 and
rescaled energy @xmath198 , for detuning @xmath199mev , @xmath4mev , @xmath5k and @xmath200 : ( a ) non - condensed ( @xmath7 , @xmath173mev ) ;
( b ) condensed ( @xmath201 , @xmath202mev ) ; ( c ) condensed ( @xmath9 , @xmath203mev ) . ]      before discussing the rrs response , let us briefly discuss a point mentioned above  the calculation of the lp blue - shift as a function of density , making use of the full density of states , and our model of saturation effects  .
this can be observed in pl experiments by the change of blue - shift as a function of intensity of the non - resonant pumping .    in the non condensed regime ,
the blue - shift of the lp in our model is a consequence of the saturation of the ( disordered ) energy levels . following the narrow band - width limit described in eq .
( [ eq : lpeup ] ) , this saturation blue - shift can be found if one has an expression for the chemical potential at a given temperature and density . in the low density limit ,
the expression for density @xmath147 as a function of @xmath126 [ eq .
( [ eq : homo2 ] ) with @xmath177 ] can be inverted in terms of elementary functions . as illustrated by fig .
[ fig : occup ] , in the low density limit , the chemical potential is far below the band edge .
thus , there are two significant contributions to density : one from low energies , in the tail of the dos , but at large occupation , and one from high energies , in the tail of the occupation , but at large dos , thus : @xmath205 here @xmath206mev is the energy which characterises the exponential decay of dos in the tail of the lifshitz states  this coefficient is extracted by an exponential fit to our numerical dos . at a temperature of @xmath207mev , the dominant contribution to the density is given by the second integral and the chemical potential increases logarithmically with the density : @xmath208 from this expression , and from the coupled oscillator expressions of the lp eq .  , we can explicitly derive the reduced rabi splitting due to saturation effects , @xmath209 and the lp blue - shift , which is thus linear in this low - density regime : @xmath210 where we have reintroduced here the density @xmath211 per unit volume .
remaining at small densities , but now considering the low temperature limit , @xmath212 , then the dominant term in eq .
( [ eq : density - assympt ] ) is the first , and so eq .  ( [ eq : appbl ] ) should be modified by replacing the thermal length @xmath213 , by a characteristic disorder length @xmath214 .
it is instructive to compare this result with that for a clean system ; in this case blue - shift of the lp has been attributed either due to coulomb interaction  @xcite or to saturation effects  @xcite ( where the expression given here is valid only in the dilute limit ) : @xmath215 because the excitonic rydberg @xmath216 can be of the same order of magnitude as the rabi splitting @xmath104 , in a clean system the two shifts can be expected to be of the same order of magnitude . considering just saturation effects , the difference between clean and dirty systems is that , at low temperature , in a clean system blue - shift depends on exciton number per square bohr radius , while in a disordered system the relevant length is that characteristic of the disorder potential , which is in general larger than the bohr radius .
finally , we wish to observe that there is a distinction between coulomb and saturation effects : coulomb interactions result in a blue - shift of both the lp and up , while saturation leads to a blue - shift of the lp and a red shift of the up , i.e. a collapse of the rabi splitting .
thus , these effects can be experimentally distinguished , and their relative magnitudes determined .    in the condensed regime ,
the equivalent of the lp is the goldstone linear mode , which by definition starts at the chemical potential .
thus , the observed blue - shift is a direct observation of chemical potential vs.  density .
figure  [ fig : vsden ] shows the variation of the spectral weight at zero momentum as a function of density , from which the energy of the zero momentum lp and up modes can be extracted . also shown , for comparison , is the chemical potential vs.density ( dashed green line ) . the locking of the lp mode to the chemical potential beyond the critical density ( red solid line ) is clearly visible , and the behaviour of the lp mode vs density allows the extraction of density from an experimental measurement of lp blue - shift .     for @xmath90 as a function of the dimensionless density @xmath147 and the rescaled energy @xmath217 , for zero detuning ( @xmath218mev ) and @xmath191k ( @xmath192 ) .
the rescaled chemical potential @xmath219 [ green ( light gray ) dashed ] , the non - condensed lower and upper polariton modes from the two coupled oscillator model [ blue ( dark gray ) dotted ] , and the critical density for condensation @xmath220 [ vertical solid red ( dark gray ) line ] are shown for comparison . ]      as described in refs .
, the rrs intensity describes the probability to inject a photon into the cavity with momentum @xmath30 , and detect a coherent photon at @xmath221 with @xmath222 .
the intensity of such a signal may be written as : @xmath223 since this involves the off - diagonal momentum space components of the green s function , it is necessary to include the off - diagonal momentum components of the inverse green s function eq .  .
however , as discussed in sec .
[ sec : photo ] , we will consider only a single polariton scattering and again neglect multiple polariton scattering ( but include all orders of exciton - disorder scattering ) .
such an approximation is reasonable for the same reasons discussed at the start of sec .
[ sec : photo ] .    in the following , we propose using the rrs signal to probe the excitation spectrum in the presence of a polariton condensate .
it is therefore necessary to be able to separate the rrs signal from the strong photoluminescence that would arise from the equilibrium state with large polariton density .
further , since a condensed polariton system may have a strong nonlinear response to an applied probe ( i.e. stimulated scattering if the probe significantly affects the population of polariton modes ) , it is necessary to use a weak rrs probe .
fortunately , the coherent nature of rrs allows exactly this : one can detect a weak rrs signal by phase sensitive measurement , as is discussed in detail in appendix  [ sec : phase - sens - detect ] .
hence , the limit on intensity of the probe is provided by the sensitivity of the ccd camera , and not by the incoherent photoluminescence background .
this therefore allows a probe sufficiently weak that only the linear rrs response is seen , and nonlinear effects can be avoided .
when non - condensed , the single polariton scattering expansion of the photon green s function in the off - diagonal terms gives @xmath224 the factors @xmath225 appearing here can be interpreted as a filter , allowing a response only when the outgoing ( or equivalently incoming for @xmath226 ) momentum has an energy close to the polariton mode at the given energy .
this means that @xmath227 and so is responsible for the ring - shaped rrs signal observed in experiments .
@xcite in contrast , the term @xmath228 describes scattering between momentum states via polariton - exciton - polariton scattering .
this term has a large variation from disorder realization to disorder realization , and is the reason for the speckle seen in rrs experiments .
this speckle , and the disorder averaged rrs intensity are shown for comparison in fig .
[ fig : excsp ] .
as the precise speckle pattern depends on the precise disorder realization , the most we can reasonably do is to describe the statistical properties of this speckle , thus the disorder averaged rrs signal shown in fig .
[ fig : rrsca ] would in experiment represent the envelope of the rrs speckle pattern .     for a single disorder realization ( gray dotted line ) and its disorder average ( black cross symbols ) versus energy for @xmath229 and a @xmath230 azimuthal angle , @xmath231 and @xmath232 . for comparison , we also plot the excitonic optical density [ green ( light gray ) plus symbols ] . ]
when condensed , the expression for the rrs intensity becomes more complicated , but can still be written in the form of a filter and scattering part : @xmath233 where the filter and scattering terms are respectively : @xmath234 as in the non - condensed case , the filter function @xmath235 restricts the allowed incoming and outgoing momenta to those for which the normal modes have the injected energy .
however , as discussed earlier , in the presence of a condensate , the normal modes supported are no longer the lp and up , but instead the new bogoliubov - like quasi - modes .    the decay of spectral weight associated with modes far below the chemical potential , as discussed in sec .
[ sec : spect ] still applies ( this can clearly be seen in fig .
[ fig : rrsca ] , fig .
[ fig : spwei ] and fig .
[ fig : plumi ] ) .
however , since the energy of the rrs signal is controlled by the energy of the incident photon , it is possible to study modes far above the chemical potential , which would have negligible equilibrium occupation , and thus negligible weight in the pl signal .
both for this reason , and because of the ability to discriminate between the coherent rrs signal and inchoherent photoluminescence , rrs provides a very powerful tool to study the interesting properties of condensed polariton systems .
in this paper we have analysed and compared different optical responses of a condensate of microcavity polaritons .
one of the conclusions of this work is that , compared to photoluminescence studies , the energy and momentum resolved resonant rayleigh scattering spectra provides many advantages in probing a condensed phase and in studying the associated coherent excitation spectrum . in particular ,
resonant rayleigh scattering , by collecting the signal at a different angle from that at which the system is probed , allows one to look directly at characteristic features of the condensed phase excitation spectrum .
such features , like the goldstone modes both above and below the chemical potential , are weak in photoluminescence emission and are likely to be masked by the strong emission from the condensate in the lowest momentum state . in contrast , rayleigh scattering ,
being coherent with the probe , can be distinguished from condensate emission .
one can ask why the accurate description of excitonic disorder considered in this paper is necessary , when the microcavity photon introduces a long length scale which averages over disorder . one might also compare this system with other excitonic systems , such as double quantum wells ,  @xcite where high density of dipole - dipole interacting excitons quickly screens the disorder .
@xcite the answer to this question is that the large ratio of exciton to photon mass means firstly that exciton density required for condensation in microcavities is far less than in double quantum wells ; and secondly that those exciton states involved in forming the thermally occupied polariton modes are strongly localised exciton states ( i.e. influenced by disorder beyond born approximation ) .
one of the consequences of an accurate treatment of disorder in microcavities is that the blue - shift of the lower polaritons due to saturation effects depends on a length scale characteristic of disorder and temperature , which can be much larger than the bohr radius , this length playing an equivalent role in the clean system .
this effect can be important in determining the polariton densities from the measured blue - shift in current experiments .
an important consequence of disorder is that the exciton - light coupling strength is characterised by a full distribution . in the non - condensed phase
, the optical density alone determines the photoluminescence response , while for resonant rayleigh scattering the averaged fourth power of the oscillator strength is required .
however , in the condensed phase one has to consider the entire distribution for each exciton energy . for energies close to and above the band edge , there is always a non - vanishing probability of an arbitrarily small oscillator strength .
for this reason , in contrast to examples like bcs superconductivity , one can show that the quasi - particle spectrum does not have a hard gap .    in conclusion , we have considered how an accurate treatment of disorder on the single - particle excitonic level , when elevated to the many body problem of interacting microcavity polaritons , leads to a variety of interesting features in various optical responses and probes of the condensed phase .
we are grateful to roland zimmermann , wolfgang langbein and ben simons for suggestions and useful discussions .
f.m.m . and m.h.s
would like to acknowledge financial support from epsrc .
j.k . would like to acknowledge financial support from the lindemann trust .
this work is supported by the eu network `` photon mediated phenomena in semiconductor nanostructures '' hprn - ct-2002 - 00298 .
in this appendix we show how the coherent nature of the resonant rayleigh signal allows a weak coherent probe to be detected in the presence of strong incoherent photoluminescence .
the rayleigh scattering probe is a perturbation , @xmath236 , which in the rotating wave approximation may be written as : @xmath237 this describes a probe at wavevector @xmath30 , frequency @xmath238 , of strength @xmath239 .
the rrs signal is the coherent scattering of this probe to other wavevectors . to isolate the part of the emission that is coherent with the probe
, one may use a homodyne measurement , interfering the emission with part of the probe signal .
this corresponds to measuring the spectrally resolved emission intensity : @xmath240 where @xmath241 \left[\hat{\psi}^{}_{{\mathbf q}}(0 ) +    a_1e^{i\phi } \right ] |n\rangle \ ; , \end{gathered}\ ] ] and where @xmath242 is the strength of the homodyne mixing , and @xmath243 is a phase delay introduced between the probe and the homodyne signal .
the states @xmath244 are the eigenstates of the system without the probe , given by @xmath245 ; @xmath246 is the free energy , for normalization .
the principle of phase sensitive detection is to vary the phase delay , and to extract the part of @xmath247 which depends on this phase delay .
since the background pl does not depend on the phase delay , this allows one to separate a small but phase dependent signal from a strong but phase independent background .
the part of the signal that depends on phase is : @xmath248 \ ; . \label{eq:4}\end{gathered}\ ] ] standard first order time - dependent perturbation theory , using the perturbation in eq .
yields : @xmath249e^{-i \omega t } + [ \ldots ] e^{i\omega t } \ ;    .
\label{eq:5}\end{gathered}\ ] ] the term written @xmath250 $ ] is similar to the first term , but with @xmath251 and @xmath252 . inserting this signal into eq .   and then into eq .
( [ eq:1a ] ) , it can be seen that such a term gives a signal at frequency @xmath253 , and so can be clearly separated and ignored . inserting the first term into eq .
, one can write : @xmath254 \frac { \langle n |    \hat{\psi}^{\dagger}_{{\mathbf q } } | m \rangle \langle m |    \hat{\psi}^{}_{{\mathbf p } } | n \rangle } { e_n - e_m - \omega - i \eta } \ ;    , \label{eq:6 } \end{gathered}\ ] ] which is the definition of the retarded green s function ,  @xcite @xmath255 .
by repeating the same analysis for the second term in eq .
( [ eq:4 ] ) , one can then write : @xmath256 where @xmath257 is the phase of the retarded green s function .
hence , the phase - dependent part of the luminescence is given by the off - diagonal in momentum space part of the retarded green s function , i.e. the rrs signal .
h. deng , d. press , s. gotzinger , g. s. solomon , r. hey , k. h. ploog , y. yamamoto , phys .
* 97 * , 146402 ( 2006 ) .
j.  kasprzak , m.  richard , s.  kundermann , a.  baas , p.  jeanbrun , j.  keeling , f.  m.  marchetti , m.  h.  szymanska , r.  andr , j.  l.  staehli , v.  savona , p.  b.  littlewood , b.  devaud , and le si dang , nature , * 443 * , 409 ( 2006 ) .
m. h. szymanska and p. b. littlewood , solid sate commun . * 124 * , 103 ( 2002 ) ; m. h. szymanska , p. b. littlewood , and b. d. simons , phys .
a * 68 * , 013818 ( 2003 ) .
a. kavokin and g. malpuech , _ cavity polaritons _ , in thin films and nanostructures vol . 32 ( elsevier , ny , 2003 ) .
j. keeling , p. r. eastham , m. h. szymanska , p. b. littlewood , phys .
lett . * 93 * , 226403 ( 2004 ) ; phys .
b * 72 * , 115320 ( 2005 )
. f. m. marchetti , b. d. simons , and p. b. littlewood , phys .
b * 70 * , 155327 ( 2004 ) ; f. m. marchetti , m. h. szymanska , p. r. eastham , b. d. simons , and p. b. littlewood , solid state commun . * 134 * , 111 ( 2005 ) .              d. sarchi and v.savona , aip conf . proc . *
893 * , 401 ( 2007 ) .
p. schwendimann and a. quattropani , phys .
b * 74 * , 045324 ( 2006 ) .
r. zimmermann , f. grosse , and e. runge , pure and appl .
chem . * 69 * , 1179 ( 1997 ) ; r. zimmermann and e. runge , phys .
( a ) * 164 * , 511 ( 1997 ) .
i. carusotto and c. ciuti , phys .
lett . * 93 * , 166401 ( 2004 ) ; c. ciuti and i. carusotto , phys .
status solidi b * 242 * , 2224 ( 2005 ) .
s. d. baranovskii and a. l. fros , fiz .
. poluprovodn .
* 12 * , 2233 ( 1978 ) [ sov .
. semicond .
* 12 * , 1328 ( 1979 ) ] .
r. zimmermann , e. runge , and v. savona , `` theory of resonant secondary emission : rayleigh scattering versus luminescence '' , in _ quantum coherence , correlation and decoherence in semiconductor nanostructures _ ( ed .
t. takagahar ) , elsevier science ( 2003 ) .
m. holland , s.j.j.m.f .
kokkelmans , m.l .
chiofalo , and r. walser , phys .
* 87 * , 120406 ( 2001 ) ; e. timmermans , k. furuya , p.w .
milonni , and a.k .
kerman , phys .
a * 285 * , 228 ( 2001 ) ; y. ohashi , a. griffin , phys .
lett . * 89 * , 130402 ( 2002 ) .