A stepping stone toward detecting gravitational wave memory: a cumulative analysis with the full $(\ell=2, m=0)$ spherical harmonic using events from GWTC-4.0 and GWTC-5.0

技術要約：重力波データにおける $(\ell= 2, m = 0)$ 球面調和関数モードの累積解析

問題提起
重力波（GW）データにおける光源パラメータの正確な推論には、測定可能なすべての信号寄与を捉える波形モデルが必要である。支配的な $(\ell= 2, m = \pm2)$ モードは十分に特性化されているが、非対称または高傾斜の系においては、サブドミナント（副次的）なモードが縮退を解く上で重要な役割を果たす。本論文は、特に $(2, 0)$ モードに焦点を当てている。このモードには、残存ブラックホールの準固有モードに関連する振動的なリングダウン成分と、非振動的な変位メモリーという、2つの明確に異なる物理的寄与が含まれている。変位メモリーは、重力放射の累積的な放出から生じ、時空の永続的な変化として現れる。このメモリーの検出は大きな課題であり、GWTC-1からGWTC-3までのカタログを用いたこれまでの観測解析では、決定的な検出には至っていない。理論的な予測によれば、現在の第2世代検出器による個別のイベントでの検出は困難であり、代わりに数百から数千のイベントをスタッキング（積み重ね）する必要がある。本研究は、現在のGWデータにおけるフル $(2, 0)$ モード（振動成分とメモリー成分の両方を含む）の観測的裏付けを評価することを目的としており、将来的な重力波メモリー検出への足掛かりとなるものである。

手法
著者らは、GWTC-4.0 カタログ（84イベント）および GWTC-5.0 から選定された高有意イベントのサブセット（6イベント）のバイナリブラックホール（BBH）イベントの系統的な再解析を行っている。解析には、準円軌道かつ非歳差運動のセクターに限定された、$(2, 0)$ モード（振動的リングダウンと変位メモリーの両方の寄与を含む）を含む時間領域フェノメノロジーモデルである IMRPhenomTHM_20 波形モデルを使用している。歳差運動を持つ候補 GW241127_061008 に対しては、歳差運動モデル IMRPhenomTPHM_20 が用いられている。

解析の核となるのは、2つの相補的なアプローチである：

1. ベイズモデル選択: 各イベントについて、$H_{20}$（$(2, 0)$ モードを含む IMRphenomTHM_20）と $H_{\emptyset}$（$(2, 0)$ モードを含まない IMRphenomTHM）の2つの仮説間のペア比較が行われる。相対的な証拠を定量化するためにベイズ因子（$B$）が計算される。個々の対数ベイズ因子は、統計的独立性を仮定して集団全体でスタッキングされ、累積的な証拠指標が導出される。
2. 信号対雑音比（SNR）のスタッキング: 著者らは、$(2, 0)$ モードの SNR を、フル、振動のみ、およびメモリーのみの寄与に分けて計算する。これらの SNR は、イベント集団にわたって二乗和の平方根（quadrature）でスタッキングされ、累積的な検出可能性を評価する。

解析には、測定されたパワースペクトル密度（PSD）と較正不確かさが組み込まれている。変位メモリー信号の非周期的な性質を扱うため、偽のアーティファクトを避けるためにフーリエ変換法 gw-foutstep が使用されている。また、ベイズ因子の SNR に対するスケーリングを検証し、パラメータ推定におけるバイアスを確認するために、インジェクション研究（数値相対論およびフェノメノロジー波形の両方を使用）も含まれている。

主な貢献と結果

• 個別イベントの証拠: GWTC-4.0 および選択された GWTC-5.0 イベントの大部分において、データは $(2, 0)$ モードの包含に関する統計的に有意な証拠を提供していない。ほとんどのベイズ因子はゼロ付近に集まっており、その不確かさはヌル値と重なっている。最も強い個別な選好は GW230814_230901 で見られ、対数ベイズ因子は $0.63 \pm 0.09$ である。これはジェフリーズ尺度によれば「実質的（substantial）」な証拠であるが、「決定的（decisive）」な証拠の閾値には達していない。
• 累積的証拠: GWTC-4.0 カタログの累積対数ベイズ因子は $\log_{10} B = 1.38 \pm 0.79$ と計算された。これはモードの存在に向けた緩やかな傾向を示しているものの、大きな不確かさがあるため、結果は決定的な証拠とはならない。
• SNR スタッキング: フル $(2, 0)$ モードの累積 SNR は、GWTC-4.0 サンプルに対して $3.53^{+0.94}_{-0.41}$ に達し、一般的に採用されている検出閾値である SNR = 3 を超えている。対照的に、メモリーのみの寄与は、約7.5ヶ月間の O4a 観測後でも $0.89^{+0.29}_{-0.11}$ に過ぎない。これは、振動成分を含めることが、メモリー単独で探索する場合と比較して、モードの検出能を大幅に向上させることを示唆している。
• 歳差運動系: IMRPhenomTPHM_20 を用いた歳差運動候補 GW241127_061008 の解析では、対数ベイズ因子は $-0.01 \pm 0.13$ であり、歳差運動を考慮しても $(2, 0)$ モードに対する追加の証拠は見られない。
• パラメータ推定: $(2, 0)$ モードの導入は、研究されたイベントの推定された光源パラメータ（質量、スピン、距離）に顕著なシフトをもたらしておらず、このモードの欠落が現在のパラメータ推定における主要な系統的バイアスの原因ではないことを示唆している。
• 検出の見通し: 観測されたベイズ因子と不確かさにパラメトリック分布を適合させることにより、著者らは、現在の検出器感度と使用された特定の波形モデルを仮定した場合、決定的な統計的証拠（$\log_{10} B - 1\sigma > 2$）を得るには、約 $N_{events} = 166^{+82}_{-55}$ 個のイベントが必要であると予測している。

意義と主張
本論文はこの解析を、重力波メモリー検出への「自然な足掛かり」として位置づけている。著者らは、$(2, 0)$ モードの振動成分が信号対雑音比を高めるため、フルモードの検出はメモリー成分単独の検出に先行して起こる可能性が高いと主張している。結果は、$(2, 0)$ モードに関する現在の観測的裏付けを現実的に評価しており、信号が最も明るいイベントには存在しているものの、依然としてサブドミナントであり、ロバストな検出には集団レベルのスタッキングが必要であることを確認している。

著者らは明示的に限界と控えめな主張を述べている：

• 本研究は、使用されたフェノメノロジーモデルにおける $(3, \pm2)$ モードの欠如などの波形系統論を考慮しておらず、これが（GW230814_230901 のケースで見られるように）ベイズ因子の推定に影響を与える可能性がある。
• 必要なイベント数（$N \approx 166$）の投影は、現在の解析における特定の仮定に基づいており、まだ広範な系統研究を取り入れていないため、楽観的であると考えられる。
• 本論文は、$(2, 0)$ モードが最も明るいイベントに対して測定可能な寄与をし始めているものの、堅牢なイベントレベルの検出は現在の地上ベースの感度では到達不能であり、将来の検出にはより大きなカタログとより広範な波形系統研究が必要になるであろうと結論付けている。