Systematic study on α-decay half-lives of uranium isotopes with a screened electrostatic barrier

Figures(4) / Tables(2)

Get Citation
Yang-Yang Xu, De-Xing Zhu, You-Tian Zou, Xi-Jun Wu, Biao He and Xiao-Hua Li. Systematic study on α decay half-lives of uranium isotopes with a screened electrostatic barrier[J]. Chinese Physics C. doi: 10.1088/1674-1137/ac7fe8
Yang-Yang Xu, De-Xing Zhu, You-Tian Zou, Xi-Jun Wu, Biao He and Xiao-Hua Li. Systematic study on α decay half-lives of uranium isotopes with a screened electrostatic barrier[J]. Chinese Physics C.  doi: 10.1088/1674-1137/ac7fe8 shu
Milestone
Received: 2022-05-24
Article Metric

Article Views(1446)
PDF Downloads(71)
Cited by(0)
Policy on re-use
To reuse of subscription content published by CPC, the users need to request permission from CPC, unless the content was published under an Open Access license which automatically permits that type of reuse.
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article

Title:
Email:

Systematic study on α-decay half-lives of uranium isotopes with a screened electrostatic barrier

  • 1. School of Nuclear Science and Technology, University of South China, Hengyang 421001, China
  • 2. School of Math and Physics, University of South China, Hengyang 421001, China
  • 3. College of Physics and Electronics, Central South University, Changsha 410083, China
  • 4. National Exemplary Base for International Sci & Tech. Collaboration of Nuclear Energy and Nuclear Safety, University of South China, Hengyang 421001, China
  • 5. Cooperative Innovation Center for Nuclear Fuel Cycle Technology & Equipment, University of South China, Hengyang 421001, China
  • 6. Key Laboratory of Low Dimensional Quantum Structures and Quantum Control, Hunan Normal University, Changsha 410081, China

Abstract: In the present work, we systematically study the α-decay half-lives of uranium (Z=92) isotopes based on the Gamow model with a screened electrostatic barrier. There are only two adjustable parameters in our model i.e. the parameter g and the screening parameter t in the Hulthen potential for considering the screened electrostatic effect of the Coulomb potential. The calculated results are in good agreement with experimental data, and the corresponding root-mean-square (rms) deviations of uranium isotopes with α transition orbital angular momentum l=0 and l=2 are 0.141 and 0.340, respectively. Moreover, we extend this model to predict α-decay half-lives of uranium isotopes whose α decay is energetically allowed or observed but not yet quantified in NUBASE2020. For comparison, the modified Hatsukawa formula (XLZ), the unified Royer formula (DZR), the universal decay law (UDL) and the Viola–Seaborg–Sobiczewski formula (VSS) are also used. The predictions are basically consistent with each other. Meanwhile, the results also indicate that N=126 shell closure is still robust at Z=92 and the spectroscopic factor Sα is almost the same for uranium isotopes with the same l.

    HTML

    I.   INTRODUCTION
    • α decay has always been a hot topic in nuclear physics because it can provide abundant nuclear structure information such as ground state [1], nuclear shell effect [2], energy levels [3], nuclear shape coexistence [46], low lying states [7] and so on. As one of the most important decay modes of superheavy nuclei [8], α decay was first explained by Rutherford in 1908 in terms of a process where a parent nucleus emits a 4He particle [9]. Following the foundation and development of quantum mechanics, Gamow [10] and Condon and Gurney [11] in 1928 intepreted the theory of α decay as a quantum penetration of α particles via tunneling.

      Up to now, based on Gamow theory, a great deal of models and/or approaches have been put forward to study α decay. The commonly used models are the cluster model [1218], the Coulomb and proximity potential model (CPPM) [19, 20], the two-potential approach (TPA) [2127], the density dependent M3Y (DDM3Y) effective interaction [2833], the shell model [34], the generalized liquid drop model (GLDM) [3537], the fission-like model [3840] and so on. These models, with their own merits and failures, have been in acceptable agreement with experimental data. Moreover, there are lots of efficient and useful empirical formulas to calculate the α-decay half-lives such as the Viola–Seaborg formula (VS) [41], the universal decay law (UDL) [42, 43], the Royer formula [44], the Hatsukawa formula [45] and Qian et al. formula (YQZR) [46, 47] while their improvements are the Viola–Seaborg–Sobiczewski formula (VSS) [48] and the modified Viola–Seaborg formula (MVS) [49], the modified universal decay law (MUDL) [50], the unitary Royer formula (DZR) [51], the modified Hatsukawa formula (XLZ) [52] and the modified Yibin et al. formula (MYQZR) [53], respectively. From the experimental aspect, the development of radioactive beams and low temperature detector technology make it possible to synthesize new superheavy nuclei [54] and search new α-emitters in artificially occurring nuclides [55, 56].

      Recently, we were intrigued by a series of reports on the synthesis of uranium isotopes. In 2015, the new neutron-deficient isotope 215U was produced in the complete-fusion reaction 180W(40Ar,5n)215U [57]. Evaporation residues recoiled from the target were separated in flight from the primary beam by the gas-filled recoil separator SHANS [58] and subsequently identified on the basis of energy–position–time correlation measurement. The α-particle energy and half-life of 215U were determined to be 8.428(30) MeV and 0.73+1.330.29 ms, respectively. In 2016, Zhang et al. identified two α-decaying states in 216U, one for the ground state and the other for the isomeric state with 8+(πh9/2π7/2) configuration. The α-decay properties for 215,216U and the systematics of 8+ isomeric state in N = 124, 126 isotones were also investigated [59, 60]. In 2021, a new α-emitting isotope 214U, produced by the fusion–evaporation reaction 182W(36Ar,4n)214U was identified by employing the gas-filled recoil separator SHANS and the recoil-α correlation technique. More precise α-decay properties of even–even nuclei 216,218U were also measured in the reactions of 40Ar and 40Ca beams with 180,182,184W targets [61].

      In addition, taking into account the electrostatic screening effect caused by the superposition of the involved charges, Budaca et al. recently proposed a simple analytical model based on the Wentzel–Kramers–Brillouin (WKB) approximation to systematically study the half-lives of proton radioactivity [62]. The results show that with the increase of proton number Z, the difference between the outer turning point radii increases, corresponding to whether the electrostatic screening effect is included in the Coulomb barrier or not. Because the penetration probability is very sensitive to outer turning point radii, it is necessary to introduce the screened electrostatic barrier, i.e. Hulthen potential [63], in the process of dealing with the decay problem. In recent years, the Hulthen potential has been widely used to study the half-lives of α decay, proton radioactivity and two-proton radioactivity [6466], and the calculated results reproduce the experimental data well.

      In 2005, Tavares et al. first evaluated the α-decay half-lives of bismuth isotopes with angular momentum l=5 carried away by the emitted α particle, based on the Gamow model [67]. Since then, it has been successfully generalized to calculate the α-decay half-lives of platinum isotopes with l=0 [68] and neptunium isotopes with l=0 and l=1 [69]. The calculations of these works are in good agreement with experimental data. Based on this model, considering the screened electrostatic barrier, in this work we systematically calculate the α-decay half-lives of uranium isotopes and explore the robustness of N = 126 shell closure at Z = 92 with the experimental data taken from the latest evaluated nuclear properties table NUBASE2020 [70].

      This article is organized as follows. In Section II, the theoretical frameworks of the Gamow model and the screened electrostatic barrier are described in detail. The calculations and discussion are presented in Section III. Finally, Section IV gives a brief summary.

    II.   THEORETICAL FRAMEWORK
    • The α-decay half-life, an important indicator of nuclear stability, is calculated by

      T1/2=ln2λ,

      (1)

      where the decay constant λ can be obtained by

      λ=ν0SαPse.

      (2)

      Here the frequency factor ν0, which represents the number of assualts on the barrier per unit of time [71, 72], is usually estimated as

      ν0=v2a=12a2Qαμ,

      (3)

      where v is the velocity of the α particle inside the parent nucleus. Sα is the spectroscopic factor (also known as α particle preformation probability at the nuclear surface) and Pse is the penetrability factor through the external barrier region [73]. They can be expressed as

      Sα=eGov,Gov=2ba2μ[V(r)Qα]dr,Pse=eGse,Gse=2cb2μ[V(r)Qα]dr,

      (4)

      where is the reduced Planck constant, μ denotes the reduced mass of α particle and daughter nucleus in the center-of-mass coordinate with md and mα being the mass of the daughter nucleus and α particle, and V(r) is the total interaction potential between the emitted α particle and daughter nucleus which is depicted in Fig. 1.

      Figure 1.  (color online) Schematic diagram of one-dimensional α-daughter nucleus interaction potential potential V(r). a, b and c are the three turning points of the potential energy while the barrier from a to b is the overlapping region and b to c is the separation region.

      In general, V(r) is composed of the nuclear potential Vn(r), the Coulomb potential Vc(r) and centrifugal potential Vl(r) by default. It can be expressed as

      V(r)=Vn(r)+Vc(r)+Vl(r).

      (5)

      Here, the Coulomb potential Vc(r) is written as

      Vc(r)=ZαZde2r,

      (6)

      where Zα and Zd represent the proton numbers of α particle and daughter nucleus, respectively, and e2 is the square of the electronic elementary charge. In the process of α decay, for the superposition of the involved charges, movement of the emitted particle which generates a magnetic field and the inhomogeneous charge distribution of the nucleus, the emitted α-daughter nucleus electrostatic potential behaves as a Coulomb potential at short distance and drops exponentially at large distance, i.e. the screened electrostatic effect [6366]. This behaviour of electrostatic potential can be described as the Hulthen type potential which is defined as

      Vh(r)=tZαZde2etr1,

      (7)

      where t is the screening parameter. The centrifugal potential Vl(r) can be expressed as

      Vl(r)=l(l+1)22μr2.

      (8)

      where l is the orbital angular momentum taken away by the emitted α particle. l=0 for the favored α decays, while l0 for the unfavored α decays. On the basis of conservation laws of angular momentum and parity, the minimum angular momentum lmin taken away by the α particle can be obtained by

      lmin={Δj,foreven  Δj  and  πp=πd,Δj+1,foreven  Δj  and  πpπd,Δj,forodd  Δj  and  πpπd,Δj+1,forodd  Δj  and  πp=πd,

      (9)

      where Δj=|jpjd| with jp, πp, πd, jd being the spin and parity values of parent and daughter nuclei, respectively.

      In addition, Gov and Gse are the Gamow factors obtained by the integral from a to b in the overlapping barrier region and from b to c in the separation barrier region in Fig. 1. Here b=Rd+Rα is the separating radius, with Rd and Rα being the radii of the daughter nucleus and the emitted α particle respectively. a=RpRα and c are the inner and outer classical turning points of potential barrier, with Rp being the radius of the parent nucleus. The turning points satisfy the conditions V(a)=V(c)=Qα. The α-decay energy Qα is calculated by

      Qα=ΔMp(ΔMd+ΔMα)+106k(ZpβZdβ) MeV,

      (10)

      where ΔMp, ΔMd and ΔMα represent the mass excess of parent nuclues, daughter nucleus and α particle, respectively, and the quantity kZβ is the total binding energy of Z electrons in the atom, where k = 8.7 eV, β=2.517 for Z60 and k=13.6 eV, β=2.408 for Z<60 [74]. Based on the above, the α-decay half-life can be expressed as

      T1/2=1.0×1022a(μQα)1/2Sα1Pse1.

      (11)

      In the overlapping region, since the α particle is still in the parent nucleus, the reduced mass μ can not be simply treated as a two-body problem. Duarte and Poenaru et al. have successfully dealed with this problem in α particle emitter 222Rn by using μ(r)=[mαmd/(mα+ md)][(ra)/(ba)]3 and V(r)= Qα+(V(b) Qα)[(ra)/(ba)]2 [38, 75]. Encouraged by these two descriptions, the power functions of both μ(r) and V(r) in the overlapping region can likewise be expressed by the following forms [67]

      μ(r)=(mαmdmα+md)(raba)p,p0 ,

      (12)

      V(r)=Qα+(V(b)Qα)(raba)q,q1,

      (13)

      with

      V(b)=Vh(b)+Vl(b)=tZαZde2etb1+l(l+1)22μb2.

      (14)

      Using Eqs. (4), (12) and (13), Gov can be obtained as

      Gov=0.4374702(ba)(1+p+q2)1{μ[2tZde2etb1+20.9008l(l+1)μb2Qα]}1/2,

      (15)

      where (1+p+q2)1 is defined as g with 0g23.

      In the separation region, the parent nucleus has separated into two free individuals, the daughter nucleus and the α particle. We can deal with the reduced mass μ(r) as the reduced mass of the final decay system μ=mαmd/(mα+md). Meanwhile, the potential energy V(r) including the Hulthen type potential and centrifugal potential can be calculated by V(r)=Vh(r)+Vl(r). From the above, Gse can be obtained as

      Gse=1.25988794Zd(μQα)1/2×F,

      (16)

      where

      F=x1/22y×ln{[x(x+2y1)]1/2+x+y(x/y)[1+(1+x/y2)1/2]1+y}+arccos{12[111/y(1+x/y2)1/2]}1/2[12y(1+x/2y1/2y)]1/2,

      (17)

      with

      x=20.9008l(l+1)μb2Qα ,y=ln(tZαZde2/Qα+1)2tb.

      (18)
    III.   RESULTS AND DISCUSSION
    • Based on the Gamow model, considering the screened electrostatic effect and introducing the Hulthen type potential, we propose an improved model to evaluate the α-decay half-lives of uranium isotopes. According to Eq. (9), we select as our database 14 nuclei with experimental α-decay half-lives of the ground-state to ground-state α transition with orbital angular momentum l=0, namely 216,218,221,222,224,226,228,229,230,232,233,234,236,238U, excluding the newly synthesized nucleus 214U, and 5 cases with l=2, namely 217,223,225,227,231U. Based on this database, using a genetic algorithm with an optimal solution of σ i.e. the deviation between the experimental data and calculated values as the objective function, we can obtain the values of the adjustable parameters t and g. In this work, σ is defined as follows

      σ=(log10Tcalc1/2log10Texp1/2)2/n,

      (19)

      where log10Texp1/2 and log10Tcalc1/2 are the logarithmic form of experimental and calculated α-decay half-lives respectively, and n is the number of nuclei involved for each case. The detailed results of the corresponding σ, t and g for the 14 α transitions with l=0 and the 5 with l=2 are listed in Table 1.

      σ/t/gl
      l = 0l = 2
      cases145
      σ0.1410.340
      t5.5239×1043.0519×108
      g0.06930.1549

      Table 1.  Standard deviations between the experimental α-decay half-lives and calculated values using our improved Gamow model for l=0 and l=2 uranium isotopes, along with the corresponding parameters t and g.

      Using our improved Gamow model and the obtained values of parameters t and g, we calculate the α-decay half-lives for these uranium isotopes, 14 with l=0 and 5 with l=2. The detailed calculations are listed in Table 2. In this table, the first five columns denote the α transition, decay energy Qα, spin–parity transformation, the minimum angular momentum taken away by the α particle, and the spectroscopic factor, respectively. The sixth column is the logarithm of the experimental α-decay half-life. The last five columns represent α-decay half-lives in logarithmic form of the uranium isotopes, calculated by our improved Gamow model, the modified Hatsukawa formula (XLZ) [52], the unified Royer formula (DZR) [51], the universal decay law (UDL) [42, 43] and the Viola–Seaborg–Sobiczewski formula (VSS) [48], respectively. As can be seen from this table, the calculated results of our model are close to the experimental values for most of the nuclei.

      α transitionQαjπpjdπlSαlgTexp1/2lgTcal1/2lgTXLZ1/2lgTDZR1/2lgTUDL1/2lgTVSS1/2
      Part I: α tansitions with l=0
      216U212Th8.370+0+00.422.652.422.352.212.182.57
      218U214Th 8.610+0+00.423.193.143.072.962.913.26
      221U217Th 9.939/2+#9/2+#00.436.186.496.196.126.365.70
      222U218Th 9.520+0+00.435.335.565.475.505.405.63
      224U220Th 8.670+0+00.423.403.393.333.243.193.42
      226U222Th 7.740+0+00.410.570.580.560.360.360.61
      228U224Th 6.840+0+00.412.752.722.722.992.942.66
      229U225Th 6.523/2+3/2+00.404.244.074.264.714.284.89
      230U226Th 6.030+0+00.406.246.316.306.616.506.19
      232U228Th 5.450+0+00.399.349.399.369.689.519.19
      233U229Th 4.955/2+5/2+00.3912.7012.5312.7013.1712.5913.12
      234U230Th 4.900+0+00.3912.8912.8612.8313.1312.9112.56
      236U232Th 4.610+0+00.3914.8714.8714.8415.1214.8714.51
      238U234Th 4.310+0+00.3917.1517.2417.2117.4617.1716.80
      Part II: α tansition with l=2
      217U213Th 8.471/2#5/220.121.711.812.431.842.471.95
      223U219Th 9.217/2+#9/2+#20.124.193.934.534.044.633.96
      225U221Th 8.055/2+7/2+#20.111.210.651.330.691.340.69
      227U223Th 7.28(3/2+)(5/2)+20.111.821.981.251.961.271.90
      231U227Th 5.625/2+#(1/2+)20.109.959.498.669.458.629.19
      Part III: other α emitters in uranium isotopes
      219U215Th 9.999/2+#(1/2)54.226.294.686.475.84
      235U231Th 4.727/25/2+116.3514.2814.8414.1314.66
      Part IV: Predictions
      214U210Th 8.530+0+00.423.282.852.772.642.613.04
      215U211Th 8.635/2#5/2#00.423.132.852.582.902.41
      220U216Th 10.330+0+00.437.347.237.377.246.58
      237U233Th 4.271/2+1/2+00.3917.5617.7318.1517.4817.98
      240U236Th 4.080+0+00.3919.2519.2119.4419.1219.64
      242U238Th 3.710+0+00.3822.8622.8122.9722.5923.08
      239U235Th 4.175/2+1/2+#20.1019.5418.5919.3018.3218.84
      241U237Th 3.867/2+#5/2+#20.1022.4221.4422.1121.0821.58
      243U239Th 3.609/2#7/2+#124.1824.5923.7224.20

      Table 2.  α-decay half-lives in logarithmic form of uranium isotopes calculated by our improved Gamow model, the modified Hatsukawa formula (XLZ), the unified Royer formula (DZR), the universal decay law (UDL) and the Viola–Seaborg–Sobiczewski formula (VSS), which are denoted as lgTcal1/2, lgTXLZ1/2, lgTDZR1/2, lgTUDL1/2 and lgTVSS1/2, respectively. The experimental α-decay half-lives, spin and parity are taken from the latest evaluated nuclear properties table NUBASE2020 [70] except for 216U and 218U, whose accurate experimental α-decay half-lives are taken from [Phys. Rev. Lett. 126, 152502 (2021)]. The Qα values are taken from the latest evaluated atomic mass table AME2020 [76, 77]. The α-decay energies and half-lives are in units of MeV and s, respectively.

      To intuitively display the accuracy of our results, in Fig. 2, we plot in logarithmic form the deviations between the experimental half-lives and the calculated values for α emitters of uranium isotopes with l=0 and l=2 using our improved Gamow model, XLZ, DZR, UDL and VSS. They are denoted as red sphere, gray square, blue pentagon, pink triangle and orange star, respectively. From this figure, we find that the differences for our improved Gamow model are basically concentrated near zero. This indicates that the calculated α-decay half-lives using our model can reproduce the experimental values well. It should be emphasized that the experimental α-decay half-lives are taken from the latest evaluated nuclear properties table NUBASE2020 [70] except for 216U and 218U, because the work [Phys. Rev. Lett. 126, 152502 (2021)] in 2021 provides us with more accurate experimental data.

      Figure 2.  (color online) Deviations in logarithmic form between the experimental and calculated α-decay half-lives for uranium isotopes with l=0 and l=2. The red sphere, gray square, blue pentagon, pink triangle and orange star represent the deviations calculated by our improved Gamow model, XLZ, DZR, UDL and VSS, respectively.

      In the following, we extend the improved Gamow model to predict the α-decay half-lives of uranium isotopes with α transition l=0 and l=2 whose α decay is energetically allowed or observed but not yet quantified in NUBASE2020. The predicted α-decay half-lives are listed in Part IV of Table 2, along with the values calculated by XLZ, DZR, UDL and VSS for comparision.

      The XLZ model was proposed by Xu et al. [52] to calculate the favored α-decay half-lives in 2022. It is expressed as

      log10T1/2=F(Z)×(AdApQα)1/2×[arccosXX(1X)]20.446+C(Z,N)+h,

      (20)

      where the coefficient F(Z)=28.274Z+2920.347/Z204.086, and C(Z,N) represents the effect of shells on α-decay half-lives. This can be written as

      C(Z,N)={1.5470.077(82Z)0.050(126N),78Z82  and  100N<126,1.3970.116(Z82)0.061(126N),82<Z90  and  110N126.

      (21)

      The last term h represents the blocking effect of unpaired nucleons, whose values for different parent nuclei are expressed as

      h={0,forevenevennuclei,hp=hn=0.2018,foroddAnuclei,hp+hn=0.4036,foroddoddnuclei.

      (22)

      The DZR model for α decay was proposed by Deng et al. [51], and can be expressed as

      log10T1/2=a+bA1/6Z+cZQα+dl(l+1)+h,

      (23)

      where A, Z, Qα and l represent the mass number, proton number, α decay energy of parent nuclei and the angular momentum taken away by the emitted α particle, respectively. By fitting the experimental data, the values of adjustable parameters are a = –26.8125, b = –1.1255, c = –1.6057, and d = 0.0513. The values of h for different α decay cases are expressed as

      hlog={0, for even-even nuclei, 0.2812, for odd Z-even N nuclei, 0.3625, for even Z-odd N nuclei, 0.7486, for odd-odd nuclei. 

      (24)

      The UDL model for α decay and cluster radioactivity modes was proposed by Qi et al. [42, 43], and can be expressed as

      log10T1/2=aZcZdAQc+bAZcZd(A1/3d+A1/3c)+c,

      (25)

      where A=AcAd/(Ac+Ad), with Ac and Ad being the mass of cluster and daughter nucleus. The constants a = 0.4314, b = –0.4087 and c = –25.7725 are determined by fitting the experimental data of both α and cluster decays.

      The VSS formula is a five-parameter formula put forward by Viola and Seaborg for calculating α-decay half-lives [41]. It can be expressed as

      log10T1/2=(aZ+B)Q1/2α+cZ+D+hlog,

      (26)

      where Z is the atomic number of the parent nucleus and a, b, c and d are equal to 1.66175, –8.5166, –0.20228 and –33.9069, respectively [48]. The last term hlog is expressed as

      hlog={0, for even-even nuclei, 0.772, for odd Z-even N nuclei, 1.066, for even Z-odd N nuclei, 1.114, for doubly odd nuclei. 

      (27)

      For l=0, there are seven nuclei i.e. 214,215,220,237,240,242U, among which 214U was newly synthesized in 2021 and its measured α-decay half-life value is 0.52+0.950.21 ms [61]. It can be well reproduced by the present model as 1.41 ms, which also checks the reliability of our improved Gamow model. Meanwhile, the predicted α-decay half-lives of 215,220,237,240,242U are 7.42×104, 4.61×108, 3.63×1017, 1.77×1019 and 7.20×1022 s, respectively. For l=2, there are two nuclei, 239,241U, with the predicted α-decay half-lives being 3.44×1019 and 2.64×1022 s. In order to test the reliability of predictions calculated by our model, we plot the relationship between log10T1/2 and Q1/2α, i.e. the Geiger–Nuttall law [78], in Fig. 3. From this figure, we can clearly see that our predicted α-decay half-lives for these eight nuclei 214,215,220,237,239,240,241,242U fit the linear relationship well. This indicates that our predicted results may be useful for future study of α-decay half-lives in newly synthesized uranium isotopes. In addition, there are two other α emitters of uranium isotopes 219U and 235U with l=5 and l=1 which can not be calculated by this work. We list them in Part III of Table 1.

      Figure 3.  (color online) Logarithmic form of α-decay half-life log10T1/2 vs Q1/2α for uranium isotopes. The line represents the Geiger–Nuttall law, black squares and red spheres denote the experimental and calculated half-lives, respectively.

      To verify whether the shell closure N = 126 is robust or not at Z = 92, in Fig. 4 we plot the logarithmic form of α-decay half-lives log10T1/2 and α-decay energy Qα against the neutron number N of the parent nuclei. From this figure, it can be clearly seen that log10T1/2 and Qα present completely opposite trends with the change of N. Particularly, when N is less than 126, log10T1/2 generally shows a slow trend to fluctuate, and drops sharply after 126 until N = 128. A similar phenomenon can be seen on Qα, where the values of Qα increase dramatically to N = 128 when the neutron number is greater than 126. The above phenomena reflect strong shell effects at N = 126.

      Figure 4.  (color online) Relationship between the logarithmic form of α-decay half-lives log10T1/2, α-decay energy Qα, spectroscopic factor Sα and neutron number N of the parent nucleus : (a) for the logarithmic form of α-decay half-lives log10T1/2, the red and blue triangle represent the experimental half-lives Texp1/2 and predicted values Tcal1/2, respectively; (b) for the α-decay energy Qα; (c) for Sα with l=0 and l = 2 in neptunium isotopes, where rhombus and pentagon represent l = 0 and l = 2 while pink and pale yellow denote experimental data and predicted values, respectively.

      Finally, we extract the spectroscopic factor Sα in two cases with l = 0 and l = 2, then plot them versus the neutron number N of the parent nuclei in Fig. 4. From this figure, we find a amazing phenomenon that the Sα values are basically the same for α transitions with the same orbital angular momentum l, with the extracted values Sα0.406 for l = 0 and Sα0.112 for l = 2. It shows that high angular momentum will hinder the spectroscopic factor and α-decay half-life [7981]. Meanwhile, its small variation trend can also reflect the robustness of N = 126 at Z = 92. The above phenomena are worthy of further study in the future.

    IV.   SUMMARY
    • In this work, based on the Gamow model, considering the effect of electrostatic screening, we systematically study the α-decay half-lives of uranium (Z = 92) isotopes. The calculated results are in reasonable agreement with the experimental data. In addition, we extend this model to predict α-decay half-lives of uranium isotopes whose α decay is energetically allowed or observed but not yet quantified in NUBASE2020, along with the predictions of XLZ, DZR, UDL and VSS for comparision. The predicted results of our model and these three formulas are close and consistent in trend. In addition, the robustness of the shell closure N = 126 is also verified at Z = 92. Finally, we find the values of spectroscopic factor Sα are basically the same for α transitions with the same orbital angular momentum l. The results of this work will prompt inquiries about nuclear structure and provide information for future experiments.

    ACKNOWLEDGEMENTS
    • We would like to thank J. G. Deng, H. M. Liu, X. Pan and H. F. Gui for useful discussion.

Reference (81)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return