[论文解读] Footprints of the QCD Crossover on Cosmological Gravitational Waves at Pulsar Timing Arrays

Gabriele Franciolini, Davide Racco+1 more|arXiv (Cornell University)|Jun 29, 2023

Cosmology and Gravitation Theories被引用 4

一句话总结

本文表明，早期宇宙中的QCD相变在脉冲星计时阵列（PTA）可探测的宇宙引力波背景低频尾部上留下了独特的光谱特征。通过将因果性限制的引力波谱建模为QCD引起的状态方程变化，作者发现该模型在贝叶斯对比中显著优于标准幂律拟合，表明标准模型物理可对nHz引力波信号的解释产生决定性影响。

ABSTRACT

Pulsar Timing Arrays (PTAs) have reported evidence for a stochastic gravitational wave (GW) background at nHz frequencies, possibly originating in the early Universe. We show that the spectral shape of the low-frequency (causality) tail of GW signals sourced at temperatures around $T\gtrsim 1$ GeV is distinctively affected by confinement of strong interactions (QCD), due to the corresponding sharp decrease in the number of relativistic species. Bayesian analyses in the NANOGrav 15 years and the previous International PTA datasets reveal a significant improvement in the fit with respect to cubic power-law spectra, previously employed for the causality tail. This suggests that the inclusion of Standard Model effects on GWs can have a potentially decisive impact on model selection.

研究动机与目标

• 研究标准模型物理（尤其是QCD相变）是否会在原初引力波背景的低频尾部留下可探测的印记。
• 通过引入QCD引起的态方程和相对论自由度的变化，改进PTA中因果性限制引力波信号的建模。
• 评估在标准幂律模型基础上引入QCD效应后，是否能显著提升对当前PTA数据的拟合效果。
• 提供一种与模型无关的框架，用于根据早期宇宙动力学和标准模型物理解释nHz引力波信号。

提出的方法

• 作者推导了在时间尺度与哈勃时间相当的早期宇宙源中产生的引力波因果性尾部（CT）的光谱形状，考虑了QCD相变对态方程的影响。
• 他们在低频区域对张量模式 $ h_k(t) $ 的演化进行建模，此时引力波对源细节不敏感，仅依赖于膨胀历史和引力波传播特性。
• 利用格点QCD计算的相对论自由度数量，他们计算了红移后的能量密度 $ \Omega_{\textsc{gw}}(f) $，表明 $ f^3 $ 幂律倾斜因QCD相变而发生改变。
• 他们使用IPTA DR2和NANOGrav 15年数据集，通过在enterprise和enterprise_extensions管道中实现PTMCMC采样，进行贝叶斯模型比较。
• 分析包括噪声参数（白噪声、红噪声、色散测量）和引力波背景参数，其先验分布基于宇宙学约束和数据驱动限制。
• 他们将新提出的QCD修正CT模型与标准幂律和分段幂律模型进行比较，通过计算贝叶斯因子评估模型偏好。

实验结果

研究问题

• RQ1 早期宇宙中的QCD相变是否会在宇宙引力波背景的低频尾部留下可探测的印记？
• RQ2 在态方程中引入QCD引起的改变，是否能提升对脉冲星计时阵列引力波数据的拟合效果，相比标准幂律模型？
• RQ3 标准模型物理对nHz随机引力波背景模型选择的定量影响是什么？
• RQ4 由于QCD相变期间相对论自由度的快速下降，因果性限制引力波信号的光谱特征如何变化？

主要发现

• 在IPTA DR2和NANOGrav 15年数据集中，将QCD相变效应纳入因果性尾部模型后，贝叶斯拟合显著优于标准立方幂律拟合。
• 对于IPTA DR2，QCD修正CT信号的后验中位数振幅为 $ \log_{10}A_{\mathrm{CT}} = -13.67^{+0.08}_{-0.07} $，光谱指数与预期的 $ f^3 $ 行为一致。
• 对于NANOGrav 15年数据，后验中位数振幅为 $ \log_{10}A_{\mathrm{CT}} = -13.85^{+0.06}_{-0.05} $，表明信号相对于噪声具有强烈偏好。
• 贝叶斯证据支持QCD修正CT模型优于标准幂律模型，表明标准模型效应可在解释当前PTA数据中发挥决定性作用。
• 分析表明，受QCD相变影响的物质主导到辐射主导的过渡，会改变高频和低频光谱指数，分别为 $ \gamma_{\text{MD}} = 4 $ 和 $ \delta_{\text{MD}} = 2 $。
• 结果表明，未来模型比较必须包含QCD效应，以避免在区分nHz引力波信号的宇宙学与天体物理起源时产生偏差。