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.

Fanaroff-Riley II low-excitation radio galaxies (FRII-LERGs) are characterized by weak nuclear excitation on pc-scales and by properties typical of powerful FRIIs (defined as high-excitation, hereafter HERGs/BLRGs) on kp-scales. Since a link between the accretion properties and the power of the produced jets is expected both from theory and observations, their nature is still debated. In this work we investigate the X-ray properties of a complete sample of 19 FRII-LERGs belonging to the 3CR catalog, exploiting Chandra and XMM-Newton archival data. We also analyze 32 FRII-HERGs/BLRGs with Chandra data as a control sample. We compared FRII-LERG and FRII-HERG/BLRG X-ray properties and optical data available in literature to obtain a wide outlook of their behavior. The low accretion rate estimates for FRII-LERGs, from both X-ray and optical bands, allow us to firmly reject the hypothesis for that they are the highly obscured counterpart of powerful FRII-HERGs/BLRGs. Therefore, at least two hypothesis can be invoked to explain the FRII-LERGs nature: (i) they are evolving from classical FRIIs because of the depletion of accreting cold gas in the nuclear region, while the extended radio emission is the heritage of a past efficiently accreting activity; (ii) they are an intrinsically distinct class of objects with respect to classical FRIs/FRIIs. Surprisingly, in this direction a correlation between accretion rates and environmental richness is found in our sample. The richer the environment, the more inefficient is the accretion. In this framework, the FRII-LERGs are intermediate between FRIs and FRII-HERGs/BLRGs both in terms of accretion rate and environment.

The $\gamma$-ray production mechanism and its localization in blazars are still a matter of debate. The main goal of this paper is to constrain the location of the high-energy emission in the blazar TXS 2013+370 and to study the physical and geometrical properties of the inner jet region on sub-pc scales. VLBI observations at 86 GHz and space-VLBI at 22 GHz allowed us to image the jet base with an angular resolution of $\sim$0.4 pc. By employing CLEAN imaging and Gaussian model-fitting, we performed a thorough kinematic analysis, which provided estimates of the jet speed, orientation, and component ejection times. Additionally, we studied the jet expansion profile and used the information on the jet geometry to estimate the location of the jet apex. VLBI data were combined with single-dish measurements to search for correlated activity between the radio and $\gamma$-ray emission. The high-resolution VLBI imaging revealed the existence of a spatially bent jet, described by moving and stationary features. New jet features are observed to emerge from the core, accompanied by flaring activity in radio bands and $\gamma$ rays. The analysis of the transverse jet width profile constrains the location of the mm core to lie $\leq$ 2 pc downstream of the jet apex, and also reveals the existence of a transition from parabolic to conical jet expansion at a distance of $\sim$54 pc from the core, corresponding to $\sim$1.5$\times$10$^{\rm 6}$ Schwarzschild radii. The cross-correlation analysis reveals a strong correlation between the radio and $\gamma$-ray data, with the 1 mm emission lagging $\sim$49 days behind the $\gamma$ rays. Based on this, we infer that the high energy emission is produced at a distance of $\sim$1 pc from the VLBI core, suggesting that the seed photon fields for the external Compton mechanism originate either in the dusty torus or in the broad-line region.

IceCube has reported a very-high-energy neutrino (IceCube-170922A) in a region containing the blazar TXS 0506+056. Correlated {\gamma}-ray activity has led to the first high-probability association of a high-energy neutrino with an extragalactic source. This blazar has been found to be in a radio outburst during the neutrino event. We have performed target-of-opportunity VLBI imaging observations at 43 GHz frequency with the VLBA two and eight months, respectively, after the neutrino event. We produced two images of TXS 0506+056 with angular resolutions of (0.2x1.1) mas and (0.2x0.5) mas, respectively. The source shows a compact, high brightness temperature core (albeit not approaching the equipartition limit) and a bright and originally very collimated inner jet. Beyond about 0.5 mas from the mm-VLBI core, the jet loses this tight collimation and expands rapidly. During the months after the neutrino event associated with this source, the overall flux density is rising. This flux density increase happens solely within the core. The core expands in size with apparent superluminal velocity during these six months so that the brightness temperature drops by a factor of three in spite of the strong flux density increase. The radio jet of TXS 0506+056 shows strong signs of deceleration and/or a spine-sheath structure within the inner 1 mas (corresponding to about 70 pc to 140 pc in deprojected distance) from the mm-VLBI core. This structure is consistent with theoretical models that attribute the neutrino and {\gamma}-ray production to interactions of electrons and protons in the highly-relativistic jet spine with external photons originating from a slower-moving jet region. Proton loading due to jet-star interactions in the inner host galaxy is suggested as the possible cause of deceleration

In March 2018, Mukherjee (2018) reported the detection by VERITAS of very-high-energy emission (VHE; > 100 GeV) from 3C 264. This is the sixth, and second most distant, radio galaxy ever detected in the TeV regime. In this article we present a radio and X-ray analysis of the jet in 3C 264. We determine the main physical parameters of the parsec-scale flow and explore the implications of the inferred kinematic structure for radiative models of this gamma-ray emitting jet. The radio data set is comprised of VLBI observations at 15 GHz from the MOJAVE program, and cover a time period of ~2 years. Through a segmented wavelet decomposition method we estimate the apparent displacement of individual plasma features; we then perform a pixel-based analysis of the stacked image to determine the jet shape. The X-ray data set includes all available observations from the Chandra, XMM, and Swift satellites, and is used, together with archival data in the other bands, to build the SED. Proper motion is mostly detected along the edges of the flow, which appears strongly limb-brightened. The apparent speeds increase as a function of distance from the core up to a maximum of ~$11.5$ c. This constrains the jet viewing angle to assume relatively small values ($\theta\lesssim10^{\circ}$). In the acceleration region, extending up to a de-projected distance of ~$4.8\times10^4$ Schwarzschild radii (~$11$ pc), the jet is collimating ($r\propto z^{0.40\pm 0.04}$), as predicted for a magnetically-driven plasma flow. By assuming that the core region is indeed magnetically dominated ($U_B/U_e>1$), the SED and the jet power can be well reproduced in the framework of leptonic models, provided that the high-energy component is associated to a second emitting region. The possibility that this region is located at the end of the acceleration zone, either in the jet layer or in the spine, is explored in the modeling.

Investigating the magnetic field structure in the innermost regions of relativistic jets is fundamental to shed light on the crucial physical processes giving rise to the jet formation, as well as to its extraordinary radiation output up to gamma-ray energies. We study the magnetic field structure of the quasar CTA 102 with 3 and 7 mm-VLBI polarimetric observations, reaching an unprecedented resolution (~50 microarcsec). We also investigate the variability and physical process occurring in the source during the observing period which coincides with a very active state of the source till high-energies. The Faraday rotation analysis between 3 and 7mm shows a gradient in rotation measure with a maximum value of ~6X10^4 rad/m^2 and intrinsic electric vector position angles (EVPAs) oriented around the centroid of the core, suggesting the presence of large-scale helical magnetic fields. Such a magnetic field structure is also visible in 7 mm images when a new superluminal component is crossing the core region. The 7mm EVPAs orientation is different when the component is exiting the core or crossing a stationary feature at ~0.1 mas. The interaction between the superluminal component and a recollimation shock at ~0.1 mas could have triggered the multi-wavelengths flares. The variability Doppler factor associated with such interaction is large enough to explain the high energy emission, as we infer from the analysis of gamma-ray and X-ray data, and it is in agreement with the Doppler factor obtained to explain the extraordinary optical flare by Raiteri et al.(2017).

Over the past few decades, our knowledge of jets produced by active galactic nuclei (AGN) has greatly progressed thanks to the development of very-long-baseline interferometry (VLBI). Nevertheless, the crucial mechanisms involved in the formation of the plasma flow, as well as those driving its exceptional radiative output up to TeV energies, remain to be clarified. Most likely, these physical processes take place at short separations from the supermassive black hole, on scales which are inaccessible to VLBI observations at centimeter wavelengths. Due to their high synchrotron opacity, the dense and highly magnetized regions in the vicinity of the central engine can only be penetrated when observing at shorter wavelengths, in the millimeter and sub-millimeter regimes. While this was recognized already in the early days of VLBI, it was not until the very recent years that sensitive VLBI imaging at high frequencies has become possible. Ongoing technical development and wide band observing now provide adequate imaging fidelity to carry out more detailed analyses. In this article we overview some open questions concerning the physics of AGN jets, and we discuss the impact of mm-VLBI studies. Among the rich set of results produced so far in this frequency regime, we particularly focus on studies performed at 43 GHz (7 mm) and at 86 GHz (3 mm). Some of the first findings at 230 GHz (1 mm) obtained with the Event Horizon Telescope are also presented.

We present for the first time Very-Long-Baseline Interferometry images of the radio galaxy Cygnus A at the frequency of $86$ $\rm GHz$. Thanks to the high spatial resolution of only ${\sim}200$ Schwarzschild radii ($R_{\bf S}$), such observations provide an extremely detailed view of the nuclear regions in this archetypal object and allow us to derive important constraints for theoretical models describing the launching of relativistic jets. A pixel-based analysis of the jet outflow, which still appears two-sided on the scales probed, was performed. By fitting Gaussian functions to the transverse intensity profiles, we could determine the jet width in the nuclear region. The base of the jets appears wide. The minimum measured transverse width of ${\sim} (227\pm98)$ $R_{\bf S}$ is significantly larger than the radius of the Innermost Stable Circular Orbit, suggesting that the outer accretion disk is contributing to the jet launching. The existence of a faster and Doppler de-boosted inner section, powered either from the rotation of the inner regions of the accretion disk or by the spinning black hole, is suggested by the kinematic properties and by the observed limb brightening of the flow.

High-resolution Very-Long-Baseline Interferometry observations of relativistic jets are essential to constrain fundamental parameters of jet formation models. At a distance of 249 Mpc, Cygnus A is a unique target for such studies, being the only Fanaroff-Riley Class II radio galaxy for which a detailed sub-parsec scale imaging of the base of both jet and counter-jet can be obtained. Observing at millimeter wavelengths unveils those regions which appear self-absorbed at longer wavelengths and enables an extremely sharp view towards the nucleus to be obtained. We performed 7 mm Global VLBI observations, achieving ultra-high resolution imaging on scales down to 90 $\mu$as. This resolution corresponds to a linear scale of only ${\sim}$400 Schwarzschild radii (for $M_{\mathrm{BH}}=2.5 \times 10^9 M_{\odot}$). We studied the kinematic properties of the main emission features of the two-sided flow and probed its transverse structure through a pixel-based analysis. We suggest that a fast and a slow layer, with different acceleration gradients, exist in the flow. The extension of the acceleration region is large (${\sim} 10^4 R_{\mathrm{S}}$), indicating that the jet is magnetically-driven. The limb brightening of both jet and counter-jet and their large opening angles ($\phi_\mathrm{J}{\sim} 10^{\circ}$) strongly favor a spine-sheath structure. In the acceleration zone, the flow has a parabolic shape ($r\propto z^{0.55\pm 0.07}$). The acceleration gradients and the collimation profile are consistent with the expectations for a jet in "equilibrium'' (Lyubarsky 2009), achieved in the presence of a mild gradient of the external pressure ($p\propto z^{-k}, k\leq2$).}

At a distance of 249 Mpc ($z$=0.056), Cygnus A is the only powerful FR II radio galaxy for which a detailed sub-parsec scale imaging of the base of both jet and counter-jet can be obtained. Observing with VLBI at millimeter wavelengths is fundamental for this object, as it uncovers those regions which appear self-absorbed or free-free absorbed by a circumnuclear torus at longer wavelengths. We performed 7 mm Global VLBI observations, achieving ultra-high resolution imaging on scales down to 90 $\mu$as. This resolution corresponds to a linear scale of only $\sim$400 Schwarzschild radii. We studied the transverse structure of the jets through a pixel-based analysis, and kinematic properties of the main emission features by modeling the interferometric visibilities with two-dimensional Gaussian components. Both jets appear limb-brightened, and their opening angles are relatively large ($\phi_\mathrm {j}\sim 10^{\circ}$). The flow is observed to accelerate within the inner-jet up to scales of $\sim$1 pc, while lower speeds and uniform motions are measured further downstream. A single component seen in the counter-jet appears to be stationary. These observational properties are explained assuming the existence of transverse gradient of the bulk Lorentz factor across the jet, consisting of a fast central spine surrounded by a slower boundary layer.

Context. The advent of Fermi is changing our understanding on the radio and gamma-ray emission in Active Galactic Nuclei. Contrary to pre-Fermi ideas, BL Lac objects are found to be the most abundant emitters in the gamma-ray band. However, since they are relatively weak radio sources, most of their parsec-scale structure and their multi-frequency properties are poorly understood and/or have not been investigated in a systematically fashion. Aims. Our main goal is to analyze the radio and gamma-ray emission properties of a sample of 42 BL Lacs selected, for the first time in the literature, with no constraint on their radio and gamma-ray flux densities/emission. Methods. Thanks to new Very Long Baseline Array observations at 8 and 15 GHz for the whole sample, we present here fundamental parameters such as radio flux densities, spectral index information, and parsec-scale structure. Moreover, we search for gamma-ray counterparts using data reported in the Second Catalog of Fermi Gamma-ray sources. Results. Parsec-scale radio emission is observed in the majority of the sources at both frequencies. Gamma-ray counterparts are found for 14/42 sources. Conclusions. The comparison between our results in radio and gamma-ray bands points out the presence of a large number of faint BL Lacs showing "non classical" properties such as low source compactness, core dominance, no gamma-ray emission and steep radio spectral indexes. A deeper multiwavelength analysis will be needed.

The advent of Fermi is changing our understanding on the radio and gamma-ray emission in Active Galactic Nuclei. In fact, contrary to previous campaigns, Fermi mission reveals that BL Lac objects are the most abundant emitters in gamma-ray band. However, since they are relatively weak sources, most of their parsec scale structure as their multifrequency properties are poorly understood and/or not systematically investigated. Our main goal is to analyse, using a multiwavelength approach, the nuclear properties of an homogeneous sample of 42 faint BL Lacs, selected, for the first time in literature, with no constraint on their radio and gamma-ray flux densities/emission. We began asking and obtaining new VLBA observations at 8 and 15 GHz for the whole sample. We derived fundamental parameters as radio flux densities, spectral index information, and parsec scale structure. Moreover, we investigated their gamma-ray emission properties using the 2LAT Fermi results. Here, we report our preliminary results on the radio and gamma-ray properties of this sample of faint BL Lacs. In the next future, we will complete the multiwavelength analysis.