We present deep Near-infrared (NIR) images of a sample of 19 intermediate-redshift ($0.3<z<1.0$) radio-loud active galactic nuclei (AGN) with powerful relativistic jets ($L_{1.4GHz} >10^{27}$ WHz$^{-1}$), previously classified as flat-spectrum radio quasars. We also compile host galaxy and nuclear magnitudes for blazars from literature. The combined sample (this work and compilation) contains 100 radio-loud AGN with host galaxy detections and a broad range of radio luminosities $L_{1.4GHz} \sim 10^{23.7} - 10^{28.3}$~WHz$^{-1}$, allowing us to divide our sample into high-luminosity blazars (HLBs) and low-luminosity blazars (LLBs). The host galaxies of our sample are bright and seem to follow the $\mu_{e}$-$R_{eff}$ relation for ellipticals and bulges. The two populations of blazars show different behaviours in the \mnuc - \mbulge plane, where a statistically significant correlation is observed for HLBs. Although it may be affected by selection effects, this correlation suggests a close coupling between the accretion mode of the central supermassive black hole and its host galaxy, that could be interpreted in terms of AGN feedback. Our findings are consistent with semi--analytical models where low--luminosity AGN emit the bulk of their energy in the form of radio jets, producing a strong feedback mechanism, and high--luminosity AGN are affected by galaxy mergers and interactions, which provide a common supply of cold gas to feed both nuclear activity and star formation episodes.

We address the highly debated issue of constraining the gamma-ray emission region in blazars from cross-correlation analysis using discrete correlation function between radio and gamma-ray light curves. The significance of the correlations is evaluated using two different approaches: simulating light curves and mixed source correlations. The cross-correlation analysis yielded 26 sources with significant correlations. In most of the sources, the gamma-ray peaks lead the radio with time lags in the range +20 and +690 days, whereas in sources 1633+382 and 3C 345 we find the radio emission to lead the gamma rays by -15 and -40 days, respectively. Apart from the individual source study, we stacked the correlations of all sources and also those based on sub-samples. The time lag from the stacked correlation is +80 days for the whole sample and the distance travelled by the emission region corresponds to 7 pc. We also compared the start times of activity in radio and gamma rays of the correlated flares using Bayesian block representation. This shows that most of the flares at both wavebands start at almost the same time, implying a co-spatial origin of the activity. The correlated sources show more flares and are brighter in both bands than the uncorrelated ones.

We explore the connection between the black hole mass and its relativistic jet for a sample of radio-loud AGN (z < 1), in which the relativistic jet parameters are well estimated by means of long term monitoring with the 14m Mets\"ahovi millimeter wave telescope and the Very Long Base-line Array (VLBA). NIR host galaxy images taken with the NOTCam on the Nordic Optical Telescope (NOT) and retrieved from the 2MASS all-sky survey allowed us to perform a detailed surface brightness decomposition of the host galaxies in our sample and to estimate reliable black hole masses via their bulge luminosities. We present early results on the correlations between black hole mass and the relativistic jet parameters. Our preliminary results suggest that the more massive the black hole is, the faster and the more luminous jet it produces.

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.

The blazar 1156+295 was active at gamma-ray energies, exhibiting three prominent flares during the year 2010. Here, we present results using the combination of broadband (X-ray through mm single dish) monitoring data and radio band imaging data at 43 GHz on the connection of gamma-ray events to the ejections of superluminal components and other changes in the jet of 1156+295. The kinematics of the jet over the interval 2007.0-2012.5 using 43 GHz Very Long Baseline Array observations, reveal the presence of four moving and one stationary component in the inner region of the blazar jet. The propagation of the third and fourth components in the jet corresponds closely in time to the active phase of the source in gamma rays. We briefly discuss the implications of the structural changes in the jet for the mechanism of gamma-ray production during bright flares. To localise the gamma-ray emission site in the blazar, we performed the correlation analysis between the 43 GHz radio core and the gamma-ray light curve. The time lag obtained from the correlation constrains the gamma-ray emitting region in the parsec-scale jet.

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.

The coexistence of Planck and Fermi satellites in orbit has enabled the exploration of the connection between the (sub-)millimeter and gamma-ray emission in a large sample of blazars. We find that the gamma-ray emission and the (sub-)mm luminosities are correlated over five orders of magnitude. However, this correlation is not significant at some frequency bands when simultaneous observations are considered. The most significant statistical correlations, on the other hand, arise when observations are quasi-simultaneous within 2 months. Moreover, we find that sources with an approximate spectral turnover in the middle of the mm-wave regime are more likely to be strong gamma-ray emitters. These results suggest a physical relation between the newly injected plasma components in the jet and the high levels of gamma-ray emission.

Spectral energy distributions (SEDs) and radio continuum spectra are presented for a northern sample of 104 extragalactic radio sources, based on the Planck Early Release Compact Source Catalogue (ERCSC) and simultaneous multifrequency data. The nine Planck frequencies, from 30 to 857 GHz, are complemented by a set of simultaneous observations ranging from radio to gamma-rays. This is the first extensive frequency coverage in the radio and millimetre domains for an essentially complete sample of extragalactic radio sources, and it shows how the individual shocks, each in their own phase of development, shape the radio spectra as they move in the relativistic jet. The SEDs presented in this paper were fitted with second and third degree polynomials to estimate the frequencies of the synchrotron and inverse Compton (IC) peaks, and the spectral indices of low and high frequency radio data, including the Planck ERCSC data, were calculated. SED modelling methods are discussed, with an emphasis on proper, physical modelling of the synchrotron bump using multiplecomponents. Planck ERCSC data also suggest that the original accelerated electron energy spectrum could be much harder than commonly thought, with power-law index around 1.5 instead of the canonical 2.5. The implications of this are discussed for the acceleration mechanisms effective in blazar shock. Furthermore in many cases the Planck data indicate that gamma-ray emission must originate in the same shocks that produce the radio emission.

We compare the gamma-ray photon flux variability of northern blazars in the Fermi/LAT First Source Catalog with 37 GHz radio flux density curves from the Metsahovi quasar monitoring program. We find that the relationship between simultaneous millimeter (mm) flux density and gamma-ray photon flux is different for different types of blazars. The flux relation between the two bands is positively correlated for quasars and does no exist for BLLacs. Furthermore, we find that the levels of gamma-ray emission in high states depend on the phase of the high frequency radio flare, with the brightest gamma-ray events coinciding with the initial stages of a mm flare. The mean observed delay from the beginning of a mm flare to the peak of the gamma-ray emission is about 70 days, which places the average location of the gamma-ray production at or downstream of the radio core. We discuss alternative scenarios for the production of gamma-rays at distances of parsecs along the length of the jet

Although the Fermi mission has increased our knowledge of gamma-ray AGN, many questions remain, such as the site of gamma-ray production, the emission mechanism, and the factors that govern the strength of the emission. Using data from a high radio band, 37 GHz, uncontaminated by other radiation components besides the jet emission, we study these questions with averaged flux densities over the the first year of Fermi operations. We look for possible correlations between the 100 MeV - 100 GeV band used by the Fermi satellite and 37 GHz radio band observed at the Aalto University Metsahovi Radio Telescope, as well as for differences between the gamma-ray emission of different AGN subsamples. We use data averaged over the 1FGL period. Our sample includes 249 northern AGN, including a complete sample of 68 northern AGN with a measured average flux density exceeding 1 Jy. We find significant correlation between both the flux densities and luminosities in gamma and radio bands. The Fermi luminosity is inversely correlated with the peak frequency of the synchrotron component of the AGN spectral energy distributions. We also calculate the gamma dominances, defined as the ratio between the gamma and radio flux densities, and find an indication that high-energy blazars are more gamma-dominated than low-energy blazars. After studying the distributions of gamma and radio luminosities, it is clear that BL Lacertae objects are different from quasars, with significantly lower luminosities. It is unclear whether this is an intrinsic difference, an effect of variable relativistic boosting across the synchrotron peak frequency range, or the result of Fermi being more sensitive to hard spectrum sources like BL Lacertae objects. Our results suggest that the gamma radiation is produced co-spatially with the 37 GHz emission, i.e., in the jet.

We compare the gamma-ray photon flux variability of northern blazars in the Fermi/LAT First Source Catalog with 37 GHz radio flux density curves from the Metsahovi quasar monitoring program. We find that the relationship between simultaneous millimeter (mm) flux density and gamma-ray photon flux is different for different types of blazars. The flux relation between the two bands is positively correlated for quasars and does not exist for BLLacs. Furthermore, we find that the levels of gamma-ray emission in high states depend on the phase of the high frequency radio flare, with the brightest gamma-ray events coinciding with the initial stages of a mm flare. The mean observed delay from the beginning of a mm flare to the peak of the gamma-ray emission is about 70 days, which places the average location of the gamma-ray production at or downstream of the radio core. We discuss alternative scenarios for the production of gamma-rays at distances of parsecs along the length the jet.

We study the jet kinematics of the blazar S5 0716+714 during its active state in 2003-2004 with multi-epoch VLBI observations. Aims. We present a kinematic analysis of the large-scale (0-12 mas) jet of 0716+714, based on the results of six epochs of VLBA monitoring at 5 GHz. Additionally, we compare kinematic results obtained with two imaging methods based on different deconvolution algorithms. The blazar 0716+714 has a diffuse large-scale jet and a very faint bright compact core. Experiments with simulated data showed that the conventional data reduction procedure based on the CLEAN deconvolution algorithm does not perform well in restoring this type of structure. This might be the reason why previous kinematic studies of this source yielded ambiguous results. In order to obtain accurate kinematics of this source, we independently applied two imaging techniques to the raw data: the conventional method, based on difference mapping, which uses CLEAN deconvolution, and the generalized maximum entropy method (GMEM) realized in the VLBImager package developed at the Pulkovo Observatory in Russia. The results of both methods give us a consistent kinematic scenario: the large-scale jet of 0716+714 is diffuse and stationary. Differences between the inner (0-1 mas) and outer (1-12 mas) regions of the jet in brightness and velocity of the components could be explained by the bending of the jet, which causes the angle between the jet direction and the line of sight to change from ~5 deg to ~11 deg. For the source 0716+714 both methods worked at the limit of their capability.

We investigate the relationship between black hole mass (MBH) and Doppler boosted emission for BL Lacertae type objects (BL Lacs) detected in the SDSS and FIRST surveys. The synthesis of stellar population and bidimensional decomposition methods allows us to disentangle the components of the host galaxy from that of the nuclear black hole in their optical spectra and images, respectively. We derive estimates of black hole masses via stellar velocity dispersion and bulge luminosity. We find that masses delivered by both methods are consistent within errors. There is no difference between the black hole mass ranges for high-synchrotron peaked BL Lacs (HBL) and low-synchrotron peaked BL Lacs (LBL). A correlation between the black-hole mass and radio, optical and X-ray luminosity has been found at a high significance level. The optical-continuum emission correlates with the jet luminosity as well. Besides, X-ray and radio emission are correlated when HBLs and LBLs are considered separately. Results presented in this work: (i) show that the black hole mass does not decide the SED shapes of BL Lacs, (ii) confirm that X-ray and optical emission is associated to the relativistic jet, and (iii) present evidence of a relation between MBH and Doppler boosted emission, which among BL Lacs may be understood as a close relation between faster jets and more massive black holes.

Radio-bright BL Lacertae objects (BLOs) are typically variable and exhibit prominent flaring. We use a sample of 24 BLOs to get a clear idea of their flaring behavior and to find possible commonalities in their variability patterns. Our goal was to compare the results given by computational time scales and the observed variability parameters determined directly from the flux curves. Also, we wanted to find out if the BLO flares adhere to the generalized shock model. We use long-term monitoring data from 4.8, 8, 14.5, 22, 37, 90 and 230 GHz. The structure function, discrete correlation function and Lomb-Scargle periodogram time scales, calculated in a previous study, are analyzed in more detail. We determine flare durations, rise and decay times, absolute and relative peak fluxes from the monitoring data. We find that BLOs demonstrate a wide range of variability behavior. BLOs include sources with fast and strong variability, such as OJ 287, PKS 1749+096 and BL Lac, but also sources with more rolling fluctuations like PKS 0735+178. The most extreme flares can last for up to 13 years or have peak fluxes of approximately 12 Jy in the observer's frame. When the Doppler boosting effect is taken into account, the peak flux of a flare does not depend on the duration of the flare. A rough analysis of the time lags and peak flux evolution indicates that BLO flares in the mm - cm wavelengths are high-peaking, i.e., are in the adiabatic stage. Thus, the results concur with the generalized shock model.

We observed the blazar \object{0716+714} with the VLBA during its active state in 2003-2004. In this paper we discuss multi-frequency analysis of the inner jet (first 1 mas) kinematics. The unprecedentedly dense time sampling allows us to trace jet components without misidentification and to calculate the component speeds with good accuracy. In the smooth superluminal jet we were able to identify and track three components over time moving outwards with relatively high apparent superluminal speeds (8.5-19.4 $c$), which contradicts the hypothesis of a stationary oscillating jet in this source. Component ejections occur at a relatively high rate (once in two months), and they are accompanied by mm-continuum outbursts. Superluminal jet components move along wiggling trajectories, which is an indication of actual helical motion. Fast proper motion and rapid decay of the components suggest that this source should be observed with the VLBI at a rate of at least once in one or two months in order to trace superluminal jet components without confusion.

The blazar sequence is a scenario in which the bolometric luminosity of the blazar governs the appearance of its spectral energy distribution. The most prominent result is the significant negative correlation between the synchrotron peak frequencies and the synchrotron peak luminosities of the blazar population. Observational studies of the blazar sequence have, in general, neglected the effect of Doppler boosting. We study the dependence of both the synchrotron peak frequency and luminosity with Doppler-corrected quantities. We determine the spectral energy distributions of 135 radio-bright AGN and find the best-fit parabolic function for the distribution to quantify their synchrotron emission. The corresponding measurements of synchrotron peak luminosities and frequencies are Doppler-corrected with a new set of Doppler factors calculated from variability data. The relevant correlations for the blazar sequence are determined for these intrinsic quantities. The Doppler factor depends strongly on the synchrotron peak frequency, the lower energy sources being more boosted. Applying the Doppler correction to the peak frequencies and luminosities annuls the negative correlation between the two quantities, which becomes positive. For BL Lacertae objects, the positive correlation is particularly strong. The blazar sequence, when defined as the anticorrelation between the peak frequency and luminosity of the synchrotron component of the spectral energy distribution, disappears when the intrinsic, Doppler-corrected values are used. It is an observational phenomenon created by variable Doppler boosting across the synchrotron peak frequency range.

Multi-frequency VLBI observations allow studies of the continuum spectrum in the different parts of the parsec scale jets of AGN, providing information on the physical properties of the plasma and magnetic fields in them. Since VLBI networks cannot be scaled, the range of spatial frequencies observed differs significantly between the different observing frequencies, which makes it difficult to obtain a broadband spectrum of the individual emission features in the jet. In this paper we discuss a model-fitting based spectral extraction method, which can significantly relieve this problem. The method uses a priori knowledge of the source structure, measured at high frequencies, to allow at lower frequencies the derivation of the sizes and flux densities of even those emission features that have mutual separations significantly less than the Rayleigh limit at the given frequency. We have successfully used this method in the analysis of 5-86 GHz VLBA data of 3C273. The spectra and sizes of several individual jet features were measured, thus allowing derivation of the magnetic flux density and the energy density of the relativistic electrons in the different parts of the jet. We discuss the results, which include e.g. a detection of a strong gradient in the magnetic field across the jet of 3C273.

We have calculated variability Doppler boosting factors, Lorentz factors, and viewing angles for a large sample of sources by using total flux density observations at 22 and 37 GHz and VLBI data. We decomposed the flux curves into exponential flares and determined the variability brightness temperatures of the fastest flares. By assuming the same intrinsic brightness temperature for each source, we calculated the Doppler boosting factors for 87 sources. In addition we used new apparent jet speed data to calculate the Lorentz factors and viewing angles for 67 sources. We find that all quasars in our sample are Doppler-boosted and that the Doppler boosting factors of BL Lacertae objects are lower than of quasars. The new Lorentz factors are about twice as high as in earlier studies, which is mainly due to higher apparent speeds in our analyses. The jets of BL Lacertae objects are slower than of quasars. There are some extreme sources with very high derived Lorentz factors of the order of a hundred. These high Lorentz factors could be real. It is also possible that the sources exhibit such rapid flares that the fast variations have remained undetected in monitoring programmes, or else the sources have a complicated jet structure that is not amenable to our simple analysis. Almost all the sources are seen in a small viewing angle of less than 20 degrees. Our results follow the predictions of basic unification schemes for AGN.

A large number of AGN have been monitored for nearly 30 years at 22, 37 and 87 GHz in Mets\"ahovi Radio Observatory. These data were combined with lower frequency 4.8, 8.0 and 14.5 GHz data from the University of Michigan Radio Astronomy Observatory, higher frequency data at 90 and 230 GHz from SEST, and supplementary higher frequency data from the literature to study the long-term variability of a large sample of AGN. Both the characteristics of individual flares from visual inspection and statistically-determined variability timescales as a function of frequency and optical class type were determined. Based on past behaviour, predictions of sources expected to exhibit large flares in 2008--2009 appropriate for study by GLAST and other instruments are made. The need for long-term data for properly understanding source behaviour is emphasised.

We have studied the flare characteristics of 55 AGN at 8 different frequency bands between 4.8 and 230 GHz. Our extensive database enables us to study the various observational properties of flares in these sources and compare our results with theoretical models. We visually extracted 159 individual flares from the flux density curves and calculated different parameters, such as the peak flux density and duration, in all the frequency bands. The selection of flares is based on the 22 and 37 GHz data from Mets\"ahovi Radio Observatory and 90 and 230 GHz data from the SEST telescope. Additional lower frequency 4.8, 8, and 14.5 GHz data are from the University of Michigan Radio Observatory. We also calculated variability indices and compared them with earlier studies. The observations seem to adhere well to the shock model, but there is still large scatter in the data. Especially the time delays between different frequency bands are difficult to study due to the incomplete sampling of the higher frequencies. The average duration of the flares is 2.5 years at 22 and 37 GHz, which shows that long-term monitoring is essential for understanding the typical behaviour in these sources. It also seems that the energy release in a flare is independent of the duration of the flare.

We present first results from a multifrequency VLBA observations of 3C273 in 2003. The source was observed simultaneously at 5.0, 8.4, 15.3, 22.2, 43.2 and 86.2 GHz, and from this multifrequency data set, spectra of 16 emission features in the parsec scale jet were carefully constructed by using a new model-fitting based method. The measured spectra and sizes of the emission features were used to calculate the magnetic field density and the energy density of the relativistic electrons in the different parts of the parsec scale jet, independent of any equipartition assumption. We measure magnetic field density of an order of 1 Gauss in the core. The magnetic energy density in the core dominates over that of the relativistic electrons, while in the downstream region our data are roughly consistent with an equipartition. A strong gradient in the magnetic field density across the jet width, coincident with a transverse velocity structure at about 1.5 mas from the core, was found: the slower superluminal component B2 on the northern side of the jet has a magnetic field density two orders of magnitude lower than the faster southern components B3 and B4.

Aims: In this paper we study whether the shock-in-jet model, widely used to explain the outbursting behaviour of quasars, can be used to explain the radio flaring behaviour of the microquasar Cygnus X-3. Method: We have used a method developed to model the synchrotron outbursts of quasar jets, which decomposes multifrequency lightcurves into a series of outbursts. The method is based on the Marscher & Gear (1985) shock model, but we have implemented the modifications to the model suggested by Bjornsson & Aslaksen (2000), which make the flux density increase in the initial phase less abrupt. We study the average outburst evolution as well as specific characteristics of individual outbursts and physical jet properties of Cyg X-3. Results: We find that the lightcurves of the February-March 1994 and September 2001 outbursts can be described with the modified shock model. The average evolution shows that instead of the expected synchrotron plateau, the flux density is still increasing during the synchrotron stage. We also find that high frequency peaking outbursts are shorter in duration than the ones peaking at lower frequencies. Finally, we show that the method can be used, complementary to radio interferometric jet imaging, for deriving the physical parameters such as the magnetic field strength and the energy density of relativistic electrons in the jet of Cyg X-3.

We present a study of variability time scales in a large sample of Active Galactic Nuclei at several frequencies between 4.8 and 230 GHz. We investigate the differences of various AGN types and frequencies and correlate the measured time scales with physical parameters such as the luminosity and the Lorentz factor. Our sample consists of both high and low polarization quasars, BL Lacertae objects and radio galaxies. The basis of this work is the 22 GHz, 37 GHz and 87 GHz monitoring data from the Metsahovi Radio Observatory spanning over 25 years. In addition,we used higher 90 GHz and 230 GHz frequency data obtained with the SEST-telescope between 1987 and 2003. Further lower frequency data at 4.8 GHz, 8 GHz and 14.5 GHz from the University of Michigan monitoring programme have been used. We have applied three different statistical methods to study the time scales: The structure function, the discrete correlation function and the Lomb-Scargle periodogram. We discuss also the differences and relative merits of these three methods. Our study reveals that smaller flux density variations occur in these sources on short time scales of 1-2 years, but larger outbursts happen quite rarely, on the average only once in every 6 years. We do not find any significant differences in the time scales between the source classes. The time scales are also only weakly related to the luminosity suggesting that the shock formation is caused by jet instabilities rather than the central black hole.

We present 37 GHz data obtained at Metsahovi Radio Observatory in 2001 December - 2005 April for a large sample of BL Lacertae objects. We also report the mean variability indices and radio spectral indices in frequency intervals 5 - 37 GHz and 37 - 90 GHz. Approximately 34 % of the sample was detected at 37 GHz, 136 BL Lacertae objects in all. A large majority of the detected sources were low-energy BL Lacs (LBLs). The variability index values of the sample were diverse, the mean fractional variability of the sample being \Delta S_2 = 0.31. The spectral indices also varied widely, but the average radio spectrum of the sample sources is flat. Our observations show that many of the high-energy BL Lacs (HBL), which are usually considered radio-quiet, can at times be detected at 37 GHz.

The IX Russian-Finnish Symposium on Radio Astronomy was held in Special astrophysical observatory RAS in Nizhnij Arkhyz, Russia on 15-20 October 2006. It was dedicated to the two-side collaboration in the field of multi-wavelength investigations of solar radio emission, studies of polar regions of the Sun, studies of stellar activity, AGNs, quasars and BL Lac objects in radio bands from millimeter to decimetre wavelengths with RATAN-600, Metsahovi 14m, RT32m radio telescope and VLBI systems. Here abstracts of all forty papers are given. The Web-site of the Symposium is http://cats.sao.ru/~satr/RFSymp/ .

We study the synchrotron flaring behaviour of the blazar 3C279 based on an extensive dataset covering 10 years of monitoring at 19 different frequencies in the radio-to-optical range. The properties of a typical outburst are derived from the observations by decomposing the 19 lightcurves into a series of self-similar events. This analysis is achieved by fitting all data simultaneously to a succession of outbursts defined according to the shock-in-jet model of Marscher & Gear (1985). We compare the derived properties of the synchrotron outbursts in 3C279 to those obtained with a similar method for the quasar 3C273. It is argued that differences in the flaring behaviour of these two sources are intrinsic to the sources themselves rather than being due to orientation effects. We also compare the start times and flux densities of our modelled outbursts with those measured from radio components identified in Very Long Baseline Interferometry (VLBI) images. We find VLBI counterparts for most of our model outbursts, although some high-frequency peaking events are not seen in the radio maps. Finally, we study the link between the appearance of a new synchrotron component and the EGRET gamma-ray state of the source at 10 different epochs. We find that an early-stage shock component is always present during high gamma-ray states, while in low gamma-ray states the time since the onset of the last synchrotron outburst is significantly longer. This statistically significant correlation supports the idea that gamma-ray flares are associated with the early stages of shock components propagating in the jet. We note, however, that the shock wave is already beyond the broad line region during the gamma-ray flaring.

We report the discovery of an extremely curved jet in the radio-loud quasar PKS2136+141. Multi-frequency Very Long Baseline Array (VLBA) images show a bending jet making a turn-around of 210 degrees in the plane of the sky, which is, to our knowledge, the largest ever observed change in the position angle of an astrophysical jet. Images taken at six different frequencies, from 2.3 to 43 GHz, reveal a spiral-like trajectory, which is likely a sign of an intrinsic helical geometry. A space-VLBI image, taken with the HALCA satellite at 5 GHz and having comparable resolution to our ground-based 15 GHz data, confirms that the bend is a frequency-independent structure. VLBA monitoring data at 15 GHz, covering eight years of observations, show knots in the jet clearly deviating from ballistic motion, which suggests that the bending may be caused by a growing helical Kelvin-Helmholtz normal mode. The jet appearance suggests a helical wave at a frequency well below the "resonant" frequency of the jet, which indicates that the wave is driven by a periodic perturbation at the base of the jet. We fit the observed structure in the source with a helical twist, and we find that a simple isothermal model with a constant wave speed and wavelength gives a good fit. The measured apparent velocities indicate some degree of acceleration along the jet, which together with an observed change in the apparent half-opening angle of the jet allow us to estimate the changes in the angle between the local jet direction and our line of sight. We suggest that the jet in PKS2136+141 is distorted by a helical Kelvin-Helmholtz normal mode externally driven into the jet (e.g. by precession), and that our line of sight falls within the opening angle of the helix cone.

PG 1553+11 was the target of a coordinated three week multiwavelength campaign during 2003 April and May. A significant X-ray flare was observed during the second half of this campaign. Although no optical flare was recorded during the X-ray campaign, optical observations obtained immediately prior to the campaign displayed a higher flux than that recorded during the campaign. An optical flare was observed a few days after the end of the X-ray campaign and may be related to the X-ray flare. Radio observations were made at three frequencies, with no significant changes in flux detected near the times of the optical and X-ray flares. The spectral energy distributions and flux ratios in different wavebands observed for this object are compared to other X-ray selected blazars to demonstrate how PG 1553+11 is an extreme member of this group.

In this first of a series of papers describing polarimetric multifrequency Very Long Baseline Array (VLBA) monitoring of 3C 273 during a simultaneous campaign with the INTEGRAL gamma-ray satellite in 2003, we present 5 Stokes I images and source models at 7 mm. We show that a part of the inner jet (1-2 milliarcseconds from the core) is resolved in a direction transverse to the flow, and we analyse the kinematics of the jet within the first 10 mas. Based on the VLBA data and simultaneous single-dish flux density monitoring, we determine an accurate value for the Doppler factor of the parsec scale jet, and using this value with observed proper motions, we calculate the Lorentz factors and the viewing angles for the emission components in the jet. Our data indicates a significant velocity gradient across the jet with the components travelling near the southern edge being faster than the components with more northern path. We discuss our observations in the light of jet precession model and growing plasma instabilities.

We have collected a large amount of multifrequency data for the objects in the Metsahovi Radio Observatory BL Lacertae sample and computed their spectral energy distributions (SED). This is the first time the SEDs of BL Lacs have been studied with a sample of over 300 objects. The synchrotron components of the SEDs were fitted with a parabolic function to determine the synchrotron peak frequency. We checked the dependence between luminosities at several frequency bands and synchrotron peak frequency to test the blazar sequence scenario, which states that the source luminosity depends on the location of the synchrotron peak. We also calculated broad band spectral indices and plotted them against each other and the peak frequency. The range of peak frequencies in our study was considerably extended compared to previous studies. There were 22 objects for which log\nu_{peak}>19. The data shows that at 5 GHz, 37 GHz and 5500 A there is negative correlation between luminosity and nu_{peak}. There is no significant correlation between source luminosity at synchrotron peak and peak frequency. Several low radio luminosity-low energy peaked BL Lacs were found. The negative correlation between broad band spectral indices and nu_{peak} is also significant, although there is substantial scatter. Therefore we find that neither alpha_{rx} nor alpha_{ro} can be used to determine the synchrotron peak of BL Lacs. On the grounds of our results we conclude that the blazar sequence scenario is not valid. In all our results the BL Lac population is continuous with no hint of the bimodality of the first BL Lac samples.

We report the discovery of an extremely curved jet in the high frequency peaking GPS quasar PKS 2136+141. Our multi-frequency VLBA images show a bending jet making a turn-around of 210 degrees on the plane of the sky, which is, to our knowledge, the largest ever observed change in a structural position angle of an extragalactic jet. Images taken at six different frequencies, from 2 to 43 GHz, beautifully reveal a spiral-like trajectory. We discuss possibilities to constrain the 3-D geometry of the source and suggest that it could be used as a testbed for models describing the bending of the relativistic jets.

We report first results of a multifrequency campaign from radio to hard X-ray energies of the prominent gamma-ray blazar 3C 279 during the first year of the INTEGRAL mission. The variable blazar was found at a low activity level, but was detected by all participating instruments. Subsequently a multifrequency spectrum could be compiled. The individual measurements as well as the compiled multifrequency spectrum are presented. In addition, this 2003 broadband spectrum is compared to one measured in 1999 during a high activity period of 3C 279.

We present results of long-term high-frequency VLBI monitoring of the relativistic jet in 3C279, consisting of 18 epochs at 22 GHz from 1991 to 1997 and 10 epochs at 43 GHz from 1995 to 1997. Three major results of this study are: apparent speeds measured for six superluminal components range from 4.8 to 7.5 c (H_{0}=70 km s^{-1} Mpc^{-1}, q_{0}=0.1), variations in the total radio flux are due primarily to changes in the VLBI core flux, and the uniform-sphere brightness temperature of the VLBI core is about 1 x 10^{13} K at 22 GHz after 1995, one of the highest direct estimates of a brightness temperature. If the variability brightness temperature measured for 3C279 by Lahteenmaki & Valtaoja is an actual value and not a lower limit, then the rest-frame brightness temperature of 3C279 is quite high and limited by inverse Compton effects rather than equipartition. The parsec-scale morphology of 3C279 consists of a bright, compact VLBI core, a jet component (C4) that moved from about 2 mas to about 3.5 mas from the core during the course of our monitoring, and an inner jet that extends from the core to a stationary component, C5, at about 1 mas from the core. Components in the inner jet are relatively short-lived, and fade by the time they reach about 1 mas from the core. The components have different speeds and position angles from each other, but these differences do not match the differences predicted by the precession model of Abraham & Carrara. Although VLBI components were born about six months prior to each of the two observed gamma-ray high states, the sparseness of the gamma-ray data prevents a statistical analysis of possible correlations.

We present multiwaveband observations of a well selected sample of 28 quasars and two radio galaxies with flat radio spectra and strong millimeter wave emission (referred to here as FSRQ's). The observations include multifrequency VLBI measurements, X-ray observations with ROSAT and submillimeter observations with the JCMT. Particularly interesting among many findings is a correlation between the X-ray to millimeter spectral index and fraction of flux density contained in the VLBI core. This tendency toward higher X-ray fluxes from sources with stronger jet emission implies that the knots in the jet are the prominent source of X-rays.

We have conducted a multiwavelength study of several radio sources within or near the error boxes of EGRET unidentified sources at mid to high Galactic latitude, under the hypothesis that the radio sources are blazars and are thus the best identification candidates for the EGRET objects. We show that one of these radio sources, PMN 0850-12, has a flux of 1.5 Jy at 22 GHz and a nearly flat spectrum up to 230 GHz and is thus very likely to be the correct identification for for the EGRET source 2EG J0852-1237.

We report results from five years of monitoring of the Einstein Cross (Q 2237+0305) with the Nordic Optical Telescope. The photometry, mainly in the R and I bands, has been performed by a PSF fitting and `cleaning' procedure, in which the four image components as well as the lens galaxy and its nucleus are iteratively removed. The resulting lightcurves exhibit several microlensing features; one event may have a timescale as short as 14 days. Variations of timescales of several years are found in all four images. This becomes even more convincing when our data are combined with data published for 1986-89. No clear high amplification event was observed during the period. A brightening of all four components during 1994 is interpreted as intrinsic variation.