Free energy of the gas of spin 1/2 fermions beyond the second order and the Stoner phase transition

By computing the free energy of a dilute spin-1/2 fermion gas up to third order in the scattering length and resumming specific infinite diagram classes, the authors demonstrate that including particle-hole ring contributions eliminates the predicted Stoner ferromagnetic phase transition, suggesting that such a transition in repulsive systems requires the inclusion of higher-order scattering parameters or is relevant primarily for cold atomic gases with artificially enhanced scattering lengths.

The Gist

The Big Picture: A Dance Floor of Tiny Particles

Imagine a giant, invisible dance floor filled with millions of tiny dancers. These dancers are fermions (a type of subatomic particle, like electrons or cold atoms). They have a very strict rule: no two dancers can occupy the exact same spot at the same time. This is known as the "Pauli Exclusion Principle."

These dancers come in two "colors": Red (spin up) and Blue (spin down).

The scientists in this paper wanted to answer a big question: If these dancers start pushing each other away (repelling), will they spontaneously sort themselves out?

• The "Normal" State (Paramagnetic): Red and Blue dancers mix randomly. The floor looks purple from a distance.
• The "Ordered" State (Ferromagnetic): The dancers get so annoyed by the pushing that they decide to segregate. All the Reds go to one side, and all the Blues go to the other. The floor becomes half red, half blue. This is called Itinerant Ferromagnetism.

For decades, physicists thought they knew exactly when this "segregation party" would start. But this paper suggests that previous calculations were missing a crucial piece of the puzzle.

---

The Story of the Calculation

1. The Old Way: Counting Handshakes (Perturbation Theory)

Imagine trying to predict how the dancers behave by counting how many times they bump into each other.

• Level 1: You count the first bump.
• Level 2: You count the second bump.
• Level 3: You count the third bump.

Previous studies calculated up to the "second bump" and found that if the dancers push hard enough, they will segregate. It looked like a clear, sudden switch (a "phase transition").

2. The New Challenge: The Third Bump

The authors of this paper decided to do the math for the third bump (the third order of interaction). This is incredibly hard because the math gets messy, like trying to untangle a knot of headphones while running.

They developed a new, systematic way to do this "knot-untangling" using Imaginary Time.

• Analogy: Imagine the dancers are moving in slow motion through a foggy dream. By calculating their paths in this "dream time," the math becomes much cleaner and easier to solve than in real time.

They successfully calculated the energy of the system when the dancers push each other three times.

3. The Twist: The "Ring" Diagrams

Here is where it gets really interesting. In physics, you don't just look at one-on-one bumps. You have to look at groups of dancers interacting in loops, called Ring Diagrams.

There are two types of rings:

1. Particle-Particle Rings: Dancers bumping into other dancers of the same type (Red hitting Red, Blue hitting Blue).
2. Particle-Hole Rings: A more complex interaction where a dancer moves, leaving an empty spot (a "hole"), and another dancer fills it.

The Surprise:

• Previous studies only looked at the Particle-Particle Rings. When they added these up (resummed them), they found the dancers did segregate. The "Ferromagnetic" state appeared!
• This paper added the Particle-Hole Rings to the mix.

The Result: When they added the second type of ring, the segregation completely disappeared.

It's as if you were predicting a riot based on people shouting at each other (Particle-Particle), but then you realized that the people were also quietly whispering and making peace deals (Particle-Hole). Once you include the whispers, the riot never happens. The dancers stay mixed (purple) no matter how hard they push.

---

Why Does This Matter?

The "Feshbach Resonance" Trick

In real-world experiments with cold atoms (like in a lab in a university), scientists use a magnetic trick called a Feshbach resonance to make the atoms push each other very hard. They do this to try to create this "Ferromagnetic" state.

• The Problem: Experiments have tried to create this state for years and failed. The atoms just don't segregate.
• The Old Explanation: Scientists thought, "Maybe the atoms are forming pairs (dimers) too quickly and ruining the experiment."
• The New Explanation (from this paper): Maybe the atoms never wanted to segregate in the first place! The previous theories that predicted segregation were missing the "Particle-Hole" interactions. When you include them, the math says the mixed state is actually the most stable one.

The "Metastable" Caveat

The authors add a small warning. Their math describes a "metastable" state.

• Analogy: Imagine a ball sitting in a shallow dip on a hill. It looks stable, but if you push it hard enough, it rolls down to a deeper valley (the true ground state, which might be a pair of atoms).
• The "Ferromagnetic" state might be a shallow dip that could exist, but it's so unstable that the atoms might just fall apart into pairs before they can organize.

The Takeaway

1. Math is Hard: Calculating how particles interact gets exponentially harder the more you look at it.
2. Missing Pieces Matter: Ignoring just one type of interaction (the "Particle-Hole" rings) led to a completely wrong prediction about the behavior of the gas.
3. The Verdict: If you have a gas of cold atoms that repel each other, and you ignore the complex "whispers" (Particle-Hole interactions), you might think they will spontaneously sort themselves into teams. But if you listen to the whole conversation, they probably won't. They'll stay mixed.

This paper doesn't just fix a number; it changes the story of what these cold atoms are actually doing, suggesting that the "Ferromagnetic" state might be a ghost that only exists in incomplete math, not in reality.

Technical Summary

1. Problem Statement

The paper addresses the theoretical prediction of the Stoner phase transition (the transition from a paramagnetic to a ferromagnetic state) in a dilute gas of non-relativistic spin-1/2 fermions interacting via a repulsive, spin-independent two-body potential.

• Context: Previous studies using effective field theory (EFT) had computed the ground-state energy density up to the second order in the gas parameter ($k_F a_0$), predicting a first-order phase transition. Later, resummations of specific "particle-particle" ring diagrams suggested a continuous transition.
• Gap: There was a lack of a systematic computation of the free energy $F$ at finite temperatures beyond the second order. Furthermore, previous resummation studies often neglected "particle-hole" ring diagrams, raising questions about the robustness of the predicted transition.
• Goal: To compute the free energy $F(T, V, N_+, N_-)$ including the complete third-order contribution $(k_F a_0)^3$ and to perform an all-orders resummation of both particle-particle and particle-hole ring diagrams to determine the stability of the ferromagnetic phase.

2. Methodology

The authors employ a systematic thermal perturbative expansion using the imaginary time formalism within the framework of Effective Field Theory (EFT).

• Formalism:
  • The system is described by a grand canonical ensemble with two chemical potentials ($\mu_+$ and $\mu_-$) for spin-up and spin-down fermions.
  • The interaction is modeled by a local potential expanded in a series of local operators (Eq. 1), where the leading term is proportional to the coupling $C_0$.
  • The free energy is derived from the thermodynamic potential $\Omega$ via Legendre transformation, utilizing a prescription that avoids the cumbersome inversion of chemical potentials by using zeroth-order chemical potentials for higher-order corrections.
• Renormalization:
  • The local interaction introduces ultraviolet (UV) divergences. These are handled by expressing the bare coupling $C_0$ in terms of the physical s-wave scattering length $a_0$ and a cutoff $\Lambda$ (Eq. 12).
  • Divergences are systematically canceled by matching the perturbative expansion to the physical scattering amplitude.
• Diagrammatic Resummation:
  • The authors identify two infinite sets of Feynman diagrams:
    1. Particle-Particle (PP) Rings: Composed of "A-blocks" (two-particle propagators).
    2. Particle-Hole (PH) Rings: Composed of "B-blocks" (particle-hole propagators).
  • They derive closed-form analytical expressions for the sum of these infinite series at finite temperature.
• Technical Resolution:
  • A key technical hurdle addressed is the handling of singularities in the third-order integrals. The authors demonstrate that naive symmetrization of integrals leads to incorrect results due to spurious singularities. They resolve this by using principal value prescriptions and carefully separating divergent and finite parts, ensuring the cancellation of UV divergences up to order $(k_F a_0)^3$.

3. Key Contributions

1. Third-Order Free Energy Derivation: The paper provides the first complete derivation of the $(k_F a_0)^3$ corrections to the free energy $F$ at finite temperature, extending previous ground-state energy calculations.
2. All-Orders Resummation: The authors derive closed-form formulas (Eqs. 47 and 60) for the resummation of both particle-particle and particle-hole ring diagrams to all orders in the gas parameter $k_F a_0$.
3. Finite Temperature Analysis: Unlike many previous works focused on $T=0$, this study explicitly calculates the temperature dependence of the free energy, allowing for the investigation of the thermal profile of the phase transition.
4. Critical Re-evaluation of the Stoner Transition: By including the particle-hole contribution, the authors fundamentally challenge the existence of the itinerant ferromagnetic phase in the specific regime of large positive scattering lengths (Feshbach resonance side).

4. Results

The numerical evaluation of the resummed free energy yields the following critical findings:

• Particle-Particle (PP) Only: When only the particle-particle ring diagrams are included (resummed), the results reproduce previous findings (e.g., Ref. [18]).
  • At $T=0$, a continuous phase transition to a ferromagnetic state occurs at a critical gas parameter $(k_F a_0)_{cr} \approx 0.79$.
  • As temperature increases, the critical value shifts to higher $k_F a_0$.
  • A "reentrant" transition to the paramagnetic state is observed at very large $k_F a_0$ (associated with the formation of bound states/dimers).
• Particle-Particle + Particle-Hole (PP + PH): When the contribution of the particle-hole ring diagrams is added to the free energy:
  • The phase transition disappears completely.
  • The free energy $F(P)$ exhibits a single global minimum at zero polarization ($P=0$) for all values of $k_F a_0$ and all temperatures studied.
  • The particle-hole contribution is of the same order of magnitude as the particle-particle contribution but has the opposite sign, effectively canceling the instability that drives the ferromagnetic transition.

5. Significance and Implications

• Theoretical Impact: The results suggest that the prediction of itinerant ferromagnetism in dilute Fermi gases, based solely on particle-particle ring resummations, is an artifact of neglecting particle-hole fluctuations. In a truly repulsive system where all scattering lengths are comparable, the transition might still occur, but in the specific regime of Feshbach resonances (where $a_0$ is artificially made very large and positive while other parameters remain small), the transition is suppressed.
• Experimental Relevance: This provides a theoretical explanation for the failure of experimental attempts to observe itinerant ferromagnetism in cold atomic gases (e.g., Refs. [24-26]). The absence of the transition in the metastable state (before dimer formation) implies that the system remains paramagnetic even at large scattering lengths.
• Metastability: The study highlights that the state being modeled is likely a metastable state on the BEC side of the Feshbach resonance. The disappearance of the transition suggests that the "textbook" mean-field expectation of ferromagnetism does not hold when higher-order fluctuations (specifically particle-hole rings) are properly accounted for in this regime.
• Future Directions: The authors conclude that further theoretical work is needed, particularly extending computations to fourth order $(k_F a_0)^4$ and including contributions from p-wave scattering lengths ($a_1$) and effective ranges ($r_0$) to fully characterize systems with truly repulsive potentials.

In summary, the paper demonstrates that the inclusion of particle-hole ring diagrams in the free energy of a Fermi gas with large s-wave scattering length completely stabilizes the paramagnetic phase, thereby eliminating the predicted Stoner ferromagnetic transition in this specific physical regime.