Search for neutrino emission from blazar $γ$-ray flares accounting for possible neutrino time delays

技術要約：考えられるニュートリノの時間遅延を考慮した、ブレーザーの$\gamma$線フレアからのニュートリノ放射の探索

問題提起
IceCubeニュートリノ観測所は、高エネルギーニュートリノ（$E_\nu \gtrsim 1$ TeV）の存在を確立しているが、拡散フラックスに寄与する特定のソースは、TXS 0506+056のような特定のイベントとブレイザーとの関連性が示唆されているにもかかわらず、依然として大部分が未特定である（例：NGC 1068などの活動銀河）。ブレイザーとIceCubeニュートリノとの相関に関する従来の探索では、通常、ニュートリノの到着と主要な電磁波（特に$\gamma$線）のフレアとの間の時間的一致を仮定してきた。しかし、理論的な検討によれば、この一致が成立しない可能性がある。もしニュートリノが$p\gamma$相互作用を通じて生成されるならば、ジェットフレームにおける陽子加速のタイムスケール（$t'_{acc}$）および陽子のエネルギー損失のタイムスケール（$t'_{p\gamma}$）は、電磁波フレアの原因となる電子のタイムスケールよりも数桁長くなる可能性がある。その結果、高エネルギーニュートリノ（$\gtrsim 100$ TeV）は、関連する$\gamma$線フレアのピークから数ヶ月から数年遅れて到着する場合がある。本論文は、このような潜在的なタイムラグを考慮した探索が欠如しているという問題に対処するものである。

手法
著者らは、IceCat-1公開カタログからの高エネルギーニュートリノのアラートと、2つの異なるブレイザー・ソース・サンプルとの間で、広範なジェットフレームの時間遅延（$t'_{delay}$）をスキャンしながら、相互照合分析を行った。

1. データセット:

   • ニュートリノ: IceTop信号が同時検知されていない、IceCat-1カタログからの348個のアラート（$E_\nu \gtrsim 100$ TeV）。
   • ソースリスト:
     • H21+: ドップラー因子（$D$）、赤方偏移（$z$）、およびFermi-LATライトカーブが既知であるMOJAVEデータベース（Homan et al. 2021）由来の294個のAGNサンプル。
     • R24+: スペクトルエネルギー分布モデリングから導出されたドップラー因子を持つ、CGRaBSカタログ（Rodrigues et al. 2024b）由来の295個のブレイザー・サンプル。これらはFermi-LATライトカーブとも相互照合されている。
   • ライトカーブ: Fermi-LATライトカーブ・リポジトリ（LCR）からの月単位でビン化された光子フラックス・ライトカーブ。

2. 統計的枠組み:

   • 空間的関連性: ニュートリノが拡大された角度不確定領域（系統的な不確定性 $\Delta = 0.78^\circ$ を含む）内に収まる場合、そのソースとブレイザーを関連付けた。
   • 重み付けスキーム: 各アラートに対し、アラートのシグナル性、大きな角度不確定性に対するペナルティ因子、および時間因子 $G_{\nu a}$ に基づいて、重み $w_\nu(t'_{delay})$ を算出した。
   • 時間因子: 時間因子は、ニュートリノの到着時刻（$t_\nu$）を、到着前のブレイザーの$\gamma$線ライトカーブの全域最大値（$t^{\gamma \max}_a$）に関連付けた。この因子は、時間差が予想される観測者フレームでの遅延（$t_{delay} = t'_{delay}(1+z)/D$）とビンサイズ内で一致する場合に最大重量を割り当て、偏差に対してガウス型ペナルティを適用した。
   • テスト統計量 (TS): グローバルなTSは、全アラートにおける個々の重みの総和として定義された。
   • 有意性検定: 実データのTSを、ニュートリノの赤経と到着時刻をランダム化した49,999個の模擬データセットから生成されたTSの分布と比較した。事前試行$p$値が計算され、その後、試行回数による補正（$t'_{delay}$のスキャン）が行われた。

主な貢献

• 新規な探索戦略: 本研究は、陽子の遅い加速やエネルギー損失によって引き起こされる普遍的なジェットフレームの時間遅延（$t'_{delay}$）を考慮して、IceCubeニュートリノとブレイザーの$\gamma$線フレアとの関連を探索する初めての試みである。
• 二重のソース・サンプル: ジェット運動学の不確定性に対する相関の堅牢性をテストするため、ドップラー因子の決定方法が異なる（輝度温度の議論 vs SEDモデリング）2つの独立したブレイザー・カタログを使用している。
• 系統的な遅延スキャン: ニュートリノ生成が$\gamma$線放射の後に遅れる可能性のある理論的範囲をカバーするため、$t'_{delay}$を1ヶ月から$\sim 3.5 \times 10^4$日までスキャンした。

結果

• 統計的有意性: 探索の結果、H21+またはR24+のいずれのソースリストにおいても、IceCat-1ニュートリノとAGNとの間に統計的に有意な相関を示す証拠は見出されなかった。
  • H21+リストについては、最小の事前試行$p$値は0.037であった（$t'_{delay} \approx 7.7 \times 10^3$ dにおいて）。
  • R24+リストについては、最小の事前試行$p$値は0.026であった（$t'_{delay} \approx 1.9 \times 10^3$ dにおいて）。
  • 試行補正後、これらの最小値はそれぞれ、帰無仮説と一致する事後試行$p$値0.14（1.06$\sigma$）および0.11（1.21$\sigma$）に対応した。
• トップ候補: 最も有意な関連性は、アラートIC120523B（3C 454.3の2010年のフレアに関連）およびIC150926A（3C 279のフレアに関連）によって駆動された。これらの特定のイベントは予想される遅延と時間的な整合性を示したが、データセット全体にわたるグローバルな相関を確立するには不十分であった。
• ドップラー因子の不一致: 最適な遅延時間は、H21とR24の異なる手法による3C 454.3および3C 279の推定ドップラー因子の差に大きく起因して、2つのソースリスト間で大きく異なった（4倍の差）。

意義と主張
著者らは、彼らの結果が、時間遅延を考慮した場合でも、検討された特定のブレイザー集団とIceCat-1ニュートリノとの間に有意な相関を支持しないことを控えめに結論付けている。彼らは、この結果が得られた理由として、互いに排他的ではないいくつかの可能性を挙げている：

1. 非効率性: ブレイザーは$\sim 100$ TeVのニュートリノの非効率なソースである可能性があり、IceCubeフラックスのわずかな割合（$\sim 0.5\%$）のみを寄与している。
2. エネルギー・スケール: もしブレイザーが効率的なニュートリノ源であるならば、その放射のピークは現在のIceCubeアラートの主要な感度範囲を越えた、より高いエネルギー（$\gtrsim 10$ PeV）にある可能性がある。
3. モデルの仮定: 普遍的なジェットフレームの時間遅延（$t'_{delay}$）という仮定は、加速メカニズムや磁場のソースごとの差異により、妥当ではない可能性がある。
4. デカップリング: 陽子加速とエネルギー散逸が、空間的または時間的に電子プロセスから分離している場合があり、その場合、ニュートリノ生成は観測された$\gamma$線フレアと因果関係を持たない。

本論文は、将来の感度の向上には、より大きなデータセット（次世代のIceCat-2やCAZカタログなど）の導入や、$\gamma$線フレアよりも遅れることが知られており、ニュートリノ生成のタイムスケールとより良く相関する可能性があるラジオ・ライトカーブの包含が必要になると主張している。