---
title: Spin-J Fermi Gas
url: https://www.emergentmind.com/topics/spin-j-fermi-gas
type: topic
---

# Spin-J Fermi Gas

A spin-\(J\) Fermi gas is a fermionic many-body system in which each particle carries an internal spin degree of freedom in \(\mathbb{C}^{2J+1}\), so the \(N\)-body wave function takes values in \(\wedge^N L^2([0,L];\mathbb{C}^{2J+1})\) and is antisymmetric under exchange of the combined space-spin variables. In the literature represented here, closely related formulations also appear as large-spin or spin-\(F\) degenerate Fermi gases with multiple hyperfine Zeeman sublevels \(m\). The presence of internal spin changes both the allowed two-body scattering channels and the structure of collective dynamics: in dilute one-dimensional systems it produces interaction corrections governed by even-wave and odd-wave scattering lengths together with an effective quantum spin chain, while in dynamical settings it leads to coherent spin-changing forward scattering, collisional damping, multicomponent hydrodynamics, and spin-sensitive collective modes [2509.00536, 1612.02618, 2209.00700].

## 1. Microscopic definition and symmetry structure

In the one-dimensional dilute problem studied explicitly as a spin-\(J\) Fermi gas, the particles occupy an interval \([0,L]\) at density
\[
\rho=\frac{N}{L},
\]
and interact through the Hamiltonian
\[
H=-\sum_{i=1}^{N}\partial_i^2+\sum_{i<j} v(x_i-x_j),
\]
with a repulsive, even, compactly supported two-body potential \(v\) [2509.00536]. The internal Hilbert space \(\mathbb{C}^{2J+1}\) is the defining spin structure, and antisymmetry is imposed on the full space-spin coordinates rather than on the spatial part alone.

A complementary large-spin formulation appears in a one-dimensional harmonically trapped degenerate Fermi gas, where the field operators \(\hat\psi_m(x)\) resolve the occupied Zeeman sublevels \(m\) of a total spin \(F\), and the Hamiltonian is
\[
\hat H=\int dx\sum_m \hat\psi_m^\dagger(x)\left[\frac{\hbar^2\nabla^2}{2M}+\frac12M\omega^2 x^2+Qm^2\right]\hat\psi_m(x) +\frac12\int dx\sum_{klmn}U_{klmn}\hat\psi_k^\dagger\hat\psi_m^\dagger\hat\psi_n\hat\psi_l.
\]
The interaction matrix elements are decomposed into total-spin scattering channels,
\[
U_{klmn}=\sum_{S,M}g_S\langle km|SM\rangle\langle SM|ln\rangle,\qquad g_S = 2\hbar\omega_\perp a_S,
\]
so the multicomponent structure is explicit already at the two-body level [1612.02618].

The opposite extreme is the fully spin-polarized or spinless limit \(q=1\), realized in three dimensions by \(N\) indistinguishable fermions in a periodic box \(\Lambda=[-L/2,L/2]^3\) with Hamiltonian
\[
H_N=\sum_{j=1}^N (-\Delta_{x_j})+\sum_{1\le j<k\le N}V(x_j-x_k),
\]
where \(V\ge 0\) is repulsive [2402.17558]. This case is not a multicomponent gas, but it is the natural reference point for understanding what spin degrees of freedom add or remove. In particular, it isolates the effect of exchange antisymmetry when no opposite-spin channel is available.

## 2. Scattering channels and the dilute limit

For a one-dimensional spin-\(J\) Fermi gas, the central two-body distinction is between even-wave spatial symmetry and odd-wave spatial antisymmetry. Because the total fermionic wave function must be antisymmetric under exchange, the spin symmetry of a pair fixes which spatial channel contributes: antisymmetric spin implies a symmetric spatial part and therefore the even-wave scattering length \(a_e\), while symmetric spin implies an antisymmetric spatial part and therefore the odd-wave scattering length \(a_o\) [2509.00536]. The associated two-body energies on a finite interval are
\[
\frac{4}{R-a_e}, \qquad \frac{4}{R-a_o}.
\]

The dilute regime is the low-density regime in which the interaction range and scattering lengths are small compared with the interparticle spacing. In the one-dimensional spin-\(J\) problem, the theorem assumes
\[
\rho R_0 \ll 1,\qquad \rho |a_e|\ll 1,\qquad \rho a_o \ll 1,
\]
with \(R_0\) the interaction range [2509.00536]. In this regime the leading energy remains the free Fermi-gas energy, while the first correction is controlled by two-body scattering lengths.

A common misconception is that the first interaction correction in a dilute Fermi gas is always \(s\)-wave. The fully spin-polarized three-dimensional case shows otherwise. There, opposite-spin \(s\)-wave scattering is absent, the Pauli principle suppresses same-spin contact, and the leading interaction effect is governed by the \(p\)-wave channel instead [2402.17558]. For \(V\in L^1\) nonnegative, radial, and compactly supported, and for \(ak_F\) sufficiently small and \(N\) large enough, the ground-state energy satisfies
\[
\frac{E_N}{N}\ge \frac{3}{5}k_F^2+\frac{3}{5\pi}a\,k_F^3 +O\!\left((a k_F)^{3+3/10}\,|\log(a k_F)|\right)+O(N^{-1/3}),
\]
with
\[
k_F=(6\pi^2\rho)^{1/3}.
\]
The \(p\)-wave scattering length \(a\) is defined through the zero-energy scattering equation
\[
x\Delta \varphi_0+2\nabla \varphi_0+xV(1-\varphi_0)=0,
\]
with \(\varphi_0(x)\to 0\) as \(|x|\to\infty\), and outside the support of \(V\),
\[
\varphi_0(x)=\frac{a}{|x|}.
\]
The identity
\[
12\pi a=\int_{\mathbb{R}^3}V(x)\,|x|\,(1-\varphi_0(x))\,dx
\]
is the link between scattering theory and the energy correction [2402.17558].

The contrast is sharp: in spinful systems with \(q\ge 2\), opposite-spin particles can interact at leading order via the \(s\)-wave scattering length, whereas in the fully spin-polarized case the leading correction is much smaller and \(p\)-wave dominated [2402.17558]. This difference is one of the main symmetry-based distinctions between single-component and spin-\(J\) gases.

## 3. Ground-state energy and effective spin chains in one dimension

The main asymptotic result for the one-dimensional dilute spin-\(J\) Fermi gas is
\[
E_{J}(N,L)= N\frac{\pi^2}{3}\rho^2\left(1+2\rho \left[a_e+(a_o-a_e)\epsilon^{\mathrm{LS}_J}\right] +\mathcal{O}\!\left((\rho \max(R_0,a_o,a_e))^{6/5}+N^{-2/3}\right)\right).
\]
Here \(N\frac{\pi^2}{3}\rho^2\) is the free one-dimensional Fermi-gas ground-state energy, and the correction of order \(\rho^3\) is encoded by
\[
a_e+(a_o-a_e)\epsilon^{\mathrm{LS}_J}.
\]
The theorem is rigorous for scalar potentials, while the upper bound is also extended to matrix-valued interactions [2509.00536].

The quantity \(\epsilon^{\mathrm{LS}_J}\) is the ground-state energy per site of the effective Lai–Sutherland spin chain
\[
H_{\mathrm{LS}}=\sum_{i=1}^{N} P_S^{i,i+1},
\]
where \(P_S^{i,i+1}\) projects onto the symmetric spin subspace of particles \(i\) and \(i+1\) [2509.00536]. Thus the dilute Fermi-gas energy depends on spin through the expectation of the local projection \(P_S\): the gas lowers its energy by arranging neighboring spins so as to favor the channel with the smaller scattering length.

For spin-\(\tfrac12\), the effective chain becomes the antiferromagnetic Heisenberg chain. Using
\[
2P_S^{i,i+1}=\frac{3}{2}+2\,S_i\cdot S_{i+1},
\]
the ground-state energy per site is
\[
\epsilon^{\mathrm{LS}_{1/2}}=1-\ln 2,
\]
and the explicit asymptotic formula becomes
\[
E(N,L)=N\frac{\pi^2}{3}\rho^2\left(1+2\rho\ln(2)\,a_e+2\rho(1-\ln(2))\,a_o +\mathcal{O}\!\left((\rho \max(R_0,a_o,a_e))^{6/5}+N^{-2/3}\right)\right)
\]
[2509.00536].

A different strong-coupling route leads to another effective spin chain. In a harmonically trapped one-dimensional SU(2) Fermi gas with strong repulsion \(g\to\infty\), the many-body wave function separates into a fixed antisymmetric orbital Slater determinant and spin-sector amplitudes \(a_P\), and to first order in \(1/g\) these amplitudes are governed by
\[
H_s = \left(E_F - \sum_i^{N-1} J_i\right)\mathbb{1} + \sum_{i=1}^{N-1} J_i P_{i,i+1},
\]
with
\[
J_i = \frac{1}{g}\int dx_1\cdots dx_N\, \delta(x_i-x_{i+1})\,\theta(x_1<\cdots<x_N)\left|\frac{\partial \Psi_A}{\partial x_i}\right|^2
\]
[2107.12075]. The sites of this chain are the ordered particle indices rather than fixed lattice positions. This suggests that, in one dimension, both the dilute limit and the strong-repulsion limit reorganize the spinful gas into an effective quantum magnet, although the precise exchange couplings and symmetry content differ.

## 4. Spin dynamics, relaxation, and transport

The nonequilibrium dynamics of a large-spin one-dimensional degenerate Fermi gas can be formulated through the Wigner function
\[
W_{mn}(x,p)=\int dy\,\langle \hat\psi_m^\dagger(x+y/2)\hat\psi_n(x-y/2)\rangle,
\]
which obeys the quantum Boltzmann equation
\[
\frac{d}{dt} W(x,p)+\partial_0 W(x,p)+\frac{i}{\hbar}\left[Q S_z^2+V^\text{mf}(x), W(x,p)\right] -\frac12\left\{\partial_x V^\text{mf}(x),\partial_p W(x,p)\right\}=C[W(x,p)],
\]
with
\[
\partial_0 = \frac{p}{M}\partial_x - M\omega^2 x\,\partial_p
\]
[1612.02618]. In this formulation, the commutator term produces coherent spin dynamics from the quadratic Zeeman term and the mean field, while the collision integral \(C[W]\) produces damping and redistribution.

The coherent process is spin-changing forward scattering, represented by the mean-field potential
\[
V^\text{mf}_{mn}(x)=\int dp \sum_{kl} U_{klnm} W_{kl}(x,p),
\]
or equivalently
\[
V^\text{mf}_{ij}(x)=2\sum_{kl}\int dp\, U_{ijlk} W_{kl}(x,p).
\]
The paper states explicitly that mean-field effects scale linearly with density, \(\omega_\text{mf}\sim ng\), whereas in one dimension the collision rate has sub-linear growth with density and therefore grows slower with increasing density than mean-field interactions [1612.02618]. Because of this, dissipative effects become relatively more important at lower density in one dimension, opposite to the usual three-dimensional intuition.

At low temperature, Pauli blocking is essential. In the large-spin kinetic theory it enters through the shielding factor
\[
S_{ijkl}(x,p,q)=\left(\delta_{ij}-W_{ij}(x,p)\right)\left(\delta_{kl}-W_{kl}(x,q)\right)-W_{ij}(x,p)W_{kl}(x,q),
\]
and in the homogeneous two-component normal gas it appears through the Fermi-liquid transport coefficients [1612.02618, 1012.1607]. For a homogeneous two-component gas in the normal phase, the spin current obeys Fick’s law
\[
\mathbf j=-D\nabla(n_\uparrow-n_\downarrow),
\]
and the diffusion coefficient has the high-temperature asymptotics
\[
D=\frac{9\pi^{3/2}\hbar}{32\sqrt 2\,m}\times 
\begin{cases}
\displaystyle \frac{1}{(k_Fa)^2}\left(\frac{T}{T_F}\right)^{1/2}, & T\ll T_a,\\[1.2ex]
\displaystyle \left(\frac{T}{T_F}\right)^{3/2}, & T\gg T_a,
\end{cases}
\]
while at low temperature Landau Fermi-liquid theory gives
\[
D = \frac{24\hbar}{\pi^3 m} \frac{1+F_0^a}{C_1^2+C_3^2/2-C_1C_3} \left(\frac{T_F}{T}\right)^2.
\]
By connecting the low- and high-temperature asymptotics, the minimum value is estimated as
\[
D_{\min}\sim \frac{9\pi^{3/2}\hbar}{32\sqrt2\,m}\approx 1.1\,\frac{\hbar}{m}
\]
[1012.1607].

In an elongated harmonic trap, spin drag is governed by the transport equation
\[
\gamma\,\bar n_0(z)\alpha_0(z)+k_B T\,\partial_z\!\left(G(\bar n_0(z))\,\partial_z\alpha_0(z)\right)=0,
\]
where \(G\) is the spin conductance and \(\gamma\) the spin damping rate [1307.6395]. In the collisionless limit,
\[
G \simeq k\,\frac{n_0(0)R_{\rm th}^2}{m\,\Gamma_0},
\]
with \(k=15.87\) for a Maxwellian gas, while in the hydrodynamic limit
\[
G \simeq \frac{2\pi n_0(0)R_{\rm th}^2}{m\,\Gamma_0}\ln\!\left(\frac{\Gamma_0}{\omega_\perp}\right)
\]
[1307.6395]. In polarized gases, spin and heat transport also mix through spin Seebeck and spin Peltier effects, and in the classical Boltzmann solution the Seebeck coefficient changes sign around
\[
\lambda/|a| \simeq 3.62
\]
in the \(L=1\) approximation [1208.6580].

## 5. Collective modes, susceptibility, and magnetic response

Collective spin motion provides a direct probe of how interactions suppress or enhance spin degrees of freedom. Near unitarity, the trapped two-component Fermi gas has a spin-dipole mode generated by the out-of-phase operator \(\sum_i z_{i,\uparrow}-\sum_i z_{i,\downarrow}\), and the sum-rule estimate is
\[
\omega_{\rm SD}^2\le \frac{N}{m\int d^3\bm{r}\, z^2\,\chi(r,T)},
\]
where
\[
\chi(r,T)=\lim_{h\rightarrow 0}\frac{n_{\uparrow}(r,T)-n_{\downarrow}(r,T)}{h}
\]
is the local spin susceptibility [1909.11425]. In the noninteracting limit,
\[
\omega_{\rm SD}=\omega_{\rm tr},
\]
but in the superfluid state the mode is enhanced because spin-singlet Cooper pairing suppresses the spin response, and in strongly interacting gases this enhancement occurs even above the superfluid transition temperature due to strong pairing correlations [1909.11425].

A different susceptibility regime appears in a nearly ferromagnetic Fermi gas with positive scattering length. There the relevant sum-rule formula in local density approximation is
\[
\omega_{SD}^2= \frac{N}{m\int d\mathbf r\, z^2\,\chi(n)},
\]
so as the magnetic susceptibility increases, the spin-dipole frequency decreases [1007.4504]. The susceptibility is
\[
\chi=\left(V\,\frac{\partial^2 A}{\partial (N_\uparrow-N_\downarrow)^2}\right)^{-1},
\]
and in Landau Fermi-liquid theory
\[
\frac{\chi_0}{\chi}=(1+F_0^a)\frac{m}{m^*}.
\]
Using ab initio Monte Carlo input, the paper reports a ferromagnetic instability around
\[
k_F a \simeq 0.83
\]
[1007.4504]. Thermal spin fluctuations satisfy
\[
\frac{\Delta (N_\uparrow-N_\downarrow)^2}{N} = k_B T\,\frac{\chi(T)}{n},
\]
whereas quantum fluctuations scale as
\[
\frac{\Delta (N_\uparrow-N_\downarrow)^2}{N} = 2\alpha \left(\frac{12}{\pi^4 N}\right)^{1/3}\ln(C N^{1/3})
\]
[1007.4504].

The contrast between these two settings is instructive. In the paired unitary gas, suppression of \(\chi\) enhances \(\omega_{\rm SD}\); near ferromagnetism, enhancement of \(\chi\) softens \(\omega_{\rm SD}\) [1909.11425, 1007.4504]. The spin-dipole mode is therefore not fixed by the trap alone; it is controlled by the microscopic channel in which interactions reorganize the spin sector.

The dynamical spin-flip response gives a further diagnostic of magnetic order. Within random-phase approximation for a repulsively interacting two-component gas, the transverse susceptibility is
\[
\chi_{+-}(\omega,\mathbf{q})= \frac{\Pi_{+-}(\omega,\mathbf{q})}{1+g\,\Pi_{+-}(\omega,\mathbf{q})},
\]
and in the ferromagnetic phase the response contains both a Stoner particle-hole continuum and a collective magnon mode defined by
\[
1+g\,\Pi_{+-}(\Omega_q,\mathbf{q})=0
\]
[1105.5997]. At small \(q\), the magnon lies below the continuum and is undamped at \(T=0\); in a trap it survives as a pronounced resonance in the spin-flip spectrum [1105.5997].

## 6. Pairing, synthetic gauge fields, and metastable branches

Spinful Fermi gases also support nonequilibrium transitions and metastable branches that have no analogue in a spinless gas. In a one-dimensional attractive Fermi gas on a ring with a time-dependent synthetic gauge flux acting on the spin sector,
\[
\Phi(t)=Ft,
\]
a spin-depairing transition occurs: the system tunnels from a spin-gapped paired state to an unpaired state carrying spin current [1110.1927]. After analytic continuation to a complex twist \(\Phi=i\Psi\), the critical imaginary momentum is
\[
b_{\mathrm{cr}}=\sinh^{-1}|u|,
\]
and the Landau-Dykhne transition probability is
\[
P=e^{-\pi F_{\mathrm{th}}/F},
\]
with threshold field
\[
F_{\mathrm{th}}=\frac{2}{\pi}\int_0^{b_{\mathrm{cr}}} \Delta_{\text{spin}}(b)\,\frac{d\Psi}{db}\,db
\]
[1110.1927]. The paired phase becomes more robust at lower filling, because \(\Psi_{\mathrm{cr}}\) increases as filling decreases.

Another attractive one-dimensional branch is the fermionic super-Tonks-Girardeau state of the spin-\(\tfrac12\) interacting Fermi gas. The homogeneous Gaudin–Yang Hamiltonian is
\[
H=-\frac{\hbar^2}{2m}\sum_{i=1}^N \frac{\partial^2}{\partial x_i^2} +g_{1d}\sum_{i<j}\delta(x_i-x_j),
\]
with
\[
g_{1d}=\frac{\hbar^2 c}{m}=-\frac{2\hbar^2}{m a_{1d}}, \qquad \gamma=\frac{c}{\rho}.
\]
For \(c<0\), the true ground state contains bound singlet pairs, but the fermionic super-Tonks-Girardeau state is the lowest real solution of the Bethe equations: a highly excited, gas-like, metastable state composed of unpaired fermions rather than bound pairs [1005.0461]. It can be prepared by a sudden quench from strong repulsion to strong attraction, and the overlap with the metastable branch becomes close to 1 in the strong-coupling regime. The paper estimates
\[
P(N,M,\gamma)\approx 1-\frac{a(N,M)}{\gamma^2},\qquad a(N,M)\le 64\pi^2 N M^2,
\]
and identifies the lowest breathing mode as an experimental signature, with a maximum
\[
\omega^2/\omega_x^2 \approx 4.3
\]
[1005.0461].

Spin imbalance produces a different instability. In a weakly attractive spin-imbalanced Fermi gas, the most favorable correlated cluster is not necessarily a Cooper pair: the optimal composition satisfies
\[
\frac{N_\uparrow}{N_\downarrow}=\frac{\nu_\uparrow}{\nu_\downarrow},
\]
so the particle-number ratio in the instability matches the ratio of the densities of states at the two Fermi surfaces [1712.09847]. In the balanced limit, this reduces to the ordinary Cooper problem.

Synthetic spin-orbit coupling provides a spectroscopic route to the dressed spin structure. In a Raman-coupled degenerate \(^{6}\)Li Fermi gas, spin-injection spectroscopy reconstructs both the energy-momentum dispersion and the spin composition of the dressed bands. The single-particle Hamiltonian can be written as
\[
\mathcal{H} = \frac{\hbar^2 \mathbf{k}^2}{2m} + \frac{\hbar^2 Q}{2m}\sigma_z q + \frac{\hbar \Omega_R}{2}\sigma_x + \frac{\delta}{2}\sigma_z + \frac{E_R}{4},
\]
revealing the spin-orbit gap opened by Raman coupling [1205.3483]. For energies within this gap, the system acts as a spin diode, and adding a radiofrequency coupling produces a spin-orbit coupled lattice with multiple gaps and spinful band structure [1205.3483].

Taken together, these results show that the term “spin-\(J\) Fermi gas” denotes more than a simple increase in internal degeneracy. It identifies a class of fermionic systems in which exchange antisymmetry, channel-resolved scattering, and multicomponent spin dynamics reorganize the many-body problem into effective spin chains, nested quasiparticle hydrodynamics, spin-sensitive collective modes, and nonequilibrium transitions whose form depends decisively on dimensionality, polarization, and pairing channel.

Source: https://www.emergentmind.com/topics/spin-j-fermi-gas