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.

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 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%.

High-energy neutrinos are unique messengers of the high-energy universe, tracing the processes of cosmic ray acceleration. This paper presents analyses focusing on time-dependent neutrino point-source searches. A scan of the whole sky, making no prior assumption about source candidates, is performed, looking for a space and time clustering of high-energy neutrinos in data collected by the IceCube Neutrino Observatory between 2012 and 2017. No statistically significant evidence for a time-dependent neutrino signal is found with this search during this period, as all results are consistent with the background expectation. Within this study period, the blazar 3C 279, showed strong variability, inducing a very prominent gamma-ray flare observed in 2015 June. This event motivated a dedicated study of the blazar, which consists of searching for a time-dependent neutrino signal correlated with the gamma-ray emission. No evidence for a time-dependent signal is found. Hence, an upper limit on the neutrino fluence is derived, allowing us to constrain a hadronic emission model.

The High Altitude Water Cerenkov (HAWC) and IceCube observatories, through the Astrophysical Multimessenger Observatory Network (AMON) framework, have developed a multimessenger joint search for extragalactic astrophysical sources. This analysis looks for sources that emit both cosmic neutrinos and gamma rays that are produced in photohadronic or hadronic interactions. The AMON system is running continuously, receiving subthreshold data (i.e., data that are not suited on their own to do astrophysical searches) from HAWC and IceCube, and combining them in real time. Here we present the analysis algorithm, as well as results from archival data collected between 2015 June and 2018 August, with a total live time of 3.0 yr. During this period we found two coincident events that have a false-alarm rate (FAR) of <1 coincidence yr-1, consistent with the background expectations. The real-time implementation of the analysis in the AMON system began on 2019 November 20 and issues alerts to the community through the Gamma-ray Coordinates Network with an FAR threshold of <4 coincidences yr-1.

We present a high-energy neutrino event generator, called LeptonInjector, alongside an event weighter, called LeptonWeighter. Both are designed for large-volume Cherenkov neutrino telescopes such as IceCube. The neutrino event generator allows for quick and flexible simulation of neutrino events within and around the detector volume, and implements the leading Standard Model neutrino interaction processes relevant for neutrino observatories: neutrino-nucleon deep-inelastic scattering and neutrino-electron annihilation. In this paper, we discuss the event generation algorithm, the weighting algorithm, and the main functions of the publicly available code, with examples. Program summary: Program Titles: LeptonInjector and LeptonWeighter CPC Library link to program files: https://doi.org/10.17632/662gkpjfd9.1 Developer's repository links: https://github.com/icecube/LeptonInjector and https://github.com/icecube/LeptonWeighter Licensing provisions: GNU Lesser General Public License, version 3. Programming Language: C++11 External Routines: • Boost • HDF5 • nuflux (https://github.com/icecube/nuflux) • nuSQuIDS (https://github.com/arguelles/nuSQuIDS) • Photospline (https://github.com/icecube/photospline) • SuiteSparse (https://github.com/DrTimothyAldenDavis/SuiteSparse) Nature of problem: LeptonInjector: Generate neutrino interaction events of all possible topologies and energies throughout and around a detector volume. LeptonWeighter: Reweight Monte Carlo events, generated by a set of LeptonInjector Generators, to any desired physical neutrino flux or cross section. Solution method: LeptonInjector: Projected ranges of generated leptons and the extent of the detector, in terms of column depth, are used to inject events in and around the detector volume. Event kinematics follow distributions provided in cross section files. LeptonWeighter: Event generation probabilities are calculated for each Generator, which are then combined into a generation weight and used to calculate an overall event weight.

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 data set is particularly well suited. For many of the signal candidates considered, this analysis is the most sensitive of any experiment to date. 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.