the bcs to bec crossover problem of a fermi gas with an adjustable attractive interaction has been investigated theoretically for quite some time @xcite . for low temperatures
the gas is superfluid and , in the case of s - wave interactions , it exhibits a smooth crossover from the well known bcs regime of weakly bound cooper pairs to the bec regime of tightly bound bosonic dimers with a residual repulsive interaction @xcite . in recent years , this crossover has been realized experimentally using ultracold fermi gases in optical traps , where the interaction can be tuned using feshbach resonances @xcite . in the experimentally relevant case of
so called broad feshbach resonances , which is in principle always realized in the dilute limit , the physical properties of the homogeneous gas at equal densities for both spin components are described by only two dimensionless parameters : the interaction strength @xmath3 and the temperature @xmath4 .
here , @xmath5 is the @xmath6-wave scattering length which fully characterizes interactions in the dilute , ultracold limit , while the scales for length and energy are determined by the fermi wave number @xmath7 and the fermi energy @xmath8 , respectively , where @xmath9 is the particle density .
a particularly interesting regime is located near the feshbach resonance , where the scattering length @xmath5 is infinite . at this point and , more generally , in the so - called _ unitary regime _ where @xmath10 , the dimensionless interaction parameter @xmath11 disappears from the problem .
all thermodynamic quantities are therefore universal functions of the dimensionless temperature @xmath4 @xcite . on a microscopic level ,
the unitary gas exhibits a particular kind of scale invariance , similar to a gas with purely inverse square two - particle interactions @xcite .
more generally , as shown by nikolic and sachdev @xcite , universality is not restricted to the unitary regime .
it is tied to the fact that the unitary balanced gas at zero density is an unstable fixpoint with only three relevant perturbations . since there is no small expansion parameter
, the unitary regime is the most challenging one from a theoretical point of view .
in addition , it is in fact precisely this regime which is accessible experimentally ( see e.g.  the recent review articles by ketterle and zwierlein @xcite , by bloch _
et al . _
@xcite and by giorgini _
et al . _
@xcite ) .    in a recent paper @xcite ,
we have presented a field theoretic approach for the thermodynamics of the bcs to bec crossover , which is based on the formalism developed by luttinger - ward @xcite and dedominicis - martin @xcite . in the following
, this approach for the homogeneous gas is applied to calculate the thermodynamic properties of the trapped fermi gas , using a local density approximation .
we compare our results with a recent experiment by luo _
et al . _
@xcite and with recent theories @xcite .
in particular , we provide results for the entropy as a function of temperature , which allows to do reliable thermometry for the trapped unitary gas and also gives a precise value for the critical temperature and the associated entropy per particle .
in addition , we show that starting well in the superfluid regime , much lower temperatures and entropies can be reached by removing atoms from the edge of the cloud , which carry most of the entropy .
in our previous paper we have calculated the thermodynamic quantities for the homogeneous system . at a given particle density @xmath9 , these are the internal energy per particle @xmath12 and the entropy per particle @xmath13 . the fermi wave number @xmath7 and the fermi energy @xmath8 can be used as scale factors in order to make the thermodynamic quantities dimensionless .    in the presence of an external confining potential @xmath14 ,
the particle density @xmath15 is non - uniform . within a local density approximation ( lda ) ,
thermodynamic quantities like pressure @xmath16 , chemical potential @xmath17 , entropy per particle @xmath18 or the internal energy per particle @xmath19 are then also spatially varying , being determined by the corresponding equilibrium values in the uniform system evaluated at the local density @xmath15 .
the local density approximation neglects the dependence of thermodynamic properties on density gradients . at zero temperature
, it is essentially a zeroth order semiclassical approximation @xcite .
it is valid as long as the local fermi wave number @xmath20 times the oscillator length @xmath21 defined by the characteristic frequency @xmath22 of the confining potential is much larger than one . except near the edge of the cloud , where the density approaches zero ,
this condition is well justified for most of the experiments because typical fermi energies @xmath23 are of the order of several khz , while the trapping frequencies are around @xmath24 or smaller .
specifically , near the trap center @xmath25 , where @xmath26 is the total particle number in a trap .
as will be shown below , the finite size corrections to the ground state energy are of relative order @xmath27 in a harmonic trap .
they are therefore negligible , at least for global observables , for the typical particle numbers in experiment , where @xmath28 @xcite . this conclusion is also supported by a recent comparison of lda with a numerical solution of the bogoliubov - degennes equations @xcite .    within the lda ,
the global thermodynamic quantities of the trapped fermi gas are obtained by integrating over the whole space . specifically , we define @xmath29 \ , \qquad   \label{eq : total_energy } \\ s & = & \int d^3r\ n(\mathbf{r})\ s(\mathbf{r } ) \ , \label{eq : entropy}\end{aligned}\ ] ] as the particle number , the potential energy , the internal energy , the total energy , and the total entropy . here
the particle density acts like a distribution function to define averages over the trap .    in an optical trap where the laser intensity profiles are gaussian functions ,
the confining potential @xmath14 is given by an anisotropic gaussian function .
close to the center of the trap this potential can be approximated by an anisotropic harmonic function @xmath30 where @xmath31 is the mass of the atoms and @xmath32 , @xmath33 , @xmath34 are the harmonic oscillator frequencies in the three spatial directions . in the following ,
we use the harmonic potential and neglect the anharmonic terms . since in lda @xmath17 is the chemical potential relative to the potential @xmath14 , in thermal equilibrium the total chemical potential @xmath35 is constant which implies the condition @xmath36 this equation together with the requirement of a constant temperature @xmath37 determines the spatial dependence of all local thermodynamic quantities . the particle density
@xmath15 implies a local fermi wave number @xmath20 and a local fermi energy @xmath38 . as a consequence the dimensionless parameters @xmath39 and
@xmath40 depend on the local position in the trap .
it is convenient to define the weighted radial coordinate @xmath41 by @xmath42 where @xmath43 is the average harmonic frequency . with this definition
the confining potential acquires the simple form @xmath44 . as a consequence ,
the anisotropy of the trap becomes irrelevant because all local quantities @xmath45 etc .
only depend on the weighted radial coordinate @xmath41 . rewriting the space integrals in eqs .  - in terms of @xmath41 , this applies also to all thermodynamic quantities derived from the local density . a convenient measure for the overall length scale in a trap
is provided by the thomas - fermi radius @xmath46 of the confined _ non - interacting _ fermi gas at zero temperature with a given total particle number @xmath26 .
similarly , as a characteristic scale for the energy we define the corresponding fermi energy @xmath47 of the non - interacting gas .
the thermodynamic quantities can be made dimensionless by considering the ratios @xmath48 , @xmath49 , and @xmath50 .
these ratios depend on the space coordinate @xmath51 only implicitly via the dimensionless parameters @xmath39 and @xmath40 .
a particular case is the unitary gas at zero temperature , where both parameters vanish identically @xmath52 , @xmath53 . using standard thermodynamic relations , it is straightforward to show that all thermodynamic quantities can be expressed in terms of a _ single _ dimensionless parameter , the so - called bertsch parameter @xmath54 @xcite in the form @xmath55 these relations hold both for an ideal fermi gas , where @xmath56 and also at the unitarity point , where @xmath54 has a nontrivial value @xcite . in our previous work @xcite
, we have calculated both ratios independently and obtained @xmath57 and @xmath58 .
the first ratio implies a bertsch parameter @xmath59 , while the second ratio implies @xmath60 .
these two values differ by about @xmath61 .
consequently , the relation @xmath62 , which is valid both for an ideal and a unitary fermi gas , is satisfied up to an error of @xmath61 .    inserting the harmonic potential into eq .   and using the weighted radial coordinate
, we obtain the chemical potential @xmath63 \label{eq : chemical_potential_2}\ ] ] where @xmath64 is the thomas fermi radius of the unitary gas at zero temperature .
the dimensionless ratios imply similar functional forms for the other quantities , e.g.   @xmath65^{1/2}$ ] for the local fermi wavevector .
these expressions are valid only for @xmath66 , because the particle density @xmath15 is nonzero only in this case and zero otherwise . evidently , the thomas - fermi radius @xmath64 of the unitary fermi gas is the only parameter which determines the spatial dependence of the thermodynamic quantities .
it is related to the thomas - fermi radius @xmath67 of the ideal fermi gas by the bertsch parameter via @xmath68^{1/4 } = \xi^{1/4 } \ .
\label{eq : thomas_fermi_radius}\ ] ] using the dimensionless ratio @xmath57 of our numerical calculations @xcite we obtain the result @xmath69 .
this value agrees very well with the most recent field theoretic result @xmath70 for the bertsch parameter , obtained from a borel resummation of an expansion around the upper critical dimension four , carried out to three loop order @xcite .
it is somewhat smaller , however , than the result @xmath71 found by using the bertsch parameter @xmath72 that follows from variational monte carlo calculations @xcite or from a theory that includes the gaussian fluctuations around the bcs mean field , extended to arbitrary coupling , where @xmath73 @xcite .
next we insert the local internal energy per particle @xmath19 , the potential @xmath74 , and the local particle density @xmath15 into eqs .  - in order to calculate the energies of the unitary fermi gas in the harmonic trap .
we thus obtain the dimensionless ratios @xmath75^{1/2 }   \nonumber\\   & = & 0.449 \ , \label{eq : potential_energy_3 } \\
2\ , \frac{u}{n e_f } & = & \frac{3}{4}\ , \frac{u(\mathbf{0})}{e_f } = \frac{3}{4}\ , \frac{u(\mathbf{0})}{\varepsilon_f(\mathbf{0 } ) } \left [ \frac{\mu(\mathbf{0})}{\varepsilon_f(\mathbf{0 } ) } \right]^{-1/2 } \nonumber\\ & = & 0.440 \
, \label{eq : internal_energy_3 } \\
\frac{e}{n e_f } & = & \frac{u+e_\mathrm{pot}}{n e_f } = 0.444 \ .
\label{eq : total_energy_3}\end{aligned}\ ] ] the explicit numbers are obtained by inserting the ratios @xmath57 and @xmath58 of our numerical calculation @xcite .    for a harmonic potential @xmath14
it is well known that the internal energy @xmath76 and the potential energy @xmath77 are related to each other by the virial theorem @xmath78 . as shown by thomas _
et al . _
@xcite , this theorem also holds for the interacting fermi gas in the unitary regime . as a consequence , the results of eqs .
- should be equal . by using the dimensionless ratios of the homogeneous system we can express the results in terms of the bertsch parameter @xmath54 according to @xmath79 in practice , our resulting numbers differ by about @xmath61 .
this difference is related to the fact , that the exact relation @xmath62 for the unitary gas is satisfied only with an error of @xmath61 in our theory .
the result for the ground state energy in the trap is based on using lda and provides the exact leading order contribution in the limit @xmath80 .
the question of how large the subleading corrections to this result are has been addressed by son and wingate @xcite . using a gradient expansion of the effective field theory describing the low energy physics of the fermionic superfluid at unitarity , they have determined the @xmath81-corrections to the density response , which is equal to the uniform compressibility @xmath82 at @xmath83 .
these corrections give rise to an additional contribution @xcite @xmath84   \ .
\label{eq : energies_5}\ ] ] to the ground state energy , which contains the two dimensionless coefficients @xmath85 and @xmath86 which appear beyond the leading coefficient @xmath54 of the uniform system in an expansion up to second order in gradients .
the leading correction to lda is thus of relative order @xmath27 ( formally it is @xmath87 ) and describes a _ curvature _ instead of the naively expected _ surface _
correction , which would scale like @xmath88 .
the absence of a surface correction also apppears for an ideal fermi gas and is a peculiarity of the harmonic confinement @xcite .
the prefactor of the @xmath27 correction in has been determined from an expansion around dimension four by rupak and schfer @xcite and is @xmath89 for an isotropic trap . for experimentally relevant particle numbers @xmath28 , the beyond lda relative corrections to the ground - state energy
are therefore only around @xmath90 and thus are clearly below the present experimental accuracy .
the final results of this section are the thomas - fermi radius @xmath69 and the total energy @xmath91 of the unitary fermi gas at zero temperature in a harmonic trap , which depend only on the bertsch parameter @xmath59 of the homogeneous system by eqs .   and .
these results will be compared with experiments and other theories in the next section , where we discuss the situation at finite temperatures .
for nonzero temperatures , the local thermodynamic quantities are space dependent via the scale factors @xmath20 , @xmath38 , and also via the dimensionless temperature @xmath40 . in the unitary regime
the dimensionless interaction @xmath92 is constant . since the effective temperature increases towards the edge of the cloud
, the local particle density @xmath15 does not follow the simple @xmath93^{3/2}$ ] law valid at @xmath2 and eq .
must be solved numerically in order to obtain the detailed density profile @xmath45 as a function of the weighted radial coordinate @xmath41 .     of the unitary fermi gas in a harmonic trap as a function of the radius @xmath41 for several temperatures @xmath94 ( black solid ) , @xmath95 ( brown short dashed ) , @xmath96 (
red long dashed ) , @xmath97 ( orange short dashed ) , @xmath98 ( green long dashed ) , @xmath99 ( blue dot - dashed ) , @xmath100 ( magenta doubledot - dashed ) , and @xmath101 ( turquoise dot - doubledashed ) .
the length scale is @xmath102.,scaledwidth=44.0% ]    using the results of our previous numerical calculations @xcite for the homogeneous system , we obtain the density profiles @xmath103 for several values of the temperature @xmath37 which are shown in fig .
[ fig : density_profiles ] .
the black solid line is the density profile for zero temperature .
it is nearly identical to the expected zero temperature thomas - fermi profile .
differences are due to the limited numerical accuracy , which are , however , much smaller than the thickness of the line . at zero temperature , the fermi gas is superfluid in the whole trap .
the colored and dot - dashed lines represent the density profile for nonzero temperatures .
for very low temperatures , changes occur only close to the surface of the cloud where @xmath104 . for this reason ,
the brown short dashed line is visible only for @xmath105 while otherwise it is nearly on top of the black solid line .
the red long dashed line is still very close to the black solid line . at a position in space
@xmath51 the fermi gas may be normal fluid or superfluid if the local dimensionless temperature @xmath106 is above or below the critical value @xmath107 . for the brown , red , orange , and green dashed lines there exists a respective weighted radius @xmath108 so that @xmath109 . in these cases
the fermi gas is superfluid in the inner region of the trap where @xmath110 and normal fluid in the outer region where @xmath111 . for the blue , magenta , and turquoise dot - dashed lines the dimensionless temperature is @xmath112 for all positions in space @xmath51 so that the fermi gas is normal fluid in the whole trap .
since @xmath113 has its lowest value for @xmath114 , the dimensionless temperature at the center of the trap @xmath115 determines the superfluid transition of the confined fermi gas .
for @xmath116 the fermi gas is completely normal fluid , while for @xmath117 there exists a superfluid region close to the center .    in a homogeneous gas ,
the normal to superfluid transition is a continuous phase transition of the 3d xy type along the complete bcs to bec crossover , because the broken symmetry is that associated with a complex scalar order parameter .
by contrast , our approach @xcite predicts a weak first - order superfluid transition because the superfluid phase of the luttinger - ward theory does not smoothly connect with the normal - fluid phase at a single critical temperature @xmath118 . as a result , there are two slightly different critical temperatures @xmath119 and @xmath120 . the upper value @xmath119 is defined by the condition that for temperatures above this value the superfluid order parameter vanishes , while the lower value @xmath120 is defined by the temperature below which the normal - fluid phase is no longer stable .
fortunately , the difference between both temperatures , which should vanish in an exact theory , is rather small over essentially the whole bcs to bec crossover .
in particular , at unitarity , the upper and lower values for @xmath121 are @xmath122 and @xmath123 , which is within the present numerical uncertainties in the determination of the critical temperature of the unitary gas .
indeed , our critical temperature agrees very well with the most precise calculations of @xmath121 so far by quantum monte carlo calculations , which give @xmath124 for the uniform gas at unitarity @xcite , a value that has been confirmed very recently @xcite .
the existence of two different critical temperatures leads to a multivaluedness of thermodynamic quantities , which is an artifact of the first order nature of the superfluid transition within our theory . in order to avoid multivalued local density or entropy profiles ,
we have connected the normal and superfluid branches with a kink at the point , where they are closest , thus providing an optimal approximation to the exact continuous profiles in a theory which properly accounts for the continuous nature of the transition in the infinite system .
this point is related to the upper value of the dimensionless critical temperature @xmath119 .
hence , in fig .  [ fig : density_profiles ] the brown , red , orange , and green dashed lines have kinks located at @xmath125 , @xmath126 , @xmath127 , @xmath128 , respectively .
the inner , superfluid branches of these lines show a bulge in the center of the trap , which is appreciable , in particular for the green curve . the formation of a small bulge upon condensation , which first appears near the trap center , has indeed been observed experimentally @xcite , although the effect in axially integrated profiles is rather small .
previous theoretical results for the density profile have been obtained by bulgac _
et al . _
@xcite within a monte - carlo approach . for zero temperature ,
their result is very close to our black solid line in fig .
[ fig : density_profiles ] .
they obtain the thomas - fermi radius @xmath129 @xcite . this is larger than our value @xmath130 because of the larger value @xmath131 of the bertsch parameter . for nonzero temperatures ( the same as those in fig .
[ fig : density_profiles ] ) , bulgac _ et al .
_  obtain single valued density profiles which agree qualitatively with our results but differ somewhat quantitatively .
they also observe a superfluid bulge which , however , is smaller .     of the unitary fermi gas in a harmonic trap as a function of the radius @xmath41 for several temperatures .
the curves are related to those in fig .
[ fig : density_profiles ] .
an additional curve is included for @xmath132 ( black short dashed).,scaledwidth=44.0% ]    a very interesting quantity is the local entropy per particle @xmath133 . from our numerical calculations
@xcite of the homogeneous system , and within lda , the resulting profiles @xmath134 of the local entropy per _ volume _ are shown in fig .
[ fig : entropy_profiles ] for several values of the temperature @xmath37 .
the colors of the curves are related to those in fig .
[ fig : density_profiles ] . for high temperatures where the fermi gas is completely normal fluid ( blue , magenta , and turquoise dot - dashed curves )
, the entropy density is distributed over the whole trap with a maximum at the center . for low temperatures ( orange , red , and brown dashed curves )
the center part is superfluid but the outer part is normal fluid .
consequently , in these cases the entropy density is minimum in the center of the trap but maximum at a nonzero radius .
for very low temperatures ( black dashed line ) , the main contribution of the entropy is located close to the surface of the atom cloud .
the width of this surface layer decreases for decreasing temperature and eventually shrinks to zero in the zero temperature limit ( black solid line ) .    in fig .
[ fig : entropy_profiles ] we have again eliminated multivalued regions which are an artifact of the first - order nature of the transition in our theory for the homogeneous system by using the same criterion as in the density profiles of fig .  [
fig : density_profiles ] . indeed , since the normal to superfluid transition is not associated with a latent heat , the local entropy density will be single valued . an interesting observation that is evident from fig .
[ fig : entropy_profiles ] , is that sufficiently below the onset of superfluidity in a trapped fermi gas , most of the entropy is located in the normal region near the cloud edge .
this observation suggests an efficient way to lower the entropy further by removing atoms in the boundary layer and simultaneously readjusting the trapping potential such that the now smaller system has a radius close to that beyond which atoms have been removed .
this idea is - of course - in the same spirit than the standard evaporation cooling in the normal state @xcite , however it is much more efficient .
indeed , consider starting with a dimensionless temperature @xmath135 ( red long dashed curve ) or - equivalently - an entropy per particle @xmath136 which is close to that reached in current experiments .
removing about @xmath137 percent of the particles in the shell beyond @xmath138 , will lower the entropy per particle by a factor @xmath139 to @xmath140 and the dimensionless temperature by a factor @xmath141 to @xmath142 . in turn , for an initial temperature @xmath143 ( blue dot - dashed curve ) where the cloud is a normal gas , removing the same amount of particles will reduce the entropy only by a factor @xmath144 . removing atoms in the outer shell
repeatedly thus provides an effective tool to reach very low entropies and temperatures . in practice , for atoms confined in an optical dipole trap , this may be achieved by lowering the depth of the optical trap , as was done e.g.   previously in experiments realizing a condensate of fermionic dimers @xcite .    from the local profiles
@xmath15 , @xmath19 , and @xmath18 , it is straightforward to calculate the the total energy @xmath1 and the global entropy @xmath145 by evaluating the integrals in eqs .   and
. in particular , we may eliminate the local value @xmath146 of the dimensionless temperature in the trap center , to obtain the function @xmath147 .
our result for the ultracold fermi gas in the unitary regime is shown in fig .
[ fig : entropy ] as blue - green - red solid line . to distinguish the superfluid and the normal - fluid parts of the curve
, the solid line is shown in blue or red color , respectively .
apparently , this curve is continuous and there is no particular feature , which indicates the superfluid transition .
in fact , since @xmath148 is just the inverse temperature , this behavior is expected not only in a trap , where none of the thermodynamic functions exhibits a singularity , but even in the homogeneous gas , because the superfluid transition is continuous . as mentioned before
, however , our theory predicts a weak first - order transition . the solid line is therefore multivalued in the intervals @xmath149 and @xmath150 , which is indicated by the green section of the line . within these intervals
the superfluid transition is located . in practice , evidently , the multivaluedness is so tiny that it can hardly be seen in fig .
[ fig : entropy ] . in the multivalued region
the blue , green , and red branch of the solid line therefore lie on top of teach other .
bulgac _ et al .
_  have obtained @xmath0 from their monte - carlo calculation @xcite , which , in fig .
[ fig : entropy ] , is shown as orange dashed line .
clearly , the agreement with our theory is nearly perfect for high temperatures well above the superfluid transition .
however , for lower temperatures close to and below the transition , the results differ .
the deviations are largest for zero temperature where @xmath151 .
they are essentially due to the fact that the ground - state energy @xmath152 in both approaches differ . indeed , from eq .
we obtain @xmath153 while the corresponding value of bulgac _ et al . _
@xmath154 is larger because of the larger value @xmath131 of the bertsch parameter . apart from the deviations in the limit of zero temperature
, the theories also differ in their predictions of the behavior near the superfluid transition temperature .
in particular , bulgac _ et al . _
@xcite have calculated the critical values @xmath155 and @xmath156 at the superfluid transition .
these results are considerably larger than our predictions @xmath157 and @xmath158 , whose uncertainty is due to the multivaluedness near @xmath159 . as a result , the value of the critical temperature @xmath160 of the unitary fermi gas in a trap in the theory of ref .
is considerably larger than our prediction @xmath161 .     as a function of the total energy @xmath1 for the unitary fermi gas in a harmonic trap : present theory ( blue - green - red solid ) , nsr theory of hu _ et al . _
@xcite ( magenta dot - dashed ) monte - carlo simulation of bulgac _ et al . _
@xcite ( orange dashed ) , and experimental data of luo _ et al . _
@xcite ( green data points).,scaledwidth=44.0% ]    now , as pointed out above , our result for the critical temperature of the uniform gas agrees rather well with the most precise numerical calculations of this quantity @xcite . in addition , it is also consistent with the recent calculations of bulgac _ et al . _
@xcite , which indicate that @xmath162 is less or equal to @xmath163 , again in very good agreement with our numbers .
this underlines that our self - consistent , conserving theory of the bcs - bec crossover @xcite provides a quantitatively reliable description of the thermodynamics of a balanced fermi gas near unitarity , despite the problems with the first order nature of the transition in our approach and the absence of a small expansion parameter . within lda , which provides a rather accurate description in the relevant regime of particle numbers @xmath164
, our critical temperature @xmath165 for the trapped gas is therefore expected to be close to the exact result .
this is consistent with a very recent analysis of the experimental data of ref .  , which indicates a critical temperature very close to this value @xcite .
previous , much higher values of the critical temperature for both the homogeneous or the trapped gas that were obtained in different extensions of the theory by nozires and schmitt - rink @xcite and also in numerical calculations @xcite are clearly ruled out .
an alternative approach to the bcs to bec crossover problem has been developed by hu _
et al . _
@xcite for the homogeneous system and applied to the ultracold fermi gas in a trap @xcite .
this theory is an extension of the approach by nozires and schmitt - rink ( nsr ) @xcite to the superfluid region at low temperatures .
while the order parameter @xmath166 is determined by the standard mean - field gap equation , the chemical potential @xmath167 is calculated by a particle - density equation , which includes condensed and noncondensed bound pairs .
the extended nsr approach is in fact a limiting case of our theory in which the full , self - consistently determined green functions are replaced by their zeroth order form obtained in bcs theory .
et al . _
@xcite have calculated @xmath147 and obtain a result which is shown as magenta dot - dashed line in fig .
[ fig : entropy ] .
again , the agreement is very good for high temperatures in the normal fluid region while , however , for low temperatures in the superfluid region there are deviations .
et al . _
@xcite obtain the ground - state energy @xmath168 which is related to the bertsch parameter @xmath73 .
moreover , similar to the theory of bulgac _ et al . _
@xcite their critical temperature @xmath169 and entropy @xmath170 are considerably larger than our values .
more recently , hu _ et al . _
@xcite have published a self - consistent result for the entropy @xmath0 which they call `` @xmath171 approximation '' .
this method is equivalent to our approach for the homogeneous system @xcite and gives results that are nearly identical to the prediction of our present theory including the critical temperature @xmath172 .    while in the experimental setup
the correct trap potential is gaussian , all the curves shown in fig .
[ fig : entropy ] are calculated for the harmonic potential .
slight differences would occur in the high - temperature regime where the energies @xmath1 and the entropies @xmath145 are large .
however , it is important that all curves are calculated for the same potential so that they all converge to a single line in the high - temperature limit .
the entropy versus the total energy has been measured experimentally for ultracold @xmath173 atoms in an optical trap by luo _
et al . _
the data are shown as green points with horizontal and vertical error bars .
clearly the data agree with all theories in the high - temperature regime .
the slightly larger entropy values of the experimental data for large energies may be due to the fact that the experiment is performed for a gaussian potential while the theoretical curves are calculated for a harmonic potential .
apparently , for low temperatures and low entropies the experimental data agree better with the theories of bulgac _ et al . _
@xcite and of hu _ et al .
_  @xcite than with our theory .
in particular , an extrapolation to zero entropy gives a ground - state energy @xmath174 , considerably higher than our value @xmath175 .
it is difficult , however , to quantify the error involved in such an extrapolation because the determination of the entropy from a comparison with its ideal fermi gas limit reached after an adiabatic ramp to magnetic fields @xmath176 far on the bcs side of the crossover becomes increasingly difficult as @xmath145 approaches zero .
as pointed out above , a quite sensitive parameter which distinguishes previous theories from our present one is the critical temperature of the superfluid transition and the associated value of the entropy and energy .
unfortunately , it is difficult to determine the critical temperature and thus also the corresponding value of the entropy from measurements of @xmath0 .
based on the excellent agreement of our value for @xmath159 and the bertsch parameter @xmath54 with the most precise numerical or field - theoretical results and the fact that all thermodynamic relations are obeyed at the few percent level , it is likely that our present theory gives a reliable description of the thermodynamics of the trapped unitary gas , despite the fact that , superficially , the agreement with the data shown in fig .
[ fig : entropy ] is not as good as those of previous theories .
the most important parameter for thermodynamic properties is - of course - the temperature , which unfortunately can not be measured directly . in principle , this is possible from the density profile @xmath15 which changes as a function of temperature . for balanced gases , however , this method is not very reliable .
indeed , for low temperatures , the density converges to the simple thomas - fermi profile . as shown in fig .
[ fig : density_profiles ] , the deviations from such a profile at finite temperatures are extremely small for the inner part of the atom cloud ( see the red and brown low - temperature curves , which are nearly identical to the black zero - temperature curve ) . only close to the surface @xmath177 , small variations with the temperature are observed .
therefore , the signal to noise ratio will be small in the interesting regime below @xmath178 .    a quantity which is much more sensitive to temperature is the entropy density , shown in fig .
[ fig : entropy_profiles ] . at low temperatures
, it is peaked near the surface of the atom cloud .
unfortunately , the local entropy density @xmath179 is not accessible experimentally .
we therefore consider the total entropy @xmath145 , which has been measured by @xcite , as discussed above .    in fig .
[ fig : entropy - temperature ] , we plot the total entropy @xmath145 in a trap as a function for the temperature @xmath37 in units of the fermi temperature @xmath180 of the trapped non - interacting gas . as in fig .  [
fig : entropy ] , we indicate the superfluid and the normal - fluid region by blue and red color of the curve , respectively . again , a tiny multivalued region is observed close to the superfluid transition which is indicated by the green section of the line . however , the three branches in the multivalued region are nearly on top of each other . consequently , the entropy @xmath181 is effectively continuous at the superfluid transition as expected . using the curve in fig .
[ fig : entropy - temperature ] , the measurements of the entropy by luo _
et al . _
@xcite therefore allow to reliably infer the related temperatures @xmath37 .     as a function of the temperature @xmath37 for the unitary fermi gas in a harmonic trap .
the blue , red , and green sections represent the superfluid , normal - fluid and multivalued branches of the curve , respectively.,scaledwidth=44.0% ]    alternatively , we may consider the excess internal energy density @xmath182 n(r)$ ] where @xmath183 is the ground - state internal energy per particle at zero temperature .
the related profiles are similar to those shown in fig .
[ fig : entropy_profiles ] . for low temperatures
the excess internal energy density is peaked close to the surface of the atom cloud .
integrating over the space we obtain the excess internal energy @xmath184 which is related to the excess total energy by the virial theorem @xmath185 .
moreover , the total energy @xmath1 is related to the mean square radius of the atom cloud @xmath186 according to @xmath187 . since @xmath186 can be measured in a model - independent way @xcite , the total energy @xmath188 as a function of the temperature @xmath37 provides an alternative method to determine the temperature of the interacting fermi gas .
we obtain a curve which is similar like fig .
[ fig : entropy - temperature ] .
however , the major drawback of this method is that it requires knowledge of the ground state energy @xmath152 , which must be determined with sufficient accuracy . as evident from fig .
[ fig : entropy ] , there are significant discrepancies in the ground state energy @xmath152 for different experiments and theories .
in the superfluid regime at very low temperatures the fermionic quasiparticles are frozen out and do not contribute to thermodynamic quantities because they have an energy gap which is related to the binding energy of the cooper pairs
. however , the spontaneous symmetry breaking implies a gapless goldstone mode which is the bogoliubov - anderson mode .
this mode propagates with a constant velocity @xmath189 like phonons . for this reason ,
the low temperature behavior of the entropy density and the internal energy density is ruled by the well known stefan - boltzmann formulas @xmath190 n(\mathbf{r } ) & = &   2\ , \frac{\sigma(\mathbf{r})}{c(\mathbf{r})}\ , t^4   \label{eq : stefan_boltzmann_energy}\end{aligned}\ ] ] for phonons with one polarization degree of freedom where @xmath191 ^ 2)$ ] is the stefan - boltzmann factor .
since the sound velocity @xmath192 depends on the local particle density , @xmath193 and @xmath192 are space - dependent parameters .
eliminating the temperature we obtain a relation between the local entropy per particle @xmath18 and the local internal energy per particle @xmath19 , which can be written in the form @xmath194^{3/4 } \ .
\label{eq : local_asymptotic_formula}\end{aligned}\ ] ] here @xmath195 is the local fermi velocity , and @xmath196 is the ground - state energy per particle . the ratio @xmath197 is constant and related to the bertsch parameter @xmath59 .    in our previous publication
@xcite we have calculated all thermodynamic quantities for the homogeneous system . hence , also the function @xmath198 is available , so that we can check the asymptotic formula for the homogeneous system .
it turns out that the resulting exponent at low energy is indeed equal to @xmath199 .
however the amplitude , which is determined fully by the sound velocity @xmath189 , gives @xmath200 for the unitary fermi gas .
this is about a factor of two larger than the expected value @xmath201 from the bertsch parameter .
this discrepancy indicates , that the accuracy of our theory at very low temperatures is not sufficient to extract a reliable value of the sound velocity from the entropy .    in the trapped case ,
unfortunately , asymptotic formulas like - do not hold for the total entropy @xmath145 and the total energy @xmath1 . to see this , we integrate eqs .   and over the whole space .
using eqs .
and we obtain well defined results @xmath145 and @xmath202 for the left hand sides , respectively .
however , while the temperature @xmath37 is constant , the space dependence on the right - hand sides arise from the factor @xmath203^{-3 }   \sim [ v_f(\mathbf{r})]^{-3 } \sim [ n(\mathbf{r})]^{-1}$ ] .
this factor is minimum at the center of the trap but maximum at the surface of the atom cloud .
thus , the main contribution of the integral arises from the surface of the atom cloud ( see also fig .  [ fig : entropy_profiles ] ) . here
the fermi gas is a normal fluid , and the dimensionless temperature @xmath204 is large , so that the stefan - boltzmann formulas are not applicable .
consequently , an asymptotic formula like can not be derived for @xmath145 and @xmath205 of the whole trap .     of the function @xmath0 where @xmath1 is two times the potential energy ( green dashed ) , two times the internal energy ( red dot - dashed ) , and one times the total energy ( blue solid ) .
here @xmath152 is the ground - state energy for @xmath2 and @xmath151.,scaledwidth=44.0% ]    empirically , it has been found in the experiments by luo _
_  @xcite that the total entropy at low energies varies with an effective power law @xmath206^b   \label{eq : asymptotic_formula}\ ] ] with an exponent @xmath207 . in order to check
whether such a behavior is consistent with a microscopic theory , we determine the exponent @xmath208 by logarithmic differentiation of the blue - green - red solid curve @xmath147 in fig .  [ fig : entropy ] . the result is shown as blue solid curve in fig .
[ fig : exponents ] and confirms that @xmath0 in a trap indeed exhibits a power law behavior in the regime near the ground state .
we have adjusted the ground - state energy @xmath152 by fine tuning in order to obtain a well defined exponent in the limit @xmath209 .
the resulting exponent @xmath210 is surprisingly close to the value @xmath211 of the local asymptotic formula but differs from the value @xmath212 inferred from the experimental fit .    in fig .
[ fig : exponents ] the superfluid transition is located at the position @xmath213=-0.65 $ ] .
the left part of the blue solid curve corresponds to the superfluid region . here
the exponent is nearly constant up to the superfluid transition . even though the stefan - boltzmann formulas are not valid , the exponent @xmath210 is remarkably close to the theoretical value @xmath211 .
the right part of the blue solid curve corresponds to the normal fluid region . here
the logarithmic derivative of the entropy with respect to energy decreases monotonically with increasing energy @xmath1 and thus no well defined exponent can be attributed to the high energy part of the curve .
the virial theorem implies the energy relations @xmath214 . in order to check the validity of the virial theorem we consider the entropy functions @xmath215 and @xmath216 and calculate the related exponents @xmath217 and @xmath218 , which are shown in fig .
[ fig : exponents ] as green dashed line and as red dot - dashed line , respectively .
these lines should be compared with the blue solid line , which represents the exponent @xmath219 . in the normal - fluid region
the virial theorem is well satisfied , because the right parts of the curves are lying nearly on top of each other where the small deviations are numerical errors .    in the superfluid region
the virial theorem is satisfied only approximately because of the modification of the theory in order to have a gapless bogoliubov - anderson mode ( see sec .
ii.j in ref .  ) . as a consequence ,
the left parts of the curves deviate from each other .
we find three different values for the , exponents which are @xmath220 , @xmath221 , and @xmath222 .
these results are related to the three different ground - state energies , , and , respectively .
based on our previous results for the bcs to bec crossover problem in a homogeneous gas @xcite , we have calculated density and entropy profiles in a trap within a local density approximation . in addition , we have determined the total entropy @xmath145 and energy @xmath1 in the unitary regime and have compared our results with both experiment and recent theories . for temperatures above the superfluid transition temperature ,
very good agreement is obtained .
however , the value of the critical temperature and the behavior at very low temperatures differ appreciably from the experimental estimates and their theoretical analysis in earlier work .    first of all , our value for the bertsch parameter @xmath54 and thus the ground - state energy @xmath152 is about @xmath223 smaller than the results obtained from variational monte carlo calculations or from gaussian fluctuation theories around the bcs ansatz for the ground state .
while our value @xmath224 agrees well with the most precise results so far obtained from an @xmath225 expansion @xcite , it differs from those obtained from variational monte carlo calculations @xcite , or from those including gaussian fluctuations around the bcs mean field @xcite , where @xmath72 or @xmath73 , respectively .
given the uncertainty in determining @xmath54 experimentally ( which requires an extrapolation to zero temperature ) it is clearly important for theory to make precise predictions for @xmath54 which do not rely on approximations that are apparently limiting all present results .
in view of the fundamental importance of this parameter in the context of strongly interacting fermi gases , progress here would be highly desirable .    as a second point ,
our values for both the critical temperature and value of entropy at @xmath159 are appreciably lower than those obtained in the theories of bulgac _ et al . _
@xcite and of hu _ et al . _
@xcite and also those inferred from the original analysis of the experimental data @xcite .
now , as is evident from fig .
[ fig : entropy ] , a measurement of the function @xmath0 does not provide a sensitive measure of the critical temperature .
quantitative results for the critical temperature of the unitary gas have been obtained by shin _
et al . _
they rely on using gases with a finite imbalance @xmath226 , which have always a fully polarized outer shell in a trap .
since a single species ultracold fermi gas is noninteracting , its temperature can be reliably determined from cloud profiles .
an extrapolation back to zero imbalance gives a critical temperature at unitarity which is close to the value predicted both in our theory and in monte - carlo calculations by burovski _
et al . _
@xcite and bulgac _ et al . _
as pointed out in section iv , there is now evidence that our result @xmath227 for the critical temperature of the unitary gas in a trap is rather precise . together with the entropy - temperature curve shown in fig .
[ fig : entropy - temperature ] , this would allow doing precise thermometry for balanced gases , that has not been possible so far .    finally , we have shown that for temperatures of order @xmath228 which have been reached experimentally @xcite , an efficient way of further cooling the gas is possible by removing the high entropy outer part of the atomic cloud and readjusting the confining potential .
this method is similar in spirit than the standard evaporative cooling idea , but potentially much more efficient .
it might open the avenue to reach regimes in which the entropy per particle is much less than @xmath229 , a necessary condition for realizing many of the nontrivially ordered states that are in principle accessible with ultracold fermions @xcite .
we would like to thank a.  bulgac , p.d .
drummond , j.e .
drut and j.  thomas for useful discussions and for making their results available for comparison with our theory . w.z .
is grateful for the hospitality as a visitor at the mit - harvard center for ultracold atoms and for a number of helpful discussions with m.  zwierlein .
this work was supported by the dfg forschergruppe 801 on ultracold quantum gases .
w. ketterle and m.w .
zwierlein , in : _ proceedings of the international school of physics `` enrico fermi '' , course clxiv , varenna , 20 - 30 june 2006 _ , edited by m. inguscio , w. ketterle , and c. salomon ( ios press , amsterdam 2008 ) .
the possibility to treat the unitary fermi gas problem by a systematic expansion around the upper and lower critical dimensions @xmath230 and @xmath231 was pointed out originally by z. nussinov and s. nussinov , phys .
a * 74 * , 053622 ( 2006 ) [ arxiv : cond - mat/0410597 ] , and has been carried out by y. nishida and d. t. son , phys .
lett . * 97 * , 050403 ( 2006 ) .
w. ketterle and n.j .
van druten , in _ advances in atomic , molecular and optical physics _
* 36 * ( 1997 ) , eds .
b. bederson and h. walther . for
an introduction see also chapter vi of the lectures notes by j. dalibard available from http://www.phys.ens.fr/@xmath232dalibard/notes_de_cours/ dea_atomes_froids.pdf h. hu , x.j .
liu , and p.d .
drummond , phys .
a * 77 * , 061605(r ) ( 2008 ) .