Extracting the speed of sound of QCD from transverse momentum fluctuations

This paper extracts the speed of sound in the quark-gluon plasma from ATLAS transverse momentum fluctuation data in ultra-central Pb+Pb collisions by correcting for detection biases and hadronization noise, yielding a result of $c_s/c=0.496\pm 0.008$ that aligns perfectly with first-principles lattice QCD calculations.

The Gist

Imagine you are a chef trying to figure out the "stiffness" of a giant, invisible soup made of the universe's most fundamental ingredients. This soup is called Quark-Gluon Plasma (QGP), and it only exists for a tiny fraction of a second when you smash two heavy atoms (like lead) together at nearly the speed of light.

The scientists in this paper are like detectives trying to measure the speed of sound in this cosmic soup. Why? Because the speed of sound tells us how "stiff" or "squishy" the material is. If sound travels fast, the material is stiff (like a steel rod); if it's slow, it's squishy (like Jell-O).

Here is how they did it, explained in simple terms:

1. The Experiment: Smashing Atoms

Think of the Large Hadron Collider (LHC) as a giant particle cannon. They fire two lead nuclei at each other. When they collide head-on (called "ultra-central" collisions), they create a tiny, super-hot drop of this plasma soup.

The scientists looked at the debris flying out of this explosion. Specifically, they measured two things:

• How many particles came out? (The "multiplicity").
• How fast were they moving sideways? (The "transverse momentum").

2. The Theory: The "Crowded Room" Analogy

The researchers had a clever idea based on a simple physical principle: Density and Temperature.

Imagine a crowded room where people are dancing.

• If you pack more people into the same size room, the room gets hotter and more chaotic. The people (particles) start moving faster because they are bumping into each other more.
• In the atom collision, if you get more particles out of the same size "soup drop," it means the soup was denser and hotter.

The paper argues that in these perfect, head-on collisions, the size of the soup drop is roughly the same every time. So, if you see more particles, the soup must have been hotter. And the hotter the soup, the faster sound travels through it.

By measuring how much the speed of the particles increases as the number of particles increases, they can calculate the speed of sound.

3. The Problem: The "Blind" Detector

Here is the catch. The detector (ATLAS) isn't perfect. It's like a security camera that is too slow to catch people walking slowly.

• The detector misses the "slow" particles (low momentum).
• This creates a bias: The data looks different than reality because the slowest particles are invisible.

If you just looked at the raw data, your calculation of the speed of sound would be wrong. It's like trying to guess the average speed of traffic in a city by only counting the race cars and ignoring the sedans and trucks.

4. The Solution: "De-blurring" the Picture

The authors didn't just ignore the missing data; they built a mathematical "lens" to fix it. They used three main tricks:

• The "Variance" Trick: They didn't just look at the average speed of particles; they looked at how much the speeds fluctuated from event to event. Even though the detector misses slow particles, the pattern of the missing particles tells them exactly how to correct the average.
• The "Noise" Filter: When atoms break apart into particles (a process called hadronization), it's a bit like popping popcorn. There's random noise. The scientists developed a method to mathematically "de-blur" this noise, separating the random popping from the actual physics of the soup.
• The "Impact Parameter" Correction: Sometimes, even in "head-on" collisions, the atoms don't hit perfectly dead center. The scientists had to account for these tiny wobbles to ensure they were only looking at the perfect collisions.

5. The Result: A Perfect Match

After all these corrections, they calculated the speed of sound ($c_s$) in the plasma.

• Their Result: They found the speed of sound is about 0.5 times the speed of light.
• The "Wow" Factor: This number matches perfectly with calculations made by supercomputers using Lattice QCD (a method that solves the laws of physics from the ground up, without any experimental data).

Why Does This Matter?

It's like finding a fingerprint at a crime scene and having it match the suspect's DNA perfectly.

• The DNA: The theoretical predictions from supercomputers (Lattice QCD).
• The Fingerprint: The real-world data from the ATLAS detector.

The fact that they match so perfectly tells us two huge things:

1. Our understanding of the universe is correct: We really do understand how the strong nuclear force works, even at temperatures 100,000 times hotter than the center of the sun.
2. The soup is "stiff": The Quark-Gluon Plasma is not a gas; it behaves like a nearly perfect liquid with a specific "stiffness" that we can now measure.

In summary: The scientists took messy, incomplete data from a particle collider, used clever math to "fix" the missing pieces and remove the noise, and discovered that the speed of sound in the early universe's soup matches our best theoretical predictions perfectly. It's a victory for both experimental physics and theoretical physics.

Technical Summary

Here is a detailed technical summary of the paper "Extracting the speed of sound of QCD from transverse momentum fluctuations".

1. Problem Statement

The primary objective of this research is to extract a quantitative value for the speed of sound ($c_s$) in the Quark-Gluon Plasma (QGP) created in ultra-relativistic heavy-ion collisions.

• Theoretical Context: In ultra-central Pb+Pb collisions (impact parameter $b \approx 0$), the mean transverse momentum per particle, $\langle [p_T] \rangle$, is predicted to increase with charged particle multiplicity ($N_{ch}$) due to thermalization and density fluctuations. The slope of this increase is directly related to $c_s^2$.
• The Challenge: Previous attempts to extract $c_s$ from this relationship suffered from significant systematic biases:
  1. Detector Acceptance: Experiments like ATLAS cannot detect particles with low transverse momentum ($p_T < 0.5$ GeV/c). This missing fraction distorts the relationship between the measured mean $\langle [p_T] \rangle$ and the true fluid properties.
  2. Statistical Fluctuations: The hadronization process introduces Poissonian statistical noise (fluctuations in the number of detected particles relative to the hydrodynamic entropy) that "blurs" the underlying hydrodynamic signal.
  3. Initial State Correlations: The extraction relies on the assumption that initial entropy ($S$) and transverse radius ($R^2$) fluctuations are uncorrelated at $b=0$. If they are correlated, the simple linear relation breaks down.

2. Methodology

The authors employ a multi-step analytical framework combining hydrodynamic response theory with rigorous statistical unfolding to correct for experimental biases.

A. Hydrodynamic Response Model

• The authors use a simplified effective description where final-state observables (multiplicity $N$ and mean $[p_T]$) are functions of initial-state quantities: total entropy $S$ and mean squared transverse radius $R^2$.
• They derive a linear relationship between relative fluctuations:
  $$\frac{\delta [p_T]}{\langle [p_T] \rangle} \approx c_s^2 \left( \frac{\delta S}{\langle S \rangle} - \frac{3}{2} \frac{\delta R^2}{\langle R^2 \rangle} \right)$$
• Crucially, they assume that at $b=0$, $S$ and $R^2$ are uncorrelated, allowing $c_s^2$ to be extracted from the slope of $\langle [p_T] \rangle$ vs. $N$.

B. Correction for Detector Acceptance (Low-$p_T$ Cuts)

• The authors introduce a correction factor using the observable $v_0(p_T)$, which describes how the normalized $p_T$ spectrum changes with temperature.
• They define two dimensionless correction factors, $C_A$ and $D_A$, derived from integrating $v_0(p_T)$ over the detector acceptance range ($p_T > 0.5$ GeV/c).
• These factors modify the relationship between measured fluctuations ($\delta N_A, \delta [p_T]_A$) and initial fluctuations, introducing a coupling between multiplicity and radius fluctuations that must be disentangled using the variance of $[p_T]$.

C. Unfolding Hadronization Noise (Deblurring)

• The transition from continuous hydrodynamic entropy to discrete hadron counts involves Poisson fluctuations.
• The authors treat the measured multiplicity $N_{ch}$ as a noisy version of the true hydrodynamic multiplicity $N_A$.
• Using Bayesian inference, they reconstruct the distribution of $N_A$ given a fixed $N_{ch}$.
• They apply a Next-to-Leading Order (NLO) approximation to "deblur" the moments of $[p_T]$, effectively removing the statistical noise introduced by the finite number of particles and reconstruction efficiency ($\epsilon \approx 0.73$). This correction is found to be substantial (renormalizing $c_s^2$ by $\approx 27\%$).

D. Impact Parameter Fluctuations

• The analysis accounts for the fact that "ultra-central" events still have a small spread in impact parameter $b$.
• They reconstruct the distribution of $b$ at fixed $N_{ch}$ and average the hydrodynamic predictions over this distribution to match experimental data.

3. Key Contributions

1. Systematic Bias Correction: The paper provides the first comprehensive method to correct for the combined effects of low-$p_T$ acceptance cuts and hadronization statistical noise in the extraction of $c_s$.
2. Utilization of Variance: The authors demonstrate that extracting $c_s$ from the mean $\langle [p_T] \rangle$ alone is insufficient when acceptance cuts are present; the variance of $[p_T]$ is required to disentangle the effects of entropy fluctuations from radius fluctuations.
3. Bayesian Deblurring: They introduce a rigorous statistical procedure to unfold Poisson fluctuations, showing that neglecting this leads to a $\sim 19-27\%$ underestimation of $c_s^2$.
4. Reverse Engineering Initial States: By fitting the data, they infer the relative standard deviations of the initial entropy ($S$) and radius ($R^2$), finding them to be smaller than some standard initial-state model predictions.

4. Results

• Speed of Sound: After applying all corrections to ATLAS data (Pb+Pb at $\sqrt{s_{NN}} = 5.02$ TeV), the authors extract:
  $$c_s/c = 0.496 \pm 0.008$$
  at an effective temperature of $T_{eff} = 221 \pm 13$ MeV.
• Agreement with Lattice QCD: This value is in perfect agreement with first-principles calculations from Lattice QCD (HotQCD and BMW collaborations), which predict a peak in $c_s$ near the crossover temperature.
• Initial State Fluctuations: The analysis reveals that the relative standard deviation of the initial entropy is $\approx 2.77\%$, and the radius is $\approx 3.09\%$. The correlation between $S$ and $R^2$ at $b=0$ is confirmed to be negligible, validating the core assumption of the extraction method.
• Sensitivity Analysis: The authors show that the result is robust against variations in the fit range and initial-state model parameters, with the largest uncertainty arising from the fit range and acceptance correction factors.

5. Significance

• Validation of QGP Properties: The result provides a robust, model-independent confirmation that the QGP behaves as a nearly perfect fluid with a speed of sound consistent with the equation of state derived from Lattice QCD.
• Methodological Benchmark: The paper establishes a new standard for analyzing fluctuations in heavy-ion collisions. It highlights that "self-correlations" and detector acceptance effects must be treated with high precision to avoid significant systematic errors.
• Future Applications: The framework developed here can be applied to other fluctuation observables (such as skewness) and can be adapted for data from other experiments (ALICE, CMS) which have different $p_T$ acceptance cuts, potentially reducing uncertainties further.

In summary, this work successfully bridges the gap between complex experimental data and fundamental QCD theory by rigorously correcting for experimental artifacts, yielding a precise measurement of the speed of sound in the quark-gluon plasma.