We present the results of nine years of the blazar observing programme at the RATAN-600 radio telescope (2005-2014). The data were obtained at six frequency bands (1.1, 2.3, 4.8, 7.7, 11.2, 21.7 GHz) for 290 blazars, mostly BL Lacs. In addition, we used data at 37 GHz obtained quasi-simultaneously with the Metsahovi radio observatory for some sources. The sample includes blazars of three types: high-synchrotron peaked (HSP), low-synchrotron peaked (LSP), and intermediate-synchrotron peaked (ISP). We present several epochs of flux density measurements, simultaneous radio spectra, spectral indices and properties of their variability. The analysis of the radio properties of different classes of blazars showed that LSP and HSP BL Lac blazars are quite different objects on average. LSPs have higher flux densities, flatter spectra and their variability increases as higher frequencies are considered. On the other hand, HSPs are very faint in radio domain, tend to have steep low frequency spectra, and they are less variable than LSPs at all frequencies. Another result is spectral flattening above 7.7 GHz detected in HSPs, while an average LSP spectrum typically remains flat at both the low and high frequency ranges we considered.

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.

Mrk 501 is one of the brightest blazars at TeV energies and has been extensively studied since its first VHE detection in 1996. Our goal is to characterize in detail the source gamma-ray emission, together with the radio-to-X-ray emission, during the non-flaring (low) activity, which is less often studied than the occasional flaring (high) activity. We organized a multiwavelength (MW) campaign on Mrk 501 between March and May 2008. This multi-instrument effort included the most sensitive VHE gamma-ray instruments in the northern hemisphere, namely the imaging atmospheric Cherenkov telescopes MAGIC and VERITAS, as well as Swift, RXTE, the F-GAMMA, GASP-WEBT, and other collaborations and instruments. Mrk 501 was found to be in a low state of activity during the campaign, with a VHE flux in the range of 10%-20% of the Crab nebula flux. Nevertheless, significant flux variations were detected with various instruments, with a trend of increasing variability with energy. The broadband spectral energy distribution during the two different emission states of the campaign can be adequately described within the homogeneous one-zone synchrotron self-Compton model, with the (slightly) higher state described by an increase in the electron number density. This agrees with previous studies of the broadband emission of this source during flaring and non-flaring states. We report for the first time a tentative X-ray-to-VHE correlation during a low VHE activity. Although marginally significant, this positive correlation between X-ray and VHE, which has been reported many times during flaring activity, suggests that the mechanisms that dominate the X-ray/VHE emission during non-flaring-activity are not substantially different from those that are responsible for the emission during flaring activity.

We present a new catalogue of the RATAN-600 multi-frequency measurements for BL Lac objects. The purpose of this catalogue is to compile the BL Lac multi-frequency data that is acquired with the RATAN-600 simultaneously at several frequencies. The BL Lac objects emit a strongly variable and polarized non-thermal radiation across the entire electromagnetic spectrum from radio to gamma-rays and represent about one percent of known AGNs. They belong to the blazar population and differ from other blazars featureless optical spectrum, which sometimes have absorption lines, or have weak and narrow emission lines. One of the most effective ways of studying the physics of BL Lacs is the use of simultaneous multi-frequency data. The multi-frequency broadband radio spectrum was obtained simultaneously with an accuracy of up to 1-2 minutes for four to six frequencies: 1.1, 2.3, 4.8, 7.7, 11.2, and 21.7 GHz. The catalogue is based on the RATAN-600 observations and on the data from: equatorial coordinate and redshift, R-band magnitude, synchrotron peak frequency, SED classes and object type, literature. The present version of the catalogue contains RATAN-600 flux densities measurements over nine years (2006-2014), radio spectra at different epochs, and their parameters for more than 300 BL Lacs objects and candidates. The BL Lacs list is constantly updated with new observational data of RATAN-600.

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.

The discovery of rapidly variable Very High Energy (VHE; E > 100 GeV) gamma-ray emission from 4C +21.35 (PKS 1222+216) by MAGIC on 2010 June 17, triggered by the high activity detected by the Fermi Large Area Telescope (LAT) in high energy (HE; E > 100 MeV) gamma-rays, poses intriguing questions on the location of the gamma-ray emitting region in this flat spectrum radio quasar. We present multifrequency data of 4C +21.35 collected from centimeter to VHE during 2010 to investigate the properties of this source and discuss a possible emission model. The first hint of detection at VHE was observed by MAGIC on 2010 May 3, soon after a gamma-ray flare detected by Fermi-LAT that peaked on April 29. The same emission mechanism may therefore be responsible for both the HE and VHE emission during the 2010 flaring episodes. Two optical peaks were detected on 2010 April 20 and June 30, close in time but not simultaneous with the two gamma-ray peaks, while no clear connection was observed between the X-ray an gamma-ray emission. An increasing flux density was observed in radio and mm bands from the beginning of 2009, in accordance with the increasing gamma-ray activity observed by Fermi-LAT, and peaking on 2011 January 27 in the mm regime (230 GHz). We model the spectral energy distributions (SEDs) of 4C +21.35 for the two periods of the VHE detection and a quiescent state, using a one-zone model with the emission coming from a very compact region outside the broad line region. The three SEDs can be fit with a combination of synchrotron self-Compton and external Compton emission of seed photons from a dust torus, changing only the electron distribution parameters between the epochs. The fit of the optical/UV part of the spectrum for 2010 April 29 seems to favor an inner disk radius of <6 gravitational radii, as one would expect from a prograde-rotating Kerr black hole.

We present the results of a series of radio, optical, X-ray and gamma-ray observations of the BL Lac object S50716+714 carried out between April 2007 and January 2011. The multi-frequency observations were obtained using several ground and space based facilities. The intense optical monitoring of the source reveals faster repetitive variations superimposed on a long-term variability trend at a time scale of ~350 days. Episodes of fast variability recur on time scales of ~ 60-70 days. The intense and simultaneous activity at optical and gamma-ray frequencies favors the SSC mechanism for the production of the high-energy emission. Two major low-peaking radio flares were observed during this high optical/gamma-ray activity period. The radio flares are characterized by a rising and a decaying stage and are in agreement with the formation of a shock and its evolution. We found that the evolution of the radio flares requires a geometrical variation in addition to intrinsic variations of the source. Different estimates yield a robust and self-consistent lower limits of \delta > 20 and equipartition magnetic field B_eq > 0.36 G. Causality arguments constrain the size of emission region \theta < 0.004 mas. We found a significant correlation between flux variations at radio frequencies with those at optical and gamma-rays. The optical/GeV flux variations lead the radio variability by ~65 days. The longer time delays between low-peaking radio outbursts and optical flares imply that optical flares are the precursors of radio ones. An orphan X-ray flare challenges the simple, one-zone emission models, rendering them too simple. Here we also describe the spectral energy distribution modeling of the source from simultaneous data taken through different activity periods.

We report on the detection of a very rapid TeV gamma-ray flare from BL Lacertae on 2011 June 28 with the Very Energetic Radiation Imaging Telescope Array System (VERITAS). The flaring activity was observed during a 34.6-minute exposure, when the integral flux above 200 GeV reached $(3.4\pm0.6) \times 10^{-6} \;\text{photons}\;\text{m}^{-2}\text{s}^{-1}$, roughly 125% of the Crab Nebula flux measured by VERITAS. The light curve indicates that the observations missed the rising phase of the flare but covered a significant portion of the decaying phase. The exponential decay time was determined to be $13\pm4$ minutes, making it one of the most rapid gamma-ray flares seen from a TeV blazar. The gamma-ray spectrum of BL Lacertae during the flare was soft, with a photon index of $3.6\pm 0.4$, which is in agreement with the measurement made previously by MAGIC in a lower flaring state. Contemporaneous radio observations of the source with the Very Long Baseline Array (VLBA) revealed the emergence of a new, superluminal component from the core around the time of the TeV gamma-ray flare, accompanied by changes in the optical polarization angle. Changes in flux also appear to have occurred at optical, UV, and GeV gamma-ray wavelengths at the time of the flare, although they are difficult to quantify precisely due to sparse coverage. A strong flare was seen at radio wavelengths roughly four months later, which might be related to the gamma-ray flaring activities. We discuss the implications of these multiwavelength results.

We present simultaneous Planck, Swift, Fermi, and ground-based data for 105 blazars belonging to three samples with flux limits in the soft X-ray, hard X-ray, and gamma-ray bands. Our unique data set has allowed us to demonstrate that the selection method strongly influences the results, producing biases that cannot be ignored. Almost all the BL Lac objects have been detected by Fermi-LAT, whereas ~40% of the flat-spectrum radio quasars (FSRQs) in the radio, soft X-ray, and hard X-ray selected samples are still below the gamma-ray detection limit even after integrating 27 months of Fermi-LAT data. The radio to sub-mm spectral slope of blazars is quite flat up to ~70GHz, above which it steepens to <\alpha>~-0.65. BL Lacs have significantly flatter spectra than FSRQs at higher frequencies. The distribution of the rest-frame synchrotron peak frequency (\nupS) in the SED of FSRQs is the same in all the blazar samples with <\nupS>=10^13.1 Hz, while the mean inverse-Compton peak frequency, <\nupIC>, ranges from 10^21 to 10^22 Hz. The distributions of \nupS and of \nupIC of BL Lacs are much broader and are shifted to higher energies than those of FSRQs and strongly depend on the selection method. The Compton dominance of blazars ranges from ~0.2 to ~100, with only FSRQs reaching values >3. Its distribution is broad and depends strongly on the selection method, with gamma-ray selected blazars peaking at ~7 or more, and radio-selected blazars at values ~1, thus implying that the assumption that the blazar power is dominated by high-energy emission is a selection effect. Simple SSC models cannot explain the SEDs of most of the gamma-ray detected blazars in all samples. The SED of the blazars that were not detected by Fermi-LAT may instead be consistent with SSC emission. Our data challenge the correlation between bolometric luminosity and \nupS predicted by the blazar sequence.

Planck's all sky surveys at 30-857 GHz provide an unprecedented opportunity to follow the radio spectra of a large sample of extragalactic sources to frequencies 2-20 times higher than allowed by past, large area, ground-based surveys. We combine the results of the Planck Early Release Compact Source Catalog (ERCSC) with quasi-simultaneous ground-based observations, as well as archival data, at frequencies below or overlapping Planck frequency bands, to validate the astrometry and photometry of the ERCSC radio sources and study the spectral features shown in this new frequency window opened by Planck. The ERCSC source positions and flux density scales are found to be consistent with the ground-based observations. We present and discuss the spectral energy distributions (SEDs) of a sample of "extreme" radio sources to illustrate the richness of the ERCSC for the study of extragalactic radio sources. Variability is found to play a role in the unusual spectral features of some of these sources.

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 report on TeV gamma-ray observations of the blazar Mrk 421 (redshift of 0.031) with the VERITAS observatory and the Whipple 10m Cherenkov telescope. The excellent sensitivity of VERITAS allowed us to sample the TeV gamma-ray fluxes and energy spectra with unprecedented accuracy where Mrk 421 was detected in each of the pointings. A total of 47.3 hrs of VERITAS and 96 hrs of Whipple 10m data were acquired between January 2006 and June 2008. We present the results of a study of the TeV gamma-ray energy spectra as a function of time, and for different flux levels. On May 2nd and 3rd, 2008, bright TeV gamma-ray flares were detected with fluxes reaching the level of 10 Crab. The TeV gamma-ray data were complemented with radio, optical, and X-ray observations, with flux variability found in all bands except for the radio waveband. The combination of the RXTE and Swift X-ray data reveal spectral hardening with increasing flux levels, often correlated with an increase of the source activity in TeV gamma-rays. Contemporaneous spectral energy distributions were generated for 18 nights, each of which are reasonably described by a one-zone SSC model.

We present observations of a major outburst at centimeter, millimeter, optical, X-ray, and gamma-ray wavelengths of the BL Lacertae object AO 0235+164. We analyze the timing of multi-waveband variations in the flux and linear polarization, as well as changes in Very Long Baseline Array (VLBA) images at 7mm with 0.15 milliarcsecond resolution. The association of the events at different wavebands is confirmed at high statistical significance by probability arguments and Monte-Carlo simulations. A series of sharp peaks in optical linear polarization, as well as a pronounced maximum in the 7 mm polarization of a superluminal jet knot, indicate rapid fluctuations in the degree of ordering of the magnetic field. These results lead us to conclude that the outburst occurred in the jet both in the quasi-stationary "core" and in the superluminal knot, both parsecs downstream of the supermassive black hole. We interpret the outburst as a consequence of the propagation of a disturbance, elongated along the line of sight by light-travel time delays, that passes through a standing recollimation shock in the core and propagates down the jet to create the superluminal knot. The multi-wavelength light curves vary together on long time-scales (months/years), but the correspondence is poorer on shorter time-scales. This, as well as the variability of the polarization and the dual location of the outburst, agrees with the expectations of a multi-zone emission model in which turbulence plays a major role in modulating the synchrotron and inverse Compton fluxes.

We combine time-dependent multi-waveband flux and linear polarization observations with sub-milliarcsecond-scale polarimetric images at lambda=7mm of the BL Lacertae-type blazar OJ287 to locate the gamma-ray emission in prominent flares in the jet of the source >14pc from the central engine. We demonstrate a highly significant correlation between the strongest gamma-ray and millimeter-wave flares through Monte-Carlo simulations. The two reported gamma-ray peaks occurred near the beginning of two major mm-wave outbursts, each of which is associated with a linear polarization maximum at millimeter wavelengths. Our Very Long Baseline Array observations indicate that the two mm-wave flares originated in the second of two features in the jet that are separated by >14 pc. The simultaneity of the peak of the higher-amplitude gamma-ray flare and the maximum in polarization of the second jet feature implies that the gamma-ray and mm-wave flares are co-spatial and occur >14 pc from the central engine. We also associate two optical flares, accompanied by sharp polarization peaks, with the two gamma-ray events. The multi-waveband behavior is most easily explained if the gamma-rays arise from synchrotron self-Compton scattering of optical photons from the flares. We propose that flares are triggered by interaction of moving plasma blobs with a standing shock. The gamma-ray and optical emission is quenched by inverse Compton losses as synchrotron photons from the newly shocked plasma cross the emission region. The mm-wave polarization is high at the onset of a flare, but decreases as the electrons emitting at these wavelengths penetrate less polarized regions.

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.

Narrow-Line Seyfert 1 (NLS1) class of active galactic nuclei (AGNs) is generally radio-quiet, but a small percent of them are radio-loud. The recent discovery by Fermi/LAT of high-energy gamma-ray emission from 4 NLS1s proved the existence of relativistic jets in these systems. It is therefore important to study this new class of gamma-ray emitting AGNs. Here we report preliminary results about the observations of the July 2010 gamma-ray outburst of PMN J0948+0022, when the source flux exceeded for the first time 10^-6 ph cm^-2 s^-1 (E > 100 MeV).

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.

The well-studied VHE (E >100 GeV) blazar Mrk 501 was observed between March and May 2008 as part of an extensive multiwavelength observation campaign including radio, optical, X-ray and VHE gamma-ray instruments. Mrk 501 was in a low state of activity during the campaign, with a low VHE flux of about 20% the Crab Nebula flux. Nevertheless, significant flux variations could be observed in X-rays as well as $\gamma$-rays. Overall Mrk 501 showed increased variability when going from radio to gamma-ray energies. The broadband spectral energy distribution during the two different emission states of the campaign was well described by a homogeneous one-zone synchrotron self-Compton model. The high emission state was satisfactorily modeled by increasing the amount of high energy electrons with respect to the low emission state. This parameterization is consistent with the energy-dependent variability trend observed during the campaign.

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.

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