1.School of Nuclear Science and Technology, University of South China, Hengyang 421001, China 2.College of Physics and Electronics, Central South University, Changsha 410083, China 3.Cooperative Innovation Center for Nuclear Fuel Cycle Technology & Equipment, University of South China, Hengyang 421001, China 4.Key Laboratory of Low Dimensional Quantum Structures and Quantum Control, Hunan Normal University, Changsha 410081, China Received Date:2020-06-29 Available Online:2021-02-15 Abstract:In the present work, a two-parameter empirical formula is proposed, based on the Geiger-Nuttall law, to study two-proton ($ 2p $) radioactivity. Using this formula, the calculated $ 2p $ radioactivity half-lives are in good agreement with the experimental data as well as with calculated results obtained by Goncalves et al. [Phys. Lett. B774, 14 (2017)] using the effective liquid drop model (ELDM), Sreeja et al. [Eur. Phys. J. A 55, 33 (2019)] using a four-parameter empirical formula, and Cui et al. [Phys. Rev. C101: 014301 (2020)] using a generalized liquid drop model (GLDM). In addition, this two-parameter empirical formula is extended to predict the half-lives of 22 possible $ 2p $ radioactivity candidates with $ 2p $ radioactivity released energy $ Q_{2p} > 0 $, obtained from the latest evaluated atomic mass table AME2016. The predicted results are highly consistent with those obtained using other theoretical models such as the ELDM, GLDM and four-parameter empirical formula.
HTML
--> --> -->
II.THEORETICAL FRAMEWORKIn 1911, Geiger and Nuttall found there is a phenomenological relationship between the $ \alpha $ decay half-life $ T_{1/2} $ and the decay energy $ Q_{\alpha} $. This relationship is the so-called Geiger-Nuttall (G-N) law. It is expressed as:
where a and b represent the two isotopic chain–dependent parameters of this formula. Later, the G-N law was widely applied to study the half-lives of $ \alpha $ decay [38, 46-48], cluster radioactivity [49-51] and proton radioactivity [52-54]. However, relative to $ \alpha $ decay and cluster radioactivity, the proton radioactivity half-life is more sensitive to the centrifugal barrier. This means that the linear relationship between the half-life of the proton radioactivity and the released energy $ Q_p $ only exists for proton-radioactive isotopes with the same orbital angular momentum l taken away by the emitted proton [44, 52, 54]. Similarly, the $ 2p $ radioactivity half-life may also depend strongly on the $ 2p $ radioactivity released energy $ Q_{2p} $ and the orbital angular momentum l taken away by the two emitted protons. Recently, considering the contributions of $ Q_{2p} $ and the orbital angular momentum l to the $ 2p $ radioactivity half-life, Sreeja et al. put forward a four-parameter empirical formula to study the $ 2p $ radioactivity half-lives, which is expressed as [36]
where a = 0.1578, b = 1.9474, $ c = -1.8795 $, and $ d = -24.847 $ denote the adjustable parameters, which are obtained by fitting the calculated results of the ELDM [33]. Their calculated results can reproduce the known experimental data well. In our previous work [44], considering the contributions of the daughter nuclear charge $ Z_d $ and the orbital angular momentum l taken away by the emitted proton to the proton radioactivity half-life, we proposed a two-parameter empirical formula for a new G-N law for proton radioactivity. This formula is written as:
where $ a_{\beta} = 0.843 $ and $ b_{\beta} = -27.194 $ are the fitted parameters. The exponent on the orbital angular momentum l taken away by the emitted proton, $ \beta $ , is 1, which is obtained by fitting 44 experimental data points of proton radioactivity in the ground state and isomeric state. Combined with the work from Sreeja et al. [36, 54] and our previous work [44], it is interesting to examine whether or not a two-parameter form of the empirical formula is suitable to investigate $ 2p $ radioactivity. In this work, because there are no experimental data for $ 2p $ radioactive nuclei with orbital angular momentum $ l \ne 0 $, we choose the experimental data for true $ 2p $ radioactive nuclei (19Mg, 45Fe, 48Ni, 54Zn and 67Kr) with l = 0, and the predicted $ 2p $ radioactivity half-lives of 7 nuclei with $ l \ne 0 $ (1 case with l = 1, 4 cases with l = 2 and 2 cases with l = 4) are extracted from Goncalves et al. [33]. First, for the $ \beta $ value describing the effect of l on the $ 2p $ radioactivity half-life, we choose the $ \beta $ value corresponding to the smallest standard deviation $ \sigma $ between the database and the calculated $ 2p $ radioactivity half-lives as the optimal value, with $ \beta $ varying from 0.1 to 0.5. The relationship between the $ \sigma $ and $ \beta $ values is shown in Fig. 1. It is clear that $ \sigma $ is smallest when $ \beta $ is equal to 0.25. Comparing with the $ \beta $ value of Eq. (3) reflecting the effect of l on the proton radioactivity half-life, this $ \beta $ value is smaller. The reason may be that the reduced mass $ \mu $ of a proton-radioactive nucleus is smaller than that of a $ 2p $-radioactive nucleus, leading to the contribution of the centrifugal barrier to the half-life of the $ 2p $-radioactive nucleus being smaller. Correspondingly, the values of parameters a and b are given as: Figure1. (color online) Relationship between the standard deviation $\sigma$ and the value of β.
$ a = 2.032,\,\;\; b = 26.832, $
(4)
Then, we can obtain a final formula, which can be written as:
III.RESULTS AND DISCUSSIONThe primary aim of this work is to verify the feasibility of using Eq. (5) to investigate $ 2p $ radioactivity. The calculated logarithmic half-lives of $ 2p $-radioactive nuclei are listed in the seventh column of Table 1. Meanwhile, for comparison, the calculated results using GLDM, ELDM and a four-parameter empirical formula are shown in the fourth to sixth column of this table, respectively. In Table 1, the first three columns denote the $ 2p $-radioactive nucleus, the experimental $ 2p $ radioactivity released energy $ Q_{2p} $ and the logarithmic experimental $ 2p $ radioactivity half-life $ {\rm{log}}_{10}T_{1/2}^{\rm{exp}} $, respectively. For quantitative comparisons between the calculated $ 2p $ radioactivity half-lives using our empirical formula and the experimental results, the last column gives the logarithm of errors between the experimental $ 2p $ radioactivity half-lives and those calculated using our empirical formula $ {\rm{log_{10}}}HF = {\rm{log}}_{10}{{T}_{1/2}^{\,\rm{exp}}} -{\rm{log}}_{10}{{T}_{1/2}^{\,\rm{cal}}} $. From this table, it can be seen that for the true $ 2p $-radioactive nuclei 19Mg, 45Fe, 48Ni, 54Zn and 67Kr ($ Q_p < 0 $, $ Q_{2p} > 0 $), most values of $ {\rm{log_{10}}}HF $ are between -1 and 1. Particularly, for the cases of 48Ni, with $ Q_{2p} $ = 1.290, and 45Fe, with $ Q_{2p} $ = 1.154, the values of $ {\rm{log_{10}}}HF $ are 0.07 and 0.24, indicating our calculated results can reproduce the experimental data well. As for the sequential or pseudo-$ 2p $-radioactive nuclei 6Be, 12O and 16Ne ($ Q_p > 0 $, $ Q_{2p} > 0 $), the values of $ {\rm{log_{10}}}HF $ for 6Be and 16Ne are relatively large. Likewise, the differences between the experimental data and the calculated $ 2p $ radioactivity half-lives using GLDM, ELDM and the four-parameter empirical formula are more than three orders of magnitude. This may be due to the limitations of the early experimental equipment, resulting in the measured decay widths of these $ 2p $ radioactivity nuclei not being accurate enough. It would be helpful to measure the experimental $ 2p $ half-lives of these nuclei again in the future. In the case of 12O, the values of $ {\rm{log_{10}}}HF $ are small, implying that our formula may also be suitable for studying pseudo-$ 2p $-radioactive nuclei which have relatively accurate experimental data.
Table1.Comparison of the experimental data for $ 2p $-radioactive nuclei with different theoretical models (GLDM, ELDM, the four-parameter empirical formula of Ref. [36] and our empirical formula. Experimental data are taken from the corresponding references.
To further test the feasibility of our empirical formula, we also use Eq. (5) to predict the $ 2p $ radioactivity half-lives of 22 nuclei with $ 2p $ radioactivity released energy $ Q_{2p} > 0 $. The $ Q_{2p} $ values are taken from the latest evaluated atomic mass table AME2016 and shown in the second column of Table 2. In this table, the first and third columns give the $ 2p $ radioactivity candidates and the angular momentum l taken away by the two emitted protons, respectively. For a benchmark, the predicted results using GLDM, ELDM and the four-parameter empirical formula, extracted from Refs. [37], [33] and [36] respectively, are also listed in this table. We can clearly see that for $ l\ne 0 $, the predicted results using our empirical formula are closer to those predicted using ELDM than those predicted using the four-parameter empirical formula. Most of the predicted results are of the same order of magnitude. As an example, in the cases of 28Cl (60As), the predicted $ 2p $ radioactivity half-lives using ELDM, the four-parameter empirical formula and our empirical formula are $ -12.95 $ ($ -8.68 $), $ -14.52 $ ($ -10.84 $) and $ -12.46 $ ($ -8.33 $), respectively. This implies that our empirical formula is also suitable for studying nuclei with orbital angular momentum $ l\neq0 $. In the case of $ l = 0 $, the predicted $ 2p $ radioactivity half-lives using our empirical formula are in good agreement with those from GLDM and ELDM. To further demonstrate the significant correlation between the $ 2p $ radioactivity half-lives $ T_{1/2} $ and the $ 2p $ radioactivity released energies $ Q_{2p} $, based on Eq. (5), we plot the quantity $ [{\rm{log_{10}}}{T}_{1/2} + 26.832]/(Z_{d}^{0.8}+l^{\,0.25}) $ as a function of $ Q_{2p}^{-1/2} $ in Fig. 2. In this figure, there is an obvious linear dependence of $ {\rm{log}}_{10}T_{1/2} $ on $ Q_{2p} $$ ^{-1/2} $, while the contributions of charge number $ Z_{d} $ and orbital angular momentum l on the $ 2p $ radioactivity half-lives are removed.
Table2.Comparison of calculated $ 2p $ radioactivity half-lives using GLDM, ELDM, the four-parameter empirical formula from Ref. [36] and our empirical formula. The $ 2p $ radioactivity released energy $ Q_{2p} $ and orbital angular momentum l taken away by the two emitted protons are taken from Ref. [33].
Figure2. (color online) The linear relationship between the quantity $[{\rm{log_{10}}}{T}_{1/2} + 26.832]/(Z_{d}^{0.8}+l^{\,0.25})$ and $Q_{2p}$ for the database used to fit the parameters of Eq. (5).