Learning objectives / 学习目标
- 从饱和、表面和库仑能建立液滴模型的尺度关系与用途。
- 逐项解释 SEMF，使用一套自洽系数估算质量并分析残差。
- 由 \(2j+1\) 简并度、强自旋—轨道分裂和 Pauli 填充说明幻数与基态 \(J^\pi\)。
- 用 \(J(J+1)\) 转动谱、近谐振动谱、\(E_{4^+}/E_{2^+}\) 与 \(B(E2)\) 识别集体性。
- 根据质量、低能能级、跃迁、电四极矩或反应任务选择模型并写明边界。
Prerequisites / 先修知识
- SEMF 趋势、稳定谷与分离能。
- 短程饱和、自旋—轨道与有效相互作用。
- 角动量耦合 与定态微扰。
Scope and conventions / 范围与约定
- 模型是对特定自由度和精度目标的压缩；任何模型预测都必须附上可观测量、参数来源和适用域。
- 液滴与 SEMF 描述平滑系统学；壳层模型描述离散单粒子/价核子结构；集体模型描述相关核子的形变转动和表面振动。
- 本课 SEMF 例题固定使用一套明示系数，绝不与 NNDC 的另一拟合集混用。
- 幻数 \(2,8,20,28,50,82,126\) 是稳定区经典闭壳层系统学；远离稳定线时壳隙可演化，不能写成永不改变的自然常数。
- \(E_{4^+}/E_{2^+}\) 与 \(B(E2)\) 是诊断指标而非单独的模型证明；过渡核、\(\gamma\) 软核和形状共存需更完整分析。
Notation and modeling contract / 符号与建模契约
| Symbol | Meaning | SI unit |
|---|---|---|
| \(B_{\rm LD}\) | liquid-drop/SEMF binding estimate | J |
| \(a_v,a_s,a_c,a_a\) | volume, surface, Coulomb, asymmetry coefficients | J |
| \(\delta\) | pairing contribution | J |
| \(n\ell j\) | dimensionless single-particle spectroscopic labels; corresponding angular-momentum operators have unit \(\mathrm{J\,s}\) | 1 |
| \(\mathcal I\) | collective moment of inertia | \(\mathrm{kg\,m^2}\) |
| \(\beta_2\) | quadrupole deformation parameter | 1 |
| \(B(E2)\) | reduced electric quadrupole transition probability | \(\mathrm{C^2m^4}\) |
1. Liquid-drop model / 液滴模型
Knowledge ID: P10-U02-L02-K01
液滴模型把核看成近似饱和、近似不可压缩、具有表面张力且带电的有限量子流体。体积能来自内部核子的充分近邻结合；表面核子缺少近邻，产生表面代价；质子库仑排斥随尺寸增大而增强。它自然给出 \(R\propto A^{1/3}\)、\(B\sim A\) 的主趋势，并适合讨论质量系统学、形变能与裂变的宏观部分。
“液滴”不是说核子是经典分子或有确定轨迹。Pauli 原理、不对称能、配对、壳修正和量子隧穿都超出简单经典液滴。核形变时，体积近似守恒而表面积与库仑能变化；裂变势垒就是表面代价与库仑获益、再加壳修正共同形成的能量景观。
图 1：内部与表面核子环境不同；圆点不表示经典固定位置。
图 2：裂变形变能可有势垒与次级谷；细节不能由纯液滴唯一给定。
2. Semi-empirical mass formula / 半经验质量公式
Knowledge ID: P10-U02-L02-K02
SEMF 把液滴趋势写成
\[B(A,Z)=a_vA-a_sA^{2/3}-a_c\frac{Z(Z-1)}{A^{1/3}} -a_a\frac{(A-2Z)^2}{A}+\delta(A,Z).\]
体积项正比于 \(A\)；表面项按面积 \(A^{2/3}\) 扣除；库仑项近似均匀带电球；不对称项惩罚过大的 \(N-Z\)，可由质子、中子分别填充费米态理解；配对项使偶偶核额外束缚、奇奇核相对抬高。不同版本使用 \(Z^2\) 或 \(Z(Z-1)\)，配对可取 \(A^{-1/2}\)、\(A^{-3/4}\) 等，系数必须与形式配套。
SEMF 是拟合到质量系统学的经验模型，不是从裸核力精确推出的定理。MIT 22.101 指出它在大部分中间质量区可给约百分之一量级结合能，但轻核、壳闭合、形变和滴线会出现结构性残差。NNDC NuDat 甚至提供 \([B-B_{\rm LDM}]/A\) 图，残差正是识别壳结构的入口。
图 3：SEMF 各项是能量账本，不应从不同参数集拼装。
图 4：系统性残差说明平滑液滴项遗漏了量子结构。
3. Shell model and magic numbers / 壳层模型与幻数
Knowledge ID: P10-U02-L02-K03
壳层模型把每个核子先置于其他核子产生的平均场中，再加入价核子间剩余相互作用。中心平均场给 \(n\ell\) 轨道，强自旋—轨道项把 \(j=\ell+1/2\) 与 \(j=\ell-1/2\) 显著分裂；每个 \(j\) 轨道可容纳 \(2j+1\) 个同类核子。按能量填充后出现大壳隙，经典质子或中子幻数为
\[2,\ 8,\ 20,\ 28,\ 50,\ 82,\ 126.\]
闭壳层通常表现更高第一激发能、更小集体性、分离能跃变与额外稳定。双幻核如 \(^{16}\mathrm O\)、\(^{40}\mathrm{Ca}\)、\(^{208}\mathrm{Pb}\) 是重要基准。但幻数不是由简谐振子单独得到：强自旋—轨道分裂对 28、50、82、126 至关重要。
极端单粒子模型用最后一个未配对核子预测奇 \(A\) 核 \(J^\pi\) 与磁矩，邻近闭壳层时尤其有效。远离闭壳层，价核子组态混合、配对和集体形变增大；远离稳定线，单粒子能量与壳隙可随中子—质子比变化。现代壳层计算必须声明模型空间、有效相互作用和有效算符。
图 5：自旋—轨道分裂重排能级并产生关键大壳隙。
图 6：幻数标记闭壳层；远离稳定线时壳隙可以演化。
4. Collective rotations and vibrations / 集体转动与振动
Knowledge ID: P10-U02-L02-K04
开壳层核的许多价核子可相干运动，形成形变平均场及转动带，或围绕平衡形状作表面振动。轴对称偶偶核的理想刚性转子基带满足
\[E_J=\frac{\hbar^2}{2\mathcal I}J(J+1), \qquad J=0,2,4,\ldots.\]
因此 \(E_{4^+}/E_{2^+}=10/3\approx3.33\)。理想谐 quadrupole vibrator 的两声子多重态约在第一 \(2^+\) 的两倍能量，常给 \(E_{4^+}/E_{2^+}\approx2\)。NNDC NuDat 用该比值作系统学指标，但明确它只是图表量；真实核常介于二者之间，且会有非谐性、\(\gamma\) 软性和带混合。
\(B(E2;0^+\to2^+)\) 测量电四极集体跃迁强度。常用尺度下
\[\beta_2=\frac{4\pi}{3ZR_0^2}\sqrt{\frac{B(E2)}{e^2}}, \qquad R_0=1.2A^{1/3}\ \mathrm{fm}.\]
NNDC B(E2) 数据库提供评估的偶偶核推荐值与不确定度；NuDat 还提醒由 \(B(E2)\) 推得的 \(\beta_2\) 在大不确定度或限值情形不可当作精确点值。实验室静态四极矩、内禀 \(Q_0\)、\(B(E2)\) 和 \(\beta_2\) 之间的关系依转动模型与约定。
图 7：理想转动带的能级间距随 \(J\) 增大。
图 8：理想谐振子多声子简并会被非谐性和相互作用劈裂。
图 9：\(\beta_2\) 参数化 quadrupole 形变；它不是直接拍摄到的几何边界。
5. Model domains and complementary descriptions / 各模型的适用范围与互补性
Knowledge ID: P10-U02-L02-K05
同一核可在不同问题中使用不同模型。总质量的平滑趋势先用 SEMF；闭壳层附近的低能 \(J^\pi\) 用壳层模型；强 \(E2\) 跃迁和 \(J(J+1)\) 带用集体模型；裂变形变能用液滴加壳修正；高精度轻核可用 ab initio；重核全局基态可用核密度泛函。选择依据是自由度、能量/长度尺度、目标可观测量和允许误差，而不是“哪个图像最像真实照片”。
壳层与集体性也非互斥：微观上，许多价核子组态的相干叠加可产生集体转动和振动；形变平均场又会重组单粒子轨道。液滴残差揭示壳层，壳修正又重塑裂变势垒。模型互补性要求用同一批数据交叉检验，而不是在每个数据点后随意更换参数。
图 10：以问题和可观测量选模型，而不是给核贴一个永久模型标签。
Worked Examples / 完整例题
Worked Example 1 — Surface fraction scaling / 表面比例标度
系统与参考系： 比较静止球形 \(A=16\) 与 \(A=128\) 核，质心系。
模型： 液滴均匀密度，表面层厚度不随 \(A\) 变；表面核子比例 \(f_s\propto A^{-1/3}\)。
坐标： 径向球坐标。
\[\frac{f_s(128)}{f_s(16)}=\left(\frac{128}{16}\right)^{-1/3} =8^{-1/3}=\frac12.\]
因此在该标度模型中，\(A=128\) 核的表面核子比例约为 \(A=16\) 的一半。
SI/单位检查： 比值无量纲；\(A\to\infty\) 时表面比例趋零，符合体积主导极限。真实轻核密度弥散和壳/团簇结构会破坏简单比例。
Worked Example 2 — SEMF estimate for \(^{56}\)Fe / SEMF 估算
系统与参考系： 静止 \(^{56}_{26}\mathrm{Fe}\) 核，质心系；球形液滴标量能量与坐标轴方向无关。
模型与系数： 本例固定 \(a_v=15.75\)、\(a_s=17.8\)、\(a_c=0.711\)、\(a_a=23.7\ \mathrm{MeV}\)，偶偶配对 \(\delta=+34A^{-3/4}\ \mathrm{MeV}\)；使用 \(Z(Z-1)\)。
已知： \(A=56\)、\(Z=26\)。
各项为
\[E_v=882.0,\quad E_s=260.5,\quad E_C=120.8,\quad E_a=6.77, \quad\delta=1.66\ \mathrm{MeV}.\]
\[B_{\rm SEMF}=882.0-260.5-120.8-6.77+1.66=495.6\ \mathrm{MeV}.\]
\[B/A=8.85\ \mathrm{MeV},\qquad B=7.94\times10^{-11}\ \mathrm J.\]
检查： 各扣除项符号为负、偶偶配对符号为正；若令表面、库仑、不对称与配对系数趋零，则回到 \(B=a_vA\)。结果接近实验约 \(492.3\ \mathrm{MeV}\)，偏差约 \(0.7\%\)；这正是经验模型精度而非“验证公式精确”。若改用 \(Z^2\) 或另一配对指数，必须重新拟合系数。
Worked Example 3 — Shell capacity and \(^{41}\)Ca / 壳容量与基态预测
系统与参考系： \(^{40}\mathrm{Ca}\) 惰性双闭壳层核芯加一个中子形成静止 \(^{41}\mathrm{Ca}\)，量子化轴 \(+z\)。
模型： 极端单粒子壳层模型，忽略核芯激发与组态混合。
轨道： 下一个中子轨道取 \(1f_{7/2}\)，\(\ell=3\)、\(j=7/2\)。
\[2j+1=2\left(\frac72\right)+1=8.\]
因此该轨道可从 \(N=20\) 填至 \(N=28\)。单个价中子给
\[J=j=\frac72,\qquad \pi=(-1)^\ell=-1,\]
即预测 \(J^\pi=7/2^-\). 最大投影的 SI 角动量为 \((7/2)\hbar=3.69\times10^{-34}\ \mathrm{J\,s}\)。
检查： 简并度无量纲，角动量单位是 \(\hbar\) 或 \(\mathrm{J\,s}\)；若 \(\ell\) 偶则宇称应翻为正。预测与评估基态相符，但磁矩和激发谱需要有效 \(g\) 因子与剩余相互作用。
Worked Example 4 — Rotational band and moment of inertia / 转动带
系统与参考系： 一个轴对称偶偶形变核的内禀轴转动，实验室系观察基带；取内禀对称轴为 \(z'\)、实验室量子化轴为 \(+z\)，且 \(E_{2^+}=80.0\ \mathrm{keV}\)。
模型： 理想刚性转子 \(E_J=\hbar^2J(J+1)/(2\mathcal I)\)，忽略离心伸长和带混合。
由 \(J=2\) 得
\[\mathcal I=\frac{3\hbar^2}{E_{2^+}}.\]
\(E_{2^+}=80.0\times10^3(1.602176634\times10^{-19})=1.2817\times10^{-14}\ \mathrm J\)，故
\[\mathcal I=\frac{3(1.054571817\times10^{-34})^2}{1.2817\times10^{-14}} =2.60\times10^{-54}\ \mathrm{kg\,m^2}.\]
\[E_{4^+}=\frac{10}{3}E_{2^+}=266.7\ \mathrm{keV}=4.272\times10^{-14}\ \mathrm J.\]
检查： \(E_4/E_2=3.33\) 符合理想转子；\(\mathcal I\to\infty\) 时能级间距趋零。真实核的配对会使 \(\mathcal I\) 偏离刚性均匀椭球值。
Worked Example 5 — Deformation from \(B(E2)\) / 由跃迁强度求形变
系统与参考系： 偶偶 \(A=150\)、\(Z=60\) 核的 \(0^+\to2^+\) 跃迁，实验室系测得 \(B(E2)=3.00\ e^2\mathrm{b^2}\)；内禀对称轴定义为 \(z'\)，实验室量子化轴为 \(+z\)。
模型： 轴对称转动模型与 \(R_0=1.20A^{1/3}\ \mathrm{fm}\)；忽略三轴性和高阶形变。
\[R_0=1.20(150)^{1/3}=6.38\ \mathrm{fm},\]
且 \(1\ \mathrm b=100\ \mathrm{fm^2}\)，所以
\[\sqrt{B(E2)/e^2}=\sqrt{3.00\times10^4}\ \mathrm{fm^2}=173.2\ \mathrm{fm^2}.\]
\[\beta_2=\frac{4\pi(173.2)}{3(60)(6.38)^2}=0.297.\]
SI 中 \(B(E2)=7.70\times10^{-94}\ \mathrm{C^2m^4}\)。
检查： \(\beta_2\) 无量纲且取正的形变幅度约 0.3，表示显著 quadrupole 集体性；\(B(E2)\to0\) 时 \(\beta_2\to0\)。若 \(B(E2)\) 只有上限，则 \(\beta_2\) 也只能报上限，不能报精确形状。
Lesson synthesis / 全课综合
液滴模型解释平滑质量、半径和形变能；SEMF 把这些趋势量化；壳层模型解释幻数、\(J^\pi\) 和离散能级；集体模型解释强 \(E2\)、转动带和振动多重态。它们在自由度和分辨率上不同，却通过残差、价核子相干与壳修正相互连接。模型的价值在可检验预测与已声明误差，不在把类比当作实体照片。
Formula-and-decision table / 公式与决策表
| Task | Relation | Use when | Boundary |
|---|---|---|---|
| smooth binding | SEMF five-term expression | global trend and rough masses | shells, deformation, drip lines need corrections |
| orbital capacity | \(2j+1\) | fill identical protons or neutrons | ordering depends on effective field |
| valence parity | \(\pi=(-1)^\ell\) | dominant single-particle configuration | configuration mixing can change simple assignment |
| rotational energies | \(E_J=\hbar^2J(J+1)/(2\mathcal I)\) | well-deformed axial band | transitional and centrifugal effects |
| rotor diagnostic | \(E_{4^+}/E_{2^+}=10/3\) | ideal even-even rotor | ratio alone is not proof |
| deformation estimate | \(\beta_2=4\pi[ B(E2)/e^2]^{1/2}/(3ZR_0^2)\) | stated collective convention | uncertainty and triaxiality matter |
Model boundaries / 模型边界
- 液滴模型描述平滑宏观能量，不给单个能级的精确 \(J^\pi\)。
- SEMF 的形式与系数成套；经验百分之一质量精度不延伸到壳残差和分离能细节。
- 经典幻数是稳定区系统学；极端中子/质子比可发生壳演化。
- 极端单粒子模型在闭壳层附近最好；开壳层需剩余相互作用与大模型空间。
- \(E_4/E_2\approx3.33\) 和大 \(B(E2)\) 支持转动性，但形状共存、三轴性和 \(\gamma\) 软性可产生复杂谱。
- \(\beta_2\) 是在指定集体模型下从电磁数据提取的参数，不是无模型几何测量。
- ab initio、壳层、DFT 与集体模型之间的参数和算符不能无条件转移。
Misconception clinic / 误区诊所
- 液滴模型不表示核子像经典液体分子沿确定轨道运动。
- SEMF 是半经验拟合，不是精确普适核 Hamiltonian。
- 幻数并非简单三维简谐振子全部给出；强自旋—轨道分裂至关重要。
- 闭壳层“特别稳定”不等于绝不发生任何衰变。
- 转动带不是整个核作为刚体的经典小球旋转；它是量子集体态序列。
- 壳层与集体模型不互相排斥，集体态可由许多壳层组态相干叠加形成。
Practice and retrieval / 练习与提取
- [K01] 液滴模型中为何有体积项和表面项？
答案 1
内部核子近邻数近似固定，给 \(A\) 标度；表面核子缺少近邻，数量按面积约 \(A^{2/3}\)，形成能量扣除。
- [K01] 液滴模型为何能讨论裂变却不擅长预测奇 \(A\) 基态自旋？
答案 2
裂变主要涉及整体形状的表面—库仑能；奇 \(A\) 自旋取决于离散价核子轨道和组态耦合。
- [K02] SEMF 中哪一项产生偶奇交错？
答案 3
配对项 \(\delta(A,Z)\)；偶偶为正、奇奇为负、奇 \(A\) 常取零，具体形式依参数集。
- [K02] 为什么不能把采用 \(Z^2\) 的库仑系数放进 \(Z(Z-1)\) 公式？
答案 4
系数已随选定函数形式整体拟合；改变变量形式却不重拟合会系统性双重改变库仑贡献。
- [K03] \(j=9/2\) 轨道可容纳多少个同类核子？
答案 5
\(2j+1=10\) 个；质子与中子各自按 Pauli 原理填充。
- [K03] 列出经典幻数并指出 28 以后为何需要强自旋—轨道分裂。
答案 6
\(2,8,20,28,50,82,126\)；中心振子简并度不能给出后续正确大壳隙，\(j=\ell\pm1/2\) 的强分裂重排能级。
- [K04] 理想转子若 \(E_{2^+}=120\ \mathrm{keV}\)，求 \(E_{4^+}\)。
答案 7
\(E_4=(10/3)E_2=400\ \mathrm{keV}\)。
- [K04] \(E_4/E_2\approx2\) 更接近哪一理想极限？
答案 8
谐 quadrupole vibrator 极限；真实核还应结合多声子谱与 \(B(E2)\) 检查。
- [K05] 若任务是预测闭壳层邻近奇核 \(J^\pi\)，首选哪个模型？
答案 9
壳层/价核子模型，并用评估能级检验；SEMF 只提供平滑质量趋势。
- [K05] 若任务是描述重形变核的强 \(E2\) 带，为什么仅用极端单粒子模型不足？
答案 10
强跃迁来自许多价核子的相干集体运动，需转动/集体模型或能产生同等集体性的微观多组态计算。
Exit ticket / 离堂检验
- 给出一个液滴模型成功但壳层模型更细的观测，以及一个相反的例子。
- 由一组 \(0^+,2^+,4^+\) 能量判断更接近转子或振子，并列出还需核对的第二项数据。
- 用一句话解释“模型互补”与“随意换模型拟合每个点”的区别。
References and source notes / 参考资料与来源说明
以下直接来源已于 2026-08-30 实际打开并核读：
- MIT OpenCourseWare 22.02, Lecture notes, Chapter 1. Introduction to Nuclear Physics. https://ocw.mit.edu/courses/22-02-introduction-to-applied-nuclear-physics-spring-2012/d0d046f78c917f107d925f11ac862ae4_MIT22_02S12_lec_ch1.pdf
- MIT OpenCourseWare 22.02, Lecture notes, Chapter 5. Nuclear Structure. https://ocw.mit.edu/courses/22-02-introduction-to-applied-nuclear-physics-spring-2012/5ac279a77468f9041a4d213442c3cb6a_MIT22_02S12_lec_ch5.pdf
- MIT OpenCourseWare 8.701, Lecture 9.5 Shell Model. https://ocw.mit.edu/courses/8-701-introduction-to-nuclear-and-particle-physics-fall-2020/resources/mit8_701f20_lec9-5/
- National Nuclear Data Center, Brookhaven National Laboratory, NuDat 3 User Guide. https://www.nndc.bnl.gov/nudat3/guide/
- National Nuclear Data Center, Brookhaven National Laboratory, Reduced Transition Probabilities, or B(E2; 0+ to 2+) Values. https://www.nndc.bnl.gov/be2/
- M. Verpelli and R. Capote Noy, Update of RIPL Nuclear Levels Segment, INDC(NDS)-0702 (2015), IAEA Nuclear Data Section. https://www-nds.iaea.org/publications/indc/indc-nds-0702/
- IAEA Atomic Mass Data Center, Atomic Mass Evaluation — AME 2020. https://www-nds.iaea.org/amdc/
Source-to-content map
| Source | Used for | Knowledge ID / Section | Boundary or convention |
|---|---|---|---|
| MIT 22.02 Chapter 1 | liquid drop and SEMF terms | K01, K02 | introductory parameterization |
| MIT 22.02 Chapter 5 | spin-orbit shells, magic numbers, valence spin and collective extension | K03–K05 | simplified shell ordering and extreme-particle limits |
| MIT 8.701 Lecture 9.5 | shell-model and magic-number overview | K03 | lecture overview, not a modern full-space calculation |
| NNDC NuDat guide | \(E_4/E_2\), \(\beta_2\), LDM residual and source provenance | K02, K04, K05 | diagnostic charts; uncertainties and limits matter |
| NNDC B(E2) | recommended quadrupole transition probabilities | K04; Example 5 | evaluated electromagnetic data, model needed for \(\beta_2\) |
| IAEA RIPL levels update | evaluated level-segment provenance | K03–K05 | data library, not a universal structure Hamiltonian |
| AME2020 | mass benchmark for SEMF residuals | K02 | evaluated masses supersede model estimates for precision |
复习与拓展
- 前课 Nuclear Interaction 说明有效两体/三体作用如何成为各模型参数的微观背景。
- 下一课 Decay Law and Statistics 将把结构与 \(Q\) 值转入随机衰变率描述。
- 后续 gamma 衰变、裂变和反应课程会分别复用壳层选择定则、集体 \(B(E2)\) 与液滴形变势垒。