双中子星并合与千新星

概述#

2017 年 8 月 17 日是天文学史上最重要的日子之一。这一天，人类第一次同时接收到来自同一天体事件的引力波、伽马射线、紫外/可见/红外光和X 射线/射电波——共五个电磁波段的信号，加上引力波，共六种"信使"。这个事件就是 GW170817，来自约 $130$ 百万光年（$40$ Mpc）外的一对双中子星的并合。

它证实了物理学家半个世纪来的预言：中子星并合是宇宙中重元素（金、铂、铀等）r 过程核合成的主要场所，并产生一种被称为千新星（kilonova）的多波段电磁瞬变——其亮度比普通新星亮约 $1000$ 倍。

千新星这个名字容易让人以为它是缩小版超新星。能源是新合成核素的衰变热，不是镍-56 主导的几周光变；看错能源，就会用超新星模板去减它，把 r 过程信号减没。它也不是“看见了金线”的同义词：先有光变与不透明度，后有锶的吸收证认。

常见误解是把金、铂写成超新星的必然产物。铁以上多数稳定核走的是快中子俘获，超新星风是否够格，争论了几十年。GW170817 把双中子星并合从候选名单推到了证据最硬的现场，但没有证明它是宇宙史上唯一的 r 过程工厂。机制与核素账本见千新星与 r 过程。

六种信使叠在同一个时间戳上，靠的不是运气，而是定位。三台引力波探测器把天区压到可搜的几十平方度，光学巡天才有机会在约十一小时内找到 AT2017gfo。引力波探测器网络与时域天文学在这里是同一条链，也是多信使天文学真正落地的那一次。

双中子星：从诞生到并合#

双中子星系统（binary neutron star, BNS）通过以下路径形成：

1. 大质量双星中的第一颗恒星演化完成，通过超新星爆炸形成中子星。超新星爆炸往往给中子星一个"踢速"（natal kick），可能解散双星，但部分系统保持束缚
2. 第二颗恒星演化为巨星，通过物质转移和公共包层演化阶段（common envelope phase），双星轨道收紧，第二颗超新星形成第二颗中子星
3. 双中子星系统发射引力波，轨道能量损失导致轨道衰减，最终在约 $1$–$10^{10}$ 年后并合

公共包层是这条路径最窄的关口。第二颗星膨胀后把第一颗中子星吞进去，轨道摩擦把双星拖近，否则引力波辐射来不及在宇宙年龄内把它们并合。超新星踢速则是另一道筛：踢得太猛，系统散掉；踢得刚好，才留下能并合的中子星双星。通道细节见双星演化。

1974 年，拉塞尔·赫尔斯（Russell Hulse）和约瑟夫·泰勒（Joseph Taylor）发现的 PSR B1913+16（赫尔斯-泰勒脉冲星）是第一颗已知的双中子星系统，其轨道衰减率与广义相对论的引力波辐射预言完全吻合（Weisberg & Taylor 2004），为 Hulse 和 Taylor 带来了 1993 年诺贝尔物理学奖——这是引力波的第一个间接证据。

赫尔斯-泰勒脉冲星并没有在我们眼皮底下并合。它证明的是轨道在按广义相对论预言的速率收缩，引力波在被直接听到之前，已经写进了脉冲到达时刻。GW170817 把同一本能量账从间接变成直接：这次是波形本身，而不只是轨道周期的导数。

GW170817：信号链#

引力波信号#

2017 年 8 月 17 日 12:41:04 UTC，LIGO Hanford、LIGO Livingston 和 Virgo 三台探测器几乎同时接收到持续约 $100$ 秒的引力波信号——双中子星旋近的"啁啾"（chirp）信号。波形分析给出：

• 两颗中子星的啁啾质量（chirp mass）$\mathcal{M}_c = 1.186^{+0.004}_{-0.002}\,M_\odot$（源参考系；探测器参考系约 $1.1975\,M_\odot$）
• 总质量约 $2.74^{+0.04}_{-0.01}\,M_\odot$
• 各自质量约 $1.17$–$1.60\,M_\odot$
• 天空定位面积约 $28\,\text{deg}^2$（90% 置信）

（Abbott et al. 2017，Physical Review Letters，119，161101）

潮汐形变是波形里多出来的那一笔。两颗中子星在并合前互相拉扯，星体越容易变形，波形越早偏离点质量预言，这直接约束中子星内部状态方程。半径仍有模型依赖，但“太软撑不住、太硬形变太小”的两端都被这次事件收紧了。

Virgo 的加入把天区从单台探测器的上千平方度压到约 $28$ 平方度。没有第三台，光学望远镜无法在一夜之内把候选星系搜完。多信使的第一层联合，发生在引力波网络内部。

伽马射线暴#

引力波结束 $1.7$ 秒后，Fermi 和 INTEGRAL 卫星探测到一个微弱的短伽马射线暴 GRB 170817A，持续约 $2$ 秒，能量约 $5 \times 10^{46}$ erg（Goldstein et al. 2017，The Astrophysical Journal Letters，848，L14）——比典型短 GRB 弱约 $1000$ 倍，可能因为喷流方向偏离视线（"偏轴观测"）。这是首次将短伽马射线暴与双中子星并合直接联系起来，证实了此前 20 年的理论预言：短 GRB 来自中子星并合（Eichler et al. 1989；Paczynski 1991）。

这段延迟本身就是物理结果。伽马射线若在并合后由喷流产生，引擎启动、喷流冲出周围物质都需要时间；把已测到的约 $1.7$ 秒放到约 $40$ Mpc 的路程上，引力波传播速度与光速的相对差被压到 $10^{-15}$ 量级（Abbott et al. 2017，The Astrophysical Journal Letters，848，L13）。这几乎排除了一大批要求引力波显著慢于光的修正引力。

偏轴不是失败，是几何。喷流有结构：核心亮、翅膀弱。视线若从核心擦过，伽马射线会弱两三个数量级，光学千新星却几乎各向同性。所以引力波才是无偏触发器，短暴样本则严重偏向正对观测者的那一小撮。机制见伽马射线暴与吸积盘与相对论喷流。

千新星：AT2017gfo#

引力波探测后约 $11$ 小时，多个光学望远镜在宿主星系 NGC 4993（距离约 $40$ Mpc）中发现了一个迅速变亮的瞬变源 AT2017gfo（也称为 SSS17a）。它的光变曲线与普通超新星完全不同：

• 蓝色成分（blue kilonova）：峰值在近紫外，约 $1$ 天后达到峰值，随后在约 $3$ 天内迅速变红并暗淡——对应轻 r 过程元素（如锶 Sr、钡 Ba）的不透明度较低的喷出物
• 红色成分（red kilonova）：峰值在近红外，约 $3$–$5$ 天后仍保持较亮，持续时间更长——对应重 r 过程元素（如金 Au、镧系元素 Lanthanides）的高不透明度（$\kappa \sim 10$–$100\,\text{cm}^2/\text{g}$）喷出物

这种双成分结构正是理论预言的千新星特征（Metzger et al. 2010，Monthly Notices of the Royal Astronomical Society，406，2650；Roberts et al. 2011）。

光学发现来自智利拉斯坎帕纳斯的斯沃普 1 米望远镜。Coulter 等人在引力波警报后约 $10.9$ 小时，在 NGC 4993 里找到迅速变亮又迅速变红的源 SSS17a，后由 IAU 命名为 AT2017gfo（Coulter et al. 2017，Science，358，1556）。它不是超新星：几天内就从蓝变红，光变时标比 Ia 型短一个数量级。

Kasen 等人用双成分不透明度模型拟合这次光变：轻 r 过程物质（质量数 $A < 140$）给出早期蓝光，含镧系的重 r 过程物质给出持续近红外（Kasen et al. 2017，Nature，551，80）。镧系的开壳层电子让不透明度比铁族高一个到两个数量级，光子被困得更久，峰值被推到红外、被推迟到数天之后。这正是“千新星比新星亮约千倍、却比超新星快而暗”的物理来源。

千新星与 r 过程核合成#

r 过程：宇宙的重元素工厂#

快中子俘获过程（rapid neutron capture process，r 过程）在极短时间内（$\sim 1$ 秒）向种子核添加大量中子，随后 $\beta$ 衰变形成稳定的重元素，是合成铁以上大多数重核（$A > 80$，包括 Kr、Sr、Ba、La、Au、Pt、Pb、U、Th 等）的主要途径。

r 过程需要极端的中子密度（$n_n > 10^{20}\,\text{cm}^{-3}$）和高温，可能的场所长期有争议：

• 核心坍缩超新星（CCSN）：中微子驱动的风中可能有 r 过程条件，但模拟表明产量受限
• 双中子星并合（BNS merger）：并合喷出的中子星物质（中子比例极高，$Y_e < 0.3$）天然提供 r 过程条件——这是目前证据最充分的 r 过程主场所

GW170817 的千新星喷出物质量约 $0.04$–$0.05\,M_\odot$，重元素产率与银河系中 r 过程元素丰度比较，表明双中子星并合是宇宙中大部分金、铂等重元素的来源（Drout et al. 2017，Science，358，1570；Pian et al. 2017，Nature，551，67）。

你手上的金戒指，很可能是某次亿年前的双中子星并合在宇宙中锻造出来的。

电子分数 $Y_e$ 决定产物站在哪一座峰。动力学喷出物往往极端丰中子，走到第三峰（金、铂、铀）；盘风和中微子照射会抬高 $Y_e$，停在第一峰（锶、钇、锆）。同一次并合可以同时喷出两种成分，所以 AT2017gfo 才会先蓝后红。把银河系的金全部算在双中子星头上，还需要并合率和每次产量的乘积对得上；单次事件给的是存在性证明，不是全宇宙的账本。

极贫金属星大气里已经有铕和钍，说明宇宙很早就能做 r 过程。若双中子星并合有很长的延迟时标，早期产量可能不够，某些特殊超新星或磁旋转塌缩是否分担，仍是未关死的通道。

元素识别：锶的确认#

千新星光谱中具体识别哪些元素是一大挑战——高温、高速（喷出物速度约 $0.1$–$0.3\,c$）导致谱线严重展宽和混叠。Watson et al. (2019，Nature，574，497) 分析了 VLT/X-Shooter 的 AT2017gfo 光谱，在约 $350$–$800$ nm 处识别出锶（Sr，$Z=38$）的吸收特征，这是 r 过程第一峰元素的直接谱线证据。此后，分析还在红外谱中识别出可能的钡（Ba）和镧系元素特征。

锶被认出来，金线并没有被单独认出来。镧系原子结构太复杂，实验室跃迁数据不完整，高速展宽又把线混成一片。Watson 等人能指认 Sr II，是因为第一峰元素的离子谱相对干净，而且正好落在 X-Shooter 的光学窗口。光谱证实的是“发生了中子俘获并合成了第一峰”，红外光变证实的是“还需要高不透明度的重核”。两句话不是同一强度的证据，不能把锶的吸收线写成“看见了金”。

并合产物：最终形成了什么？#

GW170817 并合后形成的是黑洞还是超大质量中子星？这仍是未解之谜。从质量约束看：总质量约 $2.74\,M_\odot > 2.08\,M_\odot$（最重已知中子星），超过了大多数 EOS 的最大中子星质量，可能立即坍缩为黑洞，或先形成短暂旋转支持的超大质量中子星（hypermassive neutron star，HMNS），随后在约 $10$–$100$ ms 内坍缩。偏轴短 GRB 的产生需要一个稳定存在足够时间的中心引擎（至少持续约 $2$ 秒），支持了先形成 HMNS 后坍缩的图像。

余辉和短暴还要求中心引擎活过约两秒，这对“并合后立刻变成黑洞、再也喷不出东西”的图像不友好。更像是先形成超大质量中子星，靠转动撑一会儿，再塌缩。塌缩时刻写进状态方程：最大中子星质量越高，残骸能撑得越久。电磁和引力波在这里查的是同一颗星，只是一个看潮汐，一个看引擎寿命。

波形里的潮汐形变和余辉里的引擎寿命必须对得上同一套物态。若潮汐说物质偏软、电磁却要求残骸撑过两秒，那不是两种互相独立的成功测量，而是模型还没闭合。目前还没有把两端同时钉死的第二场事件。

哈勃常数的独立测量#

GW170817 还提供了测量宇宙膨胀率哈勃常数$H_0$ 的独立方法：引力波信号给出"绝对"距离，宿主星系 NGC 4993 的退行速度已知，两者之比即 $H_0$：

$$H_0 = 70^{+12}_{-8}\,\text{km/s/Mpc} \quad \text{(Abbott et al. 2017, Nature, 551, 85)}$$

这一结果介于 CMB 测量（$H_0 \approx 67.4$ km/s/Mpc，Planck Collaboration）和造父变星/超新星阶梯（$H_0 \approx 73$ km/s/Mpc，Riess et al.）之间，有望在未来多事件积累后对哈勃张力问题提供独立输入。

这是亮标准汽笛：引力波给绝对距离，宿主星系给红移，两者相除就是 $H_0$。单次事件的误差条从约 $62$ 跨到约 $82$，正好盖住 Planck 的低值和造父变星阶梯的高值，裁决不了哈勃张力。它的价值在于方法独立于宇宙距离阶梯，不经过造父变星和 Ia 型超新星的零点传递。原理见标准汽笛与引力波宇宙学。

后续用余辉几何改进倾角约束，可以把误差收窄，但在积累更多双中子星事件之前，这一条路仍是证明可行，而不是终审。

为什么这很重要#

GW170817 开创了多信使天文学的新纪元：同一事件被引力波和五个电磁波段覆盖。它同时实现了：确认短 GRB 来源、证实重元素 r 过程核合成场所、测量哈勃常数、约束中子星状态方程。单一事件对物理学、天文学和核物理的同时推进，在科学史上罕有先例。

一次事件同时服务核物理、高能天体物理和宇宙学，是因为它把无法在地面复制的密度、无法在实验室看到的核素、无法用梯子校准的距离，叠在同一个时间戳上。它也暴露了下一阶段的瓶颈：没有第二场同样完整的多信使并合，统计仍然是样本量为 1。时域巡天和引力波网络必须一起变得更快，锶之后的元素鉴定才有第二份光谱可做。

跨域连接#

• 引力波探测器LIGO与LISA：三台探测器联网把天区压到几十平方度，望远镜才可能在十几个小时内找到那个瞬变源。推论是：多信使的成败首先是定位问题——只有一台探测器时天区会大到无法搜完，所以"联合观测"的第一层含义是引力波网络内部的联合。
• 相对论多普勒：相对论喷流的辐射被强烈前向集束，视线稍偏就让观测亮度掉几个数量级。推论是："弱"不等于本征弱，同一次并合若正对我们就是一个普通短暴；集束效应还意味着大多数并合的伽马信号我们根本收不到，因此引力波才是无偏的触发器。
• 能量守恒：轨道能量以引力波形式流失，周期就必须按特定速率缩短，而第一个已知的双中子星系统实测的轨道衰减率正好对上这条预言。推论是：引力波在被直接听到之前几十年就已经被"称"过；两类证据查的是同一本能量账，只是一个盯着源、一个盯着波。
• 因果推断：单次事件无法重复，因果结论只能靠时间顺序与巧合概率——引力波、伽马暴与光学瞬变在同一小块天区按理论预言的先后出现，偶然重合的概率极低。推论是：这类单一事件推断的可信度不来自样本量，而来自预言的先验精确度，任一环节出现在预言之外整条链就得重估。
• 概率：能并合的双中子星必须连过好几道筛：第一次超新星的踢速不把双星拆散、公共包层把轨道收紧、第二次超新星同样不拆散。推论是：每道筛的通过率都是模型参数，所以理论并合率横跨好几个数量级；观测到的率反过来成了对这条演化链最有力的约束。

参考文献#

• Abbott, B.P. et al. (2017). GW170817: Observation of gravitational waves from a binary neutron star inspiral. Physical Review Letters, 119, 161101.
• Abbott, B.P. et al. (2017). A gravitational-wave standard siren measurement of the Hubble constant. Nature, 551, 85–88.
• Abbott, B.P. et al. (2017). Gravitational waves and gamma-rays from a binary neutron star merger: GW170817 and GRB 170817A. The Astrophysical Journal Letters, 848, L13.
• Goldstein, A. et al. (2017). An ordinary short gamma-ray burst with extraordinary implications. The Astrophysical Journal Letters, 848, L14.
• Coulter, D.A. et al. (2017). Swope Supernova Survey 2017a (SSS17a), the optical counterpart to a gravitational wave source. Science, 358, 1556–1558.
• Drout, M.R. et al. (2017). Light curves of the neutron star merger GW170817/SSS17a. Science, 358, 1570–1574.
• Kasen, D. et al. (2017). Origin of the heavy elements in binary neutron-star mergers from a gravitational-wave event. Nature, 551, 80–84.
• Watson, D. et al. (2019). Identification of strontium in the merger of two neutron stars. Nature, 574, 497–500.

延伸阅读#

• Metzger, B.D. (2019). Kilonovae. Living Reviews in Relativity, 23, 1. — 千新星理论与观测的权威综述（持续更新）
• Cowan, J.J. et al. (2021). Origin of the heaviest elements: The rapid neutron-capture process. Reviews of Modern Physics, 93, 015002. — r 过程核合成的系统综述
• Margutti, R. & Chornock, R. (2021). First multimessenger observations of a neutron star merger. Annual Review of Astronomy and Astrophysics, 59, 155–202.