[论文解读] Updated Air-Shower $X_{ m max}$ Moment Parametrizations for UHECR Composition with Latest Hadronic Interaction Models

Carmelo Evoli, Igor Vaiman+2 more|arXiv (Cornell University)|Feb 20, 2026

Astrophysics and Cosmic Phenomena被引用 0

一句话总结

论文为 Xmax 的均值和方差的更新参数化，以及一个广义 Gumbel 模型，用于描述 UHECR 大气 shower 的完整 Xmax 分布，使用三种现代高能强子相互作用模型，以改进组分分析。

ABSTRACT

The mass composition of ultra-high-energy cosmic rays (UHECRs) is commonly inferred from the first two moments of the depth of shower maximum, $X_{ m max}$, measured by fluorescence and hybrid detectors. Such analyses require fast and accurate mappings between the moments of $X_{ m max}$ and those of the logarithmic mass, $\ln A$, based on realistic air-shower simulations. In this work we provide updated parametrizations of the $X_{ m max}$ moments and distributions for air showers initiated by nuclei from proton to iron, simulated with CONEX for three state-of-the-art hadronic interaction models: Epos LHC-R, Sibyll 2.3e, and QGSJet-III-01. We parametrize the mean depth $\langle X_{ m max} angle$ and the variance $σ^2(X_{ m max})$ as functions of energy and mass. For the variance we compare a second-order polynomial model with an exponential model. In addition, we model the full $X_{ m max}$ distributions with a three-parameter generalized Gumbel function. The Gumbel parameters are fitted using an unbinned likelihood and are validated by comparing the implied mean and variance with the raw CONEX samples and with the moment parametrizations. Across the full energy range considered, residuals between the parametrizations (or the Gumbel representation) and the simulations are at the level of a few g cm$^{-2}$ for the mean and a few (g cm$^{-2}$)$^2$ for the variance, making these parametrizations suitable for precision UHECR composition studies and forward-folding analyses of $X_{ m max}$ distributions.

研究动机与目标

• 提供 Xmax 时刻与 ln A 时刻在 UHECR 成分分析中的快速、准确映射。
• 使用在多种高能相互作用模型下的 CONEX 模拟更新 Xmax 时刻参数化。
• 用广义 Gumbel 函数刻画完整的 Xmax 分布并与模拟结果进行验证。
• 评估 Xmax 观测量在模型上的差异，以量化成分推断中的系统不确定性。

提出的方法

• 对 H、He、N、Si、Fe 初等粒子在 log10(E/ eV)=17.5–20.5、天顶角 60 度的条件下运行 CONEX 大气 shower 模拟。
• 将均值 Xmax 及其方差参数化为能量和质量的函数（式 4，方差有两种形式：二阶多项式（式 5）和指数式（式 6））。
• 使用 chi-squared 最小化对每个 HIM（EPOS LHC-R、Sibyll 2.3e、QGSJet III-01）拟合参数并比较残差。
• 用三参广义 Gumbel 函数（式 20）刻画 Xmax 分布，其参数 mu、sigma、lambda 依赖能量与质量（式 21–23）。
• 进行非分箱似然拟合（式 27）以确定 21 个系数，对 mu、sigma、lambda 的物理性约束和跨模型验证（图 6–9）。
• 提供数据产物（矩量表、分布和参数文件）在 Zenodo 仓库中供公众使用。

实验结果

研究问题

• RQ1 如何将 Xmax 矩量作为能量与质量的函数在当前高能相互作用模型下进行参数化？
• RQ2 广义 Gumbel 描述对 Xmax 分布在不同初粒子类型和能量下的函数形式与精度如何？
• RQ3 不同的高能相互作用模型如何影响 Xmax 的均值和方差，对 UHECR 组分推断有何影响？
• RQ4 能否将 Xmax 矩量参数化反演得到在模型一致的方式下的 ln A 矩量？
• RQ5 这些参数化在向前卷积和跨模型跨实验比较中的鲁棒性如何？

主要发现

• 作为能量和 ln A 的函数的 Xmax 均值和 Var(Xmax) 参数化在不同模型下的残差只有数 g cm-2 的数量级（均值），以及数 (g cm-2)^2 的数量级（方差）。
• 测试了两种方差模型：二阶多项式（式 5）和指数形式（式 6）；指数形式在 EPOS LHC-R 和 Sibyll 2.3e 上通常给出更小的残差，但存在模型相关差异。
• 在初粒子类型和能量范围内，Xmax 分布可被三参广义 Gumbel 函数很好描述，CDF 的残差在几个百分点等级。
• 带能量与质量依赖的广义 Gumbel 参数 mu、sigma、lambda（式 21–23）提供一个内部一致的表示，能够通过式 29–30 重现 CONEX 矩。
• 在不同模型之间比较时，质子（Protons）的 Xmax 最深（EPOS LHC-R），铁核的排序则改变；方差存在模型相关差异，影响质量推断。
• 所有数据产物与拟合参数公开可在 Zenodo 仓库中获取，便于使用与复现。