Lightening Gravity-Mediated Dark Matter

We revisit the scenario of a massive spin-2 particle as the mediator for communicating between dark matter of arbitrary spin and the Standard Model. Taking the general couplings of the spin-2 particle in the effective theory, we discuss the thermal production mechanisms for dark matter with various channels and the dark matter self-scattering. For WIMP and light dark matter cases, we impose the relic density condition and various experimental constraints from direct and indirect detections, precision measurements as well as collider experiments. We show that it is important to include the annihilation of dark matter into a pair of spin-2 particles in both allowed and forbidden regimes, thus opening up the consistent parameter space for dark matter. The UV complete models of the spin-2 mediator are presented in the context of the warped extra dimension and compared to the simplified models.


Introduction
Dark matter (DM) is a complete mystery in particle physics and cosmology, although its presence can be unambiguously inferred from galaxy rotation curves, gravitational lensing, Cosmic Microwave Background as well as large-scale structures, etc. There are null results in searching dark matter beyond gravitational interactions from various direct and indirect detection experiments, thus, in particular, a lot of parameter space for Weakly Interacting Massive Particles (WIMPs) has been ruled out [1][2][3].
The nature of dark matter is still an open question. To this, it is very important to pin down the production mechanisms for dark matter in the early universe. For instance, WIMP dark matter relies on the freeze-out process under which the DM relic density is determined in terms of weak interaction and weak-scale DM mass. Thus, this has motivated specific target materials and technologies in the direct searches for WIMP for more than three decades. New production mechanisms such as for Feebly Interacting Massive Particles (FIMPs) [4], Strongly Interacting Massive Particles (SIMPs) [5,6] and forbidden dark matter [7,8], etc, can motivate different target materials and new technologies to get access to sub-GeV DM masses and/or feeble interactions. It is known that light dark matter with sub-GeV mass can have large self-interactions to solve potentially small-scale problems at galaxies [9,10] and it may also call for new dynamics in the dark sector [11] to get the DM self-interactions velocity-dependent for galaxy clusters such as Bullet cluster [12].
Moreover, dark matter is known to be neutral under electromagnetism, so it is conceivable to communicate between dark matter and the Standard Model (SM) through messenger or mediator particles. Thus, the simplified models for dark matter with mediator particles have drawn a lot of attention, providing an important guideline for direct and indirect detections of dark matter as well as collider experiments [13,14].
In this article, we consider a massive spin-2 particle as the mediator for dark matter of arbitrary spin, which couples to the SM particles and dark matter through the energymomentum tensor, as originally proposed by one of us and collaborators [15]. This scenario has been dubbed "Gravity-mediated dark matter", due to the similarity to the way that the massless graviton interacts with the SM. The spin-2 mediator stems from a composite state in conformal field theories or a Kaluza-Klein(KK) graviton in a gravity dual with the warped extra dimension [15][16][17]. There are other works on the spin-2 mediated dark matter in similar frameworks [19]. We treat the interactions of the massive spin-2 particle in the effective theory with general couplings to the SM and dark matter and discuss the general production mechanisms for WIMP dark matter and light dark matter in this scenario.
We discuss various channels of dark matter interactions in the presence of the spin-2 mediator: direct 2 → 2 annihilations, 2 → 2 allowed and forbidden channels into a pair of spin-2 mediators, 3 → 2 assisted annihilations as well as DM self-scattering. We not only make a complete analysis of the DM-nucleon elastic scattering in the presence of quark and gluon couplings by extending the previous results in Ref. [18] but also provide new results for the DM-electron elastic scattering. We check the consistency between the correct relic density and various experimental constraints, such as direct detection, precision measurement of muon g − 2, meson decays and collider experiments, in both WIMP and light dark matter cases. We also introduce two benchmark models with the warped extra dimension as an UV completion of the spin-2 mediator, such as the Randall-Sundrum(RS) model [20] and the clockwork model [54]. Then, we discuss the impacts of heavier KK gravitons on the aforementioned processes for dark matter.
The paper is organized as follows. We begin with a brief description of our setup for the spin-2 mediator and its interactions. Then, we determine the DM relic density from various annihilation channels and discuss the self-scattering process for dark matter. Next, we consider the DM-nucleon elastic scattering for WIMP and the DM-electron elastic scattering for light dark matter and provide various direct and indirect constraints on those dark matter models. We continue to show two benchmark models with the warped extra dimension and discuss how the DM processes can be modified due to extra resonances. Finally, conclusions are drawn. There are three appendices dealing with the details on DMnucleon scattering amplitudes, decay widths of spin-2 particles as well as the KK sums.

The setup
We consider the effective interactions of a massive spin-2 field, G µν , to the SM particles as well as dark matter with arbitrary spin, in the following [15], where B µν , W µν , g µν are the strength tensors for U (1) Y , SU (2) L , SU (3) C gauge fields, respectively, ψ is the SM fermion, H is the Higgs doublet, and Λ is the dimensionful parameter for spin-2 interactions. Here, we note that c i (i = 1, 2, 3), c ψ , and c H are dimensionless couplings for the KK graviton. Depending on the spin of dark matter, s = 0, 1 2 , 1, denoted as S, χ and X, the energy-momentum tensor for dark matter, T DM µν , is given, respectively, by In the later discussion, we focus on the couplings of the spin-2 mediator to quarks, leptons and massless gauge bosons in the SM, as well as dark matter couplings. We treat those SM to mediator couplings to be independent parameters, but be universal for simplicity as well as unitarity consideration.

Dark matter annihilations and self-scattering
In this section, we discuss the Boltzmann equations for determining the relic density of dark matter and show the details for the cross sections for 2 → 2 direct annihilations. In particular, we obtain for the first time the new results for 2 → 2 forbidden channels, 3 → 2 assisted annihilations, and DM self-scattering.
First, we consider the Boltzmann equations for the relic density of real scalar dark matter S or vector dark matter X, given bẏ Similarly, for Dirac fermion dark matter χ, the corresponding Boltzmann equation for n DM = n χ + nχ iṡ Henceforth, we assume that the spin-2 particle is in thermal equilibrium with the SM plasma during the freeze-out, so we can take n G = n eq G , which is the number density in thermal equilibrium.

Direct annihilations
We focus on the cases with relatively light WIMP dark matter and light dark matter below the WW threshold, which annihilate dominantly into the SM fermions or massless gauge bosons.
If dark matter is heavier than the W W threshold, we can also take into account the DM annihilations into the electroweak sector, as shown in Ref. [15,16], allowing for smaller couplings of the spin-2 mediator to the SM particles for a correct relic density. In this work, however, for WIMP dark matter, we take the spin-2 mediator couplings to the SM quarks and gluons to be nonzero in simplified models. For consistency of gauge-invariant couplings, we choose c 1 = c 2 = c H = 0 in the electroweak sector and c l = 0 for SM leptons in the discussion for WIMP. On the other hand, for light dark matter below the W W threshold, we keep all the spin-2 mediator couplings to the SM to be nonzero.
In the case when dark matter is heavier than the spin-2 mediator, dark matter can also annihilate directly into a pair of spin-2 particles, reducing the dark matter abundance further together with the direct annihilations into the SM.
In the case where 2 → 2 annihilation channels are dominant, the Boltzmann equations, (3.1) or (3.2), becomeṅ Then, the relic density for WIMP dark matter is given by

Scalar dark matter
The annihilation cross section for scalar dark matter into a pair of SM fermions, SS → ψψ, is given [15][16][17] by where N c is the number of colors for the SM fermion ψ, and Γ G is the width of the spin-2 particle. Thus, the annihilation of scalar dark matter into the SM fermions becomes d-wave suppressed, so scalar dark matter is not constrained by indirect constraints from cosmic rays and Cosmic Microwave Background (CMB) measurements [15,16]. When m S > m G , scalar dark matter can also annihilate into a pair of spin-2 particles through the t/u-channels [15][16][17][18], becoming dominant due to sizable spin-2 couplings to dark matter. Then, the corresponding annihilation cross section is given, as follows, For light dark matter, the DM annihilations into photons or gluons are relevant. For sub-GeV dark matter, the DM annihilations into mesons must be considered instead of those into gluons. Then, for scalar dark matter, the annihilation cross sections into a pair of massless gauge bosons [15] are For 2m S 1.5 GeV, instead of the annihilation into a gluon pair, we should consider the annihilation cross section of scalar dark matter into a meson pair, as follows, where c π c q in the limit of small momenta of produced pions, because the chiral perturbation theory takes in. We also need to include the annihilation of scalar dark matter into charged pions and kaons, if kinematically allowed.

Fermion dark matter
The annihilation cross section for fermion dark matter, χχ → ψψ, is given [15][16][17] by Thus, the annihilation of fermion dark matter into the SM fermions becomes p-wave suppressed. Then, similarly to the case of scalar dark matter, fermion dark matter is not constrained by indirect constraints from cosmic rays and CMB measurements [15,16]. When m χ > m G , fermion dark matter also annihilates into a pair of spin-2 particles through to the t/u-channels [15][16][17], as follows, . Then, the resulting annihilation cross section is s-wave, so it becomes dominant in determining the relic density for fermion dark matter.
For light fermion dark matter, the annihilation cross sections into a pair of massless gauge bosons and a pair of mesons [15] are (3.14) For 2m χ 1.5 GeV, we need to include the annihilation channel into a pion pair by Similarly, the annihilation of fermion dark matter into charged pions and kaons, if kinematically allowed, should be also included.

Vector dark matter
The annihilation cross section for vector dark matter, XX → ψψ, is given [15][16][17] by Thus, the annihilation of vector dark matter into quarks becomes s-wave. In this case, smaller spin-2 mediator couplings to the SM quarks or vector dark matter can be consistent with the correct relic density, as compared to the other cases. In this case, the CMB measurement for recombination era can rule out the vector dark matter mass below 100 GeV, if the relic density is determined solely by the direction annihilation into the SM particles. But, indirect detection signals from the annihilation of vector dark matter are promising [15,16].
For m X > m G , vector dark matter also annihilates into a pair of spin-2 particles through the t/u-channels [15][16][17][18], as follows, For light vector dark matter, the annihilation cross sections into a pair of massless gauge bosons and a pair of mesons [15] are For 2m X 1.5 GeV, we also need to include the annihilation into a pion pair by (3.20) Similarly, the annihilation of vector dark matter into charged pions and kaons, if kinematically allowed, should be also included.

Forbidden channels
When dark matter is lighter than the spin-2 mediator, but their masses are comparable, that is, m DM m G , the annihilation of dark matter into a pair of spin-2 particles is forbidden at zero temperature, but it is kinematically allowed due to the tail of the Boltzmann distribution of dark matter at finite temperature, making the so called forbidden channels relevant for determining the DM abundance. In this subsection, we consider the forbidden channels in association with the spin-2 mediator.
In the case when the forbidden channels are dominant, the Boltzmann equations, (3.1) and (3.2), becomeṅ where the forbidden annihilation cross sections are given by Here, ∆ G = (m G − m DM )/m DM , and g DM is the number of degrees of freedom of dark matter, g DM = 1, 4, 3, for real scalar, Dirac fermion and vector dark matter, respectively. Here, we have used the detailed balance condition for forbidden channels. Moreover, for m DM < m G , the cross sections for the inverse annihilation channels are given by As a result, the relic density for forbidden dark matter [8] is given by There is a Boltzmann suppression factor in the effective annihilation cross sections for forbidden channels, so we would need larger couplings of dark matter to the spin-2 mediator for the correct relic density, as compared to the case with allowed 2 → 2 channels for m DM > m G .

Gravity-mediated 3 → 2 processes
Scalar dark matter can annihilate by SSS → SG, which can be dominant over the forbidden channels, SS → GG, for m S < m G < 2m S . Similarly, the 3 → 2 processes for fermion dark matter (χχχ → χG) and vector dark matter (XXX → XG) can be important for m χ < m G < 2m χ and m X < m G < 2m X , respectively. Thus, we choose m DM < m G < 2m DM in order for the 3 → 2 processes to be kinematically open and for the hidden sector 2 → 2 annihilations to be forbidden. In this subsection, we consider the assisted 3 → 2 channels with the spin-2 mediator for the first time.
When the 3 → 2 annihilation processes are dominant, the Boltzmann equation (3.1) becomesṅ Here, the corresponding 3 → 2 annihilation cross sections for scalar and fermion dark matter are As a result, the relic density for SIMP dark matter [6,8] is given by We note that the 3 → 2 annihilation cross sections are highly suppressed for perturbative couplings in most of the parameter space, so they are sub-dominant in determining the relic density, as compared to the previously discussed 2 → 2 annihilation channels. Therefore, we don't consider the SIMP option in the later discussion.

Dark matter self-scattering
Spin-2 mediator can also mediate the self-scattering process of dark matter, in particular, for fermion and vector dark matter, for which there is no renormalizable interaction for self-scattering. We can take the gravity-mediated processes to be dominant for dark matter self-scattering and consider the interplay between relic density condition and small-scale problems in galaxies.
For scalar dark matter, the self-scattering cross section for SS → SS, divided by DM mass, is in the Born approximation (3.33) For fermion dark matter, the self-scattering cross section from χχ → χχ and χχ → χχ (and its complex conjugate), divided by DM mass are similarly given by (3.34) Finally, for vector dark matter, the self-scattering cross section for XX → XX, divided by DM mass, is given by We note that for both scalar and fermion dark matter, the DM self-scattering cross section little depends on the DM velocity. In the case of scalar dark matter, there is an s-channel contribution with the spin-2 mediator too, but it is velocity-suppressed by the overall factor. On the other hand, for vector dark matter, the DM self-scattering cross section could be enhanced at a particular DM velocity due to the s-channel resonance [22], so it would be possible to accommodate the velocity-dependent self-interaction, being compatible with galaxy clusters such as Bullet Cluster [12].

Detection of dark matter and mediator couplings
We give the phenomenological discussion on the spin-2 mediator for DM-nucleon elastic scattering, DM-electron elastic scattering, g − 2 of leptons, meson decays and the direct production at colliders. We present for the first time the complete discussion of DM-nucleon scattering in the presence of both quark and gluon couplings and DM-electron scattering as well as the relevance of unitarity at colliders.

DM-nucleon elastic scattering
The scattering amplitude between DM and SM particles through the spin-2 mediator [18] is written in the limit of a small momentum transfer, as follows, is the traceless part of energy-momentum tensor given byT First, the elastic scattering amplitude between dark matter and nucleon [18] is given by For direct detection experiments, we can consider only the contributions from quarks and gluons in a nucleon, as follows, Then, we get the trace part in the effective theory for three quark flavors (u, d, s) and gluons as where scale anomalies from light quarks and gluons are separately taken into account. Moreover, the traceless part (twist-2 operators) for five quark flavors (u, d, s, c, b) and gluons is given by As a result, the nuclear matrix elements for the trace part become where f N T q , f T G are the mass fractions of light quarks and gluons in a nucleon, respectively, and f T G = 1− q=u,d,s f N T q . Here, we used the RG invariant quantity, N (p)|α S G µν G µν |N (p) = − 8π 9 f T G m N , which is obtained in the effective theory for three quark flavors. For the universal spin-2 couplings with c q = c g , we obtain the standard results for On the other hand, the nuclear matrix elements for the traceless part [23] are where q(2),q(2) and G(2) are the second moments of the parton distribution functions(PDFs) of quark, antiquark and gluon, respectively, The mass fractions are f p Tu = 0.023, f p T d = 0.032 and f p Ts = 0.020 for a proton and f n Tu = 0.017, f n T d = 0.041 and f n Ts = 0.020 for a neutron [23]. On the other hand, the second moments of PDFs are calulated at the scale µ = m Z using the CTEQ parton distribution as G(2) = 0.48, u(2) = 0.22,ū(2) = 0.034, d(2) = 0.11,d(2) = 0.036, s(2) =s(2) = 0.026, c(2) =c(2) = 0.019 and b(2) =b(2) = 0.012 [23].
There, using the results in the appendix A, the total cross section for spin-independent elastic scattering between dark matter and nucleus [18] is given by is the reduced mass of the DM-nucleus system and m A is the target nucleus mass, Z, A are the number of protons and the atomic number, respectively, and the nucleon form factors are given by the same formula for all the spins of dark matter as (4.14) where DM = χ, S, X for fermion, scalar and vector dark matter, respectively. Here, as compared to our previous work [18], we have included the twist-2 gluon operator at tree level as well as loop effects from heavy quarks and gluons in the trace part.

DM-electron elastic scattering
For light dark matter below GeV scale, the DM-nucleon elastic scattering loses the sensitivity for dark matter searches because of the low threshold of the nucleon recoil energy. Then, the DM-electron elastic scattering is relevant for direct detection [6]. The corresponding cross sections relevant for direct detection are independent of the spin of dark matter, given by where we assumed that m DM m e in the second line. Moreover, the graviton mediator should make dark matter remain in kinetic equilibrium [5,6] during the freeze-out. In this case, independent of the spin of dark matter, the momentum relaxation rate for the kinetic equilibrium of light dark matter is dominated by Then, the kinetic equilibrium of dark matter can be achieved during the freeze-out, as far as γ DM−e > H · m DM T in the case of WIMP dark matter where H is the Hubble expansion parameter, and γ DM−e > H · m DM T 2 in the case of SIMP dark matter [6,11].

Lepton g − 2 from the spin-2 mediator
When the spin-2 mediator couples to leptons, it gives an extra contribution to the anomalous magnetic moment of leptons, as follows [24], where A(y) is a monotonically decreasing function, given by For m G m l , the loop function A(x) is approximated [25] to rendering the (g − 2) l almost independent of the spin-2 mediator mass, as follows, We note that the deviation of the anomalous magnetic moment of muon between experiment and SM values is given [26,27] by which is a 3.6σ discrepancy from the SM [27]. Furthermore, there is a 2.4σ discrepancy reported between the SM prediction for the anomalous magnetic moment of electron and the experimental measurements [28,29], as follows,
The decay width of a down-type quark q 1 decaying into another down-type quark q 2 and G is given for m G < m q 1 with m q 2 = 0 [33], as follows, where V f 1 and V f 2 are the CKM matrix elements,

Mediator production at colliders
The massive spin-2 particle can be produced singly from gluon fusion or quark/anti-quark scattering at the LHC, decaying into the SM particles or a pair of dark matter. Moreover, in intensity beam or linear colliders, we may also constrain non-universal lepton and photon couplings by the photon energy distribution from e + e − → γ G. First, we obtain the squared amplitude for e + e − → γ G, as follows, cos θ). Therefore, for c γ = c e , the squared amplitude behaves like |M| 2 ∼ s Λ 2 for s m 2 G [34,35], which is expected from the dimension-5 interactions  Parameter space for m G /Λ vs m DM for WIMP dark matter. The relic density is satisfied in red solid, blue dashed and orange dotted lines for fermion, scalar and vector dark matter, respectively. The gray region is excluded by XENON1T and the light blue region is excluded by ATLAS dijet searches. We have taken the universal spin-2 mediator couplings to the SM and dark matter.
for the spin-2 mediator, − 1 Λ G µν T µν . However, for c γ = c e , the squared amplitude becomes |M| 2 ∼ s 3 m 4 G Λ 2 , which shows that the violation of unitarity at a lower energy. A similar phenomenon was observed in the QCD process, qq → g G [34,35], for which c g = c q would give rise to a similar dependence of the corresponding squared amplitude on the center of mass energy.
For c γ = c e , the production cross section for e + e − → γ G with unpolarized electron and positron is given by Thus, the angular differential cross section becomes independent of s for s m 2 G , as expected from the behavior of the squared amplitude. A similar conclusion can be drawn also for qq → g G at the LHC. The above result will be used for imposing the bounds from invisible and visible searches at BaBar in Fig. 8 of the next section.

Bounds on WIMP
Dijet and dilepton searches at the LHC can constrain relatively heavy spin-2 resonances [36].
Although not sensitive enough, the ISR photon or jet + heavy dijet resonances might be interesting to constrain non-universal quark and gluon couplings by the jet p T distribution from qq → g G at LHC and future hadron colliders [34]. Direct detection bounds from XENON1T [1], LUX [2], PandaX [3], etc, are most stringent for weak-scale or heavier dark matter.
For weak-scale spin-2 resonances, the LHC dijet searches are not sensitive due to the large QCD background. Then, dijet resonance + ISR photon [37] or jet [38,39] searches can constrain this case. In the presence of dark matter coupling to the spin-2 resonance, the invisible decay of the spin-2 particle with mono-jet of mono-photon is also promising [13,14,40,41].
In Figs. 1 and 2, we depict the parameter space for m G /Λ vs m DM in the former and m DM vs m G in the latter, satisfying the correct relic density, in red solid, blue dashed and orange dotted lines for fermion, scalar and vector dark matter, respectively. We took the universal couplings of spin-2 mediator to all the SM quarks and gluons, as well as to dark matter. We have excluded the light blue region by the bounds from dijet resonance + ISR photon [37] or jet [38,39] searches, and the gray region by the bound on DM-nucleon spin-independent cross section from the direct detection experiment in XENON1T [1].
We find from Fig. 1 that for weak-scale spin-2 mediator, the relic density region below m DM < m G /2 is disfavored by ATLAS dijet bounds. The XENON1T bound becomes stronger above m DM > m G /2, leaving only the region above m DM 200 GeV or larger masses unconstrained due to the dominance of DM DM → GG channels. But, in this case, the spin-2 mediator produced from the DM annihilation can decay into the SM particles, so the indirect detection experiments from cosmic rays such as positrons, anti-protons and gamma-rays can constrain those large mass regions [15]. In Fig. 2, XENON1T rules out the non-resonance regions below m DM 200 GeV or 160 GeV for the mediator scale, Λ/c q = 3, 5 TeV, but leaves the resonance regions with m G = 2m DM untouched.

Bounds on light dark matter
In the case of light dark matter, we would need a light spin-2 mediator in order to make the annihilation cross section of dark matter sufficiently large. In this case, monophoton + leptons at BaBar [42], and missing energy at BaBar [43], Belle-2 [44,45], LHCb (for m G > 10 GeV) [46] as well as beam dump experiments such as E137 in SLAC [47], N64 in CERN SPS [48], etc, can be important to constrain the light spin-2 mediator couplings, in particular, the couplings to leptons and dark matter. There are also direct detection bounds on DM-electron scattering from XENON10 [49], DarkSide-50 [50], Sensei experiments [51], etc.
For a light spin-2 mediator, we can consider the bounds from γ + missing energy [43] or leptons [42] at BaBar experiment. For the former case, the cosine of the scattering angle of the photon in the center of mass frame was chosen to | cos θ * γ | < 0.6, and the center of mass energy was √ s = 10.58 GeV. Then, we get the limit on the lepton couplings for m G < 8 GeV from invisible and visible searches at BaBar, respectively, as follows, Here, we assumed BR(G → DM DM) = 1 in the former and BR(G → ll) = 1 in the latter. So, in general, the above bounds scale up by 1/ √ BR. The above limits, in particular, from the invisible searches, will be improved by a factor of three in the lepton couplings in Belle-2 experiment [44,45].
We remark that if we took non-universal couplings by c γ = c e , the above bounds from BaBar would become stronger, due to the growth of the corresponding cross section.
Moreover, if the spin-2 mediator is much lighter than K-meson or B-meson, we can approximate the above partial decay rate of a flavor-changing down-type quark from eq.(4.25) to Therefore, from the current limits on the invisible decays of K + or B + , we can put the bound on the quark couplings as As a result, the bounds on quark couplings from meson decays are relatively weaker than those on lepton couplings from BaBar as will be shown in the above. When the spin-2 mediator is heavier than mesons but dark matter is light enough, mesons can still decay invisibly into a pair of dark matter [32]. But, in this case, the bounds on quark couplings become much weaker because of the phase-space suppression for three-body decays of mesons. In Figs. 3, 4 and 5, we show the parameter space for light dark matter below the GeV scale mass satisfying the correct relic density, in c e /Λ vs m G in the former and c DM /Λ vs m G in the latter two. For Fig. 3, we took m DM < m G such that dark matter annihilates only into the SM particles, not into a pair of spin-2 mediators. In this case, we find that the graviton couplings to the SM particles satisfying the correct relic density would be strongly constrained by BaBar and other intensity experiments, except the region near the resonance. On the other hand, for Figs. 4 and 5, we took m DM > m G for which dark matter can annihilate into a pair of spin-2 mediators. In this case, even for a small graviton coupling to the SM particles, for instance, for Λ/c e = 10 TeV or 100 TeV in Figs. 4 or 5, for which the current experimental constraints are satisfied, we can achieve the correct relic density in a wide range of parameter space for dark matter coupling and spin-2 mediator mass. We have also checked in Figs. 3, 4 and 5 that the DM self-scattering cross sections in the parameter space explaining the relic density are much below σ self /m DM = 1 cm 2 /g, the Bullet cluster bound [12]. We note that the difference between DM (c DM /Λ) and lepton couplings (c e /Λ) can be explained by the localization of dark matter and leptons in different positions of the extra dimension. For instance, in RS model, light dark matter can be localized on the IR brane with a small IR scale whereas the SM leptons are localized towards the UV brane [15].
In Figs. 6 and 7, we present the relic density as a function of the mass difference, ∆ G ≡ (m G − m DM )/m DM , with forbidden channels included. These plots illustrate the role of the forbidden channels in determining the relic density for the spin-2 mediator slightly heavier than dark matter. In this case, the annihilation of dark matter into a pair of spin-2 mediators is possible only at a nonzero temperature, thus leading to a Boltzmann suppression factor for the corresponding annihilation cross section. For each of Figs. 6 and 7, we have chosen m DM = 1, 10 GeV on left and right. We took Λ/c DM = 10 GeV for both, and Λ/c e = 10 TeV, 100 TeV for Figs. 6 and 7, respectively.
We find that the correct relic density for vector dark matter can be obtained with smaller couplings to the spin-2 mediator and sub-GeV DM massses, due to a mild phasespace suppression for m G m DM . On the other hand, for scalar or fermion dark matter,  dark matter masses should be about 10 GeV or larger for the correct relic density being consistent with perturbativity, due to significant phase-space suppressions for m G m DM . The forbidden channels are s-wave but get suppressed as the velocity of dark matter decreases in the later stage of the universe and in local galaxies. Thus, the forbidden channels are safe from the indirect bounds from cosmic rays or CMB measurements. In particular, it is remarkable that sub-GeV vector dark matter with m DM m G can be consistent with both the relic density and indirect detection bounds, being compatible with perturbativity.
In Fig. 8 Λ/c DM = 100 GeV on right. We note that for both plots of Fig. 8, the DM self-scattering cross sections in the parameter space of our interest are well below the Bullet cluster bound.
In the left plot of Fig. 8, the spin-2 mediator can decay dominantly into a pair of dark matter in most of the parameter space satisfying the relic density shown in red solid, blue dashed and orange dotted lines for fermion, scalar and vector dark matter, respectively. So, the bound from invisible searches at BaBar applies to the whole parameter space below m G = 8 GeV, excluding the relic density region for scalar dark matter below m G = 0.8 GeV but less constraining the counterparts for fermion or vector dark matter. The future Belle-2 results [44,45] could improve the limits or probe the larger portion of the relic density regions. We also show the (g − 2) µ favored region in green and orange at 1σ and 2σ levels, respectively, but it is excluded by BaBar. In the same plot, we show the gray contours for DM-electron scattering cross section with σ DM−e = 10 −44 , 10 −48 cm 2 , but most of the parameter space survives the current direct detection bounds on light dark matter, such as XENON10, DarkSide-50, Sensei experiments. We note that as shown in the results, (4.33) and (4.34), the bounds from K + → π + + G or B + → K + + G with G → invisible are much weaker than BaBar invisible searches, so they are not shown in Fig. 8.
On the other hand, in the right plot of Fig. 8, as shown in Figs. 4 and 5, we don't need large graviton couplings to the SM particles in the region with m DM > m G , because dark matter can annihilate directly into a pair of spin-2 mediators. Therefore, the relic density can be determined almost independent of the graviton couplings to the SM particles, so a lot of parameter space for the correct relic density can be compatible with the current experiments. In this case, the spin-2 mediator decays only into the SM particles, so mono-photon + leptons at BaBar applies, limiting the lepton couplings to the spin-2 mediator. In the same plot, we also show the gray contours for DM-electron scattering cross section with σ DM−e = 10 −48 , 10 −52 cm 2 , so most of the parameter space is unconstrained by direct detection yet.

UV completions for spin-2 mediators
We can regard the spin-2 mediator as the first Kaluza-Klein(KK) mode of graviton from the warped extra dimension. In this case, there are heavier Kaluza-Klein(KK) modes of graviton, which can be summed up to modify the DM processes, such as DM annihilation and scattering.
In this section, we first summarize the KK graviton masses and couplings for two benchmark models with the warped extra dimension and discuss the effects of the heavier KK modes in determining the relic density, the direct detection bounds as well as the direct production of KK gravitons at colliders, in order.

Spin-2 mediator from the warped extra dimension
The KK modes of graviton in Randall-Sundrum(RS) model [20] are spaced almost equally. So, if dark matter is lighter than almost twice the mass of the first KK mode, the heavier KK modes would not change much our discussion with the first KK mode only. Otherwise, we need to include the heavy KK resonances explicitly. On the other hand, in the 5D continuum limit of the clockwork model, so called the linear dilaton model [21,[52][53][54], the KK modes of graviton are almost degenerate with a mass gap from the zero mode, challenging for experimental tests [55,56]. So, it is crucial to include the heavier KK modes in the DM processes in this case.
Suppose that m n are KK graviton masses, and c DM,n , c SM,n are the couplings of the nth KK mode to dark matter and the SM, respectively, and depending on the localization in the extra dimension. Here, dark matter and the SM particles can be localized on the IR brane, in which case dark matter has sizable couplings to the SM particles. But, when the SM particles are localized away from the IR brane, we can just rescale c SM,n to small values.
In the case where dark matter and the SM particles are localized on the IR brane, the KK graviton couplings and KK graviton masses are given by Here, for RS model, m G = x 1 k RS e −k RS πR with k RS being the AdS curvature scale, and x n are the zeros of J 1 (x n ) = 0, i.e. x n = 3.83, 7.02, 10.17, 13.32 for n = 1, 2, 3, 4, which can be approximated to x n = (n + 1/4) + O(n −1 ) for n 1, and R is the radius of the warped extra dimension. For CW model, m G = k CW with k CW being the 5D curvature scale. Moreover, the overall suppression scale for massive graviton couplings is where M P , M 5 are the 4D and 5D Planck masses, respectively, and the relations between them were used in the second equality in each line. Therefore, the KK graviton mass and the KK graviton coupling can be chosen independently, attributed to the choice of the 5D curvature scale (k RS or k CW ) and the radius of the extra dimension R. We note that the ratio of the first KK graviton mass to the suppression scale are given by m G The model dependence of the widths of heavier KK gravitons is discussed in appendix B. The effects of KK modes of graviton on dark matter physics were discussed in the context of the RS model [15] and the continuum clockwork model [58]. It would be also interesting to generalize the above discussion to the case with more general warped geometries [57].
In this section, we focus on the minimal graviton interactions to catch the qualitative effects of KK modes on dark matter annihilation and scattering. When there are multiple graviton interactions, for instance, quartic interactions between dark matter and graviton [58], and cubic graviton self-interactions [55], dark matter can also annihilate into a pair of different KK gravitons, for instance, DM DM → G n G m , if kinematically allowed, contributing to the total annihilation cross section. But, we postpone a complete analysis with higher order graviton vertices to a future work.

Dark matter annihilations
First, the KK modes contribute to the s-channels of dark matter annihilating into the SM particles by where A s is the resonance-independent factors in the cross section. Then, using eqs. (C.1) and (C.6) in appendix C, we get the modified s-channel cross sections of scalar dark matter annihilating into a pair of the SM fermions, whose masses are ignored, as follows:  Finally, for vector dark matter with m X m G , using eqs. (C.5) and (C.14), the annihilation cross sections for XX → G n G n are also summed up to be We comment on the contributions of the heavier KK gravitons for dark matter annihilations. For RS model, the contributions of dark matter annihilating into a pair of heavier KK gravitons are numerically small, as compared to the contribution with the first KK graviton only, although there is a sizable effect about 10% in the case of fermion dark matter due to the smaller inverse power dependence on the KK masses. On the other hand, for CW model, the higher KK mode contributions can be important for k CW πR 1, numerically, if k CW πR 7.2 for scalar and vector dark matter and k CW R 5.0 for fermion dark matter. Therefore, we can allow for a smaller coupling for the first KK graviton in the case for m DM m G , as compared to those shown in Fig. 1 for WIMP and in Figs. 4 and 5 for light dark matter. But, we should also remark that the KK sums must be cut off at a finite n for which m n = m 2 G + n 2 /R 2 = m DM , etc, so the actual contributions from the heavier KK modes are smaller.

Dark matter scatterings
The contributions of KK gravitons to the t-channels of DM-nucleon scattering and DM self-scattering cross sections are given, respectively, by where A t is the factor independent of the KK graviton propagator in the cross section, and SM stands for nucleon for WIMP dark matter or electron for light dark matter. Similarly, the KK modes contribute similarly to the t-channels of DM-electron scattering for direct detection and kinetic equilibrium, with a similar approximate KK graviton propagator for small momentum transfer. We note that in the case of DM self-scattering, the t-channel contributions are dominant in the Born limit, so the above discussion on the t-channels would be sufficient. First, for the DM-nucleus scattering in direct detection, using eqs. (C.3) and (C.10), we only have to replace the effective nucleon couplings in eqs. (4.12) and (4.13) by the sum of KK modes, as follows,  Moreover, the momentum relaxation rate for kinetic equilibrium in eq. (4.16) becomes Finally, for the DM self-scattering, the corresponding t-channel cross sections in the Born limit are also modified due to the KK modes, as follows, As a consequence, for RS model, the contributions of the heavier KK modes to the t-channel scattering cross sections for dark matter are about 3.4 larger than the one of the first KK mode only. For CW model, the contributions from the heavier KK modes depend on the warp factor, that is, they can be important for k CW R 1.1, independent of the spins of dark matter. Therefore, in both models, we can make the direct detection bounds less stringent on the couplings of the first KK graviton by including the heavier KK modes for the t-channel scattering processes.

KK graviton productions
Each of heavier KK modes of graviton can be also singly produced with a sufficiently large center of mass energy at LHC, with similar signatures as for the first KK graviton. However, in clockwork model, the KK graviton masses can be almost degenerate, namely, the mass difference between the n + 1-th and n-the KK graviton masses is given by ∆m n ≡ m n+1 − m n = m G (2n + 1)/(2(kR) 2 ) m G for kR 1. In this case, almost continuum KK gravitons can be produced simultaneously, leading to the photon or lepton energy spectrum of periodic shape [55,56].
As we discussed in Section 4.4, another smoking-gun signal for the spin-2 mediator would be through e + e − → γ G or qq → g G, which could identify the signatures of spin-2 mediator couplings. For s m 2 G , the heavier KK modes can be also produced at the LHC. In RS model, the KK graviton masses are well separated, so we could search for the heavier KK modes as for the first KK graviton as we discussed in Section 4. On the other hand, in clockwork model, almost continuum KK gravitons could be produced against mono-jet, decaying visibly or invisibly, so the resulting experimental signatures could be significantly different from those in the effective theory only with a single spin-2 mediator case.

Conclusions
We have explored the general production mechanisms for WIMP and sub-GeV scale light dark matter with arbitrary spin in the scenarios of gravity-mediated dark matter. The spin-2 mediator interactions of dark matter as well as SM particles are constrained by direct and direct detections, precision measurements and collider experiments. We showed that the parameter space where dark matter annihilates dominantly into the SM fermions is disfavored, due to direct detection and LHC dijet bounds for weak-scale WIMP case, and mono-photon searchers at BaBar experiments for light dark matter. On the other hand, we found that when dark matter annihilates dominantly into a pair of spin-2 particles in both allowed and forbidden regimes, the model is consistent with current bounds from direct detection and collider experiments. In particular, light dark matter with forbidden channels is not constrained by current indirect detection and CMB measurements.
A DM-nucleon scattering amplitudes

Scalar dark matter
From the results, we obtain the scattering amplitude between fermion dark matter and nucleon as follows, (2)) + c g G(2)

(A.2)
On the other hand, the partial decay rates of the invisible decays of the KK graviton [15] are also given by For RS model, the heavier KK modes of graviton couple to the SM particles with the same strength as for the one for the first KK graviton, so we only have to replace the graviton mass by those for the heavier KK modes in the above formulas. Thus, the narrow width approximation holds for the heavier KK modes.
For CW model, the couplings of the KK modes of graviton are level-dependent, such as c SM(DM),n = (k CW R) n/(m n R) for the SM(DM) particles localized on the IR brane. Thus, the partial decay widths of the KK gravitons scale by the overall factor. For instance, the decay rate of the nth KK graviton G n into a gluon pair becomes etc. The overall factor, n 2 mn m G , is approximated to n 2 for k CW R 1, so the partial decay widths of heavier KK gravitons get enhanced, as compared to the case in RS model with the same coupling for the lightest KK graviton.

Randall-Sundrum model
The KK sum relevant for the s-channels in RS model is in narrow width approximation where J 1 (x n ) = 0 and we used The KK sum relevant for the t/u-channels in RS model is given by where we used J n (σ)