A Stellar Role Reversal: Multiple Features in the Mass and Mass Ratio Distributions of Merging Binary Black Holes from Stable Mass Transfer

技術要約：恒星の役割の逆転：安定質量移動に由来する合体連星ブラックホールにおける質量および質量比分布の多重な特徴

問題提起
重力波（GW）イベントのカタログが増加しているにもかかわらず、合体する連星ブラックホール（BBH）集団の起源は依然として不透明である。観測によってBBHの質量関数（例：$10 M_\odot$付近と$35 M_\odot$付近のピーク）や質量比分布に特徴があることが明らかになっているが、測定可能な特性の縮退や、恒星進化および連星相互作用の不確実性により、異なる形成チャネルの寄与を解明することは困難である。本研究では、すべてのロッシュ・ローブ・オーバーフロー（RLOF）フェーズが安定であり、共通外層（Common Envelope; CE）チャネルを除外した「安定質量移動（Stable Mass Transfer; SMT）のみ」のチャネルを具体的に調査する。著者らは、質量移動の降着効率（$\beta$）に関する仮定が、合体するBBHの主星質量および質量比分布にどのように明確な特徴を刻み込むかを決定することを目的としている。

手法
本研究は、バイナリ・ポピュレーション・シンセシス（BPS）と解析的フレームワークを組み合わせた二重のアプローチを採用している。

1. ポピュレーション・シンセシス: 著者らは、cosmicコード（v4.1.0）を用いて、合体するBBHの宇宙論的集団をシミュレートする。

   • 初期条件: シミュレーションは、30の金属量（$10^{-4}$から$0.03$）、クランパ（Kroupa 2001）の初期質量関数（$0.08$から$100 M_\odot$）、および連星率0.7をカバーするグリッドを用いる。
   • 物理学: モデルには、金属量依存の星風（Vink et al. の処方）、潮汐円化、およびドナーの進化状態（Case A, B, C, および BC）に基づく特定の質量移動安定性基準が含まれる。
   • コンパクト天体の形成: ブラックホール（BH）質量は、「遅延型」超新星処方（Fryer etal. 2012）を用いて割り当てられ、最小BH質量は$3 M_\odot$とする。対不安定超新星（PISN）は、最大BH質量を$44 M_\odot$に制限するようにモデル化される。誕生時のキックは対数正規分布から抽出され、質量降着（fallback）によって抑制される。
   • 質量移動: 基準モデルでは、恒星降着体に対しては保存的質量移動（$\beta=1$）、コンパクト天体降着体に対してはエディントン限界降着（$\beta \approx 0$）を仮定する。
   • 宇宙論: cosmergeを介してIllustrisTNGシミュレーションに基づく星形成史（SFH）を適用し、$10^6$個の合体イベントのカタログを生成し、ハッブル時間内に合体するSMTのみのシステムにフィルタリングする。

2. 解析的フレームワーク: 著者らは、降着効率（$\beta_1$）、臨界質量比（$q_{crit}$）、および初期質量比（$q_{ZAMS}$）を、結果としての最小BH質量に結びつけるために、解析的モデル（van Son et al. 2022）を適応させている。修正点には、可変的なコア質量比（$f_{core}$）および、cosmicシミュレーションにより良く一致させるための特定の超新星質量放出関数が含まれる。

主要な貢献

• サブポピュレーションの特定: 本論文は、SMTのみのチャネルが、自然に二つの観測的に区別可能なサブポピュレーション、すなわち質量比逆転（Mass Ratio Reversed; MRR）集団と非MRR集団を生み出すことを示している。
• 分離のメカニズム: 著者らは、これらのサブポピュレーションの境界を、異なる進化段階（Case A, B, C）における質量移動の安定性基準に直接マッピングしている。異なるドナー型の安定性基準が、主星質量–質量比平面において相補的な領域を切り出すことを示している。
• パラメータ感度: 本研究は、降着効率（$\beta_1$）が合体率および質量分布に与える影響を系統的に探索し、この不確実なパラメータの変化が結果としてのBBH集団をどのように変化させるかを理解するための枠組みを提供している。

結果

1. 質量分布における明確な特徴:

   • MRRサブポピュレーション: 高い主星質量（$20-40 M_\odot$）と高い質量比（$q_{BBH} > 0.75$、0.9付近にピーク）によって特徴付けられる。これらのシステムは、主に第一回質量移動（MT1）の間にCase A質量移動を行う。
   • 非MRRサブポピュレーション: 低い主星質量（$10 M_\odot$付近にピーク）と、より広い質量比分布（$q_{BBH} \sim 0.3 - 0.7$）によって特徴付けられる。これらのシステムは、主にMT1におけるCase BまたはC質量移動を介して形成される。
   • 分離: 基準モデルでは、$\sim 20 M_\odot$における主星質量のカットオフが、これら二つのサブポピュレーションを効果的に分離する。

2. 降着効率（$\beta_1$）の役割:

   • $\beta_1$を増加させると（より保存的な質量移動）、安定進化のためのパラメータ空間が縮小するため、総合体率が減少する。
   • $\beta_1$を高めることは、MRRシステムの割合を増加させ、両方のサブポピュレーションの主星質量分布をより高質量側へシフトさせる。
   • $\beta_1$が増加するにつれて、MRRサブポピュレーションはより顕著になり、より高い質量まで広がるが、非MRRのピークは$10 M_\odot$付近に留まりつつも狭くなる。

3. 境界としての安定性基準:

   • MRR集団と非MRR集団の間の遷移は、初期質量比（$q_{ZAMS}$）と特定の質量移動ケースに対する安定性限界（$q_{crit}$）の相互作用によって支配される。
   • Case A質量移動は最も厳格な安定性基準を持ち、安定性を維持するために高い$q_{Z_AMS}$と高い初期質量を必要とし、これが高質量MRRピークを形成する。
   • Case C質量移動は、より低い$q_{ZAMS}$とより低い質量を許容し、低質量非MRRピークに寄与する。

意義と主張
著者らは、安定質量移動の物理学が、区別可能で明確な特徴を重力波集団の中に作り出すと主張している。具体的には、SMTのみのチャネルは、主星質量の二峰性と質量比の二分性を予測しており、これは現在の、および将来の重力波検出器の能力と一致している。本研究は、観測されたGWTC-5カタログの特徴（例：低質量ピークおよび高質量・高比率システム）が、これら二つの異なるSMTサブポピュレーションのシグネチャーである可能性を示唆している。

本論文は、現在のモデルでは、縮退のために観測されたイベントを特定のチャネルに決定的に割り当てることはできないと控えめに結論付けているが、特定された特徴（$\sim 20 M_\odot$での分裂および高$q$のMRR集団）はテスト可能な予測として機能する。本研究は、恒星半径の処方、星風による質量損失、および固定された安定性基準の簡略化に関する制限事項を認め、これらの要因が予測される特徴の正確な位置をシフトさせる可能性があるとしている。