Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Symmetries and the Hilbert Space of Large N Extended States
Next Article in Special Issue
Re-Examination of the Effect of Pairing Gaps on Gamow–Teller Strength Distributions and β-Decay Rates
Previous Article in Journal
Revisiting a Core–Jet Laboratory at High Redshift: Analysis of the Radio Jet in the Quasar PKS 2215+020 at z = 3.572
Previous Article in Special Issue
Integral Fluxes of Neutrinos and Gamma-Rays Emitted from Neighboring X-ray Binaries
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Systematic Study of Two-Neutrino Double Electron Capture

1
Faculty of Mathematics, Physics and Informatics, Comenius University in Bratislava, 842 48 Bratislava, Slovakia
2
International Centre for Advanced Training and Research in Physics, P.O. Box MG-12, RO-077125 Măgurele, Romania
3
“Horia Hulubei” National Institute of Physics and Nuclear Engineering, 30 Reactorului, P.O. Box MG-6, RO-077125 Măgurele, Romania
4
Department of Physics, University of Bucharest, 405 Atomiştilor, P.O. Box MG-11, RO-077125 Măgurele, Romania
5
Institute of Experimental and Applied Physics, Czech Technical University in Prague, 110 00 Prague, Czech Republic
*
Author to whom correspondence should be addressed.
Universe 2024, 10(2), 98; https://doi.org/10.3390/universe10020098
Submission received: 11 January 2024 / Revised: 14 February 2024 / Accepted: 16 February 2024 / Published: 17 February 2024

Abstract

:
In this paper, we update the phase-space factors for all two-neutrino double electron capture processes. The Dirac–Hartree–Fock–Slater self-consistent method is employed to describe the bound states of captured electrons, enabling a more realistic treatment of atomic screening and more precise binding energies of the captured electrons compared to previous investigations. Additionally, we consider all s-wave electrons available for capture, expanding beyond the K and L 1 orbitals considered in prior studies. For light atoms, the increase associated with additional captures compensates for the decrease in decay rate caused by the more precise atomic screening. However, for medium and heavy atoms, an increase in the decay rate, up to 10% for the heaviest atoms, is observed due to the combination of these two effects. In the systematic analysis, we also include capture fractions for the first few dominant partial captures. Our precise model enables a close examination of low Q-value double electron capture in 152Gd, 164Er, and 242Cm, where partial K K captures are energetically forbidden. Finally, with the updated phase-space values, we recalculate the effective nuclear matrix elements and compare their spread with those associated with 2 ν β β decay.

1. Introduction

Exploring second-order weak interaction nuclear transitions represents a crucial channel in addressing unresolved challenges within modern physics. The detection of neutrinoless double beta decay (DBD) holds particular significance, as it would confirm the Majorana nature of neutrinos [1] (meaning they are their own antiparticles), and potentially provide insights into the mechanism generating their mass, as well as into charge-parity (CP) violation in the lepton sector [2,3,4,5,6,7,8]. Moreover, it would mark the first instance of the detection of a lepton number-violating process in laboratory. These types of beyond the Standard Model (SM) processes might play a significant role in certain leptogenesis mechanisms, which may elucidate the prevalence of matter over antimatter in the Universe [9,10].
Currently, the most studied modes of DBD involve the emission of two electrons from neutron-rich candidate isotopes: the two-neutrino mode ( 2 ν β β ) and the neutrinoless mode ( 0 ν β β ). While the latter remains experimentally elusive [11,12,13,14,15,16,17,18], the former serves as a playground to test different hypotheses and the predictive strengths of the nuclear structure models [19,20,21], and to constrain various parameters associated with physics beyond the SM [22,23,24,25,26,27,28,29,30,31,32]. A comprehensive and recent review on the latter topic can be found in [33].
DBD of proton-rich isotopes is also possible and is currently in the early stages of experimental exploration [34]. There are three such DBD modes [35,36]: the double-positron emitting ( 0 ν / 2 ν β + β + ) mode, the single electron capture with coincident positron emission ( 0 ν / 2 ν EC β + ) mode, and the double electron capture ( 0 ν / 2 ν ECEC ) mode. Experimental measurements of these modes have received less attention due to their longer half-lives compared to the 2 ν β β decay mode, mostly due to their lower Q-values. However, their distinct signatures might help in their detection, especially if coincidence trigger logic can be employed [34]. There are some positive indications of the 2 ν ECEC mode for 130Ba and 132Ba from geochemical measurements [37,38,39] and for 78Kr [40,41]. Recently, the XENON1T collaboration reported the direct observation of the 2 ν ECEC mode in 124Xe [42,43].
From the theoretical point of view, the pioneer estimates of the decay rates for 2 ν ECEC have been obtained by Primakoff and Rosen [44,45], followed up by the ones made by Vergados [46] and by Kim and Kubodera [47]. In these early stages, the predictions were based on a non-relativistic treatment for the captured electrons. Doi and Kotani were the first to present a detailed theoretical formulation for the process and to include the relativistic effects [35]. Still, the calculations are based on the atomic structure of a point-like nucleus, for which the electron bound states are known analytically. At the same time, Boehm and Vogel presented the results for some selected 2 ν ECEC cases, but not many details of the calculations are provided [48]. Recently, the treatment of the atomic screening was improved using a Thomas–Fermi model for an atomic cloud [49,50,51]. The diffuse nuclear surface effects have also been accounted for through a realistic charge distribution inside the nucleus [50,51].
Although 2 ν ECEC, by its very nature, stands at the interface between nuclear and atomic physics, some aspects regarding the atomic structure calculations have been overlooked or treated simplistically in previous investigations. In this paper, we employ the Dirac–Hartree–Fock–Slater (DHFS) self-consistent framework, providing an enhanced description of atomic screening for the systems involved in the decay and a more rigorous estimation of the binding energies for the captured electrons. We mention that the DHFS atomic potential includes the finite nuclear size and diffuse nuclear surface effects. Furthermore, we extend our analysis to include all s-wave electrons available for capture, surpassing the K and L 1 orbitals considered in prior studies.
We demonstrate that, for light atoms undergoing the 2 ν ECEC process, the decrease in the decay rate resulting from more realistic screening is balanced out by the increase associated with opening the capture from higher orbitals than K and L 1 . However, this balance is not observed for medium and heavy atoms, where the latter improvement in our model leads to a considerable increase in the decay rate, up to 10% for the heaviest atoms. For all cases undergoing 2 ν ECEC, we provide the capture fractions for the first few dominant channels. We specifically address the low Q-value 2 ν ECEC transitions of 152Gd, 164Er, and 242Cm, for which the capture of both K shell electrons is energetically forbidden. Finally, with the updated phase-space values, we reexamine the effective nuclear matrix elements and compare their spread with those associated with 2 ν β β .

2. Formalism for the Two-Neutrino Double Electron Capture

2.1. Energetics

We investigate the two neutrino double electron capture ( 2 ν ECEC ) process,
2 e + ( A , Z + 2 ) ( A , Z ) + 2 ν e ,
in which the initial nucleus ( A , Z + 2 ) captures two atomic electrons, changing its atomic number by two units and emitting two neutrinos. If the electrons are captured from x and y atomic orbitals, we denote the process as 2 ν xy . To keep the notation simple, we use the X-ray notation for the orbitals, i.e., x and y = K , L 1 , M 1 , N 1 , Following the 2 ν xy process, the final atomic system, while electrically neutral, remains in an excited atomic state, with vacancies in the orbitals x and y from which the electrons were initially captured. Consequently, the 2 ν ECEC process is followed by atomic de-excitation, involving a sequence of X-ray emissions and Auger electrons released from the outer shells. Figure 1 provides a schematic illustration of two channels of the 2 ν ECEC process in   124 Xe. The top panel shows that both electrons are captured from the K orbital ( 2 ν K K ), while the bottom panel shows that they are captured from L 1 orbital ( 2 ν L 1 L 1 ). It is worth noting that the atomic de-excitation results in different energy deposition in the detector depending on whether the capture is from the K or L 1 shells. This is due to the fact that the atomic relaxation following L 1 L 1 capture emits less energy compared to the K K capture.
We adopt the following energy balance,
M gs ( A , Z + 2 ) M gs ( A , Z ) = ω x + ω y + M xy ( A , Z ) M gs ( A , Z ) ,
which indicates that the difference in the atomic masses of the initial and final systems in ground states, usually denoted as Q = M gs ( A , Z + 2 ) M gs ( A , Z ) , is shared between neutrinos with energies ω x and ω y , and the atomic relaxation of the final atom. The mass M xy ( A , Z ) corresponds to an atomic excited state of the final neutral system, with holes in the shell x and y from where the electrons were captured in the initial atom, and it can be written as
M xy ( A , Z ) = E f + Z m e B xy ( Z )
where E f is the nuclear mass (energy) of the final nucleus, Z m e is the rest energy of the final atomic cloud ( m e is the mass of the electron), and B xy ( Z ) is the total electron binding energy of the final system with a configuration with two holes. We adopt the units in which = c = 1 . One can see that the atomic relaxation energy can be written as a difference in the total electron binding energies of the final atom with gs and xy configurations,
R xy = B gs ( Z ) B xy ( Z ) .
The total energy of the emitted neutrinos can be written as
ω x + ω y = Q R xy
and, making the usual approximation that B xy ( Z ) B gs ( Z + 2 ) = t x t y , we finally obtain [35,49,51]
ω x + ω y = E i E f + 2 m e t x t y = E i E f + e x + e y ,
where t x and e x = m e t x are, respectively, the binding and the total energy for the electrons in shell x .

2.2. Half-Life

Following the separation of the decay rate from [35,49], the inverse half-life of the 2 ν ECEC process from 0 + ground state of the initial nucleus to the 0 + ground state of the final nucleus is given by
T 1 / 2 2 ν ECEC 1 = g A 4 M 2 ν ECEC 2 G K edge 2 ν ECEC ,
where g A is the axial coupling constant, G K edge 2 ν ECEC is the phase-space factor (PSF) discussed later, and M 2 ν ECEC , the nuclear matrix element (NME), is given by
M 2 ν ECEC = m e A ˜ M GT 2 ν ECEC g V g A 2 M F 2 ν ECEC .
Here, A ˜ = 1.12 A 1 / 2 in MeV and g V is the vector coupling constant. This choice of A ˜ reproduces the average excitation energies in the intermediate nucleus for a wide range of isotopes [53]. Nevertheless, there are reasons supporting the use of a single value for both, as discussed in [53]. The double Gamow–Teller (GT) and Fermi (F) transition matrix elements, in closure approximation, are M GT 2 ν ECEC = 0 f + j , k τ j τ k σ j σ k 0 i + and M F 2 ν ECEC = 0 f + j , k τ j τ k 0 i + [49]. Here, τ j , k is the isospin-lowering operator, transforming a proton into a neutron and σ j , k is the nucleon spin operator.
It is important to highlight that the factorization of the half-life in Equation (7) assumes the closure approximation [35]. Consequently, the PSFs outlined below should be paired with NMEs computed within the same approximation. However, there have been various studies computing NMEs beyond the closure approximation and overlooking the constraint above [54,55,56]. Despite the lack of consistency, the predictions for half-lives are notably close to experimental values. In this study, we adhere to the closure approximation for comparison purposes. It is worth mentioning that an alternative factorization of the 2 ν ECEC half-life, beyond the closure approximation, will be presented elsewhere.
In this paper, we consider that captures can occur from any occupied s 1 / 2 orbitals of the initial atom, i.e., K , L 1 , M 1 , The capture from other orbitals is suppressed, firstly because the higher orbital angular momentum of these orbitals reduces the capture probability, and secondly, smaller NMEs are associated with these captures [35]. Under these assumptions, we found that the PSF expression is
G K edge 2 ν ECEC = m e ( G F V u d m e 2 ) 4 16 π 3 2 A ˜ 3 m e 2 1 m e 3 x , y = K E I E F + e x + e y > 0 edge B x 2 B y 2 × 0 E I E F + e x + e y K n , xy 2 + L n , xy 2 + K n , xy L n , xy ω x 2 ω y 2 d ω x = x , y = K E I E F + e x + e y > 0 edge G 2 ν xy .
where G F is the Fermi coupling constant and V u d the first element of the Cabibbo–Kobayashi–Maskawa (CKM) matrix. The sums are running over all occupied orbitals of the initial atomic system from where the electrons can be captured. The “edge” label refers to the last such orbital. The localization probability for an electron from shell x on the nuclear surface can be written as [35,49],
B x 2 = 1 4 π m e 3 g x 2 ( R ) + f x 2 ( R ) ,
in terms of the large- and small-component radial wave functions describing the bound state, g x ( r ) and f x ( r ) , respectively, evaluated on the nuclear surface R = 1.2 A 1 / 3 . The factors K n , xy and L n , xy can be written as
K n , xy = 1 e x + ω x + E n E i + 1 e y + ω y + E n E i L n , xy = 1 e x + ω y + E n E i + 1 e y + ω x + E n E i
where E n is a suitably chosen excitation energy for the intermediate nucleus, ( A , Z + 1 ) . In the actual calculations, the energy difference E n E i = A ˜ Q / 2 + m e .
For comparison with the results of the previous investigations, we define also the PSF for captures from K and L 1 orbitals only,
G K L 1 2 ν ECEC = G 2 ν K K + 2 G 2 ν K L 1 + G 2 ν L 1 L 1
which is a particularization of Equation (9). The factor of 2 stems from the equal contributions of 2 ν xy and 2 ν yx captures in the total decay rate for any x y orbitals.

2.3. Electron Bound States Description

Accurate computation of the PSFs for any 2 ν ECEC process requires precise atomic structure calculations for the atomic systems involved in the process. For this purpose, we employed the Dirac–Hartree–Fock–Slater (DHFS) self-consistent framework. The nuclear, electronic, and exchange components of the DHFS potential, as well as the convergence of the self-consistent method, were described in detail in [57,58], based on the RADIAL subroutine package [59], which we also used in our calculations. In a relativistic framework and for a spherically symmetric potential of the atomic systems, the large- and small-component radial wave functions from Equation (10) obey the following system of coupled differential equations
d d r + κ + 1 r g n κ ( r ) ( e n κ V DHFS ( r ) + m e ) f n κ ( r ) = 0 , d d r κ 1 r f n κ ( r ) + ( e n κ V DHFS ( r ) m e ) g n κ ( r ) = 0 .
In previous sections, the shells referred to as x or y can be uniquely identified by the state n κ , where n represents the principal quantum number and κ corresponds to the relativistic quantum number.
In Table 1, we present a comparison between experimental binding energies and values predicted by different models for the selected isotopes (Kr, Xe, and Ba) undergoing the 2 ν ECEC process. The model proposed by [35], assuming a uniform charge distribution for the nucleus with no screening due to the atomic charge, fails to reproduce the experimental data. Models incorporating atomic screening and diffuse nuclear surface effects, as seen in [50,51], using a Thomas–Fermi screening function and a realistic proton density inside the nucleus, respectively, provide improved results but still yield binding energies significantly larger than the experimental ones. Notably, the DHFS self-consistent method employed in this work successfully reproduces experimental binding energies within a 1% margin.

3. Results and Discussions

The PSFs for captures from K and L 1 shells only, G K L 1 2 ν ECEC , calculated for the 2 ν ECEC cases previously investigated, are shown in Table 2. Here, we have assumed the same Q-values as in [49]. Our results are consistently lower by about 5% than the values reported in [49], with a couple of exceptions observed in   92 Mo and   180 W. Our results generally fall within 70% of the values reported in [50,51], with a noticeable dependency on A: for low mass numbers, our PSFs are lower than those in [50,51], but as the mass number increases, the factors in this work approach and even exceed the values reported in that reference, suggesting a discrepancy in the underlying models. In contrast, the G K L 1 2 ν ECEC values derived in this work exhibit a consistent 10–20% reduction compared to the ones provided in [35], which could be attributed to the fact that, in [35], the screening effect is not accounted for. We chose not to display the G K L 1 2 ν ECEC values from [48], since it does not employ the same definition of the PSF due to a different separation of the decay rate.
In Table 3, we present the total PSFs G K edge 2 ν ECEC , as well as the capture fractions for selected shell pairs, for all nuclei for which 2 ν ECEC is energetically possible, and EC is energetically forbidden. When comparing with one of the most precise calculations available [49] (column four of Table 2), one can see the following trend. For light atoms, the more rigorous treatment of atomic screening balances out the effect of accounting for all κ = 1 shells. The only exception is   40 Ca, where the PSF value computed using all κ = 1 shells is higher by about 7.5% than the one in [49]. The difference is induced by the interplay between the binding energies and Q-value in the PSF integral (see Equation (9)). For medium and heavy atoms, there is an increase in the decay rate compared to results from [49], almost linear in Z, reaching about 10% for the heaviest atoms.
The low Q-value transitions of   152 Gd,   164 Er and   242 Cm exhibit an interesting behavior. In these cases, both the K K and the K L 1 captures are energetically forbidden, and the highest contribution to the total PSF is given by the L 1 L 1 , M 1 M 1 and L 1 M 1 captures. The values of the PSFs and capture fractions for these low Q-value transitions are given in Table 4.
The inverse of the PSF is proportional to the half-life for each transition (see Equation (7)). This quantity is plotted in Figure 2 for all nuclei shown in Table 3 versus atomic number. One can see a decrease in the Q-value as the atomic number increases.
Finally, we investigate the effective matrix elements for the 2 ν ECEC process of   78 Kr,   124 Xe,   130 Ba and   132 Ba. These can be obtained as
M eff 2 ν ECEC = 1 T 1 / 2 2 ν ECEC G K edge 2 ν ECEC = m e A ˜ g A 2 M GT 2 ν ECEC g V 2 M F 2 ν ECEC .
using the experimental half-lives and the PSFs from Table 3. The effective NME values are presented in Table 5. It should be noted that, in the case of   132 Ba, there is also an experimental half-life measurement [38]. The authors state that the obtained value is tentative, and indeed, we obtain an effective NME value that is more than one order of magnitude larger than the others. Hence, we employ the half-life limit from [37]. Future measurements might clarify the situation, but we note that, in general, geochemical measurements tend to underestimate the half-lives, leading to overestimated effective NMEs. This may be a consequence of the difficulty in identifying the relevant production channel of the final atom [34]. Another interesting hypothesis is the variation in the weak interaction constant with time [62,63,64].
In Figure 3, we compare the values from Table 5 with the corresponding 2 ν β β effective nuclear matrix elements. These are defined similarly to the ECEC effective NMEs, but replacing isospin-lowering with isospin-rising operators in the definition of the Gamow–Teller and Fermi matrix elements. We note that effective matrix elements from both processes span similar ranges.

4. Uncertainties and Further Improvements

We first perform a sensitivity analysis on the input parameters, namely the Q-value, the average energy of the excited 1 + states E N and the nuclear radius. The relative variations in the PSFs due to each parameter are summarized in Table 6, with the Q-value having the highest influence. However, we note that modern determinations of this quantity are highly accurate and the absolute values of the PSFs are practically insensitive to the choice of Q-value. The sensitivity of the PSFs with respect to E N and the nuclear radius are much smaller, as noted previously in [49].
Another source of uncertainties is related to the choice of the DHFS model for atomic structure computations. The accuracy of the binding energies obtained in the DHFS approach is around 1% relative to the experimental values [58]. The binding energies enter in the PSFs definition through the integration limit, through the energies of the emitted neutrinos, and through the K N and L N factors, always as additive quantities to the Q-value. Consequently, this uncertainty has a negligible effect as it translates to O ( 10 4 ) relative variation in the Q-value of the process. Low Q-value transitions are an exception, and the uncertainty in the binding energies is dominant. The values of the wave functions on the nuclear surface can also introduce some uncertainty in the PSFs. This effect was studied in [58] by comparing the Coulomb amplitudes (proportional to the wave function values on the nuclear surface) in the DHFS model with the ones obtained through more complex models. The study found that Coulomb amplitudes obtained in the DHFS framework agree within 0.25% with ones obtained through the more refined Dirac–Hartree–Fock (DHF) model for the 1 s 1 / 2 and 2 s 1 / 2 shells and for atomic numbers above 20. Considering this value as uncertainty within our model, we evaluate the relative error in the PSFs due to the wave functions alone to be 1%, since the wave functions enter the PSF definition to the fourth power. We note that this uncertainty is a systematic one, the Coulomb amplitudes in the DHFS model always overestimating the ones obtained in the DHF model.
Besides the uncertainties discussed above, one may expect variations in PSFs from improving the overall model of the 2 ν ECEC process. Firstly, the summed neutrino energy determination has been obtained through a usual approximation (see Equation (6)), but it can be more precisely determined through Equation (5). Consequences are most important in low Q-value transitions. For example, the DHFS framework predicts R K L 1 = 55.05 keV for   152 Sm. Consequently, the 2 ν K L 1 process of   152 Gd is energetically allowed (in contrast to the results shown in Table 3), as it was also found in [67,68,69] when studying the resonant neutrinoless double electron capture. Secondly, the Pauli blocking of the decay of the innermost nucleon states is not accounted for in our model. This aspect can be improved by averaging the bound electron wave function, weighted with a realistic nuclear charge distribution [65,70]. From preliminary results, we expect an increase of a few % in the PSFs due to the Pauli blocking effect.

5. Conclusions

In this paper, we have conducted a systematic study of atoms undergoing the 2 ν ECEC process. Despite the fact that 2 ν ECEC, by its very nature, resides at the interface between nuclear and atomic physics, certain aspects related to atomic structure calculations have been either overlooked or treated simplistically in previous investigations. In our model, we introduced two enhancements to the atomic part calculations of double electron capture transitions. Firstly, by employing the DHFS self-consistent framework to describe the initial atomic system, we refined both the distortion of bound wave functions due to atomic clouds and the estimation of binding energies for captured electrons. Secondly, we extended our consideration beyond previous studies, where only K and L 1 orbitals were included, by also allowing captures from the outer orbitals of the initial atoms.
With our improved model, we have updated the phase-space values for all atoms undergoing the 2 ν ECEC process. For light atoms, we observed almost no differences compared to a previous model, where a simplified atomic screening and only captures from K and L 1 orbitals were considered. This lack of difference is associated with a cancellation of the decrease in the decay rate from a more precise screening treatment with the increase associated with the possibility of captures occurring from higher orbitals. In contrast, for medium and heavy atoms, we observed an increase in the decay rate, almost linear in Z, reaching about 10% for the heaviest atoms. In the systematic study, we also provided capture fractions for the first few dominant partial channels and separately addressed the low Q-value 2 ν ECEC transitions of   152 Gd,   164 Er, and   242 Cm, for which the K K capture is energetically forbidden. Finally, we demonstrated that the effective nuclear matrix elements for both 2 ν ECEC processes and 2 ν β β decays span similar ranges.

Author Contributions

Conceptualization, O.N., S.G., S.S. and F.Š.; methodology, O.N. and S.G.; software, S.G.; validation, O.N., S.G., S.S. and F.Š.; formal analysis, S.S. and F.Š.; investigation, O.N. and S.G.; resources, S.S. and F.Š.; data curation, O.N. and S.G.; writing—original draft preparation, O.N. and S.G.; writing—review and editing, O.N., S.G., S.S. and F.Š.; visualization, O.N. and S.G.; supervision, S.S. and F.Š.; project administration, F.Š.; funding acquisition, S.S. and F.Š. All authors have read and agreed to the published version of the manuscript.

Funding

O.N., S.G. and S.S. acknowledge support from the Romanian Ministry of Research, Innovation, and Digitalization through Project No. PN 23 08 641 04 04/2023. F.Š. acknowledges support from the Slovak Research and Development Agency under Contract No. APVV-22-0413 and by the Czech Science Foundation (GAČR), project No. 24-10180S.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Majorana, E. Teoria simmetrica dell’elettrone e del positrone. II Nuovo C. 1937, 14, 171–184. [Google Scholar] [CrossRef]
  2. Schechter, J.; Valle, J.W.F. Neutrino decay and spontaneous violation of lepton number. Phys. Rev. D 1982, 25, 774–783. [Google Scholar] [CrossRef]
  3. Sujkowski, Z.; Wycech, S. Neutrinoless double electron capture: A tool to search for Majorana neutrinos. Phys. Rev. C 2004, 70, 052501. [Google Scholar] [CrossRef]
  4. Pascoli, S.; Petcov, S.; Schwetz, T. The absolute neutrino mass scale, neutrino mass spectrum, Majorana CP-violation and neutrinoless double-beta decay. Nucl. Phys. B 2006, 734, 24–49. [Google Scholar] [CrossRef]
  5. Bilenky, S.M.; Giunti, C. Neutrinoless double-beta decay: A probe of physics beyond the Standard Model. Int. J. Mod. Phys. A 2015, 30, 1530001. [Google Scholar] [CrossRef]
  6. Vergados, J.D.; Ejiri, H.; Šimkovic, F. Neutrinoless double beta decay and neutrino mass. Int. J. Mod. Phys. E 2016, 25, 1630007. [Google Scholar] [CrossRef]
  7. Girardi, I.; Petcov, S.; Titov, A. Predictions for the Majorana CP violation phases in the neutrino mixing matrix and neutrinoless double beta decay. Nucl. Phys. B 2016, 911, 754–804. [Google Scholar] [CrossRef]
  8. Šimkovic, F. Neutrino masses and interactions and neutrino experiments in the laboratory. Physics-Uspekhi 2021, 64, 1238. [Google Scholar] [CrossRef]
  9. Fukugita, M.; Yanagida, T. Barygenesis without grand unification. Phys. Lett. B 1986, 174, 45–47. [Google Scholar] [CrossRef]
  10. Buchmüller, W.; Peccei, R.; Yanagida, T. Leptogenesis as the origin of matter. Annu. Rev. Nucl. Part. Sci. 2005, 55, 311–355. [Google Scholar] [CrossRef]
  11. Ebert, J.; Fritts, M.; Gehre, D.; Gößling, C.; Hagner, C.; Heidrich, N.; Klingenberg, R.; Kröninger, K.; Nitsch, C.; Oldorf, C.; et al. Results of a search for neutrinoless double-β decay using the COBRA demonstrator. Phys. Rev. C 2016, 94, 024603. [Google Scholar] [CrossRef]
  12. Barabash, A.S.; Belli, P.; Bernabei, R.; Cappella, F.; Caracciolo, V.; Cerulli, R.; Chernyak, D.M.; Danevich, F.A.; d’Angelo, S.; Incicchitti, A.; et al. Final results of the Aurora experiment to study 2β decay of 116Cd with enriched 116CdWO4 crystal scintillators. Phys. Rev. D 2018, 98, 092007. [Google Scholar] [CrossRef]
  13. 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 Dataset. Phys. Rev. Lett. 2019, 123, 161802. [Google Scholar] [CrossRef] [PubMed]
  14. Agostini, M. et al. [GERDA Collaboration] Final Results of GERDA on the Search for Neutrinoless Double-β Decay. Phys. Rev. Lett. 2020, 125, 252502. [Google Scholar] [CrossRef] [PubMed]
  15. Azzolini, O.; Beeman, J.W.; Bellini, F.; Beretta, M.; Biassoni, M.; Brofferio, C.; Bucci, C.; Capelli, S.; Caracciolo, V.; Cardani, L.; et al. Final Result on the Neutrinoless Double Beta Decay of 82Se with CUPID-0. Phys. Rev. Lett. 2022, 129, 111801. [Google Scholar] [CrossRef]
  16. Adams, D.Q. et al. [CUORE Collaboration] New Direct Limit on Neutrinoless Double Beta Decay Half-Life of 128Te with CUORE. Phys. Rev. Lett. 2022, 129, 222501. [Google Scholar] [CrossRef] [PubMed]
  17. Abe, S. et al. [KamLAND-Zen Collaboration] Search for the Majorana Nature of Neutrinos in the Inverted Mass Ordering Region with KamLAND-Zen. Phys. Rev. Lett. 2023, 130, 051801. [Google Scholar] [CrossRef]
  18. Arnquist, I.J. et al. [Majorana Collaboration] Final Result of the Majorana Demonstrator’s Search for Neutrinoless Double-β Decay in 76Ge. Phys. Rev. Lett. 2023, 130, 062501. [Google Scholar] [CrossRef]
  19. Saakyan, R. Two-Neutrino Double-Beta Decay. Annu. Rev. Nucl. Part. Sci. 2013, 63, 503–529. [Google Scholar] [CrossRef]
  20. Barabash, A. Precise Half-Life Values for Two-Neutrino Double-β Decay: 2020 Review. Universe 2020, 6, 159. [Google Scholar] [CrossRef]
  21. Barabash, A. Double Beta Decay Experiments: Recent Achievements and Future Prospects. Universe 2023, 9, 290. [Google Scholar] [CrossRef]
  22. Deppisch, F.F.; Graf, L.; Rodejohann, W.; Xu, X.J. Neutrino self-interactions and double beta decay. Phys. Rev. D 2020, 102, 051701. [Google Scholar] [CrossRef]
  23. Deppisch, F.F.; Graf, L.; Šimkovic, F. Searching for New Physics in Two-Neutrino Double Beta Decay. Phys. Rev. Lett. 2020, 125, 171801. [Google Scholar] [CrossRef]
  24. Barabash, A.; Dolgov, A.; Dvornický, R.; Šimkovic, F.; Smirnov, A. Statistics of neutrinos and the double beta decay. Nucl. Phys. B 2007, 783, 90–111. [Google Scholar] [CrossRef]
  25. Bolton, P.D.; Deppisch, F.F.; Gráf, L.; Šimkovic, F. Two-neutrino double beta decay with sterile neutrinos. Phys. Rev. D 2021, 103, 055019. [Google Scholar] [CrossRef]
  26. Agostini, M.; Bossio, E.; Ibarra, A.; Marcano, X. Search for light exotic fermions in double-beta decays. Phys. Lett. B 2021, 815, 136127. [Google Scholar] [CrossRef]
  27. Díaz, J.S. Limits on Lorentz and CPT violation from double beta decay. Phys. Rev. D 2014, 89, 036002. [Google Scholar] [CrossRef]
  28. Niţescu, O.; Ghinescu, S.; Stoica, S. Lorentz violation effects in 2νββ decay. J. Phys. G Nucl. Part. Phys. 2020, 47, 055112. [Google Scholar] [CrossRef]
  29. Ghinescu, S.A.; Niţescu, O.; Stoica, S. Investigation of the Lorentz invariance violation in two-neutrino double-beta decay. Phys. Rev. D 2022, 105, 055032. [Google Scholar] [CrossRef]
  30. Arnold, R. et al. [NEMO Collaboration] Limits on different majoron decay modes of 100Mo and 82Se for neutrinoless double beta decays in the NEMO-3 experiment. Nucl. Phys. A 2006, 765, 483–494. [Google Scholar] [CrossRef]
  31. Arnold, R.; Augier, C.; Barabash, A.S.; Basharina-Freshville, A.; Blondel, S.; Blot, S.; Bongrand, M.; Boursette, D.; Brudanin, V.; Busto, J.; et al. Detailed studies of 100Mo two-neutrino double beta decay in NEMO-3. Eur. Phys. J. C 2019, 79, 440. [Google Scholar] [CrossRef]
  32. Arnold, R. et al. [NEMO-3 Collaboration] Search for Neutrinoless Quadruple-β Decay of 150Nd with the NEMO-3 Detector. Phys. Rev. Lett. 2017, 119, 041801. [Google Scholar] [CrossRef] [PubMed]
  33. Bossio, E.; Agostini, M. Probing beyond the standard model physics with double-beta decays. J. Phys. G Nucl. Part. Phys. 2023, 51, 023001. [Google Scholar] [CrossRef]
  34. Belli, P.; Bernabei, R.; Caracciolo, V. Status and Perspectives of 2ϵ, ϵβ+ and 2β Decays. Particles 2021, 4, 241–274. [Google Scholar] [CrossRef]
  35. Doi, M.; Kotani, T. Neutrino Emitting Modes of Double Beta Decay. Prog. Theor. Phys. 1992, 87, 1207–1231. [Google Scholar] [CrossRef]
  36. Doi, M.; Kotani, T. Neutrinoless Modes of Double Beta Decay. Prog. Theor. Phys. 1993, 89, 139–159. [Google Scholar] [CrossRef]
  37. Barabash, A.S.; Saakyan, R.R. Experimental limits on 2β+, Kβ+, and 2K processes for 130Ba and on 2K capture for 132Ba. Phys. At. Nucl. 1996, 59, 179–184. [Google Scholar]
  38. Meshik, A.P.; Hohenberg, C.M.; Pravdivtseva, O.V.; Kapusta, Y.S. Weak decay of 130Ba and 132Ba: Geochemical measurements. Phys. Rev. C 2001, 64, 035205. [Google Scholar] [CrossRef]
  39. Pujol, M.; Marty, B.; Burnard, P.; Philippot, P. Xenon in Archean barite: Weak decay of 130Ba, mass-dependent isotopic fractionation and implication for barite formation. Geochim. Cosmochim. Acta 2009, 73, 6834–6846. [Google Scholar] [CrossRef]
  40. Gavrilyuk, Y.M.; Gangapshev, A.M.; Kazalov, V.V.; Kuzminov, V.V.; Panasenko, S.I.; Ratkevich, S.S. Indications of 2ν2K capture in 78Kr. Phys. Rev. C 2013, 87, 035501. [Google Scholar] [CrossRef]
  41. Ratkevich, S.S.; Gangapshev, A.M.; Gavrilyuk, Y.M.; Karpeshin, F.F.; Kazalov, V.V.; Kuzminov, V.V.; Panasenko, S.I.; Trzhaskovskaya, M.B.; Yakimenko, S.P. Comparative study of the double-K-shell-vacancy production in single- and double-electron-capture decay. Phys. Rev. C 2017, 96, 065502. [Google Scholar] [CrossRef]
  42. Aprile, E. et al. [XENON Collaboration] Observation of two-neutrino double electron capture in 124Xe with XENON1T. Nature 2019, 568, 532–535. [Google Scholar] [CrossRef] [PubMed]
  43. Aprile, E. et al. [XENON Collaboration] Double-weak decays of 124Xe and 136Xe in the XENON1T and XENONnT experiments. Phys. Rev. C 2022, 106, 024328. [Google Scholar] [CrossRef]
  44. Primakoff, H.; Rosen, S.P. Double beta decay. Rep. Prog. Phys. 1959, 22, 121. [Google Scholar] [CrossRef]
  45. Primakoff, H.; Rosen, S.P. A Corrigendum to the Article ‘Double Beta Decay’ by H. Primakoff and S. P. Rosen in Reports on Progress in Physics. Proc. Phys. Soc. 1961, 78, 464. [Google Scholar] [CrossRef]
  46. Vergados, J. Lepton-violating ββ, β+β+ decays, (e, e+) conversion and double electron capture in gauge theories. Nucl. Phys. B 1983, 218, 109–144. [Google Scholar] [CrossRef]
  47. Kim, C.W.; Kubodera, K. Double-positron decay, positron-emitting K capture, and double K capture. Phys. Rev. D 1983, 27, 2765–2768. [Google Scholar] [CrossRef]
  48. Boehm, F.; Vogel, P. Physics of Massive Neutrinos; Cambridge University Press: Cambridge, UK, 1992. [Google Scholar]
  49. Kotila, J.; Iachello, F. Phase space factors for β+β+ decay and competing modes of double-β decay. Phys. Rev. C 2013, 87, 024313. [Google Scholar] [CrossRef]
  50. Mirea, M.; Pahomi, T.; Stoica, S. Values of the phase space factors involved in double beta decay. Rom. Rep. Phys. 2015, 67, 872–889. [Google Scholar]
  51. Stoica, S.; Mirea, M. Phase Space Factors for Double-Beta Decays. Front. Phys. 2019, 7, 12. [Google Scholar] [CrossRef]
  52. Martínez-Lema, G. et al. [NEXT Collaboration] Sensitivity of the NEXT experiment to Xe-124 double electron capture. J. High Energy Phys. 2021, 2021, 1029–8479. [Google Scholar] [CrossRef]
  53. Haxton, W.; Stephenson, G. Double beta decay. Prog. Part. Nucl. Phys. 1984, 12, 409–479. [Google Scholar] [CrossRef]
  54. Suhonen, J. Double beta decays of 124Xe investigated in the QRPA framework. J. Phys. G Nucl. Part. Phys. 2013, 40, 075102. [Google Scholar] [CrossRef]
  55. Pirinen, P.; Suhonen, J. Systematic approach to β and 2νββ decays of mass A = 100 − 136 nuclei. Phys. Rev. C 2015, 91, 054309. [Google Scholar] [CrossRef]
  56. Coello Pérez, E.; Menéndez, J.; Schwenk, A. Two-neutrino double electron capture on 124Xe based on an effective theory and the nuclear shell model. Phys. Lett. B 2019, 797, 134885. [Google Scholar] [CrossRef]
  57. Niţescu, O.; Stoica, S.; Šimkovic, F. Exchange correction for allowed β decay. Phys. Rev. C 2023, 107, 025501. [Google Scholar] [CrossRef]
  58. Sevestrean, V.A.; Niţescu, O.; Ghinescu, S.; Stoica, S. Self-consistent calculations for atomic electron capture. Phys. Rev. A 2023, 108, 012810. [Google Scholar] [CrossRef]
  59. Salvat, F.; Fernández-Varea, J.M. RADIAL: A Fortran subroutine package for the solution of the radial Schrödinger and Dirac wave equations. Comput. Phys. Commun. 2019, 240, 165–177. [Google Scholar] [CrossRef]
  60. Lotz, W. Electron Binding Energies in Free Atoms*. J. Opt. Soc. Am. 1970, 60, 206–210. [Google Scholar] [CrossRef]
  61. Wang, M.; Huang, W.; Kondev, F.; Audi, G.; Naimi, S. The AME 2020 atomic mass evaluation (II). Tables, graphs and references*. Chin. Phys. C 2021, 45, 030003. [Google Scholar] [CrossRef]
  62. Barabash, A.S. Does the weak interaction constant depend on time? J. Exp. Theor. Phys. Lett. 1998, 68, 1–6. [Google Scholar] [CrossRef]
  63. Barabash, A.S. Possible evidence for a time variation of the weak-interaction constant from double-beta-decay experiments. Phys. At. Nucl. 2000, 63, 1210–1213. [Google Scholar] [CrossRef]
  64. Barabash, A.S. Possible Evidence of Time Variation of Weak Interaction Constant from Double Beta Decay Experiments. Astrophys. Space Sci. 2003, 283, 607–612. [Google Scholar] [CrossRef]
  65. Kotila, J.; Iachello, F. Phase-space factors for double-β decay. Phys. Rev. C 2012, 85, 034316. [Google Scholar] [CrossRef]
  66. Yan, X.; Cheng, Z.; Abdukerim, A.; Bo, Z.; Chen, W.; Chen, X.; Cheng, C.; Cui, X.; Fan, Y.; Fang, D.; et al. Searching for Two-Neutrino and Neutrinoless Double Beta Decay of 134Xe with the PandaX-4T Experiment. arXiv 2023, arXiv:2312.15632. [Google Scholar]
  67. Eliseev, S.; Roux, C.; Blaum, K.; Block, M.; Droese, C.; Herfurth, F.; Kluge, H.J.; Krivoruchenko, M.I.; Novikov, Y.N.; Minaya Ramirez, E.; et al. Resonant Enhancement of Neutrinoless Double-Electron Capture in 152Gd. Phys. Rev. Lett. 2011, 106, 052504. [Google Scholar] [CrossRef] [PubMed]
  68. Krivoruchenko, M.; Šimkovic, F.; Frekers, D.; Faessler, A. Resonance enhancement of neutrinoless double electron capture. Nucl. Phys. A 2011, 859, 140–171. [Google Scholar] [CrossRef]
  69. Blaum, K.; Eliseev, S.; Danevich, F.A.; Tretyak, V.I.; Kovalenko, S.; Krivoruchenko, M.I.; Novikov, Y.N.; Suhonen, J. Neutrinoless double-electron capture. Rev. Mod. Phys. 2020, 92, 045007. [Google Scholar] [CrossRef]
  70. Šimkovic, F.; Dvornický, R.; Vogel, P. Muon capture rates: Evaluation within the quasiparticle random phase approximation. Phys. Rev. C 2020, 102, 034301. [Google Scholar] [CrossRef]
Figure 1. A simplified schematic representation of the 2 ν ECEC process in   124 Xe. Top panels: both electrons are captured from K shell. Bottom panels: both electrons are captured from L shell. Left panels: two bound electrons from K (top) or L (bottom) shell are captured by the initial nucleus   124 Xe. Right panels: the atomic relaxation of the final neutral atomic system   124 Te, with two holes in K (top) or L (bottom) shell. The de-excitation is performed via X-ray emissions and Auger electrons from outer shells. The figure is reproduced from [52] by adding the schematic representation of the 2 ν L 1 L 1 process in   124 Xe.
Figure 1. A simplified schematic representation of the 2 ν ECEC process in   124 Xe. Top panels: both electrons are captured from K shell. Bottom panels: both electrons are captured from L shell. Left panels: two bound electrons from K (top) or L (bottom) shell are captured by the initial nucleus   124 Xe. Right panels: the atomic relaxation of the final neutral atomic system   124 Te, with two holes in K (top) or L (bottom) shell. The de-excitation is performed via X-ray emissions and Auger electrons from outer shells. The figure is reproduced from [52] by adding the schematic representation of the 2 ν L 1 L 1 process in   124 Xe.
Universe 10 00098 g001
Figure 2. Theoretical half-life scaled by g A 4 M 2 ν ECEC 2 as function of atomic number for nuclei presented in Table 3. The color scale indicates the Q-value for each transition.
Figure 2. Theoretical half-life scaled by g A 4 M 2 ν ECEC 2 as function of atomic number for nuclei presented in Table 3. The color scale indicates the Q-value for each transition.
Universe 10 00098 g002
Figure 3. Absolute values of the effective nuclear matrix elements for the measured nuclei undergoing 2 ν ECEC process or 2 ν β β decay. The values for 2 ν ECEC processes and 2 ν β β decay of   134 Xe are the ones from Table 5. The other 2 ν β β values are the ones reported in [20], computed using the PSFs from [65]. Vertical bars indicate uncertainties derived from the ones of the measured half-lives. Arrows indicate that the corresponding point is an upper limit for the effective nuclear matrix element.
Figure 3. Absolute values of the effective nuclear matrix elements for the measured nuclei undergoing 2 ν ECEC process or 2 ν β β decay. The values for 2 ν ECEC processes and 2 ν β β decay of   134 Xe are the ones from Table 5. The other 2 ν β β values are the ones reported in [20], computed using the PSFs from [65]. Vertical bars indicate uncertainties derived from the ones of the measured half-lives. Arrows indicate that the corresponding point is an upper limit for the effective nuclear matrix element.
Universe 10 00098 g003
Table 1. The binding energies, t x , in units of eV, for all occupied s 1 / 2 orbitals of the neutral atoms, Kr, Xe and Ba, undergoing 2 ν ECEC process. The experimental data are sourced from [60], our values are obtained with the DHFS self-consistent method, and values from prior studies are taken from [35] and [50,51].
Table 1. The binding energies, t x , in units of eV, for all occupied s 1 / 2 orbitals of the neutral atoms, Kr, Xe and Ba, undergoing 2 ν ECEC process. The experimental data are sourced from [60], our values are obtained with the DHFS self-consistent method, and values from prior studies are taken from [35] and [50,51].
Shell ( n j ) t x  [35] t x  [50,51] t x  [This Work] t x  [Experiment] [60]
(eV)(eV)(eV)(eV)
Kr
K ( 1 s 1 / 2 ) 17 , 936 17 , 700 14 , 280 14 , 327 ± 2
L 1 ( 2 s 1 / 2 ) 4507 3100 1902 1927 ± 2
M 1 ( 3 s 1 / 2 ) 278 292 ± 2
N 1 ( 4 s 1 / 2 ) 27.46 27.51 ± 0.1
Xe
K ( 1 s 1 / 2 ) 41 , 340 39 , 400 34 , 556 34 , 565 ± 2
L 1 ( 2 s 1 / 2 ) 10 , 424 7800 5417 5452 ± 2
M 1 ( 3 s 1 / 2 ) 1122 1149 ± 2
N 1 ( 4 s 1 / 2 ) 208 213 ± 2
O 1 ( 5 s 1 / 2 ) 23.63 23.40 ± 0.1
Ba
K ( 1 s 1 / 2 ) 44 , 610 42 , 400 37 , 450 37 , 442 ± 2
L 1 ( 2 s 1 / 2 ) 11 , 293 8500 5961 5991 ± 2
M 1 ( 3 s 1 / 2 ) 1274 1293 ± 2
N 1 ( 4 s 1 / 2 ) 256 254 ± 2
O 1 ( 5 s 1 / 2 ) 38 31 ± 2
P 1 ( 6 s 1 / 2 ) 4.64 5.21 ± 0.1
Table 2. Comparison of the G K L 1 2 ν ECEC (in units of 10 24 yr   1 ) values from [35,49,50,51], with values obtained through the method explained in this study. The Q-values cited in [49] were utilized. In cases where these values were not available, the Q was calculated using atomic masses provided in [61].
Table 2. Comparison of the G K L 1 2 ν ECEC (in units of 10 24 yr   1 ) values from [35,49,50,51], with values obtained through the method explained in this study. The Q-values cited in [49] were utilized. In cases where these values were not available, the Q was calculated using atomic masses provided in [61].
NucleusQ G K L 1 2 ν ECEC  [35] G K L 1 2 ν ECEC  [49] G K L 1 2 ν ECEC  [50,51] G K L 1 2 ν ECEC  [This Work]
(MeV) 10 24 yr 1 10 24 yr 1 10 24 yr 1 10 24 yr 1
  36 Ar0.43259 2.900 × 10 4 4.168 × 10 4
  40 Ca0.19351 1.250 × 10 5 1.020 × 10 5 1.314 × 10 5
  50 Cr1.1688 4.220 × 10 1 2.380 × 10 1 4.161 × 10 1
  54 Fe0.6798 4.690 × 10 2 3.021 × 10 2 4.553 × 10 2
  58 Ni1.926317.0015.309.90014.79
  64 Zn1.09481.4101.0301.364
  74 Se1.2091695.6563.4105.454
  78 Kr2.8463774.0660.0410.0637.0
  84 Sr1.7993.6064.6290.58
  92 Mo1.651208.082.32128.2
  96 Ru2.71451 2.740 × 10 3 2.400 × 10 3 1.450 × 10 3 2.328 × 10 3
  102 Pd1.172746.0042.0944.64
  106 Cd2.77539 6.220 × 10 3 5.410 × 10 3 4.299 × 10 3 5.269 × 10 3
  108 Cd0.27204 2.070 × 10 2 6.820 × 10 2 1.975 × 10 2
  112 Sn1.91982 1.150 × 10 3 869.7 1.120 × 10 3
  120 Te1.71481888.0840.3866.3
  124 Xe2.8654 2.020 × 10 4 1.720 × 10 4 1.510 × 10 4 1.685 × 10 4
  126 Xe0.9246.1060.5944.98
  130 Ba2.619 1.630 × 10 4 1.500 × 10 4 1.477 × 10 4 1.464 × 10 4
  132 Ba0.84439.1061.9838.12
  136 Ce2.37853 1.580 × 10 4 1.250 × 10 4 1.222 × 10 4 1.224 × 10 4
  138 Ce0.69818.4034.4717.92
  144 Sm1.78259 5.150 × 10 3 6.436 × 10 3 5.055 × 10 3
  152 Gd0.0557 1.120 × 10 2 5.989 × 10 7
  156 Dy2.012 1.760 × 10 4 2.208 × 10 4 1.734 × 10 4
  158 Dy0.284 1.830 × 10 1 3.191 1.751 × 10 1
  162 Er1.844 1.810 × 10 4 1.500 × 10 4 2.008 × 10 4 1.470 × 10 4
  164 Er0.02507 8.300 × 10 3 7.392 × 10 11
  168 Yb1.40927 4.710 × 10 3 7.872 × 10 3 4.647 × 10 3
  174 Hf1.0988 1.580 × 10 3 3.432 × 10 3 1.563 × 10 3
  180 W0.1432 1.560 × 10 3 1.478 1.321 × 10 3
  184 Os1.453 1.290 × 10 4 2.422 × 10 4 1.275 × 10 4
  190 Pt1.384 1.290 × 10 4 2.815 × 10 4 1.285 × 10 4
  196 Hg0.82821.0 3.587 × 10 3 815.8
Table 3. Values of the G K edge 2 ν ECEC (in units of 10 24 yr   1 ) and of the capture fractions for various shell pairs. The Q-values are computed using atomic masses provided in [61].
Table 3. Values of the G K edge 2 ν ECEC (in units of 10 24 yr   1 ) and of the capture fractions for various shell pairs. The Q-values are computed using atomic masses provided in [61].
NucleusQ G K edge 2 ν ECEC KK KL 1 KM 1 KN 1 L 1 L 1
(MeV) 10 24 yr 1 (%)(%)(%)(%)(%)
  36 Ar0.4326 4.244 × 10 4 83.4014.231.620.61
  40 Ca0.1935 1.347 × 10 5 81.3815.442.070.160.73
  50 Cr1.1705 4.301 × 10 1 81.3515.392.180.120.73
  54 Fe0.68076 4.717 × 10 2 80.4116.012.350.160.80
  58 Ni1.926415.2280.4515.972.370.150.79
  64 Zn1.095021.40779.7316.502.480.140.85
  74 Se1.209245.64878.6417.032.740.300.92
  78 Kr2.84767661.978.6416.922.810.350.91
  84 Sr1.7897794.1677.8317.332.980.440.96
  92 Mo1.65044133.876.6917.893.250.561.04
  96 Ru2.7145 2.437 × 10 3 76.6417.853.310.601.04
  102 Pd1.2034753.7075.1618.773.590.681.17
  106 Cd2.77539 5.540 × 10 3 75.7118.313.550.691.11
  108 Cd0.27179 2.122 × 10 2 65.3924.875.121.012.26
  112 Sn1.91981 1.181 × 10 3 74.7818.823.720.761.18
  120 Te1.73558975.874.0619.203.870.831.24
  124 Xe2.85674 1.756 × 10 4 74.2219.003.880.861.22
  126 Xe0.9177847.3671.6620.624.280.951.48
  130 Ba2.6237 1.570 × 10 4 73.5719.314.000.921.27
  132 Ba0.8440740.8770.4521.264.501.031.59
  136 Ce2.37853 1.303 × 10 4 72.8819.674.140.971.33
  138 Ce0.6959419.0368.3722.434.861.141.82
  144 Sm1.7824 5.417 × 10 3 71.1120.674.491.071.50
  150 Gd1.28728 1.301 × 10 3 69.2921.724.811.161.69
  152 Gd0.05567 1.547 × 10 6 7.40 × 10 4 1.04 × 10 3 38.56
  154 Dy3.31234 2.490 × 10 5 71.2920.484.511.091.47
  156 Dy2.00595 1.837 × 10 4 70.1921.144.691.141.59
  158 Dy0.28282 2.028 × 10 1 42.3135.469.042.266.01
  162 Er1.84696 1.601 × 10 4 69.2921.644.871.191.68
  164 Er0.02508 8.606 × 10 9 0.87
  168 Yb1.40936 5.052 × 10 3 67.4222.715.201.271.90
  174 Hf1.09994 1.724 × 10 3 64.9424.065.621.392.20
  180 W0.14323 2.522 × 10 3 5.94 × 10 5 17.728.082.3434.73
  184 Os1.45289 1.399 × 10 4 65.0223.905.641.432.17
  190 Pt1.40132 1.519 × 10 4 63.8324.505.841.502.32
  196 Hg0.81859916.256.4928.257.041.843.38
  212 Rn1.71019 1.400 × 10 5 61.3625.586.271.682.63
  214 Rn0.15031 1.376 × 10 2 1.071.380.5239.52
  218 Ra1.42814 6.758 × 10 4 58.5826.926.721.833.03
  224 Th1.16975 2.802 × 10 4 54.6828.747.342.033.64
  230 U0.75251 2.393 × 10 3 42.8033.669.162.595.87
  236 Pu0.4569395.1616.5139.2512.333.6413.66
  242 Cm0.08682 1.243 × 10 3 9.66
  252 Fm0.782 8.121 × 10 3 34.1136.0710.343.057.96
  258 No1.051 7.188 × 10 4 41.7733.619.282.736.14
Table 4. Values of the G K edge 2 ν ECEC (in units of 10 24 yr   1 ) and of the capture fractions for the dominant shell pairs in case of low Q-value transitions. The Q-values are computed using atomic masses provided in [61].
Table 4. Values of the G K edge 2 ν ECEC (in units of 10 24 yr   1 ) and of the capture fractions for the dominant shell pairs in case of low Q-value transitions. The Q-values are computed using atomic masses provided in [61].
NucleusQ G K edge 2 ν ECEC L 1 L 1 L 1 M 1 L 1 N 1 L 1 O 1 M 1 M 1 M 1 N 1 M 1 O 1 N 1 N 1 N 1 O 1
(MeV) 10 24 yr 1 (%)(%)(%)(%)(%)(%)(%)(%)(%)
  152 Gd0.05567 1.547 × 10 6 38.5635.9410.121.847.564.170.750.570.21
  164 Er0.02508 8.606 × 10 9 0.8726.5912.062.3328.1120.493.793.621.33
  242 Cm0.08682 1.243 × 10 3 9.6633.4913.914.1616.5512.493.652.321.35
Table 5. Values or limits of the experimental half-life M eff 2 ν ECEC . The Q-values are computed using atomic masses provided in [61]. The bottom row corresponds to the 2 ν β β decay of   134 Xe. For this case, the value of the effective NME was obtained from the lower half-life limit, and the PSF computed following [65] ( 0.225 × 10 22 yr   1 using Q = 0.8258 MeV).
Table 5. Values or limits of the experimental half-life M eff 2 ν ECEC . The Q-values are computed using atomic masses provided in [61]. The bottom row corresponds to the 2 ν β β decay of   134 Xe. For this case, the value of the effective NME was obtained from the lower half-life limit, and the PSF computed following [65] ( 0.225 × 10 22 yr   1 using Q = 0.8258 MeV).
Nucleus T 1 / 2 2 ν ECEC M eff 2 ν ECEC
(yr)
  78 Kr 9 . 2 2.9 + 5.7 × 10 21 [40] 0 . 457 0.098 + 0.095
1 . 9 0.8 + 1.3 × 10 22 [41] 0 . 318 0.073 + 0.100
  124 Xe ( 1.1 ± 0.2 ) × 10 22 [43] 0 . 072 0.006 + 0.008
  130 Ba ( 2.2 ± 0.5 ) × 10 21 [38] 0 . 170 0.017 + 0.023
( 6.0 ± 1.1 ) × 10 20 [39] 0 . 326 0.026 + 0.035
  132 Ba 2.2 × 10 21 [37]≤3.335
Nucleus T 1 / 2 2 ν β β M eff 2 ν β β
(yr)
  134 Xe 2.8 × 10 22 [66]≤0.398
Table 6. Summary of uncertainties in the PSFs due to input parameters.
Table 6. Summary of uncertainties in the PSFs due to input parameters.
Source δ G K edge 2 ν ECEC / G K edge 2 ν ECEC
Q-value(5–7) ×   δ Q / Q
E N 0.07 ×   δ E N / E N
Nuclear radius0.2 ×   δ R / R
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Niţescu, O.; Ghinescu, S.; Stoica, S.; Šimkovic, F. A Systematic Study of Two-Neutrino Double Electron Capture. Universe 2024, 10, 98. https://doi.org/10.3390/universe10020098

AMA Style

Niţescu O, Ghinescu S, Stoica S, Šimkovic F. A Systematic Study of Two-Neutrino Double Electron Capture. Universe. 2024; 10(2):98. https://doi.org/10.3390/universe10020098

Chicago/Turabian Style

Niţescu, Ovidiu, Stefan Ghinescu, Sabin Stoica, and Fedor Šimkovic. 2024. "A Systematic Study of Two-Neutrino Double Electron Capture" Universe 10, no. 2: 98. https://doi.org/10.3390/universe10020098

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