We report on Chandra observations of the TeV emitting High Mass X-ray Binary LS 5039, for a total exposure of ~70ks, using the ACIS-S camera in Continuos Clocking mode to search for a possible X-ray pulsar in this system. We did not find any periodic or quasi-periodic signal in the 0.3-0.4 and 0.75-0.9 orbital phases, and in a frequency range of 0.005-175 Hz. We derived an average pulsed fraction 3sigma upper limit for the presence of a periodic signal of ~15% (depending on the frequency and the energy band), the deepest limit ever reached for this object. If the X-ray emission of LS 5039 is due (at least in part) to a rotational powered pulsar, the latter is either spinning faster than ~5.6 ms, or having a beam pointing away from our line of sight, or contributing to ~15% of the total X-ray emission of the system in the orbital phases we observed.

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.

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.

We report on a 95ks Chandra observation of the TeV emitting High Mass X-ray Binary LSI +61 303, using the ACIS-S camera in Continuos Clocking mode to search for a possible X-ray pulsar in this system. The observation was performed while the compact object was passing from phase 0.94 to 0.98 in its orbit around the Be companion star (hence close to the apastron passage). We did not find any periodic or quasi-periodic signal (at this orbital phase) in a frequency range of 0.005-175 Hz. We derived an average pulsed fraction 3 sigma upper limit for the presence of a periodic signal of ~10% (although this limit is strongly dependent on the frequency and the energy band), the deepest limit ever reached for this object. Furthermore, the source appears highly variable in flux and spectrum even in this very small orbital phase range, in particular we detect two flares, lasting thousands of seconds, with a very hard X-ray spectrum with respect to the average source spectral distribution. The X-ray pulsed fraction limits we derived are lower than the pulsed fraction of any isolated rotational-powered pulsar, in particular having a TeV counterpart. In this scenario most of the X-ray emission of LSI +61 303 should necessarily come from the interwind or inner-pulsar wind zone shock rather than from the magnetosphere of the putative pulsar. Furthermore, we did not find evidence for the previously suggested extended X-ray emission (abridged).

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.

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.

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.

LS I +61 303 has been recently detected as a periodic gamma-ray source by the Major Atmospheric Imaging Cerenkov (MAGIC) telescope. A distinctive orbital correlation of the $\gamma$-ray emission was found. This work shows that the range of uncertainties yet at hand in the orbital elements of the binary system LSI as well as in the possible assumptions on the stellar wind of the optical companion play a non-negligible role in the computation of opacities to high energy processes leading to $\gamma$-ray predictions. The geometry influence on the propagation and escape of $\gamma$-ray photons is explored. With this study at hand, we analyse the results of a pulsar wind zone model for the production of gamma-rays and compare it with recent MAGIC observations.

Several gamma-ray binaries have been recently detected by the High-Energy Stereoscopy Array (H.E.S.S.) and the Major Atmospheric Imaging Cerenkov (MAGIC) telescope. In two cases, their nature is unknown, since a final observational feature for a black hole or a pulsar compact object companion is still missing. One such system is LS 5039. Here we present results from a model (it includes a detailed account of the system geometry, the angular dependence of processes such as Klein-Nishina inverse Compton and gamma-gamma absorption, and a Monte Carlo simulation of cascading) of the high energy phenomenology of LS 5039 in which it is assumed that the companion object is a pulsar rotating around an O6.5V star in the 3.9 days orbit. We show that the H.E.S.S. phenomenology at all scales (spectra along the orbit in both broad and short phase-bins and lightcurve) is described within this model. We focus on presenting predictions for photons with lower energies (for E>1 GeV), subject to test in the forthcoming months with the GLAST satellite, and we also present predictions for future observations with high-energy arrays, such as H.E.S.S. II. Both set of predictions go beyond the description of the current data, and could provide a high-energy determination of the system's composition.

Several $\gamma$-ray binaries have been recently detected by the High-Energy Stereoscopy Array (H.E.S.S.) and the Major Atmospheric Imaging Cerenkov (MAGIC) telescope. In at least two cases, their nature is unknown. In this paper we aim to provide the details of a theoretical model of close $\gamma$-ray binaries containing a young energetic pulsar as compact object, earlier presented in recent Letters. This model includes a detailed account of the system geometry, the angular dependence of processes such as Klein-Nishina inverse Compton and $\gamma\gamma$ absorption in the anisotropic radiation field of the massive star, and a Monte Carlo simulation of leptonic cascading. We present and derive the used formulae and give all details about their numerical implementation, particularly, on the computation of cascades. In this model, emphasis is put in the processes occurring in the pulsar wind zone of the binary, since, as we show, opacities in this region can be already important for close systems. We provide a detailed study on all relevant opacities and geometrical dependencies along the orbit of binaries, exemplifying with the case of LS 5039. This is used to understand the formation of the very high-energy lightcurve and phase dependent spectrum. For the particular case of LS 5039, we uncover an interesting behavior of the magnitude representing the shock position in the direction to the observer along the orbit, and analyze its impact in the predictions. We show that in the case of LS 5039, the H.E.S.S. phenomenology is matched by the presented model, and explore the reasons why this happens while discussing future ways of testing the model.

At least one massive binary system containing an energetic pulsar, PSR B1259-63/SS2883, has been recently detected in the TeV gamma-rays by the HESS telescopes. These gamma-rays are likely produced by particles accelerated in the vicinity of the pulsar and/or at the pulsar wind shock, in comptonization of soft radiation from the massive star. However, the process of gamma-ray production in such systems can be quite complicated due to the anisotropy of the radiation field, complex structure of the pulsar wind termination shock and possible absorption of produced gamma-rays which might initiate leptonic cascades. In this paper we consider in detail all these effects. We calculate the gamma-ray light curves and spectra for different geometries of the binary system PSR B1259-63/SS2883 and compare them with the TeV gamma-ray observations. We conclude that the leptonic IC model, which takes into account the complex structure of the pulsar wind shock due to the aspherical wind of the massive star, can explain the details of the observed gamma-ray light curve.

LS 5039 is one of a handful of X-ray binaries that have been recently detected at high-energy $\gamma$-rays, in this case, by the High-Energy Stereoscopy Array (H.E.S.S.). The nature of this system is unknown: both a black hole and a pulsar have been invoked as possible compact object companions. Here we work with a model of the high energy phenomenology of the system in which it is assumed that the companion object is a pulsar rotating around an O6.5V star in the $\sim 3.9$ days orbit. The model assumes two different sets of power-law spectral parameters of the interacting primary leptons corresponding to the two orbital phase intervals defined by H.E.S.S. as having different gamma-ray spectra and very-high-energy (VHE) cutoffs. We show the H.E.S.S. phenomenology is completely explained by this model. We present predictions for photons with lower energies (for $E>1 $ GeV), subject to test in the forthcoming months with the GLAST satellite. We find that GLAST is able to judge on this model within one year.

Under the assumption that LS 5039 is a system composed by a pulsar rotating around an O6.5V star in a $\sim 3.9$ day orbit, we present the results of a theoretical modeling of the high energy phenomenology observed by the High Energy Stereoscopy Array (H.E.S.S.). This model (including detailed account of the system geometry, Klein-Nishina inverse Compton, $\gamma$-$\gamma$ absorption, and cascading) is able to describe well the rich observed phenomenology found in the system at all timescales, both flux and spectrum-wise.