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.

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 .

The Large Area Telescope (LAT) aboard the Fermi spacecraft routinely observes high-energy emission from gamma-ray bursts (GRBs). Here we present the second catalog of LAT-detected GRBs, covering the first 10 yr of operations, from 2008 to 2018 August 4. A total of 186 GRBs are found; of these, 91 show emission in the range 30-100 MeV (17 of which are seen only in this band) and 169 are detected above 100 MeV. Most of these sources were discovered by other instruments (Fermi/GBM, Swift/BAT, AGILE, INTEGRAL) or reported by the Interplanetary Network (IPN); the LAT has independently triggered on four GRBs. This catalog presents the results for all 186 GRBs. We study onset, duration, and temporal properties of each GRB, as well as spectral characteristics in the 100 MeV-100 GeV energy range. Particular attention is given to the photons with the highest energy. Compared with the first LAT GRB catalog, our rate of detection is significantly improved. The results generally confirm the main findings of the first catalog: the LAT primarily detects the brightest GBM bursts, and the high-energy emission shows delayed onset as well as longer duration. However, in this work we find delays exceeding 1 ks and several GRBs with durations over 10 ks. Furthermore, the larger number of LAT detections shows that these GRBs not only cover the high-fluence range of GBM-detected GRBs but also sample lower fluences. In addition, the greater number of detected GRBs with redshift estimates allows us to study their properties in both the observer and rest frames. Comparison of the observational results with theoretical predictions reveals that no model is currently able to explain all results, highlighting the role of LAT observations in driving theoretical models.

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×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 γ-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.

A majority of the γ-ray emission from star-forming galaxies is generated by the interaction of high-energy cosmic rays with the interstellar gas and radiation fields. Starforming galaxies are expected to contribute to both the extragalactic γ-ray background and the IceCube astrophysical neutrino flux. Using roughly 10 years of γ-ray data taken by the Fermi Large Area Telescope, in this study we constrain the γ-ray properties of star-forming galaxies. We report the detection of 11 bona-fide γ-ray emitting galaxies and 2 candidates. Moreover, we show that the cumulative γ-ray emission of belowthreshold galaxies is also significantly detected at ∼5 σ confidence. The γ-ray luminosity of resolved and unresolved galaxies is found to correlate with the total (8-1000 µm) infrared luminosity as previously determined. Above 1 GeV, the spectral energy distribution of resolved and unresolved galaxies is found to be compatible with a power law with a photon index of ≈ 2.2 − 2.3. Finally, we find that star-forming galaxies account for roughly 5 % and 3 % of the extragalactic γ-ray background and the IceCube neutrino flux, respectively.

We present a multi-wavelength study of the radio-loud narrow line Seyfert 1 galaxy (NLSy1), 1H 0323+342, detected by Fermi Gamma Ray Space Telescope. Multi-band light curves show many orphan X-ray and optical flares having no corresponding γ-ray counterparts. Such anomalous variability behavior can be due to different locations of the emission region from the central source. During a large flare, γ-ray flux doubling time scale as small as ∼ 3 hours is noticed. We built spectral energy distribution (SED) during different activity states and modeled them using an one-zone leptonic model. The shape of the optical/UV component of the SEDs is dominated by accretion disk emission in all the activity states. In the X-ray band, significant thermal emission from the hot corona is inferred during quiescent and first flaring states, however, during subsequent flares, non-thermal jet component dominates. The γ-ray emission in all the states can be well explained by inverse-Compton scattering of accretion disk photons reprocessed by the broad line region. The source showed violent intra-night optical variability, coinciding with one of the high γ-ray activity states. An analysis of the overall X-ray spectrum fitted with an absorbed power-law plus relativistic reflection component hints for the presence of Fe K-α line and returns a high black hole spin value of a=0.96 ± 0.14. We argue that 1H 0323+342 possesses dual characteristics, akin to flat spectrum radio quasars (FSRQs) as well as radio-quiet NLSy1s, though at a low jet power regime compared to powerful FSRQs.