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.

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.

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.

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.

Among the four known transneptunian dwarf planets, Haumea is an exotic, very elongated, and fast rotating body. In contrast to the other dwarf planets, its size, shape, albedo, and density are not well constrained. Here we report results of a multi-chord stellar occultation, observed on 2017 January 21. Secondary events observed around the main body are consistent with the presence of a ring of opacity 0.5, width 70 km, and radius 2,287$_{-45}^{+75}$ km. The Centaur Chariklo was the first body other than a giant planet to show a ring system and the Centaur Chiron was later found to possess something similar to Chariklo's rings. Haumea is the first body outside the Centaur population with a ring. The ring is coplanar with both Haumea's equator and the orbit of its satellite Hi'iaka. Its radius places close to the 3:1 mean motion resonance with Haumea's spin period. The occultation by the main body provides an instantaneous elliptical limb with axes 1,704 $\pm$ 4 km x 1,138 $\pm$ 26 km. Combined with rotational light-curves, it constrains Haumea's 3D orientation and its triaxial shape, which is inconsistent with a homogeneous body in hydrostatic equilibrium. Haumea's largest axis is at least 2,322 $\pm$ 60 km, larger than thought before. This implies an upper limit of 1,885 $\pm$ 80 kg m$^{-3}$ for Haumea's density, smaller and less puzzling than previous estimations, and a geometric albedo of 0.51 $\pm$ 0.02, also smaller than previous estimations. No global N$_2$ or CH$_4$ atmosphere with pressures larger than 15 and 50 nbar (3-$\sigma$ limits), respectively, is detected.

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.

On 28th January 2018, the large Trans-Neptunian Object (TNO) 2002TC302 occulted a m$_v= $15.3 star with ID 130957813463146112 in the Gaia DR2 stellar catalog. 12 positive occultation chords were obtained from Italy, France, Slovenia and Switzerland. Also, 4 negative detections were obtained near the north and south limbs. This represents the best observed stellar occultation by a TNO other than Pluto, in terms of the number of chords published thus far. From the 12 chords, an accurate elliptical fit to the instantaneous projection of the body, compatible with the near misses, can be obtained. The resulting ellipse has major and minor axes of 543 $\pm$ 18 km and 460 $\pm$ 11 km, respectively, with a position angle of 3 $\pm$ 1 degrees for the minor axis. This information, combined with rotational light curves obtained with the 1.5m telescope at Sierra Nevada Observatory and the 1.23m telescope at Calar Alto observatory, allows us to derive possible 3D shapes and density estimations for the body, based on hydrostatic equilibrium assumptions. The effective area equivalent diameter is $\sim$ 84 km smaller than the radiometrically derived diameter using thermal data from Herschel and Spitzer Space Telescopes. This might indicate the existence of an unresolved satellite of up to $\sim$ 300 km in diameter, to account for all the thermal flux, although the occultation and thermal diameters are compatible within their error bars given the considerable uncertainty of the thermal results. The existence of a potential satellite also appears to be consistent with other ground-based data presented here. From the effective occultation diameter combined with H$_V$ measurements we derive a geometric albedo of 0.147 $\pm$ 0.005, which would be somewhat smaller if 2002TC302 has a satellite. The best occultation light curves do not show any signs of ring features or any signatures of a global atmosphere.

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 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.

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.

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.

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.

Spectropolarimetry enables us to measure the geometry and chemical structure of the ejecta in supernova explosions, which is fundamental for the understanding of their explosion mechanism(s) and progenitor systems. We collected archival data of 35 Type Ia Supernovae (SNe Ia), observed with FORS on the Very Large Telescope at 127 epochs in total. We examined the polarization of the Si II $\lambda$6355 $\AA$ line (p$_{\rm Si II}$) as a function of time which is seen to peak at a range of various polarization degrees and epochs relative to maximum brightness. We reproduced the $\Delta$m$_{15}$-p$_{\rm Si II}$ relationship identified in a previous study, and show that subluminous and transitional objects display polarization values below the $\Delta$m$_{15}$-p$_{\rm Si II}$ relationship for normal SNe Ia. We found a statistically significant linear relationship between the polarization of the Si II $\lambda$6355 $\AA$ line before maximum brightness and the Si II line velocity and suggest that this, along with the $\Delta$m$_{15}$-p$_{\rm Si II}$ relationship, may be explained in the context of a delayed-detonation model. In contrast, we compared our observations to numerical predictions in the $\Delta$m$_{15}$-v$_{\rm Si II}$ plane and found a dichotomy in the polarization properties between Chandrasekhar and sub-Chandrasekhar mass explosions, which supports the possibility of two distinct explosion mechanisms. A subsample of SNe display evolution of loops in the $q$-$u$ plane that suggests a more complex Si structure with depth. This insight, which could not be gleaned from total flux spectra, presents a new constraint on explosion models. Finally, we compared our statistical sample of the Si II polarization to quantitative predictions of the polarization levels for the double-detonation, delayed-detonation, and violent-merger models.

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}.

The pointing system of the prototype of the Large Size Telescope (LST-1) for the Cherenkov Telescope Array observatory, should ensure mapping of the gamma-ray image of a point-like source in the Cherenkov camera to the sky coordinates with a precision better than 14 arcseconds. Detailed studies of the telescope deformations are performed in order to disentangle different deformations and quantify their contributions to the miss-pointing, to learn how to correct for them, and finally how to design the system for offline and online pointing corrections. The LST-1 pointing precision system consist of several devices mounted at the center of the dish: Starguider Camera (SG), Camera Displacement Monitor (CDM), two inclinometers, four distance meters, and an Optical Axis Reference Laser (OARL), working together with the LEDs mounted in a circle around the Cherenkov camera. The online pointing corrections are based on a bending model as currently done by existing IACTs. The offline corrections will be performed combining measurements done by the SG and CDM cameras. SG will provide the position of the Cherenkov camera center with respect to the sky coordinates with a precision of 5 arcseconds, while CDM will provide the deviation of the telescope optical axis defined by the OARL spots with respect to the Cherenkov camera center with a precision better than 5 arcseconds. Laboratory measurements on dedicated test benches showed that the required pointing precision can be achieved for SG, CDM and inclinometer.

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 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.

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.

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.

Superluminous supernovae (SLSNe) are at least $\sim$5 times more luminous than common supernovae (SNe). Especially hydrogen-poor SLSN-I are difficult to explain with conventional powering mechanisms. One possible scenario that might explain such luminosities is that SLSNe-I are powered by an internal engine, such as a magnetar or an accreting black hole. Strong magnetic fields or collimated jets can circularly polarize light. In this work, we measured circular polarization of two SLSNe-I with the FOcal Reducer and low dispersion Spectrograph (FORS2) mounted at the ESO's Very Large Telescope (VLT). PS17bek, a fast evolving SLSN-I, was observed around peak, while OGLE16dmu, a slowly evolving SLSN-I, was observed 100 days after maximum. Neither SLSN shows evidence of circularly polarized light, however, these non-detections do not rule out the magnetar scenario as the powering engine for SLSNe-I. We calculate the strength of the magnetic field and the expected circular polarization as a function of distance from the magnetar, which decreases very fast. Additionally, we observed no significant linear polarization for PS17bek at four epochs, suggesting that the photosphere near peak is close to spherical symmetry.

Highly reddened type Ia Supernovae (SNe Ia) with low total-to-selective visual extinction ratio values, $R_V$, also show peculiar linear polarization wavelength dependencies with peak polarizations at short wavelengths ($\lambda_{max} \lesssim 0.4 \mu m$). It is not clear why sightlines to SNe Ia display such different continuum polarization profiles from interstellar sightlines in the Milky Way with similar $R_V$ values. We investigate polarization profiles of a sample of Galactic stars with low $R_V$ values, along anomalous extinction sightlines, with the aim to find similarities to the polarization profiles that we observe in SN Ia sightlines. We undertook spectropolarimetry of 14 stars, and used archival data for three additional stars, and run dust extinction and polarization simulations to infer a simple dust model that can reproduce the observed extinction and polarization curves. Our sample of Galactic stars with low $R_V$ values and anomalous extinction sightlines displays normal polarization profiles with an average $\lambda_{max} \sim 0.53 {\mu m}$, and is consistent within 3$\sigma$ to a larger coherent sample of Galactic stars from literature. Despite the low $R_V$ values of dust towards the stars in our sample, the polarization curves do not show any similarity to the continuum polarization curves observed towards SNe Ia with low $R_V$ values. There is a correlation between the best-fit Serkowski parameters $K$ and $\lambda_{max}$, but we did not find any significant correlation between $R_V$ and $\lambda_{max}$. Our simulations show that the $K-\lambda_{max}$ relationship is an intrinsic property of polarization. Furthermore, we have shown that in order to reproduce polarization curves with normal $\lambda_{max}$ and low $R_V$ values, a population of large (a $\geq 0.1 \mu m$) interstellar silicate grains must be contained in the dust's composition.

The centaur 95P/(2060) Chiron is showing comet-like activity since its discovery, but the mass-loss mechanisms triggering its activity remained unexplained. Although the collision rates in the centaur region are expected to be very low, and impacts are thought not to be responsible for the mass-loss, since the recent indications that Chiron might possess a ring similar to Chariklo's, and assuming that there is debris orbiting around, the impact triggered mass-loss mechanism should not be excluded as a possible cause of its activity. From time series observations collected on Calar Alto Observatory in Spain between 2014 and 2016, we found that the photometric scatter in Chiron's data is larger than a control star's scatter, indicating a possible microactivity, possibly caused by debris falling back to Chiron's surface and lifting small clouds of material. We also present rotational light curves, and measurements of Chiron's absolute magnitudes, that are consistent with the models supporting the presumption that Chiron possesses rings. By co-adding the images acquired in 2015, we have detected a $ \sim $5 arcsec long tail, showing a surface brightness of 25.3 mag(V)/arcsec$^{2}$.

The lines-of-sight to highly reddened SNe Ia show peculiar continuum polarization curves, growing toward blue wavelengths and peaking at $\lambda_{max} \lesssim 0.4 \mu m$, like no other sight line to any normal Galactic star. We examined continuum polarization measurements of a sample of asymptotic giant branch (AGB) and post-AGB stars from the literature, finding that some PPNe have polarization curves similar to those observed along SNe Ia sight lines. Those polarization curves are produced by scattering on circumstellar dust. We discuss the similarity and the possibility that at least some SNe Ia might explode during the post-AGB phase of their binary companion. Furthermore, we speculate that the peculiar SNe Ia polarization curves might provide observational support to the core-degenerate progenitor model.

Superluminous supernovae (SLSNe) are the most luminous supernovae in the universe. They are found in extreme star-forming galaxies and are probably connected with the death of massive stars. One hallmark of very massive progenitors would be a tendency to explode in very dense, UV-bright, and blue regions. In this paper we investigate the resolved host galaxy properties of two nearby hydrogen-poor SLSNe, PTF~11hrq and PTF~12dam. For both galaxies \textit{Hubble Space Telescope} multi-filter images were obtained. Additionally, we performe integral field spectroscopy of the host galaxy of PTF~11hrq using the Very Large Telescope Multi Unit Spectroscopic Explorer (VLT/MUSE), and investigate the line strength, metallicity and kinematics. Neither PTF~11hrq nor PTF~12dam occurred in the bluest part of their host galaxies, although both galaxies have overall blue UV-to-optical colors. The MUSE data reveal a bright starbursting region in the host of PTF~11hrq, although far from the SN location. The SN exploded close to a region with disturbed kinematics, bluer color, stronger [OIII], and lower metallicity. The host galaxy is likely interacting with a companion. PTF~12dam occurred in one of the brightest pixels, in a starbursting galaxy with a complex morphology and a tidal tail, where interaction is also very likely. We speculate that SLSN explosions may originate from stars generated during star-formation episodes triggered by interaction. High resolution imaging and integral field spectroscopy are fundamental for a better understanding of SLSNe explosion sites and how star formation varies across their host galaxies.

We reduced ESO's archival linear spectropolarimetry data (4000-9000\AA) of 6 highly polarized and 8 unpolarized standard stars observed between 2010 and 2016, for a total of 70 epochs, with the FOcal Reducer and low dispersion Spectrograph (FORS2) mounted at the Very Large Telescope. We provide very accurate standard stars polarization measurements as a function of wavelength, and test the performance of the spectropolarimetric mode (PMOS) of FORS2. We used the unpolarized stars to test the time stability of the PMOS mode, and found a small ($\leq$0.1%), but statistically significant, on-axis instrumental polarization wavelength dependency, possibly caused by the tilted surfaces of the dispersive element. The polarization degree and angle are found to be stable at the level of $\leq$0.1% and $\leq$0.2 degrees, respectively. We derived the polarization wavelength dependence of the polarized standard stars and found that, in general, the results are consistent with those reported in the literature, e.g. Fossati et al. (2007) who performed a similar analysis using FORS1 data. The re-calibrated data provide a very accurate set of standards that can be very reliably used for technical and scientific purposes. The analysis of the Serkowski parameters revealed a systematic deviation from the width parameter $K$ reported by Whittet et al. (1992). This is most likely explained by incorrect effective wavelengths adopted in that study for the R and I bands.

It is well known that some Main Belt asteroids show comet-like features. A representative example is the first known Main Belt comet 133P/(7968) Elst-Pizarro. If the mechanisms causing this activity are too weak to develop visually evident comae or tails, the objects stay unnoticed. We are presenting a novel way to search for active asteroids, based on looking for objects with deviations from their expected brightnesses in a database. Just by using the MPCAT-OBS Observation Archive we have found five new candidate objects that possibly show a type of comet-like activity, and the already known Main Belt comet 133P/(7968) Elst-Pizarro. Four of the new candidates, (315) Constantia, (1026) Ingrid, (3646) Aduatiques, and (24684) 1990 EU4, show brightness deviations independent of the object's heliocentric distance, while (35101) 1991 PL16 shows deviations dependent on its heliocentric distance, which could be an indication of a thermal triggered mechanism. The method could be implemented in future sky survey programmes to detect outbursts on Main Belt objects almost simultaneously with their occurrence.

Near-Earth asteroid 2012 DA14 made its closest approach on February 15, 2013, when it passed at a distance of 27,700 km from the Earth's surface. It was the first time an asteroid of moderate size was predicted to approach that close to the Earth, becoming bright enough to permit a detailed study from ground-based telescopes. Asteroid 2012 DA14 was poorly characterized before its closest approach. We acquired data using several telescopes on four Spanish observatories: the 10.4m Gran Telescopio Canarias (GTC) and the 3.6m Telescopio Nazionale Galileo (TNG), both in the El Roque de los Muchachos Observatory (ORM, La Palma); the 2.2m CAHA telescope, in the Calar Alto Observatory (Almeria); the f/3 0.77m telescope in the La Hita Observatory (Toledo); and the f/8 1.5m telescope in the Sierra Nevada Observatory (OSN, Granada). We obtained visible and near-infrared color photometry, visible spectra and time-series photometry. Visible spectra together with color photometry of 2012 DA14 show that it can be classified as an L-type asteroid, a rare spectral type with a composition similar to that of carbonaceous chondrites. The time-series photometry provides a rotational period of 8.95 +- 0.08 hours after the closest approach, and there are indications that the object suffered a spin-up during this event. The large amplitude of the light curve suggests that the object is very elongated and irregular, with an equivalent diameter of around 18m. We obtain an absolute magnitude of H_R = 24.5 +- 0.2, corresponding to H_V = 25.0 +- 0.2. The GTC photometry also gives H_V = 25.29 +- 0.14. Both values agree with the value listed at the Minor Planet Center shortly after discovery. From the absolute photometry, together with some constraints on size and shape, we compute a geometric albedo of p_V = 0.44 +- 0.20, which is slightly above the range of albedos known for L-type asteroids (0.082 - 0.405).

From CCD observations of a fixed and large star field that contained the binary TNO Orcus, we have been able to derive high-precision relative astrometry and photometry of the Orcus system with respect to background stars. The RA residuals of an orbital fit to the astrometric data revealed a periodicity of 9.7+-0.3 days, which is what one would expect to be induced by the known Orcus companion. The residuals are also correlated with the theoretical positions of the satellite with regard to the primary. We therefore have revealed the presence of Orcus' satellite in our astrometric measurements. The photocenter motion is much larger than the motion of Orcus around the barycenter, and we show here that detecting some binaries through a carefully devised astrometric technique might be feasible with telescopes of moderate size. We also analyzed the system's mid-term photometry to determine whether the rotation could be tidally locked to the satellite's orbital period. We found that a photometric variability of 9.7+-0.3 days is clear in our data, and is nearly coincident with the orbital period of the satellite. We believe this variability might be induced by the satellite's rotation. There is also a slight hint for an additional small variability in the 10 hr range that was already reported in the literature. This short-term variability would indicate that the primary is not tidally locked and therefore the system would not have reached a double synchronous state. Implications for the basic physical properties of the primary and its satellite are discussed. From angular momentum considerations we suspect that the Orcus satellite might have formed from a rotational fission. This requires that the mass of the satellite would be around 0.09 times that of the primary, close to the value that one derives by using an albedo of 0.12 for the satellite and assuming equal densities for both objects.