Dirac mode localization in QCD near the crossover temperature

Using staggered fermions on the lattice, the study demonstrates that localized low-lying Dirac eigenmodes emerge in QCD at a temperature ($T_{\mathrm{loc}}$) between 155 and 158 MeV, which aligns precisely with the pseudocritical crossover temperature determined by chiral condensate and light-quark susceptibility.

The Gist

Imagine the universe as a giant, bubbling pot of soup. For most of its history, this soup was so hot that the ingredients (quarks and gluons) were free-floating, swimming around chaotically. This is called the Quark-Gluon Plasma. As the universe cooled down, this soup "froze" into solid chunks, binding the quarks together into particles like protons and neutrons. This is the world we live in today.

The big question physicists have is: Exactly when and how does this transition happen? It's not like water freezing at exactly 0°C; it's more like a gradual shift from a liquid to a gel.

This paper, by Giordano, Kovács, and Pittler, investigates a very specific, microscopic clue to pinpoint exactly when this transition occurs. They look at the behavior of "Dirac eigenmodes," which is a fancy way of saying: how the quantum waves of the quarks move around inside this hot soup.

Here is the breakdown using simple analogies:

1. The "Crowded Room" vs. The "Empty Hall"

Think of the quarks' quantum waves as people walking through a building.

• At Low Temperatures (The Solid World): The building is a giant, empty hall. The people (waves) can walk freely from one end to the other without bumping into anything. They are delocalized. They are spread out everywhere.
• At High Temperatures (The Plasma): Suddenly, the building gets filled with chaotic, moving obstacles (fluctuating gauge fields). Now, the people can't walk freely. They get stuck in small corners or specific rooms. They are localized. They are trapped.

The authors are looking for the exact moment the building fills up enough that the people stop walking freely and start getting stuck.

2. The "Mobility Edge" (The Security Guard)

In the middle of the spectrum of these waves, there is a theoretical "security guard" called the Mobility Edge.

• Below this guard (low energy waves), the waves are stuck (localized).
• Above this guard (high energy waves), the waves are free (delocalized).

The researchers wanted to find the temperature where this security guard first appears. Before this temperature, there is no guard; everyone is free. After this temperature, the guard shows up, trapping the low-energy waves.

3. The "Fingerprint" of Chaos

How do they know if the waves are stuck or free without watching every single one? They use a statistical trick called Random Matrix Theory.

• Imagine listening to the footsteps of the people.
• If they are walking freely in a large hall, their footsteps have a very specific, rhythmic pattern (like a well-organized marching band).
• If they are stuck in small rooms, their footsteps are random and chaotic (like a crowd at a mosh pit).

By analyzing the "rhythm" of the quantum waves in their computer simulations, the authors could tell exactly when the pattern switched from "organized marching" to "chaotic mosh pit."

4. The Big Discovery

The team ran massive computer simulations of the early universe soup at different temperatures. They found:

• Below 155 MeV: The waves are all free. No "mobility edge" exists. The system is in the "hadronic" (solid) phase.
• Between 155 MeV and 158 MeV: Suddenly, the "mobility edge" appears! The low-energy waves get trapped.
• Above 158 MeV: The trapping effect gets stronger.

The "Aha!" Moment:
This temperature range (155–158 MeV) is exactly the same temperature where other famous measurements say the transition happens.

• It matches the temperature where the "chiral condensate" (a measure of how tightly quarks are bound) changes.
• It matches the temperature where the "susceptibility" (how easily the system reacts to changes) peaks.

Why Does This Matter?

For a long time, physicists have suspected that two major changes in the universe happen at the same time:

1. Deconfinement: Quarks stop being stuck in protons and start floating freely.
2. Chiral Symmetry Restoration: A specific type of quantum symmetry is "broken" in the cold world but "restored" in the hot world.

It was a mystery why these two things happened at the exact same temperature. This paper suggests a deep, geometric link: The moment the quarks stop being able to roam freely (localization) is the exact same moment the universe switches phases.

The Bottom Line

The authors didn't just guess; they found a "geometric" switch in the math of the universe that flips at the exact same temperature as the thermodynamic switch. It's like finding that the moment a crowd stops dancing and starts sitting down is the exact same moment the music stops.

They have confirmed that the localization of quantum waves is a reliable, precise thermometer for the birth of the Quark-Gluon Plasma, and it agrees perfectly with everything else we know about the early universe.

Technical Summary

1. Problem Statement

The nature of the finite-temperature transition in Quantum Chromodynamics (QCD) at vanishing chemical potential is an analytic crossover rather than a sharp phase transition. Consequently, there is no single, sharply defined critical temperature ($T_c$) where thermodynamic properties change singularly.

• The Challenge: Standard order parameters, such as the chiral condensate (related to chiral symmetry restoration) and the Polyakov loop (related to deconfinement), are only approximate in real-world QCD (with physical quark masses). While they indicate a transition region, they cannot unambiguously define a specific temperature separating the hadronic and quark-gluon plasma (QGP) phases.
• The Hypothesis: The authors investigate whether the localization properties of low-lying Dirac eigenmodes can serve as a gauge-invariant, geometric observable that displays a singular behavior at a well-defined temperature ($T_{loc}$). Previous studies suggested that localized low modes appear at high temperatures, separated from delocalized bulk modes by a "mobility edge" in the spectrum. The goal is to determine if $T_{loc}$ coincides with the pseudocritical temperatures ($T_{pc}$) derived from thermodynamic observables.

2. Methodology

The study employs Lattice QCD simulations using rooted staggered fermions with $N_f = 2+1$ flavors at physical quark masses.

• Simulation Setup:

  • Action: Tree-level Symanzik improved Wilson gauge action.
  • Smearing: Two steps of stout smearing ($\rho=0.15$) applied to gauge fields in the fermion determinant and spectrum calculation to improve discretization effects.
  • Ensembles: Simulations performed at temperatures ranging from $150$ to $184$ MeV.
  • Finite Volume/Spacing Control:
    • Spatial aspect ratios $N_s/N_t = 6, 8, 10$ were used to control finite-volume effects.
    • Lattice spacings corresponding to temporal extents $N_t = 6, 8, 10$ were used at $T=165$ MeV to control discretization effects.
  • Algorithm: The Krylov-Schur algorithm was used to compute the low-lying positive eigenvalues of the staggered Dirac operator.

• Analysis Technique (Spectral Statistics):

  • Instead of directly calculating the Inverse Participation Ratio (IPR), the authors utilized spectral statistics, which are numerically more efficient and robust.
  • Unfolding: The spectrum was unfolded to remove system-specific density variations, mapping eigenvalues $\lambda_i$ to $x_i$ such that the mean level spacing is unity.
  • Level Spacing Distribution: The probability distribution of unfolded level spacings ($s$) was analyzed.
    • Delocalized regime: Follows Random Matrix Theory (RMT) statistics (Wigner surmise for the unitary class).
    • Localized regime: Follows Poisson statistics.
  • Mobility Edge Determination: The authors monitored the integrated unfolded level spacing distribution ($I_{s_0}$) at a specific threshold $s_0 \approx 0.508$.
    • $I_{s_0}^{RMT} \approx 0.117$
    • $I_{s_0}^{Poisson} \approx 0.398$
    • $I_{s_0}^{crit} \approx 0.1966$ (Critical value at the mobility edge).
  • The mobility edge ($\lambda_c$) is identified as the spectral point where $I_{s_0}$ crosses the critical value.

3. Key Contributions

1. Direct Measurement: This is the first direct measurement of the localization temperature ($T_{loc}$) in QCD, rather than relying on extrapolations from higher temperatures.
2. Gauge-Invariant Geometric Transition: The work establishes a "geometric" transition defined by the appearance of localized modes, providing a sharp definition for the onset of the high-temperature phase that is independent of thermodynamic ambiguities.
3. Resolution of Systematic Effects: The study rigorously addresses systematic errors inherent to staggered fermions, specifically:
   • Taste Symmetry: Demonstrated that working with sufficiently large volumes ($N_s/N_t \ge 8$) washes out the effects of taste multiplets on spectral statistics.
   • Renormalization: Confirmed that the ratio of the mobility edge to the bare quark mass ($\lambda_c/m_{ud}$) is a renormalization-group invariant quantity with a non-zero continuum limit.

4. Results

• Localization Temperature ($T_{loc}$):
  • At $T = 150$ MeV (below $T_{pc}$): All low-lying modes are delocalized ($I_{s_0} \approx I_{s_0}^{RMT}$). No mobility edge is found.
  • At $T = 158$ MeV (above $T_{pc}$): Localized low modes appear. The mobility edge is clearly detected, and $I_{s_0}$ crosses the critical value near the origin.
  • At $T = 155$ MeV: The system is in a state where $I_{s_0}$ remains close to the RMT value, indicating no mobility edge yet.
  • Conclusion: The critical temperature for localization lies in the range:
    $$155 \text{ MeV} \le T_{loc} \le 158 \text{ MeV}$$
• Comparison with Thermodynamics:
  • This range is in excellent agreement with the pseudocritical temperatures determined from:
    • Renormalized chiral condensate: $T_{pc} = 155(4)$ MeV.
    • Light-quark chiral susceptibility: $T_{pc} = 157(4)$ MeV.
    • Static quark entropy: $T_{pc} \approx 153$ MeV.
• Continuity: Extrapolation of the mobility edge $\lambda_c$ as a function of temperature suggests it vanishes smoothly as $T \to T_{loc}$ from above, consistent with the analytic nature of the crossover.

5. Significance

• Unification of Phenomena: The results provide strong evidence that chiral symmetry restoration and deconfinement (in the loose sense applicable to physical QCD) are driven by the same microscopic mechanism. This mechanism is reflected in the behavior of the low Dirac modes.
• Geometric Definition of Deconfinement: The study validates the idea that the appearance of localized low Dirac modes can serve as a precise, gauge-invariant definition for the transition to the deconfined phase.
• Theoretical Insight: It bridges the gap between the ordering of the Polyakov loop (which opens a pseudogap in the Dirac spectrum) and the localization of eigenmodes. The appearance of the mobility edge signals that the Polyakov loop has ordered sufficiently to localize the low modes.
• Methodological Benchmark: The paper sets a standard for using spectral statistics to study phase transitions in lattice QCD, demonstrating how to control taste-breaking artifacts and finite-volume effects to extract precise critical temperatures.

In summary, the paper confirms that the "localization temperature" of Dirac eigenmodes coincides with the thermodynamic pseudocritical temperature, reinforcing the deep connection between the spectral properties of the Dirac operator and the macroscopic phases of QCD.