Mass Composition Studies of the Highest Energy Cosmic Rays

The determination of the mass composition of the highest energy cosmic rays is one of the greatest challenges in cosmic ray experiments. The highest energy cosmic rays are only detected indirectly because of their very low flux. Using the atmosphere as a large target, Air Fluorescence Detectors are capable of tracing the evolution of the size of the Extensive Air Shower through the atmosphere (the shower longitudinal profile). The analysis of the characteristics of the detected longitudinal profiles is currently the most reliable way for extracting some information about the primary cosmic ray mass composition. In this proceeding, I will describe in some detail the Pierre Auger elongation rate studies, and I will show the potential for mass composition studies using the surface and the fluorescence detectors information as part of a single analysis. The interpretation of the current data with regard to mass composition, relies heavily on high energy hadron interaction models. Using standard hadron interaction models, the data suggest that the composition becomes lighter up to about 2 $\times$ 10$^{18}$ $eV$ and above that it becomes heavier again. This apparent change in the mass composition at 2 $\times$ 10$^{18}$ $eV$ seems to be correlated with a spectrum index change in the observed energy spectrum.

💡 Research Summary

The paper presents a comprehensive study of the mass composition of the highest‑energy cosmic rays using data from the Pierre Auger Observatory (PAO). Because the flux of particles above 10¹⁸ eV is extremely low, direct detection is impossible; instead, the atmosphere serves as a giant calorimeter. When a primary cosmic‑ray particle strikes the atmosphere it initiates an extensive air shower (EAS) that develops longitudinally as a cascade of hadrons, muons and an electromagnetic component. The depth in the atmosphere at which the number of charged particles (or the deposited energy) reaches its maximum, denoted Xmax, is the most sensitive observable for inferring the primary mass. Light nuclei (protons) have a larger Xmax (deeper in the atmosphere) because each nucleon carries the full primary energy, while heavy nuclei (e.g., iron) fragment into many lower‑energy sub‑showers that reach their maximum earlier (smaller Xmax).

PAO combines a surface detector (SD) array of 1600 water‑Cherenkov stations with four fluorescence detector (FD) sites, each comprising six telescopes equipped with 440 PMTs. Hybrid events—those simultaneously recorded by the FD and at least one SD station—allow a precise reconstruction of the shower geometry (angular resolution ≈ 0.6°) and consequently a high‑precision measurement of Xmax (average resolution ≈ 20 g cm⁻² after quality cuts). The longitudinal profile is obtained as dE/dX versus slant depth and fitted with the four‑parameter Gaisser‑Hillas function, yielding Xmax, the maximum deposited energy, and shape parameters λ and X₀.

A major systematic effect stems from the limited elevation field of view of the FD (≈ 2°–30°). Showers whose Xmax lies outside this window are either missed or produce a biased Xmax distribution. To mitigate this, the authors introduce two geometry‑dependent depth limits, X_low and X_up, which define the portion of the shower axis that is observable. By studying the dependence of the mean Xmax on these limits in real data, they identify flat regions where the measurement is unbiased and set energy‑dependent cuts on X_low and X_up. Only events with Xmax “bracketed” (i.e., lying within the observed profile) and with an observed grammage length > 320 g cm⁻² are retained, ensuring an Xmax uncertainty below 40 g cm⁻².

Monte‑Carlo simulations using several high‑energy hadronic interaction models (QGSJET01, QGSJETII‑03, SIBYLL 2.1, EPOS 1.6) are processed through the full detector simulation and reconstruction chain. The reconstructed mean Xmax values agree with the true simulated values for pure proton, pure iron, and mixed (50 %/50 %) compositions, confirming that the anti‑bias cuts effectively remove geometric biases. The minimal energy at which an unbiased ⟨Xmax⟩ can be measured is just above 10¹⁷ eV.

The core result is the energy dependence of the measured mean Xmax. A simple linear fit ⟨Xmax⟩ = A + D₁₀·log₁₀E yields an elongation rate D₁₀ = 54 ± 2 g cm⁻² per decade, but the χ²/ndf = 24/13 indicates that a single straight line does not adequately describe the full energy range (probability < 3 %). A broken‑line fit provides a much better description (probability ≈ 63 %). The break occurs at ≈ 2 × 10¹⁸ eV. Below this energy the elongation rate is 71 ± 5 g cm⁻² decade⁻¹, while above it drops to 40 ± 4 g cm⁻² decade⁻¹.

Interpreting these rates with the predictions of the hadronic models shows that, up to ~2 × 10¹⁸ eV, the data are compatible with a light composition (proton‑dominated). Above this energy the trend moves toward heavier nuclei, consistent with an increasing iron fraction. This transition coincides with the “ankle” feature observed in the cosmic‑ray energy spectrum, suggesting a possible link between the spectral change and a shift in source populations or acceleration mechanisms.

Systematic uncertainties—stemming from atmospheric transmission, fluorescence yield, energy scale, and the model‑dependent “missing energy” correction—are estimated to be ≤ 15 g cm⁻² at low energies and ≤ 12 g cm⁻² at high energies. The interpretation remains model‑dependent; different interaction models predict different absolute ⟨Xmax⟩ values, so the exact composition fractions cannot be uniquely determined without further constraints.

In summary, the Auger hybrid analysis, with rigorous anti‑bias selections and detailed MC validation, provides a robust measurement of the elongation rate and reveals a composition transition near 2 × 10¹⁸ eV. The findings underscore the importance of improving hadronic interaction models and increasing statistics to resolve the remaining ambiguities in the mass composition of ultra‑high‑energy cosmic rays.