Angular power spectrum analysis on current and future high-energy neutrino data

# Angular power spectrum analysis on current and future high-energy neutrino data

## Abstract

Astrophysical neutrino events have been measured in the last couple of years, which show an isotropic distribution, and the current discussion is their astrophysical origin. We use both isotropic and anisotropic components of the diffuse neutrino data to constrain the contribution of a broad number of extra-galactic source populations to the observed neutrino sky. We simulate up-going muon neutrino events by applying statistical distributions for the flux of extragalactic sources, and by Monte Carlo method we exploit the simulation for current and future IceCube, IceCube-Gen2 and KM3NeT exposures. We aim at constraining source populations by studying their angular patterns, for which we assess the angular power spectrum. We leave the characteristic number of sources () as a free parameter, which is roughly the number of neutrino sources over which the measured intensity is divided. With existing two-year IceCube data, we can already constrain very rare, bright sources with 100. This can be improved to with IceCube-Gen2 and KM3NeT with ten-year exposure, constraining the contribution of BL Lacs (). On the other hand, we can constrain weak sources with large number densities, like starburst galaxies (), if we measure an anisotropic neutrino sky with future observations.

a]Ariane Dekker a,b]and Shin’ichiro Ando

\affiliation

[a]GRAPPA Institute, University of Amsterdam, 1098 XH Amsterdam, Netherlands \affiliation[b]Kavli Institute for the Physics and Mathematics of the Universe, University of Tokyo, Kashiwa, Chiba 277-8583, Japan

## 1 Introduction

High-energy neutrinos with energies between tens of TeV and multiple PeV have been detected with IceCube at South Pole [1, 2]. Distribution of these neutrino events are isotropic, suggesting towards dominant contributions from extragalactic sources [3, 4]. Moreover, the observed neutrino flux is consistent with the extragalactic gamma-ray flux and the ultrahigh-energy cosmic ray flux, which also suggests that neutrinos and high-energy cosmic-rays are produced by the same sources [5, 6]. Theoretical models predict the production of neutrinos within or in the surroundings of cosmic ray accelerator through the decay of charged pions, whereby the pions are produced by hadronuclear collision [7, 8, 9, 10] or photo-hadronic interactions [11, 12, 13] with cosmic rays. Gamma-rays are produced alongside through the decay and production of neutral pions. The two channels result in different ratios between gamma-ray and neutrino flux and combining these observations will therefore offer insight in the production mechanism, as well as an hint for possible neutrino candidate sources. The Fermi-LAT data showed that the gamma-ray point-source emission is dominated by blazars [14, 15, 16]. Indeed, very recently, temporal and positional associations of IceCube events with a blazar, TXS 0506+056, has been found at 3–3.5 level [17, 18]. Even though these could still be caused by atmospheric neutrino background, if confirmed, the blazars will be established as a source of high-energy neutrinos through cosmic ray acceleration followed by hadronic interactions.

However, it is possible that more than one source population gives significant contributions to the measured flux of the high-energy neutrinos. In fact, Ref. [19] pointed out that the blazar flares as was identified by Refs. [17, 18] were responsible up to 10% — only a sub-dominant fraction of the total intensity of the IceCube neutrinos [1, 2]. Therefore, alternative sources are being studied for temporal and positional coincidence with neutrino signals, whereby muon track events are favorable probes due to their good angular resolution and large statistics [20]. Nevertheless, besides XS 0506+056 no other point-source was found, which sets constraints on the contribution of source populations, e.g., gamma-ray bursts [21], core-collapse supernovae [22], active galactic nuclei [23], starburst galaxies [24] and galaxy clusters and groups [25]. Moreover, also a galactic contribution is predicted, which is found to be relatively small and we will thus not consider this in our analysis [26, 27].

Another method in order to reveal sources that give dominant contribution to the measured diffuse flux, is by studying angular clustering of registered events [28, 29, 30] and is shown to be more efficient than the conventional point-source searches [20] for individually bright, rare source populations such as blazars [30]. References [31, 29] set, for instance, stringent limits on the number density of neutrino sources, which has the advantage that a wide range of neutrino source populations is discussed without using prior theoretical models. As complementary method, one-point fluctuation analysis [32] maximises information obtained in single pixels, which helps to constrain both extragalactic and Galactic source contributions [33].

In this work, we study the angular clustering of current and future up-going muon neutrino events, observed by IceCube, IceCube-Gen2 and KM3NeT, by performing an angular power spectrum analysis. Instead of using an equal flux for all sources, we apply statistical distributions based on cosmological considerations, to assign a flux to each extragalactic neutrino source. We use observables as the intensity and isotropy that were measured by IceCube [34]. By performing Monte Carlo simulations, we test the sensitivity of this analysis to a broad range of potential source classes. The work extends fully analytic treatment of Ref. [30] by adopting more realistic simulations of the expected sky map of neutrino events of given exposure. By repeating these Monte Carlo simulations, we make error estimates more robust by taking cross correlations between different multipole bins fully into account. This turns out to yield a significant effect on projections, which weakens sensitivity with the angular power spectrum approach compared with the earlier estimates of Ref. [30] based on analytic calculations. We use the characteristic number of sources () to assess the sensitivities, which is roughly the number of neutrino sources over which the measured intensity is divided, and bright sources like BL Lacs have therefore smaller values () than numerous, weak sources like starburst galaxies (). We are still able to already test the case of a characteristic number of sources of 100 with the existing two-year IceCube data, and this can be improved to with IceCube-Gen2 and KM3NeT with ten-year exposure.

The paper is organized as follows. In Sec. 2 we discuss the data sample used in our analysis, consisting of 21 high-energy up-going muon events from 2 years of IceCube exposure. Furthermore, we discuss the full-sky map simulation, where we applied statistical distributions for the astrophysical flux and the flux model from Ref. [34] for the atmospheric flux. In Sec. 3 we discuss the angular power spectrum analysis. In Sec. 4 we show the results using the data sample, where we can already constrain very bright and rare sources, and the results using future neutrino data. Subsequently, we discuss the implications of the result with respect to known source populations in Sec. 5, and we find significant constraints for sources as BL Lacs and starburst galaxies. Finally, we conclude in Sec. 6.

## 2 Half-sky map of neutrino events

In the following, we discuss the data sample and the sky map of up-going high-energy muon neutrinos, detectable by current and future neutrino telescopes IceCube, IceCube-Gen2 and KM3NeT. Muon neutrinos have a good angular resolution due to the long tracks in the detectors [20]. Moreover, up-going events traverse the Earth and thereby shield off background events from cosmic ray muons, allowing the interaction vertex to happen outside the detector and consequently enlarging the effective area [3].

### 2.1 Data sample

We obtain the data sample from two years of IceCube exposure, in which 35000 neutrino events were detected from the Northern hemisphere for GeV [34]. Most of these events are atmospheric neutrinos, produced when cosmic rays interact with the Earth’s atmosphere, and only the events with the highest energies are most probable of extra-terrestrial origin. For this reason, a threshold on the energy proxy of the events is set with  GeV in the case of IceCube, resulting in 21 muon detections in the two year dataset. The high-energy upgoing muon events are downloaded from https://icecube.wisc.edu/science/data/HE_NuMu_diffuse. The energy spectra of the astrophysical component is taken to be the best-fit to the upgoing muon neutrino data [34]:

 dΦνdEν=1.7×10−18 GeV−1 cm−2 s−1 sr−1(Eν100 TeV)−2.2, (1)

found by assuming an equal flavor composition and equal neutrino and antineutrino ratio. Furthermore, we need a description for the background, which are mainly conventional atmospheric neutrinos produced in the decay of pions and kaons [35]. We use the flux model taken from the above URL, showing a much softer spectrum.

### 2.2 Astrophysical flux distribution

We assume that the source flux distribution of neutrino sources is described by a broken power-law:

 dNsdF=N⋆×⎧⎪ ⎪⎨⎪ ⎪⎩(FF⋆)−α,F⋆

where is the total number of neutrino sources, is the flux of each source, the slopes () are fixed to (2.5, 1.5), is the characteristic flux at the break, and is the characteristic number of sources at (see Ref. [30] for further details). The slope at high-flux regime () is generically true for any sources distributed homogeneously in the local volume where expansion of the Universe can be neglected. We fix the lower flux limit at . It is justified since we are interested in the angular clustering of neutrino events, which is dominated only by sources with large fluxes.

The mean of the flux distribution is related to the measured neutrino intensity , Eq. (1), and can be found through

 ∫dFFdNsdF=4πΦν. (3)

We can analytically calculate the left-hand side using Eq. (2) in order to obtain the characteristic flux, , as

 F⋆=4πΦνN⋆η, (4)

where . The source flux distribution is thus shaped by dividing the measured neutrino intensity over , and we leave as a free parameter for the rest of the paper.

The all-sky map of the neutrino flux is obtained by randomly distributing the total number of sources, found by integrating Eq. (2) over all fluxes, and assigning each source a random flux following the source flux distribution for fixed . Additionally, we consider the case of an infinite number of sources, the isotropic case, where we obtain the flux sky map by equally dividing the measured intensity over all pixels on the sky.

### 2.3 Neutrino sky map simulation

The number of neutrino detections per given energy range from single sources with flux is calculated by multiplying the flux by the exposure (the effective area multiplied by the live time) of the neutrino detector and integrating over the neutrino energy:

 Nν=∫∞Eμ,thdEμ∫∞0dEνdFdEνdE(Eν,Eμ,θ)dEμ, (5)

where is the differential exposure to yield muon energy proxy from the incident neutrino with energy , which also depends on the zenith angle () due to the Earth absorption. The exposure data from IceCube is obtained from the previously mentioned URL. In the case of IceCube, we set a threshold on the muon energy proxy of  GeV, while in the case of the KM3NeT exposure we set a threshold on the muon neutrino energy of  TeV. The KM3NeT effective area used in this analysis is at trigger level, and no cuts were applied in order to reduce background. Before analyzing the neutrino event map, we mask all pixels with , where zenith angle represents the angle between the neutrino detectors’ normal vector and the incoming neutrino direction.

Figure 1 shows three examples of neutrino sky maps by using IceCube exposure, whereby the masked region is shown in dark red. The top panel shows the 21 neutrino events from 2 years of IceCube observations, which are isotropically distributed. The middle and bottom panels illustrate the angular clustering behavior by changing , showing and , respectively, where the exposure is enhanced to 200 years for better interpretation. The sky map with has brighter sources, which can be seen by the maximum number of events coming from one pixel being 65. While having around the same total number of neutrino events, the sky map for shows a maximum of 3 events per pixel due to the sources being very faint. The measured distribution of the events can thus be used to constrain the contribution of bright sources to the intensity due to their large fluctuations, and we will assess the parameter in this analysis.

## 3 Method

### 3.1 Auto-correlation angular power spectrum

The anisotropy of the neutrino events can be assessed through the angular power spectrum (APS), which represents the fluctuations as a function of the angular scale. Hereby, the neutrino event sky map is expanded into spherical harmonics as follows:

 N(θ,ϕ)=∑ℓmaℓmYℓm(θ,ϕ), (6)

where is the neutrino events at their latitude () and longitude (), and are the spherical harmonic functions. The APS is described by averaging the expansion coefficients over the sky:

 Cℓ=12ℓ+1ℓ∑m=−ℓ|aℓm|2. (7)

We computed the APS using anafast tool from the software package HEALPix [36].

### 3.2 Analysis

We perform the simulations for 50 logarithmically-spaced values of , in the range between 10 and , and additionally for the isotropic case of . A set of (IceCube, IceCube-Gen2) and (KM3NeT) Monte Carlo simulations are performed for each parameterization of the sky map. We adopt the angular resolution for the neutrino detectors as for IceCube, for IceCube-Gen2, and for KM3NeT. These resolutions allow a maximum multipole moment of , 329, and 1400, respectively. The specifications are summarized in Table 1.

The APS of MC simulations is illustrated in Fig. 2 in the case of with two years of IceCube exposure, showing the and containment bands as well as the APS of the observed sky map (solid line). The figure clearly shows that, in the case of , the probability distribution of cannot be well represented by the Gaussian, being highly skewed towards higher values of . In order to see this more systematically, and more importantly to discuss detectability of sources using the APS, we adopt the following two analysis methods.

distribution. The first method that we approach is by distributing the mean APS for various values of and compare their differences in shape and normalizations. We therefore take the unweighted average of , denoted as , over the multipole range , where the first 50 multipole moments are removed to reduce uncertainties at large angular scales. We discuss for , , , and . This method assumes no prior neutrino sky map information and therefore allows for a direct comparison with the APS of the observed neutrino sky.

-value. With the second approach, we fit the APS of the observed sky map to the APS of the simulated sky maps in order to test if is compatible, and we use all multipole information. We consider the following , obtained for each :

 χ2(Cℓ)=∑ℓℓ′(Cℓ−Cmeanℓ)(Covℓℓ′)−1(Cℓ′−Cmeanℓ′), (8)

where is the simulated APS, is the mean value at each , and is the covariance matrix, where the latter two are obtained from the set of simulations. For each , we calculate the probability density function (PDF) of , , as well as and . We then compare the value of obtained from the APS of the observed sky map, , and quantify the probability of obtaining the same or more extreme values (towards either greater or smaller direction) of for each , noted as the -value, as follows:

In the following, we look for constraints at 95% confidence level (CL), which is equivalent to .

## 4 Results

Following the methodology summarized in the previous section, we analyze the APS of and simulations for each parameterization of using IceCube, IceCube-Gen2 and KM3NeT specifications.

### 4.1 ¯Cℓ distributions

#### IceCube

The top left panel of Fig. 3 shows the result of the distributions for , , and using two years of IceCube exposure. The distributions for smaller deviate from Gaussian showing a tail towards larger . The obtained from the observed sky map is also illustrated, as the vertical blue line, which has the value  sr. Comparing with the exclusion values in Table 2, we can exclude with 95%CL.

#### IceCube-Gen2

Fluctuations caused by bright sources are expected to become dominant if we increase the exposure, and we expect the distributions to start showing different shapes. Indeed, in the case for ten years of IceCube-Gen2 exposure, as shown in the top right panel of Fig. 3, we find that the distributions are significantly separated from each other, allowing to constrain when comparing with future observations. Upper and lower limits on are obtained with 90%, 95% and 99% CL, and are represented in Table 3. If the measured are beyond the ranges shown in the table, one can exclude the model with at a given CL.

#### KM3NeT

The bottom panel of Fig, 3 shows the distributions using ten years of KM3NeT exposure, where we find that the distributions are again significantly distinguishable from each other in order to constrain when comparing with future . In table 4, we present the corresponding exclusion limits on . The values of are much smaller than those for IceCube and IceCube-Gen2 due to the higher detector resolution, and thus finer pixel size.

### 4.2 Correlation between multipoles

In order to estimate the -value from the analysis (see, Eq. 8), we need to evaluate the covariance matrix . The covariance matrix is a matrix where the diagonal part indicates the standard deviation at each multipole, while the off-diagonal elements represent the correlation between different multipole moments. Especially if the events are dominated by small number of sources, , it is likely that the different multipoles are highly correlated with each other, showing non-negligible off-diagonal components (e.g., [37]). We obtain the covariance matrix

 ρℓℓ′≡Covℓℓ′√CovℓℓCovℓ′ℓ′, (10)

of the APS for each using the Monte Carlo simulations. We find that the correlation depends on and , as illustrated in Fig. 4, which shows the correlation matrix using two years of IceCube exposure in the cases of (top left), (top right) and (bottom). We notice that the isotropic case does not show any correlation between different multipole moments, while decreasing number of sources show an increase in correlation. This implies that, especially at small angular scales, we see a deviation between different , and thus shows the importance of a small angular resolution using this analysis. It also shows the importance of treating the full covariance matrix rather than the diagonal components, as doing the latter will introduce systematic bias of the obtained constraints especially for small values of , which highlight an important difference from the earlier work [30].

### 4.3 p-value analysis

#### IceCube

We tested the parameterizations by fitting the APS of the simulated sky map to the APS of the observed sky map with 21 detections from two years of IceCube observations [34], and the result is presented in Fig. 5 as a red solid curve. Moreover, the blue and light-blue shaded regions show the expected and containment bands, respectively, from the analysis using the APS of the simulated isotropic neutrino sky as mock data, and the green horizontal line shows the exclusion limit with . Following this exclusion line, we find a lower limit of at 95% CL, which is well within the exclusion region using the simulated isotropic sky. This exclusion limit corresponds to very rare and bright sources, which are indeed not expected to dominate since the distribution of the 21 events are consistent with an isotropic expectation.

If we still measure an isotropic sky with 10 years of IceCube exposure, we can exclude with 95% CL, where the region is the exclusion region.

#### IceCube-Gen2

IceCube-Gen2 is expected to have an effective area of a factor of ten times larger than IceCube, and we assessed its sensitivity using ten years of exposure. We first assumed that we keep on measuring an isotropic sky in the future, with mock data , and find the -values presented in the top left panel of Fig. 6, where the shaded regions are the and containment bands from the analyzed APS, and the horizontal green line represents the exclusion limit of . We find an exclusion of bright sources below , where the range represents the band.

If, however, we do measure clustering of events in the future due to bright sources, the contribution of less bright sources can be constrained and the exclusion tendencies will change. The middle and bottom left panels of Fig. 6 show the results by assuming and , respectively, where the shaded areas and green line represent the aforementioned values. Indeed, in these cases we can significantly exclude the contribution of weak sources with large large , as well as very bright sources. Assuming that , we can exclude source populations with larger than , as well as very rare, bright sources with below with CL, where the range represents the band. Assuming , we can constrain rare, bright source populations even further, showing an exclusion band of , while we also find an upper exclusion limit with -values larger than .

#### KM3NeT

The same approach is applied by using ten years of KM3NeT exposure, from which the results are presented in the right panels of Fig. 6. The top right panel shows the case when assuming an isotropic neutrino sky, and we find that sources below can be excluded at 95% CL, where the range represents the band. The cases by assuming and are shown in the middle and bottom right panels of Fig. 6, respectively. The situation shows a lower exclusion band of and an upper exclusion band of . Whereas assuming results in a lower exclusion band of and an upper exclusion limit above . Comparing the results of IceCube-Gen2 and KM3NeT we find similar exclusion trends, possibly due to the factor 10 larger effective area of IceCube-Gen2 while a factor 10 improved angular resolution of KM3NeT. Therefore, these two future detectors, IceCube-Gen2 and KM3NeT, have similar sensitivities from the APS measurements, confirming theoretical finding of Ref. [30].

#### Case of fractional contributions

Up to now, we assumed that a single population of the sources characterized with fully contributes to the measured neutrino intensity. However, this might not be the case and different source populations can dominate the neutrino sky. We therefore also consider smaller contributions with respect to the isotropic intensity. In particular, we define the total intensity as , where is the intensity due to population and is that from the isotropic component (), and study the exclusion patterns for the fractions and . We perform Monte Carlo simulations using 10 years of IceCube-Gen2 exposure. Figure 7 shows the results with (blue), (green) and (red). We find no significant exclusion in the case of , which is consistent with previous work that predict a blazar contribution of up to . In the case of , we do find a stringent exclusion limit of , where the range represents the band. Since the results of IceCube-Gen2 and KM3NeT show similar exclusion trends (Fig. 6), we also expect similar results for the latter.

## 5 Discussion

BL Lacs, flat-spectrum radio quasars (FSRQs), radio galaxies, and starburst galaxies are generally thought to be the main contributers of the measured intensity [38, 7, 39, 40, 29], and we therefore relate the obtained exclusion values to these source classes. Reference [30] found the () parameterization for BL Lacs, and starburst galaxies of ( and (), respectively, and it is also argued that the typical value for radio galaxies is 10. The parameterizations are obtained from gamma-ray observations by assuming a linear relation between the gamma-ray and neutrino luminosity, , which can be estimated if neutrinos are produced through hadronic interactions. Lepto-hadronic sources show a more complicated luminosity relation, e.g., for BL Lacs [41] or for FSRQs [42]. We thus note that a more realistic parameterization for lepto-hadronic source classes could give stronger constraints. In fact, Ref. [43] mentions that blazars produce high-energy neutrinos for less than through hadronic interactions.

Comparing from the aforementioned source classes with the exclusion limits obtained from analyzing the APS using two years of IceCube exposure, we find that the limits are not significant yet. Increasing the exposures to 10 years for the neutrino detectors IceCube, IceCube-Gen2 and KM3NeT, shows more promising results. If we assume an isotropic sky in the future, shown in the top panels of Fig. 6, we find that BL Lacs are excluded as dominant contributers to the intensity with CL. The non-detection of clustered events does not constrain source classes with large , which produce neutrinos typically in cosmic ray reservoir models like starburst galaxies. However, if, in the future we observe clustered events, we show that our analysis can constrain large populations as well. Indeed, the middle and bottom panels in Fig. 6 show respectively the cases assuming a and neutrino sky. The contribution of starburst galaxies () is significantly constrained in both cases, and radio galaxies () are well constrained in the case of a neutrino sky with .

In order to compare the results with more physical values, we convert to the luminosity using Eq. (23) of Ref. [30], assuming that all sources have equal luminosity. We use the relation between the luminosity and the local number density of neutrino sources from Ref. [31], which is derived from the observed diffuse flux with a neutrino spectrum, and by assuming no redshift evolution. The result is shown in Fig. 8. The blue band represents the region that explains the observed diffuse flux, where the upper limit (blue) represents a contribution of to the diffuse flux, and the lower limit (red) a contribution of . The black dashed lines show the exclusion limits, taken from the results in Fig. 7, whereby the width represents the band. The gray region is excluded when no clustering of events is found in future. BL Lacs lie in the exclusion region, and moreover, lie outside the blue band representing the diffuse flux. Also galaxy clusters (GC) lie above the blue band, and inside the exclusion band, and are thus too bright to explain well the isotropic neutrino emission. FSRQs lie inside the diffuse flux band, however, due to their small number density and large luminosity, a significant clustering of events is expected, and they are therefore testable sources with APS. Low-luminosity AGNs (LL AGN), Fanaroff-Riley galaxies of type II (FR-II) and starburst galaxies (SBG) are consistent with the isotropic neutrino emission. These sources could thus be observed as point-sources in the near future with IceCube-Gen2 and KM3NeT. FR-I will be difficult to observe due to their low luminosity and number density, and they will hardly contribute to the high-energy neutrino emission.

We tested that our results are robust against varying the parameters , and . However, the results will probably change by varying the flux normalization and spectral index, for which we applied ; see Eq. (1). A spectral hardening at high energies has been suggested [44], and since the spectral index depends on the neutrino production model, this could indicate multiple production mechanisms or source populations [44]. If the spectral index would harden with better statistics in the future, we expect to have stronger constraints.

## 6 Conclusion

We have derived constraints on , by analyzing the angular power spectrum for the current and future high-energy up-going muon neutrino data, detectable by IceCube, IceCube-Gen2 and KM3NeT. The constraints were established by applying statistical distributions to the flux of the individual extragalactic sources. Our finding, by using two years of IceCube observations [34] with only 21 high-energy up-going muon events, is already a constraint on the characteristic source number with CL, excluding the dominant contribution from very bright source populations. Comparing this result with known source classes is not effective yet, however, with 10 years of future neutrino data using IceCube, IceCube-Gen2 and KM3NeT exposure, we expect to significantly constrain BL Lacs and FSRQs if the distribution remains isotropic. In particular, 10 years of IceCube-Gen2 observations excludes with CL characteristic source number less than , and KM3NeT less than , where the range represents the band and where BL Lacs are found to have . On the other hand, by observing bright sources in the future we can also find constraints on weak source classes with large number densities, such as starburst galaxies (), which could still be the case with current isotropic measurements. The angular power spectrum analysis on future neutrino data has been found to be a powerful probe to understand what astrophysical sources are dominating the neutrino sky, and in particular to predict what source classes will be observable with future neutrino telescopes, illustrated in Fig. 8 for various sources.

\acknowledgments

We thank Aart Heijboer for helpful discussions, especially for specification of KM3NeT. SA acknowledges support by JSPS KAKENHI Grant Numbers JP17H04836, JP18H04340, and JP18H04578.

### References

1. IceCube Collaboration, M. G. Aartsen et al., “Evidence for High-Energy Extraterrestrial Neutrinos at the IceCube Detector,” Science 342 (2013) 1242856, arXiv:1311.5238 [astro-ph.HE].
2. IceCube Collaboration, M. G. Aartsen et al., “Observation of High-Energy Astrophysical Neutrinos in Three Years of IceCube Data,” Phys. Rev. Lett. 113 (2014) 101101, arXiv:1405.5303 [astro-ph.HE].
3. IceCube Collaboration, M. G. Aartsen et al., “Observation and Characterization of a Cosmic Muon Neutrino Flux from the Northern Hemisphere using six years of IceCube data,” Astrophys. J. 833 no. 1, (2016) 3, arXiv:1607.08006 [astro-ph.HE].
4. IceCube Collaboration, M. G. Aartsen et al., “The IceCube Neutrino Observatory - Contributions to ICRC 2015 Part II: Atmospheric and Astrophysical Diffuse Neutrino Searches of All Flavors,” in Proceedings, 34th International Cosmic Ray Conference (ICRC 2015): The Hague, The Netherlands, July 30-August 6, 2015. 2015.
5. M. Ahlers and F. Halzen, “Opening a New Window onto the Universe with IceCube,” Prog. Part. Nucl. Phys. 102 (2018) 73–88, arXiv:1805.11112 [astro-ph.HE].
6. E. Waxman, “IceCube’s Neutrinos: The beginning of extra-Galactic neutrino astrophysics?,” in Proceedings, 9th Rencontres du Vietnam: Windows on the Universe: Quy Nhon, Vietnam, August 11-17, 2013, pp. 161–168. 2013.
7. A. Loeb and E. Waxman, “The Cumulative background of high energy neutrinos from starburst galaxies,” JCAP 0605 (2006) 003, arXiv:astro-ph/0601695 [astro-ph].
8. K. Murase, S. Inoue, and S. Nagataki, “Cosmic Rays Above the Second Knee from Clusters of Galaxies and Associated High-Energy Neutrino Emission,” Astrophys. J. 689 (2008) L105, arXiv:0805.0104 [astro-ph].
9. K. Fang and K. Murase, “Linking High-Energy Cosmic Particles by Black Hole Jets Embedded in Large-Scale Structures,” Nature Phys. 14 no. 4, (2018) 396–398, arXiv:1704.00015 [astro-ph.HE].
10. K. Murase, M. Ahlers, and B. C. Lacki, “Testing the Hadronuclear Origin of PeV Neutrinos Observed with IceCube,” Phys. Rev. D88 no. 12, (2013) 121301, arXiv:1306.3417 [astro-ph.HE].
11. K. Murase, D. Guetta, and M. Ahlers, “Hidden Cosmic-Ray Accelerators as an Origin of TeV-PeV Cosmic Neutrinos,” Phys. Rev. Lett. 116 no. 7, (2016) 071101, arXiv:1509.00805 [astro-ph.HE].
12. X. Rodrigues, A. Fedynitch, S. Gao, D. Boncioli, and W. Winter, “Neutrinos and Ultra-High-Energy Cosmic-Ray Nuclei from Blazars,” Astrophys. J. 854 no. 1, (2018) 54, arXiv:1711.02091 [astro-ph.HE].
13. C. Guépin and K. Kotera, “Can we observe neutrino flares in coincidence with explosive transients?,” Astron. Astrophys. 603 (2017) A76, arXiv:1701.07038 [astro-ph.HE].
14. M. Lisanti, S. Mishra-Sharma, L. Necib, and B. R. Safdi, “Deciphering Contributions to the Extragalactic Gamma-Ray Background from 2 GeV to 2 TeV,” Astrophys. J. 832 no. 2, (2016) 117, arXiv:1606.04101 [astro-ph.HE].
15. Fermi-LAT Collaboration, F. Acero et al., “Fermi Large Area Telescope Third Source Catalog,” Astrophys. J. Suppl. 218 no. 2, (2015) 23, arXiv:1501.02003 [astro-ph.HE].
16. Fermi-LAT Collaboration, M. Ackermann et al., “2FHL: The Second Catalog of Hard Fermi-LAT Sources,” Astrophys. J. Suppl. 222 no. 1, (2016) 5, arXiv:1508.04449 [astro-ph.HE].
17. IceCube Collaboration, “Neutrino emission from the direction of the blazar TXS 0506+056 prior to the IceCube-170922A alert,” Science 361 no. 6398, (2018) 147–151.
18. Liverpool Telescope, Swift/NuSTAR, MAGIC, H.E.S.S., AGILE, Kiso, VLA/17B-403, INTEGRAL, Kapteyn, Subaru, HAWC, Fermi-LAT, ASAS-SN, VERITAS, Kanata, IceCube Collaboration, “Multimessenger observations of a flaring blazar coincident with high-energy neutrino IceCube-170922A,” Science 361 no. 6398, (2018) eaat1378.
19. K. Murase, F. Oikonomou, and M. Petropoulou, “Blazar Flares as an Origin of High-Energy Cosmic Neutrinos?,” arXiv:1807.04748 [astro-ph.HE].
20. IceCube Collaboration, M. G. Aartsen et al., “All-sky search for time-integrated neutrino emission from astrophysical sources with 7 years of IceCube data,” arXiv:1609.04981 [astro-ph.HE].
21. IceCube Collaboration, M. G. Aartsen et al., “Constraints on minute-scale transient astrophysical neutrino sources,” Submitted to: Phys. Rev. Lett. (2018) , arXiv:1807.11492 [astro-ph.HE].
22. A. Esmaili and K. Murase, “Constraining high-energy neutrinos from choked-jet supernovae with IceCube high-energy starting events,” arXiv:1809.09610 [hep-ph].
23. K. Murase, “Active Galactic Nuclei as High-Energy Neutrino Sources,” pp. 15–31. 2017.
24. K. Bechtol, M. Ahlers, M. Di Mauro, M. Ajello, and J. Vandenbroucke, “Evidence against star-forming galaxies as the dominant source of IceCube neutrinos,” arXiv:1511.00688 [astro-ph.HE].
25. F. Zandanel, I. Tamborra, S. Gabici, and S. Ando, “Neutrinos from Clusters of Galaxies and Radio Constraints,” PoS ICRC2015 (2016) 1117, arXiv:1509.01367 [astro-ph.HE].
26. IceCube Collaboration, M. G. Aartsen et al., “Constraints on Galactic Neutrino Emission with Seven Years of IceCube Data,” arXiv:1707.03416 [astro-ph.HE].
27. F. Halzen, A. Kheirandish, and V. Niro, “Prospects for Detecting Galactic Sources of Cosmic Neutrinos with IceCube: An Update,” Astropart. Phys. 86 (2017) 46–56, arXiv:1609.03072 [astro-ph.HE].
28. IceCube Collaboration, M. G. Aartsen et al., “Searches for small-scale anisotropies from neutrino point sources with three years of IceCube data,” Astropart. Phys. 66 (2015) 39–52, arXiv:1408.0634 [astro-ph.HE].
29. K. Murase and E. Waxman, “Constraining High-Energy Cosmic Neutrino Sources: Implications and Prospects,” Phys. Rev. D94 no. 10, (2016) 103006, arXiv:1607.01601 [astro-ph.HE].
30. S. Ando, M. R. Feyereisen, and M. Fornasa, “How bright can the brightest neutrino source be?,” Phys. Rev. D95 no. 10, (2017) 103003, arXiv:1701.02165 [astro-ph.HE].
31. P. Mertsch, M. Rameez, and I. Tamborra, “Detection prospects for high energy neutrino sources from the anisotropic matter distribution in the local universe,” JCAP 1703 no. 03, (2017) 011, arXiv:1612.07311 [astro-ph.HE].
32. M. R. Feyereisen, I. Tamborra, and S. Ando, “One-point fluctuation analysis of the high-energy neutrino sky,” JCAP 1703 (2017) 057, arXiv:1610.01607 [astro-ph.HE].
33. M. R. Feyereisen, D. Gaggero, and S. Ando, “One-point fluctuation analysis of IceCube neutrino events outlines a significant unassociated isotropic component and constrains the Galactic contribution,” Phys. Rev. D97 no. 10, (2018) 103017, arXiv:1711.02820 [astro-ph.HE].
34. IceCube Collaboration, M. G. Aartsen et al., “Evidence for Astrophysical Muon Neutrinos from the Northern Sky with IceCube,” Phys. Rev. Lett. 115 no. 8, (2015) 081102, arXiv:1507.04005 [astro-ph.HE].
35. T. K. Gaisser, Cosmic rays and particle physics. 1990.
36. K. M. Gorski, E. Hivon, A. J. Banday, B. D. Wandelt, F. K. Hansen, M. Reinecke, and M. Bartelman, “Healpix – a framework for high resolution discretization, and fast analysis of data distributed on the sphere,” astro-ph/0409513. https://arxiv.org/abs/astro-ph/0409513.
37. S. Ando, A. Benoit-Lévy, and E. Komatsu, “Angular power spectrum of galaxies in the 2MASS Redshift Survey,” Mon. Not. Roy. Astron. Soc. 473 no. 4, (2018) 4318–4325, arXiv:1706.05422 [astro-ph.CO].
38. IceCube Collaboration, M. G. Aartsen et al., “The contribution of Fermi-2LAC blazars to the diffuse TeV-PeV neutrino flux,” Astrophys. J. 835 no. 1, (2017) 45, arXiv:1611.03874 [astro-ph.HE].
39. I. Tamborra, S. Ando, and K. Murase, “Star-forming galaxies as the origin of diffuse high-energy backgrounds: Gamma-ray and neutrino connections, and implications for starburst history,” JCAP 1409 (2014) 043, arXiv:1404.1189 [astro-ph.HE].
40. D. Hooper, “A Case for Radio Galaxies as the Sources of IceCube’s Astrophysical Neutrino Flux,” JCAP 1609 no. 09, (2016) 002, arXiv:1605.06504 [astro-ph.HE].
41. F. Tavecchio and G. Ghisellini, “High-energy cosmic neutrinos from spine-sheath BL Lac jets,” Mon. Not. Roy. Astron. Soc. 451 no. 2, (2015) 1502–1510, arXiv:1411.2783 [astro-ph.HE].
42. K. Murase, Y. Inoue, and C. D. Dermer, “Diffuse Neutrino Intensity from the Inner Jets of Active Galactic Nuclei: Impacts of External Photon Fields and the Blazar Sequence,” Phys. Rev. D90 no. 2, (2014) 023007, arXiv:1403.4089 [astro-ph.HE].
43. F. Krauß, K. Deoskar, C. Baxter, M. Kadler, M. Kreter, M. Langejahn, K. Mannheim, P. Polko, B. Wang, and J. Wilms, “Fermi-LAT counterparts of IceCube neutrinos above 100 TeV,” arXiv:1810.08482 [astro-ph.HE].
44. M. Kowalski and I. Collaboration, “Neutrino astronomy with icecube and beyond,” Journal of Physics: Conference Series 888 no. 1, (2017) 012007. http://stacks.iop.org/1742-6596/888/i=1/a=012007.
You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters