We report the discovery and panchromatic followup observations of the young Type Ic supernova, SN2020oi, in M100, a grand design spiral galaxy at a mere distance of $14$ Mpc. We followed up with observations at radio, X-ray and optical wavelengths from only a few days to several months after explosion. The optical behaviour of the supernova is similar to those of other normal Type Ic supernovae. The event was not detected in the X-ray band but our radio observation revealed a bright mJy source ($L_{\nu} \approx 1.2 \times 10^{27} {\rm erg\,s}^{-1} {\rm Hz}^{-1}$). Given, the relatively small number of stripped envelope SNe for which radio emission is detectable, we used this opportunity to perform a detailed analysis of the comprehensive radio dataset we obtained. The radio emitting electrons initially experience a phase of inverse Compton cooling which leads to steepening of the spectral index of the radio emission. Our analysis of the cooling frequency points to a large deviation from equipartition at the level of $\epsilon_e/\epsilon_B \gtrsim 200$, similar to a few other cases of stripped envelope SNe. Our modeling of the radio data suggests that the shockwave driven by the SN ejecta into the circumstellar matter (CSM) is moving at $\sim 3\times 10^{4}\,{\rm km\,s}^{-1}$. Assuming a constant mass-loss from the stellar progenitor, we find that the mass-loss rate is $\dot{M} \approx 1.4\times 10^{-4}\,{M}_{\odot}\,{\rm yr}^{-1}$, for an assumed wind velocity of $1000\,{\rm km\,s}^{-1}$. The temporal evolution of the radio emission suggests a radial CSM density structure steeper than the standard $r^{-2}$.

We present the discovery with Keck of the extremely infrared (IR) luminous transient AT 2017gbl, coincident with the Northern nucleus of the luminous infrared galaxy (LIRG) IRAS 23436+5257. Our extensive multi-wavelength follow-up spans ~900 days, including photometry and spectroscopy in the optical and IR, and (very long baseline interferometry) radio and X-ray observations. Radiative transfer modelling of the host galaxy spectral energy distribution shows the presence of a hitherto undetected dust obscured active galactic nucleus (AGN). The optical and near-IR spectra show broad 2000 km/s hydrogen, He I and O I emission features that decrease in flux over time. Radio imaging shows a fast evolving compact source of synchrotron emission spatially coincident with AT 2017gbl. We infer a lower limit for the radiated energy of 7.3 x 10^50 erg from the IR photometry. An extremely energetic supernova would satisfy this budget, but is ruled out by the radio counterpart evolution. Instead, we propose AT 2017gbl is related to an accretion event by the central supermassive black hole, where the spectral signatures originate in the AGN broad line region and the IR photometry is consistent with re-radiation by polar dust. Given the fast evolution of AT 2017gbl, we deem a tidal disruption event (TDE) of a star a more plausible scenario than a dramatic change in the AGN accretion rate. This makes AT 2017gbl the third TDE candidate to be hosted by a LIRG, in contrast to the so far considered TDE population discovered at optical wavelengths and hosted preferably by post-starburst galaxies.

We present an overview and description of the eMERLIN Galaxy Evolution survey (eMERGE) Data Release 1 (DR1), a large program of high-resolution 1.5 GHz radio observations of the GOODS-N field comprising $\sim140$ hours of observations with eMERLIN and $\sim40$ hours with the Very Large Array (VLA). We combine the long baselines of eMERLIN (providing high angular resolution) with the relatively closely-packed antennas of the VLA (providing excellent surface brightness sensitivity) to produce a deep 1.5 GHz radio survey with the sensitivity ($\sim 1.5\mu$Jy beam$^{-1}$), angular resolution ($0.2"$--$0.7"$) and field-of-view ($\sim15' \times 15'$) to detect and spatially resolve star-forming galaxies and AGN at $z\gtrsim 1$. The goal of eMERGE is to provide new constraints on the deep, sub-arcsecond radio sky which will be surveyed by SKA1-mid. In this initial publication, we discuss our data analysis techniques, including steps taken to model in-beam source variability over a $\sim20$ year baseline and the development of new point spread function/primary beam models to seamlessly merge eMERLIN and VLA data in the $uv$ plane. We present early science results, including measurements of the luminosities and/or linear sizes of $\sim500$ galaxes selected at 1.5 GHz. In combination with deep Hubble Space Telescope observations, we measure a mean radio-to-optical size ratio of $r_{\rm eMERGE}/r_{\rm HST}\sim1.02\pm0.03$, suggesting that in most high-redshift galaxies, the $\sim$GHz continuum emission traces the stellar light seen in optical imaging. This is the first in a series of papers which will explore the $\sim$kpc-scale radio properties of star-forming galaxies and AGN in the GOODS-N field observed by eMERGE DR1.

The quest for binary and dual supermassive black holes (SMBHs) at the dawn of the multi-messenger era is compelling. Detecting dual active galactic nuclei (AGN) -- active SMBHs at projected separations larger than several parsecs -- and binary AGN -- probing the scale where SMBHs are bound in a Keplerian binary -- is an observational challenge. The study of AGN pairs (either dual or binary) also represents an overarching theoretical problem in cosmology and astrophysics. The AGN triggering calls for detailed knowledge of the hydrodynamical conditions of gas in the imminent surroundings of the SMBHs and, at the same time, their duality calls for detailed knowledge on how galaxies assemble through major and minor mergers and grow fed by matter along the filaments of the cosmic web. This review describes the techniques used across the electromagnetic spectrum to detect dual and binary AGN candidates and proposes new avenues for their search. The current observational status is compared with the state-of-the-art numerical simulations and models for formation of dual and binary AGN. Binary SMBHs are among the loudest sources of gravitational waves (GWs) in the Universe. The search for a background of GWs at nHz frequencies from inspiralling SMBHs at low redshifts, and the direct detection of signals from their coalescence by the Laser Interferometer Space Antenna in the next decade, make this a theme of major interest for multi-messenger astrophysics. This review discusses the future facilities and observational strategies that are likely to significantly advance this fascinating field.

We report deep radio observations of nearby Type Ia Supernovae (SNe Ia) with the electronic Multi-Element Radio Linked Interferometer Net-work (e-MERLIN), and the Australia Telescope Compact Array (ATCA). No detections were made. With standard assumptions for the energy densities of relativistic electrons going into a power-law energy distribution, and the magnetic field strength (epsilon_e = epsilon_B = 0.1), we arrive at the upper limits on mass-loss rate for the progenitor system of SN 2013dy (2016coj, 2018gv, 2018pv, 2019np), to be less than 12 (2.8,1.3, 2.1, 1.7)EE(-8) solar masses per year (for a wind speed of 100 km/s). To SNe 2016coj, 2018gv, 2018pv and 2019np we add radio data for 17 other nearby SNe Ia, and model their non-detections. With the same model as described, all 21 SNe Ia have mass-loss rates less than 4EE(-8) solar masses per year (for a wind speed of 100 km/s). We compare those limits with the expected mass loss rates in different single-degenerate progenitor scenarios. We also discuss how information on epsilon_e and epsilon_B can be obtained from late observations of SNe Ia and the youngest SN Ia remnant detected in radio, G1.9+0.3, as well as stripped-envelope core-collapse SNe. We highlight SN 2011dh, and argue for epsilon_e approximately equal to 0.1 and epsilon_B approximately equal to 0.0033. Finally, we discuss strategies to observe at radio frequencies to maximize the chance of detection, given the time since explosion, the distance to the supernova and the telescope sensitivity.

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.

Water fountains are evolved stars showing early stages of collimated mass loss during transition from the asymptotic giant branch, providing valuable insight into the formation of asymmetric planetary nebulae. We report the results of multi-epoch VLBI observations, which determine the spatial and three-dimensional kinematic structure of H2O masers associated with the water fountain IRAS 18113-2503. The masers trace three pairs of high-velocity (~150-300 km/s) bipolar bow shocks on a scale of 0.18'' (~2000 au). The expansion velocities of the bow shocks exhibit an exponential decrease as a function of distance from the central star, which can be explained by an episodic, jet-driven outflow decelerating due to drag forces in a circumstellar envelope. Using our model, we estimate an initial ejection velocity ~840 km/s, a period for the ejections ~10 yr, with the youngest being ~12 yr old, and an average envelope density within the H2O maser region n(H2) ~ 10^6 cm^(-3). We hypothesize that IRAS 18113-2503 hosts a binary central star with a separation of ~10 au, revealing novel clues about the launching mechanisms of high-velocity collimated outflows in water fountains.

We report a serendipitous hard X-ray (3-24 keV) detection of the blazar LBQS 1319+0039. The X-ray spectrum is consistent with powerlaw emission with a photon index of $\Gamma = 1.78$. The observed emission and published redshift imply an observed 2-10 keV luminosity of $3.4\times10^{45}$ erg s$^{-1}$. This is consistent with extrapolations from a previous soft X-ray detection, thus there is no evidence for strong variability.

We report spatially-resolved [CII]$\lambda 158$ $\mu$m observations of HE 0433-1028, which is the first detection of a nearby luminous AGN (redshift 0.0355) with FIFI-LS onboard the airborne observatory SOFIA. We compare the spatially-resolved star formation tracers [CII], as provided by our SOFIA observations, and H$\alpha$ from MUSE optical integral-field spectroscopy. We find that the [CII] emission is mainly matching the extended star formation as traced by the extinction-corrected H$\alpha$ line emission but some additional flux is present. While a larger sample is needed to statistically confirm our findings and investigate possible dependencies on AGN luminosity and star formation rate, our study underlines the necessity of collecting a spatially-resolved optical-FIR dataset for nearby AGNs, and shows that it is technically feasible to collect such datasets with FIFI-LS onboard SOFIA.

We compute the local rate of events similar to GRB 170817A, which has been recently found to be associated with a kilonova (KN) outburst. Our analysis finds an observed rate of such events of R$_{KN}\sim 352^{+810}_{-281}$ Gpc$^{-3}$yr$^{-1}$. After comparing at their face values this density of sGRB outbursts with the much higher density of Binary Neutron Star (BNS) mergers of 1540$^{+3200}_{-1220}$ Gpc$^{-3}$yr$^{-1}$, estimated by LIGO-Virgo collaboration, one can conclude, admittedly with large uncertainty that either only a minor fraction of BNS mergers produces sGRB/KN events or the sGRBs associated with BNS mergers are beamed and observable under viewing angles as large as $\theta$ $\leq$ $40^{\circ}$. Finally we provide preliminary estimates of the number of sGRB/KN events detected by future surveys carried out with present/future ground-based/space facilities, such as LSST, VST, ZTF, SKA and THESEUS.

We report on an optical (SDSS) and X-ray (XMM) study of an optically selected sample of four dual AGN systems at projected separations of 30--60~kpc. All sources are detected in the X-ray band (0.3-10 keV); seven objects are optically identified as Seyfert, while one source, optically classified as a LINER, is likely powered by accretion in virtue of its relatively high X-ray luminosity (1.2$\times10^{41}$ erg/s). Six of the eight objects are obscured in X-rays with N$_{\rm H} \geq$ 10$^{23}$ cm$^{-2}$; three of these, whose X-ray spectrum is dominated by a reflection component, are likely Compton-thick (N$_{\rm H} \geq$ 10$^{24}$ cm$^{-2}$). This finding is in agreement with the hypothesis that galaxy encounters are effective in driving gas inflow toward the nuclear region, thus increasing the obscuration. We compare the absorption properties in our dual AGN with those in larger samples observed in X-rays but selected in different ways (optical, IR and hard X-rays). We find that the obscured (N$_{\rm H} \geq$ 10$^{22}$ cm$^{-2}$) AGN fraction within the larger sample is 84$\pm$4 per cent (taking into account the 90 per cent error on the N$_{\rm H}$ measure) up to large pair separations ($\sim$100~kpc). This is statistically higher than the fraction of obscured AGN in isolated galaxies found in X-ray surveys.

ThunderKAT is the image-plane transients programme for MeerKAT. The goal as outlined in 2010, and still today, is to find, identify and understand high-energy astrophysical processes via their radio emission (often in concert with observations at other wavelengths). Through a comprehensive and complementary programme of surveying and monitoring Galactic synchrotron transients (across a range of compact accretors and a range of other explosive phenomena) and exploring distinct populations of extragalactic synchrotron transients (microquasars, supernovae and possibly yet unknown transient phenomena) - both from direct surveys and commensal observations - we will revolutionise our understanding of the dynamic and explosive transient radio sky. As well as performing targeted programmes of our own, we have made agreements with the other MeerKAT large survey projects (LSPs) that we will also search their data for transients. This commensal use of the other surveys, which remains one of our key programme goals in 2016, means that the combined MeerKAT LSPs will produce by far the largest GHz-frequency radio transient programme to date.

After changing optical AGN type from 1.9 to 1 in 1984, the AGN Mrk 1018 recently reverted back to its type 1.9 state. Our ongoing monitoring now reveals that the AGN has halted its dramatic dimming, reaching a minimum around October 2016. The minimum was followed by an outburst rising with $\sim$0.25 U-band mag/month. The rebrightening lasted at least until February 2017, as confirmed by joint Chandra and Hubble observations. Monitoring was resumed in July 2017 after the source emerged from sunblock, at which point the AGN was found only $\sim$0.4 mag brighter than its minimum. The intermittent outburst was accompanied by the appearance of a red wing asymmetry in broad-line shape, indicative of an inhomogeneous broad-line region. The current flickering brightness of Mrk 1018 following its rapid fading either suggests that the source has reignited, remains variable at a low level, or may continue dimming over the next few years. Distinguishing between these possibilities requires continuous multiwavelength monitoring.

We present a radio light curve of supernova (SN) 2014C taken with the Arcminute Microkelvin Imager (AMI) Large Array at 15.7 GHz. Optical observations presented by Milisavljevic et al. demonstrated that SN 2014C metamorphosed from a stripped-envelope Type Ib SN into a strongly interacting Type IIn SN within 1 year. The AMI light curve clearly shows two distinct radio peaks, the second being a factor of 4 times more luminous than the first peak. This double bump morphology indicates two distinct phases of mass-loss from the progenitor star with the transition between density regimes occurring at 100-200 days. This reinforces the interpretation that SN 2014C exploded in a low density region before encountering a dense Hydrogen-rich shell of circumstellar material that was likely ejected by the progenitor prior to the explosion. The AMI flux measurements of the first light curve bump are the only reported observations taken within ~50 to ~125 days post-explosion, before the blast-wave encountered the Hydrogen shell. Simplistic synchrotron self-absorption (SSA) and free-free absorption (FFA) modelling suggest that some physical properties of SN 2014C, such as the mass-loss rate, are consistent with the properties of other Type Ibc and IIn SNe. However, our single frequency data does not allow us to distinguish between these two models, which implies they are likely too simplistic to describe the complex environment surrounding this event. Lastly, we present the precise radio location of SN 2014C obtained with eMERLIN, which will be useful for future VLBI observations of the SN.

We present new IRAM 30m spectroscopic observations of the $\sim88$ GHz band, including emission from the CCH (n=1-0) multiplet, HCN (1-0), HCO+ (1-0), and HNC (1-0), for a sample of 58 local luminous and ultraluminous infrared galaxies from the Great Observatories All-sky LIRG Survey (GOALS). By combining our new IRAM data with literature data and Spitzer/IRS spectroscopy, we study the correspondence between these putative tracers of dense gas and the relative contribution of active galactic nuclei (AGN) and star formation to the mid-infrared luminosity of each system. We find the HCN (1-0) emission to be enhanced in AGN-dominated systems ($\langle$L'$_{HCN (1-0)}$/L'$_{HCO^+ (1-0)}\rangle=1.84$), compared to composite and starburst-dominated systems ($\langle$L'$_{HCN (1-0)}$/L'$_{HCO^+ (1-0)}\rangle=1.14$, and 0.88, respectively). However, some composite and starburst systems have L'$_{HCN (1-0)}$/L'$_{HCO^+ (1-0)}$ ratios comparable to those of AGN, indicating that enhanced HCN emission is not uniquely associated with energetically dominant AGN. After removing AGN-dominated systems from the sample, we find a linear relationship (within the uncertainties) between $\log_{10}$(L'$_{HCN (1-0)}$) and $\log_{10}$(L$_{IR}$), consistent with most previous findings. L'$_{HCN (1-0)}$/L$_{IR}$, typically interpreted as the dense gas depletion time, appears to have no systematic trend with L$_{IR}$ for our sample of luminous and ultraluminous infrared galaxies, and has significant scatter. The galaxy-integrated HCN (1-0) and HCO+ (1-0) emission do not appear to have a simple interpretation, in terms of the AGN dominance or the star formation rate, and are likely determined by multiple processes, including density and radiative effects.

In this paper we present an overview of the MAGNA (Multiple AGN Activity) project aiming at a comprehensive study of multiple supemassive black hole systems. With the main goal to characterize the sources in merging systems at different stages of evolution, we selected a sample of objects optically classified as multiple systems on the basis of emission line diagnostics and started a massive multiband observational campaign. Here we report on the discovery of the exceptionally high AGN density compact group SDSS~J0959+1259. A multiband study suggests that strong interactions are taking place among its galaxies through tidal forces, therefore this system represents a case study for physical mechanisms that trigger nuclear activity and star formation. We also present a preliminary analysis of the multiple AGN system SDSS~J1038+3921.}

We have carried out a high-precision astrometric analysis of two very-long-baseline-interferometry (VLBI) epochs of observation of the 13 extragalactic radio sources in the complete S5 polar cap sample. The VLBI epochs span a time baseline of 10 years and enable us to achieve precisions in the proper motions of the source cores up to a few micro-arcseconds per year. The observations were performed at 14.4 GHz and 43.1 GHz, and enable us to estimate the frequency core-shifts in a subset of sources, for which the spectral-index distributions can be computed. We study the source-position stability by analysing the changes in the relative positions of fiducial source points (the jet cores) over a decade. We find motions of 0.1-0.9 mas among close-by sources between the two epochs, which imply drifts in the jet cores of approximately a few tens of micro-as per year. These results have implications for the standard Active Galactic Nucleus (AGN) jet model (where the core locations are supposed to be stable in time). For one of our sources, 0615+820, the morphological and spectral properties in year 2010, as well as the relative astrometry between years 2000 and 2010, suggest the possibility of either a strong parsec-scale interaction of the AGN jet with the ISM, a gravitational lens with ~1 mas diameter, or a resolved massive binary black hole.

We report the discovery that the known `changing look' AGN Mrk 1018 has changed spectral type for a second time. New VLT-MUSE data taken in 2015 as part of the Close AGN Reference Survey (CARS) shows that the AGN has returned to its original Seyfert 1.9 classification. The CARS sample is selected to contain only bright type 1 AGN, but Mrk 1018's broad emission lines and continuum, typical of type 1 AGN, have almost entirely disappeared. We use spectral fitting of the MUSE spectrum and previously available spectra to determine the drop in broad line flux and the Balmer decrement. We find that the broad line flux has decreased by a factor of 4.75+/- 0.5 in H{\alpha} since an SDSS spectrum was taken in 2000. The Balmer decrement has not changed significantly implying no enhanced reddening with time, but the remaining broad lines are more asymmetric than those present in the type 1 phase. We posit that the change is due to an intrinsic drop in flux from the accretion disk rather than variable extinction or a tidal disruption event.

We recently discovered that the active galactic nucleus (AGN) of Mrk 1018 has changed optical type again after 30 years as a type 1 AGN. Here we combine Chandra, NuStar, Swift, Hubble Space Telescope and ground-based observations to explore the cause of this change. The 2-10keV flux declines by a factor of ~8 between 2010 and 2016. We show with our X-ray observation that this is not caused by varying neutral hydrogen absorption along the line-of-sight up to the Compton-thick level. The optical-UV spectral energy distributions are well fit with a standard geometrically thin optically thick accretion disc model that seems to obey the expected $L\sim T^4$ relation. It confirms that a decline in accretion disc luminosity is the primary origin for the type change. We detect a new narrow-line absorber in Lya blue-shifted by ~700km/s with respect to the systemic velocity of the galaxy. This new Lya absorber could be evidence for the onset of an outflow or a companion black hole with associated gas that could be related to the accretion rate change. However, the low column density of the absorber means that it is not the direct cause for Mrk 1018's changing-look nature.

We present a multi-wavelength study of a newly discovered compact group (CG), SDSS J0959+1259, based data from XMM-Newton, SDSS and the Calar Alto optical imager BUSCA. With a maximum velocity offset of 500 km s$^{-1}$, a mean redshift of 0.035, and a mean spatial extension of 480 kpc, this CG is exceptional in having the highest concentration of nuclear activity in the local Universe, established with a sensitivity limit L$_{X}>4\times $10$^{40}$ erg s$^{-1}$ in 2--10 keV band and R-band magnitude $M_R < -19$. The group is composed of two type-2 Seyferts, one type-1 Seyfert, two LINERs and three star forming galaxies. Given the high X-ray luminosity of LINERs which reaches $\sim 10^{41}$ erg s$^{-1}$, it is likely that they are also accretion driven, bringing the number of active nuclei in this group to to 5 out of 8 (AGN fraction of 60\%). The distorted shape of one member of the CG suggests that strong interactions are taking place among its galaxies through tidal forces. Therefore, this system represents a case study for physical mechanisms that trigger nuclear activity and star formation in CGs.

The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI). The Spanish SKA White Book is the result of the coordinated effort of 120 astronomers from 40 different research centers. The book shows the enormous scientific interest of the Spanish astronomical community in the SKA and warrants an optimum scientific exploitation of the SKA by Spanish researchers, if Spain enters the SKA project.

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of active galactic nuclei. For a summary, we refer to the paper.

The SKA will be a transformational instrument in the study of our local Universe. In particular, by virtue of its high sensitivity (both to point sources and diffuse low surface brightness emission), angular resolution and the frequency ranges covered, the SKA will undertake a very wide range of astrophysical research in the field of nearby galaxies. By surveying vast numbers of nearby galaxies of all types with $\mu$Jy sensitivity and sub-arcsecond angular resolutions at radio wavelengths, the SKA will provide the cornerstone of our understanding of star-formation and accretion activity in the local Universe. In this chapter we outline the key continuum and molecular line science areas where the SKA, both during phase-1 and when it becomes the full SKA, will have a significant scientific impact.

Aims. We want to study the amount of molecular gas in a sample of nearby early-type galaxies (ETGs) which host low-luminosity Active Galactic Nuclei (AGN). We look for possible differences between the radio-loud (RL) and radio-quiet (RQ) AGN. Methods. We observed the CO(1-0) and CO(2-1) spectral lines with the IRAM 30m and NRO 45m telescopes for eight galaxies. They belong to a large sample of 37 local ETGs which host both RQ and RL AGN. We gather data from the literature for the entire sample. Results. We report the new detection of CO(1-0) emission in four galaxies (UGC0968, UGC5617, UGC6946, and UGC8355) and CO(2-1) emission in two of them (UGC0968 and UGC5617). The CO(2-1)/CO(1-0) ratio in these sources is $\sim0.7\pm0.2$. Considering both the new observations and the literature, the detection rate of CO in our sample is 55 $\pm$ 9%, with no statistically significant difference between the hosts of RL and RQ AGNs. For all the detected galaxies we converted the CO luminosities into the molecular masses, $M_{H_2}$, that range from 10$^{6.5}$ to 10$^{8.5}$ M$_{\odot}$, without any statistically significant differences between RL and RQ galaxies. This suggests that the amount of molecular gas does not likely set the radio-loudness of the AGN. Furthermore, despite the low statistical significance, the presence of a weak trend between the H$_{2}$ mass with various tracers of nuclear activity (mainly [O III] emission line nuclear power) cannot be excluded.

We report deep EVN and eMERLIN observations of the Type Ia SN 2014J in the nearby galaxy M 82. Our observations represent, together with JVLA observations of SNe 2011fe and 2014J, the most sensitive radio studies of Type Ia SNe ever. By combining data and a proper modeling of the radio emission, we constrain the mass-loss rate from the progenitor system of SN 2014J to $\dot{M} \lesssim 7.0\times 10^{-10}\, {\rm M_{\odot}\, yr^{-1}}$ (3-$\sigma$; for a wind speed of $100\, {\rm km s^{-1}}$). If the medium around the supernova is uniform, then $n_{\rm ISM} \lesssim 1.3 {\rm cm^3}$ (3-$\sigma$), which is the most stringent limit for the (uniform) density around a Type Ia SN. Our deep upper limits favor a double-degenerate (DD) scenario--involving two WD stars--for the progenitor system of SN 2014J, as such systems have less circumstellar gas than our upper limits. By contrast, most single-degenerate (SD) scenarios, i.e., the wide family of progenitor systems where a red giant, main-sequence, or sub-giant star donates mass to a exploding WD, are ruled out by our observations. Our estimates on the limits to the gas density surrounding SN 2011fe, using the flux density limits from Chomiuk et al. (2012), agree well with their results. Although we discuss possibilities for a SD scenario to pass observational tests, as well as uncertainties in the modeling of the radio emission, the evidence from SNe 2011fe and 2014J points in the direction of a DD scenario for both.

We present near-infrared and optical photometry, plus optical spectroscopy of two stripped-envelope supernovae (SNe) 2010O and 2010P that exploded in two different components of an interacting luminous infrared galaxy Arp 299 within only a few days of one another. SN 2010O is found to be photometrically and spectroscopically similar to many normal Type Ib SNe and our multiwavelength observations of SN 2010P suggest it to be a Type IIb SN. No signs of clear hydrogen features or interaction with the circumstellar medium are evident in the optical spectrum of SN 2010P. We derive estimates for the host galaxy line-of-sight extinctions for both SNe, based on both light-curve and spectroscopic comparison finding consistent results. These methods are also found to provide much more robust estimates of the SN host galaxy reddening than the commonly used empirical relations between extinction and equivalent width of Na I D absorption features. The SN observations also suggest that different extinction laws are present in different components of Arp 299. For completeness, we study high-resolution pre-explosion images of Arp 299 and find both SNe to be close to, but not coincident with, extended sources that are likely massive clusters. A very simple model applied to the bolometric light curve of SN 2010O implies a rough estimate for the explosion parameters of $E_{\mathrm{k}} \approx 3 \times 10^{51}$ erg, $M_{\mathrm{ej}} \approx 2.9$ M$_{\odot}$ and $M_{\mathrm{Ni}} \approx 0.16$ M$_{\odot}$.

The Javalambre-Physics of the Accelerated Universe Astrophysical Survey (J-PAS) is a narrow band, very wide field Cosmological Survey to be carried out from the Javalambre Observatory in Spain with a purpose-built, dedicated 2.5m telescope and a 4.7 sq.deg. camera with 1.2Gpix. Starting in late 2015, J-PAS will observe 8500sq.deg. of Northern Sky and measure $0.003(1+z)$ photo-z for $9\times10^7$ LRG and ELG galaxies plus several million QSOs, sampling an effective volume of $\sim 14$ Gpc$^3$ up to $z=1.3$ and becoming the first radial BAO experiment to reach Stage IV. J-PAS will detect $7\times 10^5$ galaxy clusters and groups, setting constrains on Dark Energy which rival those obtained from its BAO measurements. Thanks to the superb characteristics of the site (seeing ~0.7 arcsec), J-PAS is expected to obtain a deep, sub-arcsec image of the Northern sky, which combined with its unique photo-z precision will produce one of the most powerful cosmological lensing surveys before the arrival of Euclid. J-PAS unprecedented spectral time domain information will enable a self-contained SN survey that, without the need for external spectroscopic follow-up, will detect, classify and measure $\sigma_z\sim 0.5\%$ redshifts for $\sim 4000$ SNeIa and $\sim 900$ core-collapse SNe. The key to the J-PAS potential is its innovative approach: a contiguous system of 54 filters with $145\AA$ width, placed $100\AA$ apart over a multi-degree FoV is a powerful "redshift machine", with the survey speed of a 4000 multiplexing low resolution spectrograph, but many times cheaper and much faster to build. The J-PAS camera is equivalent to a 4.7 sq.deg. "IFU" and it will produce a time-resolved, 3D image of the Northern Sky with a very wide range of Astrophysical applications in Galaxy Evolution, the nearby Universe and the study of resolved stellar populations.

We report the results obtained for the AGN Markarian 421 by model-fitting the data in the visibility plane, studing the proper motion of jet components, the light curve, and the spectral index of the jet features. We compare the radio data with optical light curves obtained at the Steward Observatory, considering also the optical polarization information. Mrk 421 has a bright nucleus and a one-sided jet extending towards the north-west for a few parsecs. The model-fits show that brightness distribution is well described using 6-7 circular Gaussian components, four of which are reliably identified at all epochs; all components are effectively stationary except for component D, at ~0.4 mas from the core, whose motion is however subluminal. Analysis of the light curve shows two different states, with the source being brighter and more variable in the first half of 2011 than in the second half. The highest flux density is reached in February. A comparison with the optical data reveals an increase of the V magnitude and of the fractional polarization simultaneous with the enhancement of the radio activity.

We report on observations of the nova-like cataclysmic variable AE Aqr performed by MAGIC. The observations were part of a quasi-simultaneous multi-wavelength campaign carried out between 2012 May and June covering the optical, UV, X-ray and gamma-ray ranges. MAGIC conducted the campaign and observed the source during 12 hours. The other instruments involved were KVA, Skinakas, and Vidojevica in the optic and Swift in the X-ray. We also used optical data from the AAVSO. The goals were to: monitor the variability of the source at different wavelengths, perform gamma-ray studies coincident with the highest states of the source at the other wavelengths, and confirm or rule out previous claims of detection of very-high-energy emission from this object. We report on a search for steady TeV emission during the whole observation, for variable TeV emission coincident with the highest optical and X-ray states and periodic TeV emission at the 33.08 s rotation period (30.23 mHz rotation frequency) of the white dwarf and its first harmonic (60.46 mHz rotation frequency). These are the first observations under good weather conditions performed by the present generation of IACTs for this object.

Luminous and Ultra-Luminous Infrared Galaxies (U/LIRGs) do also radiate copious amounts of radio emission, both thermal (free-free) and non-thermal (mainly synchrotron). This is very handy since, unlike optical and infra-red observations, radio is not obscured by the ubiquitous dust present in U/LIRGs, which allows a direct view of the ongoing activity in the hearts of those prolific star-forming galaxies. Here, I first justify the need for this high-angular resolution radio studies of local U/LIRGs, discuss the energy budget and the magnetic field, as well as IC and synchrotron losses in U/LIRGs, and present some selected results obtained by our team on high-angular resolution radio continuum studies of U/LIRGs. Among other results, I show the impressive discovery of an extremely prolific supernova factory in the central ~150 pc of the galaxy Arp 299-A (D=45 Mpc) and the monitoring of a large number of very compact radio sources in it, the detection and precise location of the long-sought AGN in Arp 299-A. A movie showing the appearance and disappearance of new SNe in Arp 299A can be found at http://www.iaa.es/~torres/research/arp299a.html. All those results demonstrate that very-high angular resolution studies of nearby U/LIRGs are of high relevance for the comprehension of both local and high-z starbursting galaxies.

High-resolution radio observations are ideal for constraining the value of physical parameters in the inner regions of active-galactic-nucleus jets and complement results on multiwavelength (MWL) observations. This study is part of a wider multifrequency campaign targeting the nearby TeV blazar Markarian 421 (z=0.031), with observations in the sub-mm (SMA), optical/IR (GASP), UV/X-ray (Swift, RXTE, MAXI), and gamma rays (Fermi-LAT, MAGIC, VERITAS). We investigate the jet's morphology and any proper motions, and the time evolution of physical parameters such as flux densities and spectral index. The aim of our wider multifrequency campaign is to try to shed light on questions such as the nature of the radiating particles, the connection between the radio and gamma-ray emission, the location of the emitting regions and the origin of the flux variability. We consider data obtained with the Very Long Baseline Array (VLBA) over twelve epochs (one observation per month from January to December 2011) at 15 GHz and 24 GHz. We investigate the inner jet structure on parsec scales through the study of model-fit components for each epoch. The structure of Mrk 421 is dominated by a compact (~0.13 mas) and bright component, with a one-sided jet detected out to ~10 mas. We identify 5-6 components in the jet that are consistent with being stationary during the 12-month period studied here. Measurements of the spectral index agree with those of other works: they are fairly flat in the core region and steepen along the jet length. Significant flux-density variations are detected for the core component. From our results, we draw an overall scenario in which we estimate a viewing angle 2{\deg} < theta < 5{\deg} and a different jet velocity for the radio and the high-energy emission regions, such that the respective Doppler factors are {\delta}r ~3 and {\delta}h.e. ~14.

We estimate the fraction of core-collapse supernovae (CCSNe) that remain undetected by optical SN searches due to obscuration by large amounts of dust in their host galaxies. This effect is especially important in luminous and ultraluminous infrared galaxies, which are locally rare but dominate the star formation at redshifts of z~1-2. We perform a detailed investigation of the SN activity in the nearby luminous infrared galaxy Arp 299 and estimate that up to 83% of the SNe in Arp 299 and in similar galaxies in the local Universe are missed by observations at optical wavelengths. For rest-frame optical surveys we find the fraction of SNe missed due to high dust extinction to increase from the average local value of ~19% to ~38% at z~1.2 and then stay roughly constant up to z~2. It is therefore crucial to take into account the effects of obscuration by dust when determining SN rates at high redshift and when predicting the number of CCSNe detectable by future high-z surveys such as LSST, JWST, and Euclid. For a sample of nearby CCSNe (distances 6-15 Mpc) detected during the last 12 yr, we find a lower limit for the local CCSN rate of 1.5 +0.4/-0.3 x 10^-4 yr^-1 Mpc^-3, consistent with that expected from the star formation rate. Even closer, at distances less than ~6 Mpc, we find a significant increase in the CCSN rate, indicating a local overdensity of star formation caused by a small number of galaxies that have each hosted multiple SNe.

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

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

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

The high star formation rates of luminous infrared galaxies (LIRGs) make them ideal places for core-collapse supernova (CCSN) searches. At radio frequencies, free from dust extinction, it is possible to detect compact components within the innermost LIRG nuclear regions, such as SNe and SN remnants, as well as AGN buried deep in the LIRG nuclei. We studied the LIRG IC883 aiming at: (i) investigating its (circum-)nuclear regions using the e-EVN at 5GHz, and e-MERLIN at 6.9GHz, complemented by archival VLBI data; (ii) detecting at radio frequencies the two recently reported circumnuclear SNe 2010cu and 2011hi, which were discovered by near-IR (NIR) adaptive optics observations of IC883; and (iii) further investigating the nature of SN2011hi at NIR by means of observations with Gemini-North. The circumnuclear regions traced by e-MERLIN at 6.9GHz have an extension of ~1kpc, and show a striking double-sided structure, which very likely corresponds to a warped rotating ring, in agreement with previous studies. Our e-EVN observations at 5GHz and complementary archival VLBI data at 5GHz and 8.4GHz, reveal the presence of various milliarcsec compact components in the nucleus of IC883. A single compact source, an AGN candidate, dominates the emission at both nuclear and circumnuclear scales, as imaged with the e-EVN and e-MERLIN, respectively. The other milliarcsec components are very suggestive of ongoing nuclear CCSN activity. Our e-EVN observations also resulted in upper limits to the radio luminosity of the two SNe in IC883 recently discovered at NIR. We refine the classification of SN2011hi as a Type IIP SN according to our latest Gemini-North epoch from 2012, in agreement with a low-luminosity radio SN nature. We estimate a CCSN rate lower limit of 1.1_{-0.6}^{+1.3} yr^{-1} for the entire galaxy, based on three nuclear radio SNe and the circumnuclear SNe 2010cu and 2011hi. (abridged)

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

We present a preliminary analysis of new high resolution radio observations of the nearby TeV blazar Markarian 421 (z=0.031). This study is part of an ambitious multifrequency campaign, with observations in sub-mm (SMA), optical/IR (GASP), UV/X-ray (Swift, RXTE, MAXI), and gamma rays (Fermi-LAT, MAGIC, VERITAS). In this manuscript we consider only data obtained with the Very Long Baseline Array (VLBA) at seven epochs (one observation per month from January to July 2011) at 15 and 23.8 GHz. We investigate the inner jet structure on parsec scales through the study of model-fit components for each epoch. We identified 5-6 components which are consistent with being stationary during the 6-month period reported here. The aim is to try to shed light on questions such as the nature of radiating particles, the connection between radio and gamma-ray emission, the location of emitting regions and the origin of the flux variability.

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

We review the effects of source size in interferometric observations and focus on the cases of very compact sources. If a source is extremely compact and/or weak (so it is not possible to detect signature of source structure in the visibilities) we describe a test of hypothesis that can be used to set a strong upper limit to the size of the source. We also estimate the minimum possible size of a source whose structure can still be detected by an interferometer (i.e., the maximum theoretical over-resolution power of an interferometer), which depends on the overall observing time, the compactness in the array distribution, and the sensitivity of the receivers. As a result, and depending on the observing frequency, the over-resolution power of forthcoming ultra-sensitive arrays, like the Square Kilometer Array (SKA), may allow us to study details of sources at angular scales down to a few micro-as.

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

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

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

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

The nuclear region of the Luminous Infra-red Galaxy Arp 299-A hosts a recent ($\simeq 10$ Myr), intense burst of massive star formation which is expected to lead to numerous core-collapse supernovae (CCSNe). Previous VLBI observations, carried out with the EVN at 5.0 GHz and with the VLBA at 2.3 and 8.4 GHz, resulted in the detection of a large number of compact, bright, non-thermal sources in a region $\lsim$150 pc in size. We aim at establishing the nature of all non-thermal, compact components in Arp 299-A, as well as estimating its core-collapse supernova rate. We use multi-epoch European VLBI Network (EVN) observations taken at 5.0 GHz to image with milliarcsecond resolution the compact radio sources in the nuclear region of Arp 299-A. We also use one single-epoch 5.0 GHz Multi-Element Radio Linked Interferometer Network (MERLIN) observation to image the extended emission in which the compact radio sources --traced by our EVN observations-- are embedded. Twenty-six compact sources are detected, 8 of them are new objects not previously detected. The properties of all detected objects are consistent with them being a mixed population of CCSNe and SNRs. We find clear evidence for at least two new CCSNe, implying a lower limit to the CCSN rate of $\nu_{\rm SN}\gsim$0.80 SN/yr indicating that the bulk of the current star formation in Arp 299-A is taking place in the innermost $\sim 150$ pc. Our MERLIN observations trace a region of diffuse, extended emission which is co-spatial to the region where all compact sources are found. From this diffuse, non-thermal radio emission we obtain an independent estimate for the core-collapse supernova rate, which is in the range $\nu_{\rm SN}=0.40$ - 0.65 SN/yr, roughly in agreement with previous estimates and our direct estimate of the CCSN rate from the compact radio emission.

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

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

We report the discovery of two consecutive supernovae (SNe), 2010cu and 2011hi, located at 0.37" (180 pc) and 0.79" (380 pc) projected distance respectively from the centre of the K-band nucleus of the luminous infrared galaxy IC 883. The SNe were discovered in an ongoing near-infrared K-band search for core-collapse SNe in such galaxies using the ALTAIR/NIRI adaptive optics system with laser guide star at the Gemini-North Telescope. These are thus the closest SNe yet discovered to a LIRG nucleus in optical or near-infrared wavelengths. The near-infrared light curves and colours of both SNe are consistent with core-collapse events. Both SNe seem to suffer from relatively low host galaxy extinction suggesting that regardless of their low projected galactocentric distances, they are not deeply buried in the nuclear regions of the host galaxy.

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

The nearby spiral galaxy M81 harbors in its core a Low-Luminosity AGN (LLAGN), and appears closely related to the more distant and powerful AGNs seen in quasars and radio galaxies. The intrinsic size of this object is unknown due to scattering, and it has shown a core-jet morphology with weak extended emission rotating with wavelength. The proximity of M\,81 (D=3.63 Mpc) allows a detailed investigation of its nucleus to be made. The nucleus is four orders of magnitude more luminous than the Galactic centre, and is therefore considered a link between SgrA* and the more powerful nuclei of radio galaxies and quasars. Our main goal was to determine the size of M81* at a shorter wavelength thus directly testing whether the frequency-size dependent law Theta propto nu^(-0.8) was still valid for wavelengths shorter than 1cm. In addition, we also aimed to confirm the rotation of the source as a function of frequency. We observed the continuum 7mm radio emission of M81* using the Very Long Baseline Array on Sep 13, 2002, using nearby calibrators to apply their interferometric observables to the target source, to increase the chances of detection. The source was detected on all baselines and hybrid mapping was possible. We present the first 7mm VLBI image of the core of M81*, which represents the highest resolution image ever of this source. Modeling the interferometric visibilities with two Gaussian functions sets constraints on the angular size of its core down to 38 microarcseconds, corresponding to a maximum (projected) linear size of 138AU, and shows extended emission towards the NE with a position angle of ~50deg. A fit of one Gaussian elliptical function yields a position angle of 28+/-8 degrees for its elongated, compact structure. Combining the 7mm size with earlier measurements at other frequencies we determine a frequency-size dependence of Theta propto nu^(-0.88+/-0.04).

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

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

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

We report on the VLBI detection of supernova SN2011dh at 22GHz using a subset of the EVN array. The observations took place 14 days after the discovery of the supernova, thus resulting in a VLBI image of the youngest radio-loud supernova ever. We provide revised coordinates for the supernova with milli-arcsecond precision, linked to the ICRF. The recovered flux density is a factor 2 below the EVLA flux density reported by other authors at the same frequency and epoch of our observations. This discrepancy could be due to extended emission detected with the EVLA or to calibration problems in the VLBI and/or EVLA observations.

(See the complete and formatted abstract in the paper). We report on VLBI monitoring of the low-luminosity AGN (LLAGN) in M81 at 1.7, 2.3, 5, and 8.4GHz. These observations are phase-referenced to the supernova SN1993J (located in the same galaxy) and cover from late 1993 to late 2005. The source consists at all frequencies of a slightly resolved core and a small jet extension towards the north-east direction (position angle of ~65 degrees) in agreement with previous publications. We find that the position of the intensity peak in the images at 8.4GHz is very stable in the galactic frame of M81 (proper motion upper limit about 0.010 mas per year). We confirm previous reports that the peaks at all frequencies are systematically shifted among them, possibly due to opacity effects in the jet as predicted by the standard relativistic jet model. We use this model to estimate the magnetic field in the jet and the mass of the central black hole. We obtain a black-hole mass of ~2.e+7 solar masses, comparable to estimates previously reported, but the magnetic fields obtained are thousands of times lower than previous estimates. We find that the positions of the cores at 1.7, 2.3, and 5GHz are less stable than that at 8.4GHz and evolve systematically, shifting southward at a rate of several tens of micro-arcsec per year. The evolution in the jet orientation seems to be related to changes in the inclination of the cores at all frequencies. These results can be interpreted as due to a precessing jet. The jet precession also seems to be related to a flare in the flux densities at 5.0 and 8.4GHz, which lasts ~4 years. A continued monitoring of the flux density and the jet structure evolution in this LLAGN will be necessary to further confirm our jet precession model.

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

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

Aims. The interaction of two isotropic stellar winds is studied in order to calculate the free-free emission from the wind collision region. The effects of the binary separation and the wind momentum ratio on the emission from the wind-wind interaction region are investigated. Methods. We developed a semi-analytical model for calculating the thermal emission from colliding stellar winds. Assuming radiative shocks for the compressed layer, which are expected in close binaries, we obtained the emission measure of the thin shell. Then, we computed the total optical depth along each line of sight to obtain the emission from the whole configuration. Results. Here, we present predictions of the free-free emission at radio frequencies from analytic, radiative shock models in colliding wind binaries. It is shown that the emission from the wind collision region mainly arises from the optically thick region of the compressed layer and scales as ~ D^{4/5}, where D is the binary separation. The predicted flux density from the wind collision region becomes more important as the frequency increases, showing higher spectral indices than the expected 0.6 value from the unshocked winds. We also investigate the emission from short-period WR+O systems calculated with our analytic formulation. In particular, we apply the model to the binary systems WR 98 and WR 113 and compare our results with the observations. Our theoretical results are in good agreement with the observed thermal spectra from these sources.

We present a study of the CCSN rate in nuclei A and B1 of the luminous infrared galaxy Arp299, based on 11 years of Very Large Array monitoring of their radio emission at 8.4 GHz. Significant variations in the nuclear radio flux density can be used to identify the CCSN activity in the absence of high-resolution very long baseline interferometry observations. In the case of the B1-nucleus, the small variations in its measured diffuse radio emission are below the fluxes expected from radio supernovae, thus making it well-suited to detect RSNe through flux density variability. In fact, we find strong evidence for at least three RSNe this way, which results in a lower limit for the CCSN rate of 0.28 +/- 0.16 per year. In the A-nucleus, we did not detect any significant variability and found a SN detection threshold luminosity which allows only the detection of the most luminous RSNe known. Our method is basically blind to normal CCSN explosions occurring within the A-nucleus, which result in too small variations in the nuclear flux density, remaining diluted by the strong diffuse emission of the nucleus itself. Additionally, we have attempted to find near-infrared counterparts for the earlier reported RSNe in the Arp299 nucleus A, by comparing NIR adaptive optics images from the Gemini-N telescope with contemporaneous observations from the European VLBI Network. However, we were not able to detect NIR counterparts for the reported radio SNe within the innermost regions of nucleus A. While our NIR observations were sensitive to typical CCSNe at 300 mas from the centre of the nucleus A, suffering from extinction up to A_v~15 mag, they were not sensitive to such highly obscured SNe within the innermost nuclear regions where most of the EVN sources were detected. (abridged)

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

The extragalactic background light (EBL) is of fundamental importance both for understanding the entire process of galaxy evolution and for gamma-ray astronomy. However, the overall spectrum of the EBL between 0.1 and 1000 microns has never been determined directly neither from observed luminosity functions (LFs), over a wide redshift range, nor from any multiwavelength observation of galaxy spectral energy distributions (SEDs). The evolving, overall spectrum of the EBL is derived here utilizing a novel method based on observations only. The changing fractions of quiescent galaxies, star-forming galaxies, starburst galaxies and active galactic nucleus (AGN) galaxies from redshift 0.2 to 1 are estimated, and two alternative extrapolations of SED types to higher redshifts are considered. This allows calculation of the evolving EBL. The EBL uncertainties in our modelling based directly on the data are quantified, and their consequences for attenuation of very-high-energy gamma-rays due to pair production on the EBL are discussed. It is concluded that the EBL seems well constrained from the UV to the mid-IR at an intensity level roughly matching galaxy count data. Independent efforts from IR and gamma-ray astronomy are needed in order to reduce the uncertainties in the far-IR.

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

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

Radio supernovae (RSNe) are weak and rare events. Their typical maximum radio luminosities are of the order of only $10^{27}$\,erg\,s$^{-1}$\,Hz$^{-1}$. There are, however, very few cases of relatively bright (and/or close) RSNe, from which the expansion of the shock and the radio light curves at several frequencies have been monitored covering several years. Applying the standard model of radio emission from supernovae, it is possible to relate the defining parameters of the modelled expansion curve to those of the modelled light curves in a simple algebraic way, by assuming an evolution law for the magnetic field and for the energy density of the population of synchrotron-emitting electrons. However, cooling mechanisms of the electrons may affect considerably this connection between light curves and expansion curve, and lead to wrong conclusions on the details of the electron acceleration and/or on the CSM radial density profile. In this paper, we study how electron cooling modifies the flux-density decay rate of RSNe for a set of plausible/realistic values of the magnetic field and for different expansion regimes. We use these results to estimate the magnetic fields of different RSNe observed to date and compare them to those obtained by assuming energy equipartition between particles and magnetic fields. For some of the best monitored RSNe, for which deceleration measurements, optically thin spectral index, and power-law time decay have been observed (SN\,1979C, SN\,1986J, SN\,1993J, and SN\,2008iz), we find self-consistent solutions for the index of the power-law circumstellar density profile ($s=2$ for all cases), the index of the power-law relativistic electron population (rather steep values, $ p = 2.3 - 3.0$) and the initial magnetic field (ranging from $\sim 20$ to $> 100$\,G).

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

The European VLBI Network (EVN) provides us with the necessary sensitivity and angular resolution to study the nuclear and circumnuclear regions in Luminous and Ultraluminous Infrared Galaxies. The high Star Formation Rates (SFR) inferred for these galaxies implies both the presence of a high number of massive stars and a dense surrounding medium. Therefore, bright radio SNe are expected to occur. With the aim of estimating the SFR in ULIRGs by means of Core Collapse supernova (CCSN) detections, we started an observing campaign with the EVN on a small sample of the brightest and farthest ULIRGs in the local Universe. We present here our results from three epochs of quasi-simultaneous observations with the EVN at 6 and 18 cm towards one of the objects in our sample: IRAS 23365+3604.

We present partial results from our monitoring of the nuclear region of the starburst galaxy IC 694 (=Arp 299-A) at radio wavelengths, aimed at discovering recently exploded CCSNe, as well as to determine their rate of explosion, which carries crucial information on star formation rates and starburst scenarios at work. Two epochs of eEVN observations at 5.0 GHz, taken in 2008, revealed the presence of a rich cluster of compact radio emitting sources in the central 150 pc of the nuclear starburst in Arp 299A. The large brightness temperatures observed for the compact sources indicate a non-thermal origin for the observed radio emission, implying that most, if not all, of those sources were young radio supernovae (RSNe) and supernova remnants (SNRs). More recently, contemporaneous EVN observations at 1.7 and 5.0 GHz taken in 2009 have allowed us to shed light on the compact radio emission of the parsec-scale structure in the nucleus of Arp 299-A. Namely, our EVN observations have shown that one of the compact VLBI sources, A1, previously detected at 5.0 GHz, has a flat spectrum between 1.7 and 5.0 GHz and is the brightest source at both frequencies. The morphology, radio luminosity, spectral index and ratio of radio-to-X-ray emission of the A1-A5 region allowed us to identify A1-A5 with long-sought AGN in Arp 299-A. This finding may suggest that both starburst and AGN are frequently associated phenomena in mergers. Finally, we also note that component A0, identified as a young RSN, exploded at the mere distance of two parsecs from the putative AGN in Arp 299-A, which makes this supernova one of the closest to a central supermassive black hole ever detected.

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

We have performed a complete re-calibration and re-analysis of all the available VLBI observations of supernova SN1993J, following an homogeneous and well-defined methodology. Observations of SN1993J at 69 epochs, spanning 13 years, were performed by two teams, which used different strategies and analysis tools. The results obtained by each group are similar, but their conclusions on the supernova expansion and the shape and evolution of the emitting region differ significantly. From our analysis of the combined set of observations, we have obtained an expansion curve with unprecedented time resolution and coverage. We find that the data from both teams are compatible when analyzed with the same methodology. One expansion index ($m_3 = 0.87 \pm 0.02$) is enough to model the expansion observed at 1.7\,GHz, while two expansion indices ($m_1 = 0.933\pm0.010$ and $m_2 = 0.796\pm0.005$), separated by a break time, $t_{br} = 390\pm30$ days, are needed to model the data, at frequencies higher than 1.7\,GHz, up to day 4000 after explosion. We thus confirm the wavelength dependence of the size of the emitting region reported by one of the groups. We also find that all sizes measured at epochs later than day 4000 after explosion are systematically smaller than our model predictions. We estimate the fractional shell width ($0.31 \pm 0.02$, average of all epochs and frequencies) and the level of opacity to the radio emission by the ejecta. We find evidence of a spectral-index radial gradient in the supernova shell, which is indicative of a frequency-dependent ejecta opacity. Finally, we study the distribution and evolution of the azimuthal anisotropies (hot spots) found around the radio shell during the expansion. These anisotropies have intensities of $\sim 20$% of the mean flux density of the shell, and appear to systematically evolve during the expansion.

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

The extragalactic background light (EBL) is of fundamental importance both for understanding the entire process of galaxy evolution and for gamma-ray astronomy, but the overall spectrum of the EBL between 0.1-1000 microns has never been determined directly from galaxy spectral energy distribution (SED) observations over a wide redshift range. The evolving, overall spectrum of the EBL is derived here utilizing a novel method based on observations only. This is achieved from the observed evolution of the rest-frame K-band galaxy luminosity function up to redshift 4 (Cirasuolo et al. 2010), combined with a determination of galaxy SED-type fractions. These are based on fitting SWIRE templates to a multiwavelength sample of about 6000 galaxies in the redshift range from 0.2 to 1 from the All-wavelength Extended Groth Strip International Survey (AEGIS). The changing fractions of quiescent galaxies, star-forming galaxies, starburst galaxies and AGN galaxies in that redshift range are estimated, and two alternative extrapolations of SED-types to higher redshifts are considered. This allows calculation of the evolution of the luminosity densities from the UV to the IR, the evolving star formation rate density of the universe, the evolving contribution to the bolometric EBL from the different galaxy populations including AGN galaxies and the buildup of the EBL. Our EBL calculations are compared with those from a semi-analytic model, from another observationally-based model and observational data. The EBL uncertainties in our modeling based directly on the data are quantified, and their consequences for attenuation of very high energy gamma-rays due to pair production on the EBL are discussed. It is concluded that the EBL is well constrained from the UV to the mid-IR, but independent efforts from infrared and gamma-ray astronomy are needed in order to reduce the uncertainties in the far-IR.

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

Context: The dusty nuclear regions of luminous infra-red galaxies (LIRGs) are heated by either an intense burst of massive star formation, an active galactic nucleus (AGN), or a combination of both. Disentangling the contribution of each of those putative dust-heating agents is a challenging task, and direct imaging of the innermost few pc can only be accomplished at radio wavelengths, using very high-angular resolution observations. Aims: We observed the nucleus A of the interacting starburst galaxy Arp 299, using very long baseline interferometry (VLBI) radio observations at 1.7 and 5.0 GHz. Our aim was to characterize the compact sources in the innermost few pc region of Arp 299-A, as well as to detect recently exploded core-collapse supernovae. Methods: We used the European VLBI Network (EVN) to image the 1.7 and 5.0 GHz compact radio emission of the parsec-scale structure in the nucleus of Arp 299-A with milliarcsecond resolution. Results: Our EVN observations show that one of the compact VLBI sources, A1, previously detected at 5.0 GHz, has a flat spectrum between 1.7 and 5.0 GHz and is the brightest source at both frequencies. Our 1.7 GHz EVN image shows also diffuse, low-surface brightness emission extending westwards from A1 and displays a prominent core-jet structure. Conclusions: The morphology, radio luminosity, spectral index and ratio of radio-to-X-ray emission of the A1-A5 region is consistent with a low-luminosity AGN (LLAGN), and rules out the possibility that it is a chain of young radio supernovae (RSNe) and supernova remnants (SNRs). We therefore conclude that A1-A5 is the long-sought AGN in Arp 299-A. This finding may suggest that both starburst and AGN are frequently associated phenomena in mergers.

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

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

We report on a simultaneous modelling of the expansion and radio light curves of SN1993J. We have developed a simulation code capable of generating synthetic expansion and radio light curves of supernovae by taking into consideration the evolution of the expanding shock, magnetic fields, and relativistic electrons, as well as the finite sensitivity of the interferometric arrays used in the observations. Our software successfully fits all the available radio data of SN 1993J with an standard emission model for supernovae extended with some physical considerations, as an evolution in the opacity of the ejecta material, a radial drop of the magnetic fields inside the radiating region, and a changing radial density profile of the circumstellar medium beyond day 3100 after explosion.

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

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

We present and analyse integral-field observations of six type-II QSOs with z=0.3-0.4, selected from the Sloan Digital Sky Survey (SDSS). Two of our sample are found to be surrounded by a nebula of warm ionized gas, with the largest nebula extending across 8" (40 kpc). Some regions of the extended nebulae show kinematics that are consistent with gravitational motion, while other regions show relatively perturbed kinematics: velocity shifts and line widths too large to be readily explained by gravitational motion. We propose that a ~20 kpc x20 kpc outflow is present in one of the galaxies. Possible mechanisms for triggering the outflow are discussed. In this object, we also find evidence for ionization both by shocks and the radiation field of the AGN.

We present long-slit spectroscopy, continuum and [OIII]5007 imaging data obtained with the Very Large Telescope and the Gran Telescopio Canarias of the type 2 quasar SDSS J0123+00 at z=0.399. The quasar lies in a complex, gas-rich environment. It appears to be physically connected by a tidal bridge to another galaxy at a projected distance of ~100 kpc, which suggests this is an interacting system. Ionized gas is detected to a distance of at least ~133 kpc from the nucleus. The nebula has a total extension of ~180 kpc. This is one of the largest ionized nebulae ever detected associated with an active galaxy. Based on the environmental properties, we propose that the origin of the nebula is tidal debris from a galactic encounter, which could as well be the triggering mechanism of the nuclear activity. SDSS J0123+00 demonstrates that giant, luminous ionized nebulae can exist associated with type 2 quasars of low radio luminosities, contrary to expectations based on type 1 quasar studies.

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

We present first pc-scale radio imaging of the radio-quiet candidate binary black hole system SDSS J1536+0441. The observations were carried out by the European VLBI Network at the frequency of 5 GHz and allowed to image SDSS J1536+0441 with a resolution of about 10 mas (50 pc). Two compact radio cores are detected at the position of the kpc-scale components VLA-A and VLA-B, proving the presence of two compact active nuclei with radio luminosity about 10^{40} erg/s, thus ruling out the possibility that the two radio sources are both powered by one 0.1 pc binary black hole. From a comparison with published 8.5 GHz flux densities we derived an estimate of the radio spectral index of the two pc-scale cores. Both cores have flat or inverted spectral index and, at least for the case of VLA-A, we can rule out the possibility that synchrotron self-absorption is responsible for the inverted radio spectrum. We suggest that thermal free-free emission from an X-ray heated disk wind may be powering the radio emission in VLA-A.

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

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

Context. The central kiloparsec of many local Luminous Infra-red Galaxies are known to host intense bursts of massive star formation, leading to numerous explosions of core-collapse supernovae (CCSNe). However, the dust-enshrouded regions where those supernovae explode hamper their detection at optical and near-infrared wavelengths. Aims. We investigate the nuclear region of the starbust galaxy IC 694 (=Arp 299-A) at radio wavelengths, aimed at discovering recently exploded CCSNe, as well as to determine their rate of explosion, which carries crucial information on star formation rates, the initial mass function and starburst scenarios at work. Methods. We use the electronic European VLBI Network to image with milliarcsecond resolution the 5.0 GHz compact radio emission of the innermost nuclear region of IC 694. Results. Our observations reveal the presence of a rich cluster of 26 compact radio emitting sources in the central 150 pc of the nuclear starburst in IC 694. The large brightness temperatures observed for the compact sources indicate a non-thermal origin for the observed radio emission, implying that most, if not all, of those sources are young radio supernovae (RSNe) and supernova remnants (SNRs). We find evidence for at least three relatively young, slowly-evolving, long-lasting RSNe (A0, A12, and A15) which appear to display unusual properties; suggesting that the conditions in the local circumstellar medium (CSM) play a significant role in determining the radio behaviour of expanding SNe. Their radio luminosities are typical of normal RSNe, which result from the explosion of Type IIP/b and Type IIL SNe. All these results yield support for a recent (less than 10-15 Myr) instantaneous starburst in the innermost regions of IC 694.

Core-collapse supernovae are the endproducts of massive stars, and yield radio events whose brightness depends on the intensity of the interaction experienced by the supernova ejecta with the circumstellar presupernova wind material. The fact that CCSNe are intrinsically radio supernovae --albeit with a huge range of different radio powers-- and hence unaffected by dust absorption, together with the high resolution and high sensitivity provided by current VLBI arrays, has been exploited to directly image the radio brightness structure of CCSNe in nearby (D <= 20 Mpc) galaxies. This has allowed to gain insight into the physics of both CCSNe and of the circumstellar medium (CSM) with which they interact. In addition, ultra-high-resolution, ultra-high-sensitivity radio observations of CCSNe in Luminous and Ultra-Luminous Infrared Galaxies (LIRGs and ULIRGs, respectively) in the local Universe, can be used to directly measure of the current CCSN rate and star formation rate. In this contribution, I give a brief overview of VLBI observations made of some CCSNe in nearby galaxies, and then present some of the most relevant results obtained with high-resolution radio observations of (U)LIRGs in the local Universe, aimed at directly detecting CCSNe via their radio emission, and thus determine their CCSN and star formation rates, independently of models. This is of particular relevance, in view of the fact that our estimates of star formation (and CCSN) rates in high-z starburst galaxies relies on standard relationships between far-infrared luminosity and star-formation rate. In particular, I will present recently obtained results with the e-EVN on the nuclear region of Arp 299-A.

We present quasi-simultaneous, multi-frequency VLA observations at 4.8, 8.4, and 22.5 GHz, of a sample of 13 Wolf Rayet (WR) stars, aimed at disentangling the nature of their radio emission and the possible detection of a non-thermal behavior in close binary systems. We detected 12 stars from our sample, for which we derived spectral information and estimated their mass loss rates. From our data, we identified four thermal sources (WR 89, 113, 138, and 141), and three sources with a composite spectrum (similar contribution of thermal and non-thermal emission; WR 8, 98, and 156). On the other hand, from the comparison with previous observations, we confirm the non-thermal spectrum of one (WR 105), and also found evidence of a composite spectrum for WR 79a, 98a, 104, and 133. Finally, we discuss the possible scenarios to explain the nature of the emission for the observed objects.

We have studied the growth of the shell-like radio structure of supernova SN1993J in M81 from September 1993 through October 2003 with very-long-baseline interferometry (VLBI) observations at the wavelengths of 3.6, 6, and 18cm. For this purpose, we have developed a method to accurately determine the outer radius (R) of any circularly symmetric compact radio structure like SN1993J. The source structure of SN1993J remains circularly symmetric (with deviations from circularity under 2%) over almost 4000 days. We characterize the decelerated expansion of SN 1993J through approximately day 1500 after explosion with an expansion parameter $m= 0.845\pm0.005$ ($R \propto t^{m}$). However, from that day onwards the expansion is different when observed at 6 and 18cm. Indeed, at 18cm, the expansion can be well characterized by the same $m$ as before day 1500, while at 6cm the expansion appears more decelerated, and is characterized by another expansion parameter, $m_{6}= 0.788\pm0.015$. Therefore, since about day 1500 on, the radio source size has been progressively smaller at 6cm than at 18cm. These findings are in stark contrast to previous reports by other authors on the details of the expansion. In our interpretation the supernova expands with a single expansion parameter, $m= 0.845\pm0.005$, and the 6cm results beyond day 1500 are due to physical effects, perhaps also coupled to instrumental limitations. Two physical effects may be involved: (a) a changing opacity of the ejecta to the 6cm radiation, and (b) a radial decrease of the magnetic field in the emitting region. (Long abstract cut. Please, read full abstract in manuscript).

We report on 1.6 GHz Very-Long-Baseline-Interferometry (VLBI) observations of supernova SN 1979C made on 18 November 2002. We derive a model-dependent supernova size. We also present a reanalysis of VLBI observations made by us on June 1999 and by other authors on February 2005. We conclude that, contrary to our earlier claim of strong deceleration in the expansion, SN 1979C has been undergoing almost-free expansion ($m = 0.91\pm0.09$; $R \propto t^m$) for over 25 years.

We present the results of an eight-year long monitoring of the radio emission from the Luminous Infrared Galaxy (LIRG) NGC 7469, using 8.4 GHz Very Large Array (VLA) observations at 0.3'' resolution. Our monitoring shows that the late time evolution of the radio supernova SN 2000ft follows a decline very similar to that displayed at earlier times of its optically thin phase. The late time radio emission of SN 2000ft is therefore still being powered by its interaction with the presupernova stellar wind, and not with the interstellar medium (ISM). Indeed, the ram pressure of the presupernova wind is \rho_w v_w^2 \approx 7.6E-9 dyn/cm^2, at a supernova age of approximately 2127 days, which is significantly larger than the expected pressure of the ISM around SN 2000ft. At this age, the SN shock has reached a distance r_{sh \approx 0.06 pc, and our observations are probing the interaction of the SN with dense material that was ejected by the presupernova star about 5820 years prior to its explosion. From our VLA monitoring, we estimate that the swept-up mass by the supernova shock after about six years of expansion is \approx 0.29 M_sun, assuming an average expansion speed of the supernova of 10000 km/s. We also searched for recently exploded core-collapse supernovae in our VLA images. Apart from SN 2000ft (S_\nu \approx 1760 microJy at its peak, corresponding to 1.1E28 erg/s/Hz, we found no evidence for any other radio supernova (RSN) more luminous than \approx 6.0E26 erg/s/Hz, which suggests that no other Type IIn SN has exploded since 2000 in the circumnuclear starburst of NGC 7469.

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

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

We present high resolution 240 and 607 MHz GMRT radio observations, complemented with 74 MHz archival VLA radio observations of the Ophiuchus cluster of galaxies, whose radio mini-halo has been recently detected at 1400 MHz. We also present archival Chandra and XMM-Newton data of the Ophiuchus cluster. Our observations do not show significant radio emission from the mini-halo, hence we present upper limits to the integrated, diffuse non-thermal radio emission of the core of the Ophiuchus cluster. The XMM-Newton observations can be well explained by a two-temperature thermal model with temperatures of ~=1.8 keV and ~=9.0 keV, respectively, which confirms previous results that suggest that the innermost central region of the Ophiuchus cluster is a cooling core. We also used the XMM-Newton data to set up an upper limit to the (non-thermal) X-ray emission from the cluster. The combination of available radio and X-ray data has strong implications for the currently proposed models of the spectral energy distribution (SED) from the Ophiuchus cluster. In particular, a synchrotron+IC model is in agreement with the currently available data, if the average magnetic field is in the range (0.02-0.3) microG. A pure WIMP annihilation scenario can in principle reproduce both radio and X-ray emission, but at the expense of postulating very large boost factors from dark matter substructures, jointly with extremely low values of the average magnetic field. Finally, a scenario where synchrotron and inverse Compton emission arise from PeV electron-positron pairs (via interactions with the CMB), can be ruled out, as it predicts a non-thermal soft X-ray emission that largely exceeds the thermal Bremsstrahlung measured by INTEGRAL.

We report on phase-referenced 23 GHz Very-Long-Baseline-Interferometry (VLBI) observations of the type IIb supernova SN 2008ax, made with the Very Long Baseline Array (VLBA) on 2 April 2008 (33 days after explosion). These observations resulted in a marginal detection of the supernova. The total flux density recovered from our VLBI image is 0.8$\pm$0.3 mJy (one standard deviation). As it appears, the structure may be interpreted as either a core-jet or a double source. However, the supernova structure could be somewhat confused with a possible close by noise peak. In such a case, the recovered flux density would decrease to 0.48$\pm$0.12 mJy, compatible with the flux densities measured with the VLA at epochs close in time to our VLBI observations. The lowest average expansion velocities derived from our observations are $(1.90 \pm 0.30) \times 10^5$ km s$^{-1}$ (case of a double source) and $(5.2 \pm 1.3) \times 10^4$ km s$^{-1}$ (taking the weaker source component as a spurious, close by, noise peak, which is the more likely interpretation). These velocities are 7.3 and 2 times higher, respectively, than the maximum ejecta velocity inferred from optical-line observations.

We show preliminary results from a sample of Luminous and Ultra-Luminous Infrared Galaxies (LIRGs and ULIRGs, respectively) in the local universe, obtained from observations using the Very Large Array (VLA), the Multi-Element Radio Link Interferometer Network (MERLIN), and the European VLBI Network (EVN). The main goal of our high-resolution, high-sensitivity radio observations is to unveil the dominant gas heating mechanism in the central regions of local (U)LIRGs. The main tracer of recent star-formation in (U)LIRGs is the explosion of core-collapse supernovae (CCSNe), which are the endproducts of the explosion of massive stars and yield bright radio events. Therefore, our observations will not only allow us to answer the question of the dominant heating mechanism in (U)LIRGs, but will yield also the CCSN rate and the star-formation rate (SFR) for the galaxies of the sample.

The radio source M81* at the core of the nearby spiral galaxy M81 is a low-luminosity active galactic nucleus. The close distance of 3.63Mpc allows its morphology to be studied in great detail. Here we present preliminary results from continuum 7 mm VLBI observations of its core, using phase-referencing techniques. These observations set constrains on the size of M81* at this frequency and enable us to test the frequency dependence on its physical properties.

We report the discovery of a confirmed supernova (SN) and a supernova-candidate in near-infrared images from the ALTAIR/NIRI adaptive optics system on the Gemini-North Telescope and NICMOS on the Hubble Space Telescope. The Gemini images were obtained as part of a near-infrared K-band search for highly-obscured SNe in the nuclear regions of luminous infrared galaxies. SN 2008cs apparent in the Gemini images is the first SN discovered using laser guide star adaptive optics. It is located at 1500 pc projected distance from the nucleus of the luminous infrared galaxy IRAS 17138-1017. The SN luminosity, JHK colors and light curve are consistent with a core-collapse event suffering from a very high host galaxy extinction of 15.7 +- 0.8 magnitudes in V-band which is to our knowledge the highest yet measured for a SN. The core-collapse nature of SN 2008cs is confirmed by its radio detection at 22.4 GHz using our Very Large Array observations 28 days after the SN discovery, indicating a prominent interaction of the SN ejecta with the circumstellar medium. An unconfirmed SN apparent in the NICMOS images from 2004 is located in the same galaxy at 660 pc projected distance from the nucleus and has a lower extinction.

We have compared the optical emission line ratios of type 2 quasars from Zakamska et al. with standard AGN photoionization model predictions, Seyfert 2s, HII galaxies, and narrow line FRII radio galaxies. Moderate to high ionization narrow line radio galaxies and Seyfert 2s are indistinguishable from type 2 quasars based on their optical line ratios. The standard AGN photoionization models, widely discussed for other type 2 AGNs, can reproduce successfully the loci and trends of type 2 quasars in some of the main diagnostic diagrams. These models are not exempt of problems and the discrepancies with the data are the same encountered for other type 2 AGNs. As for these, realistic models must take into account a range of cloud properties, as widely demonstrated in the literature. We also study the possible impact of stellar photoionization.

Aims: The TeV BL Lac object Markarian 501 is a complex, core dominated radio source, with a one sided, twisting jet on parsec scales. In the present work, we attempt to extend our understanding of the source physics to regions of the radio jet which have not been accessed before. Methods: We present new observations of Mrk 501 at 1.4 and 86 GHz. The 1.4 GHz data were obtained using the Very Large Array (VLA) and High Sensitivity Array (HSA) in November 2004, in full polarization, with a final r.m.s. noise of 25 microJy/beam in the HSA total intensity image; the 86 GHz observations were performed in October 2005 with the Global Millimeter VLBI Array (GMVA), providing an angular resolution as good as 110 x 40 microarcseconds. Results: The sensitivity and resolution provided by the HSA make it possible to detect the jet up to ~700 milliarcseconds (corresponding to a projected linear size of ~500 pc) from its base, while the superior resolution of the 86 GHz GMVA observations probes the innermost regions of the jet down to ~200 Schwarzschild radii. The brightness temperature at the jet base is in excess of 6e10 K. We find evidence of limb brightening on physical scales from <1 pc to ~40 pc. Polarization images and fits to the trend of jet width and brightness vs. distance from the core reveal a magnetic field parallel to the jet axis.

High-resolution radio observations of nearby starburst galaxies have shown that the distribution of their radio emission consists of a compact (<150 pc), high surface brightness, central radio source immersed in a low surface brightness circumnuclear halo. This radio structure is similar to that detected in bright Seyferts galaxies like NGC 7469 or Mrk 331, which display clear circumnuclear rings. While the compact, centrally located radio emission in these starbursts might be generated by a point-like source (AGN), or by the combined effect of multiple radio supernovae and supernova remnants (e.g., the evolved nuclear starburst in Arp~220), it seems well established that the circumnuclear regions of those objects host an ongoing burst of star-formation (e.g., NGC 7469; Colina et al. 2001, Alberdi et al. 2006). Therefore, high-resolution radio observations of Luminous Infra-Red Galaxies (LIRGs) in our local universe are a powerful tool to probe the dominant dust heating mechanism in their nuclear and circumnuclear regions. In this contribution, we show results obtained from VLA-A, MERLIN, and EVN (VLBI) radio observations of the galaxies NGC 7469 (D~70 Mpc) and IRAS 18293-3413 (D ~ 79 Mpc), where two extremely bright radio supernovae have been found. High-resolution studies of these and other LIRGs would allow us to determine the core-collapse supernova rate in them, as well as their star-formation rate.

We report on the first wide-field, high-precision astrometric analysis of the 13 extragalactic radio sources of the complete S5 polar cap sample at 15.4 GHz. We describe new algorithms developed to enable the use of differenced phase delays in wide-field astrometric observations and discuss the impact of using differenced phase delays on the precision of the wide-field astrometric analysis. From this global fit, we obtained estimates of the relative source positions with precisions ranging from 14 to 200 $\mu$as at 15.4 GHz, depending on the angular separation of the sources (from $\sim$1.6 to $\sim$20.8 degrees). These precisions are $\sim$10 times higher than the achievable precisions using the phase-reference mapping technique.

We report on 8.4GHz Very Long Baseline Interferometry (VLBI) observations of the type II-P supernova SN2004et in the spiral galaxy NGC 6946, made on 20 February 2005 (151 days after explosion). The Very Large Array (VLA) flux density was 1.23$\pm$0.07 mJy, corresponding to an isotropic luminosity at 8.4GHz of (4.45$\pm$0.3)$\times10^{25}$ erg s$^{-1}$ Hz$^{-1}$ and a brightness temperature of (1.3$\pm$0.3)$\times10^{8}$ K. We also provide an improved source position, accurate to about 0.5 mas in each coordinate. The VLBI image shows a clear asymmetry. From model fitting of the size of the radio emission, we estimate a minimum expansion velocity of 15,700$\pm$2,000 km s$^{-1}$. This velocity is more than twice the expected mean expansion velocity estimated from a synchrotron self-absorbed emission model, thus suggesting that synchrotron self-absorption is not relevant for this supernova. With the benefit of an optical spectrum obtained 12 days after explosion, we favor an emission model which consists of two hot spots on an underlying expanding shell of width comparable to that of SN 1993J.

We present new continuum VLA observations of the nearby Sy 1.5 galaxy NGC 5033, made at 4.9 and 8.4 GHz on 8 April 2003. Combined with VLA archival observations at 1.4 and 4.9 GHz made on 7 August 1993, 29 August 1999, and 31 October 1999, we sample the galaxy radio emission at scales ranging from the nuclear regions (<~ 100 pc) to the outer regions of the disk (~ 40 kpc). The high-resolution VLA images show a core-jet structure for the Sy 1.5 nucleus. While the core has a moderately steep non-thermal radio spectrum (S_\nu ~ \nu^\alpha; \alpha_{1.5}^{4.9} \approx -0.4), the inner kpc region shows a steeper spectrum (\alpha_{1.5}^{8.4} \approx -0.9). This latter spectrum is typical of galaxies where energy losses are high, indicating that the escape rate of cosmic ray electrons in NGC 5033 is low. The nucleus contributes little to the total 1.4 GHz radio power of NGC 5033 and, based on the radio to far-infrared (FIR) relation, it appears that the radio and far-infrared emission from NGC 5033 are dominated by a starburst that during the last 10 Myr produced stars at a rate of ~2.8 M_sun/yr, yielding a supernova (type Ib/c and II) rate of 0.045 #/yr. This supernova rate corresponds to about 1 SN event every 22 yr. Finally, from our deep 8.4 GHz VLA-D image, we suggest the existence of a radio spur in NGC 5033, which could have been due to a hot superbubble formed as a consequence of sequential supernova explosions occurring during the lifetime of a giant molecular cloud.

We report on first results of an ongoing effort to image a small sample of high-redshift radio galaxies (HzRGs) with milliarcsecond (mas) resolution, using very-long-baseline interferometry (VLBI) techniques. Here, we present 1.7 and 5.0 GHz VLBA observations of B3 J2330+3927, a radio galaxy at z=3.087. Those observations, combined with 8.4 GHz VLA-A observations, have helped us interpret the source radio morphology, and most of our results have already been published (Perez-Torres & De Breuck 2005). In particular, we pinpointed the core of the radio galaxy, and also detected both radio lobes, which have a very asymmetric flux density ratio, R>11. Contrary to what is seen in other radio galaxies, it is the radio lobe furthest from the nucleus which is the brighest. Almost all of the Ly-alpha emission is seen between the nucleus and the furthest radio lobe, which is also unlike all other radio galaxies. The values of radio lobe distance ratio, and flux density ratio, as well as the fraction of core emission make of B3 J2330+3927 an extremely asymmetric source, and challenges unification models that explain the differences between quasars and radio galaxies as due to orientation effects.

We report on 1.7 and 5.0 GHz observations of the z=3.087 radio galaxy B3 J2330+3927, using the Very Long Baseline Array (VLBA), and archival 1.4 and 8.4 GHz Very Large Array (VLA) data. Our VLBA data identify a compact, flat spectrum (\alpha_{1.7 GHz}^{5 GHz} = -0.2 +/- 0.1; S_\nu ~ \nu^\alpha) radio component as the core. The VLA images show that the fraction of core emission is very large (f_c \approx 0.5 at 8.4 GHz), and reveal a previously undetected, very faint counterjet, implying a radio lobe flux density ratio R >= 11 and a radio lobe distance ratio Q \approx 1.9. Those values are much more common in quasars than in radio galaxies, but the optical/near-IR spectra show a clear type II AGN for B3 J2330+3927, confirming that it is indeed a radio galaxy. Unlike all other radio galaxies, the bright Ly-\alpha emitting gas is located towards the furthest radio arm. We argue against environmental and relativistic beaming effects being the cause of the observed asymmetry, and suggest this source has intrinsically asymmetric radio jets. If this is the case, B3 J2330+3927 is the first example of such a source at high redshift, and seems to be difficult to reconcile with the unified model, which explains the differences between quasars and radio galaxies as being due to orientation effects.

We report on 1.6 and 5.0 GHz observations of the ultraluminous infrared galaxy (ULIRG) Mrk 273, using the European VLBI Network (EVN) and the Multi-Element Radio-Linked Interferometer Network (MERLIN). We also make use of published 1.4 GHz VLBA observations of Mrk 273 by Carilli & Taylor (2000). Our 5 GHz images have a maximum resolution of 5-10 mas, which corresponds to linear resolutions of 3.5-7 pc at the distance of Mrk 273, and are the most sensitive high-resolution radio observations yet made of this ULIRG. Component N1, often pinpointed as a possible AGN, displays a steep spectral index ($\alpha = 1.2 \pm 0.1; S_\nu \propto \nu^{-\alpha}$); hence it is very difficult to reconcile with N1 being an AGN, and rather suggests that the compact nonthermal radio emission is produced by an extremely high luminous individual radio supernova (RSN), or a combination of unresolved emission from nested supernova remnants (SNR), luminous RSNe, or both. Component N2 is partly resolved out into several compact radio sources --none of which clearly dominates-- and a region of extended emission about 30 pc in size. The integrated spectral index of this region is flat ($\alpha = 0.15 \pm 0.1$), which can be interpreted as due to a superposition of several unresolved components, e.g., RSNe or SNRs, whose radio emission peaks at different frequencies and is partially free-free absorbed. The overall extended radio emission from component N is typical of nonthermal, optically thin radio emission ($\alpha = 0.8 \pm 0.1$), and its 1.4 GHz luminosity ($L_{1.4 \rm GHz} = (2.2 \pm 0.1)\times 10^{23} $ WHz$^{-1}$) is consistent with being produced by relativistic electrons diffused away from supernova remnants in an outburst.

We report on 8.4 GHz VLBI observations of SN2001gd in the spiral galaxy NGC5033 made on 26 June 2002 and 8 April 2003. Our data nominally suggests a relatively strong deceleration for the expansion of SN2001gd, but we cannot dismiss the possibility of a free supernova expansion. From our VLBI observations on 8 April 2003, we inferred a minimum total energy in relativistic particles and magnetic fields in the supernova shell of E_min =(0.3-14) 10^{47} ergs, and a corresponding equipartition average magnetic field of B_min = (50--350) mG. We also present multiwavelength VLA measurements of SN2001gd, which are well fit by an optically thin, synchrotron spectrum, partially absorbed by thermal plasma. We obtain a supernova flux density of (1.02 +/- 0.05) mJy at the observing frequency of 8.4 GHz for the second epoch, which results in an isotropic radio luminosity of (6.0 +/- 0.3) * 10^{36} ergs between 1.4 and 43.3 GHz, at an adopted distance of 13.1 Mpc. Finally, we report on an XMM-Newton X-ray detection of SN2001gd on 18 December 2002. The supernova X-ray spectrum is consistent with optically thin emission from a soft component (associated with emission from the reverse shock) at a temperature around 1 keV. The observed flux corresponds to an isotropic X-ray luminosity of L_X = 1.4 +/- 0.4 * 10^{39} ergs/s in the (0.3-5) keV band. We suggest that both radio and X-ray observations of SN2001gd indicate that a circumstellar interaction similar to that displayed by SN1993J in M81 is taking place.

We give an overview of circumstellar interaction in young Type II supernovae, as seen through the eyes of very-long-baseline interferometry (VLBI) observations. The resolution attained by such observations (best than 1 mas) is a powerful tool to probe the interaction that takes place after a supernova goes off. The direct imaging of a supernova permits, in principle, to estimate the deceleration of its expansion, and to obtain information on the eject and circumstellar density profiles, as well as estimates of the magnetic field intensity and relativistic particle energy density in the supernova. Unfortunately, only a handful of radio supernovae are close and bright enough as to permit their study with VLBI. We present results from our high-resolution observations of the nearby Type II radio supernovae SN1986J and SN2001gd.

We report on the status of our S5 polar cap astrometry program. Since 1997 we have observed all the 13 radio sources of the complete S5 polar cap sample at the wavelengths of 3.6 cm, 2 cm and 7 mm. Images of the radio sources at 3.6 and 2 cm have already been published reporting morphological changes. Preliminary astrometric analyses have been carried out at three frequencies with precisions in the relative position determination ranging from 80 to 20 microarcseconds. We report also on the combination of our phase-delay global astrometry results with the microarcsecond-precise optical astrometry that will be provided by future space-based instruments.

We observed the thirteen extragalactic radio sources of the complete S5 polar cap sample at 15.4 GHz with the Very Long Baseline Array, on 27 July 1999 (1999.57) and 15 June 2000 (2000.46). We present the maps from those two epochs, along with maps obtained from observations of the 2 cm VLBA survey for some of the sources of the sample, making a total of 40 maps. We discuss the apparent morphological changes displayed by the radio sources between the observing epochs. Our VLBA observations correspond to the first two epochs at 15.4 GHz of a program to study the absolute kinematics of the radio source components of the members of the sample, by means of phase delay astrometry at 8.4 GHz, 15.4 GHz, and 43 GHz. Our 15.4 GHz VLBA imaging allowed us to disentangle the inner milliarcsecond structure of some of the sources, thus resolving components that appeared blended at 8.4 GHz. For most of the sources, we identified the brightest feature in each radio source with the core. These identifications are supported by the spectral index estimates for those brightest features, which are in general flat, or even inverted. Most of the sources display core-dominance in the overall emission. We find that three of the sources have their most inverted spectrum component shifted with respect to the origin in the map, which approximately coincides with the peak-of-brightness at both 15.4 GHz and 8.4 GHz.

The high resolution obtained through the use of VLBI gives an unique opportunity to directly observe the interaction of an expanding radio supernova with its surrounding medium. We present here results from our VLBI observations of the young supernovae SN 1979C, SN 1986J, and SN 2001gd.

We present the results of monitoring the radio emission from the Type IIb supernova SN 2001gd between 2002 February 8 and 2002 October 28. Most of the data were obtained using the Very Large Array at the five wavelengths of $\lambda \lambda$1.3 cm (22.4 GHz), 2.0 cm (14.9 GHz), 3.6 cm (8.44 GHz), 6.2 cm (4.86 GHz), and 21 cm (1.4 GHz). Observations were also made with Giant Meterwave Radio Telescope at $\lambda$21 cm (1.4 GHz). The object was discovered optically well after maximum light, making any determination of the early radio evolution difficult. However, subsequent observations indicate that the radio emission has evolved regularly in both time and frequency and is well described by the SN shock/circumstellar medium interaction model.

We present VLA radio continuum measurements of SN1993J in M81 at the frequencies of 0.32 (P-band), 1.3 and 1.7 (L-band), 4.9 (C-band), 8.5 (X-band), and 14.9 (U-band) GHz carried out on December 17 and 21, 2000, about 2820 days after the supernova explosion. We find that a power-law spectrum, free-free absorbed by an homogeneous, or clumpy, distribution of ionized gas yields the best fit to the radio data. A combined homogeneous-clumpy model is not favored, but neither totally excluded. This result contrasts with the modeling of the early ($t \lessim $230 days) radio emission from SN1993J, where a mixture of homogeneous and clumpy absorbers appeared to be necessary to adequately describe the behavior of the light curves. The radio spectrum of supernova SN1993J between 0.32 and 14.9 GHz is well characterized by $\alpha = -0.67 \pm 0.02 $ ($S_\nu \propto \nu^\alpha$), typical of an optically thin radio supernova. A fit to the radio spectra of SN1993J from $\sim$ 70 up to 2820 days shows that the observed spectral index of SN1993J has been slowly evolving since $t\sim$1000 days on, with the observed spectral index changing from $\alpha \approx -1$ to $\alpha=-0.67$. This spectral evolution seems to suggest that, in addition to the radiative (synchrotron) losses, adiabatic cooling and ionization (Coulomb) losses at the lowest frequencies might be contributing significantly to the integrated electron spectrum.

We report here on 5 GHz global very-long-baseline interferometry (VLBI) observations of SN 1986J, 16 yr after its explosion. We obtained a high-resolution image of the supernova, which shows a distorted shell of radio emission, indicative of a deformation of the shock front. The angular size of the shell is $\sim4.7 {\rm mas}$, corresponding to a linear size of $\sim6.8 \times 10^{17} {\rm cm}$ for a distance of 9.6 Mpc to NGC 891. The average speed of the shell has decreased from $\sim$7400 \kms in 1988.74 down to about $6300 {\rm km s^{-1}}$ in 1999.14, indicative of a mild deceleration in the expansion of SN 1986J. Assuming a standard density profile for the progenitor wind ($\rho_{\rm cs} \propto r^{-s}, s=2$), the swept-up mass by the shock front is $\sim$2.2\Ms. This large swept-up mass, coupled with the mild deceleration suffered by the supernova, suggests that the mass of the hydrogen-rich envelope ejected at explosion was $\ga12$ \Ms. Thus, the supernova progenitor must have kept intact most of its hydrogen-rich envelope by the time of explosion, which favours a single, massive star progenitor scenario. We find a flux density for SN 1986J of $\sim$7.2 mJy at the observing frequency of 5 GHz, which results in a radio luminosity of $\sim1.4 \times 10^{37} {\rm erg s^{-1}}$ for the frequency range $10^7$--$10^{10}$ Hz ($\alpha =-0.69; S_\nu \propto \nu^{\alpha})$. We detect four bright knots that delineate the shell structure, and an absolute minimum of emission, which we tentatively identify with the centre of the supernova explosion. If this is the case, SN 1986J has then suffered an asymmetric expansion. We suggest that this asymmetry is due to the collision of the supernova ejecta with an anisotropic, clumpy (or filamentary) medium.

SN1993J is to date the radio supernova whose evolution has been monitored in greatest detail and the one which holds best promise for a comprehensive theoretical-observational analysis. The shell-like radio structure of SN1993J has expanded in general accord with models of shock excited emission, showing almost circular symmetry for over 8 years, except for a bright feature at the south-eastern region of the shell that has been observed at every epoch. The spectrum of SN1993J has flattened from alpha =-1 to alpha =-0.67 (S_(\nu) propto nu**(alpha)). The decelerated expansion can be modeled well with a single slope but apparently better with two slopes. There are also intriguing hints of structure in the expansion curve. The results by the two VLBI groups carrying out this research show general agreement, but also some differences. A comparison of the optical and VLBI results about the details of the deceleration show some discrepancies.

We observed SN1979C in M100 on 4 June 1999, about twenty years after explosion, with a very sensitive four-antenna VLBI array at the wavelength of 18cm. The distance to M100 and the expansion velocities are such that the supernova cannot be fully resolved by our Earth-wide array. Model-dependent sizes for the source have been determined and compared with previous results. We conclude that the supernova shock was initially in free expansion for 6 +/- 2 yrs and then experienced a very strong deceleration. The onset of deceleration took place a few years before the abrupt trend change in the integrated radio flux density curves. We estimate the shocked swept-up mass to be about 1.6 solar masses, assuming a standard density profile for the CSM. Such a swept-up mass for SN1979C suggests a mass of the hydrogen-rich envelope ejected at explosion no larger than $M_{\rm env} \sim$ 0.9 solar masses. If SN1979C originated in a binary star the low value of M_env suggests that the companion of the progenitor star stripped off most of the hydrogen-rich envelope mass of the presupernova star prior to the explosion.