Next Article in Journal
New Scenarios of High-Energy Particle Collisions Near Wormholes
Next Article in Special Issue
Double Beta Decay to Excited States of Daughter Nuclei
Previous Article in Journal
Jeans Instability of Dissipative Self-Gravitating Bose–Einstein Condensates with Repulsive or Attractive Self-Interaction: Application to Dark Matter
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Neutrino-Mass Sensitivity and Nuclear Matrix Element for Neutrinoless Double Beta Decay

Research Center for Nuclear Physics, Osaka University, Osaka 567-0047, Japan
Universe 2020, 6(12), 225; https://doi.org/10.3390/universe6120225
Submission received: 30 October 2020 / Revised: 23 November 2020 / Accepted: 25 November 2020 / Published: 27 November 2020
(This article belongs to the Special Issue Nuclear Issues for Neutrino Physics)

Abstract

:
Neutrinoless double beta decay (DBD) is a useful probe to study neutrino properties such as the Majorana nature, the absolute neutrino mass, the CP phase and the others, which are beyond the standard model. The nuclear matrix element (NME) for DBD is crucial to extract the neutrino properties from the experimental transition rate. The neutrino-mass sensitivity, i.e., the minimum neutrino-mass to be measured by the DBD experiment, is very sensitive to the DBD NME. Actually, the NME is one of the key elements for designing the DBD experiment. Theoretical evaluation for the DBD NME, however, is very hard. Recently experimental studies of charge-exchange nuclear and leptonic reactions have shown to be used to get single- β NMEs associated with the DBD NME. Critical discussions are made on the neutrino-mass sensitivity and the NME for the DBD neutrino-mass study and on the experimental studies of the single- β NMEs and nuclear structures associated with DBD NMEs.

1. Introduction

Neutrinoless double beta decay (DBD), which violates the lepton-number (L) conservation law by Δ L = 2, is beyond the standard electro-weak model. DBD is very powerful for studying neutrino properties such as the neutrino nature (Dirac or Majorana), the neutrino mass-scale and the mass-hierarchy, the neutrino phase, the possible right-handed weak interaction and the others, which are beyond the electro-weak standard model. Experimental and theoretical studies of DBDs have extensively been made for decades as given in review articles [1,2,3,4,5] and references therein.
The DBD transition rate is very sensitive to the DBD neutrino response (square of the DBD nuclear matrix element NME). Hereafter DBD and NME are used for the neutrinoless DBD and the neutrinoless DBD NME unless specified as two-neutrino DBD and two-neutrino DBD NME. The accurate value for the DBD NME is crucial for extracting the neutrino properties of particle physics interest from the DBD rate, if it is observed. The DBD neutrino-mass sensitivity, i.e., the minimum neutrino mass to be measured, depends on the NME. Thus one needs the NME as well to design the DBD experiment/detector. The DBD nuclear responses are discussed in review articles [6,7,8] and references therein, and the DBD NMEs are in reviews [9,10,11].
The neutrinoless DBD transition involves several transition modes such as the light neutrino-mass mode, the heavy neutrino-mass mode and others beyond the standard model [1,4]. The transition rate T 0 ν for the α -mode transition is expressed by a product of the nuclear physics factor K N ( α ) and the neutrino physics factor f ν ( α ) as
T 0 ν = K N × | f ν ( α ) | 2 , K N = g A 4 G 0 ν | M 0 ν | 2 ,
where g A = 1.27 is the axial-vector weak coupling for a free nucleon in units of the vector coupling g V , G 0 ν is the phase space volume and M 0 ν is the neutrinoless DBD NME.
The phase space factor G 0 ν is a kinematic factor, which depends on the β β energy E and the atomic number Z. It is proportional approximately to E 5 and increases as Z increases. It also depends a little on the mode α . The phase space has been well calculated [12,13,14,15].
The NME M 0 ν reflects the internal nuclear structure, and depends on the mode α of the DBD process. In case of the light neutrino-mass mode, the transition rate depends on the neutrino effective mass, and the effective mass depends on the absolute neutrino mass, the mass hierarchy and the neutrino CP phases. Thus one needs the accurate value for the NME in order to study these neutrino properties.
The DBD event-rate T 0 ν is extremely rare and the signal energy E is very low. In the case of the DBD caused by the light Majorana neutrino with the effective mass m e f f 10 meV, the rate is of the order of 10 29 per year (y). Then one needs the DBD isotopes of the order of 10 29 , i.e., N 10 tons of the DBD isotope mass to get a few signals for a 5 y exposure. The DBD signal energy E is around a few MeV, and thus is the same as the energy of typical background β γ rays. Therefore, one needs a very high-sensitivity low-background detector with the DBD isotope mass of the order of N = 10 tons.
Recent works on the neutrino oscillations show that the neutrinos are massive particles [16,17,18]. Then the neutrinoless DBD with the light neutrino-mass mode is very interesting to study the neutrino mass and the hierarchy. Actually, the DBD rate is proportional to the square of the effective neutrino mass m e f f , which depends on the CP phases and the mass hierarchy. So we discuss in the present report mainly the neutrino-mass mode DBD.
The DBD neutrino-mass sensitivity, which is defined as the minimum neutrino mass to be measured by the DBD experiment, is inversely proportional to the NME M 0 ν . The DBD NME is very sensitive to all kinds of nucleonic and non-nucleonic correlations, nuclear medium effects and effective weak-couplings. Accordingly, accurate evaluations for the NME by using existing theoretical nuclear models are indeed very hard. The DBD neutrino-mass sensitivity is discussed in the review [8] and also in a recent article [19].
Accordingly, experimental studies of the DBD NME have been made by using various kinds of nuclear and leptonic reactions. Recently charge-exchange nuclear reactions and muon capture reactions have been shown to be very useful to study the single- β NMEs and the nuclear structures associated with the DBD NME. They provide useful nuclear parameters to be used for evaluating the DBD NME.
The present report aims to discuss briefly the neutrino-mass sensitivity and the NME for the DBD experiment to access the neutrino masses of the current interest and the experimental ways to study single- β NMEs and nuclear structures associated with the DBD NME to help the theoretical evaluations for the DBD NME.
In Section 2, we discuss the DBD transition mode, the DBD NME for the neutrino-mass mode, and the DB neutrino-mass sensitivity. Experimental ways to study the single- β NMEs associated with the DBD NME and nuclear structures relevant to the DBD are discussed in Section 3, and the quenching of the axial-vector weak coupling is discussed in Section 4. Finally, concluding remarks and discussions are given in Section 5.

2. DBD Transition Mode, DBD Neutrino-Mass Sensitivity and DBD NME

2.1. Neutrinoless DBD and Effective Neutrino Mass

The neutrinoless DBD, which requires the Majorana nature (neutrino = anti-neutrino) of the neutrino, is due to several modes ( α ), the light neutrino-mass mode, the heavy neutrino-mass mode, the SYSY particle mode, the right-handed weak interaction mode, the Majoron-boson mode, and others. They are discussed as given in reviews [20,21,22] and references therein, and also in recent reviews [1,3,4,5] and references therein. The neutrinoless DBD is indeed a sensitive probe to study these neutrino properties and weak interactions beyond the standard model. The neutrinoless DBD modes beyond the standard model, together the two neutrino DBD within the standard model, are schematically shown in Figure 1.
In the case of the light neutrino-mass mode, a light Majorana neutrino is exchanged between the 2 neutrons in the DBD nucleus. The transition rate T 0 ν is expressed by using the effective neutrino mass m e f f [1,4,8] as
T 0 ν = g A 4 G 0 ν | M 0 ν | 2 | m e f f | 2 ,
where T 0 ν is written by using the halflife t 1 / 2 as T 0 ν =ln2/ t 1 / 2 , G 0 ν is the phase space volume for the light mass mode, M 0 ν is the NME for the light neutrino-mass mode, and m e f f is the effective neutrino mass, respectively. The effective neutrino mass is expressed as
m e f f = | i | U e i | 2 m i e i α i |
where U e i is the mixing coefficient, m i is the ith eigen-mass and α i is the phase, respectively. In cases of the normal mass-hierarchy (NH) of m 1 m 2 m 3 and the inverted mass-hierarchy (IH) of m 3 m 1 m 2 , the effective masses are given approximately as
| m e f f | | C 13 2 S 12 2 ( Δ m s 2 ) + S 13 2 ( Δ m a 2 ) e 2 i α 2 | NH ,
| m e f f | | C 13 2 ( Δ m a 2 ) [ 1 sin 2 2 θ 12 sin 2 α 12 ] 1 / 2 | IH ,
where C i j and S i j are the mixing coefficients, θ i j is the mixing angle, α 12 = α 2 α 1 is the Majorana phase difference, Δ m s 2 is the solar neutrino mass-square difference, and Δ m a 2 is the atmospheric neutrino mass-square difference [1,4,5,7,8].
The mixing parameters, the mixing angles and the mass-square differences in Equations (4) and (5) are known from the neutrino oscillation experiments. Then the effective mass to be studied is given as a function of the minimum neutrino mass m 0 = m 1 and the Majorana phase α 2 in the case of NH and that of the minimum mass m 0 = m 3 and the Majorana phase difference α 2 α 1 in the case of IH, respectively. The regions of the effective neutrino-masses are shown as a function of m 0 in Figure 2. Here the constrains from the cosmological bound of m 0 9 meV for IH and m 0 26 meV for NH are also shown in Figure 2.
The upper and lower bounds correspond to the effective masses for α 2 = 0 and α 2 = π /2 in the case of NH, and for α 12 = 0 and α 12 = π /2 in the case of IH. Thus the experimental value for the effective mass depends on the mass hierarchy, the minimum mass and the Majorana phases.
Then one gets m e f f 50 18 meV, depending on α 12 = 0 − π /2, and m 0 9 meV in the case of the IH, and m e f f 15 meV and m 0 26 meV in the case of the NH. The effective mass in the case of IH is plotted as a function of α 12 in Figure 3. If the measured effective mass is larger than 18 meV, one concludes the IH mass spectrum and gets the Majorana phase α 12 . If that is smaller than 4 meV, one concludes the NH mass spectrum and gets constrains on the minimum mass m 0 and the Majorana phase of α 2 .

2.2. DBD Neutrino-Mass Sensitivity and DBD Nuclear Matrix Element

The DBD transition rate is rewritten for practical use in terms of the number of the DBD signals per ton (t) of the DBD isotopes per year (y). It is expressed by using the nuclear sensitivity m n as [1,8]
T 0 ν = ( m e f f m n ) 2 , m n = k M 0 ν , k = 7.8 meV A 1 / 2 g A 2 ( G 0 ν ) 1 / 2 ,
where T 0 ν is the transition rate in units of per t y, k is the nuclear kinematic factor defined as the nuclear sensitivity m n in the case of M 0 ν = 1, and G 0 ν is the phase space volume in units of 10 14 /y, and A is the mass number of the DBD nucleus. m n corresponds to the effective neutrino mass to give the DBD rate of 1/t y. The factor k reflects the external kinematic factor and 1 / M 0 ν does the internal nuclear-structure one, and the nuclear sensitivity m n is given by the product of them.
The nuclear kinematic factor k and the nuclear sensitivity m n for typical DBD nuclei of current interests are shown in Table 1. It is interesting to note that the kinematic factors k for most nuclei are around 40 meV near the upper bound of the IH mass, while the values for 76 Ge and 150 Nd are, respectively, a factor of two larger and smaller than 40 meV because of the smaller and larger phase space factors. The m n = 40 meV sensitivity means one DBD signal per year per ton of the DBD isotopes in the case of M 0 ν = 1 and m e f f = 40 meV, and four signals in the case of M 0 ν = 2.
The DBD neutrino-mass sensitivity m m is derived from the condition that the number of the DBD signals is well above the fluctuation of the number of background (BG) signals. Then one gets
T 0 ν N T η ϵ δ 0 ( B N T ) 1 / 2 ,
where N is the DBD isotope mass in units of ton (t), η is the enrichment, ϵ , is the signal efficiency, T is the exposure time in units of year (y), B is the BG rate per t y, and δ 0 is a constant around two in the case of the 90% confidence level for the DBD signal identification. The DBD neutrino-mass sensitivity is given as [1,8]
m m = m n d , d = d 0 × η 1 / 2 ϵ 1 / 2 ( N T / B ) 1 / 4 ,
where d is the detector sensitivity and d 0 is around δ 0 1.4.
The DBD neutrino-mass sensitivity depends on the nuclear sensitivity m n = k / M 0 ν , and thus linearly on M 0 ν , but weakly on N and B by the power of 1/4. Thus the 30 % less M 0 ν is equivalent to an order of magnitude more isotope mass N or less BG rate B. Thus one needs an accurate value for the NME M 0 ν to design the DBD detector with the given neutrino-mass sensitivity.
The DBD NME is expressed as [7,8]
M 0 ν = ( g A e f f g A ) 2 [ M G T 0 ν + M T 0 ν ] + ( g V e f f g A ) 2 M F 0 ν ,
where M G T 0 ν , M T 0 ν and M F 0 ν are the Gamow–Teller (GT), Fermi (F) and tensor (T) NMEs, respectively, and g A e f f and g V e f f are the effective axial-vector and vector couplings in units of the vector coupling g V for a free nucleon. The GT, T and F NMEs are given as
M G T 0 ν = i < t ± σ h G T ( r 12 , E i ) σ t ± > ,
M T 0 ν = i < t ± h T ( r 12 , E i ) S 12 t ± > ,
M F 0 ν = i < t ± h F ( r 12 , E i ) t ± > ,
where the sum is for all the intermediate states (i), t and σ are the isospin and spin operators, and h α is the neutrino potential with α = GT, T, F, and r 12 is the distance between the two neutrons involved in the virtual-neutrino exchange.
The NMEs of M G T 0 ν , M T 0 ν and M F 0 ν are evaluated as given in Equations (10)–(12) by using appropriate theoretical models. Then the g A e f f and g V e f f stand for the quenching (re-normalization) effects due to such nucleonic and non-nucleonic correlations and nuclear medium effects that are not explicitly included in the model.
The neutrino potential is approximately given in terms of 1/ r 12 p with p being the virtual-neutrino momentum of the order of 50–150 MeV/c. Then the angular momentum is around l = (1–4) . Accordingly, the GT and T NMEs are very sensitive to nuclear spin isospin and momentum correlations and the F NME to the isospin and momentum correlations. In the case of a typical DBD run for T = 5 y by using a typical detector with N = 1 t, B = 1/t y, η = 0.9 and ϵ = 0.6, one gets the detector sensitivity d 1.5 and the neutrino mass sensitivities m m 1.5 m n .
The DBD neutrino-mass sensitivities m m in the cases of the NMEs of M 0 ν = 1, 2 and the BG rates of B = 1/t y and 0.01 /t y are shown as a function of the exposure of N T in Figure 4. Then, one needs exposures of N T 3, 16 and 256 t y in order to access the effective mass of m m = 20 meV in cases of the NMEs of M 0 ν = 3, 2 and 1, respectively. One needs almost a 1000 ton year exposure to access the NM mass in the case of M 0 ν = 2 and B = 0.01/t y. The NME is indeed a key parameter for the neutrino mass study by DBD experiments.
Experimentally one gets the value for or the upper limit on M 0 ν × m e f f from the measured value for or the upper limit on the DBD rate, i.e., the value for or the lower limit on the halflife. The present limits on the DBD halflives are discussed in the recent conference report [23]. Quite stringent limits on the DBD rates, ln2/ t 1 / 2 , have been obtained for 76 Ge [24,25], 130 Te [26], 136 Xe [27,28] and others. Figure 5 shows the allowed regions for M 0 ν × m e f f in the case of the half life limit of T 1 / 2 10 26 y for 76 Ge and 136 Xe, which is close to the present limits [25,28]. Therefore, experimental groups are trying to improve their detector sensitivities to access the IH neutrino-mass by increasing the isotope mass N and/or lowering the BG rate B. The improvement factor does depend much on the NME.
In fact, the DBD nucleus of the current interest is a medium-heavy nucleus with the mass number A 100. It is a complex many-body hadron (nucleon, meson, and isobar) system, and thus the DBD NME is very sensitive to all kinds of nucleonic and non-nucleonic correlations, nuclear deformations and nuclear medium effects. Therefore, accurate theoretical evaluations for the DBD NME M 0 ν , including g A e f f and g V e f f , are in practice very hard. They depend much on the models and the interaction parameters used for the models. Accordingly, the evaluated values for the DBD NMEs scatter over an order of magnitude, depending on the model, the g A e f f and other parameters used in the model [7,8]. Here it is noted that the current neutrino-mass limits are derived from the observed limits on the DBD halflives by using the existing minimum and maximum theoretical NMEs, which are sensitive to the models and parameters used in the calculations.

3. Experimental Approaches to DBD NMEs

DBD NMEs are very important for DBD neutrino studies, but accurate theoretical evaluations for them are very hard, and there are no direct ways to measure accurately the DBD NME M 0 ν . Recently various experimental studies have been made on single β -NMEs and nuclear structures associated with the DBD NME to provide useful information to help the theoretical evaluations for the DBD NMEs.

3.1. Experimental Studies of Two-Neutrino DBD NMEs

Here we note the 2-neutrino DBD (2 ν β β ) NMEs for the transition of Z A X Z + 2 A X, where A is the mass number of the DBD nucleus and Z and Z + 2 are the atomic numbers of the initial and final nuclei, respectively. The transition rate T 2 ν is expressed as
T 2 ν = g A 4 G 2 ν | M 2 ν | 2 ,
where G 2 ν is the phase space volume, M 2 ν is the 2-neutrino NME. The NME is given by
M 2 ν = i M G T ( i ) × M G T + ( i ) Δ ( i ) ,
where M G T ( i ) and M G T + ( i ) are, respectively, the single- β NME for Z A X Z + 1 A X and the single- β + NME for Z + 2 A X Z + 1 A X via the ith state in the intermediate nucleus, and Δ ( i ) is the energy denominator expressed as Δ ( i ) =Q i +Q β β /2 with Q i and Q β β being the EC (e-capture) Q value for the ith intermediate state and the DBD Q value.
Actually, the 2-neutrino DBD NMEs have been extensively studied experimentally and theoretically, as given in recent reviews [8,11]. Some of the theoretical works are QRPA model calculations [29,30,31], shell-model calculations [32,33,34,35], IBM-model calculations [36], and effective field theory calculations [37]. Here we show the 2-neutrino DBD NMEs evaluated empirically by using experimental single- β NMEs for low-lying (fermi-surface) quasi-particle (FSQP) states [38]. Then the single- β ± NMEs are expressed as
M G T ± ( i ) = k Q ± M G T Q P ± ( i ) ,
where M G T Q P ± ( i ) is the single- β ± GT NME for the ith FSQP, and k Q ± = ( g A Q P e f f / g A ) ± stands for the quenching (re-normalization) coefficient for the β ± due to the nucleonic and non-nucleonic correlations and nuclear medium effects, which are not included in the simple QP model. The quenching coefficient k Q ± is derived as the ratio of the experimental M G T ± ( i ) and the FSQP model M G T Q P ± ( i ) by assuming the same phase for both β + and β NMEs.
The NMEs M 2 ν for DBD nuclei of the current interests are derived by using the FSQP model with the experimental g A Q P e f f / g A and Δ ( i ) [4,8,38]. They reproduce well the NME M 2 ν derived from the observed rate T 2 ν as shown in Figure 6.
On the other hand, the neutrinoless DBD NME M 0 ν includes the neutrino potential h k ( r 12 ) to rink the single- β ± NMEs, as expressed in Equations (10)–(12). Among the three DBD modes, the GT and F mode NMEs of M G T 0 ν and M F 0 ν are major components. The axial-vector single- β ± NMEs and the vector single- β ± NMEs are then associated, respectively, with the M G T 0 ν and M F 0 ν via the neutrino potential. Accordingly, the single- β ± NMEs are used to help theoretical evaluations for M G T 0 ν and M F 0 ν .
Actually, single- β ± NNEs and nuclear structures associated with the DBD NMES are experimentally studied by means of nuclear charge-exchange reactions (CERs), leptonic (muon and neutrino) CERs and photo-nuclear reactions.

3.2. Nuclear Charge-Exchange Reactions

Nuclear CERs are used to study single- β ± NMEs in wide energy and momentum regions relevant to DBD. The high energy-resolution ( 3 He,t) CER at RCNP Osaka is a powerful probe to study axial-vector NMEs for individual states in DBD nuclei. Axial-vector states are preferentially excited by using the 0.42 GeV 3 He beam. GT NMEs with J π = 1 + have been measured, as shown in Figure 7, and used to evaluate the two-neutrino DBD NMEs [38,39,40,41].
The spin-dipole (SD) NME with J π = 2 is one of the major components of the neutrinoless DBD NMEs. Recently the SD NMEs were derived for the first time from the CER cross sections [42,43]. They are found to be consistent with the FSQP model NMEs and the NMEs derived from the β -decay f t values, as shown in Figure 8. Here it is noted that the momentum dependence of the NMEs was studied by measuring the cross sections in a wide angular range of θ = 0–4 deg, corresponding to the momentum transfer of p = 20–120 MeV/c relevant to the DBD momentum region [44]. The observed GT and SD responses (squares of the GT and SD NMEs) are found to be constant over the wide momentum range, and thus one can use both the responses for the single- β decays and the two-neutrino DBDs at the low momentum region and those for the CERs at the medium momentum region for evaluating the neutrinoless DBD responses [43,44].
The CER energy spectra show weak discrete peaks at the low energy region and the strong GT and SD giant resonances at the medium energy region, as shown in Figure 7. Note that the GT and SD strengths for low-lying states are much reduced with respect to the single quasi-particle strengths due to the destructive interference with the strong GT and SD GRs [6,7,8].

3.3. Leptonic (Muon and Neutrino) Charge-Exchange Reactions

Low-momentum negative muons stopped in a target are trapped into the lowest orbit of the target atom. Then, they are captured into the nucleus via the leptonic CER of μ + Z + 2 A X ν μ + Z + 1 A X * . This is analogous to the anti-neutrino nuclear interaction of ν ¯ μ + Z + 2 A X μ + + Z + 1 A X * . Muon capture reactions and nuclear structures are reviewed in [45]. The μ CER transfers the medium energy and medium momentum relevant to DBD. Thus the muon CER is used to study the NMEs relevant to DBD in the same energy and momentum regions. The muon CER is schematically shown in Figure 9.
The muon-capture reaction on 100 Mo shows a giant resonance around E 12 MeV [46], as shown in Figure 9. The observed GR is reproduced by the pnQRPA calculation [47]. The pnQRPA calculations for the capture rates for 100 Mo and other medium-heavy nuclei are much smaller than the empirical rates [48], suggesting some quenching coefficient of g A e f f / g A 0.4 [47,49], as in the ν NMEs studied by the nuclear β ± decays [50,51]. On the other hand the QRPA calculations for medium heavy nuclei are in accordance with the empirical rates, suggesting no quenching for g A e f f [52]. Since the capture rate is sensitive to the strength distribution as a function of the excitation energy, it is interesting to compare the theoretical and experimental distributions.
DBD NMEs are studied directly by using intense electron-neutrino beams obtained from high-intensity GeV-proton accelerators as the ORNL 1 GeV SNS and the J-PARC 3 GeV booster-synchrotron [53,54]. The neutrino energy spectrum is shown in Figure 10.
The neutrino energy of around 10–45 MeV is just appropriate for studying the NMEs for the medium energy neutrinos relevant to DBD. The expected neutrino intensity is of the order of 10 15 per sec. Since the neutrino cross-section is of the order of 10 40 cm 2 , one needs multi-ton scale DBD isotopes as used for DBD experiments.

3.4. Photo-Nuclear Reactions

Electromagnetic (EM) interactions are used to study DBD NMEs. The charged-current weak NME is derived from the isovector component of the EM NME. The electric and magnetic EM NMEs correspond to the axial-vector and vector weak NMEs, respectively. The special case is the EM NME for the isobaric analog state (IAS), as shown in Figure 11. The EM and weak NMEs are related as [8,55,56]
< f | g W m β | i > = g W g E M ( 2 T ) 1 / 2 < f | g E M m γ | I A S > ,
where g W and g E M are the weak and EM couplings, respectively, T is the isospin of the initial state and m β is the weak transition operator analogous to the EM one m γ . The β -NME analogous to the E1- γ NME was obtained for 141 Ce. The obtained NME is found to be quenched by a factor of 0.3 with respect to the QP NME [8,55,57].
NMEs for |f> →|i> are studied by measuring neutrons following photo-nuclear reactions of ( γ , n) on |f> via IAS of |i> as illustrated in Figure 11. The photo-nuclear reaction cross-section is used to get the absolute value for the NME, and the angular correlation of the ( γ , n) from the polarized photon is used to get the multi-polarity (E1,M1 or E2) of the NME [56]. Here the polarized photon is obtained from the polarized laser photon scattered off the GeV electron (Figure 12).

4. Quenching Coefficient for Weak Coupling and DBD NME

Axial-vector weak transition rates for low-lying states are quenched (reduced) with respect to simple quasi-particle (QP) model evaluations, as discussed since the 1960s [6,7]. Similarly, magnetic gamma-transition rates and magnetic moments for low lying states are also quenched (reduced) with respect to QP model evaluations. These are due to the nucleonic and non-nucleonic spin-isospin correlations and others, as discussed in review articles [6,8,57,58,59].
The quenching coefficient is conventionally expressed in terms of the effective axial-vector coupling g A e f f / g A . It is derived from the experimental and theoretical single- β NMEs as
g A e f f g A = M EXP M MODEL ,
where M EXP and M MODEL are the experimental and theoretical (model) NMEs, respectively. Here we use the pnQRPA model NME M QR , which includes the spin-isospin nucleonic correlations.
The quenching coefficients to be discussed are for medium-heavy nuclei of the DBD interest. The coefficients for axial vector GT and SD NMEs derived from the low-momentum (1–3 MeV/c) β -decay rates are around g A e f f / g A 0.5 [50,51] as shown in Figure 13. The ( 3 He,t) CERs show that the GT and SD NMEs remain constant over the wide momentum region of p = 30–100 MeV/c relevant to the DBD momentum [44], and the coefficients for GT and SD transitions are around 0.5. The muon CER covers the similar momentum region of p 80 MeV/c. The muon capture rates on 100 Mo and other DBD nuclei show the similar quenching coefficient around 0.4–0.5 [47,49], but the recent QRPA calculations suggest no quenching [52]. The magnetic hexadecapole (HD) NMEs, which are mainly isovector component, are also quenched by a coefficient around 0.33 with respect to the multi quasi-particle phonon (MQPP) model [60].
In the present work, we discussed mainly the quenching coefficient with respect to QRPA. The quenching of the axial-vector NMEs have been discussed in other model calculations as discussed in [8]. Among them, IBM-2 [61] gives similar quenching coefficients as in the present work. It is of great interest to study experimentally and theoretically how the quenching coefficients are universal or depend on the momentum and multi-polarity in medium heavy nuclei. It is noted that the quenching coefficient g A e f f / g A reflects some nuclear core effects such as non-nucleonic (isobar) and nuclear medium effects, which manifest in medium and heavy nuclei.
The neutrinoless DBD NMEs for 130 Te are plotted as a function of the quenching coefficient g A e f f / g A in Figure 14. Here we use the GT, F, and F NMEs evaluated by pnQRPA [62] and the quenching coefficient g V e f f / g V = 1 for the vector coupling. The GT and T NMEs decrease as the quenching coefficient decreases, while the F NME stays constant. In the case of the experimental coefficient of g A e f f / g A = 0.5, one gets M 0 ν = 1.8. Similarly, one gets M 0 ν = 2.4 for 76 Ge. These NMEs are based on the pnQRPA NMEs [62] and the experimental g A e f f , and thus include uncertainties, which should be further investigated.
Here it is noted that 2-neutrino DBD NMEs (particularly M + in Equation (14)) are sensitive to the particle-particle coupling g p p because the low-momentum GT proton→neutron transitions are blocked by the neutron excess [51]. However, neutrinoless DBD NMEs involve mainly medium-momentum and medium-multipole transitions and thus the NMEs are not sensitive to g p p .

5. Concluding Remarks and Discussions

The neutrinoless DBD is a very sensitive and powerful probe for studying neutrino properties of particle physics interests beyond the standard model. Extensive experimental and theoretical studies of the DBDs are under progress to study the neutrino masses in IH and NH. The neutrino-mass sensitivity m m depends on the nuclear sensitivity m n and the detector sensitivity d, and the nuclear sensitivity depends on the phase space factor and the nuclear matrix element NME. Actually, the m m is very sensitive to the NME, which is one of the key elements of the high-sensitivity DBD experiments to search for the very small ν mass. There are, however, no direct experimental ways to measure the DBD NME, and accurate theoretical evaluations for the NME, including the weak coupling constant g A , are very hard.
The high energy-resolutions nuclear CER of ( 3 He,t) and the leptonic CER of ( μ , ν μ ) at RCNP Osaka have recently been shown to be useful for experimental studies of the DBD NMEs. These reaction data provide the single- β NMEs and the nuclear structures associated with the DBD NMEs.
The single- β NMEs in the wide momentum region of the DBD interest are found to be quenched uniformly by the coefficient g A e f f / g A 0.5 with respect to the pnQRPA NMEs due to such non-nucleonic correlations and nuclear medium effects that are not explicitly included in the pnQRPA. Using this quenching coefficient, the DBD NMEs are evaluated as M 0 ν 2.4 and 1.8 for the DBD nuclei of 76 Ge and 130 Te, respectively. The halflives in the case of the IH mass of m e f f = 20 meV are 1.5 × 10 28 y and 0.4 × 10 28 y for 76 Ge and 130 Te, respectively.
The single- β NMEs relevant to the DBD NMEs are studied also by neutrino nuclear reactions using neutrino beams and photo nuclear reactions using polarized photon beams as well. Nuclear structures and NMEs relevant to the DBD NMEs may be studied by double charge-exchange reactions using light and heavy ion beams [63] at INFN-LNS Catania, RCNP Osaka, RIKEN Wakoh and others as discussed in [8].
The quenching problem has been discussed mainly on axial-vector NMEs, and thus the quenching coefficient is discussed in terms of the effective (re-normalized) axial-vector weak coupling g A e f f . In fact, medium-heavy DBD nuclei are a very complex many-body hadron (nucleon/mason/isobar) system, but most models are such simplified nucleon-based models that do not include non-nucleonic correlations, higher-shell orbits and/or nuclear medium effects. Accordingly, these effects are incorporated into the g A e f f . In case of the pnQRPA with nucleonic correlations, the g A e f f stands for the non-nucleonic correlations and nuclear medium effects, which may reduce the axial-vector NMEs.
One of the important non-nucleonic correlations are the isobar nucleon-hole contribution Δ N 1 , which is a kind of the quark spin-isospin polarization of the nuclear medium. Experimentally one may test the non-nucleonic effect by comparing the summed axial-vector strength with the nucleon-based sum-rule limit as discussed in [7,8].
Actually, the DBD strength for the grand-state transition is only a small fraction of the order of 10 5 of the total DBD strength summed over all final states up around to 60 MeV, and thus accurate evaluations for them are interesting by taking into accounts exactly all possible effects such as non-nucleonic (meson, isobars), nuclear medium effects, momentum-dependent interactions, and others that affect the DBD NME. Here it is interesting to find that the recent calculation, which includes two-body currents, reproduces the large GT rates in simple double-magic nuclei [64].
The vector weak coupling g V may be also modified due to non-nucleonic correlations and nuclear medium effects, which are not explicitly included in the models. In fact, nuclear models used to evaluate vector-type DBD NMEs are not perfect, and thus NMEs M F 0 ν depend more or less on the models and the interaction parameters used in the models. So it is interesting to study experimentally vector-type forbidden single- β NMEs and to compare them with theoretical NMEs to see if one needs some re-normalization coefficient g V e f f / g V defined as the ratio of the experimental to theoretical vector-type NMEs [6]. Experimentally, extraction of the vector NMEs in non-unique forbidden β -decays, however, is hard. One way is to study analogous γ transitions [55].
The nuclear re-normalization effect is a small fraction of the order of 10 6 of the total DBD strength, and thus not significant in the case of the super-allowed Fermi transition for the isobaric analog state with the large strength near the sum rule limit [65], but may not be negligible in the ground-state transitions associated with the isospin, spin, and nuclear-medium changes and the tiny strength of the order of 10 5 of the total strength.
Perspectives of experimental studies of single- β NMEs and nuclear structures associated with the DBD NMEs are discussed elsewhere [66].

Funding

This research received no external funding.

Acknowledgments

The author thanks C. Agodi for valuable discussions and encouragements during the course of the present work, and the nuclear and leptonic CER colleagues for the collaboration works.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Ejiri, H. Double beta decays and neutrino masses. J. Phys. Soc. Jpn. 2005, 74, 2101–2127. [Google Scholar] [CrossRef] [Green Version]
  2. Avignone, F.; Elliott, S.; Engel, J. Double beta decay, Majorana neutrino, and neutrino mass. Rev. Mod. Phys. 2008, 80, 481–516. [Google Scholar] [CrossRef] [Green Version]
  3. Ejiri, H. Double β-decays and neutrino nuclear responses. Prog. Part. Nucl. Phys. 2010, 54, 249–257. [Google Scholar] [CrossRef]
  4. Vergados, J.; Ejiri, H.; Šimkovic, F. Theory of neutrinoless double-β-decay. Rep. Prog. Phys. 2012, 75, 106301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Vergados, J.; Ejiri, H.; Šimkovic, F. Neutrinoless double β-decay and neutrino mass. Int. J. Mod. Phys. E 2016, 25, 1630007. [Google Scholar] [CrossRef] [Green Version]
  6. Ejiri, H.; Fujita, J.I. Effective coupling constants for beta and gamma transitions in medium heavy nuclei. Phys. Rep. 1978, 38, 85–131. [Google Scholar] [CrossRef]
  7. Ejiri, H. Nuclear spin isospin responses for low-energy neutrinos. Phys. Rep. 2000, 338, 265–351. [Google Scholar] [CrossRef]
  8. Ejiri, H.; Suhonen, J.; Zuber, K. Neutrino nuclear responses for astro-neutrinos, single β-decays, and double β-decays. Phys. Rep. 2019, 797, 1–102. [Google Scholar] [CrossRef]
  9. Suhonen, J.; Civitarese, O. Weak interaction and nuclear structure aspect of nuclear double beta decay. Phys. Rep. 1998, 300, 123–214. [Google Scholar] [CrossRef]
  10. Suhonen, J.; Civitarese, O. Double-beta decay nuclear matrix elements in the pnQRPA framework. J. Phys. G Nucl. Part. Phys. 2012, 39, 035105. [Google Scholar] [CrossRef]
  11. Engel, J.; Menéndez, J. Status and future of nuclear matrix elements for neutrinoless double β-decay: A review. Rep. Prog. Phys. 2017, 60, 046301. [Google Scholar] [CrossRef] [PubMed]
  12. Kolita, J.; Iachello, F. Phase space factors for double-β decay. Phys. Rev. 2012, 85, 034316. [Google Scholar]
  13. Kolita, J.; Iachello, F. Phase space factors for β+β+ decay and competing models of double-β decay. Phys. Rev. 2013, 87, 024313. [Google Scholar]
  14. Stoica, S.; Mirea, M. New calculations for phase space factors involved in double-β decay. Phys. Rev. 2013, 88, 037303. [Google Scholar]
  15. Stoica, S.; Mirea, M. Phase space factors for double-beta decays. Front. Phys. 2019, 17, 12. [Google Scholar] [CrossRef]
  16. Fukuda, Y.; Hayakawa, T.; Ichihara, E.; Inoue, K.; Ishihara, K.; Ishino, H.; Itow, Y.; Kajita, T.; Kameda, J.; Kasuga, S.; et al. (Super Kamiokande collaboration). Evidence for oscillation of atmospheric neutrinos. Phys. Rev. Lett. 1998, 81, 1562–1567. [Google Scholar] [CrossRef] [Green Version]
  17. Ahmad, Q.R.; Observ, W.; Washington, U.; Allen, R.C.; Andersen, T.C.; Anglin, J.D.; Buhler, G.; Barton, J.C.; Beier, J.W.; Bercovitch, M.; et al. Measurement of the rate of νe+d→p+p+e- interactions produced by 8B solar neutrinos at the Sudbury Neutrino Observatory. Phys. Rev. Lett. 2001, 87, 071301. [Google Scholar] [CrossRef] [Green Version]
  18. Eguchi, K.; Enomoto, S.; Furuno, K.; Goldman, J.; Hanada, H.; Ikeda, H.; Ikeda, K.; Inoue, K.; Ishihara, K.; Itoh, W.; et al. (KamLAND collaboration). First results from KamLAND: Evidence for reactor antineutrino disappearance. Phys. Rev. Lett. 2003, 90, 021802. [Google Scholar] [CrossRef] [Green Version]
  19. Ejiri, H. Nuclear and detector sensitivities for neutrinoless double β-decay experiments. Adv. High Energy Phys. 2020. submitted. [Google Scholar]
  20. Doi, M.; Kotani, T.; Takasugi, E. Double beta decay and Majorana neutrino. Prog. Theor. Phys. 1985, 83, 1–175. [Google Scholar] [CrossRef] [Green Version]
  21. Faessler, A.; Simkovic, F. Double beta decay. J. Phys. G Nucl. Part. Phys. 1998, 24, 2139–2178. [Google Scholar] [CrossRef]
  22. Vergados, J.D. The neutrinoless double beta decay from a modern perspective. Phys. Rep. 2002, 361, 1–56. [Google Scholar] [CrossRef] [Green Version]
  23. Detwiler, J. Future neutrino-less double-beta decay experiments. Proc. Neutrino 2020. [Google Scholar] [CrossRef] [Green Version]
  24. Alvis, S.I.; Arnquist, I.J.; Avignone, F.T.; Barabash, A.S.; Barton, C.J.; Basu, V.; Bertrand, F.E.; Bos, B.; Busch, M.; Buuck, M.; et al. Search for neutrinoless double-β decay in 76Ge with 26 kgy of exposure from the Majorana demonstrator. Phys. Rev. C 2019, 100, 025501. [Google Scholar] [CrossRef] [Green Version]
  25. Agostini, M.; Bakalyarov, A.M.; Balata, M.; Barabanov, I.; Baudis, L.; Bauer, C.; Bellotti, E.; Belogurov, S.; Bettini, A.; Bezrukov, L.; et al. Improved limit on neutrinoless double β-decay in 76Ge from GERDA Phase II. Phys. Rev. Lett. 2018, 120, 132503. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Adams, D.Q.; Alduino, C.; Alfonso, K.; Avignone, F.T.; Azzolini, O.; Bari, G.; Bellini, F.; Benato, G.; Biassoni, M.; Branca, A.; et al. Improved limit on neutrinoless double-beta decay in 130Te with CUORE. Phys. Rev. Lett. 2019, 124, 122501. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Gando, A. Search for Majorana neutrinos near the inverted mass hierarchy region with KamLAND-Zen. Phys. Rev. Lett. 2016, 117, 109903. [Google Scholar] [CrossRef] [Green Version]
  28. Anton, G.; Badhrees, I.; Barbeau, P.S.; Beck, D.; Belov, V.; Bhatta, T.; Breidenbach, M.; Brunner, T.; Cao, G.F.; Cen, W.R.; et al. Search for neutrinoless double-β decay with the complete EXO-200 data set. Phys. Rev. Lett. 2019, 123, 161802. [Google Scholar] [CrossRef] [Green Version]
  29. Pirinen, P.; Suhonen, J. Systematic approach to β and 2νββ decays of mass A = 100–136 nuclei. Phys. Rev. 2015, 91, 054309. [Google Scholar]
  30. Simkovic, F.; Dvornicky, R.; Stefanik, D.; Faessler, A. Improved description of the 2νββ-decay and possibility to determine the effective axial-vector coupling constant. Phys. Rev. 2018, 97, 034315. [Google Scholar] [CrossRef] [Green Version]
  31. Gando, A.; Gando, Y.; Hachiya, T.; Minh, M.H.; Hayashida, S.; Honda, Y.; Hosokawa, K.; Ikeda, H.; Inoue, K.; Ishidoshiro, K.; et al. Precision analysis of the 136Xe two-neutrino ββ spectrum in KamLAND-Zen and its impact on the quenching of nuclear matrix elements. Phys. Rev. Lett. 2019, 122, 192501. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Caurier, E.; Nowacki, F.; Poves, A. Shell model description of the ββ decay of 136Xe. Phys. Lett. 2012, 711, 62–64. [Google Scholar] [CrossRef] [Green Version]
  33. Coraggio, L.; De Angelis, L.; Fukui, T.; Gargano, A.; Itaco, N.; Nowacki, F. Renormalization of the Gamow-Teller operator within the realistic shell model. Phys. Rev. 2019, 100, 014316. [Google Scholar] [CrossRef] [Green Version]
  34. Coello Perez, E.A.; Menendez, J.; Schwenk, A. Two-neutrino double electron capture on 136Xe based on effective theory and the nuclear shell model. Phys. Lett. 2019, 797, 134885. [Google Scholar] [CrossRef]
  35. Kostensalo, J.; Suhonen, J. Consistent large-scale shell-model analysis of the two-neutrino ββ and single β bramchings in 48Ca and 96Zr. Phys. Lett. 2020, 802, 135192. [Google Scholar] [CrossRef]
  36. Yoshida, N.; Iachello, F. Two-neutrino double-β decay in the interacting boson fermion model. Prog. Theor. Exp. Phys. 2013, 2013, 043D01. [Google Scholar] [CrossRef] [Green Version]
  37. Coello Perez, E.A.; Menendez, J.; Schwenk, A. Gamow-teller and doubleβ decays of heavy nuclei with an effective theory. Phys. Rev. 2018, 98, 045501. [Google Scholar]
  38. Ejiri, H. Fermi surface quasi particle model nuclear matrix elements for two neutrino double beta decays. J. Phys. Nucl. Part. Phys. 2017, 44, 15201. [Google Scholar] [CrossRef] [Green Version]
  39. Thies, J.H.; Frekers, D.; Adachi, T.; Dozono, M.; Ejiri, H.; Fujita, H.; Fujiwara, M.; Grewe, E.W.; Hatanaka, K.; Heinrichs, P.; et al. The (3He,t) reaction on 76Ge, and double-β decay matrix element. Phys. Rev. 2012, 86, 014304. [Google Scholar]
  40. Thies, J.H.; Puppe, P.; Adachi, T.; Dozono, M.; Ejiri, H.; Frekers, D.; Fujita, H.; Fujita, Y.; Fujiwara, M.; Grewe, E.W.; et al. High-resolution 100Mo (3He,t) charge-exchange experiment and the impact on double-β decays and neutrino charged-current reactions. Phys. Rev. 2012, 86, 44309. [Google Scholar]
  41. Puppe, P.; Lennarz, A.; Adachi, T.; Akimune, H.; Ejiri, H.; Frekers, D.; Fujita, H.; Fujita, Y.; Fujiwara, M.; Grewe, E.W.; et al. High resolution (3He,t) experiment on the double-β decaying nuclei 128Te and 130Te. Phys. Rev. 2012, 86, 044603. [Google Scholar]
  42. Ejiri, H.; Frekers, D. Spin dipole nuclear matrix elements for double beta decay nuclei by charge-exchange reactions. J. Phys. G Nucl. Part. Phys. 2016, 43, 11LT01. [Google Scholar] [CrossRef] [Green Version]
  43. Akimune, H.; Ejiri, H.; Hattori, F.; Agodi, C.; Alanssari, M.; Cappuzzello, F.; Carbone, D.; Cavallaro, M.; Diel, F.; Douma, C.A.; et al. Spin-dipole nuclear matrix element for double beta decay of 76Ge by the (3He,t) charge-exchange reaction. J. Phys. G Nucl. Part. Phys. 2020, 47, 05LT01. [Google Scholar] [CrossRef]
  44. Ejiri, H. Axial-vector weak coupling at medium momentum for astro neutrinos and double beta decays. J. Phys. G Nucl. Part Phys. 2019, 46, 125202. [Google Scholar] [CrossRef] [Green Version]
  45. Measday, D.F. The nuclear physics of muon capture. Phys. Rep. 2001, 354, 243–409. [Google Scholar] [CrossRef]
  46. Hashim, I.; Ejiri, H.; Shima, T.; Takahisa, K.; Sato, A.; Kuno, Y.; Ninomiya, K.; Kawamura, N.; Miyake, Y. Muon capture reaction on 100Mo to study neutrino nuclear responses for double-β decays and neutrinos of astro-physics origins. Phys. Rev. 2018, 97, 014617. [Google Scholar]
  47. Jokiniemi, L.; Suhonen, J.; Ejiri, H.; Hashim, I. Pining down the strength function for ordinary muon capture on 100Mo. Phys. Lett. 2019, 794, 143–147. [Google Scholar] [CrossRef]
  48. Primakoff, H. Theory of muon capture. Rev. Mod. Phys. 1959, 31, 802–822. [Google Scholar] [CrossRef]
  49. Jokiniemi, L.; Suhonen, J. Muon-capture strength functions in intermediate nuclei of 0∋ββ decays. Phys. Rev. 2019, 100, 014619. [Google Scholar]
  50. Ejiri, H.; Soukouti, N.; Suhonen, J. Spin-dipole nuclear matrix elements for double beta decays and astro-neutrinos. Phys. Lett. 2014, 729, 27–32. [Google Scholar] [CrossRef] [Green Version]
  51. Ejiri, H.; Suhonen, J. GT neutrino-nuclear responses for double beta decays and astro-neutrinos. J. Phys. G 2015, 42, 055201. [Google Scholar] [CrossRef]
  52. Simkovic, F.; Dvornieky, R.; Vogel, P. Muon capture rates: Evaluation within the quasiparticle random phase approximation. Phys. Rev. 2020, 102, 034301. [Google Scholar] [CrossRef]
  53. Plasil, F. ORLaND—Oak Ridge Laboratory for Neutrino Detectors. Acta Phys. Hung. A Heavy Ion Phys. 1999, 15, 381–390. [Google Scholar] [CrossRef]
  54. Ejiri, H. Neutrino studies in nuclei and intense neutrino sources. Nucl. Instr. Methods Phys. Res. 2003, 503, 276–278. [Google Scholar] [CrossRef]
  55. Ejiri, H.; Richard, P.; Ferguson, S.; Heffner, R.; Perry, R.D. Electric dipole transition from the f7/2 isobaric analogue state to the 2d5/2 ground state in 141Pr. Phys. Rev. Lett. 1968, 21, 373–376. [Google Scholar] [CrossRef]
  56. Ejiri, H.; Titov, A.; Bosewell, M.; Yang, A. Neutrino nuclear response and photonuclear reactions. Phys. Rev. 2013, 88, 054610. [Google Scholar] [CrossRef] [Green Version]
  57. Ejiri, H. Gamma Ray and Electron Spectroscopy in Nuclear Physics; Clarendon Press: London, UK, 1989. [Google Scholar]
  58. Suhonen, J. Impact of the quenching of gA on the sensitivity of 0νββ experiments. Phys. Rev. 2017, 96, 05501. [Google Scholar] [CrossRef] [Green Version]
  59. Ejiri, H. Nuclear matrix elements for β and ββ decays and quenching of the weak coupling gA in QRPA. Front. Phys. 2019. [Google Scholar] [CrossRef]
  60. Jokiniemi, L.; Suhonen, J.; Ejiri, H. Magnetic hexadecapole γ transitions and neutrino-nuclear responses in medium heavy nuclei. Adv. High Energy Phys. 2016, 2016, 8417598. [Google Scholar] [CrossRef] [Green Version]
  61. Barea, J.; Korita, J.; Iachello, F. Nuclear matrix elements for double-β decay. Phys. Rev. 2013, 87, 014315. [Google Scholar]
  62. Jokiniemi, L.; Ejiri, H.; Frekers, D.; Suhonen, J. Neutrinoless ββ nuclear matrix elements using isovector spin-dipole Jπ=2 data. Phys. Rev. 2018, 98, 024608. [Google Scholar] [CrossRef] [Green Version]
  63. Cappuzzello, F.; Cavallaro, M.; Agodi, C.; Bondì, M.; Carbone, D.; Cunsolo, A.; Foti, A. Heavy ion double charge exchange reactions: A tool toward 0νββ nuclear matrix elements. Eur. Phys. J. 2015, 51, 145. [Google Scholar] [CrossRef] [Green Version]
  64. Gysbers, P.; Hagen, G.; Holt, J.D.; Jansen, G.R.; Morris, T.D.; Navrátil, P.; Papenbrock, T.; Quaglioni, S.; Schwenk, A.; Stroberg, S.R.; et al. Discrepancy between experimental and theoretical β decay rates resolved from first principles. Nat. Phys. 2018, 15, 428–431. [Google Scholar] [CrossRef]
  65. Hardy, J.C.; Towner, I.S. The measurement and interpretation of super allowed 0+→0+ nuclear β decay. J. Phys. G Nucl. Part Phys. 2014, 41, 114004. [Google Scholar] [CrossRef] [Green Version]
  66. Ejiri, H. Perspectives on neutrino nuclear-response studies for double beta decays and astro neutrinos. Front. Phys. 2021. to be published. [Google Scholar]
Figure 1. Schematic diagrams for 2 neutrino double beta decay (DBD) (A) and for neutrinoless DBDs with a Majorana neutrino exchange (B), the SUSY particle exchange (C) and the Majoron emission (D). DBD processes with 2 neutrons → 2 protons in the DBD nucleus are shown. d and u are d and u quarks in the 2 neutrons and the 2 protons, respectively [1,3].
Figure 1. Schematic diagrams for 2 neutrino double beta decay (DBD) (A) and for neutrinoless DBDs with a Majorana neutrino exchange (B), the SUSY particle exchange (C) and the Majoron emission (D). DBD processes with 2 neutrons → 2 protons in the DBD nucleus are shown. d and u are d and u quarks in the 2 neutrons and the 2 protons, respectively [1,3].
Universe 06 00225 g001
Figure 2. Effective neutrino masses m e f f in the cases of normal mass-hierarchy (NH) and inverted mass-hierarchy (IH) as a function of the minimum mass m 0 and the mass constrains from the cosmological bound of m i 110 meV.
Figure 2. Effective neutrino masses m e f f in the cases of normal mass-hierarchy (NH) and inverted mass-hierarchy (IH) as a function of the minimum mass m 0 and the mass constrains from the cosmological bound of m i 110 meV.
Universe 06 00225 g002
Figure 3. Effective neutrino masses m e f f in the case of IH as a function of the phase difference α 2 α 1 .
Figure 3. Effective neutrino masses m e f f in the case of IH as a function of the phase difference α 2 α 1 .
Universe 06 00225 g003
Figure 4. Double beta decay (DBD) neutrino-mass sensitivities m m in the cases of the NMEs M 0 ν = 1, 2 and B = 1/t y and 0.01/t y as a function of the exposure N T . Here k = 40 meV, η = 0.9 and ϵ = 0.6 are assumed. The arrows with IH and NH are the mass regions for the IH and NH cases.
Figure 4. Double beta decay (DBD) neutrino-mass sensitivities m m in the cases of the NMEs M 0 ν = 1, 2 and B = 1/t y and 0.01/t y as a function of the exposure N T . Here k = 40 meV, η = 0.9 and ϵ = 0.6 are assumed. The arrows with IH and NH are the mass regions for the IH and NH cases.
Universe 06 00225 g004
Figure 5. Effective neutrino mass m e f f and NME M 0 ν . IH: the IH effective mass. Two lines show upper limits in the case of the halflife limits of 10 26 y for 76 Ge and 136 Xe.
Figure 5. Effective neutrino mass m e f f and NME M 0 ν . IH: the IH effective mass. Two lines show upper limits in the case of the halflife limits of 10 26 y for 76 Ge and 136 Xe.
Universe 06 00225 g005
Figure 6. Two-neutrino DBD NMEs (triangles ) derived from the observed 2-neutrino DBD rates and the (fermi-surface) quasi-particle (FSQP) model NMEs (squares) based on the experimental quenching coefficients for single- β GT NMEs [8,38].
Figure 6. Two-neutrino DBD NMEs (triangles ) derived from the observed 2-neutrino DBD rates and the (fermi-surface) quasi-particle (FSQP) model NMEs (squares) based on the experimental quenching coefficients for single- β GT NMEs [8,38].
Universe 06 00225 g006
Figure 7. Energy spectrum of the ( 3 He,t) charge-exchange reactions (CER) on 76 Ge. The spin-dipole (SD) (blue) and Gamow–Teller (GT) (red/yellow) peaks are well separated, as shown in the red circle [39].
Figure 7. Energy spectrum of the ( 3 He,t) charge-exchange reactions (CER) on 76 Ge. The spin-dipole (SD) (blue) and Gamow–Teller (GT) (red/yellow) peaks are well separated, as shown in the red circle [39].
Universe 06 00225 g007
Figure 8. SD NMEs for DBD nuclei. The SD NMEs M C E (SD) (blue triangle) measured by the ( 3 He,t) CERs and the NMEs M β (SD) (red square) derived from the β decay f t values are plotted against the FSQP NMEs M F S Q P (SD).
Figure 8. SD NMEs for DBD nuclei. The SD NMEs M C E (SD) (blue triangle) measured by the ( 3 He,t) CERs and the NMEs M β (SD) (red square) derived from the β decay f t values are plotted against the FSQP NMEs M F S Q P (SD).
Universe 06 00225 g008
Figure 9. Schematic diagram of muon CER and the transition rate as a function of excitation energy.
Figure 9. Schematic diagram of muon CER and the transition rate as a function of excitation energy.
Universe 06 00225 g009
Figure 10. Energy spectrum for neutrinos from decays of π + and μ + produced by the SNS GeV-proton nuclear interactions.
Figure 10. Energy spectrum for neutrinos from decays of π + and μ + produced by the SNS GeV-proton nuclear interactions.
Universe 06 00225 g010
Figure 11. Schematic diagram of β decay from |i> to |f> and γ transition from |f> to the isobaric analog state (IAS) state |IAS> = T |i> with T being the isospin lowering operator. They are analogous with each other.
Figure 11. Schematic diagram of β decay from |i> to |f> and γ transition from |f> to the isobaric analog state (IAS) state |IAS> = T |i> with T being the isospin lowering operator. They are analogous with each other.
Universe 06 00225 g011
Figure 12. Angular distributions of photo-neutrons via the isobaric analog states from polarized photons on 100 Mo [8,56].
Figure 12. Angular distributions of photo-neutrons via the isobaric analog states from polarized photons on 100 Mo [8,56].
Universe 06 00225 g012
Figure 13. Quenching coefficients g A e f f / g A and the spin J for GT(J = 1), SD (J = 2) and HD (J = 4) NMEs. M EXP and M QR are the experimental and pnQRPA model NMEs, respectively.
Figure 13. Quenching coefficients g A e f f / g A and the spin J for GT(J = 1), SD (J = 2) and HD (J = 4) NMEs. M EXP and M QR are the experimental and pnQRPA model NMEs, respectively.
Universe 06 00225 g013
Figure 14. Neutrinoless DBD NMEs as a function of the quenching coefficient g A e f f / g A for 130 Te. M 0 ν : GT + T + F. M 0 ν (GT) + M 0 ν (T): GT + T. M 0 ν (F): F.
Figure 14. Neutrinoless DBD NMEs as a function of the quenching coefficient g A e f f / g A for 130 Te. M 0 ν : GT + T + F. M 0 ν (GT) + M 0 ν (T): GT + T. M 0 ν (F): F.
Universe 06 00225 g014
Table 1. Nuclear sensitivities for DBD nuclei of the current interest. Shown are the isotope, its natural abundance (N.a.), the Q value, and the nuclear sensitivity. k: the nuclear kinetic factor, m n * : nuclear sensitivity in the case of M 0 ν = 2, and  m n * * : nuclear sensitivity in the case of M 0 ν = 3.
Table 1. Nuclear sensitivities for DBD nuclei of the current interest. Shown are the isotope, its natural abundance (N.a.), the Q value, and the nuclear sensitivity. k: the nuclear kinetic factor, m n * : nuclear sensitivity in the case of M 0 ν = 2, and  m n * * : nuclear sensitivity in the case of M 0 ν = 3.
NuclideN.a.(%) Q Value (keV) k meV m n * meV m n * * meV
76 Ge7.44203980.840.426.7
82 Se8.73299740.020.013.3
100 Mo9.63303434.817.411.6
116 Cd7.49281436.01812
130 Te33.8252840.220.113.4
136 Xe8.9245840.220.113.4
150 Nd5.64337119.89.96.6
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ejiri, H. Neutrino-Mass Sensitivity and Nuclear Matrix Element for Neutrinoless Double Beta Decay. Universe 2020, 6, 225. https://doi.org/10.3390/universe6120225

AMA Style

Ejiri H. Neutrino-Mass Sensitivity and Nuclear Matrix Element for Neutrinoless Double Beta Decay. Universe. 2020; 6(12):225. https://doi.org/10.3390/universe6120225

Chicago/Turabian Style

Ejiri, Hiroyasu. 2020. "Neutrino-Mass Sensitivity and Nuclear Matrix Element for Neutrinoless Double Beta Decay" Universe 6, no. 12: 225. https://doi.org/10.3390/universe6120225

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