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.

The quasar J1918+4937 was recently suggested to harbour a milliparsec-separation binary supermassive black hole (SMBH), based upon modeling the narrow spike in its high-cadence Kepler optical light curve. Known binary SMBHs are extremely rare, and the tight constraints on the physical and geometric parameters of this object are unique. The high-resolution radio images of J1918+4937 obtained with very long baseline interferometry (VLBI) indicate a rich one-sided jet structure extending to 80 milliarcseconds. Here we analyse simultaneously-made sensitive 1.7- and 5-GHz archive VLBI images as well as snapshot 8.4/8.7-GHz VLBI images of J1918+4937, and show that the appearance of the wiggled jet is consistent with the binary scenario. We develop a jet structural model that handles eccentric orbits. By applying this model to the measured VLBI component positions, we constrain the inclination of the radio jet, as well as the spin angle of the jet emitter SMBH. We find the jet morphological model is consistent with the optical and radio data, and that the secondary SMBH is most likely the jetted one in the system. Furthermore, the decade-long 15-GHz radio flux density monitoring data available for J1918+4937 are compatible with a gradual overall decrease in the the total flux density caused by a slow secular change of the jet inclination due to the spin-orbit precession. J1918+4937 could be an efficient high-energy neutrino source if the horizon of the secondary SMBH is rapidly rotating.

Compact symmetric objects (CSOs) are radio-emitting active galactic nuclei (AGNs) typically with a double-lobed radio structure confined to within 1 kpc. CSOs represent the earliest evolutionary phase of jetted AGNs. Some of them may eventually evolve into large-scale extended double sources, while others stall within the host galaxy and die out, depending on the longevity of nuclear activity, the jet power, and parameters of the surrounding galactic environment. Studying CSOs is a useful tool for understanding the evolution of the galaxies and the interactions between the jets and the medium of the host galaxy. Based on milliarcsec-resolution imaging observations using very long baseline interferometry (VLBI), it is not always straightforward to distinguish between a compact double-lobed or a core-jet structure. The quasar J1110+4817 was considered a CSO candidate in the literature earlier, but because of the lack of clear evidence, it could not be securely classified as a CSO. Here we present a comprehensive analysis of archival multi-frequency VLBI observations combined with accurate Gaia optical astrometric information. Lower-frequency VLBI images reveal an extended radio feature nearly perpendicular to the main structural axis of the source, apparently emanating from the brighter northern feature, that is rare among the known CSOs. While the presence of a binary AGN system cannot be fully excluded, the most plausible explanation is that J1110+4817 is a CSO.

Many low-luminosity active galactic nuclei (AGNs) contain a compact radio core which can be observed with high angular resolution using very long baseline interferometry (VLBI). Combining arcsec-scale structural information with milliarcsec-resolution VLBI imaging is a useful way to characterise the objects and to find compact cores on parsec scales. VLBI imaging could also be employed to look for dual AGNs when the sources show kpc-scale double symmetric structure with flat or inverted radio spectra. We observed five such sources at redshifts 0.36 < z < 0.58 taken from an optically selected sample of Type 2 quasars with the European VLBI Network (EVN) at 1.7 and 5 GHz. Out of the five sources, only one (SDSS J1026-0042) shows a confidently detected compact VLBI core at both frequencies. The other four sources are marginally detected at 1.7 GHz only, indicating resolved-out radio structure and steep spectra. Using first-epoch data from the ongoing Karl G. Jansky Very Large Array Sky Survey, we confirm that indeed all four of these sources have steep radio spectra on arcsec scale, contrary to the inverted spectra reported earlier in the literature. However, the VLBI-detected source, SDSS J1026-0042, has a flat integrated spectrum. Radio AGNs that show kpc-scale symmetric structures with truly flat or inverted spectra could still be promising candidates of dual AGNs, to be targeted with VLBI observations in the future.

The population of intermediate-mass black holes (IMBHs) in nearby dwarf galaxies plays an important "ground truth" role in exploring black hole formation and growth in the early Universe. In the dwarf elliptical galaxy SDSS J090613.77+561015.2 (z=0.0465), an accreting IMBH has been revealed by optical and X-ray observations. Aiming to search for possible radio core and jet associated with the IMBH, we carried out very long baseline interferometry (VLBI) observations with the European VLBI Network (EVN) at 1.66 GHz. Our imaging results show that there are two 1-mJy components with a separation of about 52 mas (projected distance 47 pc) and the more compact component is located within the 1-sigma error circle of the optical centroid from available Gaia astrometry. Based on their positions, elongated structures and relatively high brightness temperatures, as well as the absence of star-forming activity in the host galaxy, we argue that the radio morphology originates from the jet activity powered by the central IMBH. The existence of the large-scale jet implies that violent jet activity might occur in the early epochs of black hole growth and thus help to regulate the co-evolution of black holes and galaxies.

We present the observational results from the 43-GHz Very Long Baseline Array (VLBA) observations of 124 compact radio-loud active galactic nuclei (AGNs) that were conducted between 2014 November and 2016 May. The typical dimensions of the restoring beam in each image are about 0.5 mas $\times$ 0.2 mas. The highest resolution of 0.2 mas corresponds to a physical size of 0.02 pc for the lowest redshift source in the sample. The 43-GHz very long baseline interferometry (VLBI) images of 97 AGNs are presented for the first time. We study the source compactness on milli-arcsec (mas) and sub-mas scales, and suggest that 95 sources in our sample are suitable for future space VLBI observations. By analyzing our data supplemented with other VLBA AGN surveys from literature, we find that the core brightness temperature increases with increasing frequency below a break frequency ~ 7 GHz, and decreases between ~7--240~GHz but increases again above~240 GHz in the rest frame of the sources. This indicates that the synchrotron opacity changes from optically thick to thin. We also find a strong statistical correlation between radio and $\gamma$-ray flux densities. Our correlation is tighter than those in literature derived from lower-frequency VLBI data, suggesting that the $\gamma$-ray emission is produced more co-spatially with the 43-GHz VLBA core emission. This correlation can also be extrapolated to the un-beamed AGN population, implying that a universal $\gamma$-ray production mechanism might be at work for all types of AGNs.

The quest for binary and dual supermassive black holes (SMBHs) at the dawn of the multi-messenger era is compelling. Detecting dual active galactic nuclei (AGN) -- active SMBHs at projected separations larger than several parsecs -- and binary AGN -- probing the scale where SMBHs are bound in a Keplerian binary -- is an observational challenge. The study of AGN pairs (either dual or binary) also represents an overarching theoretical problem in cosmology and astrophysics. The AGN triggering calls for detailed knowledge of the hydrodynamical conditions of gas in the imminent surroundings of the SMBHs and, at the same time, their duality calls for detailed knowledge on how galaxies assemble through major and minor mergers and grow fed by matter along the filaments of the cosmic web. This review describes the techniques used across the electromagnetic spectrum to detect dual and binary AGN candidates and proposes new avenues for their search. The current observational status is compared with the state-of-the-art numerical simulations and models for formation of dual and binary AGN. Binary SMBHs are among the loudest sources of gravitational waves (GWs) in the Universe. The search for a background of GWs at nHz frequencies from inspiralling SMBHs at low redshifts, and the direct detection of signals from their coalescence by the Laser Interferometer Space Antenna in the next decade, make this a theme of major interest for multi-messenger astrophysics. This review discusses the future facilities and observational strategies that are likely to significantly advance this fascinating field.

Blazars are a sub-class of quasars with Doppler boosted jets oriented close to the line of sight, and thus efficient probes of supermassive black hole growth and their environment, especially at high redshifts. Here we report on Very Long Baseline Interferometry observations of a blazar J0906+6930 at z = 5.47, which enabled the detection of polarised emission and measurement of jet proper motion at parsec scales. The observations suggest a less powerful jet compared with the general blazar population, including lower proper motion and bulk Lorentz factor. This coupled with a previously inferred high accretion rate indicate a transition from an accretion radiative power to a jet mechanical power based transfer of energy and momentum to the surrounding gas.While alternative scenarios could not be fully ruled out, our results indicate a possibly nascent jet embedded in and interacting with a dense medium resulting in a jet bending.

High-resolution observations of high-redshift ($z>4$) radio quasars offer a unique insight into jet kinematics at early cosmological epochs, as well as constraints on cosmological model parameters. Due to the general weakness of extremely distant objects and the apparently slow structural changes caused by cosmological time dilation, only a couple of high-redshift quasars have been studied with parsec-scale resolutions, and with limited number of observing epochs. Here we report on very long baseline interferometry (VLBI) observations of a high-redshift blazar J1430+4204 ($z=4.72$) in the 8 GHz frequency band at five different epochs spanning 22 years. The source shows a compact core--jet structure with two jet components being identified within 3 milli-arcsecond (mas) scale. The long time span and multiple-epoch data allow for the kinematic studies of the jet components. That results in a jet proper motion of $\mu {\rm (J1)}$ = 0.017$\pm$0.002~mas\,yr$^{-1}$ and $\mu({\rm J2})$=0.156$\pm$0.015~mas\,yr$^{-1}$, respectively. For the fastest-moving outer jet component J2, the corresponding apparent transverse speed is $19.5 \pm 1.9 \,c$. The inferred bulk jet Lorentz factor $\Gamma = 14.6 \pm 3.8$ and viewing angle $\theta = 2.2^{\circ} \pm 1.6^{\circ}$ indicate highly relativistic beaming. The Lorentz factor and apparent proper motion are the highest measured to date among the $z>4$ jetted radio sources, while the jet kinematics is still consistent with the cosmological interpretation of quasar redshifts.

We applied image stacking on empty-field Faint Images of the Radio Sky at Twenty-Centimeters (FIRST) survey maps centred on optically identified high-redshift quasars at $z\geq4$ to uncover the hidden $\mu$Jy radio emission in these active galactic nuclei (AGN). The median stacking procedure for the full sample of $2229$ optically identified AGN uncovered an unresolved point source with an integrated flux density of 52$~\mu$Jy, with a signal-to-noise ratio $\sim10$. We co-added the individual image centre pixels to estimate the characteristic monochromatic radio power at $1.4~$GHz considering various values for the radio spectral index, revealing a radio population with $P_\mathrm{1.4GHz}\sim10^{24}~$W Hz$^{-1}$. Assuming that the entire radio emission originates from star-forming (SF) activity in the nuclear region of the host galaxy, we obtained an upper limit on the characteristic star formation rate, $\sim4200~$M$_\odot~$yr$^{-1}$. The angular resolution of FIRST images is insufficient to distinguish between the SF and AGN origin of radio emission at these redshifts. However, a comparison with properties of individual sources from the literature indicates that a mixed nature is likely. Future very long baseline interferometry radio observations and ultra-deep Square Kilometre Array surveys are expected to be sensitive enough to detect and resolve the central $1-10~$kpc region in the host galaxies, and thus discriminate between SF and AGN related emission.

In the hierarchical structure formation model, galaxies grow through various merging events. Numerical simulations indicate that the mergers can enhance the activity of the central supermassive black holes in the galaxies. Pfeifle et al. (2019a) identified a system of three interacting galaxies, J0849$+$1114, and provided multiwavelength evidence of all three galaxies containing active galactic nuclei. The system has substantial radio emission, and with high-resolution radio interferometric observation we aimed to investigate its origin, whether it is related to star formation or to one or more of the active galactic nuclei in the system. We performed high-resolution continuum observation of J0849$+$1114 with the European Very Long Baseline Interferometry Network at $1.7$ GHz. We detected one compact radio emitting source at the position of the easternmost nucleus. Its high brightness temperature and radio power indicate that the radio emission originates from a radio-emitting active galactic nucleus. Additionally, we found that significant amount of flux density is contained in $\sim 100$ milliarcsec-scale feature related to the active nucleus.

The largest Key Science Program of the RadioAstron space VLBI mission is a survey of active galactic nuclei (AGN). The main goal of the survey is to measure and study the brightness of AGN cores in order to better understand the physics of their emission while taking interstellar scattering into consideration. In this paper we present detection statistics for observations on ground-space baselines of a complete sample of radio-strong AGN at the wavelengths of 18, 6, and 1.3 cm. Two-thirds of them are indeed detected by RadioAstron and are found to contain extremely compact, tens to hundreds of $\mu$as structures within their cores.

Numerical integration of orbit trajectories for a large number of initial conditions and for long time spans is computationally expensive. Semi-analytical methods were developed to reduce the computational burden. An elegant and widely used method of semi-analytically integrating trajectories of objects subject to atmospheric drag was proposed by King-Hele (KH). However, the analytical KH contraction method relies on the assumption that the atmosphere density decays strictly exponentially with altitude. If the actual density profile does not satisfy the assumption of a fixed scale height, as is the case for Earth's atmosphere, the KH method introduces potentially large errors for non-circular orbit configurations. In this work, the KH method is extended to account for such errors by using a newly introduced atmosphere model derivative. By superimposing exponentially decaying partial atmospheres, the superimposed KH method can be applied accurately while considering more complex density profiles. The KH method is further refined by deriving higher order terms during the series expansion. A variable boundary condition to choose the appropriate eccentricity regime, based on the series truncation errors, is introduced. The accuracy of the extended analytical contraction method is shown to be comparable to numerical Gauss-Legendre quadrature. Propagation using the proposed method compares well against non-averaged integration of the dynamics, while the computational load remains very low.

Recently, Rakshit et al. (2018) reported the discovery of SDSS J103024.95$+$551622.7, a radio-loud narrow-line Seyfert 1 galaxy having a $\sim 100$ kpc scale double-lobed radio structure. Here we analyse archival radio interferometric imaging data taken with the Very Large Array (VLA) at 5 GHz, and with the Very Long Baseline Array (VLBA) at 4.3 and 7.6 GHz. Two hotspots and a compact core are detected with the VLA at arcsec scale, while a single milliarcsec-scale compact radio core is seen with the highest resolution VLBA observations. The Fermi Large Area Telescope did not detect $\gamma$-ray emission at the position of this source. In the mid-infrared, the Wide-field Infrared Survey Explorer satellite light curve, covering more than 7 years and including the most recent data points, hints on flux density variability at 3.4 $\mu$m. Our findings support the notion that this source is a young version of Fanaroff-Riley type II radio galaxies.

The 4MOST instrument is a multi-object spectrograph that will address Galactic and extragalactic science cases simultaneously by observing targets from a large number of different surveys within each science exposure. This parallel mode of operation and the survey nature of 4MOST require some distinct 4MOST-specific operational features within the overall operations model of ESO. The main feature is that the 4MOST Consortium will deliver, not only the instrument, but also contractual services to the user community, which is why 4MOST is also described as a facility. This white paper concentrates on information particularly useful to answering the forthcoming Call for Letters of Intent.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

The binary neutron star merger event GW170817 was detected through both electromagnetic radiation and gravitational waves. Its afterglow emission may have been produced by either a narrow relativistic jet or an isotropic outflow. High spatial resolution measurements of the source size and displacement can discriminate between these scenarios. We present Very Long Baseline Interferometry observations, performed 207.4 days after the merger, using a global network of 32 radio telescopes. The apparent source size is constrained to be smaller than 2.5 milliarcseconds at the 90% confidence level. This excludes the isotropic outflow scenario, which would have produced a larger apparent size, indicating that GW170817 produced a structured relativistic jet. Our rate calculations show that at least 10% of neutron star mergers produce such a jet.

Active galactic nuclei are the most luminous persistent (non-transient, even if often variable) objects in the Universe. They are bright in the entire electromagnetic spectrum. Blazars are a special class where the jets point nearly to our line of sight. Because of this special geometry and the bulk relativistic motion of the plasma in the jet, their radiation is enhanced by relativistic beaming. The majority of extragalactic objects detected in gamma-rays are blazars. However, finding their counterparts in other wavebands could be challenging. Here we present the results of our 5-GHz European VLBI Network (EVN) observation of the radio source J1331+2932, a candidate blazar found while searching for possible gamma-ray emission from the stellar binary system DG CVn (Loh et al. 2017). The highest-resolution radio interferometric measurements provide the ultimate tool to confirm the blazar nature of a radio source by imaging compact radio jet structure with Doppler-boosted radio emission, and give the most accurate celestial coordinates as well.

Apparently periodic optical variations of the luminous high-redshift (z=2.06) quasar PSO J334.2028+1.4075 led Liu et al. (2015) to interpret the variability as the orbital period of a binary supermassive black hole (SMBH) residing in a single circumbinary accretion disk. The proposed orbital separation was around 0.006 pc, and the possible inspiral time about 7 yr in the rest frame of the quasar. Such objects would be of high interest as the difficult-to-find end products of binary SMBH evolution, and potential sources of low-frequency gravitational waves. However, extending the time baseline of the variability study, Liu et al. (2016) later found that the periodicity of PSO J334.2028+1.4075 does not remain persistent. Foord et al. (2017) did not find evidence for the binary active galactic nucleus scenario based on Chandra X-ray observations. The object has also been studied in detail in the radio (Mooley et al. 2018) with the Karl G. Jansky Very Large Array (VLA) and the Very Long Baseline Array (VLBA), revealing a lobe-dominated quasar at kpc scales, and possibly a precessing jet, which might retain PSO J334.2028+1.4075 as a binary SMBH candidate. Here we report on our 1.7-GHz observation with the European VLBI Network (EVN) which complements the VLBA data taken at higher frequencies, and discuss the current knowledge about the nature of this interesting object.

Volonteri et al. (2011) found that the number of radio-loud quasars above redshift 4 calculated from the luminosity function (based upon Swift/BAT observations) is much smaller than the number estimated from the known high-redshift beamed sources, blazars, assuming that for every beamed source with a Lorentz factor of $\Gamma$, statistically $2 \Gamma^2$ non-beamed sources should exist. To explain the missing misaligned (non-beamed) population of high-redshift sources, they proposed various explanations, involving heavy optical obscuration and significantly different Lorentz factors at early cosmological epochs. Our EVN observations targeting high-redshift ($z>4$) blazar candidates revealed 3 sources not showing relativistic beaming, but rather kpc-scale double structures. These three sources have significant radio emission resolved out with the EVN, while they are compact on $\sim 5-10$ arcsec scale. Our dual-frequency ($1.5$ and $5$ GHz) e-MERLIN observations of these three sources revealed a rich morphology, bending jets, and hot spots with possible sites of interaction between the jets and the surrounding medium at intermediate scales.

Previous spectral energy distribution modeling based on XMM-Newton X-ray observation of the classical double-lobed radio galaxy 3C 411 left the possibility open for the presence of a blazar-like core. We investigated this scenario by characterizing the radio brightness distribution in the inner ~10-pc region of the source. We applied the very long baseline interferometry (VLBI) technique at four different frequencies from 1.7 to 7.6 GHz. We analyzed archival data from the Very Long Baseline Array (VLBA) taken in 2014, and data from the European VLBI Network (EVN) obtained in 2017. The VLBI images reveal pc-scale extended structure in the core of 3C 411 that can be modeled with multiple jet components. The measured core brightness temperatures indicate no Doppler enhancement that would be expected from a blazar jet pointing close to the line of sight. While there is no blazar-type core in 3C 411, we found indication of flux density variability. The overall morphology of the source is consistent with a straight jet with ~50{\deg} inclination angle at all scales from pc to kpc.

The search and study of close pairs of supermassive black holes (SMBHs) is important in the study of galaxy mergers which can possibly trigger active galactic nucleus (AGN) activity, and in the context of their evolution into the gravitational wave emitting regime. The quasar SDSS~J1048+0055 was identified as a SMBH binary candidate based on the observed double-peaked \OIII$\lambda\lambda$4959,5007 emission lines and two distinct radio components separated by $\sim 20$~pc \citep{2004ApJ...604L..33Z}. To ascertain the binary nature of this source, we analyzed multi-frequency, multi-epoch very long baseline interferometry (VLBI) data to investigate its pc-scale radio properties. The source shows double components with the western feature being brighter than the eastern one. This brighter component has a brightness temperature of $\sim 10^{10}$~K, spectral index of $\alpha = -0.09 \pm 0.09$ (flat) and is indicative of mildly relativistic beaming. In contrast, the faint component has a lower brightness temperature of $\sim 10^{8-9}$~K and steep spectrum. These clues are consistent with a core--jet structure, moreover, the apparent separation speed between the two components is much higher than the expected orbital motion in a binary SMBH. Thus the present study excludes the association of the two VLBI components with the cores of a SMBH binary, although the SMBH binary possibility (e.g., a pair of radio-loud and radio-quiet AGNs) is not fully ruled out. In the single active galactic nucleus (AGN) scenario, the double-peaked optical emission lines can originate from the jet interacting with the narrow-line region as indicated by a change in the jet direction at $\sim$ 140 pc.

Recently, Saxena et al. (2018) reported the discovery of a possible radio galaxy, J1530$+$1049 at a redshift of z=5.72. We observed the source with the European Very Long Baseline Interferometry Network at $1.7$ GHz. We detected two faint radio features with a separation of $\sim 400$ mas. The radio power calculated from the VLA flux density by Saxena et al. (2018), and the projected source size derived from our EVN data place J1530$+$1049 among the medium-sized symmetric objects (MSOs) which are thought to be young counterparts of radio galaxies (An and Baan 2012). Thus, our finding is consistent with a radio galaxy in an early phase of its evolution as proposed by Saxena et al. (2018).

BL Lac objects are known to have compact jets inclined to our line of sight at a small angle, showing prominent radio emission. Two radio-weak BL Lac candidates with no counterparts in current radio surveys were recently reported by Massaro et al. (2017). Both sources were selected as candidate low-energy counterparts of unassociated Fermi {\gamma}-ray sources. We carried out very long baseline interferometry (VLBI) observations with the European VLBI Network (EVN) at 5 GHz to explore their radio properties at milli-arcsecond (mas) scale. One target, J1410+7405, is clearly detected with the EVN. Its measured 5-GHz flux density, 2.4 mJy, is consistent with recent interferometric measurements with the Karl G. Jansky Very Large Array, suggesting that the radio emission is confined to the inner < ~10 mas region. J1410+7405 is therefore identified as a radio-loud jetted active galactic nucleus, and its brightness temperature exceeds ~10^9 K. Its properties are similar to those of other {\gamma}-ray detected BL Lac objects. On the other hand, the second target, J0644+6031, remains undetected with the EVN with a 6-{\sigma} brightness upper limit of 0.12 mJy/beam. This source is thus radio-quiet, confirming its peculiarity, or possibly questioning its BL Lac nature.

Dual active galactic nuclei (DAGN) and supermassive black hole binaries (SMBHBs) at kpc and pc-scale separations, respectively, are expected during stages of galaxy merger and evolution. Their observational identification can address a range of areas of current astrophysics frontiers including the final parsec problem and their contribution towards the emission of low-frequency gravitational waves. This has however been difficult to achieve with current spectroscopy and time domain strategies. Very long baseline interferometry (VLBI) as a method of directly imaging radio structures with milli-arcsecond (mas) and sub-mas resolutions is introduced as a possible means of detecting DAGN and SMBHBs. We motivate its usage with expected observational signatures and cite some studies from literature to illustrate its current status, and present an updated list of candidates imaged with high-resolution radio observations. We then recall some shortcomings of the method with possible solutions and discuss future directions, relevant to large surveys with the upcoming Square Kilometer Array and future space VLBI missions.

The IceCube instrument detected a high-energy cosmic neutrino event on 2017 September 22 (IceCube_170922A, IceCube Collaboration 2018), which the electromagnetic follow-up campaigns associated with the flaring $\gamma$-ray blazar TXS 0506$+$056 (e.g., Padovani et al., 2018). We investigated the mid-infrared variability of the source by using the available single exposure data of the WISE satellite at $3.4$ and $4.6\mu$m. TXS 0506$+$056 experienced a $\sim 30$% brightening in both of these bands a few days prior to the neutrino event. Additional intraday infrared variability can be detected in 2010. Similar behaviour seen previously in $\gamma$-ray bright radio-loud AGN has been explained by their jet emission (e.g., Jiang et al. 2012).

Context. Radio-loud active galactic nuclei in the early Universe are rare. The quasars J0906+6930 at redshift z=5.47 and J2102+6015 at z=4.57 stand out from the known sample with their compact emission on milliarcsecond (mas) angular scale with high (0.1-Jy level) flux densities measured at GHz radio frequencies. This makes them ideal targets for very long baseline interferometry (VLBI) observations. Aims. By means of VLBI imaging we can reveal the inner radio structure of quasars and model their brightness distribution to better understand the geometry of the jet and the physics of the sources. Methods. We present sensitive high-resolution VLBI images of J0906+6930 and J2102+6015 at two observing frequencies, 2.3 and 8.6 GHz. The data were taken in an astrometric observing programme involving a global five-element radio telescope array. We combined the data from five different epochs from 2017 February to August. Results. For one of the highest redshift blazars known, J0906+6930, we present the first-ever VLBI image obtained at a frequency below 8 GHz. Based on our images at 2.3 and 8.6 GHz, we confirm that this source has a sharply bent helical inner jet structure within ~3 mas from the core. The quasar J2102+6015 shows an elongated radio structure in the east-west direction within the innermost ~2 mas that can be described with a symmetric three-component brightness distribution model at 8.6 GHz. Because of their non-pointlike mas-scale structure, these sources are not ideal as astrometric reference objects. Our results demonstrate that VLBI observing programmes conducted primarily with astrometric or geodetic goals can be utilized for astrophysical purposes as well.

Most of the radio-loud narrow-line Seyfert 1 (NLS1) galaxies resemble compact steep-spectrum sources. However, the extremely radio-loud ones show blazar-like characteristics, like flat radio spectra, compact radio cores, substantial variability and high brightness temperatures. These objects are thought to be similar to blazars as they possess relativistic jets seen at small angle to the line of sight. This claim has been further supported by the Fermi satellite discovery of gamma-ray emission from a handful of these sources. Using the Wide-Field Infrared Survey Explorer (WISE) data, we analyzed the mid-infrared variability characteristics of $42$ radio-loud NLS1 at $3.4$ and $4.6\,\mu$m. We found that $27$ out of the studied $42$ sources showed variability in at least one of the two infrared bands. In some cases, significant changes in the infrared colors can alter the location of the source in the WISE color-color diagram which might lead to different classification. More than $60$% of the variable sources also showed variability within a $1-1.5$ day interval. Such short time scales argue for a compact emission region like those associated with the jets. This connection is further strengthened by the fact that the brightest $\gamma$-ray emitters of the sample ($6$ sources), all showed short time scale infrared variability.

We present a multi-calibrator solution, i.e. MultiView, to achieve accurate astrometry on the level of the thermal noise at low VLBI frequencies dominated by ionospheric residuals. We demonstrate on L-band VLBA observations how MultiView provides superior astrometry to conventional phase referencing techniques (Rioja et al. 2017). We also introduce a new trial method to detect antenna based systematic errors in the observations (Orosz et al. 2017}. All presented methods and results are based on our recent papers (Orosz et al. 2017; Rioja et al. 2017).

To date, PMN J2134-0419 (at a redshift z=4.33) is the second most distant quasar known with a milliarcsecond-scale morphology permitting direct estimates of the jet proper motion. Based on two-epoch observations, we constrained its radio jet proper motion using the very long baseline interferometry (VLBI) technique. The observations were conducted with the European VLBI Network (EVN) at 5 GHz on 1999 November 26 and 2015 October 6. We imaged the central 10-pc scale radio jet emission and modeled its brightness distribution. By identifying a jet component at both epochs separated by 15.86 yr, a proper motion of mu=0.035 +- 0.023 mas/yr is found. It corresponds to an apparent superluminal speed of beta_a=4.1 +- 2.7 c . Relativistic beaming at both epochs suggests that the jet viewing angle with respect to the line of sight is smaller than 20 deg, with a minimum bulk Lorentz factor Gamma=4.3. The small value of the proper motion is in good agreement with the expectations from the cosmological interpretation of the redshift and the current cosmological model. Additionally we analyzed archival Very Large Array observations of J2143-0419 and found indication of a bent jet extending to ~30 kpc.

By discovering a narrow-line emitter embedded in an extended Ly-alpha nebula, Husemann et al. (2018) found that the luminous radio-quiet quasar LBQS 0302-0019 is a dual active galactic nucleus (AGN) system. By combining available archival radio interferometric data taken at 8.4 GHz with the Very Large Array (VLA), we found that the weak 0.4-mJy radio source in this system positionally coincides with the unobscured optical quasar while its newly discovered obscured companion has no detectable radio emission (<0.045 mJy/beam, 3-sigma image noise level). However, the non-detection of this object with the VLA does not rule out the presence of a radio-quiet AGN there.

The Fermi Large Area Telescope revealed that the extragalactic gamma-ray sky is dominated by blazars, active galactic nuclei (AGN) whose jet is seen at very small angle to the line of sight. To associate and then classify the gamma-ray sources, data have been collected from lower frequency surveys and observations. The gamma-ray source 3FGL J1323.0+2942 is associated with the radio source 4C+29.48 and classified as a blazar of unknown type, lacking optical spectrum and redshift. The higher-resolution radio data showed that 4C+29.48 comprises three bright radio-emitting features located within a ~1'-diameter area. We (re-)analyzed archival Very Large Array and unpublished very long baseline interferometry (VLBI) observations conducted by the Very Long Baseline Array and the European VLBI Network of 4C+29.48. We also collected data form optical, infrared and X-ray surveys. The northernmost complex of 4C+29.48 contains a blazar with a high brightness temperature compact core and a steep-spectrum jet feature. The blazar is positionally coincident with an optical source at a redshift of 1.142. Its mid-infrared colors also support its association with a gamma-ray emitting blazar. The two other radio complexes have steep radio spectra and do not have optical or infrared counterparts in currently available surveys. Based on the radio morphology, they are unlikely to be related to the blazar. We discuss the possibilities whether the two radio features are lobes of a radio galaxy, or gravitationally lensed images of a background source. We propose to associate 3FGL J1323.0+2942 in subsequent versions of the Fermi catalog with the blazar residing in northernmost complex. We suggest naming this source J1323+2941A to avoid misinterpretation arising from the fact that the coordinates of the currently listed radio counterpart 4C+29.48 is closer to a most probably unrelated radio source. (abridged)

Radio-loud active galactic nuclei (AGNs), hosting powerful relativistic jet outflows, provide an excellent laboratory for studying jet physics. Very long baseline interferometry (VLBI) enables high-resolution imaging on milli-arcsecond (mas) and sub-mas scales, making it a powerful tool to explore the inner jet structure, shedding light on the formation, acceleration and collimation of AGN jets. In this paper, we present Very Long Baseline Array (VLBA) observations of ten radio-loud AGNs at 43 and 86~GHz, which were selected from the {\it Planck} catalogue of compact sources and are among the brightest in published VLBI images at and below 15 GHz. The image noise levels in our observations are typically 0.3 mJy beam$^{-1}$ and 1.5 mJy beam$^{-1}$ at 43 and 86 GHz, respectively. Compared with the VLBI data observed at lower frequencies from the literature, our observations with higher resolution (the highest resolution up to 0.07 mas at 86 GHz and 0.18 mas at 43 GHz) and at higher frequencies detected new jet components at sub-parsec scales, offering valuable data for studies of the physical properties of innermost jets. These include compactness factor of the radio structure (the ratio of core flux density to total flux density), and core brightness temperature ($T_{\rm b}$). In all these sources, the compact core accounts for a significant fraction ($> 60\%$) of the total flux density. Their correlated flux density at the longest baselines is higher than 0.16 Jy. The compactness of these sources make them good phase calibrators of mm-wavelength ground-based and space VLBI.

Narrow-line Seyfert 1 galaxies (NLS1) are an intriguing subclass of active galactic nuclei. Their observed properties indicate low central black hole mass and high accretion rate. The extremely radio-loud NLS1 sources often show relativistic beaming and are usually regarded as younger counterparts of blazars. Recently, the object SDSS J110006.07+442144.3 was reported as a candidate NLS1 source. The characteristics of its dramatic optical flare indicated its jet-related origin. The spectral energy distribution of the object was similar to that of the gamma-ray detected radio-loud NLS1, PMN J0948+0022. Our high-resolution European Very Long Baseline Interferometry Network observations at 1.7 and 5 GHz revealed a compact core feature with a brightness temperature of >~ 10^(10) K. Using the lowest brightness temperature value and assuming a moderate Lorentz factor of ~9 the jet viewing angle is <~ 26 deg. Archival Very Large Array data show a large-scale radio structure with a projected linear size of ~150 kpc reminiscent of double-sided morphology.

Radio-loud high-redshift quasars (HRQs), although only a few of them are known to date, are crucial for the studies of the growth of supermassive black holes (SMBHs) and the evolution of active galactic nuclei (AGN) at early cosmological epochs. Radio jets offer direct evidence of SMBHs, and their radio structures can be studied with the highest angular resolution using Very Long Baseline Interferometry (VLBI). Here we report on the observations of three HRQs (J0131-0321, J0906+6930, J1026+2542) at z>5 using the Korean VLBI Network and VLBI Exploration of Radio Astrometry Arrays (together known as KaVA) with the purpose of studying their pc-scale jet properties. The observations were carried out at 22 and 43 GHz in 2016 January among the first-batch open-use experiments of KaVA. The quasar J0906+6930 was detected at 22 GHz but not at 43 GHz. The other two sources were not detected and upper limits to their compact radio emission are given. Archival VLBI imaging data and single-dish 15-GHz monitoring light curve of J0906+6930 were also acquired as complementary information. J0906+6930 shows a moderate-level variability at 15 GHz. The radio image is characterized by a core-jet structure with a total detectable size of ~5 pc in projection. The brightness temperature, 1.9x10^{11} K, indicates relativistic beaming of the jet. The radio properties of J0906+6930 are consistent with a blazar. Follow-up VLBI observations will be helpful for determining its structural variation.

High-redshift quasars are important to study galaxy and active galactic nuclei (AGN) evolution, test cosmological models, and study supermassive black hole growth. Optical searches for high-redshift sources have been very successful, but radio searches are not hampered by dust obscuration and should be more effective at finding sources at even higher redshifts. Identifying high-redshift sources based on radio data is, however, not trivial. Here we report on new multi-frequency Giant Metrewave Radio Telescope (GMRT) observations of eight z>4.5 sources previously studied at high angular resolution with very long baseline interferometry (VLBI). Combining these observations with those from the literature, we construct broad-band radio spectra of all 30 z>4.5 sources that have been observed with VLBI. In the sample we found flat, steep and peaked spectra in approximately equal proportions. Despite several selection effects, we conclude that the z>4.5 VLBI (and likely also non-VLBI) sources have diverse spectra and that only about a quarter of the sources in the sample have flat spectra. Previously, the majority of high-redshift radio sources were identified based on their ultra-steep spectra (USS). Recently a new method has been proposed to identify these objects based on their megahertz-peaked spectra (MPS). Neither method would have identified more than 18% of the high-redshift sources in this sample. More effective methods are necessary to reliably identify complete samples of high-redshift sources based on radio data.

We report on the measurement of the trigonometric parallaxes of 1612 MHz hydroxyl masers around two asymptotic giant branch stars, WX Psc and OH138.0+7.2, using the NRAO Very Long Baseline Array with in-beam phase referencing calibration. We obtained a 3-sigma upper limit of <=5.3 mas on the parallax of WX Psc, corresponding to a lower limit distance estimate of >~190 pc. The obtained parallax of OH138.0+7.2 is 0.52+/-0.09 mas (+/-18%), corresponding to a distance of 1.9(+0.4,-0.3) kpc, making this the first hydroxyl maser parallax below one milliarcsecond. We also introduce a new method of error analysis for detecting systematic errors in the astrometry. Finally, we compare our trigonometric distances to published phase-lag distances toward these stars and find a good agreement between the two methods.

In hierarchical structure formation scenarios, merging galaxies are expected to be seen in different phases of their coalescence. Simulations suggest that simultaneous activity of the supermassive black holes (SMBHs) in the centres of the merging galaxies may be expected at kpc-scale separations. Currently, there are no direct observational methods which allow the selection of a large number of such dual active galactic nuclei (AGN) candidates. SDSS J142507.32+323137.4 was reported as a promising candidate source based on its optical spectrum. Here we report on our sensitive e-MERLIN observations performed at $1.6$ and at $5$ GHz, which show that the optical spectrum of the source can be more straightforwardly explained with jet-cloud interactions instead of the dual AGN scenario.

Blazars are active galactic nuclei (AGN) whose relativistic jets point nearly to the line of sight. Their compact radio structure can be imaged with very long baseline interferometry (VLBI) on parsec scales. Blazars at extremely high redshifts provide a unique insight into the AGN phenomena in the early Universe. We observed four radio sources at redshift $z>4$ with the European VLBI Network (EVN) at 1.7 and 5 GHz. These objects were previously classified as blazar candidates based on X-ray observations. One of them, J2134$-$0419 is firmly confirmed as a blazar with our VLBI observations, due to its relativistically beamed radio emission. Its radio jet extended to $\sim$10 milli-arcsec scale makes this source a promising target for follow-up VLBI observations to reveal any apparent proper motion. Another target, J0839+5112 shows a compact radio structure typical of quasars. There is evidence for flux density variability and its radio "core" has a flat spectrum. However, the EVN data suggest that its emission is not Doppler-boosted. The remaining two blazar candidates (J1420+1205 and J2220+0025) show radio properties totally unexpected from radio AGN with small-inclination jet. Their emission extends to arcsec scales and the Doppler factors of the central components are well below 1. Their structures resemble that of double-lobed radio AGN with large inclination to the line of sight. This is in contrast with the blazar-type modeling of their multi-band spectral energy distributions. Our work underlines the importance of high-resolution VLBI imaging in confirming the blazar nature of high-redshift radio sources.

The observation and imaging of hundreds or thousands of radio sources with the technique of very long baseline interferometry (VLBI) is a computationally intensive task. However, these surveys allow us to conduct statistical investigations of large source samples, and also to discover new phenomena or types of objects. The field of view of these high-resolution VLBI imaging observations is typically a few arcseconds at cm wavelengths. For practical reasons, often a much smaller fraction of the field, the central region is imaged only. With an automated process we imaged the ~1.5-arcsec radius fields around more than 1000 radio sources, and found a variety of extended radio structures. Some of them are yet unknown in the literature.

The arrival of the Square Kilometer Array (SKA) will revitalise all aspects of Very Long Baseline Interferometry (VLBI) astronomy at the lower frequencies. In the last decade there have been huge strides towards routinely achieving high precision VLBI astrometry at frequencies dominated by the tropospheric contributions, most notably at 22GHz, using advanced phase referencing techniques. Nevertheless to increase the capability for high precision astrometric measurements at low radio frequencies (<8GHz) an effective calibration strategy of the systematic ionospheric propagation effects that is widely applicable is required. Observations at low frequencies are dominated by distinct direction dependent ionospheric propagation errors, which place a very tight limit on the angular separation of a suitable phase referencing calibrator. The MultiView technique holds the key to the compensation of atmospheric spatial-structure errors, by using observations of multiple calibrators and 2-D interpolation. In this paper we present the first demonstration of the power of MultiView using three calibrators, several degrees from the target, along with a comparative study of the astrometric accuracy between MultiView and phase-referencing techniques. MultiView calibration provides an order of magnitude improvement in astrometry with respect to conventional phase referencing, achieving ca. 100 micro-arcseconds astrometry errors in a single epoch of observations, effectively reaching the thermal noise limit. MultiView will achieve its full potential with the enhanced sensitivity and multibeam capabilities of SKA and the pathfinders, which will enable simultaneous observations of the target and calibrators. Our demonstration indicates that the 10 micro-arcseconds goal of astrometry at ~1.6GHz using VLBI with SKA is feasible using the MultiView technique.

Most of the gamma-ray emitting active galactic nuclei (AGN) are blazars, although there is still a small fraction of non-blazar AGN in the Fermi/LAT catalog. Among these misaligned gamma-ray-emitting AGN, a few can be classified as Compact Symmetric Objects (CSOs). In contrast to blazars in which gamma-ray emission is generally thought to originate from highly beamed relativistic jets, the source of gamma-ray emission in unbeamed CSOs remains an open question. The rarity of the gamma-ray emitting CSOs is a mystery as well. Here we present the radio properties of the gamma-ray CSO candidate 2234+282.

The prominent flat-spectrum radio quasar J0017+8135 (S5 0014+81) at z = 3.366 is one of the most luminous active galactic nuclei (AGN) known. Its milliarcsecond-scale radio jet structure has been studied with very long baseline interferometry (VLBI) since the 1980s. The quasar was selected as one of the original defining objects of the International Celestial Reference Frame, but left out from its current second realization (ICRF2) because of systematic long-term positional variations. Here we analyse archival 8.6- and 2.3-GHz VLBI imaging data collected at nearly 100 different epochs during more than 20 years, to obtain information about the kinematics of jet components. Because of the cosmological time dilation, extensive VLBI monitoring data are essential to reveal changes in the jet structure of high-redshift AGN. In the case of J0017+8135, the data can be described with a simple kinematic model of jet precession with a 12-year periodicity in the observer's frame.

High-redshift radio-loud quasars are used to, among other things, test the predictions of cosmological models, set constraints on black hole growth in the early universe and understand galaxy evolution. Prior to this paper, 20 extragalactic radio sources at redshifts above 4.5 have been imaged with very long baseline interferometry (VLBI). Here we report on observations of an additional ten z>4.5 sources at 1.7 and 5 GHz with the European VLBI Network (EVN), thereby increasing the number of imaged sources by 50%. Combining our newly observed sources with those from the literature, we create a substantial sample of 30 z>4.5 VLBI sources, allowing us to study the nature of these objects. Using spectral indices, variability and brightness temperatures, we conclude that of the 27 sources with sufficient information to classify, the radio emission from one source is from star formation, 13 are flat-spectrum radio quasars and 13 are steep-spectrum sources. We also argue that the steep-spectrum sources are off-axis (unbeamed) radio sources with rest-frame self-absorption peaks at or below GHz frequencies and that these sources can be classified as gigahertz peaked-spectrum (GPS) and megahertz peaked-spectrum (MPS) sources.

The quasar PG 1302$-$102 is believed to harbour a supermassive binary black hole (SMBBH) system. Using the available 15 GHz and $2-8$ GHz, multi-epoch Very Long Baseline Array data, we constrain the pc-scale jet properties based on the inferred mean proper motion, including a bulk Lorentz factor $\geq 5.1 \pm 0.8$, jet inclination angle $\leq (11.4 \pm 1.7)$ degrees, projected position angle $= 31.8$ degrees, intrinsic half opening angle $\leq (0.9 \pm 0.1)$ degrees and a mean $2-8$ GHz spectral index of 0.31. A general relativistic helical jet model is presented and applied to predict quasi-periodic oscillations of $\sim$ 10 days, power law power spectrum shape and a contribution of up to $\sim$ 53 percent to the observed variable core flux density. The model is used to make a case for high resolution, moderately sampled, long duration radio interferometric observations to reveal signatures due to helical knots and distinguish them from those due to SMBBH orbital activity including a phase difference $\sim \pi$ and an amplitude ratio (helical light curve amplitude/SMBBH light curve amplitude) of $0.2-3.3$. The prescription can be used to identify helical kinematic signatures from quasars, providing possible candidates for further studies with polarization measurements. It can also be used to infer promising SMBBH candidates for the study of gravitational waves if there are systematic deviations from helical signatures.

A bright optical flare was detected in the high-redshift ($z=2.133$) quasar CGRaBS J0809+5341 on 2014 April 13. The absolute magnitude of the object reached $-30.0$ during the flare, making it the brightest one (in flaring stage) among all known quasars so far. The 15 GHz flux density of CGRaBS J0809+5341 monitored in the period from 2008 to 2016 also reached its peak at the same time. To reveal any structural change possibly associated with the flare in the innermost radio structure of the quasar, we conducted a pilot very long baseline interferometry (VLBI) observation of CGRaBS J0809+5341 using the European VLBI Network (EVN) at 5 GHz on 2014 November 18, about seven months after the prominent optical flare. Three epochs of follow-up KaVA (Korean VLBI Network and VLBI Exploration of Radio Astrometry Array) observations were carried out at 22 and 43 GHz frequencies from 2015 February 25 to June 4, with the intention of exploring a possibly emerging new radio jet component associated with the optical flare. However, these high-resolution VLBI observations revealed only the milliarcsecond-scale compact "core" that was known in the quasar from earlier VLBI images, and showed no sign of any extended jet structure. Neither the size, nor the flux density of the "core" changed considerably after the flare according to our VLBI monitoring. The results suggest that any putative radio ejecta associated with the major optical and radio flare could not yet be separated from the "core" component, or the newly-born jet was short-lived.

According to hierarchical structure formation models, merging galaxies are expected to be seen in different stages of their coalescence. However, currently there are no straightforward observational methods neither to select nor to confirm a large number of dual active galactic nuclei (AGN) candidates. Most attempts involve the better understanding of double-peaked narrow emission line sources, to distinguish the objects where the emission lines originate from narrow-line kinematics or jet-driven outflows from those which might harbour dual AGN. We observed four such candidate sources with the Very Long Baseline Array (VLBA) at 1.5 GHz with $\sim$ 10 milli-arcsecond angular resolution where spectral profiles of AGN optical emission suggested the existence of dual AGN. In SDSS J210449.13-000919.1 and SDSS J23044.82-093345.3, the radio structures are aligned with the optical emission features, thus the double-peaked emission lines might be the results of jet-driven outflows. In the third detected source SDSS J115523.74+150756.9, the radio structure is less extended and oriented nearly perpendicular to the position angle derived from optical spectroscopy. The fourth source remained undetected with the VLBA but it has been imaged with the Very Large Array at arcsec resolution a few months before our observations, suggesting the existence of extended radio structure. In none of the four sources did we detect two radio-emitting cores, a convincing signature of duality.

Megahertz peaked-spectrum (MPS) sources have spectra that peak at frequencies below 1 GHz in the observer's frame and are believed to be radio-loud active galactic nuclei (AGN). We recently presented a new method to search for high-redshift AGN by identifying unusually compact MPS sources. In this paper, we present European VLBI Network (EVN) observations of 11 MPS sources which we use to determine their sizes and investigate the nature of the sources with ~10 mas resolution. Of the 11 sources, we detect nine with the EVN. Combining the EVN observations with spectral and redshift information, we show that the detected sources are all AGN with linear sizes smaller than 1.1 kpc and are likely young. This shows that low-frequency colour-colour diagrams are an easy and efficient way of selecting small AGN and explains our high detection fraction (82%) in comparison to comparable surveys. Finally we argue that the detected sources are all likely compact symmetric objects and that none of the sources are blazars.

HESS J1943+213 is an unidentified TeV source that is likely a high-frequency-peaked BL Lac (HBL) object but also compatible with a pulsar wind nebula (PWN) nature. Each of these enormously different astronomical interpretations is supported by some of the observed unusual characteristics. In order to finally classify and understand this object we took a three-pronged approach, through time-domain, high angular resolution, and multi-frequency radio studies. First, our deep time-domain observations with the Arecibo telescope failed to uncover the putative pulsar powering the proposed PWN. We conclude with ~70% certainty that HESS J1943+213 does not host a pulsar. Second, long-baseline interferometry of the source with e-MERLIN at 1.5- and 5- GHz, shows only a core, a point source at ~ 1 - 100 milli-arcsecond resolution. Its 2013 flux density is about one-third lower than detected in 2011 observations with similar resolution. This radio variability of the core strengthens the HBL object hypothesis. More evidence against the PWN scenario comes, third, from the radio spectrum we compiled. The extended structure follows a power-law behavior with spectral index alpha = -0.54 +- 0.04 while the core component is flat spectrum (alpha = -0.03 +- 0.03). In contrast, the radio synchrotron emission of PWNe predicts a single power-law distribution. Overall we rule out the PWN hypothesis and conclude the source is a BL Lac object. The consistently high fraction (70%) of the flux density from the extended structure then leads us to conclude that HESS J1943+213 must be a non-classical HBL object.

Hot dust-obscured galaxies (hot DOGs) are a rare class of hyperluminous infrared galaxies identified with the Wide-field Infrared Survey Explorer (WISE) satellite. The majority of them is at high redshifts (z~2-3), at the peak epoch of star formation in the Universe. Infrared, optical, radio, and X-ray data suggest that hot DOGs contain heavily obscured, extremely luminous active galactic nuclei (AGN). This class may represent a short phase in the life of the galaxies, signifying the transition from starburst- to AGN-dominated phases. Hot DOGs are typically radio-quiet, but some of them show mJy-level emission in the radio (microwave) band. We observed four hot DOGs using the technique of very long baseline interferometry (VLBI). The 1.7-GHz observations with the European VLBI Network (EVN) revealed weak radio features in all sources. The radio is free from dust obscuration and, at such high redshifts, VLBI is sensitive only to compact structures that are characteristic of AGN activity. In two cases (WISE J0757+5113, WISE J1603+2745), the flux density of the VLBI-detected components is much smaller than the total flux density, suggesting that ~70-90 per cent of the radio emission, while still dominated by AGN, originates from angular scales larger than probed by the EVN. The source WISE J1146+4129 appears a candidate compact symmetric object, and WISE J1814+3412 shows a 5.1-kpc double structure, reminiscent of hot spots in a medium-sized symmetric object. Our observations support that AGN residing in hot DOGs may be genuine young radio sources where starburst and AGN activities coexist.

We investigate the pc-scale kinematics and kpc-scale radio morphology of the quasar PG 1302-102, which may harbour a sub-pc separation supermassive binary black hole system at its centre as inferred from optical variability. High-resolution radio interferometric measurements obtained with the Very Long Baseline Array (VLBA) in the Monitoring Of Jets in Active galactic nuclei with VLBA Experiments (MOJAVE) programme at 15 GHz at 20 epochs spanning 17 years were analysed to investigate the pc-scale radio structure. Archival observations with the Very Large Array (VLA) at 1.4 GHz and 5 GHz were obtained to study the kpc-scale morphology. We find that the pc-scale jet is inclined within ~2.2 deg to the line of sight and has a half-opening angle of about 0.2 deg. The parameters derived from the pc-scale radio jet are qualitatively consistent with those obtained from the analysis of the optical light curve of PG 1302-102. We obtain at least 0.08 for the mass ratio of the two black holes in the system. We find some indication for a helical jet structure on kpc-scale, but the directions of the inner and the extended radio jets are significantly different, obstructing a straightforward connection of the pc- and kpc-scale jets within the binary scenario.

We present quasi-simultaneous, multi-epoch radio and X-ray measurements of Holmberg II X-1 using the European VLBI Network (EVN), the Karl G. Jansky Very Large Array (VLA), and the Chandra and Swift X-ray telescopes. The X-ray data show apparently hard spectra with steady X-ray luminosities 4 months apart from each other. In the high-resolution EVN radio observations, we have detected an extended milli-arcsecond scale source with unboosted radio emission. The source emits non-thermal, likely optically thin synchrotron emission and its morphology is consistent with a jet ejection. The 9-GHz VLA data show an arcsecond-scale triple structure of Holmberg II X-1 similar to that seen at lower frequencies. However, we find that the central ejection has faded by at least a factor of 7.3 over 1.5 years. We estimate the dynamical age of the ejection to be higher than 2.1 years. We show that such a rapid cooling can be explained with simple adiabatic expansion losses. These properties of Holmberg II X-1 imply that ULX radio bubbles may be inflated by ejecta instead of self-absorbed compact jets.

HESS J1943+213 is a TeV source close to the Galactic plane proposed to be a BL Lac object. Our high resolution EVN observation failed to recover two thirds of the source flux density detected simultaneously by the WSRT. Our recent e-MERLIN observations in L and C bands show only a point source with flux density comparable to the EVN detection. Thus the structure responsible for the missing flux density has to be larger than 2". It may be related to the presumed extragalactic source (thus would have a kpc-scale size), or to the Galactic foreground material close to the line of sight to the source.

Hot dust-obscured galaxies (hot DOGs) are a rare class of hyperluminous infrared galaxies recently identified with the Wide-field Infrared Survey Explorer (WISE) satellite. The majority of the ~1000-member all-sky population should be at high redshifts (z~2-3), at the peak of star formation in the history of the Universe. This class most likely represents a short phase during galaxy merging and evolution, a transition from starburst- to AGN-dominated phases. For the first time, we observed four hot DOGs with known mJy-level radio emission using the European VLBI Network (EVN) at 1.7 GHz, in a hope to find compact radio features characteristic to AGN activity. All four target sources are detected at ~15-30 mas angular resolution, confirming the presence of an active nucleus. The sources are spatially resolved, i.e. the flux density of the VLBI-detected components is smaller than the total flux density, suggesting that a fraction of the radio emission originates from larger-scale (partly starburst-related) activity. Here we show the preliminary results of our e-EVN observations made in 2014 February, and discuss WISE J1814+3412, an object with kpc-scale symmetric radio structure, in more detail.

Few high-redshift, radio-loud quasars are known to date. The extremely luminous, radio-bright quasar, SDSS J013127.34-032100.1 was recently discovered at a redshift of $z=5.18$. We observed the source with high resolution very long baseline interferometry (VLBI) at 1.7 GHz with the European VLBI Network (EVN) and found a single compact radio component. We estimated a lower limit to the brightness temperature of the detected radio component, T_B~10^{11} K. Additionaly, when compared to archival radio data, the source showed significant flux density variation. These two findings are indicative of the blazar nature of the source.

Galaxies and supermassive black holes (SMBHs) are believed to evolve through a process of hierarchical merging and accretion. Through this paradigm, multiple SMBH systems are expected to be relatively common in the Universe. However, to date there are poor observational constraints on multiple SMBHs systems with separations comparable to a SMBH gravitational sphere of influence (<< 1 kpc). In this chapter, we discuss how deep continuum observations with the SKA will make leading contributions towards understanding how multiple black hole systems impact galaxy evolution. In addition, these observations will provide constraints on and an understanding of stochastic gravitational wave background detections in the pulsar timing array sensitivity band (nanoHz -microHz). We also discuss how targets for pointed gravitational wave experiments (that cannot be resolved by VLBI) could potentially be found using the large-scale radio-jet morphology, which can be modulated by the presence of a close-pair binary SMBH system. The combination of direct imaging at high angular resolution; low-surface brightness radio-jet tracers; and pulsar timing arrays will allow the SKA to trace black hole binary evolution from separations of a galaxy virial radius down to the sub-parsec level. This large dynamic range in binary SMBH separation will ensure that the SKA plays a leading role in this observational frontier.

Adding VLBI capability to the SKA arrays will greatly broaden the science of the SKA, and is feasible within the current specifications. SKA-VLBI can be initially implemented by providing phased-array outputs for SKA1-MID and SKA1-SUR and using these extremely sensitive stations with other radio telescopes, and in SKA2 by realising a distributed configuration providing baselines up to thousands of km, merging it with existing VLBI networks. The motivation for and the possible realization of SKA-VLBI is described in this paper.

A brief review of VLBI structures in extremely high-redshift AGN (revised/edited).

The extremely high redshift (z=5.3) radio source SDSS J102623.61+254259.5 (J1026+2542) is among the most distant and most luminous radio-loud active galactic nuclei (AGN) known to date. Its one-sided radio jet structure on milli-arcsecond (mas) and ~10-mas scales typical for blazars was first imaged at 5 GHz with very long baseline interferometry (VLBI) in 2006. Here we report on our dual-frequency (1.7 and 5 GHz) imaging observations performed with the European VLBI Network (EVN) in 2013. The prominent jet structure allows us to identify individual components whose apparent displacement can be detected over the time span of 7.33 yr. This is the first time when jet proper motions are directly derived in a blazar at z>5. The small values of up to ~0.1 mas/yr are consistent with what is expected in a relativistic cosmological model if redshift is a measure of distance. The apparent superluminal jet speeds, considered tentative because derived from two epochs only, exceed 10c for three different components along the jet. Based on modeling its spectral energy distribution, J1026+2542 is known to have its X-ray jet oriented close to the line of sight, with significant Doppler boosting and a large bulk Lorentz factor (Gamma~13). The new VLBI observations, indicating ~2.3 x 10^12 K lower limit to the core brightness temperature, are consistent with this picture. The spectral index in the core region is -0.35.

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants

Galaxies are believed to evolve through merging, which should lead to multiple supermassive black holes in some. There are four known triple black hole systems, with the closest pair being 2.4 kiloparsecs apart (the third component is more distant at 3 kiloparsecs), which is far from the gravitational sphere of influence of a black hole with mass $\sim$10$^9$ M$_\odot$ (about 100 parsecs). Previous searches for compact black hole systems concluded that they were rare, with the tightest binary system having a separation of 7 parsecs. Here we report observations of a triple black hole system at redshift z=0.39, with the closest pair separated by $\sim$140 parsecs. The presence of the tight pair is imprinted onto the properties of the large-scale radio jets, as a rotationally-symmetric helical modulation, which provides a useful way to search for other tight pairs without needing extremely high resolution observations. As we found this tight pair after searching only six galaxies, we conclude that tight pairs are more common than hitherto believed, which is an important observational constraint for low-frequency gravitational wave experiments.

The Seyfert galaxy NGC 5515 has double-peaked narrow-line emission in its optical spectrum, and it has been suggested that this could indicate that it has two active nuclei. We observed the source with high resolution Very Long Baseline Interferometry (VLBI) at two radio frequencies, reduced archival Very Large Array data, and re-analysed its optical spectrum. We detected a single, compact radio source at the position of NGC 5515, with no additional radio emission in its vicinity. The optical spectrum of the source shows that the blue and red components of the double-peaked lines have very similar characteristics. While we cannot rule out unambiguously that NGC 5515 harbours a dual AGN, the assumption of a single AGN provides a more plausible explanation for the radio observations and the optical spectrum.

Recently Nyland et al. (2012) argued that the radio emission observed in the center of the dwarf galaxy NGC 404 originates in a low-luminosity active galactic nucleus (LLAGN) powered by a massive black hole ($M\sim<10^6$ M$_{\odot}$). High-resolution radio detections of MBHs are rare. Here we present sensitive, contemporaneous Chandra X-ray, and very long baseline interferometry (VLBI) radio observations with the European VLBI Network (EVN). The source is detected in the X-rays, and shows no long-term variability. If the hard X-ray source is powered by accretion, the apparent low accretion efficiency would be consistent with a black hole in the hard state. Hard state black holes are known to show radio emission compact on the milliarcsecond scales. However, the central region of NGC 404 is resolved out on 10 milliarcsecond (0.15-1.5 pc) scales. Our VLBI non-detection of a compact, partially self-absorbed radio core in NGC 404 implies that either the black hole mass is smaller than $3^{+5}_{-2}\times10^5$ M$_{\odot}$, or the source does not follow the fundamental plane of black hole activity relation. An alternative explanation is that the central black hole is not in the hard state. The radio emission observed on arcsecond (tens of pc) scales may originate in nuclear star formation or extended emission due to AGN activity, although the latter would not be typical considering the structural properties of low-ionization nuclear emission-line region galaxies (LINERs) with confirmed nuclear activity.

We investigated the archival ground-based VLBI images of the extragalactic radio sources included in both the {\it Wilkinson Microwave Anisotropy Probe (WMAP)} and the {\it Planck} catalogues, and selected 49 bright and compact sources as potential targets for space Very Long Baseline Interferometry (VLBI) observations at mm wavelengths. These sources have a flat radio continuum spectrum between 33 and 94~GHz. They are identified as core-dominated active galactic nuclei (AGN), located at declinations above $-40\degr$, and have never been observed with ground-based VLBI at 86~GHz. The radio properties of the 49 new sources are presented. We compare this new sample with similar samples of compact AGN available from earlier studies. The new candidates, together with the existing bright compact AGN sample identified from 86-GHz ground-based VLBI imaging surveys, form a catalogue of more than 160 AGN. These could be primary targets for mm-VLBI observations on the ground, as well as for future mm-wavelength space VLBI missions such as the project with two satellites currently under study in China.

A patch of sky in the SDSS Stripe 82 was observed at 1.6 GHz with Very Long Baseline Interferometry (VLBI) using the European VLBI Network (EVN). The data were correlated at the EVN software correlator at JIVE (SFXC). There are fifteen known mJy/sub-mJy radio sources in the target field defined by the primary beam size of a typical 30-m class EVN radio telescope. The source of particular interest is a recently identified high-redshift radio quasar J222843.54+011032.2 (J2228+0110) at redshift z=5.95. Our aim was to investigate the milli-arcsecond (mas) scale properties of all the VLBI-detectable sources within this primary beam area with a diameter of 20 arcmin. The source J2228+0110 was detected with VLBI with a brightness temperature T_b>10^8 K, supporting the active galactic nucleus (AGN) origin of its radio emission, which is conclusive evidence that the source is a radio quasar. In addition, two other target sources were also detected, one of them with no redshift information. Their brightness temperature values (T_b >10^7 K) measured with VLBI suggest a non-thermal synchrotron radiation origin for their radio emission. The detection rate of 20% is broadly consistent with other wide-field VLBI experiments carried out recently. We also derived the accurate equatorial coordinates of the three detected sources using the phase-referencing technique. This experiment is an early attempt of a wide-field science project with SFXC, paving the way for the EVN to conduct a large-scale VLBI survey in the multiple-phase-centre mode.

We report on the discovery of an apparent triple radio structure hidden inside the radio bubble of the ultraluminous X-ray source Holmberg II X-1. The morphology is consistent with a collimated jet structure, which is observed to emit optically thin synchrotron radiation. The central component has a steep radio spectrum and is brighter than the outer components indicating a renewed radio activity. We estimate a minimum time-averaged jet power of 2 x 10^{39} erg/s that is associated with a time-averaged isotropic X-ray luminosity of at least 4 x 10^{39} erg/s. Our results suggest that Holmberg II X-1 is powered by a black hole of M_BH \geq 25 M_sun, that is inferred to be accreting at a high Eddington rate with intermittent radio activity.

An international consortium is presently constructing a beamformer for the Atacama Large Millimeter/submillimeter Array (ALMA) in Chile that will be available as a facility instrument. The beamformer will aggregate the entire collecting area of the array into a single, very large aperture. The extraordinary sensitivity of phased ALMA, combined with the extremely fine angular resolution available on baselines to the Northern Hemisphere, will enable transformational new very long baseline interferometry (VLBI) observations in Bands 6 and 7 (1.3 and 0.8 mm) and provide substantial improvements to existing VLBI arrays in Bands 1 and 3 (7 and 3 mm). The ALMA beamformer will have impact on a variety of scientific topics, including accretion and outflow processes around black holes in active galactic nuclei (AGN), tests of general relativity near black holes, jet launch and collimation from AGN and microquasars, pulsar and magnetar emission processes, the chemical history of the universe and the evolution of fundamental constants across cosmic time, maser science, and astrometry.

The galaxy 3C\,316 is the brightest in the radio band among the optically-selected candidates exhibiting double-peaked narrow optical emission lines. Observations with the Very Large Array (VLA), Multi-Element Remotely Linked Interferometer Network (e-MERLIN), and the European VLBI Network (EVN) at 5\,GHz have been used to study the radio structure of the source in order to determine the nature of the nuclear components and to determine the presence of radio cores. The e-MERLIN image of 3C 316 reveals a collimated coherent east-west emission structure with a total extent of about 3 kpc. The EVN image shows seven discrete compact knots on an S-shaped line. However, none of these knots could be unambiguously identified as an AGN core. The observations suggest that the majority of the radio structure belongs to a powerful radio AGN, whose physical size and radio spectrum classify it as a compact steep-spectrum source. Given the complex radio structure with radio blobs and knots, the possibility of a kpc-separation dual AGN cannot be excluded if the secondary is either a naked core or radio quiet.

Context. It will soon become possible to directly link the most accurate radio reference frame with the Gaia optical reference frame using many common extragalactic objects. It is important to know the level of coincidence between the radio and optical positions of compact active galactic nuclei (AGN). Aims. Using the best catalogues available at present, we investigate how many AGN with significantly large optical-radio positional offsets exist as well as the possible causes of these offsets. Methods. We performed a case study by finding optical counterparts to the International Celestial Reference Frame (ICRF2) radio sources in the Sloan Digital Sky Survey (SDSS) Data Release 9 (DR9). The ICRF2 catalogue was used as a reference because the radio positions determined by Very Long Baseline Interferometry (VLBI) observations are about two orders of magnitude more accurate than the optical positions. Results. We find 1297 objects in common for ICRF2 and SDSS DR9. Statistical analysis of the optical-radio differences verifies that the SDSS DR9 positions are accurate to ~55 mas in both coordinates, with no systematic offset with respect to ICRF2. We find 51 sources (~4% of the sample) for which the positional offset exceeds 170 mas (~3{\sigma}). Astrophysical explanations must exist for most of these outliers. There are 3 known strong gravitational lenses among them. Dual AGN or recoiling supermassive black holes may also be possible. Conclusions. The most accurate Gaia-VLBI reference frame link will require a careful selection of a common set of objects by eliminating the outliers. On the other hand, the significant optical-radio positional non-coincidences may offer a new tool for finding e.g. gravitational lenses or dual AGN candidates. Detailed follow-up radio interferometric and optical spectroscopic observations are encouraged to investigate the outlier sources found in this study.

Context. Blazars are powerful active galactic nuclei (AGNs) radiating prominently in the whole electromagnetic spectrum, from the radio to the X-ray and gamma-ray bands. Their emission is dominated by synchrotron and inverse-Compton radiation from a relativistic jet originating from an accreting central supermassive black hole. The object IGR J12319-0749 has recently been identified as a soft gamma-ray source with the IBIS instrument of the INTEGRAL satellite, coincident with a quasar at high redshift (z=3.12). Aims. We studied the radio structure of IGR J12319-0749 to strengthen its blazar identification by detecting a compact radio jet on the milli-arcsecond (mas) angular scale, and to measure its astrometric position accurate to mas level. Methods. We used the technique of electronic very long baseline interferometry (e-VLBI) to image IGR J12319-0749 with the European VLBI Network (EVN) at 5 GHz on 2012 June 19. Results. IGR J12319-0749 (J1231-0747) is a compact radio source, practically unresolved on interferometric baselines up to ~136 million wavelengths. The estimated brightness temperature (at least ~2 x 10^12 K) indicates that the radio emission of its jet is Doppler-boosted. The accurate position of the compact radio source is consistent with the positions measured at higher energies.

Based on its broad-band spectral energy distribution, and the X-ray spectrum in particular, the radio-loud active galactic nucleus (AGN) SDSS J102623.61+254259.5 (J1026+2542) has recently been classified as a blazar. The extremely high redshift of the source, z=5.3, makes it one of the most distant and most luminous radio-loud AGN known to date. From published 5-GHz very long baseline interferometry (VLBI) imaging data obtained in 2006, the source has a typical blazar appearance on mas scales, with a prominent one-sided jet extending to ~20 mas. We estimate the brightness temperature of J1026+2542 and find no strong evidence for Doppler boosting. The jet viewing angle is possibly at least ~20 deg. The bulk Lorentz factor and the viewing angle of the jet could reliably be determined in the near future from multi-epoch VLBI observations.

HESS J1943+213, a TeV point source close to the Galactic plane recently discovered by the H.E.S.S. collaboration, was proposed to be an extreme BL Lacertae object, though a pulsar wind nebula (PWN) nature could not be completely discarded. To investigate its nature, we performed high-resolution radio observations with the European Very Long Baseline Interferometry Network (EVN) and reanalyzed archival continuum and H {\sc i} data. The EVN observations revealed a compact radio counterpart of the TeV source. The low brightness temperature and the resolved nature of the radio source are indications against the beamed BL Lacertae hypothesis. The radio/X-ray source appears immersed in a $\sim$ 1\arcmin elliptical feature suggesting a possible galactic origin (PWN nature) for the HESS source. We found that HESS\,J1943+213 is located in the interior of a $\sim1\degr$ diameter H {\sc i} feature, and explored the possibility of they being physically related.

Aims. A long timeline kinematic study of the archetypal CSO OQ 208 sheds light on the physical properties of the most compact radio sources. Methods. Archival data from the VLBA at 15 GHz over a time span of 13.6 yr are used to investigate the kinematics of the radio source. The flux density monitoring data obtained at the Michigan 26-meter radio telescope are also used as supplementary information. Results. At 8.4 and 15 GHz, the two lobes are resolved into two sub-components, identified as hotspots. A knotty jet is linked with the NE hotspot and traces back toward the geometric center. The core is too weak to be detected. Significant flux density variation is found in the primary hotspots with the maximum level of 62% (NE1) and 19% (SW1). The peak in the flux density of NE1 leads that of SW1 by approximately 5.00 yr, suggesting that the northeast lobe is advancing and the southwest lobe is receding. This light travel difference indicates a radial distance difference between the two hotspots of 1.53 pc, which indicates an inclination angle of about 80.8 degree between the radio jet and the line of sight. The angular separation rate between NE1 and SW1 is 0.027 mas/yr (or 0.133 c). The inner jet knot moves at 0.047 mas/yr (or 0.230 c), about 3.5 times the hotspot advancing speed. Conclusions. The large viewing angle and the modest jet speed suggest a mildly relativistic jet. The jet axis is close to the plane of the sky. The separation rate and the distance between the two primary hotspots result in a kinematic age of 255$\pm$17 yr, confirming that OQ 208 is indeed a young radio source. In addition to the hotspot advancing motions, sideways motions provide evidence that the lobes are obstructed by the external interstellar medium.

The Swift/XRT light curve of the ultraluminous X-ray (ULX) source NGC 5408 X-1 was re-analyzed with two new numerical approaches, Weighted Wavelet $Z$-transform (WWZ) and CLEANest, that are different from previous studies. Both techniques detected a prominent periodicity with a time scale of $115.5\pm1.5$ days, in excellent agreement with the detection of the same periodicity first reported by Strohmayer (2009). Monte Carlo simulation was employed to test the statisiticak confidence of the 115.5-day periodicity, yielding a statistical significance of $> 99.98%$ (or $>3.8\sigma$). The robust detection of the 115.5-day quasi-periodic oscillations (QPOs), if it is due to the orbital motion of the binary, would infer a mass of a few thousand $M_\odot$ for the central black hole, implying an intermediate-mass black hole in NGC 5408 X-1.

The highest-redshift quasars are still rare and valuable objects for observational astrophysics and cosmology. They provide important constraints on the growth of the earliest supermassive black holes in the Universe, and information on the physical conditions in their environment. Among the nearly 60 quasars currently known at redshifts z>5.7, only a handful are "strong" emitters in radio continuum. These can be targets of sensitive high-resolution Very Long Baseline Interferometry (VLBI) observations to reveal their innermost structure, down to ~10 pc linear scales. We review the results of our earlier European VLBI Network (EVN) experiments on three of the most distant radio quasars known to date, and give a preliminary report on the EVN detection of a fourth one. The results obtained so far suggest that we see really young active galactic nuclei - not just in a cosmological sense but also in terms of their active life in radio.

The radio-emitting quasar SDSS J1425+3231 (z=0.478) was recently found to have double-peaked narrow [O III] optical emission lines. Based on the analysis of the optical spectrum, Peng et al. (2011) suggested that this object harbours a dual active galactic nucleus (AGN) system, with two supermassive black holes (SMBHs) separated on the kpc scale. SMBH pairs should be ubiquitous according to hierarchical galaxy formation scenarios in which the host galaxies and their central black holes grow together via interactions and eventual mergers. Yet the number of presently-confirmed dual SMBHs on kpc or smaller scales remains small. A possible way to obtain direct observational evidence for duality is to conduct high-resolution radio interferometric measurements, provided that both AGN are in an evolutionary phase when some activity is going on in the radio. We used the technique of Very Long Baseline Interferometry (VLBI) to image SDSS J1425+3231. Observations made with the European VLBI Network (EVN) at 1.7 GHz and 5 GHz frequencies in 2011 revealed compact radio emission at sub-mJy flux density levels from two components with a projected linear separation of \sim2.6 kpc. These two components support the possibility of a dual AGN system. The weaker component remained undetected at 5 GHz, due to its steep radio spectrum. Further study will be necessary to securely rule out a jet--shock interpretation of the less dominant compact radio source. Assuming the dual AGN interpretation, we discuss black hole masses, luminosities, and accretion rates of the two components, using available X-ray, optical, and radio data. While high-resolution radio interferometric imaging is not an efficient technique to search blindly for dual AGN, it is an invaluable tool to confirm the existence of selected candidates.

With Gaia, it will become possible to directly link the radio and optical reference frames using a large number of common objects. For the most accurate radio-optical link, it is important to know the level of spatial coincidence between the quasars' optical positions, and the radio positions determined by Very Long Baseline Interferometry (VLBI) observations. The "outlier" objects, for which the positions are significantly offset at the two different electromagnetic wavebands, may be of astrophysical interest as well. Here we present a case study to compare the radio positions of ~800 quasars common in the second realization of the International Celestial Reference Frame (ICRF2) and in the Sloan Digital Sky Survey Data Release 7 (SDSS DR7) catalogue. Compared to the radio ICRF2, the SDSS provides two orders of magnitude less accurate astrometric data in the optical. However, its extensive sky coverage and faint magnitude limit allow us to directly relate the positions of a large sample of radio sources. This way we provide an independent check of the overall accuracy of the SDSS positions and confirm that the astrometric calibration of the latest Data Release 8 (DR8) is poorer than that of the DR7. We find over 20 sources for which the optical and radio brightness peaks are apparently not coincident at least at the 3-sigma level of SDSS DR7 positional accuracy, and briefly discuss the possible causes, including dual active galactic nuclei.

The point source list of the Wilkinson Microwave Anisotropy Probe (WMAP) is a uniform, all-sky catalogue of bright sources with flux density measurements at high (up to 94 GHz) radio frequencies. We investigated the five-year WMAP list to compile a new catalogue of bright and compact extragalactic radio sources to be potentially studied with Very Long Baseline Interferometry at millimeter wavelengths (mm-VLBI) and Space VLBI (SVLBI). After comparing the WMAP data with the existing mm-VLBI catalogues, we sorted out the yet unexplored sources. Using the 41, 61 and 94 GHz WMAP flux densities, we calculated the spectral indices. By collecting optical identifications, lower-frequency radio flux densities and VLBI images from the literature, we created a list of objects which have not been investigated with VLBI at 86 GHz before. With total flux density at least 1 Jy and declination above -40 degree, we found 37 suitable new targets. It is a nearly 25% addition to the known mm-VLBI sources. Such objects are also potentially useful as phase-reference calibrators for the future Japanese SVLBI mission ASTRO-G at its highest observing frequency (43 GHz). The phase-referencing capability of ASTRO-G would allow long integrations and hence better sensitivity for observing faint target sources close to suitable phase calibrators in the sky.

Context: There are about 60 quasars known at redshifts z>5.7 to date. Only three of them are detected in the radio above 1 mJy flux density at 1.4 GHz frequency. Among them, J1429+5447 (z=6.21) is the highest-redshift radio quasar known at present. These rare, distant, and powerful objects provide important insight into the activity of the supermassive black holes in the Universe at early cosmological epochs, and on the physical conditions in their environment. Aims: We studied the compact radio structure of J1429+5447 on the milli-arcsecond (mas) angular scale, in order to compare the structural and spectral properties with those of other two z~6 radio-loud quasars, J0836+0054 (z=5.77) and J1427+3312 (z=6.12). Methods: We performed Very Long Baseline Interferometry (VLBI) imaging observations of J1429+5447 with the European VLBI Network (EVN) at 1.6 GHz on 2010 June 8, and at 5 GHz on 2010 May 27. Results: Based on its observed radio properties, the compact but somewhat resolved structure on linear scales of <100 pc, and the steep spectrum, the quasar J1429+5447 is remarkably similar to J0836+0054 and J1427+3312. To answer the question whether the compact steep-spectrum radio emission is a "universal" feature of the most distant radio quasars, it is essential to study more, yet to be discovered radio-loud active galactic nuclei at z>6.

There are about 50 quasars known at redshifts z>5.7 to date. Only three of them are detected in the radio (J0836+0054, z=5.77; J1427+3312, z=6.12; J1429+5447, z=6.21). The highest-redshift quasars are in the forefront of current astrophysical and cosmological research since they provide important constraints on the growth of the earliest supermassive black holes in the Universe, and on the physical conditions in their environment. These sources are indeed associated with active galactic nuclei as revealed by high-resolution Very Long Baseline Interferometry (VLBI) observations. It is still unclear whether the physical properties of the few z~6 radio quasars are in general similar to those of their lower-redshift cousins. In the case of J1427+3312, the 100-pc scale double morphology suggests a young radio source. Here we report on the recent European VLBI Network (EVN) imaging observations of J1429+5447, the most distant radio quasar known. Based on its milli-arcsecond-scale structural and spectral properties, this quasar is similar to J0836+0054 and J1427+3312. This raises the question if the compact steep-spectrum radio emission is a "universal" feature of the most distant radio quasars. It could only be answered after many more objects at z>6 are discovered and studied.

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.

Recently Jamrozy et al. (2009) identified 4C 02.27 (J0935+0204) as the first possible example of a double-double radio source which is optically identified with a quasar (i.e. not a galaxy), at the redshift of z=0.649. The overall projected angular size of the radio source reaches about 1.5', with a prominent "core" component in the centre. The two opposite pairs of radio lobes might indicate two periods of episodic activity. We report on our short exploratory 1.6-GHz Very Long Baseline Interferometry (VLBI) observations of the innermost radio structure of the quasar, conducted with the electronic European VLBI Network (e-EVN) on 2009 September 30. These revealed a milliarcsecond-scale compact source which is the base of the approaching one of the two symmetric relativistic jets currently supplying the hot spots in the inner pair of the arcsecond-scale radio lobes in 4C 02.27.

Context: Until now, there have only been seven quasars at z>4.5 whose the high-resolution radio structure had been studied in detail with Very Long Baseline Interferometry (VLBI) imaging. Aims: We almost double the number of VLBI-imaged quasars at these high redshifts with the aim of studying their redshift-dependent structural and physical properties in a larger sample. Methods: We observed five radio quasars (J0813+3508, J1146+4037, J1242+5422, J1611+0844, and J1659+2101) at 4.5<z<5 with the European VLBI Network (EVN) at 1.6 GHz on 29 October 2008 and at 5 GHz on 22 October 2008. The angular resolution achieved ranges from 1.5 to 25 milli-arcseconds (mas), depending on the observing frequency, the position angle in the sky, and the source's celestial position. Results: The sources are all somewhat extended on mas scales, but compact enough to be detected at both frequencies. With one exception of a flat-spectrum source (J1611+0844), their compact emission is characterised by a steep radio spectrum. We found no evidence of Doppler-boosted radio emission in the quasars in our sample. The radio structure of one of them (J0813+3508) is extended to ~7", which corresponds to 43 kpc projected linear size. Many of the highest redshift compact radio sources are likely to be young, evolving objects, far-away cousins of the powerful gigahertz peaked-spectrum (GPS) and compact steep-spectrum (CSS) sources that populate the Universe at lower redshifts.

Aims. We selected two radio quasars (J1036+1326 and J1353+5725) based on their 1.4-GHz radio structure, which is dominated by a bright central core and a pair of weaker and nearly symmetric lobes at ~10" angular separation. They are optically identified in the Sloan Digital Sky Survey (SDSS) at spectroscopic redshifts z>3. We investigate the possibility that their core-dominated triple morphology can be a sign of restarted radio activity in these quasars, involving a significant repositioning of the radio jet axis. Methods. We present the results of high-resolution radio imaging observations of J1036+1326 and J1353+5725, performed with the European Very Long Baseline Interferometry (VLBI) Network (EVN) at 1.6 GHz. These data are supplemented by archive observations from the Very Large Array (VLA).We study the large- and small-scale radio structures and the brightness temperatures, then estimate relativistic beaming parameters. Results. We show that the central emission region of these two high-redshift, core-dominated triple sources is compact but resolved at ~10 milli-arcsecond resolution. We find that it is not necessary to invoke large misalignment between the VLBI jet and the large-scale radio structure to explain the observed properties of the sources.

Context. The high-redshift (z = 4.72) blazar J1430+4204 produced a major radio outburst in 2005. Such outbursts are usually associated with the emergence of a new component in the inner radio jet. Aims. We searched for possible changes in the radio structure on milli-arcsecond angular scales, to determine physical parameters that characterise the relativistic jet ejected from the centre of this source. Methods. We analysed 15-GHz radio interferometric images obtained with the Very Long Baseline Array (VLBA) before and after the peak of the outburst. Results. We did not identify any significant new jet component over a period of 569 days. We estimated the Doppler factor, the Lorentz factor, and the apparent transverse speed of a putative jet component using three different methods. The likely small jet angle to the line of sight and our values of the apparent transverse speed are consistent with not detecting a new jet feature.

The high-redshift (z=4.72) blazar J1430+4204 produced an exceptional radio outburst in 2006. We analyzed 15-GHz radio interferometric images obtained with the Very Long Baseline Array (VLBA) before and after the outburst, to search for possible structural changes on milli-arcsecond angular scales and to determine physical parameters of the source.

Context: There is a possibility that the optically unidentified radio source J1218+2953 may act as a gravitational lens, producing an optical arc ~4" away from the radio position. Until now, the nature of the lensing object has been uncertain since it is not detected in any waveband other than the radio. The estimated high mass-to-light ratio could even allow the total mass of this galaxy to be primarily in the form of dark matter. In this case, J1218+2953 could be the first known example of a "dark lens". Aims: We investigate the nature of J1218+2953 by means of high-resolution radio imaging observations to determine whether there is a radio-loud active galactic nucleus (AGN) in the position of the lensing object. Methods: We report on Very Long Baseline Interferometry (VLBI) observations with the European VLBI Network (EVN) at 1.6 and 5 GHz. Results: Our images, having angular resolutions of ~1 to ~10 milli-arcseconds (mas), reveal a rich and complex radio structure extending to almost 1". Based on its radio spectrum and structure, J1218+2953 can be classified as a compact steep-spectrum (CSS) source, and as a medium-size symmetric object (MSO). The object harbours an AGN. It is also found as an X-ray source in the XMM-Newton EPIC (European Photon Imaging Cameras) instrument serendipitous source catalogue. Conclusions: Rather than being a dark lens, J1218+2953 is most likely a massive, heavily obscured galaxy in which the nuclear activity is currently in an early evolutionary stage.

The e-VLBI technique offers a unique opportunity for users to probe the milliarcsecond (mas) scale structure of unidentified radio sources, and organise quick follow-up observations in case of detection. Here we report on e-EVN results for a peculiar radio source that has been suggested to act as a gravitational lens. However the lensing galaxy has not been identified in the optical or the IR bands so far. Our goal was to look for an active galactic nucleus (AGN) in this suspected dark lens system. The results indicate strong AGN activity, and rule out the possibility that the radio source itself is gravitationally lensed.

The upper envelope of the amplitude of the VLBI visibility function usually represents the most compact structural pattern of extragalactic radio sources, in particular, the "core-jet" morphologies. By fitting the envelope to a circular Gaussian model in ~3000 parsec-scale core-jet structures, we find that the apparent angular size shows significant power-law dependence on the observing frequency (power index n = -0.95 pm 0.37). The dependence is likely to result from synchrotron self-absorption in the inhomogeneous jet and not the free-free absorption (n = -2.5), nor the simple scatter broadening (n leq -2).

We show an archival 1.4-GHz Very Large Array (VLA) image of the core-dominated triple (CDT) radio source J1628+3906, along with a recent high-resolution 5-GHz Very Long Baseline Interferometry (VLBI) image of its core. This core is resolved into two nearly symmetric components at ~10 milli-arcsecond (mas) angular scale. The inner structure is markedly misaligned with respect to the large-scale structure. This morphology can be interpreted as a possible signature of once ceased and then restarted activity of the galactic nucleus, coupled with a rapid repositioning of the central radio jet axis. In a search for compact (<5"), optically identified, and highly redshifted (z>3) quasars in the Faint Images of the Radio Sky at Twenty-centimeters (FIRST) and the Sloan Digital Sky Survey (SDSS) catalogues, we have also found two quasars associated with weaker two-sided radio lobes at ~10" scale. We present archival 1.4-GHz VLA images of these quasars (J1036+1326 and J1353+5725). The compact cores of both of these CDT sources are likely to be detectable at mas scales. Future VLBI observations at the highest resolution could reveal if their inner radio structure is similar to that of the Compact Symmetric Objects (CSOs). The characteristic directions of their arcsecond-scale and mas-scale structures are probably misaligned which could also support the "restarting and rapid repositioning" scenario.

We apply an efficient selection method to identify potential weak Very Long Baseline Interferometry (VLBI) target quasars simply using optical (SDSS) and low-resolution radio (FIRST) catalogue data. Our search is restricted to within 12" from known compact radio sources that are detectable as phase-reference calibrators for ASTRO-G at 8.4 GHz frequency. These calibrators have estimated correlated flux density >20 mJy on the longest ground-space VLBI baselines. The search radius corresponds to the primary beam size of the ASTRO-G antenna. We show that ~20 quasars with at least mJy-level expected flux density can be pre-selected as potential in-beam phase-reference targets for ASTRO-G at 8.4 GHz frequency. Most of them have never been imaged with VLBI. The sample of these dominantly weak sources offers a good opportunity to study their radio structures with unprecedented angular resolution provided by Space VLBI. The method of in-beam phase-referencing is independent from the ability of the orbiting radio telescope to do rapid position-switching manoeuvres between the calibrators and the nearby reference sources, and less sensitive to the satellite orbit determination uncertainties.

We investigate the frequency-dependent radio properties of the jet of the luminous high-redshift (z = 3.2) radio quasar PKS 1402+044 (J1405+0415) by means of radio interferometric observations. The observational data were obtained with the VLBI Space Observatory Programme (VSOP) at 1.6 and 5 GHz, supplemented by other multi-frequency observations with the Very Long Baseline Array (VLBA; 2.3, 8.4, and 15 GHz) and the Very Large Array (VLA; 1.4, 5, 15, and 43 GHz). The observations span a period of 7 years. We find that the luminous high-redshift quasar PKS 1402+044 has a pronounced "core-jet" morphology from the parsec to the kilo-parsec scales. The jet shows a steeper spectral index and lower brightness temperature with increasing distance from the jet core. The variation of brightness temperature agrees well with the shock-in-jet model. Assuming that the jet is collimated by the ambient magnetic field, we estimate the mass of the central object as ~10^9 M_sun. The upper limit of the jet proper motion of PKS 1402+044 is 0.03 mas/yr (~3c) in the east-west direction.

Context: The highest redshift quasars at z>~6 receive considerable attention since they provide strong constraints on the growth of the earliest supermassive black holes. They also probe the epoch of reionisation and serve as "lighthouses" to illuminate the space between them and the observer. The source J1427+3312 (z=6.12) has recently been identified as the first and so far the only known radio-loud quasar at z>6. Aims: We investigated the compact radio structure of J1427+3312 on milli-arcsecond (mas) angular scales, to compare it with that of the second most distant radio-loud quasar J0836+0054 (z=5.77) and with lower-redshift radio quasars in general. Methods: We observed J1427+3312 in phase-reference mode with ten antennas of the European Very Long Baseline Interferometry (VLBI) Network (EVN) at 1.6 GHz on 11 March 2007 and at 5 GHz on 3 March 2007. Results: The source was clearly detected at both frequencies. At 1.6 GHz, it shows a prominent double structure. The two components are separated by 8.3 mas, corresponding to a projected linear distance of ~160 pc. Both components with sub-mJy flux densities appear resolved. In the position of the brightest component at 1.6 GHz, we detected mas-scale radio emission at 5 GHz as well. The radio spectrum of this feature is steep. The double structure and the separation of the components of J1427+3312 are similar to those of the young (<~10^4 yr) compact symmetric objects (CSOs).

We studied the radio structure of high-redshift (z>3) quasars with VSOP at 1.6 and 5 GHz. These sources are the most distant objects ever observed with Space VLBI, at rest-frame frequencies up to ~25 GHz. Here we give an account of the observations and briefly highlight the most interesting cases and results. These observations allowed us, among other things, to estimate the mass of the central black holes powering these quasars, to identify large misalignments between the milli-arcsecond (mas) and sub-mas scale radio structures, and to detect apparent superluminal motion at sub-mas scale.

We show that as many as ~50 quasars with at least mJy-level expected flux density can be pre-selected as potential in-beam phase-reference targets for ASTRO-G. Most of them have never been imaged with VLBI. These sources are located around strong, compact calibrator sources that have correlated flux density >100 mJy on the longest VLBA baselines at 8.4 GHz. All the targets lie within 12' from the respective reference source. The basis of this selection is an efficient method to identify potential weak VLBI target quasars simply using optical and low-resolution radio catalogue data. The sample of these dominantly weak sources offers a good opportunity for a statistical study of their radio structure with unprecedented angular resolution at 8.4 GHz.

Context: The Deep Extragalactic VLBI-Optical Survey (DEVOS) aims at constructing a large sample of compact radio sources up to two orders of magnitude fainter than those studied in other Very Long Baseline Interferometry (VLBI) surveys. Optical identification of the objects is ensured by selecting them from the Sloan Digital Sky Survey (SDSS) list. Aims: While continuing to build up the DEVOS data base, we investigated how the VLBI detection rate could be enhanced by refining the initial selection criteria introduced in the first paper of this series. Methods: We observed 26 sources in two adjacent, slightly overlapping 2 deg radius fields with the European VLBI Network (EVN) at 5 GHz frequency on 2 March 2007.The phase-reference calibrator quasars were J1616+3621 and J1623+3909. The objects selected were unresolved both in the Faint Images of the Radio Sky at Twenty-centimeters (FIRST) survey catalogue and the SDSS Data Release 4. Results: We present images of milli-arcsecond (mas) scale radio structures and accurate coordinates of 24 extragalactic sources. Most of them have never been imaged with VLBI. Twenty-two compact radio sources (85% of our initial sample) are considered as VLBI detections of the corresponding optical quasars in SDSS. We found an efficient way to identify quasars as potential VLBI targets with mas-scale compact radio stucture at >1 mJy level, based only on the FIRST and SDSS catalogue data by applying simple selection criteria.

In February 1997, the Japanese radio astronomy satellite HALCA was launched to provide the space-bourne element for the VLBI Space Observatory Programme (VSOP) mission. Approximately twenty-five percent of the mission time was dedicated to the VSOP Survey of bright compact Active Galactic Nuclei (AGN) at 5 GHz. This paper, the fifth in the series, presents images and models for the remaining 140 sources not included in Paper III, which contained 102 sources. For most sources, the plots of the uv-coverage, the visibility amplitude versus uv-distance, and the high resolution image are presented. Model fit parameters to the major radio components are determined, and the brightness temperature of the core component for each source is calculated. The brightness temperature distributions for all of the sources in the VSOP AGN survey are discussed.

In February 1997 the Japanese radio astronomy satellite HALCA was launched to provide the space-borne element for the VSOP mission. HALCA provided linear baselines three-times greater than that of ground arrays, thus providing higher resolution and higher AGN brightness temperature measurements and limits. Twenty-five percent of the scientific time of the mission was devoted to the ``VSOP survey'' of bright, compact, extra-galactic radio sources at 5 GHz. A complete list of 294 survey targets were selected from pre-launch surveys, 91% of which were observed during the satellite's lifetime. The major goals of the VSOP Survey are statistical in nature: to determine the brightness temperature and approximate structure, to provide a source list for use with future space VLBI missions, and to compare radio properties with other data throughout the electromagnetic spectrum. All the data collected have now been analysed and is being prepared for the final image Survey paper. In these papers we present details of the mission, images of the sources, and some statistics of the contributions and completeness plus some statistics of the images and brightness temperatures. Full resolution images are to be found on http://www.vsop.isas.jaxa.jp/survey/publications/evn.ps.gz

Systematic surveys of astronomical objects often lead to discoveries, but always provide invaluable information for statistical studies of well-defined samples. They also promote follow-up investigations of individual objects or classes. Surveys using a yet unexplored observing wavelength, a novel technique or a new instrument are of special importance. Significantly improved observing parameters (e.g. sensitivity, angular resolution, monitoring capability) provide new insight into the morphological and physical properties of the objects studied. I give a brief overview of the important Very Long Baseline Interferometry (VLBI) imaging surveys conducted in the past. A list of surveys guides us through the developments up until the present days. I also attempt to show directions for the near future.

Based on the VLBI Space Observatory Programme (VSOP) observations at 1.6 and 5 GHz, we find that the luminous high-redshift (z=3.215) quasar PKS 1402+044 (J1405+0415) has a pronounced 'core--jet' structure. The jet shows a steeper spectral index and lower brightness temperature with the increase of the distance from the core. The variation of brightness temperature is basically consistent with the shock-in-jet model. Assuming that the jet is collimated by the ambient magnetic field, we estimate the mass of the central object as ~10^9 M_sun. The upper limit of the jet proper motion of PKS 1402+044 is 0.03 mas/yr (~3c) in the east-west direction.

We discuss the importance of the direct link between the most accurate radio and optical reference frames that will become possible with the next-generation space astrometry missions in about a decade. The positions of more than 500 active galactic nuclei that are common in the precise Very Long Baseline Interferometry (VLBI) catalogues and the Sloan Digital Sky Survey (SDSS) Data Release 4 (DR4) are compared. While obtaining an ``independent'' estimate for the SDSS coordinate accuracies, we find indications that the assumption of spatially coincident brightness peaks for the same objects in radio and optical does not hold for each object.

In 1999, the highly compact and variable BL Lac object AO 0235+164 was identified as the highest brightness temperature active galactic nucleus observed with the VLBI Space Observatory Programme (VSOP), with T_B > 5.8 x 10^{13} K. The sub-milliarcsecond radio structure of this source has been studied with dual-frequency (1.6 and 5 GHz), polarization-sensitive VSOP observations during 2001 and 2002. Here we present the results of this monitoring campaign. At the time of these observations, the source was weakly polarized and characterized by a radio core that is clearly resolved on space-ground baselines.

We present the results of the pilot observations of the Deep Extragalactic VLBI-Optical Survey (DEVOS). Our ultimate aim is to collect information on compact structures in a large sample of extragalactic radio sources (~10000 objects) up to two orders of magnitude fainter than those studied in typical imaging Very Long Baseline Interferometry (VLBI) surveys up until now. This would lead to an unprecedented data base for various astrophysical, astrometric and cosmological studies. The first global VLBI observations of the DEVOS programme were successfully conducted in May 2002. We selected sources without any spectral criterion from the Very Large Array (VLA) Faint Images of the Radio Sky at Twenty-centimeters (FIRST) catalogue, that are also detected with the Multi-Element Radio Linked Interferometer Network (MERLIN). The DEVOS pilot sample sources are in the area of the sky that is covered by the Sloan Digital Sky Survey (SDSS). We describe the sample selection and present high resolution 5-GHz radio images of the sources. Based on the results of this pilot study, we estimate the outcome of and the resources needed for a full-scale DEVOS project.

We observed SDSS J0836+0054, the most distant radio-loud quasar known at present (z=5.774) with the European Very Long Baseline Interferometry (VLBI) Network at the 5 GHz frequency on 2003 November 4. The source is compact at the milli-arcsecond (mas) angular scale with a flux density of 0.34 mJy. The observations with the Very Large Array (VLA) made on the consecutive day support a conclusion that the radio emission from SDSS J0836+0054 is essentially confined within the central 40 pc. Based on our phase-referenced VLBI observation, we obtained the astrometric position of the source accurate to ~2 mas.

We report on European VLBI Network (EVN) and Very Large Array (VLA) observations of SDSS J0836+0054, the most distant radio-loud quasar known at present z=5.774). The source, with a total flux density of 1.1 mJy, shows a compact structure at 10 mas angular resolution at 1.6 GHz frequency, with no indication of multiple images produced by gravitational lensing within a 4x4 arcsec field above the brightness level of 100 microJy/beam. More recent 5 GHz EVN observations are being analysed now. As a by-product, a two-sided mas-scale jet structure is suspected in one of the phase-reference calibrator quasars, PKS 0837+012.

The VSOP mission is a Japanese-led project to study radio sources with sub-milliarcsecond angular resolution using an orbiting 8-m telescope, HALCA and global arrays of Earth-based telescopes. Approximately 25% of the observing time has been devoted to a survey of compact AGN at 5 GHz which are stronger than 1 Jy -- the VSOP AGN Survey. This paper, the second in a series, describes the data calibration, source detection, self-calibration, imaging and modeling, and gives examples illustrating the problems specific to space VLBI. The VSOP Survey web-site which contains all results and calibrated data is described.

The VSOP mission is a Japanese-led project to study radio sources with sub-milliarcsec resolution using an orbiting 8 m telescope, HALCA, along with global arrays of Earth-based telescopes. Approximately 25% of the observing time is devoted to a survey of compact AGN which are stronger than 1 Jy at 5 GHz-the VSOP AGN Survey. This paper, the third in the series, presents the results from the analysis of the first 102 Survey sources. We present high resolution images and plots of visibility amplitude versus projected baseline length. In addition, model-fit parameters to the primary radio components are listed, and from these the angular size and brightness temperature for the radio cores are calculated. For those sources for which we were able to determine the source frame core brightness temperature, a significant fraction (53 out of 98) have a source frame core brightness temperature in excess of 10^12 K. The maximum source frame core brightness temperature we observed was 1.2 X 10^13 K. Explaining a brightness temperature this high requires an extreme amount of relativistic Doppler beaming. Since the maximum brightness temperature one is able to determine using only ground-based arrays is of the order of 10^12 K, our results confirm the necessity of using space VLBI to explore the extremely high brightness temperature regime.

We present a high resolution radio image of SDSS 0836+0054 identified recently as the most distant radio-detected quasar at a redshift of z=5.82. The observation was carried out with ten antennas of the European VLBI Network, spread from Europe to China and South Africa, at 1.6 GHz frequency on 2002 June 8. The source is detected with a total flux density of 1.1 mJy, equal to its flux density measured in the VLA FIRST survey. We found no indication of multiple images produced by gravitational lensing. The radio structure of the quasar at ~10-mas angular resolution appears somewhat resolved. It resembles the radio structure typical for lower redshift radio-loud active galactic nuclei. We obtained so far the best astrometric position of the source with an accuracy better than 8 mas, limited mainly by the structural effects in the phase-reference calibrator source.

Space very long baseline interferometry (VLBI) observations of the quasar 3C216 with the VLBI Space Observatory Programme (VSOP) reveal that the parsec scale structure of the source can be well described by compact jet models developed for interpreting the core region of radio-loud active galactic nuclei (AGN). The measured brightness temperature of T_(b) = 7.9*10^{11} K is comparable to the inverse Compton limit, from which we determine a lower limit of delta ~ 3.17 for the Doppler-boosting factor. The apparent transvers velocity of the superluminal component is beta_(app) = (3.0+-0.2) h^{-1} asuming a constant velocity, but deceleration of the jet material cannot be excluded from our data. A combination of the above values indicates that the viewing angle of the core-jet to the line-of-sight is less than 18.4 degrees, and the jet Lorentz-factor exceeds 3.16. The observed small size of the source is probably caused by both interaction with the interstellar medium, and a projection effect.

We report results of a VSOP (VLBI Space Observatory Programme) observation of a high redshift quasar PKS 2215+020 (z=3.572). The ~1 milliarcsecond resolution image of the quasar reveals a prominent `core-jet' structure on linear scales from 5/h to 300/h pc ($H_0=100*h km/(s*Mpc). The brightness temperatures and sizes of bright features identified in the jet are consistent with emission from relativistic shocks dominated by adiabatic energy losses. The jet is powered by the central black hole with estimated mass of ~4*10^9 solar masses. Comparisons with VLA and ROSAT observations indicate a possible presence of an extended radio/X-ray halo surrounding 2215+020.

AO 0235+164 is a very compact, flat spectrum radio source identified as a BL Lac object at a redshift of z=0.94. It is one of the most violently variable extragalactic objects at both optical and radio wavelengths. The radio structure of the source revealed by various ground-based VLBI observations is dominated by a nearly unresolved compact component at almost all available frequencies. Dual-frequency space VLBI observations of AO 0235+164 were made with the VSOP mission in January-February 1999. The array of the Japanese HALCA satellite and co-observing ground radio telescopes in Australia, Japan, China and South Africa allowed us to study AO 0235+164 with an unprecedented angular resolution at frequencies of 1.6 and 5 GHz. We report on the sub-milliarcsecond structural properties of the source. The 5-GHz observations led to an estimate of T_B > 5.8 x 10^{13} K for the rest-frame brightness temperature of the core, which is the highest value measured with VSOP to date.

We present very long baseline interferometry (VLBI) images of ten very high redshift (z>3) quasars at 5 GHz. The sources 0004+139, 0830+101, 0906+041, 0938+119 and 1500+045 were observed in September 1992 using a global VLBI array, while 0046+063, 0243+181, 1338+381, 1428+423 and 1557+032 were observed in October 1996 with the European VLBI Network and Hartebeesthoek, South Africa. Most of the sources are resolved and show asymmetric structure. The sample includes 1428+423, the most distant radio loud quasar known to date (z=4.72). It is barely resolved with an angular resolution of about 2.0*1.4 mas.

We discuss the ``angular size - redshift'' relation for compact radio sources distributed over a wide range of redshifts 0.01 < z < 4.73. Our study is based on a sample of 330 5 GHz VLBI contour maps taken from the literature. Unlike extended source samples, the ``angular size - redshift'' relation for compact radio sources appears consistent with the predictions of standard Friedmann world models with q_o ~ 0.5 without the need to consider evolutionary or selection effects due to a ``linear size -- luminosity'' dependence. By confining our analysis to sources having a spectral index, -0.38 < alpha < 0.18, and a total radio luminosity greater than 10^{26} W/Hz, we are able to restrict the dispersion in the ``angular size -- redshift'' relation. The best fitting regression analysis in the framework of the Friedmann-Robertson-Walker model gives the value of the deceleration parameter q_o = 0.21 +- 0.30 if there are no evolutionary or selection effects due to a ``linear size -- luminosity'', ``linear size -- redshift'' or ``linear size -- spectral index'' dependence.