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.

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.

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.

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 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.

The uncharted territory of the high-redshift (z greater than or similar to 3) universe holds the key to understanding the evolution of quasars. In an attempt to identify the most extreme members of the quasar population, that is, blazars, we have carried out a multiwavelength study of a large sample of radio-loud quasars beyond z = 3. Our sample consists of nine gamma-ray-detected blazars and 133 candidate blazars selected based on the flatness of their soft X-ray spectra (0.3-10 keV photon index <= 1.75), including 15 with Nuclear Spectroscopic Telescope Array (NuSTAR).observations. The application of the likelihood profile stacking technique reveals that the high-redshift blazars are faint gamma-ray emitters with steep spectra. The high-redshift blazars host massive black holes (< log M-BH,M- M circle dot > > 9) and luminous accretion disks (< L-disk > > 10(46) erg s(-1)). Their broadband spectral energy distributions are found to be dominated by high-energy radiation, indicating their jets are among the most luminous ones. Focusing on the sources exhibiting resolved X-ray jets (as observed with the Chandra satellite), we find the bulk Lorentz factor to be larger with respect to other z > 3 blazars, indicating faster moving jets. We conclude that the presented list of high-redshift blazars may act as a reservoir for follow-up observations, such as with NuSTAR, to understand the evolution of relativistic jets at the dawn of the universe.

An incremental version (4LAC-DR2) of the fourth catalog of active galactic nuclei (AGNs) detected by the Fermi-LAT is presented. This version is associated with the second release of the 4FGL general catalog (based on 10 years of data), where the spectral parameters, spectral energy distributions, yearly light curves, and associations have been updated for all sources. The new reported AGNs include two radio galaxies and 283 blazars. We briefly describe the properties of the new sample and outline changes affecting the previously published sample. Note: Users of this incremental release are requested to cite the original 4LAC paper (Ajello M. et al., 2020, ApJ, 892, 105).

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 . 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.

e-ASTROGAM (enhanced ASTROGAM) is a breakthrough Observatory space mission, with a detector composed by a Silicon tracker, a calorimeter, and an anticoincidence system, dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV - the lower energy limit can be pushed to energies as low as 150 keV for the tracker, and to 30 keV for calorimetric detection. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LIGO-Virgo-GEO600-KAGRA, SKA, ALMA, E-ELT, TMT, LSST, JWST, Athena, CTA, IceCube, KM3NeT, and LISA.

e-ASTROGAM (enhanced ASTROGAM) is a breakthrough Observatory space mission, with a detector composed by a Silicon tracker, a calorimeter, and an anticoincidence system, dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV - the lower energy limit can be pushed to energies as low as 150 keV for the tracker, and to 30 keV for calorimetric detection. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LIGO-Virgo-GEO600-KAGRA, SKA, ALMA, E-ELT, TMT, LSST, JWST, Athena, CTA, IceCube, KM3NeT, and LISA.

We present multi-frequency simultaneous VLBI radio observations of the flat spectrum radio quasar 1633+382 (4C 38.41) as part of the interferometric monitoring of γ -ray-bright active galactic nuclei (iMOGABA) program combined with additional observations in the radio, optical, X-rays, and γ -rays carried out during the period 2012 March–2015 August. The monitoring of this source reveals a significant long-lived increase in its activity for approximately two years in the radio bands, which correlates with a similar increase in all other bands from submillimeter to γ -rays. A significant correlation is also found between radio fluxes and simultaneous spectral indices during this period. The study of the discrete correlation function indicates time lags smaller than the uncertainties of ∼40 days among both radio bands and high-energy bands, and a time lag of ∼70 days, with γ -rays leading radio emission. We interpret this as showing that the high-energy and radio fluxes arise from different emitting regions, located at 1 ± 13 and 40 ± 13 pc from the central engine respectively.

Gamma rays reveal the Universe's history How many stars have formed in the Universe, and when did they do so? These fundamental questions are difficult to answer because there are systematic uncertainties in converting the light we observe into the total mass of stars in galaxies. The Fermi-LAT Collaboration addressed these questions by exploiting the way that gamma rays from distant blazars propagate through intergalactic space, which depends on the total amount of light emitted by all galaxies. The collaboration found that star formation peaked about 3 billion years after the Big Bang (see the Perspective by Prandini). Although this is similar to previous estimates from optical and infrared observations, the results provide valuable confirmation because they should be affected by different systematic effects. Science , this issue p. 1031 ; see also p. 995

Context. The BL Lac object S5 0716+714, a highly variable blazar, underwent an impressive outburst in January 2015 (Phase A), followed by minor activity in February (Phase B). The MAGIC observations were triggered by the optical flux observed in Phase A, corresponding to the brightest ever reported state of the source in the R-band. Aims.The comprehensive dataset collected is investigated in order to shed light on the mechanism of the broadband emission. Methods. Multi-wavelength light curves have been studied together with the broadband spectral energy distributions (SEDs). The sample includes data from Effelsberg, OVRO, Metsähovi, VLBI, CARMA, IRAM, SMA, Swift-UVOT, KVA, Tuorla, Steward, RINGO3, KANATA, AZT-8+ST7, Perkins, LX-200, Swift-XRT, NuSTAR, Fermi-LAT and MAGIC. Results. The flaring state of Phase A was detected in all the energy bands, providing for the first time a multi-wavelength sample of simultaneous data from the radio band to the very-high-energy (VHE, E > 100 GeV). In the constructed SED, the Swift-XRT+NuSTAR data constrain the transition between the synchrotron and inverse Compton components very accurately, while the second peak is constrained from 0.1 GeV to 600 GeV by Fermi+MAGIC data. The broadband SED cannot be described with a one-zone synchrotron self-Compton model as it severely underestimates the optical flux in order to reproduce the X-ray to γ-ray data. Instead we use a two-zone model. The electric vector position angle (EVPA) shows an unprecedented fast rotation. An estimation of the redshift of the source by combined high-energy (HE, 0.1 GeV < E < 100 GeV) and VHE data provides a value of z = 0.31 ± 0.02stats ± 0.05sys, confirming the literature value. Conclusions. The data show the VHE emission originating in the entrance and exit of a superluminal knot in and out of a recollimation shock in the inner jet. A shock–shock interaction in the jet seems responsible for the observed flares and EVPA swing. This scenario is also consistent with the SED modeling.

We present a catalog of sources detected above 10 GeV by the Fermi Large Area Telescope (LAT) in the first 7 years of data using the Pass 8 event-level analysis. This is the Third Catalog of Hard Fermi-LAT Sources (3FHL), containing 1556 objects characterized in the 10 GeV-2 TeV energy range. The sensitivity and angular resolution are improved by factors of 3 and 2 relative to the previous LAT catalog at the same energies (1FHL). The vast majority of detected sources (79%) are associated with extragalactic counterparts at other wavelengths, including 16 sources located at very high redshift (z > 2). Of the sources, 8% have Galactic counterparts and 13% are unassociated (or associated with a source of unknown nature). The high-latitude sky and the Galactic plane are observed with a flux sensitivity of 4.4 to 9.5 x 10(-11) ph cm(-2) s(-1), respectively (this is approximately 0.5% and 1% of the Crab Nebula flux above 10 GeV). The catalog includes 214 new gamma-ray sources. The substantial increase in the number of photons (more than 4 times relative to 1FHL and 10 times to 2FHL) also allows us to measure significant spectral curvature for 32 sources and find flux variability for 163 of them. Furthermore, we estimate that for the same flux limit of 10(-12) erg cm(-2) s(-1), the energy range above 10 GeV has twice as many sources as the range above 50 GeV, highlighting the importance, for future Cherenkov telescopes, of lowering the energy threshold as much as possible.

We present multi–frequency simultaneous VLBI radio observations of the flat spectrum radio quasar 1633+382 (4C 38.41) as part of the interferometric monitoring of gamma-ray bright active galactic nuclei (iMOGABA) program combined with additional observations in radio, optical, X–rays and γ−rays carried out between the period 2012 March – 2015 August. The monitoring of this source reveals a significant long-lived increase in its activity since approximately two years in the radio bands, which correlates with a similar increase on all other bands from sub–millimeter to γ−rays. A significant correlation is also found between radio fluxes and simultaneous spectral indices during this period. The study of the discrete correlation function (DCF) indicates time lags smaller than the ∼ 40 days uncertainties among both radio bands and also high-energy bands, and a time lag of ∼70 days, with γ−rays leading radio. We interpret that the high-energy and radio fluxes are arising from different emitting regions, located at 1± 12 and 40± 13 pc from the central engine respectively. Subject headings: galaxies: active — galaxies: jets

The TANAMI program has been observing parsec-scale radio jets of southern (declination south of − 30°) γ -ray bright AGN, simultaneously with Fermi /LAT monitoring of their γ -ray emission, via high-resolution radio imaging with Very Long Baseline Interferometry techniques. We present the radio and γ -rayproperties of the TANAMI sources based on one year of contemporaneous TANAMI and Fermi /LAT data. A large fraction (72%) of the TANAMI sample can be associated with bright γ -ray sources for this time range. Association rates differ for different optical classes with all BL Lacs, 76% of quasars, and just 17% of galaxies detected by the LAT. Upper limits were established on the γ -ray flux from TANAMI sources not detected by LAT. This analysis led to the identification of three new Fermi sources whose detection was later confirmed. The γ -ray and radio luminosities are related by L γ ∝ L r 0.89±0.04 . The brightness temperatures of the radio cores increase with the average γ -ray luminosity and the presence of brightness temperatures above the inverse Compton limit implies strong Doppler boosting in those sources. The undetected sources have lower γ /radio luminosity ratios and lower contemporaneous brightness temperatures. Unless the Fermi /LAT-undetected blazars are much γ -ray-fainter than the Fermi /LAT-detected sources, their γ -ray luminosity should not be significantly lower than the upper limits calculated here.