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.

We present (sub)millimeter imaging at 0.5'' resolution of the massive star-forming region G358.93-0.03 acquired in multiple epochs at 2 and 3 months following the recent flaring of its 6.7 GHz methanol maser emission. Using SMA and ALMA, we have discovered 14 new Class II methanol maser lines ranging in frequency from 199 GHz to 361 GHz, which originate mostly from vt=1 torsionally-excited transitions and include one vt=2 transition. The latter detection provides the first observational evidence that Class II maser pumping involves levels in the vt=2 state. The masers are associated with the brightest continuum source (MM1), which hosts a line-rich hot core. The masers present a consistent curvilinear spatial velocity pattern that wraps around MM1, suggestive of a coherent physical structure 1200 au in extent. In contrast, the thermal lines exhibit a linear pattern that crosses MM1 but at progressive position angles that appear to be a function of either increasing temperature or decreasing optical depth. The maser spectral profiles evolved significantly over one month, and the intensities dropped by factors of 3.0 to 7.2, with the vt=2 line showing the largest change. A small area of maser emission from only the highest excitation lines closest to MM1 has disappeared. There are seven additional dust continuum sources in the protocluster, including another hot core (MM3). We do not find evidence for a significant change in (sub)millimeter continuum emission from any of the sources during the one month interval, and the total protocluster emission remains comparable to prior single dish measurements.

Terzan 5 is the only Galactic globular cluster that has plausibly been detected at very-high energies by the High Energy Stereoscopic System. It has an unexpectedly asymmetric very-high-energy morphology that is offset from the cluster center, in addition to a large-scale, offset radio structure and compact diffuse X-ray emission associated with this cluster. We present new data from the Fermi Large Area Telescope on this source. We model the updated broadband spectral energy distribution, attributing this to cumulative pulsed emission from a population of embedded millisecond pulsars as well as unpulsed emission from the interaction of their leptonic winds with the ambient magnetic and soft-photon fields. In particular, our model invokes unpulsed synchrotron and inverse Compton components to model the radio and TeV data, cumulative pulsed curvature radiation to fit the Fermi data, and explains the hard Chandra X-ray spectrum via a "new" cumulative synchrotron component from electron-positron pairs within the pulsar magnetospheres that has not been implemented before. We find reasonable spectral fits for plausible model parameters. We also derive constraints on the millisecond pulsar luminosity function using the diffuse X-ray data and the Chandra sensitivity. Future higher-quality spectral and spatial data will help discriminate between competing scenarios (such as dark matter annihilation, white dwarf winds or hadronic interactions) proposed for the broadband emission as well as constrain degenerate model parameters.

Spectral lines are among the most powerful signatures for dark matter (DM) annihilation searches in very-high-energy $\gamma$-rays. The central region of the Milky Way halo is one of the most promising targets given its large amount of DM and proximity to Earth. We report on a search for a monoenergetic spectral line from self-annihilations of DM particles in the energy range from 300 GeV to 70 TeV using a two-dimensional maximum likelihood method taking advantage of both the spectral and spatial features of signal versus background. The analysis makes use of Galactic Center (GC) observations accumulated over ten years (2004 - 2014) with the H.E.S.S. array of ground-based Cherenkov telescopes. No significant $\gamma$-ray excess above the background is found. We derive upper limits on the annihilation cross section $\langle\sigma v\rangle$ for monoenergetic DM lines at the level of $\sim4\times10^{-28}$ cm$^{3}$s$^{-1}$ at 1 TeV, assuming an Einasto DM profile for the Milky Way halo. For a DM mass of 1 TeV, they improve over the previous ones by a factor of six. The present constraints are the strongest obtained so far for DM particles in the mass range 300 GeV - 70 TeV. Ground-based $\gamma$-ray observations have reached sufficient sensitivity to explore relevant velocity-averaged cross sections for DM annihilation into two $\gamma$-ray photons at the level expected from the thermal relic density for TeV DM particles.

A search for dark matter line-like signals was performed in the vicinity of the Galactic Centre by the H.E.S.S. experiment on observational data taken in 2014. An unbinned likelihood analysis was developed to improve the sensitivity to line-like signals. The upgraded analysis along with newer data extend the energy coverage of the previous measurement down to 100 GeV. The 18 h of data collected with the H.E.S.S. array allow one to rule out at 95% CL the presence of a 130 GeV line (at $l = -1.5^{\circ}, b = 0^{\circ}$ and for a dark matter profile centered at this location) previously reported in Fermi-LAT data. This new analysis overlaps significantly in energy with previous Fermi-LAT and H.E.S.S. results. No significant excess associated with dark matter annihilations was found in the energy range 100 GeV to 2 TeV and upper limits on the gamma-ray flux and the velocity weighted annihilation cross-section are derived adopting an Einasto dark matter halo profile. Expected limits for present and future large statistics H.E.S.S. observations are also given.

The inner region of the Milky Way halo harbors a large amount of dark matter (DM). Given its proximity, it is one of the most promising targets to look for DM. We report on a search for the annihilations of DM particles using $\gamma$-ray observations towards the inner 300 parsecs of the Milky Way, with the H.E.S.S. array of ground-based Cherenkov telescopes. The analysis is based on a 2D maximum likelihood method using Galactic center (GC) data accumulated by H.E.S.S. over the last 10 years (2004-2014), and does not show any significant $\gamma$-ray signal above background. Assuming Einasto and Navarro-Frenk-White DM density profiles at the GC, we derive upper limits on the annihilation cross section $\langle \sigma v\rangle$. These constraints are the strongest obtained so far in the TeV DM mass range and improve upon previous limits by a factor 5. For the Einasto profile, the constraints reach $\langle \sigma v\rangle$ values of $\rm 6\times10^{-26} cm^3s^{-1}$ in the $W^+W^-$ channel for a DM particle mass of 1.5 TeV, and $\rm 2\times10^{-26} cm^3s^{-1}$ in the $\tau^+\tau^-$ channel for 1 TeV mass. For the first time, ground-based $\gamma$-ray observations have reached sufficient sensitivity to probe $\langle \sigma v\rangle$ values expected from the thermal relic density for TeV DM particles.

Microquasars are potential $\gamma$-ray emitters. Indications of transient episodes of $\gamma$-ray emission were recently reported in at least two systems: Cyg X-1 and Cyg X-3. The identification of additional $\gamma$-ray-emitting microquasars is required to better understand how $\gamma$-ray emission can be produced in these systems. Theoretical models have predicted very high-energy (VHE) $\gamma$-ray emission from microquasars during periods of transient outburst. Observations reported herein were undertaken with the objective of observing a broadband flaring event in the $\gamma$-ray and X-ray bands. Contemporaneous observations of three microquasars, GRS 1915+105, Circinus X-1, and V4641 Sgr, were obtained using the High Energy Spectroscopic System (H.E.S.S.) telescope array and the Rossi X-ray Timing Explorer (RXTE) satellite. X-ray analyses for each microquasar were performed and VHE $\gamma$-ray upper limits from contemporaneous H.E.S.S. observations were derived. No significant $\gamma$-ray signal has been detected in any of the three systems. The integral $\gamma$-ray photon flux at the observational epochs is constrained to be $I(> 560 {\rm\ GeV}) < 7.3\times10^{-13}$ cm$^{-2}$ s$^{-1}$, $I(> 560 {\rm\ GeV}) < 1.2\times10^{-12}$ cm$^{-2}$ s$^{-1}$, and $I(> 240 {\rm\ GeV}) < 4.5\times10^{-12}$ cm$^{-2}$ s$^{-1}$ for GRS 1915+105, Circinus X-1, and V4641 Sgr, respectively. The $\gamma$-ray upper limits obtained using H.E.S.S. are examined in the context of previous Cherenkov telescope observations of microquasars. The effect of intrinsic absorption is modelled for each target and found to have negligible impact on the flux of escaping $\gamma$-rays. When combined with the X-ray behaviour observed using RXTE, the derived results indicate that if detectable VHE $\gamma$-ray emission from microquasars is commonplace, then it is likely to be highly transient.

Using the High Energy Spectroscopic System (H.E.S.S.) telescopes we have discovered a steady and extended very high-energy (VHE) gamma-ray source towards the luminous blue variable candidate LBV1806-20, massive stellar cluster Cl*1806-20, and magnetar SGR1806-20. The new VHE source, HESSJ1808-204, was detected at a statistical significance of >6sigma (post-trial) with a photon flux normalisation $(2.9 \pm 0.4_{\rm stat} \pm 0.5_{\rm sys})\times 10^{-13}$\,ph\,cm$^{-2}$\,s$^{-1}$\,TeV$^{-1}$ at 1TeV and a power-law photon index of $2.3\pm0.2_{\rm stat}\pm 0.3_{\rm sys}$. The luminosity of this source (0.2 to 10 TeV; scaled to distance d=8.7kpc) is $L_{\rm VHE}\sim1.6 \times 10^{34}(d/{8.7kpc})^2$\,erg\,s$^{-1}$. The VHE gamma-ray emission is extended and is well fit by a single Gaussian with statistical standard deviation of $0.095^\circ \pm 0.015^\circ$. This extension is similar to that of the synchrotron radio nebula G10.0-0.3, which is thought to be powered by LBV1806-20. The VHE gamma-ray luminosity could be provided by the stellar wind luminosity of LBV1806-20 by itself and/or the massive star members of Cl*1806-20. Alternatively, magnetic dissipation (e.g. via reconnection) from SGR1806-20 can potentially account for the VHE luminosity. The origin and hadronic and/or leptonic nature of the accelerated particles responsible for HESSJ1808-204 is not yet clear. If associated with SGR1806-20, the potentially young age of the magnetar (650yr) can be used to infer the transport limits of these particles to match the VHE source size. This discovery provides new interest in the potential for high-energy particle acceleration from magnetars, massive stars, and/or stellar clusters.

Galactic cosmic rays reach energies of at least a few Peta-electronvolts (1 PeV =$10^\mathbf{15}$ electron volts). This implies our Galaxy contains PeV accelerators (PeVatrons), but all proposed models of Galactic cosmic-ray accelerators encounter non-trivial difficulties at exactly these energies. Tens of Galactic accelerators capable of accelerating particle to tens of TeV (1 TeV =$10^\mathbf{12}$ electron volts) energies were inferred from recent gamma-ray observations. None of the currently known accelerators, however, not even the handful of shell-type supernova remnants commonly believed to supply most Galactic cosmic rays, have shown the characteristic tracers of PeV particles: power-law spectra of gamma rays extending without a cutoff or a spectral break to tens of TeV. Here we report deep gamma-ray observations with arcminute angular resolution of the Galactic Centre regions, which show the expected tracer of the presence of PeV particles within the central 10~parsec of the Galaxy. We argue that the supermassive black hole Sagittarius A* is linked to this PeVatron. Sagittarius A* went through active phases in the past, as demonstrated by X-ray outbursts and an outflow from the Galactic Centre. Although its current rate of particle acceleration is not sufficient to provide a substantial contribution to Galactic cosmic rays, Sagittarius A* could have plausibly been more active over the last $\gtrsim 10^{6-7}$ years, and therefore should be considered as a viable alternative to supernova remnants as a source of PeV Galactic cosmic rays.

Aims: We aim for an understanding of the morphological and spectral properties of the supernova remnant RCW~86 and for insights into the production mechanism leading to the RCW~86 very high-energy gamma-ray emission. Methods: We analyzed High Energy Spectroscopic System data that had increased sensitivity compared to the observations presented in the RCW~86 H.E.S.S. discovery publication. Studies of the morphological correlation between the 0.5-1~keV X-ray band, the 2-5~keV X-ray band, radio, and gamma-ray emissions have been performed as well as broadband modeling of the spectral energy distribution with two different emission models. Results:We present the first conclusive evidence that the TeV gamma-ray emission region is shell-like based on our morphological studies. The comparison with 2-5~keV X-ray data reveals a correlation with the 0.4-50~TeV gamma-ray emission.The spectrum of RCW~86 is best described by a power law with an exponential cutoff at $E_{cut}=(3.5\pm 1.2_{stat})$ TeV and a spectral index of $\Gamma$~$1.6\pm 0.2$. A static leptonic one-zone model adequately describes the measured spectral energy distribution of RCW~86, with the resultant total kinetic energy of the electrons above 1 GeV being equivalent to $\sim$0.1\% of the initial kinetic energy of a Type I a supernova explosion. When using a hadronic model, a magnetic field of $B$~100$\mu$G is needed to represent the measured data. Although this is comparable to formerly published estimates, a standard E$^{-2}$ spectrum for the proton distribution cannot describe the gamma-ray data. Instead, a spectral index of $\Gamma_p$~1.7 would be required, which implies that ~$7\times 10^{49}/n_{cm^{-3}}$erg has been transferred into high-energy protons with the effective density $n_{cm^{-3}}=n/ 1$ cm^-3. This is about 10\% of the kinetic energy of a typical Type Ia supernova under the assumption of a density of 1~cm^-3.

Re-observations with the H.E.S.S. telescope array of the very-high-energy (VHE) source HESS J1018-589 A coincident with the Fermi-LAT $\gamma$-ray binary 1FGL J1018.6-5856 have resulted in a source detection significance of more than 9$\sigma$, and the detection of variability ($\chi^2$/$\nu$ of 238.3/155) in the emitted $\gamma$-ray flux. This variability confirms the association of HESS J1018-589 A with the high-energy $\gamma$-ray binary detected by Fermi-LAT, and also confirms the point-like source as a new very-high-energy binary system. The spectrum of HESS J1018-589 A is best fit with a power-law function with photon index $\Gamma = 2.20 \pm 0.14_{\rm stat} \pm 0.2_{\rm sys}$. Emission is detected up to ~20 TeV. The mean differential flux level is $(2.9 \pm 0.4)\times10^{-13}$ TeV$^{-1}$ cm$^{-2}$ s$^{-1}$ at 1 TeV, equivalent to ~1% of the flux from the Crab Nebula at the same energy. Variability is clearly detected in the night-by-night lightcurve. When folded on the orbital period of 16.58 days, the rebinned lightcurve peaks in phase with the observed X-ray and high-energy phaseograms. The fit of the H.E.S.S. phaseogram to a constant flux provides evidence of periodicity at the level of 3$\sigma$. The shape of the VHE phaseogram and measured spectrum suggest a low inclination, low eccentricity system with a modest impact from VHE $\gamma$-ray absorption due to pair production ($\tau$ < 1 at 300 GeV).

QU Car is listed in cataclysmic variable star catalogues as a nova-like variable. This little-studied, yet bright interacting binary is re-appraised here in the light of new high-quality ultraviolet (UV) interstellar line data obtained with STIS on board the Hubble Space Telescope. The detection of a component of interstellar absorption at a mean LSR velocity of $-$14 km s$^{-1}$ indicates that the distance to QU Car may be $\sim$2 kpc or more -- a considerable increase on the previous lower-limiting distance of 500 pc. If so, the bolometric luminosity of QU Car could exceed $10^{37}$ ergs s$^{-1}$. This would place this binary in the luminosity domain occupied by known compact-binary supersoft X-ray sources. Even at a 500 pc, QU Car appears to be the most luminous nova-like variable known. New intermediate dispersion optical spectroscopy of QU Car spanning 3800--7000 \AA is presented. These data yield the discovery that C{\sc iv} $\lambda\lambda$5801,12 is present as an unusually prominent emission line in an otherwise low-contrast line spectrum. Using measurements of this and other lines in a recombination line analysis, it is shown that the C/He abundance as proxied by the n(C$^{4+}$)/n(He$^{2+}$) ratio may be as high as 0.06 (an order of magnitude higher than the solar ratio). Furthermore, the C/O abundance ratio is estimated to be greater than 1. These findings suggest that the companion in QU Car is a carbon star. If so, it would be the first example of a carbon star in such a binary. An early-type R star best matches the required abundance pattern and could escape detection at optical wavelengths provided the distance to QU Car is $\sim$2 kpc or more.

Optical, near-infrared, and radio observations of the BL Lac object PKS2155-304 were obtained simultaneously with a continuous UV/EUV/X-ray monitoring campaign in 1994 May. Further optical observations were gathered throughout most of 1994. The radio, millimeter, and near-infrared data show no strong correlations with the higher energies. The optical light curves exhibit flickering of 0.2-0.3 mag on timescales of 1-2 days, superimposed on longer timescale variations. Rapid variations of ~0.01 mag/min, which, if real, are the fastest seen to date for any BL Lac object. Small (0.2-0.3 mag) increases in the V and R bands occur simultaneously with a flare seen at higher energies. All optical wavebands (UBVRI) track each other well over the period of observation with no detectable delay. For most of the period the average colors remain relatively constant, although there is a tendency for the colors (in particular B-V) to vary more when the source fades. In polarized light, PKS 2155-304 showed strong color dependence and the highest optical polarization (U = 14.3%) ever observed for this source. The polarization variations trace the flares seen in the ultraviolet flux.