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.

We present our results from a mid-infrared interferometric survey targeted at the planet-forming region in the circumstellar disks around low- and intermediate-mass young stars. Our sample consist of 82 objects, including T Tauri stars, Herbig Ae stars, and young eruptive stars. Our main results are: 1) Disks around T Tauri stars are similar to those around Herbig Ae stars, but are relatively more extended once we account for stellar luminosity. 2) From the distribution of the sizes of the mid-infrared emitting region we find that inner dusty disk holes may be present in roughly half of the sample. 3) Our analysis of the silicate spectral feature reveals that the dust in the inner $\sim$1 au region of disks is generally more processed than that in the outer regions. 4) The dust in the disks of T Tauri stars typically show weaker silicate emission in the N band spectrum, compared to Herbig Ae stars, which may indicate a general difference in the disk structure. Our data products are available at VizieR, and at the following web page: http://konkoly.hu/MIDI_atlas

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.

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.

We present new Karl G. Jansky Very Large Array (VLA, 1.5 GHz) radio data for the giant elliptical galaxy IC 4296, supported by archival radio, X-ray (Chandra, XMM-Newton) and optical (SOAR, HST) observations. The galaxy hosts powerful radio jets piercing through the inner hot X-ray emitting atmosphere, depositing most of the energy into the ambient intra-cluster medium (ICM). Whereas the radio surface brightness of the A configuration image is consistent with a Fanaroff-Riley Class I (FR I) system, the D configuration image shows two bright, relative to the central region, large (~160 kpc diameter), well-defined lobes, previously reported by Killeen et al., at a projected distance r~>230 kpc. The XMM-Newton image reveals an X-ray cavity associated with one of the radio lobes. The total enthalpy of the radio lobes is ~7x10^59 erg and the mechanical power output of the jets is ~10^44 erg/s. The jets are mildly curved and possibly re-brightened by the relative motion of the galaxy and the ICM. The lobes display sharp edges, suggesting the presence of bow shocks, which would indicate that they are expanding supersonically. The central entropy and cooling time of the X-ray gas are unusually low and the nucleus hosts a warm H\alpha+[NII] nebula and a cold molecular CO disk. Because most of the energy of the jets is deposited far from the nucleus, the atmosphere of the galaxy continues to cool, apparently feeding the central supermassive black hole and powering the jet activity.

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.

We studied the disc of the unclassified B[e] star HD 50138, in order to explore its structure, and to find indications for the evolutionary status of this system, whether it is a young Herbig Be or a post-main-sequence star. Using high spatial resolution interferometric measurements from MIDI instrument (N-band) on the Very Large Telescope Interferometer, we analysed the disc size, the time-variability of the disc's thermal emission, and the spectral shape of the $10~\mu$m silicate feature. By fitting simple disc models, we determined the inclination and the mid-infrared size of the disc, confirming earlier results based on a lower number of observations. We searched for mid-infrared temporal variability of different regions of the disc, and concluded that its morphology is not experiencing significant changes over the observed epochs. We characterized the mid-infrared silicate feature by determining the feature amplitude and the $11.3/9.8~\mu$m flux ratio. The latter parameter is a good indicator of the grain size. The shape of the feature suggests the presence of crystalline silicate grains in the disc. The interferometric data revealed a strong radial trend in the mineralogy: while the disc's innermost region seems to be dominated by forsterite grains, at intermediate radii both forsterite and enstatite may be present. The outer disc may predominantly contain amorphous silicate particles. A comparison of the observed spectral shape with that of a sample of intermediate-mass stars (supergiants, Herbig Ae/Be stars, unclassified B[e] stars) implied that the evolutionary state of HD 50138 cannot be unambiguously decided from mid-IR spectroscopy.

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.

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.

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.

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.

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).

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.

Context. Protoplanetary disks show large diversity regarding their morphology and dust composition. With mid-infrared interferometry the thermal emission of disks can be spatially resolved, and the distribution and properties of the dust within can be studied. Aims. Our aim is to perform a statistical analysis on a large sample of 82 disks around low- and intermediate-mass young stars, based on mid-infrared interferometric observations. We intend to study the distribution of disk sizes, variability, and the silicate dust mineralogy. Methods. Archival mid-infrared interferometric data from the MIDI instrument on the VLTI are homogeneously reduced and calibrated. Geometric disk models are used to fit the observations to get spatial information about the disks. An automatic spectral decomposition pipeline is applied to analyze the shape of the silicate feature. Results. We present the resulting data products in the form of an atlas, containing N band correlated and total spectra, visibilities, and differential phases. The majority of our data can be well fitted with a continuous disk model, except for a few objects, where a gapped model gives a better match. From the mid-infrared size--luminosity relation we find that disks around T Tauri stars are generally colder and more extended with respect to the stellar luminosity than disks around Herbig Ae stars. We find that in the innermost part of the disks ($r \lesssim 1$~au) the silicate feature is generally weaker than in the outer parts, suggesting that in the inner parts the dust is substantially more processed. We analyze stellar multiplicity and find that in two systems (AB Aur and HD 72106) data suggest a new companion or asymmetric inner disk structure. We make predictions for the observability of our objects with the upcoming MATISSE instrument, supporting the practical preparations of future MATISSE observations of T Tauri stars.

We present a detailed analysis of 30 very long baseline interferometric observations of the BL Lac object S5 1803+784 (z=0.679), obtained between mean observational time 1994.67 and 2012.91 at observational frequency 15 GHz. The long-term behaviour of the jet ridge line reveals the jet experiences an oscillatory motion superposed on its helical jet kinematics on a time-scale of about 6 years. The excess variance of the positional variability indicates the jet components being farther from the VLBI core have larger amplitude in their position variations. The fractional variability amplitude shows slight changes in 3-year bins of the component's position. The temporal variability in the Doppler boosting of the ridge line results in jet regions behaving as flaring "radio lanterns". We offer a qualitative scenario leading to the oscillation of the jet ridge line, that utilizes the orbital motion of the jet emitter black hole due to a binary black hole companion. A correlation analysis implies composite origin of the flux variability of the jet components, emerging due to possibly both the evolving jet-structure and its intrinsic variability.

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)

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.

Context. DG Tau is a low-mass pre-main sequence star, whose strongly accreting protoplanetary disk exhibits a so-far enigmatic behavior: its mid-infrared thermal emission is strongly time-variable, even turning the 10 $\mu$m silicate feature from emission to absorption temporarily. Aims. We look for the reason for the spectral variability at high spatial resolution and at multiple epochs. Methods. We study the temporal variability of the mid-infrared interferometric signal, observed with the VLTI/MIDI instrument at six epochs between 2011 and 2014. We fit a geometric disk model to the observed interferometric signal to obtain spatial information about the disk. We also model the mid-infrared spectra by template fitting to characterize the profile and time dependence of the silicate emission. We use physically motivated radiative transfer modeling to interpret the mid-infrared interferometric spectra. Results. The inner disk (r<1-3 au) spectra exhibit a 10 $\mu$m absorption feature related to amorphous silicate grains. The outer disk (r>1-3 au) spectra show a crystalline silicate feature in emission, similar to the spectra of comet Hale-Bopp. The striking difference between the inner and outer disk spectral feature is highly unusual among T Tauri stars. The mid-infrared variability is dominated by the outer disk. The strength of the silicate feature changed by more than a factor of two. Between 2011 and 2014 the half-light radius of the mid-infrared-emitting region decreased from 1.15 to 0.7 au. Conclusions. For the origin of the absorption we discuss four possible explanations: a cold obscuring envelope, an accretion heated inner disk, a temperature inversion on the disk surface and a misaligned inner geometry. The silicate emission in the outer disk can be explained by dusty material high above the disk plane, whose mass can change with time, possibly due to turbulence in the disk.

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.

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.

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.

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.

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.

Very Long Baseline Interferometry (VLBI) allows for high-resolution and high-sensitivity observations of relativistic jets, that can reveal periodicities of several years in their structure. We perform an analysis of long-term VLBI data of the quasar S5 1928+738 in terms of a geometric model of a helical structure projected onto the plane of the sky. We monitor the direction of the jet axis through its inclination and position angles. We decompose the variation of the inclination of the inner 2 milliarcseconds of the jet of S5 1928+738 into a periodic term with amplitude of ~0.89 deg and a linear decreasing trend with rate of ~0.05 deg/yr. We also decompose the variation of the position angle into a periodic term with amplitude of ~3.39 deg and a linear increasing trend with rate of ~0.24 deg/yr. We interpret the periodic components as arising from the orbital motion of a binary black hole inspiraling at the jet base and derive corrected values of the mass ratio and separation from the accumulated 18 years of VLBI data. Then we identify the linear trends in the variations as due to the slow reorientation of the spin of the jet emitter black hole induced by the spin-orbit precession and we determine the precession period T_SO=4852+/-646 yr of the more massive black hole, acting as the jet emitter. Our study provides indications, for the first time from VLBI jet kinematics, for the spinning nature of the jet-emitting black hole.

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.

Circumstellar disks around stars older than 10 Myr are expected to be gas-poor. There are, however, two examples of old (30-40 Myr) debris-like disks containing a detectable amount of cold CO gas. Here we present ALMA and Herschel Space Observatory observations of one of these disks, around HD 21997, and study the distribution and origin of the dust and its connection to the gas. Our ALMA continuum images at 886um clearly resolve a broad ring of emission within a diameter of ~4.5 arcsec, adding HD 21997 to the dozen debris disks resolved at (sub)millimeter wavelengths. Modeling the morphology of the ALMA image with a radiative transfer code suggests inner and outer radii of ~55 and ~150 AU, and a dust mass of 0.09 M_Earth. Our data and modeling hints at an extended cold outskirt of the ring. Comparison with the morphology of the CO gas in the disk reveals an inner dust-free hole where gas nevertheless can be detected. Based on dust grain lifetimes, we propose that the dust content of this gaseous disk is of secondary origin produced by planetesimals. Since the gas component is probably primordial, HD 21997 is one of the first known examples of a hybrid circumstellar disk, a so-far little studied late phase of circumstellar disk evolution.

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 BLLac object S4 0954+65 is one of the main targets of the Urumqi monitoring program targeting IntraDay Variable (IDV) sources. Between August 2005 and December 2009, the source was included in 41 observing sessions, carried out at a frequency of 4.8 GHz. The time analysis of the collected light curves, performed by applying both a structure function analysis and a specifically developed wavelet-based algorithm, discovered an annual cycle in the variability timescales, suggesting that there is a fundamental contribution by interstellar scintillation to the IDV pattern of the source. The combined use of the two analysis methods also revealed that there was a dramatic change in the variability characteristics of the source between February and March 2008, at the starting time of a strong outburst phase. The analysis' results suggest that the flaring state of the source coincides with the appearance of multiple timescales in its light curves, indicating that changes in the structure of the relativistically moving emitting region may strongly influence the variability observed on IDV timescales.

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.

A sample of classical IntraDay Variable (IDV) and IDV candidate sources has been monitored with the Urumqi 25m telescope and the Effelsberg 100m telescope. Aim of the project is to investigate the origin of IntraDay Variability, a phenomenon which has been observed in about 30% of flat spectrum radio quasars. Simultaneous Effelsberg-Urumqi observations demonstrated that the Urumqi antenna, although relatively small in diameter, is well suitable for IDV experiments. A few Urumqi datasets, however, turned out to be affected by a spurious $\sim 24$ hours modulation, an effect which has been removed by means of a new procedure for data reduction. In about 14 months, 12 epochs of observation have been collected, for a total observing time of more than 45 days. The epochs are regularly distributed over the whole year, in order to check for the presence of systematic annual changes in the variability time scales - a crucial test for verifying the consistency of source-extrinsic models of the variability. Preliminary time-analysis of the monitored sources revealed some hint for a slowing down of the characteristic time scales of S5~0716+714, a result that, if confirmed, would be compatible with a source-extrinsic origin of the variability, in contrast to previous IDV studies. No significant modulation of the time scales has been detected for S4~0954+658.

In radio bands, the study of compact radio sources can be affected by propagation effects through the Interstellar Medium. These are usually attributed to the presence of turbulent intervening plasma along the line of sight. Here, two of such effects are presented. The line of sight of B 2005+403 passes through the heavily scattered region of Cygnus, which causes substantial angular broadening of the source images obtained at frequencies between 0.6 GHz and 8 GHz. At higher frequencies, however the intrinsic source structure shines through. Therefore, multi-frequency VLBI observations allow to study the characteristics of the intervening material, the source morphology and the "coupling" of them in forming the observed image.

The BL Lac object S5 0716+71 was target of a coordinated and global multi-frequency campaign to search for rapid and correlated variability and signatures of the inverse-Compton catastrophe. Here we present first results obtained from a combined analysis of the cm- to sub-mm observations over a period of seven days aiming at a detailed study of the intra- to inter-day variability characteristics and to obtain constraints on the variability brightness temperatures and Doppler factors comparing the radio data with the high energy emission recorded by INTEGRAL. A more detailed description of the whole cm- to sub-mm observations and our analysis/results will be presented in a forthcoming paper. Our analysis reveals the source to be in a particular short-term variability phase when compared to the past with a correlated >~4 day time scale amplitude increase of up to 35%, which is systematically more pronounced towards higher frequencies. The obtained frequency dependent variability amplitudes and time lags contradict expectations from interstellar scintillation and strongly suggest a source intrinsic origin of this inter-day variability. A 7-day spectral evolution study indicate time-variable synchrotron self-absorption and expansion of the emission region, consistent with standard models. Assuming relativistic boosting, our different estimates of the Doppler factor yield robust lower limits of D_var,IC > 5-22 using the inverse-Compton limit and D_var,eq > 8-33 using the equipartition argument. Although high, these values are in good agreement with Doppler factors obtained from recent VLBI studies and from the inverse-Compton Doppler factors D_IC > 14-16 derived with the X-ray emission seen by INTEGRAL at 3-200 keV.