A key characteristic of many active galactic nuclei (AGN) is their variability, but its origin is poorly understood, especially in the radio domain. Williams et al. (2017) reported a ~50 per cent increase in peak flux density of the AGN in the Seyfert galaxy NGC 4151 at 1.5 GHz with the e-MERLIN array. We present new high resolution e-MERLIN observations at 5 GHz and compare these to archival MERLIN observations to investigate the reported variability. Our new observations allow us to probe the nuclear region at a factor three times higher-resolution than the previous e-MERLIN study. We separate the core component, C4, into three separate components: C4W, C4E and X. The AGN is thought to reside in component C4W, but this component has remained constant between epochs within uncertainties. However, we find that the Eastern-most component, C4E, has increased in peak flux density from 19.35$\pm$1.10 to 37.09$\pm$1.86 mJy/beam, representing a 8.2 sigma increase on the MERLIN observations. We attribute this peak flux density increase to continued interaction between the jet and the emission line region (ELR), observed for the first time in a low-luminosity AGN such as NGC 4151. We identify discrete resolved components at 5 GHz along the jet axis, which we interpret as areas of jet-ELR interaction.

We present the first LOFAR observations of the radio jet in the quasar 4C+19.44 (a.k.a. PKS 1354+19) obtained with the long baselines. The achieved resolution is very well matched to that of archival Jansky Very Large Array (JVLA) observations at higher radio frequencies as well as the archival X-ray images obtained with {\it Chandra}. We found that, for several knots along the jet, the radio flux densities measured at hundreds of MHz lie well below the values estimated by extrapolating the GHz spectra. This clearly indicates the presence of spectral curvature. Radio spectral curvature has been already observed in different source classes and/or extended radio structures and it has been often interpreted as due to intrinsic processes, as a curved particle energy distribution, rather than absorption mechanisms ({ Razin-Tsytovich} effect, free-free or synchrotron self absorption to name a few). Here we discuss our results according to the scenario where particles undergo stochastic acceleration mechanisms also in quasar jet knots.

The pulsar/massive star binary system PSR B1259-63 / LS 2883 is one of the best-studied gamma-ray binaries, a class of systems whose bright gamma-ray flaring can provide important insights into high-energy physics. Using the Australian Long Baseline Array we have conducted very long baseline interferometric observations of PSR B1259-63 over 4.4 years, fully sampling the 3.4-year orbital period. From our measured parallax of $0.38\pm0.05$ mas we use a Bayesian approach to infer a distance of $2.6^{+0.4}_{-0.3}$ kpc. We find that the binary orbit is viewed at an angle of $154\pm3$ degrees to the line of sight, implying that the pulsar moves clockwise around its orbit as viewed on the sky. Taking our findings together with previous results from pulsar timing observations, all seven orbital elements for the system are now fully determined. We use our measurement of the inclination angle to constrain the mass of the stellar companion to lie in the range 15-31$M_{\odot}$. Our measured distance and proper motion are consistent with the system having originated in the Cen OB1 association and receiving a modest natal kick, causing it to have moved $\sim$8 pc from its birthplace over the past $\sim3\times10^5$ years. The orientation of the orbit on the plane of the sky matches the direction of motion of the X-ray synchrotron-emitting knot observed by the Chandra X-ray Observatory to be moving away from the system.

We report on the first sub-arcsecond (0.44 $\times$ 0.41 arcsec$\rm ^2$) angular resolution image at 150 MHz of the A-nucleus in the Luminous Infrared Galaxy Arp$\,$299, from International Low Frequency Array (LOFAR) Telescope observations. The most remarkable finding is that of an intriguing two-sided, filamentary structure emanating from A-nucleus, which we interpret as an outflow that extends up to at least 14 arcseconds from the A-nucleus in the N-S direction ($\approx$ 5 kpc deprojected size) and accounts for almost 40% of the extended emission of the entire galaxy system. We also discuss HST/NICMOS [FeII] 1.64 $\rm \mu m$ and H$\rm_2$ 2.12 $\rm \mu m$ images of Arp$\,$299-A, which show similar features to those unveiled by our 150 MHz LOFAR observations, thus giving string morphological support for the outflow scenario. Finally, we discuss unpublished NaI D spectra that confirm the outflow nature of this structure. From energetic arguments, we rule out the low-luminosity active galactic nucleus in Arp$\,$299-A as a driver for the outflow. On the contrary, the powerful, compact starburst in the central regions of Arp$\,$299-A provides plenty of mechanical energy to sustain an outflow, and we conclude that the intense supernova (SN) activity in the nuclear region of Arp299-A is driving the observed outflow. We estimate that the starburst wind can support a mass-outflow rate in the range (11-63) $\rm M_{\odot} yr^{-1}$ at speeds of up to (370 - 890) $\rm km \, s^{-1}$, and is relatively young, with an estimated kinematic age of (3 - 7) Myr. Those results open an avenue to the use of low-frequency (150 MHz), sub-arcsecond imaging with LOFAR to detect outflows in the central regions of local luminous infrared galaxies.

MWC 656 is the first known Be/black hole (BH) binary system. Be/BH binaries are important in the context of binary system evolution and sources of detectable gravitational waves because they are possible precursors of coalescing neutron star/BH binaries. X-ray observations conducted in 2013 revealed that MWC 656 is a quiescent high-mass X-ray binary (HMXB), opening the possibility to explore X-ray/radio correlations and the accretion/ejection coupling down to low luminosities for BH HMXBs. Here we report on a deep joint Chandra/VLA observation of MWC 656 (and contemporaneous optical data) conducted in 2015 July that has allowed us to unambiguously identify the X-ray counterpart of the source. The X-ray spectrum can be fitted with a power law with $\Gamma\sim2$, providing a flux of $\simeq4\times10^{-15}$ erg cm$^{-2}$ s$^{-1}$ in the 0.5-8 keV energy range and a luminosity of $L_{\rm X}\simeq3\times10^{30}$ erg s$^{-1}$ at a 2.6 kpc distance. For a 5 M$_\odot$ BH this translates into $\simeq5\times10^{-9}$ $L_{\rm Edd}$. These results imply that MWC 656 is about 7 times fainter in X-rays than it was two years before and reaches the faintest X-ray luminosities ever detected in stellar-mass BHs. The radio data provide a detection with a peak flux density of $3.5\pm1.1$ $\mu$Jy beam$^{-1}$. The obtained X-ray/radio luminosities for this quiescent BH HMXB are fully compatible with those of the X-ray/radio correlations derived from quiescent BH low-mass X-ray binaries. These results show that the accretion/ejection coupling in stellar-mass BHs is independent of the nature of the donor star.

(abridged). We outline LBCS (the LOFAR Long-Baseline Calibrator Survey), whose aim is to identify sources suitable for calibrating the highest-resolution observations made with the International LOFAR Telescope, which include baselines >1000 km. Suitable sources must contain significant correlated flux density (50-100mJy) at frequencies around 110--190~MHz on scales of a few hundred mas. At least for the 200--300-km international baselines, we find around 1 suitable calibrator source per square degree over a large part of the northern sky, in agreement with previous work. This should allow a randomly selected target to be successfully phase calibrated on the international baselines in over 50% of cases. Products of the survey include calibrator source lists and fringe-rate and delay maps of wide areas -- typically a few degrees -- around each source. The density of sources with significant correlated flux declines noticeably with baseline length over the range 200--600~km, with good calibrators on the longest baselines appearing only at the rate of 0.5 per square degree. Coherence times decrease from 1--3 minutes on 200-km baselines to about 1 minute on 600-km baselines, suggesting that ionospheric phase variations contain components with scales of a few hundred kilometres. The longest median coherence time, at just over 3 minutes, is seen on the DE609 baseline, which at 227km is close to being the shortest. We see median coherence times of between 80 and 110 seconds on the four longest baselines (580--600~km), and about 2 minutes for the other baselines. The success of phase transfer from calibrator to target is shown to be influenced by distance, in a manner that suggests a coherence patch at 150-MHz of the order of 1 degree.

We analyse new observations with the International Low Frequency Array (LOFAR) telescope, and archival data from the Multi-Element Radio Linked Interferometer Network (MERLIN) and the Karl G. Jansky Very Large Array (VLA). We model the spatially resolved radio spectrum of Arp 220 from 150 MHz to 33 GHz. We present an image of Arp 220 at 150 MHz with resolution $0.65''\times0.35''$, sensitivity 0.15 mJy beam$^{-1}$, and integrated flux density $394\pm59$ mJy. More than 80% of the detected flux comes from extended ($6''\approx$2.2 kpc) steep spectrum ($\alpha=-0.7$) emission, likely from star formation in the molecular disk surrounding the two nuclei. We find elongated features extending $0.3''$ (110 pc) and $0.9''$ (330 pc) from the eastern and western nucleus respectively, which we interpret as evidence for outflows. The extent of radio emission requires acceleration of cosmic rays far outside the nuclei. We find that a simple three component model can explain most of the observed radio spectrum of the galaxy. When accounting for absorption at 1.4 GHz, Arp 220 follows the FIR/radio correlation with $q=2.36$, and we estimate a star formation rate of 220 M$_\odot\text{yr}^{-1}$. We derive thermal fractions at 1 GHz of less than 1% for the nuclei, which indicates that a major part of the UV-photons are absorbed by dust. International LOFAR observations shows great promise to detect steep spectrum outflows and probe regions of thermal absorption. However, in LIRGs the emission detected at 150 MHz does not necessarily come from the main regions of star formation. This implies that high spatial resolution is crucial for accurate estimates of star formation rates for such galaxies at 150 MHz.

The correlation between radio spectral index and redshift has been exploited to discover high redshift radio galaxies, but its underlying cause is unclear. It is crucial to characterise the particle acceleration and loss mechanisms in high redshift radio galaxies to understand why their radio spectral indices are steeper than their local counterparts. Low frequency information on scales of $\sim$1 arcsec are necessary to determine the internal spectral index variation. In this paper we present the first spatially resolved studies at frequencies below 100 MHz of the $z = 2.4$ radio galaxy 4C 43.15 which was selected based on its ultra-steep spectral index ($\alpha < -1$; $S_{\nu} \sim \nu^{\alpha}$ ) between 365 MHz and 1.4 GHz. Using the International Low Frequency Array (LOFAR) Low Band Antenna we achieve sub-arcsecond imaging resolution at 55 MHz with VLBI techniques. Our study reveals low-frequency radio emission extended along the jet axis, which connects the two lobes. The integrated spectral index for frequencies $<$ 500 MHz is -0.83. The lobes have integrated spectral indices of -1.31$\pm$0.03 and -1.75$\pm$0.01 for frequencies $\geq$1.4 GHz, implying a break frequency between 500 MHz and 1.4 GHz. These spectral properties are similar to those of local radio galaxies. We conclude that the initially measured ultra-steep spectral index is due to a combination of the steepening spectrum at high frequencies with a break at intermediate frequencies.

The transitional millisecond pulsar binary system PSR J1023+0038 re-entered an accreting state in 2013 June, in which it bears many similarities to low-mass X-ray binaries (LMXBs) in quiescence or near-quiescence. At a distance of just 1.37 kpc, PSR J1023+0038 offers an unsurpassed ability to study low-level accretion onto a highly-magnetized compact object. We have monitored PSR J1023+0038 intensively using radio imaging with the Karl G. Jansky Very Large Array, the European VLBI Network and LOFAR, seeing rapidly variable, flat spectrum emission that persists over a period of 6 months. The flat spectrum and variability are indicative of synchrotron emission originating in an outflow from the system, most likely in the form of a compact, partially self-absorbed jet, as is seen in LMXBs at higher accretion rates. The radio brightness, however, greatly exceeds extrapolations made from observations of more vigorously accreting neutron star LMXB systems. We postulate that PSR J1023+0038 is undergoing radiatively inefficient 'propeller-mode' accretion, with the jet carrying away a dominant fraction of the liberated accretion luminosity. We confirm that the enhanced gamma-ray emission seen in PSR J1023+0038 since it re-entered an accreting state has been maintained; the increased gamma-ray emission in this state can also potentially be associated with propeller-mode accretion. Similar accretion modes can be invoked to explain the radio and X-ray properties of the other two known transitional millisecond pulsar systems XSS J12270-4859 and PSR J1824-2452I (M28I), suggesting that radiatively inefficient accretion may be an ubiquitous phenomenon amongst (at least one class of) neutron star binaries at low accretion rates.

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.

The nuclear starburst in the nearby galaxy M82 provides an excellent laboratory for understanding the physics of star formation. This galaxy has been extensively observed in the past, revealing tens of radio-bright compact objects embedded in a diffuse free-free absorbing medium. Our understanding of the structure and physics of this medium in M82 can be greatly improved by high-resolution images at low frequencies where the effects of free-free absorption are most prominent. The aims of this study are, firstly, to demonstrate imaging using international baselines of the Low Frequency Array (LOFAR), and secondly, to constrain low-frequency spectra of compact and diffuse emission in the central starburst region of M82 via high-resolution radio imaging at low frequencies. The international LOFAR telescope was used to observe M82 at 110-126MHz and 146-162MHz. Images were obtained using standard techniques from very long baseline interferometry. images were obtained at each frequency range: one only using international baselines, and one only using the longest Dutch (remote) baselines. The 154MHz image obtained using international baselines is a new imaging record in terms of combined image resolution (0.3$"$) and sensitivity ($\sigma$=0.15mJy/beam) at low frequencies ($<327$MHz). We detected 16 objects at 154MHz, six of these also at 118MHz. Four weaker but resolved features are also found: a linear (50pc) filament and three other resolved objects, of which two show a clear shell structure. We do not detect any emission from either supernova 2008iz or from the radio transient source 43.78+59.3. The images obtained using remote baselines show diffuse emission, associated with the outflow in M82, with reduced brightness in the region of the edge-on star-forming disk.

Aims. An efficient means of locating calibrator sources for International LOFAR is developed and used to determine the average density of usable calibrator sources on the sky for subarcsecond observations at 140 MHz. Methods. We used the multi-beaming capability of LOFAR to conduct a fast and computationally inexpensive survey with the full International LOFAR array. Sources were pre-selected on the basis of 325 MHz arcminute-scale flux density using existing catalogues. By observing 30 different sources in each of the 12 sets of pointings per hour, we were able to inspect 630 sources in two hours to determine if they possess a sufficiently bright compact component to be usable as LOFAR delay calibrators. Results. Over 40% of the observed sources are detected on multiple baselines between international stations and 86 are classified as satisfactory calibrators. We show that a flat low-frequency spectrum (from 74 to 325 MHz) is the best predictor of compactness at 140 MHz. We extrapolate from our sample to show that the density of calibrators on the sky that are sufficiently bright to calibrate dispersive and non-dispersive delays for the International LOFAR using existing methods is 1.0 per square degree. Conclusions. The observed density of satisfactory delay calibrator sources means that observations with International LOFAR should be possible at virtually any point in the sky, provided that a fast and efficient search using the methodology described here is conducted prior to the observation to identify the best calibrator.

Context. The radio galaxy IC 310 has recently been identified as a gamma-ray emitter based on observations at GeV energies with Fermi-LAT and at very high energies (VHE, E>100GeV) with the MAGIC telescopes. Originally classified as a head-tail radio galaxy, the nature of this object is subject of controversy since its nucleus shows blazar-like behavior. Aims. In order to understand the nature of IC 310 and the origin of the VHE emission we studied the spectral and flux variability of IC 310 from the X-ray band to the VHE gamma-ray regime. Methods. The light curve of IC 310 above 300GeV has been measured with the MAGIC telescopes from Oct. 2009 to Feb. 2010. Fermi-LAT data (2008-2011) in the 10-500GeV energy range were also analyzed. In X-ray, archival observations from 2003 to 2007 with XMM, Chandra, and Swift-XRT in the 0.5-10keV band were studied. Results. The VHE light curve reveals several high-amplitude and short-duration flares. Day-to-day flux variability is clearly present. The photon index between 120GeV and 8TeV remains at the value $\Gamma\sim2.0$ during both low and high flux states. The VHE spectral shape does not show significant variability, whereas the flux at 1TeV changes by a factor of $\sim7$. Fermi-LAT detected only eight gamma-ray events in the energy range 10GeV-500GeV in three years of observation. The measured photon index of $\Gamma=1.3\pm0.5$ in the Fermi-LAT range is very hard. The X-ray measurements show strong variability in flux and photon index. The latter varied from $1.76\pm0.07$ to $2.55\pm0.07$. Conclusion. The rapid variability measured confirms the blazar-like behavior of IC 310. The TeV emission seems to originate from scales of less than 80 Schwarzschild radii within the compact core of its FRI radio jet with orientation angle 10deg-38deg. The SED resembles that of an extreme blazar, albeit the luminosity is more than two orders of magnitude lower.

We report on broad multi-wavelength observations of the 2010-2011 periastron passage of the gamma-ray loud binary system PSR B1259-63. High resolution interferometric radio observations establish extended radio emission trailing the position of the pulsar. Observations with the Fermi Gamma-ray Space Telescope reveal GeV gamma-ray flaring activity of the system, reaching the spin-down luminosity of the pulsar, around 30 days after periastron. There are no clear signatures of variability at radio, X-ray and TeV energies at the time of the GeV flare. Variability around periastron in the H$\alpha$ emission line, can be interpreted as the gravitational interaction between the pulsar and the circumstellar disk. The equivalent width of the H$\alpha$ grows from a few days before periastron until a few days later, and decreases again between 18 and 46 days after periastron. In near infrared we observe the similar decrease of the equivalent width of Br$\gamma$ line between the 40th and 117th day after the periastron. For the idealized disk, the variability of the H$\alpha$ line represents the variability of the mass and size of the disk. We discuss possible physical relations between the state of the disk and GeV emission under assumption that GeV flare is directly related to the decrease of the disk size.

Several binary systems, composed of a star and a compact object, have been detected in the GeV-TeV range. Several systems have been observed but only a handful of sources have shown emission at those energies. Here, we present the observations conducted by MAGIC of different {\gamma}-ray binary systems. On one hand, we show the latest studies on the binary system LS I +61 303, which displays variability on different timescales. With the latest MAGIC observations, we will try to shed light on our understanding of this source, by presenting super-orbital and multi-wavelength studies. On the other hand, we show the observational results on the binary system HD 215227. This source has been proposed as a new {\gamma}-ray binary for being spatially coincident with the gamma-ray source AGL J2241+4454 detected by AGILE at E >100 GeV.

Gamma-ray binaries allow us to study physical processes such as particle acceleration up to very-high energies and gamma-ray emission and absorption with changing geometrical configurations on a periodic basis. These sources produce outflows of radio-emitting particles whose structure can be imaged with Very Long Baseline Interferometry (VLBI). We have studied the changing morphology of the gamma-ray binaries LS I +61 303 and LS 5039, and we have discovered the extended emission of PSR B1259-63 and HESS J0632+057. Based on these results, we have established the basic properties and behaviour of the radio emission of gamma-ray binaries on AU scales, and we have contributed to find characteristics that are common to all of them. Here we present the most relevant properties of each source and the general properties of gamma-ray binaries, and we describe the implications on the nature of these binary systems.

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

Motivated by the Costamante & Ghisellini (2002) predictions we investigated if the blazar 1ES 1727+502 (z=0.055) is emitting very high energy (VHE, E>100 GeV) gamma rays. We observed the BL Lac object 1ES 1727+502 in stereoscopic mode with the two MAGIC telescopes during 14 nights between May 6th and June 10th 2011, for a total effective observing time of 12.6 hours. For the study of the multiwavelength spectral energy distribution (SED) we use simultaneous optical R-band data from the KVA telescope, archival UV/optical and X-ray observations by instruments UVOT and XRT on board of the Swift satellite and high energy (HE, 0.1 GeV - 100 GeV) gamma-ray data from the Fermi-LAT instrument. We detect, for the first time, VHE gamma-ray emission from 1ES 1727+502 at a statistical significance of 5.5 sigma. The integral flux above 150 GeV is estimated to be (2.1\pm0.4)% of the Crab Nebula flux and the de-absorbed VHE spectrum has a photon index of (2.7\pm0.5). No significant short-term variability was found in any of the wavebands presented here. We model the SED using a one-zone synchrotron self-Compton model obtaining parameters typical for this class of sources.

Magnetars are an extreme, highly magnetized class of isolated neutron stars whose large X-ray luminosity is believed to be driven by their high magnetic field. In this work we study for the first time the possible very high energy gamma-ray emission above 100 GeV from magnetars, observing the sources 4U 0142+61 and 1E 2259+586. We observed the two sources with atmospheric Cherenkov telescopes in the very high energy range (E > 100 GeV). 4U 0142+61 was observed with the MAGIC I telescope in 2008 for ~25 h and 1E 2259+586 was observed with the MAGIC stereoscopic system in 2010 for ~14 h. The data were analyzed with the standard MAGIC analysis software. Neither magnetar was detected. Upper limits to the differential and integral flux above 200 GeV were computed using the Rolke algorithm. We obtain integral upper limits to the flux of 1.52*10^-12cm^-2 s^-1 and 2.7*10^-12cm^-2 s^-1 with a confidence level of 95% for 4U 0142+61 and 1E 2259+586, respectively. The resulting differential upper limits are presented together with X-ray data and upper limits in the GeV energy range.

Gamma-ray binaries allow us to study physical processes such as particle acceleration up to very-high energies and gamma-ray emission and absorption with changing geometrical configurations on a periodic basis. They produce outflows of radio-emitting particles whose structure can be imaged with Very Long Baseline Interferometry (VLBI). We present recent and new VLBI observations of PSR B1259-63, LS 5039, LS I +61 303, and HESS J0632+057. For the first three cases the results show the repeatability of their radio structures with the orbit of the binary system.

Gamma-ray binaries allow us to study physical processes such as particle acceleration up to TeV energies and VHE gamma-ray emission and absorption with changing geometrical configurations on a periodic basis. These sources produce outflows of radio-emitting particles whose structure can be imaged with VLBI. LS 5039 is a gamma-ray binary that has shown variable VLBI structures in the past. We aim to characterise the radio morphological changes of LS 5039 and discriminate if they are either repeatable or erratic. We observed LS 5039 with the VLBA at 5 GHz during five consecutive days to cover the 3.9-day orbit and an extra day to disentangle between orbital or secular variability. We also compiled the available high-resolution radio observations of the source to study its morphological variability at different orbital phases. We used a simple model to interpret the obtained images. The new observations show that the morphology of LS 5039 up to projected distances of 10 milliarcseconds changes in 24 h. The observed radio morphological changes display a periodic orbital modulation. Multifrequency and multiepoch VLBI observations confirm that the morphological periodicity is stable on timescales of years. Using a simple model we show that the observed behaviour is compatible with the presence of a young non-accreting pulsar with an outflow behind it. The morphology is reproduced for inclinations of the orbit of 60-75 deg. For masses of the companion star in the range 20-50 Msun, this range of inclinations implies a mass of the compact object of 1.3-2.7 Msun. The periodic orbital modulation of the radio morphology of LS 5039 suggests that all gamma-ray binaries are expected to show a similar behaviour. The changes in the radio structure of LS 5039 are compatible with the presence of a young non-accreting neutron star, which suggests that the known gamma-ray binaries contain young pulsars.

Context. We present the discovery of very high energy (VHE, E > 100GeV) gamma-ray emission from the BL Lac object 1ES 1215+303 by the MAGIC telescopes and simultaneous multi-wavelength data in a broad energy range from radio to gamma-rays. Aims. We study the VHE gamma-ray emission from 1ES 1215+303 and its relation to the emissions in other wavelengths. Methods. Triggered by an optical outburst, MAGIC observed the source in January-February 2011 for 20.3 hrs. The target was monitored in the optical R-band by the KVA telescope that also performed optical polarization measurements. We triggered target of opportunity observations with the Swift satellite and obtained simultaneous and quasi-simultaneous data from the Fermi Large Area Telescope and from the Mets\"ahovi radio telescope. We also present the analysis of older MAGIC data taken in 2010. Results. The MAGIC observations of 1ES 1215+303 carried out in January-February 2011 resulted in the first detection of the source at VHE with a statistical significance of 9.4 sigma. Simultaneously, the source was observed in a high optical and X-ray state. In 2010 the source was observed in a lower state in optical, X-ray, and VHE, while the GeV gamma-ray flux and the radio flux were comparable in 2010 and 2011. The spectral energy distribution obtained with the 2011 data can be modeled with a simple one zone SSC model, but it requires extreme values for the Doppler factor or the electron energy distribution.

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.

The high frequency peaked BL Lac PKS 2155-304 with a redshift of z=0.116 was discovered in 1997 in the very high energy (VHE, E >100GeV) gamma-ray range by the University of Durham Mark VI gamma-ray Cherenkov telescope in Australia with a flux corresponding to 20% of the Crab Nebula flux. It was later observed and detected with high significance by the Southern Cherenkov observatory H.E.S.S. Detection from the Northern hemisphere is difficult due to challenging observation conditions under large zenith angles. In July 2006, the H.E.S.S. collaboration reported an extraordinary outburst of VHE gamma-emission. During the outburst, the VHE gamma-ray emission was found to be variable on the time scales of minutes and with a mean flux of ~7 times the flux observed from the Crab Nebula. Follow-up observations with the MAGIC-I standalone Cherenkov telescope were triggered by this extraordinary outburst and PKS 2155-304 was observed between 28 July to 2 August 2006 for 15 hours at large zenith angles. Here we present our studies on the behavior of the source after its extraordinary flare and an enhanced analysis method for data taken at high zenith angles. We developed improved methods for event selection that led to a better background suppression. The averaged energy spectrum we derived has a spectral index of -3.5 +/- 0.2 above 400GeV, which is in good agreement with the spectral shape measured by H.E.S.S. during the major flare on MJD 53944. Furthermore, we present the spectral energy distribution modeling of PKS 2155-304. With our observations we increased the duty cycle of the source extending the light curve derived by H.E.S.S. after the outburst. Finally, we find night-by-night variability with a maximal amplitude of a factor three to four and an intranight variability in one of the nights (MJD 53945) with a similar amplitude.

Context. LS I +61 303 is a member of the select group of gamma-ray binaries: galactic binary systems that contain a massive star and a compact object, show a changing milliarcsecond morphology and a similar broad spectral energy distribution (SED) that peaks at MeV-TeV energies and is modulated by the orbital motion. The nature of the compact object is unclear in LS I +61 303, LS 5039 and HESS J0632+057, whereas PSR B1259-63 harbours a 47.74 ms radio pulsar. Aims. A scenario in which a young pulsar wind interacts with the stellar wind has been proposed to explain the very high energy (VHE, E > 100 GeV) gamma-ray emission detected from LS I +61 303, although no pulses have been reported from this system at any wavelength. We aim to find evidence of the pulsar nature of the compact object. Methods. We performed phased array observations with the Giant Metrewave Radio Telescope (GMRT) at 1280 MHz centred at phase 0.54. Simultaneous data from the multi-bit phased array (PA) back-end with a sampling time of tsamp = 128 microsec and from the polarimeter (PMT) back-end with tsamp = 256 microsec where taken. Results. No pulses have been found in the data set, with a minimum detectable mean flux density of \sim 0.38 mJy at 8-sigma level for the pulsed emission from a putative pulsar with period P >2 ms and duty cycle D = 10% in the direction of LS I +61 303. Conclusions. The detection of posible radio pulsations will require deep and sensitive observations at frequencies \sim0.5-5 GHz and orbital phases 0.6-0.7. However, it may be unfeasible to detect pulses if the putative pulsar is not beamed at the Earth or if there is a strong absorption within the binary system.

Context. The gamma-ray binary LS 5039 and the isolated pulsar PSR J1825-1446 were proposed to have been formed in the supernova remnant (SNR) G016.8-01.1. Aims. We aim to obtain the Galactic trajectory of LS 5039 and PSR J1825-1446 to find their origin in the Galaxy, and in particular to check their association with SNR G016.8-01.1 to restrict their age. Methods. By means of radio and optical observations we obtained the proper motion and the space velocity of the sources. Results. The proper motion of PSR J1825-1446 corresponds to a transverse space velocity of 690 km/s at a distance of 5 kpc. Its Galactic velocity at different distances is not compatible with the expected Galactic rotation. The velocity and characteristic age of PSR J1825-1446 make it incompatible with SNR G016.8-01.1. There are no clear OB associations or SNRs crossing the past trajectory of PSR J1825-1446. We estimate the age of the pulsar to be 80-245 kyr, which is compatible with its characteristic age. The proper motion of LS 5039 is 7.09 and -8.82 mas/yr in right ascension and declination, respectively. The association of LS 5039 with SNR G016.8-01.1 is unlikely, although we cannot to discard it. The system would have had to be formed in the association Ser OB2 (at 2.0 kpc) if the age of the system is 1.0-1.2 Myr, or in the association Sct OB3 (distance 1.5-2 kpc) for an age of 0.1-0.2 Myr. If the system were not formed close to Ser OB2, the pseudo-synchronization of the orbit would be unlikely. Conclusions. PSR J1825-1446 is a high-velocity isolated pulsar ejected from the Galaxy. The distance to LS 5039, which needs to be constrained by future astrometric missions such as Gaia, is a key parameter for restricting its origin and age.

Galaxy clusters are being assembled today in the most energetic phase of hierarchical structure formation which manifests itself in powerful shocks that contribute to a substantial energy density of cosmic rays (CRs). Hence, clusters are expected to be luminous gamma-ray emitters since they also act as energy reservoirs for additional CR sources, such as active galactic nuclei and supernova-driven galactic winds. To detect the gamma-ray emission from CR interactions with the ambient cluster gas, we conducted the deepest to date observational campaign targeting a galaxy cluster at very high-energy gamma-rays and observed the Perseus cluster with the MAGIC Cherenkov telescopes for a total of ~85 hr of effective observing time. This campaign resulted in the detection of the central radio galaxy NGC 1275 at energies E > 100 GeV with a very steep energy spectrum. Here, we restrict our analysis to energies E > 630 GeV and detect no significant gamma-ray excess. This constrains the average CR-to-thermal pressure ratio to be <= 1-2%, depending on assumptions and the model for CR emission. Comparing these gamma-ray upper limits to predictions from cosmological cluster simulations that include CRs constrains the maximum CR acceleration efficiency at structure formation shocks to be < 50%. Alternatively, this may argue for non-negligible CR transport processes such as CR streaming and diffusion into the outer cluster regions. Finally, we derive lower limits on the magnetic field distribution assuming that the Perseus radio mini-halo is generated by secondary electrons/positrons that are created in hadronic CR interactions: assuming a spectrum of E^-2.2 around TeV energies as implied by cluster simulations, we limit the central magnetic field to be > 4-9 microG, depending on the rate of decline of the magnetic field strength toward larger radii.

The very high energy (VHE) \gamma-ray source HESS J0632+057 has recently been confirmed to be a \gamma-ray binary. The optical counterpart is the Be star MWC 148, and a compact object of unknown nature orbits it every ~321 d with a high eccentricity of ~0.8. We monitored HESS J0632+057 with the stereoscopic MAGIC telescopes from 2010 October to 2011 March and detected significant VHE \gamma-ray emission during 2011 February, when the system exhibited an X-ray outburst. We find no \gamma-ray signal in the other observation periods when the system did not show increased X-ray flux. Thus HESS J0632+057 exhibits \gamma-ray variability on timescales of the order of one to two months possibly linked to the X-ray outburst that takes place about 100 days after the periastron passage. Furthermore our measurements provide for the first time the \gamma-ray spectrum down to about 140 GeV and indicate no turnover of the spectrum at low energies. We compare the properties of HESS J0632+057 with the similar \gamma-ray binary LS I +61 303, and discuss on the possible origin of the multi-wavelength emission of the source

We report on the detection of very-high energy (VHE, E>100 GeV) gamma-ray emission from NGC 1275, the central radio galaxy of the Perseus cluster of galaxies. The source has been detected by the MAGIC telescopes with a statistical significance of 6.6 sigma above 100 GeV in 46 hr of stereo observations carried out between August 2010 and February 2011. The measured differential energy spectrum between 70 GeV and 500 GeV can be described by a power law with a steep spectral index of \Gamma=-4.1+/-0.7stat+/-0.3syst, and the average flux above 100 GeV is F_{gamma}=(1.3+/-0.2stat+/-0.3syst) x 10^-11 cm^-2 s^-1. These results, combined with the power-law spectrum measured in the first two years of observations by the Fermi-LAT above 100 MeV, with a spectral index of Gamma ~ -2.1, strongly suggest the presence of a break or cut-off around tens of GeV in the NGC 1275 spectrum. The light curve of the source above 100 GeV does not show hints of variability on a month time scale. Finally, we report on the nondetection in the present data of the radio galaxy IC 310, previously discovered by the Fermi-LAT and MAGIC. The derived flux upper limit F^{U.L.}_{gamma} (>300 GeV)=1.2 x 10^-12 cm^-2 s^-1 is a factor ~ 3 lower than the mean flux measured by MAGIC between October 2009 and February 2010, thus confirming the year time-scale variability of the source at VHE.

We use 73 h of stereoscopic data taken with the MAGIC telescopes to investigate the very high-energy (VHE) gamma-ray emission of the Crab pulsar. Our data show a highly significant pulsed signal in the energy range from 50 to 400 GeV in both the main pulse (P1) and the interpulse (P2) phase regions. We provide the widest spectra to date of the VHE components of both peaks, and these spectra extend to the energy range of satellite-borne observatories. The good resolution and background rejection of the stereoscopic MAGIC system allows us to cross-check the correctness of each spectral point of the pulsar by comparison with the corresponding (strong and well-known) Crab nebula flux. The spectra of both P1 and P2 are compatible with power laws with photon indices of 4.0 \pm 0.8 (P1) and 3.42 \pm 0.26 (P2), respectively, and the ratio P1/P2 between the photon counts of the two pulses is 0.54 \pm 0.12. The VHE emission can be understood as an additional component produced by the inverse Compton scattering of secondary and tertiary e\pm pairs on IR-UV photons.

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.

The W51 complex hosts the supernova remnant W51C which is known to interact with the molecular clouds in the star forming region W51B. In addition, a possible pulsar wind nebula CXO J192318.5+140305 was found likely associated with the supernova remnant. Gamma-ray emission from this region was discovered by Fermi/LAT (between 0.2 and 50 GeV) and H.E.S.S. (>1 TeV). The spatial distribution of the events could not be used to pinpoint the location of the emission among the pulsar wind nebula, the supernova remnant shell and/or the molecular cloud. However, the modeling of the spectral energy distribution presented by the Fermi/LAT collaboration suggests a hadronic emission mechanism. We performed observations of the W51 complex with the MAGIC telescopes for more than 50 hours. The good angular resolution in the medium (few hundred GeV) to high (above 1 TeV) energies allow us to perform morphological studies. We detect an extended emission of very-high-energy gamma rays, with a significance of 11 standard deviations. We extend the spectrum from the highest Fermi/LAT energies to \sim 5 TeV and find that it follows a single power law with an index of 2.58 \pm 0.07stat \pm 0.22syst . The main part of the emission coincides with the shocked cloud region, while we find a feature extending towards the pulsar wind nebula. The possible contribution of the pulsar wind nebula, assuming a point-like source, shows no dependence on energy and it is about 20% of the overall emission. The broad band spectral energy distribution can be explained with a hadronic model that implies proton acceleration above 100 TeV. This result, together with the morphology of the source, tentatively suggests that we observe ongoing acceleration of ions in the interaction zone between supernova remnant and cloud. These results shed light on the long-standing problem of the origin of galactic cosmic rays.

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.

MAGIC is a system of two Imaging Atmospheric Cherenkov Telescopes located in the Canary island of La Palma. Since autumn 2009 both telescopes have been working together in stereoscopic mode, providing a significant improvement with respect to the previous single-telescope observations. We use observations of the Crab Nebula taken at low zenith angles to assess the performance of the MAGIC stereo system. The trigger threshold of the MAGIC telescopes is 50-60 GeV. Advanced stereo analysis techniques allow MAGIC to achieve a sensitivity as good as (0.76 +/- 0.03)% of the Crab Nebula flux in 50 h of observations above 290 GeV. The angular resolution at those energies is better than ~0.07 degree. We also perform a detailed study of possible systematic effects which may influence the analysis of the data taken with the MAGIC telescopes.

We report on the observation of $\gamma$-rays above 25\,GeV from the Crab pulsar (PSR B0532+21) using the MAGIC I telescope. Two data sets from observations during the winter period 2007/2008 and 2008/2009 are used. In order to discuss the spectral shape from 100\,MeV to 100\,GeV, one year of public {\it Fermi} Large Area Telescope ({\it Fermi}-LAT) data are also analyzed to complement the MAGIC data. The extrapolation of the exponential cutoff spectrum determined with the Fermi-LAT data is inconsistent with MAGIC measurements, which requires a modification of the standard pulsar emission models. In the energy region between 25 and 100\,GeV, the emission in the P1 phase (from -0.06 to 0.04, location of the main pulse) and the P2 phase (from 0.32 to 0.43, location of the interpulse) can be described by power laws with spectral indices of $-3.1 \pm 1.0_{stat} \pm 0.3_{syst}$ and $-3.5 \pm 0.5_{stat} \pm 0.3_{syst}$, respectively. Assuming an asymmetric Lorentzian for the pulse shape, the peak positions of the main pulse and the interpulse are estimated to be at phases $-0.009 \pm 0.007$ and $0.393 \pm 0.009$, while the full widths at half maximum are $0.025 \pm 0.008$ and $0.053 \pm 0.015$, respectively.

Compilation of the papers contributed by the MAGIC collaboration to the 32nd International Cosmic Ray Conference (ICRC 2011), which took place between August 11 and 18, 2011 in Beijing, China. The papers are sorted in 6 categories: Overview and Highlight papers; Instrument, software and techniques; Galactic sources; Extragalactic sources; Multi-wavelength and joint campaigns; Fundamental physics, dark matter and cosmic rays.

Context. After the detection of a 321-days periodicity in X-rays, HESS J0632+057 can be robustly considered a new member of the selected group of gamma-ray binaries. These sources are known to show extended radio structure at scales of milliarcseconds (mas). Aims. We present the expected extended radio emission on mas scales from HESS J0632+057. Methods. We observed HESS J0632+057 with the European VLBI Network (EVN) at 1.6 GHz in two epochs: during the January/February 2011 X-ray outburst and 30 days later. Results. The VLBI image obtained during the outburst shows a compact ~0.4 mJy radio source, whereas 30 days later the source has faded and appears extended, with a projected size of ~75 AU. The peak of the emission is displaced between runs 21+/-5 AU, which is bigger than the orbit size. The position of the radio source is compatible with the Be star MWC 148, which sets the proper motion of the binary system below 3 mas yr^-1 in each coordinate. The brightness temperature of the source is above 2 \times 10^6 K. We compare the multiwavelength properties of HESS J0632+057 with those of the previously known gamma-ray binaries. Conclusions. HESS J0632+057 displays extended and variable non-thermal radio emission. Its morphology, size, and displacement at AU scales are similar to those found in the other gamma-ray binaries, PSR B1259-63, LS 5039 and LS I +61 303, supporting a similar nature for HESS J0632+057.

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.

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.

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.

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.

The discovery of emission of TeV gamma rays from X-ray binaries has triggered an intense effort to better understand the particle acceleration, absorption, and emission mechanisms in compact binary systems. Here we present the pioneering effort of the MAGIC collaboration to understand the very high energy emission of the prototype system LS I +61 303. We report on the variable nature of the emission from LS I +61 303 and show that this emission is indeed periodic. The system shows regular outburst at TeV energies in phase phi=0.6-0.7 and detect no signal at periastron (phi~ 0.275). Furthermore we find no indication of spectral variation along the orbit of the compact object and the spectral energy distribution is compatible with a simple power law with index Gamma=2.6+-0.2_(stat)+-0.2_(sys). To answer some of the open questions concerning the emission process of the TeV radiation we conducted a multiwavelength campaign with the MAGIC telescope, XMM-Newton, and Swift in September 2007. We detect a simultaneous outburst at X-ray and TeV energies, with the peak at phase 0.62 and a similar shape at both wavelengths. A linear fit to the strictly simultaneous X-ray/TeV flux pairs provides r=0.81 -0.21 +0.06. Here we present the observations and discuss the implications of the obtained results to the emission processes in the system.