On this page
The 2018 announcement of an IceCube neutrino event correlated with a gamma-ray flare from the blazar TXS 0506+056 was framed as the opening of multi-messenger neutrino astronomy. The combination of a single 290 TeV neutrino from a specific sky direction, simultaneous enhanced gamma-ray activity from a known blazar at the same direction, and a follow-up search revealing a 3-sigma excess from the same source in IceCube’s historical data was unprecedented. For the first time, a high-energy neutrino had been associated with an identified astrophysical source.
Five years and several follow-up papers later, the picture has filled in. Multi-wavelength observations have characterised the blazar jet structure in detail. Reanalyses of the historical excess have firmed up the significance. IceCube’s continuing operations have provided additional context — including the discovery of the steady-state NGC 1068 source — that has reshaped the way TXS 0506+056 fits into the broader picture. And stacking analyses of the broader blazar population have constrained how representative TXS 0506+056 is of blazars in general.
This post is about what we know now that we didn’t in 2018: how the case has strengthened, what new constraints have emerged, and what the current consensus is on blazars as cosmic neutrino sources.
The 2018 result in retrospect
The original IceCube paper from 2018 (Aartsen et al., Science 361, 147) reported the alert event IceCube-170922A: a 290 TeV muon-neutrino arriving on 22 September 2017 from celestial coordinates within 0.5 degrees of TXS 0506+056. The blazar, a flat-spectrum BL Lac object at redshift z ≈ 0.34, was in an enhanced gamma-ray state at the time of the event, with Fermi-LAT and MAGIC reporting fluxes 3-5 times above the long-term average.
A follow-up analysis of IceCube’s historical data from 2014-2015 (Aartsen et al., Science 361, 147) found a 3-sigma excess of neutrino events from the TXS 0506+056 direction during a 110-day period — an apparent “neutrino flare” that preceded the 2017 detection by years and that did not have a corresponding well-characterised electromagnetic flare.
The combined statistical case was that the chance probability of seeing the September 2017 alert event coincident with an enhanced gamma-ray state from a known blazar was roughly 0.001-0.01 per cent, and the 2014-2015 historical excess was independent corroboration. The total significance of the TXS 0506+056 association came in at approximately 3-3.5 sigma in different ways of counting.
What the multi-wavelength observations revealed
The intensive multi-messenger follow-up of TXS 0506+056 across radio, optical, X-ray, and gamma-ray bands produced detailed characterization of the source.
Radio interferometry at high resolution (VLBI) showed TXS 0506+056’s relativistic jet pointed within a few degrees of Earth’s line of sight, with measured apparent superluminal motion typical of beamed blazar jets. The jet’s structural parameters — Lorentz factor of about 10, doppler factor of about 30 — fit the standard blazar picture.
Optical and X-ray observations during the 2017 outburst identified TXS 0506+056 as having a hybrid emission profile, neither a pure low-frequency-peaked BL Lac nor a high-frequency-peaked one but somewhere between. The X-ray emission showed signs of a two-component structure that hinted at separate hadronic and leptonic emission regions within the jet.
Gamma-ray observations by Fermi-LAT, MAGIC, and HESS during the 2017 outburst showed enhanced flux from 100 MeV to 1 TeV, consistent with synchrotron-self-Compton emission from the leptonic outer zone plus a hadronic-secondary contribution at the highest energies.
Together, the multi-wavelength picture is consistent with a two-zone hadronic-leptonic model: an inner emission region where high-energy protons are accelerated and produce neutrinos through photoproduction interactions with the local radiation field, plus an outer emission region where the electromagnetic emission is dominated by leptonic synchrotron-self-Compton. The neutrino signal therefore comes from a more compact and optically thicker region than the bulk of the gamma-ray emission.
The historical excess reanalysed
The 2014-2015 historical excess has been reanalysed with improved IceCube event reconstruction and updated point-source likelihood methods. The reanalyses have raised the significance from the original 3 sigma to approximately 4 sigma, partly through better angular resolution for the relevant events and partly through more sophisticated treatment of the background expectation.
The 2014-2015 excess remains puzzling: it was a neutrino flare with no clearly corresponding electromagnetic flare, and the source did not show the kind of multi-wavelength signature that the 2017 event was associated with. Two interpretations are commonly considered. The first is that the source’s emission region is optically thick to its own gamma rays at the relevant energies, so the gamma-ray flux is absorbed locally while the neutrinos escape. The second is that the 2014-2015 emission was from a different, structurally distinct component of the source — perhaps a longer-term hadronic accretion episode rather than a flare-like event.
Neither interpretation has been definitively confirmed. The 2014-2015 excess remains the strongest evidence for sustained neutrino emission from TXS 0506+056 over a time interval longer than a flare timescale.
NGC 1068 as context
The most important development since 2018 has been IceCube’s 2022 announcement of NGC 1068 as a steady-state neutrino source. NGC 1068 is a Seyfert galaxy at distance 14 Mpc, much closer than TXS 0506+056, and its inferred neutrino emission is continuous rather than flare-like.
The NGC 1068 detection reshapes the TXS 0506+056 interpretation in important ways. It shows that AGN can produce detectable neutrino emission without a relativistic jet pointing at Earth — NGC 1068’s jet is not beamed toward us. The neutrino emission appears to come from the dense corona around the supermassive black hole rather than from the jet.
This suggests that the broader AGN population may produce neutrinos through multiple mechanisms: jet-beamed emission as in TXS 0506+056 plus corona-driven emission as in NGC 1068. The relative contributions of the two mechanisms to the diffuse astrophysical flux is one of the active research questions, with IceCube’s continuing analyses constraining the AGN population fraction in the 5-10 per cent range — substantial but not dominant.
Blazar stacking limits
IceCube’s broader blazar-population stacking analyses have constrained how representative TXS 0506+056 is of blazars in general. The principle is the same as the GRB stacking analyses discussed in our earlier post: by summing neutrino events from the directions of all known blazars (typically several thousand from the Fermi-LAT 4FGL catalogue), one can detect a population-level signal even if individual sources are too weak.
The stacking analyses have found no significant population-level signal from blazars beyond what would be expected from background fluctuations. The resulting upper limits constrain the blazar contribution to the diffuse astrophysical flux at approximately 5-15 per cent, depending on the assumed source population and weighting scheme.
This is inconsistent with TXS 0506+056 being a typical blazar example. If most blazars produced neutrinos at the TXS 0506+056 rate, the stacking signal would be substantial; the absence of such a signal means TXS 0506+056 was either unusually bright in neutrinos or unusually visible due to favourable geometry. The TXS 0506+056 event therefore appears to be a relatively rare case rather than evidence for blazars as a dominant source class.
What this means for the field
The combined picture from TXS 0506+056 plus NGC 1068 plus the blazar stacking limits is that no single source class dominates the diffuse astrophysical flux. Blazars contribute at the few-to-ten-per-cent level. Seyfert AGN contribute another few per cent through corona-driven mechanisms similar to NGC 1068. The remaining 80-90 per cent of the flux comes from sources not yet individually identified.
Candidate explanations for the unidentified majority include starburst galaxies (consistent with the cosmic-ray-induced mechanism Loeb and Waxman proposed in 2006), low-luminosity AGN below current point-source identification thresholds, tidal disruption events, and choked GRB jets that fail to break through their host stellar envelopes. None of these has been definitively confirmed but the constraints are tightening as IceCube and KM3NeT statistics grow.
The TXS 0506+056 case retains a special place as the first identified multi-messenger neutrino source. Even if blazars are not the dominant population, the 2017 association established the methodology for combining electromagnetic and neutrino observations and provided a template for the NGC 1068 result and subsequent point-source identifications.
Why the followup matters
The progression from “interesting hint” to “established source class” in observational astrophysics typically takes a decade or more. TXS 0506+056 in 2018 was an exciting single-event coincidence; TXS 0506+056 in 2024 is a well-characterised source with multi-wavelength geometry, refined statistical significance, and broader context from the discovery of NGC 1068 and the blazar stacking limits.
The case has not weakened. The 2014-2015 historical excess has firmed up. The multi-wavelength understanding of the source’s jet structure has filled in. NGC 1068 has added an independent identified source against which TXS 0506+056 can be compared.
What has changed is the framing. TXS 0506+056 is no longer “the example of cosmic neutrino sources” but “one example among several with subtle differences.” The diffuse astrophysical flux is increasingly understood as a population sum across multiple source classes rather than as the emission of a dominant single type. This is a more nuanced and more astronomically realistic picture than the one most theorists initially constructed around the 2018 announcement.
Summary
The 2018 TXS 0506+056 association of a 290 TeV IceCube neutrino event with an enhanced gamma-ray state of the blazar TXS 0506+056 has been substantially reinforced by five years of follow-up observations and analyses. The 2014-2015 historical excess has been reanalysed at higher precision to approximately 4 sigma significance. Multi-wavelength observations have characterised the source’s two-zone jet structure consistent with hadronic acceleration in an inner region and leptonic emission in an outer region. The 2022 identification of NGC 1068 as a steady-state neutrino source has provided complementary evidence that AGN are real neutrino emitters across multiple morphological classes. IceCube’s broader blazar stacking limits have shown that TXS 0506+056 is not a representative example of the blazar population — most blazars are not detectable as individual neutrino sources at current sensitivity. The combined picture is that the diffuse astrophysical neutrino flux comes from multiple source classes contributing at the few-per-cent level each, with TXS 0506+056 and NGC 1068 as the first two identified specimens of a more diverse source population that will be characterised in detail by KM3NeT and IceCube-Gen2 over the coming decade.