Energy estimation of cosmic rays by air shower radio signals

技術要約：宇宙線による空気シャワー無線信号のエネルギー推定

問題提起
宇宙線の一次エネルギーを決定することは、その起源、加速、および伝播メカニズムを理解する上で不可欠である。宇宙線によって引き起こされる広範大気シャワー（EAS）は、一次粒子の直接的な情報を破壊してしまうが、二次粒子や無線放射は再構成への経路を提供する。先行研究（例：AERA、Tunka-Rex、LOPES）では、無線信号の横分布関数（LDF）を利用し、シャワーコアからの特定の距離における信号強度と一次エネルギーを相関させてきた。しかし、これらの手法は多くの場合、各ステーションの幾何学的配置に固有の最適距離を特定する必要があり、シャワーコアの位置に敏感であるという課題がある。イランのセムナンにある4本のアンテナ無線アレイであるSURA実験は、到着方向とコア位置の再構成には成功しているが、コア位置の変動に対する感度を最小限に抑える堅牢な一次エネルギー推定手法を必要としている。

手法
著者らは、CoREASの拡張版であるCORSIKAを用いた高忠実度シミュレーションと実験的無線信号を比較することに基づく再構成手法を提案している。本手法の核心となるステップは以下の通りである：

1. シミュレーション設定:

   • SURAの設定: シミュレーションは、SURA実験の特定のパラメータ（4本のLPDA、局所磁場成分 $B_x = 39.53 \, \mu T$, $B_z = 28.1 \, \mu T$、および高度 1130 m）に一致するように行われた。
   • 高密度参照アレイ: 1メートルの間隔を持つ12,342本のアンテナからなる仮想の高密度アレイをシミュレーションし、無線信号強度の高解像度なリファレンスを生成した。このアレイは、比較のためのベースラインとして機能する。
   • イベント生成: 一次エネルギーが $10^{17}$ eV から $7 \times 10^{17}$ eV の範囲にある陽子誘起シャワーをシミュレートした。一次到着方向は、SURAのピーク検出率に基づき、天頂角 $35^\circ$、方位角 $40^\circ$ に設定された。また、傾斜の影響をテストするために、$60^\circ$ および $65^\circ$ の天頂角についても追加のシミュレーションを行った。

2. データ較正:

   • 本研究は、2段階の較正プロセス（リアルタイムFPGAフィルタリング：ハイパス、バンドパス、ノッチフィルタ、およびポストプロセッシング：到着方向の再構成、周波数応答の適用、ウェーブレットノイズ除去、および0.2 ns分解能への信号アップサンプリング）を組み込んでいる。

3. スケール因子 ($C_{ij}$) アプローチ:

   • スケール因子 $C_{ij}$ は、高密度アレイからのシミュレーション信号の平均ピーク強度 ($\bar{A}_{sim}$) と、SURA構成からの信号のピーク強度 ($\bar{A}_{exp}$) の比として定義される：
     $$C_{ij} = \frac{\bar{A}_{sim}}{\bar{A}_{exp}}$$
   • この因子は、アンテナアレイ全体にわたる実験的（シミュレーションされたSURA）信号と高密度アレイ信号を比較する $\chi^2$ 関数を最小化することによって導出される。
   • 本手法は、絶対的な信号強度はコアの位置に依存するが、疎なSURAアレイと高密度参照アレイの間の「比（スケール因子）」は、一次エネルギーと特定の相関を示すという仮定に基づいている。

主な貢献

• コア位置への独立性: 本論文は、スケール因子 $C_{ij}$ が宇宙線シャワーのコア位置にほとんど依存しないことを示している。コアの位置がアレイの中心から離れている場合（異なるテストにおいて $\pm 30$ m および $\pm 125$ m まで）でも、スケール因子は特定のエネルギーに対して事実上一定に保たれる。
• エネルギー相関: 著者らは、スケール因子と一次エネルギーの間の明確な逆相関を確立した。一次エネルギーが増加するにつれて、スケール因子は減少する。この関係は、以下の逆一次関数によってモデル化される：
  $$C_{ij} = y_0 + \frac{a}{E}$$
  ここで、$E$ は $10^{17}$ eV で正規化された一次エネルギーであり、$y_0$ と $a$ は到着方向に依存するフィッティング定数である。
• 到着方向に対する堅牢性: 本手法は、傾斜したシャワー（天頂角 $60^\circ$ および $65^\circ$）に対してもテストされた。$C_{ij}$ とエネルギーの間の逆相関は、これらの角度においても保持されるが、フィッティングパラメータ（$a$ および $y_0$）は天頂角によってわずかに変化する。

結果

• 再構成精度: 導出されたスケール因子と適合された逆関数を用いて、シミュレートされた宇宙線の一次エネルギーを再構成した。本手法は、最大で約 11% の再構成誤差を達成した（天頂角 $35^\circ$ における $5 \times 10^{17}$ eV 時）。
• 誤差解析: 相対誤差 ($\Delta E / E$) は、特定のエネルギーレベル（例：$2 \times 10^{17}$ eV）では 2% まで低く、テストされたエネルギースペクトル全体で一般に 1% ～ 11% の範囲内に収まった。
• 傾斜シャワー: より傾斜したシャワー（$60^\circ$ および $65^\circ$ の天頂角）に対しても、本手法は有効であったが、精度は変動した。最も低い精度は、$65^\circ$ の角度における $1 \times 10^{17}$ eV および $5 \times 10^{17}$ eV で観察され、$60^\circ$ の角度における最大誤差は $7 \times 10^{17}$ eV で発生した。

意義および主張
本論文は、この手法が、無線ベースのエネルギー再構成における一般的な不確実性の原因であるシャワーコア位置への感度を低減することで、従来のアプローチに対して明確な優位性を持つと主張している。疎な実験データを高密度参照アレイを用いて正規化することにより、本手法は、各イベントに対して特定の最適距離を正確に決定することを必要とせずに、信頼性の高い一次エネルギー推定値を提供する。

著者らは、この手法が垂直および傾斜したシャワーの両方に有効であり、実際のSURA実験のデータ収集が完了次第、同実験のデータに適用可能であると結論付けている。今後の課題として、この手法を実際のSURA実験データに適用すること、および、手法の堅牢性をさらに検証するために、より遠方のコア位置や多様な到着方向を持つシャワーを含むシミュレーションを拡張することが挙げられている。本研究は既存の手法を置き換えることを目的としているのではなく、SURAアレイの幾何学的制約に特化した補完的なアプローチを提示している。