The current and upcoming generation of Very Large Volume Neutrino Telescopes---collecting unprecedented quantities of neutrino events---can be used to explore subtle effects in oscillation physics, such as (but not restricted to) the neutrino mass ordering. The sensitivity of an experiment to these effects can be estimated from Monte Carlo simulations. With the high number of events that will be collected, there is a trade-off between the computational expense of running such simulations and the inherent statistical uncertainty in the determined values. In such a scenario, it becomes impractical to produce and use adequately-sized sets of simulated events with traditional methods, such as Monte Carlo weighting. In this work we present a staged approach to the generation of binned event distributions in order to overcome these challenges. By combining multiple integration and smoothing techniques which address limited statistics from simulation it arrives at reliable analysis results using modest computational resources.

Many Galactic sources of gamma rays, such as supernova remnants, are expected to produce neutrinos with a typical energy cutoff well below 100 TeV. For the IceCube Neutrino Observatory located at the South Pole, the southern sky, containing the inner part of the Galactic plane and the Galactic Center, is a particularly challenging region at these energies, because of the large background of atmospheric muons. In this paper, we present recent advancements in data selection strategies for track-like muon neutrino events with energies below 100 TeV from the southern sky. The strategies utilize the outer detector regions as veto and features of the signal pattern to reduce the background of atmospheric muons to a level which, for the first time, allows IceCube searching for point-like sources of neutrinos in the southern sky at energies between 100 GeV and several TeV in the muon neutrino charged current channel. No significant clustering of neutrinos above background expectation was observed in four years of data recorded with the completed IceCube detector. Upper limits on the neutrino flux for a number of spectral hypotheses are reported for a list of astrophysical objects in the southern hemisphere.

The past decade has welcomed the emergence of cosmic neutrinos as a new messenger to explore the most extreme environments of the universe. The discovery measurement of cosmic neutrinos, announced by IceCube in 2013, has opened a new window of observation that has already resulted in new fundamental information that holds the potential to answer key questions associated with the high-energy universe, including: what are the sources in the PeV sky and how do they drive particle acceleration; where are cosmic rays of extreme energies produced, and on which paths do they propagate through the universe; and are there signatures of new physics at TeV-PeV energies and above? The planned advancements in neutrino telescope arrays in the next decade, in conjunction with continued progress in broad multimessenger astrophysics, promise to elevate the cosmic neutrino field from the discovery to the precision era and to a survey of the sources in the neutrino sky. The planned detector upgrades to the IceCube Neutrino Observatory, culminating in IceCube-Gen2 (an envisaged $400M facility with anticipated operation in the next decade, described in this white paper) are the cornerstone that will drive the evolution of neutrino astrophysics measurements.

After the identification of the gamma-ray blazar TXS 0506+056 as the first compelling IceCube neutrino source candidate, we perform a systematic analysis of all high-energy neutrino events satisfying the IceCube realtime trigger criteria. We find one additional known gamma-ray source, the blazar GB6 J1040+0617, in spatial coincidence with a neutrino in this sample. The chance probability of this coincidence is 30% after trial correction. For the first time, we present a systematic study of the gamma-ray flux, spectral and optical variability, and multi-wavelength behavior of GB6 J1040+0617 and compare it to TXS 0506+056. We find that TXS 0506+056 shows strong flux variability in the Fermi-LAT gamma-ray band, being in an active state around the arrival of IceCube-170922A, but in a low state during the archival IceCube neutrino flare in 2014/15. In both cases the spectral shape is statistically compatible ($\leq 2\sigma$) with the average spectrum showing no indication of a significant relative increase of a high-energy component. While the association of GB6 J1040+0617 with the neutrino is consistent with background expectations, the source appears to be a plausible neutrino source candidate based on its energetics and multi-wavelength features, namely a bright optical flare and modestly increased gamma-ray activity. Finding one or two neutrinos originating from gamma-ray blazars in the given sample of high-energy neutrinos is consistent with previously derived limits of neutrino emission from gamma-ray blazars, indicating the sources of the majority of cosmic high-energy neutrinos remain unknown.

Contributions from the IceCube Collaboration presented at the 36th International Cosmic Ray Conference, 24 July - 1 August 2019, Madison, Wisconsin, USA.

Low background searches for astrophysical neutrino sources anywhere in the sky can be performed using cascade events induced by neutrinos of all flavors interacting in IceCube with energies as low as ~1 TeV. Previously, we showed that even with just two years of data, the resulting sensitivity to sources in the southern sky is competitive with IceCube and ANTARES analyses using muon tracks induced by charge current muon neutrino interactions - especially if the neutrino emission follows a soft energy spectrum or originates from an extended angular region. Here, we extend that work by adding five more years of data, significantly improving the cascade angular resolution, and including tests for point-like or diffuse Galactic emission to which this dataset is particularly well-suited. For many of the signal candidates considered, this analysis is the most sensitive of any experiment. No significant clustering was observed, and thus many of the resulting constraints are the most stringent to date. In this paper we will describe the improvements introduced in this analysis and discuss our results in the context of other recent work in neutrino astronomy.

In order to identify the sources of the observed diffuse high-energy neutrino flux, it is crucial to discover their electromagnetic counterparts. IceCube began releasing alerts for single high-energy ($E > 60$ TeV) neutrino detections with sky localisation regions of order 1 deg 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 < 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. 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 $\sim$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 z = 0.2895. 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 SiII 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 $m \sim 22$ mag, between 1 day and 25 days after detection.

High-energy neutrino emission has been predicted for several short-lived astrophysical transients including gamma-ray bursts (GRBs), core-collapse supernovae with choked jets and neutron star mergers. IceCube's optical and X-ray follow-up program searches for such transient sources by looking for two or more muon neutrino candidates in directional coincidence and arriving within 100s. The measured rate of neutrino alerts is consistent with the expected rate of chance coincidences of atmospheric background events and no likely electromagnetic counterparts have been identified in Swift follow-up observations. Here, we calculate generic bounds on the neutrino flux of short-lived transient sources. Assuming an $E^{-2.5}$ neutrino spectrum, we find that the neutrino flux of rare sources, like long gamma-ray bursts, is constrained to <5% of the detected astrophysical flux and the energy released in neutrinos (100GeV to 10PeV) by a median bright GRB-like source is $<10^{52.5}$erg. For a harder $E^{-2.13}$ neutrino spectrum up to 30% of the flux could be produced by GRBs and the allowed median source energy is $< 10^{52}$erg. A hypothetical population of transient sources has to be more common than $10^{-5}\text{Mpc}^{-3}\text{yr}^{-1}$ ($5\times10^{-8}\text{Mpc}^{-3}\text{yr}^{-1}$ for the $E^{-2.13}$ spectrum) to account for the complete astrophysical neutrino flux.

Inelasticity--the fraction of a neutrino's energy transferred to hadrons--is a quantity of interest in the study of astrophysical and atmospheric neutrino interactions at multi-TeV energies with IceCube. In this work, a sample of contained neutrino interactions in IceCube is obtained from 5 years of data and classified as 2650 tracks and 965 cascades. Tracks arise predominantly from charged-current $\nu_{\mu}$ interactions, and we demonstrate that we can reconstruct their energy and inelasticity. The inelasticity distribution is found to be consistent with the calculation of Cooper-Sarkar et al. across the energy range from $\sim$ 1 TeV to $\sim$ 100 TeV. Along with cascades from neutrinos of all flavors, we also perform a fit over the energy, zenith angle, and inelasticity distribution to characterize the flux of astrophysical and atmospheric neutrinos. The energy spectrum of diffuse astrophysical neutrinos is well-described by a power-law in both track and cascade samples, and a best-fit index $\gamma=2.62\pm0.07$ is found in the energy range from 3.5 TeV to 2.6 PeV. Limits are set on the astrophysical flavor composition that are compatible with a ratio of $\left(\frac{1}{3}:\frac{1}{3}:\frac{1}{3}\right)_{\oplus}$. Exploiting the distinct inelasticity distribution of $\nu_{\mu}$ and $\bar{\nu}_{\mu}$ interactions, the atmospheric $\nu_{\mu}$ to $\bar{\nu}_{\mu}$ flux ratio in the energy range from 770 GeV to 21 TeV is found to be $0.77^{+0.44}_{-0.25}$ times the calculation by Honda et al. Lastly, the inelasticity distribution is also sensitive to neutrino charged-current charm production. The data are consistent with a leading-order calculation, with zero charm production excluded at $91\%$ confidence level. Future analyses of inelasticity distributions may probe new physics that affects neutrino interactions both in and beyond the Standard Model.

We report on the observation of a deficit in the cosmic ray flux from the directions of the Moon and Sun with five years of data taken by the IceCube Neutrino Observatory. Between May 2010 and May 2011 the IceCube detector operated with 79 strings deployed in the glacial ice at the South Pole, and with 86 strings between May 2011 and May 2015. A binned analysis is used to measure the relative deficit and significance of the cosmic ray shadows. Both the cosmic ray Moon and Sun shadows are detected with high statistical significance ($>10\sigma$) for each year. The results for the Moon shadow are consistent with previous analyses and verify the stability of the IceCube detector over time. This work represents the first observation of the Sun shadow with the IceCube detector. We show that the cosmic ray shadow of the Sun varies with time. These results open the possibility to study cosmic ray transport near the Sun with future data from IceCube.

The IceCube Collaboration has observed a high-energy astrophysical neutrino flux and recently found evidence for neutrino emission from the blazar TXS 0506+056. These results open a new window into the high-energy universe. However, the source or sources of most of the observed flux of astrophysical neutrinos remains uncertain. Here, a search for steady point-like neutrino sources is performed using an unbinned likelihood analysis. The method searches for a spatial accumulation of muon-neutrino events using the very high-statistics sample of about $497\,000$ neutrinos recorded by IceCube between 2009 and 2017. The median angular resolution is $\sim1^\circ$ at 1 TeV and improves to $\sim0.3^\circ$ for neutrinos with an energy of 1 PeV. Compared to previous analyses, this search is optimized for point-like neutrino emission with the same flux-characteristics as the observed astrophysical muon-neutrino flux and introduces an improved event-reconstruction and parametrization of the background. The result is an improvement in sensitivity to the muon-neutrino flux compared to the previous analysis of $\sim35\%$ assuming an $E^{-2}$ spectrum. The sensitivity on the muon-neutrino flux is at a level of $E^2 \mathrm{d} N /\mathrm{d} E = 3\cdot 10^{-13}\,\mathrm{TeV}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}$. No new evidence for neutrino sources is found in a full sky scan and in an a priori candidate source list that is motivated by gamma-ray observations. Furthermore, no significant excesses above background are found from populations of sub-threshold sources. The implications of the non-observation for potential source classes are discussed.

We present the first full-sky analysis of the cosmic ray arrival direction distribution with data collected by the HAWC 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 holds 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 $\delta_{0h} = 9.16 \times 10^{-4}$ and $\delta_{6h} = 7.25 \times 10^{-4}~(\pm0.04 \times 10^{-4})$. In addition, we infer the direction ($229.2\pm 3.5^\circ$ RA , $11.4\pm 3.0^\circ$ Dec.) 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 ^{+1.0}_{-2.0} \times 10^{-4}$.

The existence of diffuse Galactic neutrino production is expected from cosmic ray interactions with Galactic gas and radiation fields. Thus, neutrinos are a unique messenger offering the opportunity to test the products of Galactic cosmic ray interactions up to energies of hundreds of TeV. Here we present a search for this production using ten years of ANTARES track and shower data, as well as seven years of IceCube track data. The data are combined into a joint likelihood test for neutrino emission according to the KRA$_\gamma$ model assuming a 5 PeV per nucleon Galactic cosmic ray cutoff. No significant excess is found. As a consequence, the limits presented in this work start constraining the model parameter space for Galactic cosmic ray production and transport.

With the observation of high-energy astrophysical neutrinos by the IceCube Neutrino Observatory, interest has risen in models of PeV-mass decaying dark matter particles to explain the observed flux. We present two dedicated experimental analyses to test this hypothesis. One analysis uses six years of IceCube data focusing on muon neutrino 'track' events from the Northern Hemisphere, while the second analysis uses two years of 'cascade' events from the full sky. Known background components and the hypothetical flux from unstable dark matter are fitted to the experimental data. Since no significant excess is observed in either analysis, lower limits on the lifetime of dark matter particles are derived: We obtain the strongest constraint to date, excluding lifetimes shorter than $10^{28}\,$s at $90\%$ CL for dark matter masses above $10\,$TeV.

We report a quasi-differential upper limit on the extremely-high-energy (EHE) neutrino flux above $5\times 10^{6}$ 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 $10^{6}$ 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 \times 10^{6}$ and $5 \times 10^{10}$ GeV. This result constrains neutrino models predicting a three-flavor neutrino flux of $E_\nu^2\phi_{\nu_e+\nu_\mu+\nu_\tau}\simeq2\times 10^{-8}\ {\rm GeV}/{\rm cm}^2\ \sec\ {\rm sr}$ at $10^9\ {\rm 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 excluded.

We report constraints on the sources of ultra-high-energy cosmic ray (UHECR) above $10^{9}$ GeV, based on an analysis of seven years of IceCube data. This analysis efficiently selects very high energy neutrino-induced events which have deposited energies from $\sim 10^6$ GeV to above $10^{11}$ GeV. Two neutrino-induced events with an estimated deposited energy of $(2.6 \pm 0.3) \times 10^6$ GeV, the highest neutrino energies observed so far, and $(7.7 \pm 2.0) \times 10^5$ GeV were detected. The atmospheric background-only hypothesis of detecting these events is rejected at 3.6$\sigma$. The hypothesis that the observed events are of cosmogenic origin is also rejected at $>$99% CL because of the limited deposited energy and the non-observation of events at higher energy, while their observation is consistent with an astrophysical origin. Our limits on cosmogenic neutrino fluxes disfavor the UHECR sources having cosmological evolution stronger than the star formation rate, e.g., active galactic nuclei and $\gamma$-ray bursts, assuming proton-dominated UHECRs. Constraints on UHECR sources including mixed and heavy UHECR compositions are obtained for models of neutrino production within UHECR sources. Our limit disfavors a significant part of parameter space for active galactic nuclei and new-born pulsar models.

We present a search for coincidence between IceCube TeV neutrinos and fast radio bursts (FRBs). During the search period from 2010 May 31 to 2016 May 12, a total of 29 FRBs with 13 unique locations have been detected in the whole sky. An unbinned maximum likelihood method was used to search for spatial and temporal coincidence between neutrinos and FRBs in expanding time windows, in both the northern and southern hemispheres. No significant correlation was found in six years of IceCube data. Therefore, we set upper limits on neutrino fluence emitted by FRBs as a function of time window duration. We set the most stringent limit obtained to date on neutrino fluence from FRBs with an $E^{-2}$ energy spectrum assumed, which is 0.0021 GeV cm$^{-2}$ per burst for emission timescales up to \textasciitilde10$^2$ seconds from the northern hemisphere stacking search.

Neutrinos interact only very weakly, so they are extremely penetrating. However, the theoretical neutrino-nucleon interaction cross section rises with energy such that, at energies above 40 TeV, neutrinos are expected to be absorbed as they pass through the Earth. Experimentally, the cross section has been measured only at the relatively low energies (below 400 GeV) available at neutrino beams from accelerators \cite{Agashe:2014kda, Formaggio:2013kya}. Here we report the first measurement of neutrino absorption in the Earth, using a sample of 10,784 energetic upward-going neutrino-induced muons observed with the IceCube Neutrino Observatory. The flux of high-energy neutrinos transiting long paths through the Earth is attenuated compared to a reference sample that follows shorter trajectories through the Earth. Using a fit to the two-dimensional distribution of muon energy and zenith angle, we determine the cross section for neutrino energies between 6.3 TeV and 980 TeV, more than an order of magnitude higher in energy than previous measurements. The measured cross section is $1.30^{+0.21}_{-0.19}$ (stat.) $^{+0.39}_{-0.43}$ (syst.) times the prediction of the Standard Model \cite{CooperSarkar:2011pa}, consistent with the expectation for charged and neutral current interactions. We do not observe a dramatic increase in the cross section, expected in some speculative models, including those invoking new compact dimensions \cite{AlvarezMuniz:2002ga} or the production of leptoquarks \cite{Romero:2009vu}.

We present an all-sky search for muon neutrinos produced during the prompt $\gamma$-ray emission of 1172 gamma-ray bursts (GRBs) with the IceCube Neutrino Observatory. The detection of these neutrinos would constitute evidence for ultra-high energy cosmic ray (UHECR) production in GRBs, as interactions between accelerated protons and the prompt $\gamma$-ray field would yield charged pions, which decay to neutrinos. A previously reported search for muon neutrino tracks from Northern Hemisphere GRBs has been extended to include three additional years of IceCube data. A search for such tracks from Southern Hemisphere GRBs in five years of IceCube data has been introduced to enhance our sensitivity to the highest energy neutrinos. No significant correlation between neutrino events and observed GRBs is seen in the new data. Combining this result with previous muon neutrino track searches and a search for cascade signature events from all neutrino flavors, we obtain new constraints for single-zone fireball models of GRB neutrino and UHECR production.

Papers on the properties of the atmospheric and astrophysical neutrino flux submitted to the 35th International Cosmic Ray Conference (ICRC 2017, Busan, South Korea) by the IceCube Collaboration

Papers on research & development towards IceCube-Gen2, the next generation neutrino observatory at South Pole, submitted to the 35th International Cosmic Ray Conference (ICRC 2017, Busan, South Korea) by the IceCube-Gen2 Collaboration.

Papers on solar flares, supernovae, event reconstruction and education & outreach, submitted to the 35th International Cosmic Ray Conference (ICRC 2017, Busan, South Korea) by the IceCube Collaboration