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.

VLBI observations of SiO masers recover at most 40-50% of the total flux obtained by single dish observations at any spectral channel. Some previous studies seems to indicate that, at least, part of the lost flux is divided up into many weak components rather than in a large resolved emission area. Taking benefit of the high sensitivity and resolution of the HSA, we investigate the problem of the missing flux in VLBI observations of SiO maser emission at 7mm in the AGB stars and obtain a high dynamic range map of IRC+10011. We conclude that the missing flux is mostly contained in many very weak maser components.

Observations of 28SiO v=0 J=1-0 line emission (7-mm wavelength) from AGB stars show in some cases peculiar profiles, composed of a central intense component plus a wider plateau. Very similar profiles have been observed in CO lines from some AGB stars and most post-AGB nebulae and, in these cases, they are clearly associated with the presence of conspicuous axial symmetry and bipolar dynamics. We present systematic observations of 28SiO v=0 J=1-0 emission in 28 evolved stars, performed with the 40~m radio telescope of the IGN in Yebes, Spain. We find that the composite core plus plateau profiles are almost always present in O-rich Miras, OH/IR stars, and red supergiants. They are also found in one S-type Mira ($\chi$ Cyg), as well as in two semiregular variables (X Her and RS Cnc) that are known to show axial symmetry. In the other objects, the profiles are simpler and similar to those of other molecular lines. The composite structure appears in the objects in which SiO emission is thought to come from the very inner circumstellar layers, prior to dust formation. The central spectral feature is found to be systematically composed of a number of narrow spikes, except for X Her and RS Cnc, in which it shows a smooth shape that is very similar to that observed in CO emission. These spikes show a significant (and mostly chaotic) time variation, while in all cases the smooth components remain constant within the uncertainties. The profile shape could come from the superposition of standard wide profiles and a group of weak maser spikes. Alternatively, we speculate that the very similar profiles detected in objects that are axisymmetric may be indicative of the systematic presence of a significant axial symmetry in the very inner circumstellar shells around AGB stars; the presence of such symmetry would be independent of the probable weak maser effects in the central spikes.

Cyanogen (NCCN) is the simplest member of the series of dicyanopolyynes. It has been hypothesized that this family of molecules can be important constituents of interstellar and circumstellar media, although the lack of a permanent electric dipole moment prevents its detection through radioastronomical techniques. Here we present the first solid evidence of the presence of cyanogen in interstellar clouds through the detection of its protonated form toward the cold dark clouds TMC-1 and L483. Protonated cyanogen (NCCNH+) has been identified through the J=5-4 and J=10-9 rotational transitions using the 40m radiotelescope of Yebes and the IRAM 30m telescope. We derive beam averaged column densities for NCCNH+ of (8.6+/-4.4)e10 cm-2 in TMC-1 and (3.9+/-1.8)e10 cm-2 in L483, which translate to fairly low fractional abundances relative to H2, in the range (1-10)e-12. The chemistry of protonated molecules in dark clouds is discussed, and it is found that, in general terms, the abundance ratio between the protonated and non protonated forms of a molecule increases with increasing proton affinity. Our chemical model predicts an abundance ratio NCCNH+/NCCN of 1e-4, which implies that the abundance of cyanogen in dark clouds could be as high as (1-10)e-8 relative to H2, i.e., comparable to that of other abundant nitriles such as HCN, HNC, and HC3N.

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.

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.

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

An international consortium is presently constructing a beamformer for the Atacama Large Millimeter/submillimeter Array (ALMA) in Chile that will be available as a facility instrument. The beamformer will aggregate the entire collecting area of the array into a single, very large aperture. The extraordinary sensitivity of phased ALMA, combined with the extremely fine angular resolution available on baselines to the Northern Hemisphere, will enable transformational new very long baseline interferometry (VLBI) observations in Bands 6 and 7 (1.3 and 0.8 mm) and provide substantial improvements to existing VLBI arrays in Bands 1 and 3 (7 and 3 mm). The ALMA beamformer will have impact on a variety of scientific topics, including accretion and outflow processes around black holes in active galactic nuclei (AGN), tests of general relativity near black holes, jet launch and collimation from AGN and microquasars, pulsar and magnetar emission processes, the chemical history of the universe and the evolution of fundamental constants across cosmic time, maser science, and astrometry.

There is growing evidence of relativistic jets in radio-loud narrow-line Seyfert 1 (RL-NLS1) galaxies. We constrain the observational properties of the radio emission in the first RL-NLS1 galaxy ever detected in gamma-rays, PMN J0948+0022, i.e., its flux density and structure in total intensity and in polarization, its compactness, and variability. We performed three real-time e-VLBI observations of PMN J0948+0022 at 22 GHz, using a global array including telescopes in Europe, East Asia, and Australia. These are the first e-VLBI science observations ever carried out with a global array, reaching a maximum baseline length of 12458 km. The observations were part of a large multiwavelength campaign in 2009. The source is detected at all three epochs. The structure is dominated by a bright component, more compact than 55 microarcsec, with a fainter component at a position angle theta~ 35deg. Relativistic beaming is required by the observed brightness temperature of 3.4x10^11 K. Polarization is detected at a level of about 1%. The parameters derived by the VLBI observations, in addition to the broad-band properties, confirm that PMN J0948+0022 is similar to flat spectrum radio quasars. Global e-VLBI is a reliable and promising technique for future studies.

The detection of gamma-ray emission by Fermi-LAT from the radio loud Narrow Line Seyfert 1 PMN J0948+0022 (Abdo et al. 2009, ApJ 699, 976) triggered a multi-wavelength campaign between March and July 2009. Given its high compactness (Doi et al. 2006, PASJ 58, 829), inverted spectrum, and 0deg declination, the source was an ideal target to observe at 22 GHz with a Global VLBI array extending from Europe to East Asia and Australia. In order to deliver prompt results to be analysed in combination with the other instruments participating in the campaign, the observations were carried out with real time VLBI, for the first time on a Global scale. Indeed, the main results have been published just a few months after the campaign (Abdo et al. 2009, ApJ 707, 727). Here we present additional details about the e-VLBI observations.

Following the recent discovery of gamma rays from the radio-loud narrow-line Seyfert 1 galaxy PMN J0948+0022 (z=0.5846), we started a multiwavelength campaign from radio to gamma rays, which was carried out between the end of March and the beginning of July 2009. The source displayed activity at all the observed wavelengths: a general decreasing trend from optical to gamma-ray frequencies was followed by an increase of radio emission after less than two months from the peak of the gamma-ray emission. The largest flux change, about a factor of about 4, occurred in the X-ray band. The smallest was at ultraviolet and near-infrared frequencies, where the rate of the detected photons dropped by a factor 1.6-1.9. At optical wavelengths, where the sampling rate was the highest, it was possible to observe day-scale variability, with flux variations up to a factor of about 3. The behavior of PMN J0948+0022 observed in this campaign and the calculated power carried out by its jet in the form of protons, electrons, radiation and magnetic field are quite similar to that of blazars, specifically of flat-spectrum radio quasars. These results confirm the idea that radio-loud narrow-line Seyfert 1 galaxies host relativistic jets with power similar to that of average blazars.

High resolution maps of maser emission provide very detailed information on processes occurring in circumstellar envelopes of late-type stars. A particularly detailed picture of the innermost shells around AGB stars is provided by SiO masers. Considerable progress is being made to provide astrometrically aligned multi-transition simultaneous observations of these masers, which are needed to better constrain the models. In view of the large amount of high quality data available, models should now be developed to fully explain all maser characteristics together (spatial distribution, variability, etc). New generation instruments (VERA, VSOP-2), new observational techniques (frequency-phase transfer), and new models promise important improvements of our knowledge on this topic.

Oxygen-rich Asymptotic Giant Branch (AGB) stars are intense emitters of SiO and H$_2$O maser lines at 43 (J=1-0, v=1 and 2) and 22 GHz, respectively. VLBI observations of the maser emission provides a unique tool to sample the innermost layers of the circumstellar envelopes in AGB stars. Nevertheless, the difficulties in achieving astrometrically aligned v=1 and v=2 SiO maser maps have traditionally prevented a unique interpretation of the observations in terms of physical underlying conditions, which depend on the nature of the SiO pumping mechanism. We have carried out observations of the SiO and H$_2$O maser emission towards RLMi, using the astrometric capabilities of VERA. Due to the too-weak emission of the reference calibrator we had to develop a special method to accurately relate the coordinates for both transitions. We present relative astrometrically aligned v=1 and v=2 J=1-0 SiO maser maps, at multiple epochs, and discuss the astrophysical results. The incorporation of astrometric information into the maps of SiO masers challenges the weak points in the current theoretical models, which will need further refinements to address the observations results.

OH 231.8+4.2 is a well studied preplanetary nebula (pPN) around a binary stellar system that shows a remarkable bipolar outflow. To study the structure and kinematics of the inner 10-80 AU nebular regions we performed high-resolution observations of the H$_2$O 6$_{1,6}$--5$_{2,3}$ and $^{28}$SiO $v$=2, $J$=1--0 maser emissions with the Very Long Baseline Array. The absolute position of both emission distributions were recovered using the phase referencing technique, and accurately registered in HST optical images. H$_2$O maser clumps are found to be distributed in two areas of 20 mas in size spatially displaced by $\sim$60 milli-arcseconds along an axis oriented nearly north-south. SiO masers are tentatively found to be placed between the two H$_2$O maser emitting regions, probably indicating the position of the Mira component of the system. The SiO maser emission traces an inner equatorial component with a diameter of 12 AU, probably a disk rotating around the M-type star. Outwards, we detect in the H$_2$O data a pair of polar caps, separated by 80 AU. We believe that the inner regions of the nebula probably have been altered by the presence of the companion, leading to an equator-to-pole density contrast that may explain the lack of H$_2$O masers and strong SiO maser emission in the denser, equatorial regions.

The study of the innermost circumstellar layers around AGB stars is crucial to understand how these envelopes are formed and evolve. The SiO maser emission occurs at a few stellar radii from the central star, providing direct information on the stellar pulsation and on the chemical and physical properties of these regions. Our data also shed light on several aspects of the SiO maser pumping theory that are not well understood yet. We aim to determine} the relative spatial distribution of the 43 GHz and 86 GHz SiO maser lines in the oxygen-rich evolved star R Leo. We have imaged with milliarcsecond resolution, by means of Very Long Baseline Interferometry, the 43 GHz (28SiO v=1, 2 J=1-0 and 29SiO v=0 J=1-0) and 86 GHz (28SiO v=1 J=2-1 and 29SiO v=0 J=2-1) masing regions. We confirm previous results obtained in other oxygen-rich envelopes. In particular, when comparing the 43 GHz emitting regions, the 28SiO v=2 transition is produced in an inner layer, closer to the central star. On the other hand, the 86 GHz line arises in a clearly farther shell. We have also mapped for the first time the 29SiO v=0 J=1-0 emission in R Leo. The already reported discrepancy between the observed distributions of the different maser lines and the theoretical predictions is also found in R Leo.

VLBI is unique, among the space geodetic techniques, in its contribution to defining and maintaining the International Celestial Reference Frame, providing precise measurements of coordinates of extragalactic radiosources. The quest for increasing accuracy of VLBI geodetic products has lead to a deeper revision of all aspects that might introduce errors in the analysis. The departure of the observed sources from perfect, stable, compact and achromatic celestial targets falls within this category. This paper is concerned with the impact of unaccounted frequency-dependent position shifts of source cores in the analysis of dual-band S/X VLBI geodesy observations, and proposes a new method to measure them. The multi-frequency phase transfer technique developed and demonstrated by Middelberg et al. (2005) increases the high frequency coherence times of VLBI observations, using the observations at a lower frequency. Our proposed SOURCE/FREQUENCY PHASE REFERENCING method endows it with astrometric applications by adding a strategy to estimate the ionospheric contributions. Here we report on the first successful application to measure the core shift of the quasar 1038+528 A at S and X-bands, and validate the results by comparison with those from standard phase referencing techniques. In this particular case, and in general in the cm-wavelength regime, both methods are equivalent. Moreover the proposed method opens a new horizon with targets and fields suitable for high precision astrometric studies with VLBI, especially at high frequencies where severe limitations imposed by the rapid fluctuations in the troposphere prevent the use of standard phase referencing techniques.

We report the first VLBI map of the v=0 J=1-0 maser line of 29SiO in the long-period variable star IRC+10011. We have found that this maser emission is composed of multiple spots distributed in an incomplete ring, suggesting that this maser is also amplified tangentially, as already proposed in other SiO circumstellar masers. We present also VLBI maps for the 7 mm 28SiO v=1 and 2 J=1-0 and the 3 mm v=1 J=2-1 lines. The 29SiO masing region appears to be located in a layer in between the 28SiO v=1 J=1-0 and 28SiO v=1 J=2-1 lines. In addition, we confirm that the 86 GHz maser v=1 J=2-1 forms in an outer region of the circumstellar envelope compared to the other 28SiO masers studied. Finally, we discuss the possible implications of the observational results on the SiO maser pumping theory.

We present milliarsecond-resolution maps of the SiO maser emission v=1 J=2-1 (3 mm) in the bipolar post-AGB nebula OH231.8+4.2, and compare them with our previous observations of the v=2 J=1-0 line (7 mm). Our observations show that the SiO masers arise in several bright spots forming a structure elongated in a direction perpendicular to the symmetry axis of the nebula. This, and the complex velocity gradient observed, is consistent with the presence of an equatorial torus in rotation and with an infall of material towards the star.

We report the first VLBI map of the 7 mm v=0 J=1-0 maser line of the 29-SiO isotopic substitution in the long period variable star IRC+10011. We have found that this maser emission is composed of multiple features distributed in an incomplete ring, suggesting that this maser is also amplified tangentially as already proposed in other SiO circumstellar masers. We present also the results for some 28-SiO 7 mm and 3 mm maser lines. We confirm in this epoch that the 86 GHz maser v=1 J=2-1 in IRC+10011 is formed in an outer region of the circumstellar envelope compared to the other SiO masers studied. The 29-SiO masing region appears to be located in a layer in between the 28-SiO v=1 J=1-0 and 28-SiO v=1 J=2-1 lines. Finally, we discuss the possible implications of the observational results in the SiO maser pumping theory.

We present sub-milliarcsecond observations of SiO masers in the late-type stars IRC +10011 and Chi Cyg. We have used the NRAO Very Long Baseline Array (VLBA) to map the 43 GHz (v=1, 2 J=1-0) and the 86 GHz (v=1, 2 J=2-1) SiO masers. All the transitions have been imaged except the v=2 J=2-1 in IRC +10011. We report the first VLBI map of the v=1 J=2-1 28SiO maser in IRC +10011 as well as the first VLBA images of SiO masers in an S-type Mira variable, Chi Cyg. In this paper we have focused on the study of the relative spatial distribution of the different observed lines. We have found that in some cases the observational results are not reproduced by the current theoretical pumping models, either radiative or collisional. In particular, for IRC +10011, the v=1 J=1-0 and J=2-1 28SiO lines have different spatial distributions and emitting region sizes, the J=2-1 emission being located in an outer region of the envelope. For Chi Cyg, the distributions also differ, but the sizes of the masing regions are comparable. We suggest that the line overlaps between ro-vibrational transitions of two abundant molecular species, H2O and 28SiO, is a possible explanation for the discrepancies found between the observations and the theoretical predictions. We have introduced this overlapping process in the calculations of the excitation of the SiO molecule. We conclude that the line overlaps can strongly affect the excitation of SiO and may reproduce the unexpected observational results for the two sources studied.

On this paper, we present simultaneous VLBA observations of the J=1-0 and J=2-1 rotational lines in the v=1 and v=2 vibrationally excited states of SiO (at 7 and 3 mm wavelengths). We have mapped these four maser lines in the circumstellar envelopes of three AGB stars: IRC +10011, Chi Cyg and TX Cam. We study the relative spatial distribution between these maser lines. In particular, for the two v=1 transitions, we found that the J=1-0 and J=2-1 maser spot distributions are unalike, challenging all the current theoretical pumping mechanisms for SiO masers (which predict quite similar distribution for maser lines in the same vibrational excited state).

We present milliarsecond-resolution maps of the SiO maser emission $v$=1 J=2--1 (3 mm) in the bipolar post-AGB nebula OH 231.8+4.2, and compare them with our previous observations of the $v$=2 J=1--0 line (7 mm). Our observations show that the SiO masers arise in several bright spots forming a structure elongated in a direction perpendicular to the symmetry axis of the nebula. This, and the complex velocity gradient observed, is consistent with the presence of an equatorial torus in rotation and with an infall of material towards the star.

We present the first VLBI images of SiO masers in the circumstellar envelope of an S-type star obtained with the newly developed capabilities of the Very Long Baseline Array (VLBA) at 86 GHz. We combine these data with those obtained quasy simultaneously at 43 GHz. These observations provide information on the structure and dynamics of the innermost circumstellar shells, where the return of large quantities of stellar material to the interstellar medium starts. Despite of the fact that these are prelimininary results, we report that the maser emission of the $v$=1 J=1--0 and J=2--1 present a very different emission distribution, being not coincident in most case.

We performed the first VLBI observations of the SiO v=1 and v=2 J=1--0 masers in a Proto Planetary Nebula OH 231.8+4.2 (also known as the Rotten Egg nebula), at milliarcsecond resolution. Only the v=2 maser transition was detected. We detect several maser spots lying along a line which is almost perpendicular to the axis of symmetry of the Nebula. We find that all the emission is concentrated in an area of ~10 mas, indicating that the SiO masers are originated very close to the surface of the star. One of the two emission areas presents an elongated structure with a clear velocity gradient. The detected emission is consistent with a torus, or disk, in rotation with a velocity of ~6 km/s and with an infall velocity of ~10 km/s.

We have performed VLBA observations of the SiO v=1 and v=2 J=1-0 masers in two AGB stars, TX Cam and IRC +10011. We confirm the ring-like spatial distribution, previously found in several AGB objects, as well as the tangential polarization pattern, already reported for TX Cam. Both properties, that seem to be systematic in this kind of objects, are characteristic of radiatively pumped SiO masers. On the contrary, we do not confirm the previous report on the spatial coincidence between the J=1-0 v=1 and 2 masers, a result that would have argued in favor of collisional pumping. We find that both lines sometimes arise from nearby spots, typically separated by 1-2 mas, but are rarely coincident. The discrepancy with previous results is explained by the very high spatial resolution of our observations, ~ 0.5 mas, an order of magnitude better than in the relevant previously published experiment.