Twin Peaks: Resolving Features in the Binary Black Hole Mass Function with COSMIC-METISSE

問題提起
LVKコラボレーションによる重力波（GW）観測は、ブラックホール連星（BBH）の集団における特定の特性、特に主質量スペクトルにおける $10 M_\odot$ 付近のグローバルなピークと、$q \approx 0.7$ における質量比分布の特徴を明らかにしている。孤立連星進化に基づく個体群合成モデルはこれらの集団を再現できるが、その予測能力は基礎となる恒星物理学の不確理性に制限されている。具体的には、標準的な高速個体群合成コードは、SSE（Single Star Evolution）などのフィッティング公式に依存することが多く、これらは広範な金属量にわたる様々な恒星風や対流境界混合（CBM）の影響を正確に捉えられない可能性がある。これらの不確実性は、質量移動の安定性、共通外層（CE）形成、および最終的な残骸質量といった連星進化の結果へと大きく波及する。

手法
著者らは、金属量 $10^{-3} \le Z/Z_\odot \le 7$ および質量最大 $300 M_\odot$ にわたる、オープンソースコード MESA（リリース r24.08.1）を用いて計算された新しい恒星トラックのグリッドを提示している。このグリッドには、水素に富む星のための85個の質量点と、裸のヘリウム星のための42個の質量点が含まれる。著者らは2つの主要な恒星物理パラメータを変化させている：

1. 風による質量放出： 標準的な処方（Vink et al. 2001; Jager & v. d. Hucht 1988）と、より新しい観測に基づいた適合式（Krticka et al. 2025; Decin et al. 2023）の比較。
2. 対流境界混合（CBM）： 標準的な指数関数的オーバーシューティング（Herwig 2000）と、対流浸透モデル（Johnston et al. 2024）の比較。

これらのMESAトラックは METISSE（Method of Interpolation for Single Stellar Evolution）を用いて補間され、高速個体群合成コード COSMIC に投入される。このパイプライン（COSMIC-METISSE）により、入力となる恒星物理学を系統的に変化させながら、局所宇宙（$z < 0.2$）における合成BBH合体集団を生成することが可能となる。本研究では、コア・エンベロープの境界を持たないドナーは合体すると仮定する「悲観的」なCEシナリオを採用しているが、比較のために「楽観的」なシナリオもテストされている。

主な結果

• $10 M_\odot$ における「ツインピーク」： ほとんどのモデル変異（Standard, New Winds, New CBM）において、主質量スペクトルは $10 M_\odot$ のピークの両側に位置する2つの明確なサブ集団を示している：一つは $\approx 8 M_\odot$、もう一つは $\approx 13 M_\odot$ である。
  • $\approx 8 M_\odot$ のピークは、質量比逆転（MRR）を起こさない系によって支配されている。
  • $\approx 13 M_\odot$ のピークは、MRRを起こす系によって支配されている。これらのMRR系は、主に初期質量比 $q_{ZAMS} \gtrsim 0.5$ の親星から形成され、最終的な質量比が $q \approx 0.7$ 付近に集まる。
• 恒星物理への感度： 弱い風と対流浸透を組み合わせた「New Winds & CBM」モデルでは、これら2つのサブ集団が $9 M_\odot$ 付近の単一のピークへと統合される。これは、特定の物理の組み合わせによって高金属量の連星がCEフェーズを生き残り、MRRを起こすことが可能になり、MRR集団がより低質量側へシフトして非MRR集団と重なるために起こる。
• 合体率： 恒星物理の選択は、予測される局所合体率に劇的な影響を与える。対流浸透を用いるモデル（New CBM および New Winds & CBM）は、$\approx 120\text{--}168 \text{ Gpc}^{-3}\text{yr}^{-1}$ の合体率を予測しており、これは標準的なオーバーシューティングを用いるモデル（$\approx 25 \text{ Gpc}^{-3}\text{yr}^{-1}$）よりも約4〜6倍高い。この増加はCEチャネルの変調によって駆動されており、対流浸透が2度目のロッシュローブオーバーフローを遅らせることで、より多くの系がCEフェーズを生き残り、ハッブル時間内に合体することを可能にしている。
• 金属量依存性： StandardおよびNew Windsモデルは、主に極めて低い金属量（$Z \lesssim 0.01 Z_\odot$）において合体するBBHを生成する。対照的に、New Winds & CBMモデルは、太陽組成および超太陽組成の金属量においても有意な合体サブ集団を生成し、これが $10 M_\odot$ 付近の単一ピークに寄与している。
• 質量比の相関： モデルは $q \approx 0.7$ における観測された特徴を再現しており、これはMRRイベントと強く関連している。MRR系においては、主質量と質量比の間に反相関が見られ、より高い主質量を持つ星ほど優先的にMRRを起こし、$q \approx 0.7$ に落ち着く。

意義と主張
本論文は、BBH質量関数に見られる「ツインピーク」は、必ずしも別々の天体集団ではなく、相対的な強度と位置が恒星風およびCBMの処方箋に高度に敏感な、2つのサブ集団（MRRおよび非MRR）の解釈であると主張している。著者らは以下のことを実証した：

1. 補間は極めて重要である： METISSEを介してMESAトラックを補間することで、従来のSSEフィッティング公式と比較して、恒星物理学の不確実性がどのように連星集団へと波及するかをより直接的に研究できる。
2. 恒星物理が率を駆動する： CBMおよび風の取り扱いの変化は、主にCEチャネルの効率を調節することによって、予測されるBBH合体率を約6倍変化させ得る。
3. 観測との整合性： 単一のモデルがLVKの局所合体率と赤方偏移進化の両方を完全に再現することはできないが、「New Winds & CBM」モデル（観測された率で正規化した場合）は、主質量スペクトルおよび質量比分布において、$10 M_\odot$ のピークや $q \approx 0.7$ の特徴に関して、GWTC-5 PixelPopの事後分布と顕著に良く一致する。

著者らは、BBH質量関数の特定の特徴を解明するには、大規模恒星進化、特に恒星風による質量放出と対流境界混合に関する精密な制約が必要であり、今後の研究は宇宙の星形成史およびCE形成メカニズムの不確実性に対処しなければならないと結論づけている。