[论文解读] Implications of Pulsar Timing Array Results for High Frequency Gravitational Waves

Junsong Cang, Yu Gao+2 more|arXiv (Cornell University)|Sep 26, 2023

Pulsars and Gravitational Waves Research被引用 7

一句话总结

本文研究了近期脉泽定时阵列（PTA）对纳赫兹引力波背景的探测结果，提出这些信号可能源自原初曲率扰动产生的标量诱导引力波（SIGW）。若属实，所需扰动振幅将产生质量范围为 $[2\times 10^{-5}, 2\times 10^{-2}] olimits\,M_\odot$ 的原初黑洞（PBH），其并合将产生可被未来观测台如爱因斯坦望远镜和DECIGO探测的 $[10^{-3}, 10^8]$ Hz 频段的高频引力波背景。

ABSTRACT

Several pulsar timing array (PTA) experiments such as NANOGrav and PPTA recently reported evidence of a gravitational wave (GW) background at nano-Hz frequency band. This signal can originate from scalar-induced gravitational waves (SIGW) generated by the enhanced curvature perturbation. Production of SIGW is expected to be accompanied by formation of primordial black holes (PBH), which can emit GW through binary mergers. Here we perform a joint likelihood inference on PTA datasets in combination with existing limits on PBH abundance and GW density, we derive full Bayesian posteriors for PBH distribution and relevant PBH merger signal. Our results show that analysis using PTA data alone implies significant overproduction of PBHs, and accounting for current PBH limits causes visible shifts in SIGW posterior. If PTA signals are indeed of SIGW origin, the required curvature perturbation amplitude produces PBHs in a narrow mass window of $[6 imes 10^{-2}, 2 imes 10^{-1}]\ m_\odot$. Mergers of these PBHs can leave a strong GW signature in $[10^{-3}, 10^5]$ Hz frequency range, to be detectable at upcoming interferometers such as LISA, aLIGO, Einstein Telescope, DECIGO and BBO, etc. This offers a multi-frequency opportunity to further scrutinize the source of the observed PTA signal and can potentially improve current PBH constraints by up to 5 orders of magnitudes.

研究动机与目标

• 评估最近在脉冲星定时阵列（PTA）中探测到的纳赫兹引力波背景是否可能源自由增强原初曲率扰动产生的标量诱导引力波（SIGW）。
• 若PTA信号确实源于SIGW，利用PTA数据与宇宙学约束的联合似然推断，确定由此产生的原初黑洞（PBH）质量分布与丰度。
• 评估未来激光干涉仪探测器（如爱因斯坦望远镜、DECIGO和BBO）在 $[10^{-3}, 10^8]$ Hz 频段探测PBH并合产生的高频引力波信号的可探测性。
• 建立低频PTA观测与高频引力波搜寻之间的多信使联系，实现对SIGW起源假说的交叉验证。

提出的方法

• 采用对数正态曲率功率谱 $P_{\mathcal{R}}(k) = \frac{A}{\sqrt{2\pi\Delta^2}} \exp\left[-\frac{(\ln k/k_*)^2}{2\Delta^2}\right]$ 建模具有振幅 $A$、峰值波数 $k_*$ 与宽度 $\Delta$ 的原初曲率扰动。
• 应用二级扰动理论计算标量诱导引力波（SIGW）的能量密度，使用公式 $\Omega_{\rm GW} = 0.29\Omega_{\rm r}\left(\frac{106.75}{g_*}\right)^{1/3} \int\!\!\!\int \cdots F(u,v)P_{\mathcal{R}}(kv)P_{\mathcal{R}}(ku)\,du\,dv$。
• 通过 $\psi = \frac{df_{\rm bh}}{d\ln m} = \frac{f_{\rm bh}}{\sqrt{2\pi}\sigma_{\rm bh}} \exp\left[-\frac{\ln^2(m/m_{\rm c})}{2\sigma_{\rm bh}^2}\right]$ 推导PBH质量分布，参数化为峰值质量 $m_{\rm c}$ 与宽度 $\sigma_{\rm bh}$。
• 利用 $\beta(m) \propto \sqrt{\frac{2\bar{\sigma}^2}{\pi\delta_{\rm c}^2}} \exp\left(-\frac{\delta_{\rm c}^2}{2\bar{\sigma}^2}\right)$ 计算PBH形成概率，其中 $\bar{\sigma}^2(m)$ 由曲率功率谱与高斯窗函数 $W(x) = \exp(-x^2/2)$ 推导。
• 对PTA数据集（NANOGrav、PPTA、IPTA）与 $\Delta N_{\rm eff}$ 及PBH丰度约束进行联合似然推断，映射SIGW与PBH参数的后验分布。
• 估算PBH双星的共动并合率 $R$ 与由此产生的高频引力波背景，利用未来探测器的灵敏度曲线评估 $[10^{-3}, 10^8]$ Hz 频段的可探测性。

实验结果

研究问题

• RQ1 若PTA探测到的纳赫兹引力波背景源自标量诱导引力波（SIGW），则需多大的原初曲率扰动振幅？
• RQ2 此类曲率扰动谱将产生何种质量范围的原初黑洞（PBH），其占暗物质的比例是多少？
• RQ3 这些PBH并合将产生何种高频引力波背景，其峰值频率位于哪个频段？
• RQ4 由PBH并合预测的高频引力波信号是否可被未来激光干涉仪观测台（如爱因斯坦望远镜、DECIGO与BBO）探测到？
• RQ5 如何利用从纳赫兹（PTA）到MHz（未来干涉仪）的多频段引力波观测，测试PTA信号的SIGW起源？

主要发现

• 为解释PTA观测到的纳赫兹引力波背景，所需原初曲率扰动振幅将产生质量范围为 $[2 \times 10^{-5}, 2 \times 10^{-2}]\,M_\odot$ 的原初黑洞（PBH）。
• PBH质量分布可用对数正态分布良好描述，其峰值质量 $m_{\rm c}$ 与宽度 $\sigma_{\rm bh}$ 与PTA数据反演的 $P_{\mathcal{R}}$ 参数一致。
• 这些PBH的并合将产生峰值位于 $[10^{-3}, 10^8]$ Hz 频段的引力波背景，其幅度在未来的探测器敏感频段内可被探测。
• 预测的高频引力波信号位于未来探测器（如爱因斯坦望远镜、DECIGO与大爆炸观测者BBO）的灵敏度范围内，具备潜在交叉验证能力。
• 对PTA数据、$\Delta N_{\rm eff}$ 约束与PBH丰度限制的联合推断表明，SIGW情景的允许参数空间与现有天体物理与宇宙学约束一致。
• 该情景提供了一条多信使路径：单一机制（增强的曲率扰动）可同时解释低频PTA信号与可探测的高频引力波背景。