Resolving XENON excess with decaying cold dark matter

We propose a decaying cold dark matter model to explain the excess of electron recoil observed at the XENON1T experiment. In this scenario, the daughter dark matter from the parent dark matter decay easily obtains velocity large enough to saturate the peak of the electron recoil energy around 2.5 keV, and the observed signal rate can be fulfilled by the parent dark matter with a mass of order 10–200 MeV and a lifetime larger than the age of Universe. We verify that this model is consistent with experimental limits from dark matter detections, Cosmic microwave background and large scale structure experiments.


Introduction
Recently, XENON1T experiment [1], which is a dark matter (DM) direct detection facility, has reported an excess of electron recoil over the background in the 1-7 keV range with 3.5σ significance. As pointed out in Ref. [1], this excess is unlikely due to solar axion, neutrino magnetic moment or statistical uncertainties about the background. So far, the observed excess has initiated extensive investigations about potential astrophysical sources. Among other things, a cold DM is a natural candidate, which is the subject of this study.
To explain the observed excess, one has to accommodate two critical quantities -the electron recoil energy range around 2-3 keV and the electron transfer momentum range near 50 keV. Unfortunately, they conflict with a conventional cold DM, which has a velocity typically of order ∼ 10 −3 c, with c the velocity of light. A few proposals have been proposed to avoid the conflicts. In the case of elastic scattering [2][3][4][5][6][7], the cold DMs can be boosted in certain circumstances before they interact with the electrons in the xenon atoms, while in the case of inelastic scattering [8][9][10][11][12] the favored electron transfer momentum range can be realized in terms a e-mail: shuaixu@cqu.edu.cn b e-mail: sibozheng.zju@gmail.com (corresponding author) of small rest mass splitting between two different DM components.
In this study, we propose a new decaying cold DM scenario, in which the parent DM (A) decays to the daughter DM particles (B), In terms of the decay, the velocity of the daughter DM can be enhanced to be comparable with c, relative to the small velocity of the parent DM. Unlike photoelectric absorption in a decaying warm DM [13], the daughter particle elastically scatters off the electrons in the xenon atoms. We will show that without any violations of current DM (in)direct detections or cosmological measurements, this decaying cold DM model can easily resolve the XENON excess.

The model
We begin with the production of the daughter particle B due to A decay. The decay yields the following velocity and present number density of B respectively, where m A and m B refer to the parent and daughter DM mass respectively with m A > 2m B , t 0 is the age of Universe, while ρ dm,0 = 0.4 GeV/cm 3 and A denote the local DM density and the decay width of the parent DM A, respectively. We will assume that compared to the decay production the thermal production of n B is subdominant. The input parameters in Eqs. (2) and (3) are constrained as follows. Firstly, in order to yield υ B of order ∼ 0.1 c the mass ratio m A /2m B should deviate from unity in percent level, which implies that m A and 2m B are highly degenerate.
Secondly, in order to fulfill the cosmological bounds on the decaying DM both from the Cosmic microwave background (CMB) and large scale structure (LSS) experiments, the lifetime of A particle τ A = −1 A should be larger than t 0 , which can be achieved by adjusting the coupling constant between A and B, with the help of a suppression by the small β factor due to the mass degeneracy.
We will return to the cosmological constraints after we have explored the signal rate of the recoil electrons at the XENON1T.

Signal rate
According to conservations of energy and momentum in the elastic scattering process, the energy 1 transferred to electron reads [14] Consider that when m B ≤ m e the daughter DM with a large velocity is severely constrained by limits such as the effective number of neutrinos, we will focus on m B m e , under which E max e ≈ 2m e υ 2 B [2] instead. This constraint implies υ B ≥ 0.05 c in order to satisfy E max e ≥ 2.5 keV. Given a fixed value of E e , Eq. (4) determines the electron transfer momentum range q − < q < q + , with The transfer momentum range in Eq. (5) affects the signal rate of the recoil electrons discussed below in the sense that the atomic factorization factor K (E e , q) [15,16] is rather sensitive to q, which takes the maximal value K max ≈ 0.1 at q peak ≈ 50 keV for E e = 2 keV, and dramatically declines as q slightly deviates from q peak . Therefore, in order to maximize the K -factor contribution to the signal rate, we should take suitable values of m B and υ B to make sure that q peak is covered by the electron transfer momentum range in Eq. (5). Furthermore, the daughter DM-free electron scattering cross sectionσ e relies on the nature of mediator [17] which communicates the interaction between the daughter DM B and electron. From the viewpoint of effective field theory,σ e can be written as where m med is the mediator mass, g m B is the coupling between the mediator and B, and g me is the coupling between the mediator and electron. If the mediator is identified as a standard model particle, only g m B in Eq. (6) is a free parameter, the magnitude of which has to be constrained by the decay width of the standard model particle. After a handle on the "luminosity" and the DM B-electron scattering cross section, we now estimate the number of events of recoil electrons where n xe ≈ 4.2 × 10 27 /ton is the number density of xenon atoms in the detector, a 0 = 1/(α em m e ) is the Bohr radius with α em = 1/137, F(q) ≈ 1 is the DM form factor, and R s is the resolution function which accounts for the "efficiency" of the detector. We will simply take the Gaussian distribution for the reconstructed energy for numerical analysis where α(E) is the efficiency [1] and σ = a √ E e + bE e , with a = (0.310 ± 0.004) √ keV and b = 0.0037 ± 0.0003, respectively. Table 1 Five sets of benchmark values which yield the same fit as shown in Fig. 1, where the required values ofσ e can be understood as an output parameter  Figure 1 shows the fit to the reported XENON1T data [1] with five different sets of benchmark values C1 to C5 as explicitly shown in Table 1

Dark matter constraints
Now we turn to possible constraints on the dark matter particles A and B. Since the interaction in Eq. (1) yields too small annihilation cross section for A to accommodate the required thermal annihilation cross section, A has to communicate either with the Standard Model (SM) sector e.g. via the same mediator as B, or mainly with other unstable freedoms in the dark sector. In the former situation, some constraints on B as below can be placed on A as well.
With the communication between B and the SM sector as inferred from the XENON1T excess, we can at least place the following constraints.
• The daughter DM B-free electron scattering cross section σ e , extracted from the XENON1T excess, can be used to constrain the model parameters. Based on the measurements onσ e within various electron recoil energies, the light daughter DM can be probed either by the current XENON1T [18,19] or the future SuperCDMS [20] experiments.
• Similar to the DM B-electron scattering, we can also constrain m B from the annihilation cross section σ B ann (B B → e + e − ), based on the cross symmetry between the two Feynman diagrams related to these two processes. While experiments such as AMS-01 [21], AMS-02 [22] or PAMELA [23] have not yet placed viable bounds on σ B ann in the sub GeV-scale DM mass, the Planck data [24] is able to constrain m B down to ∼ 1 MeV. For illustration, we show in Fig. 2 the constraints in the specific dark photon model with the mediator identified as a new vector boson A , where m A = 1 GeV and g me = 5×10 −4 e have been adopted in the light of BaBar data [25]. In this figure, we have simultaneously shown the Super-CDMS limit [20] (in blue) without relativistic effect [26], the Planck 2015 limit [27] (in red) and the parameter space (in black curve) together with the benchmark values in Table 1. Relatively weaker XENON1T limit has been ignored. One observes that in this explicit model m B beneath ∼ 15 MeV survives.
Compared to the benchmark values in Table 1, given fixed m B one can obtain largerσ e or alternatively larger g m B g me by taking larger τ A , since they are linearly correlated to each other in d R/d E ∼σ e (t 0 /τ A ) for t 0 τ A in Eq. (7). However, an increase of τ A will simultaneously lead to linearly enhanced experimental limits in Fig. 2. These trends together imply that adjusting τ A is unable to alter the SuperCDMS sensitivity as illustrated in Fig. 2.

Cosmological constraints
The decaying DM model are constrained both by the CMB and the LSS experiments for a varying dark matter energy density with time. In our scenario, it reads from Eq. (3) Compared to the baseline CDM cosmology, the DM relic density in Eq. (9) is altered by a magnitude of order | ρ dm /ρ dm,0 |≈ (1 − 2m B m A )t/τ A < 10 −3 in the small redshift region for the benchmark values in Table 1, as a result of highly degenerate dark matter mass relation m A ≈ 2m B required by the XENON excess.
For the CMB experiment [28], it mainly affects the temperature power spectrum C T T in terms of the integrated Sachs-Wolfe effect, which relies on the cosmological evolution of Universe after the last scattering. Due to the small fraction in ρ dm given by Eq. (9) relative to what attempts to explain the Hubble tension [29][30][31][32], the effect on C T T in our scenario is negligible. For the LSS experiments, the DM power spectrum δ = δρ dm /ρ dm evolves with time as where G is the Newton's constant and H is the Hubble rate. The small fraction in ρ dm gives rise to a fraction in the DM power spectrum δ less than the order of a percent level, which is far beyond the reach of future LSS experiments such as the Dark Energy Spectroscopic Instrument [33].

Conclusions
In this study we have proposed a novel decaying cold DM scenario in which the cold parent DM A decays to the daughter particle B, with the lifetime τ A larger than the age of Universe. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .