⚛️ phenomenology

A basic model for high energy cosmic ray interactions

This paper presents a new, computationally efficient Monte Carlo generator for high-energy cosmic ray interactions based on Reggeon Field Theory, designed to offer a transparent and tunable framework for studying extensive air showers while maintaining consistency with fundamental physical conservation laws and accelerator data.

Original authors: Sergey Ostapchenko, Tanguy Pierog, Günter Sigl

Published 2026-03-16

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Sergey Ostapchenko, Tanguy Pierog, Günter Sigl

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

The Big Picture: The "Cosmic Ray Detective" Problem

Imagine the Earth is a giant, invisible shield. Every day, invisible bullets (called Cosmic Rays) from deep space smash into our atmosphere. These bullets are so energetic that when they hit an air molecule, they don't just bounce off; they explode, creating a massive shower of smaller particles (like a cosmic firework) that rains down on the ground. This is called an Extensive Air Shower (EAS).

Scientists on the ground can't catch the original bullet because it's gone. They can only see the "debris" (the shower) hitting the ground. To figure out what the original bullet was made of and how fast it was going, they have to work backward, like a detective trying to guess what a car looked like by studying the skid marks and broken glass left after a crash.

To do this, they need a simulation (a computer program) that predicts exactly how that crash happens. For decades, these simulations have been like "Black Boxes." You put data in, and you get an answer out, but no one really knows how the computer decided on that answer. If the simulation doesn't match real-world observations, scientists are stuck because they can't easily tweak the "Black Box" to fix it without breaking the laws of physics.

The Solution: The "Transparent Blueprint" (QGSb)

The authors of this paper built a new simulation tool called QGSb. Think of it not as a Black Box, but as a transparent, Lego-like blueprint.

The Goal: They wanted a model that is simple enough to understand, flexible enough to be tweaked by anyone, but still follows the strict rules of physics (like conservation of energy and momentum).
The Analogy: Imagine trying to predict how a house of cards will collapse.
- Old Models: You have a super-complex robot that simulates every single card's friction and air resistance. It's accurate, but if the prediction is wrong, you can't easily tell which card caused the problem.
- The QGSb Model: You use a simplified set of rules: "If the wind blows this hard, the top card falls." It's less detailed, but you can see exactly why the house fell, and you can easily change the wind speed to see what happens.

How It Works: The "Traffic Jam" and the "String"

The paper uses a theoretical framework called Reggeon Field Theory. Let's break that down into a traffic analogy.

1. The Pomeron (The Invisible Traffic Jam)
When two particles smash into each other at near-light speed, they don't just bounce. They create a chaotic mess of energy. In this model, the authors imagine this mess as a "traffic jam" of invisible strings called Pomerons.

Soft Pomerons: These are like slow-moving, heavy trucks. They dominate at lower energies and cause a gentle, slow increase in chaos.
Semihard Pomerons: These are like fast sports cars. At very high energies, these take over, causing a rapid explosion of new particles.
The Innovation: The QGSb model uses both types of traffic jams. It assumes that at low speeds, the trucks rule, but as you speed up, the sports cars take over. This allows the model to handle everything from slow collisions to the most violent cosmic crashes.

2. The String Fragmentation (The Spaghetti Break)
Once the "traffic jam" (Pomeron) happens, the energy has to turn into actual particles (like pions, protons, and neutrons).

The Analogy: Imagine stretching a piece of sticky taffy or a rubber band between two hands. As you pull the hands apart, the string gets thinner and thinner until it snaps.
The Physics: When the string snaps, it doesn't just disappear. It creates new pairs of particles (like a new piece of taffy forming in the middle). The model calculates exactly how many pieces of "taffy" (particles) are created and how fast they fly away.

3. The "Ghost" Pions (The Invisible Messenger)
One of the paper's specific improvements is handling a tricky process where a particle emits a "ghost" pion (a temporary, invisible particle) that bounces off the target before disappearing.

The Analogy: Imagine throwing a ball at a wall, but before it hits, it bounces off a friend standing in front of the wall. The old models often ignored this "friend." The new model explicitly calculates this bounce, which turns out to be crucial for predicting how many muons (a type of heavy electron) reach the ground.

Why Does This Matter?

1. Solving the "Muon Puzzle"
For years, scientists have been confused. Real cosmic ray showers hitting the ground have more muons than the best computer models predict. It's like the "Black Box" simulations are underestimating the debris.

The QGSb Result: Because this new model is transparent, the authors could tweak the "taffy snapping" rules to match real accelerator data. They found that their model predicts more muons than previous models, bringing the simulation closer to reality.

2. The "Tuning Knob" Advantage
The biggest strength of QGSb is that it's user-friendly.

Old Way: If a scientist wanted to change a model to fit new data, they had to hire a super-expert to rewrite the code, risking a crash.
New Way: With QGSb, a scientist can turn a "knob" (a parameter) to adjust how the strings break or how the traffic jams behave. They can test: "What if the strings break differently? Does that match the new data?" This helps them figure out if the universe is behaving differently than we thought, or if our old models were just wrong.

The Bottom Line

This paper introduces a new, simpler, and more transparent tool for simulating cosmic ray crashes. It's like giving astronomers a see-through engine instead of a sealed black box.

It explains how particles smash together using the concept of "traffic jams" (Pomerons) and "stretching strings."
It successfully predicts how many particles rain down on Earth, solving some long-standing mysteries about muon counts.
Most importantly, it allows scientists to tune the model themselves to see if the universe is behaving in ways we haven't expected yet.

It's a step toward understanding the most violent events in the universe by building a clearer, more honest map of how they work.

1. Problem Statement

The study of Ultra-High Energy Cosmic Rays (UHECRs) relies on indirect detection via Extensive Air Showers (EAS). Interpreting EAS data requires accurate Monte Carlo (MC) generators for hadronic interactions. However, existing sophisticated models (e.g., EPOS, QGSJET) are often viewed as "black boxes" by experimentalists due to their complexity and lack of transparency.

The Issue: Recent UHECR observations cannot be consistently described by current models. Experimental groups often resort to ad-hoc rescaling of EAS observables to fit data, which neglects correlations between different observables and risks violating fundamental conservation laws or conflicting with accelerator data.
The Need: A new interaction model is required that offers transparent physics, sufficient parameter freedom for retuning, and high computational efficiency, while maintaining consistency with unitarity, conservation laws, and existing accelerator data.

2. Methodology: The QGSb Model

The authors developed QGSb (Quark-Gluon Strings basic), a new MC generator based on Reggeon Field Theory (RFT). The model aims to balance simplicity with physical rigor by returning to the concepts of the original Quark-Gluon String and Dual Parton models.

Core Physics Framework

Pomeron Exchange: Hadronic collisions are treated as mediated by Pomeron exchanges. The model employs a two-component Pomeron amplitude:
1. Soft Pomeron: Corresponds to low-virtuality parton cascades with a slow energy rise.
2. Semihard Pomeron: Corresponds to high-virtuality parton cascades with a rapid energy rise.
- Note: Explicit perturbative parton cascades and nonlinear corrections are not calculated microscopically; their effects are effectively "renormalized" into the two-component amplitude.
Good-Walker (GW) Formalism: Hadrons are treated as superpositions of Fock states with different transverse sizes and couplings. This accounts for inelastic diffraction and color fluctuation effects (screening), leading to a milder energy rise in cross-sections and larger fluctuations in particle production.
Reggeon Exchange: To extend applicability to low energies ( $<100$ GeV), the model includes secondary Reggeon exchanges (using an "effective Reggeon" approach) alongside Pomeron exchanges.
Hadronization: Secondary hadrons are produced via the fragmentation of color strings stretched between constituent partons. The model distinguishes fragmentation parameters between soft and semihard Pomerons to account for different parton densities.

Key Technical Features

Cross Sections: Calculated using eikonal formalism in impact parameter space, incorporating GW decomposition for total, elastic, and inelastic cross-sections.
Diffraction:
- Low Mass: Handled via GW approach.
- High Mass: Modeled via a diffractive cut of the triple-Pomeron graph.
Pion Exchange (RRP): A specific improvement over previous models (like QGSJET-III) involves explicitly treating $t$ -channel pion exchange in Reggeon-Reggeon-Pomeron configurations. This includes rescattering of virtual pions and impact-parameter-dependent absorptive corrections, crucial for predicting muon content in EAS.
Low Energy Extension: Below 100 GeV, Reggeon-dominated processes (planar diagrams or "undeveloped cylinder" contributions) are used for particle production.

3. Key Contributions

Transparent Formalism: The model avoids complex microscopic treatments of parton cascades in favor of a phenomenological but theoretically grounded RFT approach, making it easier for external users to understand and modify.
Parameter Flexibility: The model features a set of adjustable parameters (e.g., Pomeron intercepts, coupling strengths, fragmentation parameters) that can be retuned by users to explore uncertainties or fit specific datasets without breaking fundamental constraints.
Unified Energy Range: It covers a vast energy spectrum from $\sim 10$ GeV to $10^{11}$ GeV (lab frame), seamlessly transitioning between Reggeon-dominated low-energy physics and Pomeron-dominated high-energy physics.
Improved Pion Exchange: The explicit treatment of virtual pion rescattering provides a more robust prediction for forward particle production, directly impacting EAS muon predictions.

4. Results and Validation

The model parameters were tuned using a wide variety of accelerator data (NA49, NA61/SHINE, TOTEM, LHCf, ATLAS, HARP).

Cross Sections: QGSb successfully reproduces total, elastic, and diffractive cross-sections for $pp$, $\bar{p}p$ , $\pi p$ , and $Kp$ interactions across the energy range. It also accurately predicts hadron-nucleus cross-sections using the Glauber-Gribov approach.
Particle Production:
- Forward Production: Good agreement with LHCf data for neutral pions and neutrons.
- Diffraction: Accurately describes single diffractive-like events at $\sqrt{s}=7$ TeV (TOTEM data).
- Low Energy: Reasonable agreement with NA61/SHINE and HARP data for rapidity and momentum distributions, though some excesses in forward charged pion production at small angles were noted.
Extensive Air Shower (EAS) Predictions:
- Muon Content ( $N_\mu$ ): QGSb predicts a higher muon number at sea level compared to QGSJET-II/III, approaching EPOS-LHC-R values. This is attributed to enhanced (anti)nucleon production in the forward region.
- Shower Depth ( $X_{max}$ ): The model predicts a smaller elongation rate ( $d \ln X_{max} / dE_0$ ) at the highest energies compared to previous models. This is linked to the energy-dependent suppression of elastic virtual pion scattering and the specific treatment of Reggeon contributions.
- Fluctuations: The fluctuations in $X_{max}$ ( $\sigma(X_{max})$ ) agree well with other models, constrained by precise LHC cross-section measurements.

5. Significance

Bridging the Gap: QGSb serves as a vital tool for experimentalists to test "what-if" scenarios. It allows them to retune interaction parameters to fit specific EAS data and immediately check the consistency of these modifications against accelerator constraints.
Addressing the Muon Puzzle: By predicting higher muon numbers, QGSb offers a potential explanation for the "muon excess" observed in UHECR data by the Pierre Auger Observatory, suggesting that standard models may underestimate forward baryon production.
Computational Efficiency: The simplified formalism allows for faster simulations, enabling automated tuning procedures using combined accelerator and cosmic ray data.
Astrophysical Applications: Beyond EAS, the generator is suitable for calculating secondary particle spectra (electrons, neutrinos, photons) for cosmic ray propagation studies in the interstellar medium.

In conclusion, the QGSb model provides a robust, transparent, and flexible framework for high-energy hadronic interactions, specifically designed to facilitate the interpretation of cosmic ray data while maintaining strict adherence to fundamental physical principles and accelerator constraints.