Introduction

1. Model Hamiltonian A strongly interacting Fermi gas near a wide Feshbach resonance can be well described by the following single-channel Hamiltonian: where U is the atom-atom interaction rate, εkσ= k2/(2mσ) (σ=↑,↓) with mσdenotes the mass of species-σ, and μσ is the chemical potential associated with its population. We define μ=(μ↑+μ↓)/2 and h=(μ↑-μ↓)/2 for convenience in later discussion. Under the mean-field approximation, one can obtain the thermodynamical potential: where Here we use the local density approximation to account for the trapping potential. One can solve for Δ, the superfluid order parameter (T=0) or the total gap (finite T), by minimizing Ω. The number constraints are where uk2=[Ek+(εkr-μr)]/(2Ek), vk2=[Ek-(εkr-μr)]/(2Ek), and f(x)=1/(1+ex/T), subject to N=(n↑+n↓)d3r and p=(n↑-n↓)d3r/N. mr is the reduced mass mass mismatch parameter quasi-particle dispersion relation local chemical potential Superfluid shells for trapped fermions with mass and population imbalanceG.-D. Lin, W. Yi*, and L.-M. Duan FOCUS center and MCTP, Department of Physics, University of Michigan, Ann Arbor, Michigan 4810930We map out the phase diagram of strongly interacting fermions in a potential trap with mass and population imbalance between the two components. As a unique feature distinctively different from the equal-mass case, we show that the superfluid here forms a shell structure which is not simply connected in space. Different types of normal states occupy the trap regions inside and outside this superfluid shell. We calculate the atomic density profiles, which provide an experimental signature for the superfluid shell structure. Zero-Temperature Phase Diagrams Atomic Number Distributions in a Trap Introduction The recent experimental advance in superfluidity in polarized ultracold Fermi gases has raised strong interest in studying the phase configuration of this system. The experiments have suggested a phase separation picture with a superfluid core surrounded by a shell of normal gas. This picture has been confirmed also by a number of theoretical calculations of the atomic density profiles in the trap. As Feshbach resonances between different atomic species with unequal masses have been reported, the next step is to consider the properties of a strongly interacting Fermi gas with both mass and population imbalance between the two components. There have been a few recent theoretical works in this direction, with focus on the properties of a homogeneous gas. The solid (dashed) curves are for 40K (6Li) atoms, respectively. The insets in (a) and (b) show the amplified tails of the profiles. The other parameters are (a) T=0, resonance, and p=0.2, (b) T=0.1TF, resonance, and p=0.2, (c) T=0, kFas−1=−1 (BCS side), p=−0.4, and (d) T=0, kFas−1=0.2 (BEC side), p=0.3. Fermi-Surface Mismatch and Shell Structure One can understand the shell structure introduced by unequal masses by first considering the two Fermi surfaces in k-space formed by the NON-interacting atoms of two species. The radius of the Fermi surface at a certain location then can be found to be kFσ(r)=[2mσ(μσ-r2)]1/2 (in a standard unit), which decreases faster for massive atoms as r increases. When the interaction is turned on, pairing is preferred only if the two kF-surfaces are sufficiently close. Otherwise, the local phase remains normal. This picture should be qualitatively valid even if the atoms strongly interact. SF: Superfluid state (balanced population) NMx: Normal mixture with x as majority NPx: Normal fully polarized gas of species x BP1: Breached pairing phase of type 1 (only one zero in the dispersion relation) Double-to-single Well Transition Thermo-potential Ωas a function of Δ. Upper curves correspond to more eccentric locations away from the center of the trap. [1] W. Yi and L.-M. Duan, Phys. Rev. A 73, 051602(R) (2006); 74, 013610(2006). [2] M. Iskin and C.A.R. Sá de Melo, Phys. Rev. Lett. 97, 100404 (2006); e-print cond-mat/0606624. [3] S.-T Wu, C.-H. Pao, S.-K. Yip, Phys. Rev. B 74, 224504 (2006). _______________ * Current Address: Institut für Theoretische Physik, Universität Innsbruck, Technikerstrasse 25, A-6020, Innsbruck, Austria (6Li↓-40K↑ , α=0.74 p=0.2 , T=0 , On resonance )