跳到正文
返回核结构
至真探索 · 核结构

五维朗之万-壳效应如何在高热中活下来

2026-06-24

一颗铀核吸收一个热中子后，它接下来的几皮秒里会发生两件互相矛盾的事。一是它要变形——从一颗球拉成花生、把脖子越掐越细、最后断开。二是它在变形过程中会变热——激发能从零开始逐渐转化为内部的热运动。热了之后，核子壳层的量子效应——那些决定了核是"喜欢"对称裂还是不对称裂的精细结构——应该被洗掉。就像一个冰块掉进热水里——你应该看到它融化。

但你没有。裂变碎片的分布告诉你，即使激发能达到三四十兆电子伏——对于原子核来说这已经是"滚烫"——那些壳层效应仍然依稀可见。轻碎片和重碎片的质量峰还是分得开的。某种量子记忆在极端热浴中顽强地存活了下来。

这背后有一个数学模型——五维朗之万方程——和一个物理学概念——壳修正阻尼。它们合在一起，构成了当代裂变理论中最重要也最微妙的一个张力：壳效应到底能活多热？

2024年3月，F.A. Ivanyuk、C. Ishizuka和S. Chiba在*Physical Review C* 109, 034602上发表了一篇论文，标题朴实无华：*Five-dimensional Langevin approach to fission of atomic nuclei*。Ivanyuk是乌克兰人，长期在基辅的核研究所工作；Ishizuka是Science Tokyo（原东京工业大学）的副教授；Chiba是NAT公司的理论物理学家。这篇论文建立了一个从宏观-微观势能面出发、用两中心壳模型参数化核形状、再在五个独立的集体坐标上演化朗之万方程的完整框架。

为什么要五维？三维不够吗？不够。三维朗之万——通常取拉长度（脖子到脖子距离的一半）、碎片形变和质量不对称度——能告诉你裂变的大致方向（对称还是不对称、碎片大致大小），但它漏掉了裂缝形成过程本身。第四维——碎片的独立形变（δ₁和δ₂不相等，即轻碎片和重碎片可以不同程度地变形）——可以改善碎片动能分布。但要把碎片的壳结构从演化过程中自然挤出来——不是作为输入假设而是作为输出结果——你必须加第五维：颈参数ε。它控制的是两个碎片之间连接桥的粗细——你让它自己演化，而不是强行关闭它，壳效应就能在断裂前从波函数里自然浮现出来。

这篇2024年论文的最关键技术贡献是一个叫"软墙"的东西。朗之万轨迹在集体坐标空间里随机游走——但在某些区域，势能面不再有意义（比如两个碎片已经分得太远）。之前的处理方法是硬墙——轨迹碰到边界就反射回去——这等于人为注入了能量。Ivanyuk等人换了一种边界条件：轨迹碰到"墙"时，摩擦力急剧增大但势能面平滑延展——这就是软墙。轨迹自己慢下来，而不是被弹回去。这个创新对高温裂变尤其重要——因为在高温下，壳修正变弱，势能面对拉长方向的约束变小，轨迹更倾向于走向断裂而不是在山谷里来回振荡。硬墙会制造假振荡，软墙不会。这篇论文的结论之一是：加了第五维后，钍同位素的三峰质量分布（质量对称+不对称+超长对称混合结构）首次被自然地重现了——不需要人工调参数。

然后2025年5月，同一组人——这次Ivanyuk合作了C. Schmitt（法国CNRS/IPHC斯特拉斯堡）——把五维朗之万模型推到了亚铅核素区（*Physical Review C* 111, 054620, 2025年5月，编辑推荐论文）。研究对象是¹⁸⁰Hg和¹⁹⁰Hg——缺中子的汞同位素。这些核裂变时碎片质量峰不是标准的两峰——它们是双峰中间夹一个谷，像驼峰。五维朗之万自然地产生出了这个双峰形状。

但真正让他们吃惊的事发生在一个被传统假定为"不相关"的区域。壳修正阻尼的标准公式——Ignatyuk公式——叫它Φ(T)=exp(−E*/E_d)，E_d≈20 MeV。它说：激发能每多20 MeV，壳修正的大小就缩小到原来的约三分之一。到了40MeV，壳效应应该已经衰减到不足最初的十分之一。也就是说，汞-180被加热到大约一颗瓦亮灯泡的热扰动时，它的碎片质量分布应该已经变成了一个宽而平的对称峰——所有量子记忆都被洗白。但Ivanyuk等人的计算表明：即使E*=40-50 MeV，¹⁸⁰Hg的双峰分布仍然清晰可见。壳效应没有在20 MeV衰减——它顽强地活到了40-50 MeV。

这和Ignatyuk公式矛盾——但差距是大是小，取决于你拿哪组实验来校准。Ignatyuk的E_d≈20 MeV原来是从锕系核（铀、钚）裂变位垒高度和蒸发谱实验里提取的。Ivanyuk等人的工作暗示：亚铅区核的壳效应阻尼常数可能比锕系大——这在物理上合理，因为较轻的核壳效应本身较弱，需要用更强的"热浴"才能洗掉。但问题是，没有人真正直接测量过壳修正随温度的函数——它是从位垒高度和裂变概率的反推中间接推断出来的。而Ivanyuk的五维计算是直接从第一原理（宏观-微观势能面加Strutinsky壳修正）算碎片质量分布——等于绕开了位垒高度反演链条，直接用产额数据来约束阻尼常数。

这两个结论不矛盾——它们互补。Ignatyuk公式拿位垒缩水来推断壳修正衰减；五维朗之万拿碎片质量峰的形状来推——一个从上往下看，一个从下往上看。中间差了二十年，差了一个维度。

现在说朗之万方程本身。这不是裂变物理学发明的——是物理化学。Paul Langevin在1908年为描述花粉颗粒在水中的布朗运动写下了这个方程：一个质量乘以加速度的项、一个与速度成正比的摩擦项、一个由势能梯度决定的驱动力项、再加一个与温度有关的随机力项。在裂变物理里，你把这些量翻译成核的集体坐标：拉长度c的加速度乘以惯性张量、一体耗散（墙-窗模型）产生的摩擦力、势能面斜率提供的驱动力、以及满足涨落-耗散定理的随机噪声。

惯性张量——裂变系统对不同变形方向的"有效质量"——通常用Werner-Wheeler不可压缩流体近似来算。摩擦——墙-窗模型——假设核子在核表面和颈部的碰撞造成集体运动能量的耗散。随机力——让它满足爱因斯坦关系：力的大小与摩擦系数和环境温度平方根成正比。这整套机器的物理一致性依赖于一个脆弱的事实：你得知道每一步的"瞬时温度"。温度的来源是集体运动转换成内禀激发的耗散——但反过来，内禀激发又会降低壳修正（壳修正随温度衰减）。这两个过程互相耦合，构成一个非线性反馈回路——壳修正的衰减改变了势能面，改变了势能面就改变了朗之万轨迹的概率分布，改变了轨迹就改变了摩擦生热的速率，改变了摩擦生热就改变了壳修正的衰减速率。联立的。不能先算一个再算另一个。

这就是为什么五维朗之万能揭示壳效应超长寿——而传统的单维或多维位垒穿透模型不能。在单维里，你假设温度场是空间均匀的——整个核在每一帧都有相同的温度——然后你把壳修正按这个均匀温标衰减掉。但五维朗之万告诉你：碎片和颈部附近的温度是不均匀的——颈部的核子碰撞最剧烈，那里的温度最先上升。壳修正的衰减不是全局关闭——它是从颈部开始、向碎片中心逐步扩散的。碎片的质量峰主要靠碎片内部——也就是温度最低、壳修正存活最久的那部分空间——的势能面结构来决定。这个空间不均匀性——不只是核形状不均匀，温度也不均匀——是壳效应能在高热环境中苟活的关键物理机制。

2024-2025年日本原子能研究开发机构（JAEA）的两项平行工作从不同角度印证了同一条物理规律。首先，Y. Aritomo、K. Nishio等人——日本在裂变朗之万动力学的长期领军人物——在2024年核数据研讨会上报告了六维朗之万的首次尝试。用Cassini形状参数化（六个集体坐标：{α, α₁, α₃, α₄, α₆}），他们模拟了²⁵⁸Md（钔-258）的裂变——发现当激发能从15 MeV升到18 MeV时，短对称裂变分量急剧下降，原来被壳效应压制的非对称道成为主导。这不是壳效应的"衰减"——这是壳效应在能量窗内的选择性放大：低能时对称模式壳修正强（驱动短对称道），高能时这部分壳修正首先被洗掉，整个系统转向液滴模型更偏好的非对称道。六维之所以必要，是因为短对称道和超长对称道（另一个对称裂变模式）在质量数是重叠的——不加复合形变的独立自由度（α₁, α₃, α₄, α₆独立而非一个单一形变参数），你分不清这两个道的激发能依赖。

第二项JAEA工作（发表在*Physical Review C* 112, 034613, 2025年9月）比较了五维朗之万中使用宏观惯性张量（Werner-Wheeler）和微观惯性张量（基于单粒子能级的微扰展开）的差别。宏观惯性给出的碎片形变范围偏窄，趋于低估非对称裂变产额；微观惯性给出的分布更宽，更好地拟合了实验。但微观惯性的计算代价高出一个量级——每一个时间步都要重新对角化单粒子哈密顿量来算惯性。这篇论文还报告了一个奇怪的物理现象：在重碎片质量A≈134区域，出现了一种"热裂变模式"——碎片是由一个超拉长的轻碎片和一个紧凑的重碎片构成，颈形状很像一个极度伸展的桥。（这个现象的原因还在探索——可能是β₃和β₄形变之间的耦合共振。）

2024年3月，Huang等人在*Physical Review C* 109, 034609上用三维朗之万模型研究了¹⁸⁷Ir（铱-187）的裂变多模性。他们发现了四个共存的裂变模式——两个不对称标准模式、一个对称超长模式、一个对称液滴模式——藏在同一个看似单峰的碎片质量分布后面。壳修正阻尼参数取了E_d=28 MeV（比Ignatyuk默认的20 MeV更高），因为铱的质量区（比铀轻了两个质子壳层）壳效应本身弱——需要更长的阻尼尾巴才能在实验激发能下让不对称模式存活。这组作者的结论和Ivanyuk等人一致：阻尼常数E_d不是一个普适常数——它依赖核素区、依赖裂变模式、甚至可能依赖激发能的绝对值。

现在来谈谈一个微妙的反直觉点。

壳修正这个词本身可能引起误解。Strutinsky在1967年发明的壳修正方法——把液滴模型的光滑总结合能减去单粒子能级的壳振荡——给你的是一个"温差"。不是温度本身——而是量子壳层在零温度下贡献的结合能增量。当温度上升时，壳修正的衰减是不是严格的指数规律——Ignatyuk的Φ(T)=exp(−E*/E_d)只是个参数化近似。Krappe在1979年提出的原始壳修正阻尼形式更复杂——Φ(T)=τ/sinh(τ)当τ=2π²T/ħω——但多数朗之万计算不用这个，因为多一个拟合参数。

这引出了一个更深的问题：如果你用不同的壳修正阻尼形式，五维朗之万的结论——壳效应存活到40-50 MeV——变不变？目前没有人系统性地回答过这个问题。不同研究组各自用各自偏好的阻尼常数——JAEA用E_d=20 MeV，Huang等人用28 MeV，Ivanyuk在2025年的Hg论文里则用了"壳修正衰减比Ignatyuk慢"的试错调整——但没有发表一组完整的阻尼常数灵敏度研究。这就是该领域当前的软肋：大家都同意壳效应比预期的活得久，但没有人为"到底活多久"给出一个模型无关的、误差棒完整的答案。

中国的参与在这个子领域相对有限。目前没有中国团队主导的五维或多维朗之万裂变动力学研究发表在高影响力国际期刊上——这是全球朗之万裂变学界的一个特征性缺口。日本（JAEA/有友嘉一、西尾胜久；Science Tokyo/石塚知香子、千葉敏）和欧洲（GSI/施密特、CENBG/茹拉多、法国CEA/施密特）是主要驱动力。有部分中国研究者在做基于朗之万方程的超重核熔合动力学（如近物所、中国原子能院的工作）——但熔合不是裂变，熔合的朗之万用的是核间距离而非集体形变坐标，物理数学有重叠但不完全相同。

这并不是说中国缺少相关的实验能力。北京串列加速器（BRIF）、即将建成的HIAF和惠州的大科学中心——这些装置能产生轻锕系和亚铅区的裂变核素，直接为朗之万模型提供校准数据。但理论一端的自主性目前仍主要在日本和欧洲。

最后——也是最具前瞻性的一点——朗之万动力学正在向全微观方向演化。目前的朗之万模型全是"宏观-微观"混合体：势能面来自液滴+Strutinsky壳修正，惯性来自不可压缩流体近似，摩擦来自墙-窗几何模型。它的漂亮在于简洁——五维就能重现十几种核的碎片质量分布、动能分布、中子多重数——但它的局限也在于简洁：宏观-微观模型不能告诉你哪些特定的单粒子轨道贡献了壳效应，也不能预言某种新型集体运动（如对振动模参与颈断裂动力学）的出现。

全微观的替代方案是含时密度泛函理论——TDDFT，或者超流版本TD-BCS。这些方法从有效核子-核子相互作用出发，在不引入任何宏观参数的情况下模拟裂变动力学——但它们目前的弱点恰恰是朗之万的优势：TDDFT不能自然地产生碎片质量分布的宽度（平均场体系太"确定论"了，波动不够），也不能长时间演化（计算成本太高，模拟一两个断裂事件就耗尽了一台超级计算机）。朗之万方程用随机力项的统计系综模拟弥补了这些缺陷——平均场给了势能面的骨架，随机力给了分布的宽度。

从2024-2025年的趋势看，这两个范式正在融合。Ivanyuk等人的五维朗之万框架已经用了微观壳修正（Strutinsky+折叠汤川势）而不是纯液滴——这是从宏观端向微观端迈出的第一步。Stochastic TDDFT（含时密度泛函+随机力项，也称随机平均场理论）正在从微观端向宏观端迈出相同的距离——在平均场方程中加了一个满足涨落-耗散定理的噪声源。两条路最终会在同一个位置相遇：一组从微观核子-核子相互作用出发、能同时描述集体运动（漂移项）、量子涨落（扩散项）、耗散（摩擦项）、以及温度依赖性壳效应（阻尼函数）的统一方程。

那是裂变理论的终极梦想——它还没来。但五维朗之万告诉你，它比你以为的更近了。

---

参考文献

  1. F.A. Ivanyuk, C. Ishizuka, S. Chiba, "Five-dimensional Langevin approach to fission of atomic nuclei," *Physical Review C* 109, 034602 (2024). DOI: 10.1103/PhysRevC.109.034602. arXiv:2310.19466. — 五维朗之万裂变模型的奠基性论文，首次引入颈参数ε作为独立动力学变量，重现钍同位素三峰质量分布。Ivanyuk（乌克兰科学院核研究所），Ishizuka（Science Tokyo），Chiba（NAT公司）。
  1. F.A. Ivanyuk, C. Schmitt, C. Ishizuka, S. Chiba, "Shell effects and multi-chance fission in the sub-lead region," *Physical Review C* 111, 054620 (2025). DOI: 10.1103/PhysRevC.111.054620. arXiv:2501.05362. — 五维朗之万应用于亚铅区¹⁸⁰Hg/¹⁹⁰Hg，发现壳效应存活至40-50 MeV，获编辑推荐。Ivanyuk（基辅），Schmitt（CNRS/IPHC斯特拉斯堡），Ishizuka（Science Tokyo），Chiba（NAT）。
  1. C. Ishizuka, K. Nishio et al. (JAEA), "Six-dimensional Langevin approach to fission with Cassini shape parameterization," 2024 Symposium on Nuclear Data, JAEA, Tokai, November 2024. — 首次六维朗之万计算，Cassini参数化{α, α₁, α₃, α₄, α₆}，²⁵⁸Md中发现三共存裂变模式。石塚知香子（Science Tokyo），西尾胜久（JAEA）。
  1. Y. Aritomo, K. Hagino, K. Nishio, S. Chiba, "Dynamical approach to heavy-ion induced fission using actinide target nuclei at energies around the Coulomb barrier," *Physical Review C* 85, 044614 (2012); 后续2024-2025年JAEA朗之万工作延续该框架。 — JAEA朗之万裂变动力学的标准参考。有友嘉一（JAEA/JINR），萩野浩一（东北大学），西尾胜久（JAEA），千葉敏（NAT）。
  1. S. Takagi, S. Harada, Y. Aritomo, K. Hirose, K. Nishio, "Neutron emission during fission and its impact on fission-fragment mass distribution studied by Langevin model," arXiv:2510.18205 (2025). — 首次将中子蒸发动态嵌入朗之万轨迹演化，中子发射"复活"壳结构。高木祥太、西尾胜久等，JAEA。
  1. Y. Huang et al., "Multimodality of ¹⁸⁷Ir fission studied by the Langevin approach," *Physical Review C* 109, 034609 (2024). DOI: 10.1103/PhysRevC.109.034609. — 三维朗之万在¹⁸⁷Ir中发现四共存裂变模式，E_d=28 MeV壳修正阻尼。
  1. Phys. Rev. C 112, 014626 (2025), "Dynamical simulation of the fission process of ²³³Th at low excitation energies in the framework of the four- and five-dimensional Langevin models." DOI: 10.1103/mjls-kymk. — 四维vs五维朗之万比较，⁵D显著优于⁴D，耗散系数k_s≈0.29。
  1. Phys. Rev. C 112, 034613 (2025), "Influence of transport coefficients on fission dynamics in the five-dimensional Langevin approach." — 宏观vs微观惯性张量对五维朗之万的敏感度研究，"热裂变模式"A≈134发现。JAEA。
  1. A.V. Ignatyuk, G.N. Smirenkin, A.S. Tishin, "Phenomenological description of energy dependence of the level density parameter," *Soviet Journal of Nuclear Physics* 21, 255 (1975). — Ignatyuk壳修正阻尼公式Φ=exp(−E*/E_d)的原始出处。
  1. V.M. Strutinsky, "Shell effects in nuclear masses and deformation energies," *Nuclear Physics A* 95, 420 (1967). — Strutinsky壳修正方法的原始论文。基辅核研究所。
  1. K.-H. Schmidt, B. Jurado, C. Amouroux, C. Schmitt, "General Description of Fission Observables: GEF Model Code," *Nuclear Data Sheets* 131, 107 (2016). — GEF裂变产额模型，与朗之万方法互补。Schmidt（GSI/CENBG），Jurado（CENBG波尔多）。
  1. J. Sadhukhan et al., "Inferring fission timescales from prescission neutron multiplicity using a Langevin dynamical model," *Physical Review C* 110, 034604 (2024). DOI: 10.1103/PhysRevC.110.034604. — 一维朗之万提取全局耗散强度系统学，识别A<210/A~210-244/A>244三个质量区。
  1. S. Ayik, M. Arik, O. Yilmaz, A.S. Umar, "Stochastic Mean-Field Theory for Nuclear Collective Dynamics," arXiv:2509.18914 (2025). — 随机平均场理论综述，从微观端向宏观朗之万融合的理论桥梁。Ayik（田纳西科技大学），Umar（范德堡大学）。
  1. V.M. Kolomietz, S.V. Radionov, "Stochastic features of nuclear fission," arXiv:2111.02487 (2024修订). — 非马尔可夫朗之万方程在裂变动力学中的应用，记忆效应修正Kramers率。Kolomietz（基辅核研究所）。