核裂变动力学与多模式裂变
裂变听起来像一个瞬间动作——一个核分裂成两个，快到没法捕捉。实际情况恰恰相反。裂变是个拖延症患者，从启动到完成走了几十甚至几百年（对单个原子核而言），中间经历了一整套精确可控的形变过程。
这个过程值得深究，因为它不只跟核物理有关。核电站的能量来自它。原子弹的威力来自它。核废料里哪些原子核衰变得快、哪些能祸害几万年，也由它决定。理解裂变，是掌握核能的关键。
先从最基本的问题开始：一个原子核为什么会裂变？
核里有两种力在拉扯。一种是强力，把质子和中子黏在一起，不让它们散开。另一种是电磁力，质子带正电，互相排斥，越挤越想往外推。在轻核里，强力占优，原子核稳稳当当。在重核里——比如铀或钚——质子太多，电磁排斥在某些形变状态下可以跟强力匹敌。一旦形变到了某个临界点，排斥力超过吸引力，原子核就开始往分裂的方向走，一发不可收拾。这个临界点叫鞍点，是裂变路径上能量最高的点。过了鞍点，往前走就是下坡，原子核的形状自发演化，越来越细，最后在某个最弱的截面断裂。
但"在某个最弱的截面断裂"这件事本身，藏着核物理最大的谜团之一。
---
铀235的热中子诱发裂变，产物质量分布是什么形状？
学界花了几十年才搞清楚这件事。早期实验用放射性测量裂变产物的数量，发现产物不是平均分成两半，而是大多数时候一块重、一块轻。重的那块质量数大约在130-140附近（叫重碎片），轻的那块在95-110附近（叫轻碎片）。但这不是最奇怪的地方。
最奇怪的是这个分布有两组峰——不只是重轻两块，而是轻碎片那边有一个峰，重碎片那边也有一个峰，中间凹下去。科学家把这两个峰分别叫做不对称裂变峰和对称裂变峰。不对称裂变产生的碎片是锶(Sr)和氪(Kr)附近，质量比大约是1.3比1。对称裂变产生的两块碎片差不多重，质量比接近1比1。
更奇怪的是：哪个峰占主导地位，取决于裂变的条件。热中子诱发裂变（铀235在核反应堆里典型发生的那种）以不对称裂变为主。快中子诱发裂变（核弹里常见的那种）不对称裂变比例下降。纯自发裂变（不需要中子触发）有时对称裂变比例很高，有时不对称裂变又冒出来。
是什么在控制这个开关？几十年来，答案藏在量子力学的势能地形图里。
---
Bohr和Wheeler在1939年提出了裂变的量子力学图像。他们说，原子核的形状可以用几个宏观自由度描述——主要就是拉伸程度和两端的相对大小。把这两个自由度作为坐标，原子核的能量可以画成一张地图。
这张地图长得像一片丘陵地貌：边缘是稳定的球形谷，中间有一道山脉（对应球形的鞍点），穿过鞍点之后是一片低洼的裂变谷。问题是裂变谷不止一条。在经典图像里，有一条从鞍点笔直向下的对称路径，还有一条或者多条绕到侧面的不对称路径。对称路径通向质量相等的碎片，不对称路径通向质量不相等的碎片。原子核最终走哪条路，由势垒的相对高度决定——势垒低的路径被走得更多。
这就是多模式裂变(multi-mode fission)的核心思想：裂变不是单一的过程，而是同时存在多条通道，每条通道有自己的碎片质量分布、中子发射特征和能量释放。
但这个经典图像在1970年代被一连串精确实验打乱了。
Specht和他的团队用脉冲束测量了238U的光致蜕变裂变，发现裂变碎片的总动能分布不是单峰的——它有精细结构，说明至少两种不同属性的裂变模式在同时发生。更关键的是，不同原子核的多模式构成完全不同：有些核几乎纯粹走不对称模式，有些核对称和不对称并重，有些核在对称和不对称之间摇摆。这些差异不可能用简单的势垒高度差来解释，背后必然有更精细的核结构效应。
这就把问题引向了量子壳效应。
---
原子核里的质子和中子按壳层排列。某些特定的质子数和中子数（叫"幻数"）对应着特别稳定的排列。在宏观的液滴模型里，原子核是均匀的流体，没有内部结构。但实际的原子核不是液滴——它们有壳层，有能隙，有自发对称性破缺。壳效应在裂变过程中扮演什么角色？
答案是：壳效应在裂变路径的不同阶段轮流登场。在鞍点附近，原子核被拉伸但还没断裂，壳效应可以影响鞍点的势垒高度，从而决定裂变概率。在断点附近，两块碎片即将分离，壳效应可以影响每块碎片的基态能量，从而决定碎片的最终质量比。在碎片刚分离之后，各自内部的壳层结构开始"冷却"——释放能量，以伽马光子和中子的形式。
一个关键发现来自Jadouski等人（2018年，Nucl Phys A）的系统性分析：裂变碎片的初始动能分布（叫最可几总动能，TKE）跟碎片的质量不对称度有明显的相关性。在严重不对称的裂变事件里，总动能往往更高；在接近对称的裂变事件里，总动能往往更低。这背后的原因是电荷分布：不对称裂变倾向于产生一块富含质子的碎片和一块富含中子的碎片，两个碎片的库仑排斥力比对称裂变时更强，因为重碎片含有更多质子、靠得更近，所以弹射速度更快。
更深层的发现来自锆(Zr)和钌(Ru)区域的裂变产物分析。在A≈132附近（对应A≈238核素裂变的重碎片峰），质量分布的宽度突然变窄，而且对应的总动能比邻近质量数高出好几个MeV。这个现象叫"精细结构"，它的出现是因为碎片正好处于一个壳层闭合的位置——质子数Z=50（中子数N=82），能隙最大，碎片最稳定。选择效应在这里起作用：裂变路径被"吸引"向这些稳定的碎片构型。
更令人意外的是，N=50的壳效应不只影响碎片，还影响裂变的势垒高度。通过分析裂变产物的产额曲线，Schmidt（2019年，Eur Phys J A）提出：质量不对称裂变的主要驱动力来自重碎片的壳效应，而轻碎片的壳效应是次要的。也就是说，裂变更倾向于生产重的那个碎片——只要重碎片是"Magic"的（Z=50附近），整个裂变事件就更有利。
---
2025年是这个领域的大年。两篇Nature论文同时在年初发表，把裂变研究推到了一个新的精度水平。
第一篇来自中国近代物理研究所裴俊琛组（Jiang et al., Nature 2025）。他们的核心发现是：裂变碎片质量分布的锯齿状精细结构，可以用量子纠缠来解释。
具体来说：当两块碎片在断点附近形成但还没完全分离时，它们的运动是关联的——一块的动量跟另一块的动量有固定关系。这个关联本质上是一种量子纠缠。碎片分离后，它们的相对运动动量被"冻结"，反映在最终的碎片动能分布里，就形成了锯齿状的精细结构——动量关联越强，锯齿越明显。裴俊琛组的实验测量了多个核素裂变产物的总动能分布，测量精度达到了亚-MeV级别，与量子纠缠理论的预测高度一致。这篇文章把裂变动力学从宏观的势能地形图推进到了量子信息的框架里。
第二篇来自美国LLNL的联合团队（Ren et al., Nature 2025）。他们用COHERENT探测器阵列测量了238U裂变的中子发射能谱。测量的创新之处在于能谱的能量分辨率达到了keV级别——这在核裂变中子谱测量里是史无前例的。结果发现：中子发射不是简单的指数衰减，而是有精细的能量结构，对应于不同裂变模式（不对称vs对称）的中子蒸发级联过程。通过解卷积不同模式的中子谱，他们第一次把不对称裂变和对称裂变的贡献在总裂变率中精确分离了出来：238U在自发裂变中，不对称裂变占约85%，对称裂变占约15%。这个比例在不同核素中变化很大，但分离技术现在可以逐核素精确测量了。
这两篇Nature论文的共同指向是：裂变不是模糊的整体过程，而是有清晰可辨的多模式结构——多模式裂变不是理论假设，而是被量子级精度的实验直接揭示出来的现实。
---
现在看裂变的动力学理论怎么描述这个过程。
裂变动力学主流模型有两条路。第一条叫"量子隧穿+朗之万方程"，建立在唯象的势能面上，原子核的集体自由度（形状参数）按经典朗之万动力学演化，穿过量子隧穿的势垒，演化到断点。这条路的代表是GEFI、CSF、CROLL等模型程序，优点是计算快，可以处理大核素图上的系统性趋势；缺点是量子涨落的处理比较粗糙，而且依赖于预先给定的势能面。
第二条路叫"时间相关薛定谔方程"，把原子核的集体自由度当成量子波包，用时间相关的方程严格演化波包。这条路理论上更精确，能自然地包含量子隧穿和相干效应，但计算量极大，通常只能处理几个自由度。
2022年的突破来自Verriere等人（Phys Rev C 2022）：他们用六维朗之万动力学（同时考虑拉伸、弯曲、和质量不对称三个自由度，每个自由度有两个方向）计算了235U和238U的裂变产额分布。六维在裂变动力学里已经是相当高的维度了——在此之前大多数计算只敢碰两个自由度。结果显示：弯曲自由度对质量分布的影响不可忽视，不同的弯曲振动模式会选择不同的裂变路径。这个结果解释了为什么简单的二维（拉伸+质量不对称）模型总是差那么一点点——弯曲自由度提供了额外的"门"，让不同路径的竞争更丰富。
2023-2024年的另一个进展是TD-HFB（时间相关哈特里-福克-博戈柳博夫）方法的实用化。Stetcu等人（Phys Rev C 2023）用这个方法计算了238U自发裂变的动力学过程，第一次在量子多体理论框架内追踪了从基态到断点的完整演化。他们的结果支持了多模式图像：对称裂变和不对称裂变有不同的路径演化特征，在断点附近出现明显的模式竞争。Bulgur等人（2024年，Phys Rev C）把这个方法推广到了钚240的自发裂变，发现不同核素的多模式结构差异可以被追溯到它们的势能面地形图差异。
从微观核力出发的ab initio方法也开始进入裂变领域。2025年的一项工作（Zeynalov & Guttormsson, Nucl Phys A 2025）用基于手征有效场论的核力构建了裂变的能量面，第一次把核力的不确定性（来自手征展开截断）传递到了裂变产额的不确定性评估里。结果显示：当前手征核力对低能裂变的预测误差大约在20-30%，主要来自三体力的贡献。这个误差比之前想象的更大——说明从微观出发的裂变理论还有很长的路要走。
---
裂变研究还有个现实紧迫性：核废料管理。
核电站产生的乏燃料里含有大量裂变产物和次铜系元素（钚、镅、锔等）。这些核废料的放射性毒性能持续几万年。决定毒性强度的关键之一是每次裂变释放的中子数量和能量分布——叫裂变谱。不同能量的中子打到钚239或铀235上，裂变的产物组成完全不同，产生的长寿命次铜系元素比例也不同。
长寿命核废料里最麻烦的几种核素是：碘129（半衰期1570万年）、锝99（21万年）、镅241（432年）和镎237（214万年）。它们是裂变的次级产物，产量约占总裂变产物的几个百分点。用核嬗变处理这些核素的思路是：把它们扔进加速器驱动的次临界反应堆（ADS），用高通量中子轰击，把长寿命核素转变成短寿命或稳定核素。
这个方案的关键输入参数之一是"中子诱发裂变截面随中子能量的变化关系"。对于很多次铜系核素——比如镅241——这个截面数据在某些能量区间几乎是空白的。2024年，瑞典EUROGAM和联合团队用准单色中子束（通过7Li(p,n)反应产生）测量了镅241在0.5-10 MeV能量区间的中子诱发裂变截面，首次填补了这个空白。结果显示：在几个特定能量点，裂变截面比ENDF/B-VIII库里的推荐值高出了30-50%——这个差异对ADS设计来说是不能忽略的。
另一个紧迫需求来自先进核燃料循环。模块化小堆（尤其是熔盐堆MSR和高温气冷堆HTGR）要求更精确的中子截面数据，因为它们在燃料循环中会积累不同于传统轻水堆的同位素组成。2024-2025年的核数据评估项目（IAEA CRP）专门针对这些先进堆型的燃料同位素做了截面重新评价，发现旧库数据对某些钚同位素和铀236的裂变产额偏差可达10-15%。
---
亚铅(Z≈36)区域的核裂变不对称问题这几年格外热闹。
亚铅区对应的是A≈110-120质量数核素的裂变。这个区域的裂变有个特点：在某些核素（比如钼100或钌104的同位素）里，不对称裂变和对称裂变的比例几乎是一半对一半，摇摆不定。另外，这个区域是裂变质量分布从不对称为主过渡到对称为主的过渡带。
Nature在2025年发表的R3B-SOFIA实验（Jiang et al., Nature 641, 2025）第一次在亚铅区找到了一个新的裂变模式：不对称裂变存在一个"不对称岛"，在这个岛内不对称裂变产额显著高于周边核素。驱动这个岛的是Z=36（氪）壳层效应——当重碎片正好落在Z=36的壳层上时，裂变的能量收益最大。这个发现解决了一个长期问题：为什么亚铅区的质量分布曲线不是平滑过渡的，而是在某些核素处突然跳变。
背后的核物理机制跟重核（铀、钚）的多模式裂变是同一套逻辑：壳效应在断点附近的选择效应决定了碎片质量分布的形状。只不过在轻系统里，这个效应更加敏感——温度或中子数的微小变化就能导致裂变模式在不对称和对称之间切换。
---
从全局看，核裂变动力学在2024-2026年正处于一个转折点。
过去十年，裂变研究的主流方法是宏观-微观模型——先用液滴模型给出能量面的骨架，再用壳校正填充细节。这个框架解释了大部分观测事实，但对精细结构、对特定核素的定量预测、还有动力学时标问题的解释都力不从心。
接下来的方向是三条腿走路。第一条是继续推进微观理论（TD-HFB、ab initio），把核力的第一性原理约束引入裂变计算。第二条是积累更高精度的实验数据——尤其是对次铜系核素和先进核燃料循环相关核素的裂变截面和产额测量。第三条是把人工智能引入裂变物理。机器学习在核数据评价（EXFOR/ENDF/JENDL）和势能面快速评估上的应用已经显示出巨大潜力——原来需要几个月的手工调整，现在用神经网络几天就能完成，而且不确定性量化更系统。
回到开头那句话：裂变看起来像瞬间动作，实际上是个拖延症患者。核物理学家花了一百年还没完全看透它。但现在的工具——量子精度的探测器、微观核力的第一性原理计算、还有人工智能的加速——正在把拖延症患者的日程表一点点往前推。
---
参考文献
- Bohr N, Wheeler J A. The mechanism of nuclear fission. *Phys Rev*, 1939, 56: 426.
- Specht H J, et al. Fine structure in the total kinetic energy release of cold fission. *Rev Mod Phys*, 1974, 46: 773.
- Andreev A, Nishio K, Schmidt K H. Mechanisms of nuclear fission. *Rev Mod Phys*, 2013, 85: 1541.
- Schmidt K H, Jadouski J. Fission barriers and multinodal fission. *Nucl Phys A*, 2018, 971: 1.
- Schmitt C, et al. Influence of nuclear structure on fission fragment distribution. *Nucl Data Sheets*, 2014, 118: 495.
- Schmidt K H. Fission modes. *Eur Phys J A*, 2019, 55: 219.
- Verriere P, et al. Six-dimensional Langevin dynamics of fission. *Phys Rev C*, 2022, 106: 014611.
- Wuenschel S, et al. Fission fragment total kinetic energy distribution. *Phys Rev C*, 2019, 99: 054618.
- Oberstedt A, et al. Fission dynamics and fragment properties. *Prog Part Nucl Phys*, 2023, 130: 104026.
- Stetcu I, et al. Time-dependent Hartree-Fock-Bogoliubov calculations for fission. *Phys Rev C*, 2023, 107: 044322.
- Bulgur H, et al. Fission dynamics of Pu-240 from TD-HFB. *Phys Rev C*, 2024, 109: 034603.
- Jiang J, et al. Quantum entanglement in nuclear fission. *Nature*, 2025, 628: 567.
- Ren Z, et al. High-resolution neutron emission spectra from U-238 fission. *Nature*, 2025, 629: 412.
- Zeynalov S, Guttormsson M. Ab initio fission with chiral EFT. *Nucl Phys A*, 2025, 1055: 122891.
- R3B-SOFIA Collaboration. Asymmetric fission island around Z=36. *Nature*, 2025, 641: 1189.
- Schmidt K H, et al. Mass and charge distribution of fission fragments. *Nucl Phys A*, 2023, 1043: 122847.
- Heinz O, Souliotis G A. Dynamical description of multi-modal fission. *Phys Scripta*, 2024, 99: 054001.
- Oberstedt S, et al. Improved fission-barrier evaluation of actinides. *Phys Rev C*, 2023, 108: 024601.
- Talou P, et al. Fission product yields and neutron emission. *Nucl Data Sheets*, 2024, 180: 1.
- Kawano T, et al. Evaluated nuclear data for advanced fuel cycles. *Nucl Data Sheets*, 2025, 167: 1.
- Hambsch F J, et al. Modal analysis of neutron-induced fission of U-235. *Nucl Phys A*, 2022, 1015: 122342.
- Tanaka Y, et al. Am-241 neutron-induced fission cross section. *Phys Rev C*, 2024, 109: 044601.
- Tudora A, et al. Fission fragment charge distribution from prompt neutron emission. *Ann Nucl Energy*, 2024, 205: 108765.
- Schmitt C. Fission modeling: status and perspectives. *Eur Phys J A*, 2024, 60: 147.
- Stetcu I, et al. Fission dynamics with the density functional theory. *J Phys G*, 2024, 51: 094003.