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.

RZ Psc is a young Sun-like star with a bright and warm infrared excess that is occasionally dimmed significantly by circumstellar dust structures. Optical depth arguments suggest that the dimming events do not probe a typical sight line through the circumstellar dust, and are instead caused by structures that appear above an optically thick mid-plane. This system may therefore be similar to systems where an outer disk is shadowed by material closer to the star. Here we report the discovery that RZ Psc hosts a 0.12$M_\odot$ companion at a projected separation of 23au. We conclude that the disk must orbit the primary star. While we do not detect orbital motion, comparison of the angle of linear polarization of the primary with the companion's on-sky position angle provides circumstantial evidence that the companion and disc may not share the same orbital plane. Whether the companion severely disrupts the disc, truncates it, or has little effect at all, will require further observations of both the companion and disc.

Serendipitously, we have rediscovered radio emission from the binary system Ross 867 (M4.5V) and Ross 868 (M3.5V) while inspecting archival Giant Metrewave Radio Telescope (GMRT) observations. The binary system consists of two M-dwarf stars that share common characteristics such as spectral type, astrometric parameters, age and emission at infrared, optical and X-rays frequencies. The GMRT data at 610 MHz taken on July 2011 shows that the radio emission from Ross 867 is polarized and highly variable on hour time scales with a peak flux of 10.4 $\pm$ 0.7 mJy/beam. Additionally, after reviewing archival data from several observatories (VLA, GMRT, JVLA and LOFAR), we confirm that although both stars are likely coeval, only Ross 867 has been detected, while Ross 868 remains undetected at radio wavelengths. As they have a a large orbital separation, this binary stellar system provides a coeval laboratory to examine and constrain the stellar properties linked to radio activity in M dwarfs. We speculate that the observed difference in radio activity between the dwarfs could be due to vastly different magnetic field topologies or that Ross 867 has an intrinsically different dynamo.

We present trigonometric parallax and proper motion measurements toward 22 GHz water and 6.7 GHz methanol masers in 16 high-mass star-forming regions. These sources are all located in the Scutum spiral arm of the Milky Way. The observations were conducted as part of the Bar and Spiral Structure Legacy (BeSSeL) survey. A combination of 14 sources from a forthcoming study and 14 sources from the literature, we now have a sample of 44 sources in the Scutum spiral arm, covering a Galactic longitude range from 0$^\circ$ to 33$^\circ$. A group of 16 sources shows large peculiar motions of which 13 are oriented toward the inner Galaxy. A likely explanation for these high peculiar motions is the combined gravitational potential of the spiral arm and the Galactic bar.

Sharpless 269 (S269) is one of a few HII regions in the outer spiral arm of the Milky Way with strong water maser emission. Based on data from the Very Long Baseline Interferometry (VLBI) Exploration of Radio Astrometry (VERA) array, two parallax measurements have been published, which differ by nearly $2\sigma$. Each distance estimate supports a different structure for the outer arm. Moreover, given its large Galactocentric radii, S269 has special relevance as its proper motion and parallax have been used to constrain the Galactic rotation curve at large radii. Using recent Very Long Baseline Array (VLBA) observations, we accurately measure the parallax and proper motion of the water masers in S269. We interpret the position and motion of S269 in the context of Galactic structure, and possible optical counterparts. S269's 22 GHz water masers and two close-by quasars were observed at 16 epochs between 2015 and 2016 using the VLBA. We measure an annual parallax for S269 of 0.241 $\pm$ 0.012 mas corresponding to a distance from the Sun of $4.15^{+0.22}_{-0.20}$ kpc by fitting four maser spots. The mean proper motion for S269 was estimated as $0.16\pm0.26$ mas $\rm{yr^{-1}}$ and $-0.51\pm0.26$ mas $\rm{yr^{-1}}$ for $\mu_{\alpha} \ cos \delta$ and $\mu_{\delta}$ respectively, which corresponds to the motion expected for a flat Galactic rotation curve at large radius. This distance estimate, Galactic kinematic simulations and observations of other massive young stars in the outer region support the existence of a kink in the outer arm at $l \approx$ 140 degrees. Additionally, we find more than 2,000 optical sources in the Gaia DR2 catalog within 125 pc radius around the 3D position of the water maser emission; from those only three sources are likely members of the same stellar association that contains the young massive star responsible for the maser emission (S269 IRS 2w).

With the publication of Gaia DR2, 1.3 billion stars now have public parallax and proper motion measurements. In this contribution, we compare the results for sources that have both optical and radio measurements, focusing on circumstellar masers. For these large, variable and bright AGB stars, the VLBI astrometry results can be more robust. Moreover, there are a number of applications where VLBI astrometry provides unique data for studying stellar populations and Galactic structure. The BeSSel project not only provides parallax and proper motions at much larger distances than Gaia can reach, but it also uniquely samples the spiral arms of the Galaxy. The evolved stars in the BAaDE sample can potentially constrain the dynamics and stellar content of the inner bulge and bar of the Milky Way, not reachable in the optical.

Observations of dust-enshrouded evolved stars selected from infrared catalogs requiring high positional accuracy, like infrared spectroscopy or long baseline radio interferometric observations, often require preparational observational steps determining a position with an accuracy much better than 1". Using phase-referencing observations with the Very Large Array at its highest resolution, we have compared the positions of SiO 43 GHz masers in evolved stars, assumed to originate in their infrared detected circumstellar shells, with the positions listed in the MSX, WISE, 2MASS, and Gaia catalogs. Starting from an MSX position it is, in general, simple to match 2MASS and WISE counterparts. However, in order to obtain a Gaia match to the MSX source it is required to use a two-step approach due to the large number of nearby candidates and low initial positional accuracy of the MSX data. We show that the closest comparable position to the SiO maser in our limited sample never is the MSX position. When a plausible source with a characteristic signature of an evolved star with a circumstellar shell can be found in the area, the best indicator of the maser position is provided by the Gaia position, with the 2MASS position being second-best. Typical positional offsets from all catalogs to the SiO masers are reported.

Radio observations of cool stellar systems provide unique information on their magnetic fields, high-energy processes, and chemistry. Buoyed by powerful new instruments (e.g. ALMA, JVLA, LOFAR), advances in related fields (e.g., the Gaia astrometric revolution), and above all a renewed interest in the relevant stellar astrophysics, stellar radio astronomy is experiencing a renaissance. In this splinter session, participants took stock of the present state of stellar radio astronomy to chart a course for the field's future.

Radio astrometric campaigns using VLBI have provided distances and proper motions for masers associated with young massive stars (BeSSeL survey). The ongoing BAaDE project plans to obtain astrometric information of SiO maser stars located in the inner Galaxy. These stars are associated with evolved, mass-losing stars. By overlapping optical (Gaia), infrared (2MASS, MSX and WISE) and radio (BAaDE) sources, we expect to obtain important clues on the intrinsic properties and population distribution of late-type stars. Moreover, a comparison of the Galactic parameters obtained with Gaia and VLBI can be done using radio observations on different targets: young massive stars (BeSSeL) and evolved stars (BAaDE).

Disks around brown dwarfs (BDs) are excellent laboratories to study the first steps of planet formation in cold and low-mass disk conditions. The radial-drift velocities of dust particles in BD disks are higher than in disks around more massive stars. Therefore, BD disks are expected to be more depleted in millimeter-sized grains compared to disks around T Tauri or Herbig Ae/Be stars. However, recent millimeter observations of BD disks revealed low millimeter spectral indices, indicating the presence of large grains in these disks and challenging models of dust evolution. We present 3\,mm photometric observations carried out with the IRAM/Plateau de Bure Interferometer (PdBI) of three BD disks in the Taurus star forming region, which have been observed with ALMA at 0.89\,mm. The disks were not resolved and only one was detected with enough confidence ($\sim3.5\sigma$) with PdBI. Based on these observations, we obtain the values and lower limits of the spectral index and find low values ($\alpha_{\rm{mm}}\lesssim 3.0$). We compare these observations in the context of particle trapping by an embedded planet, a promising mechanism to explain the observational signatures in more massive and warmer disks. We find, however, that this model cannot reproduce the current millimeter observations for BD disks, and multiple-strong pressure bumps globally distributed in the disk remain as a favorable scenario to explain observations. Alternative possibilities are that the gas masses in BD disk are very low ($\sim2\times10^{-3}\,M_{\rm{Jup}}$) such that the millimeter grains are decoupled and do not drift, or fast growth of fluffy aggregates.

Using the VLBA, the BeSSeL survey has provided distances and proper motions of young massive stars, allowing an accurate measure of the Galactic spiral structure. By the same technique, we are planning to map the inner Galaxy using positions and velocities of evolved stars (provided by the BAaDE survey). These radio astrometric measurements (BeSSeL and BAaDE) will be complementary to Gaia results and the overlap will provide important clues on the intrinsic properties and population distribution of the stars in the bulge.

The Bulge Asymmetries and Dynamical Evolution (BAaDE) survey will provide positions and line-of-sight velocities of ~20,000 evolved, maser bearing stars in the Galactic plane. Although this Galactic region is affected by optical extinction, BAaDE targets may have Gaia cross-matches, eventually providing additional stellar information. In an initial attempt to cross-match BAaDE targets with Gaia, we have found more than 5,000 candidates. Of these, we may expect half to show SiO emission, which will allow us to obtain velocity information. The cross-match is being refined to avoid false positives using different criteria based on distance analysis, flux variability, and color assessment in the mid- and near-IR. Once the cross-matches can be confirmed, we will have a unique sample to characterize the stellar population of evolved stars in the Galactic bulge, which can be considered fossils of the Milky Way formation.

Using trigonometric parallaxes and proper motions of masers associated with massive young stars, the Bar and Spiral Structure Legacy (BeSSeL) survey has reported the most accurate values of the Galactic parameters so far. The determination of these parameters with high accuracy has a widespread impact on Galactic and extragalactic measurements. This research is aimed at establishing the confidence with which such parameters can be determined. This is relevant for the data published in the context of the BeSSeL survey collaboration, but also for future observations, in particular from the Southern Hemisphere. In addition, some astrophysical properties of the masers can be constrained, notably the luminosity function. We have simulated the population of maser-bearing young stars associated with Galactic spiral structure, generating several samples and comparing them with the observed samples used in the BeSSeL survey. Consequently, we checked the determination of Galactic parameters for observational biases introduced by the sample selection. Galactic parameters obtained by the BeSSeL survey do not seem to be biased by the sample selection used. In fact, the published error estimates appear to be conservative for most of the parameters. We show that future BeSSeL data and future observations with Southern arrays will improve the Galactic parameters estimates and smoothly reduce their mutual correlation. Moreover, by modeling future parallax data with larger distance and, thus, greater relative uncertainties for a larger numbers of sources, we found that parallax-distance biasing is an important issue. Hence, using fractional parallax uncertainty in the weighting of the motion data is imperative. Finally, the luminosity function for 6.7 GHz methanol masers was determined, allowing us to estimate the number of Galactic methanol masers.

Transition disks are protoplanetary disks with inner depleted dust cavities and excellent candidates to investigate the dust evolution under the existence of a pressure bump. A pressure bump at the outer edge of the cavity allows dust grains from the outer regions to stop their rapid inward migration towards the star and efficiently grow to millimetre sizes. Dynamical interactions with planet(s) have been one of the most exciting theories to explain the clearing of the inner disk. We look for evidence of the presence of millimetre dust particles in transition disks by measuring their spectral index with new and available photometric data. We investigate the influence of the size of the dust depleted cavity on the disk integrated millimetre spectral index. We present the 3mm photometric observations carried out with PdBI of four transition disks: LkHa330, UXTauA, LRLL31, and LRLL67. We use available values of their fluxes at 345GHz to calculate their spectral index, as well as the spectral index for a sample of twenty transition disks. We compare the observations with two kind of models. In the first set of models, we consider coagulation and fragmentation of dust in a disk in which a cavity is formed by a massive planet located at different positions. The second set of models assumes disks with truncated inner parts at different radius and with power-law dust size distributions, where the maximum size of grains is calculated considering turbulence as the source of destructive collisions. We show that the integrated spectral index is higher for transition disks than for regular protoplanetary disks. For transition disks, the probability that the measured spectral index is positively correlated with the cavity radius is 95%. High angular resolution imaging of transition disks is needed to distinguish between the dust trapping scenario and the truncated disk case.