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.

The MAGIC Project: Contributions to ICRC 2007, Merida, Mexico. Contents pages for the Contribution on behalf of the MAGIC Collaboration to the 30th ICRC that took place in July 2007 in Merida, Mexico. The contents are in html form with clickable links to the papers that exist on the Astrophysics archive. We hope that this will make it easier to access the output of the conference in a systematic way. Comments on how useful this is/ how it could be improved should be sent to michela.demaria@iuav.it.

The fast repositioning system of the MAGIC Telescope has allowed during its first data cycle, between 2005 and the beginning of year 2006, observing nine different GRBs as possible sources of very high energy gammas. These observations were triggered by alerts from Swift, HETE-II, and Integral; they started as fast as possible after the alerts and lasted for several minutes, with an energy threshold varying between 80 and 200 GeV, depending upon the zenith angle of the burst. No evidence for gamma signals was found, and upper limits for the flux were derived for all events, using the standard analysis chain of MAGIC. For the bursts with measured redshift, the upper limits are compatible with a power law extrapolation, when the intrinsic fluxes are evaluated taking into account the attenuation due to the scattering in the Metagalactic Radiation Field (MRF).

During the last decade ground-based very high-energy gamma-ray astronomy achieved a remarkable advancement in the development of the observational technique for the registration and study of gamma-ray emission above 100 GeV. It is widely believed that the next step in its future development will be the construction of telescopes of substantially larger size than the currently used 10 m class telescopes. This can drastically improve the sensitivity of the ground-based detectors for gamma rays of energy from 10 to 100 GeV. Based on Monte Carlo simulations of the response of a single stand-alone 30 m imaging atmospheric Cherenkov telescope (IACT) the maximal rejection power against background cosmic ray showers for low energy gamma-rays was investigated in great detail. An advanced Bayesian multivariate analysis has been applied to the simulated Cherenkov light images of the gamma-ray- and proton-induced air showers. The results obtained here quantitatively testify that the separation between the signal and background images degrades substantially at low energies, and consequently the maximum overall quality factor can only be about 3.1 for gamma rays in the 10-30 GeV energy range. Various selection criteria as well as optimal combinations of the standard image parameters utilized for effective image separation have been also evaluated.

The MAGIC collaboration has studied the high peaked BL-Lac object 1ES1218+30.4 at a redshift z = 0.182, using the MAGIC imaging air Cherenkov telescope located on the Canary island of La Palma. A gamma-ray signal was observed with 6.4sigma significance. The differential energy spectrum for an energy threshold of 120GeV can be fitted by a simple power law yielding F_E(E) = (8.1+-2.1)*10^-7 (E/250GeV)^(-3.0+-0.4) TeV^-1 m^-2 s^-1. During the six days of observation in January 2005 no time variability on time scales of days was found within the statistical errors. The observed integral flux above 350GeV is nearly a factor two below the the upper limit reported by the Whipple Collaboration in 2003.

The long-duration GRB050713a was observed by the MAGIC Telescope, 40 seconds after the burst onset, and followed up for 37 minutes, until twilight. The observation, triggered by a SWIFT alert, covered energies above ~175 GeV. Using standard MAGIC analysis, no evidence for a gamma signal was found. As the redshift of the GRB was not measured directly, the flux upper limit, estimated by MAGIC, is still compatible with the assumption of an unbroken power-law spectrum extending from a few hundred keV to our energy range.

We are investigated possible correlations between the calculated arrival times of the first relativistic ions at Earth and GLE start times registered by surface monitors. The analysis is based on the arrival times and energies of the first solar ions, registered by the Solar Isotope Spectrometer (SIS) on board of the ACE satellite, and protons, registered by GOES satellites. We consider both cases when the interplanetary propagation of the first high energy ions is essentially scatter-free and the diffusion of high energy ions during propagation in the interplanetary magnetic field. We extrapolate the time-velocity and time-rigidity relationships to calculate the expected arrival times of the relativistic ions that are energetic enough to enter the atmosphere at the Aragats geographical location and produce secondary fluxes that reach the monitors.

Using the EAS size spectra measured with the MAKET ANI array on Mt. Aragats, Armenia (3200m a.s.l.) in the range of the shower size 10^5 - 10^7 for different angles-of-incidence, the EAS attenuation length has been determined applying different analysis methods. The analysis is based on a data sample of 2.5*10^6 events collected in the period of June, 97 - April, 99. The results are compared with results deduced from data of the EAS TOP and KASCADE experiments.

The quality of the MAKET-ANI detector installation in view of the uniformity of the registration efficiency is demonstrated. Based on a data sample collected by the MAKET-ANI array in the period of June 1997 - March 1999, the dependencies of the age parameter on the zenith angle and the EAS size (10^5-10^7) are studied. The variation of the age parameter with the shower size can be approximately related to the elongation rate.