We show explicitly that the ground-state density of a system of noninteracting fermions in a harmonic-oscillator potential well can be given in terms of the wave functions of the single-particle states around the Fermi level.
We present an ab initio calculation of small numbers of trapped, strongly interacting fermions using the Green's function Monte Carlo method. The ground-state energy, density profile, and pairing gap are calculated for particle numbers N=2−−22 using the parameter-free ``unitary'' interaction. Trial wave functions are taken in the form of correlated pairs in a harmonic oscillator basis. We find that the lowest energies are obtained with a minimum explicit pair correlation beyond that needed to exploit the degeneracy of oscillator states. We find that the energies can be well fitted by the expression ${a}_{\mathit{TF}}{E}_{\mathit{TF}}+\ensuremath{\Delta}\phantom{\rule{0.2em}{0ex}}\mathrm{mod}(N,2)$ where ETF is the Thomas-Fermi energy of a noninteracting gas in the trap and $\ensuremath{\Delta}$ is the pairing gap. There is no evidence of a shell correction energy in the systematics, but the density distributions show pronounced shell effects. We find the value $\ensuremath{\Delta}=0.7\ifmmode\pm\else\textpm\fi{}0.2\ensuremath{\omega}$ for the pairing gap. This is smaller than the value found for the uniform gas at a density corresponding to the central density of the trapped gas.
We consider N fermions in a two-dimensional harmonic oscillator potential interacting with a very short-range repulsive pair-wise potential. The ground-state energy of this system is obtained by performing a Thomas-Fermi as well as a self-consistent Hartree-Fock calculation. The two results are shown to agree even for a small number of particles. We next use the finite-temperature Thomas-Fermi method to demonstrate that in the local density approximation, these interacting fermions are equivalent to a system of noninteracting particles obeying the Haldane-Wu fractional exclusion statistics. It is also shown that mapping onto a system of N noninteracting quasiparticles enables us to predict the energies of the ground and excited states of the N-body system. PACS Nos.: 05.30-d, 73.20Dx