Detecting superlight dark matter with Fermi-degenerate materials

We examine in greater detail the recent proposal of using superconductors for detecting dark matter as light as the warm dark matter limit of O\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{O} $$\end{document}(keV). Detection of suc light dark matter is possible if the entire kinetic energy of the dark matter is extracted in the scattering, and if the experiment is sensitive to O\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{O} $$\end{document}(meV) energy depositions. This is the case for Fermi-degenerate materials in which the Fermi velocity exceeds the dark matter velocity dispersion in the Milky Way of ∼ 10−3. We focus on a concrete experimental proposal using a superconducting target with a transition edge sensor in order to detect the small energy deposits from the dark matter scatterings. Considering a wide variety of constraints, from dark matter self-interactions to the cosmic microwave background, we show that models consistent with cosmological/astrophysical and terrestrial constraints are observable with such detectors. A wider range of viable models with dark matter mass below an MeV is available if dark matter or mediator properties (such as couplings or masses) differ at BBN epoch or in stellar interiors from those in superconductors. We also show that metal targets pay a strong in-medium suppression for kinetically mixed mediators; this suppression is alleviated with insulating targets.


Introduction
The identity of dark matter (DM) remains one of the most important mysteries in particle physics. In order to unlock the underlying nature of the DM, we rely on theories to help design and guide experiments. The dominant theoretical paradigm of massive DM over the last three decades has been the weakly interacting massive particle (WIMP), and for good reason: the observed density of DM is then naturally obtained via a freeze-out process while simultaneously ameliorating the infamous Standard Model (SM) sensitivity to ultraviolet physics, known as the hierarchy problem.

JHEP08(2016)057
As a result of this focus on DM at the weak scale, the sensitivity of experiments to such DM has dramatically increased. Direct and indirect detection experiments have made impressive gains in improving constraints on DM interaction rates with the Standard Model.
The enormous progress in the field will allow these DM experiments to push through important benchmarks in the next five to ten years. Ton scale direct detection experiments such as LUX [1], LZ [2] and Xenon1T [3] already have sensitivity to Higgs-interacting neutralino DM. The next generations of multi-ton experiments will have sensitivity not only to tree level scattering but even to one loop processes, such as wino DM scattering off nucleons through a loop of gauge bosons [4]. Indirect detection experiments are already constraining loop-generated annihilation processes, such as thermal relic neutralino DM annihilation into photons [5][6][7]. The combination of direct and indirect detection experiments will probe much of the viable parameter space for the neutralino in the minimal supersymmetric Standard Model in the upcoming years.
At the same time, theoretical developments have emphasized that compelling models of DM may be found beyond the weak scale, especially where the dark sector is complex and displays new dynamics. Unlike in the standard picture of DM, where the dark matter is inert since the time that its density is set in the early universe, such theories give rise to astrophysical and cosmological signatures that evolve with the universe itself. They often feature dark sectors with multiple particles and new dark forces. The explosion of interest in these sectors has accompanied studies of light hidden sectors at the LHC, 'Hidden Valleys' [8], where weak scale states decay into a complex dark sector with complex new dynamics. Thus, even if supersymmetry is discovered at the LHC, it in no way decreases the motivation to look for new physics beyond the SM at a much lower scale.
Most importantly, because the masses and composition of particles in such sectors are different than in the standard WIMP paradigm, new experiments must be designed to search for these dark sectors. Examples are asymmetric DM from a hidden sector [9] (natural mass scale m X Ω X Ωp m p ≈ 5 GeV) and mirror DM [10,11], where the masses are just below the reach of current direct detection probes. Going to lower masses, wellmotivated theories generated radiatively from the weak scale naturally live in the MeV-GeV mass scale [12][13][14], as well as models where the relic density is set via strong interactions, such as SIMPs [15,16]. Experiments are already moving toward detecting these theories of light DM. To this end, the detection of smaller energy deposits in the DM interaction process is required. For example, in nuclear elastic scattering processes, the deposited energy is E D q 2 /(2m N ), where m N is the target mass and the momentum transfer q ∼ µ r v X is set by the DM-nucleus reduced mass µ r and the DM velocity v X ∼ 10 −3 . The deposited energy on, e.g., a germanium nucleus is approximately 10 eV for 1 GeV DM, while past direct detection experiments, focused on weak scale DM, were sensitive to energy deposits between 5 and 100 keV. In order to access lighter DM candidates such as asymmetric DM, SuperCDMS [17][18][19] has lowered its energy sensitivity to 300 eV and plans to go to lower energies still.
Greater sensitivity to lighter DM for a given deposited energy can be achieved via inelastic processes, such as electron ionization or excitation [20][21][22]. In this case, the process may be catalyzed if the incoming DM kinetic energy, E kin = m X v 2 X /2, exceeds the JHEP08(2016)057 binding energy. Since semi-conductors feature valence electrons with binding energies as small as a few eV, semi-conductor based experiments such as SuperCDMS may be used to detect MeV DM scattering with electrons. The primary challenge is then to make heat sensors with sufficiently good energy resolution to detect 1 eV deposits of energy on electrons. This is a development challenge that SuperCDMS is currently taking on, and a subject that we return to in section 2. To go to even lower scales, however, will require even lower thresholds. Since the threshold of a semi-conductor experiment is fundamentally limited by the ionization energy of valence electrons, a new type of technology will need to be developed. In this paper we further develop the proposal laid out in ref. [23], where superconductors were considered for accessing DM energy deposits on electrons as low as a milli-eV, translating to sensitivity to DM with mass as low as m X ∼ keV. Cosmologically, the keV mass scale is significant because it corresponds to the lower bound on the DM mass from the Lyman-α forest [24] and phase space packing [25,26]; lighter (fermionic and thermalized) DM is inconsistent with cosmological observations (though see ref. [27]).
There are three features which make superconductors good DM detectors. First, ordinary metals have vanishing ionization threshold for electrons, implying no gap and hence access (in principle) to arbitrarily low DM energy depositions. Second, metals are Fermidegenerate, meaning that the conduction electrons follow a Fermi-Dirac distribution, having a Fermi velocity which is typically quite substantial, v F ∼ 10 −2 . As we will see, this Fermi-velocity is important for extracting the entire DM kinetic energy in the scattering process. Third, when the metal becomes superconducting, a gap develops between the electrons in the Fermi sea and the states into which the electrons can scatter. The size of this gap is small (of order a meV), but it is crucial for controlling the noise. It effectively allows the decoupling of energy deposited in vibrations in the lattice (phonons, dominated by thermal noise) from energy deposited directly into an electron in a hard scatter.
The energy from DM is deposited into the detector when the DM interacts with one of the electrons in the ground state of the system, namely in a Cooper pair. When the energy deposited is larger than the Cooper pair binding energy (related to the gap), the pair is broken, and two quasiparticles are excited above the gap. These excitations are then detected via a mechanism that we describe in detail in section 2.
The outline of this paper is as follows. In section 2 we outline the basic notion of the detection method, and present two concrete detector designs. In section 3 we describe the treatment of DM scattering in a Fermi-degenerate medium. Section 4 contains various constraints on DM scattering with electrons. In section 5 we discuss several particular models: scalar and vector mediation is considered in section 5.1; a kinetically mixed hidden photon is considered in section 5.2, including in-medium effects; dipole interactions of DM are tackled in section 5.4; and milli-charged DM is discussed in section 5.3. We conclude in section 6.

Detection principle and design
We begin by presenting the underlying idea behind our proposed detection method. After establishing the basic notion for detecting O(meV) energy depositions, we present concrete experimental detector designs that could be sensitive to this energy range.

Detection principle
When searching for DM with mass heavier than ∼ 10 GeV with elastic scattering, nuclear targets have three main advantages: first, DM elastic scattering has a rate that scales as the reduced mass of the DM-target system, µ r , which suppresses the scattering rate on electrons compared to that on nucleons. Second, the maximum deposited energy in an elastic recoil off of a target at rest is which is maximized when the target has mass equal to that of the DM X. Here, m T is the target mass and v X ∼ 10 −3 the DM velocity. Thus, 100 GeV DM produces nuclear recoils of O(10 keV), but e − recoils of only O(eV). Third, backgrounds in direct detection experiments, such as Compton scattering, feature mainly an electron ionization component; thus, discriminating nuclear recoils from electromagnetic activity acts as a major discriminant for reducing backgrounds.
As the DM mass drops below the mass of the nuclear target, around 10's of GeV, eq. (2.1) indicates that the deposited energy is suppressed by m 2 X /m 2 T compared to the case of DM heavier than the target. For example, sensitivity to 1 eV nuclear recoils allows reach to 100 MeV DM. Searching for such 1 eV nuclear recoils from 100 MeV DM scatterings has motivated SuperCDMS's push to lower thresholds [17][18][19]. In addition, utilizing a lighter nuclear target, such as 4 He, is also advantageous in searching for lighter DM [28,29].
To access even lighter DM, electron targets are preferred. In this case, an energy deposition sensitivity of 1 eV corresponds to probing DM models with mass down to roughly 1 MeV. Of course, this sensitivity can only be achieved if the energy deposit exceeds the binding energy of the electron. In a xenon atom, the binding energy of the outermost electron is 12 eV, while in germanium the band gap is 0.7 eV. Thus, the binding energy of electrons in atomic targets and semi-conductors fundamentally limits access to DM candidates with mass below an MeV. To access such candidates, we need a material with a gap smaller than m X v 2 X /2, which corresponds to O(meV) energy for DM at the warm DM limit of O(keV). Metals (including metals in a superconducting phase) and superfluids are examples of materials that feature a small or no gap, and as such can be appropriate.
Thus materials with a O(meV) gap and sensitivity to O(meV) energy depositions may allow for the detection of DM at the O(keV) mass scale. However, even when it is energetically possible for DM to catalyze a reaction -when the DM kinetic energy exceeds the gap of the material -kinematics may still forbid the scattering. To see this, consider the deposited energy on a target (electron or nucleus) in terms of the momentum transfer of the process q: where v i,T is the initial velocity of the target and δ (defined to have a positive sign for bound electrons) is the gap of the system. The first term in eq. (2.2) is the usual energy deposition for elastic scattering on targets at rest, just like eq. (2.1). The third term takes JHEP08(2016)057 into account that fact that DM may catalyze an inelastic process, releasing the binding energy of the target electron. The second term is the one we wish to focus on: it is the effect of the target's initial velocity on the total amount of DM kinetic energy that can be absorbed. Even if the first term is small (e.g. ∼ µeV for a keV DM scattering on an electron), the second term may allow extraction of the entire kinetic energy of the DM. In the metal and superfluid targets we are most interested in here, the velocity of the target is a property of the target ground state, and is due to Fermi statistics. (We note that semiconductor and atomic targets feature an electron velocity of similar size.) To illustrate this point, consider two relevant limits, when the DM is heavier or lighter than the target. (For the purpose of this illustration, we neglect the band gap, assuming it is substantially smaller than the DM kinetic energy.) When the DM is much heavier than the target, the center of mass frame is approximately the DM rest frame. In this frame, the collision between the DM and the target barely changes the DM velocity, and the target initial and final state velocities have the same magnitude but are in opposite directions. In the lab frame, the target velocity changes at most from v i,T to (v i,T + 2v X ), with opposite direction. Thus the maximum energy deposition can be written as E max On the other hand, in the limit that the DM is much lighter than the target, the maximum energy deposition is obtained when the DM is fully stopped by the target. For example, a target with velocity (0, v 2 i,T − v 2 X /4, v X /2) can fully stop a DM particle with velocity (0, 0, v X ), and the momentum transfer in this case is simply the DM initial momentum, m X v X . Since the deposited energy is approximately 1 2 m X v 2 X , the experiment must have meV energy resolution in order to be sensitive to keV mass DM.
What is the typical target velocity in the (nearly) gapless materials we consider? In a metal like aluminum, the valence electrons have a Fermi momentum p F = 3 keV, giving rise to a Fermi velocity for electrons of v F 10 −2 . Note that this effect is purely due to Pauli blocking and Fermi statistics in a degenerate medium at low temperature. Superfluids such as Helium-3, where the nucleus has half-integer spin, also display Fermi degeneracy. In the case of Helium-3, the Fermi energy is E F 4 × 10 −4 eV, giving rise to a Fermi velocity of the Helium-3 nucleus of order v F 10 −6 . For a typical momentum transfer of m X v X , the second term in eq. (2.2) never dominates in Helium-3, however, and instead the scattering proceeds via the ordinary nuclear recoil process.
Note that the electron velocity does play a role in DM scattering off electrons in a semi-conductor or noble gas as well [22]. In such scatterings, for DM heavier than an MeV, the electron velocity is not a necessary ingredient to catalyze the process and extract all of the DM energy in the scattering, but nevertheless it does impact the kinematics, since the electron velocity, v i,T ∼ α Z with Z the electric charge of the nucleus, is larger than the velocity of the DM.
An alternative route towards detecting DM as light as a keV is to take advantage of the gap, with inelasticity catalyzing the scattering. This is evident from eq. (2.2): even if the first two terms are below a meV, as long as the kinetic energy of the DM exceeds the gap, the DM energy may be absorbed by exciting an electron above the gap. For the metal targets we are interested in for the rest of this paper, when the metal enters JHEP08(2016)057  Figure 1. Schematic designs for superconducting detectors that are sensitive to DM-electron scattering. Left: quasiparticles produced by a recoiling e − in a large aluminum arbsorber are collected by tungsten quasiparticle collection fins and then their energy is sensed by a TES. Right: athermal phonons produced by a recoil e − in a large tantalum absorber are collected by aluminum collection fins and then their energy is sensed by a TES. the superconducting phase, a sub-meV gap appears. Since this gap is below detectable energies for the devices we consider, we ignore the presence of the gap, and focus on elastic scattering. Also note that while the presence of the superconducting gap is not important for the scattering process itself, its existence means that athermal phonons and quasiparticles have very long lifetimes, and as such can potentially be collected before they thermalize. Thus in the systems we consider, detection of DM operates via the breaking of Cooper pairs in a superconducting target. We consider this idea in more detail next.

Detector design with milli-eV sensitivity
Our detector concept is based on collecting and concentrating long lived athermal excitations from DM interactions in a superconducting target absorber onto a small volume (and thus highly sensitive) sensor. The collection and concentration of long lived excitations is a general concept that has been a core principle of detector physics, from ionization in semiconductor CCDs to athermal phonon collection in CDMS. Here we propose that this general detection philosophy be applied in large volume (very pure, single crystal) superconductors to search for DM with mass as low as the warm DM limit of a keV using standard superconducting sensor technology that has been pushed to its ultimate theoretical sensitivity. A schematic of two proposed detector concepts for light dark matter, that we describe in greater detail through the remainder of this section, is shown in figure 1.
Detection of dark matter in such detectors is comprised of a three part process: • Dark Matter Scattering on Target Absorber and Subsequent Excitation Production.
A DM particle scatters off an e − in the target metal or superconducting absorber. In subsequent interactions, the recoil energy is converted into long lived athermal phonons and quasiparticles.

JHEP08(2016)057
• Collection of Excitations. The resulting excitations must be collected and concentrated onto a small volume (and thus very sensitive) sensor; this is typically done via 'collection fins' on the surface of the absorber that efficiently collect the energy from the excitations in the absorber.
• Measurement of Excitation Energy/Quanta. After collecting the excitations, they must be measured. In the devices we consider, this is typically done via a transition edge sensor (TES) or a microwave kinetic inductance device (MKID).
In the remainder of this section we describe in detail 2 potentially feasible detector designs; each of the following three subsections is dedicated to absorption and excitation production, collection and measurement. The theorist interested primarily in a calculation of the reach of these detectors, as well as theories of DM that could be detected with such devices, can move directly to section 3.

Excitation production in superconductors
Superconductors are excellent candidates for detection of light dark matter because the absence of any unoccupied electronic states within the superconducting band gap ∆ of the Fermi surface means that both quasiparticles near the gap edge and athermal phonons with energy below the gap are long lived, and thus potentially amenable to measurement.
The valence (conducting) electrons at low temperature in a metal are well-described by a Fermi-degenerate distribution. In a metal like aluminum, the Fermi energy is 11.7 eV, and the corresponding number density of conducting electrons is This is the reservoir of electrons that are available to become superconducting once the temperature of the metal drops below the critical temperature T c and the electrons enter the superconducting phase; it is also these electrons that are available for DM scattering in the target metal. The superconducting phase is entered when it is energetically favorable for near Fermi energy electrons to bind into pairs, known as Cooper pairs. This ground state of the system is then highly correlated, and free electrons are no longer the correct degrees of freedom of the system. (For a review of superconductivity, see e.g. ref. [30].) The binding energy 2∆ of these pairs is typically quite small, e.g. of order 0.6 meV in aluminum, and their correlation length ξ 0 = v F /(π∆) macroscopic, e.g. of order a micron in aluminum. A DM particle interacts with these valence electrons bound into Cooper pairs in the superconducting target. For the purpose of the rate calculation carried out in section 3, the important point is that as long as the energy deposited in the DM-electron scattering well-exceeds the binding energy of the Cooper pair, the DM-electron scattering rate can be approximated via energy deposit onto free electrons in a Fermi degenerate free-electron sea [30]. As the energy deposited approaches the Cooper pair binding energy, a 'coherence factor', analogous to a form factor, takes into account the impact of the coherence from the Cooper pairing phenomenon; we discuss this further in the next section.

JHEP08(2016)057
The initial interaction between the Cooper pair and the DM creates two quasiparticle excitations from the ground state. As discussed in ref. [31], for E 100 meV, thermalization of these excited quasiparticles occurs predominantly via athermal phonon production. Since the phonon phase space scales as E 2 , the produced phonon distributions tend to be hard; on average most of the excess energy of the quasiparticle is converted to a single athermal phonon quanta. Then, as long as the produced athermal phonon has E > 2∆, it will break an additional Cooper pair, and so on. At the end of the cascade process, the total kinetic energy of the recoil has been converted into ∼ 60% quasiparticle potential and kinetic energy, with the remaining energy in athermal phonons with a distribution that is strongly near the 2∆ cutoff in the superconductor [32].
Once created, what are the dynamics of these excitations? Is it possible to collect them before they thermalize? We address this important issue for quasiparticles and athermal phonons next.
1. Quasiparticle dynamics. Extremely pure single crystal aluminum is very unique in that electronic excitation scattering lengths of ∼ 1.5 mm have been measured at temperatures of 4 K (with residual resistance ratio RRR ≡ R(300K)/R(4.2K) 10 5 ) [33]. Furthermore, these measured scattering lengths should underestimate the scattering length of aluminum quasiparticles in a crystal of similar quality at dilution fridge temperatures of ∼ 6mK for two reasons: • The thermal phonon population is significantly smaller at temperatures of 6mK compared to 4 K, and thus phonon up-scattering rates, which can be nonnegligible at 4 K in extremely pure aluminum crystals, are completely suppressed.
• Quasiparticle scattering rates off of impurities are suppressed compared to those of normal electrons.
Thus, for a 5 mm aluminum crystal, quasiparticle propagation is essentially ballistic.
The lifetime of quasiparticles in very high quality single crystal aluminum has not been measured to our knowledge. A 2 ms quasiparticle lifetime has been measured in 100 nm thick aluminum MKIDs [34]. However, this is probably too conservative since the RRR of the aluminum films used was only 3.3 and the MKID lifetime was found to have a strong dependence on dislocation density. Furthermore, the surface to volume ratio in these films is much larger than those found in the 5 mm cubic absorbers that we are considering and thus any recombination/trapping on the surfaces is enhanced in this thin film device compared to what we would expect. Thus, for purpose of our calculations, we'll assume a 20ms lifetime for quasi-particles in single crystal aluminum. Together with the quasiparticle group velocity [35],

JHEP08(2016)057
(depending on the quasiparticle temperature T ), this suggests that for a 5 mm aluminum absorber, quasiparticles would bounce off the single crystal surface > 10 6 times before recombination.
2. Athermal phonon dynamics. Due to the lack of electronic states within the superconducting bandgap, elastic and inelastic scattering of athermal phonons with energy below 2∆ through electron-phonon interaction is impossible to lowest order in the absence of quasiparticles. Consequently, athermal phonon dynamics in high quality single crystal superconductors at T T c mirror those found in detector grade semiconductor crystals.
In particular, elastic scattering will be driven by isotopic and impurity scattering [36]. Since natural aluminum (Al) and tantalum (Ta) are almost entirely composed of a single isotope, elastic isotopic scattering should be negligible in these superconductors. Simultaneously, impurity scattering in these materials is minimal due to the use of very pure single crystals (float zone refining is easily implemented in single crystal metals [33]). Thus -as was the case for excited quasiparticles in aluminum -athermal phonons in aluminum and tantalum should be ballistic for O(1 cm)sized absorbers. We note that this conclusion is seemingly in conflict with athermal phonon propagation studies done in single crystal lead [37,38], however this is quite expected. 1 In the bulk of these crystals, the lifetime of the athermal phonon will be limited by phonon anharmonic decay [39], in which a phonon splits to two phonons. Since the 3 rd order elastic constants (effectively, the coupling constant for this phonon splitting process) in aluminum are similar to those found in germanium and silicon, we expect the anharmonic decay lifetime for a 4 K phonon to be of O(1 s). (The equivalent numbers for tantalum could not be found in the literature, but we assume they are of similar size.) We thus estimate more than 2 × 10 5 surface bounces for athermal phonons in a 5 mm absorber. Of course, with such a large number of potential surface bounces, phonon down-conversion at the surface is most likely the dominant thermalization process. Such processes depend critically on the exact surface preparation and are thus difficult to estimate. For example, SuperCDMS has some evidence that a 30 nm amorphous silicon layer on a germanium crystal down-converts athermal phonons every ∼ 250 bounces but on bare germanium surfaces there is only a lower limit, of 1250 bounces [40]. For the purpose of estimating athermal phonon collecting 1 The reasons are as follows: (1) The source of athermal phonons for their propagation measurements was a 30-70 K thermal hot spot in a copper film on the surface of the lead crystal. Consequently, there was also a large non-equilibrium source of quasiparticles that was shown to completely dominate scattering [38].
(2) The athermal phonon scattering measurements were done at 3 < Tc/T < 5, and thus there was a non-negligible fraction of equilibrium quasiparticle scattering that would not be present at Tc/T ∼ 100.
(3) Pb naturally has large isotopic scattering not present in aluminum. (4) Isotopic and impurity phonon scattering rates scale as the inverse Debye temperature cubed [36]. As a result, even for similar impurity levels the phonon scattering rate would be suppressed by two orders of magnitude in aluminum compared to lead.

JHEP08(2016)057
detector sensitivity in both aluminum and tantalum absborbers, we will use this bare germanium lower limit.
In summary, the quasiparticle excitations in very pure single crystal aluminum and the athermal phonon excitations in aluminum and tantalum are very likely to have characteristic lifetimes 2 ms and scattering lengths potentially much larger than a mm. As a result, they have excellent potential as target material for excitation-sensitive detector technology.
The scattering length for both quasiparticles and athermal phonon bounds the size of the absorber, which we consequently take to be of order ∼ (5 mm) 3 . A large number of these small detectors would then be placed in parallel in order to obtain large exposure.
Having established the longevity of excitations in the absorber, we move to discuss collection and concentration of these excitations.

Collection and concentration of long lived excitations
The second step of these long-lived excitation detectors is to collect and concentrate the excitations from the absorber into the much smaller sensing region. In semiconductors like germanium and silicon, this is largely trivial; since the electron-hole pair are electrically charged, they can be drifted towards the sensor region by biasing the sensing region at an appropriate voltage. Such techniques are unfortunately impossible within a superconductor due to the perfect shielding of both electric and magnetic fields. (Furthermore, athermal phonons are electrically neutral.) We follow the spirit of design concepts first laid out in superconducting tunnel junctions (STJs) [41] and later by CDMS using athermal phonons [42]. We simply allow the long lived excitations to randomly propagate within the absorber, and instrument a small fraction of the overall surface of the absorber with a material that has a high probability, f trap , of collecting/trapping the excitation upon contact.
When propagation is ballistic, an excitation will, on average, be collected after ftrap bounces, corresponding in a cubic absorber to a collection time of [43] where A absorber , V absorber , and A collect are the total absorber surface area, the total absorber volume, and the instrumented collection area respectively, with |v| the average excitation velocity. The excitation collection process competes with annihilation processes (either phonon anharmonic decay or quasiparticle recombination), and thus the average excitation collection efficiency is given by The benefit of very large excitation lifetimes (τ life ) is now clear. It allows one to achieve high excitation collection efficiencies even for very large ratios of V absorber /A collect . Since sensor sensitivity scales with the size of the sensor and thus with A collect , large excitation lifetimes allow one to simultaneously satisfy the requirements of large exposures and very low energy thresholds in optimized devices.

JHEP08(2016)057
When excitation propagation is diffusive, collection probability and collection times will be highly dependent on the location of the generated excitation. In particular, excitations generated far from an instrumented surface will have long collection times and potentially very poor collection probabilities. To minimize these issues, as we already commented above, we require that the absorber be of such a size that the excitations are ballistic; thus the absorber must be relatively small (a cm or smaller).
For quasiparticle collection, the standard collector material is either a superconductor with transition temperature T c,collect that is less than that of the absorber, T c,absorber , or simply a normal metal. In both cases, the collector material is placed in direct contact with the absorber, producing a proximitized region where the critical temperature (and thus the gap) is between that of the absorber and the collector. Within this suppressed bandgap region, a portion of the quasiparticle's original potential energy will be converted into kinetic energy, and the likelihood of inelastic phonon production will be significantly enhanced. If this occurs, the quasiparticle will then be trapped within the collection volume.
The thickness of the collection film collect is a result of optimizing two competing effects. On the one hand, one would like large f trap . Since the phonon production rate scales as the excess quasiparticle energy cubed [44], and the bandgap suppression scales roughly linearly with thickness [45], the phonon dropping rate scales as ∼3 collect . Furthermore, since the amount of time spent in the film scales as collect , we have f trap ∝ ∼4 collect . (Such large power law scalings of f trap with collect are consistent with previous experiments [46]). On the other hand, we would like to maximize the efficiency of energy collection and thus minimize the energy of the phonons that are released in the trapping process, which suggests thinner films. In the quasi-particle collection interface between 35 nm tungsten and 350 nm aluminum films in SuperCDMS devices, f trap is within the range of 10 −3 -10 −4 [47]. This suggests that f trap ∼ 0.1 could be achievable with 100-200 nm thick collection films, which leads to an estimate of 87% for f collect , if one places twelve 225 µm 2 size quasiparticle collection fins on the surface of the aluminum absorber. We estimate that the average quasiparticle potential energy remaining after sub-gap phonon emission for trapping, f E Remain , will be greater than 90%.
For athermal phonon collection, we again use a superconducting collecting film with T c,collect < T c,absorber , but which is now electrically isolated from the large-volume superconducting absorber by an insulating layer (SiNx, SiOx and Al 2 O 3 are all viable options). With this configuration, athermal phonons with energies 2∆ collect < E < 2∆ absorber will ballistically travel throughout the absorber, but will annihilate within the collector film, storing their energy in quasiparticles. To collect a large fraction of the athermal phonons, the difference in the T c between the two materials should be significant. Some viable possibilities are tantalum(absorber)/aluminum(collector) and aluminum(absorber)/titanium(collector). This greater freedom in the absorber/collector duos for athermal phonon collecting detectors is due to the fact that the requirements on the absorber in this case are much less constraining since phonon scattering lengths and lifetimes should be excellent for the vast majority of single-element superconductors.
After the athermal phonons are converted into quasiparticles within the aluminum collection fin, they diffuse until they are absorbed by the connected TES -just as in JHEP08 (2016)057 CDMS athermal phonon detectors. This diffusion processes introduces an additional energy loss mechanism due to quasiparticle trapping, that has been well studied for SuperCDMS detector geometries [47]. We estimate the quasiparticle collection efficiency to be ∼ 0.65 for the aluminum collection fin geometry used in the proposed athermal phonon detector (see table 2).

Measurement
With the athermal excitations from the target now concentrated, all that is necessary is to read them out with a sensor of the appropriate sensitivity; in essence, one can just repurpose a single-infared-photon-sensitive detector. Unfortunately, the required O(1 meV) sensitivities have not yet been achieved experimentally with any technology. For both Transition Edge Sensors (TES) and Microwave Kinetic Inductance Devices (MKID) though, such sensitivities are theoretically possible. Furthermore, in both cases, engineering solutions (though extremely challenging) have been proposed, which could allow the theoretical sensitivities to be realized. Below, we detail the engineering challenges for the TES.
The TES is a superconducting film that has been artificially stabilized through electrothermal feedback at an operating point just within the superconducting transition. Biased in this manner, very small changes in the temperature of its electronic system can produce substantial changes in resistivity, which are then measured [48]. The theoretical energy resolution (squared) of the TES is given by where S p,tot (ν) is the total noise referenced to TES input power, T c,TES is the transition temperature of the biased system and C is the heat capacity of the TES. Here α is the unitless measure of sensor sensitivity defined as T R ∂R ∂T at the TES operating point, and can take values in the range of 20-200, depending on the TES film. Finally, n Tb is the temperature scaling exponent on the thermal power which flows between the TES and the heat bath. To gain intuition into eq. (2.7), note that the energy variance for a heat capacitor C coupled to a thermal bath via conductance G is 4k B T 2 C. Recognizing that the heat capacitance scales as C = γV TES T , (2.8) where γ is the specific heat coefficient for the appropriate metal, and V TES is the volume of the TES, we find that the energy resolution scales as In table 1 we list the measured energy sensitivity of three of the most sensitive TES bolometers/calorimeters that exist today, along with their physical dimensions and operating temperatures T c,TES . What is immediately clear is that none of these devices have attempted to minimize both the TES volume and T c,TES concurrently, and thus substantial sensitivity increases are conceptually possible. For an estimate of the potential gains, we have scaled these devices using eq. (2.9) to the proposed TES geometry and operating For the bolometer of ref. [51], energy resolution is estimated as the power noise multiplied by where τ eff is the sensor fall-time. For ref. [50], energy sensitivity scalings are estimated for the device with T c = 100 mK as well as for the B-Field suppressed value of 50 mK.
temperature for our quasiparticle collection detector, which would be six 9 mK tungsten TES in parallel with dimensions of 1 µm × 24 µm × 35 nm each, for a total TES volume of 4.2 µm 3 (see table 2). The resulting scaled energy resolutions σ scale E are given in the right-most column of table 1. As is evident, O(meV) energy sensitivities seem feasible.
Unfortunately, improvements in TES sensitivity to low energy recoils via eq. (2.9) are naturally accompanied by increased sensitivity to environmental noise. Since the thermal power flow between the TES and the bath scales with a power of n Tb ∼ 5, the bias power required to keep the TES within transition also scales as T n Tb c . To give a sense of scale, the TES for our proposed quasiparticle detector is estimated to have a bias power of 8.3 × 10 −20 W, nearly 3 orders of magnitude smaller than that of current devices (table 1). Ideally, the bias power is predominantly supplied by the TES readout electronics, but this is certainly not necessarily the case. Vibrations from cooling machinery (such as pulse tube cryocoolers and turbo pumps for dry dilution refrigerator systems, or 1 K pot vibrations in wet systems) could dissipate power within the TES. Likewise, thermal radiation from poorly shielded higher temperature stages could be absorbed by the TES. Finally, electromagnetic interference (even beyond the sensor bandwidth) can be coupled into the TES via the wiring. In summary, the constraints on DC environmental power loading are ×10 3 more strict than levels currently achieved.
Both electromagnetic interference and vibrational environmental noise sources will naturally have fluctuating components within the TES sensor bandwidth as well, and thus as one decreases the fundamental thermal fluctuations between the TES and the bath, these sources could begin to dominate and suppress the TES sensitivity. Roughly, a decrease in environmental power noise by a factor of 50 from levels achieved today are required in order to meet the required TES performance specifications.
An automatic benefit of operating a TES at the low temperatures of the proposed detectors, is that the sensor fall-times (which are ∝ C/G) naturally become very long (sensor bandwidths become very short) and so match the very long excitation-collection timescales that are envisioned, of O(10 msec). Thus, the problem of bandwidth mismatch JHEP08(2016)057 between the TES sensor and the excitation collection time, which degrades the energy sensitivity of current SuperCDMS [52] and CRESST detectors [53], is naturally suppressed.
Device specifications and estimated performance for both the proposed quasiparticle and athermal phonon excitation detectors are shown in table 2. In the table, we have assumed that the detector trigger threshold is 6 times the estimated detector baseline energy resolution, σ E D . Of course, σ E D is just the TES baseline resolution (σ E TES ) divided by the efficiency factors for collecting and concentrating excitation energy in the TES that were discussed in the text.

Backgrounds
Solar neutrinos are an irreducible background, with rate (per unit mass per unit time) where F ν is the neutrino flux and n T the number density of the target. As we focus on low energy depositions in the detector, the dominant contribution comes from pp neutrinos [54,55] scattering on nuclei. This assumes that an O(1) fraction of the energy deposited in nuclear recoils is converted into quasiparticles in the detector. The rate is shown in figure 2, for a few sample nuclei. We find that for an aluminum target, the integrated neutrino background for a kg·year is less than 1 event for nuclear recoils between 1 meV and 1 eV, and is 3 events for nuclear recoils between 10 meV and 10 eV. We include the 3 background events where relevant in extracting DM limits, accordingly. The expected U/Th/K compton background of 13 event/keV/kg/year for Si detectors (a material with similar stopping power to Al) within the proposed SuperCDMS SNO-LAB cryostat, plus the background due to the beta decay of 3 H produced via cosmogenic spallation during detector fabrication at sea level (60 days assuming 125 atoms/kg/day production) are also shown in figure 2 and found to be sub-dominant [56]. At first glance this might seem surprising because minimization of radiogenic backgrounds is the primary design driver in high mass dark matter direction detection. However, there are two reasons why radiogenic backgrounds are of secondary importance for light mass dark matter detection. First, the low energy coherent neutrino scattering background from pp neutrinos is much larger than the background produced by atmospheric neutrinos within the high mass dark matter region of interest. Secondly, all of the radiogenic backgrounds (comptons, 210 Pb decay products, 3 H) have characteristic energy scales which are much larger than the light mass dark matter region of interest (< 10 eV) and thus there is very little overlap between radiogenic backgrounds and light mass dark matter recoil signals. into a superconducting state at low temperature, a ∼ meV gap opens up above the Fermi surface, blocking DM-electron scattering for energy depositions below this gap. For energy deposits well above the gap, the scattering is simply characterized by allowable momentum configurations of DM-electron scattering that are consistent with Fermi statistics and Pauli blocking. As the energy deposits drop and approach the gap, an additional factor that takes into account the presence of the superconducting gap, a so-called coherence factor -similar to a form factor -kicks in [30]; this factor depends on energy and has an effect JHEP08(2016)057 only near threshold. Since the energy thresholds we consider are always above the 0.3 meV aluminum superconducting gap, we neglect the coherence factor in what follows. We thus approximate the electrons in the superconducting target as a free Fermi-degenerate gas.

JHEP08(2016)057
The most important property of the Fermi-degenerate metal or gas to properly incorporate is the phase space suppression of Pauli blocking -in a Fermi-degenerate medium, the DM must deposit enough energy to knock an electron out of the Fermi sea and into the continuum above the Fermi surface. We closely follow the discussion in ref. [57], and reformulate their calculations for non-relativistic DM-electron scattering. We denote the 4-momentum of DM initial and final states by P 1 and P 3 , the initial and final states of the electron by P 2 and P 4 , and the momentum transfer q = (E D , q). The scattering rate for a DM particle can be estimated via with E D the deposited energy, |M| 2 the squared scattering matrix element summed and averaged over spins, and f i ( is the Fermi-Dirac distribution of the electrons at temperature T . S(E D , q) characterizes the Pauli blocking effects of the process at hand. If, for instance, the scattering converts an electron to a different final state particle which does not exhibit Pauli blocking, (1 −f 4 (E 4 )) should be dropped in S(E D , q), JHEP08(2016)057 and the integrals d 3 p 3 and d 3 p 4 in eq. (3.1) reduce to the ordinary 2-to-2 scattering phase space integral. Analytically, S(E D , q) is found to be and µ is the chemical potential, identified as E F at zero temperature.
In the limit of T → 0, we have z → +∞ and ξ → 0, yielding with Θ the Heaviside theta function. We note that this limit is only valid when both e − < µ and (e − + E D ) < µ (the small region where only one inequality is satisfied is unimportant for the rate estimation). In what follows we compute the rate numerically using the full eqs. (3.2) and (3.3) at temperature much lower than the gap, in order to capture the entire kinematic range properly. Converting d 3 p 3 to energy E D and momentum transfer q, the total interaction rate of DM-electron scattering, per unit mass per unit time, is Here f MB is the velocity distribution of DM, which we take to be a modified Maxwell-Boltzmann distribution [58], with the normalization factor with z = v esc /v 0 . We use rms velocity v 0 = 220km/s and cut-off at the escape velocity v esc = 500 km/s. ρ is the mass density of the detector material, ρ X /m X the DM local number density, with the DM mass density ρ X = 0.3 GeV/cm 3 . A typical Fermi velocity is v F = O(10 3 ) km/s v esc , leading to v rel v F in eq. (3.6). As already discussed in section 2.1, when the DM is lighter than the electron, there is always an electron configuration in the Fermi sea that can fully stop the DM, and the energy cutoff is determined by the incoming kinetic energy of the DM, yielding E D ranging from 0 to 1 2 µ X v 2 esc . In contrast, when the DM is heavier than the electron, no electron can fully stop it, and the maximal deposited energy in this case is E max

JHEP08(2016)057
Intuitively, one expects the Pauli blocking effect to suppress the total rate by ∼ E D /E F , as this indicates the relative size of the shell of electrons available for energy E D deposited in the scattering. We find that the complete Pauli blocking computation is well captured by this naive estimation when E D is small compared to Fermi energy. As E D approaches E F , more of the electrons can participate in the scattering, and the effect of Pauli blocking is suppressed. As an illustration, we study the behavior of the fraction of electrons participating in the scattering when the deposited energy is maximal, which is directly related to the DM mass. The fraction of participating electrons is plotted in figure 3, and the numerical result agrees very well with the expectation: when DM is light, the maximal energy deposition is small, and the fraction of active electrons grows linearly with DM mass. When the DM mass is O(MeV), the DM kinetic energy is comparable to Fermi energy, and the linear growth approximation fails.
The scattering cross section between DM X and free electrons with a mediator φ is given by where α i ≡ g 2 i /(4π), g i is the coupling of φ to i = e, X, µ eX is the reduced mass of the DM-electron system, and q is the three-momentum transfer in the process, determined by the kinematics of the detection process. (Here we keep only the contributions from the 3-momentum transfer q since the energy transfer in a t-channel non-relativistic scattering JHEP08(2016)057 is much smaller.) This scattering cross section is related to the matrix element squared in eq. (3.1) through We define two related reference cross sections,σ DD , corresponding to the light and heavy mediator regimes:σ where v X ∼ 10 −3 is the DM velocity. In the above, the reference momentum q ref is chosen for convenience as a typical momentum exchange. It is worth noting that for very light DM that deposits E D energy, the momentum transfer in the process can be larger than this q ref by a factor of v F /v X . In-medium effects of massless mediators, which can alter eq. (3.9), will be addressed when relevant.
To establish a sense of a number of events expected, in figure 4 we plot the differential rate per kg·year as a function of deposited energy E D , for several benchmark points. The behavior of the curves can be readily understood. As we have seen, the maximum energy deposition is controlled by the DM mass, independently from the mediator. When the mediator is heavy compared to the momentum transfer of the process, the rate is peaked at high energy depositions. This is because, when the energy deposition is well below Fermi energy, the larger the energy deposits, the larger the fraction of participating electrons. For very light mediators, the rate is dominated by the minimal momentum transfer in the process, which is controlled by the detector energy threshold.

Cosmological, astrophysical and terrestrial constraints
Having established the DM interaction rate in our proposed detectors, we now consider the types of constraints such DM is subject to in order to determine whether DM candidates consistent with all constraints are within reach.

Self-interactions
DM self-interactions bound g X via a constraint on the scattering cross section weighted by the momentum transfer, In the Born regime, where α X m X m φ , the analytic perturbative result [59,60] for attractive and repulsive forces is  which reduces in the heavy mediator limit of m φ m X v, as expected, to the contact operator form, For very light mediators in the classical regime, where m X v m φ , the solution to the classical equations of motion in repulsive and attractive potentials (see e.g. [59] and references therein) reduces to in the limit of β 1, which will always be applicable to our light (but massive) mediator case. Here we have taken Dirac DM with interactions via a vector or scalar mediator; a Majorana or real scalar DM particle would have a factor of 4 larger scattering cross-section.
Bullet-cluster constraints [61][62][63] along with recent simulations which reanalyze the constraints from halo shapes [64,65], limit the DM self-interaction cross section to be roughly depending on the relevant velocity; further details can be found e.g. in ref. [66]. The self-interaction constraints will be most relevant when discussing light mediators, where the transfer cross section is proportional to 1/v 4 . In order to be conservative, in later JHEP08(2016)057 discussions in section 5, we use σ T 10 cm 2 /g with v ∼ 10 −4 to impose an upper bound on g X . For the very light mediator regime, this roughly translates to requiring where ln β −1 varies by a factor of at most a few in the region of interest.

Kinetic decoupling
The couplings of a light mediator m φ eV to DM and to electrons are constrained by CMB measurements. This is because if DM is in kinetic equilibrium with the the photon-baryon plasma during the recombination epoch, DM density fluctuations can be washed out via Silk damping [67] and the baryon acoustic peak structure can be altered. Closely following ref. [68], we require the relaxation rate of energy transfer is slower than the expansion rate of the universe: where µ bX is the reduced mass of the baryon-DM, and λ cut is the screening length for the baryon plasma; for massive mediators, this is set by 1/m φ , while for photon exchange, it corresponds to the Debye screening mass, λ D = T /(4α EM πn e ) with n e the electron density. For light mediators, eq. (4.7) is to be evaluated at the time of recombination, T = T rec 0.26 eV. For heavy mediators having m 2 φ T rec m X , ensuring the interaction is out of equilibrium at recombination is not enough -we must ensure that DM be decoupled from the plasma when the momentum transfer is of order the mediator mass, leading toT = m 2 φ /(4µ bX ) in eq. (4.7) above. In the discussion in section 5 we compare this constraint on α e α X against others and require the strongest of them holds. Of course, kinetic decoupling need only be enforced when the DM and/or the mediator have masses below a few MeV, and at temperatures below a few MeV. The constraint is most relevant for light mediators, where it roughly reads (α X α e ) light kin. dec.
within O(1) factors, for the entire mass ranges of interest to us, where we use at recombination n e,rec = n p,rec 2.4 × 10 −39 GeV −3 . We note that in the above we have considered kinetic decoupling between the DM and the SM particles it scatters with via exchange of the mediator. In the case of a light mediator, the mediator itself (rather than the DM) may be brought into thermal equilibrium via Compton-like processes; requiring the mediator is out of equilibrium during BBN results in a constraint on g e weaker than those from stars we consider next.

Stellar emission
If DM scatters off electrons via the exchange of a mediator, the mediator may be emitted from stellar objects and generate excess cooling. A variety of emission processes can occur -amongst them Compton-like processes, bremstrahlung off electrons, Primakoff emission and plasmon decay/conversion. We now summarize the relatively model-independent constraints on the mediator coupling to electrons g e that are generally applicable to a scalar mediator; the relevant constraints on a non-kinetically mixed vector differ by an O(1) factor [69].
For mediator masses beneath a keV, we find that the dominant stellar constraint comes from Horizontal Branch (HB) star cooling. The bremstrahlung process imposes the strongest bound [69,70], which is stronger by roughly a factor of 3 than that from the Compton-like process [69]. We have verified that the loop-level induced coupling of φ to photons, which generates additional emission via the Primakoff process, yields a constraint on g e that is ∼ 3 orders of magnitude weaker than eq. (4.9). The above constraint eq. (4.9) is, in principle, applicable to all mediators with mass beneath the typical stellar temperature of ∼ 10 keV. We emphasize, however, that this is not always the case; in particular models, such as the kinetically mixed light hidden photon, the constraint on the coupling g e is much weaker. This will be addressed separately in section 5.2.
For mediator masses exceeding several ten's of keV, the relevant stellar environment becomes the hotter supernovae, though in these dense objects, stellar constraints can be lifted due to trapping effects. There can be model-dependence in which process controls trapping, but decays of the mediator allow electron couplings of order g e 10 −6 at the very least (see e.g. ref. [71]) for mediator masses above a few hundred keV. Here terrestrial experiments can play a vital role.

Terrestrial constraints
Terrestrial experiments are complementary probes to the electron-coupling constraints considered thus far. Restricting to processes that do not require coupling to additional SM particles beyond electrons, the relevant experimental bounds come from measurements of the anomalous magnetic moment of the electron (g − 2) e [72][73][74], beam dump experiments such as E137 [75], and low energy e + e − machines such as BaBar [76,77].
When the mediator is heavier than twice the electron mass, it can decay visibly to a pair of electrons; if it is also heavier than twice the DM mass, invisible decays open up as well. Although a broad region of parameter space for m φ MeV decaying either visibly or invisibly is constrained by terrestrial searches, a viable window of parameter space with fairly large couplings opens up even for mediator masses as light as a few MeV. This is due to supernova trapping effects, in a similar way to the case of the well-studied kinetically mixed hidden photon. For a sense of the size of the allowed couplings, a m φ ∼ 10 MeV

JHEP08(2016)057
[100 MeV] mediator with couplings α e ∼ few × 10 −6 [few×10 −5 ] and α X ∼ 0.1 can evade all terrestrial and astrophysical constraints (see e.g. ref. [75]). For sub-MeV DM and mediator couplings this large, vanilla DM models may be brought into thermal equilibrium with the photon plasma and affect N eff and/or BBN [78]. If the DM is a real scalar, even if the DM is brought into thermal equilibrium, it remains consistent with BBN constraints. These constraints can also be evaded in 'non-vanilla' models, for instance by varying the couplings or the masses of the system with temperature. A full exploration of such models will be detailed in a separate publication [79].

Models and results
Having discussed the generic bounds that are relevant for constraining the allowed size of the direct-detection cross section, we now move to describing several concrete models and their results. Scalar and vector mediators are addressed in section 5.1; a kinetically mixed hidden photon (including in-medium effects) is treated in section 5.2; milli-charged DM is considered in section 5.3; and dipole-interacting DM is detailed in section 5.4.

Scalar and vector mediator
We begin by considering a real scalar mediator φ, described by the potential The DM-electron t-channel scattering for a vector mediator, such as a U(1) B−L vector boson, with tree level couplings to electrons and DM, are the same as those obtained for the scalar case, in the non-relativistic limit. We consider mediators both lighter and heavier than the momentum transfer involved in the process.
Light scalar/vector mediator. The self-interaction constraints of eq. (4.5) can be combined together with the kinetic decoupling requirements of eq. (4.7) and the stellar bounds of eq. (4.9) to learn how large the scattering cross section of DM and electrons eq. (3.9) can be. In the top panel of figure 5 we plotσ light DD of eq. (3.11) for several benchmark points, labeled I-III, shown in solid colored curves. As is evident, large cross sections can be obtained even for very small couplings. This is due to the enhancement of the cross section at low momentum transfer when the mediator is light, as shown in figure 5. The presented benchmarks all obey self-interaction constraints and also ensure that DM is kinetically decoupled through the time of recombination for mediator masses m φ eV. However, the depicted benchmarks may bring the mediator into equilibrium with the SM plasma via Compton-like processes; this is not a problem for a real scalar mediator, but a light vector mediator can then contribute too many degrees of freedom to N eff [78]. Likewise, the stellar constraints of eq. (4.9) are very stringent. In vanilla models, we find that stellar cooling is too severe to allow for a detectable rate for a light real scalar or vector mediator, which constrains g e itself to be below the 10 −14 level. We learn that unless stellar bounds are somehow lifted, the direct detection experiments considered in this paper will be unable to probe scalar or vector mediators with masses below ∼ O(10) keV. Stellar and BBN JHEP08(2016)057 Figure 5. Top: direct detection cross section, eq. (3.11), for light DM scattering off electrons via a scalar or (non kinetically mixed) vector mediator, for several benchmarks. These are I: α X = 10 −15 , α e = 10 −12 ; II: α X = α e = 10 −15 ; and III: α X = 10 −15 , α e = 10 −18 . These depicted parameters obey bounds from self-interactions and decoupling at recombination for m φ eV, though stellar emission (and BBN considerations for vectors) may place strong constraints; see text for details. Bottom: direct detection cross section between light DM and electrons, for several benchmarks of heavy mediators. These are A: m φ = 1 MeV, g e = 10 −5 e, α X = 0.1; B: m φ = 10 MeV, g e = 10 −5 e, α X = 0.1; and C: m φ = 100 MeV, g e = 10 −4 e, α X = 0.1. These depicted parameters obey all terrestrial and stellar-cooling constraints, though sub-MeV DM interacting with SM through a massive mediator may be strongly constrained by BBN; see text for details. The Xenon10 electron-ionization data bounds [80] are plotted in thin dashed gray. In both panels, the black solid (dashed) curve depicts the sensitivity reach of the proposed superconducting aluminum devices, for a detector sensitivity to recoil energies between 1 meV-1 eV (10 meV-10 eV), with a kg·year of exposure. We have included only the solar neutrino background in our estimate. For comparison, the gray dot-dashed curve depicts the expected sensitivity utilizing electron ionization in a germanium target as obtained in ref. [22].

JHEP08(2016)057
constraints may be lifted, however, for instance via having the coupling g e vary with environment, or via trapping effects, in which case sizable direct detection cross sections can be accommodated [79]. A kinetically mixed vector mediator can also lift stellar constraints, as will be discussed in section 5.2.2.
Heavy scalar/vector mediator. Moving to a massive scalar or vector, we focus on m φ few MeV. In the bottom panel of figure 5 we plotσ heavy DD of eq. (3.11) for several benchmark points, labeled A-C, which survive all terrestrial and stellar cooling constraints, as outlined in sections 4.3 and 4.4. Fairly large couplings to electrons are possible despite supernova constraints due to stellar trapping effects, and beam dump constraints can be evaded by decaying invisibly to additional dark-sector particles. As mentioned earlier, for values of α e and α X as large as these benchmark points, DM and/or the mediator can be brought into thermal equilibrium with the SM plasma. BBN and Planck limits on the effective number of relativistic degrees of freedom in equilibrium N eff [78] thus place constraints on, at least, the simplest of such models. If the mediator is heavy enough that it does not contribute to N eff at BBN, if DM is a real scalar, then N eff constraints are trivially satisfied. Even if a real scalar mediator, along with a real scalar DM, is in thermal equilibrium at BBN, N eff constraints (which allow one additional fermion, or two additional real scalars) are still satisfied at 95% C.L. For DM with more degrees of freedom, these bounds can potentially be lifted; this is the case, for instance, if the couplings and/or the masses of the particles involved evolve during the thermal history of the universe. (This is much in the spirit of refs. [81][82][83].) Model building efforts along these lines, both for relaxing BBN/N eff constraints for light or heavy mediators as well as lifting stellar constraints for light mediators, are being pursued in detail elsewhere [79].
Reach. The 95% expected sensitivity reach for a kg·year of our proposed superconducting aluminum experiment for light and heavy mediators is depicted in the thick black curves of both panels of figure 5, with the dashed [solid] curves showing the expected sensitivity with a 10 meV [1 meV] operating threshold. Given that the heat sensors on the detector are likely to have a somewhat limited dynamic range, we also place an upper bound of 10 eV [1 eV] on the detectable energy. The depicted curves then correspond to 8.8 [3.7] events per kg·year, taking into account the expected 3 [< 1] neutrino background events as found in section 3. For completeness, we show the Xenon10 electron-ionization bounds [80] in the thin gray dashed curves (these are absent in the top panel, as they are orders of magnitude weaker than the displayed parameter space). We also show the projected reach curves utilizing electron ionization techniques in a semi-conductor germanium target (silicon performs similarly) as obtained in ref. [22], translated toσ DD of eq. (3.11), shown in the thick gray dot-dashed curves. For massive mediators and DM heavier than a few hundred keV, the projected reach from a germanium/silicon target is comparable to than our proposed detection method, while for lower DM masses, where electron ionization techniques lose sensitivity, superconducting devices win. Light mediators further demonstrate the strength of our proposed detectors. When the mediator is light, superconductors can out-perform electron-ionization techniques by several orders of magnitude for dark masses above several JHEP08(2016)057 hundred keV. Going to even lighter masses, superconducting detectors are uniquely staged to probe such super light DM.

Kinetically mixed U(1) D
Next we study DM scattering with electrons through a kinetically mixed dark U(1) D . We consider a hidden photon mediator A which is kinetically mixed with the ordinary electromagnetic photon, with F µν (F µν ) the (hidden) photon field strength and the kinetic mixing parameter. J µ EM and J µ DM are the dark and electromagnetic currents, respectively. Diagonalizing the kinetic terms and moving to the mass basis, the hidden photon couples to the electromagnetic current of the SM with strength g e = e . The mass of the hidden photon m A is obtained via a dark Higgs mechanism or Stuckelberg mechanism. We note that we do not include the terms involving the dark Higgs in the above Lagrangian. These will be relevant only when we consider stellar cooling processes; they do not affect the DM-electron scattering of the direct detection process. Indeed, the direct-detection scattering rate is the same regardless of the mass mechanism for the hidden U(1) D , be it via a dark Higgs or the Stuckelberg case.

Photon propagator in medium
In any model that can be written in the form of eq. (5.2), the size of the effective kinetic mixing parameter, eff , is medium-dependent: where here q = (ω, q) is the four-momentum transfer of a process and Π T,L is the in-medium polarization tensor, defined according to with T,L the transverse and longitudinal polarization vectors: Now the question is how to extract Π L,T . In the appendix, we show using Maxwell's equations that, for a non-magnetic medium, JHEP08(2016)057 whereñ = n − ik is a complex index of refraction that is related to the conductivity σ and electric permittivity ε r via [84] 1 − ε r = 1 −ñ 2 = − iσ ω . (5.8) The conductivity of a target metal, such as superconducting aluminum, differs tremendously from that of an insulating target, such as Helium. For Helium, the index of refraction is very close to unity, implying that the photon mass in Helium is negligibly small. In contrast, in a metal like aluminum, the electric permittivity can be quite large and the photon mass appreciable. The relative permittivity as a general function of ω and q is given by [85]: where λ 2 TF = 3e 2 n e /(2E F ) is the Thomas-Fermi screening length. For aluminum, λ TF 4 keV, taking E F = 11.7 eV. In figure 6 we show the real and imaginary parts of √ Π as a function of |q| with fixed values of ω. Note that the imaginary part is only non-zero in a limited range of |q| and ω where the kinematics allows a photon to be absorbed; this corresponds to ω between ω = 1 2me (2|q|p F + q 2 ) and ω = 1 2me (2|q|p F − q 2 ) (or with an overall minus sign, depending on the choice of ω). As is evident, for typical q ∼ v F ω the effective photon mass in medium is approximately ∼ keV, on the order of the Thomas-Fermi screening length, implying that for typical momentum transfers of order 10 eV for m X = 1 keV, the direct detection rate in metals is severely limited. Note that it is the Thomas-Fermi screening length (of order a few keV), and not the plasma mass (typically O(10 eV)), which is the relevant screening parameter for scattering processes, where q ω; the plasma mass becomes the relevant screening mass for processes where ω q. We learn that for a kinetically mixed hidden photon mediator, an insulating target is preferred.
We can now incorporate these in-medium effects and compute the scattering crosssection for DM off of a nucleus or an electron via the exchange of a dark U(1) D . The matrix element is given by where Here G IM,βν is the in-medium photon propagator, which can be parameterized in Lorentz gauge as [86] G IM,µν (q) = P L,µν where the projection operators are Utilizing the Ward identity, one finds that the second term of the first factor in eq. (5.11) vanishes. Further, since we are only interested in non-relativistic scattering between the DM and electron, the zeroth components of the external momenta are much larger than the spatial components. In the non-relativistic limit, we find that the leading contribution comes from the longitudinal component, with the transverse components suffering velocity suppression. Thus in the following calculation, we keep only the longitudinal component of the photon propagator, and use . (5.14) where we use the relation Π L = q 2 |q| 2 Π 00 . Plugging this back to eq. (5.11), and simplifying using the Ward identity, we find .
Combining eqs. (5.10) and (5.15) we obtain where in the above we used the non-relativistic approximation, with q = (ω, q) = (p 1 − p 3 ). Utilizing eqs. (3.1), (3.6) and (5.16), we can now compute the rate for an aluminum target. The differential rate per kg·year as a function of deposited energy is given in figure 7, for several benchmark points. (Note that for heavy mediator and very light DM, the rate is always substantially smaller than the depicted range; as a result we do not show a corresponding benchmark point.) Comparing to figure 4, we find as expected that the in-medium effects essentially modify the qualitative behavior of the light mediator into that of a massive one. The resulting projected reach of such a metal target will be reduced accordingly, increasing the desirability to find an insulating target with small gap.

Kinetically mixed stellar constraints
An upper limit on the size of the direct detection cross section arises due to constraints on the relevant couplings: DM self-interactions constrain g X ; stellar cooling bounds the size of g e ; and requiring kinetic decoupling of the DM and mediator from the SM plasma such that CMB measurements are obeyed constrains the combination of the two couplings. The selfinteractions and kinetic decoupling constraints presented in section 4 directly apply here. Stellar emission constraints for a light kinetically mixed hidden photon differ, however, from those presented above, and are largely lifted, as we now discuss.
For the hidden photon masses in our range of interest m A eV, stellar constraints have been worked out in detail in the literature [87,88]. The dark photons are emitted from the sun and horizontal branch stars either through plasmon resonance conversion, or in association with photon decays, via a Higgstrahlung process. The former process proceeds regardless of the origin of the hidden vector's mass, while the latter exists only in the case of a dark Higgs mechanism. The crucial difference between these two cases arises in the JHEP08(2016)057 small-m A limit: the Higgstrahlung process does not vanish with vanishing mass, while the plasmon resonance conversion scales with ∝ m 2 A and vanishes for a massless mediator [88]; for further details, see refs. [87,88]. In the small mass region, where m A ω p with ω p ∼ 100 keV the plasma frequency in the sun and horizontal branch stars, the (direct) emission power of dark photons per volume is governed by the emission of longitudinal modes of A , and is proportional to ∝ m 2 A ω 3 p α e . The rate for the Higgstrahlung process is governed in the small mass region by decays of transverse photons, as they are more abundant than the longitudinal plasmons. The total energy power density of dark radiation is then proportional to ∝ ω 5 p α e α X q 2 H D , with q H D denoting the dark Higgs charge under the U(1) D (relative to the DM charge). The resulting stellar constraints are found to be [88]: Resonance conversion : for mediator masses 10 −5 eV m φ eV that we consider. For even lighter mediator masses, is bound by photon-dark photon mixing through level-crossing in the CMB [91], as well as from the CROWS experiment [89,90] and measurements of deviations from Coulomb's law [91]; these are lifted for m φ 10 −14 eV, where measurements of the shape of the static magnetic field of Jupiter allows kinetic mixing as large as O(10 −2 − 1) (see e.g. refs. [90,91] for a summary of constraints).
Combined with the self-interaction constraints on α X and stellar emission constraints on the DM as well, one can identify as a function of m X and m A the strongest constraints and place a bound on the combination α e α X which enters the direct detection cross section.
In the case of a dark Higgs mechanism, assuming similar dark-charges of the dark Higgs and the DM, we find that for m A below ∼ 10 −5 eV, the cooling is dominated by the Higgstrahlung process in the entire DM mass range of interest, despite the strong suppression of g X from self-interactions. For 0.1 meV m A eV, plasmon resonance conversion dominates the cooling for light DM masses until the Higgstrahlung process takes over; for m A = meV the turnover point is m X ∼ 200 keV, and increases with increased m A . For m A eV, combined with the self-interaction bounds on g X , plasmon resonance conversion is most important for the entire m X range of interest. If a hierarchy between the dark Higgs and DM charges is present, the Higgstrahlung constraints can be relaxed accordingly. Likewise, in the Stuckelberg case, only plasmon resonance conversion is relevant, and as the hidden photon mass decreases, the stellar bounds on A are weakened, at which point the other bounds mentioned above play a role.
We note that the analysis of refs. [87,88] does not include effects of trapping and absorption in the relevant stellar objects, which can, in principle, open up parameter space above the constraints presented there. Taking into account the low density in these stellar objects compared to that of supernovae, however, we expect that trapping becomes important only for very large kinetic mixing values where other (terrestrial) observations already exclude the parameter space. Higgs mechanism for light dark U 1 D Ge Figure 8. Upper bounds on the direct detection cross section of eq. (5.18) for light DM scattering off electrons, for a light kinetically mixed hidden photon mediator obtaining its mass via a dark Higgs mechanism, for a variety of different mediator masses (solid colored curves). The expected reach of a superconducting aluminum target with sensitivity to energies between 10 meV-10 eV and 1 meV-1 eV, as well as a germanium target [22], is shown in thick dashed black, solid black and dot-dashed gray, respectively. We have included only the solar neutrino background in our estimate. The in-medium effects of a metal target do not enable DM detection via a superconducting metal in this case due to strong stellar constraints on the relevant coupling; detection of these models would require an insulating target.

Kinetically mixed results
The direct detection cross section between electrons and DM through the exchange of kinetically mixed hidden photon can now be constrained. We take into account self-interactions via eq. where in the above we have taken the photon plasma mass Π L at a typical value of ∼ keV. We consider separately the light and heavy mediator regimes.
Light mediator. The largest allowed direct-detection reference cross sectionσ DD for the Higgs case with q H D ∼ 1 is depicted in the solid colored curves of figure 8, for a variety of light mediator masses m A eV. The kink in the curves as the mass of the DM increases is due to the change in the stellar constraints as the dominant cooling mechanism evolves (factoring in self-interaction constraints on α X ) from plasmon resonance conversion emitting the A to the Higgstrahlung process, as detailed above.
If the charge of the dark Higgs is substantially smaller than that of the DM, or if the hidden photon obtains its mass through the Stuckelberg mechanism, stellar constraints on the mediator are lifted as m A → 0 as discussed below eq.  Figure 9. Top: upper bounds on the direct detection cross section, eq. (5.18), for light DM scattering off electrons via a kinetically mixed hidden photon, which obtains its mass via the Stuckelberg mechanism, for a variety of different mediator masses (solid colored curves). Constraints include stellar cooling [87], CMB [91], CROWS [89,90], measurements of Coulomb's law [91], decoupling at recombination [67,68] and self-interactions [59]. Bottom: direct detection cross section between light DM and electrons, for several benchmarks of heavy mediators (same as in figure 5). These are A: m φ = 1 MeV, g e = 10 −5 e, α X = 0.1; B: m φ = 10 MeV, g e = 10 −5 e, α X = 0.1; and C: m φ = 100 MeV, g e = 10 −4 e, α X = 0.1. These depicted parameters obey all terrestrial and astrophysical constraints, though sub-MeV DM interacting with SM through a massive mediator may be strongly constrained by BBN; see text for details. In both panels, the Xenon10 electron-ionization data bounds [80] are shown in thin dashed gray. The black solid (dashed) curve depicts the sensitivity reach of the proposed superconducting aluminum devices, for a detector sensitivity to recoil energies between 1 meV-1 eV (10 meV-10 eV), with a kg·year of exposure. We have included only the solar neutrino background in our estimate. For comparison, the gray dot-dashed curve depicts the expected sensitivity utilizing electron ionization in a germanium target as obtained in ref. [22].

JHEP08(2016)057
strongest amongst all constraints, we plot in the top panel of figure 9 the upper bound on σ DD in this case for several sample mediator masses, shown in the solid colored curves. (We note that for very light hidden photon mediators, stellar emission of DM beneath 100 keV severely suppresses the allowed cross section (see figure 10 below); for this reason, only m DM 100 keV is shown here. In this region, constraints from SN emission of the DM can be released via trapping effects, and so do not control the largest allowed cross section, which we show.) As is evident, in contrast to the dark Higgs case, large direct detection cross sections are possible for the Stuckelberg case.
Heavy mediator. In the bottom panel of figure 9 we plot several benchmark points, labeled A-C, shown in solid colored curves. These theory benchmark curves are the same as those presented in the bottom panel of figure 5, modified toσ DD here. As was the case for the massive scalar/non-kinetically mixed vector mediator, if DM is lighter than approximately 100 keV, then it must either be a real scalar or be thermally unpopulated at BBN in order to satisfy constraints. Alternatively, its couplings and/or mass can vary within the thermal history of the universe. (Similar statements hold for the very light mediator as well, shown in the cyan curve in the top panel of figure 9, for very light dark matter masses.) Reach. The 95% expected sensitivity reach for a kg·year of our proposed superconducting aluminum experiment is shown in the thick black curves in figures 8 and 9. The dashed [solid] curves show the sensitivity when operating with a 10 meV to 10 eV [1 meV to 1 eV] dynamical range. The reach of the superconducting devices for both light and heavy mediators differs from that in figure 5 because the plasma effects are important. (Note that the effectively massive behavior of the mediator in-medium results in a better reach for the 10 meV-10 eV dynamical range compared to the 1 meV-1 eV range; this is because the rate is now peaked at higher energy deposits.) As is evident from figures 8 and 9, while the superconducting metal target is not appropriate for detecting some classes of kinetically mixed light hidden photon models, it is capable of probing others. When the kinetically mixed photon obtains its mass via a dark Higgs mechanism, superconductors are not ideal DM detectors unless the dark Higgs charge is substantially suppressed compared to that of the DM. In contrast, a kinetically mixed photon with Stuckelberg mass could allow for DM detection via superconductors. The in-medium effects of the metal hurt the low-DM mass reach due to the large plasma mass of the photon, and the reach of a semi-conductor target such as germanium or silicon is comparable to the superconducting devices for masses above a few hundred keV. Below that, where semi-conductors lose sensitivity, our detectors could be sensitive to DM masses above 100 keV that scatter by a kinetically mixed U(1) D with very small Stuckelberg mass.

Milli-charged dark matter
We now analyze the reach of our method into the parameter space of milli-charged DM particles X with electromagnetic charge Q. The 'mediator' between the DM and the visible sector is simply the photon, where the DM couples to the electromagnetic current with strength Qe. In our notation, this means g e = e and g X = Qe.

JHEP08(2016)057
The potential reach of the superconductoing devices we propose can easily be translated into the Q − m X plane. Constraints on the milli-charge Q as a function of the DM mass m X have been worked out extensively in the literature. Stellar cooling from red giants (RG), white dwarfs (WD) and supernovae (SN) as well as big bang nucleosynthasis (BBN) are worked out in ref. [92]. In addition, ref. [68] considers the requirement of DM be decoupling from the plasma at the time of recombination. The possibility of charged DM being evacuated from the disk was also considered in ref. [68], though the argument leading to the constraint is not bullet-proof. In the mass range of interest, Xenon10 constraints exist as well [20,80]. These existing constraints are depicted in figure 10, along with the potential reach of our proposed method. For completeness, we show the projected reach using a semi-conductor germanium target as well [22]; silicon performs similarly. For masses above a few hundred keV, a germanium/silicon experiment can outperform superconductors, while for lower masses, where semi-conducting targets loose sensitivity, the large in-medium effects of a photon in a metal suppress the reach of superconductors into the milli-charged DM parameter space. A viable region can be probed, though the region can be broadened if stellar and/or BBN constraints are lifted. For milli-charged DM, we learn that an insulating target would perform better.

Dipole operator dark matter
It is possible that DM particles directly couple to photons through a magnetic or electric dipole operator, which for Dirac fermion DM X takes the form Both operators above induce similar physics effects in stellar cooling processes, and so can be described via an effective dipole moment, For m X O(keV), the induced DM emission processes from the Red Giant (RG) branch in globular clusters place a strong bound of [93] µ DM 3 × 10 −12 µ B [RG] , (5.21) where µ B = e 2me 300 GeV −1 is the Bohr magneton. For heavier DM, emission from White Dwarfs (WD) and supernova are relevant. WD cooling restricts [94] µ DM 5 × 10 −12 µ B [WD] , (5.22) which is comparable to the RG constraint, but applicable for a wider range of masses, m X O(MeV). For heavier DM with masses up to O(100 MeV), supernovae provide the best constraints, and the allowed range is [95] µ DM 2 × 10 −11 µ B or µ DM 2 × 10 −12 µ B [SN] . (5.23) Above, the upper limit on µ DM comes from emission considerations, while the lower limit arises from trapping effects, which kick in as the coupling between DM and ordinary matter JHEP08(2016)057 Figure 10. Constraints and reach for milli-charged DM: stellar emission bounds from red giants (red), white dwarfs (green) and supernovae (orange) [92]; BBN (blue) [92]; decoupling at recombination (gray) [68]; Xenon10 (gray thin dashed curve) [80]; evacuation from the disk (dashed magenta curve) [68]; projected reach using a germanium target (thick gray dot-dashed curve) [22]; expected reach with a superconducting aluminum device with a sensitivity to recoil energies between 10 meV-10 eV (black thick dashed curve) and 1 meV-1 eV (black thick solid curve). We have included only the solar neutrino background in our estimate.
increases and release the constraints. We note that while the analysis of ref. [95] does not include Pauli blocking effects, the above bound converts into an effective suppression scale of the dimension-five dipole operator of order ∼ 10 9 GeV, which is comparable to the constraint on the suppression scale of the axion-photon coupling from supernova cooling considerations [69]. We thus believe the order of magnitude of the constraint is valid.
The DM-electron scattering cross sections from eq. (5.19) scale as [96] dσ where the electric dipole scattering is enhanced by 1/v 2 X , and no low-velocity enhancement is present for scattering through a magnetic dipole. This is in contrast to the scalar/vector mediator cases and milli-charged DM, where the scattering cross section enjoys a lowvelocity enhancement of 1/v 4 X . Comparing the milli-charged and electric dipole cases, we JHEP08(2016)057 have roughly where Q is the milli-charge of the DM and µ eX is the reduced DM-electron mass as usual.
Taking v X ∼ 10 −3 , we have (5.26) Thus we see that for DM of 10 keV mass, an experiment which is sensitive to DM with milli-charge as small as 10 −17 can be useful in probing the unconstrained parameter space in the electric dipole operator DM scenario. The reach for magnetic dipole operator is worse since there is no velocity enhancement. Comparing to figure 10, we learn that the proposed superconducting detectors will not be sensitive to dipole DM.

Conclusions
We have explored in detail a proposal for detecting DM with Fermi-degenerate materials, focused on the case of a superconducting metal target. We computed the scattering rate of DM off of the electrons, factoring in the suppression due to Pauli blocking effects. We considered cosmological and astrophysical constraints from DM self-interactions, kinetic decoupling in the early universe as well as stellar emission, together with terrestrial constraints, such as (g − 2) of the electron and beam dump experiments. These constraints were then applied to a variety of models, such as a simplified model of a scalar or vector mediator, a kinetically mixed U(1) D and milli-charged DM. We have shown that viable regions of parameter space exist, consistent with various cosmological, astrophysical and terrestrial constraints, which are detectable with our proposed experiment. A broader range of model parameter space becomes available if stellar and/or N eff constraints on light degrees of freedom are lifted; we leave the exploration of such models for future work. We also computed in-medium effects for the kinetically mixed dark U(1) D , and found that the plasma mass of a photon in a metal substantially reduces the reach of superconducting detectors for this class of models. There are several further directions that we are pursuing. First, since the reach in a metal is reduced for a kinetically mixed dark U(1) D , other types of target materials should be examined which feature small or zero energy gap and simultaneously also small inmedium photon mass; graphene is one possibility. Second, in this paper our attention was restricted to t-channel scattering of DM off the target. We plan also to examine, however, the absorption of very low mass states via the excellent energy resolution of the experiment. Third, while we focused here on targets with a substantial initial state velocity, in order to find configurations where the entire kinetic energy of the DM can be extracted, we also note that O(meV) energy deposits on a light nucleus could also be utilized for probing light DM. We are currently pursuing helium targets as well.

JHEP08(2016)057
Over the last decades, the main focus of the DM community has been directed towards axions and the weak scale as the source of DM. As the pursuit in the search for DM expands, it is important to consider as broadly as possible what types of DM experiments can be built and what types of models, consistent with all known constraints, could be detectable. Exploiting superconducting targets is an important step along this path.

(A.3)
We now use Maxwell's equations to write this in terms of the conductivity. We know that longitudinal conductivity is defined by J = σ L E. We also know that current conservation dictates ∂ µ J µ = 0 implying ωJ 0 = q i J i (we will use an Einstein summation convention and roman letters to denote spatial indices in this appendix). We can thus write, using eq. (A.1) and the fact that R 0i A i = 0 in Coulomb gauge, We also have from Maxwell's equations in Coulomb gauge ∇ · E = −∇ 2 φ. Identifying φ = A 0 , this allows us to write

JHEP08(2016)057
We use Maxwell's equation E = − ∂ A ∂t − ∇φ = iω A − i qA 0 . This allows us to write Then from the first Maxwell equation ∇ · E = −∇ 2 φ we can further write Using 2σ T = P T ij σ ij , we learn that the transverse conductivity is Identifying Π T = [δ j − q j q q 2 ]R j , we recover the transverse part of eq. (5.7).
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.