Constraints on cosmic-ray boosted dark matter in CDEX-10

Dark matter (DM) direct detection experiments have been setting strong limits on the DM–nucleon scattering cross section at the DM mass above a few GeV, but leave large parameter spaces unexplored in the low mass region. DM is likely to be scattered and boosted by relativistic cosmic rays in the expanding universe if it can generate nuclear recoils in direct detection experiments to offer observable signals. Since low energy threshold detectors using Germanium have provided good constraints on ordinary halo GeV-scale DM, it is necessary to re-analyze 102.8 kg day data in the CDEX-10 experiment assuming that DM is boosted by cosmic rays. For the DM mass range 1 keV 1 MeV and the effective distance within 1 kpc, we reach an almost flat floor limit at cm2 for the spin-independent DM–nucleon scattering cross section, at a 90% confidence level. The CDEX-10 result is able to close the gap unambiguously in the parameter space between the MiniBooNE and XENON1T constraints, which were partially hindered by the Earth attenuation effect. We also quantitatively calculate the expected neutrino floor on searching for CRBDM in future direct detection experiments using Germanium.


I. INTRODUCTION
Compelling astrophysical and cosmological evidence points to the existence of DM [1,2], yet very little is known about particle nature of DM.If DM particles could annihilate or decay into standard model particles, observations of gamma rays, anti-protons and positrons in cosmic ray telescopes and neutrinos in the related detectors will provide clues for DM properties and even its origin indirectly [3][4][5][6].In addition, several decades have been spent on searching for DM directly with low-background detectors in the deep underground laboratories on Earth in order to turn the DM-nucleus or DM-electron recoiled energy into light, heat and ionization signals.DM direct detection experiments intend to distinguish nuclear recoils (NRs) and electronic recoils (ERs) with the noble gas liquid detectors such as PandaX [7], Xenon1T [8] and LUX [9] using liquid Xenon, DEAP-3600 [10] and DarkSide [11] using liquid Argon, while a few of complementary experiments with crystals register DM energy deposits only, such as CoGeNT [12] and CDEX [13] using Germanium.Rapid progress has been made, but we have to be confronted with null results more often than not.As a result, stringent exclusion limits have be placed in the plane of the DM mass m χ and the DM coupling strength σ χ0 , given the weakly interacting massive particle (WIMP) as a DM candidate in the standard halo model (SHM) [14,15].
In the high mass DM region, the neutrino floor caused by coherent elastic neutrino-nucleus scatterings (CEνNS) is around the corner and leads to intrinsic neutrino backgrounds for the next-generation large scale DM detectors [16][17][18][19][20][21][22].Suffering from the energy threshold in the current detector technology, a relatively large space remains unexplored in the light DM scheme.Tremendous efforts are still underway, aiming at ultra-low threshold and large-scale detectors to reach better sensitivities.Novel experimental detection strategies and a number of new theoretical ideas to probe the low mass DM have been proposed and actively discussed very recently [23], including better detection techniques to overcome the energy-recoil threshold barrier [24,25], Migdal effect to observe electrons dissociated from the atom through the nuclear scattering [26][27][28][29][30], the boosted DM from heavy-generation annihilations in a multi-component DM model in order to surpass the detector threshold limitation [31].One of the interesting topics to explore low-mass DM is the cosmic ray boosted DM (CRBDM) mechanism proposed in recent studies [32][33][34][35][36][37][38].In the galaxy, there must be an isotropic and relativistic cosmic-ray stream [39,40] consisting of protons (∼86%), helium nuclei (∼11%), some other heavier nuclei and electrons (∼3%).DM can generate nuclear recoils in direct detection experiments to offer observable signals.It is likely that DM is scattered and boosted by relativistic cosmic rays in the expanding universe.Several studies have already investigated the CRBDM and offered relatively tight constraints on light DM, including the diurnal effect [41], DM-nucleon interactions propagated by various mediators [42][43][44][45] and the so-called reverse direct detection method [46].There are cosmological constraints from Big Bang Nucleisynthesis (BBN) and the DM relic density.Previous studies show that CRBDM with large enough cross section is expected to achieve thermal equilibrium with SM particles before the BBN epoch and would affect the prediction of BBN.Consequently, current observation of BBN excludes the DM mass below a few MeV, which still leaves a viable parameter region to be probed [47].In addition, the DM relic density puts constraints on the low mass DM, but it largely depends on model assumptions [48,49].
As we noticed in the previous CRBDM studies, the attenuation effect results in the upper bound of the exclusion region for CRBDM.Both MiniBooNE and XENON1T provide lower bounds on CRBDM, despite the dependence of DM and cosmic ray distribution in the galaxy.The attenuation effect in XENON1T cuts in the exclusion region of CRBDM.Therefore, there exists a gap in the parameter space between MiniBooNE and XENON1T constraints [33].Compared with other types of detector, Germanium detectors with the lower energy threshold and better energy resolution are effective means to set limits on the scattering cross section for the GeV-scale DM in the standard WIMP detection.As a consequence, the CDEX experiment located at 2.4 km underground using the point contact high purity germanium detector might be able to close the gap due to its lower level of backgrounds and lower energy threshold compared to MiniBooNE and XENON1T.Based on the CRBDM mechanism, it must be valuable to reanalyze the 102.8 kg×day data of CDEX-10 experiment, and provide constraints on spin-independent DM-nucleon cross section for m χ < 1 GeV.This article is organized as follows.In Sec.II, we review the mechanism of CRBDM and derive the theoretical energy recoil spectrum with the Earth attenuation effect taken into account.In Sec.III, we re-analyze the CDEX-10 data in the CRBDM hypothesis.Finally, we summarize and make conclusions in Sec.IV.

II. COSMIC RAY BOOSTED DARK MATTER (CRBDM)
In a CRBDM mechanism, cosmic rays and DM exchange their energy and momentum in the process of collisions, so that the boosted DM would become almost relativistic.Then, boosted DM particles traveling toward the detector on the Earth and scattering off the target nuclei, could generate enough recoil energy which would exceed the threshold energy.The schematic diagram to illustrate acceleration and detection processes is given in Fig. 1.To study a scattering process, we need to know the information of the incident particle, the target, and the collision process.The local density of DM near the solar system is ρ χ0 ∼ 0.4 GeV/cm 3 according to the cosmological observations [2].The velocity distribution of DM in the standard halo model (SHM) [14] approximates to the Maxwell-Boltzmann distribution whose the most probable velocity v 0 is about 220 km/s.In the following, we will discuss the CRs and review the collision process between CRs and DM.

A. Boosted DM flux from CRs
CRs might originate from some of the supernova remnant (SNR) in the Galactic Disk [50].Fermi pointed out that CRs from SNR bounce constantly between the shock surface and the magnetic field behind the wave surface.To explain the highly isotropic distributions for energetic charged particles, the concept of CR diffusion was proposed, in which the Galactic magnetic field plays an important role in the process [40].To calculate the differential CR flux, two main components proton and helium nucleus are considered and marked with subscript i in the following formula.With the data from PAMELA, the fitting formula of the Local Interstellar Spectra (LIS) of CRs are parameterized [51] as follows: where R means the rigidity of each nucleus.For a particle with a proton number Z, the relationship between the rigidity R and the kinetic energy of CR particle T i is: Ti R = Ze.As we all know, the differential CR flux obeys a power law with an index ∼ −2.7.As shown in Eqn.(1), a, b, c, d, e, f, g are the numerical coefficients [51,52] and we show the coefficients of proton and helium nucleus in Table I.Considering that CRs is isotropic in the interstellar, we integrate the volume to obtain the differential CR flux dΦi dTi : dΦi dTi = 4π dΦi dR dR dTi .Compared with the relativistic velocity of CRs, DM at a velocity of ∼ 10 −3 c in the SHM can be safely treated at rest.Considering the collision of DM and CRs, although the inelastic scattering occurs due to the high energy of CRs, it only contribute to a negligible portion in the DM fluxes boosted by CRs compared to elastic scattering, and the signals of γ rays and neutrinos from inelastic scattering are beyond our topic [53].Consequently, it is reasonable for us to consider only the elastic scattering process.With the assumption that the collision is elastic and isotropic, we relate the kinetic energy of DM at T χ to that of CRs at T i by energy and momentum conservation in the following: where θ is the scattering angle in the center of momentum frame, and m i is the mass for CR particle i.According to Eqn. (2), T χ tends to decrease and is proportional to m χ if m χ m i , while DM would gain more energy if m χ is close to m i .DM gets the maximimal kinetic energy T max χ (T i ) when the CR direction is completely reverted after such a collision (θ = 0).The minimum energy of the incident particle T min i (T χ ) is given in Eqn.(3): Note that the sign +/− applies to the case with When momentum transfer occurs in the scattering process, the internal structure of the nuclei would take effect.Therefore, the cross section varies with the exchange of momentum described by the form factor G i : , where σ 0 χi is the cross section for zero momentum transfer.In the collision process with large momentum transfer, just as the case that DM is boosted by CRs, we adopt the dipole form to describe the effect of nonpoint-like structure of nucleon.Note that for the small momentum transfer, the Helm model is preferred.The dipole form factor of proton and helium is given by [54]: where Λ i is inversely proportional to the charge radius, and we set Λ p = 770 MeV for proton and Λ He = 410 MeV [55] for helium nucleus.It is generally believed that the DM couplings to the proton and nucleon are the same, so the cross section can be written as Eqn.( 5) with a dependence of the mass number A for the incident particle i: It is worth mentioning that the factor A 2 in Eqn. ( 5) is 16 for the helium nucleus and 1 for the proton, respectively.The Helium nucleus and proton make contributions at the same order of magnitude to the differential flux for DM (the helium nucleus account for 10% of CRs).Based on the above information, we start to calculate the differential flux for boosted DM.Inside a volume dV , the collision rate with energy exchange dT i for CRs and dT χ for DM is given by: where dΦi dTi is the differential flux for CRs.Without a loss of generality, we make an energy-independent assumption here with such a form factor untouched that energy transfer is uniformly distributed as . By integrating the volume and the energy T i , we obtain the differential flux for the boosted DM in Eqn.(7).In addition, the volume integral asks for the distribution of DM and CRs in the galaxy.Note that the integral needs to be divided by the factor 4πd 2 to obtain the DM flux boosted by CRs, where d is the distance between the collision point and the Earth.We introduce the effective distance D eff to represent the volumn integral, into which we put the results of different assumptions on the CR spatial distribution.The CR spatial distribution is considered as the homogeneous and spherical distribution given in the reference [33] and results in D eff = 1 kpc conservatively.While other researchers assume that CRs are uniformly distributed in a cylindrical distribution with a radius R = 10 kpc and a half-height h = 1 kpc [34,35] and results in D eff = 3.7 kpc.Also, the CR distribution can be simulated with the GALPROP code [41].Therefore, it is worth further investigations on the CR distribution in the future work and we now take two benchmark examples of D eff for analyses:  2) with θ = 0. Following the uncertainty analyses on the DM local density and the DM spatial distribution [35], we take 50% uncertainty for analyses.To further study the boosted DM, we can make a comparison between the boosted DM flux induced by CRs and the DM flux under the SHM in the same reference frame: where f (v) is the Maxwellian velocity distribution with the most probable velocity at v 0 = 220 km/s and the galactic escape velocity at v escape = 540 km/s [56].Here the velocity distribution in SHM is taken as follows: where k is the normalization factor.The numerical results for the CRBDM fluxes are compared and shown in Fig. 2.
As the DM density in space is fixed, the smaller DM mass results in the larger number density of DM.Compared with the DM flux under the SHM hypothesis, the flux intensity boosted by CRs is a small portion of the total flux.Only a small portion of CRs lose energy, so that the CRBDM mechanism can survive in the current constraints from astrophysical experiments.It provides a clue for the reverse direct detection for DM [46].

B. Attenuation
As DM particles could collide with ordinary matter, the DM flux from space to the detector would be blocked by the atmosphere and the rock, which causes energy loss or de-acceleration for DM.We simply estimate the mean free path for DM as λ = 1/σ χN n N where σ χN stands for the DM-nucleus scattering cross section.Taking the mean mass number A = 20 in the Earth crust, the mean density ρ crust = 2.7 g/cm 3 in the crust [57] and σ χ0 = 10 −30 cm 2 as parameters, we get the the mean free path: λ ∼ 0.1 km.In order to suppress backgrounds from cosmic rays, direct detection experiments are usually located deep underground.For example, China Jinping Underground Laboratory to host the CDEX experiment has the rock overburden at z ∼ 2.4 km.Then it is time to establish a model to analyze the Earth attenuation effect and the ballistic trajectory approximation is adopted.
Under the assumption of energy independence: , we obtain the energy loss for DM crossing the crust: where T max r is the maximal recoil energy between the DM particle and the nucleus N .Here the form T max r is similar to Eqn. ( 2) with θ = 0 by a replacement of m χ → m N and m i → m χ .We have to solve the differential equation in Eqn.(10) to establish the relationship between DM energy in the space and detector: T z χ = T z χ (T 0 χ ).However, recent studies [58,59] involve the form factor in the energy loss calculations and make necessary modifications on Eqn.(10), which leads to a different constraint on the parameter space of DM.In the following section, we will make a comparison of these effects in the exclusion plots.After that, we will get the attenuated DM flux spectrum as follows: Adopting the nuclear composition of the Earth crust given in the reference [57] and the crust thickness as 2400 m, we present the DM flux with/without the Earth attenuation effect in Fig. 2. The comparison of different line colors tells us the interesting features associated with the DM mass while the dashed/solid lines highlight the cases with/without the attenuation effect, respectively.Under the same assumption for the cosmic ray acceleration mechanism, we see the low mass DM gets boosted into higher flux intensity region.However, the high mass DM tends to be bounced towards the low kinetic energy part by rock as the Eartch attenuation effect is taken into account.If such a bounce on Earth for CRBDM are too strong, the DM particle might not be able to reach the detector.The energy loss of DM particles crossing the earth crust depends on the maximal recoil energy, which also depends on the DM mass.As a result, lighter DM will lead to less recoil energy of particles in the Earth crust and suffer from less attenuation.It is noted that the total number of DM should be normalized in two scenarios related the Earth attenuation in order to verify the correctness in numerical calculations.Therefore, one cannot neglect the attenuation effect in sensitivity studies.As shown in the following section, we will see the ceiling limit caused by the Earth attenuation effect in re-analyzing experimental data.

C. Recoil spectrum
To obtain the recoil spectrum from CRBDM, we focus on the elastic and isotropic collision of DM and target particle in the detector.By neglecting energy dependence, we get the recoil spectrum by integrating the kinetic energy of DM particle dT χ : where T z χ and dTχ are the kinetic energy of DM and the DM flux in the detector with depth of z, respectively.Here the σ χN is given by: σ χN = σ 0 χN F 2 (E R ) where the F (E R ) is the form factor [15].Since the momentum transfer is small in the process of DM-nucleus scattering, we adopt the Helm model which describes the nucleon distribution of the target nucleus with two parameters: the effective nuclear radius and the surface thickness [60,61].We will follow the similar expression in Eqn.(3) and Eqn. ( 2) with θ = 0 to calculate T min χ and E max R .For scintillation and ionization detectors calibrated with γ sources, the observable nuclear recoil energy E v differs from the true recoil energy E R : , where f n is called the quenching factor as a function of E R calculated by the TRIM software package [62][63][64] and F (E R ) is more convenient to make the derivation.Besides, because of a finite detector energy resolution and electronic noise, recoils at energy E R would be observed as a Gaussian distribution.Consequently, we will obtain the modified recoil spectrum: where ∆E v is the energy resolution.As a demonstration, Fig. 3 shows the energy recoil spectrum from CDEX-10 data, and the predicted recoil spectrum for DM with m χ = 1 MeV and σ χ0 .Note that the quenching factor and the energy resolution are included in the calculation of theoretical predictions.The effect of the detector efficiency has been included in the experimental data so that we can directly compare the experimental data with the theoretical expectation from Eqn. (13) [65].

D. Neutrino floor
Future DM direct detection experiments will encounter an irreducible background from solar, atmospheric and diffusive supernovae neutrinos through the CEνNS process [16,18,22,66], which will hinder our searches on CRBDM.It is a question in statistics whether DM direct detection experiments will see neutrino events and how to evaluate the impact of neutrino backgrounds on the DM detection.The neutrino floor is defined as the discovery limit: when the true DM model lies above the limit, a given experiment has a 90% probability to obtain more than 3σ DM signals.
CEνNS [67] is the dominant channel of neutrino backgrounds due to the coherent enhancement.The differential cross section of the neutrino-nucleus scattering are given by: where W Z is the weak charge, θ W is the weak mixing angle, and F (E r ) is the Helm form factor.By analogy with Eqn.(12), one can obtain the CEνNS event rates in a DM detector.Then we construct a likelihood function L σ χ0 , φ based on the Monte Carlo data, where φ stands for nuisance parameters including uncertainties of different neutrino sources.In order to obtain the discovery limit, we take the background-only case as the null hypothesis H 0 and the alternative hypothesis H 1 includes both DM signals and backgrounds.Thus, the test statistic q 0 can be written as following: where ˆ φ is the value of φ which maximizes L for H 0 , while σχ0 and ˆ φ globally maximize L for H 1 .Because q 0 asymptotically follows a χ 2 distribution with one degree of freedom according to the Wilk's theorem, the 3σ statistical significance means the test statistics of an observation q obs larger or equal to 9. Finally, we carry out 500 times Monte Carlo experiments with the same m χ and σ χ0 and apply the statistical principle to obtain the neutrino floor.and MiniBooNE (blue) [33] in the framework of CRBDM at a 90% C.L. The gray region enclosed by solid and dashed black lines stand for constraint from CDEX-10 under the ballistic trajectory approximation with and without considering the form factor, respectively.As a comparison, we also include the exclusion limit from EDELWEISS [29] based on the Migdal effect and from CMB [68] (green) as well as the gas cloud cooling [69] (orange).The neutrino floor (brown) is obtained with an exposure of 1 ton•year and a recoil energy threshold of 1 keV for the Germanium detector.

III. CONSTRAINTS WITH CDEX-10 DATA
The CDEX-10 experiment located at China Jinping Underground Laboratory(CJPL) with about 2400 m of rock overburden makes use of p-type point contact Germanium (pPCGe) to reconstruct the ionization signal generated by the collision between DM and target particles [70].With the energy threshold as low as 160 eVee, CDEX-10 improves limits on σ χ0 down to m χ of 2 GeV.To obtain the constraints, we adopt the minimum-χ 2 analysis as the statistical method.Considering that high energy γ−rays from ambient radioactivity produce flat electron-recoil backgrounds at low energy [63,71], two free and positive parameters σ χ0 and b, which characterize the flat backgrounds and the spin-independent DM-nucleon scattering cross section σ χ0 , are used to construct the statistical measure such as: where S Ex k is the experimental data with the uncertainty σ k , and S T h k stands for the theoretical predictions in each bin, respectively.Note that we add a pull term accounting for the astrophysical uncertainty, where ρ is treated as a nuisance parameter with the central value ρ = 0.4 GeV/cm 3 and the uncertainty σ ρ = 0.2 GeV/cm 3 .To assess the exclusion limit, we calculate the probability in the likelihood analysis for background only data H b with respect to the null hypothesis H 0 = H s + H b where CRBDM is taken into account [72].We assume that the data in each bin follow a Guassian distribution when we construct the likelihood function, and ∆χ 2 asymptotically follows a χ 2 distribution with two degrees of freedom.The explicit distribution on ∆χ 2 could be obtained by Monte Carlo, which needs more inputs from the detector response beyond our current study.In a likelihood analysis, the excluded parameter space is obtained by χ 2 as follows: where χ 2 min is the minimum value of χ 2 and ∆χ 2 is associated with the number of considered parameters.When the ∆χ 2 exceeds 4.61 for two degrees of freedom, we reject the hypothesis H 0 at a 90% confidence level and obtain an exclusion limit in the parameter space.To be conservative, the DM mass and maximally allowed cross section to satisfy conditions in Eqn.(17) are represented as a point in the parameter space in Fig. III.Scanning over the interested low-mass DM range, we then figure out the exclusion limit for the spin-independent cross section as a lower part.Nevertheless, the Earth attenuation effect will dilute the strong exclusion limit since the bounced CRBDM might not be able to reach the detector threshold.This caveat will offer the upper ceiling limit.The eventual constraints on the DM parameter space will be obtained by a combination of exclusion limits for the lower part from an analysis of recoil spectrum in the DM direct detection experiment and the upper part from the Earth attentuation effect.
Following this strategy, we take the CDEX-10 experiment as a demonstration and reanalyze their data to obtain constraints on the DM mass m χ and the cross section σ χ0 in the CRBDM hypothesis at a 90% C.L. with two benchmark examples: D eff = 1 kpc and D eff = 3.7 kpc as shown in Fig. 4. Compared with the Migdal effect, the CRBDM framework will reach a much lower mass region but suffers from the Earth attenuation effect.However, cosmological constraints like the cosmic microwave background (CMB) [68] and the gas cloud cooling [69] help to exclude the parameter space above the reach of DM direct detection experiments.For the DM mass range 1 keV < m χ < 1 MeV and D eff = 1 kpc, we reach the almost flat floor limit at 8.32 × 10 −30 cm 2 on spin-independent DM-nucleon scattering cross section, while this limit turns up for m χ > 1 MeV and converges towards the ceiling limit caused by Earth attenuation.Here, we emphasize that the energy threshold also plays a crucial role for the attenuation effect.The low threshold of the CDEX detector somehow compensates for the sensitivity loss caused by the attenuation, though the CDEX experiment might suffer a stronger attenuation effect due to the thicker overburden.Therefore, the ceiling of the CDEX constraint can be higher than that of the XENON1T constraint.As shown in Fig. 4, the gray region enclosed by solid and dashed black lines stand for the constraints under the ballistic trajectory approximation with and without considering the form factor, respectively.We notice that the constraint without the form factor suppression in the attenuation calculations is more conservative, except for minor difference around the right endpoint.Next we come to the question Why the form factor helps to extend the ceiling limit.DM particles with higher energies suffer less energy loss when going through the crust due to the form factor suppression.Then more DM particles can reach the underground detector and the attenuation effect is relieved.The current result from CDEX-10 is one order of magnitude better than the neutrino experiment MiniBooNE but not comparable to the XENON1T data.Due to the attenuation, there will be a narrow gap for the excluded parameter regions from MiniBooNE and XENON1T.Given a larger D eff , the floor limit for each experiment will become stronger, while the ceiling limit induced from the Earth attenuation will hardly make changes.In a word, our analysis is more conservative compared with the previous study.The effective distance can enlarge the excluded regions for MiniBooNE and XENON1T to narrow the gap or even close it.It is then valuable to make a further study on D eff in detail, especially for the uncertain CR spatial distribution.With the current CRBDM hypothesis, the CDEX-10 result can close the gap in a nice manner.This fact highlights the importance of multiple detection technologies and combined analysis to examine the CRBDM mechanism.
Distinct from the ballistic trajectory assumption we take in deriving the attenuation effect, several studies performed a Monte Carlo (MC) simulation on the scattering process in the crust in detail [58,59,[73][74][75].The MC simulation is able to calculate the deflection of DM going through the Earth crust, which results in fewer DM particles reaching the underground detector.Thus, compared with the ballistic trajectory approach considering the form factor in the attenuation calculation, the MC simulation provides a slightly more conservative constraint.In addition, in the ballistic approach without considering the form factor, there is a suppression at the higher energy range (∼ 0.1 GeV).Due to the suppression, DM can not generate enough recoil energy and be detected when the cross section climbs over the ceiling limit.Since DM with higher energies can go through the crust for the form factor suppression, both the MC method and the ballistic approach in this case lead to too many recoil events.Thus, in order to be consistent with the experimental result, it is needed to enhance the attenuation effect in terms of a larger DM cross section.
As for the smaller DM cross section, the floor limit suffers less from the Earth attenuation compared with the ceiling bound.Therefore, the differences among the floor limits derived from various approaches should be smaller than that of the ceiling limits.However, for the larger cross section or the thicker crust leading to a stronger attenuation effect, the differences become larger and more careful treatment is needed to calculate the Earth attenuation.See Appendix A for more details.The data analyses in different experiments based on the MC simulation sets complementary constraints on the parameter space and the global floor limit comes down to about 3 × 10 −32 cm 2 on spin-independent DM-nucleon scattering cross section for the DM mass below 0.1 GeV.For the upper bound, the constraint derived from the MC simulation can exclude the parameter space between the floor limit and excluded regions from CMB and the gas cloud cooling constraints.
Coherent neutrino-neucleus scatterings will become intrinsic backgrounds of future DM direct detection experiments as soon as they make non-negligible contributions to nuclear recoils in the detector.Compared to the Xenon-based experiment, Germanium detector can obtain more energy transferred from neutrinos with tiny masses and has a lower energy threshold, which leads to more neutrino backgrounds.Taking neutrino backgrounds into account, we calculate the sensitivity of future germanium detectors on CRBDM with an exposure of about 1 ton•year and an energy threshold as low as 160 eVee in case of CDEX upgrades.The corresponding result of the neutrino floor is shown at the bottom area in Fig. 4. For the DM mass range 1 keV < m χ < 1 MeV and D eff = 1 kpc, it is found that we reach the neutrino floor at the cross section of about 3 × 10 −32 cm 2 in future Germanium-based DM experiments.In other words, the large-scale Germanium-based detector with a lower threshold might observe CEνNS events and simultaneously facilitate the study of astrophysical neutrino physics.

IV. CONCLUSIONS
Given that DM could elastically scatter with target nuclei in detectors, it is straightforward to think of the collisions between relativistic cosmic rays and low-speed DM in the galaxy.In this way, CRBDM would become more energetic to make the recoil energy above the threshold in detectors, especially for DM in the low mass range.Previous studies were conducted to constrain CRBDM for various experiments including MiniBooNE and XENON1T.There exists a gap in the excluded parameter regions between MiniBooNE and XENON1T partially due to the fact that the Earth attenuation effect will dilute the experimental constraints.Such a gap strongly depends on the effective distance determined by CR models where we need further investigations on the CR spatial distribution.We have reviewed the theoretical framework for the CRBDM, including the modified DM fluxes, the Earth attenuation effect from the overburden and the energy recoil spectrum for the general detector target nuclei.Taking the CDEX-10 experiment into account, we have re-analyzed 102.8 kg×day data to constrain CRBDM, and obtained constraints on the spin-independent DM-nucleon cross section in the low mass range 1 keV < m χ < 1 GeV, surpassing the detector energy threshold limitation in a classical analysis.Migdal effect can also lead to constraining the light DM but is still not comparable to the CRBDM mechanism.We have found that the CDEX-10 result is able to close the gap unambiguiously regardless of the assumed effective distance in the parameter space between MiniBooNE and XENON1T constraints.
For future direct detection experiments, the same strategy to constrain CRBDM could be applied.Note that it is still necessary to take careful treatment of uncertainties from astrophysical and nuclear physics inputs as well as the detector response in order to give more rigorous results on CRBDM.In addition, we also quantitiviely provide the neutrino floor on searching for CRBDM in future direct detection experiments using Germanium.Our work should encourage the experimental collaborations to press forward with development of the next generation Germanium detector and joint analysis with other experiments using different detection technologies.Recently, the XENON1T collaboration reported the low-energy ER excess [76].CRBDM might be able to explain the observed ER excess in the LXe detector [77].We look forward to further examination by next-generation low-threshold DM direct detection experiments.

Figure 1 .
Figure 1.The schematic diagram for acceleration and detection processes in a CRBDM mechanism.

1 ]Figure 3 .
Figure3.The measured energy recoil spectrum in CDEX-10 with black data points plus error bars[13], a background assumption in blue dashed lines the predicted CRBDM recoil spectrum in red solid line for mχ = 1 MeV and σχ0 = 7 × 10 −30 cm 2 .Here "eVee" represents electron equivalent energy derived from a charge calibration.

2 ]Figure 4 .
Figure 4. Constraints on spin-independent cross section of DM-nucleon interactions from CDEX-10 (gray), XENON1T (purple) and MiniBooNE (blue)[33] in the framework of CRBDM at a 90% C.L. The gray region enclosed by solid and dashed black lines stand for constraint from CDEX-10 under the ballistic trajectory approximation with and without considering the form factor, respectively.As a comparison, we also include the exclusion limit from EDELWEISS[29] based on the Migdal effect and from CMB[68] (green) as well as the gas cloud cooling[69] (orange).The neutrino floor (brown) is obtained with an exposure of 1 ton•year and a recoil energy threshold of 1 keV for the Germanium detector.

Figure 5 .
Figure5.The CRBDM fluxes and theoretical spectra derived from the ballistic approach (solid line) and the MC method (dashed line).Form factor effects are included here in all panels.

Table I .
[52]meters of the analytical fits to the proton and He LIS[52].