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

Follow-up observations of GW170817 with the MAGIC telescopes

The discovery of the electromagnetic counterpart AT2017gfo and the GRB 170817A, associated to the binary neutron star merger GW170817, was one of the major advances in the study of gamma-ray bursts (GRBs) and the hallmark of the multi-messenger astronomy with gravitational waves. Another breakthrough in GRB physics is represented by the discovery of the highly energetic, teraelectronvolt (TeV) component in the GRB 190114C, possibly an universal component in all GRBs. This conclusion is also suggested by the hint of TeV emission in the short GRB 160821B and a few more events reported in the literature. The missing observational piece is the joint detection of TeV emission and gravitational waves from a short GRB and its progenitor. MAGIC observed the counterpart AT2017gfo as soon as the visibility conditions allowed it, namely from January to June 2018. These observations correspond to the maximum flux level observed in the radio and X-ray bands. The upper limits derived from TeV observations are compared with the modelling of the late non-thermal emission using the multi-frequency SED

Periastron Observations of TeV Gamma-Ray Emission from a Binary System with a 50-year Period

We report on observations of the pulsar / Be star binary system PSR J2032+4127 / MT91 213 in the energy range between 100 GeV and 20 TeV with the VERITAS and MAGIC imaging atmospheric Cherenkov telescope arrays. The binary orbit has a period of approximately 50 years, with the most recent periastron occurring on 2017 November 13. Our observations span from 18 months prior to periastron to one month after. A new, point-like, gamma-ray source is detected, coincident with the location of PSR J2032+4127 / MT91 213. The gamma-ray light curve and spectrum are well-characterized over the periastron passage. The flux is variable over at least an order of magnitude, peaking at periastron, thus providing a firm association of the TeV source with the pulsar / Be star system. Observations prior to periastron show a cutoff in the spectrum at an energy around 0.5 TeV. This result adds a new member to the small population of known TeV binaries, and it identifies only the second source of this class in which the nature and properties of the compact object are firmly established. We compare the gamma-ray results with the light curve measured with the X-ray Telescope (XRT) on board the Neil Gehrels Swift Observatory and with the predictions of recent theoretical models of the system. We conclude that significant revision of the models is required to explain the details of the emission we have observed, and we discuss the relationship between the binary system and the overlapping steady extended source, TeV J2032+4130

MAGIC observations of the nearby short GRB 160821B

Gamma-ray bursts (GRBs), the most luminous explosions in the universe, have at least two types known. One of them, short GRBs, have been thought to originate from binary neutron star (BNS) mergers. The discovery of GW170817 together with a GRB was the first and only direct proof of the hypothesis, and thus the properties of the short GRBs are poorly known yet. Aiming to clarify the underlying physical mechanisms of the short GRBs, we analyzed GRB 160821B, one of the nearest short GRBs known at z=0.162, observed with the MAGIC telescopes. A hint of a gamma-ray signal is found above 0.5 TeV at a significance of >3 sigma during observations from 24 seconds until 4 hours after the burst, as presented in the past. Recently, multi-wavelength data of its afterglow emission revealed a well-sampled kilonova component from a BNS merger, and the importance of GRB 160821B increased concerning GRB-GW studies. Accordingly, we investigated GRB afterglow models again, using the revised multi-wavelength data. We found that the straightforward interpretation with one-zone synchrotron self-Compton model from the external forward shock is in tension with the observed TeV flux, contradicting the suggestion reported previously. In this contribution we discuss the implication from the TeV observation, including alternative scenarios where the TeV emission can be enhanced. We also give a brief outlook of future GeV-TeV observations of short GRBs with imaging atmospheric Cherenkov telescopes, which could shed more light on the GRB-BNS merger relation

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

Multiwavelength monitoring of the gravitationally lensed blazar QSO B0218+357 between 2016 and 2020

QSO B0218+357 is currently the only gravitationally lensed source from which very-high-energy (VHE, &100GeV) gamma-ray emission has been detected. We report the multiwavelength monitoring observations of this source performed between 2016 and 2020 in radio interferometry, optical, X-ray and gamma-ray bands. During the monitoring individual flares and hints of enhanced states in optical, X-ray and GeV bands have been observed, and the simultaneous data taken by the MAGIC telescopes allow us to search for the VHE gamma-ray emission associated with these events

The Great Markarian 421 Flare of 2010 February: Multiwavelength Variability and Correlation Studies

We report on variability and correlation studies using multiwavelength observations of the blazar Mrk 421 during the month of 2010 February, when an extraordinary flare reaching a level of ∼27 Crab Units above 1 TeV was measured in very high energy (VHE) γ-rays with the Very Energetic Radiation Imaging Telescope Array System (VERITAS) observatory. This is the highest flux state for Mrk 421 ever observed in VHE γ-rays. Data are analyzed from a coordinated campaign across multiple instruments, including VHE γ-ray (VERITAS, Major Atmospheric Gamma-ray Imaging Cherenkov), high-energy γ-ray (Fermi-LAT), X-ray (Swift, Rossi X-ray Timing Experiment, MAXI), optical (including the GASP-WEBT collaboration and polarization data), and radio (Metsahovi, Owens Valley Radio Observatory, University of Michigan Radio Astronomy Observatory). Light curves are produced spanning multiple days before and after the peak of the VHE flare, including over several flare "decline" epochs. The main flare statistics allow 2 minute time bins to be constructed in both the VHE and optical bands enabling a cross-correlation analysis that shows evidence for an optical lag of ∼25-55 minutes, the first time-lagged correlation between these bands reported on such short timescales. Limits on the Doppler factor (δ ⪆ 33) and the size of the emission region (δ-1RB≲ 3.8 × 1013cm) are obtained from the fast variability observed by VERITAS during the main flare. Analysis of 10 minute binned VHE and X-ray data over the decline epochs shows an extraordinary range of behavior in the flux-flux relationship, from linear to quadratic to lack of correlation to anticorrelation. Taken together, these detailed observations of an unprecedented flare seen in Mrk 421 are difficult to explain with the classic single-zone synchrotron self-Compton model

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]

Performance and first measurements of the MAGIC stellar intensity interferometer

In recent years, a new generation of optical intensity interferometers has emerged, leveraging the existing infrastructure of Imaging Atmospheric Cherenkov Telescopes (IACTs). The MAGIC telescopes host the MAGIC-SII system (Stellar Intensity Interferometer), implemented to investigate the feasibility and potential of this technique on IACTs. After the first successful measurements in 2019, the system was upgraded and now features a real-time, dead-time-free, 4-channel, GPU-based correlator. These hardware modifications allow seamless transitions between MAGIC’s standard very-high-energy gamma-ray observations and optical interferometry measurements within seconds. We establish the feasibility and potential of employing IACTs as competitive optical Intensity Interferometers with minimal hardware adjustments. The measurement of a total of 22 stellar diameters are reported, 9 corresponding to reference stars with previous comparable measurements, and 13 with no prior measurements. A prospective implementation involving telescopes from the forthcoming Cherenkov Telescope Array Observatory’s Northern hemisphere array, such as the first prototype of its Large-Sized Telescopes, LST-1, is technically viable. This integration would significantly enhance the sensitivity of the current system and broaden the UV-plane coverage. This advancement would enable the system to achieve competitive sensitivity with the current generation of long-baseline optical interferometers over blue wavelengths