In 2017, the Event Horizon Telescope (EHT) Collaboration succeeded in capturing the first direct image of the center of the M87 galaxy. The asymmetric ring morphology and size are consistent with theoretical expectations for a weakly accreting supermassive black hole of mass ∼6.5 × 10⁹ M _⊙. The EHTC also partnered with several international facilities in space and on the ground, to arrange an extensive, quasi-simultaneous multi-wavelength campaign. This Letter presents the results and analysis of this campaign, as well as the multi-wavelength data as a legacy data repository. We captured M87 in a historically low state, and the core flux dominates over HST-1 at high energies, making it possible to combine core flux constraints with the more spatially precise very long baseline interferometry data. We present the most complete simultaneous multi-wavelength spectrum of the active nucleus to date, and discuss the complexity and caveats of combining data from different spatial scales into one broadband spectrum. We apply two heuristic, isotropic leptonic single-zone models to provide insight into the basic source properties, but conclude that a structured jet is necessary to explain M87’s spectrum. We can exclude that the simultaneous γ-ray emission is produced via inverse Compton emission in the same region producing the EHT mm-band emission, and further conclude that the γ-rays can only be produced in the inner jets (inward of HST-1) if there are strongly particle-dominated regions. Direct synchrotron emission from accelerated protons and secondaries cannot yet be excluded.

Context. Over the past few years, on several occasions, large, continuous rotations of the electric vector position angle (EVPA) of linearly polarized optical emission from blazars have been reported. These events are often coincident with high energy γ-ray flares and they have attracted considerable attention, since they could allow us to probe the magnetic field structure in the γ-ray emitting region of the jet. The flat-spectrum radio quasar 3C 279 is one of the most prominent examples showing this behaviour.Aims. Our goal is to study the observed EVPA rotations and to distinguish between a stochastic and a deterministic origin of the polarization variability.Methods. We have combined multiple data sets of R-band photometry and optical polarimetry measurements of 3C 279, yielding exceptionally well-sampled flux density and polarization curves that cover a period of 2008–2012. Several large EVPA rotations are identified in the data. We introduce a quantitative measure for the EVPA curve smoothness, which is then used to test a set of simple random walk polarization variability models against the data.Results. 3C 279 shows different polarization variation characteristics during an optical low-flux state and a flaring state. The polarization variation during the flaring state, especially the smooth ~360° rotation of the EVPA in mid-2011, is not consistent with the tested stochastic processes.Conclusions. We conclude that, during the two different optical flux states, two different processes govern polarization variation, which is possibly a stochastic process during the low-brightness state and a deterministic process during the flaring activity.

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.

After several years of quiescence, the blazar CTA 102 underwent an exceptional outburst in 2012 September–October. The flare was tracked from γ-ray to near-infrared (NIR) frequencies, including Fermi and Swift data as well as photometric and polarimetric data from several observatories. An intensive Glast-Agile support programme of the Whole Earth Blazar Telescope (GASP–WEBT) collaboration campaign in optical and NIR bands, with an addition of previously unpublished archival data and extension through fall 2015, allows comparison of this outburst with the previous activity period of this blazar in 2004–2005. We find remarkable similarity between the optical and γ-ray behaviour of CTA 102 during the outburst, with a time lag between the two light curves of ≈1 h, indicative of cospatiality of the optical and γ-ray emission regions. The relation between the γ-ray and optical fluxes is consistent with the synchrotron self-Compton (SSC) mechanism, with a quadratic dependence of the SSC γ-ray flux on the synchrotron optical flux evident in the post-outburst stage. However, the γ-ray/optical relationship is linear during the outburst; we attribute this to changes in the Doppler factor. A strong harder-when-brighter spectral dependence is seen both the in γ-ray and optical non-thermal emission. This hardening can be explained by convexity of the UV–NIR spectrum that moves to higher frequencies owing to an increased Doppler shift as the viewing angle decreases during the outburst stage. The overall pattern of Stokes parameter variations agrees with a model of a radiating blob or shock wave that moves along a helical path down the jet.