We present the first metal abundance profiles for a representative sample of massive clusters. Our measures extend to $R_{500}$ and are corrected for a systematic error plaguing previous outskirt estimates. Our profiles flatten out at large radii, admittedly not a new result, however the radial range and representative nature of our sample extends its import well beyond previous findings. We find no evidence of segregation between cool-core and non-cool-core systems beyond $\sim 0.3 R_{500}$, implying that, as was found for thermodynamic properties (Ghirardini et al, 2019), the physical state of the core does not affect global cluster properties. Our mean abundance within $R_{500}$ shows a very modest scatter, $< $15%, suggesting the enrichment process must be quite similar in all these massive systems. This is a new finding and has significant implications on feedback processes. Together with results from thermodynamic properties presented in a previous X-COP paper, it affords a coherent picture where feedback effects do not vary significantly from one system to another. By combing ICM with stellar measurements we have found the amount of Fe diffused in the ICM to be about ten times higher than that locked in stars. Although our estimates suggest, with some strength, that the measured iron mass in clusters is well in excess of the predicted one, systematic errors prevent us from making a definitive statement. Further advancements will only be possible when systematic uncertainties, principally those associated to stellar masses, both within and beyond $R_{500}$, can be reduced.

We present the results of deep Chandra and XMM-Newton observations of a complex merging galaxy cluster Abell 2256 (A2256) that hosts a spectacular radio relic. The temperature and metallicity maps show clear evidence of a merger between the western subcluster and the primary cluster. We detect five X-ray surface brightness edges. Three of them near the cluster center are cold fronts (CFs): CF1 is associated with the infalling subcluster; CF2 is located in the east of the primary cluster and CF3 is to the west of the primary cluster core. The other two edges at cluster outskirts are shock fronts (SFs): SF1 near the radio relic in the NW has Mach numbers derived from the temperature and the density jumps respectively of $M_T=1.62\pm0.12$ and $M_\rho=1.23\pm0.06$; SF2 in the SE has $M_T=1.54\pm0.05$ and $M_\rho=1.16\pm0.13$. In the region of the radio relic there is no evidence for the correlation between X-ray and radio substructures, from which we estimate an upper limit for the inverse Compton emission, and therefore set a lower limit on the magnetic field ($\sim$ 450 kpc from primary cluster center) of $B>1.0\ \mu$G for a single power-law electron spectrum or $B>0.4\ \mu$G for a broken power-law electron spectrum. We propose a merger scenario including a primary cluster, a subcluster, and a group. Our merger scenario accounts for the X-ray edges, diffuse radio features, and galaxy kinematics, as well as projection effects.

The origin of radio relics is usually explained via diffusive shock acceleration (DSA) or re-acceleration of electrons at/from merger shocks in galaxy clusters. The case of acceleration is challenged by the low predicted efficiency of low-Mach number merger shocks, unable to explain the power observed in most radio relics. In this Letter we present the discovery of a new giant radio relic around the galaxy cluster Abell 2249 ($z=0.0838$) using LOFAR. It is special since it has the lowest surface brightness of all known radio relics. We study its radio and X-ray properties combinig LOFAR data with uGMRT, JVLA and XMM. This object has a total power of $L_{1.4\rm GHz}=4.1\pm 0.8 \times 10^{23}$ W Hz$^{-1}$ and integrated spectral index $\alpha = 1.15\pm 0.23$. We infer for this radio relic a lower bound on the magnetisation of $B\geq 0.4\, \mu$G, a shock Mach number of $\mathcal{M}\approx 3.79$, and a low acceleration efficiency consistent with DSA. This result suggests that a missing population of relics may become visible thanks to the unprecedented sensitivity of the new generation of radio telescopes.

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

The formation of cosmic structure culminates with the assembly of galaxy clusters, a process quite different from cluster to cluster. We present the study of the structure and dynamics of the Lyra complex formed of the two clusters RXC J1825.3+3026 and CIZA J1824.1+3029, very recently studied using both X-ray and radio data. This is the first analysis based on kinematics of member galaxies. New spectroscopic data for 285 galaxies were acquired at the Italian Telescopio Nazionale Galileo and used in combination with PanSTARRS photometry. The result of our member selection is a sample of 198 galaxies. For RXCJ1825 and CIZAJ1824 we report the redshifts, z=0.0645 and z=0.0708, the first estimates of velocity dispersion, sigmav=995 and 700 km/s, and dynamical mass, M200=1.1E15 and 4E14 Msun. The past assembly of RXCJ1825 is traced by the two dominant galaxies, both aligned with the major axis of the galaxy distribution along the East-West direction, and by a minor North-East substructure. We also detect a quite peculiar high velocity field in the South-West region of the Lyra complex. This feature is likely related to a high velocity, very luminous galaxy, suggested to be the central galaxy of a group in interaction with RXCJ1825 by very recent studies based on X-ray and radio data. The redshift of the whole Lyra complex is z=0.067. Assuming that the redshift difference between RXCJ1825 and CIZAJ1824 is due to the relative kinematics, the projected distance between the cluster centers is 1.3 Mpc and the los velocity difference is 1750 km/s. A dynamical analysis of the system shows that the two clusters are likely to be gravitationally bound, in a pre-merger phase, with CIZAJ1824 in front of RXCJ1825 and going toward it. Our results corroborate a picture where the Lyra region is the place of a very complex scenario of cluster assembly.

Nearby clusters of galaxies, z<0.1, are cosmic structures still under formation. Understanding the thermodynamic properties of merging clusters can provide crucial information on how they grow in the local universe. A detailed study of the intra-cluster medium (ICM) properties of un-relaxed systems is essential to understand the fate of in-falling structures and, more generally, the virialization process. We analyzed a mosaic of XMM-Newton observations (240 ks) of the Lyra system (z=0.067) that shows a complex dynamical state. We find the main cluster RXC J1825.3+3026 to be in a late merger phase, whereas its companion CIZA J1824.1+3029 is a relaxed cool-core cluster. We estimate a mass ratio of ~1:2 for the pair. No diffuse X-ray emission is found in the region between them, indicating that these clusters are in a pre-merger phase. We found evidence of a galaxy group infalling on RXC J1825.3+3026 in an advanced state of disruption. The Southern Galaxy, one of the brightest galaxies in the Lyra complex, was very likely at the center of the infalling group. This galaxy has a gaseous corona indicating that it was able to retain some of its gas after the ram-pressure stripping of the intra-group medium. In this scenario the diffuse emission excess observed southwest of RXC J1825.3+3026 could be due to gas once belonging to the group and/or to cluster ICM dislocated by the passage of the group. Finally, we identified three high-velocity galaxies aligned between RXC J1825.3+3026 and the SG, two of these showing evidence of gas stripped from them during infall. We estimate them to be currently falling onto the main cluster at an infall velocity of ~ 3000 km/s. Our study of the Lyra complex provides important clues about the processes presiding over the virialization of massive clusters in the local Universe.

With the opening of the X and gamma--ray windows in the sixties, thanks to to sounding rockets and satellite-borne instruments, extremely energetic and violent phenomena were discovered and subsequently found to be ubiquitous in the Universe. Observations in the high energy domain are fundamental for understanding how matter is organized and behaves around black holes; unravelling how these extreme objects influence their environments on a very large scale; and finding the still elusive obscured massive objects in the centre of galaxies. Other major problems in contemporary astrophysics, such as the understanding of acceleration processes at shocks of all sizes (those of pulsar wind nebulae, supernova remnants, but also at larger scales those of Active Galactic Nuclei radio lobes) in relation to the origin of cosmic-rays, or the definitive characterization of the debated non-thermal X-ray energy content of clusters of galaxies, also requires observations at very high energies. An observatory type medium mission operating from around 1 keV to about 600 keV can provide direct insights into these major questions. The essential characteristics will be coverage of the full energy range by telescopes featuring a large throughput and arc-second resolution optics, coupled to a compact focal plane assembly, with excellent imaging resolution and spectroscopy. In addition, the mission will provide unique polarimetry measurements in the hard X-ray domain, an important new diagnostic tool at energies for which the non-thermal processes dominate. The Polarimetric High-Energy Modular Telescope Observatory (PHEMTO) is designed to have performance several orders of magnitude better than the present hard X-ray instruments. This gives to PHEMTO the improvements in scientific performance needed for a mission in the 2050 era.

Diffuse radio emission associated with the intra-cluster medium (ICM) is observed in a number of merging galaxy clusters. It is currently believed that in mergers a fraction of the kinetic energy is channeled into non-thermal components, such as turbulence, cosmic rays and magnetic fields, that may lead to the formation of giant synchrotron sources in the ICM. Studying merging galaxy clusters in different evolutionary phases is fundamental to understanding the origin of radio emission in the ICM. We observed the nearby galaxy cluster pair RXC J1825.3+3026 ($z\sim0.065$) and CIZA J1824.1+3029 ($z\sim0.071$) at 120-168 MHz with the LOw Frequency ARray (LOFAR) and made use of a deep (240 ks) XMM-Newton dataset to study the non-thermal and thermal properties of the system. RXC J1825.3+3026 is in a complex dynamical state, with a primary on-going merger in the E-W direction and a secondary later stage merger with a group of galaxies in the SW, while CIZA J1824.1+3029 is dynamically relaxed. These two clusters are in a pre-merger phase. We report the discovery of a Mpc-scale radio halo with a low surface brightness extension in RXC J1825.3+3026 that follows the X-ray emission from the cluster center to the remnant of a galaxy group in the SW. This is among the least massive systems and the faintest giant radio halo known to date. Contrary to this, no diffuse radio emission is observed in CIZA J1824.1+3029 nor in the region between the pre-merger cluster pair. The power spectra of the X-ray surface brightness fluctuations of RXC J1825.3+3026 and CIZA J1824.1+3029 are in agreement with the findings for clusters exhibiting a radio halo and the ones where no radio emission has been detected, respectively. We provide quantitative support to the idea that cluster mergers play a crucial role in the generation of non-thermal components in the ICM.

The majority of the ordinary matter in the local Universe has been heated by strong structure formation shocks and resides in a largely unexplored hot, diffuse, X-ray emitting plasma that permeates the halos of galaxies, galaxy groups and clusters, and the cosmic web. We propose a next-generation "Cosmic Web Explorer" that will permit a complete and exhaustive understanding of these unseen baryons. This will be the first mission capable to reach the accretion shocks located several times farther than the virial radii of galaxy clusters, and reveal the out-of-equilibrium parts of the intra-cluster medium which are live witnesses to the physics of cosmic accretion. It will also enable a view of the thermal, kinematic, and chemical composition of the circum-galactic medium in galaxies with masses similar to the Milky Way, at the same level of detail that $Athena$ will unravel for the virialized regions of massive galaxy clusters, delivering a transformative understanding of the evolution of those galaxies in which most of the stars and metals in the Universe were formed. Finally, the proposed X-ray satellite will connect the dots of the large-scale structure by mapping, at high spectral resolution, as much as 100% of the diffuse gas hotter than $10^6$ K that fills the filaments of the cosmic web at low redshifts, down to an over-density of 1, both in emission and in absorption against the ubiquitous cosmic X-ray background, surveying up to half of the sky over 5 years in orbit. This requires a large effective area (~10 m$^2$ at 1 keV) over a large field of view (~1 deg$^2$), a megapixel cryogenic microcalorimeter array providing integral field spectroscopy with a resolving power $E/\Delta E$ = 2000 at 0.6 keV and a spatial resolution of 5" in the soft X-ray band, and a Low-Earth orbit that ensures a stable instrumental background and high sensitivity to faint, extended emission.

We present the results of a joint XMM-Newton and NuSTAR observation (200 ks) of the galaxy cluster Abell 523 at $z=0.104$. The peculiar morphology of the cluster radio halo and its outlier position in the radio power P(1.4 GHz) - X-ray luminosity plane make it an ideal candidate for the study of radio-X-ray correlations and for the search of inverse Compton (IC) emission. We constructed thermodynamic maps derived from the XMM observations to describe in detail the physical and dynamical state of the ICM. We performed a point-to-point comparison in terms of surface brightness between the X-ray and radio emissions, to quantify their morphological discrepancies. Making use of NuSTAR's hard X-ray focusing capability, we looked for IC emission both globally and locally, after modeling the purely thermal component with a multi-temperature description. The thermodynamic maps suggest the presence of a secondary merging process that could be responsible for the peculiar radio halo morphology. This hypothesis is supported by the comparison between the X-ray and radio surface brightnesses, which shows a broad intrinsic scatter and a series of outliers. The global NuSTAR spectrum can be explained by purely thermal gas emission, and there is no convincing evidence that an IC component is needed. The $3\sigma$ upper limit on the IC flux in the 20-80 keV band is in the range $\left[2.2 - 4.0\right] \times 10^{-13} \, \mathrm{erg} \, \mathrm{s}^{-1} \, \mathrm{cm}^{-2}$, implying a lower limit on the magnetic field strength in the range $B > [0.23 - 0.31] \, \mu G$. Locally, we looked for IC emission in the central region of the cluster radio halo finding a $3\sigma$ upper limit on the 20-80 keV non-thermal flux of $3.17 \times 10^{-14} \, \mathrm{erg} \, \mathrm{s}^{-1} \, \mathrm{cm}^{-2}$, corresponding to a lower limit on the magnetic field strength of $B \gtrsim 0.81 \, \mu G$.

We present the reconstruction of hydrostatic mass profiles in 13 X-ray luminous galaxy clusters that have been mapped in their X-ray and SZ signal out to $R_{200}$ for the XMM-Newton Cluster Outskirts Project (X-COP). Using profiles of the gas temperature, density and pressure that have been spatially resolved out to (median value) 0.9 $R_{500}$, 1.8 $R_{500}$, and 2.3 $R_{500}$, respectively, we are able to recover the hydrostatic gravitating mass profile with several methods and using different mass models. The hydrostatic masses are recovered with a relative (statistical) median error of 3% at $R_{500}$ and 6% at $R_{200}$. By using several different methods to solve the equation of the hydrostatic equilibrium, we evaluate some of the systematic uncertainties to be of the order of 5% at both $R_{500}$ and $R_{200}$. A Navarro-Frenk-White profile provides the best-fit in nine cases out of 13, with the remaining four cases that do not show a statistically significant tension with it. The distribution of the mass concentration follows the correlations with the total mass predicted from numerical simulations with a scatter of 0.18 dex, with an intrinsic scatter on the hydrostatic masses of 0.15 dex. We compare them with the estimates of the total gravitational mass obtained through X-ray scaling relations applied to $Y_X$, gas fraction and $Y_{SZ}$, and from weak lensing and galaxy dynamics techniques, and measure a substantial agreement with the results from scaling laws, from WL at both $R_{500}$ and $R_{200}$ (with differences below 15%), from cluster velocity dispersions, but a significant tension with the caustic masses that tend to underestimate the hydrostatic masses by 40% at $R_{200}$. We also compare these measurements with predictions from alternative models to the Cold Dark Matter, like the "Emergent Gravity" and MOND scenarios.

The hot plasma in galaxy clusters is expected to be heated to high temperatures through shocks and adiabatic compression. The thermodynamical properties of the gas encode information on the processes leading to the thermalization of the gas in the cluster's potential well as well as non-gravitational processes such as gas cooling, AGN feedback and kinetic energy. In this work we present the radial profiles of the thermodynamic properties of the intracluster medium (ICM) out to the virial radius for a sample of 12 galaxy clusters selected from the Planck all-sky survey. We determine the universal profiles of gas density, temperature, pressure, and entropy over more than two decades in radius. We exploit jointly X-ray information from XMM and Sunyaev-Zel'dovich constraints from Planck to recover thermodynamic properties out to 2 R500. We provide average functional forms for the radial dependence of the main quantities and quantify the slope and intrinsic scatter of the population as a function of radius. We find that gas density and pressure profiles steepen steadily with radius, in excellent agreement with previous observational results. Entropy profiles beyond R500 closely follow the predictions for the gravitational collapse of structures. The scatter in all thermodynamical quantities reaches a minimum in the range [0.2-0.8] R500 and increases outwards. Somewhat surprisingly, we find that pressure is substantially more scattered than temperature and density. Our results indicate that once accreting substructures are properly excised, the properties of the ICM beyond the cooling region R > 0.3 R500) follow remarkably well the predictions of simple gravitational collapse and require little non-gravitational corrections.

Galaxy clusters are the endpoints of structure formation and are continuously growing through the merging and accretion of smaller structures. Numerical simulations predict that a fraction of their energy content is not yet thermalized, mainly in the form of kinetic motions (turbulence, bulk motions). Measuring the level of non-thermal pressure support is necessary to understand the processes leading to the virialization of the gas within the potential well of the main halo and to calibrate the biases in hydrostatic mass estimates. We present high-quality measurements of hydrostatic masses and intracluster gas fraction out to the virial radius for a sample of 12 nearby clusters with available XMM-Newton and Planck data. We compare our hydrostatic gas fractions with the expected universal gas fraction to constrain the level of non-thermal pressure support. We find that hydrostatic masses require little correction and infer a median non-thermal pressure fraction of $\sim6\%$ and $\sim10\%$ at $R_{500}$ and $R_{200}$, respectively. Our values are lower than the expectations of hydrodynamical simulations, possibly implying a faster thermalization of the gas. If instead we use the mass calibration adopted by the Planck team, we find that the gas fraction of massive local systems implies a mass bias $1-b=0.85\pm0.05$ for SZ-derived masses, with some evidence for a mass-dependent bias. Conversely, the high bias required to match Planck CMB and cluster count cosmology is excluded by the data at high significance, unless the most massive halos are missing a substantial fraction of their baryons.

The Shapley Concentration ($z\approx0.048$) covers several degrees in the Southern Hemisphere, and includes galaxy clusters in advanced evolutionary stage, groups of clusters in the early stages of merger, fairly massive clusters with ongoing accretion activity, and smaller groups located in filaments in the regions between the main clusters. With the goal to investigate the role of cluster mergers and accretion on the radio galaxy population, we performed a multi-wavelength study of the BCGs and of the galaxies showing extended radio emission in the cluster complexes of Abell 3528 and Abell 3558. Our study is based on a sample of 12 galaxies. We observed the clusters with the GMRT at 235, 325 and 610 MHz, and with the VLA at 8.46 GHz. We complemented our study with the TGSS at 150 MHz, the SUMSS at 843 MHz and ATCA at 1380, 1400, 2380, and 4790 MHz data. Optical imaging with ESO-VST and mid-IR coverage with WISE are also available for the host galaxies. We found deep differences in the properties of the radio emission of the BCGs in the two cluster complexes. The BCGs in the A3528 complex and in A3556, which are relaxed cool-core objects, are powerful active radio galaxies. They also present hints of restarted activity. On the contrary, the BCGs in A3558 and A3562, which are well known merging systems, are very faint, or quiet, in the radio band. The optical and IR properties of the galaxies are fairly similar in the two complexes, showing all passive red galaxies. Our study shows remarkable differences in the radio properties of the BGCs, which we relate to the different dynamical state of the host cluster. On the contrary, the lack of changes between such different environments in the optical band suggests that the dynamical state of galaxy clusters does not affect the optical counterparts of the radio galaxies, at least over the life-time of the radio emission.

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.

In recent years, significant progress has been made in building new galaxy clusters samples, at low and high redshifts, from wide-area surveys, particularly exploiting the Sunyaev--Zel'dovich (SZ) effect. A large effort is underway to identify and characterize these new systems with optical/NIR and X-ray facilities, thus opening new avenues to constraint cosmological models using structure growth and geometrical tests. A census of galaxy clusters sets constraints on reionization mechanisms and epochs, which need to be reconciled with recent limits on the reionization optical depth from cosmic microwave background (CMB) experiments. Future advances in SZ effect measurements will include the possibility to (unambiguously) measure directly the kinematic SZ effect, to build an even larger catalogue of galaxy clusters able to study the high redshift universe, and to make (spatially-)resolved galaxy cluster maps with even spectral capability to (spectrally-)resolve the relativistic corrections of the SZ effect.

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.

In the local Universe, the growth of massive galaxy clusters mainly operates through the continuous accretion of group-scale systems. The infalling group in Abell 2142 is the poster child of such an accreting group, and as such, it is an ideal target to study the astrophysical processes induced by structure formation. We present the results of a deep (200 ks) observation of this structure with Chandra, which highlights the complexity of this system in exquisite detail. In the core of the group, the spatial resolution of Chandra reveals the presence of a leading edge and a complex AGN-induced activity. The morphology of the stripped gas tail appears straight in the innermost 250 kpc, suggesting that magnetic draping efficiently shields the gas from its surroundings. However, beyond $\sim300$ kpc from the core, the tail flares and the morphology becomes strongly irregular, which could be explained by a breaking of the drape, e.g. because of turbulent motions. The power spectrum of surface-brightness fluctuations is relatively flat ($P_{2D}\propto k^{-2.3}$), which indicates that thermal conduction is strongly inhibited even beyond the region where magnetic draping is effective. The amplitude of density fluctuations in the tail is consistent with a mild level of turbulence with a Mach number $M_{3D}\sim0.1-0.25$. Overall, our results show that the processes leading to the thermalization and mixing of the infalling gas are slow and relatively inefficient.

XMM-Newton is the direct precursor of the future ESA ATHENA mission. A study of its particle-induced background provides therefore significant insight for the ATHENA mission design. We make use of about 12 years of data, products from the third XMM-Newton catalog as well as FP7 EXTraS project to avoid celestial sources contamination and to disentangle the different components of the XMM-Newton particle-induced background. Within the ESA R&D AREMBES collaboration, we built new analysis pipelines to study the different components of this background: this covers time behavior as well as spectral and spatial characteristics.

We present an accurate characterization of the particle background behaviour on XMM-Newton based on the entire EPIC archive. This corresponds to the largest EPIC data set ever examined. Our results have been obtained thanks to the collaboration between the FP7 European program EXTraS and the ESA R&D ATHENA activity AREMBES. We used as a diagnostic an improved version of the diagnostic which compares the data collected in unexposed region of the detector with the region of the field of view in the EPIC-MOS. We will show that the in Field-of-View excess background is made up of two different components, one associated to flares produced by soft protons and the other one to a low-intensity background. Its origin needs to be further investigated.

We show the results obtained in the FP7 European program EXTraS and in the ESA R&D ATHENA activity AREMBES aimed at a deeper understanding of the XMM-Newton background to better design the ATHENA mission. Thanks to an analysis of the full EPIC archive coupled to the information obtained by the Radiation Monitor we show the cosmic ray origin of the unfocused particle background and its anti-correlation with the solar activity. We show the first results of the effort to obtain informations about the particle component of the soft proton focused background.

A detailed characterization of the particle induced background is fundamental for many of the scientific objectives of the Athena X-ray telescope, thus an adequate knowledge of the background that will be encountered by Athena is desirable. Current X-ray telescopes have shown that the intensity of the particle induced background can be highly variable. Different regions of the magnetosphere can have very different environmental conditions, which can, in principle, differently affect the particle induced background detected by the instruments. We present results concerning the influence of the magnetospheric environment on the background detected by EPIC instrument onboard XMM-Newton through the estimate of the variation of the in-Field-of-View background excess along the XMM-Newton orbit. An important contribution to the XMM background, which may affect the Athena background as well, comes from soft proton flares. Along with the flaring component a low-intensity component is also present. We find that both show modest variations in the different magnetozones and that the soft proton component shows a strong trend with the distance from Earth.

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

We study the intra-cluster magnetic field in the poor galaxy cluster Abell 194 by complementing radio data, at different frequencies, with data in the optical and X-ray bands. We analyze new total intensity and polarization observations of Abell 194 obtained with the Sardinia Radio Telescope (SRT). We use the SRT data in combination with archival Very Large Array observations to derive both the spectral aging and Rotation Measure (RM) images of the radio galaxies 3C40A and 3C40B embedded in Abell 194. The optical analysis indicates that Abell 194 does not show a major and recent cluster merger, but rather agrees with a scenario of accretion of small groups. Under the minimum energy assumption, the lifetimes of synchrotron electrons in 3C40B measured from the spectral break are found to be 157 Myrs. The break frequency image and the electron density profile inferred from the X-ray emission are used in combination with the RM data to constrain the intra-cluster magnetic field power spectrum. By assuming a Kolmogorov power law power spectrum, we find that the RM data in Abell 194 are well described by a magnetic field with a maximum scale of fluctuations of Lambda_max=64 kpc and a central magnetic field strength of <B0>=1.5 microG. Further out, the field decreases with the radius following the gas density to the power of eta=1.1. Comparing Abell 194 with a small sample of galaxy clusters, there is a hint of a trend between central electron densities and magnetic field strengths.

We report on a spectral study at radio frequencies of the giant radio halo in A2142 (z=0.0909), which we performed to explore its nature and origin. A2142 is not a major merger and the presence of a giant radio halo is somewhat surprising. We performed deep radio observations with the GMRT at 608 MHz, 322 MHz, and 234 MHz and with the VLA in the 1-2 GHz band. We obtained high-quality images at all frequencies in a wide range of resolutions. The radio halo is well detected at all frequencies and extends out to the most distant cold front in A2142. We studied the spectral index in two regions: the central part of the halo and a second region in the direction of the most distant south-eastern cold front, selected to follow the bright part of the halo and X-ray emission. We complemented our observations with a preliminary LOFAR image at 118 MHz and with the re-analysis of archival VLA data at 1.4 GHz. The two components of the radio halo show different observational properties. The central brightest part has higher surface brightess and a spectrum whose steepness is similar to those of the known radio halos, i.e. $\alpha^{\rm 1.78~GHz}_{\rm 118~MHz}=1.33\pm 0.08$. The ridge, which fades into the larger scale emission, is broader in size and has considerably lower surface brightess and a moderately steeper spectrum, i.e. $\alpha^{\rm 1.78~GHz}_{\rm 118~MHz}\sim 1.5$. We propose that the brightest part of the radio halo is powered by the central sloshing in A2142, similar to what has been suggested for mini-halos, or by secondary electrons generated by hadronic collisions in the ICM. On the other hand, the steeper ridge may probe particle re-acceleration by turbulence generated either by stirring the gas and magnetic fields on a larger scale or by less energetic mechanisms, such as continuous infall of galaxy groups or an off-axis merger.

We characterized the population of galaxy clusters detected with the SZ effect with Planck, by measuring the cool core state of the objects in a well-defined subsample of the Planck catalogue. We used as indicator the concentration parameter Santos et al. (2008). The fraction of cool core clusters is $29 \pm 4 \%$ and does not show significant indications of evolution in the redshift range covered by our sample. We compare the distribution of the concentration parameter in the Planck sample with the one of the X-ray selected sample MACS (Mann & Ebeling, 2011): the distributions are significantly different and the cool core fraction in MACS is much higher ($59 \pm 5 \%$). Since X-ray selected samples are known to be biased towards cool cores due to the presence of their prominent surface brightness peak, we simulated the impact of the "cool core bias" following Eckert et al. (2011). We found that it plays a large role in the difference between the fractions of cool cores in the two samples. We examined other selection effects that could in principle bias SZ-surveys against cool cores but we found that their impact is not sufficient to explain the difference between Planck and MACS. The population of X-ray under-luminous objects, which are found in SZ-surveys but missing in X-ray samples (Planck Collaboration 2016), could possibly contribute to the difference, as we found most of them to be non cool cores, but this hypothesis deserves further investigation.

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

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

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

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

We investigate the relationship between X-ray cooling and star formation in brightest cluster galaxies (BCGs). We present an X-ray spectral analysis of the inner regions, 10-40 kpc, of six nearby cool core clusters (z<0.35) observed with Chandra ACIS. This sample is selected on the basis of the high star formation rate (SFR) observed in the BCGs. We restrict our search for cooling gas to regions that are roughly cospatial with the starburst. We fit single- and multi-temperature mkcflow models to constrain the amount of isobarically cooling intracluster medium (ICM). We find that in all clusters, below a threshold temperature ranging between 0.9 and 3 keV, only upper limits can be obtained. In four out of six objects, the upper limits are significantly below the SFR and in two, namely A1835 and A1068, they are less than a tenth of the SFR. Our results suggests that a number of mechanisms conspire to hide the cooling signature in our spectra. In a few systems the lack of a cooling signature may be attributed to a relatively long delay time between the X-ray cooling and the star burst. However, for A1835 and A1068, where the X-ray cooling time is shorter than the timescale of the starburst, a possible explanation is that the region where gas cools out of the X-ray phase extends to very large radii, likely beyond the core of these systems.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

In the current epoch, one of the main mechanisms driving the growth of galaxy clusters is the continuous accretion of group-scale halos. In this process, the ram pressure applied by the hot intracluster medium on the gas content of the infalling group is responsible for stripping the gas from its dark-matter halo, which gradually leads to the virialization of the infalling gas in the potential well of the main cluster. Using deep wide-field observations of the poor cluster Hydra A/A780 with XMM-Newton and Suzaku, we report the discovery of an infalling galaxy group 1.1 Mpc south of the cluster core. The presence of a substructure is confirmed by a dynamical study of the galaxies in this region. A wake of stripped gas is trailing behind the group over a projected scale of 760 kpc. The temperature of the gas along the wake is constant at kT ~ 1.3 keV, which is about a factor of two less than the temperature of the surrounding plasma. We observe a cold front pointing westwards compared to the peak of the group, which indicates that the group is currently not moving in the direction of the main cluster, but is moving along an almost circular orbit. The overall morphology of the group bears remarkable similarities with high-resolution numerical simulations of such structures, which greatly strengthens our understanding of the ram-pressure stripping process.

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

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

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

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

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

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

We want to characterize the dynamical state of galaxy clusters detected with the Sunyaev-Zeldovich (SZ) effect by Planck and compare them with the dynamical state of clusters selected in X-rays survey. We analyzed a representative subsample of the Planck SZ catalogue, containing the 132 clusters with the highest signal to noise ratio and characterize their dynamical state using as indicator the projected offset between the peak of the X-ray emission and the position of the Brightest cluster galaxy. We study the distribution of our indicator in our sample and compare it to its distribution in X-ray selected samples (HIFLUGCS, MACS and REXCESS). The distributions are significantly different and the fraction of relaxed objects is smaller in the Planck sample ($52 \pm 4 \%$) than in X-ray samples ($\simeq 74\%$) We interpret this result as an indication of different selection effects affecting X-rays (e.g. "cool core bias") and SZ surveys of galaxy clusters.

Our goal is to provide a robust estimate of the metal content of the ICM in massive clusters. We make use of published abundance profiles for a sample of ~60 nearby systems, we include in our estimate uncertainties associated to the measurement process and to the almost total lack of measures in cluster outskirts. We perform a first, albeit rough, census of metals finding that the mean abundance of the ICM within r_180 is very poorly constrained, 0.06Z_sol < Z < 0.26Z_sol, and presents no tension with expectations. Similarly, the question of if and how the bulk of the metal content in clusters varies with cosmic time, is very much an open one. A solid estimate of abundances in cluster outskirts could be achieved by combining observations of the two experiments which will operate on board Athena, the XIFU and the WFI, provided they do not fall victim to the de-scoping process that has afflicted several space observatories over the last decade.

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

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.

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

We present the analysis of a local (z = 0.04 - 0.2) sample of 31 galaxy clusters with the aim of measuring the density of the X-ray emitting gas in cluster outskirts. We compare our results with numerical simulations to set constraints on the azimuthal symmetry and gas clumping in the outer regions of galaxy clusters. We exploit the large field-of-view and low instrumental background of ROSAT/PSPC to trace the density of the intracluster gas out to the virial radius. We perform a stacking of the density profiles to detect a signal beyond r200 and measure the typical density and scatter in cluster outskirts. We also compute the azimuthal scatter of the profiles with respect to the mean value to look for deviations from spherical symmetry. Finally, we compare our average density and scatter profiles with the results of numerical simulations. As opposed to some recent Suzaku results, and confirming previous evidence from ROSAT and Chandra, we observe a steepening of the density profiles beyond \sim r500. Comparing our density profiles with simulations, we find that non-radiative runs predict too steep density profiles, whereas runs including additional physics and/or treating gas clumping are in better agreement with the observed gas distribution. We report for the first time the high-confidence detection of a systematic difference between cool-core and non-cool core clusters beyond \sim 0.3r200, which we explain by a different distribution of the gas in the two classes. Beyond \sim r500, galaxy clusters deviate significantly from spherical symmetry, with only little differences between relaxed and disturbed systems. We find good agreement between the observed and predicted scatter profiles, but only when the 1% densest clumps are filtered out in the simulations. [Abridged]

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.