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 coaccelerated 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 10 4 10 erg s 45 46 1 -´-. 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 ∼10 14 -10 18 eV.

The Cherenkov Telescope Array (CTA), the new-generation ground-based observatory for γ astronomy, provides unique capabilities to address significant open questions in astrophysics, cosmology, and fundamental physics. We study some of the salient areas of γ cosmology that can be explored as part of the Key Science Projects of CTA, through simulated observations of active galactic nuclei (AGN) and of their relativistic jets. Observations of AGN with CTA will enable a measurement of γ absorption on the extragalactic background light with a statistical uncertainty below 15% up to a redshift z=2 and to constrain or detect γ halos up to intergalactic-magnetic-field strengths of at least 0.3 pG . Extragalactic observations with CTA also show promising potential to probe physics beyond the Standard Model. The best limits on Lorentz invariance violation from γ astronomy will be improved by a factor of at least two to three. CTA will also probe the parameter space in which axion-like particles could constitute a significant fraction, if not all, of dark matter. We conclude on the synergies between CTA and other upcoming facilities that will foster the growth of γ cosmology.

The analysis and combination of data from different gamma-ray instruments involves the use of collaboration proprietary software and case-bycase methods. The effort of defining a common data format for high-level data, namely event lists and instrument response functions (IRFs), has recently started for very-high-energy gamma-ray instruments, driven by the upcoming Cherenkov Telescope Array (CTA). In this work we implemented this prototypical data format for a small set of MAGIC, VERITAS, FACT, and H.E.S.S. Crab nebula observations, and we analyzed them with the open-source gammapy software package. By combining data from Fermi-LAT, and from four of the currently operating imaging atmospheric Cherenkov telescopes, we produced a joint maximum likelihood fit of the Crab nebula spectrum. Aspects of the statistical errors and the evaluation of systematic uncertainty are also commented upon, along with the release format of spectral measurements. The results presented in this work are obtained using open-access on-line assets that allow for a long-term reproducibility of the results.

Aims. The nearby (z = 0.031) TeV blazar Mrk 421 was reported to be in a high state of flux activity since November, 2009. We aim to investigate possible changes in the physical parameters of Mrk 421 during its high state of activity using multiwavelength data. Methods. We have observed this source in the bright state using the High Altitude GAmma Ray (HAGAR) telescope array at energies above 250 GeV during February 13-19, 2010. Optical, X-ray and γ-ray archival data were also used to obtain the spectral energy distribution and light curves. Results. Mrk 421 was found to undergo one of its brightest flaring episodes on February 17, 2010 by various observations in X-rays and γ-rays. HAGAR observations during February 13-19, 2010 at energies above 250 GeV show an enhancement in the flux level, with a maximum flux of ∼7 Crab units being detected on February 17, 2010. We present the spectral energy distributions during this flaring episode and investigate the correlation of the variability in X-ray and γ-ray bands. Conclusions. Our multiwavelength study suggests that the flare detected during February 16 and 17, 2010 may have been caused by a passing shock in the jet.

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 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 γγ 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% 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 µm range relative to the studied models, yet compatible with them within systematics.

Context. The HAGAR Telescope Array at Hanle, Ladakh has been regularly monitoring the nearby blazar Mkn 421 for the past seven years. Aims. Blazars show flux variability in all timescales across the electromagnetic spectrum. While there is abundant literature characterizing the short-term flares from different blazars, comparatively little work has been carried out to study the long-term variability. We aim to study the long-term temporal and spectral variability in the radiation from Mkn 421 during 2009−2015. Methods. We quantify the variability and lognormality from the radio to very high-energy (VHE) bands, and compute the correlations between various wavebands using the z-transformed discrete correlation function. We construct the spectral energy distribution (SED) contemporaneous with HAGAR observation seasons and fit this SED with a one-zone synchrotron self-Compton model to study the spectral variability. Results. The flux is found to be highly variable across all timescales. The variability is energy dependant and is maximum in the X-ray and VHE bands. A strong correlation is found between the Fermi-LAT (gamma) and radio bands and between Fermi-LAT and optical, but no correlation is found between Fermi-LAT and X-ray. Lognormality in the flux distribution is clearly detected. This is the third blazar, following BL Lac and PKS 2155+304 to show this behaviour. The SED can be well fit by a one-zone SSC model, and variations in the flux states can be attributed mainly to changes in the particle distribution. A strong correlation is seen between the break energy γ b of the particle spectrum and the total bolometric luminosity.

The coincident detection of GW170817 in gravitational waves and electromagnetic radiation spanning the radio to MeV gamma-ray bands provided the first direct evidence that short gamma-ray bursts (GRBs) can originate from binary neutron star (BNS) mergers. On the other hand, the properties of short GRBs in high-energy gamma rays are still poorly constrained, with only ∼20 events detected in the GeV band, and none in the TeV band. GRB 160821B is one of the nearest short GRBs known at z = 0.162. Recent analyses of the multiwavelength observational data of its afterglow emission revealed an optical-infrared kilonova component, characteristic of heavy-element nucleosynthesis in a BNS merger. Aiming to better clarify the nature of short GRBs, this burst was automatically followed up with the MAGIC telescopes, starting from 24 seconds after the burst trigger. Evidence of a gammaray signal is found above ∼0.5 TeV at a significance of ∼ 3 σ during observations that lasted until 4 hours after the burst. Assuming that the observed excess events correspond to gamma-ray emission from GRB 160821B, in conjunction with data at other wavelengths, we investigate its origin in the framework of GRB afterglow models. The simplest interpretation with one-zone models of synchrotronself-Compton emission from the external forward shock has difficulty accounting for the putative TeV flux. Alternative scenarios are discussed where the TeV emission can be relatively enhanced. The role of future GeV-TeV observations of short GRBs in advancing our understanding of BNS mergers and related topics is briefly addressed.

Starburst galaxies and star-forming active galactic nuclei are among the candidate sources thought to contribute appreciably to the extragalactic gamma-ray and neutrino backgrounds. NGC 1068 is the brightest of the star-forming galaxies found to emit gamma-rays from 0.1 to 50 GeV. Precise measurements of the high-energy spectrum are crucial to study the particle accelerators and probe the dominant emission mechanisms. We have carried out 125 hr of observations of NGC 1068 with the MAGIC telescopes in order to search for gamma-ray emission in the veryhigh-energy band. We did not detect significant gamma-ray emission, and set upper limits at the 95% confidence level to the gamma-ray flux above 200 GeV f<5.1×10 −13 cm −2 s −1. This limit improves previous constraints by about an order of magnitude and allows us to put tight constraints on the theoretical models for the gamma-ray emission. By combining the MAGIC observations with the Fermi-LAT spectrum we limit the parameter space (spectral slope, maximum energy) of the cosmic ray protons predicted by hadronuclear models for the gamma-ray emission, while we find that a model postulating leptonic emission from a semi-relativistic jet is fully consistent with the limits. We provide predictions for IceCube detection of the neutrino signal foreseen in the hadronic scenario. We predict a maximal IceCube neutrino event rate of 0.07 yr −1 .