Magnetic susceptibility of a hot hadronic medium and quark degrees of freedom near the QCD cross-over point

This paper proposes a quark-meson framework incorporating temperature-dependent quark masses and anomalous magnetic moments to reconcile lattice QCD magnetic susceptibility data with theoretical models, demonstrating that quark degrees of freedom must emerge significantly below the QCD cross-over temperature (around 120 MeV) to explain the observed paramagnetism.

The Gist

The Big Picture: The "Magnetic Mood" of Hot Matter

Imagine you have a giant pot of soup. This isn't just any soup; it's the "primordial soup" of the universe, made of the tiniest building blocks of matter (quarks and gluons) that existed just after the Big Bang. As this soup cools down, it changes state, much like water turning into ice.

Physicists are trying to understand exactly how this soup behaves when you stick a giant magnet near it. Specifically, they want to know: Does the soup get attracted to the magnet (paramagnetism) or repelled by it (diamagnetism)?

This property is called Magnetic Susceptibility.

The Problem: The "Recipe" Didn't Work

For a long time, physicists used a standard recipe to predict how this soup behaves. They called it the Hadron Resonance Gas (HRG) model.

• The Analogy: Think of the HRG model as a cookbook that says, "Below a certain temperature, the soup is made entirely of whole, solid vegetables (protons, neutrons, pions). Above that temperature, the vegetables melt into a smoothie of individual ingredients (quarks)."

The scientists took this recipe and calculated the magnetic mood of the soup.

• The Result: The recipe predicted that the soup should be strongly repelled by magnets (diamagnetic) at temperatures just below the "melting point."
• The Reality Check: When they compared this to super-accurate computer simulations (called Lattice QCD, which act like a digital microscope), the recipe was wrong. The real soup wasn't repelled as strongly as the recipe said. In fact, at temperatures just below the melting point, the soup started acting like it was being attracted to the magnet.

The Mystery: Why did the "vegetable-only" recipe fail? It suggested that even before the soup fully melted into a smoothie, some of the individual ingredients (quarks) were already floating around, changing the magnetic mood.

The Investigation: Why the Old Recipe Failed

The authors, Rupam Samanta and Wojciech Broniowski, decided to investigate why the "vegetable-only" model failed. They tried to fix the recipe in a few ways:

1. Adding "Spicy" Moments: They added the real-world magnetic "personality" (magnetic moments) of the vegetables.
   • Result: It helped a little, but not enough. The soup was still too repulsive.
2. Adding "Swirls" (Loops): They considered complex interactions where particles swirl around each other (pion-vector-meson loops).
   • Result: This added a tiny bit of attraction, but it was like adding a pinch of salt to a huge pot of soup. It didn't solve the problem.
3. Checking "Squishiness" (Polarizability): They checked if the vegetables got squished by the magnet.
   • Result: Negligible.

The Conclusion: The "vegetable-only" model was fundamentally missing a key ingredient. To explain the data, the soup must contain some free-floating, magnetic ingredients (quarks) even before it fully melts.

The Solution: The "Quark-Meson" Smoothie

The authors proposed a new model: The Quark-Meson Approach.

Instead of waiting for the vegetables to fully melt, they assumed that as the soup heats up, the vegetables start to break apart into their constituent parts (quarks), but the vegetables (pions) still exist alongside them.

• The Analogy: Imagine a crowd of people (the soup).
  • Old Model: Everyone is holding hands in a tight circle (diamagnetic). If you bring a magnet, they all push away.
  • New Model: As it gets hot, some people let go of hands and start running around individually (quarks). These runners are magnetic and want to get closer to the magnet (paramagnetic).
  • The Balance: You have the circle of people pushing away (pions) and the runners pulling in (quarks). The new model calculates exactly how many runners are needed to balance the circle so that the total magnetic mood matches the computer simulation.

How They Did It (The "Fingerprint" Method)

To make their new model work, they needed to know the "weight" (mass) of these running quarks at different temperatures. They couldn't just guess.

They used a clever trick:

1. They looked at other data from the computer simulations regarding how the soup reacts to changes in "baryon number" and "strangeness" (think of these as different types of flavor or charge in the soup).
2. They adjusted the "weight" of the quarks in their model until their predictions matched these other data points perfectly.
3. The Discovery: They found that the quarks get heavier as the soup gets hotter, but they are light enough to be active and magnetic even at temperatures as low as 120 MeV (which is about 1.3 trillion degrees Kelvin!).

The Final Verdict

The paper concludes that our understanding of the "melting point" of matter needs an update.

• Old View: Below the melting point, it's just solid vegetables. Above it, it's a smoothie.
• New View: The transition is a messy overlap. Even when the soup looks like solid vegetables, it's actually a mix of vegetables and free-floating quarks. These quarks are the "secret sauce" that explains why the soup behaves magnetically the way it does.

Why does this matter?
It helps us understand the extreme conditions of the early universe and what happens inside neutron stars. It tells us that the "building blocks" of matter start behaving like individuals much earlier in the heating process than we previously thought.

Summary in One Sentence

The authors discovered that the standard model of hot matter was wrong because it ignored the fact that tiny magnetic particles (quarks) start floating around and changing the soup's magnetic personality long before the soup fully melts into a smoothie.

Technical Summary

1. Problem Statement

The central problem addressed is the discrepancy between Lattice QCD results and the Hadron Resonance Gas (HRG) model regarding the magnetic susceptibility ($\chi_B$) of strongly interacting matter near the QCD cross-over temperature ($T_c \approx 155$ MeV).

• The Discrepancy: Lattice simulations indicate that $\chi_B$ transitions from diamagnetic (negative) at low temperatures to paramagnetic (positive) as $T$ approaches and exceeds $T_c$. However, the standard HRG model, even when augmented with physical hadron magnetic moments and pion-vector-meson loops, predicts a strong diamagnetism that persists well above $T_c$, failing to match the lattice data for $T \gtrsim 120$ MeV.
• The Implication: This mismatch suggests that the HRG model, which assumes the medium consists solely of hadrons below $T_c$, is incomplete. It hints at the presence of quark degrees of freedom (specifically light, paramagnetic constituents) significantly below the cross-over temperature, which are necessary to generate the observed paramagnetism.

2. Methodology

The authors employ a multi-step approach to resolve the discrepancy:

A. Re-evaluation of the HRG Model

• They performed a comprehensive HRG calculation including a complete list of hadronic states and physical magnetic moments (anomalous magnetic moments, $\kappa$).
• They analyzed the contribution of pion–vector-meson loops (e.g., $\pi-\rho$, $\pi-\omega$) to the photon polarization tensor, which generates a paramagnetic correction.
• They evaluated the impact of hadronic magnetic polarizabilities.

B. Development of a Quark-Meson Model

To address the failure of pure HRG, the authors constructed a non-interacting quark-meson model below $T_c$.

• Degrees of Freedom: The model includes constituent quarks (up, down, strange) and mesons (pions, kaons).
• Temperature-Dependent Masses: Instead of assuming fixed constituent masses, the quark masses ($m_l$ for light quarks, $m_s$ for strange quarks) are determined dynamically. They are fitted to Lattice QCD data for baryon-baryon ($\chi_{BB}$) and baryon-strangeness ($\chi_{BS}$) susceptibilities at zero magnetic field.
• Vacuum Contributions: Unlike standard HRG, the model explicitly incorporates vacuum loop contributions (quark and meson loops) to the pressure and susceptibility. These are regularized using a Pauli-Villars scheme with a cutoff $\Lambda \approx 800$ MeV.
• Anomalous Magnetic Moments: The model includes anomalous magnetic moments for constituent quarks ($\kappa_q$), estimated from octet baryon magnetic moments using a non-relativistic quark model.

C. Calculation of Susceptibility

The total magnetic susceptibility is calculated as the sum of thermal and vacuum parts for both quarks and mesons:
$$\chi_B(T) = \sum_{q} [\chi_{q,th} + \chi_{q,vac}] + \sum_{mes} [\chi_{mes,th} + \chi_{mes,vac}]$$
The results are compared against the subtracted lattice data $\chi_B(T) - \chi_B(0)$.

3. Key Contributions and Results

A. Failure of Pure Hadronic Models

• HRG Limitations: Even with physical magnetic moments, the HRG model remains too diamagnetic for $T > 120$ MeV. The inclusion of anomalous moments only raises the curve by 5–10%, insufficient to bridge the gap with lattice data.
• Loop Effects: The contribution from pion–vector-meson loops is found to be paramagnetic but small (10–20% enhancement), agreeing qualitatively with Chiral Perturbation Theory but insufficient to explain the lattice data alone.
• Polarizabilities: The contribution from hadronic magnetic polarizabilities is found to be negligible ($\sim 0.002$).

B. Success of the Quark-Meson Approach

• Mass Fitting: By fitting the quark masses to $\chi_{BB}$ and $\chi_{BS}$ lattice data, the authors found that light quark masses at $T_c$ are significantly larger ($\sim 550$ MeV) than traditional constituent masses ($\sim 300$ MeV), aligning with quasiparticle models. The mass splitting $m_s - m_l$ remains stable at $\sim 120-130$ MeV.
• Reproducing $\chi_B$: The quark-meson model, incorporating temperature-dependent masses, vacuum terms, and quark anomalous moments, successfully reproduces the lattice $\chi_B$ data from $T \approx 135$ MeV up to $T_c$.
• Mechanism: The transition from diamagnetism to paramagnetism is driven by the thermal quark contribution, which is strongly paramagnetic. The meson (pion/kaon) contribution remains diamagnetic, preserving the low-temperature behavior.
• Role of Vacuum Terms: The vacuum quark loop contribution is crucial; without it, the model cannot satisfy the lattice constraints. The vacuum term is flat at low $T$ and decreases above $T \approx 155$ MeV as quark masses drop.

C. Anomalous Magnetic Moments

The inclusion of quark anomalous magnetic moments ($\kappa_q$) significantly enhances the paramagnetic component, improving the agreement with lattice data.

4. Significance and Conclusions

• Quark Degrees of Freedom Below $T_c$: The primary conclusion is that quark degrees of freedom must exist significantly below the cross-over temperature, extending down to $\approx 120$ MeV. The medium is not purely hadronic in this region; it contains "melted" baryons or constituent quarks that provide the necessary paramagnetism.
• Unified Description: The quark-meson framework is the first model to simultaneously describe:
  1. Lattice data for conserved charge susceptibilities ($\chi_{BB}, \chi_{BS}$).
  2. Lattice data for magnetic susceptibility ($\chi_B$).
  3. The transition from diamagnetic to paramagnetic behavior.
• Implications for Heavy Ion Collisions: The findings challenge the standard interpretation of the hadronic phase in ultra-relativistic heavy ion collisions. If quarks are relevant down to 120 MeV, the successful application of HRG to particle yields (chemical freeze-out) might imply a mechanism where quarks fragment into hadrons at a later stage, rather than the system being purely hadronic throughout the freeze-out process.
• Future Directions: The authors suggest that more uniform analyses using models like the Polyakov–Nambu–Jona-Lasinio (PNJL) model, treated at the one-meson-loop level, are necessary to fully integrate these findings with all available lattice data.

In summary, the paper demonstrates that the Hadron Resonance Gas model is insufficient to describe the magnetic properties of QCD matter near the cross-over. A quark-meson approach that includes temperature-dependent quark masses and vacuum effects is required to reconcile theoretical models with Lattice QCD results, providing strong evidence for the early emergence of quark degrees of freedom.