跳到主要内容
Sandbox Physics
宇宙相变引力波

泡壁怎样推动等离子体

Sandbox Physics ·

泡壁一边推进，一边把真空能变成等离子体的热和整体运动。有多少变成运动，取决于相变强度α\alpha，也取决于泡壁比声速快还是慢。取α=0.2\alpha=0.2、ξw=0.9\xi_w=0.9，单个泡泡的效率是κ=0.274\kappa=0.274，模板里保留的动能占比是K=0.0274K=0.0274。

相变强度与袋模型

把等离子体看成相对论性气体，压强是能量密度的三分之一，声速是cs=1/3≈0.577c_s=1/\sqrt3\approx0.577。袋模型里，假真空比破缺相多带一份恒定的能量密度ϵ\epsilon。相变强度拿它和泡壁前方的辐射能量密度比：α=ϵ/(a+T+4)\alpha=\epsilon/(a_+T_+^4)，a+a_+是定义辐射能量的常数。

泡壁以速度ξw\xi_w穿过远处静止的等离子体。除此之外没有别的尺度，所以单个膨胀泡泡周围的流动只依赖ξ=r/t\xi=r/t，是自相似的。一条剖面解出来，就描述了泡泡在碰到别的泡泡之前的每一时刻。

对应阅读：Espinosa, Konstandin, No & Servant, Energy budget of cosmological first-order phase transitions, arXiv:1004.4187v1 / JCAP 1006:028 (2010) · Caprini et al., Detecting gravitational waves from cosmological phase transitions with LISA: an update, arXiv:1910.13125v2 / JCAP (2020), doi:10.1088/1475-7516/2020/03/024

三种泡壁

比声速慢的泡壁叫爆燃。它推着前方的等离子体走，一道激波以略高于声速的速度跑在前面，激波和泡壁之间的等离子体向外流。

比某个极限速度、即 Jouguet 速度ξJ\xi_J还快的泡壁叫爆轰。前方的等离子体完全没被扰动，泡壁撞上它、拖着它走；后面跟着一道稀疏波，让流体到ξ=cs\xi=c_s时重新慢下来。

速度在csc_s和ξJ\xi_J之间时，解是混合型：前有激波，后有稀疏波。α=0.2\alpha=0.2时 Jouguet 速度是0.8280.828，泡壁速度 0.9 就是爆轰；α=0.01\alpha=0.01时它只有0.6530.653。

ξJ(α)=2α/3+α2+1/31+α\xi_J(\alpha)=\frac{\sqrt{2\alpha/3+\alpha^2}+\sqrt{1/3}}{1+\alpha}
α→0\alpha\to0时它趋于csc_s；相变越强，这条分界越往光速靠。

对应阅读：Espinosa, Konstandin, No & Servant, Energy budget of cosmological first-order phase transitions, arXiv:1004.4187v1 / JCAP 1006:028 (2010) · Caprini et al., Gravitational waves from first-order phase transitions in LISA: reconstruction pipeline and physics interpretation, arXiv:2403.03723v2 / JCAP 10 (2024) 020

多少能量推动了流体

效率κ\kappa是释放的真空能里，最终变成单个泡泡流动的整体动能的那一份。它是w v2γ2ξ2w\,v^2\gamma^2\xi^2在自相似剖面上的积分，ww是焓密度，vv是流速，γ\gamma是洛伦兹因子。Espinosa、Konstandin、No 和 Servant 在 2010 年解出了这些剖面，并给出拟合公式，说明拟合与解的偏差小于 15%。

流体动能还要换算成占总能量密度的比例。2024 年的 LISA 模板再乘上 0.6，这个系数来自多泡泡模拟：泡泡多了以后，整体运动比孤立泡泡要小。

算个例子。α=0.2\alpha=0.2、ξw=0.9\xi_w=0.9时，拟合给出κ=0.2742\kappa=0.2742，我们自己解剖面得到0.27400.2740，于是K=0.6×0.2742×0.2/1.2=0.0274K=0.6\times0.2742\times0.2/1.2=0.0274。慢泡壁算例α=0.1\alpha=0.1、ξw=0.4\xi_w=0.4，拟合给0.18770.1877，我们的解是0.18900.1890，K=0.0102K=0.0102。

K≃0.6 κ α1+αK\simeq0.6\,\frac{\kappa\,\alpha}{1+\alpha}
0.6 是 arXiv:2403.03723v2 采用的多泡泡修正。2020 年的 LISA 综述 arXiv:1910.13125v2 写的是K=κα/(1+α)K=\kappa\alpha/(1+\alpha)，没有这个系数；同样的α\alpha和ξw\xi_w，两篇文章给出的KK不一样。

对应阅读：Espinosa, Konstandin, No & Servant, Energy budget of cosmological first-order phase transitions, arXiv:1004.4187v1 / JCAP 1006:028 (2010) · Caprini et al., Gravitational waves from first-order phase transitions in LISA: reconstruction pipeline and physics interpretation, arXiv:2403.03723v2 / JCAP 10 (2024) 020

我们的剖面解，以及它在哪里停下

我们自己积分了自相似流体方程，在泡壁和激波处用衔接条件，覆盖从 0.01 到 0.5 的六个强度、从 0.2 到 0.95 的泡壁速度。离 Jouguet 速度不太近时，和已发表拟合的最大相对偏差是 7.4%，出现在α=0.01\alpha=0.01、ξw=0.2\xi_w=0.2。

有两类点没有剖面。离 Jouguet 速度 0.01 以内，我们的打靶法不收敛，这些点就空着，不用拟合值去填。泡壁慢、相变又强时，根本不存在爆燃解：ξw=0.4\xi_w=0.4时有解的最大强度是α=0.590\alpha=0.590，ξw=0.2\xi_w=0.2时是0.3960.396。超过这个界，实验室不画谱。

声壳厚度也来自剖面。模板用ξshell=∣ξw−cs∣\xi_{\rm shell}=|\xi_w-c_s|，这对亚声速爆燃和爆轰成立。混合型的声壳更薄，模板没有为它标定，所以实验室把从cs−0.05c_s-0.05到ξJ\xi_J这一段标为超出适用范围。

三种相变强度下，单个泡泡的动能效率随泡壁速度的变化
α=0.03\alpha = 0.03、0.10.1 和 0.30.3 时的效率 κ\kappa。曲线是 Espinosa 等人 2010 年发表的拟合，圆点和方块是我们在袋模型里自己解出的剖面。虚线是声速 cs=1/3c_s = 1/\sqrt{3}。靠近 Jouguet 速度、我们的打靶法不收敛的地方没有画点。 查看完整图像

对应阅读：Original benchmark tutorial materials (standard-library Python) · Espinosa, Konstandin, No & Servant, Energy budget of cosmological first-order phase transitions, arXiv:1004.4187v1 / JCAP 1006:028 (2010)

练习与解释性答案

  1. α=0.05\alpha=0.05时，速度ξw=0.7\xi_w=0.7的泡壁是爆轰吗？α=0.01\alpha=0.01呢？

    展开解释

    算ξJ\xi_J：α=0.05\alpha=0.05时是0.7300.730，α=0.01\alpha=0.01时是0.6530.653。0.7 的泡壁在前一种情况是混合型，在后一种情况是爆轰。同样的速度换了类型，是因为相变越强，泡壁要跑得越快才能不扰动前方的等离子体。

  2. 用κ=0.1545\kappa=0.1545算α=0.1\alpha=0.1、ξw=0.9\xi_w=0.9时的KK，和同一强度的慢泡壁比一比。

    展开解释

    K=0.6×0.1545×0.1/1.1=0.00843K=0.6\times0.1545\times0.1/1.1=0.00843，比ξw=0.4\xi_w=0.4时的0.01020.0102略低。快爆轰并不自动更高效：很多能量用来加热泡壁后面的等离子体。效率在 Jouguet 速度附近最高。

  3. 为什么实验室不能直接对α=0.5\alpha=0.5、ξw=0.25\xi_w=0.25用拟合公式？

    展开解释

    拟合公式照样给出一个数，可袋模型里那里不存在稳定的爆燃解；我们的求解器给出ξw=0.25\xi_w=0.25时的界是α=0.429\alpha=0.429。没有物理流动支撑的数字会误导人，实验室就给出提示、不画谱。

动手检验预测

把泡壁速度扫过三种类型

  1. 先预测α=0.2\alpha=0.2时效率曲线的峰在哪里，再把泡壁速度从 0.4 慢慢拖到 0.95。
  2. 看画面标签从“激波跑在泡壁前面”变成“稀疏波拖在泡壁后面”，中间有一段不画谱。
  3. 把泡壁速度放在 0.3，再把α\alpha调到 0.5，读一读出现的提示。

来源与计算范围

  1. Caprini et al., Gravitational waves from first-order phase transitions in LISA: reconstruction pipeline and physics interpretation, arXiv:2403.03723v2 / JCAP 10 (2024) 020

    读的是 arXiv v2 的 LaTeX 源文件：摘要与引言、第 2.2.1 节（声波双折幂律、振幅、平均泡间距与红移因子）、第 3 节和附录 A（双参数 LISA 噪声与响应），以及第 5 节开头关于参数重建的部分。噪声公式把光学测量拐点印成 2e-3 mHz；我们的资料用 2 mHz，并同时给出两种结果。图 2 只读了图注。

  2. Caprini et al., Detecting gravitational waves from cosmological phase transitions with LISA: an update, arXiv:1910.13125v2 / JCAP (2020), doi:10.1088/1475-7516/2020/03/024

    只读了 LaTeX 源的两处：声波信号一节（signal.tex 约第 180 至 260 行），以及均方根流体速度和动能占比的定义（prelim.tex 第 171 至 182 行）。只用来对比早期的激波形成时间定义。

  3. Espinosa, Konstandin, No & Servant, Energy budget of cosmological first-order phase transitions, arXiv:1004.4187v1 / JCAP 1006:028 (2010)

    读了单个自相似泡泡的流体方程、爆燃、混合与爆轰三类剖面，以及附录 B 的效率拟合和其“优于 15%”的精度说明。我们在袋模型里自己解出的剖面独立检验这些拟合。

  4. Original benchmark tutorial materials (standard-library Python)

    我们的脚本、三个固定基准、四个失败案例和预期输出。它计算已发表的参数化模板、效率拟合和噪声模型；不求解爱因斯坦方程或流体方程，也不含浏览器引擎代码、论文原文或 LISA 数据。

继续阅读

这一页直接连着哪些内容

当前页面泡壁怎样推动等离子体

实验

文章

宇宙相变引力波 · 阅读路径

从泡泡成核到 LISA 可能看到的谱

  1. 文章宇宙早期若冒过泡泡，会留下怎样的引力波
  2. 文章泡泡怎样冒出来，又怎样填满空间
  3. 当前页面泡壁怎样推动等离子体
  4. 文章从声波寿命算到今天的引力波谱
  5. 文章独立复算一份基准声波谱
  6. 文章读懂 LISA 声波模板和它的边界
这个页面有错误、看不明白，或者缺了什么？告诉我们
反馈类型

提交时会附带本页地址、Lab ID、网站语言和设备类型。请不要填写个人信息。反馈最多保存 12 个月。隐私政策