贝叶斯神经网络计算核裂变碎片产额

作者丨庞龙刚

单位丨华中师范大学

研究方向丨高能核物理、人工智能

今天介绍一篇北京大学物理系使用贝叶斯神经网络计算核裂变碎片产额的文章。这篇文章发表在 PRL 上，业内同行都很感兴趣。这里对我们大同行的文章做个解读，通过费曼学习法，研究一下什么是贝叶斯神经网络。

什么是核裂变碎片产额

一个重的原子核裂变成两个或多个轻原子核，并释放能量的过程，都可以称作核裂变。裂变产物 -- 轻原子核，叫裂变碎片。如果母核 A 有 230 个核子，子核 B 有 115 个核子，子核 C 有 115 个核子，则此裂变为对称裂变，裂变产物分布为 N=115 的地方有个高峰。如果子核 A 有 90 个核子，子核 B 有 140 个核子，则此裂变为非对称裂变。裂变产物分布在 N=90 和 N=140 处有两个峰。裂变碎片产额分布如下图左边子图所示。

▲ 文章配图例子。横轴是衰变产物中核子的个数（Mass Number），纵轴是裂变产额。

通过理论计算或模型预言重原子核的裂变碎片产额分布，并与实验数据对比，能够帮助排除不符合实验的裂变动力学模型。

如果理论结算结果与实验测量的裂变碎片产额分布一致，则说明理论可能是对的，就可拿理论去预测更多的量。

原子核的质量与裂变

这一段面向大众科普一下核物理。我们先从爱因斯坦的质能关系出发：

其中 E 是能量，m 是质量，c是光速，在自然单位制下，c=1，公式可以简化为 E = m。

这个公式告诉我们能量就是质量。核裂变要用到这个公式 -- 如果一个大原子核 A 裂变成两个小原子核 B+C，则释放的能量等于裂变前后的质量差。

但是为何 A 变成 B 和 C 的时候总质量会发生变化呢？明明 A 中的质子数等于 B 和 C 中质子数之和，A 中的中子数等于 B 和 C 中中子数之和。 

这是因为 m = E，质量也是能量。原子核的质量等于里面质子质量+中子质量+质子，中子动能+质子质子，质子中子，中子中子的相互作用势能。 

另一个简单的例子是质子质量，一个质子里有 3 个组分夸克，2 u + 1 d。一个 u 夸克质量 2.16 兆电子伏特左右，一个 d 夸克质量 4.67 兆电子伏特，总质量应该是 2.16 x 2 + 4.67 = 8.99 兆电子伏特。但测量的质子质量是 938 兆电子伏特。剩余的质量全部来自夸克，胶子的动能，以及它们之间的相互作用势能。 

大自然目前发现 4 种相互作用力，1）万有引力；2）库仑力；3）强相互作用力（核力)；4）弱相互作用力。除了万有引力不贡献核子核子之间的相互作用势能，其他三种力都有贡献。

库仑力贡献了质子质子之间的电磁相互作用势能。强相互作用贡献了核子核子相互作用的一部分。弱相互作用大家可能比较陌生，它发生在基本粒子味道改变的过程，比如中子衰变成质子的过程，中子 (夸克组分 udd) 衰变成质子（夸克组分 uud）时，一个 d 夸克变成了 u 夸克，起作用的就是弱相互作用。 

知道了爱因斯坦质量能量关系，就可以预测一下不同原子核的质量了，如果一个原子核有 Z 个质子，N 个中子，则总的核子数 A=Z+N。在宏观液滴模型下，原子核的质量大概可以写成如下形式，里面包含参数化的各种各样的能量对原子核质量的贡献（当然，下列公式没有包括微观的来自库伯对和自旋轨道耦合的壳修正项）。

▲ 参考文献：Nuclear gro