Limits on Kaluza–Klein dark matter annihilation in the Sun from recent IceCube results

We interpret recent IceCube results on searches for dark matter accumulated in the Sun in terms of the lightest Kaluza–Klein excitation (assumed here to be the Kaluza–Klein photon, B1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B^1$$\end{document}), obtaining improved limits on the annihilation rate in the Sun, the resulting neutrino flux at the Earth and on the B1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B^1$$\end{document}-proton cross-sections, for B1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B^1$$\end{document} masses in the range 30–3000 GeV. These results improve previous results from IceCube in its 22-string configuration by up to an order of magnitude, depending on mass, but also extend the results to B1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B^1$$\end{document} masses as low as 30 GeV.


Introduction
There are many astrophysical and cosmological observations that point to the existence of a dark matter component as a key constituent of the Universe. Constraints on the amount of baryons in the Universe from CMB measurements, from measurements of the abundance of primordial light elements and from searches for dark objects using microlensing have practically ruled out the possibility that dark matter consists of known Standard Model particles [1]. In the most popular picture, dark matter is composed by nonrelativistic Weakly Interacting Massive Particles (WIMPs) of yet unknown nature [2]. Among the many theories beyond the Standard Model of particle physics that predict new particles that could be viable dark matter candidates, Kaluza-Klein type models [3,4] with universal extra dimensions (UED) 1 provide a WIMP in the Kaluza-Klein photon (B 1 ), the first excitation of the scalar gauge boson in the theory [5][6][7][8][9][10]. Usually denoted as the lightest Kaluza-Klein particle (LKP), it can have a mass in the range from a few hundred GeV (limit from rare decay processes [11,12]) to a few TeV (to avoid 1 Universal in this context meaning that all the Standard Model fields are free to propagate also in the new dimension. a e-mail: cph@physics.uu.se (corresponding author) overclosing the Universe), the mass being proportional to 1/R, where R is the size of the extra dimension.
LKPs in the galactic halo will suffer the same fate as any other WIMP. Assuming a non-zero B 1 -proton scattering cross section, B 1 's in the galactic halo with Sun-crossing orbits can loose energy through interaction with the matter in the Sun and eventually sink into its core, where they would accumulate, thermalize, and annihilate into Standard Model particles [13][14][15][16], which in turn can lead to a detectable neutrino flux. Neutrino telescopes like AMANDA, IceCube, ANTARES and Baikal have searched for signatures of dark matter in this way, mainly focusing on the SUSY neutralino as WIMP candidate [17][18][19][20]. Additionally, both IceCube and ANTARES have set limits to the spin-dependent LKP-proton cross section [21][22][23]. The IceCube limits are based on the event selection and analysis searching for WIMP dark matter performed in Ref. [24].
In this letter we use the latest dark matter searches from the Sun by IceCube [18] to improve the limits on the Kaluza-Klein photon cross section with protons. Given the increase in detector size (86 strings versus 22 in the previous IceCube analysis) and lifetime (104 days versus 532 days) the results presented in this letter improve those in Ref. [21] by up to an order of magnitude. Additionally, the presence of the lowenergy subdetector DeepCore allows to lower the explored LKP mass down to 30 GeV, compared to 250 GeV, the lowest mass studied in Ref. [21]. Table 1 shows the theoretical branching ratios of the B 1 selfannihilation processes in terms of the quark splitting mass: where m q 1 is the first fermion excitation in the Kaluza-Klein theory within universal extra dimensions [6]. The mass at which the LKP can be a good dark matter candidate depends on this splitting, which drives the co-annihilation with the higher KK modes which, in turn, determines the relic abundance [25]. Among the final annihilation products, only the weakly interacting neutrinos are able to escape from the Sun and be detected in neutrino observatories on Earth. Neutrinos resulting from these annihilations are also expected to be easily distinguishable from thermonuclear reaction products because the value of B 1 mass needed for it to be a good dark matter candidate (m B 1 > O100 GeV). 2 Assuming that the capture and annihilation rates of B 1 particles in the Sun, Γ C and Γ A , have reached an equilibrium [6,26], the relationship between them can be written as Γ A = 1 2 Γ C , where the one half factor comes from the fact that one annihilation requires two captured B 1 . This means that, assuming a particular velocity distribution for dark matter in the solar neighborhood and a solar structure model, the B 1proton scattering cross section, which drives the capture, is proportional to the annihilation rate,

LKP signatures from the Sun
where the proportionality constant λ i m B 1 depends on the mass of the B 1 , and the superscript i can take the values SD (for spin-dependent cross section) or SI (for spinindependent cross section). The neutrinos produced in the core of the Sun have to be numerically propagated to a detector on Earth to predict the detectable neutrino flux per unit area and time in the detector, Φ ν , which is proportional to the annihilation rate, Φ ν = η m B 1 Γ A . Such neutrino propagation must also take into account the solar composition, neutrino interaction processes with matter (such as absorption, re-emission, decays of secondary particles into neutrinos, etc.) and neutrino oscillations.
Note that the efficiency of LKP capture by the Sun is sensitive to the low-velocity tail of the, really unknown, An extra peak at the end of the spectra produced by neutrinos that escape the Sun without interacting has been omitted from the plot for the sake of legibility, but it is considered in the calculations. Note the normalization in relative energy units LKP velocity distribution in the galaxy, f (v) (lower velocity particles fall easier below the escape velocity of the Sun after collisions with solar matter). We assume here a standard Maxwellian dark matter velocity distribution. Although there is accumulating evidence for a non-Maxwellian component in f (v) from recent Gaia and Sloan Digital Sky Survey (SDSS) data [28], the modifications to the Maxwellian distribution affect mainly the high-velocity tail, and are expected from simple kinematics to bear limited consequences for capture in the Sun. See Ref. [29] for a detailed discussion on the effects of astrophysical uncertainties on the calculation of dark matter capture in the Sun.
We simulated one million annihilation events at the core of the Sun and we propagated the neutrinos to a detector in the ice at latitude 90 • S during the austral winter using WimpSim [27], for different assumed values of m B 1 . The simulations themselves do not include any information on the capture conditions or Γ A , but they do require a solar structure model as an input for the propagation of neutrinos, in this case the one from Ref. [30]. Figure 1 shows the energy distribution of the muon (plus anti-muon) neutrino fluxes at the Earth per unit area A and per annihilation in the Sun, as a function of reduced energy z = E ν /m B 1 , for some values of m B 1 . As expected, heavier B 1 particles produce steeper profiles in energy because of absorption of high-energy neutrinos on their way out of the Sun. The convolution of those spectra with the effective area of the detector gives a prediction of the number of events μ s per number of annihilations N A , expected at the detector . The bump at m B 1 = 100 GeV in the former is due to the transition from DeepCore resolution to IceCube resolution [18]. Values calculated at the discrete mases denoted by the dots. Lines are a linear interpolation between the dots to guide the eye In order to take into account the finite angular resolution of the detector, the direction of each event has been randomly smeared using the median angular resolution of Ice-Cube shown in Fig. 2, which has been computed from the median muon neutrino energy at every m B 1 and the corresponding mean angular separation between the original neutrino and the reconstructed muon trajectories (see Ref. [18] for the explicit relation). Θ is used to spread the predicted signal (4) with a 2D Gaussian distribution centered around the Sun on the celestial sphere: where ψ is the angle from the Sun (ψ = 0), and C a normalization constant. The expected angular distribution of the signal, for a few B 1 masses as a function of angle with respect to the Sun position is shown in Fig. 5.

Data selection
We use the muon-neutrino effective area corresponding to the latest solar dark matter search from IceCube [18], shown in Fig. 3. The figure shows the effective areas for IceCube and DeepCore for muon neutrinos as a function of neutrino energy, E ν , for the case of neutrinos arriving near the horizon, since the Sun is always close to the horizon in the South Pole. For this analysis, both curves are summed and the entire detector is treated as one single array. Only muon events are considered here since the long muon tracks allow pointing Fig. 3 The combined muon and anti-muon neutrino effective area for both the DeepCore (blue) and main IceCube array (red) for horizontally arriving neutrinos as a function of the neutrino energy. Uncertainties are not shown, but they can be as large as 30%. Source: Ref. [18] back to the Sun. Muon in what follows refer to both muons and anti-muons since IceCube can not distinguish between particles and antiparticles. The predicted signal (Eq. 4) can be compared with the actual observed signal rate to compute Γ A : and Γ A is used to calculate the B 1 -proton scattering cross sections through Eq. (2). The conversion factors λ i m B 1 are obtained with the DarkSUSY software [31], and include the information on the solar structure and the dark matter velocity distribution in the solar neighborhood [26]. Figure 4 shows the expected number of events at the detector per annihilation as a function of B 1 mass calculated with Eq. (4) We use the data from Fig. 6 in Ref. [18], which we present here combined for DeepCore and IceCube in Fig. 5, in order to extract limits on the spin-dependent B 1 -proton cross section as a function of B 1 mass. The data set was obtained during 532 days of exposure between May 2011 and May 2014. As thoroughly explained in Ref. [18], several filters were applied to the data sample in order to minimize the presence of background. Data was only taken into account if measured during austral winters, when the Sun is below the horizon, in order to avoid overlap with atmospheric muons originating from cosmic-ray induced showers in the atmosphere above the detector. For this reason, only muons with upward trajectories were selected. Still, atmospheric muon neutrinos created at any declination can cross the Earth without interacting and reach the detector from below, being an irreducible background. Boosted decision trees were used to maximize signal separation and reduce background.  [24], which was the basis for the previous Kaluza-Klein analysis [21]. Values calculated at the discrete mases denoted by the dots. Lines are just a linear interpolation between the dots to guide the eye 4 Results Figure 5 shows that no statistically significant deviation from the expected background was detected in the IceCube solar analysis, so a 90% confidence level limit on the B 1 -proton cross section can be extracted from the data by setting a limit to Γ A . The amount of expected signal events per unit of time can be estimated as dμ s dt μ s τ (7) where τ is the total exposure time, 532 days in our case. We proceed with a "counting experiment", comparing the number of background events to the number of observed events extracted from Fig. 5, and construct Poisson confidence intervals for the signal strength μ s , calculated with the Neyman method using the algorithm developed in Ref. [32], including the detector systematic uncertainties from Ref. [18]. A limit on μ s is easily translated to a limit on Γ A through equation (6) and further to limits on σ S I and σ S D through equation (2). Table 2 shows the results for a series of B 1 masses. Limits at 90% confidence level on the signal strength (μ s ), the annihilation rate in the Sun (Γ A ), the muon flux at the detector above 1 GeV (Φ μ ) and the spin-independent and spin-dependent cross sections (σ S I and σ S D ) are shown, along with the median angular resolution (Δθ ν ) and mean muon energy at the detector E μ for each signal model. Figure 6 shows σ S D versus B 1 mass for the current analysis (blue dots) and the previously published analysis by Ice-Cube in the 22-string configuration (orange dots), where it  [18] can be seen that the constraints have been improved by up to one order of magnitude. The figure also shows the results from the ANTARES collaboration [23] (black curve). The shaded area shows the disfavoured mass region for the first Kaluza-Klein excitation obtained from searches for UED at the LHC [33], where a limit on 1/R (GeV) is obtained by combining several searches for events with large missing transverse momentum or monojets by ATLAS and CMS at 8 TeV and 13 TeV center of mass energies. Collider searches provide a complementary approach to indirect searches for dark matter in the form of Kaluza-Klein modes with neutrino telescopes, being competitive in different regions of the LKP mass range. Additional constraints from cosmology (that the B 1 must have a relic density compatible with the estimated dark matter density from CMB measurements) require the mass of the B 1 to be below ∼ 1.6 TeV [9,25,35]. Thus, taken all results together, the allowed parameter space for the B 1 to constitute the only component of dark matter in the Universe is currently quite restricted, but non-minimal UED models, not probed here, can still provide viable dark matter candidates [10].
From the experimental point of view, a few simplifying assumptions have been made to obtain the presented results. Uncertainties on the solar structure model and dark matter velocity distribution have been ignored, and binned data (from Fig. 5) have been used instead of a continuous sample of individual events, which limits the statistical power of the analysis. Even with these approximations, the IceCube limit presented in this letter consid- Table 2 90% confidence level upper limits on the signal, the annihilation rate in the Sun, the muon flux at the detector above 1 GeV and the spin independent and spin dependent B 1 -proton cross sections for several values of the B 1 mass. The last two columns show the median angular resolution and the muon mean energy at the detector  [21]. Black line: limits from ANTARES [23]. Limits from IceCube have been evaluated at the discrete masses denoted by the dots. Lines are just a linear interpolation between the dots to guide the eye. The shaded area represents the disfavoured region (at 95% confidence level) on the mass of the LKP from the LHC [33] erably improves over those published in Refs. [21] and [23].

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: This article is based on data presented in Ref. [19].] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, pro-vide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .