Current Cherenkov telescopes have identified a population of ultra-high-frequency peaked BL Lac objects (UHBLs), also known as extreme blazars, that exhibit exceptionally hard TeV spectra, including 1ES 0229+200, 1ES 0347−121, RGB J0710+591, 1ES 1101−232, and 1ES 1218+304. Although one-zone synchrotron-self-Compton (SSC) models have been generally successful in interpreting the high-energy emission observed in other BL Lac objects, they are problematic for UHBLs, necessitating very large Doppler factors and/or extremely high minimum Lorentz factors of the emitting leptonic population. In this context, we have investigated alternative scenarios where hadronic emission processes are important, using a newly developed (lepto-)hadronic numerical code to systematically explore the physical parameters of the emission region that reproduces the observed spectra while avoiding the extreme values encountered in pure SSC models. Assuming a fixed Doppler factor δ = 30, two principal parameter regimes are identified, where the high-energy emission is due to: (1) proton-synchrotron radiation, with magnetic fields B ∼ 1–100 G and maximum proton energies E_{p; max} ≲ 10¹⁹ eV; and (2) synchrotron emission from p–γ-induced cascades as well as SSC emission from primary leptons, with B ∼ 0.1–1 G and E_{p; max} ≲ 10¹⁷ eV. This can be realized with plausible, sub-Eddington values for the total (kinetic plus magnetic) power of the emitting plasma, in contrast to hadronic interpretations for other blazar classes that often warrant highly super-Eddington values.

While active galactic nuclei with relativistic jets have long been prime candidates for the origin of extragalactic cosmic rays and neutrinos, the BL Lac object TXS 0506+056 is the first astrophysical source observed to be associated with some confidence (∼3σ) with a high-energy neutrino, IceCube-170922A, detected by the IceCube Observatory. The source was found to be active in high-energy γ-rays with Fermi-LAT and in very-high-energy γ-rays with the MAGIC telescopes. To consistently explain the observed neutrino and multi-wavelength electromagnetic emission of TXS 0506+056, we investigate in detail single-zone models of lepto-hadronic emission, assuming co-spatial acceleration of electrons and protons in the jet, and synchrotron photons from the electrons as targets for photo-hadronic neutrino production. The parameter space concerning the physical conditions of the emission region and particle populations is extensively explored for scenarios where the γ-rays are dominated by either 1) proton synchrotron emission or 2) synchrotron-self-Compton emission, with a subdominant but non-negligible contribution from photo-hadronic cascades in both cases. We find that the latter can be compatible with the neutrino observations, while the former is strongly disfavoured due to the insufficient neutrino production rate.

In 2017, the Event Horizon Telescope (EHT) Collaboration succeeded in capturing the first direct image of the center of the M87 galaxy. The asymmetric ring morphology and size are consistent with theoretical expectations for a weakly accreting supermassive black hole of mass ∼6.5 × 10⁹ M _⊙. The EHTC also partnered with several international facilities in space and on the ground, to arrange an extensive, quasi-simultaneous multi-wavelength campaign. This Letter presents the results and analysis of this campaign, as well as the multi-wavelength data as a legacy data repository. We captured M87 in a historically low state, and the core flux dominates over HST-1 at high energies, making it possible to combine core flux constraints with the more spatially precise very long baseline interferometry data. We present the most complete simultaneous multi-wavelength spectrum of the active nucleus to date, and discuss the complexity and caveats of combining data from different spatial scales into one broadband spectrum. We apply two heuristic, isotropic leptonic single-zone models to provide insight into the basic source properties, but conclude that a structured jet is necessary to explain M87’s spectrum. We can exclude that the simultaneous γ-ray emission is produced via inverse Compton emission in the same region producing the EHT mm-band emission, and further conclude that the γ-rays can only be produced in the inner jets (inward of HST-1) if there are strongly particle-dominated regions. Direct synchrotron emission from accelerated protons and secondaries cannot yet be excluded.

Context. Multiwavelength (MWL) observations of the blazar PKS 2155-304 during two weeks in July and August 2006, the period when two exceptional flares at very high energies (VHE, E ≳ 100 GeV) occurred, provide a detailed picture of the evolution of its emission. The complete data set from this campaign is presented, including observations in VHE γ-rays (H.E.S.S.), X-rays (RXTE, Chandra, Swift XRT), optical (Swift UVOT, Bronberg, Watcher, ROTSE), and in the radio band (NRT, HartRAO, ATCA). Optical and radio light curves from 2004 to 2008 are compared to the available VHE data from this period, to put the 2006 campaign into the context of the long-term evolution of the source.Aims. The data set offers a close view of the evolution of the source on different time scales and yields new insights into the properties of the emission process. The predictions of synchrotron self-Compton (SSC) scenarios are compared to the MWL data, with the aim of describing the dominant features in the data down to the hour time scale.Methods. The spectral variability in the X-ray and VHE bands is explored and correlations between the integral fluxes at different wavelengths are evaluated. SSC modelling is used to interpret the general trends of the varying spectral energy distribution.Results. The X-ray and VHE γ-ray emission are correlated during the observed high state of the source, but show no direct connection with longer wavelengths. The long-term flux evolution in the optical and radio bands is found to be correlated and shows that the source reaches a high state at long wavelengths after the occurrence of the VHE flares. Spectral hardening is seen in the Swift XRT data. Conclusions. The nightly averaged high-energy spectra of the non-flaring nights can be reproduced by a stationary one-zone SSC model, with only small variations in the parameters. The spectral and flux evolution in the high-energy band during the night of the second VHE flare is modelled with multi-zone SSC models, which can provide relatively simple interpretations for the hour time-scale evolution of the high-energy emission, even for such a complex data set. For the first time in this type of source, a clear indication is found for a relation between high activity at high energies and a long-term increase in the low frequency fluxes.

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 10³ to 10⁴, 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.

Context. The wealth of recent data from Imaging Air Cherenkov telescopes (IACTs), ultra-high energy cosmic-ray experiments and neutrino telescopes have fuelled a renewed interest in hadronic emission models for γ-loud blazars. Aims. We explore physically plausible solutions for a lepto-hadronic interpretation of the stationary emission from high-frequency peaked BL Lac objects (HBLs). The modelled spectral energy distributions are then searched for specific signatures at very high energies that could help to distinguish the hadronic origin of the emission from a standard leptonic scenario. Methods. By introducing a few basic constraints on parameters of the model, such as assuming the co-acceleration of electrons and protons, we significantly reduced the number of free parameters. We then systematically explored the parameter space of the size of the emission region and its magnetic field for two bright γ-loud HBLs, PKS 2155-304 and Mrk 421. For all solutions close to equipartition between the energy densities of protons and of the magnetic field, and with acceptable jet power and light-crossing timescales, we inspected the spectral hardening in the multi-TeV domain from proton-photon induced cascades and muon-synchrotron emission inside the source. Very-high-energy spectra simulated with the available instrument functions from the future Cherenkov Telescope Array (CTA) were evaluated for detectable features as a function of exposure time, source redshift, and flux level. Results. A range of hadronic scenarios are found to provide satisfactory solutions for the broad band emission of the sources under study. The TeV spectrum can be dominated either by proton-synchrotron emission or by muon-synchrotron emission. The solutions for HBLs cover a parameter space that is distinct from the one found for the most extreme BL Lac objects in an earlier study. Over a large range of model parameters, the spectral hardening due to internal synchrotron-pair cascades, the “cascade bump”, should be detectable for acceptable exposure times with the future CTA for a few nearby and bright HBLs.