Classical Rotons in Cold Atomic Traps

H. Ter¸cas,1, ∗ J. T. Mendon¸ca,1, 2 and V. Guerra2 1CFIF, Instituto Superior T´ecnico, Av. Rovisco Pais 1, 1049-001 Lisboa, Portugal 2IPFN, Instituto Superior T´ecnico, Av. Rovisco Pais 1, 1049-001 Lisboa, Portugal We predict the emergence of a roton minimum in the dispersion relation of elementary excitations in cold atomic gases in the presence of diffusive light. In large magneto-topical traps, multiple- scattering of light is responsible for the collective behavior of the system, which is associated to an effective Coulomb-like interaction between the atoms. In optically thick clouds, the re-scattered light undergoes diffusive propagation, which is responsible for a stochastic short-range force acting on the atoms. We show that the dynamical competition between these two forces results on a new polariton mode, which exhibits a roton minimum. Making use of Feynman’s formula for the static structure factor, we show that the roton minimum is related to the appearance of long-range order in the system.

Since the first ideas advanced by Landau [1, 2], the fusion of light, atoms in MOTs experience a complex ef- concept of the “roton minimum” in the dispersion of the fective interaction, which we find to be associated with a collective modes of a certain physical system has played a polariton dispersion relation. This polariton mode is re- central role in the description of superfluidity. After the sult of the dynamical coupling of the density waves with success of the theory in the context of superfluid phases the fluctuations of the light intensity inside the trap. of 4He, rotons have received considerable attention since A route for the most intriguing complex behavior in then and have been identified in many different quantum large magneto-optical traps relies exactly in the multiple interacting systems. Recently, Cormack et al. [3] suggest scattering of light [6, 7]. Due to the consecutive scatter- that rotons may appear in already moderately interact- ing and re-absorption of photons, the atoms experience ing ultracold Bose gases and Kalman et al. [4] have nu- a mediated long-range interaction potential similar to merically observed their emergence in two-dimensional Coulomb system (∼ 1/r) [8, 9] and the system can there- dipolar bosonic gases. fore be regarded as a one-component trapped plasma. In In fact, the emergence of a roton minimum in the ex- a series of previous works, we have put in evidence the citation spectrum strongly depends upon the shape of important consequences of such plasma description of a the interacting potential or, equivalently, on how parti- cold atoms gas [10, 11], whereas the formal analogy and cles are correlated. The correlational origin of the roton the application of plasma physics techniques reveal to be minimum has been firstly suggested by Feynman, where important in the description of driven mechanical insta- the static structure factor S(k) is expressed in terms of bilities [12–16] or even more exciting instability phenom- the dispersion relation ω(k) as ena, like phonon lasing [17]. Moreover, in such optically thick traps, it is known that the light does not propagate ballistically, rather exhibiting a diffusive behavior [18]. k2 In this situation, the energy transport velocity v , i.e. ω(k) = ~ . (1) E 2mS(k) the velocity that accounts for the propagation of energy by the scattered wave, is smaller than c [19, 20]. Labeyrie This has an enormous implication on the interpretation et al. [21] have experimentally observed that vE can, in- of the physical properties of the system in terms of the deed, be several orders of magnitude smaller than c in the dispersion relation: the presence of a roton minimum is case of resonant light propagating in traps, already with a the signature of strong correlations system. In the limit moderate optical thickness, thus putting in evidence the of large mode softening, i.e., for rotons with zero fre- phenomenon of slow light. More recently, the diffusive arXiv:1201.0298v1 [quant-ph] 31 Dec 2011 quency, the system can develop mechanical instabilities, behavior of light has been identified as a source of dy- which can lead to interesting physics phenomena. In Ref. namical instabilities leading to the formation of photon [5], Henkel et al have suggested that the presence of a “ro- bubbles in magneto-optical traps [22]. ton zero” is at the origin of crystallization in ultracold In what follows, we consider that the dynamics of cold Rydberg gases. atoms in MOTs is described by the Vlasov equation In this letter, we describe the classical origin of a ro- ton minimum in the excitation spectrum of cold atomic ! ∂ 1 X clouds confined in magneto-optical traps (MOTs). Due + v · ∇ + F · ∇ f(r, v, t) = 0, (2) ∂t m i v to competition between long-range interactions between i the atoms and the stochastic forces associated to the dif- where f(r, v, t) is the normalized distribution function

Z Z 1 = dr dvf(r, v, t). (3) ∗Electronic address: htercas@cfif.ist.utl.pt 2

P The total force i Fi = FT + Fc accounts for both the trapping and cooling forces. There are evidences [23– β Z 25] that the density profile is approximately constant for δI = 2 dvδf, (10) large traps (typically with N ∼ 109 − 1010 atoms), which iΩ − D0k allows us to consider the system to be homogeneous and where thus to neglect the effects of the trap. The collective force can be described by a Poisson equation [9, 10] 2 2 2∇ I0 2∇ I0 β = 3 2 = D0 (11) Z n0σLτ n0 ∇ · Fc(r, t) = Qeff dvf(r, v, t), (4) is the photon inhomogeneity parameter. Combining the The pre-factor in Eq. (4) represents an effective charge latter result with Eq. (7), we finally obtain the kinetic dispersion relation Qeff = σL(σR − σL)I/c of the atoms induced by light, where σR and σL represent the scattering and absorp- tion cross sections [6–9], and I is the light intensity. For 2   Z ωp ωd 1 ∂f0 most of the experimental conditions, the scattering cross 1 = 1 + dv, (12) k2 iΩ − D k2 v − Ω/k ∂v section is larger than the absorption cross section, i.e. 0 z x σR > σL, enforcing the effective charge to be a positive where we have considered perturbations parallel to the quantity. We have showed that the positiveness of Q is eff wave-vector k = kez, for definiteness. Here, we have an essential condition for the existence of stable oscilla- defined two typical frequencies of the system. The first tions in the system (see e.g. Ref. [10]). one is associated with the oscillations of the atoms due to We now assume that the diffusive behavior of light the long-range force, corresponding to an effective plasma inside the trap can be macroscopically described by the frequency [10] diffusion equation r Q n ∂I ω = eff,0 0 . (13) − ∇ ·D∇I = 0. (5) p m ∂t The second important quantity is the rate at which the The diffusion coefficient is determined by D = `2/τ , photons scatter inside the trap, or simply the diffusion where the photon mean free pass is ` = 1/nσL, with R frequency n = n0 fdv standing for the atomic density. According to experimental results [21], the diffusion time τ can be 2 considered as independent from the atom density, so the βn0 2∇ I0 diffusion coefficient explicitly reads ωd = = D0. (14) I0 I0 We notice that this frequency depends on the scale at  −2 1 1 Z which the diffusive processes occur (micro-, meso- or D(r, t) = 2 2 2 = 2 2 2 f(r, v, t)dv . (6) σLτ n σLτ n0 macroscopic), as it depends upon the spatial scale L at which the light intensity varies. We will discuss the We now linearize the Eqs. (2), (4) and (6) by allowing macroscopic case below. fluctuations around the equilibrium values f = f0 + δf, The integral in Eq. (12) can be evaluated using the I = I0 + δI and D = D0 + δD, such that Landau prescription, according to which the full informa- tion about the initial conditions is cast if the integration path is set to pass below the pole Ω = v k. We split the  ∂  1 z + v · ∇ δf + δF · ∇ f = 0 (7) integral into two parts ∂t m c v 0

Z 1 ∂f Z 1 ∂f 0 dv = Pr 0 dv ∂ 2 2 vz − Ω/k ∂vz vz − Ω/k ∂vz δI − D0∇ δI − δD∇ I0 = 0 (8) ∂t   (15) ∂f0 +iπ dv, ∂v z vz =Ω/k Z δI ∇ · δFc = Qeff,0n0 δf dv + Qeff,0n0 , (9) where Pr stands for the Cauchy principal value. Assum- I0 ing a phase speed vph = ω/k much greater than the width where Qeff,0 = σL(σR − σL)I0/c. Assuming peri- of the distribution, such that f0 and its derivatives get odic perturbations on both space and time, such that small as vz gets large, we may expand the denominator (δf, δI, δD) ∝ exp(ik · r − iΩt), Eqs. (6) and (8) yield in (15) which, together with the relation 3

1.5 2  2 2    ωp usk ωd 1 = 2 1 + 2 1 + 2 Ω Ω iΩ − D0k p 1.0 Ω (19)  2 2 Ω ωpΩ ∂f0 Landau +iπ , k2 ∂v 0.5 damping z vz =Ω/k

√where we have defined the sound atomic speed us = 0.0 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 3vth. Separating the frequency into its real and imag- inary parts, Ω = ω + iγ, with γ  ω, we may finally k ΛD 0 write -2  2  -4 ωdD0k L 2 2 2 2

2 ω = ω + u k 1 − (20)

- p s 2 2 4 - ωp + D0k

10 6 ‰ H

p -8

Ω and  Γ -10

-12 2 2 2 ωd ωp + usk 3 1 −3/(2k2λ2 ) √ D γ = 2 2 4 − 3 3 e , (21) 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 2 ωp + D0k 8π k λD k ΛD where λD = uS/ωp is the effective Debye length. This FIG. 1: (Color online) Illustration of the real (top panel) and dispersion relation describes a quasi-particle excitation imaginary (bottom panel) parts of the polariton dispersion re- resulting from the atom-photon coupling, or a polariton. 2 lation in the macroscopic regime for D0 = 2.0λDωp. We can This is the main result of this paper and we now explicitly observe the emergence of a roton minimum for ωd = 1.9ωp show that is contains a roton minimum. (full black line) and ωd = 1.99ωp. The rotons softens the The diffusive description of light is the result of a (c) frequency to ω(krot) = 0 at the critical value ωd = 2.0ωp. macroscopic approximation which is known to hold if ab- Roton instability is illustrated here for ω = 2.2ωp (red line). sorption takes place at scales much larger than the mean The short-wavelength oscillations corresponding to kλD > 1 free path ` [18]. This is true provided the following hier- are kinematically damped. The usual plasma dispersion rela- archy for the relevant length scales tion of ref. [10] is presented here for comparison (dashed gray line). λ  `  a  L, (22)

where λ is the light wavelength and a is the size of the sys- tem. According to typical experimental conditions [21], Z Z 1 ∂f0 f0 the mean-free path is found to value ` ∼ 300 µm and 2 −1 dv = 2 dv, (16) the diffusion coefficient D ' 0.66 m s . Based on our vz − Ω/k ∂vz (v − Ω/k) 0 z previous estimates [10], the effective plasma frequency and Debye length respectively value ωp ∼ 2π × 100 Hz simply yields and λD ∼ 100 µm. Therefore, provided the identifica- 2 tion L = 2∇ I0/I0 in (14) and using the inequalities in Eq. (22) with a ∼ 1 mm [21], the diffusive approx- Z Z  2 2  1 ∂f0 k vz imation is valid if the diffusion and plasma frequencies Pr dv ' f0 1 + 3 dvxdvydvz. v − Ω/k ∂v Ω2 are of the same order, i.e. ωd ∼ ωp. This is attained if z z p (17) `d ≡ D0/ωp ∼ 1 cm is of the same order of the intensity Assuming the atomic equilibrium to be described by a variation length L, which is a reasonable condition for Maxwell distribution typical experimental scenarios. Nevertheless, by varying the number of atoms in the trap it is possible to control the optical thickness of the system and, therefore, tune the value of the diffusion coefficient D0, which makes our 1 −v2/2v2 f (v) = e th , (18) present estimates even more flexible. In extremis, it may 0 (2πv )3/2 th be also possible to attain different diffusion regimes, but this situation is out of the scope of the present paper. p with vth = kBT/m standing for the thermal speed, we In Fig. (1), it is shown that a roton minimum emerges may finally write in the dispersion relation (20) in the diffusive regime. 4

7 As the values of ωd increase (i.e., for stronger diffusion), the frequency decreases around the roton wavenumber 6 2 krot ∼ λD/`d. This feature is often referred to as mode (c) 5 softening. At the critical value ωd = 2ωp, the mode soft- ens towards zero, which is a clear manifestation of a roton 4 L k

(c) H instability mechanism. For ωd > ωd , the system enters a S 3 crystallization phase. This mechanism has been recently discussed in the literature as it can lead to the forma- 2 tion of supersolids [5]. An important remark is related 1 to the Landau damping at short wavelengths. Modes in > 0 the region kλD ∼ 1 undergo a kinematic damping. For- 0.0 0.5 1.0 1.5 2.0 tunately, rotons are possible to be excited at longer wave- k ΛD lengths (krotλD < 1), thus avoiding the Landau damping mechanism. Moreover, the onset of diffusion tends to de- FIG. 2: Static structure factor S(k) depicted for D0 = crease the damping rate (see Fig. (1)). This nourishes 2 2.0λ ωp. ωd = 0 (gray dashed line), ωd = 1.9ωp (black full hope for rotons to be experimentally observable. D line) and ωd = 1.99ωp (black dashed line). A pick emerges in A remarkable feature of the polariton spectrum in (20) 2 −1 the static structure factor around krot = λD/`d = 0.5λD . is that it exhibits a roton minimum in a three dimensional system even in the absence of strong interactions. It is PV −1 −κ0r clear from the application of Landau’s criterion that the approximated by g(r) ' 1 + c0r cos(k0 + δ0)e , present spectrum does not correspond to that of a su- where c0 is a constant and z0 = κ0 +ik0 is the pole of the perfluid, as the mode is gapped at the origin. This is function S(k)−1. We remark, however, that the PY the- a consequence of the long-range nature of the interac- ory was originally developed for hard-sphere potentials, tion between nearby atoms, and the quasi-particle mass and therefore does not describe systems with long-range is proportional to ωp. This therefore corresponds to an interactions. For that reason, we have not used it to important example where rotons are not intrinsically con- compute g(r). nected to Goldstone modes. Another important property of the classical rotons de- 1.2 scribed above is that they carry useful information about the long-range correlation of the system. By using the 1.1 extension of Feynman’s formula (1) to finite temperature 1.0 systems [27], the static structure factor S(k) is given by L r H g 0.9 k2  ω(k) u2k2 Π S(k) = ~ coth ~ ' s , (23) 2 2 0.8 2mω(k) 2kBT ω(k) krot where we have used the non-degeneracy condition kBT  0.7 0 5 10 15 20 ~ω. In this limit, the structure factor corresponds to that obtained based on a hydrodynamic treatment [28]. In r ΛD Fig. (2), we illustrate the behavior of S(k) for the same FIG. 3: Two-point correlation function g(r) depicted for D0 = parameters of Fig. (1). 2 The static two-point correlation function g(r) = 2.0λDωp. ωd = 0 (gray dashed line), ωd = 1.9ωp (black full 0 2 line) and ωd = 1.99ωp (black dashed line). The correlation hn(r)n(0)i/hn(r)ihn(0)i ' hn(r)n(r )i/n0 can then be function oscillates with the period a of 2π/k ' 4πλ−1. easily calculated provided the relation g(r) = 1 + rot D F −1 [S(k) − 1] [29], which after the integrating out the In conclusion, we have derived the dispersion relation angular variables simply reads for the atom-photon polariton in large magneto-optical traps in the presence of diffusive light. We have explic- 1 Z ∞ k sin(kr) itly computed the excitation spectrum for the particu- g(r) = 1 + [S(k) − 1] . (24) 2 lar case of a thermal atomic distribution, revealing the π 0 r emergence of a roton minimum for a set of parameters As it can be observed in Fig. (3), the appearance of a compatible with current experimental conditions. We minimum in the excitation spectrum (20) is associated have also shown that the increase of the light diffusiv- with the occurrence of long-range correlation in the sys- ity lead to mode softening around the roton minimum, tem. By inspection, one founds that the correlation func- which may eventually drive the system into a roton in- tion oscillates with the period T = 2π/krot. This feature stability situation. Using the relation between the static can be qualitatively understood in the context of Percus- structure factor and the dispersion relation, we have ex- Yevick theory [30, 31], where the correlation function is plicitly demonstrated that the roton minimum is related 5 to the emergence of long-range correlations in the system. number SFRH/BD/37452/2007. This work was partially supported by Funda¸c˜aopara a Ciˆenciae Tecnologia (FCT-Portugal), through the grant

[1] L. D. Landau, J. Phys. USSR 5, 71 (1941). [18] M. C. W. van Rossum and Th. M. Nieuwenhuizen, Rev. [2] L. D. Landau, J. Phys. USSR 11, 91 (1947). Mod. Phys. 71, 313 (1999). [3] S. C. Cormack, D. Schumayer, and D. A. W. Hutchinson, [19] M. P. van Albada, M. P. van Albada, B. A. van Tiggelen, Phys. Rev. Lett. 107, 140401 (2011). A. Lagendijk, and A. Tip, Phys. Rev. Lett. 66, 3132 [4] G. J. Kalman, P. Hartmann, K. I. Golden, A. Filinov and (1991). Z. Donk´o,Eur. Phys. Lett. 90, 55002 (2010). [20] B. A. van Tiggelen, A. Lagendijk, M. P. van Albada, and [5] N. Henkel, R. Nath, and T. Pohl, Phys. Rev. Lett. 104, A. Tip, Phys. Rev. B 45, 233 (1992). 195302 (2010). [21] G. Labeyrie, E. Vaujour, C.A.Muller, D. Delande, C. [6] J. Dalibard, Opt. Commun. 68, 203 (1988). Miniatura, D. Wilkowski, and R. Kaiser, Phys. Rev. Lett. [7] D. W. Sesko, T. G. Walker, and C. E. Wieman, J. Opt. 91, 223904 (2003). Soc. Am. B, 8 946 (1991). [22] J. T. Mendon¸ca and R. Kaiser, Phys. Rev. Lett. (in [8] T. Walker, D. Sesko, and C. Wieman, Phys. Rev. Lett. press). 64, 408 (1990). [23] G. L. Gattobigio, T. Pohl, G. Labeyrie, and R. Kaiser, [9] L. Pruvost, I. Serre, H. T. Duong, and J. Jortner, Phys. Phys. Scr. 81, 025301 (2010). Rev. A, 61, 053408 (2000). [24] T. Pohl, G. Labeyrie, and R. Kaiser, Phys. Rev. A 74, [10] J.T. Mendon¸ca,R. Kaiser, H. Ter¸cas,and J. Loureiro, 023409 (2006). Phys. Rev. A 78, 013408 (2008). [25] H. Ter¸casand J. T. Mendon¸ca, Hydrodynamic equilib- [11] J.T. Mendon¸ca,Phys. Rev. A, 81, 023421 (2010). rium and normal modes in large magneto-optical traps, [12] D. Wilkowski, J. Ringot, D. Hennequin, and J. C. Gar- (in preparation) (2011). reau, Phys. Rev. Lett. 85 1839 (2000). [26] We here consider that the total incident light intensity I [13] A. di Stefano, M. Fauquembergue, P. Verkerk, and D. results from the six laser beams used in standard MOT Hennequin, Phys. Rev. A 67, 033404 (2003); A. di Ste- configurations. fano, P. Verkerk, and D. Hennequin, Eur. Phys. J. D 30, [27] R. P. Feynman, Phys. Rev. 94, 262 (1954). 243 (2004). [28] X. Wang and A. Bhattacharje, Phys. Plasmas 4, 1077 [14] D. Hennequin, Eur. Phys. J. D 28, 135 (2004). (1997). [15] K. Kim, H.-R. Nohm and W. Jhe, Opt. Comm. 236 349 [29] S. S. Z. Ashraf, A. C. Sharma and K. N. Vyas, J. Phys. (2004). Condens. Matter 19, 306201 (2007). [16] H. Ter¸cas,J.T. Mendon¸caand R. Kaiser, Europhys. Lett. [30] J. K. Percus and G. J. Yevick, Phys. Rev. 110, 1 (1958). 89, 53001 (2010). [31] M. S. Wertheim, Phys. Rev. Lett. 10, 321 (1963). [17] J.T. Mendon¸ca,H. Ter¸cas,G. Brodin, and M. Marklund, Eur. Phys. Lett. 91, 33001 (2010).