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.

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.

We present the calibration and reduction of Event Horizon Telescope (EHT) 1.3 mm radio wavelength observations of the supermassive black hole candidate at the center of the radio galaxy M87 and the quasar 3C 279, taken during the 2017 April 5-11 observing campaign. These global very long baseline interferometric observations include for the first time the highly sensitive Atacama Large Millimeter/submillimeter Array (ALMA); reaching an angular resolution of 25 μas, with characteristic sensitivity limits of ∼1 mJy on baselines to ALMA and ∼10 mJy on other baselines. The observations present challenges for existing data processing tools, arising from the rapid atmospheric phase fluctuations, wide recording bandwidth, and highly heterogeneous array. In response, we developed three independent pipelines for phase calibration and fringe detection, each tailored to the specific needs of the EHT. The final data products include calibrated total intensity amplitude and phase information. They are validated through a series of quality assurance tests that show consistency across pipelines and set limits on baseline systematic errors of 2% in amplitude and 1°in phase. The M87 data reveal the presence of two nulls in correlated flux density at ∼3.4 and ∼8.3 Gλ and temporal evolution in closure quantities, indicating intrinsic variability of compact structure on atimescale of days, or several light-crossing times for afew billion solar-mass black hole. These measurements provide the first opportunity to image horizon-scale structure in M87.

We present the results of the Gould's Belt Distances Survey of young star-forming regions toward the Orion Molecular Cloud Complex. We detected 36 young stellar objects (YSOs) with the Very Large Baseline Array, 27 of which have been observed in at least three epochs over the course of two years. At least half of these YSOs belong to multiple systems. We obtained parallax and proper motions toward these stars to study the structure and kinematics of the Complex. We measured a distance of 388±5 pc toward the Orion Nebula Cluster, 428±10 pc toward the southern portion L1641, 388±10 pc toward NGC 2068, and roughly ∼420 pc toward NGC 2024. Finally, we observed a strong degree of plasma radio scattering toward λ Ori.

We present large-scale (∼ 2000 square arcminutes), deep (∼ 20 µJy), high-resolution (∼ 1 ) radio observations of the Ophiuchus star-forming complex obtained with the Karl G. Jansky Very Large Array at λ = 4 and 6 cm. In total, 189 sources were detected, 56 of them associated with known young stellar sources, and 4 with known extragalactic objects; the other 129 remain unclassified, but most of them are most probably background quasars. The vast majority of the young stars detected at radio wavelengths have spectral types K or M, although we also detect 4 objects of A/F/B types and 2 brown dwarf candidates. At least half of these young stars are non-thermal (gyrosynchrotron) sources, with active coronas characterized by high levels of variability, negative spectral indices, and (in some cases) significant circular polarization. As expected, there is a clear tendency for the fraction of non-thermal sources to increase from the younger (Class 0/I or flat spectrum) to the more evolved (Class III or weak line T Tauri) stars. The young stars detected both in X-rays and at radio wavelengths broadly follow a Güdel-Benz relation, but with a different normalization than the most radio-active types of stars. Finally, we detect a ∼ 70 mJy compact extragalactic source near the center of the Ophiuchus core, which should be used as gain calibrator for any future radio observations of this region.

We present the first results of the Gould's Belt Distances Survey (GOBELINS), a project aimed at measuring the proper motion and trigonometric parallax of a large sample of young stars in nearby regions using multi-epoch Very Long Baseline Array (VLBA) radio observations. Enough VLBA detections have now been obtained for 16 stellar systems in Ophiuchus to derive their parallax and proper motion. This leads to distance determinations for individual stars with an accuracy of 0.3 to a few percent. In addition, the orbits of six multiple systems were modelled by combining absolute positions with VLBA (and, in some cases, near-infrared) angular separations. Twelve stellar systems are located in the dark cloud Lynds1688; the individual distances for this sample are highly consistent with one another and yield a mean parallax for Lynds1688 of v =  7.28 0.06 mas, corresponding to a distance =  d 137.3 1.2 pc. This represents an accuracy greater than 1%. Three systems for which astrometric elements could be measured are located in the eastern streamer (Lynds 1689) and yield an estimate of v =  6.79 0.16 mas, corresponding to a distance =  d 147.3 3.4 pc. This suggests that the eastern streamer is located about 10 pc farther than the core, but this conclusion needs to be confirmed by observations of additional sources in the eastern streamer (currently being collected). From the measured proper motions, we estimate the one-dimensional velocity dispersion in Lynds 1688 to be 2.8±1.8 and 3.0±2.0km s −1 , in R.A. and decl., respectively; these are larger than, but still consistent within s 1 of, those found in other studies.

In 2017 April, the Event Horizon Telescope (EHT) observed the near-horizon region around the supermassive black hole at the core of the M87 galaxy. These 1.3 mm wavelength observations revealed a compact asymmetric ring-like source morphology. This structure originates from synchrotron emission produced by relativistic plasma located in the immediate vicinity of the black hole. Here we present the corresponding linear-polarimetric EHT images of the center of M87. We find that only a part of the ring is significantly polarized. The resolved fractional linear polarization has a maximum located in the southwest part of the ring, where it rises to the level of ∼15%. The polarization position angles are arranged in a nearly azimuthal pattern. We perform quantitative measurements of relevant polarimetric properties of the compact emission and find evidence for the temporal evolution of the polarized source structure over one week of EHT observations. The details of the polarimetric data reduction and calibration methodology are provided. We carry out the data analysis using multiple independent imaging and modeling techniques, each of which is validated against a suite of synthetic data sets. The gross polarimetric structure and its apparent evolution with time are insensitive to the method used to reconstruct the image. These polarimetric images carry information about the structure of the magnetic fields responsible for the synchrotron emission. Their physical interpretation is discussed in an accompanying publication.

Aims. We take advantage of the second data release of the Gaia space mission and the state-of-the-art astrometry delivered from very long baseline interferometry observations to revisit the structure and kinematics of the nearby Taurus star-forming region. Methods. We apply a hierarchical clustering algorithm for partitioning the stars in our sample into groups (i.e., clusters) that are associated with the various molecular clouds of the complex, and derive the distance and spatial velocity of individual stars and their corresponding molecular clouds. Results. We show that the molecular clouds are located at different distances and confirm the existence of important depth effects in this region reported in previous studies. For example, we find that the L 1495 molecular cloud is located at d = 129.9 +0.4 −0.3 pc, while the filamentary structure connected to it (in the plane of the sky) is at d = 160.0 +1.2 −1.2 pc. We report B 215 and L 1558 as the closest (d = 128.5 +1.6 −1.6 pc) and most remote (d = 198.1 +2.5 −2.5 pc) substructures of the complex, respectively. The median inter-cloud distance is 25 pc and the relative motion of the subgroups is on the order of a few km/s. We find no clear evidence for expansion (or contraction) of the Taurus complex, but signs of the potential effects of a global rotation. Finally, we compare the radial velocity of the stars with the velocity of the underlying 13 CO molecular gas and report a mean difference of 0.04 ± 0.12 km/s (with r.m.s. of 0.63 km/s) confirming that the stars and the gas are tightly coupled.