Neutrinoless Double-Electron Capture
REVIEWS OF MODERN PHYSICS, VOLUME 92, OCTOBER–DECEMBER 2020 Neutrinoless double-electron capture
K. Blaum and S. Eliseev Max-Planck-Institut f¨ur Kernphysik, Saupfercheckweg 1, 69117 Heidelberg, Germany
F. A. Danevich and V. I. Tretyak Institute for Nuclear Research of NASU, Prospekt Nauky 47, Kyiv 03028, Ukraine
Sergey Kovalenko Departamento de Ciencias Físicas, Universidad Andres Bello, Sazi´e 2212, Santiago 8320000, Chile
M. I. Krivoruchenko Institute for Theoretical and Experimental Physics, NRC “Kurchatov Institute”, B. Cheremushkinskaya 25, 117218 Moscow, Russia and National Research Centre “Kurchatov Institute”, Ploschad’ Akademika Kurchatova 1, 123182 Moscow, Russia
Yu. N. Novikov Petersburg Nuclear Physics Institute, NRC “Kurchatov Institute”, Gatchina, 188300 St. Petersburg, Russia and Department of Physics, St. Petersburg State University, 199034 St. Petersburg, Russia
J. Suhonen Department of Physics, University of Jyvaskyla, P.O. Box 35, Jyvaskyla FI-40014, Finland
(published 16 December 2020) Double-beta processes play a key role in the exploration of neutrino and weak interaction properties, and in the searches for effects beyond the standard model. During the last half century many attempts were undertaken to search for double-beta decay with emission of two electrons, especially for its neutrinoless mode 0ν2β−, the latter having still not been observed. Double-electron capture (2EC) was not yet in focus because of its in general lower transition probability. However, the rate of neutrinoless double-electron capture 0ν2EC can experience a resonance enhancement by many orders of magnitude when the initial and final states are energetically degenerate. In the resonant case, the sensitivity of the 0ν2EC process can approach the sensitivity of the 0ν2β− decay in the search for the Majorana mass of neutrinos, right-handed currents, and other new physics. An overview of the main experimental and theoretical results obtained during the last decade in this field is presented. The experimental part outlines search results of 2EC processes and measurements of the decay energies for possible resonant 0ν2EC transitions. An unprecedented precision in the determination of decay energies with Penning traps has allowed one to refine the values of the degeneracy parameter for all previously known near-resonant decays and has reduced the rather large uncertainties in the estimate of the 0ν2EC half-lives. The theoretical part contains an updated analysis of the electron shell effects and an overview of the nuclear-structure models, in which the nuclear matrix elements of the 0ν2EC decays are calculated. One can conclude that the decay probability of 0ν2EC can experience a significant enhancement in several nuclides.
DOI: 10.1103/RevModPhys.92.045007
CONTENTS III. Phenomenology of Neutrinoless 2EC 9 A. Underlying formalism 10 I. Introduction 2 B. Decay amplitude of the light Majorana neutrino II. Double-Electron Capture and Physics beyond exchange mechanism 10 the Standard Model 4 C. Comparison of 0ν2EC and 0ν2β− decay half-lives 12 0ν2 A. Quark-level mechanisms of EC 5 IV. Electron Shell Effects 13 B. Examples of underlying high-scale models 6 A. Interaction energy of electron holes 15 C. Hadronization of quark-level interactions 8 1. Electrostatic interaction 16
0034-6861=2020=92(4)=045007(61) 045007-1 © 2020 American Physical Society K. Blaum et al.: Neutrinoless double-electron capture
2. Retardation correction in the Feynman gauge 17 astrophysics, and cosmology. The fermions described by the 3. Magnetostatic interaction in the Feynman gauge 17 Dirac equation can also be electrically neutral. However, even 4. Magnetostatic interaction in the Coulomb gauge 18 in this case there is a conserved current in Diracs theory, which 5. Gauge invariance of the interaction energy ensures a constant number of particles (minus antiparticles). of electron holes 18 Majorana fermions, like photons, do not have such a con- B. Double-electron ionization potentials in Auger served current. In Dirac’s theory, particles and antiparticles are spectroscopy 19 independent, while Majorana fermions are their own anti- C. Summary 21 particles and, via CPT, it follows that they must have zero V. Nuclear Matrix Elements 22 charge. Truly neutral spin-1=2 fermions are referred to as A. Overview of the calculated nuclear matrix Majorana fermions, and those with a conserved current are 0ν2 elements in EC 22 referred to as Dirac fermions. Majorana fermions in bispinor B. Overview of the calculation frameworks 23 basis are vectors of the real vector space R4. They can also be 1. Multiple-commutator model 24 described by two-component Weyl spinors in the complex 2. Deformed quasiparticle random-phase space C2. In both representations, the superposition principle approximation 26 3. Microscopic interacting boson model 26 for Majorana fermions holds over the field of real numbers. 4. Energy-density functional method 26 Majorana fermions belong to the fundamental real represen- C. Decays of nuclides with the calculated nuclear tation of the Poincar´e group. A Dirac fermion of mass m can matrix elements 27 be represented as a superposition of two Majorana fermions of − D. Summary 28 masses m and m, respectively. VI. Status of Experimental Searches 28 The neutron has a nonvanishing magnetic moment, so it A. Experimental studies of 2EC processes 28 cannot be a pure Majorana particle. On the other hand, there B. Prospects for possible future experiments 42 is no fundamental reason to claim that it is a pure Dirac C. Neutrinoless 2EC with radioactive nuclides 43 particle. In theories with nonconservation of the baryon D. Summary 43 number, the mass eigenstates include a mixture of baryons and VII. Precise Determination of 2EC Decay Energies 43 antibaryons. At the phenomenological level, the effect is A. Basics of high-precision Penning-trap mass modeled by adding a Majoranian mass term to the spectrometry 43 effective Lagrangian. As a result, the neutron experiences B. Decay energies of 2EC transitions in virtually oscillations n ↔ n¯ whereby the nuclei decay with noncon- stable nuclides 45 servation of the baryon number (Dover, Gal, and Richard, C. Prospects for measurements of decay energies in 1983, 1985; Krivoruchenko, 1996a, 1996b; Gal, 2000; radioactive nuclides 47 Kopeliovich and Potashnikova, 2010; Phillips et al., 2016). VIII. Normalized Half-Lives of Near-Resonant Nuclides 47 Experimental limits for the period of the n ↔ n¯ oscillations A. Decays of virtually stable nuclides 47 τ 2 7 108 in the vacuum vac > . × s(Abe et al., 2015) constrain B. Decays of long-lived radionuclides 51 Δ ∼ 1 τ 0 8 10−33 the neutron Majorana mass to m = vac < . × m, IX. Conclusions 51 2 List of Symbols and Abbreviations 53 where m ¼ 939.57 MeV=c is the neutron Dirac mass. Thus, Acknowledgments 53 under the condition of nonconservation of the baryon number, References 54 Majorana’s idea on the existence of truly neutral fermions can be partially realized in relation to the neutron. In contrast, I. INTRODUCTION neutrinos can be pure Majorana fermions or pure Dirac fermions or a mixture of these two extreme cases. It is In 1937, Ettore Majorana found an evolution equation for a noteworthy that none of the variants of neutrino masses are truly neutral spin-1=2 fermion (Majorana, 1937). His work possible within the SM. The neutrino mass problem leads us to was motivated by the experimental observation of a free physics beyond the SM. Alternative examples of Majorana neutron by James Chadwick in 1932. Majorana also conjec- particles include weakly interacting dark-matter candidates tured that his equation applies to a hypothetical neutrino and Majorana zero modes in solid state systems (Elliott and introduced by Wofgang Pauli to explain the continuous Franz, 2015). − spectrum of electrons in the β decay of nuclei. In the mid- Searches for neutrinoless double-beta (0ν2β ) decay, neu- 1950s, Reines and Cowan (1956, 1959) discovered a particle trinoless double-electron capture (0ν2EC) by nuclei, and other with neutrino properties. At the time of the establishment of lepton number violating (LNV) processes provide the pos- the standard model (SM) of electroweak interactions, three sibility to shed light on the question of the nature of neutrinos families of neutrinos were known. While the neutron is a whether they are Majorana or Dirac particles. By virtue of the composite fermion consisting of quarks, neutrinos acquired black-box theorem (Schechter and Valle, 1982; Hirsch, the status of elementary particles with zero masses in the Kovalenko, and Schmidt, 2006), observation of the 0ν2β SM framework. The discovery of neutrino oscillations decay would prove that neutrinos have a finite Majorana mass. by the Super-Kamiokande Collaboration (Fukuda et al., The massive Majorana neutrinos lead to a violation of the 1998) showed that neutrinos are mixed and massive. A conservation of the total lepton number L. In the quark sector comprehensive description of the current state of neutrino of the SM, the baryon charge B is a similar quantum number. physics was given by Bilenky (2018). Vector currents of B and L are classically conserved. Left- Whether neutrinos are truly neutral fermions is one handed fermions are coupled to the SUð2Þ electroweak gauge of the fundamental questions in modern particle physics, fields W and Z0 so that vector currents of B and L,
Rev. Mod. Phys., Vol. 92, No. 4, October–December 2020 045007-2 K. Blaum et al.: Neutrinoless double-electron capture
FIG. 1. Schematic representation of neutrinoless double-beta FIG. 2. Schematic representation of neutrinoless double- decay. Two neutrons in the nucleus experience β decay accom- electron capture. Two protons in the nucleus each capture a panied by the exchange of a Majorana neutrino. Neutrons, bound electron from the electron shell and turn into two neutrons protons, electrons, and neutrinos are represented by solid lines. by the exchange of a Majorana neutrino. Notation is the same as Arrows show the flow of baryon charge of protons and neutrons in Fig. 1. and the flow of lepton charge of electrons and neutrinos. The cross denotes the Majorana neutrino mass term mL that causes the helicity flip of the intermediate neutrino and violates the lepton mechanism related to the Majorana neutrino exchange is number by two units. shown in Fig. 2. Subsequent deexcitation of the nucleus occurs via gamma-ray radiation or β decays. Deexcitation of as ’t Hooft (1976) first pointed out, are sensitive to the axial the electron shell is associated with the emission of Auger anomaly. Through electroweak instantons, this leads to non- electrons or gamma rays in a cascade formed by the filling of conservation of B and L, while the difference B − L is electron vacancies. Absent special selection rules, dipole conserved. The violating amplitude is exponentially sup- radiation dominates in x rays. Since the dipole moment pressed. Another example is given by sphaleron solutions of electrons is much higher than that of nucleons in the of classical field equations of the SM that preserve B þ L but nucleus, the deexcitation of the electron shell goes faster. For violate B and L individually, which can be relevant for atoms with a low value of Z, the Auger-electron emission is cosmological implications (White, 2016). The conservation more likely. With an increase in the atomic number, the of B − L within the SM is not supported by any fundamental radiation of x-ray photons becomes dominant. The deexcita- principles analogous to local gauge symmetry, so B and L can tion of high electron orbits is due to Auger-electron emission be broken beyond the SM explicitly. Experimental observa- for all Z. tion of the proton decay or n ↔ n¯ oscillations could prove Estimates show that the sensitivity of the 0ν2EC process to nonconservation of B, while observation of the 0ν2β− decay the Majorana neutrino mass is many orders of magnitude − or the neutrinoless double-electron capture (0ν2EC) could lower than that of the 0ν2β decay. Winter pointed out that prove nonconservation of L with constant B. Moreover, these degeneracy of the energies of the parent atom ðA; ZÞ and the processes are of interest for determining the absolute daughter atom ðA; Z − 2Þ gives rise to resonant enhance- neutrino mass scale, the type of neutrino mass hierarchy, ment of the decay. In the early 1980s, Georgi, Glashow, and and the character of CP violation in the lepton sector. Nussinov (1981) and Voloshin, Mitsel’makher, and Because of the exceptional value of LNV processes, a vast Éramzhyan (1982) also remarked on the possible resonant literature is devoted to physics of 0ν2β− decay and the enhancement of the 0ν2EC process. The resonances in 2EC underlined nuclear-structure models; see Bilenky and Petcov were considered, however, an unlikely coincidence. (1987), Suhonen (2007), Avignone, Elliott, and Engel To compensate for the low probability of the 0ν2EC process (2008), Vergados, Ejiri, and Šimkovic (2012), Raduta by a resonance effect, it is necessary to determine the energy (2015), Engel and Menx´endez (2017),andEjiri, Suhonen, difference of atoms with high accuracy. The decay probability Γ ðΔ2 þ Γ2 4Þ and Zuber (2019). is proportional to the Breit-Wigner factor f= f= , − The 0ν2β decay was first discussed by Furry (1939). where Γf is the electromagnetic decay width of the daughter The process is shown in Fig. 1. A nucleus with the mass atom and Δ ¼ MA;Z − M −2 is the degeneracy parameter − A;Z number A and charge Z experiences 0ν2β decay accom- equal to the mass difference of the parent and the daughter panied by the exchange of a Majorana neutrino between the atoms. The maximum increase in probability is achieved nucleons for Δ ¼ 0 when the decay amplitude approaches the unitary limit. Taking Δ ∼ 10 keV for the typical splitting of the ð Þ → ð þ 2Þþþ þ − þ − ð Þ A; Z A; Z e e ; 1:1 masses of the atoms and Γf ∼ 10 eV for the typical decay width of the excited electron shell of the daughter atom, one þþ where ðA; Z þ 2Þ is the doubly ionized atom in the final finds a maximum enhancement of ∼106. The degeneracy state. There are many models beyond the SM that provide parameter Δ ≲ Γ gives the half-life of a nuclide with respect − f alternative mechanisms of the 0ν2β decay, some of which are to 0ν2EC comparable to the half-life of nuclides with respect discussed in Sec. II. to 0ν2β− decay. In 1955, the related 0ν2EC process The near-resonant 0ν2EC process was analyzed in detail by Bernabeu, De Rujula, and Jarlskog (1983). They developed a − þ − þð Þ → ð − 2Þ ð Þ eb eb A; Z A; Z 1:2 nonrelativistic formalism of the resonant 0ν2EC in atoms and specified a dozen of nuclide pairs for which degeneracy is not − 0ν2 was discussed by Winter (1955a). Here eb are bound excluded. The EC process became the subject of a detailed electrons. The nucleus and the electron shell of the neutral theoretical study by Sujkowski and Wycech (2004). A list of atom ðA; Z − 2Þ are in excited states. An example of the the near-resonant 0ν2EC nuclide pairs was also provided by
Rev. Mod. Phys., Vol. 92, No. 4, October–December 2020 045007-3 K. Blaum et al.: Neutrinoless double-electron capture
Karpeshin (2008). The problem acquired an experimental 0ν2β− ones. At the same time, there is a motivation to search character: the difference between masses of the parent and for the neutrinoless ECβþ and 2βþ decays owing to the daughter atoms, i.e., Q values, known to an accuracy of about potential to clarify the possible contribution of the right- 10 keV, which is too far from the accuracy required to identify handed currents to the 0ν2β− decay rate (Hirsch et al., 1994) the unitary limit. The determination of the degeneracy and the appealing possibility of the resonant 0ν2EC proc- parameter has acquired fundamental importance. esses. The complicated effect signature expected in resonant In the 1980s, there was no well-developed technique to 0ν2EC transitions becomes an advantage: the detection of measure the masses of nuclides with relative uncertainty of several gamma quanta with well-known energies could be a about 10−9 sufficient to find resonantly enhanced 0ν2EC strong proof of the pursued effect. processes. The current state-of-the-art technique high- The previously mentioned aspects of the phase space, precision Penning-trap mass spectrometry was still in its degeneracy, abundance factors, etc., play an important role þ infancy. Its triumphal advance in the field of high-precision in determining the half-lives of the 0νECβ and 0ν2EC mass measurements on radioactive nuclides began with the processes. A further ingredient affecting the decay half-lives installation of the ISOLTRAP facility at CERN in the late are the involved nuclear matrix elements (NMEs); see 1980s (Bollen et al., 1987; Mukherjee et al., 2008; Kluge, Suhonen (2012a), Maalampi and Suhonen (2013),and 2013). Recent decades have been marked by a rising variety of Ejiri, Suhonen, and Zuber (2019). These NMEs have been high-precision Penning-trap facilities in Europe, the USA, and calculated in various nuclear-theory frameworks for a Canada (Blaum, 2006; Blaum, Dilling, and Nörtersh¨auser, number of nuclei. In this review we use these NMEs, as 2013). This led to a tremendous development of Penning-trap well as NMEs that have been calculated just for this review, mass-measurement techniques (Eliseev et al., 2007, 2013, to estimate the half-lives of those 0ν2EC transitions that are 2014; George, Blaum et al., 2007; Kretzschmar, 2007, 2013; of interest due to their possibly favorable resonance Blaum, Dilling, and Nörtersh¨auser, 2013) and made it possible conditions. to routinely carry out mass measurements on a broad variety of nuclides with a relative uncertainty of about 10−9. The mass II. DOUBLE-ELECTRON CAPTURE AND PHYSICS of the ion is determined via the measurement of its free BEYOND THE STANDARD MODEL cyclotron frequency in a pure magnetic field, the most 0ν2 precisely measurable quantity in physics. The underlying quark-level physics behind the EC process [see Eq. (1.2)] is basically the same as for the These factors motivated a new study of the near-resonant 0ν2β− 0ν2βþ 0ν βþ 0ν2EC process. A relativistic formalism for calculating , , and EC decays. In Figs. 1 and 2,we electron shell effects was developed and an updated realistic show the mechanism of exchange of light or heavy neutrinos with Majorana mass. The latter arise beyond the SM list of nuclide pairs for which the measurement of Q2EC values has high priority was compiled (Krivoruchenko et al., within the Weinberg dimension-5 effective LNV operator (Weinberg, 1979) providing conditions for the existence 2011). A significantly refreshed database of the nuclides of the processes 0ν2EC and 0ν2β−. A violation of the and their excited states is now available, 30 years after lepton number can also occur from quark-lepton effective the previous publication by Bernabeu, De Rujula, and Lagrangians of higher dimensions, corresponding to other Jarlskog (1983). An overview of the investigation of the possible mechanisms of the 0ν2EC process. The neutrinoless resonant 0ν2EC was given by Eliseev, Novikov, and Blaum 2EC can be accompanied by the emission of one or more (2012) including persistent experimental attempts to extremely light particles other than neutrinos in 2ν2EC. A search for appropriate candidates for this extraordinary well-known example is the Majoron J as the Goldstone phenomenon. boson of a spontaneously broken Uð1Þ symmetry of the The advancements of the experiments in search of the 0ν2EC L lepton number. Passing to the hadronic level one meets two process are more modest than those searching for 0ν2β− decay. − possibilities of hadronization of the quark-level underlying While the sensitivity of the 0ν2β experiments approaches − 24 26 process known from 0ν2β decay: direct nucleon and pionic half-life limits T1 2 ∼10 –10 yr, which constrains the effec- = mechanisms. Next, we will take a closer look at the tive Majorana neutrino mass of electron neutrino to previously mentioned aspects of the 0ν2EC process. j j ≲ 0 1–0 7 0ν2 mββ . . eV, the results of the best EC experiments First, the underlying quark-level mechanisms of the neu- ∼ 1019–1022 are still approximately T1=2 yr. The reasons for trinoless 2EC can be classified according to the following this difference are rather obvious: there is usually a much possible exotic final states: lower relative abundance of the isotopes of interest (typically • No exotic particles in the final state. The reaction 0ν2EC lower than 1%), and additionally a more complicated effect is shown in Eq. (1.2). • 0ν2 − þ − þð Þ→ð −2Þ þ signature due to the emission of a gamma-quanta cascade The reaction ECnJ: eb eb A;Z A;Z (instead of a clear 0ν2β− peak at the decay energy). The nJ, with n being the number of Majorons or Majoron- secondcircumstanceresults in alower detection efficiencyfor like exotic particles in the final state. the most energetic peak in a 0ν2EC energy spectrum. Both kinds of reactions can be further classified by the Furthermore, the energy of the most energetic 0ν2EC peak typical distance between particles involved in the under- is generally lower than in the 0ν2β− processes, yet the higher lying quark-lepton process, depending on the masses of the energy of a certain process the better the suppression the intermediate particles (P¨as et al., 1999, 2001; of the radioactive background. As a result the scale of the Prezeau, Ramsey-Musolf, and Vogel, 2003; Cirigliano 2EC experiments is substantially smaller than that of the et al., 2017a, 2018a), as illustrated in Fig. 3.
Rev. Mod. Phys., Vol. 92, No. 4, October–December 2020 045007-4 K. Blaum et al.: Neutrinoless double-electron capture
ð9ÞXY ¯ ¯ O1 ¼ 4ðdPXuÞðdPY uÞj; ð2:6Þ
ð9ÞXY ¯ μν ¯ O ¼ 4ðdσ P uÞðdσμνP uÞj; ð2:7Þ (a) (b) (c) 2 X X
ð9ÞXY ¯ μ ¯ FIG. 3. A decomposition of the generic ΔL ¼ 2 vertex into (a), O3 ¼ 4ðdγ PXuÞðdγμPY uÞj; ð2:8Þ (b) the long-range and (c) short-range quark-level contributions to 0ν2EC. (a) The conventional Majorana neutrino mass mecha- ð9ÞXY ¯ μ ν O ¼ 4ðdγ P uÞðu¯σμνP dÞj ; ð2:9Þ nism. The blobs in (b) and (c) denote the effective ΔL ¼ 2 4 X Y vertices. ð9ÞXY ¯ μ ¯ O5 ¼ 4ðdγ PXuÞðdPY uÞjμ; ð2:10Þ
• Long-range mechanisms with the Weinberg d ¼ 5 oper- where X; Y ¼ L; R and the leptonic currents are c c ator providing Majorana mass to the intermediate neu- j ¼ e ð1 γ5Þe, jμ ¼ e γμγ5e. The first term in Eq. (2.1) ¼ 6 trino in Fig. 3(a) and an effective d operator in the describes the SM low-energy four-fermion effective interac- upper vertex of Fig. 3(b). tion of the charged current (CC): • Short-range mechanisms with a dimension-9 effective operator in the vertex of Fig. 3(c). μ μ ¯ j ¼ νγ¯ ð1 − γ5Þe; J μ ¼ dγμð1 − γ5Þu: ð2:11Þ The effective operators in the low-energy limit origi- CC CC nate from diagrams with heavy exotic particles in the The SUð3Þ ×Uð1Þ symmetric operators in Eqs. (2.3)– internal lines. c em (2.10) are written in the mass-eigenstate basis. They originate The diagrams of the 0ν2ECnJ decays are derived from the SUð3Þ ×SUð2Þ ×Uð1Þ gauge invariant opera- from those in Fig. 3 by inserting one or more scalar c W Y Majoron lines either into the blobs of effective operators tors after the electroweak symmetry breaking; see Bonnet or into the central neutrino line of Fig. 3(a). et al. (2013), Lehman (2014), and Graesser (2017). Figures 3(a) and 3(b) are of second order in the Lagrangian (2.1). The effect of ΔL ¼ 2 is introduced in A. Quark-level mechanisms of 0ν2EC Figs. 3(a) and 3(b) by the Majorana neutrino mass term and by the d ¼ 6 effective operators (2.3)–(2.5), respectively. We now consider in more detail the short- and long-range Figure 3(a) is the conventional Majorana neutrino mass mechanisms of 0ν2EC. The corresponding diagrams for the mechanism with the contribution to the 0ν2EC amplitude 0ν2EC process are shown in Fig. 3. The blobs in Figs. 3(b) X and 3(c) represent the ΔL ¼ 2 effective vertices beyond the 2 Vαβ ∼ mββ ≡ U mν ; ð2:12Þ SM. At the low-energy scales μ ∼ 100 MeV typical for the ei i i 0ν2EC process, the blobs are essentially pointlike,having been generated by the exchange of a heavy particle with the where U is the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) M 0ν2 characteristic masses H much larger than the EC scale, mixing matrix and mββ is the effective electron neutrino ≫ μ i.e., MH . Integrating them out one finds the effective Majorana mass parameter that is well known from the analysis μ ≪ Lagrangian terms describing the vertices at the scale of 0ν2β− decay. The contribution of Fig. 3(b) is independent Λ ∼ MH for any kind of underlying high-scale physics beyond of the neutrino mass but proportional to the momentum q the SM. These vertices can be written in the following form flowing in the neutrino propagator. This is the so-called qˆ-type (P¨as et al., 1999, 2001): contribution. The following comment is in order. Figures 3(a) X and 3(b) show two possible mechanisms of both 0ν2EC and ð6Þ GF μ X ð6ÞX 0ν2β− L ¼ pffiffiffi −j J μ þ C ðμÞO ðμÞ ; ð2:1Þ (with the inverted final to initial states) processes. The ql 2 CC CC i i i first mechanism in Fig. 3(a) contributes to the amplitude of these processes with terms proportional to the effective 2 X ð9Þ G ð9ÞXY Majorana mass parameter mββ defined in Eq. (2.12). The L ¼ F CXYðμÞO ðμÞ: ð : Þ ql 2 i i 2 2 mp i;XY contribution of the second mechanism in Fig. 3(b) has no explicit dependence on mββ. This is because the upper vertex The first and second lines correspond to Figs. 3(b) and 3(c), in Fig. 3(b) breaks the lepton number into two units, as respectively. The proton mass mp is introduced to match the necessary for this process to proceed without need of the conventional notations. The complete set of the ΔL ¼ 2 ΔL ¼ 2 Majorana neutrino mass insertion into the neutrino operators for d ¼ 6 and 9 is as follows (Arbeláez et al., line. Note that mββ ¼ 0 is compatible with the neutrino 2016, 2017; González, Hirsch, and Kovalenko, 2016): oscillation data in the case of normal neutrino mass ordering. This result shows that the mechanism in Fig. 3(a), propor- ð6ÞX ¯ c O1 ¼ 4ðdPXuÞðν PLeÞ; ð2:3Þ tional to mββ, can be negligible in comparison to the ð6Þ mechanism in Fig. 3(b). Therefore, even if mββ turns out to O X ¼ 4ð¯σμν Þðνcσμν Þ ð Þ − 2 d PXu PLe ; 2:4 be small, both the 0ν2β decay and the 0ν2EC process can be
ð6ÞX μ observable due to the latter mechanism. This possibility has O ¼ 4ð¯γ Þðνcγ Þ ð Þ − 3 d μPXu PRe ; 2:5 been studied in the literature for 0ν2β decay; see P¨as et al.
Rev. Mod. Phys., Vol. 92, No. 4, October–December 2020 045007-5 K. Blaum et al.: Neutrinoless double-electron capture
Eqs. (2.5) and (2.6), leading to the contributions shown in Figs. 3(b) and 3(c). The blobs in Figs. 3(b) and 3(c) can be opened up (ultraviolet completed) in terms of all possible types of FIG. 4. The contribution of the Weinberg operator to the 0ν2EC renormalizable interactions consistent with the SM process in the flavor basis of the neutrino states. gauge invariance. These are the high-scale models, which lead to the 0ν2EC process. A list of all the possible 0ν2β− (1999, 2001), Arbeláez et al. (2016, 2017), and González, ultraviolet completions for decay was given by Hirsch, and Kovalenko (2016). Bonnet et al. (2013). The SM gauge invariant Weinberg dimension-5 effective The Wilson coefficients Ci in Eqs. (2.1) and (2.2) are operator generating the neutrino mass mechanism is given by calculable in terms of the parameters (couplings and masses) Λ ∼ (Weinberg, 1979) of a particular underlying model at the scale MH, called the “matching scale.” Note that some of CiðΛÞ may vanish. c 0ν2 ð5Þ ðL HÞðLHÞ κ To make contact with EC one needs to estimate Ci at a L ¼ κ ¼ νcν 0 0 þ μ 0ν2 W Λ Λ H H scale 0 close to the typical EC-energy scale. The 0 2 coefficients Ci run from the scale Λ down to μ0 due to SSB hH i ⟶κ νcν þ ð Þ the QCD corrections. In addition, the d ¼ 9 operators Λ : 2:13 undergo the renormalization group equation mixing with In Eq. (2.13) we mean LH ¼ L H ϵ , the singlet combination each other, leading to the mixing of the corresponding i j ij Wilson coefficients. ð2Þ LðLcÞ H of two SU W doublets and . After the electroweak The general parametrization of the 0ν2EC amplitude derived spontaneous symmetry breaking (SSB) with the Higgs vac- from Fig. 3, taking into account the leading-order QCD running uum expectation value hHi ≠ 0 neutrinos acquire a Majorana ¼ −2κh i2 Λ κ (González, Hirsch, and Kovalenko, 2016; Cirigliano et al., mass mν H = , with a dimensionless parameter. In 2018a; Ayala, Cvetic, and Gonzalez, 2020; Liao, Ma, and the flavor basis of neutrino states the contribution of the Wang, 2020), reads Weinberg operator to a LNV process such as 0ν2EC is displayed in Fig. 4. The summation of multiple insertions X3 of the Weinberg operator into the bare neutrino propagator 2 2 X X Vαβ ¼ G cos θ K β ðμ0; ΛÞC ðΛÞ entails the renormalized neutrino propagator with Majorana F C Z i i i¼1 mass mν. The operator (2.13) is unique. Other operators of X5 the effective Lagrangian are suppressed by higher powers of þ βXYðμ ΛÞ XYðΛÞ A ð Þ Λ i 0; Ci αβ: 2:14 the unification scale . The study of the neutrinoless 2EC i¼1 process and 2β− decays could be the most direct way of testing physics beyond the standard model. In terms of the The parameters βX and βXY incorporate the QCD running naive dimensional counting one can expect the dominance i i of the Wilson coefficients and the matrix elements of the of the Weinberg operator, with the dimension 5, over the μ operators in Eqs. (2.3)–(2.5) combined with j J μ and operators of dimensions 6 and 9 in Eqs. (2.5) and (2.6). CC CC the operators in Eqs. (2.6)–(2.10). The wave functions However, to set mν at eV scale one should provide a small κ ∼ 10−11 of the captured electrons with quantum numbers α and β coupling for the phenomenologically interesting A – case of Λ ∼ Oð1 TeVÞ. However, the smallness of any enter the coefficients αβ defined by Eqs. (2.20) (2.23).In dimensionless coupling requires explanation. Typically in Eq. (2.14) the summation over the different chiralities ¼ ðΛÞ this case one expects the presence of some underlying X; Y L; R is implied. The Wilson coefficients Ci physics, for example, symmetry. The situation changes with entering Eqs. (2.1) and (2.2) are linked to the matching Λ the increase of the LNV scale up to Λ ∼ 1013–14,withκ ∼ 1, scale , where they are calculable in terms of the where the contribution of the Weinberg operator to 0ν2EC Lagrangian parameters of a particular high-scale under- dominates. The final count depends on the concrete high- lying model. The decay amplitude (2.14) is supplemented scale underlying LNV model: not all operators appear in the by the overlap amplitude KZ of the electron shells of the low-energy limit and κ is a small suppression factor allowing initial and final atoms. In this review, we mainly discuss TeV-scale Λ. The latter can stem from loops or the ratio of the the light Majorana neutrino exchange mechanism of 0ν2 SSB scales in multiscale models; for a recent analysis see Fig. 3(a).The EC NMEs are currently known only for Helo, Hirsch, and Ota (2016). the Majorana neutrino exchange mechanisms coupled to In this review the mechanism of the neutrino Majorana left- and right-handed currents. Calculations of the NMEs masses is discussed in detail, for which numerical evaluation corresponding to the other long- and short-range mecha- of the neutrinoless 2EC half-lives of near-resonant nuclides nisms of Figs. 3(b) and 3(c), respectively, for all oper- with the known NMEs will be given. In the case of high- ators (2.5) and (2.6) are still in progress. dimensional operators, as well as for the d ¼ 5 mechanism with the unknown NMEs, normalized estimates are given that B. Examples of underlying high-scale models take into account the factorization of nuclear effects in the 0ν2EC amplitude. Keeping the previous comments in mind, We give three examples of popular high-scale models that we also discuss mechanisms based on the operators of can underlie the 0ν2EC process. In the low-energy limit their
Rev. Mod. Phys., Vol. 92, No. 4, October–December 2020 045007-6 K. Blaum et al.: Neutrinoless double-electron capture contribution is described by the effective Lagrangians (2.1) and/or (2.2). Left-right symmetric models.—A well-known example of a high-scale model leading to ΔL ¼ 2 processes, such as 0ν2β (a) (b) decay and 0ν2EC and the generation of Majorana mass for neutrinos, is the left-right symmetric extension of the SM. The left-right symmetric model (LRSM) is based on the gauge group G spontaneously broken via the chain
G ¼ ð3Þ ð2Þ ð2Þ ð1Þ (c) (d) (e) SU C ×SU L ×SU R ×U B−L ⇓ vR FIG. 5. Possible flavor-basis contributions to 0ν2EC within the ð3Þ ð2Þ ð1Þ LRSM. SU C ×SU L ×U Y ⇓ vSM SUð3Þ ×Uð1Þ ; ð2:15Þ mL mD C em Mν ¼ ; ð2:17Þ mT m ≡ hΔ i ≫ hΦi ≡ D R where vR R vSM are the vacuum expectation ð2Þ Δ values (VEVs) of a Higgs SU R triplet R and a Higgs with m ∼ y v and m ∼ yΦv 3 × 3 block Φ L;R L;R L;R D SM bidoublet , respectively. The bidoublet belongs to the matrices in the generation space. The matrix (2.17) is doublet representation of both SUð2Þ and SUð2Þ . There is T ν L R diagonalized to U M U ¼ diagðmν ; m Þ by an orthogonal ð2Þ Δ ≡ hΔ i i Ni also a Higgs SU L triplet L with the VEV vL L . mixing matrix Left- and right-handed leptons and quarks belong to the ð2Þ ð2Þ doublet representations of the SU L and SU R gauge ð3Þ ð2Þ ð2Þ UL UD groups, respectively. The SU C ×SU L ×SU R × U ¼ ; ð2:18Þ ð1Þ T U B−L assignments of the LRSM fields are UD UR ν i with UL;R;D 3 × 3 block matrices in the generation space. The L ð Þ ¼ ∼ ½1; 2ð1Þ; 1ð2Þ; −1 ; L R l− ν i LðRÞ neutrino mass spectrum consists of three light 1;2;3 and three heavy N1 2 3 Majorana neutrino states with the masses mν u ; ; 1;2;3 ¼ i ∼ ½3 2ð1Þ 2ð1Þ; −1 3 QLðRÞ ; ; = ; and mN1 2 3 , respectively. d ; ; i LðRÞ The possible contributions to 0ν2EC within the LRSM are þ ! Δpffiffi Δþþ shown in Fig. 5. Δ ¼ 2 ∼ ½1 3ð1Þ 1ð3Þ;2 LðRÞ þ ; ; ; Figure 5(a) shows the conventional long-range light Δ0 −pΔffiffi 2 LðRÞ Majorana neutrino exchange mechanisms with the con- 0 þ tribution shown in Eq. (2.12). Figures 5(b) and 5(e) are Φ1 Φ1 Φ ¼ ∼ ½1; 2; 2; 0 ; ð2:16Þ short-range mechanisms with two heavy right-handed Φ− Φ0 2 2 bosons WR and heavy neutrino νR or doubly charged Higgs Δþþ exchange. In the low-energy limit they reduce where i ¼ 1; 2; 3 is the generation index. Previously introduced ð9ÞRR to the effective operators O3 in Eq. (2.8) depicted VEVs are related to the VEVs of the electrically neutral in Fig. 3(c). Figures 5(c) and 5(d), containing light hΔ i ≡ hΔ0 i hΦi2 ¼hΦ0i2 þhΦ0i2 ≡ components L;R L;R , 1 2 vSM. virtual neutrinos, represent the long-range mechanism There are two charged gauge bosons WL;R and two neutral of 0ν2EC. gauge bosons ZL;R with masses of the order of In the low-energy limit the upper parts of Figs. 5(c) and 5(d) ∼ ∼ ν ¼ 6 MWR ;MZR gRvR, MWL ;MZL gLvSM. Note that in the sce- with heavy particles WR and R reduce to the d effective ð2Þ ð6ÞR ð6ÞL nario with the manifest left-right symmetry the SU L;R gauge operators O3 and O3 , respectively. Note that these ¼ couplings obey gL gR. Since the bosons WR, ZR have not contributions to the 0ν2EC amplitude depend not on the light been experimentally observed, the scale of the left-right neutrino mass mν but on its momentum q flowing in the symmetry breaking vR must be sufficiently large, above a neutrino propagator. Technically this happens because differ- “ ” few TeV. On the other hand, the VEV of the left triplet vL ent chiralities of the lepton vertices project the qˆ term out of ρ ¼ 1 ˆ ˆ must be small since it affects the SM relation , which is in the neutrino propagator PLðq þ mνÞPR ¼ PLq. On the con- good agreement with the experimental measurements setting trary, Fig. 5(a), with the same chiralities in both vertices, is ≲ 8 an upper limit vL GeV. From the scalar potential of the proportional to mν due to PLðqˆ þ mνÞPL ¼ PLmν. This is ∼ 2 LRSM follows vL vSM=vR, which satisfies the previously consistent with the fact that in the latter case the source mentioned upper limit for vR ≳ 10 TeV. of LNV is the Majorana neutrino mass mν and in the limit The spontaneous gauge symmetry breaking (2.15) gener- mν → 0 the corresponding contribution must vanish. On the ates a 6 × 6 neutrino seesaw-I mass matrix given in the basis other hand, in the former case [Figs. 5(c) and 5(d)] the LVN ðν νCÞT L; R by source is the operator in the upper vertex and mν is not needed
Rev. Mod. Phys., Vol. 92, No. 4, October–December 2020 045007-7 K. Blaum et al.: Neutrinoless double-electron capture to allow for the ΔL ¼ 2 process to proceed. These are the so-called qˆ-type contributions. K The Wilson coefficients Cn in Eqs. (2.2) and (2.14) at the matching high-energy scale Λ ∼ MR corresponding to Figs. 5(b)–5(e) are given by (a) (b) (c) X3 m M 4 FIG. 6. Possible flavor-basis contributions to 0ν2EC within 5ð Þ ∼ hΔ i → RR;NR ¼ 2 p WL Fig: b yR R C3 U ; (a),(b) leptoquark and (a),(c) RPV MSSM models. The cross in Leim M i¼1 Ni WR (c) denotes the Majorana mass of a gluino g˜ or neutralino χ. ð2:19Þ X3 M 2 ¼ λ0 ¯ þ ϵ ð Þ R;qˆ WL WRPV ijkLiQjDk iLiHU; 2:23 Fig:5ðcÞ ∼ yΦy hΦihΔ i → C ¼ U U ; R R 3 Lei D;ei M i¼1 WR where Q, D, and H conventionally denote here the chiral ð Þ U 2:20 superfields of the left-chiral electroweak doublet quarks, the right-chiral electroweak singlet down quark, and the up-type X3 L;qˆ electroweak Higgs doublet, respectively. Fig:5ðdÞ ∼ ζ yΦy hΦihΔ i → C ¼ U U tan ζ ; W R R 3 Lei D;ei W The RPV SUSY models with the interactions (2.23) con- i¼1 tribute to the long- and short-range mechanisms of the 0ν2EC ð Þ 2:21 process. The corresponding diagrams are shown in Fig. 6. Figures 6(a) and 6(b) generate the long-range qˆ-type con- X3 4 Δþþ m m M ˜ χ 2 RR; R 2 Ni p WL tribution, while Fig. 6(c) with the gluino g or neutralino Fig:5ðeÞ ∼ λΔ g hΔ i → C3 ¼ U : R R R Rei m2 M exchange entails the short-range contribution. The source i¼1 ΔR WR of the lepton number violation is located in the vertices of ð Þ 2:22 Figs. 6(a) and 6(b), while in Fig. 6(c) it is given by the Majorana mass of the neutralino mχ and/or gluino m˜ . In the Here ζ is the angle of W -W mixing. For convenience we g W L R Oð6ÞL show the correspondence of the flavor-basis diagrams in Fig. 5 low-energy limit Figs. 6(a) and 6(b) lead to the operator 1 , to the particular combinations of the parameters of the LRSM while in this limit Fig. 6(c), where all internal particles are heavy, collapses to a pointlike short-range contribution given Lagrangian (quartic ΔR coupling λΔ , gauge coupling gR, and R ¼ 9 Oð9ÞLL Oð9ÞLL VEV hΔRi), then give the corresponding Wilson coefficients. by a linear combination of d operators 1 and 2 . Leptoquark models.—Leptoquarks (LQs) are exotic scalar or vector particles coupled to lepton-quark pairs in such a way C. Hadronization of quark-level interactions that L¯ ðLQÞQ. They appear in various high-scale contexts, for example, in grand unification, extended technicolor, compos- We now comment on the calculation of the structure βX βXY iteness, etc. For a generic LQ theory all the renormalizable coefficients i , i in the amplitude (2.14) depending on interactions were specified by Buchm¨uller, R¨uckl, and Wyler the NMEs and nucleon structure. The Lagrangians (2.1) (1987). Current experimental limits (Tanabashi et al., 2018) and (2.2) can, in principle, be applied to any LNV processes allow them to be relatively light at the TeV scale. The SM with whatever hadronic states: quarks, mesons, nucleons, and gauge symmetry allows LQs to mix with the SM Higgs. This other baryons, as well as nuclei. The corresponding amplitude, mixing generates ΔL ¼ 2 interactions with the chiral structure such as that in Eq. (2.14), involves the hadronic matrix leading to the long-range qˆ-type contribution not suppressed elements of the operators (2.6) and (2.5). The Wilson coefficients CXY are calculated in terms of the parameters by the smallness of the Majorana mass mν of the light virtual i neutrino displayed in Fig. 6(a), with S or V being scalar or of the high-scale model and are independent of the low-energy vector LQs. In the low-energy limit, the upper part of Fig. 6(a) scale nonperturbative hadronic dynamics. This is the cel- with heavy LQs reduces to the pointlike vertex described by ebrated property of the operator product expansion, express- ing interactions of some high-scale renormalizable model in Oð6ÞX the operator 1 in Eq. (2.3). The chirality structure of this the form of Eqs. (2.1) and (2.2) below a certain scale μ. In the vertex combined with the SM vertex in the bottom part of case of 0ν2β− decay, 0ν2EC, and other similar nuclear ˆ 0ν2 Fig. 6(a) render a q-type contribution to EC. processes, the corresponding NMEs of the operators (2.6) — R-parity violating supersymmetric models. The TeV-scale and (2.5) are calculated in the framework of the approach supersymmetric (SUSY) models offer a natural explanation of based on a nonrelativistic impulse approximation; for a – the grand unified theory SM scale hierarchy, introducing detailed description see Doi, Kotani, and Takasugi (1985). superpartners to each SM particle, so that they form super- This implies, as the first step, reformulation of the quark-level ˜ multiplets (chiral superfields): ðq; q˜Þ, ðl; lÞ, ðg; g˜Þ, etc. Here q˜ theory in terms of the nucleon degrees of freedom, which the ˜ and l are scalar squarks and sleptons, while g˜ is a spin-1=2 existing nuclear-structure models operates with. This is the so- gluino. The SUSY framework requires at least two electro- called hadronization procedure. In the absence of a firm theory weak Higgs doublets HU and HD. A class of SUSY models, of hadronization one is left to resort upon general principles the so-called R-parity violating (RPV) SUSY models, allow and particular models. Imbedding two initial (final) quarks for LNV interactions described by a superpotential into two different protons (neutrons) is conceptually a more
Rev. Mod. Phys., Vol. 92, No. 4, October–December 2020 045007-8 K. Blaum et al.: Neutrinoless double-electron capture
nucleon-nucleon interactions should be introduced for theo- retical self-consistency even in the case of the long-range light neutrino exchange mechanism in Fig. 3(a). For a detailed description of this approach see Prezeau, Ramsey-Musolf, and Vogel (2003), Cirigliano et al. (2017a, 2017b, 2018a, 2018b, (a) (b) (c) 2019, 2020), and Graesser (2017). To conclude, neutrinoless double-electron capture 0ν2EC, 0ν2 FIG. 7. Hadron-level diagrams for EC. (a) Conventional like 0νββ decay, is a ΔL ¼ 2 lepton number violating process. two-nucleon mechanism. (b) One-pion exchange mechanism. Moreover, on the level of a nucleon subprocess it is virtually (c) Two-pion exchange mechanism. equal to 0νββ decay. Consequently, the underlying ΔL ¼ 2 physics driving both these processes is the same. There are simple option that is illustrated in Fig. 7(a). This is the many formal differences in the form of the effective operators conventional two-nucleon mechanism relying on the nucleon representing this physics at low-energy sub-GeV scales. We 0ν2 form factors as a phenomenological representation of the specified a complete basis of the EC effective operators in – nucleon structure. On the other hand, putting one initial and Eqs. (2.1) (2.10) and exemplified high-energy scale models one final quark into a charged pion while the other initial presently popular in the literature, which can be reduced to these operators in the low-energy limit. Akin to 0νββ decay quark is put into a proton and one final quark into a neutron, as 0ν2 in Fig. 7(b), we deal with the one-pion mechanism. The two- there are basically three types of mechanisms of EC shown pion mechanism, displayed in Fig. 7(c), treats all quarks to be in Fig. 3: Fig. 3(a) shows the conventional neutrino exchange mechanism with the amplitude proportional to the effective incorporated in two charged pions. In both cases the pions are Majorana neutrino mass mββ defined in Eq. (2.12), and virtual and interact with nucleons via the ordinary pseudo- ¯ Fig. 3(b) displays a neutrino exchange mechanism indepen- scalar pion-nucleon coupling Niγ5τπN. One may expect dent of the Majorana neutrino mass, when the lepton number a priori dominance of the pion mechanism for the reason violation necessary for 0ν2EC to proceed is gained from a that it extends the region of the nucleon-nucleon interaction ΔL ¼ 2 vertex. Figures 3(a) and 3(b) are both long-range due to the smallness of the pion mass leading to a long-range mechanics induced by the exchange of an extremely light potential. As a result, the suppression caused by the short- particle, a neutrino. On the other hand, Fig. 3(c) represents a range nuclear correlation can be significantly alleviated in short-range mechanism induced by the exchange of heavy comparison to the conventional two-nucleon mechanism. particles with masses much larger than the typical scale Nevertheless, one should consider all these mechanisms to (approximately a few MeV) of 0ν2EC. Despite the underlying contribute to the process in question with corresponding physics of both 0ν2EC and 0νββ decay being the same, their relative amplitudes. The latter is as yet unknown. In principle, NMEs are significantly different. We discuss the nuclear- it can be evaluated in particular hadronic models. These kinds structure aspects and atomic physics involved in the calcu- of studies are still missing in the literature and consensus on lations of the 0ν2EC NMEs in subsequent sections. Here it is the dominance of one of the two mechanisms is pending. It is a worth noting that thus far only the NMEs for the Majorana common trend to posit an analysis of one of the two neutrino exchange mechanism in Fig. 3(a) have been calcu- hadronization scenarios. Note that for the long-range contri- lated in the literature. Similar calculations for NMEs of other butions described by the effective Lagrangian (2.1) the mechanisms in Figs. 3(b) and 3(c) are still pending. previously mentioned advantage of the pion mechanism is absent and one can, in a sense, safely resort to the conven- III. PHENOMENOLOGY OF NEUTRINOLESS 2EC tional two-nucleon mechanism. The light Majorana exchange contribution to 0ν2EC, on which we focus in the rest of this Figures 1 and 2 can be combined as shown in Fig. 8. In the review, is of this kind. This limitation is explained by the fact initial state, there is an atom ðA; ZÞ. The electron lines belong that for the moment there are no nuclear matrix elements yet to the electron shells, and the proton and neutron lines belong calculated in the literature for mechanisms other than this one. to the initial and intermediate nuclei, respectively. As a result The procedure of hadronization is essentially the same as of neutrinoless double-electron capture, an atom ðA; Z − 2Þ for 0νββ decay and is described in the literature. For more is formed, generally in an excited state. In what follows, details on this approach to hadronization see Doi, Kotani, and ðA; ZÞ denotes an atom with the excited electron shell and Takasugi (1985), Faessler, Kovalenko, and Šimkovic (1998), ðA; ZÞ indicates that the nucleus is also excited. The Faessler et al. (2008), and Graf et al. (2018). intermediate atom can decay by emitting a photon or Recently another approach has been developed, which Auger electrons, but it can also experience 0ν2β− transition resorts to matching the high-scale quark-level theory to chiral and evolve back to the initial state. As a result, LNV perturbation theory. The latter is believed to provide a low- oscillations ðA; ZÞ ↔ ðA; Z − 2Þ occur in the two-level energy description of QCD in terms of nucleon and pion system. These oscillations are affected by the coupling of degrees of freedom. It is expected that the parameters of the the ðA; Z − 2Þ atom to the continuum, which eventually low-energy effective theory can be determined from exper- leads to the decay of ðA; ZÞ. The Hamiltonian of the system is imental measurements or from the lattice QCD. This approach not Hermitian because ðA; Z − 2Þ has a finite width. leads to a different picture of hadronization and numerical The LNVoscillations of atoms were discussed by Šimkovic results compared to the previously sketched conventional and Krivoruchenko (2009), Krivoruchenko et al. (2011), and approach. Contrary to the conventional approach short-range Bernab´eu and Segarra (2018). The formalism of LNV
Rev. Mod. Phys., Vol. 92, No. 4, October–December 2020 045007-9 K. Blaum et al.: Neutrinoless double-electron capture qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 2 2 Ω ¼ Vαβ þ M−. The values of λ are complex, so the norm of the states is not preserved in time. A series expansion around V ¼ 0 yields
λ ≈ þ Δ − i Γ ð Þ þ Mi M 2 i; 3:4 FIG. 8. Oscillations of atoms induced by 0ν2EC and 0ν2β transitions. The notations are the same as in Figs. 1 and 2. The i λ− ≈ M − ΔM − ðΓ − Γ Þ; ð3:5Þ intermediate atom ðA; Z − 2Þ is coupled to the continuum f 2 f i through the emission of a photon or Auger electrons. These Γ 2 channels generate a finite width f in Eq. (3.1), and they are also with ΔM ¼ðMi − MfÞΓi=Γf, Γi ¼ VαβRf, and responsible for the non-Hermitian character of Heff . Γ ¼ f ð Þ Rf ð − Þ2 þð1 4ÞΓ2 : 3:6 oscillations allows one to find a relationship between the half- Mi Mf = f life T1 2 of the initial atom ðA; ZÞ, the amplitude of neutrino- = Γ ≪ Γ Γ less double-electron capture Vαβ, and the decay width of the The initial state decays at the rate i f. The width i is intermediate atom ðA; Z − 2Þ , which has an electromagnetic maximal for complete degeneracy of the atomic masses: origin. 4V2 Γmax ¼ αβ ð Þ i Γ : 3:7 A. Underlying formalism f
The evolution of a system of mixed states, each of A simple calculation gives which may be unstable due to the coupling with the con- − tinuum, can described by an effective non-Hermitian Heff Mþ UðtÞ¼expð−iMþtÞ cosðΩtÞ − i sinðΩtÞ : ð3:8Þ Hamiltonian (Weisskopf and Wigner, 1930). In the case under Ω consideration, the Hamiltonian takes the form The decay widths of single-hole excitations of atoms are M Vαβ known experimentally and were tabulated for 10 ≤ Z ≤ 92 ¼ i ð Þ Heff ; 3:1 and principal quantum numbers 1 ≤ n ≤ 4 by Campbell and V M − i Γ αβ f 2 f Papp (2001). The width of a two-hole state αβ is represented by the sum of the widths of the one-hole states Γ ¼ Γα þ Γβ. where M and M are the masses of the initial and final atoms. f i f The deexcitation width of the daughter nucleus is much The width Γ of the final excited atom with two vacancies α f smaller and can be neglected. The values Γ are used in and β is of electromagnetic origin. The off-diagonal matrix f estimating the decay rates Γ . elements are due to a violation of lepton number conservation. i The transition amplitude from the initial to the final state for They can be chosen as real by changing the phase of one of the small time t, according to Eq. (3.8), is equal to states; thus, we set Vαβ ¼ Vαβ. The real and imaginary parts of the Hamiltonian do not commute. hfjUðtÞjii¼−iVαβt þ . ð3:9Þ We now find the evolution operator
Equation (3.9) is valid for t ≲ 1=jMþj and also over a wider ð Þ¼ ð− Þ ð Þ U t exp iHeff t : 3:2 range t ≲ 1=jΩj, given that the real part of the phase can be made to vanish via redefinition of the Hamiltonian According to Sylvester’s theorem, the function of a finite- → − ℜð Þ Heff Heff Mþ . The value of Vαβ can be evaluated by dimensional n × n matrix A is expressed in terms of the means of field-theoretical methods that allow one to find the eigenvalues λk of the matrix A, which are solutions of the ð − λÞ¼0 amplitude (3.9) from first principles. The formalism described characteristic equation det A , and a polynomial of A: in this section reproduces the results of Bernabeu, De Rujula, 0ν2 X Y λ − and Jarlskog (1983) with respect to EC decay rates. ð− Þ¼ ð− λ Þ l A ð Þ f iAt f i kt λ − λ ; 3:3 ≠ l k k l k B. Decay amplitude of the light Majorana neutrino exchange mechanism where the sum runs over 1 ≤ k ≤ n, the product runs over 1 ≤ l ≤ n, l ≠ k, and the eigenvalues are assumed to be The total lepton number violation is due to the Majorana pairwise distinct. The matrix function fð−iAtÞ evolves with masses of the neutrinos. It is assumed that the left electron ð− λ Þ the time t like the superposition of n terms f i kt with the neutrino is a superposition of three left Majorana neutrinos: matrix coefficients that are projection operators onto the kth eigenstates of A. X3 The eigenvalues of the Hamiltonian (3.1) are equal νeL ¼ UekχkL; ð3:10Þ ¼1 to λ ¼ Mþ Ω, where M ¼ðMi MfÞ=2 ∓ iΓf=4 and i
Rev. Mod. Phys., Vol. 92, No. 4, October–December 2020 045007-10 K. Blaum et al.: Neutrinoless double-electron capture
where U is the PMNS mixing matrix. In the Majorana bispinor If the captured electrons occupy the states α ≡ ðn2jlÞ1 and χ ¼ð1 2Þð1 − γ Þχ χ ¼ χ β ≡ ð 2 Þ representation, kL = 5 k and k k. The ver- n jl 2, we must take the superposition of products of tex describing the creation and annihilation of a neutrino has their wave functions: the standard form X JM JM ψ ðr1; r2Þ¼ C Ψα ðr1ÞΨβ ðr2Þ; ð3:15Þ θ αβ jαmαjβmβ mα mβ GF cos C μ mαmβ HðxÞ¼ pffiffiffi jμðxÞJ ðxÞþH:c:; ð3:11Þ 2 h where jα and jβ are the total angular momenta, mα and mβ are where θ is the Cabibbo angle. The lepton and quark charged Ψ ðrÞ C their projections on the direction of the z axis, and αmα currents are defined by Eq. (2.11). In terms of the composite Ψ ðrÞ and βmβ are the relativistic wave functions of the bound fields, the hadron charged current is given by electrons in an electrostatic mean field of the nucleus and the surrounding electrons. The identity of the fermions implies μð Þ¼¯ð Þγμð − γ Þ ð Þ ð Þ Jh x n x gV gA 5 p x ; 3:12 that the wave function of two fermions is antisymmetric; thus, the final expression for the wave function takes the form where nðxÞ and pðxÞ are the neutron and proton field ¼ 1 ¼ 1 27 operators and g and g . are the vector and axial- JM JM jαþjβ−J JM V A Ψ ðr1; r2Þ¼N αβ½ψ ðr1; r2Þ − ð−Þ ψ ðr1; r2Þ ; vector coupling constants, respectively. An effective theory αβ αβ βα could also include Δ isobars, meson fields, and their vertices ð3:16Þ for decaying into lepton pairs and interacting with nucleons pffiffiffi and each other. where N αβ equals 1= 2 for α ≠ β and 1=2 for α ¼ β. As a result of the capture of electrons, the nucleus ðA; ZÞ As a consequence of the identity CJM ¼ð−1ÞJ−2jCJM , 0þ → π jm1jm2 jm2jm1 undergoes a J transition. Conservation of total angular the wave function of two electrons with equal quantum momentum requires that the captured electron pair be in the numbers α ¼ β is symmetric under the permutation mα ↔ state J. In weak interactions, parity is not conserved; thus, it is mβ provided that their angular momenta are combined to the not required that the parity of the electron pair be correlated total angular momentum J ¼ 2j mod(2). In such a case, the with the parity of the daughter nucleus. antisymmetrization (3.16) yields zero, which means that The wave function of a relativistic electron in a central ¼ 2 potential has the form the states J j mod(2) are nonexistent. The antisymmetri- zation (3.16) of the states J ¼ 2j þ 1 mod(2) leads to a doubling of the initial wave function. To keep the norm, the fαðrÞΩα ðnÞ pffiffiffi Ψ ðrÞ¼ mα ð Þ αmα ; 3:13 additional factor 1= 2 is thus required for α ¼ β. ð ÞΩ 0 ðnÞ igα r α mα The derivation of the equation for Vαβ is analogous to the 0 0 0 corresponding derivation of the 0ν2β decay amplitude, as where α ¼ðn2jlÞ, α ¼ðn2jl Þ, and l ¼ 2j − l. The described by Bilenky and Petcov (1987). The specificity is radial wave functions are defined in agreement with Ω ðnÞ ≡ that a transition from a discrete level to a quasidiscrete level is Berestetsky, Lifshitz, and Pitaevsky (1982); αmα l l0 considered. Accordingly, the delta function expressing the Ω ðnÞ, Ωα0 ðnÞ¼Ω ðnÞ are spherical spinors in the jmα mα jmα energy conservation is replaced by a time interval that can be notation of Varshalovich, Moskalev, and Khersonskii (1988). identified with the parameter t in Eq. (3.9). We thus write Ψ ðrÞ The normalization condition for αmα is given by Z hfjUðtÞjii¼−iVαβt þ ; ð3:17Þ rΨ† ðrÞΨ ðrÞ¼δ δ ð Þ d αmα βmβ αβ mαmβ : 3:14 where UðtÞ is Dyson’s U matrix. The amplitude takes the form
Z pffiffiffi 2 X −iq⋅ðr1−r2Þ GF cos θC dq 1 e Vαβ ¼ iK mββ 2N αβ pffiffiffi dr1dr2 pffiffiffiffiffiffiffiffiffiffiffiffiffiffi Z 2 ð2πÞ3 2 þ 1 2 J M q0 μν μν ½ JM ðr r Þ ðr r Þ − ð−1Þjαþjβ−J JM ðr r Þ ðr r Þ ð Þ × Tμναβ 1; 2 NJMαβ 1; 2 Tμνβα 1; 2 NJMβα 1; 2 ; 3:18 where X JM ðr r Þ¼ JM ½Ψ¯ c ðr Þγ ð1 þ γ Þγ Ψ ðr Þ Tμναβ 1; 2 Cjαmαjβmβ αmα 1 μ 5 ν βmβ 2 ; mαmβ X h j μðr Þj ih j ν ðr Þj00i h j ν ðr Þj ih j μðr Þj00i μν JM Jh 1 n n Jh 2 JM Jh 2 n n Jh 1 N ðr1; r2Þ¼ þ : JMαβ þ − − ε þ − − ε n q0 En Mi β q0 En Mi α