中国天眼FAST（500米口径球面射电望远镜）又立一功！近日，由中国科学院国家天文台等单位科研人员组成的中国脉冲星测时阵列（CPTA）研究团队利用中国天眼FAST，探测到纳赫兹引力波存在的关键性证据，表明我国纳赫兹引力波研究已与国际同步达到领先水平。相关论文在线发表于我国天文学术期刊《天文与天体物理研究（RAA）》上。
(资料图片仅供参考)
人类终于站在了期盼已久的纳赫兹引力波宇宙观测窗口前。那么，FAST是怎么探测到引力波的？纳赫兹引力波探测与获得2017年诺贝尔奖的引力波探测有何不同？此次研究成果有什么重大意义？我们今天就来详细了解。
引力波示意图
中国天眼FAST示意图
“时空的涟漪”助寻宇宙暗物质
想知道引力波为何物，得先从一百多年前说起。
1915年，爱因斯坦发表了广义相对论，这个全新的理论告诉我们，引力是时空弯曲产生的结果——时空告诉物质如何运动，物质告诉时空如何弯曲。简言之，广义相对论里的引力不是“力”，是时空的一种弯曲效果。
我们可以把时空想象成平静的水面，有质量的物体能让时空本身发生弯曲。如果宇宙中有质量的物体发生剧烈的加速运动，如超新星爆发、两个巨大星体发生碰撞、两个致密星体并合（中子星与黑洞这类星体的碰撞与并合），就像是往平静的水面扔进了石头，让时空本身泛起“涟漪”。“涟漪”携带着波源天体激烈动荡的信息和关于引力本质的线索，以光速向远处传播，这种时空“涟漪”就是引力波。
宇宙中发光的物质非常少，大部分（95%左右）都是真正的、永恒的“黑暗”——暗物质和暗能量，不管造多么先进的望远镜，这些“黑暗”人类都看不到。而引力波给人类提供了一种探索宇宙的新方式，让我们知道在看不到的宇宙空间里发生了什么。
引力波可以直接跟踪宇宙中有质量物体的运动（不管是否可见），利用引力波观测，我们就能够捕捉到“黑暗”的蛛丝马迹。比如，2015年人类第一次探测到的引力波GW150914，就是距离地球10亿光年之外的一个遥远星系中两个恒星级黑洞发生了并合。在以前，人类能知道这个过程，是不可想象的。
引力波不但大大扩展了我们对宇宙的感知，对物理学也有重要意义。通过对引力波的探测，可以检验目前的引力理论，检验人类对时间和空间的基本认知是否正确。
探测引力波需要不同的“尺子”
探测引力波的基本原理是“引力波所经过的区域，空间的长度会被周期性地拉伸和收缩”。换言之，引力波经过时，与之垂直的平面会处于不断伸缩的状态：横向收缩，纵向就拉伸；横向拉伸，纵向就收缩。空间本身的形变会让置身其中的物体也跟着形变。一旦探测出这种形变，就能说明观测到了引力波。
也许有人会问，美国科学家不是已经借助LIGO（激光干涉引力波天文台）探测到引力波，还获得诺贝尔奖了吗？为什么还需要其他装置（如FAST）来探测呢？答案是：不同的天文事件会产生不同频率的引力波，需要不同的“尺子”来测量。引力波的频率越低，波长就越长，所以承担测量功能的“尺子”也越长。
频率在10-1000赫兹波段的引力波由恒星级双黑洞、双中子星并合产生，这些波源离地球相对较近，通过前文提到的激光干涉引力波天文台可以探测到，它的“尺子”是两条相互垂直、长度达4公里的管道，激光在管道中穿行。具体来说，在引力波的影响下，激光走过的“路途”长度会发生规律性的变化，而激光干涉仪能把引力波导致的长度变化，转变为激光干涉结果的光强变化，因此可以通过测量激光干涉结果的变化来捕捉引力波。
把同样原理的激光干涉装置放到天上，“尺子”可长达数百万公里，捕捉毫赫兹频段的引力波（由双星绕转、大质量黑洞俘获致密星产生的）。欧洲的LISA计划和中国的太极计划、天琴计划等空间引力波探测装置，都是用卫星组成干涉仪网络进行长距离干涉测量的。
那么，想要探测到宇宙中更远处由超大质量双黑洞、宇宙弦引发的更低频率引力波（也就是周期在年量级的纳赫兹引力波），就需要更长的“尺子”了。目前已知的唯一探测手段是利用大型射电望远镜观测宇宙中的脉冲星，很多人可能无法想象，中国天眼FAST探测引力波时利用的脉冲星测时阵列其实有银河系那么大。至于极低频段（10-16赫兹）的引力波探测，就要利用宇宙微波背景辐射了，比如南极BICEP2、西藏阿里观测项目等。
所以，各种探测装置分别探测不同的引力波信号，彼此互为补充，不能相互替代。
探测纳赫兹引力波有多难
对频率低至纳赫兹的引力波进行探测，可以观测更深远的宇宙，有助于天文学家理解宇宙结构的起源，探测宇宙中最大质量的天体即超大质量黑洞的增长、演化及并合过程，也有助于物理学家洞察时空的基本物理原理。
当前，想要探测纳赫兹引力波，必须利用脉冲星，准确地说，是一群脉冲星。
脉冲星是一类磁场强且高速自转的中子星，它的自转很稳定，每隔固定的时间就会发出一个脉冲信号。如果不受其他因素影响，那么我们在地球上就能稳定地收到这种信号。
脉冲星的脉冲到达地球的时间极其规律，而引力波所经过的区域，空间的相对长度会被周期性地拉伸和收缩，这自然催生了一种设想：如果一列引力波经过了地球和脉冲星中间的区域，那么脉冲星发出的信号就会发生改变，从而改变该信号到达地球的时间。这样一来就很好理解了，只要观测到这种脉冲信号到达地球时间的变化，就等于探测到了引力波。这就是利用脉冲星测时（PTA）探测引力波的基本原理。
为什么不能只观测一颗脉冲星来探测引力波呢？因为仅凭一颗脉冲星，我们不知道探测到的信号是引力波、脉冲星噪声，还是其他的假信号，必须利用多颗脉冲星相关的方法，来寻找相关信号以区分引力波和噪声。如果发现多颗脉冲星同时发生某种规律性变化，就可以证明探测到引力波了。其中，利用大型射电望远镜对一批自转极其稳定的毫秒脉冲星进行长期测时观测，是纳赫兹引力波目前已知的唯一探测手段。
纳赫兹引力波探测，在国际上是一个竞争激烈的领域，相关思想早在1983年便已提出。国际上的探测团组有北美纳赫兹引力波天文台（NANOGrav）、欧洲脉冲星测时阵列（EPTA）、澳洲帕克斯脉冲星测时阵列（PPTA），以及印度脉冲星测时阵列（InPTA）、南非脉冲星测时阵列（SAPTA）和中国脉冲星测时阵列（CPTA）。
但是，由于纳赫兹引力波引起的时空改变非常微弱，引力波的周期也达到了年量级（时空的改变以年度为单位）。对它的探测在物理尺度和时间尺度上都很“大”，不但需要造银河系尺度探测器（脉冲星测时阵列），还需要观测很长时间，才能把时间周期长达数年的信号给找出来。
在此次中国天眼FAST探测到纳赫兹引力波存在的关键性证据之前，国际上还没有一个团组得到过确切的探测结果。
“天眼”拿到了什么关键证据
2016年6月，中国科学院启动了纳赫兹引力波预研究，联合北京大学、中国科学院新疆天文台、云南天文台、上海天文台、国家授时中心、广州大学等多家相关单位组建了中国脉冲星测时阵列研究团队。2019年上半年，中国天眼FAST开展了试观测，由此拉开了探测纳赫兹引力波的序幕。
不过，脉冲星测时阵列探测纳赫兹引力波的灵敏度强烈依赖于观测时间跨度，即灵敏度随着观测时间跨度的增长而迅速增加。
NANOGrav、EPTA、PPTA利用各自的大型射电望远镜，已分别开展了长达20年的纳赫兹引力波搜寻。在这个领域，中国脉冲星测时阵列研究团队虽然是“后来者”，但以数据精度、脉冲星数量和数据处理算法上的优势，弥补了时间跨度上的差距。
中国天眼FAST是目前全球最大且最灵敏的射电望远镜，也是全球搜寻脉冲星效率最高的射电望远镜，截至目前，已发现740余颗新脉冲星。中国脉冲星测时阵列研究团队面对观测时间跨度远短于美、欧、澳3个国际团队的不利局面，充分利用FAST灵敏度高、可监测脉冲星数目多、测量精度更高的优势，长期系统地监测了57颗毫秒脉冲星，并将这些毫秒脉冲星组成了银河系尺度大小的引力波探测器，来搜寻纳赫兹引力波。团队自主开发独立数据分析软件，对中国天眼FAST收集的时间跨度3年5个月的数据进行分析研究，发现了具有纳赫兹引力波特征的四极相关信号的证据。
本次测量到的引力波特征幅度很小，约为4×10-15，这意味着，它造成的时空改变极其微弱。在距离上，1公里尺度引力波引起的扰动约为百分之一个氢原子大小；在时间上，千万年尺度上才变化1秒。也就是说，对纳赫兹引力波的探测，是在挑战人类精密测量的极限。
本次关键性证据由4个国际团队分别独立获得、独立发表、相互印证，说明我国纳赫兹引力波探测灵敏度达到了与美、欧、澳相当的水平，从而同时实现此次重大科学突破。
中国脉冲星测时阵列研究团队暂时无法确定纳赫兹波段引力波的主要物理来源，其他国际团队也无法确定。不过，这个问题将随着中国天眼FAST后续观测数据时间跨度的增加得到解决。由于中国脉冲星测时阵列研究团队现有数据时间跨度较短，所以数据时间跨度增长带来的效果会更明显，例如，如果数据时间跨度再增长3年5个月，中国脉冲星测时阵列研究团队的数据时间跨度将翻倍，而其他国际团队仅增长不到20%。
未来将建成纳赫兹引力波天文台
作为目前世界上最大、最灵敏的单口径射电望远镜，中国天眼FAST将让中国在射电天文领域领先世界20年。从跟跑到领跑，这背后，是中国不断崛起的科技力量。
未来我们还将看到中国天眼FAST进行扩展和升级，基于脉冲星测时阵列方法实现纳赫兹引力波事件的常规观测，从而建成纳赫兹引力波天文台，并开启更高灵敏度和更高分辨率的低频射电观测研究新纪元。
2003年6月29日，一个中国的博士生前往澳洲，在那里开始了他利用脉冲星探测引力波的科研生涯。20年后的今天，当年的博士生已经成为中国脉冲星测时阵列研究团队的领军人物，利用中国天眼FAST发现了纳赫兹引力波存在的关键性证据，他就是李柯伽研究员。他的经历也许又一次印证了那句话：“FAST最大的意义是把人类的视野向前扩展了一大步，也让中国的天文学家第一次站在了人类视野的最前沿。”
没有人知道中国天眼FAST还能给天文学理论带来怎样的贡献，但可以确定的是，它一定会有更多贡献。
（作者单位：中国科学院计算机网络信息中心；专业审核：中国科学院国家天文台）
