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

核裂变锯齿之谜-量子纠缠解开六十年难题

2026-06-24

1962年，洛斯阿拉莫斯的James Terrell在《物理评论》上画了一条锯齿形的曲线。他横轴画碎片质量，纵轴画碎片平均放出几个中子。对铀-235的裂变——轻碎片会吐出两三个中子，重碎片只吐一个左右；在对称裂变点（两个碎片差不多大），中子数又突然弹起来。像锯条上的尖齿。

Terrell给出了一个经典解释。他说，裂变碎片里接近幻数（N=50, Z=50, N=82）的那些核——因为有闭合的核子壳层保护——特别硬、特别圆。它们吸走不了多少变形能。它们的同伴、那个互补的碎片呢？接过了全部的变形——软、热、剧烈抖动——于是蒸发出更多中子。锯齿画出来了，它是两个碎片之间的"形变失配"造成的。

这个解释优美、自洽，用了六十年。只有一个问题：它是错的。

不是因果全错了——幻数效应导致了不对称裂变这件事是对的。错的是那条锯齿本身的起源机制。它不来自碎片自己长什么样、圆不圆、壳层闭不闭合。它来自一个更根本、更反直觉的原因：两块碎片在断裂的那一刻，它们的量子波函数来不及分开。

2024年8月，北京大学物理学院的强雨和裴俊琛——合作者包括密歇根州立大学FRIB的Kyle Godbey——在arXiv上投了一篇论文。2025年2月，它正式发表在*Physics Letters B* 861, 139248上。标题直白：*Quantum Entanglement in Nuclear Fission*——核裂变中的量子纠缠。

北京大学核物理与核技术国家重点实验室在同一天发布了新闻稿，标题是"裴俊琛课题组揭示核裂变中存在量子纠缠"。强雨是该组的博雅博士后。

他们干了什么？他们用含时哈特里-福克加BCS（TD-BCS）方法——一种从核子-核子相互作用出发模拟核动力学演化微观理论——模拟了钚-240从拉长到断裂全过程的波函数。核心创新是粒子数投影：在整个动力学演化结束后，用双粒子数投影把波函数投影到固定碎片核子数的分量上，由此算出每个碎片的激发能——也就是它能发射几个中子。

在此之前，你没法做这件事。你可以在演化过程的某一帧"冻结"形状，用经典几何变形能来判断哪个碎片更激发——但那只是模拟，不包含量子关联。强雨和裴俊琛做的是：在断裂后的任何时刻，把两碎片的联合波函数拆成"轻碎片有M个中子、N个质子，重碎片有……"这样的叠加态。这个叠加态本身就是一个纠缠态——碎片甲和碎片乙的核子数分配被波函数绑在了一起。

关键来了。如果把波函数做"纯态"处理——假设两块碎片各自独立，每个有确定的核子数——锯齿图案消失了。算出来的碎片激发能分布是平的。他们试了几种消除纠缠的方法：去掉交叉项、强制波函数坍缩到某一个核子数分量、用Schmidt分解只留最大特征值——每次消除纠缠，锯齿就消失。每次保留纠缠，锯齿就回来。

这不是相关性——这是因果性。锯齿之所以存在，是因为量子纠缠存在。

物理图像是这样的。一个铀核从静态开始，被一个中子打中——变形——拉长成花生状——颈部越来越细——最终断裂。断裂本身发生在约10⁻²²秒的时间尺度。在这个速度下，波函数根本追不上形状变化。断裂前的波函数是一个在整个拉长的核上非局域化的函数——它的核子可以同时属于"将来会成为轻碎片的那半边"和"将来会成为重碎片的那半边"。断裂后，波函数没有"决定"每个粒子到底去哪一边——它继续存在于一个叠加态，其中轻碎片"可能有Z=34个质子也可能有Z=36个质子"，重碎片相应互补。

这种叠加不是数学玩具。粒子数投影给出的分布宽度证明了它确实存在——轻碎片中质子数的分布明显宽于经典模型预测。强雨和裴俊琛的论文中，用TD-BCS+PNP（粒子数投影后处理）算出来的²⁴⁰Pu碎片激发能分布自然地重现了实验锯齿：轻碎片激发能约15-20 MeV，重碎片约8-12 MeV；对称裂变点处的激发能也在两者之间。去掉量子纠缠后，两个碎片激发能几乎相等，锯齿夷为平地。

2025年同组年轻学者尚红宇（H.Y. Shang）作为第一作者发表的跟进论文（*Nuclear Science and Techniques* 36, 211, 2025，裴俊琛通讯）进一步把分析推广到²³⁸U，并细致研究了超流配对强度如何调制纠缠度和锯齿斜率。发现配对越强，碎片越不对称（趋向更极端的质量分裂）——因为配对能抑制碎片表面的激发，放大了壳效应驱动的非对称性。弱配对或纯TDHF（完全去掉超流）算出的锯齿斜率太陡——比实验值高三到四倍。强配对使锯齿展平——接近实验斜率。这提供了一个额外的自洽性检验：实验已知裂变核的配对强度大致在什么范围，而在这个范围内，纠缠模型给出大致正确的锯齿。不需要人为调节碎片温度——那正是Terrell时代统计模型不得不做的"特设"假设。

那Terrell的经典解释错在哪？不是他观察错了——他观察的锯齿形状是真实存在的实验事实。他错把锯齿归因于碎片个体属性（圆 vs 扁、壳层闭合 vs 开放），而忽略了碎片之间的量子关联。两个碎片不是先各自决定自己的形变、然后各自蒸发出对应数量的中子。它们是先作为一个整体——一个纠缠的量子多体系统——把激发能分配给两半，然后各自带着分配到的激发能开始各自的蒸发级联。

这个区别不是哲学——它有可测量的后果。如果碎片激发能仅由自身壳层闭合决定，那么给定弹核和出射道，你预言的中子多重数应该与裂变道（不对称度）之间存在一个固定的函数关系——由碎片壳层刚性唯一确定。如果激发能由纠缠度决定，它除了依赖壳层——还依赖断裂速度、颈部动力学、核黏度——这些才是决定纠缠能"分离到什么程度"的微观参数。这意味着同一对碎片在不同能量的裂变中可能展示不同的中子发射数——因为断裂速度不同，纠缠保留程度不同。这件事在实验上似乎已经开始被看到了——虽然证据还在积累中。

还有一个被纠缠效应改变的结论是关于稀有同位素束流的。FRIB、RIBF、HIAF这些装置造不稳定核靠的是裂变——一束重离子打靶产生丰中子裂变碎片。传统理解是，给定元素的较重同位素（中子数更多）由于更靠近稳定线（较少富中子）应该产量更高。尚红宇-强雨-裴俊琛2025年的论文指出了一个抵消效应：同位素碎片的激发能随中子数增加而升高——也就是说，你需要的中子丰同位素恰恰是激发能最高的那批——它们蒸发的中子比其他同位素多，产物就比预期少。这是量子纠缠通过影响蒸发链间接调制了稀有同位素束流的强度分布。定量上，这个效应在海平面到中度束流能量范围可能贡献10-20%的产额修正——不大不小，足够需要被考虑。

还有另一条平行的科学线索没有使用"量子纠缠"这个词——但它的物理与强雨和裴俊琛的结果高度一致。

2025年，P. Morfouace（法国CEA）、J. Taieb（法国CEA）、A. Chatillon（法国CEA）等人领导的R3B合作组在*Nature* 641, 339-344 (2025)上发表了一篇重磅实验论文。他们在GSI达姆施塔特用相对论能量下的库仑诱发裂变（R3B-SOFIA装置），测量了100种缺中子核体系的裂变碎片电荷分布——其中75种从未被测量过。GSI的高能碎裂装置能产生超不稳定的核——比如汞-180——然后用电磁飞行的办法在它们衰变之前拿到碎片电荷分辨。结果发现了亚铅核素区（轻于铅的核素）一个由轻碎片Z=36（氪的同位素）质子壳层驱动的不对称裂变"岛"。这是一个全新的发现——不像锕系裂变（铀、钚）的不对称来自重碎片方的Z=52-56壳层，这里的驱动力来自Z=36的质子壳闭合和相应的形变效应。

这和量子纠缠有什么关系？R3B的实验数据表明，这些轻核的不对称裂变对中子数非常敏感——多一个中子或少一个，不对称度可能陡变。在统计模型中，这种敏感度很难复现——你需要碎片的"记忆"比你想象的好。壳层效应在快断裂中存活的机制——本质上就是波函数的非局域化，也就是纠缠。只是R3B的作者用了不同的语言来描述它。

物理的深层一致性就在这里：无论你从第一原理做TD-BCS+粒子数投影（北大团队），还是从实验数据做壳层效应的系统学（R3B/GSI），你把数据推到够精细——都撞到了同一个量子多体物理的核心：核裂变的断裂过程太暴烈、太快，波函数没有被给到时间去做选择。它把所有可能的结果打包在一起扔了出去。你要理解碎片——你必须理解纠缠。

一个迷人的开放问题是：这个纠缠有多大？强雨等人用Schmidt数（纠缠度的一种度量）来定量——它大约在2到5之间不等，取决于碎片的核子数差异和具体的断裂条件。这意味着波函数中占主导的核子数分量只有几个——不是成百上千个。这不是一种"高度纠缠"——不是多体局域化极限、不是量子混沌——而是一种由少数叠加态构成的弱纠缠。恰恰因为它弱，它对可观测量（激发能、中子数）的影响才能被分离出来。太强的纠缠——所有可能态等权叠加——你算出的激发能分布是平滑的，锯齿也会消失。

这是最微妙的一点。物理上能产生锯齿的纠缠必须"恰到好处"——不能太弱（锯齿消失），不能太强（锯齿也被抹平），只能在由壳层能量尺度和配对能强度共同决定的"窗口"内。这个窗口恰好对应自然界的裂变核。没有微调，是自然的巧合——但它确实很巧。

现在来诚实地说说局限。目前纠缠解释只在宏观-微观模型的框架内——用TD-BCS描述平均场演化，用粒子数投影恢复部分量子关联。它还不是全从头算——不是从手征有效场论的核力出发在多体希尔伯特空间做精确对角化。全从头算裂变动力学在今天是不可能的——计算复杂度高出任何现有超级计算机的能力。TD-BCS是一个最好的近似——但它仍然是一种平均场理论，不能自动包含所有纠缠。粒子数投影只恢复了由于粒子数对称性破缺导致的纠缠——不是全部。其他类型的纠缠——角动量纠缠、形状纠缠（不同碎片形变的叠加）——可能也贡献了一些，但目前还没有被定量研究。

还有一个不能回避的张力：不同理论组之间的对比基准严重不足。目前全球只有极少数组能做这种含纠缠效应的裂变动力学计算——北大裴俊琛组是其中之一，法国CEA的N. Schunck和D. Regnier组在做含高斯过程回归的TDDFT裂变，美国的G.F. Bertsch和W. Younes在做基于密度泛函的朗之万裂变。但这些组用的是不同的近似、不同的有效相互作用（SkM* vs UNEDF1 vs SLy4）、不同的粒子数投影实现方式。没有过一次联合基准测试——大家各自发表各自算出来的锯齿，互相对照。这个字段默认"不同模型算出来大致差不多"——但实际上不同组对同一对裂变碎片预言的锯齿斜率可以差30-50%。

最后一个诚实的注脚。Terrell在1962年提出的壳层硬化解释不是"全错"了——它仍然是锯齿图像的一个组成部分。壳层闭合决定了形变在碎片之间的不对称分配——这给了锯齿一个基本方向（轻碎片激发多、重碎片激发少）。裴俊琛组的发现是叠加在这层物理之上：壳层给了方向，量子纠缠给了强度。你如果把波函数里纠缠去掉——方向还在（因为形变不对称由壳层决定），但强度的精确剖面（锯齿的坡度）就消失了。两者加在一起才拼出完整的锯齿。

那么接下来的路怎么走？三条。

第一，全从头算的黎明。FRIB和未来EIC产生的数据压力会让理论社区加速发展全从头算裂变——即使只能做低维截断（比如约束在特定集体坐标上）。量子蒙特卡罗和耦合簇方法已经做到了中等核的基态能——把同一套第一原理框架推广到裂变动力学是"早晚问题"——可能还要十年。

第二，纠缠的实验证实。目前所有的"纠缠产生锯齿"证据都是理论内推的——从计算中去掉纠缠看锯齿消失，再把纠缠装回去看锯齿回来。这很强，但不是直接测量。有没有可能设计一个实验，直接探测两碎片之间的量子关联？这听起来像科幻——但在量子光学里你是可以直接测量纠缠保真度的。在核物理里很难——你在跟一个约10⁻²¹秒的暴烈事件打交道，而且它发生在你够不到的真空室深处。也许通过关联中子角分布或者关联伽马射线谱可以做间接探针。目前还没有这样的实验——这是空白。

第三，实际应用。如果量子纠缠确实影响激发能分配，它应该被写进下一代裂变产额评价数据库——ENDF/B-VIII、JENDL-5、CENDL-3.2。当前这些数据库的裂变模型基础仍是统计模型（Madland-Nix或GEF），不含纠缠修正。一个"纠缠知情"的裂变产额预测器可能将少数特定反应道（如²³⁵U热中子裂变）的中子多重数预言精度提高5-10%。对反应堆设计和核廢料嬗变——这在工程设计层面有意义——不是学术细节。

2025年，裂变物理这个看似老旧的领域正经历一场安静的范式革命。六十年来大家认为不需要量子纠缠来解释锯齿——毕竟中子发射本身是经典的蒸发过程，你只需要知道碎片一开始有多热。但2024-2025年的工作证明：碎片有多热——这个"初始条件"——不是由经典几何决定，而是由量子绑定决定。两个碎片不知道自己各自得到了什么——直到波函数被观测。每次你测量一对碎片的激发能，你都在见证一次纠缠坍缩。

核裂变不是量子纠缠的一个边角例证。它是强相互作用多体系统中可以进入的最暴烈、最复杂、最高能标的一个纠缠实验室。它和量子光学的纠缠共享同一个数学骨架——Schmidt分解、约化密度矩阵、冯·诺依曼熵——只是从一个完全不同的物理维度逼近了同一个核心原理：当你让一个量子多体系统裂成两半，它的波函数不是被"切"成两半——它是"展开"成所有可能切片方式的叠加。

这就是锯齿的量子根源。

---

参考文献

  1. Y. Qiang, J.C. Pei, K. Godbey, "Quantum Entanglement in Nuclear Fission," *Physics Letters B* 861, 139248 (2025). DOI: 10.1016/j.physletb.2025.139248. arXiv:2408.07375. — 首次揭示核裂变中存在量子纠缠，证明它是锯齿结构的关键原因。第一作者强雨（北京大学"博雅"博士后），通讯作者裴俊琛（北京大学物理学院、核物理与核技术国家重点实验室），Kyle Godbey（密歇根州立大学/FRIB）。
  1. H.Y. Shang, Y. Qiang, J.C. Pei, "Energy partition between entangled fission fragments," *Nuclear Science and Techniques* 36, 211 (2025). DOI: 10.1007/s41365-025-01782-5. arXiv:2501.07118. — 研究纠缠态碎片间的能量分配，对关联对锯齿斜率的影响，以及对稀有同位素束流产生的含义。第一作者尚红宇，通讯作者裴俊琛。
  1. J. Terrell, "Neutron Yields from Individual Fission Fragments," *Physical Review* 127, 880 (1962). DOI: 10.1103/PhysRev.127.880. — 经典锯齿曲线奠基论文，提出碎片硬化（幻数效应）解释。James Terrell，洛斯阿拉莫斯国家实验室。
  1. A.C. Wahl, A.E. Norris, R.A. Rouse, J.C. Williams, "Products from Thermal-Neutron-Induced Fission of ²³⁵U," *IAEA Symposium on Physics and Chemistry of Fission*, Vienna, p. 813 (1969). — Wahl锯齿曲线，放射化学与物理方法联合分析。A.C. Wahl，华盛顿大学。
  1. 北京大学核物理与核技术国家重点实验室，"裴俊琛课题组揭示核裂变中存在量子纠缠，" 2025年1月22日新闻稿。 https://sklnpt.pku.edu.cn/info/1122/3236.htm — 官方新闻发布，确认量子纠缠解释锯齿的发现。
  1. C. Simenel, A.S. Umar, "Formation and dynamics of fission fragments," *Physical Review C* 89, 031601(R) (2014). — TDHF裂变动力学的现代计算基础。Cédric Simenel，澳大利亚国立大学。
  1. G. Scamps, C. Simenel, D. Lacroix, "Dynamical description of the fission process using the TD-BCS theory," *Physical Review C* 92, 011602(R) (2015). — TD-BCS裂变动力学奠基性工作。G. Scamps，法国CEA。
  1. J.C. Pei, W. Nazarewicz, J.A. Sheikh, A.K. Kerman, "Fission Barriers of Compound Superheavy Nuclei," *Physical Review Letters* 102, 192501 (2009). — 裴俊琛超重核裂变位垒的早期工作。裴俊琛（北京大学/当时在田纳西大学/ORNL），W. Nazarewicz（田纳西大学/ORNL）。
  1. P. Morfouace, J. Taieb, A. Chatillon et al. (R3B Collaboration), "An asymmetric fission island driven by shell effects in light fragments," *Nature* 641, 339-344 (2025). DOI: 10.1038/s41586-025-08882-7. — GSI/R3B-SOFIA实验：100种缺中子核体系的裂变碎片电荷分布，首次发现亚铅区不对称裂变岛。第一作者P. Morfouace，法国CEA/DAM/DIF。
  1. P. Möller, A.J. Sierk, T. Ichikawa, A. Iwamoto, M. Mumpower, "Fission barriers at the end of the chart of the nuclides," *Physical Review C* 91, 024310 (2015). — Möller裂变位垒系统学，布朗形变模型(Brownian shape motion)，不对称裂变模式的壳层驱动图像。
  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裂变产额模型的系统学描述。K.-H. Schmidt（GSI/CENBG），B. Jurado（CENBG）。
  1. A. Bulgac, P. Magierski, K.J. Roche, I. Stetcu, "Induced Fission of ²⁴⁰Pu within a Real-Time Microscopic Framework," *Physical Review Letters* 116, 122504 (2016). — 首次TDDFT从头算超流裂变动力学模拟。A. Bulgac，华盛顿大学。
  1. Y. Qiang, J.C. Pei, M. Kortelainen, "Fission dynamics of compound nuclei: Pairing versus fluctuations," *Physical Review C* 103, 024612 (2021). — TD-BCS裂变动力学的先前工作，²⁴⁰Pu不对称道S1和S2的识别。
  1. 核物理中的量子纠缠效应研讨会，复旦大学，2025年6月。 https://imp.fudan.edu.cn/48/31/c48933a739377/page.htm — 裴俊琛做专题报告：核裂变中的量子纠缠与锯齿结构。