Forward hadron production in pp collisions at LHC energies from an event generator based on the color glass condensate framework

This paper utilizes the MC-CGC event generator to demonstrate that LHCb data favors HERA-constrained initial conditions for the rcBK evolution equation and that the dense-dense $k_T$ factorization framework better describes mid-rapidity particle production than the dilute-dense DHJ approach, while also providing predictions for future ALICE FoCal measurements.

The Gist

Imagine two high-speed trains (protons) smashing into each other at the LHC, the world's most powerful particle accelerator. When they collide, they don't just bounce off; they explode into a shower of smaller particles. Physicists want to understand exactly how these particles are created, especially the ones flying off at extreme angles (forward) or right in the middle of the crash.

This paper is about building a digital simulation (a "movie maker" for particle collisions) to predict what happens in these crashes, based on a specific theory called the Color Glass Condensate (CGC).

Here is a breakdown of their work using simple analogies:

1. The Theory: The "Crowded Room" vs. The "Empty Hall"

In the world of particle physics, protons are made of tiny particles called gluons.

• The Problem: When protons move very fast, they get packed with so many gluons that they act like a super-dense, sticky fluid. This is the Color Glass Condensate.
• The Analogy: Imagine a crowded concert hall (the dense target) and a few people walking in from outside (the dilute projectile).
  • Old Theory (DHJ): This theory assumes the people walking in are sparse and just bump into the crowded room. It works well if the "projectile" is empty and the "target" is packed.
  • New Theory ($k_T$ Factorization): This theory assumes both sides are packed like two crowded concert halls crashing into each other.
• The Finding: The authors found that at the LHC's massive energies, both sides are actually packed. Therefore, the "two crowded halls" model ($k_T$) does a much better job of predicting what happens in the middle of the crash than the "sparse vs. crowded" model (DHJ).

2. The Tool: The "Event Generator" (MC-CGC)

The authors built a computer program called MC-CGC. Think of this as a video game engine for physics.

• Instead of just calculating one number, it simulates thousands of individual collisions, event by event.
• It takes the messy, complex math of the CGC theory and turns it into a step-by-step recipe:
  1. The Crash: Decide how many particles hit each other.
  2. The Spray: Simulate the initial explosion of particles.
  3. The Aftermath: Simulate how these particles radiate energy and stick together to form new, stable particles (hadrons) that detectors can actually see.
• This allows them to compare their "movie" directly with real footage from the LHCb experiment.

3. The "Starting Point" Mystery (Initial Conditions)

To run their simulation, they need to know what the proton looks like before the crash. They tested three different "blueprints" (initial conditions) for this starting state:

• Blueprint A (MV): The original, standard guess.
• Blueprint B & C (MV$\gamma$ and MV$_e$): Newer, more refined guesses based on data from the HERA electron accelerator.
• The Result: When they ran the simulation and compared it to real LHC data, Blueprints B and C were the winners. They matched the real data much better. Blueprint A (the original) predicted a "flatter" distribution of particles that didn't match reality, especially at higher speeds.

4. The "Forward" vs. "Middle" View

The paper looks at particles flying in two different directions:

• Forward (The "Edge" of the crash): Here, the "sparse vs. crowded" model (DHJ) works reasonably well. The data from the LHCb experiment favors the newer blueprints (MV$\gamma$ and MV$_e$).
• Middle (The "Center" of the crash): Here, the "two crowded halls" model ($k_T$) is clearly superior. The old model fails to describe the particle speeds correctly in the center.

5. Looking Ahead: The "FoCal" Crystal Ball

The authors didn't just look at past data; they used their winning simulation to predict what a future detector called FoCal (planned for the ALICE experiment) will see.

• They predicted how many neutral particles (like $\pi^0$, $\eta$, and $\omega$ mesons) and jets (tight bundles of particles) will be produced.
• They found that measuring these particles at very high speeds could help physicists fine-tune the "blueprints" of the proton even further, reducing the uncertainty in their theories.

Summary

In short, this paper says: "We built a better simulator for particle collisions. We tested three different starting assumptions, and the newer ones fit the real-world data best. We also proved that at the LHC's high energies, both colliding protons are so dense that we need a 'two-crowded-rooms' model to understand the middle of the crash, not the old 'sparse-visitor' model. Finally, we used our best model to predict what future experiments will see."

Technical Summary

Technical Summary: Forward Hadron Production in pp Collisions at LHC Energies from an Event Generator Based on the Color Glass Condensate Framework

Problem Statement
Understanding Quantum Chromodynamics (QCD) dynamics in high-energy collisions, particularly the emergence of parton saturation and the Color Glass Condensate (CGC) regime, remains a central challenge. While evidence for saturation exists from HERA, RHIC, and LHC data, alternative mechanisms can explain some observations, necessitating more robust theoretical tools. A specific gap exists in bridging the gap between complex analytical formalisms (such as NLO corrections and multi-point correlators) and practical phenomenological applications. Furthermore, there is a need to systematically assess the sensitivity of CGC predictions to initial conditions and to determine the appropriate factorization framework (dilute-dense vs. dense-dense) for different kinematic regions in proton-proton ($pp$) collisions at LHC energies.

Methodology
The authors utilize and update MC-CGC, a Monte Carlo event generator based on the leading-order (LO) CGC framework. The generator is implemented within the PYTHIA framework to handle non-perturbative processes, including beam remnants, initial/final state radiation (ISR/FSR), and Lund string fragmentation.

Key methodological components include:

1. Factorization Schemes: The model implements two distinct approaches for single-particle production:
   • Dilute-Dense (DHJ Factorization): Describes a dilute projectile parton scattering off a dense target. This uses collinear parton distribution functions (PDFs) for the projectile and dipole amplitudes for the target.
   • Dense-Dense ($k_T$-Factorization): Used when both projectile and target are in the dense regime (small-$x$). This employs unintegrated gluon distributions derived from dipole amplitudes for both sides.
2. Initial Conditions: The study systematically compares three parameterizations for the initial condition of the running-coupling Balitsky–Kovchegov (rcBK) evolution equation:
   • The McLerran–Venugopalan (MV) model.
   • The HERA DIS-constrained MV$\gamma$ model.
   • The HERA DIS-constrained MV$e$ model.
3. Event Generation: The number of hard scatterings per event is sampled from a negative binomial distribution (NBD) to reproduce charged-particle multiplicity fluctuations. The generator handles multiple parton interactions (MPI) by restricting events to at most one quark-initiated process (with others being gluon-initiated) to simplify color reconnection, though this is noted as an area for future improvement.
4. Calibration: Model parameters, including the normalization factor $K$, the Gaussian width of the thickness function, and the transverse momentum kick width ($\sigma_w$), are tuned to reproduce LHCb data at $\sqrt{s}=7$ TeV and RHIC data.

Key Contributions and Results

• Validation against Existing Data: The updated MC-CGC model successfully describes inclusive forward single-hadron production data from LHCb ($\sqrt{s}=7$ and $13$ TeV) and RHIC ($\sqrt{s}=200$ GeV). The model reproduces charged-particle multiplicity distributions ($P(N_{ch})$), pseudo-rapidity densities ($dN/d\eta$), and transverse momentum spectra ($dN/dp_T$) down to the ultra-soft region ($p_T \approx 0.5$ GeV).
• Sensitivity to Initial Conditions: The comparison of the three rcBK initial conditions reveals that current LHCb data favor the MV$\gamma$ and MV$e$ parameterizations. The original MV model produces a flatter $p_T$ spectrum that is disfavored by experimental data. The differences between the models become more pronounced at higher transverse momentum ($p_T$) and at mid-rapidity, where the small-$x$ evolution effects are less dominant.
• Factorization Framework Comparison: A critical finding is the performance difference between the DHJ and $k_T$-factorization frameworks at mid-rapidity ($|\eta| < 2.5$) at LHC energies. In this region, where both projectile and target partons reside in the small-$x$ regime (dense-dense scattering), the $k_T$-factorization framework provides a better description of the particle production spectra than the DHJ model. The DHJ model, relying on collinear PDFs for the projectile, shows clear discrepancies with mid-rapidity data unless the normalization is artificially adjusted, which then leads to overestimation at forward rapidities.
• Predictions for FoCal: The paper presents predictions for future measurements by the ALICE Forward Calorimeter (FoCal) at $\sqrt{s}=14$ TeV. These include:
  • Multiplicity distributions and mean transverse momentum $\langle p_T \rangle$ dependence for identified neutral mesons ($\pi^0, \eta, \omega$).
  • Jet observables (pseudo-rapidity distribution, $p_T$ spectra, and jet mass distributions) reconstructed from neutral pions.
  • The results indicate that jet observables, particularly at higher $p_T$, exhibit strong sensitivity to the choice of initial conditions, offering a potential pathway to constrain the shape of the dipole amplitude.

Significance and Claims
The paper claims that the MC-CGC event generator provides a consistent and flexible tool for bridging theoretical developments in saturation physics with phenomenological studies. By incorporating essential LO saturation features into a fully exclusive event generator, the authors demonstrate that the CGC framework can describe a wide range of observables across different collision energies and rapidity regions.

The study emphasizes that while the MV$\gamma$ and MV$e$ models are currently favored by data, the choice of initial condition significantly impacts predictions in the hard $p_T$ regime. Furthermore, the work highlights the necessity of using the dense-dense ($k_T$-factorization) framework for mid-rapidity $pp$ collisions at LHC energies, suggesting that the creation of a dense CGC medium from both colliding protons is a relevant physical picture in this kinematic domain. The authors conclude that their predictions for FoCal and jet observables offer concrete targets for future experimental constraints on the dipole amplitude and the validity limits of the DHJ formalism.

The authors remain modest regarding future developments, noting that the current model is limited to LO and requires extensions to include NLO corrections, detailed multiparton interactions, and generalization to $pA$ and $AA$ collisions to fully address the complexities of saturation dynamics.