Loading [MathJax]/jax/element/mml/optable/MathOperators.js

Exploring the uncertainties in theoretical predictions of nuclear β-decay half-lives

Cited by

1. Ma, S.-B., Sheng, L.-N., Zhang, X.-H. et al. Opportunities for production and property research of neutron-rich nuclei around N = 126 at HIAF[J]. Nuclear Science and Techniques, 2024, 35(6): 97. doi: 10.1007/s41365-024-01454-w
2. Gangopadhyay, G., Lahiri, C. Recent advances in some nuclear properties relevant to the astrophysical r-process[J]. European Physical Journal: Special Topics, 2024. doi: 10.1140/epjs/s11734-024-01215-1
3. Xia, J.-G., Li, W.-F., Fang, J.-Y. et al. An empirical formula of nuclear β-decay half-lives | [原子核 β 衰变寿命经验公式][J]. Wuli Xuebao/Acta Physica Sinica, 2024, 73(6): 062301. doi: 10.7498/aps.73.20231653
4. Li, W.F., Zhang, X.Y., Niu, Y.F. et al. Comparative study of neural network and model averaging methods in nuclear β-decay half-life predictions[J]. Journal of Physics G: Nuclear and Particle Physics, 2024, 51(1): 015103. doi: 10.1088/1361-6471/ad0314
5. Li, W., Zhang, X., Niu, Z. Studies of nuclear β-decay half-lives with Bayesian neural network approach | [贝叶斯神经网络方法对原子核 β 衰变寿命的研究][J]. He Jishu/Nuclear Techniques, 2023, 46(8): 080013. doi: 10.11889/j.0253-3219.2023.hjs.46.080013
6. Chen, J., Fang, J.Y., Hao, Y.W. et al. Impact of Nuclear β-decay Half-life Uncertainties on the r-process Simulations[J]. Astrophysical Journal, 2023, 943(2): 102. doi: 10.3847/1538-4357/acaeab
7. Fang, J.Y., Chen, J., Niu, Z.M. Gross theory of β decay by considering the spin-orbit splitting from relativistic Hartree-Bogoliubov theory[J]. Physical Review C, 2022, 106(5): 054318. doi: 10.1103/PhysRevC.106.054318
8. Ma, J.-Q., Zhang, Z.-H. Improved phenomenological nuclear charge radius formulae with kernel ridge regression[J]. Chinese Physics C, 2022, 46(7): 074105. doi: 10.1088/1674-1137/ac6154
9. Ma, C.-W., Wei, X.-B., Chen, X.-X. et al. Precise machine learning models for fragment production in projectile fragmentation reactions using Bayesian neural networks[J]. Chinese Physics C, 2022, 46(7): 074104. doi: 10.1088/1674-1137/ac5efb

Figures(6) / Tables(1)

Get Citation
M. Shi, J. Y. Fang and Z. M. Niu. Exploring the uncertainties in theoretical predictions of nuclear β-decay half-lives[J]. Chinese Physics C. doi: 10.1088/1674-1137/abdf42
M. Shi, J. Y. Fang and Z. M. Niu. Exploring the uncertainties in theoretical predictions of nuclear β-decay half-lives[J]. Chinese Physics C.  doi: 10.1088/1674-1137/abdf42 shu
Milestone
Received: 2020-10-14
Article Metric

Article Views(1613)
PDF Downloads(53)
Cited by(9)
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:

Exploring the uncertainties in theoretical predictions of nuclear β-decay half-lives

  • 1. School of Mathematics and Physics, Anhui Jianzhu University, Hefei 230601, China
  • 2. Key Laboratory of Architectural Acoustic Environment of Anhui Higher Education Institutes, Hefei 230601, China
  • 3. Key Laboratory of Advanced Electronic Materials and Devices, Anhui Jianzhu University, Hefei, 230601, China
  • 4. School of Physics and Materials Science, Anhui University, Hefei 230601, China
  • 5. Institute of Physical Science and Information Technology, Anhui University, Hefei 230601, China

Abstract: Nuclear β-decay half-lives are predicted based on an empirical formula and the mass predictions from various nuclear models. It is found that the empirical formula can reproduce the nuclear β-decay half-lives well, especially for short-lived nuclei with T1/2<1 s. The theoretical half-life uncertainties from β-decay energies and the parameters of the empirical formula are further investigated. It is found that the uncertainties of the half-lives are relatively large for heavy nuclei and nuclei near the neutron-drip line. For nuclei on the r-process path, the uncertainties for those with N=126 are about one order of magnitude, which are much larger than the uncertainties for those with N=50 and 82. However, theoretical uncertainties from the parameters of the empirical formula are relatively small for the nuclei on the r-process path, which indicates that the empirical formula is very suitable for predicting the β-decay half-lives in r-process simulations.

    HTML

    I.   INTRODUCTION
    • The origin of elements in the universe is a hot topic not only in astrophysics but also in nuclear physics [1, 2]. The rapid neutron-capture process (r-process) is responsible for the production of about half of the elements heavier than iron, synthesizing elements by continuous neutron captures and β decays [3-5]. The β decay governs the r-process matter flow between neighboring isotopic chains, and hence sets the time scale of the r-process. Therefore, β decay is a key process in r-process nucleosynthesis, and has been the focus of much attention in recent r-process studies [6-9].

      There are some potential candidates for r-process sites, which can provide the conditions of high neutron density and high temperature, such as the neutrino-driven wind after the formation of a proto-neutron star in core collapse supernovae [10-12] and neutron star mergers [13, 14]. However, the r-process sites have still not been fully determined, partly due to the uncertainties in its nuclear physics inputs, including nuclear β-decay half-lives and neutron-capture rates [15]. Therefore, it is essential to obtain accurate nuclear β-decay half-lives from theoretical and experimental studies. Although the measurements of β-decay half-lives have made great progress in recent years [16, 17], many neutron-rich nuclei related to the r-process still cannot be reached by present or even next-generation experimental facilities, so theoretical predictions are essential to r-process simulations. Many approaches have been developed to study nuclear β-decay half-lives, including the gross theory [18-20] and the shell model [21]. Moreover, quasiparticle random-phase approximation (QRPA) approaches based on the finite-range droplet model (FRDM) [22-24], density functional theory (DFT) [25-27], and covariant DFT [28-32] have also been widely employed to predict nuclear β-decay half-lives. In addition, the QRPA with particle-vibration coupling has also been applied to the predictions of nuclear β-decay half-lives [33, 34]. However, there are large variations among these half-life predictions with different approaches. Therefore, it is necessary to investigate the uncertainties in theoretical predictions of nuclear β-decay half-lives.

      The shell model is only feasible for light nuclei or nuclei near the magic numbers, due to the limitation of configuration spaces. Although the QRPA approach can in principle be applied to arbitrary heavy systems, it is still very time consuming when making systematic calculations for deformed systems. Recently, an empirical formula based on the Fermi theory has been proposed to predict the nuclear β-decay half-lives [35]. The parameters in this empirical formula can be obtained by fitting to the known half-lives, and it can even deal with both shell and pairing effects, so it can describe nuclear β-decay half-lives with a high accuracy approaching that of machine learning calculations [36]. This gives us the confidence to predict the half-lives of unknown neutron-rich nuclei with this empirical formula. The β-decay energy is a key input for this empirical formula, in addition to the proton number Z and neutron number N. However, the decay energies predicted by different nuclear mass models are very different, which can certainly induce the variations in the half-life predictions.

      In this work, we will investigate the uncertainties of nuclear β-decay half-life predictions with this empirical formula by employing ten nuclear mass models. These include two microscopic mass models (the relativistic mean field (RMF) [37] and the Hartree-Fock-Bogoliubov (HFB-31) [38]), five macroscopic-microscopic mass models (the extended Thomas-Fermi plus Strutinsky integral (ETFSI-Q) [39], ETFSI-2 [40], the Koura-Tachibana-Uno-Yamada (KTUY) [41], FRDM12 [42] and the Weizsäcker-Skyrme (WS4) [43]), one macroscopic mass model (the Bethe and von Weizsäcker (BW2) [44]), and two other global mass models (Bhagwat [45] and DZ28 [46]). These nuclear mass models can provide reasonable estimates of nuclear masses, while the accuracy of their mass predictions can be further improved by the radial basis function (RBF) approach [47-49], the RBF approach with odd-even effects [50-52], and the Bayesian neural network (BNN) approach [53-56]. Therefore, it is interesting to explore the uncertainties of half-life predictions with this empirical formula. Special attention is focused on the theoretical uncertainties from nuclear β-decay energies and the parameters of the empirical formula. The corresponding results and discussion are presented in Section III. The basic formula employed in this work will be given in Section II, and a summary will be presented in Section IV.

    II.   FORMALISM
    • The nuclear β-decay half-life can be described with the well-known Fermi theory [57], where the phase-space factor is complicated and should be calculated with a numerical method. The emitted electron from β decay of a neutron-rich nucleus generally has a very large speed, and hence the phase-space factor in the Fermi theory can be approximated with the extreme-relativistic limit [35]. With the corrections from pairing and shell effects, the nuclear β-decay half-life T1/2 can be predicted with the following empirical formula [35]

      logT1/2=a6+(α2Z25a7NZA)log(Qa8δ)+a9α2Z2+13α2Z2logAαZπ+S(Z,N),

      (1)

      where α=1/137 is the fine structure constant, A=Z+N is the mass number, and δ=(1)N+(1)Z is a quantity related to the nuclear pairing effect. The last term S(Z,N) is introduced to describe the shell effect in nuclear β-decay half-lives, which is:

      S(Z,N)=a1e((N28)2+(Z20)2)/12+a2e((N50)2+(Z38)2)/43+a3e((N82)2+(Z50)2)/13+a4e((N82)2+(Z58)2)/24+a5e((N110)2+(Z70)2)/244.

      (2)

      Q in Eq. (1) is the β-decay energy, which can be determined from mass predictions of nuclear models. The nine parameters ai(i=1,2,...,9) in Eqs. (1) and (2) are optimized from a least squares fit to the experimental half-lives with the Levenberg-Marquardt method. The experimental half-lives are taken from NUBASE2016 [16], and only those nuclei with Z,N8 and T1/2106 s are retained. To study the influence of these parameters on half-life predictions, they are determined with three different data sets, i.e. those sets for nuclei with half-lives shorter than 106, 103, and 1 s, which are denoted by TS61/2,TS31/2, and TS01/2 for simplicity hereafter. There are 1009, 824, and 381 experimental half-lives in the data sets for TS61/2,TS31/2, and TS01/2, respectively.

      In order to evaluate the performance of this empirical formula with Q values for different mass models, the root-mean-square (rms) deviation of the logarithms of the calculated half-lives Tcal1/2 with respect to that of the experimental data Texp1/2 is employed in the following

      σrms(log10T1/2)=1nni=1[log10(Tcal1/2Texp1/2)]2i.

      (3)

      Here, n is the number of nuclei in a given data set from NUBASE2016 [16]. For each nucleus, the theoretical uncertainties of this empirical formula from different Q predictions are estimated with the following rms quantity:

      δrms(log10T1/2)=1m1mj=1[log10(Tcal1/2Taver1/2)]2j,

      (4)

      where Taver1/2 is the average value of Tcal1/2 from the different mass models. m is the number of nuclear mass models, which are used to determine Q values. To get a reasonable estimation of uncertainty, we only calculate δrms(log10T1/2) for those nuclei whose Q values can be obtained by at least 6 mass models.

    III.   RESULTS AND DISCUSSION
    • We study the accuracy of the empirical formula for different nuclei with different half-lives. The rms deviations σrms(log10T1/2) defined in Eq. (3) are shown in Table 1 for various nuclear mass models. It is clear that this formula can reproduce the known half-lives T1/2<106 s within one order of magnitude for all mass models considered here. The BW2 and WS4 models give the largest and smallest σrms(log10T1/2), which are 0.94 and 0.66, respectively. This accuracy is similar to, and even better than, that obtained by the QRPA approach based on FRDM [23], whose σrms(log10T1/2)=0.82. It is also found that σrms(log10T1/2) gradually decreases for the nuclear sets with shorter half-lives, which indicates that this formula better reproduces the experimental data of nuclei with shorter half-lives. The best accuracy is achieved with the WS4 mass prediction, whose σrms(log10T1/2)=0.22 for the nuclear set with T1/2<1 s, i.e. it can reproduce the experimental half-lives within about a factor of 2.

      Model σrms(log10TS61/2) σrms(log10TS31/2) σrms(log10TS01/2)
      Bhagwat 0.69 0.41 0.24
      BW2 0.94 0.58 0.30
      DZ28 0.74 0.42 0.23
      ETFSI-2 0.73 0.45 0.23
      ETFSI-Q 0.79 0.57 0.26
      FRDM2012 0.72 0.44 0.25
      HFB-31 0.81 0.51 0.27
      KTUY 0.67 0.43 0.26
      RMF 0.83 0.55 0.30
      WS4 0.66 0.39 0.22

      Table 1.  The rms deviations σrms(log10T1/2) of nuclear β-decay half-lives with respect to the experimental data for various nuclear mass models. The rms deviations for TS6, TS3, and TS0 are calculated with respect to the experimental data of nuclear sets with T1/2<106 s, T1/2<103 s, and T1/2<1 s, respectively.

      To check the extrapolation ability of the empirical formula, σrms(log10T1/2) of TS01/2 for the nuclear set with T1/2<106 s are shown by the right-hand bars in Fig. 1. For comparison, the corresponding results for TS61/2 are shown by the left-hand bars in Fig. 1. Since the known half-lives with 1s<T1/2<106s are not the fitting data of TS01/2, it is found that the accuracies of TS01/2 are worse than those of TS61/2 for the nuclear set with T1/2<106 s. However, they are still reasonable, with σrms(log10T1/2) between 0.81 and 1.08 for various nuclear mass models. This indicates that the empirical formula has reliable extrapolation ability even for long-lived nuclei, whose half-lives are generally difficult to describe with microscopic nuclear models.

      Figure 1.  (color online) The rms deviations σrms(log10T1/2) for the nuclear set with T1/2<106 s. For each nuclear model, the left- and right-hand bars correspond to the results for TS61/2 and TS01/2, respectively.

      As in Fig. 1, we show σrms(log10T1/2) of TS61/2 for the nuclear set with T1/2<1 s in Fig. 2, and the corresponding results for TS01/2 are also shown for comparison. Clearly, TS61/2 still reproduces the experimental data better for the short-lived nuclei, whose σrms(log10T1/2) for the nuclear set with T1/2<1 s is only about half of its value for nuclear set with T1/2<106 s. The accuracy of TS61/2 is slightly larger than that of TS01/2, but it is still reasonable to make half-life predictions. For example, the accuracy of TS61/2 for the WS4 model is similar to those of the TS01/2 for BW2 and RMF models. Therefore, both TS61/2 and TS01/2 will be used to evaluate the uncertainties of half-lives for those unknown neutron-rich nuclei in the following, and their differences are employed to describe the uncertainties of half-life predictions from the parameters of the empirical formula.

      Figure 2.  (color online) Same as Fig. 1, but the rms deviations σrms(log10T1/2) are for the nuclear set with T1/2<1 s.

      To investigate the predictive ability of the empirical formula in detail, we show the predicted β-decay half-lives of the N=50,82, and 126 isotones with TS61/2 and TS01/2 in Fig. 3. Clearly, the empirical formula reproduces the experimental data well for both TS61/2 and TS01/2. TS61/2 predicts longer half-lives than TS01/2 for nuclei with T1/2 s, and shorter half-lives than T_{1/2}^{\rm S0} for nuclei with T_{1/2} \lesssim 0.01 s, while they predict similar half-lives for nuclei with T_{1/2}\sim 0.01 s. Furthermore, it is found that the uncertainties in the half-lives are relatively large for N = 50 and 82 isotones with T_{1/2} \lesssim 0.01 s and T_{1/2} \gtrsim 1 s, while the uncertainties remain large for all N = 126 isotones shown here. The large uncertainties can be understood from the sensitivity of T_{1/2} to the Q value and the large uncertainties of the Q values. For the long-lived nuclei, their Q values are so small that tiny differences in Q values would induce large uncertainties in the half-life predictions. For the short-lived nuclei, the uncertainties in the Q values become large since they are far from the known region, which leads to the large uncertainties in the half-lives. For N = 126 isotones, only a few long-lived nuclei are known, so the uncertainties in the Q values are always large for the unknown nuclei and hence induce the large uncertainties in the half-lives. In addition, it should be pointed out that the uncertainties for the r-path nuclei at N = 126 are about one order of magnitude, which is relatively larger than those for r-path nuclei at N = 50 and 82 . Therefore, future experiments on r-process-path nuclei around N = 126 would help to reduce the uncertainties in their half-lives and further improve our understanding of the r-process. Similar to Fig. 3, the half-lives of Z = 28,\; 50 , and 82 isotopes are shown in Fig. 4. The conclusion found in Fig. 3 is also valid for the Z = 28,\; 50 , and 82 isotopes, which further testify to the reliability of the empirical formula to study nuclear \beta -decay half-lives.

      Figure 3.  (color online) The logarithm of nuclear \beta -decay half-lives for N = 50, \;82,\;126 isotones. The red-hatched and olive-hatched bands correspond to the error bands of T_{1/2}^{\rm S0} and T_{1/2}^{\rm S6} , respectively, which are evaluated with \sigma _{\rm{rms}}(\log_{10} T_{1/2}) in Eq. (3). For comparison, the experimental half-lives are shown by the filled circles.

      Figure 4.  (color online) Same as Fig. 3, but for Z = 28,\; 50 , and 82 isotopes.

      In order to systematically study the theoretical uncertainties in the calculation of nuclear \beta -decay half-lives with the empirical formula, the rms deviations \delta_{\rm{rms}}(\log_{10} T_{1/2}) defined in Eq. (4) for T_{1/2}^{\rm S0} and T_{1/2}^{\rm S6} are shown in Fig. 5(a) and (b), respectively. Clearly, the uncertainties in the \beta -decay half-lives from nuclear mass uncertainties for T_{1/2}^{\rm S0} are generally smaller than those for T_{1/2}^{\rm S6} . The heavy nuclei and the nuclei near the neutron-drip line are generally far away from the known region, so different models usually predict very different nuclear masses and hence lead to large uncertainties in the \beta -decay energies. From Eq. (1), it is clear that the nuclear \beta -decay half-lives are very sensitive to decay energies, so the uncertainties in the half-lives would be relatively larger for the heavy nuclei and the nuclei near the neutron-drip line, which is also clearly shown in Fig. 5. For the nuclei around Z = 70 and N = 160 , their \beta -decay energies suddenly drop from about 20 MeV to about 5 MeV for the ETFSI-Q model. However, they remain about 20 MeV for all other models considered here. This sudden drop of decay energy leads to the sudden increase of half-lives in the ETFSI-Q model. Therefore, the predicted half-lives from the ETFSI-Q model are very different from those of other models for these nuclei, which further leads to large uncertainty around Z = 70 and N = 160 in both subfigures of Fig. 5.

      Figure 5.  (color online) Uncertainties of nuclear \beta -decay half-lives \delta_{\rm{rms}}(\log_{10} T_{1/2}) for T_{1/2}^{\rm S0} and T_{1/2}^{\rm S6} . To guide the eye, the r-process path calculated with the RMF mass model [37] and the stable nuclei are shown by the black lines and black filled squares, respectively.

      The differences in the average values of \beta -decay half-lives between T_{1/2}^{\rm S0} and T_{1/2}^{\rm S6} are shown in Fig. 6. A remarkable odd-even staggering is found for the differences in \beta -decay half-lives, and the differences for the even-even nuclei are generally larger than the odd-odd nuclei. This implies the importance of pairing correlation on the accurate predictions of \beta -decay half-lives. It is interesting to find that the differences in \beta -decay half-lives are generally smaller for nuclei around the r-process path, so it would be relatively reliable to make r-process calculations using the \beta -decay half-life predictions of this empirical formula. The data table for the \beta -decay half-life predicted in this work is available at xsx.ahu.edu.cn/T12Data.

      Figure 6.  (color online) Same as Fig. 5, but for the differences in the average values of \beta -decay half-lives between T_{1/2}^{\rm S0} and T_{1/2}^{\rm S6} .

    IV.   SUMMARY
    • In this work, the uncertainties in theoretical predictions of nuclear \beta -decay half-lives have been explored with the latest empirical formula. Special attention is focused on the theoretical uncertainties from nuclear \beta -decay energies and the parameters of the empirical formula. It is found that the uncertainties in the half-lives are relatively large for heavy nuclei and nuclei near the neutron-drip line. For the nuclei on the r-process path, the uncertainties for those with N = 126 are about one order of magnitude, which is considerably larger than those with N = 50 and 82 . The parameters in the empirical formula were determined for three different data sets with half-lives shorter than 10^{6} , 10^{3} , and 1 s, with the resulting empirical formulas denoted T_{1/2}^{\rm S6} , T_{1/2}^{\rm S3} , and T_{1/2}^{\rm S0} , respectively. The differences between the predictions of T_{1/2}^{\rm S6} and T_{1/2}^{\rm S0} were employed to describe the uncertainties in the half-life predictions from the parameters of the empirical formula. The empirical formula reproduces the experimental data well for both T_{1/2}^{\rm S6} and T_{1/2}^{\rm S0} , especially for those short-lived nuclei with T_{1/2}< 1 s. T_{1/2}^{\rm S6} predicts longer half-lives than T_{1/2}^{\rm S0} for nuclei with T_{1/2} \gtrsim 0.01 s and shorter half-lives than T_{1/2}^{\rm S0} for nuclei with T_{1/2} \lesssim 0.01 s, while they predict similar half-lives for nuclei with T_{1/2}\sim 0.01 s. Since the half-lives of nuclei on the r-process path are typically around 0.01 s, the uncertainties from the parameters of the empirical formula are relatively small for these r-path nuclei. This indicates that the empirical formula is very suitable for \beta -decay half-life predictions in r-process simulations.

    ACKNOWLEDGMENTS
    • We would like to express our gratitude to N. N. Zhou, B. R. Yu, and Z. B. Pei for fruitful discussions.

Reference (57)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return