From cosmological simulations to binary black hole mergers: The impact of using analytical star formation history models on gravitational-wave source populations

技術要約：宇宙論的シミュレーションから連星ブラックホール合体へ

問題提起
連星ブラックホール（BBH）合体の観測は、宇宙論的な時間スケールにおける大質量星の進化に関する重要な知見を提供する。しかし、観測された合体特性をその前駆体となる恒星へと結びつけるには、宇宙の星形成史（SFRH）および化学組成の進化に関する正確なモデルが必要である。現在、この分野では、宇宙論的シミュレーションや観測から得られた複雑な履歴のプロキシとして、金属量依存的な宇宙星形成率密度、$S(Z, z)$ に対する解析的なフィッティング手法が頻繁に用いられている。これらの解析的な処方箋は、多くの場合制約が不十分であり、推論されるBBH集団への影響も十分に理解されていない。さらに、シミュレーションの解像度およびボックスサイズが、導出される $S(Z, z)$ およびそれに続くBBH合体率に与える影響についても、系統的な定量化が行われていない。

手法
著者らは、IllustrisTNGスイート（具体的には TNG50-1、TNG100-1、TNG300-1）を用いて、ガスセルのデータから直接、シミュレーションに基づく $S(Z, z)$ を導出している。これらは解像度とボックスサイズが異なるものである。彼らは、全星形成率密度 $S(z)$ に関する Madau & Dickinson (2014) の関数形式と、歪んだ対数正規分布を組み合わせた、広く用いられている解析的フィッティングモデル（L. A. C. van Son et al. 2023に基づく）と、これらの直接的なシミュレーション結果を比較している。

異なる $S(Z, z)$ 入力がBBH集団に与える影響を評価するため、著者らは両方のモデル（シミュレーション由来および解析的）を、高速連星個体数合成コードである COMPAS (Compact Object Mergers: Population Astrophysics and Statistics) に結合させている。彼らは「宇宙積分」の枠組みを用いて、赤方偏移に伴う合体率と質量分布を計算している。本研究では特に以下の点を調査している：

1. 解析的フィッティングが、異なるTNG解像度における $S(Z, z)$ の形状と規格化をどの程度忠実に再現できているか。
2. 解析的フィッティングとフルシミュレーションデータの両方を用いた場合における、BBH合体率および一次質量分布 ($M_{BH,1}$) の偏差。
3. シミュレーションの解像度（高解像度のTNG50-1から低解像度のTNG300-1まで）がこれらの結果に与える影響。

主な結果
本研究は、$S(Z, z)$ に対する解析的な処方箋が、BBH集団をモデリングする際に重大な系統誤差を導入することを見出した：

• $S(Z, z)$ における不一致： 解析的フィッティングは、シミュレーションベースの金属量分布に見られる主要な構造的特徴を捉え損ねている。具体的には、シミュレーションに存在する高金属量領域（$Z \gtrsim 3 Z_\odot$）における二次的な「バンプ（隆起）」と、平坦化した低金属量テールを捉えられていない。このフィッティングは、高赤方偏移（$z \gtrsim 6$）かつ高金属量における星形成率を過大評価する傾向がある一方で、低金属量および極めて高い赤方偏移における寄与を過小評価しており、これはシミュレーションの解像度に依存している。
• 合体率への影響： 高赤方偏移（$z \gtrsim 6$）において、解析的フィッティングは、シミュレーションベースの合体率と比較して、10倍から$10^4$倍の範囲でBBH合体率を過大評価する。この誤差の大きさは、シミュレーションの解像度が低下するにつれて増大する。低赤方偏移（$z \approx 0$）においては、フィッティングから導かれる局所的な合体率は、概ねシミュレーションベースの合体率と因子 $\sim 1.5$ 以内で一致しているが、TNG50-1では $\sim 15$ Gpc$^{-3}$ yr$^{-1}$ の過大評価が見られる。
• 質量分布への影響： 解析的フィッティングの使用は、BBHの一次質量分布に偽の構造（スプリアスな特徴）を導入する。最も顕著な点として、$z \lesssim 2$ において、フィッティングは一次質量分布の $8 M_\odot$ 付近に人工的な過剰（バンプ）を生じさせるが、これはフルシミュレーションベースの $S(Z, z)$ を用いた場合には存在しない。この特徴は安定質量移動チャネルに起因しており、高金属量の星形成ピークを正しくモデル化できていないことに由来する。逆に、フィッティングは $M_{BH,1} \gtrsim 35 M_\odot$ において、高質量成分（$33 M_\odot$ および $45 M_\odot$ 付近）を3〜5倍の係数で抑制してしまう。
• 解像度依存性： シミュレーションの解像度は、$S(z)$ の規格化とピーク位置に大きく影響する。高解像度のシミュレーション（TNG50-1）は、低解像度の実行（$z \approx 2.5$）と比較して、より高い星形成率を示し、わずかに低い赤方偏移（$z \approx 2$）でピークに達する。解析的フィッティングの精度は、シミュレーションの解像度が低くなるにつれて低下し、特に高赤方偏移の勾配を再現する能力が低下する。

意義および主張
著者らは、解析的フィッティングは計算上の便宜を提供するものの、宇宙論的シミュレーションから導出された複雑な基礎となる $S(Z, z)$ の信頼できるプロキシとしては不十分であると結論付けている。これは、特に重力波観測を解釈する際に重要である。本研究は、単純な解析モデルが、BBHの合体率や質量分布に直接影響を与える重要な物理的特徴（高金属量のバンプや特定の低金属量テールなど）を捉えきれないことを実証している。

これらの不一致は、現在の解析的フィッティングの関数形式が、シミュレーションに見られる多峰性の金属量分布を収容できないために生じる。したがって、これらのフィッティングに依存することは、BBHの形成チャネルや進化の歴史に関する誤った推論を導く可能性がある。著者らは、重力波データの解釈において偽の構造を導入することを避けるために、宇宙論的シミュレーション、解析的フィッティング、および個体数合成を注意深く統合することを主張している。彼らは、解析的フィッティングが完全に無用であると主張しているのではなく、次世代の検出器が高赤方偏移を探索していく中で、その限界を明らかにし、高忠実度のシミュレーションデータに対して検証を行う必要性を強調している。