Sterile neutrinos are possible dark matter candidates. We examine here possible detection mechanisms, assuming that the neutrino has a mass of about 50 keV and couples to the ordinary neutrino. Even though this neutrino is quite heavy, it is nonrelativistic with a maximum kinetic energy of 0.1 eV. Thus new experimental techniques are required for its detection. We estimate the expected event rate in the following cases: (i) measuring electron recoil in the case of materials with very low electron binding; (ii) low temperature crystal bolometers; (iii) spin induced atomic excitations at very low temperatures, leading to a characteristic photon spectrum; (iv) observation of resonances in antineutrino absorption by a nucleus undergoing electron capture; (v) neutrino induced electron events beyond the end point energy of beta decaying systems, for example, in the tritium decay studied by KATRIN.
1. Introduction
There exists evidence for existence of dark matter in almost all scales, from the dwarf galaxies, galaxies, and cluster of galaxies, with the most important ones being the observed rotational curves in the galactic halos; see, for example, the review [1]. Furthermore cosmological observations have provided plenty of additional evidence, especially the recent WMAP [2] and Planck [3] data.
In spite of this plethora of evidence, it is clearly essential to directly detect such matter in the laboratory in order to unravel its nature. At present there exist many such candidates, called Weakly Interacting Massive Particles (WIMPs). Some examples are the LSP (Lightest Supersymmetric Particle) [4–11], technibaryon [12, 13], mirror matter [14, 15], and Kaluza-Klein models with universal extra dimensions [16, 17]. Among other things these models predict an interaction of dark matter with ordinary matter via the exchange of a scalar particle, which leads to a spin independent interaction (SI) or vector boson interaction and therefore to a spin dependent (SD) nucleon cross section.
Since the WIMPs are expected to be extremely nonrelativistic, with average kinetic energy 〈T〉≈50keV(mWIMP/100GeV), they are not likely to excite the nucleus, even if they are quite massive, mWIMP>100GeV. Therefore they can be directly detected mainly via the recoiling of a nucleus, first proposed more than 30 years ago [18]. There exists a plethora of direct dark matter experiments with the task of detecting WIMP event rates for a variety of targets such as those employed in XENON10 [19], XENON100 [20], XMASS [21], ZEPLIN [22], PANDA-X [23], LUX [24], CDMS [25], CoGENT [26], EDELWEISS [27], DAMA [28, 29], KIMS [30], and PICASSO [31, 32]. These consider dark matter candidates in the multi-GeV region.
Recently, however, an important dark matter particle candidate of the Fermion variety in the mass range of 10–100 keV, obtained from galactic observables, has arisen [33–35]. This scenario produces basically the same behavior in the power spectrum (down to Mpc scales) with that of standard ΛCDM cosmologies, by providing the expected large-scale structure [36]. In addition, it is not too warm; that is, the masses involved are larger than m = 1–3 keV to be in conflict with the current Lyα forest constraints [37] and the number of Milky Way satellites [38], as in standard ΛWDM cosmologies. In fact an interesting viable candidate has been suggested, namely, a sterile neutrino in the mass region of 48–300 keV [33–35, 39–43], but most likely around 50 keV. For a recent review, involving a wider range of masses, see the white paper [44].
The existence of light sterile neutrinos had already been introduced to explain some experimental anomalies like those claimed in the short baseline LSND and MiniBooNE experiments [45–47], the reactor neutrino deficit [48], and the Gallium anomaly [49, 50], with possible interpretations discussed, for example, in [51, 52] as well as in [53, 54] for sterile neutrinos in the keV region. The existence of light neutrinos can be expected in an extended see-saw mechanism involving a suitable neutrino mass matrix containing a number of neutrino singlets not all of which being very heavy. In such models it is not difficult to generate more than one sterile neutrino, which can couple to the standard neutrinos [55]. As it has already been mentioned, however, the explanation of cosmological observations requires sterile neutrinos in the 50 keV region, which can be achieved in various models [33, 56].
In the present paper we will examine possible direct detection possibilities for the direct detection of these sterile neutrinos. Even though these neutrinos are quite heavy, their detection is not easy. Since like all dark matters candidates move in our galaxies with no relativistic velocities, with average value about 10-3 c, and with energies about 0.05 eV, not all of them can be deposited in the detectors. Therefore the standard detection techniques employed in the standard dark matter experiments like those mentioned above are not applicable in this case. Furthermore, the size of the mixing parameter of sterile neutrinos with ordinary neutrinos is crucial for detecting sterile neutrinos. Thus our results concerning the expected event rates will be given in terms of this parameter.
The paper is organized as follows. In Section 2 we study the option on neutrino-electron scattering. In Section 3 we consider the case of low temperature bolometers. In Section 4 the possibility of neutrino induced atomic excitations is explored. In Section 5 we will consider the antineutrino absorption on nuclei, which normally undergo electron capture, and finally in Section 6 the modification of the end point electron energy in beta decay, for example, in the KATRIN experiment [57], is discussed. In Section 7, we summarize our conclusions.
2. The Neutrino-Electron Scattering
The sterile neutrino as dark matter candidate can be treated in the framework of the usual dark matter searches for light WIMPs except that its mass is very small. Its velocity follows a Maxwell-Boltzmann (MB) distribution with a characteristic velocity about 10-3 c. Since the sterile neutrino couples to the ordinary electron neutrino it can be detected in neutrino-electron scattering experiments with the advantage that the neutrino-electron cross section is very well known. Both the neutrino and the electron can be treated as nonrelativistic particles. Furthermore we will assume that the electrons are free, since the WIMP energy is not adequate to ionize an the atom. Thus the differential cross section is given by (1)dσ=1υCν2gV2+gA2GF222d3pν′2π3d3pe2π32π4δpν-pν′-peδpν22mν-p′ν22mν-pe22me,where Cν2 is the square of the mixing of the sterile neutrino with the standard electron neutrino νe and GF=Gcosθc, where G=1.1664×10-5 GeV-2 denotes the Fermi weak coupling constant and θc≃13° is the Cabibbo angle [58]. The integration over the outgoing neutrino momentum is trivial due to the momentum δ function yielding (2)dσ=1υCν2gV2+gA2GF22212π2d3peδpeυξ-pe22μr,where ξ=p^e·p^ν,0≤ξ≤1, υ is the WIMP velocity, and μr is the WIMP-electron reduced mass, μr≈mν. The electron energy T is given by (3)T=pe22me=2mν2meυξ2⟹0≤T≤2mν2meυesc2,where υesc is the maximum WIMP velocity (escape velocity). Integrating (2) over the angles, using the δ function for the ξ integration we obtain (4)dσ=Cν21υgV2+gA2GF22212πpe2dpe1peυ⟹dσ=Cν21υ2gV2+gA2GF22212πmedT.
We are now in a position to fold the velocity distribution assuming it to be MB with respect to the galactic center: (5)fυ′=1πυ03e-υ′/υ02.In the local frame, assuming that the sun moves around the center of the galaxy with velocity υ0=220km/s, υ′=υ+υ0z^, we obtain (6)fly,ξ=1πυ03e-1+y2+2yξ,y=υυ0,where ξ is now the cosine of the angle between the WIMP velocity υ and the direction of the sun’s motion. Eventually we will need the flux so we multiply with the velocity υ before we integrate over the velocity. The limits of integration are between υmin and υesc. The velocity is given via (3); namely, (7)υ=2meT2mνξ⟹υmin=2meT2mνWe find it convenient to express the kinetic energy T in units of T0=2(mν2/me)υ02. Then (8)ymin=x,x=TT0.Thus (9)υdσdT=1υ01T0me16πCν2GF2gV2+gA2∫xyescdyy2πe-1+y2∫-11dξe-2yξ.These integrals can be done analytically to yield(10)υdσdT=1υ0me16πCν2GF2gV2+gA2gx,gx=12erf1-x+erfx+1+erfc1-yesc+erfcyesc+1-2,where erf is the error function and erfc(x) is its complement. The function g(x) characterizes the spectrum of the emitted electrons and is exhibited in Figure 1 and it is without any particular structure, which is the case in most WIMP searches. For a 50 keV sterile neutrino we find that(11)T0=2mνme22.23210-6mec2≈5.0×10-3 eVTmax=T0yesc2=5×10-32.842≈0.04 eVT=1.6T0=8.0×10-3 eV.
The shape of the spectrum of the emitted electrons in sterile neutrino-electron scattering.
Now dT=T0dx. Thus (12)υσυ0=1υ02meT016πCν2GF2gV2+gA2∫0yesc2dxgx=1.43mν28πCν2GF2gV2+gA2,where (13)∫0yesc2dxgx=1.43.
It is clear that with this amount of energy transferred to the electron it is not possible to eject an electron out of the atom. One therefore must use special materials such that the electrons are loosely bound. It has recently been suggested that it is possible to detect even very light WIMPS, much lighter than the electron, utilizing Fermi-degenerate materials like superconductors [59]. In this case the energy required is essentially the gap energy of about 1.5kTc, which is in the meV region; that is, the electrons are essentially free. In what follows, we assume the values(14)gA=1,gV=1+4sin2θW=1.92,GF2=5.02×10-44 cm2/MeV2while Cν2 is taken as a parameter and will be discussed in Section 7. Thus we obtain(15)υσυ0=3.47×10-47Cν2 cm2.The neutrino particle density is (16)Nν=ρmν=0.3 GeV/cm350×10-6 GeV=6×103 cm-3while the neutrino flux (17)Φν=ρmνυ0=1.32×1011 cm-2s-1,where ρ=0.3 GeV/cm3, being the dark matter density. Assuming that the number of electron pairs in the target is 2×NA=2×1024 we find that the number of events per year is(18)Φνυσυ02×NA=2.89×10-4Cν2 y-1.
The authors of [59] are perhaps aware of the fact that the average energy for very light WIMPS is small and as we have seen above a small portion of it is transferred to their system. With their bolometer detector these authors probably have a way to circumvent the fact that a small amount of energy will be deposited, about 0.4 eV in a year for NA≈1024. Perhaps they may manage to accumulate a larger number of loosely bound electrons in their target.
3. Sterile Neutrino Detection via Low Temperature Bolometers
Another possibility is to use bolometers, like the CUORE detector exploiting Low Temperature Specific Heat of Crystalline TeO1302 at low temperatures. The energy of the WIMP will now be deposited in the crystal, after its interaction with the nuclei via Z-exchange. In this case the Fermi component of interaction with neutrons is coherent, while that of the protons is negligible. Thus the matrix element becomes (19)ME=GF22NgV,N=number of neutrons in the nucleus.
A detailed analysis of the frequencies of TeO1302 can be found [60]. The analysis involved crystalline phases of tellurium dioxide: paratellurite α-TeO2, tellurite β-TeO2, and the new phase -TeO2, recently identified experimentally. Calculated Raman and IR spectra are in good agreement with available experimental data. The vibrational spectra of α and β-TeO2 can be interpreted in terms of vibrations of TeO2 molecular units. The α-TeO2 modes are associated with the symmetry D4 or 422, which has 5 irreducible representations, two 1-dimensional representations of the antisymmetric type indicated by A1 and A2, two 1-dimensional representations of the symmetric types B1 and B2, and one 2-dimensional representation, usually indicated by E. They all have been tabulated in [60]. Those that can be excited must be below the Debye frequency which has been determined [61] and found to be quite low: (20)TD=232±7K°⟹ωD=0.024 eV.This frequency is smaller than the maximum sterile neutrino energy estimated to be Tmax=0.11 eV. Those frequency modes of interest to us are given in Table 1. The differential cross section is, therefore, given by(21)dσ=1υCν2N2gV2GF222∑k=18∑ni=0Nkd3pν′2π3d3q2π32π4F2q2δpν-pν′-qδpν22mν-p′ν22mν-ni+12ωk,where Ni will be specified below and q is the momentum transferred to the nucleus. The momentum transfer is small and the form factor F2(q2) can be neglected.
The frequency modes below the Debye temperature for α-TeO2 obtained from Table VIII of [60] (for notation see text).
νi=ωi2πcm-1
52
124
128
152
157
176
177
179
Symmetry
B1
E
B1
A1
B2
A2
E
B1
ωi(eV)
0.006
0.015
0.016
0.019
0.019
0.022
0.022
0.022
Ni
16
6
6
5
5
4
4
4
Emax(i)(meV)
106
100
102
103
107
98
99
100
In deriving this formula we tacitly assumed a coherent interaction between the WIMP and several nuclei, thus creating a collective excitation of the crystal, that is, a phonon or few phonons. This of course is a good approximation provided that the energy transferred is small, of a few tens of meV. We see from Table 1 that the maximum allowed energy is small, around 100 meV. We find that, if we restrict the maximum allowed energy by a factor of 2, the obtained results are reduced only by a factor of about 10%. We may thus assume that this approximation is good.
Integrating over the nuclear momentum we get (22)dσ=1υCν2N2gV2GF22212π2∑k=18∑ni=0Nkd3pν′δpν22mν-p′ν22mν-ni+12ωkdσ=1υCν2N2gV2GF22212π2∑k=18∑ni=0Nkd3pν′δpν22mν-p′ν22mν-ni+12ωk;performing the integration using the δ function we get (23)σ=1υCν2N2gV2GF2221πmν2mν∑k=18∑ni=1NkEν-ni+12ωkσ=υ0υCν2N2gV2GF2221πmν2∑k=18∑ni=0Nky2-ni+1/2ωkT1,where T1=(1/2)mνυ02,y=υ/υ0
Folding with the velocity distribution we obtain (24)υσ=υ0Cν2N2gV2GF2221πmν2∑k=18∑ni=0NkIni,ωk,Ini,ωk=∫yminyescdyfni,ωky,fni,ωky=y2-ni+1/2ωkT1ye-1-y2sinh2y,ymin=ni+1/2ωkT1.We see that we have the constraint imposed by the available energy; namely,(25)Nk=IPyesc2T1ωk-12,where IP[x]= integer part of x. We thus find Nk listed in Table 1. The functions fni,ωk(y) are exhibited in Figure 2. The relevant integrals are In(ω1)=(1.170,0.972,0.785,0.621) for n=0,1,2,3, In(ω2)=(1.032,0.609), for n=0,1, In(ω3)=(1.025,0.592), for n=0,1, and I0(ωk)=(0.979,0.970,0.934,0.932,0.929) for k=4,…,8. Thus we obtain a total of 17.8. The event rate takes with a target of mass mt which takes the form (26)R=ΦνCν2N2gV2GF218πmtAmpmν217.8.If we restrict the maximum allowed energy to half of that shown in Table 1 by a factor of two, we obtain 15.7 instead of 17.8.
(a) The function fn,ω1(y), exhibited as a function of y, associated with the mode ν1=52cm-1 for n=0,1,2,3 increasing downwards. (b) The functions fn,ω2(y) associated with ν2=124cm-1 for n=0,1 and f0,ω4 for ν4=157cm-1, n=0, and f0,ω6 for ν6=176cm-1, n=0, exhibited as a function of y, for thick solid, solid, dashed, and dotted lines, respectively. For definitions see text.
For a TeO1302 target (N=78) of 1 kg of mass get (27)R=1.7×10-6Cν2 per kg-s=51Cν2 per kg-yThis is much larger than that obtained in the previous section, mainly due to the neutron coherence arising from the Z-interaction with the target (the number of scattering centers is approximately the same, 4.5×1024). In the present case, however, targets can be larger than 1 kg. Next we are going to examine other mechanisms, which promise a better signature.
4. Sterile Neutrino Detection via Atomic Excitations
We are going to examine the interesting possibility of excitation of an atom from a level |j1,m1〉 to a nearby level |j2,m2〉 at energy Δ=E2-E1, which has the same orbital structure. The excitation energy has to be quite low; that is, (28)Δ≤12mνυesc2=1250×1032.8422.23210-6=0.11 eV.The target is selected so that the two levels |j1,m1〉 and |j2,m2〉 are closer than 0.11 eV. This can result from the splitting of an atomic level by the magnetic field so that they can be connected by the spin operator. The lower one |j1,m1〉 is occupied by electrons but the higher one |j2,m2〉 is completely empty at sufficiently low temperature. It can be populated only by exciting an electron to it from the lower one by the oncoming sterile neutrino. The presence of such an excitation is monitored by a tuned laser which excites such an electron from |j2,m2〉 to a higher state |j3,m3〉, which cannot be reached in any other way, by observing its subsequent decay by emitting photons.
Since this is a one-body transition the relevant matrix element takes the form (29)MEj1,m1;j2,m22=gV2δj1,j2δm1,m2+gA2Cl,j1,m1,j2,m22(in the case of the axial current we have gA=1 and we need evaluate the matrix element of σν·σe and then square it and sum and average over the neutrino polarization). (30)Cl,j1,m1,j2,m2=nlj2m2σnlj1m1=j1m1,1m2-m1∣j2m22j1+132l+16l12j1l12j2011expressed in terms of the Clebsch-Gordan coefficient and the nine-j symbol. It is clear that in the energy transfer of interest only the axial current can contribute to excitation.
The cross section takes the form (31)dσ=1υCν2GF222MEj1,m1;j2,m22d3pν′2π3d3pA2π32π4δpν-pν′-pAδEν-Δ-Eν′.Integrating over the atom recoil momentum, which has negligible effect on the energy, and over the direction of the final neutrino momentum and energy via the δ function we obtain (32)σ=1υCν2GF222MEj1,m1;j2,m221πEν-Δ2Eν-Δ-mνmν=1υCν2GF222MEj1,m1;j2,m221πmν22T1mνfy,ΔT1fy,ΔT1=y2-ΔT11/2,T1=12mνυ02,where we have set E-Δ=mν+T1-Δ≈mν.
Folding the cross section with the velocity distribution from a minimum Δ/T1 to yesc we obtain (33)υσυ0=Cν2GF222mν2MEj1,m1;j2,m221πgΔT1gΔT1=2π∫Δ/T1yescdyy2y2-ΔT11/2e-1+y2sinh2yy.Clearly the maximum excitation energy that can be reached is Δmax=2.842T0=0.108 eV. The function gΔ/T1 is exhibited in Figure 3.
The function gΔ/T1 for sterile neutrino scattering by an atom as a function of the excitation energy in eV.
Proceeding as in Section 2 and noting that for small excitation energy gΔ/T1≈1.4 we find(34)R=1.8×10-2Cν21ACl,j1,m1,j2,m22 kg-y.The expected rate will be smaller after the angular momentum factor Cl,j1,m1,j2,m2 is included (see Appendix A). Anyway, leaving aside this factor, which can only be determined after a specific set of levels is selected, we see that the obtained rate is comparable to that expected from electron recoil (see (18)). In fact for a target with A=100 we obtain (35)R=1.8×10-4Cν2Cl,j1,m1,j2,m22 kg-y.This rate, however, decreases as the excitation energy increases (see Figure 3). In the present case, however, we have two advantages.
The characteristic signature of photons spectrum is following the deexcitation of the level |j3,m3〉 mentioned above. The photon energy can be changed if the target is put in a magnetic field by a judicious choice of |j3,m3〉.
The target now can be much larger, since one can employ a solid at very low temperatures. The ions of the crystal still exhibit atomic structure. The electronic states probably will not carry all the important quantum numbers as their corresponding neutral atoms. One may have to consider exotic atoms (see Appendix B) or targets which contain appropriate impurity atoms in a host crystal, for example, chromium in sapphire.
In spite of this it seems very hard to detect such a process, since the expected counting rate is very low.
5. Sterile Neutrino Capture by a Nucleus Undergoing Electron Capture
This is essentially the process: (36)ν¯+eb+AN,Z⟶AN+1,Z-1∗involving the absorption of a neutrino with the simultaneous capture of a bound electron. It has already been studied [62] in connection with the detection of the standard relic neutrinos. It involves modern technological innovations like the Penning Trap Mass Spectrometry (PT-MS) and the Microcalorimetry (MC). The former should provide an answer to the question of accurately measuring the nuclear binding energies and how strong the resonance enhancement is expected, whereas the latter should analyze the bolometric spectrum in the tail of the peak corresponding to L-capture to the excited state in order to observe the relic antineutrino events. They also examined the suitability of 157Tb for relic antineutrino detection via the resonant enhancement to be considered by the PT-MS and MC teams. In the present case the experimental constraints are expected to be less stringent since the sterile neutrino is much heavier.
Let us measure all energies from the ground state of the final nucleus and assume that Δ is the mass difference of the two neutral atoms. Let us consider a transition to the final state with energy Ex. The cross section for a neutrino (here as well as in the following we may write neutrino, but it is understood that we mean antineutrino) of given velocity υ and kinetic energy Eν is given by (37)σEν=Cν21υMEExnuc2ϕe2GF222d3pA2π32π4δpA-pνδEν+mν+Δ-Ex-b,where pA is the recoiling nucleus momentum. Integrating over the recoil momentum using the δ function we obtain (38)σEν=Cν22π1υMEExnuc2ϕe2GF222δEν+mν+Δ-Ex-b.We note that since the oncoming neutrino has a mass, the excited state must be higher than the highest excited state at Ex′=Δ-b. With indicating by ϵ=Ex-Ex′ the above equation can be written as (39)σEν=Cν22π1υMEExnuc2ϕe2GF222δEν+mν-ϵ.Folding it with the velocity distribution as above we obtain (40)υσEν=Cν22πMEExnuc2ϕe2GF222∫0yescdyy22πe-1+y2sinh2yyδmν+12mνυ02y2-ϵor using the delta function(41)υσEν=2πCν2MEExnuc2ϕe2mνυ02GF222FX,FX=2πe-1+X2sinh2X,X=1υ02ϵmν-1.As expected the cross section exhibits resonance behavior though the normalized function F(X) as shown in Figure 4. It is, of course, more practical to exhibit the function F(X) as a function of the energy ϵ. This is exhibited in Figure 5. From this figure we see that the cross section resonance is quite narrow. We find that the maximum occurs at ϵ=mν(1+2.8×10-7)= 50 keV + 0.014 eV and has a width Γ=mν(1+9.1×10-7)-mν(1+0.32×10-7)≈ 0.04 eV. So for all practical purposes it is a line centered at the neutrino mass. The width may be of some relevance in the special case whereby the excited state can be determined by atomic deexcitations at the sub-eV level, but it will not show up in the nuclear deexcitations.
The cross section exhibits resonance behavior. Shown is the resonance properly normalized as a function of X=(1/υ0)2ϵ/mν-1. The width is Γ=1.49 and the location of the maximum is at 1.03.
The cross section exhibits resonance behavior. Shown is F(X) as a function of (ϵ/mν-1).
If there is a resonance in the final nucleus at the energy Ex=ϵ+(Δ-b) with a width Γ then perhaps it can be reached even if ϵ is a bit larger than mν; for example, ϵ=mν+Γ/2. The population of this resonance can be determined by measuring the energy of the deexcitation γ-ray, which should exceed by ϵ the maximum observed in ordinary electron capture.
For antineutrinos having zero kinetic energy the atom in the final state has to have an excess energy Δ-(b-mν) and this can only happen if this energy can be radiated out via photon or phonon emission. The photon emission takes place either as atomic electron or nuclear level transitions. In the first case photon energies are falling in the eV-keV energy region and this implies that only nuclei with a very small Δ-value could be suitable for this detection. In the second case, there should exist a nuclear level that matches the energy difference Ex=Δ-(b-mν) and therefore the incoming antineutrino has no energy threshold. Moreover, spontaneous electron capture decay is energetically forbidden, since this is allowed for Ex<Δ-(b+mν).
As an example we consider the capture of a very low energy ν¯ by the 65157Tb nucleus: (42)ν¯+e-+Tb65157⟶Gd64165∗taking the allowed transitions from the ground state (3/2+) of parent nucleus, Tb65157, to the first excited 5/2+ state of the daughter nucleus Gd64157. The spin and parity of the nuclei involved obey the relations ΔJ=1, ΠfΠi=+1, and the transition is dubbed as allowed. The nuclear matrix element ME can be written as (43)ME2=gAgV2GT2,where gA=1.2695 and gV=1 are the axial and vector coupling constants, respectively. The nuclear matrix element is calculated using the microscopic quasi-particle-phonon (MQPM) model [63, 64] and it is found to be |ME|2=0.96. The experimental value of first excited 5/2+ is at 64 keV [65] while that predicted by the model is at 65 keV. The Δ-value is ranging from 60 to 63 keV [65].
For K-shell electron capture where 〈ϕe〉2=(αZ/π)me3 (1s capture) with binding energy b=50.24 keV, the velocity averaged cross section takes the value (44)συ=8.98×10-46Cν2 cm2and the event rate we expect for mass mt=1 kg is (45)R=8.98×10-46Cν2×6103×6.0231023×mtA×9.28×1017y-1=19Cν2y-1.
The life time of the source should be suitable for the experiment to be performed. If it is too short, the time available will not be adequate for the execution of the experiment. If it is too long, the number of counts during the data taking will be too small. Then one will face formidable backgrounds and/or large experimental uncertainties.
The source should be cheaply available in large quantities. Clearly a compromise has to be made in the selection of the source. One can be optimistic that such adequate quantities can be produced in Russian reactors. The nuclide parameters relevant to our work can be found in [66] (see also [67]), summarized in Table 2.
Nuclides relevant for the search of the keV sterile neutrinos in the electron capture process. We give the life time T1/2, the Q-value, the electron binding energy Bi for various captures, and the value of Δ=Q-Bi. For details see [66].
Nuclide
T1/2
EC transition
Q (keV)
Bi (keV)
Bj (keV)
Q-Bi (keV)
Tb157
71 y
3/2+→3/2-
60.04(30)
K: 50.2391(5)
LI: 8.3756(5)
9.76
Ho163
4570 y
7/2-→5/2-
2.555(16)
MI: 2.0468(5)
NI: 0.4163(5)
0.51
Ta179
1.82 y
7/2+→9/2+
105.6(4)
K: 65.3508(6)
LI: 11.2707(4)
40.2
Pt193
50 y
1/2-→3/2+
56.63(30)
LI: 13.4185(3)
MI: 3.1737(17)
43.2
Pb202
52 ky
0+→2-
46(14)
LI: 15.3467(4)
MI: 3.7041(4)
30.7
Pb205
13 My
5/2-→1/2+
50.6(5)
LI: 15.3467(4)
MI: 3.7041(4)
35.3
Np235
396 d
5/2+→7/2-
124.2(9)
K: 115.6061(16)
LI: 21.7574(3)
8.6
6. Modification of the End Point Spectra of β Decaying Nuclei
The end point spectra of β decaying nuclei can be modified by the reaction involving sterile (anti)neutrinos: (46)ν+AN,Z⟶AN-1,Z+1+e-or (47)ν¯+AN,Z⟶AN+1,Z-1+e+.This can be exploited in on ongoing experiments, for example, in the tritium decay: (48)ν+H13⟶He23+e-.The relevant cross section is (49)σEν=Cν21υMEExnuc2GF222d3pA2π3d3pe2π32π4δpν-pA-peδEν+Δ-Ee,where Δ is the atomic mass difference. Integrating over the nuclear recoil momentum and the direction of the electron momentum we get (50)σEν=Cν21υMEExnuc2GF2221πEePe,where (51)Ee=mν+12mνυ2+Δ+me,Pe=Ee2-me2.
Folding the cross section with the velocity distribution we find (52)συ=Cν2GF2222π3/2∫0yescdyfy,where (53)fy=ME2ysinh2yEePee-1+y2FZf,Eewith (54)y=υυ0.The Fermi function, F(Zf,Ee), encapsulates the effects of the Coulomb interaction for a given lepton energy Ee and final state proton number Zf. The function f(y) is exhibited in Figure 6.
The shape of f(y) for the decay of 3H, where the atomic mass difference between3H and 3He is taking Δ=18.591keV [69].
In transitions happening inside the same isospin multiplet (Jπ→Jπ,J≠0) both the vector and axial form factors contribute and in this case the nuclear matrix element ME(Ex) can be written as (55)ME2=F2+gAgV2GT2,where gA=1.2695 and gV=1 are the axial and vector coupling constants, respectively. In case of H3 target we adopt 〈F〉2=0.9987 and 〈GT〉2=2.788 from [68]. Thus |ME|2=5.49.
Thus the velocity averaged cross section takes the value (56)συ=3.44×10-46Cν2 cm2and the expected event rate becomes (57)R=3.44×10-46Cν2×6103×6.0231023×mtA×9.28×1017y-1.For a mass of the current KATRIN target, that is, about 1 gr, we get (58)R=0.380Cν2y-1.
It is interesting to compare the neutrino capture rate(59)Rν=συρmν=3.44×10-46Cν2×6103×9.28×1017=1.91×10-24Cν2y-1with that of beta decay process 3H→3He+e-+ν¯j, whose rate Rβ is given by(60)Rβ=GF22π3∫meWopeEeFZ,EeME2EνpνdEe,where Wo is the maximal electron energy or else beta decay endpoint(61)Wo=Kend+mewith(62)Kend=m3H-me2-m3He+mν22m3H≃Δ=18.591 keV,the electron kinetic energy at the endpoint, and (63)me≈510.99891013 keVmH3≈2808920.820523 keVmHe3≈2808391.219324 keV.Masses mH3 and mHe3 are nuclear masses [58, 69, 70]. The calculation of (60) gives Rβ=0.055y-1. The ratio of Rν to corresponding beta decay Rβ is very small.(64)Rν=0.034·10-21Cν2Rβ.The situation is more optimistic in a narrow interval Wo-δ<Ee<Wo near the endpoint. As an example, we consider an energy resolution δ=0.2 eV close to the expected sensitivity of the KATRIN experiment [57]. Then the ratio of the event rate Rβ(δ=0.2 eV) to that of neutrino capture Rν gives (65)Rν=5.75·10-9Cν2Rβδ=0.2 eV.
In Figure 7 we present the ratio of the event rate decay rate of Rβ(δ) for the beta decay compared with the neutrino capture rate Rν as a function of the energy resolution δ in the energy region Wo-δ<Ee<Wo.
Ratio of decay rates Rν/Rβ (in units of Cν2) as a function of energy resolution δ near the endpoint.
Moreover, the electron kinetic energy Ke due to neutrino capture process (48) is(66)Ke=Eν+Kend>mν+18.591 keV;this means that the electron in the final state has a kinetic energy of at least mν above the corresponding beta decay endpoint energy. There is no reaction induced background there, but, unfortunately, the ratio obtained above is much lower than the expected KATRIN sensitivity.
7. Discussion
In the present paper we examined the possibility of direct detection of sterile neutrinos of a mass 50 keV, in dark matter searches. This depends on finding solutions to two problems. The first is the amount of energy expected to be deposited in the detector and the second one is the expected event rate. In connection with the energy we have seen that, even though these neutrinos are quite heavy, their detection is not easy, since like all dark matters candidates move in our galaxies with not relativistic velocities, 10-3 c on the average, and with energies about 0.05 eV, not all of which can be deposited in the detectors. Thus the detection techniques employed in the standard dark matter experiments, like those looking for heavy WIMP candidates, are not applicable in this case.
We started our investigation by considering neutrino-electron scattering. Since the energy of the sterile neutrino is very small one may have to consider systems with very small electron binding, for example, electron pairs in superconductors, which are limited to rather small number of electron pairs. Alternatively one may use low temperature bolometers, which can be larger in size resulting in a higher expected event rate. These experiments must be able to detect very small amount of energy.
Then we examined more exotic options by exploiting atomic and nuclear physics. In atomic physics we examined the possibility of spin induced excitations. Again to avoid background problems the detector has a crystal operating at low temperatures. Then what matters is the atomic structure of the ions of the crystal or of suitably implanted impurities. The rate in this case is less than that obtained in the case of bolometers, but one may be able to exploit the characteristic feature of the spectrum of the emitted photons.
From the nuclear physics point of view, we consider the antineutrino absorption on an electron capturing nuclear system leading to a fine resonance in the (N+1,Z-1) system, centered 50 keV above the highest excited state reached by the ordinary electron capture. The deexcitation of this resonance will lead to a very characteristic γ ray. Finally the sterile neutrino will lead to ν+A(N,Z)→e-+A(N-1,Z+1) reaction. The produced electrons will have a maximum energy which goes beyond the end point energy of the corresponding β decay essentially by the neutrino mass. The signature is less profound than in the case of antineutrino absorption.
Regarding the event rate, as we have mentioned before, it is proportional to the coupling of the sterile neutrino to the usual electron neutrino indicated above as Cν2. This parameter is not known. In neutrino oscillation experiments a value of Cν2≈10-2 has been employed. With such a value our results show that the 50 keV neutrino is detectable in the experiments discussed above. This large value of Cν2 is not consistent, however, with a sterile 50 kev neutrino. In fact such a neutrino would have a life time [71] of 2×105 y, much shorter than the age of the universe. A cosmologically viable sterile 50 keV neutrino is allowed to couple to the electron neutrino with coupling Cν2<1.3×10-7. Our calculations indicate that such a neutrino is not directly detectable with experiments considered in this work. The results, however, obtained for the various physical processes considered in this work, can be very useful in the analysis of the possible experimental searches of lighter sterile neutrinos in the mass range of 1–10 keV.
The angular momentum coefficients entering single particle transitions are shown in (A.1) and (B.4).
Equation (A.1). The coefficients (Cl,j1,m1,j2,m2)2 connect via the spin operator a given initial state |i〉=nl,j1,m1 with all possible states |f〉=|nl,j2,m2〉, for l=0,1. Note s-states are favored. (A.1)lj1m1j2m2Cl,j1,m1,j2,m22012-1212122,iflj1m1j2m2Cl,j1,m1,j2,m22112-12121229112-1232-3243112-1232-1289112-123212491121232-12491121232128911212323243132-3232-1223132-1232128913212323223.
B. Exotic Atomic Experiments
As we have mentioned the atomic experiment has to be done at low temperatures. It may be difficult to find materials exhibiting atomic structure at low temperatures. It amusing to note that one may be able to employ at low temperatures some exotic materials used in quantum technologies (for a recent review see [72]) like nitrogen-vacancy (NV), that is, materials characterized by spin S=1, which in a magnetic field allow transitions between m=0,m=1 and m=-1. These states are spin symmetric. Antisymmetry requires the space part to be antisymmetric, that is, a wave function of the form (B.1)ψ=ϕnl2rL=odd, S=1J=L-1,L,L+1.Of special interest are (B.2)ψ=ϕnl2r3PJ,ϕnl2r3FJ.Then the spin matrix element takes the form (B.3)L3J2m2σL3J1m1=12J2+1J1m1,1m2-m1∣J2m2L3J2σL3J1,L=P,F.The reduced matrix elements are given in (B.5), as well as the full matrix element 〈PJ2m23σPJ1m13〉2 of the most important component.
Equation (B.4). The same as in equation (A.1), the coefficients Cl,j1,m1,j2,m22 for l=2 are(B.4)iflj1m1j2m2Cl,j1,m1,j2,m22232-3232-12625232-3252-5285232-3252-351625232-3252-12425232-123212825232-1252-322425212-1252-122425232-12521212252321232326252321252-12122523212521224252321252322425232325212425232325232162523232525285252-5252-3225252-3252-121625252-1252121825252125232162525232525225.
Equation (B.5). The coefficients are 〈3PJ2σ3PJ1〉, 〈FJ23σFJ13〉, and 〈PJ2m23σPJ1m13〉2. For the notation see text.(B.5)J1J2PJ23σPJ130123111212562252J1J2FJ23σFJ1322-103232533376343244152J1m1J2m2PJ2m23σ3PJ1m12001m2291-11-1161-1101610100101116111116
Disclosure
Permanent address of John D. Vergados is as follows: University of Ioannina, 451 10 Ioannina, Greece.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors are indebted to Professor Marco Bernasconi for useful suggestions in connection with the phonon excitations of low temperature bolometers.
UllioP.KamioknowskiM.Velocity distributions and annual-modulation signatures of weakly-interacting massive particles200120010310.1088/1126-6708/2001/03/049SpergelD.BeanR.DoréO.Three-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Implications for Cosmology20071702AdeP. A. R.ArnaudM.AshdownM.Planck 2013 results. XXXI. Consistency of the Planck data201457110.1051/0004-6361/201423743BottinoA.de AlfaroV.FornengoN.MoralesA.PuimedonJ.ScopelS.Exploring the supersymmetric parameter space by direct search for WIMPs19974021-211312110.1016/S0370-2693(97)00390-0ArnowittR.NathP.Event rates in dark matter detectors for neutralinos including constraints from b-->s gamma decay199574234592459510.1103/PhysRevLett.74.4592ArnowittR.NathP.Predictions of neutralino dark matter event rates in minimal supergravity unification19965432374238410.1103/PhysRevD.54.2374BednyakovV.Klapdor-KleingrothausH.KovalenkoS.Possible constraints on SUSY-model parameters from direct dark matter search199432915910.1016/0370-2693(94)90509-6EllisJ.RoszkowskiL.Supergravity dark matter199228325226010.1016/0370-2693(92)90016-WEllisJ.FloresR. A.Elastic supersymmetric relic-nucleus scattering revisited199126325926610.1016/0370-2693(91)90597-JEllisJ.FloresR. A.Prospects for neutralino detection with a 73Ge + 76Ge detector19933001-217518210.1016/0370-2693(93)90767-CEllisJ.FloresR. A.Implications of LEP on laboratory searches for dark matter neutralinos1993400253610.1016/0550-3213(93)90396-7NussinovS.Some estimates of interaction in matter of neutral technibaryons made of colored constituents19922791-211111610.1016/0370-2693(92)91849-5GudnasonS. B.KouvarisC.SanninoF.Dark matter from new technicolor theories20067409500810.1103/PhysRevD.74.095008FootR.LewH.VolkasR. R.A model with fundamental improper spacetime symmetries19912721-2677010.1016/0370-2693(91)91013-LFootR.Mirror & hidden sector dark matter in the light of new CoGeNT data20117031713ServantG.TaitT. M. P.Is the lightest Kaluza–Klein particle a viable dark matter candidate?20036501-239141910.1016/S0550-3213(02)01012-XOikonomouV.VergadosJ.MoustakidisC. C.Direct detection of dark matter rates for various wimps20077731-21942GoodmanM. W.WittenE.Detectability of certain dark-matter candidates198531123059306310.1103/PhysRevD.31.3059AngleJ.AprileE.ArneodoF.Search for light dark matter in XENON10 data2013107051301AprileE.ArisakaK.ArneodoF.Dark matter results from 100 live days of XENON100 data2011107131302AbeaK.HosakaaJ.IidaT.Distillation of liquid xenon to remove krypton2009314290296GhagC.AkimovD. Y.AraújoH. M.Performance of the veto detector incorporated into the ZEPLIN-III experiment20113527686NiK.2011Beijing, China10.1142/9789814449373_0030MallingD. C.After LUX: The LZ Programhttps://arxiv.org/abs/1110.0103D. S. Akeribetal.Limits on spin-dependent WIMP-nucleon interactions from the Cryogenic Dark Matter Search200673, 011102AalsethC. E.BarbeauP. S.BowdenN. S.Cabrera-PalmerB.ColaresiJ.CollarJ. I.DazeleyS.de LurgioP.DrakeG.FastJ. E.FieldsN.GreenbergC. H.HossbachT. W.KeillorM. E.KephartJ. D.MarinoM. G.MileyH. S.MillerM. L.OrrellJ. L.RadfordD. C.ReynaD.RobertsonR. G. H.TalagaR. L.TenchO.Van WechelT. D.WilkersonJ. F.YocumK. M.Results from a search for light-mass dark matter with a p-type point contact germanium detector201110613130110.1103/PhysRevLett.106.131301ArmengaudE.AugierC.BenoîtA.Final results of the EDELWEISS-II WIMP search using a 4-kg array of cryogenic germanium detectors with interleaved electrodes20117025329335BernabeiR.BelliP.CappellaF.CerulliR.DaiC. J.d’AngeloA.HeH. L.IncicchittiA.KuangH. H.MaJ. M.MontecchiaF.NozzoliF.ProsperiD.ShengX. D.YeZ. P.First results from DAMA/LIBRA and the combined results with DAMA/NaI200856333335510.1140/epjc/s10052-008-0662-yBelliP.BernabeiR.BottinoA.Observations of annual modulation in direct detection of relic particles and light neutralinos201184, 055014LeeH. S.Limits on Interactions between Weakly Interacting Massive Particles and Nucleons Obtained with CsI(Tl) Crystal Detectors200799, 091301ArchambaultS.AubinF.AugerM.Dark matter spin-dependent limits for WIMP interactions on 19F by PICASSO20096822185192ArchambaultS.AubinF.AugerM.New Insights into Particle Detection with Superheated Liquids201113043006ArgüellesC.MavromatosN.RuedaJ.RuffiniR.The role of self-interacting right-handed neutrinos in galactic structure201620160403803810.1088/1475-7516/2016/04/038RuffiniR.ArguellesC. R.RuedaJ. A.On the core-halo distribution of dark matter in galaxies20154511622ArguellesR.RuffiniR.RuedaJ. A.Novel constraints on fermionic dark matter from galactic observableshttps://arxiv.org/abs/1606.07040BoyarskyA.RuchayskiyO.ShaposhnikovM.The role of sterile neutrinos in cosmology and astrophysics200959135937910.1146/annurev.nucl.010909.083654BoyarskyA.LesgourguesJ.RuchayskiyO.VielM.Realistic sterile neutrino dark matter with keV mass does not contradict cosmological bounds2009102, 201304TollerudE. J.BullockJ. S.StrigariL. E.WillmanB.Hundreds of milky way satellites? Luminosity bias in the satellite luminosity function200868812772892-s2.0-5584908506810.1086/592102CamposM. D.RodejohannW.Testing keV sterile neutrino dark matter in future direct detection experiments201694910.1103/PhysRevD.94.095010LiY.zhongX. Z.ShuL.Direct detection of the cosmic neutrino background including light sterile neutrinos10106924LiY. F.zhong XingZ.Captures of hot and warm sterile antineutrino dark matter on EC-decaying Ho-163 nuclei20110800610.1088/1475-7516/2011/08/006GiorgioA.The Waning of the WIMP? A Review of Models, Searches, and Constraintshttps://arxiv.org/abs/1703.07364LiaoW.201082073001AdhikariR.AgostiniM.Anh KyN.A White Paper on keV sterile neutrino Dark Matter20172017Aguilar-ArevaloA.Evidence for neutrino oscillations from the observation of ¯νe appearance in a ¯νμ beam200164, 112007Aguilar-ArevaloA. A.BrownB. C.BugelL.Improved Search for ¯νμ→¯νe Oscillations in the MiniBooNE Experiment20131106KatoriT.ConradJ. M.Beyond Standard Model Searches in the MiniBooNE Experiment201520151936297110.1155/2015/362971MuellerT. A.LhuillierD.FallotM.Improved predictions of reactor antineutrino spectra201183505461510.1103/PhysRevC.83.054615AbdurashitovJ. N.Measurement of the response of a Ga solar neutrino experiment to neutrinos from a 37Ar source200673, 045805KoppJ.MachadoP. A. N.MaltoniM.SchwetzT.Sterile neutrino oscillations: the global picture20131305GiuntiC.LavederM.LiY. F.LongH. W.Pragmatic view of short-baseline neutrino oscillations201388710.1103/PhysRevD.88.073008MentionG.FechnerM.LasserreT. H.Reactor antineutrino anomaly201183, 073006AsakaT.BlanchetS.ShaposhnikovM.The νMSM, dark matter and neutrino masses19856314151156KusenkoA.Sterile neutrinos: The dark side of the light fermions20094811-2128HeurtierL.TeresiD.Dark matter and observable lepton flavor violation2016941210.1103/PhysRevD.94.125022MavromatosN. E.PilaftsisA.Anomalous Majorana Neutrino Masses from Torsionful Quantum Gravity201286124038KATRIN collaborationKATRIN: A next generation tritium beta decay experiment with sub-eV sensitivity for the electron neutrino masshttps://arxiv.org/abs/hep-ex/0109033BeringerJ.Review of Particle Physics201286, article 010001010001HochbergY.PyleM.ZhaoY.ZurekK. M.Detecting superlight dark matter with Fermi-degenerate materials201620168, article no. 572-s2.0-8498127626710.1007/JHEP08(2016)057CeriottiM.PietrucciF.BernasconiM.2006731010.1103/PhysRevB.73.104304BarucciM.BrofferioC.GiulianiA.GottardiE.PeroniI.VenturaG.Measurement of Low Temperature Specific Heat of Crystalline TeO2 for the Optimization of Bolometric Detectors20011235-63033142-s2.0-034786893010.1023/A:1017555615150VergadosJ.NovikovY. N.Prospects of detection of relic antineutrinos by resonant absorption in electron capturing nuclei201441, article 125001125001ToivanenJ.SuhonenJ.Microscopic quasiparticle-phonon description of odd-A Xe isotopes19952111ToivanenJ.SuhonenJ.Microscopic quasiparticle-phonon description of odd-mass 127−133Xe isotopes and their β decay19985731237124510.1103/physrevc.57.1237HelmerR. G.Nuclear data sheets for A = 157200410345657822-s2.0-1294424932010.1016/j.nds.2004.11.007FilianinP. E.BlaumK.EliseevS. A.GastaldoL.NovikovY. N.ShabaevV. M.TupitsynI. I.VergadosJ.On the keV sterile neutrino search in electron capture201441910.1088/0954-3899/41/9/0950040950042-s2.0-84907202500VergadosJ. D.NovikovY. N.Exploring new features of neutrino oscillations with very low energy monoenergetic neutrinos20108391-2120SchiavillaR.StoksV. G. J.GlöckleW.Weak capture of protons by protons1998581263AudiG.WapstraA.ThibaultC.The Nubase evaluation of nuclear and decay properties200272913128BeardenJ. A.BurrA. F.Reevaluation of X-ray atomic energy levels196739112514210.1103/revmodphys.39.125DolgovA. D.HansenS. H.Massive sterile neutrinos as warm Dark Matter2002163339344DohertyM. W.MansonN. B.DelaneyP.JelezkoF.WrachtrupJ.HollenbergL. C. L.The nitrogen-vacancy colour centre in diamond20135281145