曹颖逾 郭建友

(安徽大学物理与材料科学学院, 合肥 230601)

1 引 言

众所周知, 原子核电荷半径是描述原子核性质最基本的物理量之一. 核力的性质、核物质密度以及壳结构等方面的信息可以通过研究核电荷半径得到. 核电荷半径测量的精确程度是验证原子核微观模型和理论的重要依据.

随着科技的发展, 实验仪器越来越精密, 技术和测量方法也越来越完善, 对物理量的测量也更加精确. 核电荷半径可通过高能电子在原子核上的散射[1]、原子光谱超精细结构的同位素移动[2]、 µ− 原子激发X 射线谱[3,4]、原子光谱[5,6]来获得. 在文献[7]中收集了909 个基态均方根核电荷半径, 这些数据中包含了远离b 稳定线的核素的核半径. 这表明核电荷半径的测量已经扩展到远离b 稳定线的区域.

在核电荷半径的理论计算和实验测量等方面,人们做出了大量工作. 早期, 人们通常使用简单的单参数公式 Rc=r0A1/3. 在文献[8,9]中, Nerlo 和Pomorski 提出包含同位旋的两参数公式与三参数公式, 指出同位旋与核电荷半径的关系. 文献[10]提出含Casten 因子项和d 项的四参数公式和五参数公式, 这两个公式给出了核电荷半径与壳效应和奇偶摆动的关系. 由于以 A1/3律的单参数公式缺乏物理量的解释, 曾谨言[11,12]提出了以 Z1/3律的单参数公式. 张双全等[13]提出含有同位旋的 Z1/3律公式, 该公式很好地描述了远离b 稳定线的核电荷半径. 另外壳层效应对核电荷半径的影响[14]、镜像核电荷半径的差异现象[15]和重核及其同位素半径的研究[16]等仍然是热点课题.

本文首先对上述公式进行拟合, 比较公式之间的优劣. 再考虑到原子核电四极矩与形变的关系,在三参数公式基础上提出加入电四极矩项的核电荷半径新公式. 最后对新公式进行改进, 分析新公式的优点和不足.

2 原子核电荷半径常用公式

核电荷半径通常由下面的等效半径表示:

其中r0是核电荷半径常数,A为质量数, 拟合后得核电荷半径常数r0= 1.2269 fm, 拟合图像如图1所示, 在轻核和重核区偏离较大, 中间区域符合较好. 由(2)式得到的核电荷半径理论值和实验值的均方根偏差为

通过实验数据和图1 的分析表明r0不是一个常数,对于轻核,r0较大, 到了重核r0逐渐减小. 这种系统性的偏离表示A1/3律有某些物理内容的缺失. 为了解决A1/3律在轻核与重核的偏离, 曾谨言[11,12]依据实验数据提出电荷半径的Z1/3律:

用最新数据拟合得到r0= 1.6394 fm, 均方根偏差s= 0.0939 fm, 拟合图像如图1 中右图所示. 对比A1/3律(2)式,σ减小了0.0285 fm, 这说明了Z1/3律比A1/3律描述得更好. 在文献[8,9]中提到两参数公式:

图1 Rc =r0A1/3 与 Rc =r0Z1/3的拟合曲线(左图是 Rc =r0A1/3的拟合曲线, 右图是 Rc =r0Z1/3 的拟合曲线)Fig. 1. The fitting curve of the Eqs. (2) and (4).(The left picture is the fitting curve of the Eq. (2) and the right picture is the fitting curve of the Eq. (4)).

式中有一个同位旋项 (N −Z)/A, 加入此项的目的是将远离 β 稳定线的核素的核电荷半径变化趋势更好地描述, 拟合(5)式得到拟合参数r0=1.2827 fm,a= 0.2700, 均方根偏差s= 0.0855 fm.为了更好地描述Z≽8 的偶偶核元素, (5)式加入1/A项, 得到三参数公式[9]:

拟合(6)式得到参数r0= 1.2331 fm,a= 0.1461,

b= 2.3301, 均方根偏差s= 0.0510 fm. (6)式是目前描述电荷半径变化较好的公式之一, 也是应用较为广泛的公式. 对比于(5)式和(6)式, 文献[13]提出如下公式:

其中N∗表示 β稳定线上核素的中子数,N/Z是与同位旋相关的物理量, 拟合(7)式得到r0=1.6312 fm,a= 0.0627, 均 方 根 半 径s=0.0618 fm.

为了更好地描述原子核电荷半径的变化规律,文献[10]提出以下两个公式:(8)式和(9)式中考虑到壳效应对核电荷半径的影响, 引入P项,P代表 Casten 因子[17,18], 实质意义是质子与中子的相互作用, 质子与中子的相互作用是导致原子核集团运动和形变的重要因素之一[19].P具体的表达式是NpNn/(Np+Nn) , Casten 因子广泛地应用于相位转移、形变、激发态性质、转动惯量等方面[20−23].δ项是描述奇偶摆动现象, 在同位素链中偶偶核的电荷半径比相邻的两个奇偶核的电荷半径的平均值大, 奇奇核的电荷半径比相邻的两个奇偶核的电荷半径的平均值小[10]. 可粗略表示为: 偶偶核的d= 1, 奇奇核的d= –1, 其他情况的d= 0. 对(8)式拟合, 得到拟合参数r0=1.2262 fm,a= 0.1473,b= 1.9876,c= 0.3993,均方根偏差s= 0.0273 fm. 对(9)式 拟合, 得到拟合参数r0= 1.2262 fm,a= 0.1472,b= 1.9952,c= 0.4030,d= 0.1093, 均方根偏差s= 0.0266 fm.(8)式和(9)式的均方根偏差都很小, 可以较好地描述核半径的变化.

3 形变与核电荷半径的关系

核的形变与核半径密切相关, 其中原子核的电四极矩是反映原子核电荷分布偏离球对称程度的物理量[24]. 电四极矩也是理解原子核内物质分布[25,26]、检验核模型[27]、核子与核子的相互作用的基本观测量之一[28,29]. 原子核电四极矩的研究引起了许多核物理的热点课题. 包括其中的奇异核粒子的研究[30]、原子核形状共存的研究[31]、以及壳演化的研究. 经过多年的研究, 人们对电四极矩有了大致的了解, 也开发了一系列的计算方法. 但是由于核结构较为复杂, 所以电四极矩的理论计算和实验测量都存在困难, 目前的实验测量和理论计算仍不尽完善. 电四极矩值随着核子数有着周期性变化, 计算出的理论值与实验值相比有较大的偏差[32].

电四极矩可由电势的多极展开得到, 原子核的内禀电四极矩有积分表示方法与张量表示方法, 具体计算过程见文献[33], 得到的内禀电四极矩表示为

其中e是质子的电荷元, 代入积分上下限, 得,Z为原子核的核电荷数.

假设与椭球体等体积球体的半径为R, 定义形变参量, 由等体积性得到半轴长则形变参量与内禀电四极矩的关系为

由(11) 式可计算出形变参量ε, 计算结果表明, 大多数原子核的形变参量都很小, 说明大多数核是偏离不大的非球形核. 在文献[33]中提出电四极矩与内禀电四极矩的关系为

当总自旋I= 0 或I=1/2 时, 电四极矩Q一定为0. 有了理论的支持, 再从文献[34]中取出378 个基态电四极矩, 去除中子数N≼8和质子数Z≼8的核素, 剩余368 个基态电四极矩. 这些电四极矩有正值与负值, 其中电四极矩的大小代表原子核偏离球形的程度, 正负号代表偏离的方向. 把电四极矩取绝对值, 表示成Q∗, 并在三参数公式(6) 式的基础上, 加上一个Q∗/A的项, 得到

拟合(13)式得到参数r0= 1.2221 fm,a= 0.1350,b= 2.4698,c= 0.8976, 均方根偏差s= 0.0397 fm.拟合图像如图2 所示, 相比图1, 图2 在轻核区和重核区符合更好.

与应用广泛的(6)式拟合的结果相比, 均方根偏差下降了, 说明用(13)式能够较好地描述核电荷半径的变化趋势.

图2 的拟合曲线图Fig. 2. The fitting curve of the Eq. (13).

由(12)式可知,Q还与总自旋I有关, 继续从文献[34]中取出对应的368 个核素的总自旋I, 反求内禀电四极矩Q0, 并取绝对值, 表示成, 替换(13)式中的Q∗, 得到

拟合(14)式, 得到参数r0= 1.2220 fm,a= 0.1410,b= 2.4200,c= 0.3643, 均方根偏差s= 0.0372 fm.拟合图像如图3 所示, 图3 在重核的一些区域表现得更平滑.

再取(9)式中能反映奇偶摆动现象的δ项, 将其加入(14)式中得到(15)式:

图3 的拟合曲线图Fig. 3. The fitting curve of the Eq. (14).

拟合(15)式,得到参数r0= 1.2223 fm,a=0.1421,b= 2.4577,c= 0.3660,d= 0.1705, 均方根偏差s= 0.0369 fm. 拟合图像如图4 所示.

为了看出(13)式—(15)式在具体的同位素链中拟合的好坏程度, 选取两个满壳附近的核Ba 和Fr, 与两个远离满壳附近的核Ho 和Lu. 把(13)式—(15)式计算的理论的核电荷半径值与实验值进行对比, 结果如图5 所示. 可以看出, 在实验曲线上的一些跳跃和转折, (13)式并不能很好地体现. 这些跳跃和转折现象与奇偶摆动、壳效应、

图4 的拟合曲线图Fig. 4. The fitting curve of the Eq. (15).

图5 Ba 和Fr, Ho 和Lu 四个同位素链核电荷半径的实验值与(13)式—(15)式计算的核电荷半径理论值的对比Fig. 5. The experimental values of the nuclear charge radii of Ba and Fr, Ho and Lu isotopic chains are compared with the theoretical values calculated by Eqs. (13)–(15).

形变相关. 把电四极矩项转变为内禀电四极矩项之后, 均方根偏差下降, 说明含有内禀电四极矩的(14)式与实验曲线符合得更好, 且(14)式在跳跃与转折的地方有了些许改善. 在加入了奇偶摆动项之后, 在跳跃与转折的地方符合得更好, 但是在一些没有出现跳跃与转折的地方, (15)式却表现出来跳跃与转折, 这表明仅仅加入这些项还是不足以描述核电荷半径复杂的变化.

为了从整体上看出理论值与实验值的差距, 计算了核电荷半径的实验值与(13)式—(15)式计算的理论值的差值 ∆R,做出 ∆R随质子数Z变化,结果如图6 所示. 可以看出, 整体趋势上(15)式得出的结果比(13)式和(14)式更加集中, 也印证了(15)式得出的均方根偏差最小. (13)式计算的 ∆R相对比较分散, 但在Z= 10 的轻核区, (13)式有少数核的 ∆R比(14)式和(15)式的 ∆R小. 在图6中还可以看到, 传统幻数Z= 28, 50, 82 所对应的∆R接近于零, 也印证了电四极矩与壳效应相关.对于轻核区, (13)式—(15)式的拟合结果不是很好, 计算的理论值和实验值差值较大, 这说明仅仅通过几个参数来精准地描述核电荷半径的变化趋势是很困难的. 以后还需要挖掘出蕴含在更深层的物理内容来对其进行补充. 最后, 将本文所提及的公式列在表1 中, 可以直观地看到各公式拟合的参数和精确程度.

图6 计算368 个核电荷半径的实验值分别与(13)式—(15)式计算的理论值的差值. (从上到下依次为核电荷半径的实验值与(13)式的差值图, 核电荷半径的实验值与(14)式的差值图, 核电荷半径的实验值与(15)式的差值图)Fig. 6. The difference between the experimental value of 368 nuclear charge radii and the theoretical value calculated by Eqs.(13)–(15) , respectively.

表1 各种核电荷半径公式Table 1. The mentioned equations for nuclear charge radius Rc .

4 结论与展望

本文对已有的核电荷半径实验数据和已经提出的核电荷半径公式进行验证. 并在三参数的基础上加入了电四极矩项, 发现计算结果的均方根偏差下降. 在已有的电四极矩数据的基础上, 又求出内禀电四极矩. 替换拟合公式中原本的电四极矩项,发现在质子数Z= 28, 50, 82 附近元素的实验值与公式计算的理论值相差非常小, 说明了内禀电四极矩项的加入, 很好地反映了壳效应的变化趋势,继续加入奇偶摆动项之后, 进行拟合并计算均方根偏差, 发现均方根又减小了. 相比于已经提出的公式, 新提出的公式更好地反映了原子核形变、壳效应和奇偶摆动等变化趋势, 在计算精度上也有所提高, 这将为今后的实验提供有益的参考.