[论文解读] Primordial gravitational waves in the nano-Hertz regime and PTA data -- towards solving the GW inverse problem

Eric Madge, Enrico Morgante+5 more|arXiv (Cornell University)|Jun 26, 2023

Pulsars and Gravitational Waves Research参考文献 186被引用 6

一句话总结

本文提出了适用于纳赫兹频段原初引力波的简单、极简模型，涵盖相变、宇宙弦、轴子动力学及标量涨落，以解释近期脉冲星时标阵列（PTA）对随机引力波背景的暗示。通过结合PTA数据与宇宙学约束进行贝叶斯拟合，识别出可行的参数区域，并为未来的天体物理与实验室实验提供可检验的预测，推动了引力波反问题的解决。

ABSTRACT

In recent years, several pulsar timing array collaborations have reported first hints for a stochastic gravitational wave background at nano-Hertz frequencies. Here we elaborate on the possibility that this signal comes from new physics that leads to the generation of a primordial stochastic gravitational wave background. We propose a set of simple but concrete models that can serve as benchmarks for gravitational waves sourced by cosmological phase transitions, domain wall networks, cosmic strings, axion dynamics, or large scalar fluctuations. These models are then confronted with pulsar timing data and with cosmological constraints. With only a limited number of free parameters per model, we are able to identify viable regions of parameter space and also make predictions for future astrophysical and laboratory tests that can help with model identification and discrimination.

研究动机与目标

• 识别能够解释脉冲星时标阵列（PTA）在纳赫兹频段所暗示的随机引力波背景的简单、有物理动机的原初引力波模型。
• 将这些模型同时与PTA数据及宇宙学约束（包括CMB、原初核合成与中微子退耦）进行对比。
• 通过减少模型复杂度与参数数量，提升贝叶斯模型比较的可靠性，避免过拟合，特别是针对一阶相变模型。
• 通过未来天体物理与实验室测试的预测，实现对竞争性引力波源（如超大质量双黑洞与新物理）的区分。
• 通过识别与当前数据一致的模型以及可被未来数据释放排除或证实的模型，推进引力波反问题的解决。

提出的方法

• 构建了来自宇宙学一阶相变、域壁、宇宙弦、轴子动力学及标量涨落的引力波基准模型，每种模型均仅引入最少的新场与参数。
• 使用既有的理论框架计算引力波能量密度谱：相变采用气泡壁动力学，流体/声波效率因子，标量涨落则通过随机背景生成机制。
• 利用enterprise与ceffyl软件套件，对NANOGrav 12.5年与IPTA DR2数据集进行贝叶斯参数拟合，对log10(Tp/GeV)、log10(α)与log10(β/H*)采用均匀先验。
• 应用来自CMB、BBN与中微子退耦的宇宙学约束，以限制可行参数空间，尤其针对MeV量级的相变。
• 使用Mukhanov-Sasaki方程计算标量功率谱，并与近似慢滚结果对比，以验证峰值振幅与谱指数。
• 分析引力波信号的频率依赖性，评估模型区分潜力，并识别出既匹配PTA数据又满足宇宙学约束的区域。

实验结果

研究问题

• RQ1 是否能够通过简单、极简的原初引力波模型，在不过拟合且不违反宇宙学约束的前提下，解释当前的PTA数据？
• RQ2 哪些引力波源机制——相变、宇宙弦、轴子动力学或标量涨落——能够产生与观测到的nHz信号相容的随机背景？
• RQ3 PTA信号的频率谱是否已具备区分不同原初引力波源的能力？
• RQ4 在施加宇宙学约束后，各类模型的可行参数区域为何？不同PTA合作组织之间的差异如何？
• RQ5 未来哪些天体物理或实验室实验最能有效区分竞争性模型？

主要发现

• NANOGrav 12.5年数据倾向于MeV量级的相变温度，其引力波谱的高频部分与红噪声共同过程相匹配；而IPTA DR2则偏好GeV量级的渗滤温度。
• 对于一阶相变，存在可行的参数区域，其中α < 10，Tp/GeV ∈ [10−4, 102]，β/H* ∈ [1, 1000]，且真空与流体壳层贡献在适当的效率因子下均可拟合数据。
• 具有全局或规范对称性的轴子及轴子样粒子模型，可通过域壁或宇宙弦网络产生与PTA数据一致的引力波信号。
• 由暴胀涨落引起的标量诱导引力波仅在谱指数与振幅被精细调节时才能匹配CMB观测，但此类模型通常产生振幅过低或波段过小的峰值，难以在PTA频段探测到。
• 具有Ne ≈ 60埃夫拉德与PTA频段内峰值的基准点，其振幅与谱指数均过低；而匹配CMB数据的模型其峰值过小且位于远小于PTA探测波段的尺度上。
• NANOGrav与IPTA数据的组合（简单平均）与单一来源模型一致，但各数据集偏好不同的参数区域，凸显了联合分析与未来数据释放的必要性。