Refining the Greisen Profile for Low-Energy Cosmic Gamma-Rays: Quantifying Deviations Across Altitudes and Zenith Angles

This study refines the classical Greisen profile by introducing a zenith-angle-dependent empirical correction to the shower age parameter, demonstrating that the modified model reduces deviations in particle number predictions for low-energy (20–800 GeV) cosmic gamma-ray showers to below 4.7% compared to CORSIKA simulations, thereby offering a more accurate and computationally efficient tool for high-altitude observatories like HAWC and CONDOR.

The Gist

Imagine the Earth's atmosphere as a giant, invisible forest. When a high-energy cosmic gamma-ray (a tiny, super-fast particle from deep space) hits the top of this forest, it doesn't just stop. Instead, it crashes into an air molecule and triggers a massive chain reaction, creating a "shower" of billions of new particles that rain down toward the ground. Scientists call this an Extensive Air Shower.

To understand these showers, scientists need a map. For decades, they've used a famous map called the Greisen Profile. Think of this map like an old, classic recipe for baking a cake. It works perfectly for huge, high-energy cakes (ultra-high energy particles), but when you try to use it to bake smaller, more delicate cakes (low-energy particles between 20 and 800 GeV), the recipe starts to get things wrong. It predicts the cake will be bigger or smaller than it actually is.

The Problem: The "Old Recipe" vs. Reality

The authors of this paper noticed that the old Greisen recipe was struggling in two specific situations:

1. High Altitudes: Observatories like HAWC and the proposed CONDOR array sit on top of very high mountains (5,000 to 5,900 meters up). The air there is thin, like being on a high peak of a mountain. In thin air, the "cake" (the particle shower) develops differently than it does at sea level.
2. Angles: Sometimes the cosmic rays don't fall straight down (like rain in a storm); they come in at a slant (like rain driven by wind). The old recipe didn't account well for this slant, especially when the angle was steep.

Because of these factors, the old recipe was off by as much as 12.5%. In science, that's a big error. It's like a GPS telling you to turn left when you actually need to turn right, or a weather app saying it will be sunny when it's actually pouring rain.

The Solution: A "Modified Recipe"

The team created a Modified Greisen Profile. Think of this as taking the old recipe and adding a few "tweaks" or "adjustments" to the ingredients. Specifically, they adjusted a variable called the "shower age" (which tracks how old or mature the particle shower is) to account for:

• The fact that thin air at high altitudes changes how fast the shower grows.
• The extra distance particles have to travel when coming in at a slant.
• The energy lost by particles as they bump into air molecules (ionization losses), which the old recipe ignored for low-energy particles.

The Test: The "Gold Standard" Simulation

To see if their new recipe worked, they didn't just guess. They compared it against CORSIKA, which is like a super-powerful, high-definition video game simulator. CORSIKA simulates every single particle interaction in the atmosphere with extreme detail. It is the "gold standard" or the "truth" that scientists trust.

They ran over one million simulations (1,008,000 to be exact) covering different energies, different mountain heights, and different angles. Then, they compared the predictions of the Old Recipe and the New Recipe against the simulator's results.

The Results: A Much Better Map

The results were clear:

• The Old Recipe: Still made mistakes, with errors ranging up to 12.5%. It was getting the size of the particle shower wrong, especially at steep angles.
• The New Recipe: The modified profile was much more accurate. It kept the errors down to less than 4.7%.

In everyday terms, if the old map said a mountain was 1,000 meters high but it was actually 1,125 meters, the new map says it's 1,047 meters. That's a huge improvement in precision.

Why This Matters

The paper concludes that this new, modified formula is a better tool for scientists working at high-altitude observatories.

• It's Fast: Unlike the super-powerful simulator (CORSIKA), which takes a lot of computer power and time to run, this new formula is a simple math equation. It's like using a quick mental math trick instead of a supercomputer.
• It's Accurate: Because it accounts for the thin air and the angles, it gives a much truer picture of what is happening when cosmic rays hit the atmosphere.

This allows astronomers to better understand the energy and origin of cosmic gamma rays, helping them study the universe's most energetic events with a sharper, more reliable lens. The paper does not claim this changes medical treatments or future technologies; it strictly improves the mathematical tools used to study cosmic rays right now.

Technical Summary

Technical Summary: Refining the Greisen Profile for Low-Energy Cosmic Gamma-Rays

Problem Statement
Extensive air showers (EAS) initiated by cosmic gamma rays are fundamental to astroparticle physics, serving as probes for primary energy and interaction mechanisms. While the classical Greisen formalism provides a compact analytical description of electromagnetic cascade development, its accuracy is limited in the low-energy regime (20–800 GeV). At these energies, ionization losses become significant relative to radiative processes, and particle multiplicity is reduced, leading to discrepancies between the analytical model and physical reality. Furthermore, existing models often lack systematic validation across the specific conditions relevant to high-altitude observatories (5000–5900 m), such as the HAWC array and the proposed CONDOR array. In these environments, reduced atmospheric density and varying zenith angles (0°–40°) alter cascade dynamics, shifting shower maxima and affecting particle numbers in ways the classical profile does not fully capture.

Methodology
The study employs a comparative approach, validating both the classical Greisen profile and a proposed modified version against high-fidelity Monte Carlo simulations.

• Simulation Framework: The authors utilized CORSIKA (COsmic Ray SImulations for KAscade) to simulate gamma-ray induced air showers. The simulations covered primary energies from 20 to 800 GeV, four specific altitudes (5000 m, 5300 m, 5600 m, and 5900 m), and zenith angles ranging from 0° to 40° in 2° increments. A total of 1,008,000 simulation runs were performed to ensure statistical robustness, utilizing the EGS4 electromagnetic interaction module and the US Standard Atmosphere adjusted for high-altitude density gradients.
• Analytical Models:
  • Classical Profile: Based on Approximation A (radiative processes) and B (ionization losses), the classical profile calculates the number of particles $N(t)$ as a function of atmospheric depth $t$ and shower age $s$.
  • Modified Profile: The authors introduce an empirical correction to the shower age parameter $s$ within the Greisen formalism. This modification incorporates a constant $c'$ (0.99745) into the logarithmic term of the shower age equation to better account for low-energy ionization losses and atmospheric density variations.
• Quantification: Both profiles were fitted to the CORSIKA data using a scaling parameter $\alpha$. The deviation from the simulation was quantified as $|\alpha - 1| \times 100$, and the goodness of fit was assessed using the reduced chi-squared statistic ($\chi^2/\text{n.d.f.}$).

Key Contributions

1. Systematic Validation: The paper provides a comprehensive validation of the Greisen formalism specifically for the 20–800 GeV energy range, a regime where fewer particles reach the ground due to reduced multiplicity and increased ionization losses.
2. Empirical Modification: The study proposes a modified Greisen profile that adjusts the shower age term to better reflect the accelerated particle absorption observed in the post-maximum phase of low-energy showers.
3. Altitude and Zenith Dependence: The work explicitly quantifies deviations across a range of high altitudes (5000–5900 m) and zenith angles (0°–40°), addressing the specific geometric and atmospheric challenges faced by high-altitude observatories.

Results
The comparison between the analytical profiles and CORSIKA simulations yields the following quantitative findings:

• Deviation Reduction: The modified Greisen profile consistently achieves deviations in particle numbers below 4.7% across all tested altitudes and zenith angles. In contrast, the classical profile exhibits deviations of up to 12.5%, particularly at larger zenith angles (e.g., 40°) where path-length effects are amplified.
• Fit Quality: The modified profile demonstrates lower $\chi^2/\text{n.d.f.}$ values (ranging from 0.02 to 0.16) compared to the classical profile, indicating a statistically superior fit to the simulation data.
• Specific Performance: At 5300 m and a 40° zenith angle, the classical profile shows an $\alpha$ value of 1.1258 (a 12.58% overestimation), whereas the modified profile yields an $\alpha$ of 1.0463 (a 4.63% deviation).
• Shower Maximum: The modified profile more accurately captures the shift in the shower maximum toward lower atmospheric depths caused by reduced air density at high altitudes.

Significance
The paper claims that the modified Greisen profile offers a computationally efficient alternative to resource-intensive Monte Carlo simulations, making it a practical tool for high-altitude observatories like HAWC and the CONDOR array. By reducing deviations to below 4.7%, the model enhances the accuracy of particle number predictions, which is critical for:

• Reconstructing primary gamma-ray energies.
• Improving detection efficiency.
• Discriminating gamma-ray signals from hadronic backgrounds.

The authors conclude that this refined analytical framework provides a physically consistent description of shower development under high-altitude conditions, supporting studies of cosmic accelerators and dark matter searches. They note, however, that the model relies on a simplified atmospheric description and does not currently account for seasonal variations in density or temperature. Future work is suggested to extend the model to higher energies and validate it against observational datasets.