On January 14, 2019, the Major Atmospheric Gamma Imaging Cherenkov telescopes detected GRB 190114C above 0.2 TeV, recording the most energetic photons ever observed from a gamma-ray burst. We use this unique observation to probe an energy dependence of the speed of light in vacuo for photons as predicted by several quantum gravity models. Based on a set of assumptions on the possible intrinsic spectral and temporal evolution, we obtain competitive lower limits on the quadratic leading order of speed of light modification.

Extreme High-frequency BL~Lacs (EHBL) feature their synchrotron peak of the broadband spectral energy distribution (SED) at $\nu_{\rm s} \geq $10$^{17}$\,Hz. The BL~Lac object 1ES~2344+514 was included in the EHBL family because of its impressive shift of the synchrotron peak in 1996. During the following years, the source appeared to be in a low state without showing any extreme behaviours. In August 2016, 1ES~2344+514 was detected with the ground-based $\gamma$-ray telescope FACT during a high $\gamma$-ray state, triggering multi-wavelength (MWL) observations. We studied the MWL light curves of 1ES~2344+514 during the 2016 flaring state, using data from radio to VHE $\gamma$ rays taken with OVRO, KAIT, KVA, NOT, some telescopes of the GASP-WEBT collaboration at the Teide, Crimean, and St. Petersburg observatories, \textit{Swift}-UVOT, \textit{Swift}-XRT, \textit{Fermi}-LAT, FACT and MAGIC. With simultaneous observations of the flare, we built the broadband SED and studied it in the framework of a leptonic and an hadronic model. The VHE $\gamma$-ray observations show a flux level of 55\% of the Crab Nebula flux above 300\,GeV, similar to the historical maximum of 1995. The combination of MAGIC and \textit{Fermi}-LAT spectra provides an unprecedented characterization of the inverse-Compton peak for this object during a flaring episode. The $\Gamma$ index of the intrinsic spectrum in the VHE $\gamma$-ray band is $2.04\pm0.12_{\rm stat}\pm0.15_{\rm sys}$. We find the source in an extreme state with a shift of the position of the synchrotron peak to frequencies above or equal to $10^{18}$\,Hz

Context. It has become evident that one-zone synchrotron self-Compton models are not always adequate for very-high-energy (VHE) gamma-ray emitting blazars. While two-component models are performing better, they are difficult to constrain due to the large number of free parameters. Aims. In this work, we make a first attempt to take into account the observational constraints from Very Long Baseline Interferometry (VLBI) data, long-term light curves (radio, optical, and X-rays) and optical polarisation to limit the parameter space for a two-component model and test if it can still reproduce the observed spectral energy distribution (SED) of the blazars. Methods. We selected five TeV BL Lac objects based on the availability of VHE gamma-ray and optical polarisation data. We collected constraints for the jet parameters from VLBI observations. We evaluated the contributions of the two components to the optical flux by means of decomposition of long-term radio and optical light curves as well as modeling of the optical polarisation variability of the objects. We selected eight epochs for these five objects, based on the variability observed at VHE gamma rays, for which we constructed the SEDs that we then modeled with a two-component model. Results. We found parameter sets which can reproduce the broadband SED of the sources in the framework of two-component models considering all available observational constraints from VLBI observations. Moreover, the constraints obtained from the long-term behavior of the sources in the lower energy bands could be used to determine the region where the emission in each band originates. Finally, we attempted to use optical polarisation data to shed new light on the behavior of the two components in the optical band. Our observationally constrained two zone model allows explanation of the entire SED from radio to VHE with two co-located emission regions.

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 hours 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[J_{\text{ann}}({0.5^{\circ}})/$ GeV$^{2}$ cm$^{-5}] = 19.35 \pm 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 $\tau^{+}\tau^{-}$ final state, where a cross section for annihilation down to $\langle \sigma_{\text{ann}} v \rangle = 3.05 \times 10^{-24}$ cm$^{3}$ s$^{-1}$ is excluded.

We report on a multi-band variability and correlation study of the TeV blazar Mrk 421 during an exceptional flaring activity observed from 2013 April 11 to 2013 April 19. The study uses, among others, data from GASP-WEBT, Swift, NuSTAR, Fermi-LAT, VERITAS, and MAGIC. The large blazar activity, and the 43 hours of simultaneous NuSTAR and MAGIC/VERITAS observations, permitted variability studies on 15 minute time bins, and over three X-ray bands (3-7 keV, 7-30 keV and 30-80 keV) and three very-high-energy (>0.1 TeV, hereafter VHE) gamma-ray bands (0.2-0.4 TeV, 0.4-0.8 TeV and >0.8 TeV). We detected substantial flux variations on multi-hour and sub-hour timescales in all the X-ray and VHE gamma-ray bands. The characteristics of the sub-hour flux variations are essentially energy-independent, while the multi-hour flux variations can have a strong dependence on the energy of the X-ray and the VHE gamma rays. The three VHE bands and the three X-ray bands are positively correlated with no time-lag, but the strength and the characteristics of the correlation changes substantially over time and across energy bands. Our findings favour multi-zone scenarios for explaining the achromatic/chromatic variability of the fast/slow components of the light curves, as well as the changes in the flux-flux correlation on day-long timescales. We interpret these results within a magnetic reconnection scenario, where the multi-hour flux variations are dominated by the combined emission from various plasmoids of different sizes and velocities, while the sub-hour flux variations are dominated by the emission from a single small plasmoid moving across the magnetic reconnection layer.

Imaging atmospheric Cherenkov telescope (IACT) arrays record images from air showers initiated by gamma rays entering the atmosphere, allowing astrophysical sources to be observed at very high energies. To maximize IACT sensitivity, gamma-ray showers must be efficiently distinguished from the dominant background of cosmic-ray showers using images from multiple telescopes. A combination of convolutional neural networks (CNNs) with a recurrent neural network (RNN) has been proposed to perform this task. Using CTLearn, an open source Python package using deep learning to analyze data from IACTs, with simulated data from the upcoming Cherenkov Telescope Array (CTA), we implement a CNN-RNN network and find no evidence that sorting telescope images by total amplitude improves background rejection performance.

Deep convolutional neural networks (DCNs) are a promising machine learning technique to reconstruct events recorded by imaging atmospheric Cherenkov telescopes (IACTs), but require optimization to reach full performance. One of the most pressing challenges is processing raw images captured by cameras made of hexagonal lattices of photo-multipliers, a common layout among IACT cameras which topologically differs from the square lattices conventionally expected, as their input data, by DCN models. Strategies directed to tackle this challenge range from the conversion of the hexagonal lattices onto square lattices by means of oversampling or interpolation to the implementation of hexagonal convolutional kernels. In this contribution we present a comparison of several of those strategies, using DCN models trained on simulated IACT data.

CTLearn is a new Python package under development that uses the deep learning technique to analyze data from imaging atmospheric Cherenkov telescope (IACT) arrays. IACTs use the Cherenkov light emitted from air showers, initiated by very-high-energy gamma rays, to form an image of the longitudinal development of the air shower on the camera plane. The spatial, temporal, and calorimetric information of the originating high-energy particle is then recorded electronically. The sensitivity of IACTs to astrophysical sources depends strongly on the efficient rejection of the background of much more numerous cosmic-ray showers. CTLearn includes modules for running machine learning models with TensorFlow, using pixel-wise camera data as input. Its high-level interface provides a configuration-file-based workflow to drive reproducible training and prediction. We illustrate the capabilities of CTLearn by presenting some results using IACT simulated data.

The search for Dark Matter (DM) has great potential to reveal physics beyond the Standard Model. As such, searches for evidence of DM particles are being carried out using a wide range of techniques, such as direct searches for DM particles, searches for DM produced with colliders, and indirect searches for the Standard Model annihilation products of DM. Dwarf spheroidal galaxies (dSphs) are excellent targets for indirect Dark Matter searches due to their relatively high DM content and negligible expected astrophysical background. A collaboration was formed to maximise the sensitivity of DM searches towards dSphs by combining for the first time dSph data from three imaging air Cherenkov telescope (IACT) arrays: HESS, MAGIC, and VERITAS; the Fermi-LAT satellite, and the water Cherenkov detector HAWC. Due to the diverse nature of the instruments involved, each experiment will analyse their individual datasets from multiple targets and then the results will be combined at the likelihood level. For consistency of the likelihoods across the five experiments, a common approach is used to treat the astrophysical factor (J-Factor) for each target and an agreed set of annihilation channels are considered. We also agree on a common statistical approach and treatment of instrumental systematic uncertainties. The results are presented in terms of constraints on the velocity-weighted cross section for DM self-annihilation as a function of the DM particle mass.

Major construction and initial-phase operation of a second-generation gravitational-wave detector KAGRA has been completed. The entire 3-km detector is installed underground in a mine in order to be isolated from background seismic vibrations on the surface. This allows us to achieve a good sensitivity at low frequencies and high stability of the detector. Bare-bones equipment for the interferometer operation has been installed and the first test run was accomplished in March and April of 2016 with a rather simple configuration. The initial configuration of KAGRA is named {\it iKAGRA}. In this paper, we summarize the construction of KAGRA, including the study of the advantages and challenges of building an underground detector and the operation of the iKAGRA interferometer together with the geophysics interferometer that has been constructed in the same tunnel.