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.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

The AGB star $\pi^{1}$ Gruis has a known companion (at a separation of ~400 AU) which cannot explain the strong deviations from the spherical symmetry of the CSE. Recently, hydrodynamic simulations of mass transfer in closer binary systems have successfully reproduced the spiral-shaped CSEs found around a handful of sources. There is growing evidence for an even closer, undetected companion complicating the case of $\pi^{1}$ Gruis further. The improved spatial resolution allows for the investigation of the complex circumstellar morphology and the search for imprints on the CSE of the third component. We have observed the 12CO J=3-2 line emission from $\pi^{1}$ Gruis using both the compact and extended array of Atacama Large Millimeter/submillimeter Array (ALMA). The interferometric data has furthermore been combined with data from the ALMA total power (TP) array. The imaged brightness distribution has been used to constrain a non-local, non-LTE 3D radiative transfer model of the CSE. The high-angular resolution ALMA data have revealed the first example of a source on the AGB where both a faster bipolar outflow and a spiral pattern along the orbital plane can be seen in the gas envelope. The spiral can be traced in the low- to intermediate velocity, from 13 to 25 km s$^{-1}$, equatorial torus. The largest spiral-arm separation is $\approx$5".5 and consistent with a companion with an orbital period of $\approx$330 yrs and a separation of less than 70 AU. The kinematics of the bipolar outflow is consistent with it being created during a mass-loss eruption where the mass-loss rate from the system increased by at least a factor of 5 during 10-15 yrs. The spiral pattern is the result of an undetected companion. The bipolar outflow is the result of a rather recent mass-loss eruption event.

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.

This paper presents a survey of microwave front-end receivers installed at radio telescopes throughout the World. This unprecedented analysis was conducted as part of a review of front-end developments for Italian radio telescopes, initiated by the Italian National Institute for Astrophysics in 2016. Fifteen international radio telescopes have been selected to be representative of the instrumentation used for radio astronomical observations in the frequency domain from 300 MHz to 116 GHz. A comprehensive description of the existing receivers is presented and their characteristics are compared and discussed. The observing performances of the complete receiving chains are also presented. An overview of on-going developments illustrates and anticipates future trends in front-end projects to meet the most ambitious scientific research goals.

OH231.8+4.2 is an archetypal pre-planetary nebulae (pPN). It is a binary system surrounded by bipolar nebula. Some years ago the authors extensively studied it and performed several VLBI observations from which they obtained mas-resolution maps of the SiO (7 mm) and H2O (1.3 cm) maser emissions. H2O masers were found to be distributed in two areas along the symmetry axis of the nebulae oriented nearly north-south delineating a bipolar outflow and their astrometric positions were accurately measured. SiO masers, indicating the position of the Mira component of the binary system, form a structure perpendicular to the axis of the nebulae. The general picture of the source looked satisfactory, except for the relative position of the two masers. Surprisingly, SiO masers, were tentatively placed 250 mas away (370 AU) from the apparent center of the outflow. Using the ALMA we observed the SiO maser emission at 86 GHz and accurately derived the position of the Mira component. Combining our previous VLBA data and our new ALMA observations we are now able to give a more complete and detailed description of the inner part of this amazing pPN.

We have performed a study of the characteristics of the circumstellar environment of the binary object HD101584, that provides information on a likely evolutionary scenario. We have obtained and analysed ALMA observations, complemented with observations using APEX, of a large number of molecular lines. An analysis of the spectral energy distribution has also been performed. Emissions from 12 molecular species (not counting isotopologues) have been observed, and most of them mapped with angular resolutions in the range 0.1" to 0.6". Four circumstellar components are identified: i) a central compact source of size 0.15", ii) an expanding equatorial density enhancement (a flattened density distribution in the plane of the orbit) of size 3", iii) a bipolar high-velocity outflow (150 km/s), and iv) an hourglass structure. The outflow is directed almost along the line of sight. There is evidence of a second bipolar outflow. The mass of the circumstellar gas is 0.5[D/1 kpc]^2 Msun, about half of it lies in the equatorial density enhancement. The dust mass is 0.01[D/1 kpc]^2 Msun, and a substantial fraction of this is in the form of large-sized, up to 1 mm, grains. The estimated kinetic age of the outflow is 770[D/1 kpc] yr. The kinetic energy and the scalar momentum of the accelerated gas are estimated to be 7x10^(45)[D/1 kpc]^2 erg and 10^(39)[D/1 kpc]^2 g cm/s, respectively. We provide good evidence that the binary system HD101584 is in a post-common-envelope-evolution phase, that ended before a stellar merger. Isotope ratios combined with stellar mass estimates suggest that the primary star's evolution was terminated already on the first red giant branch (RGB). Most of the energy required to drive the outflowing gas was probably released when material fell towards the companion.

We report results from a deep polarization imaging of the nearby radio galaxy 3C$\,$84 (NGC$\,$1275). The source was observed with the Global Millimeter VLBI Array (GMVA) at 86$\,$GHz at an ultra-high angular resolution of $50\mu$as (corresponding to 250$R_{s}$). We also add complementary multi-wavelength data from the Very Long Baseline Array (VLBA; 15 & 43$\,$GHz) and from the Atacama Large Millimeter/submillimeter Array (ALMA; 97.5, 233.0, and 343.5$\,$GHz). At 86$\,$GHz, we measure a fractional linear polarization of $\sim2$% in the VLBI core region. The polarization morphology suggests that the emission is associated with an underlying limb-brightened jet. The fractional linear polarization is lower at 43 and 15$\,$GHz ($\sim0.3-0.7$% and $<0.1$%, respectively). This suggests an increasing linear polarization degree towards shorter wavelengths on VLBI scales. We also obtain a large rotation measure (RM) of $\sim10^{5-6}~{\rm rad/m^{2}}$ in the core at $\gtrsim$43$\,$GHz. Moreover, the VLBA 43$\,$GHz observations show a variable RM in the VLBI core region during a small flare in 2015. Faraday depolarization and Faraday conversion in an inhomogeneous and mildly relativistic plasma could explain the observed linear polarization characteristics and the previously measured frequency dependence of the circular polarization. Our Faraday depolarization modeling suggests that the RM most likely originates from an external screen with a highly uniform RM distribution. To explain the large RM value, the uniform RM distribution, and the RM variability, we suggest that the Faraday rotation is caused by a boundary layer in a transversely stratified jet. Based on the RM and the synchrotron spectrum of the core, we provide an estimate for the magnetic field strength and the electron density of the jet plasma.

(abridged) Very long baseline interferometry (VLBI) observations at 86$\,$GHz (wavelength, $\lambda = 3\,$mm) reach a resolution of about 50 $\mu$as, probing the collimation and acceleration regions of relativistic outflows in active galactic nuclei. To extend the statistical studies of compact extragalactic jets, a large global 86 GHz VLBI survey of 162 radio sources was conducted in 2010-2011 using the Global Millimeter VLBI Array. The survey data attained a typical baseline sensitivity of 0.1 Jy and a typical image sensitivity of 5 mJy/beam, providing successful detections and images for all of the survey targets. For 138 objects, the survey provides the first ever VLBI images made at 86 GHz. Gaussian model fitting of the visibility data was applied to represent the structure of the sources. The Gaussian model-fit-based estimates of brightness temperature ($T_\mathrm{b}$) at the jet base (core) and in moving regions (jet components) downstream from the core were compared to the estimates of $T_\mathrm{b}$ limits made directly from the visibility data, demonstrating a good agreement between the two methods. The apparent brightness temperature estimates for the jet cores in our sample range from $2.5 \times 10^{9}\,$K to $ 1.3\times 10^{12}\,$K. A population model with a single intrinsic value of brightness temperature, $T_\mathrm{0}$, is applied to reproduce the observed $T_\mathrm{b}$ distribution. It yields $T_\mathrm{0} = (3.77^{+0.10}_{-0.14}) \times 10^{11}\,$K for the jet cores, implying that the inverse Compton losses dominate the emission. In the jet components, $T_\mathrm{0} =(1.42^{+0.16}_{-0.19})\times 10^{11}\,$K is found, slightly higher than the equipartition limit of $\sim5\times 10^{10}\,$K expected for these jet regions. For objects with sufficient structural detail detected, the adiabatic energy losses dominate the observed changes of $T_\mathrm{b}$ along the jet.

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 results from very long baseline interferometric (VLBI) observations of the supermassive black hole in the Galactic center, Sgr A*, at 1.3 mm (230 GHz). The observations were performed in 2013 March using six VLBI stations in Hawaii, California, Arizona, and Chile. Compared to earlier observations, the addition of the APEX telescope in Chile almost doubles the longest baseline length in the array, provides additional {\it uv} coverage in the N-S direction, and leads to a spatial resolution of $\sim$30 $\mu$as ($\sim$3 Schwarzschild radii) for Sgr A*. The source is detected even at the longest baselines with visibility amplitudes of $\sim$4-13% of the total flux density. We argue that such flux densities cannot result from interstellar refractive scattering alone, but indicate the presence of compact intrinsic source structure on scales of $\sim$3 Schwarzschild radii. The measured nonzero closure phases rule out point-symmetric emission. We discuss our results in the context of simple geometric models that capture the basic characteristics and brightness distributions of disk- and jet-dominated models and show that both can reproduce the observed data. Common to these models are the brightness asymmetry, the orientation, and characteristic sizes, which are comparable to the expected size of the black hole shadow. Future 1.3 mm VLBI observations with an expanded array and better sensitivity will allow a more detailed imaging of the horizon-scale structure and bear the potential for a deep insight into the physical processes at the black hole boundary.

M87 is one of the nearest radio galaxies with a prominent jet extending from sub-pc to kpc-scales. Because of its proximity and large mass of the central black hole, it is one of the best radio sources to study jet formation. We aim at studying the physical conditions near the jet base at projected separations from the BH of $\sim7-100$ Schwarzschild radii ($R_{\rm sch}$). Global mm-VLBI Array (GMVA) observations at 86 GHz ($\lambda=3.5\,$mm) provide an angular resolution of $\sim50\mu$as, which corresponds to a spatial resolution of only $7~R_{\rm sch}$ and reach the small spatial scale. We use five GMVA data sets of M87 obtained during 2004--2015 and present new high angular resolution VLBI maps at 86GHz. In particular, we focus on the analysis of the brightness temperature, the jet ridge lines, and the jet to counter-jet ratio. The imaging reveals a parabolically expanding limb-brightened jet which emanates from a resolved VLBI core of $\sim(8-13) R_{\rm sch}$ size. The observed brightness temperature of the core at any epoch is $\sim(1-3)\times10^{10}\,$K, which is below the equipartition brightness temperature and suggests magnetic energy dominance at the jet base. We estimate the diameter of the jet at its base to be $\sim5 R_{\rm sch}$ assuming a self-similar jet structure. This suggests that the sheath of the jet may be anchored in the very inner portion of the accretion disk. The image stacking reveals faint emission at the center of the edge-brightened jet on sub-pc scales. We discuss its physical implication within the context of the spine-sheath structure of the jet.

The S-type asymptotic giant branch (AGB) star $\pi^{1}$ Gruis has a known companion at a separation of $\approx$400 AU. The envelope structure, including an equatorial torus and a fast bipolar outflow, is rarely seen in the AGB phase and is particularly unexpected in such a wide binary system. Therefore a second, closer companion has been suggested, but the evidence is not conclusive. The new ALMA $^{12}$CO and $^{13}$CO $J$=3-2 data, together with previously published $^{12}$CO $J$=2-1 data from the Submillimeter Array (SMA), and the $^{12}$CO $J$=5-4 and $J$=9-8 lines observed with Herschel/Heterodyne Instrument for the Far-Infrared (HIFI), is modeled with the 3D non-LTE radiative transfer code SHAPEMOL. The data analysis clearly confirms the torus-bipolar structure. The 3D model of the CSE that satisfactorily reproduces the data consists of three kinematic components: a radially expanding torus with velocity slowly increasing from 8 to 13 km s$^{-1}$ along the equator plane; a radially expanding component at the center with a constant velocity of 14 km s$^{-1}$; and a fast, bipolar outflow with velocity proportionally increasing from 14 km s$^{-1}$ at the base up to 100 km s$^{-1}$ at the tip, following a linear radial dependence. The results are used to estimate an average mass-loss rate during the creation of the torus of 7.7$\times$10$^{-7}$ M$_{\odot}$ yr$^{-1}$. The total mass and linear momentum of the fast outflow are estimated at 7.3$\times$10$^{-4}$ M$_{\odot}$ and 9.6$\times$10$^{37}$ g cm s$^{-1}$, respectively. The momentum of the outflow is in excess (by a factor of about 20) of what could be generated by radiation pressure alone, in agreement with recent findings for more evolved sources. The best-fit model also suggests a $^{12}$CO/$^{13}$CO abundance ratio of 50. Possible shaping scenarios for the gas envelope are discussed

The CO(J=3-2) emission from the CSE of the binary S-type AGB star W Aql has been observed at subarcsecond resolution using ALMA. The aim of this paper is to investigate the wind properties of the AGB star and to analyse how the known companion has shaped the CSE. The average mass-loss rate during the creation of the detected CSE is estimated through modelling, using the ALMA brightness distribution and previously published single-dish measurements as observational constraints. The ALMA observations are presented and compared to the results from a 3D smoothed particle hydrodynamics (SPH) binary interaction model with the same properties as the W Aql system and with two different orbital eccentricities. Three-dimensional radiative transfer modelling is performed and the response of the interferometer is modelled and discussed. The estimated average mass-loss rate of W~Aql agrees with previous results. The size of the emitting region is consistent with photodissociation models. The CO(J=3-2) emission is dominated by a smooth component overlayed with two weak arc patterns with different separations. The larger pattern is predicted by the binary interaction model with separations of 10" and therefore likely due to the known companion. It is consistent with a binary orbit with low eccentricity. The smaller separation pattern is asymmetric and coincides with the dust distribution, but the separation timescale (200 yrs) is not consistent with any known process of the system. The separation of the known companions of the system is large enough to not have a very strong effect on the circumstellar morphology. The density contrast across the envelope of a binary with an even larger separation will not be easily detectable, even with ALMA, unless the orbit is strongly asymmetric or the AGB star has a much larger mass-loss rate.

Aims. Our goal is to characterize the intermediate age, detached shell carbon star U Antliae morphologically and physically in order to study the mass-loss evolution after a possible thermal pulse. Methods. High spatial resolution ALMA observations of unprecedented quality in thermal CO lines allow us to derive first critical spatial and temporal scales and constrain modeling efforts to estimate mass-loss rates for both the present day as well as the ejection period of the detached shell. Results. The detached shell is remarkably thin, overall spherically symmetric, and shows a barely resolved filamentary substructure possibly caused by instabilities in the interaction zone of winds with different outflow velocities. The expansion age of the detached shell is of the order of 2700 years and its overall width indicates a high expansion-velocity and high mass-loss period of only a few hundred years at an average mass-loss rate of $\approx$ 10$^{-5}$ $M_\odot$ yr$^{-1}$. The post-high-mass-loss-rate-epoch evolution of U Ant shows a significant decline to a substantially lower gas expansion velocity and a mass-loss rate amounting to 4$\times$10$^{-8}$ $M_\odot$ yr$^{-1}$, at present being consistent with evolutionary changes as predicted for the period between thermal pulses.

In September 2016, the microquasar Cygnus X-3 underwent a giant radio flare, which was monitored for 6 days with the Medicina Radio Astronomical Station and the Sardinia Radio Telescope. Long observations were performed in order to follow the evolution of the flare on a hourly scale, covering six frequency ranges from 1.5 GHz to 25.6 GHz. The radio emission reached a maximum of 13.2 +/- 0.7 Jy at 7.2 GHz and 10 +/- 1 Jy at 18.6 GHz. Rapid flux variations were observed at high radio frequencies at the peak of the flare, together with rapid evolution of the spectral index: alpha steepened from 0.3 to 0.6 within 5 hours. This is the first time that such fast variations are observed, giving support to the evolution from optically thick to optically thin plasmons in expansion moving outward from the core. Based on the Italian network (Noto, Medicina and SRT) and extended to the European antennas (Torun, Yebes, Onsala), VLBI observations were triggered at 22 GHz on five different occasions, four times prior to the giant flare, and once during its decay phase. Flux variations of 2-hour duration were recorded during the first session. They correspond to a mini-flare that occurred close to the core ten days before the onset of the giant flare. From the latest VLBI observation we infer that four days after the flare peak the jet emission was extended over 30 mas.

We present near-infrared interferometry of the carbon-rich asymptotic giant branch (AGB) star R Sculptoris. The visibility data indicate a broadly circular resolved stellar disk with a complex substructure. The observed AMBER squared visibility values show drops at the positions of CO and CN bands, indicating that these lines form in extended layers above the photosphere. The AMBER visibility values are best fit by a model without a wind. The PIONIER data are consistent with the same model. We obtain a Rosseland angular diameter of 8.9+-0.3 mas, corresponding to a Rosseland radius of 355+-55 Rsun, an effective temperature of 2640+-80 K, and a luminosity of log L/Lsun=3.74+-0.18. These parameters match evolutionary tracks of initial mass 1.5+-0.5 Msun and current mass 1.3+-0.7 Msun. The reconstructed PIONIER images exhibit a complex structure within the stellar disk including a dominant bright spot located at the western part of the stellar disk. The spot has an H-band peak intensity of 40% to 60% above the average intensity of the limb-darkening-corrected stellar disk. The contrast between the minimum and maximum intensity on the stellar disk is about 1:2.5. Our observations are broadly consistent with predictions by dynamic atmosphere and wind models, although models with wind appear to have a circumstellar envelope that is too extended compared to our observations. The detected complex structure within the stellar disk is most likely caused by giant convection cells, resulting in large-scale shock fronts, and their effects on clumpy molecule and dust formation seen against the photosphere at distances of 2-3 stellar radii.

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.

The closest ever fly-by of the Martian moon Phobos, performed by the European Space Agency's Mars Express spacecraft, gives a unique opportunity to sharpen and test the Planetary Radio Interferometry and Doppler Experiments (PRIDE) technique in the interest of studying planet - satellite systems. The aim of this work is to demonstrate a technique of providing high precision positional and Doppler measurements of planetary spacecraft using the Mars Express spacecraft. The technique will be used in the framework of Planetary Radio Interferometry and Doppler Experiments in various planetary missions, in particular in fly-by mode. We advanced a novel approach to spacecraft data processing using the techniques of Doppler and phase-referenced very long baseline interferometry spacecraft tracking. We achieved, on average, mHz precision (30 {\mu}m/s at a 10 seconds integration time) for radial three-way Doppler estimates and sub-nanoradian precision for lateral position measurements, which in a linear measure (at a distance of 1.4 AU) corresponds to ~50 m.

Supermassive black holes (SMBH) are essential for the production of jets in radio-loud active galactic nuclei (AGN). Theoretical models based on Blandford & Znajek extract the rotational energy from a Kerr black hole, which could be the case for NGC1052, to launch these jets. This requires magnetic fields of the order of $10^3\,$G to $10^4\,$G. We imaged the vicinity of the SMBH of the AGN NGC1052 with the Global Millimetre VLBI Array and found a bright and compact central feature, smaller than 1.9 light days (100 Schwarzschild radii) in radius. Interpreting this as a blend of the unresolved jet bases, we derive the magnetic field at 1 Schwarzschild radius to lie between 200 G and ~80000 G consistent with Blandford & Znajek models.

A test of a cornerstone of general relativity, the gravitational redshift effect, is currently being conducted with the RadioAstron spacecraft, which is on a highly eccentric orbit around Earth. Using ground radio telescopes to record the spacecraft signal, synchronized to its ultra-stable on-board H-maser, we can probe the varying flow of time on board with unprecedented accuracy. The observations performed so far, currently being analyzed, have already allowed us to measure the effect with a relative accuracy of $4\times10^{-4}$. We expect to reach $2.5\times10^{-5}$ with additional observations in 2016, an improvement of almost a magnitude over the 40-year old result of the GP-A mission.

Thermal pulses are fundamental to the chemical evolution of AGB stars and their circumstellar envelopes. A further consequence of thermal pulses is the formation of detached shells of gas and dust around the star. We aim to determine the physical properties of the detached gas shell around R Sculptoris, in particular the shell mass and temperature, and to constrain the evolution of the mass-loss rate during and after a thermal pulse. We analyse CO(1-0), CO(2-1), and CO(3-2) emission, observed by. The spatial resolution of the ALMA data allows us to separate the detached shell emission from the extended emission inside the shell. We perform radiative transfer modelling of both components to determine the shell properties and the post-pulse mass-loss properties. The ALMA data show a gas shell with a radius of 19.5" expanding at 14.3km/s. The different scales probed by the ALMA Cycle 0 array show that the shell must be entirely filled with gas, contrary to the idea of a detached shell. The comparison to single-dish spectra and radiative transfer modelling confirms this. We derive a shell mass of 4.5e-3 Msun with a temperature of 50K. Typical timescales for thermal pulses imply a pulse mass-loss rate of 2.3e-5 Msun/yr. For the post-pulse mass-loss rate, we find evidence for a gradual decline of the mass-loss rate, with an average value of 1.6e-5 Msun/yr. The total amount of mass lost since the last thermal pulse is 0.03 Msun, a factor four higher compared to classical models, with a sharp decline in mass-loss rate immediately after the pulse. We find that the mass-loss rate after a thermal pulse has to decline more slowly than generally expected from models of thermal pulses. This may cause the star to lose significantly more mass during a thermal pulse cycle, which affects the chemical evolution of the star and the interstellar medium.

Stars on the asymptotic giant branch (AGB) lose substantial amounts of matter, to the extent that they are important for the chemical evolution of, and dust production in, the universe. The mass loss is believed to increase gradually with age on the AGB, but it may also occur in the form of bursts, possibly related to the thermal pulsing phenomenon. Detached, geometrically thin, CO shells around carbon stars are good signposts of brief and intense mass ejection. We aim to put further constraints on the physical properties of detached CO shells around AGB stars. The photodissociation of CO and other carbon-bearing species in the shells leads to the possibility of detecting lines from neutral carbon. We have therefore searched for the CI($^3P_1-\,^3P_0$) line at 492 GHz towards two carbon stars, S Sct and R Scl, with detached CO shells of different ages, about 8000 and 2300 years, respectively. The CI($^3P_1-\,^3P_0$) line was detected towards R Scl. The line intensity is dominated by emission from the detached shell. The detection is at a level consistent with the neutral carbon coming from the full photodissociation of all species except CO, and with only limited photoionisation of carbon. The best fit to the observed $^{12}$CO and $^{13}$CO line intensities, assuming a homogeneous shell, is obtained for a shell mass of about 0.002 $M_\odot$, a temperature of about 100 K, and a CO abundance with respect to H$_2$ of 10$^{-3}$. The estimated CI/CO abundance ratio is about 0.3 for the best-fit model. However, a number of arguments point in the direction of a clumpy medium, and a viable interpretation of the data within such a context is provided.

We report about a 230 GHz very long baseline interferometry (VLBI) fringe finder observation of blazar 3C 279 with the APEX telescope in Chile, the phased submillimeter array (SMA), and the SMT of the Arizona Radio Observatory (ARO). We installed VLBI equipment and measured the APEX station position to 1 cm accuracy (1 sigma). We then observed 3C 279 on 2012 May 7 in a 5 hour 230 GHz VLBI track with baseline lengths of 2800 M$\lambda$ to 7200 M$\lambda$ and a finest fringe spacing of 28.6 micro-arcseconds. Fringes were detected on all baselines with SNRs of 12 to 55 in 420 s. The correlated flux density on the longest baseline was ~0.3 Jy/beam, out of a total flux density of 19.8 Jy. Visibility data suggest an emission region <38 uas in size, and at least two components, possibly polarized. We find a lower limit of the brightness temperature of the inner jet region of about 10^10 K. Lastly, we find an upper limit of 20% on the linear polarization fraction at a fringe spacing of ~38 uas. With APEX the angular resolution of 230 GHz VLBI improves to 28.6 uas. This allows one to resolve the last-photon ring around the Galactic Center black hole event horizon, expected to be 40 uas in diameter, and probe radio jet launching at unprecedented resolution, down to a few gravitational radii in galaxies like M 87. To probe the structure in the inner parsecs of 3C 279 in detail, follow-up observations with APEX and five other mm-VLBI stations have been conducted (March 2013) and are being analyzed.

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

We study the circumstellar evolution of the binary HD101584, consisting of a post-AGB star and a low-mass companion, which is most likely a post-common-envelope-evolution system. We used ALMA observations of the 12CO, 13CO, and C18O J=2-1 lines and the 1.3mm continuum to determine the morphology, kinematics, masses, and energetics of the circumstellar environment. The circumstellar medium has a bipolar hour-glass structure, seen almost pole-on, formed by an energetic jet, about 150 km/s. We conjecture that the circumstellar morphology is related to an event that took place about 500 year ago, possibly a capture event where the companion spiraled in towards the AGB star. However, the kinetic energy of the accelerated gas exceeds the released orbital energy, and, taking into account the expected energy transfer efficiency of the process, the observed phenomenon does not match current common-envelope scenarios. This suggests that another process must augment, or even dominate, the ejection process. A significant amount of material resides in an unresolved region, presumably in the equatorial plane of the binary system.

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 have mapped the CO(3-2) line emission around the Mira AB system at 0.5 resolution using the Atacama Large Millimeter/submillimeter Array (ALMA). The CO map shows amazing complexity. The circumstellar gas has been shaped by different dynamical actors during the evolution of the system and several morphological components can be identified. The companion is marginally resolved in continuum emission and is currently at 0.487$\pm$0.006 separation. In the main line component, centered on the stellar velocity, spiral arcs around Mira A are found. The spiral appears to be relatively flat and oriented in the orbital plane. An accretion wake behind the companion is clearly visible and the projected arc separation is of order 5''. In the blue wing of the line emission, offset from the main line, several large ($\sim$5-10''), opposing arcs are found. We tentatively suggest that this structure is created by the wind of Mira B blowing a bubble in the expanding envelope of Mira A.

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.

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 v=1 and v=2 J=1-0 (43 GHz), and v=1 J=2-1 (86 GHz) SiO masers are intense in AGB stars and have been mapped using VLBI showing ring-like distributions. Those of the v=1, v=2 J=1-0 masers are similar, but the spots are rarely coincident, while the v=1 J=2-1 maser arises from a well separated region farther out. These relative locations can be explained by models tools that include the overlap of two IR lines of SiO and H2O. The v=3 J=1-0 line is not directly affected by any line overlap and its spot structure and position, relative to the other lines, is a good test to the standard pumping models. We present single-dish and simultaneous VLBI observations of the v=1, v=2, and v=3 J=1-0 maser transitions of 28SiO in several AGB stars. The spatial distribution of the SiO maser emission in the v=3 J=1-0 transition from AGB stars is systematically composed of a series of spots that occupy a ring-like structure. The overall ring structure is extremely similar to that found in the other 43 GHz transitions and is very different from the structure of the v=1 J=2-1 maser. The positions of the individual spots of the different 43 GHz lines are, however, very rarely coincident, which in general is separated by about 0.3 AU (between 1 and 5 mas). These results are very difficult to reconcile with standard pumping models, which predict that the masers of rotational transitions within a given vibrational state require very similar excitation conditions, while the transitions of different vibrational states should appear in different positions. However, models including line overlap tend to predict v=1, v=2, v=3 J=1-0 population inversion to occur under very similar conditions, while the requirements for v=1 J=2-1 appear clearly different, and are compatible with the observational results.

[abridged] The 12CO/13CO ratio is often used as a measure of the 12C/13C ratio in the circumstellar environment, carrying important information about the stellar nucleosynthesis. External processes can change the 12CO and 13CO abundances, and spatially resolved studies of the 12CO/13CO ratio are needed to quantify the effect of these processes on the globally determined values. Additionally, such studies provide important information on the conditions in the circumstellar environment. The detached-shell source R Scl, displaying CO emission from recent mass loss, in a binary-induced spiral structure as well as in a clumpy shell produced during a thermal pulse, provides a unique laboratory for studying the differences in CO isotope abundances throughout its recent evolution. We observed both the 12CO(J=3-2) and the 13CO(J=3-2) line using ALMA. We find significant variations in the 12CO/13CO intensity ratios and consequently in the abundance ratios. The average CO isotope abundance ratio is at least a factor three lower in the shell (~19) than that in the present-day (<300 years) mass loss (>60). Additionally, variations in the ratio of more than an order of magnitude are found in the shell itself. We attribute these variations to the competition between selective dissociation and isotope fractionation in the shell, of which large parts cannot be warmer than ~35 K. However, we also find that the 12CO/13CO ratio in the present-day mass loss is significantly higher than the 12C/13C ratio determined in the stellar photosphere from molecular tracers (~19). The origin of this discrepancy is still unclear, but we speculate that it is due to an embedded source of UV-radiation that is primarily photo-dissociating 13CO. This radiation source could be the hitherto hidden companion. Alternatively, the UV-radiation could originate from an active chromosphere of R Scl itself....

A multi-transition survey of HCN (sub-) millimeter line emission from a large sample of AGB stars of different chemical type is presented. The data are analysed and circumstellar HCN abundances are estimated. The sample stars span a large range of properties such as mass-loss rate and photospheric C/O-ratio. The analysis of the new data allows for more accurate estimates of the circumstellar HCN abundances and puts new constraints on chemical models. In order to constrain the circumstellar HCN abundance distribution a detailed non-LTE excitation analysis, based on the Monte Carlo method, is performed. Effects of line overlaps and radiative excitation from dust grains are included. The median values for the derived abundances of HCN (with respect to H2) are 3x10-5, 7x10-7 and 10-7 for carbon stars (25 stars), S-type AGB stars (19 stars) and M-type AGB stars (25 stars), respectively. The estimated sizes of the HCN envelopes are similar to those obtained in the case of SiO for the same sample of sources and agree well with previous results from interferometric observations, when these are available. We find that there is a clear dependence of the derived circumstellar HCN abundance on the C/O-ratio of the star, in that carbon stars have about two orders of magnitude higher abundances than M-type AGB stars, on average. The derived HCN abundances of the S-type AGB stars have a larger spread and typically fall in between those of the two other types, however, slightly closer to the values for the M-type AGB stars. For the M-type stars, the estimated abundances are much higher than what would be expected if HCN is formed in thermal equilibrium. However, the results are also in contrast to predictions from recent non-LTE chemical models, where very little difference is expected in the HCN abundances between the various types of AGB stars.

We present the results of SiO maser observations at 43GHz toward two AGB stars using the VLBA. Our preliminary results on the relative positions of the different J=1-0 SiO masers (v=1,2 and 3) indicate that the current ideas on SiO maser pumping could be wrong at some fundamental level. A deep revision of the SiO pumping models could be necessary.

The asymptotic giant branch star R Sculptoris is surrounded by a detached shell of dust and gas. The shell originates from a thermal pulse during which the star undergoes a brief period of increased mass loss. It has hitherto been impossible to constrain observationally the timescales and mass-loss properties during and after a thermal pulse - parameters that determine the lifetime on the asymptotic giant branch and the amount of elements returned by the star. Here we report observations of CO emission from the circumstellar envelope and shell around R Sculptoris with an angular resolution of 1.3 arcsec. What was hitherto thought to be only a thin, spherical shell with a clumpy structure, is revealed to contain a spiral structure. Spiral structures associated with circumstellar envelopes have been seen previously, from which it was concluded that the systems must be binaries. Using the data, combined with hydrodynamic simulations, we conclude that R Sculptoris is a binary system that underwent a thermal pulse approximately 1800 years ago, lasting approximately 200 years. About 0.003 Msun of mass was ejected at a velocity of 14.3 km s-1 and at a rate approximately 30 times higher than the prepulse mass-loss rate. This shows that approximately 3 times more mass is returned to the interstellar medium during and immediately after a pulse than previously thought.

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.

New SiS multi-transition (sub-)millimetre line observations of a sample of AGB stars with varying photospheric C/O-ratios and mass-loss rates are presented. A combination of low- and high-energy lines are important in constraining the circumstellar distribution of SiS molecules. A detailed radiative transfer modelling of the observed SiS line emission is performed, including the effect of thermal dust grains in the excitation analysis. We find that the circumstellar fractional abundance of SiS in these environments has a strong dependence on the photospheric C/O-ratio as expected from chemical models. The carbon stars (C/O>1) have a mean fractional abundance of 3.1E-6, about an order of magnitude higher than found for the M-type AGB stars (C/O<1) where the mean value is 2.7E-7. These numbers are in reasonable agreement with photospheric LTE chemical models. SiS appears to behave similar to SiO in terms of photodissociation in the outer part of the circumstellar envelope. In contrast to previous results for the related molecule SiO, there is no strong correlation of the fractional abundance with density in the CSE, as would be the case if freeze-out onto dust grains were important. However, possible time-variability of the line emission in the lower J transitions and the sensitivity of the line emission to abundance gradients in the inner part of the CSE may mask a correlation with the density of the wind. There are indications that the SiS fractional abundance could be significantly higher closer to the star which, at least in the case of M-type AGB stars, would require non-equilibrium chemical processes.

We have performed high spatial resolution observations of SiO line emission for a sample of 11 AGB stars using the ATCA, VLA and SMA interferometers. Detailed radiative transfer modelling suggests that there are steep chemical gradients of SiO in their circumstellar envelopes. The emerging picture is one where the radial SiO abundance distribution starts at an initial high abundance, in the case of M-stars consistent with LTE chemistry, that drastically decreases at a radius of ~1E15 cm. This is consistent with a scenario where SiO freezes out onto dust grains. The region of the wind with low abundance is much more extended, typically ~1E16 cm, and limited by photodissociation. The surpisingly high SiO abundances found in carbon stars requires non-equilibrium chemical processes.

The nature of the mechanism responsible for producing the spectacular, geometrically thin, spherical shells found around some carbon stars has been an enigma for some time. Based on extensive radiative transfer modelling of both CO line emission and dust continuum radiation for all objects with known detached molecular shells, we present compelling evidence that these shells show clear signs of interaction with a surrounding medium. The derived masses of the shells increase with radial distance from the central star while their velocities decrease. A simple model for interacting winds indicate that the mass-loss rate producing the faster moving wind has to be almost two orders of magnitudes higher (~10^-5 solar masses per year) than the slower AGB wind (a few 10^-7 solar masses per year) preceding this violent event. At the same time, the present-day mass-loss rates are very low indicating that the epoch of high mass-loss rate was relatively short, on the order of a few hundred years. This, together with the number of sources exhibiting this phenomenon, suggests a connection with He-shell flashes (thermal pulses). We report the detection of a detached molecular shell around the carbon star DR Ser, as revealed from new single-dish CO (sub-)millimetre line observations. The properties of the shell are similar to those characterising the young shell around U Cam.

High angular resolution Australia Telescope Compact Array (ATCA) observations of SiO `thermal' millimetre line emission towards the two oxygen-rich, low mass loss rate AGB stars R Dor and L2 Pup are presented. In both cases the emission is resolved with an overall spherical symmetry. Detailed radiative transfer modelling of the SiO line emission has been performed, and the comparison between observations and models are conducted in the visibility plane, maximizing the sensitivity. The excitation analysis suggests that the abundance of SiO is as high as 4E-5 in the inner part of the wind, close to the predicted values from stellar atmosphere models. Beyond a radius of ~ 1E15 cm the SiO abundance is significantly lower, about 3E-6, until it decreases strongly at a radius of about 3E15 cm. This is consistent with a scenario where SiO first freezes out onto dust grains, and then eventually becomes photodissociated by the interstellar UV-radiation field. In these low expansion velocity sources the turbulent broadening of the lines plays an important role in the line formation. Micro-turbulent velocity widths in the range 1.1-1.5 km/s result in a very good reproduction of the observed line shapes even if the gas expansion velocity is kept constant. This, combined with the fact that the SiO and CO lines are well fitted using the same gas expansion velocity (to within 5-10%), suggest that the envelope acceleration occurs close to the stellar photosphere, within 20-30 stellar radii.

We present radio-interferometric observations of HCN J=1-0 line emission from the carbon star R Scl, obtained with the interim 3-mm receivers of the Australia Telescope Compact Array. The emission is resolved into a central source with a Gaussian FWHM of ~1 arcsec, which we identify as the present mass loss envelope. Using a simple photodissociation model and constraints from single-dish HCN spectra, we argue that the present mass-loss rate is low, ~2 x 10^{-7} Msol/yr, supporting the idea that R Scl had to experience a brief episode of intense mass loss in order to produce the detached CO shell at ~10 arcsec radius inferred from single-dish observations. Detailed radiative transfer modelling yields an abundance of HCN relative to H_2, of ~10^{-5} in the present-day wind. There appears to be a discrepancy between model results obtained with higher transition single-dish data included and those from the J=1-0 interferometer data alone, in that the interferometer data suggest a smaller envelope size and larger HCN abundance than the single-dish data. The lack of HCN in the detached shell is consistent with the rapid photodissociation of HCN into CN as it expands away from the star.

An extensive radiative transfer analysis of circumstellar SiO `thermal' radio line emission from a large sample of M-type AGB stars has been performed. The sample contains 18 irregulars of type Lb (IRV), 7 and 34 semiregulars of type SRa and SRb (SRV), respectively, and 12 Miras. New observational material, which contains data in several SiO rotational lines in the ground vibrational state, is presented for the IRV/SRVs. The detection rate was about 60%. SiO fractional abundances have been determined through radiative transfer modelling. The abundance distribution of the IRV/SRV sample has a median value of 8E-6. The high mass loss rate Miras have a bimodal abundance distribution, a low abundance group (on average 4E-7) and a high abundance group (on average 7E-6). The derived SiO abundances are in all cases well below the abundance expected from stellar atmosphere chemistry, on average by a factor of ten. In addition, there is a trend of decreasing SiO abundance with increasing mass loss rate. This is interpreted in terms of depletion of SiO molecules by the formation of silicate grains in the circumstellar envelopes, with an efficiency which increases with the mass loss rate. The SiO and CO radio line profiles differ in shape. In general, the dominating part of the SiO profiles are narrower than the CO profiles, but they have low-intensity wings which cover the full velocity range of the CO profile. This is interpreted (as has been done also by others) as due to the influence of gas acceleration in the region which produces a significant fraction of the SiO line emission. Finally, a number of sources which have peculiar CO line profiles are discussed from the point of view of their SiO line properties.