基于NANOGrav 15年数据的早期宇宙随机引力波背景对有效相对论粒子数ΔNeff的约束研究

《COMPUTERS and EDUCATION》：Disentangling the origins of the NANOGrav signal: Early Universe models and Δ N eff bounds

本文针对NANOGrav合作组发现的纳赫兹随机引力波背景信号，系统研究了早期宇宙多种物理机制（包括暴胀、相变、宇宙弦等）产生的原初引力波对有效相对论粒子数ΔNeff的贡献。通过贝叶斯分析发现，运行谱指数模型在CMB-S4精度下可产生16σ显著度的ΔNeff信号，为利用未来CMB实验区分不同早期宇宙模型提供了新途径。该研究发表于《COMPUTERS and EDUCATION》，为理解早期宇宙物理开辟了多信使探测新窗口。

　　

近年来，脉冲星计时阵列（PTA）探测技术取得了突破性进展。2023年，北美纳赫兹引力波天文台（NANOGrav）合作组基于15年的观测数据，首次发现了与随机引力波背景（SGWB）特征一致的Hellings-Downs角关联信号，这一发现开启了纳赫兹引力波天文学的新纪元。更令人振奋的是，中国脉冲星计时阵列（CPTA）、欧洲脉冲星计时阵列（EPTA）和帕克斯脉冲星计时阵列（PPTA）也相继独立报道了类似信号，共同指向可能存在一个随机引力波背景。如果这一信号确实源自宇宙学起源，它将为我们打开一扇探索早期宇宙物理的独特窗口。

然而，一个关键问题随之而来：这个被探测到的随机引力波背景究竟来自何种早期宇宙物理过程？是暴胀期间产生的原初引力波，还是早期宇宙相变过程中产生的引力波，或者是宇宙弦等拓扑缺陷的遗迹？不同的起源机制预测的引力波能谱形状各不相同，从而对宇宙的能量成分产生不同的贡献。特别是，这些引力波会对宇宙的有效相对论粒子数ΔN_eff产生额外贡献，而ΔN_eff是描述早期宇宙中相对论性粒子种类的关键物理量。

为了回答这一问题，研究人员在《COMPUTERS and EDUCATION》上发表了最新研究成果，系统研究了六种不同的早期宇宙模型：标准暴胀模型、具有运行谱指数的暴胀模型、再加热模型、标量诱导引力波（SIGW）、相变产生的引力波（包括气泡碰撞和声波激发两种机制）以及亚稳态宇宙弦。研究团队采用贝叶斯分析方法，利用NANOGrav 15年数据对各个模型的参数空间进行了全面探索，并计算了每种模型预测的随机引力波背景对ΔN_eff的贡献。

研究的关键创新点在于，它不仅关注随机引力波背景本身的性质，还深入分析了这些引力波对宇宙微波背景（CMB）观测量的影响，特别是对ΔN_eff的贡献。通过将脉冲星计时阵列的探测与宇宙微波背景的观测相结合，研究为区分不同的早期宇宙模型提供了新的判别标准。随着西蒙斯天文台（SO）和CMB-S4等下一代宇宙微波背景实验的推进，对ΔN_eff的测量精度将大幅提高，这将使我们可以通过测量ΔN_eff来间接限制早期宇宙的物理过程。

研究团队采用了几项关键技术方法：首先，他们使用贝T3分析框架处理NANOGrav 15年数据；其次，利用PTArcade软件包中的"ceffyl"模式进行随机引力波背景能谱的建模；第三，采用贝叶斯推断方法对各个模型的参数空间进行采样；最后，基于getdist工具进行后验分布的分析和可视化。特别值得注意的是，研究考虑了前14个频率区间的数据，以避开脉冲星本身的高频噪声。

研究结果揭示了不同早期宇宙模型对ΔN_eff贡献的显著差异：

标准暴胀模型的分析显示，张量谱指数n_t = 1.79^+0.33_-0.38，张标比log₁₀r = -8.3^+3.0_-2.6，对应的ΔN_eff贡献为log₁₀ΔN_eff = -2.72^+0.27_-0.32。在当前CMB实验精度下，这一贡献的显著性约为1.5σ，而未来西蒙斯天文台（SO）的观测有望将显著性提升至6.4σ。

运行谱指数模型表现出最强的探测前景，其参数约束为n_t = 2.48^+1.20_-0.56，log₁₀r = -14.4^+4.7_-11.0，运行系数α_t = -0.09^+0.10_-0.12。该模型预测的ΔN_eff贡献在当前精度下就达到4σ显著性，CMB-S4时代更是可达16σ，显示出极强的可探测性。

再加热模型的研究发现，再加热温度log₁₀T_re = 3.4^+3.1_-8.1 GeV，状态方程参数w_re = 0.43^+0.17_-0.29，对应的ΔN_eff贡献在当前限制下显著性较低（0.8σ），但SO实验有望将其提升至3.5σ。

相变模型的研究分为气泡碰撞和声波激发两种机制。气泡碰撞模型的log₁₀α_?? = 0.29^+0.65_-0.29，log₁₀T_?? = -0.71^+0.36_-0.47 GeV，log₁₀H_??R_?? = -0.36^+0.32_-0.16；声波激发模型对应的参数为log₁₀α_?? = 0.24^+0.46_-0.56，log₁₀T_?? = -1.71^+0.38_-0.43 GeV，log₁₀H_??R_?? = -0.71^+0.26_-0.43。两种机制对ΔN_eff的贡献均较为微弱，探测前景相对有限。

标量诱导引力波（SIGW）模型的分析给出了严格的参数约束：振幅log₁₀A_s = -2.42^+0.15_-0.17，谱指数n_s = 1.47^+0.38_-0.45，运行系数α = 0.0281^+0.0087_-0.0011，运行之运行β = 0.00995^+0.00076_-0.00088。该模型预测的ΔN_eff贡献极为微小，log₁₀ΔN_eff = -3.894^+0.031_-0.035，在当前和未来的实验精度下都难以探测。

宇宙弦模型的研究显示，弦张力参数log₁₀(Gμ) = -5.44^+0.76_-0.39，单极子与弦张力比值的平方根√κ = 7.887^+0.039_-0.11。该模型对ΔN_eff的贡献为log₁₀ΔN_eff = -1.78^+0.098_-0.210，在当前实验精度下显著性约为1.6σ，SO实验有望提升至6.5σ。

研究的结论部分强调了多信使天文学在理解早期宇宙物理中的重要性。通过将脉冲星计时阵列的引力波探测与宇宙微波背景的精确测量相结合，我们可以对早期宇宙的不同物理过程进行有效区分。特别值得注意的是，运行谱指数模型显示出极强的探测前景，这为未来的观测提供了明确的目标。

此外，研究还指出，不同早期宇宙模型对ΔN_eff的贡献存在数量级上的差异，这为利用下一代CMB实验区分这些模型提供了可能。例如，运行谱指数模型预测的ΔN_eff贡献比标准暴胀模型大一个数量级以上，这种差异在未来高精度实验中将是明显可区分的。

这项研究的重要意义在于，它首次系统地将NANOGrav发现的随机引力波背景信号与早期宇宙物理联系起来，并定量评估了不同模型对ΔN_eff的贡献。这不仅为理解当前观测数据提供了理论框架，也为未来的观测实验指明了方向。随着更多脉冲星计时阵列数据的积累和下一代CMB实验的开展，我们有望最终揭示纳赫兹引力波背景的起源，从而深化对早期宇宙物理的理解。

研究的讨论部分进一步指出，尽管当前数据还无法明确区分不同的早期宇宙模型，但已经排除了某些参数空间。例如，标量诱导引力波模型对ΔN_eff的贡献过小，即使在未来实验中也难以探测；而运行谱指数模型则显示出极强的探测前景，这提示我们应该在未来实验中重点关注这类模型。

最后，研究强调了持续监测脉冲星计时阵列信号的重要性。随着观测时间的延长和更多脉冲星的加入，PTA的灵敏度将进一步提高，有望在更宽的频率范围内刻画引力波能谱的形状。结合CMB对ΔN_eff的精确测量，我们将能够构建早期宇宙的更完整图像，从而回答宇宙极早期的一些基本物理问题。