On 2017 August 17 a binary neutron star coalescence candidate (later designated GW170817) with merger time 12:41:04 UTC was observed through gravitational waves by the Advanced LIGO and Advanced Virgo detectors. The Fermi Gamma-ray Burst Monitor independently detected a gamma-ray burst (GRB 170817A) with a time delay of ∼ 1.7 s with respect to the merger time. From the gravitational-wave signal, the source was initially localized to a sky region of 31 deg2 at a luminosity distance of 40 − 8 + 8 Mpc and with component masses consistent with neutron stars. The component masses were later measured to be in the range 0.86 to 2.26 M ⊙ . An extensive observing campaign was launched across the electromagnetic spectrum leading to the discovery of a bright optical transient (SSS17a, now with the IAU identification of AT 2017gfo) in NGC 4993 (at ∼ 40 Mpc ) less than 11 hours after the merger by the One-Meter, Two Hemisphere (1M2H) team using the 1 m Swope Telescope. The optical transient was independently detected by multiple teams within an hour. Subsequent observations targeted the object and its environment. Early ultraviolet observations revealed a blue transient that faded within 48 hours. Optical and infrared observations showed a redward evolution over ∼10 days. Following early non-detections, X-ray and radio emission were discovered at the transient’s position ∼ 9 and ∼ 16 days, respectively, after the merger. Both the X-ray and radio emission likely arise from a physical process that is distinct from the one that generates the UV/optical/near-infrared emission. No ultra-high-energy gamma-rays and no neutrino candidates consistent with the source were found in follow-up searches. These observations support the hypothesis that GW170817 was produced by the merger of two neutron stars in NGC 4993 followed by a short gamma-ray burst (GRB 170817A) and a kilonova/macronova powered by the radioactive decay of r-process nuclei synthesized in the ejecta.

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.

Previous detections of individual astrophysical sources of neutrinos are limited to the Sun and the supernova 1987A, whereas the origins of the diffuse flux of high-energy cosmic neutrinos remain unidentified. On 22 September 2017, we detected a high-energy neutrino, IceCube-170922A, with an energy of ~290 tera-electron volts. Its arrival direction was consistent with the location of a known γ-ray blazar, TXS 0506+056, observed to be in a flaring state. An extensive multiwavelength campaign followed, ranging from radio frequencies to γ-rays. These observations characterize the variability and energetics of the blazar and include the detection of TXS 0506+056 in very-high-energy γ-rays. This observation of a neutrino in spatial coincidence with a γ-ray-emitting blazar during an active phase suggests that blazars may be a source of high-energy neutrinos.

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.

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 .