NANOGrav与纳赫兹引力波背景

概述#

NANOGrav（North American Nanohertz Observatory for Gravitational Waves，北美纳赫兹引力波天文台）是美国和加拿大天文学家利用 Arecibo 射电望远镜（2020 年前）和绿岸望远镜（GBT）长期监测毫秒脉冲星计时的合作计划，自 2004 年起持续运行。

它做的不是一次闪光式捕捉，而是把银河系里几十颗毫秒脉冲星当成分布式时钟，年复一年记下脉冲到达时间。仪器是射电望远镜，方法叫脉冲星计时阵列；NANOGrav 只是北美这一支，不是全球唯一的阵列。灵敏度来自时间基线有多长、脉冲星对有多少，不来自某一次「听见了波」。

2023 年 6 月 29 日，NANOGrav 发布了基于 15 年数据集（包含 68 颗毫秒脉冲星）的分析结果，以超过 $3\sigma$ 的显著性报告发现赫尔宾斯-道-里曼（Hellings-Downs）角相关性——引力波背景最关键的"指纹"（Agazie et al. 2023，The Astrophysical Journal Letters，951，L8）。同期，欧洲（EPTA）、澳大利亚（PPTA）和中国（CPTA）的三个独立 PTA 团队也发布了一致的结果，共同构成了纳赫兹引力波天文学的历史性时刻。

那篇关键论文的标题用的是 evidence（证据），不是 detection（探测确认）。Hellings-Downs 角相关被看见了，但同一份数据换一种统计，显著性仍会在大约三到四个标准差之间移动。把它写成终审，会把尚未走完的检验说成已经结案。

NANOGrav 的发展历史#

从 Arecibo 到 GBT#

NANOGrav 的核心观测工具是：

Arecibo 望远镜（305 米口径，波多黎各）：世界最大的单口径射电望远镜之一，直到 2020 年 12 月坍塌前，是 NANOGrav 灵敏度最高的观测平台。其极高灵敏度使得对毫秒脉冲星的每次脉冲计时能达到约 $100$–$200$ ns 的残差精度（最好的脉冲星）。

绿岸望远镜（GBT）（$100$ 米×$110$ 米，美国西弗吉尼亚）：Arecibo 坍塌后成为 NANOGrav 的主力望远镜。

CHIME（2018 年加入）：拓展了南天覆盖，并提供了多个新毫秒脉冲星的日常监测。

阿雷西博坍塌等于把灵敏度最高的那只耳朵关掉。此后绿岸必须承担原来由两台望远镜分担的监测，月度采样不能断，否则最低频端会立刻变盲。CHIME 的贡献主要是日常监测与发现新星，不是把十五年时间基线一次性补齐。甚大阵后来也加入计时，用来填阿雷西博留下的天空与频率空隙。

数据集的演化#

从 2020 年的 12.5 年数据集起，NANOGrav 就已探测到所有脉冲星共同的低频红噪声，但当时尚未明确看到 Hellings-Downs 角相关——这是区分引力波背景和其他共同误差（如时钟误差）的关键。

从共同红噪声到 Hellings-Downs，中间隔着「公共谱」与「角相关」两道门槛。所有脉冲星共享的低频红噪声，也可以来自时钟误差或太阳系星历，并不自动等于引力波。15 年数据的关键进步，是把两两脉冲星的相关函数拿去对照 Hellings-Downs 曲线。

时间基线多了约三年，星表多了二十一颗。低频端仍由「观测了多少年」决定，角相关的采样密度则由「有多少对」决定。两件事不能互相替代，所以 15 年不是把 12.5 年的同一张图再画一遍。

15年数据集：Hellings-Downs 相关的探测#

信号特征#

2023 年数据分析的主要结果（Agazie et al. 2023，The Astrophysical Journal Letters，951，L8）：

功率谱：在频率 $f \approx 1/T$（约 $1/15$ yr$^{-1}$，即约 $2$ nHz）到约 $10$–$20$ nHz 的范围内，探测到一个红噪声功率谱，可用幂律描述：

$$S(f) = \frac{A_{GWB}^2}{12\pi^2} \left(\frac{f}{f_{ref}}\right)^{-\gamma} f_{ref}^{-3}$$

拟合给出振幅 $\log_{10}(A_{GWB}) \approx -14.7$（在 $f_{ref} = 1/\text{yr}$ 处的应变振幅），谱指数 $\gamma \approx 3.2$（对应的空间应变功率谱斜率）。

Hellings-Downs 相关：计时星表有 68 颗毫秒脉冲星，角相关分析用其中 67 颗的两两组合（两千余对），计算不同角间距 $\theta$ 的计时残差相关性，结果与 Hellings-Downs 曲线的吻合显著性约 $2$–$4\sigma$（不同统计方法）。

论文把「引力波背景加 Hellings-Downs」的模型，与「只有各星独立噪声」相比，贝叶斯因子超过 $10^{14}$。与「有共同谱、但星与星之间不相关」（常记作 CURN）相比，贝叶斯因子约 200 到 1000，取决于谱怎么建模。用打乱星间相关的零假设去做背景分布，得到 $p \approx 10^{-3}$，大约三个标准差；频率派的最优统计给出 $p$ 约 $5\times 10^{-5}$ 到 $1.9\times 10^{-4}$，大约 $3.5$ 到 $4$ 个标准差。这些数字都来自同一篇 Agazie et al. 2023，它们描述的是证据强度，不是发现阈值已经被跨过。

计时数据集包含 68 颗毫秒脉冲星；Hellings-Downs 分析实际用了其中 67 颗的两两相关。假定特征应变谱按双星旋近的 $f^{-2/3}$，参考频率每年一周处的应变振幅中位数为 $2.4\times 10^{-15}$（90% 可信区间 $+0.7/-0.6$）。这与上文 $\log_{10}(A_{GWB})\approx -14.7$ 同量级，不是另一套互相打架的振幅。

Hellings 与 Downs 在 1983 年给出的，是各向同性引力波背景应有的角相关形状，不是 NANOGrav 为 2023 年临时发明的判据。共同谱只能说明「大家一起抖」；只有角相关随天球夹角按这条曲线先降后升、还会变号，才把抖指向时空应变，而不是指向钟差。

与超大质量黑洞双星模型的比较#

最自然的解释是宇宙中无数超大质量黑洞双星（SMBBH）发射的引力波叠加形成随机背景。对于圆轨道、引力波驱动衰减的 SMBBH，理论预言功率谱斜率 $\gamma = 13/3 \approx 4.33$（Phinney 2001）。

观测值 $\gamma \approx 3.2$ 比理论预言更平坦，意味着谱形在高频端偏平（相对于纯 GW 驱动的 SMBBH 预言）。可能的解释：

1. 环境效应：SMBBH 的轨道衰减除引力波外，还受到星系中心恒星（"硬化"）和气体的影响，可能在低频端（长周期）加快衰减，改变谱形
2. SMBBH 参数分布：若 SMBBH 质量函数、红移分布、偏心率分布偏离假设，谱斜率会变化
3. 不同来源：谱的平坦化也可能部分来自宇宙相变、宇宙弦等非 SMBBH 信号

最可能的源是尚未并合的超大质量黑洞双星背景，因为星系并合会自然留下大量仍在旋近的核。最后一秒差距问题问的是：双星如何从约一秒差距继续收紧，才能进入引力波主导的衰减。谱比纯引力波驱动更平，既可能是恒星与气体在核区帮了忙，也可能是偏心率或质量函数与最简模型不同。目前的数据还不能在这些解释里做终审。

把谱的平坦化直接读成「发现了宇宙弦」或「发现了一阶相变」，会把尚未指定来源的背景过早写成新物理。双星种群仍是目前最省额外假设的读法，只因为星系并合史会自然填满这个频段，不是因为其他来源已被实验关掉。

其他PTA团队的同期结果#

EPTA（欧洲脉冲星计时阵列）#

使用欧洲多台射电望远镜（WSRT/MeerKAT、Effelsberg、Lovell、Nançay、SRT），数据跨度约 $25$ 年，$24$ 颗毫秒脉冲星。DR2 数据集同样发现一致的引力波背景证据（Antoniadis et al. 2023，Astronomy & Astrophysics，678，A50），且与 NANOGrav 信号特征高度一致。

欧洲阵的时间基线比北美更长，脉冲星数目更少。更长的基线有利于最低频，更少的对则让角相关采样更稀。两套数据对上，说明信号不像是某一台望远镜的独特系统误差。

PPTA（帕克斯脉冲星计时阵列）#

使用澳大利亚 Parkes 64 米望远镜（2020 年后主要使用 MeerKAT 阵列增强），DR3 数据集结果（Reardon et al. 2023，The Astrophysical Journal Letters，951，L6）同样显示 Hellings-Downs 相关的证据，振幅与其他团队一致。

南天与北天看到同一量级的背景，是对「这只是北美钟差」的一次独立否证。Parkes 与 MeerKAT 的硬件、电离层与射电频率干扰环境都与绿岸不同，系统误差不容易以同一套角相关形态重现。

CPTA（中国脉冲星计时阵列）#

使用中国 FAST 500 米望远镜，对 57 颗脉冲星的约 $5$ 年数据进行分析（Xu et al. 2023，Research in Astronomy and Astrophysics，23，075024）。尽管数据时间跨度较短，FAST 的极高灵敏度使得单颗脉冲星计时精度极高（约 $50$–$100$ ns），在较短时间内已达到可与其他 PTA 竞争的灵敏度。

五年基线在纳赫兹最低频仍然偏短。CPTA 目前的力量是单星精度和新的天空覆盖，不是已经比别人看到更低频的功率。它提供的是另一套几乎独立的系统误差，而不是把 NANOGrav 的十五年再测一遍。

IPTA（国际脉冲星计时阵列）#

IPTA 是 NANOGrav、EPTA、PPTA 的联合体，合并数据可进一步提升灵敏度。对四个团队数据的联合分析将是下一步关键，因为它同时加厚时间基线与天空覆盖，而不是把四篇通稿里的显著性做平均。

四个团队用四套望远镜、四套噪声模型，方向一致，这才是「证据」站得住的理由。下一步真正加权重的，是把北美、欧洲、澳大利亚的时间基线叠到同一批星上。单阵的钟差或星历误差，很难在三套硬件里以同一套角相关形态重现；联合分析要回答的是这条曲线会不会在合并之后更稳，而不是再报一个更大的标题数字。

科学意义#

超大质量黑洞双星演化#

如果信号确认为 SMBBH 背景，它将是：

• 第一次直接探测到超大质量黑洞双星在引力波驱动下的轨道演化（而非单次并合事件）
• 对 SMBBH 质量函数、并合率和宇宙学演化的约束——反映了星系并合历史的"引力波档案"
• 对"最后秒差距问题"（last-parsec problem，SMBBH 如何从 $\sim 1$ pc 进一步收紧以发射引力波）的约束

「如果确认」这四个字不能省。背景是海量不可分辨源的叠加，没有一张宿主星系的光学照片可以一锤定音。要把统计结论钉成个案，还需要某一对足够近、足够重的双星在背景上冒出可分辨的周期性，并最好找到电磁对应体。那是下一阶段的目标，不是 2023 年已经交付的产品。

早期宇宙物理（如果不是 SMBBH）#

若信号的谱形不符合 SMBBH 预期，可能暗示宇宙早期的新物理：

• 宇宙弦：网络拓扑缺陷的振荡产生特定幂律引力波背景（$\gamma \neq 13/3$）
• 一阶宇宙相变：如果宇宙在冷却时经历了某个一阶相变（电弱、QCD 或新物理相变），碰撞的相变气泡产生引力波，谱型与 SMBBH 不同
• 原初引力波：暴胀期间放大的张量扰动，通常在 nHz 频段太弱，但某些非标准暴胀模型可能有显著信号

要把峰放到纳赫兹，相变的特征能量必须落在很窄的窗口。标准模型里的电弱与 QCD 转变并不是那种强一阶相变。原初引力波的标准慢滚预言在这一频段通常太弱；要让它被看见，需要非标准谱倾斜，并会在微波背景的 B 模式等窗口留下痕迹。解释顺序应当是：先看双星种群能否讲圆，再考虑把剩余部分分给早期宇宙。

对哈勃张力的潜在贡献#

若 SMBBH 背景的振幅和谱与对宇宙 SMBBH 质量函数的精确模型结合，未来可以约束宇宙合并率随时间的演化，进而约束宇宙学参数（$H_0$、$\Omega_m$ 等），提供独立于 CMB 和造父变星的新方法。

用纳赫兹背景去约束哈勃常数，目前只是前景，不是已经有的第三条尺子。振幅还缠着双星种群、星系并合率与环境耗散，参数简并很大。把它写成已经能裁决哈勃张力，会把尚未展开的推断说成测量。

未来展望#

FAST 的主导角色：FAST 每年对约 $60$–$80$ 颗毫秒脉冲星进行计时，灵敏度比 Arecibo 高 $3$–$5$ 倍，未来 $5$–$10$ 年的数据积累将大幅提升 Hellings-Downs 相关的置信度，并开始分辨引力波背景的各向异性（不同天区背景强度不同）。

SKA 的变革：平方千米阵列建成后，PTA 灵敏度提升约一个数量级，预计可以：

• 将 Hellings-Downs 相关置信度提升到 $> 10\sigma$
• 探测单个近距 SMBBH 系统的持续引力波
• 可能分辨背景各向异性，重建引力波背景的"天图"

阿雷西博已经不在。北美一侧要靠绿岸、甚大阵和未来的新望远镜续测。平方千米阵列若按计划在 2030 年代形成核心阵，才谈得上把角相关推到远高于目前证据的水平。单源连续波比随机背景更苛刻：那要求某一对超大质量黑洞足够近、足够重，才能在背景上冒出可分辨的周期。各向异性则是问背景亮不亮随天区变。两件事都还在未来，不能倒过来当成 2023 年已经完成的发现。

为什么这很重要#

NANOGrav 的 2023 年结果是引力波天文学的第二次大革命——在完全不同的频率窗口（nHz vs. 百 Hz）和完全不同的物理来源（宇宙背景 vs. 单次并合）上实现。LIGO 打开了"高频"窗口，PTA 打开了"低频"窗口，而 LISA 将在未来打开"中频"毫赫兹窗口。三个窗口合在一起，形成覆盖 $10^{-9}$–$10^3$ Hz 的完整引力波天文学频谱，每个窗口都有其独特的宇宙信息。

地面激光干涉仪看见的是可指认的并合波形；NANOGrav 看见的是多年积累的相关函数。窗口不同，物理不同，确认的门槛也应不同。2023 年应当写进教科书的是 Hellings-Downs 证据，以及「最可能来自超大质量黑洞双星背景」这一工作假设。后续要看的是角相关是否随时间基线稳定上升，以及谱形是否继续指向双星，而不是过早把宇宙弦写进标准图像。

跨域连接#

• 星系合并与相互作用：观测到的谱指数比纯引力波驱动的圆轨道双星预言更平，把星系核里恒星与气体对双星轨道的作用加进来会改变低频端的衰减速率。推论是：谱形不是黑洞的性质，而是"黑洞怎么被星系送到一起"的性质，纳赫兹背景因此变成对星系并合史的一次测量。
• 元分析与证据综合：同一份数据换一种统计方法，角相关的吻合显著性就在两到四个标准差之间移动；而四个团队用四套独立数据、四条独立管线得到了方向一致的信号。推论是：判据不在数字大小，而在系统误差是否独立——这就是社区称其为"证据"而非"发现"的理由。
• 随机过程：随机背景是海量不可分辨源的叠加，只能用功率谱和两点相关来描述；单个双星则是一条相位连贯的窄带信号。推论是：两者的判别方式完全不同，背景靠角相关，单源靠周期性——所以找到一个电磁对应体，是把统计结论钉成个案的唯一途径。
• 相变与临界现象：如果宇宙早期某次相变是一阶的，碰撞的相变气泡会辐射引力波，谱形与黑洞双星背景不同。推论是：这扇窗口同时是一台对撞机替代品，能约束远超地面加速器能量的早期物理；但目前的谱形既容得下黑洞双星，也容得下宇宙弦或相变，来源之争没有定论。
• 蒙特卡洛方法：要从几十颗星的残差里同时解出逐星噪声、公共谱与角相关，是一个高维贝叶斯采样问题，报出来的显著性本身也是采样估计的产物。推论是：换一组噪声先验，数字就会移动，这也是它与对撞机那种阈值不可直接对比的根本原因。

参考文献#

• Agazie, G. et al. (NANOGrav Collaboration) (2023). The NANOGrav 15-year data set: Evidence for a gravitational-wave background. The Astrophysical Journal Letters, 951, L8.
• Agazie, G. et al. (NANOGrav Collaboration) (2023). The NANOGrav 15 yr Data Set: Observations and Timing of 68 Millisecond Pulsars. The Astrophysical Journal Letters, 951, L9.
• Hellings, R.W. & Downs, G.S. (1983). Upper limits on the isotropic gravitational radiation background from pulsar timing analysis. The Astrophysical Journal Letters, 265, L39.
• Antoniadis, J. et al. (EPTA Collaboration) (2023). The second data release from the European Pulsar Timing Array. Astronomy & Astrophysics, 678, A50.
• Reardon, D.J. et al. (2023). Search for an isotropic gravitational-wave background with the Parkes Pulsar Timing Array. The Astrophysical Journal Letters, 951, L6.
• Xu, H. et al. (2023). Searching for the nano-Hertz stochastic gravitational-wave background with the Chinese Pulsar Timing Array. Research in Astronomy and Astrophysics, 23, 075024.
• Phinney, E.S. (2001). A practical theorem on gravitational wave backgrounds. The Astrophysical Journal Letters, 554, L37.

延伸阅读#

• Burke-Spolaor, S. et al. (2019). The astrophysics of nanohertz gravitational waves. The Astronomy and Astrophysics Review, 27, 5.
• Sesana, A. et al. (2016). Measuring the expansion of the universe with gravitational-wave standard sirens. Physical Review Letters, 116, 231102.
• Arzoumanian, Z. et al. (NANOGrav Collaboration) (2020). The NANOGrav 12.5-year data set: Search for an isotropic stochastic gravitational-wave background. The Astrophysical Journal Letters, 905, L34.