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:2506.20891 · SciPost Physics Core

Energetic and structural properties of two-dimensional trapped mesoscopic Fermi gases

Emma Laird, Brendan C. Mulkerin, Jia Wang, Matthew J. Davis

Published Apr 9, 2026• 0 citations
arXivHTMLPDF

Abstract

We theoretically investigate equal-mass spin-balanced two-component Fermi gases in which pairs of atoms with opposite spins interact via a short-range isotropic model potential. We probe the distinction between two-dimensional and quasi-two-dimensional harmonic confinement by tuning the effective range parameter within two-dimensional scattering theory. Our approach, which yields numerically exact energetic and structural properties, combines a correlated Gaussian basis-set expansion with the stochastic variational method. For systems containing up to six particles, we: 1) Present the ground- and excited-state energy spectra; 2) Study non-local correlations by analysing the one- and two-body density matrices, extracting from these the occupation numbers of natural orbitals, the momentum distributions of atoms and pairs, and the molecular `condensate fraction’; 3) Study local correlations by computing the radial and pair distribution functions. This paper extends current theoretical knowledge on the properties of trapped few-fermion systems as realised in state-of-the-art cold-atom experiments.

DOI
10.21468/scipostphyscore.9.2.021

Related papers

Trapped two-component Fermi gases with up to six particles: Energetics, structural properties, and molecular condensate fraction

D. Blume, K. M. Daily

We investigate small equal-mass two-component Fermi gases under external spherically symmetric confinement in which atoms with opposite spins interact through a short-range two-body model potential. We employ a non-perturbative microscopic framework, the stochastic variational approach, and determine the system properties as functions of the interspecies s -wave scattering length a s , the orbital angular momentum L of the system, and the numbers N 1 and N 2 of spin-up and spin-down atoms (with N 1 − N 2 = 0 or 1 and N ⩽ 6 , where N = N 1 + N 2 ). At unitarity, we determine the energies of the five- and six-particle systems for various ranges r 0 of the underlying two-body model potential and extrapolate to the zero-range limit. These energies serve as benchmark results that can be used to validate and assess other numerical approaches. We also present structural properties such as the pair distribution function and the radial density. Furthermore, we analyze the one-body and two-body density matrices. A measure for the molecular condensate fraction is proposed and applied. Our calculations show explicitly that the natural orbitals and the momentum distributions of atomic Fermi gases approach those characteristic for a molecular Bose gas if the s -wave scattering length a s , a s > 0 , is sufficiently small.

arXiv:1008.3191•Comptes Rendus Physique•Jan 1, 2011•28 citations

Energetics and structural properties of trapped two-component Fermi gases

J. von Stecher, Chris H. Greene, D. Blume

Using two different numerical methods, we study the behavior of two-component Fermi gases interacting through short-range $s$-wave interactions in a harmonic trap. A correlated Gaussian basis-set expansion technique is used to determine the energies and structural properties, i.e., the radial one-body densities and pair distribution functions, for small systems with either even or odd $N$, as functions of the $s$-wave scattering length and the mass ratio $\ensuremath{\kappa}$ of the two species. Particular emphasis is put on a discussion of the angular momentum of the system in the BEC-BCS crossover regime. At unitarity, the excitation spectrum of the four-particle system with total angular momentum $L=0$ is calculated as a function of the mass ratio $\ensuremath{\kappa}$. The results are analyzed from a hyperspherical perspective, which offers unique insight into the problem. Additionally, fixed-node diffusion Monte Carlo calculations are performed for equal-mass Fermi gases with up to $N=30$ atoms. We focus on the odd-even oscillations of the ground-state energy of the equal-mass unitary system having up to $N=30$ particles, which are related to the excitation gap of the system. Furthermore, we present a detailed analysis of the structural properties of these systems.

arXiv:0801.2747•Physical Review A•Apr 16, 2008•84 citations

When does a Fermi puddle become a Fermi sea? Emergence of pairing in two-dimensional trapped mesoscopic Fermi gases

Emma Laird, Brendan C. Mulkerin, Jia Wang +1 more

Pairing lies at the heart of superfluidity in fermionic systems. Motivated by recent experiments in mesoscopic Fermi gases, we study up to six fermionic atoms with equal masses and equal populations in two different spin states, confined in a quasi-two-dimensional harmonic trap. We couple a stochastic variational approach with the use of an explicitly correlated Gaussian basis set, which enables us to obtain highly accurate energies and structural properties. Utilising two-dimensional two-body scattering theory with a finite-range Gaussian interaction potential, we tune the effective range to model realistic quasi-two-dimensional scattering. We calculate the excitation spectrum, pair correlation function, and number of pairs as a function of increasing attractive interaction strength. For up to six fermions in the ground state, we find that opposite spin and momentum pairing is maximised well below the Fermi surface in momentum space. By contrast, corresponding experiments on twelve fermions have found that pairing is maximal at the Fermi surface and strongly suppressed beneath [M. Holten et al., Nature 606, 287-291 (2022)]. This suggests that the Fermi sea — which acts to suppress pairing at low momenta via Pauli blocking — emerges in the transition from six to twelve particles.

arXiv:2408.17015•SciPost Physics•Dec 11, 2024•1 citation

Coupled-pair approach for strongly interacting trapped fermionic atoms

C. J. Bradly, Brendan C. Mulkerin, A. Martin +1 more

We present a coupled-pair approach for studying few-body physics in harmonically trapped ultracold gases. The method is applied to a two-component Fermi system of $N$ particles. A stochastically variational Gaussian expansion method is applied, focusing on optimization of the two-body correlations present in the strongly interacting, or unitary, limit. The ground-state energy of the four-, six- and eight-body problem with equal spin populations is calculated with high accuracy and minimal computational effort. We also calculate the structural properties of these systems and discuss their implication for the many-body ultracold gas and other few-body calculations.

arXiv:1405.3744•Physical Review A•Aug 19, 2014•7 citations

Trapped unitary two-component Fermi gases with up to ten particles

Xiangyu Yin, D. Blume

The properties of two-component Fermi gases with zero-range interactions are universal. We use an explicitly correlated Gaussian basis set expansion approach to investigate small equal-mass two-component Fermi gases under spherically symmetric external harmonic confinement. At unitarity, we determine the ground-state energy for systems with up to 10 particles interacting through finite-range two-body potentials for both even and odd numbers of particles. We extrapolate the energies to the zero-range limit using a novel scheme that removes the linear and, in some cases, also the quadratic dependence of the ground-state energies on the two-body range. Our extrapolated zero-range energies are compared with results from the literature. We also calculate the two-body Tan contact and structural properties.

arXiv:1504.04414•Physical Review A•Jul 6, 2015•23 citations

Pairing correlations in a trapped one-dimensional Fermi gas

Stephen Kudla, Dominique Gautreau, Daniel E. Sheehy

We use a BCS-type variational wave function to study attractively interacting quasi-one-dimensional fermionic atomic gases, motivated by cold-atom experiments that access the one-dimensional regime using an anisotropic harmonic trapping potential (with trapping frequencies ${\ensuremath{\omega}}_{x}={\ensuremath{\omega}}_{y}\ensuremath{\gg}{\ensuremath{\omega}}_{z}$) that confines the gas to a cigar-shaped geometry. To handle the presence of the trap along the $z$ direction, we construct our variational wave function from the harmonic oscillator Hermite functions, which are the eigenstates of the single-particle problem. Using an analytic determination of the effective interaction among harmonic oscillator states along with a numerical solution of the resulting variational equations, we make specific experimental predictions for how pairing correlations would be revealed in experimental probes, such as the local density and the momentum correlation function.

arXiv:1404.4081•Physical Review A•Apr 9, 2015•5 citations
We theoretically investigate equal-mass spin-balanced two-component Fermi gases in which pairs of atoms with opposite spins interact via a short-range isotropic model potential. We probe the distinction between two-dimensional and quasi-two-dimensional harmonic confinement by tuning the effective range parameter within two-dimensional scattering theory. Our approach, which yields numerically exact energetic and structural properties, combines a correlated Gaussian basis-set expansion with the stochastic variational method. For systems containing up to six particles, we: 1) Present the ground- and excited-state energy spectra; 2) Study non-local correlations by analysing the one- and two-body density matrices, extracting from these the occupation numbers of natural orbitals, the momentum distributions of atoms and pairs, and the molecular `condensate fraction’; 3) Study local correlations by computing the radial and pair distribution functions. This paper extends current theoretical knowledge on the properties of trapped few-fermion systems as realised in state-of-the-art cold-atom experiments.
We investigate small equal-mass two-component Fermi gases under external spherically symmetric confinement in which atoms with opposite spins interact through a short-range two-body model potential. We employ a non-perturbative microscopic framework, the stochastic variational approach, and determine the system properties as functions of the interspecies s -wave scattering length a s , the orbital angular momentum L of the system, and the numbers N 1 and N 2 of spin-up and spin-down atoms (with N 1 − N 2 = 0 or 1 and N ⩽ 6 , where N = N 1 + N 2 ). At unitarity, we determine the energies of the five- and six-particle systems for various ranges r 0 of the underlying two-body model potential and extrapolate to the zero-range limit. These energies serve as benchmark results that can be used to validate and assess other numerical approaches. We also present structural properties such as the pair distribution function and the radial density. Furthermore, we analyze the one-body and two-body density matrices. A measure for the molecular condensate fraction is proposed and applied. Our calculations show explicitly that the natural orbitals and the momentum distributions of atomic Fermi gases approach those characteristic for a molecular Bose gas if the s -wave scattering length a s , a s > 0 , is sufficiently small.
Using two different numerical methods, we study the behavior of two-component Fermi gases interacting through short-range sss-wave interactions in a harmonic trap. A correlated Gaussian basis-set expansion technique is used to determine the energies and structural properties, i.e., the radial one-body densities and pair distribution functions, for small systems with either even or odd NNN, as functions of the sss-wave scattering length and the mass ratio $\ensuremath{\kappa}$ of the two species. Particular emphasis is put on a discussion of the angular momentum of the system in the BEC-BCS crossover regime. At unitarity, the excitation spectrum of the four-particle system with total angular momentum L=0L=0L=0 is calculated as a function of the mass ratio $\ensuremath{\kappa}$. The results are analyzed from a hyperspherical perspective, which offers unique insight into the problem. Additionally, fixed-node diffusion Monte Carlo calculations are performed for equal-mass Fermi gases with up to N=30N=30N=30 atoms. We focus on the odd-even oscillations of the ground-state energy of the equal-mass unitary system having up to N=30N=30N=30 particles, which are related to the excitation gap of the system. Furthermore, we present a detailed analysis of the structural properties of these systems.
Pairing lies at the heart of superfluidity in fermionic systems. Motivated by recent experiments in mesoscopic Fermi gases, we study up to six fermionic atoms with equal masses and equal populations in two different spin states, confined in a quasi-two-dimensional harmonic trap. We couple a stochastic variational approach with the use of an explicitly correlated Gaussian basis set, which enables us to obtain highly accurate energies and structural properties. Utilising two-dimensional two-body scattering theory with a finite-range Gaussian interaction potential, we tune the effective range to model realistic quasi-two-dimensional scattering. We calculate the excitation spectrum, pair correlation function, and number of pairs as a function of increasing attractive interaction strength. For up to six fermions in the ground state, we find that opposite spin and momentum pairing is maximised well below the Fermi surface in momentum space. By contrast, corresponding experiments on twelve fermions have found that pairing is maximal at the Fermi surface and strongly suppressed beneath [M. Holten et al., Nature 606, 287-291 (2022)]. This suggests that the Fermi sea — which acts to suppress pairing at low momenta via Pauli blocking — emerges in the transition from six to twelve particles.
We present a coupled-pair approach for studying few-body physics in harmonically trapped ultracold gases. The method is applied to a two-component Fermi system of NNN particles. A stochastically variational Gaussian expansion method is applied, focusing on optimization of the two-body correlations present in the strongly interacting, or unitary, limit. The ground-state energy of the four-, six- and eight-body problem with equal spin populations is calculated with high accuracy and minimal computational effort. We also calculate the structural properties of these systems and discuss their implication for the many-body ultracold gas and other few-body calculations.
The properties of two-component Fermi gases with zero-range interactions are universal. We use an explicitly correlated Gaussian basis set expansion approach to investigate small equal-mass two-component Fermi gases under spherically symmetric external harmonic confinement. At unitarity, we determine the ground-state energy for systems with up to 10 particles interacting through finite-range two-body potentials for both even and odd numbers of particles. We extrapolate the energies to the zero-range limit using a novel scheme that removes the linear and, in some cases, also the quadratic dependence of the ground-state energies on the two-body range. Our extrapolated zero-range energies are compared with results from the literature. We also calculate the two-body Tan contact and structural properties.
We use a BCS-type variational wave function to study attractively interacting quasi-one-dimensional fermionic atomic gases, motivated by cold-atom experiments that access the one-dimensional regime using an anisotropic harmonic trapping potential (with trapping frequencies ${\ensuremath{\omega}}_{x}={\ensuremath{\omega}}_{y}\ensuremath{\gg}{\ensuremath{\omega}}_{z}$) that confines the gas to a cigar-shaped geometry. To handle the presence of the trap along the zzz direction, we construct our variational wave function from the harmonic oscillator Hermite functions, which are the eigenstates of the single-particle problem. Using an analytic determination of the effective interaction among harmonic oscillator states along with a numerical solution of the resulting variational equations, we make specific experimental predictions for how pairing correlations would be revealed in experimental probes, such as the local density and the momentum correlation function.