We collect data at well‐sampled frequencies from the radio to the γ‐ray range for the following three complete samples of blazars: the Slew survey, the 1‐Jy samples of BL Lacs and the 2‐Jy sample of flat‐spectrum radio‐loud quasars (FSRQs). The fraction of objects detected in γ‐rays (E ≳ 100 MeV) is ∼ 17, 26 and 40 per cent in the three samples respectively. Except for the Slew survey sample, γ‐ray detected sources do not differ either from other sources in each sample, or from all the γ‐ray detected sources, in terms of the distributions of redshift, radio and X‐ray luminosities or of the broad‐band spectral indices (radio to optical and radio to X‐ray). We compute average spectral energy distributions (SEDs) from radio to γ‐rays for each complete sample and for groups of blazars binned according to radio luminosity, irrespective of the original classification as BL Lac or FSRQ. The resulting SEDs show a remarkable continuity in that (i) the first peak occurs in different frequency ranges for different samples/luminosity classes, with most luminous sources peaking at lower frequencies; (ii) the peak frequency of the γ‐ray component correlates with the peak frequency of the lower energy one; (iii) the luminosity ratio between the high and low frequency components increases with bolometric luminosity. The continuity of properties among different classes of sources and the systematic trends of the SEDs as a function of luminosity favour a unified view of the blazar phenomenon: a single parameter, related to luminosity, seems to govern the physical properties and radiation mechanisms in the relativistic jets present in BL Lac objects as well as in FSRQs. The general implications of this unified scheme are discussed while a detailed theoretical analysis, based on fitting continuum models to the individual spectra of most γ‐ray blazars, is presented in a separate paper.

The phenomenology of gamma-ray bright blazars can be accounted for by a sequence in the source power and intensity of the diffuse radiation field surrounding the relativistic jet. Correspondingly, the equilibrium particle distribution peaks at different energies. This leads to a trend in the observed properties: an increase of the observed power corresponds to: 1) a decrease in the frequencies of the synchrotron and inverse Compton peaks; 2) an increase in the ratio of the powers of the high and low energy spectral components. Objects along this sequence would be observationally classified respectively as high frequency BL Lac objects, low frequency BL Lac objects, highly polarized quasars and lowly polarized quasars. The proposed scheme is based on the correlations among the physical parameters derived in the present paper by applying to 51 gamma ray loud blazars two of the most accepted scenarios for the broad band emission of blazars, namely the synchrotron self--Compton and external Compton models, and explains the observational trends presented by Fossati et al. (1998) in a companion paper, dealing with the spectral energy distributions of all blazars. This gives us confidence that our scheme applies to all blazars as a class.Comment: 25 pages, 11 figures, uses mn.sty and psfig.tex. Accepted for publication in MNRA

Following the detection of strong TeV γ-ray flares from the BL Lac object 1ES 1959+650 with the Whipple 10 m Cherenkov telescope on May 16 and 17, 2002, we performed intensive Target of Opportunity (ToO) radio, optical, X-ray and TeV γ-ray observations from -2correlation properties. Although the X-ray and γ-ray fluxes seemed to be correlated in general, we found an "orphan" γ-ray flare that was not accompanied by an X-ray flare. While we detected optical flux variability with the Boltwood and Abastumani observatories, the data did not give evidence for a correlation between the optical flux variability with the observed X-ray and γ-ray flares. Within statistical errors of about 0.03 Jy at 14.5 GHz and 0.05 Jy at 4.8 GHz, the radio fluxes measured with the University of Michigan Radio Astrophysical Observatory (UMRAO) stayed constant throughout the campaign; the mean values agreed well with the values measured on May 7 and June 7, 2002 at 4.9 GHz and 15 GHz with the Very Large Array (VLA), and, at 4.8 GHz with archival flux measurements. After describing in detail the radio, optical, X-ray and γ-ray light curves and Spectral Energy Distributions (SEDs) we present initial modeling of the SED with a simple Synchrotron Self-Compton (SSC) model. With the addition of another TeV blazar with good broadband data, we consider the set of all TeV blazars to begin to look for a connection of the jet properties to the properties of the central accreting black hole thought to drive the jet. Remarkably, the temporal and spectral X-ray and γ-ray emission characteristics of TeV blazars are very similar, even though the masses estimates of their central black holes differ by up to one order of magnitude.

We present multiwavelength data for twelve blazars observed from 2008-2010 as part of an ongoing optical-infrared photometric monitoring project. Sources were selected to be bright, southern (δ < 20 • ) blazars observed by the Fermi Gamma-Ray Space Telescope, with daily and weekly gamma-ray fluxes made available from the start of the Fermi mission. Light curves are presented for the twelve blazars in BVRJK (0.4 to 2.2 µm) at near-daily cadence. We find that optical and infrared fluxes are well correlated in all sources, with no measured lag between bands. Gamma-ray bright flat spectrum radio quasars (FSRQs) in our sample have optical/infrared emission correlated with gamma-rays consistent with inverse Compton-scattering leptonic models for GeV emission. In FSRQs, the variability amplitude decreases towards optical/IR wavelengths, consistent with the presence of a thermal emission component from the accretion disk varying on significantly longer timescales than the jet synchrotron emission. In BL Lac objects, variability is mainly constant across wavelengths, consistent with a weak or radiatively inefficient disk. Also consistent with this picture, FSRQs have redder optical-infrared colors when they are brighter, while BL Lac objects show no such trend. Several objects show complicated color-magnitude behavior: AO 0235+164 appears in two different states depending on whether it is gammaray bright or not. OJ 287 and 3C 279 show some hysteresis tracks in their -2color-magnitude diagrams. Individual flares may be achromatic or otherwise depart from the trend, suggesting different jet components becoming important at different times. We present a time-dependent spectral energy distribution of the bright FSRQ 3C 454.3 during December 2009 flare, which is well fit by an external Compton model in the bright state, although day to day changes in the course of the flare pose challenges to a simple one-zone model. All data from the SMARTS monitoring program are publicly available on our website.

Abstract.We have considered all blazars observed in the X-ray band and for which the slope of the X-ray spectrum is available. We have collected 421 spectra of 268 blazars, including 12 archival unpublished ASCA spectra of 7 blazars whose analysis is presented here. The X-ray spectra of blazars show trends as a function of their power, confirming that the blazar overall energy distribution can be parameterized on the basis of one parameter only, i.e. the bolometric luminosity. This is confirmed by the relatively new hard (2-10 keV) X-ray data. Our results confirm the idea that in low power objects the X-ray emission mechanism is the synchrotron process, dominating both the soft and the hard X-ray emissions. Low energy peaked BL Lac objects are intermediate, often showing harder spectra in the hard X-ray band, suggesting that the synchrotron process dominates in the soft band, with the inverse Compton process dominating at high energies. The most powerful objects have X-ray spectra that are flat both in the soft and in the hard band, consistent with a dominating inverse Compton component.

We present a detailed analysis of week-long simultaneous observations of the blazar Mrk 421 at 2Y60 keV X-rays (RXTE ) and TeV -rays ( Whipple and HEGRA) in 2001. Accompanying optical monitoring was performed with the Mt. Hopkins 48 inch telescope. The unprecedented quality of this data set enables us to establish the existence of the correlation between the TeV and X-ray luminosities, and also to start unveiling some of its characteristics, in particular its energy dependence and time variability. The source shows strong variations in both X-ray and -ray bands, which are highly correlated. No evidence of an X-ray/-ray interband lag is found on the full week data set, with P 3 ks. A detailed analysis of the March 19 flare, however, reveals that data are not consistent with the peak of the outburst in the 2Y4 keV X-ray and TeV band being simultaneous. We estimate a 2:1 AE 0:7 ks TeV lag. The amplitudes of the X-ray and -ray variations are also highly correlated, and the TeV luminosity increases more than linearly with respect to the X-ray one. The high degree of correlation lends further support to the standard model in which a unique electron population produces the X-rays by synchrotron radiation and the -ray component by inverse Compton scattering. However, the finding that for the individual best observed flares the -ray flux scales approximately quadratically with respect to the X-ray flux poses a serious challenge to emission models for TeV blazars, as it requires rather special conditions and/or fine tuning of the temporal evolution of the physical parameters of the emission region. We briefly discuss the astrophysical consequences of these new findings in the context of the competing models for the jet emission in blazars.

Abstract. We performed an observational program with the X-ray satellite BeppoSAX, to study objects with extreme synchrotron peak frequencies (ν peak > 1 keV). Of the seven sources observed, five revealed or confirmed their extreme nature. Four showed peak frequencies in the range 1-5 keV, while one (1ES 1426+428) displayed a flat power law spectrum (αx = 0.92 ± 0.04) which locates its synchrotron peak at or above 100 keV. This is the third source of this type ever found, after Mkn 501 and 1ES 2344+514. In the context of the whole blazar class, the broad band properties of these objects confirm the scenario of a synchrotron peak smoothly spanning the IR -X-ray range, which explains the multi-frequency properties of the blazar class. Our data also confirm the large ν peak variability which seems to characterize this class of sources, compared with lower ν peak objects. Given the high synchrotron peak energies, which flag the presence of high energy electrons, these extreme BL Lacs are also good candidates for TeV emission, and therefore good probes of the IR background.

Mkn 421 was repeatedly observed with BeppoSAX in 1997-1998. The source showed a very rich phenomenology, with remarkable spectral variability. This is the second of two papers presenting the results of a thorough temporal and spectral analysis of all the data available to us, focusing in particular on the flare of April 1998, which was simultaneously observed also at TeV energies. The spectral analysis and correlations are presented in this paper, while the data reduction and timing analysis are the content of the companion paper. The spectral evolution during the flare has been followed over few ks intervals, allowing us to detect for the first time the peak of the synchrotron component shifting to higher energies during the rising phase, and then receding. This spectral analysis nicely confirms the delay of the flare at the higher energies, which in Paper I we quantified as a hard lag of a few ks. Furthermore, at the highest energies, evidence is found of variations of the inverse Compton component. The spectral and temporal information obtained challenge the simplest models currently adopted for the (synchrotron) emission and most importantly provide clues on the particle acceleration process. A scenario accounting for all the observational constraints is discussed, where electrons are injected at progressively higher energies during the development of the flare, and the achromatic decay is ascribed to the source light crossing time exceeding the particle cooling timescales.