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.

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 discovery of extraterrestrial very-high-energy neutrinos by the IceCube collaboration has launched a quest for the identification of their astrophysical sources. Gamma-ray blazars have been predicted to yield a cumulative neutrino signal exceeding the atmospheric background above energies of 100 TeV, assuming that both the neutrinos and the γ-ray photons 2 are produced by accelerated protons in relativistic jets. Since the background spectrum falls steeply with increasing energy, the individual events with the clearest signature of being of an extraterrestrial origin are those at PeV energies. Inside the large positional-uncertainty fields of the first two PeV neutrinos detected by IceCube, the integrated emission of the blazar population has a sufficiently high electromagnetic flux to explain the detected IceCube events, but fluences of individual objects are too low to make an unambiguous source association.Here, we report that a major outburst of the blazar PKS B1424−418 occurred in temporal and positional coincidence with the third PeV-energy neutrino event (IC 35) detected by IceCube. Based on an analysis of the full sample of γ-ray blazars in the IC 35 field and assuming a photo-hadronic emission model, we show that the long-term average γ-ray emission of blazars as a class is in agreement with both the measured all-sky flux of PeV neutrinos and the spectral slope of the IceCube signal. The outburst of PKS B1424−418 has provided an energy output high enough to explain the observed PeV event, indicative of a direct physical association.The neutrino excess detected by IceCube comprises 37 events with energies between 30 TeV and 2 PeV, rejecting a purely atmospheric origin at a significance of 5.7 standard deviations [1][2][3] .These events show a broad distribution across both hemispheres of the sky consistent with an extragalactic source population. Due to the very steep background of atmospheric neutrinos, events at PeV energies are best suited for attempting to establish associations with individual blazars. In the first two years of observations, IceCube detected two events with about 1 PeV of deposited energy 1, 2 (IC 14, and IC 20; dubbed 'Bert' and 'Ernie') the diffuse neutrino flux due to the integrated emission of AGN in a given large field at a given time, as well as the maximum possible neutrino flux associated with an individual object of the sample.Blazars are radio-loud AGN with jets oriented close to the line of sight. This substantially increases the apparent brightness of these objects owing to the Doppler boosting of the emission from the relativistically moving emission zones. A direct association of a PeV-neutrino with an * A different analysis of IceCube muon neutrinos finds an excess signal also from the northern sky 14 . † http://pulsar.sternwarte.uni-erlangen.de/tanami 4 individual γ-ray blazar would have the important implication that a sizeable fraction of their observed γ-ray emission must be due to hadronic decays, and that blazar jets are also sources of ultra-high-energy cos...