Conjectures are made for the ground state energy of a large spin 1/2 Fermion system trapped in a 1D harmonic trap with delta function interaction. States with different spin J are separately studied. The Thomas–Fermi method is used as an effective test for the conjecture.
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 study the ground state energy of a system of N fermions with two spin states in the large N limit. The particles are placed in an inhomogeneous trapping potential and interact via scaled interactions. We study a dilute limit where the range of the interaction potential is much smaller than the typical inter-particle distance. We show that the energy per particle converges to the Thomas-Fermi energy of the system, with a perturbative term corresponding tot he interaction and exhibiting the scattering length of the potential. The proof is decomposed into two bounds. First, we construct an appropriate test-state to prove the upper bound. Then, we prove the lower bound by the Dyson lemma, which allows us to regularize the interaction potential, and several semi-classical tools.