The calculation of the characteristic function of the signal fluctuations due to clustered astrophysical sources is performed in this paper. For the typical case of power-law differential number counts and two-point angular correlation function, we present an extension of Zolotarev's theorem that allows us to compute the cumulants of the logarithm of the absolute value of the intensity. As a test, simulations based on recent observations of radio galaxies are then carried out, showing that these cumulants can be very useful for determining the fundamental parameters defining the number counts and the correlation. If the angular correlation scale of the observed source population is known, the method presented here is able to obtain estimators of the amplitude and slope of the power-law number counts with mean absolute errors that are one order of magnitude better than previous techniques, that did not take into account the correlation. Even if the scale of correlation is not well known, the method is able to estimate it and still performs much better than if the effect of correlations is not considered.

The statistical analysis of large sample of strong lensing events can be a powerful tool to extract astrophysical and/or cosmological valuable information. However, the number of such events is still relatively low, mostly because of the lengthily observational validation process on individual events. In this work we propose a new methodology with a statistical selection approach in order to increase by a factor of $\sim 5$ the number of such events. Although the methodology can be applied to address several selection problems, it has particular benefits in the case of the identification of strongly lensed galaxies: objectivity, minimal initial constrains in the main parameter space, preservation of the statistical properties. The proposed methodology is based on the Bhattacharyya distance as a measure of the similarity between probability distributions of properties of two different cross-matched galaxies. The particular implementation for the aim of this work is called SHALOS and it combines the information of four different properties of the pair of galaxies: angular separation, luminosity percentile, redshift and optical/sub-mm flux density ratio. The SHALOS method provided a ranked list of strongly lensed galaxies. The number of candidates for the final associated probability, $P_{tot}>0.7$, is 447 with an estimated mean amplification factor of 3.12 for an halo with a typical cluster mass. Additional statistical properties of the SHALOS candidates, as the correlation function or the source number counts, are in agreement with previous results indicating the statistical lensing nature of the selected sample.

We present a final description of the data-processing pipeline for the Planck, Low Frequency Instrument (LFI), implemented for the 2018 data release. Several improvements have been made with respect to the previous release, especially in the calibration process and in the correction of instrumental features such as the effects of nonlinearity in the response of the analogue-to-digital converters. We provide a brief pedagogical introduction to the complete pipeline, as well as a detailed description of the important changes implemented. Self-consistency of the pipeline is demonstrated using dedicated simulations and null tests. We present the final version of the LFI full sky maps at 30, 44, and 70 GHz, both in temperature and polarization, together with a refined estimate of the Solar dipole and a final assessment of the main LFI instrumental parameters.

This paper presents the Planck Multi-frequency Catalogue of Non-thermal (i.e. synchrotron-dominated) Sources (PCNT) observed between 30 and 857 GHz by the ESA Planck mission. This catalogue was constructed by selecting objects detected in the full mission all-sky temperature maps at 30 and 143 GHz, with a signal-to-noise ratio (S/N)>3 in at least one of the two channels after filtering with a particular Mexican hat wavelet. As a result, 29400 source candidates were selected. Then, a multi-frequency analysis was performed using the Matrix Filters methodology at the position of these objects, and flux densities and errors were calculated for all of them in the nine Planck channels. The present catalogue is the first unbiased, full-sky catalogue of synchrotron-dominated sources published at millimetre and submillimetre wavelengths and constitutes a powerful database for statistical studies of non-thermal extragalactic sources, whose emission is dominated by the central active galactic nucleus. Together with the full multi-frequency catalogue, we also define the Bright Planck Multi-frequency Catalogue of Non-thermal Sources PCNTb, where only those objects with a S/N>4 at both 30 and 143 GHz were selected. In this catalogue 1146 compact sources are detected outside the adopted Planck GAL070 mask; thus, these sources constitute a highly reliable sample of extragalactic radio sources. We also flag the high-significance subsample PCNThs, a subset of 151 sources that are detected with S/N>4 in all nine Planck channels, 75 of which are found outside the Planck mask adopted here. The remaining 76 sources inside the Galactic mask are very likely Galactic objects.

We estimate the average fractional polarisation at 143, 217 and 353 GHz of a sample of 4697 extragalactic dusty sources by applying stacking technique. The sample is selected from the second version of the Planck Catalogue of Compact Sources at 857 GHz, avoiding the region inside the Planck Galactic mask (fsky ~ 60 per cent). We recover values for the mean fractional polarisation at 217 and 353 GHz of (3.10 \pm 0.75) per cent and (3.65 \pm 0.66) per cent, respectively, whereas at 143 GHz we give a tentative value of (3.52 \pm 2.48) per cent. We discuss the possible origin of the measured polarisation, comparing our new estimates with those previously obtained from a sample of radio sources. We test different distribution functions and we conclude that the fractional polarisation of dusty sources is well described by a log-normal distribution, as determined in the radio band studies. For this distribution we estimate {\mu}_{217GHz} = 0.3 \pm 0.5 (that would correspond to a median fractional polarisation of {\Pi}_{med} = (1.3 \pm 0.7) per cent) and {\mu}_{353GHz} = 0.7 \pm 0.4 ({\Pi}_{med} = (2.0 \pm 0.8) per cent), {\sigma}_{217GHz} = 1.3 \pm 0.2 and {\sigma}_{353GHz} = 1.1 \pm 0.2. With these values we estimate the source number counts in polarisation and the contribution given by these sources to the CMB B-mode angular power spectrum at 217, 353, 600 and 800 GHz. We conclude that extragalactic dusty sources might be an important contaminant for the primordial B-mode at frequencies > 217 GHz.

In this work we apply the stacking technique to estimate the average fractional polarisation from 30 to 353 GHz of a primary sample of 1560 compact sources - essentially all radio sources - detected in the 30 GHz Planck all-sky map and listed in the second version of the Planck Catalogue of Compact Sources (PCCS2). We divide our primary sample in two subsamples according to whether the sources lay (679 sources) or not (881 sources) inside the sky region defined by the Planck Galactic mask (fsky ~ 60 per cent) and the area around the Magellanic Clouds. We find that the average fractional polarisation of compact sources is approximately constant (with frequency) in both samples (with a weighted mean over all the channels of 3.08 per cent outside and 3.54 per cent inside the Planck mask). In the sky region outside the adopted mask, we also estimate the {\mu} and {\sigma} parameters for the log-normal distribution of the fractional polarisation, finding a weighted mean value over all the Planck frequency range of 1.0 for {\sigma} and 0.7 for {\mu} (that would imply a weighted mean value for the median fractional polarisation of 1.9 per cent).

The Second Planck Catalogue of Compact Sources is a catalogue of sources detected in single-frequency maps from the full duration of the Planck mission and supersedes previous versions of the Planck compact source catalogues. It consists of compact sources, both Galactic and extragalactic, detected over the entire sky. Compact sources detected in the lower frequency channels are assigned to the PCCS2, while at higher frequencies they are assigned to one of two sub-catalogues, the PCCS2 or PCCS2E, depending on their location on the sky. The first of these catalogues covers most of the sky and allows the user to produce subsamples at higher reliabilities than the target 80% integral reliability of the catalogue. The PCCS2E contains sources detected in sky regions where the diffuse emission makes it difficult to quantify the reliability of the detections. Both the PCCS2 and PCCS2E include polarization measurements, in the form of polarized flux densities, or upper limits, and orientation angles for all seven polarization-sensitive Planck channels. The improved data-processing of the full-mission maps and their reduced noise levels allow us to increase the number of objects in the catalogue, improving its completeness for the target 80 % reliability as compared with the previous versions, the PCCS and ERCSC catalogues.

The Planck Catalogue of Compact Sources (PCCS) is the catalogue of sources detected in the first 15 months of Planck operations, the "nominal" mission. It consists of nine single-frequency catalogues of compact sources, both Galactic and extragalactic, detected over the entire sky. The PCCS covers the frequency range 30--857\,GHz with higher sensitivity (it is 90% complete at 180 mJy in the best channel) and better angular resolution (from ~33' to ~5') than previous all-sky surveys in this frequency band. By construction its reliability is >80% and more than 65% of the sources have been detected at least in two contiguous Planck channels. In this paper we present the construction and validation of the PCCS, its contents and its statistical characterization.

(abridged for arXiv) We make use of the Planck all-sky survey to derive number counts and spectral indices of extragalactic sources -- infrared and radio sources -- from the Planck Early Catalogue (ERCSC) at 100 to 857GHz. Our sample contains, after the 80% completeness cut, between 122 and 452 and sources, with flux densities above 0.3 and 1.9Jy at 100 and 857GHz, over about 31 to 40% of the sky. Using Planck HFI, all the sources have been classified as either dust-dominated or synchrotron-dominated on the basis of their spectral energy distributions (SED). Our sample is thus complete, flux-limited and color-selected to differentiate between the two populations. We find an approximately equal number of synchrotron and dusty sources between 217 and 353GHz; at 353GHz or higher (or 217GHz and lower) frequencies, the number is dominated by dusty (synchrotron) sources, as expected. For most of the sources, the spectral indices are also derived. We provide for the first time counts of bright sources from 353 to 857GHz and the contributions from dusty and synchrotron sources at all HFI frequencies in the key spectral range where these spectra are crossing. The observed counts are in the Euclidean regime. The number counts are compared to previously published data (earlier Planck, Herschel, BLAST, SCUBA, LABOCA, SPT, and ACT) and models taking into account both radio or infrared galaxies. We derive the multi-frequency Euclidean level and compare it to WMAP, Spitzer and IRAS results. The submillimetre number counts are not well reproduced by current evolution models of dusty galaxies, whereas the millimetre part appears reasonably well fitted by the most recent model for synchrotron-dominated sources. Finally we provide estimates of the local luminosity density of dusty galaxies, providing the first such measurements at 545 and 857GHz.

In this paper we describe the state-of-the art status of multi-frequency detection techniques for compact sources in microwave astronomy. From the simplest cases where the spectral behaviour is well-known (i.e. thermal SZ clusters) to the more complex cases where there is little a priori information (i.e. polarized radio sources) we will review the main advances and the most recent results in the detection problem.

The data reported in Planck's Early Release Compact Source Catalogue (ERCSC) are exploited to measure the number counts (dN/dS) of extragalactic radio sources at 30, 44, 70, 100, 143 and 217 GHz. Due to the full-sky nature of the catalogue, this measurement extends to the rarest and brightest sources in the sky. At lower frequencies (30, 44, and 70 GHz) our counts are in very good agreement with estimates based on WMAP data, being somewhat deeper at 30 and 70 GHz, and somewhat shallower at 44 GHz. Planck's source counts at 143 and 217 GHz join smoothly with the fainter ones provided by the SPT and ACT surveys over small fractions of the sky. An analysis of source spectra, exploiting Planck's uniquely broad spectral coverage, finds clear evidence of a steepening of the mean spectral index above about 70 GHz. This implies that, at these frequencies, the contamination of the CMB power spectrum by radio sources below the detection limit is significantly lower than previously estimated.

We present results on spectral index distributions, number counts, redshift distribution and other general statistical properties of extragalactic point sources in the NEWPS5 sample L\'opez-Caniego et al. (2007). The flux calibrations at all the WMAP channels have been reassessed both by comparison with ground based observations and through estimates of the effective beam areas. The two methods yield consistent statistical correction factors. A search of the NED has yielded optical identifications for 89% of sources in the complete sub-sample of 252 sources with S/N>5 and S>1.1 Jy at 23 GHz; 5 sources turned out to be Galactic and were removed. The NED also yielded redshifts for 92% of the extragalactic sources at |b|>10deg. Their distribution was compared with model predictions; the agreement is generally good but a possible discrepancy is noted. Using the 5 GHz fluxes from the GB6 or PMN surveys, we find that 76% of the 191 extragalactic sources with S_23GHz>1.3,Jy can be classified as flat-spectrum sources between 5 and 23 GHz. A spectral steepening is observed at higher frequencies: only 59% of our sources are still flat-spectrum sources between 23 and 61 GHz and the average spectral indexes steepen from <alpha_5^23>= 0.01\pm 0.03 to <alpha_41^61>= 0.37\pm 0.03. We think, however, that the difference may be due to a selection effect. The source number counts have a close to Euclidean slope and are in good agreement with the predictions of the cosmological evolution model by De Zotti et al. (2005). The observed spectral index distributions were exploited to get model-independent extrapolations of counts to higher frequencies. The risks of such operations are discussed and reasons of discrepancies with other recent estimates are clarified.

The possible contributions of the various classes of extragalactic sources (including, in addition to the canonical radio sources, GHz Peaked Spectrum sources, advection-dominated sources, starburst galaxies, high-redshift proto-spheroidal galaxies) to the arcminute scale fluctuations measured by the CBI, BIMA, and ACBAR experiments are discussed. At 30 GHz, fluctuations due to radio sources undetected by ancillary low-frequency surveys may be higher than estimated by the CBI and BIMA groups. High-redshift dusty galaxies, whose fluctuations may be strongly enhanced by the effect of clustering, could contribute to the BIMA excess signal, and dominate at 150 GHz (the ACBAR frequency). Moreover, in the present data situation, the dust emission of these high-redshift sources set an unavoidable limit to the detection of primordial CMB anisotropies at high multipoles, even at frequencies as low as $\simeq 30$ GHz. It is concluded that the possibility that the excess power at high multipoles is dominated by unsubtracted extragalactic sources cannot be ruled out. On the other hand, there is room for a contribution from the Sunyaev-Zeldovich effect within clusters of galaxies, with a density fluctuation amplitude parameter $\sigma_8$ consistent with the values preferred by current data.

We present predictions on the angular power spectrum of CMB fluctuations due to extragalactic point sources (EPS) by using a method for simulating realistic 2D distributions of clustered EPS. Both radio and far--IR selected source populations are taken into account. To analyze different clustering scenarios, we exploit angular power spectra of EPS, $P(k)$, estimated either by data coming from currently available surveys or by means of theoretical predictions. By adopting the source number counts predicted by the Toffolatti et al. (1998) evolution model -- capable of accounting well for the available data at radio cm wavelengths -- we are able to reproduce current data on the two--point angular correlation functions, $w(\theta)$, of radio sources. We can confirm that the detection of primordial CMB anisotropies is not hampered by undetected clustered sources at frequencies $\leq 150-200$ GHz. On the other hand, our current findings show that at higher frequencies the clustering signal could severely reduce the detectability of intrinsic CMB anisotropies, thus confirming previous theoretical predictions. We also show that unsubtracted EPS can account for the excess signal at high multipoles detected by recent CMB anisotropy experiments. Moreover, the additional power due to the clustering of sources gives rise to a small but not negligible contribution to the same excess signal. As a final result, we also present an example of a currently feasible {\it realistic map} of EPS at 70 GHz, by taking into account data on bright detected sources as well as the previously quoted model for number counts.

All the analyses of Cosmic Microwave Background (CMB) temperature maps up--to--date show that CMB anisotropies follow a Gaussian distribution. On the other hand, astrophysical foregrounds which hamper the detection of the CMB angular power spectrum, are not Gaussian distributed on the sky. Therefore, they should give a sizeable contribution to the CMB bispectrum. In fact, the first year data of the Wilkinson Microwave Anisotropy Probe (WMAP) mission have allowed the {\it first} detection of the extragalactic source contribution to the CMB bispectrum at 41 GHz and, at the same time, much tighter limits than before to non--Gaussian primordial fluctuations. In view of the above and for achieving higher precision in current and future CMB measurements of non--Gaussianity, in this paper we discuss a comprehensive assessment of the bispectrum due to either uncorrelated and clustered extragalactic point sources in the whole frequency interval around the CMB intensity peak. Our calculations, based on current cosmological evolution models for sources, show that the reduced angular bispectrum due to point sources, $b_{ps}$, should be detectable in all WMAP and Planck frequency channels. We also find agreement with the results on $b_{ps}$ at 41 GHz coming from the analysis of the first year WMAP data. Moreover, by comparing $b_{ps}$ with the primordial reduced CMB bispectrum, we find that only the peak value of the primordial bispectrum (which appears at $l\simeq 200$) results greater than $b_{ps}$ in a frequency window around the intensity peak of the CMB. The amplitude of this window basically depends on the capability of the source detection algorithms (i.e., on the achievable flux detection limit, $S_{lim}$, for sources).

Although the Planck Surveyor mission is optimized to map the cosmic microwave background anisotropies, it will also provide extremely valuable information on astrophysical phenomena. We review our present understanding of Galactic and extragalactic foregrounds relevant to the mission and discuss on one side, Planck's impact on the study of their properties and, on the other side, to what extent foreground contamination may affect Planck's ability to accurately determine cosmological parameters. Planck's multifrequency surveys will be unique in their coverage of large areas of the sky (actually, of the full sky); this will extend by two or more orders of magnitude the flux density interval over which mm/sub-mm counts of extragalactic sources can be determined by instruments already available (like SCUBA) or planned for the next decade (like the LSA-MMA or the space mission FIRST), which go much deeper but over very limited areas. Planck will thus provide essential complementary information on the epoch-dependent luminosity functions. Bright radio sources will be studied over a poorly explored frequency range where spectral signatures, essential to understand the physical processes that are going on, show up. The Sunyaev-Zeldovich effect, with its extremely rich information content, will be observed in the direction of a large number of rich clusters of Galaxies. Thanks again to its all sky coverage, Planck will provide unique information on the structure and on the emission properties of the interstellar medium in the Galaxy. At the same time, the foregrounds are unlikely to substantially limit Planck's ability to measure the cosmological signals. Even measurements of polarization of the primordial Cosmic Microwave background fluctuations appear to be feasible.

We simulate Planck observations by adopting a detailed model of the microwave sky including monopole, dipole, anisotropies of the cosmic microwave background (CMB) and galactic and extragalactic foregrounds. We estimate the impact of main beam optical aberrations on CMB anisotropy measurements in presence of extragalactic source fluctuations and we discuss the main implications for the Planck telescope design. By analysing the dipole pattern, we quantify the Planck performance in the determination of CMB spectral distortion parameters in presence of foreground contaminations.

Confusion noise due to extragalactic sources is a fundamental astrophysical limitation for experiments aimed at accurately determining the power spectrum of the Cosmic Microwave Background (CMB) down to arcmin angular scales and with a sensitivity $\Delta T/T \simeq 10^{-6}$. At frequencies $\lsim 200-300$ GHz, the most relevant extragalactic foreground hampering the detection of intrinsic CMB anisotropies is constituted by radio loud Active Galactic Nuclei (AGN), including ``flat--spectrum'' radiogalaxies, quasars, BL-LACs and blazars. We review our present understanding of astrophysical properties, spectra, and number counts of the above classes of sources. We also study the angular power spectrum of fluctuations due both to Poisson distributed and clustered radio sources and give preliminary predictions on the power spectrum of their polarized components. Furthermore, we discuss the capabilities of future space missions (NASA's MAP, Bennett et al. 1995; ESA's Planck Surveyor, Bersanelli et al. 1996) in studying bright radio sources over an almost unexplored frequency interval where spectral signatures, essential for the understanding of the physical processes, show up.