We present a complete analysis of Fermi Large Area Telescope (LAT) data of GRB 221009A, the brightest gamma-ray burst (GRB) ever detected. The burst emission above 30 MeV detected by the LAT preceded, by 1 s, the low-energy (<10 MeV) pulse that triggered the Fermi Gamma-Ray Burst Monitor (GBM), as has been observed in other GRBs. The prompt phase of GRB 221009A lasted a few hundred seconds. It was so bright that we identify a bad time interval of 64 s caused by the extremely high flux of hard X-rays and soft gamma rays, during which the event reconstruction efficiency was poor and the dead time fraction quite high. The late-time emission decayed as a power law, but the extrapolation of the late-time emission during the first 450 s suggests that the afterglow started during the prompt emission. We also found that high-energy events observed by the LAT are incompatible with synchrotron origin, and, during the prompt emission, are more likely related to an extra component identified as synchrotron self-Compton (SSC). A remarkable 400 GeV photon, detected by the LAT 33 ks after the GBM trigger and directionally consistent with the location of GRB 221009A, is hard to explain as a product of SSC or TeV electromagnetic cascades, and the process responsible for its origin is uncertain. Because of its proximity and energetic nature, GRB 221009A is an extremely rare event.

A 2.1-year periodic oscillation of the gamma-ray flux from the blazar PG 1553+113 has previously been tentatively identified in almost 7 year of data from the Fermi Large Area Telescope. After 15 years of Fermi sky-survey observations, doubling the total time range, we report >7 cycle gamma-ray modulation with an estimated significance of 4 sigma against stochastic red noise. Independent determinations of oscillation period and phase in the earlier and the new data are in close agreement (chance probability <0.01). Pulse timing over the full light curve is also consistent with a coherent periodicity. Multiwavelength new data from Swift X-Ray Telescope, Burst Alert Telescope, and UVOT, and from KAIT, Catalina Sky Survey, All-Sky Automated Survey for Supernovae, and Owens Valley Radio Observatory ground-based observatories as well as archival Rossi X-Ray Timing Explorer satellite-All Sky Monitor data, published optical data of Tuorla, and optical historical Harvard plates data are included in our work. Optical and radio light curves show clear correlations with the gamma-ray modulation, possibly with a nonconstant time lag for the radio flux. We interpret the gamma-ray periodicity as possibly arising from a pulsational accretion flow in a sub-parsec binary supermassive black hole system of elevated mass ratio, with orbital modulation of the supplied material and energy in the jet. Other astrophysical scenarios introduced include instabilities, disk and jet precession, rotation or nutation, and perturbations by massive stars or intermediate-mass black holes in polar orbit.

The nearby elliptical galaxy M87 contains one of the only two supermassive black holes whose emission surrounding the event horizon has been imaged by the Event Horizon Telescope (EHT). In 2018, more than two dozen multi-wavelength (MWL) facilities (from radio to gamma-ray energies) took part in the second M87 EHT campaign. The goal of this extensive MWL campaign was to better understand the physics of the accreting black hole M87*, the relationship between the inflow and inner jets, and the high-energy particle acceleration. Understanding the complex astrophysics is also a necessary first step towards performing further tests of general relativity. The MWL campaign took place in April 2018, overlapping with the EHT M87* observations. We present a new, contemporaneous spectral energy distribution (SED) ranging from radio to very high energy (VHE) gamma-rays, as well as details of the individual observations and light curves. We also conduct phenomenological modelling to investigate the basic source properties. We present the first VHE gamma-ray flare from M87 detected since 2010. The flux above 350 GeV has more than doubled within a period of about 36 hours. We find that the X-ray flux is enhanced by about a factor of two compared to 2017, while the radio and millimetre core fluxes are consistent between 2017 and 2018. We detect evidence for a monotonically increasing jet position angle that corresponds to variations in the bright spot of the EHT image. Our results show the value of continued MWL monitoring together with precision imaging for addressing the origins of high-energy particle acceleration. While we cannot currently pinpoint the precise location where such acceleration takes place, the new VHE gamma-ray flare already presents a challenge to simple one-zone leptonic emission model approaches, and emphasises the need for combined image and spectral modelling.

The recent detection of extended γ-ray emission around middle-aged pulsars is interpreted as inverse-Compton scattering of ambient photons by electron-positron pairs escaping the pulsar wind nebula, which are confined near the system by unclear mechanisms. This emerging population of γ-ray sources was first discovered at TeV energies and remains underexplored in the GeV range. To address this, we conducted a systematic search for extended sources along the Galactic plane using 14 years of Fermi-LAT data above 10 GeV, aiming to identify a number of pulsar halo candidates and extend our view to lower energies. The search covered the inner Galactic plane (| l|≤ 100^∘, | b|≤ 1^∘) and the positions of known TeV sources and bright pulsars, yielding broader astrophysical interest. We found 40 such sources, forming the Second Fermi Galactic Extended Sources Catalog (2FGES), most with 68 1.0^∘ and relatively hard spectra with photon indices below 2.5. We assessed detection robustness using field-specific alternative interstellar emission models and by inspecting significance maps. Noting 13 sources previously known as extended in the 4FGL-DR3 catalog and five dubious sources from complex regions, we report 22 newly detected extended sources above 10 GeV. Of these, 13 coincide with H.E.S.S., HAWC, or LHAASO sources; six coincide with bright pulsars (including four also coincident with TeV sources); six are associated with 4FGL point sources only; and one has no association in the scanned catalogs. Notably, six to eight sources may be related to pulsars as classical pulsar wind nebulae or pulsar halos.

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 33 erg s −1 , 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 present an incremental version (4FGL-DR4, for Data Release 4) of the fourth Fermi-LAT catalog containing 7194 gamma-ray sources. Based on the first 14 years of science data in the energy range from 50 MeV to 1 TeV, it uses the same analysis methods as the 4FGL-DR3 catalog did for 12 years of data, with only a few improvements. The spectral parameters, spectral energy distributions, light curves and associations are updated for all sources. We add four new extended sources and modify two existing ones. Among the 6658 4FGL-DR3 sources, we delete 14 and change the localization of 10, while 28 are newly associated and two associations were changed. We add 546 point sources, among which 8 are considered identified and 228 have a plausible counterpart at other wavelengths. Most are just above the detection threshold, and 14 are transient sources below the detection threshold that can affect the light curves of nearby sources.

A measurement of MeV-GeV linear polarization would provide information about the emission mechanisms of gamma-ray sources, especially pulsars and AGN (Active Galaxy Nuclei). The $\textit{Fermi}$ Large Area Telescope ($\textit{Fermi}$-LAT) detects gamma rays between $20\,MeV$ and $\sim 300\,GeV$, using their conversion into an electron-positron pair, with the linear polarization affecting the azimuthal direction of the pair. Measuring this azimuthal direction is extremely challenging due to the multiple scattering of the electron-positron pair in the detector. We present the current effort to extract the polarization information from the $\textit{Fermi}$-LAT data, including recent developments in the simulation implementation and the event reconstruction and selection. This allows us to estimate the expected sensitivity to linear polarization of the $\textit{Fermi}$-LAT, and we apply this to the Vela pulsar.

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 0 ; 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 0 + 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 55 erg) and second highest isotropic-equivalent luminosity ( L γ ,iso = 9.9 × 10 53 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.

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 co-located 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 $\leq 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 6 undetected in radio. Overall, >235 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 $\dot E$ decreases to its observed minimum near $10^{33}$ erg s$^{-1}$, 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.

Reliable neutron star mass measurements are key to determining the equation-of-state of cold nuclear matter, but these are rare. "Black Widows" and "Redbacks" are compact binaries consisting of millisecond pulsars and semi-degenerate companion stars. Spectroscopy of the optically bright companions can determine their radial velocities, providing inclination-dependent pulsar mass estimates. While inclinations can be inferred from subtle features in optical light curves, such estimates may be systematically biased due to incomplete heating models and poorly-understood variability. Using data from the Fermi Large Area Telescope, we have searched for gamma-ray eclipses from 49 spider systems, discovering significant eclipses in 7 systems, including the prototypical black widow PSR B1957$+$20. Gamma-ray eclipses require direct occultation of the pulsar by the companion, and so the detection, or significant exclusion, of a gamma-ray eclipse strictly limits the binary inclination angle, providing new robust, model-independent pulsar mass constraints. For PSR B1957$+$20, the eclipse implies a much lighter pulsar ($M_{\rm psr} = 1.81 \pm 0.07\,M_{\odot}$) than inferred from optical light curve modelling.

On October 9, 2022, the brightest GRB ever detected, also called the "Brightest Of All Time" (or BOAT), was detected by the Fermi Large Area Telescope (LAT), providing an unprecedented rich dataset that will likely help to understand the emission mechanisms underlying this class of extra-galactic objects. Some of the unique features of GRB221009A include: 1) high-energy (up to 100 MeV) emission observed during the first emission episode (precursor) at the GBM trigger time 𝑇 0 ; 2) Several GeV events, including a 100 GeV photon, observed during the brightest emission episode; and 3) a very long-lived high-energy extended emission phase covering tens of kiloseconds, with a peculiar 400 GeV photon observed at almost 10 hours after 𝑇 0 . Here we review all these observations, focusing especially on the efforts that the Fermi LAT collaboration made in order to understand and recover part of the data affected by the extremely high soft gamma-ray noise during the brightest emission episode at 𝑇 0 + 200 s.

The Fermi Large Area Telescope (LAT) light curve repository (LCR) is a publicly available, continually updated library of gamma-ray light curves of variable Fermi-LAT sources generated over multiple timescales. The Fermi-LAT LCR aims to provide publication-quality light curves binned on timescales of 3 days, 7 days, and 30 days for 1525 sources deemed variable in the source catalog of the first 10 years of Fermi-LAT observations. The repository consists of light curves generated through full likelihood analyses that model the sources and the surrounding region, providing fluxes and photon indices for each time bin. The LCR is intended as a resource for the time-domain and multi-messenger communities by allowing users to quickly search LAT data to identify correlated variability and flaring emission episodes from gamma-ray sources. We describe the sample selection and analysis employed by the LCR and provide an overview of the associated data access portal.

The third release of the fourth general catalog of Fermi-LAT sources (4FGL-DR3), based on 12 years of data between 50 MeV and 1 TeV, is presented. Improvements in the analysis method relative to the original 4FGL catalog and new features are reviewed. The 4FGL-DR3 includes similar or equal to 900 more sources than the previous release (4FGL-DR2, obtained with 10 years of data) and similar or equal to 1500 more sources than the 4FGL. About 43% of the new sources are associated with counterparts at other wavelengths, which are mostly blazar candidates. The properties of the global set of unassociated sources reported in the catalog are discussed, with particular attention to those lying close to the Galactic plane. A population of unassociated sources that do not fit in with already known classes of gamma-ray emitters is emphasized.

We present an incremental version (4FGL-DR3, for Data Release 3) of the fourth Fermi-LAT catalog of gamma-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). 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 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.

We report the discovery of PSR J1555-2908, a 1.79 ms radio and gamma-ray pulsar in a 5.6 hr binary system with a minimum companion mass of 0.052 M-circle dot. This fast and energetic (<(E) over dot> = 3 x 10(35) erg s(-1)) millisecond pulsar was first detected as a gamma-ray point source in Fermi Large Area Telescope (LAT) sky survey observations. Guided by a steep-spectrum radio point source in the Fermi error region, we performed a search at 820 MHz with the Green Bank Telescope that first discovered the pulsations. The initial radio pulse timing observations provided enough information to seed a search for gamma-ray pulsations in the LAT data, from which we derive a timing solution valid for the full Fermi mission. In addition to the discovery and timing of radio and gamma-ray pulsations, we searched for X-ray pulsations using NICER but no significant pulsations were detected. We also obtained time-series r-band photometry that indicates strong heating of the companion star by the pulsar wind. Material blown off the heated companion eclipses the 820 MHz radio pulse during inferior conjunction of the companion for approximate to 10% of the orbit, which is twice the angle subtended by its Roche lobe in an edge-on system.

After large galaxies merge, their central supermassive black holes are expected to form binary systems. Their orbital motion should generate a gravitational wave background (GWB) at nanohertz frequencies. Searches for this background use pulsar timing arrays, which perform long-term monitoring of millisecond pulsars at radio wavelengths. We used 12.5 years of Fermi Large Area Telescope data to form a gamma-ray pulsar timing array. Results from 35 bright gamma-ray pulsars place a 95% credible limit on the GWB characteristic strain of 1.0 x 10(-14) at a frequency of 1 year(-1). The sensitivity is expected to scale with t(obs), the observing time span, as t(obs)(-13/6). This direct measurement provides an independent probe of the GWB while offering a check on radio noise models.

An incremental version of the fourth catalog of active galactic nuclei (AGNs) detected by the Fermi-Large Area Telescope is presented. This version (4LAC-DR3) derives from the third data release of the 4FGL catalog based on 12 years of E>50 MeV gamma-ray data, where the spectral parameters, spectral energy distributions (SEDs), yearly light curves, and associations have been updated for all sources. The new reported AGNs include 587 blazar candidates and four radio galaxies. We describe the properties of the new sample and outline changes affecting the previously published one. We also introduce two new parameters in this release, namely the peak energy of the SED high-energy component and the corresponding flux. These parameters allow an assessment of the Compton dominance, the ratio of the Inverse-Compton to the synchrotron peak luminosities, without relying on X-ray data.

The 559 Hz black-widow pulsar PSR J1555−2908, originally discovered in radio, is also a bright gamma-ray pulsar. Timing its pulsations using 12 yr of Fermi-Large Area Telescope gamma-ray data reveals long-term variations in its spin frequency that are much larger than is observed from other millisecond pulsars. While this variability in the pulsar rotation rate could be intrinsic “timing noise,” here we consider an alternative explanation: the variations arise from the presence of a very-low-mass third object in a wide multiyear orbit around the neutron star and its low-mass companion. With current data, this hierarchical-triple-system model describes the pulsar’s rotation slightly more accurately than the best-fitting timing noise model. Future observations will show if this alternative explanation is correct.

The third release of the fourth general catalog of Fermi-LAT sources (4FGL-DR3), based on 12 years of data between 50 MeV and 1 TeV, is presented. Improvements in the analysis method relative to the original 4FGL catalog and new features are reviewed. The 4FGL-DR3 includes ~ 900 more sources than the previous release (4FGL-DR2, obtained with 10 years of data) and ~ 1500 more sources than the 4FGL. About 43% of the new sources are associated with counterparts at other wavelengths, which are mostly blazar candidates. The properties of the global set of unassociated sources reported in the catalog are discussed, with particular attention to those lying close to the Galactic plane. A population of unassociated sources that do not fit in with already known classes of gamma-ray emitters is emphasized.

We present the first Fermi - Large Area Telescope (LAT) solar flare catalog covering the 24 th solar cycle. This catalog contains 45 Fermi -LAT solar flares (FLSFs) with emission in the gamma-ray energy band (30 MeV - 10 GeV) detected with a significance greater than 5 sigma over the years 2010-2018. A subsample containing 37 of these flares exhibit delayed emission beyond the prompt-impulsive hard X-ray phase with 21 flares showing delayed emission lasting more than two hours. No prompt-impulsive emission is detected in four of these flares. We also present in this catalog the observations of GeV emission from 3 flares originating from Active Regions located behind the limb (BTL) of the visible solar disk. We report the light curves, spectra, best proton index and localization (when possible) for all the FLSFs. The gamma-ray spectra is consistent with the decay of pions produced by >300 MeV protons. This work contains the largest sample of high-energy gamma-ray flares ever reported and provides the unique opportunity to perform population studies on the different phases of the flare and thus allowing to open a new window in solar physics.