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 highenergy 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 and high-energy neutrino) emissions. Accounting properly for electromagnetic 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 highenergy 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 ultra-high-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.

The recent IceCube discovery of 0.1 − 1 PeV neutrinos of astrophysical origin opens up a new era for high-energy astrophysics. Although there are various astrophysical candidate sources, a firm association of the detected neutrinos with one (or more) of them is still lacking. A recent analysis of plausible astrophysical counterparts within the error circles of IceCube events showed that likely counterparts for nine of the IceCube neutrinos include mostly BL Lacs, among which Mrk 421. Motivated by this result and a previous independent analysis on the neutrino emission from Mrk 421, we test the BL Lac-neutrino connection in the context of a specific theoretical model for BL Lac emission. We model the spectral energy distribution (SED) of the BL Lacs selected as counterparts of the IceCube neutrinos using a one-zone leptohadronic model and mostly nearly simultaneous data. The neutrino flux for each BL Lac is self-consistently calculated, using photon and proton distributions specifically derived for every individual source. We find that the SEDs of the sample, although different in shape and flux, are all well fitted by the model using reasonable parameter values. Moreover, the modelpredicted neutrino flux and energy for these sources are of the same order of magnitude as those of the IceCube neutrinos. In two cases, namely Mrk 421 and H 1914-194, we find a suggestively good agreement between the model prediction and the detected neutrino flux. Our predictions for all the BL Lacs of the sample are in the range to be confirmed or disputed by IceCube in the next few years of data sampling.

Abstract:27 Since the discovery of the Van Allen radiation belts over 50 years ago, an explanation for Text:41

We present a time-dependent approach to the one-zone hadronic model in the case where the photon spectrum is produced by ultrarelativistic protons interacting with soft photons that are produced from protons and low magnetic fields. Assuming that protons are injected at a certain rate in a homogeneous spherical volume containing a magnetic field, the evolution of the system can be described by five coupled kinetic equations, for protons, electrons, photons, neutrons, and neutrinos. Photopair and photopion interactions are modelled using the results of Monte-Carlo simulations and, in particular, from the SOPHIA code for the latter. The coupling of energy losses and injection introduces a self-consistency in our approach and allows the study of the comparative relevancy of processes at various conditions, the efficiency of the conversion of proton luminosity to radiation, the resulting neutrino spectra, and the effects of time variability on proton injection, among other topics. We present some characteristic examples of the temporal behaviour of the system and show that this can be very different from the one exhibited by leptonic models. Furthermore, we argue that, contrary to the wide-held belief, there are parameter regimes where the hadronic models can become quite efficient. However, to keep the free parameters at a minimum and facilitate an in-depth study of the system, we have only concentrated on the case where protons are injected; i.e., we did not consider the effects of a co-accelerated leptonic component.

We investigate the origin of high-energy emission in blazars within the context of the leptohadronic one-zone model. We find that γ-ray emission can be attributed to synchrotron radiation either from protons or from secondary leptons produced via photohadronic processes. These possibilities imply differences not only in the spectral energy distribution (SED) but also in the variability signatures, especially in the Xand γ-ray regime. Thus, the temporal behavior of each leptohadronic scenario can be used to probe the particle population responsible for the high-energy emission as it can give extra information not available by spectral fits. In the present work we apply these ideas to the non-thermal emission of Mrk 421, which is one of the best monitored TeV blazars. We focus on the observations of March 2001, since during that period Mrk 421 showed multiple flares that have been observed in detail both in Xrays and γ-rays. First, we obtain pre-flaring fits to the SED using the different types of leptohadronic scenarios. Then, we introduce random-walk type, small-amplitude variations on the injection compactness or on the maximum energy of radiating particles and follow the subsequent response of the radiated photon spectrum. For each leptohadronic scenario, we calculate the X-ray and γ-ray fluxes and investigate their possible correlation. Whenever the 'input' variations lead, apart from flux variability, also to spectral variability, we present the resulting relations between the spectral index and the flux, both in X-rays and γ-rays. We find that proton synchrotron models are favoured energetically but require fine tuning between electron and proton parameters to reproduce the observed quadratic behaviour between X-rays and TeV γ−rays. On the other hand, models based on pion-decay can reproduce this behaviour in a much more natural way.