宇宙中的脉冲星不仅是“灯塔”,还是宇宙中最精确的“时钟”。脉冲星的脉冲来源于其自转。脉冲星质量大约是1.4个太阳质量,而半径仅为10公里左右。由于质量大,半径小,外部的扰动就很难有效地改变脉冲星的自转特性。实际观测发现,脉冲星的转动非常稳定,特别是毫秒脉冲星的长期转动稳定性甚至比现在投入使用的国际原子钟组性能还要优良。
这些散布在银河系中的、转动稳定的毫秒脉冲星即构成了可以通过射电天文观测,进而可以读出“时间”的高精密时钟。其良好的时间稳定性提供了脉冲星作为基本物理实验室的基础。目前已经测得的、最精确的脉冲星每亿年误差小于1秒。
射电天文观测如何去获得脉冲星时间的“读数”呢?这种技术在天文学上叫做脉冲星测时观测。它包括几个步骤:首先,人们用无线电望远镜接收脉冲星的无线电信号。
为了消除色散的影响,我们需要重新对齐不同频率的脉冲信号到达地球的时间。对于较亮的脉冲星,这个时候就能够像Hewish和Bell那样看到脉冲了。然而大部分脉冲星的信号实在是太弱了,即使消除了色散的影响,大部分情况下,脉冲信号还是不可见的。接下来让脉冲信号变得明显起来的办法叫做“周期折叠”,即按照脉冲的周期把数据分段,然后再每段数据叠加起来。这样由于将有脉冲的地方相叠加,脉冲信号就能变得明显起来。
但是仅仅获得脉冲到达地球的时间还不足以用来提取我们需要的物理信息。地球在自转,所以不同时刻,望远镜朝向天空的位置是不一样的。为了扣除地球的自转效应,我们先要计算出望远镜到地球中心的距离,把观测到的脉冲信号到达望远镜的时间归算到脉冲信号到达地心的时间。类似地,地球还在围绕太阳作公转。因此还需要进一步计算地球相对太阳系质心的位置,再把脉冲的达到时间归算到太阳系质心。
最后,在知道脉冲星相对太阳系的运动之后,我们还要进一步把太阳系质心时间折算到相对脉冲星平动的参考者那里。这样的脉冲到达时间就能表征脉冲星真实自转信息了。
脉冲星自转相当稳定。对于孤立的脉冲星来说,它的自转可以通过一个很简单的多项式模型来进行描述。测量到的脉冲时间和这个多项式模型的差异很小,但是由于测量误差和其它信号的存在,这个差异并不为零。我们把这个差异叫做“测时残差”。
测时残差包含了额外的物理信息。凡是那些不在我们上述描述中的、会影响脉冲到达时间的其它所有物理过程都会在测时残差里留下踪迹。另外,如上述描述中的物理参数不够准确,那么这些信息也会留在测时残差里边。
如果脉冲星在双星系统里边,上述过程得到的是在脉冲双星质心参考者观测到的脉冲到达时间。脉冲星和伴星还在相互绕转,进一步利用双星的轨道信息,我们把脉冲到达时间再次归算到相对脉冲星平动的观测者那里。在这个过程中,通过测量脉冲星到达时间的额外变化,就能获得脉冲星伴星的信息。例如首个系外行星系统PSR B1257+12就是这样发现的。
脉冲星双星系统和太阳系看似很像,但是脉冲星双星系统的问题要复杂得多。
图3描绘出了著名的脉冲双星系统B1913+16和J0737-3039A/B两个系统的轨道大小,我们按照相应比例画出了太阳的大小。可以看到,在这两个系统中,两个1.4倍太阳质量的脉冲星在与太阳大小相当的区域内作相互绕转。这样的系统由于引力很强,必须考虑广义相对论引起的效应才能得到准确的运动规律。
例如,众所周知的水星轨道由于广义相对论效应引起的近星点进动量为每百年43角秒,而J0737则达到了每年16度!这种相对论效应明显的系统有个学名叫做“相对论性双星”,而1993年的诺贝尔物理学奖和这种双星系统密切相关。
根据广义相对论,相对论性双星系统相互绕转之时将会辐射引力波而带走系统的能量。进而致使两颗星之间的距离会相应减小,而相互绕转的频率会升高。
通过脉冲星测时的方法,测量脉冲双星的轨道周期随时间变化的规律,再同广义相对论引力波带走能量而预言的轨道周期变化进行比较,即可检验广义相对论是否正确。在上世纪六十年代初期,科学界认为引力波探测是一项不可能完成的任务。然而脉冲双星系统B1913+16的发现及其轨道周期变化的测量间接地证明了引力波的存在、引力波会带走能量,同时也证明了引力波带走能量的方式与广义相对论预言一致。
不过,通过脉冲双星对引力波的测量仍是间接测量,因为我们测量的是脉冲星轨道周期的变化,进而推断出系统能量的减少,再而证明引力波存在的。所以脉冲星双星实验是间接推断了“引力波的存在”。或多或少受脉冲双星工作的影响,到了上世纪七十年代,科学界部份敏锐的研究人员已经意识到引力波探测也许不是遥不可及得了。
截至目前,双星系统对引力波的间接测量和引力检验的精度比上世纪七十年代已经高了不少,目前测量值和广义相对论预言的偏差小于0.1%。
现在让我们从引力波的间接探测回到脉冲星测时问题。根据前文介绍的内容,我们已经知道了三件事情:一、脉冲星会发出非常准确的周期性脉冲信号;二、我们能够利用望远镜来读出脉冲到达地球的时间;三、测时残差包含了其它的物理信息。那么哪些物理信息能够被我们观测到呢?我们又怎么提取这些信号?我们将通过两个前沿的应用例子给予说明:脉冲星时间-空间标准及引力波直接探测。
目前我们用的时间是从1972年1月1日开始通行的国际协调时。它是用原子钟定义的时间加上地球自转观测修正构成的。那么,国际协调时是如何确定的呢?首先,人们通过加权平均地球上运行的大约500个原子钟的时间构成国际原子时。在这个基础上,人们再通过闰秒修正获得国际协调时。在应用国际协调时的时候,因为每个原子钟的数据无法实时进行汇总,人们并不能实时地获得加权时间。精确时间测量仍然需要进行事后改正。
如果我们在脉冲星测时观测过程中没有考虑国际原子时的变化的话,那么脉冲星计时残差里边就会留下相应的信号。反过来,如果国际原子时本身有误差的话,那么我们也应当能够从脉冲计时残差里边看到这个误差。问题在于,每颗脉冲星也有自己的不稳定性,我们如何区分测时残差究竟是来源于脉冲星自身的不稳定性,还是来自国际原子时的误差呢?读者也许会自然而然地想到:我们同时观测了多颗脉冲星。
如果没有“阴谋论”的原因的话,可以预计,每颗脉冲星的自转特性是完全独立的,一颗脉冲星自转变化的时候,不会影响其它脉冲星的自转变化。所以如果测时残差来源于脉冲星自身,那么每颗星的变化应该是不一样的。但是如果测时残差来源于国际原子时的不稳定性,情况就不一样了。
不妨假设由于某些原因,国际原子时提前了一秒。
那么我们用这个国际原子时作为基准去测量脉冲星的脉冲到达时间,将会观察到所有脉冲星的到达时间都推迟了一秒钟,也即所有脉冲星的测时残差都会出现一秒的跳变。按照这个讨论,如果我们把脉冲星测时残差中公共的那部分提取出来,那么这部分信号就可以认为是国际原子时的系统误差。如能进一步再利用这个信号,对国际原子时进行修正,即完成了脉冲星钟和国际原子时之间的校正。
这种对多颗脉冲星进行测时观测并综合起来分析处理的技术叫做脉冲星测时阵列。
除了时间标准,我们还能够利用脉冲星建立空间框架。我们已经知道,在进行脉冲星测时观测的时候需要准确地知道地球在太阳系中的位置。如果我们对这个知识掌握的不够准确,那么就会有可观测的误差。因此我们可以在数据中搜寻太阳系中未知的小行星,或者对一些参数非常不确定的系统进行限制。在X射线波段也能看到脉冲星。
利用类似的技术,可以测量X射线卫星的位置,即实现X射线卫星的脉冲星导航。这种方式适合深空探测。因为在远离地球的地方,传统的GPS卫星定位或者干涉定位精度会逐渐变差,而使用脉冲星则可以在从太阳系到整个银河系的尺度上构造基本的时间和空间框架。不过地面设备支撑仍然必不可少,因为还需要知道地球相对于这个框架的坐标。
脉冲星测时阵列还能够直接探测引力波。
这里脉冲星被当成了标准的脉冲信号发生器来使用,其发出的脉冲信号穿过银河系的距离来到地球。如果宇宙背景中存在引力波,就会改变银河系的时空结构,从而改变信号的到达时间。广义相对论预言,引力波是时空的涟漪,而且是“横波模式”。这也就是说,在垂直于波传播方向上的空间被拉伸或者压缩。同时引力波还具有一种特性,叫做“无迹”。
这并不是说引力波不留下痕迹,而是指引力波压缩空间某个方向的时候,同时会拉伸相应的垂直方向。这样自由观测者围成的空间体积近似不变。由于引力波的“无迹”性,两颗处于垂直位置的脉冲星的脉冲信号具有“相反”的特点,也就是一颗星的脉冲到达地球如果变早了,那么另外一颗星的脉冲到达地球时间就会推迟。类似的推理不难得知,两颗位置相似的脉冲星的脉冲信号是同相的;而位置相反的两颗脉冲星脉冲信号也是同相的。
数学上有一种描述信号同相或相反的运算,叫做相关运算。完全一致的信号,计算出来的相关系数为1;如果是完全反相的信号,则计算出来的系数为-1。细致的计算表明:来自宇宙四面八方的引力波背景会导致每对脉冲星的信号都变得有相关性,而且这个相关性和脉冲星的之间的角距离是有关系的。目前国际上有大型射电天文望远镜的国家联合起来组成了国际脉冲星测时阵列合作项目。
其中包括欧洲的EPTA,澳大利亚的PPTA,以及美国的NANOGRAV。他们分别依托欧洲5个100米级望远镜、澳大利亚Parkes望远镜、美国Green Bank和Arecibo望远镜等国际上几个口径最大的望远镜展开工作。目前,国际脉冲星测时阵列的灵敏度已经逼近了做出“发现”的水平,极有可能在几年之内成功地探测到纳赫兹波段的引力波。
我国现已建成了世界上最大的500米单口径望远镜FAST,新疆在建世界最大的全可动高频110米口径望远镜QTT,还有可能建造世界最大的干涉阵列ChinaART。这些望远镜由于具有很大的接收面积,可以有效提高脉冲星测时的精度。可以预计,如果这些望远镜能够及早完成调试并投入使用,则脉冲星测时阵列探测引力波的领域将有重大的突破。