Next Article in Journal
Oscillation of Emden–Fowler-Type Differential Equations with Non-Canonical Operators and Mixed Neutral Terms
Previous Article in Journal
Controlling I-V Hysteresis in Al/Pt Bilayer Symmetric SQUIDs at Millikelvin Temperatures
Previous Article in Special Issue
Multiple Multi-Orbit Pairing Algebras in Nuclei
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Correlation of Neutrinoless Double-β Decay Nuclear Matrix Element with E2 Strength

1
School of Physics and Astronomy, Sun Yat-sen University, Zhuhai 519082, China
2
Sino-French Institute of Nuclear Engineering and Technology, Sun Yat-sen University, Zhuhai 519082, China
*
Author to whom correspondence should be addressed.
Symmetry 2023, 15(2), 552; https://doi.org/10.3390/sym15020552
Submission received: 20 January 2023 / Revised: 13 February 2023 / Accepted: 16 February 2023 / Published: 19 February 2023

Abstract

:
We explore the correlation of the neutrinoless double- β decay nuclear matrix element (NME) with electric quadrupole (E2) strength in the framework of the Hamiltonian-based generator-coordinate method, which is a configuration-mixing calculation of symmetry-restored intrinsic basis states. The restoration of symmetries that are simultaneously broken in the mean-field states allows us to compute the structural and decay properties associated with wave functions characterized by good quantum numbers. Our calculations show a clear anti-correlation between the neutrinoless double- β decay NME and the transition rate of the collective quadrupole excitation from the ground state in response to artificial changes of the quadrupole–quadrupole interaction. The anti-correlation is more remarkable in the decay from a weakly deformed parent nucleus to a more deformed grand-daughter nucleus. This interrelation may provide a way to reduce the uncertainty of the nuclear matrix element.

1. Introduction

Large-scale experiments in search of neutrinoless double-beta ( 0 ν β β ) decays are likely to determine whether neutrinos are Majorana fermions, and hence reveal corresponding new physics beyond the Standard Model of electroweak interactions. However, the estimation of the rate of 0 ν β β decay, which is crucial for a definitive choice and quantity of candidate nuclei required in these sophisticated experiments, depends on reliable descriptions of nuclear matrix elements (NMEs) [1]. In addition, if this extremely rare decay is observed, the measured decay rates can provide essential information for determining the absolute neutrino mass scale and mass hierarchy, but it relies on whether one can obtain an accurate description of the corresponding NMEs [1].
There are several nuclear structure methods for calculations of the 0 ν β β decay matrix elements. The most used ones are the shell model (SM) [2,3,4,5,6,7,8,9,10], the interacting boson model (IBM) [11,12,13], the quasiparticle random phase approximation (QRPA) [14,15,16,17,18,19,20,21,22,23,24,25], the generator coordinate method (GCM) based on energy density functional (EDF) [26,27,28,29,30] or effective Hamiltonian [31,32,33,34], the in-medium generator coordinate method (IM-GCM) [35], and the coupled-cluster (CC) method [36]. At present, NMEs predicted by these nuclear models differ by factors of 3 to 4 for most candidate nuclides, and up to 8 in the case of decay from 48 Ca to 48 Ti [37,38,39]. These large uncertainties preclude an efficient plan for experiments [37]. Reducing the uncertainty in the matrix elements thus becomes an urgent need in the nuclear structure community.
One way to improve the accuracy of 0 ν β β NME calculations is searching for their correlations with certain nuclear structure observables. Since the observables can be determined experimentally, it could help to pin down the values of 0 ν β β NMEs if there are strong correlations between them. It is found that some certain collective correlations influence significantly on the calculated 0 ν β β NMEs. In particular, the matrix elements would be suppressed when the ground states of the grandparent and grand-daughter nuclei have different intrinsic deformations [40]. This suppression was originally investigated with axial quadrupole collectivity [26,28,29,41], and later extended to non-axial quadrupole [32] and octupole correlations [30]. This leads to an interesting question: whether there is a correlation between 0 ν β β NME and the observable in associate with quadrupole collectivity, e.g., the reduced E 2 transition probability B ( E 2 ) .
Another interesting feature, which implies a potential correlation between 0 ν β β NME and E2 strength, is the Gamow-Teller (GT) transition strength. On one hand, since the GT transition replaces a neutron with a proton or the other way around, its strength of the grandparent to intermediate nuclei, combined with the strength of the grand-daughter to intermediate nuclei, is particularly related to double- β decay. On the other hand, recent work has demonstrated a remarkable anti-correlation between the calculated Gamow-Teller strength and the reduced transition probability of the lowest collective E2 excitation [42]. Though the connection between the GT strength and 0 ν β β NME is not straightforward, it is quite natural to expect a similar correlation of 0 ν β β decay NME with the transition rate of the low-lying collective E2 excitation.
Therefore it would be of particular interest to modify the quadrupole–quadrupole term in the effective Hamiltonian and to examine how the 0 ν β β NME changes are correlated with the resultant E2 strength. It is worth noting that a statistical analysis for the 0 ν β β NME of 48 Ca has been proposed recently by adding random contributions to three sets of p f -shell effective Hamiltonians, i.e., the FPD6, GXPF1A, and KB3G interactions [43]. It gives some hints that the B ( E 2 ) values in 48 Ti are correlated with the neutron occupation probabilities, and hence, indirectly influence the 48 Ca-Ti 0 ν β β NME. It is natural to presume that, instead of random contributions, a detailed analysis of NME by adding a specific term directly related to quadrupole collectivity would be of great importance.
In this work, we propose an analysis of the 0 ν β β NMEs and their correlations with the reduced collective quadrupole transition probabilities B ( E 2 ; 0 + 2 + ) . It is obtained by adding a different amount of quadrupole–quadrupole contributions to the corresponding effective Hamiltonians, and then perform the Hamiltonian-based GCM [32,33] calculation using the modified effective Hamiltonians. We apply our analysis to 48 Ca−Ti, 76 Ge−Se, 124 Sn−Te, 130 Te−Xe, and 136 Xe−Ba, which covers light-, medium-, and heavy-mass candidate nuclei associated with the 0 ν β β decay. We consider only the standard light left-handed Majorana neutrino exchange mass mechanism in this work, since it is the simplest and most studied mechanism of the 0 ν β β decay process.

2. The Model

Owing to the closure approximation, we can express the 0 ν β β decay matrix element in terms of a two-body 0 ν β β transition operator between the ground-state wave functions of the grandparent and grand-daughter nuclei. If the decay is originated from the exchange of light Majorana neutrinos with usual left-handed weak interaction, the NME can be obtained by:
M 0 ν = M GT 0 ν g V 2 g A 2 M F 0 ν + M T 0 ν
where GT, F, and T refer to the Gamow-Teller, Fermi and tensor parts of the matrix elements. g V and g A are the vector and axial coupling constants, which are taken as g V = 1 and g A = 1.254 , respectively. Our wave functions are modified at short distances by using a Jastrow-type short-range correlation (SRC) function in the parametrization of CD-Bonn [19]. A more detailed expression of the matrix element can be found in Ref. [18].
To calculate the 0 ν β β matrix element in Equation (1), a key feature is the description of ground-state wave function of grandparent and grand-daughter nuclei, i.e., | I and | F , which can be given by our GCM calculations. Firstly, we employ a shell-model effective Hamiltonian ( H eff ) . In an isospin scheme, H eff can be written in terms of one- and two-body operators as the following:
H eff = a ϵ a n ^ a + a b , c d J T V J T ( a b ; c d ) T ^ J T ( a b ; c d ) ,
where ϵ a and V J T ( a b ; c d ) stand for the one-body and two-body matrix elements of the nucleon–nucleon interaction, respectively. n ^ a denotes the number operator for the spherical single-particle orbit a labelled with quantum numbers ( n a , l a , j a ) , and
T ^ J T ( a b ; c d ) = M T z A J M T T z ( a b ) A J M T T z ( c d )
is the scalar two-body density operator for nucleon pairs occupied the spherical orbitals a , b and c , d , which are coupled to the quantum numbers J , M , T , and T z .
Changing adequately the effective Hamiltonian used in the calculations, one can increase or decrease the quadrupole correlations in grand-parent, grand-daughter, or both nuclei. In this manner, we can gauge the correlation of the decays with the E2 transition probabilities. We have artificially changed the quadrupole strength adding an extra quadrupole–quadrupole term, namely λ Q ^ · Q ^ , to the effective interaction. Increasing the values of coupling constant λ would increase the deformation and hence would enhance the E2 transition strength.
The next step in our GCM calculation is generating a reference state set | Φ ( q ) . It consists of quasiparticle vacua constrained to different expectation values q i = O i for a set of collective operators O i . In this work, we take the following operators O i :
O 1 = Q 20 , O 2 = Q 22 , O 3 = 1 2 ( P 0 + P 0 ) , O 4 = 1 2 ( S 0 + S 0 ) ,
where
Q 2 M = a r a 2 Y a 2 M , P 0 = 1 2 l l ^ [ c l c l ] 000 L = 0 , J = 1 , T = 0 , S 0 = 1 2 l l ^ [ c l c l ] 000 L = 0 , J = 0 , T = 1 ,
with M standing for the projection of angular momentum on the z-axis, a as a label of nucleons, and the brackets indicating the coupling of orbital angular momentum, spin, and isospin to different values, each of which has zero z-projection. The operator c l in Equation (5) creates a particle occupying the single-particle orbital with an orbital angular momentum l, and l ^ 2 l + 1 . In addition, the operator P 0 ( S 0 ) creates a correlated isoscalar (isovector) proton–neutron pair in the single-particle level l. To include the proton–neutron pairing effect, we start from a Bogoliubov transformation that mixes protons and neutrons in the quasiparticle creation and annihilation operators, i.e., (schematically),
α u p c p + v p c p + u n c n + v n c n .
In the practical calculations, the full equations should sum over single-particle states in the valence space, and hence we replace each of the coefficients u s and v s in Equation (6) with matrices U and V, which can be found in Ref. [44].
We then solve the constrained Hartree-Fock-Bogoliubov (HFB) equations for the effective Hamiltonian with linear constraints:
H = H eff λ N ( N N N ) λ Z ( N Z Z ) i λ i ( O i q i ) ,
where the N N and N Z signify the neutron and proton number operators, while λ N and λ Z stand for corresponding Lagrange multipliers. The sum over i includes the quadrupole operators Q 20 and Q 22 , with the addition of isoscalar or isovector proton–neutron pairing operator in Equation (4). The λ i represents the Lagrange multipliers that constrain the expectation values of these operators to specified quantities of q i . We solve the HFB equation repeatedly. In each time, the HFB vacuum is constrained to a different mesh point in the space of q i .
After we get a set of HFB vacua with different amounts of axial quadrupole deformation, triaxial quadrupole deformation, and isoscalar (or isovector) proton–neutron pairing amplitude, the GCM state can be composed of a linear superposition of the projected HFB vacua, given by
| Ψ Z N σ J = K , q f q σ J K | J M K ; Z N ; q ,
where | J M K ; Z N ; q P ^ M K J P ^ Z P ^ N | Φ ( q ) . The P ^ s are so-called projection operators which project quasiparticle vacua onto definitive angular momentum J and its z-component M, proton number Z, and neutron number N [45]. f q σ J K are the weight functions, where σ is simply a enumeration index. They can be taken as variational parameters and thus be computed by solving the Hill-Wheeler-Griffin equation [45]:
K , q H K K J ( q ; q ) E σ J N K K J ( q ; q ) f q σ J K = 0 ,
where the Hamiltonian kernel H K K J ( q ; q ) and the norm kernel N K K J ( q ; q ) are given by:
H K K J ( q ; q ) = Φ ( q ) | H eff P ^ K K J P ^ Z P ^ N | Φ ( q ) , N K K J ( q ; q ) = Φ ( q ) | P ^ K K J P ^ Z P ^ N | Φ ( q ) .
To solve the Hill-Wheeler-Griffin equation, we diagonalize the norm kernel N first:
K , q N K K J q ; q u K k J q = n k J u K k J q .
The nonzero eigenvalues n k J and corresponding eigenvectors u K k J q can be used to construct a set of orthonormal basis called “natural states”, which are defined by
| k J = K , q u K k J q n k J | J M K ; Z N ; q .
The Hamiltonian can be diagonalized in the space of these natural states, and the Hill-Wheeler-Griffin equation thus can be transformed into a normal eigenvalue problem. Then we can obtain the wave functions of GCM states | Ψ N Z σ J (see details in Refs. [46,47]). With the lowest J = 0 GCM states as ground states of the grandparent and grand-daughter nuclei, we can finally compute the 0 ν β β decay matrix element M 0 ν in Equation (1).

3. Calculations and Discussion

We start by using our GCM employing an effective Hamiltonian in a valence model space. For 48 Ca−Ti, we employ two sets of effective interactions to show the effect from taking account for a larger model space. One is the KB3G interaction restricted in one-major s d shell [48]. The other is the SDPFMU-DB interaction in the two-major s d p f shell [9]. For 76 Ge−Se, we set the valence space in the so-called f 5 p g 9 space, which comprises the 0 f 5 / 2 , 1 p 3 / 2 , 1 p 1 / 2 , and 0 g 9 / 2 orbits. We have employed the GCN2850 Hamiltonian that are fine tuned for relevant 0 ν β β decay candidate isotopes [49]. For the calculations of 124 Sn−Te, 130 Te−Xe, and 136 Xe−Ba, we use a recently proposed shell-model effective Hamiltonian called SVD Hamiltonian [50] for the j j 55 -shell configuration space that comprises the 0 g 7 / 2 , 1 d 5 / 2 , 1 d 3 / 2 , 2 s 1 / 2 , and 0 h 11 / 2 orbitals. It has been shown that these effective interactions account successfully for the low-lying level spectra, electromagnetic and Gamow-Teller transitions, and quadrupole deformations for these nuclei, respectively [7,8,9,49,50,51]. Since the effect on 0 ν β β decay NME from triaxial deformation is found to be minor [32], and most of the candidate isotopes show no theoretical or experimental evidence of being triaxially deformed, we neglect the triaxial quadrupole moment q 2 Q 22 in the choice of collective generator coordinates for simplicity. Instead, our GCM calculations use axial quadrupole moment q 1 Q 20 , as well as the proton–neutron pairing parameters q 3 1 / 2 P 0 + P 0 and q 4 1 / 2 S 0 + S 0 as generator coordinates.
We first examine the reliability of our Hamiltonian-based GCM calculation for the structural properties of the five pairs of 0 ν β β decay candidate nuclei of 48 Ca−Ti, 76 Ge−Se, 124 Sn−Te, 130 Te−Xe, and 136 Xe−Ba, which are compared with the calculated results given by the interacting shell model (ISM) calculations based on the same effective interaction without the additional λ Q ^ · Q ^ term in Figure 1. In principal, if the Hamiltonian-based GCM and the ISM itself employ the same effective interaction in the same model space, the ISM result can therefore be taken as the “exact” solution and an adequate benchmark, because it diagonalizes the Hamiltonian exactly. A great agreement between GCM and ISM calculations for all the candidate nuclei has been shown, except for the doubly magical nucleus 48 Ca in which the GCM calculation presents a slight overestimation in B ( E 2 ; 0 1 + 2 1 + ) value. The lowest excited states in 48 Ca should mainly be of particle-hole excitation, the description of which requires the incorporation of non-collective configurations.
Before investigating the correlation of the neutrinoless double- β decay NME with the low-lying electric quadrupole (E2) strength in the 0 + 2 + transitions, we make one more test. Figure 2 shows the linear anti-correlation between the excitation energies of the first 2 + states and the corresponding reduced E2 transition probabilities B ( E 2 ; 0 + 2 + ) . As the quadrupole–quadrupole interaction becomes stronger in the effective Hamiltonian by increasing the coupling constant λ in the extra λ Q ^ · Q ^ term, the 2 + energies are suppressed, while the B ( E 2 ; 0 1 + 2 1 + ) values are significantly enhanced. Qualitatively, we see a similar evolution for 48 Ti within one- and two-major-shell calculations, which employ the KB3G and SDPFMU-DB interactions, respectively. Lowered excitation energies, combining with the enhanced B ( E 2 ) values, indicate that the quadrupole collectivity of low-lying states is changed smoothly along with the modification of the additional λ Q ^ · Q ^ term, as we expect.
Figure 3 plots the changes of 0 ν β β decay NMEs as a function of the B ( E 2 ; 0 + 2 + ) values for 48 Ti. Though it is not explicitly marked in the figure, the larger constant λ always corresponds to the larger calculated B ( E 2 ; 0 + 2 + ) values. A remarkable anti-correlation can be obtained between the 0 ν β β decay NMEs and the B ( E 2 ) values. The same relation can be seen in both one- and two-shell calculations. The only difference is the 0 ν β β decay NME for 48 Ca→Ti decay calculated in the s d p f shell is slightly larger than that in the p f shell, which is in accordance with the large-scale SM calculations [9] with the omission of some cross-shell excitations. Nevertheless, the downward slope of the anti-correlation between the 0 ν β β decay NMEs and the B ( E 2 ) values are almost the same in the the p f -shell and the s d p f -shell calculations. The result suggests that enlargement of the model space may not dramatically affect the anti-correlation between NME and quadrupole collectivity, though gradual but continual effect arising from the addition of successive shells should not be ruled out.
We then turn to heavier candidate nuclei, e.g., 124 Sn−Te, 130 Te−Xe, and 136 Xe−Ba pairs. The anti-correlations are clearly presented in these candidate pairs, as shown in Figure 4. Similarly, larger constant λ leads to larger calculated B ( E 2 ; 0 + 2 + ) value, and results in smaller 0 ν β β decay NME. Since 0 ν β β decay occurs from the ground state of grand-parent nucleus to that of the grand-daughter nucleus, one may ask that what the correlation would be between 0 ν β β decay NME and the quadrupole collectivity of grand-parent nuclei. We thus plot the changes of 0 ν β β decay NMEs as a function of the B ( E 2 ; 0 + 2 + ) values for both the grand-parent nucleus 136 Xe and the grand-daughter nucleus 136 Ba in the right panel of Figure 4. Robust anti-correlations are exhibited in both 136 Xe and 136 Xe, though the downward slope is steeper in the grand-parent nucleus. This is due to the fact that 136 Xe has 82 neutrons, which is a magical number, resulting in the characteristics of a spherical nucleus. Meanwhile, 136 Ba is moderately deformed. The additional quadrupole–quadrupole interaction can barely increase the quadrupole deformation for the ground state of 136 Xe owing to the robust N = 82 shell closure, but can remarkably enhance the quadrupole collectivity for 136 Ba. A similar situation occurs in 48 Ca−Ti, 124 Sn−Te, and 130 Te−Xe, since in all these candidate nuclei pairs, grand-parent nuclei are spherical or weakly deformed, while grand-daughter nuclei have relatively larger deformation.
One may ask whether the anti-correlation between the 0 ν β β decay NMEs and the E2 transition strength merely result from the increasing difference of quadrupole collectivity between near-spherical grand-parent and deformed grand-daughter nuclei. To answer this question, we apply our analysis to 76 Ge−Se candidate nuclei, since, in this pair, both grand-parent and grand-daughter nuclei are well-deformed. The results are presented in Figure 5. To illustrate the effect from enhanced quardupole-collectivity difference between grand-parent and grand-daughter nuclei, we apply two different analyses, namely, (i) the λ Q ^ · Q ^ term is added in the effective Hamiltonian used in the GCM calculations for both 76 Ge and 76 Se, (ii) the λ Q ^ · Q ^ term is added only for grand-daughter nucleus 76 Se. Apparently, the latter case would increase the difference of quadrupole collectivity between 76 Ge and 76 Se. It is shown that, though the enhanced deformation difference between 76 Ge and 76 Se leads to a steeper downward slope of the 0 ν β β decay NME against the reduced E2 transition probability in the latter analysis, the anti-correlation between them still exists in the first analysis, in which the quadrupole collectivity is enhanced almost equally in 76 Ge and 76 Se.
To show the thorough interrelation between the 0 ν β β decay NMEs and the E2 strengths, Figure 6 summarizes our calculated results. As noted, the anti-correlation of 0 ν β β decay NMEs with the E2 strengths exists for investigated candidate nuclei in a universal way. The downward slopes are similar in both the p f -shell and s d p f -shell case of 48 Ca−Ti, as well as in the case of 48 Ca−Ti. The anti-correlated trend is even more drastic in the heavier candidates 124 Sn−Te, 130 Te−Xe, and 136 Xe−Ba. We also plot the adopted values of B ( E 2 ; 0 + 2 + ) from experiments [52], shown as vertical shades. The width of the shade displays the uncertainty of measured reduced E2 transition probability. It can be seen that the calculated B ( E 2 ; 0 + 2 + ) value without the artificial quadrupole–quadrupole interaction, which is denoted by the starting point of each curve, shows non-negligible discrepancy when compared with adopted values. It should be noted that, the more accurate correlation between 0 ν β β decay NME and E2 strength can be determined if we can obtain a better description of reduced E2 transition probability, and hence potentially help us to reduce the uncertainty of the 0 ν β β decay NME.

4. Conclusions

In this paper, we apply a detailed analysis of the interrelation between 0 ν β β decay NMEs and E2 transition rates B ( E 2 ; 0 + 2 + ) , using the Hamiltonian-based GCM built on the configuration mixing of symmetry-restored intrinsic basis states. The calculated 0 ν β β decay NMEs are clearly anti-correlated with the calculated transition rates of the collective quadrupole excitation from the ground state to the first 2 + state in response to artificial changes of quadrupole–quadrupole contributions in effective Hamiltonians. The anti-correlation is more remarkable in the decay from a spherical or weakly deformed grand-parent nucleus than a more deformed grand-daughter nucleus. Therefore, we conclude that a reliable description of the reduced collective E2 transition probability would be useful for reducing the uncertainties of the 0 ν β β decay NME.
Another interesting study would be the correlation of 0 ν β β decay NME with E2 strength for 150 Nd−Sm decay process. It is because 150 Nd−Sm is the only candidate nuclei pair of which the grandparent 150 Nd has larger quadrupole deformation compared with the grand-daughter 150 Sm. However, to compute the NME of the heaviest candidate 150 Nd−Sm with Hamiltonian-based GCM, the first issue that must be grappled with is what to use for the effective Hamiltonian. Since 150 Nd consists of 60 protons and 90 neutrons, the calculation should be carried out in the valence space between the proton and neutron 50 to 126 shell closures. The effective Hamiltonian for this extremely large model space is beyond the current capability. To overcome this difficulty, an implementation of multiple-shell valence-space Hamiltonians derived from non-perturbative ab initio methods, such as in-medium similarity renormalization group (IM-SRG) [53] or CC [54] method, would be of great importance. Recently, the valence-space Hamiltonians derived from both IM-SRG [55,56] and CC method [57] has been used to explore the open-shell nuclei in one single shell. The extension of the ab initio valence-space Hamiltonian to multiple shells is still in progress.
Finally, it should be mentioned that we only discuss the NME of 0 ν β β decay from the ground-state of grandparent nucleus to the ground-state of grand-daughter nucleus here. Actually, the decay to the low-lying excited states of the grand-daughter nucleus should also be taken into account if it is allowed energetically. The NME of this process is expected to be strongly suppressed by the phase-space factor in the standard light left-handed Majorana neutrino exchange mechanism [58,59]. Nevertheless, it may considerably contribute to the NME in the non-standard mechanism. Further investigation of 0 ν β β decay NME to the lowest 2 1 + state of grand-daughter nucleus within the framework of Hamiltonian-based GCM would be a desirable next step in future work.

Author Contributions

Conceptualization, C.J., C.Y. and J.Y.; methodology, C.J.; numerical calculations and analyses, C.J.; writing—original draft preparation, C.J.; writing—review and editing, C.Y. and J.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant Nos. 12275369 and 12141501) and the Fundamental Research Funds for the Central Universities, Sun Yat-sen University (Grant No. 22qntd3101). C. X. Yuan acknowledges the support from the Guangdong Major Project of Basic and Applied Basic Research 2021B0301030006.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Acknowledgments

We thank C. W. Johnson for fruitful discussions.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
0 ν β β Neutrinoless double- β decay
NMENuclear matrix element
E2Electric quadrupole
ISMInteracting shell model
IBMInteracting boson method
QRPAQuasiparticle random phase approximation
GCMGenerator-coordinate method
EDFEnergy density functional
IM-GCMIn-medium generator-coordinate method
CCCoupled-cluster method
HFBHartree-Fock-Bogliubov
IM-SRGIn-medium similarity renormalization group
GTGamow-Teller
SRCShort range correlation

References

  1. Avignone, F.T.; Elliott, S.R.; Engel, J. Double beta decay, Majorana neutrinos, and neutrino mass. Rev. Mod. Phys. 2008, 80, 481–516. [Google Scholar] [CrossRef] [Green Version]
  2. Caurier, E.; Menéndez, J.; Nowacki, F.; Poves, A. Influence of Pairing on the Nuclear Matrix Elements of the Neutrinoless ββ Decays. Phys. Rev. Lett. 2008, 100, 052503. [Google Scholar] [CrossRef] [Green Version]
  3. Menéndez, J.; Poves, A.; Caurier, E.; Nowacki, F. Disassembling the nuclear matrix elements of the neutrinoless ββ decay. Nucl. Phys. A 2009, 818, 139–151. [Google Scholar] [CrossRef] [Green Version]
  4. Horoi, M.; Stoica, S. Shell model analysis of the neutrinoless double-β decay of 48Ca. Phys. Rev. C 2010, 81, 024321. [Google Scholar] [CrossRef] [Green Version]
  5. Horoi, M. Shell model analysis of competing contributions to the double-β decay of 48Ca. Phys. Rev. C 2013, 87, 014320. [Google Scholar] [CrossRef] [Green Version]
  6. Horoi, M.; Brown, B.A. Shell-Model Analysis of the 136Xe Double Beta Decay Nuclear Matrix Elements. Phys. Rev. Lett. 2013, 110, 222502. [Google Scholar] [CrossRef] [Green Version]
  7. Neacsu, A.; Horoi, M. Shell model studies of the 130Te neutrinoless double-β decay. Phys. Rev. C 2015, 91, 024309. [Google Scholar] [CrossRef] [Green Version]
  8. Horoi, M.; Neacsu, A. Shell model predictions for 124Sn double-β decay. Phys. Rev. C 2016, 93, 024308. [Google Scholar] [CrossRef] [Green Version]
  9. Iwata, Y.; Shimizu, N.; Otsuka, T.; Utsuno, Y.; Menéndez, J.; Honma, M.; Abe, T. Large-Scale Shell-Model Analysis of the Neutrinoless ββ Decay of 48Ca. Phys. Rev. Lett. 2016, 116, 112502. [Google Scholar] [CrossRef] [Green Version]
  10. Sen’kov, R.A.; Horoi, M. Shell-model calculation of neutrinoless double-β decay of 76Ge. Phys. Rev. C 2016, 93, 044334. [Google Scholar] [CrossRef] [Green Version]
  11. Barea, J.; Kotila, J.; Iachello, F. Limits on Neutrino Masses from Neutrinoless Double-β Decay. Phys. Rev. Lett. 2012, 109, 042501. [Google Scholar] [CrossRef] [Green Version]
  12. Barea, J.; Kotila, J.; Iachello, F. Nuclear matrix elements for double-β decay. Phys. Rev. C 2013, 87, 014315. [Google Scholar] [CrossRef] [Green Version]
  13. Barea, J.; Kotila, J.; Iachello, F. 0νββ and 2νββ nuclear matrix elements in the interacting boson model with isospin restoration. Phys. Rev. C 2015, 91, 034304. [Google Scholar] [CrossRef] [Green Version]
  14. Šimkovic, F.; Pantis, G.; Vergados, J.D.; Faessler, A. Additional nucleon current contributions to neutrinoless double β decay. Phys. Rev. C 1999, 60, 055502. [Google Scholar] [CrossRef] [Green Version]
  15. Rodin, V.A.; Faessler, A.; Šimkovic, F.; Vogel, P. Uncertainty in the 0νββ decay nuclear matrix elements. Phys. Rev. C 2003, 68, 044302. [Google Scholar] [CrossRef] [Green Version]
  16. Kortelainen, M.; Civitarese, O.; Suhonen, J.; Toivanen, J. Short-range correlations and neutrinoless double beta decay. Phys. Lett. B 2007, 647, 128–132. [Google Scholar] [CrossRef] [Green Version]
  17. Kortelainen, M.; Suhonen, J. Improved short-range correlations and 0νββ nuclear matrix elements of 76Ge and 82Se. Phys. Rev. C 2007, 75, 051303. [Google Scholar] [CrossRef] [Green Version]
  18. Šimkovic, F.; Faessler, A.; Rodin, V.; Vogel, P.; Engel, J. Anatomy of the 0νββ nuclear matrix elements. Phys. Rev. C 2008, 77, 045503. [Google Scholar] [CrossRef] [Green Version]
  19. Šimkovic, F.; Faessler, A.; Müther, H.; Rodin, V.; Stauf, M. 0νββ-decay nuclear matrix elements with self-consistent short-range correlations. Phys. Rev. C 2009, 79, 055501. [Google Scholar] [CrossRef] [Green Version]
  20. Šimkovic, F.; Faessler, A.; Vogel, P. 0νββ nuclear matrix elements and the occupancy of individual orbits. Phys. Rev. C 2009, 79, 015502. [Google Scholar] [CrossRef] [Green Version]
  21. Mustonen, M.T.; Engel, J. Large-scale calculations of the double-β decay of 76Ge,130Te,136Xe, and 150Nd in the deformed self-consistent Skyrme quasiparticle random-phase approximation. Phys. Rev. C 2013, 87, 064302. [Google Scholar] [CrossRef] [Green Version]
  22. Šimkovic, F.; Rodin, V.; Faessler, A.; Vogel, P. 0νββ and 2νββ nuclear matrix elements, quasiparticle random-phase approximation, and isospin symmetry restoration. Phys. Rev. C 2013, 87, 045501. [Google Scholar] [CrossRef] [Green Version]
  23. Faessler, A.; González, M.; Kovalenko, S.; Šimkovic, F. Arbitrary mass Majorana neutrinos in neutrinoless double beta decay. Phys. Rev. D 2014, 90, 096010. [Google Scholar] [CrossRef] [Green Version]
  24. Hyvärinen, J.; Suhonen, J. Nuclear matrix elements for 0νββ decays with light or heavy Majorana-neutrino exchange. Phys. Rev. C 2015, 91, 024613. [Google Scholar] [CrossRef]
  25. Fang, D.L.; Faessler, A.; Šimkovic, F. 0νββ-decay nuclear matrix element for light and heavy neutrino mass mechanisms from deformed quasiparticle random-phase approximation calculations for 76Ge,82Se,130Te,136Xe, and 150Nd with isospin restoration. Phys. Rev. C 2018, 97, 045503. [Google Scholar] [CrossRef] [Green Version]
  26. Rodríguez, T.R.; Martínez-Pinedo, G. Energy Density Functional Study of Nuclear Matrix Elements for Neutrinoless ββ Decay. Phys. Rev. Lett. 2010, 105, 252503. [Google Scholar] [CrossRef] [Green Version]
  27. Vaquero, N.L.; Rodríguez, T.R.; Egido, J.L. Shape and Pairing Fluctuation Effects on Neutrinoless Double Beta Decay Nuclear Matrix Elements. Phys. Rev. Lett. 2013, 111, 142501. [Google Scholar] [CrossRef] [Green Version]
  28. Rodríguez, T.R.; Martínez-Pinedo, G. Neutrinoless ββ decay nuclear matrix elements in an isotopic chain. Phys. Lett. B 2013, 719, 174–178. [Google Scholar] [CrossRef] [Green Version]
  29. Yao, J.M.; Song, L.S.; Hagino, K.; Ring, P.; Meng, J. Systematic study of nuclear matrix elements in neutrinoless double-β decay with a beyond-mean-field covariant density functional theory. Phys. Rev. C 2015, 91, 024316. [Google Scholar] [CrossRef] [Green Version]
  30. Yao, J.M.; Engel, J. Octupole correlations in low-lying states of 150Nd and 150Sm and their impact on neutrinoless double-β decay. Phys. Rev. C 2016, 94, 014306. [Google Scholar] [CrossRef] [Green Version]
  31. Hinohara, N.; Engel, J. Proton-neutron pairing amplitude as a generator coordinate for double-β decay. Phys. Rev. C 2014, 90, 031301. [Google Scholar] [CrossRef] [Green Version]
  32. Jiao, C.F.; Engel, J.; Holt, J.D. Neutrinoless double-β decay matrix elements in large shell-model spaces with the generator-coordinate method. Phys. Rev. C 2017, 96, 054310. [Google Scholar] [CrossRef] [Green Version]
  33. Jiao, C.F.; Horoi, M.; Neacsu, A. Neutrinoless double-β decay of 124Sn, 130Te, and 136Xe in the Hamiltonian-based generator-coordinate method. Phys. Rev. C 2018, 98, 064324. [Google Scholar] [CrossRef] [Green Version]
  34. Jiao, C.F.; Johnson, C.W. Union of rotational and vibrational modes in generator-coordinate-type calculations, with application to neutrinoless double-β decay. Phys. Rev. C 2019, 100, 031303. [Google Scholar] [CrossRef] [Green Version]
  35. Yao, J.M.; Bally, B.; Engel, J.; Wirth, R.; Rodríguez, T.R.; Hergert, H. Ab Initio Treatment of Collective Correlations and the Neutrinoless Double Beta Decay of 48Ca. Phys. Rev. Lett. 2020, 124, 232501. [Google Scholar] [CrossRef]
  36. Novario, S.; Gysbers, P.; Engel, J.; Hagen, G.; Jansen, G.R.; Morris, T.D.; Navrátil, P.; Papenbrock, T.; Quaglioni, S. Coupled-Cluster Calculations of Neutrinoless Double-β Decay in 48Ca. Phys. Rev. Lett. 2021, 126, 182502. [Google Scholar] [CrossRef]
  37. Engel, J.; Menéndez, J. Status and future of nuclear matrix elements for neutrinoless double-beta decay: A review. Rep. Prog. Phys. 2017, 80, 046301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Vogel, P. Nuclear structure and double beta decay. J. Phys. G Nucl. Part. Phys. 2012, 39, 124002. [Google Scholar] [CrossRef] [Green Version]
  39. Yao, J.; Meng, J.; Niu, Y.; Ring, P. Beyond-mean-field approaches for nuclear neutrinoless double beta decay in the standard mechanism. Prog. Part. Nucl. Phys. 2022, 126, 103965. [Google Scholar] [CrossRef]
  40. Caurier, E.; Nowacki, F.; Poves, A. Nuclear-structure aspects of the neutrinoless ββ-decays. Eur. Phys. J. A 2008, 36, 195–200. [Google Scholar] [CrossRef] [Green Version]
  41. Fang, D.L.; Faessler, A.; Rodin, V.; Šimkovic, F. Neutrinoless double-β decay of deformed nuclei within quasiparticle random-phase approximation with a realistic interaction. Phys. Rev. C 2011, 83, 034320. [Google Scholar] [CrossRef] [Green Version]
  42. Zelevinsky, V.; Auerbach, N.; Loc, B.M. Nuclear structure features of Gamow-Teller excitations. Phys. Rev. C 2017, 96, 044319. [Google Scholar] [CrossRef] [Green Version]
  43. Horoi, M.; Neacsu, A.; Stoica, S. Statistical analysis for the neutrinoless double-β-decay matrix element of 48Ca. Phys. Rev. C 2022, 106, 054302. [Google Scholar] [CrossRef]
  44. Goodman, A.L. Advances in Nuclear Physics. In Advances in Nuclear Physics; Negele, J.V., Vogt, E., Eds.; Plenum Press: New York, NY, USA, 1979; Volume 11, p. 263. [Google Scholar]
  45. Ring, P.; Schuck, P. The Nuclear Many-Body Problem; Springer: Berlin, Germany, 1980. [Google Scholar]
  46. Rodríguez, T.R.; Egido, J.L. Triaxial angular momentum projection and configuration mixing calculations with the Gogny force. Phys. Rev. C 2010, 81, 064323. [Google Scholar] [CrossRef]
  47. Yao, J.M.; Meng, J.; Ring, P.; Vretenar, D. Configuration mixing of angular-momentum-projected triaxial relativistic mean-field wave functions. Phys. Rev. C 2010, 81, 044311. [Google Scholar] [CrossRef] [Green Version]
  48. Poves, A.; Sánchez-Solano, J.; Caurier, E.; Nowacki, F. Shell model study of the isobaric chains A = 50, A = 51 and A = 52. Nucl. Phys. A 2001, 694, 157–198. [Google Scholar] [CrossRef] [Green Version]
  49. Gniady, A.; Caurier, E.; Nowacki, F. (unpublised).
  50. Qi, C.; Xu, Z.X. Monopole-optimized effective interaction for tin isotopes. Phys. Rev. C 2012, 86, 044323. [Google Scholar] [CrossRef] [Green Version]
  51. Caurier, E.; Martínez-Pinedo, G.; Nowacki, F.; Poves, A.; Zuker, A.P. The shell model as a unified view of nuclear structure. Rev. Mod. Phys. 2005, 77, 427–488. [Google Scholar] [CrossRef] [Green Version]
  52. Pritychenko, B.; Birch, M.; Singh, B.; Horoi, M. Tables of E2 transition probabilities from the first 2+ states in even–even nuclei. At. Data Nucl. Data Tables 2016, 107, 1–139. [Google Scholar] [CrossRef] [Green Version]
  53. Bogner, S.K.; Hergert, H.; Holt, J.D.; Schwenk, A.; Binder, S.; Calci, A.; Langhammer, J.; Roth, R. Nonperturbative Shell-Model Interactions from the In-Medium Similarity Renormalization Group. Phys. Rev. Lett. 2014, 113, 142501. [Google Scholar] [CrossRef] [Green Version]
  54. Hagen, G.; Papenbrock, T.; Hjorth-Jensen, M.; Dean, D.J. Coupled-cluster computations of atomic nuclei. Rep. Prog. Phys. 2014, 77, 096302. [Google Scholar] [CrossRef] [Green Version]
  55. Stroberg, S.R.; Hergert, H.; Holt, J.D.; Bogner, S.K.; Schwenk, A. Ground and excited states of doubly open-shell nuclei from ab initio valence-space Hamiltonians. Phys. Rev. C 2016, 93, 051301. [Google Scholar] [CrossRef] [Green Version]
  56. Stroberg, S.R.; Calci, A.; Hergert, H.; Holt, J.D.; Bogner, S.K.; Roth, R.; Schwenk, A. Nucleus-Dependent Valence-Space Approach to Nuclear Structure. Phys. Rev. Lett. 2017, 118, 032502. [Google Scholar] [CrossRef] [Green Version]
  57. Jansen, G.R.; Schuster, M.D.; Signoracci, A.; Hagen, G.; Navrátil, P. Open sd-shell nuclei from first principles. Phys. Rev. C 2016, 94, 011301. [Google Scholar] [CrossRef] [Green Version]
  58. Tomoda, T. 0+→2+ neutrinoless ββ decay of 76Ge. Nucl. Phys. A 1988, 484, 635–646. [Google Scholar] [CrossRef]
  59. Fang, D.L.; Faessler, A. Nuclear matrix elements for the 0νββ(0+→2+) decay of 76Ge within the two-nucleon mechanism. Phys. Rev. C 2021, 103, 045501. [Google Scholar] [CrossRef]
Figure 1. (Color online) Comparison of the calculated properties of low-lying states between our Hamiltonian-based GCM and SM for the 0 ν β β decay candidate nuclei of 48 Ca−Ti, 76 Ge−Se, 124 Sn−Te, 130 Te−Xe, and 136 Xe−Ba, including (left panel) excitation energies of 2 1 + states E x ( 2 1 + ) and (right panel) E2 transition strengths B ( E 2 ; 0 1 + 2 1 + ) . The SM results are taken from Refs. [7,8,9,49].
Figure 1. (Color online) Comparison of the calculated properties of low-lying states between our Hamiltonian-based GCM and SM for the 0 ν β β decay candidate nuclei of 48 Ca−Ti, 76 Ge−Se, 124 Sn−Te, 130 Te−Xe, and 136 Xe−Ba, including (left panel) excitation energies of 2 1 + states E x ( 2 1 + ) and (right panel) E2 transition strengths B ( E 2 ; 0 1 + 2 1 + ) . The SM results are taken from Refs. [7,8,9,49].
Symmetry 15 00552 g001
Figure 2. (Color online) Calculated excitation energies of the first 2 + state of 48 Ti is linearly anti-correlated with the calculated reduced E2 transition probability B ( E 2 ; 0 + 2 + ) . The calculations are performed in the p f shell (left panel) and s d p f shell (right panel), respectively.
Figure 2. (Color online) Calculated excitation energies of the first 2 + state of 48 Ti is linearly anti-correlated with the calculated reduced E2 transition probability B ( E 2 ; 0 + 2 + ) . The calculations are performed in the p f shell (left panel) and s d p f shell (right panel), respectively.
Symmetry 15 00552 g002
Figure 3. (Color online) Calculated 0 ν β β decay NMEs of 48 Ca−Ti are anti-correlated with the calculated reduced E2 transition probabilities B ( E 2 ; 0 + 2 + ) in both p f -shell calculation (left panel) and s d p f -shell calculation (right panel).
Figure 3. (Color online) Calculated 0 ν β β decay NMEs of 48 Ca−Ti are anti-correlated with the calculated reduced E2 transition probabilities B ( E 2 ; 0 + 2 + ) in both p f -shell calculation (left panel) and s d p f -shell calculation (right panel).
Symmetry 15 00552 g003
Figure 4. (Color online) Similar to Figure 3, but for candidate nuclei 124 Sn−Te, 130 Te−Xe (left panel) and 136 Xe−Ba (right panel), respectively.
Figure 4. (Color online) Similar to Figure 3, but for candidate nuclei 124 Sn−Te, 130 Te−Xe (left panel) and 136 Xe−Ba (right panel), respectively.
Symmetry 15 00552 g004
Figure 5. (Color online) Calculated 0 ν β β decay NMEs of 76 Ge−Se changing against the calculated reduced E2 transition probabilities B ( E 2 ; 0 + 2 + ) . The calculation adding λ Q ^ · Q ^ term for both 76 Ge and 76 Se is denoted by red diamonds and a dash line, while the calculation adding λ Q ^ · Q ^ term for only 76 Se is shown as cyan circles and a solid line.
Figure 5. (Color online) Calculated 0 ν β β decay NMEs of 76 Ge−Se changing against the calculated reduced E2 transition probabilities B ( E 2 ; 0 + 2 + ) . The calculation adding λ Q ^ · Q ^ term for both 76 Ge and 76 Se is denoted by red diamonds and a dash line, while the calculation adding λ Q ^ · Q ^ term for only 76 Se is shown as cyan circles and a solid line.
Symmetry 15 00552 g005
Figure 6. (Color online) Calculated 0 ν β β decay NMEs of all investigated candidate nuclei changing against the calculated reduced E2 transition probabilities B ( E 2 ; 0 + 2 + ) . The vertical shades represent the adopted values of B ( E 2 ; 0 + 2 + ) .
Figure 6. (Color online) Calculated 0 ν β β decay NMEs of all investigated candidate nuclei changing against the calculated reduced E2 transition probabilities B ( E 2 ; 0 + 2 + ) . The vertical shades represent the adopted values of B ( E 2 ; 0 + 2 + ) .
Symmetry 15 00552 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Jiao, C.; Yuan, C.; Yao, J. Correlation of Neutrinoless Double-β Decay Nuclear Matrix Element with E2 Strength. Symmetry 2023, 15, 552. https://doi.org/10.3390/sym15020552

AMA Style

Jiao C, Yuan C, Yao J. Correlation of Neutrinoless Double-β Decay Nuclear Matrix Element with E2 Strength. Symmetry. 2023; 15(2):552. https://doi.org/10.3390/sym15020552

Chicago/Turabian Style

Jiao, Changfeng, Cenxi Yuan, and Jiangming Yao. 2023. "Correlation of Neutrinoless Double-β Decay Nuclear Matrix Element with E2 Strength" Symmetry 15, no. 2: 552. https://doi.org/10.3390/sym15020552

APA Style

Jiao, C., Yuan, C., & Yao, J. (2023). Correlation of Neutrinoless Double-β Decay Nuclear Matrix Element with E2 Strength. Symmetry, 15(2), 552. https://doi.org/10.3390/sym15020552

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop