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 ∼102 s from the northern hemisphere stacking search.

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 ∼ 1 ∘ at 1 TeV and improves to ∼ 0. 3 ∘ 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 ∼ 35 % assuming an E- 2 spectrum. The sensitivity on the muon-neutrino flux is at a level of E2dN/dE=3·10-13TeVcm-2s-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.

Adopting the Standard Halo Model (SHM) of an isotropic Maxwellian velocity distribution for dark matter (DM) particles in the Galaxy, the most stringent current constraints on their spin-dependent scattering cross-section with nucleons come from the IceCube neutrino observatory and the PICO-60 C 3F 8 superheated bubble chamber experiments. The former is sensitive to high energy neutrinos from the self-annihilation of DM particles captured in the Sun, while the latter looks for nuclear recoil events from DM scattering off nucleons. Although slower DM particles are more likely to be captured by the Sun, the faster ones are more likely to be detected by PICO. Recent N-body simulations suggest significant deviations from the SHM for the smooth halo component of the DM, while observations hint at a dominant fraction of the local DM being in substructures. We use the method of Ferrer et al. (JCAP 1509: 052, 2015) to exploit the complementarity between the two approaches and derive conservative constraints on DM-nucleon scattering. Our results constrain σSD≲ 3 × 10 - 39cm 2 (6 × 10 - 38cm 2) at ≳ 90 % C.L. for a DM particle of mass 1 TeV annihilating into τ+τ- (bb¯) with a local density of ρDM=0.3GeV/cm3. The constraints scale inversely with ρDM and are independent of the DM velocity distribution.

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 (229.°2 ± 3.°5 R.A., 11.°4 ± 3.°0 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 δN ∼ -3.97+1.0-2.0 × 10-4.

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 expected distributions of observables 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.

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 2010 May and 2011 May the IceCube detector operated with 79 strings deployed in the glacial ice at the South Pole, and with 86 strings between 2011 May and 2015 May. 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σ) 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 make it possible to study cosmic ray transport near the Sun with future data from IceCube.

We present two searches for IceCube neutrino events coincident with 28 fast radio bursts (FRBs) and 1 repeating FRB. The first improves on a previous IceCube analysis - searching for spatial and temporal correlation of events with FRBs at energies greater than roughly 50 GeV - by increasing the effective area by an order of magnitude. The second is a search for temporal correlation of MeV neutrino events with FRBs. No significant correlation is found in either search; therefore, we set upper limits on the time-integrated neutrino flux emitted by FRBs for a range of emission timescales less than one day. These are the first limits on FRB neutrino emission at the MeV scale, and the limits set at higher energies are an order-of-magnitude improvement over those set by any neutrino telescope.

In multi-messenger astronomy, rapid investigation of interesting transients is imperative. As an observatory with a 4π steradian field of view, and ∼99% uptime, the IceCube Neutrino Observatory is a unique facility to follow up transients, as well as to provide valuable insights for other observatories and inform their observational decisions. Since 2016, IceCube has been using low-latency data to rapidly respond to interesting astrophysical events reported by the multi-messenger observational community. Here, we describe the pipeline used to perform these followup analyses, and provide a summary of the 58 analyses performed as of July 2020. We find no significant signal in the first 58 analyses performed. The pipeline has helped inform various electromagnetic observation strategies, and has constrained neutrino emission from potential hadronic cosmic accelerators.

The measurement of diffuse PeV gamma-ray emission from the Galactic plane would provide information about the energy spectrum and propagation of Galactic cosmic rays, and the detection of a pointlike source of PeV gamma-rays would be strong evidence for a Galactic source capable of accelerating cosmic rays up to at least a few PeV. This paper presents several unbinned maximum-likelihood searches for PeV gamma-rays in the Southern Hemisphere using 5 yr of data from the IceTop air shower surface detector and the in-ice array of the IceCube Observatory. The combination of both detectors takes advantage of the low muon content and deep shower maximum of gamma-ray air showers and provides excellent sensitivity to gamma-rays between ∼0.6 and 100 PeV. Our measurements of pointlike and diffuse Galactic emission of PeV gamma-rays are consistent with the background, so we constrain the angle-integrated diffuse gamma-ray flux from the Galactic plane at 2 PeV to 2.61 × 10-19 cm-2 s-1 TeV-1 at 90% confidence, assuming an E -3 spectrum, and we estimate 90% upper limits on pointlike emission at 2 PeV between 10-21 and 10-20 cm-2 s-1 TeV-1 for an E -2 spectrum, depending on decl. Furthermore, we exclude unbroken power-law emission up to 2 PeV for several TeV gamma-ray sources observed by the High Energy Spectroscopic System and calculate upper limits on the energy cutoffs of these sources at 90% confidence. We also find no PeV gamma-rays correlated with neutrinos from IceCube's high-energy starting event sample. These are currently the strongest constraints on PeV gamma-ray emission.

The results of a 3+1 sterile neutrino search using eight years of data from the IceCube Neutrino Observatory are presented. A total of 305 735 muon neutrino events are analyzed in reconstructed energy-zenith space to test for signatures of a matter-enhanced oscillation that would occur given a sterile neutrino state with a mass-squared differences between 0.01 and 100  eV^{2}. The best-fit point is found to be at sin^{2}(2θ_{24})=0.10 and Δm_{41}^{2}=4.5  eV^{2}, which is consistent with the no sterile neutrino hypothesis with a p value of 8.0%.