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.

The two fundamental assumptions of the standard cosmological model - that the initial fluctuations are statistically isotropic and Gaussian - are rigorously tested using maps of the cosmic microwave background (CMB) anisotropy from the Planck satellite. Deviations from isotropy have been found and demonstrated to be robust against component separation algorithm, mask choice and frequency dependence. Many of these anomalies were previously observed in the WMAP data, and are now confirmed at similar levels of significance (about 3 sigma). However, we find little evidence for non-Gaussianity, with the exception of a few statistical signatures that seem to be associated with specific anomalies. In particular, we find that the quadrupole-octopole alignment is also connected to a low observed variance of the CMB signal. A power asymmetry is now found to persist to scales corresponding to about l=600, and can be described in the low-l regime by a phenomenological dipole modulation model. However, any primordial power asymmetry is strongly scale-dependent and does not extend to arbitrarily small angular scales. Finally, it is plausible that some of these features may be reflected in the angular power spectrum of the data, which shows a deficit of power on similar scales. Indeed, when the power spectra of two hemispheres defined by a preferred direction are considered separately, one shows evidence for a deficit in power, while its opposite contains oscillations between odd and even modes that may be related to the parity violation and phase correlations also detected in the data. Although these analyses represent a step forward in building an understanding of the anomalies, a satisfactory explanation based on physically motivated models is still lacking.

We establish for the first time heuristic correlations between harmonic space phase information and higher order statistics. Using the spherical full-sky maps of the cosmic microwave background as an example we demonstrate that known phase correlations at large spatial scales can gradually be diminished when subtracting a suitable best-fit (Bianchi-) template map of given strength. The weaker phase correlations lead in turn to a vanishing signature of anisotropy when measuring the Minkowski functionals and scaling indices in real-space and comparing them with surrogate maps being free of phase correlations. Those investigations can open a new road to a better understanding of signatures of non-Gaussianities in complex spatial structures by elucidating the meaning of Fourier phase correlations and their influence on higher order statistics.

Planck data have been used to provide stringent new constraints on cosmic strings and other defects. We describe forecasts of the CMB power spectrum induced by cosmic strings, calculating these from network models and simulations using line-of-sight Boltzmann solvers. We have studied Nambu-Goto cosmic strings, as well as field theory strings for which radiative effects are important, thus spanning the range of theoretical uncertainty in strings models. We have added the angular power spectrum from strings to that for a simple adiabatic model, with the extra fraction defined as $f_{10}$ at multipole $\ell=10$. This parameter has been added to the standard six parameter fit using COSMOMC with flat priors. For the Nambu-Goto string model, we have obtained a constraint on the string tension of $G\mu/c^2 < 1.5 x 10^{-7}$ and $f_{10} < 0.015$ at 95% confidence that can be improved to $G\mu/c^2 < 1.3 x 10^{-7}$ and $f_{10} < 0.010$ on inclusion of high-$\ell$ CMB data. For the abelian-Higgs field theory model we find, $G\mu_{AH}/c^2 < 3.2 x 10^{-7}$ and $f_{10} < 0.028$. The marginalized likelihoods for $f_{10}$ and in the $f_{10}$--$\Omega_b h^2$ plane are also presented. We have also obtained constraints on $f_{10}$ for models with semi-local strings and global textures for which $G\mu/c^2 < 1.1 x 10^{-6}$. We have made complementarity searches for the specific non-Gaussian signatures of cosmic strings, calibrating with all-sky Planck resolution CMB maps generated from networks of post-recombination strings. We have obtained upper limits on the string tension at 95% confidence of $G\mu/c^2 < 8.8 x 10^{-7}$ using modal bispectrum estimation and $G\mu/c^2 < 7.8 x 10^{-7}$ for real space searches with Minkowski functionals. These are conservative upper bounds because only post-recombination string contributions have been included in the non-Gaussian analysis.

In order to test for non-Gaussianities with respect to scale-dependencies we use so-called surrogate maps, in which possible phase correlations of the Fourier phases of the original WMAP data and simulations, respectively, are destroyed by applying a shuffling scheme to the maps. A statistical comparison of the original maps with the surrogate maps then allows to test for the existence of higher order correlations in the original maps, also and especially on well-defined Fourier modes. Using Minkowski functionals as an image analysis technique we calculate the deviation between the original data and 500 surrogates for different hemispheres in the sky and find ecliptic hemispherical asymmetries between northern and southern ecliptic sky. We find strong deviations from Gaussianity in the WMAP data when considering the low-l range with l = [2,20]. The analysis technique of the scaling indices leads to a slightly lower deviation. Although the underlying foreground reduction methods of the maps differ from each other, we find similar results for the WMAP 7yr ILC map and the WMAP 7yr (needlet-based) NILC map in the low-l range. Our results point once more to a cosmological nature of the signal. For a higher l range with l = [120,300] the results differ between the two image analysis techniques and between the two maps which makes an intrinsic nature of the signal on this l range less likely. When we decrease the size of the analysed sky regions for the low-l study, we do not find signatures of NG in the northern sky. In the south we find individual spots which show deviations from Gaussianity. In addition, we investigate non-Gaussian CMB simulations that depend on the f_NL parameter of the local type. These simulations cannot account for the detected signatures on the low-l range.

We demonstrate the feasibility to generate surrogates by Fourier-based methods for an incomplete data set. This is performed for the case of a CMB analysis, where astrophysical foreground emission, mainly present in the Galactic plane, is a major challenge. The shuffling of the Fourier phases for generating surrogates is now enabled by transforming the spherical harmonics into a new set of basis functions that are orthonormal on the cut sky. The results show that non-Gaussianities and hemispherical asymmetries in the CMB as identified in several former investigations, can still be detected even when the complete Galactic plane (|b| < 30{\deg}) is removed. We conclude that the Galactic plane cannot be the dominant source for these anomalies. The results point towards a violation of statistical isotropy.