[論文レビュー] 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の平均と分散の更新されたパラメトリゼーションと、UHECR大気 shower の完全な Xmax 分布を説明する一般化 Gumbel モデルを提供し、3つの現代のハドロン相互作用モデルを用いて組成分析を改善する。

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 大気ショーアーを実行する。
• エネルギーと質量の関数として平均 Xmax とその分散をパラメータ化する（式 4、分散の2つのオプション：多項式式 Eq. 5 と指数式 Eq. 6）。
• 各 HIM（EPOS LHC-R、Sibyll 2.3e、QGSJet III-01）について最小二乗適合でパラメータをフィットし、残差を比較する。
• エネルギーと質量に依存する μ、σ、λ の3-パラメータの一般化 Gumbel 関数（式 20）で Xmax 分布をモデル化する（式 21–23）。
• 未ビン化尤度フィットを実施して（式 27）、21個の係数を決定し、 μ、σ、λ の物理性制約とモデル間検証（図 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 の残差を達成。
• 2つの分散モデルを検討：二次多項式（Eq. 5）と指数形式（Eq. 6）；指数形式は一般に EPOS LHC-R および Sibyll 2.3e に対して残差が小さくなるが、モデル依存の差異がある。
• 原始種とエネルギーに対する Xmax 分布は三パラメータの一般化 Gumbel 関数で良く説明され、CDF の残差は数％レベル。
• エネルギーと質量に依存する μ、σ、λ の一般化 Gumbel パラメータ（式 21–23）は、式 29–30 によって CONEX モーメントを再現する内部的に一貫した表現を提供。
• モデル間で比較すると、平均 Xmax は EPOS LHC-R で陽子が最も深く、鉄で順序が変化する；分散にはモデル依存の差が見られ、質量推定に影響する。
• すべてのデータ製品とフィットパラメータは、利便性と再現性のために Zenodo リポジトリで公開提供されている。