Combined KM3NeT/ARCA and ANTARES search for compact neutrino sources

技術要約：KM3NeT/ARCAとANTARESによるコンパクト・ニュートリノ源の結合探索

問題と目的
本研究は、点状の天体起源である高エネルギー・ニュートリノの探索に取り組むものである。主な目的は、2つの異なるニュートリノ検出器、すなわちANTARES望遠鏡（2008年から2022年までフランスのトゥーロン沖で稼働）とKM3NeT/ARCA検出器（現在シチリア沖で建設中）のデータを結合する可能性を示すことである。両装置は、南天および銀河系中心部の最適な観測視野を提供している。本解析は、スタンドアロンの解析と比較して感度と発見能力を向上させるために、ビン化された尤度法（binned likelihood approach）を用いた結合探索の枠組みを確立することを目的としている。

手法
解析では、以下の4つの異なる成分からなる結合データセットを利用する。

1. ANTARES トラック・データセット： 4541日の生存時間（livetime）。非重複のトラックおよびシャワーのデータセットに分割される。
2. ANTARES シャワー・データセット： 4541日の分割に含まれる。
3. ARCA19 トラック： 19本の検出ライン（48.4日分）を用いたKM3NeT/ARCAのトラック型イベント。
4. ARCA21 トラック： 21本の検出ライン（287.4日分）を用いたKM3NeT/ARCAのトラック型イベント。

定義： 「トラック（Tracks）」とはミューオン生成を伴うイベントを指し、「シャワー（Showers）」は全ニュートリノ・フレーバーによる中性流相互作用、および電子ニュートリノによる荷電流相互作用によって誘発される。

統計的枠組み
解析には、ビン化された尤度法を採用している。観測されたイベントのヒストグラム（$N$）は、基準フラックス（$\Phi_0 = 10^{-8} \text{ GeV cm}^{-2} \text{ s}^{-1}$）に基づくスケーラブルな信号推定値（$S$）と、背景事象（$B$）と比較される。尤度関数は以下のように定義される：
$$\ln L = \sum_{i}^{N_{bins}} N_i \ln(B_i + \zeta S_i) - (B_i + \zeta S_i)$$
ここで、$\zeta$ は信号強度パラメータである。検定統計量 $\lambda$ は、最良適合信号強度（$\hat{\zeta}$）と背景のみの仮説（$\zeta=0$）の間の対数尤度比である。

• 発見ポテンシャル： 50%の擬似実験において$5\sigma$の超過を達成するために必要な信号強度（$\zeta_{5\sigma}$）として定義される。
• 感度： ネイマン法を用いて計算された、$\zeta$ に関する中央値の90%信頼限界（CL）上限として定義される。

背景モデル化と系統誤差
背景モデリングに関する系統誤差に対処するため、赤緯帯にわたって4つの独立した解析ケースを実施した。

• ケース1： ARCAの背景は完全にデータ駆動型であり、ANTARESの生成はデータ駆動型、評価はモンテカルロ（MC）を用いる。
• ケース2： ケース1と同様であるが、ANTARESの生成背景データをランダム化（一様な方位角分布）している。
• ケース3： すべてのデータセットにおいて、生成と評価の両方で背景推定にDEFTアルゴリズムを使用する。
• ケース4： ARCAはデータ駆動型であり、ANTARESは完全にMCベースである。

これら4つのケースにおける結果の広がりが、系統誤差の範囲を定義する。ヒストグラムは、角度距離（$\log_{10} \alpha$）、再構成エネルギー（$\log_{10} E$）、およびANTARESトラックについては角度誤差の推定値（$\log_{10} \beta$）の観点で構築される。

主要な結果
解析では、$E^{-2}$ エネルギー・スペクトルと1:1:1 のフレーバー比を仮定している。

• 性能向上： ANTARESのデータが結合解析の結果を大きく左右しているが、ARCAデータとの結合により、ANTARES単独の解析と比較して性能が約10%向上する。
• 候補源 MG3 J225517+2409： このソースは、$(1.2 \pm 0.4) \cdot 10^{-4}$ という最小の事前確率（pre-trial p-value）を示し、これは $(3.69 \pm 0.09)\sigma$ の有意度に相当する。
  • 最良適合フラックス：$E^2\phi = (1.39 \pm 0.02) \cdot 10^{-8} \text{ GeV cm}^{-2} \text{ s}^{-1}$。
  • 事後確率（Post-trial probability）：二項公式により $0.013 \pm 0.004$ ($(2.26 \pm 0.13)\sigma$) と算出された。シミュレーションに基づくアプローチでは、事後p値は $0.018 \pm 0.006$ ($(2.12 \pm 0.14)\sigma$) であった。
  • 1255個の天体からの初期選択とそれに続くカタログ削減を考慮した簡略化された事後推定では、このような候補が観測される確率は $\approx 0.1$ である。
• その他の候補： 解析では、3C403 ($p \approx 1.5 \cdot 10^{-3}$) や TXS 0506+056 を含む既知の他のソースも評価し、106個のカタログに対して90% CL 上限を提示した。
• スカイマップ： MG3 J225517+2409 付近のスカイマップは、1つのARCA21イベントがソースの位置を含む誤差ボックスを持っていることを示しており、その信号尤度は背景のみのシナリオよりも6.3倍高い。

意義と主張
本論文は、ANTARESとKM3NeT/ARCAの間の結合解析の実現可能性とポテンシャルを示すことを主張している。現在の結果は成熟したANTARESのデータセットに支配されているものの、ARCAデータの統合が感度の測定可能な向上をもたらすことを確立している。著者らは、KM3NeT/ARCAの規模と露出が増加し続けるにつれ、数年以内にANTARESの性能レベルに達すると予想され、将来の発見において結合データ解析がますます重要になることを指摘している。本研究は、複数の背景モデリング戦略を通じて系統誤差を扱う枠組みを提供し、幅広い天体候補に対する上限値を設定している。