Modeling the time-dependent polarization of blazars

Linear polarization is an extremely valuable observational tool for probing the dynamic physical conditions of blazar jets. Some patterns are seen in the data, suggestive of order that can be explained by shock waves and helical magnetic field components. However, much disorder is apparent, which implies that turbulence plays a major role as well, especially in the fluctuations of flux and polarization, and perhaps particle acceleration. Here, we present some actual flux and polarization versus time data, plus simulations of model jets. We analyze the output of the simulations in a manner that can be compared with observational data. The results suggest that the ratio of turbulent to ordered magnetic fields varies with time.AST-1615796 - National Science Foundation; NASA; NNX14AQ58G; NNX15AR45

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

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

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

LOFAR observations of 4C+19.44. On the discovery of low frequency spectral curvature in relativistic jet knots

We present the first LOFAR observations of the radio jet in the quasar 4C+19.44 (a.k.a. PKS 1354+19) obtained with the long baselines. The achieved resolution is very well matched to that of archival Jansky Very Large Array (JVLA) observations at higher radio frequencies as well as the archival X-ray images obtained with {\it Chandra}. We found that, for several knots along the jet, the radio flux densities measured at hundreds of MHz lie well below the values estimated by extrapolating the GHz spectra. This clearly indicates the presence of spectral curvature. Radio spectral curvature has been already observed in different source classes and/or extended radio structures and it has been often interpreted as due to intrinsic processes, as a curved particle energy distribution, rather than absorption mechanisms ({ Razin-Tsytovich} effect, free-free or synchrotron self absorption to name a few). Here we discuss our results according to the scenario where particles undergo stochastic acceleration mechanisms also in quasar jet knots.Comment: 13 pages, 4 tables, 4 figures, pre-proof version, published on the Astrophysical Journal (Harris, et al. 2019 ApJ, 873, 21

Behaviour of the Blazar CTA 102 during two giant outbursts

Blazar CTA 102 underwent exceptional optical and high-energy outbursts in 2012 and 2016-2017. We analyze its behaviour during these events, focusing on polarimetry as a tool that allows us to trace changes in the physical conditions and geometric configuration of the emission source close to the central black hole. We also use Fermi gamma-ray data in conjunction with optical photometry in an effort to localize the origin of the outbursts.AST-1615796 - Boston Universit

MOJAVE: Monitoring of Jets in AGN with VLBA Experiments. VII. Blazar Jet Acceleration

We discuss acceleration measurements for a large sample of extragalactic radio jets from the MOJAVE program which studies the parsec-scale jet structure and kinematics of a complete, flux-density-limited sample of Active Galactic Nuclei (AGN). Accelerations are measured from the apparent motion of individual jet features or "components" which may represent patterns in the jet flow. We find that significant accelerations are common both parallel and perpendicular to the observed component velocities. Parallel accelerations, representing changes in apparent speed, are generally larger than perpendicular acceleration that represent changes in apparent direction. The trend for larger parallel accelerations indicates that a significant fraction of these changes in apparent speed are due to changes in intrinsic speed of the component rather than changes in direction to the line of sight. We find an overall tendency for components with increasing apparent speed to be closer to the base of their jets than components with decreasing apparent speed. This suggests a link between the observed pattern motions and the underlying flow which, in some cases, may increase in speed close to the base and decrease in speed further out; however, common hydro-dynamical processes for propagating shocks may also play a role. About half of the components show "non-radial" motion, or a misalignment between the component's structural position angle and its velocity direction, and these misalignments generally better align the component motion with the downstream emission. Perpendicular accelerations are closely linked with non-radial motion. When observed together, perpendicular accelerations are usually in the correct direction to have caused the observed misalignment.Comment: 17 pages, 11 figures, 1 table, accepted by the Astrophysical Journa