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.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

We present RadioAstron Space VLBI imaging observations of the BL Lac object S5 0716+71 made on January 3-4 2015 at a frequency of 22 GHz (wavelength $\lambda=1.3$ cm). The observations were made in the framework of the AGN Polarization Key Science Program. The source was detected on projected space-ground baselines up to 70 833 km (5.6 Earth diameters) for both, parallel hand and cross-hand interferometric visibilities. We have used these detections to obtain a full-polarimetric image of the blazar at an unprecedented angular resolution of 24 $\mu$as, the highest for this source to date. This enabled us to estimate the size of the radio core to be $<12\times5~\mu$as and to reveal a complex structure and a significant curvature of the blazar jet in the inner 100 $\mu$as, which is an indication that the jet viewing angle lies inside the opening angle of the jet conical outflow. Fairly highly (15%) linearly polarized emission is detected in a jet region of 19 $\mu$as in size, located 58 $\mu$as downstream from the core. The highest brightness temperature in the source frame is estimated to be $>2.2\times10^{13}$ K for the blazar core. This implies that the inverse Compton limit must be violated in the rest frame of the source, even for the largest Doppler factor $\delta\thicksim25$ reported for 0716+714.

Detailed studies of relativistic jets in active galactic nuclei (AGN) require high-fidelity imaging at the highest possible resolution. This can be achieved using very long baseline interferometry (VLBI) at radio frequencies, combining worldwide (global) VLBI arrays of radio telescopes with a space-borne antenna on board a satellite. We present multiwavelength images made of the radio emission in the powerful quasar S5 0836+710, obtained using a global VLBI array and the antenna Spektr-R of the RadioAstron mission of the Russian Space Agency, with the goal of studying the internal structure and physics of the relativistic jet in this object. The RadioAstron observations at wavelengths of 18cm, 6cm, and 1.3cm are part of the Key Science Program for imaging radio emission in strong AGN. The internal structure of the jet is studied by analyzing transverse intensity profiles and modeling the structural patterns developing in the flow. The RadioAstron images reveal a wealth of structural detail in the jet of S5 0836+710 on angular scales ranging from 0.02mas to 200mas. Brightness temperatures in excess of $10^{13}$\,K are measured in the jet, requiring Doppler factors of $\ge 100$ for reconciling them with the inverse Compton limit. Several oscillatory patterns are identified in the ridge line of the jet and can be explained in terms of the Kelvin-Helmholtz (KH) instability. The oscillatory patterns are interpreted as the surface and body wavelengths of the helical mode of the KH instability. The interpretation provides estimates of the jet Mach number and of the ratio of the jet to the ambient density, which are found to be $M_\mathrm{j}\approx 12$ and $\eta\approx 0.33$. The ratio of the jet to the ambient density should be conservatively considered an upper limit because its estimate relies on approximations.

The origin of non-baryonic Dark Matter remains elusive despite ongoing sensitive searches for heavy, thermally produced dark matter particles. Recently, it has been shown, that non-thermally produced vector bosons (sometimes called hidden photons) related to a broken U(1) gauge symmetry are among the possible WISP (weakly interacting slim particles) dark matter candidates. The WISP Dark Matter eXperiment (WISPDMX) is the first direct hidden photon dark matter search experiment probing the particle masses within the 0.8-2.07 $\mu$eV range with four resonant modes of a tunable radio frequency cavity and down to 0.4 $n$eV outside of resonance. In this paper, we present the results from the first science run of WISPDMX comprising 22000 broadband spectra with a 500 MHz bandwidth and a 50 Hz spectral resolution, obtained during 10-second integrations made at each individual tuning step of the measurements. No plausible dark matter candidate signal is found, both in the individual spectra reaching minimum detectable power of $8^{-19}$ Watt and in the averaged spectrum of all the measurements with the minimum detectable power of $5^{-22}$ Watt attained for a total of 61 \hours of data taking. Using these spectra, we derive upper limits on the coupling constant of the hidden photon at the levels of $10^{e-13}$ for the resonant frequency ranges and $10^{-12}$ for broadband mass range 0.2-2.07 $\mu$eV, and steadily increasing at masses below 0.2$\mu$eV.

The astrophysical agenda of the 21st century requires a very sharp view of celestial objects. High angular resolution studies are essential for fundamental studies of a broad variety of astrophysical phenomena ranging from relativistic physics of black holes, their gravitational and electromagnetic imprints, violent transient processes, including those producing detectable gravitational waves, birth and evolution of planetary systems. Over the past decades, radio astronomy made huge leap in achieving ground-breaking angular resolution measured in tens of microarcseconds (one tenth of nanoradian and better). Recently a global Event Horizon Telescope (EHT) collaboration obtained first direct images of the shadow of a super-massive black hole in the nucleus of the active galaxy M87. These observations were conducted at 230 GHz. The two first generation Space Very Long Baseline Interferometry (VLBI) missions, VSOP/HALCA led by the Japan Aerospace Exploration Agency (JAXA) and RadioAstron led by the Russia Roscosmos State Corporation and Russia Academy of Sciences, achieved the highest angular resolution at frequencies from 0.3 to 22 GHz in observations conducted in the period 1997 - 2019. The next step in advancing high angular resolution radio astronomy is in combining high frequency (millimeter and sub-millimeter wavelengths) and interferometric baselines exceeding the Earth diameter. The present THEZA White Paper describes a combination which would unify technology developments in giga-/tera-hertz instrumentation and space-borne radio astronomy. The current preprint version of the THEZA White Paper is slightly re-formatted and edited comparing to the official submitted version.

We investigate the relationship between the variable gamma-ray emission and jet properties in the blazar 3C 279, by combining the Fermi/LAT data spanning a period of eight years with concurrent radio measurements made at multiple epochs with VLBA at 15 and 43 GHz within the MOJAVE and VLBA-BU monitoring programs. The aim of this paper is to compare the flux variability of the different components found in the VLBA observations, to the variability in the gamma-rays. This analysis helps to investigate whether any of the jet components can be associated with the gamma-ray variability. Through Spearman rank correlation we found that the gamma-ray variability is correlated with a particular region (feature B, in the MOJAVE images) downstream from the observed base (core) of the jet. This jet component is, therefore, a likely location at which an important fraction of the variable gamma-ray emission is produced. We also calculated the average proper motion of the component with respect to the VLBA core and found that it moves at an apparent superluminal velocity of $(3.70 \pm 0.35)c$, implying that one of the gamma-ray emission zones is not stationary. This jet component is also found between 6.86 mas and 8.68 mas, which translates to a distance from the radio core of at least 42 pc.

A number of extragalactic jets show periodic structures at different scales that can be associated with growing instabilities. The wavelengths of the developing instability modes and their ratios depend on the flow parameters, so the study of those structures can shed light on jet physics at the scales involved. In this work, we use the fits to the jet ridgeline obtained from different observations of S5 B0836$+$710 and apply stability analysis of relativistic, sheared flows to derive an estimate of the physical parameters of the jet. Based on the assumption that the observed structures are generated by growing Kelvin-Helmholtz (KH) instability modes, we have run numerical calculations of stability of a relativistic, sheared jet over a range of different jet parameters. We have spanned several orders of magnitude in jet-to-ambient medium density ratio, and jet internal energy, and checked different values of the Lorentz factor and shear layer width. This represents an independent method to obtain estimates of the physical parameters of a jet. By comparing the fastest growing wavelengths of each relevant mode given by the calculations with the observed wavelengths reported in the literature, we have derived independent estimates of the jet Lorentz factor, specific internal energy, jet-to-ambient medium density ratio and Mach number. We obtain a jet Lorentz factor $\gamma \simeq 12$, specific internal energy of $\varepsilon \simeq 10^{-2}\,c^2$, jet-to-ambient medium density ratio of $\eta\approx 10^{-3}$, and an internal (classical) jet Mach number of $M_\mathrm{j}\approx 12$. We also find that the wavelength ratios are better recovered by a transversal structure with a width of $\simeq 10\,\%$ of the jet radius. This method represents a powerful tool to derive the jet parameters in all jets showing helical patterns with different wavelengths.

Imaging relativistic jets in active galactic nuclei (AGN) at angular resolution significantly surpassing that of the ground-based VLBI at centimetre wavelengths is one of the key science objectives of the RadioAstron space-VLBI mission. There are three RadioAstron imaging key science programs that target both nearby radio galaxies and blazars, with one of the programs specifically focusing on polarimetry of the jets. The first images from these programs reach angular resolution of a few tens of microarcseconds and reveal unprecedented details about the jet collimation profile, magnetic field configuration, and Kelvin-Helmholtz instabilities along the flow in some of the most studied AGN (3C 84, BL Lac, 3C 273, S5 0836+710). Here we give an overview of the goals and strategy of these three ongoing programs, highlight their early results, and discuss the challenges of space-VLBI imaging.

The results of MHD numerical simulations of the formation and development of magnetized jets are presented. Similarity criteria for comparisons of the results of laboratory laser experiments and numerical simulations of astrophysical jets are discussed. The results of laboratory simulations of jets generated in experiments at the Neodim laser installation are presented.

The apparent position of jet base (core) in radio-loud active galactic nuclei changes with frequency because of synchrotron self-absorption. Studying this `core shift` effect enables us to reconstruct properties of the jet regions close to the central engine. We report here results from core shift measurements in AGNs observed with global VLBI at 2 and 8 GHz at epochs from 1994 to 2016. Our sample contains 40 objects observed at least 10 times during that period. The core shift is determined using a new automatic procedure introduced to minimize possible biases. The resulting multiple epoch measurements of the core position are employed for examining temporal variability of the core shift. We argue that the core shift variability is a common phenomenon, as established for 33 of 40 AGNs we study. Our analysis shows that the typical offsets between the core positions at 2 and 8 GHz are about 0.5 mas and they vary in time. Typical variability of the individual core positions is about 0.3 mas. The measurements show a strong dependence between the core position and its flux density, suggesting that changes in both are likely related to the nuclear flares injecting denser plasma into the flow. We determine that density of emitting relativistic particles significantly increases during these flares, while relative magnetic field changes less and in the opposite direction.

Results of a pilot project with the participation of the "Kvazar-KVO" radio interferometry array in observations carried out with the European VLBI Network are presented. The aim of the project was to conduct and analyze multi-frequency (1.7, 2.3, 5.0, 8.4 GHz) observations of the parsec-scale jets of 24 active galactic nuclei. Three observing sessions were successfully carried out in October 2008. Maps of the radio intensity distributions have been constructed in all four frequencies using phase referencing. A method for measuring the frequency-dependent shift of the position of the VLBI core by applying relative astrometry to observations of close triplets of radio sources has been developed. The fundamental possibility of detecting core shifts in ultra-compact sources for which traditional methods based on the achromatic positions of optically thin regions of the jet are not suitable is demonstrated. The conditions for successful measurement of this shift are discussed; these are determined by the closeness of the calibrator used, the effective resolution of the system, the quality of the filling of the $uv$ plane, the relative orientations of the jets in the triplets, and the brightnesses of the sources.

(abridged) Very long baseline interferometry (VLBI) observations at 86$\,$GHz (wavelength, $\lambda = 3\,$mm) reach a resolution of about 50 $\mu$as, probing the collimation and acceleration regions of relativistic outflows in active galactic nuclei. To extend the statistical studies of compact extragalactic jets, a large global 86 GHz VLBI survey of 162 radio sources was conducted in 2010-2011 using the Global Millimeter VLBI Array. The survey data attained a typical baseline sensitivity of 0.1 Jy and a typical image sensitivity of 5 mJy/beam, providing successful detections and images for all of the survey targets. For 138 objects, the survey provides the first ever VLBI images made at 86 GHz. Gaussian model fitting of the visibility data was applied to represent the structure of the sources. The Gaussian model-fit-based estimates of brightness temperature ($T_\mathrm{b}$) at the jet base (core) and in moving regions (jet components) downstream from the core were compared to the estimates of $T_\mathrm{b}$ limits made directly from the visibility data, demonstrating a good agreement between the two methods. The apparent brightness temperature estimates for the jet cores in our sample range from $2.5 \times 10^{9}\,$K to $ 1.3\times 10^{12}\,$K. A population model with a single intrinsic value of brightness temperature, $T_\mathrm{0}$, is applied to reproduce the observed $T_\mathrm{b}$ distribution. It yields $T_\mathrm{0} = (3.77^{+0.10}_{-0.14}) \times 10^{11}\,$K for the jet cores, implying that the inverse Compton losses dominate the emission. In the jet components, $T_\mathrm{0} =(1.42^{+0.16}_{-0.19})\times 10^{11}\,$K is found, slightly higher than the equipartition limit of $\sim5\times 10^{10}\,$K expected for these jet regions. For objects with sufficient structural detail detected, the adiabatic energy losses dominate the observed changes of $T_\mathrm{b}$ along the jet.

Understanding the launching, acceleration, and collimation of jets powered by active galactic nuclei remains an outstanding problem in relativistic astrophysics. This is partly because observational tests of jet formation models suffer from the limited angular resolution of ground-based very long baseline interferometry that has thus far been able to probe the transverse jet structure in the acceleration and collimation zone of only two sources. Here we report radio interferometric observations of 3C 84 (NGC 1275), the central galaxy of the Perseus cluster, made with an array including the orbiting radio telescope of the RadioAstron mission. The obtained image transversely resolves the edge-brightened jet in 3C 84 only 30 microarcseconds from the core, which is ten times closer to the central engine than what has been possible in previous ground-based observations, and it allows us to measure the jet collimation profile from ~ 100 to ~10000 gravitational radii from the black hole. The previously found, almost cylindrical jet profile on scales larger than a few thousand r_g is now seen to continue at least down to a few hundred r_g from the black hole and we find a broad jet with a transverse radius larger than about 250 r_g at only 350 r_g from the core. If the bright outer jet layer is launched by the black hole ergosphere, it has to rapidly expand laterally on scales smaller than 100 r_g. If this is not the case, then this jet sheath is likely launched from the accretion disk.

Clues to the physical conditions in radio cores of blazars come from measurements of brightness temperatures as well as effects produced by intrinsic opacity. We study the properties of the ultra compact blazar AO 0235+164 with RadioAstron ground-space radio interferometer, multi-frequency VLBA, EVN and single-dish radio observations. We employ visibility modeling and image stacking for deriving structure and kinematics of the source, and use Gaussian process regression to find the relative multi-band time delays of the flares. The multi-frequency core size and time lags support prevailing synchrotron self absorption. The intrinsic brightness temperature of the core derived from ground-based VLBI is close to the equipartition regime value. In the same time, there is evidence for ultra-compact features of the size of less than 10 $\mu$as in the source, which might be responsible for the extreme apparent brightness temperatures of up to $10^{14}$ K as measured by RadioAstron. In 2007--2016 the VLBI components in the source at 43 GHz are found predominantly in two directions, suggesting a bend of the outflow from southern to northern direction. The apparent opening angle of the jet seen in the stacked image at 43 GHz is two times wider than that at 15 GHz, indicating a collimation of the flow within the central 1.5 mas. We estimate the Lorentz factor $\Gamma = 14$, the Doppler factor $\delta=21$, and the viewing angle $\theta = 1.7^\circ$ of the apparent jet base, derive the gradients of magnetic field strength and electron density in the outflow, and the distance between jet apex and the core at each frequency.

RadioAstron is a 10 m orbiting radio telescope mounted on the Spektr-R satellite, launched in 2011, performing Space Very Long Baseline Interferometry (SVLBI) observations supported by a global ground array of radio telescopes. With an apogee of about 350 000 km, it is offering for the first time the possibility to perform {\mu}as-resolution imaging in the cm-band. We present observations at 22 GHz of 3C 273, performed in 2014, designed to reach a maximum baseline of approximately nine Earth diameters. Reaching an angular resolution of 0.3 mas, we study a particularly low-activity state of the source, and estimate the nuclear region brightness temperature, comparing with the extreme one detected one year before during the RadioAstron early science period. We also make use of the VLBA-BU-BLAZAR survey data, at 43 GHz, to study the kinematics of the jet in a 1.5-year time window. We find that the nuclear brightness temperature is two orders of magnitude lower than the exceptionally high value detected in 2013 with RadioAstron at the same frequency (1.4x10^13 K, source-frame), and even one order of magnitude lower than the equipartition value. The kinematics analysis at 43 GHz shows that a new component was ejected 2 months after the 2013 epoch, visible also in our 22 GHz map presented here. Consequently this was located upstream of the core during the brightness temperature peak. These observations confirm that the previously detected extreme brightness temperature in 3C 273, exceeding the inverse Compton limit, is a short-lived phenomenon caused by a temporary departure from equipartition. Thus, the availability of interferometric baselines capable of providing {\mu}as angular resolution does not systematically imply measured brightness temperatures over the known physical limits for astrophysical sources.

M87 is the best source in which to study a jet at high resolution in gravitational units because it has a very high mass black hole and is nearby. The angular size of the black hole is second only to Sgr A*, which does not have a strong jet. The jet structure is edge brightened with a wide opening angle base and a weak counterjet. We have roughly annual observations for 17 years plus intensive monitoring at three week intervals for a year and five day intervals for 2.5 months made with the Very Long Baseline Array (VLBA) at 43 GHz. The inner jet shows very complex dynamics, with apparent motions both along and across the jet. Speeds from zero to over 2c are seen, with acceleration observed over the first 3 milli-arcseconds. The counterjet decreases in brightness much more rapidly than the main jet, as is expected from relativistic beaming in an accelerating jet oriented near the line-of-sight. Details of the structure and dynamics are discussed. The roughly annual observations show side-to-side motion of the whole jet with a characteristic time scale of about 9 years.

Very long baseline interferometry (VLBI) imaging of radio emission from extragalactic jets provides a unique probe of physical mechanisms governing the launching, acceleration, and collimation of relativistic outflows. The two-dimensional structure and kinematics of the jet in M\,87 (NGC\,4486) have been studied by applying the Wavelet-based Image Segmentation and Evaluation (WISE) method to 11 images obtained from multi-epoch Very Long Baseline Array (VLBA) observations made in January-August 2007 at 43 GHz ($\lambda = 7$ mm). The WISE analysis recovers a detailed two-dimensional velocity field in the jet in M\,87 at sub-parsec scales. The observed evolution of the flow velocity with distance from the jet base can be explained in the framework of MHD jet acceleration and Poynting flux conversion. A linear acceleration regime is observed up to $z_{obs} \sim 2$\,mas. The acceleration is reduced at larger scales, which is consistent with saturation of Poynting flux conversion. Stacked cross correlation analysis of the images reveals a pronounced stratification of the flow. The flow consists of a slow, mildly relativistic layer (moving at $\beta \sim 0.5\,c$), associated either with instability pattern speed or an outer wind, and a fast, accelerating stream line (with $\beta \sim 0.92$, corresponding to a bulk Lorentz factor $\gamma \sim 2.5$). A systematic difference of the apparent speeds in the northern and southern limbs of the jet is detected, providing evidence for jet rotation. The angular velocity of the magnetic field line associated with this rotation suggests that the jet in M87 is launched in the inner part of the disk, at a distance $r_0 \sim 5\, R_\mathrm{s}$ from the central engine. The combined results of the analysis imply that MHD acceleration and conversion of Poynting flux to kinetic energy play the dominant roles in collimation and acceleration of the flow in M\,87.

Supermassive black holes (SMBH) are essential for the production of jets in radio-loud active galactic nuclei (AGN). Theoretical models based on Blandford & Znajek extract the rotational energy from a Kerr black hole, which could be the case for NGC1052, to launch these jets. This requires magnetic fields of the order of $10^3\,$G to $10^4\,$G. We imaged the vicinity of the SMBH of the AGN NGC1052 with the Global Millimetre VLBI Array and found a bright and compact central feature, smaller than 1.9 light days (100 Schwarzschild radii) in radius. Interpreting this as a blend of the unresolved jet bases, we derive the magnetic field at 1 Schwarzschild radius to lie between 200 G and ~80000 G consistent with Blandford & Znajek models.

Jets of compact radio sources are highly relativistic and Doppler boosted, making studies of their intrinsic properties difficult. Observed brightness temperatures can be used to study the intrinsic physical properties of the relativistic jets, and constrain models of jet formation in the inner jet region. We aim to observationally test such inner jet models. The very long baseline interferometry (VLBI) cores of compact radio sources are optically thick at a given frequency. The distance of the core from the central engine is inversely proportional to the frequency. Under the equipartition condition between the magnetic field energy and particle energy densities, the absolute distance of the VLBI core can be predicted. We compiled the brightness temperatures of VLBI cores at various radio frequencies of 2, 8, 15, and 86~GHz. We derive the brightness temperature on sub-parsec scales in the rest frame of the compact radio sources. We find that the brightness temperature increases with increasing distance from the central engine, indicating that the intrinsic jet speed (the Lorentz-factor) increases along the jet. This implies that the jets are accelerated in the (sub-)parsec regions from the central engine.

Inverse Compton cooling limits the brightness temperature of the radiating plasma to a maximum of $10^{11.5}$ K. Relativistic boosting can increase its observed value, but apparent brightness temperatures much in excess of $10^{13}$ K are inaccessible using ground-based very long baseline interferometry (VLBI) at any wavelength. We present observations of the quasar 3C273, made with the space VLBI mission RadioAstron on baselines up to 171,000 km, which directly reveal the presence of angular structure as small as 26 $\mu$as (2.7 light months) and brightness temperature in excess of $10^{13}$ K. These measurements challenge our understanding of the non-thermal continuum emission in the vicinity of supermassive black holes and require a much higher Doppler factor than what is determined from jet apparent kinematics.

Velocity measurements made from multiple-epoch astronomical images of evolving objects with optically thin continuum emission (e.g. as relativistic jets or expanding supernova shells) may be confused as a result of the overlap of semi-transparent features moving at different speeds. Multi-scale wavelet decomposition can be effectively applied to identify and track such overlapping features, provided that their respective structural responses can be separated over the spatial scales used for the decomposition. We developed a new method that combines the stacked cross-correlation with the wavelet-based image segmentation and evaluation (WISE) technique of decomposition of two-dimensional structures, to separate and track dominant spatial responses of overlapping evolving features. The method is tested on a set of simulated images of a stratified relativistic jet, demonstrating the robust detection of both the faster spine and the slower sheath speeds. The method is applied to mutliple-epoch images from the MOJAVE survey, revealing two different superluminal streams inside the jet in 3C 273 and the acceleration of the flow in 3C 120. The method can be applied to densely monitored objects with composite structural evolution such as the parsec-scale jet in M87 or heavily resolved expanding supernova shells.

We present the first polarimetric space VLBI imaging observations at 22 GHz. BL Lacertae was observed in 2013 November 10 with the RadioAstron space VLBI mission, including a ground array of 15 radio telescopes. The instrumental polarization of the space radio telescope is found to be within 9%, demonstrating the polarimetric imaging capabilities of RadioAstron at 22 GHz. Ground-space fringes were obtained up to a projected baseline distance of 7.9 Earth's diameters in length, allowing us to image the jet in BL Lacertae with a maximum angular resolution of 21 $\mu$as, the highest achieved to date. We find evidence for emission upstream of the radio core, which may correspond to a recollimation shock at about 40 $\mu$as from the jet apex, in a pattern that includes other recollimation shocks at approximately 100 $\mu$as and 250 $\mu$as from the jet apex. Polarized emission is detected in two components within the innermost 0.5 mas from the core, as well as in some knots 3 mas downstream. Faraday rotation analysis, obtained from combining RadioAstron 22 GHz and ground-based 15 GHz and 43 GHz images, shows a gradient in rotation measure and Faraday corrected polarization vector as a function of position angle with respect to the core, suggesting that the jet in BL Lacertae is threaded by a helical magnetic field. The intrinsic de-boosted brightness temperature in the unresolved core exceeds $3\!\times\!10^{12}$ K, suggesting at the very least departure from equipartition of energy between the magnetic field and radiating particles.

Polarization of radio emission in extragalactic jets at a sub-milliarcsecond angular resolution holds important clues for understanding the structure of the magnetic field in the inner regions of the jets and in close vicinity of the supermassive black holes in the centers of active galaxies. Space VLBI observations provide a unique tool for polarimetric imaging at a sub-milliarcsecond angular resolution and studying the properties of magnetic field in active galactic nuclei on scales of less than 10^4 gravitational radii. A space VLBI observation of high-redshift quasar TXS 0642+449 (OH 471), made at a wavelength of 18 cm (frequency of 1.6 GHz) as part of the Early Science Programme (ESP) of the RadioAstron} mission, is used here to test the polarimetric performance of the orbiting Space Radio Telescope (SRT) employed by the mission, to establish a methodology for making full Stokes polarimetry with space VLBI at 1.6 GHz, and to study the polarized emission in the target object on sub-milliarcsecond scales. Polarization leakage of the SRT at 18 cm is found to be within 9 percents in amplitude, demonstrating the feasibility of high fidelity polarization imaging with RadioAstron at this wavelength. A polarimetric image of 0642+449 with a resolution of 0.8 mas (signifying an ~4 times improvement over ground VLBI observations at the same wavelength) is obtained. The image shows a compact core-jet structure with low (~2%) polarization and predominantly transverse magnetic field in the nuclear region. The VLBI data also uncover a complex structure of the nuclear region, with two prominent features possibly corresponding to the jet base and a strong recollimation shock. The maximum brightness temperature at the jet base can be as high as 4*10^13 K.

Opacity-driven shifts of the apparent VLBI core position with frequency (the "core shift" effect) probe physical conditions in the innermost parts of jets in active galactic nuclei. We present the first detailed investigation of this effect in the brightest gamma-ray blazar 3C454.3 using direct measurements from simultaneous 4.6-43 GHz VLBA observations, and a time lag analysis of 4.8-37 GHz lightcurves from the UMRAO, CrAO, and Metsahovi observations in 2007-2009. The results support the standard Konigl model of jet physics in the VLBI core region. The distance of the core from the jet origin r_c(nu), the core size W(nu), and the lightcurve time lag DT(nu) all depend on the observing frequency nu as r_c(nu)~W(nu)~ DT(nu)~nu^-1/k. The obtained range of k=0.6-0.8 is consistent with the synchrotron self-absorption being the dominating opacity mechanism in the jet. The similar frequency dependence of r_c(nu) and W(nu) suggests that the external pressure gradient does not dictate the jet geometry in the cm-band core region. Assuming equipartition, the magnetic field strength scales with distance r as B = 0.4(r/1pc)^-0.8 G. The total kinetic power of electron/positron jet is about 10^44 ergs/s.

Context. Images of spatially resolved astrophysical objects contain a wealth of morphological and dynamical information, and effective extraction of this information is of paramount importance for understanding the physics and evolution of these objects. Algorithms and methods employed presently for this purpose (such as, for instance, Gaussian model fitting) often use simplified approaches for describing the structure of resolved objects. Aims. Automated (unsupervised) methods for structure decomposition and tracking of structural patterns are needed for this purpose, in order to be able to deal with the complexity of structure and large amount of data involved. Methods. A new Wavelet-based Image Segmentation and Evaluation (WISE) method is developed for multiscale decomposition, segmentation, and tracking of structural patterns in astronomical images. Results. The method is tested against simulated images of relativistic jets and applied to data from long-term monitoring of parsec- scale radio jets in 3C 273 and 3C 120. Working at its coarsest resolution, WISE reproduces exceptionally well the previous results of model fitting evaluation of the structure and kinematics in these jets. Extending the WISE structure analysis to the finer scales provides the first robust measurements of two-dimensional velocity fields in these jets and indicates that the velocity fields are likely to reflect the evolution of Kelvin-Hemlholtz instability developing the flow.

Radio-loud AGN typically manifest powerful relativistic jets extending up to millions of light years and often showing superluminal motions organised in a complex kinematic pattern. A number of physical models are still competing to explain the jet structure and kinematics revealed by radio images using the VLBI technique. Robust measurements of longitudinal and transverse velocity field in the jets would provide crucial information for these models. This is a difficult task, particularly for transversely resolved jets in objects like 3C 273 and M87. To address this task, we have developed a new technique for identifying significant structural patterns (SSP) of smooth, transversely resolved flows and obtaining a velocity field from cross-correlation of these regions in multi-epoch observations. Detection of individual SSP is performed using the wavelet decomposition and multiscale segmentation of the observed structure. The cross-correlation algorithm combines structural information on different scales of the wavelet decomposition, providing a robust and reliable identification of related SSP in multi -epoch images. The algorithm enables recovering structural evolution on scales down to 0.25 full width at half maximum (FWHM) of the image point spread function (PSF). We present here examples of applying this algorithm to obtain the first detailed transverse velocity fields and to study the kinematic evolution in the parsec-scale jets in 3C 273 and M87.

One of the most intriguing results from the gamma-ray instruments in orbit has been the detection of powerful flares from the Crab Nebula. These flares challenge our understanding of pulsar wind nebulae and models for particle acceleration. We report on the portion of a multiwavelength campaign using Keck, HST, and Chandra concentrating on a small emitting region, the Crab's inner knot, located a fraction of an arcsecond from the pulsar. We find that the knot's radial size, tangential size, peak flux, and the ratio of the flux to that of the pulsar are correlated with the projected distance of the knot from the pulsar. A new approach, using singular value decomposition for analyzing time series of images, was introduced yielding results consistent with the more traditional methods while some uncertainties were substantially reduced. We exploit the characterization of the knot to discuss constraints on standard shock-model parameters that may be inferred from our observations assuming the inner knot lies near to the shocked surface. These include inferences as to wind magnetization, shock shape parameters such as incident angle and poloidal radius of curvature, as well as the IR/optical emitting particle enthalpy fraction. We find that while the standard shock model gives good agreement with observation in many respects, there remain two puzzles: (a) The observed angular size of the knot relative to the pulsar--knot separation is much smaller than expected; (b) The variable, yet high degree of polarization reported is difficult to reconcile with a highly relativistic downstream flow.

The radio galaxy 3C84 is a representative of gamma-ray-bright misaligned active galactic nuclei (AGN) and one of the best laboratories to study the radio properties of subparsec scale jets. We discuss here the past and present activity of the nuclear region within the central 1pc and the properties of subparsec-sized components C1, C2 and C3. We compare these results with the high resolution space-VLBI image at 5GHz obtained with the RadioAstron satellite and we shortly discuss the possible correlation of radio emission with the gamma-ray emission.

DiFX is a correlator for VLBI data based on the FX architecture (first Fourier transform and then cross-multiply). DiFX is a free licensed software written in C++, developed and maintened by an international group of programmers. A new DiFX version (dra) has been developed at Max-Planck-Institut f\"ur Radioastronomie (MPIfR), in order to manage the correlation of a space-based antenna with ground stations. The dra version is running on the High Performance Computer cluster (HPC) in Bonn, and it is used for the data processing of the three AGN imaging RadioAstron Key Science Projects ongoing, based at the MPIfR.

Jet ejection by accreting black holes is a mass invariant mechanism unifying stellar and supermassive black holes (SMBHs) that should also apply for intermediate-mass black holes (IMBHs), which are thought to be the seeds from which SMBHs form. We present the detection of an off-nuclear IMBH of $\sim$5 $\times$ 10$^{4}$ M$_\odot$ located in an unusual spiral arm of the galaxy NGC 2276 based on quasi-simultaneous \textit{Chandra} X-ray observations and European VLBI Network (EVN) radio observations. The IMBH, NGC2276-3c, possesses a 1.8 pc radio jet that is oriented in the same direction as large-scale ($\sim$650 pc) radio lobes and whose emission is consistent with flat to optically thin synchrotron emission between 1.6 GHz and 5 GHz. Its jet kinetic power ($4 \times 10^{40}$ erg s$^{-1}$) is comparable to its radiative output and its jet efficiency ($\geq$ 46\%) is as large as that of SMBHs. A region of $\sim$300 pc along the jet devoid of young stars could provide observational evidence of jet feedback from an IMBH. The discovery confirms that the accretion physics is mass invariant and that seed IMBHs in the early Universe possibly had powerful jets that were an important source of feedback.

In recent years cosmology has undergone a revolution, with precise measurements of the microwave background radiation, large galaxy redshift surveys, and the discovery of the recent accelerated expansion of the Universe using observations of distant supernovae. In this light, the SKA enables us to do an ultimate test in cosmology by measuring the expansion rate of the Universe in real time. This can be done by a rather simple experiment of observing the neutral hydrogen (HI) signal of galaxies at two different epochs. The signal will encounter a change in frequency imprinted as the Universe expands over time and thus monitoring the drift in frequencies will provide a real time measure of the cosmic acceleration. Over a period of 12 years one would expected a frequency shift of the order of 0.1 Hz assuming a standard Lambda-CDM cosmology. Based on the sensitivity estimates of the SKA and the number counts of the expected HI galaxies, it is shown that the number counts are sufficiently high to compensate for the observational uncertainties of the measurements and hence allow a statistical detection of the frequency shift. [abstract abridged]

Dark Matter (DM) is a fundamental ingredient of our Universe and of structure formation, and yet its nature is elusive to astrophysical probes. Information on the nature and physical properties of the WIMP (neutralino) DM (the leading candidate for a cosmologically relevant DM) can be obtained by studying the astrophysical signals of their annihilation/decay. Among the various e.m. signals, secondary electrons produced by neutralino annihilation generate synchrotron emission in the magnetized atmosphere of galaxy clusters and galaxies which could be observed as a diffuse radio emission (halo or haze) centered on the DM halo. A deep search for DM radio emission with SKA in local dwarf galaxies, galaxy regions with low star formation and galaxy clusters (with offset DM-baryonic distribution, like e.g. the Bullet cluster) can be very effective in constraining the neutralino mass, composition and annihilation cross-section. For the case of a dwarf galaxy, like e.g. Draco, the constraints on the DM annihilation cross-section obtainable with SKA1-MID will be at least a factor $\sim 10^3$ more stringent than the limits obtained by Fermi-LAT in the $\gamma$-rays. These limits scale with the value of the B field, and the SKA will have the capability to determine simultaneously both the magnetic field in the DM-dominated structures and the DM particle properties. The optimal frequency band for detecting the DM-induced radio emission is around $\sim 1$ GHz, with the SKA1-MID Band 1 and 4 important to probe the synchrotron spectral curvature at low-$\nu$ (sensitive to DM composition) and at high-$\nu$ (sensitive to DM mass).

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.

The brightness temperature is an effective parameter that describes the physical properties of emitting material in astrophysical objects. It is commonly determined by imaging and modeling the structure of the emitting region and estimating its flux density and angular size. Reliable approaches for visibility-based estimates of brightness temperature are needed for interferometric experiments in which poor coverage of spatial frequencies prevents successful imaging of the source structure, for example, in interferometric measurements made at millimeter wavelengths or with orbiting antennas. Such approaches can be developed by analyzing the relations between brightness temperature and visibility amplitude and its r.m.s. error. A method is introduced for directly calculating the lower and upper limits of the brightness temperature from visibility measurements. The visibility-based brightness temperature estimates are shown to agree well with the image-based estimates obtained in the 2\,cm MOJAVE survey and the 3\,mm CMVA survey, with good agreement achieved for interferometric measurements at spatial frequencies exceeding $\approx 2\times 10^8$. The method provides an essential tool for constraining brightness temperature in all interferometric experiments with poor imaging capability.

AGN reverberate when the broad emission lines respond to changes of the ionizing thermal continuum emission. Reverberation measurements have been commonly used to estimate the size of the broad-line region (BLR) and the mass of the central black hole. However, reverberation mapping studies have been mostly performed on radio-quiet sources where the contribution of the jet can be neglected. In radio-loud AGN, jets and outflows may affect substantially the relation observed between the ionizing continuum and the line emission. To investigate this relation, we have conducted a series of multiwavelength studies of radio-loud AGN, combining optical spectral line monitoring with regular VLBI observations. Our results suggest that at least a fraction of the broad-line emitting material can be located in a sub-relativistic outflow ionized by non-thermal continuum emission generated in the jet at large distances (> 1pc) from the central engine of AGN. This finding may have a strong impact on black hole mass estimates based on measured widths of the broad emission lines and on the gamma-ray emission mechanisms.

We present new Jansky Very Large Array (VLA) radio images of the Crab Nebula at 5.5 GHz, taken at two epochs separated by 6 days about two months after a gamma-ray flare in 2012 July. We find no significant change in the Crab's radio emission localized to a region of <2 light-months in radius, either over the 6-day interval between our present observations or between the present observations and ones from 2001. Any radio counterpart to the flare has a radio luminosity of <~ $2 \times 10^{-4}$ times that of the nebula. Comparing our images to one from 2001, we do however find changes in radio brightness, up to 10% in amplitude, which occur on decade timescales throughout the nebula. The morphology of the changes is complex suggesting both filamentary and knotty structures. The variability is stronger, and the timescales likely somewhat shorter, nearer the centre of the nebula. We further find that even with the excellent uv~coverage and signal-to-noise of the VLA, deconvolution errors are much larger than the noise, being up to 1.2% of peak brightness of the nebula in this particular case.

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

Several observational studies and numerical simulations suggest that mergers must contribute to the evolution of galaxies; however, the role that they play is not yet fully understood. In this paper we study a sample of 52 double nucleus disk galaxies that are considered as candidates for a minor merger event. The luminosity of each of the nuclei and their relative separation are derived from a multi-component photometric fit of the galaxies in the Sloan Digital Sky Survey optical images. We find that the nuclei in most of the sources have projected separations $\leq$ 4 kpc. The ratio of nuclear luminosities indicates that most of the systems are likely in the coalescence stage of a major merger. This is supported by the existence of a single galaxy disk in 65% of the systems studied and the finding of a correlation between nuclear luminosity and host luminosity for the single-disk systems: those sources fitted with as single disk are in a more evolved stage of the merger and present an enhancement of the nuclear luminosity compared to the double-disk systems, as expected from simulations of galaxy mergers. Finally, we identify a sample of 19 double nucleus disk galaxies in which the two nuclei are physically separated by $\leq1$ kpc and constitute thus a sample of sub-kpc binary AGN candidates.

Dark sectors, consisting of new, light, weakly-coupled particles that do not interact with the known strong, weak, or electromagnetic forces, are a particularly compelling possibility for new physics. Nature may contain numerous dark sectors, each with their own beautiful structure, distinct particles, and forces. This review summarizes the physics motivation for dark sectors and the exciting opportunities for experimental exploration. It is the summary of the Intensity Frontier subgroup "New, Light, Weakly-coupled Particles" of the Community Summer Study 2013 (Snowmass). We discuss axions, which solve the strong CP problem and are an excellent dark matter candidate, and their generalization to axion-like particles. We also review dark photons and other dark-sector particles, including sub-GeV dark matter, which are theoretically natural, provide for dark matter candidates or new dark matter interactions, and could resolve outstanding puzzles in particle and astro-particle physics. In many cases, the exploration of dark sectors can proceed with existing facilities and comparatively modest experiments. A rich, diverse, and low-cost experimental program has been identified that has the potential for one or more game-changing discoveries. These physics opportunities should be vigorously pursued in the US and elsewhere.

We present very long baseline interferometry (VLBI) observations of the milliarcsecond-scale radio structure of the supernova remnant SNR 4449$-$1 in the galaxy NGC 4449. This young and superluminous remnant was observed at 1.6 GHz ($\lambda = 18$\,cm) with the European VLBI Network. The observations confirm earlier identifications of this object with a supernova remnant (SNR) while revealing a somewhat different morphology compared with the structure reported by Bietenholz et al. from VLBI observations at 1.4 GHz. This difference is discussed here in the context of structural sensitivity of both observations. The 1.6 GHz image yields accurate estimates of the size (0.0422 arcsec $\times$ 0.0285 arcsec and 0.8 $\times$ 0.5 pc) and age ($\sim$55 yr) of SNR 4449$-$1. With a total flux of 6.1 $\pm$ 0.6 mJy measured in the VLBI image, the historical lightcurve of the source can be well represented by a power-law decay with a power index of $-$1.19 $\pm$ 0.07. The SNR exhibits a decline rate of the radio emission of 2.2$%$ $\pm$ 0.1$%$ yr$^{-1}$ and a radio luminosity of 1.74 $\times$ 10$^{35}$ erg s$^{-1}$.

The most extreme ultraluminous X-ray sources (ULXs), with $L_\mathrm{X} > 5 \times 10^{40}$ erg s$^{-1}$, are amongst the best candidates for hosting intermediate-mass black holes (IMBHs) in the haloes of galaxies. Jet radio emission is expected from a sub-Eddington accreting IMBH in the low/hard (radio bright) state. In a search for such IMBH jet radio emission, we have observed with the Very Large Array (VLA) at 5 GHz a sample of seven extreme ULXs whose X-ray properties indicate they are in the hard state. Assuming they remain in this state, the non-detection of radio emission for six of the target sources allows us to constrain their black hole mass to the IMBH regime, thus ruling out a supermassive black hole nature. For the extreme ULX in the galaxy NGC\,2276, we detect extended radio emission formed by two lobes of total flux density 1.43 $\pm$ 0.22 mJy and size $\sim$650 pc. The X-ray counterpart is located between the two lobes, suggesting the presence of a black hole with jet radio emission. The radio luminosity allows us to constrain the black hole mass of this source to the IMBH regime; hence, the extreme ULX in NGC\,2276 could be the first detection of extended jet radio emission from an IMBH. The radio emission could also possibly come from a radio nebula powered by the ULX with a minimum total energy of $5.9 \times 10^{52}$ erg, thus constituting the most powerful and largest ULX radio nebula ever observed.

The origin of the high X-ray luminosities of most ultraluminous X-ray sources (ULXs) still remains poorly understood. Most of the scenarios proposed to explain their nature assume that ULXs are powered by accretion on to a black hole (BH). In this case, the detection of compact radio emission and the location of the ULXs in the Fundamental Plane (X-ray versus radio-luminosity plane) can provide an estimate of the ULX BH mass and help address the question of their physical nature. We present the results of a high-resolution (very long baseline interferometry) radio observational campaign aimed at detecting and studying compact radio emission in four ULXs with known radio counterparts. We find that one of the targets (N4559-X4) was previously misclassified: its low X-ray luminosity indicates that the source is not a ULX. No milliarcsec-scale radio emission is detected for N4559-X4 nor for the ULXs N4490-X1 and N5194-X2, for which upper limits on the radio luminosities are estimated. These limits argue strongly against the presence of supermassive BHs in these three systems. For N4559-X4, the low X-ray luminosity and the ratio of the radio and X-ray luminosities suggest the presence of an X-ray binary. Compact radio emission is detected for the ULX N5457-X9 within its \textit{Chandra} positional error, making N5457-X9 a potential intermediate-mass BH with steady jet emission.

Measurements in the radio regime embrace a number of effective approaches for WISP searches, often covering unique or highly complementary ranges of the parameter space compared to those explored in other research domains. These measurements can be used to search for electromagnetic tracers of the hidden photon and axion oscillations, extending down to ~10^-19 eV the range of the hidden photon mass probed, and closing the last gaps in the strongly favoured 1-5 micro-eV range for axion dark matter. This provides a strong impetus for several new initiatives in the field, including the WISP Dark Matter eXperiment (WISPDMX) and novel conceptual approaches for broad-band WISP searches in the 0.1-1000 micro-eV range.

The temporal and spatial spectral evolution of the jets of AGN can be studied with multi-frequency, multi-epoch VLBI observations. The combination of both, morphological and spectral parameters can be used to derive source intrinsic physical properties such as the magnetic field and the non-thermal particle density. In the first two papers of this series, we analyzed the single-dish light curves and the VLBI kinematics of the blazar CTA 102 and suggested a shock-shock interaction between a traveling and a standing shock wave as a possible scenario to explain the observed evolution of the component associated to the 2006 flare. In this paper we investigate the core-shift and spectral evolution to test our hypothesis of a shock-shock interaction. We used 8 multi-frequency VLBA observations to analyze the temporal and spatial evolution of the spectral parameters during the flare. We observed CTA 102 between May 2005 and April 2007 using the VLBA at six different frequencies spanning from 2 - 86 GHz. After the calibrated VLBA images were corrected for opacity, we performed a detailed spectral analysis. From the derived values we estimated the magnetic field and the density of the relativistic particles. The detailed analysis of the opacity shift reveals that the position of the jet core is proportional to nu^-1 with some temporal variations. The value suggests possible equipartition between magnetic field energy and particle kinetic energy densities at the most compact regions. From the variation of the physical parameters we deduced that the 2006 flare in CTA 102 is connected to the ejection of a new traveling feature (t=2005.9) and the interaction between this shock wave and a stationary structure around 0.1 mas from the core. The source kinematics together with the spectral and structural variations can be described by helical motions in an over-pressured jet.

Common extensions of the Standard Model of particle physics predict the existence of a "hidden" sector that comprises particles with a vanishing or very weak coupling to particles of the Standard Model (visible sector). For very light (m < 10^-14 eV) hidden U(1) gauge bosons (hidden photons), broad-band radio spectra of compact radio sources could be modified due to weak kinetic mixing with radio photons. Here, search methods are developed and their sensitivity discussed, with specific emphasis on the effect of the coherence length of the signal, instrumental bandwidth, and spectral resolution. We conclude that radio observations in the frequency range of 0.03--1400 GHz probe kinetic mixing of ~10^-3 of hidden photons with masses down to ~10^-17 eV. Prospects for improving the sensitivity with future radio astronomical facilities as well as by stacking data from multiple objects are discussed.

The 2006 radio flare in the blazar CTA102 is the largest ever reported for this source. The analysis of the single--dish light curves revealed a possible shock-shock interaction as a scenario for the 2006 outburst. In order to confirm this hypothesis we analyzed eight multi-frequency (2 GHz - 86 GHz) VLBI observations covering the 2006 flare. We detected a radio component ejected around 2005.9 and a stationary feature around 0.1 mas away from the core. We also performed a core-shift and spectral analysis on the fully calibrated VLBI maps. The frequency dependent position of the core slightly deviates from the typical $\nu^{-1}$-dependence during the flaring state. This behaviour may reflect a mismatch between magnetic energy density and energy density of the relativistic particles. From the core-shift corrected maps we obtained the spectral parameters and computed the temporal and spatial evolution of the jet intrinsic parameters such as the magnetic field and the particle density. The overall picture of the source could be best interpreted as an over-pressured jet perturbed by traveling shock waves.

The Square Kilometre Array (SKA) is the most ambitious radio telescope ever planned. With a collecting area of about a square kilometre, the SKA will be far superior in sensitivity and observing speed to all current radio facilities. The scientific capability promised by the SKA and its technological challenges provide an ideal base for interdisciplinary research, technology transfer, and collaboration between universities, research centres and industry. The SKA in the radio regime and the European Extreme Large Telescope (E-ELT) in the optical band are on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI) and have been recognised as the essential facilities for European research in astronomy. This "White Paper" outlines the German science and R&D interests in the SKA project and will provide the basis for future funding applications to secure German involvement in the Square Kilometre Array.

We report the detection of a statistically significant flare-like event in the Mg II\lambda 2800\AA\ emission line of 3C 454.3 during the outburst of autumn 2010. The highest levels of emission line flux recorded over the monitoring period (2008 - 2011) coincide with a superluminal jet component traversing through the radio core. This finding crucially links the broad-emission line fluctuations to the non-thermal continuum emission produced by relativistically moving material in the jet and hence to the presence of broad-line region clouds surrounding the radio core. If the radio core were located at several parsecs from the central black hole then our results would suggest the presence of broad-line region material outside the inner parsec where the canonical broad-line region is envisaged to be located. We briefly discuss the implications of broad-emission line material ionized by non-thermal continuum on the context of virial black hole mass estimates and gamma-ray production mechanisms.

We make use of VLBI observations of the radio jet in the quasar S5 0836+710 at different frequencies and epochs to study its properties. The jet shows helical structure at all frequencies. The ridge-line of the emission in the jet coincides at all frequencies and epochs, within the errors. We conclude that the helicity is a real, physical structure. Small differences between epochs reveal wave-like motion of the ridge-line transversal to the jet propagation axis. These transversal motions are measured to be superluminal. This unphysical result could correspond to a possible small amplitude oscillation of the ridge-line at the radio-core and to large errors in the determination of the positions. In addition, higher resolution images at 15 GHz show that the ridge-line does not coincide exactly with the centre of the radio jet. At arc-second scales, this powerful jet shows non-collimated, irregular structure and a lack of a hot-spot. Following this collection of evidence, we conclude that the ridge-line could be related to a pressure maximum within the jet cross-section, associated with the observed helical pattern that could lead to jet disruption at longer scales. We discuss the implications regarding the nature of the instability and the jet flow itself.

We present results from our analysis of Chandra X-ray Observatory, W. M. Keck Observatory, and Karl G. Jansky Very Large Array (VLA) images of the Crab Nebula that were contemporaneous with the gamma-ray flare of 2011 April. Despite hints in the X-ray data, we find no evidence for statistically significant variations that pinpoint the specific location of the flares within the Nebula. The Keck observations extend this conclusion to the "inner knot", i.e., the feature within an arcsecond of the pulsar. The VLA observations support this conclusion. We also discuss theoretical implications of the gamma-ray flares and suggest that the most dramatic gamma-ray flares are due to radiation-reaction-limited synchrotron emission associated with sudden, dissipative changes in the current system sustained by the central pulsar.

We performed multifrequency multiepoch Very Long Baseline Array (VLBA) observations of the blazar CTA 102 during its 2006 radio flare, the strongest ever reported for this source. These observations provide an excellent opportunity to investigate the evolution of the physical properties of blazars, especially during these flaring events. We want to study the kinematic changes in the source during the strong radio outburst in April 2006 and test the assumption of a shock-shock interaction. This assumption is based on the analysis and modeling of the single-dish observations of CTA\,102 (Paper I). In this paper we study the kinematics of CTA 102 at several frequencies using VLBI observations. From the modeled jet features we derived estimates for the evolution of the physical parameters, such as the particle density and the magnetic field. Furthermore, we combined our observations during the 2006 flare with long-term VLBA monitoring of the source at 15 GHz and 43 GHz. We cross-identified seven features throughout our entire multifrequency observations and find evidence of two possible recollimation shocks around 0.1 mas (deprojected 18 pc at a viewing angle of 2.6 degrees) and 6.0 mas (deprojected 1 kpc) from the core. The 43 GHz observations reveal a feature ejected at epoch $t_\mathrm{ej}=2005.9\pm0.2$, which could be connected to the 2006 April radio flare. Furthermore, this feature might be associated with the traveling component involved in the possible shock-shock interaction, which gives rise to the observed double peak structure in the single-dish light curves reported in Paper I.

The remarkable stability of extragalactic jets is surprising, given the reasonable possibility of the growth of instabilities. In addition, much work in the literature has invoked this possibility in order to explain observed jet structures and obtain information from these structures. For example, it was recently shown that the observed helical structures in the jet in S5 0836+710 could be associated with helical pressure waves generated by Kelvin-Helmholtz instability. Our aim is to resolve the arc-second structure of the jet in the quasar S5 0836+710 and confirm the lack of a hot-spot (reverse jet-shock) found by present observing arrays, as this lack implies a loss of jet collimation before interaction with the intergalactic medium. In this work, we use an observation performed in 2008 using EVN and MERLIN. The combined data reduction has provided a complete image of the object at arc-second scales. The lack of a hot-spot in the arc-second radio structure is taken as evidence that the jet losses its collimation between the VLBI region and the region of interaction with the ambient medium. This result, together with the previous identification of the helical structures in the jet with helical pressure waves that grow in amplitude with distance, allow us to conclude that the jet is probably disrupted by the growth of Kelvin-Helmholtz instability. This observational evidence confirms that the physical parameters of jets can be extracted using the assumption that instability is present in jets and can be the reason for many observed structures. Interestingly, the observed jet is classified as a FRII object in terms of its luminosity, but its large-scale morphology does not correspond to this classification. The implications of this fact are discussed.

We have investigated a frequency-dependent shift in the absolute position of the optically thick apparent origin of parsec-scale jets ("core shift" effect) to probe physical conditions in ultra-compact relativistic outflows in AGN. We used multi-frequency Very Long Baseline Array (VLBA) observations of 191 sources carried out in 12 epochs in 2006 within the MOJAVE program. The observations were performed at 8.1, 8.4, 12.1, and 15.4 GHz. We implemented a method of determining the core shift vector based on (i) image registration by two-dimensional normalized cross-correlation and (ii) model-fitting the source brightness distribution to take into account a non-zero core component offset from the phase center. The 15.4-8.1, 15.4-8.4, and 15.4-12.1 GHz core shift vectors are derived for 163 sources, and have median values of 0.128, 0.125, and 0.088 mas, respectively, compared to the typical measured errors of 0.050, 0.051, 0.035 mas. The effect occurs predominantly along the jet direction, with departures smaller than 45 deg from the median jet position angle in over 80% of the cases. Despite the moderate ratio of the observed frequencies (<2), core shifts significantly different from zero (>2sigma) are detected for about 55% of the sources. These shifts are even better aligned with the jet direction, deviating from the latter by less than 30 deg in over 90% of the cases. There is an indication that the core shift decreases with increasing redshift. Magnetic fields in the jet at a distance of 1 parsec from the central black hole, calculated from the obtained core shifts, are found to be systematically stronger in quasars (median B1~0.9 G) than those in BL Lacs (median B1~0.4 G). We also constrained the absolute distance of the core from the apex of the jet at 15 GHz as well as the magnetic field strength in the 15 GHz core region.

We present new spectroscopic identifications of 12 X-shaped radio galaxies and use the spectral data to derive starburst histories and masses of the nuclear supermassive black holes in these galaxies. The observations were done with the 2.1-m telescope of the Observatorio Astron\'omico Nacional at San Pedro M\'artir, M\'exico. The new spectroscopic results extend the sample of X-shaped radio galaxies studied with optical spectroscopy. We show that the combined sample of the X-shaped radio galaxies has statistically higher black-hole masses and older episodes of star formation than a control sample of canonical double-lobed radio sources with similar redshifts and luminosities. The data reveal enhanced star-formation activity in the X-shaped sample on the timescales expected in galactic mergers. We discuss the results obtained in the framework of the merger scenario.

We have used the European VLBI Network (EVN) to observe a sample of Lyman Break Analogs (LBAs), nearby (z < 0.3) galaxies with properties similar to the more distant Lyman Break Galaxies (LBGs). The study of LBGs may help define the feed-back relationship between black holes (BHs) and their host galaxies. Previous VLA observations have shown that the kpc-scale radio emission from LBAs is dominated by starbursts. The main targets of this VLBI experiment were selected because they possessed emission-line properties between starbursts and Type 2 (obscured) AGN. Eight targets (three star-forming LBAs, four composite LBAs, and one Type 1 AGN) were observed at 5 GHz, four of which were also observed at 1.7 GHz. One star-forming LBA and one composite LBA were detected above 5 \sigma at 1.7 GHz (only), while the AGN was detected at 5 GHz. In both LBAs, the radio luminosity (LR) exceeds that expected from supernovae (remnants) based on a comparison with Arp220, Arp229A and Mrk273, by factors of 2 - 8. The composite LBA exhibits a compact core emitting around 10% of the VLA flux density. The high Tb of 3.5E7 K and excess core L_R with respect to the L_R/L_X relation of radio-quiet AGN indicate that this LBA possesses an obscured AGN (MBH ~ E5-E7 M_sun). While weak AGN may co-exist with the starbursts as shown in at least one of the LBAs, their contribution to the total radio flux is fairly minimal. Our results show that the detection of such weak AGN presents a challenge at radio, X-ray and optical emission-line wavelengths at z ~ 0.2, indicating the great difficulties that need to be overcome in order to study similar processes at high redshift when these types of galaxies were common.

We review the effects of source size in interferometric observations and focus on the cases of very compact sources. If a source is extremely compact and/or weak (so it is not possible to detect signature of source structure in the visibilities) we describe a test of hypothesis that can be used to set a strong upper limit to the size of the source. We also estimate the minimum possible size of a source whose structure can still be detected by an interferometer (i.e., the maximum theoretical over-resolution power of an interferometer), which depends on the overall observing time, the compactness in the array distribution, and the sensitivity of the receivers. As a result, and depending on the observing frequency, the over-resolution power of forthcoming ultra-sensitive arrays, like the Square Kilometer Array (SKA), may allow us to study details of sources at angular scales down to a few micro-as.

Helical structures are common in extragalactic jets. They are usually attributed in the literature to periodical phenomena in the source (e.g., precession). In this work, we use VLBI data of the radio-jet in the quasar S5 0836+710 and hypothesize that the ridge-line of helical jets like this corresponds to a pressure maximum in the jet and assume that the helically twisted pressure maximum is the result of a helical wave pattern. For our study, we use observations of the jet in S5 0836+710 at different frequencies and epochs. The results show that the structures observed are physical and not generated artificially by the observing arrays. Our hypothesis that the observed intensity ridge-line can correspond to a helically twisted pressure maximum is confirmed by our observational tests. This interpretation allows us to explain jet misalignment between parsec and kiloparsec scales when the viewing angle is small, and also brings us to the conclusion that high-frequency observations may show only a small region of the jet flow concentrated around the maximum pressure ridge-line observed at low frequencies. Our work provides a potential explanation for the apparent transversal superluminal speeds observed in several extragalactic jets by means of transversal shift of an apparent core position with time.

The analysis of the radio light curves of the blazar CTA102 during its 2006 flare revealed a possible interaction between a standing shock wave and a traveling one. In order to better understand this highly non-linear process, we used a relativistic hydrodynamic code to simulate the high energy interaction and its related emission. The calculated synchrotron emission from these simulations showed an increase in turnover flux density, $S_{m}$, and turnover frequency, $\nu_{m}$, during the interaction and decrease to its initial values after the passage of the traveling shock wave.

Aims: 3C345 was recently identified as a gamma-ray emitter, based on the first 20 months of Fermi-LAT data and optical monitoring. In this paper, a connection between the gamma-ray and optical variability of 3C345 and properties of its parsec-scale radio emission is investigated. Methods: The Fermi-LAT data of 3C345, covering an energy range of 0.1-300 GeV, were combined with 32 Very Long Baseline Array observations of the object made at 43.2 GHz in the period of January 2008 - March 2010. Results: The VLBA data reveal morphology and kinematics of the flow on scales of up to ~5 milliarcseconds (mas; deprojected linear distances of 380 parsecs). The brightness temperature, T_b(r), measured along the jet first decreases with distance proportional to r^-(0.95 +/-0.69) and later exhibits a break at ~0.3 mas, with T_b(r) proportional to r^-(4.11 +/-0.85) at larger separations. Variations of the gamma-ray, optical and parsec-scale radio emission show a similar long-term trend persistent during the entire VLBA monitoring period. The gamma-ray and optical variations on shorter time scales are related to structural changes in the jet on scales of ~0.3 mas (~23 parsecs, deprojected), with the gamma-ray and optical flares possibly related to the evolution of four distinct superluminal components identified in the flow. Conclusions: The observations indicate that both the quiescent and flaring components of the gamma-ray emission are produced in a region of the jet of ~23 pc in extent. This region may mark the Compton-loss dominated zone of the flow and its large extent may favor the synchrotron self-Compton mechanism for gamma-ray production in the relativistic jet of the quasar 3C345.

The effect of a frequency dependent shift of the VLBI core position (known as the "core shift") was predicted more than three decades ago and has since been observed in a few sources, but often within a narrow frequency range. This effect has important astrophysical and astrometric applications. To achieve a broader understanding of the core shift effect and the physics behind it, we conducted a dedicated survey with NRAO's Very Long Baseline Array (VLBA). We used the VLBA to image 20 pre-selected sources simultaneously at nine frequencies in the 1.4-15.4 GHz range. The core position at each frequency was measured by referencing it to a bright, optically thin feature in the jet. A significant core shift has been successfully measured in each of the twenty sources observed. The median value of the core shift is found to be 1.21 mas if measured between 1.4 and 15.4 GHz, and 0.24 mas between 5.0 and 15.4 GHz. The core position, r, as a function of frequency, n, is found to be consistent with an r n^-1 law. This behavior is predicted by the Blandford & Koenigl model of a purely synchrotron self-absorbed conical jet in equipartition. No systematic deviation from unity of the power law index in the r(n) relation has been convincingly detected. We conclude that neither free-free absorption nor gradients in pressure and/or density in the jet itself and in the ambient medium surrounding the jet play a significant role in the sources observed within the 1.4-15.4 GHz frequency range. These results support the interpretation of the parsec-scale core as a continuous Blandford-Koenigl type jet with smooth gradients of physical properties along it.

For more than 15 years, since the days of the Energetic Gamma-Ray Experiment Telescope (EGRET) on board the Compton Gamma-Ray Observatory (CGRO; 1991-2000), it has remained an open question why the prominent blazar 3C 345 was not reliably detected at gamma-ray energies <=20 MeV. Recently a bright gamma-ray source (0FGL J1641.4+3939/1FGL J1642.5+3947), potentially associated with 3C 345, was detected by the Large Area Telescope (LAT) on Fermi. Multiwavelength observations from radio bands to X-rays (mainly GASP-WEBT and Swift) of possible counterparts (3C 345, NRAO 512, B3 1640+396) were combined with 20 months of Fermi-LAT monitoring data (August 2008 - April 2010) to associate and identify the dominating gamma-ray emitting counterpart of 1FGL J1642.5+3947. The source 3C 345 is identified as the main contributor for this gamma-ray emitting region. However, after November 2009 (15 months), a significant excess of photons from the nearby quasar NRAO 512 started to contribute and thereafter was detected with increasing gamma-ray activity, possibly adding flux to 1FGL J1642.5+3947. For the same time period and during the summer of 2010, an increase of radio, optical and X-ray activity of NRAO 512 was observed. No gamma-ray emission from B3 1640+396 was detected.

"High-resolution", or "long-baseline", science with the SKA and its precursors covers a broad range of topics in astrophysics. In several research areas, the coupling between improved brightness sensitivity of the SKA and a sub-arcsecond resolution would uncover truly unique avenues and opportunities for studying extreme states of matter, vicinity of compact relativistic objects, and complex processes in astrophysical plasmas. At the same time, long baselines would secure excellent positional and astrometric measurements with the SKA and critically enhance SKA image fidelity at all scales. The latter aspect may also have a substantial impact on the survey speed of the SKA, thus affecting several key science projects of the instrument.

We report on very long baseline interferometry (VLBI) observations of the radio emission from the inner region of the Crab Nebula, made at 1.6 GHz and 5 GHz after a recent high-energy flare in this object. The 5 GHz data have provided only upper limits of 0.4 milli-Jansky (mJy) on the flux density of the pulsar and 0.4 mJy/beam on the brightness of the putative flaring region. The 1.6 GHz data have enabled imaging the inner regions of the nebula on scales of up to ~40". The emission from the inner "wisps" is detected for the first time with VLBI observations. A likely radio counterpart (designated "C1") of the putative flaring region observed with Chandra and HST is detected in the radio image, with an estimated flux density of $0.5\pm 0.3$\,mJy and a size of 0.2-0.6". Another compact feature ("C2") is also detected in the VLBI image closer to the pulsar, with an estimated flux density of 0.4 +- 0.2 mJy and a size smaller than 0{\farcs}2. Combined with the broad-band SED of the flare, the radio properties of C1 yield a lower limit of ~0.5 mG for the magnetic field and a total minimum energy of 1.2*10^41 ergs vested in the flare (corresponding to using about 0.2% of the pulsar spin-down power). The 1.6 GHz observations provide upper limits for the brightness (0.2 mJy/beam) and total flux density (0.4 mJy) of the optical Knot 1 located at 0.6" from the pulsar. The absolute position of the Crab pulsar is determined, and an estimate of the pulsar proper motion is obtained.

Context: The blazar CTA 102 (z=1.037) underwent a historical radio outburst in April 2006. This event offered a unique chance to study the physical properties of the jet. Aims: We used multifrequency radio and mm observations to analyze the evolution of the spectral parameters during the flare as a test of the shock-in-jet model under these extreme conditions. Methods: For the analysis of the flare we took into account that the flaring spectrum is superimposed on a quiescent spectrum. We reconstructed the latter from archival data and fitted a synchrotron self-absorbed distribution of emission. The uncertainties of the derived spectral parameters were calculated using Monte Carlo simulations. The spectral evolution is modeled by the shock-in-jet model, and the derived results are discussed in the context of a geometrical model (varying viewing angle) and shock-shock interaction. Results: The evolution of the flare in the turnover frequency-turnover flux density plane shows a double peak structure. The nature of this evolution is dicussed in the frame of shock-in-jet models. We discard the generation of the double peak structure in the turnover frequency-turnover flux density plane purely based on geometrical changes (variation of the Doppler factor). The detailed modeling of the spectral evolution favors a shock-shock interaction as a possible physical mechanism behind the deviations from the standard shock-in-jet model.

We present results from our studies of radio emission from selected Ultraluminous X-ray (ULX) sources, using archival Giant Metrewave Radio Telescope (GMRT) data and new European VLBI Network (EVN) observations. The GMRT data are used to find possible faint radio emission from ULX sources located in late-type galaxies in the Chandra Deep Fields. No detections are found at 235, 325 and 610 MHz, and upper limits on the radio flux densities at these frequencies are given. The EVN observations target milliarcsecond-scale structures in three ULXs with known radio counterparts (N4449- X1, N4088-X1, and N4861-X2). We confirm an earlier identification of the ULX N4449-X1 with a supernova remnant and obtain the most accurate estimates of its size and age. We detect compact radio emission for the ULX N4088-X1, which could harbour an intermediate mass black hole (IMBH) of 10^5 M\odot accreting at a sub-Eddington rate. We detect a compact radio component in the ULX N4861-X2, with a brightness temperature > 10^6 K and an indication for possible extended emission. If the extended structure is confirmed, this ULX could be an HII region with a diameter of 8.6 pc and surface brightness temperature \geq 10^5 K. The compact radio emission may be coming from a ~ 10^5 M\odot black hole accreting at 0.1M_Edd.

It has been suggested that the X-shaped morphology observed in some radio sources can reflect either a recent merger of two supermassive black holes (SMBHs) or the presence of a second active black hole in the galactic nucleus. These scenarios are tested by studying the relationship between the black hole mass, radio and optical luminosity, starburst history, and dynamic age of radio lobes in a sample of 29 X-shaped radio galaxies drawn from a list of 100 X-shaped radio source candidates identified from the FIRST survey. The same relationships are also studied in a control sample consisting of 36 radio-loud active nuclei with similar redshifts and optical and radio luminosities. The X-shaped objects are found to have statistically higher black hole masses and older starburst activity compared to the objects from the control sample. Implications of these findings are discussed for the black hole merger scenario and for the potential presence of active secondary black holes in post-merger galaxies.

Results of processing of data of a VLBI experiment titled RAPL01 are presented. These VLBI observations were made on 4th February, 2010 at 6.28 cm between the 100-m antenna of the Max Planck Institute (Effelsberg, Germany), Puschino 22-m antenna (Astro Space Center (ASC), Russia), and two 32-m antennas of the Istituto di Radioastronomia di Bologna (Bologna, Italy) in Noto and Medicina. 2 well-known sources, 3C84 (0316+413), and BL Lac (2200+420) were included in the schedule of observations. Each of them was observed during 1 hour at all the stations. The Mark-5A registration system was used at 3 European antennae. The alternative registration system known as RDR (RADIOASTRON Data Recorder) was used in Puschino. The Puschino data were recorded in format RDF (RADIOASTRON Data Format). Two standard recording modes designed as 128-4-1 (one bit), and 256-4-2 (two bit) were used in the experiment. All the Mark-5A data from European antennae were successfully converted into the RDF format. Then, the correlation function was estimated at the ASC software correlator. A similar correlation function also was estimated at the Bonn correlator. The Bonn correlator reads Mark5A data, the RDF format was converted into Mark5B format before correlation. The goal of the experiment was to check the functioning and data analysis of the ground based radio telescopes for the RADIOASTRON SVLBI mission

3C 345 is one of the archetypical active galactic nuclei, showing structural and flux variability on parsec scales near a compact unresolved radio core. During the last 2 years, the source has been undergoing a period of high activity visible in the broad spectral range, from radio through high-energy bands. We have been monitoring parsec-scale radio emission in 3C 345 during this period at monthly intervals, using the VLBA at 15, 24, and 43 GHz. Our radio observations are compared with gamma-ray emission detected by Fermi-LAT in the region including 3C 345 (1FGL J1642.5+3947). Three distinct gamma-ray events observed in this region are associated with the propagation of relativistic plasma condensations inside the radio jet of 3C 345. We report on evidence for the gamma-rays to be produced in a region of the jet of up to 40 pc (de-projected) in extent. This suggests the synchrotron self-Compton process as the most likely mechanism for production of gamma-rays in the source.

The radio light curve and spectral evolution of the blazar CTA 102 during its 2006 outburst can be rather well explained by the standard shock-in-jet model. The results of a pixel-to-pixel spectral analysis of multi-frequency VLBI images, together with kinematics derived from the MOJAVE survey lead to the picture of an over-pressured jet with respect to the ambient medium. The interaction of a traveling shock wave with a standing one is a possible scenario which could explain the observed spectral behaviour

The blazar CTA 102 underwent a major radio flare in April 2006. We used several 15 GHz VLBI observations from the MOJAVE program to investigate the influence of this extreme event on jet kinematics. The result of modeling and analysis lead to the suggestion of an interaction between traveling and standing shocks 0.2 mas away from the VLBI core.

Apparent evolution of relativistic flows as traced by radio emission results from a combination of several factors related to propagation of relativistic blobs or shocks, velocity, density and pressure stratification of the underlying flow, plasma instability and (possibly also) phase and time travel effect. This combination can create an intricate and chaotic patterns of the observed morphological changes in radio emission, which complicates the analysis and interpretation of kinematic and physical properties of the jet plasma. Recent studies have indicated that slow and quasi-stationary patterns in jets are most likely formed by plasma instabilities while faster, superluminally moving patterns are related to highly relativistic plasma condensations produced by the nuclear flares. Some of the stationary patterns may also be related to recollimation shocks or locations where strong non-thermal continuum is produced in jets. Similarities and differences of the AGN and XRB jets in this respect are reviewed.

(abridged) Relativistic jets, formed in the vicinity of central supermassive black holes in AGN, show ample evidence connecting them to physical conditions in the accretion disc and broad-line region. The jets are responsible for a large fraction of non-thermal continuum emission (particularly during powerful flares), which makes understanding their physics an important aspect of studies of blazars characterised by profound flaring activity arising from extremely compact regions. Imaging and polarimetry of radio emission on milliarcsecond scales provided by very long baseline interferometry (VLBI) offers a range of possibilities for studying ultra-compact regions in relativistic jets and relating them to main manifestations of the blazar activity in AGN. Simultaneous monitoring of optical/high energy variability and evolution of parsec-scale radio structures yields arguably the most detailed picture of the relation between acceleration and propagation of relativistic flows and non-thermal continuum generation in blazars. These effects are reviewed and discussed in the context of deriving accurate and self-consistent models for central regions of blazars.

Single-zone synchrotron self-Compton and external Compton models are widely used to explain broad-band Spectral Energy Distributions (SEDs) of blazars from infrared to gamma-rays. These models bear obvious similarities to the homogeneous synchrotron cloud model which is often applied to explain radio emission from individual components of parsec-scale radio jets. The parsec-scale core, typically the brightest and most compact feature of blazar radio jet, could be the source of high-energy emission. We report on ongoing work to test this hypothesis by deriving the physical properties of parsec-scale radio emitting regions of twenty bright Fermi blazars using dedicated 5-43 GHz VLBA observations and comparing these parameters to results of SED modeling.

The mechanism and the region of generation of variable continuum emission are poorly understood for radio-loud AGN because of a complexity of the nuclear region. High-resolution radio VLBI (very long baseline interferometry) observations allow zooming into a subparsec-scale region of the jet in the radio-loud galaxy 3C 390.3. We combined the radio VLBI and the optical data covering the time period of 14 years to look for a link between optical flares and parsec-scale jet. We identify two stationary and nine moving radio features in the innermost subparsec-scale region of the jet. All nine ejections are associated with optical flares. We found a significant correlation (at a confidence level of $>99.99%$) between the ejected jet components and optical continuum flares. The epochs at which the moving knots pass through the location of a stationary radio feature and the optical light curve reaches the maximum are correlated. The radio events follow the maxima of optical flares with the mean time delay of $0.10\pm0.04$ years. This correlation can be understood if the variable optical continuum emission is generated in the innermost part of the jet. A possible mechanism of the energy release is the ejection of knots of high-energy electrons that are accelerated in the jet flow and generate flares of synchrotron continuum emission in the wide range of frequencies from radio to X-ray bands. In this scenario, the beamed optical continuum emission from the jet and counterjet ionizes a gas in a subrelativistic outflow surrounding the jet, which results in a formation of two outflowing conical regions with broad emission lines (in addition to the conventional broad line region around the central nucleus).

We report a relation between radio emission in the inner jet of the Seyfert galaxy 3C 120 and optical continuum emission in this galaxy. Combining the optical variability data with multi-epoch high-resolution very long baseline interferometry observations reveals that an optical flare rises when a superluminal component emerges into the jet and its maxima is related to the passage of such component through the location a stationary feature at a distance of ~1.3 parsecs from the jet origin. This indicates that a significant fraction of the optical continuum produced in 3C 120 is non-thermal and it can ionize material in a sub-relativistic wind or outflow. We discuss implications of this finding for the ionization and structure of the broad emission line region, as well as for the use of broad emission lines for determining black hole masses in radio-loud AGN.

The spatially resolved broad-band spectroscopy with Very Long Baseline Interferometry (VLBI) is one of the few methods that can probe the physical conditions inside blazar jets. We report on measurements of the magnetic field strength in parsec-scale radio structures of selected bright Fermi blazars, based on fitting the synchrotron spectrum to VLBA images made at seven frequencies in a 4.6 -- 43.2 GHz range. Upper limits of B <= 10^-2 -- 10^2 G (observer's frame) could be placed on the magnetic field strength in 13 sources. Hard radio spectra (-0.5 <= a <= +0.1, S_nu ~ nu^a) observed above the synchrotron peak may either be an indication of a hard energy spectrum of the emitting electron population or result from a significant inhomogeneity of the emitting region.

The Square Kilometre Array (SKA) will be operating at the time when several new large optical, X-ray and Gamma-ray facilities are expected to be working. To make SKA both competitive and complementary to these large facilities, thorough design studies are needed, focused in particular on imaging performance of the array. One of the crucial aspects of such studies is the choice of the array configuration, which affects substantially the resolution, rms noise, sidelobe level and dynamic range achievable with the SKA. We present here a quantitative assessment of the effect of the array configuration on imaging performance of the SKA, introducing the spatial dynamic range (SDR) and a measure of incompleteness of the Fourier domain coverage ($\Delta u/u$) as prime figures of merit.

The 16th magnitude quasar 3C 345 (redshift z=0.5928) shows structural and emission variability on parsec scales around a compact unresolved radio core. For the last three decades it has been closely monitored with very long baseline interferometry (VLBI), yielding a wealth of information about the physics of relativistic outflows and dynamics of the central regions in AGN. We present here preliminary results for the long-term jet evolution, based on the 15 GHz monitoring data collected by the MOJAVE survey and various other groups over the last ~14 years and combined with data from earlier VLBI observations of 3C 345 which started in 1979. We discuss the trajectories, kinematics, and flux density evolution of enhanced emission regions embedded in the jet and present evidence for geometrical (e.g. precession) and physical (e.g. relativistic shocks and plasma instability) factors determining the morphology and dynamics of relativistic flows on parsec scales.

The nature of X-shaped sources is a matter of considerable debate: it has even been proposed that they provide evidence for black hole mergers$ / $spin reorientation, and therefore constrain the rate of strong gravitational wave events (Merritt & Ekers 2002). Based on morphological and spectral characteristics of these sources, currently a strong contender to explain the nature of these sources is the `alternative' model of Lal & Rao (2007), in which these sources consist of two pairs of jets, which are associated with two unresolved AGNs. Detailed morphological and spectral results on milliarcsecond-scales (mas) provide a crucial test of this model, and hence these sources are excellent candidates to study on mas; {\it i.e.}, to detect he presence/absence of double nuclei/AGNs, signs of helical/disrupted jets, thereby, to investigate spatially resolved/unresolved binary AGN systems and providing clues to understanding the physics of merging of AGNs on mas. We conducted a systematic study of a large sample of known X-shaped, comparison FR II radio galaxies, and newly discovered X-shaped candidate sources using Giant Metrewave Radio Telescope and Very Large Array at several radio frequencies. In our new observations of `comparison' FR II radio galaxies we find that almost all of our targets show standard spectral steepening as a function of distance from the hotspot. However, one source, 3C 321, has a low-surface-brightness extension that shows a flatter spectral index than the high-surface-brightness hotspots$ / $lobes, as found in `known' X-shaped sources.

In extragalactic jets, the apparent position of the bright/narrow end (the core) depends on the observing frequency, owing to synchrotron self-absorption and external absorption. The effect must be taken into account in order to achieve unbiased results from multi-frequency VLBI data on AGN jets. Multi-frequency core shift measurements supplemented by other data enable estimating the absolute geometry and a number of fundamental physical properties of the jets and their environment. We have previously measured the shift between 13 and 3.6 cm in a sample of 29 AGNs to range between 0 and 1.4 mas. In these proceedings, we present and discuss first results of our follow-up study using VLBA between 1.4 and 15.4 GHz.

Variations of the optical continuum emission in the radio galaxy 3C 390.3 are compared to the properties of radio emission from the compact, sub-parsec-scale jet in this object. We showed that very long-term variations of optical continuum emission (>10 years) is correlated with the radio emission from the base of the jet located above the disk, while the optical long-term variations (1-2 years) follows the radio flares from the stationary component in the jet with time delay of about 1 year. This stationary feature is most likely to be a standing shock formed in the continuous relativistic flow seen at a distance of ~0.4 parsecs from the base of the jet. To account for the correlations observed we propose a model of the nuclear region of 3C 390.3 in which the beamed continuum emission from the jet and counterjet ionizes material in a subrelativistic outflow surrounding the jet. This results in the formation of two conical regions with double-peaked broad emission lines (in addition to the conventional broad line region around the central nucleus) at a distance ~0.6 parsecs from the central engine.

Results of high-resolution simultaneous multi-frequency 8.1-15.4 GHz VLBA polarimetric observations of relativistic jets in active galactic nuclei (the MOJAVE-2 project) are analyzed. We compare characteristics of VLBI features with jet model predictions and test if adiabatic expansion is a dominating mechanism for the evolution of relativistic shocks in parsec-scale AGN jets. We also discuss magnetic field configuration, both predicted by the model and deduced from electric vector position angle measurements.

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

The apparent position of the origin (core) of extragalactic radio jets shifts with the observing frequency, owing to synchrotron self-absorption and external absorption. One of the largest core shifts was detected by us in the quasar 0850+581 between 2 and 8 GHz. We have followed this up recently by a dedicated VLBA experiment at 5, 8, 15, 24, and 43 GHz. First results from this study enabled estimating the absolute geometry and physical conditions in the parsec-scale jet origin.

Context. The nature of X-shaped sources is a matter of considerable debate in the literature: it has even been proposed that they provide evidence for black-hole-mergers/spin-reorientation, and therefore constrain the rate of strong gravitational wave events. Aim. To explore the nature of these X-shaped radio galaxies. Method. We conduct a systematic study of a large sample of known and newly discovered X-shaped sources along with a comparison sample. We used the Giant Metrewave Radio Telescope with resolution of $\sim6^{\prime\prime}$ to $\sim15^{\prime\prime}$ at 610 MHz and 240 MHz in the dual-frequency mode. Preliminary Result. Based on our careful analysis and estimation of the possible systematic errors, the known X-shaped sources divide into the following three categories: (i) the wings have flatter spectral indices than the active lobes, (ii) the wings and the active lobes have comparable spectral indices, and (iii) the wings have steeper spectral indices than the active lobes. In addition, based on our preliminary analysis, one sample source from our comparison sample shows a spectral index result belonging to category (i). Future. Milliarcsecond scale imaging will be conducted on some of these X-shaped sources to investigate if they are examples of binary AGN systems and thereby understand the nature of these sources.

The Square Kilometre Array is going to become operational at the time when several new large optical, X-ray and Gamma-ray telescopes are expected to be working. The main drive for building the SKA is a significant improvement of sensitivity that would widen the general scope of the centimetre-wavelength radio science. To this end, a thorough design studies should be made, in order to ensure that the SKA becomes a competitive counterpart of the facilities at other wavebands. To quantify imaging performance of the SKA configurations, we are implementing figures of merit based on spatial dynamic range (SDR). This work is focused on extensive numerical tests of the analytical, SDR--based figures of merit derived in the SKA Memo 38 (A. Lobanov). Here, we present our preliminary results.

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.

We present results of simultaneous dual-frequency (2 GHz and 8 GHz) very long baseline interferometry (VLBI) observations of 12 active galactic nuclei with prominent jets. Spectral properties of the jets and evolution of their brightness temperature are discussed. Measured sizes and brightness temperatures of VLBI features are found to be consistent with emission from relativistic shocks dominated by adiabatic losses. Physical scenarios with different magnetic field orientation in the jets are discussed.

The apparent position of the "core" in a parsec-scale radio jet (a compact, bright emitting region at the narrow end of the jet) depends on the observing frequency, owing to synchrotron self-absorption and external absorption. While providing a tool probing physical conditions in the vicinity of the core, this dependency poses problems for astrometric studies using compact radio sources. We investigated the frequency-dependent shift in the positions of the cores (core shift) observed with very long baseline interferometry (VLBI) in parsec-scale jets. We discuss related physics, as well as its effect on radio astrometry and the connection between radio and optical positions of astrometric reference objects. We searched for the core shift in a sample of 277 radio sources imaged at 2.3 GHz (13 cm) and 8.6 GHz (4 cm) frequency bands using VLBI observations made in 2002 and 2003. The core shift was measured by referencing the core position to optically thin jet features whose positions are not expected to change with frequency. We present here results for 29 selected active galactic nuclei (AGN) with bright distinct VLBI jet features that can be used in differential measurements and that allow robust measurements of the shift to be made. In these AGN, the magnitude of the measured core shift between 2.3 and 8.6 GHz reaches 1.4 mas, with a median value for the sample of 0.44 mas. Nuclear flares result in temporal variability of the shift. An average shift between the radio (4 cm) and optical (6000 Angstrom) bands is estimated to be approximately 0.1 mas, and it should be taken into account in order to provide the required accuracy of the radio-optical reference frame connection. This can be accomplished with multi-frequency VLBI measurements... (abridged)

We present results from a large 86GHz global VLBI survey of compact radio sources. The main goal of the survey is to increase by factors of 3--5 the total number of objects accessible for future 3-mm VLBI imaging. The survey observations reach the baseline sensitivity of 0.1Jy and image sensitivity of better than 10 mJy/beam. The total of 127 compact radio sources have been observed. The observations have yielded images for 109 sources, extending the database of the sources imaged at 86GHz with VLBI observation by a factor of 5, and only 6 sources have not been detected. The remaining 12 objects have been detected but could not be imaged due to insufficient closure phase information. Radio galaxies are less compact than quasars and BL Lacs on sub-milliarcsecond scale. Flux densities and sizes of core and jet components of all imaged sources have been estimated using Gaussian model fitting. From these measurements, brightness temperatures have been calculated, taking into account resolution limits of the data. The cores of 70% of the imaged sources are resolved. The core brightness temperatures of the sources peak at $\sim 10^{11}$ K and only 1% have brightness temperatures higher than $10^{12}$ K. Cores of Intraday Variable (IDV) sources are smaller in angular size than non-IDV sources, and so yield higher brightness temperatures.

We present the first results on the analysis of the structures observed in the jet of the quasar 0836+710. We obtain the ridge lines of the jet at different epochs and several frequencies. We interpret the oscillatory structures obtained as waves that can be attached to the growth of instabilities. We explain how to derive information on the nature and origin of these structures by fitting together the ridge lines at different epochs and frequencies. Finally we show the predictive power of this approach: by generating an artificial wave and applying the corresponding relativistic and projection effects we show that apparent changes in the jet direction in the inner regions of jets can be attached to the transversal motion of structures.

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 aim to test the model proposed to explain the correlation between the flux density at 15 GHz of a stationary component in the parsec-scale jet and the optical continuum emission in the radio galaxy 3C~390.3. In the model, the double-peaked emission from 3C~390.3 is likely to be generated both near the disk and in a rotating subrelativistic outflow surrounding the jet, due to ionization of the outflow by the beamed continuum emission from the jet. This scenario is chosen since broad-emission lines are observed to vary following changes in the inner radio jet. For recent epochs we have imaged and modelled the radio emission of the inner jet of 3C~390.3, which was observed with very long baseline interferometry at 15 GHz, 22 GHz and 43 GHz, to image the inner part of the parsec-scale jet, locate the exact region where the bulk of the continuum luminosity is generated and search for the mechanism that drives the double-peaked profile emission. We present the preliminary results of testing the model using data from 11 years of active monitoring of 3C~390.3.

Accurate alignment of the radio and optical celestial reference frames requires detailed understanding of physical factors that may cause offsets between the positions of the same object measured in different spectral bands. Opacity in compact extragalactic jets (due to synchrotron self-absorption and external free-free absorption) is one of the key physical phenomena producing such an offset, and this effect is well-known in radio astronomy ("core shift"). We have measured the core shifts in a sample of 29 bright compact extragalactic radio sources observed using very long baseline interferometry (VLBI) at 2.3 and 8.6 GHz. We report the results of these measurements and estimate that the average shift between radio and optical positions of distant quasars would be of the order of 0.1-0.2 mas. This shift exceeds positional accuracy of GAIA and SIM. We suggest two possible approaches to carefully investigate and correct for this effect in order to align accurately the radio and optical positions. Both approaches involve determining a Primary Reference Sample of objects to be used for tying the radio and optical reference frames together.

We present results from a large global VLBI(Very Long Baseline Interferometry) survey of compact radio sources at 86 GHz which started in October 2001. The main goal of the survey is to increase the total number of objects accessible for future 3mm-VLBI imaging by factors of 3-5. The survey data reach the baseline sensitivity of 0.1 Jy, and image sensitivity of better than 10 mJy/beam. To date, a total of 127 compact radio sources have been observed. The observations have yielded images for 109 sources, and only 6 sources have not been detected. Flux densities and sizes of core and jet components of all detected sources have been measured using Gaussian model fitting. From these measurements, brightness temperatures have been estimated, taking into account resolution limits of the data. Here, we compare the brightness temperatures of the cores and secondary jet components with similar estimates obtained from surveys at longer wavelengths (e.g. 15 GHz). This approach can be used to study questions related to mechanisms of initial jet acceleration (accelerating or decelerating sub-pc jets?) and jet composition (electron-positron or electron-proton plasma?).

VLBI observations of relativistic outflows (jets) in galactic nuclei, complemented with detailed studies made in other spectral domains, have become an effective tool for investigating the physics of nuclear regions in galaxies. High-resolution radio observations access directly the regions where the jets are formed, and trace their evolution and interaction with the nuclear environment. The emission properties, dynamics, and evolution of jets in AGN are intimately connected to the characteristics of the central supermassive black hole, accretion disk, and broad-line region (BLR) in active galaxies. Large VLBI surveys (15GHz VLBA survey, MOJAVE) and dedicated monitoring programmes follow systematically the evolution of several hundreds of relativistic jets. These observations, combined with optical and X-ray studies, yield arguably the most detailed picture of the galactic nuclei. Recent results from studies of the nuclear regions in several active galaxies with prominent outflows are reviewed in this contribution.

We study the most central regions of AGN jets with an angular resolution of tens of micro-arcseconds using VLBI at millimeter wavelengths (mm-VLBI). We present and discuss a new 86 GHz VLBI survey of compact radio sources. We show new high dynamic range images of two nearby radio galaxies (3C 120 and M87). In M87 the size of the compact VLBI core (the jet base) is < 15 Schwarzschild radii. Future mm-VLBI observations at 1 mm and shorter wavelengths should lead to images of galactic and extragalactic radio sources with a spatial resolution down to a few Schwarzschild radii of the central super massive black holes. To achieve this, the participation of large and sensitive millimeter and sub-millimeter telescopes in VLBI is essential. Owing to their high sensitivity existing and planned phased interferometers like the IRAM Plateau de Bure interferometer, CARMA, the SMA and ALMA could play an important role.

Studying transversal structure in extragalactic jets is crucial for understanding their physics. The Japanese led space VLBI project VSOP has offered arguably the best opportunity for such studies, by reaching baseline lengths of up to 36,000 km and resolving structures down to an angular size of 0.3 mas at 5 GHz. VSOP observations of the jet in 0836+710 at 1.6 and 5 GHz have enabled tracing the radial structure of the flow on scales from 2 mas to 200 mas and determining the wavelengths of individual oscillatory modes responsible for the formation of the structure observed. We conclude that these modes are produced by Kelvin-Helmholtz instability in a sheared relativistic flow. Our results point towards the stratification of the jet and the growth of different modes at different jet radii. We also discuss the implications of the driving frequency on the physics of the active nucleus of the quasar.

Compact relativistic jets in active galactic nuclei offer an effective tool for investigating the physics of nuclear regions in galaxies. The emission properties, dynamics, and evolution of jets in AGN are closely connected to the characteristics of the central supermassive black hole, accretion disk and broad-line region in active galaxies. Recent results from studies of the nuclear regions in several active galaxies with prominent outflows are reviewed in this contribution.

Studying the internal structure of extragalactic jets is crucial for understanding their physics. The Japanese-led space VLBI project VSOP has presented an opportunity for such studies, by reaching baseline lengths of up to 36,000 km and resolving structures down to an angular size of $\approx 0.3$ mas at 5 GHz. VSOP observations of the jet in 0836+710 at 1.6 and 5 GHz have enabled tracing of the radial structure of the flow on scales from 2 mas to 200 mas along the jet and determination of the wavelengths of individual oscillatory modes responsible for the formation of the structure observed. We apply linear stability analysis to identify the oscillatory modes with modes of Kelvin-Helmholtz instability that match the wavelengths of the structures observed. We find that the jet structure in 0836+710 can be reproduced by the helical surface mode and a combination of the helical and elliptic body modes of Kelvin-Helmholtz instability. Our results indicate that the jet is substantially stratified and different modes of the instability grow inside the jet at different distances to the jet axis. The helical surface mode can be driven externally, and we discuss the implications of the driving frequency on the physics of the active nucleus in 0836+710.

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.

In this paper, we present new VLBA observations of the radio jet in the quasar S5 0836+710 at 8 and 22 GHz. The identification of the ridge lines allow us to interpret the jet structure in terms of Kelvin-Helmholtz instability. Combined with previous epochs of VLBA and VSOP data at 1.6 and 5 GHz, these new observations will allow us to study the evolution of the instabilities in the jet. We have detected signatures of possible jet disruption in the jet due to the growth of instabilities, which points towards a possible morphological classification of the source as an FRI, contrary to what previously thought in terms of luminosity criteria.

Nuclear activity in galaxies is closely connected to galactic mergers and supermassive black holes (SBH). Galactic mergers perturb substantially the dynamics of gas and stellar population in the merging galaxies, and they are expected to lead to formation of supermassive binary black holes (BBH) in the center of mass of the galaxies merged. A scheme is proposed here that connects the peak magnitude of the nuclear activity with evolution of a BBH system. The scheme predicts correctly the relative fractions of different types of active galactic nuclei (AGN) and explains the connection between the galactic type and the strength of the nuclear activity. It shows that most powerful AGN should result from mergers with small mass ratios, while weaker activity is produced in unequal mergers. The scheme explains also the observed lack of galaxies with two active nuclei, which is attributed to effective disruption of accretion disks around the secondary in BBH systems with masses of the primary smaller than ~10^10 solar masses.

Past years have brought an increasingly wider recognition of the ubiquity of relativistic outflows (jets) in galactic nuclei, which has turned jets into an effective tool for investigating the physics of nuclear regions in galaxies. A brief summary is given here of recent results from studies of jets and nuclear regions in several active galaxies with prominent outflows.

Over the last five decades, AGN studies have produced a number of spectacular examples of synergies and multifaceted approaches in astrophysics. The field of AGN research now spans the entire spectral range and covers more than twelve orders of magnitude in the spatial and temporal domains. The next generation of astrophysical facilities will open up new possibilities for AGN studies, especially in the areas of high-resolution and high-fidelity imaging and spectroscopy of nuclear regions in the X-ray, optical, and radio bands. These studies will address in detail a number of critical issues in AGN research such as processes in the immediate vicinity of supermassive black holes, physical conditions of broad-line and narrow-line regions, formation and evolution of accretion disks and relativistic outflows, and the connection between nuclear activity and galaxy evolution.

Relativistic outflows represent one of the best-suited tools to probe the physics of AGN. Numerical modelling of internal structure of the relativistic outflows on parsec scales provides important clues about the conditions and dynamics of the material in the immediate vicinity of the central black holes in AGN. We investigate possible causes of the structural patterns and regularities observed in the parsec-scale jet of the well-known quasar 3C 273. We present here the results from a 3D relativistic hydrodynamics numerical simulation based on the parameters given for the jet by Lobanov & Zensus (2001), and one in which the effects of jet precession and the injection of discrete components have been taken into account. We compare the model with the structures observed in 3C 273 using very long baseline interferometry and constrain the basic properties of the flow. We find growing perturbation modes in the simulation with similar wavelengths to those observed, but with a different set of wave speeds and mode identification. If the observed longest helical structure is produced by the precession of the flow, longer precession periods should be expected. Our results show that some of the observed structures could be explained by growing Kelvin-Helmholtz instabilities in a slow moving region of the jet. However, we point towards possible errors in the mode identification that show the need of more complete linear analysis in order to interpret the observations. We conclude that, with the given viewing angle, superluminal components and jet precession cannot explain the observed structures.

The Square Kilometer Array will be operating at the same time with several new large optical, X-ray and Gamma-ray facilities currently under construction or planned. Fostering synergies in astrophysical research made across different spectral bands presents a compelling argument for designing the SKA such that it would offer imaging capabilities similar to those of other future telescopes. Imaging capabilities of the SKA are compared here with those of the major future astrophysical facilities.

The possibility of radiative effects induced by the Lorentz and CPT non-invariant interaction term for fermions in the Standard Model Extension is investigated. In particular, electron-positron photo-production and photon emission by electrons and positrons are studied. The rates of these processes are calculated in the Furry picture. It is demonstrated that the rates obtained in the framework of the model adopted strongly depend on the polarization states of the particles involved. As a result, ultra-relativistic particles produced should occupy states with a preferred spin orientation, i.e., photons have the sign of polarization opposite to the sign of the effective potential, while charged particle are preferably in the state with the helicity coinciding with the sign of the effective potential. This leads to evident spatial asymmetries which may have certain consequences observable at high energy accelerators, and in astrophysical and cosmological studies.

Relativistic outflows represent one of the best-suited tools to probe the physics of AGN. Numerical modelling of internal structure of the relativistic outflows on parsec scales provides important clues about the conditions and dynamics of the material in the immediate vicinity of the central black holes in AGN. We investigate possible causes of the structural patterns and regularities observed in the parsec-scale jet of the well known quasar 3C 273. We compare the model with the radio structure observed in 3C 273 on parsec scales using very long baseline interferometry (VLBI) and constrain the basic properties of the flow. Our results show that Kelvin-Helmholtz instabilities are the most plausible mechanism to generate the observed structures.

We present observational evidence for the link between variability of the radio emission of the relativistic jet and optical continuum emission in the radio galaxy 3C 390.3 which indicates that the source of variable non-thermal continuum radiation is located in the innermost part of the relativistic jet. We suggest that the continuum emission from the jet ionizes material in a subrelativistic outflow surrounding the jet, which results in a formation of a conical region with broad emission lines at a distance of about 0.4 parsecs from the central nucleus.

The ``central engine'' of AGN is thought to be powered by accretion on a central nucleus believed to be a super-massive black hole. The localization and exact mechanism of the energy release in AGN are still not well understood. We present observational evidence for the link between variability of the radio emission of the compact jet, optical and X-ray continua emission and ejections of new jet components in the radio galaxy 3C 390.3. The time delays between the light curves of the individual jet components and the light curve of the optical continuum are estimated by using minimization methods and the discret correlation function. We find that the variations of the optical continuum are correlated with radio emission from a stationary feature in the jet. This correlation indicates that the source of variable non-thermal continuum radiation is located in the innermost part of the relativistic jet. We suggest that the continuum emission from the jet and counterjet ionizes material in a subrelativistic outflow surrounding the jet, which results in a formation of two conical regions with broad emission lines (in addition to the conventional broad line region around the central nucleus) at a distance more or equal to 0.4 parsecs from the central engine. Implications for modeling of the broad-line regions are discussed.

The luminous high-redshift (z=2.17) quasar S50836+710 has been observed in October 1997 with the VSOP at 1.6 GHz and 5 GHz. We report here a previously unpublished image made from the data at 1.6 GHz and compare the structure of a relativistic jet in this quasazr at the two frequencies. We present a spectral index image tracing spectral properties of the jet up to ~40 milliarcsecond distance from the nucleus. The curved jet ridge line observed in the images and the spectral index distribution can be described by Kelvin-Helmholtz instability developing in a relativistic outflow with a Mach number of ~6. In this description, the overall ridge line of the jet is formed by the helical surface mode of Kelvin-Helmholtz instability, while areas of flatter spectral index embedded into the flow correspond to pressure enhancements produced by the elliptical surface mode of the instability. An alternative explanation involving a sequence of slowly dissipating shocks cannot be ruled out at this point.

Radio spectroscopy offers a number of tools for studying a large variety of astrophysical phenomena, ranging from stars and their environment to interstellar and intergalactic medium, active galactic nuclei (AGN) and distant quasars. Main targets of extragalactic radio spectroscopy are atomic and molecular material in galaxies, HII regions, and maser emission originating in the dense, circumnuclear regions. These studies cover all galactic types and span an impressive range of angular scales and distances. Molecular emission, hydrogen absorption and maser lines have become the tools of choice for making an assessment of physical conditions in the nuclear regions of galaxies. In this contribution, some of the recent advances in the aforementioned fields will be reviewed and discussed in connection with future radio astronomical facilities.

We have used VLBA fringe visibility data obtained at 15 GHz to examine the compact structure in 250 extragalactic radio sources. For 171 sources in our sample, more than half of the total flux density seen by the VLBA remains unresolved on the longest baselines. There are 163 sources in our list with a median correlated flux density at 15 GHz in excess of 0.5 Jy on the longest baselines. For about 60% of the sources, we have at least one observation in which the core component appears unresolved (generally smaller than 0.05 mas) in one direction, usually transverse to the direction into which the jet extends. BL Lacs are on average more compact than quasars, while active galaxies are on average less compact. Also, in an active galaxy the sub-milliarcsecond core component tends to be less dominant. IDV sources typically have a more compact, more core-dominated structure on sub-milliarcsecond scales than non-IDV sources, and sources with a greater amplitude of intra-day variations tend to have a greater unresolved VLBA flux density. The objects known to be GeV gamma-ray loud appear to have a more compact VLBA structure than the other sources in our sample. This suggests that the mechanisms for the production of gamma-ray emission and for the generation of compact radio synchrotron emitting features are related. The brightness temperature estimates and lower limits for the cores in our sample typically range between 10^11 and 10^13 K, but they extend up to 5x10^13 K, apparently in excess of the equipartition brightness temperature, or the inverse Compton limit for stationary synchrotron sources. The largest component speeds are observed in radio sources with high observed brightness temperatures, as would be expected from relativistic beaming (abridged).

A method is introduced to derive resolution criteria for various a priori defined templates of brightness distribution fitted to represent structures and objects in astronomical images. The method is used for deriving criteria for the minimum and maximum resolvable sizes of a template. The minimum resolvable size of a template is determined by the ratio of (SNR-1)/SNR, and the maximum detectable size is determined by the ratio of 1/SNR. Application of these criteria is discussed in connection to data from filled-aperture instruments and interferometers, accounting for different aperture shapes and the effects of Fourier sampling, tapering, apodization and visibility weighting. Specific resolution limits are derived for four different templates of brightness distribution: (1) two-dimensional Gaussian, (2) optically thin spherical shell, (3) disk of uniform brightness, and (4) ring. The limiting resolution for these templates changes with SNR similarly to the quantum limit on resolution. Practical application of the resolution limits is discussed in two examples dealing with measurements of maximum brightness temperature in compact relativistic jets and assessments of morphology of young supernova remnants.

We present a hydrodynamical 3D simulation of the relativistic jet in 3C273, in comparison to previous linear perturbation analysis of Kelvin-Helmholtz instability developing in the jet. Our aim is to assess advantages and limitations of both analytical and numerical approaches and to identify spatial and temporal scales on which the linear regime of Kelvin-Helmholtz instability can be applied in studies of morphology and kinematics of parsec-scale jets.

The dynamic evolution of binary systems of supermassive black holes (SMBH) may be a key factor affecting a large fraction of the observed properties of active galactic nuclei (AGN) and galaxy evolution. Different classes of AGN can be related in general to four evolutionary stages in a binary SMBH: 1) early merger stage; 2) wide pair stage; 3) close pair stage; and 4) pre-coalescence stage. This scheme can explain a variety of properties of AGN: radio and optical luminosity differences between different classes of AGN, long-term and short-term variability, quasi-periodic nuclear flares, recurrent formation of relativistic outflows in AGN and their apparent morphology and kinematics.

Radio loud active galactic nuclei present a remarkable variety of signs indicating the presence of periodical processes possibly originating in binary systems of supermassive black holes, in which orbital motion and precession are ultimately responsible for the observed broad-band emission variations, as well as for the morphological and kinematic properties of the radio emission on parsec scales. This scenario, applied to the quasar 3C345, explains the observed variations of radio and optical emission from the quasar, and reproduces the structural variations observed in the parsec-scale jet of this object. The binary system in 3C345 is described by two equal-mass black holes with masses of approximately 7*10^8 solar masses separated by ~0.33 pc and orbiting with a period of ~480 yr. The orbital motion induces a precession of the accretion disk around the primary black hole, with a period of ~2570 yr. The jet plasma is described by a magnetized, relativistic electron-positron beam propagating inside a wider and slower electron-proton jet. The combination of Alfven wave perturbations of the beam, the orbital motion of the binary system and the precession of the accretion disk reproduces the variability of the optical flux and evolution of the radio structure in 3C345. The timescale of quasi-periodic flaring activity in 3C345 is consistent with typical disk instability timescales. The present model cannot rule out a small-mass orbiter crossing the accretion disk and causing quasi-periodic flares.

The quasi-classical theory of the spin light of neutrino ($SL\nu$) in background matter, accounting for the neutrino polarization, is developed. The neutrino transitions $\nu_{L}\to \nu_{R}$ and $\nu_{R}\to \nu_{L}$ rates in matter are calculated. It is shown that the $SL\nu$ in matter leads to the neutrino conversion from active $\nu_{L}$ to sterile $\nu_{R}$ states (neutrino self-polarization effect in matter).

NGC 1052 offers the possibility to study the obscuring torus around a supermassive black hole, predicted by the standard model of active galactic nuclei, over a wide range of wavelengths from the radio to the X-ray regime. We present a detailed VLBI study of the parsec-scale structure of the ``twin-jet'' system in NGC 1052 in both total and polarized intensity and at multiple frequencies. We report the detection of linearly polarized emission from the base of the eastern jet at 5 GHz. While the radio spectrum in this region might be still consistent with synchrotron self absorption, the highly inverted spectrum of the western jet base represents a clear sign of pronounced free-free absorption in a circumnuclear torus. We observe an abrupt change of the brightness temperature gradient at a distance of $\sim 0.2$ pc to 0.3 pc from the central engine. This might provide an observational signature of the edge of the central torus, where the transition from an external pressure-dominated jet regime to a more or less freely expanding jet takes place. We determine the absorbing column density towards the western jet core to be $\sim 2.2 \times 10^{22}$ cm$^{-2}$ in good agreement with the values derived from various X-ray observations. This suggests that the nuclear X-ray emission and the jet emission imaged by VLBI originate on the same scales.

We have observed the QSO 3C345 with VSOP at 1.6GHz and 5GHz and with the VLBA at 22GHz and 43GHz, with 24 epochs in total, from 1997.40 to 1999.69. This study, combined with previous observations, reveals periodic changes in the parsec-scale jet. This quasi-periodicity can be found in the initial position angles, the trajectories and the flux density evolution of the jet components, with a period of about 9 years. This period should reflect the dynamics of the material in the central region of 3C345 and it is possibly related to a putative binary black hole in the center of this object.

The wealth of high quality data now available on the M87 jet inspired us to carry out a detailed analysis of the plasma physical conditions in the jet. In a companion paper (Lobanov, Hardee & Eilek, this proceedings) we identify a double-helix structure within the jet, and apply Kelvin-Helmholtz stability analysis to determine the physical state of the jet plasma. In this paper we treat the jet as a test case for in situ particle acceleration. We find that plasma turbulence is likely to exist at levels which can maintain the energy of electrons radiating in the radio to optical range, consistent with the broadband spectrum of the jet.

Here we summarize our results from a detailed multi-frequency study of the QSO 3C309.1 based on the Very Long Baseline Array (VLBA) observations made in mid 1998. From our images, we find a curved jet extending up to 100 milliarcseconds (mas) to the east at low frequencies with two main components, A and B. A preliminary astrometric analysis (Ros and Lobanov 2001) provides a determination of the core position at different frequencies by phase-referencing to a nearby radio source, QSO S5 1448+76. The changes of the core position with frequency suggest high opacity close to the core caused by synchrotron self-absorption. Due to the large astrometric uncertainties we cannot draw any conclusions about the values of the opacity gradients at high frequencies. We believe that a detailed analysis of the frequency depedence of the core position will reveal the profile of the matter distribution in the broad line region, as was initially suggested by Lobanov (1998).

We present results of a VLBA multi-frequency study of the pc-scale twin jet in NGC1052. We observed this object at epoch 1998.99 with the VLBA at 5, 8.4, 22 and 43 GHz both in total and linearly polarized intensity. The spectral analysis confirms the necessity of a free-free absorbing medium, obscuring the innermost part of both jets. At 5 GHz we found a compact linearly polarized emission region at the base of the eastern jet with a degree of polarization of 1.5%. At higher frequencies there is no evidence for polarization in our data. A core shift analysis constrains the position of the central engine to ~0.03 pc. The shift rates of the apparent core position with frequency confirm the strong influence of free-free absorption in conjunction with steep pressure gradients at the bases of both jets.

The "core" of a radio source is believed to mark the frequency-dependent location where the optical depth to synchrotron self--absorption tau_s~1. The frequency dependence can be used for derive physical conditions of the radio emitting region and the ambient environment near the central engine of the radio source. In order to test and improve the models to derive this information, we made multi-frequency dual-polarization observations of 3C309.1 in 1998.6, phase-referenced to the QSOs S5 1448+76 and 4C72.20 (S5 1520+72) (4d26' and 1d49' away, respectively). We present here preliminary results from these observations: total and polarized intensity maps, spectral information deduced from these images, and the relative position of 3C309.1 with respect to S5 1448+76 at different frequencies. Finally, we discuss briefly the observed shift of the core position.

The quasar 3C345 is one of the best examples of extragalactic jet with superluminal components moving on helical paths. To get better constraints on the jet physics, we have imaged the source at the highest possible resolution, using mm-VLBI and space-VLBI observations. We observed 3C345 at 5 different epochs so far (1 epoch with 3-mm ground VLBI, and 4 epochs with space VLBI at 6cm and 18cm). These measurements expand significantly our ongoing monitoring of 3C345. In this contribution, we discuss the properties of the jet and its moving components in the immediate vicinity of the core (<3mas).

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.

So far all relativistically boosted jets with superluminal motion have only been detected in typical radio galaxies with early type host galaxies. We have now discovered superluminal motion in the Seyfert I galaxy III Zw 2, classified as a spiral. Superluminal motion was first inferred from the spectral evolution of the source and then confirmed by VLBI observations. The lower limit for the apparent expansion speed is 1.25 +/- 0.09 c. The fact that the spectral and spatial evolution are closely linked demonstrates that we are dealing with real physical expansion. Prior to this rapid expansion we have seen a period of virtually no expansion with an expansion speed less than 0.04 c. Since III Zw 2 is also part of a sample of so called radio-intermediate quasars (RIQ), it confirms earlier predictions of superluminal motion for this source, based on the argument that RIQs could be relativistically boosted jets in radio-weak quasars and Seyfert galaxies.

We have studied the parsec-scale structure of the quasar 3C345 at \lambda\lambda 1.3, 2, 3.6, and 6 cm, using total intensity and linear polarization images from 3 epochs of VLBA observations made in 1995.84, 1996.41, and 1996.81. The images show an unresolved "core" and a radio jet extending over distances up to 25 milliarcseconds. At 1.3 and 2 cm, images of the linearly polarized emission reveal an alignment of the electric vector with the extremely curved inner section of the jet. At 6 cm, the jet shows strong fractional polarization (m~15%), with the electric vector being oriented perpendicular to the jet. This is consistent with the existence of a transition from the shock-dominated to the plasma interaction-dominated emission regimes in the jet, occurring at a distance of ~1.5 milliarcseconds from the "core".

The long-term evolution of the synchrotron emission from the parsec-scale jet in the quasar 3C345 is analysed, on the basis of multi-frequency monitoring with very long baseline interferometry (VLBI) and covering the period 1979-1994. We combine the model fits from 44 VLBI observations of 3C345 made at 8 different frequencies between 2.3 and 100GHz. We calculate the turnover frequency, turnover flux density, integrated 4-25GHz flux and 4-25GHz luminosity of the core and the moving features. The core has an estimated mean luminosity of 7.1+/-3.5 * 10^42 erg/s; the estimated total luminosity of 3C345 on parsec scales is ~3*10^43 erg/s (about 1% of the observed luminosity of the source between the radio to infrared regimes). The spectral changes in the core can be reconciled with a shock or dense plasma condensation traveling through the region where the jet becomes optically thin. We are able to describe the evolution of the core spectrum by a sequence of 5 flare-like events characterized by an exponential rise and decay of the particle number density of the material injected into the jet. The same model is also capable of predicting the changes in the flux density observed in the core. The flares occur approximately every 3.5-4 years, roughly correlating with appearances of new moving features in the jet.

We present radio observations at three frequencies and contemporaneous optical monitoring of the peculiar BL Lac object AO 0235+164. During a three-week campaign with the VLA we observed intraday variability in this source and found a distinct peak which can be identified throughout the radio frequencies and tentatively connected to the R-band variations. This event is characterized by unusual properties: its strength increases, and its duration decreases with wavelength, and it peaks earlier at 20cm than at 3.6 and 6cm. We discuss several generic models (a `standard' shock-in-jet model, a precessing beam, free-free-absorption in a foreground screen, interstellar scattering, and gravitational microlensing), and explore whether they can account for our observations. Most attempts at explaining the data on 0235+164 require an extremely small source size, which can be reconciled with the $10^{12}$ K inverse Compton limit only when the Doppler factor of the bulk flow is of order 100. However, none of the models is completely satisfactory, and we suggest that the observed variability is due to a superposition of intrinsic and propagation effects.

We monitored the superluminal QSO 3C 345 at three epochs during a one-year period in 1995--1996, observing with the VLBA at 22, 15, 8.4, and 5 GHz. We imaged the radio source both in total and in polarized intensity. In the images at 5 and 8.4 GHz, the jet emission is traced up to 20 milliarcseconds (mas) from the jet core. In the 15 and 22 GHz images, we identify several enhanced emission regions moving at apparent speeds of 5c. Images of the linear polarized emission show predominantly an alignment of the electric vector with the extremely curved jet along the inner part of the high frequency jet. At 5 GHz, the jet shows remarkably strong fractional polarization (m~15%) with the electric vector perpendicular to the jet orientation.

Broad band spectra of total and compact--scale radio emission from blazars, used in combination with kinematic information inferred from VLBI monitoring programs, can be applied successfully to better constrain the models of compact jets. We discuss here a ``hands-on'' approach for tying together the kinematic and spectral properties of radio emission from blazars.

III Zw 2 is a spiral galaxy with an optical spectrum and faint extended radio structure typical of a Seyfert galaxy, but also with an extremely variable, blazar-like radio core. We have now discovered a new radio flare where the source has brightened more than twenty-fold within less than two years. A broad-band radio spectrum between 1.4 and 666 GHz shows a textbook-like synchrotron spectrum peaking at 43 GHz, with a self-absorbed synchrotron spectral index +2.5 at frequencies below 43 GHz and an optically thin spectral index -0.75 at frequencies above 43 GHz. The outburst spectrum can be well fitted by two homogenous, spherical components with equipartition sizes of 0.1 and 0.2 pc at 43 and 15 GHz, and with magnetic fields of 0.4 and 1 Gauss. VLBA observations at 43 GHz confirm this double structure and these sizes. Time scale arguments suggest that the emitting regions are shocks which are continuously accelerating particles. This could be explained by a frustrated jet scenario with very compact hotspots. Similar millimeter-peaked spectrum (MPS) sources could have escaped our attention because of their low flux density at typical survey frequencies and their strong variability.

We report results on the QSO 3C345 obtained from VLBA observations carried out on epochs 1995.84, 1996.41, and 1996.81, at 22, 15, 8.4 and 5 GHz. Total and polarized images VLBI images were obtained. We observed superluminal motions of components C8 and C7 with respect to the core component D. We imaged also the remarkably complex polarization structure near the core.

The luminous high-redshift (z=2.17) quasar S5 0836+710 has been observed at 5GHz with the VSOP. We compare the properties of three images obtained from the observation: a low-resolution ground array image (dynamic range 4600:1), a full-resolution VSOP image (900:1), and an image made with only the space baselines (200:1). The space baselines alone are sufficient for a reliable recovery of the source structure, within the limits of the achieved spatial sampling of the visibility data. The curved jet ridge line observed in the images can be described by Kelvin-Helmholtz instabilities developing in a relativistic outflow with the Mach number of about 6. This description holds on the scales of up to 700h^-1 pc, and is shown to be consistent with variable apparent speeds observed in the jet.

We discuss a technique for mapping the synchrotron turnover frequency distribution using nearly simultaneous, multi--frequency VLBI observations. The limitations of the technique arising from limited spatial sampling and frequency coverage are investigated. The errors caused by uneven spatial sampling of typical multi--frequency VLBA datasets are estimated through numerical simulations, and are shown to be of the order of 10%, for pixels with the deconvolution SNR~7. The fitted spectral parameters are corrected for the errors due to limited frequency coverage of VLBI data. First results from mapping the turnover frequency distribution in 3C345 are presented.

Broad-band (gamma to radio) variations of the flux density were observed in the first half of 1992 in the luminous high redshift (z = 2.172) quasar S5 0836+710. VLBI monitoring observations during 1993 -- 1996 performed at 86 GHz, 22 GHz, 15 GHz, and 8 GHz show the ejection of a new jet component, which most probably is directly related to a quasi simultaneous gamma-, X-ray, optical flaring activity which was observed in February 1992. During the period 1992 -- 1993 the flaring propagated through the radio spectrum. From several quasi-simultaneous radio spectra taken during this phase of activity, we determine the time evolution of the spectral turnover of the radio spectrum in the S_m-\nu_m diagram. The data indicate a correlation of the jet activity with the variability of the broad-band electromagnetic spectrum of the source. The observational findings are discussed in the framework of relativistic shock models.

We study the properties of ultracompact jets in several prominent radio sources (Cygnus A, 3C 309.1, 3C 345, 3C 395, 4C 39.25, and 1038+528 A), using the frequency dependence of observed position of the optically thick jet core. Frequency dependent offsets of the core positions are used for calculating the luminosities, magnetic fields, and geometrical properties of the jets. Opacity effects are studied, for both synchrotron self-absorption and external free-free absorption. Pressure and density gradients in the jets and in surrounding ambient medium are shown to be a primary factor determining the observed properties of ultracompact jets. We discuss possible applications of opacity effects to studying the conditions existing in central regions of active galactic nuclei.