Search for transient optical counterparts to high-energy IceCube neutrinos with Pan-STARRS1

In order to identify the sources of the observed diffuse high-energy neutrino flux, it is crucial to discover their electromagnetic counterparts. To increase the sensitivity of detecting counterparts of transient or variable sources by telescopes with a limited field of view, IceCube began releasing alerts for single high-energy (E-v > 60 TeV) neutrino detections with sky localisation regions of order 1 degrees radius in 2016. We used Pan-STARRS1 to follow-up five of these alerts during 2016-2017 to search for any optical transients that may be related to the neutrinos. Typically 10-20 faint m(ip1) less than or similar to 22.5 mag) extragalactic transients are found within the Pan-STARRS1 footprints and are generally consistent with being unrelated field supernovae (SNe) and AGN. We looked for unusual properties of the detected transients, such as temporal coincidence of explosion epoch with the IceCube timestamp, or other peculiar light curve and physical properties. We found only one transient that had properties worthy of a specific follow-up. In the Pan-STARRS1 imaging for IceCube-160427A (probability to be of astrophysical origin of similar to 50%), we found a SN PS16cgx, located at 10.0' from the nominal IceCube direction. Spectroscopic observations of PS16cgx showed that it was an H-poor SN at redshift z = 0.2895 +/- 0.0001. The spectra and light curve resemble some high-energy Type Ic SNe, raising the possibility of a jet driven SN with an explosion epoch temporally coincident with the neutrino detection. However, distinguishing Type Ia and Type Ic SNe at this redshift is notoriously difficult. Based on all available data we conclude that the transient is more likely to be a Type Ia with relatively weak Sin absorption and a fairly normal rest-frame r-band light curve. If, as predicted, there is no high-energy neutrino emission from Type Ia SNe, then PS16cgx must be a random coincidence, and unrelated to the IceCube-160427A. We find no other plausible optical transient for any of the five IceCube events observed down to a 5 sigma limiting magnitude of mip1 approximate to 22 mag, between 1 day and 25 days after detection

Differential limit on the extremely-high-energy cosmic neutrino flux in the presence of astrophysical background from nine years of IceCube data

We report a quasidifferential upper limit on the extremely-high-energy (EHE) neutrino flux above 5×106  GeV based on an analysis of nine years of IceCube data. The astrophysical neutrino flux measured by IceCube extends to PeV energies, and it is a background flux when searching for an independent signal flux at higher energies, such as the cosmogenic neutrino signal. We have developed a new method to place robust limits on the EHE neutrino flux in the presence of an astrophysical background, whose spectrum has yet to be understood with high precision at PeV energies. A distinct event with a deposited energy above 106  GeV was found in the new two-year sample, in addition to the one event previously found in the seven-year EHE neutrino search. These two events represent a neutrino flux that is incompatible with predictions for a cosmogenic neutrino flux and are considered to be an astrophysical background in the current study. The obtained limit is the most stringent to date in the energy range between 5×106 and 2×1010  GeV. This result constrains neutrino models predicting a three-flavor neutrino flux of E2νϕνe+νμ+ντ≃2×10−8  GeV/cm2 sec sr at 109  GeV. A significant part of the parameter space for EHE neutrino production scenarios assuming a proton-dominated composition of ultra-high-energy cosmic rays is disfavored independently of uncertain models of the extragalactic background light which previous IceCube constraints partially relied on

Efficient propagation of systematic uncertainties from calibration to analysis with the SnowStorm method in IceCube

Efficient treatment of systematic uncertainties that depend on a large number of nuisance parameters is a persistent difficulty in particle physics and astrophysics experiments. Where low-level effects are not amenable to simple parameterization or re-weighting, analyses often rely on discrete simulation sets to quantify the effects of nuisance parameters on key analysis observables. Such methods may become computationally untenable for analyses requiring high statistics Monte Carlo with a large number of nuisance degrees of freedom, especially in cases where these degrees of freedom parameterize the shape of a continuous distribution. In this paper we present a method for treating systematic uncertainties in a computationally efficient and comprehensive manner using a single simulation set with multiple and continuously varied nuisance parameters. This method is demonstrated for the case of the depth-dependent effective dust distribution within the IceCube Neutrino Telescope.M.G. Aartsen, M. Ackermann, J. Adams, J.A. Aguilar, M. Ahlers ... Gary C. Hil

A search for neutrino point-source populations in 7 yr of IceCube data with neutrino-count statistics

The presence of a population of point sources in a data set modifies the underlying neutrino-count statistics from the Poisson distribution. This deviation can be exactly quantified using the non-Poissonian template fitting technique, and in this work we present the first application of this approach to the IceCube high-energy neutrino data set. Using this method, we search in 7 yr of IceCube data for point-source populations correlated with the disk of the Milky Way, the Fermi bubbles, the Schlegel, Finkbeiner, and Davis dust map, or with the isotropic extragalactic sky. No evidence for such a population is found in the data using this technique, and in the absence of a signal, we establish constraints on population models with source-count distribution functions that can be described by a power law with a single break. The derived limits can be interpreted in the context of many possible source classes. In order to enhance the flexibility of the results, we publish the full posterior from our analysis, which can be used to establish limits on specific population models that would contribute to the observed IceCube neutrino flux

All-sky measurement of the anisotropy of cosmic rays at 10 TeV and mapping of the local interstellar magnetic field

We present the first full-sky analysis of the cosmic ray arrival direction distribution with data collected by the High-Altitude Water Cherenkov and IceCube observatories in the northern and southern hemispheres at the same median primary particle energy of 10 TeV. The combined sky map and angular power spectrum largely eliminate biases that result from partial sky coverage and present a key to probe into the propagation properties of TeV cosmic rays through our local interstellar medium and the interaction between the interstellar and heliospheric magnetic fields. From the map, we determine the horizontal dipole components of the anisotropy δ 0h = 9.16 × 10−4 and δ 6h = 7.25 × 10−4 (±0.04 × 10−4). In addition, we infer the direction (229fdg2 ± 3fdg5 R.A., 11fdg4 ± 3fdg0 decl.) of the interstellar magnetic field from the boundary between large-scale excess and deficit regions from which we estimate the missing corresponding vertical dipole component of the large-scale anisotropy to be ${\delta }_{N}\sim -{3.97}_{-2.0}^{+1.0}\times {10}^{-4}$