Abstract:The origin of high-energy (HE) astrophysical neutrinos has remained an elusive hot topic in the field of HE astrophysics for the past decade. Apart from a handful of individual associations, the vast majority of HE neutrinos arise from unknown sources. While there are theoretically motivated candidate populations, such as blazars -- a subclass of active galactic nuclei with jets pointed toward our line of sight -- they have not been convincingly linked to HE neutrino production yet. Here, we perform a spatio-temporal association analysis between a sample of blazars (from the CGRaBS catalog) in the radio and optical bands and the most up-to-date IceCube HE neutrino catalog. We find that if the IceCube error regions are enlarged by 1$^ in quadrature, to account for unknown systematic errors at a maximal level, a spatio-temporal correlation between the multiwavelength light curves of the CGRaBS blazars and the IceCube HE neutrinos is hinted at, least at a 2.17sigma significance level. On the other hand, when the IceCube error regions are taken as their published values, we do not find any significant correlations. A discrepancy in the blazar-neutrino correlation strengths, when using such minimal and enlarged error region scenarios, was also obtained in a recent study by the IceCube collaboration. In our study, this difference arises because several flaring blazars -- coinciding with a neutrino arrival time -- happen to narrowly miss the published 90 -likelihood error region of the nearest neutrino event. For all of the associations driving our most significant correlations, the flaring blazar is much less than 1$^ away from the published error regions. Therefore, our results indicate that the question of the blazar-neutrino connection is highly sensitive to the reconstruction of the neutrino error regions, whose reliability is expected to improve with the next generation of neutrino observatories.

Anomalous cosmic-ray correlations revisited with a complete full-sky sample of BL Lac type objects

Anomalous Cosmic-ray Correlations Revisited with a Complete Full-sky Sample of BL Lac Type Objects

Testing the Correlation of Ultra-High Energy Cosmic Rays with High Redshift Sources

A Large-Area Cross-Correlation Study of High Galactic Latutude Soft and Hard X-ray Skies

Techniques for Searching for Anisotropy in the Arrival Directions of Ultra-High Energy Cosmic Rays

Association of the IceCube neutrinos with blazars in the CGRaBS sample

The correlations between BL Lacs and ultra-high-energy cosmic rays deflected by using different GMF models

Anisotropy and Corotation of Galactic Cosmic Rays

Are There Correlations in the HAWC and IceCube High Energy Skymaps Outside the Galactic Plane?

High-energy Neutrino Source Cross-correlations with Nearest Neighbor Distributions

Search for correlations of high-energy neutrinos and ultra-high energy cosmic rays

Honing cross-correlation tools for inference on ultra-high-energy cosmic-ray composition

An isotropic full-sky sample of optically selected blazars

A Cross-Correlation Study of High-energy Neutrinos and Tracers of Large-Scale Structure

Properties of Low Radio Luminosity RGB BL Lac Objects

Cosmological Implications of BL Lac Objects

Cosmic ray anisotropies to 5 PeV

Investigating the correlations between IceCube high-energy neutrinos and Fermi-LAT $γ$-ray observations. II

Anti-correlation Between the Extended Power and the Core Dominance Parameter of BL Lac Objects