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.

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

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

The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI). The Spanish SKA White Book is the result of the coordinated effort of 120 astronomers from 40 different research centers. The book shows the enormous scientific interest of the Spanish astronomical community in the SKA and warrants an optimum scientific exploitation of the SKA by Spanish researchers, if Spain enters the SKA project.

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

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