Abstract:We present a strong hint of a connection between high energy $\gamma$-ray emitting blazars, very high energy neutrinos, and ultra high energy cosmic rays. We first identify potential hadronic sources by filtering $\gamma$-ray emitters %from existing catalogs that are in spatial coincidence with the high energy neutrinos detected by IceCube. The neutrino filtered $\gamma$-ray emitters are then correlated with the ultra high energy cosmic rays from the Pierre Auger Observatory and the Telescope Array by scanning in $\gamma$-ray flux ($F_{\gamma}$) and angular separation ($\theta$) between sources and cosmic rays. A maximal excess of 80 cosmic rays (42.5 expected) is found at $\theta\leq10^{\circ}$ from the neutrino filtered $\gamma$-ray emitters selected from the second hard {\it Fermi}-LAT catalogue (2FHL) and for $F_\gamma\left(>50\:\mathrm{GeV}\right)\geq1.8\times10^{-11}\:\mathrm{ph}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}$. The probability for this to happen is $2.4 \times 10^{-5}$, which translates to $\sim 2.4 \times 10^{-3}$ after compensation for all the considered trials. No excess of cosmic rays is instead observed for the complement sample of $\gamma$-ray emitters (i.e. not in spatial connection with IceCube neutrinos). A likelihood ratio test comparing the connection between the neutrino filtered and the complement source samples with the cosmic rays favours a connection between neutrino filtered emitters and cosmic rays with a probability of $\sim1.8\times10^{-3}$ ($2.9\sigma)$ after compensation for all the considered trials. The neutrino filtered $\gamma$-ray sources that make up the cosmic rays excess are blazars of the high synchrotron peak type. More statistics is needed to further investigate these sources as candidate cosmic ray and neutrino emitters.

A nearby FR I type radio galaxy 3C 120 as a neutrino emitting candidate

On the possible contributions of two nearby blazars to the NGC 4151 neutrino hotspot

Neutrino flares of radio blazars observed from TeV to PeV

GB6 J2113+1121: A multi-wavelength flaring gamma-ray blazar temporally and spatially coincident with the neutrino event IceCube-191001A

GB6 J2113+1121: A Multiwavelength Flaring Γ-Ray Blazar Temporally and Spatially Coincident with the Neutrino Event IceCube-191001A

Connecting blazars with ultra high energy cosmic rays and astrophysical neutrinos

Awakening of A Blazar at Redshift 2.7 Temporally Coincident with Arrival of Cospatial Neutrino Event IceCube-201221A

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

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

A two-zone blazar radiation model for "orphan" neutrino flares

Hunting for Neutrino Emission from Multifrequency Variable Sources

PKS 2254+074: A Blazar in Likely Association with the Neutrino Event IceCube-190619A

Search for Gev and X-Ray Flares Associated with the Icecube Track-Like Neutrinos

Star-forming galaxies as the origin of diffuse high-energy backgrounds: Gamma-ray and neutrino connections, and implications for starburst history

Blazar -- IceCube neutrino association revisited

Search for GeV gamma-ray flares associated with IceCube track-like neutrinos with Fermi Large Area Telescope

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

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

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

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

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