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 a measurement of the Hubble constant made using geometric distance measurements to megamaser-hosting galaxies. We have applied an improved approach for fitting maser data and obtained better distance estimates for four galaxies previously published by the Megamaser Cosmology Project: UGC 3789, NGC 6264, NGC 6323, and NGC 5765b. Combining these updated distance measurements with those for the maser galaxies CGCG 074-064 and NGC 4258, and assuming a fixed velocity uncertainty of 250 km s −1 associated with peculiar motions, we constrain the Hubble constant to be H 0 = 73.9±3.0 km s −1 Mpc −1 independent of distance ladders and the cosmic microwave background. This best value relies solely on maser-based distance and velocity measurements, and it does not use any peculiar velocity corrections. Different approaches for correcting peculiar velocities do not modify H 0 by more than ±1σ, with the full range of best-fit Hubble constant values spanning 71.8-76.9 km s −1 Mpc −1 . We corroborate prior indications that the local value of H 0 exceeds the early-Universe value, with a confidence level varying from 95-99% for different treatments of the peculiar velocities.

The Event Horizon Telescope (EHT) is a very long baseline interferometry (VLBI) array that comprises millimeter-and submillimeter-wavelength telescopes separated by distances comparable to the diameter of the Earth. At a nominal operating wavelength of ∼1.3 mm, EHT angular resolution (λ/D) is ∼25 μas, which is sufficient to resolve nearby supermassive black hole candidates on spatial and temporal scales that correspond to their event horizons. With this capability, the EHT scientific goals are to probe general relativistic effects in the strong-field regime and to study accretion and relativistic jet formation near the black hole boundary. In this Letter we describe the system design of the EHT, detail the technology and instrumentation that enable observations, and provide measures of its performance. Meeting the EHT science objectives has required several key developments that have facilitated the robust extension of the VLBI technique to EHT observing wavelengths and the production of instrumentation that can be deployed on a heterogeneous array of existing telescopes and facilities. To meet sensitivity requirements, high-bandwidth digital systems were developed that process data at rates of 64gigabit s −1 , exceeding those of currently operating cm-wavelength VLBI arrays by more than an order of magnitude. Associated improvements include the development of phasing systems at array facilities, new receiver installation at several sites, and the deployment of hydrogen maser frequency standards to ensure coherent data capture across the array. These efforts led to the coordination and execution of the first Global EHT observations in 2017 April, and to event-horizon-scale imaging of the supermassive black hole candidate in M87.

The star S2 orbiting the compact radio source Sgr A* is a precision probe of the gravitational field around the closest massive black hole (candidate). Over the last 2.7 decades we have monitored the star’s radial velocity and motion on the sky, mainly with the SINFONI and NACO adaptive optics (AO) instruments on the ESO VLT, and since 2017, with the four-telescope interferometric beam combiner instrument GRAVITY. In this Letter we report the first detection of the General Relativity (GR) Schwarzschild Precession (SP) in S2’s orbit. Owing to its highly elliptical orbit (e = 0.88), S2’s SP is mainly a kink between the pre-and post-pericentre directions of motion ≈±1 year around pericentre passage, relative to the corresponding Kepler orbit. The superb 2017−2019 astrometry of GRAVITY defines the pericentre passage and outgoing direction. The incoming direction is anchored by 118 NACO-AO measurements of S2’s position in the infrared reference frame, with an additional 75 direct measurements of the S2-Sgr A* separation during bright states (“flares”) of Sgr A*. Our 14-parameter model fits for the distance, central mass, the position and motion of the reference frame of the AO astrometry relative to the mass, the six parameters of the orbit, as well as a dimensionless parameter fSP for the SP (fSP = 0 for Newton and 1 for GR). From data up to the end of 2019 we robustly detect the SP of S2, δϕ ≈ 12′ per orbital period. From posterior fitting and MCMC Bayesian analysis with different weighting schemes and bootstrapping we find fSP = 1.10 ± 0.19. The S2 data are fully consistent with GR. Any extended mass inside S2’s orbit cannot exceed ≈0.1% of the central mass. Any compact third mass inside the central arcsecond must be less than about 1000 M⊙.

The highly elliptical, 16-year-period orbit of the star S2 around the massive black hole candidate Sgr A✻ is a sensitive probe of the gravitational field in the Galactic centre. Near pericentre at 120 AU ≈ 1400 Schwarzschild radii, the star has an orbital speed of ≈7650 km s−1, such that the first-order effects of Special and General Relativity have now become detectable with current capabilities. Over the past 26 years, we have monitored the radial velocity and motion on the sky of S2, mainly with the SINFONI and NACO adaptive optics instruments on the ESO Very Large Telescope, and since 2016 and leading up to the pericentre approach in May 2018, with the four-telescope interferometric beam-combiner instrument GRAVITY. From data up to and including pericentre, we robustly detect the combined gravitational redshift and relativistic transverse Doppler effect for S2 of z = Δλ / λ ≈ 200 km s−1/c with different statistical analysis methods. When parameterising the post-Newtonian contribution from these effects by a factor f , with f = 0 and f = 1 corresponding to the Newtonian and general relativistic limits, respectively, we find from posterior fitting with different weighting schemes f = 0.90 ± 0.09|stat ± 0.15|sys. The S2 data are inconsistent with pure Newtonian dynamics.

As part of the Megamaser Cosmology Project, here we present a new geometric distance measurement to the megamaser galaxy NGC 5765b. Through a series of very long baseline interferometry observations, we have confirmed the water masers trace a thin, sub-parsec Keplerian disk around the nucleus, implying an enclosed mass of 4.55 ± 0.40 × 107 M ⊙. Meanwhile, from single-dish monitoring of the maser spectra over two years, we measured the secular drifts of maser features near the systemic velocity of the galaxy with rates between 0.5 and 1.2 km s−1 yr−1. Fitting a warped, thin-disk model to these measurements, we determine a Hubble Constant H 0 of 66.0 ± 6.0 km s−1 Mpc−1 with an angular-diameter distance to NGC 5765b of 126.3 ± 11.6 Mpc. Apart from the distance measurement, we also investigate some physical properties related to the maser disk in NGC 5765b. The high-velocity features are spatially distributed into several clumps, which may indicate the existence of a spiral density wave associated with the accretion disk. For the redshifted features, the envelope defined by the peak maser intensities increases with radius. The profile of the systemic masers in NGC 5765b is smooth and shows almost no structural changes over the two years of monitoring time, which differs from the more variable case of NGC 4258.

We present a 0.16% precise and 0.27% accurate determination of R0, the distance to the Galactic center. Our measurement uses the star S2 on its 16-year orbit around the massive black hole Sgr A* that we followed astrometrically and spectroscopically for 27 years. Since 2017, we added near-infrared interferometry with the VLTI beam combiner GRAVITY, yielding a direct measurement of the separation vector between S2 and Sgr A* with an accuracy as good as 20 μas in the best cases. S2 passed the pericenter of its highly eccentric orbit in May 2018, and we followed the passage with dense sampling throughout the year. Together with our spectroscopy, in the best cases with an error of 7 km s−1, this yields a geometric distance estimate of R0 = 8178 ± 13stat. ± 22sys. pc. This work updates our previous publication, in which we reported the first detection of the gravitational redshift in the S2 data. The redshift term is now detected with a significance level of 20σ with fredshift = 1.04 ± 0.05.

We report the detection of continuous positional and polarization changes of the compact source SgrA* in high states (“flares”) of its variable near-infrared emission with the near-infrared GRAVITY-Very Large Telescope Interferometer (VLTI) beam-combining instrument. In three prominent bright flares, the position centroids exhibit clockwise looped motion on the sky, on scales of typically 150 μas over a few tens of minutes, corresponding to about 30% the speed of light. At the same time, the flares exhibit continuous rotation of the polarization angle, with about the same 45(±15) min period as that of the centroid motions. Modelling with relativistic ray tracing shows that these findings are all consistent with a near face-on, circular orbit of a compact polarized “hot spot” of infrared synchrotron emission at approximately six to ten times the gravitational radius of a black hole of 4 million solar masses. This corresponds to the region just outside the innermost, stable, prograde circular orbit (ISCO) of a Schwarzschild–Kerr black hole, or near the retrograde ISCO of a highly spun-up Kerr hole. The polarization signature is consistent with orbital motion in a strong poloidal magnetic field.