Time-domain phenomenological multipolar waveforms for aligned-spin binary black holes in elliptical orbits

技術要約：整列スピンを持つ楕円軌道連星ブラックホールに対する時間領域フェノメノロジカル多重極波形

問題の所在
地上設置型の重力波検出器は、数多くのコンパクト連星合体を特定してきたが、そのほとんどは孤立連星進化によって形成される準円軌道（QC）と一致している。しかし、高密度天体環境における力学的相互作用やコザイ・リドフ機構を通じて形成された連星は、検出器のバンドに入るまでにかなりの軌道離心率を保持している可能性がある。このような系の天体物理学的起源を理解するためには、この離心率を検出・特徴付けることが極めて重要である。最先端の準円軌道整列スピン波形モデルは存在するが、高次モードやスピン効果を含む離心系に対する正確な時間領域フェノメノロジカル・モデルは限られている。既存の離心系モデルは、有効一体系（EOB）形式や数値相対論（NR）サロゲートに依存することが多いが、これらは計算コストが高かったり、パラメータ空間のカバー範囲が制限されていたりする（例：非スピンまたは質量比が近い連星に限定されている）。さらに、離心系連星に関する公開されているNRシミュレーションは乏しく、新しいモデルの開発や検証を困難にしている。

手法
著者らは、離心系かつ整列スピンを持つ連星ブラックホール（BBH）のための新しい時間領域フェノメノロジカル・モデルであるIMRPhenomTEHMを導入する。このモデルは、正確な準円軌道・整列スピンモデルであるIMRPhenomTHMに、離心系ポスト・ニュートン（PN）補正を加えることで構築されている。

• モデル構築: コアとなる戦略は、IMRPhenomTHMの各モードの振幅と位相に、離心系PN補正を加えることである。モデルには、支配的な $(2, \pm 2)$ 球面調和モードに加え、劣支配的なモード $(2, \pm 1)$, $(3, \pm 3)$, $(4, \pm 4)$, $(5, \pm 5)$ が含まれる。
• 離心系ダイナミクス: モデルは、非スピンおよびスピン補正（スピン・オービットおよびスピン・スピン）を含む、完全な3PN軌道平均ダイナミクスを組み込んでいる。軌道の進化は、エネルギーおよび角運動量フラックスを接続するバランス方程式から導出された、平均化された（secular）進化方程式によって駆動される。
• 座標系: 著者らは、修正調和（MH）PN座標と有効一体系（EOB）座標の両方で、この平均化された進化方程式を実装している。EOB座標がNRシミュレーションとの一致度が高く、波形ベースの離心率の定義（$e_{GW}$）ともより密接に対応していることを見出し、EOBをデフォルトとして採用している。
• 円軌道化の仮定: 本モデルは、連星が合体までに円軌道化することを仮定している。したがって、合体およびリングダウンの局面では離心系補正は消失し、基礎となるIMRPhenomTHMモデルへのスムーズな遷移を保証する。
• 位相の再定義: ゲージ依存性を扱うために、著者らは位相の再定義を検討した。しかし、NRで較正されたQCベースラインとの一貫性を維持するため、デフォルトのモデルでは特定の位相シフトを適用せず、基礎となるモデルのNR較正が、必要な調整をすでに考慮していると仮定している。
• 計算効率: モデルは、軌道進化のための結合微分方程式を解くための高度に最適化されたアプローチを利用している。QC周波数の進化を特定のグリッド上で補間し、平均化された成分と振動成分を分離することで、モデルはEOBベースの離心系モデルと比較して大幅な高速化を実現している。

主な貢献

1. 初の時間領域フェノメノロジカル離心系モデル: IMRPhenomTEHMは、高次多重極を含む、離心系整列スピンBBHのための最初の時間領域フェノメノロジカル波形モデルとして提示される。
2. 多重極およびスピンの包含: モデルは最大$l=5$までのモードを含み、最大3PNのスピン・オービットおよびスピン・スピン相互作用を含む（ただし、QCベースラインとの一貫性のためにデフォルトでは2PNを使用）。
3. NRによる検証: モデルは、Simulating eXtreme Spacetimes (SXS) カタログからの28個の公開離心系NRシミュレーションを用いて検証されている。直接的な離心系NRデータセットへの較正なしでも、ほとんどのケースで不忠実度（ミスマッチ）は2%を下回っている。
4. 堅牢性とパラメータ空間: モデルは、広いパラメータ空間（質量比 $q \in [1, 20]$、スピン $\chi_i \in [-0.995, 0.995]$、および10 Hzでの離心率 $e < 0.4$）にわたって堅牢であることが示されている。
5. 計算性能: IMRPhenomTEHMは、低〜中程度の全質量において、最先端の離心系EOBモデル（SEOBNRv5EHMおよびTEOBResumDali）よりも約1桁高速であることが示されており、パラメータ推定（PE）研究に非常に適している。

結果

• 準円軌道極限: モデルはQC極限を忠実に再現しており、パラメータ空間全体にわたるIMRPhenomTHMに対するミスマッチは$10^{-5}$を下回っており、離心系拡張によるPE研究へのバイアスがないことを保証している。
• NR比較: 28個のSXS NRシミュレーションとの比較において、IMRPhenomTEHMは一般に2%未満のミスマッチを達成している。高い初期離心率（$e_0 \sim 0.3$）では、離心系展開されたPN式の限界に起因してミスマッチがわずかに増加する。
• NRインジェクションに対するパラメータ推定（PE）: 3つのSXSシミュレーションを用いたゼロノイズインジェクション研究において、モデルは注入されたパラメータ（質量、スピン、離心率、平均アノマリー）を90%信頼区間内で正常に回収している。高次モードの包含により、光度距離および傾斜角の回収率が向上している。
• 実イベント解析: モデルはGW150914およびGW190521の再解析に適用された。
  • GW150914: 結果は文献と一致しており、非離心、非スピン、質量が同程度の連星を支持している。ベイズ因子は、わずかにQC仮説を支持している。
  • GW190521: モデルはSEOBNRv5EHMを用いた従来の研究と一致する結果を出している。中央値の離心率はゼロではないものの、事後分布は主に情報量を欠いており、これはこの高質量イベントにおける離心率の証拠の欠如と一致している。本研究は、IMRPhenopTEHMが高質量における波形生成の制約による、他のモデルで見られるような人工的な離心率の事前分布カットを回避できることを強調している。
• タイミング・ベンチマーク: $M=10 M_\odot$ かつ $e=0.2$ の連星に対して、IMRPhenomTEHMは（高次モードを含めて）約1秒で波形を生成する。これに対し、SEOBNRv5EHMは約8秒、TEOBResumDaliは約$10^4$秒を要する。

意義および主張
論文は、IMRPhenomTEHMが、整列スピンを持ち、10 Hzでの離心率が$e < 0.4$であるブラックホール連星の、信頼性が高く計算効率的な記述を提供すると主張している。その主な意義は、離心系連星の検出が増加することが予想される次世代の重力波観測ランに向けた、大規模なベイズ推論研究を促進できる点にある。高速で多重極を含む、低速なEOBモデルに代わる選択肢を提供することで、現在の最先端モデルにおける計算上の制約を受けることなく、離心系の系統的な探索を可能にする。著者らは、現在の精度は既存の検出器には十分であるが、今後の課題として、振幅における離心系補正の取り扱いの改善（$O(e^6)$ 展開を超えたもの）および、スピン歳差運動の組み込みの可能性を挙げている。