Testing modified gravity with the eccentric neutron star--black hole merger GW200105

技術要約：離心軌道を持つ中性子星–ブラックホール合体 GW200105 による修正重力理論の検証

問題提起
重力波（GW）検出は、強重力場かつ高度に動的な領域における一般相対性理論（GR）をテストするためのユニークな実験室を提供している。しかし、現在のテストは通常、準円軌道を仮定している。最近の解析において、中性子星–ブラックホール（NSBH）合体 GW200105 の解析から、軌道の離心率に関する強い証拠が特定された。著者らは、GR をテストするために使用される波形モデルにおいてこの離心率を無視することは、「偽陽性」、すなわち、モデル化されていない軌道形態に起因する、GR からの（見かけ上の）逸脱を引き起こす可能性があると主張している。さらに、離心率によって誘発される複雑な波形構造は、複数のフーリエ調和関数を励起するため、準円モデルよりも高い精度で修正重力理論を探索できる可能性を秘めている。

手法
著者らは、離心・歳差運動を行う波形モデル pyEFPE を拡張し、以下の3つの特定の修正重力理論によって誘発される主要なポスト・ニュートン（PN）補正を組み込んだ。

1. ブランズ・ディッケ（BD）重力： スカラー場がリッチスカラーに結合する。
2. アインシュタイン・ディラトン・ガウス・ボネ（EdGB）重力： 無質量スカラー場がガウス・ボネ項に結合する。
3. 動的チャーン・サイモンズ（dCS）重力： 無質量スカラー場がポントリャギン密度に結合する。

本研究では、準ケプラー（QK）形式とPN展開を用いて離心バイナリの動力学を記述する pyEFPE フレームワークを利用している。著者らは、エネルギーおよび角運動量の損失に関する進化方程式に、主要な次数の GR 超の補正を組み込んだ。具体的には以下の通りである：

• BD および EdGB 理論については、補正は同じ PN 次数で現れ、コンパクト天体の感度に依存する双極子放射項を介して、PN 展開パラメータ（$y$）および離心率（$e^2$）の変化率を修正する。
• dCS 重力については、補正はより高い PN 次数（GR に対して 2PN）で現れ、スピン・スピンおよびモノポール・四重極結合によってソースされるため、バイナリのスピン成分に依存する。

著者らは、GWTC-3 カタログに含まれる公開されている GW200105 のストレイン・データに対し、Bilby ライブラリと Dynesty ネステッド・サンプリングを用いて、ベイズ・パラメータ推定を実施した。彼らは以下の2つのシナリオを比較している：

1. 離心モデル： 修正重力補正を含む、完全な離心・歳差運動波形を使用。
2. 準円モデル： 同様の修正重力拡張を備えた標準的な準円波形（IMRPhenomXPHM）を使用。

知見を検証するため、著者らは pyEFPE で生成された GW200105 に似た離心信号を用い、準円モデルで回収するインジェクション・リカバリー研究を行った。

主な結果

• GR テストへの離心率の影響： 準円モデルを用いて GW200105 を解析すると、PN 係数のパラメータ化された偏差に関する周辺化された事後分布は、GR からの明らかな偏差を示した（ほとんどの PN 次数において 90% 信頼区間から外れている）。しかし、インジェクション研究によれば、これらの偏差は、離心率がモデル化されていない場合に真の GR 信号を回収した際の結果と一致している。これは、この偏差が GR の破れではなく、波形モデルによるアーティファクトであることを裏付けている。
• EdGB 重力に対する制約： 離心率を含めることで、EdGB 重力の境界が大幅に厳しくなる。結合パラメータの 90% 上限は、準円モデルの $\lesssim 5.4$ km に対し、離心モデルでは $\alpha^{1/2}_{\text{EdGB}} \lesssim 2.38$ km である。これは約 2.3 倍の改善を示す。
• Brans-Dicke 重力に対する制約： 同様に、Brans-Dicke 結合パラメータの下限は、$\omega_{\text{BD}} \gtrsim 1.78$（準円）から $\omega_{\text{BD}} \gtrsim 3.5$（離心）へと改善し、約 2 倍の改善が見られた。著者らは、事後分布における GR の明らかな違反（$\omega_{\text{BD}} \to \infty$ における支持の欠如）は、事前分布に起因するアーティファクトであり、対数一様事前分布を使用することで消失することを指摘している。
• dCS 重力に対する制約： dCS 重力については、有意な制約もバイアスも得られなかった（両モデルとも $\alpha^{1/2}_{\text{dCS}} \lesssim 42$ km）。これは、GW200105 システムのスピン成分が低いことに起因する。dCS 補正はスピンによってソースされるため、有意なスピンの欠如により、このシステムは当該理論に対して感度が低い。

意義および主張
本論文は、離心バイナリの GW 信号を用いて修正重力理論をテストするベイズ解析を行った最初の研究であると主張している。主要な意義は、軌道の離心率をモデリングすることが、重力テストの正確性において極めて重要であることを示した点にある。著者らは、離心率を考慮しないことは、GR からの偏差に関する誤った結論を導くだけでなく、代替理論に対する感度を低下させることも示している。離心率を組み込むことで、本研究は EdGB および BD 理論に対して大幅に厳格な制約を導き出した。

著者らは dCS の制約については謙虚な姿勢を保ち、システムの低スピンが意味のある境界設定を妨げていることを認めている。彼らは、検出器の感度が向上するにつれ（例：アインシュタイン・テレスコープのような第三世代検出器）、残留離心率の検出可能性が高まり、準円バイナリから得られるものよりも数桁強力な制約が得られる可能性があると結論づけている。本研究は、準円軌道を仮定する現在の GR テスト・パイプラインにおける潜在的な落とし穴を浮き彫りにしており、将来の離心イベントの解析には、系統的なバイアスを避けるために離心補正を含める必要があることを示唆している。