We report Atacama-Large-millimeter-Array (ALMA) Band 6 full-polarization observations of the lensed blazar PKS1830-211 during its record-breaking radio and gamma-ray flare in the spring of 2019. The observations were taken close to the peak of the gamma activity and show a clear difference in polarization state between the two time-delayed images. The leading image has a fractional polarization about three times lower than the trailing image, implying that significant depolarization occurred during the flare. In addition, we observe clear intra-hour variability of the polarization properties between the two lensed images, with a quasi-linear increase of the differential electric-vector position angle at a rate of about two degrees per hour, associated with changes of the relative fractional polarization of ~10%. This variability, combined with the lower polarization close to the peak of gamma activity, is in agreement with models of magnetic turbulence to explain polarization variability in blazar jets. Finally, the comparison of results from the full and differential polarization analysis confirms that the differential polarization technique (Marti-Vidal et al. 2016) can provide useful information on the polarization state of sources like gravitationally lensed radio-loud quasars.

Deuterium and lithium are light elements of high cosmological and astrophysical importance. In this work we report the first detection of deuterated molecules and a search for lithium hydride, 7LiH, at redshift z=0.89 in the spiral galaxy intercepting the line of sight to the quasar PKS1830-211. We used ALMA to observe several submillimeter lines of ND, NH2D, and HDO, and their related isotopomers NH2, NH3, and H2^18O, in absorption against the southwest image of the quasar, allowing us to derive XD/XH abundance ratios. The absorption spectra mainly consist of two distinct narrow velocity components for which we find remarkable differences. One velocity component shows XD/XH abundances that is about 10 times larger than the primordial elemental D/H ratio, and no variability of the absorption profile during the time span of our observations. [...] The second component has XD/XH abundances that are 100 times larger than the primordial D/H ratio, a deepening of the absorption by a factor of two within a few months, and a rich chemical composition, with relative enhancements of N2H+, CH3OH, SO2, and complex organic molecules. We therefore speculate that this component is associated with the analog of a Galactic dark cloud, while the first component is likely more diffuse. Our search for the 7LiH (1--0) line was unsuccessful and we derive an upper limit 7LiH/H2 = 4 x 10^-13 (3sigma) in the z=0.89 absorber toward PKS1830-211. Besides, with ALMA archival data, we could not confirm the previous tentative detections of this line in the z=0.68 absorber toward B0218+357; we derive an upper limit 7LiH/H2 = 5 x 10^-11 (3sigma), although this is less constraining than our limit toward PKS1830-211. We conclude that, as in the Milky Way, only a tiny fraction of lithium nuclei are possibly bound in LiH in these absorbers at intermediate redshift.

In April 2019, the Event Horizon Telescope (EHT) collaboration revealed the first image of the candidate super-massive black hole (SMBH) at the centre of the giant elliptical galaxy Messier 87 (M87). This event-horizon-scale image shows a ring of glowing plasma with a dark patch at the centre, which is interpreted as the shadow of the black hole. This breakthrough result, which represents a powerful confirmation of Einstein's theory of gravity, or general relativity, was made possible by assembling a global network of radio telescopes operating at millimetre wavelengths that for the first time included the Atacama Large Millimeter/ submillimeter Array (ALMA). The addition of ALMA as an anchor station has enabled a giant leap forward by increasing the sensitivity limits of the EHT by an order of magnitude, effectively turning it into an imaging array. The published image demonstrates that it is now possible to directly study the event horizon shadows of SMBHs via electromagnetic radiation, thereby transforming this elusive frontier from a mathematical concept into an astrophysical reality. The expansion of the array over the next few years will include new stations on different continents - and eventually satellites in space. This will provide progressively sharper and higher-fidelity images of SMBH candidates, and potentially even movies of the hot plasma orbiting around SMBHs. These improvements will shed light on the processes of black hole accretion and jet formation on event-horizon scales, thereby enabling more precise tests of general relativity in the truly strong field regime.

We present a detailed description of the special procedures for calibration and quality assurance of Atacama Large Millimeter/submillimeter Array (ALMA) observations in Very Long Baseline Interferometry (VLBI) mode. These procedures are required to turn the phased ALMA array into a fully calibrated VLBI station. As an illustration of these methodologies, we present full-polarization observations carried out with ALMA as a phased array at 3mm (Band 3) and 1.3mm (Band 6) as part of Cycle-4. These are the first VLBI science observations conducted with ALMA and were obtained during a 2017 VLBI campaign in concert with other telescopes worldwide as part of the Global mm-VLBI Array (GMVA, April 1-3) and the Event Horizon Telescope (EHT, April 5-11) in ALMA Bands 3 and 6, respectively.

Polarization from active galactic nuclei is interpreted as a signpost of the role of magnetic fields in the launch and collimation of their relativistic radio jets. Here, we report the detection of a clear polarization signal from ALMA observations of the gravitationally lensed quasar PKS 1830-211 at submillimeter wavelengths (Band 9, 650 GHz). Applying a differential-polarimetry technique to the two compact lensed images of the quasar, we estimate a fractional polarization of ~5% for one lensed image, while the other appears nearly unpolarized, which implies that the polarization activity varies on a timescale of a few weeks. With additional ALMA Band 7 and 8 (between 300-500 GHz) concomitant data, we constrain a Faraday rotation of a few $10^5$ rad m$^{-2}$. We also observe flux-density variability of ~10% within one hour in Band 9. This work illustrates that a differential analysis can extract high-accuracy information (flux-density ratio and polarimetry) free of calibration issues from resolved sources in the submillimeter domain.

Context. With the outstanding spatial resolution and sensitivity of the Atacama Large Millimeter/sub-millimeter Array (ALMA), molecular gas other than the abundant CO can be observed and resolved in circumstellar envelopes (CSEs) around evolved stars, such as the binary S-type Asymptotic Giant Branch (AGB) star W Aquilae. Aims. We aim to constrain the chemical composition of the CSE and determine the radial abundance distribution, the photospheric peak abundance, and isotopic ratios of a selection of chemically important molecular species in the innermost CSE of W Aql. The derived parameters are put into the context of the chemical evolution of AGB stars and are compared with theoretical models. Methods. We employ one-dimensional radiative transfer modeling - with the accelerated lambda iteration (ALI) radiative transfer code - of the radial abundance distribution of a total of five molecular species (CS, SiS, 30SiS, 29SiO and H13CN) and determine the best fitting model parameters based on high-resolution ALMA observations as well as archival single-dish observations. The additional advantage of the spatially resolved ALMA observations is that we can directly constrain the radial profile of the observed line transitions from the observations. Results. We derive abundances and e-folding radii for CS, SiS, 30SiS, 29SiO and H13CN and compare them to previous studies, which are based only on unresolved single-dish spectra. Our results are in line with previous results and are more accurate due to resolution of the emission regions.

We report on the first VLBI image of the M81 nucleus (M81*) at a wavelength of 3.4 mm, obtained with the source-frequency phase-referencing (SFPR) technique. Thanks to the SFPR calibration, the coherent integration time could be eventually increased by more than an order of magnitude, which enabled the detection of fringes at the level of 45 mJy beam-1 with a dynamic range higher than 130:1. This paves the way toward future mm/sub-mm VLBI observations of weaker sources. From the analysis of the M81* visibilities, a core size of ~50 uas at 3.4 mm was estimated. This follows the power-law relationship with wavelength lambda^0.88, reported previously at lower frequencies. These results constrain the core size (at 3.4 mm) to a minimum of ~80 Schwarzschild radii of M81*.

The Atacama Millimeter/submillimeter Array (ALMA) Phasing Project (APP) has developed and deployed the hardware and software necessary to coherently sum the signals of individual ALMA antennas and record the aggregate sum in Very Long Baseline Interferometry (VLBI) Data Exchange Format. These beamforming capabilities allow the ALMA array to collectively function as the equivalent of a single large aperture and participate in global VLBI arrays. The inclusion of phased ALMA in current VLBI networks operating at (sub)millimeter wavelengths provides an order of magnitude improvement in sensitivity, as well as enhancements in u-v coverage and north-south angular resolution. The availability of a phased ALMA enables a wide range of new ultra-high angular resolution science applications, including the resolution of supermassive black holes on event horizon scales and studies of the launch and collimation of astrophysical jets. It also provides a high-sensitivity aperture that may be used for investigations such as pulsar searches at high frequencies. This paper provides an overview of the ALMA Phasing System design, implementation, and performance characteristics.

The APerture SYNthesis SIMulator is a simple interactive tool to help the students visualize and understand the basics of the Aperture Synthesis technique, applied to astronomical interferometers. The users can load many different interferometers and source models (and also create their own), change the observing parameters (e.g., source coordinates, observing wavelength, antenna location, integration time, etc.), and even deconvolve the interferometric images and corrupt the data with gain errors (amplitude and phase). The program is fully interactive and all the figures are updated in real time. APSYNSIM has already been used in several interferometry schools and has got very positive feedback from the students.

We use Atacama Large Millimeter/submillimeter Array Band 5 science verification observations of the red supergiant VY CMa to study the polarization of SiO thermal/masers lines and dust continuum at ~1.7 mm wavelength. We analyse both linear and circular polarization and derive the magnetic field strength and structure, assuming the polarization of the lines originates from the Zeeman effect, and that of the dust originates from aligned dust grains. We also discuss other effects that could give rise to the observed polarization. We detect, for the first time, significant polarization (~3%) of the circumstellar dust emission at millimeter wavelengths. The polarization is uniform with an electric vector position angle of $\sim8^\circ$. Varying levels of linear polarization are detected for the J=4-3 28SiO v=0, 1, 2, and 29SiO v=0, 1 lines, with the strongest polarization fraction of ~30% found for the 29SiO v=1 maser. The linear polarization vectors rotate with velocity, consistent with earlier observations. We also find significant (up to ~1%) circular polarization in several lines, consistent with previous measurements. We conclude that the detection is robust against calibration and regular instrumental errors, although we cannot yet fully rule out non-standard instrumental effects. Emission from magnetically aligned grains is the most likely origin of the observed continuum polarization. This implies that the dust is embedded in a magnetic field >13 mG. The maser line polarization traces the magnetic field structure. The magnetic field in the gas and dust is consistent with an approximately toroidal field configuration, but only higher angular resolution observations will be able to reveal more detailed field structure. If the circular polarization is due to Zeeman splitting, it indicates a magnetic field strength of ~1-3 Gauss, consistent with previous maser observations.

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.

Current mm/submm interferometers, like the Atacama Large mm/submm Array (ALMA), use receivers that register the sky signal in a linear polarization basis. In the case of observations performed in full-polarization mode (where the cross-correlations are computed among all the polarization channels) it is possible to reconstruct the full-polarization brightness distribution of the observed sources, as long as a proper calibration of delay offsets and leakage among polarization channels can be performed. Observations of calibrators, preferably with some linear polarization, with a good parallactic angle coverage are usually needed for such a calibration. In principle, dual-polarization observations only allow us to recover the Stokes $I$ intensity distribution of the sources, regardless of the parallactic angle coverage of the observations. In this paper, we present a novel technique of dual differential polarimetry that makes it possible to obtain information related to the full-polarization brightness distribution of the observed sources from dual-polarization observations. This technique is inspired in the Earth-rotation polarization synthesis and can be applied even to sources with complex structures.

The Gamma-ray BL Lac object OJ 287 is known to exhibit inner-parsec "jet-wobbling", high degrees of variability at all wavelengths and quasi-stationary features including an apparent (~100 deg) position angle change in projection on the sky plane. Sub-50 micro-arcsecond resolution 86 GHz observations with the global mm-VLBI array (GMVA) supplement ongoing multi-frequency VLBI blazar monitoring at lower frequencies. Using these maps together with cm/mm total intensity and Gamma-ray observations from Fermi/LAT from 2008-2014, we aimed to determine the location of Gamma-ray emission and to explain the inner-mas structural changes. Observations with the GMVA offer approximately double the angular resolution compared with 43 GHz VLBA observations and allow us to observe above the synchrotron self-absorption peak frequency. The jet was spectrally decomposed at multiple locations along the jet. From this we derived estimates of the magnetic field. How the field decreases down the jet allowed an estimate of the distance to the jet apex and an estimate of the magnetic field strength at the jet apex and in the broad line region. Combined with accurate kinematics we attempt to locate the site of Gamma-ray activity, radio flares and spectral changes. Strong Gamma-ray flares appeared to originate from either the "core" region, a downstream stationary feature, or both, with Gamma-ray activity significantly correlated with radio flaring in the downstream quasi-stationary feature. Magnetic field estimates were determined at multiple locations along the jet, with the magnetic field found to be >1.6 G in the "core" and >0.4 G in the downstream quasi-stationary feature. We therefore found upper limits on the location of the "core" as >6.0 pc from the jet apex and determined an upper limit on the magnetic field near the jet base of the order of thousands of Gauss.

We discuss about how the relative flux densities among the images of gravitationally-lensed active galactic nuclei, AGN, can be used to study the intrinsic AGN variability with high accuracy. Multi-frequency monitoring observations of resolved gravitational lenses can allow us to detect signals of very weak variability and also provide information about the jet opacity and structure. As an example, we investigate the variability of the flux-density ratio between the two lensed images of the blazar B0218+357, using dual-frequency cm-wave observations. Similar to our previously reported submm-wave observations of the lensed blazar PKS1830-211, we observe a clear chromatic variability, starting short before an increase in the flux-density of the blazar. The evolution of the flux-density ratios between the blazar images shows a more clear and rich structure than that of the mere lightcurves of each individual image. The accuracy in the ratio measurements is allowing us to see variability episodes in the blazar that are weaker than the natural scatter in the absolute flux-density measurements. A simple opacity model in the jet is used to consistently explain the difference between the flux-density-ratio evolution at the two frequencies.

Active galactic nuclei (AGN) host some of the most energetic phenomena in the Universe. AGN are thought to be powered by accretion of matter onto a rotating disk that surrounds a supermassive black hole. Jet streams can be boosted in energy near the event horizon of the black hole and then flow outward along the rotation axis of the disk. The mechanism that forms such a jet and guides it over scales from a few light-days up to millions of light-years remains uncertain, but magnetic fields are thought to play a critical role. Using the Atacama large mm/submm array (ALMA), we have detected a polarization signal (Faraday rotation) related to the strong magnetic field at the jet base of a distant AGN, PKS1830-211. The amount of Faraday rotation (rotation measure) is proportional to the magnetic field strength along the line of sight times the density of electrons. Although it is impossible to precisely infer the magnetic fields in the region of Faraday rotation, the high rotation measures derived suggest magnetic fields of at least tens of Gauss (and possibly considerably higher) on scales of the order of light days (0.01 pc) from the black hole.

Heterodyne receivers register the sky signal on either a circular polarization basis (where it is split into left-hand and right-hand circular polarization) or a linear polarization basis (where it is split into horizontal and vertical linear polarization). We study the problem of interferometric observations performed with telescopes that observe on different polarization bases, hence producing visibilities that we call "mixed basis" (i.e., linear in one telescope and circular in the other). We present novel algorithms for the proper calibration and treatment of such interferometric observations and test our algorithms with both simulations and real data. The use of our algorithms will be important for the optimum calibration of forthcoming observations with the Atacama Large mm/submm Array (ALMA) in very-long-baseline interferometry (VLBI) mode. Our algorithms will also allow us to optimally calibrate future VLBI observations at very high data rates (i.e., wide bandwidths), where linear-polarization feeds will be preferable at some stations, to overcome the polarimetric limitations due to the use of quarter-wave plates.

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.

The Atacama Large mm-submm Array (ALMA) is, by far, the most sensitive mm/submm telescope in the World. The ALMA Phasing Project (APP) will allow us to phase-up all the ALMA antennas and use them as one single VLBI station. This will be a key component of the Event Horizon Telescope (EHT), a Global VLBI array at millimeter wavelengths. A problem in the APP is the calibration and conversion of the polarization channels. Most VLBI stations record their signals in a circular basis, but the ALMA receivers record in a linear basis. The strategy that will be followed in the phased-ALMA VLBI observations will be to correlate in "mixed" basis (i.e., linear versus circular) and convert the visibilities to a pure circular basis after the correlation. We have developed an algorithm to perform such a polarization conversion of the VLBI visibilities. In these proceedings, we present the basics of the PolConvert algorithm and discuss on the polarization conversion in the general case were single dishes (besides phased arrays) record with linear receivers in VLBI observations. We show some results of PolConvert applied to realistic simulations, as well as a test with real VLBI observations at 86\,GHz between the Onsala radiotelescope (recording in linear basis) and the Effelsberg radiotelescope (recording in circular basis).

The Global mm-VLBI Array (GMVA) is a network of 14 3\,mm and 7\,mm capable telescopes spanning Europe and the United States, with planned extensions to Asia. The array is capable of sensitive maps with angular resolution often exceeding 50\,$\mu$as. Using the GMVA, a large sample of prominent $\gamma$-ray blazars have been observed approximately 6 monthly from later 2008 until now. Combining 3\,mm maps from the GMVA with near-in-time 7\,mm maps from the VLBA-BU-BLAZAR program and 2\,cm maps from the MOJAVE program, we determine the sub-pc morphology and high frequency spectral structure of $\gamma$-ray blazars. The magnetic field strength can be estimated at different locations along the jet under the assumption of equipartition between magnetic field and relativistic particle energies. Making assumptions on the jet magnetic field configuration (e.g. poloidal or toroidal), we can estimate the separation of the mm-wave "core" and the jet base, and estimate the strength of the magnetic field there. The results of this analysis show that on average, the magnetic field strength decreases with a power-law $B \propto r^{-n}$, $n=0.3 \pm 0.2$. This suggests that on average, the mm-wave "core" is $\sim 1-3$\,pc downstream of the de-projected jet apex and that the magnetic field strength is of the order $B_{\rm{apex}} \sim 5-20$\,kG, broadly consistent with the predictions of magnetic jet launching (e.g. via magnetically arrested disks (MAD)).

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.

Stacking as a tool for studying objects that are not individually detected is becoming popular even for radio interferometric data, and will be widely used in the SKA era. Stacking is typically done using imaged data rather than directly using the visibilities (the uv-data). We have investigated and developed a novel algorithm to do stacking using the uv-data. We have performed exten- sive simulations comparing to image-stacking, and summarize the results of these simulations. Furthermore, we disuss the implications in light of the vast data volume produced by the SKA. Having access to the uv-stacked data provides a great advantage, as it allows the possibility to properly analyse the result with respect to calibration artifacts as well as source properties such as size. For SKA the main challenge lies in archiving the uv-data. For purposes of robust stacking analysis, it would be strongly desirable to either keep the calibrated uv-data at least in an aver- age form, or implement a stacking queue where stacking positions could be provided prior to the observations and the uv-stacking is done almost in real time.

Adding VLBI capability to the SKA arrays will greatly broaden the science of the SKA, and is feasible within the current specifications. SKA-VLBI can be initially implemented by providing phased-array outputs for SKA1-MID and SKA1-SUR and using these extremely sensitive stations with other radio telescopes, and in SKA2 by realising a distributed configuration providing baselines up to thousands of km, merging it with existing VLBI networks. The motivation for and the possible realization of SKA-VLBI is described in this paper.

We present a comprehensive spectral and morphological analysis of the remnant of Supernova (SN) 1987A with the Australia Telescope Compact Array (ATCA) and the Atacama Large Millimeter/submillimeter Array (ALMA). The non-thermal and thermal components of the radio emission are investigated in images from 94 to 672 GHz ($\lambda$ 3.2 mm to 450 $\mu$m), with the assistance of a high-resolution 44 GHz synchrotron template from the ATCA, and a dust template from ALMA observations at 672 GHz. An analysis of the emission distribution over the equatorial ring in images from 44 to 345 GHz highlights a gradual decrease of the east-to-west asymmetry ratio with frequency. We attribute this to the shorter synchrotron lifetime at high frequencies. Across the transition from radio to far infrared, both the synchrotron/dust-subtracted images and the spectral energy distribution (SED) suggest additional emission beside the main synchrotron component ($S_{\nu}\propto\nu^{-0.73}$) and the thermal component originating from dust grains at $T\sim22$ K. This excess could be due to free-free flux or emission from grains of colder dust. However, a second flat-spectrum synchrotron component appears to better fit the SED, implying that the emission could be attributed to a pulsar wind nebula (PWN). The residual emission is mainly localised west of the SN site, as the spectral analysis yields $-0.4\lesssim\alpha\lesssim-0.1$ across the western regions, with $\alpha\sim0$ around the central region. If there is a PWN in the remnant interior, these data suggest that the pulsar may be offset westward from the SN position.

Carbon radio recombination lines (RRLs) at low frequencies (<=500 MHz) trace the cold, diffuse phase of the interstellar medium, which is otherwise difficult to observe. We present the detection of carbon RRLs in absorption in M82 with LOFAR in the frequency range of 48-64 MHz. This is the first extragalactic detection of RRLs from a species other than hydrogen, and below 1 GHz. Since the carbon RRLs are not detected individually, we cross-correlated the observed spectrum with a template spectrum of carbon RRLs to determine a radial velocity of 219 +- 9 km/s . Using this radial velocity, we stack 22 carbon-{\alpha} transitions from quantum levels n = 468-508 to achieve an 8.5 sigma detection. The absorption line profile exhibits a narrow feature with peak optical depth of 0.003 and FWHM of 31 km/s. Closer inspection suggests that the narrow feature is superimposed on a broad, shallow component. The total line profile appears to be correlated with the 21 cm H I line profile reconstructed from H I absorption in the direction of supernova remnants in the nucleus. The narrow width and centroid velocity of the feature suggests that it is associated with the nuclear starburst region. It is therefore likely that the carbon RRLs are associated with cold atomic gas in the direction of the nucleus of M82.

The Global mm-VLBI Array (GMVA) is the highest angular resolution imaging interferometer currently available as a common user facility. It is capable of angular resolutions on the order of 40 microarcseconds. Currently 14 stations in the United States and Europe participate in global 3 mm VLBI observations. The GMVA is used for continuum and spectroscopic imaging, probing the central regions of active galaxies and the origin of jets as these regions are typically not observable at longer wavelengths due to synchrotron self-absorption. In early 2012, fringes were detected to the three stations of the Korean VLBI Network (KVN), opening the possibility of extending the baseline coverage of the VLBI array to the East. In these proceedings, we will present recent images from a monitoring program of gamma-ray blazars using the GMVA, including the sources 3C454.3 and 0235+164, and an update of its current status and abilities.

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 the first results of an ALMA spectral survey of strong absorption lines for common interstellar species in the z=0.89 molecular absorber toward the lensed blazar PKS1830-211. The dataset brings essential information on the structure and composition of the absorbing gas in the foreground galaxy. In particular, we find absorption over large velocity intervals (gtrsim 100 km/s) toward both lensed images of the blazar. This suggests either that the galaxy inclination is intermediate and that we sample velocity gradients or streaming motions in the disk plane, that the molecular gas has a large vertical distribution or extraplanar components, or that the absorber is not a simple spiral galaxy but might be a merger system. The number of detected species is now reaching a total of 42 different species plus 14 different rare isotopologues toward the SW image, and 14 species toward the NE line-of-sight. The abundances of CH, H2O, HCO+, HCN, and NH3 relative to H2 are found to be comparable to those in the Galactic diffuse medium. Of all the lines detected so far toward PKS1830-211, the ground-state line of ortho-water has the deepest absorption. We argue that ground-state lines of water have the best potential for detecting diffuse molecular gas in absorption at high redshift.

The analysis of astronomical interferometric data is often performed on the images obtained after deconvolution of the interferometer's point spread function (PSF). This strategy can be understood (especially for cases of sparse arrays) as fitting models to models, since the deconvolved images are already non-unique model representations of the actual data (i.e., the visibilities). Indeed, the interferometric images may be affected by visibility gridding, weighting schemes (e.g., natural vs. uniform), and the particulars of the (non-linear) deconvolution algorithms. Fitting models to the direct interferometric observables (i.e., the visibilities) is preferable in the cases of simple (analytical) sky intensity distributions. In this paper, we present UVMULTIFIT, a versatile library for fitting visibility data, implemented in a Python-based framework. Our software is currently based on the CASA package, but can be easily adapted to other analysis packages, provided they have a Python API. We have tested the software with synthetic data, as well as with real observations. In some cases (e.g., sources with sizes smaller than the diffraction limit of the interferometer), the results from the fit to the visibilities (e.g., spectra of close by sources) are far superior to the output obtained from the mere analysis of the deconvolved images. UVMULTIFIT is a powerful improvement of existing tasks to extract the maximum amount of information from visibility data, especially in cases close to the sensitivity/resolution limits of interferometric observations.

Supernova (SN) explosions are crucial engines driving the evolution of galaxies by shock heating gas, increasing the metallicity, creating dust, and accelerating energetic particles. In 2012 we used the Atacama Large Millimeter/Submillimeter Array to observe SN 1987A, one of the best-observed supernovae since the invention of the telescope. We present spatially resolved images at 450um, 870um, 1.4mm, and 2.8mm, an important transition wavelength range. Longer wavelength emission is dominated by synchrotron radiation from shock-accelerated particles, shorter wavelengths by emission from the largest mass of dust measured in a supernova remnant (>0.2Msun). For the first time we show unambiguously that this dust has formed in the inner ejecta (the cold remnants of the exploded star's core). The dust emission is concentrated to the center of the remnant, so the dust has not yet been affected by the shocks. If a significant fraction survives, and if SN 1987A is typical, supernovae are important cosmological dust producers.

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.

The launching mechanism of the jets of active galactic nuclei is observationally poorly constrained, due to the large distances to these objects and the very small scales (sub-parsec) involved. In order to better constrain theoretical models, it is especially important to get information from the region close to the physical base of the jet, where the plasma acceleration takes place. In this paper, we report multi-epoch and multi-frequency continuum observations of the z=2.5 blazar PKS1830-211 with ALMA, serendipitously coincident with a strong $\gamma$-ray flare reported by Fermi-LAT. The blazar is lensed by a foreground z=0.89 galaxy, with two bright images of the compact core separated by 1". Our ALMA observations individually resolve these two images (although not any of their substructures), and we study the change of their relative flux ratio with time (four epochs spread over nearly three times the time delay between the two lensed images) and frequency (between 350 and 1050 GHz, rest-frame of the blazar), during the $\gamma$-ray flare. In particular, we detect a remarkable frequency-dependent behaviour of the flux ratio, which implies the presence of a chromatic structure in the blazar (i.e., a core-shift effect). We rule out the possiblity of micro- and milli-lensing effects and propose instead a simple model of plasmon ejection in the blazar's jet to explain the time and frequency variability of the flux ratio. We suggest that PKS1830-211 is likely one of the best sources to probe the activity at the base of a blazar's jet at submillimeter wavelengths, thanks to the peculiar geometry of the system. The implications of the core-shift in absorption studies of the foreground z=0.89 galaxy (e.g., constraints on the cosmological variations of fundamental constants) are discussed.

SN1993J in M81 is the best studied young radio-luminous supernova in the Northern Hemisphere. We recently reported results from the analysis of a complete set of VLBI observations of this supernova at 1.7, 2.3, 5.0, and 8.4 GHz, covering a time baseline of more than one decade. Those reported results were focused on the kinematics of the expanding shock, the particulars of its evolving non-thermal emission, the density profile of the circumstellar medium, and the evolving free-free opacity by the supernova ejecta. In the present paper, we complete our analysis by performing a search for any possible signal from a compact source (i.e., a stellar-mass black hole or a young pulsar nebula) at the center of the expanding shell. We have performed a stacking of all our VLBI images at each frequency, after subtraction of our best-fit shell model at each epoch, and measured the peak intensity in the stacked residual image. Given the large amount of available global VLBI observations, the stacking of all the residual images allows us to put upper limits to the eventual emission of a putative compact central source at the level of $\sim102$ $\mu$Jy at 5 GHz (or, more conservatively, $\sim192$ $\mu$Jy, if we make a further correction for the ejecta opacity) and somewhat larger at other wavelengths.

We report spectroscopic and imaging observations of rotational transitions of cold CO and SiO in the ejecta of SN1987A, the first such emission detected in a supernova remnant. In addition to line luminosities for the CO J=1-0, 2-1, 6-5, and 7-6 transitions, we present upper limits for all other transitions up to J=13-12, collectively measured from the Atacama Large Millimeter Array (ALMA), the Atacama Pathfinder EXperiment (APEX), and the Herschel Spectral and Photometric Imaging REceiver (SPIRE). Simple models show the lines are emitted from at least 0.01 solar masses of CO at a temperature > 14 K, confined within at most 35% of a spherical volume expanding at ~ 2000 km/s. Moreover, we locate the emission within 1'' of the central debris. These observations, along with a partial observation of SiO, confirm the presence of cold molecular gas within supernova remnants and provide insight into the physical conditions and chemical processes in the ejecta. Furthermore, we demonstrate the powerful new window into supernova ejecta offered by submillimeter observations.

In a recent publication, we reported results of a multi-frequency VLBI campaign of observations of the Active Galactic Nucleus (AGN) in galaxy M\,81, phase-referenced to the supernova SN\,1993J. We were able to extract precise information on the relative astrometry of the AGN radio emission at different epochs and frequencies. We found strong evidence of precession in the AGN jet (i.e., a systematic evolution in the jet inclination at each frequency) coupled to changes in the overall flux density at the different frequencies. In these proceedings, we summarise the main contents of our previous publication and we report on (preliminary) new results from our follow-up VLBI observations, now phase-referenced to the young supernova SN2008iz. We also briefly discuss how these results match the picture of our previously-reported precession model.

We make use of VLBI observations of the radio jet in the quasar S5 0836+710 at different frequencies and epochs to study its properties. The jet shows helical structure at all frequencies. The ridge-line of the emission in the jet coincides at all frequencies and epochs, within the errors. We conclude that the helicity is a real, physical structure. Small differences between epochs reveal wave-like motion of the ridge-line transversal to the jet propagation axis. These transversal motions are measured to be superluminal. This unphysical result could correspond to a possible small amplitude oscillation of the ridge-line at the radio-core and to large errors in the determination of the positions. In addition, higher resolution images at 15 GHz show that the ridge-line does not coincide exactly with the centre of the radio jet. At arc-second scales, this powerful jet shows non-collimated, irregular structure and a lack of a hot-spot. Following this collection of evidence, we conclude that the ridge-line could be related to a pressure maximum within the jet cross-section, associated with the observed helical pattern that could lead to jet disruption at longer scales. We discuss the implications regarding the nature of the instability and the jet flow itself.

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 the development of a semi-automatic pipeline for the calibration of 86 GHz full-polarization observations performed with the Global Millimeter-VLBI array (GMVA) and describe the calibration strategy followed in the data reduction. Our calibration pipeline involves non-standard procedures, since VLBI polarimetry at frequencies above 43 GHz is not yet well established. We also present, for the first time, a full-polarization global-VLBI image at 86 GHz (source 3C 345), as an example of the final product of our calibration pipeline, and discuss the effect of instrumental limitations on the fidelity of the polarization images. Our calibration strategy is not exclusive for the GMVA, and could be applied on other VLBI arrays at millimeter wavelengths. The use of this pipeline will allow GMVA observers to get fully-calibrated datasets shortly after the data correlation.

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.

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

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 present first results from the ongoing radio monitoring of SN 2008iz in M82. The VLBI images reveal a shell-like structure with circular symmetry, which expands in a self-similar way. There is strong evidence of a compact component with a steep spectrum at the center of the shell. The expansion curve obtained from our VLBI observations is marginally decelerated (m = 0.89) and can be modeled simultaneously with the available radio light curves. While the results of this simultaneous fitting are not conclusive (i.e. different combinations of values of the magnetic field, CSM density profile, and electron energy distribution, provide fits to the available data with similar quality), additional observations should allow a more robust and detailed modeling.

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.

Modern radio interferometers sensitive to low frequencies will make use of wide-band detectors. For such wide bandwidths, dispersive atmospheric effects introduce variations in the fringe delay which change through the band of the receivers. These undesired dispersive effects must be estimated and calibrated with the highest precision. We studied the achievable precision in the estimate of the ionospheric dispersion and the dynamic range of the correlated fringes for different distributions of sub-bands in low-frequency and wide-band interferometric observations. Our study is focused on the case of sub-bands with a bandwidth much narrower than that of the total covered spectrum (case of LOFAR). We computed the uncertainty of the ionospheric delay, the delay ambiguity, and the dynamic range of the fringes using four different kinds of sub-band distributions: constant spacing between sub-bands, random spacings, spacings based on a power-law distribution, and spacings based on Golomb rulers (sets of integers whose sets of differences have non-repeated elements). For a large number of sub-bands ($> 20$, depending on the delay window) spacings based on Golomb rulers give the most precise estimates of dispersive effects and the highest fringe dynamic ranges. Spacings based on the power-law distribution give similar results, although better than those with the Golomb rulers for smaller number of sub-bands. Random distributions result in large fringe dynamic ranges, but the estimate of dispersive effects is worse. A constant spacing of sub-bands results in very bad fringe dynamic ranges, but good estimates of ionospheric dispersion. Combining all the results, the power-law distribution gives the best compromise between homogeneity in the bandwidth sampling, precision in the estimate of ionospheric effects, dynamic range of the correlated fringes, and group-delay ambiguity.

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 present observations of the recently discovered supernova 2008iz in M82 with the VLBI High Sensitivity Array at 22 GHz, the Very Large Array at frequencies of 1.4, 4.8, 8.4, 22 and 43 GHz, and the Chandra X-ray observatory. The supernova was clearly detected on two VLBI images, separated by 11 months. The source shows a ring-like morphology and expands with a velocity of ~23000 km/s. The most likely explosion date is in mid February 2008. The measured expansion speed is a factor of ~2 higher than expected under the assumption that synchrotron self-absorption dominates the light curve at the peak, indicating that this absorption mechanism may not be important for the radio emission. We find no evidence for an asymmetric explosion. The VLA spectrum shows a broken power law, indicating that the source was still optically thick at 1.4 GHz in April 2009. Finally, we report upper limits on the X-ray emission from SN 2008iz and a second radio transient recently discovered by MERLIN observations.

We report on a set of 5 GHz Urumqi observations of the galaxy M82, made between August 2005 and May 2009. From the resulting flux densities, we detect a strong flare, starting in March or April and peaking in June 2008. We identify this flare with supernova SN 2008iz. The time sampling of the radio light curve allows us to obtain information on the precursor mass-loss rate, the strength of the magnetic field in the radiating region, the explosion date, and the deceleration of the expanding shock. We also check the possible contribution of Synchrotron Self Absorption (SSA) to the radio light curve and compare our model with other observations of the supernova at 22 GHz.

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.

Phase self-calibration (or selfcal) is an algorithm often used in the calibration of interferometric observations in astronomy. Although a powerful tool, this algorithm presents strong limitations when applied to data with a low signal-to-noise ratio. We analyze the artifacts that the phase selfcal algorithm produces when applied to extremely noisy data. We show how the phase selfcal may generate a spurious source in the sky from a distribution of completely random visibilities. This spurious source is indistinguishable from a real one. We numerically and analytically compute the relationship between the maximum spurious flux density generated by selfcal from noise and the particulars of the interferometric observations. Finally, we present two simple tests that can be applied to interferometric data for checking whether a source detection is real or whether the source is an artifact of the phase self-calibration algorithm.

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 the results of a second epoch of 18cm global Very Long-Baseline Interferometry (VLBI) observations, taken on 23 February 2001, of the central kiloparsec of the nearby starburst galaxy Messier 82. These observations further investigate the structural and flux evolution of the most compact radio sources in the central region of M82. The two most compact radio objects in M82 have been investigated (41.95+575 and 43.31+592). Using this recent epoch of data in comparison with our previous global VLBI observations and two earlier epochs of European VLBI Network observations we measure expansion velocities in the range of 1500-2000km/s for 41.95+575, and 9000-11000km/s for 43.31+592 using various independent methods. In each case the measured remnant expansion velocities are significantly larger than the canonical expansion velocity (500km/s) of supernova remnants within M82 predicted from theoretical models. In this paper we discuss the implications of these measured expansion velocities with respect to the high density environment that the SNR are expected to reside in within the centre of the M82 starburst.

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