Next Article in Journal
Neutron Star in Quantized Space-Time
Next Article in Special Issue
Science with the ASTRI Mini-Array: From Experiment to Open Observatory
Previous Article in Journal
VLBI Analysis of a Potential High-Energy Neutrino Emitter Blazar
Previous Article in Special Issue
Gamma-ray Bursts: 50 Years and Counting!
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Highlights of the Magic Florian Goebel Telescopes in the Study of Active Galactic Nuclei

by
Marina Manganaro
and
Dijana Dominis Prester
*,†
Faculty of Physics, University of Rijeka, Radmile Matejčić 2, 51000 Rijeka, Croatia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Universe 2024, 10(2), 80; https://doi.org/10.3390/universe10020080
Submission received: 8 December 2023 / Revised: 25 January 2024 / Accepted: 26 January 2024 / Published: 6 February 2024
(This article belongs to the Special Issue Recent Advances in Gamma Ray Astrophysics and Future Perspectives)

Abstract

:
The MAGIC (Major Atmospheric Gamma-ray Imaging Cherenkov) Florian Goebel telescopes are a system of two Cherenkov telescopes located on the Canary island of La Palma (Spain), at the Roque de Los Muchachos Observatory, which have been operating in stereo mode since 2009. Their low energy threshold (down to 15 GeV) allows the investigation of Active Galactic Nuclei (AGNs) in the very-high-energy (VHE, E > 100 GeV) gamma-ray range with a sensitivity up to the redshift limit of the existing IACT (Imaging Atmospheric Cherenkov Telescopes) systems. The MAGIC telescopes discovered 36 extragalactic objects emitting VHE gamma-rays and performed comprehensive studies of galaxies and their AGNs, also in a multi-wavelength (MWL) and multi-messenger (MM) context, expanding the knowledge of our Universe. Here, we report on the highlights achieved by the MAGIC collaboration since the beginning of their operations.

1. Introduction

The main targets of high-energy extragalactic astrophysics are Active Galactic Nuclei (AGNs), galaxies powered by a central supermassive black hole ( M 10 6 M , where M is the mass of the black hole and M 2 × 10 30 Kg is the solar mass) emitting non-thermal photons across the electromagnetic spectrum. AGNs whose radiation reaches the gamma-ray energy range account for only 10% of the AGNs observed in our Universe, and they exhibit collimated jets of relativistic plasma. The presence of the jet is closely linked to gamma-ray emission, and when such a jet is oriented in the direction of the observer, AGNs are called blazars and are further subdivided into BL Lac-type objects (BL Lacs) and Flat Spectrum Radio Quasars (FSRQs). AGNs can show variability in all energy bands, and due to the many studies and observations that have been conducted in each band, categorizing such objects can be challenging. When AGNs emit from the radio to the VHE (very-high-energy, E > 100 GeV) gamma-ray band, they present a double-peaked spectral energy distribution (SED) (e.g., [1]) that can be studied to identify and predict the emission mechanism of the radiation. The first bump, which extends between radio and X-ray energies, is universally attributed to synchrotron radiation emanating from a population of electrons rotating in the magnetic field of the AGN, while the second bump, which can reach VHE gamma-ray energies, can be explained, depending on the position and on its shape, by different theoretical models which can be fully leptonic or including an hadronic component. In the synchrotron self-Compton model (SSC, [2,3]), the high-energy bump is attributed to the inverse Compton scattering of lower-energy photons by relativistic electrons and is therefore often referred to as the IC bump. The SSC model is a leptonic model that successfully describes the MWL SED of most blazars, but in some cases it is necessary to consider an external component responsible for the emission of high-energy (HE, 100 MeV < E < 100 GeV) or VHE gamma rays, as in leptonic EC (external Compton) models (i.e., [4]). In the case of FSRQs, for example, the external radiation fields can be the accretion disk, broad-line region (BLR) or a dusty torus. In hadronic models, the high-energy gamma rays can also be produced by relativistic protons, by proton-synchrotron radiation or by photopion production. For more details on the theoretical models, see for example a review of the theoretical challenges in the study of AGNs [5,6] for a review of the results of high energy gamma-ray observations in the context of radiative mechanisms and models.
The subdivision of blazars into the subclasses BL Lacs and FSRQs is based on the absence or presence of emission/absorption lines in their optical spectra. BL Lacs have faint lines while the optical spectra of FSRQs have strong resolved lines that suggest the presence of radiatively efficient accretion disks and that can be measured providing a value for the luminosity and consequently the mass of their black-hole [7]. BL Lacs can be further categorized according to the position of the peak of their synchrotron bump, which can be measured by their broadband SED. They are defined as LBL (low-energy peaked), HBL (high-energy peaked) or IBL (intermediate) if their synchrotron peaks are in the submillimeter to infrared energy bands, in the ultraviolet to X-ray energy bands, or in the middle between the two aforementioned cases, respectively.
An effective recent categorization consists of dividing AGNs into two main classes, jetted and non-jetted AGNs [8]. This approach is very effective as it decouples the study of AGNs from features attributed to a single energy band and projects them in a multi-wavelength (MWL) and multi-messenger (MM) context. It is therefore of utmost importance to collect simultaneous data in a MWL context and investigate the time scale of the variability at different energies, correlations between the energy bands, polarization and periodicities, in order to derive a theoretical interpretation of the dataset and shed light on the emission scenarios of AGNs. In this context, the most important goals are:
  • to determine the location of the gamma-ray and VHE gamma-ray emission region in blazars;
  • to understand the composition of blazar jets and the particle population responsible for the observed high-energy emission;
  • to study the temporal evolution of MWL SEDs;
  • the search for a comprehensive explanation for the correlations and anti-correlations observed between different wavebands;
  • the investigation of the causes of variability.
The observation of AGNs has always been one of the priorities in the physics programme of the MAGIC Florian Goebel telescopes (Figure 1).
An overview of the latest highlights of the MAGIC Florian Goebel telescopes, which concern not only AGNs but also other targets of interest in this field (such as binary star systems, gamma-ray bursts and pulsars), can be found in [9].
In this work, Section 2 is devoted to MAGIC observations of AGNs, first presenting the peculiarities and characteristics of the MAGIC system (see Section 2.1) and then the observational strategies used (see Section 2.2). Section 3 collects the results highlights obtained with the MAGIC telescopes concerning the study of AGNs and Section 4 gives a brief overview of the results of other IACTs (Imaging Atmospheric Cherenkov Telescopes). In Section 5 we summarise the results presented.

2. AGNs Studies with MAGIC

2.1. MAGIC Characteristics

The MAGIC collaboration began its activities and data taking with a single telescope, now called MAGIC I, which has been in operation at the Roque de Los Muchachos Observatory on the Canary Island of La Palma since 2004. Since 2009, a second telescope, MAGIC II, has been part of the system and the two telescopes operate simultaneously (stereo mode). MAGIC telescopes are IACTs, ground-based instruments that detect VHE gamma rays thanks to the extensive air showers they generate in the atmosphere. The Cherenkov light emitted by the charged particles in the extensive air showers is captured by the segmented mirror dishes of the IACTs and imaged by their cameras. The IACT technique is described in detail in [10,11,12]. The MAGIC Florian Goebel telescopes were built with the aim of lowering the energy threshold for the detection of VHE gamma rays, in order to explore more distant sources in this energy range and thus study earlier parts of the Universe. At the time MAGIC I was built, the existing IACTs operated with an energy threshold ≥ 300 GeV [13,14]. A lower energy threshold was also necessary to close the gap to the HE gamma-ray detectors, as the latter reach an energy of a few tenths of a GeV. The connection to HE gamma-ray data constrains the HE bump of broadband SEDs and helps to effectively determine the emission scenarios of jetted AGNs. This goal has been successfully achieved, making MAGIC the IACT system with a lower energy threshold, 50 GeV, but capable of going down to 15 GeV using a dedicated stereoscopic analog trigger (SumTriggerII, see [15]), as was the case in MAGIC’s recent work on Geminga observations [16]. Important innovations were introduced in the design and construction of the telescopes compared to the other IACTs. Techniques known in accelerator particle physics experiments were used, such as fast electronics and automatic control of the instruments, as well as computers and networks capable of recording and reconstructing large volumes of data and performing interrelations. The main characteristics of MAGIC telescopes (see for details [17]) are as follows:
  • Active mirror surface of 236 m2, made of square elements 49.5 × 49.5 cm or 99 × 99 cm; f/D (focal length to diameter ratio) = 1.03;
  • Support frame made of reinforced carbon fibre tubes (<70 tons);
  • Approximately hexagonal camera with a diameter of 1.05 m, with 1039 PMTs of 1 inch diameter each (some PMTs of 2 inches diameter in the MAGIC I camera, see [17] for details); all PMTs have an effective quantum efficiency of 25 to 35%, depending on the wavelength; The camera is kept as light as possible and is held by an aluminium support arch, stiffened by a net of thin steel cables;
  • The maximum repositioning speed is more than 7 degrees per second, which means that the telescopes can be pointed to any point on the observable sky in less than 25 s;
  • Analogue signals are transmitted from the camera to the counting house via optical fibres; only the amplifiers and laser diode modulators for the transmission are located in the camera housing;
  • Digitization is performed by the Domino Ring Sampler (DRS4) chip with a sampling frequency of 1.64 GHz to use the timing information in the pulse.
The performance of the MAGIC Florian Goebel telescopes is described in detail in [17,18].

2.2. Observational Strategies

Sophisticated observation strategies are used to make the observation time of the MAGIC telescopes efficient. The aim is to maximise the observing time assigned to the various scientific projects within the collaboration and to fully exploit the potential of the telescopes. The MAGIC observational cycles have a duration of one year and the allocation of observing time is based on the submission of observational proposals, which can also be carried out by external scientists (https://magic.mpp.mpg.de/public/magicop/, accessed on 12 January 2024). As far as the observation of AGNs is concerned, MAGIC strategies can be summarised in two types of approaches:
  • ToO (Target of Opportunity) observations;
  • monitoring with short or long cadence.
MAGIC has a long history of ToO observations that began with the start of operations [19]. ToOs are quickly organised to respond to alerts reporting increased flux from other instruments, for example from optical or HE gamma-ray telescopes. Collaboration with other instruments and facilities is key to the prompt issuance of ToO observations and the resulting rapid re-planning of targets. The collection of MWL simultaneous data allows the target sources to be investigated in a broadband context, and in the case of a non-detection, the upper limits obtained by MAGIC can still help to describe the HE part of the broadband SED and facilitate the theoretical interpretation. Most of MAGIC’s detections in the VHE gamma-ray range, listed in Table 1, were made within the ToO MAGIC programme. Especially in the case of rapid variability, the timely response to alerts and the rapid rescheduling play a fundamental role.
The monitoring of specific targets has also yielded several results, as indicated in Section 3.8. The cadence of observation may vary depending on the characteristics of the target. Some AGNs are bright enough in the VHE gamma-ray range to be detected even when in a quiescent or low state, such as Mrk 421 and Mrk 501, allowing long-term monitoring with regularly repeated quick snapshots of observations, others are faint in the VHE gamma-ray range and require longer time windows for each scheduled observation. In any case, MWL campaigns are foreseen and organised in collaboration with other telescopes and facilities to obtain simultaneous data.
The observational proposals also consider observations at high zenith [20,21] and observations under moon conditions [22] for specific targets.
To reduce systematic uncertainties arising from atmospheric conditions, the MAGIC collaboration operates an elastic LIDAR (LIght Detection And Ranging) system [23]. This enables data corrections based on atmospheric transmission profiles [24].

3. Results

3.1. MAGIC Discoveries in the VHE Gamma-Ray Range

Since the beginning of operations, MAGIC has detected many AGNs emitting in the VHE gamma-ray range, 26 BL Lacs, six FSRQs (see Section 3.3), one radio galaxy (see Section 3.6) and three blazars (no final classification yet). Table 1 lists the AGNs discovered by MAGIC in the VHE gamma-ray range, ordered by their redshift.
Table 1. List of AGNs which were discovered to emit VHE gamma rays by the MAGIC telescopes. Data retrieved from TeVCat 1.
Table 1. List of AGNs which were discovered to emit VHE gamma rays by the MAGIC telescopes. Data retrieved from TeVCat 1.
NameTypeRedshiftDate of AnnouncementReferences
RGB J2042+244HBL0.1042019.11[25]
Mrk 180HBL0.0452006.09[26]
TXS 0210+515HBL0.0492019.01[25]
1ES 2037+521HBL0.0532016.10[25]
1ES 1727+502HBL0.0552011.11[27]
2WHSP J073326.7+515354HBL0.0652018.04[28]
1ES 1741+196HBL0.0842011.08[29]
B2 1811+31IBL0.1172020.10[30]
B3 2247+381HBL0.11872010.10[31]
TXS 1515-273HBL0.12842019.02[32]
1ES 1215+303HBL0.1312011.01[33]
RX J1136.5+6737HBL0.13422014.04[34]
1RXS J081201.8+023735HBL0.17212021.02([35] (video))
MAGIC J2001+435IBL0.17392010.07[36]
1ES 1218+304HBL0.1822006.05[37]
IC 310AGN (radio galaxy)0.01892010.03[38]
RBS 0723HBL0.1982014.01[25]
1ES 1011+496HBL0.2122007.09[39,40,41,42]
MS 1221.8+2452HBL0.2182013.05[43]
RGB J0136+391HBL>0.272012.07[44]
H 1722+119HBL 2013.05[45]
1ES 0647+250HBL>0.292010.07[46]
PKS 1413+135Blazar0.247 < z < 0.5 [47]2022.01[48]
S5 0716+714IBL0.26 [49]2008.04[50,51]
OT 081LBL0.3222016.07[52]
TXS 0506+056Blazar0.33652017.10[53,54]
S2 0109+22IBL0.362015.07[55]
S4 0954+65Blazar0.36942015.02[56]
PKS 1222+216FSRQ0.4322010.06[57]
1ES 0033+595HBL0.4672011.10[58]
GB6 J1058+2817BL Lac (class unclear)0.4793 [59]2021.04[60]
3C 279FSRQ0.53622008.06[61,62,63]
B2 1420+32FSRQ0.6822020.01[64]
TON 0599FSRQ0.72472017.12[65]
PKS 1441+25FSRQ0.9392015.04[66]
QSO B0218+357FSRQ0.9542014.07[67,68]
1 The TeVCat online source catalog, http://tevcat.uchicago.edu; (accessed on 12 January 2024).

3.2. Sources at High Redshift and EBL Studies

VHE gamma rays can interact with EBL photons by pair production. This process leads to an attenuation of the flux of VHE photons, which becomes more significant with increasing distance. For this reason, IACTs have a limit to the redshift they can achieve. This depends on the EBL properties and the energy of the VHE gamma rays emitted by the source under investigation. This limit is represented by the so-called gamma-ray horizon, which is described in detail in [69,70]. Theoretical predictions made in [69] clearly show the importance of a low energy threshold for IACTs to reach sources that are as far away as possible considering EBL attenuation. Given these predictions and the measurements of the EBL made with UV and the mid-infrared telescopes [71], the detection of AGNs near redshift 1 was considered extremely difficult. Nevertheless, thanks to their particularly low energy threshold, the MAGIC telescopes succeeded in detecting to AGNs close to redshift 1 in the VHE gamma-ray range, thus pushing the limits of the gamma-ray horizon. These two discoveries were particularly important for the study of EBL with IACTs. Both sources were detected in 2014 and were both FSRQs. Only recently an even more distant FSRQ, OP 313 with a redshift of 0.99, was detected by LST-1 (Large-Sized Telescope [72,73]).

3.2.1. QSO B0218+357

QSO B0218+357, with a redshift of z = 0.94, was also the first gravitationally lensed blazar detected in the VHE gamma ray range, triggered by increased activity in the HE gamma-ray range. Its discovery by MAGIC occurred during the arrival time of the delayed component of the emission [67].
The MWL SED challenged a simple leptonic model and a two-zone leptonic model with an external component (emission region located inside or outside the broad line region) was used to interpret the broadband emission. This model is described in [74] (third scenario). After this peculiar detection, a MWL campaign was organized to collect more information about the emission scenario, even if during a quiescent state. Also in this case, a model with an external Compton component and two zones was necessary to interpret the data [68]. It is interesting to note the different MWL SEDs observed in the first detection case (Figure 2a) and in the quiescent state monitoring (Figure 2b).

3.2.2. PKS 1441+25

PKS 1441+25 is a FSRQ with a redshift of z = 0.94 and was detected by MAGIC during a flaring state, triggered by high activity in the HE gamma-ray range. Its MWL SED could also not be interpreted by a simple leptonic model and an external Compton was considered, with the location of the emitting region explained as originating from the jet outside the broad line region [66]. With these distant sources, it was possible to study and measure the EBL for the first time at such a high redshift with VHE gamma-ray data, constraining the EBL density between 0.21 and 1.13 μ m, as reported in [75].

3.3. Flat Spectrum Radio Quasars

FSRQs are highly luminous blazars that show strong emission lines in their optical spectra. This feature indicates the presence of a radiatively efficient accretion disk [6]. The first classification of blazars based on their gamma-ray emission, where one of the properties of FSRQs is their soft gamma index ( Γ > 2 ), is the famous blazar sequence concept presented in [76]. Subsequent discoveries in the VHE gamma-ray range provided more material to investigate the classification of blazars and consolidate aspects of the blazar sequence while critically reconsidering and extending some of the initial considerations. A full overview on these topics can be found in [77]. Only 10 FSRQs have been confidently detected in the VHE gamma-ray range, six of which were discovered by MAGIC. MAGIC detected in 2006 the absolute first FSRQs in VHE energy, 3C 279, at a redshift of z = 0.53 [61,62]. This result was of utmost importance as it demonstrated the ability of IACTs to reach regions of the Universe where EBL attenuation is not negligible. Moreover, the emission of FSRQs in the VHE gamma-ray range leads to an MWL SED that exhibits a high Compton Dominance (CD, ratio between the so-called IC peak and the synchrotron peak). This latter property leads to a more sophisticated modelling approach with respect to BL Lacs. The high CD observed in the FSRQs detected in VHE gamma rays, was explored for the first time in detail with 3C 279. As shown in Figure 3, the MWL SED exhibits a high CD and models with an external component were required to describe the data set [62].
The second FSRQ detected by MAGIC was PKS 1222+216 in 2010 [57], and the VHE gamma-ray emission was found to be variable on a timescale of 10 min. This was an important result that made it possible to constrain the size of the gamma-ray emission region and to strongly constrain the emission models of blazar jets. The detection of PKS 1222+216 by MAGIC in fact challenged the models of jet emission and opened the discussion between different possible scenarios (far dissipation [78]; recollimation [79]; small compact emission regions in the jet [80]; reconfinement shocks [81]; relativistic filamentation [82]). The most distant FSRQs detected by MAGIC, already mentioned in Section 3.2, were detected in 2014 (QSO B0218+357) and 2015 (PKS 1441+25). After them, MAGIC detected another FSRQ in 2017, TON 0599, which lies at a redshift 0.72 [65]. This source also shows a MWL SED of complex interpretation which is the subject of a paper in preparation for the MAGIC collaboration. In 2020, MAGIC detected a new FSRQ, B2 1420+32 [64], adding another piece of the puzzle to the interpretation of such rare and powerful sources in the VHE gamma-ray sky.

3.4. Transitional Blazars

In some cases the usual categorization of blazars based on the properties of their optical spectra is complicated. BL Lacs are supposed to present an optical spectrum with very weak or absent emission lines. Instead in some cases (even for BL Lacertae), the optical spectra measured in different activity states can vary, making the categorization uncertain. Some blazars can be detected in the VHE gamma-ray range even when they are in quiescent state, while others can only be detected during flaring states of activity. Thus, it happens that some sources categorized as BL Lacs based on their optical spectra, suddenly show strong VHE gamma-ray emission in a flaring state that cannot be explained by the simplest leptonic models, which are often precise enough to describe the MWL emission of BL Lacs. Another categorization is related to the position of the synchrotron peak, which can also vary depending on the activity of the source. These sources could be called transitional sources [83], although it is not clear whether their behavior is simply related to exceptional gamma-ray activity in or to a real change in properties over time. To better understand these transitional AGNs, it is important to obtain a long-term dataset in MWL. MAGIC observed and detected for the first time in the VHE gamma-ray range two AGNs categorized as BL Lacs, but from the MWL dataset collected simultaneously with the MAGIC observations, they were found to have properties close to those of FSRQs. The first, S4 0954+65, was observed by MAGIC in 2015 and studied in an MWL context in [56]. It was found that the MWL SED of S4 0954+65 is reproduced by a leptonic model typical of FSRQs with an external Compton component (dusty torus). Similarly, the blazar OT 081, observed and detected in 2016 by MAGIC together with H.E.S.S. shows features of FSRQs rather than BL Lacs, in particular the high CD [52].

3.5. Extreme Sources

The BL Lacs that exhibit a synchrotron peak at unusually high energies, above 10 17 Hz, belong to the EHBL class (extreme HBLs, [84]). EHBLs are expected to be very faint in the VHE gamma-ray range and very difficult to detect by IACTs. Nevertheless, their study is appealing as they can be used for testing the gamma-ray propagation at high energy and in particular to probe EBL and derive limits on the IGMF (inter galactic magnetic field, see [85] for a review). Ref. [86] reports in detail on the progress made in studying such AGNs. MAGIC has devoted many observations to the hunting and study of EHBLs.
In 2018, MAGIC observed and detected the EHBL 2WHSP J073326.7+515354 for the first time in the VHE gamma-ray range [28]. The observations were scheduled as part of MAGIC’s hunting strategy for such extreme sources, and this source was chosen primarily because of its high synchrotron peak frequency ( ν s y n c h = 10 17.9 Hz). The successful strategy led to the VHE detection and to an in-depth study of the MWL properties of the SED. The theoretical interpretation of the broadband emission was challenging as expected for this type of AGNs, and four different theoretical models were used to test this particular case. The scenario that best described the dataset was a spine-layer two-zone leptonic model, as described in [87]. Many other EHBLs were observed by MAGIC as part of the EHBL hunting programme, and are collected in [25]. This work also includes the archetypal EHBL 1ES 0229+20 and the results of its observation from 2013 to 2017. Also included in [25] are three EHBLs detected for the first time in the VHE gamma-ray range: 1ES 2037+521, RBS 0723, and TXS 0210+515. Ref. [86] describes and emphasises the existence of two different types of EHBLs: extreme-synchrotron sources (synchrotron peak energy = h ν s y n c h 1 keV) and extreme-TeV sources (gamma-ray peak energy = h ν γ 1 TeV). This difference corresponds to a hard spectrum in the soft X-ray band (photon index Γ X r a y < 2 ), or in the TeV band ( Γ γ < 2 ). Blazars belonging to the HBL class can take on extreme properties and exhibit extreme behavior. This has been established on two occasions by MAGIC observations, for Mrk 501 [88] and for 1ES 2344+514 [89,90]. EHBLs, and in particular extreme-TeV sources, are perfect candidates for the study of IGMF, since their high energy emission and hard spectrum can be used to constrain the presence of cascades in the IGMF. For this reason, an in-depth study using gamma-ray observations (by MAGIC, H.E.S.S., VERITAS and the Fermi-LAT telescopes) of the archetypal EHBL 1ES 0229+20 was performed with the specific aim of detecting or constraining the IGMF-dependent secondary produced during the propagation of TeV gamma rays through the intergalactic medium. The results presented in [91] set robust limits consisting of a lower bound of B > 1.8 × 10 17 G for the long-correlation-length IGMF and B > 10 14 G for an IGMF of cosmological origin.

3.6. Black-Hole Lightening: IC310

A very interesting result from MAGIC concerns the radio galaxy IC 310, which is powered by a supermassive black hole ( M = 3 × 10 8 M ). In 2012, a high amount of activity of the source was observed in the VHE gamma-rays, characterized by a very fast variability (doubling time scales faster than 4.8 min). This result challenged the existing theoretical models of the variability and suggested a new interpretation of the sub-horizon scale variability consisting of a pulsar-like particle acceleration by the electric field across a magnetospheric gap at the base of the radio jet [92]. This interpretation is shown in Figure 4.

3.7. Markarian Galaxies at VHE Gamma-Rays: Mrk 421 and Mrk 501

Two AGNs emitting in the VHE gamma-ray range, Mrk 421 and Mrk 501, belonging to the category of Markarians (AGNs named after Benjamin Markarian, who brought them into attention because of their uncommon excess of emission in the UV band), have been monitored by MAGIC since the beginning of the operations. Those close-by blazars (redshift z = 0.031 and z = 0.034, respectively) are perfectly suited to study the mechanism of acceleration and broadband emission of blazars for the following reasons: their proximity implies a negligible effect of EBL absorption on their VHE gamma-ray spectra and also very precise VLBI (Very Large Baseline Interferometry) studies that can be used to constrain the emission scenarios; their brightness allows monitoring in different activity states (quiescent, flaring and intermediate states) and consequently a study of the temporal evolution of their broadband SEDs; their variability can be used to break degeneracies between emission models.
MWL campaigns were organized to obtain simultaneous MWL data from Mrk 421 and Mrk 501. This has led to several works over the years describing the broadband emission of such powerful blazars over the years, accompanied by detailed studies of the broadband correlations and variability in different activity states, as shown in Table 2 and Table 3.
All publications by MAGIC (in collaboration with many other instruments) on Mrk 421 and Mrk 501, together with the main results and the corresponding references, are listed in Table 2 and Table 3. Another Markarian, Mrk 180, was discovered in the VHE gamma-rays by MAGIC in 2006 [26] during an optical outburst. This source is very faint has only been observed in the VHE gamma-ray range by MAGIC on this occasion.

3.8. Long-Term Monitoring Campaigns

Long-term monitoring of specific AGN targets has been a successful strategy that has yielded many results.
For example, the giant radio galaxy M 87 has been studied by MAGIC over the years allowing in-depth characterization of its quiescent states [115,116], and investigation of its flare activity [117,118]. Since M 87 is a very important target for radio observation, a collaboration in MWL with many telescopes, including EHT (Event Horizon Telescope, ref. [119]), enabled an in-depth MWK study of this close radio galaxy [116,120].
The blazar PG 1553+113 was also the subject of long monitoring campaigns in MAGIC. This blazar has been observed by MAGIC from the beginning of the operations in mono mode [121,122] and as part of MWL monitoring campaigns. Ref. [123] reports on 5 years of observations and the associated results. Recently, this source, already known for its variability in different energy bands, has been found to exhibit a quasi-periodicity in gamma-rays [124,125,126], making its monitoring even more important to shed light on the complex MWL emission mechanism.
The BL Lac object 1ES 0647+250 was discovered in 2012 as emitting VHE gamma-rays [127] by MAGIC. Following this finding, a monitoring campaign was launched to collect MWL data and investigate this HBL in detail. The campaign resulted in an in-depth study of the correlations between the wavebands and the characterization of the broadband emission in four different activity states [46]. Long-term variability was found in the X-ray and VHE gamma-rays band, as well as correlations between radio, optical and HE gamma-ray fluxes.
The blazar 1ES 1959+650 was detected by MAGIC in 2004 [128]. Blazar 1ES 1959+650 is a nearby AGN (z = 0.048) and an HBL that showed high activity in various energy bands in 2015, especially in the optical, but also in the gamma-ray range. The results of MAGIC and MWL monitoring of this source in recent years are summarised in [129]. These include the detection of intra-night flux variability in the VHE gamma-ray range and the interpretation of the broadband emission within a one-zone SSC model.

3.9. Multi-Messenger Studies

The first MM event in which photons and neutrinos were detected simultaneously from an astrophysical source was the famous explosion of supernova 1987A [130].
In the following years, the development of neutrino experiments made it easier to detect neutrinos from astrophysical objects and to identify their origin. This enabled a complex system of alerts in cooperation with other observatories and telescopes at different photon energies. The efforts of the astronomical community were rewarded when in 2018 when the blazar TXS 0506+056 was found to be flaring in many wavelengths in coincidence with the neutrino event IceCube-170922A. This important detection was made possible by the coordination between the instruments and the prompt response to the IceCube alert. MAGIC was involved in the MM observations and in this discovery, reported in [53]. Several theoretical models have been tested to explain the MM SED and in [54] the dataset is interpreted using a one-zone lepto-hadronic scenario in which the gamma-ray emission is produced by the inverse Compton up-scattering of external photons (see Figure 5).
Thanks to the close collaboration between neutrino experiments and other instruments across the electromagnetic spectrum, MM observations are carried out regularly to be prepared for the next MM event of a blazar.

4. Broader Context: Other IACTs

Since MAGIC began the operations, other IACTs have also been active in the VHE gamma-ray field and have collected crucial results on AGNs and other targets. The present work is devoted only to the highlights of MAGIC, but it is important to consider the broader context of telescopes observing in a very similar energy range.
The VERITAS (Very Energetic Radiation Imaging Telescope Array System, ref. [131]) telescope array, which has been active since 2007, has observed and studied several AGNs as part of a dedicated AGN observing program [132]. The energy range of the VERITAS array extends from 85 GeV to 30 TeV. Scientific highlights of VERITAS are presented in [133] and listed on the VERITAS webpage (https://veritas.sao.arizona.edu/the-science-of-veritas/veritas-results; accessed on 12 January 2024). The VERITAS array is located in the Northern Hemisphere, like MAGIC, so many targets are common, allowing for joint publications and scientific collaboration in several cases (e.g., [88,97,100,103]).
The H.E.S.S. (High Energy Stereoscopic System) telescope array has been in operation since 2003. It is located in the Southern Hemisphere and has also studied several AGNs in the VHE gamma-ray range, from 30 GeV to 30 TeV. The highlights of their observations of AGNs are presented for example in [134,135]. A full list of scientific results can be found on their website (https://www.mpi-hd.mpg.de/HESS/pages/publications/, accessed on 12 January 2024). MAGIC has also collaborated with H.E.S.S. in the study of AGNs (e.g., [52,136]).
The FACT (First G-APD Cherenkov Telescope [137]), located at the same site of MAGIC, has been in operation since 2011 and its scientific program focuses mainly on the observations of extragalactic objects and in particular on the long-term monitoring of AGNs [138]. FACT observations can trigger follow-up observations with other IACTS with higher sensitivity, such as MAGIC, H.E.S.S. and VERITAS [88,89,104,113,139].
The LST-1 [140], which was inaugurated in 2018, is located at the same site of MAGIC (ORM observatory). This telescope, together with three others of the same design to be put in operation in the coming years, will form the group of LSTs [72] that will be part of the future CTAO (Cherenkov Telescope Array Observatory [141]). LST-1 has already achieved results in the observation of the gamma-ray source LHAASO J2108+5157 [142]. As for the study of AGNs, it has recently discovered the FSRQ OP 313 [73] in the VHE gamma-ray range. Collaboration with MAGIC in the observation of AGNs and other targets is planned, and joint data analysis will be possible thanks to a dedicated analysis pipeline [143].

5. Conclusions

The study of extragalactic sources has been one of the main targets of the MAGIC Florian Goebel telescopes for the last 20 years. The low energy threshold, which was the main goal behind their constructions (see Section 2.1), allowed a close connection to the other gamma-ray experiments and a smooth reconstruction of the high-energy bump of the MWL SEDs of many AGNs. In addition, it enabled the detection and discovery in the VHE gamma-ray range of two FSRQs close to redshift 1 (QSO B0218+357, Section 3.2.1 and PKS 1441+25, Section 3.2.2), at the limit of the IACT detection capacity, allowing the EBL to be studied for the first time at this distance with VHE gamma-ray measurements. Remarkably, of the only 10 FSRQs detected to date in the VHE gamma-ray range, six were discovered by MAGIC (See Section 3.3).
Thanks to the prompt response to MWL alerts under a dedicated ToO program, MAGIC discovered 36 AGNs in the VHE gamma-ray range (see Section 3.1) and characterized several types of AGNs in detail over the years, focusing on MWL and MM studies, which are of utmost importance for the identification of the broadband emission scenario. MAGIC was involved in the first MM observation of a neutrino blazar, TXS 0506+056, in which VHE gamma-rays were detected in spatial coincidence with the neutrino observed by IceCube (Section 3.9). The latter result suggested that blazars can be sources of high energy neutrinos and stimulated the development of theoretical models that can reproduce the MM SED.
The monitoring of AGNs such as Mrk 421 and Mrk 501, PG 1553+113, M 87, 1ES 0647+250, 1ES 1959+650 enabled long-term studies of correlations between different energy bands and variability (Section 3.7 and Section 3.8).
The results obtained by MAGIC during its long operation lifetime will certainly be complemented and extended by the new generation of IACTs—the CTAO (Cherenkov Telescope Array Observatory [141]).

Author Contributions

Conceptualization, M.M. and D.D.P.; Investigation, M.M. and D.D.P.; Writing—original draft preparation, M.M.; writing—review and editing, D.D.P. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Croatian Science Foundation (HrZZ) Project IP-2022-10-4595 and the University of Rijeka grant number uniri-prirod-18-48.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
AGNActive Galactic Nucleus
AGNsActive Galactic Nuclei
EBLExtragalactic Background Light
LIDARLIght Detection And Ranging
MAGICMajor Atmospheric Gamma-ray Imaging Cherenkov (telescopes)
IACTsImaging Atmospheric Cherenkov Telescopes
HEHigh-energy
VHEVery-high-energy
MWLMulti-wavelength
MMMulti-messenger
SEDSpectral Energy Distribution
FSRQFlat Spectrum Radio Quasar
BL LacsBL Lacertae type objects
f.o.vField of view
MJDModified Julian Date
ToOTarget of Opportunity
SSCSynchrotron self-Compton

References

  1. Ghisellini, G.; Righi, C.; Costamante, L.; Tavecchio, F. The Fermi blazar sequence. Mon. Not. R. Astron. Soc. 2017, 469, 255–266. [Google Scholar] [CrossRef]
  2. Konigl, A. Relativistic jets as X-ray and gamma-ray sources. Astrophys. J. 1981, 243, 700–709. [Google Scholar] [CrossRef]
  3. Maraschi, L.; Ghisellini, G.; Celotti, A. A Jet Model for the Gamma-Ray—Emitting Blazar 3C 279. Astrophys. J. Lett. 1992, 397, L5–L9. [Google Scholar] [CrossRef]
  4. Dermer, C.D.; Schlickeiser, R. Model for the High-Energy Emission from Blazars. Astrophys. J. 1993, 416, 458. [Google Scholar] [CrossRef]
  5. Böttcher, M. Progress in Multi-Wavelength and Multi-Messenger Observations of Blazars and Theoretical Challenges. Galaxies 2019, 7, 20. [Google Scholar] [CrossRef]
  6. Madejski, G.G.; Sikora, M. Gamma-Ray Observations of Active Galactic Nuclei. Annu. Rev. Astron. Astrophys. 2016, 54, 725–760. [Google Scholar] [CrossRef]
  7. Vestergaard, M.; Peterson, B.M. Determining Central Black Hole Masses in Distant Active Galaxies and Quasars. II. Improved Optical and UV Scaling Relationships. Astrophys. J. 2006, 641, 689–709. [Google Scholar] [CrossRef]
  8. Padovani, P. On the two main classes of active galactic nuclei. Nat. Astron. 2017, 1, 0194. [Google Scholar] [CrossRef]
  9. Loporchio, S.; MAGIC Collaboration. Highlights of the very-high-energy gamma-ray sky as seen by MAGIC. Nucl. Instrum. Methods Phys. Res. 2023, 1055, 168441. [Google Scholar] [CrossRef]
  10. Bose, D.; Chitnis, V.R.; Majumdar, P.; Acharya, B.S. Ground-based gamma-ray astronomy: History and development of techniques. Eur. Phys. J. Spec. Top. 2022, 231, 3–26. [Google Scholar] [CrossRef]
  11. de Naurois, M.; Mazin, D. Ground-based detectors in very-high-energy gamma-ray astronomy. Comptes Rendus Phys. 2015, 16, 610–627. [Google Scholar] [CrossRef]
  12. Aharonian, F.; Buckley, J.; Kifune, T.; Sinnis, G. High energy astrophysics with ground-based gamma ray detectors. Rep. Prog. Phys. 2008, 71, 096901. [Google Scholar] [CrossRef]
  13. Pühlhofer, G.; Bolz, O.; Götting, N.; Heusler, A.; Horns, D.; Kohnle, A.; Lampeitl, H.; Panter, M.; Tluczykont, M.; Aharonian, F.; et al. The technical performance of the HEGRA system of imaging air Cherenkov telescopes. Astropart. Phys. 2003, 20, 267–291. [Google Scholar] [CrossRef]
  14. Rovero, A.C.; Buckley, J.H.; Fleury, P.; Jiang, Y.; Pare, E.; Sarazin, X.; Urban, M.; Weekes, T.C. Calibration of the Whipple atmospheric Čerenkov telescope. Astropart. Phys. 1996, 5, 27–34. [Google Scholar] [CrossRef]
  15. Dazzi, F.; Schweizer, T.; Ceribella, G.; Corti, D.; Dettlaff, A.; Garcia, J.R.; Hafner, D.; Herranz, D.; Lopez-Moya, M.; Mariotti, M.; et al. The Stereoscopic Analog Trigger of the MAGIC Telescopes. IEEE Trans. Nucl. Sci. 2021, 68, 1473–1486. [Google Scholar] [CrossRef]
  16. MAGIC Collaboration; Acciari, V.A.; Ansoldi, S.; Antonelli, L.A.; Arbet Engels, A.; Asano, K.; Baack, D.; Babić, A.; Baquero, A.; Barres de Almeida, U.; et al. Detection of the Geminga pulsar with MAGIC hints at a power-law tail emission beyond 15 GeV. Astron. Astrophys. 2020, 643, L14. [Google Scholar] [CrossRef]
  17. Aleksić, J.; Ansoldi, S.; Antonelli, L.A.; Antoranz, P.; Babic, A.; Bangale, P.; Barceló, M.; Barrio, J.A.; Becerra González, J.; Bednarek, W.; et al. The major upgrade of the MAGIC telescopes, Part I: The hardware improvements and the commissioning of the system. Astropart. Phys. 2016, 72, 61–75. [Google Scholar] [CrossRef]
  18. Aleksić, J.; Ansoldi, S.; Antonelli, L.A.; Antoranz, P.; Babic, A.; Bangale, P.; Barceló, M.; Barrio, J.A.; Becerra González, J.; Bednarek, W.; et al. The major upgrade of the MAGIC telescopes, Part II: A performance study using observations of the Crab Nebula. Astropart. Phys. 2016, 72, 76–94. [Google Scholar] [CrossRef]
  19. Mazin, D.; Lindfors, E. Successful ToO triggers on the extragalactic sources with the MAGIC telescope. arXiv 2008, arXiv:0709.1694. [Google Scholar] [CrossRef]
  20. Aleksić, J.; Alvarez, E.A.; Antonelli, L.A.; Antoranz, P.; Asensio, M.; Backes, M.; Barres de Almeida, U.; Barrio, J.A.; Bastieri, D.; Becerra González, J.; et al. High zenith angle observations of PKS 2155-304 with the MAGIC-I telescope. Astron. Astrophys. 2012, 544, A75. [Google Scholar] [CrossRef]
  21. MAGIC Collaboration; Acciari, V.A.; Ansoldi, S.; Antonelli, L.A.; Arbet Engels, A.; Baack, D.; Babić, A.; Banerjee, B.; Barres de Almeida, U.; Barrio, J.A.; et al. MAGIC very large zenith angle observations of the Crab Nebula up to 100 TeV. Astron. Astrophys. 2020, 635, A158. [Google Scholar] [CrossRef]
  22. Ahnen, M.L.; Ansoldi, S.; Antonelli, L.A.; Arcaro, C.; Babić, A.; Banerjee, B.; Bangale, P.; Barres de Almeida, U.; Barrio, J.A.; Becerra González, J.; et al. Performance of the MAGIC telescopes under moonlight. Astropart. Phys. 2017, 94, 29–41. [Google Scholar] [CrossRef]
  23. Fruck, C.; Gaug, M.; Hahn, A.; Acciari, V.; Besenrieder, J.; Dominis Prester, D.; Dorner, D.; Fink, D.; Font, L.; Mićanović, S.; et al. Characterizing the aerosol atmosphere above the Observatorio del Roque de los Muchachos by analysing seven years of data taken with an GaAsP HPD-readout, absolutely calibrated elastic LIDAR. Mon. Not. R. Astron. Soc. 2022, 515, 4520–4550. [Google Scholar] [CrossRef]
  24. Schmuckermaier, F.; Gaug, M.; Fruck, C.; Moralejo, A.; Hahn, A.; Dominis Prester, D.; Dorner, D.; Font, L.; Mićanović, S.; Mirzoyan, R.; et al. Correcting Imaging Atmospheric Cherenkov Telescope data with atmospheric profiles obtained with an elastic light detecting and ranging system. Astron. Astrophys. 2023, 673, A2. [Google Scholar] [CrossRef]
  25. Acciari, V.A.; Ansoldi, S.; Antonelli, L.A.; Engels, A.A.; Asano, K.; Baack, D.; Babić, A.; Banerjee, B.; Barres de Almeida, U.; Barrio, J.A.; et al. New Hard-TeV Extreme Blazars Detected with the MAGIC Telescopes. Astrophys. J. Suppl. Ser. 2020, 247, 16. [Google Scholar] [CrossRef]
  26. Albert, J.; Aliu, E.; Anderhub, H.; Antoranz, P.; Armada, A.; Asensio, M.; Baixeras, C.; Barrio, J.A.; Bartko, H.; Bastieri, D.; et al. Discovery of Very High Energy γ-Rays from Markarian 180 Triggered by an Optical Outburst. Astrophys. J. Lett. 2006, 648, L105–L108. [Google Scholar] [CrossRef]
  27. Aleksić, J.; Antonelli, L.A.; Antoranz, P.; Asensio, M.; Backes, M.; Barres de Almeida, U.; Barrio, J.A.; Becerra González, J.; Bednarek, W.; Berger, K.; et al. Discovery of very high energy gamma-ray emission from the blazar 1ES 1727+502 with the MAGIC Telescopes. Astron. Astrophys. 2014, 563, A90. [Google Scholar] [CrossRef]
  28. MAGIC Collaboration; Acciari, V.A.; Ansoldi, S.; Antonelli, L.A.; Arbet Engels, A.; Baack, D.; Babić, A.; Banerjee, B.; Barres de Almeida, U.; Barrio, J.A.; et al. Testing emission models on the extreme blazar 2WHSP J073326.7+515354 detected at very high energies with the MAGIC telescopes. Mon. Not. R. Astron. Soc. 2019, 490, 2284–2299. [Google Scholar] [CrossRef]
  29. Ahnen, M.L.; Ansoldi, S.; Antonelli, L.A.; Antoranz, P.; Arcaro, C.; Babic, A.; Banerjee, B.; Bangale, P.; Barres de Almeida, U.; Barrio, J.A.; et al. MAGIC detection of very high energy γ-ray emission from the low-luminosity blazar 1ES 1741+196. Mon. Not. R. Astron. Soc. 2017, 468, 1534–1541. [Google Scholar] [CrossRef]
  30. Blanch, O. Detection of very-high-energy gamma-ray emission from B2 1811+31 with the MAGIC telescopes. Astron. Telegr. 2020, 14090, 1–11. [Google Scholar]
  31. Aleksić, J.; Alvarez, E.A.; Antonelli, L.A.; Antoranz, P.; Asensio, M.; Backes, M.; Barrio, J.A.; Bastieri, D.; Becerra González, J.; Bednarek, W.; et al. Discovery of VHE γ-ray emission from the BL Lacertae object B3 2247+381 with the MAGIC telescopes. Astron. Astrophys. 2012, 539, A118. [Google Scholar] [CrossRef]
  32. Acciari, V.A.; Ansoldi, S.; Antonelli, L.A.; Arbet Engels, A.; Artero, M.; Asano, K.; Baack, D.; Babić, A.; Baquero, A.; Barres de Almeida, U.; et al. First detection of VHE gamma-ray emission from TXS 1515-273, study of its X-ray variability and spectral energy distribution. Mon. Not. R. Astron. Soc. 2021, 507, 1528–1545. [Google Scholar] [CrossRef]
  33. Aleksić, J.; Alvarez, E.A.; Antonelli, L.A.; Antoranz, P.; Ansoldi, S.; Asensio, M.; Backes, M.; Barres de Almeida, U.; Barrio, J.A.; Bastieri, D.; et al. Discovery of VHE γ-rays from the blazar 1ES 1215+303 with the MAGIC telescopes and simultaneous multi-wavelength observations. Astron. Astrophys. 2012, 544, A142. [Google Scholar] [CrossRef]
  34. Mirzoyan, R. Discovery of Very High Energy Gamma-Ray Emission from BL Lac object RX J1136.5+6737 by the MAGIC Telescopes. Astron. Telegr. 2014, 6062, 1. [Google Scholar]
  35. Ventura, S.; Prandini, E.; Fallah Ramazani, V.; Jormanainen, J.; Bonnoli, G.; Da Vela, P. Very-High-Energy Detection of the Extreme Blazar 1RXS J081201.8+023735 with the MAGIC Telescopes. In Proceedings of the 43rd COSPAR Scientific Assembly, Sydney, Australia, 28 January–4 February 2021; Volume 43, p. 1281. [Google Scholar]
  36. Aleksić, J.; Ansoldi, S.; Antonelli, L.A.; Antoranz, P.; Babic, A.; Bangale, P.; Barres de Almeida, U.; Barrio, J.A.; Becerra González, J.; Bednarek, W.; et al. First broadband characterization and redshift determination of the VHE blazar MAGIC J2001+439. Astron. Astrophys. 2014, 572, A121. [Google Scholar] [CrossRef]
  37. Albert, J.; Aliu, E.; Anderhub, H.; Antoranz, P.; Armada, A.; Asensio, M.; Baixeras, C.; Barrio, J.A.; Bartelt, M.; Bartko, H.; et al. Discovery of Very High Energy Gamma Rays from 1ES 1218+30.4. Astrophys. J. Lett. 2006, 642, L119–L122. [Google Scholar] [CrossRef]
  38. Aleksić, J.; Antonelli, L.A.; Antoranz, P.; Babic, A.; Barres de Almeida, U.; Barrio, J.A.; Becerra González, J.; Bednarek, W.; Berger, K.; Bernardini, E.; et al. Rapid and multiband variability of the TeV bright active nucleus of the galaxy IC 310. Astron. Astrophys. 2014, 563, A91. [Google Scholar] [CrossRef]
  39. Albert, J.; Aliu, E.; Anderhub, H.; Antoranz, P.; Armada, A.; Baixeras, C.; Barrio, J.A.; Bartko, H.; Bastieri, D.; Becker, J.K.; et al. Discovery of Very High Energy γ-Rays from 1ES 1011+496 at z = 0.212. Astrophys. J. Lett. 2007, 667, L21–L24. [Google Scholar] [CrossRef]
  40. Ahnen, M.L.; Ansoldi, S.; Antonelli, L.A.; Antoranz, P.; Babic, A.; Banerjee, B.; Bangale, P.; Barres de Almeida, U.; Barrio, J.A.; Becerra González, J.; et al. MAGIC observations of the February 2014 flare of 1ES 1011+496 and ensuing constraint of the EBL density. Astron. Astrophys. 2016, 590, A24. [Google Scholar] [CrossRef]
  41. Ahnen, M.L.; Ansoldi, S.; Antonelli, L.A.; Antoranz, P.; Babic, A.; Banerjee, B.; Bangale, P.; Barres de Almeida, U.; Barrio, J.A.; Becerra González, J.; et al. Multiwavelength observations of the blazar 1ES 1011+496 in Spring 2008. Mon. Not. R. Astron. Soc. 2016, 459, 2286–2298. [Google Scholar] [CrossRef]
  42. Aleksić, J.; Ansoldi, S.; Antonelli, L.A.; Antoranz, P.; Arcaro, C.; Babic, A.; Bangale, P.; Barres de Almeida, U.; Barrio, J.A.; Becerra González, J.; et al. Insights into the emission of the blazar 1ES 1011+496 through unprecedented broadband observations during 2011 and 2012. Astron. Astrophys. 2016, 591, A10. [Google Scholar] [CrossRef]
  43. Cortina, J. Discovery of Very High Energy Gamma-Ray Emission from MS1221.8+2452 with the MAGIC telescopes. Astron. Telegr. 2013, 5038, 1. [Google Scholar]
  44. Mazin, D.; MAGIC Collaboration. Highlights from the MAGIC telescopes. In High Energy Gamma-Ray Astronomy: 5th International Meeting on High Energy Gamma-Ray Astronomy; Aharonian, F.A., Hofmann, W., Rieger, F.M., Eds.; American Institute of Physics Conference Series; AIP Publishing: Melville, NY, USA, 2012; Volume 1505, pp. 186–193. [Google Scholar] [CrossRef]
  45. Ahnen, M.L.; Ansoldi, S.; Antonelli, L.A.; Antoranz, P.; Babic, A.; Banerjee, B.; Bangale, P.; Barres de Almeida, U.; Barrio, J.A.; Becerra González, J.; et al. Investigating the peculiar emission from the new VHE gamma-ray source H1722+119. Mon. Not. R. Astron. Soc. 2016, 459, 3271–3281. [Google Scholar] [CrossRef]
  46. MAGIC Collaboration; Acciari, V.A.; Aniello, T.; Ansoldi, S.; Antonelli, L.A.; Arbet Engels, A.; Arcaro, C.; Artero, M.; Asano, K.; Baack, D.; et al. Long-term multi-wavelength study of 1ES 0647+250. Astron. Astrophys. 2023, 670, A49. [Google Scholar] [CrossRef]
  47. Readhead, A.C.S.; Ravi, V.; Liodakis, I.; Lister, M.L.; Singh, V.; Aller, M.F.; Blandford, R.D.; Browne, I.W.A.; Gorjian, V.; Grainge, K.J.B.; et al. The Relativistic Jet Orientation and Host Galaxy of the Peculiar Blazar PKS 1413+135. Astrophys. J. 2021, 907, 61. [Google Scholar] [CrossRef]
  48. Blanch, O.; Sitarek, J.; Striskovic, J. First detection of very-high-energy gamma-ray emission from PKS1413+135 with the MAGIC telescopes. Astron. Telegr. 2022, 15161, 1. [Google Scholar]
  49. Dorigo Jones, J.; Johnson, S.D.; Muzahid, S.; Charlton, J.; Chen, H.W.; Narayanan, A.; Schaye, J.; Wijers, N.A. Improving blazar redshift constraints with the edge of the Ly α forest: 1ES 1553+113 and implications for observations of the WHIM. Mon. Not. R. Astron. Soc. 2022, 509, 4330–4343. [Google Scholar] [CrossRef]
  50. Anderhub, H.; Antonelli, L.A.; Antoranz, P.; Backes, M.; Baixeras, C.; Balestra, S.; Barrio, J.A.; Bastieri, D.; Becerra González, J.; Becker, J.K.; et al. Discovery of very High Energy γ-Rays from the Blazar S5 0716+714. Astrophys. J. Lett. 2009, 704, L129–L133. [Google Scholar] [CrossRef]
  51. MAGIC Collaboration; Ahnen, M.L.; Ansoldi, S.; Antonelli, L.A.; Arcaro, C.; Baack, D.; Babić, A.; Banerjee, B.; Bangale, P.; Barres de Almeida, U.; et al. Multi-wavelength characterization of the blazar S5 0716+714 during an unprecedented outburst phase. Astron. Astrophys. 2018, 619, A45. [Google Scholar] [CrossRef]
  52. Manganaro, M.; Seglar-Arroyo, M.; Becerra-González, J.; Sanchez, D.; Cerruti, M.; Tavecchio, F.; Fallah-Ramazani, V.; Agudo, I.; Ciprini, S.; Filippenko, A.V.; et al. MAGIC and H.E.S.S. detect VHE gamma rays from the blazar OT081 for the first time: A deep multiwavelength study. In Proceedings of the 37th International Cosmic Ray Conference, Berlin, Germany, 12–23 July 2021; p. 815. [Google Scholar]
  53. IceCube Collaboration; Aartsen, M.G.; Ackermann, M.; Adams, J.; Aguilar, J.A.; Ahlers, M.; Ahrens, M.; Al Samarai, I.; Altmann, D.; Andeen, K.; et al. Multimessenger observations of a flaring blazar coincident with high-energy neutrino IceCube-170922A. Science 2018, 361, eaat1378. [Google Scholar] [CrossRef]
  54. Ansoldi, S.; Antonelli, L.A.; Arcaro, C.; Baack, D.; Babić, A.; Banerjee, B.; Bangale, P.; Barres de Almeida, U.; Barrio, J.A.; Becerra González, J.; et al. The Blazar TXS 0506+056 Associated with a High-energy Neutrino: Insights into Extragalactic Jets and Cosmic-Ray Acceleration. Astrophys. J. Lett. 2018, 863, L10. [Google Scholar] [CrossRef]
  55. MAGIC Collaboration; Ansoldi, S.; Antonelli, L.A.; Arcaro, C.; Baack, D.; Babić, A.; Banerjee, B.; Bangale, P.; Barres de Almeida, U.; Barrio, J.A.; et al. The broad-band properties of the intermediate synchrotron peaked BL Lac S2 0109+22 from radio to VHE gamma-rays. Mon. Not. R. Astron. Soc. 2018, 480, 879–892. [Google Scholar] [CrossRef]
  56. MAGIC Collaboration; Ahnen, M.L.; Ansoldi, S.; Antonelli, L.A.; Arcaro, C.; Baack, D.; Babić, A.; Banerjee, B.; Bangale, P.; Barres de Almeida, U.; et al. Detection of the blazar S4 0954+65 at very-high-energy with the MAGIC telescopes during an exceptionally high optical state. Astron. Astrophys. 2018, 617, A30. [Google Scholar] [CrossRef]
  57. Aleksić, J.; Antonelli, L.A.; Antoranz, P.; Backes, M.; Barrio, J.A.; Bastieri, D.; Becerra González, J.; Bednarek, W.; Berdyugin, A.; Berger, K.; et al. MAGIC Discovery of Very High Energy Emission from the FSRQ PKS 1222+21. Astrophys. J. Lett. 2011, 730, L8. [Google Scholar] [CrossRef]
  58. Aleksić, J.; Ansoldi, S.; Antonelli, L.A.; Antoranz, P.; Babic, A.; Bangale, P.; Barres de Almeida, U.; Barrio, J.A.; Becerra González, J.; Bednarek, W.; et al. Discovery of very high energy γ-ray emission from the blazar 1ES 0033+595 by the MAGIC telescopes. Mon. Not. R. Astron. Soc. 2015, 446, 217–225. [Google Scholar] [CrossRef]
  59. Massaro, F.; Masetti, N.; D’Abrusco, R.; Paggi, A.; Funk, S. Optical Spectroscopic Observations of Blazars and γ-Ray Blazar Candidates in the Sloan Digital Sky Survey Data Release Nine. Astron. J. 2014, 148, 66. [Google Scholar] [CrossRef]
  60. Blanch, O. Detection of very-high-energy gamma-ray emission from GB6 J1058+2817 with the MAGIC telescopes. Astron. Telegr. 2021, 14506, 1. [Google Scholar]
  61. MAGIC Collaboration; Albert, J.; Aliu, E.; Anderhub, H.; Antonelli, L.A.; Antoranz, P.; Backes, M.; Baixeras, C.; Barrio, J.A.; Bartko, H.; et al. Very-High-Energy gamma rays from a Distant Quasar: How Transparent Is the Universe? Science 2008, 320, 1752. [Google Scholar] [CrossRef]
  62. Aleksić, J.; Antonelli, L.A.; Antoranz, P.; Backes, M.; Barrio, J.A.; Bastieri, D.; Becerra González, J.; Bednarek, W.; Berdyugin, A.; Berger, K.; et al. MAGIC Observations and multiwavelength properties of the quasar 3C 279 in 2007 and 2009. Astron. Astrophys. 2011, 530, A4. [Google Scholar] [CrossRef]
  63. Aleksić, J.; Ansoldi, S.; Antonelli, L.A.; Antoranz, P.; Babic, A.; Bangale, P.; Barres de Almeida, U.; Barrio, J.A.; Becerra González, J.; Bednarek, W.; et al. MAGIC observations and multifrequency properties of the flat spectrum radio quasar 3C 279 in 2011. Astron. Astrophys. 2014, 567, A41. [Google Scholar] [CrossRef]
  64. MAGIC Collaboration; Acciari, V.A.; Ansoldi, S.; Antonelli, L.A.; Arbet Engels, A.; Artero, M.; Asano, K.; Baack, D.; Babić, A.; Baquero, A.; et al. VHE gamma-ray detection of FSRQ QSO B1420+326 and modeling of its enhanced broadband state in 2020. Astron. Astrophys. 2021, 647, A163. [Google Scholar] [CrossRef]
  65. Mirzoyan, R. Detection of very-high-energy gamma-ray emission from the FSRQ Ton 0599 with the MAGIC telescopes. Astron. Telegr. 2017, 11061, 1. [Google Scholar]
  66. Ahnen, M.L.; Ansoldi, S.; Antonelli, L.A.; Antoranz, P.; Babic, A.; Banerjee, B.; Bangale, P.; Barres de Almeida, U.; Barrio, J.A.; Bednarek, W.; et al. Very High Energy γ-Rays from the Universe’s Middle Age: Detection of the z = 0.940 Blazar PKS 1441+25 with MAGIC. Astrophys. J. Lett. 2015, 815, L23. [Google Scholar] [CrossRef]
  67. Ahnen, M.L.; Ansoldi, S.; Antonelli, L.A.; Antoranz, P.; Arcaro, C.; Babic, A.; Banerjee, B.; Bangale, P.; Barres de Almeida, U.; Barrio, J.A.; et al. Detection of very high energy gamma-ray emission from the gravitationally lensed blazar QSO B0218+357 with the MAGIC telescopes. Astron. Astrophys. 2016, 595, A98. [Google Scholar] [CrossRef]
  68. Acciari, V.A.; Ansoldi, S.; Antonelli, L.A.; Arbet Engels, A.; Artero, M.; Asano, K.; Baack, D.; Babić, A.; Baquero, A.; Barres de Almeida, U.; et al. Multiwavelength study of the gravitationally lensed blazar QSO B0218+357 between 2016 and 2020. Mon. Not. R. Astron. Soc. 2022, 510, 2344–2362. [Google Scholar] [CrossRef]
  69. Blanch, O.; Martinez, M. Exploring the gamma ray horizon with the next generation of gamma ray telescopes. Part 1: Theoretical predictions. Astropart. Phys. 2005, 23, 588–597. [Google Scholar] [CrossRef]
  70. Blanch, O.; Martinez, M. Exploring the gamma-ray horizon with the next generation of gamma-ray telescopes. Part 2: Extracting cosmological parameters from the observation of gamma-ray sources. Astropart. Phys. 2005, 23, 598–607. [Google Scholar] [CrossRef]
  71. Domínguez, A.; Primack, J.R.; Rosario, D.J.; Prada, F.; Gilmore, R.C.; Faber, S.M.; Koo, D.C.; Somerville, R.S.; Pérez-Torres, M.A.; Pérez-González, P.; et al. Extragalactic background light inferred from AEGIS galaxy-SED-type fractions. Mon. Not. R. Astron. Soc. 2011, 410, 2556–2578. [Google Scholar] [CrossRef]
  72. Pierro, F.D. Status of the Large-Sized Telescope of the Cherenkov Telescope Array. J. Phys. Conf. Ser. 2023, 2429, 012020. [Google Scholar] [CrossRef]
  73. Cortina, J.; CTAO LST Collaboration. First detection of VHE gamma-ray emission from FSRQ OP 313 with LST-1. Astron. Telegr. 2023, 16381, 1. [Google Scholar]
  74. Tavecchio, F.; Becerra-Gonzalez, J.; Ghisellini, G.; Stamerra, A.; Bonnoli, G.; Foschini, L.; Maraschi, L. On the origin of the γ-ray emission from the flaring blazar PKS 1222+216. Astron. Astrophys. 2011, 534, A86. [Google Scholar] [CrossRef]
  75. Acciari, V.A.; Ansoldi, S.; Antonelli, L.A.; Arbet Engels, A.; Baack, D.; Babić, A.; Banerjee, B.; Barres de Almeida, U.; Barrio, J.A.; Becerra-Gonzalez, J.; et al. Measurement of the extragalactic background light using MAGIC and Fermi-LAT gamma-ray observations of blazars up to z = 1. Mon. Not. R. Astron. Soc. 2019, 486, 4233–4251. [Google Scholar] [CrossRef]
  76. Fossati, G.; Maraschi, L.; Celotti, A.; Comastri, A.; Ghisellini, G. A unifying view of the spectral energy distributions of blazars. Mon. Not. R. Astron. Soc. 1998, 299, 433–448. [Google Scholar] [CrossRef]
  77. Prandini, E.; Ghisellini, G. The Blazar Sequence and Its Physical Understanding. Galaxies 2022, 10, 35. [Google Scholar] [CrossRef]
  78. Tavecchio, F.; Ghisellini, G.; Bonnoli, G.; Ghirlanda, G. Constraining the location of the emitting region in Fermi blazars through rapid γ-ray variability. Mon. Not. R. Astron. Soc. 2010, 405, L94–L98. [Google Scholar] [CrossRef]
  79. Marscher, A.P. Relativistic jets and the continuum emission in QSOs. Astrophys. J. 1980, 235, 386–391. [Google Scholar] [CrossRef]
  80. Ghisellini, G.; Tavecchio, F. Rapid variability in TeV blazars: The case of PKS2155-304. Mon. Not. R. Astron. Soc. 2008, 386, L28–L32. [Google Scholar] [CrossRef]
  81. Nalewajko, K.; Sikora, M. A structure and energy dissipation efficiency of relativistic reconfinement shocks. Mon. Not. R. Astron. Soc. 2009, 392, 1205–1210. [Google Scholar] [CrossRef]
  82. Frederiksen, J.T.; Haugbølle, T.; Medvedev, M.V.; Nordlund, Å. Radiation Spectral Synthesis of Relativistic Filamentation. Astrophys. J. Lett. 2010, 722, L114–L119. [Google Scholar] [CrossRef]
  83. Ghisellini, G.; Tavecchio, F.; Foschini, L.; Ghirlanda, G. The transition between BL Lac objects and flat spectrum radio quasars. Mon. Not. R. Astron. Soc. 2011, 414, 2674–2689. [Google Scholar] [CrossRef]
  84. Costamante, L.; Ghisellini, G.; Giommi, P.; Tagliaferri, G.; Celotti, A.; Chiaberge, M.; Fossati, G.; Maraschi, L.; Tavecchio, F.; Treves, A.; et al. Extreme synchrotron BL Lac objects. Stretching the blazar sequence. Astron. Astrophys. 2001, 371, 512–526. [Google Scholar] [CrossRef]
  85. Durrer, R.; Neronov, A. Cosmological magnetic fields: Their generation, evolution and observation. Astron. Astrophys. Rev. 2013, 21, 62. [Google Scholar] [CrossRef]
  86. Biteau, J.; Prandini, E.; Costamante, L.; Lemoine, M.; Padovani, P.; Pueschel, E.; Resconi, E.; Tavecchio, F.; Taylor, A.; Zech, A. Progress in unveiling extreme particle acceleration in persistent astrophysical jets. Nat. Astron. 2020, 4, 124–131. [Google Scholar] [CrossRef]
  87. Ghisellini, G.; Tavecchio, F.; Chiaberge, M. Structured jets in TeV BL Lac objects and radiogalaxies. Implications for the observed properties. Astron. Astrophys. 2005, 432, 401–410. [Google Scholar] [CrossRef]
  88. Ahnen, M.L.; Ansoldi, S.; Antonelli, L.A.; Arcaro, C.; Babić, A.; Banerjee, B.; Bangale, P.; Barres de Almeida, U.; Barrio, J.A.; Becerra González, J.; et al. Extreme HBL behavior of Markarian 501 during 2012. Astron. Astrophys. 2018, 620, A181. [Google Scholar] [CrossRef]
  89. MAGIC Collaboration; Acciari, V.A.; Ansoldi, S.; Antonelli, L.A.; Arbet Engels, A.; Babić, A.; Banerjee, B.; Barres de Almeida, U.; Barrio, J.A.; Becerra González, J.; et al. An intermittent extreme BL Lac: MWL study of 1ES 2344+514 in an enhanced state. Mon. Not. R. Astron. Soc. 2020, 496, 3912–3928. [Google Scholar] [CrossRef]
  90. Abe, H.; Abe, S.; Acciari, V.A.; Agudo, I.; Aniello, T.; Ansoldi, S.; Antonelli, L.A.; Arbet Engels, A.; Arcaro, C.; Artero, M.; et al. Multi-year characterisation of the broad-band emission from the intermittent extreme BL Lac 1ES~2344+514. arXiv 2023, arXiv:2310.03922. [Google Scholar] [CrossRef]
  91. Acciari, V.A.; Agudo, I.; Aniello, T.; Ansoldi, S.; Antonelli, L.A.; Arbet Engels, A.; Artero, M.; Asano, K.; Baack, D.; Babić, A.; et al. A lower bound on intergalactic magnetic fields from time variability of 1ES 0229+200 from MAGIC and Fermi/LAT observations. Astron. Astrophys. 2023, 670, A145. [Google Scholar] [CrossRef]
  92. Aleksić, J.; Ansoldi, S.; Antonelli, L.A.; Antoranz, P.; Babic, A.; Bangale, P.; Barrio, J.A.; González, J.B.; Bednarek, W.; Bernardini, E.; et al. Black hole lightning due to particle acceleration at subhorizon scales. Science 2014, 346, 1080–1084. [Google Scholar] [CrossRef]
  93. Albert, J.; Aliu, E.; Anderhub, H.; Antoranz, P.; Armada, A.; Asensio, M.; Baixeras, C.; Barrio, J.A.; Bartko, H.; Bastieri, D.; et al. Observations of Markarian 421 with the MAGIC Telescope. Astrophys. J. 2007, 663, 125–138. [Google Scholar] [CrossRef]
  94. Aleksić, J.; Anderhub, H.; Antonelli, L.A.; Antoranz, P.; Backes, M.; Baixeras, C.; Balestra, S.; Barrio, J.A.; Bastieri, D.; Becerra González, J.; et al. MAGIC TeV gamma-ray observations of Markarian 421 during multiwavelength campaigns in 2006. Astron. Astrophys. 2010, 519, A32. [Google Scholar] [CrossRef]
  95. Abdo, A.A.; Ackermann, M.; Ajello, M.; Baldini, L.; Ballet, J.; Barbiellini, G.; Bastieri, D.; Bechtol, K.; Bellazzini, R.; Berenji, B.; et al. Fermi Large Area Telescope Observations of Markarian 421: The Missing Piece of its Spectral Energy Distribution. Astrophys. J. 2011, 736, 131. [Google Scholar] [CrossRef]
  96. Aleksić, J.; Ansoldi, S.; Antonelli, L.A.; Antoranz, P.; Babic, A.; Bangale, P.; Barres de Almeida, U.; Barrio, J.A.; Becerra González, J.; Bednarek, W.; et al. Unprecedented study of the broadband emission of Mrk 421 during flaring activity in March 2010. Astron. Astrophys. 2015, 578, A22. [Google Scholar] [CrossRef]
  97. Aleksić, J.; Ansoldi, S.; Antonelli, L.A.; Antoranz, P.; Babic, A.; Bangale, P.; Barres de Almeida, U.; Barrio, J.A.; Becerra González, J.; Bednarek, W.; et al. The 2009 multiwavelength campaign on Mrk 421: Variability and correlation studies. Astron. Astrophys. 2015, 576, A126. [Google Scholar] [CrossRef]
  98. Baloković, M.; Paneque, D.; Madejski, G.; Furniss, A.; Chiang, J.; Ajello, M.; Alexander, D.M.; Barret, D.; Blandford, R.D.; Boggs, S.E.; et al. Multiwavelength Study of Quiescent States of Mrk 421 with Unprecedented Hard X-Ray Coverage Provided by NuSTAR in 2013. Astrophys. J. 2016, 819, 156. [Google Scholar] [CrossRef]
  99. Ahnen, M.L.; Ansoldi, S.; Antonelli, L.A.; Antoranz, P.; Babic, A.; Banerjee, B.; Bangale, P.; Barres de Almeida, U.; Barrio, J.A.; Becerra González, J.; et al. Long-term multi-wavelength variability and correlation study of Markarian 421 from 2007 to 2009. Astron. Astrophys. 2016, 593, A91. [Google Scholar] [CrossRef]
  100. Abeysekara, A.U.; Archambault, S.; Archer, A.; Benbow, W.; Bird, R.; Buchovecky, M.; Buckley, J.H.; Bugaev, V.; Cardenzana, J.V.; Cerruti, M.; et al. A Search for Spectral Hysteresis and Energy-dependent Time Lags from X-Ray and TeV Gamma-Ray Observations of Mrk 421. Astrophys. J. 2017, 834, 2. [Google Scholar] [CrossRef]
  101. Acciari, V.A.; Ansoldi, S.; Antonelli, L.A.; Asano, K.; Babić, A.; Banerjee, B.; Baquero, A.; de Almeida, U.B.; Barrio, J.A.; Becerra González, J.; et al. Multiwavelength variability and correlation studies of Mrk 421 during historically low X-ray and γ-ray activity in 2015-2016. Mon. Not. R. Astron. Soc. 2021, 504, 1427–1451. [Google Scholar] [CrossRef]
  102. Acciari, V.A.; Ansoldi, S.; Antonelli, L.A.; Arbet Engels, A.; Baack, D.; Babić, A.; Banerjee, B.; Barres de Almeida, U.; Barrio, J.A.; Becerra González, J.; et al. Unraveling the Complex Behavior of Mrk 421 with Simultaneous X-Ray and VHE Observations during an Extreme Flaring Activity in 2013 April. Astrophys. J. Suppl. Ser. 2020, 248, 29. [Google Scholar] [CrossRef]
  103. Abeysekara, A.U.; Benbow, W.; Bird, R.; Brill, A.; Brose, R.; Buchovecky, M.; Buckley, J.H.; Christiansen, J.L.; Chromey, A.J.; Daniel, M.K.; et al. The Great Markarian 421 Flare of 2010 February: Multiwavelength Variability and Correlation Studies. Astrophys. J. 2020, 890, 97. [Google Scholar] [CrossRef]
  104. MAGIC Collaboration; Acciari, V.A.; Ansoldi, S.; Antonelli, L.A.; Arbet Engels, A.; Artero, M.; Asano, K.; Babić, A.; Baquero, A.; Barres de Almeida, U.; et al. Investigation of the correlation patterns and the Compton dominance variability of Mrk 421 in 2017. Astron. Astrophys. 2021, 655, A89. [Google Scholar] [CrossRef]
  105. Aleksić, J.; Antonelli, L.A.; Antoranz, P.; Backes, M.; Baixeras, C.; Barrio, J.A.; Bastieri, D.; Becerra González, J.; Bednarek, W.; Berdyugin, A.; et al. Search for an extended VHE γ-ray emission from Mrk 421 and Mrk 501 with the MAGIC Telescope. Astron. Astrophys. 2010, 524, A77. [Google Scholar] [CrossRef]
  106. Albert, J.; Aliu, E.; Anderhub, H.; Antoranz, P.; Armada, A.; Baixeras, C.; Barrio, J.A.; Bartko, H.; Bastieri, D.; Becker, J.K.; et al. Variable Very High Energy γ-Ray Emission from Markarian 501. Astrophys. J. 2007, 669, 862–883. [Google Scholar] [CrossRef]
  107. Anderhub, H.; Antonelli, L.A.; Antoranz, P.; Backes, M.; Baixeras, C.; Balestra, S.; Barrio, J.A.; Bastieri, D.; Becerra González, J.; Becker, J.K.; et al. Simultaneous Multiwavelength Observation of Mkn 501 in a Low State in 2006. Astrophys. J. 2009, 705, 1624–1631. [Google Scholar] [CrossRef]
  108. Abdo, A.A.; Ackermann, M.; Ajello, M.; Allafort, A.; Baldini, L.; Ballet, J.; Barbiellini, G.; Baring, M.G.; Bastieri, D.; Bechtol, K.; et al. Insights into the High-energy γ-ray Emission of Markarian 501 from Extensive Multifrequency Observations in the Fermi Era. Astrophys. J. 2011, 727, 129. [Google Scholar] [CrossRef]
  109. Acciari, V.A.; Arlen, T.; Aune, T.; Beilicke, M.; Benbow, W.; Böttcher, M.; Boltuch, D.; Bradbury, S.M.; Buckley, J.H.; Bugaev, V.; et al. Spectral Energy Distribution of Markarian 501: Quiescent State Versus Extreme Outburst. Astrophys. J. 2011, 729, 2. [Google Scholar] [CrossRef]
  110. Furniss, A.; Noda, K.; Boggs, S.; Chiang, J.; Christensen, F.; Craig, W.; Giommi, P.; Hailey, C.; Harisson, F.; Madejski, G.; et al. First NuSTAR Observations of Mrk 501 within a Radio to TeV Multi-Instrument Campaign. Astrophys. J. 2015, 812, 65. [Google Scholar] [CrossRef]
  111. Aleksić, J.; Ansoldi, S.; Antonelli, L.A.; Antoranz, P.; Babic, A.; Bangale, P.; Barres de Almeida, U.; Barrio, J.A.; Becerra González, J.; Bednarek, W.; et al. Multiwavelength observations of Mrk 501 in 2008. Astron. Astrophys. 2015, 573, A50. [Google Scholar] [CrossRef]
  112. Ahnen, M.L.; Ansoldi, S.; Antonelli, L.A.; Antoranz, P.; Babic, A.; Banerjee, B.; Bangale, P.; Barres de Almeida, U.; Barrio, J.A.; Becerra González, J.; et al. Multiband variability studies and novel broadband SED modeling of Mrk 501 in 2009. Astron. Astrophys. 2017, 603, A31. [Google Scholar] [CrossRef]
  113. MAGIC Collaboration; Acciari, V.A.; Ansoldi, S.; Antonelli, L.A.; Babić, A.; Banerjee, B.; Barres de Almeida, U.; Barrio, J.A.; Becerra González, J.; Bednarek, W.; et al. Study of the variable broadband emission of Markarian 501 during the most extreme Swift X-ray activity. Astron. Astrophys. 2020, 637, A86. [Google Scholar] [CrossRef]
  114. Abe, H.; Abe, S.; Acciari, V.A.; Agudo, I.; Aniello, T.; Ansoldi, S.; Antonelli, L.A.; Arbet-Engels, A.; Arcaro, C.; Artero, M.; et al. Multimessenger Characterization of Markarian 501 during Historically Low X-Ray and γ-Ray Activity. Astrophys. J. Suppl. Ser. 2023, 266, 37. [Google Scholar] [CrossRef]
  115. Aleksić, J.; Alvarez, E.A.; Antonelli, L.A.; Antoranz, P.; Asensio, M.; Backes, M.; Barrio, J.A.; Bastieri, D.; Becerra González, J.; Bednarek, W.; et al. MAGIC observations of the giant radio galaxy M 87 in a low-emission state between 2005 and 2007. Astron. Astrophys. 2012, 544, A96. [Google Scholar] [CrossRef]
  116. MAGIC Collaboration; Acciari, V.A.; Ansoldi, S.; Antonelli, L.A.; Arbet Engels, A.; Arcaro, C.; Baack, D.; Babić, A.; Banerjee, B.; Bangale, P.; et al. Monitoring of the radio galaxy M 87 during a low-emission state from 2012 to 2015 with MAGIC. Mon. Not. R. Astron. Soc. 2020, 492, 5354–5365. [Google Scholar] [CrossRef]
  117. Albert, J.; Aliu, E.; Anderhub, H.; Antonelli, L.A.; Antoranz, P.; Backes, M.; Baixeras, C.; Barrio, J.A.; Bartko, H.; Bastieri, D.; et al. Very High Energy Gamma-Ray Observations of Strong Flaring Activity in M87 in 2008 February. Astrophys. J. Lett. 2008, 685, L23. [Google Scholar] [CrossRef]
  118. Abramowski, A.; Acero, F.; Aharonian, F.; Akhperjanian, A.G.; Anton, G.; Balzer, A.; Barnacka, A.; Barres de Almeida, U.; Becherini, Y.; Becker, J.; et al. The 2010 Very High Energy γ-Ray Flare and 10 Years of Multi-wavelength Observations of M 87. Astrophys. J. 2012, 746, 151. [Google Scholar] [CrossRef]
  119. Event Horizon Telescope Collaboration; Akiyama, K.; Alberdi, A.; Alef, W.; Asada, K.; Azulay, R.; Baczko, A.K.; Ball, D.; Baloković, M.; Barrett, J.; et al. First M87 Event Horizon Telescope Results. II. Array and Instrumentation. Astrophys. J. Lett. 2019, 875, L2. [Google Scholar] [CrossRef]
  120. EHT MWL Science Working Group; Algaba, J.C.; Anczarski, J.; Asada, K.; Baloković, M.; Chandra, S.; Cui, Y.Z.; Falcone, A.D.; Giroletti, M.; Goddi, C.; et al. Broadband Multi-wavelength Properties of M87 during the 2017 Event Horizon Telescope Campaign. Astrophys. J. Lett. 2021, 911, L11. [Google Scholar] [CrossRef]
  121. Albert, J.; Aliu, E.; Anderhub, H.; Antoranz, P.; Baixeras, C.; Barrio, J.A.; Bartko, H.; Bastieri, D.; Becker, J.K.; Bednarek, W.; et al. MAGIC observations of PG 1553+113 during a multiwavelength campaign in July 2006. Astron. Astrophys. 2009, 493, 467–469. [Google Scholar] [CrossRef]
  122. Aleksić, J.; Anderhub, H.; Antonelli, L.A.; Antoranz, P.; Backes, M.; Baixeras, C.; Balestra, S.; Barrio, J.A.; Bastieri, D.; Becerra González, J.; et al. Simultaneous multi-frequency observation of the unknown redshift blazar PG 1553+113 in March-April 2008. Astron. Astrophys. 2010, 515, A76. [Google Scholar] [CrossRef]
  123. Aleksić, J.; Alvarez, E.A.; Antonelli, L.A.; Antoranz, P.; Asensio, M.; Backes, M.; Barrio, J.A.; Bastieri, D.; Becerra González, J.; Bednarek, W.; et al. PG 1553+113: Five Years of Observations with MAGIC. Astrophys. J. 2012, 748, 46. [Google Scholar] [CrossRef]
  124. Ackermann, M.; Ajello, M.; Albert, A.; Atwood, W.B.; Baldini, L.; Ballet, J.; Barbiellini, G.; Bastieri, D.; Becerra Gonzalez, J.; Bellazzini, R.; et al. Multiwavelength Evidence for Quasi-periodic Modulation in the Gamma-Ray Blazar PG 1553+113. Astrophys. J. Lett. 2015, 813, L41. [Google Scholar] [CrossRef]
  125. Covino, S.; Landoni, M.; Sandrinelli, A.; Treves, A. Looking at Blazar Light-curve Periodicities with Gaussian Processes. Astrophys. J. 2020, 895, 122. [Google Scholar] [CrossRef]
  126. Peñil, P.; Domínguez, A.; Buson, S.; Ajello, M.; Otero-Santos, J.; Barrio, J.A.; Nemmen, R.; Cutini, S.; Rani, B.; Franckowiak, A.; et al. Systematic Search for γ-Ray Periodicity in Active Galactic Nuclei Detected by the Fermi Large Area Telescope. Astrophys. J. 2020, 896, 134. [Google Scholar] [CrossRef]
  127. De Lotto, B.; Magic Collaboration. The MAGIC telescopes: Performance, results and future perspectives. J. Phys. Conf. Ser. 2012, 375, 052021. [Google Scholar] [CrossRef]
  128. Albert, J.; Aliu, E.; Anderhub, H.; Antoranz, P.; Armada, A.; Asensio, M.; Baixeras, C.; Barrio, J.A.; Bartko, H.; Bastieri, D.; et al. Observation of Very High Energy Gamma-Ray Emission from the Active Galactic Nucleus 1ES 1959+650 Using the MAGIC Telescope. Astrophys. J. 2006, 639, 761–765. [Google Scholar] [CrossRef]
  129. MAGIC Collaboration; Acciari, V.A.; Ansoldi, S.; Antonelli, L.A.; Arbet Engels, A.; Baack, D.; Babić, A.; Banerjee, B.; Barres de Almeida, U.; Barrio, J.A.; et al. Broadband characterisation of the very intense TeV flares of the blazar 1ES 1959+650 in 2016. Astron. Astrophys. 2020, 638, A14. [Google Scholar] [CrossRef]
  130. Bionta, R.M.; Blewitt, G.; Bratton, C.B.; Casper, D.; Ciocio, A.; Claus, R.; Cortez, B.; Crouch, M.; Dye, S.T.; Errede, S.; et al. Observation of a neutrino burst in coincidence with supernova 1987A in the Large Magellanic Cloud. Phys. Rev. Lett. 1987, 58, 1494–1496. [Google Scholar] [CrossRef] [PubMed]
  131. Weekes, T.C.; Badran, H.; Biller, S.D.; Bond, I.; Bradbury, S.; Buckley, J.; Carter-Lewis, D.; Catanese, M.; Criswell, S.; Cui, W.; et al. VERITAS: The Very Energetic Radiation Imaging Telescope Array System. Astropart. Phys. 2002, 17, 221–243. [Google Scholar] [CrossRef]
  132. Benbow, W. Highlights from the VERITAS AGN Observation Program. Int. Cosm. Ray Conf. 2019, 36, 632. [Google Scholar] [CrossRef]
  133. McGrath, C. VERITAS Highlights 2022. J. Phys. Conf. Ser. 2023, 2429, 012015. [Google Scholar] [CrossRef]
  134. Wagner, R. Observations of AGN at very-high energy gamma rays with the H.E.S.S. telescopes. In Proceedings of the Active Galactic Nuclei: What’s in a Name? (AGN 2016), Garching, Germany, 27 June–1 July 2016; p. 58. [Google Scholar] [CrossRef]
  135. Taylor, A.M.; Sanchez, D.; Cerruti, M. Extragalactic Observations with HESS: Past and Future. arXiv 2017, arXiv:1708.00775. [Google Scholar] [CrossRef]
  136. H. E. S. S. Collaboration; Abdalla, H.; Adam, R.; Aharonian, F.; Ait Benkhali, F.; Angüner, E.O.; Arcaro, C.; Armand, C.; Armstrong, T.; Ashkar, H.; et al. H.E.S.S. and MAGIC observations of a sudden cessation of a very-high-energy γ-ray flare in PKS 1510-089 in May 2016. Astron. Astrophys. 2021, 648, A23. [Google Scholar] [CrossRef]
  137. Anderhub, H.; Backes, M.; Biland, A.; Boller, A.; Braun, I.; Bretz, T.; Commichau, S.; Commichau, V.; Domke, M.; Dorner, D.; et al. FACT—The first Cherenkov telescope using a G-APD camera for TeV gamma-ray astronomy. Nucl. Instrum. Methods Phys. Res. 2011, 639, 58–61. [Google Scholar] [CrossRef]
  138. Dorner, D.; Arbet-Engels, A.; Baack, D.; Balbo, M.; Biland, A.; Bretz, T.; Buss, J.; Eisenberger, L.; Elsaesser, D.; Hildebrand, D.; et al. FACT—Highlights from more than Eight Years of Unbiased TeV Monitoring. In Proceedings of the 37th International Cosmic Ray Conference, Online. 18 March 2022; p. 851. [Google Scholar] [CrossRef]
  139. Romoli, C.; Chakraborty, N.; Dorner, D.; Taylor, A.M.; Blank, M. Flux Distribution of Gamma-Ray Emission in Blazars: The Example of Mrk 501. Galaxies 2018, 6, 135. [Google Scholar] [CrossRef]
  140. Abe, H.; Abe, K.; Abe, S.; Aguasca-Cabot, A.; Agudo, I.; Alvarez Crespo, N.; Antonelli, L.A.; Aramo, C.; Arbet-Engels, A.; Arcaro, C.; et al. Observations of the Crab Nebula and Pulsar with the Large-sized Telescope Prototype of the Cherenkov Telescope Array. Astrophys. J. 2023, 956, 80. [Google Scholar] [CrossRef]
  141. Actis, M.; Agnetta, G.; Aharonian, F.; Akhperjanian, A.; Aleksić, J.; Aliu, E.; Allan, D.; Allekotte, I.; Antico, F.; Antonelli, L.A.; et al. Design concepts for the Cherenkov Telescope Array CTA: An advanced facility for ground-based high-energy gamma-ray astronomy. Exp. Astron. 2011, 32, 193–316. [Google Scholar] [CrossRef]
  142. Abe, S.; Aguasca-Cabot, A.; Agudo, I.; Alvarez Crespo, N.; Antonelli, L.A.; Aramo, C.; Arbet-Engels, A.; Artero, M.; Asano, K.; Aubert, P.; et al. Multiwavelength study of the galactic PeVatron candidate LHAASO J2108+5157. Astron. Astrophys. 2023, 673, A75. [Google Scholar] [CrossRef]
  143. Abe, H.; Abe, K.; Abe, S.; Acciari, V.A.; Aguasca-Cabot, A.; Agudo, I.; Alvarez Crespo, N.; Aniello, T.; Ansoldi, S.; Antonelli, L.A.; et al. Performance of the joint LST-1 and MAGIC observations evaluated with Crab Nebula data. Astron. Astrophys. 2023, 680, A66. [Google Scholar] [CrossRef]
Figure 1. The MAGIC Florian Goebel telescopes at the Roque de Los Muchachos Observatory, Spain. Image credit: Chiara Righi.
Figure 1. The MAGIC Florian Goebel telescopes at the Roque de Los Muchachos Observatory, Spain. Image credit: Chiara Righi.
Universe 10 00080 g001
Figure 2. (a) Broadband SED of QSO B0218+357 during the 2014 flare. Green markers show historical data while the black solid line represents the assumed two-zone model. The red markers represent the MWL dataset. Reprinted with permission from [67]. (b) MWL SED of QSO B0218+357 during the monitoring campaign (red markers) compared to the 2014 flare (black markers). Reprinted with permission from [68] (Figure 11).
Figure 2. (a) Broadband SED of QSO B0218+357 during the 2014 flare. Green markers show historical data while the black solid line represents the assumed two-zone model. The red markers represent the MWL dataset. Reprinted with permission from [67]. (b) MWL SED of QSO B0218+357 during the monitoring campaign (red markers) compared to the 2014 flare (black markers). Reprinted with permission from [68] (Figure 11).
Universe 10 00080 g002
Figure 3. MWL SED of 3C 279 in 23 February 2006. The blue solid line represents the modelling under the assumption of that the external component is inside the BLR, and the red solid line when the external component is outside the BLR. Reprinted with permission from [62], Figure 8.
Figure 3. MWL SED of 3C 279 in 23 February 2006. The blue solid line represents the modelling under the assumption of that the external component is inside the BLR, and the red solid line when the external component is outside the BLR. Reprinted with permission from [62], Figure 8.
Universe 10 00080 g003
Figure 4. (a) Light curve of IC 310, observed by the MAGIC telescopes in 12–13 November 2012. Reprinted from [92], Figure S4. (b) Emission scenario for the origin of the highly variable gamma-ray emission observed in IC 310 during the 2012 activity: the black sphere represents the rotating black hole with its event horizon accreting plasma from the center of IC 310. The apple-shaped blue form represents the ergosphere surrounding the black hole. The magnetosphere is shown in red and the polar vacuum gap regions in yellow. In the gaps, the electric field of the magnetosphere has a component that runs parallel to the magnetic field accelerating particles to ultra-relativistic energies. The observed gamma rays arise from inverse-Compton scattering and copious pair production due to interactions with low-energy thermal photons from the plasma accreted by the black hole. Reprinted with permission from [92], Figure S5.
Figure 4. (a) Light curve of IC 310, observed by the MAGIC telescopes in 12–13 November 2012. Reprinted from [92], Figure S4. (b) Emission scenario for the origin of the highly variable gamma-ray emission observed in IC 310 during the 2012 activity: the black sphere represents the rotating black hole with its event horizon accreting plasma from the center of IC 310. The apple-shaped blue form represents the ergosphere surrounding the black hole. The magnetosphere is shown in red and the polar vacuum gap regions in yellow. In the gaps, the electric field of the magnetosphere has a component that runs parallel to the magnetic field accelerating particles to ultra-relativistic energies. The observed gamma rays arise from inverse-Compton scattering and copious pair production due to interactions with low-energy thermal photons from the plasma accreted by the black hole. Reprinted with permission from [92], Figure S5.
Universe 10 00080 g004
Figure 5. MM SEDs of the blazar TXS 0506+056. Panel (a) shows the MM SED for the enhanced VHE gamma-ray state, while panel (b) shows the MM SED for the lower state. Panel (c) shows a comparison with archival data, while panels (d,e) show results for different values of the maximum proton energy. Reprinted with permission from [54].
Figure 5. MM SEDs of the blazar TXS 0506+056. Panel (a) shows the MM SED for the enhanced VHE gamma-ray state, while panel (b) shows the MM SED for the lower state. Panel (c) shows a comparison with archival data, while panels (d,e) show results for different values of the maximum proton energy. Reprinted with permission from [54].
Universe 10 00080 g005
Table 2. List of the main results obtained by MAGIC on Mrk 421.
Table 2. List of the main results obtained by MAGIC on Mrk 421.
Mrk 421
Observational PeriodMain ResultsTheor. ModelRef.
November 2004–April 2005 γ -ray/X-ray corr.,one-zone SSC[93]
IC peak ∼100 GeV
22–30 April + 14 June 2006intra-night var.leptonic[94]
(29 April, ∼36 min)
5 August 2008–12 March 2010MWL SED in aone-zone SCC,[95]
quiescent stateproton-synch
characterization
March 2010 γ -ray/X-ray corr.,one-zone SCC,[96]
γ -ray and X-ray var.two-zones SSC
January–June 2009quiescent stateone-zone SCC[97]
characterization,
X-ray harder-
when-brighter,
γ -ray/X-ray corr.,
optical/X-ray
anti-corr.
January–March 2013 γ -ray/X-ray corr.,one-zone SCC,[98]
double-bumpedsuggestion of
frac. var.,multi-zone leptonic
low state
characterization
March 2007–June 2009X-ray/soft X-raysuggested SSC,[99]
corr., frac. var.or generic
increasing withhadronic
energy, differentscenarios
levels of activity
28 April–4 May 2014X-ray spectrumone-zone SSC[100]
variability
November 2014–June 2016X-ray and γ -ray [101]
harder-whensuggesting that the
brighter, double-emission is powered
bumped frac. var.,by a multiplicative
X-ray/ γ -ray/ corr.,process
VHE intra-night
var. (27 January +
12 March 2015)
11–19 April 2013intra-night var.magnetic reconnection[102]
of X-ray and VHEin a multi-zone
γ -ray bands,
VHE γ -ray/X-ray corr.scenario
February 2010limits on theone-zone SSC[103]
Doppler factorexcluded
and size of the
emission region,
time-lagged
corr. optical/VHE
December 2016–June 2017VHE/X-ray corr.,two-zone[104]
orphan γ -rayleptonic
activity,
intra-night VHE
var., UV/X-ray
anti-corr.
December 2007–February 2009upper limits onpossible constraints[105]
extended emissionon EGMF
Table 3. List of the main results obtained by MAGIC on Mrk 501.
Table 3. List of the main results obtained by MAGIC on Mrk 501.
Mrk 501
Observational PeriodMain ResultsTheor. ModelRef.
May–July 2005VHE intra-night var.,one-zone SSC[106]
spectra hardening when
increasing flux,
var. increasing
with energy
July 2006low state in VHEone-zone SSC[107]
steep VHE photon index
spectral hardening
with flux (VHE)
15 March–1 August 2009low activityone-zone SSC[108]
characterization
March 2009quiescent state [109]
characterization,one-zone SSC
X-ray peak shift of two
orders of magnitude
1 April–10 August 2013hard X-ray var. [110]
on hour timescales,one-zone SSC
five MWL SEDs
March–May 2008low stateone-zone SSC[111]
characterization,
hint of X-ray-to-VHE
correlation
15 March–1 August 2009frac. var. increasingtwo-zones SSC[112]
with energy, flaring
activity coincident
with EVPA rotation
(1 May)
March–July 2012hard X-ray and VHEone-zone SSC[88]
spectral indexes,
extreme behaviour,
VHE/X-ray corr.,
frac. var. increasing
with energy
16–31 July 2014frac. var. increasing [113]
with energy,
VHE/X-ray corr.,
narrow feature in the
VHE spectrum at 3 TeV
(19 July)
February 2017–December 2020X-ray/VHE corr.,one-zone leptonic,[114]
HE/radio corr.,two-zone leptonic,
but also hadronic
and lepto-hadronic
are considered
May and April 2008upper limits onpossible constraints[105]
extended emissionon EGMF
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Manganaro, M.; Dominis Prester, D. Highlights of the Magic Florian Goebel Telescopes in the Study of Active Galactic Nuclei. Universe 2024, 10, 80. https://doi.org/10.3390/universe10020080

AMA Style

Manganaro M, Dominis Prester D. Highlights of the Magic Florian Goebel Telescopes in the Study of Active Galactic Nuclei. Universe. 2024; 10(2):80. https://doi.org/10.3390/universe10020080

Chicago/Turabian Style

Manganaro, Marina, and Dijana Dominis Prester. 2024. "Highlights of the Magic Florian Goebel Telescopes in the Study of Active Galactic Nuclei" Universe 10, no. 2: 80. https://doi.org/10.3390/universe10020080

APA Style

Manganaro, M., & Dominis Prester, D. (2024). Highlights of the Magic Florian Goebel Telescopes in the Study of Active Galactic Nuclei. Universe, 10(2), 80. https://doi.org/10.3390/universe10020080

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop