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.

The RadioAstron space radio telescope provides a unique opportunity to study the extreme brightness temperatures ($\mathrm{T_B }$) in AGNs with unprecedented long baselines of up to 28 Earth diameters. Since interstellar scintillation (ISS) may affect the visibilities observed with space VLBI (sVLBI), a complementary ground based flux density monitoring of the RadioAstron targets, which is performed near in time to the VLBI observation, could be beneficial. The combination/comparison with the sVLBI data can help to unravel the relative influence of source intrinsic and ISS induced effects, which in the end may alter the conclusions on the $\mathrm{T_B }$ measurements from sVLBI. Since 2013, a dedicated monitoring program has been ongoing to observe the ISS of RadioAstron AGN targets with a number of radio telescopes. Here we briefly introduce the program and present results from the statistical analysis of the Effelsberg monitoring data. We discuss the possible effects of ISS on $\mathrm{T_B }$ measurements for the RadioAstron target B0529+483 as a case study.

LS I +61$^{\circ}$303 is an X-ray binary with a radio outburst every ~27 days. Previous studies of the stellar system revealed radio microflares superimposed on the large radio outburst. We present here new radio observations of LS I +61$^{\circ}$303 at 2.2 GHz with the Westerbork Synthesis Radio Telescope (WSRT). Using various timing analysis methods we find significant Quasi-Periodic Oscillations (QPOs) of 55 minutes stable over the duration of 4 days. We also use archival data obtained from the Suzaku satellite at X-ray wavelengths. We report here for the first time significant X-ray QPOs of about 2 hours present over the time span of 21 hours. We compare our results with the previously reported QPO observations and we conclude that the QPOs seem to be associated with the radio outburst, independent of the amplitude of the outburst. Finally, the different QPO time-scales are discussed in the context of magnetic reconnection.

The high brightness temperatures, $T_\mathrm{b}\gtrsim 10^{13}$ K, detected in several active galactic nuclei by RadioAstron space VLBI observations challenge theoretical limits. Refractive scattering by the interstellar medium may affect such measurements. We quantify the scattering properties and the sub-mas scale source parameters for the quasar B0529+483. Using RadioAstron correlated flux density measurements at 1.7, 4.8, and 22 GHz on projected baselines up to 240,000 km we find two characteristic angular scales in the quasar core, about 100 $\mu$as and 10 $\mu$as. Some indications of scattering substructure are found. Very high brightness temperatures, $T_\mathrm{b}\geq 10^{13}$ K, are estimated at 4.8 GHz and 22 GHz even taking into account the refractive scattering. Our findings suggest a clear dominance of the particle energy density over the magnetic field energy density in the core of this quasar.

The JUICE mission will perform detailed measurements of the properties of the Galilean moons, with a nominal in-system science-mission duration of about 3.5 years. Using both the radio tracking data, and (Earth- and JUICE-based) optical astrometry, the dynamics of the Galilean moons will be measured to unprecedented accuracy. This will provide crucial input to the determination of the ephemerides and physical properties of the system, most notably the dissipation in Io and Jupiter. The data from Planetary Radio Interferometry and Doppler Experiment (PRIDE) will provide the lateral position of the spacecraft in the ICRF. In this article, we analyze the relative quantitative influence of the JUICE-PRIDE observables to the determination of the ephemerides of the Jovian system and the associated physical parameters. We perform a covariance analysis for a broad range of mission and system characteristics. We analyze the influence of VLBI data quality, observation planning, as well as the influence of JUICE orbit determination quality. This provides key input for the further development of the PRIDE observational planning and ground segment development. Our analysis indicates that the VLBI data are especially important for constraining the dynamics of Ganymede and Callisto perpendicular to their orbital planes. Also, the use of the VLBI data makes the uncertainty in the ephemerides less dependent on the error in the orbit determination of the JUICE spacecraft itself. Furthermore, we find that optical astrometry data of especially Io using the JANUS instrument will be crucial for stabilizing the solution of the normal equations. Knowledge of the dissipation in the Jupiter system cannot be improved using satellite dynamics obtained from JUICE data alone, the uncertainty in Io's dissipation obtained from our simulations is similar to the present level of uncertainty.

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.

The accuracy of orbit determination has a strong impact on the scientific output of the Space VLBI mission RadioAstron. The aim of this work is to improve the RadioAstron orbit reconstruction by means of sophisticated dynamical modelling of its motion in combination with multi-station Doppler tracking of the RadioAstron spacecraft. The improved orbital solution is demonstrated using Doppler measurements of the RadioAstron downlink signal and by correlating VLBI observations made by RadioAstron with ground-based telescopes using the enhanced orbit determination data. Orbit determination accuracy has been significantly improved from ~ 600 m in 3D position and ~ 2 cm/s in 3D velocity to several tens of metres and mm/s, respectively.

Aims. The phase scintillation of the European Space Agency's (ESA) Venus Express (VEX) spacecraft telemetry signal was observed at X-band (\lambda = 3.6 cm) with a number of radio telescopes of the European VLBI Network (EVN) in the period 2009-2013. Methods. We found a phase fluctuation spectrum along the Venus orbit with a nearly constant spectral index of -2.42 +/-0.25 over the full range of solar elongation angles from 0{\deg} to 45{\deg}, which is consistent with Kolmogorov turbulence. Radio astronomical observations of spacecraft signals within the solar system give a unique opportunity to study the temporal behaviour of the signal's phase fluctuations caused by its propagation through the interplanetary plasma and the Earth's ionosphere. This gives complementary data to the classical interplanetary scintillation (IPS) study based on observations of the flux variability of distant natural radio sources. Results. We present here our technique and the results on IPS. We compare these with the total electron content (TEC) for the line of sight through the solar wind. Finally, we evaluate the applicability of the presented technique to phase-referencing Very Long Baseline Interferometry (VLBI) and Doppler observations of currently operational and prospective space missions.

We present observations of a rare, rapid, high amplitude Extreme Scattering Event toward the compact BL-Lac AO 0235+164 at 6.65 GHz. The ESE cloud is compact; we estimate its diameter between 0.09 and 0.9 AU, and is at a distance of less than 3.6 kpc. Limits on the angular extent of the ESE cloud imply a minimum cloud electron density of ~ 4 x 10^3 cm^-3. Based on the amplitude and timescale of the ESE observed here, we suggest that at least one of the transients reported by Bower et al. (2007) may be attributed to ESEs.

The evidence for refractive interstellar scintillation (RISS) being the main cause for rapid intraday variations (Intraday Variability, IDV) in Quasars and BL Lacs has recently become stronger. If IDV is still a complex composition of extrinsic and source intrinsic effects, the intrinsic part of the IDV pattern should show up in the millimeter and sub-millimeter regime due to the frequency dependence of RISS. Hence, observations at higher frequencies are essential in order to exclude RISS as the sole cause of IDV. Here we report on our new attempt to search for rapid variations at much higher frequencies. In addition, the possibility of a direct detection of the postulated scattering screen in front of IDV sources will be discussed. Our recent line observations towards a few IDV sources lead to the first detection of a high latitude molecular cloud in front of an intraday variable radio core.

Intraday variability (IDV) of compact extragalactic radio sources is a complex phenomenon and shows a wavelength dependent mixture of refractive interstellar scintillation (RISS) (dominant at long cm-wavelengths) and source-intrinsic effects (dominant at shorter wavelengths). Detailed investigations of individual sources and new high frequency observations suggest a source-intrinsic contribution to the IDV pattern at least in some sources. However, the sizes of intraday variable sources at cm-wavelength are typically smaller than the scattering size set by the ISM in our galaxy and scintillation must be present, too. We present new IDV observations in different regimes (from cm to sub-mm wavelengths) and show how such multi-frequency study can be used as powerful instrument to describe different aspects of IDV.

There is new evidence which identifies seasonal changes of the variability time scale in intraday variable sources with refractive interstellar scintillation effects. Such a RISS-model takes the annual change of the Earth velocity relative to the scattering medium into account. In September 1998 we found a remarkable prolongation of the variability time scale in the IDV-source 0917+624 with only small variations in flux density during a period of 5 days. This was explained as a seasonal effect, in which the velocity vector of the Earth and the interstellar medium nearly canceled. In order to further investigate the applicability of the model for 0917+624, we performed an Effelsberg 6cm-flux monitoring program over the course of one year. Since September 2000, the source appears to be remarkably inactive and yet (May 2001), no return to its normal, faster and stronger variability pattern is observed. Here, our observational results and a possible explanation for the current quiescence are presented.

Extreme Scattering Events (ESEs) are dramatic variations of the flux density at Gigahertz frequencies caused by ray path distortions within an isolated inhomogeneity ("plasma lens") in the interstellar medium. These events are characterized by a deep flux density minimum in the light curve with, in some cases, surrounding maxima. The variability time scales range from weeks to months. These phenomena show a strong frequency dependence, in which the variability amplitudes increase with wavelength. During an Intraday Variability (IDV) monitoring project (March 2000), a feature resembling an ESE-like event appeared in the variable light curve of 0954+658, however with a time scale of less than two days. We will discuss this effect and its implications for a better description of the interstellar medium.