-
The
$ 2p $ radioactive half-life$ T_{1/2} $ can be calculated using the decay constant λ and expressed as$ T_{1/2} = \frac{ {\rm ln}2}{\lambda} $
(1) with
$ \lambda = S_{2p} \nu P . $
(2) Here,
$ S_{2p} = G_0^2[A/(A-2) ]^{2n}{\rm{\rho}}^2 $ denotes the spectroscopic factor of$ 2p $ radioactivity estimated using the cluster overlap approximation [20] with$ G_0^2 = (2n)!/[2^{2n}(n!)^2] $ [51], where n$ \approx \; (3Z)^{1/3}-1 $ is the average principal proton oscillator quantum number [52], A and Z denote the mass number and proton number of the parent nucleus, respectively,$ {\rm{\rho}}^2 = 0.015 $ denotes the proton overlap function, which is determined by a$ {\rm{\chi}}^2 $ optimization from the experimental half-lives of$ ^{19}{\rm{Mg}} $ ,$ ^{45}{\rm{Fe}} $ ,$ ^{48}{\rm{Ni}} $ , and$ ^{54}{\rm{Zn}} $ [16], and ν is the assault frequency, which can be calculated by the oscillation frequency ω and expressed as [53]$ \nu = \frac{\omega}{2\pi} = \frac{\Bigg(2n_r + l +\dfrac{3}{2}\Bigg)\hbar}{1.2\pi \mu R_n^2} = \frac{\Bigg(G+\dfrac{3}{2}\Bigg)\hbar}{1.2\pi \mu R_{00}^2} , $
(3) where
$ R_n = \sqrt{\dfrac{3}{5}}R_{00} $ is the nucleus root-mean-square (rms) radius with$ R_{00} = 1.24A^{1/3} \left(1+ \dfrac{1.646}{A}- 0.191\dfrac{A-2Z}{A}\right) $ [54], and$ G = 2n_r + l $ is the principal quantum number, with$ n_r $ and l being the radial and angular momentum quantum numbers, respectively. For$ 2p $ radioactivity, we set G = 4 or 5 corresponding to the$ 4\hbar \omega $ or$ 5\hbar \omega $ oscillator shells, depending on the individual two-proton emitters. The reduced mass of the decaying nuclear system is given by$ \mu = m_d m_{2p} /( m_d + m_{2p} ) \simeq 938.3 \times 2 \times (A-2) $ /A MeV/$ c^2 $ , with$ m_d $ and$ m_{2p} $ being the mass of daughter nucleus and the two emitted protons, respectively. Finally,$\hbar $ is the reduced Planck constant.The quantity P given in Eq. (2) is the penetration probability of the two emitted protons crossing the barrier. It can be calculated by the semi-classical WKB approximation [55, 56] and written as
$ P = {\rm{exp}} \Bigg \{-\frac{2}{\hbar} \int_{R_t}^{R_{{\rm out}}} \sqrt{2\mu[V (r)-Q_{2p}]} {\rm{dr}} \Bigg \} \ , $
(4) where
$ V(r) $ is the diproton-daughter nucleus interaction potential,$ R_{\rm out} $ is the outer turning point of potential barrier satisfying the condition$ V(R_{\rm out}) = Q_{2p} $ , and$ R_t = R_d + R_{2p} $ is the distance of the touching configuration, where$ R_d $ and$ R_{2p} $ are the radii of daughter nucleus and the two emitted protons, respectively. They can be calculated using [57]$ R = 1.28A^{1/3} - 0.76 + 0.8A^{-1/3}. $
(5) In general, the diproton-daughter nucleus electrostatic potential is by default Coulomb type, given as
$ V_{\rm C} = {Z_dZ_{2p}e^2}/{r}, $
(6) where
$ Z_d $ and$ Z_{2p} $ are the proton numbers of the daughter nucleus and the two emitted protons, respectively, and$ e^2 $ is the square of the electronic elementary charge. However, in the$ 2p $ radioactivity process, for the superposition of the involved charges, due to the magnetic field generated by the movement of the two emitted protons and the inhomogeneous charge distribution of the nucleus, the diproton-daughter nucleus electrostatic potential behaves as a Coulomb potential at a short distance and drops exponentially at large distance, i.e., the screened electrostatic effect [38, 58]. This behavior of electrostatic potential can be described using the Hulthen type potential and is defined as$ V_{\rm H} = \frac{aZ_dZ_{2p}e^2}{{\rm e}^{ar}-1}, $
(7) where a is the screening parameter. In this framework, the total diproton-daughter nucleus interaction potential
$ V(r) $ , shown in Fig. 1, is given byFigure 1. (color online) Schematic plot of the potential energy
$V(r)$ for$^{67}{\rm{Kr}}$ as a function of the distance between the centers of the daughter nucleus and the two emitted protons.$ \ V(r) = \left\{\begin{array}{*{20}{l}} -V_0, &0 \leqslant r \leqslant R_t ,\\ V_{\rm H}(r) + V_l(r), &r \geqslant R_t ,\\ \end{array}\right. $
(8) where
$ V_0 $ is the depth of square well and$ V_l $ (r) is the centrifugal potential.As
$ l (l + 1) \to \left(l + \dfrac{1}{2}\right)^2 $ is a necessary correction for one-dimensional problems [59], the centrifugal potential is chosen as the Langer modified centrifugal potential in this study. It can be expressed as$ V_l(r) = \frac{\hbar^2(l+\frac{1}{2})^2}{2\mu r^2}, $
(9) where l is the orbital angular momentum taken away by the two emitted protons. It can be obtained by the spin and parity conservation laws.
Obviously, the outer turning point of potential barrier
$ R_{\rm out} $ given in Eq. (4) is an important quantity in evaluating the half-life of$ 2p $ radioactivity. To obtain an analytic expression for$ R_{\rm out} $ , the centrifugal term should be approximated as [60]$ \frac{1}{r^2} \approx \frac{a^2}{({\rm e}^{ar}-1)^2} $
(10) with the value of a being small. Under this assumption, the outer turning point of the potential barrier
$ R_{\rm out} $ can be expressed as$ R_{\rm out} = \frac{1}{a}{\rm{ln}}\left[ \frac{2{\rm{V}}_1}{\sqrt{{{\rm{V}}_0}^2+4{\rm{V}}_1{\rm{Q}}_{{\rm{2p}}}}-{\rm{V}}_0}+1\right] $
(11) with
$ V_0 = 2ae^2Z_1,\;\;\;\; \ V_1 = \frac{a^2\hbar^2\left(l+\dfrac{1}{2}\right)^2}{2\mu} . $
(12) We define the Gamow factor as
$ G_{se} $ = -$ \log{P} $ for the outer potential barrier region. It can be given by the definite integral as [36]$ G_{se} = \frac{1}{a}[I_1(r)+ I_2(r)] \vert^{R_{\rm out}}_{R_t}. $
(13) These two terms,
$ I_1(r) $ and$ I_2(r) $ , have the following expressions: [36, 37]$ \begin{aligned}[b] I_1(r) =& -\sqrt{V_1x^2 + V_0x - Q_{2p}} \\ &+ \sqrt{Q_{2p}}\arcsin{\left[\frac{xV_0 - 2Q_{2p}}{x\sqrt{4Q_{2p}V_1 + {V_0}^2}}\right]} \\ &- \frac{V_0}{2\sqrt{V_1}} \ln{\left[2\sqrt{\left(V_1x^2 + V_0x - Q_{2p} \right) } + V_0 +2V_1x \right]} , \end{aligned} $
(14) $ \begin{aligned}[b] I_2(r) =& \sqrt{V_1y^2 + U_0y - U_1} \\ &- \sqrt{U_1}\arctan{\frac{yU_0-2U_1}{2\sqrt{U_1(V_1y^2 + U_0y - U_1)}}} \\& + \frac{U_0}{2\sqrt{V_1}} \ln\left[ 2\sqrt{V_1(V_1y^2+U_0y-U_1)}\right.\\&+U_0+2V_1y\Bigg], \end{aligned} $
(15) where the following notations are employed:
$ \begin{aligned}[b]& x = ({\rm e}^{ar} - 1)^{-1},\;\;\; \ y = 1 + ({\rm e}^{ar} - 1)^{-1}, \\& U_0 = V_0 - 2V_1,\;\;\; \ U_1 = Q_{2p} + V_0 -V_1. \end{aligned} $
(16) -
In this study, we systematically investigate the
$ 2p $ radioactive half-lives of nuclei near or beyond the proton drip line. The experimental$ 2p $ radioactive half-lives and corresponding released energies are taken from the corresponding references of$ 2p $ radioactivity. With the Hulthen potential for the electrostatic barrier being considered, this model contains one adjustable parameter, i.e., the screening parameter a, which is obtained by fitting the true$ 2p $ radioactive nuclei of$ ^{19} {\rm{Mg}} $ ,$ ^{45} {\rm{Fe}} $ ,$ ^{48} {\rm{Ni}} $ ,$ ^{54} {\rm{Zn}} $ , and$ ^{67} {\rm{Kr}} $ , amounting to 10 experimental datasets. The standard deviation σ, describing the difference between the experimental data and calculated data, can be expressed as$ \sigma = \sqrt{\frac{1}{10}\sum\limits_{i = 1}^{10}\left[\log{\left(\frac{T^{\rm{cacl}}_{1/2_i}}{T^{\rm{expt}}_{1/2_i}}\right)}\right]^2} . $
(17) Here,
$ {T^{\rm{expt}}_{1/2_i}} $ and$ {T^{\rm{cacl}}_{1/2_i}} $ represent the experimental and calculated$ 2p $ radioactive half-lives for the i-th nucleus, respectively. Through minimizing σ, the screening parameter is determined as$ a = 2.9647\times 10^{-2} \;{\rm{fm}}^{-1} $ . The value of a is small, though it observably impacts the classical outer turning point$ R_{\rm out} $ obtained by Eq. (11), with the$ 2p $ radioactive half-life being sensitive to$ R_{\rm out} $ . For intuitively displaying the screening effect, in Fig. 2, we illustrate the difference between$ R^{\rm C}_{\rm out} $ -$ R^{\rm H}_{\rm out} $ against the proton number of the daughter nucleus for different$ 2p $ radioactivity released energies,$ Q_{2p} $ . Here,$ R^{\rm C}_{\rm out} $ and$ R^{\rm H}_{\rm out} $ represent the$ R_{\rm out} $ value calculated using pure Coulomb and Hulthen potentials, respectively, while the centrifugal potential contribution is not considered. From this figure, we can clearly see that the smaller$ 2p $ radioactivity released energy$ Q_{2p} $ and larger proton number of the daughter nucleus$ Z_d $ increase the difference between$ R^{\rm C}_{\rm out} $ and$ R^{\rm H}_{\rm out} $ .Figure 2. (color online) Difference between
$R^{\rm C}_{\rm out}$ and$R^{\rm H}_{\rm out}$ obtained by V(r) =$Q_{2p}$ when only considering the Coulomb potential. For calculating$R^{\rm C}_{\rm out}$ , the Coulomb potential is taken as the potential of a uniformly charged sphere, i.e.,$V_{\rm C}$ (r) =${Z_dZ_{2p}e^2}/{r}$ , while for$R^{\rm H}_{\rm out}$ , the Coulomb potential is taken as the Hulthen potential, i.e.,$V_{\rm H}({{r}})$ =$\dfrac{aZ_dZ_{2p}e^2}{{\rm e}^{ar}-1}$ , with a =$2.9647\times10^{-2} \;{\rm{fm}}^{-1}$ .For a more intuitive illustration of the relationship between the screening effect and
$ Q_{2p} $ , taking the true$ 2p $ radioactive nucleus$ ^{45} {\rm{Fe}} $ as an example, we plot the difference of$ R^{\rm C}_{\rm out} $ -$ R^{\rm H}_{\rm out} $ against$ Q_{2p} $ in Fig. 3. From this figure, we can see an approximately linear relation between$ R^{\rm C}_{\rm out} $ -$ R^{\rm H}_{\rm out} $ and$ Q_{2p} $ , indicating a strong correlation between the screening effect and$ 2p $ radioactivity released energy.Figure 3. (color online) Difference between
$R^{\rm C}_{\rm out}$ and$R^{\rm H}_{\rm out}$ plotted against the 2p radioactivity released energy$Q_{2p}$ for$^{45}{\rm{Fe}}$ .In the following, using this model and the value of parameter a, we systematically calculate the
$ 2p $ radioactive half-lives of nuclei using the experimental data, including true and not true$ 2p $ radioactive nuclei. The detailed results are listed in Table 1. For comparison, the Gamow-like model [21], generalized liquid drop model [24], Sreeja formula [29], and Liu formula [30] are also employed. In this table, the first three columns represent the$ 2p $ radioactive parent nucleus,$ 2p $ radioactivity released energy, and the logarithmic form of experimental$ 2p $ radioactive half-lives (Expt), respectively. The last five columns represent the logarithmic form of theoretical$ 2p $ radioactive half-lives, which are calculated in this study and with the Gamow-like model (GLM) [21], generalized liquid drop model (GLDM) [24], Sreeja formula (Sreeja) [29], and Liu formula (Liu) [30], respectively. In order to intutively provide comparisons of the experimental$ 2p $ radioactive half-lives with the calculated values, we present the decimal logarithm deviations between the experimental$ 2p $ radioactive half-lives and the calculated values in Fig. 4. From this figure, it can be seen that the decimal logarithm deviations between the experimental$ 2p $ radioactive half-lives and the values calculated in this study are within$ \pm 1 $ . However, for$ ^{54} {\rm{Zn}}$ with$ Q_{2p} $ = 1.280 MeV, the decimal logarithm deviations between the experimental$ 2p $ radioactive half-lives and the values calculated using the theoretical approaches of the Gamow-like model, generalized liquid drop model, Sreeja formula, Liu formula, and this study approach approximately an order of magnitude. The result suggests that this experimental data may not be accurate enough for theoretical comparisons. Furthermore, in order to compare the reproducibility of our model with other theoretical models and/or formulas for$ 2p $ radioactive half-lives, we calculate the standard deviation σ using Eq. (17), with the results listed in Table 2. From this table, we find that the minimum σ obtained by our study is 0.736. Compared to the best result of previous studies, the value of σ is reduced by$ \dfrac{0.846-0.736}{0.846} = 13.0 $ %. It means that our study can be treated as a new and effective tool to investigate$ 2p $ radioactivity.Nuclei $Q_{2p}$ /MeV$\log_{10}T_{1/2}$ /sExpt This study GLM GLDM Sreeja Liu 6Be 1.371 [8] −20.30 [8] −19.86 −19.70 −19.37 −21.95 −23.81 12O 1.638 [61] >−20.20 [61] −17.70 −18.04 −19.71 −18.47 −20.17 1.820 [3] −20.94 [3] −18.03 −18.30 −19.46 −18.79 −20.52 1.790 [62] −20.10 [62] −17.98 −18.26 −19.43 −18.74 −20.46 1.800 [63] −20.12 [63] −18.00 −18.27 −19.44 −18.76 −20.48 16Ne 1.330 [3] −20.64 [3] −15.47 −16.23 −16.45 −15.94 −17.53 1.400 [64] −20.38 [64] −15.71 −16.43 −16.63 −16.16 −17.77 19Mg 0.750 [12] −11.40 [12] −10.58 −11.46 −11.79 −10.66 −12.03 45Fe 1.100 [11] −2.40 [11] −2.32 −2.09 −2.23 −1.25 −2.21 1.140 [10] −2.07 [10] −2.67 −2.58 −2.71 −1.66 −2.64 1.154 [15] −2.55 [15] −2.78 −2.74 −2.87 −1.80 −2.79 1.210 [65] −2.42 [65] −3.24 −3.37 −3.50 −2.34 −3.35 48Ni 1.290 [66] −2.52 [66] −2.55 −2.59 −2.62 −1.61 −2.59 1.350 [15] −2.08 [15] −3.00 −3.21 −3.24 −2.13 −3.13 54Zn 1.280 [67] −2.76 [67] −1.31 −0.93 −0.87 −0.10 −1.01 1.480 [13] −2.43 [13] −2.81 −3.01 −2.95 −1.83 −2.81 67Kr 1.690 [14] −1.70 [14] −0.95 −0.76 −1.25 0.31 −0.58 Table 1. Comparisons between the experimental
$2p$ radioactive half-lives and the calculated values using five different theoretical models and/or formulas. The experimental$2p$ radioactive half-lives$\log_{10}T_{1/2}^{\rm{expt}}$ and corresponding$2p$ released energies$Q_{2p}$ are extracted from the different references.Figure 4. (color online) Deviations between the experimental
$2p$ radioactive half-lives and the calculated values using different theoretical models and/or formulas for the true and not true$2p$ radioactive nuclei.Model This study Gamow GLDM Sreeja Liu σ 0.736 0.846 0.852 1.221 0.851 Table 2. Standard deviation σ between the experimental data and the calculated data using different theoretical models and/or formulas for true
$ 2p$ radioactivity.Finally, we extend this model to predict the half-lives of possible
$ 2p $ radioactive nuclei with$ Q_{2p} > 0$ , which are extracted from the evaluated atomic mass table AME2016 [68, 69]. The predicted results with other predictions using the Gamow-like model, GLDM, and two empirical formulas of Sreeja and Liu are all listed in Table 3. This table is similar to Table 1, except the experimental$ 2p $ radioactive half-lives denoted as "Expt" are replaced by the orbital angular momentum denoted by l in the third column. Based on Table 3, we plot the logarithmic form of the predicted$ 2p $ radioactive half-lives for these possible$ 2p $ radioactive nuclei in Fig. 5. From this figure, we observe that our predicted results are in agreement with the results from the Gamow-like model, generalized liquid drop model, Sreeja formulas, and Liu formulas. In addition, we plot the predicted$ 2p $ radioactive half-lives$ \log_{10}T_{1/2}^{\rm{Pre}} $ against$ ( Z_d^{0.8} + l^{0.25} ) Q_{2p}^{-1/2} $ , i.e., the new Geiger-Nuttall law for$ 2p $ radioactivity, [30] in Fig. 6. The figure depicts an approximately straight line, which indicates that our predicted results are reliable.Nuclei $ Q_{2p}$ /MeVl $ \log_{10}T_{1/2}$ /sThis study GLM GLDM Sreeja Liu 22Si 1.283 0 −12.17 −13.25 −13.30 −12.30 −13.74 26S 1.755 0 −12.82 −13.92 −14.59 −12.71 −14.16 34Ca 1.474 0 −8.99 −10.10 −10.71 −8.65 −9.93 36Sc 1.993 0 −10.79 −12.00 − −10.30 −11.66 38Ti 2.743 0 −12.70 −13.84 −14.27 −11.93 −13.35 39Ti 0.758 0 −1.91 −0.91 −1.34 −0.28 −1.19 40V 1.842 0 −8.97 −10.15 − −8.46 −9.73 42Cr 1.002 0 −2.87 −2.65 −2.88 −1.78 −2.76 47Co 1.042 0 −1.13 −0.42 − 0.21 −0.69 56Ga 2.443 0 −7.41 −8.57 − −6.42 −7.61 58Ge 3.732 0 −11.10 −12.32 −13.10 −9.53 −10.85 59Ge 2.102 0 −5.41 −6.31 −6.97 −4.44 −5.54 61As 2.282 0 −5.78 −6.76 − −4.74 −5.58 10N 1.300 1 −16.76 −17.36 − −20.04 −18.59 28Cl 1.965 2 −11.78 −13.11 − −14.52 −12.46 32K 2.077 2 −11.13 −12.49 − −13.46 −11.55 57Ga 2.047 2 −4.94 −5.91 − −5.22 −4.14 60As 3.492 4 −7.88 −9.40 − −10.84 −8.33 Table 3. Comparison of the predicted half-lives for possible
$ 2p$ radioactive nuclei whose$ 2p$ radioactivity is energetically allowed or observed but not yet quantified in NUBASE2016 [70].Figure 5. (color online) Predicted
$2p$ radioactive half-lives using different theoretical models and/or formulas for possible$2p$ radioactive nuclei whose$2p$ radioactivity is energetically allowed or observed but not yet quantified in NUBASE2016 [70].Figure 6. Predicted
$2p$ radioactive half-lives$\log_{10}T_{1/2}^{\rm{Pre}}$ plotted against$(Z_d^{0.8} + l^{0.25})Q_{2p}^{-1/2}$ , i.e., the new Geiger-Nuttall law for$2p$ radioactivity [30].
Systematic study of two-proton radioactivity with a screened electrostatic barrier
- Received Date: 2021-05-31
- Available Online: 2021-10-15
Abstract: In this study, a phenomenological model is proposed based on Wentzel-Kramers-Brillouin (WKB) theory and applied to investigate the two-proton (