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

No abstract

In this paper, we have selected a sample of 64 teraelectronvolt blazars, with redshift, from those classified in the fourth Fermi Large Area Telescope source catalog a) . We have obtained the values of the relevant physical parameters by performing a log-parabolic fitting of the average-state multiwavelength spectral energy distributions. We estimate the range of the radiation zone parameters, such as the Doppler factor (D), the magnetic field strength (B), the radiative zone radius (R) and the peak Lorentz factor (γ p ) of nonthermal electrons. Here, we show that (1) there is a strong linear positive correlation between the intrinsic synchrotron peak frequency and the intrinsic inverse Compton scattering (ICs) peak frequency among different types of blazars; (2) if radio bands are excluded, the spectral index of each band is negatively correlated with the intrinsic peak frequency; (3) there is a strong linear negative correlation between the curvature at the peak and the intrinsic peak frequency of the synchrotron bump, and a weak positive correlation between the curvature at the peak and the intrinsic peak frequency of the ICs bump; (4) there is a strong linear positive correlation between the intrinsic ICs peak luminosity and intrinsic γ-ray luminosity and between the intrinsic ICs peak frequency and peak Lorentz factor;(5) there is a strong negative linear correlation between log B and log γ p ; and (6) there is no correlation between log R and log γ p .

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.

In this paper, multi-wavelength data are compiled for a sample of 1425 Fermi blazars to calculate their spectral energy distributions (SEDs). A parabolic function, log(νF ν ) = P 1 (logν − P 2 ) 2 + P 3 , is used for SED fitting. Synchrotron peak frequency (logν p ), spectral curvature (P 1 ), peak flux (ν p F νp ), and integrated flux (νF ν ) are successfully obtained for 1392 blazars (461 flat spectrum radio quasarsFSRQs, 620 BL Lacs-BLs and 311 blazars of uncertain type-BCUs, 999 sources have known redshifts). Monochromatic luminosity at radio 1.4 GHz, optical R band, X-ray at 1 keV and γ-ray at 1 GeV, peak luminosity, integrated luminosity and effective spectral indexes of radio to optical (α RO ), and optical to X-ray (α OX ) are calculated. The "Bayesian classification" is employed to logν p in the rest frame for 999 blazars with available redshift and the results show that 3 components are enough to fit the logν p distribution, there is no ultra high peaked subclass. Based on the 3 components, the subclasses of blazars using the acronyms of Abdo et al. (2010a) are classified, and some mutual correlations are also studied. Conclusions are finally drawn as follows: (1) SEDs are successfully obtained for 1392 blazars. The fitted peak frequencies are compared with common sources from samples available ( Sambruna et al. 1996, Nieppola et al. 2006, 2008, Abdo et al. 2010a. (2) -2 -peak sources (ISPs) if 14.0 < log ν p (Hz) ≤ 15.3, and high synchrotron peak sources (HSPs) if log ν p (Hz) > 15.3. (3) γ-ray emissions are strongly correlated with radio emissions. γ-ray luminosity is also correlated with synchrotron peak luminosity and integrated luminosity. (4) There is an anti-correlation between peak frequency and peak luminosity within the whole blazar sample. However, there is a marginally positive correlation for HBLs, and no correlations for FSRQs or LBLs. (5) There are anti-correlations between the monochromatic luminosities (γ-ray and radio bands) and the peak frequency within the whole sample and BL Lacs. (6) The optical to X-ray (α OX ) and radio to optical (α RO ) spectral indexes are strongly anti-correlated with peak frequency (log ν p ) within the whole sample, but the correlations for subclasses of FSRQs, LBLs, and HBLs are different.

We present the first Fermi-Large Area Telescope (LAT) solar flare catalog covering the 24th solar cycle. This catalog contains 45 Fermi-LAT solar flares (FLSFs) with emission in the γ-ray energy band (30 MeV–10 GeV) detected with a significance of ≥5σ over the years 2010–2018. A subsample containing 37 of these flares exhibits 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 observations of GeV emission from three flares originating from active regions located behind the limb of the visible solar disk. We report the lightcurves, spectra, best proton index, and localization (when possible) for all FLSFs. The γ-ray spectra are consistent with the decay of pions produced by >300 MeV protons. This work contains the largest sample of high-energy γ-ray flares ever reported and provides a unique opportunity to perform population studies on the different phases of the flare and thus allowing a new window in solar physics to be opened.

The Crab Nebula is a bright source of gamma-rays powered by the Crab Pulsar's rotational energy, through the formation and termination of a relativistic electron-positron wind. We report the detection of γ-rays from this source with energies from 5 × 10−4 to 1.1 petaelectronvolts (PeV), with a spectrum showing gradual steepening over three energy decades. The ultra-high-energy photons imply the presence of a PeV electron accelerator (a pevatron) in the nebula, with an acceleration rate exceeding 15% of the theoretical limit. We constrain the pevatron's size between 0.025 and 0.1 pc, and magnetic field ≈110 μG. The production rate of PeV electrons, 2.5 × 1036 erg s−1, constitutes 0.5% of the pulsar spin-down luminosity, although we cannot exclude a contribution of PeV protons to the production of the highest energy γ-rays.

A sub-array of the Large High Altitude Air Shower Observatory (LHAASO), KM2A is mainly designed to observe a large fraction of the northern sky to hunt for γ-ray sources at energies above 10 TeV. Even though the detector construction is still underway, half of the KM2A array has been operating stably since the end of 2019. In this paper, we present the KM2A data analysis pipeline and the first observation of the Crab Nebula, a standard candle in very high energy γ-ray astronomy. We detect γ-ray signals from the Crab Nebula in both energy ranges of 10 100 TeV and 100 TeV with high significance, by analyzing the KM2A data of 136 live days between December 2019 and May 2020. With the observations, we test the detector performance, including angular resolution, pointing accuracy and cosmic-ray background rejection power. The energy spectrum of the Crab Nebula in the energy range 10-250 TeV fits well with a single power-law function dN/dE = (1.13 0.05 0.08 ) 10 (E/20 TeV) cm s TeV . It is consistent with previous measurements by other experiments. This opens a new window of γ-ray astronomy above 0.1 PeV through which new ultrahigh-energy γ-ray phenomena, such as cosmic PeVatrons, might be discovered.