
DM incoming from the Galactic halo, which is usually as-
sumed for direct detection experiments. The other is the
thermalized DM component. This thermalized compo-
nent exists as once DM enters the Earth, it can thermal-
ize, and become captured and bound to the Earth. For
sufficiently large DM-SM scattering cross sections (larger
than about 10−35 cm2), the DM rapidly thermalizes and
is said to be in local thermal equilibrium with the sur-
rounding SM matter. In this case, the DM radial profile
within the Earth, nχ, is dominantly governed by the dif-
ferential equation [36]
∇nχ
nχ
+ (κ+ 1) ∇T
T+mχg
T=Φ
nχDχN
R2
⊕
r2,(1)
where Tis the Earth’s radial temperature profile at po-
sition r,R⊕is Earth’s radius, mχis the DM mass, gis
gravitational acceleration, and Φ is the incoming flux of
DM particles from the Galactic halo. DχN ∼λvth and
κ∼ −1/[2(1+mχ/mSM)3/2] are diffusion coefficients [36],
with λthe DM mean free path, vth the DM thermal ve-
locity, and mSM the SM target mass. The DM density
profile is normalized by enforcing that its volume integral
equals the total number of particles expected within the
Earth [36].
Solving Eq. (1) for nχ(r) reveals that this thermalized
population of DM can be significantly more abundant at
the Earth’s surface than the incoming halo DM particles.
For DM masses around a GeV, the local DM density can
be as high as about ∼1014 cm−3. However, as this popu-
lation is thermalized within the Earth, its velocity is low.
We approximate the thermalized DM velocity distribu-
tion as a truncated Maxwell-Boltzmann distribution,
fχ(v) = 1
N0
e−(v/vth )2Θ(vesc −v),(2)
where N0normalizes the distribution, and v2
th =
8Tχ/πmχwith Tχ≃300 K. This velocity would require
thresholds of E<
∼0.05 eV for conventional detection
techniques. This is much lower than the reach of typ-
ical direct detection experiments, and so requires new
techniques to be detected. Our assumption of DM being
at room temperature of ∼300 K is reasonable, as even at
the largest cross sections considered the mean free path
is much larger than the size of our devices, such that DM
is not expected to thermalize with the device itself.
For halo DM, in Eq. (2)vth is replaced by the
average DM velocity in the halo v0= 230 km/s. In
this case, the relative velocity between the Earth and
DM also becomes important. Hence, for halo DM we
use the boosted velocity v→v+v⊕in Eq. (2), where
|v⊕|= 240 km/s is the Earth’s velocity in the galactic
rest frame. The halo DM density is assumed to be
0.4 GeV cm−3. We now show that quantum devices
are highly sensitive to DM with low energy depositions
through their power measurements, which includes both
the thermalized DM population, as well as light halo DM.
Scattering Rate & Energy Deposition.— As a DM
particle with velocity vscatters in the detector and trans-
fers momentum q, it deposits an amount of energy
ωq=q·v−q2
2mχ
=Ef−Ei.(3)
As a result, the target makes a transition from |i⟩to |f⟩.
For such low energy depositions, the momentum trans-
ferred is comparable to the inverse size of nuclear wave-
function in a detector crystal, and the inter-atomic forces
become important. Hence, lattice vibrations or phonon
excitations will be used to compute the DM scattering
rate. The total rate per unit target mass can be written
as [37,38]
Γ = πσχN nχ
ρTµ2Zd3vfχ(v)Zd3q
(2π)3F2
med(q)S(q, ωq) (4)
Here, fχ(v) is DM velocity distribution, ρTis the tar-
get density, σχN is the DM-nucleon scattering cross sec-
tion, µis the reduced mass of the DM-nucleon system,
Fmed(q) is a form-factor that depends on the mediator
(we assume Fmed(q) = 1), and S(q, ωq) is the dynamic
structure factor containing the detector response to DM
scattering and depends on the crystal structure of the
target material.
To compute DM scattering rates, we follow
Refs. [39,40] and use the publicly available code
DarkELF. We modify DarkELF in two main ways.
Firstly, we update the local DM density and DM veloc-
ity input to be that described in the previous section, for
halo or thermalized DM as appropriate. Secondly, the
code was developed only for materials with two atoms
per primitive cell, which is the smallest crystal unit.
Thus, we adapt it for materials like Al which has only
one atom in its primitive cell.
Detection Mechanisms and Materials.— Detect-
ing light halo DM or the captured DM population of low
thermal energy demands use of low threshold quantum
sensors that can detect ∼ O(10) meV energy deposition.
Such sensors are usually designed using superconducting
materials, which have small energy gaps [41–44]. Alu-
minum (Al) is a widely used superconductor for such a
purpose and its characterization data is readily available.
Such a small amount of energy transfer is not sufficient
for nuclear recoil or electronic ionization, however DM
can excite collective modes, such as phonons in the ma-
terial, resulting in an excess power. For example, in one
experimental setup, a bias circuit stabilizes the absorber
material at its transition temperature Tc, where its re-
sistance is very sensitive to any energy deposition. The
total power deposited in the detector by DM in the form
of phonons is
PDM =ϵZdω ω dΓ
dω ,(5)
where ϵis an efficiency factor that depends on the ex-
perimental setup. We will use this to calculate excess
2