6.2 壳模型
6.2 壳模型
本节摘要：核子不是无差异地泡在核里，而是按独立粒子轨道的能级逐格填充，构成壳层，幻数正是闭合壳层的户口登记。壳模型把核子之间的强作用近似成一个平均势场，再用自旋轨道耦合把能级掰弯，得到与实验吻合的幻数序列与基态自旋。本节用 Python 数值重现壳模型能级顺序与幻数演化，并以此解释为何氧-16、钙-40、铅-208 都特别稳。
学习目标
阅读完本节，你应当能够：
- 说明壳模型“平均势场 + 单粒子轨道”的基本假设
- 从能量填充顺序导出幻数 2、8、20、28、50、82、126
- 讲解自旋轨道耦合如何修正能级并造出大间隔
- 用壳层填充解释某核素的基态自旋宇称
- 举出闭合壳层核素在核结构与核数据里的突出地位
“核子也有座位表”
液滴模型把核当一锅粥，可实验无数次指出：当核子数是 2、8、20、28、50、82、126 时，特别稳。这强烈暗示核子里存在类似原子的“能级座位”。核物理的答案——壳模型——假定每个核子在一个平均势场中独立运动，逐个填入能级，形成闭合壳层时的幻数。强大作用被“平均化”成势场，剩下的才是个体的差别。
从简谐势到真实势与自旋轨道耦合
先看最朴素的谐振子势，能级填充 2、6、12、20…，前几个还能对，到后面就对不上幻数了。真正的魔法来自自旋轨道耦合：轨道角动量 l 与自旋 s 咬合，把 j=l+1/2 与 j=l−1/2 的能级分开，前者整体下沉。这样一来 1f7/2 沉到合适位置，闭合壳层恰好把幻数从 20 一步步推到 28、50、82、126。这就是迈耶与延森的获奖洞见。
用 Python 排出壳层填充顺序
# 壳模型能级填充顺序(计入自旋轨道耦合), 标出每一闭合壳层停止点(幻数) seq = [ ("1s1/2",2),("1p3/2",4),("1p1/2",2), ("1d5/2",6),("2s1/2",2),("1d3/2",4), ("1f7/2",8), # 沉下来 → 幻数28 ("2p3/2",4),("1f5/2",6),("2p1/2",2),("1g9/2",10), # 幻数50 ("2d5/2",6),("1g7/2",8),("3s1/2",2),("2d3/2",4),("1h11/2",12), # 幻数82 ("1h9/2",10),("3p3/2",4),("3p1/2",2),("2f7/2",8),("1i13/2",14), # 幻数126~... ] acc=0 print("累积填充(停顿点即幻数):") for name, occ in seq: acc += occ print(f" {name:8s} +{occ:2d} → 累计 {acc:3d}")
逐步累加下来的“台阶”正是幻数序列的来历。你会在计数时清楚地看到：1f7/2 这 8 个粒子正好让累计停在 28，而 1h11/2 等又把累计推到 50、82、126。跑一遍这段代码，幻数就不再是背诵口诀，而是“累计填充正好落在格子分界”的必然。
壳模型怎么报自旋宇称
核的基态自旋宇称常常由最外层空穴的那个没配对的核子决定。比如氧-17（8个质子+9中子）的基态自旋由第 9 个中子落进 1d5/2 而给出 5/2+；而闭合壳层的氧-16 自旋为 0+、特别紧。这套规则让谱学家不用查表就能预判一大批核素的量子标签，也解释了为何幻数核高稳、难激发。
图：壳模型的能级台阶与幻数
真实案例：为何氧-16 与铅-208 是“核界顶流”
氧-16（8p+8n）全闭合、基态 0+，宇宙丰度极高；铅-208（82p+126n）双幻数闭合双壳层，放射性极低、稳定异常，被当作放射实验里干净的“背景参考”。壳模型把这些稳定性翻译成“外填空、深填实、激发需花大钱”，也让核工程挑选屏蔽与结构材料时优先考虑这类不易被激发的幻数核。
壳模型的边界
壳模型对单粒子激发、幻数附近核态极为成功，可对同质异能态、强集体变形却力不从心——那是下一节集体模型的舞台。模型还抽象掉了成对关联与整体形变，需要人们用配位体模型与微观理论去补。总之壳模型告诉我们“谁是坐在座位上的粒子”，下一节则问“整个屋子会不会集体变形”。
与前后章节的关联
壳模型的幻数在第2章已有定性预告，在第7章谱学里又用它标定能级；它还和第8章核医学的稳定核素选择、核安全的结构材料选用直接挂钩。理解壳模型，等于理解了原子核“微观秩序感”的来源。
本节要点回顾
- 平均势场 + 单粒子轨道是壳模型的两块基石
- 自旋轨道耦合把能级劈裂下沉，凑出与实验一致的幻数
- 幻数 2、8、20、28、50、82、126 出自逐步累计填充的台阶
- 基态自旋宇称常由最外层未配对核子决定，幻数核多为 0+
- 氧-16、铅-208 双闭合壳层高稳，是核结构与工程选材的宠儿
一件件核有的“刚性”，有的软得能扭成椭圆——这不是壳模型失灵，而是整个核在集体变形，下一节讲集体模型。