The powerful jets of blazars have been historically considered as likely sites of high-energy cosmic-ray acceleration. However, the particulars of the launched jet and the locations of leptonic and hadronic jet loading remain unclear. In the case when leptonic and hadronic particle injection occur jointly, a temporal correlation between synchrotron radiation and neutrino production is expected. We use a first catalog of millimeter wavelength (95–225 GHz) blazar light curves from the Atacama Cosmology Telescope for a time-dependent correlation with 12 yr of muon neutrino events from the IceCube South Pole Neutrino Observatory. Such millimeter emission traces activity of the bright jet base, which is often self-absorbed at lower frequencies and potentially gamma-ray opaque. We perform an analysis of the population, as well as analyses of individual, selected sources. We do not observe a significant signal from the stacked population. TXS 0506+056 is found as the most significant, individual source, though this detection is not globally significant in our analysis of selected active galactic nuclei. Our results suggest that the majority of millimeter-bright blazars are neutrino dim. In general, it is possible that many blazars have lighter, leptonic jets, or that only selected blazars provide exceptional conditions for neutrino production.

The IceCube Neutrino Observatory has observed extragalactic astrophysical neutrinos with an apparently isotropic distribution. Only a small fraction of the observed astrophysical neutrinos can be explained by known sources. Neutrino production is thought to occur in energetic environments that are ultimately powered by the gravitational collapse of dense regions of the large-scale mass distribution in the universe. Whatever their identity, neutrino sources likely trace this large-scale mass distribution. The clustering of neutrinos with a tracer of the large-scale structure may provide insight into the distribution of neutrino sources with respect to redshift and the identity of neutrino sources. We implement a two-point angular cross correlation of the Northern sky track events with an infrared galaxy catalog derived from the Wide-field Infrared Survey Explorer (WISE) and Two Micron All Sky Survey (2MASS) source catalogs, which trace the nearby large-scale structure. No statistically significant correlation is found between the neutrinos and this infrared galaxy catalog. We find that ≤54% of the diffuse muon neutrino flux can be attributed to sources correlated with the galaxy catalog with 90% confidence. Additionally, when assuming that the neutrino source comoving density evolves following a power law in redshift, dN_s/dV ∝ (1 + z)^k, we find that sources with negative evolution, in particular k < −1.75, are disfavored at the 90% confidence level.

The discovery of joint sources of high-energy neutrinos and gravitational waves has been a primary target for the LIGO, Virgo, KAGRA, and IceCube observatories. The joint detection of high-energy neutrinos and gravitational waves would provide insight into cosmic processes, from the dynamics of compact object mergers and stellar collapses to the mechanisms driving relativistic outflows. The joint detection of multiple cosmic messengers can also elevate the significance of the common observation even when some or all of the constituent messengers are subthreshold, i.e., not significant enough to declare their detection individually. Using data from the LIGO, Virgo, and IceCube observatories, including subthreshold events, we searched for common sources of gravitational waves and high-energy neutrinos during the third observing run of the Advanced LIGO and Advanced Virgo detectors. Our search did not identify significant joint sources. We derive constraints on the rate densities of joint sources. Our results constrain the isotropic neutrino emission from gravitational-wave sources for very high values of the total energy emitted in neutrinos (>1052–1054 erg).

We report improved measurements of the all flavor astrophysical neutrino spectrum with IceCube by combining complementary neutrino samples in two independent analyses. Both analyses show evidence of a harder spectrum at energies below ∼30 TeV compared to higher energies where the spectrum is well characterized by a power law. The spectrum is better described by a log parabola or a broken power law, the latter being the preferred model. Both, however, reject a single power law over an energy range 5 TeV-10 PeV with a significance >4σ, providing new constraints on properties of cosmic neutrino sources.

IceCube recently reported the observation of TeV neutrinos from the nearby Seyfert galaxy NGC 1068, and the corresponding neutrino flux is significantly higher than the upper limit implied by observations of GeV–TeV gamma rays. This suggests that neutrinos are produced near the supermassive black hole, where the radiation density is high enough to obscure gamma rays. We use a set of muon neutrinos with interaction vertices inside the detector, which have good sensitivity to sources in the southern sky, from IceCube data recorded between 2011 and 2021. We then search for individual and collective neutrino signals from 14 Seyfert galaxies in the southern sky selected from the Swift Burst Alert Telescope AGN Spectroscopic Survey. Using the correlations between keV X-rays and TeV neutrinos predicted by disk–corona models, and assuming production characteristics similar to NGC 1068, a collective neutrino signal search reveals an excess of 6.7 ^+4.0_-3.2 events, which is inconsistent with background expectations at the 3σ level of significance. In this Letter, we present new independent evidence that Seyfert galaxies contribute to the extragalactic flux of high-energy neutrinos.

We present GollumFit, a framework designed for performing binned-likelihood analyses on neutrino telescope data. GollumFit incorporates model parameters common to any neutrino telescope and also model parameters specific to the IceCube Neutrino Observatory. We provide a high-level overview of its key features and how the code is organized. We then discuss the performance of the fitting in a typical analysis scenario, highlighting the ability to fit over tens of nuisance parameters. We present some examples showing how to use the package for likelihood minimization tasks. This framework uniquely incorporates the particular model parameters necessary for neutrino telescopes, and solves an associated likelihood problem in a time-efficient manner.

The IceCube South Pole Neutrino Observatory has discovered the presence of a diffuse astrophysical neutrino flux at energies of TeV and beyond using neutrino induced muon tracks and cascade events from neutrino interactions. We present two analyses sensitive to neutrino events in the energy range 1 TeV to 10 PeV, using more than 10 years of IceCube data. Both analyses consistently reject a neutrino spectrum following a single power-law with a significance of >4σ in favor of a broken power law. We describe the methods implemented in the two analyses, the spectral constraints obtained, and the validation of the robustness of the results. Additionally, we report the detection of a muon neutrino in the medium energy starting events sample, or MESE, with an energy of 11.4+2.46-2.53 PeV, the highest energy neutrino observed by IceCube to date. The results presented here show insights into the spectral shape of astrophysical neutrinos, which has important implications for inferring their production processes in a multimessenger picture.

The advent of multimessenger astronomy has allowed for new types of source searches by neutrino detectors. We present the results of the search for 0.5–100 GeV astrophysical neutrinos detected with IceCube and emitted from compact binary mergers detected by the LIGO, Virgo, and KAGRA interferometers from their first run of observation (O1) to the end of the first part of the fourth (O4a). An innovative approach is used to lower the energy threshold to 0.5 GeV and to search for an excess of GeV neutrinos in time coincidence with astrophysical transient events. Furthermore, we use a statistical combination of all observations, a binomial test, to search for a subpopulation of neutrino emitters. No significant excess was found from the studied mergers, with a best post-trial 𝑝-value of 40%, and there is currently no hint of a population of GeV neutrino emitters found in the IceCube data (post-trial 𝑝−value=81%).

The search for sources of high-energy astrophysical neutrinos can be significantly advanced through a multimessenger approach, which seeks to detect the γ-rays that accompany neutrinos as they are produced at their sources. Multimessenger observations have so far provided the first evidence for a neutrino source, illustrated by the joint detection of the flaring blazar TXS 0506+056 in high-energy (E > 1 GeV) and very-high-energy (VHE; E > 100 GeV) γ-rays in coincidence with the high-energy neutrino IceCube-170922A, identified by IceCube. Imaging atmospheric Cherenkov telescopes (IACTs), namely FACT, H.E.S.S., MAGIC, and VERITAS, continue to conduct extensive neutrino target-of-opportunity follow-up programs. These programs have two components: follow-up observations of single astrophysical neutrino candidate events (such as IceCube-170922A), and observation of known γ-ray sources after the identification of a cluster of neutrino events by IceCube. Here we present a comprehensive analysis of follow-up observations of high-energy neutrino events observed by the four IACTs between 2017 September (after the IceCube-170922A event) and 2021 January. Our study found no associations between γ-ray sources and the observed neutrino events. We provide a detailed overview of each neutrino event and its potential counterparts. Furthermore, a joint analysis of all IACT data is included, yielding combined upper limits on the VHE γ-ray flux.

In the IceCube Neutrino Observatory, a signal of astrophysical neutrinos is obscured by backgrounds from atmospheric neutrinos and muons produced in cosmic-ray interactions. IceCube event selections used to isolate the astrophysical neutrino signal often focus on the morphology of the light patterns recorded by the detector. The analyses presented here use the new IceCube Enhanced Starting Track Event Selection (ESTES), which identifies events likely generated by muon–neutrino interactions within the detector geometry, focusing on neutrino energies of 1–500 TeV with a median angular resolution of 1.4. Selecting for starting-track events filters out not only the atmospheric-muon background but also the atmospheric-neutrino background in the southern sky. This improves IceCube’s muon–neutrino sensitivity to southern-sky neutrino sources, especially for Galactic sources that are not expected to produce a substantial flux of neutrinos above 100 TeV. In this work, the ESTES sample was applied for the first time to search for astrophysical sources of neutrinos, including a search for diffuse neutrino emission from the Galactic plane. No significant excesses were identified from any of the analyses; however, constraining limits are set on the hadronic emission from TeV gamma-ray Galactic plane objects and models of the diffuse Galactic plane neutrino flux.