We analyze the parsec-scale jet kinematics from 2007 June to 2013 January of a sample of γ-ray bright blazars monitored roughly monthly with the Very Long Baseline Array at 43 GHz. In a total of 1929 ...images, we measure apparent speeds of 252 emission knots in 21 quasars, 12 BL Lacertae objects (BLLacs), and 3 radio galaxies, ranging from 0.02c to 78c; 21% of the knots are quasi-stationary. Approximately one-third of the moving knots execute non-ballistic motions, with the quasars exhibiting acceleration along the jet within 5 pc (projected) of the core, and knots in BLLacs tending to decelerate near the core. Using the apparent speeds of the components and the timescales of variability from their light curves, we derive the physical parameters of 120 superluminal knots, including variability Doppler factors, Lorentz factors, and viewing angles. We estimate the half-opening angle of each jet based on the projected opening angle and scatter of intrinsic viewing angles of knots. We determine characteristic values of the physical parameters for each jet and active galactic nucleus class based on the range of values obtained for individual features. We calculate the intrinsic brightness temperatures of the cores, , at all epochs, finding that the radio galaxies usually maintain equipartition conditions in the cores, while ∼30% of measurements in the quasars and BLLacs deviate from equipartition values by a factor >10. This probably occurs during transient events connected with active states. In the Appendix, we briefly describe the behavior of each blazar during the period analyzed.
Abstract
We analyze the parsec-scale jet kinematics from 2007 June to 2018 December of a sample of
γ
-ray bright blazars monitored roughly monthly with the Very Long Baseline Array (VLBA) at 43 GHz ...under the VLBA-BU-BLAZAR program. We implement a novel piecewise linear fitting method to derive the kinematics of 521 distinct emission knots from a total of 3705 total intensity images in 22 quasars, 13 BL Lacertae objects, and 3 radio galaxies. Apparent speeds of these components range from 0.01
c
to 78
c
, and 18.6% of knots (other than the “core”) are quasi-stationary. One-fifth of moving knots exhibit nonballistic motion, with acceleration along the jet within 5 pc of the core (projected) and deceleration farther out. These accelerations occur mainly at locations coincident with quasi-stationary features. We calculate the physical parameters of 273 knots with statistically significant motion, including their Doppler factors, Lorentz factors, and viewing angles. We determine the typical values of these parameters for each jet and the average for each subclass of active galactic nuclei. We investigate the variability of the position angle of each jet over the 10 yr of monitoring. The fluctuations in position of the quasi-stationary components in radio galaxies tend to be parallel to the jet, while no directional preference is seen in the components of quasars and BL Lacertae objects. We find a connection between
γ
-ray states of blazars and their parsec-scale jet properties, with blazars with brighter 43 GHz cores typically reaching higher
γ
-ray maxima during flares.
Time-variable polarization is an extremely valuable observational tool to probe the dynamical physical conditions of blazar jets. Since 2008, we have been monitoring the flux and linear polarization ...of a sample of gamma-ray bright blazars at optical frequencies. Some of the observations were performed on nightly or intra-night time-scales in four optical bands, providing information on the frequency and time dependence of the polarization. The observed behavior is similar to that found in simulations of turbulent plasma in a relativistic jet that contains a standing shock and/or a helical background magnetic field. Similar simulations predict the characteristics of X-ray synchrotron polarization of blazars that will be measured in the future by the Imaging X-ray Polarimetry Explorer (IXPE).
ABSTRACT Blazars exhibit flares across the electromagnetic spectrum. Many γ-ray flares are highly correlated with flares detected at optical wavelengths; however, a small subset appears to occur in ...isolation, with little or no variability detected at longer wavelengths. These "orphan" γ-ray flares challenge current models of blazar variability, most of which are unable to reproduce this type of behavior. We present numerical calculations of the time-variable emission of a blazar based on a proposal by Marscher et al. to explain such events. In this model, a plasmoid ("blob") propagates relativistically along the spine of a blazar jet and passes through a synchrotron-emitting ring of electrons representing a shocked portion of the jet sheath. This ring supplies a source of seed photons that are inverse-Compton scattered by the electrons in the moving blob. The model includes the effects of radiative cooling, a spatially varying magnetic field, and acceleration of the blob's bulk velocity. Synthetic light curves produced by our model are compared to the observed light curves from an orphan flare that was coincident with the passage of a superluminal knot through the inner jet of the blazar PKS 1510-089. In addition, we present Very Long Baseline Array polarimetric observations that point to the existence of a jet sheath in PKS 1510-089, thus providing further observational support for the plausibility of our model. An estimate of the bolometric luminosity of the sheath within PKS 1510-089 is made, yielding . This indicates that the sheath within PKS 1510-089 is potentially a very important source of seed photons.
Blazars exhibit flares across the entire electromagnetic spectrum. Many γ-ray flares are highly correlated with flares detected at longer wavelengths; however, a small subset appears to occur in ...isolation, with little or no correlated variability at longer wavelengths. These "orphan" γ-ray flares challenge current models of blazar variability, most of which are unable to reproduce this type of behavior. MacDonald et al. have developed the Ring of Fire model to explain the origin of orphan γ-ray flares from within blazar jets. In this model, electrons contained within a blob of plasma moving relativistically along the spine of the jet inverse-Compton scatter synchrotron photons emanating off of a ring of shocked sheath plasma that enshrouds the jet spine. As the blob propagates through the ring, the scattering of the ring photons by the blob electrons creates an orphan γ-ray flare. This model was successfully applied to modeling a prominent orphan γ-ray flare observed in the blazar PKS 1510−089. To further support the plausibility of this model, MacDonald et al. presented a stacked radio map of PKS 1510−089 containing the polarimetric signature of a sheath of plasma surrounding the spine of the jet. In this paper, we extend our modeling and stacking techniques to a larger sample of blazars: 3C 273, 4C 71.01, 3C 279, 1055+018, CTA 102, and 3C 345, the majority of which have exhibited orphan γ-ray flares. We find that the model can successfully reproduce these flares, while our stacked maps reveal the existence of jet sheaths within these blazars.
We present the first polarimetric space very long baseline interferometry (VLBI) imaging observations at 22 GHz. BL Lacertae was observed in 2013 November 10 with the RadioAstron space VLBI mission, ...including a ground array of 15 radio telescopes. The instrumental polarization of the space radio telescope is found to be less than 9%, demonstrating the polarimetric imaging capabilities of RadioAstron at 22 GHz. Ground-space fringes were obtained up to a projected baseline distance of 7.9 Earth diameters in length, allowing us to image the jet in BL Lacertae with a maximum angular resolution of 21 mu as, the highest achieved to date. We find evidence for emission upstream of the radio core, which may correspond to a recollimation shock at about 40 mu as from the jet apex, in a pattern that includes other recollimation shocks at approximately 100 and 250 mu as from the jet apex. Polarized emission is detected in two components within the innermost 0.5 mas from the core, as well as in some knots 3 mas downstream. Faraday rotation analysis, obtained from combining RadioAstron 22 GHz and ground-based 15 and 43 GHz images, shows a gradient in rotation measure and Faraday-corrected polarization vector as a function of position angle with respect to the core, suggesting that the jet in BL Lacertae is threaded by a helical magnetic field. The intrinsic de-boosted brightness temperature in the unresolved core exceeds 3 x 10 super(12) K, suggesting, at the very least, departure from equipartition of energy between the magnetic field and radiating particles.
Blazars whose synchrotron spectral energy distribution (SED) peaks at X-ray energies need to accelerate electrons to energies in the >100 GeV range in relativistic plasma jets at distances of parsecs ...from the central engine. Compton scattering by the same electrons can explain high luminosities at very high photon energies (>100 GeV) from the same objects. Turbulence combined with a standing conical shock can accomplish this. Such a scenario can also qualitatively explain the level and variability of linear polarization observed at optical frequencies in these objects. Multi-wavelength polarization measurements, including those at X-ray energies by the Imaging X-ray Polarimetry Explorer (IXPE), find that the degree of polarization is several times higher at X-ray than at optical wavelengths, in general agreement with the turbulence-plus-shock picture. Some detailed properties of the observed polarization can be naturally explained by this scenario, while others pose challenges that may require modifications to the model.
A possible source of Delta *g-ray photons observed from the jets of blazars is inverse Compton scattering by relativistic electrons of infrared seed photons from a hot, dusty torus in the nucleus. We ...use observations from the Spitzer Space Telescope to search for signatures of such dust in the infrared spectra of four Delta *g-ray bright blazars, the quasars 4C 21.35, CTA102, and PKS 1510--089, and the BL Lacertae object ON231. The spectral energy distribution (SED) of 4C 21.35 contains a prominent infrared excess indicative of dust emission. After subtracting a non-thermal component with a power-law spectrum, we fit a dust model to the residual SED. The model consists of a blackbody with temperature ~1200 K, plus a much weaker optically thin component at ~660 K. The total luminosity of the thermal dust emission is 7.9 ? 0.2 X 1045 erg s--1. If the dust lies in an equatorial torus, the density of infrared photons from the torus is sufficient to explain the Delta *g-ray flux from 4C 21.35 as long as the scattering occurs within a few parsecs of the central engine. We also report a tentative detection of dust in the quasar CTA102, in which the luminosity of the infrared excess is 7 ? 2 X 1045 erg s--1. However, in CTA102 the far-infrared spectra are too noisy to detect the 10 Delta *mm silicate feature. Upper limits to the luminosity from thermal emission from dust in PKS 1510--089, and ON231, are 2.3 X 1045, and 6.6 X 1043 erg s--1, respectively. These upper limits do not rule out the possibility of inverse Compton upscattering of infrared photons to Delta *g-ray energies in these two sources. The estimated covering factor of the hot dust in 4C 21.35, 22%, is similar to that of non-blazar quasars; however, 4C 21.35 is deficient in cooler dust.
Context. Investigating the magnetic field structure in the innermost regions of relativistic jets is fundamental to understanding the crucial physical processes giving rise to jet formation, as well ...as to their extraordinary radiation output up to γ-ray energies. Aims. We study the magnetic field structure of the quasar CTA 102 with 3 and 7 mm VLBI polarimetric observations, reaching an unprecedented resolution (∼50 μas). We also investigate the variability and physical processes occurring in the source during the observing period, which coincides with a very active state of the source over the entire electromagnetic spectrum. Methods. We perform the Faraday rotation analysis using 3 and 7 mm data and we compare the obtained rotation measure (RM) map with the polarization evolution in 7 mm VLBA images. We study the kinematics and variability at 7 mm and infer the physical parameters associated with variability. From the analysis of γ-ray and X-ray data, we compute a minimum Doppler factor value required to explain the observed high-energy emission. Results. Faraday rotation analysis shows a gradient in RM with a maximum value of ∼6 × 104 rad m−2 and intrinsic electric vector position angles (EVPAs) oriented around the centroid of the core, suggesting the presence of large-scale helical magnetic fields. Such a magnetic field structure is also visible in 7 mm images when a new superluminal component is crossing the core region. The 7 mm EVPA orientation is different when the component is exiting the core or crossing a stationary feature at ∼0.1 mas. The interaction between the superluminal component and a recollimation shock at ∼0.1 mas could have triggered the multi-wavelength flares. The variability Doppler factor associated with such an interaction is large enough to explain the high-energy emission and the remarkable optical flare occurred very close in time.
ABSTRACT BL Lac, the eponymous blazar, flared to historically high levels at millimeter, infrared, X-ray, and gamma-ray wavelengths in 2012. We present observations made with Herschel, Swift, NuSTAR, ...Fermi, the Submillimeter Array, CARMA, and the VLBA in 2012-2013, including three months with nearly daily sampling at several wavebands. We have also conducted an intensive campaign of 30 hr with every-orbit observations by Swift and NuSTAR, accompanied by Herschel, and Fermi observations. The source was highly variable at all bands. Time lags, correlations between bands, and the changing shapes of the spectral energy distributions can be explained by synchrotron radiation and inverse Compton emission from nonthermal seed photons originating from within the jet. The passage of four new superluminal very long baseline interferometry knots through the core and two stationary knots about 4 pc downstream accompanied the high flaring in 2012-2013. The seed photons for inverse Compton scattering may arise from the stationary knots and from a Mach disk near the core where relatively slow-moving plasma generates intense nonthermal radiation. The 95 spectral energy distributions obtained on consecutive days form the most densely sampled, broad wavelength coverage for any blazar. The observed spectral energy distributions and multi-waveband light curves are similar to simulated spectral energy distributions and light curves generated with a model in which turbulent plasma crosses a conical shock with a Mach disk.