Chemical potential differentials in the QCD phase diagram from heavy-ion isobar collisions

This paper utilizes Bayesian thermal analysis of hadron yields from STAR Ru+Ru and Zr+Zr isobar collisions to precisely extract chemical potential differentials in the QCD phase diagram, thereby validating these collisions as a high-precision probe for four-dimensional QCD thermodynamics against lattice-QCD and Chiral Mean Field model predictions.

The Gist

Imagine the universe as a giant, chaotic kitchen where particles are constantly cooking, mixing, and changing states. Sometimes, under extreme heat and pressure, these particles melt into a soupy state called "quark-gluon plasma." Physicists want to understand exactly how this soup behaves, but it's incredibly hard to taste the soup directly because it changes so fast.

This paper is like a team of master chefs and detectives trying to figure out the exact recipe of that soup by looking at the leftovers (the particles that survive the explosion) after a collision. Here is the story of what they did, explained simply:

1. The Experiment: A Tale of Two Twins

The scientists used a giant particle accelerator (the RHIC) to smash heavy atoms together. Usually, smashing two identical atoms is like hitting a drum with two identical hammers. But this time, they used two very specific "twins":

• Twin A (Ruthenium): Has 44 protons and 52 neutrons.
• Twin B (Zirconium): Has 40 protons and 56 neutrons.

They have the same total weight (96 parts), but Twin A is slightly "more positive" (more protons) than Twin B. It's like comparing two identical backpacks where one has a few extra heavy coins in the pocket. The scientists wanted to see how the "soup" inside the collision reacted to this tiny difference in the coins.

2. The Problem: The Noise in the Signal

When they smashed these twins apart, they looked at the particles that flew out. They wanted to measure the "chemical potential," which is a fancy physics word for the pressure or drive of different types of charges (like baryons, electric charge, and strangeness) inside the soup.

The problem? When they measured the twins separately, the difference was so small that the "static noise" of the experiment hid the answer. It was like trying to hear a whisper in a hurricane. The uncertainty was too high to say for sure if the twins produced different results.

3. The Solution: The "Double-Check" Trick

To fix this, the team used a clever statistical trick called a Bayesian analysis. Instead of measuring the twins separately, they looked at the difference between them directly.

Think of it like this: If you want to know the exact weight difference between two nearly identical apples, you don't weigh them on two different scales (which might be slightly off). You put them on a balance scale together. The errors cancel out, and you see the tiny difference clearly.

By comparing the "net charge" (the total positive minus negative particles) of the Ruthenium crash against the Zirconium crash, they could isolate the tiny shift caused by the extra protons. This reduced the "noise" and let them see the signal clearly.

4. The Findings: Mapping the Terrain

The results showed that even a tiny change in the number of protons (about 9% difference) caused a measurable shift in the "chemical pressure" of the soup.

• The Map: They created a 4-dimensional map of the QCD phase diagram (a map of how matter behaves under extreme conditions).
• The Arrow: They found that changing the proton count pushes the system in a specific direction on this map. It's like pushing a boat slightly off course; the water reacts in a predictable way.
• The Ratios: They calculated how the "baryon pressure" changes relative to the "charge pressure" and "strangeness pressure." It's like figuring out that if you add a little more sugar, the cake rises a specific amount relative to how much it spreads.

5. Checking Against Theory: The Recipe Books

The scientists then compared their experimental "leftovers" with two different theoretical "recipe books" (models) that try to predict how this soup should behave:

1. Lattice QCD (BQS): A method based on supercomputer calculations from first principles.
2. Chiral Mean Field (mCMF): An effective model that treats particles like interacting waves.

The Verdict:

• Both recipe books got the direction of the shift right (they agreed on which way the arrow pointed).
• The "Lattice" book was better at predicting how the "baryon" pressure changed relative to "charge."
• The "Mean Field" book was better at predicting how "strangeness" changed relative to "charge."
• Neither book was perfect; there were still small discrepancies, suggesting there are still some missing ingredients (like specific types of particles) in the theoretical recipes.

Why This Matters

This paper is a breakthrough because it proves that by using these "isobar" twins (atoms with the same weight but different proton counts), scientists can now measure the properties of the quark-gluon plasma with much higher precision than before.

It's like upgrading from a blurry photo to a high-definition image. They have successfully mapped out how the fundamental forces of nature respond to tiny changes in the composition of matter, bridging the gap between what we see in particle colliders and what we know about the extreme conditions inside neutron stars.

In short: They used a clever comparison trick to turn a blurry whisper into a clear shout, revealing exactly how the universe's most extreme matter reacts to a tiny change in its recipe.

Technical Summary

Technical Summary: Chemical Potential Differentials in the QCD Phase Diagram from Heavy-Ion Isobar Collisions

Problem Statement
Understanding the thermodynamics of Quantum Chromodynamics (QCD) matter under extreme conditions requires characterizing the chemical potentials associated with conserved charges: baryon number ($\mu_B$), electric charge ($\mu_Q$), and strangeness ($\mu_S$). While thermal models successfully extract the freeze-out temperature ($T_{\text{Chem}}$) and individual chemical potentials from hadron yields, previous analyses of heavy-ion collisions often imposed simplifying assumptions, such as $\mu_Q = 0$ or a fixed charge fraction. Consequently, the interdependence of these potentials, their mutual response (susceptibilities), and their differentials (e.g., $d\mu_B/d\mu_Q$) in dynamically evolving systems have remained largely unexplored. Furthermore, theoretical calculations, such as Lattice-QCD, are currently limited to specific regions of the phase diagram (vanishing $\mu_B$ or high $T$), necessitating effective models or Taylor expansions to bridge gaps.

Methodology
This work presents a Bayesian thermal analysis of hadron yields from the STAR experiment's isobar collisions: Ruthenium-96 ($^{96}_{44}\text{Ru}$) and Zirconium-96 ($^{96}_{40}\text{Zr}$). These systems share the same mass number ($A=96$) but differ in proton/neutron content, resulting in distinct charge fractions ($Y_Q \approx 0.46$ for Ru and $Y_Q \approx 0.42$ for Zr).

1. Bayesian Framework: The authors utilize the statistical model THERMUS within a Bayesian framework. A Gaussian-process emulator serves as a surrogate model to compare experimental data (particle yields, yield ratios, and net-charge differences) against theoretical predictions.
2. Handling Systematic Uncertainties: To address the issue where independent analyses of Ru and Zr yield differences consistent with zero due to large uncertainties, the study employs a method proposed in Ref. [25]. This involves analyzing the net-charge difference ($\Delta Q$) between the two isobar systems using double ratios of particle yields ($R^2_i$). This approach cancels out most systematic uncertainties.
3. Baryon Stopping: The analysis incorporates baryon stopping effects, derived from the STAR measurement of $\langle B \rangle / \Delta Q$, by introducing a factor $\alpha = 1/1.84$ in the charge-fraction constraint.
4. Theoretical Comparison: The extracted experimental differentials are compared against two theoretical frameworks:
   • BQS: A Taylor-expanded Lattice-QCD equation of state around $\mu_B = \mu_Q = \mu_S = 0$, reliable up to $\mu_B/T \approx 2.5$.
   • mCMF: An improved Chiral Mean Field effective model including baryon octet and decuplet, capable of describing the entire multidimensional phase diagram.
5. Statistical Propagation: Uncertainties are propagated via Monte Carlo sampling. For ratios of small differences (where denominators fluctuate near zero), Fieller's theorem is applied to obtain statistically meaningful confidence intervals.

Key Results

• Extraction of Differentials: The simultaneous Bayesian fit successfully extracts the chemical potential differences ($\Delta \mu_i = \mu_i^{\text{Zr}} - \mu_i^{\text{Ru}}$) and their ratios. The results show that the small deviation in charge fraction ($Y_Q$) between the isobars generates measurable shifts in the freeze-out conditions.
• Experimental Values: For the 0–10% centrality class at $\sqrt{s_{NN}} = 200$ GeV, the analysis yields:
  • $T_{\text{Chem}} \approx 158$ MeV (common to both systems).
  • $\Delta \mu_B \approx 0.14$ MeV, $\Delta \mu_S \approx 0.12$ MeV, and $\Delta \mu_Q \approx -0.16$ MeV.
  • Key ratios include $\Delta \mu_B / \Delta \mu_Q \approx -0.97$ and $\Delta \mu_S / \Delta \mu_Q \approx -0.77$.
• Theoretical Agreement:
  • Both BQS and mCMF reproduce the sign and approximate magnitude of the experimental differentials and ratios.
  • The BQS expansion shows the closest quantitative agreement for the baryon-to-charge ratio ($\Delta \mu_B / \Delta \mu_Q$).
  • The mCMF model better reproduces the strangeness-to-charge ratio ($\Delta \mu_S / \Delta \mu_Q$).
  • The largest deviations between theory and experiment occur in the baryon-to-strangeness channel ($\Delta \mu_B / \Delta \mu_S$).
• Phase Diagram Trajectories: The study maps the trajectory of the freeze-out points in the multidimensional QCD phase diagram. The directionality of the shift from Ru to Zr reflects the correlated response of conserved charges. Theoretical derivatives ($d\mu_i/d\mu_j$) computed at $\mu_B \approx 20$ MeV differ significantly from those at $\mu_B \to 0$, highlighting enhanced sensitivity near vanishing net baryon density.

Significance and Claims
The paper claims to provide the first extraction of chemical potentials and their differentials in relativistic heavy-ion collisions using a full Bayesian analysis that simultaneously fits hadron yields and isobar yield differences. By isolating the net charge through isobar comparisons, the authors substantially reduce systematic uncertainties, enabling a precise determination of the multidimensional QCD thermodynamics at $\sqrt{s_{NN}} = 200$ GeV.

The work establishes isobar collisions as a precision probe for the QCD phase diagram, demonstrating that controlled variations in isospin can reveal the response of the system to changes in conserved charges. Furthermore, it creates a framework connecting heavy-ion phenomenology with first-principles Lattice-QCD and effective models, while also linking these findings to constraints relevant for neutron star matter (via the mCMF model's applicability to high $\mu_B$/low $T$ regions). The authors note that while current models capture the general trends, discrepancies in specific ratios (particularly $\Delta \mu_B / \Delta \mu_S$) may indicate missing resonances in effective descriptions, suggesting a path for future refinements.