Eccentric or circular? A reanalysis of binary black hole gravitational wave events for orbital eccentricity signatures

技術要約：離心軌道か、円軌道か？ 軌道離心率のシグネチャに関する連星ブラックホール重力波イベントの再解析

問題提起
連星ブラックホール（BBH）の形成チャネルは、重力波天文学における中心的な未解決問題である。孤立連星進化は、地上検出器の感度帯域に進入するまでに準円軌道（quasi-circular; QC）を形成することが予想されるが、動的な形成シナリオ（例：球状星団や銀河核におけるもの）は、有意な軌道離心率を持つ合体を生じさせることがある。離心率の検出は、動的な形成を示す「決定的な証拠（smoking gun）」となる。しかし、これまでの重力波イベントにおける離心率シグネチャの特定に関する試みは、計算上の制約により、リウェイティング手法や機械学習への依存、あるいは非離心パラメータへの解析の限定といった制約を受けてきた。さらに、パラメータ推定（PE）において離心率を無視することは、チャープ質量や有効スピンなどの推定されるソース特性にバイアスをもたらす可能性がある。

手法
本研究は、Advanced LIGOおよびAdvanced Virgo検出器の第1次から第3次観測ラン（O1–O3）で検出された17件のBBHイベントの再解析を提示する。解析には、楕円軌道にある整列スピンBBHのための新しいフェノメノロジー的時系列マルチポーラ波形モデルであるIMRPhenomTEHMモデルを使用する。このモデルは、IMRPhenomTHM（QC整列スピンモデル）を拡張し、2つの離心パラメータ、すなわち離心率（$e$）と平均近点角（$l$）を組み込み、完全な3PN平均軌道力学および$O(e^6)$までの離心補正を含んでいる。

主要な手法構成要素は以下の通りである：

• モデル比較： 著者らは、IMRPhenomTEHMモデルを、基礎となるQCモデル（IMRPhenomTHM）および、特定のイベントについては歳差運動モデル（IMRPhenhemTPHMおよびNRSur7dq4）と比較することで、離心率とスピン歳差運動を区別している。
• 事前分布（Priors）： 参照周波数10 Hzにおける離心率の2つの異なる事前分布、すなわち一様事前分布（$e \in [0, 0.5]$）と対数一様事前分布（$e \in [10^{-4}, 0.5]$）をテストしている。
• ベイズ推論： bilbyフレームワークとdynestyネステッドサンプラーを用い、後験分布および、円軌道仮説に対する離心率の証拠を定量化するためのベイズ因子（$B_{E/QC}$）を算出している。
• 系統誤差およびグリッチ緩和： グリッチが含まれることが知られているイベントGW200129について、著者らは、生のGWOSCデータ、gw_subtractで処理されたデータ、および（3つの異なるグリッチ・ドローを用いた）BayesWaveで処理されたデータを用いて解析を行い、データ処理が離心率の回収に与える影響を評価している。
• 高質量イベント： マージングリッチ（merger-ringdown）のみが観測される高質量イベントについては、サンプラー設定（ライブポイントを$n_{live}=2000$に増加）と波形の系統誤差の影響を調査している。

主な結果
再解析の結果、離心率の証拠を示唆する4つのイベントが特定されたが、その証拠の強さは様々である：

1. GW200129： このイベントは、最も強い離心率の証拠を示している。

   • 非デグリッチ（non-deglitched）のGWOSCデータを用いた場合、$\log_{10}$ ベイズ因子は一様事前分布で $\log_{10} B_{E/QC} \approx 5.14$、対数一様事前分布で $4.19$ となり、離心率仮説を著しく支持している。
   • gw_subtractデータを用いた解析も、強い支持を示している（$\log_{10} B_{E/QC} \approx 4.00$ および $3.35$）。
   • BayesWaveによる緩和（特に、動的な特徴への支持が最も低いBayesAドロー）を行った場合でも、離心率仮説は整列スピンおよび歳差運動を含むQCシナリオの両方に対して優位である（$\log_{10} B_{E/QC} \approx 1.30$）。
   • 回復された離心率は事前分布間で一致している（gw_subtractデータにおいて $e_{GW}^{10Hz} \approx 0.26$）。
   • 離心モデルの後験分布は、整列スピンのQCモデルよりも歳差運動サロゲートであるNRSur7dq4とよく一致しており、これは単純なスピン歳差運動よりも動的な効果を好むことを示唆している。

2. GW200208_22： このイベントは、事前分布の選択に大きく依存して、離心率への中程度の支持を示している。

   • 一様事前分布では、$\log_{10}$ ベイズ因子は約 $1.14$ であり、離心率仮説を支持している。
   • 対数一様事前分布では、支持は減少するが（$\log_{10} B_{E/QC} \approx 0.49$）、離心率仮説はQC歳差スピンモデルよりも依然として好ましい。
   • 離心率の導入はバイナリパラメータの制約を大幅に強め、推定される全質量を（QCモデルの $\sim 174 M_\odot$ から）$\sim 65 M_\odot$ へと減少させた。

3. GW190701 および GW190929： これらの高質量イベントは、離心率に対して軽微または決定力のない支持しか示していない。

   • 一様事前分布は離心率への選好を示しているが（例：GW190929では $\log_{10} B_{E/QC} \approx 0.56$）、対数一様事前分布はこの証拠を抑制する。
   • 著者らは、これらの高質量イベントについては信号の持続時間が短く、現在の離心モデルはマージ時における円軌道化を仮定していることを指摘している。したがって、マージ付近の波形の系統誤差と離心率を区別することは困難であり、強い主張は支持されない。

意義および主張
本論文は、IMRPhenomTEHMの計算効率が高いことにより、これまで計算コストによって制限されていた大規模なデータセットへの離心波形モデルのルーチンな適用が可能になったと主張している。本研究は以下の重要な知見を強調している：

• バイアスの緩和： 離心率を無視することは、特に全質量の過小評価や有効スピンの誤認など、パラメータ推定における系統的なバイアスを引き起こす可能性がある。
• 動的な形成： GW200129（および程度は低いがGW200208_22）における離心率の堅牢な検出は、高密度な恒星環境におけるような動的な形成チャネルの証拠を提供する。
• 縮退： 軌道離心率とスピン歳差運動の間には顕著な縮退が存在する。著者らは、両者を区別するには一般的な歳差運動を含む離心モデルが必要であるが、それらは現在のフレームワークではまだ完全に実装されていないと述べている。
• 今後の展望： 著者らは、検出器の感度が向上し、BBH集団の多様性が増すにつれ、離心率を含む（そして最終的には歳差運動を含む離心）モデルを組み込むことが、正確な集団研究や形成チャネルの特定に不可欠になると強調している。彼らは、自らの結果がGW200129における軌道離心率の存在を支持し、動的な起源という解釈を強化するものであると明言している一方で、現在のモデルはマージ時の円軌道化を仮定しているため、高質量で持続時間の短い信号に対して強い主張を行うことには限界があることも認めている。