arXiv:quant-ph/9612027v1 7 Dec 1996 fara-pc Frisea.” “Fermi visualization real-space of explicit a distribution an of provide spatial would the cloud of trapped uni- the observation w same As below, the variables. see of by will scaling described suitable after be different the functions, all versal with that can gases find numbers We trapped particle [7] harmonically spin elsewhere. of trapped discussed properties two-or-more be with . will gases of states of function a properties trapped as The a of gas, heat Fermi specific spin-polarized the and potential, chemical sea. the Fermi be- compact crossover the gradual and limit a classical undergoes the gas tween Fermi phase regimes trapped a classical a and case degenerate than Bose the rather the separates principle transition in exclusion While Pauli interactions. the this classical atom-atom however, to fermions, a due for to is [5] temperature; relative effect same Bose expanded the confined the are at the both gas gases [6] , of Fermi low realization and At simple pro- system. traps particularly Fermi harmonic a confined of potential controlled, vide quadratic a The in systems setting. quantum studying for limit. in nucle- weak dilute are predom- the and effects by whose metals interact interactions gases short-range and atomic inantly however, atoms nuclei, in in approxima- electrons ons zeroth-order Unlike well-understood good and a tion. non- is old the gas where an Fermi systems interacting is precur- familiar many gas possible are there Fermi temperatures. a problem; lower ideal as at The condensate and Fermi [4] right paired quan- own a degenerate to its a sor Fermi in as trapped system both attention, of tum discontinuities merits behavior the also the as gases bosons, dramatic with as associated not di- Al- perhaps strikingly [1–3] the condensation. though Bose-Einstein as of such observation effects, rect study statistical quantitative quantum and of manipulation the for portunities ecluaetesailadmmnu distributions, momentum and spatial the calculate We laboratory excellent an provide gases atomic Trapped op- exciting provide gases atomic degenerate Trapped esuytepoete fasi-oaie em a nahar Bose a dilute a in of gas properties m Fermi corresponding and the spatial spin-polarized with the a compared for forms of Universal approximation. properties Fermi) the study We .INTRODUCTION I. † oe etrfrCmlxSses rnesUiest,Wal University, Brandeis Systems, Complex for Center Volen ∗ eateto hsc,Uiest fClfri,Berkeley California, of University Physics, of Department ..Butts D.A. rpe em gases Fermi Trapped Dcme ,1996) 6, (December ∗ e , ..Rokhsar D.S. , 1 h equation the ihasnl-atceHamiltonian single-particle potential, a harmonic with symmetric azimuthally an in moving nt o nrisls hnteFrienergy Fermi the than less is energies factor for occupation unity Fermi-Dirac the temperature zero At where hc estecaatrsi nryo h tmccloud. atomic the of energy characteristic the sets which 0 h igepril eeso q 1 r familiar: are (1) eq. of levels ne- single-particle approximation below). The (see first them between a interactions to short-ranged glect the may we Since fermions, spi polarized respectively. of close-approach prohibits directions, principle axial exclusion Pauli and radial the fe.()te ie [6] gives then (4) eq. of rae hntelvlsaig( spacing level the temperatures For than greater suppressed. been has energy zero-point est fstates with of continuum a density as spectrum single-particle discrete this where N , mnu itiuin r acltd n h results the and calculated, are distributions omentum H gas. econsider We h hmclpotential chemical The ( r ,adzr tews.Asrihfradintegration straightforward A otherwise. zero and ), , p n I.EEG N EGHSCALES LENGTH AND ENERGY III. oi rp sn h eilsia (Thomas- semiclassical the using trap, monic ω = ) x r , ∗ n and ǫ , † n y 2 x and , M 1 I EST FSTATES OF DENSITY II. ,n y ω [ N N p ,n E z A94720-7300 CA , g x 2 z F ( n = hm A02254 MA tham, + pnplrzdfrin fmass of fermions spin-polarized = E ¯ = z ¯ = Z = ) p r o-eaieitgr,adthe and integers, non-negative are hω y 2 λω hω + e r β r 2 r g [ µ ( p [6 n λ ( E ( z 2 r h rpfeunisin frequencies trap the are x (¯ E ,N T, λN − E + ] hω + ) µ k 2 dE ) B r ] n Mω 1 + 1 ) T sgvnipiil by implicitly given is ) y / 3 3 2 . + , ≫ . r 2 λn [ ¯ hω x 2 z ] ,w a treat may we ), + , y E 2 F + ≡ λ µ 2 z ( T 2 ] (4) (3) (1) (5) (2) M , n- = The characteristic size [6] of the trapped degenerate is given by RF , the excursion of a classical particle with total energy EF in the trap potential:

R [2E /Mω2]1/2 = (48Nλ)1/6σ . (6) F ≡ F r r 1/2 where σr = (¯h/Mωr) is the radial width of the Gaus- sian of the trap. For large N, the width of the degenerate Fermi cloud is much greater than the quantum length σr, and the Fermi energy is much greater than the level spacing of the trap, due to the Pauli ex- clusion induced “repulsion” between fermions. [8,9] Similarly, we may define a characteristic wavenumber, KF , which is determined by the momentum of a free particle of energy EF :

2 1/2 1/6 −1 KF [2MEF /¯h ] = (48Nλ) σ (7) vs. ≡ r FIG. 1. The temperature. Both axes 3 1/3 are scaled by the Fermi energy, which results in a universal = (48Nλ/RF ) . (8) curve that applies to all harmonically trapped Fermi gases.

From eq. (8) we see that KF is roughly the reciprocal of the typical interparticle spacing in the gas. In eqs. (9) and (10), the particle number N enters 6 As an example, consider spin-polarized Li. The radial µ(T,N) only through the Fermi energy EF . This result frequency of this fermionic isotope of lithium in the TOP holds generally for all temperatures, as can be seen by −1 trap of ref. [1] would be ωr = 3800 sec ; the character- casting eq. (4) in dimensionless form by scaling E, µ, and istic ground state length σr is 1.6 µm. The TOP trap 1/β by EF . (The same conclusion holds for any density has an intrinsic axial/radial ratio λ = √8. For N = 105 of states of the form g(E) = AEb, for constant A, b.) atoms, the radius RF is 25 µm; the typical interparticle Figure 1 is therefore a universal curve in the sense that spacing 1/KF is 0.1 µm. The Fermi temperature for this it applies to harmonically trapped Fermi gases containing gas would be 3.5 µK. any number of particles. The specific heat [10] per particle of the trapped Fermi gas, 1/N ∂U/∂T , is shown in figure 2. It is a mono- |N IV. CHEMICAL POTENTIAL AND SPECIFIC tonic function of temperature. For low temperatures the 2 HEAT VS. TEMPERATURE specific heat per particle is π kB (kBT/EF ); at high tem- peratures C = 3NkB, as we approach the equipartition For general temperature, the chemical potential µ must limit for particles in a harmonic potential. be determined numerically using eq. (4). We can find analytic expressions, however, in the limits of high and low temperature. For low temperature (kB T EF ) the chemical potential is given by the Sommerfeld≪ expansion:

2 π kBT 2 µ(T,N)= EF [1 ( ) ]. (9) − 3 EF The third and higher-order terms in the Sommerfeld se- ries vanish since the density-of-states is a quadratic func- tion of energy. At high temperatures (i.e., in the classical limit k T E ), we find B ≫ F

kBT 3 µ(T,N)= kBT ln[6( ) ]. (10) − EF

Numerical results for µ(T )/EF are compared with these two limiting forms in figure 1. Evidently the FIG. 2. The heat capacity vs. temperature. The heat ca- low temperature approximation is quantitatively accu- pacity is scaled by kBN and the temperature by EF . The rate below kBT/EF 0.5, and the classical expression classical result is shown by the dotted line. holds above k T/E ∼ 0.6. B F ∼

2 k (r)3 V. SEMICLASSICAL (THOMAS-FERMI) n(r; T =0)= F . (16) APPROXIMATION 6π2 Note that n(r) vanishes for ρ > R , where ρ is the effec- Since the exact eigenstates of the harmonic potential F tive distance are well-known, the properties of a harmonically trapped can in principle be found directly by summing ρ [x2 + y2 + λ2z2]1/2. (17) over these states. It is useful, however, to have approxi- ≡ mate forms for various observables that can be computed Combining eqs. (15) and (16), we find [6,12] directly in the large N limit, where the exact sums be- come unwieldy. Nλ 8 ρ2 n(r; T =0)= [1 ]3/2 (18) In the “semiclassical” or Thomas-Fermi approx- 3 2 2 RF π − RF imation, [11] the state of each atom is labeled by a po- sition r and a wavevector k, which can be viewed as the for ρ R . The cloud encompasses an ellipsoid with ≤ F centers of a wavepacket state. The energy of the parti- diameter 2RF in the xy plane, and diameter 2RF /λ along cle is simply the corresponding value of the Hamiltonian; the z axis. This aspect ratio is the same as that of a the density of states in the six-dimensional phase space classical gas in the same potential, since the Boltzmann −3 2 2 (r, k) is (2π) , where sums over states are replaced by distribution is proportional to exp[ Mω ρ /2kBT ]. integrals over phase space. These semiclassical approxi- − mations are valid in the limit of large N, as discussed in the Appendix. VII. MOMENTUM DISTRIBUTION AT ZERO In the semiclassical limit, the number density in phase TEMPERATURE space is One way to characterize the state of a trapped gas is to 1 1 r k allow a rapid adiabatic expansion and then measure the p( , ; T, µ)= 3 β(H(r,h¯k)−µ) . (11) (2π) e +1 velocity distribution by time-of-flight spectroscopy. [1] For the degenerate , the observed anisotropy of The chemical potential is given implicitly by the require- this velocity distribution is dramatic evidence for quan- ment tum statistical effects. The semiclassical momentum dis- tribution for a degenerate Fermi gas at zero-temperature N = d3r d3k p(r, k; T, µ). (12) Z is simply 1 It follows from the correspondence principle that the n˜(k; T =0)= d3r Θ(k (r) k ), (19) Thomas-Fermi calculation of µ(T,N) using eq. (12) re- (2π)3 Z F − | | produces the exact result obtained from eq. (4); this is easily confirmed for the harmonic oscillator, since the two where Θ(kF (r) k ) is the unit step function. The inte- − | | integrals are related by a simple change-of-variables. gral (19) is the real-space volume within which the local After computing µ(T,N), it is straightforward to cal- Fermi wavevector exceeds k : | | culate the spatial and momentum distribution functions N 8 k 2 k 3/2 n˜( ; T =0)= 3 2 [1 | |2 ] , (20) n(r; T )= d3k p(r, k; T, µ) (13) KF π − KF Z where the maximum occupied wavenumber KF was de- n˜(k; T )= d3r p(r, k; T, µ). (14) Z fined in eq. (7). Note that kF (r =0)= KF . Despite the spatial anisotropy of the trap, the momen- tum distribution of the degenerate Fermi gas is isotropic. This isotropy is a general feature of trapped Fermi gases, VI. SPATIAL DISTRIBUTION AT ZERO independent of trap potential, [13] since from eq. (19) we TEMPERATURE see thatn ˜(k) depends only on the magnitude of k. The spatial and momentum distributions (18, 20) both At zero temperature, we may define a “local” Fermi have the same functional form, because is a quadratic function of both position and momentum.H In this sense, wavenumber kF (r) by the distribution (18) can be viewed as a Fermi sea in 2 2 r ¯h kF (r) real-space. The anisotropy of n( ) is due to the unequal + V (r)= EF , (15) spring constants, while the isotropy ofn ˜(k) is due to 2M 2 2 2 the isotropy of mass. That is, px, py, and pz enter the where V (r) is the trap potential. The density n(r) is Hamiltonian with the same coefficient, while x2, y2, and then simply the volume of the local Fermi sea in k-space, z2 need not. multiplied by the density of states (2π)−3:

3 VIII. NUMERICAL RESULTS

In the semiclassical approximation, the spatial and mo- mentum distributions are easily determined numerically for any temperature as described in section V. As with the chemical potential, an appropriate scaling of these two distributions yields a universal form for all harmon- ically trapped Fermi gases when plotted vs. the scaled variables ρ/RF and k /KF , respectively. Figure 3 shows the| | scaled density vs. scaled distance for kBT/EF of 0, 0.25, 0.5, 0.75, and 1. At low tem- peratures, the density is close to their zero temperature form, with a thin evaporated “atmosphere” of thickness RF (kB T/EF ) surrounding a degenerate liquid “core.” ∼ In the classical limit, the density approaches a Gaus- sian in ρ, with a width given by the equipartition theo- rem: ρ2 =3R2 (k T/E ). As shown in figure 3, this ac- FIG. 4. The mean-square variation of the size of the cloud F B F interpolates between a low-temperature degenerate regime curatelyh i describes the density distribution for k T/E = B F and a high-temperature limit that is well-described by the 1. The evolution of the density profile from its low- equipartition theorem. temperature Fermi form (18) to the classical limit can be tracked by calculating the mean-square excursion ρ2 , which is shown in figure 4 in the dimensionless formh i 2 2 IX. PERTURBATIONS ρ /RF vs. kBT/EF . This is again a universal curve forh i all harmonically trapped Fermi gases. As noted above for zero temperature, the momentum What happens if the potential is not perfectly har- r distribution at non-zero temperature has the same form monic? We may treat δV ( ) as a perturbation. Here we as the spatial distribution, since momentum and position focus our attention on the T = 0 case. From eq. (15), a both enter the single particle Hamiltonian quadratically. change in trap potential shifts the local Fermi wavenum- 3 ber by Thus the scaled momentum distribution KF n˜(k)/N vs. k /KF is also given by figure 3. Similarly, figure 4 illus- M | | 2 2 trates the scaled mean-square momentum k /K vs. δkF (r)= [δEF δV (r)], (21) F ¯hk (r) − k /K as well. h i F | | F where δEF is the change in Fermi energy. From eq. (16), the corresponding change in density is

Mk (r) δn(r)= F [δE δV (r)], (22) 2π2¯h F − where the Fermi energy is adjusted to make d3r δn(r) vanish: R

3 d rδV (r)kF (r) δEF = 3 . (23) R d rkF (r) R Interactions between particles (which are small for spin-aligned gases) can be treated as a perturbation δV (r) = Un(r), where U is a parameter related to the p-wave scattering length.

FIG. 3. The universal spatial and momentum density dis- X. COMPARISON WITH THE BOSE GAS tributions for kBT/EF = 0 (bold), 0.25, 0.5, 0.75, and 1.0. The classical result for kB T/EF = 1 is shown as a dashed line. The interacting Bose gases of refs. [1] and [2] are in the Thomas-Fermi regime [14]. Since the gases remain dilute, two-body scattering may be treated by a delta-function pseudopotential of strength U = 4π¯h2a/M, where a is the s-wave scattering length. When the dimensionless 3 parameter UN/¯hωλσr is large (as is appropriate for the

4 experiments of refs. [1] and [2]) the density profile of the XII. APPENDIX: VALIDITY OF THE interacting Bose gas is [14] SEMICLASSICAL APPROXIMATION

R2 ρ2 n (r)= B 1 , (24) The semiclassical approximation can be safely applied B 2U R2  − B  to an inhomogeneous Fermi gas of density n(r) if we can imagine partitioning the system into cells of linear di- with maximum radius mension ℓ such that the following two conditions are si- multaneously met: 15λUN 1/5 RB = ( ) . (25) (1) The number of particles in a cell is much greater 4π than unity, so that a local Fermi sea may be envisioned: Note that the characteristic radius scales more slowly with particle number for the Fermi gas (N 1/6) than for n(r)ℓ3 1. (26) ≫ the interacting Bose gas (N 1/5); similarly, the Fermi en- ergy scales as N 1/3 while the zero temperature chemical (2) The variation of the trap potential across the cell potential of the Bose gas varies more rapidly, as N 2/5. (ℓ V ) must be small compared with the local Fermi en- ∇ 2 2 The axial/radial aspect ratio for both classical and de- ergyh ¯ kF (r) /2M, so that within a cell the potential generate trapped gases is λ, since in every case the den- energy is nearly constant. At low temperature, this con- sities are functions of ρ only. The velocity (momentum) dition becomes distributions, however, can be quite different: for classi- ¯h2 cal and Fermi gases the velocity distribution is isotropic, ℓMω2r [6π2n(r)]2/3, (27) whilen ˜(k) for a zero temperature Bose gas is the square ≪ 2M of the Fourier transform of nB(r). This is notably where we have used eq. (16). anisotropic. Note that as N increases,p both RF and RB Combining eqs. (26) and (27), we see that to be able increase, but the widths of the respective momentum dis- to choose a (possibly r-dependent) cell size ℓ that simul- tributions go in opposite directions: KF increases with taneously satisfies these two conditions, the number of N, while the typical momentum of a particle in a trapped particles per quantum volume must satisfy Bose condensate decreases with particle number, since r KB 1/RB by the uncertainty principle. n(r)σ3 , (28) It∼ is amusing to compare the interatomic repulsion in ≫ σ a Bose gas with the effective repulsion experienced by 1/2 fermions due to the Pauli exclusion principle. Equat- where σ is as before the quantum length (¯h/Mω) and ing the characteristic Fermi and Bose radii (6, 25) we we have omitted factors of order unity. see that, crudely speaking, the spatial distribution of At low temperatures, the semiclassical density given by eq. (18) scales as N/R3 N 1/2/σ3 near the ori- a degenerate Fermi gas is mimicked by that of a Bose F ∼ gas interacting via an effective “Pauli pseudopotential” gin, so the Thomas-Fermi approximation is always self- U E (R3 /N), which is the characteristic energy consistent at the center of the trap for large N. (This eff F F r multiplied∼ by the volume per particle. Equivalently, the can be confirmed at = 0 by direct summation of the effective scattering length brought about by the Pauli squares of the simple harmonic oscillator eigenfunctions −1 up to energy EF .) Near the periphery of the cloud, how- principle is aeff KF , i.e., the inter-particle spacing. (This is quite expected,∼ since the inter-particle spacing ever, the density becomes small, and the approximation is the only appropriate length in the ideal Fermi gas.) is not valid. It is easy to show that semiclassical treat- The use of such an effective interaction is limited by the ment fails within a shell at the periphery of the cloud, whose thickness δR 1/K σN −1/6 vanishes in the fact that (a) the momentum distributions of the Fermi ∼ F ∼ and Bose gases remain quite different and (b) the gas is limit of large N. Within this shell only the exponential not dilute with respect to the exclusion-induced “inter- tails of a few single particle states contribute to the den- actions” since K a is of order unity. sity; this is analogous to the corresponding region of the F eff Bose gas. [15]

XI. ACKNOWLEDGEMENTS

We thank Mark Kasevich, David Weiss, and Ike Silvera for useful and interesting discussions. This work was sup- ported by the National Science Foundation under grant [1] M.H. Anderson, J.R. Ensher, M.R. Matthews, C.E. NSF-DMR-91-57414, and the Committee on Research at Weiman, and E.A. Cornell, Science 269, 198 (1995). UC Berkeley. Work at Brandeis was supported by the [2] K.B. Davis, M.O. Mewes, M.R. Andrews, N.J. van Sloan Center for Theoretical Neurobiology. Druten, D.S. Durfee, D.M. Kurn, and W. Ketterle, Phys. Rev. Lett. 75, 3969 (1995).

5 [3] C.C. Bradley, C.A. Sackett, J.J. Tollett, and R.G. Hulet, Phys. Rev. Lett. 75, 3969 (1995). [4] H.T.C. Stoof, M. Houbiers, C.A. Sackett, and R.G. Hulet, Phys. Rev. Lett., 76 10 (1996). [5] V.V. Goldman, I.F. Silvera, and A.J. Leggett, Phys. Rev. B 24, 2870 (1981). [6] I.F. Silvera and J.T.M. Walraven, J. Appl. Phys. 52, 2304 (1981). [7] D.A. Butts and D.S. Rokhsar, unpublished. [8] O.O. Tursunov and O.V. Zhirov, Phys. Lett. 222, 110 (1994). [9] B.L.G. Bakker, M.I. Polikarpov, and A.I. Veselov, quant- ph/9511009 (1995). [10] Note that the trapped gas does not occupy a container of externally fixed size, and so we do not introduce separate constant-volume and constant-pressure specific heats. As the gas is heated, it expands against the harmonic po- tential, doing work; the specific heat of the harmonically trapped system is therefore expected to be larger than that of a comparable gas-in-a-box at fixed volume. Since the harmonic restoring force grows linearly with the size of the cloud, the effective pressure it exerts will also vary. [11] L.H. Thomas, Proc. Cambridge Phil. Soc., 23, 542 (1927); E. Fermi, Z. Physik, 48, 73 (1928). [12] I.F. Silvera, Physica 109-110B, 1499 (1982). [13] A more familiar case is the ideal Fermi gas in a Lx ×Ly × Lz rectangular box – the spatial distribution is clearly anisotropic, but the momentum distribution isn’t, since the density of states in k-space is uniformly (2π)−3. [14] G.Baym and C.J. Pethick, Phys. Rev. Lett. 76, 6 (1996). [15] S. Stringari, F. Dafolvo, and L. Pitaevski, preprint.

6