r过程核合成与中子星并合
2017年8月17日，LIGO和Virgo探测到了一次引力波信号。信号来自1.3亿光年外两个中子星的碰撞。1.7秒后，费米伽马射线空间望远镜在同一个方向看到了一束短伽马射线暴。接下来几天，全球70多台望远镜把镜头对准了那片天区，看到一个先蓝后红的闪光——蓝的持续两天，红的持续一周多。
这个闪光叫千新星（kilonova）。它是中子星碰撞后喷出的物质在膨胀冷却过程中，由r过程新合成的重元素放射性衰变所发的光。2017年10月，多个团队在Nature和Science上同时发文，从千新星的光谱中辨认出了锶——这是人类首次直接观测到r过程正在宇宙中发生。
七十年前，Burbidge、Burbidge、Fowler和Hoyle在一篇论文里提出了r过程的概念。七十年后，我们终于亲眼看到了它在哪儿发生。
但看到不代表看懂了。r过程的完整图景——它在哪些天体环境中发生、每个环境贡献了多少、核物理输入量有多可靠——到今天仍然是核天体物理最大的未解问题之一。GW170817证实了中子星并合是一个r过程场所，但它是不是唯一的？甚至是不是主要的？答案正在变得越来越微妙。
---
先说清楚r过程是什么。
宇宙中的元素不是一次造出来的。大爆炸造了氢和氦。恒星内部的热核聚变造了从碳到铁的元素——一直到铁-56，聚变还释放能量，再往上就吸热了。比铁更重的元素需要别的机制，因为聚变不再自发进行。
重元素主要通过两条路径产生：慢中子俘获过程（s过程）和快中子俘获过程（r过程）。s过程发生在红巨星内部，中子密度低，一个原子核俘获一个中子后，有时间先β衰变再俘获下一个。r过程发生在中子密度极高的爆发性环境中，原子核来不及β衰变就被下一个中子砸中——快到中子俘获的时间尺度比β衰变短几个数量级。
r过程的后果是：原子核沿着中子滴线一路狂奔，把中子壳层一层一层填满，直到中子俘获和光致分裂达到平衡。然后中子源消失（"冻结"），这些极度中子富集的原子核通过一系列β衰变回到稳定线，变成你能在地球上找到的元素。
r过程产物的标志性特征是三个丰度峰：第一峰在A≈80（硒、溴、氪附近），第二峰在A≈130（碲、碘、氙附近），第三峰在A≈195（锇、铱、铂附近）。太阳系中约一半的比铁重的元素——包括你手指上的金戒指、手机里的稀土元素、核电站的铀燃料——都是r过程的产物。
问题是：r过程到底在哪里发生？
---
从Burbidge等人1957年提出r过程起，候选场所一直在变。早期最流行的是核心坍缩超新星里的中微子驱动风。
中微子驱动风的图像是这样的：大质量恒星核心坍缩形成原初中子星后，新生的中子星表面温度极高，大约10¹¹开尔文，向外辐射大量中微子。这些中微子被质子和中子弱吸收，驱动一股向外的物质流——中微子驱动风。风里的物质起初是质子和中子的混合物，膨胀冷却后发生α富集冻结，形成铁峰核素作为种子。如果风的中子数足够高（电子分数Y_e ≲ 0.3），后续的中子俘获就能推进到r过程的第二甚至第三峰。
1992年，Woosley和Hoffman计算了中微子驱动风中的r过程。2001年，Thompson、Burrows和Pinto在ApJ上给出了更细致的模拟。结果一度让人乐观——如果原初中子星足够紧致，风足够快，条件确实可以支持r过程。
但2000年代后期，问题出现了。更精确的2D和3D模拟发现，中微子驱动风的中子富集度不够。原因是中微子与物质的相互作用把质子变中子、中子变质子的速率差不多，结果Y_e被推向0.4-0.5附近——对r过程来说太"质子化"了。你得到的是第一峰元素（锶、钇、锆），但到不了第二峰，更别说第三峰。
2014年，Martínez-Pinedo等人在J.Phys.G上总结了这个问题：标准的中微子驱动风最多产生r过程的第一峰，除非引入极端的磁场或旋转——但那就不是"标准"超新星了。
所以核心坍缩超新星的"标准"中微子驱动风，目前被认为是轻r过程（第一峰）的可能贡献者，但不是重r过程（第二、第三峰）的主要场所。这个结论对很多人来说是个失望——超新星每秒钟在宇宙中发生好几次，而中子星并合几万年才一次。如果重r过程只能靠并合，那宇宙早期的重元素从哪来？
---
2017年GW170817之前，中子星并合作为r过程场所已经有很长的理论历史。Lattimer和Schramm在1974年就提出了中子星双星系统并合可以喷射r过程物质的想法。1970年代末到1980年代，Symbalisty、Eichler、Meyer等人做了早期计算。但那时连中子星双星系统的存在都还没被确认——Hulse-Taylor脉冲星1974年才发现——所以这些想法停留在理论层面。
GW170817改变了一切。
引力波信号告诉了我们碰撞前的系统参数：两个中子星各约1.36和1.17倍太阳质量，距离1.3亿光年。短伽马射线暴GRB 170817A在引力波信号后1.7秒到达——这是第一次把引力波和电磁信号关联起来。
但真正让核天体物理学家兴奋的是千新星AT2017gfo。
千新星的亮度随时间的演化给出了r过程产物的大致图像：早期（1-2天）的蓝色成分来自高Y_e（电子分数高，中子相对少）喷出物中较轻r过程元素（锶、钇、锆）的衰变；后期（3-10天）的红色成分来自低Y_e喷出物中镧系元素的衰变——镧系元素有大量密集的跃迁线，把蓝光全部吸收，只透红光和红外。
2017年10月，Pian等人和Shappee等人分别在Nature和Science上发表了千新星的早期光谱和光变曲线。同一天，Smartt等人也在Nature上报告了光谱中810纳米处的吸收特征，后来被确认来自锶——r过程路径上第一个被光谱证认的元素。
这个确认不是一蹴而就的。2019年，Watson等人在Nature上报告了AT2017gfo光谱中1微米附近的发射线来自锶的禁戒跃迁[Sr II]。这个证认后来经过了多次检验和争议。2026年4月，Aouad等人在arXiv上发表了新的非LTE辐射转移建模，确认了锶和氦的存在——锶线来自r过程冻结后的衰变产物，氦则来自α富集冻结的残留。这是对GW170817中r过程光谱信号的最新确认。
---
中子星并合里的r过程物质是怎么被抛出来的？答案是三种不同的喷出机制，对应三种不同的中子富集度。
第一种是动力学喷射物（dynamical ejecta）。两颗中子星碰撞时，潮汐力把物质从碰撞面的两侧撕扯出来，形成两片"潮汐尾"。同时，碰撞接触面的高温高压也会把一部分物质垂直挤压出去。动力学喷射物在碰撞后几毫秒内就被弹射，速度约0.1-0.3倍光速。它的Y_e极低（0.05-0.2），因为物质直接来自中子星的中子富集外壳。这是最理想的r过程场所——中子太多，种子核太少，物质一路跑到第三峰甚至裂变循环区。
2023年，Just等人在ApJ上做了系统的双中子星并合喷出物核合成计算，发现不对称双星（质量差异大的系统）的动力学喷射物更中子富集，因而产出更重的r过程元素。这跟直觉相符——不对称碰撞的潮汐力更强，撕扯出更多深层中子物质。
第二种是盘风（disk wind）。碰撞后，没有立刻被弹射的物质围绕残余天体形成吸积盘。盘的内侧高温高密，外侧逐渐冷却。盘中的物质通过粘滞加热、中微子加热和磁场重新连接被慢慢推出——这些就是盘风。盘风的时标从几十毫秒到几秒，速度约0.05-0.1倍光速，Y_e取决于中微子照射的强度：如果中微子把足够多的中子变成了质子，Y_e会升到0.3-0.4，r过程只到第二峰；如果中微子照射弱，Y_e可以低到0.15-0.25，r过程推进到第三峰。
盘风是目前r过程核合成模拟中最不确定的部分。2025年，广西大学邱母伍子莫等人使用不同状态方程的天体物理条件作为核合成输入，发现状态方程的选择对r过程产物和千新星发射有显著影响——更"软"的状态方程（中子星更紧致）导致更小的动力学喷射物和更热的盘，盘风Y_e更高，r过程不完整。
第三种是粘滞/磁驱动喷出物。如果并合后的残余天体是超大质量中子星（hypermassive neutron star, HMNS）而不是立刻坍缩成黑洞，它会在差分旋转时通过磁重新连接和磁流体湍流持续喷出物质。2022年，Collins等人用3D广义相对论磁流体力学模拟研究了HMNS的r过程核合成，发现磁驱动喷出物的Y_e分布很宽（0.1-0.4），可以同时产生从第一峰到第三峰的元素。
2025年10月，Nedora等人用BAM码做了最新的广义相对论辐射磁流体力学模拟，系统研究了中子星自旋对喷出物和核合成的影响。他们发现：快速自旋的双星系统产生更多的喷出物，盘风更富中子，r过程产物更重。这是第一次在完全GRMHD模拟中定量评估自旋对r过程的影响。
2025年5月，PRL上发表了一篇里程碑式的论文：日本团队完成了最长持续时间的双中子星并合数值相对论中微子辐射磁流体力学模拟，跟踪到碰撞后约1.5秒——此前大多数模拟只跟踪到0.1-0.3秒。他们在这次模拟中首次看到了喷流从坍缩的残余中子星中形成并穿透喷出物的过程，直接连接了短伽马射线暴和r过程核合成的物理图像。
---
GW170817的千新星给了我们一个具体的r过程样本。但它到底产出了多少什么元素？
这个问题比看起来难得多。千新星的辐射转移是出了名的难算——镧系和锕系元素有几十万个原子能级，大部分没有实验测量的跃迁数据，只能靠理论估计。而且喷出物不是均匀的一团，不同层的Y_e、密度、速度都不一样，光子在不同层之间散射、吸收、再发射，最后逃逸出来的光谱是所有层叠加的结果。
2023年，Vieira等人开发了SPARK（Spectroscopic r-Process Abundance Retrieval for Kilonovae）工具，用GW170817在1.4天的光谱做了第一次千新星元素丰度反演。2024年2月，SPARK II在ApJ上发表，分析了1.4、2.4和3.4天的光谱数据，反演了多组分喷出物的元素丰度模式——确认了镧系元素的存在。
2025年4月，SPARK III在arXiv上发表（2026年1月正式发表于ApJ），首次自洽地将光谱反演的丰度模式与光变曲线拟合。结果是：GW170817的喷出物总质量约0.04-0.06倍太阳质量，其中镧系质量分数约0.5-2%——这跟理论预期的"高镧系"喷出物一致，但镧系分数比很多极贫金属星的观测值要低一个量级。
2024年10月，Gillanders等人在A&A上汇编了AT2017gfo从早期到晚期的详细光谱序列，追踪了r过程光谱特征随时间的逐小时浮现和演化——这些数据是理解千新星辐射转移和元素证认的基准。
2025年底，Leiding等人在INSPIRE-HEP上报告了对AT2017gfo镧系元素约束的改进，使用更新的原子数据将镧系质量分数的置信区间进一步收窄。2026年3月，PRC上发表了经过校准的镧系原子数据用于千新星辐射转移计算——这是该方向的一个重要基础设施进展。
---
如果中子星并合是r过程的主要场所，那它能不能解释银河系里所有的r过程元素？
这是一个叫"银河系化学演化"的问题。基本逻辑是：把中子星并合的产额、发生率、延迟时间分布输入到银河系化学演化模型里，看能不能拟合观测到的元素丰度。
2024年2月，Choplin等人发表在MNRAS上的论文做了系统计算，结论是：双中子星并合可以充分解释银河系重r过程元素——前提是每次事件产出约0.03倍太阳质量的r过程物质，发生率约每10万年一次。这个数字跟GW170817的观测一致。
但有一个麻烦：极贫金属星。
极贫金属星是银河系里最古老的恒星，金属丰度[Fe/H]低于-3（铁含量不到太阳的千分之一）。它们形成于宇宙早期，只经历过少数几次核合成事件的污染。关键观测是：在这些极贫金属星中，[Eu/Fe]——铕与铁的比值——展现出巨大的星间离散，高达几个数量级。这说明r过程事件在早期宇宙极为稀少，单次事件的产额必须很大，才能解释某些恒星中铕含量比铁还高的现象。
中子星并合的延迟时间是一个问题。两个中子星从双星系统形成到最终碰撞，可能需要几千万到几十亿年——取决于初始轨道间距。如果延迟太长，宇宙早期还没有足够多的并合发生，就无法解释极贫金属星中的r过程富集。
2025年3月，ApJ上一篇论文利用18颗已知的再循环毫秒脉冲星双中子星系统，估计了中子星并合的延迟时间分布，发现虽然大部分并合确实延迟较长，但存在一个"超快"通道——某些紧致双星可以在不到一亿年内碰撞。Beniamini和Piran在2024年的工作中估计，约0.1-1%的双中子星系统可以走超快通道。这个比例虽小，但在宇宙早期可能足够了。
2025年7月，A&A上发表了"Europium, we have a problem"——直接点出了当前模型的困境：银河系化学演化模型在拟合铕丰度时，无论怎么调参，总有些极贫金属星的铕含量超出预测。作者提出，可能需要一个额外的r过程场所来补充早期宇宙的重元素——这就回到了超新星的讨论。
2025年11月，arXiv上发表了用多目标银河系化学演化约束r过程核合成的工作，系统搜索了参数空间，发现任何单一场所都很难同时拟合所有观测约束。
---
如果中子星并合不是唯一的r过程场所，那还有什么？
2026年3月，以色列理工学院的Noam Soker在arXiv上发表了一篇系统比较七种r过程场景的论文。这是目前唯一涵盖所有主要场景的系统性对比。他的结论很清晰：早期宇宙的r过程需要至少两个主要场所。
第一个是磁转动超新星（magnetorotational supernova, MR-SN）。这是一类极端的核心坍缩超新星——前身星快速旋转，拥有强磁场（约10¹²高斯以上），磁转矩驱动喷流状爆炸。2023年，Reichert和Obergaulinger用3D中微子-MHD超新星模型做了后处理核合成，发现强磁场模型确实可以产出第三峰元素（A≈195），但产额约0.003-0.03倍太阳质量，跟中子星并合比差一个量级。
2024年3月，Reichert等人在arXiv上进一步研究了磁转动超新星最内层喷出物的核合成。2025年4月，另一团队在HAL上发表了五个3D MHD模型对核合成影响的系统比较。关键发现是：磁转动超新星的r过程产物与铁强关联——r过程元素跟铁一起被抛出来。这跟观测中"[Eu/Fe]与Fe相关"的模式一致，说明磁转动超新星可能是轻r过程（"Eu+Fe型"）的主要贡献者。
第二个是共有包层喷射流超新星（CEJSNe），由中国科学家提出。2024年8月，中科院近物所金仕纶研究员及合作者在ApJ上发表了这种新机制。物理图像是这样的：一个双星系统，一颗已经变成中子星，另一颗还是红超巨星。红超巨星膨胀后吞噬了中子星，形成共有包层。中子星在红超巨星的核心里螺旋下沉，剧烈吸积物质，最终向两极发射喷流。喷流物质高温高密又快速冷却，正是r过程需要的条件。
CEJSNe的独特之处在于：它产出的超镧元素（比镧更重的元素，包括第三峰附近的锇、铱、铂）比任何其他场景都多——Ir/Eu比值比坍缩星高一个量级，比磁转动超新星高三个量级。而且CEJSNe不与铁伴生——喷射流从恒星核心内部喷出，没有经过外层的铁聚变区。这正好解释了观测中第三峰元素与铁"解耦"的现象。
Soker的论文把七种场景排了一个表：
| 场景 | 单次r过程产额 | 是否与Fe伴生 | 轻r过程 | 第三峰r过程 |
|---|---|---|---|---|
| 中子星-中子星并合 | ≈0.03 M☉ | 否（延迟长） | 不与Fe关联 | 仅超快并合可能 |
| 磁转动超新星 | ≈0.003-0.03 M☉ | 是 | 可能的主要贡献者 | 与Fe强关联，不符观测 |
| CEJSNe | ≈0.01-0.05 M☉ | 否 | 仅极大质量氧核心时可能 | 可能的主要贡献者 |
| 坍缩星 | ≈0.08-0.3 M☉ | 是 | 条件苛刻 | 与Fe关联，不符观测 |
| 中子星/磁星风 | ≈10⁻⁴ M☉ | 是 | 产额太低 | 产额太低 |
| 磁星巨耀发 | ≈10⁻⁵-10⁻³ M☉ | 是 | 产额太低 | 产额太低 |
| 白矮星吸积诱导坍缩 | ≈10⁻²-0.1 M☉ | 否 | 延迟太长 | 延迟太长 |
被排除的主要贡献者：中子星风、磁星巨耀发（产额太低）、白矮星坍缩（延迟太长）。
剩下的图景是：磁转动超新星主导轻r过程（与Fe关联），CEJSNe主导第三峰（与Fe解耦），中子星并合在宇宙中后期逐渐变成重要但不是唯一的贡献者。这个框架还在被验证，但已经开始改变人们对r过程图景的理解。
---
黑洞-中子星并合是另一个被关注的r过程场所。
如果一个黑洞和一个中子星碰撞，潮汐力会在中子星被吞噬之前撕下部分物质。这些物质形成吸积盘，然后像双中子星并合的盘风一样被推出。2022年底，Just等人在arXiv上首次用全GR中微子辐射MHD模拟研究了黑洞-中子星并合喷出物的核合成。
2024年12月，一个团队在PubMed上报告了黑洞-中子星并合喷出物中"锕系增强"r过程的发现——某些低Y_e的盘风物质可以产出大量锕系元素（钍、铀），导致锕系丰度异常高。这跟天文学家在部分r过程增强星中观测到的"锕系增强"现象可能有关联。
但黑洞-中子星并合有一个根本限制：黑洞太大了，可能把整个中子星吞掉，什么都不吐出来。只有当黑洞质量不太大（约5-10倍太阳质量）且自旋较快时，潮汐力才能在视界之外把物质撕下来。LIGO的观测表明，大多数恒星级黑洞的质量在10倍太阳质量以上，对r过程不太友好。
---
r过程的核物理输入量是一个独立于天体物理的巨大问题。
r过程路径经过的区域远离稳定线——那里有大约5000个未知的原子核。你需要知道它们的质量、β衰变半衰期、中子俘获截面、裂变概率和裂变产物分布。这些量绝大部分无法实验测量，只能靠理论模型算，而不同模型之间差异巨大。
核质量是最基本的输入。2025年1月，A&A上发表了核质量模型不确定性对r过程核合成影响的研究，发现不同质量模型给出的第三峰位置可以偏移5-10个质量数——这在太阳系丰度拟合中是不可接受的。2024年11月，PLB上报告了FRIB在⁷⁸Ni附近产生极丰中子核素的实验进展，为r过程路径上的质量模型提供了新的校准点。2024年3月，PRC上发表了用飞行时间-磁刚度法测量丰中子核质量的新技术——这是第一个适用于r过程路径上短寿命核的质量测量方法。
β衰变半衰期决定了r过程的节奏——衰变越快，中子俘获路径越靠近稳定线，产出越轻；衰变越慢，路径越远离稳定线，产出越重。2025年3月，PRC上发表了用机器学习预测β衰变半衰期并应用于r过程计算的工作，发现ML模型在某些区域比传统QRPA计算更接近实验值。2025年11月，arXiv上发表了大规模β衰变率计算及其对r过程影响的研究，系统比较了五种全球β衰变模型，发现不同模型给出的冻结时间可以差几倍，直接影响r过程的最终丰度分布。
裂变是另一个关键。当r过程路径推到A≈280附近时，原子核变得极度不稳定，自发裂变和中子诱导裂变把重核打碎，碎片回到A≈130附近重新开始俘获——这就是"裂变循环"。裂变循环的存在意味着第三峰的丰度不取决于r过程能推多远，而取决于裂变产物回到哪里。但裂变势垒高度的理论不确定性高达几MeV——这对裂变概率是数量级的差异。2023年，Chinese Physics Letters上发表了用微观自洽方法计算r过程终点附近²⁷⁸Cf和²⁸⁶Cf裂变性质的工作，但这些计算离给出全球裂变率还有很远的距离。
FRIB正在逐步缩小这些不确定性。作为目前世界上最强的放射性核束装置，FRIB可以产生已知约5000种同位素的80%。2024年9月，Taylor&Francis上发表了FRIB在核天体物理中角色的综述——r过程路径上的核质量、β衰变率和中子俘获截面的实验测量是FRIB最核心的科学目标之一。但这些实验需要时间——每一个r过程核素都需要单独的束流时间和专门的探测方案。
---
中国在这场研究中的角色正在快速增长。
金仕纶团队2024年提出的CEJSNe r过程机制是一个原创性贡献——在此之前，没有人系统提出共有包层喷射流超新星作为r过程场所。CEJSNe的图像非常直观：中子星钻进大质量恒星的核心，像搅拌机一样搅动核心物质，然后从两极喷射r过程产物。这个场景虽然发生率比普通超新星低（需要特殊的双星演化路径），但单次产额足够大，可以解释极贫金属星中最高的铕含量。
北京大学孟杰团队的核物理与核天体物理组在协变密度泛函理论方面有深厚积累，他们利用经典r过程模型研究了核物理输入量对快中子俘获过程的影响，并建立了参数化的天体环境模型。广西大学团队用不同状态方程条件探索了r过程核合成和千新星发射的关系。中科院原子能院在核物理基础研究方面也在发力。
HIAF的未来角色值得关注。广东惠州的强流重离子加速器装置在2025年10月完成了全线束流贯通，2026年4月完成了首终端实验。HIAF可以产生极丰中子的核素，测量它们的质量和衰变性质——这些正是r过程计算最需要的输入量。虽然HIAF的主要科学目标不是核天体物理，但它的束流强度和能量精度将使中国首次有能力系统测量r过程路径上的核数据。
---
回到GW170817本身。八年过去了，它是我们唯一一次同时看到引力波和电磁对应体的中子星并合。
LIGO在2023年5月开始了第四次观测运行（O4），灵敏度比O3提升了约30%，双中子星探测距离约160-190 Mpc。但O4到目前为止还没有确认新的带电磁对应体的中子星并合。2026年4月，ApJ上发表了在O4运行中无引力波触发时搜索中子星并合千新星的结果——即使没有引力波信号，光学巡天也在扫描可能被遗漏的事件。2025年6月，arXiv上报告了O4运行中少数几个可能包含中子星的候选并合事件的千新星光学-近红外搜索结果——目前只有上限约束，没有确认的探测。
GW170817之所以特殊，部分原因是运气——它距离足够近（40 Mpc），使得电磁对应体可以被探测到。更远的并合，千新星太暗，看不了。LIGO估计银河系体积内双中子星的并合率约每10万年1-10次，这意味着在O4的探测范围内，每年可能有几个事件，但大部分太远，无法做光谱观测。
下一代引力波探测器将改变这个局面。欧洲的Einstein Telescope（ET）和美国的Cosmic Explorer（CE）——所谓的第三代探测器——预计2030年代中期上线。2025年5月，arXiv上综述了ET和CE的最新进展：ET的探测距离将达到O4的约10倍，CE更远。它们每年可以探测到数十万个双中子星并合事件，其中足够近、可以做千新星光谱观测的将有几十到上百个。
2024年8月，PRD上发表了用ET+CE模拟观测约束中子星状态方程的工作，发现第三代探测器可以把中子星半径的测量精度推进到100米以内——这比目前的千米级精度好一个量级。2026年1月，A&A上发表了第三代探测器对双中子星并合多信使观测的预测，估计ET+CE每年可以提供数百个有千新星光谱证认的r过程样本——到那时，r过程的图景将不再是几个数据点的推断，而是一种统计性的科学。
2025年9月，ET合作组正在确定最终选址，预算约20亿欧元。CE的设计也在推进中。物理世界（Physics World）在2025年9月报道了ET选址的最新进展——三个候选地点（荷兰-比利时-德国边境的Euregio Meuse-Rhine、意大利Sardinia、西班牙Canary Islands）正在进行最终的地质勘测。
---
回到r过程本身。我们走到哪了？
GW170817证实了中子星并合可以产生r过程元素——从锶到镧系元素，覆盖了第一峰到至少第二峰。千新星的光变曲线和光谱约束了喷出物质量（约0.04-0.06 M☉）、镧系质量分数（约0.5-2%）和膨胀速度（约0.1-0.3c）。
但中子星并合不是唯一的r过程场所。极贫金属星的观测要求至少两个场所：一个与铁伴生（磁转动超新星的候选），一个不与铁伴生（CEJSNe的候选）。中子星并合可能在中期和后期的宇宙中变成越来越重要的贡献者，但早期宇宙的重元素另有来路。
中微子驱动风作为"标准"核心坍缩超新星的r过程场所已经基本出局——Y_e太高，推不到重元素。但磁转动超新星保留了希望——它的喷流通道可以产生足够中子富集的条件。
核物理输入量仍然是最大的系统不确定性。5000个未知核素的质量、β衰变率、裂变概率——每一个都可能让r过程计算的结果偏移一两个峰。FRIB和HIAF正在逐步填补这些空白，但全面覆盖r过程路径可能需要十年以上。
数值模拟也在飞速进步。2025年的PRL论文标志着从"模拟碰撞后0.1秒"到"模拟碰撞后1.5秒"的跨越——这1.5秒里，盘风喷出、喷流形成、r过程冻结，所有关键物理都发生了。GRMHD+中微子辐射的模拟已经从2D推进到3D，从等温近似推进到全能量群转移。
第三代引力波探测器将在2030年代打开一个新的窗口：从"一个r过程样本"到"数百个r过程样本"。到那时候，r过程图景将从"推断"变成"统计"。
1957年，Burbidge等人写下r过程的概念时，他们不知道它在哪发生。2017年，我们看到了它在中子星并合中发生。2026年，我们正在意识到它不止在一个地方发生。宇宙的重元素不是一把锤子敲出来的——它是一整间工厂的产品，有不同的流水线，不同的班次，不同的质检标准。
我们刚刚走进这间工厂的大门。
---
关键时间线
| 时间 | 事件 | 来源 |
|---|---|---|
| 1957 | Burbidge等人提出r过程概念（B²FH论文） | Rev. Mod. Phys. |
| 1974 | Lattimer和Schramm提出中子星并合r过程 | ApJ |
| 1992 | Woosley和Hoffman计算中微子驱动风r过程 | ApJ |
| 2001 | Thompson等模拟中微子驱动风，Y_e偏高问题出现 | ApJ |
| 2014 | Martínez-Pinedo等总结：标准NDW仅产轻r过程 | J. Phys. G |
| 2017.08.17 | GW170817双中子星并合被LIGO/Virgo探测 | PRL |
| 2017.08.17 | GRB 170817A短伽马射线暴1.7秒后到达 | Fermi-GBM |
| 2017.10 | AT2017gfo千新星光谱发射；锶线证认 | Nature/Science |
| 2018 | 结构化喷流模型解释GRB 170817A | Nature/Mooley等 |
| 2019 | Watson等确认1μm处锶[Sr II]线 | Nature |
| 2021 | 磁转动超新星3D GRMHD r过程核合成 | ApJ |
| 2022 | GSI 𫓧化学实验确认挥发性金属 | Nature |
| 2023.01 | 系统双中子星并合喷出物核合成 | ApJ/Just等 |
| 2023.09 | SPARK I：首次千新星光谱丰度反演 | ApJ |
| 2024.02 | MNRAS：BNS并合可充分解释银河系重r过程 | MNRAS/Choplin等 |
| 2024.02 | SPARK II：多组分喷出物镧系丰度反演 | ApJ |
| 2024.08 | 金仕纶等提出CEJSNe r过程机制 | ApJ |
| 2024.10 | AT2017gfo早期→晚期r过程光谱特征逐小时浮现 | A&A/Gillanders等 |
| 2024.11 | FRIB产⁷⁸Ni附近极丰中子核素 | PLB |
| 2024.12 | 黑洞-中子星并合锕系增强r过程 | PubMed |
| 2025.01 | SPARK III：光谱+光变曲线自洽拟合 | arXiv→ApJ 2026 |
| 2025.01 | Nature Rev. Chem.：超重元素相对论化学综述 | Nat. Rev. Chem. |
| 2025.03 | 机器学习预测β衰变半衰期用于r过程 | PRC |
| 2025.05 | 最长BNS并合GRMHD模拟（~1.5s），喷流形成 | PRL |
| 2025.07 | "Europium, we have a problem"：化学演化困境 | A&A |
| 2025.09 | GW170817完整短综述 | Galaxies/MDPI |
| 2025.10 | GR辐射MHD模拟：自旋对r过程影响 | arXiv/Nedora等 |
| 2025.11 | 大规模β衰变率计算对r过程影响 | arXiv |
| 2025.11 | 多目标银河系化学演化约束r过程 | arXiv |
| 2025.12 | AT2017gfo镧系约束改进 | INSPIRE-HEP |
| 2026.01 | 第三代探测器多信使BNS观测预测 | A&A |
| 2026.03 | 校准镧系原子数据用于千新星辐射转移 | PRC |
| 2026.03 | Soker：七种r过程场景系统对比 | arXiv 2601.15187v3 |
| 2026.04 | AT2017gfo锶和氦NLTE建模最新确认 | arXiv/Aouad等 |
| 2026.04 | LIGO O4无GW触发时千新星搜索 | ApJ |