Long-term multi-wavelength study of 1ES 0647+250

The BL Lac object 1ES 0647+250 is one of the few distant $\gamma$-ray emitting blazars detected at very high energies (VHE, $\gtrsim$100 GeV) during a non-flaring state. It was detected with the MAGIC telescopes during its low activity in the years 2009-2011, as well as during three flaring activities in the years 2014, 2019 and 2020, with the highest VHE flux in the latter epoch. An extensive multi-instrument data set was collected within several coordinated observing campaigns throughout these years. We aim to characterise the long-term multi-band flux variability of 1ES 0647+250, as well as its broadband spectral energy distribution (SED) during four distinct activity states selected in four different epochs, in order to constrain the physical parameters of the blazar emission region under certain assumptions. We evaluate the variability and correlation of the emission in the different energy bands with the fractional variability and the Z-transformed Discrete Correlation Function, as well as its spectral evolution in X-rays and $\gamma$ rays. Owing to the controversy in the redshift measurements of 1ES 0647+250 reported in the literature, we also estimate its distance in an indirect manner through the comparison of the GeV and TeV spectra from simultaneous observations with Fermi-LAT and MAGIC during the strongest flaring activity detected to date. Moreover, we interpret the SEDs from the four distinct activity states within the framework of one-component and two-component leptonic models, proposing specific scenarios that are able to reproduce the available multi-instrument data.Comment: 20 pages, 7 figures. Accepted in A&A. Corresponding authors: Jorge Otero-Santos; Daniel Morcuende; Vandad Fallah Ramazani; Daniela Dorner; David Paneque (mailto: [email protected]

A detailed study of the very-high-energy Crab pulsar emission with the LST-1

Context: There are currently three pulsars firmly detected by imaging atmospheric Cherenkov telescopes (IACTs), two of them reaching TeV energies, challenging models of very-high-energy (VHE) emission in pulsars. More precise observations are needed to better characterize pulsar emission at these energies. The LST-1 is the prototype of the Large-Sized Telescope, that will be part of the Cherenkov Telescope Array Observatory (CTAO). Its improved performance over previous IACTs makes it well suited for studying pulsars. Aims: To study the Crab pulsar emission with the LST-1, improving and complementing the results from other telescopes. These observations can also be used to characterize the potential of the LST-1 to study other pulsars and detect new ones. Methods: We analyzed a total of $\sim$103 hours of gamma-ray observations of the Crab pulsar conducted with the LST-1 in the period from September 2020 to January 2023. The observations were carried out at zenith angles less than 50 degrees. A new analysis of the Fermi-LAT data was also performed, including $\sim$14 years of observations. Results: The Crab pulsar phaseogram, long-term light-curve, and phase-resolved spectra are reconstructed with the LST-1 from 20 GeV to 450 GeV for P1 and up to 700 GeV for P2. The pulsed emission is detected with a significance of 15.2$σ$. The two characteristic emission peaks of the Crab pulsar are clearly detected (>10$σ$), as well as the so-called bridge emission (5.7$σ$). We find that both peaks are well described by power laws, with spectral indices of $\sim$3.44 and $\sim$3.03 respectively. The joint analysis of Fermi-LAT and LST-1 data shows a good agreement between both instruments in the overlapping energy range. The detailed results obtained in the first observations of the Crab pulsar with LST-1 show the potential that CTAO will have to study this type of sources

Observations of the Crab Nebula and Pulsar with the Large-sized Telescope Prototype of the Cherenkov Telescope Array

The Cherenkov Telescope Array (CTA) is a next-generation ground-based observatory for gamma-ray astronomy at very high energies. The Large-Sized Telescope prototype (LST-1) is located at the CTA-North site, on the Canary Island of La Palma. LSTs are designed to provide optimal performance in the lowest part of the energy range covered by CTA, down to ≃20 GeV. LST-1 started performing astronomical observations in 2019 November, during its commissioning phase, and it has been taking data ever since. We present the first LST-1 observations of the Crab Nebula, the standard candle of very-high-energy gamma-ray astronomy, and use them, together with simulations, to assess the performance of the telescope. LST-1 has reached the expected performance during its commissioning period—only a minor adjustment of the preexisting simulations was needed to match the telescope’s behavior. The energy threshold at trigger level is around 20 GeV, rising to ≃30 GeV after data analysis. Performance parameters depend strongly on energy, and on the strength of the gamma-ray selection cuts in the analysis: angular resolution ranges from 0.°12-0.°40, and energy resolution from 15%-50%. Flux sensitivity is around 1.1% of the Crab Nebula flux above 250 GeV for a 50 hr observation (12% for 30 minutes). The spectral energy distribution (in the 0.03-30 TeV range) and the light curve obtained for the Crab Nebula agree with previous measurements, considering statistical and systematic uncertainties. A clear periodic signal is also detected from the pulsar at the center of the Nebula

Multiwavelength study of the galactic PeVatron candidate LHAASO J2108+5157

Context. Several new ultrahigh-energy (UHE) γ-ray sources have recently been discovered by the Large High Altitude Air Shower Observatory (LHAASO) collaboration. These represent a step forward in the search for the so-called Galactic PeVatrons, the enigmatic sources of the Galactic cosmic rays up to PeV energies. However, it has been shown that multi-TeV γ-ray emission does not necessarily prove the existence of a hadronic accelerator in the source; indeed this emission could also be explained as inverse Compton scattering from electrons in a radiation-dominated environment. A clear distinction between the two major emission mechanisms would only be made possible by taking into account multi-wavelength data and detailed morphology of the source. Aims. We aim to understand the nature of the unidentified source LHAASO J2108+5157, which is one of the few known UHE sources with no very high-energy (VHE) counterpart. Methods. We observed LHAASO J2108+5157 in the X-ray band with XMM-Newton in 2021 for a total of 3.8 hours and at TeV energies with the Large-Sized Telescope prototype (LST-1), yielding 49 hours of good-quality data. In addition, we analyzed 12 years of Fermi-LAT data, to better constrain emission of its high-energy (HE) counterpart 4FGL J2108.0+5155. We used naima and jetset software packages to examine the leptonic and hadronic scenario of the multi-wavelength emission of the source. Results. We found an excess (3.7σ) in the LST-1 data at energies E > 3 TeV. Further analysis of the whole LST-1 energy range, assuming a point-like source, resulted in a hint (2.2σ) of hard emission, which can be described with a single power law with a photon index of Σ = 1.6 ± 0.2 the range of 0.3 - 100 TeV. We did not find any significant extended emission that could be related to a supernova remnant (SNR) or pulsar wind nebula (PWN) in the XMM-Newton data, which puts strong constraints on possible synchrotron emission of relativistic electrons. We revealed a new potential hard source in Fermi-LAT data with a significance of 4σ and a photon index of Σ = 1.9 ± 0.2, which is not spatially correlated with LHAASO J2108+5157, but including it in the source model we were able to improve spectral representation of the HE counterpart 4FGL J2108.0+5155. Conclusions. The LST-1 and LHAASO observations can be explained as inverse Compton-dominated leptonic emission of relativistic electrons with a cutoff energy of 100-30+70 TeV. The low magnetic field in the source imposed by the X-ray upper limits on synchrotron emission is compatible with a hypothesis of a PWN or a TeV halo. Furthermore, the spectral properties of the HE counterpart are consistent with a Geminga-like pulsar, which would be able to power the VHE-UHE emission. Nevertheless, the lack of a pulsar in the neighborhood of the UHE source is a challenge to the PWN/TeV-halo scenario. The UHE γ rays can also be explained as π0 decay-dominated hadronic emission due to interaction of relativistic protons with one of the two known molecular clouds in the direction of the source. Indeed, the hard spectrum in the LST-1 band is compatible with protons escaping a shock around a middle-aged SNR because of their high low-energy cut-off, but the origin of the HE γ-ray emission remains an open question

Constraints on axion-like particles with the Perseus Galaxy Cluster with MAGIC

Axion-like particles (ALPs) are pseudo-Nambu-Goldstone bosons that emerge in various theories beyond the standard model. These particles can interact with high-energy photons in external magnetic fields, influencing the observed gamma-ray spectrum. This study analyzes 41.3 hrs of observational data from the Perseus Galaxy Cluster collected with the MAGIC telescopes. We focused on the spectra the radio galaxy in the center of the cluster: NGC 1275. By modeling the magnetic field surrounding this target, we searched for spectral indications of ALP presence. Despite finding no statistical evidence of ALP signatures, we were able to exclude ALP models in the sub-micro electronvolt range. Our analysis improved upon previous work by calculating the full likelihood and statistical coverage for all considered models across the parameter space. Consequently, we achieved the most stringent limits to date for ALP masses around 50 neV, with cross sections down to $g_{a\gamma} = 3 \times 10^{-12}$ GeV$^{-1}$.Comment: 25 pages, 10 figures, accepted for publication in Physics of the Dark Univers

MAGIC observations provide compelling evidence of hadronic multi-TeV emission from the putative PeVatron SNR G106.3+2.7

Context. Certain types of supernova remnants (SNRs) in our Galaxy are assumed to be PeVatrons, capable of accelerating cosmic rays (CRs) to ∼ PeV energies. However, conclusive observational evidence for this has not yet been found. The SNR G106.3+2.7, detected at 1- 100 TeV energies by different γ-ray facilities, is one of the most promising PeVatron candidates. This SNR has a cometary shape, which can be divided into a head and a tail region with different physical conditions. However, in which region the 100 TeV emission is produced has not yet been identified because of the limited position accuracy and/or angular resolution of existing observational data. Additionally, it remains unclear as to whether the origin of the γ-ray emission is leptonic or hadronic. Aims. With the better angular resolution provided by new MAGIC data compared to earlier γ-ray datasets, we aim to reveal the acceleration site of PeV particles and the emission mechanism by resolving the SNR G106.3+2.7 with 0.1 resolution at TeV energies. Methods. We observed the SNR G106.3+2.7 using the MAGIC telescopes for 121.7 h in total - after quality cuts - between May 2017 and August 2019. The analysis energy threshold is ∼0.2 TeV, and the angular resolution is 0.07-0.1. We examined the γ-ray spectra of different parts of the emission, whilst benefitting from the unprecedented statistics and angular resolution at these energies provided by our new data. We also used measurements at other wavelengths such as radio, X-rays, GeV γ-rays, and 10 TeV γ-rays to model the emission mechanism precisely. Results. We detect extended γ-ray emission spatially coincident with the radio continuum emission at the head and tail of SNR G106.3+2.7. The fact that we detect a significant γ-ray emission with energies above 6.0 TeV from only the tail region suggests that the emissions above 10 TeV detected with air shower experiments (Milagro, HAWC, Tibet ASγ and LHAASO) are emitted only from the SNR tail. Under this assumption, the multi-wavelength spectrum of the head region can be explained with either hadronic or leptonic models, while the leptonic model for the tail region is in contradiction with the emission above 10 TeV and X-rays. In contrast, the hadronic model could reproduce the observed spectrum at the tail by assuming a proton spectrum with a cutoff energy of ∼1 PeV for that region. Such high-energy emission in this middle-aged SNR (4-10 kyr) can be explained by considering a scenario where protons escaping from the SNR in the past interact with surrounding dense gases at present. Conclusions. The γ-ray emission region detected with the MAGIC telescopes in the SNR G106.3+2.7 is extended and spatially coincident with the radio continuum morphology. The multi-wavelength spectrum of the emission from the tail region suggests proton acceleration up to ∼PeV, while the emission mechanism of the head region could either be hadronic or leptonic

MAGIC detection of GRB 201216C at z = 1.1

Gamma-ray bursts (GRBs) are explosive transient events occurring at cosmological distances, releasing a large amount of energy as electromagnetic radiation over several energy bands. We report the detection of the long GRB 201216C by the MAGIC telescopes. The source is located at z = 1.1 and thus it is the farthest one detected at very high energies. The emission above 70 GeV of GRB 201216C is modelled together with multiwavelength data within a synchrotron and synchrotron self-Compton (SSC) scenario. We find that SSC can explain the broad-band data well from the optical to the very-high-energy band. For the late-time radio data, a different component is needed to account for the observed emission. Differently from previous GRBs detected in the very-high-energy range, the model for GRB 201216C strongly favours a wind-like medium. The model parameters have values similar to those found in past studies of the afterglows of GRBs detected up to GeV energies

MAGIC observations provide compelling evidence of the hadronic multi-TeV emission from the putative PeVatron SNR G106.3+2.7

The SNR G106.3+2.7, detected at 1--100 TeV energies by different $\gamma$-ray facilities, is one of the most promising PeVatron candidates. This SNR has a cometary shape which can be divided into a head and a tail region with different physical conditions. However, it is not identified in which region the 100 TeV emission is produced due to the limited position accuracy and/or angular resolution of existing observational data. Additionally, it remains unclear whether the origin of the $\gamma$-ray emission is leptonic or hadronic. With the better angular resolution provided by these new MAGIC data compared to earlier $\gamma$-ray datasets, we aim to reveal the acceleration site of PeV particles and the emission mechanism by resolving the SNR G106.3+2.7 with 0.1$^\circ$ resolution at TeV energies. We detected extended $\gamma$-ray emission spatially coincident with the radio continuum emission at the head and tail of SNR G106.3+2.7. The fact that we detected a significant $\gamma$-ray emission with energies above 6.0 TeV from the tail region only suggests that the emissions above 10 TeV, detected with air shower experiments (Milagro, HAWC, Tibet AS$\gamma$ and LHAASO), are emitted only from the SNR tail. Under this assumption, the multi-wavelength spectrum of the head region can be explained with either hadronic or leptonic models, while the leptonic model for the tail region is in contradiction with the emission above 10 TeV and X-rays. In contrast, the hadronic model could reproduce the observed spectrum at the tail by assuming a proton spectrum with a cutoff energy of $\sim 1$ PeV for the tail region. Such a high energy emission in this middle-aged SNR (4--10 kyr) can be explained by considering the scenario that protons escaping from the SNR in the past interact with surrounding dense gases at present.Comment: 13 pages, 7 figures, Accepted for publication in A&

A lower bound on intergalactic magnetic fields from time variability of 1ES 0229+200 from MAGIC and Fermi/LAT observations

Extended and delayed emission around distant TeV sources induced by the effects of propagation of gamma rays through the intergalactic medium can be used for the measurement of the intergalactic magnetic field (IGMF). We search for delayed GeV emission from the hard-spectrum TeV blazar 1ES 0229+200 with the goal to detect or constrain the IGMF-dependent secondary flux generated during the propagation of TeV gamma rays through the intergalactic medium. We analyze the most recent MAGIC observations over a 5 year time span and complement them with historic data of the H.E.S.S. and VERITAS telescopes along with a 12-year long exposure of the Fermi/LAT telescope. We use them to trace source evolution in the GeV-TeV band over one-and-a-half decade in time. We use Monte Carlo simulations to predict the delayed secondary gamma-ray flux, modulated by the source variability, as revealed by TeV-band observations. We then compare these predictions for various assumed IGMF strengths to all available measurements of the gamma-ray flux evolution. We find that the source flux in the energy range above 200 GeV experiences variations around its average on the 14 years time span of observations. No evidence for the flux variability is found in 1-100 GeV energy range accessible to Fermi/LAT. Non-detection of variability due to delayed emission from electromagnetic cascade developing in the intergalactic medium imposes a lower bound of B>1.8e-17 G for long correlation length IGMF and B>1e-14 G for an IGMF of the cosmological origin. Though weaker than the one previously derived from the analysis of Fermi/LAT data, this bound is more robust, being based on a conservative intrinsic source spectrum estimate and accounting for the details of source variability in the TeV energy band. We discuss implications of this bound for cosmological magnetic fields which might explain the baryon asymmetry of the Universe.Comment: 10 pages, 5 figures, accepted to A&A. Corresponding authors: Ievgen Vovk, Paolo Da Vela (mailto:[email protected]) and Andrii Neronov (mailto:[email protected]

Multi-year characterisation of the broad-band emission from the intermittent extreme BL Lac 1ES~2344+514

The BL Lac 1ES 2344+514 is known for temporary extreme properties (e.g., a shift of the synchrotron SED peak energy $\nu_{synch,p}$ above 1keV). While those extreme states were so far observed only during high flux levels, additional multi-year observing campaigns are required to achieve a coherent picture. Here, we report the longest investigation of the source from radio to VHE performed so far, focusing on a systematic characterisation of the intermittent extreme states. While our results confirm that 1ES 2344+514 typically exhibits $\nu_{synch,p}>$1keV during elevated flux periods, we also find periods where the extreme state coincides with low flux activity. A strong spectral variability thus happens in the quiescent state, and is likely caused by an increase of the electron acceleration efficiency without a change in the electron injection luminosity. We also report a strong X-ray flare (among the brightest for 1ES 2344+514) without a significant shift of $\nu_{synch,p}$. During this particular flare, the X-ray spectrum is among the softest of the campaign. It unveils complexity in the spectral evolution, where the common harder-when-brighter trend observed in BL Lacs is violated. During a low and hard X-ray state, we find an excess of the UV flux with respect to an extrapolation of the X-ray spectrum to lower energies. This UV excess implies that at least two regions contribute significantly to the infrared/optical/ultraviolet/X-ray emission. Using the simultaneous MAGIC, XMM-Newton, NuSTAR, and AstroSat observations, we argue that a region possibly associated with the 10 GHz radio core may explain such an excess. Finally, we investigate a VHE flare, showing an absence of simultaneous variability in the 0.3-2keV band. Using a time-dependent leptonic modelling, we show that this behaviour, in contradiction to single-zone scenarios, can instead be explained by a two-component model.Comment: Accepted for publication in Astronomy & Astrophysic