解析式恒星形成历史模型如何影响双黑洞并合预测？

用模拟恒星形成检验解析拟合

论文要回答的问题是：常用来描述宇宙恒星形成的简单解析公式，是否足以预测双黑洞（BBH）并合。这里的恒星形成历史指金属丰度相关的恒星形成率密度 S(Z, z)，即在给定金属丰度 Z 和红移 z 下，单位体积、单位时间内的恒星形成量。金属丰度之所以重要，是因为它影响大质量恒星的演化，进而决定它们留下的黑洞。

作者直接从三个质量分辨率和盒子大小不同的 IllustrisTNG 宇宙学模拟中取出 S(Z, z)，分别是 TNG50-1、TNG100-1 和 TNG300-1，并用一个解析形式分别拟合每个模拟：S(z) 采用 Madau–Dickinson 型函数，金属丰度分布采用偏斜对数正态分布。随后，他们把完整的模拟 S(Z, z) 和拟合版本都输入双星族合成程序 COMPAS，并且全程使用同一套 COMPAS 模拟，使两者的差异来自恒星形成输入，而不是双星物理。

这一设计只隔离出一个建模选择。比较的对象不是模拟与观测，而是详细的模拟恒星形成历史与用来替代它的简化拟合。

拟合抓不住金属丰度分布的形状

解析拟合在 10⁻¹ ≲ Z/Z☉ ≲ 10⁰ 附近与模拟吻合最好，但在两端都失效。模拟的 S(Z) 在 Z/Z☉ ∼ 1 附近有一个峰，在 Z/Z☉ ≈ 2–5 附近还有第二个凸起；拟合无法刻画这个第二凸起，只在 Z/Z☉ ≈ 1–2 处给出一个单峰，峰的位置随分辨率变化。在 Z/Z☉ ≳ 6 以上，拟合把 S(Z) 高估最多约 10–1000 倍。

低金属丰度尾部同样被错误表示，而且方式取决于所用模拟。TNG50-1 在 10⁻³ ≲ Z/Z☉ ≲ 10⁻² 之间出现斜率变化，TNG100-1 和 TNG300-1 则在 Z/Z☉ < 10⁻² 以下变平。作者指出，更好的拟合需要更陡的高金属丰度截断，以及至少多一个斜率变化的低金属丰度尾部。

在红移方向上，拟合在 z ≲ 0.5 处系统性地高估 S(z)，最多约 1.5 倍，并且无法再现模拟中 z ≈ 6–8 附近的转折点。在 z ≈ 14 处，拟合对三个模拟都高估 S(z)：TNG300-1 约 10 倍，TNG50-1 和 TNG100-1 约 2 倍。作者部分归因于优化过程更重视恒星形成率峰值附近的高值区域，对高红移这类低速率区域权重较小。

高红移并合率被放大数个量级

把拟合结果输入 COMPAS 后，预测的双黑洞并合率随之改变。在 z ≈ 0 处，使用完整模拟 S(Z, z) 时，三个模拟给出的本地并合率约为 29 至 59 Gpc⁻³ yr⁻¹；使用解析拟合时约为 28 至 74 Gpc⁻³ yr⁻¹。对 TNG50-1，拟合把本地并合率高估约 15 Gpc⁻³ yr⁻¹；对 TNG100-1 和 TNG300-1，模拟与拟合给出的本地并合率在 ±2 Gpc⁻³ yr⁻¹ 内一致。

更大的偏差出现在高红移。论文报告，解析 S(Z, z) 方案会把 z ≳ 6 处的双黑洞并合率高估最多 10–10⁴ 倍，具体倍数取决于模拟分辨率；红移越高、分辨率越低或盒子越大，这一效应越明显。拟合还会改变并合率峰值出现的位置：分辨率降低、盒子增大时，峰值移向更高红移。

本地并合率的比较与高红移的比较回答的是不同问题，不能混为一谈。一个拟合即使对三个模拟中的两个匹配了 z ≈ 0 的归一化，仍可能在 z ≳ 6 处严重偏离。

质量分布中出现虚假的 8 M☉ 特征

拟合不只改变归一化，也改变主质量分布的形状。使用解析拟合时，在 z_merger = 0.2 的主质量分布中，约 8 M☉ 附近出现人为的过量，在 z_merger ≲ 2 时最多可达 7 倍；而使用完整的模拟 S(Z, z) 时没有这一特征。论文认为该特征可能主要来自公共包层通道。与此同时，拟合压制了高质量特征：在 z_merger = 0.2 处，M₁ ≳ 20 M☉ 以上的高质量并合率随主质量增加而被逐步压制，包括 33 M☉ 和 45 M☉ 附近峰值的变平。

作者把原因追溯到金属丰度分布。由于拟合把 3 Z☉ 附近的第二个凸起抹平，它把金属丰度分布的全局极大值放在比模拟更高的金属丰度处，这一偏差随后传递到预测的并合族群中。

有一点保持稳定：对三个模拟而言，主质量分布中各特征的位置在不同红移之间、以及拟合与模拟 S(Z, z) 之间都保持一致。论文由此得出结论：双星建模假设更强地决定特征出现在哪里，而 S(Z, z) 影响总体速率和特征的相对显著程度。

分辨率主要改变归一化而非形状

选择哪一个 IllustrisTNG 模拟，主要改变结果的归一化。TNG50-1 质量分辨率最高、盒子最小，在所有红移上给出的恒星形成率都最大；TNG300-1 则被压低约 2 倍。本地并合率遵循同样的排序：使用模拟 S(Z, z) 时，从 TNG300-1 的约 29 Gpc⁻³ yr⁻¹ 到 TNG50-1 的约 59 Gpc⁻³ yr⁻¹。在并合率峰值附近，TNG50-1 比 TNG300-1 高约 3 倍；在 z ≳ 6 处，差距最多可达两个数量级。

论文对这里的“分辨率”含义保持谨慎。分辨率更高并不自动意味着恒星形成历史更可靠，因为亚网格方案仍然存在，而更小的盒子会增大宇宙方差并限制统计量。作者还指出，在这三个模拟中，分辨率与盒子大小是同时变化的，要分离两者的影响需要更多模拟，他们把这留待后续工作。

与观测到的 GWTC-4 主质量分布的比较也有局限。模型低估了 M_BH,1 ≳ 35 M☉ 处的速率，作者认为这可能需要在本文所建模的孤立双星演化之外，引入动力学环境中的层级并合等额外形成通道。

拟合可用与不可用的边界

论文给出的实际分界在于形状与归一化。TNG50-1 的解析拟合最能准确表示其 S(Z, z) 的形状；而在分辨率更低、体积更大时，速率更低，拟合与模拟的分歧也最大。拟合再现并合率演化最准确的红移区间是：TNG50-1 为 z ≈ 1–4，TNG100-1 为 z ≲ 2，TNG300-1 为 z ≲ 1.5。

因此，选择恒星形成模型是一个随红移变化的判断，而不是单一的“能用或不能用”。一个校准到匹配本地速率的拟合，仍可能扭曲高红移速率和质量特征的显著程度，而这种扭曲的大小取决于由哪个模拟提供底层的 S(Z, z)。

作者明确提出的优先事项是改进 S(Z, z) 建模，而不是直接放弃解析拟合。他们呼吁检验在 IllustrisTNG 中看到的虚假特征是否也出现在其他宇宙学模拟和由观测导出的恒星形成历史中；如果确实存在，那么转向详细恒星形成历史或更新拟合的理由就更充分。在此之前，论文最清楚的警示是：一个拟合能匹配 z ≈ 0 的并合率，并不能证明它的高红移预测或质量分布形状可信。

本文引用的文献