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.

Supermassive black holes with masses of millions to billions of solar masses are commonly found in the centers of galaxies. Astronomers seek to image jet formation using radio interferometry, but still suffer from insufficient angular resolution. An alternative method to resolve small structures is to measure the time variability of their emission. Here, we report on gamma-ray observations of the radio galaxy IC 310 obtained with the MAGIC telescopes revealing variability with doubling time scales faster than 4.8 min. Causality constrains the size of the emission region to be smaller than 20\% of the gravitational radius of its central black hole. We suggest that the emission is associated with pulsar-like particle acceleration by the electric field across a magnetospheric gap at the base of the radio jet.

HESS J1857+026 is an extended TeV gamma-ray source that was discovered by H.E.S.S. as part of its Galactic plane survey. Given its broadband spectral energy distribution and its spatial coincidence with the young energetic pulsar PSR J1856+0245, the source has been put forward as a pulsar wind nebula (PWN) candidate. MAGIC has performed follow-up observations aimed at mapping the source down to energies approaching 100 GeV in order to better understand its complex morphology. HESS J1857+026 was observed by MAGIC in 2010, yielding 29 hours of good quality stereoscopic data that allowed us to map the source region in two separate ranges of energy. We detected very-high-energy gamma-ray emission from HESS J1857+026 with a significance of $12 \sigma$ above $150$ GeV. The differential energy spectrum between $100$ GeV and $13$ TeV is well described by a power law function $dN/dE = N_0(E/1\textrm{TeV})^{-\Gamma}$ with $N_0 = (5.37 \pm0.44_{stat} \pm1.5_{sys}) \times 10^{-12} (\textrm{TeV}^{-1} \textrm{cm}^{-2}$ $\textrm{ s}^{-1})$ and $\Gamma = 2.16\pm0.07_{stat} \pm0.15_{sys}$, which bridges the gap between the GeV emission measured by Fermi-LAT and the multi-TeV emission measured by H.E.S.S.. In addition, we present a detailed analysis of the energy-dependent morphology of this region. We couple these results with archival multi-wavelength data and outline evidence in favor of a two-source scenario, whereby one source is associated with a PWN, while the other could be linked with a molecular cloud complex containing an HII region and a possible gas cavity.

It has been claimed that the nova-like cataclysmic variable (CV) AE Aquarii (AE Aqr) is a very-high-energy (VHE, $E>$100 GeV) source both on observational and theoretical grounds. We aim to search for VHE gamma-ray emission from AE Aqr during different states of the source at several wavelengths to confirm or rule out previous claims of detection of gamma-ray emission from this object. We report on observations of AE Aqr performed by MAGIC. The source was observed during 12 hours as part of a multiwavelength campaign carried out between May and June 2012 covering the optical, X-ray, and gamma-ray ranges. Besides MAGIC, the other facilities involved were the KVA, Skinakas, and Vidojevica telescopes in the optical and Swift in X-rays. We calculated integral upper limits coincident with different states of the source in the optical. We computed upper limits to the pulsed emission limiting the signal region to 30% of the phaseogram and we also searched for pulsed emission at different frequencies applying the Rayleigh test. AE Aqr was not detected at VHE energies during the multiwavelength campaign. We establish integral upper limits at the 95\% confidence level for the steady emission assuming the differential flux proportional to a power-law function d\phi/dE \propto E^{-Gamma}, with a Crab-like photon spectral index of Gamma=2.6. The upper limit above 200 GeV is 6.4\times10^{-12} cm^{-2}s^{-1} and above 1 TeV is 7.4\times10^{-13} cm^{-2}s^{-1}. We obtained an upper limit for the pulsed emission of 2.6\times10^{-12} cm^{-2}s^{-1} for energies above 200 GeV. Applying the Rayleigh test for pulsed emission at different frequencies we did not find any significant signal. Our results indicate that AE Aqr is not a VHE gamma-ray emitter at the level of emission previously claimed. We have established the most constraining upper limits for the VHE gamma-ray emission of AE Aqr.

The pulsar wind nebula (PWN) 3C 58 is one of the historical very-high-energy (VHE; E>100 GeV) gamma-ray source candidates. It is energized by one of the highest spin-down power pulsars known (5% of Crab pulsar) and it has been compared to the Crab Nebula due to their morphological similarities. This object was previously observed by imaging atmospheric Cherenkov telescopes (Whipple, VERITAS and MAGIC), although not detected, with an upper limit of 2.4% Crab Unit (C.U.) at VHE. It was detected by Fermi-LAT with a spectrum extending beyond 100 GeV. We analyzed 81 hours of 3C 58 data taken with the MAGIC telescopes and we detected VHE gamma-ray emission with a significance of 5.7 sigma and an integral flux of 0.65% C.U. above 1 TeV. The differential energy spectrum between 400 GeV and 10 TeV is well described by a power-law function d\phi/dE=f_0(E/1TeV)^{-Gamma} with f_0=(2.0\pm0.4_{stat}\pm0.6_{sys})\times10^{-13}cm^{-2}s^{-1}TeV^{-1} and Gamma=2.4\pm0.2_{stat}\pm0.2_{sys}. The skymap is compatible with an unresolved source. We report the first significant detection of PWN 3C 58 at TeV energies. According to our results 3C 58 is the least luminous VHE gamma-ray PWN ever detected at VHE and the one with the lowest flux at VHE to date. We compare our results with the expectations of time-dependent models in which electrons up-scatter photon fields. The best representation favors a distance to the PWN of 2 kpc and Far Infrared (FIR) comparable to CMB photon fields. If we consider an unexpectedly high FIR density, the data can also be reproduced by models assuming a 3.2 kpc distance. A low magnetic field, far from equipartition, is required to explain the VHE data. Hadronic contribution from the hosting supernova remnant (SNR) requires unrealistic energy budget given the density of the medium, disfavoring cosmic ray acceleration in the SNR as origin of the VHE gamma-ray emission.

We present a study of the very high energy (VHE; E > 100 GeV) gamma-ray emission of the blazar PKS 1424+240 observed with the MAGIC telescopes. The primary aim of this paper is the multiwavelength spectral characterization and modeling of this blazar, which is made particularly interesting by the recent discovery of a lower limit of its redshift of z > 0.6 and makes it a promising candidate to be the most distant VHE source. The source has been observed with the MAGIC telescopes in VHE gamma rays for a total observation time of ~33.6 h from 2009 to 2011. The source was marginally detected in VHE gamma rays during 2009 and 2010, and later, the detection was confirmed during an optical outburst in 2011. The combined significance of the stacked sample is ~7.2 sigma. The differential spectra measured during the different campaigns can be described by steep power laws with the indices ranging from 3.5 +/- 1.2 to 5.0 +/- 1.7. The MAGIC spectra corrected for the absorption due to the extragalactic background light connect smoothly, within systematic errors, with the mean spectrum in 2009-2011 observed at lower energies by the Fermi-LAT. The absorption-corrected MAGIC spectrum is flat with no apparent turn down up to 400 GeV. The multiwavelength light curve shows increasing flux in radio and optical bands that could point to a common origin from the same region of the jet. The large separation between the two peaks of the constructed non-simultaneous spectral energy distribution also requires an extremely high Doppler factor if an one zone synchrotron self-Compton model is applied. We find that a two-component synchrotron self-Compton model describes the spectral energy distribution of the source well, if the source is located at z~0.6.

We present the results of stereoscopic observations of the satellite galaxy Segue 1 with the MAGIC Telescopes, carried out between 2011 and 2013. With almost 160 hours of good-quality data, this is the deepest observational campaign on any dwarf galaxy performed so far in the very high energy range of the electromagnetic spectrum. We search this large data sample for signals of dark matter particles in the mass range between 100 GeV and 20 TeV. For this we use the full likelihood analysis method, which provides optimal sensitivity to characteristic gamma-ray spectral features, like those expected from dark matter annihilation or decay. In particular, we focus our search on gamma-rays produced from different final state Standard Model particles, annihilation with internal bremsstrahlung, monochromatic lines and box-shaped signals. Our results represent the most stringent constraints to the annihilation cross-section or decay lifetime obtained from observations of satellite galaxies, for masses above few hundred GeV. In particular, our strongest limit (95% confidence level) corresponds to a ~500 GeV dark matter particle annihilating into tau+tau-, and is of order <sigma v> ~ 1.2x10^{-24} cm^3 s^{-1} - a factor ~40 above the <sigma v> thermal value.