The current generation of ground-based imaging atmospheric Cherenkov telescopes (IACTs) operate in the very-high-energy (VHE) domain from ~100 GeV to ~100 TeV. They use electronic digital trigger ...systems to discern the Cherenkov light flashes emitted by extensive air showers (EASs), from the overwhelming light of the night sky (LoNS) background. Near the telescope energy threshold, the number of emitted Cherenkov photons by gamma-ray-induced EASs is comparable to the fluctuations of the LoNS and the photon distribution at the Cherenkov-imaging camera plane becomes patchy. This results in a severe loss of effectiveness of the digital triggers based on combinatorial logic of thresholded signals. A stereoscopic analog trigger system has been developed for improving the detection capabilities of the Major Atmospheric Gamma-ray Imaging Cherenkov (MAGIC) telescopes at the lowest energies. It is based on the analog sum of the photosensor electrical signals. In this article, the architectural design, technical performances, and configuration of this stereoscopic analog trigger, dubbed " Sum-Trigger-II ," are described.
Context.
Extended and delayed emission around distant TeV sources induced by the effects of propagation of
γ
ray s through the intergalactic medium can be used for the measurement of the ...intergalactic magnetic field (IGMF).
Aims.
We search for delayed GeV emission from the hard-spectrum TeV
γ
-ray emitting blazar 1ES 0229+200, with the goal of detecting or constraining the IGMF-dependent secondary flux generated during the propagation of TeV
γ
rays through the intergalactic medium.
Methods.
We analysed the most recent MAGIC observations over a 5 year time span, and complemented 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 used them to trace source evolution in the GeV–TeV band over a decade and a half. We used Monte Carlo simulations to predict the delayed secondary
γ
-ray flux, modulated by the source variability, as revealed by TeV-band observations. We then compared these predictions for various assumed IGMF strengths to all available measurements of the
γ
-ray flux evolution.
Results.
We find that the source flux in the energy range above 200 GeV experiences variations around its average on the 14-year time span of observations. No evidence for the flux variability is found in the 1 − 100 GeV energy range accessible to
Fermi
/LAT. The non-detection of variability due to delayed emission from electromagnetic cascade developing in the intergalactic medium imposes a lower bound of
B
> 1.8 × 10
−17
G for the long-correlation-length IGMF and
B
> 10
−14
G for an IGMF of 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 that might explain the baryon asymmetry of the Universe.
Context.
Markarian 501 (Mrk 501) is a very high-energy (VHE) gamma-ray blazar located at
z
= 0.034, which is regularly monitored by a wide range of multi-wavelength instruments, from radio to VHE ...gamma rays. During a period of almost two weeks in July 2014, the highest X-ray activity of Mrk 501 was observed in ∼14 years of operation of the
Neil Gehrels Swift
Gamma-ray Burst Observatory.
Aims.
We characterize the broadband variability of Mrk 501 from radio to VHE gamma rays during the most extreme X-ray activity measured in the last 14 years, and evaluate whether it can be interpreted within theoretical scenarios widely used to explain the broadband emission from blazars.
Methods.
The emission of Mrk 501 was measured at radio with Metsähovi, at optical–UV with KVA and
Swift
/UVOT, at X-ray with
Swift
/XRT and
Swift
/BAT, at gamma ray with
Fermi
-LAT, and at VHE gamma rays with the FACT and MAGIC telescopes. The multi-band variability and correlations were quantified, and the broadband spectral energy distributions (SEDs) were compared with predictions from theoretical models.
Results.
The VHE emission of Mrk 501 was found to be elevated during the X-ray outburst, with a gamma-ray flux above 0.15 TeV varying from ∼0.5 to ∼2 times the Crab nebula flux. The X-ray and VHE emission both varied on timescales of 1 day and were found to be correlated. We measured a general increase in the fractional variability with energy, with the VHE variability being twice as large as the X-ray variability. The temporal evolution of the most prominent and variable segments of the SED, characterized on a day-by-day basis from 2014 July 16 to 2014 July 31, is described with a one-zone synchrotron self-Compton model with variations in the break energy of the electron energy distribution (EED), and with some adjustments in the magnetic field strength and spectral shape of the EED. These results suggest that the main flux variations during this extreme X-ray outburst are produced by the acceleration and the cooling of the high-energy electrons. A narrow feature at ∼3 TeV was observed in the VHE spectrum measured on 2014 July 19 (MJD 56857.98), which is the day with the highest X-ray flux (>0.3 keV) measured during the entire
Swift
mission. This feature is inconsistent with the classical analytic functions to describe the measured VHE spectra (power law, log-parabola, and log-parabola with exponential cutoff) at more than 3
σ
. A fit with a log-parabola plus a narrow component is preferred over the fit with a single log-parabola at more than 4
σ
, and a dedicated Monte Carlo simulation estimated the significance of this extra component to be larger than 3
σ
. Under the assumption that this VHE spectral feature is real, we show that it can be reproduced with three distinct theoretical scenarios: (a) a pileup in the EED due to stochastic acceleration; (b) a structured jet with two-SSC emitting regions, with one region dominated by an extremely narrow EED; and (c) an emission from an IC pair cascade induced by electrons accelerated in a magnetospheric vacuum gap, in addition to the SSC emission from a more conventional region along the jet of Mrk 501.
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.
Aims.
We measure the Crab Nebula
γ
-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
γ
-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
γ
-ray signal.
Context. PKS 1510–089 is a flat spectrum radio quasar strongly variable in the optical and GeV range. To date, very high-energy (VHE, > 100 GeV) emission has been observed from this source either ...during long high states of optical and GeV activity or during short flares. Aims. We search for low-state VHE gamma-ray emission from PKS 1510–089. We characterize and model the source in a broadband context, which would provide a baseline over which high states and flares could be better understood. Methods. PKS 1510–089 has been monitored by the MAGIC telescopes since 2012. We use daily binned Fermi-LAT flux measurements of PKS 1510–089 to characterize the GeV emission and select the observation periods of MAGIC during low state of activity. For the selected times we compute the average radio, IR, optical, UV, X-ray, and gamma-ray emission to construct a low-state spectral energy distribution of the source. The broadband emission is modeled within an external Compton scenario with a stationary emission region through which plasma and magnetic fields are flowing. We also perform the emission-model-independent calculations of the maximum absorption in the broad line region (BLR) using two different models. Results. The MAGIC telescopes collected 75 hr of data during times when the Fermi-LAT flux measured above 1 GeV was below 3 × 10−8 cm−2 s−1, which is the threshold adopted for the definition of a low gamma-ray activity state. The data show a strongly significant (9.5σ) VHE gamma-ray emission at the level of (4.27 ± 0.61stat) × 10−12 cm−2 s−1 above 150 GeV, a factor of 80 lower than the highest flare observed so far from this object. Despite the lower flux, the spectral shape is consistent with earlier detections in the VHE band. The broadband emission is compatible with the external Compton scenario assuming a large emission region located beyond the BLR. For the first time the gamma-ray data allow us to place a limit on the location of the emission region during a low gamma-ray state of a FSRQ. For the used model of the BLR, the 95% confidence level on the location of the emission region allows us to place it at a distance > 74% of the outer radius of the BLR.
We report on the detection of flaring activity from the Fanaroff-Riley I radio galaxy NGC 1275 in very-high-energy (VHE, E > 100 GeV) gamma rays with the Major Atmospheric Gamma Imaging Cherenkov ...(MAGIC) telescopes. The observations were performed between 2016 September and 2017 February, as part of a monitoring programme. The brightest outburst, with ∼1.5 times the Crab Nebula flux above 100 GeV (C.U.), was observed during the night between 2016 December 31 and 2017 January 1. The flux is fifty times higher than the mean flux previously measured in two observational campaigns between 2009 October and 2010 February and between 2010 August and 2011 February. Significant variability of the day-by-day light curve was measured. The shortest flux-doubling timescale was found to be of (611 ± 101) min. The spectra calculated for this period are harder and show a significant curvature with respect to the ones obtained in the previous campaigns. The combined spectrum of the MAGIC data during the strongest flare state and simultaneous data from the Fermi-LAT around 2017 January 1 follows a power law with an exponential cutoff at the energy (492 ± 35) GeV. We further present simultaneous optical flux density measurements in the R-band obtained with the Kungliga Vetenskaps Akademien (KVA) telescope and investigate the correlation between the optical and gamma-ray emission. Due to possible internal pair-production, the fast flux variability constrains the Doppler factor to values that are inconsistent with a large viewing angle as observed in the radio band. We investigate different scenarios for the explanation of fast gamma-ray variability, namely emission from magnetospheric gaps, relativistic blobs propagating in the jet (mini-jets), or an external cloud (or star) entering the jet. We find that the only plausible model to account for the luminosities here observed would be the production of gamma rays in a magnetospheric gap around the central black hole, only in the eventuality of an enhancement of the magnetic field threading the hole from its equipartition value with the gas pressure in the accretion flow. The observed gamma-ray flare therefore challenges all the discussed models for fast variability of VHE gamma-ray emission in active galactic nuclei.
Context. The BL Lac object S5 0716+714, a highly variable blazar, underwent an impressive outburst in January 2015 (Phase A), followed by minor activity in February (Phase B). The MAGIC observations ...were triggered by the optical flux observed in Phase A, corresponding to the brightest ever reported state of the source in the R-band. Aims.The comprehensive dataset collected is investigated in order to shed light on the mechanism of the broadband emission. Methods. Multi-wavelength light curves have been studied together with the broadband spectral energy distributions (SEDs). The sample includes data from Effelsberg, OVRO, Metsähovi, VLBI, CARMA, IRAM, SMA, Swift-UVOT, KVA, Tuorla, Steward, RINGO3, KANATA, AZT-8+ST7, Perkins, LX-200, Swift-XRT, NuSTAR, Fermi-LAT and MAGIC. Results. The flaring state of Phase A was detected in all the energy bands, providing for the first time a multi-wavelength sample of simultaneous data from the radio band to the very-high-energy (VHE, E > 100 GeV). In the constructed SED, the Swift-XRT+NuSTAR data constrain the transition between the synchrotron and inverse Compton components very accurately, while the second peak is constrained from 0.1 GeV to 600 GeV by Fermi+MAGIC data. The broadband SED cannot be described with a one-zone synchrotron self-Compton model as it severely underestimates the optical flux in order to reproduce the X-ray to γ-ray data. Instead we use a two-zone model. The electric vector position angle (EVPA) shows an unprecedented fast rotation. An estimation of the redshift of the source by combined high-energy (HE, 0.1 GeV < E < 100 GeV) and VHE data provides a value of z = 0.31 ± 0.02stats ± 0.05sys, confirming the literature value. Conclusions. The data show the VHE emission originating in the entrance and exit of a superluminal knot in and out of a recollimation shock in the inner jet. A shock–shock interaction in the jet seems responsible for the observed flares and EVPA swing. This scenario is also consistent with the SED modeling.