A study of charged-particle multiplicity distribution in high energy p-O collisions

This study investigates charged-particle multiplicity distributions in high-energy p-O collisions across various energies and pseudorapidity intervals using Pythia and $k_T$-factorization approaches, revealing that the choice of initial oxygen nucleus configuration (specifically $\alpha$-cluster versus Woods-Saxon models) and the theoretical formalism significantly influence the resulting multiplicity profiles, particularly at large multiplicities and higher pseudorapidities.

The Gist

The Big Picture: Smashing Tiny Oranges

Imagine you are a physicist trying to understand how the universe works at its most fundamental level. To do this, you take two tiny particles and smash them together at nearly the speed of light. Usually, scientists smash heavy "lead" atoms together because they are big and messy, creating a super-hot soup of energy called the Quark-Gluon Plasma (think of it as a cosmic smoothie made of the universe's building blocks).

But recently, scientists started smashing Oxygen atoms (which are much smaller) into protons. Why? Because Oxygen is like a "mini-lead." It's small enough to test if the rules of physics change when you shrink the size of the collision.

This paper asks a simple question: Does the way the Oxygen atom is built inside matter when we smash it?

The Two Ways to Build an Oxygen Atom

The researchers realized that we don't actually know exactly how the 16 particles inside an Oxygen atom are arranged. They tested two different "blueprints" for the atom:

1. The "Smooth Ball" (Woods-Saxon Model): Imagine the Oxygen atom is like a fluffy cotton ball or a cloud. The particles are spread out evenly and smoothly throughout the sphere. There are no distinct lumps; it's just a continuous cloud of stuff.
2. The "Lego Tetrahedron" (Alpha-Cluster Model): Imagine the Oxygen atom is actually made of four tiny, hard Lego bricks (called alpha particles) glued together in a pyramid shape (a tetrahedron). The particles are clumped tightly into these four distinct groups, with empty space in between.

The Experiment: The Grand Smash

The researchers used two different computer programs to simulate smashing a proton into an Oxygen atom using these two blueprints.

• Program A (Pythia/Angantyr): Think of this as a "classic" simulation. It treats the collision like a game of billiards where particles bounce off each other, break apart, and reassemble into new particles. It's very detailed but relies on standard rules of particle physics.
• Program B (kT-factorization): Think of this as a "modern, high-speed" simulation. It focuses on the chaotic, turbulent nature of the collision, assuming that at these high speeds, the particles behave more like a fluid wave than individual billiard balls.

They ran these simulations at different energies (how hard they smashed) and looked at the results in different "viewing angles" (pseudorapidity).

The Results: What They Found

1. The Shape of the Atom Changes the Outcome
When they smashed the "Smooth Ball" Oxygen vs. the "Lego Tetrahedron" Oxygen, the results were surprisingly different.

• The Analogy: Imagine throwing a dart at a fluffy cloud versus throwing a dart at a pyramid of four hard rocks.
  • If you hit the cloud, the dart passes through gently, creating a moderate amount of debris.
  • If you hit the pyramid, and you hit one of the hard rocks, it shatters violently, creating a massive explosion of debris.
• The Finding: The "Lego" model produced significantly more particles (a bigger explosion) in the "tail" of the results—meaning in the rare, extreme collisions where the proton hit a dense cluster directly. The "Smooth Ball" model was much more consistent and less explosive. This tells us that how the nucleus is organized inside matters a lot for the final explosion.

2. The Two Computer Programs Disagree
The two simulation programs (Pythia and kT-factorization) gave very different answers.

• Pythia showed a "bumpy" result with a peak and a valley (like a rollercoaster).
• kT-factorization showed a much smoother curve.
• The Takeaway: This suggests that our current understanding of the physics is incomplete. We don't know yet which computer program is telling the truth. We need real experimental data from the Large Hadron Collider (LHC) to see which model matches reality.

3. The "Universal Rule" Holds Up
The researchers checked something called KNO Scaling.

• The Analogy: Imagine you have a bucket of popcorn. If you pop it at low heat, you get a few kernels. If you pop it at high heat, you get a lot. KNO scaling is a mathematical rule that says: If you divide the number of popped kernels by the average number, the shape of the distribution looks the same, no matter how hot the heat is.
• The Finding: Even though the Oxygen atom is small and the collisions are complex, this "Universal Rule" still worked! The shape of the particle explosion remained consistent across different energy levels. This is a comforting sign that the fundamental laws of physics are stable, even in these tiny systems.

4. Two Types of Explosions
Finally, they tried to fit the data using a "Double Negative Binomial Distribution" (a fancy math term).

• The Analogy: Think of the debris from the crash as coming from two different sources:
  1. Soft Process: A gentle puff of smoke (low energy, common).
  2. Semi-Hard Process: A violent shattering of glass (high energy, rare).
• The Finding: The data fit perfectly when they assumed there were two types of events happening at once. This confirms that particle production isn't just one random thing; it's a mix of gentle interactions and violent shattering.

The Bottom Line

This paper is like a detective story trying to figure out the internal structure of an Oxygen atom by watching what happens when it gets smashed.

• Main Discovery: The internal "architecture" of the Oxygen atom (whether it's a smooth cloud or clumpy Lego bricks) changes the outcome of the crash significantly.
• Mystery: Our two best computer simulations disagree on exactly what happens, meaning we need more real-world data to solve the puzzle.
• Good News: Despite the chaos, the universe still follows some beautiful, universal patterns (KNO scaling), and the explosions are a mix of gentle puffs and violent shatters.

This research helps us understand not just how Oxygen behaves, but also how high-energy cosmic rays (particles from space) hit our atmosphere, helping us build better models for everything from particle physics to weather prediction!

Technical Summary

1. Problem Statement

The study addresses the need to understand particle production mechanisms in small collision systems, specifically proton-Oxygen (p-O) collisions at the Large Hadron Collider (LHC). While heavy-ion collisions (e.g., Pb-Pb) are well-studied for Quark-Gluon Plasma (QGP) formation, smaller systems like p-O offer a unique testing ground for:

• Nuclear Geometry Effects: Investigating how the internal structure of the Oxygen nucleus (specifically, whether it behaves as a smooth continuous distribution or a clustered structure of $\alpha$-particles) influences final-state observables.
• Theoretical Frameworks: Comparing two distinct approaches to modeling high-energy hadronic interactions: the standard Pythia (Angantyr) generator (based on collinear factorization and Monte Carlo event generation) and the $k_T$-factorization approach (based on the Color Glass Condensate (CGC) effective field theory, which accounts for non-linear gluon saturation effects).
• Multiplicity Dynamics: Characterizing the charged-particle multiplicity distributions across various center-of-mass energies ($\sqrt{s}$) and pseudorapidity ($\eta$) intervals to test scaling laws (KNO) and phenomenological models (Double Negative Binomial Distribution).

2. Methodology

The authors employed a comparative simulation approach using two primary theoretical frameworks and two distinct nuclear geometry models.

A. Theoretical Frameworks

1. Pythia (Angantyr) Model:
   • Utilizes the Angantyr extension of Pythia 8.313 to handle nuclear collisions without QGP formation.
   • Relies on collinear factorization for perturbative processes (hard scattering, Initial/Final State Radiation) and the Lund string fragmentation model for hadronization.
   • Uses the Glauber formalism with the Woods-Saxon distribution to sample nucleon positions and calculate nucleon-nucleon interactions.
   • Incorporates Color Reconnection and Multi-Parton Interactions (MPI).
2. $k_T$-factorization Approach:
   • Based on the Color Glass Condensate (CGC) framework, which addresses high parton density and saturation at small-$x$.
   • Uses Unintegrated Gluon Distributions (UGDs) derived from the running coupling Balitsky-Kovchegov (rcBK) equation.
   • Implements event-by-event fluctuations of the saturation scale ($Q_s$) using a log-normal distribution to account for geometric fluctuations.
   • Convolved with the KKP Fragmentation Function to obtain hadronic level results.

B. Nuclear Geometry Models

The study explicitly compares two descriptions of the Oxygen-16 ($^{16}\text{O}$) nucleus:

1. Woods-Saxon (WS) Model: A continuous, smooth nuclear charge density distribution (standard approximation for heavy nuclei).
2. $\alpha$-Cluster (AC) Model: A discrete model where the 16 nucleons are arranged as four $\alpha$-clusters (tetrahedral geometry) based on Extended Quantum Molecular Dynamics (EQMD). The density within clusters is described by a three-parameter Fermi (3pF) distribution to account for non-Gaussian central curvature.

C. Simulation Parameters

• Collisions: p-O at $\sqrt{s} = 2.36, 5.02, 7.0, 9.62, 13.0$ TeV.
• Pseudorapidity Intervals: $|\eta| < 0.5, 1.0, 2.0, 3.0$.
• Analysis: Calculation of charged-particle multiplicity ($N_{\text{charged}}$) distributions, average multiplicity ($\langle N_{\text{charged}} \rangle$), KNO scaling verification, and fitting with a Double Negative Binomial Distribution (DNBD).

3. Key Contributions

• Systematic Comparison of Nuclear Geometries: The paper provides a rigorous quantitative comparison of how $\alpha$-clustering versus a smooth Woods-Saxon distribution affects charged-particle multiplicity in p-O collisions, a topic previously underexplored compared to p-Pb or O-O collisions.
• Cross-Framework Validation: It offers a direct comparison between the widely used Pythia/Angantyr generator and the more theoretical $k_T$-factorization/CGC approach for p-O systems, highlighting where they converge and diverge.
• Scaling and Parameterization: The study validates the applicability of KNO scaling and the Double NBD model to p-O collisions across a wide energy range, elucidating the contributions of soft and semi-hard processes.

4. Key Results

A. Impact of Nuclear Geometry (AC vs. WS)

• Multiplicity Differences: The geometric description of the nucleus significantly impacts the multiplicity distribution.
  • At low multiplicities, the AC and WS models yield similar results (statistical uncertainty dominates).
  • At high multiplicities (the tail of the distribution) and higher pseudorapidities, the models diverge significantly.
  • The AC model (clustered) generally produces a different tail behavior compared to the WS model because the clustered structure creates high local density regions, leading to different initial-state fluctuations. The inversion of the curves (where AC surpasses or falls below WS) shifts depending on the pseudorapidity window.

B. Pythia vs. $k_T$-factorization

• Structural Differences: The two frameworks produce distinct multiplicity distributions.
  • Pythia exhibits a characteristic "peak-valley" structure at low $N_{\text{charged}}$, which is absent in the $k_T$-factorization results.
  • At intermediate multiplicities, the models agree better at central rapidity ($|\eta| < 0.5$) but diverge at forward rapidity ($|\eta| < 3.0$).
  • Tail Divergence: Significant divergence is observed in the high-multiplicity tail for both rapidity windows, suggesting that the underlying physics mechanisms (perturbative vs. saturation-driven) treat extreme fluctuations differently.

C. Average Multiplicity and Energy Evolution

• $\langle N_{\text{charged}} \rangle$ increases with both center-of-mass energy and pseudorapidity for both models.
• The energy growth is non-abrupt and consistent with expectations.
• The comparison allows for an indirect probe of longitudinal particle production and Bjorken-$x$ evolution.

D. KNO Scaling and DNBD Fitting

• KNO Scaling: The study confirms that KNO scaling holds for p-O collisions. When distributions are scaled by the average multiplicity ($z = N/\langle N \rangle$), data from different energies collapse onto a single universal curve for both Pythia and $k_T$-factorization, indicating no violation of universality in the observed range.
• Double NBD Fit: The charged-particle multiplicity distributions are well-described by a Double Negative Binomial Distribution (DNBD).
  • The first component (NBD 1) fits the low-multiplicity region (Poissonian-like, soft processes).
  • The second component (NBD 2) fits the high-multiplicity tail (semi-hard processes).
  • This confirms that the observed multiplicity is a superposition of distinct event classes (soft and semi-hard) rather than a single homogeneous process.

5. Significance

• Nuclear Structure Probing: The results demonstrate that global inclusive observables like charged-particle multiplicity can serve as a tool to distinguish between different nuclear geometry models (clustered vs. smooth) in small systems. This is crucial for interpreting future LHC Run 3 and Run 4 data involving Oxygen beams.
• Model Constraints: The observed discrepancies between Pythia and $k_T$-factorization, particularly in the tails of the distribution and at forward rapidities, provide critical constraints for theoretical models. Experimental data will be required to determine which framework (or combination thereof) best describes the underlying physics of p-O collisions.
• Cosmic Ray Physics: Improved understanding of p-O interactions at high energies aids in refining hadronic interaction models used in atmospheric shower simulations for cosmic ray physics, where interactions with atmospheric nuclei occur at energies beyond fixed-target capabilities.
• Universality: The confirmation of KNO scaling and the success of the DNBD fit reinforce the idea that despite different initial conditions and theoretical approaches, the statistical properties of particle production in high-energy collisions follow universal patterns driven by soft and semi-hard dynamics.

In conclusion, this work establishes that the initial geometric configuration of the Oxygen nucleus and the choice of theoretical formalism (collinear vs. $k_T$-factorization) have profound and measurable effects on charged-particle multiplicity, particularly in the high-multiplicity regime and forward rapidity, offering new avenues for probing QCD dynamics in small collision systems.