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.

We consider cosmologies in which a dark-energy scalar field interacts with cold dark matter. The growth of perturbations is followed beyond the linear level by means of the time-renormalization-group method, which is extended to describe a multi-component matter sector. Even in the absence of the extra interaction, a scale-dependent bias is generated as a consequence of the different initial conditions for baryons and dark matter after decoupling. The effect is enhanced significantly by the extra coupling and can be at the 2-3 percent level in the range of scales of baryonic acoustic oscillations. We compare our results with N-body simulations, finding very good agreement.

We describe the data processing pipeline of the Planck Low Frequency Instrument (LFI) data processing centre (DPC) to create and characterize full-sky maps based on the first 15.5 months of operations at 30, 44 and 70 GHz. In particular, we discuss the various steps involved in reducing the data, starting from telemetry packets through to the production of cleaned, calibrated timelines and calibrated frequency maps. Data are continuously calibrated using the modulation induced on the mean temperature of the cosmic microwave background radiation by the proper motion of the spacecraft. Sky signals other than the dipole are removed by an iterative procedure based on simultaneous fitting of calibration parameters and sky maps. Noise properties are estimated from time-ordered data after the sky signal has been removed, using a generalized least square map-making algorithm. A destriping code (Madam) is employed to combine radiometric data and pointing information into sky maps, minimizing the variance of correlated noise. Noise covariance matrices, required to compute statistical uncertainties on LFI and Planck products, are also produced. Main beams are estimated down to the -20 dB level using Jupiter transits, which are also used for the geometrical calibration of the focal plane.

We aim to summarize the current state of knowledge regarding Galactic Faraday rotation in an all-sky map of the Galactic Faraday depth. For this we have assembled the most extensive catalog of Faraday rotation data of compact extragalactic polarized radio sources to date. In the map making procedure we use a recently developed algorithm that reconstructs the map and the power spectrum of a statistically isotropic and homogeneous field while taking into account uncertainties in the noise statistics. This procedure is able to identify some rotation angles that are offset by an integer multiple of pi. The resulting map can be seen as an improved version of earlier such maps and is made publicly available, along with a map of its uncertainty. For the angular power spectrum we find a power law behavior with a power law index of -2.14 for a Faraday sky where an overall variance profile as a function of Galactic latitude has been removed, in agreement with earlier work. We show that this is in accordance with a 3D Fourier power spectrum P(k) proportional to k^-2.14 of the underlying field n_e times B_r under simplifying geometrical and statistical assumptions.

The scientific performance of the Planck Low Frequency Instrument (LFI) after one year of in-orbit operation is presented. We describe the main optical parameters and discuss photometric calibration, white noise sensitivity, and noise properties. A preliminary evaluation of the impact of the main systematic effects is presented. For each of the performance parameters, we outline the methods used to obtain them from the flight data and provide a comparison with pre-launch ground assessments, which are essentially confirmed in flight.

We describe the processing of data from the Low Frequency Instrument (LFI) used in production of the Planck Early Release Compact Source Catalogue (ERCSC). In particular, we discuss the steps involved in reducing the data from telemetry packets to cleaned, calibrated, time-ordered data (TOD) and frequency maps. Data are continuously calibrated using the modulation of the temperature of the cosmic microwave background radiation induced by the motion of the spacecraft. Noise properties are estimated from TOD from which the sky signal has been removed using a generalized least square map-making algorithm. Measured 1/f noise knee-frequencies range from 100mHz at 30GHz to a few tens of mHz at 70GHz. A destriping code (Madam) is employed to combine radiometric data and pointing information into sky maps, minimizing the variance of correlated noise. Noise covariance matrices required to compute statistical uncertainties on LFI and Planck products are also produced. Main beams are estimated down to the approx -10dB level using Jupiter transits, which are also used for geometrical calibration of the focal plane.

One of the deepest and most long-standing mysteries in physics has been the huge discrepancy between the observed vacuum density and our expectations from theories of high energy physics, which has been dubbed the Old Cosmological Constant problem. One proposal to address this puzzle at the semi-classical level is to decouple quantum vacuum from space-time geometry via a modification of gravity that includes an incompressible fluid, known as Gravitational Aether. In this paper, we discuss classical predictions of this theory along with its compatibility with cosmological and experimental tests of gravity. We argue that deviations from General Relativity (GR) in this theory are sourced by pressure or vorticity. In particular, the theory predicts that the gravitational constant for radiation is 33% larger than that of non-relativistic matter, which is preferred by (most) cosmic microwave background (CMB), Lyman-Alpha forest, and Lithium-7 primordial abundance observations, while being consistent with other cosmological tests at ~2-sigma level. It is further shown that all Parametrized Post-Newtonian (PPN) parameters have the standard GR values aside from the anomalous coupling to pressure, which has not been directly measured. A more subtle prediction of this model (assuming irrotational aether) is that the (intrinsic) gravitomagnetic effect is 33% larger than GR prediction. This is consistent with current limits from LAGEOS and Gravity Probe B at ~2-sigma level.

We derive a method to reconstruct Gaussian signals from linear measurements with Gaussian noise. This new algorithm is intended for applications in astrophysics and other sciences. The starting point of our considerations is the principle of minimum Gibbs free energy which was previously used to derive a signal reconstruction algorithm handling uncertainties in the signal covariance. We extend this algorithm to simultaneously uncertain noise and signal covariances using the same principles in the derivation. The resulting equations are general enough to be applied in many different contexts. We demonstrate the performance of the algorithm by applying it to specific example situations and compare it to algorithms not allowing for uncertainties in the noise covariance. The results show that the method we suggest performs very well under a variety of circumstances and is indeed qualitatively superior to the other methods in cases where uncertainty in the noise covariance is present.

We present a first application of the recently proposed LITMUS test for magnetic helicity, as well as a thorough study of its applicability under different circumstances. In order to apply this test to the galactic magnetic field, the newly developed critical filter formalism is used to produce an all-sky map of the Faraday depth. The test does not detect helicity in the galactic magnetic field. To understand the significance of this finding, we made an applicability study, showing that a definite conclusion about the absence of magnetic helicity in the galactic field has not yet been reached. This study is conducted by applying the test to simulated observational data. We consider simulations in a flat sky approximation and all-sky simulations, both with assumptions of constant electron densities and realistic distributions of thermal and cosmic ray electrons. Our results suggest that the LITMUS test does indeed perform very well in cases where constant electron densities can be assumed, both in the flat-sky limit and in the galactic setting. Non-trivial distributions of thermal and cosmic ray electrons, however, may complicate the scenario to the point where helicity in the magnetic field can escape detection.

If the accelerated expansion of the Universe at the present epoch is driven by a dark energy scalar field, there may well be a non-trivial coupling between the dark energy and the cold dark matter (CDM) fluid. Such interactions give rise to new features in cosmological structure growth, like an additional long-range attractive force between CDM particles, or variations of the dark matter particle mass with time. We have implemented these effects in the N-body code GADGET-2 and present results of a series of high-resolution N-body simulations where the dark energy component is directly interacting with the cold dark matter. As a consequence of the new physics, CDM and baryon distributions evolve differently both in the linear and in the nonlinear regime of structure formation. Already on large scales a linear bias develops between these two components, which is further enhanced by the nonlinear evolution. We also find, in contrast with previous work, that the density profiles of CDM halos are less concentrated in coupled dark energy cosmologies compared with LCDM, and that this feature does not depend on the initial conditions setup, but is a specific consequence of the extra physics induced by the coupling. Also, the baryon fraction in halos in the coupled models is significantly reduced below the universal baryon fraction. These features alleviate tensions between observations and the LCDM model on small scales. Our methodology is ideally suited to explore the predictions of coupled dark energy models in the fully non-linear regime, which can provide powerful constraints for the viable parameter space of such scenarios.

Under the assumption that the variations of parameters of nature and the current acceleration of the universe are related and governed by the evolution of a single scalar field, we show how information can be obtained on the nature of dark energy from observational detection of (or constraints on) cosmological variations of the fine structure constant and the proton-to-electron mass ratio. We also comment on the current observational status, and on the prospects for improvements with future spectrographs such as ESPRESSO and CODEX.

The cosmological backreaction arises when one directly averages the Einstein equations to recover an effective Robertson-Walker cosmology, rather than assuming a background a priori. While usually discussed in the context of dark energy, strictly speaking any cosmological model should be recovered from such a procedure. We apply the Buchert averaging formalism to linear Robertson-Walker universes containing matter, radiation and dark energy and evaluate numerically the discrepancies between the assumed and the averaged behaviour, finding the largest deviations for an Einstein-de Sitter universe, increasing rapidly with Hubble rate to a 0.01% effect for h=0.701. For the LCDM concordance model, the backreaction is of the order of Omega_eff~4x10^-6, with those for dark energy models being within a factor of two or three. The impacts at recombination are of the order of 10^-8 and those in deep radiation domination asymptote to a constant value. While the effective equations of state of the backreactions in Einstein-de Sitter, concordance and quintessence models are generally dust-like, a backreaction with an equation of state w_eff<-1/3 can be found for strongly phantom models.

We show that growing neutrino models, in which the growing neutrino mass stops the dynamical evolution of a dark energy scalar field, lead to a substantial neutrino clustering on the scales of superclusters. Nonlinear neutrino lumps form at redshift z ~ 1 and could partially drag the clustering of dark matter. If observed, large scale non-linear structures could be an indication for a new attractive force stronger than gravity and mediated by the 'cosmon' dark energy scalar field.

The presence of dark energy at high redshift influences both the cosmic sound horizon and the distance to last scattering of the cosmic microwave background. We demonstrate that through the degeneracy in their ratio, early dark energy can lie hidden in the CMB temperature and polarization spectra, leading to an unrecognized shift in the sound horizon. If the sound horizon is then used as a standard ruler, as in baryon acoustic oscillations, then the derived cosmological parameters can be nontrivially biased. Fitting for the absolute ruler scale (just as supernovae must be fit for the absolute candle magnitude) removes the bias but decreases the leverage of the BAO technique by a factor 2.

A growing neutrino mass can stop the dynamical evolution of a dark energy scalar field, thus explaining the 'why now' problem. We show that such models lead to a substantial neutrino clustering on the scales of superclusters. Nonlinear neutrino lumps form at redshift z \sim 1 and could partially drag the clustering of dark matter. If observed, large scale non-linear structures could be an indication for a new attractive 'cosmon force' stronger than gravity.

We reformulate the averaged Einstein equations in a form suitable for use with Newtonian gauge linear perturbation theory and track the size of the modifications to standard Robertson-Walker evolution on the largest scales as a function of redshift for both Einstein de-Sitter and Lambda CDM cosmologies. In both cases the effective energy density arising from linear perturbations is of the order of 10^-5 the matter density, as would be expected, with an effective equation of state w ~ -1/19. Employing a modified Halofit code to extend our results to quasilinear scales, we find that, while larger, the deviations from Robertson-Walker behaviour remain of the order of 10^-5.

Einstein's theory of general relativity, which contains a universal value of the Planck mass, has been so far successfully invoked to explain gravitational dynamics from sub-millimeter scales to the scale of the cosmological horizon. However, one may envisage that in alternative theories of gravity, the effective value of the Planck mass (or Newton's constant), which quantifies the coupling of matter to metric perturbations, can run on the cosmological horizon scale. In this letter, we study the consequences of a glitch in the Planck mass from sub-horizon to super-horizon scales. We first give three examples of models that naturally exhibit this feature, and then show that current cosmological observations severely constrain this glitch to less than 1.2%. This is the strongest constraint to date, on natural (i.e. non-fine-tuned) deviations from Einstein gravity on the cosmological horizon scale.

The first three years of observation of the Wilkinson Microwave Anisotropy Probe (WMAP) have provided the most precise data on the anisotropies of the cosmic microwave background (CMB) to date. We investigate the impact of these results and their combination with data from other astrophysical probes on cosmological models with a dynamical dark energy component. By considering a wide range of such models, we find that the constraints on dynamical dark energy are significantly improved compared to the first year data.

We propose a novel parameterization of the dark energy density. It is particularly well suited to describe a non-negligible contribution of dark energy at early times and contains only three parameters, which are all physically meaningful: the fractional dark energy density today, the equation of state today and the fractional dark energy density at early times. As we parameterize Omega_d(a) directly instead of the equation of state, we can give analytic expressions for the Hubble parameter, the conformal horizon today and at last scattering, the sound horizon at last scattering, the acoustic scale as well as the luminosity distance. For an equation of state today w_0 < -1, our model crosses the cosmological constant boundary. We perform numerical studies to constrain the parameters of our model by using Cosmic Microwave Background, Large Scale Structure and Supernovae Ia data. At 95% confidence, we find that the fractional dark energy density at early times Omega_early < 0.06. This bound tightens considerably to Omega_early < 0.04 when the latest Boomerang data is included. We find that both the gold sample of Riess et. al. and the SNLS data by Astier et. al. when combined with CMB and LSS data mildly prefer w_0 < -1, but are well compatible with a cosmological constant.