position and efficiently discriminates between nuclear recoil (NR) and electron recoil (ER) signatures,
etc. The NR backgrounds arise from neutrons, neutrino-nucleus interactions, etc., whereas the ER
background arise and from β-decays, γproduced by radioactivity, neutrino-electron interactions, etc.
Experiments like Xenon, LZ, and PandaX are exploring possible DM events in the presence of these
backgrounds.
Recently, LZ collaboration has published its first result [23]. The experiment is situated at 4850
ft underground in the Davis Cavern at the Sanford Underground Research Facility (SURF) in Lead,
South Dakota, USA. The total mass of LXe is 10 t, out of which only the inner fiducial 5.5 t is used
for DM searches to reduce the backgrounds. With 60 live days of data, LZ has reached the current
strongest constraint 6 ×10−48cm2at DM mass 30 GeV. Compared to the previous strongest bound,
this is 6.7 and 1.7 times better at DM mass ∼30 GeV and ∼1000 GeV, respectively.
While the LZ result focuses on the searches for WIMP-like DM, in this paper, we take a complemen-
tary approach to investigate scenarios of sub-GeV DM interacting with nucleons via spin-independent
(SI) interactions. Non-relativistic sub-GeV DM, typically moving with velocity ∼10−3, will not be
able to impart enough energy to produce an observable NR in the LZ experiment. However, an en-
ergetic sub-GeV DM may produce sufficiently large NR. One of the simplest ways to produce such
boosted DM is to consider the interaction between high-energy cosmic-rays (CRs) and DM, known as
CR boosted DM (CRDM), proposed for the first time in Ref. [24] for nuclear scattering and Ref. [25]
for electron scattering. Further this technique has received considerable attention [26–66]. These
boosted DM particles reach the underground detector with much higher energy which helps to over-
come the energy threshold although with much lower flux. Even with this lowered flux, it is possible
to probe new regions of DM-nucleon scattering cross-section, since the bounds for sub-GeV DM using
other techniques are weak. The paradigm of CRDM premises only on the assumption of DM-nuclear
interactions, which is also true for many DD experiments. A large class of particle physics models
predicts such interaction for sub-GeV DM [67–72].
Knowledge of the CR spectrum is an important ingredient in computing CRDM flux. The direct
CR flux measurements (PAMELA [73], AMS-02 [74,75], CREAM-I [76], etc.) are done with balloons
and satellite detectors near the top or outside the atmosphere. This has been used as input CR flux
in Ref. [24]. However above 100 TeV CR fluxes are small hence direct measurements are not a
feasible choice. In this case, CR is measured indirectly through the air shower induced by it. We
utilize the parametric fit of CR flux measurement (obtained by combining direct and indirect CR
flux measurements) given in Ref. [77] as the input CR flux. Then we explore the signature of the
CR-induced DM in the LZ experiment. We find a factor ∼2 improvement compared to previous limit
of XENON1T near DM mass ∼1 MeV. We also present the projections of the upcoming XENONnT,
LZ, and Darwin in probing the DM-nucleon cross-section for sub-GeV DM. We find that there can be
a factor ∼10 improvement for Darwin compared to the current LZ limits.
The paper is organized as follows. In Sec. 2, we briefly sketch the CRDM framework. In Sec. 3, we
present limits from LZ and future xenon-based experiments. We conclude in Sec. 4.
2