Extreme Blazars Research Articles

The variety of extreme blazars in the AstroSat view

Context. Among the blazar class, extreme blazars have exceptionally hard intrinsic X-ray/TeV spectra, and extreme peak energies in their spectral energy distribution (SED). Observational evidence suggests that the non-thermal emission from extreme blazars is typically non-variable. All these unique features present a challenging case for blazar emission models, especially regarding those sources with hard TeV spectra. Aims. We aim to explore the X-ray and GeV observational features of a variety of extreme blazars, including extreme-TeV, extreme-synchrotron (extreme-Syn), and regular high-frequency-peaked BL Lac objects (HBLs). Furthermore, we aim to test the applicability of various blazar emission models that could explain the very hard TeV spectra. Methods. We conducted a detailed spectral analysis of X-ray data collected with AstroSat and Swift-XRT, along with quasi-simultaneous γ-ray data from Fermi-LAT, for five sources: 1ES 0120+340, RGB J0710+591, 1ES 1101−232, 1ES 1741+196, and 1ES 2322−409. We took three approaches to modelling the SEDs: (1) a steady-state one-zone synchrotron-self-Compton (SSC) code, (2) another leptonic scenario of co-accelerated electrons and protons on multiple shocks applied to the extreme-TeV sources only (e–p co-acceleration scenario), and (3) a one-zone hadro-leptonic (ONEHALE) code. The latter code is used twice to explain the γ-ray emission process: proton synchrotron and synchrotron emission of secondary pairs. Results. Our X-ray analysis provides well-constrained estimates of the synchrotron peak energies for both 1ES0120+340 and 1ES1741+196. These findings categorise these latter objects as extreme-synchrotron sources, as they consistently exhibit peak energies above 1 keV in different flux states. The multi-epoch X-ray and GeV data reveal spectral and flux variabilities in RGB J0710+591 and 1ES 1741+196, even on timescales of days to weeks. As anticipated, the one-zone SSC model adequately reproduces the SEDs of regular HBLs but encounters difficulties in explaining the hardest TeV emission. Hadronic models offer a reasonable fit to the hard TeV spectrum, though with the trade-off of requiring extreme jet powers. On the other hand, the lepto-hadronic scenario faces additional challenges in fitting the GeV spectra of extreme-TeV sources. Finally, the e–p co-acceleration scenario naturally accounts for the observed hard electron distributions and effectively matches the hardest TeV spectrum of RGB J0710+591 and 1ES 1101−232.

Extreme blazars: the result of unstable recollimated jets?

ABSTRACT Extreme BL Lacs (EHBL) form a subclass of blazars, which challenge standard emission scenarios. In a recent study it has been argued that their peculiar properties can be explained if emitting electrons are accelerated in a series of oblique shocks induced by the recollimation of the relativistic jet. However, new 3 D simulations of recollimated, weakly magnetized jets reveal that, in correspondence with the first recollimation shock, the flow develops a rapidly growing instability, becomes highly turbulent and decelerates, effectively hampering the formation of the multiple shock structure routinely observed in 2 D simulations. Building on these new findings, we propose here a revised scenario for EHBL, in which the emission is produced by electrons accelerated at the recollimation shock and subsequently further energized through stochastic acceleration in the turbulent downstream flow. We apply a simple version of this scenario to the prototypical EHBL 1ES 0229-200, showing that the spectral energy distribution (SED) can be satisfactorily reproduced with standard values of the main physical parameters.

Absorption Features in Sub-TeV Gamma-Ray Spectra of BL Lac Objects

Abstract The production site of gamma rays in blazars is closely related to their interaction with the photon fields surrounding the active galactic nucleus. In this paper, we discuss an indirect method that may help to unveil the presence of ambient structures in BL Lac objects through the analysis of their gamma-ray spectrum. Gamma rays, passing through structures at different distances from the black hole, interact via γ γ pair production with the corresponding photon fields and produce absorption features in their spectral energy distribution. An interaction with a putative broad-line region may reduce the gamma-ray flux only if its production site were very close to the central engine. However, if jet photons interact with a bath of optical-UV seed photons produced by a narrow-line region extended over the parsec scale, the consequent γ γ process may cause absorption features detectable at a few hundreds of GeV. The detection of such absorption features is facilitated in sources with spectra reaching TeV energies, and specifically HBLs and EHBLs (extreme blazars) may represent exceptional probes to investigate this topic. We discuss recent observations of an extreme blazar named 2WHSP J073326.7+515354 (or PGC 2402248), which shows evidence of such an absorption feature in its gamma-ray spectrum and narrow emission lines in the optical spectrum, suggesting the presence of narrow-line regions in its large-scale environment. Finally, we discuss how sub-TeV absorption features in the spectra of BL Lac objects may affect their broadband modeling and eventually represent a powerful diagnostic tool to constrain the gamma-ray production site and the jet environment.

Electron-proton co-acceleration on relativistic shocks in extreme-TeV blazars

Aims. The multi-wavelength emission from a newly identified population of ‘extreme-TeV’ blazars, with Compton peak frequencies around 1 TeV, is difficult to interpret with standard one-zone emission models. Large values of the minimum electron Lorentz factor and quite low magnetisation values seem to be required. Methods. We propose a scenario where protons and electrons are co-accelerated on internal or recollimation shocks inside the relativistic jet. In this situation, energy is transferred from the protons to the electrons in the shock transition layer, leading naturally to a high minimum Lorentz factor for the latter. A low magnetisation favours the acceleration of particles in relativistic shocks. Results. The shock co-acceleration scenario provides additional constraints on the set of parameters of a standard one-zone lepto-hadronic emission model, reducing its degeneracy. Values of the magnetic field strength of a few mG and minimum electron Lorentz factors of 103 to 104, required to provide a satisfactory description of the observed spectral energy distributions of extreme blazars, result here from first principles. While acceleration on a single standing shock is sufficient to reproduce the emission of most of the extreme-TeV sources we have examined, re-acceleration on a second shock appears needed for those objects with the hardest γ-ray spectra. Emission from the accelerated proton population, with the same number density as the electrons but in a lower range of Lorentz factors, is strongly suppressed. Satisfactory self-consistent representations were found for the most prominent representatives of this new blazar class.

Modeling Gamma-Ray SEDs and Angular Extensions of Extreme TeV Blazars from Intergalactic Proton-Initiated Cascades in Contemporary Astrophysical EGMF Models

The intrinsic spectra of some distant blazars known as “extreme TeV blazars” have shown a hint at an anomalous hardening in the TeV energy region. Several extragalactic propagation models have been proposed to explain this possible excess transparency of the Universe to gamma-rays starting from a model which assumes the existence of so-called axion-like particles (ALPs) and the new process of gamma-ALP oscillations. Alternative models suppose that some of the observable gamma-rays are produced in the intergalactic cascades. This work focuses on investigating the spectral and angular features of one of the cascade models, the Intergalactic Hadronic Cascade Model (IHCM) in the contemporary astrophysical models of Extragalactic Magnetic Field (EGMF). For IHCM, EGMF largely determines the deflection of primary cosmic rays and electrons of intergalactic cascades and, thus, is of vital importance. Contemporary Hackstein models are considered in this paper and compared to the model of Dolag. The models assumed are based on simulations of the local part of large-scale structure of the Universe and differ in the assumptions for the seed field. This work provides spectral energy distributions (SEDs) and angular extensions of two extreme TeV blazars, 1ES 0229+200 and 1ES 0414+009. It is demonstrated that observable SEDs inside a typical point spread function of imaging atmospheric Cherenkov telescopes (IACTs) for IHCM would exhibit a characteristic high-energy attenuation compared to the ones obtained in hadronic models that do not consider EGMF, which makes it possible to distinguish among these models. At the same time, the spectra for IHCM models would have longer high energy tails than some available spectra for the ALP models and the universal spectra for the Electromagnetic Cascade Model (ECM). The analysis of the IHCM observable angular extensions shows that the sources would likely be identified by most IACTs not as point sources but rather as extended ones. These spectra could later be compared with future observation data of such instruments as Cherenkov Telescope Array (CTA) and LHAASO.

Optical spectral characterization of the gamma-ray blazars S4 0954+65, TXS 1515−273, and RX J0812.0+0237

ABSTRACT The study of gamma-ray blazars is usually hindered due to the lack of information on their redshifts and on their low-energy photon fields. This information is key to understand the effect on the gamma-ray absorption due to either extragalactic background light and/or intrinsic absorption and emission processes. All this information has also an impact on the determination of the location of the emitting region within the relativistic jets. In this work, a new optical spectroscopic characterization is presented for three gamma-ray blazars: S4 0954+65, TXS 1515−273, and RX J0812.0+0237. For all the three targets, the redshift determination is successful, and for the first time in the case of TXS 1515−273 and RX J0812.0+0237. Their classification as BL Lac type is confirmed based on these new optical spectra. For S4 0954+65 (z = 0.3694 ± 0.0011), an estimation on the disc, broad-line region, and torus luminosities is performed based on the observed optical emission lines. The results from this study are compatible with the nature of S4 0954+65 as a transitional blazar. In the case of TXS 1515−273 (z = 0.1281 ± 0.0004), although its optical spectrum is dominated by the continuum emission from the jet, applying the pPXF technique, the stellar population can be unveiled and is compatible with an old and metallic population. It is also the case of RX J0812.0+0237 (z = 0.1721 ± 0.0002). Moreover, this work confirms that the optical spectrum from RX J0812.0+0237 is compatible with an extreme blazar classification.

The strange case of the transient HBL blazar 4FGL J1544.3−0649

ABSTRACT We present a multifrequency study of the transient γ-ray source 4FGL J1544.3−0649, a blazar that exhibited a remarkable behaviour raising from the state of an anonymous mid-intensity radio source, never detected at high energies, to that of one of the brightest extreme blazars in the X-ray and γ-ray sky. Our analysis shows that the averaged γ-ray spectrum is well described by a power law with a photon index of 1.87 ± 0.04, while the flux above 100 MeV is (8.0 ± 0.9) × 10−9 photon cm−2 s−1, which increases during the active state of the source. The X-ray flux and spectral slope are both highly variable, with the highest 2–10 keV flux reaching (1.28 ± 0.05) × 10−10 erg cm−2 s−1. On several observations, the X-ray spectrum hardened to the point implying as SED peak moving to energies larger than 10 keV. As in many extreme blazars the broad-band spectral energy distribution can be described by a homogeneous one-zone synchrotron-self-Compton leptonic model. We briefly discuss the potential implications for high-energy multimessenger astrophysics in case the dual behaviour shown by 4FGL J1544.3−0649 does not represent an isolated case, but rather a manifestation of a so far unnoticed relatively common phenomenon.

Hadronic X-Ray Flares from Blazars

Abstract The detection of a high-energy neutrino from the flaring blazar TXS 0506+056 and the subsequent discovery of a neutrino excess from the same direction have strengthened the hypothesis that blazars are cosmic neutrino sources. The lack, however, of γ-ray flaring activity during the latter period challenges the standard scenario of correlated γ-ray and high-energy neutrino emission in blazars. We propose instead that TeV–PeV neutrinos are produced in coincidence with X-ray flares that are powered by proton synchrotron radiation. In this case, neutrinos are produced by photomeson interactions of protons with their own synchrotron radiation, while MeV to GeV γ-rays are the result of synchrotron-dominated electromagnetic cascades developed in the source. Using a time-dependent approach, we find that this “pure hadronic flaring” hypothesis has several interesting consequences. The X-ray flux is a good proxy for the all-flavor neutrino flux, while certain neutrino-rich X-ray flares may be dark in GeV–TeV γ-rays. Lastly, hadronic X-ray flares are accompanied by an equally bright MeV component that is detectable by proposed missions like e-ASTROGAM and AMEGO. We applied this scenario to the extreme blazar 3HSP J095507.9+355101 which has been associated with IceCube-200107A while undergoing an X-ray flare. We showed that the number of muon and anti-muon neutrinos above 100 TeV during hadronic flares can be up to ∼3–10 times higher than the expected number in standard leptohadronic models. Still, frequent hadronic flaring activity is necessary for explaining the detected neutrino event IceCube-200107A.

Comprehensive Multimessenger Modeling of the Extreme Blazar 3HSP J095507.9+355101 and Predictions for IceCube

Abstract 3HSP J095507.9+355101 is an extreme blazar that has been possibly associated with a high-energy neutrino (IceCube-200107A) detected 1 day before the blazar was found to undergo a hard X-ray flare. We perform a comprehensive study of the predicted multimessenger emission from 3HSP J095507.9+355101 during its recent X-ray flare, but also in the long term. We focus on one-zone leptohadronic models, but we also explore alternative scenarios: (i) a blazar-core model, which considers neutrino production in the inner jet, close to the supermassive black hole; (ii) a hidden external-photon model, which considers neutrino production in the jet through interactions with photons from a weak broad line region; (iii) a proton-synchrotron model, where high-energy protons in the jet produce γ-rays via synchrotron; and (iv) an intergalactic cascade scenario, where neutrinos are produced in the intergalactic medium by interactions of a high-energy cosmic-ray beam escaping the jet. The Poisson probability to detect a single muon neutrino in 10 years from 3HSP J095507.9+355101 with the real-time IceCube alert analysis is ∼1% (3%) for the most optimistic one-zone leptohadronic model (the multi-zone blazar-core model). Meanwhile, detection of a single neutrino during the 44-day-long high X-ray flux-state period following the neutrino detection is 0.06%, according to our most optimistic leptohadronic model. The most promising scenarios for neutrino production also predict strong intrasource γ-ray attenuation above ∼100 GeV. If the association is real, then IceCube-Gen2 and other future detectors should be able to provide additional evidence for neutrino production in 3HSP J095507.9+355101 and other extreme blazars.

3HSP J095507.9+355101: A flaring extreme blazar coincident in space and time with IceCube-200107A

The uncertainty region of the highly energetic neutrino IceCube200107A includes 3HSP J095507.9+355101 (z = 0.557), an extreme blazar, which was detected in a high, very hard, and variable X-ray state shortly after the neutrino arrival. Following a detailed multiwavelength investigation, we confirm that the source is a genuine BL Lac. This new detection differs from TXS 0506+056, which is thus far the first source associated with IceCube neutrinos, and is considered a “masquerading” BL Lac. As in the case of TXS 0506+056, 3HSP J095507.9+355101 is also way off the so-called blazar sequence. We consider 3HSP J095507.9+355101 a possible counterpart to the IceCube neutrino. Finally, we discuss some theoretical implications in terms of neutrino production.

Optical spectral characterization of the TeV extreme blazar 2WHSP J073326.7+515354

ABSTRACT The emission from the relativistic jets in blazars usually outshines their host galaxies, challenging the determination of their distances and the characterization of the stellar population. The situation becomes more favourable in the case of the extreme blazars (EHBLs), for which the bulk of the emission of the relativistic jets is emitted at higher energies, unveiling the optical emission from the host galaxy. The distance determination is fundamental for the study of the intrinsic characteristics of the blazars, especially to estimate the intrinsic gamma-ray spectra distorted due to the interaction with the extragalactic background light. In this work, we report on the properties of 2WHSP J073326.7+515354 host galaxy in the optical band, which is one of the few EHBLs detected at TeV energies. We present the first measurement of the distance of the source, z = 0.065 04 ± 0.000 02 (velocity dispersion $\sigma =237 \pm 9\, \mathrm{km s^{-1}}$). We also perform a detailed study of the stellar population of its host galaxy. We find that the mass-weighted mean stellar age is $11.72\pm 0.06\, \mathrm{Gyr}$ and the mean metallicity [M/H] = 0.159 ± 0.016. In addition, a morphological study of the host galaxy is also carried out. The surface brightness distribution is modelled by a composition of a dominant classical bulge (Re = 3.77 ± 1 arcsec or equivalently 4.74 kpc) plus an unresolved source which corresponds to the active nucleus. The black hole mass is estimated using both the mass relation with the velocity dispersion and the absolute magnitude from the bulge yielding comparable results: $(4.8\pm 0.9)\times 10^8$ and $(3.7\pm 1.0)\times 10^8\, \mathrm{ M}_{\odot }$, respectively.

Progress in unveiling extreme particle acceleration in persistent astrophysical jets

The most powerful persistent accelerators in the Universe are jetted active galaxies. Blazars, galaxies whose jets are directed towards Earth, dominate the extragalactic gamma-ray sky. Still, most of the highest-energy particle accelerators likely elude detection. These extreme blazars, whose radiated energy can peak beyond 10 TeV, are ideal targets to study particle acceleration and radiative processes, and may provide links to cosmic rays and astrophysical neutrinos. The growing number of extreme blazars observed at TeV energies has been critical for the emergence of gamma-ray cosmology, including measurements of the extragalactic background light, tight bounds on the intergalactic magnetic field, and constraints on exotic physics at energies inaccessible with human-made accelerators. Tremendous progress has been achieved over the past decade, which bodes well for the future, particularly with the deployment of the Cherenkov Telescope Array.

Testing emission models on the extreme blazar 2WHSP J073326.7+515354 detected at very high energies with the MAGIC telescopes

ABSTRACT Extreme high-energy-peaked BL Lac objects (EHBLs) are an emerging class of blazars. Their typical two-hump-structured spectral energy distribution (SED) peaks at higher energies with respect to conventional blazars. Multiwavelength (MWL) observations constrain their synchrotron peak in the medium to hard X-ray band. Their gamma-ray SED peaks above the GeV band, and in some objects it extends up to several TeV. Up to now, only a few EHBLs have been detected in the TeV gamma-ray range. In this paper, we report the detection of the EHBL 2WHSP J073326.7+515354, observed and detected during 2018 in TeV gamma rays with the MAGIC telescopes. The broad-band SED is studied within an MWL context, including an analysis of the Fermi-LAT data over 10 yr of observation and with simultaneous Swift-XRT, Swift-UVOT, and KVA data. Our analysis results in a set of spectral parameters that confirms the classification of the source as an EHBL. In order to investigate the physical nature of this extreme emission, different theoretical frameworks were tested to model the broad-band SED. The hard TeV spectrum of 2WHSP J073326.7+515354 sets the SED far from the energy equipartition regime in the standard one-zone leptonic scenario of blazar emission. Conversely, more complex models of the jet, represented by either a two-zone spine-layer model or a hadronic emission model, better represent the broad-band SED.

Fermi-LAT Stacking Analysis Technique: An Application to Extreme Blazars and Prospects for their CTA Detection

Abstract We present a likelihood profile stacking technique based on the Fermi-Large Area Telescope (LAT) data to explore the γ-ray characteristics of Fermi-LAT undetected astrophysical populations. The pipeline is applied to a sample of γ-ray unresolved extreme blazars, i.e., sources with the highest synchrotron peak frequencies ( ), and we report a cumulative γ-ray detection with more than 32σ confidence for 2 degrees of freedom. Comparing the generated stacked γ-ray spectrum with the sensitivity limits of the upcoming Cherenkov Telescope Array (CTA), we find that the Fermi-LAT undetected population of such extreme blazars, on average, may remain well below the CTA detection threshold due to their faintness and extragalactic background light (EBL) absorption. However, γ-ray detected blazars belonging to the same class are promising candidates for CTA observations. The EBL-corrected stacked spectra of these sources do not show any softening up to 1 TeV. This finding suggests the inverse Compton peak of extreme blazars lies above 1 TeV, thus indicating a hard intrinsic TeV spectrum. Our analysis also predicts that at 100 GeV, at least ∼10% of the diffuse extragalactic γ-ray background originates from the γ-ray undetected extreme blazars. These results highlight the effectiveness of the developed stacking technique to explore the uncharted territory of γ-ray undetected astrophysical objects.

A multiwavelength view of BL Lacs neutrino candidates

The discovery of high-energy astrophysical neutrinos by IceCube kicked off a new line of research to identify the electromagnetic counterparts producing these neutrinos. Among the extragalactic sources, active galactic nuclei (AGN), and in particular Blazars, are promising candidate neutrino emitters. Their structure, with a relativistic jet pointing to the Earth, offers a natural accelerator of particles and for this reason a perfect birthplace of high energy neutrinos. A good characterisation of the spectral energy distribution (SED) of these sources can improve the understanding of the physical composition of the source and the emission processes involved. Starting from our previous works in which we assumed a correlation between the $\gamma$-ray and the neutrino flux of the BL Lacs of the 2FHL catalogue (detected by Fermi above 50GeV), we select those BL Lac in spatial correlation with the IceCube events. We obtain a sample of 7 sources and we start an observational campaign to have a better characterisation of the synchrotron peak. During the analysis of the data a new source has been added because of its position inside the angular uncertainty of a muon track event detected by IceCube. This source, namely TXS0506+056, was in a high-state during the neutrino event and we will consider it as benchmark to check the proprieties of the other sources of the sample during the related neutrino detection. We obtain a better characterisation of the SED for the sources of our sample. A prospective extreme Blazar, a very peculiar low synchrotron peak (LSP) source with a large separation of the two peaks and a \textit{twin} of TXS0506+056 come up. We also provide the $\gamma$-ray light curve to check the trend of the sources around the neutrino detection but no clears patterns are in common among the sources.

HESS J1943+213: An Extreme Blazar Shining through the Galactic Plane

Abstract HESS J1943+213 is a very high energy (VHE; >100 GeV) γ-ray source in the direction of the Galactic plane. Studies exploring the classification of the source are converging toward its identification as an extreme synchrotron BL Lac object. Here we present 38 hr of VERITAS observations of HESS J1943+213 taken over 2 yr. The source is detected with a significance of ∼20 standard deviations, showing a remarkably stable flux and spectrum in VHE γ-rays. Multifrequency Very Long Baseline Array (VLBA) observations of the source confirm the extended, jet-like structure previously found in the 1.6 GHz band with the European VLBI Network and detect this component in the 4.6 and 7.3 GHz bands. The radio spectral indices of the core and the jet and the level of polarization derived from the VLBA observations are in a range typical for blazars. Data from VERITAS, Fermi-LAT, Swift-XRT, the FLWO 48″ telescope, and archival infrared and hard X-ray observations are used to construct and model the spectral energy distribution (SED) of the source with a synchrotron self-Compton model. The well-measured γ-ray peak of the SED with VERITAS and Fermi-LAT provides constraining upper limits on the source redshift. Possible contribution of secondary γ-rays from ultra-high-energy cosmic-ray-initiated electromagnetic cascades to the γ-ray emission is explored, finding that only a segment of the VHE spectrum can be accommodated with this process. A variability search is performed across X-ray and γ-ray bands. No statistically significant flux or spectral variability is detected.

The extreme blazar AO 0235+164 as seen by extensive ground and space radio observations

Clues to the physical conditions in radio cores of blazars come from measurements of brightness temperatures as well as effects produced by intrinsic opacity. We study the properties of the ultra compact blazar AO 0235+164 with RadioAstron ground-space radio interferometer, multi-frequency VLBA, EVN and single-dish radio observations. We employ visibility modeling and image stacking for deriving structure and kinematics of the source, and use Gaussian process regression to find the relative multi-band time delays of the flares. The multi-frequency core size and time lags support prevailing synchrotron self absorption. The intrinsic brightness temperature of the core derived from ground-based VLBI is close to the equipartition regime value. In the same time, there is evidence for ultra-compact features of the size of less than 10 $\mu$as in the source, which might be responsible for the extreme apparent brightness temperatures of up to $10^{14}$ K as measured by RadioAstron. In 2007--2016 the VLBI components in the source at 43 GHz are found predominantly in two directions, suggesting a bend of the outflow from southern to northern direction. The apparent opening angle of the jet seen in the stacked image at 43 GHz is two times wider than that at 15 GHz, indicating a collimation of the flow within the central 1.5 mas. We estimate the Lorentz factor $\Gamma = 14$, the Doppler factor $\delta=21$, and the viewing angle $\theta = 1.7^\circ$ of the apparent jet base, derive the gradients of magnetic field strength and electron density in the outflow, and the distance between jet apex and the core at each frequency.

Testing models of extragalactic γ-ray propagation using observations of extreme blazars in GeV and TeV energy ranges

We briefly review contemporary extragalactic γ-ray propagation models. It is shown that the Extragalactic Magnetic Field (EGMF) strength and structure are poorly known. Strict lower limits on the EGMF strength in voids are of order 10−17 − 10−20 G, thus allowing a substantial contribution of a secondary component generated by electromagnetic cascades to the observable spectrum. We show that this “electromagnetic cascade model” is supported by data from imaging Cherenkov telescopes and the Fermi LAT detector.

LOWER BOUND ON THE COSMIC TeV GAMMA-RAY BACKGROUND RADIATION

ABSTRACT The Fermi gamma-ray space telescope has revolutionized our understanding of the cosmic gamma-ray background radiation in the GeV band. However, investigation on the cosmic TeV gamma-ray background radiation still remains sparse. Here, we report the lower bound on the cosmic TeV gamma-ray background spectrum placed by the cumulative flux of individual detected extragalactic TeV sources including blazars, radio galaxies, and starburst galaxies. The current limit on the cosmic TeV gamma-ray background above 0.1 TeV is obtained as 2.8 × 10−8(E/100 GeV)−0.55 exp(−E/2100GeV)[GeV cm−2 s−1 sr−1] < E 2 dN/dE < 1.1 × 10−7(E/100 GeV)−0.49 [GeV cm−2 s−1 sr−1], where the upper bound is set by requirement that the cascade flux from the cosmic TeV gamma-ray background radiation can not exceed the measured cosmic GeV gamma-ray background spectrum. Two nearby blazars, Mrk 421 and Mrk 501, explain ∼70% of the cumulative background flux at 0.8–4 TeV, while extreme blazars start to dominate at higher energies. We also provide the cumulative background flux from each population, i.e., blazars, radio galaxies, and starburst galaxies which will be the minimum requirement for their contribution to the cosmic TeV gamma-ray background radiation.

Extreme blazars as counterparts of IceCube astrophysical neutrinos

We explore the correlation of $\gamma$-ray emitting blazars with IceCube neutrinos by using three very recently completed, and independently built, catalogues and the latest neutrino lists. We introduce a new observable, namely the number of neutrino events with at least one $\gamma$-ray counterpart, $N_{\nu}$. In all three catalogues we consistently observe a positive fluctuation of $N_{\nu}$ with respect to the mean random expectation at a significance level of $0.4 - 1.3$ per cent. This applies only to extreme blazars, namely strong, very high energy $\gamma$-ray sources of the high energy peaked type, and implies a model-independent fraction of the current IceCube signal $\sim 10 - 20$ per cent. An investigation of the hybrid photon -- neutrino spectral energy distributions of the most likely candidates reveals a set of $\approx 5$ such sources, which could be linked to the corresponding IceCube neutrinos. Other types of blazars, when testable, give null correlation results. Although we could not perform a similar correlation study for Galactic sources, we have also identified two (further) strong Galactic $\gamma$-ray sources as most probable counterparts of IceCube neutrinos through their hybrid spectral energy distributions. We have reasons to believe that our blazar results are not constrained by the $\gamma$-ray samples but by the neutrino statistics, which means that the detection of more astrophysical neutrinos could turn this first hint into a discovery.