⚛️ phenomenology

System-size dependence of strangeness production from p+p to Pb+Pb: quantitative tests of the $K^+/\pi^+$ horn

This study quantitatively evaluates various theoretical frameworks against $K^+/\pi^+$ data across the full NA61/SHINE system-size ladder, revealing that no single model explains the entire dataset but that a combination of core-corona and canonical statistical models successfully describes the observed step-like enhancement of strangeness production and the "horn" structure in heavy-ion collisions.

Original authors: Neeraj, Amal Sarkar

Published 2026-08-31

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Neeraj, Amal Sarkar

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

When physicists smash atomic nuclei together at nearly the speed of light, they are trying to recreate a state of matter that existed just moments after the Big Bang. In these collisions, the protons and neutrons that usually make up atomic nuclei melt apart, freeing their internal building blocks—quarks and gluons—to roam freely in a hot, dense soup known as a quark-gluon plasma. A key signature of this transformation is how the collision produces particles containing strange quarks. In the early universe, this strange matter was abundant, and scientists have long suspected that if they could create enough of it in the lab, they would see a sudden, dramatic spike in its production as the system grows large enough to sustain the plasma. This spike, often called a "horn" because of its shape on a graph, has been seen in the heaviest collisions, but for decades, researchers have struggled to understand exactly when and why it appears, or if it is a true sign of a new phase of matter or just a complex side effect of the collision's energy.

A new study by researchers at the Indian Institute of Technology Mandi takes a fresh look at this puzzle by changing the variable they examine. Instead of just increasing the energy of the collisions, they varied the size of the atomic nuclei being smashed together. They analyzed data from a massive range of experiments, starting with simple collisions between single protons and moving up through collisions involving light nuclei like beryllium, intermediate ones like argon and scandium, and finally the heaviest targets like lead and gold. By comparing these different system sizes at the same collision energies, the team could isolate the effect of size from the effect of energy, asking a simple but profound question: does the strange matter production jump suddenly once a certain size is reached, or does it grow smoothly and gradually?

The researchers tested this question against five different theoretical models, each representing a different idea about how the collision works. Some models assumed the collision was purely a game of billiard balls, where particles bounce off each other without ever forming a new state of matter. Others assumed the collision created a hot, dense core surrounded by a thinner, less active outer layer. Still others relied on complex statistical rules about how particles behave in a confined space. The team ran their data through all these models to see which one could accurately predict the number of strange particles produced in every single collision size they studied.

The results were clear and surprising: no single model could explain the entire range of data. The models that worked well for the smallest collisions failed miserably for the largest ones, and vice versa. The most striking discovery was that the production of strange particles does not grow smoothly as the nuclei get bigger. Instead, the data shows a sharp, step-like jump in production between the collisions of beryllium nuclei and those of argon and scandium. Below this size, the production of strange particles looks just like what happens in simple proton collisions. But once the system reaches the size of argon and scandium, the production suddenly leaps up to match the levels seen in the heaviest lead collisions. This jump happens at a specific size, suggesting that a new physical mechanism turns on only when the colliding system is large enough to form a dense, thermalized core.

This finding helps rule out several popular explanations. The study shows that models based purely on the movement of hadrons (particles like protons and neutrons) cannot generate this sudden jump, even if they include the creation of a quark-gluon plasma. The data also contradicts the idea that the strange matter production is a smooth, continuous process that scales gently with the number of particles involved. Instead, the evidence points to a threshold effect, where the system must reach a critical size before the strange particles can be produced efficiently. The researchers calculated that this critical size corresponds to a system with an effective mass number of about eighteen, a value that sits right between the light beryllium and the intermediate argon-scandium collisions.

To ensure this jump was real and not an artifact of the measurement, the team also looked at a specific ratio involving negative kaons and pions. This comparison allowed them to cancel out the effects of the collision's energy and focus purely on the system size. The analysis confirmed that the jump is a genuine physical phenomenon, with a statistical certainty that makes it highly unlikely to be a random fluctuation. The study concludes that the formation of this dense core is a geometric process, dependent on the size of the nuclei, while the full chemical balance of the strange particles takes even larger systems to achieve.

The work leaves a clear path forward for future experiments. The researchers predict that if new data is collected from collisions involving xenon and lanthanum nuclei, the results will clearly distinguish between the model that sees a gradual transition and the one that sees a sharp onset. Until that data arrives, the current findings stand as a robust demonstration that the birth of this exotic state of matter is not a smooth slide but a sudden step, triggered the moment the colliding system becomes large enough to sustain it.

Technical Summary: System-size dependence of strangeness production from p+p to Pb+Pb

Problem Statement
The non-monotonic structure in the $K^+/\pi^+$ excitation function, known as the "horn," observed in central Pb+Pb collisions at the CERN SPS is a well-established feature often interpreted as a signature of the transition to a deconfined Quark-Gluon Plasma (QGP). However, the physical origin of this structure remains ambiguous, as multiple distinct theoretical frameworks (grand-canonical statistical hadronization, strangeness over-saturation, and microscopic transport models) can reproduce the energy dependence without invoking a sharp phase transition. While the energy-axis dependence is degenerate, the system-size dependence ( $A$ -dependence) at fixed collision energy ( $\sqrt{s_{NN}}$ ) offers a potential discriminator. Previous studies have noted a step-like enhancement in strangeness production between light systems (Be+Be) and intermediate systems (Ar+Sc), termed the "onset of fireball," but existing model comparisons along the system-size axis have been fragmentary, lacking a uniform quantitative assessment across the full NA61/SHINE system-size ladder (p+p, Be+Be, Ar+Sc, Xe+La, Pb+Pb) and STAR Au+Au data.

Methodology
The authors perform a comprehensive, quantitative comparison of five distinct model classes against experimental $K^+/\pi^+$ data spanning the complete system-size ladder. The models represent mutually exclusive physical hypotheses regarding the origin of system-size dependence:

SMASH-3.3: A purely hadronic transport model (null hypothesis) with no partonic degrees of freedom.
PHSD: A partonic transport model incorporating Dynamical Quasi-Particle Model (DQPM) degrees of freedom to test if deconfined matter alone suffices.
Glauber Core-Corona (CC): A geometric superposition model where a thermalized "core" (grand-canonical HRG, $\gamma_s=1$ ) coexists with a peripheral "corona" (modeled by SMASH output). The core fraction $f_{core}$ is determined by a Glauber collision-number threshold.
Canonical Ensemble (CE) Suppression: A statistical model testing if exact strangeness conservation in a finite correlation volume ( $V_0$ ) accounts for suppression in small systems.
Statistical Model of the Early Stage (SMES): A two-phase statistical model predicting a sharp onset of deconfinement.

The analysis utilizes a uniform $\chi^2$ definition, incorporating statistical and systematic uncertainties from published spectra (NA49, NA61/SHINE, STAR BES-I). The study distinguishes between final published data and preliminary Xe+La data (digitized from Ref. [19]), using the latter only for visual comparison and specific falsifiable predictions, excluding them from model calibration and ranking logic. Key metrics include $\chi^2/\text{ndf}$ , out-of-sample predictions, and leave-one-out cross-validation.

Key Contributions and Results

Model Hierarchy and Failure of Single Frameworks: No single model describes the full range from p+p to Pb+Pb.
- Transport Models: Neither SMASH nor PHSD generates the "horn" structure. Both exhibit monotonic system-size scaling, failing to reproduce the step-like enhancement observed between Be+Be and Ar+Sc. PHSD yields higher strangeness than SMASH but remains featureless in its excitation function.
- Core-Corona (CC): This model provides the best description for intermediate systems (Ar+Sc) and heavy systems (Pb+Pb) when calibrated on Pb+Pb data ( $\chi^2/\text{ndf} = 1.69$ ). It predicts a transition in the core fraction from $f_{core} \approx 0.22$ (Be+Be) to $f_{core} \approx 0.58$ (Ar+Sc). However, it fails for Au+Au (collider geometry) and overcorrects Be+Be.
- SMES: The SMES model with canonical strangeness conservation (SMES-CE) excels for heavy systems (Au+Au, Pb+Pb) with zero fitted parameters but significantly overpredicts intermediate systems like Ar+Sc.
- Canonical Ensemble: A global fit for the correlation volume $V_0$ yields a value ( $0.20 \text{ fm}^3$ ) far below the literature expectation ( $\sim 7 \text{ fm}^3$ ) and fails to simultaneously describe small and large systems, indicating that canonical suppression alone is insufficient.
Quantification of the "Onset of Fireball": The data exhibit a distinct step-like enhancement in $K^+/\pi^+$ between Be+Be and Ar+Sc. Within the CC framework, this is quantified as a transition in the core fraction. A sigmoid fit to the core fraction yields a critical effective mass number $A_c \approx 18$ , characterizing the midpoint of this step. This value is model-dependent (sensitive to the corona baseline, e.g., SMASH vs. PHSD) and relies on preliminary Xe+La data for constraint; without Xe+La, the uncertainty diverges.
Two-Step Strangeness Saturation ( $\gamma_s$ ): The extracted strangeness saturation factor $\gamma_s$ reveals a two-step structure rather than a single sigmoid transition.
1. Geometric Onset: A rapid rise from p+p ( $\gamma_s \approx 0.22$ ) to Ar+Sc ( $\gamma_s \approx 0.65$ ), corresponding to the formation of a dense geometric core.
2. Thermodynamic Equilibration: A further increase from Ar+Sc to Pb+Pb ( $\gamma_s \approx 0.95$ ), corresponding to the approach to full chemical equilibrium.
  This structure suggests that geometric core formation and thermodynamic equilibration are governed by distinct mechanisms at different system-size scales.
Double Ratio Analysis: The construction of the double ratio $R = (K^+/\pi^+)/(K^-/\pi^-)$ cancels the system-independent baryon chemical potential ( $\mu_B$ ) contribution. $R$ exhibits a step between light and heavy systems at $\sqrt{s_{NN}} \approx 7.6$ GeV, rising from $\approx 1.8$ to $\approx 2.8$ . This step has a global significance of $3.7\sigma$ across five SPS energies, providing evidence for genuine strangeness enhancement beyond the pair-production baseline and $\mu_B$ effects.
Falsifiable Prediction: The paper identifies a divergence of up to 27% between the CC and SMES predictions for the Xe+La system at the horn energy. Finalized Xe+La spectra are expected to discriminate between gradual equilibration (CC) and a sharp onset (SMES).

Significance
The paper establishes that the system-size axis breaks the degeneracy of mechanisms explaining the $K^+/\pi^+$ horn. The primary finding is that the mechanism of strangeness production changes with system size: no single framework captures the full ladder. The "onset of fireball" is identified as a step-like enhancement distinct from the energy-axis onset of deconfinement. The work provides a quantitative benchmark for future models, demonstrating that while geometric core-corona superposition explains intermediate systems and statistical models explain heavy systems, neither alone suffices. The identification of a two-step structure in $\gamma_s$ and the $3.7\sigma$ step in the double ratio offer robust, model-independent evidence for a system-size-driven enhancement mechanism that operates preferentially in systems with effective mass numbers $A_{eff} \gtrsim 40$ . The study concludes that the critical system size $A_c \approx 18$ is a characterization of the gap between hadronic transport and data within the CC framework, rather than a universal geometric constant, and highlights the necessity of final Xe+La data and light-ion runs (e.g., O+O) to further constrain these onset mechanisms.

More like this