The nearby radio galaxy M87 is a prime target for studying black hole accretion and jet formation^1,2. Event Horizon Telescope observations of M87 in 2017, at a wavelength of 1.3 mm, revealed a ring-like structure, which was interpreted as gravitationally lensed emission around a central black hole³. Here we report images of M87 obtained in 2018, at a wavelength of 3.5 mm, showing that the compact radio core is spatially resolved. High-resolution imaging shows a ring-like structure of $${8.4}_{-1.1}^{+0.5}$$8.4 − 1.1 + 0.5 Schwarzschild radii in diameter, approximately 50% larger than that seen at 1.3 mm. The outer edge at 3.5 mm is also larger than that at 1.3 mm. This larger and thicker ring indicates a substantial contribution from the accretion flow with absorption effects, in addition to the gravitationally lensed ring-like emission. The images show that the edge-brightened jet connects to the accretion flow of the black hole. Close to the black hole, the emission profile of the jet-launching region is wider than the expected profile of a black-hole-driven jet, suggesting the possible presence of a wind associated with the accretion flow.

Context. The nearby elliptical galaxy M87 contains one of only two supermassive black holes whose emission surrounding the event horizon has been imaged by the Event Horizon Telescope (EHT). In 2018, more than two dozen multi-wavelength (MWL) facilities (from radio to γ-ray energies) took part in the second M87 EHT campaign. Aims. The goal of this extensive MWL campaign was to better understand the physics of the accreting black hole M87*, the relationship between the inflow and inner jets, and the high-energy particle acceleration. Understanding the complex astrophysics is also a necessary first step towards performing further tests of general relativity. Methods. The MWL campaign took place in April 2018, overlapping with the EHT M87* observations. We present a new, contemporaneous spectral energy distribution (SED) ranging from radio to very high-energy (VHE) γ-rays as well as details of the individual observations and light curves. We also conducted phenomenological modelling to investigate the basic source properties. Results. We present the first VHE γ-ray flare from M87 detected since 2010. The flux above 350 GeV more than doubled within a period of ≈36 hours. We find that the X-ray flux is enhanced by about a factor of two compared to 2017, while the radio and millimetre core fluxes are consistent between 2017 and 2018. We detect evidence for a monotonically increasing jet position angle that corresponds to variations in the bright spot of the EHT image. Conclusions. Our results show the value of continued MWL monitoring together with precision imaging for addressing the origins of high-energy particle acceleration. While we cannot currently pinpoint the precise location where such acceleration takes place, the new VHE γ-ray flare already presents a challenge to simple one-zone leptonic emission model approaches, and it emphasises the need for combined image and spectral modelling.

The Event Horizon Telescope (EHT) observation of M87^∗ in 2018 has revealed a ring with a diameter that is consistent with the 2017 observation. The brightest part of the ring is shifted to the southwest from the southeast. In this paper, we provide theoretical interpretations for the multi-epoch EHT observations for M87^∗ by comparing a new general relativistic magnetohydrodynamics model image library with the EHT observations for M87^∗ in both 2017 and 2018. The model images include aligned and tilted accretion with parameterized thermal and nonthermal synchrotron emission properties. The 2018 observation again shows that the spin vector of the M87^∗ supermassive black hole is pointed away from Earth. A shift of the brightest part of the ring during the multi-epoch observations can naturally be explained by the turbulent nature of black hole accretion, which is supported by the fact that the more turbulent retrograde models can explain the multi-epoch observations better than the prograde models. The EHT data are inconsistent with the tilted models in our model image library. Assuming that the black hole spin axis and its large-scale jet direction are roughly aligned, we expect the brightest part of the ring to be most commonly observed 90 deg clockwise from the forward jet. This prediction can be statistically tested through future observations.

Context. Very long baseline interferometry (VLBI) observations at 86 GHz (wavelength, λ = 3 mm) reach a resolution of about 50 μas, probing the collimation and acceleration regions of relativistic outflows in active galactic nuclei (AGN). The physical conditions in these regions can be studied by performing 86 GHz VLBI surveys of representative samples of compact extragalactic radio sources. Aims. To extend the statistical studies of compact extragalactic jets, a large global 86 GHz VLBI survey of 162 compact radio sources was conducted in 2010–2011 using the Global Millimeter VLBI Array (GMVA). Methods. The survey observations were made in a snapshot mode, with up to five scans per target spread over a range of hour angles in order to optimize the visibility coverage. The survey data attained a typical baseline sensitivity of 0.1 Jy and a typical image sensitivity of 5 mJy beam⁻¹, providing successful detections and images for all of the survey targets. For 138 objects, the survey provides the first ever VLBI images made at 86 GHz. Gaussian model fitting of the visibility data was applied to represent the structure of the observed sources and to estimate the flux densities and sizes of distinct emitting regions (components) in their jets. These estimates were used for calculating the brightness temperature (T_b) at the jet base (core) and in one or more moving regions (jet components) downstream from the core. These model-fit-based estimates of T_b were compared to the estimates of brightness temperature limits made directly from the visibility data, demonstrating a good agreement between the two methods. Results. The apparent brightness temperature estimates for the jet cores in our sample range from 2.5 × 10⁹ K to 1.3 × 10¹² K, with the mean value of 1.8 × 10¹¹ K. The apparent brightness temperature estimates for the inner jet components in our sample range from 7.0 × 10⁷ K to 4.0 × 10¹¹ K. A simple population model with a single intrinsic value of brightness temperature, T₀, is applied to reproduce the observed distribution. It yields T₀ = (3.77_−0.14^+0.10) × 10¹¹ K for the jet cores, implying that the inverse Compton losses dominate the emission. In the nearest jet components, T₀ = (1.42_−0.19^+0.16) × 10¹¹ K is found, which is slightly higher than the equipartition limit of ∼5 × 10¹⁰ K expected for these jet regions. For objects with sufficient structural detail detected, the adiabatic energy losses are shown to dominate the observed changes of brightness temperature along the jet.

We report three epochs of polarized images of M87* at 230 GHz using data from the Event Horizon Telescope (EHT) taken in 2017, 2018, and 2021. The baseline coverage of the 2021 observations is significantly improved through the addition of two new EHT stations: the 12 m Kitt Peak Telescope and the Northern Extended Millimetre Array (NOEMA). All observations result in images dominated by a bright, asymmetric ring with a persistent diameter of 43.9 ± 0.6 μas, consistent with expectations for lensed synchrotron emission encircling the apparent shadow of a supermassive black hole. We find that the total intensity and linear polarization of M87* vary significantly across the three epochs. Specifically, the azimuthal brightness distribution of the total intensity images varies from year to year, as expected for a stochastic accretion flow. However, despite a gamma-ray flare erupting in M87 quasi-contemporaneously to the 2018 observations, the 2018 and 2021 images look remarkably similar. The resolved linear polarization fractions in 2018 and 2021 peak at ∼5%, compared to ∼15% in 2017. The spiral polarization pattern on the ring also varies from year to year, including a change in the electric vector position angle helicity in 2021 that could reflect changes in the magnetized accretion flow or an external Faraday screen. The improved 2021 coverage also provides the first EHT constraints on jet emission outside the ring, on scales of ≲1 mas. Overall, these observations provide strong proof of the reliability of the EHT images and probe the dynamic properties of the horizon-scale accretion flow surrounding M87*.