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

Research on convergence of the nuclear matrix elements for 2νββ decays

Cited by

1. Gao, Y., Di, Z., Gao, S. General mass formulas for charged Kerr-AdS black holes[J]. Physica Scripta, 2024, 99(9): 095022. doi: 10.1088/1402-4896/ad6fff
2. Sood, A., Ali, M.S., Singh, J.K. et al. Photon orbits and phase transition for Letelier AdS black holes immersed in perfect fluid dark matter[J]. Chinese Physics C, 2024, 48(6): 065109. doi: 10.1088/1674-1137/ad361f
3. Pokhrel, R., Dey, T.K. Charged AdS black holes in presence of string cloud and Cardy-Verlinde formula[J]. Nuclear Physics B, 2024. doi: 10.1016/j.nuclphysb.2024.116508
4. Pokhrel, R., Sherpa, K.P., Dey, T.K. Dissipative Force on an External Quark in AdS Gauss-Bonnet Gravity with String Cloud[J]. Springer Proceedings in Physics, 2024. doi: 10.1007/978-3-031-69146-1_47
5. Pokhrel, R., Sherpa, K.P., Dey, T.K. Holographic study of drag on a probe quark in Reissner-Nordstrom AdS black hole with Gauss-Bonnet gravity and cloud of string[J]. Proceedings of SPIE - The International Society for Optical Engineering, 2024. doi: 10.1117/12.3041616
6. Ndongmo, R., Mahamat, S., Tabi, C.B. et al. Thermodynamics of non-linear magnetic-charged AdS black hole surrounded by quintessence, in the background of perfect fluid dark matter[J]. Physics of the Dark Universe, 2023. doi: 10.1016/j.dark.2023.101299
7. Li, X.-P., Zhang, L.-C., Ma, Y.-B. et al. Thermodynamic quantities and phase transitions of five-dimensional de Sitter hairy spacetime[J]. Chinese Physics C, 2023, 47(10): 105102. doi: 10.1088/1674-1137/ace8f5
8. Alipour, M.R., Sadeghi, J., Shokri, M. WGC and WCCC of black holes with quintessence and cloud strings in RPS space[J]. Nuclear Physics B, 2023. doi: 10.1016/j.nuclphysb.2023.116184
9. López, L.A., Pedraza, O. Effects of quintessence on scattering and absorption sections of black holes[J]. Indian Journal of Physics, 2023, 97(1): 285-294. doi: 10.1007/s12648-022-02373-5
10. Aounallah, H., El Moumni, H., Khalloufi, J. et al. Insight into the microscopic structure of a quintessential black hole from the quantization concept[J]. International Journal of Modern Physics A, 2022, 37(8): 2250036. doi: 10.1142/S0217751X22500361
11. Qu, Y., Tao, J., Wu, J. New Gedanken experiment on RN-AdS black holes surrounded by quintessence[J]. European Physical Journal C, 2022, 82(2): 185. doi: 10.1140/epjc/s10052-022-10120-7
12. Liang, J., Mu, B., Wang, P. Joule-Thomson expansion of lower-dimensional black holes[J]. Physical Review D, 2021, 104(12): 124003. doi: 10.1103/PhysRevD.104.124003
13. Yin, R., Liang, J., Mu, B. Joule–Thomson expansion of Reissner–Nordström-Anti-de Sitter black holes with cloud of strings and quintessence[J]. Physics of the Dark Universe, 2021. doi: 10.1016/j.dark.2021.100884
14. Liang, J., Lin, W., Mu, B. Joule–Thomson expansion of the torus-like black hole[J]. European Physical Journal Plus, 2021, 136(11): 1169. doi: 10.1140/epjp/s13360-021-02119-y
15. Mustafa, G., Hussain, I. Radial and circular motion of photons and test particles in the Schwarzschild black hole with quintessence and string clouds[J]. European Physical Journal C, 2021, 81(5): 419. doi: 10.1140/epjc/s10052-021-09195-5
16. Yin, R., Liang, J., Mu, B. Stability of horizon with pressure and volume of d-dimensional charged AdS black holes with cloud of strings and quintessence[J]. Physics of the Dark Universe, 2021. doi: 10.1016/j.dark.2021.100831

Figures(4) / Tables(4)

Get Citation
Yue Hu, Hantao Li and Chengwei Dong. Research of convergence of the nuclear matrix elements for 2νββ decays[J]. Chinese Physics C. doi: 10.1088/1674-1137/abba13
Yue Hu, Hantao Li and Chengwei Dong. Research of convergence of the nuclear matrix elements for 2νββ decays[J]. Chinese Physics C.  doi: 10.1088/1674-1137/abba13 shu
Milestone
Received: 2020-07-11
Article Metric

Article Views(1521)
PDF Downloads(80)
Cited by(16)
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:

Research on convergence of the nuclear matrix elements for 2νββ decays

Abstract: In this work, the characteristics of 2νββ decays for six nuclei (36Ar, 46Ca, 48Ca, 50Cr, 70Zn, and 136Xe) in a mass range from A=36 to A=136 are studied within the nuclear shell model (NSM) framework. Calculations are presented for the half-lives, nuclear matrix elements (NMEs), phase space factors (G2ν), and convergence of the NMEs. The theoretical results agree well with the experimental data. In addition, we predict the half-lives of 2νββ decays for four nuclei. We focus on the convergence of the NMEs by analyzing the number of contributing intermediate 1+ states (NC) for the nuclei of interest. We assume that NC is safely determined when the accumulated NMEs saturate 99.7% of the final calculated magnitude. From the calculations of the involved nuclei, we discover a connection between NC and the total number of intermediate 1+ states (NT). According to the least squares fit, we conclude that the correlation is NC=(10.8±1.2)×N(0.29±0.02)T.

    HTML

    1.   Introduction
    • Two-neutrino double beta (2νββ) decay is an interesting process (a second-order weak interaction) in the standard model, which was proposed by Mayer [1] in 1935. It is a significant subject in nuclear and particle physics that can be used to test models of the weak interaction [2, 3]. The 2νββ decays have ultra-long half-lives of over 1018 years [4, 5], which is at least one billion times longer than the age of the universe itself. To date, 90 nuclei are believed to undergo 2νββ decay, of which only 14 have been measured [6-10]. The complementary experimental information from related 2νββ decay can improve the quality of nuclear structure models [11]. In addition, research on 2νββ decay can be used to calculate the nuclear matrix elements (NMEs) of 0νββ decay as a promising source of 0νββ decay detection [12]. It is well known that the absolute mass of the neutrino is still a mystery [13, 14]. If 0νββ decay could be detected, it could settle this problem when combined with the NMEs governed by the theory [15, 16]. Thus, whether experimental or theoretical, research into 2νββ decay is extremely significant.

      The role of accurate NMEs in 2νββ decay cannot be underestimated [17]. Although the NMEs between the initial and final nuclei ground states have been investigated using several approaches [18-24], they are still notoriously difficult to calculate, as shown in Eq. (1) [25]:

      M2ν=n0f+iσ(i)τ±(i)1+n1+niσ(i)τ±(i)0i+[12Q+E(1+n)Mi]/me+1.

      (1)

      Here, the transition operators are the ordinary Gamow-Teller operators; E(1+n)Mi denotes the energy discrepancy of the nth intermediate 1+ state and the parent nucleus; and the sum i extends over all the intermediate 1+ states [25], which poses a challenge in calculations for many nuclei.

      In recent years, research has been renewed on the single state dominance (SSD) for 2νββ decay, which suggests that the lowest intermediate 1+ state dominates the decay [26, 27]. To date, evidence of the SSD for 82Se has been established by the CUPID-0 collaboration [28]. One of its proven benefits is that the corresponding NMEs can be obtained through single β decay measurements [29]. However, the SSD of many nuclei violates this benefit, with the result that the relevant NMEs cannot be easily determined experimentally. It would hopefully be simpler to calculate only the lowest 1+ wave function instead of all wave functions in the theoretical description [29]. It is thus of great importance to discuss the number of least intermediate 1+ states that saturate the final magnitude.

      In our work, calculations for six nuclei (36Ar, 46Ca, 48Ca, 50Cr, 70Zn, and 136Xe) in a mass range from A=36 to A=136 are performed within the nuclear shell model (NSM) framework. Calculations are presented for the half-lives, NMEs, G2ν, and convergence of the NMEs. In addition, we predict the half-lives of the 2νββ decays for four nuclei. Relatively short half-lives are predicted for the nuclei 46Ca and 70Zn, providing a reference for future experimental detection. Convergence of the NMEs for 2νββ decays is systematically discussed by analyzing the number of contributing intermediate 1+ states (NC) for the nuclei of interest. A unified criterion for analyzing NC involved in 2νββ decay is proposed and adopted. According to the calculations, the general law for convergence of the NMEs for 2νββ decays is then presented.

    2.   Overview of the theoretical structure

      2.1.   Formalism for 2νββ decays

    • The 2νββ decays can be classified into two types: double β decay (2ν2β) and the decays on the double β+ side, i.e., 2ν2β+, 2νβ+EC, and 2ν2EC (where EC refers to electron capture) [30]. The formulas used for various types of 2νββ decay transitions differ, which can be seen in Refs. [13, 30, 31]. The associated half-lives, in connection with the phase space factors (Gα2ν) and the NME (|Mα2ν|), can be expressed by Eq. (2) [30]:

      1/Tα1/2=Gα2ν|Mα2ν|2, 

      (2)

      where α=2β, 2β+, β+EC, and 2EC are different types of 2νββ decays [13, 30, 31]. The related NME description of the full β strength functions of both the parent and grand nuclei is given [2]. Although formulas for Gα2ν and |Mα2ν| are supplied by the theory, the calculation for |Mα2ν| is more difficult than that for Gα2ν because it requires the precise nuclear many-body wave function [32]. The NME is defined in Eq. (1): accurate knowledge of the ground-state for the parent and grand-daughter nuclei, along with all the 1+ states of the intermediate nucleus, clearly plays an important role in evaluations of the NMEs [16].

      Detailed expressions for various modes of Gα2ν can be found in Ref. [33]. The expression is an integral in the phase space of lepton variables; the corresponding values depend on the coupling constant gA and can be calculated exactly [25, 33]. There is a global effort to address the challenge of the gA factor. This study applies the commonly used value of gA=1.27, which is consistent with that published in Refs. [13, 17, 34].

    • 2.2.   Model spaces, effective interactions, and quenching factors

    • The calculations for six nuclei (36Ar, 46Ca, 48Ca, 50Cr, 70Zn, and 136Xe) are computed in full 1ω shell model spaces without any truncation. Columns 1 and 2 in Table 1 list the nuclei of interest and their shell model spaces, in which the calculations are conducted. To improve the reliability of the calculations, different effective interactions (EIs) are adopted and listed in Column 3. For the 2νββ decay of 36Ar, the specific interactions CWH and W are adopted, where the CWH EI is determined by an iterative least squares fit to experimental level energy data [35]; W EI refers to the two-body matrix elements derived from a linear fit [35, 36]. We use the shell model empirical interaction GXFP1A for the fp shell, which provides 195 two-body matrix elements and 4 single-particle energy subsets through a fit to 699 energy observations within the mass interval A=4766 [37, 38]. Another interaction, KB3, is used to compute the lowest wave function from the full Hamiltonian, with parts of the two-nucleon EI turned off in proper order [39, 40]. It is reasonable to adopt these EIs for the calculations because the concerned nuclei are not far from these regions. We also use the large-scale shell-model interaction SN100PN [41] in the case of 136Xe, which is obtained from the free nucleon-nucleon CD-Bonn interaction. The interaction is composed of four parts: neutron-neutron, neutron-proton, proton-proton, and Coulomb repulsion between the separate protons [42].

      NucleusShell model spaceEINTN.calq
      36Arsd(0d5/2,0d3/2,1s1/2)CWH, W66660.77[35]
      46Cafp(1f7/2,2p3/2,1f5/2,2p1/2)GXFP1A, KB3236123610.74[45]
      48Cafp(1f7/2,2p3/2,1f5/2,2p1/2)GXFP1A, KB3947020000.74[45]
      50Crfp(1f7/2,2p3/2,1f5/2,2p1/2)GXFP1A, KB338393220000.74[45]
      70Znfp(1f7/2,2p3/2,1f5/2,2p1/2)GXFP1A, KB31857120000.68[46]
      136Xegdsh(0g7/2,1d5/2,1d3/2,2s1/2,0h11/2)SN100PN1664220000.45[16]

      Table 1.  Details of the NSM calculations for 2νββ decays.

      The calculations of the wave functions for the involved nuclei use the diagonalization computer program NuShellX [43]. Columns 4 and 5 provide the total number of intermediate 1+ states (NT) and the number of calculated intermediate 1+ states (N.cal), respectively. One of the most prominent advantages of the NSM is that it considers all the many-body correlations of orbits near the Fermi surface [44], but in theory, the expected total intensity is systematically greater than the observed total intensity [16, 35]. One reliable method of overcoming this problem is to correct the calculated strengths with quenching factors (q), which are tabulated in Column 6. In this article, we apply the widely used values q=0.77 for A=1640, q=0.74 for A=4050, and q=0.68 for A=6080 [35, 45, 46]. For the 136Xe nucleus, q=0.45 [16] is adopted.

      This research is inseparable from the precise description of wave functions of intermediate states. To test the reliability of the wave functions of the NSM, we compare the first few theoretical energy levels and the γ-decay half-lives of 36Cl (the intermediate nucleus for the 2νββ decay of 36Ar) in Fig. 1 together with the available experimental information. The theoretical results are obtained with W interaction. As depicted in Fig. 1, the energy levels of the low excited states and most γ-decay half-lives are in good agreement with the experimental data. Thus, within the applicable range, the calculations of the NSM are relatively accurate.

      Figure 1.  Comparison of the theoretical spectrum of 36Cl and the associated γ-decay half-lives with the experimental data. The excitation energies of the states are given in units of MeV and the γ-decay half-lives in picoseconds.

    3.   Results and discussion
    • We calculate the properties of 2νββ decays for six nuclei (36Ar, 46Ca, 48Ca, 50Cr, 70Zn, and 136Xe) with Q values adopted from AME2016 [47]. The calculated half-lives, G2ν, and NMEs are presented in Table 2, together with other results, either from recent experiments [10, 48-51] or theoretical calculations [17, 44, 52, 53]. Most of the calculated results do not differ greatly from the experimental data. For the half-lives of the 36Ar and 50Cr nuclei, our results clearly differ from previous results [52]: the reason is that prevoius works adopted gA=1.

      NucleusQ value /keVG2ν/yr1|M2ν|.1st|M2ν|half-life /yrEI
      this workothers
      36Ar (2ν2EC)432.591.28×10270.002540.1175.68×10281.4×1029[52]CWH
      0.02290.1513.43×1028W
      46Ca (2ν2β)988.41.24×10220.05640.05632.54×10241.7×1024[53]GXFP1A
      0.03050.04683.67×1024KB3
      48Ca (2ν2β)4268.084.17×10170.02190.02503.84×1019(6.4+0.0070.006±+0.0120.009)×1019[10]GXFP1A
      0.01830.02264.68×1019(4.3+2.41.1±1.4)×1019[48]KB3
      50Cr (2ν2EC)1169.61.25×10240.02940.1235.32×1025>1.3×1018[49]GXFP1A
      0.05150.1245.21×10252.1×1026[52]KB3
      70Zn (2ν2β)997.13.24×10220.06770.1381.61×1023[50]GXFP1A
      0.005770.08674.10\times 10^{23}1.4(5)\times 10^{23}[17]KB3
      ^{136}Xe (2\nu2\beta^{-})2457.84.69\times 10^{-18}0.007240.01866.16\times 10^{20}\left( 2.165\pm 0.016\pm 0.059\right) \times 10^{21}[51]SN100PN
      \left( 2.11\pm 0.25\right) \times 10^{21}[44]

      Table 2.  Summary of the calculations for the 2\nu\beta\beta decays of six nuclei.

      Recently, evidence of the SSD for ^{82} Se has been established by the CUPID-0 collaboration [28]. If the SSD is affirmed, then the corresponding NMEs can be obtained through single \beta decay measurements [29]. Therefore, it is of great significance to analyze the validity of the SSD. Calculated results show that for the 2\nu2\beta^{-} decays of ^{46} Ca and ^{48} Ca, the first intermediate 1^{+} state contributes significantly to the final NME. However, there is important canceling in the accumulations of the NMEs, as can be seen in Figs. 7 and 10 in Ref. [52]. In these cases, we cannot conclude that the SSD is satisfied. For other nuclei, the sum of the NMEs gains less than 50% of their final NMEs from the first intermediate 1^{+} state. Therefore, we cannot conclude that the SSD is satisfied for the related nuclei.

      The key of the theoretical research on 2\nu\beta\beta decays is the calculation of the NMEs. Thus, we present the calculated NMEs together with the Q values in Fig. 2. All the data are from Table 2. The figure illustrates the lack of a significant relationship between the NMEs and the Q values. The calculated results from different EIs are consistent with each other, except in the case of ^{70} Zn, for which the NME using GXFP1A EI is 1.6 times that using KB3.

      Figure 2.  (color online) Comparison of the NMEs for different EIs.

      Convergence of the NMEs for 2\nu\beta\beta decays is concluded by analyzing the number of contributing intermediate 1^{+} states ( N_{\rm{C}} ) for the concerned nuclei. We assume that a 2\nu\beta\beta decay NME is well converged when the accumulated NMEs saturate 99.7% of the final magnitude. This is illustrated in Fig. 3, where for the ^{70}{\rm{Zn}}\longrightarrow ^{70}{\rm{Ge}} decay for GXFP1A EI, points A (5th state), B (103rd state), and C (214th state) all satisfy this criterion. Only after point C has satisfied the criterion, indicating that the NME has no distinct twists, is the accumulation believed to be saturated. We thus assume that after point C (214th state), the NME is reliably converged. Accordingly, N_{\rm{C}} is determined naturally. This criterion is used throughout this article.

      Figure 3.  (color online) Accumulations of the NMEs for the ^{70}{\rm{Zn}}\longrightarrow^{70}{\rm{Ge}} decay to ground state.

      A detailed analysis of the convergence of NMEs based on the above criterion is presented in Table 3 for the nuclei of interest. Notice that the actual contributing excitation energy for ^{136} Xe decay is smaller than for the other concerned nuclei (see Column 5), which suggests that the level density is large.

      NucleusN_{\rm{T}}N_{\rm{C}}N_{\rm{C}}/N_{\rm{T}}Ex/MeVEI
      ^{36}Ar663451.5%33.6CWH
      3248.5%36.5W
      ^{46}Ca23611305.51%24.3GXFP1A
      1114.70%22.4KB3
      ^{48}Ca94701501.58%23.4GXFP1A
      1831.93%24.7KB3
      ^{50}Cr3839324170.109%22.5GXFP1A
      4230.110%23.7KB3
      ^{70}Zn185712141.15%23.0GXFP1A
      1650.888%25.4KB3
      ^{136}Xe1664210096.06%12.7SN100PN

      Table 3.  Details of the convergence of NMEs according to the criterion that the accumulated NMEs should saturate 99.7% of the final magnitude.

      From the calculations for the involved nuclei, we discover a connection between the total number of intermediate 1^{+} states ( N_{\rm{T}} ) and the number of contributing states ( N_{\rm{C}} ). A linear analysis for a general convergence law of the NMEs is depicted in Fig. 4. The unfilled circles represent the CWH and GXFP1A EIs, and the solid circles are for the W, KB3, and SN100PN EIs. All the data are from Table 3. We use a solid line for the fit using five nuclei ( ^{36} Ar, ^{46} Ca, ^{48} Ca, ^{50} Cr, and ^{70} Zn) and express the correlation between N_{\rm{C}} and N_{\rm{T}} by Eq. (3):

      Figure 4.  (color online) Correlation between the convergence proportion (N_{\rm{C}}/N_{\rm{T}}) in a \log_{10} frame and N_{\rm{T}} in \log_{10} frame. The solid line is a linear fit including all the nuclei except ^{136}Xe, and the dashed line is for all the nuclei of interest. The unfilled and solid circles represent the results from different EIs.

      N_{\rm{C}} = \left( 10.8\pm 1.2\right) \times N_{\rm{T}}^{\left( 0.29\pm 0.02\right) }, \

      (3)

      with a linear correlation coefficient of 0.996, which is very close to 1. The overall trend of the convergence properties is in good agreement with the correlation, except in the case of the ^{136} Xe nucleus, which deviates slightly from the solid line. The reason may be that the level density of ^{136} Cs (the intermediate nucleus for 2\nu\beta\beta decay of ^{136} Xe) is obviously larger than that of the others. For the other nuclei, the excitation energy exceeds 20 MeV up to the 300th intermediate 1^{+} state, whereas the excitation energy of the 1009th intermediate 1^{+} state for ^{136} Cs is only 12.7 MeV, as listed in Table 3. Thus, the 2\nu\beta\beta decay of ^{136} Xe is unusual. If the ^{136} Xe nucleus is taken into account, the correlation is expressed by Eq. (4),

      N_{\rm{C}} = \left( 10.7\pm 2.1\right) \times N_{\rm{T}}^{\left( 0.32\pm 0.08\right) }, \

      (4)

      with a linear correlation coefficient of 0.884. Hopefully, the number of least intermediate 1^{+} states needed for the precise calculation of NMEs can be obtained from Eq. (3).

      To check the applicability of Eqs. (3) and (4), two new transitions are performed. One is the ^{64}{\rm{Zn}}\longrightarrow ^{64}{\rm{Ni}} decay computed in the fpg model space with JUN45 and JJ44B interactions. Details of the calculated results are presented in Table 4. The calculated values of N_{\rm{C}} are almost identical for the different interactions within the theoretical range of the two equations.

      NucleusShell model spaceN_{\rm{T}}N_{\rm{C}}Theoretical\,{N_{\rm{C}}}EI
      Eq. (3)Eq. (4)
      ^{64}Znfpg12378182167\pm 32219\pm 66JUN45
      fpg12378186167\pm 32219\pm 66JJ44B
      ^{68}Nifp18571205187\pm 34249\pm 73GXFP1A
      fpg42895241238\pm 40326\pm 89JUN45

      Table 4.  Summary of the convergence of the NMEs of new nuclei.

      Another transition of interest is the ^{68}{\rm{Ni}}\longrightarrow ^{68}{\rm{Zn}} decay in the fp and fpg model spaces. Although this decay does not occur in nature, the related calculations could be used to further discuss the applicability of Eqs. (3) and (4) in different model spaces. In addition, the calculated values of N_{\rm{C}} for the different model spaces are within the theoretical range of the two equations. Thus, Eqs. (3) and (4) are feasible within the applicable range. We hope our research is useful for shell model calculations in 2\nu\beta\beta decay.

    4.   Summary
    • In this work, the general features and convergence of the NMEs for the 2\nu\beta\beta decays of six nuclei in a mass range from A = 36 to A = 136 are studied under the framework of the nuclear shell model (NSM). Different effective interactions (EIs) are adopted to improve the reliability of the calculations. Calculations are presented for the half-lives, NMEs, G_{2\nu} , and convergence of the NMEs. Most of the results are very close to the experimental data. In addition, we predict the half-lives of 2\nu\beta\beta decays for four nuclei. The nuclei ^{46} Ca and ^{70} Zn have relatively short half-lives, which will hopefully be probed by future experiments.

      Convergence of the NMEs for 2\nu\beta\beta decays is discussed systematically by analyzing the number of contributing intermediate 1^{+} states ( N_{\rm{C}} ) for the concerned nuclei. We propose and adopt the criterion that the accumulated NMEs should saturate 99.7% of the final calculated magnitude. From the calculations for the nuclei of interest, we discover a connection between N_{\rm{C}} and the total number of intermediate 1^{+} states ( N_{\rm{T}} ). The least squares fit is N_{\rm{C}} = \left( 10.8\pm 1.2\right) \times N_{\rm{T}}^{\left( 0.29\pm 0.02\right) } . It is hoped that the number of least intermediate 1^{+} states needed for the precise calculations of NMEs can be obtained from this equation.

      The authors are grateful to Prof. B. A. Brown of Michigan State University for providing us with the computer program NuShellX.

Reference (53)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return