本节摘要:把原子核比作一小滴不可压缩的带电液滴,用表面张力、库仑排斥与对称能的竞争,就能得到半经验质量公式,定量解释结合能随质量数的平滑走向。液滴模型虽不能交代能级细结构,却能把“重核何以不稳定”“裂变何以重核才能发生”讲得清清楚楚。本节用 Python 数值实现贝特-魏茨泽克公式,画出结合能曲线并找铁峰。
阅读完本节,你应当能够:
原子核每核子结合能约 8 MeV、波动平稳,这与液体里每个分子周边邻居数大致恒定非常相似。于是物理学家把核当一小滴不可压缩、带电的液滴:总体积正比于 A,表面有一层张力,质子间的正电荷互相排斥,还有一条让质子中子不想一边倒的对称倾向。这套“液滴”直觉后来长成半经验质量公式,是核物理最早也最耐用的整体模型。
贝特-魏茨泽克公式把原子核结合能拆成五项:体积项(与 A 成正比,体现短程饱和)、表面项(修正表面粒子邻居少、能量略低)、库仑项(质子间排斥,负贡献,∝Z²/A^(1/3))、对称项(惩罚质子中子失衡,∝(N−Z)²/A)与配对项(偶偶核更稳)。系数由拟合大量实验质量定出。它不说每个能级的秘密,却在宏观走向上高度准确。
import numpy as np # 半经验质量公式(贝特-魏茨泽克) 结合能项, 单位 MeV aV=15.8; aS=18.3; aC=0.714; aSym=23.2; aP=12.0 def binding(A, Z): N=A-Z BE = aV*A - aS*A**(2/3) - aC*Z*Z*A**(-1/3) \ - aSym*(N-Z)**2/A # 配对项 if (A%2)==0 and (Z%2)==0: BE += aP*A**(-1/2) elif (A%2)==0 and (Z%2)==1: BE -= aP*A**(-1/2) return BE # 对每个 A, 找一个结合能最大的 Z(沿稳定线) peakBE=0.0; peakA=0 for A in range(2, 250, 2): bestZ=max(range(0, A+1), key=lambda z: binding(A, z)) per=binding(A,bestZ)/A if per>peakBE: peakBE=per; peakAx=A print(f"每核子结合能最大出现在 A≈{peakAx} (实验约56, 铁), 峰值≈{peakBE:.2f} MeV/nuc") for A in (16,56,120,235): Z=max(range(0,A+1), key=lambda z:binding(A,z)) print(f"A={A:3d} 结合能/核子={binding(A,Z)/A:.2f} MeV Z={Z}")
跑出来的曲线会清楚呈现一个高台:轻核结合能快升、铁附近登顶、重核缓缓回落。这个“中段最高的山脊”全局解释了许多宏观性质,也预示了重核裂变放能、轻核聚变放能的去处。
裂变必须是重核的偏好,用液滴模型能直接解释:库仑项随 Z² 增加,而表面项让碎片表面增多、代价变大。当库仑驱动力大于表面束缚,核就倾向于分裂成两块以降低总能量——这正是铀、钚等重核得以裂变的力学根源。模型还给出与裂变相关的“库仑临界”,说明越重越可能裂、越轻越安稳。虽然它看不见量子壳层,但在“整体稳定性”这一层已相当应手。

液滴模型擅长宏观:结合能大势、裂变能、同量异位素抛物线、滴线位置。它漏掉微观:不产生能级图、说不清幻数与自旋宇称——这些要交给下一节的壳模型。工程上它仍是估裂变能、算液滴形变的快速工具,很多核设备初算就用它。
一座反应堆一年多烧掉几十吨铀,其中能真正被榨出的是那点每核子 8 MeV 的游离差值。液滴模型把“裂变一块燃料能放出多少”算得足够好,工程师才能据此设计燃料装载量与热电转换。放射性年龄测定、同位素产量核算等也都复用其结合能结论。可以说,液滴模型是核工程做“能量预算”的头号计算器。
液滴模型是第2章质量亏损的模型化提炼,为第5章裂变/聚变提供力学判据,与第4章衰变能量同源。它和第6章其他模型共同构成“既要宏观又要微观”的多面观测。理解它,就握住了核物理里“先看大势、再看细节”的思维起点。
液滴模型看到的是“整个核的平均脸”,下一节换一副更细腻的镜子——壳模型,专门对准一个一个核子的能级座位。