# Tan’s universal contact and collective oscillations of strongly interacting Fermi gases

###### Abstract

We study strongly interacting two component Fermi gas near a Feshbach resonance. By using a ground state energy functional constructed based on asymptotic limits and Monte Carlo calculations, we calculate the contact, structure factor, and collective oscillation frequencies in the BCS-BEC crossover region. The calculated contact and structure factor show excellent agreements with recent experiments. We compare these results with a standard mean-field theory and find that the contact is proportional to the square of superfluid order parameter. Further, we present the chemical potential and the polytropic index in terms of homogenous energy and the contact.

## I I. Introduction

As a result of easy controllability enabled by present day laser technology, experimentalists were able to achieve remarkable progress in trapping and cooling of atomic gases towards degenerate temperatures revCG (). The strongly interacting Fermi gas is one of the richest systems studied by the cold atom community. Although harmonically trapped quantum gases are very dilute systems, most of their properties are governed by the interaction between particles. At these characteristic densities at ultra-cold temperatures, only isotropic and short-range s-wave scattering between particles can take place. This scattering can be characterized by a single parameter, the s-wave scattering length . Experimentally, the scattering length can be tuned by using magnetically tuned Feshbach resonance (FB). In the unitarity limit, where is tuned to positive or negative infinity, the system is strongly interacting. At low temperatures and in the attractive interaction regime where , atoms form Bardeen, Cooper, and Schreifer (BCS) pairs so that the ground state is a BCS superfluid. In the repulsive interaction regime where , atomic potential supports two-body molecular bound state in vacuum. In this limit the ground state corresponds to a Bose-Einstein condensation (BEC) of these molecules. In between these two states, there is a smooth crossover when changes its sign by passing through . At zero temperature, the only relevant dimensionless parameter for physical properties is , where is the Fermi wavevector. Because of the short range nature of the interaction between atoms in a dilute gas, the range of interatomic potential is much smaller than the average inter-particle spacing.

The challenging problem concerning the system properties is how to find the ground state energy . The problem does not have an analytical solution so one has to rely on numerical solutions such as Monte Carlo simulations. In general, the solutions at the limiting cases; such as can be obtained from the asymptotic behavior. At unitarity where , must drop out of the problem so there are only two energy scales available; Fermi energy and the temperature . Here is the mass of the atom and is the atom density in the system. Thus, dimensional analysis follows that the average single particle energy must have the form , where is a dimensionless function. The universal constant has been calculated theoretically tbeta () and is in good agreement with experiments ebeta (). The universal many-body parameter for a unitary gas and for an ideal Fermi gas. The same argument holds for other thermodynamic quantities at unitarity. As the physical properties are independent of microscopic details, by studying a unitary Fermi gas, one learns about the equation of state of a generic strongly interacting gas, such as high density nuclear matter found in the center of neutron stars.

There is, however, another form of ”universality” valid not only at , but also in the entire BCS-BEC crossover region. This universality originated from the fact that , i.e, the range of the interaction is much smaller than the inter-particle distance. Using generalized functions to deal with the singularities associated with the zero-range potential, Shina Tan recently derived a number of universal relations for Fermi gases in the BCS-BEC crossover region tan (). These exact universal relations (Tan relations) connect microscopic properties to thermodynamic quantities. The connection is made through a single quantity termed contact (), which is defined as the high momentum () tail of the momentum distribution ;

(1) |

Remarkably, this contact carries with it all the properties and many-body physics of the system. The Tan relations are applicable for broad situations: homogeneous or trapped systems, many-body or few-body systems, superfluid or normal phases, and finite or zero temperatures. They were later rederived by using a renormalization scheme in the quantum field theoretical framework braaten (), by using a lattice model to regularize the singularity werner (), by using a nonlocal quantum field theory zhang (), and by using a Schrodinger formalism combescot (). For a ground state of a 4-fermion system, the contact has been calculated numerically for different values of interactions and found that the Tan relations are consistent within the systematic errors blume (). The temperature dependence of the universal contact has been investigated by several authors using a high-temperature quantum virial expansion method huc (), ab-initio results palestini (), and a combined mean field-numerical method yu (). Concurrent with the preparation of this manuscript ref. vic () presented complementary BCS mean field results with the argument that the contact is the conjugate of the inverse scattering length. Very recently, the Tan relations were experimentally verified at University of Colorado, Boulder, USA jin () and at Swinburne University of Technology, Melbourne, Australia chris ().

In this paper we investigate the strongly interacting Fermi gases at zero temperature focusing on the universal behavior. By using an energy functional constructed by Kim et al. kim (), first we calculate the contact in the entire BCS-BEC crossover region. Then we calculate the structure factor from the Tan relations. Further, we find that the contact calculated from usual mean field theory is proportional to the square of the superfluid order parameter and shows very good agreement with recent experimental data and the energy functional method. In addition, we generalize the theory to trapped gases using local density approximation and polytropic equation of state. In harmonically trapped systems, we calculate the contact, energies, and collective mode frequencies. We use both hydrodynamic theory and sum rule approach and provide detail discussion about the validity of sum rules for crossover Fermi systems. All the calculated physical quantities are presented in terms of the contact and homogenous energy.

The paper is organized as follows. In the following section, we introduce the energy functional and calculate the contact for a homogenous gas as a function of . Then we present our calculation on structure factor. In section III, we use BCS mean field theory to calculate the contact. In section IV, we generalize our theory to a harmonically trapped Fermi gas. The section V is devoted to discuss the collective oscillation frequencies and derive various oscillation frequencies. Finally, in section V, we provide a discussion and our summary.

## Ii II. The energy functional and the contact

We consider a two component Fermi gas with both spin components equally populated at zero temperature. We neglect the interaction between same spin state, but the interaction between opposite spin states is considered in the form . The total number density and the Fermi wavevector are related through . The ground state energy per particle is determined by the dimensionless interaction parameter . As changes from to to , the ground state smoothly changes from BCS superfluid state to BEC state.

In the BCS limit , the energy per particle can be expanded in powers of :

(2) |

(3) |

where the universal constant is estimated to be tbeta (). In the BEC limit , the energy per particle can be expanded in powers of :

(4) |

where is the boson-boson scattering length bs (). Here the boson is the diatomic molecule formed out of two fermions in the BEC regime. The leading term is the total binding energy density for molecules and the second term is the energy density of molecular BEC. Using the Monte Carlo results at unitarity and asymptotic results at BEC and BCS limits, the energy per particle in the entire BCS-BEC regime has been constructed from a Pade parametrization by Kim et al. kim (); kim2 (). In the BCS regime where ,

(5) |

and in the BEC regime ,

(6) |

The coefficients are found to be , , , and . The molecular energy is . This functional has shown a very good agreement with constrained variational approximation cva () and correctly produces the limiting behavior.

Adiabatic sweep theorem; one of the Tan’s universal relations, connects the contact per unit volume (contact density) and the energy of the system

(7) |

The adiabatic sweep theorem gives connection between the change in total energy density due to adiabatic change in the scattering length and the contact density. As the contact measures the probability of two antiparallel spins being close together, the adiabatic theorem links the short-range behavior to thermodynamics braaten (). By using the energy functional , we calculate the contact density in the BCS-BEC region as a function of . The calculated dimensionless contact density is shown in FIG 1. Finite temperature experimental data from ref. jin () are shown in filled circles. While the black circles represent radio frequency measurements, the gray circles represent the momentum distribution measurements. Notice that the data in figure 2 of ref. jin () shows a different scaling because of their data represents the contact per particle . As can be seen, our zero temperature model gives an excellent agreement with finite temperature experimental data on the BCS side of the resonance. The experiment is done at very low temperature , where is the Fermi temperature. Our calculation indicates that the low temperature contact is not very sensitive to the temperature. The tiny bump at the BEC side of the unitarity in our calculation is an artifact originated from the Pade approximation used to construct the energy functional. Alternatively, one can use the energy functional proposed by Manini et al. manini () to calculate the contact. Even though, Manini et al.’s energy functional is more accurate on the BEC side around , it is discontinuous at . Therefore, the contact shows a large jump at unitarity. In section III below, we compare the contact calculated in the present section with the BCS mean field theory. Later in section IV, we calculate the contact of harmonically trapped gas.

### ii.1 Structure factor

Tan also derived an universal relation for the large momentum behavior of the spin-antiparallel static structure factor . This relation originated from the short-range behavior of the pair correlation function . Taking the limit , and Fourier transforming of gives Tan’s static structure factor

(8) |

This is the large behavior, however, must be smaller than the inverse of range of the interaction potential. In FIG. 2, we present the results for in the BCS-BEC crossover region. The filled circles are the experimental data from Ref. chris () and peng (). The thin line is the theoretical calculation based on a random phase approximation (RPA) from reference hu () and peng (). Except at the BEC side of the unitarity, our model captures the universal power-law () tail of the static structure factor and shows an excellent agreement with experiments chris () and recent theory based on RPA method hu ().

## Iii III. The contact from mean field theory

We restrict ourselves to the wide Feshbach resonance, so that we do not need to explicitly consider the closed channel molecule contribution. The fermions of different hyperfine states and interact through a short-range effective potential . The system of atoms is then described by the Hamiltonian , where and contain kinetic and interaction part of the Hamiltonian defined as,

(9) | |||||

The field operators obey the usual fermionic anticommutation rules, and describe the annihilation of a fermion at position in the hyperfine state . The parameters and are the mass and the chemical potential of the atoms. Using the usual BCS mean-field decoupling bcs (), the zero temperature grand potential is given by

(10) | |||

The BCS-Bogoliubov excitation spectrum of the atoms is given by , where is the kinetic energy and is the superfluid order parameter. Notice that we have already eliminated the ultraviolet divergences originated from the nature of short-range potential. This was done by regularizing the interaction term according to the relation

(11) |

The gap equation is obtained by the minimization of the grand potential density with respect to the superfluid order parameter,

(12) |

The fermions number density is obtained by the variation of grand potential density with respect to the chemical potential,

(13) |

Finally, evaluating from the gap equation, we find the contact density

(14) |

By converting the sum into an integral over the momentum, we numerically solve both gap equation and the number equation simultaneously to find the superfluid order parameter . The mean field contact density is shown in FIG.3. As can be seen, the mean field dimensionless contact shows a reasonable agreement with our energy functional calculation carried out in the previous section. The inset shows the comparison between BCS mean field theory and the experiments. At the BCS side of the unitarity [], the contact calculated from the energy functional method (0.0951) gives a remarkable agreement with the experimental value (0.0951). However, the BCS mean field result (0.1096) is off by at . Obviously, while the contact goes to zero in the weak coupling limit (), it approaches to infinity in the strong coupling limit ().

## Iv IV. Harmonically trapped Fermions

Using the short-range nature of the interaction, Tan derived an additional universal relation for fermions in a harmonic trap;

(15) |

Here , , and are the kinetic, harmonic trapping potential, and interaction energies respectively. This generalized virial theorem was first experimentally verified at unitarity by Thomas et al. thomas () and later by Stewart et al jin () for a range of interaction strengths. In this section, we make use of generalized virial theorem in local density approximation to calculate the contact density in a trapped system.

### iv.1 Chemical potential and polytropic index

In terms of energy functional given above , the chemical potential of the Fermi gas can be calculated by using Gibbs-Duhem relation ma ()

(16) |

Using this relation, we find the chemical potential , where is the dimensionless contact density defined above. The calculated chemical potential as a function of interaction is shown in Fig 4.

The polytropic index relates the pressure and the density of the gas as . Alternatively, relates the chemical potential as . We calculate this effective polytropic index as the logarithmic derivative of the chemical potential . In the BCS-BEC crossover region, we find that the effective polytropic index reads

(17) |

As one expects, we find at both non interacting limit and the unitarity limit heisel (). As we have included the molecular energy in our formalism, we find in the deep BEC limit. Without inclusion of the molecular contribution, we find in the interaction dominated BEC regime. The effective polytropic index in BCS-BEC crossover region is plotted in FIG. 5. In the crossover region, the effective polytropic index which excludes the molecular binding energy varies between .

### iv.2 Contact density in a trapped system

In order to use the generalized virial theorem to find contact density in a harmonic trap, we first calculate the total trapped energy in the trap. Assuming the Fermi gas is locally homogenous, we use local density approximation and the effective polytropic index to calculate the trapped energy. Within local density approximation, the inhomogeneity is taken into account by using spatially varying chemical potential as , where and are the trapping frequency and the chemical potential at the center of the trap. For example, the total number of atoms, can be converted into integration over using . Assuming , we find

(18) |

Similarly, the sum of the kinetic and interaction energies of the trapped system and trapping potential energy can be converted into the integral over . Evaluating the integral, we find . Combining all the energy components, we find the total trapped energy ,

(19) |

The calculated total trapped energy in the BCS-BEC crossover region is shown in FIG. 6. Finally, the ratio of the contact density to the interaction strength in tapped Fermi gas is calculated from the Tan’s generalized virial theorem in Eq. (15),

(20) |

## V V. Collective oscillations

### v.1 Hydrodynamic theory

The dynamic of the gas can be described by the time-dependent nonlinear Schrodinger equation, known as the generalized Gross-Pitaevskii equation pethick (); ma ()

(21) |

where is the superfluid wave function. In terms of superfluid density and velocity , Eq. (21) can be converted into continuity equation and Euler equation,

(22) |

and

(23) |

where is the quantum pressure term. In Thomas-Fermi approximation (TF), we assume that the gas is locally uniform so that the quantum pressure term can be neglected. This TF regime can be realized in a trap in the limit . It is interesting to note that once we drop the quantum pressure term these equations are classical so that these hydrodynamic equations are applicable for both Fermi and Bose superfluids. However, the hydrodynamic equations are sensitive to quantum corrections, statistics, and dimensionality through equation of state which enters through the density dependent local chemical potential . This hydrodynamic description is valid as long as the collisional relaxation time is much smaller than the inverse of the collective oscillation frequencies.

For the ground state and . In order to study the collective oscillations above the ground state, we linearized the hydrodynamic equations by substituting and . After neglecting the higher order terms and keeping only the linear terms and then by combining the continuity and Euler equations, one finds the linearized version of the hydrodynamic theory

(24) |

Taking and , this equation can be converted into a linear eigenvalue problem. In the TF regime, the continuity and Euler equations admit the ground state density , where the Thomas-Fermi radius . Re-scaling the length, density, and oscillation frequency by introducing , , and respectively, Eq. (24) can be written in dimensionless form,

(25) |

For a spherically symmetric trap, the angular momentum and its projection in the z-axis m are good quantum numbers. Following the Ref. heisel () and Ref. yin (), the solution for the density fluctuation have the form , where . The function satisfies the hypergeometric differential equation,

(26) | |||

The solutions of this standard equation is hypergeometric functions heisel (); yin (); ma () and the eigenvalues are given by

(27) |

with radial quantum number and angular momentum quantum number . Using the polytropic index calculated in the previous section, all the modes of collective oscillation frequencies can be presented in terms of the contact density from this energy eigenvalues. For example, monopole mode frequency or the breathing mode frequency is given when and . Zero temperature monopole mode frequency of the two component Fermi gas in a spherically symmetric trap within the BCS-BEC crossover region is shown in FIG. 8. The frequencies are given as a function of the dimensionless interaction parameter . Our model predicts at the weak coupling BCS limit. At the weakly repulsive BEC limit, . These results are consistent with the sum rule approach present in the next subsection. At unitarity, as required by universality. The other mode frequencies, such as second breathing mode (, ) and quadrupole modes () can be readily calculated from Eq. (27).

### v.2 Sum rule approach

In this section, we use the sum rules approach in the linear response theory to express rigorous upper bounds to the energies of the collective oscillations at zero temperature. The response function at complex frequency is written in terms of the imaginary part of the response function as

(28) |

Expanding the response function in terms of , the moment expansion of is given by

(29) |

where the pth moment is defined as

(30) |

In general, one has to calculate the response function associated with the mode excitation operator to evaluate the moments. However, these moments can be computed without calculating the response function if one uses the sum rule approach sumrule (). A major advantage of sum rule approach is that the collective mode frequencies can be calculated without the full solution of the eigenstates of the Hamiltonian. One disadvantage is that this method provides the rigorous upper bounds for the lowest collective mode frequency excited by the excitation operator. The excitation energies of the collective modes are given by the ratio or . When the response function is almost exhausted by a single mode (that is, the state is highly collective), then these ratio give the exact collective mode frequency. We estimate the excitation energies of the collective modes from the ratio lc ()

(31) |

Using the completeness relations, the moments can be written as

(32) |

where represents the commutator between the operators A and B, is the mode excitation operator, and is the Hamiltonian of the system

(33) | |||||

where . The excitation operators corresponding to the dipole, monopole and quadrupole modes are , , and respectively. By evaluating the commutators above, the collective mode frequencies for the dipole, monopole, and quadrupole modes are lc ()

(34) | |||||

Notice that the dipole mode is not effected by the interaction. This mode corresponds to the oscillation of the center of mass of the atomic cloud. The calculation of monopole collective mode frequency from the first and third sum rules involves the evaluation of interaction energy of the trapped system. The calculation of this interaction energy for a crossover Fermi gas, especially at unitarity is not trivial. Indeed, the use of third sum rule for crossover fermions is problematic. For example, consider the equation of motion of the operator . Using Heisenberg equation of motion, one finds and . For a spherically harmonic oscillator potential, the three terms represent the average kinetic energy , the harmonic oscillator potential energy , and the average mean-field interaction energy . In equilibrium is time independent, so . Then summing over all three components, this follows the original virial theorem,

(35) |

Now, let’s combine this original virial theorem in Eq. (35) with the Tan’s generalized virial theorem given in Eq. (15). For both non-interacting fermions and unitary fermions, the interaction energy turns out to be zero at unitarity. In other words, within our sum rule approach, the third sum rule (together with original virial theorem) leads to an incorrect interaction energy at unitarity. Regardless to this incorrect interaction energy, the combination of these two virial theorems [Eqs. (15) and (35)] into the sum rules produces the correct breathing mode frequency. In any case, we avoid using Eq. (35), instead we use compressibility sum rule to evaluate the breathing mode frequency. Contrast to the other sum rules, the compressibility sum rule

(36) |

cannot be written as a set of commutators. By completing the Fourier transform over the frequency (or using the Kramers-Kronig relation), one finds that the compressibility sum rule is the static response function, where . Using for the breathing mode frequency, we find

(37) |

Defining , where and using , we find

(38) |

The lowest order breathing mode frequency for spherically symmetric trap is then given by . This is exactly the same results predicted by the hydrodynamic theory in previous subsection. This indicates that the rigorous upper bounds for crossover Fermi gases calculated from the sum rules are very tight. In other words, the lowest order breathing mode frequency of a trapped Fermi gas is highly collective.

### v.3 Collective modes of a Fermi gas in an axially symmetric trap

In an axially symmetric trap with trapping frequency , the monopole modes frequencies can be calculated easily using the polytropic index . Generalizing the hydrodynamic theory for an axially symmetric trap, the lowest frequency axial monopole mode and the radial monopole mode and respectively heisel (); ma (). These expressions are valid only for a cloud with aspect ratio . The other limit where , the monopole modes are and heisel (); ma (). The experimental results are available only for the case cmex1 (); cmex2 (). These two mode frequencies in the BCS-BEC crossover region are shown in FIG. 9. For the limiting cases when and , the axial mode frequencies are and and the radial mode frequencies are and lc ().

## Vi VI. Summary

By using two methods; a ground state energy functional constructed based on asymptotic limits and Monte Carlo calculations and the BCS mean field theory, we calculate the contact of a two component Fermi gas near a Feshbach resonance. The calculated contact which encapsulates all the many body physics in the BCS-BEC crossover region shows excellent agreement with the recent experiments. Within the mean field theory, we find that the contact is proportional to the square of the superfluid order parameter. We then investigate the large momentum behavior of the structure factor using a Tan’s universal relation and the contact. Our structure factor calculation in the BCS-BEC crossover region also shows an excellent agreement with recent experiments and a RPA based theoretical calculation.

We combine the Tan’s generalized virial theorem and local density approximation to derive the contact density in harmonically trapped fermions. Using the polytropic equation of state, we calculated the trapped energies and the collective mode frequencies of a Fermi gas trapped in a harmonic oscillator potential. We present all these quantities in terms of the homogenous energy density and the homogenous contact density.

## Vii Acknowledgments

We are grateful to John Gaebler for sending us their experimental data for the contact. We thank Hui Hu for sharing with us their experimental data and RPA based theoretical results. We acknowledge Andrew Snyder for very enlightening discussions and critical comments on the manuscript.

## References

- (1) For a review on ultra-cold gasses see for example, S. Giorgini, L. P. Pitaevskii, and S. Stringari, Rev. Mod. Phys. 80, 1215 (2008) and I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 80, 885 (2008).
- (2) J. Carlson, S.-Y. Chang, V. R. Pandharipande, and K. E. Schmidt, Phys. Rev. Lett. 91, 050401 (2003).
- (3) K. M. OHara, S. L. Hemmer, M. E. Gehm, S. R. Granade, and J. E. Thomas, Science 298, 2179 (2002).
- (4) Shina Tan, Ann, Phys. 323, 2952 (2008); Shina Tan, Ann, Phys. 323, 2971 (2008); Shina Tan, Ann, Phys. 323, 2987 (2008).
- (5) E. Braaten and L. Platter, Phys. Rev. Lett. 100, 205301 (2008); E. Braaten and L. Platter, Phys. Rev. A 78, 053606 (2008).
- (6) F. Werner, Phys. Rev. A 78, 025601 (2008); F. Werner, L. Tarruell, Y. Castin, Eur. Phys. J. B 68, 401 (2009); F. Werner, Y. Castin, http://arxiv.org/abs/1001.0774.
- (7) S. Zhang and A. J. Leggett, Phys. Rev. A 79, 023601 (2009).
- (8) R. Combescot, F. Alzetto, and X. Leyronas, Phys. Rev. A 79, 053640 (2009).
- (9) D. Blume and K. M. Daily, Phys. Rev. A 80, 053626 (2009).
- (10) H. Hu, X. Liu, P. D. Drummond, New J. Phys. 13 035007 (2011).
- (11) F. Palestini, A. Perali, P. Pieri, and G. C. Strinati, Phys. Rev. A 82, 021605(R) (2010).
- (12) Z. Yu, G. M. Bruun, and G. Baym, Phys. Rev. A 80, 023615 (2009).
- (13) V. Romero-Rochin, http://arxiv.org/abs/1012.0236.
- (14) J. T. Stewart, J. P. Gaebler, T. E. Drake, and D. S. Jin, Phys. Rev. Lett. 104, 235301 (2010).
- (15) E. D. Kuhnle, H. Hu, X.-J. Liu, P. Dyke, M. Mark, P. D. Drummond, P. Hannaford, C. J. Vale, Phys. Rev. Lett. 105, 070402 (2010).
- (16) Y. E. Kim and A. L. Zubarev, Phys. Rev. A 70, 033612 (2004).
- (17) D. S. Petrov, C. Salomon, and G. V. Shlyapnikov, Phys. Rev. Lett. 93, 090404 (2004).
- (18) Y. E. Kim and A. L. Zubarev, Phys. Lett. A 327, 397 (2004).
- (19) H. Heiselberg, Phys. Rev. Lett. 93, 040402 (2004).
- (20) Peng Zou, Eva D. Kuhnle, Chris J. Vale, Hui Hu, Phys. Rev. A 82, 061605(R) (2010).
- (21) H. Hu, X.-J. Liu and P. D. Drummond, Euro physics Letters 91, 20005 (2010).
- (22) A. J. Leggett, in Modern Trends in the Theory of Condensed Matter, edited by A. Pekalski and R. Przystawa (Springer-Verlag, Berlin, 1980); D. M. Eagles, Phys. Rev. 186, 456 (1969).
- (23) For example see T. N. De Silva and E. J. Mueller, Phys. Rev. A 72, 063614 (2005); O. Bohigas, A. M. Lane, and J. Martorell, Phys. Rep. 51, 267 (1979); E. Lipparini and S. Stringari, Phys. Rep. 175, 103 (1989).
- (24) J. E. Thomas, J. Kinast, and A. Turlapov, Phys. Rev. Lett. 95, 120402 (2005).
- (25) C. J. Pethick and H. Smith, BoseEinstein Condensation in Dilute Gases, Cambridge University Press, Cambridge, 2001
- (26) H. Heiselberg, Phys. Rev. A 63, 043606 (2001).
- (27) J. Yin and Y. Ma, Phys. Rev. A 74, 013609 (2006).
- (28) N. Manini and L. Salasnich, Phys. Rev. A 71, 033625 (2005).
- (29) J. Yin and Y. Ma and G. Huang, Phys. Rev. A 74, 013609 (2006).
- (30) J. Kinast, S. L. Hemmer, M. E. Gehm, A. Turlapov, and J. E. Thomas, Phys. Rev. Lett. 92, 150402 (2004).
- (31) M. Bartenstein, A. Altmeyer, S. Riedl, S. Jochim, C. Chin, J. H. Denschlag, and R. Grimm, Phys. Rev. Lett. 92, 203201 (2004).
- (32) L. Vichi and S. Stringari, Phys. Rev. A 60, 4734 (1999); S. Stringari, Europhys. Lett. 65, 749, (2004).