We present a newly developed time-dependent three-dimensional multi-zone hadronic blazar emission model. By coupling a Fokker-Planck based lepto-hadronic 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 $\gamma$-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.

In this study, we use a one-zone leptonic and a lepto-hadronic model to investigate the multi-wavelength emission and the prominent flare of the flat spectrum radio quasar 3C 454.3 in Nov 2010. We perform a parameter study with both models to obtain broadband fits to the spectral energy distribution of 3C 454.3. Starting with the baseline parameters obtained from the fits, we then investigate different flaring scenarios for both models to explain an extreme outburst and spectral hardening of 3C 454.3 that occurred in Nov 2010. We find that the one zone lepto-hadronic model can successfully explain both the broadband multi-wavelength spectral energy distribution and light curves in the optical R, Swift XRT and Fermi gamma-ray bandpasses for 3C 454.3 during quiescence and the peak of the Nov. 2010 flare. We also find that the one-zone leptonic model produces poor fits to the broadband spectra in the X-ray and high energy gamma-ray band passes for the Nov. 2010 flare.

The blazar 3C 279 exhibited a giant $\gamma$-ray outburst in 2013 December and 2014 April. Apart from the very fast $\gamma$-ray flux variability, the spectral nature of the flares were also found to vary significantly between these two flaring events. A prominent curvature in the $\gamma$-ray spectrum was noticed in 2014 April flare, on the other hand, the 2013 December displayed an extreme hardening of the spectrum. These observations, thus, put strong constraints on our understanding of the underlying particle acceleration mechanisms.

The blazar 3C 279 exhibited twin $\gamma$-ray flares of similar intensity in 2013 December and 2014 April. In this work, we present a detailed multi-wavelength analysis of the 2013 December flaring event. Multi-frequency observations reveal the uncorrelated variability patterns with X-ray and optical-UV fluxes peaking after the $\gamma$-ray maximum. The broadband spectral energy distribution (SED) at the peak of the $\gamma$-ray activity shows a rising $\gamma$-ray spectrum but a declining optical-UV flux. This observation along with the detection of uncorrelated variability behavior rules out the one-zone leptonic emission scenario. We, therefore, adopt two independent methodologies to explain the SED: a time dependent lepto-hadronic modeling and a two-zone leptonic radiative modeling approach. In the lepto-hadronic modeling, a distribution of electrons and protons subjected to a randomly orientated magnetic field produces synchrotron radiation. Electron synchrotron is used to explain the IR to UV emission while proton synchrotron emission is used to explain the high energy $\gamma$-ray emission. A combination of both electron synchrotron self Compton emission and proton synchrotron emission is used to explain the X-ray spectral break seen during the later stage of the flare. In the two-zone modeling, we assume a large emission region emitting primarily in IR to X-rays and $\gamma$-rays to come primarily from a fast moving compact emission region. We conclude by noting that within a span of 4 months, 3C 279 has shown the dominance of a variety of radiative processes over each other and this reflects the complexity involved in understanding the physical properties of blazar jets in general.

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) $\gamma$-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 {\it Nuclear Spectroscopic Telescope Array} ({\it NuSTAR}) and {\it Swift} X-Ray Telescope (XRT). In this work, we present a detailed variability analysis of {\it NuSTAR} and {\it Swift}-XRT observations of Mrk 421 during this flaring episode. We obtained the shortest flux doubling time of 14.01$\pm$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 {\it 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 \ge 0.73 \delta_1^{-2/3} \, \nu_{19}^{1/3}$~G, where $\delta_1$ is the Doppler factor in units of 10, and $\nu_{19}$ is the characteristic X-ray synchrotron frequency in units of $10^{19}$~Hz.

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 lightcurves in the optical, X-ray, high energy (HE: E > 100 MeV) and very-high-energy (VHE: E > 100 GeV) gamma-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 gamma-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 high energy protons in parameter regimes where pion and muon synchrotron emission is non-negligible.

In this paper, we discuss the light-curve features of various flaring scenarios in a time-dependent leptonic model for low-frequency-peaked blazars. The quasar 3C273 is used as an illustrative example. Our code takes into account Fermi-II acceleration and all relevant electron cooling terms, including the external radiation fields generally found to be important in the modeling of the SEDs of FSRQs, as well as synchrotron self absorption and gamma-gamma pair-production. General parameters are constrained through a fit to the average spectral energy distribution (SED) of the blazar by numerically solving the time-dependent Fokker-Planck equation for the electron evolution in a steady-state situation. We then apply perturbations to several input parameters (magnetic field, particle injection luminosity, acceleration time scale) to simulate flaring events and compute time-dependent SEDs and light curves in representative energy bands (radio, optical, X-rays, gamma-rays). Time lags between different bands are evaluated using a discrete cross correlation analysis. We find that Fermi-II acceleration has a significant effect on the distributions and that flaring events caused by increased acceleration efficiency of the Fermi II process will produce a correlation between the radio, optical and gamma-ray bandpasses, but an anti-correlation between these three bandpasses and the X-ray band, with the X-rays lagging behind the variations in other bands by up to several hours.

The blazar AO 0235+164 (z = 0.94) has been one of the most active objects observed by Fermi Large Area Telescope (LAT) since its launch in Summer 2008. In addition to the continuous coverage by Fermi, contemporaneous observations were carried out from the radio to {\gamma} -ray bands between 2008 September and 2009 February. In this paper, we summarize the rich multi-wavelength data collected during the campaign (including F-GAMMA, GASP- WEBT, Kanata, OVRO, RXTE, SMARTS, Swift, and other instruments), examine the cross-correlation between the light curves measured in the different energy bands, and interpret the resulting spectral energy distributions in the context of well-known blazar emission models. We find that the {\gamma} -ray activity is well correlated with a series of near-IR/optical flares, accompanied by an increase in the optical polarization degree. On the other hand, the X-ray light curve shows a distinct 20 day high state of unusually soft spectrum, which does not match the extrapolation of the optical/UV synchrotron spectrum. We tentatively interpret this feature as the bulk Compton emission by cold electrons contained in the jet, which requires an accretion disk corona with an effective covering factor of 19% at a distance of 100 Rg . We model the broadband spectra with a leptonic model with external radiation dominated by the infrared emission from the dusty torus.

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. 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. 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 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. Our analysis demonstrates that the jet geometry can play an extremely important role in the BL Lacertae flux and spectral variability.