Dark matter-neutrino interaction in light of collider and neutrino telescope data

We study the DM-neutrino interaction in the framework of simplified model. The phenomenology of such an interaction are derived. We also investigate the bound on DM-neutrino interaction from the LHC and neutrino telescopes. We find that for the case of a scalar dark matter, the LHC gives a stronger bound on dark matter annihilation cross-section than the neutrino telescopes. However, for the fermionic dark matter case the neutrino telescopes bounds are more stringent for dark matter mass, ≳ 200 MeV. In the case of lower DM mass, the neutrino telescopes provide better bounds for a light mediator, while the collider bounds are better for a heavy mediator. Possible UV completions of the simplified model are briefly discussed.


Introduction
The LHC has played an important role in probing the interaction between dark matter (DM) and the Standard Model (SM) particles. It provides an alternative to the direct and indirect DM searches in probing the DM parameter space. The advantage of the LHC search is its independency of the cosmological assumptions. Because DM is invisible to the LHC detectors, its LHC signature involves a visible object X recoiling against the invisible DM. The object X depends on the details of the DM interaction with the SM sector. One of the earlier considered scenario is the interaction of DM with the light quarks in the effective field theory (EFT) framework. In this scenario, the DM recoils against the initial state radiation (ISR) coming from the quarks. The strongest bounds comes from the ISR jet, hence monojet signature [1][2][3][4][5][6][7][8], while the monophoton signatures also provide a competitive bound [9][10][11][12]. It was soon realized that the EFT framework is not applicable at the LHC energy scale where the typical momentum transfer in DM-SM interaction is comparable to the suppression scale of the EFT. In order to avoid working with a specific model, simplified models [13][14][15][16][17][18][19][20][21][22][23][24][25][26][27][28][29][30] are introduced to correctly taking into account the kinematics of DM interactions. The simplified models introduce mediators between DM and the SM sector. This further enriches collider signatures of the dark sector. In some regions of parameter space, the search for the mediators provides a better bound than the search for DM alone.

JHEP06(2018)026
to be made to compare the LHC and the cosmological bounds. In the t-channel mediator framework, the DM-neutrino-mediator interaction is parametrized by a single coupling. Here no additional assumptions are needed to compare the cosmological and the LHC constraints. We will use this t-channel mediator assumption for the remainder of this paper.
We make an assumption that DM, χ, a mediator, φ (or ψ for a fermionic mediator), and the SM fields are the only light degrees of freedom. The rest of new physics particles are heavy and decoupled. In this simplified scenario, the renormalizable interaction between DM and a neutrino is given by yψP L νχ + h.c., fermionic DM, scalar DM, (2.1) where y is the DM-neutrino coupling constant. For simplicity we assume that y is real. We also assume that the coupling is flavor democratic. For χ to be stable, we must have m χ ≤ m φ(ψ) . 1 In principle, DM can be its own antiparticle if χ is a Majorana fermion or a real scalar. However, for definiteness, in this work we'll assume that DM is not its own antiparticle. The cosmological constraints for the scalar case has been considered in ref. [52].
In this setup, we'll allow for the possibility that χ and φ(ψ) are not their own antiparticles. Generically, both χ and φ(ψ) could be part of an electroweak multiplet. However, in this work we'll leave their electroweak charge assignment unspecified except that they're neutral under U (1) EM . We'll briefly revisit the topic of their possible charge assignments and UV completions in section 5. In the rest of this section, we explore the phenomenological consequences of the DMneutrino interaction in eq. (2.1).

Neutrinos signal from DM annihilation
In our simplified scenario of DM-neutrino interaction given in eq. (2.1), a pair of DM particle can annihilate into a pair of neutrinos, see figure 1. The DM self annihilation cross-section reads fermionic DM, scalar DM, where v rel is the DM relative velocity. Note in the above expression we have summed over the 3 neutrino species. The neutrinos produced by DM annihilation will each carries energy m χ . They can be looked for at neutrino telescopes such as IceCube, KamLAND and Super-Kamiokande. The flux of a particular neutrino (or an antineutrino) flavor from the DM self-annihilation JHEP06(2018)026 observed on Earth can be written as [53] where f is a lepton family index, R sc = 8.5 kpc is the Solar radius [54] and ρ sc is the DM density at the Solar radius. The factor of 1/3 comes from the assumption that each neutrino flavor is equally likely to be produced. Note in the case where DM is its own antiparticle, the factor 1/4 in the above equation becomes 1/2. The neutrinos from the annihilation point will oscillate in flight on the way to the detector. With the present neutrino oscillation parameters, we expect the flavor ratio of the neutrino at the detection point to be ν e : ν µ : ν τ 1 : 1 : 1 [55]. The factor J avg represent the average line-of-sight integral over the whole sky where θ is the angle subtended from the DM annihilation point to the Galactic Center, ρ(r) is the DM density profile as a function of distance from the Galactic Center and l is the distance from the annihilation point to the Earth. The limit of integration depends on the size of the DM halo and is given by l max = R 2 halo − R 2 sc sin 2 θ + R sc cos θ. In our calculation, we take R halo = 100 kpc.
The dark matter profile can be parametrized as where ρ 0 is the density normalization deduced from the DM density at the Solar radius (ρ sc ), r s is the scale factor and α, β, γ, δ are the shape parameters. In this work, we choose the Navaro, Frank and White (NFW) [57] and the Burkert [58] profiles as representative DM profiles. The NFW profile represents a cuspy profile while the Burkert profile represents profile with a core. The value of the parameters for both the NFW and the Burkert profiles are given in  Table 1. Parameters for the NFW and the Burkert DM density profiles [56].

DM-neutrino scattering
The DM-neutrino interaction in eq. (2.1) induce DM-neutrino scattering via Feynman diagram shown in figure 2. For the fermionic DM case, the DM-neutrino cross-section is given by where we have sum over the 3 neutrino species in the final state. The cross-section for the scalar DM case is the same up to a trivial replacement m φ → m ψ . This scattering process leaves an imprint on CMB and matter power spectra which affects the large-scale structure of the universe. The most stringent cosmological constraint on DM-neutrino interaction comes from the matter power spectrum, σ(χν → χν) < 10 −45 mχ GeV cm 2 [49]. It translates to the bound on the DM-neutrino coupling as a function of the DM and the messenger mass. Figure 3 shows such a bound for the benchmark case of messenger mass of 0.01, 0.05 and 0.1 GeV respectively. Since the χ-ν scattering cross-section scales roughly as 1/m 4 φ(ψ) , the bound on the coupling becomes weak at large messenger mass as can be seen from the plot. Similarly, one can translate the bound on χ-ν cross-section into the bound on DM annihilation cross-section, see figure 9 and 10.

New decay channels for Z and W
The DM-neutrino interaction in eq.  detector. Hence, the Z → νχφ decay contributes to the invisible decay width of the Z, 2 while the W → χφ contributes to the mono-lepton decay width of the W . The differential partial decay width for Z → νχφ is given by , v ew is the electroweak vev and m Z is the Z boson mass. The limits on the kinematic variables are The current invisible branching ratio of the Z boson is (20.00 ± 0.06)% [60]. The differential partial decay rate for W → χφ is given by where m W is the W boson mass and (2.13) The variables p φ , in this case, are given by 14) The current mono-lepton branching ratio of the W boson is (32.72 ± 0.30)% [60]. For the case of scalar DM, the expression for the differential partial decay width of the Z and the W boson can be obtained by an obvious substitution φ → ψ.
Moreover, for a hadron collider such as the LHC, the Z → νχφ(ψ) and W → χφ(ψ) decay channels inevitably lead to the mono-jet signal and the mono-lepton signal respectively. These collider signatures are actively being investigated at the LHC. We'll discuss the collider bound in details in the next section.

LHC bounds
At the LHC, the neutrinos are produced copiously from the W and Z decays. If the neutrinos interact with the dark sector according to eq. (2.1), the DM and the mediator will also be produced at the LHC, as shown in figure 4. While both the DM and the mediator are invisible at the LHC, this interaction could change the kinematics of W and Z productions. In the rest of this section we will assume the fermionic DM case for definiteness. The scalar DM case can be obtained with a trivial replacement φ → ψ for the mediator.
The LHC signature for the dark sector production via W -boson is a lepton and missing energy from φ and χ, dubbed monolepton. The cuts imposed by the ATLAS or CMS force the W to be off-shell. Hence the distribution of the lepton momentum can be approximated as dσ pp→ φχ dp ∝ p s m 2 where p is the lepton momentum, s is the center-of-mass energy and m φχ is the invariant mass of the φ and χ. The background for this process is the production of lepton and neutrino. In this case, the lepton momentum distribution is Taking into account the parton distribution function, the lepton momentum distributions are shown in figure 5. From the plot, one can see that the lepton momentum for the signal tends to be harder than the background. If φ and χ are produced via Z-boson, they need an object to recoil against to show a LHC signature. The most promising search channel for this scenario is the monojet channel. However we found that this channel is not as sensitive as the monolepton channel. The monojet channel is a dirtier channel for the LHC environment. Hence commonly ATLAS and CMS puts strong cuts on the jet p T hence the signal acceptance is reduced. For this paper, we only consider monolepton channel. Combining the monolepton channel with monojet channel does not change the results appreciably.
We base our analysis on the 13 TeV, 36.1 fb −1 ATLAS monolepton search [45]. In the electron channel, the ATLAS analysis requires exactly one electron with p T > 65 GeV and |η| < 1.37 or 1.52 < |η| < 2.47. The requirement for the muon channel is exactly one muon with p T > 55 GeV with |η| < 2.5. The missing energy has to be greater than 65 (55) GeV in the electron(muon) channel. In both channels, the events are vetoed if there is an additional lepton with p T > 20 GeV. The signal and the backgrounds are discriminated by the transverse mass distributions. The transverse mass variable is defined   Both of the background and the signal events are simulated using MadGraph 5 [61], followed by matching, parton shower and hadronization using Pythia 6 [62]. The detector simulation and pileup are simulated using Delphes 3 [63]. The average number of pileup in the simulation is 23. The signal model is generated using Feynrules 2.0 [64]. The comparison for the ATLAS and our simulation for the background process W → ν are shown in figure 6. The signal transverse mass distribution compared with the total background is given in figure 7.
ATLAS collaboration bin their results according to the transverse mass. There are seven bins, starting from m T =130(110) GeV for electron(muon) channel to 7 TeV. To get the bounds, we compared the signal and the background in each bin. The likelihood of observing the signal in a bin is given by   where b i is the number of predicted background in the bin, n i is the number of observed events in the bin, s i is the number of predicted signal in the bin which is a function of the coupling y. The function P i (ξ) is introduced to incorporate the uncertainty in each bin and take a log-normal form [65] where σ i is the systematics error for the corresponding bin. The chi-squared value is given by We combine the bounds from both of the electron and muon channels. The 90% confidence level (CL) bounds are shown in figure 8. There, we consider the bounds for both the fermionic and the scalar dark matter cases with various dark matter masses. For the benchmark points, the mediator mass are chosen to be 0.1 GeV, 1 GeV, 50 GeV and 100 GeV. We also show bounds from the invisible Z decay width. From the plot, one can see that the LHC bounds become more competitive for a higher mediator mass, where the Z decay is starting to become kinematically non-favorable. These bounds can be compared with the indirect detection bounds discussed in the next section.

Indirect detection bounds
The DM-neutrino interaction gives rise to a pair annihilation of DM particles into a pair of neutrinos. These neutrinos, in principle, can be observed at the neutrino telescopes on Earth. For example, the IceCube experiment searches for neutrinos from DM annihilation in the Milky Way for DM mass range between 30 GeV and 10 TeV [66], with the updated IceCube analysis covers DM mass range between 10 GeV and 1 TeV [67].

JHEP06(2018)026
These two IceCube studies place an upper limit on the DM self annihilation cross-section, σv O(10 −23 ) cm 3 /s. Similar analysis has been carried out by the Super-Kamiokande (SK) experiment for DM mass ranging from 1 GeV to 10 TeV [68], resulting in the upper bound σv O(10 −24 − 10 −23 ) cm 3 /s. Additionally, the KamLAND experiment has searched for neutrinos from DM annihilation for a low DM mass region from 8.3 MeV to 31.8 MeV and places an upper limit of σv O(10 −24 ) cm 3 /s [69]. In additional to the dedicated neutrinos from DM annihilation searches, one can also derive the upper bound on DM annihilation cross-section from other existing neutrino searches. In particular, we will focus on reinterpreting the SK supernova relic neutrino search [70] and the SK atmospheric neutrino measurement [71] in the context of DM annihilation.
Additionally, neutrino telescopes data could be used to bound the DM-neutrino crosssection, which translates to the constraint on DM-neutrino coupling. Ref. [72] obtained an upper bound on DM-neutrino coupling from analysing the IceCube high-energy neutrino data set [73]. Such a bound is relevant for high DM mass and/or high messenger mass and can be complimentary to the bounds obtained in this work.

Bounds from Super-Kamiokande supernova relic neutrino searches
A supernova explosion throughout the history of the universe release most of its energy in the form of neutrinos. These neutrinos, referred to as the supernova relic neutrino (SRN), should remain in existence today. The SRN signal can be detected at the SK detector via the inverse beta decay process, whereν e scatters off a free proton inside the SK detector into a neutron and a positron. These SRN data can be reinterpreted to derive a constraint on DM annihilation cross-section for DM mass in the range of 10 ∼ 130 MeV. Ref. [55] carried out such an analysis on the first 1497 days of SK data (SK-I) [74]. In this subsection, we give the updated analysis by including data from the second (SK-II) and the third (SK-III) data taking period of the SK detector [70]. Our analysis differs from that of ref. [55] in the treatment of the SK energy resolution and signal efficiency which we describe below. Also in the case of a fermionic DM, we take DM to be a Dirac particle instead of a Majorana fermion. As a result, our bounds come out to be similar to the one obtained in ref. [55] even though we use roughly twice as much data.
The SK SRN searches detect e + from the inverse beta decay process in the energy range 18-82 MeV. The detected e + is then binned into 16 four-MeV energy bins. To take into account the energy resolution of the SK detector, we introduce the energy resolution function where E is the actual e + energy, E is the detected energy and σ is the SK detector resolution which is given in [75]. 3

JHEP06(2018)026
The expected number of e + from DM self-annihilation in the energy bin i is given by where N SK = 1.5 × 10 33 is the number of proton in the SK fiducial volume, T is the total running time, is the minimum (maximum) actual e + energy, E is the detector efficiency which is a function of positron energy, and σν ep is the inverse beta decay cross-section. In our computation, we take σν ep from ref. [76]. In order to derive the upper limit on DM self annihilation cross-section from the SK data, we introduce the likelihood function for each energy bin as in eq. (3.3). Then we construct the χ 2 function, as in eq. (3.5). Note in this case the sum is over all energy bins and over the three SK data taking periods. Finally, the 90% CL upper limit on DM self-annihilation cross-section is determined from ∆χ 2 = 2.71.

Bounds from Super-Kamiokande atmospheric neutrino measurements
The SK experiment measures the atmospheric neutrino fluxes in the energy range from 100 MeV to 100 GeV. The SK detector is sensitive to both the electron-and the muonneutrinos. However, the detector cannot distinguish neutrinos from antineutrinos, thus SK reports the fluxes of ν e +ν e and ν µ +ν µ separately. The latest SK atmospheric neutrino measurements, which include SK I-IV data set, is reported in ref. [71]. In these measurements, the observed neutrino fluxes as well as their uncertainties are reported.
The SK atmospheric neutrino measurements can be reinterpret to derive a bound on DM annihilation cross-section for DM mass between 100 MeV and 100 GeV. We use the method of least square to determine the bound. The χ 2 function is defined as where Φ ν e(µ) ,i atm is the predicted atmospheric electron (muon) neutrino flux for bin i, Φ ν e(µ) ,i obs is the observed atmospheric electron (muon) neutrino flux at bin i and σ ν e(µ) ,i is the uncertainty associated with bin i. This includes both the uncertainty in the measurement (statistical and systematics) as well as the uncertainty in the predicted atmospheric neutrino flux. The differential neutrinos flux per solid unit angle from DM annihilation in bin i, dΦ ν e(µ) ,i DM /dΩ, according to eq. (2.3) can be is written as where E min and E max is the minimum and maximum energy of each bin respectively and the factor of 2 come from combining both the neutrinos and the antineutrinos. Note that in the case where DM is its own antiparticle, the factor of 1/4 in the cross-section becomes 1/2. We take the value of the predicted atmospheric neutrino flux from the HKKM simulation [77]. We estimate the uncertainty in the predicted atmospheric neutrino flux by comparing it against the Bartol [78] and FLUKA [79] simulations.   Figure 9. The 90% CL bounds from various neutrino telescope observations compared with the collider bounds. The gray line, labelled KamLAND, is the (rescaled) official KamLAND bound [69]. The green line, labelled SK, is the official SK bound [68]. The brown line, labelled IceCube, is the official IceCube bound [67]. The orange, dot-dashed line, labelled SK-SRN, is the bound obtained from recasting the SK SRN searches discussed in section 4.1. The magenta, dot-dashed line, labelled SK-νatm, is the bound obtained from recasting the SK atmospheric neutrino measurements discussed in section 4.2. The shaded yellow region is excluded by perturbativity constraint on the DM-neutrino coupling, y, and the constraint m χ ≤ m φ(ψ) . The shaded cyan region is the region that gives Ω χ h 2 > 0.1186 simulated using Micromegas 5.0 [80]. The blue dotted lines are the cosmological bounds from the DM-neutrino scattering discussed in section 2.2. For the collider bounds (dashed lines), we only show strongest bound between the LHC and the invisible Z decay width bounds. We subject the collider bounds to y ≤ 4π and m χ ≤ m φ(ψ) . In these plots, the dark matter profile is taken to be the NFW profile.

Comparison of bounds
In this section we present all bounds obtained from the previous two subsections. Figures 9  and 10 show the bounds for both the fermionic and scalar DM cases in the m χ -σv plane assuming the NFW and Burkert DM profile respectively. Figures 11 and 12 display the same bounds in the m χ -m φ(ψ) plane. The bounds are compared against those obtained from the LHC and Z decay width measurement in the context of simplified model of section 2. Additionally, we also show the official bounds from IceCube [67], Super-Kamiokande [68] and KamLAND [69] collaborations. Instead of picking a specific DM profile, the KamLAND collaboration uses a nominal value of J avg = 5. Thus, we rescale the KamLAND bounds using an appropriate J avg factor of each DM profile. On the other hand, the Super-Kamiokande collaboration did the analysis for the NFW profile only. Hence, for the Burkert profile, we rescale the official SK bound with the appropriate J avg factor.
From the plots we can see that in the case of the fermionic DM, the official bounds from SK and IceCube together with the bound from recast SK atmospheric neutrino analysis, which covers the mass region m χ 200 MeV, are stronger than the collider bounds. In the case of low dark matter mass with light mediator, the SK and KamLAND bounds are more constraining than the collider bounds. However, for low DM mass and heavy mediator, colliders provide stronger bounds.  Figure 10. The 90% CL bounds from various neutrino telescope observations compared with the collider bounds. The gray line, labelled KamLAND, is the (rescaled) official KamLAND bound [69]. The green line, labelled SK, is the (rescaled) official SK bound [68]. The brown line, labelled IceCube, is the official IceCube bound [67]. The orange, dot-dashed line, labelled SK-SRN, is the bound obtained from recasting the SK SRN searches discussed in section 4.1. The magenta, dotdashed line, labelled SK-νatm, is the bound obtained from recasting the SK atmospheric neutrino measurements discussed in section 4.2. The shaded yellow region is excluded by perturbativity constraint on the DM-neutrino coupling, y, and the constraint m χ ≤ m φ(ψ) . The shaded cyan region is the region that gives Ω χ h 2 > 0.1186 simulated using Micromegas 5.0 [80]. The blue dotted lines are the cosmological bounds from the DM-neutrino scattering discussed in section 2.2. For the collider bounds (dashed lines), we only show strongest bound between the LHC and the invisible Z decay width bounds. We subject the collider bounds to y ≤ 4π and m χ ≤ m φ(ψ) . In these plots, the dark matter profile is taken to be the Burkert profile.
For the case of scalar DM, its annihilation cross-section is p-wave suppressed. As a result, the collider bounds are more constraining than the neutrino telescope bounds for a low relative DM velocity. In figures 9b, 10b, and 12 we show the bounds in the case of v rel = 10 −3 . For different value of v rel , one can rescale the collider bound with a factor of v rel 10 −3 2 . Note in particular that the collider bound in figures 9b and 10b for the case of m ψ = 50 GeV is slightly stronger than the case of m ψ = 100GeV for low DM mass. This is because of the invisible Z decay bound gives a stronger constraint than a LHC bound.

Comments on the mediator and possible UV completions
So far we have remained agnostic about the UV completion of the interaction in eq. (2.1). In this section, we will briefly discuss the role of the mediator, in a UV complete model.
It is possible to embed the interaction in eq. (2.1) in a UV complete model such that the mediator is an SM singlet. In this case, the DM (χ), must be part of an electroweak doublet. For a concrete example, consider the supersymmetric Standard Model with righthanded neutrinos. We could have the DM being the neutralino while the mediator (φ) is the right-handed sneutrino.
Another   Figure 11. The 90% CL bounds from various neutrino telescope observations compared with the collider bounds in the m χ -m φ plane. The gray line, labelled KamLAND, is the (rescaled) official KamLAND bound [69]. The green line, labelled SK, is the official SK bound [68]. The brown line, labelled IceCube, is the official IceCube bound [67]. The orange, dot-dashed line, labelled SK-SRN, is the bound obtained from recasting the SK SRN searches discussed in section 4.1. The magenta, dot-dashed line, labelled SK-νatm, is the bound obtained from recasting the SK atmospheric neutrino measurements discussed in section 4.2. The blue dotted lines are the cosmological bounds from the DM-neutrino scattering discussed in section 2.2. The dashed lines are the bounds obtained by taking the strongest bounds between the LHC mono-lepton and the invisible Z decay width. All the regions below the lines are excluded by the respective observations/experiments. The shaded cyan region is the region that gives Ω χ h 2 > 0.1186 as simulated by Micromegas 5.0 [80]. The shaded yellow region is the region which m χ > m φ . In these plots, the dark matter profile is taken to be the NFW profile.
case of a fermionic DM, the scalar mediator resembles the a slepton of the supersymmetric Standard Model. Hence this scenario is constrained by the LHC slepton searches. The latest LHC result puts the lower bound on the mass of charged partner of the mediator (m φ ± ) to be around 500 GeV for m φ ± − m χ 60 GeV [81]. Since we consider the case of m φ 100 GeV, to satisfy this bound, the mass of the neutral and charged component of the doublet must be split. This could be achieved by introducing an interaction with the Higgs doublet such as L ⊃ λ HΦ (H † Φ)(Φ † H), where Φ is the electroweak doublet containing the mediator. Alternatively, one can avoid the bound on the mass of the charged partner of the mediator in the case of the compressed spectra, m φ ± ∼ m φ m χ + 60 GeV. In this case the leptons coming from the decay of the charged partner will be too soft to pass the search cuts.
We could instead have a fermionic mediator being part of an electroweak doublet, Ψ. The mediator must be vector-like to avoid anomalies, ie.
L ⊃ yΨ R L L χ + MΨ R Ψ L + h.c., (5.1) where L L is the left-handed lepton doublet. Ψ L and Ψ R are the left-and right-handed fermion doublets with the same charge assignment as the SM left-handed lepton doublet. As in the scalar case, there exists a bound from the charged partner of the mediator. The current LHC bounds on a vector-like lepton is m ψ > 176 GeV [82].    [69]. The green line, labelled SK, is the official SK bound [68]. The brown line, labelled IceCube, is the official IceCube bound [67]. The orange, dot-dashed line, labelled SK-SRN, is the bound obtained from recasting the SK SRN searches discussed in section 4.1. The magenta, dot-dashed line, labelled SK-νatm, is the bound obtained from recasting the SK atmospheric neutrino measurements discussed in section 4.2. The blue dotted lines are the cosmological bounds from the DM-neutrino scattering discussed in section 2.2. The dashed lines are the bounds obtained by taking the strongest bounds between the LHC mono-lepton and the invisible Z decay width. All the regions below the lines are excluded by the respective observations/experiments. The shaded cyan region is the region that gives Ω χ h 2 > 0.1186 as simulated by Micromegas 5.0 [80]. The shaded yellow region is the region which m χ > m ψ . In these plots, the dark matter profile is taken to be the NFW profile.
Lastly, both the DM and the mediator could be singlets under the SM gauge groups. The dark sector communicates with the SM sector via mixing with a right handed neutrino of the see-saw mechanism. Several bounds regarding the right handed neutrino in this scenario are discussed in [83,84].

Conclusions and discussions
DM interactions are among one of the least understood interactions in nature. There is hope that the LHC can provide insight into DM interactions.
In this work, we consider the case where DM only interact with the SM through a neutrino and a mediator. We study the interaction in the framework of the simplified model, see eq. (2.1). In this scenario, there are two viable options to probe the DM interaction. First, one can probe such an interaction using the mono-lepton and the monojet signal at the LHC. These experimental signatures arise from the W and Z production. The W and Z decays, via the off-shell neutrino, to DM and the mediator as shown in figure 4. By reinterpreting ATLAS mono-lepton result [45] as well as the invisible Z width measurement [60], we derive the bounds on DM-neutrino coupling shown in figure 8.
Moreover, one can probe the DM-neutrino interaction using the neutrino signal coming from DM annihilation in the galaxy. Here, we consider both the official bounds from a JHEP06(2018)026 dedicated neutrino from DM annihilation search at IceCube [66,67], KamLAND [69] and SK [68], as well as our reinterpretation of SK data from SRN search, see section 4.1, and the atmospheric neutrino measurement, see section 4.2.
We found that in the case of fermionic DM, neutrino telescopes give the best constraint for m χ 200 MeV. For lower DM mass, neutrino telescopes continue to provide the best constraint for a light mediator. However, for low DM mass with a heavy mediator, the LHC provides the better bounds, see figure 9a and 10a. For the case of scalar DM, we found that the LHC bounds are generally more stringent than the bound from the neutrino telescope for all DM mass considered in this work, i.e. 1 MeV < m χ < 100 GeV, see figure 9b and 10b. This can be understood by noting that the DM annihilation cross-section for the scalar DM is velocity suppressed as shown in eq. (2.2). Hence the LHC bound on DM coupling, when interpreted in term of the bound on DM annihilation cross-section, get suppressed by a small relative DM velocity.