The IceCube report of a ∼ 3.5σ excess of 13 ± 5 neutrino events in the direction of the blazar TXS 0506+056 in 2014-2015 and the 2017 detection of a high-energy neutrino event, IceCube-170922A, during a gamma-ray flare from the same blazar, have revived the interest in scenarios for neutrino production in blazars. We perform comprehensive analyses on the long-term electromagnetic emission of TXS 0506+056 using optical, X-ray, and gamma-ray data from the All-Sky Automated Survey for Supernovae, the Neil Gehrels Swift Observatory, Monitor of All-sky X-ray Image, and the Fermi Large Area Telescope. We also perform numerical modeling of the spectral energy distributions (SEDs) in four epochs prior to 2017 with contemporaneous gamma-ray and lower-energy (optical and/or X-ray) data. We find that the multi-epoch SEDs are consistent with a hybrid leptonic scenario, where the gamma-rays are produced in the blazar zone via external inverse Compton scattering of accelerated electrons, and high-energy neutrinos are produced via the photomeson production process of co-accelerated protons. The multi-epoch SEDs can be satisfactorily explained with the same jet parameters and variable external photon density and electron luminosity. Using the maximal neutrino flux derived for each epoch, we put an upper limit of ∼0.4-2 on the muon neutrino number in 10 years of IceCube observations. Our results are consistent with the IceCube-170922A detection, which can be explained as an upper fluctuation from the average neutrino rate expected from the source, but in strong tension with the 2014-2015 neutrino flare.

3HSP J095507.9+355101 is an extreme blazar that has been possibly associated with a high-energy neutrino (IceCube-200107A) detected 1 day before the blazar was found to undergo a hard X-ray flare. We perform a comprehensive study of the predicted multimessenger emission from 3HSP J095507.9+355101 during its recent X-ray flare, but also in the long term. We focus on one-zone leptohadronic models, but we also explore alternative scenarios: (i) a blazar-core model, which considers neutrino production in the inner jet, close to the supermassive black hole; (ii) a hidden external-photon model, which considers neutrino production in the jet through interactions with photons from a weak broad line region; (iii) a proton-synchrotron model, where high-energy protons in the jet produce γ-rays via synchrotron; and (iv) an intergalactic cascade scenario, where neutrinos are produced in the intergalactic medium by interactions of a high-energy cosmic-ray beam escaping the jet. The Poisson probability to detect a single muon neutrino in 10 years from 3HSP J095507.9+355101 with the real-time IceCube alert analysis is ∼1% (3%) for the most optimistic one-zone leptohadronic model (the multi-zone blazar-core model). Meanwhile, detection of a single neutrino during the 44-day-long high X-ray flux-state period following the neutrino detection is 0.06%, according to our most optimistic leptohadronic model. The most promising scenarios for neutrino production also predict strong intrasource γ-ray attenuation above ∼100 GeV. If the association is real, then IceCube-Gen2 and other future detectors should be able to provide additional evidence for neutrino production in 3HSP J095507.9+355101 and other extreme blazars.

Hadronic supercriticalities are radiative instabilities that appear when large amounts of energy are stored in relativistic protons. When the proton energy density exceeds some critical value, a runaway process is initiated resulting in the explosive transfer of the proton energy into electron-positron pairs and radiation. The runaway also leads to an increase of the radiative efficiency, namely the ratio of the photon luminosity to the injected proton luminosity. We perform a comprehensive study of the parameter space by investigating the onset of hadronic supercriticalities for a wide range of source parameters (I.e. magnetic field strengths of 1 G-100 kG and radii of 1011-1016 cm) and maximum proton Lorentz factors (103-109). We show that supercriticalities are possible for the whole range of source parameters related to compact astrophysical sources, like gamma-ray bursts and cores and jets of active galactic nuclei. We also provide an in-depth look at the physical mechanisms of hadronic supercriticalities and show that magnetized relativistic plasmas are excellent examples of non-linear dynamical systems in high-energy astrophysics.