# Observation of the scissors mode and superfluidity of a trapped Bose-Einstein condensed gas

## Abstract

We report the observation of the scissors mode of a Bose-Einstein condensed gas of <sup>87</sup>Rb atoms in a magnetic trap, which gives direct evidence of superfluidity in this system. The scissors mode of oscillation is excited by a sudden rotation of the anisotropic trapping potential. For a gas above $`T_c`$ (normal fluid) we detect the occurrence of oscillations at two frequencies, with the lower frequency corresponding to the rigid body value of the moment of inertia. Well below $`T_c`$ the condensate oscillates at a single frequency, without damping, as expected for a superfluid.

The relationship between Bose-Einstein condensation (BEC)  and superfluidity has been studied extensively in liquid helium , but only recently has it been possible to examine it in condensates of dilute alkali metal vapours . Helium below its critical point is described by a two-fluid model and the liquid is endowed with a new degree of freedom, namely the relative motion between the normal fluid and the superfluid. “This degree of freedom is the essence of the transport phenomena in He II known collectively as superfluidity” . In the dilute alkali metal condensates various phenomena which imply the occurrence of superfluidity have been observed e.g. collective modes of excitation  and demonstrations of the coherence of the wavefunction . However in a recent theoretical paper D. Guéry-Odelin and S. Stringari  describe how the superfluidity of a trapped BEC may be demonstrated directly and we report the results of such an experiment. Guéry-Odelin and Stringari analyse the so-called scissors mode in which the atomic cloud oscillates with respect to the symmetry axis of the confining potential and they point out that the scissors mode has been used in nuclear physics to demonstrate the superfluidity of neutron and proton clouds in deformed nuclei .

The full theoretical analysis of the scissors mode is given in  and we only outline the key points here. The starting point is a BEC in an anisotropic harmonic potential with three different frequencies $`\omega _x\omega _y<\omega _z`$. The scissors mode may be initiated by a sudden rotation of the trapping potential through a small angle as indicated in Fig.1. In the subsequent motion, the cloud is not deformed provided that the change in the potential is too small to excite shape oscillations. For a thermal gas both rotational and irrotational fluid flow occur in the scissors mode and the normal fluid is predicted to exhibit two frequencies corresponding to those forms of motion. For the BEC there is only irrotational flow because of its single valued wavefunction and therefore it only exhibits one frequency, which is different from either of the frequencies observed for the thermal cloud.

In our experiment the trapping potential is created by a time-averaged orbiting potential (TOP) trap  which is a combination of a static quadrupole field, of gradient $`B_Q^{^{}}`$ in the radial direction, and a time-varying field

$`𝐁_T\left(t\right)=B_r\left(\mathrm{cos}\mathrm{\Omega }t𝐞_x+\mathrm{sin}\mathrm{\Omega }t𝐞_y\right)+B_z\mathrm{cos}\mathrm{\Omega }t𝐞_z.`$   
The term $`𝐁_z\left(t\right)=B_z\mathrm{cos}\mathrm{\Omega }t𝐞_z`$ is additional to the usual field of amplitude $`B_r`$ rotating in the xy-plane. The effect of the additional term $`𝐁_z\left(t\right)`$ is to tilt the plane of the locus of $`𝐁=0`$ by an angle $`\xi =\mathrm{tan}^1(B_z/B_r)`$ with respect to the xy-plane. This causes the symmetry axes of the potential to rotate through an angle $`\varphi \frac{2}{7}\xi `$ in the xz-plane (this analytic result is only valid for $`\xi ^21`$). It is only necessary to consider the xz-plane since the absorption image of the cloud is projected onto this plane. Tilting the locus of $`𝐁=0`$ reduces the oscillation frequency in the z direction from its value when $`B_z=0`$. Thus simply switching on $`𝐁_z(t)`$ also changes the cloud shape and so excites quadrupole mode oscillations. To avoid this we first adiabatically modify the usual TOP trap to a tilted trap and then quickly change $`𝐁_z(t)`$ to $`𝐁_z(t)`$. This procedure rotates the symmetry axes of the trap potential by $`2\varphi `$ without affecting the trap oscillation frequencies (Fig.1).

Our apparatus for producing BEC of <sup>87</sup>Rb is described in  and only a few relevant features are outlined below. We trap and cool the atoms in a small glass cell that forms part of a differentially pumped system with two magneto-optical traps (MOT), one of which is a pyramidal configuration of mirrors inside the vacuum system . The Helmholtz coils which create the oscillating fields along the x, y and z directions are driven by audio amplifiers and it was found to be extremely important to filter out high frequency noise, especially in the coils driving the field along the z-axis . The <sup>87</sup>Rb atoms in the $`F=2,m_F=2`$ state, were probed by a 10 $`\mu `$s pulse of laser light propagating along the y-axis, resonant with the $`F=2,m_F=2F^{}=3,m_F^{}=3`$ transition. The probe pulse was synchronized to the point in the rotation when the magnetic field points along the y-axis and this probing scheme is not affected by $`𝐁_z\left(t\right)`$. We are restricted to probing at times which are multiples of the rotation period of the field (143 $`\mu `$s) but this is very much smaller than the period of the oscillations in the trap, as it must be for the TOP trap to work.

The following experimental procedure was used to excite the scissors mode both in the thermal cloud and in the BEC. Laser cooled atoms were loaded into the magnetic trap and after evaporative cooling the trap frequencies were $`\omega _x=90\pm 0.2`$ Hz and $`\omega _z=\sqrt{8}\omega _x`$. The trap was then adiabatically tilted in $`1`$ s by linearly ramping $`𝐁_z\left(t\right)`$ to its final value, corresponding to $`\varphi =3.6^{}`$ and a reduction of the trap frequency $`\omega _z`$ by $`2\%`$. Suddenly reversing the sign of $`𝐁_z\left(t\right)`$ in less than 100 $`\mu `$s excites the scissors mode in the trapping potential with its symmetry axes now tilted by $`\varphi `$, as described above. The initial orientation of the cloud with respect to the new axis is $`\theta _0=2\varphi `$, so this angle is the expected amplitude of the oscillations (Fig.1). The angle of the cloud was extracted from a 2-dimensional Gaussian fit of the absorption profiles such as those shown in Fig.2.

For the observation of the thermal cloud the atoms were evaporatively cooled to 1 $`\mu `$K which is about 5 times $`T_c`$, the temperature at which quantum degeneracy is observed. At this stage there were $`10^5`$ atoms remaining with a peak density of $`n_0210^{12}`$ cm<sup>-3</sup>. The scissors mode was then excited and pictures of the atom cloud in the trap were taken after a variable delay. The results of many runs are presented in Fig.3(a) showing the way the thermal cloud angle changes with time. The model used to fit this evolution is the sum of two cosines, oscillating at frequencies $`\omega _1`$ and $`\omega _2`$. From the data we deduce $`\omega _1/2\pi =338.5\pm 0.8`$ Hz and $`\omega _2/2\pi =159.1\pm 0.8`$ Hz. These values are in very good agreement with the values $`339\pm 3`$ Hz and $`159\pm 2`$ Hz predicted by theory ; which correspond to $`\omega _1=\omega _z+\omega _x`$ and $`\omega _2=\omega _z\omega _x`$. We measured $`\omega _x`$ and $`\omega _z`$ by observing the center of mass oscillations of a thermal cloud in the untilted TOP trap and calculating the modification of these frequencies caused by the tilt. The amplitudes of the two cosines were found to be the same, showing that the energy is shared equally between the two modes of oscillation.

To observe the scissors mode in a Bose-Einstein condensed gas, we carry out the full evaporative cooling ramp to well below the critical temperature, where no thermal cloud component is observable, leaving more than $`10^4`$ atoms in a pure condensate. After exciting the scissors mode we allow the BEC to evolve in the trap for a variable time and then use the time-of-flight technique to image the condensate 15 ms after releasing it from the trap. The repulsive mean-field interactions cause the cloud to expand rapidly when the confining potential is switched off, so that its spread is much greater than the initial size. The aspect ratio of the expanded cloud is opposite to that of the original condensate in the trap, so that the long axis of the time-of-flight distribution is at $`90^{}`$ to that of the thermal cloud as shown in Fig.2. However this difference in the orientation does not affect the amplitude of the angle of oscillation. The scissors mode in the condensate is described by an angle oscillating at a single frequency $`\omega _c`$

$`\theta \left(t\right)=\varphi +\theta _0\mathrm{cos}\left(\omega _ct\right)`$ (1)   
Figure 3 (b) shows some of the data obtained by exciting the scissors mode in the condensate. Consistent data, showing no damping, was recorded for times up to 100 ms. From an optimized fit to all the data for the function in Eq.1 we find a frequency of $`\omega _c/2\pi =265.6\pm 0.8`$ Hz which agrees very well with the predicted frequency of $`265\pm 2`$ Hz from $`\omega _c=\sqrt{\omega _x^2+\omega _z^2}`$. The aspect ratio of the time-of-flight distribution is constant throughout the data run confirming that there are no shape oscillations and that the initial velocity of a condensate (proportional to $`\dot{\theta }`$) does not have a significant effect.

These observations of the scissors mode clearly demonstrate the superfluidity of Bose-Einstein condensed rubidium atoms in the way predicted by Guéry-Odelin and Stringari . Direct comparison of the thermal cloud and BEC under the same trapping conditions shows a clear difference in behaviour between the irrotational quantum fluid and a classical gas. Another distinction is the lack of damping in the superfluid. The damping of the classical gas is not apparent from our data because it is not sufficiently dense and the time between collisions is many oscillation periods. However the damping in a thermal cloud can be calculated using the Direct Simulation Monte Carlo method , for a mean density and temperature that are roughly the same as those of the BEC. (Since the condensate and the thermal cloud have different spatial distributions the comparison is only approximate.) The results of such a numerical calculation for a thermal cloud at a temperature of 90 nK and a peak density of $`210^{14}cm^3`$ are shown in Fig.4 and the plot shows that the scissors mode of a classical gas is strongly damped under these conditions. Note that even for such high densities both frequency components still occur in the thermal cloud, giving behaviour that is clearly different from the undamped, single frequency oscillation observed for the BEC. The results of this numerical simulation show that the amplitude of the lower frequency component is much smaller than that of the higher frequency one. The higher frequency tends towards $`\omega _c`$, the same frequency as the condensate, when the density increases so that the hydrodynamic regime is reached (where there are many collisions per oscillation period). However in this regime the damping is so strong that only a few oscillation periods would be observed as shown in .

In the near future we plan to measure the frequency of the scissors mode at finite temperatures in the condensate i.e. where a trapped thermal cloud is present in addition to the BEC. Under these conditions the scissors mode should be damped  in a similar way to the quadrupole oscillations at finite temperature . There is still not good agreement between the observed change in frequency of the m=0 oscillation mode of a condensate at finite temperature and theoretical predictions  and therefore using the scissors mode as an alternative means of probing the interaction between the condensate and thermal cloud is important. We have found that the measurement of the angle of orientation of the cloud in the scissors mode is quite robust and could be applied in situations close to $`T_c`$ where only a small proportion of the atoms are condensed. In this way the scissors mode should enable studies of the relative motion of the thermal cloud and condensate corresponding to motion of the superfluid through a normal fluid.

This work was supported by the EPSRC and the TMR program (No. ERB FMRX-CT96-0002). O.M. Maragò acknowledges the support of a Marie Curie Fellowship, TMR program (No. ERB FMBI-CT98-3077).