A unifying view of the spectral energy distributions of blazars

We collect data at well-sampled frequencies from the radio to the gamma-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 gamma-rays (E greater than or similar to 100 MeV) is similar to 17, 26 and 40 per cent in the three samples respectively. Except for the Slew survey sample, gamma-ray detected sources do not differ either from other sources in each sample, or from all the gamma-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 gamma-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 gamma-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 gamma-ray blazars, is presented in a separate paper.