Probing Strange-Quark Hadronization via (Multi-)Strange Hadron Multiplicity Distributions in Small Collision Systems with ALICE

The ALICE collaboration presents the first event-by-event measurement of probability distributions for strange-hadron multiplicities in proton-proton collisions at 5.02 TeV, providing a novel test of strangeness production mechanisms that extends beyond mean yields to probe large imbalances between strange and non-strange content.

The Gist

Imagine a high-energy particle collision as a massive, chaotic dance party where subatomic particles are the guests. Usually, when physicists study these parties, they just count the total number of guests who leave (the "average yield"). But in this new study, the ALICE collaboration decided to look at the party in a much more detailed way: they counted exactly how many guests of a specific, rare type showed up in each individual party.

Here is a breakdown of what they found, using simple analogies:

1. The Rare Guests: "Strange" Particles

In the world of particle physics, there are "strange" particles (like $K^0_S$, $\Lambda$, $\Xi$, and $\Omega$). Think of these as the rare, exotic guests at the dance party—maybe they wear a specific hat or have a unique dance move that makes them stand out.

For a long time, scientists thought these rare guests only showed up in huge, crowded parties (heavy-ion collisions). They believed that if you had a small party (like a proton-proton collision), you wouldn't see many of them. However, ALICE discovered that even in small parties, if the room gets crowded enough (high "multiplicity"), these rare guests start showing up in surprising numbers.

2. The New Method: Counting Every Single Guest

Previously, scientists would just say, "On average, we saw 5 rare guests per party." This new study is different. They used a new technique to count exactly how many rare guests were in each specific event.

• The Old Way: Looking at the average attendance.
• The New Way: Checking the guest list for every single party to see if Party A had 0 rare guests, Party B had 2, and Party C had a whopping 7.

This allowed them to see the "tails" of the distribution—those rare events where the number of strange particles was extremely high or extremely low, which average numbers usually hide.

3. The Big Discovery: It's Not Just About the Rare Guests

The researchers wanted to know: Why do these rare guests show up more when the party gets crowded? Is it just because there are more people overall, or is there a specific rule about how they mix?

To test this, they looked at ratios. Imagine comparing the number of rare guests to the number of regular guests.

• The Surprise: They found that even when they compared groups that had the same number of "strange" ingredients (a balanced ratio), the results still changed depending on how crowded the party was.
• The Analogy: Imagine you are baking cookies.
  • Scenario A: You have a limited supply of chocolate chips (strange quarks) and a huge supply of flour (light quarks).
  • Scenario B: You have a huge supply of both.
  • The study found that in a crowded party (high multiplicity), the "flour" (light quarks) becomes so abundant that it changes how the "chocolate chips" (strange quarks) get used. Instead of the chocolate chips clumping together to make a giant chocolate bar (a heavy particle like an $\Omega$), they get spread out to make many smaller chocolate chip cookies (lighter particles like $K^0_S$).

4. The "Coalescence" Picture

The paper suggests a "quark coalescence" model. Think of this like a game of musical chairs where particles are trying to pair up to form stable groups (hadrons).

• In a small, quiet room (low multiplicity), there aren't many "light" partners available. So, the "strange" particles are forced to stick together with each other to form heavy, multi-strange groups.
• In a huge, crowded room (high multiplicity), there are plenty of "light" partners. The strange particles get "distracted" and pair up with the light ones instead, creating more light particles and fewer heavy ones.

5. Testing the Models (The Simulations)

The scientists compared their real-world data against computer simulations (like video game physics engines) to see which one got the rules right.

• Model A (PYTHIA 8 Monash): This model was like a simulator that didn't understand the crowd dynamics at all. It failed to predict the results.
• Model B (EPOS LHC): This model got some things right but overestimated how many heavy particles would form.
• Model C (PYTHIA 8 with "Color Reconnection" + Ropes): This was the winner. It's like a simulator that finally understood that when the room gets crowded, the "strings" connecting the particles get tangled and rearrange (color reconnection), changing how they pair up. This model matched the real data best.

Summary

In simple terms, this paper proves that in high-energy collisions, the way particles form isn't just about how many particles are created. It's about how they mix. When the collision is very "busy," the abundance of common particles changes the rules, forcing rare particles to behave differently than they do in quiet collisions. The best computer models to explain this are the ones that account for these complex interactions between different types of particles.

Technical Summary

Technical Summary: Probing Strange-Quark Hadronization via (Multi-)Strange Hadron Multiplicity Distributions in Small Collision Systems with ALICE

Problem and Motivation
Strangeness enhancement, defined as the increased relative production of strange hadrons in heavy-ion collisions compared to proton–proton (pp) interactions, is a signature historically associated with Quark–Gluon Plasma (QGP) formation. However, recent ALICE observations at the LHC indicate that strange-hadron-to-pion yield ratios rise with increasing charged-particle multiplicity across all system sizes, from pp to p–Pb and Pb–Pb collisions, independent of collision energy. This suggests that event activity, rather than system size alone, drives strangeness production. While mean yields have been extensively studied, the underlying mechanisms governing hadronization in dense final states remain unclear. Specifically, there is a need to probe events with large imbalances between strange and non-strange content to distinguish between competing production mechanisms, such as color reconnection, rope hadronization, and modified cluster fragmentation, which have been implemented in Monte Carlo generators to explain high-density behaviors.

Methodology
To address these questions, the ALICE collaboration analyzed pp collisions at $\sqrt{s} = 5.02$ TeV using a novel event-by-event particle counting technique. The study focused on the probability distribution, $P(n_S)$, of producing a specific number ($n$) of strange particles of the same species ($K^0_S$, $\Lambda/\bar{\Lambda}$, $\Xi^-/\bar{\Xi}^+$, and $\Omega^-/\bar{\Omega}^+$) within a single event.

The analysis proceeded through the following steps:

1. Signal Extraction: Candidates were selected based on topological criteria for weak decays and particle identification. Signal and combinatorial backgrounds were separated by fitting invariant-mass distributions in bins of transverse momentum ($p_T$) and event multiplicity. Signal shapes were modeled using a double-sided Crystal Ball function, while backgrounds were parameterized with polynomials or Gaussians.
2. Event-by-Event Probability Reconstruction: Each candidate was assigned a signal probability weight ($w_S$). The algorithm evaluated all possible signal–background configurations for the $N$ candidates in an event, summing probabilities from the case where all candidates are background to the case where all are signal. This reconstructed the true strange-particle multiplicity distribution for each event.
3. Unfolding and Correction: Detector effects were corrected using a one-dimensional Bayesian unfolding procedure based on realistic simulations. Corrections were also applied for trigger inefficiencies at both the event and particle levels.
4. Yield and Ratio Calculation: From the corrected $P(n_S)$ distributions, the average production yield of $n$-particle multiplets ($\langle Y_{nS} \rangle$) was calculated using combinatorial factors. These yields were used to construct yield ratios with balanced strangeness content ($\Delta S = 0$) to isolate effects not driven solely by net strangeness differences.

Key Results

• Multiplicity Distributions: The measured $P(n_S)$ distributions for all species show that the probability of observing $n > 0$ strange hadrons increases with charged-particle multiplicity. The separation between low- and high-multiplicity classes grows with $n$, indicating a stronger-than-linear growth in the high-multiplicity tail. These distributions are well-described by Negative Binomial Distributions (NBD).
• Multiplet Yields: The average yields of multiple strange hadrons ($\langle Y_{nS} \rangle$ for $n > 1$) exhibit a more than linear increase with charged-particle multiplicity. Comparisons with Monte Carlo generators (PYTHIA 8 Monash 2013, PYTHIA 8 QCD-CR + Ropes, and EPOS LHC) reveal that standard models struggle to describe the absolute values and multiplicity evolution as the number of strange particles per event increases.
• Yield Ratios ($\Delta S = 0$):
  • Ratios such as $\langle Y_{n\Lambda} \rangle / \langle Y_{nK^0_S} \rangle$ increase with multiplicity for $n > 1$.
  • Conversely, ratios involving multi-strange baryons in the numerator and single-strange mesons in the denominator (e.g., $\langle Y_{1\Omega^-} \rangle / \langle Y_{3K^0_S} \rangle$) decrease with multiplicity, despite the heavier mass and higher baryon number of the numerator.
  • These trends suggest that neither mass nor baryon number alone dictates the evolution; rather, the availability of light quarks relative to strange quarks plays a critical role. A naive quark-coalescence picture explains this: in low-multiplicity events, light quark scarcity favors the formation of multi-strange baryons (like $\Omega^-$), while high-multiplicity events, with abundant light quarks, favor the formation of kaons.
• Model Performance: Among the tested models, PYTHIA 8 QCD-CR + Ropes provides the best overall agreement, reproducing both the absolute scale and the multiplicity evolution of the $\Delta S = 0$ ratios. This highlights the importance of color reconnection mechanisms involving light quarks in capturing the observed hadronization trends.

Significance and Claims
This work extends the study of strangeness production beyond average yields by employing event-by-event counting to probe the full probability distribution of strange particle production. The authors claim that these measurements provide a new, stringent test bench for hadronization models, particularly for events with extreme imbalances between strange and non-strange content.

The paper asserts that the observed multiplicity dependence of balanced yield ratios contains non-trivial information that cannot be attributed solely to strangeness enhancement. Instead, the data suggests that the interplay between strange- and light-quark availability is a key factor in the hadronization process. While the production rate of strange quarks remains an open issue (as indicated by the failure of models to fully describe multiplet yields), the color-reconnection mechanism appears to successfully capture the main trends in how produced strange quarks hadronize. These results establish a new experimental approach that significantly constrains the dynamics of strange-particle production in small collision systems.