Improving low-latency multi-messenger follow-up of neutron star-black hole mergers with mode-by-mode filtering

技術要約：モード別フィルタリングによる中性子星–ブラックホール合体における低遅延マルチメッセンジャー追跡の改善

問題提起
中性子星–ブラックホール（NSBH）合体に関する迅速なパラメータ推定（PE）は、電磁波（EM）追跡観測の指針を決定する上で極めて重要である。BAYESTARのような現在の低遅延パイプラインは、通常、支配的な四重極重力波調和関数 $(\ell, |m|) = (2, 2)$ のみをモデル化している。NSBHのような非対称な連星系（質量比 $q \gtrsim 5\text{--}10$）において、このアプローチは、光度距離（$d_L$）、傾斜角（$\iota$）、および固有の連星パラメータ間の強い縮退を残してしまう。四重極のみの波形では、$d_L$ は $\iota$ と強く相関し、質量比は有効スピンと相関する。これらの縮退は、位置特定体積を広げ、二次天体の性質（NSかBHか）に関する制約を弱め、傾斜角の推定を劣化させる。これらは、EM観測の可能性やホスト銀河との関連付けに関する判断を妨げる要因となる。高次モード（HMs）を含む完全な確率論的サンプリングは、数時間から数日を要するため、低遅延アプリケーションには計算コストが高すぎる。

手法
著者らは、$(2, 2)$、$(3, 3)$、および $(4, 4)$ の各モードの信号対雑音比（SNR）時系列に対してモード別のフィルタリングを行う、低遅延PE手法を提案している。この手法は、cogwheel コードの低遅延版である flywheel ライブラリを利用しており、尤度に単純に寄与するパラメータ（例：$d_L$ および参照位相）については解析的または半解析的に周辺化を行い、残りの外的なパラメータについては適応的重要度サンプリングを使用する。

アルゴリズムは、主に以下の3つの入力を受け取る：

1. $\langle d | h_{\ell m} \rangle$：各モードテンプレートのSNR時系列。
2. $\langle h_{\ell m} | h_{\ell' m'} \rangle$：モード共分散行列。
3. 想定される天体物理学的分布から導出されたモードSNR比（$R_{\ell m}$）の事前分布サンプル。

本手法は、外的な連星パラメータ、およびモード振幅比に影響を与える固有パラメータのサブセットについて周辺化を行う。光度距離の積分は事前計算されたテーブルの補間によって評価され、参照位相の積分は台形公式による求積によって行われ、その他の積分（傾斜角、空の方向、偏光）は適応的重要度サンプリングを通じて行われる。事後分布サンプルは、積分の際に使用される準モンテカルロ・サンプルを重要度再サンプリングすることによって生成される。

主な結果
本手法は、2つのデータセットを用いて検証された：LIGO-VirgoネットワークにおけるO5設計感度でのシミュレーションされたNSBH検出、および過去に検出されたNSBHイベントの公開データである。

• シミュレーションされたO5検出: SNR $\ge 8$ の約 $1.6 \times 10^4$ 個のシミュレーションイベントに適用した結果、HMsの導入は $(2, 2)$ のみの解析と比較して制約を大幅に改善した：

  • 距離と体積: 光度距離の90%信用区間の幅、および位置特定の共動体積を、しばしば2倍以上の因子で減少させた。GW190814のようなイベントでは、共動体積が約4.2倍減少した。
  • 傾斜角と視線方向: HMsが $d_L$–$\iota$ の縮退を打破し、フェイスオン/オフ系の視線方向の制約を最大5倍改善した。これは、相対論的ジェットからのプロンプト・ガンマ線バースト（GRB）の検出確率に直接影響する。
  • 二次天体の質量: ソースフレームにおける二次天体の質量（$m_{src,2}$）の制約が、特に高質量系において改善された。高質量系では、二次天体が中性子星である確率（$m_{src,2} \le 2.7 M_\odot$）が最大50%増加し、キロノバ追跡戦略を支援した。
  • 計算コスト: 高次モードの周辺化には、単一のCPUコアあたり1イベントにつき $0.13 \pm 0.03$ 秒を要し、さらに1000個の事後分布サンプルを抽出するのに $0.11$ 秒を要する。これは既存の四重極のみのパイプラインと同等である。

• 実データの検証:

  • GW190814: 高SNR（$\sim 25$）、高質量比（$q \sim 9$）のイベントであるGW190814は、明確な高次モードのシグネチャーを示した。本手法を実データに適用したところ、$d_L$ の事後分布を1.6倍、傾斜角を2倍、位置特定体積を2.1倍狭め、低遅延の事後分布を完全なベイズPEの結果に近づけた。
  • その他のNSBHイベント: 他の3つの確実なNSBH検出（GW200105, GW200115, GW230529）は、より低いSNR（$\sim 11\text{--}14$）およびより低い質量を持ち、顕著な改善は見られなかった。これは、それらの副次的モード成分が弱いことと一致している。

意義と主張
本論文は、高次モードを迅速なPEに含めることで、現在の四重極のみの手法と同等の計算コストで、位置特定および二次天体質量の測定を大幅に改善できると主張している。著者らは、これらの改善が以下のEM追跡戦略において極めて重要であると述べている：

1. カタログからのホスト銀河のより効率的な特定。
2. ジェット放出（GRB/X線）の検出可能性を評価するための視線方向のより良い決定。
3. 潮汐破壊のポテンシャルに基づくキロノバ検出可能性の評価の向上。
4. 適切な追跡プロトコルを起動するための、より確信を持ったソース分類（NS vs BH）。

著者らは、本手法が今後の観測ラン（具体的にはO5）において、現在の四重極のみの低遅延探索を一般化するために展開可能な状態にあると結論付けている。また、現在の手法は非歳差運動連星を仮定しているが、歳差運動を取り入れることは将来的な拡張として自然なステップであると述べている。このアプローチは、モード別のSNR時系列とそれらの共分散のみを必要とするため、パイプラインに依存しない設計となっている。