Systematic study on the proton radioactivity of spherical proton emitters

Dong-Meng Zhang ·Lin-Jing Qi ·De-Xing Zhu ·Yang-Yang Xu ·Song Luo ·Xiao-Hua Li,2,3,4

Abstract In this study, based on a two-potential approach, we systematically investigated the proton radioactivity half-lives of spherical proton emitters with 69 ≤Z ≤81 from the ground and/or isomeric state, choosing the nuclear potential to be a modified Woods-Saxon potential that contains the isospin effect of the daughter nucleus.It was found that the calculated half-lives could reproduce the experimental data well.Furthermore, we extended this model to predict the half-lives of 17 protonemitting candidates whose radioactivity is energetically allowed or observed but not yet quantified in NUBASE2020.For comparison, the unified fission model, Coulomb potential and proximity potential model, universal decay law for proton emission, and new Geiger-Nuttall law were also used.All the predicted results are consistent with each other.

Keywords Proton radioactivity ·Spherical nuclei ·Two-potential approach ·Modified Woods-Saxon potential ·Half-lives

1 Introduction

Proton radioactivity, the spontaneous process in which a nucleus with a large excess of protons transforms into the stable state by emitting an individual proton, was first observed by Jackson et al.[1, 2] from an isomeric state of53Co.In the early 1980s, Hofmann et al.[3] and Klepper et al.[4] further discovered proton emission from the nuclear ground states of151Lu and147Tm, respectively.With the advancement of diverse infrastructures and radioactive beam installations, more than 40 proton emitters have been provenly illustrated in the proton number range 51 ≤Z≤83 from the ground state or low-lying isomeric state during the last decades [5-10].Proton radioactivity can be used as an effective probe to identify nuclei near the proton drip line and explore diverse nuclear information, such as the released energy, half-lives, branching ratio, and wave function of the parent nucleus [11-14].Therefore, the study of proton radioactivty of proton-rich nuclei far from theβ-stability line has become an extremely hot topic in nuclear physics.

Proton radioactivity shares the same decay mechanism asαdecay [15-20], two-proton radioactivity [21-29], and cluster radioactivity [30, 31], that is, barrier penetration.These processes can be dealt with by using the Wentzel-Kramers-Brillouin (WKB) approximation.To date, many models and/or formulae have been proposed to calculate the half-life of proton radioactivity, which can be divided into the following two categories: The first are theoretical models, including the effective interactions of density-dependent M3Y [32,33], single-folding model [21], unified fission model (UFM)[34, 35], Gamow-like model [36], Coulomb and proximity potential model [37-39], generalized liquid-drop model [7,40-42], and modified two-potential approach [11, 43, 44].The second are empirical formulae proposed in the form of the Geiger-Nuttall (G-N) law [45], such as the universal decay law for proton emission (UDLP) [46], the new Geiger-Nuttall law (NG-N) [47], the four-parameter empirical formula of Zhang and Dong [12], and the screened decay law for proton emission [48].Calculations using these theoretical methods provide excellent estimation of the half-life of proton radioactivity and improve our understanding of the proton radioactivity phenomenon.

For proton radioactivity, it is well known that the emitted particle-nucleus interaction potential is extremely crucial to improving the accuracy of half-life calculations.In 2008,based on the Skyrme-energy-density-function approach, taking into account the isospin effect of the nuclear potential,Wang et al.[49, 50] proposed a modified Woods-Saxon potential to describe a large number of heavy-ion elastic or quasi-elastic scattering and fusion reactions.Later, it was widely used to studyαdecay and cluster radioactivity [51,52].Whether this modified potential can be extended to investigate proton radioactivity or not is an interesting topic.To this end, based on the two-potential approach within the modified Woods-Saxon potential (TPA-MWS), we systematically studied the half-lives of 32 spherical proton emitters.The results indicate that our calculations are consistent with the experimental data.In addition, we extended this model to predict the half-lives of 17 proton-emitting candidates whose radioactivity is energetically allowed or observed but not yet quantified in NUBASE2020.

The remainder of this paper is organized as follows:In Sect.2, the theoretical framework of the two-potential approach and two different empirical formulae are introduced in detail.The results and discussion are presented in Sect.3.Finally, a brief summary is provided in Sect.4.

2 Theoretical framework

2.1 Two-potential approach

The proton radioactivity half-life,T1∕2, which is an important indicator of nuclear stability, is determined by

whereℏis the reduced Planck constant, and Γ represents the proton radioactivity width.In the framework of the twopotential approach [53-55], Γ can be expressed as

whereμ=mpmd∕(mp+md) ≈938.272×Ad∕(Ad+Ap)MeV∕c2represents the reduced mass of the emitted proton and daughter nucleus, withmp,md,Ap, andAdas the masses of the emitted proton and daughter nucleus, and the mass numbers of the emitted proton and daughter nucleus, respectively.Spis the formation probability of the emitted protondaughter system, which is related to the probability that the orbit of the emitted proton is empty in the daughter nucleus[7].In this study, based on Ref.[34], we choseSp=1.Fis the normalized factor describing the collision frequency of the emitted proton in the potential barrier.It satisfies the condition

P, the penetration probability of the emitted proton penetrating the barrier, is calculated by using the semiclassical WKB approximation.It can be expressed as

where the experimental data of the mass excess ΔM,ΔMd,and ΔMpare taken from the latest atomic mass table NUBASE2020 [56], representing the mass excesses of the parent and daughter nuclei and emitted proton, respectively.The termk(Zε-Zεd) denotes the screening effect of atomic electrons, withZandZdbeing the proton number of the parent and daughter nuclei, forZ ＜60 ,k=13.6 eV,ε=2.408 ,and forZ≥60 ,k=8.7 eV,ε=2.517 [57, 58].

In this study, the total interaction potentialV(r) was composed of the nuclear potentialVN(r) , Coulomb potentialVC(r) , and centrifugal potentialVl(r) (see Fig.1).It can be expressed asForVN(r) , we chose the modified Woods-Saxon potential[49, 50], which can be parameterized as

Fig.1 Variation in the total V(r) and nuclear VN(r) potentials for the parent nucleus 170 Au as a function of the nuclear radius.(Color figure online)

Here,R0is the effective nuclear radius, which can be expressed as

whereId=(Nd-Zd)∕AdandIp=(Np-Zp)∕Apare the isospin asymmetries of the daughter nucleus and emitted proton withZd,NdandZp,Npbeing the proton and neutron number of the daughter nucleus and emitted proton,respectively.

VC(r) is the Coulomb potential, which is taken as the potential of a uniformly charged sphere.It can be expressed as

wheree2≈1.44 MeV·fm denotes the square of the electronic elementary charge [61].Rdenotes the sharp radius.For convenience, it can be taken as the effective nuclear radiusR0of the modified Woods-Saxon potential [62].

For the centrifugal potentialVl(r) , we chose the Langer modified form becausel(l+1)→(l+1∕2)2is an essential correction for one-dimensional problems [63].It can be expressed as

Here,lis the angular momentum removed by the emitted proton, which satisfies the spin-parity conservation laws.It can be expressed as

where Δj=|j-jd-jp| , withj,π,jd,πd,jp,πpdenoting the spin and parity values of the parent and daughter nuclei and emitted proton, respectively.

2.2 Two different empirical formulae

2.2.1 Universal decay law for proton emission

Based on the R-matrix theory of the cluster decay process and the microscopic mechanism of charged-particle emission, Qi et al.[64, 65] proposed a universal decay law that improves the accuracy of cluster radioactivity half-lives.Because proton radioactivity shares the same decay mechanism asαdecay and cluster radioactivity, considering the influence of the centrifugal barrier, this formula was later extended to investigate proton radioactivity [46].It is expressed as

2.2.2 New Geiger-Nuttall law

In the form of the G-N law [45], considering the contributions of the orbital angular momentumland the charge of the daughter nucleusZd, Chen et al.[47] proposed a twoparameter empirical formula to study proton radioactivity in a unified way, which is expressed as

where the parametersa=0.843 andb=-27.194 were determined by fitting 44 experimental data points of the proton radioactivity half-lives from the ground state and/or isomeric state.

3 Results anddiscussion

In this study, we systematically investigated the half-lives of 32 spherical proton emitters based on the TPA-MWS.For comparison, the UFM [35], Coulomb potential and proximity potential model with the Guo-2013 formalism(CPPM) [66], UDLP [46] and NG-N [47] were also used.The detailed numerical results are listed in Table 1.In this table, the first three columns contain the proton emitter, proton radioactivity released energyQp, and angular momentumlremoved by the emitted proton.The last six columns contain the logarithmic form of the experimental half-lives and theoretical half-lives calculated using three models and two different empirical formulae denoted as Exp, TPA-MWS,UFM, CPPM, UDLP, and NG-N.As shown in Table 1, the calculated proton radioactivity half-lives obtained using our model can well reproduce the experimental data.To visually present their deviations, we plotted the differences between the experimental proton radioactivity half-lives and the calculated values using these five models and/or formulae in logarithmic form in Fig.2.From this figure, we can see that the deviations are in the range of ±1.Overall, all of the calculated proton radioactivity half-lives were consistent with the experimental data.Nevertheless, it should be noted that the logarithmic deviation between the experimental and calculated half-lives was large, corresponding to the special case177Tlm.Its deviations in logarithmic form were -0.889,-0.802,-1.176,-0.881 , and -0.826 using TPA-MWS,UFM, CPPM, UDLP, and NG-N, respectively.Exploring the causes of such a large deviation is crucial for further investigation of proton radioactivity.As we know, the selection of the proton radioactivity released energyQpand angular momentumlis important in proton radioactivity.In our previous study, we discussed that the effect of angular momentum on the half-life is larger than the released energy [67].The result indicated that the calculated half-life of proton radioactivity for177Tlmwas closer to the experimental data whenl= 6 [67].Therefore, we chosel=6 and calculated its logarithm proton radioactivity half-life as -2.779 with a deviation of 0.567, which has a higher consistency with the experimental data than the calculated value of -4.235 when the value oflis set as 5.In addition, the proton number of177Tlmis close to the magic numberZ= 82, and the proton shell effect may be another reason for this phenomenon.

Fig.2 Deviations between the experimental proton radioactivity half-lives and the calculated values using three models and two different empirical formulae for spherical nuclei.The rhombi, circles,triangles, inverted triangles, and squares correspond to the results obtained using the UFM, CPPM, UDLP, NG-N, and TPA-MWS,respectively.(Color figure online)

To obtain further insight into the agreement between the experimental data and calculated one, the root-mean-square deviationσwas used to estimate the calculated capabilities of the above five models and/or formulae for the proton radioactivity half-lives.It can be defined as

Fig.3 Decimal logarithm of the experimental half-lives plotted against -log10P for proton radioactivity.The slope and intercept of the fitted solid line for each case are shown in the corresponding figure.(Color figure online)

Table 1 Comparison of experimental proton radioactivity half-lives with the calculated values using different theoretical models and/or formulae.The symbol m denotes the first isomeric state.The experimental proton radioactivity half-lives, spin,and parity are taken from Ref.[56].The released energy is given by Eq.5, with the exception of the Qp value for 159Re, 161Rem , 164Ir, and 177Irm ,which are taken from Ref.[6]

where log10and log10Tdenote the experimental proton radioactivity half-life and the calculated one in logarithmic form for thei-th nucleus, respectively.For comparison,the specific results ofσfor the five theoretical models and/or formulae are listed in Table 2.From this table, we can clearly see that the standard deviationσobtained using TPAMWS (σTPA-MWS=0.341 ) was smaller than that of theothers.This means that the present model can better reproduce the experimental proton radioactivity half-lives than the others for spherical emitters.

Table 2 Standard deviations σ between the experimental and calculated half-lives using different theoretical models and/or formulae

A single line of the universal curve forαdecay and cluster radioactivities could be manifested by plotting the decimal logarithm of the half-life versus the negative decimal logarithm of the penetration probability through the barrier [68, 69].Whether this linear correlation can be extended to investigate proton radioactivity is an interesting topic [69].Therefore, we plotted the correlation between the decimal logarithm of the experimental halflife against the negative decimal logarithm of penetrability in Fig 3.The penetrability was calculated theoretically inthe framework of the semiclassical WKB approximation.As shown, the decimal logarithm of the experimental halflife had a clear linear dependence on the negative decimal logarithm of the penetrability, which implies that our calculations and the present methods are valid and reliable.

Table 3 Comparison of the predicted half-lives for possible radioactivity candidates whose proton radioactivity is energetically allowed or observed but not yet quantified in NUBASE2020.The symbols m and n denote the first and second isomeric states,respectively.The released energy is given by Eq.5

Fig.4 Linear relationship between the quantities log10and(+l) based on the NG-N.(Color figure online)

Considering the good agreement between the experimental half-lives and those calculated using the present model, we extended our model to predict the half-lives of proton radioactivity in the proton number region 51 ≤Z≤83.For comparison, the UFM, CPPM, UDLP,and NG-N were also used.The detailed predictions are listed in Table 3.In this table, the first three columns contain the possible proton radioactivity candidates, released energyQp, and orbital angular momentuml.The last five columns contain the predicted proton radioactivity halflives obtained using the five models and/or formulae.From this table, we can clearly see that the magnitude of all results using our model is consistent with those calculated using the other four models and/or formulae, validating the predictive power of this method.Taking116La as an example, the predictions were -10.256, -10.117, -9.963,-9.000, and -9.126.To further test the reliability of the predictions, we obtained a linear relationship between the quantities log10based on the NG-N[47] in Fig.4.The results show that the predictions are consistent with the NG-N, which indicates that our predictions of the proton radioactivity half-lives are valid and reliable.These predictions may provide valuable information for searching for new nuclides with proton radioactivity in future.

Furthermore, recent studies have shown that nuclear deformation plays an important role in proton radioactivity[37, 44, 48, 69].Barrier penetration can be affected by the deformation of daughter nuclei during the proton emission process [70].The deformation effect can result in the reduction of Coulomb and centrifugal barriers, which can influence the penetrability and reduce the theoretical half-lives[11].In addition, some methods can be used to calculate the proton radioactivity half-lives of deformed proton emitters.In 2016, based on the WKB tunneling approximation, Qian and Ren et al.[44] proposed a simple formula including the deformation term to evaluate the proton radioactivity half-lives of the deformed nuclei.In 2022, the deformation dependence of the screened decay law for proton radioactivity was proposed by Budaca et al.[48] to systematically study proton radioactivity.Inspired by this, we will also discuss the influence of the deformation effect on the estimation of proton radioactivity half-lives in more detail in future work.

4 Summary

In summary, we present a systematic study of the proton radioactivity half-lives for 32 spherical proton emitters based on the TPA-MWS, which contains the isospin effect.It was found that the calculated results could reproduce the experimental data well.In addition, we extended this model to predict the proton radioactivity half-lives for 17 possible candidates.The predicted results are in reasonable agreement with those obtained using the UFM, CPPM, UDLP,and NG-N.This study may provide useful and reliable information for predicting the existence of new proton emitters above the proton drip line, which may be detected in future.

Author contributionsAll authors contributed to the study conception and design.Material preparation, data collection and analysis were performed by Lin-Jing Qi, De-Xing Zhu, Yang-Yang Xu, Song Luo and Xiao-Hua Li.The first draft of the manuscript was written by Dong-Meng Zhang, and all authors commented on previous versions of the manuscript.All authors read and approved the final manuscript.Data availability statementThe data that support the findings of this study are openly available in Science Data Bank at https://www.doi.org/10.57760/sciencedb.j00186.00058 and http://resolve.pid21.cn/31253.11.sciencedb.j00186.00058.