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 twostage 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.
['black hole physics', 'galaxies: individual: M87', 'galaxies: jets', 'techniques: high angular resolution', 'techniques: image processing', 'techniques: interferometric', 'Astrophysics - Astrophysics of Galaxies', 'Astrophysics - Instrumentation and Methods for Astrophysics', 'General Relativity and Quantum Cosmology']
First M87 Event Horizon Telescope Results. IV. Imaging the Central Supermassive Black Hole
{'No Header': '<!-- image -->', 'First M87 Event Horizon Telescope Results. IV. Imaging the Central Supermassive Black Hole': 'The Event Horizon Telescope Collaboration ( See the end matter for the full list of authors. )  \nReceived 2019 February 11; revised 2019 March 5; accepted 2019 March 6; published 2019 April 10', 'Abstract': "We present the /uniFB01 rst 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 /uni03BC 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 /uniFB01 rst 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.  \nKey words: black hole physics -galaxies: individual ( M87 ) -galaxies: jets -techniques: high angular resolution -techniques: image processing -techniques: interferometric", '1. Introduction': "Since the discovery of the /uniFB01 rst astrophysical jet apparently connected to its nucleus /uniF0A0 ( Curtis 1918 ) , the giant elliptical galaxy M87 in the Virgo cluster has been intensively studied with imaging observations. M87 ' s nuclear gas and stellar dynamics, as traced by optical and infrared ( IR ) spectroscopy, suggest the presence of a nuclear supermassive black hole ( SMBH ) of mass M BH /uniF0A0 ∼ /uniF0A0 ( 3.3 -6.2 ) /uniF0A0 × /uniF0A0 10 9 M e ( Macchetto et al. 1997;Gebhardt& Thomas 2009;Gebhardtetal.2011; Walsh et al. 2013 ) .Thishigh mass, combined with its proximity ( D 16.8 Mpc; = Blakeslee et al. 2009;Birdetal.2010; Cantiello et al. 2018;seealsoEHT Collaboration et al. 2019e, hereafter Paper VI ) , implies that the nuclear black hole candidate in M87 ( hereafter referred to as M87 ) has an event horizon subtending the second-largest known angular size after Sagittarius A * ( Sgr A * ) in the Galactic Center.  \nUnlike Sgr /uniF0A0 A * , M87 hosts a powerful kpc-long jet that is bright in the radio, optical, and X-ray bands ( e.g., Owen et al. 1989; Sparks et al. 1996; Perlman et al. 1999;Marshalletal.2002 ) . Weak emission east of the very long baseline interferometry ( VLBI ) core from the expected counter-jet has also been detected in high-frequency VLBI images ( Walker et al. 2018 ) . Material moves down the approaching jet with a maximum apparent speed of c 6 ~ ( Biretta et al. 1999 ) . On pc and sub-pc scales, VLBI observations show the jet to be edge-brightened and parabolic in shape ( Reid et al. 1989; Dodson et al. 2006; Kovalev et al. 2007; Asada & Nakamura 2012;Hadaetal.2013; Nakamura & Asada 2013; Asada et al. 2016 ) with a characteristic progressive acceleration downstream ( Asada et al. 2014; Mertens et al. 2016; Britzen et al. 2017;Hadaetal.2017;Walkeretal.2018;  \n<!-- image -->  \nOriginal content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution of this work must maintain attribution to the author ( s ) and the title of the work, journal citation and DOI.  \nKim et al. 2018a ) . High-frequency astrometric VLBI measurements reveal a frequency-dependent shift of the radio core ( from optical depth effects ) , which asymptotically converges to ∼ 40 microarcseconds ( /uni03BC as ) east of the 7 mm core /uniF0A0 ( Hada et al. 2011 ) ; this indicates that the jet is launched in the vicinity of the central black hole ( e.g., Nakamura et al. 2018 ) residing within the central ∼ 100 /uni03BC as. The high mass and relative proximity of M87 provides an opportunity to image this black hole and jet-launching region on event-horizon scales; however, accessing these scales with ground-based VLBI requires observations with microarcsecond resolution at a wavelength of /uniF088 1 mm.  \nTo this end, we have developed the Event Horizon Telescope ( EHT ) , a global ad hoc VLBI array operating at 1.3 mm wavelength ( EHT Collaboration et al. 2019b, hereafter Paper II ) . With its longest baselines spanning nearly the diameter of the Earth, the synthesized beam size of the EHT array is approximately 20 /uni03BC as. For M87, the EHT beam size corresponds to 3 -5 R s,where the Schwarzschild radius R GM c 2 s BH 2 = subtends 3.9 -7.3 /uni03BC as for the black hole mass range and distance given above. Thus, the EHT can potentially resolve general relativistic effects associated with the SMBH in M87, most notably the ' shadow ' cast by the black hole on the bright surrounding emission ( Bardeen 1973; Luminet 1979; Falcke et al. 2000 ) . This shadow is expected to be encircled by a bright ring at the radius of the lensed photon sphere, with a diameter between approximately 4.8 and 5.2 R s for a maximally spinning black hole ( viewed face-on ) and a nonspinning ( i.e., Schwarzschild ) black hole, respectively ( Bardeen 1973; Johannsen & Psaltis 2010 ) . For M87, the expected shadow diameter is 19 -38 /uni03BC as. Physical models and general relativistic magnetohydrodynamic ( GRMHD ) simulations show that Doppler-boosted emission from rapidly rotating material near the black hole can result in substantial image brightness asymmetry very near the ring ( EHT Collaboration et al. 2019d, hereafter Paper V ) .  \nEarly EHT observations in 2009 and 2012 detected compact emission with an FWHM size of approximately 40 /uni03BC as ( Doeleman et al. 2012; Akiyama et al. 2015 ) . However, because of their limited interferometric baseline coverage, these early experiments could not synthesize an image of M87, leaving considerable uncertainty about the nature of the detected emission.  \nIn 2017 April, the EHT conducted an observing campaign using eight stations in six geographic sites. 105 The EHT observed M87 on four days ( April 5, 6, 10, and 11 ) , interleaved with observations of other targets. Notably, these observations included 37 telescopes of the Atacama Large Millimeter / submillimeter Array ( ALMA ) coherently combined to act as a single 73 m diameter telescope ( Matthews et al. 2018 ) . The addition of baselines to ALMA signi /uniFB01 cantly increases the sensitivity of the entire EHT array. Additional details of the EHT instrument are given in Paper II; details of the 2017 observations, correlation, and calibration are given in EHT Collaboration et al. ( 2019c, hereafter Paper III ) .  \nWe generated images of M87 from the 2017 EHT data in two stages. In the /uniFB01 rst stage, our aim was to compare the results of four independent imaging teams. Each team used their collective judgment to produce an image, while remaining blind to the others ' work. In the second stage, our aim was to minimize reliance on expert judgment, and we instead utilized three imaging pipelines to perform large searches over the imaging parameter spaces. We then selected /uniFB01 ducial imaging parameters for M87 from these searches based on the performance of the parameters when reconstructing images from synthetic data. In both stages, the reconstructed images were consistently dominated by a single feature: a bright ring with a ∼ 40 /uni03BC as diameter and azimuthal brightness asymmetry, consistent with the shadow of a ∼ 6.5 /uniF0A0 × /uniF0A0 10 9 M e SMBH in the nucleus of M87.  \nThis Letter presents details of this imaging procedure and analysis; our primary images of M87 are shown in Figures 11 and 15. We begin, in Section 2, by providing a brief review of interferometric measurements and imaging techniques. In Section 3, we describe the EHT observations of M87 in 2017. In Section 4, we estimate image properties for M87 such as total compact /uniFB02 ux density and size of the emission region using a nonimaging analysis of the EHT data in combination with results from VLBI imaging of M87 at longer wavelengths. In Section 5,we present the /uniFB01 rst reconstructed EHT images, generated by four separate imaging teams working independently. In Section 6,we derive /uniFB01 ducial imaging parameters for three different imaging pipelines by performing large parameter surveys and testing on synthetic data. In Section 7, we show the corresponding /uniFB01 ducial images from each pipeline, and we assess their properties and uncertainties in the image and visibility domains. In Section 8,we perform validation tests of the /uniFB01 ducial images via comparisons of gain solutions with the simultaneously observed target 3C 279 and by exploring the dependence of reconstructed images on the participation and calibration of each site. In Section 9,wegivea quantitative analysis of the image and ring properties. In Section 10, we summarize our results.", '2. Background': 'Every antenna i in an interferometric array records the incoming complex electric /uniFB01 eld as a function of time, frequency, and polarization: Ei ( t , ν , P ) . Every pair of sites ( i , j ) then gives an  \ninterferometric visibility, de /uniFB01 ned as the complex cross-correlation between their recorded electric /uniFB01 elds,  \nV t P P E t P E t P , , , , , , , . 1 ij i j 1 2 1 2 * n n n =á ñ ( ) ( ) ( ) ( )  \nIn practice, radio telescopes record data in dual circular feeds, right circular polarization ( RCP ) and left circular polarization ( LCP ) , or in an orthogonal linear basis. The set of four possible cross-correlations among the two recorded polarizations at the two sites then provides information about the four Stokes parameters ( see, e.g., Roberts et al. 1994 ) . Because this Letter is focused on monochromatic and total intensity imaging, we will suppress the frequency and polarimetric parameters for the remainder of our discussion.  \nBy the van Cittert -Zernike theorem ( van Cittert 1934; Zernike 1938 ) , the visibility measured with an ideal interferometer is related to the brightness distribution on the sky I ( x , y ) via a simple Fourier transform. The interferometer samples a spatial frequency of the image given by the vector baseline b ij joining the sites, projected orthogonal to the line of sight and measured in wavelengths: u b ij ij , l = ^ , where λ /uniF0A0 = /uniF0A0 c / ν . For a pair of telescopes at /uniFB01 xed locations on the Earth, their projected baseline coordinates ( u , v ) trace an ellipse as the Earth rotates, sampling a range of spatial frequencies ( Thompson et al. 2017 ) . Generalizing to an arbitrary projected baseline u uv , , = ( ) the ideal visibility  is given by  \nu v e I x y dx dy , , . 2 iux vy 2  = p -+ /uni222C ( ) ( ) ( ) ( )  \nIn this expression, x and y are angular coordinates on the sky in radians. Equation ( 2 ) assumes that the source radiation is spatially incoherent and in the far/uniFB01 eld limit. It also assumes that the sky brightness distribution is only non-zero within both a radian and the primary beams of the individual telescopes. All of these criteria are satis /uniFB01 ed for EHT observations. 106 In the remainder of this section, we describe the speci /uniFB01 c interferometric data products used for imaging from EHT data, and we discuss how these data products are affected by practical measurement limitations ( Section 2.1 ) . We then summarize the primary methods used to reconstruct images from interferometric data ( Section 2.2 ) . For a comprehensive discussion of the EHT data reduction pipeline, see Paper III.', '2.1. VLBI Data Products and Their Uncertainties': "In practice, each measured visibility is contaminated by thermal noise, station-dependent errors, and baseline-dependent errors. However, some combinations of visibilities ( so-called ' closure ' quantities ) can be used to construct data products that are independent of station-dependent errors. We now discuss these various data products and de /uniFB01 ne how they will be used to assess consistency between a trial image I and the data.", '2.1.1. Interferometric Visibilities': "The fundamental interferometric data product is the complex visibility, V ij ( see Equation ( 1 )) . Each visibility is complexvalued and obeys a conjugation relationship under baseline reversal, V V ij ji * = . Thus, an N el-element array samples at most N el ( N el -1 ) / 2 different visibilities at each time.  \nVisibility measurements have a simple thermal error budget. Neglecting errors in calibration ( Section 2.1.2 ) , the measured visibility V ij on a baseline u ij is given by the true visibility u ij ij   = ( ) plus baseline-dependent thermal noise ij  : Vij ij ij   = + . For perfect coherent averaging, the thermal noise on measurements obeys a complex normal distribution, with equal variance on the real and imaginary parts. The thermal noise standard deviation, σ ij , is related to the system equivalent /uniFB02 ux densities ( SEFDs ) of the two sites, the bandwidth /uni0394 ν , and the coherent integration time /uni0394 t :  \nt 1 SEFD SEFD 2 , 3 ij i j Q s h n = ' DD ( )  \nwhere η Q is a factor that accounts for sensitivity loss from quantization of the incident electric /uniFB01 eld. The EHT utilizes standard 2-bit recording, for which η Q /uniF0A0 = /uniF0A0 0.88 ( Thompson et al. 2017 ) .  \nWe quantify agreement between a trial image I and a set of measured visibilities using the mean squared standardized residual. Speci /uniFB01 cally, for a measured set of N V visibilities, we de /uniFB01 ne  \nI N V V 1 2 , 4 vis 2 V 2 2 å c s = -( ) /uni2223 ˆ /uni2223 ( )  \nwhere V ˆ denotes a model visibility corresponding to the trial image, and the sum ranges over all measured visibilities ( which may include different baselines, times, and frequencies ) . The quantity vis 2 c does not formally correspond to a reduced χ 2 , and its interpretation requires caution. For example, we do not include a correction for the effective number of image degrees of freedom, and the results can be sensitive to data processing procedures such as self-calibration. Nevertheless, we expect 1 vis 2 c » for /uniFB01 nal imaging results.  \nWe also use a mean squared standardized residual to quantify agreement between a trial image and a set of measured visibility amplitudes:  \nI N A A 1 , 5 amp 2 V 2 2 å c s = -( ) /uni2223 ˆ /uni2223 ( )  \nwhere A V = ˆ /uni2223 ˆ /uni2223 denotes a model visibility amplitude, and  \nA V V 6 2 2 s s = -Q -(/uni2223 /uni2223 ) (/uni2223 /uni2223 ) ( )  \nis a measured visibility amplitude after correcting for an upward bias from thermal noise. Here, the Heaviside Θ -function ensures that the debiased amplitudes are always real and positive.", '2.1.2. Closure Quantities': "In addition to thermal noise, measured visibilities have systematic errors on both their amplitudes and phases. These errors can often be factorized into multiplicative station-based ' gains, ' which are complex-valued and may vary in both time and frequency. In general, these gains act as linear transformations of the incident polarized electric /uniFB01 elds at each site; here, we focus on a simpler model using scalar electric /uniFB01 elds. These complex systematic errors modify the relationship between a measurement V ij and the corresponding ideal visibility ij  as  \nV g g , 7 ij i j ij ij *   » + ( )  \nwhere ò ij is the thermal noise. Equation ( 7 ) is an approximation, as it neglects cross-talk between orthogonal receivers, the nonfactorizability of the station-dependent bandpass after frequency averaging, the non-factorizability of station-dependent temporal /uniFB02 uctuations after coherent averaging, and other imperfections. We discuss these imperfections and their implications for the EHT error budget in Section 4.2 ( for a comprehensive discussion, see Paper III ) . 107  \nThe site dependence of the gains allows for the construction of data products that are insensitive to these systematic errors. Such quantities are called ' closure ' quantities. For instance, the sum of three visibility phases ψ /uniF0A0 = /uniF0A0 arg ( V ) around a directed triangle i -j -k of sites is insensitive to the phase of the complex station-based gains gi . Namely, neglecting thermal noise,  \n, 8 ijk ij jk ki ij jk ki C, y y y y = + + = Y + Y + Y ( )  \nwhere arg  Y = ( ) is the phase of the ideal visibility. This phase ijk C, y is called the closure phase of the i -j -k triangle ( Jennison 1958 ) . Because phase stability is challenging for millimeter VLBI, closure phase is a useful measurement of intrinsic source phase. For an N el-element interferometer, there are N el ( N el -1 )( N el -2 ) / 6 closure phase measurements, but only ( N el -1 )( N el -2 ) / 2 of these are non-redundant.  \nWhen the signal-to-noise ratio ( S / N ) of each component visibility is at least ∼ 3, closure phase errors are well approximated as Gaussian and depend only on the three component S / N values S ( e.g., Rogers et al. 1995 ) :  \nS S S . 9 ijk ij jk ki , 2 2 2 C s » + + y ---( )  \nIn this expression, C s y is the closure phase standard deviation in radians.  \nWe compute CP 2 c statistics to quantify image-data agreement analogously to Equation ( 4 ) :  \nI N e e 1 , 10 i i CP 2 2 2 C C C C å c s = -y y y y ( ) /uni2223 /uni2223 ( ) ˆ  \nwhere C y ˆ denotes a model closure phase from the trial image, N C y is the number of closure phases in the sum, and the sum ranges over all measurements ( which may include different triangles, times, and frequencies ) .  \nSimilar to closure phase, certain combinations of four visibility amplitudes on a quadrangle ijk ℓ are insensitive to gain amplitudes g i /uni2223 /uni2223 . For example, neglecting thermal noise,  \nA AA A A , 11 ijkℓ ij kℓ ik jℓ ij kℓ ik jℓ C,   = = ( )  \nwhere   = /uni2223 /uni2223 . Such combinations A ijkℓ C, are called closure amplitudes ( Twiss et al. 1960 ) . While N el sites give at most N el ( N el -1 )( N el -2 )( N el -3 ) / 8 closure amplitudes, only N el ( N el -3 ) / 2 of these are non-redundant. Also, while the noise statistics of closure amplitudes are approximately Gaussian if each of the four amplitudes has high S / N, the distribution is highly non-Gaussian when any of the measurements is low S / N, especially those in the denominator ( L. Blackburn et al. 2019, in preparation; A. E. Broderick et al. 2019, in preparation ) .To  \nmitigate this behavior and to symmetrize the properties of numerator and denominator terms, we use the logarithm of the closure amplitudes, A ln C , for analysis. We compute log CA 2 c statistics as  \nI N A A 1 ln ln , 12 log CA 2 lnA C C 2 lnA 2 C C å c s = -( ) /uni2223 ˆ /uni2223 ( )  \nwhere the uncertainty on the log closure amplitude is, in the high-S / N limit,  \nS S S S . 13 A ijkℓ ij kℓ ik jℓ ln , 2 2 2 2 C s » + + + ----( )  \nNote that closure quantities are generally not independent, and their uncertainties are only approximately Gaussian.", '2.2. Interferometric Imaging Methods': 'Because an interferometer incompletely samples the visibility domain, the inverse problem of determining an image from a measured set of visibilities is ill-posed. Consequently, reconstructed images are not unique -they always require information, assumptions, or constraints beyond the interferometer measurements. One strong imaging constraint is image positivity. A second strong constraint is a restricted /uniFB01 eld of view ( FOV ) for the reconstructed image ( i.e., the gridded FOV ) . 108 The choice of FOV must be made with care, as incorrect restrictions can result in false image structure. Imaging algorithms can additionally impose or favor other physically motivated properties related to the image ( e.g., image smoothness ) or to the instrument ( e.g., a maximal resolution of reconstructed features ) .  \nImaging algorithms can be broadly categorized into two methodologies: inverse modeling and forward modeling. The former includes deconvolution methods such as CLEAN ( Högbom 1974; Clark 1980 ) , while the latter includes regularized maximum likelihood ( RML ) methods such as the classical maximum entropy method ( MEM; e.g., Narayan & Nityananda 1986 ) . We now discuss each of these in more depth, with an emphasis on aspects and algorithms that are most relevant for EHT imaging.', '2.2.1. Imaging via Inverse Modeling': "Conventional inverse modeling approaches to imaging typically begin with an inverse Fourier transform of the sampled visibilities ( giving the so-called ' dirty image ' ) . They then deconvolve the effects associated with the limited baseline coverage ( giving the so-called ' dirty beam ' ) . For both VLBI and connected-element radio interferometry, the standard reconstruction algorithm in this category is CLEAN ( e.g., Högbom 1974; Schwarz 1978; Clark 1980; Schwab 1984; Cornwell et al. 1999; Cornwell 2008 ) .  \nThe classical CLEAN algorithm ( Högbom 1974; Clark 1980 ) models the image as a collection of point sources and determines the locations and /uniFB02 ux densities of these point sources iteratively. After reaching a speci /uniFB01 ed stopping criterion, CLEAN typically convolves the many-point-source image model with a ' clean beam. ' This beam is obtained by  \nmatching a Gaussian to the central component of the dirty beam, and it approximates the point-spread function of the interferometric data. However, it is the unconvolved pointsource model that is compared with the data to assess goodness of /uniFB01 t; the beam-convolved image displayed will not /uniFB01 t the measured visibility amplitudes.  \nLack of absolute phase information and a priori calibration uncertainties in EHT measurements require multiple consecutive iterations of CLEAN followed by self-calibration, a routine that solves for station gains to maximize consistency with visibilities of a speci /uniFB01 ed trial image ( Wilkinson et al. 1977; Readhead et al. 1980; Cornwell & Wilkinson 1981; Pearson & Readhead 1984; Cornwell & Fomalont 1999 ) . For the CLEAN imaging in this Letter, we used the DIFMAP software ( Shepherd 1997, 2011 ) .", '2.2.2. Imaging via Forward Modeling': "Forward modeling approaches to imaging usually represent the image as an array of pixels and only require a Fourier transform of this array to evaluate consistency between the image and data. These methods can easily integrate nonlinear image-data consistency measures and constraints on the reconstructed image properties, such as sparsity or smoothness. While these types of imaging algorithms have been developed for decades ( e.g., Frieden 1972; Gull & Daniell 1978; Cornwell & Evans 1985; Narayan & Nityananda 1986; Briggs 1995; Wiaux et al. 2009a, 2009b ) and are commonly utilized in optical interferometry ( e.g., Buscher 1994; Baron et al. 2010; Thiébaut 2013; Thiébaut & Young 2017 ) , they are used less frequently than CLEAN for radio interferometry. However, forward modeling methods have been intensively developed for the EHT ( e.g., Honma et al. 2014; Bouman et al. 2016, 2018; Chael et al. 2016, 2018; Ikeda et al. 2016; Akiyama et al. 2017a, 2017b; Johnson et al. 2017; Kuramochi et al. 2018 ) .  \nThese methods are often derived using probabilistic arguments, although they do not typically produce results that have formal probabilistic interpretations. The general approach in this type of imaging is to /uniFB01 nd the image I that minimizes a speci /uniFB01 ed objective function,  \nJI I SI . 14 D D R R data terms 2 regularizers å å ac b = -( ) ( ) ( ) ( )  \nIn this expression, each D 2 c is a goodness-of/uniFB01 t function corresponding to the data term D ( Section 2.1 ) , and each S R is a regularization term corresponding to the regularizer R ( Appendix A ) . The ' hyperparameters ' α D and β R control the relative weighting of the data and regularization terms. Placing too much weight on the regularization terms will result in /uniFB01 nal reconstructed images that are not acceptably compatible with the data ( i.e., the image will not achieve 1 D 2 c ~ for all data terms ) . Because this function is analogous to the log-likelihood of a posterior probability function, we refer to this category of imaging methods as RML methods.  \nRegularizers explored in VLBI include image entropy /uniF0A0 ( e.g., Narayan & Nityananda 1986 ) , smoothness /uniF0A0 ( e.g., Bouman et al. 2016; Chael et al. 2016; Kuramochi et al. 2018 ) , /uniF0A0 and sparsity in the image or its gradient domain /uniF0A0 ( e.g., Wiaux et al. 2009a, 2009b; Honma et al. 2014; Akiyama et al. 2017b ) . Regularizing functions can be combined; e.g., simultaneously favoring both sparsity and smoothness through regularization can  \nmitigate limitations of using only one or the other /uniF0A0 ( e.g., Akiyama et al. 2017a ) .  \nNote that while regularization is an important choice in RML methods, other decisions such as data preparation, the optimization methodology, and the FOV of the reconstructed image can be equally important in in /uniFB02 uencing the /uniFB01 nal image. Unlike CLEAN, no /uniFB01 nal restoring beam is required when using RML, and these methods often achieve a modest degree of ' super resolution ' ( i.e., /uniFB01 ner than the nominal diffraction limit, u max q l ~ /uni2223 /uni2223 ) .  \nFor the EHT, RML methods mitigate some dif /uniFB01 culties of CLEAN reconstructions through increased /uniFB02 exibility in the data products used for imaging. For instance, RML methods can directly /uniFB01 t to robust data products such as closure quantities ( e.g., Buscher 1994; Baron et al. 2010; Bouman et al. 2016; Chael et al. 2016, 2018; Akiyama et al. 2017a, 2017b ) . However, RML methods often introduce other dif /uniFB01 culties, such as requiring the solution to a complex nonlinear optimization problem ( minimizing Equation ( 14 )) . As with CLEAN, RML methods sometimes employ an iterative imaging loop, alternately imaging and then self-calibrating the data. For the RML reconstructions in this Letter, we used two open-source software libraries that have been developed speci /uniFB01 cally for the EHT: eht-imaging 109 ( Chael et al. 2016, 2018, 2019 ) and SMILI 110 ( Akiyama et al. 2017a, 2017b, 2019 ) . Appendix A gives de /uniFB01 nitions of the regularizers used throughout this Letter, as implemented in these libraries.", '3.1. EHT Observations and Data': "The EHT observed M87 with seven stations at /uniFB01 ve geographic sites on 2017 April 5, 6, 10, and 11. The participating facilities were the phased Atacama Large Millimeter / submillimeter Array ( ALMA ) and Atacama Path /uniFB01 nder Experiment telescope ( APEX ) in the Atacama Desert in Chile, the James Clerk Maxwell Telescope ( JCMT ) and the phased Submillimeter Array ( SMA ) on Maunakea in Hawai ' i, the Arizona Radio Observatory Sub-Millimeter Telescope ( SMT ) on Mt. Graham in Arizona, the IRAM 30 m ( PV ) telescope on Pico Veleta in Spain, and the Large Millimeter Telescope Alfonso Serrano ( LMT ) on Sierra Negra in Mexico. These observations of M87 were interleaved with other targets ( e.g., the quasar 3C 279 ) , some of which were visible to an eighth EHT station, the South Pole Telescope ( SPT ) .  \nData were recorded in two polarizations and two frequency bands. All sites except ALMA and the JCMT recorded dual circular polarization ( RCP and LCP ) . ALMA recorded dual linear polarization that was later converted to a circular basis via PolConvert ( Martí-Vidal et al. 2016; Matthews et al. 2018; Goddi et al. 2019 ) , and the JCMT recorded a single circular polarization ( the recorded polarization varied from day to day ) . 111 All sites recorded two 2 GHz bands centered on 227.1 and 229.1 GHz ( henceforth, low and high band, respectively ) . Paper II provides details on the setup, equipment, and station upgrades leading up to the 2017 observations.  \nPaper III outlines the correlation, calibration, and validation of these data. In particular, the data reduction utilized the sensitive baselines to ALMA to estimate and correct for stable instrumental phase offsets, RCP -LCP delays, and stochastic phase variations within scans. After these corrections, the data have suf /uniFB01 cient phase stability to coherently average over scans. The data were also amplitude calibrated using station-speci /uniFB01 c measurements; stations with an intra-site partner ( i.e., ALMA, APEX, SMA, and JCMT ) were then ' network calibrated ' to further improve the amplitude calibration accuracy and stability via constraints among redundant baselines. The /uniFB01 nal networkcalibrated data sets were frequency averaged per band and coherently averaged in 10 s intervals before being used for our imaging analysis. All data presented and analyzed in this work are Stokes /uniF0A0 I ( or pseudo I ) visibilities processed via the EHT -Haystack Observatory Processing System ( HOPS ) pipeline ( Paper III; Blackburn et al. 2019 ) . Information about accessing SR1 data and the software used for analysis can be found on the Event Horizon Telescope website ' s data portal. 112", '3.2. Data Properties': 'Figure 1 shows the baseline ( u , v ) coverage for EHT observations of M87. The shortest baselines in the EHT are intra-site ( i.e., the SMA and JCMT are separated by 0.16 km; ALMA and APEX are separated by 2.6 km ) . These intra-site baselines are sensitive to arcsecond-scale structure ( see Section 4.3 ) . In contrast, our longest baselines ( joining the SMA or JCMT to PV ) are sensitive to microarcsecond-scale structure. Baseline coverage on individual days ( bottom panels of Figure 1 ) is comparable for April /uniF0A0 5, 6, and 11 ( 18, 25, and 22 scans, respectively ) . However, April /uniF0A0 10 had signi /uniFB01 cantly less coverage, with only 7 scans.  \nFigure 2 ( left panel ) shows the S / N as a function of baseline length for M87 on April /uniF0A0 11, after coherently averaging scans. The split in S / N distributions at various baseline lengths is due to the sharp difference in sensitivity for the co-located Atacama sites ALMA and APEX. The right panel of Figure 2 shows the visibility amplitude ( correlated /uniFB02 ux density ) for M87 on April 11 after amplitude and network calibration.  \nThere is a prominent secondary peak in the networkcalibrated visibility amplitudes between two deep minima ( \' nulls \' ) ,the /uniFB01 rst at ∼ 3.4 G λ and the second at ∼ 8.3 G λ .The amplitudes along the secondary peak are weakly dependent on baseline position angle, suggesting some degree of source symmetry, and their overall trends are consistent for all days ( see Paper III, Figure 13 ) . However, evidence for source anisotropy can be seen at the location of the /uniFB01 rst null, where the east -west oriented Hawai \' i -LMT baseline gives signi /uniFB01 -cantly lower amplitudes than the north -south oriented ALMA -LMT baseline at the same projected baseline length ( see also Paper VI ) . This anisotropy is further supported by multiple measurements of non-zero closure phase ( Figure 3 ) . The majority of notable low-amplitude outliers across days are due to reduced performance of the JCMT or the LMT on a select number of scans. Despite the amplitudes of these data being low, the derived closure quantities remain stable ( Paper III ) .  \nSimilarly, most closure quantities for M87 are broadly consistent across all days, although day-to-day variations are signi /uniFB01 cant for some sensitive closure combinations involving  \nFigure 1. Top panels: aggregate baseline coverage for EHT observations of M87, combining observations on all four days. The left panel shows short-baseline coverage, comprised of ALMA interferometer baselines and intra-site EHT baselines ( SMA -JCMT and ALMA -APEX ) . These short baselines probe angular scales larger than 0.1". The right panel shows long-baseline coverage, comprised of all inter-site EHT baselines. These long baselines span angular scales from 25 to 170 /uni03BC as. Each point denotes a single scan, which range in duration from 4 to 7 minutes. Bottom panels: the full baseline coverage on M87 for each observation. In all panels, the dashed circles show baseline lengths corresponding to the indicated fringe spacings ( 0.2" for the upper-left panel; 25 and 50 /uni03BC as for the remaining panels ) .  \n<!-- image -->  \nFigure 2. Left panel: S / N as a function of projected baseline length for EHT observations of M87 on April /uniF0A0 11. Each point denotes a visibility amplitude coherently averaged over a full scan ( 4 -7 minutes ) . Points are colored by baseline. Right panel: visibility amplitudes ( correlated /uniFB02 ux density ) as a function of projected baseline length after a priori and network calibration. The amplitudes are corrected for upward bias from thermal noise ( Equation ( 6 )) . Error bars denote ± 1 σ uncertainty from thermal noise and do not include expected uncertainties in the a /uniF0A0 priori calibration ( see Paper III and Section 4.1 ) .  \n<!-- image -->  \nlong baselines to PV or to the Hawai \' i stations. Figure 3 shows examples of closure phases for various triangles and levels of variability: a \' trivial \' triangle including co-located sites ( ALMA -APEX -SMT, left panel ) that is expected to be consistent with zero, a non-trivial and mostly non-variable  \ntriangle ( ALMA -LMT -SMT, center panel ) with largely persistent structure across all days, and a non-trivial triangle ( LMT -SMA -SMT, right panel ) showing intrinsic source structure evolution in M87 between the two sets of observations on April 5, 6 and 10, 11.  \nFigure 3. Selected closure phases from coherently averaged visibilities on three triangles as a function of Greenwich Mean Sidereal Time ( GMST ) using data from all four days. Error bars denote ± 1 σ uncertainties from thermal noise. The trivial ALMA -APEX -SMT triangle ( left panel ) has closure phases near zero on all days, as expected because this triangle includes an intra-site baseline. Deviations from zero arise from a combination of thermal and systematic errors ( Paper III ) . The ALMA -LMT -SMT triangle ( middle panel ) shows persistent structure across all days, while the large LMT -SMA -SMT triangle ( right panel ) shows source evolution between the /uniFB01 rst two days and last two days.  \n<!-- image -->', '4. Pre-imaging Considerations': 'In this section, we analyze the measured M87 visibilities directly to assess what conclusions are supported by the data irrespective of choices made in the imaging procedures. This non-imaging analysis provides estimates of the quality of the a /uniF0A0 priori calibration, the level of non-closing systematic errors, and the degree of temporal variability in the source. It also guides the generation of realistic synthetic data sets ( Section 6.1 ) and motivates the choice of imaging parameters ( e.g., the FOV of reconstructed images ) used to de /uniFB01 ne parameter ranges in the imaging parameter surveys ( Section 6.2 ) .', '4.1. Expected Amplitude Calibration Limitations': 'The amplitude calibration error budget is determined from uncertainties on individual measurements of station performance, as described in Paper III. The error budget is only representative of global station performance and is not speci /uniFB01 ed for individual measurements. While this procedure is adequate for stations with stable performance and weather during the observing run, the error budget may be underestimated for stations with variable performance. The SMT, PV, SMA, JCMT, APEX, and ALMA stations are well characterized either through years of studies or via network calibration. More recent additions to the EHT ( the LMT and the SPT ) may have more variable behavior, as their observing systems are not yet well exercised and because they do not have suf /uniFB01 ciently close sites to permit network calibration.  \nSpeci /uniFB01 cally for M87, the LMT is the most undercharacterized station. The LMT began observing M87 in the evening, when the dish is still affected by thermal gradients in the back-up structure or panel distortions from daytime heating, both of which are signi /uniFB01 cant for open-air telescopes. These effects are common for sensitive millimeter-wave dishes and cause surface instability. In addition, evening conditions are inadequate for accurate pointing and focusing of the telescope ( particularly on weaker sources ) , leading to performance that can vary substantially across scans and from day to day. A defocused dish can measure persistently low amplitudes on baselines to that station between focus attempts ( typically every one to two hours ) . Changes in telescope pointing can cause amplitudes to /uniFB02 uctuate signi /uniFB01 cantly from scan to scan ( from the telescope moving to and from the source ) and between pointing attempts ( typically every half-hour ) . Issues in telescope focus  \ncan also lead to uncertainties in the a priori calibration for other sources observed during the same time period, such as 3C 279. However, pointing errors for 3C 279 are expected to be less severe, as it is bright enough for the LMT to point directly onsource prior to VLBI scans. Thus, the corrections needed for the LMT are expected to better match the a priori amplitude error budget during observations of 3C 279 ( mostly correcting for focus errors ) than during observations of M87 ( correcting for both focus and pointing errors ) .  \nIn Section 8.1 and Appendix F, we compare estimated residual gains for the SMT and the LMT from imaging M87 and 3C 279. In Paper VI, we compare these results with the estimated residual gains when /uniFB01 tting parametric models to the interferometric data.', '4.2. Estimates of Non-closing Systematic Errors': 'Apart from the thermal errors, which are well understood and can be evaluated from /uniFB01 rst principles, there are several sources of additional systematic error present in the EHT data. These include polarimetric leakage ( proportional to the magnitude of station-dependent leakage terms and the source fractional polarization ) , residual polarimetric-gain calibration offsets, loss of amplitude with long coherent averaging, the effects of biased S / N estimators, and the limitations of a /uniF0A0 linear error propagation model for closure quantities.  \nThese systematic errors can affect both visibilities and closure quantities. While these effects are small, in a /uniF0A0 very high S / N regime systematic uncertainties may dominate over thermal uncertainties. In Paper III, we quanti /uniFB01 ed the magnitude of systematic uncertainties in the EHT 2017 data both with a series of statistical tests on the distributions of trivial closure quantities and with cross-validation tests of data products across different frequency bands and polarizations. For M87, the characteristic magnitude of systematic errors was found to be only a fraction of the thermal errors, even for heavily averaged data. In the case of 3C 279, which has signi /uniFB01 cant intrinsic polarization, systematic errors can dominate over the thermal errors. In both cases the characteristic magnitudes are small: less than 2 ° for closure phases, and less than 4% for closure amplitudes.', '4.3. Constraints on the Total Compact Flux Density': "EHT baselines sample angular scales in M87 that span nearly /uniFB01 ve orders of magnitude. The largest gap in coverage  \noccurs between the arcsecond scales sampled by intra-site baselines ( ALMA -APEX and SMA -JCMT ) and the microarcsecond scales sampled by inter-site baselines. The intra-site baselines are insensitive to detailed structure on microarcsecond scales; their measurements are consistent with an unresolved point source having ∼ 1.2 Jy of /uniFB02 ux density ( i.e., the ' core ' of M87, corresponding to the bright, upstream end of the jet ) in addition to faint compact features in the largescale jet of M87, such as ' Knot A ' ( Owen et al. 1989 ) . Long inter-site baselines resolve out much of this emission, and are sensitive to only the compact central source. The shortest intersite EHT baseline is SMT -LMT, which has a fringe spacing of 139 -166 /uni03BC as for M87. Thus, there are two relevant estimates of the total /uniFB02 ux density: the total /uniFB02 ux density seen in the core on arcsecond scales, F tot /uniF0A0 ≈ /uniF0A0 1.2 Jy, and the total compact /uniFB02 ux density, F cpct /uniF0A0 /uni0084 /uniF0A0 F tot, which is the relevant quantity for the integrated /uniFB02 ux density of a reconstructed image with an FOV limited to only a few hundred microarcseconds.  \nIn Appendix B,wederiveaseriesofconstraintson F cpct using the 2017 EHT data and quasi-simultaneous observations at longer wavelengths. These constraints depend on the effective size of the compact emission region ( expressed as an equivalent Gaussian FWHM ) . We also compare these constraints with the estimates of source size and compact /uniFB02 ux density reported in Doeleman et al. ( 2012 ) and Akiyama et al. ( 2015 ) from observations of M87 in 2009 and 2012, respectively. These previous results used a threesite array, which included the Combined Array for Research in Millimeter-wave Astronomy ( CARMA; in California ) in addition to the SMT, JCMT, and SMA. The addition of CARMA provided a short inter-site baseline ( SMT -CARMA, with u 0.6 G l ~ /uni2223 /uni2223 ) , while the two remaining baselines, SMT -JCMT / SMA and CARMA -JCMT / SMA, sampled u 2.2 3.5 G l < < /uni2223 /uni2223 .  \nTaken together, the EHT constraints give F 0.64 Jy cpct 0.08 0.39 = -+ ; combined with the multi-wavelength constraints, they give F 0.66 Jy cpct 0.10 0.16 = -+ . For both these estimates, the central value gives the maximum a posteriori estimate and the uncertainties give the 95% con /uniFB01 dence interval after marginalizing over the source size. Because all of the EHT constraints use the SMT -LMT baseline, they correspond to the source size projected along that direction. However, because the visibility amplitudes are highly symmetric, we expect these estimates to apply at all position angles ( see Section 3 ) .Weinferavalueof F cpct that is signi /uniFB01 cantly lower than what was found in observations in 2009 and 2012. 113  \nAs there is a range of possible values for the total compact /uniFB02 ux density, we do not enforce a common value for all reconstructed images of M87. Instead, before evaluating image-data consistency through χ 2 metrics, we add a large Gaussian component ( FWHM > 500 /uni03BC as ) to each reconstruction, which only affects visibilities on the intra-site baselines. The integrated /uniFB02 ux density of this Gaussian component is set so that the total /uniFB02 ux density of the modi /uniFB01 ed image matches the intra-site visibility amplitudes.  \nWe emphasize that the constraints on the size and total /uniFB02 ux density of the compact emission were derived by making rigid assumptions about the station gains, by combining constraints among all four observing days ( over which M87 may have  \nTable 1 Metrics of EHT Angular Resolution for the 2017 Observations of M87  \nNote. Because the EHT coverage is fairly symmetric, to avoid asymmetries introduced by restoring beams and to homogenize the images among epochs, we adopt a circular Gaussian restoring beam with 20 /uni03BC as FWHM for all CLEAN reconstructions.  \nvaried in F cpct and size ) , and by imposing short-baseline approximations for the visibility function. Thus, while we will use these constraints to guide our imaging and analysis, we will not strictly enforce them.", '4.4. Image Resolution and Degrees of Freedom': "The diffraction-limited resolution of an interferometer depends upon its /uniFB01 nest fringe spacing: u f min max q = q /uni2223 /uni2223 , where f θ is a scaling coef /uniFB01 cient. The resolution of reconstructed images also depends upon the baseline coverage and sensitivity and is commonly quanti /uniFB01 ed using procedures to /uniFB01 t the central peak of the point-spread function ( or ' dirty beam ' ) . Table 1 gives several common representations of the EHT resolution for observations of M87.  \nThe maximum number of independent degrees of freedom in a reconstructed image is then N dof fov min 2  q q ( ) , where θ fov is the FOV of the image. Of course, the image may only be nonzero in a small fraction of the full FOV, reducing N dof correspondingly. We emphasize that the number of degrees of freedom in a reconstructed image is not determined by the number of pixels in the image. Likewise, for a constrained FOV, it is not possible to arbitrarily /uniFB01 t ( or over /uniFB01 t ) interferometric data by making the pixel resolution increasingly /uniFB01 ne because the baseline correlation length for measurements in the visibility domain is determined by the FOV rather than by the pixel resolution. Consequently, a reconstructed image and its effective number of degrees of freedom is sensitive to the FOV of the image, but not to the chosen pixel size.", '4.5. Image Conventions': 'Throughout this Letter, we present images using their equivalent brightness temperatures de /uniFB01 ned by the Rayleigh -Jeans law: T I c k b 2 2 2 = n n , where I ν is the speci /uniFB01 c intensity, c is the speed of light, ν is the observing frequency, and k is the Boltzmann constant ( e.g., Rybicki & Lightman 1979 ) . We use brightness temperature rather than the standard radio convention of /uniFB02 ux density per beam ( e.g., Jy / beam ) because our images are spatially resolved and because RML methods do not have a natural associated beam. However, we emphasize that brightness temperature does not necessarily correspond to any  \nphysical temperature of the radio-emitting plasma. The radio spectrum of M87 is not a blackbody, and its 230 GHz emission is from synchrotron radiation ( Paper V ) . Finally, for visualization of our images, we use perceptually uniform colormaps from the ehtplot 114 library.  \nThroughout this Letter, we restore CLEAN images using a circular Gaussian beam with 20 /uni03BC as FWHM, which is comparable to the geometric mean of the principal axes of the CLEAN beam ( see Table 1 ) . Unless otherwise stated, any image that has been restored will include the restoring beam in the lower right. Also, for consistency with RML methods but in contrast with standard practice, our presented CLEAN images do not include the residual image, corresponding to the inverse Fourier transform of gridded residual visibilities. In Section 7.2, we describe characteristics of the residual images.', '5. First M87 Images from Blind Imaging': 'VLBI images are sensitive to choices made in the imaging and self-calibration process. Choices required in using any imaging method include deciding which data are used ( e.g., low and / or high band, /uniFB02 agging ) , specifying the self-calibration procedures, and /uniFB01 xing the reconstructed image FOV. In addition, imaging methods also require choices that are particular to their assumptions and methodology. For CLEAN, these choices include choosing a set of CLEAN windows and a data-weighting scheme. For RML methods, choices include the selection of which data and regularizer terms and weights to use in the objective function ( Equation ( 14 )) . With this abundance of user input, it can be dif /uniFB01 cult to assess what image properties are reliable from a given imaging method.  \nThe dangers of false con /uniFB01 dence and collective con /uniFB01 rmation bias are magni /uniFB01 ed for the EHT because the array has fewer sites than typical VLBI arrays, there are no previous VLBI images of any source at 1.3 mm wavelength, and there are no comparable black hole images on event-horizon scales at any wavelength. To minimize the risk of collective bias in /uniFB02 uencing our /uniFB01 nal images, in our /uniFB01 rst stage of analysis we reconstructed images of M87 in four independent imaging teams.', '5.1. Imaging Procedure and Team Structure': "We subdivided our /uniFB01 rst M87 imaging efforts into four separate imaging teams. The teams were blind to each others ' work, prohibited from discussing their imaging results and even from discussing aspects of the data that might in /uniFB02 uence imaging ( e.g., which stations or data might be of poor quality ) . No restrictions were imposed on the data pre-processing or imaging procedures used by each team. Teams 1 and 2 focused on RML methods, while Teams 3 and 4 primarily used CLEAN. In addition to independently imaging M87, teams also independently imaged other sources observed by the EHT in 2017 to test the blind imaging procedure.  \nBlind imaging procedures have long been used to reduce the risk of group bias. Prior to the 2017 observations, we organized a series of ' imaging challenges ' that used synthetic data to assess how conventional and newly developed imaging algorithms would perform for the EHT ( Bouman 2017 ) . 115 Reconstructing images independently in  \nthese challenges helped us identify which image features were likely intrinsic, and which were likely to be spurious. To compare EHT 2017 results among teams while keeping submissions blind, we built a website that allowed users to independently upload images and automatically compare them to the ground truth images and submissions from other users ( Bouman 2017 ) .", '5.2. First M87 Imaging Results': 'The /uniFB01 rst M87 imaging analysis used an early-release engineering data set ( ER4; Paper III ) . These data had a priori and network calibration applied but did not have calibrated relative RCP -LCP gains. Consequently, each team imaged the data using only parallel-hand products ( i.e., RCP -RCP or LCP -LCP ) to approximate total intensity. The April /uniF0A0 11 data set was selected for the /uniFB01 rst comparison, as it had the best coverage for the M87 / 3C 279 pair and the most stable a /uniF0A0 priori amplitude calibration ( especially for the LMT ) .  \nThe imaging teams worked on the data independently, without communication, for seven weeks, after which teams submitted images to the image comparison website using LCP data ( because the JCMT recorded LCP on April 11 ) . After ensuring image consistency through a variety of blind metrics ( including normalized cross-correlation, Equation ( 15 )) , we compared the independently reconstructed images from the four teams.  \nFigure 4 shows these /uniFB01 rst four images of M87. All four images show an asymmetric ring structure. For both RML teams and both CLEAN teams, the ring has a diameter of approximately 40 /uni03BC as, with brighter emission in the south. In contrast, the ring azimuthual pro /uniFB01 le, thickness, and brightness varies substantially among the images. Some of these differences are attributable to different assumptions about the total compact /uniFB02 ux density and systematic uncertainties ( see Table 2 ) .  \nThe initial blind imaging stage indicated that the image of M87 is dominated by a ∼ 40 /uni03BC as ring. The ring persists across the imaging methods. Next, we moved to a second, non-blind imaging stage that focused on exploring the space of acceptable images for each method. The independent team structure was only used for the /uniFB01 rst stage of imaging; the remainder of this Letter will categorize results by imaging methodology.', '6. Determining Imaging Parameters via Surveys on Synthetic Data': 'To explore the dependence of the reconstructed images on imaging assumptions and impartially determine a combination of /uniFB01 ducial imaging parameters, we introduced a second stage of image production and analysis: performing scripted parameter surveys for three imaging pipelines. To objectively evaluate the /uniFB01 delity of the images reconstructed by our surveys -i.e., to select imaging parameters that were independent of expert judgment -we performed these surveys on synthetic data from a suite of model images as well as on the M87 data. The synthetic data sets were designed to have properties that are similar to the EHT M87 visibility amplitudes ( e.g., prominent amplitude nulls ) . This suite of synthetic data allowed us to test the scripted reconstructions with knowledge of the corresponding ground truth images and, thereby, select /uniFB01 ducial imaging parameters for each method. These /uniFB01 ducial parameters were selected to perform well across a variety of source  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \nFigure 4. The /uniFB01 rst EHT images of M87, blindly reconstructed by four independent imaging teams using an early, engineering release of data from the April /uniF0A0 11 observations. These images all used a single polarization ( LCP ) rather than Stokes /uniF0A0 I , which is used in the remainder of this Letter. Images from Teams 1 and 2 used RML methods ( no restoring beam ) ; images from Teams 3 and 4 used CLEAN ( restored with a circular 20 /uni03BC as beam, shown in the lower right ) . The images all show similar morphology, although the reconstructions show signi /uniFB01 cant differences in brightness temperature because of different assumptions regarding the total compact /uniFB02 ux density ( see Table 2 ) and because restoring beams are applied only to CLEAN images.  \n<!-- image -->  \nTable 2 Image Properties and Data Consistency Metrics for the First M87 Images ( See Figure 4 )  \nNote. Data metrics are shown as originally computed on April /uniF0A0 11 data ( using 10 s averaged engineering data with LCP ) and using the data from the /uniFB01 rst EHT science release ( scan-averaged, Stokes I ) when 0%, 1% and 10% systematic error has been included. Teams /uniF0A0 2 -4 chose to exclude the intra-site baselines in their imaging. However, for consistency with our later χ 2 values computed from science release data, we include these baselines when computing χ 2 after adding an extended component to these images containing the missing /uniFB02 ux density.  \nstructures, including sources without the prominent ring observed in our images of M87.  \nWe emphasize that the ensemble of results from these parameter surveys do not correspond to a posterior distribution of reconstructed images. Our surveys are coarse-grained and do not completely explore the choices in the imaging process. Nonetheless, they identify regions of imaging parameter space that consistently produce faithful image reconstructions on synthetic data, and they help us identify which features of our reconstructions are consistent and which features vary with speci /uniFB01 c parameter choices.', '6.1. Synthetic Source Models and Data': 'To create a testing suite of synthetic data, we selected simple geometric models that have corresponding visibility amplitudes that are similar to those observed in M87 ( Figure 2 ) . The primary data properties that we used to de /uniFB01 ne similarity are ( 1 ) /uniF0A0 a large decrease in /uniFB02 ux density on baselines between 0 and 1G λ , indicating extended structure, ( 2 ) /uniF0A0 visibility nulls at ∼ 3.4 and ∼ 8.3 G λ , and ( 3 ) /uniF0A0 a high secondary peak between the nulls at ∼ 6G λ , which recovers ∼ 15% of the total compact /uniFB02 ux density.  \nWe selected four models with distinct compact morphologies that each reproduce these features of the M87 data. The four models are ( 1 ) a tapered ring with 44 /uni03BC as ring diameter, ( 2 ) a tapered crescent of the same diameter with its brightest point oriented directly south, ( 3 ) a tapered disk with 70 /uni03BC as diameter, and ( 4 ) two different circular Gaussian components separated by 32.3 /uni03BC as at a position angle of 292 ° . To ensure rough consistency and compatibility with the M87 parameters estimated in Section 4, we adopted a total compact /uniFB02 ux density of 0.6 Jy for all these simple geometric models. Note that none of the synthetic EHT data sets generated from these simple models reproduces all features seen in the M87 data. For example, the ring and disk models both have point symmetry, so all their closure phases are either 0 ° or 180 ° .  \nTo simulate the effects of a large-scale jet on our data ( which only signi /uniFB01 cantly affects intra-site visibilities ) , we added a three-component Gaussian model that approximates the inner M87 jet at 3 mm ( e.g., Kim et al. 2018a ) . The jet also has 0.6 Jy of total /uniFB02 ux density, giving a total image /uniFB02 ux density in each case ( compact + jet ) of 1.2 Jy. To produce non-closing systematic errors from polarimetric leakage, we also included linear polarization in each model. For additional details on these simulated models and data, see Appendix C.1. Figure 5 shows these model images.  \nWe generated synthetic data from each image using the eht-imaging software library. The synthetic data were produced with the baseline coverage and sensitivity of the EHT on all four days of the 2017 observations. Station-based errors were added in a Jones matrix formalism ( Thompson et al. 2017; see Appendix C.2 ) . To simulate a lack of absolute phase  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \nFigure 5. The four simple geometric models and synthetic data sets used in the parameter surveys ( see Appendix C for details ) . Top: linear scale images, highlighting the compact structure of the models. Middle: logarithmic scale images, highlighting the larger-scale jet added to each model image. Bottom: one realization of simulated visibility amplitudes corresponding to the April /uniF0A0 11 observations of M87. We indicate the conventions for cardinal direction and position angle used throughout this Letter on the upper-right panel. Note that east is oriented to the left, and position angles are de /uniFB01 ned east of north.  \n<!-- image -->  \ncalibration, a random station phase is adopted for each scan. Independent station-based gains were applied to each visiblity, with two components of each gain factor drawn from normal distributions: one term that was constant over the observation, and one that varied randomly among scans. We captured the increased uncertainty and variability in station gains at the LMT ( Section 4.1 ) by giving these gain terms larger variances. Time-independent polarimetric leakage terms were drawn from a complex Gaussian distribution with 5% standard deviation, motivated by previous estimates of the polarization leakage at EHT sites ( Johnson et al. 2015 ) . Identical gain, phase, and leakage contamination was applied to the high- and low-band data. Figure 5 shows four examples of visibility amplitudes generated using this procedure. Appendix C.2 provides full details about the synthetic data generation procedure.  \nFurthermore, to test whether or not our results are sensitive to speci /uniFB01 c choices made in the synthetic data generation, we compared our generated synthetic data and reconstructed images to results from another VLBI data simulator, MeqSilhouette+rPICARD ( Blecher et al. 2017; Janssen et al. 2019 ) . MeqSilhouette+rPICARD takes a more physical approach to mm-VLBI synthetic data generation ( with added corruption based on, e.g., measured weather parameters and  \nantenna pointing offsets ) and it uses the full CASA -based EHT reduction pipeline for calibration. Despite the differences in approach, both synthetic data pipelines yield comparable results ( see Appendix C.3 and Figure 31 ) . We use ehtimaging for all further synthetic data in this Letter.', '6.2. Imaging Pipelines and Parameter Survey Space': 'To survey the space of imaging parameters relevant to each imaging software package and to test their effects on reconstructions of simulated data, we designed scripted imaging pipelines in three software packages: DIFMAP , eht-imaging , and SMILI . Each pipeline has some choices that are /uniFB01 xed ( e.g., the convergence criterion, the pixel size, etc. ) but takes additional parameters ( e.g., the regularizer weights, the total compact /uniFB02 ux density ) as arguments. We then reconstructed images from all M87 and synthetic data sets using all possible parameter combinations on a coarse grid in the space of these parameters. We chose large ranges for each parameter, deliberately including values that we expected to produce poor reconstructions.  \nThe parameter choices and the values surveyed vary among the three pipelines, and the pipelines also differ in which frequency  \nbands were used. However, a consistent feature of the three pipelines is that none of them do any data /uniFB02 agging, even of intrabaselines. 116 We now describe each pipeline in detail.', '6.2.1. DIFMAP': "While CLEAN was implemented by hand in our /uniFB01 rst stage of analysis ( Section 4 ) , a scripted version of CLEAN was developed in DIFMAP to carry out the parameter search in this second stage. This script has /uniFB01 ve parameters taken as arguments: the total assumed compact component /uniFB02 ux density, the stop condition, the relative weight correction factor for ALMA in the self-calibration process, the diameter of the cleaning region, and a parameter controlling how the ( u , v ) grid weights are scaled by the visibility errors. CLEAN imaging with DIFMAP was performed separately for the low- and highband data. The DIFMAP images we present in this Letter are derived from the low-band data averaged to 10 s in time.  \nThe DIFMAP script begins by gridding the visibilities in the ( u , v ) plane, after which the direct Fourier transform, or ' dirty image, ' and the point-spread function, or ' dirty beam, ' are computed. A uniform ( u , v ) density weighting function is used for gridding the visibilities. This choice gives equal weight to each gridded spatial frequency containing a measurement, irrespective of sample density. The visibility weights in ( u , v ) gridding can be further multiplied by σ κ ,where σ is the thermal noise, with our parameter survey exploring 2 k =-, -1, and 0.  \nIn addition to the ( u , v ) weighting, one or more antennas can be downweighted in the self-calibration process. When /uniFB01 tting for complex antenna gains, the self-calibration algorithm of DIFMAP uses weights that are proportional to σ -2 . In the case of EHT data, it is particularly helpful to downweight data on baselines to phased ALMA -which has a much higher S / N than other stations -so that these ALMA baselines do not dominate the corrections to phases and amplitudes. In the parameter survey, we have scaled the ALMA weights by a factor ranging from 0.01 to 1.0. All other station weights were kept at 1.0 ( i.e., at their original values ) .  \nA disk-shaped set of CLEAN windows, or imaging mask, aligned with the previously found geometrical center of the underlying emission structure ( e.g., the ring for M87 ) and with a speci /uniFB01 ed diameter is then supplied; these windows de /uniFB01 ne the area on the image where the CLEAN algorithm searches for point sources. We limited the cleaning windows to image only the compact structure ( /uni0084 100 /uni03BC as ) in order to prevent CLEAN from adding larger-scale emission features that are poorly constrained by the lack of short EHT baselines. The emission structure common to all early trials ( Section 5 ) was used as a guide for setting the diskshaped window for scripted reconstructions, but different diameters from 40 to 100 /uni03BC as were tested in the parameter search.  \nTwo main stopping criteria are implemented in the DIFMAP script to control the amount of cleaning; imaging may stop ( 1 ) when the desired total compact /uniFB02 ux density is reached, or ( 2 ) when the relative decrease in the rms of the residual image over the image noise estimated from the visibilities dips below a given threshold. In the parameter survey, two cases were tested: ( a ) criterion 1 alone, and ( b ) a combination of criteria 1 and 2, such that cleaning will stop whenever one of the two  \ncriteria is reached /uniFB01 rst. In both cases, the assumed total compact /uniFB02 ux density was varied from 0.5 to 0.8 Jy. The threshold for the relative decrease in the residual image rms in the criterion 2 was always /uniFB01 xed to 10 -4 .  \nAfter the script /uniFB01 nishes a round of cleaning, it introduces a large ( /uniF089 2 mas ) Gaussian component that recovers the extended missing /uniFB02 ux density and does not affect the compact emission structure. By /uniFB01 tting the most extended structure, the Gaussian component also acts as an ' anchor ' for the intra-site baselines, minimizing the variations in the antenna gains during the amplitude self-calibration loops.  \nThe initial image, produced by iterative cleaning and phase self-calibration, is used to amplitude self-calibrate the LMT ( with a solution interval equal to the scan length ) , which is known to have poor a priori amplitude calibration ( see Section 4.1 ) . In this step, the gains of the other antennas in the array are /uniFB01 xed. The initial model is then removed and further iterations of cleaning loops are performed until a new model that recovers all the compact /uniFB02 ux density is obtained. This round is followed by self-calibration at all sites in both phase and amplitude. The loop of cleaning and self-calibration is then repeated with progressively shorter solution intervals in the amplitude self-calibration step, from an initial solution equal to the whole duration of the observations down to a minimum solution interval of 10 s. For consistency with the outputs of the eht-imaging and SMILI pipelines, the /uniFB01 nal CLEAN image omits the large Gaussian component.  \nFigure 6 shows a 2D slice of the DIFMAP parameter survey results from varying the weights on ALMA data in the selfcalibration step and the diameter of the circular mask of CLEAN windows on synthetic data. We show the results for synthetic data of a crescent model and for M87 data from April 11. It is important to note that although DIFMAP images are displayed using a restoring beam of 20 /uni03BC as ( FWHMshown on /uniFB01 gures ) ,asis standard in CLEAN imaging, the underlying δ -function CLEAN components are used for computing χ 2 so as to avoid modifying the reconstructed amplitudes through Gaussian convolution.", '6.2.2. eht-imaging': 'The eht-imaging parameter survey used a template imaging script that takes seven parameters as arguments. These are the total assumed compact component /uniFB02 ux density, the fractional systematic error on the measured visibilities, the FWHM of the circular Gaussian initial image ( also used as a prior for maximum entropy regularization ) , and the relative weights of four regularization terms: MEM, the ℓ 1 norm ( ℓ 1 ) , total variation ( TV ) , and total squared variation ( TSV ) . For details on each of these regularization terms, see Appendix A.  \nThe eht-imaging script begins by preparing the data for imaging. After scan-averaging the data and combining the high- and low-band data, it rescales the amplitude of zero baseline measurements to account for the unresolved extended components. It then performs an initial self-calibration of the LMT by using only the LMT -SMT baseline and adjusting only the LMT amplitude gains to match a Gaussian image with a total /uniFB02 ux density of 0.6 Jy and a size of 60 /uni03BC as, based on the constraints determined in Section 4.  \nThe script also incorporates additional systematic uncertainty into the error budget during imaging. Namely, for each visibility, a speci /uniFB01 ed fraction of the visibility amplitude is added in quadrature to the thermal noise; this in /uniFB02 ated thermal noise is then used to estimate uncertainties on all data products,  \nFigure 6. Selection of the DIFMAP ( CLEAN ) parameter survey results on real and synthetic data with April /uniF0A0 11 EHT baseline coverage. A 2D slice of the 5D parameter space is displayed, corresponding to different diameters of the circular mask and the data weight on ALMA in self-calibration. All other parameters are kept constant ( Compact Flux /uniF0A0 = /uniF0A0 0.5 Jy, κ = -1, Stop Condition /uniF0A0 = /uniF0A0 Flux Reached ) . The left panel shows results of the parameter search on the Crescent synthetic data, while the right panel shows reconstructions for the same parameters on M87 data. Images that meet the threshold for the Top Set are outlined in green ( see Section 6.3.1 ) .  \n<!-- image -->  \nthereby accounting for small, non-closing systematic errors from, e.g., polarimetric leakage. We surveyed values of 0%, 1%, 2%, and 5% in computing this additional systematic uncertainty. Note that these levels of systematic uncertainty do not necessarily correspond with the estimated level of nonclosing errors presented in Paper III; however, additional systematic uncertainty beyond the a /uniF0A0 priori estimate can help in avoiding over /uniFB01 tting to a few sensitive baselines ( e.g., those to ALMA ) during image optimization.  \nThe /uniFB01 nal pre-imaging step is to inverse taper the visibility amplitudes. Speci /uniFB01 cally, we divide the visibility amplitudes and their estimated error by the Fourier transform of a 5 /uni03BC as FWHM Gaussian, thereby raising the amplitudes on long baselines while preserving the S / N. At the end of the script, after producing an image from the inverse-tapered data, we convolve the reconstructed image with a 5 /uni03BC as FWHM Gaussian to ensure that it will /uniFB01 t the original, untapered data. This procedure enforces a limiting angular resolution of 5 /uni03BC as in the /uniFB01 nal reconstructed images.  \nAfter the data have been prepared, the eht-imaging script proceeds to imaging, using a 64 /uniF0A0 × /uniF0A0 64 pixel grid with a 128 /uni03BC as FOV. The script alternates between solving for the optimal image under the current data constraints and self-calibrating the complex station gains ( Equation ( 7 )) . The /uniFB01 rst imaging iteration uses only visibility amplitudes and closure quantities. Because the amplitudes are uncalibrated in this iteration, systematic error tolerance is added on top of the ( potentially in /uniFB02 ated ) thermal uncertainty for the visibility amplitudes, using the systematic error budget provided in Paper III for all sites  \nexcept LMT, which has its budget increased by 15%. After phase self-calibration, complex visibilities are used for optimization in combination with closure quantities.  \nBetween imaging iterations, the image is convolved with a circular Gaussian with a FWHM corresponding to the nominal array resolution, u 25 as max l m » /uni2223 /uni2223 . This procedure aids convergence to a global minimum by moving intermediate images away from local minima. In each iteration, we increase the weight on data terms relative to the regularizers and reduce the systematic noise tolerance for amplitude calibration uncertainties. The /uniFB01 nal reconstructed image is convolved with the same 5 /uni03BC as Gaussian that was used for the inverse taper, ensuring consistency with the original data while imposing a constraint that no features in the reconstructed image can be /uniFB01 ner than 5 /uni03BC as.  \nFigure 7 shows a 2D slice of the eht-imaging parameter survey results from varying the weights on the MEM and the TV regularizers. In this /uniFB01 gure, we show the results for both synthetic data from a crescent model and for EHT M87 data from April /uniF0A0 11.', '6.2.3. SMILI': "In the SMILI imaging pipeline, we reconstructed images using low-band EHT data and utilized weightedℓ 1 ( ℓ 1 w ) , TV, and TSV regularizers ( see Appendix A ) . ℓ 1 w favors sparsity in the image domain, using a circular Gaussian prior image as a ' soft mask, ' which increasingly penalizes pixel intensities farther from the image center. In contrast, TV favors sparsity in the gradient domain ( leading to piecewise-smooth images ) , and TSV favors overall image smoothness. The SMILI search  \nFigure 7. Selection of the eht-imaging ( RML ) parameter survey results on real and synthetic data with April /uniF0A0 11 EHT baseline coverage. A 2D slice of the 7D parameter space is displayed, corresponding to different weights on the MEM and TV regularizers. All other parameters are kept constant ( Compact /uniF0A0 Flux /uniF0A0 = /uniF0A0 0.6 Jy, Initial / MEM /uniF0A0 FWHM /uniF0A0 = /uniF0A0 40 /uni03BC as, Systematic /uniF0A0 Error /uniF0A0 = /uniF0A0 1%, TSV /uniF0A0 = /uniF0A0 0, and ℓ 1 /uniF0A0 = /uniF0A0 0 ) . The left panel shows results of the parameter search on the Crescent synthetic data, while the right panel shows reconstructions for the same parameters on M87 data. Images that meet the threshold for the Top Set are outlined in green. Note that the upper-left reconstruction has no regularization; it is produced by enforcing only image positivity and a constrained FOV.  \n<!-- image -->  \nadopts an FOV of 128 /uni03BC as ( matching the eht-imaging pipeline ) . The survey included /uniFB01 ve values for the FWHM of the prior image for ℓ 1 w , ranging from 40 to 80 /uni03BC as. The pipeline also surveyed /uniFB01 ve values for the total compact /uniFB02 ux density, between 0.4 and 0.8 Jy.  \nPrior to imaging, the script coherently averages the visibilities over scans and then pre-calibrates the LMT gains using an assumed Gaussian model, as in the procedure described above for the eht-imaging script. The script then adds systematic uncertainties in quadrature to the thermal noise of complex visibilities to account for small, non-closing errors. The survey included options for 0% ( no systematic error ) and 1% systematic error, which correspond to the fraction of each visibility amplitude that is added in quadrature to the reported thermal noise σ values.  \nSimilar to the eht-imaging and DIFMAP scripts, the SMILI imaging procedure is iterative, with four stages of imaging and self-calibration. Reconstructions at each stage begin with a circular Gaussian with an FWHM of 20 /uni03BC as as the initial image. After the /uniFB01 rst imaging attempt in each stage, subsequent initializations use the previously obtained image convolved with the initial Gaussian. We perform three imaging cycles for each self-calibration stage, self-calibrating only to the /uniFB01 nal reconstructed image in each of the three cycles.  \nIn the /uniFB01 rst two self-calibration stages, the imaging step uses closure data ( closure amplitudes and phases ) and visibility amplitudes. To account for errors in the amplitude calibration, during these stages the visibility amplitudes have additional 5%  \nand 30% uncertainties added in quadrature to the thermal noise ( on top of the systematic error tolerance noted above ) for nonLMT and LMT baselines, respectively. For the /uniFB01 nal two selfcalibration stages, the imaging uses closure data and complex visibilities.  \nFigure 8 shows a 2D slice of the SMILI parameter survey results from varying the weights on the TSV regularizer and the prior radius. We show the results for synthetic data from a crescent model and for M87 data from April /uniF0A0 11. Although the SMILI and eht-imaging pipelines have a number of similarities, these pipelines use different regularizers, calibrate and constrain data differently, and were implemented using independent software libraries.", '6.3. Selecting Imaging Parameters': "We derived two principal outputs from each parameter survey: a ' Top Set ' of parameter combinations that each produce acceptable images on the synthetic data, and a single combination of best-performing, ' /uniFB01 ducial ' parameters for each pipeline. In this section, we describe our ranking procedure to /uniFB01 nd acceptable parameter combinations, the possible limitations of this procedure, and the resulting acceptable parameter combinations for each imaging pipeline.", '6.3.1. Ranking Imaging Parameters': "To determine the Top Sets, we designed a two-step selection procedure. The /uniFB01 rst step requires that all reconstructions of M87  \nFigure 8. Selection of the SMILI ( RML ) parameter survey results on real and synthetic data with April /uniF0A0 11 EHT baseline coverage. A 2D slice of the 6D parameter space is displayed, corresponding to varying the diameter of the soft mask region and the weight on the TSV regularizer. All other parameters are kept constant ( Compact Flux /uniF0A0 = /uniF0A0 0.5 Jy, Systematic Error /uniF0A0 = /uniF0A0 0%, TV /uniF0A0 = /uniF0A0 0, ℓ 1 1 W = ) . The left panel shows results of the parameter search on the Crescent synthetic data, while the right panel shows reconstructions for the same parameters on M87 data. Images that meet the threshold for the Top Set are outlined in green.  \n<!-- image -->  \nwith a given combination of parameters are consistent with their underlying data; the second step requires that the reconstructions of simulated data are suf /uniFB01 ciently similar to their corresponding ground truth images. Parameter combinations that meet both criteria de /uniFB01 ne the Top Set, and the parameter combination in this set that gives the best performance on the simulated data is used to de /uniFB01 ne the /uniFB01 ducial images of M87 for each day.  \nEach parameter survey spans a large range of parameter space, including combinations that are poorly adapted to the M87 data ( e.g., very large regularizer weights in the RML methods ) .To /uniFB01 lter parameter combinations that lead to /uniFB01 nal images with poor /uniFB01 ts to the data, we calculate χ 2 values on closure phase and log closure amplitudes using scan-averaged low-band data. For the χ 2 computation, we only construct closure quantities from averaged visibilities that have synchronized start times. Any combinations with 2 2  c on any day of EHT observations are excluded from the Top Set. However, the DIFMAP imaging pipeline uses 10 s averages ( while we estimate χ 2 values using scan averages ) , and the pipeline downweights the extremely high S / N baselines to ALMA during the iterative self-calibration process. Consequently, DIFMAP reconstructions tend to fail this strict χ 2 requirement unless a systematic uncertainty is included in the error budget. As a result, before performing this initial χ 2 /uniF0A0 < /uniF0A0 2 thresholding on DIFMAP reconstructions, we include 10% systematic uncertainty in quadrature to the thermal uncertainty of every visibility.  \nNext, to determine the /uniFB01 delity of the different reconstructions on synthetic data, we compute the normalized cross-correlation ρ NX between each surviving image and the corresponding ground  \nFigure 9. Normalized cross-correlation, ρ NX, of the four simulated images compared before and after convolution with a circular Gaussian kernel with FWHM α . The vertical dashed line shows the nominal diffraction-limited resolution of the EHT ( see Table 1 ) . For each source, we invert these curves to obtain from a given ρ NX an equivalent blurring kernel size α . The Top Set reconstructions are de /uniFB01 ned as those with this α smaller than the nominal EHT resolution.  \n<!-- image -->  \ntruth image. For a given reconstruction X and its ground truth image Y , the normalized cross-correlation is  \nX Y N X X Y Y , 1 , 15 i i i X Y NX å r ss = -á ñ -á ñ ( ) ( )( ) ( )  \nFigure 10. Cross-validation of the imaging parameter selection procedure. In each of the left four columns, we show reconstructed images for the simple geometric source models. These reconstructions do not use the /uniFB01 ducial imaging parameters identi /uniFB01 ed by the full training set; instead, we selected the imaging parameters for each geometric source model after excluding that particular model from the parameter selection process. For example, in the disk reconstructions, the parameters were selected by assessing reconstructions of only the ring, crescent, and double source models. Thus, the selected parameters vary among these four columns, but we can verify that the training sets do not overly constrain the outcomes. In the /uniFB01 fth column, we show reconstructions of a GRMHD snapshot ( Paper V ) using the /uniFB01 ducial M87 parameters selected from all four geometric models. That is, the script and parameters used to produce these GRMHD image reconstructions are identical to those used to produce our /uniFB01 ducial M87 images ( shown in Figure 11 ) . Because the GRMHD snapshot has a substantially higher peak brightness than the reconstructions, its column has been scaled to the peak brightness of the eht-imaging reconstruction.  \n<!-- image -->  \nwhere the sum is over all N pixels in the two images, X á ñ and Y á ñ are the mean pixel values in the images, and σ X and σ Y are the standard deviations of the pixel values in the two images. The normalized cross-correlation quanti /uniFB01 es image similarity, where a value of 1 implies perfect correlation between images, 0 implies no correlation, and -1 implies perfect anticorrelation. Because the absolute visibility phase is unknown and, thus, the image centroid for reconstructions is not /uniFB01 xed, we compute ρ NX for each possible discrete shift between the images and select the largest value. Before computing ρ NX for DIFMAP images, the CLEAN components are convolved with a restoring beam with FWHM of 20 /uni03BC as.  \nWe choose the median value of ρ NX across all days as the representative value for a given combination of imaging parameters on a particular model. Taking the median allows us to incorporate days with different baseline coverage in the scoring procedure while preventing the poorest reconstructions ( e.g., those on April 10 ) from dominating the outcome.  \nWhile ρ NX is a useful measure of relative agreement among images, it does not provide an absolute score of reconstruction quality. For instance, Figure 9 shows that convolving the four geometric source models with the same circular Gaussian of FWHM α produces ρ NX values ( computed against the original image ) that are different depending on the underlying model. Thus, to better standardize the interpretation of ρ NX for  \nreconstructions of different underlying images, we invert the curves in Figure 9 to compute, for each value of ρ NX, an effective kernel size α . That is, we quantify the /uniFB01 delity of a reconstruction in terms of an effective blurring kernel FWHM α that would give the same ρ NX as we compute from the reconstruction.  \nWe compute a characteristic value α char for each combination of imaging parameters by taking the mean α across the four sources, and we then rank the parameters according to this characteristic value. The top-ranked combination of parameters is then denoted as the ' /uniFB01 ducial parameters ' of that method. We also select a Top Set of parameter combinations that obtain an α char smaller than the nominal EHT resolution. That is, the Top Set includes all imaging parameter combinations with u 1 char max a < /uni2223 /uni2223 . A selection of the Top Set reconstructions is highlighted in green for the DIFMAP , eht-imaging , and SMILI parameter searches in Figures 6 -8, respectively.", '6.3.2. Parameter Selection Considerations': 'In Figure 10, we test whether or not our parameter selection procedure leads to reconstructed images that are unduly similar to the training set. In this test, we repeat the parameter selection process after withholding a speci /uniFB01 ed geometric source model from the training set. For example, in Figure 10, the parameter combinations used to produce the images of the disk model were determined by selecting the best-performing parameters on only the ring, crescent, and double. Despite excluding the disk from the training set, all three pipelines produce an image with a disk morphology that does not resemble any of the other images in the training set.  \nA common feature of all four geometric models in the training set is that they restrict the compact emission to a region spanning /uniF088 70 /uni03BC as. One possible concern is that the resultant Top Set parameter combinations may then fail to reconstruct images of sources that have compact components extending over larger regions. However, such bright compact components would introduce variations over time in the visibilities that could not be /uniFB01 t with an overly restricted FOV. Thus, our requirement that Top Set reconstructions of M87 have χ 2 /uniF0A0 < /uniF0A0 2 ensures that they are not missing bright extended components.  \nOther common features in all our training set models are the lack of small-scale structure ( the /uniFB01 nest features in each image is /uni0085 10 /uni03BC as ) and a /uniFB01 xed total compact /uniFB02 ux density of 0.6. We test that our /uniFB01 ducial scripts produce accurate results for an image with small-scale structure using using a GRMHD snapshot image ( see Figure 10 ) , and we explore systematic uncertainties related to the total compact /uniFB02 ux density in Appendix H.', '6.3.3. Results of Parameter Selection': 'In Table 3, we summarize the speci /uniFB01 cations and outcomes of the three parameter surveys. In this table, for each parameter value, we list the fraction of parameter combinations in the Top Set that include that value. We also provide the total number of parameter combinations surveyed, the number of these combinations in the Top Set, and the /uniFB01 ducial imaging parameters for each pipeline. Table 4 quanti /uniFB01 es the /uniFB01 delity of the /uniFB01 ducial and Top Set images reconstructed from simulated data by reporting the equivalent blurring kernel α obtained with each method.  \nThe distributions in Table 3 show that the parameter survey eliminates regions of parameter space from consideration based  \nTable 3 Parameters and Their Surveyed Values for the DIFMAP ( Top ) , eht-imaging ( Middle ) , and SMILI ( Bottom ) Pipelines  \nNote. For instance, the eht-imaging survey produced 5 /uniF0A0 × /uniF0A0 3 /uniF0A0 × /uniF0A0 4 /uniF0A0 × /uniF0A0 5 4 /uniF0A0 = 37,500 parameter combinations, and 1572 of these passed the criteria for inclusion in the Top Set. Below each parameter value we specify the fraction of the Top Set parameter combinations that include that value ( see Section 6.3.1 ) . Boxed parameters are those corresponding to the /uniFB01 ducial images. Note that the /uniFB01 ducial parameters are determined by identifying the parameter combination that jointly performs best on all the synthetic data sets; the /uniFB01 ducial parameters do not necessarily correspond with the parameters that have the largest share in the Top Set.  \non our selection criteria. For instance, in the eht-imaging results, reconstructions using large hyperparameters ( > 10 ) on the TV and TSV regularizers or those using a large initial Gaussian ( FWHM /uniF0A0 > /uniF0A0 50 /uni03BC as ) perform poorly regardless of the other parameter values used. These parameters are not represented in the Top Set.  \nIn addition, differences in the procedures for each pipeline lead to different /uniFB01 ducial parameters. In some cases, these are  \nFigure 11. Fiducial images of M87 on all four observed days from each of the three imaging pipelines. CLEAN images ( from DIFMAP ) are shown after convolution with a 20 /uni03BC as beam; eht-imaging and SMILI results have no restoring beam applied. Different selected /uniFB01 ducial imaging parameters ( e.g., compact /uniFB02 ux ) result in different peak brightness temperatures for each method, as indicated by the unique color bars for each row.  \n<!-- image -->  \nTable 4  \nEquivalent Blurring Values α from the Fiducial and Top Set Images from Each Pipeline on Each of the Four Simple Geometric Models in the Training Set  \nNote. For each pipeline, we give α for the /uniFB01 ducial image, as well as the mean and standard deviation of α across the Top Set.  \nexpected from differences in the underlying methods. For instance, the prior and mask sizes are treated differently in different scripts: a hard mask was used for DIFMAP , a softer mask for SMILI , and only weak constraints on compactness were used for eht-imaging . However, other choices are unexpected. For example, the top-ranked parameter combinations ( i.e., the /uniFB01 ducial parameters ) for the DIFMAP and SMILI pipelines use a total compact /uniFB02 ux density of 0.5 Jy, favoring this value over the true total compact /uniFB02 ux density of 0.6 Jy. However, these preferences are mild, and the Top /uniF0A0 Set from each pipeline includes images with all values of total compact /uniFB02 ux density surveyed, from 0.4 to 0.8 Jy ( see Table 3 ) .We explore the dependence of the assumed compact /uniFB02 ux density in Appendix H.', '7. Fiducial M87 Images: Properties and Uncertainties': 'Having determined /uniFB01 ducial imaging parameters via tests on synthetic data ( Section 6.3.1 ) , we now show the corresponding /uniFB01 ducial M87 images from each imaging method for each day of the 2017 EHT observations. We also assess the consistency of the /uniFB01 ducial images with the data, discuss other visibilitydomain properties of these images, and estimate image uncertainties resulting from imaging parameter choices.', '7.1. Fiducial M87 Images': 'Figure 11 shows the /uniFB01 ducial images from each day of EHT observations and each imaging method. These /uniFB01 ducial images are broadly consistent; all twelve images have a prominent,  \nTable 5 Closure Quantity χ 2 Values for the Fiducial M87 Images and χ 2 Statistics ( Mean and Standard Deviation ) for Top Set Images  \nNote. For each, values are shown after including systematic uncertainties of 0%, 1%, and 10%. Inclusion in the top set requires χ 2 /uniF0A0 < /uniF0A0 2 for both closure phases and amplitudes with 0% systematic uncertainty for eht-imaging and SMILI , and with 10% systematic uncertainty for DIFMAP . See Section 2.1 for de /uniFB01 nitions and limitations of these metrics.  \nasymmetric ring feature, approximately 40 /uni03BC as in diameter with enhanced brightness toward the south. Each image has a prominent depression in the center of the ring and a peak brightness temperature along the ring of ∼ 10 10 K. Because they are restored with a 20 /uni03BC as beam, the DIFMAP images have broader structure, leading to a weaker central depression and lower peak brightness temperatures. For all pipelines, the /uniFB01 nite resolution of the reconstructions results in brightness temperatures that are only lower limits on the maximum intrinsic brightness temperature.  \nAs expected from tests on synthetic data ( Section 6 ) , the details of the reconstructions differ between the imaging methods. Although the RML images ( eht-imaging and SMILI ) are generally more similar to each other than to the DIFMAP images, SMILI reconstructions are fainter inside and outside the ring than the eht-imaging reconstructions from the strong sparsity regularization. Based on the results of the parameter survey, M87 data was reconstructed with a total compact /uniFB02 ux density of 0.5 Jy for DIFMAP and SMILI , and with 0.6 Jy for eht-imaging . While 0.5 Jy is slightly below the range of acceptable values for F cpct identi /uniFB01 ed in Section 4, the constraints presented in Section 4 made strict assumptions about the gains that are not rigidly enforced in imaging, and our parameter surveys using synthetic data demonstrate that some imaging methods produce higher /uniFB01 delity ( i.e., higher ρ NX ) images when reconstructing an image with less compact /uniFB02 ux density than the true value. We explore the effects of the assumed compact /uniFB02 ux density further in Appendix H.  \nIn Table 5, we provide χ 2 values for the 12 /uniFB01 ducial images as well as χ 2 statistics across the Top Sets. We compute these values after including varying amounts of systematic uncertainty ( 0%, 1%, and 10% of visibility amplitudes ) added in quadrature to the thermal uncertainty. This additional uncertainty accounts for possible systematic errors ( e.g., from polarimetric leakage ) and also ensures that the statistics quantify overall image-data consistency without being overwhelmed by the extremely sensitive ALMA baselines. For the largest value, 10%, the systematic uncertainty allows modest tolerance for source structure that is absent from the reconstructed images ( e.g., from extended structure in the M87 jet outside the masked region or the FOV of the reconstructed image ) .  \nFigures 12 and 13 compare self-calibrated visibility amplitudes and closure phases on three selected triangles with the corresponding model visibility amplitudes and closure phases computed from each April 11 /uniFB01 ducial image. As indicated in the χ 2 values listed in Table 5, the two RML methods produce /uniFB01 ducial images that are consistent with the data to within roughly their thermal uncertainties, and the χ 2 values have little scatter across the Top Sets. The DIFMAP /uniFB01 ducial image also shows agreement with the data, but requires modest systematic tolerance for comparable image-data consistency ( because of the different time averaging used in DIFMAP imaging the ALMA downweighting ) , and the χ 2 values show more scatter across the Top Set. All three /uniFB01 ducial images /uniFB01 t the nulls in visibility amplitudes at 3.4 and 8.3 G λ .  \nFigure 12. Top: visibility amplitudes of self-calibrated data on April 11 as a function of baseline length compared with corresponding values for the three /uniFB01 ducial images. In each panel, visibilities are plotted after self-calibrating the data to the indicated /uniFB01 ducial image; thus, the visibility amplitudes differ across the three panels. For each image, we add an additional large-scale Gaussian component before self-calibration to account for the excess /uniFB02 ux density seen on intra-site baselines ( see Section 4.3 ) . Bottom: visibility amplitude residuals for each /uniFB01 ducial image, after normalizing by the thermal noise with 1% systematic uncertainty added in quadrature.  \n<!-- image -->  \nFigure 13. Closure phases plotted as function of GMST on three selected triangles from the April 11 observations. The solid lines indicate the corresponding closure phase curves from the /uniFB01 ducial images produced by the three pipelines.  \n<!-- image -->  \nOur images of M87 do not show extended emission outside the observed ring, although a jet is prominent in VLBI images of M87 at longer wavelengths. Likewise, while the geometric models in the training set ( Section 6.1 ) include emission from a simple extended jet model, none of the Top Set reconstructions from simulated EHT data recover this extended component nor are their recovered compact structure affected by the missing jet components. Furthermore, our /uniFB01 ducial images in all three pipelines do not recover the more realistic faint jet emission in the GRMHD simulation reconstructions in Figure 10. All of these results reinforce the expectation ( see Section 4.3 ) that the current EHT array lacks the short-baseline coverage and dynamic range necessary to recover a faint, extended jet.  \nIn Figure 14,weshowthe /uniFB01 ducial images from the three pipelines on April /uniF0A0 11 after blurring each of them separately to achieve a common, conservative resolution. Speci /uniFB01 cally, we restore the DIFMAP image with the same 20 /uni03BC as beam used throughout this Letter. For eht-imaging and SMILI , we use the circular Gaussian convolution kernels that maximize the normalized cross-correlation between the restored RML image and the restored DIFMAP image ( 17.1 /uni03BC as for eht-imaging ; 18.6 /uni03BC as for SMILI ) . The three blurred images have a consistent ring diameter ( ≈ 40 /uni03BC as ) and the same overall asymmetry ( oriented south ) . Figure 15 shows the simple average of the three equivalently blurred /uniFB01 ducial images for each day. We adopt the images in  \nFigure 14. Fiducial images of M87 on April /uniF0A0 11 from our three separate imaging pipelines after restoring each to an equivalent resolution. The eht-imaging and SMILI images have been restored with 17.1 and 18.6 /uni03BC as FWHM Gaussian beams, respectively, to match the resolution of the DIFMAP reconstruction restored with a 20 /uni03BC as beam.  \n<!-- image -->  \nFigure 16 shows the visibility amplitude and phase for each of the three April /uniF0A0 11 /uniFB01 ducial images as a function of vector baseline. Note that no restoring beam is required for CLEAN in this visibility-domain analysis. Each image produces nulls in the visibility amplitude near the SMA -SMT baseline, consistent with the observed amplitudes ( see Figure 2 ) . The visibility phase shows rapid swings at these nulls. The visibilities of the images from the different pipelines are most similar near the EHT measurements, as expected. On longer baselines than those sampled by the EHT, the DIFMAP image produces notably higher visibility amplitudes than those of the  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \nFigure 15. Averages of the three /uniFB01 ducial images of M87 for each of the four observed days after restoring each to an equivalent resolution, as in Figure 14. The indicated beam is 20 /uni03BC as ( i.e., that of DIFMAP , which is always the largest of the three individual beams ) .  \n<!-- image -->  \nFigure 15 as a conservative representation of our /uniFB01 nal M87 imaging results.  \nThe /uniFB01 ducial images from each pipeline ( Figure 11 ) , as well as the conservative, blurred, pipeline-averaged images ( Figure 15 ) provide some evidence for evolution in the ring structure between April 5, 6 and 10, 11. We discuss this evolution in more detail in Appendix E ( Figure 33 ) . Some change in the image structure between April 5, 6 and 10, 11 is necessitated by the variations seen in the underlying closure phases ( Paper III ) .We /uniFB01 nd more variation in the image pairs that are separated by larger intervals, suggesting that these variations are intrinsic. However, we cannot unambiguously associate the observed variability with coherent evolution of a speci /uniFB01 c image feature.  \neht-imaging and SMILI images, as expected from the fact that the DIFMAP image is fundamentally a collection of point sources.', '7.2. Image Uncertainties': "Measuring the variation in images produced in a parameter survey Top Set allows us to evaluate image uncertainties due to the explored imaging choices. Figure 17 shows uncertainties related to imaging assumptions from the largest Top Set ( that of the eht-imaging parameter survey ) on April /uniF0A0 11 data.  \nReconstructed image uncertainties are concentrated in the regions with enhanced brightness, notably in the three ' knots ' in the lower half of the ring ( Figure 17; top panel ) . These are also the regions that show the most variation among different imaging methods ( Appendix I compares their azimuthal pro /uniFB01 les ) . Visibility-domain modeling provides another method to assess image structure. In Paper VI,weexplore /uniFB01 tting simple crescent models to the data. For instance, a crescent with a brightness gradient and blurring reproduces the north -south asymmetry in images without additional azimuthal structure ( the ' blurred and slashed with LSG ' crescent of Paper VI ) . However, this model gives CP 2 c between 3.2 and 11.5 and log CA 2 c between 2.2 and 6.6 for different days and bands when assuming 0% systematic error ( compare with Table 5 ) . Adding additional degrees of freedom in the form of three elliptical Gaussian components to the crescent  \nFigure 16. Visibility-domain representation of the three /uniFB01 ducial pipeline images for April /uniF0A0 11. The top panels show visibility amplitude as a function of vector baseline; the bottom panels show visibility phase as a function of vector baseline. Amplitude contours are spaced logarithmically, with three contours per decade. Thicker contours correspond to 10 and 100 mJy; the lowest contour is at 10 1 / 3 /uniF0A0 ≈ /uniF0A0 2.2 mJy, although the images produce nulls in the visibility domain. Plotted points show the April /uniF0A0 11 EHT coverage.  \n<!-- image -->  \nyields better /uniFB01 ts, similar to those found with imaging. These Gaussian components sometimes reside near or on the crescent, affecting the azimuthal structure ( see Appendix /uniF0A0 B.2.2inPaperVI ) . Thus, our current image reconstructions and visibility-domain analyses are not able to con /uniFB01 rm or reject the reality of the ' knots ' seen in images, and these features should be interpreted with caution.  \nIn the visibility domain, the uncertainties are concentrated at the origin, re /uniFB02 ecting the spread of total compact /uniFB02 ux density choices present in the Top /uniF0A0 Set ( Figure 17; bottom panel ) . After normalizing the total /uniFB02 ux densities of the Top Set images, the visibility domain uncertainties at the origin are zero; the uncertainties peak on baselines at ∼ 2G λ , then fall gradually with increasing baseline length. The fractional uncertainties are the most straightforward to interpret, as they are low ( /uniF088 10% ) in regions near EHT coverage and grow to values near or exceeding unity for baselines that are separated by more than a few G λ from an EHT measurement.  \nThe strong correlation between the visibility-domain fractional uncertainties and gaps in EHT coverage suggests that imaging uncertainties are dominated by imaging parameter choices rather than by thermal or systematic noise. That is, the /uniFB01 ducial image /uniFB01 delity is primarily limited by the sparse EHT baseline coverage, rather than by sensitivity or calibration accuracy. For example, the standard deviation of visibility  \namplitudes across the Top /uniF0A0 Set for M87 is ∼ 50 mJy on baselines at ∼ 2G λ , falling to ∼ 10 mJy at ∼ 8G λ . In contrast, the median thermal noise of non-ALMA baselines is 7 mJy and of ALMA baselines is 0.7 mJy.  \nIn the image domain, the median rms variation is 1.6 /uniF0A0 × /uniF0A0 10 8 K for the Top Set images of M87, and 1.3 /uniF0A0 × /uniF0A0 10 8 K for the normalized Top Set images. Because the peak image brightness is ∼ 10 10 K, these variations indicate that our images achieve a dynamic range of ∼ 50. Across the image, the minimum fractional uncertainty from imaging parameter choices is 8.7% across the Top Set images without /uniFB02 ux density normalization, and 7.9% with normalization. For comparison, the rms of the DIFMAP residual images is signi /uniFB01 cantly smaller: ∼ 3 /uniF0A0 × /uniF0A0 10 7 K across all days for the /uniFB01 ducial parameters. Thus, adding in the residual images would have negligible effect on our /uniFB01 nal reconstructions, especially relative to image differences between images from different imaging methods and from different parameter combinations using a single method.", '8. Image Validation': 'Having determined /uniFB01 ducial images of M87 from each imaging method on each observing day, we now perform additional validation tests to assess their reliability. First, in Section 8.1, we compare the residual telescope gains estimated for M87 with those for the calibrator source 3C 279 to  \nFigure 17. Uncertainties related to imaging assumptions from the eht-imaging parameter search on April /uniF0A0 11 data. The top row shows the mean image and associated uncertainties in the image domain, while the bottom row shows visibility amplitudes of the mean image and their uncertainties in the visibility domain. From left to right, panels correspond to the mean Top Set image, the standard deviation of the Top Set reconstructions, the standard deviation after normalizing each reconstructed image to the total /uniFB02 ux density of the mean image ( 0.62 Jy ) , and the standard deviation divided by the mean. In the image domain panels, contours are those of the mean image, at 10%, 20%, 40%, and 80% of the peak. In the visibility domain panels, contours are equally spaced logarithmically with three contours per decade; the thicker contours are 10 and 100 mJy ( left three ) and 10% ( right ) .  \n<!-- image -->  \ndetermine whether or not the signi /uniFB01 cant variation seen in the inferred LMT residual gains are consistent between the two sources. Next, in Section 8.2, we perform data resampling tests, using identical scripted reconstructions on data sets that exclude individual sites, to assess how sensitive the image morphology is to the loss of different combinations of data. Finally, in Section 8.3, we test whether or not the time variability seen in the /uniFB01 ducial images is related to differences in the sampled ( u , v ) coverage across the days.', '8.1. Gain Validation with the Calibrator 3C 279': 'The 2017 EHT observations of M87 were interleaved with those of the active galactic nucleus 3C 279. We have reconstructed images of 3C 279 using the three software packages used for M87 imaging in order to assess the consistency of these 3C 279 images with the M87 reconstructions in terms of the inferred time-variable residual gains. Because the sources are nearby on the sky ( separation 19 =  ) , the inferred time-variable residual gains at each site for the same day should be similar. We do not use the observations of 3C 279 to derive gain corrections for M87; instead, we compare the derived gains on both sources after independent imaging as a post hoc consistency test.  \nFigure 18 shows the aggregate baseline coverage for the interleaved EHT observations of 3C 279 in 2017 April. While the SPT could not observe M87, it participated in the observations of 3C 279, viewing it at an elevation of ∼ 6 ° ( i.e., a relative air mass of ∼ 10 ) . The addition of the SPT signi /uniFB01 cantly improves the north -south resolution of the array. Figure 19 shows the April 11 visibility amplitudes from 3C 279 after a priori and network calibration. From ALMA  \nFigure 18. Aggregate baseline coverage for EHT observations of 3C 279 in 2017 April. The coverage of M87 is shown in light gray for comparison.  \n<!-- image -->  \ninterferometric measurements, the total /uniFB02 ux density for 3C 279 ( 8 -10 Jy ) is nearly an order of magnitude higher than that of M87 /uniF0A0 ( Goddi et al. 2019 ) . As expected for a source with bright, linear jet features, the 3C 279 amplitudes vary strongly with baseline position angle and have a more complex structure on long baselines than M87 ( Figure 2 ) .  \nFigure 19. Visibility amplitudes of 3C 279 on April 11 as a function of projected baseline length after a priori and network calibration. The amplitudes are corrected for upward bias from thermal noise ( Equation ( 6 )) . Points are colored by baseline as in Figure 18.  \n<!-- image -->  \nWe imaged 3C 279 using both traditional CLEAN ( DIFMAP ) and RML ( eht-imaging and SMILI ) methods. Because of pronounced differences between the sources in different characteristics ( total /uniFB02 ux density, FOV, compact structure morphology ) , we do not use the same /uniFB01 ducial scripts derived for M87 on the 3C 279 data. Furthermore, as a consequence of having very few short baselines in the EHT array, we have found imaging 3C 279 to be more dif /uniFB01 cult than M87 because of its more extended structure. A full description of 3C 279 imaging procedures and results will be presented separately.  \nFigure 20 shows reconstructed images of 3C 279 from all three methods using data from April 11. The source exhibits two compact and bright features with a separation of ∼ 100 /uni03BC as: a primary component extended in the northwest to southeast direction, and a secondary component perpendicular to the /uniFB01 rst. This secondary component extends from the core in the direction of the 3C 279 jet observed at lower frequencies ( see, e.g., Lister et al. 2016; Jorstad et al. 2017 ) .  \nIn Figure 21, we compare the interleaved multiplicative station gains for M87 and 3C 279 on April 5. The station gains were derived via self-calibration to the /uniFB01 ducial images of M87 from the three imaging pipelines ( Section 7 ) and on the three images of 3C 279 ( Figure 20 ) . In Table 6, we present the median gains for the two sources and compare them to the expected a priori visibility amplitude error budget, assuming nominal pointing and focus ( Paper III ) .  \nThe gains derived from images produced with the three imaging methods are consistent. For all stations except the LMT, the derived gains are time-variable, but close to unity. Issues with the performance of the LMT ( Section 4,andPaperIII ) result in a larger median gain correction at that station, in line with the gain corrections derived directly from the visibility amplitudes via the use of crossing tracks on M87 described in Appendix B. Furthermore, the interleaved LMT gain curve has large variations with time on both sources. In particular, the M87 gains at the LMT have large, single-scan excursions not seen in the 3C 279 gain trends. These excursions are likely due to poor pointing on M87, which is nearly an order of magnitude fainter than 3C 279.  \nIn Table 14 of Appendix F, we present the median gains for the two sources, compared to the expected a priori visibility amplitude error budget, across all days. Figure 34 in  \nAppendix F compares the gain variations recovered from imaging M87 and 3C 279 across all days for the SMT and LMT stations. In some cases, the absolute gain is not identical across the different sources and imaging pipelines. This variation can be partly ascribed to the large uncertainty in the total /uniFB02 ux density of 3C 279. However, in all cases the relative gain trends with time are broadly consistent, except for the occasional large excursions seen in the LMT M87 gains.  \nThe consistency of the derived gain variations between the two sources and the different imaging methods for all stations on all days, particularly for the large corrections at LMT, provides con /uniFB01 dence that these corrections are not the result of imaging artifacts or missing structure in M87 reconstructions. In Paper VI we show that derived gains from /uniFB01 tting parametric models are similar to those derived from imaging.', '8.2. Image Dependence on Sites': "Our comparisons of reconstructed images on different days show that the primary image features in M87 are not sensitive to small changes in baseline coverage ( Figure 11 ) . We now explore how more severe changes in baseline coverage affect imaging. In particular, we assess whether or not the image structure is closely tied to the data from any one site.  \nThe top row of Figure 22 shows reconstructed images of M87 using the /uniFB01 ducial eht-imaging parameters after excluding all data on baselines to one of the /uniFB01 ve EHT sites that observed the source. As shown in Figure 1, removing some sites entirely from the array does not lead to a signi /uniFB01 cant decrease in the length of the longest baseline in either the north -south or east -west direction. The notable exceptions are Chile ( ALMA and APEX ) and Pico Veleta ( PV ) in Spain. The Chilean sites are part of all baselines with north -south lengths v 2.5 G l > /uni2223 /uni2223 , while PV provides the only baselines in the northeast direction with lengths longer than 5 G λ . Thus, the loss of either of these sites signi /uniFB01 cantly elongates the synthesized beam and impedes the reconstruction of the ring feature.  \nThe presence of a ring in our reconstructed images of M87 is resilient to several site losses, even though the /uniFB01 ducial imaging parameters were not chosen using a reduced ( u , v ) coverage. The ring is successfully reconstructed when we entirely omit both Hawai ' i stations ( SMA and JCMT ) , even though these stations are required to sample the two visibility nulls ( at ∼ 3.4 and ∼ 8.3 G λ ; see Section 2.1 ) . The /uniFB01 ducial scripts also reconstruct a ring when LMT data are entirely omitted, but the ring orientation is reversed. Unsurprisingly, the site losses that result in the poorest reconstructions of the ring are those that most degrade the synthesized beam, namely Chile and PV.  \nWe also tested whether or not severe calibration errors could have affected imaging. To do so, we applied our /uniFB01 ducial scripts after omitting calibrated visibility amplitudes to a speci /uniFB01 ed site while still retaining closure phases and closure amplitudes. The lower panels of Figure 22 show this test for M87 reconstructions using the eht-imaging pipeline ( for similar results with SMILI , see Appendix D ) . 117 Comparing the top and bottom row of Figure 22, it is apparent that even sites with poor amplitude calibration ( e.g., the LMT ) can provide essential constraints through closure quantities. For this reason, we do not /uniFB02 ag data with poor amplitude calibration before imaging. We also show an image that uses only closure quantities for the reconstruction, therefore using no calibration information for  \nFigure 20. Representative images of 3C 279 from the April /uniF0A0 11 EHT observations produced using DIFMAP , eht-imaging , and SMILI . To simplify visual comparisons and display the images at similar resolutions, the images are restored with circular Gaussian beams of 20, 17.1, and 18.6 /uni03BC as FWHM, respectively.  \n<!-- image -->  \nTable 6 5  \nMedian, 25th, and 75th Percentile Residual Gain Corrections for M87 and 3C 279 on April /uniF0A0  \nNote. These gains were derived via self-calibration ( with no systematic error included ) . The error budget on a /uniF0A0 priori calibration is derived in Paper III. Note that the median gain corrections for ALMA, APEX, SMA, and JCMT can reasonably be much smaller than this error budget because network calibration has already been applied. The variation in the recovered gains among pipelines is partly due to the large uncertainty in the total /uniFB02 ux density ( between 8 and 10 Jy ) and total compact /uniFB02 ux density of 3C 279.  \nFigure 21. Multiplicative residual station gains for the SMT ( left ) and LMT ( right ) derived from the 3C 279 images ( Figure 20 ) and /uniFB01 ducial M87 images ( Figure 14 ) from the three imaging pipelines on April 5. Gains for the /uniFB01 ducial images of M87 are shown in red; those for 3C 279 are shown in blue. The particularly large excursions on the LMT M87 gains are likely due to poor pointing. Note that LMT could not observe 3C 279 before 2 h 30 UTC.  \n<!-- image -->  \nany site. This image demonstrates that closure quantities are suf /uniFB01 cient to recover the diameter and asymmetry of the ring in the /uniFB01 ducial images.", '8.3. Tests of the Observed Time Variability': "Our /uniFB01 ducial images of M87 ( Figure 11 ) over the four observing days show an evolution in the ring structure between  \nthe two observations on April 5, 6 and the two observations on April 10, 11 ( Figure 33 ) . However, the data used to reconstruct these images do not identically sample the ( u , v ) plane ( Figure 1 ) . To verify whether or not the observed evolution in the /uniFB01 ducial M87 images is intrinsic, it is necessary to con /uniFB01 rm that the evolution is not an artifact imprinted on the reconstructions by the changing ( u , v ) coverage sampled by the EHT on the different observing days.  \nTo this end, we constructed reduced data sets from the April 6 and 11 data with ' overlapping ' baseline coverage. Speci /uniFB01 -cally, we /uniFB02 ag all ( u , v ) points along a particular baseline track that are not common to the April 6 and 11 observations ( i.e., we /uniFB02 ag points on April 6 with no measurement on April 11 within 0.01 G λ , and vice versa ) . The number of ( 10 s averaged ) visibilities is reduced from ∼ 18,000 to ∼ 5,000 in the overlapping data set. We then used the eht-imaging script with the identi /uniFB01 ed /uniFB01 ducial parameters ( Section 6.2.2 ) to reconstruct images on both days with these subsampled data sets.  \nFigure 23 compares the reconstructions made from these overlapping ( u , v ) data to our original images from the full data sets. The results show no signi /uniFB01 cant variation in the morphology due to the changing ( u , v ) coverage. In particular, the enhanced southwest brightness in the April /uniF0A0 11 image  \nFigure 22. Example reconstructions of M87 on April /uniF0A0 11 after omitting visibilities to each geographical site. The top row shows reconstructions that exclude all baselines to the indicated site ( i.e., mimicking an observation without that site ) ; the bottom row shows reconstructions that exclude all visibility amplitudes on baselines to the indicated site, while still including closure phases and closure amplitudes ( i.e., exploring the limit of arbitrarily poor amplitude calibration ) . The lowerleft panel shows an image reconstructed after excluding visibility amplitudes to all sites ( i.e., using only closure quantities ) . These images were reconstructed using the eht-imaging pipeline with its /uniFB01 ducial parameters ( see Section 6.3.3 ) . The ellipse in each panel shows the corresponding synthesized beam with uniform weighting, but the image is not convolved with this beam. While the /uniFB01 ducial parameters used are not necessarily optimal for imaging these modi /uniFB01 ed data sets, the reconstructed images of M87 are rather resilient to the loss of most sites, especially if the closure information is preserved, demonstrating that our results are insensitive to calibration assumptions. Figure 32 in Appendix D shows corresponding results for the DIFMAP and SMILI pipelines.  \n<!-- image -->  \nFigure 23. Validating the observed variation in the reconstructed images between April 6 and 11. Top-left panels: /uniFB01 ducial eht-imaging image reconstructions from observations on April 6 and 11 ( Section 7.1 ) . Top-right panels: reconstructions using the same /uniFB01 ducial script and parameters, but using data only on the ( u , v ) points which overlap on April 6 and 11. The reconstructions using only the overlapping ( u , v ) coverage show similar time evolution ( a shift in brightness to the southwest; see Appendix E ) to the reconstructions from the original data sets, indicating that the observed variation in images is intrinsic. Bottom panels: plots of two selected closure phases ( LMT -SMA -SMT and ALMA -LMT -SMA ) that clearly identify intrinsic structural evolution between April 6 and 11. The colored points identify the closure phases of the scan-averaged data, and the phases shown as solid and dotted colored lines are taken from the reconstructions of full and overlapping coverage data, respectively.  \n<!-- image -->  \nrelative to the April /uniF0A0 6 image is still prominent in the overlapping coverage reconstructions. This test shows that the reconstructed evolution is not due, at least solely, to the changes in ( u , v ) coverage between days, although each image is still affected by the limited ( u , v ) coverage. The bottom panels of Figure 23 show selected closure phases on April /uniF0A0 6 and 11 that are signi /uniFB01 cantly different; reconstructed images with both the full and overlapping data sets recover this signi /uniFB01 cant evolution of closure phases. We discuss the observed time variability further in Appendix E.", '9. Image Analysis': 'To assess the consistency of our images with each other and determine which speci /uniFB01 c image features are most reliable, we de /uniFB01 ne and measure several properties that characterize asymmetric rings: the diameter d , width w , orientation angle η , azimuthal brightness asymmetry A , and fractional central brightness f C. This approach is also motivated by GRMHD simulations of M87, which predict an image dominated by a bright and nearly circular ring morphology across a large range of parameters ( Paper V ) . In Section 9.1, we de /uniFB01 ne these quantities, explain how they are measured, and give procedures to estimate their uncertainties. In Section 9.2, we evaluate biases and systematic errors in these measurements using tests with synthetic data. In Section 9.3, we apply our results to estimate properties of the /uniFB01 ducial images of M87 presented in Section 7. In Paper VI, we explore model /uniFB01 tting in more detail, both through visibility-domain /uniFB01 ts of geometrical models and through image-domain analysis of the reconstructions from each imaging pipeline.', '9.1. Ring Parameter De /uniFB01 nitions': 'To identify a ring in a given reconstructed image, we search for the central point from which radial pro /uniFB01 les are peaked at a similar distance. From a candidate ring center position ( x , y ) , we sample the linearly interpolated image equally in azimuthal angle θ between 0 ° and 360 ° and in radius r between 0 and 50 /uni03BC as to obtain a transformed image: I ( r , θ ; x , y ) . 118 Then, for each radial pro /uniFB01 le at /uniFB01 xed θ , we identify the distance r pk at which the angular pro /uniFB01 le assumes its peak brightness. We de /uniFB01 ne the ring radius at a given position, r xy , pk ¯ ( ) as the mean of these peak distances:  \nr x y I r x y r x y r x y ; , argmax , ; , , , ; , . 16 r pk pk pk 0,2 q q q = =á ñ q p Î ( ) [ ( )] ¯ ( ) ( ) ( ) [ ]  \nTo estimate an associated uncertainty, we use the standard deviation xy , r s ( ) ¯ of the r pk ( θ ; x , y ) values. For a similar  \nmethod using Gaussian /uniFB01 ts to the radial pro /uniFB01 les, see Johannsen et al. ( 2016 ) .  \nTo /uniFB01 nd the ring center ( x 0, y 0 ) , we search over ( x , y ) and identify the position that minimizes the normalized radial peak dispersion:  \nx y x y r x y , argmin , , . 17 r xy 0 0 pk , s = /uni23A1 /uni23A3 /uni23A2 /uni23A4 /uni23A6 /uni23A5 ( ) ( ) ¯ ( ) ( ) ¯ ( )  \nThe measured diameter d is then twice r pk ¯ measured from the identi /uniFB01 ed center ( x 0, y 0 ) ,  \nd r x y 2 , , 18 pk 0 0 = ¯ ( ) ( )  \nand the associated uncertainty σ d in the diameter is  \nx y 2 , . 19 d r 0 0 s s = ( ) ( ) ¯  \nAlthough we have designed our analysis tools to speci /uniFB01 cally search for circular features, σ d / d can be interpreted as a measure of the circularity of the identi /uniFB01 ed feature in an individual image. See Paper VI for further discussion of the circularity of these images and model /uniFB01 ts to the visibility data.  \nTo avoid spurious detections when searching for the center location ( Equation ( 17 )) ,we /uniFB01 rst blur the image with a 2 /uni03BC as FWHM Gaussian ( approximately the pixel size of the original reconstructions ) and threshold the search image below 5% of the peak brightness. We also restrict the range of allowed diameters to 10 -100 /uni03BC as. We use the original ( unconvolved and unthresholded ) image for all subsequent analysis. Because the remaining features are all estimated using the /uniFB01 xed center coordinates estimated by this procedure, we suppress the center coordinates for the remainder of our discussion, and we denote the unwrapped, centered image as I ( r , θ ) .  \nAfter identifying the ring center and diameter, we measure the remaining characteristic features. The ring width w is determined by measuring the FWHM of each radial slice at constant θ and taking the mean. To avoid bias in the measurement from the resampled image I ( r , θ ) having a nonzero mean /uniFB02 oor value outside the ring, we subtract the value I Ir 50 as, floor max m q =á = ñ q ( ) from each radial pro /uniFB01 le before computing the FWHM:  \nw I r I FWHM , . 20 floor q =á -ñ [ ( ) ] ( )  \nWe estimate an uncertainty σ w using the standard deviation of the set of FWHMs. Note that the measured width is dependent upon both the intrinsic width of the source and the /uniFB01 nite resolution of the array. Thus, w should be viewed only as an upper limit on the intrinsic ring width, and it is upward biased by the application of a restoring beam ( e.g., for DIFMAP reconstructions ) . We explore this bias further in Appendix G.  \nWe de /uniFB01 ne the ring orientation angle η ( measured east of north ) by considering the angular pro /uniFB01 les I ( r , θ ) at /uniFB01 xed r .At each radius, we /uniFB01 nd the argument of the /uniFB01 rst angular mode ( m = 1 ) of the angular pro /uniFB01 le and then de /uniFB01 ne the overall orientation angle η as the circular mean of these angles over the ring width; from r d w 2 in = -( ) to r out /uniF0A0 = /uniF0A0 ( d + w ) / 2. That is,  \nI e d Arg . 21 i r r r 0 2 , in out ò h q q = p q Î /uni23A1 /uni23A3 /uni23A2 /uni23A4 /uni23A6 /uni23A5 ( ) ( ) [ ]  \nSimilarly, we de /uniFB01 ne the associated uncertainty ση as the circular standard deviation of the angle measurements for r r r , in out Î [ ] .  \nTo characterize the degree of azimuthal asymmetry in a ring, we compute the normalized amplitude of the /uniFB01 rst angular mode for radii between r in and r out and again take the mean. That is,  \nA I e d I d . 22 i r r r 0 2 0 2 , in out ò ò q q qq = p q p Î ( ) ( ) ( ) [ ]  \nThe associated uncertainty σ A is the standard deviation of the asymmetry at each r r r , in out Î [ ] . The asymmetry A takes values in the range from 0 to 1, with 0 corresponding to perfect azimuthal symmetry and 1 corresponding to a delta function concentrating all of the /uniFB02 ux density at a single orientation angle. The crescent models de /uniFB01 nedinSection6.1 have A /uniF0A0 = /uniF0A0 0.23, with angular brightness pro /uniFB01 les I A 1 2 cos q q h µ + -( ) ( ) .  \nFinally, we de /uniFB01 ne a fractional central brightness ( or inverse contrast ratio ) f C as the ratio of the mean brightness interior to the ring to the mean brightness around the ring. To de /uniFB01 ne the interior brightness, we average over a disk of radius 5 /uni03BC as. That is,  \nf Ir Id , 2, . 23 r C , 0,5 as 0,2 q q = á ñ á ñ q m q p Î Î ( ) ( ) ( ) [ ] [ ]  \nWe /uniFB01 nd that this statistic has an extremely large scatter across the Top Sets, primarily because the interior brightness can become arbitrarily low. Thus, our imaging methods only securely identify an upper limit on f C.', '9.2. Tests with Synthetic Data Reconstructions': 'We tested the analysis methods described in Section 9.1 on synthetic data from a crescent model with asymmetry parameter A /uniF0A0 = /uniF0A0 0.23 oriented at a orientation angle η /uniF0A0 = /uniF0A0 150 ° ( Appendix C ) , and on the GRMHD simulation image shown in Figure 10. For both sources, we measured the ring features de /uniFB01 ned in Section 9.1 for every image reconstructed with parameter combinations in the M87 Top Sets ( see Section 6 ) . While some of these images may not meet the χ 2 /uniF0A0 < /uniF0A0 2 constraint imposed for M87 images, using the M87 Top Set allows us to assess the expected variation resulting from imaging choices irrespective of image-data consistency.  \nThe results of this analysis for the diameter ( d ) , width ( w ) , and orientation angle ( η ) are displayed in Figure 24. The points shown correspond to the measured quantities from the /uniFB01 ducial reconstructions. These measurements have two distinct sources of error: the intrinsic measurement uncertainty on each quantity from a single image, and the uncertainty in the quantity from varying imaging parameters across the Top Set. We combine these two sources of error in the error bars in Figure 24 by taking the quadrature sum of the measurement uncertainty from the /uniFB01 ducial image and the median absolute deviation of the parameter estimates across the Top Set.  \nIn both simulated data tests, the ring diameter d is the most accurately recovered quantity; however, the diameter measurement is correlated with the ring width and is biased downward by several microarcseconds when the image is blurred ( see Appendix G ) . In the simulated crescent, the diameter of the unblurred model is 44 /uni03BC as; the value measured using our approach on the ground truth image is pushed down to 43 /uni03BC as because of the Gaussian convolution ( Appendix G ) . Taking the median across all days, the diameters extracted from the  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \nFigure 24. Measurements of ring features on /uniFB01 ducial reconstructions of a crescent model ( top row ) and a GRMHD simulation snapshot ( bottom row ) from images reconstructed with the /uniFB01 ducial parameters for three imaging pipelines. From left to right, panels show the measured ring diameter d , width w , and orientation angle η . The DIFMAP results are shown for images restored with both a 10 /uni03BC as ( cyan ) and a 20 /uni03BC as ( blue ) Gaussian beam. The eht-imaging ( red ) and SMILI ( green ) are shown for the unblurred images. Solid lines indicate the ground truth values of the three quantities in the crescent model, and the photon ring diameter of the GRMHD simulation. Dashed lines indicate the measured values from the ground truth images, and dotted lines indicate the measured values from the ground truth images convolved with a 20 /uni03BC as FWHM circular Gaussian beam. The error bars are computed as the quadrature sum of the measurement uncertainty from the /uniFB01 ducial images ( Section 9.1 ) and the median absolute deviation of the parameter estimates across the Top Set.  \n<!-- image -->  \n/uniFB01 ducial DIFMAP crescent model reconstructions have a median value of 40 /uniF0A0 ± /uniF0A0 2 /uni03BC as when restored with a 20 /uni03BC as FWHM beam. When they are instead restored with a smaller 10 /uni03BC as beam, the DIFMAP results become more accurate ( compared to the model values ) , with a measured diameter of 43 /uniF0A0 ± /uniF0A0 2 /uni03BC as. The RML crescent model reconstructions measure a median diameter of 42 /uniF0A0 ± /uniF0A0 1 /uni03BC as.  \nTheGRMHDimagethatweuse ( Paper V ) has a lensed photon ring with a diameter of 9.79 M / D /uniF0A0 ≈ /uniF0A0 35.5 /uni03BC as; 119 we recover a value of 36.5 /uni03BC as from the ground truth image, indicating that even with a perfect image reconstruction, /uniFB01 ne-scale substructure and extended /uniFB02 ux in an image can bias the ring diameter measurement away from the photon ring value. The ring diameter measured from the eht-imaging and SMILI GRMHD /uniFB01 ducial reconstructions ( Figure 24 ) has a median value of 38 /uniF0A0 ± /uniF0A0 2 /uni03BC as across all four observed days. However, for this image, the DIFMAP results are strongly dependent upon the chosen restoring beam. The DIFMAP reconstructions restored with a 20 /uni03BC as beam lack prominent rings, leading to poor recovery of the ring diameter ( 26 ± 20 /uni03BC as ) . When blurred with a smaller 10 /uni03BC as beam, the DIFMAP results align with the RML methods but have larger uncertainties due to larger scatter across the Top Set ( 37 ± 11 /uni03BC as ) .  \nAs shown in Figure 24, the measured orientation angles from the crescent model reconstructions using the RML imaging pipelines have a median value of 141 ° /uniF0A0 ± /uniF0A0 5 ° , which is moderately discrepant from the true value of 150 ° . The values from the DIFMAP images are even more discrepant and are unstable among days. The GRMHD image has no de /uniFB01 ned a priori orientation angle; the measured orientation angles from  \nall pipelines are more stable than for the simple crescent model, and are consistent with the value measured from the blurred simulation image ( 105 ° ) . Both of these cases indicate that our procedure may underestimate uncertainties on η .  \nBoth the crescent and GRMHD models have intrinsic widths that are much narrower than the resolution of the EHT. As expected from the 20 /uni03BC as restoring beam, the DIFMAP reconstructions lead to a larger measurement of ring width than the RML reconstructions. When a 10 /uni03BC as FWHM beam is used to restore the DIFMAP images, the measured width aligns with the RML results. Nonetheless, all reconstructions have widths that are systematically biased upward from the true values, and the extracted ring widths from image-domain /uniFB01 tting can at best be viewed as upper limits.  \nSimilarly, the brightness depression contrast ratio f C is highly sensitive to the image resolution and particular imaging choices. In particular, we /uniFB01 nd that the Top Sets from the RML imaging pipelines have an extremely large scatter in f C. When blurred to the same 20 /uni03BC as resolution, both RML and CLEAN methods give measurements of f C that are consistently in the range 0.2 -0.5 ( Appendix G ) . Consequently, we can only determine an upper limit for f C.  \nFor all the imaging pipelines, we /uniFB01 nd that the scatter in ring diameters across the Top Set ( typically /uniF088 1 /uni03BC as ) is subdominant to the intrinsic measurement uncertainty estimated from a single image ( typically 1 -2 /uni03BC as ) . Thus, choices made in the imaging process do not signi /uniFB01 cantly affect the measured ring diameter. In contrast, the other measured features have error budgets that are more evenly divided between intrinsic uncertainty in a single image and the scatter across the Top Set.', '9.3. Results for M87': "Table 7 lists the values of all ring parameters measured for the /uniFB01 ducial M87 reconstructions for each day from the three  \nTable 7  \nDiameter d , Width w , Orientation Angle η , Asymmetry A , and Floor-to-ring Contrast Ratio f C Measured from the Fiducial M87 Images for Each Day from Each Method ( Section 7.1 )  \nNote. The DIFMAP results are presented for a restoring beam of 20 /uni03BC as. The estimated diameters and orientation angles are generally consistent across different imaging pipelines on all four days, but the measured width, contrast, and asymmetry differ between imaging pipelines.  \nparameter surveys. As in the previous section, for the /uniFB01 ducial images the uncertainties are computed by adding the scatter in the measured quantities across the Top Sets in quadrature to the intrinsic measurement uncertainty.  \nFigure 25 plots the measured diameter, width, and orientation angle from each method over all four days, and Figure 26 visualizes these parameters. Across all days, the DIFMAP reconstructions recover a median diameter of 40 /uniF0A0 ± /uniF0A0 3 /uni03BC as when restored with a 20 /uni03BC as beam, and 44 /uniF0A0 ± /uniF0A0 5 /uni03BC as when restored with a 10 /uni03BC as beam. The RML methods recover a median diameter of 41 /uniF0A0 ± /uniF0A0 2 /uni03BC as. For each imaging pipeline, there is slight upward trend in the diameter over time, increasing by ≈ 2 /uni03BC as from the /uniFB01 rst to the last day of the observing campaign. However, this trend is well within the estimated uncertainty of the diameter measurements. The measured ring diameters are consistent among all imaging pipelines and largely unchanged from the /uniFB01 rst April 11 images produced from early-release engineering data ( Section 5 ) , for which we measure a median diameter of 41 /uniF0A0 ± /uniF0A0 3 /uni03BC as. The other parameters are less consistent from these /uniFB01 rst images to the /uniFB01 nal /uniFB01 ducial images selected by the parameter survey, indicating more sensitivity in these parameters to the data quality and the imaging method.  \nAs in the synthetic data results presented in Figure 24,the beam-convolved DIFMAP reconstructions produce larger measured widths than the RML methods. The width measurements become consistent when the DIFMAP images are restored with a smaller 10 /uni03BC as beam, but the ring width remains limited by the resolution of the EHT. Thus, from these imaging results, we can only /uniFB01 rmly conclude that the ring width is /uni0084 20 /uni03BC as.  \nIn the reconstructions from all three imaging methods, there is a counterclockwise trend in the orientation angle from April 5 to 11, which is consistent with the apparent shift in brightness along the ring remarked on in Section 7.1. However, we emphasize that this ≈ 20 ° change could be the result of spurious azimuthal structure, and we have seen that our orientation angle uncertainties may be underestimated ( see Figure 24 ) . Even if this shift in angle is physical, it does not necessarily indicate continuous motion or a /uniFB02 ow direction associated with the black  \nhole in M87, and it is opposite to the inferred rotation direction for the large-scale jet ( Walker et al. 2018 ) .  \nIdentifying a ring in the /uniFB01 ducial images allows us to ' unwrap ' the brightness distribution around the ring and display the image in I ( r , θ ) space. Figure 27 shows these unwrapped ring pro /uniFB01 les of the /uniFB01 ducial images for all pipelines and days. In Appendix I,we present more diagnostic /uniFB01 gures to characterize the radial and azimuthal structure of the recovered rings.  \nEven under the assumption that the ring in our images corresponds to the lensed photon ring of an SMBH ( i.e., d ≈ ( 5.0 ± 0.2 ) R s / D ; Bardeen 1973; Johannsen & Psaltis 2010 ) , translating the estimated ring diameter to a Schwarzschild radius -or, equivalently, to a black hole mass via R s /uniF0A0 = /uniF0A0 GM BH / c 2 -requires evaluating biases and uncertainties related to the unresolved image structure, the black hole vector spin, and the estimated distance to M87. We perform this analysis separately, in Paper VI.", '10. Summary and Conclusions': "We have presented imaging analysis and results for EHT observations of M87 in 2017 April. These observations achieved an angular resolution ( synthesized beam size ) of ∼ 20 /uni03BC as, and they are the /uniFB01 rst EHT observations with suf /uniFB01 cient sensitivity and baseline coverage to reconstruct images. Images can provide agnostic results, potentially identifying entirely unexpected features. However, a multitude of imaging techniques and parameters can lead to results that depend on speci /uniFB01 c imaging choices. Consequently, we adopted a staged imaging approach to control and evaluate potential biases.  \nWe designed our /uniFB01 rst imaging stage to minimize shared bias among different groups of imaging experts and to test whether or not their independent inferences were consistent. In this stage, described in Section 5, we separately imaged M87 in four imaging teams. These imaging teams did not communicate about their methods or /uniFB01 ndings while analyzing the M87 data. Each team relied upon the judgment of its members to select strategies for data /uniFB02 agging, calibration, and imaging, using whatever software  \nFigure 25. Measured ring properties on the /uniFB01 ducial images of M87 produced with all three imaging pipelines. From left to right, panels show the measured ring diameter d , width w , and orientation angle η . The DIFMAP results are shown for images restored with both a 10 /uni03BC as ( cyan ) and a 20 /uni03BC as ( blue ) Gaussian beam. The eht-imaging ( red ) and SMILI ( green ) are shown for the unblurred images. The three imaging pipelines produce consistent measurements of the ring diameter across all days. The measured orientation angles indicate a modest shift between April 5, 6 and 10, 11. The error bars are computed in the same way as in Figure 24.  \n<!-- image -->  \nFigure 26. Summary of the estimated ring properties overlaid on the April 11 /uniFB01 ducial images from each imaging pipeline. Our procedure estimates the ring diameter d , width w , orientation angle η , and fractional central brightness f C, as well as the asymmetry A ( not shown ) . In each panel, the magenta cross indicates the location of peak ring brightness, and the dashed green circle shows the region used to de /uniFB01 ne interior brightness for f C.  \n<!-- image -->  \nFigure 27. Unwrapped ring pro /uniFB01 les of the /uniFB01 ducial images from April /uniF0A0 5to11 ( top to bottom ) and for the three imaging pipelines ( left to right ) . The columns are each scaled as in Figure 26. The estimated radius d / 2 is shown with a horizontal line, with dotted lines denoting the associated uncertainty ( Section 9.1 ) . Horizontal dashed lines at ( d ± w ) / 2 show the measured ring width. Vertical blue lines give the orientation angle η and its uncertainty. The magenta cross marks the peak brightness in each reconstruction.  \n<!-- image -->  \nand techniques that they judged to be useful. The /uniFB01 nal images from the teams were blindly submitted for comparison. These /uniFB01 rst images of M87 all showed a prominent ring of ∼ 40 /uni03BC as diameter with enhanced brightness in the south. The four teams made different assumptions about the total compact /uniFB02 ux density of M87, varying from ∼ 0.4 to 0.9 Jy. The images also showed differences in the azimuthal ring structure.  \nWe designed our second imaging stage to select imaging parameters objectively and to assess how parameter choices correspond to uncertainties in the /uniFB01 nal images. In this stage, described in Section 6, we developed scripted imaging pipelines for one established ( DIFMAP ) and two EHT-tailored imaging libraries ( eht-imaging and SMILI ) . For each method, we conducted a parameter survey, exploring ∼ 10 3 -10 4 parameter combinations for each imaging pipeline. We then reconstructed images from synthetic data sets and compared them to the corresponding ground truth images. To determine which parameter combinations in the survey produced acceptable images, we imposed an image-data consistency requirement for reconstructions of M87, as well as an image /uniFB01 delity requirement for reconstructions of synthetic data.  \nThis procedure identi /uniFB01 ed a ' Top Set ' of parameter choices as well as a single ' /uniFB01 ducial ' combination of parameters for each of the three surveys. When applied to M87 data, the Top Set reconstructions provide a measure of uncertainty corresponding to choices made in the imaging procedure. In addition, the /uniFB01 ducial image ( Figure 11 ) acts as a single representative image from each imaging pipeline that has been identi /uniFB01 ed based on objective performance, rather than expert judgment.  \nOur reconstructed images of M87 are consistently dominated by a ring of roughly 40 /uni03BC as diameter, with enhanced brightness in the south. The ring is present in both imaging stages, across all imaging pipelines, and all observing days. The residual timedependent station gains derived from our images are consistent between the interleaved EHT observations of M87 and the bright quasar 3C 279. Our reconstructed images of M87 show variation between the observations on April 5 and 11, as expected from the variations seen in the underlying data ( Paper III ) .However, we cannot unambiguously associate the variability with evolution of a speci /uniFB01 c image feature.  \nThese /uniFB01 rst images from the EHT achieve the highest angular resolution in the history of ground-based VLBI. The images reveal the central engine of M87 and are dominated by a ∼ 40 /uni03BC as diameter ring with pronounced azimuthal asymmetry. Such an asymmetry is expected for relativistic rotation of emitting material near a black hole. For an assumed distance of 16.8 Mpc, the ring diameter is consistent with the expected ' shadow ' of a ∼ 6.5 /uniF0A0 × /uniF0A0 10 9 M e black hole.  \nThe authors of this Letter thank the following organizations and programs: the Academy of Finland ( projects 274477, 284495, 312496 ) ; the Advanced European Network of E-infrastructures for Astronomy with the SKA ( AENEAS ) project, supported by the European Commission Framework Programme Horizon 2020 Research and Innovation action under grant agreement 731016; the Alexander von Humboldt Stiftung; the Black Hole Initiative at Harvard University, through a grant ( 60477 ) from the John Templeton Foundation; the China Scholarship Council; Comisión Nacional de Investigación Cientí /uniFB01 ca y Tecnológica ( CONICYT, Chile, via PIA ACT172033, Fondecyt 1171506, BASAL AFB-170002, ALMA-conicyt 31140007 ) ; Consejo Nacional de Ciencia y Tecnología ( CONACYT, Mexico,  \nprojects 104497, 275201, 279006, 281692 ) ; the Delaney Family via the Delaney Family John A. Wheeler Chair at Perimeter Institute; Dirección General de Asuntos del Personal AcadémicoUniversidad Nacional Autónoma de México ( DGAPA-UNAM, project IN112417 ) ; the European Research Council Synergy Grant ' BlackHoleCam: Imaging the Event Horizon of Black Holes ' ( grant 610058 ) ; the Generalitat Valenciana postdoctoral grant APOSTD / 2018 / 177; the Gordon and Betty Moore Foundation ( grants GBMF-3561, GBMF-5278 ) ; the Istituto Nazionale di Fisica Nucleare ( INFN ) sezione di Napoli, iniziative speci /uniFB01 che TEONGRAV; the International Max Planck Research School for Astronomy and Astrophysics at the Universities of Bonn and Cologne; the Jansky Fellowship program of the National Radio Astronomy Observatory ( NRAO ) ; the Japanese Government ( Monbukagakusho: MEXT ) Scholarship; the Japan Society for the Promotion of Science ( JSPS ) Grant-in-Aid for JSPS Research Fellowship ( JP17J08829 ) ; JSPS Overseas Research Fellowships; the Key Research Program of Frontier Sciences, Chinese Academy of Sciences ( CAS, grants QYZDJ-SSW-SLH057, QYZDJ-SSWSYS008 ) ; the Leverhulme Trust Early Career Research Fellowship; the Max-Planck-Gesellschaft ( MPG ) ; the Max Planck Partner Group of the MPG and the CAS; the MEXT / JSPS KAKENHI ( grants 18KK0090, JP18K13594, JP18K03656, JP18H03721, 18K03709, 18H01245, 25120007 ) ; the MIT International Science and Technology Initiatives ( MISTI ) Funds; the Ministry of Science and Technology ( MOST ) of Taiwan ( 105-2112-M-001-025-MY3, 106-2112-M-001-011, 106-2119-M-001-027, 107-2119-M-001017, 107-2119-M-001-020, and 107-2119-M-110-005 ) ; the National Aeronautics and Space Administration ( NASA, Fermi Guest Investigator grant 80NSSC17K0649 ) ; the National Institute of Natural Sciences ( NINS ) of Japan; the National Key Research and Development Program of China ( grant 2016YFA0400704, 2016YFA0400702 ) ; the National Science Foundation ( NSF, grants AST-0096454, AST-0352953, AST-0521233, AST-0705062, AST-0905844, AST-0922984, AST-1126433, AST-1140030, DGE-1144085, AST-1207704, AST-1207730, AST-1207752, MRI-1228509, OPP-1248097, AST-1310896, AST-1312651, AST-1337663, AST-1440254, AST-1555365, AST-1715061, AST-1614868, AST-1615796, AST-1716327, OISE-1743747, AST-1816420 ) ; the Natural Science Foundation of China ( grants 11573051, 11633006, 11650110427, 10625314, 11721303, 11725312, 11873028, 11873073, U1531245, 11473010 ) ; the Natural Sciences and Engineering Research Council of Canada ( NSERC, including a Discovery Grant and the NSERC Alexander Graham Bell Canada Graduate Scholarships-Doctoral Program ) ; the National Youth Thousand Talents Program of China; the National Research Foundation of Korea ( grant 2015R1D1A1A01056807, the Global PhD Fellowship Grant: NRF2015H1A2A1033752, and the Korea Research Fellowship Program: NRF-2015H1D3A1066561 ) ; the Netherlands Organization for Scienti /uniFB01 c Research ( NWO ) VICI award ( grant 639.043.513 ) and Spinoza Prize SPI 78-409; the New Scienti /uniFB01 c Frontiers with Precision Radio Interferometry Fellowship awarded by the South African Radio Astronomy Observatory ( SARAO ) , which is a facility of the National Research Foundation ( NRF ) ,an agency of the Department of Science and Technology ( DST ) of South Africa; the Onsala Space Observatory ( OSO ) national infrastructure, for the provisioning of its facilities / observational support ( OSO receives funding through the Swedish Research Council under grant 2017-00648 ) the Perimeter Institute for Theoretical Physics ( research at Perimeter Institute is supported by the Government of Canada through the Department of  \nInnovation, Science and Economic Development Canada and by the Province of Ontario through the Ministry of Economic Development, Job Creation and Trade ) ; the Russian Science Foundation ( grant 17-12-01029 ) ; the Spanish Ministerio de Economía y Competitividad ( grants AYA2015-63939-C2-1-P, AYA2016-80889-P ) ; the State Agency for Research of the Spanish MCIU through the ' Center of Excellence Severo Ochoa ' award for the Instituto de Astrofísica de Andalucía ( SEV-2017-0709 ) ;the Toray Science Foundation; the US Department of Energy ( USDOE ) through the Los Alamos National Laboratory ( operated by Triad National Security, LLC, for the National Nuclear Security Administration of the USDOE ( Contract 89233218CNA000001 )) ; the Italian Ministero dell ' Istruzione Università e Ricerca through the grant Progetti Premiali 2012-iALMA ( CUP C52I13000140001 ) ; ALMA North America Development Fund; Chandra TM6-17006X. This work used the Extreme Science and Engineering Discovery Environment ( XSEDE ) , supported by NSF grant ACI-1548562, and CyVerse, supported by NSF grants DBI0735191, DBI-1265383, and DBI-1743442. XSEDE Stampede2 resource at TACC was allocated through TG-AST170024 and TGAST080026N. XSEDE JetStream resource at PTI and TACC was allocated through AST170028. The simulations were performed in part on the SuperMUC cluster at the LRZ in Garching, on the LOEWE cluster in CSC in Frankfurt, and on the HazelHen cluster at the HLRS in Stuttgart. This research was enabled in part by support provided by Compute Ontario ( http: // computeontario.ca ) , Calcul Quebec ( http: // www.calculquebec.ca ) and Compute Canada ( http: // www.computecanada.ca ) . We thank the staff at the participating observatories, correlation centers, and institutions for their enthusiastic support. This Letter makes use of the following ALMA data: ADS / JAO.ALMA # 2016.1.01154.V. ALMA is a partnership of the European Southern Observatory ( ESO; Europe, representing its member states ) , NSF, and National Institutes of Natural Sciences of Japan, together with National Research Council ( Canada ) , Ministry of Science and Technology ( MOST; Taiwan ) , Academia Sinica Institute of Astronomy and Astrophysics ( ASIAA; Taiwan ) , and Korea Astronomy and Space Science Institute ( KASI; Republic of Korea ) , in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, Associated Universities, Inc. ( AUI ) / NRAO, and the National Astronomical Observatory of Japan ( NAOJ ) .TheNRAO is a facility of the NSF operated under cooperative agreement by AUI. APEX is a collaboration between the Max-Planck-Institut für Radioastronomie ( Germany ) , ESO, and the Onsala Space Observatory ( Sweden ) . The SMA is a joint project between the SAO and ASIAA and is funded by the Smithsonian Institution and the Academia Sinica. The JCMT is operated by the East Asian Observatory on behalf of the NAOJ, ASIAA, and KASI, as well as the Ministry of Finance of China, Chinese Academy of Sciences, and the National Key R&D Program ( No. 2017YFA0402700 ) of China. Additional funding support for the JCMT is provided by the Science and Technologies Facility Council ( UK ) and participating universities in the UK and Canada. The LMT project is a joint effort of the Instituto Nacional de Astró /uniFB01 sica, Óptica, y Electrónica ( Mexico ) and the Universit The IRAM 30-m telescope on Pico Veleta, Spain is operated by IRAM and supported by CNRS ( Centre National de la Recherche Scienti /uniFB01 que, France ) , MPG ( Max-Planck-Gesellschaft, Germany ) and IGN ( Instituto Geográ /uniFB01 co Nacional, Spain ) . The SMT is operated by the Arizona Radio Observatory, a part of the Steward Observatory of the University of Arizona, with /uniFB01 nancial support of operations from the State of Arizona and /uniFB01 nancial  \nsupport for instrumentation development from the NSF. Partial SPT support is provided by the NSF Physics Frontier Center award ( PHY-0114422 ) to the Kavli Institute of Cosmological Physics at the University of Chicago ( USA ) , the Kavli Foundation, and the GBMF ( GBMF-947 ) . The SPT hydrogen maser was provided on loan from the GLT, courtesy of ASIAA. The SPT is supported by the National Science Foundation through grant PLR-1248097. Partial support is also provided by the NSF Physics Frontier Center grant PHY-1125897 to the Kavli Institute of Cosmological Physics at the University of Chicago, the Kavli Foundation and the Gordon and Betty Moore Foundation grant GBMF 947. This work made use of Global Millimeter VLBI Array ( GMVA ) data ( project MM07B; PI: A. Marscher ) . The GMVA consists of telescopes operated by the MPIfR, IRAM, Onsala, Metsahovi, Yebes, the Korean VLBI Network, the Green Bank Observatory and the VLBA. The VLBA is an instrument of the NRAO, which is a facility of the National Science Foundation operated by Associated Universities, Inc. This research has made use of data obtained with 12 radio telescopes from the East Asian VLBI Network ( EAVN ) : 2 stations from the Chinese VLBI Network ( the Tianma radio telescope operated by Shanghai Astronomical Observatory of Chinese Academy of Sciences ( CAS ) and the Nanshan radio telescope operated by Xinjiang Astronomical Observatory of CAS ) , 2 stations from the Japanese VLBI Network ( the Hitachi station operated by Ibaraki University and the Kashima station operated by the National Institute of Information and Communications Technology ) , the Korean VLBI Network operated by the Korea Astronomy and Space Science Institute ( KASI ) , the VLBI Exploration of Radio Astrometry ( VERA ) operated by the National Astronomical Observatory of Japan ( NAOJ ) , and the Nobeyama 45-meter radio telescope operated by NAOJ. We are grateful to KVN staff who helped to operate the array and to correlate the data for quasi-simulataneous multi-wavelength observations of M87. The KVN is a facility operated by KASI. The KVN operations are supported by KREONET ( Korea Research Environment Open NETwork ) , which is managed and operated by KISTI ( Korea Institute of Science and Technology Information ) . The EHTC has received generous donations of FPGA chips from Xilinx Inc., under the Xilinx University Program. The EHTC has bene /uniFB01 ted from technology shared under open-source license by the Collaboration for Astronomy Signal Processing and Electronics Research ( CASPER ) . The EHT project is grateful to T4Science and Microsemi for their assistance with Hydrogen Masers. This research has made use of NASA ' s Astrophysics Data System. We gratefully acknowledge the support provided by the extended staff of the ALMA, both from the inception of the ALMA Phasing Project through the observational campaigns of 2017 and 2018. We would like to thank A. Deller and W. Brisken for EHT-speci /uniFB01 c support with the use of DiFX. We acknowledge the signi /uniFB01 cance that Maunakea, where the SMA and JCMT EHT stations are located, has for the indigenous Hawaiian people.  \nFacility: EHT.  \nSoftware: eht-imaging , DIFMAP , SMILI .", 'Appendix A Regularizer De /uniFB01 nitions': "Here we give de /uniFB01 nitions for the regularization terms used in the RML methods that we employ in this Letter. In particular, these terms have normalization factors applied that differ among  \nthe imaging software libraries ( eht-imaging and SMILI ) .The normalization factors are generally selected so that the associated regularization hyperparameter is of order unity, and to ensure that the regularization strength is independent of, e.g., the chosen pixel size and FOV of the reconstructed image.  \nTotal Flux Density Regularizer : the total /uniFB02 ux density regularizer favors reconstructed images that have total /uniFB02 ux density near a speci /uniFB01 ed value F :  \nS I F 1 . 24 i i tot flux 2 å z =--/uni239B /uni239D /uni239C /uni239E /uni23A0 /uni239F ( )  \nIn eht-imaging , the normalization factor is  \nF . 25 2 z = ( )  \nIn SMILI , the normalization factor is  \nfF , 26 F 2 z = ( ) ( )  \nwhere f F is the target fractional error and is /uniFB01 xed to be 10 -2 for imaging in this Letter.  \nRelative Entropy : entropy ( MEM ) favors pixel-to-pixel similarity to a ' prior image ' with pixel values P i ( in this work, the prior is always taken to be a circular Gaussian with a variable FWHM ) . The corresponding regularization term is  \nS I I P 1 log . 27 i i i i MEM å z =-/uni239B /uni239D /uni239C /uni239E /uni23A0 /uni239F ( )  \nIn eht-imaging , the normalization factor ζ for MEM is  \nF , 28 z = ( )  \nwhere F is the user-speci /uniFB01 ed target total /uniFB02 ux density. In this Letter, SMILI reconstructions do not use MEM regularization.  \nℓ 1 Norm : the ℓ 1 norm favors image sparsity. In ehtimaging , it takes the form  \nS I 1 . 29 ℓ i i 1 å z =-/uni2223 /uni2223 ( )  \nwhere the normalization factor ζ is simply  \nF . 30 z = ( )  \nIn SMILI , a generalized weighted ℓ 1 norm ( ℓ 1 w ) is used, given by  \nS I P 1 , 31 ℓ i i i s 1 w  å z =-+ /uni2223 /uni2223 /uni2223 /uni2223 ( )  \nwhere Pi is a speci /uniFB01 ed prior image, and ò s is /uniFB01 xed to be 10 -12 Jy / pixel. The total scaling factor ζ is given by  \nz  \n=  \nå  \ni  \nP  \ni  \ns  \nP  \ni  \n/uni2223  \n/uni2223  \n+  \n/uni2223  \n/uni2223  \nFor the absolute value operator ... /uni2223 /uni2223 , SMILI adopts a smoothedapproximation operator, which is differentiable everywhere, given by  \nx x , 33 a 2  » + /uni2223 /uni2223 ( )  \nwhere ò a is /uniFB01 xed to be 10 -12 Jy / pixel. We use a circular Gaussian with a speci /uniFB01 ed FWFM for the prior image.  \nTotal Variation : TV favors piecewise-smooth images with /uniFB02 at regions separated by sharp edges  \nS I I I I 1 . 34 ℓm ℓ m ℓm ℓm ℓm TV , 1, , 2 , 1 , 212 å z =--+ -+ + [( ) ( ) ] ( )  \n.  \n32  \n(  \n)  \nIn this equation, the two sums are taken over the two image dimensions and the image pixels I ℓ , m are now indexed by their position ( l , m ) in a 2D grid. This regularizer is not differentiable everywhere. In eht-imaging , the normalization factor ζ for TV is  \nF , 35 beam z q = D Q ( )  \nwhere q D is the pixel size and Θ beam is the beam size.  \nOn the other hand, SMILI adopts a smoothed TV that is differentiable everywhere:  \nS I I I I 1 , 36 ℓm ℓ m ℓm ℓm ℓm TV , 1, , 2 , 1 , 2 1 2  å z = --+ -+ + + [( ) ( ) ] ( )  \nwhere ò is set to be 10 -12 Jy / pixel for imaging in this Letter. The normalization factor for SMILI is  \nF . 37 z = ( )  \nTotal Squared Variation : TSV favors images with smooth edges:  \nS I I I I 1 . 38 ℓm ℓ m ℓm ℓm ℓm TSV , 1, , 2 , 1 , 2 å z =--+ -+ + [( ) ( ) ] ( )  \nIn eht-imaging , the normalization factor for TSV is  \nF . 39 2 beam 4 z q = D Q /uni239B /uni239D /uni239C /uni239E /uni23A0 /uni239F ( )  \nIn SMILI , the normalization factor for TSV is  \nF N , 40 2 pix z = ( )  \nwhere N pix is the total number of pixels in the reconstructed image.", 'Appendix B Derivation of Constraints on the Total Compact Flux Density of M87': 'In this section, we use a combination of constraints directly accessible from EHT observations on long baselines ( Section B.1 ) in combination with constraints from multi-wavelength observations ( Section B.2 ) to constrain the compact ( /uniF088 100 /uni03BC as ) 230 GHz /uniFB02 ux density of M87 during 2017 April.', 'B.1. Constraints from EHT Observations': "As discussed in Section 4.3, there are two relevant estimates of the total /uniFB02 ux density of M87: the total /uniFB02 ux density seen in the core on arcsecond scales, F tot /uniF0A0 ≈ /uniF0A0 1.2 Jy, and the total compact /uniFB02 ux density, F cpct /uniF0A0 /uni0084 /uniF0A0 F tot, which is the integrated /uniFB02 ux density for a reconstructed image with an FOV limited to only ∼ 100 /uni03BC as. Intra-site EHT baselines are sensitive to the former, while inter-site baselines are sensitive to the latter.  \nThe shortest inter-site EHT baseline is SMT -LMT, which has a fringe spacing of ( 139 -166 ) /uni03BC as for M87. However, using this baseline to estimate F cpct requires assumptions about the compact source structure and the amplitude calibration accuracy. In particular, while SMT gains are well constrained, the use of provisional equipment at the LMT led to severe uncertainty and variability in LMT gains, as is evident by anomalously large variations in amplitudes on LMT baselines ( see Figure 2 and Paper III ) . Nevertheless, we can still obtain useful limits that are  \nindependent of source structure. To proceed, we will consider four separate constraints, each of which requires different assumptions about the data and image. Moreover, we will show that constraints on F cpct are degenerate with constraints on the overall source size. For reference, the values of F tot estimated using ALMA interferometric data are 1.13, 1.14, 1.17, and 1.21 /uniF0A0 Jy on our four observing days, April 5, 6, 10, and 11. Because of the short spacing combined with the relatively small difference in F tot for adjacent days, we will analyze each pair of consecutive observations separately, allowing for the possibility of source evolution in the four-day observing gap.  \nIn the following discussion, we will denote complex station gains by, e.g., g LMT, and measured complex visibilities by, e.g., V LMT SMT -. Apart from thermal noise, the measured visibilities are then related to the true visibilities LMT SMT  -via  \nV g g . 41 LMT SMT LMT SMT LMT SMT *  = ( ) --  \nIn expressions using visibility amplitudes, we will assume that these are suitably debiased for the contribution of thermal noise ( Equation ( 6 )) .  \nConstraint 1 . Visibility amplitudes are maximum at zero baseline length. In addition, a priori calibration may overestimate station sensitivity ( e.g., if there are unaccounted pointing errors ) , but it should not underestimate sensitivity; i.e., g 1 i  /uni2223 /uni2223 for every station i , and every visibility amplitude should only underestimate the true amplitude ( e.g., V LMT SMT LMT SMT   /uni2223 /uni2223 /uni2223 /uni2223 --) , apart from thermal noise. Combining these properties, we derive a lower limit to F cpct by identifying the maximum LMT -SMT visibility amplitude measured for each of the four observing days: 0.34, 0.23, 0.43, and 0.43 /uniF0A0 Jy, respectively. Thus, F cpct /uniF0A0 > /uniF0A0 0.34 Jy on the /uniFB01 rst two days and F cpct /uniF0A0 > /uniF0A0 0.43 Jy on the latter two days.  \nThese constraints can be strengthened under the additional assumption that the LMT -SMT baseline partially resolves the compact emission region. In this case, it is convenient to express the visibility amplitudes in terms of an equivalent Gaussian visibility function,  \nu V I I e ; , . 42 G 0 0 u 2 4ln2 q = -pq ( ) ( ) ( /uni2223 /uni2223)  \nIn this expression, I 0 is the total /uniFB02 ux density of the Gaussian source, and θ is its FWHM in radians. Visibility amplitudes will match this generic Gaussian form until a source is at least moderately resolved, providing a characteristic image FWHM that does not make speci /uniFB01 c assumptions about image structure. 120 Because LMT -SMT amplitudes are up to ∼ 20% -40% of F tot ( and F cpct may be considerably less than F tot ) ,weexpectLMT -SMT amplitudes to lie within this regime.  \nCombining these properties, we obtain a modi /uniFB01 ed constraint on total compact /uniFB02 ux density that is dependent upon the compact source size, θ cpct. For any LMT -SMT visibility amplitude on a baseline u ,  \nF V e . 43 cpct LMT SMT u cpct 2 4ln2 > p q /uni2223 /uni2223 ( ) -( /uni2223 /uni2223 )  \nConstraint 2 . Our second constraint comes from the observation that the visibility amplitudes fall steeply with  \nbaseline length on short baselines. By comparing the relative visibility amplitudes of simultaneous ALMA -LMT and SMT -LMT measurements, we obtain an estimate of the relative source visibility amplitudes that only depends on the station gains for ALMA and the SMT:  \nV V g g . 44 ALMA LMT SMT LMT ALMA SMT ALMA LMT SMT LMT   = /uni2223 /uni2223 /uni2223 /uni2223 /uni2223 /uni2223 /uni2223 /uni2223 /uni2223 /uni2223 /uni2223 /uni2223 ( ) ----  \nBoth the SMT and ALMA are stable and have wellcharacterized performance. Moreover, ALMA visibility amplitiudes have absolute calibration to within a few percent via ' network calibration, ' using the short baseline to APEX. Consequently, we can make the approximation  \nV V . 45 ALMA LMT SMT LMT ALMA LMT SMT LMT   » /uni2223 /uni2223 /uni2223 /uni2223 /uni2223 /uni2223 /uni2223 /uni2223 ( ) ----  \nFor M87, we /uniFB01 nd values for V V ALMA LMT SMT LMT /uni2223 /uni2223 /uni2223 /uni2223 --ranging from 0.18 to 0.34 on the /uniFB01 rst two days, and from 0.17 to 0.28 on the last two days.  \nFor each measurement, we can estimate a minimum image size by /uniFB01 nding the equivalent Gaussian that achieves the same quotient on the corresponding baselines. This criterion is motivated by the property that long baselines will tend to measure more /uniFB02 ux density than extrapolated by the equivalent Gaussian on short baselines, re /uniFB02 ecting small-scale power from image substructure. This criterion also makes the simplifying assumption that the source has isotropic size; this assumption is both supported by our data ( see Section 3.2 ) and is not especially problematic because the ALMA -LMT and SMT -LMT baselines sample similar orientation angles ( see Figure 1 ) .  \nWe thereby estimate a minimum source size of θ /uniF0A0 /uni0085 /uniF0A0 38 /uni03BC as for the /uniFB01 rst two days, and θ /uniF0A0 /uni0085 /uniF0A0 41 /uni03BC as for the second two days. Note that these estimates are independent of the total ( and total compact ) /uniFB02 ux density. These minimum size estimates can be applied to tighten the limits from Constraint 1, via Equation ( 43 ) .Usingthese minimum source size estimates, we obtain F cpct /uniF0A0 > /uniF0A0 0.42 Jy on the /uniFB01 rst two days and F cpct /uniF0A0 > /uniF0A0 0.56 Jy on the last two days.  \nConstraint 3 . For our third constraint, we again use the properties that LMT -SMT visibility amplitudes fall within the universal regime described via a Gaussian visibility amplitude function, and that LMT -SMT measurements will only underpredict true visibility amplitudes. Combining these properties, each measured amplitude V LMT SMT /uni2223 /uni2223 -on an associated baseline u determines an upper limit on the image size θ cpct. This upper limit occurs when the required total compact /uniFB02 ux density must exceed F tot, giving  \nu F V 2 ln2 ln . 46 cpct tot LMT SMT  q p /uni2223 /uni2223 /uni2223 /uni2223 ( ) -  \nWe thereby /uniFB01 nd 82 as  q m for the /uniFB01 rst two days, and θ /uniF0A0 /uni0084 /uniF0A0 77 /uni03BC as for the last two days.  \nConstraint 4 . Our /uniFB01 nal constraint utilizes information from different baselines that sample nearly the same vector baseline u , though not simultaneously. Speci /uniFB01 cally, the SMT -PV and LMT -PV tracks cross ( see Figure 1 ) , sampling a single baseline near ( u , v ) /uniF0A0 = /uniF0A0 ( 5.9, 0.9 ) G λ at times t 1 and t 2, respectively. Under the assumption of stable source structure, the true source visibility amplitudes at the crossing point must  \nTable 8 SMT -LMT Amplitudes before and after Applying a Correction Due to Crossing Baseline Tracks for Each Observing Day  \nNote. Non-simultaneous SMT -PV and LMT -PV visibilities at a crossing baseline location ( Near 6 G λ ; see Figure 1 ) determine the LMT SEFD at the time of the LMT -PV measurement, yielding correction factors on a priori estimates of LMT SEFDs ranging from 2.4 to 54.6. Substituting these measurements into Equation ( 48 ) gives an estimate of the LMT -SMT visibility amplitude, shown as bold values in the table. Stated uncertainties on original measurements include only thermal noise, but those on the corrected LMT -SMT amplitudes also account for uncertainties in the LMT and SMT SEFDs. On all four days, the corrected LMT -SMT visibility amplitude is approximately 400 mJy. Note that these corrections require a substantial adjustment to the LMT SEFD. In Section 8.1 we show that these large gain corrections are consistent with those recovered through imaging.  \nbe equal, providing an estimate of the relative station gains at their corresponding times:  \nV V g g g g g g . 47 t t t t t t t t t t SMT PV LMT PV SMT PV LMT PV SMT PV LMT PV SMT LMT 1 2 1 2 1 1 2 2 1 2   = » /uni2223 /uni2223 /uni2223 /uni2223 /uni2223 /uni2223 /uni2223 /uni2223 ( ) ----  \nThe approximation makes the assumptions of stable source structure between t 1 and t 2, negligible change in the visibility function on the two baselines sampled near the crossing point, and a stable gain amplitude between t 1 and t 2 at PV.  \nIf we further assume that the SMT gain is near unity ( the a priori error budget estimates g 1.00 0.04 SMT =  /uni2223 /uni2223 ) , then Equation ( 47 ) provides an estimate of the LMT gain amplitude at time t 2. We can then apply this estimate to a simultaneous measurement of the SMT -LMT amplitude at time t 2 to estimate SMT LMT  /uni2223 /uni2223 -:  \ng g V V g g V V V V V . 48 t t t t t t t t t t t t LMT SMT SMT 1 LMT 1 LMT SMT SMT PV SMT SMT LMT PV LMT SMT SMT PV LMT SMT LMT PV 2 2 2 2 1 1 2 2 2 1 2 2  = » » --/uni2223 /uni2223 /uni2223 /uni2223 /uni2223 /uni2223 /uni2223 /uni2223 /uni2223 /uni2223 /uni2223 /uni2223/uni2223 /uni2223/uni2223 /uni2223 /uni2223 /uni2223 /uni2223 /uni2223/uni2223 /uni2223 /uni2223 /uni2223 ( ) --------  \nHowever, because the gain of the LMT is substantially time variable, extrapolating the gain estimated in Equation ( 47 ) to other times ( even nearby times ) is not reliable. To summarize,  \nthis /uniFB01 nal approximation assumes the following: that the SMT -PV measurement at time t 1 samples a baseline close to the LMT -PV measurement at t 2, that the source visibility on that baseline is stable between t 1 and t 2, that the PV gain amplitude is stable between t 1 and t 2, and that the SMT gain amplitude is near unity at both times.  \nWe derive an uncertainty on t LMT SMT 2  /uni2223 /uni2223 -by summing errors in quadrature from all these effects in addition to the nominal thermal noise on each of the three visibility amplitudes. We assume a conservative uncertainty in the relative PV gain amplitude between t 1 and t 2 of 10% ( the a priori error budget for absolute gain calibration of PV is 5% ) , and we assume 3.5% uncertainty in the SMT gain amplitude ( added independently for t 1 and t 2 ) . Because we require mutual visibility of SMT and LMT at t 2, the nearest suitable points are not necessarily near the crossing point. 121 Consequently, the nearest suitable SMT -PV and LMT -PV baselines are displaced by 0.5-0.8 G λ . Empirically, we estimate up to a 20 mJy difference from this displacement ( which occurs at u 6G l » /uni2223 /uni2223 ) ; the resultant uncertainty on the estimate of t LMT SMT 2  /uni2223 /uni2223 -depends on this amplitude as a fraction of the SMT -PV and LMT -PV amplitudes.  \nTable 8 shows a triplet of visibility amplitudes and corresponding estimate of t LMT SMT 2  /uni2223 /uni2223 -for each of the four observing days. With no LMT gain correction, the measurements range from 55 -247 mJy, while after the LMT gain  \nFigure 29. Joint constraints on the size and total compact /uniFB02 ux density of 1.3 mm emission in M87. Numbered constraints ( described in Appendix B.1 ) correspond to: ( 1 ) the maximum LMT -SMT amplitude, ( 2 ) the LMT -SMT / LMT -ALMA amplitude ratio, ( 3 ) LMT -SMT amplitudes coupled with the requirement F cpct /uniF0A0 /uni0084 /uniF0A0 F tot, and ( 4 ) the SMT -LMT amplitude after LMT gain correction using the LMT -PV and SMT -PV crossing tracks in the ( u , v ) plane. Estimates on F cpct and F tot ( with ± 1 σ uncertainties ) from extrapolating VLBI observations at 3 mm are shown in brown ( see Appendix B.2 ) . The orange circle shows the estimates from earlier EHT observations by Doeleman et al. ( 2012 ) and Akiyama et al. ( 2015 ) .  \n<!-- image -->  \nFigure 28. M87 visibility amplitudes measured on the SMT -LMT baseline, as a function of baseline length. Colored circles show measurements and their associated thermal noise ( ± 1 σ ) after a /uniF0A0 priori calibration; colored diamonds show the measurements after the SMT -PV / LMT -PV crossing track correction. Uncertainties in corrected amplitudes also account for propagated a /uniF0A0 priori calibration uncertainties in the derived correction factor. Correction factors are only applied if the nearest SMT -PV / LMT -PV baseline pair are located within 1G λ . SMT -LMT amplitudes with no associated pair meeting this criterion are colored gray and have plotted uncertainties that only account for thermal noise. All uncorrected amplitudes ( circles ) have scatter far exceeding their indicated thermal noise, demonstrating that amplitude calibration uncertainties are the dominant source of error.  \n<!-- image -->  \ncorrection the measurements range from 349 -408 mJy. The corrections also all raise the SMT -LMT amplitude, as is expected. Some of these points represent the most signi /uniFB01 cant LMT amplitude losses in the entire experiment. However, Figure 28 shows the corrected SMT -LMT amplitudes for all points with a suitable crossing track correction, demonstrating that the corrected amplitudes are broadly consistent between points with large and small derived LMT correction factors.  \nIn summary, constraints 1 -3 lead to the conclusion that 0.42 /uniF0A0 /uni0084 /uniF0A0 F cpct /uniF0A0 /uni0084 /uniF0A0 1.14 Jy on April 5 and 6, and that 0.56 /uniF0A0 /uni0084 /uniF0A0 F cpct /uniF0A0 /uni0084 /uniF0A0 1.21 Jy on April 10 and 11. We also estimate that the compact source size has an equivalent Gaussian FWHM between 38 and 82 /uni03BC as on the /uniFB01 rst two days and between 41 and 77 /uni03BC as on the last two days. Constraint /uniF0A0 4 provides an absolute estimate of the LMT -SMT amplitude at a single point for each day, which allows us to provide a joint constraint on source size and F cpct. Figure 29 summarizes these constraints. Combining these constraints and marginalizing over the unknown source size, we obtain F 0.64 Jy cpct 0.08 0.39 = -+ .", 'B.2. Constraints from Quasi-simultaneous Multi-wavelength Observations': 'In addition to the non-imaging /uniFB02 ux density constraints derived from EHT data ( Section B.1 ) , we can also utilize VLBI observations at longer wavelengths to estimate expected properties at 1.3 mm. Because both the compact and extended emission from M87 are variable, these observations must be quasi-simultaneous with the EHT observations to give precise constraints. Our approach assumes that the total compact emission of M87 at 230 GHz F tot originates from the combination of ( i ) a compact emitter with unknown /uniFB02 ux density F cpct and spectral index at 230 GHz, and ( ii ) a larger and diffuse ( milliarcsecond-scale ) jet, which can be characterized by its total /uniFB02 ux density F jet /uniF0A0 = /uniF0A0 F tot /uniF0A0 -/uniF0A0 F cpct and a steep spectral index of α jet /uniF0A0 ∼ /uniF0A0 -( 0.7 -1.0 ) , as measured at cm- and mm-wavelengths ( with typical errors of ∼ 0.3; see Hovatta et al. 2014; Hada et al. 2016 ) . Then, by measuring F jet and F tot at lower frequencies and extrapolating them to 230 GHz using a power-law model, we can estimate the total compact /uniFB02 ux density F cpct at 230 GHz: F cpct /uniF0A0 = /uniF0A0 F tot /uniF0A0 -/uniF0A0 F jet.  \nTable 9 Quasi-simultaneous Lower-frequency VLBI Flux Density Measurements of M87 within a 0.5 /uniF0A0 × /uniF0A0 0.5 mas 2 Box around the Peak Flux Density  \nIn parallel with the EHT campaign, we densely monitored M87 at 22 and 43 GHz with the East Asian VLBI Network ( EAVN; see Wajima et al. 2016; Hada et al. 2017; and An et al. 2018 for more details ) , which included 12 stations from Japan, Korea, and China ( Mizusawa, Iriki, Ogasawara, Ishigaki, Hitachi, Kashima, Nobeyama, Yonsei, Ulsan, Tamna, Tianma, and Nanshan ) . For each epoch, a continuous 5 -7 hr track was observed. The angular resolution was as /uniFB01 ne as 0.55 mas at 22 GHz and 0.63 mas at 43 GHz. At both frequencies, the total /uniFB02 ux densities of the nuclear region ( i.e., F tot = F jet + F cpct ) were measured using a large window of 0.5 /uniF0A0 × /uniF0A0 0.5 mas 2 around the peak of the image and summing up the /uniFB02 ux densities of the CLEAN components within the window. From this, we obtain total /uniFB02 ux densities of F tot,22 /uniF0A0 ≈ /uniF0A0 1.3 Jy and F tot,43 /uniF0A0 ≈ /uniF0A0 1.1 -1.2 Jy at 22 and 43 GHz, respectively.  \nWe also make use of quasi-simultaneously obtained data sets of M87 from Korean VLBI Network ( KVN ) observations at 22, 43, 86, and 129 GHz on 2017 April 19 ( see Lee et al. 2016 and Kim et al. 2018b for the monitoring program, details of data analysis, typical angular resolutions, and discussions about the extended jet /uniFB02 ux ) . At 22 and 43 GHz, the correlated /uniFB02 ux densities were comparable to the results of the EAVN observations. At 86 and 129 GHz, we adopted the same window of 0.5 /uniF0A0 × /uniF0A0 0.5 mas 2 and obtained F tot,86 /uniF0A0 ≈ /uniF0A0 1.1 Jy and F tot,129 /uniF0A0 ≈ /uniF0A0 0.9 Jy, respectively. These measurements and their uncertainties are summarized in Table 9.  \nWe /uniFB01 t a power-law model to the milliarcsecond-scale /uniFB02 ux densities ( i.e., F tot tot n µ a + ) and /uniFB01 nd α tot /uniF0A0 = /uniF0A0 -( 0.14 ± 0.04 ) .By extrapolating this model to 230 GHz, we estimate F tot,230 /uniF0A0 = /uniF0A0 0.91 /uniF0A0 ± /uniF0A0 0.07 Jy. While this value is slightly lower than the ALMA interferometric measurements during the EHT observations, it only includes the total /uniFB02 ux density within the inner 0.5 /uniF0A0 × /uniF0A0 0.5 mas 2 , which is nearly three orders of magnitude below the resolution limit of ALMA.  \nNext, we estimate the compact /uniFB02 ux density F cpct,86 at 86 GHz using previous GMVA observations of M87 ( Kim et al. 2018a ) . To proceed, we compute the /uniFB02 ux densities within a much smaller window of 60 /uniF0A0 × /uniF0A0 60 /uni03BC as 2 around the peak of the GMVA images from 2004 to 2015, but excluding the 2009 epoch when the core of the jet showed exceptionally high /uniFB02 ux densities. The values of F cpct,86 are relatively stable from epoch to epoch, and therefore we adopt a value of F cpct,86 /uniF0A0 = /uniF0A0 0.55 /uniF0A0 ± /uniF0A0 0.06 Jy, where the mean and uncertainty correspond to the mean and standard deviation of the measured /uniFB02 ux densities. We can then estimate F jet at 86 GHz in 2017 April as F jet,86 /uniF0A0 = /uniF0A0 F tot,86 /uniF0A0 -/uniF0A0 F cpct,86. By propagating the errors in F tot,86 and F cpct,86, we obtain F jet,86 /uniF0A0 = /uniF0A0 0.57 /uniF0A0 ± /uniF0A0 0.18 Jy at 86 GHz.  \nFigure 30. Multi-wavelength VLBI spectrum of M87 in 2017 and expected /uniFB02 ux densities at 230 GHz ( see Section B.2 for details ) . The black line and shaded region gives the model and its ± 1 σ uncertainties for F tot.  \n<!-- image -->  \nFinally, we compute the expected F cpct,230 by estimating F jet,230 and subtracting it from the previously estimated F tot,230 /uniF0A0 = /uniF0A0 0.91 /uniF0A0 ± /uniF0A0 0.07 Jy. As explained above, we assume that F jet,230 follows a power law that can be extrapolated up to 230 GHz ( i.e., F jet,230 jet n µ a + ) . For a range of α jet /uniF0A0 = -( 0.7 -1.0 ) ,we /uniFB01 nd F jet,230 /uniF0A0 = /uniF0A0 0.25 /uniF0A0 ± /uniF0A0 0.08 Jy. This value gives F cpct,230 /uniF0A0 = /uniF0A0 0.67 /uniF0A0 ± /uniF0A0 0.10 Jy within the central ∼ 60 /uniF0A0 × /uniF0A0 60 /uni03BC as 2 . Figure 30 summarizes the spectrum of M87 at lower frequencies and the expected /uniFB02 ux densities at 230 GHz. We emphasize that this estimate for F cpct,230 represents a characteristic value, and that the exact value may have varied slightly over the EHT observing window.  \nWe note that this value for F cpct,230 is consistent with the constraints derived in Section B.1 using EHT measurements. Combining them, we obtain F 0.66 Jy cpct,230 0.10 0.16 = -+ .', 'Appendix C Synthetic Data': 'To test and calibrate the parameters that we use in our /uniFB01 nal image reconstructions from imaging parameter surveys ( Section 6 ) , we generate synthetic data from different geometrical models and GRMHD simulations using two independent pipelines ( as introduced in Section 6.1 ) . We consider four geometric source models with parameters chosen such that their visibility data on EHT baselines match salient features of the observations, in particular the visibility nulls at ∼ 3.4 and ∼ 8.3 G λ .', 'C.1. Geometric Models': 'We use four geometric models with different compact structure ( displayed in Figure 5 ) . The model parameters in each case were chosen such that their visibility data on EHT baselines match salient features of the observations ( e.g., the visibility null at ∼ 3.4 G λ ) .  \nAll of our models have 0.6 Jy of /uniFB02 ux density in a compact component and 0.6 Jy in an extended jet that is heavily resolved by EHT inter-site baselines. The large-scale jet model is common to all images and consists of three elliptical Gaussian components aligned to images of the M87 jet at 3 mm ( e.g., Kim et al. 2018a ) . The jet parameters are displayed in Table 10.  \nWe used four different models for the compact /uniFB02 ux density. These are as follows.  \nTable 10 Gaussian Parameters of the Large-scale Jet Model used in Synthetic Data Generation  \n- 1. A delta-function ring with radius r 0 /uniF0A0 = /uniF0A0 22 /uni03BC as, convolved with a circular Gaussian with FWHM of 10 /uni03BC as.\n- 2. A crescent composed of a delta-function ring with radius r 0 /uniF0A0 = /uniF0A0 22 /uni03BC as and a dipolar angular intensity pro /uniFB01 le  \nIr I A r r r , 1 2 cos 2 , 49 0 0 0 q q h d p = + --( ) [ ( )] ( ) ( )  \nwhere I 0 /uniF0A0 = /uniF0A0 0.55 Jy, the asymmetry parameter A /uniF0A0 = /uniF0A0 0.23 ( see Section 9.1 ) , and η is the orientation angle east of north. The crescent is then also convolved with a 10 /uni03BC as Gaussian.  \n- 3. A uniform disk with radius r /uniF0A0 = /uniF0A0 35 /uni03BC as, convolved with a circular Gaussian with FWHM of 10 /uni03BC as.\n- 4. Two circular Gaussian components, each with a FWHM of 20 /uni03BC as. The /uniFB01 rst component is located at the origin and has a total /uniFB02 ux density of 0.27 Jy, while the second is located at /uni0394 R.A. /uniF0A0 = /uniF0A0 30 /uni03BC as and /uni0394 decl. /uniF0A0 = /uniF0A0 -12 /uni03BC as and has a total /uniFB02 ux density of 0.33 Jy.  \nBoth the jet and compact structure are polarized in order to produce non-closing systematic errors from polarimetric leakage. The jet components all have a constant fractional polarization of 20% and a random position angle that varies across the image with a coherence length of 400 /uni03BC as. The compact components all have a constant fractional polarization of 40% and a random polarization position angle with a coherence length of 5 /uni03BC as.', 'C.2. Synthetic Data Generation with eht-imaging': 'The eht-imaging library produces synthetic visibilities from model images by computing their Fourier transform, adding random thermal noise ( Equation ( 3 )) , and corrupting the data with systematic station-based effects using the Jones matrix formalism ( Thompson et al. 2017 ) . Off-diagonal terms in the Jones matrices mix the measured polarizations, which in the absence of noise and systematic error would correspond simply to the Fourier transform of the image polarizations by Equation ( 2 ) .  \nWe take circular polarizations as our primary basis, and consider a single baseline ij . The Fourier transform of the image gives the uncorrupted visibility of each polarization ( RR ij ¢ , LL ij ¢ , LR ij ¢ , RL ij ¢ ) , which are assembled into a 2 /uniF0A0 × /uniF0A0 2 correlation matrix:  \nV RR RL LR LL . 50 ij ij ij ij ij ¢= ¢ ¢ ¢ ¢ /uni239B /uni239D /uni239C /uni239C /uni239E /uni23A0 /uni239F /uni239F ( )  \nThe thermal noise standard deviation on each baseline and polarization is given by Equation ( 3 ) and assembled into a matrix ij s of the same form.  \nTable 11 Field Rotation Parameters for the EHT Stations  \nIncluding the effects of systematic and thermal noise, the simulated visibility matrix is V ij :  \nV J V J , 51 ij i ij j ij  s = ¢ + ( ) ( ) †  \nwhere J i and J j are the station Jones matrices. Each Jones matrix has the form  \nJ g e d g e d g g , 52 i R R R i L L L = j j -/uni239B /uni239D /uni239C /uni239C /uni239E /uni23A0 /uni239F /uni239F ( )  \nwhere g R, g L are the complex gain terms, d R and d L are the constant station d -terms, and j is a term from /uniFB01 eld rotation. We typically expect g /uniF0A0 < /uniF0A0 1 due to losses in telescope sensitivity.  \nThe complex gain terms include an absolute amplitude gain offset g /uni2223 /uni2223 , and a random atmospheric phase f . For the purposes of this Letter, we set g R /uniF0A0 = /uniF0A0 g L because the atmosphere is not signi /uniFB01 cantly birefringent at millimeter wavelengths and because differential gains between right and left polarizations are removed in the calibration process ( Paper III ) :  \ng g g t e . 53 i t R L = = f /uni2223 ( )/uni2223 ( ) ( )  \nThe phase error f is drawn from a uniform distribution once per scan; it is kept stable within scans to mimic the ad hoc phase calibration procedure. The amplitude gain of a site at each time has a time-stable component g 1 and a time-varying component g 2 ( t ) :  \ngt g g t 1 . 54 1 2 1 2 = -+ ( ) /uni2223 /uni2223 /uni2223 ( )/uni2223 ( )  \nBoth g 1 and g 2 are sampled from normal distributions with zero mean: g 0, g 1 1  s ~ ( ) and g t 0, g 2 2  s ~ ( ) ( ) . To simulate its higher median gain error and variability, for the LMT we set 0.6 g 1 s = and 0.5 g 2 s = . For all other sites we set 0.15 g 1 s = and 0.05 g 2 s = . g 2. Like the atmospheric phase, we /uniFB01 xed the gain errors so that g 2 ( t ) changes between scans but is kept stable across scans.  \nThe complex d terms, d R and d L, are stationary in time and drawn from a complex normal distribution with zero mean and a standard deviation of 0.05 in each of the real and imaginary parts ( motivated by the estimates in Johnson et al. 2015 ) . The mean amplitude of the leakage terms is then 6.3%. The /uniFB01 eld rotation phase term j has three possible contributions depending upon the receiver mount:  \nf f , 55 el el par par off j q y j = + + ( )  \nwhere θ el is the elevation angle, ψ par is the parallactic angle, and j off is a constant offset. Cassegrain mounts have f par /uniF0A0 = /uniF0A0 1 and f el /uniF0A0 = /uniF0A0 0. Nasmyth mounts have f par /uniF0A0 = /uniF0A0 1 and f el /uniF0A0 = /uniF0A0 ± 1, depending on the handedness. The EHT station /uniFB01 eld rotation parameters are listed in Table 11.  \nAfter being simulated with Equation ( 51 ) , the synthetic data were network calibrated and converted to Stokes I , following the procedure used for the EHT data pipeline in Paper III.', 'C.3. Synthetic Data Generation with MeqSilhouette and rPICARD': 'An mm-VLBI synthetic data pipeline built on the synthetic data generator MeqSilhouette ( Blecher et al. 2017 ) and the EHT CASA calibration pipeline rPICARD ( Janssen et al. 2019 ) has been developed ( F. Roelofs et al. 2019, in preparation ) . Like eht-imaging , this pipeline is able to produce synthetic data sets with statistics that closely match real VLBI data, but it is done using an a priori physics-based approach instead of an a posteriori data-based approach.  \nMeqSilhouette ( Blecher et al. 2017 ) is a synthetic VLBI data generator that predicts visibilities from both parametric and non-parametric sky models using MeqTrees ( Noordam & Smirnov 2010 ) and WSClean ( Offringa et al. 2014 ) , respectively, and applies physically motivated propagation and instrumental effects. Uncorrupted frequency-resolved visibilities are generated from a full-Stokes sky model, antenna locations, and a real observing schedule. Elevation-dependent mean tropospheric delays, opacity, and sky temperatures are simulated from input weather parameters ( ground pressure P g,  \nground temperature T g, and precipitable water vapor ) using the ATM software ( Pardo et al. 2001 ) . In addition, turbulent phases are simulated based on input atmospheric coherence times t c and assuming a thin scattering screen model characterized by Kolmogorov turbulence. Station-dependent thermal noise drawn from a Gaussian distribution is included and contains contributions from both the receiver and the atmosphere. Scanbased antenna pointing errors are added using the primary beam sizes at 230 /uniF0A0 GHz and observatory-measured pointing uncertainties. Full polarization simulation of time-variable sources is also possible ( as will be relevant to Sgr A * observations ) , with the option of introducing station-based, frequency-dependent polarization leakage and complex gain errors.  \nIn the MeqSilhouette+rPICARD pipeline ( F. Roelofs et al. 2019, in preparation ) , the full observation and calibration process is simulated by running the synthetic data generated by MeqSilhouette through the CASA VLBI calibration pipeline rPICARD ( Janssen et al. 2019 ) . In this calibration process, fringe /uniFB01 tting, amplitude calibration, and time and frequency averaging are done in the same way as for the real EHT data.  \nFigure 31 shows a comparison between synthetic data generated from the same model image using eht-imaging and  \nFigure 31. Comparison of the synthetic data from the eht-imaging and MeqSilhouette data generation libraries. Low-band data were generated from the same source image. The reconstructed images, visibility amplitudes, and S / Ns from the data generated with the MeqSilhouette+rPICARD pipeline are broadly consistent with those generated with eht-imaging , despite the different approaches of the software packages. Points are colored by baseline as in Figure 18. Images are reconstructed using the same script in eht-imaging using the /uniFB01 ducial parameters from Section 6.  \n<!-- image -->  \nThe Astrophysical Journal Letters, 875:L4 ( 52pp ) , 2019 April 10  \nTable 13  \nNote. The parameters PWV, P g, and T g are median values for the EHT 2017 campaign as logged in VLBImonitor ( Paper II ) .  \nthe MeqSilhouette+rPICARD simulation software. For the MeqSilhouette simulations, input weather parameters ( Table 12 ) were based on site measurements logged in VLBImonitor ( Paper II ) during the EHT campaign. Coherence times, SEFDs, and pointing offsets were based on measurements from individual EHT stations. A leakage of 5% was added to the visibilities with the /uniFB01 eld rotation parameters listed in Table 11, and an overall gain factor of 0.6 was set for the LMT.  \nThe reconstructed images ( generated with the M87 /uniFB01 ducial parameters of the eht-imaging imaging pipeline ) , visibility amplitudes, and S / Ns from the data generated with the MeqSilhouette and rPICARD pipeline are broadly consistent with those generated with eht-imaging , despite the different approaches used in the two software packages. Table 13 shows the normalized cross-correlation ρ NX and equivalent blurring α ( given in /uni03BC as ) , for images made from the eht-imaging and MeqSilhouette+rPICARD data simulation software. Each software package was used to generate three different realizations of visibilities from the same underlying crescent geometric model ( see Figure 5 ) . The ρ NX and α computed from reconstructions of this data show broad consistency. This consistency reinforces our con /uniFB01 dence in both  \nImage Fidelity Metrics, Normalized Cross-correlation ρ NX, and Equivalent Blurring α ( given in /uni03BC as ) , for Images Made from Data Simulated using eht-imaging and MeqSilhouette+rPICARD  \nNote. Images from three sets of visibilities of the crescent ( + jet ) model, with different realizations of noise ( samples 1 -3 ) , were generated for each data simulation software package using the /uniFB01 ducial parameters from the ehtimaging imaging pipeline. The ρ NX and α /uniFB01 delity metrics are insensitive to both the particular noise realization and the simulation software. Thus, our selected Top Sets and /uniFB01 ducial parameters are insensitive to the particular noise realization and the simulation software.  \npipelines, particularly in corroborating that the image /uniFB01 delity is primarily determined by the limited baselines coverage rather than data sensitivity or calibration limitations.', 'Appendix D Image Dependence on Sites for SMILI and DIFMAP Pipelines': 'In Section 8.2, we use the eht-imaging pipeline with /uniFB01 ducial parameters to show that the ring structure in reconstructed images is resilient to the loss of most sites. Figure 32 shows the corresponding reconstructions from the other two pipelines on April /uniF0A0 11. The ring structures produced by all three pipelines are similarly sensitive to dropping most individual sites, with Chile and PV being the most critical. The image degradation due to the loss of all data to a site tends to be larger for RML approaches, as these methods use weaker constraints on the allowed FOV for image /uniFB02 ux.  \nFigure 32. Example reconstructions of M87 on April /uniF0A0 11 after omitting visibilities to each geographical site for the DIFMAP and SMILI pipelines. Following Figure 22, the top row of each pipeline shows reconstructions that exclude all baselines to the indicated site, while the bottom row shows reconstructions that exclude all visibility amplitudes on baselines to the indicated site but retain closure phases and closure amplitudes. The lower-left panel shows an image reconstructed after excluding visibility amplitudes to all sites ( i.e., using only closure quantities ) . Because DIFMAP uses only visibilities, we only show the case of excluding baselines. The ellipse in each panel shows the corresponding synthesized beam with uniform weighting; only the DIFMAP images are restored with these beams.  \n<!-- image -->', 'Appendix E Interday Variability in M87 Images': 'The closure phases of M87 show variation across the four observations ( Paper III ) . Similarly, our /uniFB01 ducial images show corresponding variation in image structure between April 5, 6 and 10, 11 ( Figure 11 ) . These differences are also apparent in the averaged /uniFB01 ducial images restored to a common, conservative resolution ( see Figure 15 ) . In Section 8.3, we demonstrated that the image variability is not simply from changes in the baseline coverage. We now assess this image variation.  \nFigure 33 shows differences of reconstructed image pairs on different days. To suppress features that are speci /uniFB01 c to imaging parameter choices, we show the beam-convolved difference of two mean images, where the mean is taken over all corresponding Top Set images for the speci /uniFB01 ed pipeline and  \nday. For all pipelines, but especially for eht-imaging and SMILI , the differences are larger across the multi-day gap than across adjacent days, as expected if these differences re /uniFB02 ect source evolution.  \nFor both eht-imaging and SMILI , the difference images indicate a shift of the /uniFB02 ux from the left to the right side of the ring ( east to west ) over the observations, which is consistent with the slight increase in η over the days ( Section 9 ) . In contrast, DIFMAP difference images are dominated by largerscale morphology changes. For all pipelines, the beamconvolved temporal differences are nearly an order of magnitude larger than the median standard deviation of differences seen across eht-imaging Top Set images ( 1.6 × 10 8 K; see Figure 17 ) . Nevertheless, in all cases, we cannot unambiguously associate the variability with smooth motion or evolution of a speci /uniFB01 c image feature.  \nFigure 34 compares the interleaved multiplicative station gains for M87 and 3C 279 for each day; Table 14 presents the gain statistics for the two sources and compares them to the  \n<!-- image -->  \nFigure 33. Variations seen in reconstructed images between days. Each row shows the beam-convolved difference of the aligned mean images from the Top Set of the corresponding parameter survey. All images are shown using the same color scale. Positive values ( red ) indicate an increase in brightness at the later date, and negative values ( blue ) indicate a decrease in brightness at the later date.', 'Appendix F Multi-day Gain Comparisons with 3C 279': 'In Section 8, we compared the residual amplitude gains for the LMT and SMT on April /uniF0A0 5 derived from separately imaging M87 and 3C 279. We now expand this comparison to all four days on which these sources were observed.  \nexpected amplitude calibration error budget ( Paper III ) . The trends of inferred residual gain amplitudes are consistent among the imaging pipelines, and those of M87 are similar to those of 3C 279. The inferred gains for M87 at the LMT have large excursions that are not seen in the 3C 279 gain trends; these excursions are likely due to poor pointing and tracking. These large LMT correction factors ( up to g 6 LMT 1 » -/uni2223 /uni2223 ) are also consistent with those estimated using the SMT -PV and LMT -PV crossing tracks ( see Table 8 and Figure 28 ) .  \nFigure 34. Derived residual gains for the SMT ( left ) and LMT ( right ) using self-calibration to images of M87 ( red ) and 3C 279 ( blue ) . The M87 images are the /uniFB01 ducial images from each pipeline; the 3C 279 images were reconstructed separately, using adapted imaging scripts. The particularly large excursions on the LMT M87 gains are often due to poor pointing. For instance, excursions at ≈ 6 UTC are from dif /uniFB01 culties tracking the source during transit ( ∼ 81 ° elevation ) . Note the different ranges shown for each observing day.  \n<!-- image -->  \nTable 14 Residual Gain Corrections Derived for M87 and 3C 279 for Observations on All Days ( April 5, 6, 10, 11 )  \nTable 14 ( Continued )  \nNote. Numbers indicate the median of g 1 /uni2223 /uni2223 , with stated uncertainties corresponding to the 25th and 75th percentiles. The multiplicative gains on the visibilities for each station are computed by self-calibrating the network-calibrated data sets to the reconstructions after rescaling the intra-site baseline amplitudes to match the total /uniFB02 ux of the reconstructions. The percentage deviation from unity can be compared to the expectations from the station-based a priori error budget on the visibility amplitudes derived from Paper III. The SPT is omitted from gain comparisons on the /uniFB01 rst two days as it did not observe 3C 279 on April 5, and it observed only a single 3C 279 scan on April /uniF0A0 6.', 'Appendix G Finite Resolution Bias on Ring Parameters': "In addition to uncertainties from thermal noise, systematic noise, and algorithmic imaging assumptions, estimated image properties are necessarily limited by the image resolution. Image structure at scales /uniFB01 ner than the diffraction-limited resolution can bias properties such as the magnitude and location of maximum image brightness.  \nAs a simple example, we will consider the /uniFB01 nite resolution bias for a thin ring. Written in cylindrical coordinates ( r , θ ) ,the image of an in /uniFB01 nitesimally thin ring of diameter d takes the form  \nI r d d r d , ; 1 2, 56 ring q p d = -( ) ( ) ( )  \nwhere the image is normalized to have a total /uniFB02 ux density of unity. Next, suppose that this image is convolved with a circular Gaussian kernel having FWHM α . The resulting blurred ring is given by  \nI r d d rdrd r d r r rr r d I rd , ; , 4ln2 2 exp 4ln2 2 cos 4ln2 exp 4ln2 4 4ln2 , ring 2 2 2 2 2 2 2 2 2 0 2 ò q a pa qd a q pa a a = ¢ ¢ ¢ ¢ -' -+ ¢ -¢ ¢ = -+ ' /uni239C /uni239F /uni23A1 /uni23A3 /uni23A2 /uni23A4 /uni23A6 /uni23A5 /uni23A1 /uni23A3 /uni23A2 /uni23A4 /uni23A6 /uni23A5 /uni239B /uni239D /uni239E /uni23A0 ( ) ( ) ( ) ( )  \nwhere I 0 ( x ) in the /uniFB01 nal line is a modi /uniFB01 ed Bessel function of the /uniFB01 rst kind. In the limit α /uniF0A0 = /uniF0A0 d , we can simplify this expression using the identity I x 1 e x x 0 2 1 x  = + p /uni23A1 /uni23A3 /uni23A4 /uni23A6 ( ) ( ) ( Abramowitz & Stegun 1972 ) . Substituting r d 2 D = -and only keeping  \nterms to leading order in /uni0394 / d and α / d , we obtain  \nI r d d d e ; , 4ln2 1 1 . 57 ring 3 4ln2 2 a p a » -D a -D /uni23A1 /uni23A3 /uni23A2 /uni23A4 /uni23A6 /uni23A5 ( ) ( ) ( )  \nThus, when α /uniF0A0 = /uniF0A0 d , the FWHM of the blurred ring is approximately equal to α . More generally, convolving a ring that has a Gaussian pro /uniFB01 le of intrinsic FWHM w true with a circular Gaussian kernel with FWHM α then gives an effective width  \nw w . 58 meas true 2 2 a » + ( )  \nLikewise, we can use Equation ( 57 ) to estimate the radius of peak brightness after blurring; to leading order in α / d , r d 2 d pk 1 8ln2 2 » -a . Thus, letting d true denote the angular diameter of the unconvolved ring and d meas denote the angular diameter of peak brightness of the convolved ring, we obtain  \nd d 1 , 59 d true meas 1 4ln2 2 meas » -a ( ) ( )  \nwhere the approximation is accurate to leading order in α / d meas.  \nAs a concrete example, choosing parameters that are similar to our M87 measurements in Section 9.3,anin /uniFB01 nitesimally thin ring with true diameter d /uniF0A0 = /uniF0A0 44 /uni03BC as and width α /uniF0A0 = /uniF0A0 15 /uni03BC as will have its measured diameter biased downward by ≈ 2 /uni03BC as.  \nWhile the simple Gaussian convolution assumed in this example only crudely approximates the effects of /uniFB01 nite resolution on reconstructed images, it indicates that for images dominated by a thin ring, estimated diameters ( widths ) will be biased downward ( upward ) by the /uniFB01 nite resolution of an image  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \nFigure 35. From left to right, measured diameter d , width w , orientation angle η , and central brightness ratio f C measured from the April /uniF0A0 11 /uniFB01 ducial images blurred with circular Gaussian kernels of increasing FWHM α . The solid lines indicate the measured value, and the shaded regions give the 1 σ uncertainty as de /uniFB01 ned in Section 9.1 ( these uncertainties do not include a contribution from scatter across the Top Set ) . The dashed line on the /uniFB01 rst panel shows the prediction of Equation ( 59 ) for the measured diameter d meas as a function of the blur kernel assuming d true /uniF0A0 = /uniF0A0 42 /uni03BC as. The dashed line in the second panel shows the FWHM of a 15 /uni03BC as Gaussian convolved with the kernel of FWHM α ( Equation ( 58 )) . Because the DIFMAP images are fundamentally composed of point sources, the measurements from these images become highly uncertain when α /uniF0A0 /uniF088 /uniF0A0 10 /uni03BC as.  \n<!-- image -->  \nreconstruction. For the /uniFB01 ducial image measurements in Section 9, the diameter bias from this /uniFB01 nite-resolution effect would be maximum in the case of an in /uniFB01 nitesimal intrinsic ring, with α then approximately corresponding to the measured width. Thus, the /uniFB01 nite-resolution diameter bias is at most a few /uni03BC as.  \nIn Figure 35, we explore the dependence of estimated ring parameters on the size of an applied restoring beam. Using the /uniFB01 ducial M87 images on April /uniF0A0 11 from all three imaging pipelines, we show parameters computed after convolving each image with a circular Gaussian kernel having FWHM α in the range 0 -20 /uni03BC as. The DIFMAP images are fundamentally composed of point sources ( ' CLEAN ' components ) . When these point sources are restored with α /uniF0A0 /uniF088 /uniF0A0 10 /uni03BC as, the CLEAN components still appear in the convolved image as individual point sources. As a result, the feature extraction methods of Section 9, which assume a smooth ring structure, have large uncertainties when extracting parameters on these images. This uncertainty is particularly apparent in the orientation angle measurement. The RML images, by contrast, have a /uniFB01 nite width and smooth structure even at α /uniF0A0 = /uniF0A0 0. As a result, their parameters vary smoothly and have similar uncertainties at all values of α .  \nThe leftmost panel of Figure 35 shows that the dependence of the measured diameter on α closely follows Equation ( 59 ) . In addition, the /uniFB01 ducial images have a scatter of approximately 1 /uni03BC as across the different imaging methods, but an uncertainty of ≈ 3 -4 /uni03BC as across the varying α . Thus, when extracting physical parameters from the measured diameter, it is important to take into consideration the additional bias and uncertainty induced by this effect.  \nFor small values of the restoring beam α , Figure 35 shows that the measured width w of the RML reconstructions follows the prediction of Equation ( 58 ) . For larger α /uniF0A0 > /uniF0A0 15 /uni03BC as, the kernel size approaches the ring radius, and higher-order effects become important ( i.e., contributions from the opposite side of the ring in the convolved width ) . Because it is intrinsically built of point sources that are not con /uniFB01 ned to a δ -function in radius, the DIFMAP image does not follow this simple prediction for the increase in width with blurring kernel size, and it instead increases more rapidly to converge with the RML result at α /uniF0A0 ≈ /uniF0A0 20 /uni03BC as. For the RML reconstructions ( and for DIFMAP  \nwith α > 10 /uni03BC as ) , the measured orientation angle is relatively unaffected by the Gaussian convolution.  \nOf all the parameters de /uniFB01 ned in Section 9.1, the fractional central brightness f C between the average ring center brightness and the rim varies the most with resolution. In the absence of a restoring beam, both SMILI and DIFMAP produce rings with practically zero brightness in the ring center; as a result of this near-zero /uniFB02 oor, f C is extremely small ( f C < 10 -5 ) . As convolution with a /uniFB01 nite Gaussian kernel /uniFB01 lls in the center of the ring, f C increases rapidly with α , by several orders of magnitude, for these images. By contrast, because they include inverse-tapering of the data and a /uniFB01 nal blurring with a 5 /uni03BC as Gaussian, the eht-imaging /uniFB01 ducial reconstructions always have non-zero central brightness. All three imaging methods give f C /uniF0A0 ≈ /uniF0A0 0.3 for α /uniF0A0 = /uniF0A0 20 /uni03BC as. Because f C is sensitive to resolution and imaging methodology, we can at only securely conclude that f C /uniF0A0 /uniF088 /uniF0A0 0.3 for M87.", 'Appendix H Image Dependence on the Assumed Total Compact Flux Density': 'While all three parameter surveys explored different choices of total compact /uniFB02 ux density from 0.4 to 0.8 Jy, they were tested against simulated images that all had F cpct /uniF0A0 = /uniF0A0 0.6. For M87, the true value of F cpct is unknown, adding additional uncertainty that is not explored systematically in the parameter surveys or Top Sets. In Figure 36, we explore how the estimated ring parameters of M87 from reconstructed images depend upon the assumed total compact /uniFB02 ux density. In this analysis, we /uniFB01 x all other imaging parameters to their /uniFB01 ducial values.  \nNot all values of F cpct produce acceptable images, with two primary problems. First, a reconstructed image may not be in acceptable agreement with the underlying data. Second, a reconstructed image may have corresponding station gains from self-calibration that are signi /uniFB01 cantly different than a priori expectations. The SMT and LMT gains are most sensitive to F cpct because the SMT -LMT baseline is the only inter-site EHT baseline that weakly resolves M87. To de /uniFB01 ne acceptable images in this test, we require two conditions to be satis /uniFB01 ed ( similar to the requirements for the Top Sets; Section 6.3.1 ) .First,werequirethat , 2 CP 2 log CA 2 c c < with no systematic uncertainty for eht-  \nFigure 38 shows the angular pro /uniFB01 les along the ring from the three imaging pipelines on April /uniF0A0 11 data. While Figure 37 showed variability in the radial pro /uniFB01 les from the individual /uniFB01 ducial images, Figure 38 considers variability in the azimuthal structure across the entire Top Set of images. In each panel, the solid line represents the /uniFB01 ducial value of the angular pro /uniFB01 le, and the bands show variability across the entire Top Set.  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \nFigure 36. Measured diameter d , width w , orientation angle η , and central brightness ratio f C of M87, measured from image reconstructions with varying total compact /uniFB02 ux density, F cpct. All other imaging parameters were set to the /uniFB01 ducial parameters of the corresponding pipeline. DIFMAP values were measured after restoring with a 10 /uni03BC as FWHM Gaussian beam. The solid lines indicate the measured value, and the shaded regions give the ± 1 σ uncertainty as de /uniFB01 ned in Section 9.1 ( these uncertainties do not include a contribution from scatter across the Top Set ) . The plotted values are dark in the range of F cpct that produces images that have ( 1 ) /uniF0A0 , 2 CP 2 log CA 2 c c < ( with 5% systematic uncertainty for DIFMAP ) , and ( 2 ) /uniF0A0 a corresponding self-calibration solution with g 0.9 median 1 1.1 SMT   ( /uni2223 /uni2223 ) . The vertical dashed lines indicate the 0.56 -0.82 Jy constraint on the total compact /uniFB02 ux density from pre-imaging considerations ( Section 4 ) ; the vertical solid lines indicate the 0.4 -0.8 Jy range explored in the parameter surveys ( Section 6 ) .  \n<!-- image -->  \nimaging /uniF0A0 and SMILI /uniF0A0 and 5% systematic uncertainty for DIFMAP . Second, we require that the median SMT gain amplitude satis /uniFB01 es g 0.9 median 1 1.1 SMT   ( /uni2223 /uni2223 ) ( we only assess the images based upon g SMT /uni2223 /uni2223 because the LMT a /uniF0A0 priori calibration is comparatively poor; see Appendix F ) .For eht-imaging /uniF0A0 and SMILI , the SMT gain constraint gives 0.40 Jy /uniF0A0 /uni0084 /uniF0A0 F cpct /uniF0A0 /uni0084 /uniF0A0 0.81 Jy and 0.40 Jy /uniF0A0 /uni0084 /uniF0A0 F cpct /uniF0A0 /uni0084 /uniF0A0 0.85 Jy, respectively; all reconstructions across these ranges satisfy the χ 2 criterion. For DIFMAP ,theSMT gain constraint gives 0.43 Jy /uniF0A0 /uni0084 /uniF0A0 F cpct /uniF0A0 /uni0084 /uniF0A0 0.91 Jy; however, only images with 0.18 Jy /uniF0A0 /uni0084 /uniF0A0 F cpct /uniF0A0 /uni0084 /uniF0A0 0.56 Jy satisfy the χ 2 criterion. Thus, for DIFMAP , the range of images meeting both criteria is 0.43 Jy /uniF0A0 /uni0084 /uniF0A0 F cpct /uniF0A0 /uni0084 /uniF0A0 0.56 Jy. For all pipelines, the acceptable range of F cpct extends lower than that derived in Section 4 because we do not strictly enforce g 1 1 SMT  /uni2223 /uni2223 , instead enforcing the weaker constraint that the g median 1 0.9 SMT  ( /uni2223 /uni2223 ) .  \nFor DIFMAP , the ring diameter is relatively insensitive to F cpct, although the azimuthal variations σ d increase signi /uniFB01 -cantly for F cpct /uniF0A0 /uniF089 /uniF0A0 0.55. For eht-imaging and SMILI , larger values of F cpct lead to broader rings with slightly lower diameters, as expected from the /uniFB01 nite resolution bias discussed in Appendix G. Across the range of F cpct identi /uniFB01 ed in Appendix B, the most signi /uniFB01 cant variations are in the orientation angle η estimated with the SMILI pipeline, and in the central brightness ratio f C, which increases steeply with F cpct for all pipelines.', 'Appendix I Image Radial and Azimuthal Pro /uniFB01 les': "In this section, we present additional radial and azimuthal pro /uniFB01 les of the /uniFB01 ducial M87 images. These pro /uniFB01 les use the ring centroids determined in Section 9.  \nFigure 37 shows radial pro /uniFB01 les taken across the rings identi /uniFB01 ed in the three April 11 /uniFB01 ducial images, with the DIFMAP image restored by the nominal 20 /uni03BC as beam. For each image, we identi /uniFB01 ed the orientation angle η , then divided the ring in two by the line perpendicular to the measured orientation angle. On each half-ring, we compute the median angular pro /uniFB01 le ( solid lines in Figure 37 ) , as well as the 25% -75% and 0% -100% percentile ranges of variation. Figure 37 shows that the peak-to-peak ring diameters measured from the three imaging methods are broadly consistent, as indicated by  \nour measurements in Figure 25. The overall sense of ring asymmetry recovered from each imaging pipeline is also consistent; the ring is always brighter in the south.  \nHowever, the shape of the median radial pro /uniFB01 les differ among the different imaging methods due to the different assumptions made in the imaging process. The DIFMAP reconstructions in Figure 37 are restored with a 20 /uni03BC as beam, so they produce wider, shallower radial pro /uniFB01 les. In contrast, the SMILI images tend to zero out low brightness regions, and the /uniFB01 ducial images have a near-zero brightness depression in the center of the ring. The eht-imaging results are of a higher resolution than the DIFMAP results, but they have a less dramatic brightness depression that the SMILI images as a result of the choice to include a 5 /uni03BC as maximum resolution in the eht-imaging script and the use of entropy regularization.  \nThe DIFMAP angular pro /uniFB01 le in Figure 38 is smooth and shallow, and it is distinct from those of the RML methods due to its 20 /uni03BC as restoring beam. The dashed lines in the ehtimaging and SMILI panels show the angular pro /uniFB01 les from the /uniFB01 ducial image blurred to match the DIFMAP resolution. When blurred, the angular pro /uniFB01 les from the eht-imaging and SMILI images better match the broad DIFMAP pro /uniFB01 le, but they still differ in the measured orientation angle and the position of the brightest location on the ring.  \nThe angular pro /uniFB01 les for the unblurred eht-imaging and SMILI reconstructions are similar, with the 0.1 Jy difference in total /uniFB02 ux density in these reconstructions manifesting in an overall lower pro /uniFB01 le for the SMILI reconstruction. The ' knot ' features in the eht-imaging and SMILI reconstructions show signi /uniFB01 cant variation across the Top Sets. This suggests that these ring features are highly sensitive to imaging parameter choices, and that they are potentially artifacts of the limited ( u , v ) coverage. In general, the azimuthal structure  \nFigure 37. 1D radial brightness pro /uniFB01 les of the three /uniFB01 ducial M87 images on April /uniF0A0 11. For each image, radial pro /uniFB01 les in the semi-circle centered on η are plotted with negative values of r , and radial pro /uniFB01 les through the opposing semi-circle centered on η /uniF0A0 + /uniF0A0 180 ° are plotted with positive r . The solid curves show the median pro /uniFB01 le over the corresponding semi-circle, the darker band shows the 25th to 75th percentile range, and the lighter band shows the full range of pro /uniFB01 les in the /uniFB01 ducial images.  \n<!-- image -->  \nFigure 38. 1D angular pro /uniFB01 les of M87 on April /uniF0A0 11 from the three imaging methods. For each method, the solid line shows the angular pro /uniFB01 le obtained from the /uniFB01 ducial image, the darker band shows the 25th to 75th percentile range across the Top Set, and the lighter band shows the full Top Set range. For the RML methods, the dashed line shows the angular pro /uniFB01 le from the /uniFB01 ducial image blurred to the resolution of the DIFMAP image. The ' knot ' features in the unblurred eht-imaging and SMILI reconstructions show signi /uniFB01 cant variation across the Top Set, suggesting that these features are sensitive to imaging parameter choices.  \n<!-- image -->  \nin the reconstructed images from all three imaging pipelines is more variable than the radial structure, making measurements of the orientation angle intrinsically more uncertain than measurements of the ring diameter.", 'ORCID iDs': 'Kazunori Akiyama https: // orcid.org / 0000-0002-9475-4254 Antxon Alberdi https: // orcid.org / 0000-0002-9371-1033 Rebecca Azulay https: // orcid.org / 0000-0002-2200-5393 Anne-Kathrin Baczko https: // orcid.org / 0000-00033090-3975  \nMislav Balokovi ć https: // orcid.org / 0000-0003-0476-6647 John Barrett https: // orcid.org / 0000-0002-9290-0764  \nLindy Blackburn https: // orcid.org / 0000-0002-9030-642X Katherine L. Bouman https: // orcid.org / 0000-00030077-4367  \nGeoffrey C. Bower https: // orcid.org / 0000-0003-4056-9982 Christiaan D. Brinkerink https: // orcid.org / 0000-00022322-0749  \nRoger Brissenden https: // orcid.org / 0000-0002-2556-0894 Silke Britzen https: // orcid.org / 0000-0001-9240-6734 Avery E. Broderick https: // orcid.org / 0000-00023351-760X  \nDo-Young Byun https: // orcid.org / 0000-0003-1157-4109 Andrew Chael https: // orcid.org / 0000-0003-2966-6220 Chi-kwan Chan https: // orcid.org / 0000-0001-6337-6126 Shami Chatterjee https: // orcid.org / 0000-0002-2878-1502  \nIlje Cho https: // orcid.org / 0000-0001-6083-7521 Pierre Christian https: // orcid.org / 0000-0001-6820-9941 John E. Conway https: // orcid.org / 0000-0003-2448-9181 Geoffrey B. Crew https: // orcid.org / 0000-0002-2079-3189 Yuzhu Cui https: // orcid.org / 0000-0001-6311-4345 Jordy Davelaar https: // orcid.org / 0000-0002-2685-2434 Mariafelicia De Laurentis https: // orcid.org / 0000-00029945-682X  \nRoger Deane https: // orcid.org / 0000-0003-1027-5043 Jessica Dempsey https: // orcid.org / 0000-0003-1269-9667 Gregory Desvignes https: // orcid.org / 0000-0003-3922-4055 Jason Dexter https: // orcid.org / 0000-0003-3903-0373 Sheperd S. Doeleman https: // orcid.org / 0000-00029031-0904  \nRalph P. Eatough https: // orcid.org / 0000-0001-6196-4135 Heino Falcke https: // orcid.org / 0000-0002-2526-6724 Vincent L. Fish https: // orcid.org / 0000-0002-7128-9345 Raquel Fraga-Encinas https: // orcid.org / 0000-00025222-1361  \nJosé L. Gómez https: // orcid.org / 0000-0003-4190-7613 Peter Galison https: // orcid.org / 0000-0002-6429-3872 Charles F. Gammie https: // orcid.org / 0000-0001-7451-8935 Boris Georgiev https: // orcid.org / 0000-0002-3586-6424 Roman Gold https: // orcid.org / 0000-0003-2492-1966 Minfeng Gu ( /uni987E /uni654F /uni5CF0 ) https: // orcid.org / 0000-0002-4455-6946 Mark Gurwell https: // orcid.org / 0000-0003-0685-3621 Kazuhiro Hada https: // orcid.org / 0000-0001-6906-772X Ronald Hesper https: // orcid.org / 0000-0003-1918-6098 Luis C. Ho ( /uni4F55 /uni5B50 /uni5C71 ) https: // orcid.org / 0000-0001-6947-5846 Mareki Honma https: // orcid.org / 0000-0003-4058-9000 Chih-Wei L. Huang https: // orcid.org / 0000-0001-5641-3953 Shiro Ikeda https: // orcid.org / 0000-0002-2462-1448 Sara Issaoun https: // orcid.org / 0000-0002-5297-921X David J. James https: // orcid.org / 0000-0001-5160-4486 Michael Janssen https: // orcid.org / 0000-0001-8685-6544 Britton Jeter https: // orcid.org / 0000-0003-2847-1712 Wu Jiang ( /uni6C5F /uni609F ) https: // orcid.org / 0000-0001-7369-3539 Michael D. Johnson https: // orcid.org / 0000-0002-4120-3029 Svetlana Jorstad https: // orcid.org / 0000-0001-6158-1708 Taehyun Jung https: // orcid.org / 0000-0001-7003-8643 Mansour Karami https: // orcid.org / 0000-0001-7387-9333 Ramesh Karuppusamy https: // orcid.org / 0000-00025307-2919  \nTomohisa Kawashima https: // orcid.org / 0000-00018527-0496  \nGarrett K. Keating https: // orcid.org / 0000-0002-3490-146X Mark Kettenis https: // orcid.org / 0000-0002-6156-5617 Jae-Young Kim https: // orcid.org / 0000-0001-8229-7183 Junhan Kim https: // orcid.org / 0000-0002-4274-9373 Motoki Kino https: // orcid.org / 0000-0002-2709-7338 Jun Yi Koay https: // orcid.org / 0000-0002-7029-6658 Patrick M. Koch https: // orcid.org / 0000-0003-2777-5861 Shoko Koyama https: // orcid.org / 0000-0002-3723-3372 Michael Kramer https: // orcid.org / 0000-0002-4175-2271 Carsten Kramer https: // orcid.org / 0000-0002-4908-4925 Thomas P. Krichbaum https: // orcid.org / 0000-00024892-9586  \nTod R. Lauer https: // orcid.org / 0000-0003-3234-7247 Sang-Sung Lee https: // orcid.org / 0000-0002-6269-594X Yan-Rong Li ( /uni674E /uni5F66 /uni8363 ) https: // orcid.org / 0000-0001-5841-9179 Zhiyuan Li ( /uni674E /uni5FD7 /uni8FDC ) https: // orcid.org / 0000-0003-0355-6437 Michael Lindqvist https: // orcid.org / 0000-0002-3669-0715  \nKuo Liu https: // orcid.org / 0000-0002-2953-7376 Elisabetta Liuzzo https: // orcid.org / 0000-0003-0995-5201 Laurent Loinard https: // orcid.org / 0000-0002-5635-3345 Ru-Sen Lu ( /uni8DEF /uni5982 /uni68EE ) https: // orcid.org / 0000-0002-7692-7967 Nicholas R. MacDonald https: // orcid.org / 0000-00026684-8691  \nJirong Mao ( /uni6BDB /uni57FA /uni8363 ) https: // orcid.org / 0000-0002-7077-7195 Sera Markoff https: // orcid.org / 0000-0001-9564-0876 Daniel P. Marrone https: // orcid.org / 0000-0002-2367-1080 Alan P. Marscher https: // orcid.org / 0000-0001-7396-3332 Iván Martí-Vidal https: // orcid.org / 0000-0003-3708-9611 Lynn D. Matthews https: // orcid.org / 0000-0002-3728-8082 Lia Medeiros https: // orcid.org / 0000-0003-2342-6728 Karl M. Menten https: // orcid.org / 0000-0001-6459-0669 Yosuke Mizuno https: // orcid.org / 0000-0002-8131-6730 Izumi Mizuno https: // orcid.org / 0000-0002-7210-6264 James M. Moran https: // orcid.org / 0000-0002-3882-4414 Kotaro Moriyama https: // orcid.org / 0000-0003-1364-3761 Monika Moscibrodzka https: // orcid.org / 0000-00024661-6332  \nCornelia Müller https: // orcid.org / 0000-0002-2739-2994 Hiroshi Nagai https: // orcid.org / 0000-0003-0292-3645 Neil M. Nagar https: // orcid.org / 0000-0001-6920-662X Masanori Nakamura https: // orcid.org / 0000-00016081-2420  \nRamesh Narayan https: // orcid.org / 0000-0002-1919-2730 Iniyan Natarajan https: // orcid.org / 0000-0001-8242-4373 Chunchong Ni https: // orcid.org / 0000-0003-1361-5699 Aristeidis Noutsos https: // orcid.org / 0000-0002-4151-3860 Héctor Olivares https: // orcid.org / 0000-0001-6833-7580 Daniel C. M. Palumbo https: // orcid.org / 0000-00027179-3816  \nUe-Li Pen https: // orcid.org / 0000-0003-2155-9578 Dominic W. Pesce https: // orcid.org / 0000-0002-5278-9221 Oliver Porth https: // orcid.org / 0000-0002-4584-2557 Ben Prather https: // orcid.org / 0000-0002-0393-7734 Jorge A. Preciado-López https: // orcid.org / 0000-00024146-0113  \nHung-Yi Pu https: // orcid.org / 0000-0001-9270-8812 Venkatessh Ramakrishnan https: // orcid.org / 0000-00029248-086X  \nRamprasad Rao https: // orcid.org / 0000-0002-1407-7944 Alexander W. Raymond https: // orcid.org / 0000-00025779-4767  \nBart Ripperda https: // orcid.org / 0000-0002-7301-3908 Freek Roelofs https: // orcid.org / 0000-0001-5461-3687 Eduardo Ros https: // orcid.org / 0000-0001-9503-4892 Mel Rose https: // orcid.org / 0000-0002-2016-8746 Alan L. Roy https: // orcid.org / 0000-0002-1931-0135 Chet Ruszczyk https: // orcid.org / 0000-0001-7278-9707 Benjamin R. Ryan https: // orcid.org / 0000-0001-8939-4461 Kazi L. J. Rygl https: // orcid.org / 0000-0003-4146-9043 David Sánchez-Arguelles https: // orcid.org / 0000-00027344-9920  \nMahito Sasada https: // orcid.org / 0000-0001-5946-9960 Tuomas Savolainen https: // orcid.org / 0000-0001-6214-1085 Lijing Shao https: // orcid.org / 0000-0002-1334-8853 Zhiqiang Shen ( /uni6C88 /uni5FD7/uni5F3A ) https: // orcid.org / 0000-00033540-8746  \nDes Small https: // orcid.org / 0000-0003-3723-5404 Bong Won Sohn https: // orcid.org / 0000-0002-4148-8378 Jason SooHoo https: // orcid.org / 0000-0003-1938-0720  \nFumie Tazaki https: // orcid.org / 0000-0003-0236-0600 Paul Tiede https: // orcid.org / 0000-0003-3826-5648 Remo P. J. Tilanus https: // orcid.org / 0000-0002-6514-553X Michael Titus https: // orcid.org / 0000-0002-3423-4505 Kenji Toma https: // orcid.org / 0000-0002-7114-6010 Pablo Torne https: // orcid.org / 0000-0001-8700-6058 Sascha Trippe https: // orcid.org / 0000-0003-0465-1559 Ilse van Bemmel https: // orcid.org / 0000-0001-5473-2950 Huib Jan van Langevelde https: // orcid.org / 0000-00020230-5946 Daniel R. van Rossum https: // orcid.org / 0000-0001-7772-6131 John Wardle https: // orcid.org / 0000-0002-8960-2942 Jonathan Weintroub https: // orcid.org / 0000-00024603-5204 Norbert Wex https: // orcid.org / 0000-0003-4058-2837 Robert Wharton https: // orcid.org / 0000-0002-7416-5209 Maciek Wielgus https: // orcid.org / 0000-0002-8635-4242 George N. Wong https: // orcid.org / 0000-0001-6952-2147 Qingwen Wu ( /uni5434 /uni5E86 /uni6587 ) https: // orcid.org / 0000-0003-4773-4987 André Young https: // orcid.org / 0000-0003-0000-2682 Ken Young https: // orcid.org / 0000-0002-3666-4920 Ziri Younsi https: // orcid.org / 0000-0001-9283-1191 Feng Yuan ( /uni8881 /uni5CF0 ) https: // orcid.org / 0000-0003-3564-6437 J. Anton Zensus https: // orcid.org / 0000-0001-7470-3321 Guangyao Zhao https: // orcid.org / 0000-0002-4417-1659 Shan-Shan Zhao https: // orcid.org / 0000-0002-9774-3606 Joseph R. Farah https: // orcid.org / 0000-0003-4914-5625 Daniel Michalik https: // orcid.org / 0000-0002-7618-6556 Andrew Nadolski https: // orcid.org / 0000-0001-9479-9957 Rurik A. Primiani https: // orcid.org / 0000-0003-3910-7529 Paul Yamaguchi https: // orcid.org / 0000-0002-6017-8199', 'References': "- Abramowitz, M., & Stegun, I. A. 1972, Handbook of Mathematical Functions\n- Akiyama, K., Ikeda, S., Pleau, M., et al. 2017a, AJ, 153, 159\n- Akiyama, K., Kuramochi, K., Ikeda, S., et al. 2017b, ApJ, 838, 1\n- Akiyama, K., Lu, R.-S., Fish, V. L., et al. 2015, ApJ, 807, 150\n- Akiyama, K., Tazaki, F., Moriyama, K., et al. 2019, SMILI: Sparse Modeling Imaging Library for Interferometry, Zenodo, doi:10.5281 / zenodo.2616725 An, T., Sohn, B. W., & Imai, H. 2018, NatAs, 2, 118\n- Asada, K., & Nakamura, M. 2012, ApJL, 745, L28\n- Asada, K., Nakamura, M., Doi, A., Nagai, H., & Inoue, M. 2014, ApJL, 781, L2\n- Asada, K., Nakamura, M., & Pu, H.-Y. 2016, ApJ, 833, 56\n- Bardeen, J. M. 1973, in Black Holes ( Les Astres Occlus ) , ed. C. DeWitt & B. S. DeWitt ( New York: Gordon and Breach ) , 215  \nBaron, F., Monnier, J. D., & Kloppenborg, B. 2010, Proc. SPIE, 7734, 77342I  \nBird, S., Harris, W. E., Blakeslee, J. P., & Flynn, C. 2010, A&A, 524, A71  \nBiretta, J. A., Sparks, W. B., & Macchetto, F. 1999, ApJ, 520, 621  \nBlackburn, L., Chan, C-K., Crew, G. B., et al. 2019, arXiv:1903.08832  \n- Blakeslee, J. P., Jordán, A., Mei, S., et al. 2009, ApJ, 694, 556\n- Blecher, T., Deane, R., Bernardi, G., & Smirnov, O. 2017, MNRAS, 464, 143\n- Bouman, K. L. 2017, PhD thesis, Massachusetts Institute of Technology\n- Bouman, K. L., Johnson, M. D., Dalca, A. V., et al. 2018, IEEE Transactions on Computational Imaging, 4, 512\n- Bouman, K. L., Johnson, M. D., Zoran, D., et al. 2016, in The IEEE Conf. Computer Vision and Pattern Recognition ( New York: IEEE ) , 913 Briggs, D. S. 1995, BAAS, 27, 1444\n- Britzen, S., Fendt, C., Eckart, A., & Karas, V. 2017, A&A, 601, A52\n- Buscher, D. F. 1994, in IAU Symp. 158, Very High Angular Resolution Imaging, ed. J. G. Robertson & W. J. Tango ( Dordrecht: Kluwer ) , 91  \nCantiello, M., Blakeslee, J. P., Ferrarese, L., et al. 2018, ApJ, 856, 126  \n- Chael, A., Bouman, K., Johnson, M., et al. 2019, eht-imaging: v1.1.0: Imaging interferometric data with regularized maximum likelihood, Zenodo, doi:10. 5281 / zenodo.2614016\n- Chael, A. A., Johnson, M. D., Bouman, K. L., et al. 2018, ApJ, 857, 23\n- Chael, A. A., Johnson, M. D., Narayan, R., et al. 2016, ApJ, 829, 11  \n- Shepherd, M. 2011, Difmap: Synthesis Imaging of Visibility Data, Astrophysics Source Code Library, ascl:1103.001\n- Shepherd, M. C. 1997, in ASP Conf. Ser. 125, Astronomical Data Analysis Software and Systems VI, ed. G. Hunt & H. Payne ( San Francisco, CA: ASP ) , 77\n- Sparks, W. B., Biretta, J. A., & Macchetto, F. 1996, ApJ, 473, 254\n- Thiébaut, É. 2013, in New Concepts in Imaging: Optical and Statistical Models, EAS Publications Ser. /uniF0A0 59, ed. D. Mary, C. Theys, & C. Aime ( Les Ulis: EDP Sciences ) , 157\n- Thompson, A. R., Moran, J. M., & Swenson, G. W., Jr. 2017, Interferometry and Synthesis in Radio Astronomy ( 3rd ed.; Hoboken, NJ: Wiley )\n- Twiss, R. Q., Carter, A. W. L., & Little, A. G. 1960, Obs, 80, 153 van Cittert, P. H. 1934, Phy, 1, 201  \n- Wajima, K., Hagiwara, Y., An, T., et al. 2016, in Proc. Workshop Frontiers in Radio Astronomy and FAST Early Sciences Symposium 2015, ed.\n- Walker, R. C., Hardee, P. E., Davies, F. B., Ly, C., & Junor, W. 2018, ApJ,\n- Walsh, J. L., Barth, A. J., Ho, L. C., & Sarzi, M. 2013, ApJ, 770, 86\n- Wiaux, Y., Jacques, L., Puy, G., Scaife, A. M. M., & Vandergheynst, P. 2009a,\n- Wiaux, Y., Puy, G., Boursier, Y., & Vandergheynst, P. 2009b, MNRAS, 400, 1029 Wilkinson, P. N., Readhead, A. C. S., Purcell, G. H., & Anderson, B. 1977,  \nL. Qain & D. Li ( San Francisco, CA: ASP ) , 81 855, 128 MNRAS, 395, 1733 Natur, 269, 764 Thiébaut, É., & Young, J. 2017, JOSAA, 34, 904 Zernike, F. 1938, Phy, 785, 5  \nThe Event Horizon Telescope Collaboration  \nKazunori Akiyama 1,2,3,4 , Antxon Alberdi 5 , Walter Alef 6 , Keiichi Asada 7 , Rebecca Azulay 8,9,6 , Anne-Kathrin Baczko 6 , David Ball 10 , Mislav Balokovi ć 4,11 , John Barrett 2 , Dan Bintley 12 , Lindy Blackburn 4,11 , Wilfred Boland 13 , Katherine L. Bouman 4,11,14 , Geoffrey C. Bower 15 , Michael Bremer 16 , Christiaan D. Brinkerink 17 , Roger Brissenden 4,11 , Silke Britzen 6 , Avery E. Broderick 18,19,20 , Dominique Broguiere 16 , Thomas Bronzwaer 17 , Do-Young Byun 21,22 , John E. Carlstrom 23,24,25,26 , Andrew Chael 4,11 , Chi-kwan Chan 10,27 , Shami Chatterjee 28 , Koushik Chatterjee 29 , Ming-Tang Chen 15 , Yongjun Chen ( /uni9648 /uni6C38 /uni519B ) 30,31 , Ilje Cho 21,22 , Pierre Christian 10,11 , John E. Conway 32 , James M. Cordes 28 , Geoffrey B. Crew 2 , Yuzhu Cui 33,34 , Jordy Davelaar 17 , Mariafelicia De Laurentis 35,36,37 , Roger Deane 38,39 , Jessica Dempsey 12 , Gregory Desvignes 6 , Jason Dexter 40 , Sheperd S. Doeleman 4,11 , Ralph P. Eatough 6 , Heino Falcke 17 , Vincent L. Fish 2 , Ed Fomalont 1 , Raquel Fraga-Encinas 17 , William T. Freeman 41,42 , Per Friberg 12 , Christian M. Fromm 36 , José L. Gómez 5 , Peter Galison 4,43,44 , Charles F. Gammie 45,46 , Roberto García 16 , Olivier Gentaz 16 , Boris Georgiev 19,20 , Ciriaco Goddi 17,47 , Roman Gold 36 , Minfeng Gu ( /uni987E /uni654F /uni5CF0 ) 30,48 , Mark Gurwell 11 , Kazuhiro Hada 33,34 , Michael H. Hecht 2 , Ronald Hesper 49 , Luis C. Ho ( /uni4F55 /uni5B50 /uni5C71 ) 50,51 , Paul Ho 7 , Mareki Honma 33,34 , Chih-Wei L. Huang 7 , Lei Huang ( /uni9EC4 /uni78CA ) 30,48 , David H. Hughes 52 , Shiro Ikeda 3,53,54,55 , Makoto Inoue 7 , Sara Issaoun 17 , David J. James 4,11 , Buell T. Jannuzi 10 , Michael Janssen 17 , Britton Jeter 19,20 , Wu Jiang ( /uni6C5F /uni609F ) 30 , Michael D. Johnson 4,11 , Svetlana Jorstad 56,57 , Taehyun Jung 21,22 , Mansour Karami 18,19 , Ramesh Karuppusamy 6 , Tomohisa Kawashima 3 , Garrett K. Keating 11 , Mark Kettenis 58 , Jae-Young Kim 6 , Junhan Kim 10 , Jongsoo Kim 21 , Motoki Kino 3,59 , Jun Yi Koay 7 , Patrick M. Koch 7 , Shoko Koyama 7 , Michael Kramer 6 , Carsten Kramer 16 , Thomas P. Krichbaum 6 , Cheng-Yu Kuo 60 , Tod R. Lauer 61 , Sang-Sung Lee 21 , Yan-Rong Li ( /uni674E /uni5F66 /uni8363 ) 62 , Zhiyuan Li ( /uni674E /uni5FD7 /uni8FDC ) 63,64 , Michael Lindqvist 32 , Kuo Liu 6 , Elisabetta Liuzzo 65 , Wen-Ping Lo 7,66 , Andrei P. Lobanov 6 , Laurent Loinard 67,68 , Colin Lonsdale 2 , Ru-Sen Lu ( /uni8DEF /uni5982 /uni68EE ) 30,6 , Nicholas R. MacDonald 6 , Jirong Mao ( /uni6BDB /uni57FA /uni8363 ) 69,70,71 , Sera Markoff 29,72 , Daniel P. Marrone 10 , Alan P. Marscher 56 , Iván Martí-Vidal 32,73 , Satoki Matsushita 7 , Lynn D. Matthews 2 , Lia Medeiros 10,74 , Karl M. Menten 6 , Yosuke Mizuno 36 , Izumi Mizuno 12 , James M. Moran 4,11 , Kotaro Moriyama 33,2 , Monika Moscibrodzka 17 , Cornelia Müller 6,17 , Hiroshi Nagai 3,34 , Neil M. Nagar 75 , Masanori Nakamura 7 , Ramesh Narayan 4,11 , Gopal Narayanan 76 , Iniyan Natarajan 39 , Roberto Neri 16 , Chunchong Ni 19,20 , Aristeidis Noutsos 6 , Hiroki Okino 33,77 , Héctor Olivares 36 , Tomoaki Oyama 33 , Feryal Özel 10 , Daniel C. M. Palumbo 4,11 , Nimesh Patel 11 , Ue-Li Pen 18,78,79,80 , Dominic W. Pesce 4,11 , Vincent Piétu 16 , Richard Plambeck 81 , Aleksandar PopStefanija 76 , Oliver Porth 36,29 , Ben Prather 45 , Jorge A. Preciado-López 18 , Dimitrios Psaltis 10 , Hung-Yi Pu 18 , Venkatessh Ramakrishnan 75 , Ramprasad Rao 15 , Mark G. Rawlings 12 , Alexander W. Raymond 4,11 , Luciano Rezzolla 36 , Bart Ripperda 36 , Freek Roelofs 17 , Alan Rogers 2 , Eduardo Ros 6 , Mel Rose 10 , Arash Roshanineshat 10 , Helge Rottmann 6 , Alan L. Roy 6 , Chet Ruszczyk 2 , Benjamin R. Ryan 82,83 , Kazi L. J. Rygl 65 , Salvador Sánchez 84 , David Sánchez-Arguelles 52,85 , Mahito Sasada 33,86 , Tuomas Savolainen 6,87,88 , F. Peter Schloerb 76 , Karl-Friedrich Schuster 16 , Lijing Shao 6,51 , Zhiqiang Shen ( /uni6C88 /uni5FD7/uni5F3A ) 29,30 , Des Small 58 , Bong Won Sohn 21,22,89 , Jason SooHoo 2 , Fumie Tazaki 33 , Paul Tiede 18,19,20 , Remo P. J. Tilanus 17,47,90 , Michael Titus 2 , Kenji Toma 91,92 , Pablo Torne 6,84 , Tyler Trent 10 , Sascha Trippe 93 , Shuichiro Tsuda 33 , Ilse van Bemmel 58 , Huib Jan van Langevelde 58,94 , Daniel R. van Rossum 17 , Jan Wagner 6 , John Wardle 95 , Jonathan Weintroub 4,11 , Norbert Wex 6 , Robert Wharton 6 , Maciek Wielgus 4,11 , George N. Wong 45 , Qingwen Wu ( /uni5434 /uni5E86 /uni6587 ) 96 , André Young 17 , Ken Young 11 , Ziri Younsi 97,36 , Feng Yuan ( /uni8881 /uni5CF0 ) 30,48,98 , Ye-Fei Yuan ( /uni8881 /uni4E1A /uni98DE ) 99 , J. Anton Zensus 6 , Guangyao Zhao 21 , Shan-Shan Zhao 17,63 , Ziyan Zhu 44 , Joseph R. Farah 11,100,4 , Zheng Meyer-Zhao 7,101 , Daniel Michalik 102,103 , Andrew Nadolski 46 , Hiroaki Nishioka 7 , Nicolas Pradel 7 , Rurik A. Primiani 104 , Kamal Souccar 76 , Laura Vertatschitsch 11,104 , and Paul Yamaguchi 11  \n- 1 National Radio Astronomy Observatory, 520 Edgemont Rd, Charlottesville, VA 22903, USA\n- 2 Massachusetts Institute of Technology Haystack Observatory, 99 Millstone Road, Westford, MA 01886, USA\n- 3 National Astronomical Observatory of Japan, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan\n- 4 Black Hole Initiative at Harvard University, 20 Garden Street, Cambridge, MA 02138, USA  \nInstituto de Astrofísica de Andalucía  \n-  \nCSIC, Glorieta de la Astronomía s  \n6  \n5  \n/  \nn, E-18008 Granada, Spain  \n- Max-Planck-Institut für Radioastronomie, Auf dem Hügel 69, D-53121 Bonn, Germany 7 Institute of Astronomy and Astrophysics, Academia Sinica, 11F of Astronomy-Mathematics Building, AS / NTU No. 1, Sec. 4, Roosevelt Rd, Taipei 10617, Taiwan, R.O.C.  \nDepartament d  \n'  \nAstronomia i Astrofísica, Universitat de València, C. Dr. Moliner 50, E-46100 Burjassot, València, Spain  \n9  \nObservatori Astronòmic, Universitat de València, C. Catedrático José Beltrán 2, E-46980 Paterna, València, Spain  \nSteward Observatory and Department of Astronomy, University of Arizona, 933 N. Cherry Ave., Tucson, AZ 85721, USA  \n- 11 Center for Astrophysics | Harvard & Smithsonian, 60 Garden Street, Cambridge, MA 02138, USA 12 East Asian Observatory, 660 N. A ' ohoku Pl., Hilo, HI 96720, USA  \nNederlandse Onderzoekschool voor Astronomie  \n14  \n(  \nNOVA  \n)  \n, PO Box 9513, 2300 RA Leiden, The Netherlands  \nCalifornia Institute of Technology, 1200 East California Boulevard, Pasadena, CA 91125, USA  \nInstitute of Astronomy and Astrophysics, Academia Sinica, 645 N. A  \n'  \nohoku Place, Hilo, HI 96720, USA  \n13  \n15  \n16  \nInstitut de Radioastronomie Millimétrique, 300 rue de la Piscine, F-38406 Saint Martin d  \n'  \nHères, France  \n17  \nDepartment of Astrophysics, Institute for Mathematics, Astrophysics and Particle Physics  \n(  \nIMAPP  \n)  \n,  \nRadboud University, P.O. Box 9010, 6500 GL Nijmegen, The Netherlands  \nPerimeter Institute for Theoretical Physics, 31 Caroline Street North, Waterloo, ON, N2L 2Y5, Canada  \n18  \nDepartment of Physics and Astronomy, University of Waterloo, 200 University Avenue West, Waterloo, ON, N2L 3G1, Canada  \n20  \nWaterloo Centre for Astrophysics, University of Waterloo, Waterloo, ON N2L 3G1 Canada  \nKorea Astronomy and Space Science Institute, Daedeok-daero 776, Yuseong-gu, Daejeon 34055, Republic of Korea  \n22  \nUniversity of Science and Technology, Gajeong-ro 217, Yuseong-gu, Daejeon 34113, Republic of Korea  \nKavli Institute for Cosmological Physics, University of Chicago, 5640 South Ellis Avenue, Chicago, IL 60637, USA  \n21  \n23  \n24  \nDepartment of Astronomy and Astrophysics, University of Chicago, 5640 South Ellis Avenue, Chicago, IL 60637, USA  \n- 25 Department of Physics, University of Chicago, 5720 South Ellis Avenue, Chicago, IL 60637, USA\n- 26 Enrico Fermi Institute, University of Chicago, 5640 South Ellis Avenue, Chicago, IL 60637, USA 27 Data Science Institute, University of Arizona, 1230 N. Cherry Ave., Tucson, AZ 85721, USA\n- 28 Cornell Center for Astrophysics and Planetary Science, Cornell University, Ithaca, NY 14853, USA\n- 29 Anton Pannekoek Institute for Astronomy, University of Amsterdam, Science Park 904, 1098 XH, Amsterdam, The Netherlands  \nShanghai Astronomical Observatory, Chinese Academy of Sciences, 80 Nandan Road, Shanghai 200030, People  \n'  \ns Republic of China  \n31  \nKey Laboratory of Radio Astronomy, Chinese Academy of Sciences, Nanjing 210008, People  \n'  \ns Republic of China  \nDepartment of Space, Earth and Environment, Chalmers University of Technology, Onsala Space Observatory, SE-439 92 Onsala, Sweden  \n32  \n33  \nMizusawa VLBI Observatory, National Astronomical Observatory of Japan, 2-12 Hoshigaoka, Mizusawa, Oshu, Iwate 023-0861, Japan  \n(  \n)  \n- INFN Sez. di Napoli, Compl. Univ. di Monte S. Angelo, Edi /uniFB01 cio G, Via Cinthia, I-80126, Napoli, Italy\n- Department of Astronomical Science, The Graduate University for Advanced Studies SOKENDAI , 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan 35 Dipartimento di Fisica ' E. Pancini ' , Universitá di Napoli ' Federico II ' , Compl. Univ. di Monte S. Angelo, Edi /uniFB01 cio G, Via Cinthia, I-80126, Napoli, Italy 36 Institut für Theoretische Physik, Goethe-Universität Frankfurt, Max-von-Laue-Straße 1, D-60438 Frankfurt am Main, Germany 37\n- 38 Department of Physics, University of Pretoria, Lynnwood Road, Hat /uniFB01 eld, Pretoria 0083, South Africa\n- 39 Centre for Radio Astronomy Techniques and Technologies, Department of Physics and Electronics, Rhodes University, Grahamstown 6140, South Africa 40 Max-Planck-Institut für Extraterrestrische Physik, Giessenbachstr. 1, D-85748 Garching, Germany\n- 41 Department of Electrical Engineering and Computer Science, Massachusetts Institute of Technology,\n- 32-D476, 77 Massachussetts Ave., Cambridge, MA 02142, USA\n- 42 Google Research, 355 Main St., Cambridge, MA 02142, USA\n- 43 Department of History of Science, Harvard University, Cambridge, MA 02138, USA  \n44  \nDepartment of Physics, Harvard University, Cambridge, MA 02138, USA  \n- 45 Department of Physics, University of Illinois, 1110 West Green St, Urbana, IL 61801, USA  \nDepartment of Astronomy, University of Illinois at Urbana-Champaign, 1002 West Green Street, Urbana, IL 61801, USA  \n- 47 Leiden Observatory -Allegro, Leiden University, P.O. Box 9513, 2300 RA Leiden, The Netherlands  \nKey Laboratory for Research in Galaxies and Cosmology, Chinese Academy of Sciences, Shanghai 200030, People  \n'  \ns Republic of China  \n48  \nNOVA Sub-mm Instrumentation Group, Kapteyn Astronomical Institute, University of Groningen, Landleven 12, 9747 AD Groningen, The Netherlands  \n50  \nDepartment of Astronomy, School of Physics, Peking University, Beijing 100871, People  \n'  \ns Republic of China  \n- 51 Kavli Institute for Astronomy and Astrophysics, Peking University, Beijing 100871, People ' s Republic of China 52  \nInstituto Nacional de Astrofísica, Óptica y Electrónica. Apartado Postal 51 y 216, 72000. Puebla Pue., México  \n53  \nThe Institute of Statistical Mathematics, 10-3 Midori-cho, Tachikawa, Tokyo, 190-8562, Japan  \nDepartment of Statistical Science, The Graduate University for Advanced Studies  \n55  \n(  \nSOKENDAI  \n)  \n, 10-3 Midori-cho, Tachikawa, Tokyo 190-8562, Japan  \nKavli Institute for the Physics and Mathematics of the Universe, The University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa, 277-8583, Japan  \n56  \nInstitute for Astrophysical Research, Boston University, 725 Commonwealth Ave., Boston, MA 02215, USA  \nAstronomical Institute, St. Petersburg University, Universitetskij pr., 28, Petrodvorets,198504 St. Petersburg, Russia  \n58  \nJoint Institute for VLBI ERIC  \n(  \nJIVE  \n)  \n, Oude Hoogeveensedijk 4, 7991 PD Dwingeloo, The Netherlands  \nKogakuin University of Technology & Engineering, Academic Support Center, 2665-1 Nakano, Hachioji, Tokyo 192-0015, Japan  \n60  \nPhysics Department, National Sun Yat-Sen University, No. 70, Lien-Hai Rd, Kaosiung City 80424, Taiwan, R.O.C  \n61  \nNational Optical Astronomy Observatory, 950 North Cherry Ave., Tucson, AZ 85719, USA  \n- 62 Key Laboratory for Particle Astrophysics, Institute of High Energy Physics, Chinese Academy of Sciences,  \n19B Yuquan Road, Shijingshan District, Beijing, People ' s Republic of China  \n- 63 School of Astronomy and Space Science, Nanjing University, Nanjing 210023, People ' s Republic of China  \nKey Laboratory of Modern Astronomy and Astrophysics, Nanjing University, Nanjing 210023, People  \n'  \ns Republic of China  \n- 65 Italian ALMA Regional Centre, INAF-Istituto di Radioastronomia, Via P. Gobetti 101, 40129 Bologna, Italy\n- 66 Department of Physics, National Taiwan University, No.1, Sect.4, Roosevelt Rd., Taipei 10617, Taiwan, R.O.C\n- 67 Instituto de Radioastronomía y Astrofísica, Universidad Nacional Autónoma de México, Morelia 58089, México 68 Instituto de Astronomía, Universidad Nacional Autónoma de México, CdMx 04510, México\n- 69 Yunnan Observatories, Chinese Academy of Sciences, 650011 Kunming, Yunnan Province, People ' s Republic of China  \nCenter for Astronomical Mega-Science, Chinese Academy of Sciences, 20A Datun Road, Chaoyang District, Beijing, 100012, People  \n'  \ns Republic of China  \n71  \nKey Laboratory for the Structure and Evolution of Celestial Objects, Chinese Academy of Sciences, 650011 Kunming, People  \n'  \ns Republic of China  \n34  \n49  \n54  \n70  \n30  \n19  \n59  \n8  \n10  \n46  \n64  \n57  \n- 72 Gravitation Astroparticle Physics Amsterdam ( GRAPPA ) Institute, University of Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands 73 Centro Astronómico de Yebes ( IGN ) , Apartado 148, E-19180 Yebes, Spain 74  \nDepartment of Physics, Broida Hall, University of California Santa Barbara, Santa Barbara, CA 93106, USA  \n- 75 Astronomy Department, Universidad de Concepción, Casilla 160-C, Concepción, Chile\n- 76 Department of Astronomy, University of Massachusetts, 01003, Amherst, MA, USA  \nDepartment of Astronomy, Graduate School of Science, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan  \n78  \n79  \nCanadian Institute for Theoretical Astrophysics, University of Toronto, 60 St. George Street, Toronto, ON M5S 3H8, Canada  \nDunlap Institute for Astronomy and Astrophysics, University of Toronto, 50 St. George Street, Toronto, ON M5S 3H4, Canada  \n80  \nCanadian Institute for Advanced Research, 180 Dundas St West, Toronto, ON M5G 1Z8, Canada  \n81  \nRadio Astronomy Laboratory, University of California, Berkeley, CA 94720, USA  \nCCS-2, Los Alamos National Laboratory, P.O. Box 1663, Los Alamos, NM 87545, USA  \n82  \n- 83 Center for Theoretical Astrophysics, Los Alamos National Laboratory, Los Alamos, NM, 87545, USA  \nInstituto de Radioastronomía Milimétrica, IRAM, Avenida Divina Pastora 7, Local 20, E-18012, Granada, Spain  \n85  \nConsejo Nacional de Ciencia y Tecnología, Av. Insurgentes Sur 1582, 03940, Ciudad de México, México  \nHiroshima Astrophysical Science Center, Hiroshima University, 1-3-1 Kagamiyama, Higashi-Hiroshima, Hiroshima 739-8526, Japan  \n87  \nAalto University Department of Electronics and Nanoengineering, PL 15500, FI-00076 Aalto, Finland  \n88  \nAalto University Metsähovi Radio Observatory, Metsähovintie 114, FI-02540 Kylmälä, Finland  \nDepartment of Astronomy, Yonsei University, Yonsei-ro 50, Seodaemun-gu, 03722 Seoul, Republic of Korea  \n89  \n90  \n93  \nNetherlands Organisation for Scienti  \n91  \n/uniFB01  \nc Research  \n(  \nNWO  \n)  \n, Postbus 93138, 2509 AC Den Haag, The Netherlands  \nFrontier Research Institute for Interdisciplinary Sciences, Tohoku University, Sendai 980-8578, Japan  \n92 Astronomical Institute, Tohoku University, Sendai 980-8578, Japan  \nDepartment of Physics and Astronomy, Seoul National University, Gwanak-gu, Seoul 08826, Republic of Korea  \n94  \nLeiden Observatory, Leiden University, Postbus 2300, 9513 RA Leiden, The Netherlands  \n95  \nPhysics Department, Brandeis University, 415 South Street, Waltham, MA 02453, USA  \n- 96 School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, 430074, People ' s Republic of China  \n97  \nMullard Space Science Laboratory, University College London, Holmbury St. Mary, Dorking, Surrey, RH5 6NT, UK  \nSchool of Astronomy and Space Sciences, University of Chinese Academy of Sciences, No. 19A Yuquan Road, Beijing 100049, People  \n'  \ns Republic of China  \n99  \n102  \nAstronomy Department, University of Science and Technology of China, Hefei 230026, People  \n'  \ns Republic of China  \n100  \nUniversity of Massachusetts Boston, 100 William T, Morrissey Blvd, Boston, MA 02125, USA  \n101  \nASTRON, Oude Hoogeveensedijk 4, 7991 PD Dwingeloo, The Netherlands  \nScience Support Of  \n/uniFB01  \nce, Directorate of Science, European Space Research and Technology Centre  \nKeplerlaan 1, 2201 AZ Noordwijk, The Netherlands  \nUniversity of Chicago, 5640 South Ellis Avenue, Chicago, IL 60637, USA  \n103  \n104 Systems and Technology Research, 600 West Cummings Park, Woburn, MA 01801, USA  \n98  \n77  \n86  \n84  \n(  \nESA  \nESTEC  \n/  \n)  \n,"}