Accueil > Physique > Trous Noirs :l’image de M87 > M87-2 Array et Instrumentation
M87-2 Array et Instrumentation
Publié le 15 avril 2019, dernière mise-à-jour le 1er mai 2019, 1 visites, visites totales.
II. Array and Instrumentation
The Event Horizon Telescope (EHT) is a very long baseline interferometry (VLBI) array that comprises millimeter- and submillimeter-wavelength telescopes separated by distances comparable to the diameter of the Earth. At a nominal operating wavelength of 1.3 mm, EHT angular resolution (λ/D) is 25 μas, which is sufficient to resolve nearby supermassive black hole candidates on spatial and temporal scales that correspond to their event horizons. With this capability, the EHT scientific goals are to probe general relativistic effects in the strong-field regime and to study accretion and relativistic jet formation near the black hole boundary. In this Letter we describe the system design of the EHT, detail the technology and instrumentation that enable observations, and provide measures of its performance. Meeting the EHT science objectives has required several key developments that have facilitated the robust extension of the VLBI technique to EHT observing wavelengths and the production of instrumentation that can be deployed on a heterogeneous array of existing telescopes and facilities. To meet sensitivity requirements, high-bandwidth digital systems were developed that process data at rates of 64 gigabit s−1, exceeding those of currently operating cm-wavelength VLBI arrays by more than an order of magnitude. Associated improvements include the development of phasing systems at array facilities, new receiver installation at several sites, and the deployment of hydrogen maser frequency standards to ensure coherent data capture across the array. These efforts led to the coordination and execution of the first Global EHT observations in 2017 April, and to event-horizon-scale imaging of the supermassive black hole candidate in M87.
1. Introduction
It is generally accepted that active galactic nuclei (AGNs) are powered by accretion onto supermassive black holes (SMBHs ; Heckman & Best 2014). These central engines are powerful actors on the cosmic stage, with roles in galactic evolution, star formation, mergers, and particle acceleration as evidenced by relativistic jets that both dynamically influence and redistribute matter on galactic scales (Blandford et al. 2018). Inflowing material typically obscures the event horizons of these black hole candidates, but it is in this extreme environment of the black hole boundary that strong-field effects of general relativity become evident and the accretion and outflow processes that govern black hole feedback on galactic scales originate (Ho 2008). Imaging black holes on scales that resolve these effects and processes would enable new tests of general relativity and the extraordinarily detailed study of core AGN physics. Realization of this goal requires a specialized instrument that does two things. It must have the ultra-high angular resolution required to resolve the nearest SMBH candidates, and it must operate in a range of the electromagnetic spectrum where light streams unimpeded from the innermost accretion region to telescopes on Earth. Achieving these specifications is the primary objective of the Event Horizon Telescope (EHT) : a very long baseline interferometry (VLBI) array of millimeter (mm) and submillimeter (submm) wavelength facilities that span the globe, creating a telescope with an effective Earth-sized aperture (Doeleman et al. 2009a).
While the EHT is uniquely designed for the imaging of SMBHs, other pioneering instruments are capable of probing similar angular scales for other purposes. Operating in the infrared, the GRAVITY interferometer delivers relative astrometry at the 10 micro-arcsecond (μas) level and has provided evidence for relativistic motion of material in close proximity to Sgr A* (GRAVITY Collaboration et al. 2018a). These infrared observations are an important and parallel probe of the spacetime surrounding Sgr A*, but cannot be used to make spatially resolved images of the black hole candidate because the interferometer intrinsic resolution is only 3 mas. The RadioAstron satellite, used as an orbiting element of combined Earth-Space VLBI arrays, is also capable of 10 μas angular resolution (Kardashev et al. 2013), but it operates at longer radio wavelengths that cannot penetrate the self-absorbed synchrotron plasma that surrounds the event horizon. In many ways the EHT is complementary to the Laser Interferometer Gravitational-Wave Observatory (LIGO) facility, which has detected the gravitational wave signatures from merging stellar-mass black holes (Abbott et al. 2016). LIGO and the EHT observe black holes that differ in mass by factors of $10^4\mbox—10^7$ ; LIGO events are transient, while the EHT carries out long-term studies of its main targets.
This Letter is one in a sequence of manuscripts that describes the first EHT results. The full sequence includes an abstract Letter with a summary of results (EHT Collaboration et al. 2019a, hereafter Paper I), this Letter with a description of the array and instrumentation (Paper II), a description of the data pipeline and processing (EHT Collaboration et al. 2019b, hereafter Paper III), a description of imaging techniques (EHT Collaboration et al. 2019c, hereafter Paper IV), and theoretical analyses of astrophysical and physics results (EHT Collaboration et al. 2019d, hereafter Paper V ; EHT Collaboration et al. 2019e, hereafter Paper VI, respectively).
1.1. EHT Science Goals
The scientific goals of the EHT array are to resolve and detect strong general relativistic signatures expected to arise on event-horizon scales. The best-known effect is that a black hole, surrounded by an optically thin luminous plasma, should exhibit a ’’silhouette’’ or ’’shadow’’ morphology : a dim central region delineated by the lensed photon orbit (Falcke et al. 2000). The apparent size of the photon orbit, described not long after Schwarzschild’s initial solution was published (Hilbert 1917 ; von Laue 1921), defines a bright ring or crescent shape that was calculated for arbitrary spin by Bardeen (1973), first imaged through simulations by Luminet 1979, and subsequently studied extensively (Chandrasekhar 1983 ; Takahashi 2004 ; Broderick & Loeb 2006). The size and shape of the resulting shadow depends primarily on the mass of the black hole, and only very weakly on its spin and the observing orientation. For a non-spinning black hole, the shadow diameter is equal to $\sqrt27$ Schwarzschild radii (Rs = 2GM/c2). Over all black hole spins and orientations, the shadow diameter ranges from 4.8 to 5.2 Rs (Bardeen 1973 ; Johannsen & Psaltis 2010). Because any light that crosses the photon orbit from outside will eventually reach the event horizon, use of the term ’’horizon scale’’ will hereafter be understood to mean the size of the shadow and the lensed photon orbit. Detection of the shadow via lensed electromagnetic radiation would provide new evidence for the existence of SMBHs by confining the masses of EHT targets to within their expected photon orbits. A more detailed study of the precise shape of the photon orbit can be used to test the validity of general relativity on horizon scales (Johannsen & Psaltis 2010). Full polarimetric imaging can similarly be used to map magnetic field structure near the event horizon, placing important constraints on modes of accretion and the launching of relativistic jets (Broderick & Loeb 2009 ; Johnson et al. 2015 ; Chael et al. 2016 ; Akiyama et al. 2017 ; Gold et al. 2017).
Separate signatures, potentially offering a more sensitive probe of black hole spin, are the timescales of dynamical processes at horizon scales. A characteristic timescale for such processes is given by the orbital period of test particles at the innermost stable circular orbit (ISCO), which depends sensitively on the black hole spin (Bardeen et al. 1972). Monitoring VLBI observables to track the orbital dynamics of inhomogeneities in the accretion flow can thus be used to probe the spacetime, and potentially spin, of the black hole (Broderick & Loeb 2006 ; Doeleman et al. 2009b ; Fish et al. 2009 ; Fraga-Encinas et al. 2016 ; Medeiros et al. 2017 ; Roelofs et al. 2017).
1.2. Target Sources and Confirmation of Horizon-scale Structure
The bright radio core of M87 and Sgr A* (Table 1) are the primary EHT targets, as the combination of their estimated mass and proximity make them the two most suitable sources for studying SMBH candidates at horizon-scale resolution (Narayan & McClintock 2008 ; Johannsen et al. 2012). With bolometric luminosities well below the Eddington limit (Di Matteo et al. 2000 ; Baganoff et al. 2003), both the nucleus of M87 and Sgr A* are representative of the broad and populous class of low-luminosity AGN (LLAGN). AGN spend most of their time in this low state after prior periods of high accretion (Ho 2008), and LLAGN share many characteristics of stellar black hole emission in X-ray binaries, which exhibit episodic emission and jet production (Narayan & McClintock 2008). Thus, the prospects for applying results from horizon-scale observations of the central region of M87 and Sgr A* across a broad range of astrophysical contexts are excellent.
table1
aSgr A* : $\alpha _\rmJ2000.0=17^\rmh
45^\rmm
40\buildrel\rms\over. 0409,\,\delta _\rmJ2000.0=-29^\circ 00^\prime 28\buildrel\prime\prime\over. 118$ (10) ; M87 : $\alpha _\rmJ2000.0=12^\rmh
30^\rmm
49\buildrel\rms\over. 4234,\,\delta _\rmJ2000.0=12^\circ 23^\prime 28\buildrel\prime\prime\over. 044$ (11). bThe shadow diameter is within the range 4.8–5.2 Rs depending on black hole spin and orientation to the observer’s line of sight (Johannsen & Psaltis 2010). cBrightness temperatures are reported for an observing frequency of 230 GHz. d $P_\mathrmISCO$ range is given in the case of maximum spin for both prograde (shortest) and retrograde (longest) orbits (Bardeen et al. 1972). eMass accretion rates $\dotM$ are estimated from measurements of Faraday rotation imparted by material in the accretion flow around the black hole.
References. (1) GRAVITY Collaboration et al. (2018a), (2) Reid et al. (2014), (3) Lu et al. (2018), (4) Marrone et al. (2007), (5) Walsh et al. (2013), (6) Gebhardt et al. (2011), (7) Blakeslee et al. (2009), EHT Collaboration et al. (2019e), (8) Akiyama et al. (2015), (9) Kuo et al. (2014), (10) Reid & Brunthaler (2004), (11) Lambert & Gontier (2009).
Within the past decade, VLBI observations at a wavelength of 1.3 mm have confirmed the existence of structure on the scale of the shadow in both the nucleus of M87 and Sgr A*. M87 (Virgo A) is thought to harbor a black hole in the range of $3.3\mbox—6.2\times 10^9$ solar mass (M⊙) (Gebhardt et al. 2011 ; Walsh et al. 2013), when scaled to a distance of 16.8 Mpc (Paper VI). Structure in M87 measuring 40 μas in extent has been resolved with 1.3 mm VLBI, corresponding to a total extent of 5.5 Schwarzschild radii at the upper end of the mass range (Doeleman et al. 2012 ; Akiyama et al. 2015). For Sgr A*, the black hole candidate at the Galactic center with a presumed mass of $4\times 10^6\,M_\odot $ (Balick & Brown 1974 ; Ghez et al. 2008 ; Genzel et al. 2010 ; GRAVITY Collaboration et al. 2018a), the 1.3 mm emission has been measured to have a size of 3.7 Rs (Doeleman et al. 2008 ; Fish et al. 2011). More recently, full polarimetric VLBI observations at 1.3 mm wavelength have revealed ordered and time-variable magnetic fields within Sgr A* on horizon scales (Johnson et al. 2015), and extension to longer baselines has confirmed compact structure on 3 Rs scales (Lu et al. 2018). These results, obtained with three- and four-site VLBI arrays consisting of the former Combined Array for Research in Millimeter-wave Astronomy (CARMA) in California, the Submillimeter Telescope (SMT) in Arizona, the James Clerk Maxwell Telescope (JCMT) and Submillimeter Array (SMA) facilities on Maunakea in Hawaii, and the Atacama Pathfinder Experiment (APEX) telescope in Chile, demonstrated that direct imaging of emission structures near the event horizon of SMBH candidates is possible in principle. For comparison and perspective, the closest approach of the orbiting stars used to determine the mass of Sgr A* is 1400 Rs (Gravity Collaboration et al. 2018b).
1.3. Array Architecture and Context
To realize these fundamental science goals, our international collaboration has engineered the EHT to move beyond the detection of horizon-scale structure and achieve the required imaging and time-domain sampling capability.
One of the key enabling technologies behind the EHT observations has been the development of high-bandwidth (wideband) VLBI systems that compensate to some degree for the generally smaller telescope apertures at millimeter and submillimeter wavelengths. The first detections of horizon-scale structure followed directly from deployment of new digital VLBI backend and recording instrumentation, custom-built for mm-wavelength observations that achieved a recording rate of 4 gigabit s−1 ( Gbps)136 (Doeleman et al. 2008). Continued development led to an increased recording rate of 16 Gbps (Whitney et al. 2013). Adoption of industry-standard high-speed data protocols, increased hard disk storage capacity, and flexible field programmable gate array (FPGA) computational fabric have enabled the EHT to reach data throughputs of 64 Gbps (Vertatschitsch et al. 2015), or 32 times the maximum recording rate and corresponding bandwidth, offered by open access VLBI facilities at longer wavelengths (e.g., the NRAO Very Long Baseline Array ; Napier et al. 1994).
In addition to the increased sensitivity provided by such large data recording rates, several factors, some engineered and some serendipitous, have converged to enable spatially and temporally resolved observations of black hole candidates by the EHT. By temporal resolution we mean that increased resolution, sensitivity, and baseline coverage allows the EHT to detect and spatially resolve horizon-scale time-variable structures, which would otherwise only be studied through light-curve analysis and light-crossing time assumptions. A list of the key enabling factors is as follows.
1.
Angular resolution : the angular resolution of Earth-diameter VLBI baselines at wavelengths of 1.3 mm can resolve the lensed photon orbits of Sgr A* and M87 (Table 1 : about 50 μas and 38 μas, respectively).
2.
Fourier coverage : using Earth-rotation aperture synthesis, the number of existing and planned mm/submm wavelength telescopes provides a sufficient sampling of VLBI baseline lengths and orientations to produce images with horizon-scale resolution (Fish et al. 2014 ; Honma et al. 2014 ; Lu et al. 2014 ; Ricarte & Dexter 2015 ; Bouman et al. 2016 ; Chael et al. 2016).
3.
Atmospheric transparency : at the required mm/submm observing wavelengths, the Earth’s atmosphere at high-altitude sites is reliably transparent enough that global VLBI arrays can be formed for long-duration observations (Thompson et al. 2017).
4.
Optically thin accretion : for both Sgr A* and M87 the spectral energy density of the accretion flow begins to turn over at mm wavelengths, allowing photons from deep within the gravitational potential well to escape, presuming a synchrotron emission mechanism (see Broderick & Loeb 2009 ; Genzel et al. 2010, for M87 and Sgr A*, respectively).
5.
Interstellar scattering : radio images of Sgr A* are blurred due to interstellar scattering by free electrons (Lo et al. 1998 ; Shen et al. 2005 ; Bower et al. 2006 ; Lu et al. 2011 ; Bower et al. 2015 ; Johnson et al. 2018 ; Psaltis et al. 2018). This blurring decreases with wavelength as λ2 and becomes sub-dominant for wavelengths of 1.3 mm and shorter, where observations enable direct access to intrinsic structures in close proximity to the event horizon (Doeleman et al. 2008 ; Issaoun et al. 2019).
1.4. Current EHT Array
Figure 1 shows a map of the EHT array. In 2017 April, the EHT carried out global observations with an array of eight telescopes (see Table 2) that included the Atacama Large Millimeter/submillimeter Array (ALMA) for the first time. A purpose-built system electronically combined the collecting area of 37 × 12 m diameter ALMA dishes (see Appendix A.1) : the equivalent of adding a 70 m dish to the EHT array. Other participating telescopes were APEX, JCMT, SMA, SMT, the Large Millimeter Telescope Alfonso Serrano (LMT), the Pico Veleta 30 m telescope (PV), and the South Pole Telescope (SPT). Operating in the 1.3 mm window in full polarimetric mode and with an aggregate bandwidth of 8 GHz, the resulting increase in sensitivity above the first horizon-scale detections was nearly an order of magnitude (e.g., Section 3.8). For observations with phased-ALMA in 2018, the EHT added an additional facility (the 12 m diameter Greenland Telescope (GLT)) and doubled the aggregate bandwidth to its nominal target of 16 GHz.
figure1
Map of the EHT. Stations active in 2017 and 2018 are shown with connecting lines and labeled in yellow, sites in commission are labeled in green, and legacy sites are labeled in red. Nearly redundant baselines are overlaying each other, i.e., to ALMA/APEX and SMA/JCMT. Such redundancy allows improvement in determining the amplitude calibration of the array (Paper III).
Table 2. EHT Station Information
aGeocentric coordinates with X pointing to the Greenwich meridian, Y pointing 90° away in the equatorial plane (eastern longitudes have positive Y), and positive Z pointing in the direction of the North Pole. This is a left-handed coordinate system. Elevations are relative to the GRS 80 ellipsoid (Moritz 2000). bArray coordinates indicate the phasing center used in 2017. c2017 April position : the ice sheet at the South Pole moves at a rate of about 10 m yr−1. Effects of this slow drift are removed during VLBI correlation. dNOEMA : Northern Extended Millimeter Array ; KP 12 m : Kitt Peak 12 m. These stations have not participated in global EHT observations and their coordinates are approximate.
American Astronomical Society logo
American Astronomical Society logo
iop-2016.png
iop-2016.png
A publishing partnership
First M87 Event Horizon Telescope Results. II. Array and Instrumentation
The Event Horizon Telescope Collaboration, Kazunori Akiyama1,2,3,4, Antxon Alberdi5, Walter Alef6, Keiichi Asada7, Rebecca Azulay8,9,6, Anne-Kathrin Baczko6, David Ball10, Mislav Baloković4,11, John Barrett2Show full author list
Published 2019 April 10 • © 2019. The American Astronomical Society.
The Astrophysical Journal Letters, Volume 875, Number 1
Focus on the First Event Horizon Telescope Results
Download Article PDF DownloadArticle ePub
Figures
Tables
References
Download PDF Download ePub
Turn on MathJax
Get permission to re-use this article
Share this article
Share this content via email
Share on Facebook
Share on Twitter
Share on Google+
Share on CiteULike
Share on Mendeley
Article information
Abstract
The Event Horizon Telescope (EHT) is a very long baseline interferometry (VLBI) array that comprises millimeter- and submillimeter-wavelength telescopes separated by distances comparable to the diameter of the Earth. At a nominal operating wavelength of 1.3 mm, EHT angular resolution (λ/D) is 25 μas, which is sufficient to resolve nearby supermassive black hole candidates on spatial and temporal scales that correspond to their event horizons. With this capability, the EHT scientific goals are to probe general relativistic effects in the strong-field regime and to study accretion and relativistic jet formation near the black hole boundary. In this Letter we describe the system design of the EHT, detail the technology and instrumentation that enable observations, and provide measures of its performance. Meeting the EHT science objectives has required several key developments that have facilitated the robust extension of the VLBI technique to EHT observing wavelengths and the production of instrumentation that can be deployed on a heterogeneous array of existing telescopes and facilities. To meet sensitivity requirements, high-bandwidth digital systems were developed that process data at rates of 64 gigabit s−1, exceeding those of currently operating cm-wavelength VLBI arrays by more than an order of magnitude. Associated improvements include the development of phasing systems at array facilities, new receiver installation at several sites, and the deployment of hydrogen maser frequency standards to ensure coherent data capture across the array. These efforts led to the coordination and execution of the first Global EHT observations in 2017 April, and to event-horizon-scale imaging of the supermassive black hole candidate in M87.
Export citation and abstract BibTeX RIS
Related links
cc-by
Original 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.
1. Introduction
It is generally accepted that active galactic nuclei (AGNs) are powered by accretion onto supermassive black holes (SMBHs ; Heckman & Best 2014). These central engines are powerful actors on the cosmic stage, with roles in galactic evolution, star formation, mergers, and particle acceleration as evidenced by relativistic jets that both dynamically influence and redistribute matter on galactic scales (Blandford et al. 2018). Inflowing material typically obscures the event horizons of these black hole candidates, but it is in this extreme environment of the black hole boundary that strong-field effects of general relativity become evident and the accretion and outflow processes that govern black hole feedback on galactic scales originate (Ho 2008). Imaging black holes on scales that resolve these effects and processes would enable new tests of general relativity and the extraordinarily detailed study of core AGN physics. Realization of this goal requires a specialized instrument that does two things. It must have the ultra-high angular resolution required to resolve the nearest SMBH candidates, and it must operate in a range of the electromagnetic spectrum where light streams unimpeded from the innermost accretion region to telescopes on Earth. Achieving these specifications is the primary objective of the Event Horizon Telescope (EHT) : a very long baseline interferometry (VLBI) array of millimeter (mm) and submillimeter (submm) wavelength facilities that span the globe, creating a telescope with an effective Earth-sized aperture (Doeleman et al. 2009a).
While the EHT is uniquely designed for the imaging of SMBHs, other pioneering instruments are capable of probing similar angular scales for other purposes. Operating in the infrared, the GRAVITY interferometer delivers relative astrometry at the 10 micro-arcsecond (μas) level and has provided evidence for relativistic motion of material in close proximity to Sgr A* (GRAVITY Collaboration et al. 2018a). These infrared observations are an important and parallel probe of the spacetime surrounding Sgr A*, but cannot be used to make spatially resolved images of the black hole candidate because the interferometer intrinsic resolution is only 3 mas. The RadioAstron satellite, used as an orbiting element of combined Earth-Space VLBI arrays, is also capable of 10 μas angular resolution (Kardashev et al. 2013), but it operates at longer radio wavelengths that cannot penetrate the self-absorbed synchrotron plasma that surrounds the event horizon. In many ways the EHT is complementary to the Laser Interferometer Gravitational-Wave Observatory (LIGO) facility, which has detected the gravitational wave signatures from merging stellar-mass black holes (Abbott et al. 2016). LIGO and the EHT observe black holes that differ in mass by factors of $10^4\mbox—10^7$ ; LIGO events are transient, while the EHT carries out long-term studies of its main targets.
This Letter is one in a sequence of manuscripts that describes the first EHT results. The full sequence includes an abstract Letter with a summary of results (EHT Collaboration et al. 2019a, hereafter Paper I), this Letter with a description of the array and instrumentation (Paper II), a description of the data pipeline and processing (EHT Collaboration et al. 2019b, hereafter Paper III), a description of imaging techniques (EHT Collaboration et al. 2019c, hereafter Paper IV), and theoretical analyses of astrophysical and physics results (EHT Collaboration et al. 2019d, hereafter Paper V ; EHT Collaboration et al. 2019e, hereafter Paper VI, respectively).
1.1. EHT Science Goals
The scientific goals of the EHT array are to resolve and detect strong general relativistic signatures expected to arise on event-horizon scales. The best-known effect is that a black hole, surrounded by an optically thin luminous plasma, should exhibit a ’’silhouette’’ or ’’shadow’’ morphology : a dim central region delineated by the lensed photon orbit (Falcke et al. 2000). The apparent size of the photon orbit, described not long after Schwarzschild’s initial solution was published (Hilbert 1917 ; von Laue 1921), defines a bright ring or crescent shape that was calculated for arbitrary spin by Bardeen (1973), first imaged through simulations by Luminet 1979, and subsequently studied extensively (Chandrasekhar 1983 ; Takahashi 2004 ; Broderick & Loeb 2006). The size and shape of the resulting shadow depends primarily on the mass of the black hole, and only very weakly on its spin and the observing orientation. For a non-spinning black hole, the shadow diameter is equal to $\sqrt27$ Schwarzschild radii (Rs = 2GM/c2). Over all black hole spins and orientations, the shadow diameter ranges from 4.8 to 5.2 Rs (Bardeen 1973 ; Johannsen & Psaltis 2010). Because any light that crosses the photon orbit from outside will eventually reach the event horizon, use of the term ’’horizon scale’’ will hereafter be understood to mean the size of the shadow and the lensed photon orbit. Detection of the shadow via lensed electromagnetic radiation would provide new evidence for the existence of SMBHs by confining the masses of EHT targets to within their expected photon orbits. A more detailed study of the precise shape of the photon orbit can be used to test the validity of general relativity on horizon scales (Johannsen & Psaltis 2010). Full polarimetric imaging can similarly be used to map magnetic field structure near the event horizon, placing important constraints on modes of accretion and the launching of relativistic jets (Broderick & Loeb 2009 ; Johnson et al. 2015 ; Chael et al. 2016 ; Akiyama et al. 2017 ; Gold et al. 2017).
Separate signatures, potentially offering a more sensitive probe of black hole spin, are the timescales of dynamical processes at horizon scales. A characteristic timescale for such processes is given by the orbital period of test particles at the innermost stable circular orbit (ISCO), which depends sensitively on the black hole spin (Bardeen et al. 1972). Monitoring VLBI observables to track the orbital dynamics of inhomogeneities in the accretion flow can thus be used to probe the spacetime, and potentially spin, of the black hole (Broderick & Loeb 2006 ; Doeleman et al. 2009b ; Fish et al. 2009 ; Fraga-Encinas et al. 2016 ; Medeiros et al. 2017 ; Roelofs et al. 2017).
1.2. Target Sources and Confirmation of Horizon-scale Structure
The bright radio core of M87 and Sgr A* (Table 1) are the primary EHT targets, as the combination of their estimated mass and proximity make them the two most suitable sources for studying SMBH candidates at horizon-scale resolution (Narayan & McClintock 2008 ; Johannsen et al. 2012). With bolometric luminosities well below the Eddington limit (Di Matteo et al. 2000 ; Baganoff et al. 2003), both the nucleus of M87 and Sgr A* are representative of the broad and populous class of low-luminosity AGN (LLAGN). AGN spend most of their time in this low state after prior periods of high accretion (Ho 2008), and LLAGN share many characteristics of stellar black hole emission in X-ray binaries, which exhibit episodic emission and jet production (Narayan & McClintock 2008). Thus, the prospects for applying results from horizon-scale observations of the central region of M87 and Sgr A* across a broad range of astrophysical contexts are excellent.
Table 1. Assumed Physical Properties of Sgr A* and M87 Used to Establish Technical Goalsa
Sgr A* M87
Black Hole Mass M (M⊙) $4.1\times 10^6$ (1) (3.3–$6.2)\times 10^9$ (5), (6)
Distance D (pc) $8.34\times 10^3$ (2) $16.8\times 10^6$ (7)
Schwarzschild Radius Rs (μas) 9.7 3.9–7.3
Shadow Diameterb $D_\mathrmsh$ (μas) 47–50 19–38
Brightness Temperaturec TB (K) $3\times 10^9$ (3) 1010 (8)
Period ISCOd $P_\mathrmISCO$ 4–54 minutes 2.4–57.7 days
Mass Accretion Ratee $\dotM$ (M⊙ yr−1) 10−9–10−7 (4) $\lt 10^-3$ (9)
Notes.
aSgr A* : $\alpha _\rmJ2000.0=17^\rmh
45^\rmm
40\buildrel\rms\over. 0409,\,\delta _\rmJ2000.0=-29^\circ 00^\prime 28\buildrel\prime\prime\over. 118$ (10) ; M87 : $\alpha _\rmJ2000.0=12^\rmh
30^\rmm
49\buildrel\rms\over. 4234,\,\delta _\rmJ2000.0=12^\circ 23^\prime 28\buildrel\prime\prime\over. 044$ (11). bThe shadow diameter is within the range 4.8–5.2 Rs depending on black hole spin and orientation to the observer’s line of sight (Johannsen & Psaltis 2010). cBrightness temperatures are reported for an observing frequency of 230 GHz. d $P_\mathrmISCO$ range is given in the case of maximum spin for both prograde (shortest) and retrograde (longest) orbits (Bardeen et al. 1972). eMass accretion rates $\dotM$ are estimated from measurements of Faraday rotation imparted by material in the accretion flow around the black hole.
References. (1) GRAVITY Collaboration et al. (2018a), (2) Reid et al. (2014), (3) Lu et al. (2018), (4) Marrone et al. (2007), (5) Walsh et al. (2013), (6) Gebhardt et al. (2011), (7) Blakeslee et al. (2009), EHT Collaboration et al. (2019e), (8) Akiyama et al. (2015), (9) Kuo et al. (2014), (10) Reid & Brunthaler (2004), (11) Lambert & Gontier (2009).
Download table as : ASCIITypeset image
Within the past decade, VLBI observations at a wavelength of 1.3 mm have confirmed the existence of structure on the scale of the shadow in both the nucleus of M87 and Sgr A*. M87 (Virgo A) is thought to harbor a black hole in the range of $3.3\mbox—6.2\times 10^9$ solar mass (M⊙) (Gebhardt et al. 2011 ; Walsh et al. 2013), when scaled to a distance of 16.8 Mpc (Paper VI). Structure in M87 measuring 40 μas in extent has been resolved with 1.3 mm VLBI, corresponding to a total extent of 5.5 Schwarzschild radii at the upper end of the mass range (Doeleman et al. 2012 ; Akiyama et al. 2015). For Sgr A*, the black hole candidate at the Galactic center with a presumed mass of $4\times 10^6\,M_\odot $ (Balick & Brown 1974 ; Ghez et al. 2008 ; Genzel et al. 2010 ; GRAVITY Collaboration et al. 2018a), the 1.3 mm emission has been measured to have a size of 3.7 Rs (Doeleman et al. 2008 ; Fish et al. 2011). More recently, full polarimetric VLBI observations at 1.3 mm wavelength have revealed ordered and time-variable magnetic fields within Sgr A* on horizon scales (Johnson et al. 2015), and extension to longer baselines has confirmed compact structure on 3 Rs scales (Lu et al. 2018). These results, obtained with three- and four-site VLBI arrays consisting of the former Combined Array for Research in Millimeter-wave Astronomy (CARMA) in California, the Submillimeter Telescope (SMT) in Arizona, the James Clerk Maxwell Telescope (JCMT) and Submillimeter Array (SMA) facilities on Maunakea in Hawaii, and the Atacama Pathfinder Experiment (APEX) telescope in Chile, demonstrated that direct imaging of emission structures near the event horizon of SMBH candidates is possible in principle. For comparison and perspective, the closest approach of the orbiting stars used to determine the mass of Sgr A* is 1400 Rs (Gravity Collaboration et al. 2018b).
1.3. Array Architecture and Context
To realize these fundamental science goals, our international collaboration has engineered the EHT to move beyond the detection of horizon-scale structure and achieve the required imaging and time-domain sampling capability.
One of the key enabling technologies behind the EHT observations has been the development of high-bandwidth (wideband) VLBI systems that compensate to some degree for the generally smaller telescope apertures at millimeter and submillimeter wavelengths. The first detections of horizon-scale structure followed directly from deployment of new digital VLBI backend and recording instrumentation, custom-built for mm-wavelength observations that achieved a recording rate of 4 gigabit s−1 ( Gbps)136 (Doeleman et al. 2008). Continued development led to an increased recording rate of 16 Gbps (Whitney et al. 2013). Adoption of industry-standard high-speed data protocols, increased hard disk storage capacity, and flexible field programmable gate array (FPGA) computational fabric have enabled the EHT to reach data throughputs of 64 Gbps (Vertatschitsch et al. 2015), or 32 times the maximum recording rate and corresponding bandwidth, offered by open access VLBI facilities at longer wavelengths (e.g., the NRAO Very Long Baseline Array ; Napier et al. 1994).
In addition to the increased sensitivity provided by such large data recording rates, several factors, some engineered and some serendipitous, have converged to enable spatially and temporally resolved observations of black hole candidates by the EHT. By temporal resolution we mean that increased resolution, sensitivity, and baseline coverage allows the EHT to detect and spatially resolve horizon-scale time-variable structures, which would otherwise only be studied through light-curve analysis and light-crossing time assumptions. A list of the key enabling factors is as follows.
1.
Angular resolution : the angular resolution of Earth-diameter VLBI baselines at wavelengths of 1.3 mm can resolve the lensed photon orbits of Sgr A* and M87 (Table 1 : about 50 μas and 38 μas, respectively).
2.
Fourier coverage : using Earth-rotation aperture synthesis, the number of existing and planned mm/submm wavelength telescopes provides a sufficient sampling of VLBI baseline lengths and orientations to produce images with horizon-scale resolution (Fish et al. 2014 ; Honma et al. 2014 ; Lu et al. 2014 ; Ricarte & Dexter 2015 ; Bouman et al. 2016 ; Chael et al. 2016).
3.
Atmospheric transparency : at the required mm/submm observing wavelengths, the Earth’s atmosphere at high-altitude sites is reliably transparent enough that global VLBI arrays can be formed for long-duration observations (Thompson et al. 2017).
4.
Optically thin accretion : for both Sgr A* and M87 the spectral energy density of the accretion flow begins to turn over at mm wavelengths, allowing photons from deep within the gravitational potential well to escape, presuming a synchrotron emission mechanism (see Broderick & Loeb 2009 ; Genzel et al. 2010, for M87 and Sgr A*, respectively).
5.
Interstellar scattering : radio images of Sgr A* are blurred due to interstellar scattering by free electrons (Lo et al. 1998 ; Shen et al. 2005 ; Bower et al. 2006 ; Lu et al. 2011 ; Bower et al. 2015 ; Johnson et al. 2018 ; Psaltis et al. 2018). This blurring decreases with wavelength as λ2 and becomes sub-dominant for wavelengths of 1.3 mm and shorter, where observations enable direct access to intrinsic structures in close proximity to the event horizon (Doeleman et al. 2008 ; Issaoun et al. 2019).
1.4. Current EHT Array
Figure 1 shows a map of the EHT array. In 2017 April, the EHT carried out global observations with an array of eight telescopes (see Table 2) that included the Atacama Large Millimeter/submillimeter Array (ALMA) for the first time. A purpose-built system electronically combined the collecting area of 37 × 12 m diameter ALMA dishes (see Appendix A.1) : the equivalent of adding a 70 m dish to the EHT array. Other participating telescopes were APEX, JCMT, SMA, SMT, the Large Millimeter Telescope Alfonso Serrano (LMT), the Pico Veleta 30 m telescope (PV), and the South Pole Telescope (SPT). Operating in the 1.3 mm window in full polarimetric mode and with an aggregate bandwidth of 8 GHz, the resulting increase in sensitivity above the first horizon-scale detections was nearly an order of magnitude (e.g., Section 3.8). For observations with phased-ALMA in 2018, the EHT added an additional facility (the 12 m diameter Greenland Telescope (GLT)) and doubled the aggregate bandwidth to its nominal target of 16 GHz.
Figure 1.
Figure 1. Map of the EHT. Stations active in 2017 and 2018 are shown with connecting lines and labeled in yellow, sites in commission are labeled in green, and legacy sites are labeled in red. Nearly redundant baselines are overlaying each other, i.e., to ALMA/APEX and SMA/JCMT. Such redundancy allows improvement in determining the amplitude calibration of the array (Paper III).
Download figure :
Standard image High-resolution image Export PowerPoint slide
Table 2. EHT Station Information
Facility Diameter Location Xa Ya Za Latitude Longitude Elevationa
(m) (m) (m) (m) (m)
Facilities that Participated in the 2017 Observations
ALMAb 12 ($\times 54$) and 7 ($\times 12$) Chile $2225061.3$ $-5440061.7$ $-2481681.2$ $-23^\circ 01^\prime 45\buildrel\prime\prime\over. 1$ $-67^\circ 45^\prime 17\buildrel\prime\prime\over. 1$ 5074.1
APEX 12 Chile $2225039.5$ $-5441197.6$ $-2479303.4$ $-23^\circ 00^\prime 20\buildrel\prime\prime\over. 8$ $-67^\circ 45^\prime 32\buildrel\prime\prime\over. 9$ 5104.5
JCMT 15 Hawaii, USA $-5464584.7$ $-2493001.2$ $2150654.0$ $+19^\circ 49^\prime 22\buildrel\prime\prime\over. 2$ $-155^\circ 28^\prime 37\buildrel\prime\prime\over. 3$ 4120.1
LMT 50 Mexico $-768715.6$ $-5988507.1$ $2063354.9$ $+18^\circ 59^\prime 08\buildrel\prime\prime\over. 8$ $-97^\circ 18^\prime 53\buildrel\prime\prime\over. 2$ 4593.3
PV 30 m 30 Spain $5088967.8$ $-301681.2$ $3825012.2$ $+37^\circ 03^\prime 58\buildrel\prime\prime\over. 1$ $-3^\circ 23^\prime 33\buildrel\prime\prime\over. 4$ 2919.5
SMAb 6 ($\times 8$) Hawaii, USA $-5464555.5$ $-2492928.0$ $2150797.2$ $+19^\circ 49^\prime 27\buildrel\prime\prime\over. 2$ $-155^\circ 28^\prime 39\buildrel\prime\prime\over. 1$ 4115.1
SMT 10 Arizona, USA $-1828796.2$ $-5054406.8$ $3427865.2$ $+32^\circ 42^\prime 05\buildrel\prime\prime\over. 8$ $-109^\circ 53^\prime 28\buildrel\prime\prime\over. 5$ 3158.7
SPTc 10 Antarctica 809.8 −816.9 $-6359568.7$ $-89^\circ 59^\prime 22\buildrel\prime\prime\over. 9$ $-45^\circ 15^\prime 00\buildrel\prime\prime\over. 3$ 2816.5
Facilities Joining EHT Observations in 2018 and Later
GLT 12 Greenland $541547.0$ $-1387978.6$ $6180982.0$ $+76^\circ 32^\prime 06\buildrel\prime\prime\over. 6$ $-68^\circ 41^\prime 08\buildrel\prime\prime\over. 8$ 89.4
NOEMAd 15 ($\times 12$) France 4524000.4 468042.1 4460309.8 $+44^\circ 38^\prime 01\buildrel\prime\prime\over. 2$ $+5^\circ 54^\prime 24\buildrel\prime\prime\over. 0$ 2617.6
KP 12 md 12 Arizona, USA $-1995954.4$ $-5037389.4$ $3357044.3$ $+31^\circ 57^\prime 12\buildrel\prime\prime\over. 0$ $-111^\circ 36^\prime 53\buildrel\prime\prime\over. 5$ 1894.5
Facilities Formerly Participating in EHT Observations
CARMA 10.4, 6.1 $(\times 8)$ California, USA $-2397378.6$ $-4482048.7$ $3843513.2$ $+37^\circ 16^\prime 49\buildrel\prime\prime\over. 4$ $-118^\circ 08^\prime 29\buildrel\prime\prime\over. 9$ 2168.9
CSO 10 Hawaii, USA $-5464520.9$ $-2493145.6$ $2150610.6$ $+19^\circ 49^\prime 20\buildrel\prime\prime\over. 9$ $-155^\circ 28^\prime 31\buildrel\prime\prime\over. 9$ 4107.2
Notes.
aGeocentric coordinates with X pointing to the Greenwich meridian, Y pointing 90° away in the equatorial plane (eastern longitudes have positive Y), and positive Z pointing in the direction of the North Pole. This is a left-handed coordinate system. Elevations are relative to the GRS 80 ellipsoid (Moritz 2000). bArray coordinates indicate the phasing center used in 2017. c2017 April position : the ice sheet at the South Pole moves at a rate of about 10 m yr−1. Effects of this slow drift are removed during VLBI correlation. dNOEMA : Northern Extended Millimeter Array ; KP 12 m : Kitt Peak 12 m. These stations have not participated in global EHT observations and their coordinates are approximate.
Download table as : ASCIITypeset image
The sections that follow describe the specifications and characteristics of the array (Section 2), the EHT instrumentation deployed (Section 3), the observing strategy (Section 4), correlation, calibration, and detection (Section 5), and future enhancements (Section 6).
2. EHT Specifications and Characteristics
Extending the VLBI technique to wavelengths of 1.3 mm presents technical challenges. Heterodyne receivers exhibit greater noise, and the required stability of atomic frequency standards is higher than that typically specified for VLBI at longer wavelengths. Early 1.3 mm wavelength VLBI experiments and observations in the 1990s succeeded in making first detections of AGNs and Sgr A* on modest-length baselines (≤1100 km ; Padin et al. 1990 ; Greve et al. 1995 ; Krichbaum et al. 1998). Following this pioneering work (see Doeleman & Krichbaum 1999 ; Boccardi et al. 2017, for summaries), 1.3 mm VLBI efforts over the next decade focused on order-of-magnitude bandwidth expansion as a means to boost sensitivity. This development led to observations in 2007 with enough resolution to resolve emission on the scale of the event horizon of Sgr A* by using a three-station array with telescopes in Hawaii, California, and Arizona. Motivated by this detection, an EHT array capable of imaging strong general relativistic features was planned based on organizing coordinated observations on a network of mm-wavelength observatories (Doeleman et al. 2009a). These sites and their general characteristics are given in Table 2, which lists current and planned EHT sites as well as now-decommissioned observatories used in prior experiments. For the most part, the EHT consists of pre-existing telescopes that perform single-dish or connected-element astronomy observations during most of the year, but which required modifications or upgrades (in some cases significant ones) to carry out VLBI. The GLT, for example, was commissioned primarily to image M87 (Inoue et al. 2014), and a heterodyne receiver for the SPT was built specifically for EHT observations (Kim et al. 2018a). Details of how each site was modified for EHT work are given in the Appendix.
Technical specifications for the EHT were adopted to unify VLBI recordings across the heterogeneous array by establishing a common frequency configuration, polarization configuration, and sampling rate (Tilanus et al. 2013 ; Marrone et al. 2014). When ALMA participates in EHT observations, it is by far the most sensitive site, so the overall sensitivity of the EHT array is optimized on a per-bandwidth basis when all other sites can match the recorded frequency bands at ALMA. The EHT converged on a scheme that matches ALMA specifications : two 4 GHz sidebands in each of two polarizations for the 1.3 and 0.87 mm receiver bands, which could be realized through feasible modifications and enhancements at most sites. The resulting global array geometry and sensitivity is well matched to the science goals.
2.1. Angular Resolution
Imaging a black hole shadow requires several resolution elements in each direction across the lensed innermost photon orbit, in addition to field-of-view coverage that extends beyond the feature. The longest baselines of the EHT (e.g., South Pole to Arizona, Hawaii, or Spain) provide nominal angular resolutions of $\lambda /D\simeq 25\ $μas in the 1.3 mm wavelength band. Regularized maximum likelihood (RML) imaging methods (see Paper IV) typically achieve angular resolutions that are better than the nominal figure by factors of 2–3 (Narayan & Nityananda 1986). For the EHT case in particular, RML methods have been extensively tested using realistic and synthetic interferometric data to set the optimal resolution of the array (Honma et al. 2014 ; Bouman et al. 2016 ; Chael et al. 2016 ; Akiyama et al. 2017 ; Chael et al. 2018 ; Kuramochi et al. 2018). This results in an anticipated effective EHT angular resolution at 1.3 mm of 20 μas, yielding between about 36 resolution elements across a 120 × 120 μas field of view. For Sgr A* and M87 this field of view is expected to encompass the dim interior, bright annulus of the photon orbit, and sufficient spatial extent to extract the shadow feature (e.g., Mościbrodzka et al. 2009 ; Dexter et al. 2010). These considerations, and the expectation based on simulations that the EHT instrument and array could achieve such resolution, were important factors in specifying the EHT architecture. Further tests using EHT observations of quasar calibrator sources (e.g., 3C 279) obtained in the 2017 April observations demonstrate robust structural agreement between RML methods and more traditional CLEAN-based radio imaging techniques. Similar comparisons and results on M87, one of the EHT primary targets, can be found in Paper IV.
2.2. Sensitivity
Maximizing detections across sparse Fourier coverage is essential for high-fidelity image reconstruction. This is especially true for Sgr A* because its structure is scatter broadened by interstellar effects resulting in reduced VLBI visibility amplitudes on the longest baselines.
To maximize the number of interferometric detections across the array, the EHT uses a two-stage approach to fringe detection. First, detections are found on baselines from all stations to ALMA within the atmospheric coherence time. Next, these ALMA detections are used to remove the effects of phase fluctuations (due to atmospheric turbulence) above the non-ALMA stations, allowing coherent integration on non-ALMA baselines for intervals up to many minutes, thereby boosting the signal-to-noise ratio (S/N) to recover the full baseline coverage of the array (Paper III). Thus, the EHT sensitivity specification corresponds to the requirement that for baselines connecting each EHT site to ALMA, Sgr A* and M87 are detected with a typical signal-to-noise that enables this atmospheric phase correction with acceptable loss. In the usual case where the VLBI observing scan duration greatly exceeds the atmospheric coherence time, the loss due to noise in the phase-correction algorithm is $\sim e^
\rm-(S/N)
^2/2$, where the S/N is the signal-to-noise on the ALMA baselines. To ensure negligible loss, we specify S/N > 3 for EHT baselines to ALMA within an integration time, $T_\rmint
$, where Tint is less than the atmospheric coherence time. We note that the expected change in interferometric phase due to source structure over the imaged field of view (Section 2.1) would be less than a few degrees over typical VLBI scan lengths of a few minutes.
The S/N of a VLBI signal on a single baseline between stations is
Equation (1)
where ηQ is the digital loss due to sampling the received signal at each antenna with finite precision ($\eta _Q=0.88$ for 2-bit samples), $\rm\Delta \nu $ is the bandwidth of the recording, SEFDi is the system equivalent flux density (SEFD)137 for station i, and $S_\rmcor
$ is the expected correlated flux on the baseline between stations 1 and 2. $T_\rmint
$ is the integration interval of the VLBI signal. It is typically much less than the atmospheric coherence time (Tcoh), or the coherent integration time beyond which the VLBI visibility138 signal decreases by 10 % due to phase fluctuations imposed by turbulence in the troposphere (typically from a few to 20 s). On the weakest baseline involving ALMA, ALMA–SPT, with an estimated flux for Sgr A* of 0.1 Jy (see APEX-CARMA baseline in Lu et al. 2018) and with Δν = 4 GHz, an integration time of 3 s yields an S/N of about 12 (see Table 3 for parameters), which far exceeds the EHT S/N specification for detections on ALMA baselines. Figure 2 shows that more than 75% of scans on baselines in the 2017 array that include ALMA had Tcoh greater than 10 s.
Figure 2. Atmospheric coherence time for EHT ALMA baselines. VLBI observations consist of scans, each 3–7 minutes long (see Figure 12). Each scan in the 2017 EHT observation campaign was divided into sequential coherent integration intervals, which were then incoherently averaged to obtain an estimate of VLBI signal amplitude per scan. As the duration of the coherent interval decreases, the loss in amplitude due to atmospheric phase fluctuations decreases. Solid lines represent the 25th percentile of the coherence distribution for each baseline : 75% of the data exhibit higher time coherence than indicated by the line. The heavy black dashed line shows the 25th percentile of the coherence distribution for all baselines. The light dashed and dotted lines show the 95 % and 90 % coherence limits for the baseline with lowest coherence, which is from PV to ALMA (largely due to low-elevation PV scans on Sgr A* during the observations). The number of scans used to compute the coherence curves are shown next to each station in parenthesis.
Table 3. EHT Station Sensitivities
aSEFD = 2k$T_\rmsys
^* $/$\eta _\rmA
A_\rmgeom
$ with $T_\rmsys
^* $ the effective system temperature and $A_\rmgeom
$ the geometric collecting area of the telescope. The effective system temperature lists representative values for the 2017 Sgr A* and M87 observations, i.e., for median elevations of the sources and typical atmospheric opacities at each facility. It includes effects from elevation-dependent gains and phasing efficiency relevant for some of the telescopes. bWith σrms-ALMA the thermal noise for the station’s baseline with ALMA using Equation (1) for a bandwidth of 4 GHz and integration time of 10 s (see Figure 2). cThe median number of phased antennas in 2017 for the SMA was 6 × 6 m antennas and for ALMA was 37 × 12 m antennas. If all ALMA antennas are phased that would be equivalent to a single 88 m antenna. dThe LMT was upgraded from 32.5 m to a full 50 m diameter telescope in early 2018. At SPT the aperture was under-illuminated in 2017 with an effective diameter of 6 m.
2.3. Fourier Coverage
An interferometer like the EHT samples the Fourier transform of the image on the sky. By correlating the data obtained from N stations, $N(N-1)/2$ spatial frequencies are measured. As the Earth rotates, those spatial frequencies form tracks in the Fourier plane (i.e., the (u, v) plane) to produce a sparsely sampled Fourier transform of the sky image. With the sensitivity requirements established, one can estimate the baseline coverage that results after the two-stage detection process above (Section 2.2) is followed. Figure 3 shows that the measured spatial frequencies that satisfy the EHT signal-to-noise specification when ALMA is included in the array result in near-full coverage on all baselines in the array, thereby maximizing imaging potential. At the outset of EHT build-out, this full coverage goal served as a design goal based on imaging simulations of synthetic data (see Section 2.1). Subsequent observations in 2017 April confirm that this (u, v) coverage is sufficient to image horizon-scale features (Paper IV).
Figure 3. Expected EHT Fourier space coverage on Sgr A*. The left panel shows both detections (red) and non-detections (gray) of Sgr A* in the 2013 EHT campaign. Participating telescopes were : APEX, CARMA, JCMT, SMA, and SMT. The dashed circles mark baselines with a fringe spacing equal to 50 μas (approximately the diameter of the shadow of the SMBH candidate Sgr A*) and 25 μas. The two remaining panels show simulated EHT observations in 2020 : (1) without ALMA and (2) with ALMA. Specifications to determine baseline detections shown are those detailed in Section 2.2. These figures emphasize the benefit of including ALMA in the array : its high sensitivity allows detections for SgrA* on all observed baselines. Because each EHT site requires at least one strong baseline to identify an interferometric fringe and to correct for residual delays, rates, and phase wander, ALMA significantly extends the Fourier coverage and sensitivity even for non-ALMA baselines. Coverage shown in the two right panels corresponds to an array including Chile (APEX, ALMA), Mexico (LMT), France (NOEMA), Spain (PV), Hawaii (SMA, JCMT), Arizona (SMT), and South Pole (SPT). The corresponding baseline coverage of the 2017 observations is shown in Figure 11.
2.4. Time Resolution
Characteristic timescales that affect the evolution of structures in horizon-scale images include the light-crossing time ($T_\rmlight
=GM/c^3$) and the period of the ISCO (PISCO), both of which scale with mass. With an assumed mass of $4.1\times 10^6M_\odot $ for Sgr A*, $T_\rmlight
=20\,\rms$ and PISCO ranges from 4 minutes for a prograde orbit around a maximally spinning black hole to 54 minutes for a retrograde orbit around a maximally spinning black hole. If a mass of $6.2\times 10^9M_\odot $ is assumed for M87, $T_\rmlight
=8.5\,\rmh$ and PISCO ranges from 4.5 to 58 days for the same orbits. All of these timescales, while important for imaging and modeling time-variable structures, far exceed the coherence times of the atmosphere, which set the integration intervals and sensitivity requirements for the EHT. Hence, the EHT must already track VLBI observables on considerably finer timescales than Tlight or PISCO for the primary cosmic targets, M87 and Sgr A*.
2.5. Frequency Configuration
The instrumentation at EHT sites is designed to match ALMA’s bandwidth and intermediate frequency (IF) ranges. ALMA antennas are equipped with dual-polarization, sideband-separating receivers with four IF outputs (Section 3.1). The total bandwidth that can be transported back from each ALMA antenna is 16 GHz, thus resulting in a maximum IF bandwidth of 4 GHz per output. The sky frequencies for both the lower- and upper-sideband must overlap perfectly from station to station, which fixes both the local oscillator (LO) frequency and IF frequency ranges at each observatory (Tilanus et al. 2013). Table 4 shows the frequency configurations of the EHT in both atmospheric transmission windows around 230 and 345 GHz.
Table 4. EHT Frequency Configurations
The main considerations for selecting the chosen LO frequencies are (Marrone et al. 2014) as follows.
1.
The tuning range of the receivers at participating facilities.
2.
Atmospheric transmission.
3.
The avoidance of Galactic 12CO : $\nu \,\sim $ 230.3–230.8 GHz and $\nu \,\sim $ 345.4–346.1 GHz within the VLBI observing bands.
4.
To a lesser degree : avoidance of Galactic 13CO : $\nu \,\sim $ 220.2–220.6 GHz and $\nu \,\sim $ 330.3–330.9 GHz within the VLBI observing bands.
5.
Access to 12CO spectral lines in an extended tuning range of the receivers above 9 GHz when observing in the 1.3 mm band, or below 4 GHz when observing in the 0.87 mm band.
6.
Access to maser lines within the VLBI band, e.g., the SiO maser line at 215.596 GHz ($v=1,\,J=5\to \,4$).
7.
Performance of existing quarter-wave plates used to observe circular polarization.
For ALMA’s 230 GHz band, the IF band previously was restricted to a lower limit of 5 GHz (now 4.5 GHz), which resulted in a common IF range across the telescopes of 5–9 GHz, whereas the common IF range for the 345 GHz band is 4–8 GHz.
3. Instrumentation
A schematic of the EHT’s VLBI signal chain at single-dish telescopes is shown in Figure 4. The front end is typically a dual-polarization sideband-separating receiver in the 1.3 mm or 0.87 mm bands, often the product of a joint development project between the EHT and the telescope facility. These efforts are described in the Appendix. A hydrogen maser provides a frequency reference standard of sufficient stability for mm-VLBI (Section 3.2), and is used to phase lock all analog systems as well as digital sampling clocks throughout the signal chain.
Figure 4. System diagram at a single-dish EHT site with 64 Gbps capability. Dual-polarization (left and right circular polarization (LCP and RCP)) and upper- and lower-sideband (USB and LSB) analog IF signals are sent from the receiver. The receiver local oscillator is locked to the station hydrogen maser 10 MHz reference. Block downconverters (BDCs) mix these signals to baseband. R2DBEs sample the analog signals and distribute packetized data to Mark 6 recorders over a 10 GbE network. All timing is locked to a 10 MHz maser reference and synchronized with a pulse-per-second (PPS) Global Positioning System (GPS) signal. The components are controlled through a 1 GbE network. See Figure 8 for a photograph.
Dedicated block downconverters (BDCs ; Section 3.3) mix IF bands coming from the receivers to baseband. Each 4 GHz wide IF band is split and downconverted into two 2 GHz wide sections at baseband. High-bandwidth digital backends (DBEs ; Section 3.4) are used to sample two 2 GHz baseband signals each. The sampled data are put into the VLBI Data Interchange Format (VDIF ; Whitney et al. 2009) and transmitted to a Mark 6 recorder (Section 3.5) via two 10 Gbps Ethernet (GbE) small form-factor pluggable (SFP+) interfaces.139 The 16 Gbps output from each DBE matches the recording rate of a single Mark 6 VLBI recorder that records in parallel to 32 hard drives (Whitney et al. 2013). Four DBE—Mark 6 systems are used to sample and record the four 4 GHz wide IF bands coming from the two receiver sidebands and dual polarizations for an aggregate data rate of 64 Gbps.
For operations at high-altitude, helium-filled hermetically sealed hard drives are used, avoiding the need to build pressure chambers around the recorders. The hard drives collectively accumulated over 10,000 hr of recording at up to 5100 m altitude without a disk failure during the 2017 observations. The full 64 Gbps signal chain thus uses 128 hard drives in parallel, each with a capacity of 6–10 TB for a total storage capacity of about 1 PB. For the stations on the periphery of the network that cannot participate in all scans, 1 PB supports a typical VLBI observing campaign of about 6 days, while stations in the center need a double set of modules, or about 2 PB. Across the array, an EHT observing campaign involving all stations in Table 3 would require a total of about 15 PB of data storage.
Connected-element (sub)millimeter interferometers have more complex systems that use phased-array processors to sum signals from multiple telescopes into a single signal with much greater sensitivity. These processors are integrated with their digital correlators and output the summed signal as VDIF data packets to Mark 6 recorders as at single-dish stations.
The following sections describe each subsystem in more detail.
3.1. Receivers
The past three decades have seen the development and widespread use of heterodyne receivers in the millimeter and submillimeter bands based on superconductor–insulator–superconductor (SIS) junctions (e.g., Phillips et al. 1981 ; Maier 2009 ; Carter et al. 2012 ; Tong et al. 2013 ; Kerr et al. 2014). Over this period, instantaneous bandwidths increased by more than a factor of 30, while noise temperatures decreased by an order of magnitude. Improvements in receiver and antenna reflector technology have combined with the increased recording rates to lay the foundations for a millimeter wavelength VLBI array that is capable of observing targets with a flux density below 1 Jy. In the receivers, high-frequency radiation from the sky is mixed with a pure tone (an LO signal derived from a high-stability frequency reference) and downconverted to an IF using photon-assisted tunneling to transport current across the SIS junction’s energy gap. Presently installed SIS-based receivers at observatories typically support IF bandwidths of 4–12 GHz, with higher bandwidth systems becoming available at selected sites.
Heterodyne receivers have two observing bands with sky frequencies centered at $\nu _\mathrmLO-\nu _\mathrmIF$ (LSB) and $\nu _\mathrmLO+\nu _\mathrmIF$ (USB), respectively, from which the signals are folded on top of one another (double-sideband (DSB) receiver). With the use of an RF 90° hybrid, two mixers, and an IF 90° hybrid, modern millimeter receivers can separate the sidebands into two distinct IF channels (sideband-separating (2SB) receiver). The SIS mixers employed at these wavelengths use single-polarization rectangular waveguides that couple to a single polarization of the incident radiation through planar antennas mounted in rectangular waveguides. Dual-polarization receivers employ two independent orthogonal signal chains. Polarizations are split using a wire grid or a waveguide-based orthomode transducer (OMT). The latter can be inserted in front of the mixer blocks and cryogenically cooled, resulting in reduced ohmic losses. Some of the EHT stations use OMTs with circular polarizers, but the majority use quarter-wave plates in front of the receivers. ALMA records in dual orthogonal linear polarization, which is converted to circular polarization in post-processing as described in Section 5.
3.2. Hydrogen Maser Frequency Standards
As with connected-element arrays, VLBI relies on the fundamental ability of radio band receivers to accurately preserve the phase of detected cosmic radiation. For connected-element arrays, a common LO, derived from a station frequency reference, is used for all antenna receivers. On long VLBI baselines, sharing a common frequency reference is currently not technically feasible, and each VLBI site must use its own clock. This practice imposes a strict requirement on the stability of the frequency standards used for VLBI. Generally, one requires $\omega \tau \sigma _y(\tau )\ll 1$, where ω is the observing frequency in rad s−1, $\sigma _y(\tau )$ is the square root of the Allan variance (the Allan standard deviation), and τ is the integration time. This limit keeps rms phase fluctuations well below 1 rad (Rogers & Moran 1981).
Because of their excellent stability on timescales that match VLBI integrations (about 10 s, Figure 5), hydrogen masers are used almost without exception as VLBI frequency references worldwide. At $\sigma _y(10\,\rms)=1.5\times 10^-14$, VLBI observations at wavelengths down to 0.87 mm can be carried out with coherence losses ($\mathrmLoss\simeq 1-e^-(\omega \tau \sigma _y(\tau ))^2$) limited to $\lt 5$ % on a baseline where stations at either end are equipped with similarly stable masers (Doeleman et al. 2011). The use of hydrogen masers ensures that phase fluctuations in the VLBI signal due to atmospheric turbulence are typically the dominant source of coherence loss in EHT observations (Rogers et al. 1984 ; see also Figure 2).
Figure 5. Top panel : measured Allan standard deviation, $\sigma _y(\tau )$, of the hydrogen maser at the LMT compared to a precision quartz oscillator. The open red points are the manufacturer specifications of the hydrogen maser when compared to another maser. At a 1 s integration time, the quartz oscillator and maser have similar stability, so the measurements indicate that the maser is meeting its specifications as installed at a coherence time of 1 s. The flattening of the $\sigma _y(\tau )$ curve beyond 10 s is due to the decreased stability of the quartz crystal. The extrapolation of Allan deviation from short integration times to 10 s indicates that the maser meets the stability goal of $\sigma _y(10\,\rms)=1.5\times 10^-14$. Bottom panel : the long-term drift of the maser at the SMA compared to GPS, measured by differencing the 1 PPS ticks from the maser and local GPS receiver. The vertical width of the trace is due to variable ionospheric and tropospheric delays of the GPS signal (including the excursion near hour 200), while the long-term trend represents the frequency error of the maser. The drift measured from this plot, and its effects on the fringe visibility, are removed during VLBI correlation.
The hydrogen masers used in the EHT are monitored for stability using two tests. The first compares the stability of the 10 MHz frequency reference output of the maser with that of a high-precision quartz oscillator, known to have a similar performance as a maser on 1 s timescales. This in situ measurement confirms nominal maser operation at 1 s integration times, and extrapolation then yields an estimate of maser performance at 10 s (see Figure 5). A second test is to routinely monitor the offset between 1 PPS signals derived from the maser and a GPS receiver. GPS stability on 1 s timescales is poor due to variations induced by the ionosphere, but on longer timescales ($\gt 10^4$ s) GPS precision, which is referenced to the average of many atomic clocks, will exceed that of any single maser at an EHT location. An observed linear drift (Figure 5) over month to year timescales confirms long-term maser stability and also provides an instrumental delay-rate that is required for later interferometric detection (see Section 5).
3.3. Block Downconverters
A BDC design was developed to mix selected pieces of the incoming IF band from the receiver (in the range 4–9 GHz) down to baseband (0–2 GHz) for sampling. For observing at 1.3 mm, bandpass filters select 5–7 GHz and 7–9 GHz bands that are mixed against a LO at 7 GHz to convert to baseband. Similarly, for observing at 0.87 mm, 4–6 GHz and 6–8 GHz bands are mixed against an LO at 6 GHz to convert to baseband. The LO is phase locked to a 10 MHz reference from the station frequency standard. Coaxial relays in the output stage select which set of baseband outputs is connected to the data acquisition system, and the LO of the unused set is muted to avoid interference. The downconverter has two duplicate chains of filters and mixers for processing two of the four receiver IF channels, and two such downconverters are used in a complete system. A simplified schematic of a BDC is shown in Figure 6.
Figure 6. Simplified schematic diagram of the BDC. Only one of the two (identical) polarization channels is shown, and several intermediate amplification stages are omitted. Two local oscillators maser-locked and tuned to 6 GHz and 7 GHz are used to mix all the filtered IF bands between 4 and 9 GHz to baseband (0–2 GHz). For 230 GHz operation, a 5–9 GHz IF band is split with output of 5–7 GHz (’’low-band’’) and 7–9 GHz (’’high-band’’). For 345 GHz operation, a 4–8 GHz IF band is split with an output of 4–6 GHz (’’low-band’’) and 6–8 GHz (’’high-band’’). The 6 and 7 GHz local oscillator signals are also sent to the other polarization channel. The attenuation and band selection can be set using a control panel or remotely via an Ethernet connection.
The maximum conversion gain of each channel is about 23 dB, and it can be adjusted downward to about −8 dB in 0.5 dB steps to match the output power of the BDC with the input dynamic range of the digitizers. At a nominal output power of −7 dBm, the BDC operates in the linear regime over the entire programmable range of the attenuators. An 8-bit controller provides control of the synthesizers, digital attenuators, and coaxial relays, and interfaces to a keypad and display unit for manual operation, as well as to an Ethernet port for remote control.
3.4. Wideband VLBI Digital Backend
Coherent received station data in modern VLBI systems are recorded digitally. The VLBI instrument that digitizes and formats the analog received signal for recording is termed the digital backend (DBE). The high bandwidths that enhance the EHT sensitivity require proportionately fast digital sampling speeds. The timing of the samples is an implicit time stamping of the data, so the sampler clock must be timed with maser stability and precision.
For several generations of instrumentation the EHT has based its digital hardware on open-source technology shared by the Collaboration for Astronomy Signal Processing and Electronics Research (CASPER ;140 Hickish et al. 2016). The open-source hardware currently in use includes five gigasample-per-second (Gsps) analog-to-digital converter (ADC) boards, based on an integrated circuit that interleaves four cores, each with a maximum sample rate of 1.25 Gsps. The ADC circuit design and hardware testing is documented in Jiang et al. (2014). The compute hardware platform is the CASPER second-generation Reconfigurable Open Access Computing Hardware (ROACH2). The ROACH2 uses an FPGA as its digital signal processing engine.141
The current EHT DBE instrument was developed in 2014, and is called the ROACH2 Digital Back End (R2DBE ; Vertatschitsch et al. 2015). The R2DBE samples two analog channels at 4.096 Gsps, each with 8-bit resolution, i.e., two 2.048 GHz wide bands, before re-quantizing the samples with 2-bit resolution, packing the data in VDIF format, and then transmitting it over two 10 GbE connections to the Mark 6 recorder.142 Each R2DBE transmits two 8 Gbps data streams to a Mark 6 for recording. Accurate timing and synchronization is achieved by referencing the ROACH2 clocks to the maser via an external EHT-developed 2.048 GHz synthesizer.
The per-channel signal processing flow within the R2DBE is shown in Figure 7. Gain, offset, and sampler clock phase mismatches between the cores of the ADC produce spurious artifacts in the spectrum of the digitized signal, most prominently a spur at one quarter of the sample rate frequency that is caused by interleaving imperfection of the quad-core ADC. A calibration routine is performed at the start of each observation to tune the distribution of 8-bit samples for each core in offset, gain, and phase, removing these artifacts (Patel et al. 2014). A digital power meter is implemented in the firmware to calculate input power every millisecond, which is useful for pointing and calibration measurements.
Figure 7. Functional diagram of the R2DBE. A pair of ADCs sample two channels of 2.048 GHz Nyquist bandwidth IF bands, typically representing two antenna polarizations. After requantization of the samples to 2 bits, the data are distributed in VDIF packets over a 10 GbE network to Mark 6 recorders. The maximum FPGA clock speed is too slow to process the ADC output stream rate of 4096 megasamples per second in series, so 16 ADC samples are transferred from the sampler in parallel on each FPGA clock cycle and are processed in parallel though the FPGA, which is clocked at a sixteenth of the sample clock, or 256 MHz. The FPGA and sampler clocks are locked to the maser reference frequency. In addition to the sample clock and dual IF inputs, each ADC circuit board is also equipped with a low-frequency synchronization input. On one of the ADC boards, the synchronization input is connected to GPS PPS. An internal PPS for VDIF time stamping is synthesized from the 2.048 GHz maser-referenced clock and synchronized to the GPS real-time PPS time-tick once at the beginning of an observation. The internal PPS drift relative to GPS real-time PPS (Figure 5) is measured by a counter and is available to be read from a register over Ethernet, as well as copied into the VDIF data-header, as it is a prior parameter for VLBI correlation.
An important design consideration is the optimal number of bits per sample recorded onto storage media at the telescopes. This is determined by trading off the sensitivity increase realized through sampling more bandwidth at higher precision (see Equation (1)) with the cost of media required to store the data (for a detailed explanation, see Thompson et al. 2017). For the EHT, a maximum aggregate bandwidth of 16 GHz is set by the characteristics of ALMA receivers. When bandwidth is limited in this way, one examines the increase in sensitivity per additional unit media. For 1-bit recording (two-level sampling), the digital efficiency (ηQ in Equation (1)) is 0.64. Moving to 2-bit recording (four-level sampling) increases ηQ to 0.88 at the expense of doubling the required recording media, for a 38 % sensitivity increase per additional unit media—this is approximately equivalent to doubling the front-end bandwidth and keeping 1-bit sampling. Tripling the media cost and recording 3 bits (eight-level sampling) delivers only a 25.5 % increase in sensitivity per additional unit media. The 2-bit sampling scheme was chosen as it ensures ample margin for sensitivity requirements while minimizing recording media, which dominates the cost of VLBI correlation (Deller et al. 2011).
In a 2-bit system, the noise voltage thresholds of the four sampling levels should be properly set to maximize digital processing efficiency. A proper setting ensures that all four levels are optimally populated for a given input voltage signal from the telescope receiver system (Cooper 1970 ; Thompson et al. 2017). The processing efficiency is not a sensitive function of these thresholds, but the statistics of the level populations can vary significantly as elevation and weather condition changes cause fluctuations in the receiver power arriving at the sampler. For this reason, a calibration of the sampler threshold setting are performed periodically during calibration scans to maintain optimal settings.
3.5. High-speed Data Recorders
The Mark 6 VLBI Data System (Whitney et al. 2013) is a packet recording system used to store the R2DBE output streams to hard drives. The recorder captures the two 8 Gbps streams from the R2DBE using commercial 10 Gbps Ethernet network interface cards. It strips the internet packet headers and stores the payload containing VDIF data frames. For sustained recording at 16 Gbps, each Mark 6 recorder writes the data in time slices across 32 hard drives with a round-robin algorithm. The disks are mounted in groups of eight in four removable modules for ease of handling and shipping. For the EHT, four such recorders are configured in parallel to achieve an aggregate data capture rate of 64 Gbps.
Mark 6 recorders and modules are commercially available and in use in other VLBI applications. Several modifications were required for EHT use at high altitude, as hardware specifications typically state 10,000 ft as the maximum operating altitude. The low ambient air density necessitates sealed, helium-filled hard drives, both for the system disk for recorders and also in all modules onto which data are recorded. In addition, the Mark 6 interior was modified to direct a high-volume airflow onto the CPU, network, and data interface cards, which were found to be sensitive to overheating at altitude.
A photograph of the EHT VLBI backend (BDCs, R2DBEs, and Mark 6s) at the PV 30 m telescope is shown in Figure 8. Recorders used for earlier EHT campaigns included the Mark 5C system, which was capable of capturing only a 4 Gbps data stream onto two modules. The bandwidth of the EHT backend since 2004 is plotted in Figure 9 and represents a 64-fold increase, corresponding to an eight-fold improvement in sensitivity.
Figure 8. EHT digital VLBI backend as installed at the Institut de Radioastronomie Millimétrique (IRAM) PV 30 m telescope in Spain. The upper portion of the right-hand side rack holds the four R2DBE units. Two block downconverters are installed near the middle. The VLBI backend computer is mounted on the bottom right. The rack on the left and the lower portion of the rack on the right hold the four Mark 6 recorders with four disk modules each, providing a total of about 1 PB in data storage.
Figure 9. Recording rate of EHT observations over time. As of 2018, EHT stations record at 64 Gbps, equivalent to a doubling of recorded bandwidth every two years for over a decade (blue curve). The high bandwidth, a result of linking EHT instrumentation to industry trends and commodity electronics, is a crucial component of the EHT’s sensitivity, enabling detections on long baselines, providing resilience of the network against poor weather and low-elevation targets, and allowing detections from all stations to ALMA within short atmospheric coherence times (see Section 2.2).
3.6. Phased Arrays
To use the total collecting area available at connected-element interferometers that participate as single stations in the EHT VLBI network (ALMA, SMA, and, in the future, the Northern Extended Millimeter Array (NOEMA) these arrays coherently add the signal received from the target source by each antenna and record as if from a single antenna. Practical constraints on the maximum data rate that can be recorded at each site require that this summation be performed in real time.
Forming a coherent sum requires correcting for deterministic delays, such as geometric and known instrumental delays, as well as non-deterministic delays, which at EHT frequencies are significantly affected by the distribution of water vapor in the atmosphere. The atmospheric delay varies over time, antenna location, and direction so that accurate compensation can only be achieved through the use of in situ calibration methods.
Phasing systems were developed for the SMA and ALMA (see Appendix), and these observatories participated as phased arrays since 2011 and 2017, respectively. A phasing system for NOEMA is in the process of being implemented and commissioned for VLBI. Figure 10 shows the SMA phasing efficiency, i.e., the fraction of beamformed power compared to ideal phasing, estimated using Equation (2), achieved over the course of several scans during one night of the 2017 EHT campaign. For most of the scans, the efficiency is well above 0.9 (see inset). Typical phasing performance of the ALMA array can be found in Matthews et al. (2018).
Figure 10. Phasing efficiency of the SMA on various sources over the course of one observing schedule. The inset shows the cumulative distribution of phasing efficiency of the SMA over the entire 2017 EHT observing campaign.
3.7. Setup and Verification
As part of the integration of the new VLBI systems in 2015, most station backends were configured at Haystack Observatory prior to shipment to the telescopes. New equipment still typically passes through Haystack for initial inspection, check-out, and configuration. Each site takes responsibility for the installation of the equipment and ensuring that the connections to facility hardware work to specification and are not affected by factors such as telescope motion.
Before each observing campaign, each site goes through a comprehensive setup and verification procedure, which includes the completion of a checklist by site operational staff. The procedure verifies operation of the BDC, R2DBE, and Mark 6 systems, and the exact frequency and lock status of LOs, which set the exact observed sky band, and the ADC clock that time stamps the data. The check also includes the coherence and drift of the hydrogen maser time standard as this clock rate measurement is needed during correlation for fringe-stopping.
Due to the remote nature of many of its sites and the large data volume recorded, the EHT lacks real-time verification of fringes. Shipping data from remote stations to the central correlator and processing takes several days at minimum, and many months in the worst case (from the South Pole ; see Appendix A.11). If a key system fails, it is possible for a site to take data that never result in fringes, making careful testing and retesting of subsystems throughout the observation absolutely crucial. Data from brief observations (10–30 s) can be transferred from sites with fast internet connections for near real-time fringe verification. While possible, this requires robust data transfer connections.
The most complete in situ full-system check consists of the injection of a test tone at a known frequency in front of the receiver. A short data segment is recorded, and the autocorrelations and zero-baseline cross-correlations are inspected to verify that the test tone appears at the correct frequency and with the correct profile. Further, the tone in the baseband is mixed down to 10 kHz using a third LO and compared to a 10 kHz tone derived from the maser 10 MHz reference. The two should be phase locked when examined on an oscilloscope, which verifies the phase stability of the whole system. This test is exquisitely sensitive to small disturbances in the system. At the JCMT an equivalent test is done by verification of local connected-element fringes with the SMA.
Results and logs from the setup and verification are centrally archived and available to the correlator centers for the interpretation of station issues and for fault diagnosis when data quality issues emerge.
3.8. The Array
The EHT included eight observing facilities in 2017. Three additional facilities have since joined (GLT) or will soon join (KP 12 m and NOEMA), and two facilities (CARMA and the CSO) have participated in past EHT observations but have now been decommissioned. Properties of these facilities are summarized in Tables 2 and 3, and the Appendix. The baseline Fourier coverage provided by the array in 2017 for M87 is shown in Figure 11. Table 5 shows the actual performance of the array on scans of M87 during the 2017 science observations. The scan-averaged thermal noise within one 2 GHz frequency band was significantly better than 1 mJy on most baselines to ALMA, and generally a few mJy on baselines excluding ALMA. These achieved sensitivities, combined with the realized (u, v) coverage in Figure 11, confirm that the EHT has met its essential specification in providing an array that can image features on the scale of an SMBH shadow.
Figure 11. Aggregate EHT baseline coverage for M87 over four nights of observing with the 2017 array. Only detections are shown. The dashed circles show baseline lengths corresponding to fringe spacings of 25 and 50 μas. See Paper III for details.
Table 5. Median Thermal Noise (mJy) for Observations in 2017
Note. Median scan-averaged thermal noise per baseline in mJy for all M87 detections in 2017 April. Entries for the SPT reflect the median thermal noise on 3C 279, since the SPT cannot observe M87 due to its geographic location. Scan durations were three to seven minutes. For more details, see Paper III.
4. Observing
EHT science observing campaigns are scheduled for March or April when Sgr A* and M87 are night-time sources and the weather tends to be best on average over all sites. In addition, ALMA tends to be in a more compact configuration with better prospects for including more antennas in a phased array for increased sensitivity. Test and commissioning runs are scheduled a few months prior to campaigns to ensure that VLBI equipment, which may be dormant for months at a time, is operational.
4.1. Weather
Weather is an important consideration for VLBI observations at millimeter wavelengths. Most of the EHT observatories are located in the northern hemisphere, and those stations have especially large weather variations between seasons. Opacity and turbulence are typically lowest at night and in winter and early spring. At many sites, and particularly the connected-element arrays, the reduced atmospheric turbulence during night-time hours is essentially required. To protect against inclement weather, the EHT uses flexible observing with windows that are about twice as long as the intended number of observing nights. Within the window, a few hours before the start of observing each night, weather conditions are reviewed at all sites and a decision is made whether or not to observe that night. When an EHT night is not triggered, it is often possible for the time to be used for other observing programs at an observatory.
4.2. Scheduling
The process of scheduling starts with the list of approved target sources, time allocations at EHT telescopes and ALMA, and the dates of the observing window (Section 4.1). Schedule construction has to satisfy multiple requirements :
1.
a high total on-source time for each target,
2.
a wide range of parallactic angles sampled for polarimetry,
3.
long baseline tracks on all sources within the limited number of observing days,
4.
randomized scan lengths and scan start times on Sgr A* for periodicity analysis,
5.
good baseline coverage in the (u, v) plane, and
6.
regular gaps for telescope pointing and calibration.
AGN sources are chosen as calibrators based on their brightness, compactness, and proximity to the target sources. The calibration sources are used for multiple purposes, including fringe finding, bandpass calibration, and polarization calibration.
Observing blocks are created for each of the sources, and these blocks are merged into tracks, each corresponding to a single night of observing. Depending on time allocations and Fourier-coverage needs, multiple blocks on the same source may appear in the same track or in different tracks. Each observatory has different needs for overhead because the EHT is an inhomogeneous array. Overhead accommodates time for pointing, focus, and primary and secondary flux observations at all sites, plus phase-up time and array-polarimetric calibrations at the phased arrays. A typical EHT schedule (Figure 12) records VLBI scans with a duty cycle of approximately 50%, with a substantial fraction of those scans on targets strong enough to use for array calibration.
Figure 12. EHT 2017 observing schedule for M87 and 3C 279, covering one day of observations (April 6). Empty rectangles represent scans that were scheduled, but were not observed successfully due to weather, insufficient sensitivity, or technical issues. The filled rectangles represent scans corresponding to detections available in the final data set. Scan durations vary between 3 and 7 minutes, as reflected by the width of each rectangle, and scans are separated by periods of time used at each site for local pointing adjustments and calibration measurements. For this schedule, 3C 279 was the calibrator source for the primary target, M87. For information on other observing days, see Paper III.
4.3. Monitoring and VLBI Backend Control
The EHT developed a centralized monitoring tool called VLBImonitor to visualize observing status, collect ancillary calibration and weather data from each site, and provide real-time and predicted weather information from meteorological services. The information collected, logged, and displayed includes atmospheric and local weather conditions, observatory metadata such as telescope coordinates and on-source status, system temperature measurements, opacity measurements, digital backend and recorder state information, and comments from on-site observers and operators. Communication with the metadata server (by the software clients at each site or via a web interface) uses the JSON-RPC (remote procedure call) protocol over HTTPS. A ’’masterlist’’ defines metadata parameters that are accepted (white-listed) by the server and all their properties (e.g., data type, measurement cadence, a function that evaluates if the current value is valid or invalid, and units).
The VLBImonitor software143 was introduced for the 2017 observations. In 2018, a VLBI backend computer and network were added at the sites to provide a common monitoring and control platform at each station. Both monitoring and control remain in active development. At present, the EHT deploys specialist teams at each station, which, in practice, limits the observing window to about 12 days. It is a long-term objective for the EHT to both increase the length of this window and to conduct VLBI observations without the need for on-site specialists, after an initial setup and verification by local and remote experts. A critical component toward this goal is the implementation of comprehensive remote monitoring of the stations and VLBI equipment, as well as remote control of the VLBI backend.
5. Correlation and Calibration
The recorded data modules at all sites are separated by frequency band, with the ’’low-band’’ shipped to the VLBI correlator at MIT Haystack Observatory in Westford, Massachusetts, USA, and the ’’high-band’’ to the correlator at Max-Planck-Institut für Radioastronomie (MPIfR) in Bonn, Germany. Correlation is performed using the Distributed FX (DiFX) software correlator (Deller et al. 2011) running on clusters of more than 1000 compute cores at each site, and is split between the two sites to speed processing and allow cross-checks. At least as many Mark 6 playback units are needed at each correlator as there are stations in the EHT. The Mark 6 playback units at the MIT correlator are connected via 40 Gbps data links. A 100 Gbps network switch then delivers data to the processing nodes using 25 Gbps links. At MPIfR the internode communication, which includes the Mark 6 playback units, is realized via 56 Gbps connections, exceeding the maximum playback rate of the Mark 6 units of 16 Gbps.
Each 2 GHz observing band is correlated independently, with multiple passes required to correlate the full 4 or 8 GHz in an experiment. The correlation coefficients between pairs of antennas are calculated after correcting for an a priori clock model (Earth rotation, instrumental delays, and clock offsets and drift rates). All sites except ALMA record left and right circular polarizations (L/R) producing cross-correlations in the standard circular basis (LL, LR, RL, RR). ALMA antennas, however, are natively dual-linear polarization (X/Y), so a linear-to-circular polarization conversion is performed on ALMA baselines after VLBI correlation. PolConvert (Martí-Vidal et al. 2016) is a software routine developed for this purpose as part of the program to phase the ALMA array for VLBI operation (Matthews et al. 2018 ; Goddi et al. 2019). This routine forms linear combinations of the cross-polarization-basis correlator products (XR, YR) and (XL, YL) with ±90° phase shifts introduced to the visibility phase to produce complex visibilities in the circular-polarization basis, for each integration time. The input correlator products to PolConvert are equalized in gain and phase, and the X–Y channel phase differences are removed before polarization conversion. This X–Y channel equalization is performed using standard calibration techniques for connected-element ALMA operation (for details see Paper III).
To reduce data volumes and increase S/N, the cross-products computed in the correlator are averaged both in time and frequency. These complex cross-products are then fitted for fringe-delay (linear change in phase versus frequency) and fringe-rate (linear change in phase versus time), which are then removed. Residual fringe-delay and fringe-rate are due to several factors. Largest among these are residual clock drifts, instrumental electronic delays, and atmospheric phase fluctuations that can result in rapidly varying fringe-rate residuals. For the compact, narrow-field EHT continuum targets, delay and rate variability due to intrinsic source structure is expected to produce much smaller residuals. The averaging time and bandwidth in the correlator is thus set to ensure that any coherence losses due to delay or rate variations are negligible, or equivalently that such variations can be tracked both in time and frequency. For EHT observations, the typical averaging bandwidth is 0.5 MHz and the averaging time is 0.4 s, allowing residual delays up to ±1 μs and residual rates up to ±2.5 Hz. These settings enable the isolation and identification of instrumental spectral features, and also the ability to track atmospheric phase variations. This averaging in both time and frequency, results in a decimation of data volumes by more than a factor of 1 million.
After the initial correlation, the data are further processed through a pipeline that results in final data products for use in imaging, time-domain analyses, and modeling (Blackburn et al. 2019 ; Janssen et al. 2019). During fringe detection, the excess in correlated signal power due to the source on a VLBI baseline is identified, and the complex Fourier component of source brightness distribution is measured. As outlined in Section 2.2, high signal-to-noise detections on baselines to ALMA can be used to remove atmospheric phase fluctuations on non-ALMA baselines. Phase stabilizing the array in this way allows coherent integration of the VLBI signal on non-ALMA baselines beyond the atmospheric coherence time. Figure 13 demonstrates the process on real EHT data during relatively poor weather conditions, confirming the general approach and basis for the specifications in Section 2. Precise estimates of the observed systematic errors during the 2017 April EHT campaign are detailed in Paper III.
Figure 13. Fringes detected on EHT baselines before and after phase steering using the ALMA phased array as a reference station. Shown here are residual fringe-rate spectra, where the peak occurs at a single fringe-rate that is consistent with linear clock drift over the entire scan. Stochastic atmospheric phase variations result in fringe-rate variation over time, spreading the signal in the plot on the left. To recover the total fringe amplitude, baselines to ALMA are used to phase reference the array at a short timescale. Data shown are from a VLBI scan on a quasar calibrator (3C 279) obtained during EHT observations in 2017 April. The data are taken at low elevation, and at several sites during the afternoon and early evening when the atmosphere is often unstable. Baselines to ALMA are able to phase steer the remaining sites in the array at an effective timescale of 1–3 s depending on S/N. At SMA due to the mid-afternoon local observation, phasing is particularly challenging. This leads to an amplitude efficiency of about 90 % on SMA baselines, as well as increased measurement thermal noise for the APEX-SMA spectrum with a full-scan S/N of only 10.
Subsequent calibration steps convert the correlation coefficients derived from fringe detection to correlated flux density (Jy). This is accomplished through use of a priori information about station sensitivities, and by the application of self-calibration techniques to correct for variations in telescope gain over time and frequency. In cases where two EHT telescopes are in close geographical proximity (e.g., ALMA-APEX and SMA-JCMT), additional calibration constraints can be derived from the resulting baseline redundancy and using only general assumptions about the observed source (Paper III).
6. Future Developments
The EHT array is continuing to develop. With the ability to image SMBHs on horizon scales now confirmed (Papers I ; III ; III ; V ; VI), the focus of EHT development will shift to enabling observations that can refine constraints on fundamental black hole properties, processes of black hole accretion and outflow, and tests of general relativity. This depends on achieving higher angular resolution, enhancing image fidelity, and enabling dynamic imaging of time-variable phenomena. Higher resolution will allow more detailed studies and modeling of sub-horizon structures as well as sensitive tests for asymmetries in shadow features. Greater image fidelity will bring fainter emission near the horizon into focus for the study of accretion and jet processes, and it will enable a sensitive comparison across imaging epochs, which is especially germane for M87 with its dynamical timescale of days to weeks. For Sgr A*, the light-crossing time of 20 s requires a dynamic approach to image reconstruction, with the potential of observing the near real-time evolution of a black hole.
The planned and in-progress addition of telescope facilities at new geographic sites will improve the (u, v), or Fourier, coverage, and thus imaging fidelity. Over the course of the next two years, the EHT expects to add two more facilities : a beamformed NOEMA in France, which, once completed, will be the equivalent of an approximately 50 m dish, and a 12 m diameter dish on Kitt Peak in Arizona. A newly realized and important consequence of designing high-bandwidth systems is that adding telescopes with modest diameters (sime6 m) creates VLBI baselines with sufficient sensitivity to detect the primary EHT targets. An expansion of the EHT can, therefore, include the possibility of deploying numerous smaller apertures, enabling not only improved image quality but also snapshot capability, which is a precursor to constructing black hole movies (Bouman et al. 2018 ; Johnson et al. 2018). Additional plans to enhance the data throughput of VLBI backend and phased array systems are underway, linked to new generations of wideband millimeter and submillimeter receivers and industry trends that will allow for rapid design and implementation. With the inclusion of more sites and wider bandwidths, the computational requirements for VLBI correlation will increase : linearly with bandwidth, and as the square of the number of stations. The exploration of scalable approaches to address future EHT correlation needs is underway (Gill et al. 2019). The combination of these efforts is aimed at the substantial expansion of aperture plane sampling of the EHT and improved imaging.
To significantly improve angular resolution requires development in different directions. Planned extension of the EHT to operation at 0.87 mm wavelength, a standard receiving band at many facilities, would increase the angular resolution of the array by 40%. VLBI tests at 0.87 mm are underway and this capability is expected at a subset of EHT sites over the next 3–5 yr (see Table 4). An alternate approach to increased angular resolution is to deploy EHT antennas in space where baseline length is not limited to the diameter of the Earth, potentially allowing horizon-scale imaging of additional SMBH candidates (Johannsen et al. 2012). Space platforms also offer the possibility of rapidly sampling the Fourier plane, thereby opening the potential for dynamical imaging of black hole accretion and outflow processes (D. Palumbo et al. 2019, in preparation ; F. Roelofs et al. 2019, in preparation ; M. Shea et al. 2019, in preparation).
EHT observation strategies continue to be refined. With further development of remote monitoring and control tools, the EHT can explore triggering observations during the best conditions throughout the year outside the current March–April window. One consequence of this flexibility would be that observations for Sgr A* and M87 could be optimized separately in a given observing cycle, instead of grouped together as they are currently. Distributing the observations over a larger portion of the year also increases the likelihood that the EHT would detect emission transients or flaring on intermediate timescales should they occur, and affords, in general, the opportunity to study M87 on timescales that correspond to the expected ISCO period of order one to several weeks (Table 1). Not all EHT sites may be able to participate in such flexible campaigns, but even a subset of the array, especially if augmented with many smaller dishes, could provide useful observations for variability studies.
7. Conclusion
The goal of the EHT project is to observe SMBHs with spatial and temporal resolution that permits imaging on the scale of the lensed photon orbit, and the study of dynamics on commensurate light-crossing timescales. Leading up to the detection of event-horizon-scale structure in Sgr A* in 2008, and now a decade afterward, the development of EHT instrumentation proceeded through a series of systems that supported increasingly ambitious observations. These observing campaigns led to precursor scientific results that motivated a strategy of building an imaging array : one with baselines long enough to resolve horizon-scale structure, that included enough geographical sites to sufficiently sample the Fourier plane, and had the sensitivity required to detect compact, low flux density targets. Technology maturation, fueled by industry-driven trends, was focused primarily on improving the bandwidth of observations by increasing digitization and recording rates. Over a period of 10 yr, EHT data capture rates increased from 4 to 64 Gbps, allowing for observations that take full advantage of the capabilities of modern millimeter receivers. The parallel development of phased array technologies and investment in site-specific infrastructure (see the Appendix) have now led to a global EHT array with the capability to address its science goals.
Had the EHT project needed to build a VLBI array from the ground up, the cost would have been prohibitive. By developing new systems and infrastructure that allowed for the use of existing telescopes and facilities, most of which were not originally conceived to operate as VLBI elements, the EHT has created a purpose-built array with unique capability at modest cost.
Appendix : Specific Details of Participating Facilities in EHT Observations
One of the technical challanges of the EHT array is that it is assembled from existing telescope facilities that differ to a greater or lesser extent in their technical characteristics and readiness for VLBI operations. As a result, EHT instrumentation installed at each site requires customization. This appendix provides information on the specifics of each site in turn, not general to the array overall.
A.1. ALMA (5100 m Altitude)
The ALMA, located on the Chajnantor plain in Chile, is an international partnership between Europe, the United States, Canada, Japan, South Korea, Taiwan, and Chile (Wootten & Thompson 2009). The connected-element interferometer consists of fifty 12 m and twelve 7 m antennas, and is supplemented with four 12 m total-power antennas. Reflector surface accuracy is better than 25 μm rms.
The ALMA Phasing Project (APP) was an international effort to produce an antenna phasing system leveraging the superb (sub)millimeter capabilities and large collecting area of ALMA for VLBI applications (Doeleman 2010 ; Matthews et al. 2018). The APP added to the hardware and software already available at the ALMA site for connected-element interferometry the additional components necessary for coherently summing the ALMA antennas and recording VLBI data products. The rubidium clock that was formerly used to generate the reference 5 MHz LO tone was replaced with a hydrogen maser. Phasing interface cards were added to the ALMA Baseline Correlator to serve as the VLBI backend. A fiber link system was built and deployed to transport the phased-sum signal from the ALMA high-elevation site to the ALMA Operations Support Facility at an elevation of about 2900 m, where a set of Mark 6 VLBI recorders was installed. Numerous software enhancements were also required, including the implementation of an ALMA VLBI Observing Mode (VOM) and a phase solver to calculate the phases needed to adjust each of the ALMA antennas to allow coherent summation of their signals. The ALMA phasing system treats the phased sum as another input to the Baseline Correlator, so the phased-sum signal can be correlated with individual antenna signals. Standard ALMA single-field interferometric data are generated as a matter of course during VOM operation, and ALMA-only data products are archived as usual. When all antennas are included in the phased sum, ALMA is equivalent to an effective 88 m diameter aperture for VLBI operations. Details are described in Matthews et al. (2018) and Goddi et al. (2019).
A.2. APEX (5100 m Altitude)
APEX is a 12 m antenna located on the Chajnantor plain in Chile about 2 km from the ALMA site (Güsten et al. 2006). It was built by VERTEX Antennentechnik from the ALMA prototype development and inaugurated in 2005. It was outfitted and operated by the Max-Planck-Institut für Radioastronomie in Bonn (50%), the Onsala Space Observatory (23%), and the European Southern Observatory (27%) as the first telescope on the Chajnantor plain. The first VLBI fringes were obtained with APEX in 2012 at 230 GHz (Wagner et al. 2015).
The reflector surface accuracy is 17 μm rms, which gives good efficiency at frequencies exceeding 1 THz, and the telescope is correspondingly equipped with receivers spanning 200 GHz–1.4 THz. The surface panels and subreflector were replaced in early 2018 to improve the surface to 10 μm rms and so further raise the efficiency at the highest frequencies. At 230 GHz the elevation gain curve is flat to about 3% due to the high structural rigidity of the antenna. Pointing accuracy is better than 2 arcsec rms, and tracking accuracy is typically below 1 arcsec rms depending on wind conditions. The beam FWHM is 27.1 arcsec at 228.1 GHz.
The maser 10 MHz reference is transported to the front-end synthesizers and backend equipment on four low temperature-coefficient, double-shielded 85 m length cables. The cables are wrapped in thermal insulation and firmly fixed to support structures in the receiver cabin to minimize phase instabilities resulting from temperature changes or telescope movement. No round-trip phase stabilization is used. The front-end synthesizers (first LO and tone generation) each have their own frequency reference cable from the maser. References for the backend equipment are supplied through a distributor because the frequencies generated are much lower than for the front-end. The R2DBE and another VLBI backend called the DBBC are located in the antenna close to the receiver (15 m cable length) for low loss transmission of the dual-polarization dual-sideband 4–12 GHz receiver IF output. The Mark 6 recorders are located in a control container about 50 m from the telescope. The backends transmit eight 8 Gbps VDIF-format data streams on fiber with short-range 10 GbE 850 nm SFP+ fiber transceivers to the recorders.
A.3. GLT (100 m Altitude Currently, 3200 m Altitude Planned)
The GLT (Inoue et al. 2014) was originally the 12 m ALMA-North America prototype antenna (Mangum et al. 2006) located at the Very Large Array (VLA) site near Socorro, New Mexico. The antenna was awarded to the Smithsonian Astrophysical Observatory (SAO) in 2011, representing the partnership of the SAO and Academia Sinica Institute of Astronomy (ASIAA). The antenna was retrofitted and rebuilt for operation in the extreme arctic conditions of northern Greenland. Presently located at a temporary site on Thule Air Base, Greenland, the telescope is equipped with a new set of submillimeter receivers operating at 86, 230, and 345 GHz. In 2018, the GLT dish surface underwent iterative adjustments to improve surface accuracy, eventually achieving better than 40 μm rms. The pointing accuracy is about 3−5 arcsec during 230 GHz observations. The telescope is planned to be relocated to Greenland Summit, where the atmosphere provides excellent transparency for submm and THz observations (Matsushita et al. 2017).
The GLT has a complete VLBI backend. In 2018 January, fringes were detected between the GLT and ALMA, and in 2018 April, the GLT participated in EHT observations (Chen et al. 2018). The maser and associated equipment are located in a dedicated trailer approximately 40 m from the telescope. The maser 10 MHz is fed through a phase-stable coaxial cable to an adjacent synthesizer that generates the first LO (18–31 GHz) reference for the receivers. This first LO reference along with a separate maser reference of 100 MHz are combined and transmitted to the antenna through a single-mode fiber and used to phase-lock all the associated instruments in the receiver cabin (Kubo et al. 2018). The BDCs and R2DBEs are located in the receiver cabin and are interconnected using short segments of phase-stable coaxial cables. The Mark 6 recorders are located in the VLBI trailer approximately 35 m from the telescope. Communications between the R2DBEs and the Mark 6 recorders are via a 10 GbE optical network. The current system supports a throughput of 64 Gbps to the recorders.
A.4. JCMT (4100 m Altitude)
The JCMT is located on Maunakea in Hawai’i. It was dedicated in 1987 and was operated by a consortium (the United Kingdom, Canada, and the Netherlands) until 2015. Since then, the East Asia Observatory has operated the JCMT with funding from Taiwan, China, Japan, Korea, the United Kingdom, and Canada. The reflector is a 15 m Cassegrain design and has a surface accuracy of about 25 μm rms. The telescope is enclosed in a co-rotating dome and is protected by a 0.3 mm thick GoreTexTM windblind that is not removed for observing. The atmospheric opacity is measured with a line-of-sight 183 GHz radiometer.
The telescope participated in the VLBI observations that in 2007 detected the horizon-scale emission in Sgr A* (Doeleman et al. 2008) and has successfully participated in all EHT observing since then. This was done in cooperation with the SMA, which generated the phase-stable frequency reference signals and recorded the IF, utilizing the extended SMA (eSMA) infrastructure that enables the JCMT to join the SMA array as a connected element. From the 2019 VLBI campaign the JCMT will use a locally installed VLBI backend for recording, and will be equipped with a new receiver outfitted with ALMA 230 and 345 GHz mixers at 1.3 and 0.87 mm built by ASIAA.
A.5. Institut de Radioastronomie Millimétrique : 30 m (PV, 2900 m Altitude)
The IRAM 30 m telescope is located on Pico Veleta in the Spanish Sierra Nevada near Granada. The telescope was built in the early 1980s (Baars et al. 1987). The first VLBI fringes at 1.3 mm were obtained in 1995 (Greve et al. 1995) with one antenna of the Plateau de Bure interferometer, and since 1996 the PV 30 m has participated in the regular 3 mm VLBI CMVA/GMVA biannual runs.
The telescope surface has been adjusted to an accuracy of about $55\,\mu \rmm$ rms using phase-coherent holography with geostationary satellites. This yields a 230 GHz aperture efficiency of 47% at the optimal elevation of 50°. The gain elevation curve drops to about 35 % at elevations of 20° and 80°. The blind pointing accuracy is about 3 arcsec. The tracking accuracy is better than 1 arcsec for wind speeds below about 10 m s−1. Active thermal control of the telescope backstructure allows for efficient operation around the clock, and a de-icing system enables quick startup after winter storms. A 225 GHz radiometer continuously monitors the sky opacity at a fixed azimuth. Since 2018 October, the 30 m telescope is connected via a 1 Gbps fiber link to IRAM/Granada and the Spanish research network RedIRIS.
The PV 30 m is currently equipped with two heterodyne receivers. One of these, the Eight MIxer Receiver (EMIR) installed in 2009, has eight sideband-separating (2SB) mixers. The EMIR has 8 GHz of bandwidth per sideband covering the atmospheric windows at 3, 2, 1.3, and 0.87 mm wavelengths in dual-polarization. The EHT uses the EMIR 1.3 and 0.87 mm bands. Conversion from linear to circular polarization is done by quarter-wave plates, which are installed in a carousel to quickly switch between non-VLBI and VLBI observations. The IF signal from the receiver cabin is transported to a temperature-controlled backend room over 100 m coaxial cables. A new hydrogen maser was installed in 2018 January and then used for the GMVA and EHT sessions in 2018 April. For the EHT, four R2DBEs deliver 64 Gbps (two polarizations, two sidebands) to four Mark 6 recorders (see Figure 8).
A.6. IRAM : NOEMA (2600 m Altitude)
NOEMA is located in the French Alps and is the largest millimeter-wave facility in the northern hemisphere. With twelve 15 meter antennas, 2SB receivers with 8 GHz IF per sideband and polarization, and baselines up to 1.7 km, NOEMA will operate at 10 times the sensitivity and 4 times the spatial resolution of its predecessor, the former Plateau de Bure Interferometer (Guilloteau et al. 1992).
The NOEMA Phase 1 extension program was completed in 2018 September with the delivery of four additional antennas to the existing six, each with a surface accuracy of 35 μm rms. The dual-polarization 2SB receivers operate in the 3, 2, and 1.3 mm wavelength bands, with noise temperatures of, respectively, 35 K, 40 K, and 50 K in each band. Signal processing is done using an FPGA-based FFX correlator with flexible configurations modes. Beginning in 2019, NOEMA will have a phased-array mode with an initial processing bandwidth of 2 × 2 GHz (16 Gbps), eventually growing to 128 Gbps. The NOEMA Phase 2 extension program (2019) includes the construction of two additional antennas, the delivery of a second correlator for dual-band operation (2 × 32 GHz), a receiver upgrade to perform full-array observations up to 373 GHz (0.8 mm), and a baseline extension to 1.7 km for high-fidelity imaging down to a spatial resolution of 0.1 arcsec.
A.7. Kitt Peak 12 m (KP, 1900 m Altitude)
The telescope on Kitt Peak was acquired by the University of Arizona in 2013 March. The 12 m dish was a prototype antenna (Mangum et al. 2006) made for the ALMA project by the Alcatel/European Industrial Engineering consortium. It was relocated from the VLA site in New Mexico to the enclosure of the previous NRAO 12 m telescope on Kitt Peak, Arizona, in 2013 November. The telescope is operated by the Arizona Radio Observatory (ARO), a division of the Steward Observatory of the University of Arizona. First light was observed in 2014 October. The surface was adjusted to an rms figure error of 16 μm in 2018 September using photogrammetry. A 225 GHz radiometer measures atmospheric opacity.
VLBI observations will make use of a newly developed multi-band receiver (1.3–4 mm), which incorporates ALMA Band 3 and Band 6 mixers for the 3 and 1.3 mm bands. All bands sample both polarizations, and the 1.3 mm band will be converted to circular polarizations through a quartz quarter-wave plate. The receiver is under commissioning and the telescope is scheduled to participate in the next scheduled EHT observations. The master reference for this site is located within the domed enclosure adjacent to the computer room that houses the VLBI backend electronics.
A.8. LMT (4600 m Altitude)
The LMT is situated at the summit of Volcán Sierra Negra in Central Mexico. It is jointly operated by the Instituto Nacional de Astrofísica, Óptica y Electrónica (INAOE) and the University of Massachusetts (UMass). In 2017, two additional rings of precision surface panels were installed to enlarge the LMT primary mirror from 32.5 to 50 m. The primary mirror is an active surface comprised of 180 segments steered by 720 actuators to correct thermal and mechanical distortions. The large collecting area and a central geographical location with respect to other EHT sites make the LMT particularly important for array sensitivity and imaging fidelity.
Initial technical work to develop 3 mm wavelength VLBI at the LMT was the product of a multi-year collaboration between UMass, INAOE, SAO, the Massachusetts Institute of Technology Haystack Observatory, the Universidad Nacional Autónoma de México (UNAM), and NRAO. The first 3 mm observations were carried out in 2013 April with the facility Redshift Search Receiver (Erickson et al. 2007), a broadband dual-polarization, dual-beam receiver that acted as the front end. Those observations were performed with a VLBI backend using a high-precision quartz oscillator frequency reference running at 2 Gbps recording rate. Subsequent integration of a hydrogen maser and further testing in 2014 led to first VLBI science observations at 3 mm wavelength in 2015 (Ortiz-León et al. 2016).
In 2014, an interim special-purpose VLBI 1.3 mm wavelength receiver was assembled using mixers from the CARMA observatory and resourced through a collaboration including the institutes above and the University of California, Berkeley. This receiver was a single-pixel, dual-polarization, double-sideband receiver. It was deployed in 2015 for commissioning along with a 32 Gbps VLBI system (Figure 4), both of which were used for the 2017 April EHT observations. In early 2018, a new dual-polarization, sideband-separating 230 GHz receiver, using ALMA mixers, was delivered by UMass and installed along with an upgrade of the VLBI system to 64 Gbps capability.
The hydrogen maser at the LMT is located in a temperature-stable environment deep in the apex cone. The maser reference is brought to a 10 MHz distribution system in the VLBI backend room. The backend room is located one level below the receiver cabin and houses the 2.048 GHz synthesizer, the 1 PPS distribution system, BDCs, R2DBEs, and the Mark 6 recorders.
A.9. SMA (4100 m Altitude)
The SMA is an eight-element radio interferometer located atop Maunakea operated by SAO and ASIAA (Ho et al. 2004). The 6 m dishes are configurable with baselines as long as 509 m, which produces a synthesized beam of about 1 arcsec width at 230 GHz, and sub-arcsec width at 345 GHz. The SMA has participated in EHT observations since 2006. In 2006 and 2007, the SMA provided a maser reference tuning signal to the Caltech Submillimeter Observatory and the JCMT, respectively. In 2009, the SMA contributed collecting area to the EHT through the Phased Array Recording Instrument for Galactic Event Horizon Studies (PhRInGES), a 4 Gbps phased-array system.
During EHT observations, the SMA uses its own correlator, known as the SMA Wideband Astronomical ROACH2 Machine (SWARM), instead of the R2DBE. SWARM is a combined correlator and phased array that was initially deployed in 2014 and built out to its full 32 GHz bandwidth by 2017 (Primiani et al. 2016). SWARM supports VLBI through a built-in beamformer based on a signal-noise decomposition of the correlator output while observing the target source (Young et al. 2016). Phasing the array requires tracking all sources of delay, including fluctuations in water vapor concentration in the atmosphere. The SWARM phasing system is equipped with a real-time phasing solver that continually updates the beamforming weights over the course of the observation to compensate for these variable delays, which manifest as variable phase errors in each antenna. Because the phased array capability is used to observe sources that are unresolved on baselines within the array, the corrective beamformer weights can be computed by extracting from the correlator output that contribution associated with a point-like source. Furthermore, as the weights are applied to the signal from each antenna before computing cross-correlations between antenna pairs, the solution obtained from the correlator output for a particular integration period can also be used to calculate the average phasing efficiency over that same period. In addition to solving and applying corrective phasing in real time, the system produces an estimate of the phasing efficiency, ηphgr,
Equation (2)
where wi is the complex-valued weight applied to antenna $\rmi$. See Young et al. (2016)144 for a more detailed discussion of the phased array and performance assessment thereof.
While the SWARM correlator architecture lends itself conveniently to phased-array processing, it also presents several challenges for integrating the instrument into a VLBI array. First, the beamformer processing of 2 GHz bands is spread across eight ROACH2 devices. These parallel data streams must be collected and formatted in real time in order to record to the Mark 6, which is done using a custom FPGA device called the SWARM Digital Backend (SDBE). The data are received on four of the eight 10 GbE ports on the SDBE. The packets are time-stamped, the frequency domain samples are quantized from 4-bit complex to 2-bit complex, the packets are formatted with VDIF headers, and then transported over Ethernet to the Mark 6. As the beamformer data packets are relatively small, several of these packets are bundled into each UDP packet to reduce the interrupt rate on the Mark 6 and prevent packet loss. The second complication to using SWARM in VLBI is that the data are in complex frequency domain representation as opposed to the time-domain format used by the majority of stations in the EHT, and the data are sampled at a different rate than the single-dish EHT stations : 4.576 Gsps compared to 4.096 Gsps at single-dish sites. Prior to the correlation of SWARM data with data from other sites, an offline inversion and resampling program known as the Adaptive Phased-array and Heterogeneous Interpolating Downsampler for SWARM (APHIDS) is used to preprocess the SMA data. Development of a real-time preprocessing system to be integrated into the SDBE is currently in progress. The SMA VLBI pipeline is illustrated in Figure 14.
Figure 14. Block diagram demonstrating the VLBI data pipeline, from SWARM to correlatable data in VDIF format. The SDBE is integrated with SWARM and does the real-time processing necessary to interface with the on-site Mark 6 during an observation. After observing, the data is preprocessed offline in APHIDS prior to correlation with data from other sites.
A.10. SMT (3200 m Altitude)
The SMT, located on Mount Graham, Arizona, was built by the Max-Planck-Institut für Radioastronomie in 1993, and is now operated and maintained by the ARO (Baars et al. 1999). Work to incorporate the telescope in mm-wavelength VLBI experiments began in the late 1990s (Doeleman et al. 2002 ; Krichbaum et al. 2002). In 2007, the SMT participated in the observations that detected horizon-scale emission in Sgr A* (Doeleman et al. 2008). The 10 m telescope is mounted in a co-rotating enclosure that minimizes the effects of insolation and wind during observations. The telescope surface accuracy is 15 μm rms and the telescope pointing is typically accurate to better than 1.5 arcsec rms. Due to weak diurnal variations in water vapor and the protection of the telescope enclosure, the SMT can observe 24 hr per day. Atmospheric opacity is monitored at this site with a tipping 225 GHz radiometer that is mounted to the rotating dome.
The SMT is equipped with dual-polarization heterodyne receivers covering the atmospheric windows from 200 to 700 GHz. For VLBI observations, the 1.3 mm receiver system uses prototype ALMA sideband-separating mixers, with polarization splitting achieved via a room-temperature wire grid, and circular polarization induced through a quartz quarter-wave plate. The SMT also has a double-sideband receiver for observations at 0.87 mm ; a replacement sideband-separating receiver for this band is in development at the University of Arizona. The 64 Gbps EHT backend hardware is mounted in the receiver cabin. The master reference for the SMT is a hydrogen maser mounted at the base of the telescope pier in a temperature-stabilized closet.
A.11. SPT (2800 m Altitude)
The SPT is a 10 m diameter mm-wavelength telescope located at the geographic South Pole (Carlstrom et al. 2011) and is operated by a collaboration led by the University of Chicago. The off-axis Gregorian design is optimized for observations of the cosmic microwave background with wide-field, multi-color bolometer cameras. The primary reflector surface accuracy is 20 μm rms.
EHT observations with the SPT are managed by the University of Arizona. The standard observing campaigns occur during the austral fall when the South Pole is inaccessible. VLBI setup and testing is performed by a two-person SPT winter-over team, with remote supervision from Arizona. EHT observing at the SPT employs a dual-frequency 230/345 GHz VLBI receiver (Kim et al. 2018a), which was constructed at the University of Arizona. The receiver was first installed in late 2014 and first operated in 2015 January (Kim et al. 2018b). The receiver was reinstalled in an upgraded receiver cabin in 2016 in preparation for the 2017 EHT campaign. To illuminate the VLBI receiver, an alternate secondary mirror is installed ahead of prime focus to form a Cassegrain telescope, which relays the light to an ellipsoidal tertiary mirror atop the receiver. The full VLBI signal chain includes a hydrogen maser and coherence monitoring loops to account for the unusual range of environmental temperatures at the site. Mark 6 data modules are shipped after flights resume in the austral spring season, around November 1.
M87-3 Data Processing and Calibration
Footnotes
136 The bandwidth of an EHT observation is typically expressed as a recording rate in gigabit-per-second, or Gbps. The recording rate in Gbps is four times the recorded bandwidth in GHz, a factor of two coming from the need to take samples at a rate of twice the bandwidth (Nyquist rate), and another factor of two because each sample is 2 bits.
137 The SEFD of a radio telescope is the total system noise represented in units of equivalent incident flux density above the atmosphere (see Paper III, Equation (3)).
138 The visibility is the complex two-point correlation of the electric field measured on a VLBI baseline (see Thompson et al. 2017, Chapter 1).
140 For more information on CASPER, please see https://casper.berkeley.edu/.
141 See https://www.xilinx.com/support/documentation/data_sheets/ds150.pdf.
142 An alternative DBE that is planned to be deployed at some EHT sites is the third-generation Digital BaseBand Converter (DBBC3 ; Tuccari et al. 2014).
143 https://bitbucket.org/vlbi
144 See also the related SMA Memo 163 at https://www.cfa.harvard.edu/sma/memos/163.pdf.