We have collected a large amount of multifrequency data for objects in the Metsähovi Radio Observatory BL Lacertae sample and computed their spectral energy distributions (SED) in the log ν -log νF -representation. 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, ν peak . 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 ν peak . The range of ν peak in our study was considerably extended compared to previous studies. There were 22 objects for which log ν peak >19. The data shows that at 5 GHz, 37 GHz, and 5500 Å, there is negative correlation between luminosity and ν peak , whereas in X-rays the correlation turns slightly positive. There is no significant correlation between source luminosity at synchrotron peak and ν peak . Several low radio luminosity-low energy peaked BL Lacs were found. The negative correlation between broad band spectral indices and ν peak is also significant, although there is substantial scatter. Therefore we find that neither α rx nor α 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.

Context. 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. Aims. 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. Methods. 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. Results. 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. Conclusions. 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.

We report on TeV γ -ray observations of the blazar Mrk 421 (redshift of 0.031) with the VERITAS observatory and the Whipple 10 m Cherenkov telescope. The excellent sensitivity of VERITAS allowed us to sample the TeV γ -ray fluxes and energy spectra with unprecedented accuracy where Mrk 421 was detected in each of the pointings. A total of 47.3 hr of VERITAS and 96 hr of Whipple 10 m data were acquired between 2006 January and 2008 June. We present the results of a study of the TeV γ -ray energy spectra as a function of time and for different flux levels. On 2008 May 2 and 3, bright TeV γ -ray flares were detected with fluxes reaching the level of 10 Crab. The TeV γ -ray data were complemented with radio, optical, and X-ray observations, with flux variability found in all bands except for the radio wave band. The combination of the Rossi X-ray Timing Explorer and Swift X-ray data reveal spectral hardening with increasing flux levels, often correlated with an increase of the source activity in TeV 1 The Astrophysical Journal, 738:25 (19pp), 2011 September 1 Acciari et al.γ -rays. Contemporaneous spectral energy distributions were generated for 18 nights, each of which are reasonably described by a one-zone synchrotron self-Compton model.

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 γ-ray bands, with additional 5 GHz flux-density limits to ensure a good probability of a Planck detection. We compare our results to those of a companion paper presenting simultaneous Planck and multi-frequency observations of 104 radio-loud northern active galactic nuclei selected at radio frequencies. While we confirm several previous results, our unique data set allows 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 the Fermi Large Area Telescope (LAT), whereas 30% to 40% of the flat-spectrum radio quasars (FSRQs) in the radio, soft X-ray, and hard X-ray selected samples are still below the γ-ray detection limit even after integrating 27 months of Fermi-LAT data. The radio to sub-millimetre spectral slope of blazars is quite flat, with α ∼ 0 up to about 70 GHz, above which it steepens to α ∼ −0.65. The BL Lacs have significantly flatter spectra than FSRQs at higher frequencies. The distribution of the rest-frame synchrotron peak frequency (ν of BL Lacs are much broader and are shifted to higher energies than those of FSRQs; their shapes strongly depend on the selection method. The Compton dominance of blazars, defined as the ratio of the inverse Compton to synchrotron peak luminosities, ranges from less than 0.2 to nearly 100, with only FSRQs reaching values larger than about 3. Its distribution is broad and depends strongly on the selection method, with γ-ray selected blazars peaking at ∼7 or more, and radio-selected blazars at values close to 1, thus implying that the common assumption that the blazar power budget is largely dominated by high-energy emission is a selection effect. A comparison of our multi-frequency data with theoretical predictions shows that simple homogeneous SSC models cannot explain the simultaneous SEDs of most of the γ-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 ν S peak predicted by the blazar sequence.

We combine time-dependent multi-waveband flux and linear polarization observations with submilliarcsecond-scale polarimetric images at λ = 7 mm of the BL Lacertae-type blazar OJ287 to locate the γ-ray emission in prominent flares in the jet of the source > 14 pc from the central engine. We demonstrate a highly significant correlation between the strongest γ-ray and millimeter-wave flares through Monte-Carlo simulations. The two reported γ-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 γ-ray flare and the maximum in polarization of the second jet feature implies that the γ-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 γ-ray events. The multi-waveband behavior is most easily explained if the γ-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 γ-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.

Context. The quasar 3C 279 is among the most extreme blazars in terms of luminosity and variability of flux at all wavebands. Its variations in flux and polarization are quite complex and therefore require intensive monitoring observations at multiple wavebands to characterise and interpret the observed changes. Aims. In this paper, we present radio-to-optical data taken by the WEBT, supplemented by our VLBA and RXTE observations, of 3C 279. Our goal is to use this extensive database to draw inferences regarding the physics of the relativistic jet. Methods. We assemble multifrequency light curves with data from 30 ground-based observatories and the space-based instruments SWIFT (UVOT) and RXTE, along with linear polarization vs. time in the optical R band. In addition, we present a sequence of 22 images (with polarization vectors) at 43 GHz at resolution 0.15 milliarcsec, obtained with the VLBA. We analyse the light curves and polarization, as well as the spectral energy distributions at different epochs, corresponding to different brightness states. Results. We find that the IR-optical-UV continuum spectrum of the variable component corresponds to a power law with a constant slope of −1.6, while in the 2.4-10 keV X-ray band it varies in slope from −1.1 to −1.6. The steepest X-ray spectrum occurs at a flux minimum. During a decline in flux from maximum in late 2006, the optical and 43 GHz core polarization vectors rotate by ∼300• . Conclusions. The continuum spectrum agrees with steady injection of relativistic electrons with a power-law energy distribution of slope −3.2 that is steepened to −4.2 at high energies by radiative losses. The X-ray emission at flux minimum comes most likely from a new component that starts in an upstream section of the jet where inverse Compton scattering of seed photons from outside the jet is important. The rotation of the polarization vector implies that the jet contains a helical magnetic field that extends ∼20 pc past the 43 GHz core.

Aims. 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. Methods. 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ähovi 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. Results. 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 compare the γ-ray photon flux variability of northern blazars in the Fermi/LAT First Source Catalog with 37 GHz radio flux density curves from the Metsähovi quasar monitoring program. We find that the relationship between simultaneous millimeter (mm) flux density and γ-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 γ-ray emission in high states depend on the phase of the high frequency radio flare, with the brightest γ-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 γ-ray emission is about 70 days, which places the average location of the γ-ray production at or downstream of the radio core. We discuss alternative scenarios for the production of γ-rays at distances of parsecs along the length of the jet.