Chasing the neutrino blazar candidates II: SED modeling with hadronic model

技術的サマリー：ニュートリノ・ブレーザー候補の追跡 II：ハドロンモデルによる SED モデリング

問題提起
ブレーザーは高エネルギーニュートリノ源の有望な候補と考えられているが、そのニュートリノ放射の物理的起源は依然として不確かである。レプトンシナリオ（逆コンプトン散乱）はブレーザーのスペクトルエネルギー分布（SED）をモデル化する標準的な枠組みであるが、高エネルギーニュートリノの検出は、相対論的ハドロンが重要な役割を果たしていることを示唆している。著者らの先行研究（論文 I）では、少数ショット学習アプローチを用いて 199 のニュートリノ・ブレーザー候補（NBC）を同定した。本研究は、ハドロン的枠組みにおける広帯域 SED モデリングを通じてこれらの候補の物理的特性を調査し、特に最大可能なニュートリノ出力を推定するとともに、ブレーザージェット内の粒子組成を制約することで、その研究を拡張することを目的としている。

手法
著者らは、十分なマルチウェーブ長カバレッジに基づいて元の 199 個から選定された 103 個の NBC と、確認済みの源 TXS 0506+056 について、広帯域 SED を構築した。モデリングには、以下の主要な仮定と手順を用いた 1 領域・時間依存ハドロンコード（AM3）が採用された：

• ハドロン優位性： ニュートリノ放射の上限を推定するため、モデルは高エネルギー（HE）および極高エネルギー（VHE）光子が主に $p\gamma$ 相互作用および関連する電磁カスケードを通じて生成されると仮定する。レプトン的逆コンプトン（IC）散乱からの寄与は強く抑制される。
• 物理的設定： 放射領域は、磁場 $B$ 中をバルクローレンツ因子 $\Gamma$ で運動する球状のブロブとしてモデル化される。非熱的電子と陽子がべき乗スペクトルで注入される。
• モデル化された過程： このモデルは、シンクロトロン放射、光パイオン生成（$p\gamma \to n\pi^+$ または $p\pi^0$）、ベッテ・ハイテラー対生成（$p\gamma \to pe^+e^-$）、光子 - 光子対生成、およびニュートリノ生成のためのパイオンとミューオンの崩壊を考慮する。高エネルギーガンマ線には銀河系外背景光（EBL）による減衰が適用される。
• フィッティング戦略： 著者らは PySwarms パッケージを介した粒子群最適化（PSO）アルゴリズムを用い、ラジオからガンマ線帯にわたるマルチウェーブ長データに対して、磁場強度、領域サイズ、粒子注入光度、スペクトル指数などの 9 つの自由パラメータをフィッティングした。1 領域モデルにおけるシンクロトロン自己吸収の制限により、$10^{11}$ Hz 以下のラジオデータは除外された。許容可能な SED 表現を確保するため、「目視によるフィッティング」による精査が適用された。
• 相関分析： 赤方偏移依存性を制御しつつ、ニュートリノ光度（$L_\nu$）と電磁気的光度（ラジオ、光学、X 線、ガンマ線）および中心エンジン特性（ブラックホール質量 $M_{BH}$ と降着円盤光度 $L_{disk}$）の間の関係を調査するために、部分相関分析が実施された。

主要な貢献と結果

1. SED モデリングとパラメータ制約： 本研究は 103 個の NBC に対して 9 つの主要パラメータを成功裡に制約した。モデリングは、磁場強度が高くなる傾向がある（66% の源で $\log B \ge 1.0$）ことを明らかにしており、これは電子エネルギーをシンクロトロン放射を通じて効率的に散逸させ、IC 放射を抑制して高エネルギーのヒルムを陽子が支配することを可能にするものである。相対論的陽子が運ぶエネルギー（$P_p$）は、電子のエネルギー（$P_e$）および磁場のエネルギー（$P_B$）の両方を凌駕しており、平均比は $\log(P_p/P_B) \simeq 1.67$ である。
2. MeV 帯の予測： このモデルの顕著な特徴は、MeV 帯でピークを持つ陽子シンクロトロン放射の予測である。103 個の NBC のうち 99 個において、陽子シンクロトロンピークは 0.1–100 MeV の範囲内に位置する。これらの 90 個の源において、MeV 帯の陽子シンクロトロンフラックスは、GeV 帯のベッテ・ハイテラー対生成フラックスを上回る。
3. 電磁放射との相関：
   • ガンマ線と X 線： ニュートリノ光度とガンマ線（$r=0.78$）および X 線（$r=0.52$）光度の間には、強い正の相関が見出された。これは、高エネルギー陽子がニュートリノと高エネルギー光子の両方を生成するハドロン起源を支持するものである。
   • 光学（R 帯）： ニュートリノと光学 R 帯光度の間には、弱から中程度の相関（$r=0.30$）が存在する。著者らは、これがジェットからの光学シンクロトロン光子が光パイオン生成のための種光子として機能するため、あるいは降着円盤との関連によるものである可能性を指摘している。
   • ラジオ： ニュートリノと 1.4 GHz ラジオ光度の間には、有意な相関は見出されなかった。
   • 中心エンジン： ニュートリノ光度とブラックホール質量（$M_{BH}$）または降着円盤光度（$L_{disk}$）の間には統計的に有意な相関は見出されなかったが、$L_{disk}$ に対しては視覚的に弱い正の傾向が認められた。
4. 検出の可能性： モデルによって予測された最大ニュートリノフラックスに基づくと：
   • IceCube： 3 個の NBC が潜在的に検出可能である。
   • 将来の施設： 最大 22 個の源が KM3NeT によって、45 個が NEON によって、62 個が TRIDENT によって検出可能である可能性がある。
   • サンプルからの合計ニュートリノフラックスは、現在の IceCube のスタッキング限界を超えないため、「最大ニュートリノ」シナリオは物理的に妥当であることが示唆される。

意義と主張
本論文は、その結果が将来のマルチメッセンジャー観測に対する検証可能な予測を提供し、ブレーザージェット組成に関する新たな洞察をもたらすと主張している。具体的には：

• MeV ウィンドウ： 著者らは、ハドロンモデルが標準的なレプトンモデルには存在しない明確な陽子シンクロトロンピークをこの範囲で予測するため、レプトンシナリオとハドロンシナリオを区別するための重要な観測ウィンドウとして MeV 帯を強調している。
• ジェット組成： モデリングは、これらの候補におけるジェットエネルギー収支を相対論的陽子が支配するというシナリオを支持しており、TXS 0506+056 からのマルチメッセンジャー証拠と整合的である。
• 上限： 結果はニュートリノ放射の上限を表している。実際のフラックスはより低い可能性があり、これがモデルにおいて現在の感度曲線以上にあるにもかかわらず、多くの候補がまだニュートリノ放射源として確認されていない理由を説明する。
• 将来の観測： 本研究は、これらのハドロンシナリオを検証し、ブレーザージェット内の粒子組成を制約するために、次世代 MeV ガンマ線ミッション（e-ASTROGAM、AMEGO など）および次世代ニュートリノ検出器（KM3NeT、NEON、TRIDENT）が不可欠であると強調している。

著者らは、光学 - ニュートリノ相関の起源については慎重な姿勢を維持しており、これはサンプル選択効果やジェット対円盤といった複数のメカニズムに起因する可能性があると指摘し、一部の源において超エディントン・ジェットパワーが必要となる要件については、完全な整合性を得るために多領域モデルが必要となる可能性を認めている。