跳转到内容
← 返回知识库
核物理约 19 分钟阅读

核壳层模型与液滴模型

核物理核壳层模型液滴模型幻数结合能核形状玛丽亚·迈耶延森核裂变核稳定性

关键词

核壳层模型; 液滴模型; 幻数; 结合能; 核形变; 玛丽亚·迈耶; 延森; 魏茨泽克公式; 泡利不相容; 自旋轨道耦合; 核裂变理论

原子核是一堆紧密堆叠的球吗?

标题:"原子核就是核子堆在一起的球"——核的量子结构远比这复杂

"原子核是质子和中子紧密堆积的球"——这个图像在大众科普中很普遍,却完全忽视了原子核内部丰富的量子结构。核子(质子和中子)在原子核中并非随机堆积,而是和原子中的电子一样,占据着离散的量子壳层(nuclear shells)。特定质子数或中子数(所谓"幻数":2, 8, 20, 28, 50, 82, 126)的原子核表现出异常的稳定性,类似于惰性气体原子的电子层填满时的化学稳定性——这不是巧合,而是量子力学(泡利不相容原理+自旋轨道耦合)在核尺度上的直接体现。

核子不能像弹珠那样挤成一堆,是因为它们服从泡利不相容原理。质子与中子各自占据分立的量子态,同一态至多容纳两个自旋相反的同类核子。于是核内出现壳层:能量由低到高填满,而不是连续地堆叠。

液滴图像之所以仍然有用,是因为核力程极短、近乎饱和。内部密度大体均匀,表面则像液体一样可以起伏。集体形变与裂变,正是从这种可变形的表面开始讲起的。

两种截然不同的图像——液滴模型(把核当成宏观液滴,强调集体行为)和壳层模型(把核当成量子多体系统,强调单粒子结构)——看似矛盾,实际上互补,共同构成了我们理解原子核的基础。核物理是量子力学在强相互作用束缚态中的最复杂应用之一,并非"用旧理论解决旧问题",而是至今仍有大量未解问题(例如:核力从QCD的第一原理推导、核物质状态方程在极端密度下的行为)。

两种模型回答的问题并不相同。问结合能随质量数如何平滑变化、问裂变势垒从何而来,液滴公式已经够用。问某个核素为何特别稳定、磁矩为何偏离经典预期,就必须回到单粒子能级。它们不是前后取代的理论,而是同一对象的两种分辨率。

1935年的液滴与1949年的壳层:两种图像如何各自解释核的行为

标题:两个互补的理论——魏茨泽克液滴(1935年)与迈耶-延森壳层模型(1949年)

液滴模型:1935年,卡尔·魏茨泽克(Carl Friedrich von Weizsäcker,1912—2007)提出将原子核类比为带正电荷的不可压缩液滴,写出了著名的半经验结合能公式(Bethe-Weizsäcker formula,又称魏茨泽克公式):

B(A,Z)=avA−asA2/3−acZ(Z−1)A1/3−aA(A−2Z)2A+δ(A,Z)B(A,Z) = a_v A - a_s A^{2/3} - a_c \frac{Z(Z-1)}{A^{1/3}} - a_A \frac{(A-2Z)^2}{A} + \delta(A,Z)

其中五项分别是:体积项(核力的短程饱和性,av≈15.8 MeVa_v \approx 15.8\ \text{MeV})、表面项(表面核子结合能少,as≈18.3 MeVa_s \approx 18.3\ \text{MeV})、库仑项(质子间斥力,ac≈0.714 MeVa_c \approx 0.714\ \text{MeV})、不对称项(质子中子数不等时能量增加,aA≈23.2 MeVa_A \approx 23.2\ \text{MeV})和配对项(偶偶核更稳定,δ>0\delta > 0)。液滴模型成功解释了核结合能的总体趋势和核裂变(波尔和惠勒,1939年,用液滴模型分析了铀核裂变过程)。

魏茨泽克把原子核写成带电液滴,并不是比喻先行、公式随后。体积项来自核力的短程饱和:每个核子只跟近邻相互作用,结合能大致与核子数成正比。

表面项扣除那些少了邻居的核子。库仑项计入质子之间的长程排斥。不对称项则惩罚质子、中子数目相差过大。配对项承认偶偶核更稳,那是后来壳层图像才会讲清楚的细节,在宏观公式里先被写成一个符号。

汉斯·贝特在稍后的综述里把同一套半经验公式系统化,因此文献常称贝特–魏茨泽克公式。玻尔与惠勒在一九三九年用这滴带电液体分析铀核的形变路径:表面能与库仑能此消彼长,越过势垒便裂成两块。液滴模型因此从结合能曲线走进了核裂变动力学。

壳层模型:但液滴模型无法解释幻数核($Z$ 或 $N = 2, 8, 20, 28, 50, 82, 126$)的异常稳定性、激发能隙大、磁矩特异等性质。1949年,玛丽亚·格佩特-迈耶(Maria Goeppert Mayer,1906—1972)和汉斯·延森(J. Hans D. Jensen,1907—1973)各自独立发现:如果在核势阱(有限势阱或谐振子势)之外加上强的自旋-轨道耦合(spin-orbit coupling,−l⋅s-l \cdot s 项,比原子中大约100倍),单粒子能级的简并性被部分解除,恰好在 $N$ 或 $Z = 2, 8, 20, 28, 50, 82, 126$ 处出现能隙,形成壳层闭合。迈耶和延森于1963年共同获诺贝尔物理学奖(迈耶是继居里夫人之后第二位获得诺贝尔物理学奖的女性)。

迈耶当时在芝加哥与阿贡实验室整理核素的自旋与磁矩表。延森则与奥托·哈克塞尔、汉斯·聚斯在海德堡从另一组数据得出同一结论。两边都发现:只靠谐振子或方势阱填壳,幻数会停在错误的位置。必须加上很强的自旋–轨道耦合,高角动量轨道才被劈开,闭合壳层才落到已经写在文中的那些幻数上。

一九六三年的诺贝尔奖承认的是这条共同的发现路径,而不是某一次单独的灵感。迈耶一九四八年关于闭合壳层的短讯,与一九四九年哈克塞尔、延森、聚斯的通信,以及迈耶次年对自旋–轨道组态的展开,构成同一条证据链。壳层模型因此是被核素表逼出来的,不是从势阱形状先验推出来的。

幻数2、8、20、28、50、82、126:自旋轨道耦合如何打开能隙

标题:幻数、壳层闭合与核形状

核壳层模型的单粒子势阱通常用谐振子势加修正项来近似: V(r)=12mω2r2−κℏω(2l2−⟨l2⟩N)−κμℏωl⋅sV(r) = \frac{1}{2}m\omega^2 r^2 - \kappa\hbar\omega(2\mathbf{l}^2 - \langle l^2 \rangle_{N}) - \kappa\mu\hbar\omega \mathbf{l}\cdot\mathbf{s}

自旋轨道耦合项使高轨道角动量(高 $l$)态大幅分裂:$l$ 较大的态的 $j = l+1/2$ 分量下移到下一个壳层,在正确的位置产生能隙,对应幻数。

自旋–轨道项把同一轨道角动量的两个总角动量劈开。角动量较大的那个分量被压低,有时会掉进下一个主壳,于是在特定核子数处突然出现宽能隙。能隙之上激发一个核子要付很高的代价,所以双幻核的第一激发态特别高。这不是液滴表面张力能直接给出的结果。

幻数的实验证据:

  • 双幻核($Z$ 和 $N$ 都是幻数)极端稳定:4^{4}He($Z=N=2$)、16^{16}O($Z=N=8$)、40^{40}Ca($Z=N=20$)、48^{48}Ca($Z=20, N=28$)、208^{208}Pb($Z=82, N=126$)是稳定性最强的核,其第一激发态能量特别高,核半径特别小。
  • 幻数核的自然丰度高于预期(e.g. 40^{40}Ca、56^{56}Fe、90^{90}Zr丰度相对较高)。
  • r过程元素合成在幻数处形成"峰"——$N=50$、$N=82$、$N=126$ 处的中子捕获速率减慢,产生元素丰度峰(对应Sr/Y/Zr、Ba-La-Ce、Pb/Bi-Po区域)。

幻数处的丰度峰,不只是核结构教科书里的摆设。快中子俘获经过这些壳层时截面下降,核素在那里堆积,再经β衰变回到稳定谷。中子俘获核合成因此把实验室里的壳间隙,翻译成了太阳系元素表上的峰。超新星与重元素起源与中子星合并,是这条翻译最剧烈的现场。

核形状:大多数原子核不是完美的球形。壳层部分填充时,未配对核子的各向异性相互作用可以使核变形——轴对称椭球形(长椭球"拉长"或扁椭球"扁平")。核形变用四极形变参数 β2\beta_2 描述(β2>0\beta_2 > 0:拉长,β2<0\beta_2 < 0:扁平)。变形核(如稀土区 Z≈64Z \approx 64—76,N≈88N \approx 88—112)的转动谱(rotational spectrum,EJ=ℏ22IJ(J+1)E_J = \frac{\hbar^2}{2\mathcal{I}}J(J+1),I\mathcal{I} 是转动惯量)是集体运动的典型表现,由玻尔(Aage Bohr,尼尔斯·玻尔之子)和莫特尔森(Ben Mottelson)系统研究,两人与雷恩沃特(James Rainwater)于1975年共同获诺贝尔物理学奖。

形变核的转动谱看起来像一个刚体在转,其实转动的是整个自洽势场。斯文·戈斯塔·尼尔逊在一九五五年把壳层模型搬到变形势里,单粒子能级随形变参数移动。某些核素在拉长或压扁时反而降低了能量。集体运动与单粒子运动在这里咬合:玻尔–莫特尔森的贡献,正是把液滴的形状自由度与壳层的单粒子谱写成同一套语言。

稳定岛在哪里?超重元素与"消失的幻数"

标题:"岛屿稳定性"与放射性核束——核壳层的极端区域

超重元素的"稳定岛":液滴模型预言 $Z > 100$ 的超重元素因库仑斥力应在极短时间内裂变,无法稳定存在。但核壳层模型预言:在下一个质子幻数(理论预言为 $Z = 114$ 或 $Z = 120$ 或 $Z = 126$,有争议)附近存在一个"稳定岛"(island of stability),其中超重元素的自发裂变寿命可能延长数个数量级。目前已合成到 $Z = 118$(𬊛,oganesson,2002年,JINR,杜布纳);$Z = 114$ 的寿命约0.5秒,尚未找到明显的稳定岛。更重元素的合成是当前核物理实验的前沿目标(RIKEN、GSI/FAIR、JINR)。

稳定岛是外推,不是已经上岸的陆地。液滴图像认定库仑斥力会随质子数迅速压倒表面张力,超重核应几乎立刻裂变。壳层修正却可能在下一组闭合壳附近把寿命拉长若干数量级。理论对质子幻数落在何处仍有分歧,文中已列出的几种候选都还在计算与实验之间摇摆。

已经合成的最重元素证实超重核可以短暂存在。但它们的寿命仍以秒或更短计,远不是可以装进试剂瓶的稳定核素。把稳定岛写成已经制成,是把预言地图误认成了实测海岸线。

"消失"的幻数:在远离稳定线的奇特核(富中子核或富质子核)中,"岸上"观测到的幻数可能不再成立。例如,实验发现 $N = 8$(He、Li同位素)和 $N = 20$(Na、Mg同位素,所谓"sd壳崩溃区")在富中子区域的壳层间隙缩小甚至消失,出现新的幻数(如 $N = 14, 16, 32, 34$)。这说明核壳层结构随着中子-质子比的变化而演化,是放射性核束实验的核心研究对象,涉及核力的三体效应和张量力成分。

所谓幻数消失,指的是壳间隙在极端中子质子比下被压窄,不再构成激发能谱上的明显门槛。这并不推翻岸上那些幻数,只说明单粒子势本身依赖于核子组成。三体力与张量力在富中子核里变得更刺眼。放射性核束装置要测的,正是这张随组成而变的能级图,而不是一张永恒的幻数表。

从临界质量计算到中子星合并:核结构理论的跨领域遗产

标题:从核裂变到天体物理s/r过程——核结构理论的跨领域影响

核壳层模型和液滴模型不是互斥的,而是在不同能量和尺度下各有适用范围。现代核物理采用的是二者的综合:集体模型(interacting boson model,IBM)、壳模型配对作用、以及近年兴起的手征有效场论+格点QCD方法。

现代计算不再强迫人在液滴与壳层之间单选。平均场给出变形与密度,组态相互作用再把价核子的关联补上。手征有效场论则试图从量子色动力学的对称性出发,写出核子之间的力,而不是从拟合参数倒推。第一原理到重核仍然艰难,这是未完成的推导,不是已经关闭的课题。

核裂变理论:波尔-惠勒(1939年)用液滴模型的表面张力与库仑斥力之比,计算了铀核自发裂变和诱发裂变的势垒——这直接奠定了曼哈顿计划中铀和钚的临界质量计算基础。

结合能的平滑部分由液滴公式承担,起伏部分交给壳修正。这正是结合能与质量亏损曲线在铁峰之后仍能被读懂的原因。临界质量估算吃的是液滴外推;幻数核对裂变势垒的修正,则提醒工程不要把半经验公式当成逐核真理。

元素合成的核数据:r过程(快中子俘获)在极端中子丰富环境(中子星合并,2017年LIGO/Virgo引力波+电磁对应体AT 2017gfo首次直接证实)中经过幻数壳层时的捕获截面,决定了金、铂、铀等重元素的宇宙丰度。精确的核壳层模型是预测r过程核素网络的基础。

天体物理把同一套结构数据接到极端环境。中子星外壳的核组成、合并事件里的r过程路径,都依赖于幻数在丰中子区是否还在原来的位置。致密星——白矮星与中子星内部可能出现的超流中子,也是核多体理论的外推,而不是已经采样完毕的样品。

核医学:幻数核(尤其是双幻核 208^{208}Pb)的高稳定性使其成为放射性药物(如Pb-212、Pb-203)的理想核素基础。核壳层对同质异能态(isomer)寿命的决定,也影响着医用放射性同位素(如Tc-99m,最常用的核医学诊断核素)的应用。

事实卡

  1. 卡1

    魏茨泽克半经验结合能公式(1935年)用五项(体积、表面、库仑、不对称、配对)拟合核结合能曲线,是核物理中应用最广泛的公式之一。

  2. 卡2

    核壳层模型的幻数(2, 8, 20, 28, 50, 82, 126)来自核势阱中强自旋-轨道耦合引起的能级分裂,由迈耶和延森于1949年独立发现,1963年共同获诺贝尔物理学奖。

  3. 卡3

    玛丽亚·格佩特-迈耶是继居里夫人之后第二位获得诺贝尔物理学奖的女性,以核壳层模型的发现闻名。

  4. 卡4

    远离稳定线的奇特核中,"传统"幻数(如N=20)可能消失,新幻数(如N=32,34)可能出现,是当前放射性核束实验(RIKEN、GSI)的核心研究主题。

跨域连接

  • 结合能与质量亏损:液滴公式用体积、表面、库仑、不对称、配对五项拟合整条结合能曲线,平均趋势漂亮吻合,逐核的偏差则全留给壳层效应解释。推论:配对项其实就是壳层效应在宏观公式里打的补丁——两个模型的接缝清晰可见。
  • 范式:液滴模型把核当连续介质,壳层模型把核当单粒子的量子体系,看似矛盾,实则分工——问裂变势垒用前者,问幻数用后者。推论:"哪个模型对"是个坏问题;模型的适用范围由问题的类型决定,这是科学模型观的教科书案例。
  • 元素周期表:幻数核的异常稳定与惰性气体的化学惰性同出一源——都是泡利原理下的壳层填满,一个填电子,一个填核子。推论:双幻核像惰性气体一样"不活泼"且丰度偏高,核素版图上也有着自己的稀有气体族。
  • 核工程:玻尔与惠勒用液滴模型的表面张力对库仑斥力之比算出裂变势垒,直接奠定了临界质量的估算方法。推论:半经验公式的工程价值恰恰在于敢外推——参数由已知核素定出,就能估算尚未测量的核。
  • 超重元素合成:液滴模型预言 Z 过百的核应瞬间裂变,壳层模型却预言 Z 约 114、N 约 184 附近存在"稳定岛"——两个模型的外推在超重区正面冲突。推论:合成更重的核素并测量其寿命,是对壳层模型外推能力的判决性检验——每多合成一种超重核素,就是给这场模型之争投出一票。

参考文献

  1. Weizsäcker, Carl Friedrich von. "Zur Theorie der Kernmassen." Zeitschrift für Physik, 96(7–8): 431–458, 1935.
  2. Mayer, Maria Goeppert. "On Closed Shells in Nuclei." Physical Review, 74(3): 235–239, 1948.
  3. Mayer, Maria Goeppert. "Nuclear Configurations in the Spin-Orbit Coupling Model." Physical Review, 78(1): 16–21, 1950.
  4. Bohr, Aage, and Ben R. Mottelson. Nuclear Structure, Vols. I & II. Benjamin, 1969 & 1975.
  5. Casten, Richard F. Nuclear Structure from a Simple Perspective (2nd ed.). Oxford University Press, 2000.
  6. 卢希庭. 《原子核物理》(修订版). 原子能出版社, 2001.
  7. Bohr, Niels, and John A. Wheeler. "The Mechanism of Nuclear Fission." Physical Review, 56(5): 426–450, 1939.
  8. Haxel, Otto, J. Hans D. Jensen, and Hans E. Suess. "On the \"Magic Numbers\" in Nuclear Structure." Physical Review, 75(11): 1766, 1949.
  9. Nilsson, Sven Gösta. "Binding States of Individual Nucleons in Strongly Deformed Nuclei." Kongelige Danske Videnskabernes Selskab, Matematisk-Fysiske Meddelelser, 29(16), 1955.
  10. Bethe, H. A., and R. F. Bacher. "Nuclear Physics A. Stationary States of Nuclei." Reviews of Modern Physics, 8(2): 82–229, 1936.