Accretion of Primordial Black Holes in Stellar Interiors

技术摘要：恒星内部原初黑洞的吸积作用

问题陈述
原初黑洞（PBHs）仍是极具竞争力的暗物质候选者，特别是在受限较少的亚月球质量范围（$10^{-16} M_\odot \lesssim M_{BH} \lesssim 10^{-10} M_\odot$）内。恒星为探测该质量范围内的 PBH 提供了一个独特的实验室，因为嵌入在恒星核心中的捕获 PBH 会吸积周围的气体，从而可能改变恒星的演化。先前的研究（例如 Bellinger 等人，2023）使用一维恒星演化代码对这一场景进行了建模，但依赖于源自薄吸积盘理论的假设。具体而言，它们假设了固定的辐射效率 $\eta \approx 0.08$，并将爱丁顿光度极限应用于所有质量。然而，对于该质量范围内的 PBH，邦迪半径（$r_B$）比恒星圆化半径小几个数量级，导致无法形成吸积盘。相反，吸积过程是球对称的，气体在接近史瓦西半径（$r_S$）时被压缩并加热至相对论性温度（$T \sim 10^{10}–10^{12}$ K）。在这种状态下，微观物理过程（相对论性制动辐射、对产生过程、克莱因-尼尔森散射）以及流体的光学厚度与薄盘范式有着本质的区别，因此必须对吸积率（$\dot{M}$）、辐射效率（$\eta$）和流体结构进行自洽的计算。

方法论
作者求解了带有隐式冷却源项的随时间变化的球面欧拉方程，以模拟嵌入类太阳恒星核心的 PBH 的吸积流。

• 微观物理： 冷却项 $\epsilon(\rho, T)$ 包括相对论性制动辐射（电子-离子和电子-电子）以及电子-正电子对湮灭。不透明度考虑了汤姆逊散射、自由-自由吸收以及高能光子的克莱因-尼尔森抑制。
• 数值方法： 作者利用开源 Python 包 radbondi，将流体从初始的绝热邦迪剖面演化至自洽的稳态。这种随时间演化的方法是必要的，因为在冷却主导的机制下，声速点变成了一个焦点（复特征值）而非鞍点，这使得标准的稳态 ODE 射击法失效。
• 机制处理：
  • 光学薄区域（$M_{BH} \lesssim 5 \times 10^{-13} M_\odot$）： 欧拉求解器直接计算物理辐射效率和吸积率。
  • 光学厚区域（$M_{BH} \gtrsim 5 \times 10^{-13} M_\odot$）： 邦迪球变得光学厚，进入光子俘获机制。作者采用了 Begelman (1979) 的解析处方，其中位于俘获半径（$r_{trap}$）内部局部产生的光子会被带回黑洞，而逃逸的光度则受到自调节，满足 $L_{esc} \lesssim 0.6 L_{Edd}$。
• 反馈机制： 研究自洽地模拟了辐射反馈，即沉积在邦迪球之外的光度会加热环境等离子体。这种加热通过对流输运（混合长度理论）得到缓解，从而防止了纯辐射扩散所预测的极端温度升高。

关键结果
研究识别了三种不同的球面吸积机制，并得出以下结果：

1. 三种吸积机制：

   • 热邦迪机制（Hot Bondi, $M_{BH} \lesssim 10^{-14} M_\odot$）： 冷却在动力学上可以忽略不计。流体接近绝热状态，在 $r_S$ 附近达到 $T \sim 10^{11}$ K 的温度。内层流体在达到对产生温度之前就已变为无碰撞状态，从而抑制了对产生平衡。制动辐射在冷却中占主导地位。
   • 制动辐射冷却（Bremsstrahlung Cooling, $10^{-14} \lesssim M_{BH} \lesssim 5 \times 10^{-13} M_\odot$）： 冷却变得高效，在气体达到对产生温度之前便排出了热能。流体趋向于等温状态（$T \approx T_\infty$），吸积率相对于绝热邦迪吸积率增强了约 $\sim 2–7$ 倍。
   • 光子俘获（Photon Trapping, $M_{BH} \gtrsim 5 \times 10^{-13} M_\odot$）： 邦迪球具有光学厚度。虽然局部发射的光度随质量迅速上升，但逃逸的光度被 Begelman 极限限制（$L_{esc} \lesssim 0.6 L_{Edd}$）。至关重要的是，吸积率保持在接近邦迪值的水平（$\dot{M} \propto M_{BH}^2$），因为对于流入过程不存在球面爱丁顿极限；辐射压力仅调制有效引力。

2. 辐射效率： 自洽计算得到的辐射效率在光学薄范围内为 $\eta \approx 10^{-2}$，比此前文献中假设的典型薄盘值（$\eta \approx 0.08$）低了一个数量级。效率在 $10^{-15} M_\odot$ 附近达到峰值，随后在冷却和等温极限下下降。

3. 增长轨迹： 由于不存在球面爱丁顿极限，PBH 的增长在整个光子俘获机制中保持为超指数增长（$\dot{M} \propto M_{BH}^2$）。增长的瓶颈完全由最低质量阶段缓慢的“热邦迪”阶段决定。

4. 临界初始质量： 使 PBH 能在哈勃时间内消耗掉一颗太阳质量恒星的临界初始质量为 $M_{0,crit} \sim 10^{-16} M_\odot$。这比此前基于高效率和爱丁顿瓶颈的估计值（$\sim 10^{-11} M_\odot$）低了五个数量级。

5. 磁场： 作者发现磁通量的积累可以忽略不计。太阳核心的背景磁场太弱，无法在动力学上阻碍流体，且磁重联将任何磁通量积累限制在远低于动力学重要性的水平（$P_{mag}/P_{ram} \ll 1$）。

意义与主张
本文声称提供了首个关于跨越亚月球质量范围的球面 PBH 吸积辐射效率和吸积率的自洽计算。通过使用微观物理处理取代薄盘假设，作者证明了：

• 辐射效率显著低于此前假设的水平。
• 对于球面吸积而言，不存在爱丁顿限制阶段；在质量通量方面，流体保持为超爱丁顿状态，而不受辐射压限制。
• 因此，来自恒星观测对 PBH 暗物质的约束比此前认为的要严格得多，因为初始质量低至 $10^{-16} M_\odot$ 的 PBH 也能在哈勃时间内增长并吞噬其宿主恒星。

作者指出，其工作侧重于吸积的微观物理，尚未包括对全局恒星结构的反馈作用（这需要耦合到如 MESA 之类的恒星演化代码），但他们提供了必要的 $\eta(M_{BH})$ 和 $\dot{M}(M_{BH})$ 关系，以便为未来的此类研究提供支持。他们也将关于盘形成（发生在更高质量或旋转更快的恒星中）的讨论留到了另一篇配套论文中。