A unique and often overlooked property of a source loaded with relativistic protons is that it can become supercritical, i.e. it can undergo an abrupt transition from a radiatively inefficient to a radiatively efficient state once its proton energy density exceeds a certain threshold. In this paper, we investigate the temporal variability of hadronic systems in this hardly explored regime. We show that there exists a range of proton densities that prevent the system from reaching a steady state, but drive it instead in a quasi-periodic mode. The escaping radiation then exhibits limit cycles, even if all physical parameters are held constant in time. We extend our analysis to cases where the proton injection rate varies with time and explore the variability patterns of escaping radiation as the system moves in and out from the supercritical regime. We examine the relevance of our results to the variability of the prompt gamma-ray burst emission and show that, at least on a phenomenological level, some interesting analogies exist.

Detection of the IceCube-170922A neutrino coincident with the flaring blazar TXS 0506+056, the first and only ∼3σ high-energy neutrino source association to date, offers a potential breakthrough in our understanding of high-energy cosmic particles and blazar physics. We present a comprehensive analysis of TXS 0506+056 during its flaring state, using newly collected Swift, NuSTAR, and X-shooter data with Fermi observations and numerical models to constrain the blazar’s particle acceleration processes and multimessenger (electromagnetic (EM) and high-energy neutrino) emissions. Accounting properly for EM cascades in the emission region, we find a physically consistent picture only within a hybrid leptonic scenario, with γ-rays produced by external inverse-Compton processes and high-energy neutrinos via a radiatively subdominant hadronic component. We derive robust constraints on the blazar’s neutrino and cosmic-ray emissions and demonstrate that, because of cascade effects, the 0.1-100 keV emissions of TXS 0506+056 serve as a better probe of its hadronic acceleration and high-energy neutrino production processes than its GeV-TeV emissions. If the IceCube neutrino association holds, physical conditions in the TXS 0506+056 jet must be close to optimal for high-energy neutrino production, and are not favorable for ultrahigh-energy cosmic-ray acceleration. Alternatively, the challenges we identify in generating a significant rate of IceCube neutrino detections from TXS 0506+056 may disfavor single-zone models, in which γ-rays and high-energy neutrinos are produced in a single emission region. In concert with continued operations of the high-energy neutrino observatories, we advocate regular X-ray monitoring of TXS 0506+056 and other blazars in order to test single-zone blazar emission models, clarify the nature and extent of their hadronic acceleration processes, and carry out the most sensitive possible search for additional multimessenger sources.

Blazars are a sub-category of radio-loud active galactic nuclei with relativistic jets pointing towards to the observer. They are well-known for their non-thermal variable emission, which practically extends over the whole electromagnetic spectrum. Despite the plethora of multi-wavelength observations, the issue about the origin of the γ -ray and radio emission in blazar jets remains unsettled. Here, we construct a parametric leptonic model for studying the connection between the γ -ray and radio emission in both steady-state and flaring states of blazars. Assuming that relativistic electrons are injected continuously at a fixed distance from the black hole, we numerically study the evolution of their population as it propagates to larger distances while losing energy due to expansion and radiative cooling. In this framework, γ -ray photons are naturally produced at small distances (e.g., 10 - 3 pc) when the electrons are still very energetic, whereas the radio emission is produced at larger distances (e.g., 1 pc), after the electrons have cooled and the emitting region has become optically thin to synchrotron self-absorption due to expansion. We present preliminary results of our numerical investigation for the steady-state jet emission and the predicted time lags between γ -rays and radio during flares.