Experiments designed to detect ultra-high energy (UHE) neutrinos using radio techniques are also capable of detecting the radio signals from cosmic-ray (CR) induced air showers. These CR signals are important both as a background and as a tool for calibrating the detector. The Askaryan Radio Array (ARA), a radio detector array, is designed to detect UHE neutrinos. The array currently comprises five independent stations, each instrumented with antennas deployed at depths of up to 200 meters within the ice at the South Pole. In this study, we focus on a candidate event recorded by ARA Station 2 (A2) that shows features consistent with a downward-going CR-induced air shower. This includes distinctive double-pulse signals in multiple channels, interpreted as geomagnetic and Askaryan radio emissions arriving at the antennas in sequence. To investigate this event, we use detailed simulations that combine a modern ice-impacting CR shower simulation framework, FAERIE, with a realistic detector simulation package, AraSim. We will present results for an optimization of the event topology, identified through simulated CR showers, comparing the vertex reconstruction of both the geomagnetic and Askaryan signals of the event, as well as the observed time delays between the two signals in each antenna.

We analyzed the 7.92 × 10 ^11 cosmic-ray-induced muon events collected by the IceCube Neutrino Observatory from 2011 May 13, when the fully constructed experiment started to take data, to 2023 May 12. This data set provides an up-to-date cosmic-ray arrival direction distribution in the Southern Hemisphere with unprecedented statistical accuracy covering more than a full period length of a solar cycle. Improvements in Monte Carlo event simulation and better handling of year-to-year differences in data processing significantly reduce systematic uncertainties below the level of statistical fluctuations compared to the previously published results. We confirm the observation of a change in the angular structure of the cosmic-ray anisotropy between 10 TeV and 1 PeV, more specifically in the 100–300 TeV energy range. For the first time, we analyzed the angular power spectrum at different energies. The observed variations of the power spectra with energy suggest relatively reduced large-scale features at high energy compared to those of medium and small scales. The large volume of data enhances the statistical significance at higher energies, up to the PeV scale, and smaller angular scales, down to approximately 6° compared to previous findings.

While the sources of the diffuse astrophysical neutrino flux detected by the IceCube Neutrino Observatory are still largely unknown, one of the promising methods to improve our understanding of them is investigating the potential temporal and spatial correlations between neutrino alerts and the electromagnetic radiation from blazars. We report on the multiwavelength target-of-opportunity observations of the blazar B3 2247+381, taken in response to an IceCube multiplet alert for a cluster of muon neutrino events compatible with the source location between 2022 May 20 and 2022 November 10. B3 2247+381 was not detected with VERITAS during this time period. The source was found to be in a low-flux state in the optical, ultraviolet, and gamma-ray bands for the time interval corresponding to the neutrino event, but was detected in the hard X-ray band with NuSTAR during this period. We find the multiwavelength spectral energy distribution is described well using a simple one-zone leptonic synchrotron self-Compton radiation model. Moreover, assuming the neutrinos originate from hadronic processes within the jet, the neutrino flux would be accompanied by a photon flux from the cascade emission, and the integrated photon flux required in such a case would significantly exceed the total multiwavelength fluxes and the VERITAS upper limits presented here. The lack of flaring activity observed with VERITAS, combined with the low multiwavelength flux levels, as well as the significance of the neutrino excess being at a 3 σ level (uncorrected for trials), makes B3 2247+381 an unlikely source of the IceCube multiplet. We conclude that the neutrino excess is likely a background fluctuation.

We report a search for high-energy astrophysical neutrino multiplets, detections of multiple neutrino clusters in the same direction within 30 days, based on an analysis of 11.4 yr of IceCube data. A new search method optimized for transient neutrino emission with a monthly timescale is employed, providing a higher sensitivity to neutrino fluxes. This result is sensitive to neutrino transient emission, reaching per-flavor flux of approximately 1 0 − 10 erg cm − 2 s − 1 from the Northern Sky in the energy range E ≳ 50 TeV. The number of doublets and triplets identified in this search is compatible with the atmospheric background hypothesis, which leads us to set limits on the nature of neutrino transient sources with emission timescales of one month.

Active galactic nuclei (AGN) are promising candidate sources of high-energy astrophysical neutrinos, since they provide environments rich in matter and photon targets where cosmic-ray interactions may lead to the production of gamma rays and neutrinos. We searched for high-energy neutrino emission from AGN using the Swift-BAT Spectroscopic Survey catalog of hard X-ray sources and 12 yr of IceCube muon track data. First, upon performing a stacked search, no significant emission was found. Second, we searched for neutrinos from a list of 43 candidate sources and found an excess from the direction of two sources, the Seyfert galaxies NGC 1068 and NGC 4151. We observed NGC 1068 at flux ϕ ν μ + ν ¯ μ = 4.0 2 − 1.52 + 1.58 × 1 0 − 11 TeV −1 cm −2 s −1 normalized at 1 TeV, with a power-law spectral index γ = 3.10 − 0.22 + 0.26 , consistent with previous IceCube results. The observation of a neutrino excess from the direction of NGC 4151 is at a posttrial significance of 2.9 σ . If interpreted as an astrophysical signal, the excess observed from NGC 4151 corresponds to a flux ϕ ν μ + ν ¯ μ = 1.5 1 − 0.81 + 0.99 × 1 0 − 11 TeV −1 cm −2 s −1 normalized at 1 TeV and γ = 2.83 − 0.28 + 0.35 .

The Askaryan Radio Array (ARA), located near the South Pole, is among the first experiments there designed to detect ultra-high energy neutrinos through the Askaryan effect. When such neutrinos interact in dense media like ice, they initiate particle cascades that generate coherent radio pulses. Operating in the 150-850 MHz band, ARA is deployed 80-200 meters deep in Antarctic ice, where the radio background is exceptionally low. Despite this, experiments such as ARA must still account for continuous wave (CW) signals, which can originate from anthropogenic sources, instrumental noise, and other environmental factors. These CW signals can obscure faint neutrino-induced pulses, complicating data analysis and event identification. Over the years, ARA has developed and refined techniques for CW filtering and identification, including spectral analysis, notch filtering, and phase-variance methods. These approaches exploit CW characteristics such as their narrowband nature and persistence to separate contamination from genuine impulsive events. We review the main CW identification and filtering techniques developed within the ARA collaboration and present recent improvements in adaptive, multi-stage filtering pipelines. These advances have led to faster processing, simpler operation, and more accurate CW suppression, improving data analysis quality. Their efficacy is demonstrated across all ARA stations, underscoring their role in reducing event misclassification and enhancing performance. By refining these techniques, this work not only improves the sensitivity of ARA but also highlights the importance of robust CW filtering for current and future neutrino radio detection experiments.

We summarize the motivation for (as well as the presentations from) the April, 2024 workshop held in India and focused on radiowave techniques for cosmic ray and neutrino detection.

The Radio Neutrino Observatory – Greenland (RNO-G) seeks discovery of ultra-high energy neutrinos from the cosmos through their interactions in ice. The science program extends beyond particle astrophysics to include radioglaciology and, as we show herein, solar observations, as well. Currently seven of 35 planned radio-receiver stations (24 antennas/station) are operational. These stations are sensitive to impulsive radio signals with frequencies between 80 and 700 MHz and feature a neutrino trigger threshold for recording data close to the thermal floor. RNO-G can also trigger on elevated signals from the Sun, resulting in nanosecond resolution time-domain flare data; such temporal resolution is significantly shorter than from most dedicated solar observatories. In addition to possible RNO-G solar flare polarization measurements, the Sun also represents an extremely useful above-surface calibration source.Using RNO-G data recorded during the summers of 2022 and 2023, we find signal excesses during solar flares reported by the solar-observing Callisto network and also in coincidence with ∼2/3 of the brightest excesses recorded by the SWAVES satellite. These observed flares are characterized by significant time-domain impulsivity. Using the known position of the Sun, the flare sample is used to calibrate the RNO-G absolute pointing on the radio signal arrival direction to sub-degree resolution. We thus establish the Sun as a regularly observed astronomical calibration source to provide the accurate absolute pointing required for neutrino astronomy.

One of the most promising techniques for detecting ultra-high energy neutrinos involves the use of radio antennas to observe the 10-1000 MHz radiation generated by the showers that neutrinos induce in large volumes of ice. The expected neutrino detection rates of one neutrino or less per detector station per 10 years make the characterization of backgrounds a priority. The largest natural background comes from ultra-high energy cosmic rays which are orders of magnitude more abundant than neutrinos. Particularly crucial is the understanding of geometries in which substantial energy of the cosmic-ray-induced air shower is deposited in the ice giving rise to a compact in-ice shower close to the ice surface. We calculated the radio emission of air-shower cores using the novel CORSIKA 8 code and found it to be similar to the predictions for neutrino-induced showers. For the first time, we calculated the detection rates for O(100m) deep antennas yielding 10-100 detections per year and detector station, which makes this a useful calibration source as these downward-going signals can be differentiated from neutrino-induced showers based on the signal arrival direction. However, the presence of reflection layers in the ice confuses the arrival directions, which makes this a potentially important background. We review the existing information on reflecting layers in the South Pole glacier and, for the first time, quantify the corresponding rate of reflected air-shower signals for the proposed IceCube-Gen2 radio array and discuss mitigation strategies. The reflectivity of the layers is the dominant uncertainty resulting in rate predictions of much less than one detection to several detections per year for IceCube-Gen2 if not mitigated.

We present the results of a search for 10–1000 GeV neutrinos from 2268 gamma-ray bursts (GRBs) over 8 yr of IceCube-DeepCore data. This work probes burst physics below the photosphere where electromagnetic radiation cannot escape. Neutrinos of tens of giga electronvolts are predicted in sub-photospheric collision of free-streaming neutrons with bulk-jet protons. In a first analysis, we searched for the most significant neutrino-GRB coincidence using six overlapping time windows centered on the prompt phase of each GRB. In a second analysis, we conducted a search for a group of GRBs, each individually too weak to be detectable, but potentially significant when combined. No evidence of neutrino emission is found for either analysis. The most significant neutrino coincidence is for Fermi-GBM GRB bn 140807500, with a p -value of 0.097 corrected for all trials. The binomial test used to search for a group of GRBs had a p -value of 0.65 after all trial corrections. The binomial test found a group consisting only of GRB bn 140807500 and no additional GRBs. The neutrino limits of this work complement those obtained by IceCube at tera electronvolt to peta electronvolt energies. We compare our findings for the large set of GRBs as well as GRB 221009A to the sub-photospheric neutron-proton collision model and find that GRB 221009A provides the most constraining limit on baryon loading. For a jet Lorentz factor of 300 (800), the baryon loading on GRB 221009A is lower than 3.85 (2.13) at a 90% confidence level.

Neutrino flares in the sky are searched for in data collected by IceCube between 2011 and 2021 May. This data set contains cascade-like events originating from charged-current electron neutrino and tau neutrino interactions and all-flavor neutral-current interactions. IceCube’s previous all-sky searches for neutrino flares used data sets consisting of track-like events originating from charged-current muon neutrino interactions. The cascade data set is statistically independent of the track data sets, and while inferior in angular resolution, the low-background nature makes it competitive and complementary to previous searches. No statistically significant flare of neutrino emission was observed in an all-sky scan. Upper limits are calculated on neutrino flares of varying duration from 1 hr to 100 days. Furthermore, constraints on the contribution of these flares to the diffuse astrophysical neutrino flux are presented, showing that multiple unresolved transient sources may contribute to the diffuse astrophysical neutrino flux.

The origin of high-energy galactic cosmic rays is yet to be understood, but some galactic cosmic-ray accelerators can accelerate cosmic rays up to PeV energies. The high-energy cosmic rays are expected to interact with the surrounding material or radiation, resulting in the production of gamma-rays and neutrinos. To optimize for the detection of such associated production of gamma-rays and neutrinos for a given source morphology and spectrum, a multimessenger analysis that combines gamma-rays and neutrinos is required. In this study, we use the Multi-Mission Maximum Likelihood framework with IceCube Maximum Likelihood Analysis software and HAWC Accelerated Likelihood to search for a correlation between 22 known gamma-ray sources from the third HAWC gamma-ray catalog and 14 yr of IceCube track-like data. No significant neutrino emission from the direction of the HAWC sources was found. We report the best-fit gamma-ray model and 90% CL neutrino flux limit from the 22 sources. From the neutrino flux limit, we conclude that, for five of the sources, the gamma-ray emission observed by HAWC cannot be produced purely from hadronic interactions. We report the limit for the fraction of gamma-rays produced by hadronic interactions for these five sources.

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.

We present the first three-parameter fit to a 3+1 sterile neutrino model using 7.634 years of data from the IceCube Neutrino Observatory on nu(mu) + (nu) over bar (mu) charged-current interactions in the energy range 500-9976 GeV. Our analysis is sensitive to the mass-squared splitting between the heaviest and lightest mass state (Delta m(41)(2)), the mixing matrix element connecting muon flavor to the fourth mass state (vertical bar U-mu 4 vertical bar(2)), and the element connecting tau flavor to the fourth mass state (vertical bar U-tau 4 vertical bar(2)). Predicted propagation effects in matter enhance the signature through a resonance as atmospheric neutrinos from the Northern Hemisphere traverse the Earth to the IceCube detector at the South Pole. The remaining sterile neutrino matrix elements are left fixed, with vertical bar U-e4 vertical bar(2) = 0 and delta(14) = 0, as they have a negligible effect, and delta(24) = pi is set to give the most conservative limits. The result is consistent with the no-sterile neutrino hypothesis with a probability of 4.3%. Profiling the likelihood of each parameter yields the 90% confidence levels: 2.4 eV(2) < Delta m(41)(2) < 9.6 eV(2), 0.0081 < vertical bar U-mu 4 vertical bar(2) < 0.10, and vertical bar U-tau 4 vertical bar(2) < 0.035, which narrows the allowed parameter-space for vertical bar U-tau 4 vertical bar(2). However, the primary result of this analysis is the first map of the 3+1 parameter space exploring the interdependence of Delta m(41)(2), vertical bar U-mu 4 vertical bar(2), and |U-tau 4 vertical bar(2).

This Letter presents the result of a 3+1 sterile neutrino search using 10.7 yr of IceCube data. We analyze atmospheric muon neutrinos that traverse the Earth with energies ranging from 0.5 to 100 TeV, incorporating significant improvements in modeling neutrino flux and detector response compared to earlier studies. Notably, for the first time, we categorize data into starting and throughgoing events, distinguishing neutrino interactions with vertices inside or outside the instrumented volume, to improve energy resolution. The best-fit point for a 3+1 model is found to be at sin^{2}(2θ_{24})=0.16 and Δm_{41}^{2}=3.5  eV^{2}, which agrees with previous iterations of this Letter. The result is consistent with the null hypothesis of no sterile neutrinos with a p value of 3.1%.

More than 10,000 photomultiplier tubes (PMTs) with a diameter of 80 mm will be installed in multi-PMT Digital Optical Modules (mDOMs) of the IceCube Upgrade. These have been tested and pre-calibrated at two sites. A throughput of more than 1000 PMTs per week with both sites was achieved with a modular design of the testing facilities and highly automated testing procedures. The testing facilities can easily be adapted to other PMTs, such that they can, e.g., be re-used for testing the PMTs for IceCube-Gen2. Single photoelectron response, high voltage dependence, time resolution, prepulse, late pulse, afterpulse probabilities, and dark rates were measured for each PMT. We describe the design of the testing facilities, the testing procedures, and the results of the acceptance tests.