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.

We analyse the implications of the Planck data for cosmic inflation. The Planck nominal mission temperature anisotropy measurements, combined with the WMAP large-angle polarization, constrain the scalar spectral index to be n s = 0.9603 ± 0.0073, ruling out exact scale invariance at over 5σ. Planck establishes an upper bound on the tensor-to-scalar ratio of r < 0.11 (95% CL). The Planck data thus shrink the space of allowed standard inflationary models, preferring potentials with V < 0. Exponential potential models, the simplest hybrid inflationary models, and monomial potential models of degree n ≥ 2 do not provide a good fit to the data. Planck does not find statistically significant running of the scalar spectral index, obtaining dn s /dln k = −0.0134 ± 0.0090. We verify these conclusions through a numerical analysis, which makes no slowroll approximation, and carry out a Bayesian parameter estimation and model-selection analysis for a number of inflationary models including monomial, natural, and hilltop potentials. For each model, we present the Planck constraints on the parameters of the potential and explore several possibilities for the post-inflationary entropy generation epoch, thus obtaining nontrivial data-driven constraints. We also present a direct reconstruction of the observable range of the inflaton potential. Unless a quartic term is allowed in the potential, we find results consistent with second-order slow-roll predictions. We also investigate whether the primordial power spectrum contains any features. We find that models with a parameterized oscillatory feature improve the fit by ∆χ 2 eff ≈ 10; however, Bayesian evidence does not prefer these models. We constrain several single-field inflation models with generalized Lagrangians by combining power spectrum data with Planck bounds on f NL. Planck constrains with unprecedented accuracy the amplitude and possible correlation (with the adiabatic mode) of non-decaying isocurvature fluctuations. The fractional primordial contributions of cold dark matter (CDM) isocurvature modes of the types expected in the curvaton and axion scenarios have upper bounds of 0.25% and 3.9% (95% CL), respectively. In models with arbitrarily correlated CDM or neutrino isocurvature modes, an anticorrelated isocurvature component can improve the χ 2 eff by approximately 4 as a result of slightly lowering the theoretical prediction for the < ∼ 40 multipoles relative to the higher multipoles. Nonetheless, the data are consistent with adiabatic initial conditions.

Recent detections of the starburst galaxies M82 and NGC 253 by gamma-ray telescopes suggest that galaxies rapidly forming massive stars are more luminous at gamma-ray energies compared to their quiescent relatives. Building upon those results, we examine a sample of 69 dwarf, spiral, and luminous and ultraluminous infrared galaxies at photon energies 0.1-100 GeV using 3 years of data collected by the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope (Fermi). Measured fluxes from significantly detected sources and flux upper limits for the remaining galaxies are used to explore the physics of cosmic rays in galaxies. We find further evidence for quasi-linear scaling relations between gamma-ray luminosity and both radio continuum luminosity and total infrared luminosity which apply both to quiescent galaxies of the Local Group and lowredshift starburst galaxies (conservative P-values 0.05 accounting for statistical and systematic uncertainties). The normalizations of these scaling relations correspond to luminosity ratios of log(L 0.1−100 GeV /L 1.4 GHz ) = 1.7 ± 0.1 (statistical) ± 0.2 (dispersion) and log(L 0.1−100 GeV /L 8−1000 µm ) = −4.3 ± 0.1 (statistical) ± 0.2 (dispersion) for a galaxy with a star formation rate of 1 M ⊙ yr −1 , assuming a Chabrier initial mass function. Using the relationship between infrared luminosity and gamma-ray luminosity, the collective intensity of unresolved star-forming galaxies at redshifts 0 < z < 2.5 above 0.1 GeV is estimated to be 0.4-2.4 ×10 −6 ph cm −2 s −1 sr −1 (4-23% of the intensity of the isotropic diffuse component measured with the LAT). We anticipate that ∼ 10 galaxies could be detected by their cosmic-ray induced gamma-ray emission during a 10-year Fermi mission.