We present an incremental version (4FGL-DR3, 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 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.

The fourth catalog of active galactic nuclei (AGNs) detected by the Fermi Gamma-ray Space Telescope Large Area Telescope (4LAC) between 2008 August 4 and 2016 August 2 contains objects located at high Galactic latitudes ( ). It includes 85% more sources than the previous 3LAC catalog based on 4 yr of data. AGNs represent at least 79% of the high-latitude sources in the fourth Fermi-Large Area Telescope Source Catalog (4FGL), which covers the energy range from 50 MeV to 1 TeV. In addition, gamma-ray AGNs are found at low Galactic latitudes. Most of the 4LAC AGNs are blazars (98%), while the remainder are other types of AGNs. The blazar population consists of 24% Flat Spectrum Radio Quasars (FSRQs), 38% BL Lac-type objects, and 38% blazar candidates of unknown types (BCUs). On average, FSRQs display softer spectra and stronger variability in the gamma-ray band than BL Lacs do, confirming previous findings. All AGNs detected by ground-based atmospheric Cerenkov telescopes are also found in the 4LAC.

We present a search for spatial extension in high-latitude ( ) sources in recent Fermi point source catalogs. The result is the Fermi High-Latitude Extended Sources Catalog, which provides source extensions (or upper limits thereof) and likelihood profiles for a suite of tested source morphologies. We find 24 extended sources, 19 of which were not previously characterized as extended. These include sources that are potentially associated with supernova remnants and star-forming regions. We also found extended γ-ray emission in the vicinity of the Cen A radio lobes and—at GeV energies for the first time—spatially coincident with the radio emission of the SNR CTA 1, as well as from the Crab Nebula. We also searched for halos around active galactic nuclei, which are predicted from electromagnetic cascades induced by the e ⁺ e ⁻ pairs that are deflected in intergalactic magnetic fields. These pairs are produced when γ-rays interact with background radiation fields. We do not find evidence for extension in individual sources or in stacked source samples. This enables us to place limits on the flux of the extended source components, which are then used to constrain the intergalactic magnetic field to be stronger than 3 × 10⁻¹⁶ G for a coherence length λ ≳ 10 kpc, even when conservative assumptions on the source duty cycle are made. This improves previous limits by several orders of magnitude.

A neutrino with energy ∼290 TeV, IceCube-170922A, was detected in coincidence with the BL Lac object TXS 0506+056 during enhanced gamma-ray activity, with chance coincidence being rejected at ∼3 σ level. We monitored the object in the very-high-energy (VHE) band with the Major Atmospheric Gamma-ray Imaging Cherenkov (MAGIC) telescopes for ∼41 hr from 1.3 to 40.4 days after the neutrino detection. Day-timescale variability is clearly resolved. We interpret the quasi-simultaneous neutrino and broadband electromagnetic observations with a novel one-zone lepto-hadronic model, based on interactions of electrons and protons co-accelerated in the jet with external photons originating from a slow-moving plasma sheath surrounding the faster jet spine. We can reproduce the multiwavelength spectra of TXS 0506+056 with neutrino rate and energy compatible with IceCube-170922A, and with plausible values for the jet power of . The steep spectrum observed by MAGIC is concordant with internal γγ absorption above ∼100 GeV entailed by photohadronic production of a ∼290 TeV neutrino, corroborating a genuine connection between the multi-messenger signals. In contrast to previous predictions of predominantly hadronic emission from neutrino sources, the gamma-rays can be mostly ascribed to inverse Compton upscattering of external photons by accelerated electrons. The X-ray and VHE bands provide crucial constraints on the emission from both accelerated electrons and protons. We infer that the maximum energy of protons in the jet comoving frame can be in the range ∼10 ¹⁴ – 10 ¹⁸ eV.

We present 294 pulsars found in GeV data from the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope. Another 33 millisecond pulsars (MSPs) discovered in deep radio searches of LAT sources will likely reveal pulsations once phase-connected rotation ephemerides are achieved. A further dozen optical and/or X-ray binary systems colocated with LAT sources also likely harbor gamma-ray MSPs. This catalog thus reports roughly 340 gamma-ray pulsars and candidates, 10% of all known pulsars, compared to ≤11 known before Fermi. Half of the gamma-ray pulsars are young. Of these, the half that are undetected in radio have a broader Galactic latitude distribution than the young radio-loud pulsars. The others are MSPs, with six undetected in radio. Overall, ≥236 are bright enough above 50 MeV to fit the pulse profile, the energy spectrum, or both. For the common two-peaked profiles, the gamma-ray peak closest to the magnetic pole crossing generally has a softer spectrum. The spectral energy distributions tend to narrow as the spindown power E ̇ decreases to its observed minimum near 10³³ erg s⁻¹, approaching the shape for synchrotron radiation from monoenergetic electrons. We calculate gamma-ray luminosities when distances are available. Our all-sky gamma-ray sensitivity map is useful for population syntheses. The electronic catalog version provides gamma-ray pulsar ephemerides, properties, and fit results to guide and be compared with modeling results.

We report the discovery of GRB 221009A, the highest flux gamma-ray burst (GRB) ever observed by the Fermi Gamma-ray Burst Monitor (Fermi-GBM). This GRB has continuous prompt emission lasting more than 600 s, which smoothly transitions to afterglow emission visible in the Fermi-GBM energy range (8 keV–40 MeV), and total energetics higher than any other burst in the Fermi-GBM sample. By using a variety of new and existing analysis techniques we probe the spectral and temporal evolution of GRB 221009A. We find no emission prior to the Fermi-GBM trigger time (t ₀; 2022 October 9 at 13:16:59.99 UTC), indicating that this is the time of prompt emission onset. The triggering pulse exhibits distinct spectral and temporal properties suggestive of the thermal, photospheric emission of shock breakout, with significant emission up to ∼15 MeV. We characterize the onset of external shock at t ₀ + 600 s and find evidence of a plateau region in the early-afterglow phase, which transitions to a slope consistent with Swift-XRT afterglow measurements. We place the total energetics of GRB 221009A in context with the rest of the Fermi-GBM sample and find that this GRB has the highest total isotropic-equivalent energy (E _γ,iso = 1.0 × 10⁵⁵ erg) and second highest isotropic-equivalent luminosity (L _γ,iso = 9.9 × 10⁵³ erg s^–1) based on its redshift of z = 0.151. These extreme energetics are what allowed us to observe the continuously emitting central engine of Fermi-GBM from the beginning of the prompt emission phase through the onset of early afterglow.