On the phase structure of massless many-flavour QCD with staggered fermions

This paper investigates the phase structure of massless many-flavour QCD using unimproved staggered fermions to map chiral phase boundaries in the bare lattice parameter space, concluding that the continuum chiral transition is second-order for all $N_f$ below the conformal window and proposing methods to identify the conformal onset despite lattice artifacts.

The Gist

Imagine the universe is built from tiny, invisible Lego bricks. In the world of particle physics, the most important bricks are quarks, which stick together to form protons and neutrons. The "glue" that holds them together is a force called the Strong Force (or QCD).

Usually, these quarks have a specific "flavor" (like up, down, strange, etc.). But in this paper, the scientists are playing a thought experiment: What happens if we keep adding more and more types of massless (weightless) quark flavors to the mix?

They are trying to find a "tipping point." If you add too many flavors, the rules of the game change completely. The universe stops breaking symmetry (a fancy way of saying the particles stop behaving differently at low energy) and enters a strange, scale-free state called the Conformal Window.

Here is the breakdown of their journey, using simple analogies.

1. The Goal: Finding the "Magic Number"

The scientists want to know: What is the maximum number of quark flavors ($N_f$) we can have before the laws of physics change?

• Below the limit: The universe behaves normally. Quarks break symmetry, and we get the matter we see today.
• Above the limit (The Conformal Window): The universe becomes "conformal." It loses its scale. It's like a fractal pattern that looks the same whether you zoom in or out. No matter how hot or cold it gets, the quarks never settle down into the structures we know.

The big question is: Is the magic number 7, 8, or something else?

2. The Problem: The "Foggy Mirror"

To study this, they use supercomputers to simulate the universe on a grid (a lattice). But there's a catch:

• The Grid is Rough: Computers can't simulate a perfect, smooth universe. They have to use a grid with "pixels" (lattice spacing).
• The Fog: Because the grid is rough, two different things happen that look the same to the computer:
  1. The Real Thing: The actual thermal transition (like water boiling into steam) that happens in the real universe.
  2. The Artifact: A fake transition caused just by the grid being too coarse (like a pixelated image glitching).

It's like trying to watch a movie on a very low-resolution screen. You can't tell if the character is crying because of the plot (real physics) or because the pixels are just blurry (lattice artifact).

3. The Strategy: Mapping the Terrain

The authors decided to map out the entire "terrain" of this simulation, looking at four variables:

• $\beta$ (Beta): How fine the grid is (resolution).
• $N_\tau$: The "time" steps in the simulation.
• $m$: The mass of the quarks (they are testing massless, but start with a little weight to see how it changes).
• $N_f$: The number of flavors.

They are essentially drawing a 3D map to see where the "Real Transition" lives and where the "Fake Transition" lives.

4. The Discovery: Two Different Worlds

They found two very different landscapes depending on how many flavors they used.

Scenario A: The "Normal" World ($N_f \le 7$)

Imagine a mountain range.

• On coarse grids (low resolution), the mountain looks like a sharp cliff (a First-Order Transition). It's a sudden jump.
• As they increase the resolution (make the grid finer), the cliff smooths out.
• Eventually, on a fine grid, the cliff becomes a gentle slope (a Second-Order Transition).
• Conclusion: For 7 or fewer flavors, the "cliff" was just an illusion caused by the low-resolution grid. In the real, smooth universe, the transition is a gentle slope. The physics is stable.

Scenario B: The "Strange" World ($N_f = 8$)

Now, they added the 8th flavor. The map changed drastically.

• They tried to find that gentle slope (the second-order transition) leading to the real universe.
• But it disappeared!
• Instead of a slope leading to the "Real Universe," the path hit a wall. The "Real Transition" got blocked by the "Fake Transition" (the bulk transition) before it could reach the massless limit.
• The Metaphor: Imagine trying to walk to the top of a mountain (the real universe). For 7 flavors, you can walk up a ramp. For 8 flavors, the ramp ends abruptly at a cliff edge, and you fall into a pit of "lattice artifacts." You can't reach the top.

5. The Big Conclusion: The Conformal Window?

This is the exciting part.

• If the path to the "Real Universe" is blocked for $N_f = 8$, it suggests that the Real Universe doesn't exist for 8 flavors.
• In other words, if you have 8 massless flavors, the universe never breaks symmetry. It stays in that weird, scale-free Conformal Window forever.
• The "cliff" they see in the simulation isn't a glitch; it's a sign that the physics has fundamentally changed.

Why Does This Matter?

This isn't just about quarks.

• Composite Higgs: Some theories suggest that the Higgs boson (the particle that gives us mass) isn't a fundamental brick, but a composite object made of these new, strange quarks.
• New Physics: If we can find a theory that sits right on the edge of this "Conformal Window," it could explain why the universe is the way it is, potentially solving mysteries about dark matter or why the Higgs is so light.

Summary in One Sentence

The scientists used supercomputers to map the "landscape" of particle physics and found that while 7 flavors behave normally, 8 flavors might push the universe into a strange, new state where the usual rules of matter no longer apply, hinting that the "Conformal Window" might start right at 8 flavors.

Technical Summary

1. Problem Statement

The primary objective of this work is to determine the critical number of massless fermion flavors ($N_f^*$) at which Quantum Chromodynamics (QCD) enters the conformal window. In this regime, the theory becomes chirally symmetric in the vacuum due to the emergence of an infrared fixed point (Banks-Zaks), and chiral symmetry breaking ceases to occur.

Determining $N_f^*$ via lattice simulations is challenging because:

• Simulations are performed at finite lattice spacing ($a$), finite quark mass ($am$), and finite temporal extent ($N_\tau$).
• Two distinct types of phase transitions complicate the analysis:
  1. Thermal transitions: Physical transitions related to chiral symmetry restoration at finite temperature.
  2. Bulk transitions: Unphysical transitions driven by lattice artifacts at strong coupling, which obscure the continuum physics.
• Previous studies indicated that first-order chiral transitions observed on coarse lattices disappear as the lattice spacing is reduced, suggesting a second-order transition in the continuum for $N_f \leq 7$. However, the behavior for $N_f = 8$ remains ambiguous, with conflicting evidence regarding whether it marks the onset of the conformal window.

2. Methodology

The authors employ a systematic mapping of the chiral phase boundaries within the extended bare lattice parameter space defined by:

• Gauge coupling: $\beta = 6/g^2$
• Number of flavors: $N_f$
• Bare quark mass: $am$
• Temporal lattice extent: $N_\tau$

Simulation Setup:

• Action: Standard Wilson gauge action and unimproved staggered fermions.
• Lattice Sizes: $N_\tau \times N_\sigma^3$ with aspect ratios $N_\sigma/N_\tau \in \{2, 3, 4\}$.
• Software: OpenCL framework CL2QCD and bash tool BaHaMAS for job management.
• Parameters: Simulations cover $N_f \in [6, 10]$ and $N_\tau \in [6, 16]$.

Analysis Technique:

• Order Parameter: The chiral condensate $\langle \bar{\psi}\psi \rangle$.
• Cumulants: The order of the transition is determined using generalized cumulants of the condensate distribution:
  • Skewness ($B_3$): Used to locate the (pseudo-)critical coupling $\beta_{pc}$ where $B_3 = 0$.
  • Kurtosis ($B_4$): Used to determine the order of the transition via finite-size scaling.
    • $B_4 = 1$: First-order transition.
    • $B_4 \approx 1.6044$: Second-order transition ($Z_2$ universality class).
    • $B_4 = 3$: Crossover.
• Strategy: The authors scan $\beta$ at fixed $N_f$ and $N_\tau$ for various masses $am$ to locate the transition lines and extrapolate to the chiral limit ($am \to 0$) and continuum limit ($N_\tau \to \infty$).

3. Key Contributions and Results

A. Phase Structure for $N_f \leq 7$ (Non-Conformal)

• Tricritical Point: For all $N_f \leq 7$, the authors confirm the existence of a tricritical point $(\beta_{tric}, N_{tric}^\tau)$ in the chiral limit.
• Transition Order:
  • For coarse lattices ($N_\tau < N_{tric}^\tau$), the chiral transition is first-order.
  • For finer lattices ($N_\tau \geq N_{tric}^\tau$), the transition becomes second-order ($Z_2$ universality).
• Continuum Limit: The first-order region is identified as a cutoff effect. In the continuum chiral limit ($am=0, \beta \to \infty, N_\tau \to \infty$), the transition is strictly second-order for all $N_f \leq 7$.
• Bulk Transition: A zero-temperature bulk transition separates the physical weak-coupling regime from a strong-coupling lattice artifact regime. For $N_f > 6$, this bulk transition becomes non-analytic (first-order at small masses) but terminates at a second-order critical point.

B. Phase Structure for $N_f = 8$ (Candidate for Conformal Window)

• Interference of Transitions: At $N_f = 8$, the thermal transition line is increasingly interfered with by the bulk transition as $N_\tau$ increases.
• Disappearance of Tricriticality: Unlike $N_f \leq 7$, no tricritical point is observed in the chiral limit.
  • For $N_\tau = 8$, a first-order thermal transition is seen at small masses, but it terminates on the bulk transition line at finite mass.
  • For $N_\tau \geq 10$, no first-order or second-order thermal transition is found in the accessible mass range; only a crossover remains.
• Hypothesis: The thermal transition line for $N_f=8$ appears to "end" on the bulk transition at finite mass $am > 0$ rather than extending to the chiral limit ($am=0$). This suggests that in the true chiral limit, the thermal transition is absent, consistent with the theory being conformal.

C. Proposed Identification Strategy for the Conformal Window

The authors propose a novel method to distinguish conformal from non-conformal theories using lattice simulations away from the chiral limit:

1. Non-Conformal ($N_f < N_f^*$): The thermal transition lines (crossover/second-order) for increasing $N_\tau$ converge toward a second-order line that connects to the chiral limit ($am=0$) at finite $\beta$.
2. Conformal ($N_f \geq N_f^*$): The thermal transition lines for increasing $N_\tau$ become steeper and intersect the bulk transition at finite mass ($am > 0$). Consequently, no thermal transition exists in the chiral limit ($am=0$) for any finite $N_\tau$; the thermally broken phase is disconnected from the continuum limit.

4. Significance

• Clarification of Lattice Artifacts: The work provides a coherent picture of how unphysical bulk transitions obscure the physical thermal transitions in many-flavor QCD, offering a roadmap to disentangle them.
• Evidence for $N_f=8$: The results strongly suggest that $N_f=8$ lies at or near the onset of the conformal window, as the thermal transition fails to reach the chiral limit, unlike the cases $N_f \leq 7$.
• Methodological Advance: The proposed strategy of analyzing the intersection of thermal and bulk transition surfaces in the $(N_\tau, \beta, am)$ space allows for the identification of the conformal window without requiring simulations strictly at the chiral limit, which is computationally expensive.
• Broader Impact: Understanding the conformal window is crucial for Beyond Standard Model (BSM) physics, particularly for composite Higgs scenarios and strongly coupled extensions of the Standard Model.

5. Conclusion

The paper concludes that continuum QCD with massless quarks exhibits a second-order chiral phase transition for all $N_f$ up to the onset of the conformal window. For $N_f=8$, the absence of a thermal transition connecting to the chiral limit, combined with the "hiding" of the transition behind the bulk phase at finite mass, serves as a strong indicator that $N_f=8$ may be the lower edge of the conformal window. Further simulations at larger $N_\tau$ are required to definitively confirm this scenario.