Lower Your Rates: On Claims of a Binary Black Hole Merger-Rate Crisis

技術要約：「レートを下げよ：連星ブラックホール合体率の危機に関する主張に対して」

問題提起
近年の文献では、孤立連星進化のシミュレーションは、金属量依存的な宇宙星形成史の変化を考慮したとしても、観測される局所的な連星ブラックホール（BBH）合体率を一般的に過大評価していると論じられている。特定の研究（例：Sgalletta et al. 2025; Boco et al. 2026）は、この不一致を劇的な連星恒星進化モデルの修正を必要とする証拠として解釈しており、理論を観測（具体的には $\lesssim 20$ Gpc$^{-3}$ yr$^{-1}$ のレート）と整合させるには、星形成史単独の変化を超えた大幅な変更が必要であると示唆している。しかし、シミュレートされたBBHレートに関する文献は数桁にわたる幅を持っており、限定的なモデルセットから導き出された結論が、孤立連星進化研究のより広い領域全体に一般化できるかどうかは依然として不明である。

手法
著者らは、Mandel & Broekgaarden (2022) によって提示された、孤立連星進化研究による局所的（$z \approx 0$）な合体レートのコンピレーションに対する包括的な更新および拡張を提示する。

• データ・コンピレーション： 本研究は、2026年までに発表された57の研究から得られた1490個のシミュレートされたBBH合体レートを集計している。このコンピレーションには、金属量依存的な宇宙星形成史を考慮した固有の局所合体レートを報告している研究のみが含まれている。
• 比較分析： シミュレートされたレートを、LIGO–Virgo–KAGRA（LVK）コラボレーションによる制約、具体的にはGWTC-5の区間（約27.5–49.4 Gpc$^{-3}$ yr$^{-1}$）およびそれ以前のGWTCリリースと比較する。
• パラメータの分離： レートの変動の要因を特定するために、著者らは2543組のペアとなる「関係性」をサブモデル間で構築した。これらのペアは、単一の物理的仮定またはパラメータ（例：共通外層の効率、誕生時キック、質量移動の安定性）によって異なる。このアプローチにより、連星恒星進化物理学の特定の変化が合体レートにどのように影響するかを定量化することが可能となる。
• 分類： サブモデルは、パラメータファミリー（例：誕生時キック、共通外層進化、星形成史）およびポピュレーション・シンセシス・フレームワーク（例：COMPAS, SEVN, MOBSE, StarTrack）ごとにグループ化され、クラスタリングおよび「シミュレーション・サイロ」を分析するために使用される。

主要な貢献と結果

1. 普遍的な合体率の危機の否定：
   集計された1490個のBBHサブモデルのうち、約80%がGWTC-5の区間の上限を超えるレートを示しているが、観測されたレートを再現、あるいは下回る相当数のサブセットが存在する。文献の分布は高レート側に偏っているが、これは理論的枠組みが必然的に高いレートを生み出すことを意味しない。不一致の存在とその大きさは、普遍的なものではなく、極めてモデル依存的であることが示されている。

2. 低レートへの複数の経路：
   低BBH合体レートは、強い誕生時キックや低金属量星形成の減少に特有に関連しているという最近の主張とは対照的に、本解析は、低レートが広範な物理的仮定から生じ得ることを実証している。GWTC-5の区間と一致、あるいはそれを下回るモデルは、以下の変数を変化させた際に発見される：

   • 共通外層の効率および生存基準。
   • 外層結合エネルギー（$\lambda$）。
   • 質量移動の安定性と効率。
   • 角運動量喪失の処方。
   • 残骸質量の処方。
   • 恒星風および進化トラック。
   • 初期条件（例：初期質量関数）。
     局所的なBBH合体レートは縮退した観測量である。観測との一致は、明確に異なり、かつ互いに補完し得る複数の異なる仮定の組み合わせを通じて達成され得る。

3. 「シミュレーション・サイロ」効果：
   本研究は、シミュレートされたレートが、使用されるポピュレーション・シンセシス・コードに基づいて、明確な「シミュレーション・サイロ」へとクラスター化することを特定している。SEVN、COSMIC、COMPASなどのフレームワークは、最も多くのサブモデルを提供し、主にパラメータ空間の高レート領域を占めている。したがって、特定のサブセット（例：SEVNシミュレーションに基づくもの）から導き出された結論は、孤立連星進化の一般的な特徴ではなく、その特定のフレームワークに固有の特性的なレート分布を反映している可能性がある。異なるコード間の分布の重なりは、異なるフレームワークによって同様のレートが生成され得ることを裏付けるとともに、単一のフレームワークであっても仮定を変化させることで、極めて異なるレートを生み出し得ることを示している。

4. 観測的ベンチマークへの感度：
   認識される緊張の深刻さは、選択された観測的ベンチマークに敏感である。GWTC-4の中心値（$\sim 20$ Gpc$^{-3}$ yr$^{-1}$）をベンチマークとして用いると、フルなGWTC-5の区間を用いる場合よりも、より大きな過剰生産因子が示される。観測的に推論されるレートが進化するにつれ、定量的な緊張のレベルは減少する。

意義と主張
本論文は、普遍的なBBH合体率の危機という主張は、強くモデル依存的であることを論じている。著者らは、孤立連星進化が重力波観測と矛盾していると結論付けているのではなく、むしろ、その緊張は、大規模で不完全に探索されたモデル空間の特定の領域に限定されていると主張している。

本研究の主要な意義は、孤立連星進化が根本的に欠陥があるという結論を出す前に、複数のポピュレーション・シンセシス・フレームワークにわたって、高次元のパラメータ空間の全容を探索する必要性を強調することにある。結果は、現在の不一致が孤立連星進化のパラダイム自体の失敗ではなく、現在文献を支配している特定の仮定やコードの実装を反映したものであることを示唆している。今後の解決には、合体レートだけでなく、BBH質量分布、スピン、および赤方偏移進化などの他の観測量に対しても、モデルの結果を比較することが求められる。