-
根据Bhor等[20]关于复合核的形成和衰变过程相互独立假设与 Weisskopf[21]提出的核合成反应的三阶段模型,理论上通常将超重核合成动力学过程分为俘获阶段、熔合阶段(复合核形成阶段)以及复合核的退激发阶段。于是,目标超重核的生成截面可以表述为不同分波下三个因子乘积再求和,表示为[22]
$$ \sigma (E_{\mathrm{c.m.}}) = \sum\limits_{J}\sigma_{\mathrm{c}}(E_{\mathrm{c.m.}},\,J) P_{\mathrm{CN}}(E_{\mathrm{c.m.}},\,J) W_{\mathrm{sur}} (E_{\mathrm{c.m.}},\,J)\\, $$ (1) 其中:
$E_{\mathrm{c.m.}}$ 为质心系下的入射能;$\sigma_\mathrm{c}(E_{\mathrm{c.m.}},J)$ 为双核系统的俘获截面;$P_{\mathrm{CN}}(E_{\mathrm{c.m.}},J)$ 为熔合几率;$W_{\mathrm{sur}}(E_{\mathrm{c.m.}},\,J)$ 为存活几率。$ J $ 是入射角动量,受存活几率限制,一般最大取值到30[1]。 -
当两原子核从无穷远处相互靠近时,其间包含库仑力和核力相互作用,在以双核质心距离R为函数的有效势能曲线上会形成一个“口袋”和库仑位垒。俘获截面
$\sigma_\mathrm{c}(E_{\mathrm{c.m.}},\,J)$ [16]表示弹、靶核碰撞过程中克服库仑位垒后形成双核系统的几率:$$ \sigma_{\mathrm{c}}(E_{\mathrm{c.m.}},J) = \frac{\pi \hbar^2}{2\mu E_{\mathrm{c.m.}}}(2J+1) T(E_{\mathrm{c.m.}},J),\, $$ (2) 其中
$ \mu $ 是折合质量。穿透几率$T(E_{\mathrm{c.m.}},\,J)$ 受库仑位垒高度及口袋宽度的影响,表示为$$ T = \int\cfrac{f(B)}{1+\exp \left\{-{\cfrac{2\pi}{\hbar \omega(J)}} \left[ E_{\mathrm{c.m.}}- B-{\cfrac{\hbar ^{2}}{2\mu R_{\mathrm{B}}^2(J)}}J(J+1)\right] \right\}} \, \mathrm{d}B , $$ (3) 其中
$ \omega(J) $ 是库仑位垒位置$ R_{\mathrm{B}}(J) $ 对应的曲率刚度,$ f(B) $ 是非对称高斯分布函数[16]。考虑了动力学形变的核-核相互作用势为
$$ \begin{split} V(r,{\beta_1},{\beta_2},{\theta_1},{\theta_2}) = & V_{\rm{C}}(r,{\beta_1},{\beta_2},{\theta_1},{\theta_2})+V_{\rm{N}}(r,{\beta_1},{\beta_2},{\theta_1},{\theta_2})+\\& \dfrac{1}{2}C_1(\beta_1-\beta_1^{\,0})+\dfrac{1}{2}C_2(\beta_2-\beta_2^{\,0}), \\[-14pt] \end{split} $$ (4) 其中1和2代表弹核与靶核;
$ \beta_1 $ 与$ \beta_2 $ 是弹靶的动态四级形变参数;$\beta_1^{\,0}$ 与$\beta_2^{\,0}$ 是弹靶静态四级形变参数。$\theta_{i}$ (i = 1, 2)是第$ i $ 个核的径向矢量和形变核对称轴之间的夹角,硬度参数$C_{i}$ (i = 1, 2)由液滴模型[23]给出:$C_{{{{i}}\lambda}}$ = ($\lambda-1)\Big[$ $(\lambda$ +2)$ R_0^2 $ $ \sigma $ −$ \frac{3}{2\pi} $ $\frac{Z_i^2{\rm e}^2}{R_0(2\lambda+1)}\Big]$ ($\lambda\geqslant$ 2)。库仑相互作用
$V_{\rm{C}}(r,\,{\beta_1},\,{\beta_2},\,{\theta_1},\,{\theta_2})$ 采用Wong势[24]:$$ \begin{split} V_{\rm{C}} = & \dfrac{Z_{1}Z_{2}{ e}^{2}}{R}+\left(\dfrac{9}{20\pi}\right)^{1/2}\dfrac{Z_{1}Z_{2}{ e}^{2}}{R^{3}}\displaystyle\sum\limits_{i = 1}^2 \Re_{i}^{2}\beta_{2}^{\,i}P_{2}^{}(\cos\theta_{i})+ \\& \left(\dfrac{3}{7\pi}\right)^{1/2}\dfrac{Z_{1}Z_{2}{ e}^{2}}{R^{3}}\displaystyle\sum\limits_{i = 1}^2 \Re_{i}^{2}\big[\beta_{2}^{\,i}P_{2}^{}(\cos\theta_{i})\big]^{2} \text{。} \end{split} $$ (5) 对于核势采用Woods-Saxon势[25]:
$$ V_{\rm{N}} = - \frac{V_0}{1+\exp\left\{r-\displaystyle\sum\nolimits_{i = 1}^2 R_i\left[1+(5/4\pi)^{1/2}\beta_iP_2(\cos\theta_i) \right]/a \right\}}\text{。} $$ (6) -
弹、靶核在克服库仑位垒相接触后,两核之间发生核子转移。一个核的核子全部转移给另一个核时形成复合核,这个过程可以用主方程描述[5]。
设反应中始终保持两体过程,碎片1、2的中子质子数分别是(
$ N_1 $ ,$ Z_1 $ )、($ N_2 $ ,$ Z_2 $ )。设在$ t $ 时刻,碎片1在激发能$ E_1 $ 处的分布几率为$ P $ ($ N_1 $ ,$ Z_1 $ ,$ E_1 $ ,$ t $ ),各组态分布几率为$$ \begin{split} &\frac{\mathrm{d}P(N_{1},\,Z_{1},\,E_{1},\,t)}{\mathrm{d}t} = \sum\limits_{N_{1}^\prime}W_{{N_{1},\,N_{1}^\prime }}\big[d_{{N_{1}}}P(N_{1}^\prime,Z_{1}^{},E_{1}^\prime,\,t)-\\& d_{{N_{1}^\prime}}P(N_{1},Z_{1},\,E_{1},\,t)\big]+\sum\limits_{Z_{1}^\prime}W_{{Z_{1},Z_{1}^\prime }}\big[d_{{Z_{1}}}P(N_{1}^{},\,Z_{1}^\prime,E_{1}^\prime,t)-\\& d_{{Z_{1}^\prime}}P(N_{1},\,Z_{1},\,E_{1},\,t)\big] -\big[\varLambda{^{{\rm qf}}_{N_1,\,Z_1,\,E_1,\,t}}(\varTheta)+\varLambda{^{{\rm fis}}_{N_1,\,Z_1,\,E_1,\,t}}(\varTheta)\big]\times\\& P(N_1,\,Z_1,\,E_1,\,t), \end{split} $$ (7) 式(7)中,
$W_{{N_{1},\,Z_{1}^{};\,N_{1}^{\prime},\,Z_{1}^{}}}$ 表示碎片1从态 ($ N_1^\prime $ ,$Z_1^{}$ ,$ E_1^{\prime} $ )向态 ($ N_1 $ ,$ Z_1 $ ,$ E_1 $ )跃迁的几率,并对碎片1所能取的所有状态进行求和。$d_{\rm{{N_1},\,{Z_1}}}$ 表示碎片在($ N_1 $ ,$ Z_1 $ )宏观状态时的微观维度,由单粒子哈密顿量通过微观计算[26]得到。准裂变速率$\varLambda {^{\rm qf}_{A_{1},\,E_{1},\,t}}(\varTheta)$ 与裂变速率$\varLambda {^{\rm fis}_{A_{1},\,E_{1},\,t}}(\varTheta)$ 由一维的Kramers公式[27]给出:$$ \begin{split} \varLambda{^{{{\rm qf}}}_{{{A}}_{1},\,{{E}}_{1},\,{{t}}}}(\varTheta) =& \dfrac{\omega}{2\pi \omega^{B^{\rm{qf}}}}\left( \sqrt{\left(\dfrac{\varGamma}{2\hbar}\right)^2+(\omega^{B^{\rm{qf}}})^2}-\dfrac{\varGamma}{2\hbar}\right)\times\\& \exp\left(-\dfrac{B^{\rm{qf}}(A_1)}{\varTheta(A_{1},E_{1},t)}\right), \end{split} $$ (8) 其中
$ B^{\rm{qf}} $ 是双核组态的准裂变位垒,核温度为$\varTheta(A_{1},\,E_{1},\,t)$ =$ \sqrt{\epsilon^*/a} $ ,$ \epsilon ^* $ 为局域激发能,能级密度参数$ a $ =$\frac{A}{12}~{\rm{MeV}}^{-1}$ 。$ \omega^{B^{\rm{qf}}} $ 和$ \omega $ 分别为在准裂变位垒顶部和底部的曲率刚度相对应的谐振子频率。宽度$\varGamma$ [28]常取2.8 MeV。$$\begin{split} \varLambda{^{\rm fis}_{A_{1},E_{1},t}}(\varTheta) =& \dfrac{\omega_0}{2\pi \omega^{B^{\rm{f}}}}\left( \sqrt{\left(\dfrac{\varGamma_0}{2\hbar}\right)^2+(\omega^{\rm{B^{f}}})^2}-\dfrac{\varGamma_0}{2\hbar}\right)\times\\ &\exp\left(-\dfrac{B^{\rm{f}}(A_1)}{\varTheta(A_{1},E_{1},t)}\right), \end{split} $$ (9) 相应地,
$ {B^{\rm{f}}} $ 表示裂变位垒高度,$\omega^{B^{\rm{f}}}$ 和$ \omega_0 $ 分别为在裂变位垒和基态处对应的谐振子频率,宽度$\varGamma_0$ [28]= 2 MeV。只考虑能量跃迁和核子转移的情况下,质子跃迁几率可以写成如下形式:
$$ \begin{split} W_{{Z_{1}^{},\,N_{1}^{};\,Z_{1}^\prime,\,N_{1}^{}}} =& \dfrac{\tau_{\rm{mem}}(Z_{1}^{},N_1^{},E_1^{},Z_{1}^\prime,E_{1}^\prime;t)}{\hbar^{2}d_{{Z_1,\,N_1}}{d_{Z_{1}^\prime,\,N_1^{}}}}\times\\ & \sum\limits_{ii^\prime}\left |\langle Z_{1}^\prime,N_1^{},E_{1}^\prime,i^\prime | V |Z_1^{},N_1^{},E_1^{},i \rangle \right |^2,\, \end{split} $$ (10) 其中
$ i $ 代表所有的剩余量子数。在重离子反应过程中,相对运动动能耗散转移为两核子的内部激发能,使得两核在费米面两侧各自产生价空间
$ \Delta_{\epsilon K} $ ,其内的核子会被激发或发生转移,$$ \Delta_{\rm{\epsilon K}} = \sqrt{\frac{4\epsilon_{\rm{K}^{*}}}{g_{\rm{K}}}},\, $$ (11) 式中
$ \epsilon_{{\rm{K}}^*} $ =$ \epsilon ^*\frac {A_{\rm{K}}}{A} $ ,$ g{\rm_K} $ =$ \frac{A{\rm_K}}{12} $ ,$ A_{\rm{K}} $ =$ N_{\rm{K}} $ +$ Z_{\rm{K}} $ 。在价空间,价态的个数$ N_{\rm{K}} $ =$ \Delta_{\rm{\epsilon K}} $ $ g_{\rm{K}} $ ,受激的核子数$ m_{\rm{K}} $ =$ N_{\rm{K}} $ /2。微观维度数目:$$ d({m_1},{m_2}) = \left( {\begin{array}{*{20}{c}} {{N_1}}\\ {{m_1}} \end{array}} \right)\left( {\begin{array}{*{20}{c}} {{N_2}}\\ {{m_2}} \end{array}} \right) \text{。} $$ $ \epsilon^* $ 是双核系统$ (N_1,Z_1) $ 的局域激发能,$$ \epsilon^* = E_{\rm{diss}}-\big[U(Z_1,N_1)-U(Z_{\rm{P}},N_{\rm{P}})\big]-\frac{(J-M)^2}{2{ \zeta_{{\rm{rel}}}}}-\frac{M^2}{2{ \zeta_{{\rm{int}}}}},\, $$ (12) 其中
$ E_{\rm{diss}}(t) $ 表示耗散到双核系统的能量,随时间的变化而变化。$ J $ 、$ \zeta_{\rm{rel}} $ 分别为入射道的初始相对运动角动量和转动惯量,$ M $ 、$ \zeta_{\rm{int}} $ 分别为耗散成双核系统的内禀角动量和转动惯量。$ U(Z_1,N_1 ) $ 、$ U(Z_{\rm{P}},N_{\rm{P}} ) $ 分别为双核系统中的某一双核组态和入射点对应的势能,因此,双核系统$ (N_1,Z_1) $ 区域激发能和驱动势直接相关,跃迁几率最终要落到驱动势求解上。如图1所示,双核系统相互作用势随质心距离R存在变化关系。图 1
$^{48}{\rm{Ca}}$ +$^{238}{\rm{U}}$ 相互作用势随两核中心距离$R$ 的变化曲线[29]驱动势表述为
$$ \begin{split} U(N_1,Z_1;&N_2,Z_2,R) = B(N_1,Z_1)+B(N_2,Z_2)-\\ & [B(N,Z)+V_{\rm{rot}}^\prime(J)]+V(R,J), \end{split} $$ (13) $$ V(R,J) = V_{\rm{C}}(R)+V_{\rm{N}}(R)+V_{\rm{rot}}(R,J),\, $$ (14) 其中:
$ N $ =$ N_1 $ +$ N_2 $ ,$Z=Z_1+Z_2$ ;$B(N_1,\,Z_1)$ 、$B(N_2,\,Z_2)$ 、$B(N,\,Z)$ 分别是两碎片与复合核的结合能,包括了壳修正与对能修正;$ V_{\rm{C}}(R) $ 、$ V_{\rm{N}}(R) $ 分别是库仑势(式5)与核势;$V_{\rm{rot}}^{}$ 、$ V_{\rm{rot}}^\prime $ 分别是双核系统与复合核的离心势。由于人们对核力的认识还不够清楚以及核多体处理的困难,实际计算核-核相互作用势是按具体情况作一定的近似处理。通常以考虑介质效应的核子-核子有效相互作用并对核子密度分布进行折叠计算得到,而这种相互作用又有零程和有限程的近似处理,故相互作用势的计算具有模型依赖性,且因拟合的实验对象不同,所得到的核子-核子的有效相互作用或密度分布的参数(组)亦有不同的适用范围。目前还没有准确的相互作用势(特别是核相互作用势)的计算公式或方法。尽管不同形式的相互作用势计算或许会有所差别,但这种差别更可能的是系统性的平移,而俘获截面和熔合几率的计算结果也都有系统性的不确定,但它们的趋势没有改变,故最后所得结论也不会有太大影响。
对于核势
$ V_{\rm{N}}(R) $ ,我们采用Skyrme作用力基础上的双折叠势[30]:$$ \begin{split} & U_{\rm{N}}(R) = C_0\left\{\dfrac{F_{\rm{in}}-F_{\rm{ex}}}{\rho_{00}}\bigg[\displaystyle\int \rho_1^2(r)\rho_2(r-R) \, \mathrm{d}r + \right.\\& \left. \displaystyle\int \rho_1^{}(r)\rho_2^2(r-R) \, \mathrm{d}r\bigg]+F_{\rm{ex}}^{}\displaystyle\int \rho_1(r)^{}\rho_2(r-R)^{} \, \mathrm{d}r \right\}, \, \end{split} $$ (15) 其中:
$F_{\rm{in (ex)}}^{} = f_{\rm{in (ex)}}^{} + f_{\rm{in (ex)}}^\prime \boldsymbol\cdot \frac{N_1-Z_1}{A_1} \boldsymbol\cdot \frac{N_2-Z_2}{A_2} $ 。参数选择为:$ C_0 $ = 300 MeV·fm3,$f_{\rm{in}}^{}$ = 0.09,$f_{\rm{ex}}^{}$ = −2.59,$ f_{\rm{in}}^\prime $ = 0.42,$ f_{\rm{ex}}^\prime $ = 0.54。通过分步差分的方法数值解主方程,可以得到
$ t $ 时刻,激发能为$ E^* $ 的碎片1的分布几率$ P(Z_1,N_1,E_1,t) $ 。如图2所示,在入射道左侧,随着核子质量数不断减小,驱动势在不断增大,在 B.G.点达到最大,即为双核系统势能曲面的鞍点位置。当系统演化到B.G.点左侧时,双核系统发生熔合。于是熔合几率可表示为对B.G.点左边的分布几率求和[31]:$$ {P_{{\rm{CN}}}}({E_{{\rm{c}}.{\rm{m.}}}},J) = \mathop \sum \limits_{{Z_1} = 1}^{{Z_{{\rm{BG}}}}} \mathop \sum \limits_{{N_1} = 1}^{{N_{{\rm{BG}}}}} P({Z_1},{N_1},{E_1},{\tau _{{\rm{int}}}}) \text{。} $$ (16) 双核系统必须克服内部熔合位垒
$ B_{\rm{fus}} $ 才能发生熔合,这部分能量由局域激发能提供。相应地,对复合核激发能$ E^* $ 、入射能$E_{\rm{c.m.}}$ [29]有$$ E^* = U(\eta_{\rm{i}})+B_{\rm{fus}},\, $$ (17) $$ E_{\rm{c.m.}} = Q+E^*,\, $$ (18) 其中
$ U(\eta_{\rm{i}}) $ 为入射点的驱动势,$ Q $ 为合成复合核的反应能。 -
重离子熔合形成的复合核需要通过发射中子、带电粒子等方式退激发。对于重核,只有当激发能低于中子分离能时,发射带电粒子作用才比较明显。因此,目前考虑复合核退激发的存活几率时只考虑了中子蒸发与裂变的竞争[26]。
根据复合核的统计蒸发模型,自旋为
$ J $ 、激发能为$ E^* $ 的超重复合核蒸发$ x $ 个中子的存活几率表述为$$ W_{\mathrm{sur}}(E^{*},x,J) = P(E^{*})\prod\limits_{i = 1}^x \left[\frac{\varGamma_{\mathrm{n}}^{}(E{_{i}^{*}},J)}{\varGamma_{\mathrm{n}}^{}(E{_{i}^{*}},J)+\varGamma_{\mathrm{f}}^{}(E{_{i}^{*}},J)}\right]_{i},\, $$ (19) 其中
$ P(E^{*}) $ 是蒸发$ x $ 个中子退激的实现几率,由Jackson蒸发公式[32]给出:$$ P(E^*,x) = I(\Delta{_x},2x-3)-I(\Delta{_{x+1}},2x-1) \text{。} $$ (20) $\varGamma_{\rm{n}}(E^{*}_{i},J)$ 为复合核能级的第$ i$ 个中子的蒸发宽度,$\varGamma_{\rm{f}}(E^{*}_{i},J)$ 为蒸发第$ i $ 个中子的裂变宽度,它们都是$ E{_ i^*} $ 的函数,由Bohr - Wheeler公式[20]可知:$$ \begin{split} \varGamma_{\rm{n}} = \dfrac{2s+1}{2\pi\rho(E_i^*,J)}\dfrac{2m_{\rm{n}}R^2}{\hbar^2} \displaystyle\int_{0}^{E{_i^*}-B{_i^{\,\rm{n}}}-\frac{1}{a{_i}}-\delta} \epsilon \rho(E{_i^*}-B{_i^{\rm{n}}}-\epsilon,J)\, \mathrm{d}\epsilon , \end{split}$$ (21) $$ \begin{array}{l} \varGamma_{\rm{f}} = \dfrac{1}{2\pi\rho_{\rm{f}}(E{_i^*},J)} \displaystyle \int_{0}^{E{_i^*}-B{_i^{\,\rm{f}}}-\frac{1}{a_{{{\rm f}i}}}-\delta} \dfrac{\rho_{\rm{f}}(E{_i^*}-B{_i^{\,\rm{f}}}-\epsilon,J)}{1+\exp[2\pi(\epsilon+B{_i^{\rm{f}}}-E{_i^*})/\hbar\omega]}\, \mathrm{d}\epsilon , \end{array} $$ (22) 式 (21) 中,
$ m_{\rm{n}} $ 、$ R $ 、$ s $ 分别是蒸发第$ i $ 个中子前的中子质量、复合核半径和自旋量子数;$ \hbar\omega $ = 2.2 MeV。$B{_i^{\rm{n}}}$ 是第$ i $ 个中子的分离能;$B{_i^{\rm{f}}}$ 是蒸发第i个中子前的裂变位垒,能级密度$ \rho $ 由Fermi-gas模型[33]得到。 -
我们计算了115-118号元素各反应道的蒸发剩余截面,如图3所示,各反应道除
$ ^{48} {\rm{Ca}}$ +$ ^{249} {\rm{Bk}}$ 外均在4n蒸发道处有最大值,其中图例用*标注的为相应的实验值[34]。显然,各反应道在4n蒸发道处误差范围更小,在2n、3n蒸发道处误差大致在允许范围内,利用双核模型计算的结果与实验值基本符合。 -
研究超重核生成截面的同位素依赖性是为了寻找反应体系的最佳弹靶组合。一般在靶(弹)核确定的情况下,与合成超重核的弹(靶)核同位素组成反应体系,通过对各种体系的生成截面进行研究,可以为合成超重核提供重要信息[46]。
探究影响重离子核反应生成截面大小的反应机制是筛选出最佳弹靶组合的重要理论基础之一。对此,我们根据双核模型的数值计算结果,将重点从熔合阶段、存活阶段讨论影响生成截面同位素依赖性的主要因素,以此说明超重核反应理论机制。
-
对于合成Z = 119,120超重核,除
$ ^{50} {\rm{Ti}}$ 相关的弹靶组合外,理论工作者[6, 47]对$ ^{64} {\rm{Ni}}$ 、$ ^{54} {\rm{Cr}}$ 等弹核也展开讨论。在本文中,我们选取三条弹核同位素链:$ ^{52-59} {\rm{Cr}}$ +$ ^{243} {\rm{Am}}$ $ \rightarrow $ 119、$ ^{54-62} {\rm{Mn}}$ +$ ^{243} {\rm{Am}}$ $ \rightarrow $ 120、$ ^{56-72} {\rm{Ni}}$ +$ ^{238} {\rm{U}}$ $ \rightarrow $ 120,通过其各反应道研究超重核生成截面的同位素依赖性。图9给出了各反应道最大生成截面与复合核质量数
$ A $ 的关系。需要说明的是:考虑奇偶效应,$ 3{\rm n} $ 蒸发道分别用黑、红两种颜色表示,$ 4{\rm n} $ 蒸发道用蓝色单独表示。$ ^{56-72} {\rm{Ni}}$ +$ ^{238} {\rm{U}}$ 各反应道的情况相对复杂,其中对四个弹核$^{59,\,60,\,66,\,67} {\rm{Ni}}$ +$ ^{238} {\rm{U}}$ 过小的计算结果进行舍弃。由图9 (a) (b) (c)大致可以看出,弹核包含中子数越多,蒸发剩余截面越小。针对蒸发剩余截面同弹核同位素链的依赖性,我们首先讨论俘获阶段。一般重离子核反应的实验耗时长,在核反应实验中弹(靶)核素大多采用稳定元素。因此,对
$ ^{52-59} {\rm{Cr}}$ +$ ^{243} {\rm{Am}}$ 反应系统,重点关注$ ^{52} {\rm{Cr}}$ 、$ ^{53} {\rm{Cr}}$ 、$ ^{54} {\rm{Cr}}$ 三个半衰期长的弹核。将
$ ^{52-59} {\rm{Cr}}$ +$ ^{243} {\rm{Am}}$ 各反应道的俘获截面按其对应的弹核质量数$ A $ 进行排序,得到如图9 (d)所呈现的趋势,$ ^{52-59} {\rm{Cr}}$ +$ ^{243} {\rm{Am}}$ 各反应道最大值对应的俘获截面均在一个数量级上,侧面反映出对$ ^{52-59} {\rm{Cr}}$ +$ ^{243} {\rm{Am}}$ 系统,俘获截面并不是影响其蒸发剩余截面的重要因素。针对
$ ^{54-62} {\rm{Mn}}$ +$ ^{243} {\rm{Am}}$ 反应系统,我们重点关注弹核$ ^{54} {\rm{Mn}}$ 、$ ^{55} {\rm{Mn}}$ 。从图9 (e)也可以得出与$ ^{52-59} {\rm{Cr}}$ +$ ^{243} {\rm{Am}}$ 反应系统一致的结论。相对而言,
$ ^{56-72} {\rm{Ni}}$ +$ ^{238} {\rm{U}}$ 反应系统的情况较为复杂,我们重点考虑$ ^{58} {\rm{Ni}}$ 、$ ^{61} {\rm{Ni}}$ 、$ ^{62} {\rm{Ni}}$ 、$ ^{63} {\rm{Ni}}$ 、$ ^{64} {\rm{Ni}}$ 五个半衰期长的核素。如图9 (f)所示,对于考虑的五个核素(除$ ^{58} {\rm{Ni}}$ 外),我们依然得到相似的结果,相较图9 (d) (e),$ ^{56-72} {\rm{Ni}}$ +$ ^{238} {\rm{U}}$ 整个同位素链,截面结果存在1~2个量级的差异。鉴于俘获截面并非是影响生成截面的重要因素,我们重点分析熔合和存活两个过程。
-
对
$ ^{54-62} {\rm{Mn}}$ +$ ^{243} {\rm{Am}}$ 反应系统,图10 (a)大致显示出最大值条件下,各反应道的熔合几率随弹核质量数的增加呈先指数级下降后趋于平缓的变化关系,最大可相差2个数量级。图10 (b)给出了反应系统各反应道的内部熔合位垒,其中对$ ^{54} {\rm{Mn}}$ 、$ ^{55} {\rm{Mn}}$ 两个稳定核素用阴影标注,大致可看出内部熔合位垒随弹核质量数先呈现上升趋势,后趋于平缓变化,且在弹核质量数$ A $ = 56之后有1 MeV左右的波动。一方面,随着弹核质量数
$ A $ 的增加,反应系统的质量不对称度$ \eta $ 逐渐减小,势能曲面上超重复合核演化路径更远,降低了熔合几率。另一方面,在$ ^{54} {\rm{Mn}}$ 、$ ^{55} {\rm{Mn}}$ 弹核附近的结果表明,内部熔合位垒直接影响了熔合几率。当内部熔合位垒越小时,熔合几率越大,越有利于复合核合成,进而说明了熔合几率是弹核同位素依赖性的重要因素。同时,我们给出了
$ ^{52-59} {\rm{Cr}}$ +$ ^{243} {\rm{Am}}$ 反应系统的计算结果。由图11 (a)可知,反应系统$ ^{52-59} {\rm{Cr}}$ +$ ^{243} {\rm{Am}}$ 与$ ^{54-62} {\rm{Mn}}$ +$ ^{243} {\rm{Am}}$ 趋势相同,对三个稳定核素$ ^{52} {\rm{Cr}}$ 、$ ^{53} {\rm{Cr}}$ 、$ ^{54} {\rm{Cr}}$ ,熔合几率呈显著的下降趋势;图11 (b)为内部熔合位垒随弹核质量数的变化关系,考虑到弹核同位素的奇偶效应,对于偶偶核$ ^{52} {\rm{Cr}}$ 、$ ^{54} {\rm{Cr}}$ ,前面的结论依然适用。对
$ ^{56-72} {\rm{Ni}}$ +$ ^{238} {\rm{U}}$ 反应系统,由图12 (a)可得相应的熔合几率随弹核质量数增加呈先指数级下降而后趋于平缓的变化趋势,跨越的数量级大,因此进一步印证了熔合几率是同位素依赖性的重要因素。在图12 (b)中内部熔合位垒与弹核质量数呈现出更复杂的变化关系:选取偶偶核$^{56,\,58,\,62} {\rm{Ni}}$ 可看出其内部熔合位垒逐渐增大,相应的熔合几率呈指数级下降;而对稳定偶偶核$ ^{62} {\rm{Ni}}$ 、$ ^{64} {\rm{Ni}}$ ,尽管$ ^{64} {\rm{Ni}}$ 内部熔合位垒相较低2.32 MeV,但熔合几率近似相等,结果表明:一方面,$ ^{64} {\rm{Ni}}$ +$ ^{238} {\rm{U}}$ 反应道有更低的内部熔合位垒,越容易合成超重复合核,由此成为合成Z = 120超重核的一个重要弹靶组合;另一方面,说明除了考虑内部熔合位垒,其他影响如弹靶取向效应等也需要考虑。 -
在双核模型中,影响存活几率的因素主要是中子分离能[48]与裂变位垒。其中,中子分离能可表示为
$$ B^{\rm{n}} = \Delta(Z,\,A-1)+\Delta(n)-\Delta(Z,\,A), $$ (23) $ \Delta(n) $ 表示中子的质量过剩,$\Delta(Z,\,A)$ 表示质子数为$ Z $ 、质量数为$ A $ 的原子核的质量过剩,$\Delta(Z,\,A-1)$ 同理。中子分离能越小,对相应的中子蒸发道越有利。当Z = 104左右时,裂变位垒趋近零,自发裂变寿命小于
$ {10}^{-14} $ s,按照经典液滴模型理论预言,超重核的液滴能为零。壳结构效应[2]使得超重核得以存在,几个MeV的壳修正能,增加了原子核的裂变位垒,延缓其自发裂变[5]。因此,较高的裂变位垒更有利于超重核存活。裂变位垒的大小是存活几率计算非常敏感的一个重要参数,目前它的准确描述还不清楚,只能通过约束以描述其数值大小的趋势。通常,它的经验取值是通过宏观液滴能部分减去微观壳修正部分得到,表示为$$ B^{\rm{f}} = B_{\rm{LD}}-\Delta E_{\rm{shell}},\, $$ (24) 其中:
$ B_{\rm{LD}} $ 表示液滴能;$ \Delta E_{\rm{shell}} $ 表示壳修正能。本文所关注的是超重核的合成,其液滴能部分的势垒高度按零近似,壳修正能量取自质量表FDRM2012[49],即壳修正的绝对值作为裂变位垒大小。尽管不同的模型给出壳修正能量不同,但大体趋势特征近乎一致。而存活几率的大小还有其他不确定参数的影响,比如能级密度参数、裂变粘滞系数、能级密度的集体增强因子等,这些都会影响存活几率的数值计算,但对其趋势的描述几乎不会太大改变,故通过这些参数的合理选取并固定,在能较好地再现实验数据的同时,也能给出相对可靠的理论预言。
结合图9给出的最大生成截面与弹核质量数的变化关系,我们同样对存活阶段做相关讨论。对
$ ^{54-62} {\rm{Mn}}$ +$ ^{243} {\rm{Am}}$ 反应系统,各反应道在不同中子蒸发道取得最大值。其中,$ ^{54-59} {\rm{Mn}}$ +$ ^{243} {\rm{Am}}$ 均在3n蒸发道取得最大值,图中用黑、红两种颜色表示。从图13 (a)中可知,$ ^{54-62} {\rm{Mn}}$ +$ ^{243} {\rm{Am}}$ 各反应道的存活几率相差不大,在一个数量级范围内。图13 (b) (c)给出了该反应系统各反应道对应的中子分离能与裂变位垒,对$ ^{54} {\rm{Mn}}$ (奇奇核)、$ ^{55} {\rm{Mn}}$ (奇偶核)两个稳定核素,相邻的核素$ ^{56} {\rm{Mn}}$ (奇奇核)、$ ^{57} {\rm{Mn}}$ (奇偶核)的裂变位垒稍低于$ ^{54} {\rm{Mn}}$ 、$ ^{55} {\rm{Mn}}$ ,使得$ ^{56} {\rm{Mn}}$ 、$ ^{57} {\rm{Mn}}$ 在退激发的过程中相较更容易越过裂变位垒,发生裂变,导致存活几率降低;但$ ^{56} {\rm{Mn}}$ 、$ ^{57} {\rm{Mn}}$ 的中子分离能更小,使得中子蒸发宽度增加,增加了存活几率,进而导致$ ^{56} {\rm{Mn}}$ 、$ ^{57} {\rm{Mn}}$ 的存活几率分别大于$ ^{54} {\rm{Mn}}$ 、$ ^{55} {\rm{Mn}}$ 。由此说明:存活过程中子蒸发与裂变存在竞争关系。一般考虑奇偶效应的情况下,当反应道裂变位垒越高,中子分离能越小时,超重核存活几率越大。$ ^{52-59} {\rm{Cr}}$ +$ ^{243} {\rm{Am}}$ 反应系统均在3n蒸发道处取得最大值,用黑、红两种颜色表示。图14 (a)给出了$ ^{52-59} {\rm{Cr}}$ +$ ^{243} {\rm{Am}}$ 相应各反应道最大值处存活几率随弹核质量数的变化关系,即$ ^{52-59} {\rm{Cr}}$ +$ ^{243} {\rm{Am}}$ 反应系统大致呈先呈指数级上升后平缓波动的趋势,且在A = 59处取得极大值。图14 (b) (c)表示$ ^{52-59} {\rm{Cr}}$ +$ ^{243} {\rm{Am}}$ 各反应道的中子分离能、裂变位垒与复合核质量数的关系,对稳定核素$ ^{52} {\rm{Cr}}$ (偶偶核)、$ ^{53} {\rm{Cr}}$ (偶奇核)、$ ^{54} {\rm{Cr}}$ (偶偶核),$ ^{54} {\rm{Cr}}$ +$ ^{243} {\rm{Am}}$ 的中子分离能在3n处最小,裂变位垒最高,相应的存活几率最大。在A = 58处相应的中子分离能取得极小值,却同$ ^{54} {\rm{Cr}}$ 、$ ^{56} {\rm{Cr}}$ 处的存活几率相近,原因在于裂变位垒较$ ^{54} {\rm{Cr}}$ 、$ ^{56} {\rm{Cr}}$ 小,进一步印证了中子蒸发与裂变位的竞争关系。相对而言,稳定核素$ ^{54} {\rm{Cr}}$ 有较高的裂变位垒,较小的中子分离能,更有利于合成超重复合核。对
$ ^{56-72} {\rm{Ni}}$ +$ ^{238} {\rm{U}}$ 反应系统,黑、红两种颜色表示反应道在3n蒸发道处取得最大值。以图15 (a)偶奇核为例,各反应道存活几率随弹核质量数增加呈先指数级上升后下降的趋势。图15 (b) (c)给出了各反应道弹核质量数与中子分离能、裂变位垒的关系。选取各反应道在3n蒸发道处的结果讨论,对比稳定核素$ ^{58} {\rm{Ni}}$ (偶偶核)、$ ^{62} {\rm{Ni}}$ (偶偶核),$ ^{62} {\rm{Ni}}$ +$ ^{238} {\rm{U}}$ 相较有较小的中子分离能,较高的裂变位垒,因此存活几率相较$ ^{58} {\rm{Ni}}$ +$ ^{238} {\rm{U}}$ 约高1个量级;对比$ ^{63} {\rm{Ni}}$ (偶奇核)、$ ^{65} {\rm{Ni}}$ (偶奇核)、$ ^{69} {\rm{Ni}}$ (偶奇核)、$ ^{71} {\rm{Ni}}$ (偶奇核),随着弹核包含的质量数越多,裂变位垒明显下降,增加了裂变宽度,反应系统越容易发生裂变,相应存活几率降低。由此说明:对于$ ^{56-72} {\rm{Ni}}$ +$ ^{238} {\rm{U}}$ 反应系统,裂变位垒高度同弹核同位素链具有强烈的依赖性。 -
通过4.3节的讨论,各反应系统熔合几率随弹核质量数的增加整体呈现先指数级下降后平缓的趋势,且在指数级下降阶段,大致可以得到内部熔合位垒逐步增加的结论。由此可以说明熔合几率(驱动势)直接影响弹核同位素依赖性。通过4.4节的讨论,存活过程涉及中子蒸发与裂变的竞争关系,与之对应的裂变位垒、中子分离能需要重点关注。对照三条弹核同位素链存活几率的计算结果,存活几率具有弹核同位素依赖性。
由于熔合阶段和存活阶段直接影响最终生成截面,为比较两个阶段的影响大小,我们计算了各同位素链稳定弹核对应的截面,如图16所示。
对
$^{54,\,55} {\rm{Mn}}$ +$ ^{243} {\rm{Am}}$ 、$^{52,\,53,\,54} {\rm{Cr}}$ +$ ^{243} {\rm{Am}}$ 两个反应系统,在激发能一定的情况下,相邻弹核同位素熔合几率相差1~2个量级,存活几率大致在1个量级范围内,最终生成截面相差1个量级左右,由此可见:熔合过程影响生成截面作用更大,即驱动势(内部熔合势垒)才是决定$^{54,\,55} {\rm{Mn}}$ +$ ^{243} {\rm{Am}}$ 、$^{52,\,53,\,54} {\rm{Cr}}$ +$ ^{243} {\rm{Am}}$ 最大生成截面的重要反应机制。对
$^{58,\,61,\,62,\,63,\,64} {\rm{Ni}}$ +$ ^{238} {\rm{U}}$ 反应系统,结果有些新的变化,如图16 (c) (f),$^{58,\,61,\,62,\,63,\,64} {\rm{Ni}}$ +$ ^{238} {\rm{U}}$ 随激发能的不断增大,熔合、存活过程对生成截面影响存在一个竞争关系,具体来说:在$ E^* <38$ MeV时,存活过程占主导作用,需要重点考虑中子分离能、裂变位垒等相关物理量;在$ E^* <38$ MeV时,熔合过程则占主导作用。对于不同激发能大小,计算结果也会出现大的变化,本文不再讨论。综上,系统激发能大小、驱动势(内部熔合位垒)对于一个反应系统需要重点讨论。
-
摘要: 合成Z=119,120超重核是当今各核物理实验室争相追逐的目标,理论预言可靠的弹靶组合、入射能等信息有助于超重核合成的实验设计和探测。本工作基于双核模型研究影响重离子核反应生成截面大小的反应机制,计算了
$^{50}{\rm{Ti}}$ +$^{249}{\rm{Bk}}$ 、$^{50}{\rm{Ti}}$ +$^{249}{\rm{Cf}}$ 两个弹靶组合,预测$^{50}{\rm{Ti}}$ +$^{249}{\rm{Bk}}$ 的生成截面为0.021 1 pb。考虑双核系统熔合与存活两个过程,重点关注$^{52-59}{\rm{Cr}}$ +$^{243}{\rm{Am}}$ 、$^{54-62}{\rm{Mn}}$ +$^{243}{\rm{Am}}$ 、$^{56-72}{\rm{Ni}}$ +$^{238}{\rm{U}}$ 生成截面的同位素链依赖性,研究表明熔合几率随弹核质量数呈现强烈的依赖行为,直接影响蒸发剩余截面大小。Abstract: Synthesis Z=119, 120 superheavy nuclei is the goal that nuclear physics laboratories are chasing. Theoretically, it is helpful to design experiments and detect newly synthesized superheavy nuclei by giving reliable information such as projectile-target combination, incident energy, etc. Based on dinuclear system model, we explore the reaction mechanism affecting the generation cross section of the heavy ion nuclear reaction and calculate$^{50}{\rm{Ti}}$ +$^{249}{\rm{Bk}}$ ,$^{50}{\rm{Ti}}$ +$^{249}{\rm{Cf}}$ projectile-target combination. Among them, we predicted that$^{50}{\rm{Ti}}$ +$^{249}{\rm{Bk}}$ generation cross section is 0.021 1 pb. Considering the fusion process and survival process, we focus on the isotopic dependence of the projectile isotopes, such as$^{52-59}{\rm{Cr}}$ +$^{243}{\rm{Am}}$ ,$^{54-62}{\rm{Mn}}$ +$^{243}{\rm{Am}}$ ,$^{56-72}{\rm{Ni}}$ +$^{238}{\rm{U}}$ . It is shown that the fusion probability is significantly dependent on the projectile nuclear mass number, which directly affects evaporation residual cross section. -
图 1
$^{48}{\rm{Ca}}$ +$^{238}{\rm{U}}$ 相互作用势随两核中心距离$R$ 的变化曲线[29] -
[1] 李文飞. 基于重离子输运理论对超重核合成及非对称核物质状态方程的研究[D]. 兰州: 中国科学院近代物理研究所, 2004. LI W F. On the Synthesis of the Superheavy Nuclei and Asymmetric Nuclear Matter Equation of State Based on the Transport Theory[D]. Lanzhou: Institute of Modern Physics, Chinese Academy of Sciences, 2004. (in Chinese) [2] MYERS W D, SWIATECKI W J. Nucl Phys, 1966, 81(2): 1. doi: 10.1016/s0029-5582(66)80001-9 [3] NILSSON S G, TSANG C F, SOBICZEWSKI A, et al. Phys A, 1969, 131(1): 1. doi: 10.1016/0375-9474(69)90809-4 [4] 罗亦孝. 科学通报, 2016, 61(21): 2326. doi: 10.1360/N972016-00710 LUO Y X. Chin Sci Bull, 2016, 61(21): 2326. (in Chinese) doi: 10.1360/N972016-00710 [5] 王艺澎, 郭树青, 包小军, 等. 物理学进展, 2021, 41(04): 157. doi: 10.13725/j.cnki.pip.2021.04.001 WANG Y P, GUO S Q, BAO X J, et al. Progress in Physics, 2021, 41(04): 157. (in Chinese) doi: 10.13725/j.cnki.pip.2021.04.001 [6] HOFMANN S, HEINZ S, MANN R, et al. Eur Phys J A, 2016, 52(6): 180. doi: 10.1140/epja/i2016-16180-4 [7] 卢希庭. 原子核物理[M]. 北京: 原子能出版社, 2000. LU X T. Nuclear Physics[M]. Beijing: Atomic Energy Press, 2000. (in Chinese) [8] ZHANG Z Y, MA L, GAN Z G, et al. Nucl Instr and Meth, B, 2013, 317(12): 315. doi: 10.1016/j.nimb.2013.05.062 [9] 徐瑚珊, 周小红, 肖国青, 等. 原子核物理评论, 2003, 02(2): 76. doi: 10.11804/NuclPhysRev.20.02.076 XU H S, ZHOU X H, XIAO G Q, et al. Nucl Phys Rev, 2003, 02(2): 76. (in Chinese) doi: 10.11804/NuclPhysRev.20.02.076 [10] HAMILTON J H, HOFMANN S, OGANESSIAN Y T, et al. Annu Rev Nucl Part Sci, 2013, 63: 383. doi: 10.1146/annurev-nucl-102912-144535 [11] MORITA K. Nuclear Physics A, 2015, 944: 30. doi: 10.1016/j.nuclphysa.2015.10.007 [12] ZHANG Z Y, GAN Z G, YANG H B, et al. Phys Rev Lett, 2019, 122(19): 192503. doi: 10.1103/PhysRevLett.122.192503 [13] MA L, ZHANG Z Y, GAN Z G, et al. Phys Rev Lett, 2020, 125(3): 032502. doi: 10.1103/PhysRevLett.125.032502 [14] SWIATECKI W J. Phys Scr, 1981, 24(1B): 113. doi: 10.1088/0031-8949/24/1B/007 [15] ARITOMO Y, WADA T, OHTA M, et al. Phys Rev C, 1999, 59(2): 796. doi: 10.1103/PhysRevC.59.796 [16] ZAGREBAEV V I. Phys Rev C, 2001, 64(3): 034606. doi: 10.1103/PhysRevC.64.034606 [17] ADAMIAN G G, ANTONENKO N V, SCHEID W, et al. Nucl Phys A, 1998, 633: 409. doi: 10.1016/s0375-9474(98)00124-9 [18] LI W F, WANG N, LI J F, et al. Europhys Lett, 2007, 64(6): 750. doi: 10.1209/epl/i2003-00622-0 [19] 左维, 李君清, 赵恩广. 原子核物理评论, 2006, 23(4): 8. doi: 10.3969/j.issn.1007-4627.2006.04.005 ZUO W, LI J Q, ZHAO E G. Nucl Phys Rev, 2006, 23(4): 8. (in Chinese) doi: 10.3969/j.issn.1007-4627.2006.04.005 [20] BOHR N, WHEELER J A. Phys Rev, 1939, 56(5): 426. doi: 10.1103/PhysRev.56.426 [21] WEISSKOPF V. Physical Review, 1937, 52(4): 295. doi: 10.1103/physrev.52.295 [22] LI W, WANG N, JIA F, et al. Journal of Physics G: Nuclear and Particle Physics, 2006, 32(8): 1143. doi: 10.1088/0954-3899/32/8/006 [23] BOHR A, MOTTELSON B R. Nuclear Structure: Volume II: Nuclear Deformations[M]. Singapore: World Scientific Publishing Company, 1998. [24] WONG C Y. Phys Rev lett, 1973, 31(12): 766. doi: 10.1103/PhysRevLett.31.766 [25] 徐躬耦, 王顺金. 原子核理论[M]. 北京: 高等教育出版社, 1992. XU G O, WANG S J. Nuclear Theory[M]. Beijing: Higher Education Press, 1992. (in Chinese) [26] FENG Z Q, JIN G M, FU F, et al. Nucl Phys A, 2006, 771: 50. doi: 10.10/j.nuclphysa.2006.03.002 [27] ADAMIAN G G, ANTONENKO N V, SCHEID W. Phys Rev C, 2003, 68(3): 034601. doi: 10.1103/PhysRevC.68.034601 [28] 冯兆庆. 近垒重离子熔合反应和超重核合成机制研究[D]. 兰州: 中国科学院近代物理研究所, 2007. FENG Z Q. Studies on Heavy-ion Fusion Reactions near Coulomb Barrier and Synthesis of Superheavy Nuclei[D]. Lanzhou: Institute of Modern Physics, Chinese Academy of Sciences, 2007. (in Chinese) [29] 黄明辉. 超重核合成机制-包含中子质子跃迁的二维主方程[D].兰州: 中国科学院近代物理研究所, 2008. HUANG M H. The Synthesis Mechanism of Super Heavy Nuclei- the Neutron and Proton Transfer in the Two-dimensional Master Equation[D]. Lanzhou: Institute of Modern Physics, Chinese Academy of Sciences, 2008. (in Chinese) [30] LI Q F, ZUO W, LI W F, et al. The European Physical Journal A - Hadrons and Nuclei, 2005, 24(2): 223. doi: 10.1140/epja/i2004-10138-1 [31] FENG Z Q, JIN G M, LI J Q, et al. Phys Rev C, 2007, 76(4): 044606. doi: 10.1103/PhysRevC.76.044606 [32] JACKSON J D. Can J Phys, 1956, 34(8): 767. doi: 10.1139/p56-087 [33] ADAMIAN G G, ANTONENKO N V, JOLOS R V, et al. International Journal of Modern Physics E, 1996, 5(01): 191. doi: 10.1142/S0218301396000098 [34] OGNAESSIAN Y T, UTYONKOV V K, LOBANOV Y V, et al. Physical Review C, 2006, 74(4): 044602. doi: 10.1103/PhysRevC.74.044602 [35] ZAGREBAEV V I, KARPOV A V, GREINER W. Physical Review C, 2012, 85(1): 014608. doi: 10.1103/physrevc.85.014608 [36] NASIRO A K, GIARDINA G, MANDAGLIO G, et al. Phys Rev C, 2009, 79(2): 024606. doi: 10.1103/PhysRevC.79.024606 [37] ADAMIAN G G, ANTONENKO N V, SCHEID W. Eur Phys J A, 2009, 41(2): 235. doi: 10.1140/epja/i2009-10795-4 [38] NASIROV A K, MANDAGLIO G, GIARDINA G. Phys Rev C, 2011, 84(4): 044612. doi: 10.1103/PhysRevC.84.044612 [39] LIU Z H, BAO J D. Phys Rev C, 2011, 84(3): 031602. doi: 10.1103/PhysRevC.84.031602 [40] LIU Z H, BAO J D. Phys Rev C, 2009, 80(5): 054608. doi: 10.1103/PhysRevC.80.054608 [41] WANG N, ZHAO E G, SCHEID W, et al. Phys Rev C, 2012, 85(4): 041601. doi: 10.1103/PhysRevC.85.041601 [42] ADAMIAN G G, ANTONENKO N V, DIAZ-TORRES A, et al. Eur Phys J A, 2020, 56(2): 47. doi: 10.1140/epja/s10050-020-00046-7 [43] 杨秀秀, 张根, 李静静, 等. 原子核物理评论, 2020, 37(2): 151. doi: 10.11804/NuclPhysRev.37.2020005 YANG X X, ZHANG G, LI J J, et al. Nuclear Physics Review, 2020, 37(2): 151. (in Chinese) doi: 10.11804/NuclPhysRev.37.2020005 [44] ZAGREBAEV V, GREINER W. Phys Rev C, 2008, 78(3): 034610. doi: 10.1103/PhysRevC.78.034610 [45] LIU L, SHEN C W, LI Q F, et al. Eur Phys J A, 2016, 52(2): 35. doi: 10.1140/epja/i2016-16035-0 [46] 刘硕, 刘玲. 沈阳师范大学学报: 自然科学版, 2014, 32(2): 237. doi: 10.3969/j.issn.1673-5862.2014.02.024 LIU S, LIU L. Journal of Shenyang Normal University (Natural Science Edition), 2014, 32(2): 237. (in Chinese) doi: 10.3969/j.issn.1673-5862.2014.02.024 [47] NIU F, CHEN P H, FENG Z Q. Nucl Sci Tech, 2021, 32: 103. doi: 10.1007/s41365-021-00946-3 [48] WANG N, ZHAO K, SCHEID W, et al. Phys Rev C, 2008, 77(1): 014603. doi: 10.1103/PhysRevC.77.014603 [49] MOLLER P, SIERK A J, ICHIKAWA T, et al. Atomic Data and Nuclear Data Tables, 2016, 109. doi: 10.1016/j.adt.2015.10.002 [50] 包小军. 超重核合成与衰变性质的理论研究[D]. 兰州: 兰州大学, 2016. BAO X J. Theoretical Studies of Superheavy Nuclei Synthesis and Decay Properties[D]. Lanzhou: Lanzhou University, 2016. (in Chinese)