THRILL: neuTrinos from HadRonIc fLares in bLazars
Collaborators: A. Mastichiadis (NKUA), P. Giommi (INAF), P. Padovani (ESO), T. Hovatta (University of Turku), E. Lindfors (University of Turku)
Blazars are a subclass of active galactic nuclei (AGN) with relativistic jets closely aligned to our line of sight  that are powered by accretion onto a central supermassive black hole . They are the most powerful persistent astrophysical sources of non-thermal electromagnetic radiation in the Universe, with spectral energy distributions spanning ~15 decades in energy, from radio frequencies up to high-energy γ-rays. Blazars are also characterized by flaring activity, namely fluctuations in flux away from an average value. Blazar flares are detected across the electromagnetic spectrum on a variety of timescales, ranging from years to minutes (for a review, see ).
In 2013, the IceCube neutrino telescope at the South Pole discovered a diffuse and isotropic flux of neutrinos of astrophysical origin [4,5]. The detection of a neutrino with energy >290 TeV, IceCube-170922A, in spatial and temporal coincidence with a γ-ray flare from blazar TXS 0506+056 in 2017 provided the first ~3σ neutrino source association at the time . A follow-up archival search at the position of the blazar revealed an excess of ~13 high-energy neutrinos with respect to the atmospheric background in 2014/15. This “neutrino flare” provided another ~3.5σ evidence for neutrino emission from the direction of TXS 0506+056 .
While the association of neutrino events and γ-ray flaring blazars has been getting plenty of attention, the connection has proven to be difficult to establish (see e.g. ). The most striking example is the lack of γ-ray flaring activity from TXS 0506+056 during the 2014/15 neutrino flare (e.g. [7, 9, 10]). From a theoretical perspective, this is not surprising. High-energy neutrinos are produced via inelastic collisions of relativistic protons with low-energy photons (optical or X-rays), which can also attenuate γ-rays. For sufficiently dense radiation fields, the neutrino flux may become high, but so does the opacity to γγ pair production.
This simple observation motivates the search for neutrino counterparts in other wavelengths. Recent studies have shown that AGN spatially associated with high-energy neutrinos typically have brighter parsec-scale radio emission compared to other sources in the sample . There is also evidence that radio flares (at GHz frequencies) observed close to the arrival time of a neutrino event from a blazar are unlikely to happen by chance .
Motivated by the recent progress in the field of multi-messenger astrophysics, we study alternative theoretical scenarios of neutrino production in flaring blazars. We explore the ramifications of our scenarios by comparing our model predictions against electromagnetic data obtained with satellites, like the Neil Gehrels Swift Observatory and the Fermi Gamma-Ray Telescope, and various radio/optical ground-based observatories, as well as neutrino data from the IceCube Neutrino Telescope.
We recently proposed a scenario according to which TeV-PeV neutrinos are produced in coincidence with X-ray flares that are powered by proton synchrotron radiation. In this case, neutrinos are produced by photomeson interactions of protons with their own synchrotron radiation, while MeV to GeV γ-rays are the result of synchrotron-dominated electromagnetic cascades developed in the source — see Figure 1. Using a time-dependent approach, we find that this “pure hadronic flaring” hypothesis has several interesting consequences. The X-ray flux is a good proxy for the all-flavor neutrino flux, while certain neutrino-rich X-ray flares may be dark in GeV-TeV γ-rays. Lastly, hadronic X-ray flares are accompanied by an equally bright MeV component that is detectable by proposed missions like e-ASTROGAM and AMEGO.
Figure 1: Animated figures showing the temporal evolution of the photon spectrum (grey lines) and the all-flavor neutrino spectrum (blue lines) for four indicative hadronic X-ray flares.
The research paper with our results can be found here: https://ui.adsabs.harvard.edu/abs/2020arXiv200912158M/abstract
 C. M. Urry and P. Padovani, 1995, PASP, 107  M. C. Begelman, R. D. Blandford, and M. J. Rees, 1984, Reviews of Modern Physics, 56.2  M. Boettcher, 2019, Galaxies, 7  M.G. Aartsen et al., 2013, Science, 342  M.G. Aartsen et al., 2013, Phys.Rev.Lett., 111  IceCube Collaboration, 2018a, Science, 361.6398  IceCube Collaboration, 2018b, Science, 361.6398  A. Franckowiak et al., 2020, ApJ, 893  A. Reimer et al., 2019, ApJ, 881  M. Petropoulou et al., 2020, ApJ 891  A. Plavin et al., 2020, ApJ 894  T. Hovatta et al., 2020, arXiv e-prints, arXiv:2009.10523