# Shaking-induced dynamics of cold atoms in magnetic traps

## Abstract

We describe an experiment in which cold rubidium atoms, confined in an elongated magnetic trap, are excited by transverse oscillation of the trap centre. The temperature after excitation exhibits resonance as a function of the driving frequency. We measure these resonances at several different trap frequencies. In order to interpret the experiments, we develop a simple model that incorporates both collisions between atoms and the anharmonicity of the real three-dimensional trapping potential. As well as providing a precise connection between the transverse harmonic oscillation frequency and the temperature resonance frequency, this model gives insight into the heating and loss mechanisms, and into the dynamics of driven clouds of cold trapped atoms.

###### pacs:

37.10.Gh, 37.90.+j, 02.70.Uu^{1}

## I Introduction

It is common to determine the oscillation frequencies of atoms in a magnetic trap by exciting their motion. A sudden displacement of the trap excites a dipole (centre-of-mass) oscillation of the cloud at the trap frequency, while a sudden compression or de-compression induces quadrupole (width) oscillation of the cloud. In many experiments, the trapped cloud is long and thin, with transverse frequencies of 1 kHz or above. The oscillation amplitude is then typically below the resolution of the imaging system and therefore difficult to observe directly. In that case, the oscillation may be driven, and the excitation of the cloud may be detected through the increase in its length, which results from heating. This exhibits a resonant behaviour, with a maximum temperature at a specific driving frequency, related to the natural transverse oscillation frequency.

This paper presents temperature resonances of magnetically trapped atoms, induced by dipole oscillation, with six different trap frequencies. These traps are formed on an atom chip using the field of a permanently magnetised videotape. A full characterisation of the videotape traps is carried out with particular emphasis on the effects of anharmonicity. Our experimental results are compared with a numerical model that provides valuable insight into the dynamics of the cold atoms moving in the real confining potential of the trap under the important influences of collisions and atom loss.

Several references already consider the dipolar excitation of trapped particles by shaking the trap centre Savard et al. (1997); Gehm et al. (1998); Jáuregui (2001); van Eijkelenborg et al. (1999); Kumakura et al. (2003). The first three are theory papers Savard et al. (1997); Gehm et al. (1998); Jáuregui (2001) that describe heating due to laser noise in harmonic optical traps far from resonance. Ref. Gehm et al. (1998) calculates the evolving energy distribution, accounting for atom loss, in a truncated harmonic trap. Ref. van Eijkelenborg et al. (1999) describes measurements of the resonant frequencies and laser cooling rates for ions in a shaken Penning trap. In Ref. Kumakura et al. (2003) M. Kumakura et al. investigate the excitation of neutral atoms in a shaken cloverleaf magnetic trap. They measure resonances in atom loss and temperature and discuss how different ratios of atomic temperature to effective trap depth result in either heating or cooling after shaking. They use a classical 1D equation of motion, without collisions, to illustrate some features of their experiment.

There have also been discussions of parametric excitation, i.e. modulation of the trap frequency, in various contexts. These include heating and cooling of neutral atoms in optical dipole traps Savard et al. (1997); Gehm et al. (1998); Jáuregui (2001); Gardiner et al. (2000); Friebel et al. (1998); Poli et al. (2001), one-dimensional optical lattices Jáuregui et al. (2001); Friebel et al. (1998); Jáuregui (2001) or magnetic traps Zhou et al. (2007), and measurement of oscillation frequencies for electrons in a Penning trap Paasche et al. (2002); Tommaseo et al. (2004), ions in a quadrupole trap Higaki et al. (2007) or neutral atoms in a MOT Kim et al. (2004, 2005). We do not study parametric resonance here because there is no straightforward way to modulate the frequency of our videotape trap without also shaking its position. The same is true of any permanent-magnet atom trap.

The experiments described here determine both the position and the shape of the dipolar temperature resonances. We also investigate these resonances theoretically with the aid of a three-dimensional numerical model that accounts simultaneously for anharmonicity of the trap, atom loss due to finite trap depth, interatomic collisions and evolution of the collision rate during the excitation. All these aspects of our model go beyond what has been done before. Our model yields good quantitative agreement between experiment and theory, even though we use a rather simplified collision model. We anticipate that our simplifying assumption could be exported to a number of ensemble dynamic problems in order to achieve faster simulations.

## Ii The experiment

Our experiment uses the field gradient of a magnetized videotape to trap the atoms. The tape, lying in the plane, has magnetisation , where . Above the surface, this magnetisation produces a transverse (i.e. in the plane) magnetic field of Hinds and Hughes (1999) with G. The addition of a transverse bias field, , produces an array of zeros in the total transverse field located at a height . On Taylor expanding the field around any of these zeros one finds a quadrupole field in the plane whose orientation depends on the angle of the bias field. The magnitude of the field grows linearly with the cylindrical radius according to , providing a cylindrically symmetrical, linear confining potential for weak-field seeking atoms with magnetic quantum number and -factor ( is the Bohr magneton). With the addition of a further bias field , the trapping potential for transverse displacements becomes approximately harmonic in the region close to the axis, where . In that region, atoms oscillate in the trap with a transverse frequency of

(1) |

where is the mass of the atom. We use Rb atoms in the ground state, for which . The bias field is tuneable over the range G, which varies the distance from the trap to the videotape between and m. Over this range of positions the axial bias drops from to G. These parameters give transverse harmonic oscillation frequencies in the range . The axial bias field has a minimum value near the centre of the videotape and increases roughly quadratically with to form a weak axial trap with a frequency of Hz. In summary, this combination of videotape fields and bias fields produces an array of cigar-shaped 3D traps with strong transverse confinement and weaker trapping along . For more details of videotape traps, see Hinds and Hughes (1999); Llorente García et al. (2010).

The upper graphs in Figure 1 show how the total magnetic field strength varies along the three Cartesian axes through the centre of the trap. These are calculated using a full numerical model of the apparatus, with G, directed along . The trap along repeats every m because of the periodic magnetisation of the videotape, but only one of these is shown as only one trap is used in the experiments. Along , we see an asymmetric trap, with a strong repulsive wall as the videotape is approached and an asymptote far from the tape that is equal to the total bias field strength. Along , there is the weak trapping due to the inhomogeneous axial bias field.

We load one such trap using an experimental sequence similar to that described in Llorente García et al. (2010). The temperature of the atom cloud depends on the bias field and ranges from K at G up to K at G. We then add a modulation field along , to displace the trap by a small distance along . We choose mG, which gives a shaking amplitude of nm over the range of transverse bias fields used. Zooming into the centre of the trap, the three lower graphs in Fig. 1 show the extrema of this modulation. The first two show the (equal) movements along and resulting from our particular arrangement of fields. The third shows the absence of movement along and also illustrates the modulation of the minimum field, due to the variation of with transverse displacement of the trap. Since affects the transverse trap frequency, this could cause parametric heating, but the effect is considered in Llorente García (2008) and found to be negligible.

With initially set at G, the atoms are shaken for and held for s, after which the temperature of the trapped cloud is determined by measuring its density distribution along the -direction. This is done using a CCD camera to record the absorption of resonant laser light Llorente García et al. (2010), taking into account the inhomogeneous Zeeman shift of the trapped atoms. The experiment is repeated for a range of shaking frequencies to map out a resonance curve, shown in the leftmost peak of Fig. 2(a). The next two curves are obtained in the same way with G and G. At still higher transverse bias fields ( and G), the atoms are too close to the videotape to yield clean absorption images. In these three cases the trap is moved away from the surface by lowering the bias to G over s before taking the image. This weakens the trap, thereby cooling the cloud, and it is the lower temperature in this final trap that we plot in Fig. 2(a).

These temperature resonances show that the atoms absorb energy most efficiently near a particular frequency. In the case of a one-dimensional harmonic trap, the physics would be that of a driven, weakly damped harmonic oscillator, whose Lorentzian resonance would be centred on the oscillator frequency with a width given by the collision rate. The actual widths are much greater than the collision rate, but nevertheless, motivated by this thought, we fit a Lorentzian to each resonance curve and plot the centre frequencies (the red squares) as a function of in Fig. 3. For comparison, the line in Fig. 3 shows the frequency of small transverse oscillations, given by Eq. 1, with evaluated at the centre of the trap for each value of transverse bias. The data points lie below this line because (i) the most energetic atoms move out of the region of small where the harmonic approximation of Eq. (1) is valid and (ii) the atoms are also displaced from the centre along , where the increased value of reduces the radial frequency. These effects also produce inhomogeneous broadening of the temperature resonances, making them wider for the hotter clouds at higher bias fields, as seen in Fig. 2(a).

Numerically integrating the equations of motion in the full trap potential, we have calculated the period of oscillation along the direction through the centre of the trap. The inverse of this is plotted as a function of amplitude by the solid lines in Fig. 4, the upper(lower) panel being for the case of G. At small amplitude, the frequency coincides with the harmonic approximation indicated by the dotted line. At large amplitude, the potential approaches that of linear trap, and the frequency tends correspondingly to , indicated by the dashed line. The filled circles(squares) mark the rms radius of the cloud in this direction before(after) resonant heating. These show that the atoms explore the anharmonic region of the transverse trap even before the cloud is heated, and move further into this region after heating. These frequency shifts are larger when the bias field is larger. Figure 5 shows how axial displacement from the centre of the trap lowers the harmonic frequency for small transverse oscillations. The essence of this effect is already captured in Eq. (1) through the dependence of on , but here we show the result of the full numerical model of our experiment. Again, the circles(squares) represent the rms size, this time along before(after) heating. Both of these mechanisms contribute appreciably to the inhomogeneous broadening seen in Fig. 2(a) and the lowering of the trap frequency seen in Fig. 3.

For a more quantitative understanding of the resonances it is necessary to build a dynamical model that allows the atoms to collide. In the next section we develop such a model and use it to simulate the data presented in Fig. 2(a).

## Iii Simulation

### iii.1 The role of collisions

Our goal here is to find a numerical model that is simplified as far as possible, while still reproducing the data of Fig. 2(a). We use the known currents and videotape magnetisation to determine an accurate 3D potential, , for the shaken magnetic trap. We note that the inclusion of gravity has no significant effect because these traps are so strong vertically. The cloud of approximately atoms is represented by an ensemble of 500-5000 point particles, this being sufficient to represent the average properties of the ensemble. The three initial velocity components for each particle are chosen at random from the Maxwell-Boltzmann distribution corresponding to initial temperature , which is the baseline temperature in Fig. 2(a) for that particular trap. Similarly, the initial positions are distributed with a probability density proportional to , where is the Boltzmann constant. The classical equations of motion are then integrated numerically to follow the movement of the particles.

The two curves in Fig. 6(a) show how the energy of the cloud increases with time when the trap having G is shaken at 7.2kHz. These parameters correspond to the peak of the pale blue resonance curve at kHz in Fig. 2(a). A simple model without collisions produces the green (lower) curve in Fig. 6(a). Particles are placed in the trap at time , which remains static for the first s. Then the shaking is switched on and the energy rises rapidly, increasing by over the next s. This is due to the excitation of particles whose transverse oscillation frequency is close to the drive frequency. Once they are sufficiently excited, the anharmonicity moves these particles out of resonance. Because the resonant group has been depleted, there is almost no subsequent energy increase even though the shaking continues until s. The shaking is then switched off leaving the cloud to evolve freely over the last second of the simulation. This behaviour disagrees with the experiment. In reality, the temperature doubles (though not in Fig. 2(a) because there the trap was relaxed before recording the temperature). The discrepancy is removed when we allow the simulation to redistribute the momentum through collisions.

The following very simple model of momentum redistribution is sufficient for our purpose. The thermally averaged atom-atom scattering cross-section is , where is the s-wave scattering length Burke et al. (1998) and is the temperature of the cloud. The average relative velocity between particles is Guggenheim et al. (.). Knowing the mean density of trapped atoms, , we obtain a mean collision rate per atom, . At appropriate time intervals (short compared with the inverse collision rate), the numerical integration is paused and a random number is generated for each particle to determine whether it has had a collision. Ideally, this should account for the local variations of density Bird (1994), which we ignore here in order to keep the model simple. If there is no collision, the atom continues unperturbed. Otherwise, the momentum of the particle is redirected by the collision, according to some angular distribution. For the data presented here, we used that of elastic scattering from an infinitely heavy sphere. However, we find that the results are quite insensitive to the chosen distribution and there is nothing special about this particular one. Given our experimental conditions, the average time between collisions is at the start of the excitation, and this becomes longer as the cloud heats up.

With the momentum redistribution thus incorporated, we obtain the dark blue (higher) curve in Fig. 6(a). Now, the ensemble is able to continue absorbing energy after the initial absorption because the depleted velocity group is steadily replenished through collisions. The refilling rate slows down as the atoms become more energetic, and this is responsible for the saturation of heating, seen in Fig. 6(a). In this case, the energy doubles over the s of shaking and our model approximates well the measured heating of the atom cloud.

Figure 6(b) plots the number of atoms in the trap as a function of time, with collisions (blue, lower) and without (green, higher). In the collision-free case, few atoms leave the trap because the heating is weak and because the phase space is not efficiently sampled in the absence of collisions. By contrast, when collisions are included, the cloud heats much more strongly and the energetic atoms are more easily able to find an exit route from the trap. There is competition between the heating rate due to resonant excitation and the cooling rate due to evaporation from the trap. In this example, the heating is dominant because the trap is deep in comparison with the mean energy absorbed by each atom.

Figure 7(a) shows calculated energy spectra in the absence of collisions. The ordinate is the probability density for a given energy, normalised to unity, while the abscissa shows that energy, normalised to the trap depth mK. The red curve (s), showing the initial distribution just before the trap starts to shake, has a single peak just below . After only ms of shaking (green curve, s) a deep notch appears in the distribution close to the energy of the initial peak. This shows that atoms close to that energy, having an oscillation period close to the period of the drive, are the ones excited by the shaking. Their excitation causes a second peak in the distribution at approximately . There is no further significant change in the distribution, even after s of shaking, as shown by the blue curve. This behaviour is to be compared with Figure 7(b), which shows energy distributions for the same simulated experiment when collisions are included. The initial distribution is the same, as is the notch appearing at s, but in this case continued shaking does produce additional heating because the collisions refill the velocity group that absorbs energy from the drive. This is clearly seen in the growing probability on the high-energy end of the spectrum.

### iii.2 Simulated Resonances

Figure 2(b) shows our simulated temperature resonances at each of the six different values of the bias field. After the shaking stops in the experiment, we allow the atoms to thermalise for one second before measuring the temperature. In the simulation, we simply take the final energy as a measure of the temperature that would be reached in equilibrium. The central frequencies of these simulated resonances agree very closely with the experiment, as indicated by the blue circles in Fig. 3.

The peak temperature rises are also remarkably well reproduced by the simulations, given the simplicity of the model. In particular, these model collisions redistribute momentum, but do not permit re-thermalisation of the energy distribution because the energy of a given atom is conserved in the collision. Despite that, the energy increase in the model reproduces all the measured temperature rises to well within a factor of two and that remains the case for a variety of model angular distributions.

The widths of the simulated resonant peaks are a little too large - by a factor of . In searching for an explanation we simulated the yellow (G) resonance with the initial temperature reduced from K to K. This only reduced the width by , so we do not think a temperature calibration error can explain the discrepancy between the measured and calculated widths. It could well be that our very simple collision model causes the resonances to be too broad, although we do not see a clear reason why that should be so.

Figure 8 shows the simulated resonances in atom loss for the same six bias fields used in Fig. 2. We see that the dips in atom number are on the low-frequency side of the temperature resonances indicated by arrows. This is because the energetic atoms most likely to be driven out of the trap are also those most able to explore the anharmonic regions and hence to oscillate at lower frequencies. The same behaviour has been reported by several groups for both parametric shaking Jáuregui et al. (2001); Poli et al. (2001); Roati et al. (2001); Zhou et al. (2007); Barrett (2002), and for shaking of the trap centre Kumakura et al. (2003), in agreement with the results of our simulations.

## Iv Summary and Conclusions

We have measured the temperature rise in a cigar-shaped cloud of cold atoms after shaking it sinusoidally in the transverse direction. Unlike most previous measurements, we have modulated the position of the trap, not its curvature. We have recorded temperature resonances as a function of frequency for several values of the bias field that controls the curvature of the transverse trapping potential. Essential to the interpretation of our measurements is the understanding that, at the temperatures involved, atoms oscillate with a wide range of transverse frequencies in the trap and hence that the central frequencies of the observed resonances lie below the calculated harmonic frequencies for small oscillations.

We have developed a simple numerical model that has provided a clear quantitative understanding of the driven dynamics. When compared with other models and simulations in the literature, this provides one of the most detailed and complete attempts to reproduce the observed resonances. No other models include both the collisions between particles and the anharmonicity of the trapping potential, yet both of these are shown to be essential for a full understanding of the behaviour.

###### Acknowledgements.

The authors thank Robert Nyman and Michael Trupke for useful discussions. We are indebted to the FastNet and AtomChips European networks and to the UK EPSRC and Royal Society for their funding support.### Footnotes

- preprint: APS/123-QED

### References

- T. A. Savard, K. M. O’Hara, and J. E. Thomas, Phys. Rev. A 56, R1095 (1997).
- M. E. Gehm, K. M. O’Hara, T. A. Savard, and J. E. Thomas, Phys. Rev. A 58, 3914 (1998).
- R. Jáuregui, Phys. Rev. A 64, 053408 (2001).
- M. A. van Eijkelenborg, K. Dholakia, M. E. M. Storkey, D. M. Segal, and R. C. Thompson, Optics Communications 159, 169 (1999).
- M. Kumakura, Y. Shirahata, Y. Takasu, Y. Takahashi, and T. Yabuzaki, Phys. Rev. A 68, 021401(R) (2003).
- C. W. Gardiner, J. Ye, H. C. Nagerl, and H. J. Kimble, Phys. Rev. A 61, 045801 (2000).
- S. Friebel, C. D’Andrea, J. Walz, M. Weitz, and T. W. Hänsch, Phys. Rev. A 57, R20 (1998).
- N. Poli, R. J. Brecha, G. Roati, and G. Modugno, Phys. Rev. A 65, 021401(R) (2001).
- R. Jáuregui, N. Poli, G. Roati, and G. Modugno, Phys. Rev. A 64, 033403 (2001).
- S. Zhou, Z. Xu, S. Zhou, L. Liu, and Y. Wang, Phys. Rev. A 75, 053414 (2007).
- P. Paasche, T. Valenzuela, D. Biswas, C. Angelescu, and G. Werth, Eur. Phys. J. D 18, 295 (2002).
- G. Tommaseo, P. Paasche, C. Angelescu, and G. Werth, Eur. Phys. J. D 28, 39 (2004).
- H. Higaki, K. Ito, R. Takai, K. Nakayama, W. Saiki, K. Izawa, and H. Okamoto, Hyperfine Interact 174, 77 (2007).
- K. Kim, H. R. Noh, H. J. Ha, and W. Jhe, Phys. Rev. A 69, 033406 (2004).
- K. Kim, H. R. Noh, and W. Jhe, Phys. Rev. A 71, 033413 (2005).
- E. A. Hinds and I. G. Hughes, J. Phys. D 32, R119 (1999).
- I. Llorente García, B. Darquié, E. A. Curtis, C. D. J. Sinclair, and E. A. Hinds, New J. Phys. 12, 093017 (2010).
- I. Llorente García, Ph.D. thesis, Imperial College London (2008).
- J. P. Burke, J. L. Bohn, B. D. Esry, and C. H. Greene, Phys. Rev. Lett. 80, 2097 (1998).
- E. A. Guggenheim, J. E. Mayer, and F. C. Tompkins, The International Encyclopedia of Physical Chemistry and Chemical Physics (Pergamon Press. Oxford, .).
- G. A. Bird, Molecular Gas Dynamics and the Direct Sim- ulation of Gas Flows (Clarendon Press, Oxford, 1994).
- G. Roati, W. Jastrzebski, A. Simoni, G. Modugno, and M. Inguscio, Phys. Rev. A 63, 052709 (2001).
- M. D. Barrett, Ph.D. thesis, Georgia Institute of Technology (2002).