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

We present 23 interferometric images of parsec-scale jet of the quasar PKS 1741--03 at 15, 24 and 43 GHz spanning about 13 yr. We model the images as a superposition of discrete two--dimensional elliptical Gaussian components, with parameters determined by the cross--entropy technique. All the images present a spatially unresolved component (core) and usually two or three components receding from it. The same components were found in simultaneous 24 and 43 GHz maps, showing the robustness of our model-fitting. The core-shift opacity effect between these frequencies is weak. We have identified seven components moving along straight lines at constant apparent superluminal speeds ($3.5\le\beta_\mathrm{obs}\le 6.1$), with different sky position angles ($-186\deg\le \eta\le -125\deg$). The core flux density tracks quite well the fluctuations seen in the historical single-dish light curve at 14.5 GHz, with no measurable delay. The total flux density from the moving jet components is delayed $\sim$2 yr in relation to the core light curve, roughly the same as the lag between the ejection epoch and the maximum flux density in the light curves of the jet components. We propose that there are three non-exclusive mechanisms for producing these delays. From the kinematics of the most robust jet components and the core brightness temperature, we determined the bulk Lorentz factor ($4.8\le\gamma\le 24.5$) and the jet viewing angle ($0.\deg 35 \le \theta \le 4.\deg 2$); these values agree with previous estimates from the spectral energy distribution of PKS 1741--03 and its radio variability.

Modeling VLBI ejections of nuclei of extragalactic radio sources, indicates that their nuclei contain a binary black hole system. One can derive the distance and the positions of the two black holes in the plane of the sky. We can also use the RMS of the time series of the ICRF2 survey to obtain an estimate of the structure and the size of the nuclei. We will discuss the possible problems to link VLBI observations and GAIA optical observations of radio quasars if they contain a binary black hole system.

If we assume that nuclei of extragalactic radio sources contain binary black hole systems, the two black holes can eject VLBI components in which case two families of different VLBI trajectories will be observed. Another important consequence of a binary black hole system is that the VLBI core is associated with one black hole, and if a VLBI component is ejected by the second black hole, one expects to be able to detect the offset of the origin of the VLBI component ejected by the black hole that is not associated with the VLBI core. The ejection of VLBI components is perturbed by the precession of the accretion disk and the motion of the black holes around the center of gravity of the binary black hole system. We modeled the ejection of the component taking into account the two pertubations and present a method to fit the coordinates of a VLBI component and to deduce the characteristics of the binary black hole system. Specifically, this is the ratio Tp/Tb where Tp is the precession period of the accretion disk and Tb is the orbital period of the binary black hole system, the mass ratio M1/M2, and the radius of the Binary Black Hole system Rbin. From the variations of the coordinates as a function of time of the ejected VLBI component, we estimated the inclination angle io and the bulk Lorentz factor gamma of the modeled component. We applied the method to component S1 of 1823+568 and to component C5 of 3C 279.

We present observations of the European Large-Area {\it ISO} Survey-North 1 (ELAIS-N1) at 325 MHz using the Giant Metrewave Radio Telescope (GMRT), with the ultimate objective of identifying active galactic nuclei and starburst galaxies and examining their evolution with cosmic epoch. After combining the data from two different days we have achieved a median rms noise of $\approx40 \mu$Jy beam$^{-1}$, which is the lowest that has been achieved at this frequency. We detect 1286 sources with a total flux density above $\approx270 \mu$Jy. In this paper, we use our deep radio image to examine the spectral indices of these sources by comparing our flux density estimates with those of Garn et al. at 610 MHz with the GMRT, and surveys with the Very Large Array at 1400 MHz. We attempt to identify very steep spectrum sources which are likely to be either relic sources or high-redshift objects as well as inverted-spectra objects which could be Giga-Hertz Peaked Spectrum objects. We present the source counts, and report the possibility of a flattening in the normalized differential counts at low flux densities which has so far been reported at higher radio frequencies.

We present a new method to fit the variations of both coordinates of a VLBI component as a function of time, assuming that the nucleus of the radio source contains a binary black hole system (BBH system). The presence of a BBH system produces 2 perturbations of the trajectory of the ejected VLBI components. By using only the VLBI coordinates, the problem we have to solve reduces to an astrometric problem. Knowledge of the variations of the VLBI coordinates as a function of time contains the kinematical information, thus we are able to deduce the inclination angle of the source and the bulk Lorentz factor of the ejected component. Generally, there is a family of the BBH system producing the same fit to our data. To illustrate this method, we apply it to the source 1807+784. We find that the inclination of the source is i = 5.8+-1.8 degrees and the VLBI component is ejected with a bulk Lorentz factor of 3.7+-0.3. We determine the family of the BBH system which provides the best fit, assuming at first that the masses of the 2 black holes are equal and then that the masses are different. Each family of BBH systems is characterized by Tp/Tb~1.967, where Tp and Tb are the precession period of the accretion disk of the black hole ejecting the VLBI component and the orbiting period of the BBH system.

We present the results of a multi-frequency analysis of the structural variability in the parsec-scale jet of the blazar S5 1803+784. More than 90 epochs of observations at 6 frequencies from 1.6 GHz up to 22 GHz have been combined and analyzed. We discuss an alternative jet model for the source. In contrast to previously discussed motion scenarios for S5 1803+784, we find that the jet structure within 12 mas of the core can most easily be described by seven "oscillating" jet features. We find that the parameters of jet features, such as core separation, position angle and flux density, change in a periodic way with a timescale of about 4 years. We also find evidence for a correlation between these parameters and the total flux density variability. We suggest a scenario incorporating a periodic form of motion (e.g. rotation, precession), with a non-negligible geometrical contribution to explain the observational results.

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.