Extracting freeze-out conditions in beam energy scan via functional QCD

This paper utilizes functional QCD to determine QCD freeze-out conditions by matching theoretical baryon number susceptibilities with experimental proton cumulants, successfully predicting kurtosis values that agree with data and revealing a peak structure near 5 GeV indicative of a critical end point.

The Gist

Imagine you are a detective trying to solve the mystery of how the universe was built. Specifically, you want to understand what happened a fraction of a second after the Big Bang, when the universe was a super-hot, super-dense soup of particles called quarks and gluons.

Scientists recreate this "primordial soup" by smashing heavy atoms (like gold) together at nearly the speed of light in giant machines called particle colliders. This creates a tiny, fleeting drop of this hot soup, which then cools down and freezes into ordinary matter (like protons and neutrons) that we see today.

The big question is: Exactly when and where does this "freezing" happen? And, even more excitingly, is there a hidden "Critical Point" in this process—a specific spot where the rules of physics change dramatically, like water turning into ice but with a twist?

This paper is a report from a team of theoretical physicists who used advanced math (called "Functional QCD") to answer these questions and compare their predictions with real-world experiments.

Here is the breakdown of their work using simple analogies:

1. The Problem: Comparing Apples to Oranges

The scientists have two sources of information:

• The Theory (The Apples): Their math models calculate the behavior of all baryons (a family of particles that includes protons and neutrons) in a perfect, theoretical equilibrium.
• The Experiment (The Oranges): Real-world detectors (like the STAR experiment at RHIC) can only count protons specifically. They also can't capture the perfect "equilibrium" state because the collision happens so fast it's a bit chaotic.

Usually, comparing apples to oranges is a bad idea. But the team found a clever trick. They looked at ratios (like comparing the number of apples to the number of oranges in a basket) rather than the raw numbers. They discovered that even though the "apples" (theory) and "oranges" (experiment) are different, the ratios of their fluctuations match up surprisingly well at a specific temperature and pressure.

2. The Detective Work: Finding the "Freeze-Out" Point

Think of the collision as a hot cup of coffee cooling down on a table.

• Freeze-out is the exact moment the coffee stops being a liquid and starts forming ice crystals. In the particle world, this is the moment the hot quark-gluon soup stops interacting and turns into stable protons and neutrons.
• The scientists used two specific "thermometers" (mathematical ratios of particle fluctuations) to pinpoint this moment.
  • Thermometer A is sensitive to the pressure (chemical potential).
  • Thermometer B is sensitive to the temperature.

By looking at where the "pressure reading" from the theory matched the "pressure reading" from the experiment, and doing the same for the "temperature," they found a single point where both matched. This point tells them exactly where the "freezing" happened for each collision energy.

3. The Map: Drawing the Freeze-Out Line

By doing this for many different collision energies (from very fast to slower ones), they drew a map (a curve) showing the path of the freeze-out points.

• The Result: For high-energy collisions (fast crashes), their map lines up perfectly with other methods and experimental data. It's like they found a secret path that everyone else was guessing at, but they have the GPS coordinates.
• The Surprise: When they looked at lower energies (slower crashes, around 5 GeV), their map started to diverge from the standard predictions. This is the "interesting discrepancy" mentioned in the paper.

4. The "Smoking Gun": The Kurtosis Peak

Now, for the most exciting part. The scientists predicted a specific shape for a statistical measurement called Kurtosis (which basically measures how "spiky" or "peaked" the data distribution is).

• The Analogy: Imagine a mountain range. In normal conditions, the mountains are smooth and rolling. But if you hit a "Critical Point," the mountain range might suddenly spike into a sharp, jagged peak.
• The Prediction: The team's map predicts that at a collision energy of about 5 to 6 GeV, the Kurtosis should shoot up into a sharp peak.
• Why it matters: This peak is the "smoking gun" for the Critical End Point (CEP). If the experimental data shows this peak, it proves that the QCD phase diagram has a special critical point, just like water has a critical point where liquid and gas become indistinguishable.

5. The Conclusion

The paper concludes that:

1. Their method of comparing theory and experiment works, even when comparing different types of particles.
2. They have successfully drawn a reliable map of where the "freezing" happens in the universe's early moments.
3. They predict a peak in the data around 5 GeV. This is a strong hint that the Critical End Point exists right there.

The Caveat:
The authors are humble. They admit they are still comparing "apples to oranges" (baryons vs. protons) and haven't fully accounted for the chaos of the collision (non-equilibrium effects). They say this is just the "first step" toward a complete picture. But, they have provided the first solid, equilibrium-based baseline to help experimentalists know exactly where to look for the "smoking gun" signal in their data.

In short: They built a high-tech GPS for the subatomic world, found a specific spot where the rules of physics might change, and told the experimentalists, "Go look right here; that's where the magic happens."

Technical Summary

Here is a detailed technical summary of the paper "Extracting freeze-out conditions in beam energy scan via functional QCD" by Yi Lu et al.

1. Problem Statement

The primary goal of heavy-ion collision experiments, such as the Beam Energy Scan (BES) at the Relativistic Heavy Ion Collider (RHIC), is to map the Quantum Chromodynamics (QCD) phase diagram and locate the Critical End Point (CEP). A key signature of the CEP is the non-monotonic behavior of fluctuations in conserved charges (specifically net-baryon number).

However, a significant challenge exists in comparing theoretical predictions with experimental data:

• Theory: Calculates equilibrium properties of net-baryon number fluctuations (susceptibilities, $\chi_B$).
• Experiment: Measures net-proton number fluctuations (cumulants, $C_n$) because protons are the only stable baryons detected.
• Gap: There is a lack of a quantitative, first-principles theoretical baseline that connects these two quantities across the full range of collision energies, particularly in the high baryon chemical potential ($\mu_B$) region where the CEP is expected. Previous comparisons often relied on statistical models or lattice QCD, which are limited to low $\mu_B$.

2. Methodology

The authors employ Functional QCD (fQCD) approaches (specifically Dyson-Schwinger Equations and Functional Renormalization Group methods) to compute thermodynamic quantities from first principles.

Core Strategy:

1. Freeze-out Determination: The authors define the "freeze-out" condition (the point where the system decouples from thermal equilibrium) by matching theoretical ratios of baryon number susceptibilities to experimental ratios of proton number cumulants.

   • They utilize the lowest-order ratios:
     $$\frac{\chi_2^B}{\chi_1^B} \quad \text{and} \quad \frac{\chi_3^B}{\chi_1^B}$$
   • These are equated to experimental cumulants:
     $$\frac{C_2}{C_1} \quad \text{and} \quad \frac{C_3}{C_1}$$
   • By solving these equations simultaneously for temperature ($T_f$) and baryon chemical potential ($\mu_{B,f}$) at each collision energy ($\sqrt{s_{NN}}$), they extract the freeze-out curve.

2. Handling Discrepancies:

   • The paper acknowledges that comparing baryons (theory) to protons (experiment) and ignoring non-equilibrium effects is an approximation.
   • They assume these differences approximately cancel out in the ratios of fluctuations.
   • For low collision energies ($\sqrt{s_{NN}} < 10$ GeV), where theoretical trajectories sometimes fail to intersect experimental data (indicating a mismatch between equilibrium QCD and experimental freeze-out), they use an optimization method to find the $\mu_{B,f}$ that minimizes the gap between the two extracted temperatures.

3. Parametrization:

   • The extracted points are fitted to a parametric form for the freeze-out line:
     $$T_f = T_0 \left[ 1 - \frac{\kappa_f^2}{2} \left(\frac{\mu_{B,f}}{T_0}\right)^2 - \frac{\kappa_f^4}{4} \left(\frac{\mu_{B,f}}{T_0}\right)^4 + \dots \right]$$
   • Constraints are applied based on the chiral crossover line (e.g., $T_0 \approx 157$ MeV) to ensure physical consistency.
   • A Padé fit is used to relate $\mu_{B,f}$ to collision energy $\sqrt{s_{NN}}$.

3. Key Contributions

• First Equilibrium Baseline: This work provides the first self-consistent prediction of the kurtosis (fourth-order cumulant) along the freeze-out line derived directly from functional QCD calculations at finite density.
• Bridging Theory and Experiment: It establishes a quantitative link between net-baryon susceptibility ratios (theory) and net-proton cumulants (experiment) across the entire BES energy range (200 GeV down to 3 GeV).
• Refined Freeze-out Curve: It maps out the freeze-out curve in the $(T, \mu_B)$ plane with narrowed error bars, specifically improving the precision in the critical region ($3 < \sqrt{s_{NN}} < 7.7$ GeV).
• CEP Prediction: It offers a specific prediction for the location of the kurtosis peak, which serves as a "smoking gun" signal for the CEP.

4. Key Results

• Freeze-out Curve Consistency:
  • For high collision energies ($\sqrt{s_{NN}} \gtrsim 10$ GeV), the extracted freeze-out points align well with the chiral crossover line and previous statistical model extractions.
  • For low energies ($\sqrt{s_{NN}} \lesssim 5$ GeV), the extracted points show a deviation from standard statistical models, suggesting the onset of new dynamics or phases (e.g., inhomogeneous phases or "moat regimes").
• Kurtosis Prediction:
  • The authors predict a distinct peak in the kurtosis along the freeze-out line.
  • Peak Location:
    • Upper bound of error: $\sqrt{s_{NN}} \approx 5$ GeV, with $(T_f, \mu_{B,f}) = (120, 506)$ MeV.
    • Lower bound of error: $\sqrt{s_{NN}} \approx 6$ GeV, with $(T_f, \mu_{B,f}) = (125, 432)$ MeV.
  • The peak height is lower than some previous functional QCD estimates but agrees qualitatively with the location of the CEP found in other fQCD studies.
• Agreement with Data:
  • The predicted kurtosis values show quantitative agreement with STAR BES-II data where available, despite the theoretical "apples-to-oranges" comparison (baryons vs. protons).
  • The results suggest that the lack of a critical signal in current BES-II data is consistent with the CEP being located at lower energies ($\sim 5$ GeV) than currently probed with high precision.

5. Significance

• Validation of Functional QCD: The study demonstrates that functional QCD methods are robust enough to compute thermodynamic quantities at large chemical potentials ($\mu_B/T \lesssim 4.5$) with quantitative reliability, filling a gap left by lattice QCD (which suffers from the sign problem at high $\mu_B$).
• Guidance for Future Experiments: By pinpointing the kurtosis peak around 5–6 GeV, the paper provides a specific target for future experiments (such as BES-III or the NICA program) to search for the CEP.
• Theoretical Framework: It establishes a rigorous framework for extracting freeze-out parameters without relying solely on statistical hadronization models, grounding the analysis in the fundamental thermodynamics of QCD.
• Future Directions: The authors highlight that the next critical step is to calculate proton susceptibilities directly within the functional framework to eliminate the "apples-to-oranges" approximation and to incorporate non-equilibrium effects, which are crucial near the critical point.

In summary, this paper successfully uses functional QCD to map the QCD freeze-out curve and predicts a kurtosis peak indicative of the Critical End Point at $\sqrt{s_{NN}} \approx 5$ GeV, offering a refined theoretical baseline for interpreting heavy-ion collision data.