Most of the light from blazars, active galactic nuclei with jets of magnetized plasma that point nearly along the line of sight, is produced by high-energy particles, up to around 1 TeV. Although the jets are known to be ultimately powered by a supermassive black hole, how the particles are accelerated to such high energies has been an unanswered question. The process must be related to the magnetic field, which can be probed by observations of the polarization of light from the jets. Measurements of the radio to optical polarization—the only range available until now—probe extended regions of the jet containing particles that left the acceleration site days to years earlier^1–3, and hence do not directly explore the acceleration mechanism, as could X-ray measurements. Here we report the detection of X-ray polarization from the blazar Markarian 501 (Mrk 501). We measure an X-ray linear polarization degree Π_X of around 10%, which is a factor of around 2 higher than the value at optical wavelengths, with a polarization angle parallel to the radio jet. This points to a shock front as the source of particle acceleration and also implies that the plasma becomes increasingly turbulent with distance from the shock.

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.01c to 78c, 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.

In an effort to locate the sites of emission at different frequencies and physical processes causing variability in blazar jets, we have obtained high time-resolution observations of BL Lacertae over a wide wavelength range: with the Transiting Exoplanet Survey Satellite (TESS) at 6000–10000 Å with 2 minute cadence; with the Neil Gehrels Swift satellite at optical, UV, and X-ray bands; with the Nuclear Spectroscopic Telescope Array at hard X-ray bands; with the Fermi Large Area Telescope at γ-ray energies; and with the Whole Earth Blazar Telescope for measurement of the optical flux density and polarization. All light curves are correlated, with similar structure on timescales from hours to days. The shortest timescale of variability at optical frequencies observed with TESS is ∼0.5 hr. The most common timescale is 13 ± 1 hr, comparable with the minimum timescale of X-ray variability, 14.5 hr. The multiwavelength variability properties cannot be explained by a change solely in the Doppler factor of the emitting plasma. The polarization behavior implies that there are both ordered and turbulent components to the magnetic field in the jet. Correlation analysis indicates that the X-ray variations lag behind the γ-ray and optical light curves by up to ∼0.4 day. The timescales of variability, cross-frequency lags, and polarization properties can be explained by turbulent plasma that is energized by a shock in the jet and subsequently loses energy to synchrotron and inverse Compton radiation in a magnetic field of strength ∼3 G.

We report the results of decade-long (2008–2018) γ-ray to 1 GHz radio monitoring of the blazar 3C 279, including GASP/WEBT, Fermi and Swift data, as well as polarimetric and spectroscopic data. The X-ray and γ-ray light curves correlate well, with no delay $\gtrsim 3$ h, implying general cospatiality of the emission regions. The γ-ray–optical flux–flux relation changes with activity state, ranging from a linear to a more complex dependence. The behaviour of the Stokes parameters at optical and radio wavelengths, including 43 GHz Very Long Baseline Array images, supports either a predominantly helical magnetic field or motion of the radiating plasma along a spiral path. Apparent speeds of emission knots range from 10 to 37c, with the highest values requiring bulk Lorentz factors close to those needed to explain γ-ray variability on very short time-scales. The Mg ii emission line flux in the ‘blue’ and ‘red’ wings correlates with the optical synchrotron continuum flux density, possibly providing a variable source of seed photons for inverse Compton scattering. In the radio bands, we find progressive delays of the most prominent light-curve maxima with decreasing frequency, as expected from the frequency dependence of the τ = 1 surface of synchrotron self-absorption. The global maximum in the 86 GHz light curve becomes less prominent at lower frequencies, while a local maximum, appearing in 2014, strengthens toward decreasing frequencies, becoming pronounced at ∼5 GHz. These tendencies suggest different Doppler boosting of stratified radio-emitting zones in the jet.

The quasar 3C 454.3 is a blazar known for its rapid and violent outbursts seen across the electromagnetic spectrum. Using γ-ray, X-ray, multiband optical, and very-long-baseline interferometric data we investigate the nature of two such events that occurred in 2013 and 2014 accompanied by strong variations in optical polarization, including a ∼230° electric vector position angle (EVPA) rotation. Our results suggest that a single disturbance was responsible for both flaring events. We interpret the disturbance as a shock propagating down the jet. Under this interpretation the 2013 flare originated most likely due to changes in the viewing angle caused by perhaps a bent or helical trajectory of the shock upstream of the radio core. The 2014 flare and optical polarization behavior are the result of the shock exiting the 43 GHz radio core, suggesting that shock crossings are one of the possible mechanisms for EVPA rotations.