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.

The Large Area Telescope (LAT), the primary instrument for the Fermi Gamma-ray Space Telescope (Fermi) mission, is an imaging, wide field-of-view, high-energy gamma-ray telescope, covering the energy range from 30 MeV to more than 300 GeV. We describe the performance of the instrument at the 10 yr milestone. LAT performance remains well within the specifications defined during the planning phase, validating the design choices and supporting the compelling case to extend the duration of the Fermi mission. The details provided here will be useful when designing the next generation of high-energy gamma-ray observatories.

An excess of $γ$-ray emission from the Galactic Center (GC) region with respect to predictions based on a variety of interstellar emission models and $γ$-ray source catalogs has been found by many groups using data from the {\it Fermi} Large Area Telescope (LAT). Several interpretations of this excess have been invoked. In this paper we search for members of an unresolved population of $γ$-ray pulsars located in the inner Galaxy that are predicted by the interpretation of the GC excess as being due to a population of such sources. We use cataloged LAT sources to derive criteria that efficiently select pulsars with very small contamination from blazars. We search for point sources in the inner $40^\circ\times40^\circ$ region of the Galaxy, derive a list of approximately 400 sources, and apply pulsar selection criteria to extract pulsar candidates among our source list. We performed the entire data analysis chain with two different interstellar emission models (IEMs), and found a total of 135 pulsar candidates, of which 66 were selected with both IEMs.

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.

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

We present a search for spatial extension in high-latitude (vertical bar b vertical bar > 5 degrees) 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 gamma-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 x 10(-16) G for a coherence length lambda greater than or similar to 10 kpc, even when conservative assumptions on the source duty cycle are made. This improves previous limits by several orders of magnitude.

Millisecond pulsars (MSPs) are old neutron stars that spin hundreds of times per second and appear to pulsate as their emission beams cross our line of sight. To date, radio pulsations have been detected from all rotation-powered MSPs. In an attempt to discover radio-quiet gamma-ray MSPs, we used the aggregated power from the computers of tens of thousands of volunteers participating in the Einstein@Home distributed computing project to search for pulsations from unidentified gamma-ray sources in Fermi Large Area Telescope data. This survey discovered two isolated MSPs, one of which is the only known rotation-powered MSP to remain undetected in radio observations. These gamma-ray MSPs were discovered in completely blind searches without prior constraints from other observations, raising hopes for detecting MSPs from a predicted Galactic bulge population.

Black holes with masses below approximately 10(15) g are expected to emit gamma-rays with energies above a few tens of MeV, which can be detected by the Fermi Large Area Telescope (LAT). Although black holes with these masses cannot be formed as a result of stellar evolution, they may have formed in the early universe and are therefore called primordial black holes (PBHs). Previous searches for PBHs have focused on either short-timescale bursts or the contribution of PBHs to the isotropic gamma-ray emission. We show that, in cases of individual PBHs, the Fermi-LAT is most sensitive to PBHs with temperatures above approximately 16 GeV and masses 6 x 10(11) g, which it can detect out to a distance of about 0.03 pc. These PBHs have a remaining lifetime of months to years at the start of the Fermi mission. They would appear as potentially moving point sources with gamma-ray emission that become spectrally harder and brighter with time until the PBH completely evaporates. In this paper, we develop a new algorithm to detect the proper motion of gamma-ray point sources, and apply it to 318 unassociated point sources at a high galactic latitude in the third Fermi-LAT source catalog. None of the unassociated point sources with spectra consistent with PBH evaporation show significant proper motion. Using the nondetection of PBH candidates, we derive a 99% confidence limit on the PBH evaporation rate in the vicinity of Earth, <(rho)over dot>(PBH) < 7.2 x 10(3) pc(-3) yr(-1). This limit is similar to the limits obtained with ground-based gamma-ray observatories.

of interesting and frequently exciting new physics. This has been particularly true for the energy range of the J/$ and T; however, results from the PEP/PETRA machines have also been valuable, e.g., the measurement of the B lifetime. Much effort has been expended in the search for new effects. These include searches for charged and neutral fundamental scalars, new families of quarks, new leptons, lepton-quark substructure, and supersymmetric particles, plus many other topics. In this report I will discuss a few of these experiments in some detail, some in a broader review, and unfortunately some not at all. For those that I have missed I apologize, as these experiments are necessary and rather difficult to carry out, while the fate of unrequited searches is frequently ignominy.

We present the Fermi Large Area Telescope (LAT) observations of the binary neutron star merger event GW170817 and the associated short gamma-ray burst (SGRB) GRB 170817A detected by the Fermi Gamma-ray Burst Monitor. The LAT was entering the South Atlantic Anomaly at the time of the LIGO/Virgo trigger (t(GW)) and therefore cannot place constraints on the existence of high-energy (E > 100 MeV) emission associated with the moment of binary coalescence. We focus instead on constraining high-energy emission on longer timescales. No candidate electromagnetic counterpart was detected by the LAT on timescales of minutes, hours, or days after the LIGO/Virgo detection. The resulting flux upper bound (at 95% C. L.) from the LAT is 4.5. x. 10(-10) erg cm(-2) s(-1) in the 0.1-1 GeV range covering a period from tGW. +. 1153 s to t(GW). +. 2027 s. At the distance of GRB 170817A, this flux upper bound corresponds to a luminosity upper bound of 9.7. x. 10(43) erg s(-1), which is five orders of magnitude less luminous than the only other LAT SGRB with known redshift, GRB 090510. We also discuss the prospects for LAT detection of electromagnetic counterparts to future gravitational-wave events from Advanced LIGO/Virgo in the context of GW170817/GRB 170817A.

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.

The Large Area Telescope on board the Fermi Gamma-ray Space Telescope has collected the largest ever sample of high-energy cosmic-ray electron and positron events since the beginning of its operation. Potential anisotropies in the arrival directions of cosmic-ray electrons or positrons could be a signature of the presence of nearby sources. We use almost seven years of data with energies above 42 GeV processed with the Pass 8 reconstruction. The present data sample can probe dipole anisotropies down to a level of 10^-3. We take into account systematic effects that could mimic true anisotropies at this level. We present a detailed study of the event selection optimization of the cosmic-ray electrons and positrons to be used for anisotropy searches. Since no significant anisotropies have been detected on any angular scale, we present upper limits on the dipole anisotropy. The present constraints are among the strongest to date probing the presence of nearby young and middle-aged sources.

We present a measurement of the cosmic-ray electron + positron spectrum between 7 GeV and 2 TeV performed with almost seven years of data collected with the Fermi Large Area Telescope. We find that the spectrum is well fit by a broken power law with a break energy at about 50 GeV. Above 50 GeV, the spectrum is well described by a single power law with a spectral index of 3.07 +/- 0.02 (stat +/- syst) +/- 0.04 (energy measurement). An exponential cutoff lower than 1.8 TeV is excluded at 95% CL.

The spatial extension of a gamma-ray source is an essential ingredient to determine its spectral properties as well as its potential multi-wavelength counterpart. The capability to spatially resolve gamma-ray sources is greatly improved by the newly delivered Fermi-Large Area Telescope (LAT) Pass 8 event-level analysis which provides a greater acceptance and an improved point spread function, two crucial factors for the detection of extended sources. Here, we present a complete search for extended sources located within 7 degrees from the Galactic plane, using 6 years of LAT data above 10 GeV. We find 46 extended sources and provide their morphological and spectral characteristics. This constitutes the first catalog of hard LAT extended sources, named the Fermi Galactic Extended Source Catalog, which allows a thorough study of the properties of the Galactic plane in the sub-TeV domain.

We present the Fermi Gamma-ray Burst Monitor (GBM) and Large Area Telescope (LAT) observations of the LIGO binary black hole merger (BBH) event GW170104. No candidate electromagnetic counterpart was detected by either GBM or LAT. A detailed analysis of the GBM and LAT data over timescales from seconds to days covering the Laser Interferometer Gravitational-wave Observatory (LIGO) localization region is presented. The resulting flux upper bound from the GBM is (5.2-9.4). x. 10(-7) erg cm(-2) s(-1) in the 10-1000 keV range and from the LAT is (0.2-90). x. 10(-9) erg cm(-2) s(-1) in the 0.1-1 GeV range. We also describe the improvements to our automated pipelines and analysis techniques for searching for and characterizing the potential electromagnetic counterparts for future gravitational-wave events from Advanced LIGO/Virgo.

We present the second catalog of flaring gamma-ray sources (2FAV) detected with the Fermi All-sky Variability Analysis (FAVA), a tool that blindly searches for transients over the entire sky observed by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope. With respect to the first FAVA catalog, this catalog benefits from a larger data set, the latest LAT data release (Pass 8), as well as from an improved analysis that includes likelihood techniques for a more precise localization of the transients. Applying this analysis to the first 7.4 years of Fermi observations, and in two separate energy bands 0.1-0.8 GeV and 0.8-300 GeV, a total of 4547 flares were detected with significance greater than 6s (before trials), on the timescale of one week. Through spatial clustering of these flares, 518 variable gamma-ray sources were identified. Based on positional coincidence, likely counterparts have been found for 441 sources, mostly among the blazar class of active galactic nuclei. For 77 2FAV sources, no likely gamma-ray counterpart has been found. For each source in the catalog, we provide the time, location, and spectrum of each flaring episode. Studying the spectra of the flares, we observe a harder-when-brighter behavior for flares associated with blazars, with the exception of BL Lac flares detected in the low-energy band. The photon indexes of the flares are never significantly smaller than 1.5. For a leptonic model, and under the assumption of isotropy, this limit suggests that the spectrum of freshly accelerated electrons is never harder than p similar to 2.

The region around the Galactic center (GC) is now well established to be brighter at energies of a few GeV than expected from conventional models of diffuse gamma-ray emission and catalogs of known gamma-ray sources. We study the GeV excess using 6.5 years of data from the Fermi Large Area Telescope. We characterize the uncertainty of the GC excess spectrum and morphology due to uncertainties in cosmic-ray source distributions and propagation, uncertainties in the distribution of interstellar gas in the Milky Way, and uncertainties due to a potential contribution from the Fermi bubbles. We also evaluate uncertainties in the excess properties due to resolved point sources of gamma rays. The Galactic center is of particular interest as it would be expected to have the brightest signal from annihilation of weakly interacting massive dark matter particles. However, control regions along the Galactic plane, where a dark-matter signal is not expected, show excesses of similar amplitude relative to the local background. Based on the magnitude of the systematic uncertainties, we conservatively report upper limits for the annihilation cross section as function of particle mass and annihilation channel.