Thermo-optical interactions in a dye-microcavity photon Bose–Einstein condensate

Superfluidity and Bose–Einstein condensation are usually considered as two closely related phenomena. Indeed, in most macroscopic quantum systems, like liquid helium, ultracold atomic Bose gases, and exciton-polaritons, condensation and superfluidity occur in parallel. In photon Bose–Einstein condensates realized in the dye microcavity system, thermalization does not occur by direct interaction of the condensate particles as in the above described systems, i.e. photon–photon interactions, but by absorption and re-emission processes on the dye molecules, which act as a heat reservoir. Currently, there is no experimental evidence for superfluidity in the dye microcavity system, though effective photon interactions have been observed from thermo-optic effects in the dye medium. In this work, we theoretically investigate the implications of effective thermo-optic photon interactions, a temporally delayed and spatially non-local effect, on the photon condensate, and derive the resulting Bogoliubov excitation spectrum. The calculations suggest a linear photon dispersion at low momenta, fulfilling the Landau’s criterion of superfluidity. We envision that the temporally delayed and long-range nature of the thermo-optic photon interaction offer perspectives for novel quantum fluid phenomena.


I. INTRODUCTION
For a gas sufficiently cold and dense that the thermal de Broglie wavelength exceeds the interparticle spacing, quantum statistical effects come into play.Specifically, for massive bosonic particles, Bose-Einstein condensation into a macroscopically populated ground state minimizes the free energy, as has been experimentally demonstrated in the gaseous regime with ultracold atoms more than 20 years ago 1 .
More recently, Bose-Einstein condensation has also been observed with exciton-polaritons, mixed states of matter and light, and with photons [2][3][4][5] .Unlike particles with a non-vanishing rest mass, photons usually do not show Bose-Einstein condensation.The thermal photons of blackbody radiation have no chemical potential, corresponding to a non-conserved particle number upon temperature variation.Therefore, blackbody radiation photons vanish at low temperature instead of showing a phase transition to a condensate.This difficulty was overcome in 2010 by confining photons in a dye-solution filled optical resonator made of two mirrors spaced in the micrometer regime 3,5 .The short mirror spacing effectively imprints a low-frequency cutoff for the photon gas, and in the presence of a mirror curvature the problem becomes formally equivalent to a two-dimensional system of harmonically confined massive bosons.By repeated absorption and re-emission processes, photons thermalize to the (rovibrational) temperature of the dye solution at room temperature.
Experimentally, both the thermalization of the photon gas to the dye temperature 6 and Bose-Einstein condensation of photons above a critical particle number have been observed 3,5 .For larger condensate fractions, the size of the condensate increases, which is attributed to a weak repulsive interaction mainly due to the thermo-optic effects.
So far, it is not known if the photon gas, in addition to exhibiting Bose-Einstein condensation, is a superfluid 7 .The existence of superfluidity is believed to require direct interparticle interactions, as e.g.present for polaritons, for which superfluidity has been established [8][9][10][11][12][13] .
On the theoretical side, the concept of a nonlinear photon fluid was first introduced by Brambilla et al. and Staliunas 14,15 , who used hydrodynamic equations to describe electromagnetic fields in a cavity.Chiao et al. subsequently proposed to generate a photon superfluid in a nonlinear optical cavity using the Kerr-effect and furthermore predicted sound-like modes at the low-momentum part of the Bogoliubov dispersion 16,17 .
In this paper, we examine the effect of thermo-optic interactions on the dispersion of a photon gas trapped inside an optical microcavity.The thermo-optic effect, also known as thermal lensing, was first introduced by J. Gordon while studying transient effects of the output power and the beam size upon inserting a liquid cell inside the laser cavity 18 .In a propagating configuration in stationary conditions, the temporal delay of the thermo-optic effect does not play a role.The signatures of superfluidity and non-local effects from the associated thermo-optic nonlinearity have been experimentally observed in this configuration 19,20 .In another work, Strinati and Conti have theoretically studied the stationary state of dye microcavity photon condensate subject to a non-local thermo-optic nonlinearity 21 .
We report a theoretical study on the effects of a thermo-optic nonlinearity giving rise to a temporally delayed and non-local effective photon interaction for the photon condensate in the dye-microcavity system.Assuming a plane mirror microcavity geometry, we derive the Bogoliugov dispersion for such a system subject to small perturbations.For a suitable parameter range, we find a linear dispersion at low momentum, corresponding to a phonontype dispersion.We discuss the possibility of superfluidity of the photon condensate based on such temporally delayed and long-range interactions.
In the following, Chapter II discusses some general properties of the photon gas subject to the thermo-optic interactions, and section III gives steady-state solution.In section IV we derive Bogoliubov modes of the system and the resulting elementary excitation dispersion.
Finally, section V concludes the manuscript.

II. THERMO-OPTIC INTERACTIONS IN A PHOTON GAS
We begin by discussing some general formulas describing the system of a photon gas trapped in a dye microcavity subject to thermo-optic interactions.In the experimental scheme of 3,4 , Bose-Einstein condensation of photons is achieved in a dye microcavity (see Fig. 1) with the mirror spacing in the wavelength regime.This leads to a large frequency spacing between longitudinal modes, comparable to the emission width of the dye molecules.
Therefore, to good approximation only photons of a fixed longitudinal mode order q are found in the resonator, and the two remaining transverse mode quantum numbers make the system two-dimensional.The transverse T EM 00 -mode has the lowest allowed frequency, which imposes a low-frequency cutoff.Moreover, the photon dispersion becomes quadratic, i.e. massive particle-like, and the mirror curvature imposes a harmonic confinement on the photon gas.
One can show that the photon gas in the cavity is formally equivalent to a two-dimensional gas of massive bosons with effective mass m ph = k z n 0 /c = ω c (n 0 /c) 2 , where ω c denotes the cutoff frequency, k z the longitudinal wavevector, and n 0 the refractive index of the solution.
The photon energy in the paraxial limit is: where k r denotes the transverse momentum and is the trapping potential.The trapping frequency of the harmonic potential is given by Ω = (c/n 0 )/ LR/2, where L and R are the mirror spacing and curvature, respectively.Finally E int is the effective photon interaction energy, as will be discussed below.
Thermal equilibrium of the photon gas is achieved by repeated absorption and re-emission processes by the dye molecules.For the described two-dimensional, harmonically confined system it is known that a Bose-Einstein condensate exists at finite temperature.Accounting for the two-fold polarization degeneracy of photons, one finds the critical particle number: For Ω = 2π × 3.6 × 10 10 Hz, as derived for L = 2 µm, R = 1 m, and n 0 = 1.33, at room temperature (T = 300 o K) one obtains N c ≈ 99000 for the critical particle number.
During the course of the absorption re-emission processes of the dye molecules, a small fraction of inelastic processes due to dye's finite quantum efficiency (η ≃ 95% for the case of rhodamine dye solution) causes local heating of the solvent.Due to the temperature dependence of the solution refractive index, the optical distance between the mirrors is decreased at the corresponding transverse position in the cavity.This is equivalent to a local rise of the photon gas potential.In other words, the heating with a corresponding decrease of the refractive index results in a smaller optical wavelength, hence a higher photon energy is required to locally match the mirror boundary conditions.The resulting interaction energy is: where ∆n = β∆T , with β = ∂n/∂T as the thermo-optic coefficient of the solution.
The photon gas is well described by an equilibrium state, if thermalization by coupling to the dye molecules is faster than both loss and pump processes.In this case, photons can relax towards the ground state of the harmonic trapping potential and form a Bose-Einstein condensate before they are lost through mirror transmission or inelastic processes in the dye 22,23 .Accounting for the thermo-optic effective photon interactions, the condensate dynamics can be described with a generalized time-dependent Gross-Pitaevskii equation: In the above equation ψ( r, t) denotes the slowly varying envelope of the condensate wavefunction (in the mean-field treatment) in the rotated frame of ω c .Note that the above equation is more general and captures the full 3D behavior of the wavefunction as (ψ(x, y, z, t)).
Also instead of using the effective trapping potential of eq. 1, we impose Dirichlet boundary conditions on the spherical mirror surfaces.Similar as in cold atom physics literature 24,25 , here we do not consider interaction effects of the thermal cloud, due to its much lower density.
The time evolution of the relative local temperature is determined by the heat transport equation as: The first term on the right hand side accounts for the heat diffusion and the second term for heating through inelastic processes in the dye solution.The parameters κ and C v are the thermal conductivity and volume heat capacity, respectively.Further, α in is the inelastic absorption coefficient and is related to the absorption coefficient α via the quantum efficiency of the dye η, as α in = (1 − η)α.
Moreover, the following normalization condition for the wavefunction is held: where N BEC is the total number of photons in the condensate.Table I gives relevant parameters of a typical dye-filled cavity setup with methanol solvent 23 , 26 .
List of physical parameters of a dye-filled cavity setup used in this paper.The properties of the solvent, methanol, are from. 26The value of the inelastic absorption coefficient is for 1mM solution of R6G in methanol solvent and is calculated from the experimental data in 23 , assuming a quantum yield of 95% for the dye 27 .

III. STEADY-STATE
In the steady state, the temperature is settled to a stationary value T ss ( r), and the the time evolution of the condensate wavefunction is given by ψ( r; t) = ψ ss ( r)e −iµt .The stationary forms of eq. 4 and 5 are: The mirrors (M 1,2 in Fig. 1) are macroscopically large.Moreover, their thermal conductivity exceeds that of the dye solution, so in all of the calculations we assume ∆T = 0 on the mirror surfaces.

A. Numerical solution
In the presence of a mirror curvature, as required to obtain a trapping potential, the coupled sets of eq.7 can only be solved numerically.Without a thermo-optic non-linearity (i.e. for β = 0), the ground state of the harmonic trap is the lowest energy state.Following the usual convention we normalize the chemical potential to yield µ = 0 in this interaction-less case.For earlier work investigating the steady state properties of a thermo-optic interaction in the harmonically trapped condensate inside a microcavity, see Ref. 21 .
Using a fully numerical algorithm we solved the above coupled non-linear equations to investigate the effect of non-linearity on the interacting condensate.The results for a symmetric cavity with mirrors of R = 1 m radius of curvature are shown by the solid lines in A comparison with experimental results is not straightforward since the condensate mode diameter increase reported in 6 corresponds to the accumulative stationary value observed for a pulsed pump with certain repetition rate.The values here are nevertheless smaller than the experimentally observed ones.We attribute this discrepancy to the boundary conditions used for the temperature distribution, which imposes that ∆T vanishes at the mirror surfaces.A more realistic model in the presence of the thermo-optic non-linearity needs to account for the finite thermal-conductivity of the mirrors and include the thermal properties of both mirror layers as well as the substrate.However, we expect the dynamic properties of the condensate, as discussed in the following, to be less affected by the thermal properties of the mirrors since the local properties of the heat transfer will dominate in that case.

B. Green's function approach
Aside from being computationally expensive, the fully numerical method fails to provide one with more physical insight about the problem.In this section we propose an alternative method by employing the Green's function of the heat diffusion problem, making the analysis more efficient.Moreover, this form could be used further for elementary excitation studies as will be discussed in the next sections.
The coupled eigenvalue problem of eq.7 can be efficiently solved using Green's function of heat diffusion problem when proper boundary conditions are applied.The Green's function of the heat transfer problem will be determined as: For simplicity in our analytic calculations the mirrors are assumed to be flat at z = ±L/2, and extended to infinity in the transverse plane (xy-plane).
When combined with the eigenvalue problem of the wavefunction in eq.7 one obtains the following equation for the steady-state: Equation 10 has the form of a Gross-Pitaevskii equation, while the interaction term is different from the usual contact form.Here, the interaction potential has an integral form describing a non-local interaction given by the Green's function, where the strength of the interaction potential at point r is related to the whole condensate wavefunction distribution.
Figure 3 shows the variation of the Green's function G N L ( r; 0) as a function of ρ in the transverse plane for different mirror separations.As can be seen at a fixed radial point the value of the Green's function increases as the cavity length becomes larger.This behavior implies that the temperature changes from a 2D distribution to 3D in thicker cavities.In the limit of very large mirror spacing the Green's function approaches a Coulomb-type longrange interaction given by G N L ( r; r ′ ) ≈ 1/| r − r ′ |.This behavior can be clearly observed in Fig. 3 where the thick cavity Green's function (L = 10 µm) and the asymptotic Columbic form (black dashed line) are in a good agreement.
To compare the predictions of these two approaches, we have used the generalized Gross-Pitaevskii eq. 10 to calculate the stationary features of the condensate.Unlike the numerical method however, the Green's function approach only allows for a treatment of a flat mirror geometry.The red dashed line with squares in Fig. 2(a) shows the chemical potential variation for the homogeneous problem calculated with this method.Similarly, the dashed line with circles in Fig. 2(b) shows the corresponding variation of the maximum temperature with the number of photons in the condensate.The obtained results from the Green's function approach are in approximate agreement with the numerical results for the problem with curved mirrors, and show the same trend.We point out that an exact agreement is not expected here due to the different mirror geometries, with correspondingly different condensate mode profiles.
In contrary to the case of atomic BECs where the range of interaction is strictly limited by the interaction type, the photon fluid shows a unique feature of possessing a tunable interaction range from a local form to a highly non-local gravitational type interaction in relatively thin and thick cavities, respectively.In practice, the requirement of photon Bose-Einstein condensation for a low-frequency cutoff imposes a maximum usable cavity length for corresponding experiments.

IV. SMALL PERTURBATIONS AND BOGOLIUBOV MODES
For weak perturbations the dynamics of the system can be found by assuming small fluctuations around the stationary solutions.We use the following ansatz to determine the Bogoliubov modes, where Ω denotes the frequency of the perturbations.
Inserting this ansatz into the original non-linear equations and neglecting the terms with orders higher than one in the perturbation, one can derive the following linear coupled equations for the small fluctuations u( r), v( r), and δt( r): where: To begin with, we derive a solution that neglects the temporal delay of the temperature and assumes that it follows the time dependency of the condensate density |ψ( r, t)| 2 .With this assumption the temperature change ∆T ( r, t) is given by a diffusion-type equation, as for T ss in eq. 7, which yields a non-local, but instantaneous effective interaction.Later in this section we will modify this assumption to take into account the temporal delay of the temperature distribution due to the finite heat conductivity.
For a translationally invariant problem, a plane wave ansatz of the form e i k• r would be well suited to describe the spatial part of the excitations.More suited for this problem, given that the condensate wavefunction ψ( r) is not spatially uniform, is the Green's function approach, building upon the stationary eigenfunctions discussed in the previous section.The mirrors impose Dirichlet boundary conditions and break the translational invariance along the z-axis.Therefore, we approximate the longitudinal variation of the wavefunction as 2/L sin(qπz/L), and define an effective 2D Green's function in the transverse plane: where K 0 is the 0 th -order modified Bessel function.
As can be inferred, this effective 2D potential is translationally invariant in the transverse plane, and leads to a well-defined transverse momentum k for the elementary excitation.
Therefore, the Bogoliubov dispersion of this system is properly defined for frequency Ω and the transverse momentum k, and from eq. 12 is determined as: Here ĜNL (k) is the Fourier transform of G 2D ef f ( ρ), and in the first formula |ψ ss (0)| 2 can be approximated as: Equation 15 gives the dispersion in the presence of the thermo-optic interaction.At large momenta, the first term on the right hand side dominates, yielding the usual particle-like quadratic dispersion of photons, see also eq. 1.In other words, when the momentum is large the interactions have negligible effect.
The thermo-optic interactions impact the low-momentum part of the dispersion, and the corresponding effect is predominantly determined by the function ĜNL (k), i.e. the Fourier transform of the Green's function.For wavevectors k ≤ k z = qπ/L, within the range that the paraxial limit is fulfilled, ĜNL (k) becomes almost k-independent, hence a linear tendency for low-momentum excitations is expected.
So far we have ignored the explicit dynamics of the delayed temperature given by the Ṫ -term in the heat equation.This effect can be taken into account by employing the proper Green's function of the heat equation in eq. 5 which depends on time as well as the position.
This time-dependent Greens' function modifies the aforementioned dispersion of eq. 15 to a transcendental equation for the dispersion (i.e.Ω(k)) after substituting ĜNL (k) with ĜNLD (Ω, k) as: Equation 18 together with eq. 15 gives the final results for the quasi-particle dispersion, including the effects of both non-locality and the delay in such thermo-optic interactions.As the form of the transcendental equation suggests, the delayed nature of the interaction leads to complex frequencies Ω, implying an instability of the condensate.Compared to the nonlocal case only, this is the main qualitative modification of the temporal effect.We notice this behavior can be compared with the dynamical instability predicted and observed in polariton condensates 28 .While the instability of the latter is due to the coupling of the condensate with the exciton reservoir, in this dye microcavity system the condensate instability occurs due to the thermal coupling to the dye solutions.At longer times, we expect that the thermo-optic interaction destroys the photon condensate for larger interaction strengths.This feature suggests that the photon BEC can potentially be a superfluid, a feature that can be better understood using the non-local effect of the Green's function.
In eqs.15 the effect of non-locality is implicit in ĜNL (k) and ĜNLD (k, Ω).A good physical intuition can be established by studying two extreme cases for this effect.As demonstrated in Fig. 3, the interaction range decreases in thin cavities.In such cavities the interaction ultimately reduces to a local one with 1, leading to a linear dispersion.Therefore, for a contact interaction the dispersion is separated to two distinct forms, one for free particles at large momenta (Ω ∝ k 2 ) and one for sonic modes for small momenta (Ω ∝ k).The other extreme happens for thick cavities where the fluid becomes fully 3D with a long range, gravitational-type interaction where this Green's function would lead to a constant, k-independent dispersion at low momenta.
To provide a better understanding of the dependency of the dispersion in terms of physical parameters such as L, N BEC , and α in , we define a critical momentum k critical , as the largest momentum at which the dispersion is significantly modified by the interaction.Using the above given dispersion relations this parameter is determined as: Where Ĝ is a general representation for the Green's function and could be either of ĜNL (k) for the non-local case, or ĜNLD (Ω, k) for the non-local and delayed one.
Figure 5(a)-(c) shows the behavior of the critical momentum as a function of some experimental parameters.The solid blue lines correspond to the predictions considering both delay and non-locality given by the Green's function of eq.18.For comparison, the dashed red lines are obtained when only considering the non-locality as given by the Green's function of eq.16.While the critical wavevectors have similar dependencies to the experimental parameters, the inclusion of the temporal delay decreases the values.As discussed earlier, the inclusion of the temporal delay also results in complex eigenfrequencies Ω (Fig. 4(b)), leading to condensate instabilities.
The results of Fig. 5 on the inelastic absorption coefficient, the cavity length, and the number of photons in the condensate, respectively, showing that the critical velocity increases with an increase in the inelastic absorption and photon number and decreases with increasing the cavity length.

V. CONCLUSIONS
In this work we investigated the effect of a thermo-optic interaction, a temporally delayed and spatially non-local interaction, on a photon Bose-Einstein condensate in a dye-filled microcavity system.We derived the general form of the dispersion in such systems, calculated the spectrum of Bogoliubov modes in a plane-mirror microcavity, and identified a linear scaling for the low-energy modes.At larger transverse momenta, the usual quadratic freeparticle dispersion of cavity photons is restored.The derived linear dispersion, corresponding to sonic modes, fulfills Landau's criterion for superfluidity.
We envision several experimental and theoretical follow up studies based on the reported results here.To obtain accurate quantitative predictions, it would be important to investigate the implications of a finite heat conductivity of the mirrors on both the stationary and the dynamical features of the condensate subject to the thermo-optic non-linearity.
For an experimental test of superfluidity based on the thermo-optic non-linearity a possible setup could use a photon fluid in a wedge-shaped cavity composed of two tilted flat mirrors, allowing the photon droplet to flow freely.By intentionally making a perturbation on one of the mirrors in the flowing path of the photon fluid, scattering would be observed if the fluid is dissipative.However, in the superfluid phase the photon droplet passes the perturbation without being scattered.By directly monitoring the photon condensate-defect interaction one should be able to distinguish between these two different phases.
Another follow up work could search for long-lived vortices of the photon condensate in a curved mirror microcavity.Along these lines, the implications of such a trapping potential on the Bogoliubov modes and the spectrum of elementary excitations should be studied.
Another fascinating topic for future investigations is the search for analogies between the long-range thermo-optic interaction and gravitational physics and its consequences, including possible black-hole physics .For the Green's function method calculation investigating the flat mirror problem, the quoted value for the photon number refers to the number of photons in an area πr 2 BEC .r BEC corresponds to the interaction less condensate radius of the curved mirror problem.In this way, the assumed average photon density is the same for both of the numerical and the Green's function approaches.thermo-optic effective photon interaction, given by eq. 15 and eq.18 for different photon numbers N BEC .Other parameters are L = 2 µm, α in = 0.63 m −1 , and q = 9.In (a) the black dotted line shows the free-particle dispersion for comparison.The effect of interaction in modifying the low-momentum part of dispersion from this free-particle tendency is increased as more photons are put in the condensate.At larger momenta all condensates seems to behave as free-particles which their dispersion tends to behave quadratically as free-particles (a), with zero imaginary parts (b).
Zoomed-in (c) real and (d) imaginary part of the dispersion at very low momentum.For this range the real part behaves linearly and the imaginary part is small, suggesting a superfluid behavior.both for the delay and the non-locality of the thermo-optic effective photon interaction.For comparison, the dashed red lines give the critical wavevector when only the non-locality of this interaction, following eq.16, is considered.The lower panels (d)-(f) show the behavoir of the critical velocity (v c ) for corresponding experimental parameters, where both the non-locality and the temporally delayed features of thermo-optic interaction are considered.As described in the text, by increasing the interaction strength either via more absorption (a,d) or larger photon numbers (c,f), the interaction affects quasi-particles up to larger momenta, as seen from the larger critical wavevector as well as the larger critical velocity.Panels (b,e) show that the use of a larger cavity length, which makes the interaction more non-local, weakens the interaction effect, and reduces both the critical wavevector and the critical velocity.

Fig. 2 .
Fig. 2. In the absence of the thermo-optic interaction the Gaussian condensate radius is r BEC = 1 /2m ph Ω ≃ 6 µm for the quoted values of the cavity length and mirror radius of curvature.Panel (a) shows the variation of the condensate radius from the interaction-less case ∆r and the chemical potential ∆µ as a function of number of photons in the condensate.As can be seen both of these parameters linearly increase as the photon number becomes larger, implying a larger condensate with a higher energy.This is consistent with the physical consequence of a repulsive interaction mediated by thermo-optic non-linearity.To clarify the behavior further, Fig.2(b)shows the maximum value of the temperature increase in the dye microcavity as a function of condensate photon number N BEC .As can be seen the temperature monotonically increases with increasing number of photons, giving rise to a larger change of the refractive index, hence a stronger effective photon interaction.

Figures 4 (
Figures 4(a),(b) show the real and imaginary part of the dispersion respectively for various number of photons in the condensate N BEC .The thick, dotted black line in panel (a) gives the (quadratic) free-particle dispersion in the absence of interactions.As stated above, at large momenta this free particle behavior is approached, and the imaginary part asymptotically approaches zero (Fig. 4(b)), leading to a stable condensate with a quadratic dispersion 29 .The low-momentum behavior however, noticeably deviates from the non-interacting dispersion.The difference between the interacting photon fluid dispersion and the ideal one increases for larger photon numbers.Figure4(c) and (d)show the zoomed-in real and imaginary part of Bogoliubov dispersion at very low momenta.The linear dispersion behavior of Fig.4(c), accompanied with low imaginary value as in Fig.4(d), means that the low momentum quasi-particles behave like phonons, and move with a constant velocity v c in the photon condensate.This feature suggests that the photon BEC can potentially be a (a) and (c) show that the critical wavevector increases with both inelastic absorption coefficient α in and photon number N BEC with corresponding dependency √ α in and √ N BEC , respectively.As shown in Fig. 5(b) the critical wavevector decreases for larger mirror separations L, when the interaction becomes more non-local.According to Fig. 4(c), the thermo-optic interaction leads to a linear dispersion at low momenta, indicating the existence of sonic modes in this regime, which fulfills the Landau's criterion of superfluidiy.The slope of this line at low energies determines the velocity of sound in the condensate.Figure 5(d)-(f) shows the dependency of this critical velocity v c 27 C. V. Bindhu, S. S. Harilal, V. P. N. Nampoori, and C. P. G. Vallabhan, Modern Physics Letters B 13, 0178 (1999).

2 FIG. 1 .FIG. 2 .
FIG.1.Schematics of the dye-filled microcavity setup.The mirrors M 1,2 are spherically curved with radius R, leading to a harmonic trapping potential for the photon gas.The optical cavity imposes a non-zero cut-off frequency for the trapped photons, with the longitudinal modal quantum number q.The remaining transverse modal quantum numbers make the photon gas two-dimensional.The photon gas thermalizes by repeated absorption re-emission processes to the spectral temperature of the dye molecules.

1 FIG. 3 .FIG. 4 .
FIG.3.Variation of the Green's function G N L ( r, r ′ ) vs. ρ (the radial parameter) in the transverse plane when r ′ = 0, z = 0.The effect of mirror separation and non-locality are compared for different L. As can be seen the Green's function extends further with an increases in the cavity length.This means that the interaction range is longer in thicker cavities.For large mirror spacing, the Green's function is in good agreement with the gravitational-type dependency given by 1/| r − r ′ | (black dashed line).