Multiwavelength Variability Signatures of Relativistic Shocks in Blazar Jets
Mildly-relativistic shocks that are embedded in colliding magnetohydrodynamic flows are prime sites for relativistic particle acceleration and production of strongly variable, polarized multi-wavelength emission from relativistic jet sources such as blazars and gamma-ray bursts. The principal energization mechanisms at these shocks are diffusive shock acceleration and shock drift acceleration. In recent work, we had self-consistently coupled shock acceleration and radiation transfer simulations in blazar jets in a basic one-zone scenario. These one-zone models revealed that the observed spectral energy distributions (SEDs) of blazars strongly constrain the nature of the shock layer hydromagnetic turbulence. In this paper, we expand our previous work by including full time dependence and treating two zones, one being the site of acceleration, and the second being a larger emission zone. This construction is applied to multiwavelength flares of the flat spectrum radio quasar 3C 279, fitting snap-shot SEDs and generating light curves that are consistent with observed variability timescales. We also present a generic study for the typical flaring behavior of the BL Lac object Mrk 501. The model predicts correlated variability across all wavebands, but cross-band time lags depending on the type of blazar (FSRQ vs. BL Lac), as well as distinctive spectral hysteresis patterns in all wavelength bands, from mm radio waves to gamma-rays. These evolutionary signatures serve to provide diagnostics on the competition between acceleration and radiative cooling.