Baryon fluctuation signatures of the onset of deconfinement

The paper predicts that anomalous proton number fluctuations observed in heavy-ion collisions around a center-of-mass energy of 10 GeV serve as a signature for the onset of deconfinement, offering a unified explanation for experimental data from both RHIC and the CERN SPS.

The Gist

Imagine you are a chef trying to figure out exactly when your soup changes from a thick, chunky stew into a smooth, clear broth. You can't see the transformation happening inside the pot, so instead, you look at how the ingredients are bouncing around.

This paper is about a similar "cooking" process, but the pot is a particle accelerator, the ingredients are subatomic particles, and the "soup" is the matter that existed just after the Big Bang.

Here is the story of the paper, broken down into simple concepts:

1. The Big Goal: Finding the "Switch"

Scientists smash heavy atoms (like gold or lead) together at incredibly high speeds. They want to see if they can turn normal matter (where protons and neutrons are stuck together like Lego bricks) into a new state of matter called Quark-Gluon Plasma (QGP).

Think of normal matter as a crowded dance floor where everyone is holding hands in tight groups (protons and neutrons).
Think of the QGP as a mosh pit where everyone has let go of their hands and is running around freely.

The big question is: At what exact speed (energy) does the dance floor turn into a mosh pit? The authors believe this "switch" happens around a specific energy level (about 10 GeV).

2. The Clue: Counting the "Bouncers"

How do we know the switch happened? We can't see the quarks directly. Instead, the scientists look at fluctuations.

Imagine you are counting how many people leave the dance floor every minute.

• In the "Stew" (Normal Matter): People leave in groups of three (because protons are made of three quarks). If you count the groups, the numbers are very predictable. It's like counting full boxes of eggs; you know exactly how many eggs are in there.
• In the "Broth" (QGP): The boxes have melted. Now, individual eggs (quarks) are running around. Because quarks are "fractional" (they are only 1/3 of a proton), the way they fluctuate changes completely.

The authors predict that as you increase the collision energy, the way these numbers wiggle and change will suddenly shift. It's like hearing the music change from a steady drumbeat to a chaotic jazz solo.

3. The "Net-Baryon" Mystery

The scientists focus on a specific count: the difference between protons and anti-protons (called "net-baryons").

• Low Energy (The Stew): The math says the fluctuations should be very steady. If you count the groups, the ratio of "wiggle" to "average" is always 1.
• High Energy (The Broth): Once the matter melts into quarks, the math changes. Because quarks carry only 1/3 of the charge, the "wiggle" becomes much smaller relative to the average. The ratio drops significantly.

The paper argues that if you plot this ratio against the collision energy, you should see a sharp dip or a sudden change right around that 10 GeV mark. This dip is the "signature" of the deconfinement (the moment the Lego bricks fall apart).

4. The Real-World Messiness

The authors admit that real experiments aren't as clean as their math models.

• The "Window" Problem: Detectors can't see every particle leaving the collision; they only see a slice of the action (like looking through a keyhole). The authors had to create a mathematical "correction" to guess what the whole picture looks like based on that small slice.
• The "Proton" Filter: The detectors actually count protons, not the raw quarks. Since protons are made of three quarks, the signal gets "diluted." It's like trying to hear a whisper in a noisy room; the signal is there, but it's quieter.

5. The Conclusion: A Match Made in Heaven

When the authors compared their "soup recipe" predictions with real data from experiments at CERN and the RHIC (the Relativistic Heavy Ion Collider), they found something exciting:

The data actually shows that weird dip!

The experimental data shows a strange, non-smooth change in how proton numbers fluctuate right around the energy level the authors predicted (8–12 GeV). This suggests that the "switch" from normal matter to the Quark-Gluon Plasma is indeed happening right there.

The Takeaway

This paper is like a detective story. The scientists built a theory about how the "suspect" (the phase transition) would behave. They looked at the "crime scene" (the experimental data) and found fingerprints that matched their theory perfectly.

They aren't just saying "we think this happened." They are saying, "Look at this specific pattern in the data; it's the smoking gun that proves matter changes its nature at this specific energy level." This helps us understand the fundamental rules of the universe and how the very first moments after the Big Bang worked.

Technical Summary

1. Problem Statement

The primary goal of high-energy nucleus-nucleus collision experiments is to map the QCD phase diagram, specifically identifying the transition from confined hadronic matter to deconfined Quark-Gluon Plasma (QGP). While the existence of QGP is established, the precise nature of the transition (crossover vs. first-order phase transition) and the location of the conjectured Critical Point (CP) remain open questions.

Recent experimental data from the STAR collaboration (RHIC) and NA49 (SPS) regarding proton number fluctuations in central collisions show an anomalous, non-monotonic dependence on collision energy ($\sqrt{s_{NN}}$). The authors hypothesize that this anomaly is not necessarily a signature of the Critical Point, but rather a direct consequence of the onset of deconfinement occurring around $\sqrt{s_{NN}} \approx 8\text{--}12$ GeV. The paper aims to theoretically predict the collision-energy dependence of baryon and proton number cumulants across this transition to explain the experimental observations.

2. Methodology

The authors employ a comparative modeling approach using two distinct statistical ensembles to represent the two phases of matter:

• Confined Phase (Hadron Resonance Gas - HRG): Modeled as an ideal Boltzmann gas of baryons and antibaryons. This is valid for $\sqrt{s_{NN}} \leq 8$ GeV.
• Deconfined Phase (Quark-Gluon Plasma - QGP): Modeled as an ideal gas of massless $u, d, s$ quarks and antiquarks. This is valid for $\sqrt{s_{NN}} \geq 12$ GeV.

The calculation proceeds in three sequential steps to bridge the gap between theoretical models and experimental observables:

1. Grand Canonical Ensemble (GCE) Calculation:

   • HRG: Assumes Poisson distributions for baryon ($N_b$) and antibaryon ($N_{\bar{b}}$) numbers. The net-baryon number ($B = N_b - N_{\bar{b}}$) follows a Skellam distribution.
   • QGP: Uses the ideal-gas equation of state with massless quarks. Cumulants ($\kappa_n$) are derived from derivatives of the pressure $P$ with respect to the baryon chemical potential $\mu_B$.
   • Key Theoretical Insight: In the confined phase, baryon number is carried by integer units (baryons). In the deconfined phase, it is carried by fractional units (quarks, $1/3$). This leads to a fundamental suppression of fluctuations in the QGP phase.

2. Finite Acceptance and Conservation Laws:

   • Experiments measure particles within a limited momentum acceptance, not the full system.
   • The authors apply the Sub-ensemble Acceptance Method (SAM) to account for global net-baryon number conservation. This introduces a correction factor based on $\alpha$, the fraction of baryons detected in the acceptance.
   • Formulas (Eqs. 13–15) relate GCE cumulants to acceptance-limited cumulants ($BA$).

3. Hadronization to Protons:

   • Experiments measure protons, not net-baryons.
   • The authors assume that for a fixed net-baryon number, the proton number fluctuates via binomial thinning with probability $\beta$ (the fraction of baryons that are protons).
   • Standard relations (Appendix A) are used to convert net-baryon cumulants to proton number cumulants ($pA$).

3. Key Contributions

• Theoretical Mechanism for Anomaly: The paper provides a quantitative mechanism explaining the "anomalous" energy dependence of proton fluctuations. It attributes the change in cumulant ratios to the shift in effective degrees of freedom from integer-charged baryons to fractionally-charged quarks.
• Quantitative Prediction of Cumulant Ratios: The authors derive specific predictions for the ratios of the first four cumulants ($\kappa_2/\kappa_1$, $\kappa_3/\kappa_2$, $\kappa_4/\kappa_2$) for both net-baryons and protons across the transition region.
• Integration of Conservation Laws: Unlike simple GCE predictions, this work rigorously incorporates global conservation laws and experimental acceptance effects, which are crucial for comparing theory with data from STAR and NA49.
• Unified Interpretation: The model offers a consistent interpretation of data spanning the CERN SPS, BNL RHIC, and GSI SIS energies without invoking the Critical Point.

4. Results

The study produces the following key results, visualized in Figures 1–3 and Appendix B:

• Qualitative Change at $\sqrt{s_{NN}} \approx 8\text{--}12$ GeV:

  • HRG Phase ($\leq 8$ GeV): Cumulant ratios for net-baryons are close to unity ($\kappa_n/\kappa_{n-1} \approx 1$), characteristic of a Poissonian distribution.
  • QGP Phase ($\geq 12$ GeV): Cumulant ratios drop significantly due to the fractional charge of quarks. For example, $\kappa_2/\kappa_1$ drops to $\approx 1/3$ in the ideal QGP limit.
  • Transition: A rapid, non-monotonic change occurs in the transition region.

• Suppression in Proton Observables:

  • The magnitude of the change is most pronounced for the total net-baryon number ($B$).
  • The effect is dampened for the net-baryon number in acceptance ($BA$) due to conservation laws (factor $1-\alpha$).
  • The effect is further dampened for proton number ($pA$) due to binomial thinning (factor $\beta$).
  • Despite this damping, the model predicts a distinct qualitative change (a "dip" or non-monotonic behavior) in the proton cumulant ratios around 10 GeV.

• Comparison with STAR Data:

  • When the model predictions are converted to factorial cumulants (to match STAR's reported $C_n/C_1$) and plotted against experimental data (Fig. B.4), the model successfully reproduces the non-monotonic behavior observed in the data.
  • The model captures the trend where ratios decrease as energy increases from the SPS range into the RHIC Beam Energy Scan range.

5. Significance and Implications

• Alternative to Critical Point: The results suggest that the observed anomalies in proton fluctuations may be a robust signature of the onset of deconfinement rather than the Critical Point. This challenges the exclusive focus on the CP and highlights the phase transition itself as a source of fluctuation signatures.
• Consistency Across Experiments: The model provides a unified framework that explains data from GSI SIS, CERN SPS, and RHIC, suggesting a continuous evolution of the system's degrees of freedom.
• Future Directions: The authors note that while the qualitative agreement is strong, quantitative precision requires:
  • Using experimentally determined, energy-dependent values for acceptance ($\alpha$) and proton fraction ($\beta$).
  • Better modeling of nuclear cluster production.
  • Investigating the sensitivity to the QGP equation of state and the freeze-out conditions.

In conclusion, the paper establishes that the transition from confined to deconfined matter fundamentally alters the statistical properties of baryon number fluctuations, providing a natural and consistent explanation for recent experimental anomalies in heavy-ion collisions.