Abstract

Recent polarimetric surveys of extragalactic radio sources (ERS) at frequencies  GHz are reviewed. By exploiting all the most relevant data we study the frequency dependence of polarization properties of ERS between 1.4 and 86 GHz. For flat-spectrum sources the median (mean) fractional polarization increases from 1.5% (2–2.5%) at 1.4 GHz to 2.5–3% (3–3.5%) at  GHz. Steep-spectrum sources are typically more polarized, especially at high frequencies where Faraday depolarization is less relevant. Current data suggest that at high radio frequencies  GHz the fractional polarization of ERS does not depend on total flux density and moderately increases with frequency. We estimate ERS number counts in polarization and the contribution of unresolved polarized ERS to angular power spectra. A first application is for the Planck satellite mission: we predict that only a dozen polarized ERS will be detected by the Planck LFI, and a few tens by the HFI. As for CMB power spectra, ERS should not be a strong contaminant to the CMB E-mode polarization at  GHz, but they can become a relevant constraint for the detection of the cosmological B-mode polarization if the tensor-to-scalar ratio is .

1. Introduction

The radiation emitted by extragalactic radio sources (ERS) at high radio frequencies, mainly synchrotron radiation from relativistic electrons in their jets and lobes, can be highly polarized, with an intrinsic degree of linear polarization as high as ~70–75% in homogeneous sources with an unidirectional magnetic field (see, e.g., [1, 2], for comprehensive discussions on the subject).

These have to be considered as maximum values, which could be detected only in the case of very homogeneous sources with highly aligned magnetic field lines, which are unlikely to be observed among ERS. Actually, a much lower degree of total linear polarization, , is commonly observed in ERS at cm or mm wavelengths (see e.g., [3, 4]), with only very few ERS showing a total fractional polarization, , as high as ~10% of the total flux density, .

Nevertheless, this low fractional polarization observed in ERS may constitute a problem for the detection of primordial polarization in the Cosmic Microwave Background (CMB), since ERS are the dominant polarized foreground at small angular scales and the CMB polarized signal is only a few percent of CMB temperature anisotropy. Therefore, CMB polarization studies need a careful determination of the polarized emission from foreground sources in general and of ERS in particular.

CMB polarized anisotropies can be decomposed into modes of even-parity (E-mode) and odd-parity (B-mode). As for the E-modes, generated by scalar perturbations in the primordial universe, they were first detected by [5] and all the following observations (e.g., [6, 7]) have confirmed measurements of CMB E-modes compatible with the concordance CDM model. The more elusive B-modes are generated by tensor metric perturbations, that is, “primordial” gravitational waves generated during inflation, and, according to predictions of inflationary models, with an amplitude directly proportional to the energy scale at which the inflation occurred [8, 9]. A detection of these “primordial” B-modes would provide the first real measure of the energy scale of inflation; thus, it would produce a real breakthrough in modern cosmology. But the detection of the B-modes CMB polarization—parameterized by the tensor-to-scalar ratio, —is still a great challenge. In fact, the Planck mission [10, 11] (detailed discussions of the Planck Low Frequency Instrument (LFI) and High Frequency Instrument (HFI) expected polarization capabilities and calibrations are given by [12, 13]) would be marginally able to detect tensor metric perturbations by the direct detection of the primordial CMB B-mode, but only in the case of very high values of the tensor-to-scalar ratio of primordial perturbations, 0.05–0.1 [14]. Otherwise, we will have to rely on proposed future space experiments, like [15], specifically designed to detect CMB polarization by virtue of a much higher sensitivity than current ones.

The above discussion clearly illustrates the importance of giving the best up-to-date estimate of the average/median fractional polarization as well as of the distributions of fractional polarization [16], , at least for each one of the different ERS populations [1719] which are providing relevant contributions to temperature anisotropies of the CMB.

Statistical studies of polarized emission from ERS have shown that at ~1.4–5 GHz, the fractional polarization, , increases with decreasing flux density [3, 16, 2022]. Different explanations have been proposed of this result, for example, a population change at fainter flux density [20, 22] or a changing fraction of radio-quiet AGN [3]. However, the cause of this increase in is still unknown. On the other hand, it has been shown by [16] that this increase of the degree of polarization with decreasing flux density is only observed, at cm wavelengths, for steep-spectrum ERS, that is, sources with if , and not for flat-spectrum ERS ().

Other very recent studies of the polarized emission in ERS have tried to analyze the dependence of the fractional polarization with luminosity, redshift, and the source environment. The current, still preliminary, results show no correlation between the fractional polarization and redshift, whereas a weak correlation is found between decreasing luminosity and increasing degree of polarization [23]. According to the analysis of highly polarized elliptical galaxies [24] no differences have been found in the source environments between low polarization and ultrahigh polarization sources. This result should indicate that a high polarization must be a result of intrinsic properties of ERS.

The outline of the paper is as follows: in Section 2 we briefly summarize the main processes giving rise to linear and circular polarization in ERS; in Section 3 we present current published data on the polarized emission from ERS; Section 4 is dedicated to discuss our current results on statistical properties of the polarized emission in ERS; in Section 5 we give a short review on the most recent cosmological evolution model for ERS; Section 6 presents our predictions on the contributions of ERS polarized radiation, given by unresolved sources, to the E- and B-modes; finally, in Section 7, we present our conclusions.

2. Polarized Emission from ERS

Observations as well as statistical studies of the fractional polarization of ERS are very interesting on their own, and not only because ERS constitute a major contaminant of the CMB polarized signal. In fact, they are an important tool in active galactic nuclei (AGN) research, as they give valuable information on the physical quantities which determine the characteristics of the synchrotron radiation emitted by relativistic electrons accelerated by homogeneous and/or random magnetic fields in AGN [1, 2].

On the one hand, precise measurements of the total synchrotron radiation emitted by a radio source give an estimate—under some assumption, that is, equipartition between the field and particle energy—of the total magnetic field strength. On the other hand, the degree of polarization of the radio wave provides information on the direction of the main magnetic field in the source environment and also on its degree of ordering. The magnetic field structure can, in turn, provide information on the relationship between the environments of the ERS and their properties.

As already noted, compact ERS typically show a degree of total linear polarization of a few percent of their radio total intensity. Therefore, magnetic fields in radio sources are believed to be highly inhomogeneous, or almost without ordering, although the observed nonvanishing linear polarization gives an indication of a certain degree of ordering of the field [25, 26].

The precise orientation of magnetic fields lines inside jets and lobes is still unknown, but theoretical arguments as well as observational evidence show that magnetic fields are indeed partially ordered. As for observations, these conclusions are based on measurements of the orientations of the linear polarization, revealing coherent structures across the images. From the theoretical point of view, an ordered magnetic field is expected when shocks compress an initially random field (with the field perpendicular to the jet axis) or when such initial fields are sheared to lie in a plane, with parallel to the jet axis [2730].

Circular polarization (CP), which is a common feature of quasars and BL Lacertae objects (or simply BL Lacs), commonly known as blazars, (blazar sources are jet-dominated extragalactic objects—observed within a small angle of the jet axis—in which the beamed component dominates the observed emission [31]), is preferentially generated near synchrotron self-absorbed jet cores and is detected in about 30%–50% of these sources [25, 32]. In these inhomogeneous, optically thick, synchrotron sources, the emission of the electron population at lower energies is hidden by self-absorption. These invisible electrons produce Faraday rotation and conversion, which is the most likely mechanism capable of creating the observed CP in blazar sources [26, 32]. In any case, the measured degrees of CP are generally well below the levels of linear polarization [33], and thus negligible, at least at GHz frequencies.

The change in the position angle of the linearly polarized radiation which passes through a magnetoionic medium, that is, the Faraday rotation, can be expressed by the rotation angle [26, 34] experienced by the polarization vector, where indicates the rotation measure. The rotation measure, , is the line-of-sight integral , being the electron density and the component of the magnetic field along the line-of-sight [35].

If the rotation depth is the same for all the emission volumes of the radio source, the net result is a rotation equal to of the direction of polarization, without any effect on the degree of polarization. To change the degree of polarization there must be a variation in depth, either along or transverse to the line-of-sight [26].

Faraday rotation is commonly observed towards extragalactic radio sources. It was first observed by [36], who found an RM of −60 rad/ towards the center of Centaurus A. More recently, observations of the rest-frame RM in AGN cores and jets reported values from a few hundreds to several thousands rad/m2 (e.g., [3739]). On the other hand, ERS dominated by emission from radio lobes show lower values of RMs, that is, from a few tens to hundreds of rad/m2 (e.g., [40]).

3. Polarization Data on ERS at cm and mm Wavelengths

3.1. Polarization Surveys at 1.4 GHz

A large-scale polarization catalogue of radio sources is provided by the NRAO VLA Sky Survey (NVSS) at 1.4 GHz [41]. This survey covers  sr of the sky with . The catalogue contains the flux density and the Stokes parameters and of almost discrete sources with  mJy. Extensive analyses of these data were carried out by [20] and [16]. By correlating NVSS sources with sources from the Green Bank 4.85 GHz (GB6) catalogue [42], statistical polarization properties have been derived for a subsample of  ERS with  mJy, divided into steep- and flat-spectrum sources. Steep-spectrum sources are found to increase the degree of polarization () with decreasing flux density: in fact, the median (mean) value increases from (%) for  mJy up to (%) for the faintest ERS of the sample, at  mJy. On the other hand, flat-spectrum sources show a lower fractional polarization (median ~1.3% and mean %) with no significant trend with flux density [16].

Two recent sub-mJy surveys at 1.4 GHz with polarization measurements, that is, the Dominion Radio Astrophysical Observatory (DRAO) Planck Deep Fields project [3, 23, 43] and the Australia Telescope Low-Brightness Survey (ATLBS) [22], seem to confirm the above result of a higher level of fractional polarization in steep-spectrum sources with fainter flux densities. However, there is no clear explanation for the origin of the anticorrelation between and flux density in this source population.

In Figure 1 we plot differential number counts at 1.4 GHz for the polarized intensity, computed from the surveys presented above. The continuous curve represents a fit to total number counts of AGNs from [44], whereas the dotted curve is obtained from the previous fit by assuming a constant fractional polarization of % for all ERS. This value is chosen in order to fit NVSS data; the fit is extremely good down to polarized intensities of  mJy. At fainter fluxes the predicted curve bends down whereas observational data keep a flatter shape, in agreement with a higher fractional polarization for ERS at fainter flux density levels.

3.2. Multiwavelength Samples of Polarized ERS
3.2.1. Linear Polarization in a Sample of the B3-VLA Survey

Polarization measurements of sources in the B3-VLA survey [45] were carried out by [46] at 10.5, 4.85, and 2.7 GHz using the Effelsberg 100 m telescope. Taking into account only sources with  mJy, a sample of 106 objects (out of 208) was defined with detected polarization. They found that flat-spectrum sources are significantly less polarized than steep-spectrum ones at 10.5 GHz. The latter ERS population shows a median fractional polarization that strongly depends on the frequency, from ~2% at 1.4 GHz to ~6 at 10.5 GHz, indicating a high Faraday depolarization at lower frequencies. Moreover, they notice that compact steep-spectrum sources exhibit much stronger depolarization than noncompact ones, and that sources showing larger linear size are more polarized. On the other hand, flat-spectrum sources are characterized by almost constant and low degrees of polarization (~2.5%) over the whole wavelength range considered.

3.2.2. The Australia Telescope 20 GHz (AT20G) Survey

The AT20G survey is a blind survey of the whole Southern sky at 20 GHz with follow-up observations at 4.8 and 8.6 GHz carried out with the Australia Telescope Compact Array (ATCA) [4]. The full source catalogue includes 5890 sources detected above a flux-density limit of 40 mJy. The AT20G catalogue is found to be 91% complete above 100 mJy and 79% above 50 mJy in regions south of declination . Polarization was detected at 20 GHz for 768 sources, 467 of which also have simultaneous polarization detections at 5 and/or 8 GHz (upper limits on the polarized flux density are provided for nondetections, as described in [4]). Taking into account upper limits, a median (mean) fractional polarization of 2.6% (2.7%) at 20 GHz is found for sources with  mJy, whereas an average 17% depolarization is observed at 5 GHz with respect to 20 GHz [47]. Compact sources with detected polarization are separated into flat- and steep-spectrum sources: the mean values of are 2.9% and 3.8%, respectively, thus confirming the higher fractional polarization in steep-spectrum sources.

3.2.3. VLA Polarization Measurements of WMAP Point Sources

Using the VLA, [48] carried out polarization measurements at 8.4, 22 and 43 GHz of a complete sample of ERS brighter than 1 Jy in the 5-year WMAP catalogue and with declinations north of . The sample consists of 203 objects: polarized emission was detected for 123, 169, and 167 sources at 8.4, 22, and 43 GHz, respectively, (at 8.4 GHz only a subset of 134 were observed) and 105 were detected at all the 3 frequencies. An accurate analysis of the statistical properties of the polarized intensity of the sample is done by [39], including an analysis of the correlations between the fractional polarization and the spectral indices. Here below, we summarize the main results obtained by them.(i)The distribution of the fractional polarization varies slightly as a function of the frequency. Including sources undetected in polarization, the mean fractional polarization is 3.0 and 3.5 percent at 8.4, 22, and 43 GHz.(ii)No correlation is found between the fractional polarization at 22 GHz and the intensity spectral indices, and (we remind readers that the sample, selected at 22 GHz, is dominated by flat-spectrum sources).(iii)There is a significant change in the polarization angle between 8.4 and 22 GHz. For 45 sources, the position angle satisfies the dependence, and intrinsic rotation measures in excess of 1000 were observed for a large number of them.(iv)Polarization of 71 sources was also measured at 86 GHz by [49]. The fractional polarization is typically higher at 86 GHz than at 43 GHz.

3.2.4. VLA Polarization Measurements in an ACT Survey Field

In [50], VLA observations in total intensity and polarization at 4.86, 8.46, 22.46, and 43.34 GHz are presented. The sample is selected from the AT20G survey and consists of 159 ERS, out of the almost 200 sources with flux density  mJy in a field of the Atacama Comsmology Telescope (ACT) survey. Polarized emission for about 60 sources was detected at all 4 frequencies, whereas the detections are 141, 146, 89, and 59 from low to high frequencies. Fractional polarization distributions are very similar at 5 and 8.5 GHz, whereas a trend of increasing polarization with increasing frequency is indicated at 22 and 43 GHz. In particular, a tail of strongly polarized (%) sources is observed at 43 GHz. Data at 22 and 43 GHz suggest also that the polarization fraction in steep-spectrum sources is significant higher than in flat-spectrum sources.

3.3. Other Samples of Polarized ERS at High Frequencies
(i) Reference [51] presented polarization observations of 250 (out of 258) southern sources in the complete 5-GHz 1-Jy sample of [52] by using the ATCA facilities at 18.5 GHz. Polarized flux densities were measured for 170 sources (114 flat-spectrum and 56 steep-spectrum), upper limits were set for an additional subset of 27 sources (12 flat-spectrum and 15 steep-spectrum), and 53 sources were rejected (probably extended objects). The final flat-spectrum sample is almost complete (80%), while only 49% of steep-spectrum sources have reliable detections. In the flat-spectrum sample the median fractional polarization is % and the mean %, and no sources have . The median fractional polarization for the flat-spectrum sources included in the NVSS catalogue is about a factor 2 lower than that at 18.5 GHz. However, a relevant increase in the polarization fraction at 18.5 GHz is noticed only for sources with at 1.4 GHz.(ii) Reference [49] presented a 3.5 mm polarimetric survey of ERS using the IRAM 30 m Telescope. Their sample consists of 145 flat-spectrum sources with and flux density Jy at 86 GHz. Linear polarization is detected for 76% of the sample (110 objects). They found that BL Lacs () are more strongly polarized than quasars (). This result seems to be in contradiction with the idea that quasars should be more polarized at high frequencies than BL Lacs because in the latter sources the synchrotron self-absorbed spectrum is maintained up to higher frequencies. A possible explanation provided by the authors comes from the recent evidence that the view angle of jets in quasars is smaller than that in BL Lacs [53, 54]. So, if the magnetic field is not homogeneous along the jet, a lower fractional polarization level is expected from sources better oriented to the line of sight (i.e., quasars).Moreover, for those sources with detected polarization at 15 GHz, they found that is larger than by a median factor , and about 20% of sources have the ratio larger than 4. Reference [49] suggest that this increase may be explained by a combination of two phenomena: the 86-GHz emission comes from a region with greater degree of order in the magnetic field or/and, at 15 GHz, emission is still affected by Faraday depolarization.(iii) Reference [55] looked for polarized sources in the WMAP five-year data, using a new technique named filtered fusion. They detected polarization in 13 ERS at a confidence level ≥ and polarized flux density higher than 300 mJy.

4. Statistical Properties of the Polarized Emission in ERS

In order to provide reliable estimates of ERS contamination to CMB anisotropy polarization measurements, we need to address the following main questions about polarization properties of ERS: how the fractional polarization, , varies from cm to mm wavelengths;   if the fractional polarization depends on the flux density; how polarization properties change among the different populations of ERS. In this section we address these questions on the basis of the observational data presented and discussed in the previous Section.

4.1. Flat-Spectrum Sources

We investigate more deeply the polarization of ERS observed in the AT20G survey. We consider the almost-complete sample at and, when available, we use 5- and 8-GHz measurements to separate ERS into steep- and flat-spectrum sources. In Table 1 we report statistical properties of the fractional polarization (i.e., the mean , the median , and ) for flat-spectrum sources as a function of the flux density range and the frequency. They are computed using the Survival Analysis techniques and the Kaplan-Meyer estimator as implemented in the ASURV code [56], which takes into account upper limits on the fractional polarization for estimating the above quoted quantities. In fact, when polarization is not detected, an upper limit is provided. About 10% of sources in the AT20G sample we use here have measurements only at 20 GHz (see Table 1). Although most of them should be flat-spectrum sources, for a more consistent comparison with results at 5 and 8 GHz we prefer to exclude them from the analysis.

At 20 GHz the median and the mean fractional polarization do not present any significant variation between the subsamples defined by  Jy and  Jy. At the lower frequencies, and show a moderate decrease at fainter flux densities. In particular, two-sample tests implemented in the ASURV code (e.g., the Gehan’s generalized Wilcoxon test and the Peto & Prentice generalized Wilcoxon test) yield a probability of ~10% and ~1% that the distributions of the fractional polarization for  Jy and  Jy at 5 and 8 GHz, respectively are drawn from the same parent distribution (compared to a probability of 30% at 20 GHz). At flux densities lower than 500 mJy, the high number of upper limits (%) makes our estimates unreliable.

On the other hand, a larger fractional polarization is observed as the frequency increases, with, on average, 18% of depolarization at 5 GHz with respect to 20 GHz, in agreement with results from [47]. The ASURV two-sample tests yield only a 10% probability to have the same parent distribution for at 5 and 20 GHz. For a better comparison, we show in Figure 2 the distributions of the fractional polarization discussed above as obtained by using the Kaplan-Meyer estimator. These distributions are well fitted by a log-normal distribution with and values taken from Table 1 (see, e.g., the case at 20 GHz in Figure 2). (Figure 2 clearly shows that the log-normal distribution, when averaged in each bin, is giving predictions always compatible with the observed values well inside the 1 level, except for very few bins at high levels where the statistics is very poor. We have also verified that adopting other distribution functions, for example, a truncated gaussian, we are not able to reproduce equally well the observed distributions of .)

In Figure 3 we also study the correlation of with and : no clear correlation is found between the fractional polarization and the flux density at 20 GHz (the generalized Kendall’s tau test yields a probability of no correlation). On the other hand, as expected, we see a strong correlation between the fractional polarization at 20 and 4.8 GHz (with a of no correlation). Using the Schmitt’s method from the ASURV code, we find a linear regression . This result seems to indicate that only sources with very low fractional polarization at 5 GHz (i.e., ) have a significant increase of at 20 GHz (see also Figure 3). The large offset term in the linear regression we find, however, could be also partially due to Eddington bias in the AT20G catalogue and to the large number of sources with upper limits in fractional polarization at %.

In Figure 5 we plot the mean and median values of the fractional polarization as obtained from the surveys presented in the previous section. These data allow us to cover frequencies from 1.4 to 86 GHz. Although there is large scatter in the data, we see a general increasing trend for the median fractional polarization: is ~1.5% at 1.4 GHz, around 2–2.5% in the range 5–10 GHz, and finally 2- at . The mean fractional polarization has a more linear increase with the frequency in all the samples, with varying from 2–2.5% at 1.4 GHz to 3– at GHz.

Most of the data presented in Figure 5 are coming from samples of bright sources, typically with Jy. Samples with fainter sources are the B3-VLA (mJy) and the sample discussed in [50]. The latter one provides polarization for sources with mJy. However, since at the frequencies 20 and 43 GHz the number of polarized ERS detected in the sample is less than 50 percent, we have decided to consider only flat-spectrum sources with  mJy. Spectral indices are estimated using flux densities at 4.86 and 8.46 GHz. Table 2 reports the number of detections and the corresponding values of the mean, median, and of the fractional polarization. For the two highest frequencies of the sample we estimate the median in two ways: firstly, we use the ASURV code, by taking into account the upper limits; secondly, we assume a fractional polarization ≤1% for those sources without measured polarization (values indicated in brackets in Table 2). The large spread between the two values at 43 GHz is an indication that the sample, also including upper limits, could provide biased values of statistical properties of at these frequencies.

If we compare results from the previous two samples (i.e., B3-VLA and [50]) with surveys of bright sources, we cannot find any clear evidence of higher fractional polarization in faint sources. However, larger and deeper samples of data are clearly required to settle the question.

In Figure 5 we also include the values. This quantity is important in order to estimate the angular power spectra of the polarized signal due to undetected ERS. At GHz, that is, the most interesting frequencies for CMB data analyses, sources have . In [49] optical identifications are provided for ERS in the [50] sample. Over a total of 145 objects, 107 are identified as quasars and 26 as BL Lacs: we find , 3.6% and , 4.5% for FSRQs and BL Lacs, respectively.

4.2. Steep-Spectrum Sources

The number of steep-spectrum sources with polarization measurements becomes very small at frequencies  GHz, preventing any study of the correlation between fractional polarization and flux density (see, e.g., Figure 4). In the subsample of AT20G at there are 51 steep-spectrum sources with Jy and only 25 of them have polarization detected at 20 GHz. We find that the median fractional polarization for these sources is less than 2% at all the frequencies (see Figure 5(b)). These very small values could be biased due to the small number of detected sources and to the incompleteness of the sample at faint flux densities. Nevertheless, from Figure 4 we can see a general increase of the fractional polarization between 4.8 and 20 GHz: the ratio of at 20 and 4.8 GHz is typically close to one (<1.5) for sources with , but become higher for sources weakly polarized at 4.8 GHz. In fact, the linear regression found by the Schmitt’s method implemented in the ASURV code yields .

From Figure 5 we can observe the strong increase of the fractional polarization from 1.4 GHz to 5 GHz, where Faraday depolarization is probably very relevant. The increase becomes more moderate up to 20 GHz. The large difference in the median values among different samples is perhaps related to the small samples considered and the large fraction of sources with upper limits in polarization. The largest sample of steep-spectrum sources is provided by B3-VLA (77 sources, all with detected polarization) and the Ricci et al. sample (71 sources, 15 of them with upper limits): in these samples the median fractional polarization is ~5% between 5–20 GHz and the mean varies from 5 to 6.5%.

4.3. Number Counts in Polarization of ERS at 20 GHz

We provide a first estimate of number counts in polarization, , at 20 GHz by exploiting WMAP and AT20G polarization source catalogues. For very bright sources, we use the polarized source sample detected by [55] in the 5 yr WMAP CMB anisotropy maps (see their Table 2). We exclude from number counts Fornax A, Virgo A, and Centaurus A because they are all local objects. Then, we consider the nearly-complete subsample of AT20G at and  mJy. This sample allows us to compute down to polarized fluxes of ~10 mJy. Number counts, that are reported in Table 3 (see also Figure 6), have been corrected for the estimated incompleteness of the sample (i.e., the completeness is 0.91 for mJy and 0.79 for mJy [4]). As displayed by Figure 6, number counts in total polarization, , are almost flat at 20 GHz, at least in the flux density range in which they can be estimated.

Given that an upper limit on the fractional polarization is always provided for the ERS in this sample (see [4], for more details), except for 5 very bright objects with no information on polarization, we also give a tentative estimate (Table 3, third column) of their contribution to the number counts between and by means of The sum is done over all the sources without polarization detection and is the flux density of the th object. is the probability function for the fractional polarization: we take a log-normal function with % and % (see Table 1 and Figure 3) if and zero otherwise. For the 5 objects without any polarization information, no upper limits are considered.

As shown by Table 3 and Figure 6, the contribution of sources undetected in polarization (displayed by empty squares) is generally negligible, except for bins at very faint polarization levels (i.e., less than 20 mJy), and for bins at mJy, due to the 5 sources without polarization information and flux density between 1 and 10 Jy. Only for the faintest bin is larger than the uncertainty on the number counts. Therefore, we can conclude that number counts in polarization estimated from the AT20G sample are not affected by upper limits in the polarized flux density.

5. Cosmological Evolution Models and High-Frequency Number Counts of ERS

Early evolutionary models of radio sources [17, 6062] were able to give remarkably successful fits to the majority of data coming from surveys at GHz, and down to flux densities of a few mJy. More recently, [18, 44] exploited the wealth of newly available data on luminosity functions, multifrequency source counts and redshift distributions to provide new cosmological evolution models of radio sources at frequencies  GHz and  GHz, respectively. These two models are based on the determination of the epoch-dependent luminosity functions for different source populations, starting from the local luminosity function and by adopting different luminosity evolution laws with free parameters fitted from observational data.

The predictions of high-frequency number counts of ERS provided by the above evolution models assume a simple power-law spectrum for ERS. Each source population is characterized by an “average,” fixed, spectral index, or by two spectral indices (at most). This “classical” modeling has to be considered as a first-although successful-approximation, but it gives rise to an increasing mismatch with observed high-frequency (>30 GHz) number counts currently available.

Indeed, the ERS energy spectra can be quite different from a single power-law if analyzed in large frequency intervals. Different mechanisms can be responsible for this: (a) a spectral steepening due to the more rapid energy loss of high-energy electrons with source age, that is, “electron ageing”; (b) a transition from the optically thick to the optically thin regime at high radio frequencies; (c) at different wavelengths radio emission can be dominated by different components characterized by distinct spectral behaviors [63]. In particular, a clear steepening at mm wavelengths is theoretically expected for radio flat spectra of AGN core emission [64, 65]. This steepening has been already observed in blazars (see, e.g., [66]) and has also been statistically suggested by recent analyses of different ERS samples at GHz [19, 59, 6769].

A first attempt of taking such steepening in blazar spectra into account has been done by [19]. In this work the spectral behaviour of blazars at mm wavelengths is statistically described by considering the main physical mechanisms responsible for the emission. In agreement with classical models of the synchrotron emission in the inner jets of blazars these authors interpret the high-frequency steepening observed in blazar spectra as caused, at least partially, by the transition from the optically-thick to the optically-thin regime in the observed spectra at mm wavelengths. Based on the published models of synchrotron emission from inhomogeneous, unresolved, relativistic jets [30, 65, 70], the value of the frequency at which the spectral break occurs is estimated as a function of the relevant physical parameters of AGNs: the redshift, the Doppler factor, and the linear dimension of the region (approximated as homogeneous and spherical) that is mainly responsible for the emission at the break frequency.

Recent high frequency data [57, 58], and in particular the first 1.6 Planck survey [59], have provided new important results on number counts and related statistics of ERS, and they can be used to constrain the different possible cases/models featuring a spectral break in the emission of AGN jets (see [19]). The most successful one (indicated in [19] as C2Ex, that is, the relevant emission is emitted from a more “extended” region in the inner jet of FSRQs) assumes different distributions of the break frequency for BL Lacs and FSRQs. According to this model, most of the FSRQs should show a bend in their otherwise flat spectra between 10 and 100 GHz, whereas in BL Lac spectral breaks should be typically observed at GHz (implying that the observed synchrotron radiation comes from more compact emitting regions). This dichotomy has been indeed found in the Planck ERCSC [66, 71]: almost all radio sources show very flat spectral indices at LFI frequencies, ; at HFI frequencies BL Lacs keep flat spectra () while most of FSRQs show steeper spectra, that is, . Moreover, the same model gives a remarkably good fit to all the observed data on number counts and spectral index distributions of ERS at frequencies above 5 and up to 220 GHz.

For the above reasons, hereafter we adopt the number counts provided by the “C2Ex’’ case in [19].

5.1. Number Counts of ERS in Polarization at mm Wavelengths

To assess the contamination due to undetected ERS in CMB polarization maps it is necessary to know how many sources can be found with polarized flux density (typically, a debiased estimator is however used for the polarized flux density of point sources (as, for example, , where is the uncertainty on ; see [50] for more details)) above a given flux limit . Answering this question involves the determination of source counts in polarization. This estimate is quite difficult to perform directly, by the statistical analysis of observational data, since the polarized signal is weak and many samples are usually defined by their completeness in terms only of total intensity. However, this problem can be overcome by using the source number counts in total intensity complemented by information on the statistical properties of the fractional polarization , expressed in term of some probability function [16].

Let us discuss this point in more detail. Polarization number counts can be written as where is the total number of sources with in the sample, and are the probability functions of observing in a source of flux density , a polarized intensity and a fractional polarization , respectively. Assuming that is independent of , can be determined by The probability function can be constrained from the observed distributions of the fractional polarization. In agreement with [39] we model by a log-normal distribution: where is the median of the distribution and . These formulas are strictly valid only if . However, because of the very low fractional polarization observed in ERS, the upper limit of can be effectively assumed as infinite. In Figure 2 we compared log-normal functions with the polarization distributions observed in the AT20G survey, confirming the very good fit of the model with data.

A critical point for our estimates is the variation with frequency of the fractional polarization observed in ERS. The data on ERS polarization discussed in this paper clearly show an higher fractional polarization at 10–20 GHz than at few GHz in both flat- and, more prominently, in steep-spectrum sources. At higher frequencies, data become scarce, but there are still indication of a possible further increase of the polarization fraction from 10–20 GHz to 40–90 GHz [39, 49]. For flat-spectrum sources, this increase could be due to the combination of two different effects: the polarization degree actually increases with the frequency; and BL Lacs, which are observed to be more polarized than quasars [49], become more relevant in number at higher frequencies.

In the following predictions we take into account both effects. We consider two possible values for the median and the dispersion of the log-normal distribution in order to provide a range of estimates for number counts and power spectra that could take into account uncertainties in observational data at GHz. Table 4 reports our estimates for a more “optimistic” (lower) and a more “conservative” (upper) case. (Please note that the choice of the two adjectives, “optimistic” and “conservative”, has been done under the point of view that polarized ERS do “contaminate” CMB maps in polarization and, thus, have to be removed from them.) The frequency dependence of the fractional polarization in ERS is simply modeled by means of two different sets of median and r.m.s. values at frequencies below and above about 40 GHz. Although somewhat arbitrary, this choice is motivated by the fact that this is the highest frequency at which multifrequency polarization samples of ERS are available.(i)At GHz we choose a lower and higher median and r.m.s. values of the fractional polarization based on the results displayed in Figure 5. Moreover, we require that number counts computed using values in Table 4 are (at least partially) compatible with AT20G/WMAP counts in polarization. As shown in Figure 6, the more “optimistic” case fits quite well with observational counts, especially at mJy. On the other hand, the more “conservative” case tends to overestimate number counts at lower polarized fluxes, whereas it fits particularly well with current data at high fluxes. This is not unexpected because, in this latter case, we take as median and dispersion of the fractional polarization for blazars the values provided by [49] at 90 GHz.(ii)At GHz, we make the assumption of an increase of the median fractional polarization of about 20% with respect to the corresponding cases at  GHz. This choice is not firmly constrained by current data sets; however, we have been guided by the fractional polarization levels of ERS observed by [49] at 86 GHz which fall in the middle between our present lower and upper cases.

In Figure 6 we compare the predicted number counts, , of ERS in polarized intensity at 20 and 143 GHz. By comparing the two panels of the same Figure we can appreciate that decreases by a factor at 143 GHz at Jy, in comparison with the estimated at 20 GHz. At still fainter polarized intensities, Jy, this decrement increases due to the fact that the contribution of steep-spectrum ERS, which are on average more polarized, becomes less and less relevant at high CMB frequencies. Finally, Figure 6(b) indicates a possible higher relative contribution to total polarized counts, , coming from BL Lacs with respect to FSRQs at 143 GHz. This is a direct consequence of the present analysis, which is based on the observations discussed by [49].

The good agreement with observational data at 20 GHz, as well as the reliability of our predictions on ERS number counts in flux density, , gives us confidence in making extrapolations of integral number counts, , in total linear polarization at Planck frequencies (see Table 5). It is expected that Planck LFI will be able to detect sources in polarization down to  mJy at 30 GHz and 300–400 mJy at 44 and 70 GHz, whereas Planck HFI, thanks to a better resolution and sensitivity, should reach polarized flux limits of mJy. In both cases, however, a catalogue of only a few tens of compact ERS should be provided by Planck in polarization.

The above quoted detection limits in total polarization are estimated by the performances foreseen for the Planck LFI and HFI instruments and are based on prelaunch measurements [12, 13] of the detectors calibration and capabilities. These detection limits take also into account the future application to the Planck data—corresponding to the end of the nominal mission, that is, in January 2012—of new detection techniques, specifically designed for detecting compact polarized sources in CMB maps (see, e.g., [72] for a recent discussion on the subject). These techniques have already been applied with success to WMAP 5 yr maps [55], improving on the results published by the WMAP team on the same data set.

6. Predictions on the Contribution of ERS to the CMB E- and B-Modes

As for temperature fluctuations, the analysis of CMB polarization measurements is usually made by the estimate of angular power spectra, that is, E- and B-modes spectra [8, 9]. B-mode polarization, that arises only from tensor perturbations, is expected to be extremely weak, and even for the most optimistic cases foreseen by inflationary models the rms signal is only a fraction of K, less than 1 per cent of the level of temperature anisotropies at degree scales. On the other hand, polarization of foregrounds (and in particular of extragalactic sources) is equally shared between E- and B-modes [73]. ERS are therefore expected to dominate the sky B-mode polarization at subdegree angular scales at frequencies GHz [74].

By using the statistical characterization of the fractional polarization described in the previous Sections, we are able to estimate the polarized angular power spectra given by undetected ERS in CMB anisotropy maps. First of all, we assume that ERS follow a Poisson distribution in the sky. The contribution of clustered ERS to the angular power spectrum of CMB temperature anisotropy is in fact small and can be neglected, if ERS are not subtracted down to faint flux limits (the signal due to clustered ERS becomes more relevant only at relatively low fluxes, that is,  mJy, [75, 76]).

It is well known that an ensemble of Poisson distributed point sources gives rise to a flat power spectrum of temperature fluctuations [77]. For a sample of sources with flux density below some cut-off , the amplitude of this white noise spectrum is given by where and are, respectively, the total and the differential number of sources per steradian, and is the conversion factor from brightness to temperature, that is, K/(Jy  and GHz.

As an analogy to the expression of the angular power spectrum in total intensity or CMB temperature, it is possible to define the angular power spectrum for the Stokes parameters and . Because point sources contribute, on average, equally to , and to the E-, B-mode power spectra, we assume (we generally refer to them as polarization spectra). Following the treatment given by [16], we have that where the Stokes parameter is written in terms of and of the polarization angle in the chosen reference system, . The factor 1/2 arises because of the uniform distribution of polarization angles. It is easy to demonstrate that the cross-correlation temperature-polarization spectra are null, for example, .

In (3) and (6) we have assumed that the fractional polarization is independent of the total intensity of the source. Observations for flat-spectrum sources in the flux density range mJy seem to support this assumption (see Section 4.1). We expect that this is maintained also at fainter fluxes but only if FSRQs and BL Lacs are separated into two different populations. On the other hand, this hypothesis may not be true for steep-spectrum sources: a clear anticorrelation between and is observed from data at 1.4 GHz, whereas, at higher frequencies, the lack of large samples of steep-spectrum sources does not allow to determine it. In any case, steep-spectrum ERS are giving a negligible contribution to number counts at GHz and, thus, this lack of information does not affect our current predictions.

In Figure 7 we present the results on ERS polarization power spectra for the six Planck frequencies where ERS are relevant. The value of is taken according to Table 4 for the different radio source populations. Moreover, we consider two cutoffs in flux density: , and 0.1 Jy. The former value is close to the completeness limit obtained by the Planck ERSCS for the frequency channels GHz [59]; the latter one can be seen as a (somewhat optimistic) reference value for the Planck high frequency channels, or a reference value for future experiments.

In Figure 7 we also plot the CMB power spectrum for the E-mode and for the B-mode with a tensor-to-scalar ratio , 0.01 and 0.001. In this way we can have an indication of the level of ERS subtraction required to allow the detection of a gravitational wave induced primordial CMB B-mode signal. ERS should not be a strong constraint on detecting E-mode polarization or the B-modes with . On the other hand, a primordial CMB B-mode signal corresponding to lower values requires the subtraction of ERS down to flux detection limits of ~100 mJy, which will not be easy, or even possible, with the Planck sensitivity.

7. Conclusions

In this contribution we have reviewed recently available polarimetric surveys of extragalactic radio sources at frequencies GHz. These data point out that the typical intrinsic fractional polarization of ERS is around 2–5% of the total flux density, , of the source even at frequencies as high as 20 GHz, and that in very few objects the fractional polarization is . This may be due to the low degree of uniformity of magnetic fields in the internal part of AGN jets and in lobes. Faraday depolarization is probably the cause of the large number of sources with a very low level of polarization, that is, , at GHz frequencies, which also explains the strong increase of the fractional polarization observed at GHz in those objects. This conclusion is supported by high or extreme values of rotation measure RM  rad  observed in some blazars (e.g., [37, 39]).

Moreover, we have studied how polarization properties of ERS change from cm to mm wavelengths. For flat-spectrum sources a weak but constant increase of fractional polarization is observed, with median (mean) values varying from 1.5% (2–2.5%) at 1.4 GHz, to 2–2.5% (2.5–3%) at 5–10 GHz and 2-3% (3–3.5%) at 10–40 GHz. Indications that fractional polarization in blazars could further increase above 40 GHz come from the recent works by [39, 49]. On the other hand, a significantly higher fractional polarization is typically found in steep-spectrum sources, especially at high frequencies (median is 4-5% at 10–20 GHz and mean between 5 and 6.5%). However, because of incompleteness of the samples and of the small number of steep-spectrum sources in surveys at GHz, current observations could be biased by high-polarized objects.

In general, we do not find any dependence of the fractional polarization of ERS with the flux density at high radio frequencies. However, more conclusive evidences require larger and deeper surveys. Nevertheless, an anticorrelation between and in blazars is expected when very faint flux densities are considered (see Section 4). For the flux-density ranges covered by available surveys, flat-spectrum sources are dominated by quasars (FSRQs). These objects are typically less polarized than BL Lacs [49], which become increasingly relevant at fainter flux densities and become the dominant population at  mJy [18, 78].

We also discuss a formalism to estimate ERS number counts in polarization and to predict the contribution of unresolved ERS to angular power spectra at CMB frequencies. As a first application, we attempt to predict how many polarized ERS the Planck Satellite will be able to detect in the different channels sensible to polarization measurements: we expect that only a dozen polarized ERS could be detected by Planck LFI, and a few tens at the HFI frequencies. Although the number of Planck detected sources is low, these data will allow us to study the frequency dependence of the fractional polarization in a wide range of frequencies, from 30 to 353 GHz, thus providing original and valuable information on polarization properties of ERS in the innermost regions of AGN jets.

Finally, our results on polarization power spectra demonstrate that ERS should not be a strong contaminant to the CMB E-mode polarization when observing at frequencies GHz. Moreover, it seems unlikely that ERS will have a significant impact on our ability to detect the B-mode polarization from primordial graviational waves if . On the contrary, if the cosmological B-mode signal is fainter, some strategy will be required to subtract the confusion noise produced by radio sources with flux-density Jy. At sub-mm wavelengths, where radio sources become less and less relevant, confusion noise of point sources may be dominated by dusty galaxies (e.g., [79]). Polarization of these objects should be relatively low, but they are expected to be significantly clustered. Although the level of this effect is still quite uncertain, polarization spectra from dusty galaxies could begin to dominate the one produced by ERS already at 200–300 GHz, at least at angular scales relevant for the detection of the B-modes.

Acknowledgments

The authors thank the referee, R. B. Partridge, for his insightful comments and criticisms that helped a lot in clarifying the main assumptions of the paper and also in improving its final presentation. M. Tucci acknowledges financial support from the French “Centre National d’ Études Spatiales” (CNES). L. Toffolatti acknowledges partial financial support from the Spanish Ministry of Science and Innovation (MICINN), under project AYA2010-21766-C03-01.