When surrounded by a transparent emission region, black holes are expected to reveal a dark shadow caused by gravitational light bending and photon capture at the event horizon. To image and study this phenomenon, we have assembled the Event Horizon Telescope, a global very long baseline interferometry array observing at a wavelength of 1.3 mm. This allows us to reconstruct event-horizon-scale images of the supermassive black hole candidate in the center of the giant elliptical galaxy M87. We have resolved the central compact radio source as an asymmetric bright emission ring with a diameter of 42±3 μas, which is circular and encompasses a central depression in brightness with a flux ratio 10:1. The emission ring is recovered using different calibration and imaging schemes, with its diameter and width remaining stable over four different observations carried out in different days. Overall, the observed image is consistent with expectations for the shadow of a Kerr black hole as predicted by general relativity. The asymmetry in brightness in the ring can be explained in terms of relativistic beaming of the emission from a plasma rotating close to the speed of light around a black hole. We compare our images to an extensive library of ray-traced general-relativistic magnetohydrodynamic simulations of black holes and derive a central mass of M=(6.5±0.7)×10 9 M e . Our radiowave observations thus provide powerful evidence for the presence of supermassive black holes in centers of galaxies and as the central engines of active galactic nuclei. They also present a new tool to explore gravity in its most extreme limit and on a mass scale that was so far not accessible.

We present measurements of the properties of the central radio source in M87 using Event Horizon Telescope data obtained during the 2017 campaign. We develop and fit geometric crescent models (asymmetric rings with interior brightness depressions) using two independent sampling algorithms that consider distinct representations of the visibility data. We show that the crescent family of models is statistically preferred over other comparably complex geometric models that we explore. We calibrate the geometric model parameters using general relativistic magnetohydrodynamic (GRMHD) models of the emission region and estimate physical properties of the source. We further fit images generated from GRMHD models directly to the data. We compare the derived emission region and black hole parameters from these analyses with those recovered from reconstructed images. There is a remarkable consistency among all methods and data sets. We find that >50% of the total flux at arcsecond scales comes from near the horizon, and that the emission is dramatically suppressed interior to this region by a factor >10, providing direct evidence of the predicted shadow of a black hole. Across all methods, we measure a crescent diameter of 42±3 μas and constrain its fractional width to be <0.5. Associating the crescent feature with the emission surrounding the black hole shadow, we infer an angular gravitational radius of GM/Dc 2 =3.8±0.4 μas. Folding in a distance measurement of -+ 16.8 Mpc 0.7 0.8 gives a black hole mass of =  ´ | | M M 6.5 0.2 0.7 10 stat sys 9. This measurement from lensed emission near the event horizon is consistent with the presence of a central Kerr black hole, as predicted by the general theory of relativity.

The Event Horizon Telescope (EHT) has mapped the central compact radio source of the elliptical galaxy M87 at 1.3 mm with unprecedented angular resolution. Here we consider the physical implications of the asymmetric ring seen in the 2017 EHT data. To this end, we construct a large library of models based on general relativistic magnetohydrodynamic (GRMHD) simulations and synthetic images produced by general relativistic ray tracing. We compare the observed visibilities with this library and confirm that the asymmetric ring is consistent with earlier predictions of strong gravitational lensing of synchrotron emission from a hot plasma orbiting near the black hole event horizon. The ring radius and ring asymmetry depend on black hole mass and spin, respectively, and both are therefore expected to be stable when observed in future EHT campaigns. Overall, the observed image is consistent with expectations for the shadow of a spinning Kerr black hole as predicted by general relativity. If the black hole spin and M87's large scale jet are aligned, then the black hole spin vector is pointed away from Earth. Models in our library of non-spinning black holes are inconsistent with the observations as they do not produce sufficiently powerful jets. At the same time, in those models that produce a sufficiently powerful jet, the latter is powered by extraction of black hole spin energy through mechanisms akin to the Blandford-Znajek process. We briefly consider alternatives to a black hole for the central compact object. Analysis of existing EHT polarization data and data taken simultaneously at other wavelengths will soon enable new tests of the GRMHD models, as will future EHT campaigns at 230 and 345 GHz.

We have examined the compact structure in 250 flat-spectrum extragalactic radio sources using interferometric fringe visibilities obtained with the Very Long Baseline Array ( VLBA) at 15 GHz. With projected baselines out to 440 Mk, we are able to investigate source structure on typical angular scales as small as 0.05 mas. This scale is similar to the resolution of the VLBI Space Observatory Programme data obtained on longer baselines at a lower frequency and with somewhat poorer accuracy. For 171 sources in our sample, more than half of the total flux density seen by the VLBA remains unresolved on the longest baselines. There are 163 sources in our list with a median correlated flux density at 15 GHz in excess of 0.5 Jy on the longest baselines; these will be useful as fringe finders for shortwavelength VLBA observations. The total flux densities recovered in the VLBA images at 15 GHz are generally close to the values measured around the same epoch at the same frequency with the RATAN-600 and University of Michigan Radio Astronomy Observatory telescopes.We have modeled the core of each source with an elliptical Gaussian component. For about 60% of the sources, we have at least one observation in which the core component appears unresolved (generally smaller than 0.05 mas) in one direction, usually transverse to the direction into which the jet extends. BL Lac objects are on average more compact than quasars, while active galaxies are on average less compact. Also, in an active galaxy the submilliarcsecond core component tends to be less dominant. Intraday variability (IDV) sources typically have a more compact, more core-dominated structure on sub-milliarcsecond scales than non-IDV sources, and sources with a greater amplitude of intraday variations tend to have a greater unresolved VLBA flux density. The objects known to be GeV gamma-ray-loud appear to have a more compact VLBA structure than the other sources in our sample. This suggests that the mechanisms for the production of gamma-ray emission and for the generation of compact radio synchrotron-emitting features are related. The brightness temperature estimates and lower limits for the cores in our sample typically range between 10 11 and 10 13 K, but they extend up to 5 ; 10 13 K, apparently in excess of the equipartition brightness temperature or the inverse Compton limit for stationary synchrotron sources. The largest component speeds are observed in radio sources with high observed brightness temperatures, as would be expected from relativistic beaming. Longer baselines, which can be obtained by space VLBI observations, will be needed to resolve the most compact high brightness temperature regions in these sources.

Aims. We have investigated a frequency-dependent shift in the absolute position of the optically thick apparent origin of parsec-scale jets ("core shift" effect) to probe physical conditions in ultra-compact relativistic outflows in active galactic nuclei. , and 88 μas, respectively, compared to the typical measured errors of 50, 51, 35 μas. The effect occurs predominantly along the jet direction, with departures smaller than 45 • from the median jet position angle in over 80% of the cases. Despite the moderate ratio of the observed frequencies (<2), core shifts significantly different from zero (>2σ) are detected for about 55% of the sources. These shifts are even better aligned with the jet direction, deviating from the latter by less than 30 • in over 90% of the cases. There is an indication that the core shift decreases with increasing redshift. Magnetic fields in the jet at a distance of 1 parsec from the central black hole, calculated from the obtained core shifts, are found to be systematically stronger in quasars (median B 1 ≈ 0.9 G) than those in BL Lacs (median B 1 ≈ 0.4 G). We also constrained the absolute distance of the core from the apex of the jet at 15 GHz as well as the magnetic field strength in the 15 GHz core region.

Context. Very long baseline interferometry (VLBI) imaging of radio emission from extragalactic jets provides a unique probe of physical mechanisms governing the launching, acceleration, and collimation of relativistic outflows. Aims. VLBI imaging of the jet in the nearby active galaxy M 87 enables morphological and kinematic studies to be done on linear scales down to ∼ 100 Schwarzschild radii (R s ). Methods. The two-dimensional structure and kinematics of the jet in M 87 (NGC 4486) have been studied by applying the Waveletbased Image Segmentation and Evaluation (WISE) method to 11 images obtained from multi-epoch Very Long Baseline Array (VLBA) observations made in January-August 2007 at 43 GHz (λ = 7 mm).Results. The WISE analysis recovers a detailed two-dimensional velocity field in the jet in M 87 at sub-parsec scales. The observed evolution of the flow velocity with distance from the jet base can be explained in the framework of MHD jet acceleration and Poynting flux conversion. A linear acceleration regime is observed up to z obs ∼ 2 mas. The acceleration is reduced at larger scales, which is consistent with saturation of Poynting flux conversion. Stacked cross correlation analysis of the images reveals a pronounced stratification of the flow. The flow consists of a slow, mildly relativistic layer (moving at β ∼ 0.5 c), associated either with instability pattern speed or an outer wind, and a fast, accelerating stream line (with β ∼ 0.92, corresponding to a bulk Lorentz factor γ ∼ 2.5). A systematic difference of the apparent speeds in the northern and southern limbs of the jet is detected, providing evidence for jet rotation. The angular velocity of the magnetic field line associated with this rotation suggests that the jet in M87 is launched in the inner part of the disk, at a distance r 0 ∼ 5 R s from the central engine.Conclusions. The combined results of the analysis imply that MHD acceleration and conversion of Poynting flux to kinetic energy play the dominant roles in collimation and acceleration of the flow in M 87.

We present the first Event Horizon Telescope (EHT) images of M87, using observations from April 2017 at 1.3 mm wavelength. These images show a prominent ring with a diameter of ∼40 μas, consistent with the size and shape of the lensed photon orbit encircling the "shadow" of a supermassive black hole. The ring is persistent across four observing nights and shows enhanced brightness in the south. To assess the reliability of these results, we implemented a two-stage imaging procedure. In the first stage, four teams, each blind to the others' work, produced images of M87 using both an established method (CLEAN) and a newer technique (regularized maximum likelihood). This stage allowed us to avoid shared human bias and to assess common features among independent reconstructions. In the second stage, we reconstructed synthetic data from a large survey of imaging parameters and then compared the results with the corresponding ground truth images. This stage allowed us to select parameters objectively to use when reconstructing images of M87. Across all tests in both stages, the ring diameter and asymmetry remained stable, insensitive to the choice of imaging technique. We describe the EHT imaging procedures, the primary image features in M87, and the dependence of these features on imaging assumptions.

Context. The apparent position of the "core" in a parsec-scale radio jet (a compact, bright emitting region at the narrow end of the jet) depends on the observing frequency, owing to synchrotron self-absorption and external absorption. While providing a tool probing physical conditions in the vicinity of the core, this dependency poses problems for astrometric studies using compact radio sources. Aims. We investigated the frequency-dependent shift in the positions of the cores (core shift) observed with very long baseline interferometry (VLBI) in parsec-scale jets. We discuss related physics, as well as its effect on radio astrometry and the connection between radio and optical positions of astrometric reference objects. Methods. We searched for the core shift in a sample of 277 radio sources imaged at 2.3 GHz (13 cm) and 8.6 GHz (4 cm) frequency bands using VLBI observations made in 2002 and 2003. The core shift was measured by referencing the core position to optically thin jet features whose positions are not expected to change with frequency. Results. We present here results for 29 selected active galactic nuclei (AGN) with bright distinct VLBI jet features that can be used in differential measurements and that allow robust measurements of the shift to be made. In these AGN, the magnitude of the measured core shift between 2.3 and 8.6 GHz reaches 1.4 mas, with a median value for the sample of 0.44 mas. Nuclear flares result in temporal variability of the shift.Conclusions. An average shift between the radio (4 cm) and optical (6000 Å) bands is estimated to be approximately 0.1 mas, and it should be taken into account in order to provide the required accuracy of the radio-optical reference frame connection. This can be accomplished with multi-frequency VLBI measurements yielding estimates of the core shift in the sources used for the radio reference frame and radio-optical position alignment.