The IceCube report of a $\sim 3.5\sigma$ excess of $13\pm5$ neutrino events in the direction of the blazar TXS 05056+056 in 2014-2015 and the 2017 detection of a high-energy neutrino, 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 05056+056 using optical, X-ray, and gamma-ray data from the All-Sky Automated Survey for Supernovae (ASAS-SN), the Neil Gehrels Swift Observatory (Swift), the Monitor of All-sky X-ray Image (MAXI), and the Fermi Large Area Telescope (Fermi-LAT). 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 $\sim0.4-2$ on the muon neutrino number in ten 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.

Active galactic nuclei (AGN) with relativistic jets, powered by gas accretion onto their central supermassive black hole (SMBH), are unique laboratories for studying the physics of matter and elementary particles in extreme conditions that cannot be realized on Earth. For a long time since the discovery of AGN, photons were the only way to probe the underlying physical processes. The recent discovery of a very high energy neutrino, IceCube-170922A, coincident with a flaring blazar, TXS 0506+056, provides the first evidence that AGN jets are multi-messenger sources; they are capable of accelerating hadrons to very high energies, while producing non-thermal EM radiation and high-energy neutrinos. This new era of multi-messenger astronomy, which will mature in the next decade, offers us the unprecedented opportunity to combine more than one messenger to solve some long-standing puzzles of AGN jet physics: How do jets dissipate their energy to accelerate particles? What is the jet total kinetic power? Where and how do jets produce high-energy emission and neutrinos? What physical mechanisms drive the particle acceleration?

Detection of the IceCube-170922A neutrino coincident with the flaring blazar TXS 0506+056, the first and only 3-sigma 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 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 gamma-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-100keV 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 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 prime candidate sources for the high energy neutrinos recently detected by IceCube. Being one of the brightest sources in the extragalactic X-ray and $\gamma$-ray sky as well as one of the nearest blazars to Earth, Mrk 421 is an excellent source for testing the scenario of the blazar-neutrino connection. Here, we model the spectral energy distribution of Mrk 421 during a 13-day flare in 2010 with unprecedented multi-wavelength coverage, and calculate the respective neutrino flux. We find a correlation between the $>1$ PeV neutrino and photon fluxes, in all energy bands. Using typical IceCube through-going muon event samples with good angular resolution and high statistics, we derive the mean event rate above 100 TeV ($\sim0.57$ evt/yr) and show that it is comparable to that expected from a four-month quiescent period in 2009. Due to the short duration of the flare, an accumulation of similar flares over several years would be necessary to produce a meaningful signal for IceCube. To better assess this, we apply the correlation between the neutrino and $\gamma$-ray fluxes to the 6.9 yr Fermi-LAT light curve of Mrk 421. We find that the mean event count above 1 PeV for the full IceCube detector livetime is $3.59\pm0.60$ ($2.73\pm0.38$) $\nu_\mu+\bar{\nu}_\mu$ with (without) major flares included in our analysis. This estimate exceeds, within the uncertainties, the $95\%$ ($90\%$) threshold value for the detection of one or more muon (anti-)neutrinos. Meanwhile, the most conservative scenario, where no correlation of $\gamma$-rays and neutrinos is assumed, predicts $1.60\pm0.16$ $\nu_\mu+\bar{\nu}_\mu$ events. We conclude that a non-detection of high-energy neutrinos by IceCube would probe the neutrino/$\gamma$-ray flux correlation during major flares or/and the hadronic contribution to the blazar emission.

We present two ways of modeling the spectral energy distribution of blazars in the hadronic context and discuss the predictions of each "hadronic variant" on the spectral shape, the multi-wavelength variability, the cosmic-ray flux, and the high-energy neutrino emission. Focusing on the latter, we then present an application of the hadronic model to individual BL Lacs that were recently suggested to be the counterparts of some of the IceCube neutrinos, and conclude by discussing the contribution of the whole BL Lac class to the observed neutrino background.

We present a method of constraining the properties of the $\gamma$-ray emitting region in flat spectrum radio quasars (FSRQs) in the one-zone proton synchrotron model, where the $\gamma$-rays are produced by synchrotron radiation of relativistic protons. We show that for low enough values of the Doppler factor $\delta$, the emission from the electromagnetic (EM) cascade which is initiated by the internal absorption of high-energy photons from photohadronic interactions may exceed the observed $\sim$GeV flux. We use that effect to derive an absolute lower limit of $\delta$; first, an analytical one, in the asymptotic limit where the external radiation from the broad line region (BLR) is negligible, and then a numerical one in the more general case that includes BLR radiation. As its energy density in the emission region depends on $\delta$ and the region's distance from the galactic center, we use the EM cascade to determine a minimum distance for each value of $\delta$. We complement the EM cascade constraint with one derived from variability arguments and apply our method to the FSRQ 3C 273. We find that $\delta \gtrsim 18-20$ for $B \lesssim 30$ G and $\sim$day timescale variability; the emission region is located outside the BLR, namely at $r \gtrsim 10 R_{\rm BLR} \sim 3$ pc; the model requires at pc-scale distances stronger magnetic fields than those inferred from core shift observations; while the jet power exceeds by at least one order of magnitude the accretion power. In short, our results disfavour the proton synchrotron model for the FSRQ 3C 273.

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 model-predicted 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.

The recent discovery of extragalactic PeV neutrinos opens a new window to the exploration of cosmic-ray accelerators. The observed PeV neutrino flux is close to the Waxman-Bahcall upper bound implying that gamma-ray bursts (GRBs) may be the source of ultra-high energy cosmic rays (UHECRs). Starting with the assumption of the GRB-UHECR connection, we show using both analytical estimates and numerical simulations that the observed neutrinos can originate at the jet as a result of photopion interactions with the following implications: the neutrino spectra are predicted to have a cutoff at energy $\le 10$ PeV; the dissipation responsible for the GRB emission and cosmic-ray acceleration takes place at distances $r_{\rm diss} \approx 3 \times 10^{11}-3 \times 10^{13} \ {\rm cm}$ from the central engine; the Thomson optical depth at the dissipation region is $\tau_{\rm T} \sim 1$; the jet carries a substantial fraction of its energy in the form of Poynting flux at the dissipation region, and has a Lorentz factor $\Gamma \simeq 100-500$. The non-detection of $\sim$PeV neutrinos coincident with GRBs will indicate that GRBs are either poor cosmic accelerators or the dissipation takes place at small optical depths in the jet.

We explore a one-zone hadronic model that may be able to reproduce $\gamma$-ray burst (GRB) prompt emission with a minimum of free parameters. Assuming only that GRBs are efficient high-energy proton accelerators and without the presence of an ab initio photon field, we investigate the conditions under which the system becomes supercritical, i.e. there is a fast, non-linear transfer of energy from protons to secondary particles initiated by the spontaneous quenching of proton-produced $\gamma$-rays. We first show analytically that the transition to supercriticality occurs whenever the proton injection compactness exceeds a critical value, which favours high proton injection luminosities and a wide range of bulk Lorentz factors. The properties of supercriticality are then studied with a time-dependent numerical code that solves concurrently the coupled equations of proton, photon, electron, neutron and neutrino distributions. For conditions that drive the system deep into the supercriticality we find that the photon spectra obtain a Band-like shape due to Comptonization by cooled pairs and that the energy transfer efficiency from protons to $\gamma$-rays and neutrinos is high reaching $\sim 0.3$. Although some questions concerning its full adaptability to the GRB prompt emission remain open, supercriticality is found to be a promising process in that regard.

Aims: We investigate the role of the second synchrotron self-Compton (SSC) photon generation to the multiwavelength emission from the compact regions of sources that are characterized as misaligned blazars. For this, we focus on the nearest high-energy emitting radio galaxy Centaurus A and we revisit the one-zone SSC model for its core emission. Methods: We have calculated analytically the peak luminosities of the first and second SSC components by, first, deriving the steady-state electron distribution in the presence of synchrotron and SSC cooling and, then, by using appropriate expressions for the positions of the spectral peaks. We have also tested our analytical results against those derived from a numerical code where the full emissivities and cross-sections were used. Results: We show that the one-zone SSC model cannot account for the core emission of Centaurus A above a few GeV, where the peak of the second SSC component appears. We, thus, propose an alternative explanation for the origin of the high energy ($\gtrsim 0.4$ GeV) and TeV emission, where these are attributed to the radiation emitted by a relativistic proton component through photohadronic interactions with the photons produced by the primary leptonic component. We show that the required proton luminosities are not extremely high, e.g. $\sim 10^{43}$ erg/s, provided that the injection spectra are modelled by a power-law with a high value of the lower energy cutoff. Finally, we find that the contribution of the core emitting region of Cen A to the observed neutrino and ultra-high energy cosmic-ray fluxes is negligible.

We examine the neutrino and cosmic ray spectra resulting from two models of fitting the spectral energy distribution (SED) of the blazar Mrk~421 using a self-consistent leptohadronic code. The $\gamma$-ray emission is attributed to either synchrotron radiation of ultra-high energy protons (LHs model) or to synchrotron radiation from electrons that result from photopion interactions of lower energy protons (LH$\pi$ model). Although both models succeed in fitting satisfactorily the SED, the parameter values that they use result in significantly different neutrino and cosmic-ray spectra. For the LH$\pi$ model, which requires high proton energy density, we find that the neutrino spectrum peaks at an energy $E_{\rm \nu,peak}= 3.3$ PeV which falls well within the energy range of recent neutrino observations. While at the same time its peak flux is just under the sensitivity limit of IC-40 observations, it cannot produce ultra-high energy cosmic rays. In the LHs model, on the other hand, neutrinos are far from being detectable because of their low flux and peak energy at $E_{{\rm \nu,peak}} \simeq 100$ PeV. However, the propagation of protons produced by the decay of escaping neutrons results in an ultra-high energy cosmic ray flux close to that observed by Pierre Augere, HiRes and Telescope Array at energies $E_p \simeq 30$ EeV.

We investigate the origin of high-energy emission in blazars within the context of the leptohadronic one-zone model. We find that $\gamma$-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 X- and $\gamma$-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. 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 X-rays and $\gamma$-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 $\gamma$-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 $\gamma$-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 $\gamma-$rays. On the other hand, models based on pion-decay can reproduce this behaviour in a much more natural way.

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 present some results on the radiative signatures of the one zone hadronic model. For this we have solved five spatially averaged, time-dependent coupled kinetic equations which describe the evolution of relativistic protons, electrons, photons, neutrons and neutrinos in a spherical volume containing a magnetic field. Protons are injected and lose energy by synchrotron, photopair and photopion production. We model photopair and photopion using the results of relevant MC codes, like the SOPHIA code in the case of photopion, which give accurate description for the injection of secondaries which then become source functions in their respective equations. This approach allows us to calculate the expected photon and neutrino spectra simultaneously in addition to examining questions like the efficiency and the temporal behaviour of the hadronic models.

Recent observations by the H.E.S.S. collaboration of the Galactic Centre region have revealed what appears to be gamma-ray emission from the decay of pions produced by interactions of recently accelerated cosmic rays with local molecular hydrogen clouds. Synthesizing a 3-D hydrogen cloud map from the available data and assuming a diffusion coefficient of the form kappa(E) = kappa_0(E/E0)^delta, we performed Monte Carlo simulations of cosmic ray diffusion for various propagation times and values of kappa_0 and delta. By fitting the model gamma-ray spectra to the observed one we were able to infer the value of the diffusion coefficient in that environment (kappa = 3.0 +/- 0.2 kpc^2 Myr^-1 for E = 10^12.5 eV and for total propagation time 10^4 yr) as well as the source spectrum (2.1 < gamma < 2.3). Also, we found that proton losses can be substantial, which justifies our approach to the problem.

Diffuse VHE gamma radiation from the Galactic Centre ridge observed by the H.E.S.S. telescope has been convincingly linked with the propagation of recently accelerated cosmic rays that interact with molecular hydrogen clouds during their diffusion. Through a series of time-dependent simulations of that diffusion for different propagation parameters we have obtained the most probable values of the diffusion coefficient for the Galactic Centre region. Assuming that the diffusion coefficient is of the form kappa(E) = kappa_0*(E/E_0)^delta, then for different optimal combinations of kappa_0 and delta its value is obtained for cosmic rays originating from a central point (possibly Sgr A East) 10 kyr ago.

The non-linear back reaction of accelerated cosmic rays at the shock fronts, leads to the formation of a smooth precursor with a length scale corresponding to the diffusive scale of the energetic particles. Past works claimed that shocklets could be created in the precursor region of a specific shock width, which might energize few thermal particles to sufficient acceleration and furthermore this precursor region may act as confining large angle scatterer for very high energy cosmic rays. On the other hand, it has been shown that the smoothing of the shock front could lower the acceleration efficiency. These controversies motivated us to investigate numerically by Monte Carlo simulations the particle acceleration efficiency in oblique modified shocks. The results show flatter spectra compared to the spectra of the pressumed sharp discontinuity shock fronts. The findings are in accordance with theoretical predictions, since the scattering inside the precursor confines high energy particles to further scattering, resulting in higher energies making the whole acceleration process more efficient.