Abstract:A correlation has been reported between the arrival directions of high-energy IceCube events and gamma -ray blazars classified as intermediate- and high-synchrotron-peaked BL Lacs. Subsequent studies have investigated the optical properties of these sources, compiled and analyzed public multiwavelength data, and constrained their individual neutrino emission based on public IceCube point-source data. We provide a theoretical interpretation of public multiwavelength and neutrino point source data for the 32 BL Lac objects in the sample previously associated with an IceCube alert event. We combined the individual source results to draw conclusions regarding the multimesssenger properties of the sample and the required power in relativistic protons. We performed particle interaction modeling using open-source numerical simulation software. We constrained the model parameters using a novel and unique approach that simultaneously describes the host galaxy contribution, the observed synchrotron peak properties, the average multiwavelength fluxes, and, where possible, the IceCube point source constraints. We show that a single-zone leptohadronic model can describe the multiwavelength broadband fluxes from all 32 IceCube candidates. In some cases, the model suggests that hadronic emission may contribute a considerable fraction of the gamma -ray flux. The required power in relativistic protons ranges from a few percent to a factor of ten of the Eddington luminosity, which is energetically less demanding compared to other leptohadronic blazar models in recent literature. The model can describe the 68 confidence level IceCube flux for a large fraction of the masquerading BL Lacs in the sample, including TXS 0506+056; whereas, for true BL Lacs, the model predicts a low neutrino flux in the IceCube sensitivity range. Physically, this distinction is due to the presence of photons from broad line emission in masquerading BL Lacs, which increase the efficiency of hadronic interactions. The predicted neutrino flux peaks between a few petaelectronvolt and 100 PeV and scales positively with the flux in the gigaelectronvolt, megaelectronvolt, X-ray, and optical bands. Based on these results, we provide a list of the brightest neutrino emitters, which can be used for future searches targeting the 10-100 PeV regime.

Did IceCube discover Dark Matter around Blazars?

A Two-zone Model for Blazar Emission: Implications for TXS 0506+056 and the Neutrino Event IceCube-170922A

Hadronic processes at work in 5BZB J0630-2406

Leptohadronic Multimessenger Modeling of 324 Gamma-Ray Blazars

Recent IceCube Results: Diffuse Flux, Point Sources and Dark Matter

Blazar Flares as an Origin of High-Energy Cosmic Neutrinos?

New constraints on the dark matter-neutrino and dark matter-photon scattering cross sections from TXS 0506+056

Hadronuclear interactions in AGN jets as the origin of the diffuse high-energy neutrino background

AGN constraints on neutrino-dark matter scattering

The Spectra of IceCube Neutrino (SIN) candidate sources. V. Modeling and interpretation of multiwavelength and neutrino data

Blazar -- IceCube neutrino association revisited

Hadronuclear Interpretation of a High-Energy Neutrino Event Coincident with a Blazar Flare

Searching for temporary gamma-ray dark blazars associated with IceCube neutrinos

The Spectra of IceCube Neutrino (SIN) candidate sources -- V. Modeling and interpretation of multiwavelength and neutrino data

A Self-Consistent Leptonic–hadronic Interpretation of the Electromagnetic and Neutrino Emissions from Blazar TXS 0506+056

An NGC 1068-Informed Understanding of Neutrino Emission of the Active Galactic Nucleus TXS 0506+056

Search for correlations of high-energy neutrinos detected in IceCube with radio-bright AGN and gamma-ray emission from blazars

A Search for AGN sources of the IceCube Diffuse Neutrino Flux

Study of blazars detected by Fermi-LAT as high-energy neutrino sources

Possible correlation between unabsorbed hard X-rays and neutrinos in radio-loud and radio-quiet AGN