Superposition model for energy reconstruction and mass identification in cosmic ray spectra

arXiv:2602.23637 · astro-ph.HE, astro-ph.IM · Submitted 2026-02-27 · Read on arXiv

Listen

Transcript

Introduction to the show: ident: Astrophysics Radio. Generated commentary on the latest astrophysics papers.

Vera: Today's paper: "Superposition model for energy reconstruction and mass identification in cosmic ray spectra".

Jocelyn: The gist The method introduces a novel method for reconstructing energy and logarithm mass (lnA) based on a superposition model for cosmic ray spectra How it works The reconstruction of…

Vera: First, who's behind it and why it matters.

Paper summary: Vera: So to recap this paper on "Superposition model for energy reconstruction and mass identification in cosmic ray spectra", they introduce a novel technique based on a superposition model to reconstruct both the energy and logarithm mass of cosmic rays.

Jocelyn: The core thesis is that instead of using integrated particle counts, they use particle densities derived from electromagnetic and muon detectors to get better resolution for these properties. They claim this density-based approach improves reconstruction, particularly for heavier nuclei.

Subrahmanyan: The paper is organized around describing the KM2A detector setup and then presenting the two formulas used to fit the lateral distributions of these detectors, which are then used to derive energy and log mass variables.

Vera: They claim that if the superposition model holds, for all primary compositions, the relationship between (N A E e/mu /A) versus (E/A) simplifies down to a single universal curve—the one for a proton (<ref:2602.23637#pg5>).

Jocelyn: That universality is key because it means the reconstruction formulas, specifically equations five and six become very simple when you relate (rho ED) and (rho MD) to the energy variables (<ref:2602.23637#pg5>).

Subrahmanyan: Essentially, they are showing that the physics of how these nuclei interact with the shower can be described by a single underlying relationship, which makes the reconstruction cleaner than if we had to treat every nucleus type separately.

Vera: The paper sets up this framework using rho ED and rho MD, where these densities are derived from fitting the lateral distributions of the ED and MD hits, respectively (<ref:2602.23637#pg3>).

Jocelyn: They define resolution for density as the sigma of a Gaussian function fitted to the distribution of the logarithm of that density, and they chose specific distance parameters—one hundred m for rho ED and one hundred fifty m for rho MD —to simplify things across all particle types (<ref:2602.23637#pg4>).

Subrahmanyan: So, the main claim is that by using these two density variables, they can reconstruct energy and log mass simultaneously through fitting equations one and two (<ref:2602.23637#pg3>).

Vera: And the performance evaluation shows that the bias for energy stays within ±five percent up to one hundred PeV, with better results at lower zenith angles, which is a very concrete statement about the accuracy of this new method (<ref:2602.23637#pg1>).

Jocelyn: It matters because it gives us a robust way to look at the data from experiments like LHAASO and provides reliable values for energy spectra and composition identification in that high-energy regime.

Subrahmanyan: It’s about providing a tool where the hadronic model dependencies scale with (E/A) and are nearly independent of primary composition, which is what makes this framework theoretically sound (<ref:2602.23637#pg1>).

Conclusion: Vera: So, wrapping up this discussion on "Superposition model for energy reconstruction and mass identification in cosmic ray spectra", the main thing is that this method offers a consistent way to measure both the energy and the log mass of cosmic rays using density measurements from particle detectors.

Jocelyn: The authors are showing that by treating these densities through a superposition model, they can connect them back to a single universal curve based on protons, which simplifies how we think about the composition dependence (<ref:2602.23637#pg5>).

Subrahmanyan: It’s really about giving us a reliable way to reconstruct A from the distance between points A and B in a specific coordinate system, which allows for that combined measurement of energy and log mass (<ref:2602.23637#pg1>).

Vera: The implication is that this method gives us a universal energy contour line, one that doesn't depend on the primary energy or composition, which helps us pin down those systematic uncertainties in our measurements (<ref:2602.23637#pg1>).

Jocelyn: So, when you think about it simply, this paper is giving us a better set of tools to analyze cosmic ray data from observatories like LHAASO, helping us get a clearer picture of the energy and mass distribution at these very high energies.

Subrahmanyan: It provides important input for measuring individual energy spectra because the reconstruction parameters are not overly dependent on what nucleus we assume it is when we calculate the results (<ref:2602.23637#pg1>).

Vera: The authors acknowledge that the two universal calibration lines are actually the main source of systematic uncertainty, which is a necessary piece of information for anyone using this method to understand its limits (<ref:2602.23637#pg1>).

Jocelyn: So, ultimately, this work gives us a strong framework for combining electromagnetic and muon detector data to reconstruct the fundamental properties of these high-energy cosmic rays.

Subrahmanyan: That framework is important because it helps us understand the origin of phenomena like the cosmic ray knee by giving us better inputs on energy and composition (<ref:2602.23637#pg1>).

School of Physical Science and Technology, Southwest Jiaotong University · Key Laboratory of Particle Astrophysics & Experimental Physics Division & Computing Center, Institute of High Energy Physics, Chinese Academy of Sciences · TIANFU Cosmic Ray Research Center

astro-ph.HE, astro-ph.IM

Submitted: 2026-02-27

Updated: 2026-10-08

Comments: 16 pages, 16 figures, Accepted for publication in Physical Review D

License: http://arxiv.org/licenses/nonexclusive-distrib/1.0/

Importance score: 77/100

The gist: The gist The method introduces a novel method for reconstructing energy and logarithm mass (lnA) based on a superposition model for cosmic ray spectra How it works The reconstruction of energy and

Key concepts

Superposition Model: This model approximates a composite nucleus as many independent nucleons. It suggests that the relationship between the logarithm of the nuclear mass number and energy scales follows a single, universal curve, similar to that of protons, regardless of the primary nucleus's composition.
Particle Densities ($ ho_{ED}$ and $ ho_{MD}$): Instead of using integrated particle counts in annular bands, this method uses particle densities measured by the electromagnetic (ED) and muon (MD) detectors at a fixed distance from the shower axis. This density-based approach improves resolution for energy and lnA reconstruction.
Calibration Lines: The reconstruction relies on two universal calibration lines derived from particle densities. These lines are used to simplify the relationship between measured densities ($ ho_{ED}$ and $ ho_{MD}$) and the physical parameters of energy ($E$) and logarithm mass ($ ext{ln}A$).
Lateral Distribution Fitting: The reconstruction involves fitting the lateral distribution of ED and MD hits using specific equations. Variables like $ ho_{ED}$ (for ED) and $ ho_{MD}$ (for MD) are derived from these fits, which are then used to determine the energy and lnA values.

Terminology

Core — 16× in this paper · explained in 2 episodes
Calibration — 7× in this paper · explained in 7 episodes
Monte Carlo — 4× in this paper
Large High Altitude Air Shower Observatory — 3× in this paper
Accuracy — 1× in this paper · explained in 1 episode
alpha — 1× in this paper · explained in 1 episode

Summary

The gist The method introduces a novel method for reconstructing energy and logarithm mass (lnA) based on a superposition model for cosmic ray spectra

How it works

The reconstruction of energy and lnA is based on using two universal, composition- and energy-independent calibration lines The method utilizes particle densities—measured by LHAASO’s electromagnetic and muon detectors at a fixed distance from the shower axis—rather than integrated particle counts in annular bands This density-based approach improves resolution for both energy and lnA, especially for heavy nuclei The resulting energy resolution ranges from below 5% to ∼ 15% above 1 PeV, with the best mass resolution for iron achieved being below 25% above 10 PeV

Detector Description and Simulation

The KM2A array consists of an electromagnetic particle detector (ED) array and a muon detector (MD) array The ED unit consists of four plastic scintillation tiles covered by a 5-mm-thick lead plate to absorb low-energy charged particles and convert γ-rays into electron-positron pairs, improving angular and core position resolution The MD array is composed of 1188 water Cherenkov tanks deployed on a grid with a spacing of 30 m Simulations were performed using the CORSIKA software (v77410) [28] with hadronic interaction models QGSJET-II-04, EPOS-LHC, and SIBYLL 2.3d for hadronic interactions

Reconstruction Methodology

The reconstruction procedure involves fitting the lateral distribution of ED and MD hits to derive variables used for energy and lnA reconstruction For the ED reconstruction, equation 1 is used to fit the lateral distribution, where ρem is the density of particles measured by the ED, rm is fixed at 130 m, and r is the perpendicular distance to the shower axis The MD lateral distribution is fitted with equation 2 using a maximum likelihood method similar to that employed for ED reconstruction The variables used for energy and lnA reconstruction are ρED and ρMD, which are derived from these fits

Energy and lnA Reconstruction

The two variables of ρED and ρMD, derived from the fitted lateral distributions, will be used for energy and lnA reconstruction The resolution of density is defined as the sigma of a Gaussian function fitted to the distribution of the logarithm of the density To simplify reconstruction, 100 m for density measured by ED (denoted as ρED below) and 150 m for density measured by MD (denoted as ρMD below) were chosen for all types of primary particles in all energy range and zenith angle range in this work

Superposition Model Application

The superposition model approximates a composite nucleus as a collection of independent nucleons, each carrying a fraction of the total energy The relationship between lg(N A,E e/µ /A) and lg(E/A) for a primary nucleus of mass number A reduces to lg(N P,E/A e/µ) versus lg(E/A) If the superposition model holds, for all compositions, lg(N A,E e/µ /A) versus lg(E/A) follows a single, universal curve: that of a proton This relationship is expressed as lg(ρED) − αelgA ≡ fe(E/A) (5) and lg(ρMD) − αµlgA ≡ fµ(E/A) (6)

Performance Evaluation

The performance evaluation shows that the bias of energy is within ±5% for all masses and zenith angles from 300 TeV to 100 PeV, with a minimum bias of ±2% occurring at zenith angles of 35o−40o For lnA reconstruction, the bias is within 0.3 for all masses and zenith angles from 300 TeV to 100 PeV The resolution of lnA for iron is less than 0.4 above 1 PeV at all zenith angles, with the best resolution less than 0.25 (corresponding to a 25% mass resolution) above ∼ 10 PeV The hadronic model dependencies of energy and lnA scale with lg(E/A) and are nearly independent of primary composition, as predicted by the superposition model

Conclusion

The method’s good energy and lnA performance provides important input for measuring individual energy spectra The two universal calibration lines are the main source of systematic uncertainty in the reconstruction The reconstructed parameters are used to reconstruct lnA from the distance between points A and B in a specific coordinate system This allows for the reconstruction of both energy and lnA The resulting energy contour line is universal and independent of primary energy and composition, defining the main sources of systematic uncertainty in energy reconstruction The method's good performance provides important input for measuring individual energy spectra

--- Page 1 ---

The gist The method introduces a novel method for reconstructing energy and logarithm mass (lnA) based on a superposition model for cosmic ray spectra

Detector Description and Simulation

Improvements for AI systems

Bold header: Superposition Model Integration for Energy Reconstruction

The improved AI system can reconstruct energy and logarithm mass (lnA) using two universal, composition- and energy-independent calibration lines, which are derived from fitting lateral distributions of electromagnetic particles (ED) at 100m and muons (MD) at 150m.

Bold header: Density-Based Measurement Superiority

The system can prioritize particle densities—measured by LHAASO’s electromagnetic and muon detectors at a fixed distance from the shower axis—rather than integrated particle counts in annular bands, as this approach improves resolution for both energy and lnA, especially for heavy nuclei.

Bold header: Universal Curve Application

The system can utilize the derived relationships, such as lg(ρED) − αelgA ≡ fe(E/A) (5) and lg(ρMD) − αµlgA ≡ fµ(E/A) (6), which demonstrate that for all compositions, the relationship between log density and log energy is nearly independent of primary composition, as predicted by the superposition model.

Bold header: Composition-Independent lnA Reconstruction

The system can reconstruct lnA using the relationship derived from combining universal curves, as shown in lg(ρMD/Aαµ) versus lg(ρED/Aαe) also follows a universal curve for all compositions, with a residual within about 2%.

Bold header: Hadronic Model Uncertainty Quantification

The AI system can quantify the impact of hadronic models on reconstruction by analyzing the mean differences in the reconstructed energy and lnA for EPOS-LHC and SIBYLL data relative to QGSJet-II-04, noting that these differences scale approximately linearly with lg(E/A) and are nearly independent of primary composition.