中子星研究的过往今来
|作者：俞云伟(华中师范大学天体物理研究所)
本文选自《物理》2021年第6期
摘要 中子星概念的形成既是人们对物质基本结构认识的一个自然推论，同时也是理解恒星演化的一个必然环节。自1967年通过发现射电脉冲星证实了中子星的存在以来，基于半个多世纪的多波段、多信使观测，人们已经发现了数以千记的多种类型中子星，了解了单个中子星的电磁辐射机制、中子星双星系统的相互作用以及双中子星系统的引力波辐射等等，并在多个方面为中子星的内部物质组分及其状态这一核心科学问题做出了观测限制。
关键词 中子星，恒星演化，脉冲星
1 中子星概念的形成背景
对于恒星乃至几乎所有天体而言，它们的一生始终是一部与万有引力相爱相杀的历史。从恒星诞生于星云、维持于主序核燃烧、灭亡于引力坍缩、归寂于致密天体，概莫如此。引力主宰着宇宙中所有亮暗物质的流动，而运动所导致的离心效应则是物质抵抗引力的天然法宝，这种运动既指宏观亦指微观。微观的无规运动在宏观上即表现为气体的压强。当一团足够巨大的星云在引力的作用下开始坍缩、分裂，分裂、坍缩，引力势能的释放将不断加速气体分子的无规运动，使其内能增加，逐渐导致温度升高、压强变大、坍缩减缓。在温度达到一定临界条件时，坍缩气体的核燃烧将被最终点燃，其导致的巨大压强将使引力坍缩完全停止，恒星得以形成[1，2]。此后，通过源源不断的热核聚变，恒星将长久地维持于高温的状态以保证其与引力达到平衡的气体压强。然而，伴随着不可逆的核聚变在亿万年后走向终结，恒星核心的气体压强最终将无以为继，它儿时的引力坍缩噩梦再次降临。随着恒星核心的引力坍缩，其半径将不断缩小，而物质的密度将不断增大。那么，这个过程的终点将在哪里呢？
1844年，F. Bessel发现天狼星的伴星是一颗具有太阳质量但很难被观测到的奇特天体[3]。数十年后，人们发现这个被称作“小狼”的天体具有天狼星千分之一的光度和大概两万开尔文的温度，据此它被定性为一颗白矮星。基于“小狼”的质量、光度和温度，我们容易发现白矮星的半径大概与地球相当，但其密度高达~106 g·cm-3，是一种前所未知的物质状态。1926年，基于E. Fermi和P. Dirac所发现的费米子量子统计规律，R. Fowler提出在白矮星密度下，从原子中电离出来的大量自由电子将能够形成非常强大的简并压，远高于白矮星温度下的热压强[4]。这个理论可以很好地解释具有强大引力场的白矮星为什么能够稳定存在。同时也使人们认识到，恒星核心的坍缩必将经历电子简并压的快速增长，从而可能阻碍坍缩过程的进一步发生，使其止步于一颗稳定白矮星的形成。换言之，观测到的白矮星正是恒星核心坍缩的一种自然产物。不同于普通的恒星，白矮星的支撑压取决于它的密度而非温度，因而其存在不再依赖于额外的能源供给。这个结局无疑使得恒星演化理论变得完整而美妙。然而，S. Chandrasekhar却意识到，随着密度的增大，电子气的费米能将很容易达到甚至超过它的静能量从而成为相对论性气体。在此情况下，电子简并压随密度增长的趋势将变得十分疲软，将无力支撑相应密度下的引力，其结果是白矮星的质量将存在上限，即钱德拉塞卡(Chandrasekhar)极限(1.44 M⊙)[5，6]。这一推论使得那些具有较大质量的恒星死后无所归处，引发了恒星演化理论的一次危机。
图1 恒星核心坍缩为中子星的基本过程示意图
1932年，J. Chadwick用阿尔法粒子(氦原子核)轰击铍，再用铍所产生的射线轰击氢、氮，首次发现了中子，进而使人们认识到原子核由中子和质子所组成[7，8]。恰在这一重大发现之前，L. Landau提出过一个设想，他认为比白矮星更致密的天体可能就像是一个巨大的原子核，其基本组成单元就是这种亚于原子核的未知粒子[9]。这一设想被认为是中子星概念原型的最早出处。随着中子的发现，W. Baade和F. Zwicky在1934年研究超新星现象的论文中首次明确提出了中子星这个名词，并正确地指出超新星现象应起源于大质量恒星向中子星转化的过程(图1)，该过程中释放出来的巨大引力势能正是超新星爆发的能量来源[10]。1939年，R. C. Tolman、J. R. Oppenheimer和G. M. Volkoff在广义相对论框架下，基于理想中子气简并压和引力的平衡建立了第一个定量的中子星模型(TOV方程)[11，12]，得到了中子星内部的物质分布情况以及它的质量和半径大小。此外，从观测角度考虑，人们也非常关心中子星表面的热辐射性质。通过研究它们的冷却机制和过程，人们发现中子星表面的温度将在百万年内主要维持在数十万到百万开尔文的量级[13]，因而其热辐射将出现在X射线能段。然而，由于中子星的半径只有十几公里，以当时的技术是完全不可观测的，因此人们也就很快对这种奇特天体失去了研究的兴趣。
2 中子星的多波段多信使发现
第二次世界大战后，受益于成熟的雷达技术，国际射电天文学迅速发展并在1960年代迎来了高潮。1967年，A. Hewish和他的研究生J. Bell发现了一种具有准确周期性的射电脉冲信号[14]。经过分析，人们认为这种周期性信号很可能来自于一种旋转天体的辐射(类似于灯塔的效应)，并将这种天体命名为脉冲星。更具体来看，只有当脉冲星具有像中子星那样的质量、体积和密度时，才能够成功解释观测到的射电脉冲辐射强度和发生频率。同年，在不知道发现脉冲星的情况下，F. Pacini指出，如果中子星具有很强的磁场并能够快速旋转的话，那么它们就可能发出低频的电磁波辐射，从而造成某种观测效应[15]。因此，人们认为射电脉冲星辐射的能量应主要来自于中子星的旋转能。其实在发现脉冲星之前两年，A. Hewish和S. Okoye还曾在蟹状星云中发现过一个具有高亮温度的奇特射电源[16]，后来便知道该射电源正是位于该星云中间的一颗中子星。蟹状星云和它中心的这颗中子星一起，是我国宋史中所记载的发生于1054年(宋至和元年)的壮观超新星爆发事件的遗物。无论如何，脉冲星的发现终于使中子星从一个理论猜想变成了一个可被实际观测的真实天体，无疑称得上是天文学史上的一个里程碑。脉冲星的发现者A. Hewish也因此被授予了1974年的诺贝尔物理学奖。这一重大发现重新激发了人们对于中子星的研究热情。在观测方面，脉冲漂移、脉冲缺失、周期跃变等越来越多的辐射特征被不断地揭示，其中还包括一些脉冲极度缺失的旋转射电暂现源。人们甚至认为，当前炙手可热的快速射电暴现象也很可能和中子星的射电辐射存在密切联系。在理论方面，天体物理学家们在上个世纪就已迅速建立起一整套描述射电脉冲星辐射机制的电动力学理论体系[17—19]。一般认为，中子星应具有以偶极为主的强磁场结构，如图2所示。在星体周围充斥着与星体共转的高度电离的等离子体，称为磁层，这里是各种脉冲辐射的产生之所。不过，当前的一些观测也表明，多极场结构有时候会变得非常显著，可以使其辐射性质更为丰富。
图2 脉冲星磁场结构和射电脉冲辐射示意图
1960年代同样也是X射线天文学的开启年代，以1962年R. Giacconi发现第一个宇宙X射线源Scorpius X-1为标志[20]。当时，便有人试图将此X射线源和中子星的热辐射相联系，这在一定程度上为理解1967年发现的射电脉冲星创造了理论预热的机会。而在1967年之后，中子星研究也就顺理成章地在X射线波段发展起来，使射电波段的孤立星研究扩展到X射线波段的双星系统。首先，通过检查ScorpiusX-1的X射线和光学观测结果，I. Shklovsky提出这些辐射应来自于处于吸积状态的中子星[21]。其后，R. Giacconi等人成功从X射线源Centaurus X-3中发现了一个周期为4.8 s的脉冲信号[22]，并指出这种X射线脉冲辐射是由于中子星从伴星或星际介质中吸积物质到星体表面上所导致的(图3)。在这种情况下，脉冲辐射的能量来源是吸积物质的引力势能而非中子星的旋转能。1975年，J. Grindlay和J. Heise更是从中子星X射线源中发现了两次短暂的X射线爆发(流量增加10左右)[23]。这些发现使人们对中子星双星系统产生了浓厚的兴趣。1982年，D. Backer等人发现了第一颗毫秒脉冲星[24]，每秒钟可转642次，被认为正是在双星系统中吸积加速的结果。与此同时，这种极快的旋转状态对脉冲星的质量和半径给出了极为严格的限制，进一步强化了脉冲星的中子星属性。值得注意的是，X射线脉冲星并不总是处于双星系统中，有时候也可以是孤立的，并具有旋转驱动和吸积驱动所不能解释的辐射光度。1992年，R. Duncan和C. Thompson最早在理论上研究了一类具有超朗道临界磁场的特殊中子星(称为磁陀星)[25，26]，可以为这些反常的X射线脉冲星以及其他一些软伽马射线重复暴现象提供很好的解释。在这些现象中，星体的辐射主要由磁能耗散驱动，而高度扭曲的强磁场所引发的一些不稳定性也将自然导致X射线暴等剧烈活动的频繁发生。
图3 中子星吸积伴星物质产生X射线辐射的示意图
针对中子星的X射线观测还进一步揭示，中子星对外的能量输出实际上不仅仅只通过磁层的脉冲辐射，实际上这可能只占中子星全部能量输出的一小部分。从整体上看，中子星的旋转能将主要通过以低频电磁波和等离子体(常以正负电子对为主)相耦合的坡印亭外流形式输出，并将在远大于中子星半径的尺度上逐渐转化为相对论性的中子星星风。在该星风的演变和最终击打到外围超新星抛射物的过程中，将可能产生明亮的脉冲星风云辐射。不妨让我们再次凝视蟹状星云，可以看到，在它的中间脉冲星风云的辐射正在X射线望远镜的镜头下熠熠生辉、璀璨夺目(图4)。不过，要在理论上精确重现坡印亭能流转化为相对论性星风的具体过程并不是一件容易的事情。一般相信，大尺度磁场的重联一定在其中扮演了重要的角色[27]。
再来看处于双星系统中的中子星，它们实际上也不一定总是会处于吸积的状态，很多时候也可能只是在做简单的轨道运动，观测表现为脉冲星的脉冲到达时间具有明显的轨道调制。此时，如果中子星和主序伴星的星风都很强烈的话，星风之间的相互作用就有可能造成可被观测的辐射信号。近年来，Fermi伽马射线望远镜所看到的不少伽马射线源便有可能属于这种星风相互作用系统[28]。在某些情况下，中子星的星风甚至还可能深度剥离主序伴星的包层物质乃至完全摧毁主序伴星，使其最终仅剩下一颗加速后的孤立毫秒脉冲星。更使人感兴趣的是，中子星双星系统进一步演化的结果还有可能导致双中子星系统的形成。1974年，J. Taylor和R. Hulse发现了第一个双中子星系统PSR B1913+16，其中一颗可以观测到脉冲辐射[29]。利用它的周期性信号，可以很好地限制两颗致密星绕质心公转的轨道参数。2003年，M. Burgay还首次发现了第一对双脉冲星系统PSR J0737-3039，从而可以更加精确地测定双星参数和更好地检验广义相对论效应[30]。
根据广义相对论，两个天体的相互绕转可以导致引力波辐射，辐射的强度高度依赖于系统的致密性。因此，双中子星系统被认为是宇宙中最理想的引力波辐射源之一(图5)。引力波辐射的能量来自于双星绕转轨道的引力势能。因此，随着引力波的持续辐射，双星系统的轨道半径和周期将变短。J. Taylor等人对PSR 1913+16做了持续数十年的跟踪观测，发现其轨道变化与广义相对论的预言高度一致，间接证明了引力波辐射的存在[31，32]。他们也因此获得了1993年的诺贝尔物理学奖。受此项发现的激励，从信号最强的双中子星并合事件中直接探测到引力波辐射成为了几代物理学家孜孜以求的科学目标。经过数十年的努力，2017年8月17日美国激光干涉引力波天文台(LIGO)和欧洲室女座引力波天文台(Virgo)终于首次探测到了来自于双中子星并合事件GW170817的引力波辐射[33]，完成了几代人的夙愿。在这一事件中，人们还同时观测到了与引力波信号成协的短时标伽马射线暴辐射和千新星辐射。这些多波段辐射信号对并合过程产生的抛射物性质乃至并合产物性质给出了重要的限制，开启了天文学多信使研究的新时代。
3 中子星研究的核心问题和探索途径
半个多世纪以来的天文观测不断为人们理解中子星的各种物理细节提供着重要的依据，尤其是它们的磁场、磁层结构以及产生各类辐射信号的机制，甚至是它们的大气层和内外壳层的组成成分以及局部的热核爆发过程等等。但是，时至今日，人们对于中子星最内部的物质状态却仍然知之甚少。而实际上，从中子星概念最初被提出开始，所有中子星研究的一个最终指向始终在于揭示中子星内部究竟由何种物质所组成并处于何种状态(即物态问题)。这是一个事关基本物理的核心科学问题。
1964年，M. Gell-Mann和G. Zweig提出了强子结构的夸克模型[34，35]。以此为基础建立起来的粒子物理标准模型为中子星的内部物质组分提供了更多的可能性。1984年，E.Witten等人计算发现，在很大的参数范围内由u、d、s三味夸克组成的奇异夸克物质可以是物质的真正基态，从而可能改变人们一直认为的铁56最稳定的观点[36，37]。在此设想下，中子星的物理本质也许有可能是奇异夸克星[38]或者是奇子(理论设想的由大量夸克所组成的类强子结构)星[39]，它们的基本属性同样可以和脉冲星观测相吻合。同时，奇异夸克星也可以表现出一些独特的属性，比如表面数千飞米范围内可能存在强电场，比如它们的质量和半径可以非常小等。有鉴于此，黄永锋等人提出，在一些正常质量大小的脉冲星周围可能存在奇异夸克行星[40]。无论如何，微观物理的不确定性极大地增加了中子星研究的自由度。除了传统的中子星和上述提到的奇异夸克星外，还有超子星、混合星/混杂星、中子超流、质子超导、介子凝聚和夸克色超导等一系列新的概念被陆续引入到中子星研究中，它们之间的最大区别主要反映在对星体内核组分的不同构造上(图6)[41]。因此，当我们现在再讲中子星这个名词的时候，除非加以特别说明，它实际上常常会是一个非常广义的概念，包含了很多不同的物理内涵。这种现状使得针对中子星的天文观测成为了人类探索微观物理的新途径，并具有地面实验室所不具有的特殊优势。
原则上人们可以从多个角度去设法鉴别中子星的不同物质属性，但最有力的观测判据可能还是在于确定中子星的质量上限及其质量半径关系，尽管这是一件非常困难的工作。近年，人们利用测量双星轨道的Shapiro延迟等方法，找到了数例具有较大质量的中子星，其中最大的质量达到了
近二十年来，随着时域天文研究的不断发展，人们对超新星和伽马射线暴等恒星爆发现象也有了越来越丰富的认识。1998年，戴子高和陆埮提出，一部分伽马射线暴的残留中心致密天体可能正是一颗处于极限旋转状态的磁陀星(即毫秒磁陀星)[47，48]，其理论预言受到了数年后Swift卫星观测结果的强力支持[49]，从而使中子星研究在暂现源现象中大放异彩。近年来发现的一类极为明亮的超亮超新星便是一个重要的例证[50]，它们的光度演化常常可以与毫秒磁陀星的自转能损相一致，却很难由传统的镍56衰变理论解释。更为重要的是，戴子高等人还在2006年提出，双中子星并合后的产物可能仍然是一颗大质量的中子星[51]，这一观点对于理解短伽马射线暴的诸多观测特征具有重要帮助。尤其是通过分析GW170817引力波事件中的千新星辐射，俞云伟等提出该事件的并合产物很可能是一颗大质量中子星[52]。如果事实确实如此，那就意味着中子星的极限质量可以高达 2.5 M⊙ 以上，这将对中子星的物态和起源提供极为强烈的限制。由此可见，双中子星并合产物的属性无疑将是未来引力波探测和多信使天文学研究有待解决的一个重大问题。双中子星并合产物可能是大质量中子星的设想同时也表明，宇宙中的中子星可能具有多种不同的起源[53]，不单单只来自于超新星爆发。不同起源的中子星甚至可能具有不同的物态(即有些星体可能处于亚稳态)。实际上，早在1976年R. Canal和E. Schatzman就曾指出，吸积白矮星在质量趋近于Chandrasekhar极限的时候，有可能经吸积诱导坍缩为一颗中子星[54]，而不一定如通常认为的总是导致Ia型超新星爆发。目前观测上也的确发现了不少起源未知的暂现源现象，为人们提供了不少遐想的空间。有鉴于此，我们有理由相信当前以及未来的暂现源观测能够为中子星研究打开一扇全新的窗口，对于了解中子星初生时期的状态具有至关重要的价值。
4 结束语
综上所述，中子星的研究既久远悠长，又生机盎然、热点纷呈，在它们的身上仍然存在着一系列的谜团有待探索，并为人们研究极端条件(高密度、强引力场、强电磁场)下的物理规律提供了天然的实验室。同时，现代天文观测技术的快速发展正逐渐使得这种探索达到一种前所未有的广度和深度，而其中我国自主研发的一些大科学装置也正逐渐发挥出特有的效力。2020年4月28日，我国的慧眼HXMT望远镜成功地从磁陀星SGR 1935+2154观测到了与快速射电暴FRB200428成协的X射线暴[55]，为揭示快速射电暴现象的起源做出了关键贡献，也使人们对磁陀星的性质有了全新的认识。而就在本文完成之际，我国500 m口径球面射电望远镜(FAST)也再次传来好消息，它通过开展银道面脉冲星巡天，新发现了212颗脉冲星，其中包括42颗毫秒脉冲星、16颗脉冲双星、一批最暗弱的脉冲星、一批模式变化和消零脉冲星，以及射电暂现源等[56]。这些重大的发现无疑将为中子星研究带来新的高潮。
最后，让我们回到本文开头的那一段描述，因为中子星无论如何都存在着质量上限，因此当那些最大质量的恒星在坍缩的时候，其核心的密度终将远远超过核密度。然而，此时再讨论物质的状态似乎已变得了无意义，因为所有的物质信息都将禁闭在一个无限弯曲的时空中而不可窥测。也就是说，它将最终不可避免地坍缩为一个黑洞。当然，这将是另一个波澜壮阔的新故事了。
[1] Henyey L G，Lelevier R，Levée R D. PASP，1955，67：154
[2] Hayashi C. PASJ，1961，13：450
[3] Bessel F W. MNRAS，1844，6：136
[4] Fowler R H. MNRAS，1926，87：114
[5] Chandrasekhar S. MNRAS，1931，91：456
[6] Chandrasekhar S. MNRAS，1935，95：207
[7] Chadwick J. Nature，1932，129：312
[8] Chadwick J. Proceedings of the Royal Society of London SeriesA，1932，136：692
[9] Landau L. PhysikalischeZeitschrift der Sowjetunion，1932，1：285
[10] Baade W，Zwicky F. Proceedings of the National Academy ofScience，1934，20：254
[11] Tolman R C. Physical Review，1939，55：364
[12] Oppenheimer J R，Volkoff G M. Physical Review，1939，55：374
[13] Zwicky F. Physical Review，1939，55：726
[14] Hewish A，Bell S J，Pilkington J D H et al. Nature，1968，217：709
[15] Pacini F. Nature，1968，219：145
[16] Hewish A，Okoye S E. Nature，1965，207：59
[17] Goldreich P，Julian W H. ApJ，1969，157：869
[18] Ostriker J P，Gunn J E. ApJ，1969，157：1395
[19] Ruderman M. ARA&A，1972，10：427
[20] Giacconi R，Gursky H，Paolini F R et al. Phys. Rev. Lett.，1962， 9：439
[21] Shklovsky I S. ApJ，1967，148：L1
[22] Giacconi R，Gursky H，Kellogg E et al. ApJ，1971，167：L67
[23] Grindlay J，Heise J. IAU Circ.，1975，2879
[24] Backer D C，Kulkarni S R，Heiles C et al. Nature，1982，300：615
[25] Duncan R C，Thompson C. ApJ，1992，392：L9
[26] Thompson C，Duncan R C. ApJ，1993，408：194
[27] Drenkhahn G. A&A，2002，387：714
[28] Chen A M，Takata J，Yi S X et al. A&A，2019，627：A87
[29] Taylor J H，Hulse R A，Margon B et al. IAU Circ.，1974，2704
[30] Burgay M，D’Amico N，Possenti A et al. Nature，2003，426：531
[31] Taylor J H，Weisberg J M. ApJ，1982，253：908
[32] Taylor J H，Weisberg J M. ApJ，1989，345：434
[33] Abbott B P，Abbott R，Abbott T D et al. Phys. Rev. Lett.，2017，119：161101
[34] Gell-Mann M. Physics Letters，1964，8：214
[35] Zweig G. An SU(3) Model for Strong Interaction Symmetry and its Breaking. 1964. CERN-TH-401
[36] Witten E. Phys. Rev. D，1984，30：272
[37] Farhi E，Jaffe R L. Phys. Rev. D，1984，30：2379
[38] Alcock C，Farhi E，Olinto A. ApJ，1986，310：261
[39] Lai X Y，Xu R X. in Journal of Physics Conference Series，2017，861：012027
[40] Huang Y F，Yu Y B. ApJ，2017，848：115
[41] Haensel P，Potekhin A Y，Yakovlev D G. Neutron Stars 1：Equation of State and Structure(Astrophysics and Space Science Library，Vol. 326). New York：Springer，2007
[42] Cromartie H T，Fonseca E，Ransom S M et al. Nature Astronomy，2020，4：72
[43] Guerra Chaves A，Hinderer T. Journal of Physics G NuclearPhysics，2019，46：123002
[44] Miller M C，Lamb F K，Dittmann A J et al. ApJ，2019，887：L24
[45] Kaaret P，Prieskorn Z，in′tZand J J M et al. ApJ，2007，657：L97
[46] Shternin P S，Yakovlev D G，Heinke C O et al. MNRAS，2011，412：L108
[47] Dai Z G，Lu T. A&A，1998，333：L87
[48] Dai Z G，Lu T. Phys. Rev. Lett.，1998，81：4301
[49] Zhang B. Chinese J. Astron. Astrophys.，2007，7：1
[50] Gal-Yam A. Science，2012，337：927
[51] Dai Z G，Wang X Y，Wu X F et al. Science，2006，311：1127
[52] Yu Y W，Liu L D，Dai Z G. ApJ，2018，861：114
[53] Yu Y W，Chen A，Dai Z G et al，AIP Conference Proceeding，2019，2017：020024
[54] Canal R，Schatzman E. A&A，1976，46：229
[55] Li C K，Lin L，Xiong S L et al. Nature，in preprint，2020，arXiv：2005.11071
[56] Han J L，Wang C，Wang P F et al. 2021，arXiv:2105.08460
《物理》在淘宝店和微店上线，扫码即可购买过刊和现刊。
▼往期精彩回顾▼
2.非线性光学超构表面3.降低超导体的高压4.光学超构表面中的复合相位调控5.转角石墨烯的三昧｜量子多体中的呐喊与彷徨之九6.关于加强物理资助的对话7.原本需要15万年才能解决的问题，有了它只需1秒 | 彭新华——周末讲坛8.芯片能够参与竞争的任何技术，最终都不可避免地成为失败者 | 魏少军——科学讲坛