Evidence for mass-dependent spin subpopulations in GWTC-4

技術的サマリー：GWTC-4 における質量依存性のスピンサブ集団の証拠

問題提起
合体する連星ブラックホール（BBH）系の天体物理学的起源は依然として不確実であり、さまざまな形成チャネルは、観測される質量分布とスピン分布に明確な特徴を刻み込むことが期待される。以前の集団研究では質量分布と有効スピン（$\chi_{\text{eff}}$）に特徴が特定されてきたが、$\chi_{\text{eff}}$ などの縮約された観測量は、異なる物理的スピン構成を混合させることで、それらを曖昧にする可能性がある。さらに、成分スピンの大きさ（$\chi_1, \chi_2$）の分布を推論することは技術的に困難であり、特に複数の狭い集団を支持するモデルにおいて顕著である。標準的なモンテカルロ（MC）尤度推定量は、狭い分布や集中したサブ集団を扱う際に高い分散に悩まされ、データ中の微妙な構造を隠蔽する可能性がある。本研究は、これら限界に対処するため、第 4 リーゴ・バーゴ・カグラ重力波過渡現象カタログ（GWTC-4）におけるスピン大きさの質量依存性を調査し、主質量（$m_1$）に応じて相対割合が変化する複数のスピンサブ集団を許容するモデルを調査する。

手法
著者は、階層ベイズ推論を用いて、GWTC-4 から得られた 153 の BBH 候補信号（誤報率 $< 1 \text{ 年}^{-1}$）から集団モデルのハイパーパラメータ（$\Lambda$）を推定する。

1. 集団モデル: ベースラインの質量および赤方偏移モデルは、Abac ら（2025a）による標準的な強力にモデル化された BBH 集団に従う。核心的な革新は、スピン大きさモデルにあり、単位正方形 $[0, 1] \times [0, 1]$ 上で、2 つの相関した切断された 2 次元正規分布（サブ集団 $a$ と $b$）の混合を仮定する。
2. 質量依存混合: 以前の定数混合モデルとは異なり、本研究は質量依存混合割合 $\eta(m_1)$ を導入する。この割合は、$m_s$ と $m_e$ で定義される遷移領域を跨いで、低質量値 $\eta_s$ から高質量値 $\eta_e$ へ線形に移行する。モデルは、低質量においてサブ集団 $a$ が支配的となるよう制約されている（$\eta_s \geq 0.5$）。
3. スピン傾きモデル: 分析では 2 つのスピン傾き構成を考慮する。スピン傾きの切断正規分布のピークが自由パラメータである「自由」モデルと、ピークが $\cos \theta_i = 1$（整列スピン）に制約される「固定」モデルである。標準的な分析では自由傾きモデルを使用する。
4. 尤度推定（TGMM）: 狭い集団に対する標準的な MC 推定量に内在する分散問題を軽減するため、著者は切断ガウス混合モデル（TGMM）を利用する。彼らは TGMM を基礎となる事象の事後分布と回復された注入事象に適合させる。これにより、スピン領域での半解析的積分と質量/赤方偏移領域での MC 積分が可能となり、サブ集団が狭い場合やパラメータ空間の隅に集中している場合であっても、分散を大幅に低減する。
5. 推論: ハイパー事後分布は、numpyro に実装されたノ・Uターン・サンプラー（NUTS）ハミルトニアン MC アルゴリズムを用いてサンプリングされる。

主要な貢献

• 質量依存混合: 本論文は、低速回転と高速回転の BBH の割合が主質量の関数として変化するモデルを導入し、これを制約する。これは、一定の全球的混合を仮定するものではない。
• 正則化された尤度: 集団尤度を正則化するための TGMM の適用は、標準的な推定量が分散のために見逃す可能性のある狭いサブ集団を検出するための堅牢な手法を提供する。
• 成分スピン分解能: 分析は、縮約された観測量にのみ依存するのではなく、直接成分スピン空間（$\chi_1, \chi_2$）で実行され、形成チャネルに関する特定の構成（例：非対称スピン）の検出を可能にする。

結果

• 2 つのサブ集団の証拠: データは、単一集団に対する 2 集団モデルを支持する対数ベイズ因子 $\log_{10} B = 4^{+6}_{-2}$ を伴い、2 つの明確なスピン大きさサブ集団に対する強力な証拠を提供する。
• 質量依存遷移: 質量依存混合割合に対する強力な証拠がある。遷移は $m_s \approx 36^{+19}_{-21} M_\odot$ で始まり、$m_e \approx 68^{+65}_{-31} M_\odot$ で終わる。
  • 低質量領域（$m_1 \lesssim 35 M_\odot$）: 集団は低速回転するサブ集団（$\eta_s \approx 0.95$）によって支配されている。スピンはゼロ付近に密に集中しており、平均値は $\mu_a \approx (0.07, 0.09)$、幅は $\sigma_a \approx (0.08, 0.12)$ である。データは、同一に消滅するスピン（$\chi_1 = \chi_2 = 0$）という仮説を除外しない。
  • 高質量領域（$m_1 \gtrsim 70 M_\odot$）: 高速回転するサブ集団が有意になる（$\eta_e \approx 0.17$）。この成分はより広範で、$\mu_b \approx (0.75, 0.69)$ にピークを持ち、大きな幅（$\sigma_b \approx (0.43, 0.49)$）を持つ。これは両方のブラックホールが高いスピン（$\chi \sim 0.7$）を持つ構成を支持するが、非対称スピン（例：$\chi_1 \sim 0.7, \chi_2 \sim 0$）への無視できない尾部も保持する。
• 相関: 低速回転サブ集団におけるスピン間の弱い反相関の証拠がある一方、高速回転サブ集団における相関は事前分布と一致している。
• 堅牢性: 結果は、最も極端な事象（GW231123 および GW190517）を除外するリーブ・ワン・アウト分析に対して堅牢である。定数混合割合モデルは質量依存モデルに対して不利である（$\log_{10} B \approx 1.6$）が、サブ集団の推定特性はモデル間で一貫している。
• 傾きモデルの感度: 低速回転集団の特性の推論はスピン傾きモデルに敏感である。「自由」傾きモデルは消滅するスピンを許容するのに対し、「固定」傾きモデルは低速回転のピークを $\chi \approx 0.15$ へシフトさせる。

意義と主張
著者は、その結果が BBH 形成モデルに対する新たな実証的な目標を提供すると主張する。発見は、BBH 質量スペクトル全体にわたって異なる形成チャネルが機能していることを示唆する。

1. 低質量: 効率的な角運動量輸送がブラックホールに無視できる出生スピンを残す、孤立連星進化と一致する。
2. 高質量: $\chi \sim 0.7$ の高速回転成分の存在は、階層的合体（2 世代目ブラックホール）または化学的に均一な進化を示唆する。データは、高質量系がより広範で高速回転するこの成分によって支配されるという図を支持し、これは対不安定質量ギャップを埋める可能性がある。

本論文は、有効スピン単独よりも、個々の成分スピン（$p(\chi_1, \chi_2 | m_1)$）の質量依存構造の方が、形成チャネルを区別するためのより鋭い診断を提供すると強調する。著者は、結果が高密度環境（球状星団、活動銀河核円盤）における階層的合体および化学的に均一な進化と整合的である一方で、現在のデータはこれらの特定のシナリオを一意に区別しないことに言及する。彼らは、これらの制約をさらに鋭くし、これらの爆発的事象の不確実な起源を解決するためには、イベント数の増加を伴う将来のカタログが必要であると結論づける。