Heating and phase-space decompression of evanescent-wave cooled atoms by multiple photon reabsorption

We show that multiple reabsorption of resonance-frequency photons in a cloud of evanescent-wave cooled atoms can have a significant influence on the cooling efficiency and maximum value of the atomic phase-space density. © 2003 Optical Society of America OCIS codes: (020.7010) Trapping; (020.1670) Coherent optical effects References and links 1. J. S̈oding, R. Grimm, and Yu. B. Ovchinnikov, “Gravitational laser trap for atoms with evanescent-wave cooling,”Opt. Commun.119, 652-662 (1995). 2. Yu. B. Ovchinnikov, I. Manek, and R. Grimm, “Surface Trap for Cs atoms based on Evanescent-Wave Cooling,” Phys. Rev. Lett. 79, 2225-2228 (1997). 3. M. Hammes, D. Rychtarik, V. Druzhinina, U. Moslener, I. Manek-H önninger, and R. Grimm, “Optical and evaporative cooling of caesium atoms in the gravito-optical surface trap,” J. Mod. Opt. 47, 2755-2767 (2000). 4. P. Desbiolles, M. Arndt, P. Szriftgiser, and J. Dalibard, “Elementary Sisyphus process close to a dielectric surface,” Phys. Rev. A54, 4292-4198 (1996). 5. R. J. C. Spreeuw, D. Voigt, B. T. Wolschrijn, and H. B. van Linden den Heuvell, “Creating a low-dimensional quantum gas using dark states in an inelastic evanescent-wave mirror,” Phys. Rev. A 61, 053604-1-7 (2000). 6. M. Hammes, D. Rychtarik, H.-C. N̈ agel, and R. Grimm, “Cold-atom gas at very high densities in an optical surface microtrap,” Phys. Rev. A 66, 051401-1-4 (2002). 7. M. Hammes, D. Rychtarik, B. Engeser, H.-C. N ägel, and R. Grimm, “Evanescent-Wave Trapping and Evaporative Cooling of an Atomic Gas at the Crossover to Two Dimensions,” Phys. Rev. Lett. 90, 173001-1-4 (2003). 8. H. Gauck, M. Hartl, D. Schneble, H. Schnitzler, T. Pfau, and J. Mlynek, “Quasi-2D Gas of Laser Cooled Atoms in a Planar Matter Waveguide,” Phys. Rev. Lett. 81, 5298-5301 (1998). 9. P. Domokos and H. Ritsch, “Efficient laoding and cooling in a dynamic optical evanescent-wave microtrap,” Europhys. Lett.54, 306-312 (2001). 10. M. Gorliki, S. Feron, V. Lorent, and M. Dukloy, “Interferometric approaches to atom-surface van der Waals interactions in atomic mirrors,” Phys. Rev. A 61, 013603-1-9 (2000). 11. C. Henkel, K. Mølmer, R. Kaiser, N. Vansteenkiste, C. I. Westbrook, and A. Aspect, “Diffuse atomic reflection at a rough mirror,” Phys. Rev. A55, 1160-1178 (1997). 12. V. Savalli, D. Stevens, J. Est ève, P. D. Featonby, V. Josse, N. Westbrook, C. I. Westbrook, and A. Aspect, “Specular Reflection of Matter Waves from a Rough Mirror,” Phys. Rev. Lett 88, 250404-1-4 (1997). 13. L. Khaykovich, N. Davidson, “Adiabatic focusing of cold atoms in a blue-detuned laser standing wave,” Appl. Phys. B70, 683-688 (2000). 14. S, Meneghini, V. I. Savichev, K. A. H. van Leeuwen, W. P. Schleich, “Atomic focusing and near field imaging: A combination for producing small-period nanostructures,” Appl. Phys. B 70, 675-682 (2000). 15. A. Lenef, T. D. Hammond, E. T. Smith, M. S. Chapman, R. A. Rubenstein, and D. E. Pritchard, “Rotation Sensing with an Atom Interferometer,” Phys. Rev. Lett 78, 760-763 (1997). 16. M.Özcan, “Influence of electric potentials on atom interferometers: Increased rotation sensitivity,” J. Appl. Phys. 83, 6185-6186 (1998). (C) 2003 OSA 11 August 2003 / Vol. 11, No. 16 / OPTICS EXPRESS 1827 #2485 $15.00 US Received May 20, 2003; Revised July 08, 2003 17. A. Peters, K. Y. Chung, and S. Chu, “Measurement of gravitational acceleration by dropping atoms,” Nature 400, 849-852 (1999). 18. T. Calarco, E. A. Hinds, D. Jaksch, J. Schmiedmayer, J. I. Cirac, and P. Zoller, “Quantum gates with neutral atoms: Controlling colisional interactions in time-dependent traps,” Phys. Rev. A 61, 022304-1-11 (2000). 19. M. D. Lukin, M. Fleischhauer, and R. Cote, “Dipole blockade and Quantum Information Processing in Mesoscopic Atomic Ensembles,” Phys. Rev. Lett 87, 037901-1-4 (2001). 20. M. V. Subbotin, V. I. Balykin, D. V. Laryushin, V. S. Letokhov, “Laser controlled atom waveguide as a source of ultracold atoms,” Opt. Commun. 139, 107-116 (1997). 21. H. Nha and W. Jhe, “Sisyphus cooling on the surface of a hollow-mirror atom trap,” Phys. Rev. A 56, 729-736 (1997). 22. J. Yin, Y. Zhu, Y. Wang, “Gravito-optical trap for cold atoms with doughnut-hollow-beam cooling,” Phys. Lett. A 248, 309-318 (1998). 23. I. Manek, Yu. B. Ovchinnikov and R. Grimm, “Generation of a hollow laser beam for atom trapping using an axicon,” Opt. Commun. 147, 67-70 (1998). 24. D. J. Harris and C. M. Savage, “Atomic gravitational cavities from hollow optical fibers,” Phys. Rev. A 51, 3967-3971 (1995).


Introduction
Recent progress in laser cooling and trapping of neutral atoms has opened up the possibility to study a great variety of phenomena inherent to matter at ultralow temperatures.Among the many methods of cooling, those making use of evanescent optical waves [1][2][3][4][5][6][7] are of particular interest when creating cold atomic samples in close proximity to a solid surface.Such samples may be used for fundamental studies of low-dimensional atomic gases [7][8][9], investigations of near-field properties of material surfaces [10][11][12], and also for applied techniques such as atom lithography [13,14] and acceleration sensing [15][16][17].In addition, microscopic atom traps and guides may be useful as integrated atom-optical devices for quantum information processing [18,19].
The evanescent wave created by total internal reflection of a laser beam on a dielectricvacuum interface can serve as a mirror for atoms, if the laser frequency is detuned to the blue from the atomic resonance.In such a repulsive potential, alkali atoms that comprise two ground state hyperfine levels can efficiently be made to lose kinetic energy in a so-called Sisyphus cooling process by letting the atoms repetitively bounce from the mirror potential under the effect of gravity.In this process, a low-power laser field is needed to repump the atoms back to the initial lower-energy hyperfine state in between the bounces.If the repumping beam is directed downwards, an additional slowing of the reflected atoms will take place.The combined effect of this geometric cooling process together with the Sisyphus mechanism is able to cool the atoms to a temperature corresponding to a few photon recoils [1][2][3].Depending on the number density of the atoms, their phase-space density can reach a high value and allow the sample to be further processed towards quantum degeneracy [3,6,7].
A theoretical model describing the cooling of atoms on a tilted evanescent-wave mirror was introduced by J. Söding et al. in [1] and it has since been applied to calculate the cooling parameters in a variety of trap configurations [1][2][3][20][21][22].This model describes well the cooling of atoms at low densities or in small traps.However, the temperature achieved in experiments with large atomic samples turns out to be much higher than what the model predicts [3].In this paper, we show that this discrepancy can be explained by the heating of the sample through the multiple reabsorption of resonant photons, emitted by the atoms in the repumping process.
In Section 2 of this work, we first briefly recall the theory of the evanescent-wave atom cooling and than calculate the contribution of the multiple photon reabsorption.Section 3 describes an example of the cooling of 133 Cs in a gravito-optical surface trap.In Section 4, we discuss the implication of the reabsorption effect on the cooling of atoms on an evanescent wave.

Evanescent-wave cooling and multiple photon reabsorption
Evanescent-wave cooling of atoms has been experimentally demonstrated by using a gravitooptical surface trap [1][2][3].This trap is created by two blue-detuned laser beams, one of which experiences total internal reflection on a flat dielectric-vacuum interface and the other one, being hollow, defines the horizontal extent of the trap (see Fig. 1).In general, the intensity of the evanescent wave created by total internal reflection decays exponentially above the surface with a decay length of Λ =(λ /4π)(n 2 sin 2 θ − 1) −1/2 .Here λ is the wavelength, n the index of refraction of the dielectric, and θ the angle of incidence of the reflected beam.Considering alkali atoms, we assume the evanescent wave to be tuned above the D 2 -resonance frequency by an amount δ with respect to the lower hyperfine ground state |g1 .The low-power repumping beam can be tuned into resonance with the transition from the upper hyperfine ground state |g2 to the excited state |e (in experiments with Cs it is usually the F = 4 → F ′ = 4 transition) in order to optically pump atoms that have entered |g2 back into the state |g1 .If an atom in the state |g1 enters the repulsive evanescent wave and, near the turning point, makes a transition to the state |g2 (through scattering of an evanescent-wave photon), the energy lost by the atom as it moves up the potential will be larger than that it gains when leaving the field after the reflection.Recycling the process leads to a cooling of the atomic sample at the rate [1] γ Sis ≈ 2 9 where δ hfs is the ground-state hyperfine splitting and q e is the mean branching ratio to the lower hyperfine ground state for the elastic scattering of a photon.The parameter τ c denotes the mean time between incoherent reflections.In the case of gravitational confinement of the atoms on the evanescent-wave mirror, this parameter is given by where Γ is the natural width of the atomic line, m the atom's mass, g the acceleration of gravity, and ϕ the angle between the vertical axis and the vacuum-dielectric interface.In contrast to Ref. [1], the factor sin ϕ is in the denominator and, therefore, when ϕ approaches zero, the time τ c goes to infinity, since in this case gravity can no longer return the reflected atoms back to the evanescent wave.Equation ( 2) is derived from the equation of conservation of energy, similar to that considered in Ref. [1], but containing also the potential energy due to gravity.The rate of the geometric cooling, originating from the absorption of the downwards propagating repumping photons by the |g2 -state atoms on their way up in the gravitational field, is calculated to be where q r is the mean branching ratio to the state |g1 for transitions in the repumping field starting from the state |g2 , k B the Boltzmann's constant, and T the instantaneous temperature of the atoms.We point out that, in our calculations, the rate γ geo in Eq. ( 3) turns out to be twice that given in equation ( 13) of Ref. [1].
Assuming that in each optical transition the kinetic energy of an atom increases by one recoil energy, one can derive an expression for the corresponding heating rate (see Ref. [1]): In thermal equilibrium, the heating rate is equal to the overall cooling rate, and the solution of this equation with respect to T will give the final temperature T eq .Since Eq. ( 4) does not contain the contribution of the heating caused by the reabsorption of photons emitted by the atoms in the repumping process, T eq can be considered to be correct for low atomic densities.
In order to extend the theory to include large and/or dense atomic samples, we apply the law of energy conservation and take into account multiple reabsorption events.As before, all the relevant quantities are considered in terms of their mean values per atom.We start by calculating the mean rate at which the resonant photons originally appear in the atomic cloud.In between the reflection events, the atoms spend most of the time outside the strong evanescent field and, due to repumping, they will be found in the state |g1 [1][2][3].Thus, photons that are scattered from the evanescent wave or emitted in spontaneous transitions from the excited state to the state |g2 pass the sample with a negligibly small probability to be absorbed or Raman scattered, since they are far enough from resonance with the transition |g1 →|e .As each cooling cycle ends with a transition |e →|g1 , there will be resonant photons produced with a rate R 1 = (1 − q e )/τ c .The higher the density of atoms in the cloud, the higher is the probability that these photons will be absorbed by the atoms.These photons are created in spontaneous transitions and are emitted in all directions with equal probability.Since they originally appear in a very thin layer above the evanescent-wave mirror, only a fraction of them, η in ≈ 1/2, enter the sample.The mean rate at which the resonant photons get into the sample is, therefore, R in = η in R 1 .
The total mean rate R rad of resonant-photon emissions per atom is equal to the sum of R in and R r , where R r is the additional rate due to multiple reabsorption events.The mean power of the resonance-frequency radiation emitted by each atom is, therefore, P rad = 2π hν 1 R rad , where ν 1 is the resonant-photon frequency.A part of this radiation continuously escapes the sample at a power of P out = η out P rad , where η out is smaller than 1.Since each absorption of a resonancefrequency photon is followed by emission of another such photon, energy conservation requires that the radiation power entering the sample, P in ≡ 2π hν 1 R in , has to be equal to the escaping power P out .Thus we obtain Supposing that η out is known, we now calculate the rate of heating of the atoms in the reabsorption events.This rate is given by γ r =(∆T r /T )R r , where ∆T r is the temperature change caused by optical recoils during the time 1/R r .On average, each atom in this time absorbs and emits one resonant photon, in total undergoing a mean number of 2/q r transitions between the ground and the excited state (due to the presence of the repumping field).Each of the transitions is associated with an increase of the atomic kinetic energy by one recoil energy.Thus, ∆T r =(2/q r )(2E rec /3k B ).By making use of Eq. ( 6) and writing R in and E rec in their explicit forms, we obtain γ r = 1 − q e q r (2π h) 2 3λ 2 mk B T The overall heating rate including the contribution of the resonant-photon reabsorption is then γ heat = γ heat + γ r , where γ heat is given by Eq. ( 4).Taking into account the fact that η in ≈ 1/2, we obtain γ heat ≈ 2 + 1 − q e q r η out Compared with Eq. ( 4), this equation contains the parameter η out in the denominator of the second term.Since η out is always smaller than unity, the rate γ heat is always higher than γ heat , and only if resonance-photon reabsorption can be neglected is η out ≈ 1 and γ heat ≈ γ heat .The thermal equilibrium allowing for the reabsorption effect is obtained by equating γ heat = γ Sis + γ geo .This can be written in the form showing the explicit dependence of the equilibrium temperature T eq on the parameter η out where and The only unknown quantity in Eq. ( 9) is η out ≡ P out /P rad , which may itself depend on T eq .To calculate η out , we assume that the atomic sample has a constant density equal to the peak density n 0 within the effective trap volume.Then, by modelling the emitted resonance-frequency photons as spherical waves, continuously radiated by each atom at power P rad , we can estimate the mean power P out and thus find η out .In these calculations, the attenuation of the waves with the distance of propagation has to be taken into account.The attenuation coefficient for the weak resonance-frequency radiation is given by α = σ res n 0 , where σ res ≈ 3λ 2 /2π is the atomic absorption cross section corresponding to a lifetime broadened transition.The procedure for calculating η out is outlined in Section 3 assuming the geometry of a horizontally aligned gravito-optical surface trap.

Cooling in a gravito-optical surface trap
We apply the model of Section 2 to the particular case of cooling of alkali-metal atoms in a gravito-optical surface trap.This kind of a trap is produced by separating a part of a horizontal evanescent-wave atom mirror by a vertically aligned blue-detuned hollow laser beam [2,3,23].The atoms in the trap are confined in the vertical direction by the evanescent wave and gravity, and in the horizontal directions by the surrounding potential barrier of the hollow beam.The hollow-beam diameter D is usually on the order of 1 mm and, therefore, when the atoms are cooled down, the effective vertical size of the sample, u = k B T /mg, turns out to be much smaller than D. In this case, the calculation of the parameter η out becomes particularly simple, since D can be extended to infinity.The atomic sample can thus be considered to be in the form of an infinite slab with thickness u.Applying the calculation technique described in the previous Section, the average radiated power per one atom leaving the sample through the slab surfaces is obtained in the form where z denotes the coordinate along the normal to the slab and Γ(0, αu) is the incomplete gamma function.Equation ( 12) can in fact be cast in an approximate form, which, while still being quite accurate, allows a more transparent interpretation, In the case of no absorption, αu approaches zero and P out approaches its maximum value of P rad , as it should.On the other hand, for strong absorption, αu is much larger than 1 and P out approaches zero.The attenuation coefficient α = 3λ 2 n 0 /2π depends on the atomic density n 0 = N/(uπD 2 /4) with N being the total number of atoms in the trap.The product αu is, therefore, independent of u and Eq. ( 13) can be rewritten in terms of N as where β = 6λ 2 /(πD) 2 .We remind that this equation is valid for u ≪ D. Substituting Eq. ( 14) into Eq.( 9), we obtain an expression for the equilibrium temperature T eq , which will now be a function of N and D.
As an example, we consider the cooling of a sample of 133 Cs in a trap that was recently used in experiments by M. Hammes et al. [6,7].The inner diameter D of their hollow beam was 0.8 mm and the evanescent-wave detuning was δ = 2π × 5 GHz.For 133 Cs the ground-state hyperfine splitting is δ hfs = 2π × 9.2 GHz, the spontaneous decay rate Γ = 2π × 5.3 MHz, the resonance wavelength λ = 852 nm, the mass m = 2.2 × 10 −25 kg, and the branching ratios are q e = 0.75 and q r = 0.611 [1].Neglecting first the resonant-photon reabsorption, we calculate the equilibrium temperature of the atoms with the aid of Eq. ( 5) by setting ϕ = π/2, and obtain a value of T eq = 1 µK, independent of the number of the trapped atoms.The lowest temperature which has been achieved so far in the experiments with low-density atomic samples of 133 Cs is two times higher than this calculated value [3].One of the possible reasons could be the "roughness" of the realistic evanescent-wave mirror which leads to a diffuse atomic reflection [11] and to a decrease of the geometric-cooling efficiency (see sinϕ in Eq. ( 3)).We now include the reabsorption processes and numerically solve Eq. ( 9), using Eq. ( 14) for η out .The temperature T eq calculated in this way is shown as a function of N in Fig. 2(a) ( T eq (N) for small values of N is shown in the inset).The resonance-photon reabsorption is seen to significantly affect the final temperature of the cooled atoms.The lowest possible temperature in cooling of, e.g., 10 8 atoms turns out to be almost two orders of magnitude higher than T eq given by Eq. (5).At large values of N, which provide small η out , the temperature T eq shows a linear dependence on N with a slope given by as can be seen from Eqs. ( 9) and ( 14).We note that the linear dependence has been observed in experiments of Ref. [3], although the number of atoms was only N ≤ 4 × 10 6 .If we substitute δ = 2π × 3 GHz and D = 0.52 mm as applied in Ref. [3] into Eq.( 15), we obtain a ≈ 2 × 10 −12 K.This compares well with the measured value of 1.5 × 10 −12 K.The density n 0 of the atoms scales as N/ T eq (n 0 = 4N/uπD 2 ) and it is plotted as a function of N in Fig. 2(b).The function n 0 (N) has a local maximum at N ≈ 8 × 10 6 , since at that point the relative increase of temperature, d T eq / T eq , starts to exceed the relative growth of N. When N increases further, the density n 0 slowly drops towards a constant value.This value can be determined to be n ∞ = 4mg/ak B πD 2 ≈ 3 × 10 17 .However, for N > 10 8 the height u of the sample starts to exceed the diameter D and, therefore, Eq. ( 14) cannot be used anymore.
One of the important parameters characterizing the cooling efficiency is the phase-space density.For a fully unpolarized cloud of Cs, it can be expressed in terms of n 0 and T eq as Ω =(n 0 /7)(2π h2 /mk B T eq ) 3/2 .Figure 3a shows the dependence of Ω on N for N < 2 × 10 7 .
The maximum value of Ω is reached at N ≈ 1.8 × 10 6 .We note that in Ref. [6], where the parameter values are as in our calculations, the number of cooled atoms is in fact very close to this calculated optimum value (N = 2 × 10 6 in Ref. [6]).Our simple model shows that Ω starts to drop fast beyond this optimum value as a consequence of the multiple photon reabsorption.
Our model indicates that it should be possible to increase the maximum value of Ω by making the effective size of the trap smaller, since then the resonant radiation would escape the sample more efficiently.In the trap considered, the height u of the sample is much smaller than the diameter D, when N corresponds to a high value of Ω.Therefore, changing D does not result in a significant change of the maximum achievable phase-space density (Ω max ).If, on the other hand, there would be some mechanism that would reduce the height u, the value of Ω max could be increased.As an example, we show the dependence of Ω on N for atoms acted on by an additional force f = 10 × mg (Fig. 3b).This force can be created, e.g., by a spatially inhomogeneous far-red-detuned laser field, similar to one used in Ref. [7], but which could be applied during the cooling.

Discussion
We have shown that in the cooling of an atomic sample by an optical evanescent wave, the presence of photons at the resonance frequency of the atoms leads to a crucial dependence of the temperature and phase-space density of the atomic sample on the number of atoms and on the effective size of the sample.This dependence is due to the heating caused by the process of multiple photon reabsorption.In particular, the phase-space density versus the number of atoms in the trap shows a well defined maximum.This maximum can be increased by reducing the trap size, since the trapped resonant radiation escapes smaller samples more efficiently.Of course, the trap cannot be reduced to an arbitrarily small size for a fixed number of atoms, since at higher densities, the trap lifetime becomes shorter due to inelastic interatomic collisions.The cooling results could be improved further, if, for example, the repumping mechanism was replaced with a process that would not involve spontaneously emitted resonance-frequency photons.Such repumping could in principle be realized through Raman transitions in a field which does not contain the near-resonance frequency components.
The extension of the existing theory to include large and dense atomic samples is important, since high values of the phase-space density occur at such number densities at which the photon reabsorption has already a high rate.Although the model implies cooling of atoms on a planar surface, it can readily be adopted to more complicated traps based on blue-detuned optical fields [20,21,24].The mean cooling and heating rates can in such cases be calculated by considering the atoms as a sample of ideal gas that obeys the Maxwell-Boltzmann density distribution.For traps having comparable sizes in all directions, the calculation of η out will be more complicated than that given in Section 3 for the gravito-optical surface trap.
The theoretical approach described to find the reabsorption contribution can also be applied to calculate the cooling parameters for atoms confined on the evanescent wave by a far reddetuned laser beam.In this case, the effective size of the trap and, as a result, the number of atoms which could be trapped and cooled will be limited by the beam power.In addition, heating of atoms by photon recoils and loss of them due to two-body inelastic collisions (including light-induced collisions leading to molecular states) will have a higher rate.Nevertheless, in the case of a small trap size, the phase-space density can still reach a high level.To find an optimum balance between the required beam power and the detuning, multiple scattering of the photons from the red-detuned laser beam has to be taken into account in addition to the multiple reabsorption.The contribution of the multiple scattering processes to the overall heating rate can be calculated in a similar way by making use of the law of energy conservation.

Fig. 3 .
Fig.3.Phase-space density Ω versus number of atoms N in a gravito-optical surface trap, where a) no external force is applied and b) the applied force is f = 10 × mg.