Rf spectra and pseudogap in ultracold Fermi gases across the BCS-BEC crossover from pairing fluctuation theory

By advancing pairing fluctuation theory with particle-hole fluctuations and full numerical convolution to accurately model spectral broadening and pair-hole scattering, this study quantitatively demonstrates that the pseudogap in ultracold Fermi gases emerges continuously across the BCS-BEC crossover due to pairing fluctuations, achieving strong agreement with experimental data.

The Gist

Imagine you have a crowded dance floor filled with people (the fermions). In a normal crowd, everyone dances independently, bumping into each other occasionally but mostly following their own rhythm. This is like a normal gas.

But what happens if you turn up the music and make the dancers incredibly sensitive to each other? They start pairing up, holding hands, and moving in perfect sync. This is superfluidity (or superconductivity in metals), where the whole crowd moves as one giant, coordinated wave.

Now, here is the mystery: What happens just before they fully lock into that perfect dance? Is there a "ghost" of the dance already forming? This ghost is called the pseudogap.

For decades, scientists have argued about what this pseudogap is. Is it a pre-formed dance pair waiting to happen? Or is it something else entirely, like a different kind of order?

This paper, by Li, Sun, and their team, acts like a high-tech "dance floor simulator" to settle the argument. They used a supercomputer to model ultracold atoms (which act like the dancers) and found that the pseudogap is indeed caused by pairing fluctuations—essentially, the dancers are trying to hold hands and let go, over and over again, even before the music gets loud enough for a full-blown dance party.

Here is a breakdown of their discovery using simple analogies:

1. The Problem: The "Blurry" Dance Floor

In the past, scientists tried to model this using a "perfect crystal" approach. They assumed that if pairs formed, they were solid and unchanging.

• The Analogy: Imagine taking a photo of the dancers. The old models assumed the dancers were frozen in perfect poses.
• The Reality: In the real world, dancers are jittery. They form a pair, spin, break apart, and find a new partner. This "jitter" (or fluctuation) blurs the picture. The old models ignored this blur, so their predictions didn't match real experiments.

2. The New Tool: The "Convolution" Camera

The authors developed a new mathematical method they call a "full numerical convolution."

• The Analogy: Instead of taking a frozen photo, they used a high-speed video camera that captures every wobble, every near-miss, and every momentary hand-hold.
• What it does: This method accounts for two critical things the old models missed:
  1. Finite Lifetime: Pairs don't last forever; they have a "half-life." They exist for a split second and then vanish. This makes the energy levels "fuzzy" rather than sharp.
  2. The "Hartree" Shift: When dancers crowd together, they push against each other, shifting the whole group's position. The authors calculated this "push" (the Hartree energy) accurately for the first time in this context.

3. The Journey: From Solo Dancers to a Mosh Pit

The paper studies the BCS-BEC crossover. Think of this as a slider on a mixing board:

• BCS Side (Weak Interaction): The dancers are far apart. They only pair up when the music is very specific and slow. They are like shy couples holding hands across a room.
• BEC Side (Strong Interaction): The dancers are glued together. They form tight, permanent couples (molecules) even before the music starts.
• The Unitary Point (The Middle): This is the "sweet spot" where the interaction is strongest. It's chaotic, energetic, and the hardest to predict.

The authors simulated the entire journey from the shy couples to the glued-together mosh pit.

4. The Big Discovery: The "Ghost" is Real

Their simulation showed that as you move from the shy-couple side to the glued-together side:

• The Pseudogap appears naturally: Even before the full superfluid state forms, the "ghost pairs" (the jittery, short-lived hand-holds) create a gap in the energy spectrum.
• It matches reality: When they compared their simulation to real experiments with Lithium-6 atoms (the "dancers" in a lab), the numbers matched perfectly.
• The Pair Lifetime: They found that above a certain energy level (twice the gap size), the pairs stop being distinct dancers and become a "diffusive soup." They are no longer holding hands; they are just bumping into each other in a chaotic cloud.

5. Why This Matters

This isn't just about cold atoms in a lab. The same physics applies to high-temperature superconductors (the materials that could one day power lossless power grids or levitating trains).

• The Analogy: If you understand how the "ghost pairs" form in a simple, controllable system (like the cold atoms), you can finally understand why those mysterious superconductors behave the way they do.
• The Conclusion: The paper confirms that the pseudogap is not a mysterious, competing force. It is simply the result of particles trying to pair up, failing, and trying again. It's the sound of the crowd trying to dance before the beat drops.

Summary

The authors built a sophisticated "virtual dance floor" that accounts for the jittery, short-lived nature of particle pairs. By doing so, they proved that the mysterious "pseudogap" is just the natural result of particles pairing up and breaking apart. Their model fits the experimental data perfectly, solving a decades-old debate and giving us a clearer map of how matter behaves when it's pushed to its limits.

Technical Summary

1. Problem Statement

The pseudogap phenomenon—a suppression of the density of states (DOS) in the normal state above the critical temperature ($T_c$)—is a hallmark of strongly interacting Fermi systems, observed in both high-temperature superconductors and ultracold atomic gases. However, its precise microscopic origin remains debated:

• The Debate: Is the pseudogap caused by pre-formed Cooper pairs (pairing fluctuations) or by competing orders (e.g., d-density waves, loop currents)?
• The Challenge: While ultracold Fermi gases offer a tunable platform to study the BCS-BEC crossover via Feshbach resonances, existing theoretical models often rely on analytic approximations (such as the standard pseudogap approximation) that fail to capture spectral broadening (finite quasiparticle lifetime) and Hartree energy shifts accurately. Consequently, previous theories have struggled to achieve quantitative agreement with recent high-precision momentum-resolved radio-frequency (rf) spectroscopy experiments.

2. Methodology

The authors employ and advance a Pairing Fluctuation Theory to calculate spectral functions and rf spectra across the entire BCS-BEC crossover. Key methodological advancements include:

• Incorporation of Particle-Hole Fluctuations:
  • The theory extends beyond the standard particle-particle channel by including particle-hole fluctuations (induced interactions).
  • These fluctuations renormalize the effective interaction in the particle-particle channel, effectively screening the bare pairing interaction. This is crucial for accurately shifting the $T_c$ curve and determining the effective pairing strength.
• Iterative Numerical Convolution Framework:
  • Moving beyond analytic pseudogap approximations (which assume delta-function spectral weights), the authors develop an iterative framework to compute the self-energy ($\Sigma$) and pair susceptibility ($\chi$) via full numerical convolutions in the real-frequency domain.
  • Process:
    1. Start with a BCS-like Green's function ($G$) as an initial guess.
    2. Compute the retarded pair susceptibility $\chi^R$ and the imaginary part of the T-matrix $Im t^R_{pg}$.
    3. Calculate the retarded self-energy $\Sigma^R$ using convolution integrals.
    4. Update the Green's function $G$ and the spectral function $A(k, \omega)$.
    5. Iterate until convergence (or near-convergence) to capture lifetime effects.
• Key Physical Quantities Calculated:
  • Spectral Function $A(k, \omega)$: Includes quasiparticle broadening due to finite pair lifetimes.
  • Pair Spectral Function $B(q, \Omega)$: Describes the excitation spectrum of finite-momentum pairs.
  • Hartree Energy ($\bar{E}_{Hartree}$): Calculated as the real part of the self-energy at the Fermi surface, representing a significant energy shift previously neglected in simple approximations.
  • Chemical Potential ($\mu$): Determined self-consistently with the number equation.

3. Key Contributions

1. Quantitative Agreement with Experiment: The theory achieves quantitative agreement with recent momentum-resolved microwave spectroscopy data on homogeneous $^6$Li gases (specifically Ref. [26]) across the BCS-BEC crossover, including the unitary limit.
2. Resolution of Spectral Broadening: The framework explicitly captures the spectral broadening of fermions caused by finite pair lifetimes, a feature missing in previous analytic models. This allows for the accurate extraction of quasiparticle dispersion from Energy Distribution Curves (EDCs).
3. Identification of Pair-Hole Scattering: The inclusion of particle-hole fluctuations reveals a substantial Hartree energy contribution, which acts as a shift in the chemical potential and is essential for matching experimental observations.
4. Characterization of Pair Dynamics: The study elucidates the transition of pairs from long-lived bound states to diffusive, short-lived excitations as energy increases.

4. Key Results

• Evolution of the Pseudogap:
  • The pseudogap ($\Delta$) emerges continuously as the system moves from the BCS regime toward the BEC regime.
  • The gap grows monotonically with interaction strength, evidenced by the DOS and EDCs.
  • In the BEC regime, the pseudogap persists well above $T_c$ (up to $T \approx 1.8 T_c$), whereas in the near-BCS regime, it vanishes closer to $T_c$.
• Spectral Function Characteristics:
  • BCS Regime: Shows clear particle and hole branches with back-bending behavior near the Fermi surface.
  • Unitary & BEC Regimes: The dispersion becomes S-shaped or hybridized. The spectral weight spreads significantly, and the "back-bending" point shifts, reflecting strong particle-hole mixing and a reduced effective Fermi surface.
  • Broadening: The spectral peaks exhibit significant broadening, particularly in the unitary regime, indicating short quasiparticle lifetimes.
• Pair Spectral Function ($B(q, \Omega)$):
  • At low energies ($\Omega < 2\Delta$), pairs exhibit well-defined, long-lived parabolic dispersions.
  • Diffusive Transition: Above an energy threshold of roughly $2\Delta$, the pair spectral function broadens rapidly. Pairs become diffusive and short-lived, governed by virtual binding and unbinding processes. This suggests $2\Delta$ acts as the pair-breaking energy scale.
• Thermodynamic Properties:
  • The calculated Bertsch parameter ($\xi = \mu/E_F$ at $T=0$) at unitarity is 0.364, consistent with experimental and Monte Carlo results.
  • The physical chemical potential $\mu$ shows a non-monotonic temperature dependence, peaking near $T_c$, consistent with mean-field BCS behavior but modified by strong fluctuations.

5. Significance

• Validation of Pairing Origin: The quantitative match between the calculated rf spectra/EDCs and experimental data provides strong evidence that the pseudogap in strongly interacting Fermi gases is primarily driven by pairing fluctuations (pre-formed pairs) rather than competing orders.
• Microscopic Understanding: By moving beyond the "pseudogap approximation" (which treats pairs as static), this work offers a more rigorous microscopic description of how finite pair lifetimes and Hartree shifts shape the spectral properties of the system.
• Implications for High-$T_c$ Superconductors: Since the pseudogap mechanism in ultracold gases is now quantitatively understood via pairing fluctuations, the authors argue this picture is likely applicable to high-temperature superconductors, suggesting that pre-formed pairs play a central role in the normal state of cuprates.
• Theoretical Benchmark: The developed iterative numerical framework sets a new standard for calculating spectral functions in strongly correlated Fermi systems, overcoming the limitations of previous analytic approximations.

In conclusion, this paper successfully bridges the gap between theoretical pairing fluctuation models and high-precision experimental data, establishing a robust framework for understanding the pseudogap phenomenon across the BCS-BEC crossover.