## Searches for small-scale anisotropies from neutrino point sources with three years of IceCube data

,

M. G. Aartsen , M. Ackermann b as , J. Adams , J. A. Aguilar o w , M. Ahlers ab , M. Ahrens aj , D. Altmann , T. Anderson v ap , C. Arguelles ab , T. C. Arlen ap , J. Auffenberg , X. Bai a ah , S. W. Barwick , V. Baum y ac , J. J. Beatty q,r , J. Becker Tjus , K.-H. Becker j ar , S. BenZvi ab , P. Berghaus as , D. Berley , p E. Bernardini as , A. Bernhard ae, ∗ , D. Z. Besson , G. Binder z h,g , D. Bindig ar , M. Bissok , E. Blaufuss , J. Blumenthal , D. J. Boersma a p a aq , C. Bohm aj , F. Bos , D. Bose j al , S. B¨ser , O. Botner o k aq , L. Brayeur m , H.-P. Bretz as , A. M. Brown , J. Casey , M. Casier o e m , E. Cheung , D. Chirkin p ab , A. Christov w , B. Christy , K. Clark p am , L. Classen , F. Clevermann , v t S. Coenders ae , D. F. Cowen ap,ao , A. H. Cruz Silva as , M. Danninger aj , J. Daughhetee , J. C. Davis , M. Day e q ab , J. P. A. M. de Andr´ e ap , C. De Clercq m , S. De Ridder , P. Desiati x ab , K. D. de Vries m , M. de With , i T. DeYoung ap , J. C. D´ ıaz-V´lez e ab , M. Dunkman ap , R. Eagan ap , B. Eberhardt ac , B. Eichmann , J. Eisch j ab , S. Euler aq , P. A. Evenson af , O. Fadiran ab , A. R. Fazely , A. Fedynitch , J. Feintzeig f j ab , J. Felde p , T. Feusels , K. Filimonov , C. Finley x g aj , T. Fischer-Wasels ar , S. Flis aj , A. Franckowiak , K. Frantzen , T. Fuchs , T. K. Gaisser k t t af , R. Gaior , n J. Gallagher aa , L. Gerhardt h,g , D. Gier , L. Gladstone a ab , T. Gl¨senkamp u as , A. Goldschmidt , G. Golup h m , J. G. Gonzalez af , J. A. Goodman , D. G´ra p o as D. Grant , P. Gretskov , J. C. Groh u a ap , A. Groß ae , C. Ha h,g , C. Haack , a A. Haj Ismail , P. Hallen , A. Hallgren x a aq , F. Halzen ab , K. Hanson , l D. Hebecker , D. Heereman , D. Heinen , K. Helbing k l a ar , R. Hellauer , p D. Hellwig , S. Hickford , G. C. Hill a o b , K. D. Hoffman , R. Hoffmann p ar , A. Homeier , K. Hoshina k ab,1 , F. Huang ap , W. Huelsnitz , P. O. Hulth p aj , K. Hultqvist aj , S. Hussain af , A. Ishihara , E. Jacobi n as , J. Jacobsen ab , K. Jagielski , G. S. Japaridze a d , K. Jero ab , O. Jlelati x , M. Jurkovic ae , B. Kaminsky as , A. Kappes , T. Karg v as , A. Karle ab , M. Kauer ab , J. L. Kelley ab , A. Kheirandish ab , J. Kiryluk ak , J. Kl¨s a ar , S. R. Klein h,g , J.-H. K¨hne , G. Kohnen o t ad , H. Kolanoski , A. Koob , L. K¨pke i a o ac , C. Kopper ab , S. Kopper ar , D. J. Koskinen , M. Kowalski , A. Kriesten , s k a K. Krings , G. Kroll a ac , M. Kroll , J. Kunnen j m , N. Kurahashi ab , T. Kuwabara , M. Labare , D. T. Larsen n x ab , M. J. Larson , s M. Lesiak-Bzdak ak , M. Leuermann a, ∗ , J. Leute ae , J. L¨nemann u ac , J. Madsen ai , G. Maggi m , R. Maruyama ab , K. Mase , H. S. Matis , n h

R. Maunu , F. McNally p ab , K. Meagher , M. Medici , A. Meli , T. Meures , p s x l S. Miarecki h,g , E. Middell as , E. Middlemas ab , N. Milke , J. Miller t m , L. Mohrmann as , T. Montaruli w , R. Morse ab , R. Nahnhauer as , U. Naumann ar , H. Niederhausen ak , S. C. Nowicki , D. R. Nygren , u h A. Obertacke ar , S. Odrowski , A. Olivas , A. Omairat u p ar , A. O'Murchadha , l T. Palczewski an , L. Paul , O. Penek , J. A. Pepper a ¨ a an , C. P´ erez de los Heros aq , C. Pfendner , D. Pieloth , E. Pinat , J. Posselt q t l ar , P. B. Price , G. T. Przybylski g h , J. P¨tz , M. Quinnan u a ap , L. R¨del , a a M. Rameez w , K. Rawlins , P. Redl , I. Rees c p ab , R. Reimann , M. Relich , a n E. Resconi ae , W. Rhode , M. Richman , B. Riedel t p ab , S. Robertson , b J. P. Rodrigues ab , M. Rongen , C. Rott a al , T. Ruhe , B. Ruzybayev t af , D. Ryckbosch , S. M. Saba , H.-G. Sander x j ac , J. Sandroos , M. Santander s ab , S. Sarkar s,ag , K. Schatto ac , F. Scheriau , T. Schmidt , M. Schmitz , t p t S. Schoenen , S. Sch¨neberg , A. Sch¨nwald a o j o as , A. Schukraft , L. Schulte , a k O. Schulz ae , D. Seckel af , Y. Sestayo ae , S. Seunarine ai , R. Shanidze as , M. W. E. Smith ap , D. Soldin ar , G. M. Spiczak ai , C. Spiering as , M. Stamatikos q,2 , T. Stanev af , N. A. Stanisha ap , A. Stasik , k T. Stezelberger h , R. G. Stokstad , A. St¨ßl h o as , E. A. Strahler m , R. Str¨m o aq , N. L. Strotjohann , G. W. Sullivan , H. Taavola k p aq , I. Taboada , e A. Tamburro af , A. Tepe ar , S. Ter-Antonyan , A. Terliuk f as , G. Teˇ si´ c ap , S. Tilav af , P. A. Toale an , M. N. Tobin ab , D. Tosi ab , M. Tselengidou , v E. Unger , M. Usner , S. Vallecorsa j k w , N. van Eijndhoven m , J. Vandenbroucke ab , J. van Santen ab , M. Vehring , M. Voge , M. Vraeghe , a k x C. Walck aj , M. Wallraff , Ch. Weaver a ab , M. Wellons ab , C. Wendt ab , S. Westerhoff ab , B. J. Whelan , N. Whitehorn b ab , C. Wichary , K. Wiebe a ac , C. H. Wiebusch , D. R. Williams a an , H. Wissing , M. Wolf p aj , T. R. Wood , u K. Woschnagg , D. L. Xu g an , X. W. Xu , J. P. Yanez f as , G. Yodh , y S. Yoshida , P. Zarzhitsky n an , J. Ziemann , S. Zierke , M. Zoll t a aj

- c
- a III. Physikalisches Institut, RWTH Aachen University, D-52056 Aachen, Germany b School of Chemistry &amp; Physics, University of Adelaide, Adelaide SA, 5005 Australia Dept. of Physics and Astronomy, University of Alaska Anchorage, 3211 Providence Dr., Anchorage, AK 99508, USA

d CTSPS, Clark-Atlanta University, Atlanta, GA 30314, USA e School of Physics and Center for Relativistic Astrophysics, Georgia Institute of Technology, Atlanta, GA 30332, USA

f Dept. of Physics, Southern University, Baton Rouge, LA 70813, USA

g

Dept. of Physics, University of California, Berkeley, CA 94720, USA

h

Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA

i Institut f¨r Physik, Humboldt-Universit¨t zu Berlin, D-12489 Berlin, Germany u a

| j Fakult¨t a f¨r u Physik & Astronomie,                                                                                                                                                                                                                                                                         | Ruhr-Universit¨t a Bochum, D-44780 Bochum, Germany                                                                                                                                                                                                                                                              |
|-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| k Physikalisches Institut, Universit¨t a Bonn, l Universit´ e Libre de Bruxelles, Science m Vrije Universiteit Brussel, Dienst n Dept. of Physics, Chiba                                                                                                                                                        | Nussallee 12, D-53115 Bonn, Germany Faculty CP230, B-1050 Brussels, Belgium ELEM, B-1050 Brussels, Belgium University, Chiba 263-8522, Japan                                                                                                                                                                    |
| o Dept. of Physics and Astronomy, Christchurch,                                                                                                                                                                                                                                                                 | University of Canterbury, Private Bag 4800, New Zealand                                                                                                                                                                                                                                                         |
| p Dept. of Physics, University of Maryland, College Park, MD 20742, USA                                                                                                                                                                                                                                         | p Dept. of Physics, University of Maryland, College Park, MD 20742, USA                                                                                                                                                                                                                                         |
| q Dept. of Physics and Center for Cosmology and Astro-Particle Physics, Ohio State University, Columbus, OH 43210, USA                                                                                                                                                                                          | q Dept. of Physics and Center for Cosmology and Astro-Particle Physics, Ohio State University, Columbus, OH 43210, USA                                                                                                                                                                                          |
| r Dept. of Astronomy, Ohio State University, Columbus, OH 43210, USA                                                                                                                                                                                                                                            | r Dept. of Astronomy, Ohio State University, Columbus, OH 43210, USA                                                                                                                                                                                                                                            |
| s Niels Bohr Institute, University of Copenhagen, DK-2100 Copenhagen, Denmark t Dept. of Physics, TU Dortmund University, D-44221 Dortmund, Germany u Dept. of Physics, University of Alberta, Edmonton, Alberta, Canada T6G 2E1 v Erlangen Centre for Astroparticle Physics, Friedrich-Alexander-Universit¨t a | s Niels Bohr Institute, University of Copenhagen, DK-2100 Copenhagen, Denmark t Dept. of Physics, TU Dortmund University, D-44221 Dortmund, Germany u Dept. of Physics, University of Alberta, Edmonton, Alberta, Canada T6G 2E1 v Erlangen Centre for Astroparticle Physics, Friedrich-Alexander-Universit¨t a |
| Erlangen-N¨rnberg, u D-91058 Erlangen, Germany                                                                                                                                                                                                                                                                  | Erlangen-N¨rnberg, u D-91058 Erlangen, Germany                                                                                                                                                                                                                                                                  |
| w D´partement e de physique nucl´aire e et corpusculaire, Universit´ e de Gen`ve, e CH-1211 Gen`ve, e Switzerland                                                                                                                                                                                               | w D´partement e de physique nucl´aire e et corpusculaire, Universit´ e de Gen`ve, e CH-1211 Gen`ve, e Switzerland                                                                                                                                                                                               |
| x Dept. of Physics and Astronomy, University of Gent, B-9000 Gent, Belgium y Dept. of Physics and Astronomy, University of California, Irvine, CA 92697, USA                                                                                                                                                    | x Dept. of Physics and Astronomy, University of Gent, B-9000 Gent, Belgium y Dept. of Physics and Astronomy, University of California, Irvine, CA 92697, USA                                                                                                                                                    |
| z Dept. of Physics and Astronomy, University of Kansas, Lawrence, KS 66045, USA                                                                                                                                                                                                                                 | z Dept. of Physics and Astronomy, University of Kansas, Lawrence, KS 66045, USA                                                                                                                                                                                                                                 |
| aa Dept. of Astronomy, University of Wisconsin, Madison, WI 53706, USA ab                                                                                                                                                                                                                                       | aa Dept. of Astronomy, University of Wisconsin, Madison, WI 53706, USA ab                                                                                                                                                                                                                                       |
| Dept. of Physics and Wisconsin IceCube Particle Astrophysics Center, University of                                                                                                                                                                                                                              | Dept. of Physics and Wisconsin IceCube Particle Astrophysics Center, University of                                                                                                                                                                                                                              |
| Wisconsin, Madison, WI 53706, USA                                                                                                                                                                                                                                                                               | Wisconsin, Madison, WI 53706, USA                                                                                                                                                                                                                                                                               |
| Institute of Physics, University of Mainz, Staudinger Weg 7, D-55099 Mainz, Germany ad Universit´ e de Mons, 7000 Mons, Belgium ae                                                                                                                                                                              | Institute of Physics, University of Mainz, Staudinger Weg 7, D-55099 Mainz, Germany ad Universit´ e de Mons, 7000 Mons, Belgium ae                                                                                                                                                                              |
| Technische Universit¨t a M¨nchen, u D-85748 Garching, Germany Bartol Research Institute and Dept. of Physics and Astronomy, University of                                                                                                                                                                       | Technische Universit¨t a M¨nchen, u D-85748 Garching, Germany Bartol Research Institute and Dept. of Physics and Astronomy, University of                                                                                                                                                                       |
| Delaware,                                                                                                                                                                                                                                                                                                       | Delaware,                                                                                                                                                                                                                                                                                                       |
| Newark, DE 19716, USA                                                                                                                                                                                                                                                                                           | Newark, DE 19716, USA                                                                                                                                                                                                                                                                                           |
| ag Dept. of Physics, University of Oxford, 1 Keble Road, Oxford OX1 3NP, UK ah Physics Department, South Dakota School of Mines and Technology, Rapid City, SD                                                                                                                                                  | ag Dept. of Physics, University of Oxford, 1 Keble Road, Oxford OX1 3NP, UK ah Physics Department, South Dakota School of Mines and Technology, Rapid City, SD                                                                                                                                                  |
| 57701, USA ai Dept. of Physics, University of Wisconsin, River Falls, WI 54022, USA                                                                                                                                                                                                                             | 57701, USA ai Dept. of Physics, University of Wisconsin, River Falls, WI 54022, USA                                                                                                                                                                                                                             |
| Oskar                                                                                                                                                                                                                                                                                                           | Oskar                                                                                                                                                                                                                                                                                                           |
| aj                                                                                                                                                                                                                                                                                                              | aj                                                                                                                                                                                                                                                                                                              |
| Klein Centre and Dept. of Physics, Stockholm University, SE-10691 Stockholm,                                                                                                                                                                                                                                    | Klein Centre and Dept. of Physics, Stockholm University, SE-10691 Stockholm,                                                                                                                                                                                                                                    |
| Sweden ak Dept. of Physics and Astronomy, Stony Brook University, Stony Brook, NY 11794-3800, USA al                                                                                                                                                                                                            | Sweden ak Dept. of Physics and Astronomy, Stony Brook University, Stony Brook, NY 11794-3800, USA al                                                                                                                                                                                                            |
| Dept. of Physics, Sungkyunkwan University, Suwon 440-746, Korea                                                                                                                                                                                                                                                 | Dept. of Physics, Sungkyunkwan University, Suwon 440-746, Korea                                                                                                                                                                                                                                                 |
| Dept. of Physics, University of Toronto, Toronto, Ontario, Canada, M5S 1A7                                                                                                                                                                                                                                      | Dept. of Physics, University of Toronto, Toronto, Ontario, Canada, M5S 1A7                                                                                                                                                                                                                                      |
| am of Physics and Astronomy, University of Alabama, Tuscaloosa, AL 35487, USA                                                                                                                                                                                                                                   | am of Physics and Astronomy, University of Alabama, Tuscaloosa, AL 35487, USA                                                                                                                                                                                                                                   |
| Dept.                                                                                                                                                                                                                                                                                                           | Dept.                                                                                                                                                                                                                                                                                                           |
| an Dept. of Astronomy and Astrophysics, Pennsylvania State University, University Park,                                                                                                                                                                                                                         | an Dept. of Astronomy and Astrophysics, Pennsylvania State University, University Park,                                                                                                                                                                                                                         |
| PA 16802, USA                                                                                                                                                                                                                                                                                                   | PA 16802, USA                                                                                                                                                                                                                                                                                                   |
| ao                                                                                                                                                                                                                                                                                                              | ao                                                                                                                                                                                                                                                                                                              |

aq Dept. of Physics and Astronomy, Uppsala University, Box 516, S-75120 Uppsala, Sweden ar Dept. of Physics, University of Wuppertal, D-42119 Wuppertal, Germany as DESY, D-15735 Zeuthen, Germany

## Abstract

Recently, IceCube found evidence for a diffuse signal of astrophysical neutrinos in an energy range of ∼ 60 TeV to the PeV-scale [1]. The origin of those events, being a key to understanding the origin of cosmic rays, is still an unsolved question. So far, analyses have not succeeded to resolve the diffuse signal into point-like sources. Searches including a maximumlikelihood-ratio test, based on the reconstructed directions and energies of the detected down- and up-going neutrino candidates, were also performed on IceCube data leading to the exclusion of bright point sources. In this paper, we present two methods to search for faint neutrino point sources in three years of IceCube data, taken between 2008 and 2011. The first method is an autocorrelation test, applied separately to the northern and southern sky. The second method is a multipole analysis, which expands the measured data in the northern hemisphere into spherical harmonics and uses the resulting expansion coefficients to separate signal from background. With both methods, the results are consistent with the background expectation with a slightly more sparse spatial distribution, corresponding to an underfluctuation. Depending on the assumed number of sources, the resulting upper limit on the flux per source in the northern hemisphere for an E -2 energy spectrum ranges from ∼ 1 5 . · 10 -8 GeV/cm 2 s -1 , in the case of one assumed source, to ∼ 4 · 10 -10 GeV/cm 2 s -1 , in the case of 3500 assumed sources.

Keywords: Extraterrestrial Neutrinos, Astrophysical Neutrinos, Point Sources, IceCube, 2pt-Correlation, Multipole Analysis

∗ Corresponding authors: Anna.Bernhard@tum.de, Leuermann@physik.rwth-aachen.de

1 Earthquake Research Institute, University of Tokyo, Bunkyo, Tokyo 113-0032, Japan 2 NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA

## 1. Introduction

The unsolved problem of the origin of cosmic rays is one of the biggest challenges in high-energy astrophysics. In hadronic interactions of cosmic rays with matter, high-energy neutrinos are produced in the direct environment of cosmic ray sources. Possible candidates for these sources are, for example, Gamma Ray Bursts (GRBs) [2, 3] or Active galactic Nuclei (AGNs) [4, 5]. From Fermi Acceleration, the resulting energy-dependent differential flux ( d φ d E ν ) of astrophysical neutrinos is expected to follow a power law of E -γ ν , where E ν is the neutrino energy and γ = 2+ glyph[epsilon1] is close to the spectral index of the cosmic ray production process [6]. The so-called inefficiency ( ) is often glyph[epsilon1] assumed to be 0 ≤ glyph[epsilon1] ≤ 1, while γ ≈ 2 is favored by Fermi Acceleration [6]. Neutrinos are also expected in coincidence with high-energy photons. Thus, the measurement of thousands of gamma-ray sources, like AGNs [7, 8], provides a large number of potential neutrino sources.

Since neutrinos do not experience deflections or scattering, they are ideal messengers for observing and tracing the hadronic interactions described above. The detection of high-energy neutrinos of cosmic origin would give important insights for identifying the sources of cosmic rays [9, 10].

The IceCube neutrino telescope [11] recently detected a diffuse cosmic flux of high-energy neutrino events in two years of IceCube data by searching for neutrino-induced events with an interaction vertex within the detector [1]. This permits one to reduce the dominant atmospheric muon background. In a follow-up analysis, this diffuse flux was investigated in more detail using three years of IceCube data [12]. This high energy starting events analysis (HESE) yielded 37 events compared to an expected background of 15.0 events from atmospheric muons and neutrinos. These events could provide information about potential cosmic ray astrophysical sources and motivate additional searches.

IceCube, with its large field-of-view, offers a unique opportunity to study the production and interaction of high-energy cosmic rays using neutrinos. The detector, which is located at the geographical South Pole, has a detection volume of ∼ 1 km 3 deep in the Antarctic glacier and an additional ∼ 1 km 2 surface air shower detector, called IceTop. The IceCube detector consists of digital optical modules (DOMs) [13], placed on strings deployed vertically at depths between 1450 m and 2450 m. The strings hold 60 DOMs each equipped with a photomultiplier tube and digitizing electronics to detect neutrinos by measuring Cherenkov radiation of their secondary particles [14].

A detailed description of the data acquisition system can be found in [13]. The detector was built in several stages between 2005 and 2010, such that each year several strings were added, until reaching the final configuration of 86-strings, containing more than 5000 DOMs. IceCube has completed different point source searches, including an energy-dependent likelihood point source search scanning the full sky [15], as well as searches for flaring and periodic neutrino emission [16]. Additionally there are searches for diffuse neutrino emission looking for deviations in the two dimensional distribution of energy and zenith angle [17]. Point source searches are most sensitive for finding individual sources of astrophysical neutrinos among the background of atmospheric events (neutrinos and muons from cosmic ray interactions at Earth). Diffuse searches, on the other hand, are most sensitive for detecting within this background the presence of high-energy astrophysical neutrinos throughout the sky, without identifying individual sources. In between these two scenarios is the possibility that many weak sources exist. These could contribute to the detected diffuse signal and create a small number of events clustering on the background event distribution, while the individual clusters remain too weak to be detected by the point source searches. In this paper we present two searches for such small-scale clustering.

The first is an autocorrelation test using a two-point (2-pt) correlation function performed in the northern and southern hemisphere. Since most of the signal-like events are at high energies, we extended the autocorrelation test to include the most energetic events in an additional test (2-pt HE). The second test is a multipole expansion of the skymap of neutrino arrival directions. The goal of both methods is to gain sensitivity to faint sources at unknown positions in the sky with unidentified energy spectra using a three year data-set of the partially completed detector. Both methods are complementary to previous searches and therefore an important addition to IceCube searches for a cosmic neutrino flux. A similar search for clustering was recently performed by ANTARES, details can be found in [18].

The paper is organized as follows: In section 2, the data sample used and the generation of pseudo experiments is described. The 2-pt correlation test is explained in section 3.1, and the multipole analysis in section 3.2. The performance of both analyses, decribed by the discovery potential, is given in section 4. In section 5, the experimental result is presented and exclusion limits are calculated. Systematic uncertainties are discussed in section 6. In section 7 a conclusion is drawn.

## 2. Data Sample

## 2.1. Experimental Data

The data presented was taken between 2008 and 2011. It consists of 3 different detector configurations: 40-, 59- and 79-string configurations (IC40, IC59 and IC79, respectively), where the labeling corresponds to the number of deployed strings. For the selection criteria applied [15], 108 310 and 146 047 track-like events were found in the northern and southern sky, respectively. The data are background dominated, consisting mostly of atmospheric muons in the southern sky and of atmospheric neutrinos in the northern sky where the atmospheric muon background is shielded by the Earth. At trigger level, there is a significant contribution from atmospheric muons also in the northern sky, caused by mis-reconstructed atmospheric muons near the horizon and events with two muons within the same readout window reconstructed as a single track [19]. This background is reduced by rejecting events with a poor angular resolution and events with multiple muon tracks which are identified by a clustering algorithm. Requirements on the amount of light seen in the detector are used to also reduce the atmospheric neutrino background in the northern hemisphere and result in a neutrino energy threshold of around 100 GeV. Following these guidelines, the event selection was performed separately for each detector configuration. For IC40, straight cuts were used [20, 15], while for IC59 and IC79, the event selection was performed using several Boosted Decision Trees [15]. To reduce the atmospheric muon background in the southern hemisphere, declination dependent cuts on the reconstructed muon energy were developed for each data sample. While the cuts were designed to result in a smooth distribution of the events, there are remaining disuniformities induced by the steep change of the event density with the declination if a smaller binning is used to present the data than was used in the event selection. The effects of these disuniformities are fully included in the background estimation and the calculated statistical significance of the result. For each selected event, the declination, the right ascension and the energy of the muon inside the detector were reconstructed, using the methods described in [21] and [19]. Atmospheric neutrinos yield events with reconstructed muon energies of up to ∼ 100 TeV. Above this energy, their flux is too low to make a significant contribution. The atmospheric muon background in the southern hemisphere however is significantly larger and contributes events with up to 10 times higher energies. The declination distribution of the experimental data after the event selection is shown in

Fig. 1(a).

## 2.2. Signal Simulation

The energy spectrum of candidate neutrinos is motivated by Fermi acceleration [6], and is assumed to follow a power law of E -γ with different spectral indices ( γ ). Additionally, the number of sources ( N Sou ) is varied between one and several thousand sources with respect to the wide range of possible source candidates. From full detector Monte-Carlo simulation (MC) [22], the detection efficiency and the angular resolution of these signal events can be estimated. Since the angular resolution of the detector depends on the energy spectrum, the distribution of the angular reconstruction error, i.e. the Point Spread Function (PSF), is obtained for each of the investigated energy spectra ( E -2 , E -2 25 . and E -3 ) and for each detector configuration. As an example, for an E -2 spectrum, the PSFs for IC40, IC59 and IC79 are shown in Fig. 1(b).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0634v2-figure-dc0f03c5-2633-4387-8bfd-54e5446ae098.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0634v2-figure-de38ce7d-87c9-4d8c-a961-6e00aa3841c7.jpg)

- (a) Experimental declination distribution (b) E -2 Point Spread Function from MC

Figure 1: (a) Experimental declination distribution and (b) Point Spread Function (PSF) obtained from full detector MC simulations. In (b) the median and the 90%-quantile for each distribution are shown as solid and dashed vertical lines, respectively.

The detector efficiency is a function of the energy E and the zenith angle θ of an incoming neutrino. It can be characterized by the effective area A eff ( E,θ ), which is the corresponding area of a hypothetical detector, that has 100 % efficiency for detecting neutrinos generated uniformly and omnidirectionally at the surface of the atmosphere of Earth with given energy E and zenith angle θ . For a neutrino source at a given declination in the sky, the number of detected neutrinos is proportional to a convolution of the energy-dependent

Figure 2: The solid-angle-averaged effective area of IC40, IC59, and IC79 for neutrinos from the northern hemisphere (a) and the southern hemisphere (b). The effective area is determined from Monte Carlo simulations for the final event selection and is a function of the neutrino energy E ν and the incident zenith angle (not shown here, but averaged over solid angle).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0634v2-figure-53ce5ee2-9068-406c-abf8-bf37a16f6ef1.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0634v2-figure-72123615-e2c2-4466-8120-b9dd46e23763.jpg)

effective area (evaluated at the appropriate zenith angle) and the energy spectrum of the source. As an illustration, the solid-angle-averaged effective area A eff ( E ) for the northern and southern sky are shown in Fig. 2 for the final event selection.

The declination-dependent detector acceptance, as obtained from MC simulation, is shown in Fig. 3(a) for E -2 , E -2 25 . and E -3 energy spectra. Additionally, the detector acceptance for the top 10%, 1% and 0.1% high-energy events of an E -2 energy spectrum is shown in Fig. 3(b).

For the northern hemisphere, the expected neutrino flux from the pole region decreases at high energies which is due to the Earth's declining transparency to neutrinos. Thus, for hard energy spectra and the highest-energy bins the signal is dominated by the horizon region, while for soft energy spectra and low-energy bins, the expected number of signal neutrinos from an isotropic signal is largely the same at all declinations. In addition, the declinationacceptance is influenced by declination-dependent cuts, applied to the data sample in [15]. In the southern hemisphere, the detector acceptance is comparatively much smaller due to strong cuts applied in [15] to reduce the atmospheric muon background. The observed smooth transition at the horizon is due to the decrease of atmospheric muons and the corresponding increase in signal efficiency.

(sin((mc\_Zd-90)*3.141/180)) {(((mc\_Zd)&gt;0) *(OW*MC\_Neutrino\_En**(-2)/NEvents*86400*365))}

(sin((mc\_Zd-90)*3.141/180)) {(((log10(MC\_Neutrino\_En)&gt;3.8 &amp;&amp; mc\_Zd&gt;90) || (log10(MC\_Neutrino\_En)&gt;5.8 &amp;&amp; mc\_Zd&lt;90)) *(OW*MC\_Neutrino\_En**(-2)/NEvents*86400*365))}

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0634v2-figure-fed60c6a-aa38-4bfb-8bba-4cc66e366bd8.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0634v2-figure-a6cad060-91ce-4500-8b66-2c7ea7fd9e7d.jpg)

(b) For different energy bins of

E

-

2

Figure 3: Declination distribution δ for (a) different energy spectra and (b) different energy bins of an E -2 spectrum used as detector acceptance for the signal simulation. All histograms are obtained from MC simulation of isotropic signal events.

## 2.3. Simulation of Signal and Background Neutrino Sky Maps

Using toy MC simulations, pure atmospheric background sky maps, sky maps of pure point source signal and mixed sky maps are generated. These pseudoexperiments are used to develop the analysis method and to quantify the analysis performance (discovery potential). The generation of signal and background neutrinos is described in the following.

For all maps, the total number of events, n tot , is fixed to the number of events in the experimental dataset. For the signal simulations, three signal parameters are used: the mean number of neutrinos per source ( µ ), the number of sources ( N sou ) in the sky (for the corresponding hemisphere) and the spectral index ( γ ) of the signal energy spectrum. The spatial distribution of point sources is assumed to be uniform in the full sky. In an additional signal hypothesis, it is assumed to follow the galactic plane. The number of events produced by the source is drawn from a Poisson distribution with a mean source strength of µ . All source events are distributed around the true source position according to the PSF of the specified energy spectrum to take the detector's angular resolution into account. Each generated signal event is either rejected from or accepted to the simulated sky map according to its declination using the detector's declination-dependent acceptance in a hit-and-miss procedure [23]. Thus, µ is the number of generated neutrinos per source before detector acceptance and corresponds to a flux. The expected number of neutrinos per source ( µ eff ) in the simulated map depends

on the declination, such that µ eff is given by the product of µ and the detector acceptance at the source's declination. To do this, the MC declination distribution for the specified energy spectrum is normalized such that its peak value is at 1 0 and taken to be the detector's declination acceptance. . Thus, the number of neutrinos from each source within one sky map is determined by the declination of the source and Poissonian fluctuations.

The background, consisting of atmospheric neutrinos and mis-reconstructed atmospheric muons, is generated by performing a large number of pseudoexperiments using a uniform right ascension distribution and taking into account the declination distribution from the experimental data shown in Fig. 1(a). Using this data driven approach, no additional systematic effects due to Monte-Carlo simulation are introduced.

## 3. Methods

## 3.1. Two-Point Autocorrelation Test

The autocorrelation test is based on the distribution of pairwise calculated spatial distances (Ψ ij ) between events i and j , which are compared to the background expectation. It is designed to detect an event clustering at angular scales θ comparable to the detector resolution, while no prior information about the potential sources is required. The amount of clustering is then obtained from scanning simultaneously over θ and different bins for the minimum energy E min that optimizes the sensitivity [24, 25].

With Ψ being the spatial distance between two events, the test statistic of the 2-pt analysis can be defined as a function of θ by

$$T S ( \theta, E _ { \min } ) = \frac { \text{obs. no. pairs with $\Psi_{i,j}\leq \theta,E_{i,j}\geq E_{min}$} } { \text{avg. no. bg. pairs with $\Psi_{i,j}\leq \theta,E_{i,j}\geq E_{min}$} } \,, \quad ( 1 )$$

or, more precisely, by

$$T S ( \theta, E _ { \min } ) = \frac { \sum _ { i, j \in H, i > j } \Theta ( \theta - \Psi _ { i j } ) \cdot \Theta ( E _ { i, j } - E _ { \min } ) } { \langle \sum _ { m, n \in H, m > n } \Theta ( \theta - \Psi _ { m n } ) \cdot \Theta ( E _ { m, n } - E _ { \min } ) \rangle _ { b g } } \,, \quad ( 2 )$$

where E i is the reconstructed muon energy of event i (see [19] for details on the energy reconstruction), which is a lower bound of the primary neutrino energy and E i, j = min( E, E i j ). The pairs of events i and j inside the hemisphere, H , are counted, with Θ being the Heaviside function. The background expectation enters in the denominator and is obtained by averaging

over a large number of pseudo-experiments using a uniform right ascencion distribution. E min is the minimum energy defining each energy bin. In the data sample of the northern hemisphere, events have neutrino energies between 240 GeV and 1.2 PeV. For the 10% sample, the lower energy bound increases to 6.3 TeV, while for the 1% and 0 1% samples the lower threshold . is 31.6 TeV and 120 TeV, respectively [19]. For the southern hemisphere, the energies range from 490 GeV to 8 PeV. For the 10% sample this range increases to 630 TeV and for the 1% and 0 1% samples the lower threshold lies . at 1.3 PeV and 2.7 PeV. This search is an extension of the Multi Point Source (MPS) analysis that has previously been applied to the Icecube data [26].

## 3.2. Multipole Analyses

A multipole analysis is performed as a second test. It is based on the expansion of the measured skymap into spherical harmonics which are given by

$$\lambda _ { J } ^ { m } ( \vartheta, \varphi ) = \sqrt { \frac { ( 2 \ell + 1 ) ( \ell - m )! } { 4 \pi ( \ell + m )! } } P _ { \ell } ^ { m } ( \cos ( \vartheta ) ) \exp ( i m \varphi ) \,,$$

with glyph[lscript] = 0 1 2 , , , ... , - ≤ glyph[lscript] m ≤ glyph[lscript] and P m glyph[lscript] ( x ) being the Legendre Polynomials. Here, 0 ≤ ϑ &lt; π and 0 ≤ ϕ &lt; 2 π are two free parameters of each spherical harmonic that can be connected to declination( δ ) and right ascension( α ) by ϕ = α and ϑ = π 2 -δ .

Since spherical harmonics form a complete and orthonormal system, any square-integrable function f ( ϑ, ϕ ) on a sphere can be expanded into them, such that it is expressed by a superposition of spherical harmonics and the corresponding complex expansion coefficients a m glyph[lscript] . The expansion coefficients can be obtained by solving the integral

$$a _ { \ell } ^ { m } = \int _ { 0 } ^ { 2 \pi } \mathrm d \varphi \int _ { 0 } ^ { \pi } \mathrm d \vartheta \sin ( \vartheta ) Y _ { \ell } ^ { * m } ( \vartheta, \varphi ) f ( \vartheta, \varphi ) \,.$$

In this analysis, f ( ϑ, ϕ ) is the neutrino arrival direction skymap, described by

$$f ( \vartheta, \varphi ) = \sum _ { i = 0 } ^ { n _ { t o t } } \delta ^ { D } ( \cos ( \vartheta ) - \cos ( \frac { \pi } { 2 } - \delta _ { i } ) ) \cdot \delta ^ { D } ( \varphi - \alpha _ { i } ) \,,$$

where ( δ , i α i ) are the delination and right ascension coordinates of event i and δ D is the Dirac-Delta distribution.

The expansion is performed using the HEALPix software package [27], which splits the sky into 786 432 equal sized bins corresponding to an angular size of ∼ 0 13 (mean bin radius), which is sufficiently smaller than the angular . ◦ resolution of the events.

glyph[negationslash]

For each simulated skymap, the resulting expansion coefficients with m = 0 are complex quantities where the distribution follows a two-dimensional Gaussian centered at the origin of the complex plane. For different simulated signal strengths, µ , different standard deviations σ of the two-dimensional Gaussian are observed, while σ is minimal for pure background. As an example, the complex expansion coefficient for glyph[lscript] = m = 1 is shown in Fig. 4 for 10 000 simulated sky maps of pure background and pure signal of different source strengths µ . Since the distributions are rotationally symmetric about the origin, no separation power exists within the phase of the expansion coefficients φ a ( m glyph[lscript] ), while all separation power rather is contained in their absoulte value | a m glyph[lscript] | .

The coefficients for m = 0 describe the pure zenith-dependence and are real numbers. They contain all power in the background-only case and thus are not centered at the origin of the complex plane. Due to its location at the geographic South Pole and the daily rotation of the earth, possible right ascension systematics due to the detector's geometry are averaged out. Remaining possible deviations between the simulated and true signal event distributions are mainly contained in the zenith, and thus declination, spectrum of the measured signal. Therefore, the m = 0 coefficients carry almost all of this systematic uncertainty.

Since the complex phase of the expansion coefficients ( φ ) is uniformly distributed, and thus carries no separation power between signal and background, the analysis is based on an effective power spectrum ( C eff glyph[lscript] ) defined by glyph[negationslash]

$$C _ { \ell } ^ { \text{eff} } = \frac { 1 } { 2 \ell } \sum _ { \substack { m = - \ell \\ m \neq 0 } } ^ { \ell } | a _ { \ell } ^ { m } | ^ { 2 } \,,$$

using only the absolute values of a m glyph[lscript] . Additionally, the m = 0 coefficients are omitted due to their different behavior and their strong systematic dependencies. Since at large glyph[lscript] of ∼ 100, their contribution to f ( ϑ, ϕ ) becomes very small, the corresponding loss in performance is negligible. The resulting effective power spectra for pure signal and pure background averaged over 10 000 pseudo experiments are shown in Fig. 5 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0634v2-figure-3c3968ea-331c-49f2-825e-ebc73af07d61.jpg)

/Rfractur

·

glyph[negationslash]

Figure 4: Example of the distribution of a complex expansion coefficient, a m =1 glyph[lscript] =1 , shown in the complex plane for 10 000 simulated pure background and pure signal pseudo experiments. For pure signal, two cases of µ = 10 and µ = 30 are shown. All distributions follow a two-dimensional Gaussian centered at the origin, such that all separation power is contained in the absoulte value of a m =1 glyph[lscript] =1 , but none in the phase. This is analogous for all other expansion coefficients with m = 0.

To quantify the differences in the power spectra, a test statistic D 2 eff is defined by

$$D _ { \text{eff} } ^ { 2 } = \frac { 1 } { \sum _ { \ell = 1 } ^ { \text{eff} } w _ { \ell } ^ { \text{eff} } } \sum _ { \ell = 1 } ^ { \ell _ { \max } } \left ( w _ { \ell } ^ { \text{eff} } \cdot \text{sign} _ { \ell } \frac { \left ( C _ { \ell, \exp } ^ { \text{eff} } - \langle C _ { \ell, \text{bg} } ^ { \text{eff} } \rangle \right ) ^ { 2 } } { \sigma _ { C _ { \ell, \text{bg} } ^ { \text{eff} } } ^ { 2 } } \right ) \,, \quad \quad ( 7 ) \\ \dots \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \,.$$

which is motivated by a single-sided, weighted χ 2 -test of the experimentally observed, effective power spectrum C eff glyph[lscript], exp . The spectra 〈 C eff glyph[lscript], sig 〉 and 〈 C eff glyph[lscript], bg 〉 are the mean values of the effective power spectrum for each glyph[lscript] for pure signal and pure background sky maps, respectively. For all glyph[lscript] , they are both averaged over 10 000 pseudo experiments. For the pure background sky maps, the corresponding standard deviation for all glyph[lscript] of the effective power spectrum is called σ C eff glyph[lscript], bg . The parameters w glyph[lscript] and sign glyph[lscript] are defined by

Figure 5: Effective power spectrum C eff glyph[lscript] shown for pure signal sky maps and for various values of µ for an E -2 energy spectrum. As described in section 2.3, for pure signal the number of sources N Sou is chosen such that the map contains as many neutrinos as in the experimental sample. The plot shows the averaged values for 10 000 simulated sky maps.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0634v2-figure-7fa37bc7-5af0-45cf-8493-e9a999b462e1.jpg)

$$w _ { \ell } ^ { \text{eff} } = \frac { \langle C _ { \ell, \text{sig} } ^ { \text{eff} } \rangle - \langle C _ { \ell, \text{bg} } ^ { \text{eff} } \rangle } { \sigma _ { C _ { \ell, \text{bg} } ^ { \text{eff} } } }$$

C

glyph[lscript], bg

$$\text{sign} _ { \ell } = \frac { C _ { \ell, \exp } ^ { \text{eff} } - \langle C _ { \ell, \text{bg} } ^ { \text{eff} } \rangle } { | C _ { \ell, \exp } ^ { \text{eff} } - \langle C _ { \ell, \text{bg} } ^ { \text{eff} } \rangle | } \,,$$

such that each deviation in C eff glyph[lscript] is weighted by the expected deviation in the case of a point source signal. Thus, C eff glyph[lscript] that are very sensitive to point-source signals obtain a large weight, while insensitive C eff glyph[lscript] obtain a small weight to increase sensitivity by keeping the test statistic from being dominated by statistical fluctuations on other angular scales than those relevant for pointsource searches.

Additionally, the parameter sign glyph[lscript] guarantees that only deviations in the expected direction are counted positively, while deviations in the opposite direction are counted negatively. This is a natural definition, because under-

Figure 6: Weight spectrum w glyph[lscript] shown for pure signal sky maps and for various values of µ for an E -2 energy spectrum. As described in section 2.3, for pure signal the number of sources N Sou is chosen such that the map contains as many neutrinos as in the experimental sample. The weights are calculated from averaged values of C eff glyph[lscript] for 10 000 simulated sky maps. The curves exhibit similar shapes and differ only by their total normalization.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0634v2-figure-4cc8cb85-50a4-4482-945b-88d7467523da.jpg)

fluctuations on scales, where point-sources were supposed to cause an excess, are not meant to count positively for the test statistic, since an underfluctuation does rather contradict a point-source signal on this angular scale than support it.

In Fig. 6, the weight spectrum w glyph[lscript] is shown for various source strengths µ and pure signal. For all µ , it exhibits the similar shape, while the shown spectra differ only in their total normalization. Since the test statistic is divided by the sum over all weights (Eq. 7), the test statistic does not depend on the total normalization of the weight spectrum, but only on its shape. Thus, the shown weight spectra lead to the same test statistic.

A softer spectral index, γ , would lead to a slightly poorer angular resolution. This change in angular resolution also alters the characteristic scale of the point sources. Since the weights are large for glyph[lscript] corresponding to the characteristic scale of the structure size, a slightly broader angular resolution would lead to a slight shift of the distribution to lower glyph[lscript] , although this effect

and the corresponding loss in sensitivity are small and thus are neglected in this analysis. Therefore, in the following only one set of weights, which was calculated for µ = 30 and γ = 2, is applied to experimental data.

## 4. Analysis Performance

## 4.1. Discovery Potential of the Autocorrelation Analysis

In this test, the angular scale θ is varied from 0 ◦ to 5 ◦ with a step size of 0.25 . ◦ In addition, four energy bins are used, that contain: all events, the 10 %, 1 % and 0.1 % of the most energetic events observed in data, in order to have a better background suppression. The 0.1 % sample for the northern hemisphere, for example, contains only the 100 highest-energy events. By using different energy thresholds, the discovery potential for high energy signals is improved, while the sensitivity to sources with soft energy spectra is retained. By varying the step size, the scan itself can determine the best energy and θ binnings that maximize the signal [24].

The autocorrelation analysis is performed on the data described in section 2. The test statistic is acquired from over 20 000 pseudo-experiments of randomized data. The test statistic is then fit with a Gamma distribution for the two highest energy bins, which contain the 1 % and 0.1 % of the events with the highest energy and a Gaussian for the two lowest energy bins, containing all events and the 10 % of events with the highest energy. Figure 7 shows the distribution of the number of observed pairs for pure background sky maps, a clustering scale of 0 25 . ◦ and the energy bin containing 1 % of the data with a Gaussian fit. Additionally, two signal scenarios with a uniform distribution of sources in the northern sky are shown. In order to estimate the discovery potential, the number of signal events is varied for different numbers of sources and is compared to the test statistic distribution for background. Furthermore, the thresholds for the energy bins are re-calculated in order to keep the number of events fixed in every energy bin. Given the test statistic distribution for different numbers of signal neutrinos ( n sig ), the required number of signal neutrinos, where 50 % of the cases give at least a 5 σ deviation from the observed background, is defined to be the discovery potential. Applying the solid-angle-averaged effective area ( A eff ( E ) -see section 2.2) of the detector for the three-years (one-year) sample, the resulting discovery

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0634v2-figure-3aa15df1-595e-4b51-8080-67fcc3bf2088.jpg)

/ ndf 2.35 / 16

2

c

1

p0    .1261

0.0

71

-

0

p1

5.86   0.30

-

2

p2   6.249

0.2

5

-

Figure 7: Example of the number of observed pairs in the autocorrelation analysis for simulated data sets with and without signal. The case for the angular scale θ = 0 25 . ◦ and the energy bin that contains the 1 % highest energy events is shown. The randomized data is fitted with a Gaussian and two signal scenarios with a uniform distribution of E -2 sources in the northern sky are shown. For N Sou = 50, the mean number of neutrinos per source µ was adjusted until the total number of signal events added to the full data sample (and replacing randomized data events at the corresponding declinations) was n sig = 750. Similarly for N Sou = 100, the example shown was constructed by adjusting µ per source until the number of signal events inserted in the full data sample was n sig = 1500

potential in terms of n sig can be converted to fluxes. This is done using

$$E ^ { \gamma } \frac { \text{d} \phi } { \text{d} E } = \frac { n _ { \text{sig} } } { T _ { \text{up} } \overset { \infty } { \underset { 0 } { \int } } \text{d} E \, A _ { \text{eff} } ( E ) E ^ { - \gamma } } \,,$$

where n sig is the number of signal neutrinos fixed by the signal parameters ( N Sou , µ, γ ) and T up is the detector uptime.

The discovery potential includes the correction for trial factors that come from testing different angular scales θ and different energy thresholds E min . The first signal model considered here contains sources with equal strength and a uniform distribution in the northern sky. The second scenario is a spatial distribution according to the Green catalogue [28] of SNRs in the Milky Way, and contains 274 SNRs. Thus, sources are distributed according to randomly chosen positions of SNRs of the catalogue and can vary from

50 to 200 sources. Even if this catalogue gives a biased view of Galactic SNRs, we decided to use it for limit calculations in our search. However, for future applications we will consider improved catalogues of galactic objects. Since the sources are distributed only inside the Galactic Plane, this scenario exhibits a larger clustering between the sources compared to the first one. In order to study the high energy events in more detail we perform an additional test (2-pt HE) for the 100 most energetic events in the IC79 data sample. A smaller energy binning in steps of 10 events is applied up to the most energetic 10 events observed in this sample, while the binning in θ stays the same.

## 4.2. Discovery Potential of the Multipole Analysis

The test statistic for mixed sky maps containing only a small fraction of signal events and mainly background is calculated to estimate the analysis performance. This is done by varying the number of sources for fixed values of the mean number of neutrinos per source and the spectral index. The resulting test statistic for each combination of ( N Sou , µ, γ ) is then calculated 1000 times and 10 000 times for the background distribution. The resulting test statistics for µ = 5 are shown in Fig. 8 for an E -2 energy spectrum and various values of N Sou . To facilitate comparison with the autocorrelation analysis, shown in Fig. 7, the final (average) number of signal events is also given. Note that in contrast to Fig. 7, there is no restriction on the event's energy.

The corresponding 5 σ discovery potential in terms of the number of signal neutrinos n sig is defined analogous to section 4.1 and is converted to a physical flux using Eq. 10. Further, in contrast to the autocorrelation analysis, the discovery potential of the multipole analysis is determined only for the case of uniformly distributed sources on the northern hemisphere. Furthermore, it does not use a binning in energy or angular distance, such that there is no correction for trials and the resulting pre-trial significance is also the post-trial value.

## 4.3. Comparison of Autocorrelation and Multipole Analyses

The 5 σ discovery significance for both analyses are shown in Fig. 9(a) and Fig. 9(b) for the northern and southern hemisphere, respectively, assuming an E -2 neutrino spectrum and a uniform source distribution. The significances are compared to the discovery flux of the time-integrated point source likelihood search [15] averaged for each hemisphere and to the recently found

Figure 8: Test statistic for the multipole analysis shown for various amounts of signal in the expanded sky map, given by µ = 5, γ = 2 and various values for N Sou . The average number of signal events n sig , which are added to the data in each case, is also shown. A Gaussian fit is shown for all cases and provides good agreement with each distribution.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0634v2-figure-ec60f8d5-be42-4d05-a13b-0d4cfa7c411e.jpg)

diffuse flux of astrophysical sources [12]. This is achieved by evaluating the best fit astrophysical flux for each spectral index which is converted into a flux per source by normalizing with 2 π for each hemisphere and dividing the diffuse flux by the number of sources, assuming sources of equal flux at Earth. One should note that the point source discovery flux is shown for a single source and does not include trial factors for searching many locations which, over the whole sky, would increase the discovery flux by a factor of about 2.

For an E -2 spectrum, and more than ∼ 20 sources in the northern sky, the autocorrelation analysis is able to identify a signal that the point source likelihood search would not observe, while for the multipole analysis this is the case for more than ∼ 45 sources. The large difference between these two analyses is due to the hard energy spectrum, which is easier to extract with the autocorrelation analysis since it uses energy as an additional observable. For the southern hemisphere, the autocorrelation analysis is performing better than the point source likelihood analysis for more than ∼ 10 sources, while the multipole analysis is not shown, because it was only performed on the northern hemisphere. For the autocorrelation analysis, the 5 σ discovery

flux is shown in Fig. 10(a) and Fig. 10(b) for the northern and southern hemisphere, respectively, assuming an E -2 energy spectrum and a galactic source distribution. As described above, it is also compared to the average discovery flux of the point source likelihood search. In Fig. 11(a), the 5 σ discovery flux for an E -3 neutrino spectrum with uniformly distributed sources is shown and compared to the discovery flux of the point source likelihood search for that spectrum. Above ∼ 20 sources, the multipole and the autocorrelation analysis perform better than the point source likelihood analysis. For this energy spectrum, the discovery potential of both analyses is similar, since the energy observable that is used by the autocorrelation analysis carries only little separation power between astrophysical and atmospheric neutrinos. Figure 11(b) illustrates the discovery flux for sources with an E -2 25 . neutrino spectrum for uniformly distributed sources. Again, the performance of both analyses is similar, while a discovery potential for the point source likelihood analysis is not available.

For a fit of an E -2 spectrum to the HESE data, neither analysis is able to detect the underlying source population except in the case that it consists of very few sources (the scenario for a detection in the previous point source likelihood searches). However, for fits of softer spectral indices to the HESE data like E -3 or E -2 25 . , there are a range of source populations compatible with the HESE flux that could be detected by the autocorrelation and multipole analysis, but not by the point source likelihood analysis. For all tested signal hypotheses, the declination-dependent detector acceptance was correctly taken into account. However, it should be noted that in case of a non-isotropic signal, the given discovery potential is not valid since different declinations contribute differently to the total significances (s. Fig. 3). Additionally, large-scale structures in the source distribution would increase the clustering of events, allowing both methods to detect even smaller fluxes. Thus, the sensitivity depends strongly on the source distribution and can not simply be applied to any model predictions.

The calculations presented in this paper are not limited by computation time, although this does scale differently with the number of events in each sky map. The information-based complexity, and thus the required computation time, of the multipole analysis scales with O ( n tot ), while the autocorrelation analysis scales with O ( n 2 tot ).

## 5. Results

## 5.1. Results of the Autocorrelation Test

The autocorrelation test was applied to the presented data sample of the IC40, IC59 and IC79 configurations. In Fig. 12, the observed number of pairs for the northern hemisphere are shown as a function of the clustering scale θ and compared to the background expectation. The same plots are shown for the southern hemisphere in Fig. 13. Since the data sample contains more events in the southern part of the sky, the fluctuations for the highest energy bin are smaller.

The ratios of the observed number of pairs and the background expectation of this analysis are shown in Fig. 14 for the northern hemisphere and Fig. 15 for the southern hemisphere. In both hemispheres, fewer pairs than expected were observed and a small underfluctuation is visible. The background distribution is used for the evaluation of the local p-values. The best pre-trial p-values are 0.16 for the northern and 0.055 for the southern hemisphere. Taking the trials for the different angular and energy bins into account results in a post-trial p-value of 0.84 for the northern hemisphere and 0.73 for the southern hemisphere. The result of the additional high energy test of IC79 data only (2-pt HE) gives a pre-trial p-value of 0.035. After taking the trial correction into account results into a post-trial p-value of 0.38, with one pair at 1.75 ◦ inside the 10 events energy bin. The pair of events has reconstructed (right ascension, declination) of (285.7 ,+3.1 ) and (287.2 ,+3.6 ). ◦ ◦ ◦ ◦ All results are consistent with fluctuations of the background.

## 5.2. Results of the Multipole Analysis

The method is applied to experimental data from the northern hemisphere by calculating the experimental power spectrum and the value of the test statistic. The observed value of the test statistic is shown in Fig. 16 with the background expectation. This corresponds to a -0 3 . σ deviation from background and thus an underfluctuation with a p-value of 0 63 to find a . higher value of D 2 eff .

## 5.3. Upper Limits on the neutrino flux

Using Eq. 10, the experimental value of the test statistic for both analyses can be converted into upper limits on n sig and thus on the physical flux normalizations (Φ 0 = E γ dΦ d E ) for different spectral indices ( γ ). Upper limits at the 90% confidence level are calculated based on the classical (frequentist)

approach [29]. The resulting limits on the flux normalization for the northern hemisphere and uniformly distributed E -2 neutrino sources are shown in Fig. 9(a). For comparison, additional lines are drawn for the limit of the point source likelihood analysis [15] and for the converted flux of the HESE analysis. The average upper limit of the point source likelihood analysis is 2 7 . · 10 -9 GeV cm / -2 s -1 . Figure 9(b) shows the limit for the southern hemisphere and uniformly distributed E -2 neutrino sources for the autocorrelation analysis. Again, the limit for the point source likelihood analysis and the converted HESE flux are shown for comparison and the average upper limits of the point source likelihood analysis is 5 5 . · 10 -9 GeV cm / -2 s -1 .

In Fig. 10(a) and Fig. 10(b), the limits are shown for the northern and southern hemisphere for an E -2 neutrino spectrum in the galactic plane scenario. Both are compared to the converted HESE flux and the limits from the point source likelihood analysis.

For both analyses, the choice of sampling points is different due to different approaches in determining the limit by varying the signal parameters ( N Sou , µ, γ ). While the autocorrelation analysis fixes N Sou and varies n sig by changing µ to determine the limit, the multipole analysis varies N Sou keeping µ fixed. This results in a different set of sampling points for the limit lines, which has no influence to the lines themselves or their physical interpretation. In Fig. 11(a) and Fig. 11(b), the limits of both analyses are shown for an E -3 and an E -2 25 . spectrum, respectively. They are compared to the converted HESE flux and to the limit of the point source likelihood analysis for E -3 , while for E -2 25 . no limit from the point source likelihood analysis is available. The E -2 25 . spectrum is motivated by the HESE best fit of the spectral index.

## 6. Systematic Influences

The analyses are affected by only few systematic uncertainties due to uncertainties in the background estimation and in the signal efficiency. To estimate the systematic influences we use only the IC79 data sample since it contains most of the events in the combined sample and this is a good approximation. The background estimation is based on randomized experimental data, and systematic effects in the background can only result from preexisting largescale anisotropies in the experimental sample. Therefore, no systematic effects on background estimation are introduced due to MC simulations, e.g. from assumptions on any hadronic models or on the composition of cosmic rays.

Since point source searches look for small-scale anisotropies, these must be distinguishable from possible large-scale structures in the data sample. To investigate the effect of a preexisting large-scale anisotropy in the background, the analyses are repeated including a large-scale anisotropy in the background of atmospheric neutrinos and muons for the mixed sky maps, while keeping the test statistic of the null hypothesis fixed. The large-scale anisotropy is simulated according to a measurement by Milagro [30]. Since the Milagro anisotropy is an anisotropy in cosmic rays, a possible anisotropy of atmospheric muons and neutrinos of about the same order of magnitude is expected. The corresponding systematic errors are given by the resulting shift in the sensitivity of the analyses.

Systematic uncertainties in the signal efficiency arise mainly from the DOMefficiency and optical ice properties. In this context the DOM-efficiency describes the absolute light detection efficiency of the optical modules. The variation of the ice parameters refers to the optical properties of the ice, including absorption and scattering. These systematic effects are estimated using MC simulations analogously to [15]. To estimate the systematic errors, the following three uncertainties are propagated through the analysis while calculating the sensitivity:

- 1. Variation of the DOM-efficiency by ± 10%
- 2. Variation of the Ice Parameters by ± 10%
- 3. Influence of a large-scale anisotropy in the background estimation.

The resulting effects on the sensitivities are shown in Table 1. While the uncertainties for the DOM-efficiency and absorption and scattering in the ice are the same for both analyses, different values for the uncertainties due to a large-scale anisotropy are shown.

## 7. Conclusions

Two methods to search for a small-scale anisotropy with IceCube were presented. The results of both searches are consistent with background expectations with small underfluctuation. Depending on the number of assumed sources, the resulting upper limits range from 10 -8 GeV/cm s 2 -1 for one source to 10 -9 GeV/cm s 2 -1 for 3500 E -2 neutrino sources in the northern hemisphere. Limits were also set for other assumed energy spectra, including E -3 and E -2 25 . in the northern hemisphere. Since both analyses use a data-driven background estimation they are more robust against systematic

Table 1: Systematic errors on the flux normalization for the combined data sample while varying the DOM-efficiency and absorption and scattering of the ice. The uncertainty due to a large-scale anisotropy in the background estimation is calculated separately for each analysis. The maximum uncertainty for both analysis is listed above.

|            | effect on sensitivity   | effect on sensitivity   | effect on sensitivity   | effect on sensitivity   | effect on sensitivity   | effect on sensitivity   |
|------------|-------------------------|-------------------------|-------------------------|-------------------------|-------------------------|-------------------------|
| spectrum   | DOM Eff.                | Abs. & Scatt.           | Large-Scale Aniso.      | Large-Scale Aniso.      | Combined                | Combined                |
|            |                         |                         | Multip.                 | 2-pt.                   | Multip.                 | 2-pt.                   |
| E - 2      | ± 13%                   | ± 6%                    | ± 6%                    | ± 4%                    | ± 16%                   | ± 15%                   |
| E - 2 . 25 | ± 14%                   | ± 6%                    | ± 8%                    | ± 8%                    | ± 17%                   | ± 17%                   |
| E - 3      | ± 26%                   | ± 11%                   | ± 3%                    | ± 6%                    | ± 28%                   | ± 29%                   |

uncertainties than estimations from MC simulations.

Considering the astrophysical flux previously observed in IceCube [12], a small number ( ≤ 10) of isotropically distributed sources in the northern hemisphere of very hard energy spectra, like E -2 , is excluded as it was by former IceCube analyses [15]. For softer energy spectra, the analyses presented here disfavor the observed flux to come from less than ∼ 20 sources for E -2 25 . and from less than ∼ 5000 sources for E -3 . Additionally, for sources distributed along the galactic plane in the northern hemisphere the autocorrelation limit is close to the flux predicted by HESE. In the southern hemisphere, the data sample contains predominantly atmospheric muons from cosmic ray air showers above the detector. Due to this background the autocorrelation analysis is not sensitive to a population of sources at the HESE flux level. For all these tests, the sources are assumed to have the same flux at Earth, since the true spatial flux distribution is not known for the observed astrophysical flux.

For hard energy spectra, the 2-pt correlation analysis is more sensitive than the multipole analysis since it uses the energy information as an additional variable. For soft energy spectra, the multipole analysis becomes slightly more competitive.

Acknowledgements We acknowledge the support from the following agencies: U.S. National Science Foundation-Office of Polar Programs, U.S. National Science Foundation-Physics Division, University of Wisconsin Alumni Research Foundation, the Grid Laboratory Of Wisconsin (GLOW) grid infrastructure at the University of Wisconsin - Madison, the Open Science Grid (OSG) grid infrastructure;

U.S. Department of Energy, and National Energy Research Scientific Computing Center, the Louisiana Optical Network Initiative (LONI) grid computing resources; Natural Sciences and Engineering Research Council of Canada, WestGrid and Compute/Calcul Canada; Swedish Research Council, Swedish Polar Research Secretariat, Swedish National Infrastructure for Computing (SNIC), and Knut and Alice Wallenberg Foundation, Sweden; German Ministry for Education and Research (BMBF), Deutsche Forschungsgemeinschaft (DFG), Helmholtz Alliance for Astroparticle Physics (HAP), Research Department of Plasmas with Complex Interactions (Bochum), Germany; Fund for Scientific Research (FNRS-FWO), FWO Odysseus programme, Flanders Institute to encourage scientific and technological research in industry (IWT), Belgian Federal Science Policy Office (Belspo); University of Oxford, United Kingdom; Marsden Fund, New Zealand; Australian Research Council; Japan Society for Promotion of Science (JSPS); the Swiss National Science Foundation (SNSF), Switzerland; National Research Foundation of Korea (NRF); Danish National Research Foundation, Denmark (DNRF)

## References

- [1] M. Aartsen, et al. (IceCube collaboration), Evidence for high-energy extraterrestrial neutrinos at the IceCube detector, Science 342 1242856 (2013).
- [2] P. M´ esz´ros, a Gamma-ray bursts, Rep.Prog.Phys. 69 (2006) 2259.
- [3] E. Waxman, Neutrino astronomy and gamma-ray bursts, Phil.Trans.R.Soc. A 365 (2007) 1323-1334.
- [4] F. W. Stecker, et al. , High-energy neutrinos from active galactic nuclei, Phys.Rev.Lett. 66 (1991) 2697-2700.
- [5] E. Waxman, J. Bahcall, High energy neutrinos from astrophysical sources: An upper bound, Phys.Rev. D 59 (1998) 023002.
- [6] F. Halzen, E. Zas, Neutrino fluxes from active galaxies: A modelindependent estimate, Astrophys.J. 488 (1997) 669.
- [7] M. Doert, et al. , Correlated neutrino and gamma-ray emission from active galactic nuclei - an estimation, J.P.C.S. 355 (2012) 012039.

| [8]   | W. Atwood, et al. (Fermi-LAT collaboration), The Large Area Telescope on the Fermi Gamma-Ray Space Telescope Mission, Astrophys.J. 697 (2009) 1071.                                          |
|-------|----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [9]   | J. K. Becker, High-energy neutrinos in the context of multimessenger astrophysics, Phys.Rep. 458 (2008) 173-246.                                                                             |
| [10]  | E. Waxman, J. Bahcall, High energy neutrinos from cosmological gamma-ray burst fireballs, Phys.Rev.Lett. 78 (1997) 2292-2295.                                                                |
| [11]  | A. Achterberg, et al. (IceCube collaboration), First year performance of the IceCube Neutrino Telescope, Astropart.Phys. 26 (2006) 155-173.                                                  |
| [12]  | M. Aartsen, et al. (IceCube collaboration), Observation of high- energy astrophysical neutrinos in three years of IceCube data, arXiv eprint:1405.5303 (2014).                               |
| [13]  | R. Abbasi, et al. (IceCube collaboration), The IceCube data ac- quisition system: Signal capture, digitization, and timestamping, Nucl.Instr.Meth.Phys.Res. 601 (2009) 294 - 316.            |
| [14]  | R. Abbasi, et al. (IceCube collaboration), Calibration and characteri- zation of the IceCube photomultiplier tube, Nucl.Instr.Meth.Phys.Res. 618 (2010) 139 - 152.                           |
| [15]  | M. G. Aartsen, et al. (IceCube collaboration), Search for time- independent neutrino emission from astrophysical sources with 3 yr of IceCube data, Astrophys.J. 779 (2013) 132.             |
| [16]  | M. Aartsen, et al. (IceCube collaboration), The IceCube Neutrino Ob- servatory part I: Point source searches, arXiv eprint:1309.6979 (2013).                                                 |
| [17]  | M. Aartsen, et al. (IceCube collaboration), Search for a diffuse flux of astrophysical muon neutrinos with the icecube 59-string configuration, Phys.Rev.D 89 (2014) 062007.                 |
| [18]  | S. Adri´n-Mart´nez, a ı et al. (ANTARES collaboration), Searches for clus- tering in the time integrated skymap of the ANTARES neutrino tele- scope, J. Cosmol. Astropart. Phys. 5 (2014) 1. |
| [19]  | M. Aartsen, et al. (IceCube collaboration), Energy reconstruction meth- ods in the IceCube neutrino telescope, J.Instrum. 9 (2014) P03009.                                                   |

| [20]   | R. Abbasi, et al. (IceCube collaboration), Time-integrated searches for point-like sources of neutrinos with the 40-string IceCube detector, As- trophys.J. 732 (2011) 18.   |
|--------|------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [21]   | J. Ahrens, et al. (AMANDA Collaboration), Muon track reconstruction and data selection techniques in AMANDA, Nucl.Instrum.Meth. A524 (2004) 169-194.                         |
| [22]   | N. Metropolis, S. Ulam, The monte carlo method, J.A.S.A. 44 (1949) 335-341.                                                                                                  |
| [23]   | J. Hammersley, D. Handscomb, Monte Carlo Methods, Methuen's mono- graphs on applied probability and statistics, Methuen, 1964.                                               |
| [24]   | C. B. Finley, S. Westerhoff, On the evidence for clustering in the arrival directions of AGASA's ultrahigh energy cosmic rays, Astropart.Phys. 21 (2004) 359-367.            |
| [25]   | R. Abbasi, et al. (HiRes collaboration), A search for arrival direc- tion clustering in the HiRes-I monocular data above 10**19.5 eV, As- tropart.Phys. 22 (2004) 139-149.   |
| [26]   | Y. Sestayo, E. Resconi, Uncovering neutrinos from cosmic ray factories: the multi point source method, Astropart.Phys. 44 (2013) 15-23.                                      |
| [27]   | K. Gorski, et al. , HEALPix - a framework for high resolution discretiza- tion, and fast analysis of data distributed on the sphere, Astrophys.J. 622 (2005) 759-771.        |
| [28]   | D. A. Green, A catalogue of 294 Galactic supernova remnants, ArXiv e-prints (2014).                                                                                          |
| [29]   | J. Neyman, Outline of a theory of statistical estimation based on the classical theory of probability, Phil.Trans.R.Soc. A 236 (1937) 333.                                   |
| [30]   | A. Abdo, et al. (Milagro collaboration), The large scale cosmic-ray anisotropy as observed with Milagro, Astrophys.J. 698 (2009) 2121- 2130.                                 |

(a) Discovery potential and limits for the northern sky

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0634v2-figure-d31f3785-597f-422a-8ae7-d680616c9739.jpg)

(b) Discovery potential and limits for the southern sky

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0634v2-figure-4cf294ab-35f9-4fbb-a9ce-95fa3b371397.jpg)

Figure 9: Discovery potential and upper limits for uniform E -2 neutrino sources for the autocorrelation analysis and the multipole analysis (a) on the northern hemisphere and (b) on the southern hemisphere. They are compared to the discovery potential of the point source search [15]. The yellow band corresponds to the converted flux of the HESE analysis [12].

(a) Discovery potential and limits for the northern sky

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0634v2-figure-9f0daed5-1dbd-4b14-930d-460790fc1c33.jpg)

(b) Discovery potentials and limits for the southern sky

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0634v2-figure-b0d1a3a6-667d-4409-8e5e-5e972842a111.jpg)

Figure 10: (a) The discovery potential and the upper limits for E -2 neutrino sources, distributed in the galactic plane with the autocorrelation analysis (a) for the northern hemisphere and (b) for the southern hemisphere. They are compared to the discovery potential of the point source search [15].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0634v2-figure-7647448f-b543-4e58-a062-0145fa1c57bb.jpg)

(a) Discovery potential and limits for E -3

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0634v2-figure-8f344d68-7b9c-47e4-a4c1-df395e1d9325.jpg)

- (b) Discovery potential and limits for E -2 25 .

Figure 11: (a) The discovery potential and upper limits (a) for E -3 neutrino sources and (b) for E -2 25 . neutrino sources with the autocorrelation analysis and the multipole analysis for the northern hemisphere. They are compared to the discovery potential and the upper limit of the point source search [15]. Additionally, the converted flux from the HESE analysis is shown [12].

Figure 12: Result for the autocorrelation test on the northern hemisphere as a function of the clustering scale θ . The black points refer to the observed number of pairs, while the black line represents the average number of background pairs.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0634v2-figure-cc5d2233-a5aa-4fd2-b214-100888580a8e.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0634v2-figure-6bd89a15-9741-4924-9074-82f3e9441365.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0634v2-figure-11626750-b75a-46f6-bbe6-67741f66f2e8.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0634v2-figure-5af0fb8e-7282-457c-a1c2-07b952d4ca5c.jpg)

## Southern Hemisphere

Figure 13: Result for the autocorrelation test for the south as a function of the clustering scale θ . The black points refer to the observed number of pairs, while the black line represents the average number of background pairs.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0634v2-figure-d25b07be-3f7b-4641-8204-68e0ed64287f.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0634v2-figure-150e2260-e188-4887-a266-867f3fca3b29.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0634v2-figure-b42832f5-0d4a-4917-8da8-b5631a4318f2.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0634v2-figure-528a5145-52d3-42e3-8942-fda3f2807932.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0634v2-figure-a8624d21-6b87-47fd-8699-dba38a394b74.jpg)

#obs. pairs/#avg. pairs

1.2

1% energy bin

1.1

1

0.9

0.8

0.7

0

1

0.8

0.6

0.4

0.2

0

prob.

#obs. pairs/#avg. pairs prob.

#obs. pairs/#avg. pairs prob.

1.03

1.02

1.01

1

0.99

0.98

0.97

1

0.8

0.6

0.4

0.2

0

0

2.5

2

1.5

1

0.5

0

0

1

0.8

0.6

0.4

0.2

0

0

maximum separation

]

GLYPH&lt;176&gt;

[

q

0.5

1

1.5

2

2.5

3

3.5

4

4.5

5

]

GLYPH&lt;176&gt;

[

q

maximum separation

0

0.5

1

1.5

2

2.5

3

3.5

4

4.5

5

Figure 14: Results for the anisotropy analysis on the northern hemisphere. In each plot, the upper panel shows the ratio of the number of observed pairs and the average number of background pairs, with the +1 σ and -1 σ (light blue), as well as the +2 σ and -2 σ (dark blue) contours as a function of the clustering scale θ . The lower panel shows the probability before trials. The best p-value for the northern hemisphere is 0.16 and is found at θ &lt; 4 ◦ in the highest-energy 10 % selection. The final p-value after correcting for trials is 0.84.

5

theta

0.5

1

1.5

2

2.5

3

3.5

4

4.5

0

0.5

0.5

0.5

1

1

1

10% energy bin final p-value: 0.84

1.5

2

2.5

3

3.5

1.5

maximum separation

2

2.5

3

3.5

0.1% energy bin

1.5

2

2.5

3

]

[

q

3.5

GLYPH&lt;176&gt;

4

4

4

4.5

4.5

4.5

5

theta

5

5

theta

histo7

Entries   20

Mean   2.495

RMS    1.443

Entries   20

Mean     2.64

RMS    1.474

Entries   20

Mean   2.045

RMS     1.531

Entries  100

Mean   2.147

RMS    1.576

Figure 15: Results for the anisotropy analysis for the southern hemisphere. For each plot, the upper panel shows the ratio of the number of observed pairs and the average number of background pairs, with the +1 σ and -1 σ (light blue), as well as the +2 σ and -2 σ (dark blue) contours for the clustering scale θ . The lower panel shows the probability before trials. The best p-value for the southern hemisphere is 0.055 and is found at θ &lt; 4 75 . ◦ in the highest-energy 1 % selection. The final p-value after correcting for trials is 0.73.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0634v2-figure-46cb0385-483c-4f58-bf1d-29aca3afeb98.jpg)

Figure 16: Test statistic for the multipole analysis for pure background compared to the experimentally observed value.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0634v2-figure-e518ff6b-7faf-45ff-b575-a568ee75a8ed.jpg)
IceCube Collaboration M. G. Aartsen M. Ackermann J. Adams J. A. Aguilar M. Ahlers M. Ahrens D. Altmann T. Anderson C. Arguelles T. C. Arlen J. Auffenberg X. Bai S. W. Barwick V. Baum J. J. Beatty J. Becker Tjus K. -H. Becker S. BenZvi P. Berghaus D. Berley E. Bernardini A. Bernhard D. Z. Besson G. Binder D. Bindig M. Bissok E. Blaufuss J. Blumenthal D. J. Boersma C. Bohm F. Bos D. Bose S. Böser O. Botner L. Brayeur H. -P. Bretz A. M. Brown J. Casey M. Casier E. Cheung D. Chirkin A. Christov B. Christy K. Clark L. Classen F. Clevermann S. Coenders D. F. Cowen A. H. Cruz Silva M. Danninger J. Daughhetee J. C. Davis M. Day J. P. A. M. de André C. De Clercq S. De Ridder P. Desiati K. D. de Vries M. de With T. DeYoung J. C. Díaz-Vélez M. Dunkman R. Eagan B. Eberhardt B. Eichmann J. Eisch S. Euler P. A. Evenson O. Fadiran A. R. Fazely A. Fedynitch J. Feintzeig J. Felde T. Feusels K. Filimonov C. Finley T. Fischer-Wasels S. Flis A. Franckowiak K. Frantzen T. Fuchs T. K. Gaisser R. Gaior J. Gallagher L. Gerhardt D. Gier L. Gladstone T. Glüsenkamp A. Goldschmidt G. Golup J. G. Gonzalez J. A. Goodman D. Góra D. Grant P. Gretskov J. C. Groh A. Groß C. Ha C. Haack A. Haj Ismail P. Hallen A. Hallgren F. Halzen K. Hanson D. Hebecker D. Heereman D. Heinen K. Helbing R. Hellauer D. Hellwig S. Hickford G. C. Hill K. D. Hoffman R. Hoffmann A. Homeier K. Hoshina F. Huang W. Huelsnitz P. O. Hulth K. Hultqvist S. Hussain A. Ishihara E. Jacobi J. Jacobsen K. Jagielski G. S. Japaridze K. Jero O. Jlelati M. Jurkovic B. Kaminsky A. Kappes T. Karg A. Karle M. Kauer J. L. Kelley A. Kheirandish J. Kiryluk J. Kläs S. R. Klein J. -H. Köhne G. Kohnen H. Kolanoski A. Koob L. Köpke C. Kopper S. Kopper D. J. Koskinen M. Kowalski A. Kriesten K. Krings G. Kroll M. Kroll J. Kunnen N. Kurahashi T. Kuwabara M. Labare D. T. Larsen M. J. Larson M. Lesiak-Bzdak M. Leuermann J. Leute J. Lünemann J. Madsen G. Maggi R. Maruyama K. Mase H. S. Matis R. Maunu F. McNally K. Meagher M. Medici A. Meli T. Meures S. Miarecki E. Middell E. Middlemas N. Milke J. Miller L. Mohrmann T. Montaruli R. Morse R. Nahnhauer U. Naumann H. Niederhausen S. C. Nowicki D. R. Nygren A. Obertacke S. Odrowski A. Olivas A. Omairat A. O'Murchadha T. Palczewski L. Paul Ö. Penek J. A. Pepper C. Pérez de los Heros C. Pfendner D. Pieloth E. Pinat J. Posselt P. B. Price G. T. Przybylski J. Pütz M. Quinnan L. Rädel M. Rameez K. Rawlins P. Redl I. Rees R. Reimann M. Relich E. Resconi W. Rhode M. Richman B. Riedel S. Robertson J. P. Rodrigues M. Rongen C. Rott T. Ruhe B. Ruzybayev D. Ryckbosch S. M. Saba H. -G. Sander J. Sandroos M. Santander S. Sarkar K. Schatto F. Scheriau T. Schmidt M. Schmitz S. Schoenen S. Schöneberg A. Schönwald A. Schukraft L. Schulte O. Schulz D. Seckel Y. Sestayo S. Seunarine R. Shanidze M. W. E. Smith D. Soldin G. M. Spiczak C. Spiering M. Stamatikos T. Stanev N. A. Stanisha A. Stasik T. Stezelberger R. G. Stokstad A. Stößl E. A. Strahler R. Ström N. L. Strotjohann G. W. Sullivan H. Taavola I. Taboada A. Tamburro A. Tepe S. Ter-Antonyan A. Terliuk G. Tešić S. Tilav P. A. Toale M. N. Tobin D. Tosi M. Tselengidou E. Unger M. Usner S. Vallecorsa N. van Eijndhoven J. Vandenbroucke J. van Santen M. Vehring M. Voge M. Vraeghe C. Walck M. Wallraff Ch. Weaver M. Wellons C. Wendt S. Westerhoff B. J. Whelan N. Whitehorn C. Wichary K. Wiebe C. H. Wiebusch D. R. Williams H. Wissing M. Wolf T. R. Wood K. Woschnagg D. L. Xu X. W. Xu J. P. Yanez G. Yodh S. Yoshida P. Zarzhitsky J. Ziemann S. Zierke M. Zoll
Searches for small-scale anisotropies from neutrino point sources with three years of IceCube data
Recently, IceCube found evidence for a diffuse signal of astrophysical
neutrinos in an energy range of $60\,\mathrm{TeV}$ to the $\mathrm{PeV}$-scale.
The origin of those events, being a key to understanding the origin of cosmic
rays, is still an unsolved question. So far, analyses have not succeeded to
resolve the diffuse signal into point-like sources. Searches including a
maximum-likelihood-ratio test, based on the reconstructed directions and
energies of the detected down- and up-going neutrino candidates, were also
performed on IceCube data leading to the exclusion of bright point sources. In
this paper, we present two methods to search for faint neutrino point sources
in three years of IceCube data, taken between 2008 and 2011. The first method
is an autocorrelation test, applied separately to the northern and southern
sky. The second method is a multipole analysis, which expands the measured data
in the northern hemisphere into spherical harmonics and uses the resulting
expansion coefficients to separate signal from background. With both methods,
the results are consistent with the background expectation with a slightly more
sparse spatial distribution, corresponding to an underfluctuation. Depending on
the assumed number of sources, the resulting upper limit on the flux per source
in the northern hemisphere for an $E^{-2}$ energy spectrum ranges from $1.5
\cdot 10^{-8}\,\mathrm{GeV}/(\mathrm{cm}^2 \mathrm{s})$, in the case of one
assumed source, to $4 \cdot 10^{-10} \,\mathrm{GeV}/(\mathrm{cm}^2
\mathrm{s})$, in the case of $3500$ assumed sources.