We present an incremental version 4FGLDR3 for Data Release 3 of the fourth Fermi Large Area Telescope LAT catalog of ray sources. Based on the first 12 years of science data in the energy range from 50 MeV to 1 TeV it contains 6658 sources. The analysis improves on that used for the 4FGL catalog over eight years of data more sources are fit with curved spectra we introduce a more robust spectral parameterization for pulsars and we extend the spectral points to 1 TeV. The spectral parameters spectral energy distributions and associations are updated for all sources. Light curves are rebuilt for all sources with 1 yr intervals not 2 month intervals. Among the 5064 original 4FGL sources 16 were deleted 112 are formally below the detection threshold over 12 yr but are kept in the list while 74 are newly associated 10 have an improved association and seven associations were withdrawn. Pulsars are split explicitly between young and millisecond pulsars. Pulsars and binaries newly detected in LAT sources as well as more than 100 newly classified blazars are reported. We add three extended sources and 1607 new point sources mostly just above the detection threshold among which eight are considered identified and 699 have a plausible counterpart at other wavelengths. We discuss the degreescale residuals to the global sky model and clusters of soft unassociated point sources close to the Galactic plane which are possibly related to limitations of the interstellar emission model and missing extended sources.
['Gamma-ray sources', 'Gamma-ray astronomy', 'Blazars', 'Pulsars', 'Sky surveys', 'Celestial objects catalogs', '633', '628', '164', '1306', '1464', '212', 'Astrophysics - High Energy Astrophysical Phenomena']
Incremental Fermi Large Area Telescope Fourth Source Catalog
{'No Header': 'Draft version July 4, 2022  \nTypeset using L A T E X twocolumn', 'Incremental Fermi Large Area Telescope Fourth Source Catalog': "S. Abdollahi, 1 F. Acero, 2 L. Baldini, 3 J. Ballet, 2 D. Bastieri, 4, 5, 6 R. Bellazzini, 7 B. Berenji, 8 A. Berretta, 9 E. Bissaldi, 10, 11 R. D. Blandford, 12 E. Bloom, 12 R. Bonino, 13, 14 A. Brill, 15 R. J. Britto, 16 P. Bruel, 17 T. H. Burnett, 18 S. Buson, 19 R. A. Cameron, 12 R. Caputo, 20 P. A. Caraveo, 21 D. Castro, 22, 20 S. Chaty, 2 C. C. Cheung, 23 G. Chiaro, 21 N. Cibrario, 13 S. Ciprini, 24, 25 J. Coronado-Bl'azquez, 26, 27 M. Crnogorcevic, 28 S. Cutini, 29 F. D'Ammando, 30 S. De Gaetano, 10, 11 S. W. Digel, 12 N. Di Lalla, 12 F. Dirirsa, 31 L. Di Venere, 10, 11 A. Dom'ınguez, 32 V. Fallah Ramazani, 33 S. J. Fegan, 17 E. C. Ferrara, 20, 28, 34 A. Fiori, 35 H. Fleischhack, 36, 20, 34 A. Franckowiak, 33 Y. Fukazawa, 37 S. Funk, 38 P. Fusco, 10, 11 G. Galanti, 21 V. Gammaldi, 26, 27 F. Gargano, 11 S. Garrappa, 39 D. Gasparrini, 24, 25 F. Giacchino, 24, 25 N. Giglietto, 10, 11 F. Giordano, 10, 11 M. Giroletti, 30 T. Glanzman, 12 D. Green, 40 I. A. Grenier, 2 M.-H. Grondin, 41 L. Guillemot, 42, 43 S. Guiriec, 44, 20 M. Gustafsson, 45 A. K. Harding, 46 E. Hays, 20 J.W. Hewitt, 47 D. Horan, 17 X. Hou, 48, 49, 50 G. J'ohannesson, 51, 52 C. Karwin, 53 T. Kayanoki, 37 M. Kerr, 23 M. Kuss, 7 D. Landriu, 2 S. Larsson, 54, 55, 56 L. Latronico, 13 M. Lemoine-Goumard, 41 J. Li, 57 I. Liodakis, 58 F. Longo, 59, 60 F. Loparco, 10, 11 B. Lott, 41 P. Lubrano, 29 S. Maldera, 13 D. Malyshev, 38 A. Manfreda, 3 G. Mart'ı-Devesa, 61 M. N. Mazziotta, 11 I.Mereu, 9, 29 M. Meyer, 38 P. F. Michelson, 12 N. Mirabal, 20, 62 W. Mitthumsiri, 63 T. Mizuno, 64 A. A. Moiseev, 34, 28 M. E. Monzani, 12 A. Morselli, 24 I. V. Moskalenko, 12 M. Negro, 34, 62 E. Nuss, 65 N. Omodei, 12 M. Orienti, 30 E. Orlando, 66, 12 D. Paneque, 40 Z. Pei, 5 J. S. Perkins, 20 M. Persic, 59, 67 M. Pesce-Rollins, 7 V. Petrosian, 12 R. Pillera, 10, 11 H. Poon, 37 T. A. Porter, 12 G. Principe, 60, 59, 30 S. Rain'o, 10, 11 R. Rando, 5, 4, 6 B. Rani, 68, 20, 69 M. Razzano, 3 S. Razzaque, 70 A. Reimer, 61, 12 O. Reimer, 61 T. Reposeur, 41 M. S'anchez-Conde, 26, 27 P. M. Saz Parkinson, 71, 72, 73 L. Scotton, 65 D. Serini, 11 C. Sgr'o, 7 E. J. Siskind, 74 D. A. Smith, 41, 75 G. Spandre, 7 P. Spinelli, 10, 11 K. Sueoka, 37 D. J. Suson, 76 H. Tajima, 77, 12 D. Tak, 78, 20 J. B. Thayer, 12 D. J. Thompson, 20 D. F. Torres, 79, 80 E. Troja, 20, 28 J. Valverde, 62, 20 K. Wood, 81 and 82", 'G. Zaharijas': "1 IRAP, Universit'e de Toulouse, CNRS, UPS, CNES, F-31028 Toulouse, France  \n2 Universit'e Paris Saclay and Universit'e Paris Cit'e, CEA, CNRS, AIM, F-91191 Gif-sur-Yvette, France 3 Universit'a di Pisa and Istituto Nazionale di Fisica Nucleare, Sezione di Pisa I-56127 Pisa, Italy 4 Istituto Nazionale di Fisica Nucleare, Sezione di Padova, I-35131 Padova, Italy 5 Dipartimento di Fisica e Astronomia 'G. Galilei', Universit'a di Padova, I-35131 Padova, Italy 6 Center for Space Studies and Activities 'G. Colombo', University of Padova, Via Venezia 15, I-35131 Padova, Italy 7 Istituto Nazionale di Fisica Nucleare, Sezione di Pisa, I-56127 Pisa, Italy 8 California State University, Los Angeles, Department of Physics and Astronomy, Los Angeles, CA 90032, USA 9 Dipartimento di Fisica, Universit'a degli Studi di Perugia, I-06123 Perugia, Italy 10 Dipartimento di Fisica 'M. Merlin' dell'Universit'a e del Politecnico di Bari, via Amendola 173, I-70126 Bari, Italy 11 Istituto Nazionale di Fisica Nucleare, Sezione di Bari, I-70126 Bari, Italy 12 W. W. Hansen Experimental Physics Laboratory, Kavli Institute for Particle Astrophysics and Cosmology, Department of Physics and SLAC National Accelerator Laboratory, Stanford University, Stanford, CA 94305, USA 13 Istituto Nazionale di Fisica Nucleare, Sezione di Torino, I-10125 Torino, Italy 14 Dipartimento di Fisica, Universit'a degli Studi di Torino, I-10125 Torino, Italy 15 NASA Postdoctoral Program Fellow, NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 16 Department of Physics, University of the Free State, P.O. Box 339, Bloemfontein 9300, South Africa 17 Laboratoire Leprince-Ringuet, ' Ecole polytechnique, CNRS/IN2P3, F-91128 Palaiseau, France 18 Department of Physics, University of Washington, Seattle, WA 98195-1560, USA 19 Institut fur Theoretische Physik and Astrophysik, Universitat Wurzburg, D-97074 Wurzburg, Germany 20 NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 21 INAF-Istituto di Astrofisica Spaziale e Fisica Cosmica Milano, via E. Bassini 15, I-20133 Milano, Italy 22 Harvard-Smithsonian Center for Astrophysics, Cambridge, MA 02138, USA 23 Space Science Division, Naval Research Laboratory, Washington, DC 20375-5352, USA  \n24 Istituto Nazionale di Fisica Nucleare, Sezione di Roma 'Tor Vergata', I-00133 Roma, Italy  \njean.ballet@cea.fr  \nPhilippe.Bruel@llr.in2p3.fr  \ntburnett@u.washington.edu  \nlott@cenbg.in2p3.fr", 'Fermi -LAT collaboration': "- 25 Space Science Data Center - Agenzia Spaziale Italiana, Via del Politecnico, snc, I-00133, Roma, Italy\n- 26 Instituto de F´ısica Te´orica UAM/CSIC, Universidad Aut´onoma de Madrid, E-28049 Madrid, Spain\n- 27 Departamento de F´ısica Te´orica, Universidad Aut´onoma de Madrid, E-28048 Madrid, Spain  \n28 Department of Astronomy, University of Maryland, College Park, MD 20742, USA  \n29 Istituto Nazionale di Fisica Nucleare, Sezione di Perugia, I-06123 Perugia, Italy  \n30 INAF Istituto di Radioastronomia, I-40129 Bologna, Italy  \n- 31 Astronomy and Astrophysics Research Development Department, Entoto Observatory and Research Center, Ethiopian Space Science and Technology Institute, Ethiopia  \n32 Grupo de Altas Energ´ıas, Universidad Complutense de Madrid, E-28040 Madrid, Spain  \n- 33 Ruhr University Bochum, Faculty of Physics and Astronomy, Astronomical Institute (AIRUB), D-44780 Bochum, Germany\n- 34 Center for Research and Exploration in Space Science and Technology (CRESST) and NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA  \n35 Dipartimento di Fisica 'Enrico Fermi', Universit`a di Pisa, Pisa I-56127, Italy  \n36 Catholic University of America, Washington, DC 20064, USA  \n37 Department of Physical Sciences, Hiroshima University, Higashi-Hiroshima, Hiroshima 739-8526, Japan  \n- 38 Friedrich-Alexander Universit¨at Erlangen-N¨urnberg, Erlangen Centre for Astroparticle Physics, Erwin-Rommel-Str. 1, 91058 Erlangen, Germany\n- 39 Deutsches Elektronen Synchrotron DESY, D-15738 Zeuthen, Germany  \n40 Max-Planck-Institut f¨ur Physik, D-80805 M¨unchen, Germany  \n41 Universit´e Bordeaux, CNRS, LP2I Bordeaux, UMR 5797, F-33170 Gradignan, France  \n- 42 Laboratoire de Physique et Chimie de l'Environnement et de l'Espace - Universit´e d'Orl´eans / CNRS, F-45071 Orl´eans Cedex 02, France  \n43 Station de radioastronomie de Nan¸cay, Observatoire de Paris, CNRS/INSU, F-18330 Nan¸cay, France  \n44  \nThe George Washington University, Department of Physics, 725 21st St, NW, Washington, DC 20052, USA  \n- 45 Georg-August University G¨ottingen, Institute for theoretical Physics - Faculty of Physics, Friedrich-Hund-Platz 1, D-37077 G¨ottingen, Germany  \n46 Los Alamos National Laboratory, Los Alamos, NM 87545, USA  \n47 University of North Florida, Department of Physics, 1 UNF Drive, Jacksonville, FL 32224 , USA  \n- 48 Yunnan Observatories, Chinese Academy of Sciences, 396 Yangfangwang, Guandu District, Kunming 650216, P. R. China\n- 49 Key Laboratory for the Structure and Evolution of Celestial Objects, Chinese Academy of Sciences, 396 Yangfangwang, Guandu District, Kunming 650216, P. R. China\n- 50 Center for Astronomical Mega-Science, Chinese Academy of Sciences, 20A Datun Road, Chaoyang District, Beijing 100012, P. R. China  \n51 Science Institute, University of Iceland, IS-107 Reykjavik, Iceland  \n- 52 Nordita, Royal Institute of Technology and Stockholm University, Roslagstullsbacken 23, SE-106 91 Stockholm, Sweden  \n53  \nDepartment of Physics and Astronomy, Clemson University, Kinard Lab of Physics, Clemson, SC 29634-0978, USA  \n54  \nDepartment of Physics, KTH Royal Institute of Technology, AlbaNova, SE-106 91 Stockholm, Sweden  \n55  \nThe Oskar Klein Centre for Cosmoparticle Physics, AlbaNova, SE-106 91 Stockholm, Sweden  \n56  \nSchool of Education, Health and Social Studies, Natural Science, Dalarna University, SE-791 88 Falun, Sweden  \n- 57 Department of Astronomy, School of Physical Sciences, University of Science and Technology of China, Hefei, Anhui 230026, China 58 Finnish Centre for Astronomy with ESO (FINCA), University of Turku, FI-21500 Piikii¨o, Finland  \n59 Istituto Nazionale di Fisica Nucleare, Sezione di Trieste, I-34127 Trieste, Italy  \n60 Dipartimento di Fisica, Universit`a di Trieste, I-34127 Trieste, Italy  \n- 61 Institut f¨ur Astro- und Teilchenphysik, Leopold-Franzens-Universit¨at Innsbruck, A-6020 Innsbruck, Austria\n- 62 Department of Physics and Center for Space Sciences and Technology, University of Maryland Baltimore County, Baltimore, MD 21250, USA  \n63 Department of Physics, Faculty of Science, Mahidol University, Bangkok 10400, Thailand  \n- 64 Hiroshima Astrophysical Science Center, Hiroshima University, Higashi-Hiroshima, Hiroshima 739-8526, Japan  \n65  \nLaboratoire Univers et Particules de Montpellier, Universit´e Montpellier, CNRS/IN2P3, F-34095 Montpellier, France  \n66  \nIstituto Nazionale di Fisica Nucleare, Sezione di Trieste, and Universit`a di Trieste, I-34127 Trieste, Italy  \n67  \nOsservatorio Astronomico di Trieste, Istituto Nazionale di Astrofisica, I-34143 Trieste, Italy  \n68  \nKorea Astronomy and Space Science Institute, 776 Daedeokdae-ro, Yuseong-gu, Daejeon 30455, Republic of Korea  \n69  \nDepartment of Physics, American University, Washington, DC 20016, USA  \n- 70 Centre for Astro-Particle Physics (CAPP) and Department of Physics, University of Johannesburg, PO Box 524, Auckland Park 2006, South Africa\n- 71 Santa Cruz Institute for Particle Physics, Department of Physics and Department of Astronomy and Astrophysics, University of California at Santa Cruz, Santa Cruz, CA 95064, USA  \n72 Department of Physics, The University of Hong Kong, Pokfulam Road, Hong Kong, China  \n73 Laboratory for Space Research, The University of Hong Kong, Hong Kong, China  \n74  \nNYCB Real-Time Computing Inc., Lattingtown, NY 11560-1025, USA  \n- 75 Laboratoire d'Astrophysique de Bordeaux, Universit'e de Bordeaux, CNRS, B18N, all'ee Geoffroy Saint-Hilaire, F-33615 Pessac, France  \n76 Purdue University Northwest, Hammond, IN 46323, USA  \n77 Solar-Terrestrial Environment Laboratory, Nagoya University, Nagoya 464-8601, Japan  \n78 Department of Physics, University of Maryland, College Park, MD 20742, USA  \n79 Institute of Space Sciences (ICE, CSIC), Campus UAB, Carrer de Magrans s/n, E-08193 Barcelona, Spain; and Institut d'Estudis Espacials de Catalunya (IEEC), E-08034 Barcelona, Spain  \n80 Instituci'o Catalana de Recerca i Estudis Avan¸cats (ICREA), E-08010 Barcelona, Spain  \n81 Praxis Inc., Alexandria, VA 22303, resident at Naval Research Laboratory, Washington, DC 20375, USA  \n82 Center for Astrophysics and Cosmology, University of Nova Gorica, Nova Gorica, Slovenia", 'ABSTRACT': 'We present an incremental version (4FGL-DR3, for Data Release 3) of the fourth Fermi -LAT catalog of γ -ray sources. Based on the first twelve years of science data in the energy range from 50 MeV to 1 TeV, it contains 6658 sources. The analysis improves on that used for the 4FGL catalog over eight years of data: more sources are fit with curved spectra, we introduce a more robust spectral parameterization for pulsars, and we extend the spectral points to 1 TeV. The spectral parameters, spectral energy distributions, and associations are updated for all sources. Light curves are rebuilt for all sources with 1 yr intervals (not 2 month intervals).  \nAmong the 5064 original 4FGL sources, 16 were deleted, 112 are formally below the detection threshold over 12 yr (but are kept in the list), while 74 are newly associated, 10 have an improved association, and seven associations were withdrawn. Pulsars are split explicitly between young and millisecond pulsars. Pulsars and binaries newly detected in LAT sources, as well as more than 100 newly classified blazars, are reported. We add three extended sources and 1607 new point sources, mostly just above the detection threshold, among which eight are considered identified, and 699 have a plausible counterpart at other wavelengths.  \nWe discuss the degree-scale residuals to the global sky model and clusters of soft unassociated point sources close to the Galactic plane, which are possibly related to limitations of the interstellar emission model and missing extended sources.  \nKeywords: Gamma rays: general - surveys - catalogs', '1. INTRODUCTION': 'The Fermi Large Area Telescope (LAT) has been surveying the high-energy γ -ray sky since 2008 (Atwood et al. 2009), and the LAT Collaboration has published a succession of source catalogs based on comprehensive analyses of LAT data. The fourth source catalog (4FGL, where FGL stands for Fermi Gamma-ray LAT) was derived from the analysis of the first 8 yr of LAT science data and contained 5064 sources (Abdollahi et al. 2020, hereafter the 4FGL paper, or simply 4FGL). Every FGL catalog until and including 4FGL had used a new analysis method, new calibrations, a new diffuse model, and even a new reconstruction of the events themselves (like Pass 8 for 4FGL). Therefore each successive version of the catalog was largely independent of the previous one. Results for different versions were compared only after the fact.  \nSince the development of the 4FGL, the data have been stable, no new model for the Galactic interstel-  \nlar emission was built, and the analysis method has not evolved much. So we decided to change our approach and start building incremental 4FGL versions every two years, until one major analysis or data component changes. The Data Release 2 (DR2; Ballet et al. 2020) covered 10 yr of data and contained 723 additional sources. We refer to the original 4FGL sources as Data Release 1 (DR1).  \nIn this paper, we describe the incremental procedure and apply it to 12 yr of data, creating the Data Release 3 (4FGL-DR3 or DR3 for short). The 4FGL paper remains the reference for the detailed methodology. DR2 was strictly incremental (all DR1 sources are included and the analysis procedure was exactly the same). Meanwhile we have improved a few analysis steps, and recent publications have revealed a few new extended sources. Therefore we have relaxed this condition in DR3, to allow for deleting a few sources, and accept a few improved localizations.  \nAll the information present in the DR1 catalog is updated, except the 2 month light curves. Those are very costly (in terms of CPU and disk space), and we have shown in the 4FGL paper that the 1 yr light curves capture most of the variability information. Monthly light curves, including indications of spectral variability, are now provided with similar accuracy in the Fermi -LAT Light Curve Repository 1 . Detection of transient sources on monthly time scales is covered by 1FLT 2 (Baldini et al. 2021). Transients on weekly time scales are covered by FAVA 3 (Abdollahi et al. 2017). Phase-folded spectral analysis of periodic sources (pulsars and binaries) is provided in dedicated papers.  \nThe condition for switching the spectral models to a curved representation is now less stringent, thereby increasing the fraction of curved spectra to about half. We improved the spectral parameterization of pulsars, quantified the residuals using a new method, and now include in the Catalog FITS file the peak energy (in νF ν ) for each source.  \nThe community has worked hard on the unassociated LAT sources in recent years. Optical observations allowed for classifying many blazars, and radio, optical, and γ -ray timing allowed for detecting new pulsars and binary systems. We have included all this information in the source associations. The DR2 and DR3 results provide about 900 more unassociated sources to look into. Unassociated sources close to the Galactic plane are spatially clustered and have distinct spectral characteristics. We have singled them out by a specific flag. Section 2 describes the data and the updates to the diffuse model, Section 3 describes the updates to the analysis including the new extended sources, Section 4 describes the results including a comparison to DR1, Section 5 describes the updates to the associations, and Section 6 discusses the unassociated sources, with particular emphasis on the category of soft clustered sources close to the Galactic plane.', '2.1. The LAT Data': 'The data for the 4FGL-DR3 catalog were taken during the period 2008 August 4 (15:43 UTC) to 2020 August 2 (8:33 UTC) covering twelve years. During most of this time, Fermi was operated in sky-scanning survey mode, with the viewing direction rocking north and south of  \nthe zenith on alternate orbits such that the entire sky is observed every ∼ 3 hr. Starting on 2018 March 16, the Fermi spacecraft was put in safe hold after one solar array drive became stuck. Scientific operations of the LAT were interrupted for more than three weeks, by far the longest missing time interval since 2008, and restarted on April 8 in partial sky-scanning mode 4 . The Sun more rarely enters the field of view, and during some phases of the ∼ 53-day precession period of the orbit, the entire sky is not covered every three hours. This hampers the solar observations and short-term light curves, but has little impact on the integrated sky coverage.  \nAs in 4FGL, intervals around solar flares and bright γ -ray bursts (GRBs) were excised. During the additional four years, 108 ks were cut due to two successive bright solar flares in 2017 September, and 2.4 ks around 5 new bright GRBs. The current version of the LAT data is Pass 8 P8R3 (Atwood et al. 2013; Bruel et al. 2018). For DR2 and DR3, we used the P8R3 V3 IRFs instead of P8R3 V2 used to build DR1. The relative effective areas for the point-spread function (PSF) event types were recalibrated on bright sources 5 , reducing the systematics from 20% to 5%. The energy range remains 50 MeV to 1 TeV. In 4FGL, the event types were not treated separately above 10 GeV. We now decompose the events into 4 components up to 30 GeV, resulting in 19 components in total (Table 1), for better precision but at the cost of more disk space and CPU.  \n2.2. Model for the Diffuse Gamma-Ray Background  \nTable 1. 4FGL-DR3 Summed Likelihood components  \nNote -Weused 19 components (all in binned mode) in the Summed Likelihood approach ( § 3.3) for DR3. Components in a given energy interval share the same number of energy bins, the same zenith angle selection, and the same region of interest (RoI) core, but have different pixel sizes in order to adapt to the PSF width. Each filled entry under Pixel Size corresponds to one component of the summed log-likelihood. NBins is the number of energy bins in the interval, ZMax is the zenith angle cut, Ring Width refers to the difference between the RoI core and the extraction region (see Abdollahi et al. 2020, for details).  \nFigure 1. Model for the inverse Compton emission of the Sun at 1 GeV averaged over one year, in equatorial coordinates. The top and bottom rows are, respectively, years 12 (2019-2020) and 9 (2016-2017), chosen to be after and before the solar array drive anomaly. The color bar is in sqrt scale, and the intensity unit is ph cm -2 s -1 sr -1 MeV -1 .  \n<!-- image -->  \nThe interstellar emission is a very important ingredient when characterizing the point sources, because it is structured and 90% of the photons in the Galactic plane are of diffuse origin. We used exactly the same model as in 4FGL 6 . The isotropic spectrum was updated to P8R3 V3. The model for the emission of the Sun and the Moon was also kept the same, neglecting the modulation of the emission along the solar cycle. It became neces-  \nsary, however, to generate a different effective model for each year (for the light curves). After the solar array drive anomaly during the tenth year, the coverage of the Sun changed significantly (Figure 1). For intervals of about two weeks, the Sun does not enter the field of view at all.', '3. CONSTRUCTION OF THE CATALOG': 'Most of the steps were identical to 4FGL (Abdollahi et al. 2020). As before the extended sources are taken from the literature. We only search for point sources. We use the same Test Statistic (TS), TS = 2 ln( L / L 0 ), to quantify how significantly a source emerges from the background, comparing the maximum value of the likelihood function L including the source in the model with L 0 , the value without the source. The 1 yr light curves were generated in exactly the same way.', '3.1. Extended Sources': "Several extended sources (existing and new) were the subject of publications since the 4FGL paper. They were not considered in DR2, but were considered for DR3:  \n- · We adopted the Gaussian model proposed by Devin et al. (2020) for the broad supernova remnant (SNR) G150.3+4.5, after we confirmed that it is a better fit than the disk model used in DR1.\n- · Tibaldo et al. (2018) have studied the Vela X pulsar wind nebula (PWN) and reported that the radio template was superior to the geometric templates. We found that it does not provide a better  \nfit than the DR1 model, but it eliminates six dubious point sources, so it is simpler and we adopted it. The same authors reported a spatial evolution as a function of energy, but we found that their high-energy disk is not significant when using the radio template.  \n- · Araya (2020) has shown that the cluster of DR1 sources around the SNR G279.0+1.1 is actually an extended LAT source. We adopted their ring model.\n- · HESS J1640 -465 is seen as extended by H.E.S.S., but was considered a point source in the DR1 catalog. We adopted the H.E.S.S. Gaussian template, after it was shown by Mares et al. (2021) that it is a better fit to the LAT data as well.\n- · Two papers discussed the Fermi -LAT emission from the PWN HESS J1825 -137. The morphology considered by Araya et al. (2019) was too  \ncomplex for our catalog. We tried the broader Gaussian proposed by Principe et al. (2020) but found it worsened the fit. We did change the template, however, because we noticed that the one that we used since 2FGL was not the final version of Grondin et al. (2011).  \n- · Ambrogi et al. (2019) have shown that the previous template for the SNR HB 21 was too broad, because it mistakenly included a separate point source. We adopted their smaller template.\n- · Xin et al. (2019) reported the detection with Fermi -LAT of the extended VERITAS source VER J2227+608, corresponding to the SNR G106.3+2.7 next to the much brighter pulsar PSR J2229+6114. We have included it in the DR3 model.  \nTable 2 . Extended Sources Entered or Modified in the 4FGL-DR3 Analysis  \nNote - List of updated and new sources that have been modeled as spatially extended. The Origin column gives the name of the Fermi -LAT catalog in which that extended source was first introduced (with a different template). The Extent column indicates the radius for Disk (flat disk) sources, the 68% containment radius for Gaussian sources, the inner and outer radii for Ring (flat annulus) sources, and an approximate radius for Map (external template) sources.  \nTable 2 lists the source name, origin, spatial template, and the reference for the dedicated analysis. These sources are tabulated with the point sources, with the only distinctions being that no position uncertainties are reported and their names end in e (see Appendix A). The new and modified extended sources have DataRelease set to 3. Unidentified point sources inside extended ones are indicated as 'xxx field' in the ASSOC2 column of the catalog, where 'xxx' is the name of the extended source. A total of 78 extended sources are considered in DR3.", '3.2. Detection and Localization': "The source detection procedure followed the same lines as in 4FGL. It used pointlike (Kerr 2010) and a specialized diffuse model in which the nontemplate features are estimated differently. It started from the previous list of sources (DR1 for DR2, DR2 for DR3), relocalized them over 12 yr of data, looked for peaks in the residual TS maps generated for several spectral shapes, introduced those in the model, refit, and iterated over the full procedure. The resulting list (called uw1216) contained 11,128 seed sources at TS > 10.  \nThat procedure naturally resulted in changing the positions of all sources. Since we wanted the catalog to be incremental, we forced most DR1 and DR2 sources back to their original positions (consistent with their names). This requires associating the uw1216 seed sources with the DR2 catalog. To do that, we started by defining  \nFigure 2. Comparison of the pointlike localizations of uw1216 seeds (for DR3) with those of uw8606 (for DR1). The abscissa is the ratio of TS values, whereas the ordinate is the inverse of the ratio of the areas of the error ellipses (such that larger is better). Sources in the upper right corner (materialized by the dashed lines) are stronger and much better localized in uw1216, and their positions and error ellipses were updated in DR3. The source at TS ratio around 300 is the bright γ -ray binary LS 5039, which appeared double in uw8606 because of its complex spectrum (resulting in low TS) but was well localized.  \n<!-- image -->  \npairs in which each member is the nearest counterpart of the other, within at most 1 · . We repeated that procedure after removing the pairs found at the first step, in order to handle faint sources close to much brighter ones (in which the nearest counterpart can be the bright source). We rejected 58 pairs in which the distance between the two members was larger than the largest of the two 99.9% position errors (approximately 1.5 times the 95% error) and 0 . · 15. We eliminated the uw1216 seeds in a pair with a DR2 source, and formed the input list from the DR2 sources (including 139 with no counterpart in uw1216) and the remaining uw1216 seeds.  \nIn DR2 we actually used the original positions of all DR1 sources strictly, and kept all of them. However a number of sources became much brighter in recent years than they were during the first 8 yr, so they are now much better localized. In recognition of this, we adopted for DR3 the new error ellipse when its area was at least ten times smaller and TS at least ten times larger than the same quantities in DR1. We used those very strict conditions to preserve the incremental nature of the catalog. This selected ten sources, clear outliers when looking at the distribution of ratios of ellipse areas and TS (Figure 2). All are associated with blazars (very variable sources). One of them was in uw8606 (list of seeds for DR1) but too faint to be included in DR1 (it appeared in DR2), so only nine DR1 sources were updated. They  \nhave a 4FGL counterpart in the ASSOC FGL column, but have DataRelease set to 3 because they have moved.  \nWe have also adopted the new error ellipses of one DR1 source close to the modified extended source HESS J1825 -137 and three faint soft sources (two from DR1, one from DR2) whose 95% semi-major axes R 95 were larger than 1 · in DR1 or DR2. Three DR1 sources with broad error ellipses R 95 ∼ 0 . · 4 were split into two better localized seeds inside the original error ellipse. We replaced the DR1 sources by the two seeds. The same applied to two DR2 sources. Conversely we replaced two DR1 sources (4FGL J1831.5 -0935 and J1830.2 -1005) by a single uw1216 seed compatible with PSR J1831 -0952.  \nWe applied the same procedure as in 4FGL to eliminate seeds too close to a bright source and inside extended sources. Since we changed the extended sources ( § 3.1), this resulted in deleting several DR1 sources as well, not only new seeds. More precisely, 4FGL J1640.6 -4632 was replaced by HESS J1640 -465, five DR1 and two DR2 sources were deleted in or close to Vela X, three DR1 and one DR2 sources in the new SNR G279.0+1.1. We also deleted the faint 4FGL J1814.1 -1710 next to the extended source HESS J1813 -178 (even though this one did not change) because its TS decreased from 30 in DR1 to 16 in DR3, and it was not spectrally distinct from HESS J1813 -178.  \nFollowing the catXcheck verification ( § 3.6) we have applied several modifications to fit the data better. We replaced 4FGL J0857.7 -4507 (in one of the clusters described in § 6.2) by a stronger uw1216 source about 0 . · 5 away. We replaced 4FGL J1750.0 -3849 just outside the extended source FHES J1741.6 -3917 by a stronger point source just inside it.  \nWe have looked at individual faint (TS < 25) sources with a very soft (Γ > 3 . 1) spectrum, a large error ellipse ( R 95 > 0 . · 5) or a stronger close neighbor in uw1216 and deleted three DR1 (4FGL J0533.9+2838, J0313.6 -7508, and J2326.5+8555) and three DR2 (4FGL J0750.0+7140, J1112.0+1021, and J1752.2 -3002) sources. Nearby seeds appeared in DR3 and fit the data better. Similarly, we have deleted a faint very hard (Γ < 1 . 3) DR2 source (4FGL J1830.3 -1601) that was superseded by a brighter nearby DR3 seed.  \nAs before we manually added 11 known LAT pulsars that could not be localized by the automatic procedure without phase selection. Only one of those (PSR J1731 -4744) reached TS > 25 in DR3. It appears with a NULL error ellipse. The 10 others were discarded. A total of 5029 new seeds were entered to the gtlike source characterization for DR3 in addition to the 5690 remaining DR1 and DR2 point sources.  \nTable 3. 4FGL-DR3 priors on diffuse parameters  \nNote -The Galactic interstellar emission model is modulated by a power law K ( E/E 0 ) -Γ d , in which E 0 is set to 1 GeV. K tends to be smaller than 1. This is probably because of the sources that we added since the model was devised (just before 4FGL).  \nWe reassessed the systematic corrections to localization. Up to 4FGL, they were defined globally over the entire sky, so their values were dominated by highlatitude sources (92% of associated point sources in DR1 are at | b | > 5 · ). However, confusion and strong interstellar emission make localization much more difficult in the Galactic plane. In DR2 we estimated localization systematics separately at | b | < 5 · , based on 110 pulsars and low-latitude active galactic nuclei (AGN). We found that they should indeed be increased there, to 27 '' for the absolute 95% error, and 1.37 for the systematic factor. We did not change the systematics in the highlatitude sky (25 '' and 1.06). The DR2 systematics are used again for DR3.", '3.3. Thresholding and Light Curves': 'We used the Fermi Tools 1.4.7 analysis suite 7 . Compared to the Science Tools v11r7p0 that were used in 4FGL, the main difference is the new edisp bins setting for energy dispersion, which allows for broadening of the energy interval considered in the model space. The Science Tools setting was equivalent to edisp bins = -1, not enough to fully fill the first and last data bins for the small energy bins that we use (10 per decade). We used edisp bins = -2 instead, considering two bins for energy dispersion on each side beyond the data bins. The acceptance for the PSF3 event type increases by a factor 2.4 between 50 and 100 MeV, so to account for this rapid change, we used 6 energy bins in the first component (Table 1) instead of 3 and set edisp bins to -4 accordingly.  \nThe likelihood weights were recomputed over 12 yr of data, resulting in slightly smaller weights everywhere. We reoptimized all RoIs, resulting in 1988 RoIs containing up to 10 sources in their core. The maximum of 8 sources in the RoI core enforced in 4FGL was not optimal (it led to too many small neighboring RoIs and too many iterations).  \nFor DR1 and DR2, the normalization of one of the diffuse components sometimes became very large (close to the maximum, set to 2) or very small (close to 0) in regions where that component contributed little. This was particularly obvious for the isotropic component in the Galactic ridge. This could happen because of missing features in the interstellar emission model (IEM), leading to a higher isotropic component in the fit. This was not satisfying since the isotropic emission is supposed to be approximately constant over the sky, and the stuck parameters caused trouble in the error estimates. In order to stabilize the diffuse parameters, we replaced the hard limits with Bayesian priors. We chose the priors from the distributions of the parameters (in DR2) in regions where they are dominant ( | b | < 10 · for the Galactic parameters and | b | > 30 · for the isotropic normalization). This led to the values reported in Table 3. The intrinsic reasons for the large or small values did not go away, of course, so the resulting isotropic normalizations were still larger than one in the Galactic ridge, but they rarely exceeded 1.5. The worst RoIs around ( l, b ) ∼ (315 · , 3 · ) reached 1.7.  \nA major difference with the 4FGL procedure is that we retained all previous 4FGL sources (DR1 and DR2, except those discussed in § 3.2), so that they would not be deleted from the model even if they have TS < 25. The resulting catalog contains 6659 entries, among which 1610 are new (1607 point sources plus the three extended sources discussed in § 3.1), and 181 have TS < 25 (112 DR1 and 69 DR2 sources).  \nThe light curves over 1 yr bins were recomputed (even the first 8 yr because adding new sources changes the fluxes of DR1 sources). The procedure was very similar to 4FGL, except we considered events up to 1 TeV, and we used a different model for the Sun and Moon at each year (see § 2.2). Neither change has a major effect. Recomputing light curves over 2 month bins would have been very time consuming, and we showed in the 4FGL paper that the finer time binning was not critical to detecting variability. Monthly (and shorter intervals) light curves are provided in the Fermi -LAT Light Curve Repository.  \nThe threshold on Variability Index is now 24.725 (same 1% probability of false detection, applied to 12 bins). The number of significantly variable sources in-  \nFigure 3. Ratio of the error on the Γ S PLEC parameter (Eq. 2) in DR3 to that on the asymptotic low-energy index Γ = Γ S -d/b in DR2 (diamonds). The dots show the same error ratio on the α LP parameter (Eq. 1), to gauge the effect of increased exposure (small). The large scatter at significance larger than 100 is due to the fact that those sources have free exponential index b in DR3.  \n<!-- image -->  \ncreased from 1443 to 1695 (among which 158 are DR2 or DR3 sources). The fraction of variable sources decreased because part of the variability information in 4FGL came from the 2 month light curves. The distribution of fractional variability ( Frac Variability in the FITS file) is similar to 4FGL, peaking between 50% and 90%. The blazars all have fractional variability larger than 10% except 4C +55.17 (4FGL J0957.6+5523). The flux of that source, first discussed by McConville et al. (2011), decreased by about 13% over 10 yr, not showing any flare.', '3.4. Spectral Fitting and Energy Distributions': 'The spectral analysis proceeded along the same lines as in 4FGL, with five important differences detailed below: we changed the parameterization for pulsars, we fit more sources with four parameters, we changed the TS curv threshold for using a curved spectral shape, we added new columns reporting the peak energy in νF ν , and we added one spectral bin to the spectral energy distributions (SEDs).  \nMost significantly curved sources are still fit with a lognormal function ( LogParabola under SpectrumType in the FITS table, hereafter LP):  \nd N d E = K ( E E 0 ) -α -β ln( E/E 0 ) . (1)  \nThe significantly curved pulsars are still fit with a subexponentially cutoff power law ( PLSuperExpCutoff under SpectrumType in the FITS table, hereafter PLEC). However we noticed that the parameters of the  \nparameterization used in 4FGL ( PLSuperExpCutoff2 in the Fermi Tools) were strongly correlated, particularly when the exponential index ( b in Eq. 2) was free. This hampered convergence and precluded freeing b except in the very brightest sources. To work around this, we devised a new parameterization ( PLSuperExpCutoff4 in the Fermi Tools):  \nd N d E = K ( E E 0 ) d b -Γ S exp [ d b 2 ( 1 -( E E 0 ) b )] (2)  \nd N d E = K ( E E 0 ) -Γ S -d 2 ln E E 0 -db 6 ln 2 E E 0 -db 2 24 ln 3 E E 0 for ∣ ∣ ∣ ∣ b ln E E 0 ∣ ∣ ∣ ∣ < 10 -2 , (3)  \nin which the normalization K is directly the flux density at the reference energy E 0 and the shape parameters are the spectral slope Γ S = dlog( dN/dE ) / dlog E and the spectral curvature d = d 2 log( dN/dE ) / d(log E ) 2 at E 0 .  \nThe development in Eq. 3 shows explicitly that for b = 0 the expression converges to LP with d = 2 β . For b < 0 the shape (in log) is reversed (it decreases exponentially toward low energies and as a PL toward high energies). The parameters K , Γ S , d and b appear as PLEC Flux Density , PLEC IndexS , PLEC ExpfactorS and PLEC Exp Index in the FITS table, respectively. The reference energy E 0 is chosen so that the error on K is minimal and appears as Pivot Energy in the FITS table. This parameterization proved indeed much more stable, and the correlation between parameters is considerably reduced.  \nFigure 3 shows that on average Γ S is much better defined (error ratio below 1) than the low-energy asymptotic index Γ = Γ S -d/b used in 4FGL. On average the error on α (∆ α ) decreased by only 2.3% between 10 and 12 yr, whereas ∆Γ S is on average 0.626 ∆Γ. Note that 15 pulsars had no ∆Γ at all in DR2 because Γ was stuck at 0. This problem does not exist in the new parameterization. The small error ratio at significance below 40 corresponds to PSR J0622+3749. It had Γ = 0.05 in DR2, very close to the boundary. At low significance, a fraction of the pulsars have ∆Γ in DR2 that is similar or less than ∆Γ S in DR3. The DR2 errors may have been underestimated because of the large correlation between parameters.  \nSince the new PLEC parameters describe the same mathematical quantities (spectral slope and curvature) as the LP parameters (except the curvature is 2 β ), it is interesting to compare them. Both spectral models are fit to all sources, so this information is readily available. Figure 4 shows that, as expected, the correlation between α and Γ S is quite tight. The correlation between 2 β and d is not as tight, because the curvature has larger error. On average d < 2 β , because the curvature in the  \n<!-- image -->  \nFigure 4. Left: comparison of α of the LP model and Γ S of the PLEC model ( LP Index and PLEC IndexS in the FITS file), showing that these parameters are largely similar. Right: comparison of LP β and PLEC d ( LP beta and PLEC ExpfactorS in the FITS file), showing that d is well correlated to 2 β . The dashed lines show a one-to-one correlation. The outliers have large errors on both plots (none is farther off than 2 σ ).  \n<!-- image -->  \nPLEC model depends on energy, and the pivot energy (where d is defined) is below the energy at which the two curvatures would be the same. This effect is stronger for bright sources that have lower pivot energy for the same spectral shape.  \nThe new PLEC parameterization allows a 4-parameter fit (with free b ) in fainter sources. We have applied it to all pulsars with TS > 10,000 (28, up from 6 in DR1 and DR2). At TS < 10,000 ∆ b becomes larger than 0.15, which is the natural scatter on b in the brightest pulsars, so freeing b is no longer beneficial. For all other significantly curved pulsars, b is fixed to 2/3 as in 4FGL. For comparison, the median b over the 28 pulsars with free b is 0.51, its weighted (with 1/ σ 2 ) average is 0.55, and its intrinsic dispersion is 0.16. The Small Magellanic Cloud has b = 1, as before (Abdo et al. 2010).  \nBesides 3C 454.3, the five other AGN with TS > 80,000 were modeled with PLEC and free b as well: CTA 102, Mkn 421, S5 0716+71, 3C 279 and PKS 1424 -41. The new parameterization contains the LP model (for b = 0) so there is no risk of a worse fit, but nonconvergence can occur when the model is not constrained well enough. The TS threshold for free b is higher in AGN than pulsars, because the curvature is much less in AGN, so ∆ b is larger. The main objective of fitting PLEC with free b remains to improve the modeling of the surroundings of very bright sources at low energy.  \nIn the DR1 and DR2 catalogs, the sources were represented with a curved spectral model (LP or PLEC) rather than a power law (PL) when TS curv = 2 (log L (curved spectrum) -log L (PL)) was larger than 9 (3 σ ). For DR3 we have lowered that threshold to 4 (2 σ ) for two reasons, illustrated by Figure 5:  \n- 1. The left-hand plot shows that all sources are curved to some degree, when bright enough that the curvature is significant. This is physically understandable (the PL is an idealized asymptotic case).\n- 2. The PL approximation tends to overestimate the source flux outside of the energy range in which the data are well constrained. The right-hand plot quantifies that on the sum of all sources, at three different TS curv thresholds. Between 200 MeV and 20 GeV, it makes very little difference. Highlatitude blazars (blue) are not strongly curved compared to pulsars and other Galactic sources (green and black). But above 20 GeV and below 200 MeV, the TS curv threshold makes a difference. Most conspicuously, the global spectrum in the plane outside the ridge (green) diverges toward low energies at a threshold of 9 (what was actually used in DR2). The global spectrum is better behaved when the TS curv threshold is lower.  \nThe implication is that we should choose a TS curv threshold as low as possible while preserving global convergence of the model. For DR3 we set the threshold to 4 as a compromise because the test with a threshold at 1 required manual intervention, and global convergence could not be guaranteed. This is still a significant improvement as demonstrated in Figure 5 (right).  \nThanks to the improved statistics and the lowered TS curv threshold, many more sources are considered significantly curved in DR3 than in DR1. The number of LP spectral models is now 3130 (among which 650 are new sources), while the number of PLEC models is now  \n<!-- image -->  \nFigure 5. Left: fraction of spectrally curved sources as a function of significance, outside the Galactic plane. The error bars are based on the source counts in each bin, assuming a binomial distribution. The sharp increase followed by a plateau at low significance (instead of the expected regular increase) is due to the fact that the source significance is boosted (moved to the right) for curved sources, combined with the detection threshold and a decreasing dN/d log TS . Right: impact of changing the TS curv threshold from 1 to 9 on the global source spectrum (sum of all sources, excluding the diffuse model) in three regions of the sky; high ( | b | > 10 · ) in blue, ridge ( | b | ≤ 2 · and | l | < 60 · ) in black, and plane ( | b | ≤ 10 · excluding the ridge) in green. This test was carried out on the 10 yr data set and DR2 sources. The plot starts at 10 MeV, lower than the data (50 MeV).  \n<!-- image -->  \nFigure 6. Ratio between the energy flux (0.1 to 100 GeV) in the PL approximation and that in the LP or PLEC representation, for those sources that are considered significantly curved, outside the Galactic plane. Plus signs refer to the LP model, diamonds to the PLEC model. The scatter increases toward low significance because of the larger statistical fluctuations.  \n<!-- image -->  \n258 (among which 2 are new sources). The fraction of sources represented by a curved spectral model increased from 30% to 51% overall.  \nChanging the spectral models for many sources from PL to curved also has an impact on their energy flux estimate. The PL flux is nearly always larger, due to the integration down to 100 MeV for soft sources (or up to 100 GeV for hard sources Γ < 2). This is quantified in Figure 6. Well-measured pulsars (diamonds at high significance) have on average a 50% larger flux estimate  \nin the PL model than in the PLEC model. The wellmeasured blazars (plus signs) are in the lower branch at a ratio around 1.1, because their spectra are much less curved than those of pulsars. The diamonds in the lower branch are the very bright blazars fit with PLEC (see above). A ratio less than one occurs when the peak of the curved spectrum is close to either 100 MeV or 100 GeV (the boundaries of the energy interval over which the energy flux is extracted), because the full analysis interval (50 MeV to 1 TeV) extends beyond those boundaries and maintains the PL fit below the peak of the curved model.  \nTo provide more information on curved spectra, we now report systematically in the catalog the peak energy in νF ν and its uncertainty for all sources (including those not significantly curved) and both models as (Unc )LP EPeak and (Unc )PLEC EPeak :  \nE peak (LP)= E 0 exp ( 2 -α 2 β ) (4)  \nE peak (PLEC)= E 0 ( 1 + b d (2 -Γ S ) ) 1 /b (5)  \nThe uncertainties are obtained using the covariance matrix. We set those columns to NULL when the peak energy is undefined (when Γ ≥ 2 in PLEC or curvature is upwards) or found too far from the LAT energy range ( | ln( E peak /E 0 ) | > 10). E peak (PLEC) is undefined in only 4 PLEC spectra, but in 516 LP and 1682 PL spectra (mostly because Γ ≥ 2). E peak (LP) is always defined in PLEC and LP spectra, but remains undefined in 850  \nFigure 8. Distribution of spectral fit quality over all DR3 sources (histogram). The χ 2 distributions for 4 (dashed line, good fit) and 6 (solid line, bad fit) degrees of freedom are overlaid. The vertical lines correspond to a p -value of 10 -3 . The threshold for Flag 10 in Table 4 corresponds to the dashed vertical line (for 4 dof).  \n<!-- image -->  \nFigure 7. SED of the BL Lac object PG 1218+304 (4FGL J1221.3+3010) for DR2 (purple crosses and PL fit) and DR3 (green diamonds and LP fit).  \n<!-- image -->  \nPL spectra in which β is found negative or close to 0. The uncertainty on E peak (LP) is also set to NULL when β is stuck to 1 and has no error (166 occurrences in LP spectra, 6 in PLEC spectra, and 3 in PL spectra). Even when E peak is formally defined, it should be treated with caution when Γ(PLEC) gets close to 2 or the curvature is close to 0.  \nIn the DR1 and DR2 catalogs, we provided SEDs over N bands = 7 energy bands. Thanks to the increased statistics after 12 yr (particularly valuable at high energy where the precision is limited by counts rather than background), we can provide them over 8 energy bands in DR3, extending to 1 TeV as in the 3FHL catalog (Ajello et al. 2017). We replaced the 30 - 300 GeV bin by a 30 - 100 GeV bin and a 100 GeV - 1 TeV bin. We computed those two fluxes using unbinned likelihood, as in 4FGL. In the highest energy band, 172 sources have TS > 25. An example of such a hard source is shown in Figure 7. That source has TS curv (LP) = 8.70, below the DR2 threshold of 9 but above the DR3 threshold of 4, so it is fit with the LP model in DR3. The significant differences between 1 and 10 GeV are due to the natural variability of that strong source.  \nThe systematic uncertainties 8 have improved when moving from P8R3 V2 to P8R3 V3 and are now 0.1 in bands 1 and 8 and 0.05 in bands 2-7. We compute the spectral fit quality (SFQ) as before (Eq. 5 in Abdollahi et al. 2020). Figure 8 shows how that quantity is distributed. Up to DR2, we flagged (Flag 10 in Table 4) sources with SFQ above the 99.9% quantile (10 -3 p -value) of the χ 2 distribution with N bands -2 degrees  \nof freedom (dof), considering that most sources were fit with PL spectra that have two free parameters. For DR3 that would result in 6 dof. Clearly the data do not follow the χ 2 (6) distribution, but are much closer to the χ 2 (4) distribution. This is due to two effects: first, half the sources are now fit with curved spectra (3 free parameters); second, the first ( < 100 MeV) and last ( > 100 GeV) bands do not behave as a Gaussian distribution (very large errors in the first band, and very low counts in the last band), so they contribute little to SFQ. For DR3 we flagged sources with SFQ > 18.47, the 10 -3 p -value of the χ 2 (4) distribution. The number of sources that trigger this remains small (only 16). The largest SFQ is obtained in PSR J0205+6449 at 49.6. That pulsar indeed has a known PWN (3C 58), which is very significant above the PLEC fit in the two bands above 30 GeV. The other sources flagged with Flag 10 have one spectral band off by more than 3 σ .', '3.5. Analysis Flags': "As before, we attempt to flag the sources affected by systematic errors larger than the statistical errors that we report, and which should therefore be treated with more caution. This can impact all source parameters, and even call into question the existence of the faintest flagged sources with significance less than 6 σ or 7 σ (TS < 50). Because the largest source of systematics is the underlying diffuse emission, sources close to the Galactic plane, where diffuse emission is both the strongest and most structured, are most affected.  \nTable 4. Comparison of the numbers of flagged sources between DR1 and DR3, separately for the DR1 sources ('DR1 in DR3' column) and the sources introduced in DR2 and DR3 ('Post DR1' column).  \n- a In the FITS file, the values are encoded as individual bits in the Flags column, with Flag n having value 2 ( n -1) .  \nThe flags ( Flags in the FITS file) are recalled in Table 4 (see the 4FGL paper for the detailed definitions), together with the numbers of sources flagged for each reason and their evolution since DR1. The effect of the underlying IEM was estimated by launching the procedure described in § 3.3 for a second time using the same seeds but the previous IEM (gll iem v06).  \nWe note six changes with respect to the 4FGL procedure:  \n- · For DR1, Flag 1 was applied to the comparison with both pointlike and the alternative IEM. For DR3, the comparison to pointlike was performed with no TS selection in order to apply Flag 3 to sources at TS < 25. Flag 1 was determined only by the comparison with the old IEM, after selecting TS > 25 on both sides.\n- · Flag 2 mostly relies on the comparison between the new uw1216 pointlike positions of the 4FGL sources ( § 3.2) and the 4FGL positions (uw8606), because no localization was performed over 12 yr with the old IEM. The comparison with the old IEM also flags a few sources (although the seed positions are exactly the same) when among two nearby seeds one survives in DR3 and the other one with the old IEM.\n- · All sources have been visually screened for potential deficiencies in the underlying diffuse model. The visual screening for diffuse features (Flag 6) depends on the source parameters. It was entirely redone for DR2. To be conservative, we left that flag (and the c suffix) in all the sources  \nthat had it in DR1, even if they would not have been flagged in DR2, and added the flag to all the new cases, including the sources already in DR1. For DR3 we checked only the new sources and did not change Flag 6 in the DR1 and DR2 sources. Two sources in the Large Magellanic Cloud (4FGL J0517.9 -6930c and J0535.7 -6604c) have that flag set as well.  \n- · Flag 12 still applies to PLEC sources with Γ = Γ S -d/b ≤ 0 (as well as to LP sources with β fixed to 1, as before), even though the new PLEC parameterization (Eq. 2) has no explicit limit at Γ = 0.\n- · We introduced for DR2 a new Flag (13) to flag explicitly sources from a previous version with TS < 25 in the current version. It now applies to DR2 sources as well (69 have TS < 25 in DR3).\n- · For DR3, we introduce a new Flag (14) to flag explicitly soft Galactic unassociated (SGU) sources in regions with high source density ( § 6.2). The three strongest SGUs (4FGL J1702.7 -5655 at b = -9 . 2, 4FGL J1823.3 -1340 inside HESS J1825 -137, and 4FGL J1839.4 -0553 inside HESS J1841 -055), all at TS > 800, are very clear point sources (good pulsar candidates) and are not flagged.  \nOverall the fraction of sources with any flag set increased from 23% in DR1 to 32% in DR3. Among the new DR2 and DR3 sources, 75% are flagged in the Galactic plane ( | b | < 10 · ), but only 26% are flagged above the Galactic plane.", '3.6. All-sky verification': "In order to verify that the 4FGL-DR3 catalog provides a good representation of the γ -ray sky, we have performed an independent all-sky analysis, named catXcheck , using the p -value statistic (PS) data/model deviation estimator developed by Bruel (2021). catXcheck consists of the analysis of 540 RoIs (12 × 12 · , with a pixel size of 0 . · 1) covering the whole sky. The RoI centers lie on 19 Galactic parallels, whose latitudes range from -90 · to 90 · with a 9 · step. The longitude step is 9 · for | b | ≤ 18 · , 10 · for | b | = 27 · and 36 · , and it is 12 · , 15 · , 18 · , 24 · , and 45 · for | b | = 45 · , 54 · , 63 · , 72 · , and 81 · , respectively. The difference with respect to the catalog main analysis is that we combine all PSF event types together, using all events whose zenith angle is less than 90 · to suppress contamination from Earth limb γ rays. The energy range is [100 MeV, 1 TeV], and we use 10 bins per decade.  \nThe source model of each RoI comprises the Galactic diffuse emission, the isotropic template, the Sun and Moon steady emission templates, as well as all point-like and extended sources from DR3, within 5+0 . 015( σ src -4) degrees of the RoI border, where σ src is the source significance ( Signif Avg in the FITS version of the catalog). That formula selects sources outside of the RoI up to 5 · for the faintest ones and to 15 · for the brightest ones, corresponding respectively to the 68% and 95% containment radii at 100 MeV.  \nSince we want to check the goodness of the DR3 information, the spectral parameters of the DR3 sources are fixed to their catalog values, and the only free parameters of the binned-likelihood fit are the ones of the Galactic diffuse emission (power-law correction) and isotropic (normalization) templates. After the fit, we compute the 3D map of the number of photons predicted by the source model and compute the PS map in order to detect any significant deviation between the observed and predicted 3D maps.  \nAs explained in Bruel (2021), the PS estimator is computed for each spatial pixel of a RoI, and rather than comparing the observed and predicted count spectra at that pixel, we compare spatially integrated count spectra around that pixel, with an integration radius varying with energy as p 0 ( E/ 100 MeV) -p 1 ⊕ p 2 , where ⊕ is the sum in quadrature. We use the point-source optimized integration parameterization (which is very close to the PSF 68% containment energy dependence): p 0 = 4 · , p 1 = 0 . 9 and p 2 = 0 . · 1, and we merge the energy bins to have an energy binning of 0.3 in log 10 E . PS is defined such that | PS | = -log 10 ( p -value), where the p -value is the probability that the statistical fluctuations can reach a level of deviation as large as the one observed in the  \nFigure 9. PS distribution for the three Galactic-latitude selections. The thin solid histogram corresponds to the expected 10 -| x | distribution. All histograms are normalized such that their integral is 1.  \n<!-- image -->  \ndata, under the assumption that the model represents the data. The PS sign is the sign of the sum of the count spectra residuals in units of σ . The likelihood weights are taken into account in the PS computation. We chose to use the PS map to investigate data/model deviations because, contrary to the TS map, it is sensitive to both positive and negative deviations.  \nFor one pixel, the 3 σ , 4 σ , and 5 σ thresholds correspond to PS = 2.57, 4.20 and 6.24. Since the level of correlation between pixels is negligible with p 2 = 0 . · 1, we only have to take into account the number of trials when deriving the all-sky significance thresholds. The number of RoIs and their centers have been chosen so that the overlap between neighboring RoIs ensures that, for any given direction in the sky, there is at least one RoI containing this direction and providing a distance to the RoI border larger than ∼ 1 . · 5. As a consequence, when combining the results of the 540 PS maps, we loop over the RoIs and the RoI pixels, and we only consider the pixels for which the RoI provides the largest distance to the RoI border. This allows us to avoid double counting the sky directions. The total number of considered pixels is 4,175,468 and adding its logarithm to the one-pixel significance thresholds gives the all-sky 3 σ , 4 σ , and 5 σ thresholds: PS = 9.19, 10.82, and 12.86.  \nThe PS distributions for samples of sources in three different ranges of Galactic latitude are shown in Figure 9. The high-latitude ( | b | > 35 · ) distribution follows the 10 -| x | expectation, whereas some deviations are visible in the low-latitude ( | b | < 5 · ) and mid-latitude (5 · < | b | < 35 · ) distributions. To show the locations of these deviations, we build an all-sky HEALPix (G'orski et al. 2005) map in Galactic coordinates with N side = 256 from the 540 PS maps. Figure 10 shows the re-  \nFigure 10. All-sky PS map with 'source' settings in Galactic coordinates (Mollweide projection). The color scale corresponds to -4 σ to 4 σ deviations (accounting for the number of pixels in the map).  \n<!-- image -->  \nFigure 11. All-sky PS map with 'diffuse' settings in Galactic coordinates (Mollweide projection). The color scale corresponds to -4 σ to 4 σ deviations (accounting for the number of independent pixels in the map).  \n<!-- image -->  \nll-sky PS map. Overall the data/model agreement is good over the whole sky. Using a 4 σ threshold, we find two positive and five negative deviations. The positive ones are both located inside bright extended sources: 4FGL J1923.2+1408e (W 51C) and 4FGL J0822.1 -4253e (Puppis A). They occur because the geometric templates that we use for those sources are far from perfect, and we do not allow additional point sources inside them ( § 3.2). Regarding the negative deviations, the one at a distance of ∼ 1 · from the Galactic center corresponds to only one pixel. The other four are clusters of several pixels with PS < -10 . 82 and correspond to extended deviations explaining the negative tail in the mid-latitude histogram in Figure 9. They all are close to large molecular clouds: one can be associated with Cepheus, one with Perseus, and two with Orion B.  \nSince these extended negative deviations seem to be related to the Galactic diffuse emission, we also compute the 540 PS maps with the set of PS spatial integration parameters optimized for large extended deviations ( p 0 = 5 · , p 1 = 0 . 8 and p 2 = 1 · ). With p 2 = 1 · , the level of correlation between pixels is expected to increase. Assuming that pixels within 0 . · 5 of each other are highly correlated, we divide the number of trials by 69, and the 3 σ , 4 σ , 5 σ thresholds become PS = 7.4, 9, 11. The corresponding all-sky PS map, shown in Figure 11, exhibits many deviations above 4 σ . A preliminary analysis of the negative ones shows that many of them are correlated with some components of the IEM, especially the local CO template (Galactocentric radius range 7-9 kpc), which may indicate that allowance for variation of the X CO factor among the local clouds is necessary. The positive deviations are less numerous and less significant. Like the negative ones, some are likely due to imperfections in the modeling of the Galactic diffuse emission, but some could correspond to new extended sources.  \nAlthough a complete analysis of these extended deviations is outside the scope of this paper, we try to quantify their impact on the DR3 catalog. In order to do so, we look at the 540 PS maps and count the number of pixels with | PS | > 7 . 4 and derive the expected number of DR3 sources in these regions from the average source density in the | PS | < 7 . 4 region. In the 40 RoIs centered along the Galactic plane, the fraction of negative and positive deviation pixels is 2% and 0.9%, respectively. The numbers of observed and expected DR3 sources are 15 and 24.2 in the negative deviation pixels, and 33 and 10.1 in the positive deviation pixels. In the RoIs farther from the Galactic plane, the fraction of negative and positive deviation pixels is 0.15% and 0.03%, respectively. There are no DR3 sources in the negative deviation pixels, while 8.4 are expected. We conclude that the catalog results are biased in the negative and positive deviation regions, but the effect is small and the number of missed or spurious sources is at most ∼ 20.", '4. THE 4FGL-DR3 CATALOG': 'The catalog is available online 9 , together with associated products. It contains 6659 entries (1607 new since DR1). The source designation is 4FGL JHHMM.m+DDMM . The format is similar to 4FGL (except it does not contain the 2 month light curves) but not identical. It is detailed in App. A. A new column ( DataRelease ) is set to 1 for the unchanged DR1 sources, 2 for DR2 sources,  \n9  \nSee  \n12yr catalog/.  \nFigure 12. Distributions of the 4FGL-DR3 Test Statistic at high latitude ( | b | > 10 · ) for the sources already in DR1 (dashed) and those new to DR2 or DR3 (solid). The vertical dotted line is the threshold at TS = 25. The solid histogram extends below the threshold because 69 DR2 sources have TS < 25 in DR3, including 46 at | b | > 10 · .  \n<!-- image -->  \nand 3 for DR3 sources. Sources existing in DR1 or DR2 but whose position (and therefore name) changed ( § 3.2) have DataRelease set to 3. The DR1 or DR2 names of those sources are given in the ASSOC 4FGL column of the FITS file. They are considered DR1 or DR2 sources in the text, tables, and plots.', '4.1. The 4FGL-DR1 sources': 'The TS of the 5048 remaining DR1 sources in DR3 increased on average compared to DR1. Among sources that are fit with the PL model in both catalogs (no TS boost from switching to a curved model), TS increased on average by 33% at high latitude ( | b | > 10 · ) and by 22% closer to the Galactic plane. The difference between the two is due to the lower log-likelihood weights, which partly offset the exposure increase close to the Galactic plane where diffuse emission is strong, particularly for soft sources. Even at high latitude, the TS ratio is less than expected from the exposure increase (50%). This is not due to the lower weights (the TS ratio does not depend on index) but to the selection bias (sources tend to be brighter in the interval in which they were defined) and signal-splitting with new sources. The statistical uncertainties on the parameters decreased by 15% at high latitude and 9% close to the plane.  \nThe random character of adding new data and the variable nature of the γ -ray sources inevitably lead to broadening the TS distribution, in addition to shifting it to larger values. The sharp selection threshold at TS = 25 in DR1 becomes a tail of sources at TS < 25 in DR3. Figure 12 shows the TS distribution outside the Galactic plane, where 74 DR1 sources have TS < 25.  \nOver the entire sky, 112 DR1 sources have TS < 25 in DR3, among which 15 have TS < 16, and the smallest TS is 9.8 (4FGL J1816.4 -0057).  \nThe energy flux decreased on average by 5% (0.2 σ ) between DR1 and DR3 (considering only PL sources to avoid the effect of switching to a curved model). This is again a selection bias dominated by faint sources: close to the detection threshold, DR1 sources can get either brighter or fainter in DR3, but faint sources in DR3 cannot be fainter in DR1 because they would have TS < 25. Real variability is apparent for bright sources (TS > 1000). The scatter on the energy flux ratio is 5 σ on variable DR1 sources, but only 1 σ on nonvariable ones. There is no offset and less scatter (only 0.8 σ overall) on the difference of PL indices because source fluxes vary strongly, but spectral shapes only slightly.', '4.2. The new DR2 and DR3 sources': 'The 1607 new DR2 and DR3 sources are on average very close to the detection threshold (median TS of 34). This is obvious in Figure 12. There are a few newly active bright blazars at TS > 100, and even three at TS > 800 (4FGL J0154.5 -6604, J1544.3 -0649, and J2250.0 -1250 associated with PKS 0153 -663, NVSS J154419 -064913, and PKS 2247 -131 respectively) but the majority are faint sources that rise above the detection threshold thanks to the additional four years. Outside the Galactic plane, the confusion limit is still distant, because ever-deeper exposures reach to higher energies at which the PSF is better. Among the 1042 faint (TS < 100) new sources at | b | > 10 · , only 7% have the confusion flag (Flag 5 in Table 4) set.  \nFigure 13 (left) shows that the energy flux distribution of the new sources peaks between 1 and 2 × 10 -12 erg cm -2 s -1 . The majority (140) of the 156 new DR2 and DR3 sources considered variable (mostly blazars) are outside the Galactic plane ( | b | > 10 · ). Their average flux is larger than that of the nonvariable sources. This is largely a selection bias (variability cannot be detected in the very faint sources).  \nFigure 13 (right) compares the distribution of PL indices of the new sources with that of the DR1 ones. The variable new sources are soft on average, similar to flat-spectrum radio quasars (FSRQ). The same effect is observed in the 1FLT catalog (Baldini et al. 2021), even more strongly because detections in 1FLT are based on one-month data sets. The nonvariable new sources have a peak on the soft side as well (probably similar to the variable ones, but not bright enough to be significantly variable). But they also have a shoulder on the hard side. This is because adding more exposure is more beneficial at high energies (in the count-dominated regime)  \n<!-- image -->  \nFigure 13. Left: distributions of the 4FGL-DR3 energy flux at high latitude ( | b | > 10 · ) for the sources already in DR1 (dotted) and those new to DR2 or DR3, separately for those considered variable (dashed) or not (solid). The energy flux is obtained with the best spectral model, which can be either a curved model or a power law ( § 3.4). Right: same distributions for the PL index (always obtained with the power-law model). The DR1 distribution was broadened to account for the larger uncertainties of the index values for the new sources, so the comparison is fair.  \n<!-- image -->  \nthan at low energies, and therefore helps with detecting hard sources (like BL Lac objects).', '4.3. Step-by-step from DR1 to DR3': 'To understand the effects of the analysis changes introduced in DR3 with respect to 4FGL, we have considered their effects on the DR2 data set (ten years). To that end, we started with the same seeds as the DR2 catalog, changed each element in sequence (in the order of the list below), and compared each intermediate result with the previous one.  \n- · We first switched from the Science Tools to the Fermi Tools, and from the P8R3 V2 to the P8R3 V3 data and Instrument Response Functions (changing the isotropic templates accordingly). This led to essentially no change at all, as expected.\n- · Setting edisp bins to -2 ( § 3.3) when accounting for energy dispersion had a significant effect. Prior to that, the model counts were underestimated by about 5% in the first energy band (50 to 100 MeV), so the model fit was too high by the same amount at low energy. We found 16 fewer sources with the new settings, due to a slight TS reduction in soft sources. The energy flux decreased by 2% on average. The largest effect was spectral: we found 155 more LP spectra (+ 10%). Hard PL sources (Γ < 2 . 1) were unchanged, but Γ decreased more and more for softer sources, reaching ∆Γ = -0 . 1 at Γ = 3. Note that this is specific to the catalog setting that uses many components  \nfor the summed likelihood in relatively narrow energy bands (Table 1). Analyses using few broader energy bands are not as sensitive to edisp bins .  \n- · Introducing diffuse priors had very little effect on the sources. We found 7 fewer sources, 11 fewer LP spectra. In the Galactic ridge (where the effect of the priors is strongest), the average energy flux of sources decreased by 2% (0.1 σ ), and the average PL index decreased (became harder) by 0.005 (0.06 σ ). This is completely innocuous.\n- · Lowering the TS curv threshold from 9 to 4 ( § 3.4) had an important effect on the sources, besides the collective effect illustrated in Figure 5 (right). The number of LP spectra increased by 1000 (+ 58%). The total number of sources increased by 127 (+ 2%) because more LP sources allowed more soft neighbors to reach the detection threshold. The energy fluxes of sources that switched from PL to LP decreased by 1.8 σ (well above the statistical precision) as illustrated in Figure 6. Their neighbors picked up most of that flux, so their energy fluxes increased. They also became somewhat softer in the process (by 0.02 on average). On the other hand, the neighbors of those neighbors became a little harder, so that the global average PL index did not change.  \nIn conclusion, the change with the greatest consequences was lowering the TS curv threshold. We are confident that the new setting is closer to the physical reality.', '5. ASSOCIATIONS': 'We have performed the same association procedure as in 4FGL, using the Bayesian method and likelihoodratio method (see Abdollahi et al. 2020). The procedure yielded associations for 341 new DR2 sources and 366 new DR3 sources, representing association fractions of 48% and 41% respectively. These fractions are roughly consistent with that expected from DR1 for these lowTS sources. We provide low-probability (0 . 1 < P < 0 . 8) associations for 95 new DR2 or DR3 sources and associations with Planck sources for 44 others. One major change regarding classes is that now we split the generic pulsar class used in previous catalogs into two separate classes, namely young pulsars (PSR) and millisecond pulsars (MSP). We recall that some pulsars with LATdetected pulsations do not fill the defining significance criterion of this general source catalog and are thus omitted. The full list of LAT-detected pulsars is regularly updated, 10 and the current population will be discussed in the forthcoming third LAT pulsar catalog. We also point out that associations with pulsar candidates are based on spatial coincidence only and that γ -ray pulsations from these sources are not necessarily expected to be detected soon.', '5.1. Census of new DR2/DR3 sources relative to 4FGL-DR1': "For completeness, we recapitulate below the census of associations of new DR2 sources already given in the DR2 document (Ballet et al. 2020), with some updates on blazar classes:  \n- 1. two pulsar candidates, PSR J1439 -5501 and PSR J1904 -11;\n- 2. one globular cluster (GLC), NGC 362;\n- 3. one high-mass binary (HMB), PSR B1259 -63; although detected by the LAT during its 2010 periastron passage (Abdo et al. 2011), this source is included in an FGL catalog for the first time thanks to the gain in significance from its 2017 passage (Johnson et al. 2018);\n- 4. three SNRs, 3C 397, SNR G001.4 -00.1, and SNR G003.7 -00.2; they are considered as reliable associations (as opposed to SPPs, see below) because their radio sizes are quite small (4 ' -14 ' ) and their positions well match those of the 4FGL sources;\n- 5. two star-forming regions, NGC 346, which is the brightest star-forming region in the SMC, and Sh 2 -152;  \n- 6. two galaxies, IC 678 and NGC 5380;\n- 7. one starburst galaxy, Arp 299;\n- 8. two radio galaxies, NGC 3078 and NGC 4261;\n- 9. 281 blazars, including 176 blazars of unknown type (BCUs), 66 BL Lac objects (BLL), and 39 FSRQs; 11\n- 10. 16 SPPs (SPP stands for 'SNR, Pulsar or PWN' and refers to sources of unknown nature but overlapping with known SNRs or PWNe and thus candidate members of these classes);\n- 11. 30 sources of unknown nature (UNK, | b | < 10 · sources solely associated with the likelihood-ratio method from large radio and X-ray surveys).  \nThe pulsar PSR J1757 -2421, formerly associated with 4FGL J1756.6 -2352 in DR1, is now associated with 4FGL J1757.9 -2419. The source 4FGL J1756.6 -2352 is now a SPP.  \nBelow is a census of the associations for the new DR3 sources:  \n- 1. three pulsars with detected pulsations, including two young pulsars, PSR J0922+0638 and PSR J1731 -4744 (Smith et al. 2019), and one MSP, PSR J1455 -3330;\n- 2. one MSP candidate, PSR J1957+2516;\n- 3. five GLCs, NGC 1851, NGC 5286, NGC 6205, NGC 6712, and M 54;\n- 4. two PWNe, PWN G327.1 -1.1 and PWN G54.1+0.3;\n- 5. one SNR, the Kepler SNR; the DR3 source 4FGL J1730.4 -2131 was associated by the Likelihood Ratio method with NVSS J173036 -212910; this NVSS source is actually part of the Kepler SNR (extended for the Very Large Array at 3 ' diameter), so we report the SNR as the counterpart (see also Acero et al. 2022);\n- 6. three novae, V1369 Cen (Cheung et al. 2016), V906 Car (Jean et al. 2018; Aydi et al. 2020), and YZ Ret (MGAB-V207; Li et al. 2020); the position of the DR3 source (4FGL J0358.5 -5432) is closer to that of a blazar candidate, PMN J0358 -5434, than to that of YZ Ret, but the light curve clearly points to the nova association;  \n- 7. one high-mass X-ray binary, 1RXS J172006.1 -311702 (Esposito et al. 2014); the nature of this source is debated as it was classified as a cataclysmic variable in Fortin et al. (2018);\n- 8. one radio galaxy, LEDA 55267, classified as a Fanaroff-Riley type 0 (Baldi et al. 2019) with a LAT detection reported in Paliya (2021), and one AGN, NGC 6454;\n- 9. 306 blazar candidates including 219 BCUs, 51 BLLs, and 36 FSRQs;\n- 10. 22 SPPs;\n- 11. 20 sources of unknown nature.\n- 5.2. Association/classification changes in 4FGL-DR1 sources\n- 1. A specific class (GC) has been created for the Galactic center (4FGL J1745.6 -2859), replacing its former SPP class.\n- 2. The increase of localization systematic uncertainties for sources lying at | b | < 5 · ( § 3.2) has led to the association of five extra BCUs, five SPPs, and the suppression of the association of 4FGL J1804.9 -3001 with the GLC NGC 6528. Four sources classified as unknown in DR1 are now classified as BCUs because they gained Bayesianbased associations. Two other formerly unknown sources are now associated with a pulsar and a binary star (Kleinmann's star).\n- 3. Two associations with pulsar candidates for 4FGL J1618.7 -4633 and 4FGL J1720.8 -1937 have been discarded on the realization that their estimated γ -ray luminosities would exceed their spin-down powers by factors greater than 30.\n- 4. Pulsations have been detected for 22 more DR1 sources, including 10 PSR and 12 MSP (e.g., Clark et al. 2021; Tabassum et al. 2021).\n- 5. A total of 29 associations with MSP from the West Virginia University list 12 have been added (three replacing associations with BCUs and one with a SPP), along with one association with a young pulsar candidate.\n- 6. The pulsar PSR J1909 -3744 was mistakenly associated with 4FGL J1912.2 -3636. The position of this LAT-detected pulsar was used as a seed, but  \n- the resulting 4FGL source ended up too far away from the seed to make the association plausible. The 4FGL source is now unassociated.\n- 7. The association of the MSP binary recently discovered by Wang et al. (2020) with 4FGL J0935.3+0901 has been implemented. Its class is binary (BIN).\n- 8. Following Bogdanov & Halpern (2015), 4FGL J1544.5 -1126 is a candidate transitional MSP binary (1RXS J154439.4 -112820). Its class, BCU in DR1, has been changed to 'low-mass X-ray binary' (LMB). Another candidate transitional MSPbinary, CXOU J110926.4 -650224, was found (Coti Zelati et al. 2019, 2021) in a systematic search for X-ray counterparts to unassociated sources in the preliminary LAT 8 yr point-source list (FL8Y). While the position of the X-ray source lies well within the 95% error ellipse of the FL8Y source, it is outside of the ellipse of the corresponding DR3 source 4FGL J1110.3 -6501. In consequence, 13 we have not retained this association in our tables.\n- 9. Another transitional MSP binary, 1SXPS J042749.2 -670434, displaying simultaneous optical, X-ray, and γ -ray eclipses has been identified by Kennedy et al. (2020) in 4FGL J0427.8 -6704 (previously unassociated).\n- 10. Gamma 2 Velorum (4FGL J0809.5 -4714; Mart'ıDevesa et al. 2020) has been classified as a binary system.\n- 11. We have implemented the results of Corbet et al. (2019), who discovered that 4FGL J1405.1 -6119 is a HMB.\n- 12. The source 4FGL J1832.9 -0913, associated with the TeV source HESS J1832 -093, has been classified as a HMB (see Mart'ı-Devesa & Reimer 2020).\n- 13. Swihart et al. (2021) found a candidate redback MSP binary in a follow-up observation of 4FGL J0940.3 -7610. Other redbacks have been discovered in a similar fashion in 4FGL J0407.7 -5702 (Miller et al. 2020) and 4FGL J0540.0 -7552 (Strader et al. 2021). These associations have been classified as LMBs.  \n- 14. Similarly, Li et al. (2021) found a black widow binary at the position of 4FGL J0336.0+7502, which has then also been classified as a LMB.\n- 15. Three periodic optical and X-ray sources were recently discovered positionally consistent with DR1 sources. These redback MSP candidates have been classified as BIN. They are 4FGL J0212.1+5321 (Li et al. 2016), 4FGL J0523.3 -2527 (Strader et al. 2014), and 4FGL J0838.7 -2827 (Halpern et al. 2017).\n- 16. The tentative association of 4FGL J0647.7 -4418 with the HMB RX J0648.0 -4418 called out in DR1 has been replaced by association with the BCU SUMSS J064744 -441946 following the multiwavelength investigation of Mart'ı et al. (2020).\n- 17. Following Jarvela et al. (2020), we have reclassified TXS 2116 -077 (4FGL J2118.8 -0723) as a Seyfert galaxy instead of a NLSY1.\n- 18. The latest version of the Radio Fundamental Catalog 14 has enabled the association of six previously unassociated sources and two SPPs with blazar candidates.\n- 19. Recent follow-up observations of 4FGL blazars (Pe˜na-Herazo et al. 2017, 2019, 2020, 2021a,b,c; Desai et al. 2019; Rajagopal et al. 2021) have enabled the classification of 240 former BCUs, 2 AGNs, and 2 UNKs into 214 BLLs and 30 FSRQs.\n- 20. Two FSRQs (PKS 0736 -770, TXS 1530 -131) have been reclassified as BCUs and one other (Two Micron All Sky Survey, 2MASS J02212698+2514338) as a BLL.  \nNote: A typo has been found in § 6.2 of the 4FGL paper concerning the name of the star cluster associated with the extended H ii region encompassing 4FGL J1115.1 -6118. The star cluster is NGC 3603 (and not NGC 4603).  \nTable 5 lists the class tallies of the whole DR3 catalog.", '5.3. Spectral properties': "In this section, the spectral properties for the main source classes in the whole 4FGL-DR3 sample are presented. We use these properties to shed light on the nature of the unassociated sources in the next section.  \nFigure 14. Photon-index distributions for different source classes.  \n<!-- image -->  \nAlthough some classes, most notably pulsars, show significant levels of spectral curvature, the PL photon index, Γ, is a convenient parameter enabling comparison across classes. Figure 14 compares the photon-index distributions for different established classes. Pulsars show remarkably consistent soft spectra, along with GLCs. MSP have somewhat harder spectra on average compared to PSR and GLCs. The observed difference between the MSP and GLC distributions seems to challenge the conventional wisdom that MSP dominate the GLC's γ -ray emission (e.g., Harding et al. 2005; Venter et al. 2009), and is possibly due to different ambient stellar radiation fields comptonized by the MSP electrons (Bednarek & Sitarek 2007). The spectra of SNRs and PWNe are notably harder than those of pulsars, in agreement with the observation that several SNRs and almost all PWNe are detected in the TeV domain. The distribution for blazars spans a broad range of photon index, the distributions for FSRQs and BLLs exhibiting moderate overlap (see Ajello et al. 2020). Blazars constitute the only abundant class showing a hard-spectrum component.  \nTurning to the spectral curvature, while some blazars show significant curvature, as do most pulsars, the two classes occupy different branches in the LP TS curv versus TS plane (Figure 15 top). These branches merge at the low-TS end as curvature cannot be assessed with high significance in that case. Despite the fact that PSR and MSP show slightly different spectra as can be seen from their photon-index distributions, both fall on the same branch in Figure 15 top.", '6. UNASSOCIATED SOURCES': "The variation of the photon-index distribution of unassociated sources with Galactic latitude indicates that different populations coexist (see Figure 21 of the  \nTable 5. LAT 4FGL-DR3 Source Classes  \nNote -The designation 'spp' indicates potential association with SNR or PWN. 'Unknown' are | b | < 10 · sources solely associated with the likelihood-ratio method from large radio and X-ray surveys. Designations shown in capital letters are firm identifications; lower-case letters indicate associations.  \n4FGL paper), with the spectra becoming softer on average for lower Galactic latitudes. This observation justifies treating separately the cases of high- and lowlatitude unassociated sources. We set a somewhat arbitrary limit at | b | = 10 · , with 1122 and 1018 unassociated point sources lying at higher and lower latitudes respectively.", '6.1. High-latitude unassociated sources': 'Drawing on the class-wise discriminating power of the photon-index distributions illustrated in Figure 14, the distribution for high-latitude unassociated sources is compared to that of BCUs, which are the most relevant blazar set given the low TS, in Figure 16 top. The good agreement between the two distributions supports the conclusion that most high-latitude unassociated sources are indeed blazars 15 . This conclusion is slightly qualified  \nwhen considering the loci of the unassociated sources in the TS versus LP TS curv plane (Figure 15 bottom). A significant fraction ( glyph[similarequal] 70%) of the higher-TS (TS > 100) unassociated sources fall onto the pulsar branch, strongly suggesting that they are MSP. This idea is also supported by considering the evolution of the photonindex and variability-index distributions with TS (Figure 17). While for the lower-TS bins, the distributions quite closely match those of BCUs, clear differences are observed for the higher-TS bin (100 < TS < 500), where the distributions are compatible with those expected for MSP. In conclusion, their photon indices, significance of spectral curvature, and lack of variability make many of the high-latitude, 100 < TS < 500 unassociated sources good MSP candidates.', '6.2. Low-latitude unassociated sources': "Since some AGN catalogs suffer from significant incompleteness close to the plane due to extinction, some background blazars detected through the plane could  \nFigure 15. Top: parameter LP TS curv (reflecting the spectral-curvature significance) as a function of Test Statistic for blazars and pulsars (no selection on the Galactic latitude was applied). The dotted line, selected by eye, approximately delineates the separation between the two corresponding branches. Bottom: same for the high-latitude unassociated sources. The same dotted line as in the top panel is displayed for orientation.  \n<!-- image -->  \nelude association and contribute to the unassociated sample. Based on the tallies observed at high latitude and the difference in flux limits due to the brighter diffuse emission background close to the plane, the number of | b | < 10 · blazars is estimated to be 340 ± 20, while 399 blazars are actually reported as associations in DR3. Comparing the low-latitude photon-index distribution with that obtained for higher latitudes (normalized so that the total yield of Γ < 2.4 sources match), an excess of soft sources (Γ glyph[similarequal] 2.5) is found in the low-latitude sample relative to the other. From the difference of the spectral-index distributions, this excess, which was less pronounced in 3FGL and is very likely due to a contamination from a nonblazar population, amounts to 75 ± 4 sources. In conclusion, an excess of soft-spectra sources stands out in the low-latitude blazar sample. It can be attributed to a contamination from another population, likely of Galactic origin. After subtracting this  \nFigure 16. Top: photon-index distribution of high-latitude unassociated sources compared to that of blazars of unknown type. Bottom: photon-index distribution of low-latitude unassociated sources compared to that of young pulsars.  \n<!-- image -->  \nexcess from the number of reported blazars, the remaining number of detected blazars agrees reasonably well with that extrapolated from the high-latitude population (324 vs. 340 ± 20). Few blazars with missed associations are thus expected to contribute to the low-latitude unassociated sample.  \nUnassociated sources around the Galactic plane (referred to as GUs) account for a growing fraction of the detected sources, from 12% in DR1 to 24% in the new sample. The features presented here were already present in DR1 but have become more salient with the larger population.  \nThe GU latitude distribution, shown in Figure 18, shows a narrow Gaussian-like component ( glyph[similarequal] 210 sources) centered on the plane, and two broad shoulders ( glyph[similarequal] 760 sources) extending out to about | b | = 10 · (these components will be referred to as 'spike' and 'shoulders' respectively). The spike is narrower ( σ sin( b ) =0.005) than the young-pulsar distribution ( σ sin( b ) =0.012 for sources with sin( b ) < 0 . 03). The loci of the unassociated sources close to the plane exhibit a notable clustering. Some high-density regions ('hot spots') are readily seen in Figure 19, while other areas  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \nFigure 17. Variability-index (left) and photon-index (right) distributions of high-latitude unassociated sources and blazars for different bins in Test Statistic.  \n<!-- image -->  \nFigure 18. Galactic-latitude distribution of unassociated sources around the Galactic plane compared to that of pulsars.  \n<!-- image -->  \nare mostly devoid of unassociated sources. The clustering aspect also manifests itself in the large fraction (30%) of GUs with Flag 5 (i.e., close to a brighter neighbor). High-density regions have been determined using a  \n<!-- image -->  \nHEALPix tesselation with N side = 8 produced in Galactic coordinates. The average density of GUs per pixel in the three pixel rows spanning the Galactic Plane is about 7. We define the high-density pixels as those encompassing more than 14 GUs, i.e. twice the average value.  \nThe photon-index distribution is compared to that of PSR in the lower panel of Figure 16. The distribution is remarkably narrow and contrasts sharply with that observed at higher latitude (Figure 16 top). It peaks around Γ = 2 . 5, with a tail toward lower Γ values. This tail might be due to blazars, SNRs, or PWNe (Figure 14). The high-Γ end of the distribution is peculiar. No established class of Galactic γ -ray emitters exhibits such a soft component. We will refer to sources in this soft component as SGUs.  \nA dedicated flag (Flag 14) has been devised to indicate SGUs that are in regions of relatively high source density as seen in Figure 19. Flag 14 is set for SGUs lying inside a high-density pixel, having Γ > 2 . 4, TS < 500, and with a closest neighbor (classified either as GU, SPP, or UNK) less than 2 · away. Since most SPP or UNK sources probably belong to the same population as GUs as ascertained by their similar photon-index distributions (peaking around Γ=2.5), they are also flagged if the above conditions are filled. Flag 14 is set for 446 SGUs, 54 SPPs, and 49 UNKs. 16 SGUs were already significantly present in DR1, but their number in highdensity regions has now increased by 50% (see census of Flag 14 in Table 4).  \nReturning to the spectral properties, the photon-index distribution is essentially constant across the plane, as illustrated in Figure 20 for different (not mutually exclusive) sky regions. In this figure, 'clusters' refers to the high-density pixels mentioned above. The unassociated sources preferably lie within the pulsar branch in the LP TS curv vs. TS plane (Figure 21), supporting their spectral resemblance with pulsars. Their distribution of the LP curvature parameter β is also consistent with that of pulsars. Despite these similarities, other factors distinguish GUs and pulsars. First, the photon-index distributions do not closely match. The significant clustering observed also disfavors the possibility that GUs are mostly pulsars. Moreover, the lowb unassociated sources outnumber the detected PSR by about a factor of 9, making the scenario that so many pulsars have gone undetected so far implausible. The final blow to this possibility comes from their broad latitude distribution,  \n16 We checked that the set of flagged sources depends moderately (typical > 85% overlap) on the centering of the HEALPix pixels.  \nFigure 19. Positions of unassociated sources (crosses) around the Galactic center. The background is a count map of 83-228 MeV patch photons simulated over 12 yr with a pixel size of 1 · .  \n<!-- image -->  \nFigure 21. Same as Figure 15 bottom, but for low-latitude unassociated sources.  \n<!-- image -->  \nFigure 20. Photon-index distributions of low-latitude unassociated sources in different sky regions.  \n<!-- image -->  \nwhich does not match that of PSR (nor any of the distributions of known γ -ray emitters for that matter). Older pulsars ( > 10 6 yr), which could have drifted farther off the plane over time, exhibit a wider Galactic-latitude distribution, more compatible with that observed. However, these pulsars should be undetectable in γ -rays as they would have crossed the 'pair-creation death line' (see Abdo et al. 2013).  \nMost GUs exhibit a low source-to-background ratio (30% of them having the corresponding Flag 4 set). Whether they originate from mismodeled Galactic diffuse emission becomes a valid question. This emission is not fully understood, and its modeling suffers many uncertainties. To produce the Pass 8 model 17 , the contributions of different components, H i , CO, dark neutral  \nmedium, and inverse Compton, whose properties are established at larger wavelengths, were fitted to the data in 9 rings of increasing Galactocentric distances. An ad hoc component called the patch encapsulates the contributions of nontemplate components like the Galactic center excess, Loop I, or the Fermi bubbles. Other, yetto-be-identified large scale structures may be contributing (see J'ohannesson & Porter 2021, for a search performed at high latitude), and a comprehensive analysis of the patch still needs to be carried out. Understanding the patch low-energy behavior would be of particular relevance to this work. The patch covers a limited sky area and was smoothed out not to absorb features with spatial scales less than 4 · . Looking for spatial correlations between the SGUs and features present in the different components, the best match was obtained with the patch component at low energy (see Figure 19). The SGU latitude distribution exhibits a north/south sym-  \nry and a relative extension (Figure 18) that find a good correspondence in the patch at | b | < 10 · . Comprehending the photon excess underlying the patch will certainly be key in unveiling the nature of the SGUs.  \nBelow we list some possibilities regarding the SGUs and lay out avenues for further exploring them.  \n- · As mentioned above, SGUs may result from mismodeled Galactic diffuse emission. Allowing more flexibility to the model fitting by leaving free individual diffuse components, like H i and CO, would constrain the robustness of the SGU detections, if the templates used for the model are accurate on scales that are finer than the RoI size, while inaccurate on larger scales.\n- · A possibility to be considered is that SGUs trace the large-scale regions associated with fresh injection of cosmic-rays, leading to locally harder γ -ray emission than assumed in the model components. A worthwhile test would be to investigate (e.g., by means of simulations) whether the spurious point-like sources spawned from a possible, not-accounted-for diffuse excess would share common spectral features with the observed SGUs.\n- · Some SGUs could be spatially extended if they are related to structures in the diffuse emission. This would introduce some spectral distortion when analyzing them under the assumption that they are point-like, making them softer than they really are.\n- · Some correlation with the general structures of interstellar clouds is found by visual inspection, suggesting that the emission ascribed to SGUs could be due to molecular gas unaccounted for on small scales, e.g., because of saturated CO lines. This situation is reminiscent of that encountered in the EGRET era, which led to the discovery of 'dark gas' (Grenier et al. 2005) and the suppression of about 30% of unassociated EGRET sources (Casandjian & Grenier 2008) once a modified diffuse emission model was considered. Assessing the effect of saturated CO lines using 13 COdata would shed light on this possibility.\n- · Alternative ways could be pursued to build the patch (diffuse emission that is not accounted for by standard components of the interstellar emission model) and distinguish it from point sources (e.g., Selig et al. 2015).\n- · The alternative to mismodeled diffuse emission, namely the existence of an abundant, entirely  \nnew class of sources remains highly speculative. Follow-up multiwavelength studies (as in Bruzewski et al. 2021), e.g., of regions of high SGU densities or surrounding the brightest SGUs, could help elucidate this issue.", '7. CONCLUSIONS': "We built two incremental versions (DR2 after 10 yr and DR3 after 12 yr) of the fourth Fermi -LAT source catalog, the deepest-yet in the GeV energy range. Incremental versions preserve the existence and positions of sources in the previous versions (with a few exceptions, see § 3.2), and add new sources. We will continue providing incremental versions every two years, until a major ingredient (the event reconstruction, the calibrations, or the underlying diffuse emission model) improves significantly.  \nWe applied a new method to check the quality of the global model of the sky (see § 3.6). That method indicates that the interstellar emission model is too bright in a number of local molecular clouds. Conversely, it may be too faint close to the Galactic plane where we find many clustered soft unassociated sources (see § 6.2). Understanding the nature of these sources and what fraction are actually diffuse emission will be an important challenge, which could potentially lead to major advances like the discovery of the dark gas in the EGRET data (Grenier et al. 2005).  \nThe 4FGL-DR3 catalog includes 6658 sources, 1695 of which (25.5%) are found to be significantly variable on one-year timescales, and 3388 of which (50.9%) are fit with curved spectral shapes (see § 3.4). We mark 317 (4.8%) of the sources as potentially related to imperfections in the model for Galactic diffuse emission; the character c (for 'on top of interstellar gas clump') is appended to their names (except those already marked as e for extended). An additional 1836 (27.6%) are flagged in the catalog for less serious concerns, e.g., for depending sensitively on the details of the analysis or for being close to a brighter source (see § 3.5).  \nOf the 6658 sources in the catalog, 389 (5.8%) are considered identified, based on pulsations, correlated variability, or correlated angular sizes with observations at other wavelengths. We find likely lower-energy counterparts for 4112 other sources (61.8%). The remaining 2157 sources (32.4%) are unassociated.  \nThe identified and associated sources in the 4FGLDR3 catalog include many Galactic and extragalactic source classes (see § 5). The largest Galactic source class continues to be pulsars, with 135 young, 120 millisecond γ -ray pulsars, and 37 associations to non-LAT pulsars (mostly MSP). Other Galactic source classes have con-  \nnued to grow; 35 GLCs, 43 SNRs and 19 PWNe are now associated with LAT sources. Blazars remain the largest class of extragalactic sources, with 2251 identified or associated with BL Lac objects or FSRQ active galaxies and 1492 with unclassified blazars.  \nThe Fermi LAT Collaboration acknowledges generous ongoing support from a number of agencies and institutes that have supported both the development and the operation of the LAT as well as scientific data analysis. These include the National Aeronautics and Space Administration and the Department of Energy in the United States, the Commissariat 'a l'Energie Atomique and the Centre National de la Recherche Scientifique / Institut National de Physique Nucl'eaire et de Physique des Particules in France, the Agenzia Spaziale Italiana and the Istituto Nazionale di Fisica Nucleare in Italy, the Ministry of Education, Culture, Sports, Science and Technology (MEXT), High Energy Accelerator Research Organization (KEK) and Japan Aerospace Exploration Agency (JAXA) in Japan, and the K. A. Wallenberg Foundation, the Swedish Research Council and the Swedish National Space Board in Sweden.  \nAdditional support for science analysis during the operations phase is gratefully acknowledged from the Istituto Nazionale di Astrofisica in Italy and the Centre National d' ' Etudes Spatiales in France. This work performed in part under DOE Contract DE-AC0276SF00515.  \nThis work made extensive use of the ATNF pulsar catalog 18 (Manchester et al. 2005). This research has made use of the NASA/IPAC Extragalactic Database (NED) (2019) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration, and of archival data, software, and online services provided by the ASI Science Data Center (ASDC) operated by the Italian Space Agency. We used the Manitoba SNR catalog (Ferrand & Safi-Harb 2012) to check recently published extended sources. We acknowledge the Einstein@Home project for providing new pulsar associations through the dedicated efforts of the Einstein@Home volunteers. The Einstein@Home project is supported by the NSF award 1816904.  \nSoftware: Gardian (Ackermann et al. 2012), GALPROP 19 (Porter et al. 2017), HEALPix 20 (G'orski et al. 2005), Aladin 21 , TOPCAT 22 (Taylor 2005)  \nFacility:  \nFermi -LAT", 'REFERENCES': "Abdo, A. A., Ackermann, M., Ajello, M., et al. 2010, A&A, 523, A46, doi: 10.1051/0004-6361/201014855  \n-. 2011, ApJL, 736, L11,  \ndoi: 10.1088/2041-8205/736/1/L11  \nAbdo, A. A., Ajello, M., Allafort, A., et al. 2013, ApJS,  \n208, 17, doi: 10.1088/0067-0049/208/2/17  \n- 19 http://galprop.stanford.edu\n- 20 http://healpix.jpl.nasa.gov/\n- 21 http://aladin.u-strasbg.fr/\n- 22 http://www.star.bristol.ac.uk/ ∼ mbt/topcat/  \nCorbet, R. H. D., Chomiuk, L., Coe, M. J., et al. 2019, ApJ, 884, 93, doi: 10.3847/1538-4357/ab3e32 Coti Zelati, F., Papitto, A., de Martino, D., et al. 2019, A&A, 622, A211, doi: 10.1051/0004-6361/201834835 Coti Zelati, F., Hugo, B., Torres, D. F., et al. 2021, A&A, 655, A52, doi: 10.1051/0004-6361/202141431 Desai, A., Marchesi, S., Rajagopal, M., & Ajello, M. 2019, ApJS, 241, 5, doi: 10.3847/1538-4365/ab01fc Devin, J., Lemoine-Goumard, M., Grondin, M. H., et al. 2020, A&A, 643, A28, doi: 10.1051/0004-6361/202038503 Esposito, P., Israel, G. L., Sidoli, L., et al. 2014, MNRAS, 441, 1126, doi: 10.1093/mnras/stu659 Ferrand, G., & Safi-Harb, S. 2012, Advances in Space Research, 49, 1313, doi: 10.1016/j.asr.2012.02.004 Fortin, F., Chaty, S., Coleiro, A., Tomsick, J. A., & Nitschelm, C. H. R. 2018, A&A, 618, A150, doi: 10.1051/0004-6361/201731265 G'orski, K. M., Hivon, E., Banday, A. J., et al. 2005, ApJ, 622, 759, doi: 10.1086/427976 Grenier, I. A., Casandjian, J.-M., & Terrier, R. 2005, Science, 307, 1292, doi: 10.1126/science.1106924 Grondin, M.-H., Funk, S., Lemoine-Goumard, M., et al. 2011, ApJ, 738, 42, doi: 10.1088/0004-637X/738/1/42 Halpern, J. P., Strader, J., & Li, M. 2017, ApJ, 844, 150, doi: 10.3847/1538-4357/aa7cff Harding, A. K., Usov, V. V., & Muslimov, A. G. 2005, ApJ, 622, 531, doi: 10.1086/427840 Jarvela, E., Berton, M., Ciroi, S., et al. 2020, A&A, 636, L12, doi: 10.1051/0004-6361/202037826 Jean, P., Cheung, C. C., Ojha, R., van Zyl, P., & Angioni, R. 2018, ATel, 11546, 1 J'ohannesson, G., & Porter, T. A. 2021, ApJ, 917, 30, doi: 10.3847/1538-4357/ac01c9 Johnson, T. J., Wood, K. S., Kerr, M., et al. 2018, ApJ, 863, 27, doi: 10.3847/1538-4357/aad185 Kennedy, M. R., Breton, R. P., Clark, C. J., et al. 2020, MNRAS, doi: 10.1093/mnras/staa912 Kerr, M. 2010, PhD thesis, University of Washington, ArXiv:1101.6072. https://arxiv.org/abs/1101.6072 Li, K.-L., Jane Yap, Y. X., Hui, C. Y., & Kong, A. K. H. 2021, ApJ, 911, 92, doi: 10.3847/1538-4357/abeb76 Li, K.-L., Kong, A., Aydi, E., et al. 2020, ATel, 13868, 1 Li, K.-L., Kong, A. K. H., Hou, X., et al. 2016, ApJ, 833, 143, doi: 10.3847/1538-4357/833/2/143 Manchester, R. N., Hobbs, G. B., Teoh, A., & Hobbs, M. 2005, AJ, 129, 1993, doi: 10.1086/428488 Mares, A., Lemoine-Goumard, M., Acero, F., et al. 2021, ApJ, 912, 158, doi: 10.3847/1538-4357/abef62 Mart'ı, J., S'anchez-Ayaso, E., Luque-Escamilla, P. L., et al. 2020, MNRAS, 492, 4291, doi: 10.1093/mnras/staa072", 'A. DESCRIPTION OF THE FITS VERSION OF THE 4FGL-DR3 CATALOG': 'Table 6 . LAT 4FGL-DR3 FITS format: LAT Point Source Catalog extension. Changes since 4FGL appear in Italics . The file is available at https://fermi.gsfc.nasa.gov/ssc/data/access/lat/12yr catalog/.Table 6 continued', 'Table 6 (continued)': 'Table 6 continued', '4FGL DR3': "Table 6 (continued)  \nThe FITS format version of the third release of the 4FGL catalog has seven binary table extensions. The extension LAT Point Source Catalog Extension has all of the information about the sources, and has changed to adapt to the evolutions described in § 3.4. Its format is described in Table 6. Columns (Unc )PLEC Index and (Unc )PLEC Expfactor have been replaced by (Unc )PLEC IndexS and (Unc )PLEC ExpfactorS corresponding to the new pulsar parameterization (Eq. 2). The table has 6659 rows for 6658 sources because the Crab nebula is described by two entries (the synchrotron component and the inverse Compton component) but counted as only one source. The Crab pulsar is another entry and counted as a separate source. The extended sources are indicated by an 'e' appended to their names in the main table.  \nThe extensions GTI , ExtendedSources , ROIs , Components , EnergyBounds , Hist Start describe auxiliary details and have the same format as in 4FGL."}