A measurement of the diffuse astrophysical neutrino spectrum is presented using IceCube data collected from 2011-2022 (10.3 years). We developed novel detection techniques to search for events with a contained vertex and exiting track induced by muon neutrinos undergoing a charged-current interaction. Searching for these starting track events allows us to not only more effectively reject atmospheric muons but also atmospheric neutrino backgrounds in the southern sky, opening a new window to the sub-100 TeV astrophysical neutrino sky. The event selection is constructed using a dynamic starting track veto and machine learning algorithms. We use this data to measure the astrophysical diffuse flux as a single power law flux (SPL) with a best-fit spectral index of $\gamma = 2.58 ^{+0.10}_{-0.09}$ and per-flavor normalization of $\phi^{\mathrm{Astro}}_{\mathrm{per-flavor}} = 1.68 ^{+0.19}_{-0.22} \times 10^{-18} \times \mathrm{GeV}^{-1} \mathrm{cm}^{-2} \mathrm{s}^{-1} \mathrm{sr}^{-1}$ (at 100 TeV). The sensitive energy range for this dataset is 3 - 550 TeV under the SPL assumption. This data was also used to measure the flux under a broken power law, however we did not find any evidence of a low energy cutoff.

Gamma-ray bursts (GRBs) have long been considered a possible source of high-energy neutrinos. While no correlations have yet been detected between high-energy neutrinos and GRBs, the recent observation of GRB 221009A - the brightest GRB observed by Fermi-GBM to date and the first one to be observed above an energy of 10 TeV - provides a unique opportunity to test for hadronic emission. In this paper, we leverage the wide energy range of the IceCube Neutrino Observatory to search for neutrinos from GRB 221009A. We find no significant deviation from background expectation across event samples ranging from MeV to PeV energies, placing stringent upper limits on the neutrino emission from this source.

Abstract Atmospheric muon neutrinos are produced by meson decays in cosmic-ray-induced air showers. The flux depends on meteorological quantities such as the air temperature, which affects the density of air. Competition between decay and re-interaction of those mesons in the first particle production generations gives rise to a higher neutrino flux when the air density in the stratosphere is lower, corresponding to a higher temperature. A measurement of a temperature dependence of the atmospheric $$\nu _{\mu }$$ ν μ flux provides a novel method for constraining hadronic interaction models of air showers. It is particularly sensitive to the production of kaons. Studying this temperature dependence for the first time requires a large sample of high-energy neutrinos as well as a detailed understanding of atmospheric properties. We report the significant ( $$> 10 \; \sigma $$ > 10 σ ) observation of a correlation between the rate of more than 260,000 neutrinos, detected by IceCube between 2012 and 2018, and atmospheric temperatures of the stratosphere, measured by the Atmospheric Infrared Sounder (AIRS) instrument aboard NASA’s AQUA satellite. For the observed 10 $$\%$$ % seasonal change of effective atmospheric temperature we measure a 3.5(3) $$\%$$ % change in the muon neutrino flux. This observed correlation deviates by about 2-3 standard deviations from the expected correlation of 4.3 $$\%$$ % as obtained from theoretical predictions under the assumption of various hadronic interaction models.

Core-collapse supernovae are a promising potential high-energy neutrino source class. We test for correlation between seven years of IceCube neutrino data and a catalog containing more than 1000 core-collapse supernovae of types IIn and IIP and a sample of stripped-envelope supernovae. We search both for neutrino emission from individual supernovae as well as for combined emission from the whole supernova sample, through a stacking analysis. No significant spatial or temporal correlation of neutrinos with the cataloged supernovae was found. All scenarios were tested against the background expectation and together yield an overall p -value of 93%; therefore, they show consistency with the background only. The derived upper limits on the total energy emitted in neutrinos are 1.7 × 10 48 erg for stripped-envelope supernovae, 2.8 × 10 48 erg for type IIP, and 1.3 × 10 49 erg for type IIn SNe, the latter disfavoring models with optimistic assumptions for neutrino production in interacting supernovae. We conclude that stripped-envelope supernovae and supernovae of type IIn do not contribute more than 14.6% and 33.9%, respectively, to the diffuse neutrino flux in the energy range of about [ 10 3 –10 5 ] GeV, assuming that the neutrino energy spectrum follows a power-law with an index of −2.5. Under the same assumption, we can only constrain the contribution of type IIP SNe to no more than 59.9%. Thus, core-collapse supernovae of types IIn and stripped-envelope supernovae can both be ruled out as the dominant source of the diffuse neutrino flux under the given assumptions.

We present a catalog of likely astrophysical neutrino track-like events from the IceCube Neutrino Observatory. IceCube began reporting likely astrophysical neutrinos in 2016, and this system was updated in 2019. The catalog presented here includes events that were reported in real time since 2019, as well as events identified in archival data samples starting from 2011. We report 275 neutrino events from two selection channels as the first entries in the catalog, the IceCube Event Catalog of Alert Tracks, which will see ongoing extensions with additional alerts. The Gold and Bronze alert channels respectively provide neutrino candidates with a 50% and 30% probability of being astrophysical, on average assuming an astrophysical neutrino power-law energy spectral index of 2.19. For each neutrino alert, we provide the reconstructed energy, direction, false-alarm rate, probability of being astrophysical in origin, and likelihood contours describing the spatial uncertainty in the alert's reconstructed location. We also investigate a directional correlation of these neutrino events with gamma-ray and X-ray catalogs, including 4FGL, 3HWC, TeVCat, and Swift-BAT.

This paper presents the results of a search for neutrinos that are spatially and temporally coincident with 22 unique, nonrepeating fast radio bursts (FRBs) and one repeating FRB (FRB 121102). FRBs are a rapidly growing class of Galactic and extragalactic astrophysical objects that are considered a potential source of high-energy neutrinos. The IceCube Neutrino Observatory's previous FRB analyses have solely used track events. This search utilizes seven years of IceCube cascade events which are statistically independent of track events. This event selection allows probing of a longer range of extended timescales due to the low background rate. No statistically significant clustering of neutrinos was observed. Upper limits are set on the time-integrated neutrino flux emitted by FRBs for a range of extended time windows.

Using data from the IceCube Neutrino Observatory, we searched for high-energy neutrino emission from the gravitational-wave events detected by the advanced LIGO and Virgo detectors during their third observing run. We did a low-latency follow-up on the public candidate events released during the detectors’ third observing run and an archival search on the 80 confident events reported in the GWTC-2.1 and GWTC-3 catalogs. An extended search was also conducted for neutrino emission on longer timescales from neutron star containing mergers. Follow-up searches on the candidate optical counterpart of GW190521 were also conducted. We used two methods; an unbinned maximum likelihood analysis and a Bayesian analysis using astrophysical priors, both of which were previously used to search for high-energy neutrino emission from gravitational-wave events. No significant neutrino emission was observed by any analysis, and upper limits were placed on the time-integrated neutrino flux as well as the total isotropic equivalent energy emitted in high-energy neutrinos.

The D-Egg, an acronym for “Dual optical sensors in an Ellipsoid Glass for Gen2,” is one of the optical modules designed for future extensions of the IceCube experiment at the South Pole. The D-Egg has an elongated-sphere shape to maximize the photon-sensitive effective area while maintaining a narrow diameter to reduce the cost and the time needed for drilling of the deployment holes in the glacial ice for the optical modules at depths up to 2700 m. The D-Egg design is utilized for the IceCube Upgrade, the next stage of the IceCube project also known as IceCube-Gen2 Phase 1, where nearly half of the optical sensors to be deployed are D-Eggs. With two 8-inch high-quantum efficiency photomultiplier tubes (PMTs) per module, D-Eggs offer an increased effective area while retaining the successful design of the IceCube digital optical module (DOM). The convolution of the wavelength-dependent effective area and the Cherenkov emission spectrum provides an effective photodetection sensitivity that is 2.8 times larger than that of IceCube DOMs. The signal of each of the two PMTs is digitized using ultra-low-power 14-bit analog-to-digital converters with a sampling frequency of 240 MSPS, enabling a flexible event triggering, as well as seamless and lossless event recording of single-photon signals to multi-photons exceeding 200 photoelectrons within 10 ns. Mass production of D-Eggs has been completed, with 277 out of the 310 D-Eggs produced to be used in the IceCube Upgrade. In this paper, we report the design of the D-Eggs, as well as the sensitivity and the single to multi-photon detection performance of mass-produced D-Eggs measured in a laboratory using the built-in data acquisition system in each D-Egg optical sensor module.

The understanding of novae, the thermonuclear eruptions on the surfaces of white dwarf stars in binaries, has recently undergone a major paradigm shift. Though the bolometric luminosity of novae was long thought to arise directly from photons supplied by the thermonuclear runaway, recent GeV gamma-ray observations have supported the notion that a significant portion of the luminosity could come from radiative shocks. More recently, observations of novae have lent evidence that these shocks are acceleration sites for hadrons for at least some types of novae. In this scenario, a flux of neutrinos may accompany the observed gamma rays. As the gamma rays from most novae have only been observed up to a few GeV, novae have previously not been considered as targets for neutrino telescopes, which are most sensitive at and above TeV energies. Here, we present the first search for neutrinos from novae with energies between a few GeV and 10 TeV using IceCube-DeepCore, a densely instrumented region of the IceCube Neutrino Observatory with a reduced energy threshold. We search both for a correlation between gamma-ray and neutrino emission as well as between optical and neutrino emission from novae. We find no evidence for neutrino emission from the novae considered in this analysis and set upper limits for all gamma-ray detected novae.

Dark Matter particles in the Galactic Center and halo can annihilate or decay into a pair of neutrinos producing a monochromatic flux of neutrinos. The spectral feature of this signal is unique and it is not expected from any astrophysical production mechanism. Its observation would constitute a dark matter smoking gun signal. We performed the first dedicated search with a neutrino telescope for such signal, by looking at both the angular and energy information of the neutrino events. To this end, a total of five years of IceCube's DeepCore data has been used to test dark matter masses ranging from 10~GeV to 40~TeV. No significant neutrino excess was found and upper limits on the annihilation cross section, as well as lower limits on the dark matter lifetime, were set. The limits reached are of the order of $10^{-24}$~cm$^3/s$ for an annihilation and up to $10^{27}$ seconds for decaying Dark Matter. Using the same data sample we also derive limits for dark matter annihilation or decay into a pair of Standard Model charged particles.

Beginning in 2016, the IceCube Neutrino Observatory has sent out alerts in real time containing the information of high-energy ($E \gtrsim 100$~TeV) neutrino candidate events with moderate-to-high ($\gtrsim 30$\%) probability of astrophysical origin. In this work, we use a recent catalog of such alert events, which, in addition to events announced in real-time, includes events that were identified retroactively, and covers the time period of 2011-2020. We also search for additional, lower-energy, neutrinos from the arrival directions of these IceCube alerts. We show how performing such an analysis can constrain the contribution of rare populations of cosmic neutrino sources to the diffuse astrophysical neutrino flux. After searching for neutrino emission coincident with these alert events on various timescales, we find no significant evidence of either minute-scale or day-scale transient neutrino emission or of steady neutrino emission in the direction of these alert events. This study also shows how numerous a population of neutrino sources has to be to account for the complete astrophysical neutrino flux. Assuming sources have the same luminosity, an $E^{-2.5}$ neutrino spectrum and number densities that follow star-formation rates, the population of sources has to be more numerous than $7\times 10^{-9}~\textrm{Mpc}^{-3}$. This number changes to $3\times 10^{-7}~\textrm{Mpc}^{-3}$ if number densities instead have no cosmic evolution.

Galactic PeV cosmic-ray accelerators (PeVatrons) are Galactic sources theorized to accelerate cosmic rays up to PeV in energy. The accelerated cosmic rays are expected to interact hadronically with nearby ambient gas or the interstellar medium, resulting in γ -rays and neutrinos. Recently, the Large High Altitude Air Shower Observatory (LHAASO) identified 12 γ -ray sources with emissions above 100 TeV, making them candidates for PeVatrons. While at these high energies the Klein–Nishina effect exponentially suppresses leptonic emission from Galactic sources, evidence for neutrino emission would unequivocally confirm hadronic acceleration. Here, we present the results of a search for neutrinos from these γ -ray sources and stacking searches testing for excess neutrino emission from all 12 sources as well as their subcatalogs of supernova remnants and pulsar wind nebulae with 11 yr of track events from the IceCube Neutrino Observatory. No significant emissions were found. Based on the resulting limits, we place constraints on the fraction of γ -ray flux originating from the hadronic processes in the Crab Nebula and LHAASO J2226+6057.

Recent observations of GeV gamma-rays from novae have led to a paradigm shift in the understanding of these objects. While it is now believed that shocks contribute significantly to the energy budget of novae, it is still unknown if the emission is hadronic or leptonic in origin. Neutrinos could hold the key to definitively differentiating between these two scenarios, though the energies of such particles would be much lower than are typically targeted with neutrino telescopes. IceCube's densely instrumented DeepCore sub-array provides the ability to reduce the threshold for observation from 1 TeV down to approximately 10 GeV. We will discuss recent measurements in this low energy regime, details of a new sub-TeV selection, and prospects for future searches for transient neutrino emission.

Measurements of neutrinos at and below 10 GeV provide unique constraints of neutrino oscillation parameters as well as probes of potential Non-Standard Interactions (NSI). The IceCube Neutrino Observatory's DeepCore array is designed to detect neutrinos down to GeV energies. IceCube has built the world's largest data set of neutrinos >10 GeV, making searches for NSI a computational challenge. This work describes the use of convolutional neural networks (CNNs) to improve the energy reconstruction resolution and speed of reconstructing O(10 GeV) neutrino events in IceCube. Compared to current likelihood-based methods which take seconds to minutes, the CNN is expected to provide approximately a factor of 2 improvement in energy resolution while reducing the reconstruction time per event to milliseconds, which is essential for processing large datasets.

New optical sensors called the "D-Egg" have been developed for cost-effective instrumentation for the IceCube Upgrade. With two 8-inch high quantum efficient photomultiplier tubes (PMTs), they offer increased effective photocathode area while retaining as much of the successful IceCube Digital Optical Module design as possible. Mass production of D-Eggs has started in 2020. By the end of 2021, there will be 310 D-Eggs produced with 288 deployed in the IceCube Upgrade. The D-Egg readout system uses advanced technologies in electronics and computing power. Each of the two PMT signals is digitised using ultra-low-power 14-bit ADCs with a sampling frequency of 240 MSPS, enabling seamless and lossless event recording from single-photon signals to signals exceeding 200 PE within 10 ns, as well as flexible event triggering. In this paper, we report the single photon detection performance as well as the multiple photon recording capability of D-Eggs from the mass production line which have been evaluated with the built-in data acquisition system.

IceCube, a cubic-kilometer array of optical sensors built to detect atmospheric and astrophysical neutrinos between 1 GeV and 1 PeV, is deployed 1.45 km to 2.45 km below the surface of the ice sheet at the South Pole. The classification and reconstruction of events from the in-ice detectors play a central role in the analysis of data from IceCube. Reconstructing and classifying events is a challenge due to the irregular detector geometry, inhomogeneous scattering and absorption of light in the ice and, below 100 GeV, the relatively low number of signal photons produced per event. To address this challenge, it is possible to represent IceCube events as point cloud graphs and use a Graph Neural Network (GNN) as the classification and reconstruction method. The GNN is capable of distinguishing neutrino events from cosmic-ray backgrounds, classifying different neutrino event types, and reconstructing the deposited energy, direction and interaction vertex. Based on simulation, we provide a comparison in the 1-100 GeV energy range to the current state-of-the-art maximum likelihood techniques used in current IceCube analyses, including the effects of known systematic uncertainties. For neutrino event classification, the GNN increases the signal efficiency by 18% at a fixed false positive rate (FPR), compared to current IceCube methods. Alternatively, the GNN offers a reduction of the FPR by over a factor 8 (to below half a percent) at a fixed signal efficiency. For the reconstruction of energy, direction, and interaction vertex, the resolution improves by an average of 13%-20% compared to current maximum likelihood techniques in the energy range of 1-30 GeV. The GNN, when run on a GPU, is capable of processing IceCube events at a rate nearly double of the median IceCube trigger rate of 2.7 kHz, which opens the possibility of using low energy neutrinos in online searches for transient events.