HydroGrav: Precise hydrodynamics and gravitational waves for cosmological phase transitions

“HydroGrav：用于宇宙学相变的高精度流体力学与引力波研究”技术摘要

问题陈述
探测来自宇宙学一阶相变的随机引力波背景（SGWB）是未来空间观测站（如 LISA）的主要目标。该信号的主要贡献来自于由破缺相气泡膨胀产生的声波。然而，预测由此产生的引力波谱需要对流体动力学演化进行精确建模。目前的文献主要依赖于简化的状态方程（EOS）——具体为“袋模型”（bag model）和 $\mu\nu$ 模型（改进的袋模型）——这些模型假设声速恒定（$c_s^2 = 1/3$）或在真空相中具有固定的声速。这些近似忽略了直接源自底层粒子物理有效势的温度依赖热力学性质。作者认为，在流体力学中使用简化的 EOS，同时使用基于这些简化模型推导出的拟合公式来预测引力波，会引入概念上的不一致性。此外，这些近似在确定流体动力学模式（爆燃、爆轰或混合模式）以及谱形状方面的有效性，尚未通过与全模型依赖的计算进行严格测试。

方法论
本文介绍了 HydroGrav，这是一个旨在利用直接源自有限温度有效势的精确状态方程来计算自相似流体剖面和引力波谱的 C++ 代码。该框架运行如下：

1. 有效势与热力学： 代码直接从有效势 $V_{\text{eff}}(\phi, T)$ 计算热力学量（压力 $p$、能量密度 $e$、焓 $w$、熵 $s$），而非假设为辐射主导形式。这使得能够实现全模型依赖的声速 $c_s^2(T) = dp/de$ 计算。
2. 相对论流体力学： HydroGrav 在自相似坐标系（$\xi = r/t$）中求解完美流体的相对论流体力学运动方程。它为三种流体动力学模式构建了流体剖面（速度 $v(\xi)$、焓 $w(\xi)$、$T(\xi)$）：爆燃（deflagration）、爆轰（detonation）和混合（hybrid）模式。
   • 对于简化 EOS（袋模型和 $\mu\nu$ 模型），代码使用代数匹配条件来确定冲击波前沿和壁速度。
   • 对于精确 EOS，代码采用数值寻根算法（有界二维牛顿法）来满足气泡壁和冲击波前沿处的匹配条件，因为热力学量的代数形式是未知的。
3. 引力波计算： HydroGrav 没有使用基于有限模拟校准的拟合公式，而是实现了声壳模型（Sound Shell Model, SSM）。这种半解析方法直接利用流体动力学解生成的速度和能量密度扰动剖面来构建引力波谱。这保持了微观热力学与宏观信号之间的联系，捕捉到了在破碎幂律拟合中丢失的谱特征（例如双峰结构）。
4. 验证与实现： 代码通过了针对数值求解器（黄金分割最小化、牛顿法、RK4 积分）的单元测试验证。通过使用 $Z_2$ 对称的标量单态扩展标准模型（xSM）进行了演示，其中有效势是通过 DRalgo 包利用三维有效场论（3dEFT）构建的。

核心贡献

• HydroGrav 框架： 一种新的开源 C++ 工具，通过使用精确热力学的精确状态方程，连接了粒子物理模型与引力波预测，支持两种 EOS 处理方式。
• 精确 EOS 实现： 一种稳健的数值方法，用于求解使用精确状态方程的流体剖面，避免了袋模型和 $\mu\nu$ 模型中存在的恒定声速假设。
• 声壳模型集成： 一个一致的流水线，使用相同的流体剖面进行流体动力学演化和引力波源计算，消除了复杂的流体力学与简化谱拟合之间的不匹配。
• 系统性比较： 对 xSM 参数空间进行了全面的扫描，比较了不同 EOS 模型下的流体剖面和引力波谱。

结果
作者对 xSM 模型进行了参数扫描（1,037 个点），并针对爆燃、混合和爆轰模式分析了特定的基准点：

• 流体动力学模式： $\mu\nu$ 模型通常比袋模型更准确地重现了精确 EOS 所预测的流体动力学模式（爆燃 vs 爆轰）。然而，差异依然存在；袋模型在约 5% 的采样点中预测错误，常将爆轰误判为混合或爆燃模式。$\mu\nu$ 模型表现出极高的一致性，失败率 $\lesssim 0.02\%$。
• 流体剖面： 与辐射主导假设偏离最大的部分出现在破缺相中，此时热力学量的温度依赖性导致声速偏离 $c_s^2 = 1/3$。这在温度剖面和冲击波前沿的位置上表现得最为明显。
• 引力波谱：
  • 谱形状： 声壳模型揭示了一个编码了平均气泡间距和声壳厚度的双峰结构，而单一的破碎幂律拟合公式无法捕捉到这一点。
  • 振幅与形状差异： 虽然 $\mu\nu$ 模型通常比袋模型更接近精确 EOS，但显著的峰值振幅和谱形状差异（通过 Wasserstein 距离衡量）仍然存在，特别是在混合解以及声速显著偏离辐射极限的区域。
  • 峰值振幅： 对于 xSM 模型，$\mu\nu$ 近似在约 76% 的爆燃点、50% 的混合点和 47% 的爆轰点中，能将峰值振幅捕捉到精确结果的 5% 以内。
• 信噪比 (SNR)： 利用精确 EOS，作者估计了 LISA 在 xSM 参数空间中的 SNR。该特定模型的预测信号仍低于名义检测阈值（$SNR \approx 1$），但作者指出这是对该流水线的演示，而非对该模型的最终排除，考虑到理论上的不确定性。

意义与主张
论文声称 HydroGrav 为实现“精密引力波宇宙学”提供了必要的一步，建立了一个完全模型依赖的流水线。作者强调：

1. 一致性： 使用精确 EOS 进行流体力学计算，同时保留简化的拟合公式，会产生概念上的不一致；声壳模型通过将相同的流体剖面用于动力学和信号生成，解决了这一问题。
2. 近似的有效性： 虽然 $\mu\nu$ 模型比袋模型有了显著改进，但简化的 EOS 并非普遍可靠。它们可能导致对流体动力学过程的定性描述错误（例如错误的流体动力学模式），或导致定量不同的引力波预测，特别是在声速强烈偏离辐射极限的参数空间区域。
3. 未来的必要性： 随着像 LISA 这样的观测站接近探测电弱标度相变所需的灵敏度，对流体力学进行完全模型依赖的处理将变得越来越重要，以建立观测与底层粒子物理之间的稳健联系。

作者谦虚地总结道，其工作展示了此类精确处理的能力与必要性，而非声称所考虑的特定 xSM 模型是确定可探测的，或断言简化的模型总是无效的。该框架被呈现为一种量化现有近似方法有效范围的工具。