Simulation of PBH formation in a matter-dominated universe

技术摘要：物质主导宇宙中原初黑洞形成的模拟

问题陈述
原初黑洞（PBH）的形成对主导宇宙组分的物态方程高度敏感。虽然辐射主导时期 PBH 的形成过程已得到较充分的理解，但在早期物质主导时期，该过程仍存在不确定性。在物质主导时期，压力梯度可以忽略不计，使得过密区域更容易发生坍缩。然而，形成准则受非球面动力学、不均匀性和椭率度的显著影响。以往的解析估计通常依赖于摄动近似或泽尔多维奇（Zel'dovich）近似，但在高度非线性的引力坍缩机制下，这些近似可能不再适用。研究这一机制的一个主要挑战在于如何在数值模拟中实现非球面对称性，特别是在处理会导致数值崩溃（在黑洞视界完全解析之前就发生壳层交叉奇异性）的尘埃流体时。

方法论
作者采用基于 Baumgarte-Shapiro-Shibata-Nakamura (BSSN) 形式论的全非线性数值相对论，模拟了在膨胀的物质主导背景下的 PBH 形成过程。

• 初始条件： 模拟使用的曲率扰动 $\zeta$ 函数形式包含了椭率度，从而构建了三轴配置。初始数据是利用具有振幅参数 $\mu$ 以及椭率参数 ($e, p$) 的长波梯度展开解构建的尘埃流体。
• 规范与坐标： 采用了一种新型的时间切片条件，即从基于 e-folding 数的规范过渡到在视界进入阶段采用的避免奇异性的“1+log”型切片。使用非笛卡尔坐标系以扩大模拟框的中心区域，从而在坍缩发生的区域提高分辨率。
• 物质描述： 测试了两种不同的物质模型：
  1. 尘埃流体： 一种无压流体描述。作者指出，这种描述不可避免地会导致流体密度发散的壳层交叉奇异性，从而导致数值代码崩溃。
  2. 无碰撞粒子： 为了规避流体近似的失效，作者在全数值相对论中实现了粒子单元法（PIC）。通过为粒子分配质量来重现所需的密度不均匀性，其运动受测地线方程支配。这使得模拟能够在流体元素发生交叉之前继续进行。

关键结果

• 球面对称性验证： 代码通过了针对球对称尘埃坍缩的 Lemaître-Tolman-Bondi (LTB) 解的验证。数值结果成功重现了解析解，直至视界形成点，并正确识别了局部裸奇异性（PBH 形成）与全局裸奇异性之间的阈值。
• 尘埃流体模拟（非球面）：
  • 对于三轴坍缩，由于尘埃流体模拟在视界形成之前就因壳层交叉奇异性而崩溃，除非初始振幅 $\mu$ 显著较大。
  • 研究发现，观察到视界形成之前代码崩溃所需的阈值振幅，大约是此前基于泽尔多维奇近似和箍缩猜想（hoop conjecture）的解析估计值的两倍。
  • 这表明，在流体近似有效的范围内，PBH 的形成比此前认为的更难实现。
• 无碰撞粒子模拟：
  • 当使用无碰撞粒子建模同一系统时，模拟不会因壳层交叉而崩溃。
  • 在初始振幅低至 $\mu \approx 0.045$（对于椭率 $e=0.2$）时，便观察到了黑洞形成（由视界形成指示）。
  • 该阈值比先前的解析估计值小一个数量级，也显著低于从尘埃流体模拟中得到的阈值。
  • 然而，粒子模拟在粒子相交后表现出显著的哈密顿约束违反，这表明虽然坍缩在进行，但由于粒子描述与几何结构之间的不匹配，其演化的定量准确性受到了限制。

意义与主张
本文声称展示了在物质主导时期，PBH 形成准则存在巨大的理论不确定性，这取决于物质场的物理描述：

1. 流体近似极限： 如果将流体描述中与壳层交叉相关的奇异性解释为一种物理抑制机制（或者如果流体近似在坍缩前失效），那么 PBH 形成的阈值将显著高于解析预测，使得 PBH 形成变得更加困难。
2. 粒子动力学： 如果接受无碰撞粒子描述作为单流体近似之外的有效扩展，那么 PBH 的形成将比此前预期的要容易得多（阈值更低）。
3. 理论不确定性： 粒子模拟与流体模拟之间的差异凸显了标准解析准则（泽尔多维奇近似）可能无法准确捕捉三轴坍缩中的非线性动力学。作者得出结论，阈值高度依赖于高密度下对物质场的处理方式。

本文并非旨在最终解决 PBH 丰度问题，而是旨在揭示形成阈值对底层物质模型和数值处理方式的敏感性。作者指出，未来的工作有必要精细化高密度处理方法，改进粒子模拟中的约束保持能力，并量化由此产生的 PBH 丰度和引力波特征。