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 $\ ...\nu 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, $\nu_{\rm 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 $\nu_{\rm peak}$. The range of $\nu_{\rm peak}$ in our study was considerably extended compared to previous studies. There were 22 objects for which log $\nu_{\rm peak}>$19. The data shows that at 5 GHz, 37 GHz, and 5500 ${\rm \AA}$, there is negative correlation between luminosity and $\nu_{\rm peak}$, whereas in X-rays the correlation turns slightly positive. There is no significant correlation between source luminosity at synchrotron peak and $\nu_{\rm peak}$. Several low radio luminosity-low energy peaked BL Lacs were found. The negative correlation between broad band spectral indices and $\nu_{\rm peak}$ is also significant, although there is substantial scatter. Therefore we find that neither $\alpha_{\rm rx}$ nor $\alpha_{\rm 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 address the highly debated issue of constraining the γ-ray emission region in blazars from cross-correlation analysis using discrete correlation function between radio and γ-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 γ-ray peaks lead the radio with time lags in the range +20 and +690 d, whereas in sources 1633+382 and 3C 345 we find the radio emission to lead the γ-rays by −15 and −40 d, respectively. Apart from the individual source study, we stacked the correlations of all sources and also those based on subsamples. The time lag from the stacked correlation is +80 d 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 γ-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 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.
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 present the results of a series of radio, optical, X-ray, and γ-ray observations of the BL Lac object S50716+714 carried out between April 2007 and January 2011. The multifrequency 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 on 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 γ-ray frequencies favors the synchrotron self-Compton mechanism for the production of the high-energy emission. Two major low-peaking radio flares were observed during this high optical/γ-ray activity period. The radio flares are characterized by a rising and a decaying stage and agrees 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 robust and self-consistent lower limits of δ ≥ 20 and equipartition magnetic field Beq ≥ 0.36 G. Causality arguments constrain the size of emission region θ ≤ 0.004 mas. We found a significant correlation between flux variations at radio frequencies with those at optical and γ-rays. Theoptical/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.
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.
Aims. 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 γ rays and represent about 1% 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. Methods. The multi-frequency broadband radio spectrum was obtained simultaneously with an accuracy of up to 1−2 min 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. Results. 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 of the catalogue for more than 300 BL Lacs objects and candidates. The BL Lacs list is constantly updated with new observational data of RATAN-600.
Aims. 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. Methods. 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 Metsähovi 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. Results. 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.
Aims. Since the CGRO operation in 1991–2000, one of the primary unresolved questions about the blazar γ-ray emission has been its possible correlation with the low-energy (in particular optical) ...emission. To help answer this problem, the Whole Earth Blazar Telescope (WEBT) consortium has organized the GLAST-AGILE Support Program (GASP) to provide the optical-to-radio monitoring data to be compared with the γ-ray detections by the AGILE and GLAST satellites. This new WEBT project started in early September 2007, just before a strong γ-ray detection of 0716+714 by AGILE. Methods. We present the GASP-WEBT optical and radio light curves of this blazar obtained in July–November 2007, about various AGILE pointings at the source. We construct NIR-to-UV spectral energy distributions (SEDs), by assembling GASP-WEBT data together with UV data from the Swift ToO observations of late October. Results. We observe a contemporaneous optical-radio outburst, which is a rare and interesting phenomenon in blazars. The shape of the SEDs during the outburst appears peculiarly wavy because of an optical excess and a UV drop-and-rise. The optical light curve is well sampled during the AGILE pointings, showing prominent and sharp flares. A future cross-correlation analysis of the optical and AGILE data will shed light on the expected relationship between these flares and the γ-ray events.