Gravitational-wave signatures of primordial black hole clusters and the imprint of an early dense-core collapse

技术摘要：原初黑洞星团的引力波特征与早期致密核塌缩的印记

问题陈述
理论上预计原初黑洞（PBHs）的自旋极小（$\chi \simeq 0$），且如果原初曲率涨落具有非高斯尾部，它们将以强聚类构型而非泊松分布的形式形成。当前的引力波（GW）观测（如 GWTC-5.0）表明，存在一个具有显著自旋（$\chi \sim 0.7$）的合并黑洞亚群，这可以作为等级化组装或第二代合并的潜在信号。然而，在模拟 PBH 星团时存在一个根本性的张力：产生观测到的自旋水平需要极高的初始密度（亚秒差距尺度）以促进频繁合并，但这类致密构型由于快速蒸发而无法在宇宙学时间尺度上生存。相反，秒差距尺度的星团可以生存，但产生的自旋极微。本文研究了是否可以通过两阶段动力学演化——即一个短寿命、超致密的核嵌入在一个长寿命、扩展的晕之中——来解决这一张力，并重现 PBH 星团的多信使特征。

方法论
作者使用修改后的 nbody6++gpu 代码进行直接 $N$ 体模拟，以模拟 PBH 星团的动力学演化。为了解决物理尺度巨大的差异（长度差异超过 6 个数量级，时间差异约 7 个数量级），研究被分为两个独立的模拟模块：

1. 塌缩核模块（Collapsing-Core Block）： 模拟一个包含 4,000 个 PBH 的紧凑子系统，其 Plummer 半径（$R_{\text{Plum}}$）范围为 $3 \times 10^{-6}$ 至 $10^{-4}$ pc。初始速度设为零以保证立即发生维里塌缩。测试了三种截断对数正态质量函数（2–50 $M_\odot$, 10–50 $M_\odot$, 以及 10–80 $M_\odot$），以探测高质量尾部的影响。代码经过修改，用于处理近抛物线动力学捕获合并，并利用基于数值相对论结果校准的启发式处方，根据特定无量纲能量、线性动量和轨道角动量来确定残余自旋和反冲速度。
2. 长期晕模块（Long-Term Halo Block）： 模拟一个包含 $2 \times 10^4$ 个 PBH 的秒差距尺度星团（$R_{\text{Plum}} = 10$ pc），跨越一个哈勃时间（13.8 Gyr）。该模块使用了预测的完整 PBH 能谱（峰值在 $2 M_\odot$，并在 $50 M_\odot$ 附近有一个次级特征），并包括包含或不包含质量分别为 1,000、1,500 或 2,000 $M_\odot$ 的中心中质量黑洞（IMBH）的情景。

主要贡献与结果

• 解决自旋-密度张力： 模拟证实，要在 $30\text{--}50 M_\odot$ 范围内产生均方根自旋 $\chi_{\text{rms}} \sim 0.1\text{--}0.2$，需要初始星团半径在 $10^{-5}$ pc 数量级。此类密度会导致快速解体，使得单尺度星团无法在宇宙学时间内生存。作者提出了一个“核-晕”方案，即一个致密、短寿命的核发生塌缩并合并，从而在残余物上留下自旋印记，而周围扩展的晕则得以生存。
• 核内的自旋生成： 核模块在几天到几年的时间内发生塌缩，产生爆发式的近抛物线动力学捕获（$e \to 1$, $l \ll l_{\text{crit}}$）。由于这些遭遇携带的轨道角动量较低，残余自旋相对于准圆轨道值（$\chi \simeq 0.7$）被系统性地抑制，通常产生 $\chi_{\text{rms}} \sim 0.1\text{--}0.2$。
  • 如果初始质量函数在 $50 M_\odot$ 处截断，则高于此质量的存活群体完全由合并残余物组成，导致 $\chi_{\text{rms}} \sim 0.4\text{--}0.5$。
  • 如果能谱延伸至 $80 M_\odot$，则在该范围内会存留无自旋的第一代 PBH，从而稀释平均自旋。
  • 这表明，观测数据中支持在 $\sim 45 M_\odot$ 以上具有快速自旋黑洞的现象，意味着原初质量函数在 $50 M_\odot$ 附近截断。
• 失控增长与 IMBH 形成： 每个核的实现过程都会在初始塌缩后经历失控增长，从而形成一个 IMBH。IMBH 的质量饱和在约 $0.2 M_{\text{cl},0}$（其中 $M_{\text{cl},0}$ 是初始核质量），且独立于具体的质量函数细节。
• 长期演化与双星逃逸者： 秒差距尺度的晕模块在一个哈勃时间内产生的星团内合并少于 5 次，证实了仅靠长寿命星团本身无法产生旋转亚群。主要的观测输出是“双星逃逸者”。
  • 星团内的动力学形成高度偏向于最重的成员；尽管初始质量函数峰值在 $2 M_\odot$，但逃逸者由质量 $>10 M_\odot$ 的天体主导，峰值在 $30\text{--}50 M_\odot$。
  • 逃逸者的丰度在初始分布的截断质量处急剧下降，提供了与基于自旋的诊断一致的基于质量的诊断。
  • 在星团内形成的双星如果变得足够紧密（$a \lesssim 1$ AU），会因硬化遭遇中的能量释放而被弹出。因此，星团无法保留那些在哈布尔时间内合并的双星；只有被弹出的双星才会圆化并进入 LIGO-Virgo-Kagra (LVK) 检测频段。
• 中心 IMBH 的影响： 被保留的 $2,000 M_\odot$ 中心 IMBH 通过稳定星团对抗蒸发并阻止周围 PBH 形成硬双星，几乎完全抑制了双星形成。较轻的 IMBH（$1,000\text{--}1,500 M_\odot$）经常在最初的几 Gyr 内被弹出（在 6 个实现中分别占 5/6 和 2/6），之后星团的行为表现得如同不存在 IMBH 一样。

随机引力波背景 (SGWB)
本文识别了两个成分的 SGWB：

1. 高红移核成分： 来自辐射时期塌缩核的合并事件，经红移后在 $10^{-6}\text{--}10^{-4}$ Hz 波段产生随机背景，可由 LISA 探测。峰值频率和振幅取决于初始核密度（更致密的核峰值频率更低）。
2. 低红移晕成分： 由晕的长期演化产生，其频谱随频率上升，并由双星逃逸者在 $10^{-5}$ Hz 以上主导。如果保留了质量巨大的中心 IMBH，它会通过涉及 IMBH 的宽轨道双星贡献一个低频凸起（$10^{-7}\text{--}10^{-5}$ Hz）。

意义与主张
本文声称提供了一个将 PBH 星团的自旋分布、质量分布和随机背景联系起来的一致动力学框架。它指出：

• 观测到的自旋亚群需要一个早期、紧凑的动力学阶段（即“核”），该阶段与负责当前合并探测的长寿命“晕”是不同的。
• 自旋和质量数据均独立支持原初质量函数在 $50 M_\odot$ 附近截断。
• 观测任务的分工是明确的：自旋由原初塌缩设定，而可探测的合并则由长期的晕动力学决定。
• 作者将此作为两阶段演化模型的“概念验证”。他们明确指出了一些局限性，包括理想化的初始条件（零速度）、对核心和晕模块的独立处理，以及合并处方的启发式性质。他们并未声称确定证明了 PBH 星团的存在，而是展示了这种情景如何能够同时满足多个难以用单尺度模型调和的观测约束。