# Constraints on Heavy Neutrino and SUSY Parameters Derived from the Study of Neutrinoless Double Beta Decay

Authors: Andrei Neacsu; Sabin Stoica
Journal: Advances in High Energy Physics
 (2014)

---

## Abstract

New constraints on the lepton number violating (LNV) parameters are derived from the analysis of the neutrinoless double beta (0νββ) decay in the hypothesis that this process would occur through the exchange of heavy neutrinos and/or SUSY particles. For derivation, we use new values of both phase space factors (PSFs) and nuclear matrix elements (NMEs) calculated with numerical codes developed recently, as well as the most recent experimental lifetimes. The NMEs are computed with a shell model (ShM) code for 48Ca, 76Ge, and 82Se nuclei, while at present similar ShM results are available only for the first nucleus. We compare our results with similar ones from literature, obtained with ShM, QRPA, and IBM-2 methods, and conclude that more results are needed for a relevant analysis on the validity of NMEs associated with these decay mechanisms.


---

## Body

## 1. Introduction

Neutrinoless double beta (0νββ) decay is a beyond standard model (BSM) process by which an even-even nucleus transforms into another even-even nucleus with the emission of two electrons/positrons but no antineutrinos/neutrinos in the final states. Its study is important since it would clarify the question about the lepton number violation, decide on the neutrinos character (are they distinguished or not from their antiparticles?), and give a hint on the scale of their absolute masses. Moreover, the study of 0νββ decay has a broader potential to search for other BSM phenomena, such as new neutrino flavors and exotic particles. The great interest in understanding these fundamental issues has led to many theoretical and experimental investigations of this process. The reader can find up-to-date information on these studies in several recent reviews [1–6], which also contain therein a comprehensive list of references in the domain.One of the open issues concerning the0νββ decay is related to the possible mechanisms responsible for its occurrence. The exchange of a light Majorana neutrino in the presence of left handed (LH) weak interaction is still the most popular and studied mechanism. However, other mechanisms have also been considered, for example, the exchange of heavy neutrinos [7, 8] and SUSY particles [3, 9], which could also contribute to the total 0νββ decay rate. On the other hand, constraints on heavy neutrino and SUSY parameters can be now derived at hadron colliders, as well, by analyzing same sign dilepton decay channels, which are triggered by processes that violate the lepton number conservation by two units (like 0νββ decay). Indeed, the CMS, ATLAS, and LHCb experiments at LHC include now in their data analysis the search of such channels, and the first results were already reported [10–13]. Thus, complementary information that allows the constraint of various LNV parameters can also be obtained now from high-energy experiments.For all mechanisms, the0νββ lifetimes can be expressed in a good approximation as a product of a phase space factor, a nuclear matrix element related to the nuclear structure of the parent and daughter nuclei, and a LNV parameter related to the BSM mechanism considered. These three factors can be identified in (1) and (2) from the next section. Hence, to extract reliable upper limits for the LNV parameters, we need accurate calculation of both PSFs and NMEs factors, as well as reliable measurements of the lifetimes.The largest uncertainties in the theoretical calculations for double beta decay (DBD) arise from the calculated values of the NMEs. For the LH light neutrino exchange mechanism, the NMEs are currently computed by several methods, the present most employed ones being the proton-neutron quasi-random phase approximation (pnQRPA) [14–20], interacting shell model (ISM) [21–25], interacting boson model (IBM) [26–28], projected Hartree Fock Bogoliubov (PHFB) [29], and energy density functional (EDS) [30] methods. At present, there are still large discrepancies between the NMEs values computed with different methods and by different groups, which have been widely discussed in literature (see, e.g., [5, 6]) for the light neutrino exchange mechanism. At present, for heavy neutrino and SUSY exchange mechanisms, there are fewer NME calculations, performed with QRPA [31, 32] and IBA-2 [28] and ShM methods [33, 34]. The discrepancies between the existing NME values associated with heavy neutrino or SUSY mechanisms are even larger than those in the case of LH light neutrino exchange, so there is a need for new calculations.In this paper, we report new constraints on the heavy neutrino and SUSY parameters derived from the analysis of the neutrinoless double beta decay of three experimentally interested nuclei,Ca48, Ge76, and Se82, in the hypothesis that this process would occur through the exchange of such particle. The computation of the NMEs is performed with a ShM code developed recently and presented in detail in [35, 36]. Also, we estimate the uncertainties introduced in calculations by the use of different NN interactions and SRC parameterizations and found that the calculated NME values depend significantly on these nuclear ingredients. We compare our results with other previous ones from literature performed with ShM, QRPA, and IBM-2 approaches and discuss the discrepancies. We remark that, until now, large scale ShM calculations of NMEs corresponding to these mechanisms are available only for Ca48, and they differ significantly from the IBM-2 calculations. In order to understand these discrepancies, more ShM calculations of the NMEs for these (less discussed) 0νββ decay mechanisms, including other nuclei, are needed.The values of the PSFs are taken from our recent work [37], which are very close to those reported in [38, 39]. We mention that PSFs for DBD have been calculated since long time [40, 41], but they were less discussed, being considered to be computed with enough precision. Recently, they were recalculated with improved approaches [38, 39, 42] and several discrepancies were revealed as compared to the old calculations, which should be taken into account for accurate predictions of 0νββ lifetimes and derivation of LNV parameters.The paper is organized as follows. In the next section, we present the formulae for the0νββ decay lifetimes together with the relevant expressions for the PSFs and NMEs which have to be computed for deriving the LNV parameters, for heavy neutrino and SUSY mechanisms. In Section 3, we present and discuss our results, while in the last section we formulate the conclusions of our work.


## 2. Formalism

In this section, we give a short description of the heavy neutrino and SUSY formalisms, displaying the PSF and NME formulas that have to be computed. As we already mentioned, the expressions of the0νββ lifetimes can be written as a product of three terms: a phase space, a nuclear matrix elements, and a term related to the corresponding LNV mechanism:
(1)(T1/20ν)-1=G0ν(Qββ,Z)∑k‍(|Mk0ν|2)(ηk)2,
where k is an index which denotes the mechanisms that can contribute to the 0νββ decay. In this work, we refer to the heavy neutrino and SUSY exchange mechanisms, so the lifetime reads
(2)(T1/20ν)-1=G0ν(Qββ,Z)(|MN0ν|〈ηN〉+|Mλ′0ν|〈ηλ′〉+|Mq0ν|〈ηq〉)2,
where G0ν is the phase space factor for this decay mode, depending on the energy decay Qββ and nuclear charge Z. MN,λ′,q0ν are the NMEs associated with the heavy neutrino (N) exchange mechanism and with the gluino (λ′) and squark-neutrino (q) exchange mechanisms. The last mechanisms may appear in SUSY theories with R-parity violation [3]. ηN,λ′,q are the corresponding coupling parameters associated with these mechanisms.For the0νββ decay, the expression of the PSFs reads
(3)G0ν=24gA4R2ln⁡⁡2∫mec2T0-mec2‍f11(0)w0νdϵ1,
where f11(0) are the solutions of the Dirac equation and
(4)w0ν=gA4(Gcos⁡⁡θC)416π5(mec2)2(ℏc2)(p1c)(p2c)ϵ1ϵ2,
where gA  (gA=1.25) is the axial-vector coupling constant, G=1.16637×10-5 GeV-2 is the Fermi constant, θC is the Cabbibo angle, me is the electron rest mass, while p1, p2 and ϵ1, ϵ2 are the electron momenta and energies, respectively. In (3), it is convenient to redefine the PSFs by a renormalization that eliminates the constant gA and correlates (by dividing by 4R2, R the nuclear radius) the dimension of G0ν with the NMEs which are dimensionless, such that the PSFs are expressed in [yr-1].Concerning the NMEs, their expressions can be written, in general, as a sum of three components:(5)M0ν=MGT0ν-(gVgA)2·MF0ν-MT0ν,
where MGT0ν, MF0ν, and MT0ν are the Gamow-Teller (GT), Fermi (F), and Tensor (T) components, respectively.Due to the two-body nature of the transition operator, the NMEs can be expressed as a sum of products of two-body transition densities (TBTDs) and matrix elements of the two-body transition operators for two-particle states, shortly, two-body matrix elements (TBMEs) [35, 36]. Consider
(6)Mα0ν=∑jpjp′jnjn′JπTBTD(jpjp′,jnjn′;Jπ)×〈jpjp′;Jπ∥τ-1τ-2O12α∥jnjn′;Jπ〉,
where |jj′;Jπ〉 represent the antisymmetrized two-particle states and Omnα are DBD transition operators (α=GT,F,T). Since the NN effective interaction can be treated by means of a central (single-particle) potential, the NMEs can be conveniently calculated using Moshinsky’s transformations between the relative and center of mass (CM) coordinates and the proper use of nuclear states in different coupling notations. The calculation of the matrix elements of these operators can be decomposed into products of reduced matrix elements within the spin and relative coordinates subspaces [21, 35]. The expressions of the two-body transition operators are
(7)O12GT=σ1·σ2H(r),O12F=H(r),O12T=23[σ1×σ2]2·rRH(r)C(2)(r^),
where σ are the Pauli spin operators and C(2)(r^) is the center of mass operator. The expressions for the neutrino potentials Hα in the case of heavy neutrino mechanisms are
(8)Hα(r)=2Rπmemp∫0∞jk(qr)hα(q2)q2dq,
where mp is the proton mass and jk (k=0 for α=GT,F, and k=2 for α=T) are the Bessel spherical functions. The expressions of hα(α=GT,F,T) are
(9)hF(q2)=GV2(q2),(10)hGT(q2)=GA2(q2)gA2[1-23q2q2+mπ2+13(q2q2+mπ2)2]+23GM2(q2)gA2q24mp2,(11)hT(q2)=GA2(q2)gA2[23q2q2+mπ2-13(q2q2+mπ2)2]+13GM2(q2)gA2q24mp2,
where mπ is the pion mass and
(12)GM(q2)=(μp-μn)GV(q2),
with (μp-μn)=4.71. The expressions (9)–(11) include the following nuclear effects: (i) finite nucleon size (FNS) effects through the nucleon form factors GA, GV and (ii) higher order current (HOC) effects by the second and third terms in the hGT expression and by the appearance of the hT contribution.TheGV and GA form factors which take into account the finite size of the nucleons effect are
(13)GA(q2)=gA(ΛA2ΛA2+q2)2,GV(q2)=gV(ΛV2ΛV2+q2)2.
For the vector and axial coupling constants, the majority of calculations take gV=1 and the unquenched value gA=1.25, while the values of the vector and axial vector form factors are ΛV=850MeV and ΛA=1086MeV [1], respectively.For the SUSY mechanisms, only theGT and T components of the M0ν are contributing. For the gluino exchange mechanism (SUSY1), the radial neutrino potentials Hα have a similar form as that for the heavy neutrino mechanism (9), but with different hα factors [3, 14]:
(14)Mλ′0ν=c1π(MGT1π-MT1π)+c2π(MGT2π-MT2π),
with
(15)c1π=-292fπmπ4mp3me(mu+md)gsFPgA2,c2π=118fπ2mπ4mp3me(mu+md)2gs2gA2,
where mπ is the charged-pion mass, 139.6 MeV, fπ=0.668mπ MeV, mu and md are the up and down quark masses (mu+md)=11.6MeV, mp=938.27MeV, me=0.511MeV, gs=13.4, and FP≈4.41 [3]. Computing the expressions from (15), we obtain c1π=-85.33 and c2π=365.72. Consider
(16)Hα1π(r)=-2Rπ∫0∞jk(qr)hα1π(q2)dq,Hα2π(r)=-4Rπ∫0∞jk(qr)hα2π(q2)dq,
with
(17)hα1π(q2)=GA2(q2)gA2q4/mπ41+q2/mπ2,hα2π(q2)=GA2(q2)gA2q4/mπ4(1+q2/mπ2)2.For the squark-neutrino mechanism (SUSY2),Hα has an expression similar to the light neutrino mechanism:
(18)Hα(r)=2Rπ∫0∞jk(qr)hα(q2)q+〈E〉Gα(q2)qdq,
where 〈E〉 is the energy used in the closure approximation and represents the average excitation energy of the states in the intermediate odd-odd nucleus that contribute to the decay. The hα factors have in this case the following expressions [3, 14]:
(19)hGT,T=-16mπ2me(mu+md)q2/mπ2(1+q2/mπ2)2.For computing the radial matrix elements〈nl|Hα|n′l′〉, we use the harmonic oscillator wave functions ψnl(lr) and ψn′l′(r) corrected by a factor [1+f(r)], which takes into account the SRCs induced by the nuclear interaction:
(20)ψnl(r)⟶[1+f(r)]ψnl(r).
For the correlation function, we take the functional form
(21)f(r)=-c·e-ar2(1-br2),
where a, b, and c are constants which have particular values in different parameterizations. In this work, we consider the Miller-Spencer (MS) [44], AV18 [45], and CD-Bonn [46, 47] parameterizations.


## 3. Numerical Results and Discussions

We calculate first the NMEs for the three nuclei, that is,Ca48, Ge76, and Se82, and the three mechanisms mentioned above, using the code described in detail in [35, 36]. For Ca48, we performed calculations with two different NN interactions, GXPF1A [51] and KB3G [52] interactions, while for the other two isotopes we used the JUN45 [53] NN interaction. In the case of Ca48, our model space includes the whole fp shell f7/2, p3/2, f5/2, and p1/2, while, for Ge76 and Se82, the model space is jj44, including the following orbitals f5/2, p3/2, p1/2, and g9/2. For the nuclear parameters involved in calculations, we used the values mentioned in the previous section and, for the nuclear radius, R=r0A1/3, we used r0=1.2fm. The results are presented in Tables 1, 2, and 3. One can see that there are significant differences between the results when using different NN interactions and different SRC parameterizations. Indeed, especially for the short-range exchange mechanisms, as the heavy neutrino and gluino (SUSY1) ones, one expects the results to be sensitive to the NN interaction and to the type of SRC parameterizations. The uncertainties in the NME calculated values due to the use of different NN interactions are in the range of (10–16)%. The uncertainties due to the use of different SRC parameterizations are larger; sometimes differences between results are larger than 100%, especially when comparing the NME values obtained with MS [44] parameterizations with those obtained with softer parameterizations, such as AV18 [45] and CD-Bonn [46, 47]. Further, we compare our results with other similar ones found in literature. To the best of our knowledge, explicit values of NMEs calculated with large scale ShM approaches are reported only in [34] and for one nucleus, Ca48. Our results agree within ~10% with those from [34], for all three mechanisms. At this point, it is worth mentioning that the GT and tensor contributions in (5) must have a relative opposite sign, such that the adding of the tensor contribution decreases the M0ν total value. Hence, one has to pay attention when applying formula (5) because the GT and tensor contributions can come out from numerical calculations with opposite signs. In this respect, we suspect that our NME values agree even better with those from [34]. For the other two isotopes, that is, Ge76 and Se82, we compare our (ShM) results with QRPA [31, 32, 43] and IBM-2 [28] results. For the heavy neutrino and SUSY2 mechanisms, the differences between our NME values and NME values calculated with QRPA and IBM-2 methods are significant, but they are similar (in magnitude) to those encountered usually in the case of light neutrino exchange mechanism. For the SUSY1 mechanism, our results compare differently with the QRPA results reported by Tuebingen group in [32, 43], the only ones found in literature. We note that the agreement is good when comparing with the results of [43], where the MS parameterizations of the SRCs are employed, and it is bad when comparing with those of [32], where AV18 and CD-Bonn parameterizations for SRCs are used. As a conclusion, it turns out that more results for the NMEs associated with heavy neutrino and SUSY mechanisms, performed with different methods and including more nuclei, are needed for a relevant analysis of their validity.Table 1
0
ν
β
β NMEs values for heavy neutrino mechanism and comparison to other results in the literature.
No SRC
Jastrow
AV-18
CD-BONN
48Ca (GXPF1A)
91.5
23.7
47.5
70.3
48Ca (KB3G)
106.8
29.7
56.8
82.8
48Ca [34]
52.9
75.5
48Ca [28]
16.3
46.3
76.0
76Ge (JUN45)
255.7
73.1
137.8
199.2
76Ge [31]a, [32]b
32.6a
233b
351b
76Ge [28]
48.1
107
163
82Se (JUN45)
237.3
66.1
126.9
184.5
82Se [31]a, [32]b
30.0a
226b
340b
82Se [28]
35.6
84.4
132
The indices a and b are used to indicate different calculations performed by the Tubingen group.Table 2
0
ν
β
β NMEs values for SUSY1 mechanism and comparison to other results in the literature.
No SRC
Jastrow
AV-18
CD-BONN
48Ca (GXPF1A)
643.7.2
242.4
417.2
548.2
48Ca (KB3G)
750.1
293.3
492.47
641.7
48Ca [34]
453
618
48Ca [43]
392
147
76Ge (JUN45)
1762.1
678.3
1156.2
1509.4
76Ge [43]a, [32]b
1831a
625a
587b
515b
82Se (JUN45)
1628.0
612.8
1062.5
1393.5
82Se [43]a, [32]b
1667a
583a
574b
504b
The indices a and b are used to indicate different calculations performed by the Tubingen group.Table 3
0
ν
β
β NMEs values for SUSY2 mechanism and comparison to other results in the literature.
M
0
ν
No SRC
Jastrow
AV-18
CD-BONN
48Ca (GXPF1A)
65.1
46.4
65.0
70.0
48Ca (KB3G)
72.9
51.8
72.9
78.6
48Ca [34]
81.8
86.7
76Ge (JUN45)
281.7
233.4
283.5
296.8
76Ge [32]
594
612
82Se (JUN45)
253.6
208.6
255.7
268.1
82Se [32]
578
595In Table4, we present new upper limits for the LNV parameters associated with the studied 0νββ decay mechanisms, in the hypothesis of a single dominance mechanism. The PSFs values are taken from [37] and we mention that they are very close to PSF values reported in [38, 39]. The PSFs have been calculated since long time [40, 41] but they were less discussed, being considered to be computed with enough precision. Recently, they were recalculated with improved approaches, by using exact electron Dirac wave functions, taking into account the finite nuclear size and electron screening effects [38] and, in addition, a realistic Coulomb potential [39, 42], and differences/discrepancies were found in several cases between the old and the recent PSF values. For the three isotopes considered in this work, the new PSF values differ from older calculations (e.g., from those reported in [41, 54–56]) up to 10%, justifying hence their use in a precise derivation of LNV parameters. The uncertainties in NME values reflect on the precision of derivation the LNV parameters. According to the consensus agreed upon in literature that softer SRC parameterizations are indicated, we adopt the NME values calculated with the CD-Bonn parameterization and derive the LNV parameters for the three nuclei and the three mechanisms considered. These results are presented in Table 4. For derivation, we used the most recent results for the experimental lifetimes, reported in references indicated in parenthesis. For Ca48, our limits for the LNV parameters are close to those reported in [34], while for the other two nuclei they differ from other results from literature [28, 32, 34, 43], derived with NMEs values calculated with QRPA and IBM-2 methods. As we already mentioned, information on LNV parameters for the heavy neutrino mechanism can be now extracted, as well, from the data analysis of the same sign dilepton channels at the LHC experiments. Thus, we expect that complementary information from both low- and high-energy experiments could be used in the future for better constraining of the LNV parameters.Table 4
Upper limits for Majorana neutrino mass parameters together with the other components of the0νββ decay halftimes: the experimental lifetimes lower limits, the phase space factors, and the nuclear matrix elements.
T
exp
⁡
0
ν [yr]
G
0
ν [yr−1]
M
N
0
ν
M
λ
′
0
ν
M
q
0
ν
〈
η
N
〉
〈
η
λ
′
〉
〈
η
q
〉
48Ca*
5.8
·
1
0
22
[48]
2.46
E
-
14
70.3
548.2
70
2.42
·
1
0
-
7
3.10
·
1
0
-
8
2.42
·
1
0
-
7
48Ca†
5.8
·
1
0
22 [48]
2.46
E
-
14
82.8
641.7
78.6
2.05
·
1
0
-
7
2.64
·
1
0
-
8
2.16
·
1
0
-
7
76
Ge
2.1
·
1
0
25 [49]
2.37
E
-
15
199.2
1509.4
296.8
0.14
·
1
0
-
7
0.19
·
1
0
-
8
0.10
·
1
0
-
7
82
Se
3.6
·
1
0
23 [50]
1.01
E
-
14
184.5
1393.5
268.1
0.58
·
1
0
-
7
0.76
·
1
0
-
8
0.40
·
1
0
-
7
*Denotes GXPF1A [51] effective interaction and †KB3G [52] effective interaction.


## 4. Conclusions

We report new constraints for the heavy neutrino and SUSY parameters derived from an analysis of0νββ decay of three experimentally interesting nuclei, Ca48, Ge76, and Se82. To the best of our knowledge, for the last two isotopes, no other large scale ShM results for these mechanisms are available so far. For deduction of these parameters, precise values of the NMEs and PSFs which enter the lifetime formulae are needed, as well as accurate experimental measurements of the lifetimes. The NME calculations are performed with a recent ShM code described widely in [35, 36]. The PSFs are taken from our most recent calculations reported in [37]. For SUSY mechanisms, we discuss two mechanisms associated with possible R-parity breaking in SUSY theories, a short-range one with exchange of heavy Majorana neutrino and scalar SUSY particles (gluinos and squarks and/or neutralinos and selectrons) and a long-range one involving exchange of both heavy squarks and light Majorana neutrinos—“squark-neutrino” mechanism. First, we calculated the NMEs using different nuclear ingredients. For Ca48, we use two different NN interactions, GXPF1A and KB3G, and found differences between the results within 10%. Also, we performed the calculations using three different parameterizations for the SRCs and found a significant dependence of the results on the way the SRCs are introduced in calculations. Further, we compare our (ShM) results with similar results from literature. For Ca48, our results are in good agreement with the results from [34], for all three mechanisms. For the other two isotopes, Ge76 and Se82, we compare our (ShM) NMEs with those computed by QRPA [31, 32, 43] and IBM-2 [28] methods. For the heavy neutrino and SUSY2 mechanisms, the differences between our results and the results from those references are significant but similar (up to a factor of two) to those encountered usually in literature in the case of light neutrino exchange mechanism. For the SUSY1 mechanism, our results compare differently with QRPA results, reported by Tuebingen group. There is a good agreement with the results from [43], when the MS SRC parameterization is employed, but the agreement is bad when comparing with the values from [32], when softer SRC parameterizations are employed. As a conclusion, it turns out that more results for the NMEs associated with these 0νββ decay mechanisms, alternative to the light neutrino exchange one, performed with different nuclear methods and including more nuclei, are needed for a relevant analysis of their validity. Then, we derived upper limits of the LNV parameters associated with the three mechanisms, in the hypothesis of one single dominance mechanism, using NMEs values computed with CD-Bonn parameterization. Finally, it is worth noting that information on LNV parameters can now be provided by LHC experiments at CERN from the analysis of the same sign dilepton channels. Thus, there is now the possibility to gather complementary information from both low- and high-energy experiments and use it to better constrain LNV parameters.



---

*Source: 724315-2014-06-16.xml*