Arxiv:Hep-Ph/0312040V2 21 May 2004 Nacmn Ftecpuertsi Xetdt Cu Ten- at Resonance Occur 2 to a the Expected to Comparable Favoured

Home , Double electron capture, Electron capture

arXiv:hep-ph/0312040v2 21 May 2004 nacmn ftecpuertsi xetdt cu ten- at resonance occur 2 to A the expected to comparable favoured. is release rates strongly ergy capture are the atoms of enhancement Z High culated. nry h xeiet erhn o h 0 the decay for total searching the to experiments equal The is electrons energy. two the of energies of related ns asn h ea 1,[2]. [1], decay mech- the the the of of causing regardless nonconservation anism true the remains This as number. well lepton Majorana as the neutrino unambiguously of proves nature effect the of vation inlt h aoaanurn as hl h ex- successful the a While from deduced propor- mass mass. 0 are neutrino the neutrino process of value Majorana a act such the for to Amplitudes tional fig.1). (see lyietclwt hto h eetr n h tracking the and usu- detector, is the source of the that with of calorimetric identical material the ally the categories: which two in into experiments, divided be can cay easre naohroe hslast h nuclear the to leads This may one. 0 vertex the another reaction Thus interaction in weak absorbed conjugate. one charge be def- in day its by produced then to present particle identical neutrino the Majorana is a it of is inition topics neutrino are particles If intriguing these physics. most of nature the Majorana among or Dirac the and ASnmes 34.s 46.q 23.40.Bw 14.60.Pq, 23.40.-s, numbers: PACS an estimated is a feasibility to studies encouraging. experimental decays such found for for The Candidates met nuclei. be 4considered. final to in likely states are excited conditions resonance The etdlne ieie.Rtsfrte0 the ex- for the Rates offset lifetimes. partly longer advantages pected These suppres to background. trigger random coincidence of existence the portantly, aiew rps osuyteivrepoes h radiativ the 0 process, capture inverse major electron the a double study al is neutrino-less as to attractive mass an propose neutrino its As we of of native physics. measure neutrino nature sensitive to-day’s Majorana of a challenge provide the to prove as to well thus and cay aue te datgsicuetefvual ai fth of ratio favourable the 0 include advantages Other s nature. experimental convenient a provides photon monoenergetic ν νββ ∗ e-mails: 2 h 0 The Introduction. h osblt oosretenurn-esdouble neutrino-less the observe to possibility The EC etiols obeEeto atr olt erhfor search to tool a - Capture Electron Double Neutrino-less xeieti oe eedn,temr obser- mere the dependent, model is experiment otecmeig2 competing the to β νββ lcrn.Isuiu intr sta h sum the that is signature unique Its electrons. νββ ukwfweup,[email protected] [email protected], rcs rcesvateeiso ftocor- two of emission the via proceeds process ( ,Z A, : − h xsec fmsieneutrinos massive of existence The 2) ola nttt o ula tde,H˙a6,PL-00-681 69, Ho˙za Studies, Nuclear for Institute So ltan → ν 2 ( EC ,Z A, P − atr ae n,vr im- very and, rates capture + ) 1 S tmclvldifference. level atomic ν e γ 2 2 + EC EC e h associated The . rcs r cal- are process .ukwk n .Wycech S. and Z.Sujkowski νβ − β − β the s ter- de- (1) de- ig- re d e e 1 utfll h ihrslto eurmn.Otherwise requirement. resolution 0 experiment between high the of sandwiched the kinds sheets both fulfill in thin must detectors of The form detectors. the the order in the of tons, material, nec- source of the the handling of amount in large difficulties other cessarily practical the are On there background. hand, are the the to coincidence, tracking sensitive of in The less electrons material somewhat two the surrounding. the all the counting of of detectors, purity and conditions extreme housing shielding the detector extreme to as way requiring well by only as in is The rests difficulty this background. main attain quantities The random Large used. the be suppressing can emitted. ener- material particles the the charged are up of the sums detector of automatically the gies former and The source the separate. which in experiments, ttscntu ecniee nraitceprmnsin experiments realistic in the to considered contrast be thus can states suppress to • background; trigger random coincidence overwhelming precious the a provides this • [8]; low is process • • degra- energy without material dation; source the of layers thick : • process capture electron diative ruddet h 2 back- the physical to continuous due ground dominating the from cernible h rbesivle erfrt eiw 2 [6]. - β [2] reviews to refer we involved problems the apn hntepoo nryeul h nryof energy the equals energy which 2 effect photon atomic resonant the the when is radia- interest the happens case. special of capture of theory electron point a double The of the to basis adapted the capture on tive section encourag- next discussed are are the Crude These process in experiments. radiative photon. feasible a a suggesting such ing, by away for carried [7] is estimates energy excess The ffc a togyehnete0 the enhance strongly may effect fteiooe r icse ntefia sections. final the in discussed are isotopes choice optimal the the of and questions experimental The rates. trigger. coincidence characteristic another yet offer oiatapiue o h aitv obeelectron double radiative the for amplitudes dominant − h htneiso sfloe yta fteKX-ray; K the of that by followed is emission photon the h orecnb eaaefo h detector; the from separate be can source the h hsclbcgon u otecmeig2 competing the to due background physical the o ea nrisaefvue.Dcy oexcited to Decays favoured. are energies decay low fairly from easily escapes photon monoenergetic the hr r eea xeietlavnae ftera- the of advantages experimental several are There esgett td h nes fnurn-esdouble neutrino-less of inverse the study to suggest We mltdsfrterdaiecapture. radiative the for Amplitudes transition the of reduction sizeable a is pay to price The ea,ie h etiols obeeeto capture. electron double neutrino-less the i.e. decay, νβ − P β − − ∗ eki h u pcrmwl o edis- be not will spectrum sum the in peak 1 S ββ asw Poland Warsaw, , rniin nfruaestain this situations fortunate In transition. eas The decays. νβ − β − aoaaneutrinos Majorana ea.Frtedsrpinof description the For decay. γ rniin hc follow which transitions νcc rates. νeeγ The g (b) Two electrons are captured in a virtual process p n n p which generates a final atom with two electron holes. + W- W - + e e bound This final atom radiates and one of the holes is filled. m Process (a) is standard in the single electron radia- m n + n n n + n n L m e R L m e R tive captures. The description of electron propagation in- b- b- (Z,N) (Z+2,N – 2) (Z,N) ECEC (Z – 2,N+2) volves the electron Coulomb Green’s functions [12], [13].

m m The two electron capture requires some modifications. n + n n n L m e R n L + m n R e To conserve the angular momentum, for one electron in + e bound 1S state the other one must be in a higher nS state. The - W- e W+ transition is of magnetic type. The electron wave func- n p p n tions must be antisymmetrized. Process (b) is less important in general. However, for FIG. 1. Diagrams for the 0ν double beta decay and double small photon energy Q the virtual two electron hole state electron capture processes. may be degenerate to the final atom-photon state (see below, eq. (12); see also refs. [15], [16], [17]). This leads to a singularity in Rγ . The most interesting situation capture can be factorised into weak-nuclear, weak- 0ν occurs for Q = E(2P ) E(1S) i.e. when the final photon leptonic and radiative factors [9], [11]. The expressions resonates with the 2P − 1S transitions in the final atom. for radiative factors presented below are based on the The capture rate is enhanced− up to the limit given by the theory for the radiative single electron capture process natural K-L line widths. [12], [13]. The cases (a) and (b) require different descriptions. The amplitude for 0νββ decays from the (atom + In the simplest theory of Majorana neutrino [9], [11] the nucleus) ground state 0) to a final nuclear state Nf cc | matrix element R0ν is: and final two electron state nf is 2 ββ cc G ρ σ R0ν = (0 Hw nf ,Nf ). (2) R =2 dxdy JN (x) mν hν (x y)JN (y) | | 0ν √2 Z h − i This differs from the amplitude for the double electron C Ψ(n1, x)γρ(1 γ5)γσΨ¯ (n2, y). (5) capture × − The nuclear weak currents are denoted by J and the Rcc = (¯n ,N H 0) (3) N 0ν f f | w | effect of neutrino propagation is included into the ”neu- by the electron wave functions. Heren ¯ denotes two trino potential” with hν (r) exp( qr)/r, where q is f ≈ − electron holes in the final atom (say 1S, 2S or 1S, 2P ). an average momentum carried by the neutrino (the clo- The reversed transition cannnot be tested in the same sure approximation over nuclear states is used here). The cc nucleus and R0ν does not correspond to a real physical projection on left-handed intermediate neutrino brings transition. Assuming the radiative process instead, we about the neutrino mass factor mν . The electronic part can write the corresponding amplitudes in the second of this formula contains atomic wave functions Ψ and order perturbation theory, to the leading order in the the charge conjugates ΨC. For the 0+ 0+ transi- → radiative interaction Hγ, as tions the nuclear part is reduced to Fermi and Gamow- Teller matrix elements defined by MF = 0 hν 0 and (ni,Ni Hw 0) h | | i γ MGT = 0 hν σ1σ2 0 . These enter via a combination R0 =Σni,Ni (nf ,Nf Hγ ni,Ni) | | . ν 0ν h | gV |2 i | | E E(ni) E(Ni) M = MGT ( ) MF and in this way the transition − − − gA (4) matrix element (3) is brought to the form where H is the weak interaction describing the neutrino- 2 w cc G 0ν R =2 M mν [Ψ(n1, 0)S(1, 2)Ψ(n2, 0)]A, (6) less process of fig.1. Intermediate states ni,Ni) in the 0ν √2 radiative process may involve ground and| excited nuclear states. Only the electron is considered to radiate as the which resembles the standard ββ decay expression, [9]. γ emission by the nucleus is unlikely [10], [9]. The electron wave functions are required in the nuclear We look first for the transitions which conserve the region. These are Ψ(ne, 0)û, where Ψ(ne, 0) is the ra- nuclear angular momentum (0+ 0+) as these offer the dial part of the large component andu ˆ is the Dirac largest nuclear matrix elements.→ The final photon and spinor. The spin matrix element is given by S(1, 2) = C its spin is generated by the electron radiation. Two basic u(n1)(1 + γ5)¯u (n2). In eq.(6) the electron wave func- processes are possible: tion is to be antisymmetrised. For two 1S electrons this (a) One of the electrons captured from the initial involves antisymmetric spin zero combination of spinors atomic state nini radiates while it propagates towards which compensates an ”inner antisymmetry” built into the nucleus. operator S(1, 2) by the charge conjugation.

2 To describe the radiative capture of the type (a), one At smaller Q the capture amplitude indicates an inter- function Ψ is to be replaced by some function Ψγ which esting structure. The electron propagator in eq.(7) has a takes care of the photon emission and electron propaga- pole at Q = E(1S) E(2P ) Qres. To elucidate the tion. In eq.(4), in coordinate representation, this func- atomic physics| involved− in this|≡ radiation process, ampli- tion is given by the expression tude (4) is presented in terms of the initial state

(r ni)(ni Hγ nini) (ni,Ni Hw 0) Ψγ(r,nini)=Σni | | | (7) Ψini)= 0)+Σni,Ni | | ni,Ni), (12) E E(ni) E(0) | | E E(n ) E(N ) | − − − i − i which involves the sum over continuum and discrete where E = EN (0) + Eatom(0) is the initial energy of states of the electron, i.e. the Dirac Green’s function the system composed of the nuclear part EN and the in the external Coulomb field of the nucleus (Glauber atomic component Eatom. Equation (12) describes the and Martin [12], [13]). This solution is now implemented virtual transitions of the initial state 0) under weak into eq.(4) to give the radiative amplitude interactions which include virtual two| electron capture states. Some of these states may be degenerate with the 2 γ G 0ν ~ final atom + photon state to within the natural width of R0ν =2 mν M Ψ(n1, 0)S(1, 2,~ǫ, Q)γ Ψ(n2, 0)γ √2 A the K-L line. Rate enhancements of the order of 104 106 h i − (8) may be obtained. In the single electron description by Glauber and Martin the electron propagator in eq.(7) has a pole at Q = E(1S,Z) E(2P,Z) Q (Z) where Z where the leptonic spins and the photon polarisation ǫ | − |≡ res enter the last term refers to the initial atom. In the many electron process considered here we use the energy levels for final(Z 2) √α atom (see also ref. [17]). We obtain − S(1, 2,~ǫ, Q~ )γ = 2me r ~ ~ C 0νγ Γ (2P 1S) cc 2 (u(n1)[iA(Q)~σ~ǫ Q B(Q)~ǫQ](1 + γ5)¯u (n2)). (9) Γ (Q)= → R (13) × × − [Q Q (Z 2)]2 + [Γr/2]2 | 0ν | − res − For the discussion of A(Q) and B(Q) we refer to ref. where Γr is the radiative width of the final two electron [13]. If Q m then A, B are close to unity for both the e hole atom. We use the experimental energy values for the 1S and 2S≈ electrons. The spin conservation requires the Z 2 atoms corrected for the screening differences in the two electrons to be in the symmetric spin and antysym- ionised− atoms; the combined width is 3Γr(1S)+Γr(2P ) metric space combinations. This causes cancelations at similarly as for the hypersatellite KhX-ray spectra. small Z. Otherwise the effect of symmetrisation is small. Experimental feasibility, optimal targets. The rate for radiative no-neutrino process is now The usual choice for targets in the double β decay is

0νγ 2πdk γ 2 motivated by the phase space that favors large energy Γ (Q)=Σpol δ(k Q) R0ν , (10) 5 Z (2π)32k − | | release ( roughly Qββ dependence). In contrast, the phase space dependence of the double electron capture to be summed over photon polarisation and possible elec- eq.(11) is rather weak and the radiative factors favor tron pairs. The result is small Q. Two capture situations can be considered. At Q larger than the electron mass, the ”magnetic” radia- 4 0νγ G gV 2 2 2 tive capture of 1S, 2S pair dominates. The calculated Γ = (MGT ( ) MF ) (mν /me) √2 − gA rates (Table 1) are based on values MGT 0.6/fm and ≈ Q MF /MGT 0.3 obtained in typical ββ0ν reactions in Ψ(1S, 0)Ψ(2S, 0) 2 Σ S(n ,n ,~ǫ, Q~ ) 2 (11) ≈ − × | | 2π polh 1 2 γ i a number of nuclei [18], [19]. These values give a crude estimate also for the 0νcc processes, although some cal- where the summation over spin factor S is to be taken. h i culations indicate a larger coupling, [20]. This factor takes care of the angular momentum con- At small Q the capture amplitude is given by the res- servation. The best possibility is the 1S, 2S pair cap- onant conditions. Several candidates for such a resonant ture accompanied by a magnetic photon transition. This capture may be found. Examples are given in Table II. rate is scaled by the S electron wave function at the nu- The region of practical interest is limited to Q< 80 keV. 2 3 2 cleus. One has Ψ(n1S, 0) = (Zmeα) fs /π where factor Such conditions are likely to be met in transitions to ex- fs comes from the weak singularity of Dirac atomic wave cited states of final nuclei. The rate predictions are ham- function. Averaged over the nuclear volume it becomes pered by the poor knowledge of mass differences. Be- λ−1 2 fs = (2RZmeα) with λ = 1 (Zα) and the nu- cause of the sharp atomic resonances a few keV mass 6 − clear radius R. The Z factorp arises and thus one is uncertainty results in rates changing by several orders interested in the heaviest possible atoms. There the rel- (cf. Table II). The atomic g.s. mass differences ∆M are ativistic singularity enhances the rate further.

3 TABLE I. The radiative neutrino-less capture rates from TABLE II. The resonant situations. The radiative neu- the 1S, 2S states [1/y], mν = 1 eV. trino-less capture rates per year R/y, and per year and ton of the isotope R/y · ton, m = 1 eV. The uncertainties of R/y · Atom abundance % Q[keV ] Γ(2S, 1S) ν 92 −32 ton are due to 1σ errors in the mass determination. Q - pho- 42Mo 15.84 1628(5) 2 · 10 108 −31 ton energies, Qr - the resonant energies, QX ,Qγ - energies of 48 Cd 0.88 241(7) 1 · 10 ∗ 144 −31 accompanying X-ray, γ-ray photons, E - excitations of ﬁnal 62 Sm 3.07 1730(4) 1 · 10 162 −30 nuclei, [keV]. 68 Er 0.14 1781(4) 2 · 10 180 −31 112 136 152 162 164 180 74 W 0.12 66(5) 2 · 10 Atom 50 Sn 58 Ce 64 Gd 68 Er 68 Er 74 W 196 −30 80 Hg 0.15 729(3) 4 · 10 abnd.% 1.01 0.19 0.20 0.14 1.56 0.13 + + + f.st. 02 ,1S,2P 03 ,1S,2P g.s.,1S,2P 1 ,1S,2P g.s.,2S,2P g.s.,1S,2P + + + int.st. 02 ,1S,1S 03 ,1S,1S g.s.,1S,2S 1 ,1S,1S g.s.,2S,2S g.s.,1S,1S ∆M 1919.5(4.6) 2418.9(13) 54.6(3.5) 1843.8(3.9) 23.3(3.9) 144.4(4.5) ∗ taken from refs. [23], [24] and the Q values are ∆M re- E 1870.9 2315.4 1745.5 duced by the two electron-hole excitation energies. The Q 18.0(4.6) 60(13) 39.7(3.5) 36.1(3.9) 6(4) 68.6(4.5) precision needed is well below 1 keV. Qr 23.8 33 39.7 46.6 1 56.3(4.5) The experimental signature for the double electron Q − Qr -5.8(4.6) 27(13) 0.04(3.5) -10.5(3.9) 5(4) 12.3(4.5) capture in the resonance conditions will be the KX- QX 23 32 6.5 46 6.5 55 KhX-ray coincidence and, in the case of decays to ex- Qγ 1253.4 1496.9 1665.1 617.6 818.5 80.6 cited states, the triple coincidence with the gamma rays −30 −30 −25 −28 −32 −28 h R/y 7 · 10 2 · 10 2 · 10 1 · 10 2 · 10 3 · 10 deexciting these states. The K line has a measurable en- − − − − − − − R/y · t 1 − 10 2 10 1 − 10 3 103 − 10 1 1 − 10 1 10 1 − 10 3 3 − .2 ergy shift. The experimental feasibility arguments have to include the decay rate and the cost estimates. Leav- ing the cost arguments aside and assuming 1 ton of the 2 source material with 100 % isotopic purity and the cor- ννγ 2 0νγ (2π) Γ = dLδ(Q Σki)RN Γ (kγ ) (15) respondingly larger amount of a high resolution detector Z − kγ

(e.g. a multilayered source-detector sandwich for X-ray k k k′ 2πd γ d ν d ν detection combined with the medium resolution larger γ- where the phase space element dL = 9 ′ and (2π) 8kγ kν kν 2ν 0ν ray detectors) it seems possible to design feasible exper- RN = (4πM )/(M mν ) is a dimensionless ratio of the iments for the 0ν double electron capture process. The nuclear two-neutrino and no-neutrino matrix elements. count rates expected depend strongly on the energy dif- Typical values M 2ν 1, M 0ν 1/fm follow from nu- ≈ ≈ ference with respect to the resonance value (see Table clear model calculations, [3], [11]. To obtain RS/B, ex- II). The 112Sn isotope indicated in ref. [17] as the best pression (15) is presented as an integral over the photon choice, now with the the recent mass determination seems ννγ Q energy distribution Γ = 0 W (kγ )dkγ and the back- less profitable. Much higher decay rates can be expected ground contribution becomesR for 152Gd and 164Er g.s. g.s. decays. However, there → Q are no convenient experimental signatures in these tran- ννγ Γ ND = W (kγ )dkγ . (16) sitions. Best signatures are offered in the cases of decays ZQ−D/2 to excited states (e.g. 112Sn, 136Ce or 162Er) where there are characteristic high energy γ-rays in addition to the Formula W (k ) = (Q k )3Γ0νγ (k )R2 /6(2π)4 follows γ − γ γ N kαX-rays and the resonant transitions. from eq.(15). This indicates a cubic cut-off at the end of The physical background. One advantage of the the spectrum due to the phase space. The ratio RS/B is 0νγ process is a convenient ratio of the signal to the phys- now given by the interplay of this cut-off and the energy 0νγ ical background generated by the dominant ννγ channel. dependence of Γ (kγ ). Here, we briefly estimate this ratio defined as First, we consider the magnetic type transition related to the 1S, 2S electron capture. For characteristic values 0νγ 0νγ ννγ Γ (Q) Q = 1 MeV and mν = 1 eV one obtains Γ /Γ = RS/B = ννγ (14) −5 Γ ND 5 10 . The photon energy resolution D = 3 KeV would yield· very convenient R =2 106. The case of resonant where N is the fraction of photons from the dominant S/B D electron captures is more involved· due to the rapid energy ννγ decay mode emitted into the region from the end dependence in Γ0νγ (k ) The result is plotted in Fig.2. of the spectrum Q down to Q D/2 and D is the pho- γ − The Q Qres situation is very favorable and RS/B is ton energy resolution. More details may be found in ref. ≤9 r about 10 . For Q>Qres +Γ the RS/B falls down as [8]. For easier comparison the two-neutrino radiative rate − [Q Q ] 5. At Q Q = 1 keV the ratio is still very is presented in terms of the no-neutrino radiative rate res res convenient,− R −104, but the conditions deteriorate Γ0νγ (k ), given by eq.(11) or eq.(13), as a function of S/B γ quickly with the increasing≈ energy separation. the photon energy The random background (RB) varies with the source material. Crudely, the effect/BR ratios are expected

4 [1] R.N. Mohapatra and P.B. Pal, Massive Neutrinos, World

Scientiﬁc (2001) 1010 [2] F.Boehm and P.Vogel, Physics of Massive Neutrinos, R (res) S/B 109 Cambridge Univ.Press (1987) [3] H.Ejiri, Physics Reports 338, 265 (2000).

108 [4] S. Elliott and P. Vogel, Annu. Rev. Nucl. Part. Sci.52,(2002)115 107 0νγ [5] H.V. Klapdor-Kleingrothaus, Particles and Nuclei, Let- Γ (res)

B 1 / ters, JINR, ,20(2001). H.V. Klapdor-Kleingrothaus et

R (magn) S S/B 6 586 R 10 al., Phys.Lett.B ,198 (2004) [6] Z.Sujkowski, Acta Phys.Pol.B 34,2207(2003). 105 -5 [7] Z.Sujkowski and S.Wycech, Acta Phys.Pol.B ~ Q - Qres 33,471(2002). 4 10 [8] S.Wycech and Z.Sujkowski, Acta Phys.Pol.B 35,1221(2004), nucl-th/0402103. 103 -1.0 -0.5 0.0 0.5 1.0 1.5 [9] M.Doi, T.Kotani and E.Takasugi, Prog.Th.Phys.Suppl. 83 89 Q - Q [k eV] , 1 (1985); M.Doi, T.Kotani Prog.Th.Phys. , 139 res (1993). FIG. 2. Schematic plot for the signal/physical background [10] J.D.Vergados, Nucl.Phys.B 218,109 (1983). 300 ratios: RS/B(res)- continuous line for the resonant transition [11] J.Suhonen and O.Civitarese Phys.Rep. ,123 (1998). 104 and RS/B (magn)- dashed line for the magnetic transition. [12] R.J. Glauber and P.C. Martin, Phys.Rev. ,158(1956). 0νγ 109 The rate dependence Γ (res) vs (Q − Qres) is shown in [13] P.C. Martin, R.J. Glauber, Phys.Rev. ,1307(1958). arbitrary units to indicate the resonance width. [14] T.Tomoda, Rep.Progr.Phys.54,53 (1991). [15] R.G. Winter, Phys.Rev. 100,142(1955). [16] H.M.Georgi, S.L.Glashow and S.Nussinov, Nucl.Phys.B to be improved in comparison to the calorimetric 193,297 (1981). β−β−experiments [5] by a factor between 10 and 1000, [17] J.Bernabeu, A. de Rujula and C.Jarlskog, Nucl.Phys.B depending on the X-ray and the γ-ray coincidence trig- 223,15 (1983). ger. [18] J.Suhonen Phys.Rev.C62(2000)042510(R) 64 We conclude that detecting the neutrino-less double [19] F.Simkovic et al., Phys.Rev.C (2001)035501 beta process of any kind, be it β−, β+, or EC, will mean [20] S.Stoica and H.V.Klapdor-Klaingrothaus, Eur. Phys. Journ.A 17,529 (1998). discovering the lepton number nonconservation and prov- [21] T. Fritioff et al., Eur. Phys. J. A15 (2002) 423 ing the Majorana nature of neutrinos. It will also pro- [22] NDM03 Proc. Conf., Nara, June 2003. vide a sensitive measure of the neutrino mass, thus giv- [23] V.I.Tretyak, Y.G.Zdesenko, At. Dat. Nucl. Dat. Tabl. 66 ing answers to some of the most urgent questions of to- (2002) 83 day physics. This search calls for a massive, diversified [24] G.Audi,A.H.Wapstra and C.Thibault, Nucl.Phys. experimental as well as theoretical effort. The radia- A729(2003)337. tive neutrino-less double electron capture, proposed in this work, is a valid alternative to double β± emission, with several experimental advantages. The predicted resonance rate enhancement makes this search plausible. Exact rate predictions require precise measurements of masses and likeways calculations of nuclear matrix elements. The former are feasible with modern techniques [21], [24]. The latter are strongly assumption dependent. The consensus [22] is that they can be trusted to about a factor 3. One way to minimize this uncertainty is to carry out measurements for several nuclear species, preferably on both sides of the mass parabolae. Thus the double electron capture search can be considered not only as an attractive alternative to the double β−. Assuming both kinds of experiment successful, it will also provide a much needed complementary piece of information.