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 flat-spectrum radio quasar PKS 1441+25 at a redshift of z = 0.940 is detected between 40 and 250 GeV with a significance of 25.5 {\sigma} using the MAGIC telescopes. Together with the gravitationally lensed blazar QSO B0218+357 (z = 0.944), PKS 1441+25 is the most distant very high energy (VHE) blazar detected to date. The observations were triggered by an outburst in 2015 April seen at GeV energies with the Large Area Telescope on board Fermi. Multi-wavelength observations suggest a subdivision of the high state into two distinct flux states. In the band covered by MAGIC, the variability time scale is estimated to be 6.4 +/- 1.9 days. Modeling the broadband spectral energy distribution with an external Compton model, the location of the emitting region is understood as originating in the jet outside the broad line region (BLR) during the period of high activity, while being partially within the BLR during the period of low (typical) activity. The observed VHE spectrum during the highest activity is used to probe the extragalactic background light at an unprecedented distance scale for ground-based gamma-ray astronomy.

In this work we present data from observations with the MAGIC telescopes of SN 2014J detected in January 21 2014, the closest Type Ia supernova since Imaging Air Cherenkov Telescopes started to operate. We probe the possibility of very-high-energy (VHE; $E\geq100$ GeV) gamma rays produced in the early stages of Type Ia supernova explosions. We performed follow-up observations after this supernova explosion for 5 days, between January 27 and February 2 in 2014. We search for gamma-ray signal in the energy range between 100 GeV and several TeV from the location of SN 2014J using data from a total of $\sim5.5$ hours of observations. Prospects for observing gamma-rays of hadronic origin from SN 2014J in the near future are also being addressed. No significant excess was detected from the direction of SN 2014J. Upper limits at 95$\%$ confidence level on the integral flux, assuming a power-law spectrum, d$F/$d$E\propto E^{-\Gamma}$, with a spectral index of $\Gamma=2.6$, for energies higher than 300 GeV and 700 GeV, are established at $1.3\times10^{-12}$ and $4.1\times10^{-13}$ photons~cm$^{-2}$s$^{-1}$, respectively. For the first time, upper limits on the VHE emission of a Type Ia supernova are established. The energy fraction isotropically emitted into TeV gamma rays during the first $\sim10$ days after the supernova explosion for energies greater than 300 GeV is limited to $10^{-6}$ of the total available energy budget ($\sim 10^{51}$ erg). Within the assumed theoretical scenario, the MAGIC upper limits on the VHE emission suggest that SN 2014J will not be detectable in the future by any current or planned generation of Imaging Atmospheric Cherenkov Telescopes.

We present the first detection of the nearby (z=0.084) low-luminosity BL Lac object 1ES 1741+196 in the very high energy (VHE: E$>$100 GeV) band. This object lies in a triplet of interacting galaxies. Early predictions had suggested 1ES 1741+196 to be, along with several other high-frequency BL Lac sources, within the reach of MAGIC detectability. Its detection by MAGIC, later confirmed by VERITAS, helps to expand the small population of known TeV BL Lacs. The source was observed with the MAGIC telescopes between 2010 April and 2011 May, collecting 46 h of good quality data. These observations led to the detection of the source at 6.0 $\sigma$ confidence level, with a steady flux $\mathrm{F}(> 100 {\rm GeV}) = (6.4 \pm 1.7_{\mathrm{stat}}\pm 2.6_{\mathrm{syst}}) \cdot 10^{-12}$ ph cm$^{-2}$ s$^{-1}$ and a differential spectral photon index $\Gamma = 2.4 \pm 0.2_{\mathrm{stat}} \pm 0.2_{\mathrm{syst}}$ in the range of $\sim$80 GeV - 3 TeV. To study the broad-band spectral energy distribution (SED) simultaneous with MAGIC observations, we use KVA, Swift/UVOT and XRT, and Fermi/LAT data. One-zone synchrotron-self-Compton (SSC) modeling of the SED of 1ES 1741+196 suggests values for the SSC parameters that are quite common among known TeV BL Lacs except for a relatively low Doppler factor and slope of electron energy distribution. A thermal feature seen in the SED is well matched by a giant elliptical's template. This appears to be the signature of thermal emission from the host galaxy, which is clearly resolved in optical observations.

We present an extensive study of the BL Lac object Mrk 501 based on a data set collected during the multi-instrument campaign spanning from 2009 March 15 to 2009 August 1 which includes, among other instruments, MAGIC, VERITAS, Whipple 10-m, Fermi-LAT, RXTE, Swift, GASP-WEBT and VLBA. We find an increase in the fractional variability with energy, while no significant interband correlations of flux changes are found in the acquired data set. The higher variability in the very high energy (>100 GeV, VHE) gamma-ray emission and the lack of correlation with the X-ray emission indicate that the highest-energy electrons that are responsible for the VHE gamma-rays do not make a dominant contribution to the ~1 keV emission. Alternatively, there could be a very variable component contributing to the VHE gamma-ray emission in addition to that coming from the synchrotron self-Compton (SSC) scenarios. The space of SSC model parameters is probed following a dedicated grid-scan strategy, allowing for a wide range of models to be tested and offering a study of the degeneracy of model-to-data agreement in the individual model parameters. We find that there is some degeneracy in both the one-zone and the two-zone SSC scenarios that were probed, with several combinations of model parameters yielding a similar model-to-data agreement, and some parameters better constrained than others. The SSC model grid-scan shows that the flaring activity around 2009 May 22 cannot be modeled adequately with a one-zone SSC scenario, while it can be suitably described within a two-independent-zone SSC scenario. The observation of an electric vector polarization angle rotation coincident with the gamma-ray flare from 2009 May 1 resembles those reported previously for low frequency peaked blazars, hence suggesting that there are many similarities in the flaring mechanisms of blazars with different jet properties.

Context. We present the results of a multi-year monitoring campaign of the Galactic Center (GC) with the MAGIC telescopes. These observations were primarily motivated by reports that a putative gas cloud (G2) would be passing in close proximity to the super-massive black hole (SMBH), associated with Sagittarius A*, located at the center of our galaxy. This event was expected to give astronomers a unique chance to study the effect of in-falling matter on the broad-band emission of a SMBH. Aims. We search for potential flaring emission of very-high-energy (VHE; $\geq$100 GeV) gamma rays from the direction of the SMBH at the GC due to the passage of the G2 object. Using these data we also study the morphology of this complex region. Methods. We observed the GC region with the MAGIC Imaging Atmospheric Cherenkov Telescopes during the period 2012-2015, collecting 67 hours of good-quality data. In addition to a search for variability in the flux and spectral shape of the GC gamma-ray source, we use a point-source subtraction technique to remove the known gamma-ray emitters located around the GC in order to reveal the TeV morphology of the extended emission inside that region. Results. No effect of the G2 object on the VHE gamma-ray emission from the GC was detected during the 4 year observation campaign. We confirm previous measurements of the VHE spectrum of Sagittarius A*, and do not detect any significant variability of the emission from the source. Furthermore, the known VHE gamma-ray emitter at the location of the supernova remnant G0.9+0.1 was detected, as well as the recently discovered VHE source close to the GG radio Arc.

We describe and report the status of a neutrino-triggered program in IceCube that generates real-time alerts for gamma-ray follow-up observations by atmospheric-Cherenkov telescopes (MAGIC and VERITAS). While IceCube is capable of monitoring the whole sky continuously, high-energy gamma-ray telescopes have restricted fields of view and in general are unlikely to be observing a potential neutrino-flaring source at the time such neutrinos are recorded. The use of neutrino-triggered alerts thus aims at increasing the availability of simultaneous multi-messenger data during potential neutrino flaring activity, which can increase the discovery potential and constrain the phenomenological interpretation of the high-energy emission of selected source classes (e.g. blazars). The requirements of a fast and stable online analysis of potential neutrino signals and its operation are presented, along with first results of the program operating between 14 March 2012 and 31 December 2015.

The MAGIC (Major Atmospheric Gamma-ray Imaging Cherenkov) telescopes observed the BL Lac object H1722+119 (redshift unknown) for six consecutive nights between 2013 May 17 and 22, for a total of 12.5 h. The observations were triggered by high activity in the optical band measured by the KVA (Kungliga Vetenskapsakademien) telescope. The source was for the first time detected in the very high energy (VHE, $E > 100$ GeV) $\gamma$-ray band with a statistical significance of 5.9 $\sigma$. The integral flux above 150 GeV is estimated to be $(2.0\pm 0.5)$ per cent of the Crab Nebula flux. We used contemporaneous high energy (HE, 100 MeV $ < E < 100$ GeV) $\gamma$-ray observations from Fermi-LAT (Large Area Telescope) to estimate the redshift of the source. Within the framework of the current extragalactic background light models, we estimate the redshift to be $z = 0.34 \pm 0.15$. Additionally, we used contemporaneous X-ray to radio data collected by the instruments on board the Swift satellite, the KVA, and the OVRO (Owens Valley Radio Observatory) telescope to study multifrequency characteristics of the source. We found no significant temporal variability of the flux in the HE and VHE bands. The flux in the optical and radio wavebands, on the other hand, did vary with different patterns. The spectral energy distribution (SED) of H1722+119 shows surprising behaviour in the $\sim 3\times10^{14} - 10^{18}$ Hz frequency range. It can be modelled using an inhomogeneous helical jet synchrotron self-Compton model.

Context. QSO B0218+357 is a gravitationally lensed blazar located at a redshift of 0.944. The gravitational lensing splits the emitted radiation into two components, spatially indistinguishable by gamma-ray instruments, but separated by a 10-12 day delay. In July 2014, QSO B0218+357 experienced a violent flare observed by the Fermi-LAT and followed by the MAGIC telescopes. Aims. The spectral energy distribution of QSO B0218+357 can give information on the energetics of z ~ 1 very high energy gamma- ray sources. Moreover the gamma-ray emission can also be used as a probe of the extragalactic background light at z ~ 1. Methods. MAGIC performed observations of QSO B0218+357 during the expected arrival time of the delayed component of the emission. The MAGIC and Fermi-LAT observations were accompanied by quasi-simultaneous optical data from the KVA telescope and X-ray observations by Swift-XRT. We construct a multiwavelength spectral energy distribution of QSO B0218+357 and use it to model the source. The GeV and sub-TeV data, obtained by Fermi-LAT and MAGIC, are used to set constraints on the extragalactic background light. Results. Very high energy gamma-ray emission was detected from the direction of QSO B0218+357 by the MAGIC telescopes during the expected time of arrival of the trailing component of the flare, making it the farthest very high energy gamma-ray sources detected to date. The observed emission spans the energy range from 65 to 175 GeV. The combined MAGIC and Fermi-LAT spectral energy distribution of QSO B0218+357 is consistent with current extragalactic background light models. The broad band emission can be modeled in the framework of a two zone external Compton scenario, where the GeV emission comes from an emission region in the jet, located outside the broad line region.

1ES 1011+496 $(z=0.212)$ was discovered in very high energy (VHE, E >100 GeV) $\gamma$-rays with MAGIC in 2007. The absence of simultaneous data at lower energies led to a rather incomplete characterization of the broadband spectral energy distribution (SED). We study the source properties and the emission mechanisms, probing whether a simple one-zone synchrotron-self-Compton (SSC) scenario is able to explain the observed broadband spectrum. We analyzed VHE to radio data from 2011 and 2012 collected by MAGIC, $Fermi$-LAT, $Swift$, KVA, OVRO, and Mets\"ahovi in addition to optical polarimetry data and radio maps from the Liverpool Telescope and MOJAVE. The VHE spectrum was fit with a simple power law with a photon index of $3.69\pm0.22$ and a flux above 150 GeV of $(1.46\pm0.16)\times10^{-11}$ ph cm$^{-2}$ s$^{-1}$. 1ES 1011+496 was found to be in a generally quiescent state at all observed wavelengths, showing only moderate variability from radio to X-rays. A low degree of polarization of less than 10% was measured in optical, while some bright features polarized up to 60% were observed in the radio jet. A similar trend in the rotation of the electric vector position angle was found in optical and radio. The radio maps indicated a superluminal motion of $1.8\pm0.4\,c$, which is the highest speed statistically significantly measured so far in a high-frequency-peaked BL Lac. For the first time, the high-energy bump in the broadband SED of 1ES 1011+496 could be fully characterized from 0.1 GeV to 1 TeV which permitted a more reliable interpretation within the one-zone SSC scenario. The polarimetry data suggest that at least part of the optical emission has its origin in some of the bright radio features, while the low polarization in optical might be due to the contribution of parts of the radio jet with different orientations of the magnetic field to the optical emission.

We study the multi-band variability and correlations of the TeV blazar Mrk 421 on year time scales, which can bring additional insight on the processes responsible for its broadband emission. We observed Mrk 421 in the very high energy (VHE) gamma-ray range with the Cherenkov telescope MAGIC-I from March 2007 to June 2009 for a total of 96 hours of effective time after quality cuts. The VHE flux variability is quantified with several methods, including the Bayesian Block algorithm, which is applied to data from Cherenkov telescopes for the first time. The 2.3 year long MAGIC light curve is complemented with data from the Swift/BAT and RXTE/ASM satellites and the KVA, GASP-WEBT, OVRO, and Mets\"ahovi telescopes from February 2007 to July 2009, allowing for an excellent characterisation of the multi-band variability and correlations over year time scales. Mrk 421 was found in different gamma-ray emission states during the 2.3 year long observation period. Flares and different levels of variability in the gamma-ray light curve could be identified with the Bayesian Block algorithm. The same behaviour of a quiet and active emission was found in the X-ray light curves measured by Swift/BAT and the RXTE/ASM, with a direct correlation in time. The behaviour of the optical light curve of GASP-WEBT and the radio light curves by OVRO and Mets\"ahovi are different as they show no coincident features with the higher energetic light curves and a less variable emission. The fractional variability is overall increasing with energy. The comparable variability in the X-ray and VHE bands and their direct correlation during both high- and low-activity periods spanning many months show that the electron populations radiating the X-ray and gamma-ray photons are either the same, as expected in the Synchrotron-Self-Compton mechanism, or at least strongly correlated, as expected in electromagnetic cascades.

The gamma-ray binary LS I +61$^{\circ}$303 is a well established source from centimeter radio up to very high energy (VHE; E$>$100 GeV). Its broadband emission shows a periodicity of $\sim$26.5 days, coincident with the orbital period. A longer (super-orbital) period of 1667 $\pm$ 8 days was discovered in radio and confirmed in optical and high energy (HE; E>100 MeV) gamma-ray observations. We present a four-year campaign performed by MAGIC together with archival data concentrating on a search for a long timescale signature in the VHE emission. We focus on the search for super-orbital modulation of the VHE peak and on the search for correlations between TeV emission and optical determination of the extension of the circumstellar disk. A four-year campaign has been carried out by MAGIC. The source was observed during the orbital phases when the periodic VHE outbursts have occurred ($\phi$=0.55-0.75). Additionally, we included archival MAGIC observations and data published by the VERITAS collaboration in these studies. For the correlation studies, LS I +61$^{\circ}$303 has also been observed during the orbital phases where sporadic VHE emission had been detected in the past ($\phi$=0.75-1.0). These MAGIC observations were simultaneous with optical spectroscopy from the LIVERPOOL telescope. The TeV flux of the periodical outburst in orbital phases $\phi$=0.5--0.75 was found to show yearly variability consistent with the $\sim$4.5 years long-term modulation found in the radio band. This modulation of the TeV flux can be well described by a sine function with the best fit period of $1610\pm 58$ days. The complete dataset span two super-orbital periods. There is no evidence for a correlation between the TeV emission and the mass-loss rate of the Be star but this may be affected by the strong, short timescale (as short as intra-day) variation displayed by the H$\alpha$ fluxes.

The BL Lac object 1ES 1011+496 was discovered at Very High Energy gamma-rays by MAGIC in spring 2007. Before that the source was little studied in different wavelengths. Therefore a multi-wavelength (MWL) campaign was organized in spring 2008. Along MAGIC, the MWL campaign included the Metsahovi radio observatory, Bell and KVA optical telescopes and the Swift and AGILE satellites. MAGIC observations span from March to May, 2008 for a total of 27.9 hours, of which 19.4 hours remained after quality cuts. The light curve showed no significant variability. The differential VHE spectrum could be described with a power-law function. Both results were similar to those obtained during the discovery. Swift XRT observations revealed an X-ray flare, characterized by a harder when brighter trend, as is typical for high synchrotron peak BL Lac objects (HBL). Strong optical variability was found during the campaign, but no conclusion on the connection between the optical and VHE gamma-ray bands could be drawn. The contemporaneous SED shows a synchrotron dominated source, unlike concluded in previous work based on nonsimultaneous data, and is well described by a standard one zone synchrotron self Compton model. We also performed a study on the source classification. While the optical and X-ray data taken during our campaign show typical characteristics of an HBL, we suggest, based on archival data, that 1ES 1011+496 is actually a borderline case between intermediate and high synchrotron peak frequency BL Lac objects.

The Geminga pulsar, one of the brighest gamma-ray sources, is a promising candidate for emission of very-high-energy (VHE > 100 GeV) pulsed gamma rays. Also, detection of a large nebula have been claimed by water Cherenkov instruments. We performed deep observations of Geminga with the MAGIC telescopes, yielding 63 hours of good-quality data, and searched for emission from the pulsar and pulsar wind nebula. We did not find any significant detection, and derived 95% confidence level upper limits. The resulting upper limits of 5.3 x 10^{-13} TeV cm^{-2} s^{-1} for the Geminga pulsar and 3.5 x 10^{-12} TeV cm^{-2} s^{-1} for the surrounding nebula at 50 GeV are the most constraining ones obtained so far at VHE. To complement the VHE observations, we also analyzed 5 years of Fermi-LAT data from Geminga, finding that the sub-exponential cut-off is preferred over the exponential cut-off that has been typically used in the literature. We also find that, above 10 GeV, the gamma-ray spectra from Geminga can be described with a power law with index softer than 5. The extrapolation of the power-law Fermi-LAT pulsed spectra to VHE goes well below the MAGIC upper limits, indicating that the detection of pulsed emission from Geminga with the current generation of Cherenkov telescopes is very difficult.

In February-March 2014, the MAGIC telescopes observed the high-frequency peaked BL Lac 1ES 1011+496 (z=0.212) in flaring state at very-high energy (VHE, E>100GeV). The flux reached a level more than 10 times higher than any previously recorded flaring state of the source. We present the description of the characteristics of the flare presenting the light curve and the spectral parameters of the night-wise spectra and the average spectrum of the whole period. From these data we aim at detecting the imprint of the Extragalactic Background Light (EBL) in the VHE spectrum of the source, in order to constrain its intensity in the optical band. For this we implement the method developed by the H.E.S.S. collaboration in which the intrinsic energy spectrum of the source is modeled with a simple function, and the EBL-induced optical depth is calculated using a template EBL model. The likelihood of the observed spectrum is then maximized, including a normalization factor for the EBL opacity among the free parameters. From the data collected differential energy spectra was produced for all nights of the observed period. Evaluating the changes in the fit parameters we conclude that the spectral shape for most of the nights were compatible, regardless of the flux level, which enabled us to produce an average spectrum from which the EBL imprint could be constrained. The likelihood ratio test shows that the model with an EBL density 1.07(-0.20,+0.24)_{stat+sys}, relative to the one in the tested EBL template (Dominguez et al.2011), is preferred at the 4.6 sigma level to the no-EBL hypothesis, with the assumption that the intrinsic source spectrum can be modeled as a log-parabola. This would translate into a constraint of the EBL density in the wavelength range [0.24 um,4.25 um], with a peak value at 1.4 um of F=12.27_{-2.29}^{+2.75} nW m^{-2} sr^{-1}, including systematics.

We present the first joint analysis of gamma-ray data from the MAGIC Cherenkov telescopes and the Fermi Large Area Telescope (LAT) to search for gamma-ray signals from dark matter annihilation in dwarf satellite galaxies. We combine 158 hours of Segue 1 observations with MAGIC with 6-year observations of 15 dwarf satellite galaxies by the Fermi-LAT. We obtain limits on the annihilation cross-section for dark matter particle masses between 10 GeV and 100 TeV - the widest mass range ever explored by a single gamma-ray analysis. These limits improve on previously published Fermi-LAT and MAGIC results by up to a factor of two at certain masses. Our new inclusive analysis approach is completely generic and can be used to perform a global, sensitivity-optimized dark matter search by combining data from present and future gamma-ray and neutrino detectors.

Clusters of galaxies are expected to be reservoirs of cosmic rays (CRs) that should produce diffuse gamma-ray emission due to their hadronic interactions with the intra-cluster medium. The nearby Perseus cool-core cluster, identified as the most promising target to search for such an emission, has been observed with the MAGIC telescopes at very-high energies (VHE, E>100 GeV) for a total of 253 hr from 2009 to 2014. The active nuclei of NGC 1275, the central dominant galaxy of the cluster, and IC 310, lying at about 0.6$^\circ$ from the centre, have been detected as point-like VHE gamma-ray emitters during the first phase of this campaign. We report an updated measurement of the NGC 1275 spectrum, which is well described by a power law with a photon index of $3.6\pm0.2_{stat}\pm0.2_{syst}$ between 90 GeV and 1.2 TeV. We do not detect any diffuse gamma-ray emission from the cluster and set stringent constraints on its CR population. In order to bracket the uncertainties over the CR spatial and spectral distributions, we adopt different spatial templates and power-law spectral indexes $\alpha$. For $\alpha=2.2$, the CR-to-thermal pressure within the cluster virial radius is constrained to be below 1-2%, except if CRs can propagate out of the cluster core, generating a flatter radial distribution and releasing the CR-to-thermal pressure constraint to <20%. Assuming that the observed radio mini-halo of Perseus is generated by secondary electrons from CR hadronic interactions, we can derive lower limits on the central magnetic field, $B_0$, that depend on the CR distribution. For $\alpha=2.2$, $B_0\gtrsim5-8 \mu$G, which is below the 25 $\mu$G inferred from Faraday rotation measurements, whereas, for $\alpha\lesssim2.1$, the hadronic interpretation of the diffuse radio emission is in contrast with our gamma-ray flux upper limits independently of the magnetic field strength.

We present coordinated multiwavelength observations of the bright, nearby BL Lac object Mrk 421 taken in 2013 January-March, involving GASP-WEBT, Swift, NuSTAR, Fermi-LAT, MAGIC, VERITAS, and other collaborations and instruments, providing data from radio to very-high-energy (VHE) gamma-ray bands. NuSTAR yielded previously unattainable sensitivity in the 3-79 keV range, revealing that the spectrum softens when the source is dimmer until the X-ray spectral shape saturates into a steep power law with a photon index of approximately 3, with no evidence for an exponential cutoff or additional hard components up to about 80 keV. For the first time, we observed both the synchrotron and the inverse-Compton peaks of the spectral energy distribution (SED) simultaneously shifted to frequencies below the typical quiescent state by an order of magnitude. The fractional variability as a function of photon energy shows a double-bump structure which relates to the two bumps of the broadband SED. In each bump, the variability increases with energy which, in the framework of the synchrotron self-Compton model, implies that the electrons with higher energies are more variable. The measured multi-band variability, the significant X-ray-to-VHE correlation down to some of the lowest fluxes ever observed in both bands, the lack of correlation between optical/UV and X-ray flux, the low degree of polarization and its significant (random) variations, the short estimated electron cooling time, and the significantly longer variability timescale observed in the NuSTAR light curves point toward in-situ electron acceleration, and suggest that there are multiple compact regions contributing to the broadband emission of Mrk 421 during low-activity states.

We report on simultaneous broadband observations of the TeV-emitting blazar Markarian 501 between 1 April and 10 August 2013, including the first detailed characterization of the synchrotron peak with Swift and NuSTAR. During the campaign, the nearby BL Lac object was observed in both a quiescent and an elevated state. The broadband campaign includes observations with NuSTAR, MAGIC, VERITAS, the Fermi Large Area Telescope (LAT), Swift X-ray Telescope and UV Optical Telescope, various ground-based optical instruments, including the GASP-WEBT program, as well as radio observations by OVRO, Mets\"ahovi and the F-Gamma consortium. Some of the MAGIC observations were affected by a sand layer from the Saharan desert, and had to be corrected using event-by-event corrections derived with a LIDAR (LIght Detection And Ranging) facility. This is the first time that LIDAR information is used to produce a physics result with Cherenkov Telescope data taken during adverse atmospheric conditions, and hence sets a precedent for the current and future ground-based gamma-ray instruments. The NuSTAR instrument provides unprecedented sensitivity in hard X-rays, showing the source to display a spectral energy distribution between 3 and 79 keV consistent with a log-parabolic spectrum and hard X-ray variability on hour timescales. None (of the four extended NuSTAR observations) shows evidence of the onset of inverse-Compton emission at hard X-ray energies. We apply a single-zone equilibrium synchrotron self-Compton model to five simultaneous broadband spectral energy distributions. We find that the synchrotron self-Compton model can reproduce the observed broadband states through a decrease in the magnetic field strength coinciding with an increase in the luminosity and hardness of the relativistic leptons responsible for the high-energy emission.

We have performed a search for muon neutrinos from dark matter annihilation in the center of the Sun with the 79-string configuration of the IceCube neutrino telescope. For the first time, the DeepCore sub-array is included in the analysis, lowering the energy threshold and extending the search to the austral summer. The 317 days of data collected between June 2010 and May 2011 are consistent with the expected background from atmospheric muons and neutrinos. Upper limits are therefore set on the dark matter annihilation rate, with conversions to limits on spin-dependent and spin-independent WIMP-proton cross-sections for WIMP masses in the range 20 - 5000 GeV. These are the most stringent spin-dependent WIMP-proton cross-sections limits to date above 35 GeV.

Context. In the last five years the Fermi Large Area Telescope (LAT) instrument detected GeV {\gamma}-ray emission from five novae. The GeV emission can be interpreted in terms of an inverse Compton process of electrons accelerated in a shock. In this case it is expected that protons in the same conditions can be accelerated to much higher energies. Consequently they may produce a second component in the {\gamma}-ray spectrum at TeV energies. Aims. We aim to explore the very-high-energy domain to search for {\gamma}-ray emission above 50 GeV and to shed light on the acceleration process of leptons and hadrons in nova explosions. Methods. We have performed observations with the MAGIC telescopes of the classical nova V339 Del shortly after the 2013 outburst, triggered by optical and subsequent GeV {\gamma}-ray detec- tions. We also briefly report on VHE observations of the symbiotic nova YY Her and the dwarf nova ASASSN-13ax. We complement the TeV MAGIC observations with the analysis of con- temporaneous Fermi-LAT data of the sources. The TeV and GeV observations are compared in order to evaluate the acceleration parameters for leptons and hadrons. Results. No significant TeV emission was found from the studied sources. We computed upper limits on the spectrum and night-by-night flux. The combined GeV and TeV observations of V339 Del limit the ratio of proton to electron luminosities to Lp<~0.15 Le.

The MAGIC stereoscopic system collected 69 hours of Crab Nebula data between October 2009 and April 2011. Analysis of this data sample using the latest improvements in the MAGIC stereoscopic software provided an unprecedented precision of spectral and night-by-night light curve determination at gamma rays. We derived a differential spectrum with a single instrument from 50 GeV up to almost 30 TeV with 5 bins per energy decade. At low energies, MAGIC results, combined with Fermi-LAT data, show a flat and broad Inverse Compton peak. The overall fit to the data between 1 GeV and 30 TeV is not well described by a log-parabola function. We find that a modified log-parabola function with an exponent of 2.5 instead of 2 provides a good description of the data ($\chi^2=35/26$). Using systematic uncertainties of red the MAGIC and Fermi-LAT measurements we determine the position of the Inverse Compton peak to be at (53 $\pm$ 3stat + 31syst -13syst) GeV, which is the most precise estimation up to date and is dominated by the systematic effects. There is no hint of the integral flux variability on daily scales at energies above 300 GeV when systematic uncertainties are included in the flux measurement. We consider three state- of-the-art theoretical models to describe the overall spectral energy distribution of the Crab Nebula. The constant B-field model cannot satisfactorily reproduce the VHE spectral measurements presented in this work, having particular difficulty reproducing the broadness of the observed IC peak. Most probably this implies that the assumption of the homogeneity of the magnetic field inside the nebula is incorrect. On the other hand, the time-dependent 1D spectral model provides a good fit of the new VHE results when considering a 80 {\mu}G magnetic field. However, it fails to match the data when including the morphology of the nebula at lower wavelengths.

We present constraints derived from a search of four years of IceCube data for a prompt neutrino flux from gamma-ray bursts (GRBs). A single low-significance neutrino, compatible with the atmospheric neutrino background, was found in coincidence with one of the 506 observed bursts. Although GRBs have been proposed as candidate sources for ultra-high energy cosmic rays, our limits on the neutrino flux disfavor much of the parameter space for the latest models. We also find that no more than $\sim1\%$ of the recently observed astrophysical neutrino flux consists of prompt emission from GRBs that are potentially observable by existing satellites.

A flare from the TeV blazar Mrk 421, occurring in March 2010, was observed for 13 consecutive days from radio to very high energy (VHE, E > 100 GeV) gamma-rays with MAGIC, VERITAS, Whipple, FermiLAT, MAXI, RXTE, Swift, GASP-WEBT, and several optical and radio telescopes. We model the day-scale SEDs with one-zone and two-zone synchrotron self-Compton (SSC) models, investigate the physical parameters, and evaluate whether the observed broadband SED variability can be associated to variations in the relativistic particle population. Flux variability was remarkable in the X-ray and VHE bands while it was minor or not significant in the other bands. The one-zone SSC model can describe reasonably well the SED of each day for the 13 consecutive days. This flaring activity is also very well described by a two-zone SSC model, where one zone is responsible for the quiescent emission while the other smaller zone, which is spatially separated from the first one, contributes to the daily-variable emission occurring in X-rays and VHE gamma-rays. Both the one-zone SSC and the two-zone SSC models can describe the daily SEDs via the variation of only four or five model parameters, under the hypothesis that the variability is associated mostly to the underlying particle population. This shows that the particle acceleration and cooling mechanism producing the radiating particles could be the main one responsible for the broadband SED variations during the flaring episodes in blazars. The two-zone SSC model provides a better agreement to the observed SED at the narrow peaks of the low- and high-energy bumps during the highest activity, although the reported one-zone SSC model could be further improved by the variation of the parameters related to the emitting region itself ($\delta$, $B$ and $R$), in addition to the parameters related to the particle population.

In this paper searches for flaring astrophysical neutrino sources and sources with periodic emission with the IceCube neutrino telescope are presented. In contrast to time integrated searches, where steady emission is assumed, the analyses presented here look for a time dependent signal of neutrinos using the information from the neutrino arrival times to enhance the discovery potential. A search was performed for correlations between neutrino arrival times and directions as well as neutrino emission following time dependent lightcurves, sporadic emission or periodicities of candidate sources. These include active galactic nuclei, soft $\gamma$-ray repeaters, supernova remnants hosting pulsars, micro-quasars and X-ray binaries. The work presented here updates and extends previously published results to a longer period that covers four years of data from 2008 April 5 to 2012 May 16 including the first year of operation of the completed 86-string detector. The analyses did not find any significant time dependent point sources of neutrinos and the results were used to set upper limits on the neutrino flux from source candidates.

The MAGIC telescopes are two Imaging Atmospheric Cherenkov Telescopes (IACTs) located on the Canary island of La Palma. The telescopes are designed to measure Cherenkov light from air showers initiated by gamma rays in the energy regime from around 50 GeV to more than 50 TeV. The two telescopes were built in 2004 and 2009, respectively, with different cameras, triggers and readout systems. In the years 2011-2012 the MAGIC collaboration undertook a major upgrade to make the stereoscopic system uniform, improving its overall performance and easing its maintenance. In particular, the camera, the receivers and the trigger of the first telescope were replaced and the readout of the two telescopes was upgraded. This paper (Part I) describes the details of the upgrade as well as the basic performance parameters of MAGIC such as raw data treatment, dead time of the system, linearity in the electronic chain and sources of noise. In Part II, we describe the physics performance of the upgraded system.

The high-frequency-peaked BL Lac (HBL) 1ES 0806+524 (z = 0.138) was discovered in VHE $\gamma$ rays in 2008. Until now, the broad-band spectrum of 1ES 0806+524 has been only poorly characterized, in particular at high energies. We analysed multiwavelength observations from $\gamma$ rays to radio performed from 2011 January to March, which were triggered by the high activity detected at optical frequencies. These observations constitute the most precise determination of the broad-band emission of 1ES 0806+524 to date. The stereoscopic MAGIC observations yielded a $\gamma$-ray signal above 250 GeV of $(3.7 \pm 0.7)$ per cent of the Crab Nebula flux with a statistical significance of 9.9 $\sigma$. The multiwavelength observations showed significant variability in essentially all energy bands, including a VHE $\gamma$-ray flare that lasted less than one night, which provided unprecedented evidence for short-term variability in 1ES 0806+524. The spectrum of this flare is well described by a power law with a photon index of $2.97 \pm 0.29$ between $\sim$150 GeV and 1 TeV and an integral flux of $(9.3 \pm 1.9)$ per cent of the Crab Nebula flux above 250 GeV. The spectrum during the non-flaring VHE activity is compatible with the only available VHE observation performed in 2008 with VERITAS when the source was in a low optical state. The broad-band spectral energy distribution can be described with a one-zone Synchrotron Self Compton model with parameters typical for HBLs, indicating that 1ES 0806+524 is not substantially different from the HBLs previously detected.

PG 1553+113 is a very-high-energy (VHE, $E>100\,\mathrm{GeV}$) $\gamma$-ray emitter classified as a BL Lac object. Its redshift is constrained by intergalactic absorption lines in the range $0.4<z<0.58$. The MAGIC telescopes have monitored the source's activity since 2005. In early 2012, PG 1553+113 was found in a high-state, and later, in April of the same year, the source reached its highest VHE flux state detected so far. Simultaneous observations carried out in X-rays during 2012 April show similar flaring behaviour. In contrast, the $\gamma$-ray flux at $E<100\,\mathrm{GeV}$ observed by Fermi-LAT is compatible with steady emission. In this paper, a detailed study of the flaring state is presented. The VHE spectrum shows clear curvature, being well fitted either by a power law with an exponential cut-off or by a log-parabola. A simple power-law fit hypothesis for the observed shape of the PG 1553+113 VHE $\gamma$-ray spectrum is rejected with a high significance (fit probability P=2.6 $\times 10^{-6}$). The observed curvature is compatible with the extragalactic background light (EBL) imprint predicted by current generation EBL models assuming a redshift $z\sim0.4$. New constraints on the redshift are derived from the VHE spectrum. These constraints are compatible with previous limits and suggest that the source is most likely located around the optical lower limit, $z=0.4$, based on the detection of Ly$\alpha$ absorption. Finally, we find that the synchrotron self-Compton (SSC) model gives a satisfactory description of the observed multi-wavelength spectral energy distribution during the flare.

MAGIC is a system of two Imaging Atmospheric Cherenkov Telescopes located in the Canary island of La Palma, Spain. During summer 2011 and 2012 it underwent a series of upgrades, involving the exchange of the MAGIC-I camera and its trigger system, as well as the upgrade of the readout system of both telescopes. We use observations of the Crab Nebula taken at low and medium zenith angles to assess the key performance parameters of the MAGIC stereo system. For low zenith angle observations, the standard trigger threshold of the MAGIC telescopes is ~50GeV. The integral sensitivity for point-like sources with Crab Nebula-like spectrum above 220GeV is (0.66+/-0.03)% of Crab Nebula flux in 50 h of observations. The angular resolution, defined as the sigma of a 2-dimensional Gaussian distribution, at those energies is < 0.07 degree, while the energy resolution is 16%. We also re-evaluate the effect of the systematic uncertainty on the data taken with the MAGIC telescopes after the upgrade. We estimate that the systematic uncertainties can be divided in the following components: < 15% in energy scale, 11-18% in flux normalization and +/-0.15 for the energy spectrum power-law slope.

We present the development and application of a generic analysis scheme for the measurement of neutrino spectra with the IceCube detector. This scheme is based on regularized unfolding, preceded by an event selection which uses a Minimum Redundancy Maximum Relevance algorithm to select the relevant variables and a Random Forest for the classification of events. The analysis has been developed using IceCube data from the 59-string configuration of the detector. 27,771 neutrino candidates were detected in 346 days of livetime. A rejection of 99.9999% of the atmospheric muon background is achieved. The energy spectrum of the atmospheric neutrino flux is obtained using the TRUEE unfolding program. The unfolded spectrum of atmospheric muon neutrinos covers an energy range from 100 GeV to 1 PeV. Compared to the previous measurement using the detector in the 40-string configuration, the analysis presented here, extends the upper end of the atmospheric neutrino spectrum by more than a factor of two, reaching an energy region that has not been previously accessed by spectral measurements.

A diffuse flux of astrophysical neutrinos above $100\,\mathrm{TeV}$ has been observed at the IceCube Neutrino Observatory. Here we extend this analysis to probe the astrophysical flux down to $35\,\mathrm{TeV}$ and analyze its flavor composition by classifying events as showers or tracks. Taking advantage of lower atmospheric backgrounds for shower-like events, we obtain a shower-biased sample containing 129 showers and 8 tracks collected in three years from 2010 to 2013. We demonstrate consistency with the $(f_e:f_{\mu}:f_\tau)_\oplus\approx(1:1:1)_\oplus$ flavor ratio at Earth commonly expected from the averaged oscillations of neutrinos produced by pion decay in distant astrophysical sources. Limits are placed on non-standard flavor compositions that cannot be produced by averaged neutrino oscillations but could arise in exotic physics scenarios. A maximally track-like composition of $(0:1:0)_\oplus$ is excluded at $3.3\sigma$, and a purely shower-like composition of $(1:0:0)_\oplus$ is excluded at $2.3\sigma$.

We performed a 4.5-month multi-instrument campaign (from radio to VHE gamma rays) on Mrk421 between January 2009 and June 2009, which included VLBA, F-GAMMA, GASP-WEBT, Swift, RXTE, Fermi-LAT, MAGIC, and Whipple, among other instruments and collaborations. Mrk421 was found in its typical (non-flaring) activity state, with a VHE flux of about half that of the Crab Nebula, yet the light curves show significant variability at all wavelengths, the highest variability being in the X-rays. We determined the power spectral densities (PSD) at most wavelengths and found that all PSDs can be described by power-laws without a break, and with indices consistent with pink/red-noise behavior. We observed a harder-when-brighter behavior in the X-ray spectra and measured a positive correlation between VHE and X-ray fluxes with zero time lag. Such characteristics have been reported many times during flaring activity, but here they are reported for the first time in the non-flaring state. We also observed an overall anti-correlation between optical/UV and X-rays extending over the duration of the campaign. The harder-when-brighter behavior in the X-ray spectra and the measured positive X-ray/VHE correlation during the 2009 multi-wavelength campaign suggests that the physical processes dominating the emission during non-flaring states have similarities with those occurring during flaring activity. In particular, this observation supports leptonic scenarios as being responsible for the emission of Mrk421 during non-flaring activity. Such a temporally extended X-ray/VHE correlation is not driven by any single flaring event, and hence is difficult to explain within the standard hadronic scenarios. The highest variability is observed in the X-ray band, which, within the one-zone synchrotron self-Compton scenario, indicates that the electron energy distribution is most variable at the highest energies.

Context: MWC 656 has recently been established as the first observationally detected high-mass X-ray binary system containing a Be star and a black hole (BH). The system has been associated with a gamma-ray flaring event detected by the AGILE satellite in July 2010. Aims: Our aim is to evaluate if the MWC 656 gamma-ray emission extends to very high energy (VHE > 100 GeV) gamma rays. Methods. We have observed MWC 656 with the MAGIC telescopes for $\sim$23 hours during two observation periods: between May and June 2012 and June 2013. During the last period, observations were performed contemporaneously with X-ray (XMM-Newton) and optical (STELLA) instruments. Results: We have not detected the MWC 656 binary system at TeV energies with the MAGIC Telescopes in either of the two campaigns carried out. Upper limits (ULs) to the integral flux above 300 GeV have been set, as well as differential ULs at a level of $\sim$5% of the Crab Nebula flux. The results obtained from the MAGIC observations do not support persistent emission of very high energy gamma rays from this system at a level of 2.4% the Crab flux.

Recently, IceCube found evidence for a diffuse signal of astrophysical neutrinos in an energy range of $60\,\mathrm{TeV}$ to the $\mathrm{PeV}$-scale. The origin of those events, being a key to understanding the origin of cosmic rays, is still an unsolved question. So far, analyses have not succeeded to resolve the diffuse signal into point-like sources. Searches including a maximum-likelihood-ratio test, based on the reconstructed directions and energies of the detected down- and up-going neutrino candidates, were also performed on IceCube data leading to the exclusion of bright point sources. In this paper, we present two methods to search for faint neutrino point sources in three years of IceCube data, taken between 2008 and 2011. The first method is an autocorrelation test, applied separately to the northern and southern sky. The second method is a multipole analysis, which expands the measured data in the northern hemisphere into spherical harmonics and uses the resulting expansion coefficients to separate signal from background. With both methods, the results are consistent with the background expectation with a slightly more sparse spatial distribution, corresponding to an underfluctuation. Depending on the assumed number of sources, the resulting upper limit on the flux per source in the northern hemisphere for an $E^{-2}$ energy spectrum ranges from $1.5 \cdot 10^{-8}\,\mathrm{GeV}/(\mathrm{cm}^2 \mathrm{s})$, in the case of one assumed source, to $4 \cdot 10^{-10} \,\mathrm{GeV}/(\mathrm{cm}^2 \mathrm{s})$, in the case of $3500$ assumed sources.