Abstract
In the present work, the λ mechanism (left-right weak boson exchange) and the light neutrino-exchange mechanism of neutrinoless double beta decay is studied. In particular, much attention is paid to the calculation of nuclear matrix elements for one of the neutrinoless double beta decaying isotopes 82Se. The interacting shell model framework is used to calculate the nuclear matrix element. The widely used closure approximation is adopted. The higher-order effect of the pseudoscalar term of nucleon current is also included in some of the nuclear matrix elements that result in larger Gamow-Teller matrix elements for the λ mechanism. Bounds on Majorana neutrino mass and lepton number violating parameters are also derived using the calculated nuclear matrix elements.
1 Introduction
Neutrinoless double beta decay (0νββ) is a rare second-order weak nuclear process. In this process, neutrino comes as a virtual intermediate particle when two neutron pairs decay into two proton pairs inside some even-even nuclei. Thus, it violates the lepton number by two units. The 0νββ experiment is one of the possible ways to determine the effective neutrino mass (; ; ; ; ) and can help to solve many mysteries of neutrinos, such as whether neutrinos are their own anti-particle (Majorana neutrino) or not (Dirac neutrino) (; ; ).
As lepton number conservation is not exact in most of the beyond the standard model (BSM) physics theories, many particle mechanisms of 0νββ have been proposed in different BSM theories such as light neutrino-exchange mechanism (; ), heavy neutrino-exchange mechanism (), left-right symmetric mechanism (; ), and the supersymmetric particles exchange mechanism (; ).
The decay rate for any particle mechanism of 0νββ is connected by nuclear matrix elements (NMEs) and absolute neutrino mass. These NMEs are calculated using theoretical nuclear many-body models (). Popular nuclear models are quasiparticle random phase approximation (QRPA) (), the interacting shell-model (ISM) (; ; ; ; ), the interacting boson model (IBM) (; ), the generator coordinate method (GCM) (), the energy density functional (EDF) theory (; ), and the projected Hartree-Fock Bogolibov model (PHFB) (). Other techniques includes, ab initio calculations for lower mass nuclei (A = 6–12) using variational Monte Carlo (VMC) method (; ; ).
In the present work, we focus on the left-right weak boson (WL-WR) exchange λ mechanism along with the standard light neutrino-exchange mechanism (WL − WL exchange) of the 0νββ mediated by light neutrinos (; ; ). The λ mechanism has origin in the left-right symmetric mechanism with right-handed gauge boson at the TeV scale (). Thus, it will be interesting to study how the λ mechanism can compete with the standard light neutrino-exchange mechanism when both the mechanisms co-exist. Hence, in the present work, we are eager to study the λ and light neutrino-exchange mechanisms together.
In left-right symmetric model, there is another mass independent mechanism called η mechanism which occurs through WL − WR mixing. It will be interesting to study η mechanism along with λ mechanism of 0νββ. But, η mechanism is suppressed due to WL − WR mixing as compared to λ mechanism (). Hence, in the present work, we are interested to study the mass independent λ mechanism along with the mass dependent standard light neutrino-exchange mechanism. In future studies, we will extensively explore the η mechanism of 0νββ along with other mass independent and dependent mechanisms in left-right symmetric model.
One of the motivations of the present work is to include effects of some of the revisited formalism of Ref. () on light neutrino-exchange and λ mechanism of 0νββ. The revised formalism was exploited to include the effects of the pseudoscalar term of nucleon currents. Using the revised formalism of Ref. (), the NMEs for λ, and light neutrino-exchange mechanisms of 0νββ are calculated using the QRPA model for several 0νββ decaying isotopes using closure approximation in Ref. (). Most of the NMEs relevant for λ and light neutrino-exchange mechanisms are also calculated using ISM in Ref. () using the closure approximation for different 0νββ decaying isotopes (including 82Se). In this case, some of the NMEs are calculated without including the higher-order terms (for example, pseudoscalar and weak magnetism terms) of the nucleon currents. Recently, using the revised formalism of Ref. (), we have also calculated the NMEs for 48Ca in Ref. () using the non-closure approximation and found a significant change in some of the NMEs for including the pseudoscalar term. Thus, we have tried here to include the revised higher-order effect of the pseudoscalar term of nucleon current for the λ mechanism of 0νββ of 82Se using ISM. The 0νββ of 82Se is one of the experimental interests of CUPID (; ) and NEMO-3 () experiments. Hence, it is important to study the nuclear structure aspects of 0νββ of 82Se theoretically. In recent years, one of the most important studies on light neutrino-exchange 0νββ of 82Se was performed in the ISM framework in Ref. () using the non-closure approximation. Here we focus on the λ mechanism of 0νββ of 82Se in the closure approximation using the revised nucleon current term.
This paper is organized as follows. In Section 2, the expression for decay rate and the theoretical formalism to calculate NMEs for the λ and light neutrino-exchange mechanisms of 0νββ are presented. The results and discussion are presented in Section 3. A summary of the work is given in Section 4.
2 Theoretical Framework
2.1 Decay Rate for λ Mechanism of 0νββ
If both light neutrino-exchange (WL − WL exchange) and λ mechanisms (WL − WR exchange) of 0νββ co-exist, one can write the decay rate for 0νββ as (; )where the coupling constant λ is defined as ()
The and are masses of the Standard Model left-handed WL and right-handed WR gauge bosons, respectively. The ην of Eq. 1 is an effective lepton number violating parameters for WL − WL exchange, ηλ is an effective lepton number violating parameters for WL − WR exchange, and ψ denotes the CP violating phase. These parameters are given in Ref. () as
Here, mββ is the effective Majorana neutrino mass defined by the neutrino mass eigenvalues mj and the neutrino mixing matrix elements Uej ():
The U, and T are the 3 × 3 block matrices in flavor space, which constitute a generalization of the Pontecorvo-Maki-Nakagawa-Sakata matrix, namely the 6 × 6 unitary neutrino mixing matrix (; ).
The amplitude of λ mechanism is given by ()where λ is defined earlier, GF is the Fermi constant for weak interaction, and q is the virtual Majorana neutrino momentum.
The coefficients CI (I = mm, mλ and λλ) of Eq. 1 are linear combinations of products of nuclear matrix elements and phase-space factors ().
Calculated values of phase-space factors G0i (i = 1, 2, 3, 4, 10 and 11) for different 0νββ decaying nuclei are given in Ref. ().
2.2 Nuclear Matrix Elements for λ Mechanism of 0νββ
Matrix elements required in the expression of CI are ().
The (MGT,ωGT,qGT) (MF,ωF,qF), and (MT,ωT,qT) matrix elements of the scalar two-body transition operator of 0νββ can be expressed as ()where, 〉, and 〉 are the initial and the final 0+ ground state (g.s) for 0νββ decay, respectively, and α = (GT, F, T, ν, ωGT, ωF, ωT, νω, qGT, qF, qT, 1 +, 2−), τ− is the isospin annihilation operator. The scalar two-particle transition operators of 0νββ containing spin and radial neutrino potential operators can be written aswhere, , r = r1 −r2, and r = |r| is inter nucleon distance of the decaying nucleons. The Ek is the energy of the virtual intermediate state () of 0νββ. The intermediate state is achieved when one neutron from the initial state is converted into one proton. Subsequently, from the state, another neutron is converted into another proton to achieve the final state of the 0νββ. For the present manuscript, is the 0+ g.s. of 82Se, is the 0+ g.s. of 82Kr, and are all the allowed spin-parity states of intermediate nucleus 82Br.
There are two approximations for calculating the NME, one is non-closure approximation and another is the widely used closure approximation. In non-closure approximation, the radial neutrino potential Hα(r, Ek) has explicit dependence on energy of the intermediate state 〉. In non-closure approximation, the radial neutrino potential for λ mechanism of 0νββ are is given as integral over Majorana neutrino momentum q ():where R is the radius of the parent nucleus, and the fα(q, r) factor (Appendix B) contains the form factors that incorporates the effects of finite nucleon size (FNS), and higher-order currents (HOC) of nucleons (), which is given in Appendix B of the manuscript. The Ei and Ef are the g.s. energy of the initial and final nucleus of the 0νββ decay, respectively. The non-closure approximation is computationally very challenging, because in this approximation, the NME has explicit dependence on the energy of large numbers of virtual intermediate state |k⟩ and calculating these states requires enormous computational power. Particularly, for higher mass region isotopes, some of the calculations are still beyond the reach of current generation’s high-performance computers. Fortunately, the most of the contributions on NME of 0νββ come from low lying energy states up to 10–12 MeV of the intermediate nucleus (; ). Thus, one can replace the effects of Ek with a suitable constant energy called closure energy without affecting the value of NME too much, and this approximation is known as closure approximation. In this approximation, one assumes ()and the radial neutrino potential operator of Eq. 16 becomes
In closure approximation, the 0νββ decay operators defined in Eq. 15 become
The closure approximation is widely used in literature as it eliminates the complexity of calculating a large number of virtual intermediate states (; ; ). One can find suitable values of using the method described in Ref. (), such that using closure approximation, one can get NME near to the non-closure approximation.
In the calculation of the NME of 0νββ, it is also necessary to take into account the effects of short-range correlations (SRC). A standard method to include SRC is via a phenomenological Jastrow-like function (; ; ). Including SRC effect in the Jastrow approach, one can write the NME of 0νββ defined in Eq. 14 as ()where Jastrow-type SRC function is defined as
In literature, three different SRC prametrizationparameterization are used: Miller-Spencer, Charge-Dependent Bonn (CD-Bonn), and Argonne V18 (AV18) to parametrize a, b, and c (). These parameters are chosen in such a way that the two-body wave function of two-body matrix elements (TBME) for 0νββ are still normalized. The parameters a, b, and c in different SRC parametrizations are given in Table 1.
TABLE 1
| SRC type | a | b | c |
|---|---|---|---|
| Miller-Spencer | 1.10 | 0.68 | 1.00 |
| CD-Bonn | 1.52 | 1.88 | 0.46 |
| AV18 | 1.59 | 1.45 | 0.92 |
Parameters for the short-range correlation (SRC) parametrization of Eq. 21. Values are taken from Ref. ().
This approach of using a Jastrow-like function to include the effects of SRC is extensively used in Refs. (; ; ).
2.3 The Closure Method of Nuclear Matrix Elements Calculation for 0νββ in ISM
The (MGT,ωGT,qGT) (MF,ωF,qF), and (MT,ωT,qT) matrix elements of the scalar two-body transition operator of 0νββ can be expressed as the sum over the product of the two-body transition density (TBTD) and anti-symmetric two-body matrix elements ():where, α = (F, GT, T, ωF, ωGT, ωT, qF, qGT, qT), J is the coupled spin of two decaying neutrons or two final created protons, τ− is the isospin annihilation operator, A denotes that the two-body matrix elements (TBME) (Appendix A) are obtained using anti-symmetric two-nucleon wavefunctions, and k1 stands for the set of spherical quantum numbers (n1; l1; j1) (similar definition for k2, k1′, k2′). The |i⟩ is 0+ ground state (g.s.) of the parent nucleus, and |f⟩ is the 0+ g.s of the granddaughter nucleus.
The TBTD can be expressed as ()where,andare the two particle creation and annihilation operator of rank J, respectively. Most of the available public shell model code does not provide the option to calculate TBTD directly. One of the ways is to calculate TBTD in terms of a large number of two nucleon transfer amplitudes (TNA), assuming 0νββ decay occurs through (n − 2) channel (). In (n − 2) channel of 0νββ, the TNA are calculated with a large set of intermediate states |m⟩ of the (n − 2) nucleons system, where n is the number of nucleons for the parent nucleus. In this approach, the TBTD in terms of TNA is expressed as ()where, TNA are given by
Here, Jm is the spin of the allowed states |m⟩ of intermediate nuclei. J0 is spin of |i⟩ and |f⟩. Jm = J when J0 = 0 ().
3 Results and Discussion
We have used JUN45 effective shell model Hamiltonian () of fpg model space to calculate the relevant initial, intermediate, and final nuclear states for 0νββ of 82Se. In the fpg model space, valence nucleons can occupy the orbitals 0f5/2, 1p3/2, 1p1/2, and 0g9/2. For the 0νββ decay of 82Se through (n − 2) channel, the states of allowed spin-parity of 80Se acts as intermediate states for TNA calculations. The nuclear shell model code KSHELL () was used in the calculation. For comparing some of the TNA values, NushellX@MSU () code was also used. In the present calculation, we have included the first 100 states of different allowed spin-parity of 80Se in calculating the TNA. Earlier, it was found that considering around the first 50 states is enough to get the saturated value of NME, as the most dominating contributions come from the first few initial states (; ).
We have adopted the widely used closure approximation with the closure energy ⟨E⟩ = 0.5 MeV. Earlier studies of Refs. (; ) suggested that ⟨E⟩ = 0.5 MeV is a suitable value that is close to optimal closure energy and, thus, gives NME near to the NME in the non-closure approximation. The non-closure method can give the exact value of NME, but the present study is beyond the scope of studying it. But, according to earlier results (; ), with ⟨E⟩ = 0.5 MeV, one can get NME in the closure approximation close to the NME in non-closure approximation (within 1% difference).
Different types of NMEs for light neutrino-exchange and λ mechanism of 0νββ for 82Se is shown in Table 2. Here, NMEs are calculated in different SRC parameterization schemes. All standard effects of FNS + HOC are taken care of in all calculations. It is found that the Gamow-Teller matrix elements dominate over Fermi and tensor type matrix elements. Also, it is found that the MqGT type matrix element associated with the λ mechanism is relatively large as compared to standard light neutrino-exchange Gamow-Teller matrix element MGT. This leads to the large value of total NME M1+ for λ mechanism as compared to total NME Mν for light neutrino-exchange mechanism.
TABLE 2
| SRC Type | ||||
|---|---|---|---|---|
| NME Type | None | Miller-Spencer | CD-Bonn | AV18 |
| MF | −0.633 | −0.442 | −0.674 | −0.621 |
| MGT | 3.681 | 2.536 | 3.247 | 3.068 |
| MT | −0.020 | −0.020 | −0.020 | −0.020 |
| Mν | 3.529 | 2.790 | 3.645 | 3.433 |
| MωF | −0.630 | −0.441 | −0.671 | −0.618 |
| MωGT | 3.075 | 2.453 | 3.165 | 2.986 |
| MωT | −0.020 | −0.020 | −0.020 | −0.020 |
| Mνω | 3.485 | 2.751 | 3.599 | 3.388 |
| MqF | −0.330 | −0.274 | −0.384 | −0.372 |
| MqGT | 11.667 | 10.167 | 12.538 | 12.184 |
| MqT | −0.097 | −0.097 | -0.097 | −0.097 |
| M1+ | 11.636 | 10.241 | 12.409 | 12.076 |
| M2− | 2.192 | 1.613 | 2.220 | 2.046 |
NMEs for 0νββ (light neutrino-exchange and λ mechanism) of 82Se.
This increment of MqGT type of NME, which is obtained through the new revised expression of the nucleon currents of Ref. (), is surprisingly high. It is coming through the new revised expression of the nucleon currents of Ref. () which includes the higher-order term (pseudoscalar) of the nucleon currents. In our calculation, Eq. 39 is used to calculate MqGT type NME using the revised formalism of nucleon currents of Refs. (; ).
An old equivalent expression of Eq. 39 is also found in Ref. (), which one can write using Eq. (A2c) and Eq. (A4b) of Ref. () as
Using this old value of fqGT (q, r), the MqGT type NME will be significantly smaller, as reported earlier.
Here we include the higher-order current effect of pseudoscalar term in Eq. 39 as suggested in Ref. () which is enhancing the MqGT type NME as compared to standard MGT type NME. A similar type of enhancement in MqGT type NME was also found in our earlier study for 48Ca ().
We have also decomposed the NME in terms of coupled spin-parity (Jπ) of two decaying neutrons and two created protons in the decay. Decomposed NME gives us a picture of the role of individual spin-parity on NME. The contribution of NMEs through different Jπ is shown in Figures 1–3 for different types of NME. Figure 1 examines the decomposition for MF,GT,T type matrix elements associated with light neutrino-exchange mechanism, where Figures 2, 3 examine the NME as function of Jπ for MωF,ωGT,ωT and MqF,qGT,qT type NMEs, respectively, for λ and interference mechanism. All results are presented for AV18 SRC parameterization.
FIGURE 1
FIGURE 2
FIGURE 3
For all types of NMEs, the most dominating contribution comes from 0+ states and 2+ states. The pairing effect is in play for dominating even-Jπ contributions (). The NME from 0+ and 2+ states has opposite signs and, thus, cancel the effects of each other. Other non-negligible contributions come through 4+, 3−, 5−, and 7− states.
Now we will discuss how the calculated NMEs will help to determine the bounds on Majorana neutrino mass and various lepton number violating parameters, using the lower limit on the experimental half-life of the decay. The inverse of half-life for 0νββ is given in Eq. 1. It is found that the half-life is influenced by the term CI (I = mm, mλ, λλ), lepton number violating term ην and ηλ, which are unknown, and CP-violating phase ψ. The CI are defined in Eq 7 and (8), (9), which contains mainly phase space factors and relevant NMEs. To calculate CI, we have used the improved values of phase space factors calculated in Ref. (), and for the NMEs, we have used the results of Table 2 using ISM.
The results for CI of light neutrino-exchange and λ mechanisms of 0νββ decay of 82Se and 48Ca are presented in Table 3. Here, the results for 48Ca are taken from our earlier work using the closure approximation on the λ mechanism (). It is found that values of Cmm (light neutrino-exchange) and Cλλ (λ mechanism) are similar in values, which shows the dominance of each of these mechanisms on 0νββ half-life. The interference term (Cmλ) of both the mechanisms are relatively smaller, which shows the less importance of the interference mechanism.
TABLE 3
| Quantity | 82Se | 82Se Ref. | 48Ca | 48Ca Ref. ) |
|---|---|---|---|---|
| [Years] | 2.5 × 1023 | 2.5 × 1023 | 2.0 × 1022 | 2.0 × 1022 |
| Cmm [Years]−1 | 31.21 × 10–14 | 51.3 × 10–14 | 4.06 × 10–14 | 2.33 × 10–14 |
| Cmλ [Years]−1 | 10.46 × 10–14 | −27.0 × 10–14 | 3.37 × 10–14 | −1.04 × 10–14 |
| Cλλ [Years]−1 | 36.19 × 10–14 | 150.0 × 10–14 | 5.39 × 10–14 | 10.1 × 10–14 |
| mββ [eV] | 1.83 | 1.43 | 17.92 | 23.7 |
| ηλ | 3.32 × 10–6 | 1.63 × 10–6 | 30.44 × 10–6 | 22.30 × 10–6 |
Results for half-life and bounds on neutrino mass and lepton number violating parameters. The is taken from the experimental lower limit on half-life from Ref. () for 82Se and from Ref. () for 48Ca. All results are for AV18 type SRC parameterizaionparameterization. We have assumed CP conservation (ψ =0). The results are compared with QRPA calculations for λ mechanism of Ref. ().
We have also calculated the upper bound on unknown Majorana neutrino mass (mββ) and lepton number violating parameter: the right-handed current coupling strength ηλ, using the experimental constraint on of Ref. () for 82Se and of Ref. () for 48Ca. The upper limits on mββ and ηλ are also presented in Table 3 for 82Se and 48Ca when both light neutrino-exchange and λ mechanisms co-exist. With the experimental lower limit on , the upper limits on Majorana neutrino mass (mββ) are found to be 1.83 and 17.92 eV, respectively, for 82Se and 48Ca. This difference of mββ value for 82Se and 48Ca is quite large and also found in earlier work (). With the recent progress and future prospects of new generation experiments, lower limits on will be gradually improved and thus, will improve the upper limit on mββ and also reduce the differences for different isotopes.
4 Summary
In summary, we have studied how the left-right weak boson exchange (λ) mechanism of 0νββ decay is competing with the standard light neutrino-exchange mechanism. Our interest of isotope was one of the prominent 0νββ decaying isotope 82Se. Particularly, we have calculated the NMEs for 0νββ of 82Se when both standard light neutrino-exchange and λ mechanisms co-exist. The revised formalism for nucleon currents to include the pseudoscalar term was taken care of. The nuclear shell model framework was used in the calculation, and the widely used closure approximation was adopted with suitable closure energy. Nuclear states of initial, final, and intermediate states are calculated for fpg model space with JUN45 effective shell model Hamiltonian using shell model code KSHELL. These nuclear states are used to calculate TNA, which comes in the expression of NME of 0νββ through (n − 2) decay channel. Using the calculated NMEs, we have also calculated the upper bounds on Majorana neutrino mass and lepton number violating parameters.
The results show that particularly MqGT type matrix element of λ mechanism is significantly enhanced as compared to standard MGT type NME for the inclusion of the higher-order effect of the pseudoscalar term in the nucleon current. A similar type of enhancement in MqGT type NME was also found in our earlier study for 48Ca (). The dominance of 0+ and 2+ states of neutron-neutron (proton-proton) pairs were also observed, just like earlier studies.
With the experimental lower limits on the half-life, we have used our calculated NMEs to set the upper bounds on Majorana neutrino mass (mββ). The upper limits of values of mββ are found to be 1.83 and 17.92 eV, respectively, for 82Se and 48Ca. With the new generation of experiments, the lower limit on half-life will be further improved, and thus we can expect a much more refined upper bound on mββ, which may be below 1 eV. Also, the difference for the value of mββ will be reduced.
The term CI (I = mm, mλ, λλ), which contains the phase space factors and NMEs, was also evaluated. The Cmm for light neutrino exchange and Cλλ for λ mechanism were found to be similar in values, that were larger than the term Cmλ for the interference of both the mechanisms. This shows the dominance of light neutrino exchange and the λ mechanisms over the interference mechanism. The overall dominant effect of light neutrino-exchange mechanism is observed over λ mechanism and interference of both the mechanisms for very small values of lepton number violating ηλ parameter.
In the future, it will be interesting to see the competing effect of the λ mechanism on the light neutrino-exchange mechanism and also how their contribution on 0νββ half-life will be evaluated in the current and future generation experiments.
Statements
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Author contributions
The idea of the article was originated by YI. He has also contributed to calculating the nuclear states, interpretation of results, and manuscript writing. SS is responsible for the calculation of the TNA and TBME part of the NME. He has also actively participated in preparing the manuscript. Overall, both the authors have contributed enough in to various stages of preparing the final manuscript.
Funding
YI is grateful for the funding support from JSPS KAKENHI Grant No.17K05440.
Acknowledgments
Numerical computation in this work was carried out at the Yukawa Institute Computer Facility. YI acknowledges the Tokyo Institute of Technology for allowing to use of the high-performance computing facility to perform nuclear states calculation using KSHELL. YI is also grateful to Prof. Noritaka Shimizu, CNS, the University of Tokyo, for providing the 2020 version of shell model code KSHELL.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
ArnoldR.AugierC.BakerJ.BarabashA.BroudinG.BrudaninV.et al (2005). First Results of the Search for Neutrinoless Double-Beta Decay with the NEMO 3 Detector. Phys. Rev. Lett.95, 182302. 10.1103/PhysRevLett.95.182302
2
ArnoldR.AugierC.BakalyarovA. M.BakerJ. D. (2016). Measurement of the Double-Beta Decay Half-Life and Search for the Neutrinoless Double-Beta Decay of 48Ca with the NEMO-3 Detector. Phys. Rev. D93, 112008. 10.1103/PhysRevD.93.112008
3
ArnoldR.AugierC.BarabashA. S.Basharina-FreshvilleA.BlondelS.BlotS.et al (2020). Search for the Double-Beta Decay of 82se to the Excited States of 82kr with Nemo-3. Nucl. Phys. A996, 121701. 10.1016/j.nuclphysa.2020.121701
4
AvignoneF. T.ElliottS. R.EngelJ. (2008). Double Beta Decay, Majorana Neutrinos, and Neutrino Mass. Rev. Mod. Phys.80, 481–516. 10.1103/RevModPhys.80.481
5
BareaJ.IachelloF. (2009). Neutrinoless Double-Βdecay in the Microscopic Interacting Boson Model. Phys. Rev. C79, 044301. 10.1103/PhysRevC.79.044301
6
BareaJ.KotilaJ.IachelloF. (2012). Limits on Neutrino Masses from Neutrinoless Double-β Decay. Phys. Rev. Lett.109, 042501. 10.1103/PhysRevLett.109.042501
7
BarryJ.RodejohannW. (2013). Lepton Number and Flavour Violation in Tev-Scale Left-Right Symmetric Theories with Large Left-Right Mixing. J. High Energ. Phys.2013, 1–45. 10.1007/jhep09(2013)153
8
Bhupal DevP. S.GoswamiS.MitraM. (2015). Tev-scale Left-Right Symmetry and Large Mixing Effects in Neutrinoless Double Beta Decay. Phys. Rev. D91, 113004. 10.1103/PhysRevD.91.113004
9
BrownB. A.HoroiM.Sen’kovR. A. (2014). Nuclear Structure Aspects of Neutrinoless Double-β Decay. Phys. Rev. Lett.113, 262501. 10.1103/PhysRevLett.113.262501
10
BrownB. A.RaeW. D. M. (2014). The Shell-Model Code NuShellX@MSU. Nucl. Data Sheets120, 115–118. 10.1016/j.nds.2014.07.022
11
CaurierE.MenéndezJ.NowackiF.PovesA. (2008). Influence of Pairing on the Nuclear Matrix Elements of the Neutrinoless Betabeta Decays. Phys. Rev. Lett.100, 052503. 10.1103/PhysRevLett.100.052503
12
CiriglianoV.DekensW.de VriesJ.GraesserM. L.MereghettiE.PastoreS.et al (2019). Renormalized Approach to Neutrinoless Double-β Decay. Phys. Rev. C100, 055504. 10.1103/physrevc.100.055504
13
DeppischF. F.HirschM.PäsH. (2012). Neutrinoless Double-Beta Decay and Physics beyond the Standard Model. J. Phys. G: Nucl. Part. Phys.39, 124007. 10.1088/0954-3899/39/12/124007
14
DolinskiM. J.PoonA. W. P.RodejohannW. (2019). Neutrinoless Double-Beta Decay: Status and Prospects. Annu. Rev. Nucl. Part. Sci.69, 219–251. 10.1146/annurev-nucl-101918-023407
15
EngelJ.MenéndezJ. (2017). Status and Future of Nuclear Matrix Elements for Neutrinoless Double-Beta Decay: a Review. Rep. Prog. Phys.80, 046301. 10.1088/1361-6633/aa5bc5
16
HonmaM.OtsukaT.MizusakiT.Hjorth-JensenM. (2009). New Effective Interaction Forf5pg9-Shell Nuclei. Phys. Rev. C80. 10.1103/physrevc.80.064323
17
HoroiM.NeacsuA. (2016). Analysis of Mechanisms that Could Contribute to Neutrinoless Double-Beta Decay. Phys. Rev. D93, 113014. 10.1103/PhysRevD.93.113014
18
HoroiM.NeacsuA. (2018). Shell Model Study of Using an Effective Field Theory for Disentangling Several Contributions to Neutrinoless Double-β Decay. Phys. Rev. C98, 035502. 10.1103/PhysRevC.98.035502
19
HoroiM.StoicaS. (2010). Shell Model Analysis of the Neutrinoless Double-β Decay of 48Ca. Phys. Rev. C81, 024321. 10.1103/PhysRevC.81.024321
20
IwataY.ShimizuN.OtsukaT.UtsunoY.MenéndezJ.HonmaM.et al (2016). Large-Scale Shell-Model Analysis of the Neutrinoless ββ Decay of 48Ca. Phys. Rev. Lett.116, 112502. 10.1103/PhysRevLett.116.112502
21
MenéndezJ.PovesA.CaurierE.NowackiF. (2009). Disassembling the Nuclear Matrix Elements of the Neutrinoless ββ Decay. Nucl. Phys. A818, 139–151. 10.1016/j.nuclphysa.2008.12.005
22
MohapatraR. N. (1986). New Contributions to Neutrinoless Double-Beta Decay in Supersymmetric Theories. Phys. Rev. D34, 3457–3461. 10.1103/PhysRevD.34.3457
23
MohapatraR. N.SenjanovićG. (1980). Neutrino Mass and Spontaneous Parity Nonconservation. Phys. Rev. Lett.44, 912–915. 10.1103/PhysRevLett.44.912
24
MohapatraR. N.VergadosJ. D. (1981). New Contribution to Neutrinoless Double Beta Decay in Gauge Models. Phys. Rev. Lett.47, 1713–1716. 10.1103/PhysRevLett.47.1713
25
NeacsuA.StoicaS.HoroiM. (2012). Fast, Efficient Calculations of the Two-Body Matrix Elements of the Transition Operators for Neutrinoless Double-β Decay. Phys. Rev. C86, 067304. 10.1103/physrevc.86.067304
26
PagnaniniL.AzzoliniO.BeemanJ.BelliniF.BerettaM.BiassoniM.et al (2019). “Results on Double Beta Decay of 82se with CUPID-0 Phase I,” in AIP Conference Proceedings (AIP Publishing LLC), 020019. 10.1063/1.5130980
27
PastoreS.CarlsonJ.CiriglianoV.DekensW.MereghettiE.WiringaR. B. (2018). Neutrinoless Double-β Decay Matrix Elements in Light Nuclei. Phys. Rev. C97, 014606. 10.1103/PhysRevC.97.014606
28
RathP. K.ChandraR.ChaturvediK.RainaP. K.HirschJ. G. (2010). Uncertainties in Nuclear Transition Matrix Elements for Neutrinoless ββ Decay within the Projected-Hartree-Fock-Bogoliubov Model. Phys. Rev. C82, 064310. 10.1103/PhysRevC.82.064310
29
RodejohannW. (2011). Neutrino-less Double Beta Decay and Particle Physics. Int. J. Mod. Phys. E20, 1833–1930. 10.1142/s0218301311020186
30
RodinV. A.FaesslerA.ŠimkovicF.VogelP. (2006). Assessment of Uncertainties in QRPA 0ν ββ-decay Nuclear Matrix Elements. Nucl. Phys. A766, 107–131. 10.1016/j.nuclphysa.2005.12.004
31
RodríguezT. R.Martínez-PinedoG. (2010). Energy Density Functional Study of Nuclear Matrix Elements for Neutrinoless ββ Decay. Phys. Rev. Lett.105, 252503. 10.1103/physrevlett.105.252503
32
SarkarS.IwataY.RainaP. K. (2020). Nuclear Matrix Elements for the λ Mechanism of 0ν ββ Decay of 48Ca in the Nuclear Shell-Model: Closure versus Nonclosure Approach. Phys. Rev. C102, 034317. 10.1103/PhysRevC.102.034317
33
SarkarS.KumarP.JhaK.RainaP. K. (2020). Sensitivity of Nuclear Matrix Elements of 0νββ of 48Ca to Different Components of the Two-Nucleon Interaction. Phys. Rev. C101, 014307. 10.1103/PhysRevC.101.014307
34
SchechterJ.ValleJ. W. F. (1982). Neutrinoless Double-Βdecay in SU(2)×U(1) Theories. Phys. Rev. D25, 2951–2954. 10.1103/physrevd.25.2951
35
Sen'kovR. A.HoroiM. (2013). Neutrinoless Double-β Decay of 48Ca in the Shell Model: Closure versus Nonclosure Approximation. Phys. Rev. C88, 064312. 10.1103/PhysRevC.88.064312
36
Sen’kovR. A.HoroiM.BrownB. A. (2014). Neutrinoless Double-β Decay of 82Se in the Shell Model: Beyond the Closure Approximation. Phys. Rev. C89, 054304. 10.1103/PhysRevC.89.054304
37
ShimizuN.MizusakiT.UtsunoY.TsunodaY. (2019). Thick-restart Block Lanczos Method for Large-Scale Shell-Model Calculations. Comp. Phys. Commun.244, 372–384. 10.1016/j.cpc.2019.06.011
38
ŠimkovicF.FaesslerA.MütherH.RodinV.StaufM. (2009). 0νββ-decay Nuclear Matrix Elements with Self-Consistent Short-Range Correlation. Phys. Rev. C79, 055501. 10.1103/PhysRevC.79.055501
39
ŠimkovicF.PantisG.VergadosJ. D.FaesslerA. (1999). Additional Nucleon Current Contributions to Neutrinoless Double β Decay. Phys. Rev. C60, 055502. 10.1103/PhysRevC.60.055502
40
ŠimkovicF.ŠtefánikD.DvornickỳR. (2017). The λ Mechanism of the 0νββ-Decay. Front. Phys.5, 57. 10.3389/fphy.2017.00057
41
SongL. S.YaoJ. M.RingP.MengJ. (2014). Relativistic Description of Nuclear Matrix Elements in Neutrinoless Double-β Decay. Phys. Rev. C90, 054309. 10.1103/PhysRevC.90.054309
42
ŠtefánikD.DvornickỳR.ŠimkovicF.VogelP. (2015). Reexamining the Light Neutrino Exchange Mechanism of the 0νββ Decay with Left-And Right-Handed Leptonic and Hadronic Currents. Phys. Rev. C92, 055502. 10.1103/PhysRevC.92.055502
43
TomodaT. (1991). Double Beta Decay. Rep. Prog. Phys.54, 53–126. 10.1088/0034-4885/54/1/002
44
VergadosJ. D.EjiriH.ŠimkovicF. (2012). Theory of Neutrinoless Double-Beta Decay. Rep. Prog. Phys.75, 106301. 10.1088/0034-4885/75/10/106301
45
VergadosJ. D. (1987). Neutrinoless Double β-decay without Majorana Neutrinos in Supersymmetric Theories. Phys. Lett. B184, 55–62. 10.1016/0370-2693(87)90487-4
46
VogelP. (2012). Nuclear Structure and Double Beta Decay. J. Phys. G: Nucl. Part. Phys.39, 124002. 10.1088/0954-3899/39/12/124002
47
WangX. B.HayesA. C.CarlsonJ.DongG. X.MereghettiE.PastoreS.et al (2019). Comparison between Variational Monte Carlo and Shell Model Calculations of Neutrinoless Double Beta Decay Matrix Elements in Light Nuclei. Phys. Lett. B798, 134974. 10.1016/j.physletb.2019.134974
Appendix A
One can write anti-symmetric two-body matrix elements for transition operator of 0νββ defined in Eq. 22 aswhere,
One can write in terms of 9j symbol
In terms of 6j symbol one can write
is the harmonic oscillator bracket used to convert the radial integral of neutrino potential from individual coordinate system of nucleons to relative and center of mass coordinate system of the nucleons.
Appendix B
The fα(q, r) factor of Eq. 16 can be written in terms of radial dependence, spherical Bessel function jp (qr) (p = 0, 1, 2 and 3), and FNS + HOC coupling form factors in closure approximation as ().where one can write in dipole approximation ().
μp − μn = 4.7, MV = 850 MeV, MA = 1,086 MeV mp and mπ are the mass of protons and pions (). In the present calculation, vector constant gV = 1.0 and bare axial-vector constant gA = 1.27 () are used. Both the pseudo scalar and weak magnetism terms of the nucleon currents are included in fGT,T,ωGT,ωT (q, r) factors whereas pseudo scalar term is included in fqGT,qT (q, r) factors ().
Summary
Keywords
neutrinoless double beta decay, λ mechanism, nuclear shell model, nuclear matrix element, right-handed weak boson
Citation
Iwata Y and Sarkar S (2021) Interacting Shell Model Calculations for Neutrinoless Double Beta Decay of 82Se With Left-Right Weak Boson Exchange. Front. Astron. Space Sci. 8:727880. doi: 10.3389/fspas.2021.727880
Received
20 June 2021
Accepted
30 August 2021
Published
19 November 2021
Volume
8 - 2021
Edited by
Hiroyasu Ejiri, Research Center for Nuclear Physics, Osaka University, Japan
Reviewed by
Bhupal Dev, Washington University in St. Louis, United States
J.D. Vergados, University of Ioannina, Greece
Javier Menendez, University of Barcelona, Spain
Updates
Copyright
© 2021 Iwata and Sarkar.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yoritaka Iwata, iwata_phys@08.alumni.u-tokyo.ac.jp
This article was submitted to High-Energy and Astroparticle Physics, a section of the journal Frontiers in Astronomy and Space Sciences
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.