跳转到内容
← 返回知识库
多信使天文学约 20 分钟阅读

kilonova与r过程核合成

千新星r过程核合成重元素双中子星并合中子俘获

概述

地球上所有比铁重的元素——银、金、铂、铀、钍,以及稀土元素(镧系)——都无法在恒星内部的正常核聚变中产生。它们需要极端的中子密度环境,在秒级时间内完成数十次中子俘获反应,这就是快中子俘获过程(rapid neutron capture process,r 过程)。r 过程在哪里发生?这是核天体物理学最重要的未解问题之一,持续争论了半个多世纪。双中子星并合(binary neutron star merger)和其产生的千新星(kilonova)如今被认为是主要场所——2017 年的 GW170817 事件提供了迄今最直接的观测证据。

常见误解是把 s 过程与 r 过程当成两套互不相干的元素表,或者把金、铂写成任何超新星炉里都会炼出的铁以上灰烬。真正的分界是时标:中子俘获若慢过 β 衰变,路径贴着稳定核素谷走,就是慢过程;若快过 β 衰变,路径被赶到丰中子一侧,中子流关掉后再衰变回来,才是快过程。金、铂、铀落在第三丰度峰附近,需要电子分数足够低的喷出物才能走到,并不是所有核心坍缩都能供应。光谱上目前被认出来的是锶这类第一峰核,金线本身还没有被单独指认。事件链与多信使时间戳见双中子星并合与千新星。

核合成的阶梯

氢燃烧到铁:核聚变的尽头

恒星在其一生中通过一系列核聚变反应产生能量,从氢聚变为氦,经氦聚变、碳燃烧、氧燃烧、硅燃烧,最终在铁(56^{56}Fe,$Z=26$,$A=56$)处达到核结合能的峰值——铁的每核子结合能(约 $8.8$ MeV)是所有核素中最高的(按核子数归一化的结合能曲线的顶峰)。铁核无法再进行放热聚变,核心坍缩超新星的机制正源于此。

因此,铁以上元素($Z > 26$)不能在正常核聚变中产生,需要通过中子俘获(neutron capture)在已有的种子核上添加中子,随后 β\beta 衰变增加质子数。

太初核合成只把氢、氦和极少量锂放到宇宙里,铁以上的稳定核素这条早期路径完全够不着、只能留给恒星。恒星内部的带电粒子反应在铁附近撞上库仑势垒和结合能峰值的双重墙,再往上加质子会吸热而不是放热。剩下的路只能靠不带电的中子往种子核上堆,再靠 β 衰变把中子改成质子,所以阶梯在铁处并没有断,只是从聚变换成了俘获。

慢中子俘获(s 过程)

s 过程(slow neutron capture)发生在 AGB 恒星(渐近巨星支)和大质量恒星的碳-氮-氧层中,中子密度约 108 cm−310^8\,\text{cm}^{-3},中子俘获时标(约 10310^3 年)远大于 β\beta 衰变时标,路径沿着稳定核素谷走。

s 过程产生钡(Ba)、铅(Pb)、铋(Bi)等约 50%50\% 的中等质量重元素(A∼90A \sim 90–$209$),但不能产生铀(235,238^{235,238}U)、钍(232^{232}Th)等超铀元素,也不能产生某些"中子丰富"的稳定核素(r 过程唯一能产生的核)。

渐近巨星支里的中子主要来自碳十三口袋的 13C(α,n)16O^{13}\mathrm{C}(\alpha,n)^{16}\mathrm{O},大质量恒星的氦燃烧壳层则可以靠 22Ne(α,n)25Mg^{22}\mathrm{Ne}(\alpha,n)^{25}\mathrm{Mg} 继续补中子。因为每俘获一次都要等 β 衰变,路径爬不到远离稳定线的超铀核,铀、钍因此几乎是快过程的独占产物。太阳系里钡、锶这类核往往是慢过程与快过程的混合物,不能把观测到的钡丰度整段记在慢过程账上,更不能把金写成慢过程沿途的副产品。

快中子俘获(r 过程)

r 过程需要中子密度 nn>1020n_n > 10^{20}–1024 cm−310^{24}\,\text{cm}^{-3}(比 s 过程高 $12$–$16$ 个数量级)、温度约 10910^9–101010^{10} K、时间尺度约 $0.1$–$10$ 秒。

在这种条件下,种子核(通常是铁族)以比 β\beta 衰变快得多的速率俘获中子,沿着远离稳定核素谷的"丰中子"路径快速爬升,形成高度不稳定的重核。当中子流停止后,这些不稳定核通过一系列 β\beta 衰变回到稳定核素谷,形成最终的重元素分布。r 过程的"峰"对应于幻数(neutron magic numbers,$N = 50, 82, 126$)处的特别稳定壳层,在最终分布上形成约 $A = 80$(锶、钇、锆)、$A = 130$(碲、碘、氙)和 $A = 195$(铂、金、汞)处的丰度峰。

中子幻数处的等待点把物质堆在特定质量数,冻结后再 β 衰变,这才在稳定谷上留下这三座峰。第三峰大约在质量数 195 附近,铂、金、汞都坐在这一带,所以要问金从哪来,先要问喷出物有没有够低的电子分数、够长的中子照射,好让种子核跨过第二峰走到第三峰。裂变循环还可能把更重的核打回来并改写锕系与第三峰的相对高度,这一段核数据仍然不完整,不能写成已经结算的元素配比。

千新星的物理图像

千新星(kilonova,也称 macronova)是双中子星(或黑洞-中子星)并合喷出物的光学/红外瞬变,最早由李、帕新斯基(Li & Paczynski 1998)和麦茨格(Metzger et al. 2010)在理论上预言。

喷出物不是从中子星表面“蒸发”出来的均匀气体,而是潮汐、冲击和盘风几路物质叠在一起,电子分数可以差出一截。中微子照射把中子转成质子,等于当场改写后续核合成路径,照射越强就越难让种子核跨过第二峰走到金所在的第三峰。所以金、铂的产量对并合几何和残骸寿命敏感,不能从一次事件的总喷出质量直接按太阳系比例去拆。

并合与喷出物

双中子星并合产生约 $0.01$–0.1 M⊙0.1\,M_\odot 的高速喷出物,主要有以下成分:

动力学喷出物(dynamical ejecta,∼0.001\sim 0.001–0.01 M⊙0.01\,M_\odot):

  • 由潮汐力剥离(tidal tail,较冷,极端中子丰富,Ye<0.1Y_e < 0.1)
  • 冲击加热喷出(shock-heated,在碰撞界面,YeY_e 较高,约 $0.2$–$0.4$)

星风喷出物(wind ejecta,∼0.01\sim 0.01–0.05 M⊙0.05\,M_\odot):

  • 中微子风(neutrino-driven wind,来自中心超大质量中子星或黑洞吸积盘,Ye∼0.25Y_e \sim 0.25–$0.45$)
  • 磁流体动力学风(MHD wind)和粘性驱动盘风

YeY_e(电子分数,即质子数与总重子数之比)是决定 r 过程产物的关键参数:

  • Ye<0.25Y_e < 0.25:极端丰中子,重 r 过程(镧系元素、锕系元素),κ∼10\kappa \sim 10–100 cm2/g100\,\text{cm}^2/\text{g}(高不透明度,红色千新星)
  • Ye>0.25Y_e > 0.25:轻 r 过程(锶、钡、银等),κ∼1\kappa \sim 1–5 cm2/g5\,\text{cm}^2/\text{g}(低不透明度,蓝色千新星)

潮汐尾相对冷、几乎来不及被中微子照射,电子分数可以低到让反应网络跑完重 r 过程,金和锕系主要指望这一支。碰撞界面和盘风会被中微子把中子转成质子,电子分数抬高后往往停在第一峰或第二峰,锶、钇可以很多,金却不必跟着等比例出现。同一场并合因此可以同时喷出“能做金”和“做不出金”的气体;把银河系的金全部记在平均每次并合头上,还缺少并合率与每事件产量的乘积,不能用一次光变去结算宇宙化学。Freiburghaus、Rosswog 与 Thielemann 在 1999 年用流体加核反应网证明:电子分数大约 $0.1$ 的中子星物质减压膨胀,就能在质量数 $130$ 及以上复现太阳系 r 过程图案——这是并合作为场所的理论底气,不是一张已经点名的金元素清单。

千新星的光变曲线

不同不透明度的喷出物成分导致千新星呈现出双成分光变曲线,这在 AT2017gfo(GW170817 的千新星对应体)中得到了完美体现:

蓝色成分(blue kilonova,t≲3t \lesssim 3 天):κ∼1\kappa \sim 1 cm2^2/g,喷出速度约 $0.2c$–$0.3c$,峰值在近紫外/蓝色波段(λ∼3000\lambda \sim 3000–$5000$ Å),对应中子丰度适中(Ye>0.25Y_e > 0.25)的轻 r 过程元素。

红色成分(red kilonova,t∼3t \sim 3–$10$ 天):κ∼10\kappa \sim 10–$100$ cm2^2/g(镧系元素),喷出速度约 $0.1c$,峰值在近红外(λ∼10000\lambda \sim 10000–$20000$ Å),对应极端丰中子(Ye<0.2Y_e < 0.2)的重 r 过程元素。峰值时间 tpeak∼κMej/(4πvc)t_{peak} \sim \sqrt{\kappa M_{ej}/(4\pi v c)},对于镧系元素高不透明度,峰值被推迟到约 $5$–$7$ 天。

峰值光度约 Lpeak∼1041L_{peak} \sim 10^{41}–104210^{42} erg/s,是经典新星的约 10310^3 倍(故名"千新星"),但比 Ia 型超新星约弱 $10$–$100$ 倍。

Barnes 与 Kasen 在 2013 年指出,只要喷出物里有镧系,不透明度会比铁族高出几个数量级,光变会变红、变暗、拖到大约一周,峰值被推到近红外。红千新星因此是“发生了重 r 过程、镧系在场”的辐射转移证据,仍然不是金原子被单独指认出来的指纹谱。把红外变亮写成“看见了金”,等于把不透明度约束和谱线证认当成同一件事;能源是新合成不稳定核的衰变热,不是镍-56 那条超新星衰变链。

AT2017gfo 的元素鉴定

整体谱特征

AT2017gfo 的早期光谱($t < 2$ 天)显示蓝色连续谱,有宽吸收特征(宽度对应速度约 $0.2c–0.3c$),说明光球在高速膨胀。晚期($t > 5$ 天)光谱红化,近红外增亮,符合高不透明度 r 过程元素的特征。

锶(Sr)的识别

Watson et al. (2019,Nature,574,497) 通过与原子数据比较,在约 $400$–$800$ nm 的早期光谱中识别出Sr II 的特征吸收(约 $800$–$900$ nm 处的宽吸收,多普勒蓝移对应速度约 $0.27c$)。锶($Z=38$,A≈88A \approx 88–$90$)是 r 过程第一峰元素,其存在证实了中子俘获过程在并合喷出物中确实发生。

其他元素的可能证据

Pian et al. (2017,Nature,551,67) 分析了约 ∼1\sim 1–$7$ 天的宽带光谱,拟合出需要镧系元素(Nd、Ce、La 等,$Z = 57$–$71$)的存在来解释红外光度。Smartt et al. (2017,Nature,551,75) 则提出铯(Cs,$Z=55$)可能的识别。完整的元素鉴定仍在进行中,主要障碍是镧系元素的原子结构复杂、实验室测量的跃迁数据不完整,以及多元素混叠使谱线无法单独分辨。

金、铂的开壳层和镧系一样复杂,实验室波长数据更稀,高速膨胀又把线混成一片,所以目前没有人在 AT2017gfo 里单独指认出金的吸收。能说的是:红成分要求镧系量级的不透明度,而核反应网在足够丰中子的喷出物里会同时做出第三峰,金是这条链上的预期产物,不是已经点名的谱线。第一峰被看见、第三峰被辐射转移间接要求,两句话的证据强度不一样,不能写成同一级别的发现。

r 过程的宇宙历史

银河系中重元素的化学演化给出了 r 过程宇宙历史的线索:

早期银河系问题:在银河系中金属贫乏([Fe/H]<−2\text{[Fe/H]} < -2)的古老星中,r 过程元素(如钡、铕 Eu)丰度有极大散射,暗示 r 过程来源是"稀少且产量大"的事件,而非普遍分布的超新星。双中子星并合率低但每次产量高,符合这一约束(Côté et al. 2019,The Astrophysical Journal,875,106)。

时间延迟问题:双中子星形成需要约 10810^8–101010^{10} 年(轨道衰减时标),这意味着 r 过程元素注入星际介质有相对于恒星形成的时间延迟,可能与某些金属贫化 r 过程元素丰度模式不符。这促使一些学者保留核心坍缩超新星("坍缩星",collapsar)作为额外的 r 过程场所(Siegel et al. 2019,Nature,569,241)。

延迟时标来自双星演化,不是核反应本身的快慢。第一颗超新星若把系统踢散,就永远等不到并合;公共包层若没把轨道收够紧,引力波也来不及在宇宙年龄内把两颗中子星拖到一起。所以早期银河系里已经出现的铕、钍,对“并合是唯一工厂”构成压力,并不否定 2017 年那一次事件里确实跑过 r 过程。

目前更稳妥的说法是:双中子星并合已被观测证明能够运行 r 过程,并且很可能供应了银河系中相当一部分重 r 过程核(质量数大于约 $130$ 的金、铂族等),份额却没有被钉成一个精确百分比。核心坍缩超新星、尤其是快速旋转的坍缩星,仍可能分担早期宇宙或轻 r 过程的产量,场所清单并没有在 2017 年一次性关账。Côté 等人 2019 年的化学演化计算强调:只靠并合,未必能同时满足早期散射与总产量,额外通道是否必要,仍取决于延迟时标和每次喷出物质量。

为什么这很重要

r 过程核合成是物理学最美丽的问题之一:它把核物理、中子星物理、天文观测和宇宙化学演化连接在一起。你今天佩戴的黄金,极大概率是某次亿年前的双中子星并合在宇宙中锻造出来的,随后被超新星遗迹的激波卷入星际气体云,参与太阳系的形成,最终在地球上以矿石形式被人类挖掘。这条从中子星到手指上的戒指的链条,是宇宙物质循环最动人的叙事之一。

地球上的金是第三峰 r 过程核,地质循环只是把它从地幔送到矿脉,并不能在地壳里靠化学反应新造。要解释这块金,必须同时回答核路径有没有走到 $N=126$ 等待点,以及天体喷出物够不够丰中子,缺一则只剩类比。把戒指叙事写成已经完成的元素鉴定,会把“发生了快中子俘获”说成“望远镜已经看见了金线”。

跨域连接

  • 恒星核合成:慢中子俘获沿稳定核素谷爬行,产不出铀钍这类超铀核,也产不出某些丰中子的稳定核素——这些"只有快过程能给"的核素就是判据。推论是:要证明快过程发生在哪里,不必看到全过程,只要在喷出物里认出这类核素的谱线即可。
  • 核壳层模型:最终丰度的几个峰落在中子幻数处,因为壳层闭合的核特别稳定、再俘获一个中子更难,于是丰度在那里堆积。推论是:峰的位置由核结构定死,与在哪里发生无关,所以峰位是"是否经历过快中子俘获"的判据,而峰的相对高度才反映环境的中子丰度。
  • 原子结构:认出锶靠的是它的特征吸收线,静止波长由原子能级定死,观测到的偏移则给出接近光速三成的膨胀速度。推论是:瓶颈在原子数据而不在望远镜——镧系元素能级极复杂、实验室跃迁数据不全、多元素谱线互相混叠,所以完整的元素清单至今没有列完。
  • 电子轨道与量子数:镧系元素的电子结构复杂、能级密集,形成成片的吸收,不透明度比轻元素高一两个数量级。推论是:同一次并合因此给出两段光变——低不透明度的轻元素成分几天内变蓝又变暗,高不透明度的镧系成分把峰值推迟到五到七天并落在近红外;颜色本身就是成分的读数。
  • 冰芯古气候:读银河系的化学演化和读冰芯是同一套办法:古老的贫金属恒星是不同时刻星际介质成分的封存样本,按金属丰度排序就得到一条时间序列。推论是:这些星里重元素丰度散射极大,指向"稀少但每次产量大"的来源;而并合还要等轨道衰减,注入有时间延迟,是否还需要另一类场所目前没有定论。

参考文献

  • Metzger, B.D. et al. (2010). Electromagnetic counterparts of compact object mergers powered by the radioactive decay of r-process nucleosynthesis products. Monthly Notices of the Royal Astronomical Society, 406, 2650–2662.
  • Watson, D. et al. (2019). Identification of strontium in the merger of two neutron stars. Nature, 574, 497–500.
  • Pian, E. et al. (2017). Spectroscopic identification of r-process nucleosynthesis in a double neutron-star merger. Nature, 551, 67–70.
  • Cowan, J.J. et al. (2021). Origin of the heaviest elements: The rapid neutron-capture process. Reviews of Modern Physics, 93, 015002.
  • Côté, B. et al. (2019). Neutron star mergers might not be the only source of r-process elements in the Milky Way. The Astrophysical Journal, 875, 106.
  • Freiburghaus, C., Rosswog, S. & Thielemann, F.-K. (1999). r-Process in Neutron Star Mergers. The Astrophysical Journal Letters, 525, L121–L124.
  • Barnes, J. & Kasen, D. (2013). Effect of a High Opacity on the Light Curves of Radioactively Powered Transients from Compact Object Mergers. The Astrophysical Journal, 775, 18.

延伸阅读

  • Metzger, B.D. (2019). Kilonovae. Living Reviews in Relativity, 23, 1. — 千新星物理的权威综述(持续更新)
  • Burbidge, E.M. et al. (1957). Synthesis of the elements in stars. Reviews of Modern Physics, 29, 547–650. — 核合成理论的奠基论文(B²FH 论文)
  • Lattimer, J.M. & Schramm, D.N. (1974). Black-hole–neutron-star collisions. The Astrophysical Journal Letters, 192, L145. — 早期预言中子星并合产生 r 过程元素的论文