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

读懂 LISA 声波模板和它的边界

Sandbox Physics ·

打开 arXiv:2403.03723v2，“Gravitational waves from first-order phase transitions in LISA: reconstruction pipeline and physics interpretation”（Caprini 等，已发表于 JCAP 10 (2024) 020）。这篇导读带你读它关于声波的 2.2.1 节、参数依赖图里的声波分图，以及第 5 节热力学参数的重建。另外两篇讨论声波源寿命的论文，arXiv:2003.07360v4 和 arXiv:2007.08537v3，会告诉你模板里哪些地方是一种选择，而不是一个结果。

每篇来源说的是哪一类话

本课程用到的五篇论文都已在期刊发表：arXiv:2403.03723 在 JCAP 2024，arXiv:2003.07360 在 JCAP 2020，arXiv:2007.08537 在 JCAP 2021，arXiv:1910.13125 在 JCAP 2020，arXiv:1004.4187 在 JCAP 2010。我们读的是这里注明版本号的 arXiv 稿，不是期刊排版稿。它们都没有报告观测：相变引力波背景至今没有被探测到。

内容分三类。对流体模拟的拟合，比如谱的斜率和常数AswA_{\rm sw}，是模拟范围内的数值结果。模板、寿命规则和 LISA 预测是模型预言。接近 Jouguet 速度时的谱、激波形成之后的谱，论文提出了却没有解决，属于未解决问题。我们自己的蒙特卡罗记录和剖面解，是为这门课做的合成计算。

对应阅读：Caprini et al., Gravitational waves from first-order phase transitions in LISA: reconstruction pipeline and physics interpretation, arXiv:2403.03723v2 / JCAP 10 (2024) 020 · Ellis, Lewicki & No, Gravitational waves from first-order cosmological phase transitions: lifetime of the sound wave source, arXiv:2003.07360v4 / JCAP 07 (2020) 050 · Guo, Sinha, Vagie & White, Phase transitions in an expanding universe: stochastic gravitational waves in standard and non-standard histories, arXiv:2007.08537v3 / JCAP 01 (2021) 001 · 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

2.2.1 节：先认清每个符号，再看公式

这一节先讲物理图像：流体的整体运动起初以压缩为主，泡泡相撞以后，压缩波成为引力波的源。接着引用了格点模拟和“无希格斯场”方法；后者不再求解标量场，而把不断扩张的真真空当作给定的背景，从而能精确解出流体和其中的激波。

sw_shape、sw_amplitude 和 Rstar 三个式子定义了模板。仔细读Δw=ξshell/max⁡(ξw,cs)\Delta_w=\xi_{\rm shell}/\max(\xi_w,c_s)：文中说ξshell=∣ξw−cs∣\xi_{\rm shell}=|\xi_w-c_s|对亚声速爆燃和爆轰成立，并说明这是“本文关注的情形，因为我们取ξw\xi_w接近 1”。把混合型泡壁代进同一个式子，就超出了作者说明的范围。

动能占比是K≃0.6 κα/(1+α)K\simeq0.6\,\kappa\alpha/(1+\alpha)。0.6 来自多泡泡模拟，文中称它比 2020 年综述的取法K=κα/(1+α)K=\kappa\alpha/(1+\alpha)更保守。常见误读：同样的α\alpha和ξw\xi_w下直接比较两篇文章的振幅，没注意K2K^2差了 0.36 倍。

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

参数依赖图里的声波分图

在源文件里找到标签为 fig:pars_PT_sw 的分图，它是参数依赖图的第二个分图。图注把黑色参考谱定在K=0.1K=0.1、H∗R∗=0.1H_*R_*=0.1、ξw=0.9\xi_w=0.9、T∗=1 TeVT_*=1\,\mathrm{TeV}；每组彩色曲线只改变一个参数，灰色虚线是 LISA 噪声。

我们读的是图注和公式，没有逐点读图上的曲线。从公式可以先预测每组曲线该怎么变，再去图里核对。改T∗T_*，曲线左右平移，高度不变。改ξw\xi_w，只有第二个拐点通过Δw\Delta_w移动。改H∗R∗H_*R_*，两个拐点都动，高度也变，因为H∗R∗H_*R_*和寿命都进了振幅。改KK，寿命停在 1 时曲线按K2K^2升高；激波先出现时升得慢一些。

实验室画的是同一个模板，你可以自己做这些单参数改动；模板不适用的地方，金色曲线就不出现。

对应阅读：Caprini et al., Gravitational waves from first-order phase transitions in LISA: reconstruction pipeline and physics interpretation, arXiv:2403.03723v2 / JCAP 10 (2024) 020 · Original benchmark tutorial materials (standard-library Python)

声波能响多久：三种写法

Ellis、Lewicki 和 No（arXiv:2003.07360v4）把激波形成时间估为τsh∼Lf/Uˉf\tau_{\rm sh}\sim L_f/\bar U_f，其中Uˉf2≃34κα/(1+α)\bar U_f^2\simeq\tfrac34\kappa\alpha/(1+\alpha)。他们在摘要里指出，没有零温势垒的模型和标度不变模型里，声波只响不到一个哈勃时间。2020 年的 LISA 综述把功率乘上H∗τsh=H∗R∗/K1/2H_*\tau_{\rm sh}=H_*R_*/K^{1/2}。2024 年的模板用H∗R∗/K/ΓH_*R_*/\sqrt{K/\Gamma}，并以 1 为上限。

对我们的基准算例α=0.2\alpha=0.2、β/H∗=50\beta/H_*=50、ξw=0.9\xi_w=0.9，三种写法分别给出0.2850.285、0.2470.247和0.3680.368。物理图像一样，只是流速的定义不同。Guo、Sinha、Vagie 和 White（arXiv:2007.08537v3）在辐射主导下用Υ=1−1/1+2H∗τsh\Upsilon=1-1/\sqrt{1+2H_*\tau_{\rm sh}}代替硬性上限，得0.2410.241。

激波形成以后引力波怎样继续产生，这几篇论文都没有定论。2024 年的模板此后不再把它算作声波，而是另加一项不确定的湍流贡献，配上它自己的效率参数。寿命因子因此是振幅里最大的建模不确定性之一。

两种声源寿命因子随激波形成时间的变化，标出三个基准
金线是模板用的因子 min⁡(H∗τsh,1)\min(H_*\tau_{\mathrm{sh}},1)，紫线是 Guo 等人 2021 年的膨胀压低因子 Υ=1−1/1+2H∗τsh\Upsilon = 1 - 1/\sqrt{1+2H_*\tau_{\mathrm{sh}}}。实心点和空心点把三个基准分别放到两条曲线上。H∗τsh≪1H_*\tau_{\mathrm{sh}} \ll 1 时两者一致；差别在 H∗τsh=1H_*\tau_{\mathrm{sh}} = 1 处最大，模板因子约为 Υ\Upsilon 的 2.42.4 倍。 查看完整图像

对应阅读：Ellis, Lewicki & No, Gravitational waves from first-order cosmological phase transitions: lifetime of the sound wave source, arXiv:2003.07360v4 / JCAP 07 (2020) 050 · Guo, Sinha, Vagie & White, Phase transitions in an expanding universe: stochastic gravitational waves in standard and non-standard histories, arXiv:2007.08537v3 / JCAP 01 (2021) 001 · 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 · Caprini et al., Gravitational waves from first-order phase transitions in LISA: reconstruction pipeline and physics interpretation, arXiv:2403.03723v2 / JCAP 10 (2024) 020

泡壁速度是输入，不是输出

这里所有公式里，ξw\xi_w都是自由参数。我们读的论文没有一篇从泡壁和等离子体之间的摩擦把它算出来。2020 年的综述说，它的拟合式在泡壁速度远离 Chapman–Jouguet 速度时最可靠：爆燃要明显低于csc_s，爆轰要明显高于它。接近这个速度时声壳很薄，厚度约为R∗∣ξw−cs∣/csR_*|\xi_w-c_s|/c_s，谱形公式失效。

这关系到能不能分辨参数。如果一次真实相变的泡壁速度接近声速，模板会把第二个拐点算错，而携带ξw\xi_w信息的正是这个拐点。

对应阅读：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 · Caprini et al., Gravitational waves from first-order phase transitions in LISA: reconstruction pipeline and physics interpretation, arXiv:2403.03723v2 / JCAP 10 (2024) 020

第 5 节：LISA 能分辨出什么

arXiv:2403.03723v2 第 5 节点明了关键的计数问题：双折幂律只有三个几何参数，即振幅和两个拐点；声波却依赖四个热力学参数KK、H∗R∗H_*R_*、ξw\xi_w和T∗T_*，反推必然有简并。作者表明，同时重建不止一种贡献时，简并可以解开；他们的做法是先重建几何参数，再换算成热力学参数。

引用“精度优于约 10%”这个预测之前，先读摘要和引言里写明的前提。模拟数据和重建用的是同一套信号、噪声和前景模板，模板本身的误差按构造就不存在；泡泡大到与哈勃尺度相当时的广义相对论效应被忽略；噪声模型只有两个参数。预测成立的条件是“到 LISA 升空时，某些实验和理论上的不确定性已经解决”。

读这个领域下一篇论文时可以照着核对：模板覆盖哪些泡壁类型；KK和流速用的是哪种定义；寿命怎样封顶；有没有包括前景；注入数据和拟合用的是不是同一个模板；哪些数字是模拟拟合，哪些是预测，哪些问题还没有答案。

对应阅读：Caprini et al., Gravitational waves from first-order phase transitions in LISA: reconstruction pipeline and physics interpretation, arXiv:2403.03723v2 / JCAP 10 (2024) 020 · Robson, Cornish & Liu, The construction and use of LISA sensitivity curves, arXiv:1803.01944v2 / Class. Quantum Grav. (2019), doi:10.1088/1361-6382/ab1101

练习与解释性答案

  1. 某篇论文对α=0.3\alpha=0.3、ξw=0.95\xi_w=0.95，用K=κα/(1+α)K=\kappa\alpha/(1+\alpha)给出h2Ωpeak=10−10h^2\Omega_{\rm peak}=10^{-10}。同样的输入和寿命下，2024 年的模板大约给出多少？

    展开解释

    2024 年的模板把KK乘以 0.6，K2K^2降为 0.36 倍。如果两者寿命都封顶为 1，谱峰降到约3.6×10−113.6\times10^{-11}。如果寿命由激波决定，H∗τsh∝K−1/2H_*\tau_{\rm sh}\propto K^{-1/2}会变长，部分抵消，因子约为0.36/0.6≈0.460.36/\sqrt{0.6}\approx0.46；不过这时流速定义的差别也要算进去。先对齐定义，再比较数字。

  2. 假设 LISA 把f1f_1、f2f_2和Ω2\Omega_2测得完全准确。没有额外信息时，哪个热力学组合仍然分不开？

    展开解释

    两个拐点都正比于H∗,0/(H∗R∗)∝T∗/(H∗R∗)H_{*,0}/(H_*R_*)\propto T_*/(H_*R_*)，所以f1f_1只定出比值T∗/(H∗R∗)T_*/(H_*R_*)。比值f2/f1f_2/f_1定出Δw\Delta_w，进而定出ξw\xi_w。振幅再定出K2(H∗τsw)(H∗R∗)K^2(H_*\tau_{\rm sw})(H_*R_*)。还缺一个关系：需要粒子物理模型给出T∗T_*，或者测到湍流这样的第二种源。

动手检验预测

一次只改一个参数，看模板怎样变

  1. 用默认预设，只把ξw\xi_w从 0.85 移到 0.95，再只改β/H∗\beta/H_*。实验室的滑块是α\alpha、β/H∗\beta/H_*和ξw\xi_w，所以改泡壁速度也会改变R∗R_*和κ\kappa；即便如此，第二个拐点仍比第一个移得多。
  2. 设α=0.5\alpha=0.5、ξw=0.95\xi_w=0.95，在寿命图里比较两种寿命规则。

基准复算教程资料

下载离线教程 ZIP · 阅读单位、版本与适用范围

只用 Python 标准库，为三个固定基准复算声源寿命、两个拐点频率、今天的模板谱和 LISA 信噪比，并演示四个失败案例怎样被拒绝。它计算的是已发表的参数化模板，不求解爱因斯坦方程或流体方程。

文件、字节与 SHA-256

来源与计算范围

  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. Ellis, Lewicki & No, Gravitational waves from first-order cosmological phase transitions: lifetime of the sound wave source, arXiv:2003.07360v4 / JCAP 07 (2020) 050

    读了 LaTeX 源中关于激波形成时间、用效率表示的平均流体速度，以及“强相变时声波源远早于一个哈勃时间就关闭”的论证。

  5. Guo, Sinha, Vagie & White, Phase transitions in an expanding universe: stochastic gravitational waves in standard and non-standard histories, arXiv:2007.08537v3 / JCAP 01 (2021) 001

    读了主 LaTeX 文件里辐射主导宇宙的膨胀压低因子。它和模板的寿命因子并排给出作比较，模板本身没有用它。

  6. Robson, Cornish & Liu, The construction and use of LISA sensitivity curves, arXiv:1803.01944v2 / Class. Quantum Grav. (2019), doi:10.1088/1361-6382/ab1101

    在 LaTeX 源中找到了光学测量噪声、检验质量噪声和灵敏度公式。只用来独立核对噪声曲线形状和 2 mHz 拐点，没有逐节通读。

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

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

继续阅读

这一页直接连着哪些内容

当前页面读懂 LISA 声波模板和它的边界

实验

文章

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

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

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

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