The impact of precession and higher-order multipoles for gravitational wave cosmological inference

技術要約：重力波宇宙論的推論における歳差運動および高次多重極の影響

問題提起
重力波（GW）天文学は、「質量スペクトル（またはスペクトル・サイレン）」法を用いることで、ハッブル定数（$H_0$）を独立して測定する経路を提供します。この手法は、観測された検出器フレームでの質量と、仮定された連星ブラックホール（BBH）集団のソースフレームでの質量分布を組み合わせることで、$H_0$を推定し、電磁波による対応天体の存在を必要としません。しかし、この推論は、光源の特性（距離や質量など）を推定するための理論的な波形モデルに依存しています。現在のモデルは一般相対性理論の不完全な記述であり、具体的には、多くの解析において、スピン誘起の軌道歳差運動や高次多重極モーメントを除外した簡略化されたモデルが利用されています。これらの物理効果を含めることは、縮退（特に距離と傾斜角の間）を解き、単一イベントのパラメータ推定を改善しますが、計算コストを大幅に増大させます。ここで検討される中心的な問いは、これらの複雑な物理現象の除外が、より単純で計算コストの低いモデルを使用するメリットを上回るような、$H_0$の階層的推論における系統的なバイアスを導入するかどうかです。

手法
著者らは、異なる波形モデルから導出された宇宙論的推論を比較するために、階層的ベイズフレームワークを採用しています。本研究では、主に2つのデータセットを利用しています。

1. GWTC-3の再解析： 著者らは、第3次重力波過渡現象カタログ（GWTC-3）に含まれる42個のBBH合体を再解析しています。これには、以下の4つの異なるPhenom波形モデルを用いて生成された事後分布サンプルが使用されます。
   • XAS: スピン整列、主要な四重極のみ。
   • XP: スピンの不整合（歳差運動）、主要な四重極。
   • XHM: スピン整列、高次多重極（$\ell < 5$まで）。
   • XPHM: スピンの不整合、高次多重極（最も完全なモデル）。
2. シミュレーションされた集団： 現在の結論の限界をテストするため、著者らは、高いスピン絶対値（平均 $a_i=0.7$）と優先的な非対称質量比（$\beta_q = -1.1$）を持つ、歳差運動と高次モードの影響を最大化するように設計された「ワーストケース」の80個のBBH集団をシミュレートしました。これらの条件は、歳差運動と高次モードの影響を最大化することが期待されるものです。また、天文学的に動機付けられた低スピン集団もシミュレートしています。すべての信号は、Advanced LIGOおよびVirgoネットワークの感度を表す理想的なガウス雑音の中に注入されています。

両方のデータセットに対して、bilbyパッケージとDynestyサンプラーを用いて単一イベントのパラメータ推定が行われます。これらの単一イベントの事後分布は、その後、階層的推論を行うためにicarogwパイプラインへと送られ、ハイパーパラメータ（$H_0$、質量分布のべき乗則指数、および赤方偏移の進化を含む）が推定されます。異なるモデルから導出された事後分布間の統計的な一致は、イェンセン・シャノン・ダイバージェンス（JSD）を用いて定量化されます。

主な結果

• GWTC-3の再解析： 推定された$H_0$の事後分布は、スピン歳差運動や高次多重極の有無に対して、無視できる程度の依存性を示しました。最も複雑なモデル（XPHM）とより単純なモデル（XP、XHM）の間のJSDは、それぞれ約 $4 \times 10^{-3}$ ビットおよび $5 \times 10^{-4}$ ビットでした。これらの値は、統計的有意性の閾値（$5 \times 10^{-2}$ ビット）を大きく下回り、ランダムシードによる変動と同程度です。質量分布のパラメータ（最大ブラックホール質量 $m_{max}$ や質量比のべき乗則 $\beta_q$）において軽微な差異は見られますが、$H_0$の90%確信区間はすべてのモデル間で一貫しています。
• 高スピン集団のシミュレーション： 強い歳差運動と高次多重極の効果を示すように特別に設計された集団（単一イベント解析では距離と質量比の事後分布に顕著な差が見られる条件）を分析した場合でも、階層的推論による$H_0$は波形モデルに対してほぼ無関心（アグノスティック）でした。すべてのモデルが、真の$H_0$を$1\sigma$の確信区間内で回収しています。
  • 注：質量比について： 歳差運動を無視すると、推論される質量比分布（$\beta_q$）がより対称的な連星へと偏る（バイアスがかかる）ことが観察されました。これは、単一イベントにおける縮退と一致しています。しかし、この質量分布へのバイアスは、$H_0$の推論に重大なバイアスをもたらすことはありませんでした。
  • 注：高次モードについて： 高次多重極を無視することは、イベント数が過度に多くない限り、$H_0$および質量分布の両方に無視できる影響しか与えませんでした。
• 計算効率： 80個のシミュレートされたイベントをXAS（整列スピン、四重極のみ）で解析する計算コストは、約30 CPU日であったのに対し、XPHMでは180 CPU日でした。これは、同等の$H_0$制約を得るために、6倍の計算コスト削減が可能であることを示しています。
• 波形の系統誤差： 6つの高歳差運動信号に対するPhenom（XPHM）系とSEOBNR（SEOBNRv5PHM）系の限定的な比較では、質量と距離の周辺事後分布において極めて高い一致が見られました。これは、波形ファミリーの選択が現在の$H_0$推論における制限要因ではないことを示唆しています。

意義および主張
本論文は、現在のおよび近い将来のGWイベント数（約1000イベントまで）において、質量スペクトル法による$H_0$の推論という特定の目的のために、スピン歳差運動や高次多重極を含む計算コストの高いモデルを使用することには、大きな利点はないと結論付けています。個々のGW解析に内在する統計的不確定性が、欠落した物理学によって導入される系統誤差よりも支配的であるためです。

著者らは、単一イベントのパラメータ推定における縮退を解く（距離や傾斜角の制約を強める）ためには複雑なモデルが必要ですが、データを集計して集団レベルの$H_0$を推定する際には、「正確な単一イベントの事後分布の形状」は重要ではないと主張しています。したがって、より単純で精度の低いモデル（XASなど）を使用することで、同等の宇宙論的推論を、わずかな計算コストで得ることが可能です。この結論は、次世代の検出器が数千のイベントを検出することが予想される中で、フル物理モデルの計算負荷が深刻な問題になる可能性があることを踏まえると、特に重要です。ただし、著者らは、この結論は質量スペクトル法に特有のものであることを注記しています。もし直接的な電磁波対応天体が得られる場合（「明るいサイレン」法）、正確な距離推定のために正確なスピン歳差運動モデルが重要となります。