Plasmoids -- magnetized quasi-circular structures formed self-consistently in reconnecting current sheets -- were previously considered to be the graveyards of energetic particles. In this paper, we demonstrate the important role of plasmoids in shaping the particle energy spectrum in relativistic reconnection (i.e., with upstream magnetization $\sigma_{\rm up} \gg 1$). Using two dimensional particle-in-cell simulations in pair plasmas with $\sigma_{\rm up}=10$ and $100$, we study a secondary particle energization process that takes place inside compressing plasmoids. We demonstrate that plasmoids grow in time, while their interiors compress, amplifying the internal magnetic field. The magnetic field felt by particles injected in an isolated plasmoid increases linearly with time, which leads to particle energization as a result of magnetic moment conservation. For particles injected with a power-law distribution function, this energization process acts in such a way that the shape of the injected power law is conserved, while producing an additional non-thermal tail $f(E)\propto E^{-3}$ at higher energies followed by an exponential cutoff. The cutoff energy, which increases with time as $E_{\rm cut}\propto\sqrt{t}$, can greatly exceed $\sigma_{\rm up} m_e c^2$. We analytically predict the secondary acceleration timescale and the shape of the emerging particle energy spectrum, which can be of major importance in certain astrophysical systems, such as blazar jets.

The afterglow emission from gamma-ray bursts (GRBs) is believed to originate from a relativistic blast wave driven into the circumburst medium. Although the afterglow emission from radio up to X-ray frequencies is thought to originate from synchrotron radiation emitted by relativistic, non-thermal electrons accelerated by the blast wave, the origin of the emission at high energies (HE; $\gtrsim$~GeV) remains uncertain. The recent detection of sub-TeV emission from GRB~190114C by MAGIC raises further debate on what powers the very high-energy (VHE; $\gtrsim 300$GeV) emission. Here, we explore the inverse Compton scenario as a candidate for the HE and VHE emissions, considering two sources of seed photons for scattering: synchrotron photons from the blast wave (synchrotron self-Compton or SSC) and isotropic photon fields external to the blast wave (external Compton). For each case, we compute the multi-wavelength afterglow spectra and light curves. We find that SSC will dominate particle cooling and the GeV emission, unless a dense ambient infrared photon field, typical of star-forming regions, is present. Additionally, considering the extragalactic background light attenuation, we discuss the detectability of VHE afterglows by existing and future gamma-ray instruments for a wide range of model parameters. Studying GRB~190114C, we find that its afterglow emission in the \fermi-LAT band is synchrotron-dominated.The late-time \fermi-LAT measurement (i.e., $t\sim 10^4$~s), and the MAGIC observation also set an upper limit on the energy density of a putative external infrared photon field (i.e. $\lesssim 3\times 10^{-9}\,{\rm erg\,cm^{-3}}$), making the inverse Compton dominant in the sub-TeV energies.

The central engine in long gamma-ray bursts (GRBs) is thought to be a compact object produced by the core collapse of massive stars, but its exact nature (black hole or millisecond magnetar) is still debatable. Although the central engine of GRB collapsars is hidden to direct observation, its properties may be imprinted on the accompanying electromagnetic signals. We aim to decipher the generic properties of central engines that are consistent with prompt observations of long GRBs detected by the Burst Alert Telescope (BAT) on board the Neil Gehrels Swift Observatory. Adopting a generic model for the central engine, in which the engine power and activity timescale are independent of each other, we perform Monte Carlo simulations of long GRBs produced by jets that successfully breakout from the star. Our simulations consider the dependence of the jet breakout timescale on the engine luminosity and the effects of the detector's flux threshold. The two-dimensional (2D) distribution of simulated detectable bursts in the gamma-ray luminosity versus gamma-ray duration plane is consistent with the observed one for a range of parameter values describing the central engine. The intrinsic 2D distribution of simulated collapsar GRBs peaks at lower gamma-ray luminosities and longer durations than the observed one, a prediction that can be tested in the future with more sensitive detectors. Black-hole accretors, whose power and activity time are set by the large-scale magnetic flux through the progenitor star and stellar structure, respectively, are compatible with the properties of the central engine inferred by our model.

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 $10^{11}$ cm$-10^{16}$ cm) and maximum proton Lorentz factors ($10^3-10^9$). We show that supercriticalities are possible for the whole range of source parameters related to compact astrophysical sources, like gamma-rays bursts, 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.

Tidal disruption events (TDE) have been considered as cosmic-ray and neutrino sources for a decade. We suggest two classes of new scenarios for high-energy multi-messenger emission from TDEs that do not have to harbor powerful jets. First, we investigate high-energy neutrino and gamma-ray production in the core region of a supermassive black hole. In particular, we show that about 1-100 TeV neutrinos and MeV gamma-rays can efficiently be produced in hot coronae around an accretion disk. We also study the consequences of particle acceleration in radiatively inefficient accretion flows (RIAFs). Second, we consider possible cosmic-ray acceleration by disk-driven winds or interactions between tidal streams, and show that subsequent hadronuclear and photohadronic interactions inside the TDE debris lead to GeV-PeV neutrinos and sub-GeV cascade gamma-rays. We demonstrate that these models should be accompanied by soft gamma-rays or hard X-rays, which can be used for future observational tests. Although this work aims to present models of non-jetted high-energy emission, we discuss the implications of the TDE AT2019dsg that might coincide with the high-energy neutrino IceCube-191001A, by considering the corona, RIAF, hidden sub-relativistic wind, and hidden jet models. It is not yet possible to be conclusive about their physical association and the expected number of neutrinos is typically much less than unity. We find that the most optimistic cases of the corona and hidden wind models could be consistent with the observation of IceCube-191001A, whereas jet models are unlikely to explain the multi-messenger observations.

The joint detection of GW~170817/GRB 170817 confirmed the long-standing theory that binary neutron star mergers produce short gamma-ray burst (sGRB) jets that can successfully break out of the surrounding ejecta. At the same time, the association with a kilonova provided unprecedented information regarding the physical properties (such as masses and velocities) of the different ejecta constituents. Combining this knowledge with the observed luminosities and durations of cosmological sGRBs detected by the Burst Alert Telescope (BAT) onboard the Neil Gehrels Swift Observatory, we revisit the breakout conditions of sGRB jets. Assuming self-collimation of sGRB jets does not play a critical role, we find that the time interval between the binary merger and the launching of a typical sGRB jet is $\lesssim0.1$~s. We also show that for a fraction of at least $\sim 30\%$ of sGRBs, the usually adopted assumption of static ejecta is inconsistent with observations, even if the polar ejecta mass is an order of magnitude smaller than the one in GRB 170817. Our results disfavour magnetar central engines for powering cosmological sGRBs, limit the amount of energy deposited in the cocoon prior to breakout, and suggest that the observed delay of $\sim 1.$7~s in GW 170817 /GRB 170817 between the gravitational wave and $\gamma$-ray signals is likely dominated by the propagation time of the jet to the $\gamma$-ray production site.

3HSP J095507.9+355101 is an extreme blazar to have been possibly associated with a high-energy neutrino, IceCube-200107A, which was detected one day before the blazar was observed and found to undergo a hard X-ray flare. Motivated by this observation, 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 single-zone leptohadronic models, where the broadband photon and high-energy neutrino emissions originate from same region of the jet, but we also explore alternative scenarios: (i) a blazar-core model, which considers neutrino production in the inner jet close to the accreting supermassive black hole, (ii) a hidden external-photon model, which considers neutrino production in the jet through interactions with photons from a possible weak broad line region, (iii) a one-zone proton synchrotron model, where high-energy protons produce $\gamma$-rays in the jet 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 blazar jet. We find that the Poisson probability to detect one or more neutrinos in 8 years from 3HSP J095507.9+355101 is $\sim 2\%$ ($3\%$) when considering the most optimistic one-zone leptohadronic model (the multi-zone blazar-core model), while detection of one neutrino event during the X-ray flare is much less likely. 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, assuming similar physical conditions.

The origin of high-energy emission in blazars jets (i.e., leptonic versus hadronic) has been a long-standing matter of debate. Here, we focus on one variant of hadronic models where proton synchrotron radiation accounts for the observed steady $\gamma$-ray blazar emission. Using analytical methods, we derive the minimum jet power ($P_{j,\min}$) for the largest blazar sample analyzed to date (145 sources), taking into account uncertainties of observables and jet's physical parameters. We compare $P_{j,\min}$ against three characteristic energy estimators for accreting systems, i.e., the Eddington luminosity, the accretion disk luminosity, and the power of the Blandford-Znajek process, and find that $P_{j,\min}$ is about 2 orders of magnitude higher than all energetic estimators for the majority of our sample. The derived magnetic field strengths in the emission region require either large amplification of the jet's magnetic field (factor of 30) or place the $\gamma$-ray production site at sub-pc scales. The expected neutrino emission peaks at $\sim 0.1-10$ EeV, with typical peak neutrino fluxes $\sim 10^{-4}$ times lower than the peak $\gamma$-ray fluxes. We conclude that if relativistic hadrons are present in blazar jets, they can only produce a radiatively subdominant component of the overall spectral energy distribution of the blazar's steady emission.

This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process.

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.

We report on a multi-band variability and correlation study of the TeV blazar Mrk 421 during an exceptional flaring activity observed from 2013 April 11 to 2013 April 19. The study uses, among others, data from GASP-WEBT, Swift, NuSTAR, Fermi-LAT, VERITAS, and MAGIC. The large blazar activity, and the 43 hours of simultaneous NuSTAR and MAGIC/VERITAS observations, permitted variability studies on 15 minute time bins, and over three X-ray bands (3-7 keV, 7-30 keV and 30-80 keV) and three very-high-energy (>0.1 TeV, hereafter VHE) gamma-ray bands (0.2-0.4 TeV, 0.4-0.8 TeV and >0.8 TeV). We detected substantial flux variations on multi-hour and sub-hour timescales in all the X-ray and VHE gamma-ray bands. The characteristics of the sub-hour flux variations are essentially energy-independent, while the multi-hour flux variations can have a strong dependence on the energy of the X-ray and the VHE gamma rays. The three VHE bands and the three X-ray bands are positively correlated with no time-lag, but the strength and the characteristics of the correlation changes substantially over time and across energy bands. Our findings favour multi-zone scenarios for explaining the achromatic/chromatic variability of the fast/slow components of the light curves, as well as the changes in the flux-flux correlation on day-long timescales. We interpret these results within a magnetic reconnection scenario, where the multi-hour flux variations are dominated by the combined emission from various plasmoids of different sizes and velocities, while the sub-hour flux variations are dominated by the emission from a single small plasmoid moving across the magnetic reconnection layer.

The IceCube collaboration reported a $\sim 3.5\sigma$ excess of $13\pm5$ neutrino events in the direction of the blazar TXS 0506+56 during a $\sim$6 month period in 2014-2015, as well as the ($\sim3\sigma$) detection of a high-energy muon neutrino during an electromagnetic flare in 2017. We explore the possibility that the 2014-2015 neutrino excess and the 2017 multi-messenger flare are both explained in a common physical framework that relies on the emergence of a relativistic neutral beam in the blazar jet due to interactions of accelerated cosmic rays (CRs) with photons. We demonstrate that the neutral beam model provides an explanation for the 2014-2015 neutrino excess without violating X-ray and $\gamma$-ray constraints, and also yields results consistent with the detection of one high-energy neutrino during the 2017 flare. If both neutrino associations with TXS 05065+056 are real, our model requires that (i) the composition of accelerated CRs is light, with a ratio of helium nuclei to protons $\gtrsim5$, (ii) a luminous external photon field ($\sim 10^{46}$ erg s$^{-1}$) variable (on year-long timescales) is present, and (iii) the CR injection luminosity as well as the properties of the dissipation region (i.e., Lorentz factor, magnetic field, and size) vary on year-long timescales.

The All-sky Medium Energy Gamma-ray Observatory (AMEGO) is a probe class mission concept that will provide essential contributions to multimessenger astrophysics in the late 2020s and beyond. AMEGO combines high sensitivity in the 200 keV to 10 GeV energy range with a wide field of view, good spectral resolution, and polarization sensitivity. Therefore, AMEGO is key in the study of multimessenger astrophysical objects that have unique signatures in the gamma-ray regime, such as neutron star mergers, supernovae, and flaring active galactic nuclei. The order-of-magnitude improvement compared to previous MeV missions also enables discoveries of a wide range of phenomena whose energy output peaks in the relatively unexplored medium-energy gamma-ray band.

Blazar emission models based on magnetic reconnection succeed in reproducing many observed spectral and temporal features, including the short-duration luminous flaring events. Plasmoids, a self-consistent by-product of the tearing instability in the reconnection layer, can be the main source of blazar emission. Kinetic simulations of relativistic reconnection have demonstrated that plasmoids are characterized by rough energy equipartition between their radiating particles and magnetic fields. This is the main reason behind the apparent shortcoming of plasmoid-dominated emission models to explain the observed Compton ratios of BL Lac objects. Here, we demonstrate that the radiative interactions among plasmoids, which have been neglected so far, can assist in alleviating this contradiction. We show that photons emitted by large, slow-moving plasmoids can be a potentially important source of soft photons to be then up-scattered, via inverse Compton, by small fast-moving, neighboring plasmoids. This inter-plasmoid Compton scattering process can naturally occur throughout the reconnection layer, imprinting itself as an increase in the observed Compton ratios from those short and luminous plasmoid-powered flares within BL Lac sources, while maintaining energy equipartition between radiating particles and magnetic fields.

A high-energy muon neutrino event, IceCube-170922A, was recently discovered in both spatial and temporal coincidence with a gamma-ray flare of the blazar TXS 0506+056. It has been shown, with standard one-zone models, that neutrinos can be produced in the blazar jet via hadronic interactions, but with a flux which is mostly limited by the X-ray data. In this work, we explore the neutrino production from TXS 0506+056 by invoking two physically distinct emission zones in the jet, separated by the broad line region (BLR). Using the Doppler-boosted radiation of the BLR as the target photon field, the inner zone accounts for the neutrino and gamma-ray emission via $p\gamma$ interactions and inverse Compton scattering respectively, while the outer zone produces the optical and X-ray emission via synchrotron and synchrotron self-Compton processes. The different conditions of the two zones allow us to suppress the X-ray emission from the electromagnetic cascade, and set a much higher upper limit on the muon neutrino flux (i.e., $\sim 10^{-11}\rm erg~cm^{-2}s^{-1}$) than in one-zone models. We compare, in detail, our scenario with one-zone models discussed in the literature, and argue that differentiating between such scenarios will become possible with next generation neutrino telescopes, such as IceCube-Gen2.

Inverse Compton-pair cascades are initiated when gamma-rays are absorbed on an ambient soft photon field to produce relativistic pairs, which in turn up-scatter the same soft photons to produce more gamma-rays. If the Compton scatterings take place in the deep Klein-Nishina regime, then triplet pair production ($e\gamma_b \rightarrow ee^{+}e^{-}$) becomes relevant and may even regulate the development of the cascade. We investigate the properties of pair-Compton cascades with triplet pair production in accelerating gaps, i.e., regions with an unscreened electric field. Using the method of transport equations for the particle evolution, we compute the growth rate of the pair cascade as a function of the accelerating electric field in the presence of black-body and power-law ambient photon fields. Informed by the numerical results, we derive simple analytical expressions for the peak growth rate and the corresponding electric field. We show that for certain parameters, which can be realized in the vicinity of accreting supermassive black holes at the centers of active galactic nuclei, the pair cascade may well be regulated by inverse Compton scattering in the deep Klein-Nishina regime and triplet pair production. We present indicative examples of the escaping gamma-ray radiation from the gap, and discuss our results in application to the TeV observations of radio galaxy M87.

NGC300 ULX1 is an ultra-luminous X-ray pulsar, showing an unprecedented spin evolution, from about 126 s to less than 20 s in only 4 years, consistent with steady mass accretion rate. Following its discovery we have been monitoring the system with Swift/XRT and NICER to further study its properties. We found that even though the observed flux of the system dropped by a factor of $\gtrsim$20, the spin-up rate remained almost constant. A possible explanation is that the decrease in the observed flux is a result of increased absorption of obscuring material due to outflows or a precessing accretion disk.

This is a science white paper submitted to the Astro-2020 and Plasma-2020 Decadal Surveys. The paper describes the present status and emerging opportunities in Extreme Plasma Astrophysics -- a study of astrophysically-relevant plasma processes taking place under extreme conditions that necessitate taking into account relativistic, radiation, and QED effects.

Magnetic reconnection is often invoked to explain the non-thermal radiation of relativistic outflows, including jets of active galactic nuclei (AGN). Motivated by the largely unknown plasma composition of AGN jets, we study reconnection in the unexplored regime of electron-positron-proton (pair-proton) plasmas with large-scale two-dimensional particle-in-cell simulations. We cover a wide range of pair multiplicities (lepton-to-proton number ratio $\kappa=1-199$) for different values of the all-species plasma magnetization ($\sigma=1,3$ and 10) and electron temperature ($\Theta_e\equiv kT_e/m_ec^2=0.1-100$). We focus on the dependence of the post-reconnection energy partition and lepton energy spectra on the hot pair plasma magnetization $\sigma_{e,h}$ (i.e., the ratio of magnetic to pair enthalpy densities). We find that the post-reconnection energy is shared roughly equally between magnetic fields, pairs, and protons for $\sigma_{e,h}\gtrsim 3$. We empirically find that the mean lepton Lorentz factor in the post-reconnection region depends on $\sigma, \Theta_e$, and $\sigma_{e,h}$ as $\langle \gamma_e-1\rangle \approx \sqrt{\sigma}(1+4\Theta_e)\left(1+\sigma_{e,h}/30\right)$, for $\sigma\ge1$. The high-energy part of the post-reconnection lepton energy distributions can be described by a power law, whose slope is mainly controlled by $\sigma_{e,h}$ for $\kappa \gtrsim 3-6$, with harder power laws obtained for higher magnetizations. We finally show that reconnection in pair-proton plasmas with multiplicities $\kappa \sim 1-20$, magnetizations $\sigma \sim 1-10$, and temperatures $\Theta_e \sim 1-10$ results in particle power law slopes and average electron Lorentz factors that are consistent with those inferred in leptonic models of AGN jet emission.

MeV blazars are the most luminous persistent sources in the Universe and emit most of their energy in the MeV band. These objects display very large jet powers and accretion luminosities and are known to host black holes with a mass often exceeding $10^9 M_{\odot}$. An MeV survey, performed by a new generation MeV telescope which will bridge the entire energy and sensitivity gap between the current generation of hard X-ray and gamma-ray instruments, will detect $>$1000 MeV blazars up to a redshift of $z=5-6$. Here we show that this would allow us: 1) to probe the formation and growth mechanisms of supermassive black holes at high redshifts, 2) to pinpoint the location of the emission region in powerful blazars, 3) to determine how accretion and black hole spin interplay to power the jet.

The possible association of the blazar TXS 0506+056 with a high-energy neutrino detected by IceCube holds the tantalizing potential to answer three astrophysical questions: 1. Where do high-energy neutrinos originate? 2. Where are cosmic rays produced and accelerated? 3. What radiation mechanisms produce the high-energy {\gamma}-rays in blazars? The MeV gamma-ray band holds the key to these questions, because it is an excellent proxy for photo-hadronic processes in blazar jets, which also produce neutrino counterparts. Variability in MeV gamma-rays sheds light on the physical conditions and mechanisms that take place in the particle acceleration sites in blazar jets. In addition, hadronic blazar models also predict a high level of polarization fraction in the MeV band, which can unambiguously distinguish the radiation mechanism. Future MeV missions with a large field of view, high sensitivity, and polarization capabilities will play a central role in multi-messenger astronomy, since pointed, high-resolution telescopes will follow neutrino alerts only when triggered by an all-sky instrument.

The constantly improving sensitivity of ground-based and space-borne observatories has made possible the detection of high-energy emission (X-rays and gamma-rays) from several thousands of extragalactic sources. Enormous progress has been made in measuring the continuum flux enabling us to perform imaging, spectral and timing studies. An important remaining challenge for high-energy astronomy is measuring polarization. The capability to measure polarization is being realized currently at X-ray energies (e.g. with IXPE), and sensitive gamma-ray telescopes capable of measuring polarization, such as AMEGO, AdEPT, e-ASTROGAM, etc., are being developed. These future gamma-ray telescopes will probe the radiation mechanisms and magnetic fields of relativistic jets from active galactic nuclei at spatial scales much smaller than the angular resolution achieved with continuum observations of the instrument. In this white paper, we discuss the scientific potentials of high-energy polarimetry, especially gamma-ray polarimetry, including the theoretical implications, and observational technology advances being made. In particular, we will explore the primary scientific opportunities and wealth of information expected from synergy of multi-wavelength polarimetry that will be brought to multi-messenger astronomy.

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?

The high-energy universe has revealed that energetic particles are ubiquitous in the cosmos and play a vital role in the cultivation of cosmic environments on all scales. Though they play a key role in cultivating the cosmological environment and/or enabling our studies of it, there is still much we do not know about AGNs and GRBs, particularly the avenue in which and through which they supply radiation and energetic particles, namely their jets. This White Paper is the second of a two-part series highlighting the most well-known high-energy cosmic accelerators and contributions that MeV gamma-ray astronomy will bring to understanding their energetic particle phenomena. The focus of this white paper is active galactic nuclei and gamma-ray bursts.

This Astro2020 white paper advocates for a multi-messenger approach that combines high-energy neutrino and broad multi-wavelength electromagnetic observations to study AGN during the coming decade. The unique capabilities of these joint observations promise to solve several long-standing issues in our understanding of AGN as powerful cosmic accelerators.

We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science.

Motivated by the observation of a $>290$ TeV muon neutrino by IceCube, coincident with a $\sim$6 month-long $\gamma$-ray flare of the blazar TXS 0506+056, and an archival search which revealed $13 \pm 5$ further, lower-energy neutrinos in the direction of the source in 2014-2015, we discuss the likely contribution of blazars to the diffuse high-energy neutrino intensity, the implications for neutrino emission from TXS 0506+056 based on multi-wavelength observations of the source, and a multi-zone model that allows for sufficient neutrino emission so as to reconcile the multi-wavelength cascade constraints with the neutrino emission seen by IceCube in the direction of TXS 0506+056.

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 $\gamma$-ray and radio emission in blazar jets remains unsettled. Here, we construct a parametric leptonic model for studying the connection between the $\gamma$-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, $\gamma$-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 $\gamma$-rays and radio during flares.

We present evidence that TXS 0506+056, the first plausible non-stellar neutrino source, despite appearances, is not a blazar of the BL Lac type but is instead a masquerading BL Lac, i.e., intrinsically a flat-spectrum radio quasar with hidden broad lines and a standard accretion disk. This re-classification is based on: (1) its radio and O II luminosities; (2) its emission line ratios; (3) its Eddington ratio. We also point out that the synchrotron peak frequency of TXS 0506+056 is more than two orders of magnitude larger than expected by the so-called "blazar sequence", a scenario which has been assumed by some theoretical models predicting neutrino (and cosmic-ray) emission from blazars. Finally, we comment on the theoretical implications this re-classification has on the location of the $\gamma$-ray emitting region and our understanding of neutrino emission in blazars.

The joint detection of gravitational waves (GWs) and $\gamma$-rays from a binary neutron-star (NS) merger provided a unique view of off-axis GRBs and an independent measurement of the NS merger rate. Comparing the observations of GRB170817 with those of the regular population of short GRBs (sGRBs), we show that an order unity fraction of NS mergers result in sGRB jets that breakout of the surrounding ejecta. We argue that the luminosity function of sGRBs, peaking at $\approx 2\times 10^{52}\, \mbox{erg s}^{-1}$, is likely an intrinsic property of the sGRB central engine and that sGRB jets are typically narrow with opening angles $\theta_0 \approx 0.1$. We perform Monte Carlo simulations to examine models for the structure and efficiency of the prompt emission in off axis sGRBs. We find that only a small fraction ($\sim 0.01-0.1$) of NS mergers detectable by LIGO/VIRGO in GWs is expected to be also detected in prompt $\gamma$-rays and that GW170817-like events are very rare. For a NS merger rate of $\sim 1500$ Gpc$^{-3}$ yr$^{-1}$, as inferred from GW170817, we expect within the next decade up to $\sim 12$ joint detections with off-axis GRBs for structured jet models and just $\sim 1$ for quasi-spherical cocoon models where $\gamma$-rays are the result of shock breakout. Given several joint detections and the rates of their discoveries, the different structure models can be distinguished. In addition the existence of a cocoon with a reservoir of thermal energy may be observed directly in the UV, given a sufficiently rapid localisation of the GW source.

The multi-wavelength spectral and temporal variability observed in blazars set tight constraints on current theoretical emission models. Here, we investigate the relativistic magnetic reconnection process as a source of blazar emission in which quasi-spherical plasmoids, containing relativistic particles and magnetic fields, are associated with the emission sites in blazar jets. By coupling recent two-dimensional particle-in-cell simulations of relativistic reconnection with a time-dependent radiative transfer code, we compute the non-thermal emission from a chain of plasmoids formed during a reconnection event. The derived photon spectra display characteristic features observed in both BL Lac sources and flat spectrum radio quasars, with the distinction made by varying the strength of the external photon fields, the jet magnetization, and the number of pairs per proton contained within. Light curves produced from reconnection events are composed of many fast and powerful flares that appear on excess of a slower evolving envelope produced by the cumulative emission of medium-sized plasmoids. The observed variability is highly dependent upon the orientation of the reconnection layer with respect to the blazar jet axis and to the observer. Our model provides a physically motivated framework for explaining the multi-timescale blazar variability across the entire electromagnetic spectrum.

We consider implications of high-energy neutrino emission from blazar flares, including the recent event IceCube-170922A and the 2014-2015 neutrino flare that could originate from TXS 0506+056. First, we discuss their contribution to the diffuse neutrino intensity taking into account various observational constraints. Blazars are likely to be subdominant in the diffuse neutrino intensity at sub-PeV energies, and we show that blazar flares like those of TXS 0506+056 could make <1-10 percent of the total neutrino intensity. We also argue that the neutrino output of blazars can be dominated by the flares in the standard leptonic scenario for their gamma-ray emission, and energetic flares may still be detected with a rate of <1 per year. Second, we consider multi-messenger constraints on the source modeling. We show that luminous neutrino flares should be accompanied by luminous broadband cascade emission, emerging also in X-rays and gamma-rays. This implies that not only gamma-ray telescopes like Fermi but also X-ray sky monitors such as Swift and MAXI are critical to test the canonical picture based on the single-zone modeling. We also suggest a two-zone model that can naturally satisfy the X-ray constraints while explaining the flaring neutrinos via either photomeson or hadronuclear processes.

Magnetic reconnection is invoked as an efficient particle accelerator in a variety of astrophysical sources of non-thermal high-energy radiation. With large-scale two-dimensional particle-in-cell simulations of relativistic reconnection (i.e., with magnetization $\sigma\gg1$) in pair plasmas, we study the long-term evolution of the power-law slope and high-energy cutoff of the spectrum of accelerated particles. We find that the high-energy spectral cutoff does not saturate at $\gamma_{\rm cut}\sim 4\sigma$, as claimed by earlier studies, but it steadily grows with time as long as the reconnection process stays active. At late times, the cutoff scales approximately as $\gamma_{\rm cut}\propto \sqrt{t}$, regardless of the flow magnetization and initial temperature. We show that the particles dominating the high-energy spectral cutoff reside in plasmoids, and in particular in a strongly magnetised ring around the plasmoid core. The growth of their energy is driven by the increase in the local field strength, coupled with the conservation of the first adiabatic invariant. We also find that the power-law slope of the spectrum ($p=-{\rm d}\log N/{\rm d}\log \gamma$) evolves with time. For $\sigma\gtrsim10$, the spectrum is hard at early times ($p\lesssim 2$), but it tends to asymptote to $p\sim 2$; the steepening of the power-law slope allows the spectral cutoff to extend to higher and higher energies, without violating the fixed energy budget of the system. Our results demonstrate that relativistic reconnection is a viable candidate for accelerating the high-energy particles emitting in relativistic astrophysical sources.

Frequency-dependent brightness fluctuations of radio sources, the so-called extreme scattering events (ESEs), have been observed over the last three decades. They are caused by Galactic plasma structures whose geometry and origin are still poorly understood. In this paper, we construct axisymmentric two-dimensional (2D) column density profiles for the plasma lens and explore the resulting ESEs for both point-like and extended sources. A quantity that becomes relevant is the impact parameter $b$, namely the distance of the observer's path from the lens symmetry axis. We demonstrate its effects on the shape of ESE light curves and use it for a phenomenological classification of ESEs into four main types. Three of them are unique outcomes of the 2D model and do not show a characteristic U-shaped dip in the light curve, which has been traditionally used as an identification means of ESEs. We apply our model to five well-studied ESEs and show that elongated plasma tubes or quasi-spherical clouds are favoured over plasma sheets for four of them, while the remaining one is compatible with both lens geometries.

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.

We report an optical/UV jet and counterjet in M84, previously unreported in archival HST imaging. With archival VLA, ALMA, and Chandra imaging, we examine the first well-sampled spectral energy distribution of the inner jet of M84, where we find that multiple co-spatial spectral components are required. In particular, the ALMA data reveal that the radio spectrum of all four knots in the jet turns over at approximately 100 GHz, which requires a second component for the bright optical/UV emission. Further, the optical/UV has a soft spectrum and is inconsistent with the relatively flat X-ray spectrum, which indicates a third component at higher energies. Using archival VLA imaging, we have measured the proper motion of the innermost knots at 0.9+/-0.6 and 1.1+/-0.4 c, which when combined with the low jet-to-counterjet flux ratio yields an orientation angle for the system of 74 (+9,-18) degrees. In the radio, we find high fractional polarization of the inner jet of up to 30% while in the optical no polarization is detected (< 8%). We investigate different scenarios for explaining the particular multi-component SED of the knots. Inverse Compton models are ruled out due to the extreme departure from equipartition and the unrealistically high total jet power required. The multi-component SED can be naturally explained within a leptohadronic scenario, but at the cost of very high power in relativistic protons. A two-component synchrotron model remains a viable explanation, but more theoretical work is needed to explain the origin and properties of the electron populations.

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.

We report on the X-ray and optical properties of two high-mass X-ray binary systems located in the Large Magellanic Cloud (LMC). Based on the obtained optical spectra, we classify the massive companion as a supergiant star in both systems. Timing analysis of the X-ray events collected by XMM-Newton revealed the presence of coherent pulsations (spin period $\sim$2013 s) for XMMU J053108.3-690923 and fast flaring behaviour for XMMU J053320.8-684122. The X-ray spectra of both systems can be modelled sufficiently well by an absorbed power-law, yielding hard spectra and high intrinsic absorption from the environment of the systems. Due to their combined X-ray and optical properties we classify both systems as SgXRBs: the 19$^{\rm th}$ confirmed X-ray pulsar and a probable supergiant fast X-ray transient in the LMC, the second such candidate outside our Galaxy.

PSR J2032+4127 is a young and rapidly rotating pulsar on a highly eccentric orbit around the high-mass Be star MT91 213. X-ray monitoring of the binary system over a $\sim$4000 day period with Swift has revealed an increase of the X-ray luminosity which we attribute to the synchrotron emission of the shocked pulsar wind. We use Swift X-ray observations to infer a clumpy stellar wind with $r^{-2}$ density profile and constrain the Lorentz factor of the pulsar wind to $10^5<\gamma_w<10^6$. We investigate the effects of an axisymmetric stellar wind with polar gradient on the X-ray emission. Comparison of the X-ray light curve hundreds of days before and after the periastron can be used to explore the polar structure of the wind.

Plasmoids, quasi-spherical regions of plasma containing magnetic fields and high-energy particles, are a self-consistent by-product of the reconnection process in the relativistic regime. Recent two-dimensional particle-in-cell (PIC) simulations have shown that plasmoids can undergo a variety of processes (e.g. mergers, bulk acceleration, growth, and advection) within the reconnection layer. We developed a Monte Carlo (MC) code, benchmarked with the recent PIC simulations, to examine the effects of these processes on the steady-state size and momentum distributions of the plasmoid chain. The differential plasmoid size distribution is shown to be a power law, $N(w)\propto w^{-\chi}$, ranging from a few plasma skin depths to $\sim 0.1$ of the reconnection layer's length. We demonstrate numerically and analytically that the power law slope $\chi$ is linearly dependent upon the ratio of the plasmoid acceleration and growth rates and that it slightly decreases (from $\sim 2$ to $\sim 1.3$) with increasing plasma magnetization (from 3 to 50). We perform a detailed comparison of our results with those of recent PIC simulations and briefly discuss the astrophysical implications of our findings through the representative case of flaring events from blazar jets.

Jets in long-duration $\gamma$-ray bursts (GRBs) have to drill through the collapsing star in order to break out of it and produce the $\gamma$-ray signal while the central engine is still active. If the breakout time is shorter for more powerful engines, then the jet-collapsar interaction acts as a filter of less luminous jets. We show that the observed broken power-law GRB luminosity function is a natural outcome of this process. For a theoretically motivated breakout time that scales with jet luminosity as $L^{-\chi}$ with $\chi\sim 1/3-1/2$, we show that the shape of the $\gamma$-ray duration distribution can be uniquely determined by the GRB luminosity function and matches the observed one. This analysis has also interesting implications about the supernova-central engine connection. We show that not only successful jets can deposit sufficient energy in the stellar envelope to power the GRB-associated supernovae, but also failed jets may operate in all Type Ib/c supernovae.

Type IIn supernovae (SNe), a rare subclass of core collapse SNe, explode in dense circumstellar media that have been modified by the SNe progenitors at their last evolutionary stages. The interaction of the freely expanding SN ejecta with the circumstellar medium gives rise to a shock wave propagating in the dense SN environment, which may accelerate protons to multi-PeV energies. Inelastic proton-proton collisions between the shock-accelerated protons and those of the circumstellar medium lead to multi-messenger signatures. Here, we evaluate the possible neutrino signal of type IIn SNe and compare with IceCube observations. We employ a Monte Carlo method for the calculation of the diffuse neutrino emission from the SN IIn class to account for the spread in their properties. The cumulative neutrino emission is found to be ~ 10 per cent of the observed IceCube neutrino flux above 60 TeV. Type IIn SNe would be the dominant component of the diffuse astrophysical flux, only if 4 per cent of all core collapse SNe were of this type and 20 to 30 per cent of the shock energy was channeled to accelerated protons. Lower values of the acceleration efficiency are accessible by the observation of a single type IIn SN as a neutrino point source with IceCube using up-going muon neutrinos. Such an identification is possible in the first year following the SN shock breakout for sources within 20 Mpc.

Radiatively inefficient accretion flow models have been shown to accurately account for the spectrum and luminosity observed from Sgr A* in the X-ray regime down to mm wavelengths. However, observations at a few GHz cannot be explained by thermal electrons alone but require the presence of an additional non-thermal particle population. Here, we propose a model for the origin of such a population in the accretion flow via means of a pulsar orbiting the supermassive black hole in our Galaxy. Interactions between the relativistic pulsar wind with the disc lead to the formation of a bow shock in the wind. During the pulsar's transit through the accretion disc, relativistic pairs, accelerated at the shock front, are injected into the disc. The radio-emitting particles are long-lived and remain within the disc long after the pulsar's transit. Periodic pulsar transits through the disc result in regular injection episodes of non-thermal particles. We show that for a pulsar with spin-down luminosity $L_{\rm sd} \sim 3\times 10^{35}$ erg s$^{-1}$ and a wind Lorentz factor of $\gamma_{\rm w} \sim 10^4$ a quasi-steady synchrotron emission is established with luminosities in the $1-10$ GHz range comparable to the observed one.

The flat spectrum radio quasar 3C 279 is a known $\gamma$-ray variable source that has recently exhibited minute-scale variability at energies $>100$ MeV. One-zone leptonic models for blazar emission are severely constrained by the short timescale variability that implies a very compact emission region at a distance of hundreds of Schwarzschild radii from the central black hole. Here, we investigate a hadronic scenario where GeV $\gamma$-rays are produced via proton synchrotron radiation. We also take into account the effects of the hadronically initiated electromagnetic cascades (EMC). For a $\gamma$-ray emitting region in rough equipartition between particles and kG magnetic fields, located within the broad-line region (BLR), the development of EMC redistributes the $\gamma$-ray luminosity to softer energy bands and eventually leads to broad-band spectra that differ from the observed ones. Suppression of EMC and energy equipartition are still possible, if the $\gamma$-ray emitting region is located beyond the BLR, is fast moving with Doppler factor ($>70$), and contains strong magnetic fields ($>100$ G). Yet, these conditions cannot be easily met in parsec-scale jets, thus disfavouring a proton synchrotron origin of the Fermi-LAT flare.

Ap Librae is one out of a handful of low-frequency peaked blazars to be detected at TeV $\gamma$-rays and the only one with an identified X-ray jet. Combined observations of Fermi-LAT at high energies (HE) and of H.E.S.S at very high energies (VHE) revealed a striking spectral property of Ap Librae; the presence of a broad high-energy component that extends more than nine orders of magnitude in energy and is, therefore, hard to be explained by the usual single-zone synchrotron self-Compton model. We show that the superposition of different emission components related to photohadronic interactions can explain the $\gamma$-ray emission of Ap Librae without invoking external radiation fields. We present two indicative model fits to the spectral energy distribution of Ap Librae where the VHE emission is assumed to originate from a compact, sub-pc scale region of the jet. A robust prediction of our model is VHE flux variability on timescales similar to those observed at X-rays and HE $\gamma$-rays, which can be further used to distinguish between a sub-pc or kpc scale origin of the TeV emission. We thus calculate the expected variability signatures at X-rays, HE and VHE $\gamma$-rays and show that quasi-simultaneous flares are expected, with larger amplitude flares appearing at $\gamma$-rays. We assess the detectability of VHE variability from Ap Librae with CTA, next generation of IACTs. We show that $\sim$hr timescale variability at $E_{\gamma}>0.1$ TeV could be detectable at high significance with shorter exposure times than current Cherenkov telescopes.

We present our results on the {\gamma}-ray emission from interaction-powered supernovae (SNe), a recently discovered SN type that is suggested to be surrounded by a circumstellar medium (CSM) with densities 10^7-10^12~ cm^-3. Such high densities favor inelastic collisions between relativistic protons accelerated in the SN blast wave and CSM protons and the production of {\gamma}-ray photons through neutral pion decays. Using a numerical code that includes synchrotron radiation, adiabatic losses due to the expansion of the source, photon-photon interactions, proton-proton collisions and proton-photon interactions, we calculate the multi-wavelength non-thermal photon emission soon after the shock breakout and follow its temporal evolution until 100-1000 days. Focusing on the {\gamma}-ray emission at >100 MeV, we show that this could be detectable by the Fermi-LAT telescope for nearby (<10 Mpc) SNe with dense CSM (>10^11 cm^-3).

Powerful flares from blazars with short ($\sim$ min) variability timescales are challenging for current models of blazar emission. Here, we present a physically motivated ab initio model for blazar flares based on the results of recent particle-in-cell (PIC) simulations of relativistic magnetic reconnection. PIC simulations demonstrate that quasi-spherical plasmoids filled with high-energy particles and magnetic fields are a self-consistent by-product of the reconnection process. By coupling our PIC-based results (i.e., plasmoid growth, acceleration profile, particle and magnetic content) with a kinetic equation for the evolution of the electron distribution function we demonstrate that relativistic reconnection in blazar jets can produce powerful flares whose temporal and spectral properties are consistent with the observations. In particular, our model predicts correlated synchrotron and synchrotron self-Compton flares of duration of several hours--days powered by the largest and slowest moving plasmoids that form in the reconnection layer. Smaller and faster plasmoids produce flares of sub-hour duration with higher peak luminosities than those powered by the largest plasmoids. Yet, the observed fluence in both types of flares is similar. Multiple flares with a range of flux-doubling timescales (minutes to several hours) observed over a longer period of flaring activity (days or longer) may be used as a probe of the reconnection layer's orientation and the jet's magnetization. Our model shows that blazar flares are naturally expected as a result of magnetic reconnection in a magnetically-dominated jet.

Blobs, or quasi-spherical emission regions containing relativistic particles and magnetic fields, are often assumed ad hoc in emission models of relativistic astrophysical jets, yet their physical origin is still not well understood. Here, we employ a suite of large-scale two-dimensional particle-in-cell simulations in electron-positron plasmas to demonstrate that relativistic magnetic reconnection can naturally account for the formation of quasi-spherical plasmoids filled with high-energy particles and magnetic fields. Our simulations extend to unprecedentedly long temporal and spatial scales, so we can capture the asymptotic physics independently of the initial setup. We characterize the properties of the plasmoids that are continuously generated as a self-consistent by-product of the reconnection process: they are in rough energy equipartition between particles and magnetic fields; the upper energy cutoff of the plasmoid particle spectrum is proportional to the plasmoid width w, corresponding to a Larmor radius ~0.2 w; the plasmoids grow in size at ~0.1 of the speed of light, with most of the growth happening while they are still non-relativistic (first they grow); their growth is suppressed once they get accelerated to relativistic speeds by the field line tension, up to the Alfven speed (then they go). The largest plasmoids, whose typical recurrence interval is ~2.5 L/c, reach a characteristic size w ~ 0.2 L independently of the system length L, they have nearly isotropic particle distributions and they contain the highest energy particles, whose Larmor radius is ~0.03 L. The latter can be regarded as the Hillas criterion for relativistic reconnection. We briefly discuss the implications of our results for the high-energy emission from relativistic jets and pulsar winds.

Sgr A* is an ideal target to study low-luminosity accreting systems. It has been recently proposed that properties of the accretion flow around Sgr A* can be probed through its interactions with the stellar wind of nearby massive stars belonging to the S-cluster. When a star intercepts the accretion disk, the ram and thermal pressures of the disk terminate the stellar wind leading to the formation of a bow shock structure. Here, a semi-analytical model is constructed which describes the geometry of the termination shock formed in the wind. With the employment of numerical hydrodynamic simulations, this model is both verified and extended to a region prone to Kelvin-Helmholtz instabilities. Because the characteristic wind and stellar velocities are in $\sim10^{8}$ cm s$^{-1}$ range, the shocked wind may produce detectable X-rays via thermal bremsstrahlung emission. The application of this model to the pericenter passage of S2, the brightest member of the S-cluster, shows that the shocked wind produces roughly a month long X-ray flare with a peak luminosity of $L\approx 4 \times 10^{33}$ erg s$^{-1}$ for a stellar mass-loss rate, disk number density, and thermal pressure strength of $\dot{M}_{\rm w}= 10^{-7} M_\odot\, {\rm yr}^{-1}$, $n_{\rm d} = 10^{5}$ cm$^{-3}$, and $\alpha=0.1$, respectively. This peak luminosity is comparable to the quiescent X-ray emission detected from Sgr A* and is within the detection capabilities of current X-ray observatories. Its detection could constrain the density and thickness of the disk at a distance of $\sim 3000$ gravitational radii from the supermassive black hole.

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.

Several supernovae (SNe) with an unusually dense circumstellar medium (CSM) have been recently observed at radio frequencies. Their radio emission is powered by relativistic electrons that can be either accelerated at the SN shock (primaries) or injected as a by-product (secondaries) of inelastic proton-proton collisions. We investigate the radio signatures from secondary electrons, by detailing a semi-analytical model to calculate the temporal evolution of the distributions of protons, primary and secondary electrons. With our formalism, we track the cooling history of all the particles that have been injected into the emission region up to a given time, and calculate the resulting radio spectra and light curves. For a SN shock propagating through the progenitor wind, we find that secondary electrons control the early radio signatures, but their contribution decays faster than that of primary electrons. This results in a flattening of the light curve at a given radio frequency that depends only upon the radial profiles of the CSM density and of the shock velocity, $v_0$. The relevant transition time at the peak frequency is $\sim 190 \, K_{\rm ep,-3}^{-1} A_{\rm w, 16} \beta_{0, -1.5}^{-2}\,{\rm d}$, where $A_{\rm w}$ is the wind mass-loading parameter, $\beta_0=v_0/c$ and $K_{\rm ep}$ is the electron-to-proton ratio of accelerated particles. We explicitly show that late peak times at 5 GHz (i.e., $t_{\rm pk} \gtrsim 300-1000$ d) suggest a shock wave propagating in a dense wind ($A_{\rm w} \gtrsim 10^{16}-10^{17}$ gr cm$^{-1}$), where secondary electrons are likely to power the observed peak 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 report on the first detection of X-ray dust scattered rings from the Galactic low mass X-ray binary V404 Cyg. The observation of the system with Swift/XRT on June 30 2015 revealed the presence of five concentric ring-like structures centred at the position of V404 Cyg. Follow-up Swift/XRT observations allowed a time-dependent study of the X-ray rings. Assuming that these are the result of small-angle, single X-ray scattering by dust grains along the line of sight, we find that their angular size scales as $\theta \propto\sqrt{t}$ in agreement with theoretical predictions. The dust grains are concentrated in five dust layers located at about 2.12, 2.05, 1.63, 1.50 and 1.18 kpc from the observer. These coincide roughly with locations of enhanced extinction as determined by infrared photometry. Assuming that the grain size distribution is described by a generalized Mathis-Rumpl-Nordsieck model, we find that the power-law index of the most distant cloud is $q\sim 4.4$, while $q \sim 3.5-3.7$ in all other clouds. We constrain at a $3\sigma$ level the maximum grain size of the intermediate dust layers in the range $0.16-0.20\,\mu$m and set a lower limit of $\sim 0.2\,\mu$m in the other clouds. Hints of an exponential cutoff at the angular intensity profile of the outermost X-ray ring suggest that the smallest grains have sizes $0.01 \mu{\rm m}\le \alpha_{\min} \lesssim 0.03\,\mu$m. Based on the relative ratios of dust column densities we find the highest dust concentration at $\sim 1.6$ kpc. Our results indicate a gradient in the dust properties within 1 kpc from V404 Cyg.

The inverse Compton catastrophe is defined as a dramatic rise in the luminosity of inverse Compton scattered photons. It is described by a non-linear loop of radiative processes that sets in for high values of the electron compactness and is responsible for the efficient transfer of energy from electrons to photons, predominantly through inverse Compton scatterings. We search for the conditions that drive a magnetized non-thermal source to the inverse Compton catastrophe regime and study its multi-wavelength (MW) photon spectrum. We develop a generic analytical framework and use numerical calculations as a backup to the analytical predictions. We find that the escaping radiation from a source in the Compton catastrophe regime bears some unique features. The MW photon spectrum is a broken power law with a break at $\sim m_e c^2$ due to the onset of the Klein-Nishina suppression. The spectral index below the break energy depends on the electron and magnetic compactnesses logarithmically, while it is independent of the electron power-law index ($s$). The maximum radiating power emerges typically in the $\gamma$-ray regime, at energies $\sim m_e c^2$ ($\sim \gamma_{\max} m_e c^2$ ) for $s>2$ ($s\lesssim 2$), where $\gamma_{\max}$ is the maximum Lorentz factor of the injected electron distribution. We apply the principles of the inverse Compton catastrophe to blazars and $\gamma$-ray bursts using the analytical framework we developed, and show how these can be used to impose robust constraints on the source parameters.

Blazars have been suggested as possible neutrino sources long before the recent IceCube discovery of high-energy neutrinos. We re-examine this possibility within a new framework built upon the blazar simplified view and a self-consistent modelling of neutrino emission from individual sources. The former is a recently proposed paradigm that explains the diverse statistical properties of blazars adopting minimal assumptions on blazars' physical and geometrical properties. This view, tested through detailed Monte Carlo simulations, reproduces the main features of radio, X-ray, and gamma-ray blazar surveys and also the extragalactic gamma-ray background at energies > 10 GeV. Here we add a hadronic component for neutrino production and estimate the neutrino emission from BL Lacs as a class, "calibrated" by fitting the spectral energy distributions of a preselected sample of BL Lac objects and their (putative) neutrino spectra. Unlike all previous papers on this topic, the neutrino background is then derived by summing up at a given energy the fluxes of each BL Lac in the simulation, all characterised by their own redshift, synchrotron peak energy, gamma-ray flux, etc. Our main result is that BL Lacs as a class can explain the neutrino background seen by IceCube above ~ 0.5 PeV while they only contribute ~ 10% at lower energies, leaving room to some other population(s)/physical mechanism. However, one cannot also exclude the possibility that individual BL Lacs still make a contribution at the ~ 20% level to the IceCube low-energy events. Our scenario makes specific predictions testable in the next few years.

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.

Observations of gamma-ray-bursts and jets from active galactic nuclei reveal that the jet flow is characterized by a high radiative efficiency and that the dissipative mechanism must be a powerful accelerator of non-thermal particles. Shocks and magnetic reconnection have long been considered as possible candidates for powering the jet emission. Recent progress via fully-kinetic particle-in-cell simulations allows us to revisit this issue on firm physical grounds. We show that shock models are unlikely to account for the jet emission. In fact, when shocks are efficient at dissipating energy, they typically do not accelerate particles far beyond the thermal energy, and vice versa. In contrast, we show that magnetic reconnection can deposit more than 50% of the dissipated energy into non-thermal leptons as long as the energy density of the magnetic field in the bulk flow is larger than the rest mass energy density. The emitting region, i.e., the reconnection downstream, is characterized by a rough energy equipartition between magnetic fields and radiating particles, which naturally accounts for a commonly observed property of blazar jets.

In 2012 Markarian 421 underwent the largest flare ever observed in this blazar at radio frequencies. In the present study, we start exploring this unique event and compare it to a less extreme event in 2013. We use 15 GHz radio data obtained with the Owens Valley Radio Observatory 40-m telescope, 95 GHz millimeter data from the Combined Array for Research in Millimeter-Wave Astronomy, and GeV gamma-ray data from the Fermi Gamma-ray Space Telescope. The radio light curves during the flaring periods in 2012 and 2013 have very different appearances, both in shape and peak flux density. Assuming that the radio and gamma-ray flares are physically connected, we attempt to model the most prominent sub-flares of the 2012 and 2013 activity periods by using the simplest possible theoretical framework. We first fit a one-zone synchrotron self-Compton (SSC) model to the less extreme 2013 flare and estimate parameters describing the emission region. We then model the major gamma-ray and radio flares of 2012 using the same framework. The 2012 gamma-ray flare shows two distinct spikes of similar amplitude, so we examine scenarios associating the radio flare with each spike in turn. In the first scenario, we cannot explain the sharp radio flare with a simple SSC model, but we can accommodate this by adding plausible time variations to the Doppler beaming factor. In the second scenario, a varying Doppler factor is not needed, but the SSC model parameters require fine tuning. Both alternatives indicate that the sharp radio flare, if physically connected to the preceding gamma-ray flares, can be reproduced only for a very specific choice of parameters.

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.

We present the spectral signatures of the Bethe-Heitler pair production ($pe$) process on the spectral energy distribution (SED) of blazars, in scenarios where the hard $\gamma$-ray emission is of photohadronic origin. If relativistic protons interact with the synchrotron blazar photons producing $\gamma$ rays through photopion processes, we show that, besides the $2-20$ PeV neutrino emission, the typical blazar SED should have an emission feature due to the synchrotron emission of $pe$ secondaries that bridges the gap betweeen the low-and high-energy humps of the SED, namely in the energy range 40 keV-40 MeV. We first present analytical expressions for the photopion and $pe$ loss rates in terms of observable quantities of blazar emission. For the $pe$ loss rate in particular, we derive a new approximate analytical expression for the case of a power-law photon distribution, which has an excellent accuracy with the numerically calculated exact one, especially at energies above the threshold for pair production. We show that for typical blazar parameters, the photopair synchrotron emission emerges in the hard X-ray/soft $\gamma$-ray energy range with a characteristic spectral shape and non negligible flux, which may be even comparable to the hard $\gamma$-ray flux produced through photopion processes. We argue that the expected "$pe$ bumps" are a natural consequence of leptohadronic models, and as such, they may indicate that blazars with a three-hump SED are possible emitters of high-energy neutrinos.

We apply both leptonic and leptohadronic emission scenarios for modelling the multiwavelength photon spectra and the observed variability in the optical, X-ray, and TeV gamma-ray energy bands of blazar PKS 2155-304 while being in a low state between 25 August and 6 September 2008. We consider three emission models, i.e. a one-component synchrotron self-Compton model (1-SSC), a one-zone proton synchrotron model (LHs), and a two-component SSC model (2-SSC). Only in the first scenario can the emission from the optical up to TeV gamma-rays be attributed to a single particle population from one emission region. Using a time-dependent numerical code that solves the kinetic equations for each particle species, we derived, in all cases, acceptable fits to the time-averaged SED. By imposing variations to one (or more) model parameters according to observed variability pattern in one (or more) frequencies we calculated the respective light curves and compared them with the observations. We show that the 1-SSC model cannot account for the anticorrelation observed between the X-rays and VHE gamma-rays, although it can explain the time-averaged SED. The anticorrelation can be more naturally explained by the two-component emission models. Both of them reproduce satisfactorily the optical, X-ray, and TeV variability but at the cost of additional free parameters, which from four in the 2-SSC model increase to six in the LHs model. Although the results of our time-resolved analysis do not favour one of the aforementioned models, they suggest that a two-component scenario is more adequate for the emission of PKS 2155-304 in the low state of 2008, which agrees with a recent independent analysis. This suggests that the quiescent blazar radiation might result from a superposition of the radiation from different components, while a flare might still be the result of a single component.

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.

We investigate the effects of hadronic cascades on the gamma-ray burst (GRB) prompt emission spectra in scenarios of efficient neutrino production. By assuming a fiducial GRB spectrum and a power-law proton distribution extending to ultra-high energies, we calculate the proton cooling rate and the neutrino emission produced through photopion processes. For this, we employ a numerical code that follows the formation of the hadronic cascade by taking into account non-linear feedback effects, such as the evolution of the target photon field itself due to the contribution of secondary particles. We show that in cases of efficient proton cooling and subsequently efficient high-energy neutrino production, the emission from the hadronic cascade distorts and may even dominate the GRB spectrum. Taking this into account, we constrain the allowable values of the ratio $\eta_p=L_p/L_{\gamma}$, where $L_{p}$ and $L_{\gamma}$ are the isotropic equivalent proton and prompt gamma-ray luminosities. For the highest value of $\eta_{p}$ that does not lead to the dominance of the cascading emission, we then calculate the maximum neutrino luminosity from a single burst and show that it ranges between $(0.01- 0.6)L_{p}$ and $(0.5-1.4)L_{\gamma}$ for various parameter sets. We discuss possible implications of other parameters, such as the magnetic field strength and the shape of the initial gamma-ray spectrum, on our results. Finally, we compare the upper limit on $\eta_p$ derived here with various studies in the field, and we point out the necessity of a self-consistent treatment of the hadronic emission in order to avoid erroneously high neutrino fluxes from GRB models.

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 study a mechanism for producing intrinsic broken power-law gamma-ray spectra in compact sources. This is based on the principles of automatic photon quenching, according to which, gamma-rays are being absorbed on spontaneously produced soft photons, whenever the injected luminosity in gamma-rays lies above a certain critical value. We derive an analytical expression for the critical gamma-ray compactness in the case of power-law injection. For the case where automatic photon quenching is relevant, we calculate analytically the emergent steady-state gamma-ray spectra. We perform also numerical calculations in order to back up our analytical results. We show that a spontaneously quenched power-law gamma-ray spectrum obtains a photon index 3{\Gamma}/2, where {\Gamma} is the photon index of the power-law at injection. Thus, large spectral breaks of the gamma-ray photon spectrum, e.g. $\Delta \Gamma \gtrsim 1$, can be obtained by this mechanism. We also discuss additional features of this mechanism that can be tested observationally. Finally, we fit the multiwavelength spectrum of a newly discovered blazar (PKS 0447-439) by using such parameters, as to explain the break in the gamma-ray spectrum by means of spontaneous photon quenching, under the assumption that its redshift lies in the range 0.1<z<0.24.

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.

In the present work we propose an innovative estimation method for the minimum Doppler factor and energy content of the gamma-ray emitting region of quasar 3C 279, using a standard proton synchrotron blazar model and the principles of automatic photon quenching. The latter becomes relevant for high enough magnetic fields and results in spontaneous annihilation of gamma-rays. The absorbed energy is then redistributed into electron-positron pairs and soft radiation. We show that as quenching sets an upper value for the source rest-frame gamma-ray luminosity, one has, by neccessity, to resort to Doppler factors that lie above a certain value in order to explain the TeV observations. The existence of this lower limit for the Doppler factor has also implications on the energetics of the emitting region. In this aspect, the proposed method can be regarded as an extension of the widely used one for estimating the equipartition magnetic field using radio observations. In our case, the leptonic synchrotron component is replaced by the proton synchrotron emission and the radio by the VHE gamma-ray observations. We show specifically that one can model the TeV observations by using parameter values that minimize both the energy density and the jet power at the cost of high-values of the Doppler factor. On the other hand, the modelling can also be done by using the minimum possible Doppler factor; this, however, leads to a particle dominated region and high jet power for a wide range of magnetic field values. Despite the fact that we have focused on the case of 3C 279, our analysis can be of relevance to all TeV blazars favoring hadronic modelling that have, moreover, simultaneous X-ray observations.

The hadronic model of Active Galactic Nuclei and other compact high energy astrophysical sources assumes that ultra-relativistic protons, electron-positron pairs and photons interact via various hadronic and electromagnetic processes inside a magnetized volume, producing the multiwavelength spectra observed from these sources. A less studied property of such systems is that they can exhibit a variety of temporal behaviours due to the operation of different feedback mechanisms. We investigate the effects of one possible feedback loop, where \gamma-rays produced by photopion processes are being quenched whenever their compactness increases above a critical level. This causes a spontaneous creation of soft photons in the system that result in further proton cooling and more production of \gamma-rays, thus making the loop operate. We perform an analytical study of a simplified set of equations describing the system, in order to investigate the connection of its temporal behaviour with key physical parameters. We also perform numerical integration of the full set of kinetic equations verifying not only our analytical results but also those of previous numerical studies. We find that once the system becomes `supercritical', it can exhibit either a periodic behaviour or a damped oscillatory one leading to a steady state. We briefly point out possible implications of such a supercriticality on the parameter values used in Active Galactic Nuclei spectral modelling, through an indicative fitting of the VHE emission of blazar 3C 279.

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.

Aims: We investigate photon quenching in compact non-thermal sources. This involves photon-photon annihilation and lepton synchrotron radiation in a network that can become non-linear. As a result the gamma-ray luminosity of a source cannot exceed a critical limit that depends only on the radius of the source and on the magnetic field. Methods: We perform analytic and numerical calculations that verify previous results and extend them so that the basic properties of photon quenching are investigated. Results: We apply the above to the 2006 TeV observations of quasar 3C279 and obtain the parameter space of allowed values for the radius of the emitting source, its magnetic field strength and the Doppler factor of the flow. We argue that the TeV observations favour either a modest Doppler factor and a low magnetic field or a high Doppler factor and a high magnetic field.

Aims: We investigate the behavior of the frequency-centered light curves expected within the standard model of Gamma Ray Bursts allowing the maximum electron energy to be a free parameter permitted to take low values. Methods: We solve the spatially averaged kinetic equations which describe the simultaneous evolution of particles and photons, obtaining the multi-wavelength spectra as a function of time. From these we construct the frequency-centered light curves giving emphasis in the X-ray and optical bands. Results: We show that in cases where the maximum electron energy takes low values, the produced X-ray light curves show a plateau as the synchrotron component gives its place to the Synhro Self-Compton one in the X-ray band.

Aims: Drawing an analogy with Active Galactic Nuclei, we investigate the one-zone SSC model of Gamma Ray Bursts afterglows in the presence of electron injection and cooling both by synchrotron and SSC losses. Methods: We solve the spatially averaged kinetic equations which describe the simultaneous evolution of particles and photons, obtaining the multi-wavelength spectrum as a function of time. We back up our numerical calculations with analytical solutions of the equations using various profiles of the magnetic field evolution under certain simplifying assumptions. Results: We apply the model to the afterglow evolution of GRBs in a uniform density environment and examine the impact various parameters have on the multiwavelength spectra. We find that in cases where the electron injection and/or the ambient density is high, the losses are dominated by SSC and the solutions depart significantly from the ones derived in the synchrotron standard cases.