We present 10 years of R-band monitoring data of 31 northern blazars which were either detected at very high energy (VHE) gamma rays or listed as potential VHE gamma-ray emitters. The data comprise 11820 photometric data points in the R-band obtained in 2002-2012. We analyze the light curves by determining their power spectral density (PSD) slopes assuming a power-law dependence with a single slope $\beta$ and a Gaussian probability density function (PDF). We use the multiple fragments variance function (MFVF) combined with a forward-casting approach and likelihood analysis to determine the slopes and perform extensive simulations to estimate the uncertainties of the derived slopes. We also look for periodic variations via Fourier analysis and quantify the false alarm probability through a large number of simulations. Comparing the obtained PSD slopes to values in the literature, we find the slopes in the radio band to be steeper than those in the optical and gamma rays. Our periodicity search yielded one target, Mrk 421, with a significant (p<5%) period. Finding one significant period among 31 targets is consistent with the expected false alarm rate, but the period found in Mrk~421 is very strong and deserves further consideration}.

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.

[Abridged] Context. To construct and interpret the spectral energy distribution (SED) of BL Lacertae objects, simultaneous broad-band observations are mandatory. Aims. We present the results of a dedicated multi-wavelength study of the high-frequency peaked BL Lacertae (HBL) object and known TeV emitter 1ES 2344+514 by means of a pre-organised campaign. Methods. The observations were conducted during simultaneous visibility windows of MAGIC and AGILE in late 2008. The measurements were complemented by Mets\"ahovi, RATAN-600, KVA+Tuorla, Swift and VLBA pointings. Additional coverage was provided by the ongoing long-term F-GAMMA and MOJAVE programs, the OVRO 40-m and CrAO telescopes as well as the Fermi satellite. The obtained SEDs are modelled using a one-zone as well as a self-consistent two-zone synchrotron self-Compton model. Results. 1ES 2344+514 was found at very low flux states in both X-rays and very high energy gamma rays. Variability was detected in the low frequency radio and X-ray bands only, where for the latter a small flare was observed. The X-ray flare was possibly caused by shock acceleration characterised by similar cooling and acceleration time scales. MOJAVE VLBA monitoring reveals a static jet whose components are stable over time scales of eleven years, contrary to previous findings. There appears to be no significant correlation between the 15 GHz and R-band monitoring light curves. The observations presented here constitute the first multi-wavelength campaign on 1ES 2344+514 from radio to VHE energies and one of the few simultaneous SEDs during low activity states. The quasi-simultaneous Fermi-LAT data poses some challenges for SED modelling, but in general the SEDs are described well by both applied models. The resulting parameters are typical for TeV emitting HBLs. Consequently it remains unclear whether a so-called quiescent state was found in this campaign.

Context: The blazar Markarian 421 is one of the brightest TeV gamma-ray sources of the northern sky. From December 2007 until June 2008 it was intensively observed in the very high energy (VHE, E > 100 GeV) band by the single-dish Major Atmospheric Gamma-ray Imaging Cherenkov telescope (MAGIC-I). Aims: We aimed to measure the physical parameters of the emitting region of the blazar jet during active states. Methods: We performed a dense monitoring of the source in VHE with MAGIC-I, and also collected complementary data in soft X-rays and optical-UV bands; then, we modeled the spectral energy distributions (SED) derived from simultaneous multi-wavelength data within the synchrotron self--compton (SSC) framework. Results: The source showed intense and prolonged gamma-ray activity during the whole period, with integral fluxes (E > 200 GeV) seldom below the level of the Crab Nebula, and up to 3.6 times this value. Eight datasets of simultaneous optical-UV (KVA, Swift/UVOT), soft X-ray (Swift/XRT) and MAGIC-I VHE data were obtained during different outburst phases. The data constrain the physical parameters of the jet, once the spectral energy distributions obtained are interpreted within the framework of a single-zone SSC leptonic model. Conclusions: The main outcome of the study is that within the homogeneous model high Doppler factors (40 <= delta <= 80) are needed to reproduce the observed SED; but this model cannot explain the observed short time-scale variability, while it can be argued that inhomogeneous models could allow for less extreme Doppler factors, more intense magnetic fields and shorter electron cooling times compatible with hour or sub-hour scale variability.

The quasar-type blazar 4C 38.41 (B3 1633+382) experienced a large outburst in 2011, which was detected throughout the entire electromagnetic spectrum. We present the results of low-energy multifrequency monitoring by the GASP project of the WEBT consortium and collaborators, as well as those of spectropolarimetric/spectrophotometric monitoring at the Steward Observatory. We also analyse high-energy observations of the Swift and Fermi satellites. In the optical-UV band, several results indicate that there is a contribution from a QSO-like emission component, in addition to both variable and polarised jet emission. The unpolarised emission component is likely thermal radiation from the accretion disc that dilutes the jet polarisation. We estimate its brightness to be R(QSO) ~ 17.85 - 18 and derive the intrinsic jet polarisation degree. We find no clear correlation between the optical and radio light curves, while the correlation between the optical and \gamma-ray flux apparently fades in time, likely because of an increasing optical to \gamma-ray flux ratio. As suggested for other blazars, the long-term variability of 4C 38.41 can be interpreted in terms of an inhomogeneous bent jet, where different emitting regions can change their alignment with respect to the line of sight, leading to variations in the Doppler factor \delta. Under the hypothesis that in the period 2008-2011 all the \gamma-ray and optical variability on a one-week timescale were due to changes in \delta, this would range between ~ 7 and ~ 21. If the variability were caused by changes in the viewing angle \theta\ only, then \theta\ would go from ~ 2.6 degr to ~ 5 degr. Variations in the viewing angle would also account for the dependence of the polarisation degree on the source brightness in the framework of a shock-in-jet model.

We present the results of a long M87 monitoring campaign in very high energy $\gamma$-rays with the MAGIC-I Cherenkov telescope. We aim to model the persistent non-thermal jet emission by monitoring and characterizing the very high energy $\gamma$-ray emission of M87 during a low state. A total of 150\,h of data were taken between 2005 and 2007 with the single MAGIC-I telescope, out of which 128.6\,h survived the data quality selection. We also collected data in the X-ray and \textit{Fermi}--LAT bands from the literature (partially contemporaneous). No flaring activity was found during the campaign. The source was found to be in a persistent low-emission state, which was at a confidence level of $7\sigma$. We present the spectrum between 100\,GeV and 2\,TeV, which is consistent with a simple power law with a photon index $\Gamma=2.21\pm0.21$ and a flux normalization at 300\,GeV of $(7.7\pm1.3) \times 10^{-8}$ TeV$^{-1}$ s$^{-1}$ m$^{-2}$. The extrapolation of the MAGIC spectrum into the GeV energy range matches the previously published \textit{Fermi}--LAT spectrum well, covering a combined energy range of four orders of magnitude with the same spectral index. We model the broad band energy spectrum with a spine layer model, which can satisfactorily describe our data.

Abridged: The giant radio galaxy M 87 with its proximity, famous jet, and very massive black hole provides a unique opportunity to investigate the origin of very high energy (VHE; E>100 GeV) gamma-ray emission generated in relativistic outflows and the surroundings of super-massive black holes. M 87 has been established as a VHE gamma-ray emitter since 2006. The VHE gamma-ray emission displays strong variability on timescales as short as a day. In this paper, results from a joint VHE monitoring campaign on M 87 by the MAGIC and VERITAS instruments in 2010 are reported. During the campaign, a flare at VHE was detected triggering further observations at VHE (H.E.S.S.), X-rays (Chandra), and radio (43 GHz VLBA). The excellent sampling of the VHE gamma-ray light curve enables one to derive a precise temporal characterization of the flare: the single, isolated flare is well described by a two-sided exponential function with significantly different flux rise and decay times. While the overall variability pattern of the 2010 flare appears somewhat different from that of previous VHE flares in 2005 and 2008, they share very similar timescales (~day), peak fluxes (Phi(>0.35 TeV) ~= (1-3) x 10^-11 ph cm^-2 s^-1), and VHE spectra. 43 GHz VLBA radio observations of the inner jet regions indicate no enhanced flux in 2010 in contrast to observations in 2008, where an increase of the radio flux of the innermost core regions coincided with a VHE flare. On the other hand, Chandra X-ray observations taken ~3 days after the peak of the VHE gamma-ray emission reveal an enhanced flux from the core. The long-term (2001-2010) multi-wavelength light curve of M 87, spanning from radio to VHE and including data from HST, LT, VLA and EVN, is used to further investigate the origin of the VHE gamma-ray emission. No unique, common MWL signature of the three VHE flares has been identified.

We study the non-thermal jet emission of the BL Lac object B3 2247+381 during a high optical state. The MAGIC telescopes observed the source during 13 nights between September 30th and October 30th 2010, collecting a total of 14.2 hours of good quality very high energy (VHE) $\gamma$-ray data. Simultaneous multiwavelength data was obtained with X-ray observations by the Swift satellite and optical R-band observations at the KVA-telescope. We also use high energy $\gamma$-ray (HE, 0.1 GeV-100 GeV) data from the Fermi satellite. The BL Lac object B3 2247+381 (z=0.119) was detected, for the first time, at VHE $\gamma$-rays at a statistical significance of 5.6 $\sigma$. A soft VHE spectrum with a photon index of -3.2 $\pm$ 0.6 was determined. No significant short term flux variations were found. We model the spectral energy distribution using a one-zone SSC-model, which can successfully describe our data.

We present the results of five years (2005-2009) of MAGIC observations of the BL Lac object PG 1553+113 at very high energies (VHEs, E > 100 GeV). Power law fits of the individual years are compatible with a steady mean photon index \Gamma = 4.27 $\pm$ 0.14. In the last three years of data, the flux level above 150 GeV shows a clear variability (probability of constant flux < 0.001%). The flux variations are modest, lying in the range from 4% to 11% of the Crab Nebula flux. Simultaneous optical data also show only modest variability that seems to be correlated with VHE gamma ray variability. We also performed a temporal analysis of (all available) simultaneous Fermi/LAT data of PG 1553+113 above 1 GeV, which reveals hints of variability in the 2008-2009 sample. Finally, we present a combination of the mean spectrum measured at very high energies with archival data available for other wavelengths. The mean spectral energy distribution can be modeled with a one-zone Synchrotron Self Compton (SSC) model, which gives the main physical parameters governing the VHE emission in the blazar jet.

We present very high energy (VHE, E > 100 GeV) {\gamma}-ray observations of the {\gamma}-ray binary system LS I+61 303 obtained with the MAGIC stereo system between 2009 October and 2010 January. We detect a 6.3{\sigma} {\gamma}-ray signal above 400 GeV in the combined data set. The integral flux above an energy of 300 GeV is F(E>300 GeV)=(1.4 +- 0.3stat +- 0.4syst) * 10^{-12} cm^{-2} s^{-1}, which corresponds to about 1.3% of the Crab Nebula flux in the same energy range. The orbit-averaged flux of LS I +61 303 in the orbital phase interval 0.6--0.7, where a maximum of the TeV flux is expected, is lower by almost an order of magnitude compared to our previous measurements between 2005 September and 2008 January. This provides evidence for a new low emission state in LS I +61 303. We find that the change to the low flux state cannot be solely explained by an increase of photon-photon absorption around the compact star.

We report the results of the observation of the nearby satellite galaxy Segue 1 performed by the MAGIC-I ground-based gamma-ray telescope between November 2008 and March 2009 for a total of 43.2 hours. No significant gamma-ray emission was found above the background. Differential upper limits on the gamma-ray flux are derived assuming various power-law slopes for the possible emission spectrum. Integral upper limits are also calculated for several power-law spectra and for different energy thresholds. The values are of the order of 10^{-11} ph cm^{-2}$ s^{-1} above 100 GeV and 10^{-12} ph cm^{-2} s^{-1} above 200 GeV. Segue 1 is currently considered one of the most interesting targets for indirect dark matter searches. In these terms, the upper limits have been also interpreted in the context of annihilating dark matter particles. For such purpose, we performed a grid scan over a reasonable portion of the parameter space for the minimal SuperGravity model and computed the flux upper limit for each point separately, taking fully into account the peculiar spectral features of each model. We found that in order to match the experimental upper limits with the model predictions, a minimum flux boost of 10^{3} is required, and that the upper limits are quite dependent on the shape of the gamma-ray energy spectrum predicted by each specific model. Finally we compared the upper limits with the predictions of some dark matter models able to explain the PAMELA rise in the positron ratio, finding that Segue 1 data are in tension with the dark matter explanation of the PAMELA spectrum in the case of a dark matter candidate annihilating into tau+tau-. A complete exclusion however is not possible due to the uncertainties in the Segue 1 astrophysical factor.

Since 2005, the blazar 3C 454.3 has shown remarkable flaring activity at all frequencies, and during the last four years it has exhibited more than one gamma-ray flare per year, becoming the most active gamma-ray blazar in the sky. We present for the first time the multi-wavelength AGILE, SWIFT, INTEGRAL, and GASP-WEBT data collected in order to explain the extraordinary gamma-ray flare of 3C 454.3 which occurred in November 2010. On 2010 November 20 (MJD 55520), 3C 454.3 reached a peak flux (E>100 MeV) of F_gamma(p) = (6.8+-1.0)E-5 ph/cm2/s on a time scale of about 12 hours, more than a factor of 6 higher than the flux of the brightest steady gamma-ray source, the Vela pulsar, and more than a factor of 3 brighter than its previous super-flare on 2009 December 2-3. The multi-wavelength data make a thorough study of the present event possible: the comparison with the previous outbursts indicates a close similarity to the one that occurred in 2009. By comparing the broadband emission before, during, and after the gamma-ray flare, we find that the radio, optical and X-ray emission varies within a factor 2-3, whereas the gamma-ray flux by a factor of 10. This remarkable behavior is modeled by an external Compton component driven by a substantial local enhancement of soft seed photons.

We report on TeV gamma-ray observations of the blazar Mrk 421 (redshift of 0.031) with the VERITAS observatory and the Whipple 10m Cherenkov telescope. The excellent sensitivity of VERITAS allowed us to sample the TeV gamma-ray fluxes and energy spectra with unprecedented accuracy where Mrk 421 was detected in each of the pointings. A total of 47.3 hrs of VERITAS and 96 hrs of Whipple 10m data were acquired between January 2006 and June 2008. We present the results of a study of the TeV gamma-ray energy spectra as a function of time, and for different flux levels. On May 2nd and 3rd, 2008, bright TeV gamma-ray flares were detected with fluxes reaching the level of 10 Crab. The TeV gamma-ray data were complemented with radio, optical, and X-ray observations, with flux variability found in all bands except for the radio waveband. The combination of the RXTE and Swift X-ray data reveal spectral hardening with increasing flux levels, often correlated with an increase of the source activity in TeV gamma-rays. Contemporaneous spectral energy distributions were generated for 18 nights, each of which are reasonably described by a one-zone SSC model.

The acceleration of particles up to GeV or higher energies in microquasars has been the subject of considerable theoretical and observational efforts in the past few years. Sco X-1 is a microquasar from which evidence of highly energetic particles in the jet has been found when it is in the so-called Horizontal Branch (HB), a state when the radio and hard X-ray fluxes are higher and a powerful relativistic jet is present. Here we present the first very high energy gamma-ray observations of Sco X-1 obtained with the MAGIC telescopes. An analysis of the whole dataset does not yield a significant signal, with 95% CL flux upper limits above 300 GeV at the level of 2.4x10^{-12} ph/cm^2/s. Simultaneous RXTE observations were conducted to search for TeV emission during particular X-ray states of the source. A selection of the gamma-ray data obtained during the HB based on the X-ray colors did not yield a signal either, with an upper limit of 3.4x10^{-12} ph/cm^2/s. These upper limits place a constraint on the maximum TeV luminosity to non-thermal X-ray luminosity of L_{VHE}/L_{ntX}<0.02, that can be related to a maximum TeV luminosity to jet power ratio of L_{VHE}/L_{j}<10^{-3}. Our upper limits indicate that the underlying high-energy emission physics in Sco X-1 must be inherently different from that of the hitherto detected gamma-ray binaries.

We report on the observation of the region around supernova remnant G65.1+0.6 with the stand-alone MAGIC-I telescope. This region hosts the two bright GeV gamma-ray sources 1FGL J1954.3+2836 and 1FGL J1958.6+2845. They are identified as GeV pulsars and both have a possible counterpart detected at about 35 TeV by the Milagro observatory. MAGIC collected 25.5 hours of good quality data, and found no significant emission in the range around 1 TeV. We therefore report differential flux upper limits, assuming the emission to be point-like (<0.1 deg) or within a radius of 0.3 deg. In the point-like scenario, the flux limits around 1 TeV are at the level of 3 % and 2 % of the Crab Nebula flux, for the two sources respectively. This implies that the Milagro emission is either extended over a much larger area than our point spread function, or it must be peaked at energies beyond 1 TeV, resulting in a photon index harder than 2.2 in the TeV band.

We report on the extreme gamma-ray activity from the FSRQ PKS 1510-089 observed by AGILE in March 2009. In the same period a radio-to-optical monitoring of the source was provided by the GASP-WEBT and REM. Moreover, several Swift ToO observations were triggered, adding important information on the source behaviour from optical/UV to hard X-rays. We paid particular attention to the calibration of the Swift/UVOT data to make it suitable to the blazars spectra. Simultaneous observations from radio to gamma rays allowed us to study in detail the correlation among the emission variability at different frequencies and to investigate the mechanisms at work. In the period 9-30 March 2009, AGILE detected an average gamma-ray flux of (311+/-21)x10^-8 ph cm^-2 s^-1 for E>100 MeV, and a peak level of (702+/-131)x10^-8 ph cm^-2 s^-1 on daily integration. The gamma-ray activity occurred during a period of increasing activity from near-IR to UV, with a flaring episode detected on 26-27 March 2009, suggesting that a single mechanism is responsible for the flux enhancement observed from near-IR to UV. By contrast, Swift/XRT observations seem to show no clear correlation of the X-ray fluxes with the optical and gamma-ray ones. However, the X-ray observations show a harder photon index (1.3-1.6) with respect to most FSRQs and a hint of harder-when-brighter behaviour, indicating the possible presence of a second emission component at soft X-ray energies. Moreover, the broad band spectrum from radio-to-UV confirmed the evidence of thermal features in the optical/UV spectrum of PKS 1510-089 also during high gamma-ray state. On the other hand, during 25-26 March 2009 a flat spectrum in the optical/UV energy band was observed, suggesting an important contribution of the synchrotron emission in this part of the spectrum during the brightest gamma-ray flare, therefore a significant shift of the synchrotron peak.

Very High Energy (VHE) gamma-ray emission from the flat spectrum radio quasar (FSRQ) PKS 1222+21 (4C 21.35, z=0.432) was detected with the MAGIC Cherenkov telescopes during a short observation (~0.5 hr) performed on 2010 June 17. The MAGIC detection coincides with high energy MeV/GeV gamma-ray activity measured by the Large Area Telescope (LAT) on board the Fermi satellite. The VHE spectrum measured by MAGIC extends from about 70 GeV up to at least 400 GeV and can be well described by a power law dN/dE \propto E^-Gamma with a photon index Gamma= 3.75+/-0.27stat +/-0.2syst. The averaged integral flux above 100 GeV is (4.56+/-0.46)x10^(-10) cm^-2 s^-1 (~1 Crab Nebula flux). The VHE flux measured by MAGIC varies significantly within the 30 min exposure implying a flux doubling time of about 10 min. The VHE and MeV/GeV spectra, corrected for the absorption by the extragalactic background light (EBL), can be described by a single power law with photon index 2.72+/-0.34 between 3 GeV and 400 GeV, and is consistent with emission belonging to a single component in the jet. The absence of a spectral cutoff constrains the gamma-ray emission region outside the Broad Line Region, which would otherwise absorb the VHE gamma-rays. Together with the detected fast variability, this challenges present emission models from jets in FSRQ. Moreover, the combined Fermi/LAT and MAGIC spectral data yield constraints on the density of the Extragalactic Background Light in the UV-optical to near-infrared range that are compatible with recent models.

Context. 3C 279, the first quasar discovered to emit VHE gamma-rays by the MAGIC telescope in 2006, was reobserved by MAGIC in January 2007 during a major optical flare and from December 2008 to April 2009 following an alert from the Fermi space telescope on an exceptionally high gamma -ray state. Aims. The January 2007 observations resulted in a detection on January 16 with significance 5.2 sigma, corresponding to a F(> 150 GeV) (3.8 \pm 0.8) \cdot 10^-11 ph cm^-2 s^-1 while the overall data sample does not show significant signal. The December 2008 - April 2009 observations did not detect the source. We study the multiwavelength behavior of the source at the epochs of MAGIC observations, collecting quasi-simultaneous data at optical and X-ray frequencies and for 2009 also gamma-ray data from Fermi. Methods. We study the light curves and spectral energy distribution of the source. The spectral energy distributions of three observing epochs (including the February 2006, which has been previously published in Albert et al. 2008a) are modeled with one-zone inverse Compton models and the emission on January 16, 2007 also with two zone model and with a lepto-hadronic model. Results. We find that the VHE gamma-ray emission detected in 2006 and 2007 challenges standard one-zone model, based on relativistic electrons in a jet scattering broad line region photons, while the other studied models fit the observed spectral energy distribution more satisfactorily.

Between 2004 and 2009 a sample of 28 X-ray selected high- and intermediate-frequency peaked blazars with a X-ray flux larger than 2 uJy at 1 keV in the redshift range from 0.018 to 0.361 was observed with the MAGIC telescope at energies above 100 GeV. Seven among them were detected and the results of these observations are discussed elsewhere. Here we concentrate on the remaining 21 blazars which were not detected during this observation campaign and present the 3 sigma (99.7 %) confidence upper limits on their flux. The individual flux upper limits lie between 1.6 % and 13.6 % of the integral flux from the Crab Nebula. Applying a stacking method to the sample of non-detections with a total of 394.1 hours exposure time, we find evidence for an excess with a cumulative significance of 4.9 standard deviations. It is not dominated by individual objects or flares, but increases linearly with the observation time as for a constant source with an integral flux level of ~1.5 % of that observed from the Crab Nebula above 150 GeV.

We report new observations of the intermediate-frequency peaked BL Lacertae object 3C 66A with the MAGIC telescopes. The data sample we use were taken in 2009 December and 2010 January, and comprises 2.3 hr of good quality data in stereoscopic mode. In this period, we find a significant signal from the direction of the blazar 3C 66A. The new MAGIC stereoscopic system is shown to play an essential role for the separation between 3C 66A and the nearby radio galaxy 3C 66B, which is at a distance of only $6^\prime$. The derived integral flux above $100\eh{GeV}$ is 8.3% of Crab Nebula flux and the energy spectrum is reproduced by a power law of photon index $3.64 \pm 0.39_{\rm stat} \pm 0.25_{\rm sys}$. Within errors, this is compatible with the one derived by VERITAS in 2009. From the spectra corrected for absorption by the extragalactic background light, we only find small differences between the four models that we applied, and constrain the redshift of the blazar to $z < 0.68$.

The very high energy (VHE; E > 100 GeV) blazar Markarian 501 has a well-studied history of extreme spectral variability and is an excellent laboratory for studying the physical processes within the jets of active galactic nuclei. However, there are few detailed multiwavelength studies of Markarian 501 during its quiescent state, due to its low luminosity. A short-term multiwavelength study of Markarian 501 was coordinated in March 2009, focusing around a multi-day observation with the Suzaku X-ray satellite and including {\gamma}-ray data from VERITAS, MAGIC, and the Fermi Gamma-ray Space Telescope with the goal of providing a well-sampled multiwavelength baseline measurement of Markarian 501 in the quiescent state. The results of these quiescent-state observations are compared to the historically extreme outburst of April 16, 1997, with the goal of examining variability of the spectral energy distribution between the two states. The derived broadband spectral energy distribution shows the characteristic double-peaked profile. We find that the X-ray peak shifts by over two orders of magnitude in photon energy between the two flux states while the VHE peak varies little. The limited shift in the VHE peak can be explained by the transition to the Klein-Nishina regime. Synchrotron self-Compton models are matched to the data and the implied Klein-Nishina effects are explored.

We report on the detection with the MAGIC telescopes of very high energy gamma-rays from IC 310, a head-tail radio galaxy in the Perseus galaxy cluster, observed during the interval November 2008 to February 2010. The Fermi satellite has also detected this galaxy. The source is detected by MAGIC at a high statistical significance of 7.6sigma in 20.6 hr of stereo data. The observed spectral energy distribution is flat with a differential spectral index of -2.00 \pm 0.14. The mean flux above 300 GeV, between October 2009 and February 2010, (3.1 \pm 0.5)x10^{-12} cm^{-2} s^{-1}, corresponds to (2.5 \pm 0.4)% of Crab Nebula units. Only an upper limit, of 1.9% of Crab Nebula units above 300 GeV, was obtained with the 2008 data. This, together with strong hints (>3sigma) of flares in the middle of October and November 2009, implies that the emission is variable. The MAGIC results favour a scenario with the very high energy emission originating from the inner jet close to the central engine. More complicated models than a simple one-zone SSC scenario, e.g. multi-zone SSC, external Compton or hadronic, may be required to explain the very flat spectrum and its extension over more than three orders of magnitude in energy.

Cygnus X-3 is a microquasar consisting of an accreting compact object orbiting around a Wolf-Rayet star. It has been detected at radio frequencies and up to high-energy gamma rays (above 100 MeV). However, many models also predict a very high energy (VHE) emission (above hundreds of GeV) when the source displays relativistic persistent jets or transient ejections. Therefore, detecting such emission would improve the understanding of the jet physics. The imaging atmospheric Cherenkov telescope MAGIC observed Cygnus X-3 for about 70 hours between 2006 March and 2009 August in different X-ray/radio spectral states and also during a period of enhanced gamma-ray emission. MAGIC found no evidence for a VHE signal from the direction of the microquasar. An upper limit to the integral flux for energies higher than 250 GeV has been set to 2.2 x 10-12 photons cm-2 s-1 (95% confidence level). This is the best limit so far to the VHE emission from this source. The non-detection of a VHE signal during the period of activity in the high-energy band sheds light on the location of the possible VHE radiation favoring the emission from the innermost region of the jets, where absorption is significant. The current and future generations of Cherenkov telescopes may detect a signal under precise spectral conditions.

The Major Atmospheric Gamma Imaging Cerenkov (MAGIC) telescope participated in three multiwavelength (MWL) campaigns, observing the blazar Markarian (Mkn) 421 during the nights of 2006 April 28, 29, and 2006 June 14. We analyzed the corresponding MAGIC very-high energy observations during 9 nights from 2006 April 22 to 30 and on 2006 June 14. We inferred light curves with sub-day resolution and night-by-night energy spectra. A strong gamma-ray signal was detected from Mkn 421 on all observation nights. The flux (E > 250 GeV) varied on night-by-night basis between (0.92+-0.11)10^-10 cm^-2 s^-1 (0.57 Crab units) and (3.21+-0.15)10^-10 cm^-2 s^-1 (2.0 Crab units) in 2006 April. There is a clear indication for intra-night variability with a doubling time of 36+-10(stat) minutes on the night of 2006 April 29, establishing once more rapid flux variability for this object. For all individual nights gamma-ray spectra could be inferred, with power-law indices ranging from 1.66 to 2.47. We did not find statistically significant correlations between the spectral index and the flux state for individual nights. During the June 2006 campaign, a flux substantially lower than the one measured by the Whipple 10-m telescope four days later was found. Using a log-parabolic power law fit we deduced for some data sets the location of the spectral peak in the very-high energy regime. Our results confirm the indications of rising peak energy with increasing flux, as expected in leptonic acceleration models.

Context: Part of the very high energy $\gamma$-ray radiation coming from extragalactic sources is absorbed through the pair production process on the extragalactic background light photons. Extragalactic magnetic fields alter the trajectories of these cascade pairs and, in turn, convert cosmic background photons to gamma-ray energies by inverse Compton scattering. These secondary photons can form an extended halo around bright VHE sources. Aims: We searched for an extended emission around the bright blazars Mrk 421 and Mrk 501 using the MAGIC telescope data. Methods: If extended emission is present, the angular distribution of reconstructed gamma-ray arrival directions around the source is broader than for a point-like source. In the analysis of a few tens of hours of observational data taken from Mrk 421 and Mrk 501 we used a newly developed method that provides better angular resolution. This method is based on the usage of multidimensional decision trees. Comparing the measured shapes of angular distributions with those expected from a point-like source one can detect or constrain possible extended emission around the source. We also studied the influence of different types of systematic errors on the shape of the distribution of reconstructed gamma-ray arrival directions for a point source. Results: We present upper limits for an extended emission calculated for both sources for various source extensions and emission profiles. We discuss possible constraints on the extragalactic magnetic fields strength. We obtain upper limits on the extended emission around the Mrk~421 (Mrk~501) on the level of < 5% (< 4%) of the Crab Nebula flux above the energy threshold of 300 GeV. Using these results we discuss possible constraints on the extragalactic magnetic fields strength around a few times $10^{-15}$~G.

The blazar PG 1553+113 is a well known TeV gamma-ray emitter. In this paper, we determine its spectral energy distribution using simultaneous multi-frequency data in order to study its emission processes. An extensive campaign was carried out between March and April 2008, where optical, X-ray, high-energy (HE) gamma-ray, and very-high-energy (VHE) gamma-ray data were obtained with the KVA, Abastumani, REM, RossiXTE/ASM, AGILE and MAGIC telescopes, respectively. This is the first simultaneous broad-band (i.e., HE+VHE) gamma-ray observation, though AGILE did not detect the source. We combine data to derive source's spectral energy distribution and interpret its double peaked shape within the framework of a synchrotron self compton model

Context: Gamma-ray bursts are cosmological sources emitting radiation from the gamma-rays to the radio band. Substantial observational efforts have been devoted to the study of gamma-ray bursts during the prompt phase, i.e. the initial burst of high-energy radiation, and during the long-lasting afterglows. In spite of many successes in interpreting these phenomena, there are still several open key questions about the fundamental emission processes, their energetics and the environment. Aim: Independently of specific gamma-ray burst theoretical recipes, spectra in the GeV/TeV range are predicted to be remarkably simple, being satisfactorily modeled with power-laws, and therefore offer a very valuable tool to probe the extragalactic background light distribution. Furthermore, the simple detection of a component at very-high energies, i.e. at $\sim 100$\,GeV, would solve the ambiguity about the importance of various possible emission processes, which provide barely distinguishable scenarios at lower energies. Methods: We used the results of the MAGIC telescope observation of the moderate resdhift ($z\sim0.76$) \object{GRB\,080430} at energies above about 80\,GeV, to evaluate the perspective for late-afterglow observations with ground based GeV/TeV telescopes. Results: We obtained an upper limit of $F_{\rm 95%\,CL} = 5.5 \times 10^{-11}$\,erg\,cm$^{-2}$\,s$^{-1}$ for the very-high energy emission of \object{GRB\,080430}, which cannot set further constraints on the theoretical scenarios proposed for this object also due to the difficulties in modeling the low-energy afterglow. Nonetheless, our observations show that Cherenkov telescopes have already reached the required sensitivity to detect the GeV/TeV emission of GRBs at moderate redshift ($z \lesssim 0.8$), provided the observations are carried out at early times, close to the onset of their afterglow phase.

The Perseus galaxy cluster was observed by the MAGIC Cherenkov telescope for a total effective time of 24.4 hr during 2008 November and December. The resulting upper limits on the gamma-ray emission above 100 GeV are in the range of 4.6 to 7.5 x 10^{-12} cm^{-2} s^{-1} for spectral indices from -1.5 to -2.5, thereby constraining the emission produced by cosmic rays, dark matter annihilations, and the central radio galaxy NGC1275. Results are compatible with cosmological cluster simulations for the cosmic-ray-induced gamma-ray emission, constraining the average cosmic ray-to-thermal pressure to <4% for the cluster core region (<8% for the entire cluster). Using simplified assumptions adopted in earlier work (a power-law spectrum with an index of -2.1, constant cosmic ray-to-thermal pressure for the peripheral cluster regions while accounting for the adiabatic contraction during the cooling flow formation), we would limit the ratio of cosmic ray-to-thermal energy to E_CR/E_th<3%. The upper limit also translates into a level of gamma-ray emission from possible annihilations of the cluster dark matter (the dominant mass component) that is consistent with boost factors of ~10^4 for the typically expected dark matter annihilation-induced emission. Finally, the upper limits obtained for the gamma-ray emission of the central radio galaxy NGC1275 are consistent with the recent detection by the Fermi-LAT satellite. Due to the extremely large Doppler factors required for the jet, a one-zone synchrotron self-Compton model is implausible in this case. We reproduce the observed spectral energy density by using the structured jet (spine-layer) model which has previously been adopted to explain the high-energy emission of radio galaxies.

(Abridged) OJ287 is a BL Lac object that has shown double-peaked bursts at regular intervals of ~12 yr during the last ~40 yr. We analyse optical photopolarimetric monitoring data from 2005-2009, during which the latest double-peaked outburst occurred. The aim of this study is twofold: firstly, we aim to analyse variability patterns and statistical properties of the optical polarization light-curve. We find a strong preferred position angle in optical polarization. The preferred position angle can be explained by separating the jet emission into two components: an optical polarization core and chaotic jet emission. The optical polarization core is stable on time scales of years and can be explained as emission from an underlying quiescent jet component. The chaotic jet emission sometimes exhibits a circular movement in the Stokes plane. We interpret these events as a shock front moving forwards and backwards in the jet, swiping through a helical magnetic field. Secondly, we use our data to assess different binary black hole models proposed to explain the regularly appearing double-peaked bursts in OJ287. We compose a list of requirements a model has to fulfil. The list includes not only characteristics of the light-curve but also other properties of OJ287, such as the black hole mass and restrictions on accretion flow properties. We rate all existing models using this list and conclude that none of the models is able to explain all observations. We discuss possible new explanations and propose a new approach to understanding OJ287. We suggest that both the double-peaked bursts and the evolution of the optical polarization position angle could be explained as a sign of resonant accretion of magnetic field lines, a 'magnetic breathing' of the disc.

The discovery of very high energy (VHE) gamma-ray emitting X-ray binaries has triggered an intense effort to better understand the particle acceleration, absorption, and emission mechanisms in compact binary systems, which provide variable conditions along eccentric orbits. Despite this, the nature of some of these systems, and of the accelerated particles producing the VHE emission, is unclear. To answer some of these open questions, we conducted a multiwavelength campaign of the VHE gamma-ray emitting X-ray binary LS I +61 303 including the MAGIC telescope, XMM-Newton, and Swift during 60% of an orbit in 2007 September. We detect a simultaneous outburst at X-ray and VHE bands, with the peak at phase 0.62 and a similar shape at both wavelengths. A linear fit to the simultaneous X-ray/VHE pairs obtained during the outburst yields a correlation coefficient of r=0.97, while a linear fit to all simultaneous pairs provides r=0.81. Since a variable absorption of the VHE emission towards the observer is not expected for the data reported here, the correlation found indicates a simultaneity in the emission processes. Assuming that they are dominated by a single particle population, either hadronic or leptonic, the X-ray/VHE flux ratio favors leptonic models. This fact, together with the detected photon indices, suggests that in LS I +61 303 the X-rays are the result of synchrotron radiation of the same electrons that produce VHE emission as a result of inverse Compton scattering of stellar photons.

We report results from an intensive multiwavelength campaign on the intermediate-frequency-peaked BL Lacertae object W Com (z=0.102) during a strong outburst of very high energy gamma-ray emission in June 2008. The very high energy gamma-ray signal was detected by VERITAS on 2008 June 7-8 with a flux F(>200 GeV) = (5.7+-0.6)x10^-11 cm-2s-1, about three times brighter than during the discovery of gamma-ray emission from W Com by VERITAS in 2008 March. The initial detection of this flare by VERITAS at energies above 200 GeV was followed by observations in high energy gamma-rays (AGILE, E>100 MeV), and X-rays (Swift and XMM-Newton), and at UV, and ground-based optical and radio monitoring through the GASP-WEBT consortium and other observatories. Here we describe the multiwavelength data and derive the spectral energy distribution (SED) of the source from contemporaneous data taken throughout the flare.

We report the detection by the AGILE satellite of a rapid gamma-ray flare from the powerful gamma-ray quasar PKS 1510-089, during a pointing centered on the Galactic Center region from 1 March to 30 March 2008. This source has been continuosly monitored in the radio-to-optical bands by the GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT). Moreover, the gamma-ray flaring episode triggered three ToO observations by the Swift satellite in three consecutive days, starting from 20 March 2008. In the period 1-16 March 2008, AGILE detected gamma-ray emission from PKS 1510-089 at a significance level of 6.2-sigma with an average flux over the entire period of (84 +/- 17) x 10^{-8} photons cm^{-2} s^{-1} for photon energies above 100 MeV. After a predefined satellite re-pointing, between 17 and 21 March 2008, AGILE detected the source at a significance level of 7.3-sigma, with an average flux (E > 100 MeV) of (134 +/- 29) x 10^{-8} photons cm^{-2} s^{-1} and a peak level of (281 +/- 68) x 10^{-8} photons cm^{-2} s^{-1} with daily integration. During the observing period January-April 2008, the source also showed an intense and variable optical activity, with several flaring episodes and a significant increase of the flux was observed at millimetric frequencies. Moreover, in the X-ray band the Swift/XRT observations seem to show an harder-when-brighter behaviour of the source spectrum. The spectral energy distribution of mid-March 2008 is modelled with a homogeneous one-zone synchrotron self Compton emission plus contributions from inverse Compton scattering of external photons from both the accretion disc and the broad line region. Indeed, some features in the optical-UV spectrum seem to indicate the presence of Seyfert-like components, such as the little blue bump and the big blue bump.

Index entry listing the contributed papers of the MAGIC collaboration to the 31th International Cosmic Ray Conference (ICRC 2009), July 7-15 2009, \L\'od\'z, Poland. The individual papers are sorted by subject: Overview and Highlight Papers; MAGIC-II Status and Components; Software and Analysis Techniques; Technical Developments; Scientific Results. This HTML document includes clickable links to the papers that exist on the astro-ph arXiv. We hope that this will make it easy to access the MAGIC contributions in a systematic way.

We report the gamma-ray activity from the intermediate BL Lac S5 0716+714 during observations acquired by the AGILE satellite in September and October 2007. These detections of activity were contemporaneous with a period of intense optical activity, which was monitored by GASP-WEBT. This simultaneous optical and gamma-ray coverage allows us to study in detail the light curves, time lags, gamma-ray photon spectrum, and Spectral Energy Distributions (SEDs) during different states of activity. AGILE observed the source with its two co-aligned imagers, the Gamma-Ray Imaging Detector (GRID) and the hard X-ray imager (Super-AGILE), which are sensitive to the 30 MeV-50 GeV and 18-60 keV energy ranges, respectively. Observations were completed in two different periods, the first between 2007 September 4-23, and the second between 2007 October 24-November 1. Over the period 2007 September 7-12, AGILE detected gamma-ray emission from the source at a significance level of 9.6-sigma with an average flux (E>100 MeV) of (97 +/- 15) x 10^{-8} photons cm^{-2} s^{-1}, which increased by a factor of at least four within three days. No emission was detected by Super-AGILE for the energy range 18-60 keV to a 3-sigma upper limit of 10 mCrab in 335 ksec. In October 2007, AGILE repointed toward S5 0716+714 following an intense optical flare, measuring an average flux of (47 +/- 11) x 10^{-8} photons cm^{-2} s^{-1} at a significance level of 6.0-sigma. The gamma-ray flux of S5 0716+714 detected by AGILE is the highest ever detected for this blazar and one of the most intense gamma-ray fluxes detected from a BL Lac object. The SED of mid-September appears to be consistent with the synchrotron self-Compton (SSC) emission model, but only by including two SSC components of different variabilities.

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

We have derived R-band host galaxy fluxes of 20 TeV candidate BL Lacertae objects as a function of aperture radius and FWHM. The results are given as correction tables, that list the fluxes (in mJy) of all ``contaminating'' sources (host galaxy + significant nearby objects) as a function of aperture radius and FWHM. We found that the derived fluxes depend strongly on aperture radius, but the FWHM has only a minor effect (a few percent). We also discuss the implications of our findings to optical monitoring programs and potential sources of error in our derived fluxes. During this work we have also constructed new calibration star sequences for 9 objects and present the finding charts and calibrated magnitudes.

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

We report first results of a multifrequency campaign from radio to hard X-ray energies of the prominent gamma-ray blazar 3C 279 during the first year of the INTEGRAL mission. The variable blazar was found at a low activity level, but was detected by all participating instruments. Subsequently a multifrequency spectrum could be compiled. The individual measurements as well as the compiled multifrequency spectrum are presented. In addition, this 2003 broadband spectrum is compared to one measured in 1999 during a high activity period of 3C 279.