First eccentric inspiral-merger-ringdown analysis of neutron star-black hole mergers

技術要約：中性子星・ブラックホール合体における初の離心率を持つインスパイラル・マージャー・リングダウン解析

問題と背景
コンパクト連星の重力波（GW）検出は、天体物理学的な形成チャネルの研究を可能にしてきた。孤立進化を通じて形成される恒星質量連星は、検出器のバンドに入る前に円軌道化すると予想されるが、高密度環境（例：星団）や三重系における力学的相互作用は、軌道の離心率を保持させることがある。この離心率を検出することは、形成チャネルを区別する上で極めて重要である。中性子星・ブラックホール（NSBH）イベントGW200105を用いた、先行するインスパイラルのみの離心率を持つ歳差運動波形モデル（Morras et al. 2025a）による解析は、軌道の離心率に関する強い証拠を提供した。しかし、サブドミナントな高次モードを含む、インスパイラル、マージャー、およびリングダウン（IMR）の全フェーズを組み込んだ完全な解析は、NSBHシステムに対してはまだ行われていなかった。さらに、既存のIMR離心率モデルは、離心率セクターにおける数値相対論（NR）へのキャリブレーションを欠いているか、あるいは低質量系のパラメータ推定（PE）において計算コストが高すぎる場合がある。

手法
著者らは、3つのNSBHイベント（GW200105、GW200115、GW230529）に対して、初のフルIMRパラメータ推定解析を行う。この解析では、準円軌道（QC）モデルであるIMRPhenomTHMを拡張した、整列スピンを持つ離心率モデルであるIMRPhenomTEHM波形モデルを使用している。本手法の主な特徴は以下の通りである：

• 波形モデル: IMRPhenomTEHMは、支配的な$(2, \pm 2)$モードと、サブドミナントな高次モード $(l, m) = \{(2, \pm 1), (3, \pm 3), (4, \pm 4), (5, \pm 5)\}$ を含む。これは質量比1/18までのQC NRシミュレーションにキャリブレーションされており、連星は合体時に円軌道化すると仮定している。本研究で対象とする質量比および信号対雑音比（SNR）においては、物質効果（潮汐破壊）は無視できると判断し、含まれていない。
• 推論フレームワーク: ベイズ推論は、dynestyネストサンプリングアルゴリズムを用いたbilbyフレームワークを用いて行われる。解析では、成分質量、整列スピン、基準離心率（$e_{20\text{Hz}}$）、平均アノマリー、および外部パラメータを含む13次元のパラメータ空間をサンプリングする。
• 事前分布: 2種類の離心率事前分布がテストされた：一様分布（$e \in [0, 0.5]$）と、対数一様分布（$e \in [10^{-4}, 0.5]$）である。
• 系統誤差の制御: 時間領域の離心率モデルに特有の計算上の課題、特に、尤度積分開始周波数（20 Hz）におけるすべての関連する高次モードを捉えるために、波形生成を十分に早期に開始する必要性について対処している。著者らは、開始周波数を17 Hz、20 Hz、13.3 Hzとテストし、低周波分解能とノイズ混入を評価するために信号の長さ（32秒 vs 64秒）を変化させている。

主な結果

1. GW200105（離心率の証拠）: 本解析は、非ゼロの軌道離心率の強い証拠を裏付けている。

   • 一様事前分布を用いた場合、事後分布のピークは $e_{20\text{Hz}} = 0.12^{+0.02}_{-0.03}$ にあり、QC仮説に対する離心率仮説のベイズ因子は $\log_{10} B_{E/QC} = 1.22^{+0.12}_{-0.12}$ であった。
   • 対数一様事前分布を用いた場合、主要なモードはゼロ方向にシフトするが、サブドミナントなモードは $e_{20\text{Hz}} \sim 0.11$ に残り、ベイズ因子は $\log_{10} B_{E/QC} = 0.11^{+0.11}_{-0.11}$ と緩やかな値となった。
   • 離心率の事後分布は二峰性を示し、$e \sim 0.11$ と $e \sim 0.14$ にピークを持つ。これは、離心率、成分質量、および整列スピンの間の縮退に起因しており、これらが信号の持続時間に影響を与えるためである。高い離心率を持つモードは、より低い一次質量（$m_1 \sim 8 M_\odot$）およびより反並行なスピンと相関している。
   • 開始周波数と信号の長さを変化させたテストにより、二峰性と非ゼロの離心率信号の頑健性が確認されたが、著者らは、低周波分解能の限界（32秒のデータウィンドウによるもの）が事後分布の構造に影響を与えている可能性についても言及している。

2. GW200115 および GW230529（QCとの整合性）:

   • 両イベントは、準円軌道連星と矛盾しない結果となった。
   • GW200115は $e_{20\text{Hz}} = 0.02^{+0.03}_{-0.02}$、$\log_{10} B_{E/QC} = -0.91$ を示した。
   • GW230529は $e_{20\text{Hz}} = 0.01^{+0.02}_{-0.01}$、$\log_{10} B_{E/QC} = -1.36$ を示した。
   • これらの結果は、これらの特定のイベントにおいて離心率の証拠がないことを示している。

3. パラメータの縮退: モデルに離心率を含めることで、QCのみの解析と比較して、回収されるQCパラメータ（質量およびスピン）にシフトが生じる。GW200105については、離心率によって引き起こされる信号の短縮を補償するために、離心率モデルはより低い質量を支持する傾向がある。

意義と主張
本論文は、NSBH合体における初のフルIMR整列スピン離心率パラメータ推定解析を提示していると主張している。その主な意義は以下の通りである：

• 検証: 高次モードを含む完全なIMRモデルを用いて、GW200105の離心率の主張を検証し、インスパイラルのみの解析（Morras et al. 2025a）の結果を補強したこと。
• 手法の進展: 離心率セクターにおける直接的なNRキャリブレーションを欠いているにもかかわらず、計算効率の高いIMR離心率モデルが、低質量系（NSBH）の離心率測定に適用可能であることを示したこと。
• 集団に関する知見: 確認された少なくとも一つのNSBH合体（GW200105）が、孤立連星進化ではなく、力学的チャネルを通じて形成された可能性が高いという証拠を提供した一方で、他のイベント（GW200115, GW230529）は標準的な形成チャネルと矛盾しないことを示したこと。
• 系統的課題への意識: 時間領域の離心率PEにおける計算上および系統的な課題、特に低周波分解能、信号の長さ、および離心率とスピン/質量の縮退の相互作用を浮き彫りにした。著者らは、GW200105およびGW200129（別の離心率候補）におけるデータ品質の問題（ノイズのグリッチ）について、将来の解析において慎重な評価が必要であると述べている。

著者らは、系統的な不確実性（二峰性や低周波分解能の限界など）は依然として存在するものの、本結果は重力波天文学におけるロバストな離心率測定への一歩であると結論付けている。