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 5GHz 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 AreaTelescope (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 70GHz, 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 (νpeak S) in the spectral energy distribution (SED) of FSRQs is the same in all the blazar samples with ⟨νpeak S⟩ = 1013.1 ± 0.1 Hz, while the mean inverse Compton peak frequency, ⟨νpeak IC⟩, ranges from 1021 to 1022 Hz. The distributions of νpeak S and νpeak IC 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 νpeak S predicted by the blazar sequence.

In order to determine the location of the gamma-ray emission site in blazars, we investigate the time-domain relationship between their radio and gamma-ray emission. Light curves for the brightest detected blazars from the first 3 yr of the mission of the Fermi Gamma-ray Space Telescope are cross-correlated with 4 yr of 15 GHz observations from the Owens Valley Radio Observatory 40 m monitoring programme. The large sample and long light-curve duration enable us to carry out a statistically robust analysis of the significance of the cross-correlations, which is investigated using Monte Carlo simulations including the uneven sampling and noise properties of the light curves. Modelling the light curves as red noise processes with power-law power spectral densities, we find that only one of 41 sources with high-quality data in both bands shows correlations with significance larger than 3σ (AO 0235+164), with only two more larger than even 2.25σ (PKS 1502+106 and B2 2308+34). Additionally, we find correlated variability in Mrk 421 when including a strong flare that occurred in 2012 July–September. These results demonstrate very clearly the difficulty of measuring statistically robust multiwavelength correlations and the care needed when comparing light curves even when many years of data are used. This should be a caution. In all four sources, the radio variations lag the gamma-ray variations, suggesting that the gamma-ray emission originates upstream of the radio emission. Continuous simultaneous monitoring over a longer time period is required to obtain high significance levels in cross-correlations between gamma-ray and radio variability in most blazars.

Polarization observations of the cosmic microwave background with the Cosmic Background Imager from September 2002 to May 2004 provide a significant detection of the E -mode polarization and reveal an angular power spectrum of polarized emission showing peaks and valleys that are shifted in phase by half a cycle relative to those of the total intensity spectrum. This key agreement between the phase of the observed polarization spectrum and that predicted on the basis of the total intensity spectrum provides support for the standard model of cosmology, in which dark matter and dark energy are the dominant constituents, the geometry is close to flat, and primordial density fluctuations are predominantly adiabatic with a matter power spectrum commensurate with inflationary cosmological models.

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 multiple components. Planck ERCSC data also suggest that the original accelerated electron energy spectrum could be much harder than commonly thought, with power-law indexaround 1.5 instead of the canonical 2.5. The implications of this are discussed for the acceleration mechanisms effective in blazar shocks. Furthermore in many cases the Planck data indicate that gamma-ray emission must originate in the same shocks that produce the radio emission.

Most large galaxies host supermassive black holes in their nuclei and are subject to mergers, which can produce a supermassive black hole binary (SMBHB), and hence periodic signatures due to orbital motion. We report unique periodic radio flux density variations in the blazar PKS 2131−021, which strongly suggest an SMBHB with an orbital separation of ∼0.001–0.01 pc. Our 45.1 yr radio light curve shows two epochs of strong sinusoidal variation with the same period and phase to within ≲2% and ∼10%, respectively, straddling a 20 yr period when this variation was absent. Our simulated light curves accurately reproduce the “red noise” of this object, and Lomb–Scargle, weighted wavelet Z-transform and least-squares sine-wave analyses demonstrate conclusively, at the 4.6σ significance level, that the periodicity in this object is not due to random fluctuations in flux density. The observed period translates to 2.082 ± 0.003 yr in the rest frame at the z = 1.285 redshift of PKS 2131−021. The periodic variation in PKS 2131−021 is remarkably sinusoidal. We present a model in which orbital motion, combined with the strong Doppler boosting of the approaching relativistic jet, produces a sine-wave modulation in the flux density that easily fits the observations. Given the rapidly developing field of gravitational-wave experiments with pulsar timing arrays, closer counterparts to PKS 2131−021 and searches using the techniques we have developed are strongly motivated. These results constitute a compelling demonstration that the phenomenology, not the theory, must provide the lead in this field.