Nonlinear growth and amplification of phase-transition gravitational waves induced by cosmic expansion

技术摘要：宇宙膨胀诱导的相变引力波的非线性增长与放大

问题陈述
由宇宙首次阶相变（FOPTs）产生的引力波（GWs）是 LISA、Taiji 和 TianQin 等未来观测站以及当前脉冲星计时阵列（PTA）的主要观测目标。虽然在热 FOPTs 过程中产生的声波被认为是随机引力波背景（SGWB）的主要贡献来源，但传统的理论预测通常忽略了一个关键的物理要素：宇宙膨胀。以往考虑了膨胀效应的半解析研究通常预测了信号的抑制以及峰值频率的偏移，这类研究依赖于声波线性叠加的假设。此外，这些方法往往无法解释由于辐射能量密度相对于恒定真空能量随膨胀而稀释所导致的相变强度（$\alpha$）的时间依赖性。对于缓慢相变而言，这种疏忽尤为显著，因为膨胀与成核动力学之间的相互作用可能具有非平凡的特性。

研究方法
为了解决这些局限性，作者在弗里德曼–勒梅特–罗伯逊–沃尔克（FLRW）背景下，进行了首次系统的宇宙首次阶相变三维流体动力学模拟。该研究采用了以下方法论框架：

• 流体动力学框架： 流体动力学使用无希格斯框架下的包方程状态（bag equation of state）进行建模。能量-动量张量在膨胀时空中求解，明确包含了随时间变化的标度因子 $a(t)$。
• 时间依赖的强度： 该方法将热力学量分解为均匀背景分量和相对分量。这揭示了相变强度参数 $\alpha$ 本质上是随时间变化的，由于辐射浴的稀释，$\alpha(t)$ 随标度因子的四次方增长（$\alpha(t) \propto a(t)^4$）。
• 模拟设置： 作者利用具有恒定壁速度（$v_w$）的移动边界来模拟膨胀的气泡，并采用 Kurganov-Tadmor 方案来解析冲击波前沿。模拟在 $512^3$ 网格上进行，箱体尺寸为 $L = 40v_w/\beta$，以确保气泡数量恒定。
• 引力波计算： 通过求解由流体各向异性应力驱动的线性化爱因斯坦方程来计算引力波能量密度。作者追踪了引力波能量份额 $\tilde{\Omega}^*_{GW}(t)$ 以及谱密度 $\Omega^*_{GW}(k)$ 的演化。

主要结果
模拟对宇宙膨胀对引力波产生的影响得出了几个违反直觉且具有重要意义的发现：

1. 非线性放大： 与此前预测的信号抑制相反，宇宙膨胀在引力波产生的早期阶段诱发了高度非线性的、近乎指数级的增长。这导致最终的引力波谱发生了实质性的放大，与非膨胀情景相比，其量级达到了 $O(10)$ 至 $O(100)$。
2. 增长机制： 这种放大是由整个成核过程中有效相变强度 $\alpha(t)$ 的增加驱动的。随着宇宙膨胀，潜热转化为流体动能的效率增加，从而增强了流体的均方根四维速度（$U_f$）。由于引力波信号强度与 $U_f^4$ 成比例，这导致了显著增强的信号。
3. 对相变强度的依赖性： 与中等强度（$\alpha_0 = 0.05$）相比，放大效应在初始较弱的相变（$\alpha_0 = 0.005$）中最为显著。对于弱相变，声波寿命足够长，使得指数增长阶段能够占据主导地位；而对于中等强度的相变，减小的平均气泡间距和较短的声波寿命部分缓解了这种增强效应。
4. 谱偏移： 虽然低频和高频处的谱斜率基本保持不变，但峰值频率向高频移动。这一偏移归因于在膨胀体积中有效增强的成核率导致平均气泡间距减小。
5. 壁速度的作用： 对于初始较弱的相变，放大效应在较低壁速度（爆燃机制）下更为强烈。然而，对于中等强度的相变，由于在膨胀情景下较低速度会抑制声波寿命，这种关系发生了逆转。

意义与主张
本文断言，精确模拟缓慢相变必须包含宇宙膨胀，因为忽略膨胀会导致严重低估引力波信号。作者强调，此处识别出的非线性增长机制是基于线性叠加或静态强度参数的半解析框架无法捕捉到的。

研究得出结论，这种“$O(100)$”的放大是一个由于 $\alpha(t)$ 有效增长而产生的稳健特征，它改变了相变的流体动力学过程。虽然作者承认随时间变化的壁速度（此处视为常数）可能会引入进一步的修正，但他们认为这些效应不太可能消除所观察到的实质性放大。因此，这些结果强调了在膨胀背景下进行全三维流体动力学模拟对于可靠预测来自早期宇宙的引力波（特别是用于解释来自 PTA 和空间干涉仪的潜在信号）的必要性。