Equal-spin and opposite-spin density-density correlations in the BCS-BEC crossover: Gauge Symmetry, Pauli Exclusion Principle, Wick's Theorem and Experiments

We develop a general theory of spin-dependent density-density correlations, that is valid for any temperature, interactions, dimensions and mass or population status of Fermi gases with two internal states. We use gauge invariance and the Pauli principle to establish constraints on the spin-dependent density-density correlations that are consistent with the fluctuation-dissipation and Wick’s theorem. As an example, we study the spin-dependent density-density correlations from the BCS to the Bose regime in two dimensions at zero temperature, inspired by experiments in 6Li. We show that two-particle irreducible contributions involving collective excitations, many-particle scattering and vertex corrections, are essential to describe experiments. In particular they turn out to be responsible for the emergence of an experimentally observed minimum in the opposite-spin density-density correlations.

💡 Research Summary

The paper presents a comprehensive theoretical framework for spin‑dependent density‑density correlation functions in two‑component Fermi gases that is valid across all temperatures, interaction strengths, dimensionalities, mass ratios, and population imbalances. By invoking global U(1) gauge invariance and the Pauli exclusion principle, the authors derive exact constraints on the correlation functions χσσ′(q,ω). These constraints guarantee consistency with the fluctuation‑dissipation theorem and Wick’s theorem, ensuring that the correlation matrix is real, symmetric, and respects the sign reversal required for identical‑spin pairs.

To illustrate the power of the formalism, the authors focus on a two‑dimensional gas of ^6Li atoms at zero temperature and trace the evolution of the spin‑resolved correlations from the Bardeen‑Cooper‑Schrieffer (BCS) regime to the Bose‑Einstein condensate (BEC) side. In the weak‑coupling BCS limit, a simple mean‑field description captures the formation of Cooper pairs but fails to reproduce the detailed q‑dependence of the correlations measured experimentally. The authors therefore incorporate collective excitations—Goldstone (phase) and Higgs (amplitude) modes—through a functional‑integral approach, and they supplement this with many‑body scattering processes encoded in a two‑particle‑irreducible (2PI) vertex.

The analysis reveals that opposite‑spin correlations χ↑↓ are dominated by the coupling to the Goldstone mode. This coupling generates a pronounced dip (negative minimum) in χ↑↓ at a finite momentum, exactly as observed in recent Bragg‑scattering experiments on 6Li. In contrast, same‑spin correlations χ↑↑ are largely governed by amplitude fluctuations and by vertex corrections arising from repeated particle‑particle scattering. By solving the Bethe‑Salpeter equation with a T‑matrix that includes both ladder and bubble diagrams, the authors demonstrate that these vertex corrections are essential to reproduce the overall magnitude and momentum dependence of χ↑↑.

A key result is that neglecting the 2PI contributions—i.e., using only mean‑field Green’s functions or single‑particle self‑energies—cannot generate the experimentally observed minimum in χ↑↓. Only when the full set of collective, many‑body, and vertex corrections is retained does the theory quantitatively match the data, reproducing both the position and depth of the dip. This finding underscores the necessity of treating the density response beyond the random‑phase approximation in the strongly interacting crossover regime.

The paper also discusses extensions to finite temperature and three dimensions. The same gauge‑symmetry and Pauli‑principle constraints remain operative, but thermal fluctuations introduce additional Matsubara sums and require temperature‑dependent vertex functions. The authors argue that the qualitative role of collective modes and vertex corrections persists, suggesting that the present framework can be adapted to a broad class of ultracold‑atom experiments, including imbalanced mixtures and mass‑asymmetric systems.

In summary, the work establishes a rigorous, symmetry‑based foundation for spin‑resolved density correlations, identifies the indispensable contributions of collective excitations and many‑body vertex corrections, and provides a quantitative explanation for the minimum observed in opposite‑spin correlations across the BCS‑BEC crossover. This advances our theoretical understanding of strongly correlated fermionic superfluids and offers a direct bridge between microscopic many‑body theory and precision measurements in ultracold atomic gases.