We study the extended Thomas-Fermi (ETF) density functional of the superfluid unitary Fermi gas. This functional includes a gradient term which is essential to describe accurately the surface effects of the system, in particular with a small number of atoms, where the Thomas-Fermi (local density) approximation fails. We find that our ETF functional gives density profiles which are in good agreement with recent Monte Carlo results and also with a more sophisticated superfluid density functional based on Bogoliubov-de Gennes equations. In addition, by using extended hydrodynamics equations of superfluids, we calculate the frequencies of collective surface oscillations of the unitary Fermi gas, showing that quadrupole and octupole modes strongly depend on the number of trapped atoms.
Superconductivity and superfluidity of fermions require, within the BCS theory, matching of the Fermi energies of the two interacting Fermion species. Difference in the number densities of the two species leads either to a normal state, to phase separation, or---potentially---to exotic forms of superfluidity such as Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state, Sarma state, or breached pair state. We consider ultracold Fermi gases with polarization, i.e., spin-density imbalance. We show that, due to the gases being trapped and isolated from the environment in terms of particle exchange, exotic forms of superfluidity appear as a shell around the Bardeen-Cooper-Schrieffer (BCS) superfluid core of the gas and, for large density imbalance, in the core as well. We obtain these results by describing the effect of the trapping potential by using the Bogoliubov-de Gennes equations. For comparison to experiments, we calculate also the condensate fraction, and show that, in the center of the trap, a polarized superfluid leads to a small dip in the central density difference. We compare the results to those given by local density approximation and find qualitatively different behavior.
We investigate the collective oscillations in the BCS-BEC crossover region of a trapped Fermi superfluid at T=0. We explicitly include the effect of a Feshbach resonance, which leads to a tunable pairing interaction adjusted by the threshold energy of the resonance, as well as the associated mol ecules. In a previous paper,we obtained solutions of the Bogoliubov-de Gennes coupled equations, describing a trapped Fermi superfluid ga s. Using these Bogoliubov quasi-particle excitations, we calculate the density correlation function in the Hartree-Fock-Bogoliubov generalized random phase approximation (HFB-GRPA). Previous discussions of collective modes were based on solving the T=0 equation of motion for the superfluid in a trap, ignoring the effect of single-particle excitations. We determine the frequencies of the quadrupole and the monopole modes. In the crossover region, the frequencies of these two modes are shown to be stro ngly depressed by the threshold energy of the two-particle continuum, originating from Andreev single-particle bound states (or in-gap states) localized at the edge of the trap. This effect becomes less important as one enters the BEC regime, due to the absence of these low energy fermion bound states. This suppression of collective mode frequencies by the two-particle continuum may be a useful experimental signature of the single-particle energy gap in the crossover region of trapped Fermi superfluids.
We investigate the single-particle properties at
T=0 of a trapped superfluid gas of Fermi atoms with a Feshbach resonance. A tunable pairing interaction associated with the Feshbach resonance leads to the BCS--Bose-Einstein condensate (BEC) crossover, where the character of superfluidity continuously changes from the BCS-type to a BEC of composite bosons. In this paper, we extend our previous work for a uniform superfluid Fermi gas [Y. Ohashi and A. Griffin, Phys. Rev. A 67, 063612 (2003)] to include the effect of a harmonic trap. We do not use the local density approximation (LDA), but directly solve the Bogoliubov--de Gennes (BdG) coupled equations. While our explicit numerical solutions are for a weak (narrow) Feshbach resonance, we argue that the single-particle BdG excitation spectrum will exhibit the same features for a strong (broad) Feshbach resonance. Using these equations, we find self-consistent values for the spatially dependent local density
n(r) as well as the composite BCS order parameter $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\Delta}}(\mathbf{r})$, the latter describing both the Cooper-pair and molecular condensate contributions. Using these results, we calculate the single-particle density of states in the crossover region, and from this determine the true single-particle energy gap
(Eg) of the trapped Fermi superfluid. This is associated with the in-gap (or Andreev) states in the low-density region at the edge of the trap. We calculate the laser-induced tunneling current $I(\ensuremath{\omega})$ into another hyperfine state, as measured in recent rf spectroscopy experiments. This rf spectrum gives a direct probe of the quasiparticle spectrum. We show how the high-energy part of $I(\ensuremath{\omega})$ gives information about $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\Delta}}(r=0)$ at the center of the trap (which is comparable to the Fermi energy ${\ensuremath{\epsilon}}_{F}$ in the crossover region). We show that $I(\ensuremath{\omega})$ is very dependent on the spatial profile of the pair potential $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\Delta}}(r)$ that is used. We also emphasize that the narrow ``unpaired atom" peak in the rf data gives information about
Eg and the low-energy $(⪡{\ensuremath{\epsilon}}_{F})$ in-gap states of a Fermi superfluid. While our calculations are limited to
T=0, we use them to discuss the recent data of Chin et al. and the LDA calculations of T\"orm\"a and co-workers. The LDA, while useful, can lead to an incorrect physical picture of the low-density surface region of the Fermi superfluid.
Near Feshbach resonances, n | a | 3 ≫ 1, systems of Bose and Fermi particles become strongly interacting/dense. In this unitary limit both bosons and fermions have very different properties than in a dilute gas, e.g., the energy per particle approaches a value ℏ 2 n 2/3 / m times a universal many-body constant. Calculations based upon an approximate Jastrow wavefunction can quantitatively describe recent measurements of trapped Bose and Fermi atoms near Feshbach resonances. The pairing gap between attractive fermions also scales as Δ ∼ ℏ 2 n 2/3 / m near Feshbach resonances and is a large fraction of the Fermi energy—promising for observing BCS superfluidity in traps. Pairing undergoes several transitions depending on interaction strength and the number of particles in the trap and can also be compared with pairing in nuclei.
Near Feshbach resonances, n|a|3 ≫ 1, systems of Bose and Fermi particles become strongly interacting/dense. In this unitary limit both bosons and fermions have very different properties than in a dilute gas, e.g., the energy per particle approaches a value ℏ2n2/3/m times a universal many-body constant. Calculations based upon an approximate Jastrow wavefunction can quantitatively describe recent measurements of trapped Bose and Fermi atoms near Feshbach resonances. The pairing gap between attractive fermions also scales as Δ ∼ ℏ2n2/3/m near Feshbach resonances and is a large fraction of the Fermi energy—promising for observing BCS superfluidity in traps. Pairing undergoes several transitions depending on interaction strength and the number of particles in the trap and can also be compared with pairing in nuclei.