Evidence for Three Subpopulations of Merging Binary Black Holes at Different Primary Masses

問題提起
第4次LIGO–Virgo–KAGRA重力波カタログ（GWTC-4.0）のリリースに伴い、確信を持って検出された連星ブラックホール（BBH）合体の数は150を超えた。初期の予想では、少数の検出があれば単一の主要な形成チャネルを特定するのに十分であると考えられていたが、現在のデータは、単一のメカニズムでは容易に説明できない複雑な集団を明らかにしている。スピン・ティルト（スピンの傾き）を用いたサブポピュレーション（副集団）の特定に関する過去の試みは、測定の困難さや理論的な不確実性のために課題に直面してきた。さらに、パラメータ間のグローバルな相関は存在するものの、その天体物理学的な解釈は非自明なままである。ここで対処される中心的な問題は、一次質量（$m_1$）において明確な遷移境界を持ち、質量比（$q$）およびスピン大きさ（$\chi$）の分布が異なる、明確な特性を持つBBHのサブポピュレーションが存在するかどうかである。

手法
著者らは、累積的なGWTC-4.0データ（二次的な曖昧さがあるGW190814、および波形系統誤差や極端なスピンに関する懸念があるGW231123を除く）を用いて、BBH集団をモデル化するための現象論的なアプローチを採用している。彼らは、階層的ベイズ推論を用いて、複雑さのレベルが異なるモデルを比較している。

主な手法の特徴は以下の通りである：

• パラメータ化： 一次質量分布（$p(m_1)$）は、2つのガウスピークを持つ破断型パワーローとしてモデル化される。スピン大きさ（$p(\chi)$）およびティルト（$p(\cos \theta)$）は、それぞれ、切断正規分布および等方的成分と整列成分の混合としてモデル化される。赤方偏移（$p(z)$）はパワーローに従う。
• 遷移質量： コアとなる革新性は、$q$、$\chi$、および $\cos \theta$ の分布が最大3つの遷移質量で形状を変化させることを許容することであり、これにより最大4つのサブポピュレーションを生成している。分布の関数形式は遷移を越えて同じであってもよいが、形状を制御するハイパーパラメータは異なってよいとされる。
• モデル選択： 著者らは、過学習を避けるためにオッカムの剃刀（Occam penalty）を適用しながら、モデル選択を行うためにベイズ因子を利用している。彼らは、「マルチポップ（multipop）」モデル（遷移を許容するもの）を、基準となる強くパラメータ化されたモデル（gwtc4-fiducial）および特定の機能が除去されたバリアントと比較している。

主な貢献と結果
この解析は、2つの遷移質量、$m_t^{AB} = 27.7^{+4.1}_{-3.4} M_\odot$ および $m_t^{BC} = 40.2^{+4.7}_{-3.2} M_\odot$ によって隔てられた、少なくとも3つの明確な連星ブラックホール合生のサブポピュレーションが存在するという強い証拠を提供している。「マルチポップ」モデルは、自然対数ベイズ因子12.3で、基準となるモデルよりも好ましい結果となった。

3つのサブポピュレーションは以下のように特徴付けられる：

1. サブポピュレーションA ($m_1 \lesssim 27.7 M_\odot$):

   • 質量比： ほぼ平坦な分布を示す（$\beta_0 \approx -0.8$）。
   • スピン： 小さなスピン大きさを特徴とし、主に $\chi \leq 0.5^{+0.1}_{-0.1}$ に支持される。
   • 解釈： 効率的な角運動量輸送を伴う孤立連星を含む、複数の形成チャネルからの寄与を含んでいる可能性がある。

2. サブポピュレーションB ($27.7 M_\odot \lesssim m_1 \lesssim 40.2 M_\odot$):

   • 質量比： 等質量ペアリング（$q \approx 1$）への非常に鋭い嗜好を示す（$\beta_1 \approx 3.6$）。この明確なペアリング分布は、一次質量分布における $\sim 35 M_\odot$ 付近の以前に報告されたピークを解消しており、そのピークがモデルの誤指定によるアーティファクトであったことを示唆している。
   • スピン： スピン大きさとティルトの分布は、サブポピュレーションAと統計的に区別がつかない。
   • 解釈： 著者らは、化学的均質進化（CHE）またはポピュレーションIIIの星を潜在的な起源として示唆しているが、標準的なCHEの予測とスピンの差の欠如を整合させることの難しさについても述べている。スピン・ティルト分布の明確な違いが見られないことから、動力学的形成の可能性は低いと考えられる。

3. サブポピュレーションC ($m_1 \gtrsim 40.2 M_\odot$):

   • スピン： 実質的な $\chi > 0.5$ の支持（$\chi \sim 0.7$ でピーク）を伴う、著しく広いスピン大きさ分布を特徴とする。
   • 質量比： $q \approx 0.5$ への嗜好を仮定する、暫定的な支持を示す。
   • 解釈： その特性は、広範なスピン分布と $q \approx 0.5$ への嗜好から、階層的合体（前世代のブラックホールの合体）と概ね一致している。しかし、このグループ内に小さなスピンが存在することや、スピン・ティルトにおける遷移の欠如は、高密度星団における階層的合体がこのサブポピュレーションを完全には説明できない可能性を示唆している。

意義と主張
本論文は、これらの3つのサブポピュレーションが単一の連続的な集団としてではなく、共存するものとして解釈されることがデータから強く支持されると主張している。質量比とスピン大きさの分布における鋭い遷移境界の特定は、BBHパラメータの集団分布における分離可能性の仮定が $\sim 25 M_\odot$ 以上で崩れることを意味している。

著者らは、彼らの結果が強くパラメータ化されたモデルに依存しており、モデルの誤指定の可能性もあることを強調しているが、信号の統計的な強さと、他の文献分析におけるこれらの遷移質量の独立した確認は、これらのサブポピュレーションの存在が現実であることを示唆している。彼らは、異なる形成経路が異なる質量範囲で支配的である可能性が高いと結論付けており、サブポピュレーションBはCHEまたはポピュレーションIIIの星に関連し、サブポピュレーションCは階層的合体に関連している可能性があるが、代替シナリオ（AGNディスクにおける降着など）も依然として可能であるとしている。本研究は、$10 M_\odot$ のピーク付近には遷移の証拠は見いだせなかったと明記しており、この構造については今後の調査に委ねている。