RABBITS - III. Modelling relativistic accretion discs around spinning black holes in galaxy formation simulations

这是一份关于论文《RABBITS – III. MODELLING RELATIVISTIC ACCRETION DISCS AROUND SPINNING BLACK HOLES IN GALAXY FORMATION SIMULATIONS》（RABBITS 系列第三部分：星系形成模拟中旋转黑洞周围的相对论吸积盘建模）的详细技术总结。

1. 研究背景与问题 (Problem)

• 现有模拟的局限性： 在现有的星系形成流体动力学模拟中，黑洞吸积通常采用亚网格模型（sub-grid models），最常用的是基于 Bondi-Hoyle-Lyttleton (BHL) 公式的随机吸积方案。这些方法存在显著的随机性，且通常假设吸积盘是几何薄盘但忽略相对论效应（如使用 Shakura & Sunyaev 1973 的牛顿近似），导致无法准确描述旋转黑洞附近的能量耗散、角动量传输以及吸积盘的真实结构。
• 物理过程的缺失： 经典方法难以自洽地演化黑洞的质量、自旋矢量（Spin vector）以及吸积盘的几何结构。特别是，黑洞自旋对吸积效率（Radiative efficiency）和反馈能量有决定性影响，而现有的简化模型无法捕捉这种复杂的相互作用。
• 观测对比的困难： 由于缺乏对吸积盘大小、温度分布和光谱能量分布（SED）的物理建模，模拟结果难以直接与观测到的类星体吸积盘半光半径（half-light radii）和光谱特征进行对比。

2. 方法论 (Methodology)

该研究开发并实现了一个几何薄相对论吸积盘模型，并将其集成到 RABBITS 系列的星系形成模拟中。该方法结合了三个核心组件：

1. 相对论吸积盘物理： 基于 Novikov & Thorne (1973) 和 Compère & Oliveri (2017) 的理论。
2. Ketju 积分器： 用于解析黑洞周围的小尺度动力学（Rantala et al. 2017, 2020）。
3. SPH 星系形成模拟： 使用 Gadget-3 和 Sphgal 代码进行的平滑粒子流体动力学模拟。

核心建模步骤：

• 吸积盘结构的构建（分段求解）：

  • 模型不假设单一的吸积盘状态，而是根据局部物理条件（压力主导机制和不透明度主导机制）将吸积盘划分为不同的区域（Regimes）。
  • 压力机制： 气体压力 ($p_{gas}$) 或辐射压力 ($p_{rad}$)。
  • 不透明度机制： 电子散射 ($\kappa_{es}$) 或自由 - 自由吸收 ($\kappa_{ff}$)。
  • 由此形成三种局部解：Gas-ES（气体压力 + 电子散射）、Rad-ES（辐射压力 + 电子散射）、Gas-FF（气体压力 + 自由 - 自由吸收）。
  • 此外，还结合了 Mummery & Balbus (2023) 的解析解，以平滑连接 ISCO（最内层稳定圆轨道）内部（intra-ISCO）和外部（extra-ISCO）的区域，打破了传统“零力矩”边界条件的限制。
  • 通过求解不等式确定各区域的适用范围，将这些局部解拼接成全局吸积盘结构（表面密度、温度、不透明度剖面）。

• 黑洞与吸积盘的自洽演化：

  • 角动量对齐： 利用 King et al. (2005) 的判据，通过积分计算吸积盘的总角动量，判断吸积盘与黑洞自旋是共转（co-rotating）还是逆转（counter-rotating），从而动态调整 ISCO 半径和吸积盘结构。
  • 质量与自旋演化： 黑洞通过吸积盘物质增加质量，同时吸积物质带来的角动量改变黑洞的自旋大小和方向（Bardeen-Petterson 效应）。模型显式计算了这些力矩，而非使用解析近似。
  • 吸积率限制： 吸积率受爱丁顿极限限制，并考虑了自旋依赖的辐射效率。

• 反馈机制：

  • 采用热/类星体模式反馈（Thermal/Quasar mode）。
  • 辐射效率 ($\epsilon_r$) 是自旋和吸积盘结构的函数，不仅包含 ISCO 处的结合能释放，还包含了 ISCO 处额外力矩产生的辐射。
  • 反馈能量沉积率由 $\dot{E} = \epsilon_f \epsilon_r \dot{M} c^2$ 计算，其中 $\epsilon_f$ 为耦合效率。

3. 关键贡献 (Key Contributions)

1. 首个自洽的相对论吸积盘模型： 首次在星系形成模拟中实现了基于第一性原理（First Principles）的吸积盘结构建模，消除了 Bondi 型吸积的随机性。
2. 动态演化黑洞自旋与吸积盘： 模型能够同时演化黑洞的质量、自旋矢量以及吸积盘的几何结构（包括共转和逆转情况），并显式计算 Bardeen-Petterson 效应导致的自旋进动和对齐。
3. 可观测量的直接提取： 模型能够直接输出吸积盘的半光半径（Half-light radius）、表面密度、温度剖面和光谱能量分布（SED），使得模拟结果可以直接与观测数据（如微透镜观测、反响映射）进行对比。
4. 开源工具发布： 发布了 C 语言版本（用于集成到模拟中）和 Python 版本（用于后处理和构建虚拟吸积盘），以及更新版的 Relagn 模型。

4. 主要结果 (Results)

• 参数空间探索：

  • 粘度 ($\alpha$)： 对吸积盘的区域划分影响较小，但显著影响表面密度和盘的稳定性（高粘度导致更大的盘半径）。
  • 黑洞质量 ($M_\bullet$)： 质量越大，吸积盘起始半径（光子半径）和终止半径（自引力半径）越大，且盘更倾向于由气体压力和自由 - 自由吸收主导（Gas-FF 区域扩大）。
  • 吸积率 ($\dot{M}$)： 高吸积率会导致辐射压力主导（Rad-ES 区域出现），并可能引发盘的不稳定性。
  • 自旋 ($\alpha_\bullet$) 与取向： 共转和逆转显著改变 ISCO 半径和吸积盘结构。逆转会导致更冷的盘和更小的有效吸积区域。

• 孤立星系模拟结果：

  • 吸积盘质量演化： 模拟展示了气体粒子被捕获、形成吸积盘并逐渐向内迁移的过程。吸积盘质量呈现非线性的演化特征。
  • 黑洞生长与自旋： 比较了“相干吸积”（Coherent，吸积盘角动量继承自捕获气体）和“混沌吸积”（Chaotic，随机方向）两种模式。
    • 混沌模式下，由于吸积盘尺寸较大，导致更多吸积事件被判定为共转（Co-rotating），从而使得黑洞自旋增加更快。
    • 自旋增加导致辐射效率提高，反馈增强，反而抑制了气体的进一步吸积，最终导致混沌模式下的黑洞质量略低于相干模式。
  • 半光半径预测： 模拟预测的吸积盘半光半径与观测数据（如 Morgan et al. 2010, Blackburne et al. 2011 等）相比，系统性地偏大（约几倍），这与标准薄盘理论预测偏小但观测偏大的“长期争议”一致，但模型提供了更物理的物理解释（如自引力半径的确定）。

• 光谱能量分布 (SED)：

  • 基于模型构建的温度剖面生成的 SED 与传统的 Novikov & Thorne (1973) 黑体辐射模型相比，在红外、光学和紫外波段（低频尾部）存在显著差异。
  • 新模型在盘的外围区域（低温区）与经典模型差异最大，这直接影响类星体“大蓝包”（Big Blue Bump）的形状和反响映射的时间延迟预测。

5. 意义与展望 (Significance)

• 连接微观物理与宏观演化： 该模型成功地将黑洞吸积盘尺度的相对论物理与星系尺度的流体动力学模拟连接起来，填补了亚网格模型与全解析模拟之间的空白。
• 多信使天文学的基石： 通过更准确地捕捉黑洞自旋演化、吸积效率及反馈能量，该模型为解释引力波事件（黑洞并合后的反冲速度、自旋分布）和电磁波观测（类星体光谱、光变）提供了更可靠的理论基础。
• 解决观测与理论的长期矛盾： 通过显式计算吸积盘结构，模型能够生成更真实的 SED 和尺寸预测，有助于解决长期存在的吸积盘尺寸观测值与理论预测不符的问题。
• 未来方向： 作者计划在未来工作中引入辐射效率低下的吸积态（如 ADAF）、喷流模式反馈，并将该模型扩展至双黑洞吸积系统，以完整追踪超大质量黑洞从形成到并合的演化过程。

总结： 这篇论文通过引入一个基于相对论流体动力学和局部解析解的复杂吸积盘模型，显著提升了星系形成模拟中黑洞物理的保真度，为理解黑洞与宿主星系的共同演化以及解释多信使观测数据提供了强有力的新工具。