-
考虑到
$ \cosh^{-2} $ 势的位垒穿透几率有严格的解析解,我们以此为标准来检验WKB近似和传输矩阵算法在计算势垒穿透几率时 的准确程度。$ \cosh^{-2} $ 势具体写为$ V(r) = \frac{V_0}{\cosh^2(\frac{r}{d})} $ 。其中$ V_0 $ 表示的是势垒高度,$ d $ 表示宽度参数。粒子穿透该势垒的几率[4]为其中
$ k = \sqrt{\frac{2m Q}{\hbar^2}} $ 。假设入射的粒子是质子,$ m = 938 $ MeV$ /c^2 $ 。这里取$ V_0 = 10 $ MeV,$ d = 10 $ fm。这样就可以利用传输矩阵、WKB方法以及解析公式(12)计算质子穿透势垒的几率。在本文研究中,$ V(r) $ 被划分为50 000个方势垒,因此能够保证计算的精度。相关结果显示在图2中。可以看到,当入射能量较高时,WKB方法、传输矩阵方法给出的结果都与精确的解析公式给出的结果一致。但当能量较低时,WKB方法给出的结果与精确值有较大差异,而传输矩阵方法给出的结果始终与精确解相符合,表明该方法比WKB近似方法更适合计算粒子穿透势垒的几率。WKB方法不能正确地给出入射粒子能量较低时的穿透几率,这与文献[4]中显示的结果相一致。 -
为探究WKB近似方法在计算
$ \alpha $ 粒子穿透该势垒几率时的误差,本文采用WKB近似方法以及传输矩阵方法计算了铀同位素链$ \alpha $ 粒子穿透势垒的几率,如图3所示。由图可知,$ \alpha $ 粒子穿透势垒的几率与衰变能有关。能量越大,穿透概率越大。对于铀同位素链$ \alpha $ 衰变时的穿透几率大小跨度了30个数量级,WKB近似方法和传输矩阵方法给出的势垒穿透几率都很接近。在图3中还显示了这两种方法给出的穿透几率之比,可以看到WKB近似方法给出的结果低于传输矩阵方法的结果,而且这种低估效应与衰变能有很好的依赖关系。总体上WKB近似方法给出的穿透几率比真实值低40%左右,这一结果与文献[2]中给出的结论也相似。因此在$ \alpha $ 衰变的研究中采用准确度更高的传输矩阵方法计算穿透几率十分有必要。 -
利用结团模型计算原子核
$ \alpha $ 衰变半衰期时的参数主要有势阱深度$ V_0 $ 、弥散宽度$ a $ 以及主量子数$ G $ 。文献[8]考察多个原子核$ \alpha $ 衰变后给出的一组参数设置为$ V_0 = 162.3 $ MeV,$ a = 0.40 $ fm,$ G = 20 $ (当$ N<126 $ )或者22($ N\geqslant 126 $ )。我们首先探究主量子数$ G $ 对半衰期的影响,相关结果显示在图4中。可以看到,$ G $ 越大计算得到的半衰期就越小;$ G $ 值变大2,则半衰期大概减小一个数量级。$ G = 20 $ (当$ N<126 $ )或者22($ N\geqslant 126 $ )这样的取值得到的计算值与实验值更加靠近,因此本文后面的计算都按照这种方式选取G。图5中显示的是半衰期受弥散宽度
$ a $ 的影响情况。可以看到随着$ a $ 增大半衰期减小,这与文献[9]中给出结果一致。整体上,弥散宽度$ a $ 从0.35变化到0.45 fm,计算的半衰期与实验值都能较好地符合。图6中显示的是势阱深度$ V_0 $ 对半衰期的影响,随着$ V_0 $ 的增大,半衰期增大。这是因为势阱越深,$ \alpha $ 粒子就越不容易穿透势垒,从而导致半衰期增大。考虑到势阱深度
$ V_0 $ 应该随原子核的同位旋变化[10],我们引入同位旋依赖的势阱深度参数$ V_0 = 152.5 (1+ 0.2\frac{N-Z}{N+Z}) $ MeV,重新调整了弥散宽度$ a = 0.39 $ fm。基于改进的参数计算的半衰期能与实验值更好地符合,其结果列在表1中。可以看到,大多数计算值与实验值之间的差异都小于3倍左右,表明计算的合理性。理论计算与实验数据有差异的原因可能是:(1) 原子核的形变效应还未被考虑在本次研究中。根据文献[12]的研究表明,考虑原子核的二极形变后计算得到的半衰期会减小;(2) 势阱以及原子核密度分布中的同位旋依赖性考虑得还不细致。文献[9]指出,引入弥散宽度和半径的同位旋依赖性后,计算得到的半衰期会降低,并且更加靠近实验值。文献[10]表明,在方势阱基础上考虑势阱深度对原子核质量数和质子数的依赖后,可以更好地解释实验数据。因此考虑更细致的同位旋依赖信息将对解释铀同位素链$ \alpha $ 衰变半衰期产生积极作用;(3)$ \alpha $ 粒子的形成因子的影响还没有在本次研究中考虑。通过对上述三方面的改进,可以预期理论计算将与实验值更加接近。计算结果与实验值差异最大的是$ {}^{217}{\rm{U}} $ 和$ {}^{219}{\rm{U}} $ ,理论计算值分别是实验值的20倍和0.04倍。$ {}^{217}{\rm{U}}$ 的计算值高于实验值,而${}^{219}{\rm{U}} $ 的计算值低于实验值的现象在其它理论模型的计算中也常被看到,如文献[13-15]采用机器学习方法以及密度依赖的结团模型(DDCM)同样有这一现象。文献[16]中采用推广的液滴模型(GLDM)也能看到这一现象。这两个原子核刚好跨越了$ N = 126 $ 的壳层,因此$ \alpha $ 粒子的形成因子将有较大变化;而且此时壳效应的贡献也较大,简单的相互作用势不能完全描述壳层的贡献。当理论模型中考虑的形成因子和量子效应不准确时,就会导致计算值与实验值有较大差异。此外,我们注意到AME2020给出的$ {}^{217}{\rm{U}} $ 半衰期评价值是$ 0.85\pm 0.71 $ ms,而文献[6]中给出的实验值为$ 19.3^{+13.3}_{-5.6} $ ms,文献[17]中给出的实验值为$ 15.6^{+21.3}_{-5.7} $ ms。这两个实验值与 本文计算的理论值17.22 ms非常接近。还可以明显地看到,对于偶偶核,计算的结果与实验符合得更好,计算的平均偏差$\Delta = \sqrt{\frac{1}{n}\sum_{i = 1}^n(\log_{10}T_{1/2}^{{\rm Cal},\,i}-\log_{10}T_{1/2}^{{\rm Expt},\,i})^2}$ = 0.30,表明计算值与实验值的平均偏差几乎在2倍以内;而奇A核通常相差较大,平均偏差$ \Delta $ = 0.73,这是因为目前的研究中$ \alpha $ 粒子的形成因子假定为1,这一假设对偶偶核的影响较小,但对奇A核的影响较大。考虑形成因子后,计算结果将与实验值更加靠近,相关研究正在进行中。N A ${Q}$/MeV $T^{\rm Expt}_{1/2}$ $T^{\rm Cal}_{1/2}$ $T^{\rm Cal}_{1/2}/T^{\rm Expt}_{1/2}$ 122 214 8.696 $0.52^{+0.95}_{-0.21}$ ms[7] $1.58$ ms 3.04 123 215 8.588 $1.4\pm 0.9$ ms $3.21$ ms 2.29 124 216 8.542 $6.9\pm 2.9$ ms $4.40$ ms 0.64 125 217 8.430 $0.85\pm 0.71$ ms $17.22$ ms 20.26 126 218 8.773 $0.35\pm 0.09$ ms $0.10$ ms 0.30 127 219 9.950 $60\pm 7$ $\mu {\rm{s}}$ $2.34$ $\mu {\rm{s}}$ 0.04 128 220 10.210 $60$ ns $41.4$ ns 0.69 129 221 9.889 $0.66\pm 0.14$ $\mu {\rm{s}}$ $0.21$ $\mu {\rm{s}}$ 0.31 130 222 9.481 $4.7\pm 0.7$ $\mu {\rm{s}}$ $1.78$ $\mu {\rm{s}}$ 0.38 131 223 9.158 $65\pm 12$ $\mu {\rm{s}}$ $19.50$ $\mu {\rm{s}}$ 0.30 132 224 8.633 $396\pm 17$ $\mu {\rm{s}}$ $273.07$ $\mu {\rm{s}}$ 0.69 133 225 8.009 $62\pm 4$ ms $33.89$ ms 0.55 134 226 7.702 $269\pm 6$ ms $189.32$ ms 0.70 135 227 7.235 $1.1\pm 0.1$ min $0.25$ min 0.23 136 228 6.799 $9.3\pm 0.2$ min $6.95$ min 0.75 137 229 6.473 $289\pm 2.5$ min $170.22$ min 0.59 138 230 5.993 $20.23\pm 0.02$ d $21.40$ d 1.06 139 231 5.580 $287.5\pm 6.8$ y $68.22$ y 0.24 140 232 5.356 $68.9\pm 0.4$ y $176.89$ y 2.57 141 233 4.908 $1.59\times10^5$ y $1.28\times10^5$ y 0.81 142 234 4.858 $2.45\times10^5$ y $2.87\times10^5$ y 1.17 144 236 4.571 $2.34\times10^7$ y $3.59\times10^7$ y 1.53 146 238 4.270 $4.463\times10^9$ y $9.68\times10^{9}$ y 2.17 实验值取自AME2020[11]以及最新的文献[7]。 -
最近我国科研人员测量到
$ {}^{214}{\rm{U}} $ 的$ \alpha $ 衰变的半衰期为$ 0.52^{+0.95}_{-0.21} $ ms,而表1中理论计算的结果为1.58 ms,此计算值还没有考虑$ \alpha $ 粒子形成因子的影响。文献[18]对$ \alpha $ 粒子形成因子与其衰变能的关联进行了细致研究,提出了一个计算形成因子的参数化公式,可以得到$ {}^{214} {\rm{U}}$ 中$ \alpha $ 粒子形成因子为0.024 7。在考虑这一形成因子的基础上,我们对势阱深度$ V_0 $ 和弥散宽度$ a $ 对半衰期的影响进行了分析,相关计算结果显示在图7中。可以看到半衰期随$ V_0 $ 和$ a $ 的变化规律与前文一致。在考虑形成因子的基础上,当$ V_0 $ 在125到137 MeV,$ a $ 在0.47到0.67 fm之间计算得到的半衰期都能与实验数据较好地符合。与实验中心值最接近的参数设置为$ V_0 = 129 $ MeV和$ a = 0.6 $ fm,在这组参数设置下得到的式(7)中势阱宽度参数$ R = 7.93 $ fm,这与原子核半径的经验公式$\sim1.3~A^{1/3}$ 得到的结果接近。文献[9]中给出的弥散宽度$ a $ 的经验公式为$ a = 0.54\big[1+0.787\; 95(N-Z)/(N+Z)\big] $ ,根据此公式计算得到$ a = 0.599\; 65 $ fm,这也表明我们所选参数的合理性。
Alpha Decay for Uranium Isotopes
doi: 10.11804/NuclPhysRev.39.2021083
- Received Date: 2021-08-10
- Rev Recd Date: 2021-10-19
- Available Online: 2022-06-29
- Publish Date: 2022-06-29
-
Key words:
- $\alpha$ decay /
- penetrability /
- half-live
Abstract: Both the transfer matrix approach and the Wentzel Kramers Brillouin(WKB) approximation are used to calculate the penetrability of
Citation: | Xiaoyun YE, Zichun WANG, Shuangshuang ZHANG, Yuting WANG, Yongjia WANG, Qingfeng LI. Alpha Decay for Uranium Isotopes[J]. Nuclear Physics Review, 2022, 39(2): 154-159. doi: 10.11804/NuclPhysRev.39.2021083 |