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

Among the four known transneptunian dwarf planets, Haumea is an exotic, very elongated, and fast rotating body. In contrast to the other dwarf planets, its size, shape, albedo, and density are not well constrained. Here we report results of a multi-chord stellar occultation, observed on 2017 January 21. Secondary events observed around the main body are consistent with the presence of a ring of opacity 0.5, width 70 km, and radius 2,287$_{-45}^{+75}$ km. The Centaur Chariklo was the first body other than a giant planet to show a ring system and the Centaur Chiron was later found to possess something similar to Chariklo's rings. Haumea is the first body outside the Centaur population with a ring. The ring is coplanar with both Haumea's equator and the orbit of its satellite Hi'iaka. Its radius places close to the 3:1 mean motion resonance with Haumea's spin period. The occultation by the main body provides an instantaneous elliptical limb with axes 1,704 $\pm$ 4 km x 1,138 $\pm$ 26 km. Combined with rotational light-curves, it constrains Haumea's 3D orientation and its triaxial shape, which is inconsistent with a homogeneous body in hydrostatic equilibrium. Haumea's largest axis is at least 2,322 $\pm$ 60 km, larger than thought before. This implies an upper limit of 1,885 $\pm$ 80 kg m$^{-3}$ for Haumea's density, smaller and less puzzling than previous estimations, and a geometric albedo of 0.51 $\pm$ 0.02, also smaller than previous estimations. No global N$_2$ or CH$_4$ atmosphere with pressures larger than 15 and 50 nbar (3-$\sigma$ limits), respectively, is detected.

On 28th January 2018, the large Trans-Neptunian Object (TNO) 2002TC302 occulted a m$_v= $15.3 star with ID 130957813463146112 in the Gaia DR2 stellar catalog. 12 positive occultation chords were obtained from Italy, France, Slovenia and Switzerland. Also, 4 negative detections were obtained near the north and south limbs. This represents the best observed stellar occultation by a TNO other than Pluto, in terms of the number of chords published thus far. From the 12 chords, an accurate elliptical fit to the instantaneous projection of the body, compatible with the near misses, can be obtained. The resulting ellipse has major and minor axes of 543 $\pm$ 18 km and 460 $\pm$ 11 km, respectively, with a position angle of 3 $\pm$ 1 degrees for the minor axis. This information, combined with rotational light curves obtained with the 1.5m telescope at Sierra Nevada Observatory and the 1.23m telescope at Calar Alto observatory, allows us to derive possible 3D shapes and density estimations for the body, based on hydrostatic equilibrium assumptions. The effective area equivalent diameter is $\sim$ 84 km smaller than the radiometrically derived diameter using thermal data from Herschel and Spitzer Space Telescopes. This might indicate the existence of an unresolved satellite of up to $\sim$ 300 km in diameter, to account for all the thermal flux, although the occultation and thermal diameters are compatible within their error bars given the considerable uncertainty of the thermal results. The existence of a potential satellite also appears to be consistent with other ground-based data presented here. From the effective occultation diameter combined with H$_V$ measurements we derive a geometric albedo of 0.147 $\pm$ 0.005, which would be somewhat smaller if 2002TC302 has a satellite. The best occultation light curves do not show any signs of ring features or any signatures of a global atmosphere.

Context. Red supergiant stars possess surface features and extended molecular atmospheres. Photospheric convection may be a crucial factor of the levitation of the outer atmospheric layers. However, the mechanism responsible is still poorly understood. Aims. We image the stellar surface of V602 Carinae (V602 Car) to constrain the morphology and contrast of the surface features and of the extended atmospheric layers. Methods. We observed V602 Car with the Very Large Telescope Interferometer (VLTI) PIONIER instrument (1.53-1.78 $\mathrm{\mu}$m) between May and July 2016, and April and July 2019 with different telescope configurations. We compared the image reconstructions with 81 temporal snapshots of 3D radiative-hydrodynamics (RHD) CO$^5$BOLD simulations in terms of contrast and morphology, using the Structural Similarity Index. Results. The interferometric data are compatible with an overall spherical disk of angular diameter 4.4$\pm$0.2 mas, and an extended molecular layer. In 2016, the reconstructed image reveals a bright arc-like feature toward the northern rim of the photospheric surface. In 2019, an arc-like feature is seen at a different orientation and a new peak of emission is detected on the opposite side. The contrasts of the reconstructed surface images are 11%$\pm$2% and 9%$\pm$2% for 2016 and 2019, respectively. The morphology and contrast of the two images are consistent with 3D RHD simulations, within our achieved spatial resolution and dynamic range. The extended molecular layer contributes 10--13% of the total flux with an angular diameter of 6--8 mas. It is present but not clearly visible in the reconstructed images because it is close to the limits of the achieved dynamic range. The presence of the molecular layer is not reproduced by the 3D RHD simulations. Conclusions...

Studies of fundamental parameters of very low-mass objects are indispensable to provide tests of stellar evolution models that are used to derive theoretical masses of brown dwarfs and planets. However, only objects with dynamically determined masses and precise photometry can effectively evaluate the predictions of stellar models. AB Dor C (0.090 solar masses) has become a prime benchmark for calibration of theoretical evolutionary models of low-mass young stars. One of the ambiguities remaining in AB Dor C is the possible binary nature of this star. We observed AB Dor C with the VLTI/AMBER instrument in low-resolution mode at the J, H and K bands. The interferometric observables at the K-band are compatible with a binary brown dwarf system with tentative components AB Dor Ca/Cb with a K-band flux ratio of 5$\pm$1% and a separation of 38$\pm$1 mas. This implies theoretical masses of 0.072$\pm$0.013 M$_{\rm \odot}$ and 0.013$\pm$0.001 M$_{\rm \odot}$ for each component, near the hydrogen-burning limit for AB Dor Ca, and near the deuterium-burning limit, straddling the boundary between brown dwarfs and giant planets, for AB Dor Cb. The possible binarity of AB Dor C alleviates the disagreement between observed magnitudes and theoretical mass-luminosity relationships.

The astrophysical agenda of the 21st century requires a very sharp view of celestial objects. High angular resolution studies are essential for fundamental studies of a broad variety of astrophysical phenomena ranging from relativistic physics of black holes, their gravitational and electromagnetic imprints, violent transient processes, including those producing detectable gravitational waves, birth and evolution of planetary systems. Over the past decades, radio astronomy made huge leap in achieving ground-breaking angular resolution measured in tens of microarcseconds (one tenth of nanoradian and better). Recently a global Event Horizon Telescope (EHT) collaboration obtained first direct images of the shadow of a super-massive black hole in the nucleus of the active galaxy M87. These observations were conducted at 230 GHz. The two first generation Space Very Long Baseline Interferometry (VLBI) missions, VSOP/HALCA led by the Japan Aerospace Exploration Agency (JAXA) and RadioAstron led by the Russia Roscosmos State Corporation and Russia Academy of Sciences, achieved the highest angular resolution at frequencies from 0.3 to 22 GHz in observations conducted in the period 1997 - 2019. The next step in advancing high angular resolution radio astronomy is in combining high frequency (millimeter and sub-millimeter wavelengths) and interferometric baselines exceeding the Earth diameter. The present THEZA White Paper describes a combination which would unify technology developments in giga-/tera-hertz instrumentation and space-borne radio astronomy. The current preprint version of the THEZA White Paper is slightly re-formatted and edited comparing to the official submitted version.

Context. Colliding winds in massive binaries are able to accelerate particles up to relativistic speeds as the result of the interaction between the winds of the different stellar components. HD 167971 exhibits this phenomenology which makes it a strong radio source. Aims. We aim at characterizing the morphology of the radio emission and its dependence on the orbital motion, traced independently by NIR-interferometry, of the spectroscopic binary and the tertiary component that conforms HD 167971. Methods. We analyze 2006 and 2016 very long baseline interferometric data at C and X bands. We complement our analysis with a geometrical model of the wind-wind collision region, and with an astrometric description of the system. Results. We confirm that the detected non-thermal radio emission is associated with the wind-wind collision region of the spectroscopic binary and the tertiary component in HD 167971. The wind-wind collision region changes orientation in agreement with the orbital motion of the tertiary around the spectroscopic binary. The total intensity also changes between the two observing epochs in a way inversely proportional to the separation between the SB and T, with a negative-steep spectral index typical of an optically thin synchrotron emission possibly steepened by an inverse Compton cooling effect. The wind-wind collision bow-shock shape and its position with respect to the stars indicates that the wind momentum from the spectroscopic binary is stronger than that of the tertiary. Finally, the astrometric solution derived for the stellar system and the wind-wind collision region is consistent with independent Gaia data.

With the purpose to investigate the radio emission of new ultracool objects, we carried out a targeted search in the recently discovered system VHS J125601.92$-$125723.9 (hereafter VHS 1256$-$1257); this system is composed by an equal-mass M7.5 binary and a L7 low-mass substellar object located at only 15.8\,pc. We observed in phase-reference mode the system VHS 1256$-$1257 with the Karl G. Jansky Very Large Array at $X$- and $L$- band and with the European VLBI Network at $L$-band in several epochs during 2015 and 2016. We discovered radio emission at $X$-band spatially coincident with the equal-mass M7.5 binary with a flux density of 60 $\mu$Jy. We determined a spectral index $\alpha = -1.1 \pm 0.3$ between 8 and 12 GHz, suggesting that non-thermal, optically-thin, synchrotron or gyrosynchrotron radiation is responsible for the observed radio emission. Interestingly, no signal is seen at $L$-band where we set a 3-$\sigma$ upper limit of 20 $\mu$Jy. This might be explained by strong variability of the binary or self-absorption at this frequency. By adopting the latter scenario and gyrosynchrotron radiation, we constrain the turnover frequency to be in the interval 5--8.5 GHz, from which we infer the presence of kG-intense magnetic fields in the M7.5 binary. Our data impose a 3-$\sigma$ upper bound to the radio flux density of the L7 object of 9 $\mu$Jy at 10\,GHz.

The star V766 Cen (=HR 5171A) was originally classified as a yellow hypergiant but lately found to more likely be a 27-36 Msun red supergiant (RSG). Recent observations indicated a close eclipsing companion in the contact or common-envelope phase. Here, we aim at imaging observations of V766 Cen to confirm the presence of the close companion. We used near-infrared H -band aperture synthesis imaging at three epochs in 2014, 2016, and 2017, employing the PIONIER instrument at the Very Large Telescope Interferometer (VLTI). The visibility data indicate a mean Rosseland angular diameter of 4.1+/-0.8 mas, corresponding to a radius of 1575+/-400 Rsun. The data show an extended shell (MOLsphere) of about 2.5 times the Rosseland diameter, which contributes about 30% of the H-band flux. The reconstructed images at the 2014 epoch show a complex elongated structure within the photospheric disk with a contrast of about 10%. The second and third epochs show qualitatively and quantitatively different structures with a single very bright and narrow feature and high contrasts of 20-30%. This feature is located toward the south-western limb of the photospheric stellar disk. We estimate an angular size of the feature of 1.7+/-0.3 mas, corresponding to a radius of 650+/-150 Rsun, and giving a radius ratio of 0.42+0.35/-0.10} compared to the primary stellar disk. We interpret the images at the 2016 and 2017 epochs as showing the close companion, or a common envelope toward the companion, in front of the primary. At the 2014 epoch, the close companion is behind the primary and not visible. Instead, the structure and contrast at the 2014 epoch are typical of a single RSG harboring giant photospheric convection cells. The companion is most likely a cool giant or supergiant star with a mass of 5+15/-3 Msun.

Context. Precise determination of stellar masses is necessary to test the validity of pre-main-sequence (PMS) stellar evolutionary models, whose predictions are in disagreement with measurements for masses below 1.2 Msun. To improve such a test, and based on our previous studies, we selected the AB Doradus moving group (AB Dor-MG) as the best-suited association on which to apply radio-based high-precision astrometric techniques to study binary systems. Aims. We seek to determine precise estimates of the masses of a set of stars belonging to the AB Dor-MG using radio and infrared observations. Methods. We observed in phase-reference mode with the Very Large Array (VLA) at 5 GHz and with the European VLBI Network (EVN) at 8.4 GHz the stars HD 160934, EK Dra, PW And, and LO Peg. We also observed some of these stars with the near-infrared CCD AstraLux camera at the Calar Alto observatory to complement the radio observations. Results. We determine model-independent dynamical masses of both components of the star HD 160934, A and c, which are 0.70+/-0.07 Msun and 0.45+/-0.04 Msun , respectively. We revised the orbital parameters of EK Dra and we determine a sum of the masses of the system of 1.38+/-0.08 Msun. We also explored the binarity of the stars LO Peg and PW And. Conclusions. We found observational evidence that PMS evolutionary models underpredict the mass of PMS stars by 10%-40%, as previously reported by other authors. We also inferred that the origin of the radio emission must be similar in all observed stars, that is, extreme magnetic activity of the stellar corona that triggers gyrosynchrotron emission from non-thermal, accelerated electrons.

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 add four warmer late-type supergiants to our previous spectro-interferometric studies of red giants and supergiants. V766 Cen (=HR 5171 A) is found to be a high-luminosity log(L/L_sun)=5.8+-0.4 source of Teff 4290+-760 K and radius 1490+-540 Rsun located close to both the Hayashi and Eddington limits; this source is consistent with a 40 Msun evolutionary track without rotation and current mass 27-36 Msun. It exhibits NaI in emission arising from a shell of radius 1.5 Rphot and a photocenter displacement of about 0.1 Rphot. V766 Cen shows strong extended molecular (CO) layers and a dusty circumstellar background component. This suggest an optically thick pseudo-photosphere at about 1.5 Rphot at the onset of the wind. V766 Cen is a red supergiant located close to the Hayashi limit instead of a yellow hypergiant already evolving back toward warmer Teff as previously discussed. The stars sigma Oph, BM Sco, and HD 206859 are found to have lower luminosities of about log(L/Lsun)=3.4-3.5 and Teff of 3900-5300 K, corresponding to 5-9 Msun tracks. They do not show extended molecular layers as observed for higher luminosity red supergiants of our sample. BM Sco shows an unusually strong contribution by an over-resolved circumstellar dust component. These stars are more likely high-mass red giants instead of red supergiants. This leaves us with an unsampled locus in the HR diagram corresponding to luminosities log(L/Lsun)~3.8-4.8 or masses 10-13 Msun, possibly corresponding to the mass region where stars explode as type II-P supernovae during the RSG stage. Our previously found relation of increasing strength of extended molecular layers with increasing luminosities is now confirmed to extend to double our previous luminosities and up to the Eddington limit. This might further point to steadily increasing radiative winds with increasing luminosity. [Abridged]

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.

Context. ABDoradus is the main system of the ABDoradus moving group. It is a quadruple system formed by two widely separated binaries of pre-main-sequence (PMS) stars: ABDor A/C and ABDor Ba/Bb. The pair ABDor A/C has been extensively studied and its dynamical masses have been determined with high precision, thus making of ABDor C a benchmark for calibrating PMS stellar models. If the orbit and dynamical masses of the pair ABDor Ba/Bb can be determined, they could not only play a similar role to that of ABDor C in calibrating PMS models, but would also help to better understand the dynamics of the whole ABDoradus system. Aims. We aim to determine the individual masses of the pair ABDor Ba/Bb using VLBI observations and archive infrared data, as part of a larger program directed to monitor binary systems in the ABDoradus moving group. Methods. We observed the system ABDor B between 2007 and 2013 with the Australian Long Baseline Array (LBA), at a frequency of 8.4 GHz in phase-reference mode. Results. We detected, for the first time, compact radio emission from both stars in the binary, ABDor Ba and ABDor Bb. This result allowed us to determine the orbital parameters of both the relative and absolute orbits and, consequently, their individual dynamical masses: 0.28+/-0.05 Msun and 0.25+/-0.05 Msun, respectively. Conclusions. Comparisons of the dynamical masses with the prediction of PMS evolutionary models show that the models underpredict the dynamical masses of the binary components Ba and Bb by 10-30% and 10-40%, respectively, although they all still agree at the 2-sigma level. Some of stellar models considered favour an age between 50 and 100 Myr for this system, meanwhile others predict older ages. We also discuss the evolutionary status of ABDor Ba/Bb in terms of an earlier double-double star scenario that might explain the strong radio emission detected in both components.

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.

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants

The main goal of this research is to determine the angular size and the atmospheric structures of cool giant stars and to compare them with hydrostatic stellar model atmospheres, to estimate the fundamental parameters, and to obtain a better understanding of the circumstellar environment. We conducted spectro-interferometric observations of epsilon Oct, beta Peg, NU Pav, and psi Peg in the near-infrared K band (2.13-2.47 microm), and gamma Hya (1.9-2.47 microm) with the VLTI/AMBER instrument at medium spectral resolution. To obtain the fundamental parameters, we compared our data with hydrostatic atmosphere models (PHOENIX). We estimated the Rosseland angular diameters of epsilon Oct, beta Peg, NU Pav, psi Peg, and gamma Hya. Together with distances and bolometric fluxes, we estimated radii, effective temperatures, and luminosities of our targets. In the beta Peg visibility, we observed a molecular layer of CO with a size similar to that modeled with PHOENIX. However, there is an additional slope in absorption starting around 2.3 microm. This slope is possibly due to a shell of water that is not modeled with PHOENIX. The visibility of psi Peg shows a low increase in the CO bands, compatible with the modeling of the PHOENIX model. The visibility data of the other sources show no increase in molecular bands. The spectra and visibilities predicted by the PHOENIX atmospheres agree with the spectra and the visibilities observed in our stars (except for beta Peg). This indicates that the opacity of the molecular bands is adequately included in the model, and the atmospheres of our targets have an extension similar to the modeled atmospheres. The atmosphere of beta Peg is more extended than that predicted by the model. The targets are located close to the red limits of the evolutionary tracks of the STAREVOL model.

Precise determination of dynamical masses of pre-main-sequence (PMS) stars is essential to calibrate stellar evolution models that are widely used to derive theoretical masses of young low-mass objects. Binary stars in young, nearby loose associations are particularly good candidates for this calibration since all members share a common age. Interestingly, some of these young binaries present a persistent and compact radio emission, which makes them excellent targets for astrometric VLBI studies. We aim to monitor the orbital motion of the binary system HD 160934, a member of the AB Doradus moving group. We observed HD 160934 with the Very Large Array and the European VLBI Network at 8.4 and 5 GHz, respectively. The orbital information derived from these observations was analyzed along with previously reported orbital measurements. We show that the two components of the binary, HD 160934 A and HD 160934 c, display compact radio emission at VLBI scales, providing precise information on the relative orbit. Revised orbital elements were estimated. Future VLBI monitoring of this pair should determine precise model-independent mass estimates for the A and c components, which will serve as calibration tests for PMS evolutionary models.

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.

The pre-main sequence (PMS) star ABDorA is the main component of the quadruple system ABDoradus. The precise determination of the mass and photometry of the close companion to ABDorA, ABDorC, has provided an important benchmark for calibration of theoretical evolutionary models of low-mass stars. The limiting factor to the precision of this calibration is the age of the system, as both the mass and luminosity of ABDorA and C are well monitored by other ongoing programs. In this paper we present VLTI/AMBER observations of ABDorA which provide a direct measurement of the size of this star, 0.96+/-0.06 Rsun. The latter estimate, combined with other fundamental parameters also measured for this star, allows a precise test of PMS evolutionary models using both H-R diagrams and mass-radius relationships. We have found that our radius measurement is larger than that predicted by the models, which we interpret as an evidence of the oversizing produced by the strong magnetic activity of ABDorA. Considering, at least partially, this magnetic effect, theoretical isochrones have been used to derive constraints to the age of ABDorA, favouring an age about 40-50 Myr for this system. Older ages are not completely excluded by our data.

We report on K-band VLTI/AMBER observations at medium spectral resolution ($\sim$1500) of RS Capricorni, an M6/M7III semi-regular AGB star. From the spectrally-dispersed visibilities, we measure the star diameter as a function of observing wavelength from 2.13 to 2.47 microns. We derive a Rosseland angular diameter of $7.95 \pm 0.07$ mas, which corresponds to an effective temperature of $3160 \pm 160$ K. We detect size variations of around 10% in the CO band heads, indicating strong opacity effects of CO in the stellar photosphere. We also detect a linear increase of the size as a function of wavelength, beginning at 2.29 microns. Models of the stellar atmosphere, based on the mass of the star as estimated from stellar-evolution models, predict CO-size effects about half of those observed, and cannot reproduce the linear size increase with wavelength, redward of 2.29 microns. We are able to model this linear size increase with the addition of an extended water-vapor envelope around the star. However, we are not able to fit the data in the CO bandheads. Either the mass of the star is overestimated by the stellar-evolution models and/or there is an additional extended CO envelope in the outer part of the atmosphere. In any case, neither the water-vapor envelope, nor the CO envelope, can be explained using the current models.

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.

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: Phase referencing is a standard calibration technique in radio interferometry, particularly suited for the detection of weak sources close to the sensitivity limits of the interferometers. However, effects from a changing atmosphere and inaccuracies in the correlator model may affect the phase-referenced images, leading to wrong estimates of source flux densities and positions. A systematic observational study of signal decoherence in phase referencing, and its effects in the image plane, has not been performed yet. Aims: We systematically studied how the signal coherence in Very-Long-Baseline-Interferometry (VLBI) observations is affected by a phase-reference calibration at different frequencies and for different calibrator-to-target separations. The results obtained should be of interest for a correct interpretation of many phase-referenced observations with VLBI. Methods: We observed a set of 13 strong sources (the S5 polar cap sample) at 8.4 and 15 GHz in phase-reference mode, with 32 different calibrator/target combinations spanning angular separations between 1.5 and 20.5 degrees. We obtained phase-referenced images and studied how the dynamic range and peak flux density depend on observing frequency and source separation. Results: We obtained dynamic ranges and peak flux densities of the phase-referenced images as a function of frequency and separation from the calibrator. We compared our results with models and phenomenological equations previously reported. Conclusions: The dynamic range of the phase-referenced images is strongly limited by the atmosphere at all frequencies and for all source separations. The limiting dynamic range is inversely proportional to the sine of the calibrator-to-target separation. We also find that the peak flux densities, relative to those obtained with the self-calibrated images, decrease with source separation.

Phase referencing is a standard calibration procedure in radio interferometry. It allows to detect weak sources by using quasi-simultaneous observations of closeby sources acting as calibrators. Therefore, it is assumed that, for each antenna, the optical paths of the signals from both sources are similar. However, atmospheric turbulence may introduce strong differences in the optical paths of the signals and affect, or even waste, phase referencing for cases of relatively large calibrator-to-target separations and/or bad weather. The situation is similar in wide-field observations, since the random deformations of the images, mostly caused by atmospheric turbulence, have essentially the same origin as the random astrometric variations of phase-referenced sources with respect to the phase center of their calibrators. In this paper, we present the results of a Monte Carlo study of the astrometric precision and sensitivity of an interferometric array (a realization of the Square Kilometre Array, SKA) in phase-referenced and wide-field observations. These simulations can be extrapolated to other arrays by applying the corresponding corrections. We consider several effects from the turbulent atmosphere (i.e., ionosphere and wet component of the troposphere) and also from the antenna receivers. We study the changes in dynamic range and astrometric precision as a function of observing frequency, source separation, and strength of the turbulence. We find that, for frequencies between 1 and 10 GHz, it is possible to obtain images with high fidelity, although the atmosphere strongly limits the sensitivity of the instrument compared to the case with no atmosphere. Outside this frequency window, the dynamic range of the images and the accuracy of the source positions decrease. [...] (Incomplete abstract. Please read manuscript.)

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 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 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 an astrometric analysis of the binary systems ABDorA /ABDorC and ABDorBa / ABDorBb. These two systems of well-known late-type stars are gravitationally associated and they constitute the quadruple ABDoradus system. From the astrometric data available at different wavelengths, we report: (i) a determination of the orbit of ABDorC, the very low mass companion to ABDorA, which confirms the mass estimate of 0.090Msun reported in previous works; (ii) a measurement of the parallax of ABDorBa, which unambiguously confirms the long-suspected physical association between this star and ABDorA; and (iii) evidence of orbital motion of ABDorBa around ABDorA, which places an upper bound of 0.4Msun on the mass of the pair ABDorBa / ABDorBb (50% probability). Further astrometric monitoring of the system at all possible wavelengths would determine with extraordinary precision the dynamical mass of its four components.

We expand upon the results of Close et al. 2005 regarding the young, low-mass object AB Dor C and its role as a calibration point for theoretical tracks. We present an improved spectral reduction and a new orbital solution with two additional epochs. Our improved analysis confirms our spectral type of M8 (+/- 1) and mass of 0.090+/-0.003 solar masses for AB Dor C. Comparing the results for AB Dor C with other young, low-mass objects with dynamical masses we find a general trend where current evolutionary models tend to over-predict the temperature (or under-predict the mass) for low mass stars and brown dwarfs. Given our precision, there is a ~99% chance that the mass of AB Dor C is underestimated by the DUSTY tracks in the HR diagram.

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

Nearby M-dwarfs are best suited for searches of low mass companions. VLBI phase-referencing observations with sensitive telescopes are able to detect radio star flux-densities of tenths of mJy as well as to position the star on the sky with submilliarcsecond precision. We have initiated a long-term observational program, using EVN telescopes in combination with NASA DSN dishes, to revisit the kinematics of nearby, single M dwarfs. The precision of the astrometry allows us to search for possible companions with masses down to 1 Jupiter mass. In this contribution we report preliminary results of the first observation epochs, in which we could detect some of the radio stars included in our program.

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.

We present multi-epoch 8.4 and 43 GHz Very Long Baseline Array images of the BL Lac object 0735+178. The images confirm the presence of a twisted jet with two sharp apparent bends of 90$^{\circ}$ within two milliarcseconds of the core, resembling a helix in projection. The observed twisted geometry could be the result of precession of the jet inlet, but is more likely produced by pressure gradients in the external medium through which the jet propagates. Quasi-stationary components are observed at the locations of the 90$^{\circ}$ bends, possibly produced by differential Doppler boosting. Identification of components across epochs, since the earliest VLBI observations of this source in 1979.2, proves difficult due to the sometimes large time gaps between observations. One possible identification suggests the existence of superluminal components following non--ballistic trajectories with velocities up to $11.6\pm 0.6 h_{65}^{-1} c$. However, in images obtained after mid-1995, components show a remarkable tendency to cluster near several jet positions, suggesting a different scenario in which components have remained nearly stationary in time at least since mid-1995. Comparison with the earlier published data, covering more than 19 years of observations, suggests a striking qualitative change in the jet trajectory sometime between mid-1992 and mid-1995, with the twisted jet structure with stationary components becoming apparent only at the later epochs. This would require a re-evaluation of the physical parameters estimated for 0735+178, such as the observing viewing angle, the plasma bulk Lorentz factor, and those deduced from these.

We observed the thirteen extragalactic radio sources of the S5 polar cap sample at 8.4 GHz with the Very Long Baseline Array, on 1997.93 and 1999.41. We present the maps from those two epochs and briefly discuss the morphological changes experimented by some of the radio sources in the 1.4 yr elapsed. These results correspond to the first two epochs at 8.4 GHz of a program directed 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, 15 and 43 GHz.

We present 5 GHz space-VLBI observations of the quasar pair B1342+662/B1342+663 that demonstrate the feasibility of phase-reference techniques using an antenna in space. The space-based data were taken by the satellite HALCA, of the space-VLBI mission VSOP. From residual (referenced) phases we derive an upper bound of 10 meters to the uncertainty of the spacecraft orbit reconstruction. An analysis of the phase-reference maps of the sources additionally suggests that the above mentioned uncertaintyis likely not larger than 3 meters. With errors of this magnitude, HALCA is a useful tool for astrometric studies of close pairs of radio sources.

We present a new centimeter polarimetric VLBI image of the BL Lac object 0735+178. This source exhibits one of the most pronounced curvatures observed in jets of AGNs, with two sharp apparent bends of 90 degrees within the inner 2 milliarcseconds from the core. Through the analysis of the data gathered over the past decades we study whether this curvature is produced by precession of the jet with ejection of ballistic components, or a precession in such a way that components' velocity vector are always parallel to the jet axis. These possibilities are also studied by comparison with 3D hydrodynamic relativistic simulations of precessing jets.

We present wide-field images of the quadruple gravitational lenses B1422+231 and MGJ0414+0534 obtained from global Very Long Baseline Interferometry (VLBI) observations at 8.4 GHz on 23 November 1997. We present also a lens model for MGJ0414+0534, which reproduces the core positions and flux densities of the VLBI images, combining a singular isothermal ellipsoid with external shear, and a singular isothermal sphere to represent, respectively, the main lens galaxy and its neighbor, a faint galaxy near one of the images.

The use of the phase-delay improves substantially the accuracy obtained in VLBI astrometry with respect to the group-delay observable. Recently, Ros et al. (1999) have extended the related phase-connection technique to a triangle of radio sources with relative separations up to 6.8deg (the S5 sources BL1803+784/ QSO1928+738/ BL2007+777). This technique has also been extended for separations up to 15deg in the studies of the pair of S5 radio sources QSO1150+812/ BL1803+784 (Perez-Torres et al. 2000). We are carrying out a long-term astrometric programme at 8.4, 15, and 43 GHz to determine the absolute kinematics of radio source components in the 13 members of the complete S5 polar cap sample. For each epoch, all the radio sources are observed at the same frequency for 24 hours in total. Bootstrapping techniques are used to phase-connect the data jointly for the 13 objects throughout the observations. An accurate registration of the maps of the radio sources at different frequencies will allow a study of jet components with unprecedented precision and provide spectral information. This observing scheme could be extended in the future to other regions of the sky and eventually lead to global phase-delay astrometry.

We observed the quadruple gravitationally lensed image of MGJ0414+0534 on 23 November 1997 with a global VLBI array at 8.4 GHz. We report wide-field imaging results of its four components at submilliarcsecond resolution, displaying complex core-like and jet-like extended structures. A simple model combining a singular isothermal ellipsoid to represent the main lens galaxy, external shear, and a singular isothermal sphere to represent an additional, nearby object accounts well for the core positions and flux densities of the VLBI images. This model predicts delays between the different lensed images of several weeks.

On 18 November 1993, we observed the QSO 1150+812 and the BL Lac object 1803+784, nearly 15\deg apart on the sky, with a VLBI array, simultaneously recording at 8.4 and 2.3 GHz. Using difference astrometry techniques, we determined the coordinates of 1803+784 relative to those of 1150+812 to be \Delta\alpha = 6h 7m 33.18469 arcsec +- 0.00020 arcsec \Delta\delta = -2deg 30' 25''.13557 +- 0''.00075 These standard errors contain estimated contributions from the propagation medium, the effects of source opacity, and the possible misidentification of the reference point in some of the maps. We obtained comparably accurate and consistent relative positions of the radio sources with GPS-based ionosphere corrections, thus demonstrating that dual-frequency observations are not required for state-of-the-art accuracy in VLBI astrometry. These results also demonstrate the feasibility of using single-frequency measurements for sources separated by 15$^\circ$ on the sky, and open the avenue for the application of this technique on a full-sky scale.

We report on the present status of our S5 polar cap phase-connected astrometry program. We observe 13 radio sources in the northernmost 20deg of the sky at the wavelengths of 3.6cm and 2cm, and we plan to extend the program to 0.7cm. We phase-connect jointly all our data successfully. We image the radio sources and some of them show morphological changes, in which astrometric registration is needed to determine the kinematics of the source components. We aim at unprecedented astrometric accuracy and at a check of the jet standard model at the 5-10 microarcsec/yr level.

From 43GHz VLBA observations of the pair of radio sources 1928+738 and 2007+777 we have demonstrated the feasibility of precision phase-delay differential astrometric techniques at millimeter wavelengths. For a pair of sources with 5 degrees separation and high antenna elevations, we have shown that present astrometric models and millimeter arrays are advanced enough to model the differential phase-delay to within 2 picoseconds, less than one tenth of a phase-cycle at 43GHz. The root-mean-square of the differential phase-delay residuals is dominated by the fluctuations of the atmospheric water vapor. We have determined the relative position of the observed sources with a precision twofold better than previous determinations at lower frequencies and, more importantly, largely free from ambiguous definitions of the reference point on the structure of the radio sources. Our result makes 43GHz VLBI phase-delay differential astrometry an ideal tool to study the absolute kinematics of the highly variable structures of regions near the core of extragalactic radio sources.

Global Positioning System (GPS) satellites broadcast at frequencies of 1,575.42 MHz (L1) and of 1,227.60 MHz (L2). The dispersive property of the ionosphere can be used to combine independent measurements at the two frequencies to estimate the total electron content (TEC) between a GPS receiver site and a broadcasting satellite. Such measurements, made at sites near to Very Long Baseline Interferometry (VLBI) sites, can be used to estimate the ionospheric contribution to VLBI observables. For our 1991.9 astrometric VLBI experiment in which we obtained group-delay observations in the 8.4 and 2.3 GHz bands simultaneously, we found that the GPS and VLBI determinations of the ionosphere delays agreed with root-mean-square differences below 0.15 ns for intercontinental baselines and 0.10 ns for continental ones. We also successfully applied the GPS-based procedure to reduce the ionospheric effect in phase delays used for high precision differenced astrometry at 8.4 GHz for this same experiment.

We determined the separations of the radio sources in the triangle formed by the BL Lac objects 1803+784 and 2007+777, and the QSO 1928+738 from intercontinental interferometric observations carried out in November 1991 at the frequencies of 2.3 and 8.4 GHz simultaneously. We successfully connected differenced phase delays over 7^o on the sky at 8.4 GHz at an epoch of maximum solar activity. The effects of the ionosphere on these VLBI data were mostly removed by estimates of the total electron content from observations of GPS satellites. The comparison of the estimated separation of QSO 1928+738 and BL 2007+777 with previous such estimates obtained from data at different epochs leads us to a particular alignment of the maps of QSO 1928+738 at different epochs relative to those of BL 2007+777, although with significant uncertainty. For this alignment, the jet components of QSO 1928+738 show a mean proper motion of 0.32+/-0.10 mas/yr and also suggest an identification for the position of the core of this radio source.

The VLBI technique can image compact radio sources with resolutions of the order of the milliarcsecond and can determine astrometrically relative positions with submilliarcsecond precisions. We determined the angular separations among the radio sources in the triangle 1803+784/1928+738/2007+777 with submilliarcsecond accuracy from a weighted least squares analysis of the differential phases. The comparison of the estimated separation of the pair 1928+738/2007+777 with previous such estimates obtained from data at different epochs allowed us to align the maps of 1928+738 at diferent epochs relative to those of 2007+777, showing that the jet components of 1928+738 display proper motions, and suggesting an identification for the position of the core of this radio source. We have also observed astrometrically the complete S5 polar cap sample, that among the 13 sources within 20 degrees to the celestial North Pole includes the above 3. These observations are carried out in the framework of a long-term astrometric program to determine the absolute kinematics of radio sources in the S5 complete sample.

We observed the pair of radio sources 1150+812 and 1803+784 in November 1993 with a VLBI array, simultaneously recording at 8.4 and 2.3 GHz. We determined the angular separation between the two sources with submilliarcsecond accuracy by using differential techniques. This result demonstrates the feasibility of high precision differential astrometry for radio sources separated in the sky by almost 15 degrees, and opens the avenue to its application to larger samples of radio sources.

A rarity among supernova, SN 1993J in M81 can be studied with high spatial resolution. Its radio power and distance permit VLBI observations to monitor the expansion of its angular structure. This radio structure was previously revealed to be shell-like and to be undergoing a self-similar expansion at a constant rate. From VLBI observations at the wavelengths of 3.6 and 6 cm in the period 6 to 42 months after explosion, we have discovered that the expansion is decelerating. Our measurement of this deceleration yields estimates of the density profiles of the supernova ejecta and circumstellar material in standard supernova explosion models.