Abstract

In the present work we examine the possibility of detecting electrons in light dark matter searches. These detectors are considered to be the most appropriate for detecting dark matter particles with a mass in the MeV region. We analyze theoretically some key issues involved in such detection. More specifically we consider a particle model involving WIMPs interacting with fermions via -exchange. We find that for WIMPs with mass in the electron mass range the cross section for WIMP-atomic electron scattering is affected by the electron binding. For WIMPs more than 20 times heavier than the electron, the binding affects the kinematics very little. As a result, many electrons can be ejected with energy which increases linearly with the WIMP mass, but the cross section is somewhat reduced depending on the bound state wave function employed. On the other hand for lighter WIMPs, the effect of binding is dramatic. More specifically at most 10 electrons, namely, those with binding energy below 10 eV, become available even in the case of WIMPs with a mass as large as 20 times the electron mass. Even fewer electrons contribute if the WIMPs are lighter. The cross section is, however, substantially enhanced by the Fermi function corrections, which become more important at low energies of the outgoing electrons. Thus events of 0.5–2.5 per kg-y become possible.

1. Introduction

The combined MAXIMA-1 [13], BOOMERANG [4, 5] DASI [6], and COBE/DMR Cosmic Microwave Background (CMB) observations [7] imply that the universe is flat [8] and that most of the matter in the universe is dark [9], that is, exotic. These results have been confirmed and improved by the recent WMAP [10] and Planck [11] data. Combining the data of these quite precise measurements one finds (the more recent Planck data yield a slightly different combination, , . It is worth mentioning that both the WMAP and the Plank observations yield essentially the same value of , but they differ in the value of , namely, (WMAP) and (Planck). Since any “invisible” nonexotic component cannot possibly exceed of the above [12], exotic (nonbaryonic) matter is required and there is room for cold dark matter candidates or WIMPs (Weakly Interacting Massive Particles).

Even though there exists firm indirect evidence for a halo of dark matter in galaxies from the observed rotational curves, see, for example, the review [13], it is essential to directly detect such matter in order to unravel the nature of the constituents of dark matter. The possibility of such detection, however, depends on the nature of the dark matter constituents and their interactions.

The WIMPs are expected to have a velocity distribution with an average velocity which is close to the rotational velocity of the Sun around the galaxy; that is, they are completely nonrelativistic. In fact a Maxwell-Boltzmann leads to a maximum energy transfer which is close to the average WIMP kinetic energy . Thus for GeV WIMPs this average is in the keV regime, not high enough to excite the nucleus, but sufficient to measure the nuclear recoil energy. For light dark matter particles in the MeV region, which we will also call WIMPs, the average energy that can be transferred is in the few eV regions. So this light WIMPs can be detected by measuring the electron recoil after the collision. Electrons may of course be produced by heavy WIMPs after they collide with a heavy target which results in a shakeup of the atom yielding “primordial" electron production [1416]. This approach for sufficiently heavy WIMPs and target nuclei can produce electrons energies even in the 30 keV region, with a spectrum very different from that arising after a direct WIMP-electron collision. Furthermore WIMP-electron collisions involving WIMPs with masses in the few GeV region have also recently appeared [17, 18]. In the present work, however, we will restrict ourselves in the case of light WIMPs with a mass in the region of the electron mass.

We will draw from the experience involving WIMPs in the GeV region. The event rate for such a process can be computed from the following ingredients [19]: (i) the elementary WIMP-electron cross section. The total Higgs decay width [20]  MeV determined by the LHC data imposes severe constraints on the WIMP-fermion interactions mediated with Higgs exchange. Thus we will consider the case of a WIMP which interacts with fermions via -exchange. In this case we will relate the needed WIMP-electron cross to the isovector spin dependent WIMP-nucleon cross and exploit the information available for the latter. (ii) The knowledge of the WIMP particle density in our vicinity: this is extracted from WIMP density in the neighborhood of the solar system, obtained from the rotation curves measurements. The number density of these MeV WIMPs, however, is expected to be six orders of magnitude bigger than that of the standard WIMPs due to the smaller WIMP mass involved. (iii) The WIMP velocity distribution: in the present work we will consider a Maxwell-Boltzmann (MB) distribution.

In the electron recoil experiments, like the nuclear measurements first proposed more than 30 years ago [21], one has to face the problem that the process of interest does not have a characteristic feature to distinguish it from the background. So since low counting rates are expected the background is a formidable problem. Some special features of the WIMP interaction can be exploited to reduce the background problems, such as the modulation effect: this yields a periodic signal due to the motion of the Earth around the Sun. Unfortunately this effect, also proposed a long time ago [22] and subsequently studied by many authors [19, 2331], is small in the case of nuclear recoils, but we expect it to be a bit larger in the case of the electron recoils. There has always been an interest in light WIMPs, see, for example, the recent work [32]. In fact the first direct detection limits on sub-GeV dark matter from XENON10 have recently been obtained [33]. This is encouraging, but based on our experience with standard nuclear recoil experiments to excited states [34], one has to make sure that the proper kinematics has to be used in dealing with bound electrons. Clearly the binding electron energy plays a similar role as the excitation energy of the nucleus, in determining the small fraction of the WIMP’s energy to be transferred to the recoiling system. It is therefore clear that light WIMPs are quite different in energy, mass, interacting particle, and flux. Accordingly one needs detectors capable of detecting low energy light WIMPs in the midst of formidable backgrounds, that is, detectors which are completely different from current WIMP detectors employed for heavy WIMP searches.

2. The Particle Model

The narrow decay width of the Higgs boson obtained with LHC data imposes severe constraints on the Higgs mediated WIMP-SM fermion interaction. It may thus be more favorable to look for - boson mediated interactions, which have so far been considered only in the hadron spin induced dark matter searches. One finds that the SM coupling of the to the electron takes the form , where is the Weiberg angle. On the other hand the corresponding isovector quark interaction, expected to dominate in the scattering, takes the form , which leads to at the nucleon level. It is therefore natural to relate, if possible, the elementary WIMP-electron cross section to that of the spin induced isovector elementary WIMP-nucleon interaction, since for the latter there exists more theoretical and experimental information

Indeed if both the WIMP and the target with mass are nonrelativistic particles the cross section can be cast in the form Integrating over the momenta we find where with . Thus

To this end for the elementary isovector proton cross section we writewhere is a constant. Before proceeding further let us examine available information on [35], which is summarized here: (a) supersymmetric models in which the LSP (lightest supersymmetric particle) communicates with the quarks via the exchange of the -boson [36]. In this case in the constrained minimal supersymmetric standard model one finds sizable spin induced cross sections [37, 38], and as large as  pb (the proton neutron representation was chosen, since the experiments are analyzed this way). If, however, the isoscalar contribution is negligible the proton and the neutron amplitudes are opposite and . (b) Kaluza-Klein theories in models with Universal Extra Dimensions [39] in which the dark matter is a heavy neutrino communicating with matter via -exchange. In such models the couplings are of the order of the standard model couplings and one finds  pb. (c) The WIMP is a spin 3/2 particle [40]. In this case only the isovector contribution exists, leading to  pb. (d) Experimental limits: there exist some, namely, for Xe and Xe [41] and F [4246]. From the Xe data a limit is extracted on the elementary neutron SD cross section of pb and  pb for the proton SD cross section, while from the F target a slightly smaller limit is extracted on the proton SD cross section,  pb. These limits were based on nuclear physics considerations, namely, the nuclear spin matrix elements in the proton neutron representation. This explains the difference of the two limits extracted from the Xe data. For illustration purposes we will adopt the value [35]  pb. With this value of we find that the parameter relevant for electron scattering is

3. Experimental Aspects

Light WIMPs are quite different in energy, mass, interacting particle, and flux. Accordingly one needs detectors which are completely different from current WIMP detectors for heavy WIMPs. Detectors are required to observe low energy light WIMP signals beyond/among BG signals to identify the light WIMPs. Experimental aspects to be considered for light WIMP detectors are: (i) the particle to be detected, (ii) the event rate, (iii) the signal pulse height, (iv) the background rate, and (v) the detector threshold energy. We will now examine each of these items.

(1) Particle to Be Detected. Light WIMPs are detected by observing a recoil/scattered electron in the continuum region. In case that the WIMP interaction produces an ion-electron pair, one can detect the ion and/or the electron, and/or photons associated with the ion-electron pair. If the recoil electron or the ion-electron pair energy is absorbed by the detector, one may measure the temperature change. These are similar to those from heavy WIMPs except that their energies are very different. It is noted that atomic bound electrons are not excited by the light WIMPs with a few eV.

(2) Event Rate. The cross section of  pb is an order of magnitude larger than the present XENON limit of  pb for heavy WIMPs [47]. The flux rate is around , which is larger by a factor of 50 GeV/0.5 MeV ≈ . Then the event rate for Xe detector is around /(t·y), which is also a bit larger than the present limit for 50 GeV heavy WIMPs [47].

(3) Signal Pulse Height. The electron signal energy for light WIMPs is around 0.5–1 eV. This energy is 4 orders of magnitude smaller than the Xe nuclear recoil energy of around 25 keV for the 50 GeV WIMP. The nuclear recoil signal is quenched by a factor 2–20, depending on the atomic number, in most heavy WIMP detectors. Thus the actual signal height for the light WIMP is 3 orders of magnitude smaller than that for the heavy WIMP.

(4) Background Rate. There are three types of background origins for WIMP detectors, radioactive (RI) impurities, neutrons associated with cosmic rays, and electric noises. - rays from RI impurities produce BG electron signals, which are similar to electron signals from light WIMPs as well as those encountered in the case of double decay detectors, which measure rays. BG rate for a typical future DBD (double decay) detectors is around 1/(t y keV) = /(t y eV) at a few MeV regions [48]. Then one may expect a similar BG rate in the eV region. This is 3 orders of magnitude smaller than the signal rate. Neutrons do not contribute to BGs in light WIMP detectors, although nuclear recoils from neutron nuclear reactions are most serious BGs for heavy WIMP detectors. Electric noises are most serious for light WIMP detectors because of the very low energy signals. The nuclear recoil energy from heavy WIMPs is typically a few 10 keV, and the signal pulse height is around a few keV if they are quenched, depending on the detector. This is of the same order of magnitude as electric noise levels. Thus one can search for heavy WIMPs by measuring the higher velocity component above the electric noises On the other hand, the signal height for light WIMP is far below that of typical electric noises for current heavy WIMP detector.

(5) Energy Threshold. The energy threshold for WIMP detectors is set necessarily below the WIMP signal, but just above the electric noise to be free from the noise. Then a very low energy threshold of an order of sub eV is required for light WIMP searches. This is 3-4 orders of magnitude smaller than the level around 1–3 keV for most heavy WIMP detectors [47, 49]. Germanium semiconductor detectors are widely used to study low energy neutrinos and WIMPs. The ionization energy is 0.67 eV. Thus it can be used in principle for energetic light WIMPs. In practice, their threshold of around 200 eV [50] or more is still far above the light WIMP signals. Bolometers are, in principle, low energy threshold and high energy resolution detectors, but the energy threshold of practical 10 kg-scale bolometers is orders of magnitude higher than the light WIMP signal. Thus light WIMP detectors are necessarily different types from the present heavy WIMP detectors.

It is indeed a challenge to develop light WIMP detectors with low-threshold energy of the order of eV. Since the event rate is as large as 2 /t y, one can use a small volume detector of the order of 10 kgr at low temperature. In general, electric noises are random in time. Then, coincidence measurements of two signals are quite effective in reducing electric noise signals in case that one light WIMP produces 2 or more signals. One possible detector would be an ionization- scintillation detector, where one light WIMP interaction produces one ionized ion and one electron. In case that the ionized ion traps an electron nearby and emits a scintillation photon, one may measure the primary electron in coincidence with the scintillation photon. Nuclear emulsion may be of potential interest for low energy electrons. We briefly discuss possible new detection methods in Section 7.

4. The Differential WIMP-Electron Rate

The evaluation of the rate proceeds as in the case of the standard WIMP-nucleon scattering, but we will give the essential ingredients here to establish notation. We will begin by examining the case of a free electron.

4.1. Free Electrons

The differential cross section for WIMP-electron scattering in an atom can be cast in a form containing the electron recoil energy:

Thus the differential cross section takes the formor

We find from (8) that the average energy of the electron is given byFrom (8) we also find that the fraction of the energy of the WIMP transferred to the electron, when taking for scattering at forward angles, isThis situation is exhibited in Figure 1. We thus see that this fraction attains a maximum when , that is, when the two masses are equal. Away from this value it becomes smaller. The effect is more crucial for very light WIMPs, since their average energy is much smaller. Thus for MeV WIMPs the average energy transfer is in the eV region, which is reminiscent of the standard WIMPs where GeV mass leads to an energy transfer in the keV region. The maximum energy transfer corresponds to the escape velocity which is , which leads to a value four times higher. The exact expression of the maximum electron energy will be given below.

From (8) we find thatIn other words the minimum velocity consistent with the energy transfer and the WIMP mass is constrained as above. The maximum velocity allowed is determined by the velocity distribution and it will be indicated by . From this we can obtain the differential rate per electron in a given velocity volume as follows:where is the velocity distribution of WIMPs in the laboratory frame. Integrating over the allowed velocity distributions we obtain is a crucial parameter. Before proceeding further we find it convenient to express the velocities in units of the Sun’s velocity. We should also take note of the fact the velocity distribution is given with respect to the center of the galaxy. For a M-B distribution this takes the formWe must transform it to the local coordinate system:with , a unit vector in the Sun’s direction of motion, a unit vector radially out of the galaxy in our position, and . The last term in parenthesis in (17) corresponds to the motion of the Earth around the Sun with km/s being the modulus of the Earth’s velocity around the Sun and the phase of the Earth ( around June 3rd). The above formula assumes that the motion of both the Sun around the Galaxy and of the Earth around the Sun are uniformly circular. Since is small, we expand the distribution in powers of and we keep terms up to linearity in . Then (15) can be cast in the formwhere, in the above equation, the first term in parenthesis represents the average flux of WIMPs, the second provides the scale of the elementary cross section (in the present model the elementary cross section contains an additional mass dependence), the third term gives the number of electrons available for the scattering in a target of mass containing atoms with mass number and active electrons , and the fourth is essentially the inverse of the square of the Sun’s velocity in units of (its origin has its root in (10)). Moreover, and are functionals of and , where In the above expression the Heaviside function guarantees that the required kinematical condition is satisfied. One can factor the constants out in the above equation to getwhereThe meaning of will become clear after we consider the fact that the electrons are not free but bound in the atom. Thus they are not all available for scattering; that is, . We will now estimate considering the following input: (a) the elementary cross section ; (b) the total cross section, in units of , for example, = 0.2 for a WIMP mass about the electron mass (see below); (c) the particle density of WIMPs in our vicinity: (we use the electron mass in this estimate; the correct mass dependence has been included in evaluating ); this value leads to a flux: ; (d) The number of electrons in a Kg of Xe: . Taking , that is, all electrons in Xe participating, we expect about ≈0.3 events per kg·y. Encouraged by this estimate, even though it has been obtained with a much smaller elementary cross section than previous estimates [32], we are going to proceed in evaluating the expected spectrum of the recoiling electrons.

The shape of the differential cross section as a function of the electron energy in eV for various WIMP masses is exhibited in Figure 2. The electrons are assumed to be free. Obviously, this does not exhibit any special structure.

The various atomic physics approximations previously considered involve relatively high electron energies, as, for example, in the recent works [17, 18], which are not important in our case. The obtained rate, however, can be modified substantially at low energies by including the correction of the outgoing electron wave due to the Coulomb field. In beta decay this is done via the simple Fermi function [51]:whereand is the Coulomb function represented by a confluent hypergeometric function.

After some approximations (see, e.g., Landau’s book [52]) one finds the Fermi functionThis function, which depends on the electron energy as well as the atomic parameters and , may lead to an enhancement for low energy electrons. Integrating the differential rate given by (20) over the electron spectrum we obtain the total rate:where and are the average and time dependent (modulated) cross sections, respectively, in units of , that is, of the elementary cross section. The obtained quantity is shown in Figure 3 for free electrons both without the Fermi function as well as with the Fermi function for two values and , assumed to be sort of average. We see that the effect on the total cross section is large. Thus we find that the cross sections (in units of ), for a WIMP with the mass of the electron we get values 0.2, 0.8, and 2.2 for cases (a), (b), and (c) respectively (see Figure 3), while for and we find 0.3. In our estimates we will adopt an average enhancement factor of 8 due to the Fermi function taking into account the Fermi function the above estimate becomes 2.5 events/(kg·y).

5. Effects of Binding of Electrons

The binding of the electrons comes in two ways. The first is the most obvious. A portion of the energy of the WIMP will not go to recoil, but it will be spent to release the bound electron. The second comes from the fact that the initial electron is not at rest but it has a momentum distribution, which is the Fourier transform of its wave function in coordinate space.

Since the propagator involves the - exchange with a large mass, the momentum transfer can be neglected. So the cross section for WIMP-electron scattering is proportional towhere is the electron wf in momentum space. A plane wave outgoing electron has been assumed. Deviations from plane wave will be accounted for in the usual way, that is, by the Fermi function as described in the previous section. After the integration over the momentum via the function we obtainThe integration over the magnitude of can be done using the energy conserving function and we obtainwith and a unit vector in the direction of . We thus find the important constraintThis already sets a limit on the range of the variables and of interest to experiments.

5.1. The Allowed Range of Binding Energy and Recoil Energy

The range of these quantities is not affected by the specific form of the bound electron wave function. In particular regarding the above conditions imply where is associated with zero maximum recoil energy and is exhibited as a function of WIMP mass in Figure 4. What seriously affects the detection of light WIMPs in the presence of large binding energies is the minimum WIMP velocity required to eject an electron. One finds that in order to surpass the barrier of a given binding energy the WIMP must have a minimum velocity at most and a high mass, even if the electron energy is zero. The minimum required for this purpose is exhibited in Figure 5(a). The actual value of must, of course, be larger to get a reasonable rate.

Finally we present in Figure 5(b) the maximum possible energy for outgoing electrons as a function of the WIMP mass for various binding energies. We see that the maximum energy of the electron increases due to the binding provided that the binding energy is small. Thus the threshold energy needs not be extremely low, provided that for there exist enough electrons in the atom with low enough binding energies. For heavier WIMPs the binding energy of the electron plays a minor role and the electron recoil energy increases linearly with the WIMP mass.

5.2. The Effect of the Bound Electron Wave Function

To proceed further we must evaluate the integral in (28). Firstly, we must select a suitable coordinate system, for example, one with the -axis along the initial WIMP velocity, the -axis in the direction of the outgoing WIMP, and -axis perpendicular to the plane of the other two. Actually, there is no hope for obtaining an analytic expression even for the angular integrals without further approximations. So we will assume hydrogenic wave functions. This approximation is expected to be reasonable for electrons which are not tightly bound, which are of interest to us in the present work. We concentrate on the case of states and find that the cross section is given bywhere is a complicated functional given bywhereand , the scale of momentum of the bound electron wave function. is the invariant amplitude for the process which has the formNow, sinceone can cast (35) in a form similar to the expression for free electrons, namely,or where

In order to get the differential rate one must fold the above expression with the velocity distribution in a similar fashion as in (15), that is,Note the appearance of the quantities and as a result of the electron binding. The behavior of the function as a function of the velocity and its numerical value significantly affects the obtained rates. In order to obtain exactly one must evaluate the remaining integrals numerically for each electron orbit separately. After that, in order to obtain the total rate, one must sum up over all available electrons involved in any orbit of a given atom, which is not trivial. This is currently under study, but here we will report results obtained in an approximate scheme valid for relatively low mass WIMPs, which is adequate for our purposes.

5.3. A Convenient Approximation for Light WIMPs

We have seen that the momentum of the outgoing electron is small compared to . Let us assume that This means that This quantity is small as long as For such values of we can expand the integral up to second order in the small parameter. The result, for example, for hydrogenic wave functions, is Integrating over the angles of the outgoing electron we obtainProceeding as above we obtainorwhere and are in eV and is the WIMP velocity in units of the Sun’s velocity. A similar expression with a slightly different constant is expected to hold for other electron orbitals. This function must be multiplied with the velocity distribution before proceeding with the needed integrations to obtain the rate.

6. Some Results for Bound Electrons

We will limit ourselves to eV and After integrating with the velocity distribution we obtain the electron spectra shown in Figures 68. After integrating over the energy spectrum we obtain the cross section in units of shown in Figure 9 as a function of the WIMP mass for various binding energies. It is perhaps better to show as a function of the binding energy. For only are available. For we find in units of . For larger , is exhibited in Figure 10.

We thus see that using hydrogenic wave functions we find that, for , greater , electrons with eV become available. The number of electrons with relatively small binding energy for some targets of interest is shown in Table 1. Anyway there seems to be a wide range of the parameters which leads to = 5 for a suitable atom with large . Thus for one can conservatively set    to be ≈0.15, which leads to about 0.6 event per kg-y compared to the 2.0 per kg-y we got above for lighter WIMPs without the electron binding. The Fermi function is incorporated into the results.

7. Discussion

We have seen that the use of electron detectors may be a good way for directly detecting light WIMPs in the MeV region. The electron density in our vicinity is very high, the elementary WIMP-electron cross section may be quite large, and the event rate may be further enhanced by the behavior of the Fermi function at low energies. This is fine for WIMPs heavier than 20 times the electron mass. For WIMPs less than about 20 times the mass of the electron, however, there appears to be a problem arising from the electron binding. A small fraction of electrons can be exploited, namely, those electrons with binding energies below the 10 eV range. This is reminiscent of the difficulty encountered in the inelastic WIMP nucleus scattering, whereby only very low excited states can be reached. We have seen, however, that, due to the behavior of the Fermi function in the case of very low energy electron recoils, the rates may increase up to factors of 8. So event rates of about 0.8 to 2.4 events per kg-y are possible.

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 [53]. In this case the energy required is essentially the gap energy of about , which is in the meV region; that is, the electrons are essentially free. The authors are perhaps aware of the fact that the average energy for very light WIMPs is small and not all of it can be transferred to their system (the maximum fraction occurs if the mass of the WIMP is equal to , i.e., for in Figure 1). With their detector these authors probably have a way to circumvent the fact that a small amount of energy will be deposited, if they manage to accumulate a large number of loosely bound electrons in their targets so that the obtained rates are sufficient. More recently it is claimed that even smaller energies in meV can be detected in the case of Liquid Helium [54]. The expected event rates and the total energy deposited in such essentially bolometer type detectors are currently being estimated more precisely and they will appear elsewhere.

It thus appears that light WIMPs in the MeV region can, in principle, be detected. The detection techniques and targets employed, however, may have to be different than the ones employed in standard WIMP searches.

Disclosure

The permanent address of J. D. Vergados is University of Ioannina, 451 10 Ioannina, Greece.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

J. D. Vergados is happy to acknowledge support of this work by (i) the National Experts Council of China via a “Foreign Master” grant and (ii) IBS-R017-D1-2016-a00 in the Republic of Korea. A substantial part of this work performed while J. D. Vergados was on a visit to the University of Adelaide, supported by CoEPP and the Centre for the Subatomic Structure of Matter (CSSM). He is happy to thank Professors Yannis K. Semertzidis of KAIST, Anthony Thomas of Adelaide, and Edna Cheung of Nanjing University for their hospitality. Yeuk-Kwan E. Cheung acknowledges support from the Jiangsu Ministry of Science and Technology under Contract BK20131264 and the Priority Academic Program Development for Jiangsu Higher Education Institutions (PAPD).