Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:cond-mat/0112248 · Physical Review Letters

Shell Structure and Pairing for Interacting Fermions in a Trap

H. Heiselberg, Ben R. Mottelson

Published Apr 24, 2002• 36 citations
arXivHTMLPDF

Abstract

The shell structures for weakly interacting fermions in harmonic oscillator traps at zero temperature undergo several transitions depending on the number of particles in the trap and their interaction strength. Calculations of the one and two-particle spectra give the pairing gaps, and the many particle state is constructed by the seniority scheme. For sufficiently few particles, N less than or similar to 10(4), the pairing field exceeds the mean-field splitting of single particle levels at the Fermi surface and the greater number of almost degenerate states due to the SU(3) symmetry in the harmonic oscillator shell lead to supergaps. Deformation and rotation are also discussed.

DOI
10.1103/physrevlett.88.190401

Related papers

Pairing of fermions in atomic traps and nuclei

H. Heiselberg

Pairing gaps for fermionic atoms in harmonic-oscillator traps are calculated for a wide range of interaction strengths and particle numbers, and compared to pairing in nuclei. Especially systems where the pairing gap exceeds the single-level spacing but is smaller than the shell splitting $\ensuremath{\Elzxh}\ensuremath{\omega}$ are studied, which applies to most trapped Fermi atomic systems as well as to finite nuclei. When solving the gap equation for a large trap with such multilevel pairing, one finds that the matrix elements between nearby harmonic-oscillator levels and the quasi particle energies lead to a double logarithm of the gap, and a pronounced shell structure at magic numbers. It is argued that neutron and proton pairing in nuclei belongs to the class of multilevel pairing, that their shell structure follows naturally, and that the gaps scale as $\ensuremath{\sim}{A}^{\ensuremath{-}1/3}$---all in qualitative agreement with odd-even staggering of nuclear binding energies. Pairing in large systems is related to that in the bulk limit. For large nuclei the neutron and proton superfluid gaps approach the asymptotic value in infinite nuclear matter: $\ensuremath{\Delta}\ensuremath{\simeq}1.1\mathrm{MeV}.$

arXiv:physics/0304005•Physical Review A•Nov 24, 2003•36 citations

Cooper pairing and single-particle properties of trapped Fermi gases

Georg M. Bruun, H. Heiselberg

We calculate the elementary excitations and pairing of a trapped atomic Fermi gas in the superfluid phase. The level spectra and pairing gaps undergo several transitions as the strength of the interactions between atoms and the number of atoms are varied. For weak interactions, the Cooper pairs are formed between particles residing in the same harmonic oscillator shell. In this regime, the nature of the paired state is shown to depend critically on the position of the chemical potential relative to the harmonic oscillator shells and on the size of the mean field. For stronger interactions, we find a region where pairing occur between time-reversed harmonic oscillator states in different shells also.

arXiv:cond-mat/0203130•Physical Review A•Apr 23, 2002•49 citations

Shell-Model Monte Carlo Simulations of Pairing in Few-Fermion Systems

Nikolaj Thomas Zinner, Klaus Mølmer, C. Özen +2 more

We study a trapped system of fermions with an attractive zero-range two-body interaction using the Shell-Model Monte Carlo method. The method provides {\em ab initio} results in the low $N$ limit where mean-field theory is not applicable. The energy and pairing properties are presented as functions of interaction strength, particle number, and temperature. In the interesting region where typical matrix elements of the two-body interaction are comparable to the level spacing of the trap we find large odd-even effects and signatures of shell structure. As a function of temperature, we observe the disappearance of these effects as in a phase transition.

arXiv:0803.2861•Mar 19, 2008•0 citations

Crossover between few and many fermions in a harmonic trap

Tomasz Grining, Michał Tomza, Michał Lesiuk +5 more

The properties of a balanced two-component Fermi gas in a one-dimensional harmonic trap are studied by means of the coupled-cluster method. For few fermions we recover the results of exact diagonalization, yet with this method we are able to study much larger systems. We compute the energy, the chemical potential, the pairing gap, and the density profile of the trapped clouds, smoothly mapping the crossover between the few-body and many-body limits. The energy is found to converge surprisingly rapidly to the many-body result for every value of the interaction strength. Many more particles are instead needed to give rise to the nonanalytic behavior of the pairing gap, and to smoothen the pronounced even-odd oscillations of the chemical potential induced by the shell structure of the trap.

arXiv:1507.03174•Physical Review A•Dec 3, 2015•51 citations

Unitary Fermi gas in a harmonic trap

S. Y. Chang, G. F. Bertsch

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 ${E}_{\mathit{TF}}$ 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.

arXiv:physics/0703190•Physical Review A•Aug 2, 2007•74 citations

Super-shell structure in harmonically trapped fermionic gases and its semi-classical interpretation

Magnus Ögren, Y Yu, S Åberg +2 more

It was recently shown in self-consistent Hartree-Fock calculations that a harmonically trapped dilute gas of fermionic atoms with a repulsive two-body interaction exhibits a pronounced {\it super-shell} structure: the shell fillings due to the spherical harmonic trapping potential are modulated by a beat mode. This changes the ``magic numbers'' occurring between the beat nodes by half a period. The length and amplitude of the beating mode depends on the strength of the interaction. We give a qualitative interpretation of the beat structure in terms of a semiclassical trace formula that uniformly describes the symmetry breaking U(3) $\to$ SO(3) in a 3D harmonic oscillator potential perturbed by an anharmonic term $\propto r^4$ with arbitrary strength. We show that at low Fermi energies (or particle numbers), the beating gross-shell structure of this system is dominated solely by the two-fold degenerate circular and (diametrically) pendulating orbits.

arXiv:cond-mat/0612659•Physica Scripta•Jun 28, 2006•1 citation
The shell structures for weakly interacting fermions in harmonic oscillator traps at zero temperature undergo several transitions depending on the number of particles in the trap and their interaction strength. Calculations of the one and two-particle spectra give the pairing gaps, and the many particle state is constructed by the seniority scheme. For sufficiently few particles, N less than or similar to 10(4), the pairing field exceeds the mean-field splitting of single particle levels at the Fermi surface and the greater number of almost degenerate states due to the SU(3) symmetry in the harmonic oscillator shell lead to supergaps. Deformation and rotation are also discussed.
Pairing gaps for fermionic atoms in harmonic-oscillator traps are calculated for a wide range of interaction strengths and particle numbers, and compared to pairing in nuclei. Especially systems where the pairing gap exceeds the single-level spacing but is smaller than the shell splitting $\ensuremath{\Elzxh}\ensuremath{\omega}$ are studied, which applies to most trapped Fermi atomic systems as well as to finite nuclei. When solving the gap equation for a large trap with such multilevel pairing, one finds that the matrix elements between nearby harmonic-oscillator levels and the quasi particle energies lead to a double logarithm of the gap, and a pronounced shell structure at magic numbers. It is argued that neutron and proton pairing in nuclei belongs to the class of multilevel pairing, that their shell structure follows naturally, and that the gaps scale as $\ensuremath{\sim}{A}^{\ensuremath{-}1/3}$---all in qualitative agreement with odd-even staggering of nuclear binding energies. Pairing in large systems is related to that in the bulk limit. For large nuclei the neutron and proton superfluid gaps approach the asymptotic value in infinite nuclear matter: $\ensuremath{\Delta}\ensuremath{\simeq}1.1\mathrm{MeV}.$
We calculate the elementary excitations and pairing of a trapped atomic Fermi gas in the superfluid phase. The level spectra and pairing gaps undergo several transitions as the strength of the interactions between atoms and the number of atoms are varied. For weak interactions, the Cooper pairs are formed between particles residing in the same harmonic oscillator shell. In this regime, the nature of the paired state is shown to depend critically on the position of the chemical potential relative to the harmonic oscillator shells and on the size of the mean field. For stronger interactions, we find a region where pairing occur between time-reversed harmonic oscillator states in different shells also.
We study a trapped system of fermions with an attractive zero-range two-body interaction using the Shell-Model Monte Carlo method. The method provides {\em ab initio} results in the low NNN limit where mean-field theory is not applicable. The energy and pairing properties are presented as functions of interaction strength, particle number, and temperature. In the interesting region where typical matrix elements of the two-body interaction are comparable to the level spacing of the trap we find large odd-even effects and signatures of shell structure. As a function of temperature, we observe the disappearance of these effects as in a phase transition.
The properties of a balanced two-component Fermi gas in a one-dimensional harmonic trap are studied by means of the coupled-cluster method. For few fermions we recover the results of exact diagonalization, yet with this method we are able to study much larger systems. We compute the energy, the chemical potential, the pairing gap, and the density profile of the trapped clouds, smoothly mapping the crossover between the few-body and many-body limits. The energy is found to converge surprisingly rapidly to the many-body result for every value of the interaction strength. Many more particles are instead needed to give rise to the nonanalytic behavior of the pairing gap, and to smoothen the pronounced even-odd oscillations of the chemical potential induced by the shell structure of the trap.
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−−22N=2--22N=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{E}_{\mathit{TF}}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.
It was recently shown in self-consistent Hartree-Fock calculations that a harmonically trapped dilute gas of fermionic atoms with a repulsive two-body interaction exhibits a pronounced {\it super-shell} structure: the shell fillings due to the spherical harmonic trapping potential are modulated by a beat mode. This changes the ``magic numbers'' occurring between the beat nodes by half a period. The length and amplitude of the beating mode depends on the strength of the interaction. We give a qualitative interpretation of the beat structure in terms of a semiclassical trace formula that uniformly describes the symmetry breaking U(3) →\to→ SO(3) in a 3D harmonic oscillator potential perturbed by an anharmonic term ∝r4\propto r^4∝r4 with arbitrary strength. We show that at low Fermi energies (or particle numbers), the beating gross-shell structure of this system is dominated solely by the two-fold degenerate circular and (diametrically) pendulating orbits.