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

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

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

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

Multiwavelength variability and correlation studies of Mrk 421 during historically low X-ray and γ-ray activity in 2015-2016

We report a characterization of the multiband flux variability and correlations of the nearby (z = 0.031) blazar Markarian 421 (Mrk 421) using data from Metsahovi, Swift, Fermi-LAT, MAGIC, FACT, and other collaborations and instruments from 2014 November till 2016 June. Mrk 421 did not show any prominent flaring activity, but exhibited periods of historically low activity above 1 TeV (F->1 TeV 0.1 TeV) gamma-rays, which, despite the low activity, show a significant positive correlation with no time lag. The HRkeV and HRTeV show the harder-when-brighter trend observed in many blazars, but the trend flattens at the highest fluxes, which suggests a change in the processes dominating the blazar variability. Enlarging our data set with data from years 2007 to 2014, we measured a positive correlation between the optical and the GeV emission over a range of about 60 d centred at time lag zero, and a positive correlation between the optical/GeV and the radio emission over a range of about 60 d centred at a time lag of 43(-6)(+9) d. This observation is consistent with the radio-bright zone being located about 0.2 parsec downstream from the optical/GeV emission regions of the jet. The flux distributions are better described with a lognormal function in most of the energy bands probed, indicating that the variability in Mrk 421 is likely produced by a multiplicative process

Searching for GRBs at VHE with MAGIC: the status before CTA

Gamma-Ray Bursts (GRBs) are one of the main targets for current and next generation Imaging Atmospheric Cherenkov Telescopes (IACTs). Given their transient behavior, especially in the case of their prompt emission phase, performing fast follow-up observations is challenging for IACTs, which have a narrow field of view and limited duty cycle. Despite this, MAGIC plays a major role in the search for Very High Energy (VHE, E>100 GeV) gamma-ray emission from GRBs: this is possible thanks to its fast repositioning speed, low energy threshold and high sensitivity at the lowest energies. In 2013 the MAGIC GRB automatic procedure was upgraded, increasing the number of GRBs followed in the prompt and early afterglow phases and decreasing dramatically hardware failures during fast repositioning. Currently, only GRB 190114C was firmly detected in the VHE band, while for other GRBs no significant detection was achieved. In such a case, upper limits (ULs) can give insight into the physics driving such eluding sources, especially on their emission mechanisms. In this contribution we report on the status of the GRB follow-up with MAGIC and focus on the ULs and results obtained from a sample of GRBs observed between 2013 and 2018. This GRB catalog is the result of the MAGIC well-designed and tested follow-up procedure, and it serves as a precursor of GRBs observation with the next generation IACT system, the Cherenkov Telescope Array (CTA)

Following up Transient Sources at Very High Energies with MAGIC

Several classes of sources are known to emit different messengers. Among them, transient sources are a special case, due to their serendipitous occurrence, time variability and duration on different timescales. They are associated with explosive and catastrophic events where very compact objects like neutron stars and black holes are involved. The difficulty of observing such elusive and possibly short-lasting events requires a fast reaction and a well-organized alert network between different experiments. In order to characterize them in the best possible way, instruments with a wide field of view should serve as external triggers for facilities with small sky coverage. MAGIC, as a Cherenkov telescope, belongs to the latter category. The search for transients by MAGIC is possible thanks to an automatic alert system listening to the alerts sent by the Gamma-ray Coordinate Network (GCN). In this contribution we describe the MAGIC alert system, which was designed mainly for the follow-up of Gamma-Ray Bursts in its initial conception. The alert system was recently updated in a multi-messenger context, receiving alerts also from neutrino and GW observatories. Finally we will present the MAGIC program for transient sources and how it was adapted in the current multi-wavelength and multi-messenger panorama

MAGIC observations of the diffuse γ -ray emission in the vicinity of the Galactic center

Aims. In the presence of a sufficient amount of target material, γ-rays can be used as a tracer in the search for sources of Galactic cosmic rays (CRs). Here we present deep observations of the Galactic center (GC) region with the MAGIC telescopes and use them to infer the underlying CR distribution and to study the alleged PeV proton accelerator at the center of our Galaxy. Methods. We used data from ≈100 h observations of the GC region conducted with the MAGIC telescopes over five years (from 2012 to 2017). Those were collected at high zenith angles (58-70 deg), leading to a larger energy threshold, but also an increased effective collection area compared to low zenith observations. Using recently developed software tools, we derived the instrument response and background models required for extracting the diffuse emission in the region. We used existing measurements of the gas distribution in the GC region to derive the underlying distribution of CRs. We present a discussion of the associated biases and limitations of such an approach. Results. We obtain a significant detection for all four model components used to fit our data (Sgr A∗, "Arc", G0.9+0.1, and an extended component for the Galactic Ridge). We observe no significant difference between the γ-ray spectra of the immediate GC surroundings, which we model as a point source (Sgr A∗) and the Galactic Ridge. The latter can be described as a power-law with index 2 and an exponential cut-off at around 20 TeV with the significance of the cut-off being only 2σ. The derived cosmic-ray profile hints to a peak at the GC position and with a measured profile index of 1.2 ± 0.3 is consistent with the 1/r radial distance scaling law, which supports the hypothesis of a CR accelerator at the GC. We argue that the measurements of this profile are presently limited by our knowledge of the gas distribution in the GC vicinity

Monitoring the magnetar SGR 1935+2154 with the MAGIC telescopes

The Galactic magnetar SGR 1935+2154 was associated with a bright, millisecond-timescale fast radio burst (FRB) which occured in April 2020, during a flaring episode. This was the first time an FRB was unequivocally associated with a Galactic source, and the first FRB for which the nature of the emitting source was identified. Moreover, it was the first FRB with a counterpart at another wavelength correlated in time, an atypical, hard X-ray burst, which provides clear evidence for accompanying non-thermal processes. The MAGIC Telescopes are Imaging Air Cherenkov Telescopes (IACTs) sensitive to very-high-energy (VHE, E>100 GeV) gamma rays. Located at the center of the camera lies the MAGIC Central pixel, a single fully-modified photosensor-toreadout chain to measure millisecond-duration optical signals, displaying a maximum sensitivity at a wavelength of 350 nm. This allows MAGIC to operate simultaneously both as a VHE gammaray and a fast optical telescope. The MAGIC telescopes have monitored SGR 1935+2154 in a multiwavelength campaign involving X-ray, radio and optical facilities. In this contribution, we will show the results on the search for the VHE counterpart of the first SGR-FRB source in this multiwavelength context, as well as the search for fast optical bursts with the MAGIC Central Pixel

A search for dark matter in Triangulum II with the MAGIC telescopes

We present the first results from very-high-energy observations of the dwarf spheroidal satellite candidate Triangulum II with the MAGIC telescopes from 62.4 h of good-quality data taken between August 2016 and August 2017. We find no gamma-ray excess in the direction of Triangulum II, and upper limits on both the differential and integral gamma-ray flux are presented. Currently, the kinematics of Triangulum II are affected by large uncertainties leading to a bias in the determination of the properties of its dark matter halo. Using a scaling relation between the annihilation J-factor and heliocentric distance of well-known dwarf spheroidal galaxies, we estimate an annihilation J-factor for Triangulum II for WIMP dark matter of log[Jann(0.5°)∕GeV2cm−5]=19.35±0.37. We also derive a dark matter density profile for the object relying on results from resolved simulations of Milky Way sized dark matter halos. We obtain 95% confidence-level limits on the thermally averaged annihilation cross section for WIMP annihilation into various Standard Model channels. The most stringent limits are obtained in the τ−τ+ final state, where a cross section for annihilation down to 〈σannv〉=3.05×10−24 cm3 s−1 is excluded