An intermittent extreme BL Lac: MWL study of 1ES 2344+514 in an enhanced state

Extreme high-frequency BL Lacs (EHBL) feature their synchrotron peak of the broad-band spectral energy distribution (SED) at vs ≥ 1017 Hz. The BL Lac object 1ES 2344+514 was included in the EHBL family because of its impressive shift of the synchrotron peak in 1996. During the following years, the source appeared to be in a low state without showing any extreme behaviours. In 2016 August, 1ES 2344+514 was detected with the groundbased γ -ray telescope FACT during a high γ -ray state, triggering multiwavelength (MWL) observations. We studied the MWL light curves of 1ES 2344+514 during the 2016 flaring state, using data from radio to very-high-energy (VHE) γ -rays taken with OVRO, KAIT, KVA, NOT, some telescopes of the GASP-WEBT collaboration at the Teide, Crimean, and St. Petersburg observatories, Swift-UVOT, Swift-XRT, Fermi-LAT, FACT, and MAGIC. With simultaneous observations of the flare, we built the broad-band SED and studied it in the framework of a leptonic and a hadronic model. The VHE γ -ray observations show a flux level of 55 per cent of the Crab Nebula flux above 300 GeV, similar to the historical maximum of 1995. The combination of MAGIC and Fermi-LAT spectra provides an unprecedented characterization of the inverse-Compton peak for this object during a flaring episode. The _ index of the intrinsic spectrum in the VHE γ -ray band is 2.04 ± 0.12stat ± 0.15sys.We find the source in an extreme state with a shift of the position of the synchrotron peak to frequencies above or equal to 1018 Hz.The financial support of the German BMBF and MPG, the Italian INFN and INAF, the Swiss National Fund SNF, the ERDF under the Spanish MINECO (FPA2017-87859-P, FPA2017-85668-P, FPA2017-82729-C6-2-R, FPA2017-82729-C6-6-R, FPA2017-82729-C6-5-R, AYA2015-71042-P, AYA2016-76012-C3-1-P, ESP2017-87055-C2-2-P, and FPA2017-90566-REDC), the Indian Department of Atomic Energy, the Japanese JSPS and MEXT, the Bulgarian Ministry of Education and Science, National RI Roadmap Project DO1-153/28.08.2018, and the Academy of Finland grant no. 320045 is gratefully acknowledged. This work was also supported by the Spanish Centro de Excelencia ‘Severo Ochoa’ SEV-2016-0588 and SEV-2015-0548, and Unidad de Excelencia ‘María de Maeztu’ MDM-2014-0369, by the Croatian Science Foundation (HrZZ) Project IP-2016-06-9782 and the University of Rijeka Project 13.12.1.3.02, by the DFG Collaborative Research Centers SFB823/C4 and SFB876/C3, the Polish National Research Centre grant UMO-2016/22/M/ST9/00382 and by the Brazilian MCTIC, CNPq, and FAPERJ. The FACT collaboration acknowledges the important contributions from ETH Zurich grants ETH-10.08-2 and ETH-27.12-1 as well as the funding by the Swiss SNF and the German BMBF (Verbundforschung Astro- und Astroteilchenphysik) and HAP (Helmoltz Alliance for Astroparticle Physics) are gratefully acknowledged. Part of this work is supported by Deutsche Forschungsgemeinschaft (DFG) within the Collaborative Research Center SFB 876 ‘Providing Information by Resource-Constrained Analysis’, project C3

Measurement of the Extragalactic Background Light using MAGIC and Fermi-LAT gamma-ray observations of blazars up to z = 1

We present a measurement of the extragalactic background light (EBL) based on a joint likelihood analysis of 32 gamma-ray spectra for 12 blazars in the redshift range z = 0.03 to 0.944, obtained by the MAGIC telescopes and Fermi-LAT. The EBL is the part of the diffuse extragalactic radiation spanning the ultraviolet, visible and infrared bands. Major contributors to the EBL are the light emitted by stars through the history of the universe, and the fraction of it which was absorbed by dust in galaxies and re-emitted at longer wavelengths. The EBL can be studied indirectly through its effect on very-high energy photons that are emitted by cosmic sources and absorbed via photon-photon interactions during their propagation across cosmological distances. We obtain estimates of the EBL density in good agreement with state-of-the-art models of the EBL production and evolution. The 1-sigma upper bounds, including systematic uncertainties, are between 13% and 23% above the nominal EBL density in the models. No anomaly in the expected transparency of the universe to gamma rays is observed in any range of optical depth.We also perform a wavelength-resolved EBL determination, which results in a hint of an excess of EBL in the 0.18 - 0.62 $\mu$m range relative to the studied models, yet compatible with them within systematics.Comment: Accepted by MNRA

Study of the GeV to TeV morphology of the γ Cygni SNR (G 78.2+2.1) with MAGIC and Fermi-LAT: Evidence for cosmic ray escape

Context. Diffusive shock acceleration (DSA) is the most promising mechanism that accelerates Galactic cosmic rays (CRs) in the shocks of supernova remnants (SNRs). It is based on particles scattering caused by turbulence ahead and behind the shock. The turbulence upstream is supposedly generated by the CRs, but this process is not well understood. The dominant mechanism may depend on the evolutionary state of the shock and can be studied via the CRs escaping upstream into the interstellar medium (ISM). Aims. Previous observations of the γ Cygni SNR showed a difference in morphology between GeV and TeV energies. Since this SNR has the right age and is at the evolutionary stage for a significant fraction of CRs to escape, our aim is to understand γ-ray emission in the vicinity of the γ Cygni SNR. Methods. We observed the region of the γ Cygni SNR with the MAGIC Imaging Atmospheric Cherenkov telescopes between 2015 May and 2017 September recording 87 h of good-quality data. Additionally, we analysed Fermi-LAT data to study the energy dependence of the morphology as well as the energy spectrum in the GeV to TeV range. The energy spectra and morphology were compared against theoretical predictions, which include a detailed derivation of the CR escape process and their γ-ray generation. Results. The MAGIC and Fermi-LAT data allowed us to identify three emission regions that can be associated with the SNR and that dominate at different energies. Our hadronic emission model accounts well for the morphology and energy spectrum of all source components. It constrains the time-dependence of the maximum energy of the CRs at the shock, the time-dependence of the level of turbulence, and the diffusion coefficient immediately outside the SNR shock. While in agreement with the standard picture of DSA, the time-dependence of the maximum energy was found to be steeper than predicted, and the level of turbulence was found to change over the lifetime of the SNR. © 2023 EDP Sciences. All rights reserved

First detection of VHE gamma-ray emission from TXS 1515-273, study of its X-ray variability and spectral energy distribution

We report here on the first multi-wavelength (MWL) campaign on the blazar TXS 1515-273, undertaken in 2019 and extending from radio to very-high-energy gamma rays (VHE). Up until now, this blazar had not been the subject of any detailed MWL observations. It has a rather hard photon index at GeV energies and was considered a candidate extreme high-synchrotronpeaked source. MAGIC observations resulted in the first-time detection of the source in VHE with a statistical significance of 7.6$\sigma$. The average integral VHE flux of the source is 6 $\pm$ 1% of the Crab nebula flux above 400 GeV. X-ray coverage was provided by Swift-XRT, XMMNewton, and NuSTAR. The long continuous X-ray observations were separated by $\sim$ 9 h, both showing clear hour scale flares. In the XMM-Newton data, both the rise and decay timescales are longer in the soft X-ray than in the hard X-ray band, indicating the presence of a particle cooling regime. The X-ray variability timescales were used to constrain the size of the emission region and the strength of the magnetic field. The data allowed us to determine the synchrotron peak frequency and classify the source as a flaring high, but not extreme, synchrotron peaked object. Considering the constraints and variability patterns from the X-ray data, we model the broad-band spectral energy distribution. We applied a simple one-zone model, which could not reproduce the radio emission and the shape of the optical emission, and a two-component leptonic model with two interacting components, enabling us to reproduce the emission from radio to VHE band

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]

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]

Investigating the blazar TXS 0506+056 through sharp multi-wavelength eyes during 2017-2019

The blazar TXS 0506+056 got into the spotlight of the astrophysical community in September 2017, when a high-energy neutrino detected by IceCube (IceCube-170922A) was associated at the 3 $\sigma$ level to a $\gamma$-ray flare from this source. This multi-messenger photon-neutrino association remains, as per today, the most significant one ever observed. TXS 0506+056 was a poorly studied object before the IceCube-170922A event. To better characterize its broad-band emission, we organized a multi-wavelength campaign lasting 16 months (November 2017 to February 2019), covering the radio-band (Mets\"ahovi, OVRO), the optical/UV (ASAS-SN, KVA, REM, Swift/UVOT), the X-rays (Swift/XRT, NuSTAR), the high-energy $\gamma$ rays (Fermi/LAT) and the very-high-energy (VHE) $\gamma$ rays (MAGIC). In $\gamma$ rays, the behaviour of the source was significantly different from the 2017 one: MAGIC observations show the presence of flaring activity during December 2018, while the source only shows an excess at the 4$\sigma$ level during the rest of the campaign (74 hours of accumulated exposure); Fermi/LAT observations show several short (days-to-week timescale) flares, different from the long-term brightening of 2017. No significant flares are detected at lower energies. The radio light curve shows an increasing flux trend, not seen in other wavelengths. We model the multi-wavelength spectral energy distributions in a lepto-hadronic scenario, in which the hadronic emission emerges as Bethe-Heitler and pion-decay cascade in the X-rays and VHE $\gamma$ rays. According to the model presented here, the December 2018 $\gamma$-ray flare was connected to a neutrino emission that was too brief and not bright enough to be detected by current neutrino instruments.Comment: 18 pages, 6 figures; in press in Ap

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&

Multi-messenger characterization of Mrk 501 during historically low X-ray and γ\gamma-ray activity

We study the broadband emission of Mrk 501 using multi-wavelength observations from 2017 to 2020 performed with a multitude of instruments, involving, among others, MAGIC, Fermi-LAT, NuSTAR, Swift, GASP-WEBT, and OVRO. Mrk 501 showed an extremely low broadband activity, which may help to unravel its baseline emission. Nonetheless, significant flux variations are detected at all wavebands, with the highest occurring at X-rays and very-high-energy (VHE) $\gamma$-rays. A significant correlation ($>$3$\sigma$) between X-rays and VHE $\gamma$-rays is measured, supporting leptonic scenarios to explain the variable parts of the emission, also during low activity. This is further supported when we extend our data from 2008 to 2020, and identify, for the first time, significant correlations between Swift-XRT and Fermi-LAT. We additionally find correlations between high-energy $\gamma$-rays and radio, with the radio lagging by more than 100 days, placing the $\gamma$-ray emission zone upstream of the radio-bright regions in the jet. Furthermore, Mrk 501 showed a historically low activity in X-rays and VHE $\gamma$-rays from mid-2017 to mid-2019 with a stable VHE flux ($>$0.2 TeV) of 5% the emission of the Crab Nebula. The broadband spectral energy distribution (SED) of this 2-year-long low-state, the potential baseline emission of Mrk 501, can be characterized with one-zone leptonic models, and with (lepto)-hadronic models fulfilling neutrino flux constraints from IceCube. We explore the time evolution of the SED towards the low-state, revealing that the stable baseline emission may be ascribed to a standing shock, and the variable emission to an additional expanding or traveling shock.Comment: 56 pages, 30 figures, 14 tables, submitted. Corresponding authors are L. Heckmann, D. Paneque, S. Gasparyan, M. Cerruti, and N. Sahakya

MAGIC very large zenith angle observations of the Crab Nebula up to 100 TeV

Aims. We measure the Crab Nebula gamma-ray spectral energy distribution in the 100 TeV energy domain and test the validity of existing leptonic emission models at these high energies.Methods. We used the novel very large zenith angle observations with the MAGIC telescope system to increase the collection area above 10 TeV. We also developed an auxiliary procedure of monitoring atmospheric transmission in order to assure proper calibration of the accumulated data. This employs recording optical images of the stellar field next to the source position, which provides a better than 10% accuracy for the transmission measurements.Results. We demonstrate that MAGIC very large zenith angle observations yield a collection area larger than a square kilometer. In only 56 h of observations, we detect the gamma-ray emission from the Crab Nebula up to 100 TeV, thus providing the highest energy measurement of this source to date with Imaging Atmospheric Cherenkov Telescopes. Comparing accumulated and archival MAGIC and Fermi/LAT data with some of the existing emission models, we find that none of them provides an accurate description of the 1 GeV to 100 TeV gamma-ray signal