National Radio Astronomy Observatory

The all sky surveys done by the Palomar Observatory Schmidt, the European Southern Observatory Schmidt, and
\nthe United Kingdom Schmidt, the InfraRed Astronomical Satellite, and the Two Micron All Sky Survey have
\nproven to be extremely useful tools for astronomy with value that lasts for decades. The Wide-field Infrared
\nSurvey Explorer (WISE) is mapping the whole sky following its launch on 2009 December 14. WISE began
\nsurveying the sky on 2010 January 14 and completed its first full coverage of the sky on July 17. The survey
\nwill continue to cover the sky a second time until the cryogen is exhausted (anticipated in 2010 November).
\nWISE is achieving 5σ point source sensitivities better than 0.08, 0.11, 1, and 6 mJy in unconfused regions on
\nthe ecliptic in bands centered at wavelengths of 3.4, 4.6, 12, and 22μm. Sensitivity improves toward the ecliptic
\npoles due to denser coverage and lower zodiacal background. The angular resolution is 6".1, 6".4, 6".5, and 12".0 at 3.4, 4.6, 12, and 22μm, and the astrometric precision for high signal-to-noise sources is better than 0".15.
\n

Abstract The Astropy Project supports and fosters the development of open-source and openly developed Python packages that provide commonly needed functionality to the astronomical community. A key element of the Astropy Project is the core package astropy , which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package, as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of interoperable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy Project.

The NRAO VLA Sky Survey (NVSS) covers the sky north of J2000.0 δ = -40° (82% of the celestial sphere) at 1.4 GHz. The principal data products are (1) a set of 2326 4° × 4° continuum "cubes" with three planes containing Stokes I, Q, and U images plus (2) a catalog of almost 2 × 106 discrete sources stronger than S ≈ 2.5 mJy. The images all have θ = 45'' FWHM resolution and nearly uniform sensitivity. Their rms brightness fluctuations are σ ≈ 0.45 mJy beam-1 ≈ 0.14 K (Stokes I) and σ ≈ 0.29 mJy beam-1 ≈ 0.09 K (Stokes Q and U). The rms uncertainties in right ascension and declination vary from ≲1'' for the N ≈ 4 × 105 sources stronger than 15 mJy to 7'' at the survey limit. The NVSS was made as a service to the astronomical community. All data products, user software, and updates are being released via the World Wide Web as soon as they are produced and verified.

Full list of authors: Price-Whelan, Adrian M.; Lim, Pey Lian; Earl, Nicholas; Starkman, Nathaniel; Bradley, Larry; Shupe, David L.; Patil, Aarya A.; Corrales, Lia; Brasseur, C. E.; Noethe, Maximilian; Donath, Axel; Tollerud, Erik; Morris, Brett M.; Ginsburg, Adam; Vaher, Eero; Weaver, Benjamin A.; Tocknell, James; Jamieson, William; van Kerkwijk, Marten H.; Robitaille, Thomas P.; Merry, Bruce; Bachetti, Matteo; Gunther, H. Moritz; Aldcroft, Thomas L.; Alvarado-Montes, Jaime A.; Archibald, Anne M.; Bodi, Attila; Bapat, Shreyas; Barentsen, Geert; Bazan, Juanjo; Biswas, Manish; Boquien, Mederic; Burke, D. J.; Cara, Daria; Cara, Mihai; Conroy, Kyle E.; Conseil, Simon; Craig, Matthew W.; Cross, Robert M.; Cruz, Kelle L.; D'Eugenio, Francesco; Dencheva, Nadia; Devillepoix, Hadrien A. R.; Dietrich, Jorg P.; Eigenbrot, Arthur Davis; Erben, Thomas; Ferreira, Leonardo; Foreman-Mackey, Daniel; Fox, Ryan; Freij, Nabil; Garg, Suyog; Geda, Robel; Glattly, Lauren; Gondhalekar, Yash; Gordon, Karl D.; Grant, David; Greenfield, Perry; Groener, Austen M.; Guest, Steve; Gurovich, Sebastian; Handberg, Rasmus; Hart, Akeem; Hatfield-Dodds, Zac; Homeier, Derek; Hosseinzadeh, Griffin; Jenness, Tim; Jones, Craig K.; Joseph, Prajwel; Kalmbach, J. Bryce; Karamehmetoglu, Emir; Kaluszynski, Mikolaj; Kelley, Michael S. P.; Kern, Nicholas; Kerzendorf, Wolfgang E.; Koch, Eric W.; Kulumani, Shankar; Lee, Antony; Ly, Chun; Ma, Zhiyuan; MacBride, Conor; Maljaars, Jakob M.; Muna, Demitri; Murphy, N. A.; Norman, Henrik; O'Steen, Richard; Oman, Kyle A.; Pacifici, Camilla; Pascual, Sergio; Pascual-Granado, J.; Patil, Rohit R.; Perren, Gabriel, I; Pickering, Timothy E.; Rastogi, Tanuj; Roulston, Benjamin R.; Ryan, Daniel F.; Rykoff, Eli S.; Sabater, Jose; Sakurikar, Parikshit; Salgado, Jesus; Sanghi, Aniket; Saunders, Nicholas; Savchenko, Volodymyr; Schwardt, Ludwig; Seifert-Eckert, Michael; Shih, Albert Y.; Jain, Anany Shrey; Shukla, Gyanendra; Sick, Jonathan; Simpson, Chris; Singanamalla, Sudheesh; Singer, Leo P.; Singhal, Jaladh; Sinha, Manodeep; Sipocz, Brigitta M.; Spitler, Lee R.; Stansby, David; Streicher, Ole; Sumak, Jani; Swinbank, John D.; Taranu, Dan S.; Tewary, Nikita; Tremblay, Grant R.; De Val-Borro, Miguel; Vasovic, Zlatan; Van Kooten, Samuel J.; Verma, Shresth; Cardoso, Jose Vinicius de Miranda; Williams, Peter K. G.; Wilson, Tom J.; Winkel, Benjamin; Wood-Vasey, W. M.; Xue, Rui; Yoachim, Peter; Zhang, Chen; Zonca, Andrea; Astropy Project Contributors; TARDIS Collaboration; Astropy Coordination Comm.--This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.

The Swift mission, scheduled for launch in 2004, is a multiwavelength observatory for gamma-ray burst (GRB) astronomy. It is a first-of-its-kind autonomous rapid-slewing satellite for transient astronomy and pioneers the way for future rapid-reaction and multiwavelength missions. It will be far more powerful than any previous GRB mission, observing more than 100 bursts yr-1 and performing detailed X-ray and UV/optical afterglow observations spanning timescales from 1 minute to several days after the burst. The objectives are to (1) determine the origin of GRBs, (2) classify GRBs and search for new types, (3) study the interaction of the ultrarelativistic outflows of GRBs with their surrounding medium, and (4) use GRBs to study the early universe out to z > 10. The mission is being developed by a NASA-led international collaboration. It will carry three instruments: a new-generation wide-field gamma-ray (15-150 keV) detector that will detect bursts, calculate 1'-4' positions, and trigger autonomous spacecraft slews; a narrow-field X-ray telescope that will give 5'' positions and perform spectroscopy in the 0.2-10 keV band; and a narrow-field UV/optical telescope that will operate in the 170-600 nm band and provide 0farcs3 positions and optical finding charts. Redshift determinations will be made for most bursts. In addition to the primary GRB science, the mission will perform a hard X-ray survey to a sensitivity of ~1 mcrab (~2 × 10-11 ergs cm-2 s-1 in the 15-150 keV band), more than an order of magnitude better than HEAO 1 A-4. A flexible data and operations system will allow rapid follow-up observations of all types of high-energy transients, with rapid data downlink and uplink available through the NASA TDRSS system. Swift transient data will be rapidly distributed to the astronomical community, and all interested observers are encouraged to participate in follow-up measurements. A Guest Investigator program for the mission will provide funding for community involvement. Innovations from the Swift program applicable to the future include (1) a large-area gamma-ray detector using the new CdZnTe detectors, (2) an autonomous rapid-slewing spacecraft, (3) a multiwavelength payload combining optical, X-ray, and gamma-ray instruments, (4) an observing program coordinated with other ground-based and space-based observatories, and (5) immediate multiwavelength data flow to the community. The mission is currently funded for 2 yr of operations, and the spacecraft will have a lifetime to orbital decay of ~8 yr.

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

Pulsar Tests Gravity Because of their extremely high densities, massive neutron stars can be used to test gravity. Based on spectroscopy of its white dwarf companion, Antoniadis et al. (p. 448 ) identified a millisecond pulsar as a neutron star twice as heavy as the Sun. The observed binary's orbital decay is consistent with that predicted by general relativity, ruling out previously untested strong-field phenomena predicted by alternative theories. The binary system has a peculiar combination of properties and poses a challenge to our understanding of stellar evolution.

We present the final data release of observations of λ21-cm emission from Galactic neutral hydrogen over the entire sky, merging the Leiden/Dwingeloo Survey (LDS: Hartmann & Burton 1997, Atlas of Galactic Neutral Hydrogen) of the sky north of δ = -30° with the Instituto Argentino de Radioastronomía Survey (IAR: Arnal et al. 2000, A&AS, 142, 35; and Bajaja et al. 2005, A&A, 440, 767) of the sky south of δ = -25°. The angular resolution of the combined material is HPBW ~ . The LSR velocity coverage spans the interval -450 km s-1 to +400 km s-1, at a resolution of 1.3 km s-1. The data were corrected for stray radiation at the Institute for Radioastronomy of the University of Bonn, refining the original correction applied to the LDS. The rms brightness-temperature noise of the merged database is K. Residual errors in the profile wings due to defects in the correction for stray radiation are for most of the data below a level of mK. It would be necessary to construct a telescope with a main beam efficiency of % to achieve the same accuracy. The merged and refined material entering the LAB Survey of Galactic H i is intended to be a general resource useful to a wide range of studies of the physical and structural characteristices of the Galactic interstellar environment. The LAB Survey is the most sensitive Milky Way H i survey to date, with the most extensive coverage both spatially and kinematically.

LOFAR, the LOw-Frequency ARray, is a new-generation radio interferometer constructed in the north of the Netherlands and across europe. Utilizing a novel phased-array design, LOFAR covers the largely unexplored low-frequency range from 10–240 MHz and provides a number of unique observing capabilities. Spreading out from a core located near the village of Exloo in the northeast of the Netherlands, a total of 40 LOFAR stations are nearing completion. A further five stations have been deployed throughout Germany, and one station has been built in each of France, Sweden, and the UK. Digital beam-forming techniques make the LOFAR system agile and allow for rapid repointing of the telescope as well as the potential for multiple simultaneous observations. With its dense core array and long interferometric baselines, LOFAR achieves unparalleled sensitivity and angular resolution in the low-frequency radio regime. The LOFAR facilities are jointly operated by the International LOFAR Telescope (ILT) foundation, as an observatory open to the global astronomical community. LOFAR is one of the first radio observatories to feature automated processing pipelines to deliver fully calibrated science products to its user community. LOFAR’s new capabilities, techniques and modus operandi make it an important pathfinder for the Square Kilometre Array (SKA). We give an overview of the LOFAR instrument, its major hardware and software components, and the core science objectives that have driven its design. In addition, we present a selection of new results from the commissioning phase of this new radio observatory.

We present the final data release of observations of lambda 21-cm emission from Galactic neutral hydrogen over the entire sky, merging the Leiden/Dwingeloo Survey (LDS: Hartmann & Burton, 1997) of the sky north of delta = -30 deg with the Instituto Argentino de Radioastronomia Survey (IAR: Arnal et al., 2000, and Bajaja et al., 2005) of the sky south of delta = -25 deg. The angular resolution of the combined material is HPBW ~ 0.6 deg. The LSR velocity coverage spans the interval -450 km/s to +400 km/s, at a resolution of 1.3 km/s. The data were corrected for stray radiation at the Institute for Radioastronomy of the University of Bonn, refining the original correction applied to the LDS. The rms brightness-temperature noise of the merged database is 0.07 - 0.09 K. Residual errors in the profile wings due to defects in the correction for stray radiation are for most of the data below a level of 20 - 40 mK. It would be necessary to construct a telescope with a main beam efficiency of eta_{MB} > 99% to achieve the same accuracy. The merged and refined material entering the LAB Survey of Galactic HI is intended to be a general resource useful to a wide range of studies of the physical and structural characteristices of the Galactic interstellar environment. The LAB Survey is the most sensitive Milky Way HI survey to date, with the most extensive coverage both spatially and kinematically.

CO line emission represents the most accessible and widely used tracer of the molecular ISM. This renders the translation of observed CO intensity into total H 2 gas mass critical to understanding star formation and the ISM in our Galaxy and beyond. We review the theoretical underpinning, techniques, and results of efforts to estimate this CO-to-H 2 “conversion factor,” X CO , in different environments. In the Milky Way disk, we recommend a conversion factor of [Formula: see text] cm −2 (K km s −1 ) −1 with ±30% uncertainty. Studies of other “normal galaxies” return similar values in Milky Way–like disks, but with greater scatter and systematic uncertainty. Departures from this Galactic conversion factor are both observed and expected. Dust-based determinations, theoretical arguments, and scaling relations all suggest that X CO increases with decreasing metallicity, turning up sharply below metallicity ∼1/3–1/2 solar in a manner consistent with model predictions that identify shielding as a key parameter. Based on spectral line modeling and dust observations, X CO appears to drop in the central, bright regions of some but not all galaxies, often coincident with regions of bright CO emission and high stellar surface density. This lower X CO is also present in the overwhelmingly molecular ISM of starburst galaxies, where several lines of evidence point to a lower CO-to-H 2 conversion factor. At high redshift, direct evidence regarding the conversion factor remains scarce; we review what is known based on dynamical modeling and other arguments.

The Cosmic Evolution Survey (COSMOS) is designed to probe the correlated evolution of galaxies, star formation, active galactic nuclei (AGNs), and dark matter (DM) with large-scale structure (LSS) over the redshift range z>0.5-6. The survey includes multiwavelength imaging and spectroscopy from X-ray-to-radio wavelengths covering a 2 deg2 area, including HST imaging. Given the very high sensitivity and resolution of these data sets, COSMOS also provides unprecedented samples of objects at high redshift with greatly reduced cosmic variance, compared to earlier surveys. Here we provide a brief overview of the survey strategy, the characteristics of the major COSMOS data sets, and a summary the science goals. Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which is operated by AURA, Inc., under NASA contract NAS 5-26555 also based on data collected at the Subaru Telescope, which is operated by the National Astronomical Observatory of Japan; the XMM-Newton, an ESA science mission with instruments and contributions directly funded by ESA Member States and NASA; the European Southern Observatory under Large Program 175.A-0839, Chile; Kitt Peak National Observatory, Cerro Tololo Inter-American Observatory, and the National Optical Astronomy Observatory, which are operated by the AURA, Inc., under cooperative agreement with the National Science Foundation; the National Radio Astronomy Observatory, which is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc.; and the Canada-France-Hawaii Telescope with MegaPrime/MegaCam operated as a joint project by the CFHT Corporation, CEA/DAPNIA, the NRC and CADC of Canada, the CNRS of France, TERAPIX, and the University of Hawaii.

Pulsar surveys offer a rare opportunity to monitor the radio sky for impulsive burst-like events with millisecond durations. We analyzed archival survey data and found a 30-jansky dispersed burst, less than 5 milliseconds in duration, located 3 degrees from the Small Magellanic Cloud. The burst properties argue against a physical association with our Galaxy or the Small Magellanic Cloud. Current models for the free electron content in the universe imply that the burst is less than 1 gigaparsec distant. No further bursts were seen in 90 hours of additional observations, which implies that it was a singular event such as a supernova or coalescence of relativistic objects. Hundreds of similar events could occur every day and, if detected, could serve as cosmological probes.

Abstract We present the first Event Horizon Telescope (EHT) observations of Sagittarius A* (Sgr A*), the Galactic center source associated with a supermassive black hole. These observations were conducted in 2017 using a global interferometric array of eight telescopes operating at a wavelength of λ = 1.3 mm. The EHT data resolve a compact emission region with intrahour variability. A variety of imaging and modeling analyses all support an image that is dominated by a bright, thick ring with a diameter of 51.8 ± 2.3 μ as (68% credible interval). The ring has modest azimuthal brightness asymmetry and a comparatively dim interior. Using a large suite of numerical simulations, we demonstrate that the EHT images of Sgr A* are consistent with the expected appearance of a Kerr black hole with mass ∼4 × 10 6 M ⊙ , which is inferred to exist at this location based on previous infrared observations of individual stellar orbits, as well as maser proper-motion studies. Our model comparisons disfavor scenarios where the black hole is viewed at high inclination ( i > 50°), as well as nonspinning black holes and those with retrograde accretion disks. Our results provide direct evidence for the presence of a supermassive black hole at the center of the Milky Way, and for the first time we connect the predictions from dynamical measurements of stellar orbits on scales of 10 3 –10 5 gravitational radii to event-horizon-scale images and variability. Furthermore, a comparison with the EHT results for the supermassive black hole M87* shows consistency with the predictions of general relativity spanning over three orders of magnitude in central mass.

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

Context. Measurement of the Galactic neutral atomic hydrogen (H i) column density, N H i , and brightness temperatures, T B , is of high scientific value for a broad range of astrophysical disciplines. In the past two decades, one of the most-used legacy H i datasets has been the Leiden/Argentine/Bonn Survey (LAB).

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

Abstract We report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15 yr pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. The correlations follow the Hellings–Downs pattern expected for a stochastic gravitational-wave background. The presence of such a gravitational-wave background with a power-law spectrum is favored over a model with only independent pulsar noises with a Bayes factor in excess of 10 14 , and this same model is favored over an uncorrelated common power-law spectrum model with Bayes factors of 200–1000, depending on spectral modeling choices. We have built a statistical background distribution for the latter Bayes factors using a method that removes interpulsar correlations from our data set, finding p = 10 −3 (≈3 σ ) for the observed Bayes factors in the null no-correlation scenario. A frequentist test statistic built directly as a weighted sum of interpulsar correlations yields p = 5 × 10 −5 to 1.9 × 10 −4 (≈3.5 σ –4 σ ). Assuming a fiducial f −2/3 characteristic strain spectrum, as appropriate for an ensemble of binary supermassive black hole inspirals, the strain amplitude is <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:msubsup> <mml:mrow> <mml:mn>2.4</mml:mn> </mml:mrow> <mml:mrow> <mml:mo>−</mml:mo> <mml:mn>0.6</mml:mn> </mml:mrow> <mml:mrow> <mml:mo>+</mml:mo> <mml:mn>0.7</mml:mn> </mml:mrow> </mml:msubsup> <mml:mo>×</mml:mo> <mml:msup> <mml:mrow> <mml:mn>10</mml:mn> </mml:mrow> <mml:mrow> <mml:mo>−</mml:mo> <mml:mn>15</mml:mn> </mml:mrow> </mml:msup> </mml:math> (median + 90% credible interval) at a reference frequency of 1 yr −1 . The inferred gravitational-wave background amplitude and spectrum are consistent with astrophysical expectations for a signal from a population of supermassive black hole binaries, although more exotic cosmological and astrophysical sources cannot be excluded. The observation of Hellings–Downs correlations points to the gravitational-wave origin of this signal.

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

Abstract We present the fourth Fermi Large Area Telescope catalog (4FGL) of γ -ray sources. Based on the first eight years of science data from the Fermi Gamma-ray Space Telescope mission in the energy range from 50 MeV to 1 TeV, it is the deepest yet in this energy range. Relative to the 3FGL catalog, the 4FGL catalog has twice as much exposure as well as a number of analysis improvements, including an updated model for the Galactic diffuse γ -ray emission, and two sets of light curves (one-year and two-month intervals). The 4FGL catalog includes 5064 sources above 4 σ significance, for which we provide localization and spectral properties. Seventy-five sources are modeled explicitly as spatially extended, and overall, 358 sources are considered as identified based on angular extent, periodicity, or correlated variability observed at other wavelengths. For 1336 sources, we have not found plausible counterparts at other wavelengths. More than 3130 of the identified or associated sources are active galaxies of the blazar class, and 239 are pulsars.