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.

Context. Red supergiant stars possess surface features and extended molecular atmospheres. Photospheric convection may be a crucial factor of the levitation of the outer atmospheric layers. However, the mechanism responsible is still poorly understood. Aims. We image the stellar surface of V602 Carinae (V602 Car) to constrain the morphology and contrast of the surface features and of the extended atmospheric layers. Methods. We observed V602 Car with the Very Large Telescope Interferometer (VLTI) PIONIER instrument (1.53-1.78 $\mathrm{\mu}$m) between May and July 2016, and April and July 2019 with different telescope configurations. We compared the image reconstructions with 81 temporal snapshots of 3D radiative-hydrodynamics (RHD) CO$^5$BOLD simulations in terms of contrast and morphology, using the Structural Similarity Index. Results. The interferometric data are compatible with an overall spherical disk of angular diameter 4.4$\pm$0.2 mas, and an extended molecular layer. In 2016, the reconstructed image reveals a bright arc-like feature toward the northern rim of the photospheric surface. In 2019, an arc-like feature is seen at a different orientation and a new peak of emission is detected on the opposite side. The contrasts of the reconstructed surface images are 11%$\pm$2% and 9%$\pm$2% for 2016 and 2019, respectively. The morphology and contrast of the two images are consistent with 3D RHD simulations, within our achieved spatial resolution and dynamic range. The extended molecular layer contributes 10--13% of the total flux with an angular diameter of 6--8 mas. It is present but not clearly visible in the reconstructed images because it is close to the limits of the achieved dynamic range. The presence of the molecular layer is not reproduced by the 3D RHD simulations. Conclusions...

Studies of fundamental parameters of very low-mass objects are indispensable to provide tests of stellar evolution models that are used to derive theoretical masses of brown dwarfs and planets. However, only objects with dynamically determined masses and precise photometry can effectively evaluate the predictions of stellar models. AB Dor C (0.090 solar masses) has become a prime benchmark for calibration of theoretical evolutionary models of low-mass young stars. One of the ambiguities remaining in AB Dor C is the possible binary nature of this star. We observed AB Dor C with the VLTI/AMBER instrument in low-resolution mode at the J, H and K bands. The interferometric observables at the K-band are compatible with a binary brown dwarf system with tentative components AB Dor Ca/Cb with a K-band flux ratio of 5$\pm$1% and a separation of 38$\pm$1 mas. This implies theoretical masses of 0.072$\pm$0.013 M$_{\rm \odot}$ and 0.013$\pm$0.001 M$_{\rm \odot}$ for each component, near the hydrogen-burning limit for AB Dor Ca, and near the deuterium-burning limit, straddling the boundary between brown dwarfs and giant planets, for AB Dor Cb. The possible binarity of AB Dor C alleviates the disagreement between observed magnitudes and theoretical mass-luminosity relationships.

This paper discusses the science case for a sensitive spectro-polarimetric survey of the microwave sky. Such a survey would provide a tomographic and dynamic census of the three-dimensional distribution of hot gas, velocity flows, early metals, dust, and mass distribution in the entire Hubble volume, exploit CMB temperature and polarisation anisotropies down to fundamental limits, and track energy injection and absorption into the radiation background across cosmic times by measuring spectral distortions of the CMB blackbody emission. In addition to its exceptional capability for cosmology and fundamental physics, such a survey would provide an unprecedented view of microwave emissions at sub-arcminute to few-arcminute angular resolution in hundreds of frequency channels, a data set that would be of immense legacy value for many branches of astrophysics. We propose that this survey be carried-out with a large space mission featuring a broad-band polarised imager and a moderate resolution spectro-imager at the focus of a 3.5m aperture telescope actively cooled to about 8K, complemented with absolutely-calibrated Fourier Transform Spectrometer modules observing at degree-scale angular resolution in the 10-2000 GHz frequency range. We propose two observing modes: a survey mode to map the entire sky as well as a few selected wide fields, and an observatory mode for deeper observations of regions of specific interest.

Cool evolved stars are the main source of chemical enrichment of the interstellar medium. Understanding their mass loss offers a unique opportunity to study the cycle of matter. We discuss interferometric studies and their comparison to latest state-of-the-art dynamic model atmospheres. They show broad agreement for asymptotic giant branch stars. For red supergiants, however, current models cannot explain observed extensions by far, pointing to missing physical processes in their models, and uncertainties in our general understanding of mass loss. We present ongoing imaging and time-series observations that may provide the strongest constraint and may help to identify missing dynamic processes. VLTI studies will remain the highest spatial resolution observations at ESO into the ELT era, complemented by ALMA observations. We discuss crucial improvements in both instrumental and operational areas for surface imaging of cool evolved stars in the era of the ELT.

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.