The cosmic infrared background (CIB) contains emissions accumulated over the entire history of the Universe, including from objects inaccessible to individual telescopic studies. The near-IR (~1-10 mic) part of the CIB, and its fluctuations, reflects emissions from nucleosynthetic sources and gravitationally accreting black holes (BHs). If known galaxies are removed to sufficient depths the source-subtracted CIB fluctuations at near-IR can reveal sources present in the first-stars-era and possibly new stellar populations at more recent times. This review discusses the recent progress in this newly emerging field which identified, with new data and methodology, significant source-subtracted CIB fluctuations substantially in excess of what can be produced by remaining known galaxies. The CIB fluctuations further appear coherent with unresolved cosmic X-ray background (CXB) indicating a very high fraction of BHs among the new sources producing the CIB fluctuations. These observations have led to intensive theoretical efforts to explain the measurements and their properties. While current experimental configurations have limitations in decisively probing these theories, their potentially remarkable implications will be tested in the upcoming CIB measurements with the ESA's Euclid dark energy mission. We describe the goals and methodologies of LIBRAE (Looking at Infrared Background Radiation with Euclid), a NASA-selected project for CIB science with Euclid, which has the potential for transforming the field into a new area of precision cosmology.

The correlation of weak lensing and Cosmic Microwave Anisotropy (CMB) data traces the pressure distribution of the hot, ionized gas and the underlying matter density field. The measured correlation is dominated by baryons residing in halos. Detecting the contribution from unbound gas by measuring the residual cross-correlation after masking all known halos requires a theoretical understanding of this correlation and its dependence with model parameters. Our model assumes that the gas in filaments is well described by a log-normal probability distribution function, with temperatures $10^{5-7}$K and overdensities $\xi\le 100$. The lensing-comptonization cross-correlation is dominated by gas with overdensities in the range $\xi\approx[3-33]$; the signal is generated at redshifts $z\le 1$. If only 10\% of the measured cross-correlation is due to unbound gas, then the most recent measurements set an upper limit of $\bar{T}_e\lesssim 10^6$K on the mean temperature of Inter Galactic Medium. The amplitude is proportional to the baryon fraction stored in filaments. The lensing-comptonization power spectrum peaks at a different scale than the gas in halos making it possible to distinguish both contributions. To trace the distribution of the low density and low temperature plasma on cosmological scales, the effect of halos will have to be subtracted from the data, requiring observations with larger signal-to-noise ratio than currently available.

We determine the relation between the Comptonization parameter predicted using X-ray data $Y_{C,Xray}$ and the X-ray luminosity $L_X$, both magnitudes derived from ROSAT data, with the Comptonization parameter $Y_{C,SZ}$ measured on {\it Planck} 2013 foreground cleaned Nominal maps. The 560 clusters of our sample includes clusters with masses $M\ge 10^{13}M_\odot$, one order of magnitude smaller than those used by the Planck Collaboration in a similar analysis. It also contains eight times more clusters in the redshift interval $z\le 0.3$. The prediction of the $\beta=2/3$ model convolved with the Planck antenna beam agrees with the anisotropies measured in foreground cleaned Planck Nominal maps within the X-ray emitting region, confirming the results of an earlier analysis (Atrio-Barandela et al. 2008). The universal pressure profile overestimates the signal by a 15-21\% depending on the angular aperture. We show that the discrepancy is not due to the presence of {\it cool-core} systems but it is an indication of a brake in the $L_X-M$ relation towards low mass systems. We show that relation of the Comptonization parameter averaged over the region that emits 99\% of the X-ray flux and and the X-ray luminosity is consistent with the predictions of the self-similar model. We confirm previous findings that the scaling relations studied here do not evolve with redshift within the range probed by our catalog.

Models where Dark Matter and Dark Energy interact with each other have been proposed to solve the coincidence problem. We review the motivations underlying the need to introduce such interaction, its influence on the background dynamics and how it modifies the evolution of linear perturbations. We test models using the most recent observational data and we find that the interaction is compatible with the current astronomical and cosmological data. Finally, we describe the forthcoming data sets from current and future facilities that are being constructed or designed that will allow a clearer understanding of the physics of the dark sector.

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

The Euclid space mission, designed to probe evolution of the Dark Energy, will map a large area of the sky at three adjacent near-IR filters, Y, J and H. This coverage will also enable mapping source-subtracted cosmic infrared background (CIB) fluctuations with unprecedented accuracy on sub-degree angular scales. Here we propose methodology, using the Lyman-break tomography applied to the Euclid-based CIB maps, to accurately isolate the history of CIB emissions as a function of redshift from 10 < z < 20, and to identify the baryonic acoustic oscillations (BAOs) at those epochs. To identify the BAO signature, we would assemble individual CIB maps over conservatively large contiguous areas of >~ 400 sq deg. The method can isolate the CIB spatial spectrum by z to sub-percent statistical accuracy. We illustrate this with a specific model of CIB production at high z normalized to reproduce the measured Spitzer-based CIB fluctuation. We show that even if the latter contain only a small component from high-z sources, the amplitude of that component can be accurately isolated with the methodology proposed here and the BAO signatures at z>~ 10 are recovered well from the CIB fluctuation spatial spectrum. Probing the BAO at those redshifts will be an important test of the underlying cosmological paradigm, and would narrow the overall uncertainties on the evolution of cosmological parameters, including the Dark Energy. Similar methodology is applicable to the planned WFIRST mission, where we show that a possible fourth near-IR channel at > 2 micron would be beneficial.

Our velocity relative to the rest frame of the cosmic microwave background (CMB) generates a dipole temperature anisotropy on the sky which has been well measured for more than 30 years, and has an accepted amplitude of v/c = 0.00123, or v = 369km/s. In addition to this signal generated by Doppler boosting of the CMB monopole, our motion also modulates and aberrates the CMB temperature fluctuations (as well as every other source of radiation at cosmological distances). This is an order 0.1% effect applied to fluctuations which are already one part in roughly one hundred thousand, so it is quite small. Nevertheless, it becomes detectable with the all-sky coverage, high angular resolution, and low noise levels of the Planck satellite. Here we report a first measurement of this velocity signature using the aberration and modulation effects on the CMB temperature anisotropies, finding a component in the known dipole direction, (l,b)=(264, 48) [deg], of 384km/s +- 78km/s (stat.) +- 115km/s (syst.). This is a significant confirmation of the expected velocity.

We constrain the deviation of adiabatic evolution of the Universe using the data on the Cosmic Microwave Background (CMB) temperature anisotropies measured by the {\it Planck} satellite and a sample of 481 X-ray selected clusters with spectroscopically measured redshifts. To avoid antenna beam effects, we bring all the maps to the same resolution. We use a CMB template to subtract the cosmological signal while preserving the Thermal Sunyaev-Zeldovich (TSZ) anisotropies; next, we remove galactic foreground emissions around each cluster and we mask out all known point sources. If the CMB black-body temperature scales with redshift as $T(z)=T_0(1+z)^{1-\alpha}$, we constrain deviations of adiabatic evolution to be $\alpha=-0.007\pm 0.013$, consistent with the temperature-redshift relation of the standard cosmological model. This result could suffer from a potential bias $\delta\alpha$ associated with the CMB template, that we quantify it to be $|\delta\alpha|\le 0.02$ and with the same sign than the measured value of $\alpha$, but is free from those biases associated with using TSZ selected clusters; it represents the best constraint to date of the temperature-redshift relation of the Big-Bang model using only CMB data, confirming previous results.

Planck data when combined with ancillary data provide a unique opportunity to separate the diffuse emission components of the inner Galaxy. The purpose of the paper is to elucidate the morphology of the various emission components in the strong star-formation region lying inside the solar radius and to clarify the relationship between the various components. The region of the Galactic plane covered is l=300-0-60deg where star-formation is highest and the emission is strong enough to make meaningful component separation. The latitude widths in this longitude range lie between 1deg and 2deg, which correspond to FWHM z-widths of 100-200pc at a typical distance of 6kpc. The four emission components studied here are synchrotron, free-free, anomalous microwave emission (AME), and thermal (vibrational) dust emission. These components are identified by constructing spectral energy distributions (SEDs) at positions along the Galactic plane using the wide frequency coverage of Planck (28.4-857GHz) in combination with low-frequency radio data at 0.408-2.3GHz plus WMAP data at 23-94GHz, along with far-infrared (FIR) data from DIRBE and IRAS. The free-free component is determined from radio recombination line (RRL) data. AME is found to be comparable in brightness to the free-free emission on the Galactic plane in the frequency range 20-40GHz with a width in latitude similar to that of the thermal dust; it comprises 45+/-1% of the total 28.4GHz emission in the longitude range l=300-0-60deg. The free-free component is the narrowest, reflecting the fact that it is produced by current star-formation as traced by the narrow distribution of OB stars. It is the dominant emission on the plane between 60 and 100GHz. RRLs from this ionized gas are used to assess its distance, leading to a free-free z-width of FWHM ~100pc...(abridged)

Planck has mapped the intensity and polarization of the sky at microwave frequencies with unprecedented sensitivity. We make use of the Planck 353 GHz I, Q, and U Stokes maps as dust templates, and cross-correlate them with the Planck and WMAP data at 12 frequencies from 23 to 353 GHz, over circular patches with 10 degree radius. The cross-correlation analysis is performed for both intensity and polarization data in a consistent manner. We use a mask that focuses our analysis on the diffuse interstellar medium at intermediate Galactic latitudes. We determine the spectral indices of dust emission in intensity and polarization between 100 and 353 GHz, for each sky-patch. The mean values, $1.59\pm0.02$ for polarization and $1.51\pm0.01$ for intensity, for a mean dust temperature of 19.6 K, are close, but significantly different ($3.6\,\sigma$). We determine the mean spectral energy distribution (SED) of the microwave emission, correlated with the 353 GHz dust templates, by averaging the results of the correlation over all sky-patches. We find that the mean SED increases for decreasing frequencies at $\nu < 60$ GHz, for both intensity and polarization. The rise of the polarization SED towards low frequencies may be accounted for by a synchrotron component correlated with dust, with no need for any polarization of the anomalous microwave emission. We use a spectral model to separate the synchrotron and dust polarization and to characterize the spectral dependence of the dust polarization fraction. The polarization fraction ($p$) of the dust emission decreases by $(21\pm6)$ % from 353 to 70 GHz. The decrease of $p$ could indicate differences in polarization efficiency among components of interstellar dust (e.g., carbon versus silicate grains). Our observational results provide inputs to quantify and optimize the separation between Galactic and cosmological polarization.

The Planck survey provides unprecedented full-sky coverage of the submillimetre polarized emission from Galactic dust, bringing new constraints on the properties of dust. The dust grains that emit the radiation seen by Planck in the submillimetre also extinguish and polarize starlight in the visible. Comparison of the polarization of the emission and of the interstellar polarization on selected lines of sight probed by stars provides unique new diagnostics of the emission and light scattering properties of dust. Using ancillary catalogues of interstellar polarization and extinction of starlight, we obtain the degree of polarization, $p_V$ , and the optical depth in the V band to the star, $\tau_V$. Toward these stars we measure the submillimetre polarized intensity, $P_S$, and total intensity, $I_S$, in the Planck 353 GHz channel. For those lines of sight through the diffuse interstellar medium with comparable values of the estimated column density and polarization directions close to orthogonal, we correlate properties in the submillimetre and visible to find two ratios, $R_{S/V} = (P_S/I_S)/(p_V/\tau_V)$ and $R_{P/p} = P_S/p_V$ , the latter focusing directly on the polarization properties of the aligned grain population alone. We find $R_{S/V}$ = 4.2, with statistical and systematic uncertainties 0.2 and 0.3, respectively, and $R_{P/p}$ = 5.4 MJy sr$^{-1}$, with uncertainties 0.2 and 0.3 MJy sr$^{-1}$, respectively. Our estimate of $R_{S/V}$ is compatible with predictions based on a range of polarizing dust models that have been developed for the diffuse interstellar medium. However, our estimate of $R_{P/p}$ is not compatible with predictions, which are too low by a factor of about 2.5. This more discriminating diagnostic, $R_{P/p}$, indicates that changes to the optical properties in the models of the aligned grain population are required.

We analyse the implications of the Planck data for cosmic inflation. The Planck nominal mission temperature anisotropy measurements, combined with the WMAP large-angle polarization, constrain the scalar spectral index to $n_s = 0.9603 \pm 0.0073$, ruling out exact scale invariance at over 5 $\sigma$. Planck establishes an upper bound on the tensor-to-scalar ratio of r < 0.11 (95% CL). The Planck data thus shrink the space of allowed standard inflationary models, preferring potentials with V" < 0. Exponential potential models, the simplest hybrid inflationary models, and monomial potential models of degree n > 2 do not provide a good fit to the data. Planck does not find statistically significant running of the scalar spectral index, obtaining $d n_s/d ln k = -0.0134 \pm 0.0090$. Several analyses dropping the slow-roll approximation are carried out, including detailed model comparison and inflationary potential reconstruction. We also investigate whether the primordial power spectrum contains any features. We find that models with a parameterized oscillatory feature improve the fit $\chi^2$ by ~ 10; however, Bayesian evidence does not prefer these models. We constrain several single-field inflation models with generalized Lagrangians by combining power spectrum data with bounds on $f_\mathrm{NL}$ measured by Planck. The fractional primordial contribution of CDM isocurvature modes in the curvaton and axion scenarios has upper bounds of 0.25% or 3.9% (95% CL), respectively. In models with arbitrarily correlated CDM or neutrino isocurvature modes, an anticorrelation can improve $\chi^2$ by approximatively 4 as a result of slightly lowering the theoretical prediction for the $\ell<40$ multipoles relative to the higher multipoles. Nonetheless, the data are consistent with adiabatic initial conditions.

We update the all-sky Planck catalogue of 1227 clusters and cluster candidates (PSZ1) published in March 2013, derived from Sunyaev-Zeldovich (SZ) effect detections using the first 15.5 months of Planck satellite observations. Addendum. We deliver an updated version of the PSZ1 catalogue, reporting the further confirmation of 86 Planck-discovered clusters. In total, the PSZ1 now contains 947 confirmed clusters, of which 214 were confirmed as newly discovered clusters through follow-up observations undertaken by the Planck Collaboration. The updated PSZ1 contains redshifts for 913 systems, of which 736 (~80.6%) are spectroscopic, and associated mass estimates derived from the Y_z mass proxy. We also provide a new SZ quality flag, derived from a novel artificial neural network classification of the SZ signal, for the remaining 280 candidates. Based on this assessment, the purity of the updated PSZ1 catalogue is estimated to be 94%. In this release, we provide the full updated catalogue and an additional readme file with further information on the Planck SZ detections.

Analytical $f(R)$-gravity models introduce Yukawa-like corrections to the Newtonian potential in the weak field limit. These models can explain the dynamics of galaxies and cluster of galaxies without requiring dark matter. To test the model, we have computed the pressure profile of 579 X-ray galaxy clusters assuming the gas is in hydrostatic equilibrium within the potential well of the modified gravitational potential. We have compared those profiles with the ones measured in the foreground cleaned SMICA released by the Planck Collaboration. Our results show that Extended Theories of Gravity explain the dynamics of self-gravitating systems at cluster scales and represent an alternative to dark matter haloes.

We have constructed the first all-sky map of the thermal Sunyaev-Zeldovich (tSZ) effect by applying specifically tailored component separation algorithms to the 100 to 857 GHz frequency channel maps from the Planck survey. These maps show an obvious galaxy cluster tSZ signal that is well matched with blindly detected clusters in the Planck SZ catalogue. To characterize the signal in the tSZ map we have computed its angular power spectrum. At large angular scales ($\ell < 60$), the major foreground contaminant is the diffuse thermal dust emission. At small angular scales ($\ell > 500$) the clustered Cosmic Infrared Background (CIB) and residual point sources are the major contaminants. These foregrounds are carefully modelled and subtracted. We measure the tSZ power spectrum in angular scales, $0.17^{\circ} \lesssim \theta \lesssim 3.0^{\circ}$, that were previously unexplored. The measured tSZ power spectrum is consistent with that expected from the Planck catalogue of SZ sources, with additional clear evidence of signal from unresolved clusters and, potentially, diffuse warm baryons. We use the tSZ power spectrum to obtain the following cosmological constraints: $\sigma_8(\Omega_{\mathrm{m}}/0.28)^{3.2/8.1}=0.784 \pm 0.016 (68% C.L.). Marginalized band-powers of the Planck tSZ power spectrum and the best-fit model are given. The non-Gaussianity of the Compton parameter map is further characterized by computing its 1D probability distribution function and its bispectrum. These are used to place additional independent constraints on $\sigma_{8}$.

Late stages of stellar evolution are characterized by copious mass-loss events whose signature is the formation of circumstellar envelopes (CSE). Planck multi-frequency measurements have provided relevant information on a sample of Galactic planetary nebulae (PNe) in the important and relatively unexplored observational band between 30 and 857GHz. Planck enables the assembly of comprehensive PNe spectral energy distributions (SEDs) from radio {\bf to} far-infrared frequencies. Modelling of the derived SEDs provides us with information on physical properties of CSEs and the mass content of both main components: ionised gas, traced by the free-free emission at cm--mm waves; and thermal dust, traced by the millimetre and far-IR emission. In particular, the amount of ionised gas and dust has been derived here. Such quantities have also been estimated for the very young PN CRL618, where the strong variability observed in its radio and millimetre emission has previously prevented the construction of its SED. A morphological study of the Helix Nebula has also been performed. Planck maps reveal, for the first time, the spatial distribution of the dust inside the envelope, allowing us to identify different components, the most interesting of which is a very extended component (up to 1pc) that may be related to a region where the slow expanding envelope is interacting with the surrounding interstellar medium.

Polarized emission observed by Planck HFI at 353 GHz towards a sample of nearby fields is presented, focusing on the statistics of polarization fractions $p$ and angles $\psi$. The polarization fractions and column densities in these nearby fields are representative of the range of values obtained over the whole sky. We find that: (i) the largest polarization fractions are reached in the most diffuse fields; (ii) the maximum polarization fraction $p_\mathrm{max}$ decreases with column density $N_\mathrm{H}$ in the more opaque fields with $N_\mathrm{H} > 10^{21}\,\mathrm{cm}^{-2}$; and (iii) the polarization fraction along a given line of sight is correlated with the local spatial coherence of the polarization angle. These observations are compared to polarized emission maps computed in simulations of anisotropic magnetohydrodynamical (MHD) turbulence in which we assume a uniform intrinsic polarization fraction of the dust grains. We find that an estimate of this parameter may be recovered from the maximum polarization fraction $p_\mathrm{max}$ in diffuse regions where the magnetic field is ordered on large scales and perpendicular to the line of sight. This emphasizes the impact of anisotropies of the magnetic field on the emerging polarization signal. The decrease of the polarization fraction with column density in nearby molecular clouds is well reproduced in the simulations, indicating that it is essentially due to the turbulent structure of the magnetic field: an accumulation of variously polarized structures along the line of sight leads to such an anti-correlation. In the simulations, polarization fractions are also found to anti-correlate with the angle dispersion function $\mathcal{S}$. [abridged]

This paper presents an all-sky model of dust emission from the Planck 857, 545 and 353 GHz, and IRAS 100 micron data. Using a modified black-body fit to the data we present all-sky maps of the dust optical depth, temperature, and spectral index over the 353-3000 GHz range. This model is a tight representation of the data at 5 arc min. It shows variations of the order of 30 % compared with the widely-used model of Finkbeiner, Davis, and Schlegel. The Planck data allow us to estimate the dust temperature uniformly over the whole sky, providing an improved estimate of the dust optical depth compared to previous all-sky dust model, especially in high-contrast molecular regions. An increase of the dust opacity at 353 GHz, tau_353/N_H, from the diffuse to the denser interstellar medium (ISM) is reported. It is associated with a decrease in the observed dust temperature, T_obs, that could be due at least in part to the increased dust opacity. We also report an excess of dust emission at HI column densities lower than 10^20 cm^-2 that could be the signature of dust in the warm ionized medium. In the diffuse ISM at high Galactic latitude, we report an anti-correlation between tau_353/N_H and T_obs while the dust specific luminosity, i.e., the total dust emission integrated over frequency (the radiance) per hydrogen atom, stays about constant. The implication is that in the diffuse high-latitude ISM tau_353 is not as reliable a tracer of dust column density as we conclude it is in molecular clouds where the correlation of tau_353 with dust extinction estimated using colour excess measurements on stars is strong. To estimate Galactic E(B-V) in extragalactic fields at high latitude we develop a new method based on the thermal dust radiance, instead of the dust optical depth, calibrated to E(B-V) using reddening measurements of quasars deduced from Sloan Digital Sky Survey data.

We discuss the methods employed to photometrically calibrate the data acquired by the Low Frequency Instrument on Planck. Our calibration is based on a combination of the Orbital Dipole plus the Solar Dipole, caused respectively by the motion of the Planck spacecraft with respect to the Sun and by motion of the Solar System with respect to the CMB rest frame. The latter provides a signal of a few mK with the same spectrum as the CMB anisotropies and is visible throughout the mission. In this data release we rely on the characterization of the Solar Dipole as measured by WMAP. We also present preliminary results (at 44GHz only) on the study of the Orbital Dipole, which agree with the WMAP value of the Solar System speed within our uncertainties. We compute the calibration constant for each radiometer roughly once per hour, in order to keep track of changes in the detectors' gain. Since non-idealities in the optical response of the beams proved to be important, we implemented a fast convolution algorithm which considers the full beam response in estimating the signal generated by the dipole. Moreover, in order to further reduce the impact of residual systematics due to sidelobes, we estimated time variations in the calibration constant of the 30GHz radiometers (the ones with the largest sidelobes) using the signal of a reference load. We have estimated the calibration accuracy in two ways: we have run a set of simulations to assess the impact of statistical errors and systematic effects in the instrument and in the calibration procedure, and we have performed a number of consistency checks on the data and on the brightness temperature of Jupiter. Calibration errors for this data release are expected to be about 0.6% at 44 and 70 GHz, and 0.8% at 30 GHz. (Abriged.)

On the arcminute angular scales probed by Planck, the CMB anisotropies are gently perturbed by gravitational lensing. Here we present a detailed study of this effect, detecting lensing independently in the 100, 143, and 217GHz frequency bands with an overall significance of greater than 25sigma. We use the temperature-gradient correlations induced by lensing to reconstruct a (noisy) map of the CMB lensing potential, which provides an integrated measure of the mass distribution back to the CMB last-scattering surface. Our lensing potential map is significantly correlated with other tracers of mass, a fact which we demonstrate using several representative tracers of large-scale structure. We estimate the power spectrum of the lensing potential, finding generally good agreement with expectations from the best-fitting LCDM model for the Planck temperature power spectrum, showing that this measurement at z=1100 correctly predicts the properties of the lower-redshift, later-time structures which source the lensing potential. When combined with the temperature power spectrum, our measurement provides degeneracy-breaking power for parameter constraints; it improves CMB-alone constraints on curvature by a factor of two and also partly breaks the degeneracy between the amplitude of the primordial perturbation power spectrum and the optical depth to reionization, allowing a measurement of the optical depth to reionization which is independent of large-scale polarization data. Discarding scale information, our measurement corresponds to a 4% constraint on the amplitude of the lensing potential power spectrum, or a 2% constraint on the RMS amplitude of matter fluctuations at z~2.

We describe the processing of the 531 billion raw data samples from the High Frequency Instrument (hereafter HFI), which we performed to produce six temperature maps from the first 473 days of Planck-HFI survey data. These maps provide an accurate rendition of the sky emission at 100, 143, 217, 353, 545, and 857 GHz with an angular resolution ranging from 9.7 to 4.6 arcmin. The detector noise per (effective) beam solid angle is respectively, 10, 6, 12 and 39 microKelvin in HFI four lowest frequency channel (100--353 GHz) and 13 and 14 kJy/sr for the 545 and 857 GHz channels. Using the 143 GHz channel as a reference, these two high frequency channels are intercalibrated within 5% and the 353 GHz relative calibration is at the percent level. The 100 and 217 GHz channels, which together with the 143 GHz channel determine the high-multipole part of the CMB power spectrum (50 < l <2500), are intercalibrated at better than 0.2 %.

The ESA's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14 May 2009 and has been scanning the microwave and submillimetre sky continuously since 12 August 2009. This paper gives an overview of the mission and its performance, the processing, analysis, and characteristics of the data, the scientific results, and the science data products and papers in the release. The science products include maps of the CMB and diffuse extragalactic foregrounds, a catalogue of compact Galactic and extragalactic sources, and a list of sources detected through the SZ effect. The likelihood code used to assess cosmological models against the Planck data and a lensing likelihood are described. Scientific results include robust support for the standard six-parameter LCDM model of cosmology and improved measurements of its parameters, including a highly significant deviation from scale invariance of the primordial power spectrum. The Planck values for these parameters and others derived from them are significantly different from those previously determined. Several large-scale anomalies in the temperature distribution of the CMB, first detected by WMAP, are confirmed with higher confidence. Planck sets new limits on the number and mass of neutrinos, and has measured gravitational lensing of CMB anisotropies at greater than 25 sigma. Planck finds no evidence for non-Gaussianity in the CMB. Planck's results agree well with results from the measurements of baryon acoustic oscillations. Planck finds a lower Hubble constant than found in some more local measures. Some tension is also present between the amplitude of matter fluctuations derived from CMB data and that derived from SZ data. The Planck and WMAP power spectra are offset from each other by an average level of about 2% around the first acoustic peak.

Anomalous microwave emission (AME) is believed to be due to electric dipole radiation from small spinning dust grains. The aim of this paper is a statistical study of the basic properties of AME regions and the environment in which they emit. We used WMAP and Planck maps, combined with ancillary radio and IR data, to construct a sample of 98 candidate AME sources, assembling SEDs for each source using aperture photometry on 1deg-smoothed maps from 0.408 GHz up to 3000 GHz. Each spectrum is fitted with a simple model of free-free, synchrotron (where necessary), cosmic microwave background (CMB), thermal dust, and spinning dust components. We find that 42 of the 98 sources have significant (>5sigma) excess emission at frequencies between 20 and 60 GHz. An analysis of the potential contribution of optically thick free-free emission from ultra-compact HII regions, using IR colour criteria, reduces the significant AME sample to 27 regions. The spectrum of the AME is consistent with model spectra of spinning dust. The AME regions tend to be more spatially extended than regions with little or no AME. The AME intensity is strongly correlated with the submillimetre/IR flux densities and comparable to previous AME detections in the literature. AME emissivity, defined as the ratio of AME to dust optical depth, varies by an order of magnitude for the AME regions. The AME regions tend to be associated with cooler dust in the range 14-20 K and an average emissivity index of +1.8, while the non-AME regions are typically warmer, at 20-27 K. In agreement with previous studies, the AME emissivity appears to decrease with increasing column density. The emerging picture is that the bulk of the AME is coming from the polycyclic aromatic hydrocarbons and small dust grains from the colder neutral interstellar medium phase (Abridged).

Models of $f(R)$ gravity that introduce corrections to the Newtonian potential in the weak field limit are tested at the scale of galaxy clusters. These models can explain the dynamics of spiral and elliptical galaxies without resorting to dark matter. We compute the pressure profiles of 579 galaxy clusters assuming that the gas is in hydrostatic equilibrium within the potential well of the modified gravitational field. The predicted profiles are compared with the average profile obtained by stacking the data of our cluster sample in the Planck foreground clean map SMICA. We find that the resulting profiles of these systems fit the data without requiring a dominant dark matter component, with model parameters similar to those required to explain the dynamics of galaxies. Our results do not rule out that clusters are dynamically dominated by Dark Matter but support the idea that Extended Theories of Gravity could provide an explanation to the dynamics of self-gravitating systems and to the present period of accelerated expansion, alternative to the concordance cosmological model.

We use Planck HFI data combined with ancillary radio data to study the emissivity index of the interstellar dust emission in the frequency range 10 - 353 GHz, or 3 - 0.8 mm, in the Galactic plane. We analyse the region l=20 degr - 44 degr and |b| \leq 4 degr where the free-free emission can be estimated from radio recombination line data. We fit the spectra at each sky pixel with a modified blackbody model and two spectral indices, beta_mm and beta_FIR, below and above 353 GHz respectively. We find that beta_mm is smaller than beta_FIR and we detect a correlation between this low frequency power-law index and the dust optical depth at 353 GHz, tau_353. The opacity spectral index beta_mm increases from about 1.54 in the more diffuse regions of the Galactic disk, |b| = 3 degr - 4 degr and tau_353 ~ 5 x 10^{-5}, to about 1.66 in the densest regions with an optical depth of more than one order of magnitude higher. We associate this correlation with an evolution of the dust emissivity related to the fraction of molecular gas along the line of sight. This translates into beta_mm ~ 1.54 for a medium that is mostly atomic and beta_mm ~ 1.66 when the medium is dominated by molecular gas. We find that both the Two-Level System model and the emission by ferromagnetic particles can explain the results. The results improve our understanding of the physics of interstellar dust and lead towards a complete model of the dust spectrum of the Milky Way from far-infrared to millimetre wavelengths.