Superradiant Cooling, Trapping, and Lasing of Dipole-Interacting Clock Atoms

A cold atomic gas with an inverted population on a transition coupled to a field mode of an optical resonator constitutes a generic model of a laser. For quasi-continuous operation, external pumping, trapping and cooling of the atoms is required to confine them in order to achieve enough gain inside the resonator. As inverted atoms are high-field seekers in blue detuned light fields, tuning the cavity mode to the blue side of the atomic gain transition allows for combining lasing with stimulated cavity cooling and dipole trapping of the atoms at the antinodes of the laser field. We study such a configuration using a semiclassical description of particle motion along the cavity axis. In extension of earlier work we include free space atomic and cavity decay as well as atomic dipole-dipole interactions and their corresponding forces. We show that for a proper choice of parameters even in the bad cavity limit the atoms can create a sufficiently strong field inside the resonator such that they are trapped and cooled via the superradiant lasing action with less than one photon on average inside the cavity.


Introduction
The idea of building a superradiant laser operating on an ultra-narrow optical clock transition in a cold gas has fostered the vision of implementing the optical analog of microwave clock masers [1] with a precision and accuracy improved by many orders of magnitude [2][3][4][5]. Today, a central limitation of the best optical clock implementations [6,7] is noise within the mirrors of the reference oscillators [8] that act as the flywheels locked to the atomic transition frequency. When operated on the clock transition in the bad cavity regime and at low photon numbers, superradiant lasers have been predicted to be very insensitive to these fluctuations and create an accurate and precise frequency reference [3,4,9,10].
In principle, operated at high photon numbers sufficiently above the lasing threshold, lasers do not exhibit a fundamental limit of their linewidth [11,12]. In practice, however, the operational laser linewidth is determined by technical noise in the resonator and in the active medium. Technological advances have reduced this limit down to the order of Hz [13], which has lead to a growing interest in using long-lived clock states as the gain medium in a new generation of so-called superradiant lasers [3,[14][15][16]. However, the long lifetimes, i.e. the small linewidths, of those states entail a minute dipole moment of the involved transitions, thus making it necessary to work in the strong collective coupling regime. In this domain, by means of synchronization through the cavity field [17][18][19], a large collective dipole will build up, which can provide the necessary gain. Here, the atoms do not need to be confined in a small volume as is the case with Dicke superradiance [20], but they have to couple almost equally to the cavity.
In the present manuscript we investigate a model of a superradiant laser where the gain medium is self-trapped by the cavity field it creates via stimulated emission into the resonator. At a suitably chosen detuning of the cavity above the atomic transition frequency, the atoms will also slow down and experience a cooling within their prescribed trap positions while simultaneously acting as the gain medium for the laser [21,22]. Recently, very efficient cooling has been predicted involving cavity mediated collective superradiant decay and atomic dipole-dipole interactions [23]. Since inverted atoms in a blue detuned cavity are high-field seekers they are drawn to mode antinodes and almost equal coupling can be achieved.

Model
Let us consider N identical two-level atoms confined to one-dimensional motion along the axis of a Fabry Perot cavity. At finite temperature we can assume a classical description of atomic motion along the cavity axis. All transition dipoles are assumed parallel and perpendicular to the cavity axis as in a J = 0 → J = 1 transition. The atoms couple to the cavity mode via the well known Tavis-Cummings interaction with a strength given by the cavity mode function at the atomic position, g(r i ). Given that the atomic ensemble is closely spaced, we need to take coherent dipole-dipole energy exchange (Ω i j ) and collective spontaneous emission (Γ i j ), which are both mediated by the surrounding vacuum, into account. Furthermore, we assume to create population inversion of the relevant two atomic levels via an individual transverse incoherent pump with the rate R. In practise this has to be implemented via a multistep process involving intermediate levels. Photons can leak through the cavity mirrors at a cavity loss rate κ (see Fig. 1). Fig. 1. Schematic of our model. We consider a confined ensemble of two-level atoms moving along the axis of a cavity resonantly coupled to a single mode with amplitude g(r i ). The atoms directly interact via resonant dipole-dipole coupling inducing pairwise energy exchange Ω i j and collective decay with decay rates Γ i j . A uniform transverse pump mechanism individually excites atoms at rate R, while the cavity loses photons at rate κ.
The Hamiltonian of this system in the rotating wave approximation and in the reference frame of the atoms is where a † (a) is the bosonic creation (annihilation) operator which creates (annihilates) a photon with frequency ω c in the cavity. The operators σ + i and σ − i are the atomic raising and lowering operators of the ith two-level atom with transition frequency ω a . The ith dipole couples to the cavity mode with the position-dependent coupling strength g(r i ) = g cos(k c r i ). The coupling constant is denoted by g and k c = 2π/λ c is the wave number of the cavity mode. The frequency Ω i j quantifies the resonant dipole-dipole energy transfer between atoms i and j. The detuning between the cavity resonance frequency and the atomic transition frequency is given by Dissipative processes are accounted for by the Liouvillian L in the master equation Within the Markov approximation our Liouvillian consists of three parts, namely where the individual incoherent transversal pump is characterized by the pump rate R, the cavity losses occur at the cavity decay rate κ, and the collective atomic decay is determined by the generalized spontaneous emission rates Γ i j , The resonant dipole-dipole couplings Ω i j and the collective decay rates Γ i j depend on the interatomic distances [24,25] and are given by and Here, k a = ω a /c is the wavenumber corresponding to the atomic transition frequency and Θ denotes the angle between the atomic dipoles and the distance vector between atom i and atom j.
For the time evolution of the classical variables we have for the velocity of the ith particle and the force acting on a particle is (see Appendix A: for details) Here, we defined ω r := k 2 a /(2m) as the recoil frequency, with m the mass of an atom. Note that the above equations are only valid for sufficiently slow particles. This is because we include the time dependence of the collective dipole-dipole interactions via the time-dependent atomic positions only. Thus, Doppler shifts and other effects depending on the velocity, or higher-order derivatives of the position, are neglected in the dipole-dipole coupling. Furthermore, we note that forces stemming from the collective decay are neglected here (see Appendix A: ). Additionally, since we assume classical motion, the recoil from spontaneous emission is neglected in the kinetic energy. This is probably the most drastic approximation made here.

Cooling and trapping properties
We investigate the stability of the system described above in the lasing regime by showing that the atoms are cooled and trapped within the cavity field potential created by the photons scattered from the inverted atoms. Due to the exponential scaling of the Hilbert space dimension with the number of atoms numerical methods are limited. Thus, we restrict ourselves to treating a sufficiently small system that still exhibits collective effects.
We study a system with three atoms inside the cavity. The initial state is as follows. The atoms are placed λ c /2 apart at the cavity field antinodes. Furthermore, the atoms are in the ground state and there are no photons inside the cavity. The set of initial momenta is picked from a normal distribution that depends on the recoil frequency. Namely, we always choose a normal distribution for the atomic momentum such that the average kinetic energy is constant in respect to ω r . In general, the kinetic energy of the ith particle is p 2 i /(2m). Thus, if the recoil frequency is multiplied by an arbitrary constant c (i.e. the mass is divided by c), we need to scale the momentum with 1/ √ c in order to keep the kinetic energy constant. The standard deviation p 0 of the initial momentum distribution of the atoms is chosen depending on the choice of ω r according to this relation.
To analyze the trapping and cooling properties of the system we study the time evolution of the particle positions. In Fig. 2(a) we show a case where the particles are cooled until they are cold enough to get trapped in the potential created by the cavity field. The particles distribute themselves relatively far from each other, which means that collective effects hardly play a role. As we aim to investigate collective effects as well, we restrict our calculations to particle trajectories that remain in their initially prescribed trap for the entire time evolution. We call them completely stable trajectories, see for example Fig. 2(b). We refer to Appendix C: for more details on the stability. The momentum transfer from particle 2 to particle 3, depicted in Fig. 2(a), stems from the collective dipole-dipole effects.
In order to quantitatively capture the cooling process, we study the time evolution of the particles' kinetic energy and average over 100 thermally distributed initial momenta. However, we consider the completely stable trajectories only. As evident from Fig. 2(c) (blue line) this kinetic energy of the stable trajectories oscillates very rapidly on the time scale of the cooling process and thus does not yield comparable results. Therefore, we introduce the time-averaged kinetic energyĒ kin (t), which is obtained by taking the midpoints between two adjacent extrema of the kinetic energy as seen in Fig. 2(c) (orange line). The parameter we use in order to characterize cooling or heating is which we call time-averaged relative kinetic energy. We scan over the experimentally most accessible parameters using the procedure described above for three different ω r . The thermally distributed initial kinetic energy corresponds to a normal distribution of the initial momenta. As discussed above, in order to ensure that the particles start with the same average kinetic energy for all values of ω r we scale the standard deviationp 0 of the momentum distribution. We choosep 0 = 2 k a for ω r = 0.1Γ,p 0 = 2/ √ 10 k a for ω r = 1Γ andp 0 = 2/10 k a for ω r = 10Γ. Fig. 3 shows the scan over ∆ and R as well as over g and R for all three values of ω r . Every value ofĒ rel kin (t) above 1.0 corresponds to heating and is artificially fixed to 1.0, as these are points of little interest.
We observe that only trajectories with ∆ > 0 realize cooling, which corresponds to the expected blue detuning of the cavity mode with respect to the atoms. This is due to the fact that atoms inside a cavity favour the emission of photons near the cavity resonance frequency [26]. An atom in a blue detuned cavity emits photons at a frequency higher than its transition frequency. Therefore, the atom has to exert energy in order to lose a photon into the cavity, which it does by losing kinetic energy. The atoms feel an effective friction force that is largest at the points where the cavity field is maximal (high-field seeking). As can be seen from the scans, there is an optimum for the detuning where the cooling is maximal. This is similar to the maximal force in the process of Doppler cooling. The force is also proportional to the excited state population of the atoms. The cooling is thus best when the pump is sufficiently strong to keep the atoms inverted at almost all times.
Furthermore, we can see that atoms with a larger ω r reach a lower relative kinetic energy during a fixed cooling time. On the one hand, this means that lighter particles are cooled down faster. On the other hand, their initially larger velocities make them more difficult to trap, i.e. more trajectories are unstable (see Appendix C: ). Heavier atoms (smaller ω r ) are easier to trap for the same initial kinetic energy (see Fig. 8), even though they do not cool as much during the observed time interval. Note the difference between cooling and trapping here: heavier atoms can still be trapped if they cool poorly, since even a larger kinetic energy in this case oftentimes corresponds to a relatively small velocity insufficient for the particles to climb the potential walls of the trapping potential created by the cavity field. Since they start with a lower initial velocity, however, the cooling is much slower. The inverse line of argument holds for lighter atoms: they are more difficult to trap, but if they are trapped the cooling is more efficient.
Finally, as can be seen in Figs. 3(d)-3(f), the coupling to the cavity mode should not be too large in order for the system to cool the atoms. This can be explained by the growing probability of the atoms absorbing photons from the cavity, which causes heating. Hence, the coupling strength should always be well below the cavity loss rate, such that it is much more probable for a photon to leave the cavity than to be reabsorbed. Higher pump strengths can also counteract the heating. If the atoms are pumped strongly they are inverted at almost all times and thus cannot absorb an incident cavity photon. Note that the red areas in Fig. 3 with ∆ > 0 are mainly caused by extremely slow initial atoms. In these cases the atoms do remain trapped, even though they are heated (note again the difference between cooling and trapping). They are initially so slow that the noise stemming from the cavity field causes heating inside their trap. If the atoms here started with a larger kinetic energy (temperature), they would indeed be cooled. However, they would then also be fast enough to leave their initial traps.
At this point we would like to emphasize again that we describe the system in a semiclassical treatment and we neglect the recoil arising from spontaneous emission. Since the absolute values of the particles' momenta are around k a we need to view the cooling and trapping results critically, especially for ω r = 1Γ and ω r = 10Γ. Note, though, that the recoil heating would only have an effect here since we formulate the superradiant laser regime in terms of a toy model. Specifically, the spontaneous emission rate is taken to be in the limit where Γ κ. However, it is still chosen much larger than it would be for realistic clock atoms in order to avoid numerical difficulties due to different timescales. While we choose Γ ∼ 10 −1 κ, a more realistic choice would be Γ ∼ 10 −6 κ. In that case, a spontaneous emission event is so rare that the recoil can be safely neglected. Still, the fact that we do not include recoil effects for our choice of parameters may be viewed as a rather drastic simplification in our model.

Collective cooling effects
Let us now investigate the relevance of the collective effects for the cooling process. Therefore, we set Ω i j = 0 and Γ i j = δ i j Γ, and compare this independent cooling to the collective cooling from above. In order to acquire collective effects we find that we need to extend the cooling time by orders of magnitude. In Fig. 4 we see the time evolution ofĒ rel kin for collectively interacting atoms in comparison to independent ones.
The main result from Fig. 4 is that independent atoms will always reach a lower final kinetic energy for long cooling times. In the collective case the atoms push or pull each other away from We compare the motional energy loss for collectively interacting (solid lines) and independent atoms (dashed lines) showingĒ rel kin for both cases. The independent case describes atoms far apart from each other. The parameters are N = 3, ∆ = 5Γ, g = 5Γ, κ = 10Γ and R = 8Γ for both, ω r = 0.1Γ and ω r = 1Γ. the cavity field antinodes and thus their displacement amplitude is larger. Therefore, their kinetic energy is bigger on average. Until approximately Γt = 5000 the collective line is slightly below the independent line, for ω r = 0.1Γ. The reason for this is that parts of the kinetic energy are absorbed into the dipole-dipole interaction potential. The fact that there is a minimum below the final value in the collective case stems also from the dipole-dipole interaction. The minimal temperatures in the two collective cases shown in Fig. 4 are reached at approximately the same time in units of ω r .
As we mentioned in the beginning, we restrict our considerations to N = 3 due to the exponential growth of the Hilbert space. Let us still comment on what one might expect for a larger number of atoms in terms of cooling. In [23], it has been shown that efficient cavity cooling can be achieved without direct dipole-dipole interactions between the atoms. Rather, the interactions there stem from the cavity mediated dipole coupling. These findings in combination with the result shown in Fig. 4 suggest that the limit imposed on the final kinetic energy by direct dipole-dipole coupling will be more pronounced for larger atom numbers. More precisely, the cooling without direct dipole-dipole interactions yields lower final kinetic energies with growing atom number [23]. The larger energy due to the displacement caused by direct dipole-dipole interactions should thus cause an increasing difference to non-interacting atoms.

Laser properties
After having established that the system is stable for a given set of parameters, we proceed by analyzing its lasing properties. To this end we study the cavity spectrum, the average photon number as well as the second-order correlation function. Furthermore, we look at the atomic inversion. We use the density matrices at Γt = 500, which describe a quasi-stationary final state, in order to calculate the properties mentioned above.

Laser spectrum
The laser spectrum can be calculated as the Fourier transform of the first order correlation function g (1) (τ) = a † (t + τ)a(t) . According to the Wiener-Khinchin theorem [27], we have

Fig. 5. Properties of the emitted laser light as a function of different system parameters.
We depict scans for the peak frequency shift δ 0 from the atomic resonance in (a) and (b), the laser linewidth γ in (c) and (d), the average photon number n in (e) and (f), and the bunching parameter g (2) (0) in (g) and (h). We focus on values of ∆ ≥ 5Γ, since below this threshold the particle motion shows few stable trajectories only (see Fig. 8).
The remaining parameters are the same as in Fig. 3.
Appendix B: provides details on how the spectrum is calculated in our semiclassical approximation. Most of the spectra are well described by a Lorentzian distribution (see Fig. 7). Thus, we determine the full width at half maximum (FWHM) γ and the offset to the atomic resonance frequency δ 0 . The dependency of the linewidth γ and the offset δ 0 on our scan parameters is depicted in Figs. 5(a)-5(d). We show these plots for one choice of the recoil frequency only, namely ω r = 1Γ, since they are qualitatively identical for the other two choices of ω r . The central observation from Figs. 5(a)-5(d) is that the laser offset δ 0 is much smaller than the corresponding detuning ∆ for all parameters. Mathematically, this means that the slope of the offset's dependency on the detuning (cavity pulling coefficient) is smaller than one. For a conventional laser in the good cavity regime the cavity pulling coefficient is approximately one. In our case it is roughly between 0.1 and 0.2, depending on the pump rate R. In addition, the linewidth γ does not vary much with the detuning. The significance of these two features is that the spectrum of the superradiant laser depends less on cavity fluctuations than the spectrum of a normal laser, which is the expected behaviour. Between ∆ = 5Γ and ∆ = 15Γ the frequency offset grows as expected, but for larger detunings it seems to reduce. The reason for this is that too large detunings lead to the formation of a distinct second peak approximately at the cavity resonance frequency, as shown in Fig. 6. Therefore, in these cases we determine the frequency close to the atomic resonance only, which almost vanishes.
The incoherent drive effectively broadens the atomic transition, resulting in an actual linewidth of R + Γ. We hence plot the FWHM in units of this effective atomic linewidth in Fig. 5(c) and Fig. 5(d). One can see that there are areas where the laser linewidth is even lower than the effective atomic linewidth. This is the case for high pump strengths and small detunings. The smallest laser linewidth, however, is achieved for small pump rates. We can also see that, for all stable parameters, the laser linewidth is well below the cavity linewidth, γ < 2κ. The linewidth and the offset grow with increasing pump strength, which implies that the narrowest laser spectrum featuring a low frequency shift is achieved at small pump rates just above the lasing threshold. The atom-field coupling does not affect the offset, but the linewidth grows with it.
We note that the lasing properties shown in Fig. 5 are almost identical for the three different Fig. 6. Appearance of a second maximum in the spectrum for large atom cavity detuning. For detunings larger than ∆ = 25Γ a second peak emerges at the cavity resonance to the right-hand side of the atomic peak. For ∆ = 50Γ this peak is almost completely separated. This shifts the average of the emitted intensity towards the cavity resonance. We call the offset from the atomic transition frequency δ a and the one from the cavity resonance δ c . The parameters are N = 3, ω r = 0.1Γ, g = 5Γ, κ = 10Γ and R = 10Γ for both, ∆ = 25Γ and ∆ = 50Γ.
ω r , which indicates that the laser properties do not change dramatically compared to a laser with fixed particle positions as described in [4]. To further support this statement we calculate the spectrum in the same manner as before, but for fixed atomic positions (r 1 = 0, r 2 = λ c /2 and r 3 = λ c ). Comparing the resulting spectrum to the one with moving atoms, we find an almost perfect overlap. Therefore, the atomic motion appears to merely change the effective atom-field coupling, which does not significantly alter the spectrum.

Photon number, second-order correlation and population inversion
Besides the spectrum we also calculate other characteristic quantities of a laser. Specifically, we compute the average intra-cavity photon number, and the second-order correlation function at zero time delay, Finally, the population inversion of the atoms is relevant as well. The overall excitation is given by where inversion is achieved if p e > N/2. Figs. 5(e)-5(h) depict n and g (2) (0) as functions of the scan parameters for ω r = 1Γ. The most significant feature is that we always have less than half a photon on average inside the cavity. The figure also shows that the most photons are created and the field is most coherent (g (2) (0) = 1) for small detunings, large pump strengths and large atom-field coupling. This behaviour coincides with that of a conventional (good-cavity) laser. The excited state population is always above 1.5, and we note that the overall scan of the atomic excitation is similar to that obtained for a conventional laser.
Comparing the cooling and the lasing scan, we find that the optimal lasing point does not coincide with the best cooling, specifically for the pump strength dependency. We therefore conclude that there is a certain trade-off between the optimal cooling and lasing regimes.

Conclusions
We have seen that even less than one average intra-cavity photon can be sufficient in order to accumulate excited state atoms dynamically at positions of maximal light coupling, i.e. at field mode antinodes, in the blue-detuned regime. For a sufficient pumping one can thus achieve population inversion and gain which subsequently leads to superradiant lasing. This behaviour is stable with respect to forces and heating induced by dipole-dipole interaction. The output spectrum of such a laser exhibits a very low sensitivity to cavity length fluctuations with a linewidth determined by the atomic linewidth broadened by the pump rate. We have obtained these results by means of a semiclassical model, in which we have treated the atomic states as well as the cavity field mode quantum mechanically, whereas the atomic motion has been described classically.
Overall, for sufficiently slow atoms, the atomic motion only marginally affects the operating conditions and output characteristics of such a laser. In particular, its spectral and coherence properties remain almost unchanged as long as the photon numer is low. This is a promising result for the construction of a superradiant laser, where inverted atoms are moved through the cavity by an optical lattice conveyer belt. It seems that using as many atoms as possible with a weak pump and a large bandwidth cavity is the optimal way to operate such a device. Note, that we have used a rather generic rate based spatially uniform pumping scheme. This should be refined and modeled in more detail for future considerations. (generally 3D) motion of the atoms is accounted for by the position and momentum operatorsr i andp i . The coupling to the cavity mode is determined by the component of the position along the cavity axisr c i = k c ·r i /k c . The free energy part of the cavity and the atoms is given by Note, that we perform the so-called independent bath assumption for the atoms and the cavity, i.e., the cavity decay does not affect the coupling of the atoms to the environment. Since the cavity damping does not affect the motion of the atoms directly either, we will neglect it for now. The density operator describing the internal atomic dynamics as well as the motional degrees of freedom, the cavity mode and the 3D vacuum modes ρ tot is then governed by the von Neumann equation, Essentially, the semiclassical approximation consists of two assumptions. First, we assume that there are no correlations (entanglement) between the motion and the remaining degrees of freedom. Secondly, we will assume that the motion is classical such that all expectation values factorize.
The assumption that there are no correlations between the motion and the remaining degrees of freedom amounts to setting On the one hand, the density operator ρ acf describes the state of the atomic excitation, the cavity, as well as the free-space vacuum modes. On the other hand, the motional degrees of freedom are given by ρ m (t). We now aim at finding an equation for the reduced system density operator ρ acf .
To this end, we take the partial trace, Here, we have defined the reduced Hamiltonian where we wrote Additionally, we made our second assumption of treating the motion clasically, such that f (r) ≈ f (r) for any function f . The assumptions from above constitute our semiclassical approximation. We can now proceed by eliminating the field modes. This leads to the dipole-dipole interactions among the atoms in the form of coherent energy exchange as well as collective decay. The only additional assumption one has to make to arrive at this is the Markov approximation for the atomic positions. The remaining procedure remains the same and yields [24], with H the Hamiltonian from Eq. (1). The motion of the atoms, however, is still determined by H tot via and equivalently for r i . While the velocity is given in Eq. (9), for the average force on the ith particle we have (in 1D) Note that the term proportional to the collective decay does not significantly contribute since in our system Im σ + i σ − j ≈ 0 ∀ i, j due to almost perfect phase invariance [14].

Appendix B: Calculation of the spectrum
The spectrum is given by the Fourier transform of the correlation function A common method to calculate this is to define a new density operator at time t which is then evolved up to a time t + τ. This is also known as the optical regression theorem and the essential steps are as follows. Let H tot be a Hamiltonian which describes the entire system and bath dynamics. The evolution of the system is then reversible and given by the unitary operator U(t) = exp(−iH tot t/ ). Thus, we can write the correlation function as where ρ tot is the total density operator. Upon defining a new density operatorρ tot (0) := aρ tot (t), we may write g(τ) = tr a †ρ tot (τ) .
The time evolution ofρ tot is given by the same unitary operator as before. Therefore, eliminating the bath leads to the same master equation, but for a new operatorρ = aρ. In this way it is possible to compute g(τ) from the reduced system density operator via the master equation.
One has to be careful when deriving the semiclassical master equation forρ, though. In particular, the force on the atoms is proportional to average values of system variables such as a † σ j . These have to be computed from the actual density operator ρ, rather than fromρ. Assuming, as before, that there is no entanglement between the atomic motion and the remaining degrees of freedom [Eq. (20)] is equivalent to writing U(t) ≈ U acf (t) ⊗ U m (t). We can hence writeρ where in the second step we have implicitly definedρ acf (0) := aρ acf (t) and have used the fact that a does not act on the motional degrees of freedom. It is then possible to obtain a semiclassical master equation forρ by tracing out the motion as well as the vacuum modes. However, as can be seen from Eq. (30), the motional degrees of freedom are still determined by ρ m (t + τ). Thus, we need to compute the motion up to a time t + τ in a time evolution with the density operator ρ.
Only then can we calculate the proper time evolution ofρ by using the previously calculated particle positions and obtain the correlation function g(τ).
If the detuning between the cavity and the atoms is not too large, the cavity output spectrum can be well described by a Lorentzian distribution, see Fig. 7. Fig. 7. Lorentzian fit of the normalized spectrum. The calculated spectrum is fitted with a three-parameter Lorentzian function. We call the FWHM γ and the offset to the atomic resonance frequency δ 0 . The parameters are the same as in Fig. 2(b).

Appendix C: Stability
In Fig. 8 we provide a scan of the number of completely stable trajectories, i.e. trajectories where the atoms stay in their initial trap for the whole time evolution. As mentioned in Sec. 3 we see that lighter particles with the same kinetic energy are more difficult to trap. This is simply because they have a higher initial velocity and therefore it is harder for the cavity field to keep them trapped. Also, it takes a certain amount of time for the cavity to build up a sufficiently strong field which can confine the atoms. If the atomic velocity is too large, an atom may leave its initial trap during this build-up time.