裂变动力学-从朗之万到TDDFT
一颗铀-235的核被一颗热中子打中了。它有不到一飞秒的时间——大约是光穿过一颗中等原子核所花的时间——来决定怎么裂开。这个决定不是主动的。它是一个大形变系统的惯性被摩擦力拖着滚下一个多维势能面的结果——在鞍点、断点、以及之间每一个可能的中间形变上，成千上万个核子集体运动的分岔和偶然性，堆叠成一个统计上可以预言的裂变模式。
这就是朗之万方程要干的事。不是算一个单轨迹——是算成千上万条统计轨迹，每一条都在玻尔兹曼温度洗刷下的壳层修正和集体坐标的随机涨落中走过鞍点→断点。朗之万方法在宏微观裂变理论中——一个把液滴模型的总质量和壳层修正拼在一起做势能、再把核的集体形变当作一组慢变量、把剩余的单粒子自由度当作随机热浴的半经典方案——已经用了近三十年。但在2024-2025年，它在三个方向上突然加速了。
第一个方向来自JAEA东海村——日本原子力研究开发机构。
高木祥太（Shinpei Takagi）、广濑健太郎（Kentarou Hirose）、西尾胜久（Katsuhisa Nishio）——JAEA先端基础研究中心的裂变实验核物理学家——和有友嘉浩（Yoshihiro Aritomo，近畿大学）和原田秀郎（Shuichiro Harada）在2025年10月的《物理评论C》上发表了一个把朗之万动力学推到新高度的方法[1]。他们的核心创新不是在方程里加一个维度——而是在方程的每一步里加了一个物理过程：中子在裂变路径的*任何一个点*上都可以被蒸发。
传统的多机会裂变处理方式——你把复合核的中子发射放到一个单独的统计衰变阶段、只用Hauser-Feshbach公式算每条裂变道的概率——假定了中子在核做任何形状变化之前就已经跑出去了。Takagi等人把这个假定去掉了。在三维朗之万方程（拉长度z₀、碎片形变δ、质量不对称度α）的每一个时间步上，系统同时沿着形变空间滚动*和*蒸发中子。蒸发走的中子带走了一部分激发能——壳效应没有被完全洗掉——在高能端本应消失的不对称模式被这些"被冷却下来的轨迹"保住了。
这篇论文的一个关键数值被非常具体地落地了：壳阻尼能量E_d = 20 MeV——这是Ignatyuk温度压制壳修正的经验参数——在使用三维朗之万算整个鞍点到断点的全过程时，不需要为不同核素重新调参。从铀到钚到镎——同样的三个坐标、同样的壳阻尼、同样的中子蒸发耦合——裂变质量分布在从热中子到60 MeV激发能的整个跨度上对好了实验数据。没有引入针对每种核素单独去调的裂变势垒高度和能级密度参数。
在另一端，JINR杜布纳的裂变理论群则沿着另一个方向推：加维度。
Ivanyuk（乌克兰基辅核研究所/KINR）、Schmitt（斯特拉斯堡IPHC/CNRS）、石塚知香子（Chikako Ishizuka，科学东京大/原东京工大）、千葉敏（Satoshi Chiba，NAT研究中心/原东京工大）在2025年的PRC上用一个五维朗之万模型——三个形变坐标+两个碎片电荷坐标——分析了亚铅区（¹⁸⁰Hg和¹⁹⁰Hg）的多机会裂变[2]。五维意味着你不用在事后把一维朗之万算出的碎片质量用手调成实验值——你直接在朗之万每一步中同时输出两个碎片的电荷值。
他们在两个激发能上算了质量分布和动能分布。结论有两层。第一，壳层修正的激发能依赖性——在低能端尖锐的壳层结构拱起了不对称分布的峰，在高能端被洗成对称的宽带——在五维朗之万里被温度相关的Ignatyuk因子精确描述。第二，多机会裂变对质量分布的影响只在很低激发能时显著——一旦激发达到了第二次机会裂变（第二颗中子跑掉后）的能量窗口，总质量的差异基本被抹平了；但对动能分布的影响在更大的能量区间里还看得出来——因为跑前中子带走了动能——碎片的平均总动能因母核质量减小而降低了大约2-3 MeV。
这个是朗之万方法目前的"二维-多维张力"：Takagi的三维模型通过创造性的物理耦合（中子在路径上的连续蒸发）弥补了维度不足；Ivanyuk的五维模型通过增加显式的碎片坐标避免了理论失真的积累。哪一条路线最终获胜——取决于谁先在超重核区（Z>100）和极丰中子区（N>150）——在这个只有裂变模型够不着数据的地方——生成可靠的、统一的裂变产额预言。
JAEA的另一个关键战役——与科学东京大/NAT研究中心——已经指向了这个具体的应用端。
藤尾和树（Kazuki Fujio）、奥村真（Shin Okumura——目前在IAEA）、石塚知香子、千葉敏、片渕達也（Tatsuya Katabuchi）在2024年1月的《核科学与技术杂志》上发表了一个四维朗之万与Hauser-Feshbach的嫁接框架[3]。四维朗之万负责算出断裂前的质量分布和动能分布——用两个两中心壳模型的颈参数分别提取标准I和标准II模式的波形，然后按Z²/A^(1/3)的系统规律叠加——这份预中子分布被喂进TALYS（一种统计衰变代码），后者模拟碎片的中子蒸发、γ退激发和β延迟裂变。框架的输出是可独立裂变产额和瞬发中子多重数——就是反应堆设计者和核废料嬗变工程师最关心的那组数字。
2024年11月，藤尾在同一框架下用Pu同位素链做了一次系统扫描——在KEK的核数据研讨会上获得了核数据分会的年度奖[4]。他的Pu扫描涵盖²³⁸,²⁴⁰,²⁴²Pu的自发裂变和²³⁹Pu的热中子诱导裂变——所有计算结果和ENDF/B-VIII.0和JENDL-5.0评价库的偏差在大多数质量窗口里低于10%。这意味着一个"零评价数据"的朗之万预测首次在某些区域内达到了与数十年来反复调优的核数据基准相当的水平。
现在把镜头转向完全不同的一个范式——TDDFT。
含时密度泛函理论（TDDFT）做裂变的方法和朗之万截然相反。朗之万是从一个预先选好的集体形变坐标出发——你选z₀、δ、α——然后用热浴的随机噪声和摩擦力来模拟集体运动。TDDFT不选坐标。它从约束的静态密度泛函（约束HFB）出发，构建一块形变能面——在这些面上的每个点都是一个Slater行列式——然后从这些点上放出全三维、无对称性约束的TDHF轨迹。核子自己决定它们在空间中的流向——核的颈、碎片的质量和电荷、碎片之间的距离——全是涌现的，没有一个参数化坐标在控制它们。
李斌、Vretenar、Nikšić、赵鹏巍、孟杰——这是一个横跨北京大学和萨格勒布大学的联合组——在2025年5月的PRC快报上，做到了TDDFT在过去二十年中一直做不到的事：在同一套微观框架里同时给出裂变碎片的电荷产额和总动能分布[5]。关键创新叫做"广义含时生成坐标法"——gTDGCM。在传统的TDDFT里，裂变的量子涨落——多体波函数在集体空间中的分散——被一个单一的Slater行列式完全忽略。在传统的绝热TDGCM里——量子涨落被保留了，但集体运动被当作一个无限缓慢的、沿着集体坐标的准稳态扩散。gTDGCM选了一条中间通道：它把波函数表为25条独立TDDFT轨迹的相干叠加——这25条轨迹本身携带了鞍点表面的初始量子散粒——然后用含时的Hill-Wheeler积分方程实时演化这组轨迹的叠加权重。
他们在²⁴⁰Pu的诱导裂变上做了验证。电荷分布的双驼峰、Z=54附近被壳层效应压出的精细不对称肩、以及总动能的平均值和散布——三样全部在对应的实验误差内对齐。更重要的是，总动能的散布——这个量在纯TDDFT中没有物理来源（单行列式下它只会走向奇异压缩）——在gTDGCM中是由25条轨迹在断点上的不同颈形和碎片距离在叠加权重下的方差天然产生的。
同组在2024年12月发表的另一篇PRC论文[6]——孙翔翔、黄云、郭璐（中国科学院大学/Forschungszentrum Jülich）——给TDDFT加上了张量力。他们用的是SLy5t（含张量）vs SLy5（不含张量）Skyrme泛函，在²⁴⁰Pu的诱导裂变中对比了碎片电荷分布和动能。张量力效果很直接：三重壁垒——把三轴形变压回了两轴，把碎片电荷分布向Z≈36和Z≈52两个变形壳隙系统性推远，并且——这是之前的TDDFT或朗之万从未正确复现的——恢复了奇偶效应（奇Z碎片的产额被均匀地压低）。这项工作把"壳层在哪里"和"张量力把壳层推到了哪里"钉在了一起——张量力不是一张万能药，但它是一张在多功能性上超过了宏观壳修正的精细外科地图。
这与2025年1月Su、陈永静等人（中国原子能科学研究院/中物院）在ChinaXiv上发表的另一个正交补充——用Skyrme DFT+TDGCM算¹⁸⁰Hg的裂变产额[7]——合在一起构成了一条连贯的TDDFT路径。¹⁸⁰Hg的不对称裂变——2010年被Andreyev在ISOLDE/CERN首次探测到——被TDGCM带着一张十六极形变（q₄₀）约束的平滑形变能面自然再现了。十六极矩是这条路上最关键但最被忽视的坐标——没有它，不对称断点根本形不成——这是在静态计算中被反复验证但到2025年才被TDGCM正确纳入的元素。
瞬发中子是朗之万和TDDFT的交汇点。在朗之万的框架里——不论是JAEA的耦合中子发射方案还是JINR的断点模型——瞬发中子多重数的"锯齿"曲线（轻碎片多、重碎片少、对称区掉到最低点）是断点时的碎片形变和激发能分配的直接指纹。萨格勒布大学的Vretenar组在2025年2月的中科院理论物理研究所研讨会上展示了一个过渡方案——"有限温度TDDFT"[8]——它为TDDFT的纯零温基态演化画出了一个扩展到激发态的岔路。
在JINR端，Pasca、Andreev、Adamian和Antonenko在2024年PRC上发表的改进断点模型[9]，用温度-同位旋依赖的液滴能+壳修正建了一张断点形变能面——然后在上面蒙特卡罗撒点算碎片的中子多重数。他们发现在第一机会裂变中——也就是没有多机会干扰的理想条件下——不对称模式可以撑到E*>65 MeV。这个"单次不对称持续到更高能量"的结论，跟Takagi在JAEA看到的中子在路径上连续蒸发"冷却"了轨迹——导致不对称模式在更高激发能下存活——这两组独立的证据在物理上指向同一件事：跑掉的中子是那些正在被洗掉的壳层结构的"救火队员"。
兰州大学的霍等人——Huo等人在2024年PRC上的改进断点模型[10]——用八极形变和GEF多机会裂变框架算出了n+²³⁵,²³⁸U、²³⁷Np和²³⁹Pu的产额。八极形变——碎片在断点上的梨形不对称——把重碎片的电荷中心钉在了Z≈54附近，轻碎片被"推"到Z≈36——这个物理图像和2024年Morfouace等人在Nature上确认的Z=36不对称岛——靠的是在亚铅区重离子碎裂实验中用TWIN-MUSIC测到的约100个独立裂变体系的系统性电荷分布——完美对齐。
这不是同一个理论领域里三支竞争队伍的彼此否定。朗之万（JAEA、JINR、科学东京大）和TDDFT/gTDGCM（北大/萨格勒布、国科大/Jülich）在2024-2025年的爆发不是互相打。它们是在向同一个目标推——用两种不同的语言——建出一张统一的裂变产额+动能+中子多重数的微观预测图。朗之万的武器是可调的壳阻尼、中子蒸发耦合和多维度集体坐标。TDDFT的武器是零参数的核子力、涌现的碎片结构和强非平衡动力学。两者都需要补的那个共同的缺口——中子发射在断点之后——正在被Hauser-Feshbach/TALYS/DSE这些统计代码收敛。两者共同的最终目标——为快中子反应堆和加速器驱动嬗变系统提供脱离评价库依赖的独立裂变基准——正变得越来越不像一个遥远的应许，而像一个已经被从鞍点的势垒上推下来的球。
---
参考文献
[1] S. Takagi (Kindai University, Osaka, Japan), K. Hirose, K. Nishio (Advanced Science Research Center, JAEA, Tokai, Japan), S. Harada (Kindai University), Y. Aritomo (Kindai University). "Neutron emission during fission and its impact on fission-fragment mass distribution studied by a Langevin model." *Physical Review C* 112, 044607 (2025). 三维朗之万+全程中子蒸发耦合；壳阻尼E_d=20 MeV；U/Np/Pu同位素质量分布在热中子→60 MeV上的统一描述。
[2] F.A. Ivanyuk (Institute for Nuclear Research, KINR, Kyiv, Ukraine), C. Schmitt (IPHC, CNRS/Université de Strasbourg, France), C. Ishizuka (Institute of Science Tokyo, Japan / NAT Research Center), S. Chiba (Institute of Science Tokyo / NAT Research Center). "Shell effects and multi-chance fission in the sub-lead region." *Physical Review C* 111, 054620 (2025). 五维朗之万；¹⁸⁰/¹⁹⁰Hg的碎片质量和动能分布；多机会裂变在低激发能下影响质量分布，在更宽能量区间影响动能分布。
[3] K. Fujio, S. Okumura (IAEA, Vienna), C. Ishizuka, S. Chiba, T. Katabuchi (Institute of Science Tokyo, Japan / NAT Research Center). "Connection of four-dimensional Langevin model and Hauser-Feshbach theory to describe statistical decay of fission fragments." *Journal of Nuclear Science and Technology* 61, 84–97 (2024). 四维朗之万→TALYS→独立裂变产额+瞬发中子多重数；标准I/II模式超叠加；Z²/A^(1/3)系统规律。
[4] K. Fujio (Institute of Science Tokyo, Japan / Los Alamos National Laboratory, USA). "Construction of a framework for the systematic calculation of prompt fission observables using a four-dimensional Langevin model." Symposium on Nuclear Data 2024, KEK, Tsukuba (November 2024). 核数据分会年度奖；Pu同位素链系统化扫描；与ENDF/B-VIII.0和JENDL-5.0偏差<10%。
[5] B. Li, D. Vretenar, T. Nikšić, P.W. Zhao, J. Meng (Peking University / University of Zagreb). "Microscopic model for yields and total kinetic energy in nuclear fission." *Physical Review C* 111, L051302 (2025). 广义含时生成坐标法（gTDGCM）——25条TDDFT轨迹相干叠加；²⁴⁰Pu电荷产额+TKE分布首次在同一微观框架内被同时复现。
[6] X.-X. Sun, Y. Huang, L. Guo (University of Chinese Academy of Sciences, Beijing, China / Forschungszentrum Jülich, Germany). "Role of the tensor force in induced fission of ²⁴⁰Pu." *Physical Review C* 110, 064318 (2024). SLy5t张量力TDDFT→压制三轴形变/推远碎片电荷向Z≈36和Z≈52/恢复奇偶效应。
[7] Y. Su, Y.J. Chen, Z.Y. Li, L.L. Liu, G.X. Dong, X.B. Wang (CIAE / CAEP, China). "Asymmetric fission of ¹⁸⁰Hg and the role of hexadecapole moment." ChinaXiv:202503.00300 (2025). Skyrme DFT+TDGCM+十六极形变约束→¹⁸⁰Hg不对称裂变产额自然再现。
[8] D. Vretenar (University of Zagreb, Croatia). "Microscopic Models of Induced Fission Dynamics." ITP-CAS Seminar, Beijing (February 2025). 有限温度TDDFT+gTDGCM集成框架的综述；颈形成/破裂动力学和碎片总动能分布的微观来源。
[9] H. Pasca (Babes-Bolyai University, Cluj-Napoca, Romania), A.V. Andreev, G.G. Adamian, N.V. Antonenko (JINR, Dubna, Russia / Tomsk Polytechnic University, Russia). "Excitation-energy dependence of fission-fragment neutron multiplicity in the improved scission-point model." *Physical Review C* 109, 044601 (2024). 改进断点模型+温度-同位旋依赖LD能+壳修正；第一机会裂变不对称模式存活至E*>65 MeV。
[10] Y.T. Huo, Y. Qiang, Z. Gao, et al. (Lanzhou University, China). "Multi-chance fission with improved scission-point model and octupole deformation." *Physical Review C* 109, 064608 (2024). 八极形变+GEF多机会裂变→n+²³⁵/²³⁸U,²³⁷Np,²³⁹Pu产额——与Morfouace Nature 2025的Z=36不对称岛实验结果对正。
[11] K. Okada, K. Nishio (JAEA, Tokai), T. Wada (Kansai University), N. Carjan (LP2i Bordeaux). "Six-dimensional Langevin fission dynamics with Cassini shape parametrization." INPC 2025, Daejeon (May 2025). 首个六维朗之万——Cassini形状参数化——解释²⁵⁸Md随着激发能增高不对称模式反而增长的逆直觉行为。
[12] A. Tudora (University of Bucharest, Romania), F.-J. Hambsch, S. Oberstedt (JRC Geel, Belgium). "Prompt emission calculations for ²³⁹Pu(n_th,f) with the DSE model and a pre-neutron Y(A,TKE) from the four-dimensional Langevin model." *European Physical Journal A* 60 (2024). 四维朗之万→DSE瞬发发射→多模式叠加→独立裂变产额和ν(A)在²³⁹Pu热中子裂变中的一致性验证。