We report the direct detection of the kinetic Sunyaev-Zel'dovich effect in galaxy clusters with a 3.8 sigma significance level. The measurement is performed by stacking the Planck map at 217 GHz at the positions of galaxy clusters from the Wen-Han-Liu (WHL) catalog. To avoid the cancellation of positive and negative kSZ signals, we use the large scale distribution of the SDSS galaxies to estimate the peculiar velocities of the galaxy clusters along the line-of-sight and incorporate the sign in the (velocity weighted) stacking of the kSZ signals. Using this technique, we are able to measure the kSZ signal around galaxy clusters beyond 3R500. Assuming a standard beta-model, we also find that the gas fraction within R500 is f(gas,500) = 0.13 +- 0.03 for the clusters with the mass of M500 ~ 0.9 * 10^14 Msun/h, and is slightly lower than the universal baryon fraction. We compare this result to predictions from the Magneticum cosmological hydrodynamic simulations as well as other kSZ and X-ray measurements, where most of them show a lower gas fraction than the universal value for the same mass of clusters.

Galaxy clusters are a recent cosmological probe. The precision and accuracy of the cosmological parameters inferred from these objects are affected by the knowledge of cluster physics, entering the analysis through the mass-observable scaling relations, and the theoretical description of their mass and redshift distribution, modelled by the mass function. In this work, we forecast the impact of different modelling of these ingredients for clusters detected by future optical and near-IR surveys. We consider the standard cosmological scenario and the case with a time-dependent equation of state for dark energy. We analyse the effect of increasing accuracy on the scaling relation calibration, finding improved constraints on the cosmological parameters. This higher accuracy exposes the impact of the mass function evaluation, which is a subdominant source of systematics for current data. We compare two different evaluations for the mass function. In both cosmological scenarios, the use of different mass functions leads to biases in the parameter constraints. For the $\Lambda$CDM model, we find a $1.6 \, \sigma$ shift in the $(\Omega_m,\sigma_8)$ parameter plane and a discrepancy of $\sim 7 \, \sigma$ for the redshift evolution of the scatter of the scaling relations. For the scenario with a time-evolving dark energy equation of state, the assumption of different mass functions results in a $\sim 8 \, \sigma$ tension in the $w_0$ parameter. These results show the impact, and the necessity for a precise modelling, of the interplay between the redshift evolution of the mass function and of the scaling relations in the cosmological analysis of galaxy clusters.

We report new galaxy clusters previously unknown included in the PSZ1 catalogue. The results here presented have been achieved during the second year of a 2-year observational programme, the ITP13, developed at the Roque de los Muchachos Observatory (La Palma, Spain). We characterise 75 SZ sources with low SZ significance, SZ S/N$<5.32$ by performing deep optical imaging and spectroscopy in order to associate actual galaxy clusters to the SZ Planck source. We adopt robust criteria, based on the 2D-spatial distribution, richness and velocity dispersions to confirm actual optical counterparts up to $z<0.85$. At the end of the ITP13 observational programme, we study 256 SZ sources with $Dec \geq -15^{\circ}$ (212 of them completely unknown), finding optical counterparts for 152 SZ sources. The ITP13 validation programme has allowed us to update the PSZ1 purity, which is now more refined, increasing from 72\% to 83\% in the low SZ S/N regime. Our results are consistent with the predicted purity curve for the full PSZ1 catalogue and with the expected fraction of false detections caused by the non-Gaussian noise of foreground signals. Indeed, we find a strong correlation between the number of unconfirmed sources and the thermal emission of diffuse galactic dust at 857 GHz, thus increasing the fraction of false Planck SZ detections at low galactic latitudes.

Galaxy clusters are connected at their peripheries to the large scale structures by cosmic filaments that funnel accreting material. These filamentary structures are studied to investigate both environment-driven galaxy evolution and structure formation and evolution. In the present work, we probe in a statistical manner the azimuthal distribution of galaxies around clusters as a function of the cluster-centric distance, the cluster richness, and the galaxy activity (star-forming or passive).We perform a harmonic decomposition in large photometric galaxy catalogue around 6400 SDSS clusters with masses M >1e14 solar masses, in the redshift range of 0.1< z <0.3. The same analysis is performed on the mock galaxy catalogue from the light-cone of Magneticum hydrodynamical simulation. We use the multipole analysis to quantify asymmetries in the 2-D galaxy distribution. In the inner cluster regions at R <2 R500, we confirm that the galaxy distribution traces an ellipsoidal shape, which is more pronounced for richest clusters. In the clusters' outskirts (R= [2-8]R500), filamentary patterns are detected in harmonic space with a mean angular scale m_mean= 4.2+/-0.1. Massive clusters seem to have a larger number of connected filaments than low massive ones. We also find that passive galaxies appear to better trace the filamentary structures around clusters, even if the contribution of SF ones tend to increase with the cluster-centric distance, suggesting a gradient of galaxy activity in filaments around clusters.

We studied physical properties of matter in 24,544 filaments ranging from 30 to 100 Mpc in length, identified in the Sloan Digital Sky Survey (SDSS). We stacked the Comptonization y map produced by the Planck Collaboration around the filaments, excluding the resolved galaxy groups and clusters above a mass of ~3*10^13 Msun. We detected the thermal Sunyaev-Zel'dovich signal for the first time at a significance of 4.4 sigma in filamentary structures on such a large scale. We also stacked the Planck cosmic microwave background (CMB) lensing convergence map in the same manner and detected the lensing signal at a significance of 8.1 sigma. To estimate physical properties of the matter, we considered an isothermal cylindrical filament model with a density distribution following a beta-model (beta=2/3). Assuming that the gas distribution follows the dark matter distribution, we estimate that the central gas and matter overdensity and gas temperature are overdensity = (19.0 +27.3 -12.1) and temperature = (1.2 +- 0.4)*10^6 K, which results in a measured baryon fraction of (0.080 +0.116 -0.051) * Omega_b.

We present a statistical study of the filamentary structures of the Cosmic Web in the large hydro-dynamical simulations TNG300-1 and MAGNETICUM at redshift z=0. We focus on the radial distribution of the galaxy density around filaments detected using the Discrete Persistent Structure Extractor (DisPerSE). We show that the average profile of filaments presents an excess of galaxy density (> 5 sigma) up to radial distances of 27 Mpc from the core. The relation between galaxy density and the length of filaments is further investigated showing that short (L_f < 9 Mpc) and long (L_f > 20 Mpc) filaments are two statistically different populations. Short filaments are puffier, denser and more connected to massive objects, whereas long filaments are thinner, less dense and more connected to less massive structures. These two populations trace different environments and may correspond to bridges of matter between over-dense structures (short filaments) and to the cosmic filaments shaping the skeleton of the Cosmic Web (long filaments). Through MCMC explorations, we find that the density profiles of both short and long filaments can be described by the same empirical models (GNFW, beta-model, a single and a double power law) with different and distinct sets of parameters.

Detecting the large-scale structure of the Universe from the galaxy distribution and characterising its components is of fundamental importance in astrophysics but also a difficult task to achieve. This is due to the necessity of having wide-area spectroscopic redshift surveys to accurately measure galaxy positions in space while also covering large areas of the sky, as well as to the difficulty of creating algorithms to extract cosmic web structures (e.g. filaments). Moreover, these detections will be affected by systematic uncertainties, stemming both from the characteristics of the survey used (e.g. its completeness and coverage) and from the unique properties of the specific method adopted to detect the cosmic web (i.e. the assumptions it relies on and the free parameters it may employ). For these reasons, the creation of new catalogues of cosmic web features on wide sky areas is important, as this allows users to have at their disposal a well understood sample of structures, whose systematic uncertainties have been thoroughly investigated. In this paper we present the filament catalogues created using the Discrete Persistent Structure Extractor (DisPerSE) tool in the Sloan Digital Sky Survey (SDSS) and we fully characterise them in terms of their dependence on the choice of parameters pertaining to the algorithm, as well as with respect to several systematic issues, which may arise in the skeleton due to the properties of the galaxy distribution (such as Finger-of-God redshift distortions and defects of the density field due to the boundaries of the survey). We create several fully characterised filament catalogues, which we make available to the scientific community upon request.

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.

The role played by the large-scale structures in the galaxy evolution is not quite well understood yet. In this study, we investigate properties of galaxy in the range 0.1<z<0.3 from a value-added version of the WISExSCOS catalogue around cosmic filaments detected with DisPerSE. We have fitted a profile of galaxy over-density around cosmic filaments and found a typical radius of r_m = 7.5+-0.2 Mpc. We have measured an excess of passive galaxies near the filament's spine, higher than the excess of transitioning and active galaxies. We have also detected SFR and Mstar gradients pointing towards the filament's spine. We have investigated this result and found an Mstar gradient for each type of galaxies: active, transitioning, and passive, and a positive SFR gradient for passive galaxies. We also link the galaxy properties and the gas content in the Cosmic Web. To do so, we have investigated the quiescent fraction fQ profile of galaxies around the cosmic filaments. Based on recent studies about the effect of the gas and of the Cosmic Web on galaxy properties, we have modeled fQ with a beta model of gas pressure. The slope obtained here, beta=0.54+-0.18, is compatible with the scenario of projected isothermal gas in hydrostatic equilibrium (beta=2/3), and with the profiles of gas fitted in SZ.

We construct the mean thermal Sunyaev-Zel'dovich (tSZ) Comptonization y profile around Luminous Red Galaxies (LRGs) in the redshift range 0.16 < z < 0.47 from the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) using the Planck y map. The mean central tSZ signal for the full sample is y ~ 1.8 * 10^(-7) and we detect tSZ emission out to ~30 arcmin, which is well beyond the 10 arcmin angular resolution of the y map and well beyond the virial radii of the LRGs. We compare the measured profile with predictions from the cosmo-OWLS suite of cosmological hydrodynamical simulations. This comparison agrees well for models that include feedback from active galactic nuclei (AGN), but not with hydrodynamic models without this energetic feedback mechanism. This suggests that an additional heating mechanism is required over SNe feedback and star formation to explain the y data profile. We also compare our results with predictions based on the halo model with a universal pressure profile (UPP) giving the y signal. The predicted profile is consistent with the data, but only if we account for the clustering of haloes via a two-halo term and if halo masses are estimated using the mean stellar-to-halo mass (SHM) relation of Coupon et al. (2015) or Wang et al.(2016) estimated from gravitational lensing measurements. We also discuss the importance of scatter in the SHM relation on the model predictions.

Numerical simulations and observations show that galaxies are not uniformly distributed in the universe but rather spread on a filamentary structure. In this large scale pattern, highly dense regions are linked together by bridges and walls, all of them surrounded by vast nearly-empty areas. While nodes of the network are widely studied in the literature, simulations indicate that half of the mass budget comes from a more diffuse part of the network made of filaments. In the context of recent and upcoming large galaxy surveys, it becomes essential to identify and classify features of the Cosmic Web in an automatic way to study their physical properties and the impact of the cosmic environment on galaxies and their evolution. In this work, we propose a new approach to automatically retrieve the underlying filamentary structure from a 2D or 3D galaxy distribution using graph theory and the assumption that paths linking galaxies together with the minimum total length highlight the underlying distribution. To obtain a smoothed version of this topological prior, we embed it in a Gaussian mixtures framework. In addition to a geometrical description of the pattern, a bootstrap-like estimate of these regularized minimum spanning trees allows to obtain a map characterising the frequency at which an area of the domain is crossed. Using distribution of halos derived from numerical simulations, we show that the proposed method is able to recover the filamentary pattern in a 2D or 3D distribution of points with noise and outliers robustness with few and comprehensible parameters.

The Cosmic Web is a complex network of filaments, walls and voids that represent the largest structures in the Universe. In this network, which is the direct result of structure formation, galaxy clusters occupy central positions as the nodes, connected through the filaments. In this work, we investigate the position in the Cosmic Web of one of the most known and best studied clusters of galaxies, the Coma cluster. We make use of the Sloan Digital Sky Survey Data Release 7 (SDSS DR7) Main Galaxy Sample and of the Discrete Persistent Structure Extractor (DisPerSE) to detect large scale filaments around the Coma cluster and we analyze the properties of the Cosmic Web. We study the network of filaments around Coma in a 75 Mpc radius region. We find that the Coma cluster has a median connectivity of 2.5, in agreement with measurements from clusters of similar mass in the literature, as well as with what expected from numerical simulations. It is indeed connected to 3 secure filaments which connect Coma to Abell 1367 and to several other clusters in the field. The location of these filaments in the vicinity of Coma is consistent with features detected in the X-ray, as well as the likely direction of infall of galaxies, as for example NGC 4839. The overall picture that emerges of the Coma cluster is that of a highly connected structure occupying a central position as a dense node of the Cosmic Web. We also find a tentative detection, at 2.1$\sigma$ significance, of the filaments in the SZ signal.

We search the Planck data for a thermal Sunyaev-Zel'dovich (tSZ) signal due to gas filaments between pairs of Luminous Red Galaxies (LRG's) taken from the Sloan Digital Sky Survey Data Release 12 (SDSS/DR12). We identify $\sim$260,000 LRG pairs in the DR12 catalog that lie within 6-10 $h^{-1} \mathrm{Mpc}$ of each other in tangential direction and within 6 $h^{-1} \mathrm{Mpc}$ in radial direction. We stack pairs by rotating and scaling the angular positions of each LRG so they lie on a common reference frame, then we subtract a circularly symmetric halo from each member of the pair to search for a residual signal between the pair members. We find a statistically significant (5.3$\sigma$) signal between LRG pairs in the stacked data with a magnitude $\Delta y = (1.31 \pm 0.25) \times 10^{-8}$. The uncertainty is estimated from two Monte Carlo null tests which also establish the reliability of our analysis. Assuming a simple, isothermal, cylindrical filament model of electron over-density with a radial density profile proportional to $r_c/r$ (as determined from simulations), where $r$ is the perpendicular distance from the cylinder axis and $r_c$ is the core radius of the density profile, we constrain the product of over-density and filament temperature to be $\delta_c \times (T_{\rm e}/10^7 \, {\rm K}) \times (r_c/0.5h^{-1} \, {\rm Mpc}) = 2.7 \pm 0.5$. To our knowledge, this is the first detection of filamentary gas at over-densities typical of cosmological large-scale structure. We compare our result to the BAHAMAS suite of cosmological hydrodynamic simulations (McCarthy et al. 2017) and find a slightly lower, but marginally consistent Comptonization excess, $\Delta y = (0.84 \pm 0.24) \times 10^{-8}$.

We present cosmological parameter results from the final full-mission Planck measurements of the CMB anisotropies. We find good consistency with the standard spatially-flat 6-parameter $\Lambda$CDM cosmology having a power-law spectrum of adiabatic scalar perturbations (denoted "base $\Lambda$CDM" in this paper), from polarization, temperature, and lensing, separately and in combination. A combined analysis gives dark matter density $\Omega_c h^2 = 0.120\pm 0.001$, baryon density $\Omega_b h^2 = 0.0224\pm 0.0001$, scalar spectral index $n_s = 0.965\pm 0.004$, and optical depth $\tau = 0.054\pm 0.007$ (in this abstract we quote $68\,\%$ confidence regions on measured parameters and $95\,\%$ on upper limits). The angular acoustic scale is measured to $0.03\,\%$ precision, with $100\theta_*=1.0411\pm 0.0003$. These results are only weakly dependent on the cosmological model and remain stable, with somewhat increased errors, in many commonly considered extensions. Assuming the base-$\Lambda$CDM cosmology, the inferred late-Universe parameters are: Hubble constant $H_0 = (67.4\pm 0.5)$km/s/Mpc; matter density parameter $\Omega_m = 0.315\pm 0.007$; and matter fluctuation amplitude $\sigma_8 = 0.811\pm 0.006$. We find no compelling evidence for extensions to the base-$\Lambda$CDM model. Combining with BAO we constrain the effective extra relativistic degrees of freedom to be $N_{\rm eff} = 2.99\pm 0.17$, and the neutrino mass is tightly constrained to $\sum m_\nu< 0.12$eV. The CMB spectra continue to prefer higher lensing amplitudes than predicted in base -$\Lambda$CDM at over $2\,\sigma$, which pulls some parameters that affect the lensing amplitude away from the base-$\Lambda$CDM model; however, this is not supported by the lensing reconstruction or (in models that also change the background geometry) BAO data. (Abridged)

This paper discusses the science case for a sensitive spectro-polarimetric survey of the microwave sky. Such a survey would provide a tomographic and dynamic census of the three-dimensional distribution of hot gas, velocity flows, early metals, dust, and mass distribution in the entire Hubble volume, exploit CMB temperature and polarisation anisotropies down to fundamental limits, and track energy injection and absorption into the radiation background across cosmic times by measuring spectral distortions of the CMB blackbody emission. In addition to its exceptional capability for cosmology and fundamental physics, such a survey would provide an unprecedented view of microwave emissions at sub-arcminute to few-arcminute angular resolution in hundreds of frequency channels, a data set that would be of immense legacy value for many branches of astrophysics. We propose that this survey be carried-out with a large space mission featuring a broad-band polarised imager and a moderate resolution spectro-imager at the focus of a 3.5m aperture telescope actively cooled to about 8K, complemented with absolutely-calibrated Fourier Transform Spectrometer modules observing at degree-scale angular resolution in the 10-2000 GHz frequency range. We propose two observing modes: a survey mode to map the entire sky as well as a few selected wide fields, and an observatory mode for deeper observations of regions of specific interest.

Voyage 2050 White Paper highlighting the unique science opportunities using spectral distortions of the cosmic microwave background (CMB). CMB spectral distortions probe many processes throughout the history of the Universe. Precision spectroscopy, possible with existing technology, would provide key tests for processes expected within the cosmological standard model and open an enormous discovery space to new physics. This offers unique scientific opportunities for furthering our understanding of inflation, recombination, reionization and structure formation as well as dark matter and particle physics. A dedicated experimental approach could open this new window to the early Universe in the decades to come, allowing us to turn the long-standing upper distortion limits obtained with COBE/FIRAS some 25 years ago into clear detections of the expected standard distortion signals.

This paper describes the 2018 Planck CMB likelihoods, following a hybrid approach similar to the 2015 one, with different approximations at low and high multipoles, and implementing several methodological and analysis refinements. With more realistic simulations, and better correction and modelling of systematics, we can now make full use of the High Frequency Instrument polarization data. The low-multipole 100x143 GHz EE cross-spectrum constrains the reionization optical-depth parameter $\tau$ to better than 15% (in combination with with the other low- and high-$\ell$ likelihoods). We also update the 2015 baseline low-$\ell$ joint TEB likelihood based on the Low Frequency Instrument data, which provides a weaker $\tau$ constraint. At high multipoles, a better model of the temperature-to-polarization leakage and corrections for the effective calibrations of the polarization channels (polarization efficiency or PE) allow us to fully use the polarization spectra, improving the constraints on the $\Lambda$CDM parameters by 20 to 30% compared to TT-only constraints. Tests on the modelling of the polarization demonstrate good consistency, with some residual modelling uncertainties, the accuracy of the PE modelling being the main limitation. Using our various tests, simulations, and comparison between different high-$\ell$ implementations, we estimate the consistency of the results to be better than the 0.5$\sigma$ level. Minor curiosities already present before (differences between $\ell$<800 and $\ell$>800 parameters or the preference for more smoothing of the $C_\ell$ peaks) are shown to be driven by the TT power spectrum and are not significantly modified by the inclusion of polarization. Overall, the legacy Planck CMB likelihoods provide a robust tool for constraining the cosmological model and represent a reference for future CMB observations. (Abridged)

We present measurements of the cosmic microwave background (CMB) lensing potential using the final $\textit{Planck}$ 2018 temperature and polarization data. We increase the significance of the detection of lensing in the polarization maps from $5\,\sigma$ to $9\,\sigma$. Combined with temperature, lensing is detected at $40\,\sigma$. We present an extensive set of tests of the robustness of the lensing-potential power spectrum, and construct a minimum-variance estimator likelihood over lensing multipoles $8 \le L \le 400$. We find good consistency between lensing constraints and the results from the $\textit{Planck}$ CMB power spectra within the $\rm{\Lambda CDM}$ model. Combined with baryon density and other weak priors, the lensing analysis alone constrains $\sigma_8 \Omega_{\rm m}^{0.25}=0.589\pm 0.020$ ($1\,\sigma$ errors). Also combining with baryon acoustic oscillation (BAO) data, we find tight individual parameter constraints, $\sigma_8=0.811\pm0.019$, $H_0=67.9_{-1.3}^{+1.2}\,\text{km}\,\text{s}^{-1}\,\rm{Mpc}^{-1}$, and $\Omega_{\rm m}=0.303^{+0.016}_{-0.018}$. Combining with $\textit{Planck}$ CMB power spectrum data, we measure $\sigma_8$ to better than $1\,\%$ precision, finding $\sigma_8=0.811\pm 0.006$. We find consistency with the lensing results from the Dark Energy Survey, and give combined lensing-only parameter constraints that are tighter than joint results using galaxy clustering. Using $\textit{Planck}$ cosmic infrared background (CIB) maps we make a combined estimate of the lensing potential over $60\,\%$ of the sky with considerably more small-scale signal. We demonstrate delensing of the $\textit{Planck}$ power spectra, detecting a maximum removal of $40\,\%$ of the lensing-induced power in all spectra. The improvement in the sharpening of the acoustic peaks by including both CIB and the quadratic lensing reconstruction is detected at high significance (abridged).

Using precise galaxy stellar mass function measurements in the COSMOS field we determine the stellar-to-halo mass relationship (SHMR) using a parametric abundance matching technique. The unique combination of size and highly complete stellar mass estimates in COSMOS allows us to determine the SHMR over a wide range of halo masses from $z\sim0.2$ to $z\sim5$. At $z\sim 0.2$ the ratio of stellar-to-halo mass content peaks at a characteristic halo mass $M_{\rm h} =10^{12} M_\odot$ and declines at higher and lower halo masses. This characteristic halo mass increases with redshift reaching $M_{\rm h} =10^{12.5} M_\odot$ at $z\sim2.3$ and remaining flat up to $z=4$. We considered the principal sources of uncertainty in our stellar mass measurements and also the variation in halo mass estimates in the literature. We show that our results are robust to these sources of uncertainty and explore likely explanation for differences between our results and those published in the literature. The steady increase in characteristic halo mass with redshift points to a scenario where cold gas inflows become progressively more important in driving star-formation at high redshifts but larger samples of massive galaxies are needed to rigorously test this hypothesis.

Galaxy clusters observed through the thermal Sunyaev-Zeldovich (tSZ) effect are a recent cosmological probe. The precision on the cosmological constraints is affected mainly by the current knowledge of cluster physics, which enters the analysis through the scaling relations. Here we aim to study one of the most important sources of systematic uncertainties, the mass bias, $b$. We have analysed the effects of a mass-redshift dependence, adopting a power-law parametrisation. We applied this parametrisation to the combination of tSZ number counts and power spectrum, finding a hint of redshift dependence that leads to a decreasing value of the mass bias for higher redshift. We tested the robustness of our results for different mass bias calibrations and a discrete redshift dependence. We find our results to be dependent on the clusters sample that we are considering, in particular obtaining an inverse (decreasing) redshift dependence when neglecting $z<0.2$ clusters. We analysed the effects of this parametrisation on the combination of cosmic microwave background (CMB) primary anisotropies and tSZ galaxy clusters. We find a preferred constant value of mass bias, having $(1-b) =0.62 \pm 0.05$. The corresponding value of $b$ is too high with respect to weak lensing and numerical simulations estimations. Therefore we conclude that this mass-redshift parametrisation does not help in solving the remaining discrepancy between CMB and tSZ clusters observations.

Using a thermal Sunyaev-Zel'dovich (tSZ) signal, we search for hot gas in superclusters identified using the Sloan Digital Sky Survey Data Release 7 (SDSS/DR7) galaxies. We stack a Comptonization y map produced by the Planck Collaboration around the superclusters and detect the tSZ signal at a significance of 6.4 sigma. We further search for an intercluster component of gas in the superclusters. For this, we remove the intracluster gas in the superclusters by masking all galaxy groups/clusters detected by the Planck tSZ, ROSAT X-ray, and SDSS optical surveys down to a total mass of 10^13 Msun. We report the first detection of intercluster gas in superclusters with y = (3.5 +- 1.4) * 10^(-8) at a significance of 2.5 sigma. Assuming a simple isothermal and flat density distribution of intercluster gas over superclusters, the estimated baryon density is (Omega_gas / Omega_b) * (T_e/(8*10^6 K)) = 0.067 +- 0.006 +- 0.025. This quantity is inversely proportional to the temperature, therefore taking values from simulations and observations, we find that the gas density in superclusters may account for 17 - 52 % of missing baryons at low redshifts. A better understanding of the physical state of gas in the superclusters is required to accurately estimate the contribution of our measurements to missing baryons.

Following the pioneering observations with COBE in the early 1990s, studies of the cosmic microwave background (CMB) have focused on temperature and polarization anisotropies. CMB spectral distortions - tiny departures of the CMB energy spectrum from that of a perfect blackbody - provide a second, independent probe of fundamental physics, with a reach deep into the primordial Universe. The theoretical foundation of spectral distortions has seen major advances in recent years, which highlight the immense potential of this emerging field. Spectral distortions probe a fundamental property of the Universe - its thermal history - thereby providing additional insight into processes within the cosmological standard model (CSM) as well as new physics beyond. Spectral distortions are an important tool for understanding inflation and the nature of dark matter. They shed new light on the physics of recombination and reionization, both prominent stages in the evolution of our Universe, and furnish critical information on baryonic feedback processes, in addition to probing primordial correlation functions at scales inaccessible to other tracers. In principle the range of signals is vast: many orders of magnitude of discovery space could be explored by detailed observations of the CMB energy spectrum. Several CSM signals are predicted and provide clear experimental targets, some of which are already observable with present-day technology. Confirmation of these signals would extend the reach of the CSM by orders of magnitude in physical scale as the Universe evolves from the initial stages to its present form. The absence of these signals would pose a huge theoretical challenge, immediately pointing to new physics.

We present 353 GHz full-sky maps of the polarization fraction $p$, angle $\psi$, and dispersion of angles $S$ of Galactic dust thermal emission produced from the 2018 release of Planck data. We confirm that the mean and maximum of $p$ decrease with increasing $N_H$. The uncertainty on the maximum polarization fraction, $p_\mathrm{max}=22.0$% at 80 arcmin resolution, is dominated by the uncertainty on the zero level in total intensity. The observed inverse behaviour between $p$ and $S$ is interpreted with models of the polarized sky that include effects from only the topology of the turbulent Galactic magnetic field. Thus, the statistical properties of $p$, $\psi$, and $S$ mostly reflect the structure of the magnetic field. Nevertheless, we search for potential signatures of varying grain alignment and dust properties. First, we analyse the product map $S \times p$, looking for residual trends. While $p$ decreases by a factor of 3--4 between $N_H=10^{20}$ cm$^{-2}$ and $N_H=2\times 10^{22}$ cm$^{-2}$, $S \times p$ decreases by only about 25%, a systematic trend observed in both the diffuse ISM and molecular clouds. Second, we find no systematic trend of $S \times p$ with the dust temperature, even though in the diffuse ISM lines of sight with high $p$ and low $S$ tend to have colder dust. We also compare Planck data with starlight polarization in the visible at high latitudes. The agreement in polarization angles is remarkable. Two polarization emission-to-extinction ratios that characterize dust optical properties depend only weakly on $N_H$ and converge towards the values previously determined for translucent lines of sight. We determine an upper limit for the polarization fraction in extinction of 13%, compatible with the $p_\mathrm{max}$ observed in emission. These results provide strong constraints for models of Galactic dust in diffuse gas.

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.

Recent years have brought strong observational evidences for the standard LCDM cosmological model. Cosmic microwave background (CMB) anisotropy and large scale structure (LSS) probes do not favour any extensions of the standard model. Nevertheless, in this framework, the preferred cosmological parameters may differ from probe to probe, from experiment to experiment. This is the well known case of the tension between CMB and Sunyaev Zel'dovich (SZ) galaxy clusters (GC) from Planck. In 2013, the Planck team has shown that the preferred matter content ({\Omega}M) and density fluctuation power spectrum amplitude ({\sigma}8), the two main cosmological parameters probed by the galaxy cluster number count, are different in the CMB analyses and in the SZ cluster analyses at more than 2 sigmas (a result confirmed in subsequent analyses). We present the results of our new analysis using more recent measurements of the CMB, SZ clusters and SZ power spectrum of 2016 and show that the tension on ({\Omega}M,{\sigma}8) is mostly releaved. The lower value of the reionisation optical depth and thus of {\sigma}8 in the recent Planck studies is the main reason. We also show that basic extensions of the standard model (massive neutrinos or non-lambda dark energy) do not help improving the agreement between the probes. In order to fully reconcile SZ clusters with CMB best model, the mass of the galaxy clusters derived from hydrostatic equilibrium should be 40% lower than the true mass. While current numerical simulations and weak lensing measurements agree for a mass bias of 20%, investigations are still going on to explain such disagreement on the mass bias. We show that considering a mass bias evolving with redshift or mass does not help in eliminating the discrepancy.

Star-formation activity is a key property to probe the structure formation and hence characterise the large-scale structures of the universe. This information can be deduced from the star formation rate (SFR) and the stellar mass (Mstar), both of which, but especially the SFR, are very complex to estimate. Determining these quantities from UV, optical, or IR luminosities relies on complex modeling and on priors on galaxy types. We propose a method based on the machine-learning algorithm Random Forest to estimate the SFR and the Mstar of galaxies at redshifts in the range 0.01<z<0.3, independent of their type. The machine-learning algorithm takes as inputs the redshift, WISE luminosities, and WISE colours in near-IR, and is trained on spectra-extracted SFR and Mstar from the SDSS MPA-JHU DR8 catalogue as outputs. We show that our algorithm can accurately estimate SFR and Mstar with scatters of sigma_SFR=0.38 dex and sigma_Mstar=0.16 dex for SFR and stellar mass, respectively, and that it is unbiased with respect to redshift or galaxy type. The full-sky coverage of the WISE satellite allows us to characterise the star-formation activity of all galaxies outside the Galactic mask with spectroscopic redshifts in the range 0.01<z<0.3. The method can also be applied to photometric-redshift catalogues, with best scatters of sigma_SFR=0.42 dex and sigma_Mstar=0.24 dex obtained in the redshift range 0.1<z<0.3.

Using the ${\it Planck}$ full-mission data, we present a detection of the temperature (and therefore velocity) dispersion due to the kinetic Sunyaev-Zeldovich (kSZ) effect from clusters of galaxies. To suppress the primary CMB and instrumental noise we derive a matched filter and then convolve it with the ${\it Planck}$ foreground-cleaned `${\tt 2D-ILC\,}$' maps. By using the Meta Catalogue of X-ray detected Clusters of galaxies (MCXC), we determine the normalized ${\it rms}$ dispersion of the temperature fluctuations at the positions of clusters, finding that this shows excess variance compared with the noise expectation. We then build an unbiased statistical estimator of the signal, determining that the normalized mean temperature dispersion of $1526$ clusters is $\langle \left(\Delta T/T \right)^{2} \rangle = (1.64 \pm 0.48) \times 10^{-11}$. However, comparison with analytic calculations and simulations suggest that around $0.7\,\sigma$ of this result is due to cluster lensing rather than the kSZ effect. By correcting this, the temperature dispersion is measured to be $\langle \left(\Delta T/T \right)^{2} \rangle = (1.35 \pm 0.48) \times 10^{-11}$, which gives a detection at the $2.8\,\sigma$ level. We further convert uniform-weight temperature dispersion into a measurement of the line-of-sight velocity dispersion, by using estimates of the optical depth of each cluster (which introduces additional uncertainty into the estimate). We find that the velocity dispersion is $\langle v^{2} \rangle =(123\,000 \pm 71\,000)\,({\rm km}\,{\rm s}^{-1})^{2}$, which is consistent with findings from other large-scale structure studies, and provides direct evidence of statistical homogeneity on scales of $600\,h^{-1}{\rm Mpc}$. Our study shows the promise of using cross-correlations of the kSZ effect with large-scale structure in order to constrain the growth of structure.

This paper presents the High Frequency Instrument (HFI) data processing procedures for the Planck 2018 release. Major improvements in mapmaking have been achieved since the previous 2015 release. They enabled the first significant measurement of the reionization optical depth parameter using HFI data. This paper presents an extensive analysis of systematic effects, including the use of simulations to facilitate their removal and characterize the residuals. The polarized data, which presented a number of known problems in the 2015 Planck release, are very significantly improved. Calibration, based on the CMB dipole, is now extremely accurate and in the frequency range 100 to 353 GHz reduces intensity-to-polarization leakage caused by calibration mismatch. The Solar dipole direction has been determined in the three lowest HFI frequency channels to within one arc minute, and its amplitude has an absolute uncertainty smaller than $0.35\mu$K, an accuracy of order $10^{-4}$. This is a major legacy from the HFI for future CMB experiments. The removal of bandpass leakage has been improved by extracting the bandpass-mismatch coefficients for each detector as part of the mapmaking process; these values in turn improve the intensity maps. This is a major change in the philosophy of "frequency maps", which are now computed from single detector data, all adjusted to the same average bandpass response for the main foregrounds. Simulations reproduce very well the relative gain calibration of detectors, as well as drifts within a frequency induced by the residuals of the main systematic effect. Using these simulations, we measure and correct the small frequency calibration bias induced by this systematic effect at the $10^{-4}$ level. There is no detectable sign of a residual calibration bias between the first and second acoustic peaks in the CMB channels, at the $10^{-3}$ level.

We discuss constraints on cosmic reionisation and their implications on a cosmic SFR density $\rho_\mathrm{SFR}$ model; we study the influence of key-parameters such as the clumping factor of ionised hydrogen in the intergalactic medium (IGM) $C_{H_{II}}$ and the fraction of ionising photons escaping star-forming galaxies to reionise the IGM $f_\mathrm{esc}$. Our analysis uses SFR history data coming from luminosity functions, assuming that star-forming galaxies were sufficient to lead the reionisation process at high redshift. We add two other sets of constraints: measurements of the IGM ionised fraction and the most recent result from Planck Satellite about the integrated Thomson optical depth of the Cosmic Microwave Background (CMB) $\tau_\mathrm{Planck}$. We also consider various possibilities for the evolution of these two parameters with redshift, and confront them with observational data cited above. We conclude that, if the model of a constant clumping factor is chosen, the fiducial value of $3$ often used in papers is consistent with observations; even if a redshift-dependent model is considered, the resulting optical depth is strongly correlated to $C_{H_{II}}$ mean value at $z>7$, an additional argument in favour of the use of a constant clumping factor. Besides, the escape fraction is related to too many astrophysical parameters to allow us to use a complete and fully satisfactory model. A constant value with redshift seems again to be the most likely expression: considering it as a fit parameter, we get from the maximum likelihood (ML) model $f_\mathrm{esc}=0.24\pm0.08$; with a redshift-dependent model, we find an almost constant evolution, slightly increasing with $z$, around $f_\mathrm{esc}=0.23$. Last, our analysis shows that a reionisation beginning as early as $z\geq14$ and persisting until $z\sim6$ is a likely storyline.

The Planck catalogues of SZ sources, PSZ1 and PSZ2, are the largest catalogues of galaxy clusters selected through their SZ signature in the full sky. In 2013, we started a long-term observational program at Canary Island observatories with the aim of validating about 500 unconfirmed SZ sources. In this work we present results of the initial pre-screening of possible cluster counterparts using photometric and spectroscopic data of the Sloan Digital Sky Survey DR12. Our main aim is to identify previously unconfirmed PSZ2 cluster candidates and to contribute in determination of the actual purity and completeness of Planck SZ source sample. Using the latest version of the PSZ2 catalogue, we select all sources overlapping with the SDSS DR12 footprint and without redshift information. We validate these cluster fields following optical criteria (mainly distance with respect to the Planck pointing, magnitude of the brightest cluster galaxy and cluster richness) and combining them with the profiles of the Planck Compton y-maps. Together, this procedure allows for a more robust identification of optical counterparts compared to simply cross-matching with existing SDSS cluster catalogues that have been constructed from earlier SDSS Data Releases. The sample contains new redshifts for 37 Planck galaxy clusters that were not included in the original release of PSZ2 Planck catalogue. We detect three cases as possible multiple counterparts. We show that a combination of all available information (optical images and profile of SZ signal) can provide correct associations between the observed Planck SZ source and the optically identified cluster. We also show that Planck SZ detection is very sensitive even to high-z (z>0.5) clusters. In addition, we also present updated spectroscopic information for 34 Planck PSZ1 sources (33 previously photometrically confirmed and 1 new identification).

The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration.

We identify new clusters and characterize previously unknown Planck Sunyaev-Zeldovich (SZ) sources from the first Planck catalogue of SZ sources (PSZ1). The results presented here correspond to an optical follow-up observational programme developed during approximately one year (2014) at Roque de los Muchachos Observatory, using the 2.5m Isaac Newton telescope, the 3.5m Telescopio Nazionale Galileo, the 4.2m William Herschel telescope and the 10.4m Gran Telescopio Canarias. We characterize 115 new PSZ1 sources using deep optical imaging and spectroscopy. We adopt robust criteria in order to consolidate the SZ counterparts by analysing the optical richness, the 2D galaxy distribution, and velocity dispersions of clusters. Confirmed counterparts are considered to be validated if they are rich structures, well aligned with the Planck PSZ1 coordinate and show relatively high velocity dispersion. Following this classification, we confirm 53 clusters, which means that 46% of this PSZ1 subsample has been validated and characterized with this technique. Sixty-two SZ sources (54% of this PSZ1 subset) remain unconfirmed. In addition, we find that the fraction of unconfirmed clusters close to the galactic plane (at |b|<25deg) is greater than that at higher galactic latitudes (|b|>25deg), which indicates contamination produced by radio emission of galactic dust and gas clouds on these SZ detections. In fact, in the majority of the cases, we detect important galactic cirrus in the optical images, mainly in the SZ target located at low galactic latitudes, which supports this hypothesis.

Thermal Sunyaev-Zel'dovich effect is one of the recent probes of cosmology and large scale structures. We update constraints on cosmological parameters from galaxy clusters observed by the Planck satellite in a first attempt to combine cluster number counts and power spectrum of hot gas, using the new value of the optical depth, and sampling at the same time on cosmological and scaling-relation parameters. We find that in the $\Lambda$CDM model, the addition of tSZ power spectrum provides only small improvements with respect to number counts only, leading to the $68\%$ c.l. constraints $\Omega_m = 0.32 \pm 0.02$, $\sigma_8 = 0.77\pm0.03 $ and $\sigma_8 (\Omega_m/0.3)^{1/3}= 0.78\pm0.03$ and lowering the discrepancy with CMB primary anisotropies results (updated with the new value of $\tau$) to $\simeq 1.6\, \sigma$ on $\sigma_8$. We analyse extensions to standard model, considering the effect of massive neutrinos and varying the equation of state parameter for dark energy. In the first case, we find that the addition of tSZ power spectrum helps in strongly improving cosmological constraints with respect to number counts only results, leading to the $95\%$ upper limit $\sum m_{\nu}< 1.53 \, \text{eV}$. For the varying dark energy EoS scenario, we find again no important improvements when adding tSZ power spectrum, but still the combination of tSZ probes is able in providing constraints, producing $w = -1.0\pm 0.2$. In all cosmological scenari the mass bias to reconcile CMB and tSZ probes remains low: $(1-b)\lesssim 0.66$ as compared to estimates from weak lensing and Xray mass estimate comparisons or numerical simulations.

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.

Galaxy clusters play a crucial role in constraining cosmological parameters. Cluster detection and characterization have therefore become an important field of modern cosmology. In this work, we aim at constructing reasonably representative composite luminosity function of cluster members in the WISE survey that can be used in the future as an input in matching algorithms for cluster detection. We use a sample of massive clusters from the redMaPPer catalog to match the positions of the WISE sources. We build the composite luminosity function of the WISE members in different redshift bins. We find that galaxy cluster members have a characteristic composite luminosity function, with a clear change in the slope at a given magnitude apparent $M^{*}$, which becomes fainter with increasing redshift. The best-fit bright-end slope $\beta$ is compatible with a constant value of $3.54$ with no clear trend with redshift, while the faint-end slope $\alpha$ is remarkably different, ranging from $2$ to $3$. We present the first characterization of the composite luminosity function of the WISE counterparts of redMaPPer cluster members and in this way we provide an element for building cluster detection techniques based on matching algorithms.

We present the results of optical identifications and spectroscopic redshifts measurements for galaxy clusters from 2-nd Planck catalogue of Sunyaev-Zeldovich sources (PSZ2), located at high redshifts, $z\approx0.7-0.9$. We used the data of optical observations obtained with Russian-Turkish 1.5-m telescope (RTT150), Sayan observatory 1.6-m telescope, Calar Alto 3.5-m telescope and 6-m SAO RAS telescope (Bolshoi Teleskop Alt-azimutalnyi, BTA). Spectroscopic redshift measurements were obtained for seven galaxy clusters, including one cluster, PSZ2 G126.57+51.61, from the cosmological sample of PSZ2 catalogue. In central regions of two clusters, PSZ2 G069.39+68.05 and PSZ2 G087.39-34.58, the strong gravitationally lensed background galaxies are found, one of them at redshift $z=4.262$. The data presented below roughly double the number of known galaxy clusters in the second Planck catalogue of Sunyaev-Zeldovich sources at high redshifts, $z\approx0.8$.

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.

We have performed a multi-wavelength analysis of two galaxy cluster systems selected with the thermal Sunyaev-Zel'dovich (tSZ) effect and composed of cluster pairs and an inter-cluster filament. We have focused on one pair of particular interest: A399-A401 at redshift z~0.073 seperated by 3 Mpc. We have also performed the first analysis of one lower significance newly associated pair: A21-PSZ2 G114.09-34.34 at z~0.094, separated by 4.2 Mpc. We have characterised the intra-cluster gas using the tSZ signal from Planck and, when this was possible, the galaxy optical and infra-red (IR) properties based on two photometric redshift catalogs: 2MPZ and WISExSCOS. From the tSZ data, we measured the gas pressure in the clusters and in the inter-cluster filaments. In the case of A399-A401, the results are in perfect agreement with previous studies and, using the temperature measured from the X-rays, we further estimate the gas density in the filament and find n0=4.3+-0.7x10^-4 cm-3. The optical and IR colour-colour and colour-magnitude analyses of the galaxies selected in the cluster system, together with their Star Formation Rate, show no segregation between galaxy populations, in the clusters and in the filament of A399-A401. Galaxies are all passive, early type, and red and dead. The gas and galaxy properties of this system suggest that the whole system formed at the same time and corresponds to a pre-merger, with a cosmic filament gas heated by the collapse. For the other cluster system, the tSZ analysis was performed and the pressure in the clusters and in the inter-cluster filament was constrained. However, the limited or nonexistent optical and IR data prevent us from concluding on the presence of an actual cosmic filament or from proposing a scenario.

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.

We derive and compare the fractions of cool-core clusters in the {\em Planck} Early Sunyaev-Zel'dovich sample of 164 clusters with $z \leq 0.35$ and in a flux-limited X-ray sample of 100 clusters with $z \leq 0.30$, using {\em Chandra} observations. We use four metrics to identify cool-core clusters: 1) the concentration parameter: the ratio of the integrated emissivity profile within 0.15 $r_{500}$ to that within $r_{500}$, and 2) the ratio of the integrated emissivity profile within 40 kpc to that within 400 kpc, 3) the cuspiness of the gas density profile: the negative of the logarithmic derivative of the gas density with respect to the radius, measured at 0.04 $r_{500}$, and 4) the central gas density, measured at 0.01 $r_{500}$. We find that the sample of X-ray selected clusters, as characterized by each of these metrics, contains a significantly larger fraction of cool-core clusters compared to the sample of SZ selected clusters (44$\pm$7\% vs. 28$\pm$4\% using the concentration parameter in the 0.15--1.0 $r_{500}$ range, 61$\pm$8\% vs. 36$\pm$5\% using the concentration parameter in the 40--400 kpc range, 64$\pm$8\% vs. 38$\pm$5\% using the cuspiness, and 53$\pm$7\% vs. 39$\pm$5\% using the central gas density). Qualitatively, cool-core clusters are more X-ray luminous at fixed mass. Hence, our X-ray flux-limited sample, compared to the approximately mass-limited SZ sample, is over-represented with cool-core clusters. We describe a simple quantitative model that uses the excess luminosity of cool-core clusters compared to non-cool-core clusters at fixed mass to successfully predict the observed fraction of cool-core clusters in X-ray selected samples.

The six parameters of the standard $\Lambda$CDM model have best-fit values derived from the Planck temperature power spectrum that are shifted somewhat from the best-fit values derived from WMAP data. These shifts are driven by features in the Planck temperature power spectrum at angular scales that had never before been measured to cosmic-variance level precision. We investigate these shifts to determine whether they are within the range of expectation and to understand their origin in the data. Taking our parameter set to be the optical depth of the reionized intergalactic medium $\tau$, the baryon density $\omega_{\rm b}$, the matter density $\omega_{\rm m}$, the angular size of the sound horizon $\theta_*$, the spectral index of the primordial power spectrum, $n_{\rm s}$, and $A_{\rm s}e^{-2\tau}$ (where $A_{\rm s}$ is the amplitude of the primordial power spectrum), we examine the change in best-fit values between a WMAP-like large angular-scale data set (with multipole moment $\ell<800$ in the Planck temperature power spectrum) and an all angular-scale data set ($\ell<2500$ Planck temperature power spectrum), each with a prior on $\tau$ of $0.07\pm0.02$. We find that the shifts, in units of the 1$\sigma$ expected dispersion for each parameter, are $\{\Delta \tau, \Delta A_{\rm s} e^{-2\tau}, \Delta n_{\rm s}, \Delta \omega_{\rm m}, \Delta \omega_{\rm b}, \Delta \theta_*\} = \{-1.7, -2.2, 1.2, -2.0, 1.1, 0.9\}$, with a $\chi^2$ value of 8.0. We find that this $\chi^2$ value is exceeded in 15% of our simulated data sets, and that a parameter deviates by more than 2.2$\sigma$ in 9% of simulated data sets, meaning that the shifts are not unusually large. Comparing $\ell<800$ instead to $\ell>800$, or splitting at a different multipole, yields similar results. We examine the $\ell<800$ model residuals in the $\ell>800$ power spectrum data and find that the features there... [abridged]

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.

We present the results of spectroscopic redshift measurements for the galaxy clusters from the first all-sky Planck catalogue of the Sunyaev-Zeldovich sources, that have been mostly identified by means of the optical observations performed previously by our team (Planck Collaboration, 2015a). The data on 13 galaxy clusters at redshifts from z=~0.2 to z=~0.8, including the improved identification and redshift measurement for the cluster PSZ1 G141.73+14.22 at z=0.828, are provided. The measurements were done using the data from Russian-Turkish 1.5-m telescope (RTT-150), 2.2-m Calar Alto Observatory telescope, and 6-m SAO RAS telescope (Bolshoy Teleskop Azimutalnyi, BTA).

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.

We investigate constraints on cosmic reionization extracted from the Planck cosmic microwave background (CMB) data. We combine the Planck CMB anisotropy data in temperature with the low-multipole polarization data to fit LCDM models with various parameterizations of the reionization history. We obtain a Thomson optical depth tau=0.058 +/- 0.012 for the commonly adopted instantaneous reionization model. This confirms, with only data from CMB anisotropies, the low value suggested by combining Planck 2015 results with other data sets and also reduces the uncertainties. We reconstruct the history of the ionization fraction using either a symmetric or an asymmetric model for the transition between the neutral and ionized phases. To determine better constraints on the duration of the reionization process, we also make use of measurements of the amplitude of the kinetic Sunyaev-Zeldovich (kSZ) effect using additional information from the high resolution Atacama Cosmology Telescope and South Pole Telescope experiments. The average redshift at which reionization occurs is found to lie between z=7.8 and 8.8, depending on the model of reionization adopted. Using kSZ constraints and a redshift-symmetric reionization model, we find an upper limit to the width of the reionization period of Dz < 2.8. In all cases, we find that the Universe is ionized at less than the 10% level at redshifts above z~10. This suggests that an early onset of reionization is strongly disfavoured by the Planck data. We show that this result also reduces the tension between CMB-based analyses and constraints from other astrophysical sources.

Using the Planck 2015 data release (PR2) temperature maps, we separate Galactic thermal dust emission from cosmic infrared background (CIB) anisotropies. For this purpose, we implement a specifically tailored component-separation method, the so-called generalized needlet internal linear combination (GNILC) method, which uses spatial information (the angular power spectra) to disentangle the Galactic dust emission and CIB anisotropies. We produce significantly improved all-sky maps of Planck thermal dust emission, with reduced CIB contamination, at 353, 545, and 857 GHz. By reducing the CIB contamination of the thermal dust maps, we provide more accurate estimates of the local dust temperature and dust spectral index over the sky with reduced dispersion, especially at high Galactic latitudes above $b = \pm 20{\deg}$. We find that the dust temperature is $T = (19.4 \pm 1.3)$ K and the dust spectral index is $\beta = 1.6 \pm 0.1$ averaged over the whole sky, while $T = (19.4 \pm 1.5)$ K and $\beta = 1.6 \pm 0.2$ on 21 % of the sky at high latitudes. Moreover, subtracting the new CIB-removed thermal dust maps from the CMB-removed Planck maps gives access to the CIB anisotropies over 60 % of the sky at Galactic latitudes $|b| > 20{\deg}$. Because they are a significant improvement over previous Planck products, the GNILC maps are recommended for thermal dust science. The new CIB maps can be regarded as indirect tracers of the dark matter and they are recommended for exploring cross-correlations with lensing and large-scale structure optical surveys. The reconstructed GNILC thermal dust and CIB maps are delivered as Planck products.

Parity violating extensions of the standard electromagnetic theory cause in vacuo rotation of the plane of polarization of propagating photons. This effect, also known as cosmic birefringence, impacts the cosmic microwave background (CMB) anisotropy angular power spectra, producing non-vanishing $T$--$B$ and $E$--$B$ correlations that are otherwise null when parity is a symmetry. Here we present new constraints on an isotropic rotation, parametrized by the angle $\alpha$, derived from Planck 2015 CMB polarization data. To increase the robustness of our analyses, we employ two complementary approaches, in harmonic space and in map space, the latter based on a peak stacking technique. The two approaches provide estimates for $\alpha$ that are in agreement within statistical uncertainties and very stable against several consistency tests. Considering the $T$--$B$ and $E$--$B$ information jointly, we find $\alpha = 0.31^{\circ} \pm 0.05^{\circ} \, ({\rm stat.})\, \pm 0.28^{\circ} \, ({\rm syst.})$ from the harmonic analysis and $\alpha = 0.35^{\circ} \pm 0.05^{\circ} \, ({\rm stat.})\, \pm 0.28^{\circ} \, ({\rm syst.})$ from the stacking approach. These constraints are compatible with no parity violation and are dominated by the systematic uncertainty in the orientation of Planck's polarization-sensitive bolometers.

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.

The characterization of the Galactic foregrounds has been shown to be the main obstacle in the challenging quest to detect primordial B-modes in the polarized microwave sky. We make use of the Planck-HFI 2015 data release at high frequencies to place new constraints on the properties of the polarized thermal dust emission at high Galactic latitudes. Here, we specifically study the spatial variability of the dust polarized spectral energy distribution, and its potential impact on the determination of the tensor-to-scalar ratio. We use the correlation ratio of the $C_\ell^{BB}$ angular power spectra between the 217- and 353-GHz channels as a tracer of these potential variations, computed on different high Galactic latitude regions, ranging from 80% to 20% of the sky. The new insight from Planck data is a departure of the correlation ratio from unity that cannot be attributed to a spurious decorrelation due to the cosmic microwave background, instrumental noise, or instrumental systematics. The effect is marginally detected on each region, but the statistical combination of all the regions gives more than 99% confidence for this variation in polarized dust properties. In addition, we show that the decorrelation increases when there is a decrease in the mean column density of the region of the sky being considered, and we propose a simple power-law empirical model for this dependence, which matches what is seen in the Planck data. We explore the effect that this measured decorrelation has on simulations of the BICEP2-Keck Array/Planck analysis and show that the 2015 constraints from those data still allow a decorrelation between the dust at 150 and 353GHz of the order of the one we measure. Finally we show that either spatial variation of the dust SED or of the dust polarization angle could produce decorrelations between 217- and 353-GHz data similar to those we observe in the data.

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]

This paper describes the identification, modelling, and removal of previously unexplained systematic effects in the polarization data of the Planck High Frequency Instrument (HFI) on large angular scales, including new mapmaking and calibration procedures, new and more complete end-to-end simulations, and a set of robust internal consistency checks on the resulting maps. These maps, at 100, 143, 217, and 353 GHz, are early versions of those that will be released in final form later in 2016. The improvements allow us to determine the cosmic reionization optical depth $\tau$ using, for the first time, the low-multipole $EE$ data from HFI, reducing significantly the central value and uncertainty, and hence the upper limit. Two different likelihood procedures are used to constrain $\tau$ from two estimators of the CMB $E$- and $B$-mode angular power spectra at 100 and 143 GHz, after debiasing the spectra from a small remaining systematic contamination. These all give fully consistent results. A further consistency test is performed using cross-correlations derived from the Low Frequency Instrument maps of the Planck 2015 data release and the new HFI data. For this purpose, end-to-end analyses of systematic effects from the two instruments are used to demonstrate the near independence of their dominant systematic error residuals. The tightest result comes from the HFI-based $\tau$ posterior distribution using the maximum likelihood power spectrum estimator from $EE$ data only, giving a value $0.055\pm 0.009$. In a companion paper these results are discussed in the context of the best-fit Planck $\Lambda$CDM cosmological model and recent models of reionization.

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.

We study the statistical properties of interstellar dust polarization at high Galactic latitude, using the Stokes parameter Planck maps at 353 GHz. Our aim is to advance the understanding of the magnetized interstellar medium (ISM), and to provide a model of the polarized dust foreground for cosmic microwave background component-separation procedures. Focusing on the southern Galactic cap, we examine the statistical distributions of the polarization fraction ($p$) and angle ($\psi$) to characterize the ordered and turbulent components of the Galactic magnetic field (GMF) in the solar neighbourhood. We relate patterns at large angular scales in polarization to the orientation of the mean (ordered) GMF towards Galactic coordinates $(l_0,b_0)=(70^\circ \pm 5^\circ,24^\circ \pm 5^\circ)$. The histogram of $p$ shows a wide dispersion up to 25 %. The histogram of $\psi$ has a standard deviation of $12^\circ$ about the regular pattern expected from the ordered GMF. We use these histograms to build a phenomenological model of the turbulent component of the GMF, assuming a uniform effective polarization fraction ($p_0$) of dust emission. To model the Stokes parameters, we approximate the integration along the line of sight (LOS) as a sum over a set of $N$ independent polarization layers, in each of which the turbulent component of the GMF is obtained from Gaussian realizations of a power-law power spectrum. We are able to reproduce the observed $p$ and $\psi$ distributions using: a $p_0$ value of (26 $\pm$ 3)%; a ratio of 0.9 $\pm$ 0.1 between the strengths of the turbulent and mean components of the GMF; and a small value of $N$. We relate the polarization layers to the density structure and to the correlation length of the GMF along the LOS.

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 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.