Lucio Angelo Antonelli

We present very-high-energy gamma-ray observations of two binary black hole merger candidates, GW240615_113620 and GW241125_010116, performed with the Major Atmospheric Gamma Imaging Cherenkov (MAGIC) telescopes and the first Large-Sized Telescope of the Cherenkov Telescope Array Observatory’s (CTAO LST-1). GW240615_113620 was the best localized event of the fourth observing run of the LIGO-Virgo-KAGRA gravitational waves interferometers. GW241125_010116 was temporally and spatially coincident with a sub-threshold short-duration burst detected with the Swift-Burst Alert Telescope, the Swift-X-Ray Telescope and the Einstein Probe Follow-up X-ray Telescope. We observed the two events in stereoscopic mode, taking advantage of the improved sensitivity of joint MAGIC+LST-1 observations. No detection was achieved in the GeV–TeV gamma-ray band for any of the two sources. The unfavourable observing conditions of both events posed a challenge for a standard analysis and therefore a non standard analysis was necessary for both objects. Owing to the small localization area and the association with a GRB-like burst respectively, these events represented an unprecedented opportunity to study in details the electromagnetic emission from binary black holes merger events and, in particular, we discussed two theoretical models that predict a detectable gamma-ray emission and the possible future applications.

The ASTRI-Horn telescope is a prototype Imaging Atmospheric Cherenkov Telescope (IACT) developed within the Astrofisica con Specchi a Tecnologia Replicante Italiana (ASTRI) project. It represents the first implementation of a dual-mirror Schwarzschild–Couder (SC) optical configuration for ground-based gamma-ray IACT observations. It also served as a testbed for the Small-Sized Telescopes (SSTs) of the Cherenkov Telescope Array Observatory (CTAO) that will host at least 37 SSTs at the Paranal site in Chile, and for the forthcoming ASTRI Mini-Array. Installed in 2014 at the INAF “M.C. Fracastoro” observing station on Mount Etna, the telescope is equipped with a SiPM-based Cherenkov camera and a complete control and data acquisition system. Between 2020 and 2022, after the demonstration of the validity of the SC configuration in observing gamma-ray sources, ASTRI-Horn underwent refurbishment activities that included mirror replacement, new optical coatings, camera upgrades, and software improvements. The telescope has since been re-validated through calibration campaigns and observational tests, including the detection and spectral characterization of the Crab Nebula at TeV energies. The results confirm the correct operation of the optical system, camera electronics, and data processing chain, providing a consolidated basis for the design and development of the ASTRI Mini-Array and future CTAO-SST telescopes.

We report on the average and orbital phase-resolved polarization of Cyg X-3 in the hard state (HS) during the 2023 Imaging X-ray Polarimetry Explorer (IXPE) observational campaign. We find the polarization degree of 21.4% ± 0.4% and polarization angle of 92 . ° 2 ± 0 . ° 5 , compatible within 3 σ with the first HS IXPE observation in 2022. As the observed polarization depends on both the accretion geometry and the X-ray emission mechanism, and assuming it arises from reflection from the optically thick envelope surrounding the central source, our result indicates that both are similar in the two HS epochs on year-long timescales despite the source transitioning through an intermediate state in between. We discuss time- and energy-dependent polarization properties and their implications for the geometry and stability of the accretion funnel, finding significant orbital modulation of the polarization properties.

We report on a nine-month monitoring campaign of the blazar PG 1553+113, relying on three observations carried out in 2025 with the Imaging X-ray Polarimetry Explorer (IXPE) and supported by multiwavelength facilities. The source displayed pronounced variability across the electromagnetic spectrum, with X-ray flux changes by up to a factor of ∼5 and complex evolution of the optical polarization properties, including one of the largest (exceeding 150°) and fastest rotations in the electric vector position angle (EVPA) ever recorded. This swing of the EVPA was also accompanied by a temporary drop of the optical polarization degree to nearly zero. Significant X-ray polarization was observed during the third IXPE pointing, with a polarization degree Π _X = (18.4 ± 5.8)% and Ψ _X = 74° ± 9° in the 2–8 keV band, while only upper limits were obtained in the first two epochs. The optical data show that the second IXPE observation occurred shortly after a dramatic optical polarization event characterized by a rapid EVPA swing and strong depolarization. Two possible scenarios may explain the broadband polarimetric behavior: (i) the superposition of two emitting regions with nearly orthogonal magnetic field configurations and variable relative contributions, and (ii) the interaction of a single emitting region with a shock that temporarily reorders the magnetic field. In both cases, the data support a picture in which the X-ray and optical emissions arise from closely related but not strictly cospatial regions within a dynamically evolving, magnetically structured jet.

Extremely high-peaked BL Lac objects—also named extreme blazars —are among the most energetic and persistent extragalactic accelerators in the Universe, defined by a synchrotron emission peaking above 10 ¹⁷ Hz in X-rays. Such emission is then reprocessed and produces radiation extending deeply into very-high-energy (VHE; energy E > 100 GeV) gamma-rays. Observations in this energy band—optimally investigated by the Imaging Air-Shower Cherenkov telescopes—are crucial for probing the physical processes that drive their extreme behavior. This study extends our investigation of extreme blazars in the VHE gamma-ray range, providing a second new mini-catalog of sources observed by the MAGIC telescopes. We report on the monitoring of seven targets between 2017 and 2025, including four newly observed sources and three that have been part of long-term observation campaigns, for a total of approximately 338 hr of observations. The analysis of MAGIC data reveals two new VHE detections of extreme blazars, along with three additional sources showing hints of VHE emission. Joint observations of MAGIC and the first Large-Sized Telescope (LST-1) also confirmed a new VHE extreme blazar. Our results are complemented by simultaneous multiwavelength observations in other energy bands, including optical-UV, X-rays, and high-energy gamma-rays (100 MeV < E < 100 GeV). We confirm typical behavior of extreme blazars, such as a modest variability and a “harder-when-brighter” trend in X-rays across the sample. This new set increases the population of extreme blazars, and their broadband analysis confirms the physical properties of these extreme sources.

Dwarf spheroidal galaxies (dSphs) are excellent targets for indirect dark matter (DM) searches using gamma-ray telescopes because they are thought to have high DM content and a low astrophysical background. The sensitivity of these searches is improved by combining the observations of dSphs made by different gamma-ray telescopes. We present the results of a combined search by the most sensitive currently operating gamma-ray telescopes, namely: the satellite-borne Fermi -LAT telescope; the ground-based imaging atmospheric Cherenkov telescope arrays H.E.S.S., MAGIC, and VERITAS; and the HAWC water Cherenkov detector. Individual datasets were analyzed using a common statistical approach. Results were subsequently combined via a global joint likelihood analysis. We obtain constraints on the velocity-weighted cross section 〈 σv 〉 for DM self-annihilation as a function of the DM particle mass. This five-instrument combination allows the derivation of up to 2-3 times more constraining upper limits on 〈 σv 〉 than the individual results over a wide mass range spanning from 5 GeV to 100 TeV. Depending on the DM content modeling, the 95% confidence level observed limits reach 1.5×10 ^-24 cm ³ s ^-1 and 3.2×10 ^-25 cm ³ s ^-1 , respectively, in the τ ⁺ τ ^- annihilation channel for a DM mass of 2 TeV.

RCW 86 is a supernova remnant whose origin has recently been linked to an off-center explosion within a cavity created by its progenitor star. In the southwestern region, the forward shock is thought to have reached the cavity wall, encountering diverse environmental conditions. We report on the spatially resolved X-ray polarimetric observation of RCW 86 with the Imaging X-ray Polarimetry Explorer (IXPE). In the 2–4.5 keV energy band we find no significant detection of polarization. Employing a dedicated background subtraction procedure and Bayesian spectropolarimetric fitting, we derive 99% upper limits on the polarization degree of the synchrotron component: 15% in higher-statistics regions and 30%–40% in lower-statistics regions. These upper limits on the polarization degree in several regions exclude the possibility of a strongly coherent magnetic field down to the subparsec scale, and that of a moderately coherent one on the scale of the synchrotron features as resolved by IXPE. The results indicate that the shocks in the southwestern rim of RCW 86 propagate more slowly than the unshocked ejecta at their locations, yet exceed the measured proper motion speeds. This behavior is consistent with reflected shocks occurring in tenuous regions of the shocked ejecta, distinct from regions that are radio-bright.

In 2013 April, the TeV blazar Markarian 421 underwent one of its most powerful emission outbursts recorded to date. An extensive multi-instrument campaign featuring MAGIC, VERITAS, and NuSTAR provided comprehensive very high-energy (VHE; E > 100 GeV) and X-ray coverage over nine consecutive days. The VHE flux peaked at approximately 15 times that of the Crab Nebula, with rapid variability detected on timescales down to 15 minutes in both X-ray and VHE bands. This rich data set, characterized by its dense temporal coverage and high photon statistics, offers an unparalleled opportunity to probe the broadband emission dynamics in blazars. In this work, we perform a detailed spectral analysis of the X-ray and VHE emissions on subhour timescales throughout the flare. We identify several clockwise spectral hysteresis loops in the X-rays, revealing a spectral evolution more complex than a simple harder-when-brighter trend. The VHE spectrum extends beyond 10 TeV, and its temporal evolution closely mirrors the behavior in the X-rays. Crucially, we report the first evidence of VHE spectral hysteresis occurring simultaneously with the X-ray loops. To interpret these findings, we apply a time-dependent leptonic model to 240 broadband spectral energy distributions (SEDs) binned on a 15 minute scale, allowing us to self-consistently track the particle distribution’s history. Our modeling shows that the majority of the subhour flux and spectral variations are driven by changes in the luminosity and slope of the injected electron distribution. The required variations in the electron slope are difficult to reconcile with magnetic reconnection, but they are consistent with a shock-acceleration scenario where the shock compression ratio evolves by a factor of ∼2. The model also points to a relatively stable magnetic field and emitting region size, favoring a scenario where the emission originates from a stationary feature in the jet, such as a recollimation shock. However, this scenario requires a jet Lorentz factor that significantly exceeds values from VLBI measurements to account for the high minimum electron energy implied by the lack of variability in the optical band.

The search for sources of high-energy astrophysical neutrinos can be significantly advanced through a multimessenger approach, which seeks to detect the γ -rays that accompany neutrinos as they are produced at their sources. Multimessenger observations have so far provided the first evidence for a neutrino source, illustrated by the joint detection of the flaring blazar TXS 0506+056 in high-energy ( E > 1 GeV) and very-high-energy (VHE; E > 100 GeV) γ -rays in coincidence with the high-energy neutrino IceCube-170922A, identified by IceCube. Imaging atmospheric Cherenkov telescopes (IACTs), namely FACT, H.E.S.S., MAGIC, and VERITAS, continue to conduct extensive neutrino target-of-opportunity follow-up programs. These programs have two components: follow-up observations of single astrophysical neutrino candidate events (such as IceCube-170922A), and observation of known γ -ray sources after the identification of a cluster of neutrino events by IceCube. Here we present a comprehensive analysis of follow-up observations of high-energy neutrino events observed by the four IACTs between 2017 September (after the IceCube-170922A event) and 2021 January. Our study found no associations between γ -ray sources and the observed neutrino events. We provide a detailed overview of each neutrino event and its potential counterparts. Furthermore, a joint analysis of all IACT data is included, yielding combined upper limits on the VHE γ -ray flux.

Seyfert galaxies are emerging as a promising source class of high-energy neutrinos. The Seyfert galaxies NGC 4151 and NGC 1068 have respectively come up as the most promising counterparts of a 3σ and of a 4.2σ neutrino excesses detected by IceCube in the TeV energy range. Constraining the very high energy (VHE) emission associated with the neutrino signal is crucial to unveiling the mechanism and site of neutrino production. In this work, we present the first results of the VHE observations (∼29 hours) of NGC 4151 with the MAGIC telescopes. We detected no gamma-ray excess in the direction of NGC 4151, and we derived constraining upper limits on the VHE gamma-ray flux. The integral flux upper limit (at the 95% confidence level) above 200 GeV is f = 2.3 × 10⁻¹² cm⁻² s⁻¹. Comparison of the MAGIC and IceCube measurements suggests the presence of a gamma-ray obscured accelerator, and it allowed us to constrain the gamma-ray optical depth and the size of the neutrino production site.