Time dependent spectral modeling of Markarian 421 during a violent outburst in 2010

We present the results of extensive modeling of the spectral energy distributions (SEDs) of the closest blazar (z=0.031) Markarian 421 (Mrk 421) during a giant outburst in February 2010. The source underwent rapid flux variations in both X-rays and very high energy (VHE) gamma-rays as it evolved from a low-flux state on 2010 February 13-15 to a high-flux state on 2010 February 17. During this period, the source exhibited significant spectral hardening from X-rays to VHE gamma-rays while exhibiting a "harder when brighter" behavior in these energy bands. We reproduce the broadband SED using a time-dependent multi-zone leptonic jet model with radiation feedback. We find that an injection of the leptonic particle population with a single power-law energy distribution at shock fronts followed by energy losses in an inhomogeneous emission region is suitable for explaining the evolution of Mrk 421 from low- to high-flux state in February 2010. The spectral states are successfully reproduced by a combination of a few key physical parameters, such as the maximum $\&$ minimum cutoffs and power-law slope of the electron injection energies, magnetic field strength, and bulk Lorentz factor of the emission region. The simulated light curves and spectral evolution of Mrk 421 during this period imply an almost linear correlation between X-ray flux at 1-10 keV energies and VHE gamma-ray flux above 200 GeV, as has been previously exhibited by this source. Through this study, a general trend that has emerged for the role of physical parameters is that, as the flare evolves from a low- to a high-flux state, higher bulk kinetic energy is injected into the system with a harder particle population and a lower magnetic field strength.Comment: 13 pages, 5 figures, accepted for publication in MNRA

Stochastic model of optical variability of BL Lacertae

We use optical photometric and polarimetric data of BL Lacertae that cover a period of 22 years to study the variability of the source. The long-term observations are employed for establishing parameters of a stochastic model consisting of the radiation from a steady polarized source and a number of variable components with different polarization parameters, proposed by Hagen-Thorn et al. earlier. We infer parameters of the model from the observations using numerical simulations based on a Monte Carlo method, with values of each model parameter selected from a Gaussian distribution. We determine the best set of model parameters by comparing model distributions to the observational ones using the chi-square criterion. We show that the observed photometric and polarimetric variability can be explained within a model with a steady source of high polarization, ~40%, and with direction of polarization parallel to the parsec scale jet, along with 10+-5 sources of variable polarization.Comment: 4 pages, 10 figures, published by Astronomy and Astrophysics; v2: typos correcte

The connection between the radio jet and the gamma-ray emission in the radio galaxy 3C 120

We present the analysis of the radio jet evolution of the radio galaxy 3C 120 during a period of prolonged gamma-ray activity detected by the Fermi satellite between December 2012 and October 2014. We find a clear connection between the gamma-ray and radio emission, such that every period of gamma-ray activity is accompanied by the flaring of the mm-VLBI core and subsequent ejection of a new superluminal component. However, not all ejections of components are associated with gamma-ray events detectable by Fermi. Clear gamma-ray detections are obtained only when components are moving in a direction closer to our line of sight.This suggests that the observed gamma-ray emission depends not only on the interaction of moving components with the mm-VLBI core, but also on their orientation with respect to the observer. Timing of the gamma-ray detections and ejection of superluminal components locate the gamma-ray production to within almost 0.13 pc from the mm-VLBI core, which was previously estimated to lie about 0.24 pc from the central black hole. This corresponds to about twice the estimated extension of the broad line region, limiting the external photon field and therefore suggesting synchrotron self Compton as the most probable mechanism for the production of the gamma-ray emission. Alternatively, the interaction of components with the jet sheath can provide the necessary photon field to produced the observed gamma-rays by Compton scattering.Comment: Already accepted for publication in The Astrophysical Journa

Multiwavelength Variations of 3C 454.3 during the November 2010 to January 2011 Outburst

We present multiwavelength data of the blazar 3C 454.3 obtained during an extremely bright outburst from November 2010 through January 2011. These include flux density measurements with the Herschel Space Observatory at five submillimeter-wave and far-infrared bands, the Fermi Large Area Telescope at gamma-ray energies, Swift at X-ray, ultraviolet (UV), and optical frequencies, and the Submillimeter Array at 1.3 mm. From this dataset, we form a series of 52 spectral energy distributions (SEDs) spanning nearly two months that are unprecedented in time coverage and breadth of frequency. Discrete correlation anlaysis of the millimeter, far-infrared, and gamma-ray light curves show that the variations were essentially simultaneous, indicative of co-spatiality of the emission, at these wavebands. In contrast, differences in short-term fluctuations at various wavelengths imply the presence of inhomegeneities in physical conditions across the source. We locate the site of the outburst in the parsec-scale core, whose flux density as measured on 7 mm Very Long Baseline Array images increased by 70 percent during the first five weeks of the outburst. Based on these considerations and guided by the SEDs, we propose a model in which turbulent plasma crosses a conical standing shock in the parsec-scale region of the jet. Here, the high-energy emission in the model is produced by inverse Compton scattering of seed photons supplied by either nonthermal radiation from a Mach disk, thermal emission from hot dust, or (for X-rays) synchrotron radiation from plasma that crosses the standing shock. For the two dates on which we fitted the model SED to the data, the model corresponds very well to the observations at all bands except at X-ray energies, where the spectrum is flatter than observed.Comment: Accepted for publication in Astrophysical Journal. 82 pages, 13 figure

Rapid TeV variability in Blazars as result of Jet-Star Interaction

We propose a new model for the description of ultra-short flares from TeV blazars by compact magnetized condensations (blobs), produced when red giant stars cross the jet close to the central black hole. Our study includes a simple dynamical model for the evolution of the envelope lost by the star in the jet, and its high energy nonthermal emission through different leptonic and hadronic radiation mechanisms. We show that the fragmented envelope of the star can be accelerated to Lorentz factors up to 100 and radiate effectively the available energy in gamma-rays predominantly through proton synchrotron radiation or external inverse Compton scattering of electrons. The model can readily explain the minute-scale TeV flares on top of longer (typical time-scales of days) gamma-ray variability as observed from the blazar PKS 2155-304. In the framework of the proposed scenario, the key parameters of the source are robustly constrained. In the case of proton synchrotron origin of the emission a mass of the central black hole of $M_{\rm BH}\approx 10^8 M_{\odot}$, a total jet power of $L_{\rm j} \approx 2\times 10^{47} \, \rm erg\,s^{-1}$ and a Doppler factor, of the gamma-ray emitting blobs, of $\delta\geq 40$ are required. Whilst for the external inverse Compton model, parameters of $M_{\rm BH}\approx 10^8 M_{\odot}$, $L_{\rm j} \approx 10^{46} \, \rm erg\,s^{-1}$ and the $\delta\geq 150$ are required.Comment: 25 pages, 11 figures, Submitted to Ap

Blazar 3C 454.3 in Outburst and Quiescence During 2005-2007: Two Variable Synchrotron Emission Peaks

We monitored the flaring blazar 3C 454.3 during 2005 June-July with the Spitzer Infrared Spectrograph (IRS: 15 epochs), Infrared Array Camera (IRAC: 12 epochs) and Multiband Imaging Photometer (MIPS: 2 epochs). We also made Spitzer IRS, IRAC, and MIPS observations from 2006 December-2007 January when the source was in a low state, the latter simultaneous with a single Chandra X-ray observation. In addition, we present optical and sub-mm monitoring data. The 2005-2007 period saw 3 major outbursts. We present evidence that the radio-optical SED actually consists of two variable synchrotron peaks, the primary at IR and the secondary at sub-mm wavelengths. The lag between the optical and sub-mm outbursts may indicate that these two peaks arise from two distinct regions along the jet separated by a distance of 0.07-5 pc. The flux at 5-35 microns varied by a factor of 40 and the IR peak varied in frequency from <1E13 Hz to 4E13 Hz between the highest and lowest states in 2005 and 2006, respectively. Variability was well correlated across the mid-IR band, with no measurable lag. Flares that doubled in flux occurred on a time scale of 3 days. The IR SED peak moved to higher frequency as a flare brightened, then returned to lower frequency as it decayed. The fractional variability amplitude increased with frequency, which we attribute to decreasing synchrotron-self absorption optical depth. Mid-IR flares may signal the re-energization of a shock that runs into inhomogeneities along the pre-existing jet or in the external medium. The synchrotron peak frequencies during each major outburst may depend upon both the distance from the jet apex and the physical conditions in the shocks. Variation of the Doppler parameter along a curved or helical jet is another possibility. Frequency variability of the IR synchrotron peak may have important consequences for the interpretation of the blazar sequence, and the presence of a secondary peak may give insight into jet structure.Comment: 38 pages, 15 figures, submitted to ApJS, comments welcom