Deformation and the Nuclear Matrix Elements of the Neutrinoless \beta\beta Decay

Deformation and the Nuclear Matrix Elements of the Neutrinoless Decay

J. Menéndez \atqueA. Poves Departamento de Física Teórica and IFT-UAM/CSIC, Universidad Autónoma de Madrid, E-28049, Madrid, SpainIPHC, IN2P3-CNRS/Université Louis Pasteur BP 28, F-67037, Strasbourg Cedex 2, France    E. Caurier \atqueF. Nowacki Departamento de Física Teórica and IFT-UAM/CSIC, Universidad Autónoma de Madrid, E-28049, Madrid, SpainIPHC, IN2P3-CNRS/Université Louis Pasteur BP 28, F-67037, Strasbourg Cedex 2, France

In this talk I will review the “state of the art” of the calculations of the nuclear matrix elements (NME) of the neutrinoless double beta decays () for the nuclei Ca, Ge, Se, Sn, Te, Te and Xe in the framework of the Interacting Shell Model (ISM), and compare them with the NME’s obtained using the Quasi-particle RPA approach (QRPA). I will also discuss the effect of the competition between the pairing and quadrupole correlations in the value of these NME’s. In particular I will show that, as the difference in deformation between parent and grand daughter grows, the NME’s of both the neutrinoless and the two neutrino modes decrease rapidly.

1 Introduction

The discovery of neutrino oscillations in recent experiments at Super-Kamiokande [1], SNO [2] and KamLAND [3] has changed the old conception of neutrinos by proving that they are massive particles. According to the origin of their mass, neutrinos can be either Dirac or Majorana particles, the latter case being particularly interesting since it would imply an extension to the standard model of electroweak interactions, and, being neutrinos their own antiparticles in this scenario, lepton number conservation would be broken. Besides, it happens that the best way to detect one of these violating processes and consequently establish the Majorana character of the neutrinos would be detection of the neutrinoless double beta decay ().

Double beta decay is a very slow weak process. It takes place between two even-even isobars when the single beta decay is energetically forbidden or hindered by large spin difference. Two neutrinos beta decay is a second order weak process —the reason of its low rate—, and has been measured in a few nuclei. The decay is analog but needs neutrinos to be Majorana particles. With the exception of one unconfirmed claim [4, 5], it has never been observed, and currently there is a number of experiments either taking place [6, 7, 8] or expected for the near future —see e.g. ref. [9]— devoted to detect this processes and to set up firmly the nature of neutrinos.

Furthermore, decay is also sensitive to the absolute scale of neutrino mass, and hence to the mass hierarchy —at present, only the difference between different mass eigenstates is known. Since the half-life of the decay is determined, together with the masses, by the nuclear matrix element (NME) for this process, the knowledge of these NME’s is essential to predict the most favorable decays and, once detection is achieved, to settle the neutrino mass scale and hierarchy.

Two different and complementary methods are mainly used to calculate NME’s for decays. One is the family of the quasiparticle random-phase approximation (QRPA). This method has been used by different groups and a variety of techniques is employed, with results for most of the possible emitters [10, 11, 12]. This work concerns to the alternative, the interacting shell model (ISM) [13].

In previous works [14, 15], the NME’s for the decay were calculated taking into account only the dominant terms of the nucleon current. However, in ref. [16] it was noted that the higher order contributions to the current (HOC) are not negligible and it was claimed that they could reduce up to 20%-30% the final NME’s. Subsequently, other QRPA calculations [17, 18] have also taken into account these terms, although resulting in a somewhat smaller correction. These additional nucleon current contributions have been recently included for the first time in the ISM framework [19]. In addition, the short range correlations (SRC), are now modeled either by the Jastrow prescription or by the UCOM method [20].

2 Ism QRPA Nuclear Matrix Elements

The expression for the half-life of the decay can be written as [21, 22]:


where is the averaged neutrino mass, a combination of the neutrino masses due to the neutrino mixing matrix —as we see, the neutrino mass scale is directly related to the decay rate— and is a kinematic factor —dependent on the charge, mass and available energy of the process. is the NME object of study in this work.

The kinematic factor depends on the value of the coupling constant . Therefore we have to take this into account when comparing the values of NME’s obtained with different values. In these cases we will use a NME modified as:


These ’s are directly comparable between them no matter which was the value of employed in their calculation, since they share a common factor —that of . In this sense, the translation of ’s into half-lives is transparent. The QRPA results obtained with different values are already expressed in this way by the authors of refs. [23, 12] while the results of refs. [17, 18, 24] will be translated by us into the above form when compared with other results.

We have calculated the ISM NME’s both taking the UCOM and Jastrow ansatzs for the short range correlations. The former use the correlator of the channel [25], throughout the calculation. The correlator of the other important —even— channel is very similar to this one, and it should not make much difference on this result.

In figure 1 the ISM and QRPA results for the NME’s are compared within this UCOM treatment of the SRC. The same figure but considering Jastrow type SRC was shown in ref. [19], but, inadvertently, the QRPA results from refs. [17, 18, 24] obtained with =1.0 where not transformed properly. This is corrected in figure 2. By comparing both figures, it is confirmed that there is a common trend; when the nuclei that participate in the decay have a low level of quadrupole correlations, as in the decays of Sn and Xe, both approaches agree. The QRPA in a spherical basis seems not to be able to capture the totality of the quadrupole correlations when they are strong. As these correlations tend to reduce the NME’s, the QRPA produces NME’s that are too large in Ge, Se, Te, and Te. For both ISM and QRPA the only net effect of UCOM is an increase of the Jastrow results of about 20%.

Figure 1: The neutrinoless double beta decay ’s for ISM and QRPA calculations treating the SRC with the UCOM approach. Tu07 QRPA results from ref. [23] and Jy07 results from refs. [17, 18].
Figure 2: Same as fig.1 but with Jastrow type SRC. Tu07 QRPA results from ref. [12] and Jy08 results from ref. [24].

3 The Influence of Deformation in the NME’s

An important issue regarding decay is the role of pairing and deformation. It has recently been discussed in ref. [19] that the pairing interaction favors the decay and that, consequently, truncations in seniority, not including the anti-pairing-like effect of the missing uncoupled pairs, tend to overestimate the value of the NME’s. On the other hand, the NME is also reduced when the parent and grand-daughter nuclei have different deformations [26]. Thus, we have studied the interplay between both pairing and deformation, this is, to which extent a wave function in the laboratory frame, truncated in seniority, can capture the correlations induced by the quadrupole-quadrupole part of the nuclear interaction, and its eventual influence in NME’s. There is an extra motivation to pursue this study; the possibility to carry on an experiment with Nd, which is a well deformed nuclei, decaying into Sm which is a much less deformed one.

To study the interplay between pairing, seniority truncations, and quadrupole correlations we need first to decide how to measure the quadrupole correlations of the ground state. Our choice is to refer to non energy-weighted sum rule:


The operator represents the mass quadrupole. No effective “nuclear” charges are included. Using Kr as our test bench, we proceed to compute , first with our standing effective interaction and different seniority truncations (s means the maximum seniority allowed in the wave functions of parent and grand daughter nuclei). The results are drawn in Figure 3 as the black circles labeled . We can see that at —roughly, the implicit level of seniority truncation in the spherical QRPA— some 70% of the full quadrupole correlations are incorporated in the wave function. We would like to know how this behavior evolves when more correlations are enforced in the system. For this we recalculate the ground state of Kr with a new hamiltonian that consists of the standing one plus a quadrupole-quadrupole term , whose effect will be gauged by its influence in the sum rule. To have an idea of the relevant range of values of in this nucleus and valence space, we have gone to the limit of pure quadrupole-quadrupole interaction with degenerate single particle energies, getting fm. The results for and are also shown in Figure 3 ( corresponds to in Figs. 6-9 and in Fig. 10 to . It is evident in the figure that, as we try to increase the correlations, the truncation becomes more and more ineffective. For , only 57% of the full correlations are present, and for only 50%. The situation is different for Se: while for the values of as a function of seniority are similar, albeit a bit smaller than the Kr ones, for and there is scarcely any increase of the ground state correlations. This means also that, as we increase , the “deformation” of Kr grows, whereas that of Se remains constant.

Figure 3: Quadrupole correlations in the ground state of Kr as a function of the amount of quadrupole-quadrupole interaction added to the hamiltonian and of the maximum seniority permitted in the wave functions


Figure 4: Se Kr Gamow-Teller matrix element, , as a function of the maximum seniority of the wave functions, for different values of the strength of the extra quadrupole-quadrupole interaction.

This behavior offers us the opportunity of exploring the effect of the difference in deformation between parent and grand daughter in the NME’s. To this goal, we have computed the Gamow-Teller matrix element for different values of —the amount of extra quadrupole-quadrupole interaction— and —the maximum seniority allowed in the wave function—. The results are gathered in Figure 4. For , we observe that the Gamow-Teller matrix element grows as a function of . This may seem paradoxical, but is not, because at this seniority truncation, the only effect of adding more quadrupole-quadrupole interaction is to augment the pairing content of the wave functions, thus increasing . At , remains constant as a function of meaning that the minor increase of the correlations of Kr, that we have shown in Figure 3, is barely enough to compensate the increase of at . On the contrary, the full space results are sensitive to the difference in deformation —or, to be more precise, to the difference in the level of quadrupole correlations in the ground state— between parent and grand daughter. The effect goes in the direction of reducing the value of . In the decay, doubling the quadrupole correlations in Kr, roughly halves .

Figure 5: Se Kr NME, , as a function of the difference between the mass quadrupole sum rule between Kr and Se for a large number of different values of the strength of the added quadrupole quadrupole interaction.

We can go much further in the exploration of the deformation effects, calculating the NME of the decay for initial and final states computed with different amounts of supplementary quadrupole-quadrupole interaction. As before , we measure the quadrupole correlations by means of the mass quadrupole sum rule. The results of this search are plotted in Fig. 5. We observe that the NME decreases almost linearly as the difference of the sum rules for the final and initial states increases. In fact, the maximum values of the NME are reached when this difference is close to zero. On the contrary for large differences the NME can be extremely quenched. Notice that the values in the figure range between 0.07 and 2.7, a factor of 40 span. The =0 value is 2.18, indicating that the difference in quadrupole correlations between Kr and Se is not very large.

4 0 (Unphysical) Mirror Decays: A Case Study

We have also studied the transitions between mirror nuclei in order to have a clearer view of the role of deformations in the NME’s. These transitions have the peculiarity that the wave functions of the initial and final nuclei are identical (provided Coulomb effects are neglected) and consequently the interplay of the operator and of the nuclear wave functions in the NME may be easier to understand.

We have studied four parent nuclei with six valence protons and four valence neutrons, Mg, Cr, Ge and Xe, decaying into the four grand daughter nuclei with four valence protons and six valence neutrons, Si, Fe, Se and Ba. The valence spaces considered are the -shell, the -shell, rg (1p, 1p, 0f, 0g) and rh (0g, 1d, 1d, 2s, 0h). The starting interactions are USD [27], KB3 [28], GCN28.50 and GCN50.82 [29]. The deformation of the nuclei is modified as in the previous section, and it is quantified also in the same way.

Figure 6: Ge Se NME, , as a function of the strength of the added quadrupole quadrupole interaction. Equally deformed case; the same amount of extra QQ interaction is added to Ge and Se. On the right hand y axis the pairing and quadrupole sum rules are represented, normalized so that their maximum value is 1. Note the change of scale in the x axis at .
Figure 7: Same as previous figure, but now the additional quadrupole interaction is only added to Se. The normalized overlap between the initial and final states is also included.

The results for in the case of equally deformed initial and final nuclei, are shown on Fig. 6. There we see that, as the nuclei become more deformed, the NME and the pairing content of the wave function get smaller, while the quadrupole sum rule grows. All these changes are nearly linear for reasonable deformations and then the saturation is approached more smoothly. Note that the purely quadrupole interaction ( limit) gives a NME which is about a half of the value obtaind with no additional quadrupole.

Fig. 7 shows the same quantities than Fig. 6 but now only the final nucleus has been artificially deformed by adding an extra quadrupole-quadupole term. In addition, the overlap between initial and final wave functions has been included. We see that now, the reduction of the NME is more pronounced and, what is more interesting, that it follows closely the overlap between wave functions. This means that, if we write the final wave function as: , the operator only connects the parts of the wave functions that have the same deformation among themselves.

The behavior of the NME’s with respect to the difference of deformation between parent and grand daughter is common to all the other transitions between mirror nuclei that we have studied. Therefore we can submit that this is a robust result. However, when we consider the transitions between equally deformed nuclei, the evolution of the NME’s with the deformation that we have found in A=66 is only shared by the A=110 case. When the valence space is a full major oscillator shell — and — the situation is quite different. Indeed, what is observed is that the NME does not decrease for moderate values of but remains rather constant until a point —with large deformation— where its value increases significantly —up to 50%—. This is due to the fact that, at this point, the major contribution to the NME ceases to come only form the decay of pairs coupled to , since other values like , which usually have a contribution to the NME contrary to that of , reverse sign and grow until being comparable with this contribution, thus resulting in this notorious rise of the NME. In the -shell and -shell cases, the limit is equivalent to Elliott’s SU(3) limit, and, the fact that both the initial and final nuclei should belong to the same irrep of SU(3) may be the reason of the increase of the NME, but, for the moment we have not found a formal explanation.

5 2 (Unphysical) Mirror Decays: A Case Study

Very similar conclusions may be reached for the effect of deformation in decay. For instance, the mirror nuclei diagonal and non-diagonal NME’s are represented as in the case in Figs. 8 and 9 for the transition. We see that the figures resemble very much that of the previous section, with the only exception that, in the equally deformed, case, the lowering of the NME due to deformation is more pronounced.

Figure 8: Equally deformed Ge Se NME. Note the change of scale in the x axis at .
Figure 9: The same as the previous figure, but now the only nuclei calculated with additional quadrupole interaction is the final one. The normalized overlap between initial and final states is also included.
Figure 10: Influence of deformation in the Ti Cr decay.

If we look to a non mirror —but again fictitious— transition, for instance that of Ti Cr, we get the results shown in Fig. 10. If we compare with those of ref. [26], which are again the equivalent ones for the transition, we see that there are not substantial changes. In this sense, deformation seems to affect similarly to and decays.

6 Summary

After a brief discussion of the “state of the art” results for the nuclear matrix elements of the neutrinoless double beta decay in the context of the Interacting Shell Model and of the Quasiparticle Random Phase Approximation, we have analyzed the role of the pairing correlations and the deformation in the NME’s, concluding that seniority truncations are less reliable when the quadrupole correlations are large. Since the NME’s are reduced when the deformation of parent and grand daughter is different, a bad treatment of the quadrupole correlations can lead to an artificial enhancement of the NME’s in the transitions among nuclei with unequal deformations, or in the cases of nuclei with different –and large- amounts of quadrupole correlations.

This work has been supported by a grant of the Spanish Ministry of Education and Science, FPA2007-66069, by the IN2P3-CICyT collaboration agreements, by the Spanish Consolider-Ingenio 2010 Program CPAN (CSD2007-00042), and by the Comunidad de Madrid (Spain), project HEPHACOS P-ESP-00346.


  • [1] \NAMEFukuda Y., Hayakawa T., Ichihara E., Inoue K., Ishihara K., Ishino H., Itow Y., Kajita T., Kameda J., Kasuga S., Kobayashi K., Kobayashi Y., Koshio Y., Miura M., Nakahata M., Nakayama S., Okada A., Okumura K., Sakurai N., Shiozawa M., Suzuki Y., Takeuchi Y., Totsuka Y., Yamada S., Earl M., Habig A. \atqueKearns E., \INPhys. Rev. Lett. 8119981562.
  • [2] \NAMEAhmad Q. R., Allen R. C., Andersen T. C., D.Anglin J., Barton J. C., Beier E. W., Bercovitch M., Bigu J., Biller S. D., Black R. A., Blevis I., Boardman R. J., Boger J., Bonvin E., Boulay M. G., Bowler M. G., Bowles T. J., Brice S. J., Browne M. C., Bullard T. V., Bühler G., Cameron J., Chan Y. D., Chen H. H., Chen M., Chen X. \atqueCleveland B. T., \INPhys. Rev. Lett. 892002011301.
  • [3] \NAMEEguchi K., Enomoto S., Furuno K., Goldman J., Hanada H., Ikeda H., Ikeda K., Inoue K., Ishihara K., Itoh W., Iwamoto T., Kawaguchi T., Kawashima T., Kinoshita H., Kishimoto Y., Koga M., Koseki Y., Maeda T., Mitsui T., Motoki M., Nakajima K., Nakajima M., Nakajima T., Ogawa H., Owada K., Sakabe T. \atqueShimizu I., \INPhys. Rev. Lett. 902003021802.
  • [4] \NAMEKlapdor-Kleingrothaus H. V., Dietz A., Harney H. L. \atqueKrivosheina I. V., \INMod. Phys. Lett. A1620012409.
  • [5] \NAMEKlapdor-Kleingrothaus H. V., Krivosheina I. V., Dietz A. \atqueChkvorets O., \INPhys. Lett. B5862004198.
  • [6] \NAMEArnold R., Augier C., Baker J., Barabash A., Broudin G., Brudanin V., Caffrey A. J., Caurier E., Egorov V., Errahmane K., Etienvre A. I., Guyonnet J. L., Hubert F., Hubert P., Jollet C., Jullian S., Kochetov O., Kovalenko V., Konovalov S., Lalanne D., Leccia F., Longuemare C., Lutter G., Marquet C., Mauger F., Nowacki F., Ohsumi H., Piquemal F., Reyss J. L., Saakyan R., Sarazin X., Simard L., Simkovic F., Shitov Y., Smolnikov A., Stekl L., Suhonen J., Sutton C. S., Szklarz G., Thomas J., Timkin V., Tretyak V., Umatov V., Vala L., Vanushin I., Vasilyev V., Vorobel V. \atqueVylov T., \INPhys. Rev. Lett. 952005182302.
  • [7] \NAMEArnaboldi C., Artusa D. R., Avignone F. T., Balata M., Bandac I., Barucci M., Beeman J. W., Brofferio C., Bucci C., Capelli S., Carbone L., Cebrian S., Cremonesi O., Creswick R. J., de Waard A., Farach H. A., Fiorini E., Frossati G., Guardincerri E., Giuliani A., Gorla P., Haller E. E., McDonald R. J., Morales A., Norman E. B., Nucciotti A., Olivieri E., Pallavicini M., Palmieri E., Pasca E., Pavan M., Pedretti M., Pessina G., Pirro S., Previtali E., Risegari L., Rosenfeld C., Sangiorgio S., Sisti M., Smith A. R., Torres L. \atqueVentura G., \INPhys. Rev. Lett. 952005142501.
  • [8] \NAMEBloxham T., Boston A., Dawson J., Dobos D., Fox S. P., Freer M., Fulton B. R., Gössling C., Harrison P. F., Junker M., Kiel H., McGrath J., Morgan B., Münstermann D., Nolan P., Oehl S., Ramachers Y., Reeve C., Stewart D., Wadsworth R., Wilson J. R. \atqueZuber K., \INPhys. Rev. C762007025501.
  • [9] \NAMEAvignone F. T., Elliott S. R. \atqueEngel J., \INRev. Mod. Phys. 802008481.
  • [10] \NAMESuhonen J. \atqueCivitarese O., \INPhys. Rept. 3001998123.
  • [11] \NAMERodin V. A., Faessler A., Simkovic F. \atqueVogel P., \INNucl. Phys. A7662006107.
  • [12] \NAMERodin V. A., Faessler A., Simkovic F. \atqueVogel P., \INNucl. Phys. A7932007213.
  • [13] \NAMECaurier E., Martinez-Pinedo G., Nowacki F., Poves A. \atqueZuker A. P., \INRev. Mod. Phys. 772005427.
  • [14] \NAMERetamosa J., Caurier E. \atqueNowacki F., \INPhys. Rev. C511995371.
  • [15] \NAMECaurier E., Nowacki F., Poves A. \atqueRetamosa J., \INPhys. Rev. Lett. 7719961954.
  • [16] \NAMESimkovic F., Pantis G., Vergados J. D. \atqueFaessler A., \INPhys. Rev. C601999055502.
  • [17] \NAMEKortelainen M. \atqueSuhonen J., \INPhys. Rev. C752007051303.
  • [18] \NAMEKortelainen M. \atqueSuhonen J., \INPhys. Rev. C762007024315.
  • [19] \NAMECaurier E., Menendez J., Nowacki F. \atquePoves A., \INPhys. Rev. Lett. 1002008052503.
  • [20] \NAMEFeldmeier H., Neff T., Roth R. \atqueSchnack J., \INNucl. Phys. A632199861.
  • [21] \NAMEDoi M., Kotani T., Nishiura H., Okuda K. \atqueTakasugi E., \INPhys. Lett. B1031981219.
  • [22] \NAMEDoi M., Kotani T. \atqueTakasugi E., \INProg. Theor. Phys. Suppl. 8319851.
  • [23] \NAMESimkovic F., Faessler A., Rodin V. A., Vogel P. \atqueEngel J., \INPhys. Rev. C772008045503.
  • [24] \NAMESuhonen J. \atqueKortelainen M., \INInt. J. Mod. Phys. E1720081.
  • [25] \NAMERoth R., Hergert H., Papakonstantinou P., Neff T. \atqueFeldmeier H., \INPhys. Rev. C722005034002.
  • [26] \NAMECaurier E., Nowacki F. \atquePoves A., \INEur. Phys. J. A362008195.
  • [27] \NAMEWildenthal B. H., \INProg. Part. Nucl. Phys. 1119845.
  • [28] \NAMEPoves A. \atqueZuker A., \INPhys. Rep. 701981235.
  • [29] \NAMEGniady A., Caurier E. \atqueNowacki F., to be published.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description