ALBERT

Description: 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.

Description: 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.

Description: 1ES 1959+650 is a bright TeV high-frequency-peaked BL Lac object exhibiting interesting features like “orphan” TeV flares and broad emission in the high-energy regime that are difficult to interpret using conventional one-zone Synchrotron Self-Compton (SSC) scenarios. We report the results from the Major Atmospheric Gamma Imaging Cherenkov (MAGIC) observations in 2016 along with the multi-wavelength data from the Fermi Large Area Telescope (LAT) and Swift instruments. MAGIC observed 1ES 1959+650 with different emission levels in the very-high-energy (VHE, E 〉  100 GeV) γ-ray band during 2016. In the long-term data, the X-ray spectrum becomes harder with increasing flux and a hint of a similar trend is also visible in the VHE band. An exceptionally high VHE flux reaching ∼3 times the Crab Nebula flux was measured by MAGIC on the 13 and 14 of June, and 1 July 2016 (the highest flux observed since 2002). During these flares, the high-energy peak of the spectral energy distribution (SED) lies in the VHE domain and extends up to several TeV. The spectrum in the γ-ray (both Fermi-LAT and VHE bands) and the X-ray bands are quite hard. On 13 June and 1 July 2016, the source showed rapid variations in the VHE flux within timescales of less than an hour. A simple one-zone SSC model can describe the data during the flares requiring moderate to large values of the Doppler factors (δ ≥ 30−60). Alternatively, the high-energy peak of the SED can be explained by a purely hadronic model attributed to proton-synchrotron radiation with jet power Ljet ∼ 1046 erg s−1 and under high values of the magnetic field strength (∼100 G) and maximum proton energy (∼few EeV). Mixed lepto-hadronic models require super-Eddington values of the jet power. We conclude that it is difficult to get detectable neutrino emission from the source during the extreme VHE flaring period of 2016.

Description: Aims. We aim to characterize the multiwavelength emission from Markarian 501 (Mrk 501), quantify the energy-dependent variability, study the potential multiband correlations, and describe the temporal evolution of the broadband emission within leptonic theoretical scenarios. Methods. We organized a multiwavelength campaign to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration. Results. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0.2 TeV of ∼0.5 times the Crab Nebula flux (CU) for most of the campaign. The highest activity occurred on 2012 June 9, when the VHE flux was ∼3 CU, and the peak of the high-energy spectral component was found to be at ∼2 TeV. Both the X-ray and VHE gamma-ray spectral slopes were measured to be extremely hard, with spectral indices 〈 2 during most of the observing campaign, regardless of the X-ray and VHE flux. This study reports the hardest Mrk 501 VHE spectra measured to date. The fractional variability was found to increase with energy, with the highest variability occurring at VHE. Using the complete data set, we found correlation between the X-ray and VHE bands; however, if the June 9 flare is excluded, the correlation disappears (significance 〈 3σ) despite the existence of substantial variability in the X-ray and VHE bands throughout the campaign. Conclusions. The unprecedentedly hard X-ray and VHE spectra measured imply that their low- and high-energy components peaked above 5 keV and 0.5 TeV, respectively, during a large fraction of the observing campaign, and hence that Mrk 501 behaved like an extreme high-frequency-peaked blazar (EHBL) throughout the 2012 observing season. This suggests that being an EHBL may not be a permanent characteristic of a blazar, but rather a state which may change over time. The data set acquired shows that the broadband spectral energy distribution (SED) of Mrk 501, and its transient evolution, is very complex, requiring, within the framework of synchrotron self-Compton (SSC) models, various emission regions for a satisfactory description. Nevertheless the one-zone SSC scenario can successfully describe the segments of the SED where most energy is emitted, with a significant correlation between the electron energy density and the VHE gamma-ray activity, suggesting that most of the variability may be explained by the injection of high-energy electrons. The one-zone SSC scenario used reproduces the behavior seen between the measured X-ray and VHE gamma-ray fluxes, and predicts that the correlation becomes stronger with increasing energy of the X-rays.

Description: Aims. The very high energy (VHE ≳100 GeV) γ-ray MAGIC observations of the blazar S4 0954+65, were triggered by an exceptionally high flux state of emission in the optical. This blazar has a disputed redshift of z = 0.368 or z ≥ 0.45 and an uncertain classification among blazar subclasses. The exceptional source state described here makes for an excellent opportunity to understand physical processes in the jet of S4 0954+65 and thus contribute to its classification. Methods. We investigated the multiwavelength (MWL) light curve and spectral energy distribution (SED) of the S4 0954+65 blazar during an enhanced state in February 2015 and have put it in context with possible emission scenarios. We collected photometric data in radio, optical, X-ray, and γ-ray. We studied both the optical polarization and the inner parsec-scale jet behavior with 43 GHz data. Results. Observations with the MAGIC telescopes led to the first detection of S4 0954+65 at VHE. Simultaneous data with Fermi-LAT at high energy γ-ray (HE, 100 MeV 〈 E 〈 100 GeV) also show a period of increased activity. Imaging at 43 GHz reveals the emergence of a new feature in the radio jet in coincidence with the VHE flare. Simultaneous monitoring of the optical polarization angle reveals a rotation of approximately 100°. Conclusions. The high emission state during the flare allows us to compile the simultaneous broadband SED and to characterize it in the scope of blazar jet emission models. The broadband spectrum can be modeled with an emission mechanism commonly invoked for flat spectrum radio quasars (FSRQs), that is, inverse Compton scattering on an external soft photon fieldfrom the dust torus, also known as external Compton. The light curve and SED phenomenology is consistent with an interpretation of a blob propagating through a helical structured magnetic field and eventually crossing a standing shock in the jet, a scenario typically applied to FSRQs and low-frequency peaked BL Lac objects (LBL).

Description: The mechanisms producing fast variability of the γ-ray emission in active galactic nuclei (AGNs) are under debate. The MAGIC telescopes detected a fast, very-high-energy (VHE, E  〉  100 GeV) γ-ray flare from BL Lacertae on 2015 June 15. The flare had a maximum flux of (1.5 ± 0.3) × 10−10 photons cm−2 s−1 and halving time of 26 ± 8 min. The MAGIC observations were triggered by a high state in the optical and high-energy (HE, E  〉  100 MeV) γ-ray bands. In this paper we present the MAGIC VHE γ-ray data together with multi-wavelength data from radio, optical, X-rays, and HE γ rays from 2015 May 1 to July 31. Well-sampled multi-wavelength data allow us to study the variability in detail and compare it to the other epochs when fast, VHE γ-ray flares have been detected from this source. Interestingly, we find that the behaviour in radio, optical, X-rays, and HE γ-rays is very similar to two other observed VHE γ-ray flares. In particular, also during this flare there was an indication of rotation of the optical polarization angle and of activity at the 43 GHz core. These repeating patterns indicate a connection between the three events. We also test modelling of the spectral energy distribution based on constraints from the light curves and VLBA observations, with two different geometrical setups of two-zone inverse Compton models. In addition we model the γ-ray data with the star-jet interaction model. We find that all of the tested emission models are compatible with the fast VHE γ-ray flare, but all have some tension with the multi-wavelength observations.

Description: 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.

Description: Aims. The aim of this study is to search for evidence of a common emission engine between radio giant pulses (GPs) and very-high-energy (VHE, E 〉  100 GeV) γ-rays from the Crab pulsar. Methods. We performed 16 h of simultaneous observations of the Crab pulsar at 1.4 GHz with the Effelsberg radio telescope and the Westerbork Synthesis Radio Telescope (WSRT), and at energies above 60 GeV we used the Major Atmospheric Gamma-ray Imaging Cherenkov (MAGIC) telescopes. We searched for a statistical correlation between the radio and VHE γ-ray emission with search windows of different lengths and different time lags to the arrival times of a radio GP. A dedicated search for an enhancement in the number of VHE γ-rays correlated with the occurrence of radio GPs was carried out separately for the P1 and P2 phase ranges, respectively. Results. In the radio data sample, 99444 radio GPs were detected. We find no significant correlation between the GPs and VHE photons in any of the search windows. Depending on phase cuts and the chosen search windows, we find upper limits at a 95% confidence level on an increase in VHE γ-ray events correlated with radio GPs between 7% and 61% of the average Crab pulsar VHE flux for the P1 and P2 phase ranges, respectively. This puts upper limits on the flux increase during a radio GP between 12% and 2900% of the pulsed VHE flux, depending on the search window duration and phase cuts. This is the most stringent upper limit on a correlation between γ-ray emission and radio GPs reported so far.

Description: 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.

Description: 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.