跳转到内容

kilonova与r过程核合成

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

概述

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

核合成的阶梯

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

恒星在其一生中通过一系列核聚变反应产生能量，从氢聚变为氦，经氦聚变、碳燃烧、氧燃烧、硅燃烧，最终在铁（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 过程唯一能产生的核）。

快中子俘获（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$（铂、金、汞）处的丰度峰。

千新星的物理图像

千新星（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}（低不透明度，蓝色千新星）

千新星的光变曲线

不同不透明度的喷出物成分导致千新星呈现出双成分光变曲线，这在 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$ 倍。

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$）可能的识别。完整的元素鉴定仍在进行中，主要障碍是镧系元素的原子结构复杂、实验室测量的跃迁数据不完整，以及多元素混叠使谱线无法单独分辨。

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）。

目前学界共识是：双中子星并合是宇宙中大部分（可能约 ∼80%\sim 80\%）重 r 过程元素（$A > 130$，金、铂族等）的主要来源，核心坍缩超新星（collapsar）可能提供额外贡献，特别是对轻 r 过程元素（A∼80A \sim 80–$130$）。

为什么这很重要

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

跨域连接

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

参考文献

  • 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.

延伸阅读

  • 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 过程元素的论文

正文中引用本条 · 2 · 其中 1 条来自其他领域

化学核化学与放射化学