Released energy formula for proton radioactivity based on the liquid-drop model*

Xiao Pan(潘霄),You-Tian Zou(邹有甜),Hong-Ming Liu(刘宏铭),Biao He(何彪),Xiao-Hua Li(李小华),3,4,5,∗∗and Dong Xiang(向东),3,4,∗∗

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 National Exemplary Base for International Sci & Tech.Collaboration of Nuclear Energy and Nuclear Safety,University of South China,Hengyang 421001,China

4 Cooperative Innovation Center for Nuclear Fuel Cycle Technology & Equipment,University of South China,Hengyang 421001,China

5 Key Laboratory of Low Dimensional Quantum Structures and Quantum Control,Hunan Normal University,Changsha 410081,China

Abstract In this work,based on the liquid-drop model and considering the shell correction,we propose a simple formula to calculate the released energy of proton radioactivity(Qp).The parameters of this formula are obtained by fitting the experimental data of 29 nuclei with proton radioactivity from ground state.The standard deviation between the theoretical values and experimental ones is only 0.157 MeV.In addition,we extend this formula to calculate 51 proton radioactivity candidates in region 51≤Z≤83 taken from the latest evaluated atomic mass table AME2016 and compared with the Qp calculated by WS4 and HFB-29.The calculated results indicate that the evaluation ability of this formula for Qp is inferior to WS4 while better than HFB-29.

Keywords:the released energy of proton radioactivity,liquid-drop model,shell correction

1.Introduction

With the continuous developments of radioactive beam facilities,the study of exotic nuclei far from the β-stability line has became a hot project in nuclear physics field[1–9].Exotic nuclei usually exhibit distinctive features that can provide a lot of information about the nuclear structure different from the well-known stable nuclei,such as new shell closures[10],coupling of bound state and unbound state and so on[11].Proton radioactivity,the spontaneous emission of a proton by the nucleus,is an important decay mode of exotic nuclei.In 1970,this decay mode was firstly observed in an isomeric state of53Co by Iackson et al[12,13].Subsequently,proton emissions from ground state of151Lu[14]and147Tm[15]were detected by Hofmann et al and Klepper et al.Up to now,there are 44 proton emitters decaying from their ground states or isomeric states between Z=51 and Z=83[11,16–22]being identified in experiments.

Proton radioactivity shares the similar theory of barrier penetration with different kinds of charged particles’radioactivity,such as α decay,heavy ion emission,spontaneous fission,etc[23–31].Theoretically,a great deal of methods have been proposed to deal with proton radioactivity half-life,which can be divided into two categories.One kind is theoretical models that the probability of proton penetration barrier is calculated by the Wentzel–Kramers–Brillouin approximation,such as the unified fission model[32,33],the single-folding model[17],the modified two-potential approach[20,34],the Gamow-like model[35],the effective interactions of density-dependent M3Y[36–38],the generalized liquid-drop model[19,28],the Coulomb and proximity potential model,etc[21,39,40].The other one is empirical formulas,which rely on the accurate released energy Qpand orbital angular momentum l taken away by emitted proton,such as universal decay law for proton radioactivity[41],the formula of Zhang and Dong[1],the New Geiger–Nuttall law[42],etc.Using these formulas,ones can estimate the half-life of proton radioactivity easily and rapidly.However,without precise Qp,we can not to predict the half-life of proton radioactivity for the certain nuclei whose proton radioactivity is energetically allowed or observed but not yet quantified.Therefore,obtaining a simple and accurate Qpformula is very necessary.In 2010,based on the macroscopic-microscopic model,a simple and local formula for calculating α-decay energies of superheavy nuclei was proposed by Dong et al[43].They used this formula to calculate the α-decay energy,the theoretical calculations can well reproduce the experimental data.Meanwhile,they also predicted the α-decay energy of newly synthesized superheavy nuclei and obtained positive results.Refering to the shell correction form obtained by Dong et al,we propose the Qpformula including the liquid-drop part and shell correction.The liquid-drop part is the main part of our formula,which is derived from the Bethe–Weizsäcker binding energy formula.

This article is organized as follows.In the next section,the formula of released energy for proton radioactivity is briefly described.The detailed calculations and discussion are presented in section 3.Finally,a summary is given in section 4.

2.The formula of released energy for proton radioactivity

The relationship between the released energy of proton radioactivity and binding energy can be expressed as Qp=B(A−1,Z−1)−B(A,Z).If the change of binding energy with Z and A is relatively smooth,it can be replaced by

withΔA=ΔZ=−1.From the Bethe–Weizsäcker binding energy formula

we obtain the following formula

HereΔBpis the change of pairing energy terms when the parent nuclei emits a proton.Up to now,the observed nuclei with proton radioactivity are all odd-Z nuclei.ThenΔBpcan be divided into the following categories

So,we can obtained the liquid-drop part of Qpformula as follows

However,simply considering the part of the liquid-drop in this formula,the Qpcan only reflect the average trend of binding energy.While the effect of shell correction on the Qpis very important to nuclei around shell closures.Refering to the shell correction form for α-decay energies proposed by Dong et al,the final Qpformula is expressed as

where N0is the number of neutron for the shell correction,and it will be shown in the results and discussion.

3.Results and discussion

The σ represents the deviation between the calculated value and experimental ones,in this work,which is defined as follows:

where(expt.)and(cal.)denote the experimental data and calculated ones of the ith proton emitter released energy,respectively.

Table 1.The experimental data[42]and calculated results for 29 nuclei with proton radioactivity from ground state obtained by equation(6).

At first,fitting the experimental data of Qpfor 29 nuclei with proton radioactivity from ground state listed in table 1,we obtain the best parameters for equation(5)as follows:

Using equation(5)and the above parameters,we obtain the σ=0.179 MeV calculated by equation(7).For a more intuitive displaying the deviation and attempting to find N0,we plot the difference of Qpbetween the calculated value and experimental ones in figure 1(a).From this figure we can clearly see there is a peak around N=76,meaning that N0=76 required by equation(6).

In the following,fixed the parameters a1to a5as mentioned above and let N0in the equation(6)equal to 76,fitting the experimental data of the same 29 nuclei,we obtain the a6=0.4745 and a7=−0.0686.Using these parameters and equation(6),we calculate the Qpfor these 29 nuclei.The results are listed in table 1.In this table,the first three columns represent the proton emitter,the experimental Qpand the calculated result,respecttively.The last column is the difference of Qpbetween the calculated value and experimental data.In order to intuitively survey their deviations,we plot the difference of Qpbetween the calculated value and experimental ones in figure 1(b).From table 1,it is clearly seen that the deviations of Qpare almost less than 0.2 MeV except 7 nuclei.The σ is 0.157 MeV calculated by equation(6).Comparing to equation(5),the result is slightly improved by 12.3% with considering the shell correction.From figure 1(b),we can see that the peak around N=76 is reduced.It may be the contribution by shell correction.

Finally,as an application,we extend this formula to calculate 51 proton radioactivity candidates in region 51≤Z≤83 taken from the latest evaluated atomic mass table AME2016[44,45].For comparison,the Qpare also evaluated using Weizsäcker–Skyrme(WS4)formula[46]and Hartree–Fock–Bogoliubov(HFB-29)model[47]In these two models,the Qpis calculated by[21]

Figure 1.The part(a)is the difference of proton radioactivity released energy between the calculated value and experimental one[42]using equation(5)without the shell correction,and the part(b)shows the results obtained using equation(6)with the shell correction term.

whereΔMp,ΔMdandΔMprare the mass excess of parent nuclei,daughter nuclei and emitted proton,respectively.The last itemrepresents the screening effect of the atmoic electrous[48],where k=13.6eV,ε=2.408 for Z<60 and k=8.7 eV,ε=2.517 for Z≥60[49].

Table 2.The experimental data and calculated Qp values of this work,the WS4 model and the HFB-29 model.Δ1,Δ2 andΔ3 are the deviation between experimental data and calculated value for equation(6),the WS4 model and the HFB-29 model.

Figure 2.The difference of Qp between the experimental data and calculated value for this work,WS4 and HBF-29.

The detailed results are given in table 2.In this table,the first five columns represent the proton emitter,the experimental data and the calculated Qpobtained by this work,WS4 and HFB-29,respectively.The last three column are the difference of Qpbetween the experimental data and calculated value derived by the three models mentioned above,denoted asΔ1,Δ2,Δ3.From this table,it indicates that all of the Qpin this work are positive,while in WS4 or HFB-29 the Qpof some nuclei such as119La,129Pm,145Ho,154,155Lu,171Ir,177Au are negative.To obtain further insight into the well of agreement and the systematics of results,we plot the difference of Qpbetween the experimental data and calculated value using these three models in figure 2.From this figure,we can clearly see that the best result is WS4 to reproduce the experimental data but the worst is HFB-29 in these three models.The evaluation abilityour of our work is between them.Furthermore,we calculate the standard deviation between experimental Qpand calculated ones for this work,WS4 and HFB-29 are 0.347,0.231 and 0.482 MeV,respectively.The reason may be that this work and WS4 are considered the effect of shell correction,but WS4 considers more comprehensively.Due to the uncertainty of skyme force and the pairing force that its corresponding parameters are obtained by fitting the properties of less nuclei and then extend it to the unknown regions,the results of Qpcalculated by the HFB-29 model are inferior to the first two models.

4.Summary

In summary,based on the liquid-drop model and considered shell correction,we put forward a simple formula for the Qp.The corresponding parameters are obtained by fitting the experimental data of 29 nuclei with proton radioactivity from ground state.The standard deviation between the theoretical values and experimental ones is 0.157 MeV.In order to further test the reliability of this formula,we calculate the Qpvalues of 51 proton radioactivity candidates in region 51≤Z≤83 obtained from the latest evaluated atomic mass table AME2016 and compare with the Qpvalues extracted from the WS4 and HFB-29 models.It is found that the accuracy of our formula is higher than that of the HFB-29 model,although its accuracy is lower than that of the WS4 model.This indicates our formula is suitable for the predictions of Qpto some extent.

ORCID iDs