A generalized definition of the isothermal compressibility in (2+1)-flavor QCD

This paper introduces a generalized definition of isothermal compressibility in (2+1)-flavor QCD based on conserved charge fluctuations, presenting lattice results that align with hadron resonance gas models and heavy-ion collision data while demonstrating that the compressibility on the pseudo-critical line closely resembles that of an ideal gas.

The Gist

The Big Picture: Squeezing the "Perfect Fluid"

Imagine you have a giant, invisible balloon filled with a super-hot, super-dense soup. This isn't water or air; it's a "soup" made of the fundamental building blocks of matter (quarks and gluons) that existed just a fraction of a second after the Big Bang. Scientists call this strong-interaction matter.

In heavy-ion collision experiments (like smashing gold or lead atoms together at nearly the speed of light), physicists recreate this soup for a tiny split second. They want to know: How easy is it to squeeze this soup?

In physics, the measure of how easy it is to squeeze something is called Isothermal Compressibility ($\kappa_T$).

• High compressibility: The substance is like a fluffy pillow; you can squish it down easily.
• Low compressibility: The substance is like a steel block; it resists being squished.

The Problem: Counting the Wrong Things

For a long time, scientists tried to calculate this "squishiness" by counting the total number of particles in the soup (like counting every single grain of sand in a bucket).

The Analogy: Imagine trying to measure the density of a crowd at a concert.

• The Old Way: You try to count every single person (total particles).
• The Problem: In this specific "soup" (QCD), particles are constantly popping in and out of existence, turning into their anti-particles (like a person turning into their twin and then vanishing). Because of this, the "total count" becomes a moving target that makes the math break down (it goes to infinity) when the soup is perfectly balanced.

The authors of this paper realized that trying to count the total number of people is a bad idea for this specific crowd. Instead, you should count the net difference.

• The New Way: Count the difference between "Team Quark" and "Team Anti-Quark." If there are 100 Quarks and 100 Anti-Quarks, the net number is zero. This number is stable and conserved, even if the individuals keep swapping places.

The Solution: A New Definition

The authors introduced a Generalized Definition of compressibility. Instead of asking, "How does the volume change if we add more total particles?" they asked:
"How does the volume change if we keep the fluctuations (the natural wiggles and jiggles) of the net charge constant?"

Think of it like this:

• Old Method: Trying to measure the stiffness of a trampoline by counting every single person jumping on it.
• New Method: Measuring the stiffness by watching how much the trampoline bounces when the pattern of the jumpers stays the same, even if the individual jumpers change.

The Results: It's Like an Ideal Gas

The team used Lattice QCD (a super-computer simulation that acts like a digital microscope for the subatomic world) to calculate this new "squishiness" at the exact temperature where the soup turns back into normal particles (a phase called "freeze-out").

The Discovery:
They found that at this critical temperature, the "squishiness" of this exotic, super-hot soup is almost exactly the same as a perfect, ideal gas (like air in a balloon that doesn't interact with itself).

• The Number: They calculated the value to be 13.8 (in specific physics units).
• The Comparison: This matches perfectly with data from the ALICE experiment at the Large Hadron Collider (LHC), but only after the scientists corrected the ALICE data to include neutral particles (like neutrons) that were previously ignored.

Why This Matters

1. It Validates the Models: It proves that our theoretical models (like the Hadron Resonance Gas model) are doing a great job describing the universe at the moment of "freeze-out."
2. It's "Perfectly" Simple: Even though the matter inside a heavy-ion collision is incredibly complex and strong, at the moment it cools down, it behaves just like a simple, non-interacting gas. It's as if a chaotic mosh pit suddenly organizes itself into a perfectly polite line.
3. No Critical Point Yet: They looked at different energy levels (different "squeezes") and found that the compressibility stays steady. This suggests they haven't found a "Critical Point" (a place where the soup suddenly becomes infinitely squishy) yet, which is a major goal of these experiments.

The Takeaway

The authors successfully invented a new ruler to measure the "squishiness" of the universe's earliest moments. By switching from counting "total particles" to counting "net charge fluctuations," they avoided mathematical dead-ends.

Their conclusion is comforting: The most complex matter in the universe, at the moment it forms, behaves with the simplicity of an ideal gas. It's a reminder that sometimes, the most chaotic systems find a way to be perfectly orderly.

Technical Summary

1. Problem Statement

The isothermal compressibility ($\kappa_T$) is a fundamental thermodynamic quantity describing the response of a system's volume to changes in pressure at constant temperature. In standard statistical physics, it is defined by keeping the total particle number density ($n = N/V$) constant:
$$\kappa_T = \frac{1}{n^2} \frac{\partial^2 P}{\partial \mu^2}$$
However, applying this definition to Quantum Chromodynamics (QCD) and the Quark-Gluon Plasma (QGP) presents a fundamental theoretical challenge:

• Non-conservation of Particle Number: In QCD, the number of particles (hadrons) is not a conserved quantity due to pair creation and annihilation. The only strictly conserved quantities are the net charges associated with baryon number ($N_B$), electric charge ($N_Q$), and strangeness ($N_S$).
• Divergence at $\mu=0$: If one attempts to define $\kappa_T$ by keeping net-charge densities fixed (e.g., $N_B$) in the limit of vanishing chemical potential ($\mu \to 0$), the compressibility diverges. This is because the net charge of a system with equal numbers of particles and antiparticles is zero, making the derivative ill-defined in a non-interacting Hadron Resonance Gas (HRG) model.
• Experimental Discrepancy: Experimental determinations of $\kappa_T$ in heavy-ion collisions (HIC) often rely on charged particle yields. Previous analyses frequently neglected the contribution of neutral hadrons, leading to potential underestimations of the total compressibility.

The authors aim to resolve these issues by introducing a generalized definition of isothermal compressibility that is calculable in first-principles Lattice QCD and directly comparable to experimental data without relying on ill-defined total particle numbers.

2. Methodology

A. Generalized Definition

The authors propose a generalized isothermal compressibility, $\kappa_{T, f_Q}$, defined by keeping a specific function of the net-electric charge ($f_Q$) fixed, rather than the total particle number. The definition is:
$$\kappa_{T, f_Q, r_Q, r_S} = -\frac{1}{V} \left( \frac{\partial V}{\partial P} \right)_{T, f_Q, r_Q, r_S}$$
where $r_Q = N_Q/N_B$ and $r_S = N_S/N_B$ are fixed ratios (specifically $r_S=0$ and $r_Q=1/2$ for isospin symmetry).

They consider two specific cases for $f_Q$:

1. Fixed Net Charge ($f_Q = N_Q$): This still leads to divergence at $\mu=0$.
2. Fixed Net Charge Fluctuations ($f_Q = \sigma_Q^2$): This is the primary focus. Here, $f_Q = V T \frac{\partial^2 P}{\partial \mu_Q^2} \equiv \sigma_Q^2$.
   • This choice is motivated by the observation that in a non-interacting HRG, total particle numbers are well-approximated by net-charge fluctuations.
   • Crucially, this definition yields a finite limit as chemical potentials vanish, avoiding the divergence problem.

B. Theoretical Framework

• Lattice QCD: The authors calculate the generalized compressibility using (2+1)-flavor QCD simulations. They utilize Taylor series expansions of the pressure ($P/T^4$) up to sixth order in the baryon chemical potential ($\hat{\mu}_B$) to evaluate cumulants of conserved charges ($\hat{\chi}_{ijk}^{BQS}$).
• Hadron Resonance Gas (HRG) Model: They perform parallel calculations using a non-interacting HRG model (using the QMHRG2020 hadron list) to serve as a baseline for comparison and to interpret the QCD results near the freeze-out temperature.
• Normalization: To facilitate comparison with ideal gas behavior, the results are normalized by the pressure $P$, examining the dimensionless product $P \kappa_{T, \sigma_Q^2}$.

C. Experimental Comparison

The theoretical results are compared with data from:

• ALICE Collaboration: Data from Pb-Pb collisions at $\sqrt{s_{NN}} = 2.76$ TeV ($\hat{\mu}_B \approx 0$).
• STAR Collaboration: Data from the Beam Energy Scan (BES) at RHIC covering a range of $\hat{\mu}_B$.
• Correction: The authors re-analyze the ALICE data by replacing the number of charged hadrons ($N_{ch}$) with the total number of hadrons ($N_{tot}$), accounting for neutral particles which were previously neglected in similar analyses.

3. Key Contributions

1. Resolution of the Divergence: The paper successfully defines a thermodynamic compressibility for QCD that remains finite at $\mu=0$ by fixing charge fluctuations ($\sigma_Q^2$) instead of net charge or particle number.
2. Unified Framework: It establishes a direct link between first-principles Lattice QCD calculations and experimental heavy-ion collision observables using a consistent definition of compressibility.
3. Correction of Experimental Analysis: The authors demonstrate that neglecting neutral hadrons in experimental compressibility calculations leads to a significant underestimation (by nearly a factor of 2). They provide a "rescaled" experimental value that accounts for total hadron multiplicity.
4. Validation of Ideal Gas Behavior: The study confirms that near the pseudo-critical temperature ($T_{pc}$), the strongly interacting matter behaves thermodynamically like an ideal gas regarding compressibility.

4. Key Results

• Value at $T_{pc}$: At the pseudo-critical temperature $T_{pc,0} = 156.5(1.5)$ MeV and $\hat{\mu}_B = 0$, the generalized isothermal compressibility is calculated as:
  $$\kappa_{T, \sigma_Q^2} = 13.8(1.3) \, \text{fm}^3/\text{GeV}$$
• Normalization: When normalized by the pressure ($P$), the result is:
  $$P \kappa_{T, \sigma_Q^2} \approx 1.08(10)$$
  This is remarkably close to the ideal gas value of $1$ (where $PV=NkT \implies \kappa_T = 1/P$).
• Comparison with HRG:
  • The HRG model predicts $P \kappa_{T, \sigma_Q^2} \approx 1.72$ (or $21.93 \, \text{fm}^3/\text{GeV}$).
  • The discrepancy between QCD and HRG is attributed to the strong modification of the doubly charged $\Delta$-resonance in the QCD medium, which is treated as a stable particle in the non-interacting HRG model.
• Experimental Consistency:
  • The original ALICE result (using only charged particles) was $\kappa_T \approx 27.9 \, \text{fm}^3/\text{GeV}$.
  • The rescaled ALICE result (including neutral hadrons) is $\kappa_T \approx 15.8(1.3) \, \text{fm}^3/\text{GeV}$.
  • This rescaled experimental value is consistent with the Lattice QCD result of $13.8(1.3) \, \text{fm}^3/\text{GeV}$.
• Dependence on $\hat{\mu}_B$: The compressibility shows only a weak dependence on the baryon chemical potential along the pseudo-critical line $T_{pc}(\hat{\mu}_B)$. There is no evidence of a divergence that would signal the proximity of a critical point within the studied range.

5. Significance

• Thermodynamic Characterization: The paper confirms that the QCD matter created in heavy-ion collisions at the time of hadronization (freeze-out) behaves as an "almost perfect fluid" not just in terms of viscosity, but also in its compressibility, which remains close to that of an ideal Boltzmann gas.
• Methodological Advancement: The generalized definition provides a robust tool for future studies of the QCD phase diagram, particularly at finite baryon density where standard definitions fail.
• Bridging Theory and Experiment: By correcting the treatment of neutral hadrons in experimental analyses, the paper resolves a long-standing discrepancy between lattice calculations and experimental data, validating the HRG model's utility near $T_{pc}$ while highlighting specific medium effects (like $\Delta$-resonance modification) that distinguish QCD from a simple gas of non-interacting hadrons.
• Critical Point Search: The lack of divergence in $\kappa_T$ along the pseudo-critical line up to the studied chemical potentials suggests that the system has not yet reached the vicinity of a critical point in this regime.