Search for Dark Matter Annihilation in the Galactic Center with IceCube-79

The Milky Way is expected to be embedded in a halo of dark matter particles, with the highest density in the central region, and decreasing density with the halo-centric radius. Dark matter might be indirectly detectable at Earth through a flux of stable particles generated in dark matter annihilations and peaked in the direction of the Galactic Center. We present a search for an excess flux of muon (anti-) neutrinos from dark matter annihilation in the Galactic Center using the cubic-kilometer-sized IceCube neutrino detector at the South Pole. There, the Galactic Center is always seen above the horizon. Thus, new and dedicated veto techniques against atmospheric muons are required to make the southern hemisphere accessible for IceCube. We used 319.7 live-days of data from IceCube operating in its 79-string configuration during 2010 and 2011. No neutrino excess was found and the final result is compatible with the background. We present upper limits on the self-annihilation cross-section, $\left<\sigma_\mathrm{A} v\right>$, for WIMP masses ranging from 30 GeV up to 10 TeV, assuming cuspy (NFW) and flat-cored (Burkert) dark matter halo profiles, reaching down to $\simeq 4 \cdot 10^{-24}$ cm$^3$ s$^{-1}$, and $\simeq 2.6 \cdot 10^{-23}$ cm$^3$ s$^{-1}$ for the $\nu\overline{\nu}$ channel, respectively.


Introduction
The first clear evidence for the existence of an invisible mass component in the universe was Zwicky's observation of the dynamics of the Coma galaxy cluster [1].Subsequently, a broad range of cosmological and astrophysical observations supported the existence of this dark matter (DM) on various scales, from galaxy cluster scales down to galactic scales.Measurements of galactic velocity profiles hint at invisible mass distributed beyond the visible disks [2].Galaxy cluster dynamics exhibit a similar behavior [3].
Further evidence for the existence of dark matter can be found in galaxy cluster mergers like the Bullet Cluster [4,5].Following a collision, the interstellar and intergalactic gas components as seen in X-ray observations are spatially separated from the reconstructed mass distribution.Such a separation strongly disfavors theories of modified gravity.
According to the current understanding of the formation and evolution of large-scale structures, cold (nonrelativistic), or warm dark matter is preferred over hot (relativistic) dark matter.Otherwise, the formation of the observed large-structures on time scales of the order of the age of the universe would not have been possible [6][7][8].
Though the nature of dark matter is largely unknown, some of its properties may be deduced from the above-mentioned observations.Analyses of temperature fluctuations in the Cosmic Microwave Background (CMB) by the Planck collaboration [8] yield the current best estimate for the total content of DM in the universe: Ω CDM h 2 = 0.1199 ± 0.0027, with the cold DM density parameter Ω CDM , and h = 0.673 ± 0.012 being the Hubble parameter divided by 100 km/s Mpc.
Besides inference from gravitational interaction, particle DM may also be detected indirectly.A weakly interacting massive particle (WIMP) at roughly GeVscale masses is a favorable class of DM; it naturally provides the right order of magnitude for the thermal relic abundance of DM in the early universe [9].Examples of WIMPs are neutralinos in supersymmetric extensions of the Standard Model [10], or the lightest stable excitations in Kaluza-Klein models [11].
If DM decays, or (self-)annihilates, a flux of stable final-state messenger particles, e.g.charged leptons, photons, and neutrinos, may be detected at Earth, making DM experimentally accessible by indirect searches (e.g.[12][13][14]).The neutrino is an attractive messenger particle because it propagates without absorption, and neutrino vacuum oscillations do not alter the energy and direction information.Further, no fore-or background from astrophysical objects has been identified yet.Regions of increased DM density, like massive celestial objects, dwarf galaxies, galactic halos, and the Galactic Center, provide targets to search for an increased flux of neutrinos.Due to its proximity, the Galactic Center is expected to yield the highest flux of annihilation products [15].While most of these sources would appear as (nearly) point-like sources in the sky, the Galactic Center is an extended source, and a signal from the Milky Way halo would lead to a large-scale anisotropy in neutrino arrival directions [16,17].With its 4π acceptance, the IceCube neutrino detector [18], is well-suited for DM searches from all of the abovementioned sources.
In this paper we present the results from a search for a neutrino signal from DM self-annihilation in the Galactic Center, targeting DM masses ranging from 30 GeV to 10 TeV.Due to the wide range of event topologies associated with neutrinos from this energy range, two event selections are motivated and presented.One event selection focuses on the low-mass region ranging from 30 GeV to 100 GeV, accessible through the low-energy in-fill array DeepCore (DC) [19], with the surrounding parts of IceCube used as veto.The sec-ond event selection focuses on the mass range 100 GeV-1 TeV, but extends up to 10 TeV.For this selection a larger part of the IceCube detector is defined as fiducial volume.Throughout this paper we denote the low-mass event selection as LE and the high-mass selection as HE.
2 Dark Matter Halos DM halos are considered to be gravitationally self-bound overdensities of DM particles, formed through hierarchical merging of proto-halos from primordial density fluctuations [20].There is a tension between halo profile fits to DM overdensities in N-body simulations, and fits to observational data (the cusp-core problem) [21].N-body simulations seem to favor cuspy halos, while observations of e.g.dwarf spheroidal galaxies imply a rather flat central core region.However, the inner part of the halo profile may depend on the host halo mass [22].Simple models of DM halos describe the density by a smooth spherically-symmetric function of the halocentric radius, with the maximal density at the center, and a decreasing density with increasing radius.One parametrization of such density profiles is given by (modified from [23]) with the shape parameters α, β, γ, the scale radius r s , and the mass density normalization ρ 0 , which is usually determined from the assumed local DM density, ρ local , in our solar system.We introduced the parameter δ to allow for a central core if set to 1, while δ = 0 describes a cuspy halo profile.The parametrization of Eq. ( 1) is a combination of power laws, where e.g. the power-law index γ describe the inner slope, while α and β describe the outer slope.The halo-centric distance of this transition region depends on the scale radius r s .
In view of the unresolved cusp-core problem, we present results for two halo density profiles.The widely used Navarro-Frenk-White (NFW) profile represents cuspy halos [24], and is chosen for comparability among different experimental results.The Burkert profile is chosen as representative of flat-cored profiles [25].Based on observation, the latter profile is currently favored for the Milky Way [26].Table 1 shows the parameter values for the two models used in this work.The flux of final state particles from annihilating dark matter depends on the DM mass density squared, integrated along the line-of-sight (l.o.s.) through the DM halo, and is given by the J a -factor.Following e.g.[13,15], the J a -factor is Here, the density profile along the l.o.s. is parametrized for a given angle between the l.o.s. and the direction of the center of the galaxy, Ψ .The parameters are the radius of the solar circle, R SC ≈ 8.5 kpc, and the maximal distance from the observer along the l.o.s., l max .The latter is with the assumed radius of the Milky Way, R MW ≈ 50 kpc.Typically, radii larger than the scale radius do not contribute significantly to the total value of J a .
with σ A v being the thermally averaged product of self-annihilation cross-section, σ A and WIMP velocity v. Further, m χ is the WIMP mass, dN ν /dE is the neutrino energy spectrum per annihilating WIMP pair, and J a is the DM abundance along the l.o.s..We consider several benchmark annihilation channels with 100% branching ratios (χχ → b b, W + W − , µ + µ − , τ + τ − ,ν ν) for the calculation of the neutrino energy spectrum dN ν /dE.The resulting spectra bracket the realistic annihilation neutrino energy spectra with a mixture of different annihilation branching ratios.We generated a neutrino energy spectrum from annihilating DM for a particular WIMP mass and annihilation channel using the PYTHIA8 (version 8.175) software package [27].Our PYTHIA8 simulation was set up to simulate a generic resonance with an energy of twice the WIMP mass forming only the final state particle pair (i.e.bb, W + W − , etc.) in question.Subsequent processes like hadronization and decay were simulated using the default PYTHIA8 implementations.The generic resonance ensures an isotropic decay of weak bosons, e.g. in the W + W − channel.Thus, the spin of the annihilating WIMPs is not considered and we don't assume a specific WIMP model like the lightest neutralino described by supersymmetric models [9].If the WIMP is indeed the lightest neutralino, the spin (1/2) of the WIMP would affect the generation of the neutrino energy spectra.The spin of such a WIMP would lead to fully transversely polarized W -bosons in the final state of the annihilation process, thus altering the neutrino energy spectrum [28].The differences in the differential neutrino yield compared to the isotropically decaying W bosons is about ±40%.Examples of the neutrino energy spectra used here are shown in the bottom panel of Fig. 1 for the bb, W + W − , and µ + µ − annihilation channels.
In general, neutrinos are subject to neutrino oscillations on the way from the Galactic Center to the Earth.Due to the very long baseline, we assume a relative neutrino flavor ratio of 1:1:1 at Earth.These simulations have the detector response folded in and are generated using the ANIS event generator [29] modelling the neutrino-nucleon charged and neutral current interactions via CTEQ5 [30] parton distributions for neutrino and anti-neutrino interactions.
Finally, the neutrino energy distributions are used to weight generic simulated neutrino data to DM annihilation signal.

The IceCube Neutrino Observatory
The IceCube Neutrino Observatory, situated at the geographic South Pole, consists of an in-ice detector array, IceCube, and a surface air shower detector array, Ice-Top [31], dedicated to neutrino and cosmic ray research, respectively.IceCube [32] is installed in the glacial ice at depths between 1450 m and 2450 m below the surface, instrumenting a total volume of one cubic kilometer.IceCube detects neutrinos by optical detection of Cherenkov radiation induced by secondary charged leptons which are produced in neutrino interactions in the surrounding ice or the nearby bedrock.
Construction of the IceCube detector started in the Austral summer of 2004.In January 2010, 79 detec-  tor strings with 60 digital optical modules (DOMs) on each string [32,33] were deployed.Each DOM contains a 25 cm Hamamatsu photomultiplier tube (PMT) and on-board electronics to readout and digitize the signal from the PMT [34].In December 2010, the Ice-Cube detector construction was completed.The final IceCube detector consists of 86 strings.A schematic layout of the detector is shown in Fig. 2. The 79 strings used in this analysis are marked by green and yellow markers within the outer shaded grey area.The square markers denote the additional strings constituting the completed IceCube detector.Of the 79 strings used in this analysis, 73 strings (green) have a horizontal spacing of 125 m and a vertical spacing of 17 m between DOMs.The six remaining strings (yellow) are located near the central string of IceCube.The DOMs on these strings are equipped with PMTs with a 30% increased quantum efficiency.Together with their nearest IceCube strings, these strings constitute the inner detector, DeepCore [19] (for IceCube-79).The vertical distance between DOMs is reduced to 7 m (10 m) for the bottom 50 (upper 10) DOMs.The horizontal distance between strings in DeepCore is less than 75 m.
These two densely-instrumented parts (see Fig. 2) are separated by a region with significantly reduced scattering and absorption lengths for Cherenkov pho-tons due to dust particles.It is located at a depth of about 2050 m.
For muon-neutrino events, the neutrino arrival directions are inferred from the muon arrival direction.The latter is reconstructed using a likelihood approach, based on the arrival times of photons at DOMs [35].Thus, a good understanding of the absorption and scattering of photons is necessary for direction reconstruction.The clean glacial ice is a natural medium, built up over tens of thousands of years.Thus, the optical properties exhibit a variation over the 1 km depth of the instrumented volume.A detailed description of the optical ice properties is given in [36].

Event Selection
The signal for this analysis is muons produced in chargedcurrent neutrino interactions.These muons produce tracklike event signatures in the detector, which allow for a reconstruction of the arrival direction.
At the South Pole, the Galactic Center is always 29 • above the horizon.Thus, neutrinos from the direction of the Galactic Center will appear as down-going events in IceCube.The backgrounds for this analysis are therefore down-going atmospheric muons and, at a lower rate, muons produced by neutrinos originating from cosmic-ray showers in the atmosphere.The overwhelming majority of the 2500 Hz trigger rate is due to atmospheric muons.The atmospheric neutrino background contributes to the trigger rate at the 1 mHzlevel.This event class is an irreducible background, since in the energy range of interest for this analysis the accompanying muon component of the atmospheric shower is absorbed in the ice sheet above the detector.
The approach adopted here to reduce the background is to consider only neutrinos which interact within the detector, and reject the background of penetrating (incoming) muons.In order to select events which appear to start within the detector we developed several complementary veto techniques, exploiting differences in timing and topology of background and signal events.
This work is based on two independently developed event selections, referred to as low-energy (LE) and high-energy (HE) selections or samples.There are two reasons for such energy-specific optimization.First, the efficiency of these vetoes decreases rapidly with decreasing event energy, because low-energy muons are able to traverse several string-layers without being detected.Second, the likelihood function (described in Section 6) does not use the event energy.
Though the individual event selections differ, the general selection techniques and the analysis pipelines  are very similar.First, an initial online selection is performed at the South Pole.Second, a set of cuts on the event quality, topology, and arrival directions is applied.Third, a boosted decision tree (BDT) is used to remove remaining background events, using the TMVA software package [37].The BDT is trained on a representative signal assumption for each sample.Finally, a likelihood analysis is performed, exploiting the different distributions of arrival directions of background events and events originating in dark matter annihilations in the Galactic Center.The differences between the two event selections are highlighted in the following sections.
This analysis uses data collected with IceCube in its 79-string configuration between May 31, 2010 and May 13, 2011 with a total live-time of 319.7 days of stable high-quality data.The LE sample contains 35,538 events, and the HE sample contains 293,043 events.4,706 events appear in both samples; about 13% of the LE events are in the HE sample, and about 1.6% of the HE events are in the LE sample.

Low-Energy Event Selection
The LE event selection considers events from the Deep-Core online-filter [19], and is optimized for low-mass WIMPs below 100 GeV, and thus uses the bottom part of the densely-instrumented DeepCore sub-array as fiducial volume.The remaining instrumented IceCube volume as well as the two bottom DOM layers are used as a veto.The fiducial volume is illustrated in Fig. 2, and corresponds to roughly 27 Mton of ice.
The LE selection cuts are based on experience from the IceCube-79 Solar WIMP analysis [38], which used DeepCore for the first time in low-mass WIMP searches.The signal used for BDT training are events that are fully contained in the fiducial volume, and originate in annihilations of 65 GeV WIMPs to b b-pairs in the NFW halo.The search window for the LE analysis extends to ±30 • in right ascension (α) with respect to the Galactic Center, while the declination (δ) width is asymmetric and extends from −39 The LE sample data rate at the analysis level is 1.4 mHz.

High-Energy Event Selection
The IceCube array has a trigger energy threshold of about 100 GeV.Therefore, the search for WIMPs in the mass range above a few hundred GeV benefits from the large volume of IceCube in addition to Deep-Core at the cost of a decreasing veto efficiency.The HE selection considers events from the dedicated Galactic Center online-filter and the DeepCore online-filter.
The veto for the HE event selection is defined by the upper 12 DOM layers and the two outer string layers, which roughly corresponds to 200 m and 125 m of instrumented distance, respectively.The fiducial volume for the HE selection is shown in Fig. 2. The signal assumed for BDT training are events that are events which start in the fiducial volume, and originate in annihilations of 600 GeV WIMPs to W + W − -pairs in the NFW halo.The search window around the Galactic Center is given by ±15 • in both declination and right ascension.
The HE sample data rate at the analysis level is 10 mHz.

Analysis Method and Sensitivity
A maximum likelihood analysis is performed independently on each event selection for a number of different WIMP masses ranging from 30 GeV to 10 TeV, assuming a 100% branching ratio for each tested annihilation channel.
Considering the large number of events in the two final samples, a binned likelihood method was chosen.To reduce the number of bins, event arrival directions were only considered in a search window around the Galactic Center.The search window shape and size differ slightly between the two event samples as defined in the two previous sections 5.1 and 5.2 as well as the background estimations.Due to these differences the likelihood analysis performed on the LE and HE selections is not identical.The LE likelihood has the more complicated form and is defined as a function of signal fraction, ξ, in the following way: where the binomial factor in front of the product accounts for the probability of observing n events in the search window given N total events in the event selection, and the shape term, i.e the direction, X = (δ, α), of the events is accounted for by the term f (X i , |ξ).
The binomial probability was chosen in favor of a Poisson probability since the search window covers a non negligible fraction of the declination band.The probability of an event to fall in the search window is defined as p = π s ξ + π bg (1 − ξ), where π s and π bg are the probability for a signal or a background event, respectively, to fall in the search window.Note that the signal originating in a dark matter halo is an extreme case of an extended source as it is present in the whole sky.Any background estimation based on data will be contaminated by signal.π bg is determined from the relation between the size of the search window and the size of the background estimation region, while π s is determined from simulation.The directional probability density function (pdf), f (X|ξ), in Eq. ( 5) is constructed from binned expectations of event directions for background and signal.
Figure 3 shows examples of these binned expectations.The software package HEAL-Pix [39] was used to ensure equal-area bins on the sphere for these two-dimensional pdfs.To determine the signal pdfs, IceCube neutrino simulations were used, weighting events according to Eq. ( 4) and the corresponding DM annihilation spectrum.Background pdfs were created by scrambling the right ascension of experimental data events in the final event selections.
The same reasoning regarding signal contamination, as stated above, applies to the directional pdf, f (X|ξ).In effect, the expected arrival directions of background events will depend on the signal strength.This needs to be accounted for in the directional pdf, f (X|ξ) in the likelihood, as well as in the background simulation during the construction of confidence intervals.The directional pdf is defined as where f s and f bg are the signal and background directional pdfs, respectively, f sc is a pdf describing the signal scrambled in right ascension.The signal fraction inside the search window is A different approach is used for the HE analysis.The background estimation is performed on off-source data, excluding all events within ±30 • of the Galactic Center.Therefore, any signal contamination of the background estimate is ignored, and the likelihood function from Eq. ( 5) simplifies to: where n s is the number of signal events and n bg is the expected number of background events in the search window which is given by the number of off-source events multiplied by the ratio of the on-source and off-source region sizes.The directional pdf consequently becomes: All confidence intervals are constructed using the prescription by Feldman and Cousins [40].The sensitivity is defined as the median upper limit on the number of signal events at 90% confidence level.
The final limits for each WIMP mass and annihilation channel are obtained from the sample (LE or HE) which gives the best sensitivity (w/o systematics), i.e the lowest median upper limit, for the particular 10 -1 10 0 p.d.f.value per sr WIMP mass and annihilation channel.Thus, the crossover point in the WIMP mass between the two event samples depends on the DM halo model and WIMP annihilation channel.This procedure circumvents the necessity to deal with the small overlap of both samples.Figure 4 illustrates how the two event samples contribute to the best sensitivity at different WIMP masses, in this case for WIMP annihilation to ν ν.Each sample has a WIMP mass range where it outperforms the other.
Figure 5 shows the neutrino effective area for the two event selections.Even though the effective area for the LE event selection is smaller than that of the HE event selection, the sensitivity to the number of signal events for low-mass WIMPs is better, as can be seen in Fig. 4. The reasons are the larger on-source region for the LE event selection, and a lower background event rate due to higher veto efficiency.
In order to avoid confirmation bias throughout the development of the analysis, blindness with respect to the right-ascension information was imposed by scrambling the right-ascension information of the experimental data.Only the declination information of the events was used for cut development.
Following the optimization of the event selections, meaning all cuts are fixed, and the choice of the event selection for each channel, halo, and WIMP mass (based on best sensitivity), the right ascension information was unblinded.Although the HE effective area is bigger than the LE effective area at low energies, the higher background contamination at low energies in the HE selection makes it less efficient.

Discussion of Uncertainties
The uncertainties relevant for this analysis can be categorized into two classes: -Detector systematic uncertainties impacting the signal efficiency -Astrophysical uncertainties (choice of halo model, model-specific parameter) The former are incorporated into the calculated limits, while the latter are studied to estimate model uncertainties.Both classes are discussed in the following sections.

Detector Systematics
The uncertainties in the signal efficiency are mainly governed by the uncertainties in the optical efficiency of the DOMs and the optical properties of the glacial ice, manifested in the absorption and scattering length.
To determine the effects on σ A v due to the mentioned uncertainties, the event selections and analysis were applied to sets of simulated data where the optical properties of the DOMs and the ice were changed.The optical efficiency of the DOMs was varied by ±10%.The same was done for the absorption and scattering lengths of the ice.The resulting uncertainties on σ A v generally lie in the range 10% -20% except for the lowest neutrino energies where they reach up to ≈ 70%.This is due to threshold effects where events with just enough hit DOMs to trigger the detector would fail to do so with increased absorption and scattering or decreased optical efficiency.The above-described systematic uncertainties are included into the limits by degrading the baseline results by the relative variation of the detector uncertainties with respect to the baseline, as stated above.

Astrophysical Uncertainties
The astrophysical uncertainty is studied by using two different halo profiles, and also by varying the parameters ρ local and r s within the uncertainties stated in [26], and summarized in Tab. 1. Figure 6 compares the sensitivity for WIMPs annihilating to ν ν-pairs for both profiles.The bands depict the variation of the sensitivity within each profile that arises from varying the profile parameters within the given uncertainty.
The relative variation of σ A v due to the halo profile parameter uncertainties was estimated to be 60% -100% for the LE likelihood analysis and 60% -200% for the HE likelihood analysis, and is shown in Fig. 6.However, we refrain from including uncertainties on the halo profile parameters into the limits.

Results
Table 2 shows the number of events for the two unblinded event selections.The quantities n obs and n bg are the number of measured on-source events and expected background events in the search windows, respectively.A 2σ under-fluctuation of experimental data events was observed for the LE event selection.A systematic origin of this under-fluctuation due to an uneven right-ascension exposure was excluded.
A small over-fluctuation was measured for the HE event selection.However, after applying the likelihood analysis to both event selections independently, all the resulting upper limits were smaller than their corresponding sensitivity.For the HE selection this implies that despite the over-fluctuation in the number of events, the spatial distribution of these events within the search window is incompatible with the expectation from dark matter annihilation in the halo.As an example, Fig. 7 shows the sensitivity and the observed upper limit after unblinding on σ A v for the NFW and Burkert profiles, assuming WIMP annihilation into neutrinos.In addition the ±1σ and ±2σ statistical uncertainty bands of the median upper limit are shown as green and yellow shaded areas, respectively.The contribution from the LE and HE event selection can be clearly seen through the upper limit curve (solid black line) being lower than the expected median upper limit (dashed black line) in both cases, but to a different extent, with the switchover between m χ = 100 GeV and 200 GeV for this annihilation channel.3 summarizes the upper limits for all considered annihilation channels and WIMP masses.The limits are shown separately for the two considered DM halo profiles; the NFW (top) and Burkert (bottom) profile.Figure 8 shows the sensitivity and limit for three different annihilation channels.
To compare the performance of this analysis to previous IceCube analyses and other experiments, we choose the τ + τ − annihilation channel assuming a NFW DM halo profile.The comparison is shown in Fig. 9.The black solid line shows the limit of this analysis, whereas dashed lines with markers show the limits from previous galactic halo [16,17] and dwarf spheroidal galaxies [41] analyses with IceCube.The other lines show the limits from gamma-ray experiments, in particular the limit from the dwarf spheroidal galaxy Segue 1 analysis by VERITAS [42] (dash-dotted) and MAGIC [43] (dash-dot-dotted), and the limit from the Fermi analysis of several dwarf spheroidal galaxies [44] (dashed).Also shown is the DM interpretation of the positronfraction excess reported by the PAMELA collaboration (dark gray shaded region) and the 3σ and 5σ preferred regions from the e + + e − -flux excess reported by the Fermi and H.E.S.S. collaborations as dark green and green shaded regions, respectively.All the shaded region data are taken from [45] and rescaled to a local dark matter density of ρ local = 0.471 GeV cm −3 to this DM halo profile parameter with the one considered in the other analyses.For a WIMP mass below 1TeV the present analysis improves significantly in sensitivity on previous IceCube analyses.Furthermore by using the DeepCore detector array, the self-annihilation crosssection for WIMP masses below 100 GeV is probed for the first time by IceCube.

Conclusion
We have presented limits on the cross-section on dark matter annihilation in the Galactic Center, probing down to σ A v 4 • 10 −24 cm 3 s −1 at 65 GeV WIMP mass, assuming the NFW halo profile and direct annihilation to neutrino pairs.This analysis is the first IceCube Galactic Center DM search using the nearly complete detector configuration.Further, it is the first IceCube DM search probing σ A v for WIMP masses below 100 GeV by utilizing the DeepCore infill-array of IceCube.
We have presented methods for a selection of downgoing muon neutrinos in IceCube, making the southern hemisphere accessible to low-energy neutrino searches in the energy range 10 GeV -10 TeV.These methods have been applied to create two event selections, that are optimized for neutrino signals from the direction of the Galactic Center.Based on these event selections a likelihood analysis looking for a neutrino flux from annihilating dark matter in the Galactic Center was performed, testing a number of dark matter annihilation channels at different masses.The results are compatible with the background-only hypothesis, thus upper limits on σ A v were set (c.f.Fig. 9).The limits from the low-energy selection are almost 2σ lower than their sensitivity due to an under-fluctuation in the number of background events.The limits presented here for direct annihilation to νν-pairs are model-independent and conservative upper bounds for dark matter annihilation to Standard Model final states [49]; even small branching ratios to other -more visible -species at the σ A v -level presented here would yield a detectable flux in gamma-ray experiments, or otherwise stronger constraints.Thus, these limits complement gamma-ray detection channels.
Future improvements to this analysis can be expected from improvements in the background rejection in the energy region corresponding to the highest probed WIMP masses, and the inclusion of an energy term in the likelihood function.
Long-term improvements should also be expected from possible IceCube extensions.The low-energy upgrade PINGU [50] would increase the sensitivity to lowmass WIMPs, and extend the probed mass range below 30 GeV.PINGU is a possible future in-fill array with a denser instrumentation than DeepCore.The highmass (TeV-PeV) sensitivity would benefit from a future high-energy extension, IceCube-Gen2 [51].The aim for IceCube-Gen2 is an expanded instrumented volume of the order of 10 km 3 with a larger inter-string spacing, compared to IceCube.Fig. 9 Comparison of limits from this work (IC79 GC) to other IceCube (designated by IC + "string number") searches for dark matter annihilation in self-bound structures [16,17,41].Further, photon search limits from observation of dwarf spheroidals by VERITAS [42], MAGIC [43] and Fermi [44] are shown.The grey-shaded region is a dark-matter interpretation of the positron excess reported by the PAMELA collaboration.The green-shaded regions are the 3σ and 5σ preferred regions from the e + + e − -flux excess reported by Fermi and H.E.S.S.All shaded region data taken from [45].The region data and the IC22 halo limits are rescaled to the here assumed local dark matter density of ρ local = 0.471 GeV cm −3 .The natural scale is the self-annihilation cross-section region for WIMPs to be thermal relics from the Big Bang [46].The black dotted line in the upper right part of the figure is the unitarity bound [47].We note that preliminary Galactic Center limits from the ANTARES neutrino telescope have recently been released [48].

Fig- ure 1 (
top panel) shows J a for the NFW (solid line) and Burkert (dashed line) profiles.The final differential neutrino flux from DM annihilation, dφ ν /dE, depends on the neutrino energy spectrum of the actual annihilation channel.The differential neutrino flux is

Fig. 2
Fig.2A footprint (top) and side view (bottom) of Ice-Cube in detector coordinates.The green circles mark regular IceCube DOMs with 17 m vertical spacing, while the yellow diamonds mark DeepCore DOMs with about half the regular vertical DOM spacing and higher quantum efficiency PMTs.The black squares mark DOMs which are part of the final 86-string configuration, but are not present in the here used 79-string configuration of IceCube.The purple shaded and red shaded areas illustrate the fiducial volumes used by the LE and HE event selections, respectively.

Fig. 3
Fig.3Example skymaps of the background (top) and signal (bottom) pdfs in equatorial coordinates for the LE event selection, for 100 GeV WIMPs annihilating into W + W − .The LE search window is marked with a dashed line.For illustration the search window for the HE event selection is marked with a dash-dotted line.

Fig. 4
Fig.4Median upper limits, i.e. sensitivities, at 90% C.L. (w/o systematics) for the two event selections assuming WIMP annihilation to ν ν for the NFW DM profile.The black solid line shows the combined best sensitivity for this particular annihilation channel and DM profile, considering both event selections.At 200 GeV the HE selection yields a slightly better sensitivity and is thus used here.

Fig. 5
Fig.5Neutrino effective areas as a function of energy for the two event selections.The effective areas for the LE and HE selection are shown as solid, and dotted lines, respectively.Although the HE effective area is bigger than the LE effective area at low energies, the higher background contamination at low energies in the HE selection makes it less efficient.

Fig. 6
Fig.6Uncertainties on the sensitivity due to the DM halo model parameter uncertainties for the NFW (dashed line, blue band) and Burkert (solid line, red band) DM profile, assuming WIMP annihilation to νν.The reduced width of the bands below 100 GeV is caused by different on-source regions, and thus differences in the integration of the J-factors.The dip below 100 GeV is caused by the under-fluctuation in the LE sample.

Fig. 7
Fig. 7 Sensitivity (dashed) and observed upper limit (solid, w/o systematics) at 90% C.L. for WIMPs annihilating to neutrinos assuming a NFW (top) and Burkert (bottom) DM halo profile.The statistical uncertainty on the sensitivity is shown at the 1σ (green band) and 2σ (yellow band) level.The dip below 100 GeV is caused by the under-fluctuation in the LE sample.

Fig. 8
Fig.8Sensitivity (dashed) and observed upper limits (solid) at 90% C.L., including detector systematics, for WIMPs annihilating to bb (stars), τ + τ − (squares), and directly to neutrinos (triangles) assuming a NFW DM halo profile.The shaded areas are guides for the reader's eyes connecting sensitivity and limit for a particular annihilation channel.

Table 2
Results in terms of number of events for the two event selections.The differences of observed and expected events in the on-source region, ∆n = n obs − n bg , show an under-and over-fluctuation of data for the LE and HE event selection, respectively.

Table 3
Final upper limits (including detector systematics) on the self-annihilation cross-section, σ A v , for different annihilation channels and WIMP masses, m χ , for the NFW (top) and Burkert (bottom) DM halo profiles.