Hartree shift and pairing gap in ultracold Fermi gases in the framework of low-momentum interactions

This paper employs a diagrammatic Bogoliubov many-body perturbation theory with low-momentum interactions to calculate the Hartree shift and pairing gap in ultracold Fermi gases, successfully reproducing established weak-coupling corrections and showing reasonable agreement with experimental and quantum Monte-Carlo results near the unitary regime.

The Gist

The Big Picture: A Dance Floor of Atoms

Imagine a giant, super-cold dance floor filled with two types of dancers: Spin-Up and Spin-Down. In the world of physics, these are fermions (like electrons or Lithium-6 atoms).

In a normal gas, these dancers just bump into each other randomly. But if you cool them down enough and tune their interactions just right, they start to pair up. They hold hands and dance in perfect sync. This is called superfluidity (or superconductivity in metals).

This paper is about figuring out exactly how tightly these dancers hold hands (the "pairing gap") and how much the crowd pushes or pulls on them (the "Hartree shift").

The Problem: The "Perfect" Theory vs. Reality

Physicists have a simple theory called BCS theory (named after three scientists) that predicts how these pairs form. It's like a simple map that says, "If you walk this path, you'll get to the destination."

However, in the real world, things are messy. The dancers don't just hold hands with one partner; they are constantly jostled by the crowd.

• The "Hartree Shift": Imagine the crowd pushing the dancers apart or pulling them together. This changes the energy of the dance.
• The "Pairing Gap": This is the energy required to break a pair apart. If the gap is small, the pair breaks easily. If it's large, they are stuck together.

The simple BCS map is often wrong. It overestimates how strong the pairs are. In the 1960s, scientists found that if you account for the crowd's jostling (fluctuations), the pairs are actually 50% weaker than the simple map predicted. This is a famous correction known as the Gor'kov-Melik-Barkhudarov (GMB) correction.

The Authors' Approach: A New Way to Count

The authors of this paper, Michael Urban and S. Ramanan, wanted to see if they could calculate these corrections accurately using a method called Many-Body Perturbation Theory.

Think of this method like trying to calculate the total cost of a party.

1. Level 1 (Mean Field): You just count the cost of the food. (This is the simple BCS theory).
2. Level 2 (Second Order): You add the cost of the drinks and the music.
3. Level 3 (Third Order): You add the cost of the decorations, the cleanup crew, and the tips.

The authors went up to Level 3. They wanted to see if adding these extra layers of complexity would give them the right answer.

The Twist: The "Cutoff" and the "Self-Consistent" Fix

There was a catch. In their math, they had to draw a line in the sand, called a cutoff ($\Lambda$). Imagine they are only allowed to count dancers moving slower than a certain speed.

• If they set the speed limit too high (infinity), the math breaks and gives nonsense results.
• If they set it too low, they miss important dancers.

The Old Way: They used to say, "Let's just pick a speed limit, do the math, and hope the result doesn't change much if we tweak the limit."

• Result: For the "Hartree shift" (the crowd push), this worked okay. But for the "Pairing gap" (the hand-holding strength), the math was unstable. The answer kept changing depending on where they drew the line.

The New Way (The Breakthrough): The authors realized they needed to be Self-Consistent.

• Analogy: Imagine you are trying to guess the temperature of a room.
  • Old Way: You guess the temperature, measure the air, and say, "Okay, my guess was close."
  • New Way: You guess the temperature, measure the air, realize your guess was wrong, update your guess, measure again, and repeat until your guess and the measurement match perfectly.

By forcing their math to "match itself" (making the crowd push and the hand-holding strength consistent with each other), they found that the math finally stabilized.

What Did They Find?

1. In Weak Interactions (The "Cold, Quiet" Dance):
   When the dancers don't interact much, their new method worked perfectly. It reproduced the famous 50% reduction in pairing strength (the GMB correction). The math became stable, and the "cutoff" didn't matter anymore. This is a huge success.

2. In Strong Interactions (The "Unitary" Dance):
   When the dancers interact very strongly (the "Unitary" regime, where they are as chaotic as possible), the math got messy again. Even with their new self-consistent method, the results still depended a bit on where they drew the line.

   • Why? Because at this level of chaos, simple "pairing" isn't enough. You need to account for groups of three or four dancers interacting at once (three-body forces), which their current math doesn't fully capture.

3. Comparison with Experiments:
   They compared their results with real experiments done with Lithium atoms.

   • Good News: Their results were very close to what was observed in the lab.
   • The "Temperature" Glitch: The experiments were done at a tiny, but non-zero, temperature. The authors' math was for absolute zero. They realized that if you account for the slight warmth of the lab, the experimental data matches their theory even better. It's like the dancers were slightly distracted by the heat, making the pairs look weaker than they would be in a perfect freeze.

Why Does This Matter?

This isn't just about cold atoms in a lab.

• Neutron Stars: The inside of a neutron star is essentially a giant gas of neutrons (which are fermions) interacting strongly.
• The Connection: The physics of these cold atoms is very similar to the physics of neutron stars. By perfecting the math for the atoms, the authors are building a better map for understanding the dense, super-dense matter inside stars.

Summary in a Nutshell

The authors built a sophisticated mathematical model to predict how atoms pair up in a super-cold gas.

• They discovered that to get the right answer, you can't just do the math once; you have to loop back and check your work (self-consistency).
• When they did this, they successfully predicted how the "crowd" weakens the atomic pairs in simple situations.
• In the most chaotic situations, the math is still a bit wobbly, suggesting that even more complex interactions (groups of three atoms) are needed to get the perfect answer.
• Overall, their work bridges the gap between simple theories and the messy reality of the universe, helping us understand everything from cold lab experiments to the hearts of neutron stars.

Technical Summary

1. Problem Statement

The paper addresses the theoretical challenge of accurately calculating the pairing gap ($\Delta$) and the Hartree shift (normal self-energy, $U$) in two-component ultracold Fermi gases (e.g., $^6$Li) across the BCS-BEC crossover, specifically on the BCS side.

• Limitations of Mean-Field Theory: Standard BCS mean-field theory (Hartree-Fock-Bogoliubov or HFB) significantly overestimates the pairing gap. It fails to account for medium polarization effects (screening), which are known to suppress the gap by a factor of $\approx 0.45$ in the weak-coupling limit (the famous Gor'kov-Melik-Barkhudarov or GMB correction).
• Limitations of Standard Regularization: Conventional perturbative approaches often rely on taking the momentum cutoff $\Lambda \to \infty$ to regularize contact interactions. This leads to divergent integrals requiring complex resummation of ladder diagrams to recover physical results like the GMB correction.
• Goal: The authors aim to compute these quantities up to third order in perturbation theory using a low-momentum effective interaction with a finite cutoff scaled to the Fermi momentum ($k_F$), and to implement a self-consistent scheme to recover the correct weak-coupling physics without ad-hoc resummations.

2. Methodology

A. Low-Momentum Effective Interaction

Instead of using a standard contact interaction with $\Lambda \to \infty$, the authors employ a separable interaction $V(q, q') = g F(q) F(q')$ with a finite momentum cutoff $\Lambda$.

• Scaling: The cutoff is scaled with the Fermi momentum, $\Lambda \propto k_F$.
• Form Factor: The form factor $F(q)$ is constructed to reproduce the s-wave scattering phase shifts of a contact interaction up to the cutoff $\Lambda$.
• Advantage: This approach naturally includes screening effects and higher-order corrections without needing to explicitly resum infinite ladder diagrams in the vertices, provided the cutoff is finite and scaled correctly.

B. Diagrammatic Bogoliubov Many-Body Perturbation Theory (BMBPT)

The authors use the Nambu-Gor'kov formalism to handle the superfluid state.

• Formalism: They define a two-component field operator $\Psi_k = (a_{k\uparrow}, a^\dagger_{-k\downarrow})^T$.
• Expansion: The self-energy $\Sigma$ is expanded in powers of the interaction.
  • First Order: Includes the Hartree shift and the pairing gap (HFB level).
  • Second & Third Order: Includes diagrams representing particle-hole and particle-particle fluctuations (screening).
• Diagrams: Calculations include all relevant diagrams up to third order (14 diagrams in the standard HFB expansion).

C. Self-Consistent Scheme

A critical innovation of this work is the move from a "naive" perturbative expansion (where the gap is fixed to the HFB value) to a self-consistent scheme.

• Motivation: In naive perturbation theory, the gap reduction is insufficient because the screening is calculated on top of an unphysical, large HFB gap.
• Implementation: The authors impose a condition where the total self-energy vanishes at the Fermi surface ($k_F$) for specific energies:
  • $U_{k_F} = \Sigma_{11}(k_F, 0)$
  • $\Delta_{k_F} = \Sigma_{12}(k_F, E_{k_F})$
• Iterative Solution: They solve the non-linear equation $\Delta = \Delta^{(1)} + \Delta^{(2)} + \Delta^{(3)}$ self-consistently. This effectively renormalizes the propagators used in higher-order diagrams, allowing the medium polarization to feedback into the gap calculation.

3. Key Contributions

1. Third-Order Calculation: The paper provides the first explicit calculation of the pairing gap and Hartree shift up to third order in BMBPT using Nambu-Gor'kov propagators for ultracold gases with low-momentum interactions.
2. Recovery of GMB Limit: By employing the self-consistent scheme, the authors successfully recover the Gor'kov-Melik-Barkhudarov (GMB) suppression factor ($\approx 0.45$) in the weak-coupling limit. This demonstrates that self-consistency is essential to capture the physics of medium polarization without resummation techniques.
3. Cutoff Independence Analysis: The study systematically analyzes the dependence of results on the cutoff ratio $\Lambda/k_F$. They identify a "plateau" region (typically $1.5 \lesssim \Lambda/k_F \lesssim 3$) where results become cutoff-independent, indicating convergence of the perturbative expansion.
4. Comparison with Experiments and QMC: The results are benchmarked against Quantum Monte Carlo (QMC) simulations and experimental data (radiofrequency spectroscopy) across the BCS-BEC crossover.

4. Key Results

• Weak Coupling ($1/(k_F a) \lesssim -2$):

  • The self-consistent third-order results show excellent convergence.
  • The pairing gap is reduced by $\sim 50\%$ compared to HFB, matching the GMB prediction.
  • The Hartree shift converges to the Galitskii result.
  • Results are largely independent of the cutoff within the identified plateau range.

• Intermediate Coupling ($1/(k_F a) \approx -1$):

  • Convergence is less robust; the plateau region shrinks.
  • The third-order correction remains significant, suggesting missing higher-order terms or induced many-body forces.

• Unitary Limit ($1/(k_F a) = 0$):

  • Perturbation theory fails to converge. No plateau is observed in the cutoff dependence for either the gap or the Hartree shift.
  • Despite the lack of convergence, the results (within the large uncertainty bands) are qualitatively compatible with QMC and experimental data.
  • The chemical potential $\mu$ shows surprising stability, as the cutoff dependencies of $\Delta$ and $U$ largely cancel out in the expression for $\mu$.

• Finite Temperature Effects:

  • The authors discuss discrepancies between their $T=0$ results and experimental data (which were taken at finite $T \approx 0.13 T_F$).
  • A rough estimate suggests that finite temperature effects could explain why experimental gaps appear smaller than the theoretical $T=0$ predictions, particularly near the crossover region.

5. Significance and Outlook

• Validation of Low-Momentum Interactions: The work validates the approach of using finite, scaled cutoffs to handle ultracold Fermi gases. It shows that this method can reproduce known many-body physics (like GMB) more naturally than traditional $\Lambda \to \infty$ regularization.
• Neutron Matter Applications: The framework is directly applicable to neutron matter (relevant for neutron star crusts). Since neutron matter operates in a regime similar to the weak-to-intermediate coupling of cold atoms (but never reaches the unitary limit due to effective range effects), the authors anticipate this method will yield reliable predictions for neutron pairing gaps, improving upon previous phenomenological approaches.
• Future Directions:
  • Extending the formalism to finite temperatures to directly compare with experiments.
  • Investigating induced three-body interactions that arise when the cutoff is lowered below $\Lambda \approx 2k_F$. These repulsive three-body forces are hypothesized to be necessary to further reduce the gap and improve convergence in the unitary regime.
  • Applying the method to calculate the equation of state for neutron stars.

In summary, this paper establishes a robust, self-consistent perturbative framework for calculating superfluid properties in Fermi gases, successfully bridging the gap between simple mean-field theory and complex many-body correlations, with direct implications for nuclear astrophysics.