More uses for Thermal Models

This paper proposes parameter-free combinations of particle and anti-particle yields to test thermal models and extract freeze-out chemical potentials ($\mu_B, \mu_S, \mu_Q$) from RHIC BES data, successfully validating these methods against known results while extending predictions to unmeasured (anti-)nuclei yields and providing updated parametrizations for the energy dependence of freeze-out parameters.

The Gist

Imagine you are a detective trying to figure out what happened at a massive, chaotic party that happened a billionth of a second ago. This party was a collision between two heavy atoms (like gold nuclei) smashed together at nearly the speed of light.

When these atoms crash, they create a tiny, super-hot soup of fundamental particles called Quark-Gluon Plasma. As this soup cools down, it freezes into a "soup" of ordinary particles like protons, neutrons, and strange particles. This moment when the particles stop changing into one another is called "Chemical Freeze-out."

The goal of this paper is to figure out the exact "temperature" and "pressure" (chemical potentials) of this soup at the moment it froze, without needing to do incredibly complex, computer-heavy math.

Here is the breakdown of their clever detective work, using simple analogies:

1. The Old Way vs. The New Way

The Old Way (The Global Fit):
Traditionally, scientists tried to figure out the party conditions by measuring the total number of every single type of guest (protons, pions, strange particles, etc.) and running a massive computer simulation to see what conditions would produce that exact mix. It's like trying to guess the temperature of a room by counting every single dust mote, which is hard because you need to know exactly how many dust mites were there to begin with.

The New Way (The Ratio Trick):
The authors of this paper say, "Why count the total number? Let's just look at the ratio of guests to their 'anti-guests'."

• Imagine for every Proton (a guest), there is an Anti-Proton (a guest wearing a hat backwards).
• In this hot soup, the number of "hat-backwards" guests depends entirely on how much "baryon pressure" (chemical potential) was in the room.
• By looking at the ratio of Protons to Anti-Protons, the messy details (like the total size of the room or the mass of the particles) cancel out. It's like comparing the number of left shoes to right shoes; you don't need to know how many people are in the room to know if the ratio is 1:1.

2. The "Double Ratio" Test (The Reality Check)

The authors found a brilliant shortcut. They noticed that if you take the ratio of Lambda particles to Anti-Lambdas, and divide that by the ratio of Protons to Anti-Protons, you get a very specific number.

• The Analogy: Imagine you have a recipe for a cake. If you double the sugar, the cake gets sweeter. If you double the flour, it gets denser. But if you look at the ratio of sugar to flour in two different cakes, and that ratio stays the same, you know the bakers used the same recipe.
• The Result: The paper shows that for many different particles (Protons, Lambdas, Cascades), these "double ratios" are all equal to each other. This proves that the "soup" really did reach a state of thermal equilibrium (it was a well-mixed, hot soup) and that the thermal model is working correctly. It's a quick "sanity check" before doing the heavy lifting.

3. The "Cosh" Calculator (Cracking the Code)

Once they established the ratios work, they used a specific mathematical formula (involving a function called "hyperbolic cosine," or cosh) to reverse-engineer the conditions of the soup.

• They measured the ratio of Protons (to get the Baryon pressure).
• They measured the ratio of Lambdas (to get the Strangeness pressure).
• They measured the ratio of Cascades (to get a mix of both).
• By plugging these three simple numbers into their "decoder ring" (the equations), they could instantly calculate the Temperature (T), Baryon Chemical Potential ($\mu_B$), and Strangeness Chemical Potential ($\mu_S$).

It's like having three different locks on a door. Instead of trying to pick all three locks, they found a master key (the ratios) that opens them all at once.

4. Predicting the Invisible (The Crystal Ball)

The most exciting part of the paper is that they used this method to predict things that haven't been measured yet.

• The Omega Particle: They used their calculated "soup conditions" to predict how many Omega particles (a very heavy, rare particle) should exist. When they checked the actual data, the prediction was spot on. This proved their "decoder ring" was accurate.
• Anti-Nuclei: They went a step further. They predicted how many Anti-Deuterons and Anti-Tritons (anti-nuclei made of anti-protons and anti-neutrons) should exist at low energies where we haven't measured them yet.
  • Why is this cool? It's like predicting how many snowflakes will fall in a blizzard you haven't seen yet, just by knowing the temperature and humidity of the air. They predicted the yield of Anti-Deuterons at a specific low energy (7.7 GeV) where no one had measured it before.

5. The Big Picture: Why Do We Care?

The authors updated a "map" of the conditions in these collisions.

• High Energy: At very high energies (like the LHC), the soup is hot and the "baryon pressure" is low. It's like a hot, empty room.
• Low Energy: At lower energies (like the RHIC BES experiments), the soup is cooler but has high "baryon pressure." It's like a dense, crowded room.

By mapping out exactly how the temperature and pressure change as we slow down the collisions, they are helping physicists understand the Phase Diagram of QCD (the rulebook of how matter behaves).

The Ultimate Goal:
This isn't just about particle colliders. The conditions in these low-energy, high-density collisions are very similar to the core of Neutron Stars. By understanding the "soup" in the lab, we can better understand the physics of the densest objects in the universe.

Summary

In short, this paper is a guide on how to stop counting every single particle in a cosmic crash and start using simple ratios to instantly decode the temperature and pressure of the resulting soup. It's a smarter, faster, and more elegant way to understand the fundamental building blocks of our universe, and it even lets us predict the existence of particles we haven't found yet.

Technical Summary

1. Problem Statement

In relativistic heavy-ion collisions, the extraction of chemical freeze-out parameters—specifically temperature ($T$) and chemical potentials for baryon number ($\mu_B$), strangeness ($\mu_S$), and charge ($\mu_Q$)—is conventionally performed via global thermal fits to measured hadron yields. This standard approach relies on statistical hadronization or hadron resonance gas models and is sensitive to systematic uncertainties arising from:

• The specific list of hadrons included.
• Decay treatments and feed-down corrections.
• Ensemble choices (canonical vs. grand-canonical).
• Overall volume factors and particle masses.

The authors identify a need for parameter-free methods to test the validity of thermal models and to extract chemical potential ratios ($\mu/T$) directly from experimental data with reduced systematic dependencies.

2. Methodology

The paper proposes a novel approach based on constructing specific combinations of particle and anti-particle yields that cancel out mass-dependent terms and volume factors.

A. Theoretical Framework:

• Particle/Anti-particle Ratios: In the Grand-Canonical ensemble, the ratio of particle ($n_h$) to anti-particle ($n_{\bar{h}}$) densities is given by:
  $$\frac{n_h}{n_{\bar{h}}} = \exp\left[\frac{2(B_h\mu_B + S_h\mu_S + Q_h\mu_Q)}{T}\right]$$
  This eliminates mass dependence ($m_i$) and degeneracy factors ($g_i$).
• Double Ratios: The authors define double ratios (e.g., $\frac{n_\Lambda/n_{\bar{\Lambda}}}{n_p/n_{\bar{p}}}$) to isolate specific chemical potentials. Under the approximation that $\mu_Q/T$ is small, these double ratios become equal to simple functions of $\mu_S/T$ (e.g., $\exp[4\mu_S/T]$). These serve as direct tests of local thermal equilibrium.
• The $R_h$ Variable: A new variable is introduced, defined as the ratio of the arithmetic mean to the geometric mean of particle and anti-particle yields:
  $$R_h \equiv \frac{n_h + n_{\bar{h}}}{2\sqrt{n_h n_{\bar{h}}}} = \cosh\left(\frac{B_h\mu_B + S_h\mu_S + Q_h\mu_Q}{T}\right)$$
  Crucially, $R_h$ depends only on quantum numbers and chemical potentials, independent of mass and degeneracy.

B. Inversion Strategy:
By measuring $R_p$, $R_\Lambda$, and $R_\Xi$ (for protons, lambdas, and cascades), the authors derive analytical expressions to invert the system and extract $\mu_B/T$, $\mu_S/T$, and $\mu_Q/T$ directly:

• $\mu_B/T$ and $\mu_Q/T$ are derived from combinations of $\cosh^{-1}(R_p)$, $\cosh^{-1}(R_\Lambda)$, and $\cosh^{-1}(R_\Xi)$.
• $\mu_S/T$ is isolated using the difference between $\cosh^{-1}(R_p)$ and $\cosh^{-1}(R_\Xi)$.

C. Data Sources:
The method is applied to experimental data from:

• RHIC BES Phase-1: Au+Au collisions at $\sqrt{s_{NN}} = 7.7$–39 GeV (STAR collaboration).
• AGS and SPS: Lower energy data from previous experiments.
• Validation: The extracted parameters are used to predict $R_\Omega$ (Omega baryon ratio) and light (anti-)nuclei yields ($R_d$, $R_t$) to verify consistency.

3. Key Contributions

1. Parameter-Free Extraction: Development of a method to extract $\mu_B/T$, $\mu_S/T$, and $\mu_Q/T$ directly from particle/anti-particle ratios without global fitting or reliance on specific hadron lists.
2. Thermalization Tests: Introduction of double-ratio tests (Eq. 3 and Eq. 4) that verify local thermal equilibrium. If the model holds, specific double ratios of yields (and higher-order cumulants) must be equal.
3. First Parametrization of $\mu_S$: The paper presents the first functional parametrization of the strangeness chemical potential ($\mu_S$) as a function of collision energy ($\sqrt{s_{NN}}$).
4. Prediction of Unmeasured Yields: The method is extended to predict yields of light anti-nuclei (anti-deuterons, anti-tritons) in energy regimes where they have not yet been measured.

4. Key Results

• Validation of Thermal Equilibrium:
  • Double ratios constructed from $p, \Lambda, \Xi,$ and $K$ yields are consistent with each other at high energies ($\sqrt{s_{NN}} \ge 20$ GeV), confirming the thermal model's validity.
  • Deviations are observed at lower energies, suggesting potential challenges to complete thermalization in the BES-I low-energy regime.
• Extraction of Chemical Potentials:
  • The extracted $\mu_B/T$ increases with centrality (number of participants, $\langle N_{part} \rangle$), reflecting higher baryon density.
  • $\mu_S/T$ is independent of centrality but decreases significantly with increasing beam energy, consistent with the OZI rule (strange/anti-strange pair creation).
  • $\mu_Q/T$ is found to be small and negative, consistent with the system being nearly isoscalar.
• Agreement with STAR Fits:
  • The values of $\mu_B/T$ and $\mu_S/T$ extracted via this direct method show excellent agreement with the published global thermal fit results from the STAR collaboration (using THERMUS).
• Predictive Success:
  • $R_\Omega$ Verification: The method successfully predicts $R_\Omega$ (using inputs from $p, \Lambda, \Xi$), validating the approach for multi-strange baryons.
  • Light Nuclei: The model successfully predicts anti-deuteron ($\bar{d}$) and anti-triton ($\bar{t}$) yields. Specifically, it provides a genuine prediction for $\bar{d}$ yields at $\sqrt{s_{NN}} = 7.7$ GeV, where experimental data is currently unavailable.
• New Parametrizations:
  • Updated functional forms for $T(\sqrt{s_{NN}})$, $\mu_B(\sqrt{s_{NN}})$, and the first-ever $\mu_S(\sqrt{s_{NN}})$ are provided.
  • Extrapolations are provided for energies where data is sparse or missing (e.g., AGS 4.85 GeV, SPS 6.27 GeV, and RHIC 14.5 GeV).

5. Significance

• Reduced Systematics: By utilizing ratios that cancel mass and volume dependencies, this method reduces systematic uncertainties inherent in global fits, offering a more robust check of the thermal model.
• Low-Energy Physics: As the Beam Energy Scan (BES) moves to lower energies and higher statistics (BES-II), these double-ratio tests become critical for probing the limits of thermalization and the QCD phase diagram.
• Neutron Star Connection: Precise determination of $\mu_B, \mu_S,$ and $\mu_Q$ is essential for extrapolating heavy-ion collision data to the conditions found in neutron star interiors. The new parametrizations, particularly for $\mu_S$, provide better constraints for these astrophysical models.
• Predictive Power: The ability to predict unmeasured anti-nuclei yields allows experimentalists to plan future measurements and test the coalescence vs. statistical hadronization mechanisms for light nuclei formation.

In summary, the paper demonstrates that specific combinations of particle/anti-particle yields offer a powerful, computationally simple, and systematic-reduced tool for extracting freeze-out parameters and testing the fundamental assumptions of thermal models in high-energy nuclear physics.