This white paper describes the science case for Very Long Baseline Interferometry (VLBI) and provides suggestions towards upgrade paths for the European VLBI Network (EVN). The EVN is a distributed long-baseline radio interferometric array, that operates at the very forefront of astronomical research. Recent results, together with the new science possibilities outlined in this vision document, demonstrate the EVN's potential to generate new and exciting results that will transform our view of the cosmos. Together with e-MERLIN, the EVN provides a range of baseline lengths that permit unique studies of faint radio sources to be made over a wide range of spatial scales. The science cases are reviewed in six chapters that cover the following broad areas: cosmology, galaxy formation and evolution, innermost regions of active galactic nuclei, explosive phenomena and transients, stars and stellar masers in the Milky Way, celestial reference frames and space applications. The document concludes with identifying the synergies with other radio, as well as multi-band/multi-messenger instruments, and provide the recommendations for future improvements. The appendices briefly describe other radio VLBI arrays, the technological framework for EVN developments, and a selection of spectral lines of astrophysical interest below 100 GHz. The document includes a glossary for non-specialists, and a list of acronyms at the end.

Fast radio bursts (FRBs) are bright, millisecond-duration radio transients originating from extragalactic distances. Their origin is unknown. Some FRB sources emit repeat bursts, ruling out cataclysmic origins for those events. Despite searches for periodicity in repeat burst arrival times on time scales from milliseconds to many days, these bursts have hitherto been observed to appear sporadically, and though clustered, without a regular pattern. Here we report the detection of a $16.35\pm0.15$ day periodicity (or possibly a higher-frequency alias of that periodicity) from a repeating FRB 180916.J0158+65 detected by the Canadian Hydrogen Intensity Mapping Experiment Fast Radio Burst Project (CHIME/FRB). In 38 bursts recorded from September 16th, 2018 through February 4th, 2020, we find that all bursts arrive in a 5-day phase window, and 50% of the bursts arrive in a 0.6-day phase window. Our results suggest a mechanism for periodic modulation either of the burst emission itself, or through external amplification or absorption, and disfavour models invoking purely sporadic processes.

Infrared imaging of the colliding-wind binary Apep has revealed a spectacular dust plume with complicated internal dynamics that challenges standard colliding-wind binary physics. Such challenges can be potentially resolved if a rapidly-rotating Wolf-Rayet star is located at the heart of the system, implicating Apep as a Galactic progenitor system to long-duration gamma-ray bursts. One of the difficulties in interpreting the dynamics of Apep is that the spectral composition of the stars in the system was unclear. Here we present visual to near-infrared spectra that demonstrate that the central component of Apep is composed of two classical Wolf-Rayet stars of carbon- (WC8) and nitrogen-sequence (WN4-6b) subtypes. We argue that such an assignment represents the strongest case of a classical WR+WR binary system in the Milky Way. The terminal line-of-sight wind velocities of the WC8 and WN4-6b stars are measured to be $2100 \pm 200$ and $3500 \pm 100$ km s$^{-1}$, respectively. If the mass-loss rate of the two stars are typical for their spectral class, the momentum ratio of the colliding winds is expected to be $\approx$ 0.4. Since the expansion velocity of the dust plume is significantly smaller than either of the measured terminal velocities, we explore the suggestion that one of the Wolf-Rayet winds is anisotropic. We can recover a shock-compressed wind velocity consistent with the observed dust expansion velocity if the WC8 star produces a significantly slow equatorial wind with a velocity of $\approx$530 km s$^{-1}$. Such slow wind speeds can be driven by near-critical rotation of a Wolf-Rayet star.

We report on simultaneous radio and X-ray observations of the repeating fast radio burst source FRB 180916.J0158+65 using the Canadian Hydrogen Intensity Mapping Experiment (CHIME), Effelsberg, and Deep Space Network (DSS-14 and DSS-63) radio telescopes and the Chandra X-ray Observatory. During 33 ks of Chandra observations, we detect no radio bursts in overlapping Effelsberg or Deep Space Network observations and a single radio burst during CHIME/FRB source transits. We detect no X-ray events in excess of the background during the Chandra observations. These non-detections imply a 5-$\sigma$ limit of $<5\times10^{-10}$ erg cm$^{-2}$ for the 0.5--10 keV fluence of prompt emission at the time of the radio burst and $1.3\times10^{-9}$ erg cm$^{-2}$ at any time during the Chandra observations at the position of FRB 180916.J0158+65. Given the host-galaxy redshift of FRB 180916.J0158+65 ($z\sim0.034$), these correspond to energy limits of $<1.6\times10^{45}$ erg and $<4\times10^{45}$ erg, respectively. We also place a 5-$\sigma$ limit of $<8\times10^{-15}$ erg s$^{-1}$ cm$^{-2}$ on the 0.5--10\,keV absorbed flux of a persistent source at the location of FRB 180916.J0158+65. This corresponds to a luminosity limit of $<2\times10^{40}$ erg s$^{-1}$. Using Fermi/GBM data we search for prompt gamma-ray emission at the time of radio bursts from FRB 180916.J0158+65 and find no significant bursts, placing a limit of $4\times10^{-9}$ erg cm$^{-2}$ on the 10--100 keV fluence. We also search Fermi/LAT data for periodic modulation of the gamma-ray brightness at the 16.35-day period of radio-burst activity and detect no significant modulation. We compare these deep limits to the predictions of various fast radio burst models, but conclude that similar X-ray constraints on a closer fast radio burst source would be needed to strongly constrain theory.

Fast radio bursts (FRBs) are brief, bright, extragalactic radio flashes. Their physical origin remains unknown, but dozens of possible models have been postulated. Some FRB sources exhibit repeat bursts. Though over a hundred FRB sources have been discovered to date, only four have been localised and associated with a host galaxy, with just one of the four known to repeat. The properties of the host galaxies, and the local environments of FRBs, provide important clues about their physical origins. However, the first known repeating FRB has been localised to a low-metallicity, irregular dwarf galaxy, and the apparently non-repeating sources to higher-metallicity, massive elliptical or star-forming galaxies, suggesting that perhaps the repeating and apparently non-repeating sources could have distinct physical origins. Here we report the precise localisation of a second repeating FRB source, FRB 180916.J0158+65, to a star-forming region in a nearby (redshift $z = 0.0337 \pm 0.0002$) massive spiral galaxy, whose properties and proximity distinguish it from all known hosts. The lack of both a comparably luminous persistent radio counterpart and a high Faraday rotation measure further distinguish the local environment of FRB 180916.J0158+65 from that of the one previously localised repeating FRB source, FRB 121102. This demonstrates that repeating FRBs have a wide range of luminosities, and originate from diverse host galaxies and local environments.

This paper reports observations of a 22 GHz water maser `superburst' in the G25.65+1.05 massive star forming region, conducted in response to an alert from the Maser Monitoring Organisation (M2O). Very long baseline interferometry (VLBI) observations using the European VLBI Network (EVN) recorded a maser flux density of $1.2 \times 10^{4}$ Jy. The superburst was investigated in the spectral, structural and temporal domains and its cause was determined to be an increase in maser path length generated by the superposition of multiple maser emitting regions aligning in the line of sight to the observer. This conclusion was based on the location of the bursting maser in the context of the star forming region, its complex structure, and its rapid onset and decay.

Fast radio bursts (FRBs) are an enigmatic class of extragalactic transients emitting Jy-level radio bursts in the GHz band, lasting for only a few ms. So far, some objects are known to repeat while several others are not, likely indicating multiple origins. There are many theoretical models, some predict prompt VHE or optical emission correlated with FRBs while others imply VHE afterglows hours after the FRB. To test these predictions and unravel the nature of FRB progenitors, the stereoscopic Imaging Atmospheric Cherenkov Telescopes (IACTs) system MAGIC has been participating in FRB observation campaigns since 2016. As IACTs are sensitive to Cherenkov photons in the UV/blue region of the electromagnetic spectrum and use photo-detectors with time response faster than a ms, MAGIC is also able to perform simultaneous optical observations through a dedicated system installed in the central PMT of its camera. The main challenge faced by MAGIC in searching for optical counterpart of FRBs is the presence of irreducible background optical events due to terrestrial sources. We present new results from MAGIC observations of the first repeating FRB 121102 during several MWL observation campaigns. The recently improved instrument and refined strategy to search for counterparts of FRBs in the VHE and optical bands will also be presented.

Pulsars typically exhibit radio emission in the form of narrow pulses originated from confined regions of their magnetospheres. A potential presence of magnetospherically originated emission outside this region, the so-called off-pulse emission, would challenge the existing theories. Detection of significant off-pulse emission has been reported so far from only two pulsars, B0525+21 and B2045-16, at 325 and 610 MHz. However, the nature of this newly uncovered off-pulse emission remains unclear. To probe its origin we conducted very high resolution radio observations of B0525+21 and B2045-16 with the European VLBI Network (EVN) at 1.39 GHz. Whereas the pulsed emission is detected at a level consistent with previous observations, we report absence of any off-pulse emission above $42$ and $96\ \mathrm{\mu Jy\ beam^{-1}}$ (three times the rms noise levels) for B0525+21 and B2045-16, respectively. Our stringent upper limits imply the off-pulse emission to be less than $0.4$ and $0.3\%$ of the period-averaged pulsed flux density, i.e., much fainter than the previously suggested values of $1$-$10\%$. Since the EVN data are most sensitive to extremely compact angular scales, our results suggest a non-magnetospheric origin for the previously reported off-pulse emission. Presence of extended emission that is resolved out on these milliarcsecond scales still remains plausible. In this case, we constrain the emission to arise from structures with sizes of $\sim (0.61$-$19) \times 10^3\ \mathrm{au}$ for B0525+21 and $\sim (0.48$-$8.3) \times 10^3\ \mathrm{au}$ for B2045-16. These constraints might indicate that the two pulsars are accompanied by compact bow-shock pulsar wind nebulae. Future observations probing intermediate angular scales ($\sim 0.1$-$5\ \mathrm{arcsec}$) will help in clarifying the actual origin of the off-pulse emission.

Ofek (2017) identified FIRST J141918.9+394036 (hereafter FIRST J1419+3940) as a radio source sharing similar properties and host galaxy type to the compact, persistent radio source associated with the first known repeating fast radio burst, FRB 121102. Law et al. (2018) showed that FIRST J1419+3940 is a transient source decaying in brightness over the last few decades. One possible interpretation is that FIRST J1419+3940 is a nearby analogue to FRB 121102 and that the radio emission represents a young magnetar nebula (as several scenarios assume for FRB 121102). Another interpretation is that FIRST J1419+3940 is the afterglow of an `orphan' long gamma-ray burst (GRB). The environment is similar to where most such events are produced. To distinguish between these hypotheses, we conducted radio observations using the European VLBI Network at 1.6 GHz to spatially resolve the emission and to search for millisecond-duration radio bursts. We detect FIRST J1419+3940 as a compact radio source with a flux density of $620 \pm 20\ \mathrm{\mu Jy}$ (on 2018 September 18) and a source size of $3.9 \pm 0.7\ \mathrm{mas}$ (i.e. $1.6 \pm 0.3\ \mathrm{pc}$ given the angular diameter distance of $83\ \mathrm{Mpc}$). These results confirm that the radio emission is non-thermal and imply an average expansion velocity of $(0.10 \pm 0.02)c$. Contemporaneous high-time-resolution observations using the 100-m Effelsberg telescope detected no millisecond-duration bursts of astrophysical origin. The source properties and lack of short-duration bursts are consistent with a GRB jet expansion, whereas they disfavor a magnetar birth nebula.

The binary neutron star merger event GW170817 was detected through both electromagnetic radiation and gravitational waves. Its afterglow emission may have been produced by either a narrow relativistic jet or an isotropic outflow. High spatial resolution measurements of the source size and displacement can discriminate between these scenarios. We present Very Long Baseline Interferometry observations, performed 207.4 days after the merger, using a global network of 32 radio telescopes. The apparent source size is constrained to be smaller than 2.5 milliarcseconds at the 90% confidence level. This excludes the isotropic outflow scenario, which would have produced a larger apparent size, indicating that GW170817 produced a structured relativistic jet. Our rate calculations show that at least 10% of neutron star mergers produce such a jet.

Fast Radio Bursts (FRBs) are transient sources that emit a single radio pulse with a duration of only a few milliseconds. Since the discovery of the first FRB in 2007, tens of similar events have been detected. However, their physical origin remains unclear, and a number of scenarios even larger than the number of known FRBs has been proposed during these years. The presence of repeating bursts in FRB 121102 allowed us to perform a precise localization of the source with the Very Large Array and the European VLBI Network (EVN). Optical observations with Keck, Gemini and HST unveiled the host to be a low-metallicity star-forming dwarf galaxy located at a redshift of 0.193. The EVN results showed that the bursts are co-located (within a projected separation of $< 40$ pc) to a compact and persistent radio source with a size of $< 0.7$ pc inside a star-forming region. This environment resembles the ones where superluminous supernovae (SLSNe) or long gamma-ray bursts are produced. Although the nature of this persistent source and the origin of the bursts remain unknown, scenarios considering a pulsar/magnetar energizing a young SLSN, or a system with a pulsar/magnetar in the vicinity of a massive black hole are the most plausible ones to date. More recent observations have shown that the bursts from FRB 121102 are almost 100% linearly polarized at an unexpectedly high and variable Faraday rotation measure, that has been observed to date only in vicinities of massive black holes. The bursts are thus likely produced from a neutron star in such environment, although the system can still be explained by a young neutron star embedded in a highly magnetized nebula. Upcoming interferometric searches are expected to report tens of these localizations in the coming years, unveil if this source is representative of the whole population or a particular case, and dramatically boosting the field of FRBs.

In a rare and spectacular display, two well-known massive star forming regions, W49N and G25.65+1.05, recently underwent maser 'super burst' - their fluxes suddenly increasing above 30,000 and 18,000 Jy, respectively, reaching several orders of magnitude above their usual values. In quick-response, ToO observations with the EVN, VLBA and KaVA were obtained constituting a 4 week campaign - producing a high-cadence multi-epoch VLBI investigation of the maser emission. The combination of high-resolution, polarisation and flux monitoring during the burst provides one of the best accounts, to date, of the maser super burst phenomenon, aiding their use as astrophysical tools. These proceedings contain the preliminary results of our campaign.

Aims: To investigate the extension of the very-high-energy spectral tail of the Crab pulsar at energies above 400 GeV. Methods: We analyzed $\sim$320 hours of good quality data of Crab with the MAGIC telescope, obtained from February 2007 until April 2014. Results: We report the most energetic pulsed emission ever detected from the Crab pulsar reaching up to 1.5 TeV. The pulse profile shows two narrow peaks synchronized with the ones measured in the GeV energy range. The spectra of the two peaks follow two different power-law functions from 70 GeV up to 1.5 TeV and connect smoothly with the spectra measured above 10 GeV by the Large Area Telescope (LAT) on board of the Fermi satellite. When making a joint fit of the LAT and MAGIC data, above 10 GeV, the photon indices of the spectra differ by 0.5$\pm$0.1. Conclusions: We measured with the MAGIC telescopes the most energetic pulsed photons from a pulsar to date. Such TeV pulsed photons require a parent population of electrons with a Lorentz factor of at least $5\times 10^6$. These results strongly suggest IC scattering off low energy photons as the emission mechanism and a gamma-ray production region in the vicinity of the light cylinder.

FRB 121102 is the only known repeating fast radio burst source. Here we analyze a wide-frequency-range (1-8 GHz) sample of high-signal-to-noise, coherently dedispersed bursts detected using the Arecibo and Green Bank telescopes. These bursts reveal complex time-frequency structures that include sub-bursts with finite bandwidths. The frequency-dependent burst structure complicates the determination of a dispersion measure (DM); we argue that it is appropriate to use a DM metric that maximizes frequency-averaged pulse structure, as opposed to peak signal-to-noise, and find DM = 560.57 +/- 0.07 pc/cc at MJD 57644. After correcting for dispersive delay, we find that the sub-bursts have characteristic frequencies that typically drift lower at later times in the total burst envelope. In the 1.1-1.7 GHz band, the ~ 0.5-1-ms sub-bursts have typical bandwidths ranging from 100-400 MHz, and a characteristic drift rate of ~ 200 MHz/ms towards lower frequencies. At higher radio frequencies, the sub-burst bandwidths and drift rate are larger, on average. While these features could be intrinsic to the burst emission mechanism, they could also be imparted by propagation effects in the medium local to the source. Comparison of the burst DMs with previous values in the literature suggests an increase of Delta(DM) ~ 1-3 pc/cc in 4 years, though this could be a stochastic variation as opposed to a secular trend. This implies changes in the local medium or an additional source of frequency-dependent delay. Overall, the results are consistent with previously proposed scenarios in which FRB 121102 is embedded in a dense nebula.

Gamma-ray binaries are systems composed of a massive star and a compact object that exhibit emission from radio to very high energy gamma rays. They are ideal laboratories to study particle acceleration and a variety of physical processes that vary as a function of the orbital phase. We aim to study the radio emission of the gamma-ray binary 1FGL J1018.6-5856 to constrain the emitting region and determine the peculiar motion of the system within the Galaxy to clarify its origin by analyzing an observation with the Australian Long Baseline Array at 8.4~GHz. We combined these data with the optical Gaia DR2 and UCAC4 catalogs to consolidate the astrometry information therein. 1FGL J1018.6-5856 shows compact radio emission ($< 3$ mas or $\lesssim 20$ au at $\sim 6.4$ kpc distance), implying a brightness temperature of $\gtrsim 5.6 \times 10^6$ K, and confirming its nonthermal origin. We report consistent results between the proper motion reported by Gaia DR2 and the positions obtained from the Gaia, UCAC4, and LBA data. We also determined the distance to 1FGL J1018.6-5856 to be $6.4_{-0.7}^{+ 1.7}$ kpc. Together with the radial velocity of the source we computed its three-dimensional proper and peculiar motion within the Galaxy. We obtained a peculiar motion on its regional standard of rest (RSR) frame of $|u| = 45_{-9}^{+30}$ km s$^{-1}$, with the system moving away from the Galactic plane. In the simplest scenario of a symmetric stellar core collapse we estimate a mass loss of $ 4 \lesssim \Delta M \lesssim 9$ M$_{\odot}$ during the creation of the compact object. 1FGL J1018.6-5856 exhibits compact radio emission similar to that detected in other gamma-ray binaries. We provide the first accurate peculiar motion estimations of the system and place it within the Galaxy. The obtained motion and distance excludes its physical relation with the supernova remnant G284.3-1.8.

Fast radio bursts (FRBs) are bright flashes observed typically at GHz frequencies with millisecond duration, whose origin is likely extragalactic. Their nature remains mysterious, motivating searches for counterparts at other wavelengths. FRB 121102 is so far the only source known to repeatedly emit FRBs and is associated with a host galaxy at redshift $z \simeq 0.193$. We conducted simultaneous observations of FRB 121102 with the Arecibo and MAGIC telescopes during several epochs in 2016--2017. This allowed searches for millisecond-timescale burst emission in very-high-energy (VHE) gamma rays as well as the optical band. While a total of five FRBs were detected during these observations, no VHE emission was detected, neither of a persistent nature nor burst-like associated with the FRBs. The average integral flux upper limits above 100 GeV at 95% confidence level are $6.6 \times 10^{-12}~\mathrm{photons\ cm^{-2}\ s^{-1}}$ (corresponding to luminosity $L_{\rm VHE} \lesssim 10^{45}~\mathrm{erg\ s^{-1}}$) over the entire observation period, and $1.2 \times 10^{-7}~ \mathrm{photons\ cm^{-2}\ s^{-1}}$ ($L_{\rm VHE} \lesssim 10^{49}~\mathrm{erg\ s^{-1}}$) over the total duration of the five FRBs. We constrain the optical U-band flux to be below 8.6 mJy at 5-$\sigma$ level for 1-ms intervals around the FRB arrival times. A bright burst with U-band flux $29~\mathrm{mJy}$ and duration $\sim 12$ ms was detected 4.3 s before the arrival of one FRB. However, the probability of spuriously detecting such a signal within the sampled time space is 1.5% (2.2 $\sigma$, post-trial), i.e. consistent with the expected background. We discuss the implications of the obtained upper limits for constraining FRB models.

e-ASTROGAM (enhanced ASTROGAM) is a breakthrough Observatory space mission, with a detector composed by a Silicon tracker, a calorimeter, and an anticoincidence system, dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV - the lower energy limit can be pushed to energies as low as 150 keV for the tracker, and to 30 keV for calorimetric detection. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LIGO-Virgo-GEO600-KAGRA, SKA, ALMA, E-ELT, TMT, LSST, JWST, Athena, CTA, IceCube, KM3NeT, and LISA.

Gamma-ray binaries are systems composed of a massive star and a compact object that produce emission from radio to very high energy gamma-rays. LS I +61 303 is one of the only six gamma-ray binaries discovered so far. It is thought that gamma-ray binaries contain a young highly rotating neutron star as compact object, and the emission is produced by the interaction between its relativistic pulsar wind and the stellar wind, However, in the case of LS I +61 303 a microquasar scenario is still considered and results pointing to oppose directions have been published during the last decades. Here we provide a review about the state of the art of LS I +61 303, summarizing the observed emission from radio to very high energy gamma-rays along all these years, and we discuss the proposed scenarios that can explain such emission.

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

Fast radio bursts (FRBs) are millisecond-duration, extragalactic radio flashes of unknown physical origin. FRB 121102, the only known repeating FRB source, has been localized to a star-forming region in a dwarf galaxy at redshift z = 0.193, and is spatially coincident with a compact, persistent radio source. The origin of the bursts, the nature of the persistent source, and the properties of the local environment are still debated. Here we present bursts that show ~100% linearly polarized emission at a very high and variable Faraday rotation measure in the source frame: RM_src = +1.46 x 10^5 rad m^-2 and +1.33 x 10^5 rad m^-2 at epochs separated by 7 months, in addition to narrow (< 30 mus) temporal structure. The large and variable rotation measure demonstrates that FRB 121102 is in an extreme and dynamic magneto-ionic environment, while the short burst durations argue for a neutron star origin. Such large rotation measures have, until now, only been observed in the vicinities of massive black holes (M_BH > 10^4 MSun). Indeed, the properties of the persistent radio source are compatible with those of a low-luminosity, accreting massive black hole. The bursts may thus come from a neutron star in such an environment. However, the observed properties may also be explainable in other models, such as a highly magnetized wind nebula or supernova remnant surrounding a young neutron star.

We undertook coordinated campaigns with the Green Bank, Effelsberg, and Arecibo radio telescopes during Chandra X-ray Observatory and XMM-Newton observations of the repeating fast radio burst FRB 121102 to search for simultaneous radio and X-ray bursts. We find 12 radio bursts from FRB 121102 during 70 ks total of X-ray observations. We detect no X-ray photons at the times of radio bursts from FRB 121102 and further detect no X-ray bursts above the measured background at any time. We place a 5$\sigma$ upper limit of $3\times10^{-11}$ erg cm$^{-2}$ on the 0.5--10 keV fluence for X-ray bursts at the time of radio bursts for durations $<700$ ms, which corresponds to a burst energy of $4\times10^{45}$ erg at the measured distance of FRB 121102. We also place limits on the 0.5--10 keV fluence of $5\times10^{-10}$ erg cm$^{-2}$ and $1\times10^{-9}$ erg cm$^{-2}$ for bursts emitted at any time during the XMM-Newton and Chandra observations, respectively, assuming a typical X-ray burst duration of 5 ms. We analyze data from the Fermi Gamma-ray Space Telescope Gamma-ray Burst Monitor and place a 5$\sigma$ upper limit on the 10--100 keV fluence of $4\times10^{-9}$ erg cm$^{-2}$ ($5\times10^{47}$ erg at the distance of FRB 121102) for gamma-ray bursts at the time of radio bursts. We also present a deep search for a persistent X-ray source using all of the X-ray observations taken to date and place a 5$\sigma$ upper limit on the 0.5--10 keV flux of $4\times10^{-15}$ erg s$^{-1}$ cm$^{-2}$ ($3\times10^{41}$ erg~s$^{-1}$ at the distance of FRB 121102). We discuss these non-detections in the context of the host environment of FRB 121102 and of possible sources of fast radio bursts in general.

We present optical, near- and mid-infrared imaging of the host galaxy of FRB 121102 with the Gemini North telescope, the Hubble Space Telescope and the Spitzer Space Telescope. The FRB 121102 host galaxy is resolved, revealing a bright star forming region located in the outskirts of the irregular, low-metallicity dwarf galaxy. The star forming region has a half-light radius of 0.68 kpc (0.20 arcsec), encompassing the projected location of the compact (<0.7 pc), persistent radio source that is associated with FRB 121102. The half-light diameter of the dwarf galaxy is 5 to 7 kpc, and broadband spectral energy distribution fitting indicates that it has a total stellar mass of M*~10^8 Msun. The metallicity of the host galaxy is low, 12+log10 ([O/H])=8.0+-0.1. The properties of the host galaxy of FRB 121102 are comparable to those of extreme emission line galaxies, also known to host hydrogen-poor superluminous supernovae and long-duration gamma-ray bursts. The projected location of FRB 121102 within the star forming region supports the proposed connection of FRBs with newly born neutron stars or magnetars.

The binary system AR Sco contains an M star and the only known radio-pulsing white dwarf. The system shows emission from radio to X-rays, likely dominated by synchrotron radiation. The mechanism that produces most of this emission remains unclear. Two competing scenarios have been proposed: Collimated outflows, and direct interaction between the magnetospheres of the white dwarf and the M star. The two proposed scenarios can be tested via very long baseline interferometric radio observations. We conducted a radio observation with the Australian Long Baseline Array (LBA) on 20 Oct 2016 at 8.5 GHz to study the compactness of the radio emission. Simultaneous data with the Australian Telescope Compact Array (ATCA) were also recorded for a direct comparison of the obtained flux densities. AR Sco shows radio emission compact on milliarcsecond angular scales ($\lesssim 0.02\ \mathrm{AU}$, or $4\ \mathrm{R_{\odot}}$). The emission is orbitally modulated, with an average flux density of $\approx 6.5\ \mathrm{mJy}$. A comparison with the simultaneous ATCA data shows that no flux is resolved out on mas scales, implying that the radio emission is produced in this compact region. Additionally, the obtained radio light curves on hour timescales are consistent with the optical light curve. The radio emission in AR Sco is likely produced in the magnetosphere of the M star or the white dwarf, and we see no evidence for a radio outflow or collimated jets significantly contributing to the radio emission.

Fast radio bursts, or FRBs, are transient sources of unknown origin. Recent radio and optical observations have provided strong evidence for an extragalactic origin of the phenomenon and the precise localization of the repeating FRB 121102. Observations using the Karl G. Jansky Very Large Array (VLA) and very-long-baseline interferometry (VLBI) have revealed the existence of a continuum non-thermal radio source consistent with the location of the bursts in a dwarf galaxy. All these new data rule out several models that were previously proposed, and impose stringent constraints to new models. We aim to model FRB 121102 in light of the new observational results in the active galactic nucleus (AGN) scenario. We propose a model for repeating FRBs in which a non-steady relativistic $e^\pm$-beam, accelerated by an impulsive magnetohydrodynamic (MHD)-driven mechanism, interacts with a cloud at the centre of a star-forming dwarf galaxy. The interaction generates regions of high electrostatic field called cavitons in the plasma cloud. Turbulence is also produced in the beam. These processes, plus particle isotropization, the interaction scale, and light retardation effects, provide the necessary ingredients for short-lived, bright coherent radiation bursts. The mechanism studied in this work explains the general properties of FRB 121102, and may also be applied to other repetitive FRBs. Coherent emission from electrons and positrons accelerated in cavitons provides a plausible explanation of FRBs.

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

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

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.

The precise localization of the repeating fast radio burst (FRB 121102) has provided the first unambiguous association (chance coincidence probability $p\lesssim3\times10^{-4}$) of an FRB with an optical and persistent radio counterpart. We report on optical imaging and spectroscopy of the counterpart and find that it is an extended ($0.6^{\prime\prime}-0.8^{\prime\prime}$) object displaying prominent Balmer and [OIII] emission lines. Based on the spectrum and emission line ratios, we classify the counterpart as a low-metallicity, star-forming, $m_{r^\prime} = 25.1$ AB mag dwarf galaxy at a redshift of $z=0.19273(8)$, corresponding to a luminosity distance of 972 Mpc. From the angular size, the redshift, and luminosity, we estimate the host galaxy to have a diameter $\lesssim4$ kpc and a stellar mass of $M_*\sim4-7\times 10^{7}\,M_\odot$, assuming a mass-to-light ratio between 2 to 3$\,M_\odot\,L_\odot^{-1}$. Based on the H$\alpha$ flux, we estimate the star formation rate of the host to be $0.4\,M_\odot\,\mathrm{yr^{-1}}$ and a substantial host dispersion measure depth $\lesssim 324\,\mathrm{pc\,cm^{-3}}$. The net dispersion measure contribution of the host galaxy to FRB 121102 is likely to be lower than this value depending on geometrical factors. We show that the persistent radio source at FRB 121102's location reported by Marcote et al (2017) is offset from the galaxy's center of light by $\sim$200 mas and the host galaxy does not show optical signatures for AGN activity. If FRB 121102 is typical of the wider FRB population and if future interferometric localizations preferentially find them in dwarf galaxies with low metallicities and prominent emission lines, they would share such a preference with long gamma ray bursts and superluminous supernovae.

The millisecond-duration radio flashes known as Fast Radio Bursts (FRBs) represent an enigmatic astrophysical phenomenon. Recently, the sub-arcsecond localization (~ 100mas precision) of FRB121102 using the VLA has led to its unambiguous association with persistent radio and optical counterparts, and to the identification of its host galaxy. However, an even more precise localization is needed in order to probe the direct physical relationship between the millisecond bursts themselves and the associated persistent emission. Here we report very-long-baseline radio interferometric observations using the European VLBI Network and the 305-m Arecibo telescope, which simultaneously detect both the bursts and the persistent radio emission at milliarcsecond angular scales and show that they are co-located to within a projected linear separation of < 40pc (< 12mas angular separation, at 95% confidence). We detect consistent angular broadening of the bursts and persistent radio source (~ 2-4mas at 1.7GHz), which are both similar to the expected Milky Way scattering contribution. The persistent radio source has a projected size constrained to be < 0.7pc (< 0.2mas angular extent at 5.0GHz) and a lower limit for the brightness temperature of T_b > 5 x 10^7K. Together, these observations provide strong evidence for a direct physical link between FRB121102 and the compact persistent radio source. We argue that a burst source associated with a low-luminosity active galactic nucleus or a young neutron star energizing a supernova remnant are the two scenarios for FRB121102 that best match the observed data.

Fast radio bursts are astronomical radio flashes of unknown physical nature with durations of milliseconds. Their dispersive arrival times suggest an extragalactic origin and imply radio luminosities orders of magnitude larger than any other kind of known short-duration radio transient. Thus far, all FRBs have been detected with large single-dish telescopes with arcminute localizations, and attempts to identify their counterparts (source or host galaxy) have relied on contemporaneous variability of field sources or the presence of peculiar field stars or galaxies. These attempts have not resulted in an unambiguous association with a host or multi-wavelength counterpart. Here we report the sub-arcsecond localization of FRB 121102, the only known repeating burst source, using high-time-resolution radio interferometric observations that directly image the bursts themselves. Our precise localization reveals that FRB 121102 originates within 100 mas of a faint 180 uJy persistent radio source with a continuum spectrum that is consistent with non-thermal emission, and a faint (25th magnitude) optical counterpart. The flux density of the persistent radio source varies by tens of percent on day timescales, and very long baseline radio interferometry yields an angular size less than 1.7 mas. Our observations are inconsistent with the fast radio burst having a Galactic origin or its source being located within a prominent star-forming galaxy. Instead, the source appears to be co-located with a low-luminosity active galactic nucleus or a previously unknown type of extragalactic source. [Truncated] If other fast radio bursts have similarly faint radio and optical counterparts, our findings imply that direct sub-arcsecond localizations of FRBs may be the only way to provide reliable associations.

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

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

Blazars are variable emitters across all wavelengths over a wide range of timescales, from months down to minutes. It is therefore essential to observe blazars simultaneously at different wavelengths, especially in the X-ray and gamma-ray bands, where the broadband spectral energy distributions usually peak. In this work, we report on three "target-of-opportunity" (ToO) observations of Mrk 421, one of the brightest TeV blazars, triggered by a strong flaring event at TeV energies in 2014. These observations feature long, continuous, and simultaneous exposures with XMM-Newton (covering X-ray and optical/ultraviolet bands) and VERITAS (covering TeV gamma-ray band), along with contemporaneous observations from other gamma-ray facilities (MAGIC and Fermi-LAT) and a number of radio and optical facilities. Although neither rapid flares nor significant X-ray/TeV correlation are detected, these observations reveal subtle changes in the X-ray spectrum of the source over the course of a few days. We search the simultaneous X-ray and TeV data for spectral hysteresis patterns and time delays, which could provide insight into the emission mechanisms and the source properties (e.g. the radius of the emitting region, the strength of the magnetic field, and related timescales). The observed broadband spectra are consistent with a one-zone synchrotron self-Compton model. We find that the power spectral density distribution at $\gtrsim 4\times 10^{-4}$ Hz from the X-ray data can be described by a power-law model with an index value between 1.2 and 1.8, and do not find evidence for a steepening of the power spectral index (often associated with a characteristic length scale) compared to the previously reported values at lower frequencies.

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

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

We investigate the nature of the compact, and possibly variable nuclear radio source in the centre of WISE J0716-19, the proposed host galaxy of fast radio burst, FRB 150418. We observed WISE J0716-19 at 5.0 GHz with the European VLBI Network four times between 2016 March 16 and June 2. At three epochs, we simultaneously observed the source with e-MERLIN at the same frequency. We detected a compact source in the EVN data in each epoch with a significance up to ~8sigma. The four epochs yielded consistent results within their uncertainties, for both peak surface intensity and positions. The mean values for these quantities are I_peak = (115+-9) {\mu}Jy/beam and r.a. = 07:16:34.55496(7), dec. = -19:00:39.4754(8), respectively. The e-MERLIN data provided ~3-5sigma detections, at a position consistent with those of the EVN data. The presence of emission on angular scales intermediate between the EVN and e-MERLIN is consistent with being null. The brightness temperature of the EVN core is Tb~10^8.5K, close to the value required by Akiyama & Johnson (2016) to explain the radio properties of WISE J0716-19 in terms of interstellar induced short-term variability. Our observations provide direct, independent evidence of the existence of a nuclear compact source in WISE J0716-19, a physical scenario with no evident connection with FRB 150418. However, the EVN data do not show indication of the variability observed with the VLA.

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

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

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

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

High-mass binary systems involve extreme environments that produce non-thermal emission from radio to gamma rays. Only three types of these systems are known to emit persistent gamma-ray emission: colliding-wind binaries, high-mass X-ray binaries and gamma-ray binaries. This thesis is focused on the radio emission of high-mass binary systems through interferometric observations, and we have explored several of these sources with low- and high-frequency radio observations, and very high-resolution VLBI ones. We have studied two gamma-ray binaries, LS 5039 and LS I +61 303, at low frequencies. We have obtained their light-curves and spectra, and we have determined the physical properties of their radio emitting regions. We have also studied the gamma-ray binary HESS J0632+057 through VLBI observations. A new colliding wind binary, HD 93129A, has been discovered through VLBI and optical observations. Finally, we have conducted radio observations of two sources that were candidates to be gamma-ray binaries.

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

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

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

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

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

LS I +61 303 is a gamma-ray binary that exhibits an outburst at GHz frequencies each orbital cycle of $\approx$ 26.5 d and a superorbital modulation with a period of $\approx$ 4.6 yr. We have performed a detailed study of the low-frequency radio emission of LS I +61 303 by analysing all the archival GMRT data at 150, 235 and 610 MHz, and conducting regular LOFAR observations within the Radio Sky Monitor (RSM) at 150 MHz. We have detected the source for the first time at 150 MHz, which is also the first detection of a gamma-ray binary at such a low frequency. We have obtained the light-curves of the source at 150, 235 and 610 MHz, all of them showing orbital modulation. The light-curves at 235 and 610 MHz also show the existence of superorbital variability. A comparison with contemporaneous 15-GHz data shows remarkable differences with these light-curves. At 15 GHz we see clear outbursts, whereas at low frequencies we see variability with wide maxima. The light-curve at 235 MHz seems to be anticorrelated with the one at 610 MHz, implying a shift of $\sim$ 0.5 orbital phases in the maxima. We model the shifts between the maxima at different frequencies as due to changes in the physical parameters of the emitting region assuming either free-free absorption or synchrotron self-absorption, obtaining expansion velocities for this region close to the stellar wind velocity with both mechanisms.

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

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

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

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

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

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

We have studied in detail the 0.15-15 GHz radio spectrum of the gamma-ray binary LS 5039 to look for a possible turnover and absorption mechanisms at low frequencies, and to constrain the physical properties of its emission. We have analysed two archival VLA monitorings, all the available archival GMRT data and a coordinated quasi-simultaneous observational campaign conducted in 2013 with GMRT and WSRT. The data show that the radio emission of LS 5039 is persistent on day, week and year timescales, with a variability $\lesssim 25~\%$ at all frequencies, and no signature of orbital modulation. The obtained spectra reveal a power-law shape with a curvature below 5 GHz and a turnover at $\sim0.5$ GHz, which can be reproduced by a one-zone model with synchrotron self-absorption plus Razin effect. We obtain a coherent picture for a size of the emitting region of $\sim0.85~\mathrm{mas}$, setting a magnetic field of $B\sim20~\mathrm{mG}$, an electron density of $n_{\rm e}\sim4\times10^5~{\rm cm^{-3}}$ and a mass-loss rate of $\dot M\sim5\times10^{-8}~{\rm M_{\odot} yr^{-1}}$. These values imply a significant mixing of the stellar wind with the relativistic plasma outflow from the compact companion. At particular epochs the Razin effect is negligible, implying changes in the injection and the electron density or magnetic field. The Razin effect is reported for first time in a gamma-ray binary, giving further support to the young non-accreting pulsar scenario.

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

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

Radio observations are an effective tool to discover particle acceleration regions in colliding-wind binaries, through detection of synchrotron radiation; these regions are natural laboratories for the study of relativistic particles. Wind-collision region (WCR) models can reproduce the radio continuum spectra of massive binaries that contain both thermal and non-thermal radio emission; however, key constraints for models come from high-resolution imaging. Only five WCRs have been resolved to date at radio frequencies at milliarcsec (mas) angular scales. The source HD 93129A, prototype of the very few known O2 I stars, is a promising target for study: recently, a second massive, early-type star about 50 mas away was discovered, and a non-thermal radio source detected in the region. Preliminary long-baseline array data suggest that a significant fraction of the radio emission from the system comes from a putative WCR. We sought evidence that HD 93129A is a massive binary system with colliding stellar winds that produce non-thermal radiation, through spatially resolved images of the radio emitting regions. We completed observations with the Australian Long Baseline Array (LBA) to resolve the system at mas angular resolutions and reduced archival Australia Telescope Compact Array (ATCA) data to derive the total radio emission. We also compiled optical astrometric data of the system in a homogeneous way. We reduced historical Hubble Space Telescope data and obtained absolute and relative astrometry with milliarcsec accuracy. The astrometric analysis leads us to conclude that the two stars in HD 93129A form a gravitationally bound system. The LBA data reveal an extended arc-shaped non-thermal source between the two stars, indicative of a WCR. The wind momentum-rate ratio of the two stellar winds is estimated. The ATCA data show a point source with a change in flux level ...

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

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

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

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

The number of known very high energy (VHE) blazars is $\sim\,50$, which is very small in comparison to the number of blazars detected in other frequencies. This situation is a handicap for population studies of blazars, which emit about half of their luminosity in the $\gamma$-ray domain. Moreover, VHE blazars, if distant, allow for the study of the environment that the high-energy $\gamma$-rays traverse in their path towards the Earth, like the extragalactic background light (EBL) and the intergalactic magnetic field (IGMF), and hence they have a special interest for the astrophysics community. We present the first VHE detection of 1ES\,0033+595 with a statistical significance of 5.5\,$\sigma$. The VHE emission of this object is constant throughout the MAGIC observations (2009 August and October), and can be parameterized with a power law with an integral flux above 150 GeV of $(7.1\pm1.3)\times 10^{-12} {\mathrm{ph\,cm^{-2}\,s^{-1}}}$ and a photon index of ($3.8\pm0.7$). We model its spectral energy distribution (SED) as the result of inverse Compton scattering of synchrotron photons. For the study of the SED we used simultaneous optical R-band data from the KVA telescope, archival X-ray data by \textit{Swift} as well as \textit{INTEGRAL}, and simultaneous high energy (HE, $300$\,MeV~--~$10$\,GeV) $\gamma$-ray data from the \textit{Fermi} LAT observatory. Using the empirical approach of Prandini et al. (2010) and the \textit{Fermi}-LAT and MAGIC spectra for this object, we estimate the redshift of this source to be $0.34\pm0.08\pm0.05$. This is a relevant result because this source is possibly one of the ten most distant VHE blazars known to date, and with further (simultaneous) observations could play an important role in blazar population studies, as well as future constraints on the EBL and IGMF.