Virial shocks around galaxy clusters are expected to show a cutoff in the thermal Sunyaev-Zel'dovich (SZ) signal, coincident with a leptonic ring. However, until now, leptonic virial signals were reported only in Coma and in stacked Fermi-LAT clusters, and an SZ virial shock signal was reported only in A2319. We point out that a few clusters --- presently Coma, A2319, and A2142 --- already show a sharp drop in Planck SZ pressure near the virial radius, coincident with a weak LAT $\gamma$-ray excess. These signatures are shown to trace the virial shocks of the clusters, at joint medium to high confidence levels. The electron acceleration rates inferred from $\gamma$-rays are consistent with previous measurements. The combined signal allows a separate measurement of the $\sim0.5\%$ acceleration efficiency and of the accretion rate. Lower limits of order a few are imposed on the shock Mach numbers.

We present the first combination of thermal Sunyaev-Zel'dovich (tSZ) map with a multi-frequency quality assessment of the sky pixels based on Artificial Neural Networks (ANN) aiming at detecting tSZ sources from sub-millimeter observations of the sky by Planck. We construct an adapted full-sky ANN assessment on the fullsky and we present the construction of the resulting filtered and cleaned tSZ map, MILCANN. We show that this combination allows to significantly reduce the noise fluctuations and foreground residuals compared to standard tSZ maps. From the MILCANN map, we constructed the HAD tSZ source catalog that consists of 3969 sources with a purity of 90\%. Finally, We compare this catalog with ancillary catalogs and show that the galaxy-cluster candidates in the HAD catalog are essentially low-mass (down to $M_{500} = 10^{14}$ M$_\odot$) high-redshift (up to $z \leq 1$) galaxy cluster candidates.

We present a synthetic galaxy lightcone specially designed for narrow-band optical photometric surveys. To reduce time-discreteness effects, unlike previous works, we directly include the lightcone construction in the \texttt{L-Galaxies} semi-analytic model applied to the subhalo merger trees of the {\tt Millennium} simulation. Additionally, we add a model for the nebular emission in star-forming regions, which is crucial for correctly predicting the narrow/medium-band photometry of galaxies. Explicitly, we consider, individually for each galaxy, the contribution of 9 different lines: $\rm Ly{\alpha}$ (1216\AA), \Hb (4861\AA), \Ha (6563\AA), {\oii} (3727\AA, 3729\AA), {\oiii} (4959\AA, 5007\AA), $\rm [\ion{Ne}{III}]$ (3870\AA), {\oi} (6300\AA), $\rm [\ion{N}{II}]$ (6548\AA, 6583\AA), and $\rm [\ion{S}{II}]$ (6717\AA, 6731\AA). We validate our lightcone by comparing galaxy number counts, angular clustering, and \Ha, \Hb, {\oii} and {\oiiiFd} luminosity functions to a compilation of observations. As an application of our mock lightcones, we generate catalogues tailored for J-PLUS, a large optical galaxy survey featuring 5 broad and 7 medium band filters. We study the ability of the survey to correctly identify, with a simple \textit{three filter method}, a population of emission-line galaxies at various redshifts. We show that the $4000\AA$ break in the spectral energy distribution of galaxies can be misidentified as line emission. However, all significant excess (larger than 0.4 magnitudes) can be correctly and unambiguously attributed to emission line galaxies. Our catalogues are publicly released to facilitate their use in interpreting narrow-band surveys and for quantifying the impact of line emission in broad band photometry.

The accurate determination of the galaxy cluster mass-observable relations is one of the major challenge of modern astrophysics and cosmology. We present a new statistical methodology to constrain the evolution of the mass-observable relations. Instead of measuring individual mass of galaxy clusters, we only consider large scale homogeneity of the Universe. In this case, we expect the present galaxy cluster mass function to be the same everywhere in the Universe. Using relative abundance matching, we contraint the relation between the richness, $\lambda(z)$, and the expected present mass, $M(t_0)$, of galaxy clusters. We apply this approach to the redMaPPer galaxy cluster catalogue in 10 redshift bins from $z=0.1$ to $0.6$. We found that the $\lambda(z)$-$M(t_0)$ relation is not evolving from $z=0.1$ to $0.4$, whereas it starts to significantly evolve at higher redshift. This results implies that the redMaPPer richness appears to be a better proxy for the expected present-day galaxy cluster mass than for the mass at the observational redshift. Assuming cosmology and galaxy cluster mass accretion history, it is possible to convert $M(t_0)$ to the mass at the galaxy cluster redshift $M(t_z)$. We found a significant evolution of the $\lambda(z)$-$M(t_z)$ over all the covered redshift range. Consequently, we provide a new redshift-dependent richness-mass relation for the redMaPPer galaxy cluster catalogue. This results demonstrates the efficiency of this new methodology to probe the evolution of scaling relations compared to individual galaxy cluster mass estimation.

Galaxy cluster number count has been proven to be a powerful cosmological probe. However, cosmological constraints established with galaxy cluster number count are highly dependent on the calibration of the mass-observable relations. Thanks to its nearly mass independence the specific mass accretion rate of galaxy clusters is nearly insensitive to the calibration of mass-observable relations. The study of galaxy cluster number count evolution allows to probe the galaxy cluster mass accretion history in the context of an homogenous Universe. In this paper, we use relative abundance matching to infer the galaxy cluster mass accretion rate (MAR) for $z \in [0.0,0.6[$. Then, we use the MAR to set cosmological constraints. We found that this cosmological probe is sensitive to $\sigma_8 \Omega_{\rm m}^{-0.3} H_0^{-0.2}$ whereas the galaxy cluster count is sensitive to $\sigma_8 \Omega_{\rm m}^{0.3}$. We used the second $Planck$ catalog of Sunyaev-Zel'dovich sources and we derive $\sigma_8 \Omega_{\rm m}^{-0.3} H_0^{-0.2} = 0.75 \pm 0.06$. This results is consistent with cosmological constraints derived from galaxy clusters number counts, angular power spectrum, and cosmic microwave background analyses. Therefore, the MAR is a key cosmological probe that can break the $\sigma_8$-$\Omega_{\rm m}$ degeneracy and that is not sensitive to the calibration of the mass-observable relations and does not requires a parametric form for the galaxy cluster mass-function.

Rotational transition lines of CO are one of the major tracers used to study star forming regions and Galactic structures. A large number of observations of CO rotational lines are covering the galactic plane, recently the $Planck$ collaboration release the first full sky coverage of J=1-0, J=2-1, and J=3-2 CO lines at a resolution of 10', 5', and 5' FWHM. However, the measured signal in $Planck$ detectors is integrated over large bandpass $(\Delta \nu)/\nu \simeq 0.2$. Consequently, the derived CO products are composite maps including several rotational lines. In the 100 GHz $Planck$ channel, the two main lines are the J=1-0 transitions of $^{12}$CO and $^{13}$CO at 115 and 110 GHz respectively. In the present paper, we present and applied a method to construct separate CO integrated intensity maps for the two isotopes. The measurement of the $^{13}$CO rotational transitions provide an unprecedented all-sky view of the Galaxy for this isotope.

We present the optimal reconstruction of the thermal Sunyaev-Zel'dovich (tSZ) effect signal based on the combination of a heterogeneous dataset consisting of ACT and Planck data, with different numbers of channels, angular resolutions and noise levels. We combine both datasets using two different approaches, a Matched Multi-Filter (MMF) technique and an optimised Internal Linear Combination (ILC). We show that when applying the MMF to the combination of ACT and Planck data, the size-flux degeneracy is reduced and the signal-to-noise of clusters detected with their SZ signal improves by up to a factor of three. In the case of the optimised ILC method, we show that the tSZ map is reconstructed with a resolution of $\sim 1.5$ arcmin. This is more than a factor two improvement compared with the Planck resolution, and with a very good control of noise, i.e. limited only by the intrinsic noise of the individual experiments. The combination of ACT and Planck data offers a unique opportunity to improve on the study of the pressure profiles and to study substructure in clusters through their tSZ.

Magnetized plasmas within halos of galaxies leave their footprint on the polarized anisotropies of the cosmic microwave background. The two dominant effects for astrophysical halos are Faraday rotation generating rotation of the plane of linear polarization, and Faraday conversion inducing a leakage from linear polarization to circular polarization. We revisit these sources of secondary anisotropies by computing the angular power spectra of the Faraday rotation angle and of the Faraday conversion rate by the large scale structures. To this end, we use the halo model and we pay special attention to the impact of magnetic field projections. Assuming magnetic fields of halos to be uncorrelated, we found a vanishing 2-halo term, and angular power spectra peaking at multipoles $\ell\sim10^4$. The Faraday rotation angle is dominated by the contribution of thermal electrons. For the Faraday conversion rate, we found that both thermal electrons and relativistic, non-thermal electrons contribute equally in the most optimistic case for the density and Lorentz factor of relativistic electrons, while in more pessimistic cases the thermal electrons give the dominant contribution. Assuming the magnetic field to be independent of the halo mass, the angular power spectra for both effects roughly scale with the amplitude of matter perturbations as $\sim\sigma_8^3$, and with a very mild dependence with the density of cold dark matter. Introducing a dependence of the magnetic field strength with the halo mass leads to an increase of the scaling with the amplitude of matter fluctuations, up to $\sim\sigma_8^{9.5}$ for Faraday rotation and $\sim\sigma_8^{15}$ for Faraday conversion for a magnetic field strength scaling linearly with the halo mass.

There is a lack of large samples of spectroscopically confirmed clusters and protoclusters at high redshifts, $z>$1.5. Discovering and characterizing distant (proto-)clusters is important for yielding insights into the formation of large-scale structure and on the physical processes responsible for regulating star-formation in galaxies in dense environments. The Spitzer Planck Herschel Infrared Cluster (SPHerIC) survey was initiated to identify these characteristically faint and dust-reddened sources during the epoch of their early assembly. We present Spitzer IRAC observations of 82 galaxy (proto-)cluster candidates at 1.3<$z_p$<3.0 that were vetted in a two step process: (1) using Planck to select by color those sources with the highest star-formation rates, and (2) using Herschel at higher resolution to separate out the individual red sources. The addition of the Spitzer data enables efficient detection of the central and massive brightest red cluster galaxies (BRCGs). We find that BRCGs are associated with highly significant, extended and crowded regions of IRAC sources which are more overdense than the field. This result corroborates our hypothesis that BRCGs within the Planck - Herschel sources trace some of the densest and actively star-forming proto-clusters in the early Universe. On the basis of a richness-mass proxy relation, we obtain an estimate of their mean masses which suggests our sample consists of some of the most massive clusters at z$\approx$2 and are the likely progenitors of the most massive clusters observed today.

J-PLUS is an ongoing 12-band photometric optical survey, observing thousands of square degrees of the Northern hemisphere from the dedicated JAST/T80 telescope at the Observatorio Astrof\'isico de Javalambre. T80Cam is a 2 sq.deg field-of-view camera mounted on this 83cm-diameter telescope, and is equipped with a unique system of filters spanning the entire optical range. This filter system is a combination of broad, medium and narrow-band filters, optimally designed to extract the rest-frame spectral features (the 3700-4000\AA\ Balmer break region, H$\delta$, Ca H+K, the G-band, the Mgb and Ca triplets) that are key to both characterize stellar types and to deliver a low-resolution photo-spectrum for each pixel of the sky observed. With a typical depth of AB $\sim 21.25$ mag per band, this filter set thus allows for an indiscriminate and accurate characterization of the stellar population in our Galaxy, it provides an unprecedented 2D photo-spectral information for all resolved galaxies in the local universe, as well as accurate photo-z estimates ($\Delta\,z\sim 0.01-0.03$) for moderately bright (up to $r\sim 20$ mag) extragalactic sources. While some narrow band filters are designed for the study of particular emission features ([OII]/$\lambda$3727, H$\alpha$/$\lambda$6563) up to $z < 0.015$, they also provide well-defined windows for the analysis of other emission lines at higher redshifts. As a result, J-PLUS has the potential to contribute to a wide range of fields in Astrophysics, both in the nearby universe (Milky Way, 2D IFU-like studies, stellar populations of nearby and moderate redshift galaxies, clusters of galaxies) and at high redshifts (ELGs at $z\approx 0.77, 2.2$ and $4.4$, QSOs, etc). With this paper, we release $\sim 36$ sq.deg of J-PLUS data, containing about $1.5\times 10^5$ stars and $10^5$ galaxies at $r<21$ mag.

The mass of galaxy clusters is not a direct observable, nonetheless it is commonly used to probe cosmological models. Based on the combination of all main cluster observables, that is, the X-ray emission, the thermal Sunyaev-Zel'dovich (SZ) signal, the velocity dispersion of the cluster galaxies, and gravitational lensing, the gravitational potential of galaxy clusters can be jointly reconstructed. We derive the two main ingredients required for this joint reconstruction: the potentials individually reconstructed from the observables and their covariance matrices, which act as a weight in the joint reconstruction. We show here the method to derive these quantities. The result of the joint reconstruction applied to a real cluster will be discussed in a forthcoming paper. We apply the Richardson-Lucy deprojection algorithm to data on a two-dimensional (2D) grid. We first test the 2D deprojection algorithm on a $\beta$-profile. Assuming hydrostatic equilibrium, we further reconstruct the gravitational potential of a simulated galaxy cluster based on synthetic SZ and X-ray data. We then reconstruct the projected gravitational potential of the massive and dynamically active cluster Abell 2142, based on the X-ray observations collected with XMM-Newton and the SZ observations from the Planck satellite. Finally, we compute the covariance matrix of the projected reconstructed potential of the cluster Abell 2142 based on the X-ray measurements collected with XMM-Newton. The gravitational potentials of the simulated cluster recovered from synthetic X-ray and SZ data are consistent, even though the potential reconstructed from X-rays shows larger deviations from the true potential. Regarding Abell 2142, the projected gravitational cluster potentials recovered from SZ and X-ray data reproduce well the projected potential inferred from gravitational-lensing observations. (abridged)

We present cluster counts and corresponding cosmological constraints from the Planck full mission data set. Our catalogue consists of 439 clusters detected via their Sunyaev-Zeldovich (SZ) signal down to a signal-to-noise ratio of 6, and is more than a factor of 2 larger than the 2013 Planck cluster cosmology sample. The counts are consistent with those from 2013 and yield compatible constraints under the same modelling assumptions. Taking advantage of the larger catalogue, we extend our analysis to the two-dimensional distribution in redshift and signal-to-noise. We use mass estimates from two recent studies of gravitational lensing of background galaxies by Planck clusters to provide priors on the hydrostatic bias parameter, $(1-b)$. In addition, we use lensing of cosmic microwave background (CMB) temperature fluctuations by Planck clusters as an independent constraint on this parameter. These various calibrations imply constraints on the present-day amplitude of matter fluctuations in varying degrees of tension with those from the Planck analysis of primary fluctuations in the CMB; for the lowest estimated values of $(1-b)$ the tension is mild, only a little over one standard deviation, while it remains substantial ($3.7\,\sigma$) for the largest estimated value. We also examine constraints on extensions to the base flat $\Lambda$CDM model by combining the cluster and CMB constraints. The combination appears to favour non-minimal neutrino masses, but this possibility does little to relieve the overall tension because it simultaneously lowers the implied value of the Hubble parameter, thereby exacerbating the discrepancy with most current astrophysical estimates. Improving the precision of cluster mass calibrations from the current 10%-level to 1% would significantly strengthen these combined analyses and provide a stringent test of the base $\Lambda$CDM model.

We present the joint analysis of the X-ray and SZ signals in A2319, the galaxy cluster with the highest signal-to-noise ratio in Planck maps and that has been surveyed within our XMM Cluster Outskirts Project (X-COP). We recover the thermodynamical profiles by the geometrical deprojection of the X-ray surface brightness, of the SZ comptonization parameter, and an accurate and robust spectroscopic measurements of the temperature. We resolve the clumpiness of the density to be below 20 per cent demonstrating that most of this clumpiness originates from the ongoing merger and can be associated to large-scale inhomogeneities. This analysis is done in azimuthally averaged radial bins and in eight independent angular sectors, enabling us to study in details the azimuthal variance of the recovered properties. Given the exquisite quality of the X-ray and SZ datasets, we constrain at $R_{200}$ the total hydrostatic mass, modelled with a NFW profile, with very high precision ($M_{200} = 9.76 \pm 0.16^{stat.} \pm 0.31^{syst.} \times 10^{14} M_\odot$). We identify the ongoing merger and how it is affecting differently the gas properties in the resolved azimuthal sectors. We have several indications that the merger has injected a high level of non-thermal pressure in this system: the clumping free density profile is above the average profile obtained by stacking Rosat observations; the gas mass fraction exceeds the expected cosmic gas fraction beyond $R_{500}$; the pressure profile is flatter than the fit obtained by the Planck collaboration; the entropy profile is flatter than the mean one predicted from non-radiative simulations; the analysis in azimuthal sectors has revealed that these deviations occur in a preferred region of the cluster. All these tensions are resolved by requiring a relative support of about 40 per cent from non-thermal to the total pressure at $R_{200}$.

Shocks produced by the accretion of infalling gas in the outskirt of galaxy clusters are expected in the hierarchical structure formation scenario, as found in cosmological hydrodynamical simulations. Here, we report the detection of a shock front at a large radius in the pressure profile of the galaxy cluster A2319 at a significance of $8.6\sigma$, using Planck thermal Sunyaev-Zel'dovich data. The shock is located at $(2.93 \pm 0.05) \times R_{500}$ and is not dominated by any preferential radial direction. Using a parametric model of the pressure profile, we derive a lower limit on the Mach number of the infalling gas, $\mathcal{M} > 3.25$ at 95\% confidence level. These results are consistent with expectations derived from hydrodynamical simulations. Finally, we use the shock location to constrain the accretion rate of A2319 to $\dot{M} \simeq (1.4 \pm 0.4) \times 10^{14}$ M$_\odot$ Gyr$^{-1}$, for a total mass, $M_{200} \simeq 10^{15}$ M$_\odot$.

A fraction of galaxy clusters host diffuse radio sources whose origins are investigated through multi-wavelength studies of cluster samples. We investigate the presence of diffuse radio emission in a sample of seven galaxy clusters in the largely unexplored intermediate redshift range (0.3 < z < 0.44). In search of diffuse emission, deep radio imaging of the clusters are presented from wide band (1.1-3.1 GHz), full resolution ($\sim$ 5 arcsec) observations with the Australia Telescope Compact Array (ATCA). The visibilities were also imaged at lower resolution after point source modelling and subtraction and after a taper was applied to achieve better sensitivity to low surface brightness diffuse radio emission. In case of non-detection of diffuse sources, we set upper limits for the radio power of injected diffuse radio sources in the field of our observations. Furthermore, we discuss the dynamical state of the observed clusters based on an X-ray morphological analysis with XMM-Newton. We detect a giant radio halo in PSZ2 G284.97-23.69 (z=0.39) and a possible diffuse source in the nearly relaxed cluster PSZ2 G262.73-40.92 (z=0.421). Our sample contains three highly disturbed massive clusters without clear traces of diffuse emission at the observed frequencies. We were able to inject modelled radio halos with low values of total flux density to set upper detection limits; however, with our high-frequency observations we cannot exclude the presence of RH in these systems because of the sensitivity of our observations in combination with the high z of the observed clusters.

Since the publication of the results of the Planck satellite mission in 2013, the local and early universes have been considered to be in tension in respect of the determination of amplitude of the matter density spatial fluctuations ($\sigma_8$) and the amount of matter present in the universe ($\Omega_m$). This tension can be seen as a lack of massive galaxy clusters in the local universe compared to the prediction inferred from Planck cosmic microwave background (CMB) best-fitting cosmology. In the present analysis, we perform the first detection of the cross-correlation between X-rays and CMB weak-lensing at 9.1 $\sigma$. We next combine thermal Sunyaev-Zel'dovich, X-rays, and weak-lensing angular auto and cross power spectra to determine the galaxy cluster hydrostatic mass bias. We derive $(1-b_H) = 0.70 \pm 0.05$. Considering these constraints, we observe that estimations of $\sigma_8$ in the local Universe are consistent with Planck CMB best-fitting cosmology. However, these results are in clear tension with the output of hydrodynamical simulations that favor $(1-b_H) > 0.8$.

The cosmological parameters prefered by the cosmic microwave background (CMB) primary anisotropies predict many more galaxy clusters than those that have been detected via the thermal Sunyaev-Zeldovich (tSZ) effect. This tension has attracted considerable attention since it could be evidence of physics beyond the simplest $\Lambda$CDM model. However, an accurate and robust calibration of the mass-observable relation for clusters is necessary for the comparison, which has been proven difficult to obtain so far. Here, we present new contraints on the mass-pressure relation by combining tSZ and CMB lensing measurements about optically-selected clusters. Consequently, our galaxy cluster sample is independent from the data employed to derive cosmological constrains. We estimate an average hydrostatic mass bias of $b = 0.26 \pm 0.07$, with no significant mass nor redshift evolution. This value greatly reduces the tension between the predictions of $\Lambda$CDM and the observed abundance of tSZ clusters while being in agreement with recent estimations from tSZ clustering. On the other hand, our value for $b$ is higher than the predictions from hydro-dynamical simulations. This suggests the existence of mechanisms driving large departures from hydrostatic equilibrium and that are not included in state-of-the-art simulations, and/or unaccounted systematic errors such as biases in the cluster catalogue due to the optical selection.

We present the implications for cosmic inflation of the Planck measurements of the cosmic microwave background (CMB) anisotropies in both temperature and polarization based on the full Planck survey. The Planck full mission temperature data and a first release of polarization data on large angular scales measure the spectral index of curvature perturbations to be $n_\mathrm{s} = 0.968 \pm 0.006$ and tightly constrain its scale dependence to $d n_s/d \ln k =-0.003 \pm 0.007$ when combined with the Planck lensing likelihood. When the high-$\ell$ polarization data is included, the results are consistent and uncertainties are reduced. The upper bound on the tensor-to-scalar ratio is $r_{0.002} < 0.11$ (95% CL), consistent with the B-mode polarization constraint $r< 0.12$ (95% CL) obtained from a joint BICEP2/Keck Array and Planck analysis. These results imply that $V(\phi) \propto \phi^2$ and natural inflation are now disfavoured compared to models predicting a smaller tensor-to-scalar ratio, such as $R^2$ inflation. Three independent methods reconstructing the primordial power spectrum are investigated. The Planck data are consistent with adiabatic primordial perturbations. We investigate inflationary models producing an anisotropic modulation of the primordial curvature power spectrum as well as generalized models of inflation not governed by a scalar field with a canonical kinetic term. The 2015 results are consistent with the 2013 analysis based on the nominal mission data.

The hot electrons in the intra-cluster medium produce a spectral distorsion of the cosmic microwave background (CMB) black body emission, the thermal Sunyaev-Zel'dovich effect (tSZ). This characteristic spectral distorsion is now commonly used to detect and characterize the properties of galaxy clusters. The tSZ effect spectral distorsion does not depends on the redshift, and is only slightly affected by the galaxy cluster properties via the relativistic corrections, when the electrons reach relativistic velocities. The present work proposes a linear component separation approach to extract the tSZ effect Compton parameter and relativistic corrections for next-generation CMB experiments. We demonstrated that relativistic corrections, if neglected, would induce a significant bias on galaxy cluster Compton parameter, tSZ scaling relation slope, and tSZ angular power spectrum shape measurements. We showed that tSZ relativistic corrections mapping can be achieved at high signal-to-noise ratio with a low level of contamination up to $\ell=3000$ for next generation CMB experiments. At smaller angular scales the contamination produced by infra-red emission will be a significant source of bias. Such tSZ relativistic corrections mapping enables the study of galaxy cluster temperature profile via the tSZ effect only.

The Sunyaev-Zel'dovich (SZ) effects are produced by the interaction of cosmic microwave background (CMB) photons with the ionized and diffuse gas of electrons inside galaxy clusters integrated along the line of sight. The two main effects are the thermal SZ (tSZ) produced by thermal pressure inside galaxy clusters and the kinetic SZ (kSZ) produce by peculiar motion of galaxy clusters compared to CMB rest-frame. The kSZ effect is particularly challenging to measure as it follows the same spectral behavior as the CMB, and consequently can not be separated from the CMB using spectral considerations. In this paper, we explore the feasibility of detecting the kSZ through the computation of the tSZ-CMB-CMB cross-correlation bispectrum for current and future CMB experiments. We conclude that next generation of CMB experiments will offer the possibility to detect at high S/N the tSZ-kSZ-kSZ bispectrum. This measurement will constraints the intra-cluster dynamics and the velocity field of galaxy cluster that is extremely sensitive to the growth rate of structures and thus to dark energy properties. Additionally, we also demonstrate that the tSZ-kSZ-kSZ bispectrum can be used to break the degeneracies between the mass-observable relation and the cosmological parameters to set tight constraints, up to 4\%, on the $Y-M$ relation calibration.

The Sunyaev-Zel'dovich (SZ) effect is a powerful probe of the evolution of structures in the universe, and is thus highly sensitive to cosmological parameters $\sigma_8$ and $\Omega_m$, though its power is hampered by the current uncertainties on the cluster mass calibration. In this analysis we revisit constraints on these cosmological parameters as well as the hydrostatic mass bias, by performing (i) a robust estimation of the tSZ power-spectrum, (ii) a complete modeling and analysis of the tSZ bispectrum, and (iii) a combined analysis of galaxy clusters number count, tSZ power spectrum, and tSZ bispectrum. From this analysis, we derive as final constraints $\sigma_8 = 0.79 \pm 0.02$, $\Omega_{\rm m} = 0.29 \pm 0.02$, and $(1-b) = 0.71 \pm 0.07$. These results favour a high value for the hydrostatic mass bias compared to numerical simulations and weak-lensing based estimations. They are furthermore consistent with both previous tSZ analyses, CMB derived cosmological parameters, and ancillary estimations of the hydrostatic mass bias.

The recent release of {\it Planck} data gives access to a full sky coverage of the thermal Sunyaev-Zel'dovich (tSZ) effect and of the cosmic microwave background (CMB) lensing potential ($\phi$). The cross-correlation of these two probes of the large-scale structures in the Universe is a powerful tool for testing cosmological models, especially in the context of the difference between galaxy clusters and CMB for the best-fitting cosmological parameters. However, the tSZ effect maps are highly contaminated by cosmic infra-red background (CIB) fluctuations. Unlike other astrophysical components, the spatial distribution of CIB varies with frequency. Thus it cannot be completely removed from a tSZ Compton parameter map, which is constructed from a linear combination of multiple frequency maps. We have estimated the contamination of the CIB-$\phi$ correlation in the tSZ-$\phi$ power-spectrum. We considered linear combinations that reconstruct the tSZ Compton parameter from {\it Planck} frequency maps. We conclude that even in an optimistic case, the CIB-$\phi$ contamination is significant with respect to the tSZ-$\phi$ signal itself. Consequently, We stress that tSZ-$\phi$ analyses that are based on Compton parameter maps are highly limited by the bias produced by CIB-$\phi$ contamination.

The thermal Sunyaev-Zel'dovich (tSZ) effect is produced by the interaction of cosmic microwave background (CMB) photons with the hot (a few keV) and diffuse gas of electrons inside galaxy clusters integrated along the line of sight. This effect produces a distortion of CMB blackbody emission law. This distortion law depends on the electronic temperature of the intra-cluster hot gas, $T_{e}$, through the so-called tSZ relativistic corrections. In this work, we have performed a statistical analysis of the tSZ spectral distortion on large galaxy cluster samples. We performed a stacking analysis for several electronic temperature bins, using both spectroscopic measurements of X-ray temperatures and a scaling relation between X-ray luminosities and electronic temperatures. We report the first high significance detection of the relativistic tSZ at a significance of 5.3 $\sigma$. We also demonstrate that the observed tSZ relativistic corrections are consistent with X-ray deduced temperatures. This measurement of the tSZ spectral law demonstrates that tSZ effect spectral distorsion can be used as a probe to measure galaxy cluster temperatures.

We present the first detection of the cross-correlation angular power spectrum between the thermal Sunyaev-Zel'dovich (tSZ) effect and the X-ray emission over the full sky. The tSZ effect and X-rays are produced by the same hot gas within groups and clusters of galaxies, which creates a naturally strong correlation between them that can be used to boost the joint signal and derive cosmological parameters. We computed the correlation between the ROSAT All Sky Survey in the 0.5-2 keV energy band and the tSZ effect reconstructed from six Planck all-sky frequency maps between 70 and 545 GHz. We detect a significant correlation over a wide range of angular scales. In the range $50<\ell < 2000$, the cross-correlation of X-rays to tSZ is detected at an overall significance of 28 $\sigma$. As part of our systematic study, we performed a multi-frequency modelling of the AGN contamination and the correlation between cosmic infra-red background and X-rays. Taking advantage of the strong dependence of the cross-correlation signal on the amplitude of the power spectrum, we constrained $\sigma_8 = 0.804 \pm 0.037$, where modelling uncertainties dominate statistical and systematic uncertainties. We also derived constraints on the mass indices of scaling relations between the halo mass and X-ray luminosity, L$_{500}$-M$_{500}$, and SZ signal, Y$_{500}$-M$_{500}$, $a_{\rm sz} + a_{\rm x} = 3.37 \pm 0.09$, and on the indices of the extra-redshift evolution, $\beta_{\rm sz} + \beta_{\rm x} = 0.4^{+0.4}_{-0.5}$.

Thermal Sunyaev-Zel'dovich (tSZ) effect and X-ray emission from galaxy clusters have been extensively used to constrain cosmological parameters. These constraints are highly sensitive to the relations between cluster masses and observables (tSZ and X-ray fluxes). The cross-correlation of tSZ and X-ray data is thus a powerful tool, in addition of tSZ and X-ray based analysis, to test our modeling of both tSZ and X-ray emission from galaxy clusters. We chose to explore this cross correlation as both emissions trace the hot gas in galaxy clusters and thus constitute one the easiest correlation that can be studied. We present a complete modeling of the cross correlation between tSZ effect and X-ray emission from galaxy clusters, and focuses on the dependencies with clusters scaling laws and cosmological parameters. We show that the present knowledge of cosmological parameters and scaling laws parameters leads to an uncertainties of 47\% on the overall normalization of the tSZ-X cross correlation power spectrum. We present the expected signal-to-noise ratio for the tSZ-X cross-correlation angular power spectrum considering the sensitivity of actual tSZ and X-ray surveys from {\it Planck}-like data and ROSAT. We demonstrate that this signal-to-noise can reach 31.5 in realistic situation, leading to a constraint on the amplitude of tSZ-X cross correlation up to 3.2\%, fifteen times better than actual modeling limitations. Consequently, used in addition to other probes of cosmological parameters and scaling relations, we show that the tSZ-X is a powerful probe to constrain scaling relations and cosmological parameters.

The high-significance measurement of large-scale structure signals enables testing the isotropy of the Universe. The measurement of cosmological parameters through the large-scale distribution of matter is now a mature domain. This approach is mainly limited by our knowledge of astrophysical processes that are used to observe the large-scale structure. However, when we assume that these astrophysical processes are the same across the Universe, then it is possible to tightly constrain the isotropy of cosmological parameters across the sky. Particularly the X-SZ cross-correlation has been shown to be a probe of the large scale structures that has a high signal-to-noise ratio and low bias. For this analysis, we used a localized measurement of the X-SZ cross-correlation as a test of the cosmological parameter isotropy. Using the scatter of the X-SZ cross-correlation across the sky, we derive cosmological constraints $\sigma_{8} \left( \Omega_{\rm m}/0.28 \right)^{0.34} = 0.78 \pm 0.02$ and tight isotropy constraints on the dark energy dipole $\Delta \Omega_{\Lambda} < 0.07$ at 95\% confidence level.

This paper presents a study of the ISW effect from the Planck 2015 temperature and polarization data release. The CMB is cross-correlated with different LSS tracers: the NVSS, SDSS and WISE catalogues, and the Planck 2015 lensing map. This cross-correlation yields a detection at $4\,\sigma$, where most of the signal-to-noise is due to the Planck lensing and NVSS. In fact, the ISW effect is detected only from the Planck data (through the ISW-lensing bispectrum) at $\approx 3\,\sigma$, which is similar to the detection level achieved by combining the cross-correlation signal coming from all the catalogues. The ISW signal allow us to detect $\Omega_\Lambda$ at more than $3\,\sigma$. This cross-correlation analysis is performed only with the Planck temperature data, since the polarization scales available in the 2015 release do not permit significant improvement of the CMB-LSS cross-correlation detectability. Nevertheless, polarization data is used to study the anomalously large ISW signal previously reported through the stacking of CMB features at the locations of known superstructures. We find that the current Planck polarization data do not exclude that this signal could be caused by the ISW effect. In addition, the stacking of the Planck lensing map on the locations of superstructures exhibits a positive cross-correlation with these large-scale structures. Finally, we have improved our previous reconstruction of the ISW temperature fluctuations by combining the information encoded in all the previously mentioned LSS tracers. In particular, we construct a map of the ISW secondary anisotropies and the corresponding uncertainties map, obtained from simulations. We also explore the reconstruction of the ISW anisotropies caused by the LSS traced by the 2MPZ survey by directly inverting the density field into the gravitational potential field.

The secondary cosmic microwave background (CMB) $B$-modes stem from the post-decoupling distortion of the polarization $E$-modes due to the gravitational lensing effect of large-scale structures. These lensing-induced $B$-modes constitute both a valuable probe of the dark matter distribution and an important contaminant for the extraction of the primary CMB $B$-modes from inflation. Planck provides accurate nearly all-sky measurements of both the polarization $E$-modes and the integrated mass distribution via the reconstruction of the CMB lensing potential. By combining these two data products, we have produced an all-sky template map of the lensing-induced $B$-modes using a real-space algorithm that minimizes the impact of sky masks. The cross-correlation of this template with an observed (primordial and secondary) $B$-mode map can be used to measure the lensing $B$-mode power spectrum at multipoles up to $2000$. In particular, when cross-correlating with the $B$-mode contribution directly derived from the Planck polarization maps, we obtain lensing-induced $B$-mode power spectrum measurement at a significance level of $12\,\sigma$, which agrees with the theoretical expectation derived from the Planck best-fit $\Lambda$CDM model. This unique nearly all-sky secondary $B$-mode template, which includes the lensing-induced information from intermediate to small ($10\lesssim \ell\lesssim 1000$) angular scales, is delivered as part of the Planck 2015 public data release. It will be particularly useful for experiments searching for primordial $B$-modes, such as BICEP2/Keck Array or LiteBIRD, since it will enable an estimate to be made of the lensing-induced contribution to the measured total CMB $B$-modes.

We present the most significant measurement of the cosmic microwave background (CMB) lensing potential to date (at a level of 40 sigma), using temperature and polarization data from the Planck 2015 full-mission release. Using a polarization-only estimator we detect lensing at a significance of 5 sigma. We cross-check the accuracy of our measurement using the wide frequency coverage and complementarity of the temperature and polarization measurements. Public products based on this measurement include an estimate of the lensing potential over approximately 70% of the sky, an estimate of the lensing potential power spectrum in bandpowers for the multipole range 40<L<400 and an associated likelihood for cosmological parameter constraints. We find good agreement between our measurement of the lensing potential power spectrum and that found in the best-fitting LCDM model based on the Planck temperature and polarization power spectra. Using the lensing likelihood alone we obtain a percent-level measurement of the parameter combination $\sigma_8 \Omega_m^{0.25} = 0.591\pm 0.021$. We combine our determination of the lensing potential with the E-mode polarization also measured by Planck to generate an estimate of the lensing B-mode. We show that this lensing B-mode estimate is correlated with the B-modes observed directly by Planck at the expected level and with a statistical significance of 10 sigma, confirming Planck's sensitivity to this known sky signal. We also correlate our lensing potential estimate with the large-scale temperature anisotropies, detecting a cross-correlation at the 3 sigma level, as expected due to dark energy in the concordance LCDM model.

We investigate the possible presence of diffuse radio emission in the intermediate redshift, massive cluster PLCK G285.0-23.7 (z=0.39, M_500 = 8.39 x 10^(14) M_Sun). Our 16cm-band ATCA observations of PLCK G285.0-23.7 allow us to reach a rms noise level of ~11 microJy/beam on the wide-band (1.1-3.1 GHz), full-resolution (~5 arcsec) image of the cluster, making it one of the deepest ATCA images yet published. We also re-image visibilities at lower resolution in order to achieve a better sensitivity to low-surface-brightness extended radio sources. We detect one of the lowest luminosity radio halos known at z>0.35, characterised by a slight offset from the well-studied 1.4 GHz radio power vs. cluster mass correlation. Similarly to most known radio-loud clusters (i.e. those hosting diffuse non-thermal sources), PLCK G285.0-23.7 has a disturbed dynamical state. Our analysis reveals a similarly elongated X-ray and radio morphology. While the size of the radio halo in PLCK G285.0-23.7 is smaller than lower redshift radio-loud clusters in the same mass range, it shows a similar correlation with the cluster virial radius, as expected in the framework of hierarchical structure formation.

The Planck full mission cosmic microwave background(CMB) temperature and E-mode polarization maps are analysed to obtain constraints on primordial non-Gaussianity(NG). Using three classes of optimal bispectrum estimators - separable template-fitting (KSW), binned, and modal - we obtain consistent values for the local, equilateral, and orthogonal bispectrum amplitudes, quoting as our final result from temperature alone fNL^local=2.5+\-5.7, fNL^equil=-16+\-70 and fNL^ortho=-34+\-33(68%CL). Combining temperature and polarization data we obtain fNL^local=0.8+\-5.0, fNL^equil=-4+\-43 and fNL^ortho=-26+\-21 (68%CL). The results are based on cross-validation of these estimators on simulations, are stable across component separation techniques, pass an extensive suite of tests, and are consistent with Minkowski functionals based measurements. The effect of time-domain de-glitching systematics on the bispectrum is negligible. In spite of these test outcomes we conservatively label the results including polarization data as preliminary, owing to a known mismatch of the noise model in simulations and the data. Beyond fNL estimates, we present model-independent reconstructions of the CMB bispectrum and derive constraints on early universe scenarios that generate NG, including general single-field and axion inflation, initial state modifications, parity-violating tensor bispectra, and directionally dependent vector models. We also present a wide survey of scale-dependent oscillatory bispectra, and we look for isocurvature NG. Our constraint on the local primordial trispectrum amplitude is gNL^local=(-9.0+\-7.7)x10^4 (68%CL), and we perform an analysis of additional trispectrum shapes. The global picture is one of consistency with the premises of the LambdaCDM cosmology, namely that the structure we observe today was sourced by adiabatic, passive, Gaussian, and primordial seed perturbations.[abridged]

This paper presents the Planck 2015 likelihoods, statistical descriptions of the 2-point correlations of CMB data, using the hybrid approach employed previously: pixel-based at $\ell<30$ and a Gaussian approximation to the distribution of spectra at higher $\ell$. The main improvements are the use of more and better processed data and of Planck polarization data, and more detailed foreground and instrumental models, allowing further checks and enhanced immunity to systematics. Progress in foreground modelling enables a larger sky fraction. Improvements in processing and instrumental models further reduce uncertainties. For temperature, we perform an analysis of end-to-end instrumental simulations fed into the data processing pipeline; this does not reveal biases from residual instrumental systematics. The $\Lambda$CDM cosmological model continues to offer a very good fit to Planck data. The slope of primordial scalar fluctuations, $n_s$, is confirmed smaller than unity at more than 5{\sigma} from Planck alone. We further validate robustness against specific extensions to the baseline cosmology. E.g., the effective number of neutrino species remains compatible with the canonical value of 3.046. This first detailed analysis of Planck polarization concentrates on E modes. At low $\ell$ we use temperature at all frequencies and a subset of polarization. The frequency range improves CMB-foreground separation. Within the baseline model this requires a reionization optical depth $\tau=0.078\pm0.019$, significantly lower than without high-frequency data for explicit dust monitoring. At high $\ell$ we detect residual errors in E, typically O($\mu$K$^2$); we recommend temperature alone as the high-$\ell$ baseline. Nevertheless, Planck high-$\ell$ polarization allows a separate determination of $\Lambda$CDM parameters consistent with those from temperature alone.

Context. Galaxy clusters are continuously growing through the accretion of matter in their outskirts. This process induces inhomogeneities in the gas density distribution (clumping) which need to be taken into account to recover the physical properties of the intracluster medium (ICM) at large radii. Aims. We studied the thermodynamic properties in the outskirts (R > R500) of the massive galaxy cluster Abell 2142 by combining the Sunyaev Zel'dovich (SZ) effect with the X-ray signal. Methods. We combined the SZ pressure profile measured by Planck with the XMM-Newton gas density profile to recover radial profiles of temperature, entropy and hydrostatic mass out to 2R500. We used a method that is insensitive to clumping to recover the gas density, and we compared the results with traditional X-ray measurement techniques. Results. When taking clumping into account, our joint SZ/X-ray entropy profile is consistent with the predictions from pure gravitational collapse, whereas a significant entropy flattening is found when the effect of clumping is neglected. The hydrostatic mass profile recovered using joint X-ray/SZ data agrees with that obtained from spectroscopic X-ray measurements and with mass reconstructions obtained through weak lensing and galaxy kinematics. Conclusions. We found that clumping can explain the entropy flattening observed by Suzaku in the outskirts of several clusters. When using a method insensitive to clumping for the reconstruction of the gas density, the thermodynamic properties of Abell 2142 are compatible with the assumption that the thermal gas pressure sustains gravity and that the entropy is injected at accretion shocks, with no need to evoke more exotic physics. Our results highlight the need for X-ray observations with sufficient spatial resolution, and large collecting area, to understand the processes at work in cluster outer regions.

We present results based on full-mission Planck observations of temperature and polarization anisotropies of the CMB. These data are consistent with the six-parameter inflationary LCDM cosmology. From the Planck temperature and lensing data, for this cosmology we find a Hubble constant, H0= (67.8 +/- 0.9) km/s/Mpc, a matter density parameter Omega_m = 0.308 +/- 0.012 and a scalar spectral index with n_s = 0.968 +/- 0.006. (We quote 68% errors on measured parameters and 95% limits on other parameters.) Combined with Planck temperature and lensing data, Planck LFI polarization measurements lead to a reionization optical depth of tau = 0.066 +/- 0.016. Combining Planck with other astrophysical data we find N_ eff = 3.15 +/- 0.23 for the effective number of relativistic degrees of freedom and the sum of neutrino masses is constrained to < 0.23 eV. Spatial curvature is found to be |Omega_K| < 0.005. For LCDM we find a limit on the tensor-to-scalar ratio of r <0.11 consistent with the B-mode constraints from an analysis of BICEP2, Keck Array, and Planck (BKP) data. Adding the BKP data leads to a tighter constraint of r < 0.09. We find no evidence for isocurvature perturbations or cosmic defects. The equation of state of dark energy is constrained to w = -1.006 +/- 0.045. Standard big bang nucleosynthesis predictions for the Planck LCDM cosmology are in excellent agreement with observations. We investigate annihilating dark matter and deviations from standard recombination, finding no evidence for new physics. The Planck results for base LCDM are in agreement with BAO data and with the JLA SNe sample. However the amplitude of the fluctuations is found to be higher than inferred from rich cluster counts and weak gravitational lensing. Apart from these tensions, the base LCDM cosmology provides an excellent description of the Planck CMB observations and many other astrophysical data sets.

Continuum spectra covering centimetre to submillimetre wavelengths are presented for a northern sample of 104 extragalactic radio sources, mainly active galactic nuclei, based on four-epoch Planck data. The nine Planck frequencies, from 30 to 857 GHz, are complemented by a set of simultaneous ground-based radio observations between 1.1 and 37 GHz. The single-survey Planck data confirm that the flattest high-frequency radio spectral indices are close to zero, indicating that the original accelerated electron energy spectrum is much harder than commonly thought, with power-law index around 1.5 instead of the canonical 2.5. The radio spectra peak at high frequencies and exhibit a variety of shapes. For a small set of low-z sources, we find a spectral upturn at high frequencies, indicating the presence of intrinsic cold dust. Variability can generally be approximated by achromatic variations, while sources with clear signatures of evolving shocks appear to be limited to the strongest outbursts.

(abridged) We discuss the Galactic foreground emission between 20 and 100GHz based on observations by Planck/WMAP. The Commander component-separation tool has been used to separate the various astrophysical processes in total intensity. Comparison with RRL templates verifies the recovery of the free-free emission along the Galactic plane. Comparison of the high-latitude Halpha emission with our free-free map shows residuals that correlate with dust optical depth, consistent with a fraction (~30%) of Halpha having been scattered by high-latitude dust. We highlight a number of diffuse spinning dust morphological features at high latitude. There is substantial spatial variation in the spinning dust spectrum, with the emission peak ranging from below 20GHz to more than 50GHz. There is a strong tendency for the spinning dust component near many prominent HII regions to have a higher peak frequency, suggesting that this increase in peak frequency is associated with dust in the photodissociation regions around the nebulae. The emissivity of spinning dust in these diffuse regions is of the same order as previous detections in the literature. Over the entire sky, the commander solution finds more anomalous microwave emission than the WMAP component maps, at the expense of synchrotron and free-free emission. This can be explained by the difficulty in separating multiple broadband components with a limited number of frequency maps. Future surveys (5-20GHz), will greatly improve the separation by constraining the synchrotron spectrum. We combine Planck/WMAP data to make the highest S/N ratio maps yet of the intensity of the all-sky polarized synchrotron emission at frequencies above a few GHz. Most of the high-latitude polarized emission is associated with distinct large-scale loops and spurs, and we re-discuss their structure...

Full-sky CMB maps from the 2015 Planck release allow us to detect departures from global isotropy on the largest scales. We present the first searches using CMB polarization for correlations induced by a non-trivial topology with a fundamental domain intersecting, or nearly intersecting, the last scattering surface (at comoving distance $\chi_{rec}$). We specialize to flat spaces with toroidal and slab topologies, finding that explicit searches for the latter are sensitive to other topologies with antipodal symmetry. These searches yield no detection of a compact topology at a scale below the diameter of the last scattering surface. The limits on the radius $R_i$ of the largest sphere inscribed in the topological domain (at log-likelihood-ratio $\Delta\ln{L}>-5$ relative to a simply-connected flat Planck best-fit model) are $R_i>0.97\chi_{rec}$ for the cubic torus and $R_i>0.56\chi_{rec}$ for the slab. The limit for the cubic torus from the matched-circles search is numerically equivalent, $R_i>0.97\chi_{rec}$ (99% CL) from polarisation data alone. We also perform a Bayesian search for a Bianchi VII$_h$ geometry. In the non-physical setting where the Bianchi cosmology is decoupled from the standard cosmology, Planck temperature data favour the inclusion of a Bianchi component. However, the cosmological parameters generating this pattern are in strong disagreement with those found from CMB anisotropy data alone. Fitting the induced polarization pattern for this model to Planck data requires an amplitude of $-0.1\pm0.04$ compared to +1 if the model were to be correct. In the physical setting where the Bianchi parameters are fit simultaneously with the standard cosmological parameters, we find no evidence for a Bianchi VII$_h$ cosmology and constrain the vorticity of such models to $(\omega/H)_0<7.6\times10^{-10}$ (95% CL). [Abridged]

Recent models for the large-scale Galactic magnetic fields in the literature have been largely constrained by synchrotron emission and Faraday rotation measures. We use three different but representative models to compare their predicted polarized synchrotron and dust emission with that measured by the Planck satellite. We first update these models to match the Planck synchrotron products using a common model for the cosmic-ray leptons. We discuss the impact on this analysis of the ongoing problems of component separation in the Planck microwave bands and of the uncertain cosmic-ray spectrum. In particular, the inferred degree of ordering in the magnetic fields is sensitive to these systematic uncertainties, and we further show the importance of considering the expected variations in the observables in addition to their mean morphology. We then compare the resulting simulated emission to the observed dust polarization and find that the dust predictions do not match the morphology in the Planck data but underpredict the dust polarization away from the plane. We modify one of the models to roughly match both observables at high latitudes by increasing the field ordering in the thin disc near the observer. Though this specific analysis is dependent on the component separation issues, we present the improved model as a proof of concept for how these studies can be advanced in future using complementary information from ongoing and planned observational projects.

We study the implications of Planck data for models of dark energy (DE) and modified gravity (MG), beyond the cosmological constant scenario. We start with cases where the DE only directly affects the background evolution, considering Taylor expansions of the equation of state, principal component analysis and parameterizations related to the potential of a minimally coupled DE scalar field. When estimating the density of DE at early times, we significantly improve present constraints. We then move to general parameterizations of the DE or MG perturbations that encompass both effective field theories and the phenomenology of gravitational potentials in MG models. Lastly, we test a range of specific models, such as k-essence, f(R) theories and coupled DE. In addition to the latest Planck data, for our main analyses we use baryonic acoustic oscillations, type-Ia supernovae and local measurements of the Hubble constant. We further show the impact of measurements of the cosmological perturbations, such as redshift-space distortions and weak gravitational lensing. These additional probes are important tools for testing MG models and for breaking degeneracies that are still present in the combination of Planck and background data sets. All results that include only background parameterizations are in agreement with LCDM. When testing models that also change perturbations (even when the background is fixed to LCDM), some tensions appear in a few scenarios: the maximum one found is \sim 2 sigma for Planck TT+lowP when parameterizing observables related to the gravitational potentials with a chosen time dependence; the tension increases to at most 3 sigma when external data sets are included. It however disappears when including CMB lensing.

Although infrared (IR) overall dust emission from clusters of galaxies has been statistically detected using data from the Infrared Astronomical Satellite (IRAS), it has not been possible to sample the spectral energy distribution (SED) of this emission over its peak, and thus to break the degeneracy between dust temperature and mass. By complementing the IRAS spectral coverage with Planck satellite data from 100 to 857 GHz, we provide new constraints on the IR spectrum of thermal dust emission in clusters of galaxies. We achieve this by using a stacking approach for a sample of several hundred objects from the Planck cluster sample; this procedure averages out fluctuations from the IR sky, allowing us to reach a significant detection of the faint cluster contribution. We also use the large frequency range probed by Planck, together with component-separation techniques, to remove the contamination from both cosmic microwave background anisotropies and the thermal Sunyaev-Zeldovich effect (tSZ) signal, which dominate below 353 GHz. By excluding dominant spurious signals or systematic effects, averaged detections are reported at frequencies between 353 and 5000 GHz. We confirm the presence of dust in clusters of galaxies at low and intermediate redshifts, yielding an SED with a shape similar to that of the Milky Way. Planck's beam does not allow us to investigate the detailed spatial distribution of this emission (e.g., whether it comes from intergalactic dust or simply the dust content of the cluster galaxies), but the radial distribution of the emission appears to follow that of the stacked SZ signal, and thus the extent of the clusters. The recovered SED allows us to constrain the dust mass responsible for the signal, as well as its temperature. We additionally explore the evolution of the IR emission as a function of cluster mass and redshift.

We compute and investigate four types of imprint of a stochastic background of primordial magnetic fields (PMFs) on the cosmic microwave background (CMB) anisotropies: the impact of PMFs on the CMB spectra; the effect on CMB polarization induced by Faraday rotation; the impact of PMFs on the ionization history; magnetically-induced non-Gaussianities; and the magnetically-induced breaking of statistical isotropy. Overall, Planck data constrain the amplitude of PMFs to less than a few nanogauss. In particular, individual limits coming from the analysis of the CMB angular power spectra, using the Planck likelihood, are $B_{1\,\mathrm{Mpc}}< 4.4$ nG (where $B_{1\,\mathrm{Mpc}}$ is the comoving field amplitude at a scale of 1 Mpc) at 95% confidence level, assuming zero helicity, and $B_{1\,\mathrm{Mpc}}< 5.6$ nG for a maximally helical field.For nearly scale-invariant PMFs we obtain $B_{1\,\mathrm{Mpc}}<2.0$ nG and $B_{1\,\mathrm{Mpc}}<0.9$ nG if the impact of PMFs on the ionization history of the Universe is included. From the analysis of magnetically-induced non-Gaussianity we obtain three different values, corresponding to three applied methods, all below 5 nG. The constraint from the magnetically-induced passive-tensor bispectrum is $B_{1\,\mathrm{Mpc}}< 2.8$ nG. A search for preferred directions in the magnetically-induced passive bispectrum yields $B_{1\,\mathrm{Mpc}}< 4.5$ nG, whereas the the compensated-scalar bispectrum gives $B_{1\,\mathrm{Mpc}}< 3$ nG. The analysis of the Faraday rotation of CMB polarization by PMFs uses the Planck power spectra in $EE$ and $BB$ at 70 GHz and gives $B_{1\,\mathrm{Mpc}}< 1380$ nG. In our final analysis, we consider the harmonic-space correlations produced by Alfv\'en waves, finding no significant evidence for the presence of these waves. Together, these results comprise a comprehensive set of constraints on possible PMFs with Planck data.

We present an updated description of the Planck Low Frequency Instrument (LFI) data processing pipeline, associated with the 2015 data release. We point out the places where our results and methods have remained unchanged since the 2013 paper and we highlight the changes made for the 2015 release, describing the products (especially timelines) and the ways in which they were obtained. We demonstrate that the pipeline is self-consistent (principally based on simulations) and report all null tests. For the first time, we present LFI maps in Stokes Q and U polarization. We refer to other related papers where more detailed descriptions of the LFI data processing pipeline may be found if needed.

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.

We present all-sky modelling of the high resolution Planck, IRAS, and WISE infrared (IR) observations using the physical dust model presented by Draine and Li in 2007 (DL). We study the performance and results of this model, and discuss implications for future dust modelling. The present work extends the DL dust modelling carried out on nearby galaxies using Herschel and Spitzer data to Galactic dust emission. We employ the DL dust model to generate maps of the dust mass surface density, the optical extinction Av, and the starlight intensity parametrized by Umin. The DL model reproduces the observed spectral energy distribution (SED) satisfactorily over most of the sky, with small deviations in the inner Galactic disk and in low ecliptic latitude areas. We compare the DL optical extinction Av for the diffuse interstellar medium with optical estimates for 2 10^5 quasi-stellar objects (QSOs) observed in the Sloan digital sky survey. The DL Av estimates are larger than those determined towards QSOs by a factor of about 2, which depends on Umin. The DL fitting parameter Umin, effectively determined by the wavelength where the SED peaks, appears to trace variations in the far-IR opacity of the dust grains per unit Av, and not only in the starlight intensity. To circumvent the model deficiency, we propose an empirical renormalization of the DL Av estimate, dependent of Umin, which compensates for the systematic differences found with QSO observations. This renormalization also brings into agreement the DL Av estimates with those derived for molecular clouds from the near-IR colours of stars in the 2 micron all sky survey. The DL model and the QSOs data are used to compress the spectral information in the Planck and IRAS observations for the diffuse ISM to a family of 20 SEDs normalized per Av, parameterized by Umin, which may be used to test and empirically calibrate dust models.

We present the results of approximately three years of observations of Planck Sunyaev-Zeldovich (SZ) sources with telescopes at the Canary Islands observatories as part of the general optical follow-up programme undertaken by the Planck collaboration. In total, 78 SZ sources are discussed. Deep-imaging observations were obtained for most of these sources; spectroscopic observations in either in long-slit or multi-object modes were obtained for many. We effectively used 37.5 clear nights. We found optical counterparts for 73 of the 78 candidates. This sample includes 53 spectroscopic redshift determinations, 20 of them obtained with a multi-object spectroscopic mode. The sample contains new redshifts for 27 Planck clusters that were not included in the first Planck SZ source catalogue (PSZ1).

We present the current accounting of systematic effect uncertainties for the Low Frequency Instrument (LFI) that are relevant to the 2015 release of the Planck cosmological results, showing the robustness and consistency of our data set, especially for polarization analysis. We use two complementary approaches: (i) simulations based on measured data and physical models of the known systematic effects; and (ii) analysis of difference maps containing the same sky signal ("null-maps"). The LFI temperature data are limited by instrumental noise. At large angular scales the systematic effects are below the cosmic microwave background (CMB) temperature power spectrum by several orders of magnitude. In polarization the systematic uncertainties are dominated by calibration uncertainties and compete with the CMB $E$-modes in the multipole range 10--20. Based on our model of all known systematic effects, we show that these effects introduce a slight bias of around $0.2\,\sigma$ on the reionization optical depth derived from the 70\,GHz $EE$ spectrum using the 30 and 353\,GHz channels as foreground templates. At 30\,GHz the systematic effects are smaller than the Galactic foreground at all scales in temperature and polarization, which allows us to consider this channel as a reliable template of synchrotron emission. We assess the residual uncertainties due to LFI effects on CMB maps and power spectra after component separation and show that these effects are smaller than the CMB amplitude at all scales. We also assess the impact on non-Gaussianity studies and find it to be negligible. Some residuals still appear in null maps from particular sky survey pairs, particularly at 30 GHz, suggesting possible straylight contamination due to an imperfect knowledge of the beam far sidelobes.

We test the statistical isotropy and Gaussianity of the cosmic microwave background (CMB) anisotropies using observations made by the Planck satellite. Our results are based mainly on the full Planck mission for temperature, but also include some polarization measurements. In particular, we consider the CMB anisotropy maps derived from the multi-frequency Planck data by several component-separation methods. For the temperature anisotropies, we find excellent agreement between results based on these sky maps over both a very large fraction of the sky and a broad range of angular scales, establishing that potential foreground residuals do not affect our studies. Tests of skewness, kurtosis, multi-normality, N-point functions, and Minkowski functionals indicate consistency with Gaussianity, while a power deficit at large angular scales is manifested in several ways, for example low map variance. The results of a peak statistics analysis are consistent with the expectations of a Gaussian random field. The "Cold Spot" is detected with several methods, including map kurtosis, peak statistics, and mean temperature profile. We thoroughly probe the large-scale dipolar power asymmetry, detecting it with several independent tests, and address the subject of a posteriori correction. Tests of directionality suggest the presence of angular clustering from large to small scales, but at a significance that is dependent on the details of the approach. We perform the first examination of polarization data, finding the morphology of stacked peaks to be consistent with the expectations of statistically isotropic simulations. Where they overlap, these results are consistent with the Planck 2013 analysis based on the nominal mission data and provide our most thorough view of the statistics of the CMB fluctuations to date.

The Planck catalogue of SZ sources limits itself to a significance threshold of 4.5 to ensure a low contamination rate by false cluster candidates. This means that only the most massive clusters at redshift z>0.5, and in particular z>0.7, are expected to enter into the catalogue, with a large number of systems in that redshift regime being expected around and just below that threshold. In this paper, we follow-up a sample of SZ sources from the Planck SZ catalogues from 2013 and 2015. In the latter maps, we consider detections around and at lower significance than the threshold adopted by the Planck Collaboration. To keep the contamination rate low, our 28 candidates are chosen to have significant WISE detections, in combination with non-detections in SDSS/DSS, which effectively selects galaxy cluster candidates at redshifts $z\gtrsim0.5$. By taking r- and z-band imaging with MegaCam@CFHT, we bridge the 4000A rest-frame break over a significant redshift range, thus allowing accurate redshift estimates of red-sequence cluster galaxies up to z~0.8. After discussing the possibility that an overdensity of galaxies coincides -by chance- with a Planck SZ detection, we confirm that 16 of the candidates have likely optical counterparts to their SZ signals, 13 (6) of which have an estimated redshift z>0.5 (z>0.7). The richnesses of these systems are generally lower than expected given the halo masses estimated from the Planck maps. However, when we follow a simplistic model to correct for Eddington bias in the SZ halo mass proxy, the richnesses are consistent with a reference mass-richness relation established for clusters detected at higher significance. This illustrates the benefit of an optical follow-up, not only to obtain redshift estimates, but also to provide an independent mass proxy that is not based on the same data the clusters are detected with, and thus not subject to Eddington bias.

We present a description of the pipeline used to calibrate the Planck Low Frequency Instrument (LFI) timelines into thermodynamic temperatures for the Planck 2015 data release, covering four years of uninterrupted operations. As in the 2013 data release, our calibrator is provided by the spin-synchronous modulation of the cosmic microwave background dipole, but we now use the orbital component, rather than adopting the Wilkinson Microwave Anisotropy Probe (WMAP) solar dipole. This allows our 2015 LFI analysis to provide an independent Solar dipole estimate, which is in excellent agreement with that of HFI and within $1\sigma$ (0.3% in amplitude) of the WMAP value. This 0.3% shift in the peak-to-peak dipole temperature from WMAP and a global overhaul of the iterative calibration code increases the overall level of the LFI maps by 0.45% (30 GHz), 0.64% (44 GHz), and 0.82% (70 GHz) in temperature with respect to the 2013 Planck data release, thus reducing the discrepancy with the power spectrum measured by WMAP. We estimate that the LFI calibration uncertainty is now at the level of 0.20% for the 70 GHz map, 0.26% for the 44 GHz map, and 0.35% for the 30 GHz map. We provide a detailed description of the impact of all the changes implemented in the calibration since the previous data release.

By looking at the kinetic Sunyaev-Zeldovich effect (kSZ) in Planck nominal mission data, we present a significant detection of baryons participating in large-scale bulk flows around central galaxies (CGs) at redshift $z\approx 0.1$. We estimate the pairwise momentum of the kSZ temperature fluctuations at the positions of the CGC (Central Galaxy Catalogue) samples extracted from Sloan Digital Sky Survey (DR7) data. For the foreground-cleaned maps, we find $1.8$-$2.5\sigma$ detections of the kSZ signal, which are consistent with the kSZ evidence found in individual Planck raw frequency maps, although lower than found in the WMAP-9yr W band ($3.3\sigma$). We further reconstruct the peculiar velocity field from the CG density field, and compute for the first time the cross-correlation function between kSZ temperature fluctuations and estimates of CG radial peculiar velocities. This correlation function yields a $3.0$-$3.7$$\sigma$ detection of the peculiar motion of extended gas on Mpc scales, in flows correlated up to distances of 80-100 $h^{-1}$ Mpc. Both the pairwise momentum estimates and kSZ temperature-velocity field correlation find evidence for kSZ signatures out to apertures of 8 arcmin and beyond, corresponding to a physical radius of $> 1$ Mpc, more than twice the mean virial radius of halos. This is consistent with the predictions from hydro simulations that most of the baryons are outside the virialized halos. We fit a simple model, in which the temperature-velocity cross-correlation is proportional to the signal seen in a semi-analytic model built upon N-body simulations, and interpret the proportionality constant as an "effective" optical depth to Thomson scattering. We find $\tau_T=(1.4\pm0.5)\times 10^{-4}$; the simplest interpretation of this measurement is that much of the gas is in a diffuse phase, which contributes little signal to X-ray or thermal SZ observations.

The Virgo cluster is the largest Sunyaev-Zeldovich (SZ) source in the sky, both in terms of angular size and total integrated flux. Planck's wide angular scale and frequency coverage, together with its high sensitivity, allow a detailed study of this large object through the SZ effect. Virgo is well resolved by Planck, showing an elongated structure, which correlates well with the morphology observed from X-rays, but extends beyond the observed X-ray signal. We find a good agreement between the SZ signal (or Compton paranmeter, y_c) observed by Planck and the expected signal inferred from X-ray observations and simple analytical models. Due to its proximity to us, the gas beyond the virial radius can be studied with unprecedented sensitivity by integrating the SZ signal over tens of square degrees. We study the signal in the outskirts of Virgo and compare it with analytical models and a constrained simulation of the environment of Virgo. Planck data suggest that significant amounts of low-density plasma surround Virgo out to twice the virial radius. We find the SZ signal in the outskirts of Virgo to be consistent with a simple model that extrapolates the inferred pressure at lower radii while assuming that the temperature stays in the keV range beyond the virial radius. The observed signal is also consistent with simulations and points to a shallow pressure profile in the outskirts of the cluster. This reservoir of gas at large radii can be linked with the hottest phase of the elusive warm/hot intergalactic medium. Taking the lack of symmetry of Virgo into account, we find that a prolate model is favoured by the combination of SZ and X-ray data, in agreement with predictions.

This paper presents the characterization of the in-flight beams, the beam window functions, and the associated uncertainties for the Planck Low Frequency Instrument (LFI). The structure of the paper is similar to that presented in the 2013 Planck release; the main differences concern the beam normalization and the delivery of the window functions to be used for polarization analysis. The in-flight assessment of the LFI main beams relies on measurements performed during observations of Jupiter. By stacking data from seven Jupiter transits, the main beam profiles are measured down to -25 dB at 30 and 44 GHz, and down to -30 dB at 70 GHz. It has been confirmed that the agreement between the simulated beams and the measured beams is better than 1% at each LFI frequency band (within the 20 dB contour from the peak, the rms values are 0.1% at 30 and 70 GHz; 0.2% at 44 GHz). Simulated polarized beams are used for the computation of the effective beam window functions. The error budget for the window functions is estimated from both main beam and sidelobe contributions, and accounts for the radiometer band shapes. The total uncertainties in the effective beam window functions are 0.7% and 1% at 30 and 44 GHz, respectively (at $\ell \approx 600$); and 0.5% at 70 GHz (at $\ell \approx 1000$).

The quest for a B-mode imprint from primordial gravity waves on the polarization of the cosmic microwave background (CMB) requires the characterization of foreground polarization from Galactic dust. We present a statistical study of the filamentary structure of the 353 GHz Planck Stokes maps at high Galactic latitude, relevant to the study of dust emission as a polarized foreground to the CMB. We filter the intensity and polarization maps to isolate filaments in the range of angular scales where the power asymmetry between E-modes and B-modes is observed. Using the Smoothed Hessian Major Axis Filament Finder, we identify 259 filaments at high Galactic latitude, with lengths larger or equal to 2\deg\ (corresponding to 3.5 pc in length for a typical distance of 100 pc). These filaments show a preferred orientation parallel to the magnetic field projected onto the plane of the sky, derived from their polarization angles. We present mean maps of the filaments in Stokes I, Q, U, E, and B, computed by stacking individual images rotated to align the orientations of the filaments. Combining the stacked images and the histogram of relative orientations, we estimate the mean polarization fraction of the filaments to be 11 %. Furthermore, we show that the correlation between the filaments and the magnetic field orientations may account for the E and B asymmetry and the $C_{\ell}^{TE}/C_{\ell}^{EE}$ ratio, reported in the power spectra analysis of the Planck 353 GHz polarization maps. Future models of the dust foreground for CMB polarization studies will need to take into account the observed correlation between the dust polarization and the structure of interstellar matter.

We use Planck data to detect the cross-correlation between the thermal Sunyaev-Zeldovich (tSZ) effect and the infrared emission from the galaxies that make up the the cosmic infrared background (CIB). We first perform a stacking analysis towards Planck-confirmed galaxy clusters. We detect infrared emission produced by dusty galaxies inside these clusters and demonstrate that the infrared emission is about 50% more extended than the tSZ effect. Modelling the emission with a Navarro--Frenk--White profile, we find that the radial profile concentration parameter is $c_{500} = 1.00^{+0.18}_{-0.15}$. This indicates that infrared galaxies in the outskirts of clusters have higher infrared flux than cluster-core galaxies. We also study the cross-correlation between tSZ and CIB anisotropies, following three alternative approaches based on power spectrum analyses: (i) using a catalogue of confirmed clusters detected in Planck data; (ii) using an all-sky tSZ map built from Planck frequency maps; and (iii) using cross-spectra between Planck frequency maps. With the three different methods, we detect the tSZ-CIB cross-power spectrum at significance levels of (i) 6 $\sigma$, (ii) 3 $\sigma$, and (iii) 4 $\sigma$. We model the tSZ-CIB cross-correlation signature and compare predictions with the measurements. The amplitude of the cross-correlation relative to the fiducial model is $A_{\rm tSZ-CIB}= 1.2\pm0.3$. This result is consistent with predictions for the tSZ-CIB cross-correlation assuming the best-fit cosmological model from Planck 2015 results along with the tSZ and CIB scaling relations.

We present the 8th Full Focal Plane simulation set (FFP8), deployed in support of the Planck 2015 results. FFP8 consists of 10 fiducial mission realizations reduced to 18144 maps, together with the most massive suite of Monte Carlo realizations of instrument noise and CMB ever generated, comprising $10^4$ mission realizations reduced to about $10^6$ maps. The resulting maps incorporate the dominant instrumental, scanning, and data analysis effects; remaining subdominant effects will be included in future updates. Generated at a cost of some 25 million CPU-hours spread across multiple high-performance-computing (HPC) platforms, FFP8 is used for the validation and verification of analysis algorithms, as well as their implementations, and for removing biases from and quantifying uncertainties in the results of analyses of the real data.

Planck has mapped the polarized dust emission over the whole sky, making it possible to trace the Galactic magnetic field structure that pervades the interstellar medium (ISM). We combine polarization data from Planck with rotation measure (RM) observations towards a massive star-forming region, the Rosette Nebula in the Monoceros molecular cloud, to study its magnetic field structure and the impact of an expanding HII region on the morphology of the field. We derive an analytical solution for the magnetic field, assumed to evolve from an initially uniform configuration following the expansion of ionized gas and the formation of a shell of swept-up ISM. From the RM data we estimate a mean value of the line-of-sight component of the magnetic field of about 3microG (towards the observer) in the Rosette Nebula, for a uniform electron density of about 12cm-3. The dust shell that surrounds the Rosette HII region is clearly observed in the Planck intensity map at 353 GHz, with a polarization signal significantly different from that of the local background when considered as a whole. The Planck observations constrain the plane-of-the-sky orientation of the magnetic field in the Rosette's parent molecular cloud to be mostly aligned with the large-scale field along the Galactic plane. The Planck data are compared with the analytical model, which predicts the mean polarization properties of a spherical and uniform dust shell for a given orientation of the field. This comparison leads to an upper limit of about 45degr on the angle between the line of sight and the magnetic field in the Rosette complex, for an assumed intrinsic dust polarization fraction of 4%. This field direction can reproduce the RM values detected in the ionized region if the magnetic field strength in the Monoceros molecular cloud is in the range 6.5--9microG.

Within ten nearby (d < 450 pc) Gould Belt molecular clouds we evaluate statistically the relative orientation between the magnetic field projected on the plane of sky, inferred from the polarized thermal emission of Galactic dust observed by Planck at 353 GHz, and the gas column density structures, quantified by the gradient of the column density, $N_H$. The selected regions, covering several degrees in size, are analyzed at an effective angular resolution of 10' FWHM, thus sampling physical scales from 0.4 to 40 pc in the nearest cloud. The column densities in the selected regions range from $N_H \approx 10^{21}$ to $10^{23}$ cm$^{-2}$, and hence they correspond to the bulk of the molecular clouds. The relative orientation is evaluated pixel by pixel and analyzed in bins of column density using the novel statistical tool called "Histogram of Relative Orientations". Throughout this study, we assume that the polarized emission observed by Planck at 353 GHz is representative of the projected morphology of the magnetic field in each region, i.e., we assume a constant dust grain alignment efficiency, independent of the local environment. Within most clouds we find that the relative orientation changes progressively with increasing $N_H$, from preferentially parallel or having no preferred orientation to preferentially perpendicular. In simulations of magnetohydrodynamic turbulence in molecular clouds this trend in relative orientation is a signature of Alfv\'enic or sub-Alfv\'enic turbulence, implying that the magnetic field is significant for the gas dynamics at the scales probed by Planck. We compare the deduced magnetic field strength with estimates we obtain from other methods and discuss the implications of the Planck observations for the general picture of molecular cloud formation and evolution.

The Andromeda Galaxy (M31) is one of a few galaxies that has sufficient angular size on the sky to be resolved by the Planck satellite. Planck has detected M31 in all of its frequency bands, and has mapped out the dust emission with the High Frequency Instrument, clearly resolving multiple spiral arms and sub-features. We examine the morphology of this long-wavelength dust emission as seen by Planck, including a study of its outermost spiral arms, and investigate the dust heating mechanism across M31. We find that dust dominating the longer wavelength emission ($\gtrsim 0.3\,$mm) is heated by the diffuse stellar population (as traced by 3.6$\,\mu$m emission), with the dust dominating the shorter wavelength emission heated by a mix of the old stellar population and star-forming regions (as traced by 24$\,\mu$m emission). We also fit spectral energy distributions (SEDs) for individual 5' pixels and quantify the dust properties across the galaxy, taking into account these different heating mechanisms, finding that there is a linear decrease in temperature with galactocentric distance for dust heated by the old stellar population, as would be expected, with temperatures ranging from around 22$\,$K in the nucleus to 14$\,$K outside of the 10$\,$kpc ring. Finally, we measure the integrated spectrum of the whole galaxy, which we find to be well-fitted with a global dust temperature of ($18.2\pm1.0$)$\,$K with a spectral index of $1.62\pm0.11$ (assuming a single modified blackbody), and a significant amount of free-free emission at intermediate frequencies of 20-60$\,$GHz, which corresponds to a star formation rate of around $0.12$M$_\odot\,$yr$^{-1}$. We find a $2.3\,\sigma$ detection of the presence of spinning dust emission, with a 30$\,$GHz amplitude of $0.7\pm0.3\,$Jy, which is in line with expectations from our Galaxy.

The Planck mission, thanks to its large frequency range and all-sky coverage, has a unique potential for systematically detecting the brightest, and rarest, submillimetre sources on the sky, including distant objects in the high-redshift Universe traced by their dust emission. A novel method, based on a component-separation procedure using a combination of Planck and IRAS data, has been applied to select the most luminous cold submm sources with spectral energy distributions peaking between 353 and 857GHz at 5' resolution. A total of 2151 Planck high-z source candidates (the PHZ) have been detected in the cleanest 26% of the sky, with flux density at 545GHz above 500mJy. Embedded in the cosmic infrared background close to the confusion limit, these high-z candidates exhibit colder colours than their surroundings, consistent with redshifts z>2, assuming a dust temperature of 35K and a spectral index of 1.5. First follow-up observations obtained from optical to submm have confirmed that this list consists of two distinct populations. A small fraction (around 3%) of the sources have been identified as strongly gravitationally lensed star-forming galaxies, which are amongst the brightest submm lensed objects (with flux density at 545GHz ranging from 350mJy up to 1Jy) at redshift 2 to 4. However, the vast majority of the PHZ sources appear as overdensities of dusty star-forming galaxies, having colours consistent with z>2, and may be considered as proto-cluster candidates. The PHZ provides an original sample, complementary to the Planck Sunyaev-Zeldovich Catalogue; by extending the population of the virialized massive galaxy clusters to a population of sources at z>1.5, the PHZ may contain the progenitors of today's clusters. Hence the PHZ opens a new window on the study of the early ages of structure formation, and the understanding of the intensively star-forming phase at high-z.