We introduce a new time-dependent lepto-hadronic model for blazar emission that takes into account the radiation emitted by secondary particles, such as pions and muons, from photo hadronic interactions. Starting from a baseline parameter set guided by a fit to the spectral energy distribution of the blazar 3C 279, we perform a parameter study to investigate the effects of perturbations of the input parameters to mimic different flaring events to study the resulting light curves in the optical, X-ray, high-energy (HE: > E 100 MeV), and very-high-energy > E ( 100 GeV) γ-rays as well as the neutrino emission associated with charged-pion and muon decay. We find that flaring events from an increase in the efficiency of Fermi II acceleration will produce a positive correlation between all bandpasses and a marked plateau in the HE γ-ray lightcurve. We also predict a distinctive dip in the HE lightcurve for perturbations caused by a change in the proton injection spectral index. These plateaus/dips could be a tell tale signature of hadronic models for perturbations that lead to more efficient acceleration of highenergy protons in parameter regimes where pion and muon synchrotron emission is non-negligible.

The well studied blazar Markarian 421 (Mrk 421, z=0.031) was the subject of an intensive multi-wavelength campaign when it flared in 2013 April. The recorded X-ray and very high energy (VHE, E>100 GeV) γ-ray fluxes are the highest ever measured from this object. At the peak of the activity, it was monitored by the hard X-ray focusing telescope Nuclear Spectroscopic Telescope Array (NuSTAR) and Swift X-Ray Telescope (XRT). In this work, we present a detailed variability analysis of NuSTAR and Swift-XRT observations of Mrk 421 during this flaring episode. We obtained the shortest flux doubling time of 14.01±5.03 minutes, which is the shortest hard X-ray (3−79 keV) variability ever recorded from Mrk 421 and is on the order of the light crossing time of the black hole's event horizon. A pattern of extremely fast variability events superposed on slowly varying flares is found in most of the NuSTAR observations. We suggest that these peculiar variability patterns may be explained by magnetic energy dissipation and reconnection in a fast moving compact emission region within the jet.Based on the fast variability, we derive a lower limit on the magnetic field strength of B ≥ 0.73δ

Aims. In a previous study we suggested that the broad-band emission and variability properties of BL Lacertae can be accounted for by a double synchrotron emission component with related inverse-Compton emission from the jet, plus thermal radiation from the accretion disc. Here we investigate the matter with further data extending over a wider energy range. Methods. The GLAST-AGILE Support Program (GASP) of the whole earth blazar telescope (WEBT) monitored BL Lacertae in 2008-2009 at radio, near-IR, and optical frequencies to follow its flux behaviour. During this period, high-energy observations were performed by XMM-Newton, Swift, and Fermi. We analyse these data with particular attention to the calibration of Swift UV data, and apply a helical jet model to interpret the source broad-band variability.Results. The GASP-WEBT observations show an optical flare in 2008 February-March, and oscillations of several tenths of mag on a few-day time scale afterwards. The radio flux is only mildly variable. The UV data from both XMM-Newton and Swift seem to confirm a UV excess that is likely caused by thermal emission from the accretion disc. The X-ray data from XMM-Newton indicate a strongly concave spectrum, as well as moderate (∼4-7%) flux variability on an hour time scale. The Swift X-ray data reveal fast (interday) flux changes, not correlated with those observed at lower energies. We compare the spectral energy distribution (SED) corresponding to the 2008 low-brightness state, which was characterised by a synchrotron dominance, to the 1997 outburst state, where the inverse-Compton emission was prevailing. A fit with an inhomogeneous helical jet model suggests that two synchrotron components are at work with their self inverse-Compton emission. Most likely, they represent the radiation from two distinct emitting regions in the jet. We show that the difference between the source SEDs in 2008 and 1997 can be explained in terms of pure geometrical variations. The outburst state occurred when the jet-emitting regions were better aligned with the line of sight, producing an increase of the Doppler beaming factor. Conclusions. Our analysis demonstrates that the jet geometry can play an extremely important role in the BL Lacertae flux and spectral variability. Indeed, the emitting jet is probably a bent and dynamic structure, and hence changes in the emitting regions viewing angles are likely to happen, with strong consequences on the source multiwavelength behaviour.

No abstract

We present a newly developed time-dependent three-dimensional multi-zone hadronic blazar emission model. By coupling a Fokker-Planck based leptohadronic particle evolution code 3DHad with a polarization-dependent radiation transfer code, 3DPol, we are able to study the time-dependent radiation and polarization signatures of a hadronic blazar model for the first time. Our current code is limited to parameter regimes in which the hadronic γ-ray output is dominated by proton synchrotron emission, neglecting pion production. Our results demonstrate that the time-dependent flux and polarization signatures are generally dominated by the relation between the synchrotron cooling and the light crossing time scale, which is largely independent of the exact model parameters. We find that unlike the low-energy polarization signatures, which can vary rapidly in time, the high-energy polarization signatures appear stable. As a result, future high-energy polarimeters may be able to distinguish such signatures from the lower and more rapidly variable polarization signatures expected in leptonic models.