Implications of the nanoHertz Gravitational-Wave Background for Galactic Feedback and Massive Black Hole Growth

技术摘要：纳赫兹引力波背景对星系反馈与大质量黑洞生长的启示

问题陈述
脉冲星计时阵列（PTA）的最新测量为随机纳赫兹引力波背景（GWB）提供了强有力的证据，该背景很可能由超大质量黑洞双星（SMBHBs）的非相干叠加产生。该背景的振幅对超大质量黑洞（SMBH）的质量函数敏感。然而，存在一种张力：宇宙学流体动力学模拟（如 IllustrisTNG、MillenniumTNG 和 Simba）一致低估了由 PTA 数据推断出的 GWB 振幅。这一差异表明，当前关于 SMBH 生长、种子形成以及调节其生长的反馈机制（活动星系核反馈和恒星反馈）的模型可能是不完整的。具体而言，尚不清楚星系反馈设定的变化如何改变黑洞质量函数（BHMF）的高质量端，进而影响预测的 GWB 信号。

方法论
作者采用基于 Casey-Clyde 等人（2022）的半解析框架，将模拟预测的 BHMF 与 GWB 振幅联系起来。方法论步骤如下：

1. 模拟套件：本研究分析了四个主要模拟套件：IllustrisTNG、MillenniumTNG、Simba 以及 CAMELS（基于机器学习的宇宙学与天体物理模拟）项目。
2. 盒子尺寸与收敛性：鉴于 GWB 对 BHMF 高质量端的敏感性，作者利用多种盒子尺寸（范围从 25 Mpc/h 到 500 Mpc/h）进行了收敛性研究。他们确定，小于约 35 Mpc/h（针对 TNG）和约 50 Mpc/h（针对 Simba）的盒子会遭受泊松噪声的影响，且对大质量晕（$M_{BH} > 10^9 M_\odot$）的采样不足。
3. 反馈变化：本研究通过以下方式探索星系反馈的影响：
   • 比较不同模拟代码中的基准模型。
   • 分析 CAMELS 套件中的"1P"（单参数）变化，其中控制恒星和活动星系核反馈的 23 个天体物理参数被单独改变。
   • 检查 Simba 反馈的特定变体，其中移除了活动星系核喷流模式或对其进行了修改（例如，仅“SW"模式、"SW + 辐射”模式）。
4. GWB 计算：模拟得到的 BHMF 通过爱丁顿比率分布函数（ERDF）转换为双星种群，以识别活跃类星体。该种群通过归一化因子（$N_{cal} \approx 4$）进行缩放，以考虑双星与类星体的比率，并通过质量比分布分配次级质量。将生成的双星种群进行积分，以计算特征应变谱以及在参考频率 $f = 1 \text{ yr}^{-1}$ 处的 GWB 振幅。

主要贡献

• 量化反馈影响：本文证明，恒星和活动星系核反馈模型的选择会显著改变 BHMF 的高质量端。低效或缺失的活动星系核反馈会导致大质量 SMBH 过剩，使预测的 GWB 振幅增加 2 倍（在基准变化中）至 10 倍（在极端的 CAMELS 变化中）。
• 收敛性分析：作者确立了稳健的 GWB 预测所需的最小盒子尺寸，指出在标准的 25 Mpc/h CAMELS 盒子中，$M_{BH} > 10^9 M_\odot$ 的统计量尚未收敛。
• 系统性低估：研究证实，来自 IllustrisTNG、MillenniumTNG 和 Simba 的基准模型低估了 PTA 推断的 GWB 振幅，倍数约为 2。为了与观测结果匹配，这些模型需要在 $M_{BH} > 10^8 M_\odot$ 范围内拥有 1.5 到 5.5 倍更多的 SMBH。

结果

• 基准模型：与 NANOGrav 15 年数据相比，所有标准（基准）模拟模型均低估了 GWB 振幅。这意味着当前的反馈设定对 SMBH 生长的抑制过于有效，或者种子形成发生得太晚。
• 反馈变体：减少或缺失活动星系核反馈的模拟（特别是没有喷流模式的 Simba 运行）产生的 GWB 振幅与 PTA 数据一致。然而，这些模型无法重现真实的星系种群（例如，它们未能有效淬灭恒星形成），因此缺乏物理现实性。
• 参数敏感性：在 CAMELS 1P 集合中，特定参数驱动了 GWB 振幅：
  • 对于 CAMELS-TNG，参数 $\beta$（控制低吸积和高吸积活动星系核状态之间的切换）影响最大。$\beta=0$ 的值与 PTA 数据匹配，但低估了低质量恒星质量函数。
  • 对于 CAMELS-Simba，种子质量比参数（$A_{BHseed}$）最为关键。降低种子阈值允许更早的 SMBH 生长，从而匹配 PTA 数据，但这可能会扰乱早期的恒星形成。
• 恒星反馈：虽然通常假设恒星反馈仅影响低质量黑洞，但研究发现，某些恒星反馈参数（例如，风速度归一化和热分数）也会通过调节黑洞生长可用的燃料，影响高质量 BHMF 和 GWB 振幅。

意义与主张
本文主张，PTA 测量提供了一种强大且独立的星系反馈物理探针，是对传统电磁约束的补充。模拟预测与 PTA 数据之间的不匹配表明，早期宇宙中的 SMBH 生长比当前模型所捕捉到的更为高效或发生得更早。

作者得出结论，虽然修改反馈参数或种子设定可以人为地将模拟与 PTA 数据对齐，但此类修改往往会破坏与观测到的星系统计量（如恒星质量函数）的一致性。因此，GWB 振幅与模拟预测之间的张力表明，有必要重新评估宇宙学流体动力学模拟中的 SMBH 种子形成和生长设定。研究强调，未来的模型必须同时重现真实的星系种群和观测到的 GWB 振幅，这表明当前的反馈模型可能过度抑制了 SMBH 的生长。作者指出，他们的框架依赖于简化的假设（例如，忽略环境硬化和并合时间的细节），未来的精度研究需要对双星演化进行自洽处理。