The Galactic Neutrino Sky: Predictions from Gamma-ray Source Populations

技術要約：銀河系ニュートリノの空：ガンマ線源集団からの予測

問題提起
IceCubeによって、銀河面（GP）からの高エネルギーニュートリノ放射が5.7$\sigma$の有意性で検出されており、銀河中心部に向かって顕著な過剰が示されている。総フラックスは、銀河拡散$\gamma$線放射（GDGE）から導出されるマルチメッセンジャー予測とおおむね一致しているが、信号の空間的形態は、離散的なソースからの寄与を示唆している。従来の解析では、主に（GDNEを追跡するように設計された）Fermi-$\pi^0$やKRA$\gamma$モデルといった、拡散成分のみの空間テンプレートが使用されてきた。しかし、これらの拡散のみのモデルは、観測された銀河中心部の過剰に適合させるために、再正規化を必要とすることが多い。対処される中心的な問題は、観測されたニュートリノ強度プロファイル、特に銀河中心部へのピークが、アドホックな拡散成分の再正規化を行うことなく、分解および未分解のハドロン的$\gamma$線源の空間分布によって説明できるかどうかである。

手法
著者らは、分解された$\gamma$線源、拡散放射、および未分解のソース集団の3つの成分を統合することにより、包括的な銀河ニュートリノ放射モデルを構築している。

1. 分解されたソース・テンプレート（ReGal-$\gamma$ および CTA）：

   • ReGal-$\gamma$： GeVからPeVまでのエネルギー領域にわたる4つの主要な$\gamma$線カタログ（Fermi-LAT 4FGL-DR4、H.E.S.S. Galactic Plane Survey (HGPS)、LHAASO 1LHAASO、およびHAWC 4HWC）を組み合わせた主要なテンプレート。著者らは、その$\gamma$線放射が主にレプトン的であると仮定し、銀河外ソース、パルサー（PSR）、パルサー風星雲（PWN）、およびTeVヘイローを除外している。
   • CTA テンプレート： チェレンコフ・テレスコープ・アレイ観測装置（CTAO）およびSouthern Wide-field Gamma-ray Observatory (SWGO) のシミュレーションに使用されるカタログ（gamma-cat, 3FHL, 2HWC, HGPS）に基づき、PSR、PWN、および高緯度ソースを除外して作成された独立したテンプレート。
   • ニュートリノへの変換： 各ソースについて、$\gamma$線のハドロン分率 $\chi$（$\gamma$線のうちハドロン相互作用に由来する割合）を仮定して、$\gamma$線フラックスをニュートリノフラックスに変換する。ニュートリノスペクトルは、$E_\nu^2 dN_\nu/dE_\nu \approx \frac{\chi}{2} (E_\gamma^2 dN_\gamma/dE_\gamma)|_{E_\gamma=2E_\nu}$ という関係式を用いて導出される。パラメータ $\chi$ は自由パラメータとして扱われ、テストのベースラインとして $\chi=1$ が使用される。

2. 拡散および未分解成分：

   • 拡散放射： GDNEテンプレートとしてCRINGEモデルが採用されている。これは、局所的なデータにフィットさせた宇宙線輸送方程式を解くことで導出されたものである。
   • 未分解ソース： 現在の検出閾値以下のソースを考慮するために、2つのモデルがテストされている：
     • モデルA： SNR、PWN、および$\gamma$線連星をモデリングした、Abe et al. (2024) に基づく合成集団（PWNは除外）。
     • モデルB： Schwefer et al. (2023) に基づき、PWN的な寄与の除去を近似するために全フラックスの50%にスケーリングしたもの。

3. 観測との比較：
   予測されるニュートリノフラックスは、IceCubeの有効面積（$A_{eff}$）および生存時間（livetime）と畳み込まれ、期待イベント数が生成される。解析では、結合モデル（分解 ＋ 未分解 ＋ 拡散）の経度および緯度プロファイルを、IceCubeの12年間のカスケードおよびトラック・データの、特に銀河中心部（$|l| \le 20^\circ, |b| \le 15^\circ$）と比較する。

主な貢献

• ReGal-$\gamma$ の構築： ハドロン的ニュートリノ放射の候補を特定するために特別に調整された、分解された銀河$\gamma$線源（GeV–PeV）の統一的なマルチ波長テンプレートの作成。
• 空間形態の解析： $\gamma$線カタログから推測されるニュートリノ源の空間分布が、標準的な拡散モデルとは異なる明確な形態学的シグネチャーを提供することを実証。
• ロバスト性の検証： 特定のカタログ選択によるアーティファクトではないことを保証するため、2つの独立したソーステンプレート（ReGal-$\gamma$ と CTA）および2つの異なる未分解ソースモデルを用いて結果を検証。

結果

• 空間的集中： ReGal-$\gamma$ テンプレートは、拡散のみのモデル（CRINGE, Fermi-$\pi^0$, KRA$\gamma$）と比較して、ニュートリノ強度が銀河中心部により強く集中することを予測している。これは、分解されたハドロン的ソースの密度が銀河中心部でより高いためである。
• スペクトルおよび経度フィット： 再正規化されていないCRINGE拡散モデルを、分解および未分解のソース集団（ReGal-$\gamma$ または CTA）と組み合わせたとき、モデルはIceCubeが報告しているスペクトルエネルギー分布（SED）および経度カウントプロファイルを正常に再現する。
• 再正規化の不要性： 拡散テンプレートのみを用いたフィットとは異なり、結合モデルは、放射モデルの追加的な再正規化を必要とせずに、観測された銀河中心部の過剰を説明できる。データは、平均的なハドロン分率 $\chi \sim 0.5–1$ と一致している。
• エネルギー依存性： 高エネルギー（例：100 TeV）において、ReGal-$\gamma$ テンプレートは、南半球におけるエアシャワー観測所の欠如による南天のソースの希少性を反映している。一方、スペクトル外挿に依存するCTA テンプレートは、より大きな不確実性を伴う。

意義
本論文は、銀河中心部における顕著なニュートリノ強度のピークは、拡散的な宇宙線相互作用のみではなく、離散的なハドロン的ソースの集団によって引き起こされている可能性が高いと主張している。研究結果は、従来の拡散のみの解析で見られた「過剰」が、正規化の問題ではなく、形態学的なミスマッチであったことを示唆している。観測されたニュートリノ信号を、分解および未分解の$\gamma$線源によって説明できることを確立することで、本研究は、銀河ニュートリノの空に関する自己整合的なマルチメッセンジャー・フレームワークを提供している。著者らは、エネルギー依存の経度および緯度ニュートリノ分布のさらなる測定が、これらのモデルをさらに制約し、銀河ニュートリノ源の真の空間分布を決定するために極めて重要であると結論付けている。