We investigate the role of dense Mpc-scale environments in processing molecular gas of cluster galaxies as they fall into the cluster cores. We consider $\sim20$ luminous infrared galaxies (LIRGs) in intermediate-$z$ clusters, from the Hershel Lensing Survey and the Local Cluster Substructure Survey. They include MACS J0717.5+3745 at $z=0.546$ and Abell 697, 963, 1763, and 2219 at $z=0.2-0.3$. We have performed far infrared to ultraviolet spectral energy distribution modeling of the LIRGs, which span cluster-centric distances within $r/r_{200}\simeq0.2-1.6$. We have observed the LIRGs in CO(1$\rightarrow$0) or CO(2$\rightarrow$1) with the Plateau de Bure interferometer and its successor NOEMA, as part of five observational programs carried out between 2012 and 2017. We have compared the molecular gas to stellar mass ratio $M(H_2)/M_\star$, star formation rate (SFR), and depletion time ($\tau_{\rm dep}$) of the LIRGs with those of a compilation of cluster and field star forming galaxies. The targeted LIRGs have SFR, $M(H_2)/M_\star$, and $\tau_{\rm dep}$ that are consistent with those of both main sequence (MS) field galaxies and star forming galaxies from the comparison sample. However we find that the depletion time, normalized to the MS value, increases with increasing $r/r_{200}$, with a significance of $2.8\sigma$, which is ultimately due to a deficit of cluster core LIRGs with $\tau_{\rm dep}\gtrsim\tau_{\rm dep,MS}$. We suggest that a rapid exhaustion of the molecular gas reservoirs occurs in the cluster LIRGs and is effective in suppressing their star formation. This mechanism may explain the exponential decrease of the fraction of cluster LIRGs with cosmic time. The compression of the gas in LIRGs, possibly induced by intra-cluster medium shocks, may be responsible for the short depletion timescales, observed in a large fraction of cluster core LIRGs.

The main limiting factor of cosmological analyses based on thermal Sunyaev-Zel'dovich (SZ) cluster statistics comes from the bias and systematic uncertainties that affect the estimates of the mass of galaxy clusters. High-angular resolution SZ observations at high redshift are needed to study a potential redshift or morphology dependence of both the mean pressure profile and of the mass-observable scaling relation used in SZ cosmological analyses. The NIKA2 camera is a new generation continuum instrument installed at the IRAM 30-m telescope. With a large field of view, a high angular resolution and a high-sensitivity, the NIKA2 camera has unique SZ mapping capabilities. In this paper, we present the NIKA2 SZ large program, aiming at observing a large sample of clusters at redshifts between 0.5 and 0.9, and the characterization of the first cluster oberved with NIKA2.

Clusters of galaxies, the largest bound objects in the Universe, constitute a cosmological probe of choice, which is sensitive to both dark matter and dark energy. Within this framework, the Sunyaev-Zel'dovich (SZ) effect has opened a new window for the detection of clusters of galaxies and for the characterization of their physical properties such as mass, pressure and temperature. NIKA, a KID-based dual band camera installed at the IRAM 30-m telescope, was particularly well adapted in terms of frequency, angular resolution, field-of-view and sensitivity, for the mapping of the thermal and kinetic SZ effect in high-redshift clusters. In this paper, we present the NIKA cluster sample and a review of the main results obtained via the measurement of the SZ effect on those clusters: reconstruction of the cluster radial pressure profile, mass, temperature and velocity.

We present the results of the analysis of the very massive cluster MOO J1142$+$1527 at a redshift $z = 1.2$ based on high angular resolution NIKA2 Sunyaev-Zel'dovich (SZ) and $Chandra$ X-ray data. This multi-wavelength analysis enables us to estimate the shape of the temperature profile with unprecedented precision at this redshift and to obtain a map of the gas entropy distribution averaged along the line of sight. The comparison between the cluster morphological properties observed in the NIKA2 and $Chandra$ maps together with the analysis of the entropy map allows us to conclude that MOO J1142$+$1527 is an on-going merger hosting a cool-core at the position of the X-ray peak. This work demonstrates how the addition of spatially-resolved SZ observations to low signal-to-noise X-ray data can bring valuable insights on the intracluster medium thermodynamic properties at $z>1$.

High-resolution mapping of the hot gas in galaxy clusters is a key tool for cluster-based cosmological analyses. Taking advantage of the NIKA2 millimeter camera operated at the IRAM 30-m telescope, the NIKA2 SZ Large Program seeks to get a high-resolution follow-up of 45 galaxy clusters covering a wide mass range at high redshift in order to re-calibrate some of the tools needed for the cosmological exploitation of SZ surveys. We present the second cluster analysis of this program, targeting one of the faintest sources of the sample in order to tackle the difficulties in data reduction for such faint, low-SNR clusters. In this study, the main challenge is the precise estimation of the contamination by sub-millimetric point sources, which greatly affects the tSZ map of the cluster. We account for this contamination by performing a joint fit of the SZ signal and of the flux density of the compact sources. A prior knowledge of these fluxes is given by the adjustment of the SED of each source using data from both NIKA2 and the \textit{Herschel} satellite. The first results are very promising and demonstrate the possibility to estimate thermodynamic properties with NIKA2, even in a compact cluster heavily contaminated by point sources.

The characterization of the Intra-Cluster Medium (ICM) properties of high-redshift galaxy clusters is fundamental to our understanding of large-scale structure formation processes. We present the results of a multi-wavelength analysis of the very massive cluster MOO J1142$+$1527 at a redshift $z = 1.2$ discovered as part of the Massive and Distant Clusters of WISE Survey (MaDCoWS). This analysis is based on high angular resolution $Chandra$ X-ray and NIKA2 Sunyaev-Zel'dovich (SZ) data. Although the X-ray data have only about 1700 counts, we are able to determine the ICM thermodynamic radial profiles, namely temperature, entropy, and hydrostatic mass. These have been obtained with unprecedented precision at this redshift and up to $0.7R_{500}$, thanks to the combination of high-resolution X-ray and SZ data. The comparison between the galaxy distribution mapped in infrared by $Spitzer$ and the morphological properties of the ICM derived from the combined analysis of the $Chandra$ and NIKA2 data leads us to the conclusion that the cluster is an on-going merger. We measure the hydrostatic mass profile of the cluster in four angular sectors centered on the large-scale X-ray centroid. This allows us to estimate a systematic uncertainty on the cluster total mass that characterizes both the impact of the observed deviations from spherical symmetry and of the core dynamics on the mass profile. We further combine the X-ray and SZ data at the pixel level to obtain maps of the temperature and entropy distributions averaged along the line of sight. We find a relatively low entropy core at the position of the X-ray peak and high temperature regions located on its south and west sides. The increase in ICM temperature at the location of the SZ peak is expected given the merger dynamics. (abridged)

We study the dynamical state and the integrated total mass profiles of 75 massive (M500 > 5 e+14 M_sun) SZ-selected clusters at 0.08<z< 1.1. The sample is built from the Planck catalogue, with the addition of 4 SPT clusters at z>0.9. Using XMM imaging observations, we characterise the dynamical state with the centroid shift, the concentration, and their combination, M, which simultaneously probes the core and the large scale gas morphology. Using spatially-resolved spectroscopy and assuming hydrostatic equilibrium, we derive the total integrated mass profiles. The mass profile shape is quantified by the sparsity, the ratio of M500 to M2500, the masses at density contrast 500 and 2500, respectively. We study the correlations between the various parameters and their dependence on redshift. We confirm that SZ-selected samples, thought to reflect most closely the underlying cluster population, are dominated by disturbed and non-cool core objects at all z. There is no significant evolution or mass dependence of either the cool core fraction or the centroid shift parameter. The M parameter evolves slightly with z, having a correlation coefficient of rho= -0.2 $\pm$ 0.1 and a null hypothesis p-value of 0.01. In the high mass regime considered here, the sparsity evolves minimally with redshift, increasing by 10% between z<0.2 and z>0.55, an effect significant at less than 2 sigma. In contrast, the dependence of the sparsity on dynamical state is much stronger, increasing by a factor of $\sim$60% from the 1/3 most relaxed to the 1/3 most disturbed objects, an effect significant at more than 3 sigma. This is the first observational evidence that the shape of the integrated total mass profile in massive clusters is principally governed by the dynamical state, and is only mildly dependent on redshift. We discuss the consequences for the comparison between observations and theoretical predictions.

The Planck satellite has detected cluster candidates via the Sunyaev Zel'dovich (SZ) effect, but the optical follow-up required to confirm these candidates is still incomplete, especially at high redshifts and for SZ detections at low significance. In this work we present our analysis of optical observations obtained for 32 Planck cluster candidates using ACAM on the 4.2-m William Herschel Telescope. These cluster candidates were preselected using SDSS, WISE, and Pan-STARRS images to likely represent distant clusters at redshifts $z \gtrsim 0.7$. We obtain photometric redshift and richness estimates for all of the cluster candidates from a red-sequence analysis of $r$-, $i$-, and $z$-band imaging data. In addition, long-slit observations allow us to measure the redshifts of a subset of the clusters spectroscopically. The optical richness is often lower than expected from the inferred SZ mass when compared to scaling relations previously calibrated at low redshifts. This likely indicates the impact of Eddington bias and projection effects or noise-induced detections, especially at low SZ-significance. Thus, optical follow-up not only provides redshift measurements, but also an important independent verification method. We find that 18 (7) of the candidates at redshifts $z > 0.5$ ($z > 0.8$) are at least half as rich as expected from scaling relations, thereby clearly confirming these candidates as massive clusters. While the complex selection function of our sample due to our preselection hampers its use for cosmological studies, we do provide a validation of massive high-redshift clusters particularly suitable for further astrophysical investigations.

Galaxy cluster counts in bins of mass and redshift have been shown to be a competitive probe to test cosmological models. This method requires an efficient blind detection of clusters from surveys with a well-known selection function and robust mass estimates. The Euclid wide survey will cover 15000 deg$^2$ of the sky in the optical and near-infrared bands, down to magnitude 24 in the $H$-band. The resulting data will make it possible to detect a large number of galaxy clusters spanning a wide-range of masses up to redshift $\sim 2$. This paper presents the final results of the Euclid Cluster Finder Challenge (CFC). The objective of these challenges was to select the cluster detection algorithms that best meet the requirements of the Euclid mission. The final CFC included six independent detection algorithms, based on different techniques, such as photometric redshift tomography, optimal filtering, hierarchical approach, wavelet and friend-of-friends algorithms. These algorithms were blindly applied to a mock galaxy catalog with representative Euclid-like properties. The relative performance of the algorithms was assessed by matching the resulting detections to known clusters in the simulations. Several matching procedures were tested, thus making it possible to estimate the associated systematic effects on completeness to $<3$%. All the tested algorithms are very competitive in terms of performance, with three of them reaching $>80$% completeness for a mean purity of 80% down to masses of $10^{14}$ M$_{\odot}$ and up to redshift $z=2$. Based on these results, two algorithms were selected to be implemented in the Euclid pipeline, the AMICO code, based on matched filtering, and the PZWav code, based on an adaptive wavelet approach. [abridged]

We present the first all-sky catalogue of galaxy clusters and cluster candidates obtained from joint X-ray-SZ detections using observations from the Planck satellite and the ROSAT all-sky survey (RASS). The catalogue contains 2323 objects and has been validated by careful cross-identification with previously known clusters. This validation shows that 1597 candidates correspond to already known clusters, 212 coincide with other cluster candidates still to be confirmed, and the remaining 514 are completely new detections. With respect to Planck catalogues, the ComPRASS catalogue is simultaneously more pure and more complete. Based on the validation results in the SPT and SDSS footprints, the expected purity of the catalogue is at least 84.5%, meaning that more than 365 clusters are expected to be found among the new or still to be confirmed candidates with future validation efforts or specific follow-ups.

In recent years, the outskirts of galaxy clusters have emerged as one of the new frontiers and unique laboratories for studying the growth of large scale structure in the universe. Modern cosmological hydrodynamical simulations make firm and testable predictions of the thermodynamic and chemical evolution of the X-ray emitting intracluster medium. However, recent X-ray and Sunyaev-Zeldovich effect observations have revealed enigmatic disagreements with theoretical predictions, which have motivated deeper investigations of a plethora of astrophysical processes operating in the virialization region in the cluster outskirts. Much of the physics of cluster outskirts is fundamentally different from that of cluster cores, which has been the main focus of X-ray cluster science over the past several decades. A next-generation X-ray telescope, equipped with sub-arcsecond spatial resolution over a large field of view along with a low and stable instrumental background, is required in order to reveal the full story of the growth of galaxy clusters and the cosmic web and their applications for cosmology.

The Sunyaev-Zeldovich (SZ) effect was first predicted nearly five decades ago, but has only recently become a mature tool for performing high resolution studies of the warm and hot ionized gas in and between galaxies, groups, and clusters. Galaxy groups and clusters are powerful probes of cosmology, and they also serve as hosts for roughly half of the galaxies in the Universe. In this white paper, we outline the advances in our understanding of thermodynamic and kinematic properties of the warm-hot universe that can come in the next decade through spatially and spectrally resolved measurements of the SZ effects. Many of these advances will be enabled through new/upcoming millimeter/submillimeter (mm/submm) instrumentation on existing facilities, but truly transformative advances will require construction of new facilities with larger fields of view and broad spectral coverage of the mm/submm bands.

The total mass of a galaxy cluster is one of its most fundamental properties. Together with the redshift, the mass links observation and theory, allowing us to use the cluster population to test models of structure formation and to constrain cosmological parameters. Building on the rich heritage from X-ray surveys, new results from Sunyaev-Zeldovich and optical surveys have stimulated a resurgence of interest in cluster cosmology. These studies have generally found fewer clusters than predicted by the baseline Planck LCDM model, prompting a renewed effort on the part of the community to obtain a definitive measure of the true cluster mass scale. Here we review recent progress on this front. Our theoretical understanding continues to advance, with numerical simulations being the cornerstone of this effort. On the observational side, new, sophisticated techniques are being deployed in individual mass measurements and to account for selection biases in cluster surveys. We summarise the state of the art in cluster mass estimation methods and the systematic uncertainties and biases inherent in each approach, which are now well identified and understood, and explore how current uncertainties propagate into the cosmological parameter analysis. We discuss the prospects for improvements to the measurement of the mass scale using upcoming multi-wavelength data, and the future use of the cluster population as a cosmological probe.

NIKA2 is a dual-band millimetric camera of thousands of Kinetic Inductance Detectors (KID) installed at the IRAM 30-meter telescope in the Spanish Sierra Nevada. The instrument commissioning was completed in September 2017, and NIKA2 is now open to the scientific community and will operate for the next decade. NIKA2 has well-adapted instrumental design and performance to produce high-resolution maps of the thermal Sunyaev-Zel'dovich (SZ) effect toward intermediate and high redshift galaxy clusters. Moreover, it benefits from a guaranteed time large program dedicated to mapping a representative sample of galaxy clusters via SZ and that includes X-ray follow-ups. The main expected outputs of the SZ large program are the constraints on the redshift evolution of the pressure profile and the mass-observable relation. The first SZ mapping of a galaxy cluster with NIKA2 was produced, as part of the SZ large program. We found a sizable impact of the intracluster medium dynamics on the integrated SZ observables. This shows NIKA2 capabilities for the precise characterisation of the mass-observable relation that is required for accurate cosmology with galaxy clusters.

We study the abundance of star-forming and quiescent galaxies in a sample of 21 massive clusters at 0.5<z<0.7, detected with the Planck satellite. We measure the cluster galaxy stellar mass function (SMF), which is a fundamental observable to study and constrain the formation and evolution of galaxies. Our measurements are based on homogeneous and deep multi-band photometry spanning u- to the Ks-band for each cluster and are supported by spectroscopic data from different programs. The galaxy population is separated between quiescent and star-forming galaxies based on their rest-frame U-V and V-J colours. The SMF is compared to that of field galaxies at the same redshifts, using data from the COSMOS/UltraVISTA survey. We find that the shape of the SMF of star-forming galaxies does not depend on environment, while the SMF of quiescent galaxies has a significantly steeper low-mass slope in the clusters compared to the field. We estimate the environmental quenching efficiency (f_EQ), i.e. the probability for a galaxy that would normally be star forming in the field, to be quenched due to its environment. The f_EQ shows no stellar-mass dependence in any environment, but it increases from 40% in the cluster outskirts to ~90% in the cluster centres. The radial signature of f_EQ provides constraints on where the dominant quenching mechanism operates in these clusters and on what timescale. Exploring these using a simple model based on galaxy orbits obtained from an N-body simulation, we find a clear degeneracy between both parameters. For example, the quenching process may either be triggered on a long (~3 Gyr) time scale at large radii (r~8R_500), or happen well within 1 Gyr at r<R_500. The radius where quenching is triggered is at least r_quench> 0.67R_500 (95%CL). The ICM density at this location suggests that ram-pressure stripping of the cold gas is a likely cause of quenching. [Abridged]

We present a detailed study of the integrated total hydrostatic mass profiles of the five most massive ($M^{\mathrm{SZ}}_{500} > 5 \times 10^{14}$ M$_{\odot}$) galaxy clusters selected at $z\sim1$ via the Sunyaev-Zel'dovich effect. These objects represent an ideal laboratory to test structure formation models where the primary driver is gravity. Optimally exploiting spatially-resolved spectroscopic information from XMM-Newton and Chandra observations, we used both parametric (forward, backward) and non-parametric methods to recover the mass profiles, finding that the results are extremely robust when density and temperature measurements are both available. Our X-ray masses at $R_{500}$ are higher than the weak lensing masses obtained from the Hubble Space Telescope (HST), with a mean ratio of $1.39^{+0.47}_{-0.35}$. This offset goes in the opposite direction to that expected in a scenario where the hydrostatic method yields a biased, underestimated, mass. We investigated halo shape parameters such as sparsity and concentration, and compared to local X-ray selected clusters, finding hints for evolution in the central regions (or for selection effects). The total baryonic content is in agreement with the cosmic value at $R_{500}$. Comparison with numerical simulations shows that the mass distribution and concentration are in line with expectations. These results illustrate the power of X-ray observations to probe the statistical properties of the gas and total mass profiles in this high-mass, high-redshift regime.

The complete characterization of the pressure profile of high-redshift galaxy clusters, from their core to their outskirts, is a major issue for the study of the formation of large-scale structures. It is essential to constrain a potential redshift evolution of both the slope and scatter of the mass-observable scaling relations used in cosmology studies based on cluster statistics. In this paper, we present the first thermal Sunyaev-Zel'dovich (tSZ) mapping of a cluster from the sample of the NIKA2 SZ large program that aims at constraining the redshift evolution of cluster pressure profiles and the tSZ-mass scaling relation. We have observed the galaxy cluster PSZ2 G144.83+25.11 at redshift $z=0.58$ with the NIKA2 camera, a dual-band (150 and 260 GHz) instrument operated at the IRAM 30-meter telescope. We identify a thermal pressure excess in the south-west region of PSZ2 G144.83+25.11 and a high redshift sub-millimeter point source that affect the intracluster medium (ICM) morphology of the cluster. The NIKA2 data are used jointly with tSZ data acquired by the MUSTANG, Bolocam and $Planck$ experiments in order to non-parametrically set the best constraints on the electronic pressure distribution from the cluster core ($\rm{R} \sim 0.02 \rm{R_{500}}$) to its outskirts ($\rm{R} \sim 3 \rm{R_{500}} $). We investigate the impact of the over-pressure region on the shape of the pressure profile and on the constraints on the integrated Compton parameter $\rm{Y_{500}}$. A hydrostatic mass analysis is also performed by combining the tSZ-constrained pressure profile with the deprojected electronic density profile from XMM-$Newton$. This allows us to conclude that the estimates of $\rm{Y_{500}}$ and $\rm{M_{500}}$ obtained from the analysis with and without masking the disturbed ICM region differ by 65 and 79% respectively. (abridged)

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.

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.

The combination of X-ray and SZ observations can potentially improve the cluster detection efficiency when compared to using only one of these probes, since both probe the same medium: the hot ionized gas of the intra-cluster medium. We present a method based on matched multifrequency filters (MMF) for detecting galaxy clusters from SZ and X-ray surveys. This method builds on a previously proposed joint X-ray-SZ extraction method (Tarr\'io et al. 2016) and allows to blindly detect clusters, that is finding new clusters without knowing their position, size or redshift, by searching on SZ and X-ray maps simultaneously. The proposed method is tested using data from the ROSAT all-sky survey and from the Planck survey. The evaluation is done by comparison with existing cluster catalogues in the area of the sky covered by the deep SPT survey. Thanks to the addition of the X-ray information, the joint detection method is able to achieve simultaneously better purity, better detection efficiency and better position accuracy than its predecessor Planck MMF, which is based on SZ maps only. For a purity of 85%, the X-ray-SZ method detects 141 confirmed clusters in the SPT region, whereas to detect the same number of confirmed clusters with Planck MMF, we would need to decrease its purity to 70%. We provide a catalogue of 225 sources selected by the proposed method in the SPT footprint, with masses ranging between 0.7 and 14.5 $\cdot 10^{14}$ Msun and redshifts between 0.01 and 1.2.

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.

Sub-structures in the hot gas of galaxy clusters are related to their formation history and to the astrophysical processes at play in the intracluster medium (ICM). The thermal Sunyaev-Zel'dovich (tSZ) effect is directly sensitive to the line-of-sight integrated ICM pressure, and is thus particularly adapted to study ICM sub-structures. We apply structure-enhancement filtering algorithms to high resolution tSZ observations of distant clusters, in order to search for pressure discontinuities, compressions, as well as secondary peaks in the ICM. The same filters are applied to synthetic tSZ images extracted from RHAPSODY-G hydrodynamic simulations, in order to better interpret the extracted features. We also study the noise propagation trough the filters and quantify the impact of systematic effects, point source residuals being identified as the dominant potential contaminant. In 3 of our 6 NIKA clusters we identify features at high S/N that show clear evidence for merger events. In MACSJ0717 (z=0.55), three strong pressure gradients are observed on the E, SE and W sectors, and two main peaks in the pressure distribution are identified. We observe a lack of tSZ compact structure in the cool-core cluster PSZ1G045.85 (z=0.61), and a tSZ gradient ridge dominates in the SE. In the highest z cluster, CLJ1227 (z=0.89), we detect a ~45" (360 kpc) long ridge pressure gradient associated with a secondary pressure peak in the W region. Our results show that current tSZ facilities have now reached the angular resolution and sensitivity to allow an exploration of the details of pressure sub-structures in clusters, even at high z. This opens the possibility to quantify the impact of the dynamical state on the relation between the tSZ signal and the mass of clusters, which is important when using tSZ clusters to test cosmological models. This work also marks the first NIKA cluster sample data release.

Galaxy clusters constitute a major cosmological probe. However, Planck 2015 results have shown a weak tension between CMB-derived and cluster-derived cosmological parameters. This tension might be due to poor knowledge of the cluster mass and observable relationship. As for now, arcmin resolution Sunyaev-Zeldovich (SZ) observations ({\it e.g.} SPT, ACT and Planck) only allowed detailed studies of the intra cluster medium for low redshift clusters ($z<0.2$). For high redshift clusters ($z>0.5$) high resolution and high sensitivity SZ observations are needed. With both a wide field of view (6.5 arcmin) and a high angular resolution (17.7 and 11.2 arcsec at 150 and 260 GHz), the NIKA2 camera installed at the IRAM 30-m telescope (Pico Veleta, Spain) is particularly well adapted for these observations. The NIKA2 SZ observation program will map a large sample of clusters (50) at redshifts between 0.5 and 0.9. As a pilot study for NIKA2, several clusters of galaxies have been observed with the pathfinder, NIKA, at the IRAM 30-m telescope to cover the various configurations and observation conditions expected for NIKA2.}

We investigate the evolution of the dark matter density profiles of the most massive galaxy clusters in the Universe. Using a `zoom-in' procedure on a large suite of cosmological simulations of total comoving volume of $3\,(h^{-1}\,\rm Gpc)^3$, we study the 25 most massive clusters in four redshift slices from $z\sim 1$ to the present. The minimum mass is $M_{500} > 5.5 \times 10^{14}$ M$_{\odot}$ at $z=1$. Each system has more than two million particles within $r_{500}$. Once scaled to the critical density at each redshift, the dark matter profiles within $r_{500}$ are strikingly similar from $z\sim1$ to the present day, exhibiting a low dispersion of 0.15 dex, and showing little evolution with redshift in the radial logarithmic slope and scatter. They have the running power law shape typical of the NFW-type profiles, and their inner structure, resolved to $3.8\,h^{-1}$ comoving kpc at $z=1$, shows no signs of converging to an asymptotic slope. Our results suggest that this type of profile is already in place at $z>1$ in the highest-mass haloes in the Universe, and that it remains exceptionally robust to merging activity.

We examine the reconstruction of galaxy cluster radial density profiles obtained from Chandra and XMM X-ray observations, using high quality data for a sample of twelve objects covering a range of morphologies and redshifts. By comparing the results obtained from the two observatories and by varying key aspects of the analysis procedure, we examine the impact of instrumental effects and of differences in the methodology used in the recovery of the density profiles. We find that the final density profile shape is particularly robust. We adapt the photon weighting vignetting correction method developed for XMM for use with Chandra data, and confirm that the resulting Chandra profiles are consistent with those corrected a posteriori for vignetting effects. Profiles obtained from direct deprojection and those derived using parametric models are consistent at the 1% level. At radii larger than $\sim$6", the agreement between Chandra and XMM is better than 1%, confirming an excellent understanding of the XMM PSF. We find no significant energy dependence. The impact of the well-known offset between Chandra and XMM gas temperature determinations on the density profiles is found to be negligible. However, we find an overall normalisation offset in density profiles of the order of $\sim$2.5%, which is linked to absolute flux cross-calibration issues. As a final result, the weighted ratios of Chandra to XMM gas masses computed at R2500 and R500 are r=1.03$\pm$0.01 and r=1.03$\pm$0.03, respectively. Our study confirms that the radial density profiles are robustly recovered, and that any differences between Chandra and XMM can be constrained to the $\sim$ 2.5% level, regardless of the exact data analysis details. These encouraging results open the way for the true combination of X-ray observations of galaxy clusters, fully leveraging the high resolution of Chandra and the high throughput of XMM.

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

The development of precision cosmology with clusters of galaxies requires high-angular resolution Sunyaev-Zel'dovich (SZ) observations. As for now, arcmin resolution SZ observations (e.g. SPT, ACT and Planck) only allowed detailed studies of the intra cluster medium for low redshift clusters (z<0.2). With both a wide field of view (6.5 arcmin) and a high angular resolution (17.7 and 11.2 arcsec at 150 and 260 GHz), the NIKA2 camera installed at the IRAM 30-m telescope (Pico Veleta, Spain), will bring valuable information in the field of SZ imaging of clusters of galaxies. The NIKA2 SZ observation program will allow us to observe a large sample of clusters (50) at redshifts between 0.4 and 0.9. As a pilot study for NIKA2, several clusters of galaxies have been observed with the pathfinder, NIKA, at the IRAM 30-m telescope to cover the various configurations and observation conditions expected for NIKA2.

X-ray observations show that galaxy clusters have a very large range of morphologies. The most disturbed systems which are good to study how clusters form and grow and to test physical models, may potentially complicate cosmological studies because the cluster mass determination becomes more challenging. Thus, we need to understand the cluster properties of our samples to reduce possible biases. This is complicated by the fact that different experiments may detect different cluster populations. For example, SZ selected cluster samples have been found to include a greater fraction of disturbed systems than X-ray selected samples. In this paper we determined eight morphological parameters for the Planck Early Sunyaev-Zeldovich (ESZ) objects observed with XMM-Newton. We found that two parameters, concentration and centroid-shift, are the best to distinguish between relaxed and disturbed systems. For each parameter we provide the values that allow one to select the most relaxed or most disturbed objects from a sample. We found that there is no mass dependence on the cluster dynamical state. By comparing our results with what was obtained with REXCESS clusters, we also confirm that indeed the ESZ clusters tend to be more disturbed, as found by previous studies.

We propose a method to map the temperature distribution of the hot gas in galaxy clusters that uses resolved images of the thermal Sunyaev-Zel'dovich (tSZ) effect in combination with X-ray data. Application to images from the New IRAM KIDs Array (NIKA) and XMM-Newton allows us to measure and determine the spatial distribution of the gas temperature in the merging cluster MACS J0717.5+3745, at $z=0.55$. Despite the complexity of the target object, we find a good morphological agreement between the temperature maps derived from X-ray spectroscopy only -- using XMM-Newton ($T_{\rm XMM}$) and Chandra ($T_{\rm CXO}$) -- and the new gas-mass-weighted tSZ+X-ray imaging method ($T_{\rm SZ}$). We correlate the temperatures from tSZ+X-ray imaging and those from X-ray spectroscopy alone and find that $T_{\rm SZ}$ is higher than $T_{\rm XMM}$ and lower than $T_{\rm CXO}$ by $\sim 10\%$ in both cases. Our results are limited by uncertainties in the geometry of the cluster gas, contamination from kinetic SZ ($\sim 10\%$), and the absolute calibration of the tSZ map ($7\%$). Investigation using a larger sample of clusters would help minimise these effects.

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.

Measurement of the gas velocity distribution in galaxy clusters provides insight into the physics of mergers, through which large scale structures form in the Universe. Velocity estimates within the intracluster medium (ICM) can be obtained via the Sunyaev-Zel'dovich (SZ) effect, but its observation is challenging both in term of sensitivity requirement and control of systematic effects, including the removal of contaminants. In this paper we report resolved observations, at 150 and 260 GHz, of the SZ effect toward the triple merger MACS J0717.5+3745 (z=0.55), using data obtained with the NIKA camera at the IRAM 30m telescope. Assuming that the SZ signal is the sum of a thermal (tSZ) and a kinetic (kSZ) component and by combining the two NIKA bands, we extract for the first time a resolved map of the kSZ signal in a cluster. The kSZ signal is dominated by a dipolar structure that peaks at -5.1 and +3.4 sigma, corresponding to two subclusters moving respectively away and toward us and coincident with the cold dense X-ray core and a hot region undergoing a major merging event. We model the gas electron density and line-of-sight velocity of MACS J0717.5+3745 as four subclusters. Combining NIKA data with X-ray observations from XMM-Newton and Chandra, we fit this model to constrain the gas line-of-sight velocity of each component, and we also derive, for the first time, a velocity map from kSZ data (i.e. that is model-dependent). Our results are consistent with previous constraints on the merger velocities, and thanks to the high angular resolution of our data, we are able to resolve the structure of the gas velocity. Finally, we investigate possible contamination and systematic effects with a special care given to radio and submillimeter galaxies. Among the sources that we detect with NIKA, we find one which is likely to be a high redshift lensed submillimeter galaxy.

We present a pilot X-ray study of the five most massive ($M_{500}>5 \times 10^{14} M_{\odot}$), distant (z~1), galaxy clusters detected via the Sunyaev-Zeldovich effect. We optimally combine XMM-Newton and Chandra X-ray observations by leveraging the throughput of XMM to obtain spatially-resolved spectroscopy, and the spatial resolution of Chandra to probe the bright inner parts and to detect embedded point sources. Capitalising on the excellent agreement in flux-related measurements, we present a new method to derive the density profiles, constrained in the centre by Chandra and in the outskirts by XMM. We show that the Chandra-XMM combination is fundamental for morphological analysis at these redshifts, the Chandra resolution being required to remove point source contamination, and the XMM sensitivity allowing higher significance detection of faint substructures. The sample is dominated by dynamically disturbed objects. We use the combined Chandra-XMM density profiles and spatially-resolved temperature profiles to investigate thermodynamic quantities including entropy and pressure. From comparison of the scaled profiles with the local REXCESS sample, we find no significant departure from standard self-similar evolution, within the dispersion, at any radius, except for the entropy beyond 0.7$R_{500}$. The baryon mass fraction tends towards the cosmic value, with a weaker dependence on mass than observed in the local Universe. We compare with predictions from numerical simulations. The present pilot study demonstrates the utility and feasibility of spatially-resolved analysis of individual objects at high-redshift through the combination of XMM and Chandra observations. Observations of a larger sample will allow a fuller statistical analysis to be undertaken, in particular of the intrinsic scatter in the structural and scaling properties of the cluster population. (abridged)

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 determination of the thermodynamic properties of clusters of galaxies at intermediate and high redshift can bring new insights into the formation of large-scale structures. It is essential for a robust calibration of the mass-observable scaling relations and their scatter, which are key ingredients for precise cosmology using cluster statistics. Here we illustrate an application of high resolution $(< 20$ arcsec) thermal Sunyaev-Zel'dovich (tSZ) observations by probing the intracluster medium (ICM) of the \planck-discovered galaxy cluster \psz\ at redshift $z = 0.61$, using tSZ data obtained with the NIKA camera, which is a dual-band (150 and 260~GHz) instrument operated at the IRAM 30-meter telescope. We deproject jointly NIKA and \planck\ data to extract the electronic pressure distribution from the cluster core ($R \sim 0.02\, R_{500}$) to its outskirts ($R \sim 3\, R_{500}$) non-parametrically for the first time at intermediate redshift. The constraints on the resulting pressure profile allow us to reduce the relative uncertainty on the integrated Compton parameter by a factor of two compared to the \planck\ value. Combining the tSZ data and the deprojected electronic density profile from \xmm\ allows us to undertake a hydrostatic mass analysis, for which we study the impact of a spherical model assumption on the total mass estimate. We also investigate the radial temperature and entropy distributions. These data indicate that \psz\ is a massive ($M_{500} \sim 5.5 \times 10^{14}$ M$_{\odot}$) cool-core cluster. This work is part of a pilot study aiming at optimizing the treatment of the NIKA2 tSZ large program dedicated to the follow-up of SZ-discovered clusters at intermediate and high redshifts. (abridged)

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.

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.

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]

NIKA2 (New IRAM KID Arrays) is a dual band (150 and 260 GHz) imaging camera based on Kinetic Inductance Detectors (KIDs) and designed to work at the IRAM 30 m telescope (Pico Veleta, Spain). Built on the experience of the NIKA prototype, NIKA2 has been installed at the 30 m focal plane in October 2015 and the commissioning phase is now ongoing. Through the thermal Sunyaev-Zeldovich (tSZ) effect, NIKA2 will image the ionized gas residing in clusters of galaxies with a resolution of 12 and 18 arcsec FWHM (at 150 and 260 GHz, respectively). We report on the recent tSZ measurements with the NIKA camera and discuss the future objectives for the NIKA2 SZ large Program, 300h of observation dedicated to SZ science. With this program we intend to perform a high angular resolution follow-up of a cosmologically-representative sample of clusters belonging to SZ catalogues, with redshift greater than 0.5. The main output of the program will be the study of the redshift evolution of the cluster pressure profile as well as that of the scaling laws relating the cluster global properties.

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 investigate the properties of the hot gas in four fossil galaxy systems detected at high significance in the Planck Sunyaev-Zeldovich (SZ) survey. XMM-Newton observations reveal overall temperatures of kT ~ 5-6 keV and yield hydrostatic masses M500,HE > 3.5 x 10e14 Msun, confirming their nature as bona fide massive clusters. We measure the thermodynamic properties of the hot gas in X-rays (out to beyond R500 in three cases) and derive their individual pressure profiles out to R ~ 2.5 R500 with the SZ data. We combine the X-ray and SZ data to measure hydrostatic mass profiles and to examine the hot gas content and its radial distribution. The average Navarro-Frenk-White (NFW) concentration parameter, c500 = 3.2 +/- 0.4, is the same as that of relaxed `normal' clusters. The gas mass fraction profiles exhibit striking variation in the inner regions, but converge to approximately the cosmic baryon fraction (corrected for depletion) at R500. Beyond R500 the gas mass fraction profiles again diverge, which we interpret as being due to a difference in gas clumping and/or a breakdown of hydrostatic equilibrium in the external regions. Overall our observations point to considerable radial variation in the hot gas content and in the gas clumping and/or hydrostatic equilibrium properties in these fossil clusters, at odds with the interpretation of their being old, evolved and undisturbed. At least some fossil objects appear to be dynamically young.

The hot ionized gas of the intra-cluster medium emits thermal radiation in the X-ray band and also distorts the cosmic microwave radiation through the Sunyaev-Zel'dovich (SZ) effect. Combining these two complementary sources of information through innovative techniques can therefore potentially improve the cluster detection rate when compared to using only one of the probes. Our aim is to build such a joint X-ray-SZ analysis tool, which will allow us to detect fainter or more distant clusters while maintaining high catalogue purity. We present a method based on matched multifrequency filters (MMF) for extracting cluster catalogues from SZ and X-ray surveys. We first designed an X-ray matched-filter method, analogous to the classical MMF developed for SZ observations. Then, we built our joint X-ray-SZ algorithm by combining our X-ray matched filter with the classical SZ-MMF, for which we used the physical relation between SZ and X-ray observations. We show that the proposed X-ray matched filter provides correct photometry results, and that the joint matched filter also provides correct photometry when the $F_{\rm X}/Y_{500}$ relation of the clusters is known. Moreover, the proposed joint algorithm provides a better signal-to-noise ratio than single-map extractions, which improves the detection rate even if we do not exactly know the $F_{\rm X}/Y_{500}$ relation. The proposed methods were tested using data from the ROSAT all-sky survey and from the Planck survey.

The thermal Sunyaev-Zeldovich effect (tSZ) is a powerful probe to study clusters of galaxies and is complementary with respect to X-ray, lensing or optical observations. Previous arcmin resolution tSZ observations ({\it e.g.} SPT, ACT and Planck) only enabled detailed studies of the intra-cluster medium morphology for low redshift clusters ($z < 0.2$). Thus, the development of precision cosmology with clusters requires high angular resolution observations to extend the understanding of galaxy cluster towards high redshift. NIKA2 is a wide-field (6.5 arcmin field of view) dual-band camera, operated at $100 \ {\rm mK}$ and containing $\sim 3300$ KID (Kinetic Inductance Detectors), designed to observe the millimeter sky at 150 and 260 GHz, with an angular resolution of 18 and 12 arcsec respectively. The NIKA2 camera has been installed on the IRAM 30-m telescope (Pico Veleta, Spain) in September 2015. The NIKA2 tSZ observation program will allow us to observe a large sample of clusters (50) at redshift ranging between 0.5 and 1. As a pathfinder for NIKA2, several clusters of galaxies have been observed at the IRAM 30-m telescope with the NIKA prototype to cover the various configurations and observation conditions expected for NIKA2.

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

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

The prototype of the NIKA2 camera, NIKA, is an instrument operating at the IRAM 30-m telescope, which can observe simultaneously at 150 and 260GHz. One of the main goals of NIKA2 is to measure the pressure distribution in galaxy clusters at high resolution using the thermal SZ (tSZ) effect. Such observations have already proved to be an excellent probe of cluster pressure distributions even at high redshifts. However, an important fraction of clusters host submm and/or radio point sources, which can significantly affect the reconstructed signal. Here we report on <20" resolution observations at 150 and 260GHz of the cluster MACSJ1424, which hosts both radio and submm point sources. We examine the morphology of the tSZ signal and compare it to other datasets. The NIKA data are combined with Herschel satellite data to study the SED of the submm point source contaminants. We then perform a joint reconstruction of the intracluster medium (ICM) electronic pressure and density by combining NIKA, Planck, XMM-Newton, and Chandra data, focusing on the impact of the radio and submm sources on the reconstructed pressure profile. We find that large-scale pressure distribution is unaffected by the point sources because of the resolved nature of the NIKA observations. The reconstructed pressure in the inner region is slightly higher when the contribution of point sources are removed. We show that it is not possible to set strong constraints on the central pressure distribution without accurately removing these contaminants. The comparison with X-ray only data shows good agreement for the pressure, temperature, and entropy profiles, which all indicate that MACSJ1424 is a dynamically relaxed cool core system. The present observations illustrate the possibility of measuring these quantities with a relatively small integration time, even at high redshift and without X-ray spectroscopy.

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

The all-sky Planck survey in 9 frequency bands was used to search for emission from all 274 known Galactic supernova remnants. Of these, 16 were detected in at least two Planck frequencies. The radio-through-microwave spectral energy distributions were compiled to determine the emission mechanism for microwave emission. In only one case, IC 443, is there high-frequency emission clearly from dust associated with the supernova remnant.In all cases, the low-frequency emission is from synchrotron radiation. A single power law, as predicted for a population of relativistic particles with energy distribution that extends continuously to high energies, is evident for many sources, including the Crab and PKS 1209-51/52. A decrease in flux density relative to the extrapolation of radio emission is evident in several sources. Their spectral energy distributions can be approximated as broken power laws, $S_\nu\propto\nu^{-\alpha}$, with the spectral index, alpha, increasing by 0.5-1 above a break frequency in the range 10-60 GHz. The break could be due to synchrotron losses.

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.

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

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

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

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

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

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

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

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

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

The Planck design and scanning strategy provide many levels of redundancy that can be exploited to provide tests of internal consistency. One of the most important is the comparison of the 70GHz and 100GHz channels. Based on different instrument technologies, with feeds located differently in the focal plane, analysed independently by different teams using different software, and near the minimum of diffuse foreground emission, these channels are in effect two different experiments. The 143GHz channel has the lowest noise level on Planck, and is near the minimum of unresolved foreground emission. In this paper, we analyse the level of consistency achieved in the 2013 Planck data. We concentrate on comparisons between the 70/100/143GHz channel maps and power spectra, particularly over the angular scales of the first and second acoustic peaks, on maps masked for diffuse Galactic emission and for strong unresolved sources. Difference maps covering angular scales from 8deg-15arcmin are consistent with noise, and show no evidence of cosmic microwave background structure. Including small but important corrections for unresolved-source residuals, we demonstrate agreement between 70 and 100GHz power spectra averaged over 70<l<390 at the 0.8% level, and agreement between 143 and 100GHz power spectra of 0.4% over the same l range. These values are within and consistent with the overall uncertainties in calibration given in the Planck 2013 results. We also present results based on the 2013 likelihood analysis showing consistency at the 0.35% between the 100/143/217GHz power spectra. We analyse calibration procedures and beams to determine what fraction of these differences can be accounted for by known approximations or systematic errors that could be controlled even better in the future, reducing uncertainties still further. Several possible small improvements are described...(abridged)