Quark Number Susceptibilities and Conserved Charge Fluctuations in $(2+1)$-flavor QCD with Möbius domain-wall fermions (MDWF)

This paper calculates second- and fourth-order conserved-charge fluctuations in $(2+1)$-flavor QCD using Möbius domain-wall fermions to study quark number susceptibilities across different quark masses and temperatures, comparing the results to hadron resonance gas models and Stefan–Boltzmann limits.

The Gist

The Cosmic Soup: A Recipe for the Early Universe

Imagine you are trying to understand the very first moments of the universe—a time when everything was so hot and dense that atoms couldn't exist. Instead, the universe was a "primordial soup" called Quark-Gluon Plasma (QGP). In this soup, the fundamental building blocks of matter (quarks) were swimming freely, rather than being locked inside protons and neutrons.

Scientists want to know exactly how this soup changes as it cools down. Does it turn from a liquid to a solid instantly (like water freezing into ice), or does it change smoothly (like honey thickening as it cools)? This paper uses some of the world's most powerful supercomputers to simulate this "cooling soup" and figure out its recipe.

---

1. The Tool: The Digital Microscope (Lattice QCD)

To study this, scientists can't just look through a telescope; they have to build a universe inside a computer. This is called Lattice QCD.

Think of it like a digital Minecraft world. Instead of a continuous space, the scientists create a grid (a "lattice") of points. They place the laws of physics on these points and let the simulation run to see how particles behave.

The researchers in this paper used a special, high-quality "lens" called Möbius Domain Wall Fermions.

• The Analogy: Imagine trying to photograph a fast-moving hummingbird. If your camera is cheap, the bird looks like a blurry smudge (this is what happens with older methods called "staggered fermions"). The Möbius method is like having a high-speed, ultra-HD camera that captures the bird’s wings with incredible clarity and precision, even when things are moving incredibly fast.

2. The Measurement: The "Flavor" Fluctuations

The scientists are looking at "Conserved Charge Fluctuations." That sounds intimidating, but think of it as measuring the "jitteriness" of the soup.

Imagine you have a bowl of soup containing three ingredients: salt (Baryon number), pepper (Electric charge), and spice (Strangeness).

• If the soup is a calm, organized liquid, the ingredients stay mostly where they are.
• If the soup is boiling violently or undergoing a massive change, the ingredients will suddenly "jitter" or clump together in unpredictable ways.

By measuring these "jitters" (fluctuations), scientists can tell exactly when the soup is transitioning from the chaotic Quark-Gluon Plasma back into the organized matter (protons and neutrons) that makes up our world today.

3. The Discovery: The Smooth Transition

What did they find?

The "Smooth Transition" Result:
The simulation showed that the transition isn't a sudden "snap." It’s a smooth crossover.

• The Analogy: It’s not like a light switch flipping from OFF to ON. It’s more like dimming a light bulb or a sunset. The universe didn't suddenly "become" matter; it gradually shifted from a wild, free-flowing state into a structured one.

The "Hadron Gas" Comparison:
At lower temperatures, the particles behave exactly like a "Hadron Resonance Gas"—a predictable collection of particles like protons and pions.

• The Analogy: Imagine a crowded dance floor. At high temperatures, everyone is running around wildly and bumping into each other (the Quark-Gluon Plasma). As it cools, people start to pair up and dance in specific patterns (the Hadron Gas). The researchers found that their digital simulation matched the mathematical "dance patterns" perfectly.

4. Why does this matter?

By understanding these fluctuations, we are essentially reading the "DNA of the Big Bang."

When scientists smash atoms together in giant machines like the Large Hadron Collider (LHC), they are trying to recreate this primordial soup. This paper provides the "Gold Standard" guidebook. It tells experimental scientists: "If you see this specific amount of 'jitter' in your particle detector, it means you have successfully recreated the soup of the early universe."

Summary in a Nutshell

Scientists used a super-precise digital simulation to study how the universe's "primordial soup" cooled down. They confirmed that the transition from wild quarks to stable matter was a smooth, gradual process, and they provided a high-definition map that helps us understand the very birth of everything.

Technical Summary

Technical Summary: Quark Number Susceptibilities and Conserved Charge Fluctuations in (2 + 1)-flavor QCD with Möbius Domain Wall Fermions

1. Problem Statement

A central goal in strong-interaction physics is mapping the phase structure of Quantum Chromodynamics (QCD) at finite temperature and baryon density. While it is established that the transition from the hadronic phase to the Quark-Gluon Plasma (QGP) is a smooth crossover at zero baryon chemical potential ($\mu_B=0$), the existence of a critical endpoint at higher densities remains a major open question.

Experimental probes in heavy-ion collisions (RHIC, LHC) measure event-by-event fluctuations of conserved charges (baryon number $B$, electric charge $Q$, and strangeness $S$). To interpret this data, theoretical "baselines" are required. Lattice QCD provides these via derivatives of the pressure with respect to chemical potentials. However, previous lattice studies have predominantly used staggered fermions, which suffer from "taste-symmetry breaking" at finite lattice spacing. This discretization artifact artificially increases the pion mass, potentially distorting low-temperature observables—particularly electric-charge fluctuations, which are highly sensitive to the pion sector.

2. Methodology

The JLQCD Collaboration employs Möbius Domain Wall Fermions (MDWF) to perform these calculations. MDWF offer superior control over chiral symmetry at finite lattice spacing compared to staggered fermions, providing a cleaner probe of light-hadron dynamics.

• Lattice Setup: The study uses (2+1)-flavor QCD (two light quarks, one strange quark) along a Line of Constant Physics (LCP).
• Mass Regimes:
  • Heavier-than-physical: $m_l/m_s = 1/10$ (to study lattice-spacing effects using $N_\tau = 12$ and $16$).
  • Physical point: $m_l/m_s = 1/27.4$ (to study light-quark-mass dependence using $N_\tau = 12$).
• Mathematical Approach: Since direct simulations at $\mu_B \neq 0$ are hindered by the "sign problem," the authors use the Taylor expansion method around $\mu=0$. They calculate quark number susceptibilities ($\chi_{ijk}^{uds}$) up to the fourth order.
• Numerical Techniques: To manage the high computational cost and stochastic noise of disconnected diagrams, they utilize:
  • Gaussian random sources with dilution schemes (spin and time-slice dilution).
  • Mass reweighting to correct for mistuning in the LCP for the heavier mass ensembles.

3. Key Contributions

• Chiral Symmetry Advantage: This work provides a high-precision calculation using a fermion discretization (MDWF) that preserves chiral symmetry better than the standard staggered approach, specifically targeting the resolution of the "pion sector" issue.
• Comprehensive Fluctuation Mapping: The paper computes a wide array of second- and fourth-order conserved-charge cumulants, including diagonal and off-diagonal correlations (e.g., $\chi_{BS}$, $\chi_{QS}$, $\chi_{BQ}$).
• Mass and Cutoff Assessment: It systematically evaluates how results change with the light-quark mass (moving toward the physical point) and the lattice temporal extent $N_\tau$ (assessing discretization effects).

4. Results

• Second-Order Fluctuations:
  • Below the pseudocritical temperature ($T_{pc}$), most observables (electric charge, strangeness, and off-diagonal correlations) are consistent with the Hadron Resonance Gas (HRG) model, specifically the QMHRG2020 calculation.
  • Mass Sensitivity: The electric-charge fluctuation ($\chi_Q^2$) shows the strongest dependence on the light-quark mass, confirming that the pion sector dominates this channel. In contrast, strangeness and baryon-dominated channels are much less sensitive to the light-quark mass.
  • Comparison with Staggered Fermions: At the physical pion mass, MDWF results for $\chi_Q^2$ lie closer to the HRG expectation than staggered-fermion results (HISQ/Stout), suggesting that MDWF better captures the physical pion dynamics at finite lattice spacing.
• Fourth-Order Fluctuations:
  • These are numerically more challenging and noisier. However, for the strange sector ($\chi_S^4, \chi_{QS}^{13}, \chi_{QS}^{22}$), the data show a clear trend toward the Stefan-Boltzmann (SB) limit (the free quark-gas limit) as temperature increases.
  • The QMHRG2020 baseline generally provides a better description of the low-temperature lattice data for strange-sector fourth-order cumulants than the standard PDG-HRG.
• Kurtosis Ratios: The strangeness ratio $R_{S}^{42} = \chi_S^4/\chi_S^2$ is well-resolved and shows the expected rise toward the SB limit, while the electric-charge ratio $R_{Q}^{42}$ remains qualitatively constrained due to higher statistical uncertainty.

5. Significance

This paper establishes a robust, first-principles lattice QCD baseline for conserved-charge fluctuations using a chiral-symmetric fermion action. By demonstrating that MDWF results align more closely with HRG models in the pion-sensitive sector than staggered fermions, the study provides critical evidence that discretization artifacts in previous studies may have distorted the low-temperature hadronic physics. These results are vital for the interpretation of heavy-ion collision data from RHIC and the LHC, particularly in identifying the transition from hadronic matter to the quark-gluon plasma.