The origins of the extraterrestrial neutrinos observed in IceCube have yet to be determined. In this study we perform a one-point fluctuation analysis of the six-year high-energy starting event (HESE) shower data, with fixed non-Poissonian contributions from atmospheric, Galactic and some extragalactic components, as well as an isotropic (and weakly non-Poissonian) template. In addition to the star-forming galaxies and blazars, our analysis suggests the presence of an additional isotropic component, not associated with any known class of sources, with best-fit intensity of $(2.8\pm0.2)\times 10^{-18}\,(E/100~{\rm TeV})^{-2.7\pm 0.5}$ ${\mathrm{cm}^{-2}\, \mathrm{s}^{-1}\,\mathrm{sr}^{-1}\,\mathrm{GeV}^{-1}}$. For the first time, we also consider high-energy extrapolations of several phenomenological models for the diffuse Galactic emission (tuned to both local cosmic-ray data and diffuse gamma-ray emission in the GeV-TeV domain). We demonstrate the potential of our framework in discriminating between different scenarios, with possible implications on the physics of cosmic ray transport in the TeV-PeV range.

After the discovery of extraterrestrial high-energy neutrinos, the next major goal of neutrino telescopes will be identifying astrophysical objects that produce them. The flux of the brightest source $F_{\rm max}$, however, cannot be probed by studying the diffuse neutrino intensity. We aim at constraining $F_{\rm max}$ by adopting a broken power-law flux distribution, a hypothesis supported by observed properties of any generic astrophysical sources. The first estimate of $F_{\rm max}$ comes from the fact that we can only observe one universe, and hence, the expected number of sources above $F_{\rm max}$ cannot be too small compared with one. For abundant source classes such as starburst galaxies, this one-source constraint yields a value of $F_{\rm max}$ that is an order of magnitude lower than the current upper limits from point-source searches. Then we derive upper limits on $F_{\rm max}$ assuming that the angular power spectrum is consistent with neutrino shot noise yet. We find that the limits obtained with upgoing muon neutrinos in IceCube can already be quite competitive, especially for rare but bright source populations such as blazars. The limits will improve nearly quadratically with exposure, and therefore be even more powerful for the next generation of neutrino telescopes.

We perform the first one-point fluctuation analysis of the high-energy neutrino sky. This method reveals itself to be especially suited to contemporary neutrino data, as it allows to study the properties of the astrophysical components of the high-energy flux detected by the IceCube telescope, even with low statistics and in the absence of point source detection. Besides the veto-passing atmospheric foregrounds, we adopt a simple model of the high-energy neutrino background by assuming two main extra-galactic components: star-forming galaxies and blazars. By leveraging multi-wavelength data from Herschel and Fermi, we predict the spectral and anisotropic probability distributions for their expected neutrino counts in IceCube. We find that star-forming galaxies are likely to remain a diffuse background due to the poor angular resolution of IceCube, and we determine an upper limit on the number of shower events that can reasonably be associated to blazars. We also find that upper limits on the contribution of blazars to the measured flux are unfavourably affected by the skewness of the blazar flux distribution. One-point event clustering and likelihood analyses of the IceCube HESE data suggest that this method has the potential to dramatically improve over more conventional model-based analyses, especially for the next generation of neutrino telescopes.

The one-point function (i.e., the isotropic flux distribution) is a complementary method to (anisotropic) two-point correlations in searches for a gamma-ray dark matter annihilation signature. Using analytical models of structure formation and dark matter halo properties, we compute the gamma-ray flux distribution due to annihilations in extragalactic dark matter halos, as it would be observed by the Fermi Large Area Telescope. Combining the central limit theorem and Monte Carlo sampling, we show that the flux distribution takes the form of a narrow Gaussian of `diffuse' light, with an `unresolved point source' power-law tail as a result of bright halos. We argue that this background due to dark matter constitutes an irreducible and significant background component for point-source annihilation searches with galaxy clusters and dwarf spheroidal galaxies, modifying the predicted signal-to-noise ratio. A study of astrophysical backgrounds to this signal reveals that the shape of the total gamma-ray flux distribution is very sensitive to the contribution of a dark matter component, allowing us to forecast promising one-point upper limits on the annihilation cross section. We show that by using the flux distribution at only one energy bin, one can probe the canonical cross section required for explaining the relic density, for dark matter of masses around tens of GeV.