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.

We examine a model for the observed temporal variability of powerful blazars in the $\gamma$-ray band in which the dynamics is described in terms of a stochastic differential equation, including the contribution of a deterministic drift and a stochastic term. The form of the equation is motivated by the current astrophysical framework, accepting that jets are powered through the extraction of the rotational energy of the central supermassive black hole mediated by magnetic fields supported by a so-called \emph{magnetically arrested} accretion disk. We apply the model to the $\gamma$-ray light curves of several bright blazars and we infer the parameters suitable to describe them. In particular, we examine the differential distribution of fluxes ($dN/dF_{\gamma}$) and we show that the predicted probability density function for the assumed stochastic equation naturally reproduces the observed power law shape at large fluxes $dN/dF_{\gamma} \propto F_{\gamma}^{-\alpha}$ with $\alpha>2$.

We investigate with relativistic MHD simulations the dissipation physics of BL Lac jets, by studying the synchrotron polarization signatures of particles accelerated by the kink instability in a magnetically-dominated plasma column. The nonlinear stage of the kink instability generates current sheets, where particles can be efficiently accelerated via magnetic reconnection. We identify current sheets as regions where s = J d/B is above some predefined threshold (where B is the field strength, J the current density and d the grid scale), and assume that the particle injection efficiency scales as proportional to the square of the current. X-ray emitting particles have short cooling times, so they only probe the field geometry of their injection sites. In contrast, particles emitting in the optical band, which we follow self-consistently as they propagate away from their injection sites while cooling, sample a larger volume, and so they may be expected to produce different polarimetric signatures. We find that the degree of polarization is roughly the same between X-ray and optical bands, because even the optical-emitting particles do not travel far from the current sheet where they were injected, due to lack of sufficient kink-generated turbulence. The polarization angle shows a different temporal evolution between the two bands, due to the different regions probed by X-ray and optical emitting particles. In view of the upcoming IXPE satellite, our results can help constrain whether kink-induced reconnection (as opposed to shocks) can be the source of multi-wavelength emission from BL Lacs.

Flat Spectrum Radio Quasars (FSRQ) are the most powerful blazars in the gamma-ray band. Although they are supposed to be good candidates in producing high energy neutrinos, no secure detection of FSRQs has been obtained up to now, except for a possible case of PKS B1424-418. In this work, we compute the expected flux of high energy neutrinos from FSRQs using standard assumptions for the properties of the radiation fields filling the regions surrounding the central supermassive black hole. We obtain as a result that high energy neutrinos are naturally expected from FSRQs in the sub-EeV-EeV energy range and not at PeV energies. This justifies the non-observation of neutrinos from FSRQs with the present technology, since only neutrinos below 10 PeV have been observed. We found that for a non-negligible range of the parameters the cumulative flux from FSRQs is comparable to or even exceeds the expected cosmogenic neutrino flux. This result is intriguing and highlights the importance to disentangle these point-source emissions from the diffuse cosmogenic background.

Long-duration gamma-ray bursts (GRBs) originate from ultra-relativistic jets launched from the collapsing cores of dying massive stars. They are characterised by an initial phase of bright and highly variable radiation in the keV-MeV band that is likely produced within the jet and lasts from milliseconds to minutes, known as the prompt emission. Subsequently, the interaction of the jet with the external medium generates external shock waves, responsible for the afterglow emission, which lasts from days to months, and occurs over a broad energy range, from the radio to the GeV bands. The afterglow emission is generally well explained as synchrotron radiation by electrons accelerated at the external shock. Recently, an intense, long-lasting emission between 0.2 and 1 TeV was observed from the GRB 190114C. Here we present the results of our multi-frequency observational campaign of GRB~190114C, and study the evolution in time of the GRB emission across 17 orders of magnitude in energy, from $5\times10^{-6}$ up to $10^{12}$\,eV. We find that the broadband spectral energy distribution is double-peaked, with the TeV emission constituting a distinct spectral component that has power comparable to the synchrotron component. This component is associated with the afterglow, and is satisfactorily explained by inverse Compton upscattering of synchrotron photons by high-energy electrons. We find that the conditions required to account for the observed TeV component are not atypical, supporting the possibility that inverse Compton emission is commonly produced in GRBs.

Gamma-ray bursts (GRBs) of the long-duration class are the most luminous sources of electromagnetic radiation known in the Universe. They are generated by outflows of plasma ejected at near the speed of light by newly formed neutron stars or black holes of stellar mass at cosmological distances. Prompt flashes of MeV gamma rays are followed by longer-lasting afterglow emission from radio waves to GeV gamma rays, due to synchrotron radiation by energetic electrons in accompanying shock waves. Although emission of gamma rays at even higher, TeV energies by other radiation mechanisms had been theoretically predicted, it had never been detected previously. Here we report the clear detection of GRB 190114C in the TeV band, achieved after many years of dedicated searches for TeV emission from GRBs. Gamma rays in the energy range 0.2--1 TeV are observed from about 1 minute after the burst (at more than 50 standard deviations in the first 20 minutes). This unambiguously reveals a new emission component in the afterglow of a GRB, whose power is comparable to that of the synchrotron component. The observed similarity in the radiated power and temporal behaviour of the TeV and X-ray bands points to processes such as inverse Compton radiation as the mechanism of the TeV emission, while processes such as synchrotron emission by ultrahigh-energy protons are disfavoured due to their low radiative efficiency.

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.

Aims: $\gamma$ rays can be used as a tracer in the search of sources of Galactic cosmic rays (CRs). We present deep observations of the Galactic Centre (GC) region with the MAGIC telescopes, which we use for inferring the underlying CR distribution. Methods: We observed the GC region for ${\approx}100$ hours with the MAGIC telescopes from 2012 to 2017, at high zenith angles (58-70~deg). This implies a larger energy threshold, but also an increased effective collection area compared to low zenith observations. Using new software, we derive instrument response and background models, enabling us to study the diffuse emission in the region. We use pre-existing data of the gas distribution in the GC region to derive the underlying distribution of CRs. 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 find that the diffuse component is best 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~$\sigma$. The derived cosmic-ray profile hints to a peak at the GC position, with a measured profile index of $1.2 \pm 0.3$, supporting 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.

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.

Within the classical emission model, where the emission region is placed within the broad line region (BLR), flat spectrum radio quasars (FSRQs) were believed not to emit photons with energies above few tens of GeV because of the absorption with the optical-UV photons from the BLR. However, photons with observed energies up to about $300 \, \rm GeV$ have been detected for few FSRQs, whose most iconic example is PKS 1441+25 at redshift $z = 0.94$. The most conservative explanation for these observations is that the emission occurs at distances comparable to the size of the dusty torus. In this case, absorption of high-energy gamma-ray photons for energies above $200-300 \, {\rm GeV}$ is dominated by the interaction with infrared radiation emitted by the torus. We investigate if current observational data about FSRQs in flaring state can give us information about: (i) the importance of the torus absorption and (ii) the properties of the torus i.e. its temperature and its geometry. We find that present data do not arrive at energies where the torus influence is prominent and as a result it is currently hardly possible to infer torus properties from observations. However, with dedicated simulations, we demonstrate that observations with the forthcoming Cherenkov Telescope Array (CTA) will be able to constrain the torus parameters (temperature and geometry).

We report on variability and correlation studies using multiwavelength observations of the blazar Mrk 421 during the month of February, 2010 when an extraordinary flare reaching a level of $\sim$27~Crab Units above 1~TeV was measured in very-high-energy (VHE) $\gamma$-rays with the VERITAS observatory. This is the highest flux state for Mrk 421 ever observed in VHE $\gamma$-rays. Data are analyzed from a coordinated campaign across multiple instruments including VHE $\gamma$-ray (VERITAS, MAGIC), high-energy (HE) $\gamma$-ray (Fermi-LAT), X-ray (Swift}, RXTE, MAXI), optical (including the GASP-WEBT collaboration and polarization data) and radio (Mets\"ahovi, OVRO, UMRAO). Light curves are produced spanning multiple days before and after the peak of the VHE flare, including over several flare `decline' epochs. The main flare statistics allow 2-minute time bins to be constructed in both the VHE and optical bands enabling a cross-correlation analysis that shows evidence for an optical lag of $\sim$25-55 minutes, the first time-lagged correlation between these bands reported on such short timescales. Limits on the Doppler factor ($\delta \gtrsim 33$) and the size of the emission region ($ \delta^{-1}R_B \lesssim 3.8\times 10^{13}\,\,\mbox{cm}$) are obtained from the fast variability observed by VERITAS during the main flare. Analysis of 10-minute-binned VHE and X-ray data over the decline epochs shows an extraordinary range of behavior in the flux-flux relationship: from linear to quadratic to lack of correlation to anti-correlation. Taken together, these detailed observations of an unprecedented flare seen in Mrk 421 are difficult to explain by the classic single-zone synchrotron self-Compton model.

1ES 1959+650 is a bright TeV high-frequency-peaked BL Lac object exhibiting interesting features like "orphan" TeV flares and a 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) gamma-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 13th, 14th of June and 1st 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 gamma-ray (both Fermi-LAT and VHE bands) and the X-ray bands are quite hard. On 13th June and 1st July 2016, the source showed rapid variations of 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 high 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 L_{jet}~10^{46} erg/s 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.

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) gamma-rays from the Crab pulsar. 16 hours 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 with the Major Atmospheric Gamma-ray Imaging Cherenkov (MAGIC) telescopes were performed. We searched for a statistical correlation between the radio and VHE gamma-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 gamma-rays correlated with the occurrence of radio GPs was carried out separately for the P1 and P2 phase ranges respectively. 99444 radio GPs have been detected in the radio data sample. 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 95% confidence level on an increase in VHE gamma-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 of 12% to 2900% (depending on search window duration and phase cuts) of the pulsed VHE flux. This is the most stringent upper limit on a correlation between gamma-ray emission and radio GPs reported so far.}

M87 is one of the closest (z=0.00436) extragalactic sources emitting at very-high-energies (VHE, E > 100 GeV). The aim of this work is to locate the region of the VHE gamma-ray emission and to describe the observed broadband spectral energy distribution (SED) during the low VHE gamma-ray state. The data from M87 collected between 2012 and 2015 as part of a MAGIC monitoring programme are analysed and combined with multi-wavelength data from Fermi-LAT, Chandra, HST, EVN, VLBA and the Liverpool Telescope. The averaged VHE gamma-ray spectrum can be fitted from 100GeV to 10TeV with a simple power law with a photon index of (-2.41 $\pm$ 0.07), while the integral flux above 300GeV is $(1.44 \pm 0.13) \times 10^{-12} cm^{-2} s^{-1}$. During the campaign between 2012 and 2015, M87 is generally found in a low emission state at all observed wavelengths. The VHE gamma-ray flux from the present 2012-2015 M87 campaign is consistent with a constant flux with some hint of variability ($\sim3\sigma$) on a daily timescale in 2013. The low-state gamma-ray emission likely originates from the same region as the flare-state emission. Given the broadband SED, both a leptonic synchrotron self Compton and a hybrid photo-hadronic model reproduce the available data well, even if the latter is preferred. We note, however, that the energy stored in the magnetic field in the leptonic scenario is very low suggesting a matter dominated emission region.

Aims: We aim to measure the Crab Nebula gamma-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 use the novel very large zenith angle observations with the MAGIC telescope system to increase the collection area above 10 TeV. We also develop an auxiliary procedure of monitoring atmospheric transmission in order to assure proper calibration of the accumulated data. This employs recording of 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 hr of observations, we detect the gamma-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 gamma-ray signal.

The most powerful persistent accelerators in the Universe are jetted active galaxies. Galaxies whose jets are directed towards Earth, so called blazars, dominate the extragalactic gamma-ray sky. Still, most of the highest-energy accelerators likely elude detection. These extreme blazars, whose radiated energy can peak beyond 10 TeV, are ideal targets to study particle acceleration and radiative processes, and may provide links to cosmic rays and astrophysical neutrinos. The growing number of extreme blazars observed at TeV energies has been critical for the emergence of gamma-ray cosmology, including measurements of the extragalactic background light, tight bounds on the intergalactic magnetic field, and constraints on exotic physics at energies inaccessible with human-made accelerators. Tremendous progress is expected in the decade to come, particularly with the deployment of the Cherenkov Telescope Array.

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.

Mrk501 is a very high-energy (VHE) gamma-ray blazar located at z=0.034. During a period of two weeks in July 2014, the highest X-ray activity of Mrk501 was observed in ~14 years of operation of the Neil Gehrels Swift Gamma-ray Burst Observatory. We characterize the broadband variability of Mrk501 from radio to VHE gamma rays, and evaluate whether it can be interpreted within theoretical scenarios widely used to explain the broadband emission from blazars. The temporal evolution of the most prominent and variable segments of the SED 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. 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$\sigma$. A fit with a log-parabola plus a narrow component is preferred over the fit with a single log-parabola at more than 4$\sigma$, and a dedicated Monte Carlo simulation estimated the significance of this extra component to be larger than 3$\sigma$. 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.

Surprisingly, the contribution to the cosmic opacity to UHE ($E>10^{18} \, \rm eV$) $\gamma$-rays has been systematically computed so far for the $\gamma \gamma \to e^+e^-$ and $\gamma \gamma \to e^+ e^- e^+ e^-$ processes alone. We go a step further by systematically evaluating the additional opacity brought about by other leptons and hadrons. We find that the dominant channels are those leading to the production of $\mu^{\pm}$ and hadrons (mainly $\pi^{\pm}$, $\pi^0$, $K^{\pm}$, $K^0$, $\eta$). For nearly the GZK radius, the photon survival probability becomes smaller by about a factor of two with respect to current estimates.

Very-high-energy (VHE) BL Lac spectra extending above $10 \, \rm TeV$ provide a unique opportunity for testing physics beyond the standard model of elementary particle and alternative blazar emission models. We consider the hadron beam, the photon to axion-like particle (ALP) conversion, and the Lorentz invariance violation (LIV) scenarios by analyzing their consequences and induced modifications to BL Lac spectra. In particular, we consider how different processes can provide similar spectral features (e.g. hard tails) and we discuss the ways they can be disentangled. We use HEGRA data of a high state of Markarian 501 and the HESS spectrum of the extreme BL Lac (EHBL) 1ES 0229+200. In addition, we consider two hypothetical EHBLs similar to 1ES 0229+200 located at redshifts $z=0.3$ and $z=0.5$. We observe that both the hadron beam and the photon-ALP oscillations predict a hard tail extending to energies larger than those possible in the standard scenario. Photon-ALP interaction predicts a peak in the spectra of distant BL Lacs at about $20-30 \, \rm TeV$, while LIV produces a strong peak in all BL Lac spectra around $\sim 100 \, \rm TeV$. The peculiar feature of the photon-ALP conversion model is the production of oscillations in the spectral energy distribution, so that its detection/absence can be exploited to distinguish among the considered models. The above mentioned features coming from the three models may be detected by the upcoming Cherenkov Telescope Array (CTA). Thus, future observations of BL Lac spectra could eventually shed light about new physics and alternative blazar emission models, driving fundamental research towards a specific direction.

Extreme high-energy peaked BL Lac objects (EHBLs) are an emerging class of blazars. Their typical two-hump structured spectral energy distribution (SED) peaks at higher energies with respect to conventional blazars. Multi-wavelength (MWL) observations constrain their synchrotron peak in the medium to hard X-ray band. Their gamma-ray SED peaks above the GeV band, and in some objects it extends up to several TeV. Up to now, only a few EHBLs have been detected in the TeV gamma-ray range. In this paper, we report the detection of the EHBL 2WHSP J073326.7+515354, observed and detected during 2018 in TeV gamma rays with the MAGIC telescopes. The broad-band SED is studied within a MWL context, including an analysis of the Fermi-LAT data over ten years of observation and with simultaneous Swift-XRT, Swift-UVOT, and KVA data. Our analysis results in a set of spectral parameters that confirms the classification of the source as an EHBL. In order to investigate the physical nature of this extreme emission, different theoretical frameworks were tested to model the broad-band SED. The hard TeV spectrum of 2WHSP J073326.7+515354 sets the SED far from the energy equipartition regime in the standard one-zone leptonic scenario of blazar emission. Conversely, more complex models of the jet, represented by either a two-zone spine-layer model or a hadronic emission model, better represent the broad-band SED.

Extreme high-frequency peaked BL Lac objects (EHBLs) are blazars which exhibit extremely energetic synchrotron emission. They also feature non-thermal gamma-ray emission whose peak lies in the very high-energy (VHE, E > 100 GeV) range, and in some sources exceeds 1TeV: this is the case of hard-TeV EHBLs such as 1ES 0229+200. With the aim of increasing the EHBL population, ten targets were observed with the MAGIC telescopes from 2010 to 2017, for a total of 262 h of good quality data. The data were complemented by coordinated Swift observations. The X-ray data analysis confirms that all the sources but two 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 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 and a proton synchrotron models. 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 a first mini-catalog of VHE gamma-ray and multi-wavelength observations of EHBLs.

We investigate the impact on the properties of high-energy emitting BL Lac objects of a population of electrons with an anisotropic momentum distribution. We adopt a simple phenomenological description of the anisotropy, in which the most energetic electrons have a small pitch angle and the least energetic electrons are isotropic, as proposed by Sobacchi \& Lyubarsky (2019). We explore (i) a simple model that assumes a phenomenological shape for the electron energy distribution, and (ii) a self-consistent scheme in which the electrons follow a distribution which is the result of the balance between injection and radiative losses (we include the effects of the anisotropy on the synchrotron cooling rate). Considering the BL Lac object Mkn 421 as representative of the entire class, we show that in both cases the emission can be satisfactorily reproduced under equipartition between the magnetic field and the relativistic electrons. This is in better agreement with the idea that jets are launched as Poynting dominated flows with respect to the standard isotropic scenario, which requires both a low magnetization and a low radiative efficiency to reproduce the observed emission. The hard spectrum predicted for the inverse Compton continuum at TeV energies could be used as a potential test of the anisotropic model.

The gamma-ray blazar TXS 0506+056, was discovered in VHE gamma-rays by the MAGIC telescopes in 2017 in a follow-up campaign of a high energy neutrino event IceCube-170922A (IC+Fermi+MAGIC++, Science 361, eaat1378 (2018)). Subsequent multivawelenght (MWL) observations and theoretical modeling in a frame of hadro-leptonic emission confirmed that this source could be a potential cosmic ray and neutrino emitter (MAGIC Collaboration, Ansoldi et al., (2018)). This is, by far, the most significant association between a high-energy neutrino and an astrophysical source emitting gamma rays and X-rays. TXS 0506+056 is a key object to help the astrophysics community to establish connections between high-energy neutrinos and astrophysical sources. Accurate and contemporaneous MWL spectral measurements are essential ingredients to achieve this goal. In the conference, we present the measurements from the MAGIC and MWL monitoring of this source, spanning the time period from November 2017 till February 2019. These include the lowest VHE gamma-ray emission state measured from this source so far as well as a flaring episode in December 2018.

Extreme high-energy peaked BL Lac objects (EHBLs) are blazars whose synchrotron emission peaks at exceptionally high energies, above few keV, in the hard X-ray regime. So far, only a handful of those objects has been detected at very high energy (VHE, E > 100 GeV) gamma rays by Imaging Atmospheric Cherenkov Telescopes. Very remarkably, VHE observations of some of these blazars (like 1ES 0229+200) have provided evidence of a VHE gamma-ray emission extending to several TeV, which is difficult to explain with standard, one-zone synchrotron self-Compton models usually applied to BL Lac objects. The MAGIC collaboration coordinated a multi-year, multi-wavelength observational campaign on ten targets. The MAGIC telescopes detected VHE gamma rays from four EHBLs. In this paper we focus on the source 1ES 1426+426 and its X-ray and VHE gamma-ray properties. The results of different models (synchrotron self-Compton, spine-layer, hadronic) reproducing the broadband spectral energy distribution are also presented.

Extreme high-energy peaked BL Lac objects (EHBLs) are a new emerging class of blazars. The typical two-hump structured spectral energy distribution (SED) is shifted to higher energies with respect to other more established classes of blazars. Multi-wavelength observations allow us to constrain their synchrotron peak in the medium and hard X-ray bands. Their gamma-ray emission dominates above the GeV gamma-ray band, and in some objects it extends up to several TeV (e.g. 1ES 0229+200). Their hard TeV spectrum is also interesting for the implications on the extragalactic background light indirect measurements, the intergalactic magnetic field estimate, and the possible origin of extragalactic high-energy neutrinos. Up to now, only a few objects have been studied in the TeV gamma-ray range. In this contribution, we will present the new detection of the EHBL object PGC 2402248, recently discovered in TeV gamma rays with the MAGIC telescopes. The analysis results of a set of multi-wavelength simultaneous observations up to the VHE gamma-ray band provide the broad-band SED of the blazar, which will be used to probe different emission models. Given the extreme characteristics of this blazar, constraints on the physical parameters within the framework of leptonic and hadronic models are derived.

Starburst galaxies and star-forming active galactic nuclei (AGN) 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 hours of observations of NGC 1068 with the MAGIC telescopes in order to search for gamma-ray emission in the very high energy band. We did not detect significant gamma-ray emission, and set upper limits at 95\% confidence level to the gamma-ray flux above 200 GeV f<5.1x10^{-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}.

Context. Measurements of the Extragalactic Background Light (EBL) are a fundamental source of information on the collective emission of cosmic sources. Aims. At infrared wavelengths, however, these measurements are precluded by the overwhelming dominance from Interplanetary Dust emission and the Galactic infrared foreground. Only at $\lambda > 300 \ \mu$m, where the foregrounds are minimal, has the Infrared EBL (IR EBL) been inferred from analysis of the COBE maps. The present paper aims to assess the possibility of evaluating the IR EBL from a few $\mu$m up to the peak of the emission at >100 $\mu$m using an indirect method that avoids the foreground problem. Methods. To this purpose we exploit the effect of pair-production from gamma-gamma interaction by considering the highest energy photons emitted by extragalactic sources and their interaction with the IR EBL photons. We simulate observations of a variety of low redshift emitters with the forthcoming Imaging Atmospheric Cherenkov Telescope (IACT) arrays (CTA in particular) and water Cherenkov observatories (LHAASO, HAWC, SWGO) to assess their suitability to constrain the EBL at such long wavelengths. Results. We find that, even under the most extremely favorable conditions of huge emission flares, extremely high-energy emitting blazars are not very useful for our purpose because they are much too distant (>100 Mpc the nearest ones, MKN 501 and MKN 421). Observations of more local Very High Energy (VHE) emitting AGNs, like low-redshift radio galaxies (M87, IC 310, Centaurus A), are better suited and will potentially allow us to constrain the EBL up to $\lambda \simeq 100\ \mu$m.

We examine the potential multi-messenger role of BL Lac objects emitting at high energy (so-called HBL) focusing on the limits on the energy density of soft radiation in the jet frame, a critical parameter which regulates the proton cooling and the fragmentation of ultra-high energy cosmic ray (UHECR) nuclei possibly accelerated in the jet. We show that (under the assumption that the high-energy emission bump is dominated by inverse Compton emission) the energy density of any external soft radiation field (e.g. produced by a layer surrounding the jet or in the accretion flow) cannot be larger than few times that associated to the observed synchrotron radiation produced in the emission region. Quite interestingly, the constraint that we derive is generally stronger than the limit obtained from the condition that the source is transparent to very-high energy $\gamma$-rays. Using this constraint we can derive a robust upper limit for the efficiency of the photopion reaction leading to the emission of PeV neutrinos, $f_{\pi}\lesssim 10^{-5}$, which makes HBL quite inefficient neutrino sources. For the photodisintegration of nuclei the results are more dependent on the spectral properties of the radiation field. The photodisintegration efficiency is safely below 1 (and nuclei can escape intact) for a `canonical' spectrum of the soft radiation field $\propto \nu^{-0.5}$. For radiation fields characterized by a softer spectrum and extended over a large portion of the jet, the efficiency increases and for an appreciable fraction of the sources nuclei with energies above $10^{19}$ eV might suffer significant photodisintegration.

We observed with the NuSTAR satellite 3 blazars at z>2, detected in gamma-rays by Fermi/LAT and in the soft X-rays, but not yet observed above 10 keV. The flux and slope of their X-ray continuum, together with Fermi/LAT data allows us to estimate their total electromagnetic output and peak frequency. For some of them we can study the source in different states, and investigate the main cause of the observed different spectral energy distribution. We then collected all blazars at redshift greater than 2 observed by NuSTAR, and confirm that these hard and luminous X-ray blazars are among the most powerful persistent sources in the Universe. We confirm the relation between the jet power and the disk luminosity, extending it at the high energy end.

Prompted by the increasing interest of axion-like particles (ALPs) for very-high-energy (VHE) astrophysics, we have considered a full scenario for the propagation of a VHE photon/ALP beam emitted by a BL Lac and reaching us in the light of the most up-to-date astrophysical information and for energies up to above $100 \, \rm TeV$. During its trip, the beam -- generated in a small region of a BL Lac jet -- crosses a variety of magnetic structures in very different astronomical environments: the BL Lac jet, the host elliptical galaxy, the extragalactic space and the Milky Way. We have taken an effort to model all these magnetic fields in the most realistic fashion and using a new model developed by us concerning the extragalactic magnetic field. Assuming an intrinsic spectrum with a power law exponentially truncated at a fixed cut-off energy, we have evaluated the resulting observed spectra of Markarian 501, the extreme BL Lac 1ES 0229+200 and a similar source located at $z = 0.6$ up to above $100 \, \rm TeV$. We obtain interesting results: the model with photon-ALP oscillations possesses features (spectral energy oscillatory behaviour and photon excess above $20 \, \rm TeV$) which can be tested by $\gamma$-ray observatories like CTA, HAWC, GAMMA 400, LHAASO, TAIGA-HiSCORE and HERD. In addition, our ALP can be detected in dedicated laboratory experiments like the upgrade of ALPS II at DESY, the planned IAXO and STAX experiments, as well as with other techniques developed by Avignone and collaborators.

Several observations are revealing the widespread occurrence of mildly relativistic wide-angle AGN winds strongly interacting with the gas of their host galaxy. Such winds are potential cosmic-ray accelerators, as supported by gamma-ray observations of the nearby Seyfert galaxy NGC 1068 with the Fermi gamma-ray space telescope. The non-thermal emission produced by relativistic particles accelerated by the AGN-driven wind observed in the circum-nuclear molecular disk of such galaxy is invoked to produce the gamma-ray spectrum. The AGN wind model predicts a hard spectrum that extend in the very high energy band which differs significantly from those corresponding to other models discussed in the literature, like starburst or AGN jet. We present dedicated simulations of observations through the Cherenkov Telescope Array (CTA), the next-generation ground based gamma-ray observatory, of the very high energy spectrum of the Seyfert galaxy NGC 1068 assuming the AGN wind and the AGN jet models. We demonstrate that, considering 50 hours of observations, CTA can be effectively used to constrain the two different emission models, providing important insight into the physics governing the acceleration of particles in non-relativistic AGN-driven outflows. This analysis strongly motivates observations of Seyfert and starburst galaxies with CTA in order to test source population models of the extragalactic gamma-ray and neutrino backgrounds.

The Cherenkov Telescope Array (CTA) is the major next-generation observatory for ground-based very-high-energy gamma-ray astronomy. It will improve the sensitivity of current ground-based instruments by a factor of five to twenty, depending on the energy, greatly improving both their angular and energy resolutions over four decades in energy (from 20 GeV to 300 TeV). This achievement will be possible by using tens of imaging Cherenkov telescopes of three successive sizes. They will be arranged into two arrays, one per hemisphere, located on the La Palma island (Spain) and in Paranal (Chile). We present here the optimised and final telescope arrays for both CTA sites, as well as their foreseen performance, resulting from the analysis of three different large-scale Monte Carlo productions.

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 to 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 photon-photon 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-sigma 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 $\mu$m range relative to the studied models, yet compatible with them within systematics.

A population of globular clusters (GCs) has been recently established by the Fermi-LAT telescope as a new class of GeV $\gamma$-ray sources. Leptons accelerated to TeV energies, in the inner magnetospheres of MSPs or in their wind regions, should produce $\gamma$-rays through the inverse Compton scattering in the dense radiation field from the huge population of stars. We have conducted deep observations of the globular cluster M15 with the MAGIC telescopes and used 165 hrs in order to search for $\gamma$-ray emission. A strong upper limit on the TeV $\gamma$-ray flux $<3.2\times 10^{-13}\mathrm{cm^{-2}s^{-1}}$ above 300 GeV ($<0.26\%$ of the Crab nebula flux) has been obtained. We interpret this limit as a constraint on the efficiency of the acceleration of leptons in the magnetospheres of the MSPs. We constrain the injection rate of relativistic leptons, $\eta_{\rm e}$, from the MSPs magnetospheres and their surrounding. We conclude that $\eta_{\rm e}$ must be lower than expected from the modelling of high energy processes in MSP inner magnetospheres. For leptons accelerated with the power law spectrum in the MSP wind regions, $\eta_{\rm e}$ is constrained to be much lower than derived for the wind regions around classical pulsars. These constraints are valid for the expected range of magnetic field strengths within the GC and for the range of likely energies of leptons injected from the inner magnetospheres, provided that the leptons are not removed from the globular cluster very efficiently due to advection process. We discuss consequences of these constraints for the models of radiation processes around millisecond pulsars.

The HAWC Collaboration released the 2HWC catalog of TeV sources, in which 19 show no association with any known high-energy (HE; E > 10 GeV) or very-high-energy (VHE; E > 300 GeV) sources. This catalog motivated follow-up studies by both the MAGIC and Fermi-LAT observatories with the aim of investigating gamma-ray emission over a broad energy band. In this paper, we report the results from the first joint work between HAWC, MAGIC and Fermi-LAT on three unassociated HAWC sources: 2HWC J2006+341, 2HWC J1907+084* and 2HWC J1852+013*. Although no significant detection was found in the HE and VHE regimes, this investigation shows that a minimum 1 degree extension (at 95% confidence level) and harder spectrum in the GeV than the one extrapolated from HAWC results are required in the case of 2HWC J1852+013*, while a simply minimum extension of 0.16 degrees (at 95% confidence level) can already explain the scenario proposed by HAWC for the remaining sources. Moreover, the hypothesis that these sources are pulsar wind nebulae is also investigated in detail.

The mechanisms producing fast variability of the $\gamma$-ray emission in active galactic nuclei are under debate. The MAGIC telescopes detected a fast very high energy (VHE, E$>100$ GeV) $\gamma$-ray flare from BL Lacertae on 2015 June 15. The flare had a maximum flux of $(1.5\pm 0.3)\times 10^{-10}$ photons cm$^{-2}$ s$^{-1}$ and halving time of $26\pm8$ minutes. The MAGIC observations were triggered by a high state in the optical and high energy (HE, E$>100$ MeV) $\gamma$-ray bands. In this paper we present the MAGIC VHE $\gamma$-ray data together with multiwavelength data from radio, optical, X-rays, and HE $\gamma$ rays from 2015 May 1 to July 31. Well-sampled multiwavelength data allow us to study the variability in detail and compare it to the other epochs when fast VHE $\gamma$-ray flares have been detected from this source. Interestingly, we find that the behaviour in radio, optical, X-rays and HE $\gamma$-rays is very similar to two other observed VHE $\gamma$-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 $\gamma$-ray data with the star-jet interaction model. We find that all of the tested emission models are compatible with the fast VHE $\gamma$-ray flare, but all have some tension with the multiwavelength observations.

Simulating and analysing detailed observations of astrophysical sources for very high energy (VHE) experiments, like the Cherenkov Telescope Array (CTA), can be a demanding task especially in terms of CPU consumption and required storage. In this context, we propose an innovative cloud computing architecture based on Amazon Web Services (AWS) aiming to decrease the amount of time required to simulate and analyse a given field by distributing the workload and exploiting the large computational power offered by AWS. We detail how the various services offered by the Amazon online platform are jointly used in our architecture and we report a comparison of the execution times required for simulating observations of a test source with the CTA, by a single machine and the cloud-based approach. We find that, by using AWS, we can run our simulations more than 2 orders of magnitude faster than by using a general purpose workstation for the same cost. We suggest to consider this method when observations need to be simulated, analysed, and concluded within short timescales.

SNR G24.7+0.6 is a 9.5 kyrs radio and $\gamma$-ray supernova remnant evolving in a dense medium. In the GeV regime, SNR G24.7+0.6 (3FHL\,J1834.1--0706e/FGES\,J1834.1--0706) shows a hard spectral index ($\Gamma$$\sim$2) up to $200$\,GeV, which makes it a good candidate to be observed with Cherenkov telescopes such as MAGIC. We observed the field of view of \snr\ with the MAGIC telescopes for a total of 31 hours. We detect very high energy $\gamma$-ray emission from an extended source located 0.34\degr\ away from the center of the radio SNR. The new source, named \mgc\ is detected up to 5\,TeV, and its spectrum is well-represented by a power-law function with spectral index of $2.74 \pm 0.08$. The complexity of the region makes the identification of the origin of the very-high energy emission difficult, however the spectral agreement with the LAT source and overlapping position at less than 1.5$\sigma$ point to a common origin. We analysed 8 years of \fermi-LAT data to extend the spectrum of the source down to 60\,MeV. \fermi-LAT and MAGIC spectra overlap within errors and the global broad band spectrum is described by a power-law with exponential cutoff at $1.9\pm0.5$\,TeV. The detected $\gamma$-ray emission can be interpreted as the results of proton-proton interaction between the supernova and the CO-rich surrounding.

The aim of this white paper is to discuss the observing strategies for the LSST Wide-Fast-Deep that would improve the study of blazars (emission variability, census, environment) and Fast Radio Bursts (FRBs). For blazars, these include the adoption of: i) a reference filter to allow reconstruction of a well-sampled light curve not affected by colour changes effects; ii) two snapshots/visit with different exposure times to avoid saturation during flaring states; iii) a rolling cadence to get better-sampled light curves at least in some time intervals. We also address the potential importance of Target of Opportunity (ToO) observations of blazar neutrino sources, and the advantages of a Minisurvey with a star trail cadence (see white paper by David Thomas et al.) for both the blazar science and the detection of possible very fast optical counterparts of FRBs.

The origin of the astrophysical high-energy neutrinos discovered by IceCube is currently a major mystery. The recent detection of IceCube-170922A, a $\sim$300 TeV neutrino potentially correlated with the flaring $\gamma$-ray source TXS 0506+056, directs attention toward BL Lac objects (BL Lacs), the subclass of blazars with weak emission lines. While high-energy neutrinos can be produced via photohadronic interactions between protons accelerated in their jets and ambient low-energy photons, the density of the latter in such objects had generally been thought to be too low for efficient neutrino emission. Here we consider the role of radiatively inefficient accretion flows (RIAFs), which can plausibly exist in the nuclei of BL Lacs, as the source of target photons for neutrino production. Based on simple model prescriptions for the spectra of RIAFs at different accretion rates, we find that they can be sufficienly intense to allow appreciable neutrino emission for the class of low-synchrotron-peak BL Lacs such as TXS 0506+056. In constrast, for high-synchrotron-peak BL Lacs including Mkn 421 and Mkn 501, the contribution of RIAFs is subdominant and their neutrino production efficiency can remain low, consistent with their non-detection by IceCube to date.

Aims: To investigate the extension of the very-high-energy spectral tail of the Crab pulsar at energies above 400 GeV. Methods: We analyzed $\sim$320 hours of good quality data of Crab with the MAGIC telescope, obtained from February 2007 until April 2014. Results: We report the most energetic pulsed emission ever detected from the Crab pulsar reaching up to 1.5 TeV. The pulse profile shows two narrow peaks synchronized with the ones measured in the GeV energy range. The spectra of the two peaks follow two different power-law functions from 70 GeV up to 1.5 TeV and connect smoothly with the spectra measured above 10 GeV by the Large Area Telescope (LAT) on board of the Fermi satellite. When making a joint fit of the LAT and MAGIC data, above 10 GeV, the photon indices of the spectra differ by 0.5$\pm$0.1. Conclusions: We measured with the MAGIC telescopes the most energetic pulsed photons from a pulsar to date. Such TeV pulsed photons require a parent population of electrons with a Lorentz factor of at least $5\times 10^6$. These results strongly suggest IC scattering off low energy photons as the emission mechanism and a gamma-ray production region in the vicinity of the light cylinder.

The discovery of high-energy astrophysical neutrinos by IceCube kicked off a new line of research to identify the electromagnetic counterparts producing these neutrinos. Among the extragalactic sources, active galactic nuclei (AGN), and in particular Blazars, are promising candidate neutrino emitters. Their structure, with a relativistic jet pointing to the Earth, offers a natural accelerator of particles and for this reason a perfect birthplace of high energy neutrinos. A good characterisation of the spectral energy distribution (SED) of these sources can improve the understanding of the physical composition of the source and the emission processes involved. Starting from our previous works in which we assumed a correlation between the $\gamma$-ray and the neutrino flux of the BL Lacs of the 2FHL catalogue (detected by Fermi above 50GeV), we select those BL Lac in spatial correlation with the IceCube events. We obtain a sample of 7 sources and we start an observational campaign to have a better characterisation of the synchrotron peak. During the analysis of the data a new source has been added because of its position inside the angular uncertainty of a muon track event detected by IceCube. This source, namely TXS0506+056, was in a high-state during the neutrino event and we will consider it as benchmark to check the proprieties of the other sources of the sample during the related neutrino detection. We obtain a better characterisation of the SED for the sources of our sample. A prospective extreme Blazar, a very peculiar low synchrotron peak (LSP) source with a large separation of the two peaks and a \textit{twin} of TXS0506+056 come up. We also provide the $\gamma$-ray light curve to check the trend of the sources around the neutrino detection but no clears patterns are in common among the sources.

Hadron beams are invoked to explain the peculiar properties of a subclass of BL Lac objects, the so-called extreme BL Lacs (EHBLs). This scenario predicts a quite distinctive feature for the high-energy gamma-ray spectrum of these sources, namely a hard energy tail extending up to $\sim100$ TeV. It has been proposed that the detection of this tail can offer an unambiguous way to distinguish between the hadron beam scenario and the standard one, which instead assumes gamma-ray emission from the jet strongly depleted at the highest energies ($E>30$ TeV) because of the interaction with the optical-IR cosmic radiation field. We present dedicated simulations of observations through the presently under construction Cherenkov Telescope Array (CTA) of the very-high energy spectrum of the prototypical EHBL 1ES 0229+200 assuming the two alternative models. We demonstrate that, considering 50 hours of observations from the southern site of CTA (the most sensitive at the highest energies), in the case of the hadron beam model it is possible to detect the source up to 100 TeV. This, together with the non detection of the source above 10 TeV in the standard case, ensures that CTA observations can be effectively used to unambiguously confirm or rule out the hadron beam scenario.

We report on observations of the pulsar / Be star binary system PSR J2032+4127 / MT91 213 in the energy range between 100 GeV and 20 TeV with the VERITAS and MAGIC imaging atmospheric 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 0.5 TeV. 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 (XRT) on board the Neil Gehrels \textit{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 we have observed, and we discuss the relationship between the binary system and the overlapping steady extended source, TeV J2032+4130.

Gamma-ray emitting narrow-line Seyfert 1 ($\gamma$-NLSy1) galaxies are thought to harbour relatively low-mass black holes (10$^6$-10$^8$ M$_{\odot}$) accreting close to the Eddington limit. They show characteristics similar to those of blazars, such as flux and spectral variability in the gamma-ray energy band and radio properties which point toward the presence of a relativistic jet. These characteristics make them an intriguing class of sources to be investigated with the Cherenkov Telescope Array (CTA), the next-generation ground-based gamma-ray observatory. We present our extensive set of simulations of all currently known $\gamma$-ray emitters identified as NLS1s (20 sources),investigating their detections and spectral properties, taking into account the effect of both the extra-galactic background light in the propagation of gamma-rays and intrinsic absorption components.We find that the prospects for observations of $\gamma$-NLSy1 with CTA are promising. In particular, the brightest sources of our sample, SBS 0846+513, PMN J0948+0022, and PKS 1502+036 can be detected during high/flaring states, the former two even in the case in which the emission occurs within the highly opaque central regions, which prevent $\gamma$ rays above few tens of GeV to escape. In this case the low-energy threshold of CTA will play a key role. If, on the other hand, high-energy emission occurs outside the broad line region, we can detect the sources up to several hundreds of GeV-depending on the intrinsic shape of the emitted spectrum. Therefore, CTA observations will provide valuable information on the physical conditions and emission properties of their jets.

PKS 1510-089 is a flat spectrum radio quasar strongly variable in the optical and GeV range. We search for low-state VHE gamma-ray emission from PKS 1510-089. We aim to characterize and model the source in a broad-band context, which would provide a baseline over which high states and flares could be better understood. 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 modelled within an External Compton scenario with a stationary emission region through which plasma and magnetic field are flowing. We perform also 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 hrs of data during times when the Fermi-LAT flux measured above 1 GeV was below 3x10-8cm-2s-1, which is the threshold adopted for the definition of a low gamma-ray activity state. The data show a strongly significant (9.5{\sigma}) VHE gamma-ray emission at the level of (4.27+-0.61stat)x10-12cm-2s-1 above 150GeV, a factor 80 smaller 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 broad-band emission is compatible with the EC 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% C.L. on the location of the emission region allows us to place it at the distance >74% of the outer radius of the BLR.

Fast radio bursts (FRBs) are bright flashes observed typically at GHz frequencies with millisecond duration, whose origin is likely extragalactic. Their nature remains mysterious, motivating searches for counterparts at other wavelengths. FRB 121102 is so far the only source known to repeatedly emit FRBs and is associated with a host galaxy at redshift $z \simeq 0.193$. We conducted simultaneous observations of FRB 121102 with the Arecibo and MAGIC telescopes during several epochs in 2016--2017. This allowed searches for millisecond-timescale burst emission in very-high-energy (VHE) gamma rays as well as the optical band. While a total of five FRBs were detected during these observations, no VHE emission was detected, neither of a persistent nature nor burst-like associated with the FRBs. The average integral flux upper limits above 100 GeV at 95% confidence level are $6.6 \times 10^{-12}~\mathrm{photons\ cm^{-2}\ s^{-1}}$ (corresponding to luminosity $L_{\rm VHE} \lesssim 10^{45}~\mathrm{erg\ s^{-1}}$) over the entire observation period, and $1.2 \times 10^{-7}~ \mathrm{photons\ cm^{-2}\ s^{-1}}$ ($L_{\rm VHE} \lesssim 10^{49}~\mathrm{erg\ s^{-1}}$) over the total duration of the five FRBs. We constrain the optical U-band flux to be below 8.6 mJy at 5-$\sigma$ level for 1-ms intervals around the FRB arrival times. A bright burst with U-band flux $29~\mathrm{mJy}$ and duration $\sim 12$ ms was detected 4.3 s before the arrival of one FRB. However, the probability of spuriously detecting such a signal within the sampled time space is 1.5% (2.2 $\sigma$, post-trial), i.e. consistent with the expected background. We discuss the implications of the obtained upper limits for constraining FRB models.

Clusters of galaxies are the largest known gravitationally bound structures in the Universe, with masses around $10^{15}\ M_\odot$, most of it in the form of dark matter. The ground-based Imaging Atmospheric Cherenkov Telescope MAGIC made a deep survey of the Perseus cluster of galaxies using almost 400 h of data recorded between 2009 and 2017. This is the deepest observational campaign so far on a cluster of galaxies in the very high energy range. We search for gamma-ray signals from dark matter particles in the mass range between 200 GeV and 200 TeV decaying into standard model pairs. We apply an analysis optimized for the spectral and morphological features expected from dark matter decays and find no evidence of decaying dark matter. From this, we conclude that dark matter particles have a decay lifetime longer than $\sim10^{26}$~s in all considered channels. Our results improve previous lower limits found by MAGIC and represent the strongest limits on decaying dark matter particles from ground-based gamma-ray instruments.

A multiwavelength campaign was organized 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. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0.2 TeV of $\sim$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 $\sim$3 CU, and the peak of the high-energy spectral component was found to be at $\sim$2 TeV. This study reports very hard X-ray spectra, and the hardest VHE spectra measured to date for Mrk 501. The fractional variability was found to increase with energy, with the highest variability occurring at VHE, and a significant correlation between the X-ray and VHE bands. 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 one-zone synchrotron self-Compton (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 behaviour 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.

e-ASTROGAM (enhanced ASTROGAM) is a breakthrough Observatory space mission, with a detector composed by a Silicon tracker, a calorimeter, and an anticoincidence system, dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV - the lower energy limit can be pushed to energies as low as 150 keV for the tracker, and to 30 keV for calorimetric detection. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LIGO-Virgo-GEO600-KAGRA, SKA, ALMA, E-ELT, TMT, LSST, JWST, Athena, CTA, IceCube, KM3NeT, and LISA.

A neutrino with energy of $\sim$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 $\sim 3\sigma$ level. We monitored the object in the very-high-energy (VHE) band with the MAGIC telescopes for $\sim$41 hours 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 $\sim 10^{45} - 4 \times 10^{46} {\rm erg \ s^{-1}}$. The steep spectrum observed by MAGIC is concordant with internal $\gamma\gamma$ absorption above a few tens of GeV entailed by photohadronic production of a $\sim$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 up-scattering 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 co-moving frame can be in the range $\sim 10^{14}$ to $10^{18}$ eV.

e-ASTROGAM is a concept for a breakthrough observatory space mission carrying a gamma-ray telescope dedicated to the study of the non-thermal Universe in the photon energy range from 0.15 MeV to 3 GeV. The lower energy limit can be pushed down to energies as low as 30 keV for gamma-ray burst detection with the calorimeter. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with remarkable polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous and current generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will be a major player of the multiwavelength, multimessenger time-domain astronomy of the 2030s, and provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LISA, LIGO, Virgo, KAGRA, the Einstein Telescope and the Cosmic Explorer, IceCube-Gen2 and KM3NeT, SKA, ALMA, JWST, E-ELT, LSST, Athena, and the Cherenkov Telescope Array.

The MAGIC telescopes observed S2 0109+22 in 2015 July during its flaring activity in high energy gamma rays observed by Fermi-LAT. We analyse the MAGIC data to characterise the very high energy (VHE) gamma-ray emission of S2 0109+22, which belongs to the subclass of intermediate synchrotron peak (ISP) BL Lac objects. We study the multi-frequency emission in order to investigate the source classification. Finally, we compare the source long-term behaviour to other VHE gamma-ray emitting (TeV) blazars. We performed a temporal and spectral analysis of the data centred around the MAGIC interval of observation (MJD 57225-57231). Long-term radio and optical data have also been investigated using the discrete correlation function. The redshift of the source is estimated through optical host-galaxy imaging and also using the amount of VHE gamma-ray absorption. The quasi-simultaneous multi-frequency spectral energy distribution (SED) is modelled with the conventional one-zone synchrotron self-Compton (SSC) model. MAGIC observations resulted in the detection of the source at a significance level of $5.3\,\sigma$. The VHE gamma-ray emission of S2 0109+22 is variable on a daily time scale. VHE gamma-ray luminosity of the source is lower than the average of TeV BL Lacs. The optical polarization, and long-term optical/radio behaviour of the source are different from the general population of TeV blazars. All these findings agree with the classification of the source as an ISP BL Lac object. We estimate the source redshift as $z = 0.36 \pm 0.07$. The SSC parameters describing the SED are rather typical for blazars.

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. The comprehensive dataset collected is investigated in order to shed light on the mechanism of the broadband emission. Multi-wavelength light curves have been studied together with the broadband Spectral Energy Distributions (SEDs). The data set collected spans from radio, optical photometry and polarimetry, X-ray, high-energy (HE, 0.1 GeV < E < 100 GeV) with \textit{Fermi}-LAT to the very-high-energy (VHE, E>100 GeV) with MAGIC. 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 VHE. In the constructed SED the \textit{Swift}-XRT+\textit{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 \textit{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 $\gamma$-ray data. Instead we use a two-zone model. The EVPA shows an unprecedented fast rotation. An estimation of the redshift of the source by combined HE and VHE data provides a value of $z = 0.31 \pm 0.02_{stats} \pm 0.05_{sys}$, confirming the literature value. The data show the VHE emission originating in the entrance and exit of a superluminal knot in and out 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 modelling.

The dissipation of energy flux in blazar jets plays a key role in the acceleration of relativistic particles. Two possibilities are commonly considered for the dissipation processes, magnetic reconnection -- possibly triggered by instabilities in magnetically-dominated jets -- , or shocks -- for weakly magnetized flows. We consider the polarimetric features expected for the two scenarios analyzing the results of state-of-the-art simulations. For the magnetic reconnection scenario we conclude, using results from global relativistic MHD simulations, that the emission likely occurs in turbulent regions with unstructured magnetic fields, although the simulations do not allow us to draw firm conclusions. On the other hand, with local particle-in-cell simulations we show that, for shocks with a magnetic field geometry suitable for particle acceleration, the self-generated magnetic field at the shock front is predominantly orthogonal to the shock normal and becomes quasi-parallel downstream. Based on this result we develop a simplified model to calculate the frequency-dependent degree of polarization, assuming that high-energy particles are injected at the shock and cool downstream. We apply our results to HBLs, blazars with the maximum of their synchrotron output at UV-soft X-ray energies. While in the optical band the predicted degree of polarization is low, in the X-ray emission it can ideally reach 50\%, especially during active/flaring states. The comparison between measurements in the optical and in the X-ray band made during active states (feasible with the planned {\it IXPE} satellite) are expected to provide valuable constraints on the dissipation and acceleration processes.

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 MAGIC telescopes. Observations were performed between 2016 September and 2017 February as part of a monitoring program. The brightest outburst with $\sim1.5$ times the Crab Nebula flux above 100 GeV (C.U.) was observed during the night between 2016 December 31 and 2017 January 1 (fifty times higher than the mean previously measured in two observational campaigns between 2009 and 2011). Significant variability of the day-by-day light curve was measured, the shortest flux-doubling time-scales was found to be of $(611\pm101)$ min. 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\pm35)$ GeV. Simultaneous optical flux density measurements in the R-band obtained with the KVA telescope are also presented and the correlation between the optical and gamma-ray emission is investigated. Due to possible internal pair-production, the fast flux variability constrains the Doppler factor to values which 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 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.

We study a sample of Narrow-line Seyfert 1 galaxies (NLS1) in their optical and radio features, in order to understand the differences between their radio silent, radio-loud and radio-quiet subclasses. We first show that the different redshift and mass distributions of radio-loud and -quiet NLS1s could be ascribed to observational biases. We then present a geometrical model according to which most of the different observational features of radio-loud and radio-quiet NLS1s are ascribed to the orientation of an intrinsically structured jet. We estimate the fraction of intrinsically jetted sources among NLS1s that justifies the observed radio-detected population. Noticeably, under the assumptions of the geometrical model, we derive a fraction of jetted sources significantly larger than in standard AGN.

A search for tau neutrino induced showers with the MAGIC telescopes is presented. The MAGIC telescopes located at an altitude of 2200 m a.s.l. in the Canary Island of La Palma, can point towards the horizon or a few degrees below across an azimuthal range of about 80 degrees. This provides a possibility to search for air showers induced by tau leptons arising from interactions of tau neutrinos in the Earth crust or the surrounding ocean. In this paper we show how such air showers can be discriminated from the background of very inclined hadronic showers by using Monte Carlo simulations. Taking into account the orography of the site, the point source acceptance and the event rates expected have been calculated for a sample of generic neutrino fluxes from photo-hadronic interactions in AGNs. The analysis of about 30 hours of data taken towards the sealeads to a 90\% C.L. point source limit for tau neutrinos in the energy range from $1.0 \times 10^{15}$ eV to $3.0 \times 10^{18}$ eV of about $E_{\nu_{\tau}}^{2}\times \phi (E_{\nu_{\tau}}) < 2.0 \times 10^{-4}$ GeV cm$^{-2}$ s$^{-1}$ for an assumed power-law neutrino spectrum with spectral index $\gamma$=-2. However, with 300 hours and in case of an optimistic neutrino flare model, limits of the level down to $E_{\nu_{\tau}}^{2}\times \phi (E_{\nu_{\tau}}) < 8.4 \times 10^{-6}$ GeV cm$^{-2}$ s$^{-1}$ can be expected.

Hard-TeV BL Lacs are a new type of blazars characterized by a hard intrinsic TeV spectrum, locating the peak of their gamma-ray emission in the spectral energy distribution (SED) above 2-10 TeV. Such high energies are problematic for the Compton emission, using a standard one-zone leptonic model. We study six examples of this new type of BL Lacs in the hard X-ray band with the NuSTAR satellite. Together with simultaneous observations with the SWIFT satellite, we fully constrain the peak of the synchrotron emission in their SED, and test the leptonic synchrotron self-Compton (SSC) model. We confirm the extreme nature of 5 objects also in the synchrotron emission. We do not find evidence of additional emission components in the hard X-ray band. We find that a one-zone SSC model can in principle reproduce the extreme properties of both peaks in the SED, from X-ray up to TeV energies, but at the cost of i) extreme electron energies with very low radiative efficiency, ii) conditions heavily out of equipartition (by 3 to 5 orders of magnitude), and iii) not accounting for the simultaneous UV data, which then should belong to a different emission component, possibly the same as the far-IR (WISE) data. We find evidence of this separation of the UV and X-ray emission in at least two objects. In any case, the TeV electrons must not "see" the UV or lower-energy photons, even if coming from different zones/populations, or the increased radiative cooling would steepen the VHE spectrum.