Null test of general relativity based on gravitational waves

技術要約：重力波に基づく一般相対性理論のヌルテスト

問題提起
一般相対性理論（GR）は、数多くの古典的および現代的なテストを通過してきたが、それが重力の最終理論ではないと広く信じられている。重力波（GW）天文学は、GRの妥当性の境界を探索するためのユニークなツールを提供する。既存のパラメータ化されたテスト（ppE、TIGER、FTIフレームワークなど）は、ポスト・ニュートン（PN）係数に偏差パラメータを導入することでGWデータを解析している。これらの研究は偏差パラメータが小さいことを示しているが、決定的な手法上の欠落が残っている。すなわち、「どの程度小さければ」GRを決定的に裏付けるのに十分であるのか、という点である。従来の手法は、偏差パラメータの大きさに依存することが多いが、系統的なバイアスやパラメータ間の相関が存在する場合に、データが真にGRを代替理論よりも支持しているかどうかを評価するための厳密な統計的枠組みを欠いている。

手法
著者らは、偏差パラメータを質量–質量関係へと変換することにより、既存のパラメータ化アプローチを拡張した「ヌルテスト」スキームを提案している。その手法は以下の通りである：

1. 波形モデリング: 本研究では、全インスパイラル–マージャー–リングダウン（IMR）過程をカバーするIMRPhenomXAS波形モデル（および特定のデータ解析にはLVKのIMRPhenomPv2）を利用する。インスパイラル相は、位相進化におけるPN係数（$\psi_j$ および $\psi_{jl}$）を修正することでパラメータ化されており、これらは標準的なソースパラメータ（全質量 $M$、対称質量比 $\eta$、スピン）と共に自由に変化させることができる。
2. 帰無仮説の構築: GRにおいて、特定のPN係数はソースパラメータ（$m_1, m_2$, スピン）によって厳密に決定される。著者らはこの関係を反転させた。すなわち、測定された偏差パラメータ（またはPN係数）の特定の値に対して、もしGRが真であるならば必要となる対応する質量–質量関係の曲線を描出する。
3. 交差解析:
   • 領域A: 標準的なGW解析から導出されたソース質量（$m_1, m_2$）の事後分布。
   • 領域B: 特定のPN係数（例：$\phi_{2,NS}, \phi_{6,NS}$）の事後分布から導出された、質量–質量図における「GRが可能な領域」。
   • テスト: もしGRが正しいならば、領域Aと領域Bは必ず交差しなければならない。著者らは、これらの領域が最初に接触するために必要な累積確率（信頼レベル）を計算する。
4. クリティカル確率: 「クリティカル確率」は、質量事後分布領域と、GRに整合する領域（一つの、あるいは複数のPN係数から導出されたもの）が共通の交差部を共有する最小の累積確率として定義される。クリティカル確率が高いことは、データのGRとの整合性が高い信頼レベルで成立していること（すなわち、テンションの確率が低いこと）を示し、クリティカル確率が低いことは、より低い信頼レベルでGRとのテンションが存在することを示唆する。
5. データ適用: このスキームを、LIGO-VIRGO-KAGRA（LVK）のO1/O2ランにおける5つの連星ブラックホール（BBH）イベント（GW150914, GW151226, GW170104, GW170608, GW170814）に適用する。本解析では、単一の係数、および2つ、3つ、4つ、5つのスピン独立なPN係数の組み合わせをテストする。

主な貢献

• 全IMRへの拡張: インスパイラル段階に限定されていた先行研究[32–34]とは異なり、本研究はヌルテストの概念を全インスパイラル–マージャー–リングダウン信号へと拡張し、インスパイラル段階の切り出しに伴う不確実性を回避するために検出データ全体を利用している。
• 質量–質量図フレームワーク: 本論文は、パルサー連星のタイミング解析テストから概念を借り、偏差パラメータの結果を質量–質量図へとマッピングする手法を導入している。これにより、異なるパラメータ間での整合性を直接的かつ視覚的に比較することが可能になる。
• 多パラメータ整合性: 著者らは、複数の係数を同時に扱うテストへとこれを一般化し、質量領域と複数の係数から導出されたGR領域が交差するクリティカル確率を定義している。
• 定量的閾値: 本研究は、単なる定性的な偏差の大きさに関する記述を超えて、「どの程度小さければ十分なのか」という問いに答えるための具体的な指標（クリティカル確率）を提供している。

結果

• イベントごとのテンション: 解析の結果、5つのイベント間でGRとのテンションの度合いが異なることが明らかになった。
  • GW150914 および GW170104: GRとのテンションが高い信頼レベル（$>90\%$）で示されており、これはデータが低い信頼レベルでしかGRと整合していない（あるいは逆に、テンションが有意である）ことを意味する。
  • GW151226, GW170608, および GW170814: GRとのテンションが非常に小さく（主に$<60\%$）、GRとの整合性が高いことを示している。
• SNRへの独立性: 本研究では、ネットワークの信号対雑音比（SNR）とGRからの偏差の強さとの間に、単調な相関は認められない。例えば、GW150914 (SNR=26.0) と GW170104 (SNR=13.8) は共に高いテンションを示したが、GW170814 (SNR=17.7) は低いテンションを示した。このことは、観察されたテンションが単に信号強度によって駆動されているのではないことを示唆している。
• 多パラメータ挙動: 多パラメータテストにおいて、クリティカル確率はしばしば最も制約の強い単一のパラメータ（質量領域と交差するために最大の累積確率を必要とするパラメータ）によって決定される。

意義および主張
本論文は、現在のGW観測データにはGRから逸脱する情報が含まれている可能性があると主張しているが、これらの偏差をGRの破綻の証拠として解釈することについては注意を促している。 著者らは、観察されたテンションが決定的なものではないことを明示している。彼らは、潜在的な偏差を以下のいくつかの系統的な要因に帰している：

1. 波形の系統誤差: マージャー–リングダウン領域におけるIMRPhenomPv2モデルの限界、および波形ファミリー間の差異。
2. 検出器のアーティファクト: 特に低SNRイベントに影響を与えるノイズのアーティファクトや較正の不確実性。
3. 事前分布の選択: PN偏差パラメータに対する広範な一様事前分布の使用が、広範な事後分布を生む可能性がある。
4. パラメータ間の相関: PN係数とソースパラメータ（スピン、質量比）の間の縮退。

本研究の意義は、決定的なGRの反証ではなく、補完的な診断ツールを提供することにある。著者らは、本ヌルテスト・スキームが、より正確な波形モデル、独立した解析パイプライン、およびより大規模なイベントカタログ（GWTC-1からGWTC-5への拡張）を用いたさらなる調査を促すものであると結論付けている。このスキームは、他の波形モデルへ容易に一般化でき、利用可能なすべてのGW観測データに適用できるように設計されている。