Multiplicity dependence of thermal parameters in pp collisions at $\sqrt{s}=7$ TeV from statistical hadronization fits

This study utilizes the statistical hadronization model to analyze ALICE data from 7 TeV proton-proton collisions, revealing that while chemical freeze-out temperature and volume scale predictably with multiplicity, a persistent tension between hidden and open strangeness constraints suggests that a single global freeze-out description may be insufficient to fully characterize the strange sector in high-multiplicity events.

The Gist

Imagine you are a chef trying to figure out the secret recipe for a massive, chaotic kitchen party. In this kitchen, particles (the ingredients) are flying everywhere, colliding, and turning into new particles (the dishes).

For decades, physicists have used a "Statistical Hadronization Model" (SHM) to describe how this kitchen works in huge, explosive collisions (like smashing two heavy atoms together). They found that if you wait just long enough, the ingredients settle into a perfect, balanced soup where the temperature and volume are predictable. It's like a perfectly cooked stew.

But here's the puzzle: What happens in a tiny kitchen?

This paper investigates proton-proton collisions. Think of these as tiny, quick "snack-sized" collisions compared to the "full-course dinner" of heavy-ion collisions. The big question is: Does a tiny, chaotic snack party ever behave like a perfectly balanced stew?

Here is the breakdown of what the researchers found, using simple analogies:

1. The Setup: The "Thermal" Kitchen

The researchers used a sophisticated computer program (called Thermal-FIST) to act as a "kitchen simulator." They fed it data from the ALICE experiment at the Large Hadron Collider (LHC), where protons were smashed together at 7 TeV (a very high energy).

They looked at these collisions through different "lenses" based on how many particles came out (multiplicity).

• Low Multiplicity: A quiet, small party with few guests.
• High Multiplicity: A wild, crowded rave with thousands of guests.

They tried to fit the data to a "thermal model," which assumes the particles reached a state of equilibrium (like a pot of water reaching a boil).

2. The Three Main Ingredients (Parameters)

To describe the "soup," the model needs three main numbers:

1. Temperature ($T$): How hot the kitchen is.
2. Volume ($V$): How big the pot is.
3. Strangeness Saturation ($\gamma_S$): A measure of how well "strange" ingredients (a specific type of particle called strange quarks) are mixed in.

3. The Findings: What the Data Told Them

The Temperature: The "Thermostat" Stays the Same

Analogy: Imagine a thermostat that refuses to change, no matter how many people are in the room.
Result: The temperature stayed almost exactly the same (around 155–165 MeV) whether it was a small party or a huge rave.
Meaning: This suggests that the "boiling point" for creating particles is universal. It doesn't matter if the system is small or large; the moment particles stop interacting and freeze into their final forms, they do so at the same specific temperature. This temperature is very close to the theoretical "crossover point" where matter changes from a soup of quarks into normal particles.

The Volume: The "Pot" Grows with the Crowd

Analogy: If you invite more people to a party, you naturally need a bigger room.
Result: The size of the system ($V$) grew in a straight, predictable line as the number of particles increased.
Meaning: This makes perfect sense. More activity means a larger effective space where particles are created.

The "Strange" Ingredient: The "Mixing" Gets Better

Analogy: Imagine you are trying to mix a specific, rare spice (strangeness) into a soup. In a tiny cup (small system), it's hard to mix it evenly, so you end up with clumps of spice and areas with none. In a giant pot (large system), the stirring is so vigorous that the spice is perfectly distributed.
Result: In small collisions, the "strange" particles were suppressed (not enough of them). But as the collision got more crowded (higher multiplicity), the amount of strange particles increased, approaching a perfect mix.
Meaning: In tiny systems, the rules of conservation (you can't create a particle without its partner) make it hard to produce strange particles. But in high-multiplicity collisions, the system gets big enough that it behaves like a giant pot where the "strange" spice is fully mixed in.

4. The Big Twist: The "Tension" in the Recipe

This is the most exciting part of the paper. The researchers tried to check their recipe using two different "taste testers":

1. The $\phi$-meson: A particle with "hidden" strangeness (like a spice hidden inside a dumpling).
2. The $\Omega$-baryon: A particle with "open" strangeness (like a spice sitting right on top).

The Problem: When they used the $\phi$-meson to set the recipe, the math said one thing. When they used the $\Omega$-baryon, the math said something slightly different.
The Analogy: It's like trying to bake a cake. If you taste the frosting, you think the oven is at 350°F. But if you taste the cake inside, you think the oven is at 340°F. They don't quite agree.
The Result: There is a statistically significant "tension" (a 4-sigma difference). This means that while the whole system looks like it's in thermal equilibrium, the "strange" sector might not be fully settled yet. A single, simple recipe might not be enough to describe every single type of particle in these tiny collisions.

5. The Energy Check

The researchers also calculated the "energy per particle."
Analogy: In heavy-ion collisions, the average energy per particle is like a standard serving size (about 1 GeV).
Result: In these tiny proton collisions, as the party got bigger, the energy per particle rose from about 0.85 GeV to almost 1 GeV.
Meaning: High-multiplicity proton collisions are starting to look and act more like the massive heavy-ion collisions. They are approaching the same "thermodynamic" state.

The Bottom Line

This paper tells us that even in the smallest, most chaotic particle collisions, nature is trying to act like a calm, balanced system.

• The temperature is universal.
• The volume scales with the crowd.
• The strange particles get better mixed as the crowd grows.

However, the "tension" between different types of strange particles suggests that while these tiny systems are almost perfectly balanced, they aren't quite there yet. It's like a party that is 95% organized, but the DJ (the strange sector) is still fumbling with the playlist.

This helps physicists understand that the transition from "chaos" to "order" isn't just for giant explosions; it happens even in the smallest collisions, provided there are enough guests to make the party crowded.

Technical Summary

1. Problem Statement

The Statistical Hadronization Model (SHM) has been highly successful in describing hadron production in heavy-ion collisions, suggesting that the system reaches chemical equilibrium before hadronization. However, applying thermal models to small collision systems, such as proton-proton (pp) collisions, remains conceptually challenging.

• The Core Question: Do high-multiplicity pp collisions exhibit thermal-like behavior similar to heavy-ion collisions, and can a single set of thermal parameters consistently describe hadron yields across different multiplicity classes?
• Specific Tension: Recent ALICE data show "strangeness enhancement" in pp collisions (strange hadron yields increasing faster than non-strange ones with multiplicity). It is unclear if a unified chemical freeze-out description can simultaneously account for different sectors of the strange hadron spectrum (e.g., hidden-strangeness mesons like $\phi$ vs. multi-strange baryons like $\Omega$).
• Gap in Knowledge: Previous studies have not systematically addressed the sensitivity of thermal fits to the specific choice of hadron species included, nor have they thoroughly analyzed derived thermodynamic quantities (energy density, average energy per particle) across multiplicity classes in a single framework.

2. Methodology

The study performs a systematic thermal analysis of identified hadron yields measured by the ALICE Collaboration in pp collisions at a center-of-mass energy of $\sqrt{s} = 7$ TeV.

• Framework: The analysis utilizes the Thermal-FIST framework, a modern implementation of the Hadron Resonance Gas (HRG) model.
• Data Source: Identified hadron yields ($\pi, K, p, \Lambda, \Xi, \Omega, \phi$, etc.) across various charged-particle multiplicity classes ($\langle dN_{ch}/d\eta \rangle$).
• Fit Strategy:
  • Ensemble: Grand Canonical Ensemble with vanishing chemical potentials (appropriate for LHC energies).
  • Parameters: The global fits extract three primary parameters:
    1. Chemical freeze-out temperature ($T$).
    2. Effective system volume ($V$).
    3. Strangeness saturation parameter ($\gamma_S$), which accounts for deviations from full chemical equilibrium ($\gamma_S < 1$ implies suppression).
  • Decay Chains: Full inclusion of strong and electromagnetic resonance decays; weak decays are excluded to match experimental data treatment.
• Comparative Configurations: To test the robustness of the model, two distinct fit configurations were performed for each multiplicity class:
  1. $\phi$-constrained: Includes the $\phi$ meson (hidden strangeness) but excludes the $\Omega$ baryon.
  2. $\Omega$-constrained: Includes the $\Omega$ baryon (multi-strange) but excludes the $\phi$ meson.
• Derived Quantities: Using the extracted parameters, the study calculates energy density ($\varepsilon$) and average energy per particle ($E/N$).

3. Key Contributions

• Systematic Multiplicity Scan: Provides a comprehensive extraction of thermal parameters across the full range of charged-particle multiplicities in pp collisions using a single, consistent computational framework.
• Sensitivity Analysis: Explicitly quantifies the tension between fits constrained by hidden-strangeness mesons ($\phi$) and multi-strange baryons ($\Omega$), revealing a persistent offset in $\gamma_S$.
• Thermodynamic Scaling: Validates the scaling of derived thermodynamic quantities (specifically $E/N$) in small systems against empirical freeze-out criteria established in heavy-ion collisions.
• Quantification of Equilibration: Demonstrates that while high-multiplicity pp collisions approach thermal equilibrium, the strange sector may not be fully equilibrated under a single global freeze-out condition.

4. Key Results

A. Thermal Parameters

• Temperature ($T$): Remains approximately constant across all multiplicity classes, fluctuating between 155–165 MeV. This value is consistent with the QCD pseudo-critical temperature ($T_c$) from lattice calculations, suggesting a universal chemical freeze-out condition.
• Volume ($V$): Exhibits a clear, linear increase with charged-particle multiplicity ($V \propto \langle dN_{ch}/d\eta \rangle$). This confirms that multiplicity serves as a valid proxy for the effective system size.
• Strangeness Saturation ($\gamma_S$): Shows a strong, systematic dependence on multiplicity.
  • At low multiplicity: $\gamma_S \approx 0.7$ (significant suppression).
  • At high multiplicity: $\gamma_S \to 1.0$ (approaching full equilibration).
  • This confirms the "strangeness enhancement" phenomenon as a reduction of canonical suppression with increasing system size.

B. Derived Thermodynamic Quantities

• Energy Density ($\varepsilon$): Increases with multiplicity, reflecting larger system volumes and higher particle production.
• Average Energy per Particle ($E/N$): Increases from $\sim 0.85$ GeV at low multiplicity to $\sim 0.99$ GeV at high multiplicity. The value remains close to 1 GeV, supporting the empirical freeze-out criterion observed in heavy-ion collisions.

C. The $\Omega$-$\phi$ Tension (Critical Finding)

A systematic comparison between the two fit configurations revealed a significant discrepancy:

• $\gamma_S$ Offset: Fits constrained by the $\Omega$ baryon consistently yield lower $\gamma_S$ values than those constrained by the $\phi$ meson.
• Statistical Significance: The mean difference is $\langle \Delta \gamma_S \rangle = 0.063 \pm 0.015$, corresponding to a $\sim 4\sigma$ systematic offset.
• Implication: The $\Omega$ yield (scaling as $\gamma_S^3$) is highly sensitive to the saturation parameter. The inability of a single parameter set to simultaneously reproduce the yields of both $\phi$ and $\Omega$ suggests that a unified global freeze-out description may not fully capture the dynamics of the strange sector in small systems.

D. Particle-to-Pion Ratios

• Ratios like $\Lambda/\pi$, $\Xi/\pi$, and $\Omega/\pi$ show a hierarchy where enhancement increases with strangeness content.
• The model reproduces the qualitative trend but struggles to quantitatively match all strange species simultaneously, particularly in the multi-strange sector.

5. Significance and Conclusion

This work provides strong evidence that high-multiplicity proton-proton collisions at the LHC exhibit thermodynamic properties remarkably similar to those of heavy-ion collisions, specifically regarding the constancy of freeze-out temperature and the scaling of volume and energy density.

However, the study highlights a critical limitation: the strange sector in small systems may not be in complete chemical equilibrium describable by a single set of parameters. The persistent tension between fits constrained by hidden-strangeness mesons and multi-strange baryons suggests that the mechanisms driving strangeness production in pp collisions are more complex than a simple grand-canonical thermal model can describe.

Conclusion: While the statistical hadronization model provides a robust baseline for understanding multiplicity dependence in pp collisions, the observed discrepancies indicate that future models must account for non-equilibrium effects or distinct freeze-out conditions for different hadronic sectors in small collision systems.