#
Tum-Hep 912/13
High-energy neutrino signals from the Sun

in dark matter scenarios

with internal bremsstrahlung

###### Abstract

We investigate the prospects to observe a high energy neutrino signal from dark matter annihilations in the Sun in scenarios where the dark matter is a Majorana fermion that couples to a quark and a colored scalar via a Yukawa coupling. In this minimal scenario, the dark matter capture and annihilation in the Sun can be studied in a single framework. We find that, for small and moderate mass splitting between the dark matter and the colored scalar, the two-to-three annihilation plays a central role in the calculation of the number of captured dark matter particles. On the other hand, the two-to-three annihilation into gives, despite its small branching fraction, the largest contribution to the neutrino flux at the Earth at the highest energies. We calculate the limits on the model parameters using IceCube observations of the Sun and we discuss their interplay with the requirement of equilibrium of captures and annihilations in the Sun and with the requirement of thermal dark matter production. We also compare the limits from IceCube to the limits from direct detection, antiproton measurements and collider searches.

## 1 Introduction

The search for the particles which are presumed to be produced by dark matter annihilations is one of the primary goals in Astroparticle Physics (see [1, 2] for reviews on particle dark matter). This search is however hindered by the low fluxes expected from annihilations and by the existence of large, and still poorly understood, astrophysical backgrounds. Therefore, in order to detect a signal from annihilations, it is necessary to devise search strategies that mildly rely on the modelling of the backgrounds. Up to now, the most promising strategies to identify a signal from dark matter annihilations consist in i) the search for features in the measured energy spectrum of the products of annihilation and ii) the search for an excess of events following a spatial distribution with the morphology expected from an annihilation signal.

Both strategies have been successfully applied to the search for annihilation signals in gamma-rays [3]. Several spectral features have been identified [4, 5, 6, 7, 8, 9, 10] and searched for in the gamma-ray sky [11, 12, 13, 14]. Besides, the characteristic morphology of the signal expected from dark matter annihilations has been searched for in the Galactic Center region [15, 16], galaxy clusters [17, 18] and, notably, dwarf galaxies [19, 20, 21, 22]. Both search strategies already lead to very stringent limits on the annihilation cross section which rule out some models. Pursuing these two promising strategies is, however, challenging when searching for antimatter particles produced in dark matter annihilations, since the propagation of the charged cosmic rays in the galaxy practically erases any spectral feature at production and any directional information, yielding a rather featureless spectrum and a fairly isotropic distribution of events. As a result, the limits from antimatter searches are typically weaker than in gamma-rays (see, however, the recent analyses [23, 24]). Lastly, the search for dark matter annihilation using neutrinos as messengers shares many similarities with gamma-rays: both are electrically neutral particles and therefore the propagation in the galaxy does not significantly distort the energy spectrum nor the direction of these particles. Nevertheless, searches with neutrinos as messengers are limited by the small neutrino detection rate and by the poor energy and angular resolution of neutrino telescopes, thus making the search for features in the neutrino spectrum challenging. Indeed, the current limits on the annihilation cross section for annihilations into [25, 26] are significantly weaker than the limits for annihilations into . Interestingly, if dark matter particles interact with nucleons, they can be captured inside the Sun [27] and subsequently annihilate, producing a high-energy neutrino flux in the direction of the center of the Sun [28, 29]. The very characteristic spatial morphology of the signal then opens the possibility of probing the scattering cross section of dark matter particles with nucleons. Conversely, the observation of an excess of high-energy neutrinos in the (time-dependent) direction of the Sun would be a clear signal of dark matter annihilations. The same phenomenon of capture might also occur in the Earth and other celestial bodies.

In this paper we will concentrate in the possibility of detecting a high-energy neutrino flux originated in the annihilations of dark matter particles that have been previously captured inside the Sun or inside the Earth. Most analyses postulate that the dark matter particles had been captured and that their population is constant today, due to the equilibration between the processes of capture and annihilation. In this way, it is possible to set strong limits on the capture cross section from the non-observation of a neutrino flux originated in annihilations into, for example, or [30]. However, this requires that the dark matter particle must couple to neutrinos or to weak gauge bosons, but also sizably to quarks to allow their capture. This scenario can be justified in concrete realizations although does not necessarily hold in general.

We will then discuss the capture and the annihilation in a single Particle Physics framework. More concretely, we will consider a toy model where the dark matter particle is a Majorana fermion that couples to a quark and a scalar via a Yukawa coupling. The coupling to the quark leads to the capture of dark matter particles inside the Sun or the Earth. On the other hand, and as is well known, in this toy model the annihilation rate into a quark-antiquark pair is helicity- and velocity- suppressed. Hence, in the equilibration the internal bremsstrahlung processes , with a gauge boson, play a major role and, in fact, determine the minimal value of the coupling that allows equilibration inside the Sun or the Earth. Moreover, these processes determine the energy spectrum of neutrinos produced in the annihilation and accordingly the detection rate at the Earth (neutrino signals from two-to-three annihilations in leptophilic dark matter models have been considered previously in [31, 32]). The parameter space of this toy model is further constrained by the requirement of reproducing the correct dark matter abundance from thermal production, as well as from direct, indirect and collider dark matter searches. Concretely, we will consider, following [33, 34], the limits from the XENON100 experiment [35], the PAMELA measurements on the antiproton-to-proton fraction [36], as well as collider limits from the CMS experiment [37]. All these limits, together with the requirement of building up a population of dark matter particles inside the Sun or the Earth, will then allow us to assess the prospects to observe a high energy neutrino signal from the Sun or the Earth in neutrino telescopes, concretely in IceCube.

The paper is organized as follows. In Section 2 we present the toy model under consideration as well as the properties relevant for our analysis. Then, in Section 3 we describe the process of capture of dark matter particles in the Sun and in the Earth and we present the constraints on the parameter space of our toy model from the requirement of equilibration between capture and annihilation. In Section 4 we calculate the neutrino flux from the Sun, emphasizing the role of the internal bremsstrahlung process in generating a high-energy component. In Section 5 we present the limits on the parameter space which stem from the negative searches of a high-energy neutrino excess in the direction of the Sun in IceCube and we compare our limits to those from other experiments. Lastly, in Section 6 we present our conclusions.

## 2 Toy model with internal bremsstrahlung

We consider a toy model where the dark matter particle is a Majorana fermion , singlet under the Standard Model gauge group, which interacts with a quark and a scalar via a Yukawa interaction with coupling constant . The Lagrangian of the model is given by:

(1) |

Here, is the Standard Model Lagrangian, while and are the parts of the Lagrangian involving just the new fields and and which are given, respectively, by:

(2) |

Here denotes the usual covariant derivative. Lastly, the interaction term in the Lagrangian is given by

(3) |

where is the Standard Model Higgs boson. Here we have assumed for concreteness that the quark is right-handed, hence the quantum numbers of the scalar under are , where is the right-handed quark electric charge (which in this case coincides with minus the hypercharge).

This toy model has been thoroughly discussed in relation to indirect dark matter searches. The lowest order annihilation process is helicity- and velocity- suppressed, resulting in an annihilation cross section which is far below the reach of present and foreseeable indirect dark matter search experiments. Interestingly, the helicity suppression is lifted by the associated emission of a spin-one boson [7, 8].^{1}^{1}1The emission of a Higgs boson also lifts the helicity suppression, as recently pointed out in [38]. The cross section, however, strongly depends on the parameters of the scalar potential, therefore we will neglect these processes in our analysis for simplicity. More specifically, the annihilation processes , with being a gluon [8, 39, 40], a -boson [41, 42, 43, 44, 45, 46, 47] (see also [48, 49]) or a photon [7, 8, 9] can dominate over the two-to-two process, provided the scalar that mediates the interaction is close in mass to the dark matter particle [42]. Moreover, in this regime a sharp spectral feature arises in the gamma-ray spectrum which, if observed, would constitute a smoking gun for dark matter detection. ^{2}^{2}2Sharp spectral features from internal bremsstrahlung also arise in the annihilation of scalar dark matter particles [50, 51, 52]. Therefore, new opportunities open for indirect dark matter searches using antimatter, gamma-rays and, as we will discuss in this paper, neutrinos.

The relative importance of the various annihilation channels is illustrated in Fig. 1, where we show the branching ratios as a function of the dark matter mass for the exemplary cases , namely a scenario with large degeneracy, and . In these plots, and in the rest of the paper, we neglect the quark masses in the calculation of the rates of the two-to-three processes. It is worth recalling that when , which is the range of masses relevant for neutrino telescopes, the ratios of cross sections of the two-to-three processes are fairly independent of the mass splitting and take the values [41]:

(4) |

for couplings to right-handed up-type quarks and

(5) |

for right-handed down-type quarks. The numerical values in these formulas have been obtained evaluating, for illustration, the strong coupling constant at the scale .

The dark matter interaction to a quark in Eq. (3) not only induces dark matter annihilations but also the elastic scattering of dark matter particles with a nucleus in a detector, thus opening the possibility of observing signals also in direct dark matter search experiments and, as we will see in Section 3, of capturing dark matter particles inside the Sun or the Earth. The scattering cross sections for this model were calculated in [53, 54, 55]. The spin dependent part reads

(6) |

where is the reduced mass of the proton-dark matter system and parametrizes the quark spin content of the proton for the light quark to which the dark matter couples; for heavy quarks, . Besides, the spin independent cross section with protons reads:

(7) |

where

(8) |

with an analogous expression for the neutron-dark matter cross section (note that we do not assume isospin invariance). Here, is a function of the masses of the respective quark, the dark matter particle and the scalar mediator which we take from [54, 56]. Moreover, () parametrizes the scalar quark (gluon) current in the proton/neutron; we use the values and uncertainties from [55]. Lastly, are the second moments of the particle distribution function of the quark inside the nucleon, given in [53]. Note that both the spin independent and the spin dependent cross sections get enhanced when the dark matter particle and the intermediate scalar are close in mass, which is precisely the region of the parameter space relevant for indirect dark matter searches in this scenario [33]. The expressions for the cross sections, Eqs. (6,7), on the other hand, become divergent when and therefore it is necessary to include close to the resonance the width of the exchanged scalar. However, and in order to avoid the model dependence introduced by the modelling of the scalar decay, we will restrict our analysis to values of the model parameters sufficiently away from the resonance, concretely we will only consider .

## 3 Dark matter capture in the Sun and in the Earth

In our toy model we have postulated that the dark matter particle couples to a quark, therefore, a dark matter particle traversing the Sun could scatter-off a nucleus in the solar interior and lose energy. After subsequent scatterings, the dark matter particles eventually sink to the solar core where they accumulate. At the same time, the dark matter population in the solar core is depleted by annihilations and – in case of light dark matter particles – by thermal evaporation. The time evolution of the number of dark matter particles in the solar core is then described by the following differential equation [57]:

(9) |

where is the capture rate, the annihilation constant and the evaporation constant.

For dark matter masses above the evaporation term can be safely neglected [56, 58]. Then, after solving Eq.(9), one finds an annihilation rate as a function of time given by

(10) |

where

(11) |

The annihilation rate reaches a maximum when . In this regime, captures and annihilations are in equilibrium and the annihilation rate is determined only by the capture rate: . It is important to note that, despite the Earth and the Sun have been capturing dark matter particles since their formation approximately ago, for some choices of the model parameters equilibrium might have not been reached. Whether dark matter particles are or not in equilibrium in the Sun or the Earth is determined by the values of the capture rate and the annihilation constant, cf. Eq. (11). The capture rate can be calculated from the scattering cross section of the process dark matter-nucleon and relies on assumptions on the density and velocity distributions of dark matter particles in the Solar System, as well as on the composition and density distribution of the interior of the Sun or the Earth. In our analysis we will use DarkSUSY to evaluate the capture rate from the scattering cross sections in Eqs. (6,7), adopting for concreteness a local dark matter density [59] and a homogeneous Maxwell-Boltzmann dark matter distribution with velocity dispersion [60]. On the other hand, the annihilation constant for the Sun is approximately given by [56]

(12) |

while for the Earth by

(13) |

where the annihilation cross sections for are given in [41].

Both the capture rate and the annihilation constant depend on the fourth power of the Yukawa coupling , thus allowing to define a lower limit on the coupling constant from requiring equilibration. Taking for concreteness the case of the Sun, the equilibration condition reads , where . More specifically, we define the equilibrium coupling constant from the requirement that the argument of the hyperbolic tangent in Eq.(10) is equal to 1 for , namely

(14) |

where and are, respectively, the capture and annihilation constants for . This implies that the ratio of annihilation rate and capture rate can be cast as:

(15) |

It follows from this equation that for values of the coupling larger than the equilibrium value, the annihilation rate, and hence the signal strength, reaches exponentially fast the maximal value. On the other hand, when the coupling is smaller than the equilibrium value, the annihilation rate exponentially decreases with . For example, for , the signal strength drops to about 0.4% of the maximal (equilibrium) signal strength.

We show in Fig. 2 as red solid lines the minimal values of the coupling constant which are required to achieve equilibration between capture and annihilation for dark matter particles that couple to the right-handed up- (left plot) or bottom-quarks (right plot) as a function of the dark matter mass, for the mass ratios , 1.1 and 2.0. These lines scale with the mass as for the spin independent dominated capture (as is the case of couplings to the bottom quark) and as for the spin dependent dominated capture (as is the case of couplings to the up-quark).
^{3}^{3}3These scalings follow from the dependence of the annihilation cross section with the mass , which implies an annihilation rate , as well as the dependences of the capture cross sections and , which translate into a capture rate and for the spin dependent and the spin independent capture, respectively. Inserting these dependences in Eq. (14) one obtains and for the spin dependent and the spin independent dominated captures. Below these lines the annihilation rate becomes very suppressed, due to the small number of particles trapped inside the Sun, and therefore observing a high energy neutrino signal from the Sun becomes very challenging, if not impossible in practice.

We also show, for comparison, the values of the coupling calculated in [33], that lead to the correct dark matter relic density. Note the existence of a lower limit on the dark matter mass from the requirement of reproducing the correct relic abundance when the scalar is highly degenerate in mass with the dark matter particle, and which is due to the impact of coannihilations in the mass degenerate regime (for details, see [33],[34]). Concretely, for couplings to up-quarks this lower limit is roughly 1 TeV and 150 GeV for and 1.1, respectively. Furthermore, for thermally produced dark matter particles, there are upper and lower limits on the dark matter mass from the requirement of equilibration between captures and annihilations in the Sun. Namely, , and for , 1.1 and 2.0, respectively. For masses not fulfilling these conditions, the number of dark matter particles captured in the Sun is suppressed compared to the equilibrium values and, accordingly, the annihilation rate. Interestingly, for dark matter particles coupling only to bottom quarks, the coupling required to thermally produce the dark matter particles is smaller than the lower limit from the equilibration condition Eq. (14) (this can be understood from the fact that ). Hence considerably higher couplings are needed in order to achieve equilibrium and therefore in this scenario the commonly used approximation does not hold but instead the annihilation rate is always smaller than this value.

The values of the coupling that lead to equilibration in the Earth are shown in Fig. 2, bottom plots. As apparent from the plots, within our framework a thermal relic can never be in equilibrium in the Earth. Therefore, the prospects for the detection of a neutrino flux from the center of the Earth are poor and hence we will restrict our discussion in the following to neutrino fluxes from dark matter annihilations in the Sun.

## 4 High energy neutrino flux from the Sun

The quarks and gauge bosons produced in the annihilations will hadronize and decay, eventually producing a neutrino flux. The dense medium inside the Sun where the annihilations take place significantly affects the energy spectrum of neutrinos, since the particles produced in the annihilations may lose a substantial fraction of their energy due to electromagnetic and strong interactions with the medium before decaying. We have simulated the energy loss of particles in the interior of the Sun following [61]. We first calculate the primary spectrum of particles produced in the annihilation with PYTHIA8.1 [62] interfaced with CalcHEP [63, 64]. We then assume that the muons and light hadrons, such as pions and kaons, are stopped in the Sun before decaying and hence produce neutrinos with energies below 1 GeV, which we neglect in our analysis (the low energy neutrino flux from dark matter annihilations can also be used to probe dark matter models, as discussed in [65],[66]). On the other hand, taus and heavy hadrons, such as charmed or beauty hadrons, decay in flight after losing a fraction of their kinetic energy. We have simulated the energy loss and subsequent decay in flight of these particles in an event-by-event basis following the chain of scatterings they undergo inside the Sun from the first scattering to the last (contrary to [67] which follows the chain of scatterings from the last to the first; the difference between both approaches turns out to be, however, negligible for our analysis). In addition to the neutrinos produced in hadronic and leptonic decays there is a component in the spectrum of hard neutrinos produced by the prompt decay of the boson from the annihilation . This process, despite having a small branching ratio, can be an important source of high energy neutrinos from the Sun and can even be the dominant contribution to the flux at the highest energies, as we will show below.

Neutrinos produced in dark matter annihilations then propagate from the solar interior to the surface, undergoing flavour oscillations and scatterings off solar matter, and then from the solar surface to a neutrino telescope at the Earth, undergoing flavour oscillations in vacuum. To calculate the fluxes at the Earth, we use WimpSim [67], a publicly available Monte Carlo code that simulates, in an event-by-event basis, the propagation of neutrinos from the production point in the Solar center to the Earth, including the effect of neutrino interactions in matter and three-flavour neutrino oscillations; for the calculation we have adopted the most recent best fit values for the neutrino oscillation parameters [68].

We show in Fig. 3 the differential flux at Earth of muon neutrinos and antineutrinos produced in the annihilation in the Sun of dark matter particles with mass which couple only to right-handed up-quarks (top panels) or to right-handed bottom-quarks (bottom panels), for the exemplary cases with mass degeneracy parameter (left panels) and (right panels). In the scenario where the dark matter particle couples to up-quarks, the contribution to the neutrino flux from annihilations into is negligible, since this annihilation channel produces mostly light hadrons that, as explained above, do not contribute to the high energy neutrino flux and besides the branching fraction for this process is helicity suppressed. The annihilation into is not helicity suppressed and gives a larger contribution to the total flux, although still suppressed since this channel also produces mostly light hadrons. More relevant is then the annihilation into , which has the largest branching fraction (even when is sizable, cf. Fig. 1) and moreover hadronizes producing a significant amount of heavy hadrons, which decay in flight contributing to the high energy neutrino flux. Finally, the annihilation channel into has, despite the fairly small branching fraction, a notable impact on the high energy neutrino flux, due to the hard neutrinos produced in the decay . In fact, for , this channel is the dominant source of neutrinos when . Therefore, this channel will be particularly important for IceCube, which is mostly sensitive to high energy neutrinos. Note also that in the limit , which is the relevant one for IceCube, the total neutrino spectrum is fairly insensitive to the degeneracy parameter, as long as . For small and moderate degeneracy parameters, the gauge bosons produced by the internal bremsstrahlung have similar energy spectra and, accordingly, also the neutrino specta. Furthermore, as discussed in [41], when the cross sections for the two-to-three processes are in fixed relations which depend only on the couplings of the gauge bosons to the final fermions. As a result, the total neutrino energy spectrum depends, for a given dark matter mass, very mildly on the degeneracy parameter.

The conclusions for the case of a dark matter particle coupling to the up-quark also apply for the coupling to the bottom-quark. In this case, the annihilation cross section into is not as suppressed as into , due to the larger bottom quark mass, although for large dark matter masses the helicity suppression is still very strong and the dominant annihilation channels are the two-to-three processes. As apparent from the plot, also for the case of a dark matter particle coupling to the b-quark, the largest contribution to the neutrino spectrum at the highest energies is the annihilation .

## 5 Limits

The muon (anti-)neutrino flux produced in the dark matter annihilations in the Sun generates a (anti-)muon signal which can be detected in a neutrino telescope. For a given neutrino spectrum we calculate the induced number of (anti-)muon events in IceCube following the approach of [69], using the effective area presented in [70]. The background consists of muons induced by the atmospheric neutrino flux and is given, together with the actual data, in [30]. For our analysis, it is necessary to choose a cut on the angle between the reconstructed muon direction and the Sun; to optimize the search we use, for each dark matter mass and mass ratio, the angle that gives the best constraint under a background only hypothesis. More details of our method of calculating the limits from the IceCube data can be found in appendix A.

The non-observation in IceCube-79 of an excess of events with respect to the expectations from the atmospheric background then allows to exclude regions of the parameter space of the model. We show in Fig. 4 as a light red band the limits on the coupling as a function of the dark matter mass, in a toy model where the dark matter particle only couples to the right-handed up-quarks for three exemplary choices of the degeneracy parameter, , 1.1 and 2.0. The width of the band brackets the sensitivity of the upper limit to systematic uncertainties in the determination of the capture rate, which were estimated to be () for spin dependent (spin independent) capture [70], as well as to uncertainties in the determination of the astrophysical parameter , the nuclear parameters , and the second moments of the quark PDFs. The values of the corresponding uncertainties are chosen as in [33]. We also show for comparison the limits derived in [33] on the coupling from the non-observation of an excess in the cosmic antiproton-to-proton fraction measured by PAMELA (grey band), and which in this model stem mostly from the two-to-three process , or from the non-observation of a dark matter signal in the XENON100 experiment (dark blue band). Again, the width of the band brackets the sensitivity of the upper limit to the astrophysical and nuclear uncertainties. As apparent from the plot, IceCube-79 gives limits which are competitive, and for small degeneracies better, than the PAMELA antiproton limits. On the other hand, the IceCube-79 limits are worse than the XENON100 limits, due to the fact that this toy model has both spin dependent and spin independent interactions with matter, the latter being strongly constrained by direct search experiments.

We also show in the plot the region of the parameter space where, for this toy model, captures and annihilations are in equilibrium in the Sun and therefore the annihilation rate reaches its maximum value. The shaded blue region then corresponds to model parameters where the neutrino signal from the Sun is very suppressed due to the inefficient dark matter capture and therefore observing an annihilation signal from the Sun becomes very challenging. As can be seen from the plot, IceCube-79, and especially XENON100, already probe a fairly large region of the parameter space where a signal of annihilations from the Sun could be found, especially for scenarios with large mass degeneracy. The XENON1T experiment will close in on the parameter space even further, concretely, the factor of increase in sensitivity projected in the mass range of interest for our analysis [71], will translate into an improvement in the limit on the coupling by a factor of . Note also that XENON100 rules out the values of the couplings necessary to achieve equilibriation inside the Earth, as can be checked from comparing Fig. 2, bottom plot, with Fig. 4. Therefore, the observation of a high energy neutrino signal from the center of the Earth is very unlikely in this scenario.

The limits of the coupling can be straightforwardly translated into limits of the spin dependent cross section, as commonly presented in searches for neutrinos from dark matter annihilations in the Sun. The result is shown in Fig. 5 for our exemplary cases , 1.1 and 2.0. We show the limit on the spin dependent cross section that follows from assuming annihilations in the Sun only into as well as the limit that results from considering also the two-to-three annihilations. As apparent from the plot, including the three body annihilations in the Sun dramatically improves the limits on this scenario. We also show for completeness the experimental limits on the spin dependent cross section of this model from the COUPP [72] and the XENON100 [35] experiments. The lines always correspond, where applicable, to the choice of astrophysical and nuclear parameters giving the most conservative constraints. It should be stressed that the XENON100 limit in this plot refers to the limit on the spin dependent cross section from the non-observation of events in the XENON100 experiment in this specific model, which produces spin dependent and spin independent dark matter-nucleon interactions. The limits then differs from those reported in [73, 74], where it was implicitly assumed a model with just spin dependent interactions.

The limits for the scenario where the dark matter particle couples just to the right-handed bottom-quark are shown in Fig. 6. As explained at the end of Section 2, we restrict our analysis to , therefore we only show in this case the limits for and , for which the previous condition is automatically fulfilled provided . Due to the suppressed spin dependent scattering rate, and as is apparent from the plots, for couplings to bottom-quarks the IceCube limits are rather weak, weaker than the antiproton limits and sensibly weaker than the XENON100 limits. It also follows from the plot that the XENON100 limits rule out already a large part of the parameter space where equilibration in the Sun is possible; the XENON1T experiment will improve the limit in the coupling by a factor , if no signal is found.

In deriving these limits we have not made any assumption on how the dark matter particle was produced. It is then interesting to analyze the limits under the well motivated assumption that the dark matter particle was produced thermally in the early Universe. This fixes one of the three parameters of the toy model, , or . The value of the coupling constant required to produce the correct dark matter abundance is shown in Figs. 4,6 as a dashed line and lies well below the IceCube-79 limits, although close to, and for some choices of parameters above, the XENON100 limit.

We further investigate in Fig. 7 the complementarity of the various search strategies to probe this scenario, under the assumption that the observed dark matter abundance was thermally produced. We have calculated, following [34] (see also [75, 76, 77]), the excluded regions of the parameter space, spanned by the dark matter mass and the mass splitting between and , from the XENON100 results (light green) and from the LHC searches of a colored scalar by the CMS collaboration employing the -analysis based on 11.7 at 8 TeV center of mass energy [37] (light red). The present limits for the model with couplings to just right-handed up- (bottom-) quarks are shown in the upper (lower) left plots, together with the regions of the parameter space where it is impossible to generate thermally the observed relic abundance of dark matter particles (light gray) and the region where the perturbative calculation is not valid (dark gray). We also show in the plot the lines of constant (shown in black) and the line that marks the boundary of the parameter space where equilibration in the Sun is not attained, (shown in red). Note that indirect search experiments with antiprotons and with high energy neutrinos are not yet sensitive to probe thermally produced dark matter particles, cf. Figs.4,6, hence these experiments do not rule out any point in this parameter space. For couplings to up-quarks, XENON100 is the experiment that gives the strongest constraints on the parameter space at low mass splittings, and in particular rules out the possibility of observing a signal from the Sun for when ; for larger mass splittings, it is the CMS experiment the one giving the strongest constraints. Future experiments, notably the XENON1T experiment, will cover a much larger region of the parameter space where a dark matter signal from the Sun could be observed, as apparent from Fig. 7, upper right plot.

For couplings to bottom-quarks, limits from present experiments are very weak, as can be seen in Fig. 7, lower left plot. In particular, the region of the parameter space where equilibration in the Sun takes place, around and , is barely constrained by the XENON100 experiment. Interestingly, the upcoming XENON1T experiment will be able to cover the whole region where equilibration is predicted to take place in the Sun, cf. Fig. 7, lower right plot. The non-observation of a dark matter signal at XENON1T will then have important implications for the feasibility of observing a high energy neutrino signal from the Sun in this concrete scenario.

## 6 Conclusions

The observation of an excess of high energy neutrinos in the direction of the Sun would be a clear signature of the annihilation of dark matter particles captured in the solar interior. In this paper we have investigated the prospects to observe this signature in a minimal scenario where the dark matter is a Majorana fermion that couples to a right-handed quark, concretely an up-quark or a bottom-quark, and a colored scalar via a Yukawa coupling. This scenario then allows to study the capture and the annihilation in a single Particle Physics framework. We have found that, for small and moderate mass splitting between the dark matter and the colored scalar, the most important channel in the calculation of the number of captured dark matter particles in the Sun is the two-to-three annihilation into . We have determined, for a given dark matter and colored scalar mass, the minimal value of the Yukawa coupling necessary to attain equilibrium between capture and annihilation. Upon comparing with the value of the coupling necessary to thermally produce the observed dark matter abundance, we conclude that for a thermal relic equilibration is never reached in the interior of Sun if the dark matter particle only couples to bottom-quarks. Equilibration in the interior of the Earth is not attained neither for couplings to up-quarks nor to bottom-quarks. Therefore, the observation of a high-energy neutrino flux from the Earth is, although not precluded in theory, very challenging in practice due to the very suppressed annihilation rate.

We have also calculated the neutrino flux produced in the dark matter annihilations including the two-to-three annihilations into a quark-antiquark pair and a gauge boson. We have found that the neutrino flux at the Earth is, close to the kinematical endpoint, dominated by the annihilations into , despite the rather small branching fraction in this channel. We have then derived an upper limit on the Yukawa coupling from the non-observation at IceCube of an excess of high-energy neutrino events in the direction of the Sun. For couplings to up-quarks, and small mass degeneracies between the dark matter particle and the colored scalar, the IceCube limits are stronger than the limits on the same model parameters from the PAMELA measurements of the cosmic antiproton-to-proton fraction. On the other hand, for couplings to bottom-quarks the IceCube limits are weaker, due to the inefficient capture of dark matter particles in the Sun in this scenario. Both for couplings to up- and to bottom-quarks, the direct detection limits on the coupling from XENON100 are significantly stronger than the IceCube limits and already exclude regions of the parameter space where equilibration between capture and annihilation is reached inside the Sun.

We have investigated in detail the limits on the parameter space under the assumption that the observed dark matter abundance was thermally generated, calculating the excluded regions from present searches of new physics at LHC and from XENON100. For couplings to bottom-quarks, the present limits are rather weak, while for couplings to up-quarks, they already exclude some regions of the parameter space. Future direct detection experiments, such as LUX and notably XENON1T will continue closing in on the the possibility of observing a high energy neutrino flux from the Sun. In particular, in the scenario where the dark matter couples just to right-handed bottom-quarks, XENON1T will cover the whole region of the parameter space where equilibration is reached in the solar interior. Therefore, the non-observation of a dark matter signal at XENON1T will then make very unlikely the possibility of observing a high energy neutrino signal from the Sun in this scenario.

## Note Added

After the completion of this work, the LUX experiment released new limits on the spin independent dark matter-nucleon scattering cross section [78], which are approximately a factor of three better than the XENON100 limits in the range of masses of interest for our work. In view of the new results, our limits on the coupling constant are improved by approximately 30%.

## Acknowledgements

We are grateful to Sergio Palomares-Ruiz, Miguel Pato and Stefan Vogl for useful discussions. This work was partially supported by the DFG cluster of excellence “Origin and Structure of the Universe,” the Universität Bayern e.V., the TUM Graduate School and the Studienstiftung des Deutschen Volkes.

## Appendix A Calculation of limits from the IceCube data

In this appendix, we shortly review our method of converting a (anti-)neutrino flux from dark matter annihilations in the Sun into a number of (anti-)muon events in IceCube as well as the statistical method of extracting limits from the data, closely following [69].

The number of expected muon signal events, , for a given (anti-)neutrino flux is obtained by [69]

(16) |

Here, is the live-time of the 79-string analysis of IceCube given in [70] and is the effective area for the detection of neutrinos/antineutrinos in IceCube, encapsulating all information about the efficiency of the IceCube detector as well as the scattering cross sections of neutrinos and the muon range in rock and ice. We use the total effective area for the IceCube 79-string configuration given in [70] and we convert it into separate effective areas for neutrinos and antineutrinos using the method described in [79]. Besides, is the angular loss factor, defined as the probability that the angle of the reconstructed muon track, originated from an incoming neutrino with energy , with respect to the direction of the Sun is smaller than a certain cut angle , namely

(17) |

This assumes a 2D Gaussian point spread function of the muon track directions (see [69] for more details); we extract the corresponding median angle between the neutrino and the muon track from [70].

We then use the standard method for calculating upper limits on , based on a hypothesis test using the likelihood ratio

(18) |

For a definition of the likelihood functions and the corresponding -values of the hypothesis test we refer to [69]. The number of observed events and the number of off-source measured background events as a function of the angle between the muon direction and the Sun are given in [30].

For our analysis we choose the cut angle such that the limit on the annihilation cross section is optimized. To this end, we use the number of neutrino events observed by IceCube in directions away from the Sun, which are presumably of atmospheric origin, and we make the plausible assumption that the background of atmospheric neutrinos in the direction of the Sun is identical to the background in these directions. We then calculate, for given and , the limit on the annihilation cross section for eight different cut angles between and (corresponding to the the bins in used in [30]) from the requirement that the dark matter signal would not produce an excess over the atmospheric neutrino background with a significance larger than 95% C.L. In this step of the calculation, we use the off-source measured background events instead of the real data (i.e. the data in the direction of the sun); hence our procedure to determine the optimal cut angle is not biased by the actual measurements in the direction of the Sun, which might receive an exotic component. Lastly, we select the cut angle that gives the best limit under this “background-only” hypothesis, and use it to calculate the limit on the annihilation cross section from the actual number of neutrino events measured by the IceCube collaboration in the direction of the Sun.

The optimal cut angle as a function of the dark matter mass is shown in Fig. 8 for for couplings to up- and to bottom-quarks. The differences in the spectra for couplings to up- or to bottom-quarks shown in Fig. 3 get translated into different choices for the optimal cut angle in Fig. 8, due to the decrease of the average scattering angle between a neutrino and the muon produced in a charged-current interaction with the energy. Therefore, in a model with a neutrino spectrum peaked at high-energies, as in our toy model with coupling to up-quarks, small cut angles are favored, as in this case most of the produced muons have directions within a small angle around the Sun; choosing a larger cut angle would increase the number of expected background events and hence lead to worse constraints. Conversely, in a model where the neutrino spectrum has an important soft component, as in our toy model with coupling to bottom-quarks, more signal events are found at larger angles and hence larger cut angles are favored, otherwise too many signal events would be cut, leading again to a worse constraint.

As a final comment, we note that the published 79-string data of IceCube [30] actually consists of three data sets, named “Winter high”, “Winter low” and “Summer low”. These are three non-overlapping data sets, each with its own distribution of , and as well as its own effective area. The three event selections thereby correspond to different sets of experimental cuts, each optimized for a different signal model; in particular the “Summer low” data set consists of downward going muons which were detected in DeepCore, using the surrounding IceCube strings as an active veto against atmospheric muons. While it would be in principle possible to use all three data sets in a combined likelihood for calculating the constraints, we choose the simpler (and conservative) strategy of calculating for each point in the parameter space three independent constraints using the three event selections and then choose the most constraining one. In our model, this turns out to be almost always the “Winter high” data set, except for dark matter masses below 200 GeV, where the “Winter low” or “Summer low” event selection can be relevant.

## References

- [1] G. Bertone, D. Hooper, and J. Silk, Particle dark matter: Evidence, candidates and constraints, Phys.Rept. 405 (2005) 279–390, [hep-ph/0404175].
- [2] L. Bergstrom, Nonbaryonic dark matter: Observational evidence and detection methods, Rept.Prog.Phys. 63 (2000) 793, [hep-ph/0002126].
- [3] T. Bringmann and C. Weniger, Gamma Ray Signals from Dark Matter: Concepts, Status and Prospects, Phys.Dark Univ. 1 (2012) 194–217, [arXiv:1208.5481].
- [4] M. Srednicki, S. Theisen, and J. Silk, Cosmic Quarkonium: A Probe of Dark Matter, Phys.Rev.Lett. 56 (1986) 263.
- [5] S. Rudaz, Cosmic production of quarkonium?, Phys.Rev.Lett. 56 (1986) 2128.
- [6] L. Bergstrom and H. Snellman, Observable monochromatic photons from cosmic photino annihilation, Phys.Rev. D37 (1988) 3737–3741.
- [7] L. Bergstrom, Radiative processes in dark matter photino annihilation, Phys.Lett. B225 (1989) 372.
- [8] R. Flores, K. A. Olive, and S. Rudaz, Radiative processes in LSP annihilation, Phys.Lett. B232 (1989) 377–382.
- [9] T. Bringmann, L. Bergstrom, and J. Edsjo, New Gamma-Ray Contributions to Supersymmetric Dark Matter Annihilation, JHEP 0801 (2008) 049, [arXiv:0710.3169].
- [10] A. Ibarra, S. Lopez Gehler, and M. Pato, Dark matter constraints from box-shaped gamma-ray features, JCAP 1207 (2012) 043, [arXiv:1205.0007].
- [11] Fermi-LAT Collaboration, M. Ackermann et. al., Search for Gamma-ray Spectral Lines with the Fermi Large Area Telescope and Dark Matter Implications, arXiv:1305.5597.
- [12] H.E.S.S. Collaboration, A. Abramowski et. al., Search for photon line-like signatures from Dark Matter annihilations with H.E.S.S, Phys.Rev.Lett. 110 (2013) 041301, [arXiv:1301.1173].
- [13] T. Bringmann, X. Huang, A. Ibarra, S. Vogl, and C. Weniger, Fermi LAT Search for Internal Bremsstrahlung Signatures from Dark Matter Annihilation, arXiv:1203.1312.
- [14] A. Ibarra, H. M. Lee, S. Lopez Gehler, W.-I. Park, and M. Pato, Gamma-ray boxes from axion-mediated dark matter, JCAP 1305 (2013) 016, [arXiv:1303.6632].
- [15] LAT Collaboration, M. Ackermann et. al., Fermi LAT Search for Dark Matter in Gamma-ray Lines and the Inclusive Photon Spectrum, Phys.Rev. D86 (2012) 022002, [arXiv:1205.2739].
- [16] H.E.S.S.Collaboration, A. Abramowski et. al., Search for a Dark Matter annihilation signal from the Galactic Center halo with H.E.S.S, Phys.Rev.Lett. 106 (2011) 161301, [arXiv:1103.3266].
- [17] M. Ackermann, M. Ajello, A. Allafort, L. Baldini, J. Ballet, et. al., Constraints on Dark Matter Annihilation in Clusters of Galaxies with the Fermi Large Area Telescope, JCAP 1005 (2010) 025, [arXiv:1002.2239].
- [18] HESS Collaboration, A. Abramowski et. al., Search for Dark Matter Annihilation Signals from the Fornax Galaxy Cluster with H.E.S.S, Astrophys.J. 750 (2012) 123, [arXiv:1202.5494].
- [19] Fermi-LAT Collaboration, M. Ackermann et. al., Dark Matter Constraints from Observations of 25 Milky Way Satellite Galaxies with the Fermi Large Area Telescope, arXiv:1310.0828.
- [20] A. Geringer-Sameth and S. M. Koushiappas, Exclusion of canonical WIMPs by the joint analysis of Milky Way dwarfs with Fermi, Phys.Rev.Lett. 107 (2011) 241303, [arXiv:1108.2914].
- [21] MAGIC Collaboration, J. Aleksic et. al., Searches for Dark Matter annihilation signatures in the Segue 1 satellite galaxy with the MAGIC-I telescope, JCAP 1106 (2011) 035, [arXiv:1103.0477].
- [22] HESS Collaboration, A. Abramowski et. al., H.E.S.S. constraints on Dark Matter annihilations towards the Sculptor and Carina Dwarf Galaxies, Astropart.Phys. 34 (2011) 608–616, [arXiv:1012.5602].
- [23] L. Bergstrom, T. Bringmann, I. Cholis, D. Hooper, and C. Weniger, New limits on dark matter annihilation from AMS cosmic ray positron data, Phys.Rev.Lett. (2013) [arXiv:1306.3983].
- [24] A. Ibarra, A. S. Lamperstorfer, and J. Silk, Dark matter annihilations and decays after the AMS-02 positron measurements, arXiv:1309.2570.
- [25] IceCube collaboration, R. Abbasi et. al., Search for Neutrinos from Annihilating Dark Matter in the Direction of the Galactic Center with the 40-String IceCube Neutrino Observatory, arXiv:1210.3557.
- [26] IceCube Collaboration, M. Aartsen et. al., The IceCube Neutrino Observatory Part IV: Searches for Dark Matter and Exotic Particles, arXiv:1309.7007.
- [27] W. H. Press and D. N. Spergel, Capture by the sun of a galactic population of weakly interacting massive particles, Astrophys.J. 296 (1985) 679–684.
- [28] J. Silk, K. A. Olive, and M. Srednicki, The Photino, the Sun and High-Energy Neutrinos, Phys.Rev.Lett. 55 (1985) 257–259.
- [29] A. Gould, Resonant Enhancements in WIMP Capture by the Earth, Astrophys.J. 321 (1987) 571.
- [30] IceCube collaboration, M. Aartsen et. al., Search for dark matter annihilations in the Sun with the 79-string IceCube detector, Phys.Rev.Lett. 110 (2013) 131302, [arXiv:1212.4097].
- [31] N. F. Bell, A. J. Brennan, and T. D. Jacques, Neutrino signals from electroweak bremsstrahlung in solar WIMP annihilation, JCAP 1210 (2012) 045, [arXiv:1206.2977].
- [32] K. Fukushima, Y. Gao, J. Kumar, and D. Marfatia, Bremsstrahlung signatures of dark matter annihilation in the Sun, Phys.Rev. D86 (2012) 076014, [arXiv:1208.1010].
- [33] M. Garny, A. Ibarra, M. Pato, and S. Vogl, Closing in on mass-degenerate dark matter scenarios with antiprotons and direct detection, JCAP 1211 (2012) 017, [arXiv:1207.1431].
- [34] M. Garny, A. Ibarra, M. Pato, and S. Vogl, Internal bremsstrahlung signatures in light of direct dark matter searches, arXiv:1306.6342.
- [35] XENON100 Collaboration, E. Aprile et. al., Dark Matter Results from 225 Live Days of XENON100 Data, Phys.Rev.Lett. 109 (2012) 181301, [arXiv:1207.5988].
- [36] PAMELA Collaboration, O. Adriani et. al., PAMELA results on the cosmic-ray antiproton flux from 60 MeV to 180 GeV in kinetic energy, Phys.Rev.Lett. 105 (2010) 121101, [arXiv:1007.0821].
- [37] CMS Collaboration, S. Chatrchyan et. al., Search for supersymmetry in hadronic final states with missing transverse energy using the variables AlphaT and b-quark multiplicity in pp collisions at 8 TeV, Eur.Phys.J. C73 (2013) 2568, [arXiv:1303.2985].
- [38] F. Luo and T. You, Enhancement of Majorana Dark Matter Annihilation Through Higgs Bremsstrahlung, arXiv:1310.5129.
- [39] M. Drees, G. Jungman, M. Kamionkowski, and M. M. Nojiri, Neutralino annihilation into gluons, Phys.Rev. D49 (1994) 636–647, [hep-ph/9306325].
- [40] V. Barger, W.-Y. Keung, H. E. Logan, and G. Shaughnessy, Neutralino annihilation to q anti-q g, Phys.Rev. D74 (2006) 075005, [hep-ph/0608215].
- [41] M. Garny, A. Ibarra, and S. Vogl, Dark matter annihilations into two light fermions and one gauge boson: General analysis and antiproton constraints, JCAP 1204 (2012) 033, [arXiv:1112.5155].
- [42] M. Garny, A. Ibarra, and S. Vogl, Antiproton constraints on dark matter annihilations from internal electroweak bremsstrahlung, JCAP 1107 (2011) 028, [arXiv:1105.5367].
- [43] P. Ciafaloni, M. Cirelli, D. Comelli, A. De Simone, A. Riotto, et. al., On the Importance of Electroweak Corrections for Majorana Dark Matter Indirect Detection, JCAP 1106 (2011) 018, [arXiv:1104.2996].
- [44] P. Ciafaloni, M. Cirelli, D. Comelli, A. De Simone, A. Riotto, et. al., Initial State Radiation in Majorana Dark Matter Annihilations, JCAP 1110 (2011) 034, [arXiv:1107.4453].
- [45] P. Ciafaloni, D. Comelli, A. De Simone, A. Riotto, and A. Urbano, Electroweak Bremsstrahlung for Wino-Like Dark Matter Annihilations, JCAP 1206 (2012) 016, [arXiv:1202.0692].
- [46] N. F. Bell, J. B. Dent, A. J. Galea, T. D. Jacques, L. M. Krauss, et. al., W/Z Bremsstrahlung as the Dominant Annihilation Channel for Dark Matter, Revisited, Phys.Lett. B706 (2011) 6–12, [arXiv:1104.3823].
- [47] N. F. Bell, J. B. Dent, T. D. Jacques, and T. J. Weiler, Dark Matter Annihilation Signatures from Electroweak Bremsstrahlung, Phys.Rev. D84 (2011) 103517, [arXiv:1101.3357].
- [48] M. Kachelriess, P. Serpico, and M. A. Solberg, On the role of electroweak bremsstrahlung for indirect dark matter signatures, Phys.Rev. D80 (2009) 123533, [arXiv:0911.0001].
- [49] P. Ciafaloni, D. Comelli, A. Riotto, F. Sala, A. Strumia, et. al., Weak Corrections are Relevant for Dark Matter Indirect Detection, JCAP 1103 (2011) 019, [arXiv:1009.0224].
- [50] T. Toma, Internal Bremsstrahlung Signature of Real Scalar Dark Matter and Consistency with Thermal Relic Density, Phys.Rev.Lett. 111 (2013) 091301, [arXiv:1307.6181].
- [51] F. Giacchino, L. Lopez-Honorez, and M. H. G. Tytgat, Scalar Dark Matter Models with Significant Internal Bremsstrahlung, arXiv:1307.6480.
- [52] C. Garcia-Cely and A. Ibarra, Novel Gamma-ray Spectral Features in the Inert Doublet Model, JCAP 1309 (2013) 025, [arXiv:1306.4681].
- [53] J. Hisano, K. Ishiwata, and N. Nagata, Gluon contribution to the dark matter direct detection, Phys.Rev. D82 (2010) 115007, [arXiv:1007.2601].
- [54] M. Drees and M. Nojiri, Neutralino - nucleon scattering revisited, Phys.Rev. D48 (1993) 3483–3501, [hep-ph/9307208].
- [55] J. R. Ellis, K. A. Olive, and C. Savage, Hadronic Uncertainties in the Elastic Scattering of Supersymmetric Dark Matter, Phys.Rev. D77 (2008) 065026, [arXiv:0801.3656].
- [56] G. Jungman, M. Kamionkowski, and K. Griest, Supersymmetric dark matter, Phys.Rept. 267 (1996) 195–373, [hep-ph/9506380].
- [57] K. Griest and D. Seckel, Cosmic Asymmetry, Neutrinos and the Sun, Nucl.Phys. B283 (1987) 681.
- [58] G. Busoni, A. De Simone, and W.-C. Huang, On the Minimum Dark Matter Mass Testable by Neutrinos from the Sun, JCAP 1307 (2013) 010, [arXiv:1305.1817].
- [59] R. Catena and P. Ullio, A novel determination of the local dark matter density, JCAP 1008 (2010) 004, [arXiv:0907.0018].
- [60] M. Pato, L. Baudis, G. Bertone, R. Ruiz de Austri, L. E. Strigari, et. al., Complementarity of Dark Matter Direct Detection Targets, Phys.Rev. D83 (2011) 083505, [arXiv:1012.3458].
- [61] S. Ritz and D. Seckel, Detailed Neutrino Spectra From Cold Dark Matter Annihilations in the Sun, Nucl.Phys. B304 (1988) 877.
- [62] T. Sjostrand, S. Mrenna, and P. Z. Skands, A Brief Introduction to PYTHIA 8.1, Comput.Phys.Commun. 178 (2008) 852–867, [arXiv:0710.3820].
- [63] A. Pukhov, E. Boos, M. Dubinin, V. Edneral, V. Ilyin, et. al., CompHEP: A Package for evaluation of Feynman diagrams and integration over multiparticle phase space, hep-ph/9908288.
- [64] A. Pukhov, CalcHEP 2.3: MSSM, structure functions, event generation, batchs, and generation of matrix elements for other packages, hep-ph/0412191.
- [65] C. Rott, J. Siegal-Gaskins, and J. F. Beacom, New Sensitivity to Solar WIMP Annihilation using Low-Energy Neutrinos, Phys.Rev. D88 (2013) 055005, [arXiv:1208.0827].
- [66] N. Bernal, J. Martin-Albo, and S. Palomares-Ruiz, A novel way of constraining WIMPs annihilations in the Sun: MeV neutrinos, JCAP 1308 (2013) 011, [arXiv:1208.0834].
- [67] M. Blennow, J. Edsjo, and T. Ohlsson, Neutrinos from WIMP annihilations using a full three-flavor Monte Carlo, JCAP 0801 (2008) 021, [arXiv:0709.3898].
- [68] Particle Data Group, J. Beringer et. al., Review of Particle Physics (RPP), Phys.Rev. D86 (2012) 010001.
- [69] IceCube Collaboration, P. Scott et. al., Use of event-level neutrino telescope data in global fits for theories of new physics, JCAP 1211 (2012) 057, [arXiv:1207.0810].
- [70] M. Danninger, PhD thesis, Stockholm University.
- [71] XENON1T collaboration, E. Aprile, The XENON1T Dark Matter Search Experiment, arXiv:1206.6288.
- [72] COUPP Collaboration, E. Behnke et. al., First Dark Matter Search Results from a 4-kg CFI Bubble Chamber Operated in a Deep Underground Site, Phys.Rev. D86 (2012) 052001, [arXiv:1204.3094].
- [73] M. Garny, A. Ibarra, M. Pato, and S. Vogl, On the spin-dependent sensitivity of XENON100, Phys.Rev. D87 (2013), no. 5 056002, [arXiv:1211.4573].
- [74] XENON100 Collaboration, E. Aprile et. al., Limits on spin-dependent WIMP-nucleon cross sections from 225 live days of XENON100 data, Phys.Rev.Lett. 111 (2013) 021301, [arXiv:1301.6620].
- [75] H. An, L.-T. Wang, and H. Zhang, Dark matter with -channel mediator: a simple step beyond contact interaction, arXiv:1308.0592.
- [76] Y. Bai and J. Berger, Fermion Portal Dark Matter, arXiv:1308.0612.
- [77] A. DiFranzo, K. I. Nagao, A. Rajaraman, and T. M. P. Tait, Simplified Models for Dark Matter Interacting with Quarks, arXiv:1308.2679.
- [78] LUX Collaboration, D. Akerib et. al., First results from the LUX dark matter experiment at the Sanford Underground Research Facility, arXiv:1310.8214.
- [79] C. A. Arguelles and J. Kopp, Sterile neutrinos and indirect dark matter searches in IceCube, JCAP 1207 (2012) 016, [arXiv:1202.3431].