A neutrino with energy ∼290 TeV, IceCube-170922A, was detected in coincidence with the BL Lac object TXS 0506+056 during enhanced gamma-ray activity, with chance coincidence being rejected at ∼3 ...level. We monitored the object in the very-high-energy (VHE) band with the Major Atmospheric Gamma-ray Imaging Cherenkov (MAGIC) telescopes for ∼41 hr from 1.3 to 40.4 days after the neutrino detection. Day-timescale variability is clearly resolved. We interpret the quasi-simultaneous neutrino and broadband electromagnetic observations with a novel one-zone lepto-hadronic model, based on interactions of electrons and protons co-accelerated in the jet with external photons originating from a slow-moving plasma sheath surrounding the faster jet spine. We can reproduce the multiwavelength spectra of TXS 0506+056 with neutrino rate and energy compatible with IceCube-170922A, and with plausible values for the jet power of . The steep spectrum observed by MAGIC is concordant with internal γγ absorption above ∼100 GeV entailed by photohadronic production of a ∼290 TeV neutrino, corroborating a genuine connection between the multi-messenger signals. In contrast to previous predictions of predominantly hadronic emission from neutrino sources, the gamma-rays can be mostly ascribed to inverse Compton upscattering of external photons by accelerated electrons. The X-ray and VHE bands provide crucial constraints on the emission from both accelerated electrons and protons. We infer that the maximum energy of protons in the jet comoving frame can be in the range ∼1014 - 1018 eV.
ABSTRACT
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–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 that 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 γγ 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σ upper bounds, including systematic uncertainties, are between 13 per cent and 23 per cent 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\mathrm{ m}$ range relative to the studied models, yet compatible with them within systematics.
Abstract
It is widely believed that the bulk of the Galactic cosmic rays is accelerated in supernova remnants (SNRs). However, no observational evidence of the presence of particles of PeV energies ...in SNRs has yet been found. The young historical SNR Cassiopeia A (Cas A) appears as one of the best candidates to study acceleration processes. Between 2014 December and 2016 October, we observed Cas A with the MAGIC telescopes, accumulating 158 h of good quality data. We derived the spectrum of the source from 100 GeV to 10 TeV. We also analysed ∼8 yr of Fermi-LAT to obtain the spectral shape between 60 MeV and 500 GeV. The spectra measured by the LAT and MAGIC telescopes are compatible within the errors and show a clear turn-off (4.6σ) at the highest energies, which can be described with an exponential cut-off at
$E_c = 3.5(^{+1.6}_{-1.0})_{{\rm stat}} (^{+0.8}_{-0.9})_{{\rm sys}}$
TeV. The gamma-ray emission from 60 MeV to 10 TeV can be attributed to a population of high-energy protons with a spectral index of ∼2.2 and an energy cut-off at ∼10 TeV. This result indicates that Cas A is not contributing to the high energy (∼PeV) cosmic ray sea in a significant manner at the present moment. A one-zone leptonic model fails to reproduce by itself the multiwavelength spectral energy distribution. Besides, if a non-negligible fraction of the flux seen by MAGIC is produced by leptons, the radiation should be emitted in a region with a low magnetic field (B⪅180 μG) like in the reverse shock.
Abstract
Extreme high-frequency-peaked BL Lac objects (EHBLs) are blazars that exhibit extremely energetic synchrotron emission. They also feature nonthermal gamma-ray emission whose peak lies in the ...very high-energy (VHE,
E
> 100 GeV) range, and in some sources exceeds 1 TeV: this is the case for hard-TeV EHBLs such as 1ES 0229+200. With the aim of increasing the EHBL population, 10 targets were observed with the MAGIC telescopes from 2010 to 2017, for a total of 265 hr of good-quality data. The data were complemented by coordinated
Swift
observations. The X-ray data analysis confirms that all but two sources are EHBLs. The sources show only a modest variability and a harder-when-brighter behavior, typical for this class of objects. At VHE gamma-rays, three new sources were detected and a hint of a signal was found for another new source. In each case, the intrinsic spectrum is compatible with the hypothesis of a hard-TeV nature of these EHBLs. The broadband spectral energy distributions (SEDs) of all sources are built and modeled in the framework of a single-zone, purely leptonic model. The VHE gamma-ray-detected sources were also interpreted with a spine–layer model and a proton synchrotron model. The three models provide a good description of the SEDs. However, the resulting parameters differ substantially in the three scenarios, in particular the magnetization parameter. This work presents the first mini catalog of VHE gamma-ray and multiwavelength observations of EHBLs.
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.
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.
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 observations of the pulsar/Be star binary system PSR J2032+4127/MT91 213 in the energy range between and with the Very Energetic Radiation Imaging Telescope Array and Major Atmospheric ...Gamma Imaging 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 . 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 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 that we have observed, and we discuss the relationship between the binary system and the overlapping steady extended source, TeV J2032+4130.