|作者：裴俊琛 强雨 乔春源
(北京大学物理学院 核物理与核技能国家重点实行室)
本文选自《物理》2022年第11期
择要 核裂变的发现深刻地影响了人类社会。核裂变的研究还在不断深入，一方面核裂变有新的应用需求，另一方面核裂变是一个复杂的量子多体动力学过程。比年来，核裂变理论和实行研究有很大希望，人们对核裂变几率、裂变产物和裂变机制都获得了新的熟悉，这有助于澄清一些唯象模子的经验假设。此外，机器学习的应用为发掘使用不精确不完备的核数据提供了可能。等待未来更精确更自洽的核裂变理论可以更好地支撑应用创新。
关键词 核裂变，先辈核能，裂变机制
1. 引 言
1939年2月，Meitner与Frisch初次展现了铀原子核像液滴一样发生了分裂，并用fission这个词来形貌核裂变。更紧张的是，他们基于玻尔的液滴模子估算出一次核裂变会释放约200 MeV的能量。现实上fission一词最早是指生物学中的细胞分裂。核裂变释放的能量是云云之巨大，很快就引起了科学家们的极大爱好。1942年12月，费米在芝加哥大学实现了可控的链式核裂变反应，开启了和平使用原子能的时代。1945年7月，美国乐成爆炸了第一颗原子弹，深刻地改变了人类汗青。核裂变的发现是一个曲折的传奇故事，一些大科学家曾与之失之交臂，它生动地展现了科学熟悉积累到一定程度后灵光一现的思想突破。
核裂变的发现至今已经80多年，它深刻地影响了人类社会。人们猜测地球内核就是一个巨大的核裂变反应堆，不停保持着人类生存的温暖。核裂变一方面会释放巨大的能量造福人类，另一方面假如控制不好会带来劫难性的影响。这需要我们进一步研究核裂变，更精准地熟悉核裂变，更好地使用核裂变。北京大学胡济民先生所著的《核裂变物理》一书对核裂变研究举行了全面体系的阐述。比年来，核裂变的理论和实行研究取得了显著希望，产生了一些新的熟悉，这为核裂变的应用带来了新的可能。
2. 为什么我们还要研究核裂变？
由于地球上U储量有限，发展先辈的可持续、更安全、更清洁的核裂变能将越来越紧张。目前主流的核能是压水堆，此中裂变产生的中子经过慢化后变成能量很低的热中子。第四代先辈核能(图1(a))的一个主流方向是快中子堆，快中子堆无需中子慢化剂，可以更紧凑。快中子堆可以通过增殖反应将U变成易裂变的Pu，将铀资源的使用率从1%提升到60%。同时快中子堆大幅度地减少了核废料的放射性寿命。相比于压水堆，发展快中子堆需要更精确的、中子能量连续的核裂变数据。而目前国际上主要核数据库的核裂变产物的产额只有热中子、0.5 MeV与14 MeV三个能量点的评价数据。更精确的核模子与核数据也有助于设计更细密紧凑的专用核动力(船用、月基、空间，图1(b))和更好的支撑国防研究。
图1 核裂变的新应用，包括先辈核能(a)、空间核动力(b)、同位素电池(c)、天体环境下 R-过程中的核裂变(d)、同位素药物(e)、超重元素的合成(f)、反应堆中微子的研究(g)等
除了使用核裂变释放的巨大能量以外，裂变产物核的循环使用将是一个巨大的机遇。通过核裂变产生的Mo可以获得Tc，Tc是用于核医学诊断的紧张同位素，已经有广泛成熟的应用(图1(e))。英国科学家使用核废料长期放射性的特点，通过钻石包裹制成了能够稳定供电两千多年的核电池，创新性地实现了变废为宝。使用反应堆还可以生产Pu，已经将其制成同位素电池用于中国的火星车和月球车，但是其生产还很昂贵(图1(c))。反应堆内的核反应网络十分复杂，这也为实现先辈核能提供了新的可能。通过核裂变可以产生数百种同位素核，大部门裂变产物核是不稳定的。通过加速核裂变碎片可以形成放射性束流，国际上新一代放射性束流装置的主要科学目标是研究极度条件下的奇特核物质，如美国的有数同位素束流装置(FRIB)，中国的强流重离子加速器装置(HIAF)等，这将极大地扩展核物理的研究范围。长命命放射性核素在反应堆内会大量积累起来，对反应堆设计、核废料处置惩罚、裂变产物循环使用都十分关键。
核裂变对一些紧张的根本题目，比如超重新元素的合成、宇宙中元素的演化过程、反应堆中微子的研究，也不可或缺。实行上熔合反应合成的超重核处于高引发态，它的存活概率取决于中子蒸发与裂变之间的竞争。实行上合成超重核极其困难，往往一年才观测到1—2个事例，需要可靠的理论指引(图1(f))。在双中子星并和与超新星发作的喷射物中，会发生R-过程快中子俘获反应，从而产生重元素(图1(d))。地球上的铀、钚都劈头于天体环境下的R-过程，但是R-过程到了极度丰中子超重核区由于裂变而停止，他们裂变的产物又循环参与R-过程从而显著地影响终极宇宙元素的丰度。此外，在核反应堆中裂变产物的β衰变会产生大量的反中微子(图1(g))，对其能谱的观测将展现一个基本物理题目，即是否存在第4种中微子——惰性中微子。这些新的应用和根本研究都依赖更可靠的核裂变的几率与产物产额。
图2 核裂变过程示意图。此中横轴是随着裂变核拉长的形变增加，红线表示裂变位垒，曲率表示位势的弯曲程度
核裂变虽然是一个老题目，但是从微观角度看，核裂变是一个极其复杂的非平衡非绝热的量子多体动力学过程，如图2所示。传统的唯象裂变模子通过引入一些参数，对实行数据较多的核区能较好地形貌，但是无助于深刻明白核裂变以及外推到实行很难到达的核区。原则上微观核裂变理论可以自洽地形貌多种裂变观丈量，但是微观模子离应用需求的精度尚有一定的间隔。比年来，随着超等盘算机的发展，科学家对核裂变机制获得了一些新熟悉。发展能形貌多种裂变观丈量，包括碎片产额、释放动能、释放γ光子数、释放中子数、裂变几率与裂变截面等观丈量的综合可靠的微观裂变理论是一个紧张科学目标，可以更深刻地明白核裂变过程，并对很难精确丈量的核数据和空缺核数据提供关键的增补。此外人工智能与机器学习的应用可以帮助我们更好地模拟核裂变和挖掘核数据。比年来实行上提供了亘古未有的精确的裂变观丈量，为进一步验证、束缚和发展新的裂变理论提供了机遇。
3. 核裂变几率
核裂变的寿命或者裂变几率是一个关键的观丈量。重核的自发裂变是裂变位垒的量子隧穿过程，这是一个十分缓慢的过程。裂变位垒是指原子核结合能随着核形状拉长而变化的曲线或多维曲面。原子核的多维团体形状空间由于量子壳效应而呈现复杂的裂变位垒。当原子核处于高引发态，量子效应(对关联、壳效应)渐渐消失，裂变几率可以由统计模子形貌，裂变寿命为10—10s。随着引发能增加，裂变寿命先是下降很陡，到了高引发时变成缓慢下降，所以裂变的机制是能量相干的。
传统的玻尔—惠勒统计模子，也叫过渡态理论，在核裂变寿命的盘算中有广泛的应用，但依赖较多唯象参数。为了形貌裂变位垒的能量相干性，唯象模子通常引入一个参数，来形貌壳修正能随引发能增加而指数衰减的因子。而高引发态的裂变位垒可以通过微观的有限温度的能量密度泛函盘算给出。微观盘算可以自洽地思量量子效应随温度增加而渐渐消失的过程。德国重离子研究中心(GSI)与日本的理化学研究所(RIKEN)通过冷熔合合成了107至113号元素，俄罗斯的杜布纳联合核子研究所(Dubna)通过热熔合合成了114至118号元素。我们的微观盘算结果表明，超重复合核的裂变位垒随引发能增加而减小的因子在差别核区十分差别。热熔合区比冷熔合区的裂变位垒下降要慢，在高引发时仍有一定的位垒，这是唯象模子所没有预料到的。相对于冷熔合所呈现的趋势来说，通过热熔合生成超重元素有很大的截面，这在当时是一个让人狐疑和质疑的题目。我们在2009年的理论工作对澄清和推动后续热熔合实行做出了贡献。此外微观盘算可以形貌位垒形状随能量的变化，热熔合区复合核在高引发时的位谷曲率比一般重核的曲率要小4—5倍，而裂变几率正比于位谷曲率。总之基于微观盘算，热熔合超重复合核有显著的存活概率，展现了114—118号超重元素合成的关键因素。我们理论预言合成119、120号新元素的存活概率与118号相似。此外，我们通过微观盘算的复合核裂变位垒可以解释兰州近代物理所的丰质子重核区的新核素实行。
基于微观盘算的裂变位垒的能量相干性也可以解释实行观测的裂变产物分布的能量相干性。随着引发能增加，裂变模式从不对称裂变渐渐演化到对称裂变。如铀、钚的第二个裂变位垒的高度在引入反射不对称形变后可以低沉2.5 MeV，说明低引发时由不对称裂变主导，但是它们的差异随着引发能增加而渐渐减少。
在热浴环境下，核裂变几率可以用虚自由能法(ImF)来盘算。该方法在化学反应中有广泛的应用。我们推导了玻尔—惠勒模子与ImF方法之间的联系，发现它们的主要区别是位垒与位谷的能级密度参数的区别。当处于极高的温度时，量子壳效应消失，玻尔—惠勒模子的裂变几率与虚自由能法只相差一个因子，即。统计模子十分依赖位谷与位垒的能级密度。原子核的能级密度与核的形状和引发能有关，微观盘算能级密度还是一个很大的挑衅。随着复合核引发能的增加，可能出现发射中子后的裂变。这一定程度上反映了高温核物质有较大的耗散系数和粘滞性，会延缓裂变动力学过程。通过丈量裂变前中子发射多重数可以推测出高引发核的耗散系数在增加。
4. 核裂变产物
核裂变会形成很多差别碎片产物的组合，这些裂变产物核的产额分布是非常紧张的裂变观丈量。此外，裂变碎片会携带很大的动能，这是核能释放的主要形式。裂变碎片处于引发态，会迅速释放中子而冷却。差别裂变碎片释放的中子数、携带的动能和角动量也差别。往往轻质量碎片释放两个中子，而重碎片释放一个中子，这与我们的预想不一样。裂变碎片会通过β衰变形成终极的累积产额。裂变产物的产额与其他裂变后的观丈量是关联在一起的。
理论上形貌裂变产物的产额分布主要基于多维裂变位能面的形状演化，比如基于多维朗之万方程求解可以公道地形貌裂变碎片的质量分布。与花粉在水中的布朗运动相似，在朗之万方程中，核裂变是原子核在团体形变空间的缓慢演化，而核的单粒子运动作为随机配景在快速的变化，这是一个经典的动力学方程，思量了涨落—耗散效应。此外基于微观裂变位垒的时间相干的生成坐标法(TD-GCM)也能大抵形貌裂变产物的分布。TD-GCM主要是基于裂变位能势的驱动，原则上微观盘算多维位能面可以更公道地形貌裂变产额，但是盘算量极大。这些基于静态的裂变位垒的形状演化本质上是绝热动力学，不能自洽思量碎片引发。基于时间相干的密度泛函(TD-DFT)可以形貌裂变的非绝热非平衡的动力学过程。TD-DFT是基于微观的单粒子波函数随时间的自洽演化，不需要盘算裂变位能面。TD-DFT能自洽盘算多种裂变观丈量的平均性子，但TD-DFT长期存在的一个题目是无法给出足够展宽的分布，这是因为TD-DFT缺乏团体自由度的涨落。为了解决这一题目，我们提出了裂变过程中单粒子能级随机跃迁的图像。在TD-DFT中，单粒子运动与团体运动是交织在一起的，随着有用温度增加，随机跃迁的效应增加，经过长时间的演化累积而得到有展宽的分布。此前国际上提出了随机平均场模子来形貌产额分布，是基于很大的初始涨落，但是这与自发裂变抵牾。
图3 (a)基于Brosa模子对裂变产物分布的形貌；(b)微观计算的Pu核裂变动力学演化路径(此中单元b表示10 m)
唯象的Brosa模子从裂变产物的质量分布出发，认为存在两种不对称的裂变模式(图3(a))，并认为差别的不对称裂变模式的劈头是受到裂变位垒的影响。Brosa模子还通过颈部随机断裂来形貌裂变观丈量的展宽，脖子越长分布越宽。Brosa模子可以公道地解释裂变产额分布、总动能分布、中子发射多重数之间的关联，是物理直觉的很大乐成，但是不停缺乏微观理论的支持。我们的结果展现了动力学涨落效应正是Brosa模子中的S1、S2两种不对称裂变模式的劈头。如图3(b)所示，随着涨落增加，长脖子S2裂变道的成分在增加，这与实行是一致的。这两种模式的裂变路径相似，不大可能是静态位垒的影响。
图4 (a)基于贝叶斯机器学习对U的不完备裂变产额的评价；(b)对U裂变碎片Xe的产额—能量关系的评价
比年来人工智能与机器学习在很多学科中的交织应用获得了很大关注。实行上丈量的裂变产物的产额往往是不完备的，有误差或存在分歧。特殊是能量相干的裂变产额，在中子入射能量处于2 MeV与14 MeV之间的数据比力希罕。在这种配景下，我们提出基于贝叶斯机器学习来学习增补缺失的裂变产额数据(图4(a))，展示了机器学习的现实应用代价和上风。此外通过输入碎片的电荷奇偶信息，以及在学习中引入负值处罚等，将物理信息和物理束缚与机器学习举行了实行结合。近来我们提出通过数据融合来更好地评价不完备、有分歧、有误差的核裂变产额。在反应堆中，裂变产物Xe有很大的中子吸收截面，是反应堆的“毒物”，会显著低沉反应堆运行功率，其产额的评价很紧张。图4(b)是我们基于贝叶斯机器学习对Xe产额的评价。当一种裂变数据在某些能区很希罕时，它与别的数据在其他能区的关联有助于这种数据的推断。数据融合可以思量潜在的、高维的、非局域的关联，可以给出综合的误差传播，可以发掘出不精确的裂变数据的最大代价，有望形成新的核数据评价方法。
5. 核裂变机制
核裂变是一个极其复杂的非平衡非绝热的量子多体动力学过程，裂变后碎片之间存在量子胶葛。更深入地熟悉裂变机制有助于发展精确的裂变理论。TD-DFT理论最得当研究核裂变的微观机制。比年来，超等盘算机的应用为微观裂变动力学的发展提供了很好的机遇，使我们有可能澄清或者更新一些唯象的裂变模子的经验图像。
在TD-DFT的根本上，对关联是最紧张的剩余相互作用，对裂变机制有紧张影响。人们熟悉到对关联相称于裂变的“润滑剂”，可以加速裂变过程。静态对关联也会导致裂变位垒的稍微低沉，这会显著减少自发裂变的寿命。当核的密度缓慢变化时，动力学对关联有快速的涨落。当对关联非常大时，核体系形成一个超流的团体态，涨落效应被压制。在TD-DFT盘算中，当对关联弱时，可能出现三分裂核裂变，而对关联强时则是二分裂。当对关联很弱时，裂变路径往往走短脖子的S1裂变道，与实行不符。当复合核处于高温引发时，对关联会很快衰减，涨落也更迅速。
核物质的耗散系数或者粘滞系数到底有多大呢？在著名的沥青滴漏实行中，由于沥青有很大的粘滞性，形成一个脖子拉伸很长的液滴，约10年才滴出一滴，如图5(a)所示。我们认为核物质的粘滞性比沥青小，但是比水大。通过把微观动力学裂变路径映射到经典动力学方程，可以提取出形状相干的耗散系数。这个盘算需要提取出动力学的裂变位势，相比于静态位垒，动力学的裂变断点更远，碎片间的库仑能更小。这意味着相比于非绝热裂变，绝热盘算的总动能会显著偏大，这也验证了非绝热裂变的公道性。我们的结果表明，裂变过程中耗散系数也在发生变化，一般在2×10—4×10s，这有助于束缚唯象模子的耗散参数。耗散系数随着引发能增加而增加，随对关联增加而减小。随着引发能增加，由于很强的耗散，裂变动力学演化时间越来越长。这时涨落效应成为裂变的主要驱动机制，这与涨落—耗散定理是一致的。很强的耗散和粘滞性将导致一个拉长的裂变颈部构型(图5(b))，从而使库仑能低沉，导致裂变释放的总动能减少。由于能量守恒，总动能减小，导致碎片的引发能增加，从而发射更多中子。微观TD-DFT盘算能公道地解释裂变机制的能量相干性。
图5 (a)沥青滴漏实行展示出有拉长脖子的液滴；(b)微观TD-DFT盘算给出的Pu裂变的断点构型
裂变分裂成的两个碎片之间存在很强的动力学胶葛，主导着碎片之间的能量分配、核子数分配、角动量关联等。由于断裂的碎片很快飞开，以至于部门的动力学胶葛还来不及塌缩。图6分别展示了差别裂变碎片发射中子数的平均值随碎片质量变化的分布、碎片中的中子与质子的平均比值，以及差别碎片所携带的平均角动量的分布。碎片发射中子的几率主要由碎片的引发能决定。TD-DFT盘算得到的轻碎片引发能大于重碎片引发能，且它们的差异随引发能增加而减少。现实上碎片发射中子的多重数分布是一个型的锯齿结构(图6(a))，其能量分配机制尚有待微观理论的解释。裂变是颈部随机断裂还是量子胶葛主导呢？唯象Brosa模子认为是随机断裂主导的。基于碎片之间的胶葛，两个碎片的核子数分配和能量分配具有不确定性。微观盘算通过粒子数投影也能获得有展宽的分布。比年来法国实行组取得了很大的希望。实行上可以获得全部碎片的产额分布，此中碎片的中子与质子的比值也具有锯齿结构(图6(b))，但与中子多重数的锯齿结构相反，这为进一步熟悉裂变断裂机制提供了观丈量。此外近来实行上获得了裂变碎片的角动量分布(图6(c))，这是一种新的观丈量，引起了理论上很大的关注。美国一些理论组很快对裂变碎片角动量分布提出了多个解释。需要指出的是，美国在裂变理论模子方面有长期积累的上风。裂变碎片角动量的获得是与断裂模式(比如脖子的扭曲或弯曲断裂)有关呢，还是断裂后获得的？碎片的角动量分布也有相似的锯齿结构，可能主要由能量分配机制主导。原则上TD-DFT可以形貌多种裂变观丈量之间的关联，但是还需额外思量超越平均场的效应，终极形成一个综合、自洽、可靠的微观裂变理论。
图6 (a)实行上观测的裂变碎片的平均发射中子数的分布；(b)碎片中平均的中子/质子比；(c)碎片所携带的平均角动量
6. 总结与展望
核裂变是一团强关联的量子物质分裂成两块的独特的量子动力学过程。核裂变的发现至今已经有80多年，但是核裂变过程非常复杂，对它的熟悉尚有待进一步深入。裂变过程中既有单粒子自由度，也有团体自由度、团体自由度，尚有涨落—耗散效应等交织在一起。裂变断裂前的脖子构型对裂变后观丈量有紧张影响。裂变断点既有随机性，也存在碎片之间的量子胶葛。断裂前体系的裂变位垒、能级密度、耗散系数也具有能量相干性，此外对关联对裂变机制有紧张影响。我们看到微观的TD-DFT可以乐成地解释裂变机制，有助于澄清或更新一些唯象的裂变模子的图像。基于BBGKY框架，进一步思量更高阶的关联动力学，可以更现实地形貌裂变。随着盘算本领的增加，思量高阶关联的裂变动力学将是一个紧张方向。目前微观核裂变理论基于有用核力，存在一定的不确定性，包括裂变位垒的预言也存在误差。从现实核力出发，发展从头盘算(ab initio)核结构是核物理的前沿方向。基于ab initio盘算重核裂变过程还很遥远，但是可以为发展更精确的有用核力和有用哈密顿量提供指引。
无疑，核裂变的研究有很强的应用配景。为了应对气候变化，先辈核能的发展将受到更大的器重。发展更精确可靠的核裂变理论，对升级核能应用十分紧张，对一些巨大的根本研究也很关键，比如超重元素合成、天体环境中的R-过程、中微子研究等。核裂变的研究既是迷人的量子多体题目，也有很强的交织应用需求。随着超等盘算、机器学习、量子盘算的应用以及实行装置的发展，核裂变的根本研究迎来了新的机遇，将为核裂变应用提供新的线索。比年来，美国、法国等在核裂变根本研究方面取得系列希望，中国在核裂变的根本研究方面尚有很大的发展空间和前景。
<hr/>参考文献
[1] Meitner L，Frisch O R. Nature(London)，1939，143：239
[2] 胡济民. 核裂变物理学. 北京：北京大学出书社，2014
[3] Kolos K，Sobes V，Vogt R et al. Phys. Rev. Research，2022，4：021001
[4] 核物理与等离子体物理—学科前沿及发展战略 . 北京：科学出书社，2017
[5] http://hiaf.impcas.ac.cn/
[6] 周善贵. 物理，2014，43(12)：817
[7] Eichler M，Arcones A，Kelic A et al. Astrophys. J.，2015，808：30
[8] Mention G，Fechner M，Lasserre Th et al. Phys. Rev. D，2011，83：073006
[9] Bender M，Bernard R，Bertsch G et al. J. Phys. G，2020，47：113002
[10] Bohr N，Wheeler J A. Phys. Rev.，1939，56：426
[11] Pei J C，Nazarewicz W，Sheikh J A et al. Phys. Rev. Lett.，2009，102：192501
[12] Zhu Y，Pei J C. Phys. Rev. C，2016，94：024329
[13] Qiao C Y，Pei J C. Phys. Rev. C，2022，106：014608
[14] Ma L et al. Phys. Rev. C，2022，106：034316
[15] Pei J C，Zhu Y. Nucl. Phys. Rev.，2017，34(1)：87
[16] Liu L L，Wu X Z，Chen Y J et al. Phys. Rev. C，2019，99：044614
[17] Randrup J，Moller P，Sierk A J. Phys. Rev. C，2011，84：034613
[18] Regnier D，Dubray N，Schunck N et al. Phys. Rev. C，2016，93：054611
[19] Zhao J，Niki T，Vretenar D et al. Phys. Rev. C，2019，99：014618
[20] Lu B N，Zhao J，Zhao E G et al. Phys. Rev. C，2014，89：014323
[21] Negele J W，Koonin S E，Mller P et al. Phys. Rev. C，1978，17：1098
[22] Bulgac A，Jin S，Roche K J et al. Phys. Rev. C，2019，100：034615
[23] Qiang Y，Pei J C，span >Stevenson P D. Phys. Rev. C，2021，103：L031304
[24] Lacroix D，Ayik S. Eur. Phys. J. A，2014，50：95
[25] Brosa U，Grossmann S，Müller A. Phys. Rep.，1990，197：167
[26] Wang Z A，Pei J C，Liu Y et al. Phys. Rev. Lett.，2019，123：122501
[27] Qiao C Y，Pei J C，Wang Z A et al. Phys. Rev. C，2021，103：034621
[28] Wang Z A，Pei J C. Phys. Rev. C，2021，104：064608
[29] Wang Z A，Pei J C，Chen Y J et al. Phys. Rev. C，2022，106：L021304
[30] Qiang Y，Pei J C. Phys. Rev. C，2021，104：054604
[31] Schmidt K H，Estienne M，Fallot M et al. EPJ Conf.，2021，256：00015
[32] Martin J F et al. Phys. Rev. C，2021，104：044602
[33] Ramos D et al. Phys. Rev. Lett.，2019，123：092503
[34] Wilson J N et al. Nature，2021，590：566
[35] Verriere M，Schunck N，Regnier D. Phys. Rev. C，2021，103：054602
<hr/>转载内容仅代表作者观点
不代表中科院物理所立场
如需转载请联系原公众号
<hr/>来源：中国物理学会期刊网
编辑：草莓熊 |
