Evaluating the Contribution of Active Galactic Nuclei to the Diffuse High-Energy Neutrino Flux

技術要約：高エネルギーニュートリノの拡散フラックスに対する活動銀河核（AGN）の寄与の評価

問題提起
2013年のIceCubeによる高エネルギー天体物理学的ニュートリノの拡散フラックスの検出に続き、このフラックスの起源は高エネルギー天体物理学における主要な疑問であり続けている。ブレイザーであるTXS 0506+56やセイファート銀河NGC 1068といった特定の源は特定されているが、拡散フラックスの大部分はいまだ未同定である。活動銀河核（AGN）は、宇宙線の加速とニュートリノ生成の有望な候補である。しかし、より小さなデータセット（1〜3年）を用いた先行研究では、ガンマ線が明るいAGNがこのフラックスに大きく寄与している証拠は見つからず、これはニュートリノ生成がガンマ線が遮蔽される光学的に厚い環境で起こっている可能性を示唆している。本研究は、10年間のIceCubeデータを用い、様々なAGN集団（具体的にはX線が明るいものおよびガンマ線が明るいブレイザーおよび非ブレイザー）の拡散ニュートリノフラックスへの寄与を体系的に評価することを目的とする。

手法
著者らは、10年間の公開されているIceCubeミューオン・トラック・ニュートリノ候補イベント（2008年〜2018年）を解析した。解析には、ソース間のニュートリノ放射の差異を明示的に許容することで、従来のスタック型頻度論的探索とは異なる、基礎となるソース集団の特性を推論するためのモデルベースの統計的枠組みを用いた。

• ソースカタログ:
  • X線が明るいAGN: 70ヶ月間のSwift-BAT硬X線サーベイに基づく732個の非ブレイザーAGN。
  • ガンマ線が明るいAGN: Fermi 4LAC-DR3カタログに基づく3,339個のブレイザーおよび64個の非ブレイザーAGN。非変動ソースのサブセットも解析対象とした。
• 仮説: 電磁気学的特性に基づいてニュートリノフラックス（$F_\nu$）を予測するために、3つのスケーリング仮説をテストした。
  1. X線スケーリング: $F_\nu \propto$ 内在的X線フラックス（軟X線および硬X線バンド）。
  2. ガンマ線スケーリング: $F_\nu \propto$ ガンマ線フラックス。
  3. 幾何学的スケーリング: $F_\nu \propto 1/D_L^2$ （光度距離の逆二乗）。
• 統計的枠組み: ソース集団からの総イベント数に関する尤度関数を構築した。決定的な違いとして、本モデルはソース間のフラックスの変動（$\delta$）を考慮するために、$\delta$を自由パラメータとして扱う対数正規確率密度関数を組み込んでいる。テスト統計量（TS）は、帰無仮説（信号なし）に対する対数尤度の変化の2倍として定義された。
• パラメータ: ニュートリノスペクトル指数（$\alpha$）を2.5および3.0とし、変動パラメータ$\delta$を0から2.5の範囲で変化させた。カタログ化されたサンプルからソースクラス全体へと制約を外挿するために、完備性係数が適用された。

主な結果

• ガンマ線が明るいブレイザー: 解析の結果、ガンマ線が明るいブレイザーからのニュートリノ放射の証拠は見られなかった。本研究は厳格な上限値を設定しており、スペクトル指数および変動性の仮定に応じて、これらのソースが全拡散ニュートリノフラックスの14.5%から15.6%（2$\sigma$信頼区間）以上を占めることはできないと結論付けている。
• ガンマ線が明るい非ブレイザーAGN: 同様に、ガンマ線が明るい非ブレイザーAGNからのニュートリノ放射に関する統計的に有意な証拠は見出されなかった。本研究は、有意な寄与の可能性を排除できないが、正の信号は検出されなかった。
• X線が明るい非ブレイザーAGN (Swift-BAT):
  • 有意性: この集団に対して、ニュートリノ放射に対する統計的に有意な選好が見られた。最も高い有意性は、スペクトル指数 $\alpha=3.0$ における硬X線スケーリングで4.18$\sigma$（3.98$\sigma$）であった。
  • NGC 1068による支配: この統計的有意性は、近傍のセイファート銀河NGC 1068によって強く駆動されている。NGC 1068を除外すると、有意性は1.35$\sigma$に低下する。
  • フラックスの寄与: NGC 1068が支配的であるものの、このクラスの総寄与に関する上限値は堅牢である。著者らは、X線が明るい非ブレイザーAGNが、2$\sigma$信頼区間で拡散フラックスの**11.2%から100%**の間を占めると推定している。
  • 個別のソース: NGC 1068（4.9$\sigma$）に加え、解析ではいくつかの近傍のX線が明るいセイファート型AGNからのニュートリノ放射のプレトリアル証拠を特定した。これには、SWIFT J1041.4-1740（2.6$\sigma$）、SWIFT J0202.4+6824A/B（2.6$\sigma$）、SWIFT J0744.0+2914（2.6$\sigma$）、NGC 4151（2.5$\sigma$）、およびNGC 3079（2.5$\sigma$）が含まれる。
  • 相関: IceCubeニュートリノとSwift-BATカタログとの間に4.2$\sigma$の相関が特定されたが、これはNGC 1068に支配されている。
• ソース間の変動: Swift-BAT集団の最良適合モデルは、大きなソース間変動（$\delta \sim 1.7$）を示しており、これはニュートリノフラックスが平均的にはX線フラックスに比例するものの、個々のソース（NGC 1068のような）は平均的な予測から大きく逸脱（因子 $\sim e^{2\delta} \approx 30$）し得ることを示唆している。

意義と主張
本論文は、ガンマ線が明るいブレイザーは高エネルギーニュートリノの拡散フラックスの主要な源ではないことが否定された一方で、X線が明るくガンマ線が遮蔽されたAGN（特にセイファート銀河）が有力な候補クラスであることを結論付けている。NGC 1068からのニュートリノ検出およびSwift-BATカタログとの相関は、これらのシステムにおけるニュートリノ生成が、ガンマ線が吸収される超巨大ブラックホール周囲のコロナのような、光学的に厚い環境で起こっていることを示唆している。

著者らは、これらの結果が「ガンマ線遮蔽AGNが拡散フラックスの重要な寄与者であるという新たな議論を強化する」と述べている。彼らは、個々のソースに関する証拠（トライアル係数を克服できていない点）は決定的ではないものの、近傍のX線が明るいセイファート型AGNにおける高有意なソースの集まりは、より大きなソース集団を示唆していると強調している。本研究は、拡散フラックスの起源を解明したと主張しているのではなく、実行可能なソースクラスを、遮蔽された非ブレイザーAGNへと大幅に絞り込んだものである。論文は、これらのソースを決定的に特定するためには、将来の大型体積ニュートリノ検出器（例：IceCube-Gen2）が必要になる可能性があると断じている。