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 aim to investigate the overall properties of the ICRF3 with the help of the Gaia Data release 2 (Gaia DR2). This could serve as an external check of the quality of the ICRF3. The radio source positions of the ICRF3 catalog were compared with the Gaia DR2 positions of their optical counterparts at G < 18.7. Their properties were analyzed in terms of the dependency of the quoted error on the number of observations, on the declination, and the global difference, the latter revealed by means of expansions in the vector spherical harmonics. The ICRF3 S/X-band catalog shows a more smooth dependency on the number of observations than the ICRF1 and ICRF2, while the K and X/Ka-band yield a dependency discrepancy at the number of observations of about 50. The rotation of all ICRF catalogs show consistent results, except for the X-component of the X/Ka-band which arises from the positional error in the non-defining sources. No significant glides were found between the ICRF3 S/X-band component and Gaia DR2. However, the K- and X/Ka- band frames show a dipolar deformation in Y-component of +50{\mu}as and several quadrupolar terms of 50{\mu}as in an absolute sense. A significant glide along Z-axis exceeding 200 {\mu}as in the X/Ka-band was also reported. These systematics in the ICRF catalog are shown to be less dependent on the limiting magnitude of the Gaia sample when the number of common sources is sufficient (> 100). The ICRF3 S/X-band catalog shows improved accuracy and systematics at the level of noise floor. But the zonal errors in the X/Ka-band should be noted, especially in the context of comparisons of multi-frequency positions for individual sources.

Geodetic Very Long Baseline Interferometry (VLBI) data are capable of measuring the light deflection caused by the gravitational field of the Sun and large planets with high accuracy. The parameter $\gamma$ of the parametrized Post-Newtonian (PPN) formalism estimated using observations of reference radio sources near the Sun should be equal to unity in the general relativity. We have run several VLBI experiments tracking reference radio sources from 1 to 3 degrees from the Sun. The best formal accuracy of the parameter $\gamma$ achieved in the single-session mode is less than 0.01 percent, or better than the formal accuracy obtained with a global solution included all available observations at arbitrary elongation from the Sun. We are planning more experiments starting from 2020 using better observing conditions near the minimum of the Solar activity cycle.

We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.

Context. We highlight the capabilities of the geodetic VLBI technique to test General relativity in the classical astrometric style, i.e., measuring the deflection of light in the vicinity of the Sun. Aims. In previous studies, the parameter was estimated by global analyses of thousands of geodetic VLBI sessions. Here we estimate from a single session where the Sun has approached two strong reference radio sources 0229+131 and 0235+164 at an elongation angle of 1-3 degrees. Methods. The AUA020 VLBI session of 1 May 2017 was designed to obtain more than 1000 group delays from the two radio sources. The Solar corona effect was effectively calibrated with the dual-frequency observations even at small elongation from the Sun. Results. We obtained with a precision better than what is obtained through global analyses of thousands of standard geodetic sessions over decades. Current results demonstrate that the modern VLBI technology is capable of establishing new limits on observational test of General Relativity.

The second Gaia data release (DR2), contains very precise astrometric and photometric properties for more than one billion sources, astrophysical parameters for dozens of millions, radial velocities for millions, variability information for half a million of stellar sources and orbits for thousands of solar system objects. Before the Catalogue publication, these data have undergone dedicated validation processes. The goal of this paper is to describe the validation results in terms of completeness, accuracy and precision of the various Gaia DR2 data. The validation processes include a systematic analysis of the Catalogue content to detect anomalies, either individual errors or statistical properties, using statistical analysis, and comparisons to external data or to models. Although the astrometric, photometric and spectroscopic data are of unprecedented quality and quantity, it is shown that the data cannot be used without a dedicated attention to the limitations described here, in the Catalogue documentation and in accompanying papers. A particular emphasis is put on the caveats for the statistical use of the data in scientific exploitation.

The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration.

Before the publication of the Gaia Catalogue, the contents of the first data release have undergone multiple dedicated validation tests. These tests aim at analysing in-depth the Catalogue content to detect anomalies, individual problems in specific objects or in overall statistical properties, either to filter them before the public release, or to describe the different caveats of the release for an optimal exploitation of the data. Dedicated methods using either Gaia internal data, external catalogues or models have been developed for the validation processes. They are testing normal stars as well as various populations like open or globular clusters, double stars, variable stars, quasars. Properties of coverage, accuracy and precision of the data are provided by the numerous tests presented here and jointly analysed to assess the data release content. This independent validation confirms the quality of the published data, Gaia DR1 being the most precise all-sky astrometric and photometric catalogue to-date. However, several limitations in terms of completeness, astrometric and photometric quality are identified and described. Figures describing the relevant properties of the release are shown and the testing activities carried out validating the user interfaces are also described. A particular emphasis is made on the statistical use of the data in scientific exploitation.

Lorentz symmetry violations can be described by an effective field theory framework that contains both General Relativity and the Standard Model of particle physics called the Standard-Model extension (SME). Recently, post-fit analysis of Gravity Probe B and binary pulsars lead to an upper limit at the $10^{-4}$ level on the time-time coefficient $\bar s^{TT}$ of the pure-gravity sector of the minimal SME. In this work, we derive the observable of Very Long Baseline Interferometry (VLBI) in SME and then we implement it into a real data analysis code of geodetic VLBI observations. Analyzing all available observations recorded since 1979, we compare estimates of $\bar s^{TT}$ and errors obtained with various analysis schemes, including global estimations over several time spans and with various Sun elongation cut-off angles, and with analysis of radio source coordinate time series. We obtain a constraint on $\bar s^{TT}=(-5\pm 8)\times 10^{-5}$, directly fitted to the observations and improving by a factor 5 previous post-fit analysis estimates.

Lorentz symmetry violations can be parametrized by an effective field theory framework that contains both General Relativity and the Standard Model of particle physics, called the Standard-Model Extension or SME. We consider in this work only the pure gravitational sector of the minimal SME. We present new constraints on the SME coefficients obtained from lunar laser ranging, very long baseline interferometry, and planetary motions.

Some models of the expanding Universe predict that the astrometric proper motion of distant radio sources embedded in space-time are non-zero as the radial distance from observer to the source grows. Systematic proper motion effects would produce a predictable quadrupole pattern on the sky that could be detected using Very Long Baseline Interferometry (VLBI) technique. This quadrupole pattern can be interpreted either as an anisotropic Hubble expansion, or as a signature of the primordial gravitational waves in the early Universe. We present our analysis of a large set of geodetic VLBI data spanning 1979--2015 to estimate the dipole and quadrupole harmonics in the expansion of the vector field of the proper motions of quasars in the sky. The analysis is repeated for different redshift zones.

The goal of this presentation is to report the latest progress in creation of the next generation of VLBI-based International Celestial Reference Frame, ICRF3. Two main directions of ICRF3 development are improvement of the S/X-band frame and extension of the ICRF to higher frequencies. Another important task of this work is the preparation for comparison of ICRF3 with the new generation optical frame GCRF expected by the end of the decade as a result of the Gaia mission.

Modeling of the trajectories of VLBI components ejected by the nucleus of 1928+738 shows the VLBI jet contains three families of trajectories, i.e. VLBI components are ejected from three different origins. The fit of components C1, C6 and C8 indicates that the nucleus of 1928+738 contains two binary black hole systems. The first binary black hole system is associated with the stationary components Cg and CS and is characterized by a radius $R_{bin,1} \approx 0.220$ mas; both black holes ejected VLBI components quasi regularly between 1990 and 2010. The second binary black hole system is not associated with stationary components and is characterized by a radius $R_{bin,2} \approx 0.140$ mas; it ejected only three VLBI components between 1994 and 1999. The two black hole sytems are separated by $\approx 1.35$ mas. We briefly discuss the consequences of the existence of binary black holes systems in radio quasars to make the link between radio quasars and GAIA.

In this paper we outline several problems related to the realization of the international celestial and terrestrial reference frames ICRF and ITRF at the millimeter level of accuracy, with emphasis on ICRF issues. The main topics considered are: analysis of the current status of the ICRF, mutual impact of ICRF and ITRF, and some considerations for future ICRF realizations.

We propose new estimates of the secular aberration drift, mainly due to the rotation of the Solar System about the Galactic center, based on up-to-date VLBI observations and and improved method of outlier elimination. We fit degree-2 vector spherical harmonics to extragalactic radio source proper motion field derived from geodetic VLBI observations spanning 1979-2013. We pay particular attention to the outlier elimination procedure to remove outliers from (i) radio source coordinate time series and (ii) the proper motion sample. We obtain more accurate values of the Solar system acceleration compared to those in our previous paper. The acceleration vector is oriented towards the Galactic center within 7 deg. The component perpendicular to the Galactic plane is statistically insignificant. We show that an insufficient cleaning of the data set can lead to strong variations in the dipole amplitude and orientation, and statistically biased results.

This report exists to provide high-level guidance for the strategic and engineering development of Data Management and Preservation plans for 'Big Science' data. Although the report's nominal audience is therefore rather narrow, we intend the document to be of use to other planners and data architects who wish to implement good practice in this area. For the purposes of this report, we presume that the reader is broadly persuaded (by external fiat if nothing else) of the need to preserve research data appropriately, and that they have both sophisticated technical support and the budget to support developments. The goal of the document is not to provide mechanically applicable recipes, but to allow the user to develop and lead a high-level plan which is appropriate to their organisation. Throughout, the report is informed where appropriate by the OAIS reference model.

While analyzing decades of very long baseline interferometry (VLBI) data, we detected the secular aberration drift of the extragalatic radio source proper motions caused by the rotation of the Solar System barycenter around the Galactic center. Our results agree with the predicted estimate to be 4-6 micro arcseconds per year ({\mu}as/yr) towards {\alpha} = 266\circ and {\delta} = -29\circ. In addition, we tried to detect the quadrupole systematics of the velocity field. The analysis method consisted of three steps. First, we analyzed geodetic and astrometric VLBI data to produce radio source coordinate time series. Second, we fitted proper motions of 555 sources with long observational histories over the period 1990-2010 to their respective coordinate time series. Finally, we fitted vector spherical harmonic components of degrees 1 and 2 to the proper motion field. Within the error bars, the magnitude and the direction of the dipole component agree with predictions. The dipole vector has an amplitude of 6.4 \pm 1.5 {\mu}as/yr and is directed towards equatorial coordinates {\alpha} = 263\circ and {\delta} = -20\circ. The quadrupole component has not been detected. The primordial gravitational wave density, integrated over a range of frequencies less than 10-9 Hz, has a limit of 0.0042 h-2 where h is the normalized Hubble constant is H0/(100 km s-1).

Many of the compact extragalactic radio sources that are used as fiducial points to define the celestial reference frame are known to have proper motions detectable with long-term geodetic/astrometric Very Long Baseline Interferometry (VLBI) measurements. These changes can be as high as several hundred micro-arcseconds per year for certain objects. When imaged with VLBI at milli-arcsecond (mas) angular resolution, these sources (radio-loud active galactic nuclei) typically show structures dominated by a compact, often unresolved "core" and a one-sided "jet". The positional instability of compact radio sources is believed to be connected with changes in their brightness distribution structure. For the first time, we test this assumption in a statistical sense on a large sample, not only for individual objects. We investigate a sample of 62 radio sources for which reliable long-term time series of astrometric positions as well as detailed 8-GHz VLBI brightness distribution models are available. We compare the characteristic direction of their extended jet structure and the direction of their apparent proper motion. We present our data and analysis method, and conclude that there is indeed a correlation between the two characteristic directions. However, there are cases where the ~1-10-mas scale VLBI jet directions are significantly misaligned with respect to the apparent proper motion direction.

The International Celestial Reference Frame (ICRF, Ma et al. 1998) is currently the best realization of a quasi-inertial reference system. It is based on more than 10 years of cumulated geodetic and astrometric VLBI observations of compact extragalactic objects at centimetric wavelengths. In the perspective of the realization of an accurate optical counterpart of the ICRF using future space astrometry missions like GAIA or SIM, this paper investigates the consistency of celestial reference frames realized through the same subset of compact extragalactic radio sources at optical wavelengths. Celestial reference frames realized in radio wavelengths with the VLBA Calibrator Survey (VCS) data and in optical wavelengths with the Sloan Digital Sky Survey (SDSS) data (DR3 quasar catalogue and DR5) are compared in terms of radio-optical distances between the common sources, global rotation of the axes and offset of the equator.

The new generation of atomic clocks will reach unprecedented uncertainties in frequency of $10^{-18}$. In order to prepare space missions such as ACES, we compute all relativistic frequency shifts detectable during this mission in the case of a clock aboard the International Space Station.