Critical velocity for superfluid flow across the BEC-BCS crossover

Critical velocity for superfluid flow across the BEC-BCS crossover


Critical velocities have been observed in an ultracold superfluid Fermi gas throughout the BEC-BCS crossover. A pronounced peak of the critical velocity at unitarity demonstrates that superfluidity is most robust for resonant atomic interactions. Critical velocities were determined from the abrupt onset of dissipation when the velocity of a moving one dimensional optical lattice was varied. The dependence of the critical velocity on lattice depth and on the inhomogeneous density profile was studied.

03.75.Kk, 03.75.Lm, 03.75.Ss

The recent realization of the BEC-BCS crossover in ultracold atomic gases Gior07 () allows one to study how bosonic superfluidity transforms into fermionic superfluidity. The critical velocity for superfluid flow is determined by the low-lying excitations of the superfluid. For weakly bound fermions, the (Landau) critical velocity is proportional to the binding energy of the pairs, which increases monotonically along the crossover into the BEC regime. However, the speed of sound, which sets the critical velocity for phonon excitations, is almost constant in the BCS regime, but then decreases monotonically on the BEC side, since the strongly bound molecules are weakly interacting. At the BEC-BCS crossover, one expects a rather narrow transition from a region where excitation of sound limits superfluid flow, to a region where pair breaking dominates. In this transition region, the critical velocity is predicted to reach a maximum Sens06 (); Stri06 (); Stri07 (). This makes the critical velocity one of the few quantities which show a pronounced peak across the BEC-BCS crossover in contrast to the chemical potential, the transition temperature deMelo_PhysRevLett.71.3202 (), the speed of sound zwer07 (); jose07 () and the frequencies of shape oscillations Grim07 (), which all vary monotonically.

In this paper, we report the first study of critical velocities across the BEC-BCS crossover, where a Feshbach resonance allows the magnetic tuning of the atomic interactions, and find that superfluid flow is most robust near the resonance. Our observation of a pronounced maximum of the critical velocity is in agreement with the predicted crossover between the two different mechanisms for dissipation.

Critical velocities have been determined before in atomic BECs perturbed by a stirring beam Rama99 (); Onof00 (); Rama01 () as well as by a 1D moving optical lattice ingu05 (). In both cases, the inhomogeneous density of the harmonically trapped sample had to be carefully accounted for in order to make quantitative comparisons to theory. Here we mitigate this problem by probing only the central region of our sample with a tightly focused moving lattice formed from two intersecting laser beams. For decreasing lattice depths, the critical velocity increases and, at very small depths, approaches a value which is in agreement with theoretical predictions.

In our experiments, we first create a superfluid of Li pairs according to the procedure described in previous work chin06 (). Forced evaporative cooling of an even mixture of the two lowest hyperfine states is performed at a magnetic field of 822 G, on the BEC side of a broad Feshbach resonance centered at G. This results in a nearly pure Bose-Einstein condensate of pairs in a cross optical dipole trap with harmonic trapping frequencies Hz. The Fermi energy of the system is kHz. To form the moving lattice, we focus two phase-locked 1064 nm laser beams to intersect at the sample with an angle of (see Figure 1). The resulting 1D lattice has a spatial period of m. A frequency difference between the two beams of causes the lattice to move with velocity . The beams have waists of 20 m and 60 m respectively, and address a relatively homogeneous region at the center of the cloud which has Thomas-Fermi radii m. The minimum density at the position of the waist is 42% of the central density.

The lattice which necessarily varies in depth across the sample, is characterized by its peak depth specified in units of or the recoil energy kHz, where is the molecular mass. The lattice depth is calibrated using Kapitza-Dirac scattering. Due to the inhomogeneity of the lattice, the uncertainty is . The lattice depths explored in this work are sufficiently small such that motion induced in the laboratory frame is negligible, in contrast to mun07 ().

Lattice Velocity (mm/s)

Figure 1: Onset of dissipation for superfluid fermions in a moving optical lattice. (inset) Schematic of the experiment in which two intersecting laser beams produced a moving optical lattice at the center of an optically trapped cloud (trapping beams not shown). Number of fermion pairs which remained in the condensate after being subjected to a deep optical lattice for 500 ms, moving with velocity , at a magnetic field of 822 G ()) An abrupt onset of dissipation occurred above a critical velocity , which we identify from a fit to Equation 1.

The lattice moving at a constant velocity is adiabatically ramped up and held for a time up to 2 s, after which the lattice is ramped down and all confinement is switched off. As in previous work chin06 (), a fast magnetic field ramp is used to reduce strong interactions in order to probe the center-of-mass momentum distribution of the pairs. Subsequently, absorption imaging is done on the atomic resonance line at 730 G. A bimodal fit reveals the number of pairs remaining in the condensate , providing a measure of the heating incurred during application of the moving lattice.

Figure 1 illustrates the characteristic dependence of dissipation on the velocity of the moving lattice. At low velocities, the sample is unaffected. Above some critical velocity , dissipation sets in abruptly. We determine from a fit of to the intersection of two lines with slopes 0 and :


The critical velocity which we obtain from this procedure is consistent for a large range of hold times, varying by less than 15% when the hold time is changed by a factor of 20. We explore the BEC-BCS crossover by adiabatically ramping the magnetic field to different values after evaporation and repeating the measurement as before. The crossover is parameterized by the interaction parameter , where is the Fermi wavevector and is the B-field dependent s-wave scattering length bart05 (). Again, we observe a threshold for dissipation.

Interaction Parameter Magnetic Field (Gauss)

Figure 2: Critical velocities throughout the BEC-BCS crossover. A pronounced maximum was found at resonance. Data is shown for a deep lattice, held for t=500 ms. The solid line is a guide to the eye.

Figure 2 shows the measured critical velocity throughout the BEC-BCS crossover. The maximum near resonance is consistent with the picture of a crossover between two different types of excitation, as discussed in the introduction, and confirms that superfluidity is most robust on resonance.

To illuminate the role of the inhomogeneous density distribution, we performed experiments in which the entire sample was perturbed by a uniform lattice. Lattice beams with m waists probed a more tightly confined sample of pairs, with spatial extent Rm. The onset of dissipation seen in Figure 3 is still striking, but now loss is observed at much lower lattice velocities, in spite of a larger Fermi energy kHz. Moreover, the onset of dissipation is slightly more gradual. When the magnetic field was varied across the Feshbach resonance, we again found a maximum of the critical velocity near resonance. The lowering of the critical velocity due to the inhomogeneous density profile is expected, since at lower density, both the speed of sound and (on the BCS side) the pairing energy decrease. Although the critical velocity should approach zero in the low density wings of the cloud, we still observe a sudden onset of dissipation at a finite velocity, similar to studies in Ref. Rama99 (); Onof00 (); Rama01 (), where a laser beam pierced through the whole condensate, but in contrast to studies reported in ingu05 ().

In the limit of vanishing perturbation, the critical velocity should be given by the Landau criterion. In Figure 4 we address the effects of a finite lattice potential in the original lattice configuration, as depicted in Figure 1. The critical velocity is shown to be a decreasing function of , saturating in the limit of low lattice depth (). This behavior is consistent with numerical simulation iane06 (); Stri07 (). Measurements at the smallest lattice depths had large uncertainties, as the hold time required to observe a heating effect of the lattice approached the natural lifetime of our sample. For this reason we studied the field dependence (Figure 2) at an intermediate lattice depth, where was more well defined.

Lattice Velocity (mm/s)

Figure 3: Effects of density inhomogeneity on the critical velocity. A configuration in which the lattice beams (m) were larger than the trapped sample (m) results in loss in the condensate number at significantly lower velocity. Data is shown for a deep optical lattice held for 200 ms at a magnetic field of 822 G.

For comparison with theory we reference the local Fermi velocity at the trap center mm/s, where is the Fermi velocity of a non-interacting gas at the trap center, and is a universal parameter characterizing unitarity limited interactions Carl03 (); Astr04 (); Carl05 (). For vanishing lattice depth, the observed critical velocity at unitarity approaches . If we use the local Fermi velocity at the waist of the lattice, we obtain . The difference between these values indicates the uncertainty due to residual density inhomogeneity. The local speed of sound in a Fermi gas at unitarity is


The critical velocity for pair breaking is


with Carl03 (); Carl05 () and . These two values should provide approximate upper bounds to the critical velocity at unitarity Sens06 (); Stri06 (). It seems natural that the combination of both excitation mechanisms lowers the critical velocity further. Within these uncertainties, and those of the density, the theoretical predictions agree with the experimental results.

Up until now, we have deferred a discussion of how the moving lattice couples to the excitations. In a pure system at zero temperature, one would expect the excitation spectrum to exhibit discrete resonances, where the perturbation couples only to modes with the -vector of the lattice. On the other hand, at finite temperature, it is possible that the lattice drags along thermal atoms which are point-like perturbations and can create excitations at all -vectors. Our observation that the dissipation sets in at a certain threshold velocity and increases monotonically with velocity is consistent with the participation of the thermal component.

Lattice Depth

Figure 4: Critical velocities at different lattice depths. The results show to be a decreasing function of lattice depth . In the limit of low , converges to a maximum value of 0.25 . Data was taken near resonance, at 822 G () for hold times t = 250 ms, 500 ms, 1000 ms, 2000 ms (squares, diamonds, circles, triangles).

We further elucidated the role of thermal excitations, by varying the temperature. Gradually reducing the trap depth from to , during exposure to a lattice moving above , will suppress the accumulation of a thermal component. The lifetime in this case exceeded that for a sample held at a fixed depth of either or . For Bose-Einstein condensates, theoretical papers emphasized the role of the thermal component in the Landau damping process in a moving lattice griff04 (); niku05 (). This was confirmed qualitatively in an experiment at Florence ingu05 () in which the lifetime of the sample was drastically improved by eliminating the thermal atoms.

In our experiments, the clouds heated up during the exposure to the moving lattice. Figure 5 shows the increase in the number of thermal atoms and the loss in the total number of atoms due to evaporative cooling. In an idealized model, where density is fixed, constant dissipation would result in a linear decrease in the number of atoms due to evaporative cooling. Our data show an accelerated decrease, possibly reflecting increased dissipation due to the increasing fraction of thermal atoms. However, an accurate model should include the change in density (and therefore critical velocity) during the exposure time. Additional impurity atoms (e.g. sodium atoms) could cause dissipation even at zero temperature and would allow more controlled studies of the dissipation mechanism. Unpaired atoms in clouds with population imbalance may not play this role because of phase separation effects shin06 ().

Hold Time (ms)

Figure 5: Number of pairs which remained in the condensate (filled circles) and thermal component (open circles) after being held in a deep optical lattice moving with velocity mm/s for a variable hold time. The thermal component shows a linear increase (dashed line), whereas showed an accelerated loss, and is fit to a quadratic function (solid line).

Another possible dissipation mechanism in a lattice is the creation of two excitations through a dynamical or modulational instability. Such an instability Wu_PhysRevA.64.061603 () occurs already for weakly interacting particles moving through a lattice with momentum , when they collide and scatter into states with momenta , analogous to optical parametric generation camp06 (). This process is energetically possible only above , where is the Bragg momentum which defines the edge of the Brillouin zone. This corresponds to a velocity mm/s for fermion pairs (and twice this value for single atoms). Since the highest critical velocities we observe are slightly below this threshold, and strongly decrease already for relatively small , it is very unlikely, that dynamical instabilities play a role in our experiments. Moreover, such instabilities should be strongly modified by Pauli blocking. For our ratio of local Fermi momentum to the Bragg momentum of 0.9, the first band is nearly full in the center of the cloud. For Bose-Einstein condensates, it has been recently predicted Altm05 () and experimentally shown mun07 () that strong interactions can lower the threshold for the dynamical instability, close to the Mott-insulator transition. The range of 1-D lattice depths explored here () is far from the 1D Mott-insulator regime. We have observed the loss of coherence which typically accompanies the superfluid to Mott insulation transition to occur only beyond .

In conclusion, we have used a novel optical lattice geometry to determine critical velocities in the BEC-BCS crossover without the complications of strong density inhomogeneity. This configuration could be applied to studies in atomic Bose gases which so far have been limited by the inhomogeneous density Rama99 (); Onof00 (); Rama01 (); ingu05 (). In addition, it would be interesting to study dynamical instabilities for fermions and the role of Pauli blocking. The authors would like to thank Aviv Keshet for experimental assistance. This research has been supported by the NSF and the Office of Naval Research.


  1. Present address: Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138
  2. Website:


  1. S. Giorgini, L. P. Pitaevskii, and S. Stringari, preprint condmat/0706.3360 .
  2. R. Sensarma, M. Randeria, and T.-L. Ho, Phys. Rev. Lett. 96, 090403 (2006).
  3. R. Combescot, M. Y. Kagan, and S. Stringari, Phys. Rev. A 74, 042717 (2006).
  4. A. Spuntarelli, P. Pieri, and G. C. Strinati, pre-print /condmat/0705.2658 .
  5. C. A. R. Sá de Melo, M. Randeria, and J. R. Engelbrecht, Phys. Rev. Lett. 71, 3202 (1993).
  6. R. Haussmann, W. Rantner, S. Cerrito, and W. Zwerger, Phys. Rev. A 75, 023610 (2007).
  7. J. Joseph et al., Phys. Rev. Lett. 98, 170401 (2007).
  8. A. Altmeyer et al., Phys. Rev. Lett. 98, 040401 (2007) 98, 040401 (2007).
  9. C. Raman et al., Phys. Rev. Lett. 83, 2502 (1999).
  10. R. Onofrio et al., Phys. Rev. Lett. 85, 2228 (2000).
  11. C. Raman et al., J. Low Temp. Phys. 122, 99 (2001).
  12. L. De Sarlo et al., Phys. Rev. A 72, 013603 (2005).
  13. J. Chin et al., Nature 443, 961 (2006).
  14. J. Mun et al., preprint condmat/0706.3946 .
  15. M. Bartenstein et al., Phys. Rev. Lett. 94, 103201 (2005).
  16. S. Ianeselli, C. Menotti, and A. Smerzi, J. Phys. B 39, S135 (2006).
  17. J. Carlson, S.-Y. Chang, V. R. Pandharipande, and K. E. Schmidt, Phys. Rev. Lett. 91, 050401 (2003).
  18. G. E. Astrakharchik, J. Boronat, J. Casulleras, and S. Giorgini, Phys. Rev. Lett. 93, 200404 (2004).
  19. J. Carlson and S. Reddy, Phys. Rev. Lett. 95, 060401 (2005).
  20. S. Tsuchiya and A. Griffin, Phys. Rev. A 70, 023611 (2004).
  21. S. Konabe and T. Nikuni, J. Phys. B 39, S101 (2005).
  22. Y. Shin et al., Phys. Rev. Lett. 97, 030401 (2006).
  23. B. Wu and Q. Niu, Phys. Rev. A 64, 061603 (2001).
  24. G. K. Campbell et al., Phys. Rev. Lett. 96, 020406 (2006).
  25. E. Altman et al., Phys. Rev. Lett. 95, 020402 (2005).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description