Assessing the imprint of eccentricity in GW signatures using two independent waveform models

技術的概要：重力波信号における軌道離心率の痕跡の評価

問題提起
合体するコンパクト連星からの重力波（GW）信号には、その軌道特性と固有特性に関する情報が含まれている。軌道離心率は、形成チャネルを区別する重要な指標であり、特に、孤立した連星進化（準円軌道であることが期待される）と、高密度の恒星環境における動的な組み立て（測定可能な離心率を保持する可能性がある）を区別するものであるが、この痕跡を検出することは困難である。過去の分析では、潜在的な離心率を持つ候補事象が特定されてきたが、その結果はしばしば分析設定、波形モデル、データ条件付け（例えば、グリッチの軽減）に対して敏感であった。さらに、異なる波形モデルを使用するとパラメータ推論に大幅な差異が生じる可能性があり、ソースパラメータの信頼性のある識別を複雑にしている。本研究は、軌道離心率と整列したスピンを考慮した最先端の波形モデルを用いて、多数の GW 事象群に対して離心率を包括的かつ一貫して評価する必要性に対応するものである。

手法
著者は、LIGO-Virgo-KAGRA（LVK）ネットワークの O3 および O4a 観測ラン中に特定された 162 の確実な GW ソースに対してパラメータ推論を行った。サンプルには、偽警報率（FAR）が年間 1 回未満のすべての事象に加え、暫定的な離心率の兆候を示す 11 の追加 O3 事象が含まれる。

• 推論エンジン: 本研究では、RIFT（Rapid Inference Framework for Transients）パラメータ推論エンジンを使用し、複雑な波形の評価に伴う高い計算コストを軽減するために、周辺尤度を計算し事後分布を効率的に構築する 2 段階の反復プロセスを採用した。
• 波形モデル: 離心率を評価するために、2 つの独立して実装された最先端の波形モデル、すなわちSEOBNRv5EHMとTEOBResumS-Daliが採用された。両モデルとも、利用可能な高次モードを用いて、軌道離心率と整列したコンパクト天体のスピンを考慮している。
• 分析設定:
  • スピンは軌道角運動量と整列しているものとして制約された（"z 事前分布"）。
  • サンプリング周波数は、RIFT の内部離散時間サンプリングに対応するため、少なくとも 4096 Hz（一部事象では 8192 Hz）に設定された。
  • 離心率（$e$）と平均近点角（$\ell$）は、一様事前分布を持つパラメータとして含まれた。GWTC-3 事象については、事前分布範囲を $e \in [0, 0.9]$ とし（最終事後分布ではモデルの有効性を反映するため $e \in [0, 0.5]$ に切り捨て）、GWTC-4 事象については事前分布を $e \in [0, 0.5]$ とした。
  • 開始周波数は、事象の特性と過去の分析の選択に基づき変動した（GWTC-4 は 13.33 Hz、より古い事象は 20 Hz または 10 Hz）。
  • 雑音パワースペクトル密度（PSD）の推定値は、各事象に対して Bayeswave を用いて独立して生成された。
• 比較指標: 本研究では、2 つの波形モデル間の結果、および以前に発表された準円軌道（QC）の結果との比較を行った。一貫性は、周辺事後分布間のジェンセン・シャノン（JS）ダイバージェンスを用いて定量化された。ベイズ因子は、サベージ・ディッキー密度比を用いて計算され、離心率仮説に対する準円軌道仮説に対する支持度を定量化した。

主な貢献

1. 包括的な集団分析: 本研究は、単一の分析フレームワーク（RIFT）内で、2 つの異なる高精度な離心率波形モデルを用いて、162 事象からなる大規模で一貫したサンプルに対する離心率のパラメータ推論を初めて体系的に提供した。
2. モデル一貫性の評価: 2 つの独立した波形実装を採用することで、波形の系統誤差に対する離心率推論の堅牢性を厳密に検証した。
3. 系統誤差の特定: 離心率および他の固有パラメータに関する結論が波形モデル間で著しく異なる特定の事象を特定し、波形の系統誤差が顕著なパラメータ空間の領域（特に高質量系）を浮き彫りにした。
4. 候補事象の再評価: 以前に特定された離心率候補（例：GW200129、GW190521）および新たな候補を再評価し、更新されたベイズ因子と事後分布を提供した。

結果

• 一般的な一貫性: 162 の事象の大部分において、SEOBNRv5EHM と TEOBResumS-Dali の両モデルは、 chirp 質量、質量比、有効スピン、および離心率の事後分布が強く一致しており、非常に一貫した結論を導き出した。ほとんどの事象は、離心率がゼロであることと一致している。
• 離心率候補:
  • GW200129 065458: 両波形モデルとも離心率に対する肯定的な証拠を見出した。SEOBNRv5EHM は $\log_{10} B \approx 3.12$、$e \approx 0.12$ を与え、TEOBResumS-Dali は $\log_{10} B \approx 3.59$、$e \approx 0.11$ を与えた。ただし、著者はこの事象の結論が分析設定（サンプリング周波数、基準周波数、グリッチ軽減など）に敏感であることを指摘しており、これは過去の文献と一致する。
  • GW231001 140220: TEOBResumS-Dali は離心率に modest なピーク（$e \approx 0.36$、$\log_{10} B \approx 0.30$）を示したが、SEOBNRv5EHM は有意な証拠を見出さなかった（$\log_{10} B \approx -0.12$）。著者は、TEOB の特徴の一部を、高質量系における波形の系統誤差に起因すると帰属している。
  • GW231123 135430: TEOBResumS-Dali は離心率に対して強い支持を示した（$\log_{10} B \approx 6.05$、$e \approx 0.46$）が、分布は事前分布の上限に対して「レールに乗り」ているように見えた。これに対し、SEOBNRv5EHM は有意な証拠を見出さなかった（$\log_{10} B \approx 0.04$）。著者は、この事象は著しい波形の系統誤差を示しており、その不一致に基づいて確認された離心率候補とはみなさないと結論づけた。
• 系統誤差: いくつかの高質量事象（GW230814、GW231001、GW231123）において、モデル間で離心率および有効スピンパラメータに関して著しい不一致（JS ダイバージェンス > 0.2）が観測された。TEOBResumS-Dali は、高質量系において $e \gtrsim 0.5$ の高い離心率に無視できない確率を割り当てる傾向があり、これは著者によって波形の系統誤差である可能性が示唆されている。
• 他の研究との比較: 本研究の結論は、GW190521 074359 や GW190701 203306 などの事象においてより強い離心率の証拠を見出した以前の報告としばしば異なっていた。著者は、これらの差異を分析設定、波形モデル、データ条件付け（例えば、デグリッチング戦略）の相違に起因すると帰属している。

意義
本論文は、一貫したフレームワークを用いて、現在の GW 事象カタログに対する離心率パラメータ推定の包括的なカタログを確立した。主な意義は、2 つの独立した最先端の離心率波形モデルが、ほとんどの連星ブラックホール合体の準円軌道的性質について概ね一致していることを示しつつ、特定の高質量ケースにおいて著しい系統誤差が生じ得ることを実証した点にある。著者は、よく研究され分析設定に敏感な GW200129 065458 を除き、サンプル内の他の連星ブラックホール候補は離心率に対する有意な証拠を示さないと結論づけた。この研究は、潜在的な離心率の痕跡を解釈する際に、分析設定と波形の系統誤差を評価することの決定的な重要性を強調し、モデルや分析選択のアーティファクトに過ぎない可能性のある特徴を過剰に解釈しないよう警告している。その結果は、コンパクト連星の形成チャネルの将来の集団レベルの研究のための堅牢な基準を提供する。