We introduce a new mathematical tool (a direction-dependent probe) to analyse the randomness of purported isotropic Gaussian random fields on the sphere. We apply the probe to assess the full-sky cosmic microwave background (CMB) temperature maps produced by the {\it Planck} collaboration (PR2 2015 and PR3 2018), with special attention to the inpainted maps. To study the randomness of the fields represented by each map we use the autocorrelation of the sequence of probe coefficients (which are just the full-sky Fourier coefficients $a_{\ell,0}$ if the $z$ axis is taken in the probe direction). If the field is {isotropic and Gaussian} then the probe coefficients for a given direction should be realisations of uncorrelated scalar Gaussian random variables. We introduce a particular function on the sphere (called the \emph{AC discrepancy}) that accentuates the departure from Gaussianity and isotropy. We find that for some of the maps, there are many directions for which the departures are significant, especially near the galactic plane. We also study the effect of varying the highest multipole used to calculate the AC discrepancy from the initial value of $1500$ to $2500$. In the case of Commander 2015, the AC discrepancy now exhibits antipodal "blobs" well away from the galactic plane. Finally, we look briefly at the non-inpainted Planck maps, for which the computed AC discrepancy maps have a very different character, with features that are global rather than local. For the particular case of the non-inpainted 2018 \texttt{SEVEM} map (which has visible equatorial pollution), we model with partial success the observed behaviour by an isotropic Gaussian random field added to a non-random needlet-like structure located near the galactic centre.

Cosmic microwave background radiation (CMB) is critical to the understanding of the early universe and precise estimation of cosmological constants. Due to the contamination of thermal dust noise in the galaxy, the CMB map that is an image on the two-dimensional sphere has missing observations, mainly concentrated on the equatorial region. The noise of the CMB map has a significant impact on the estimation precision for cosmological parameters. Inpainting the CMB map can effectively reduce the uncertainty of parametric estimation. In this paper, we propose a deep learning-based variational autoencoder --- CosmoVAE, to restoring the missing observations of the CMB map. The input and output of CosmoVAE are square images. To generate training, validation, and test data sets, we segment the full-sky CMB map into many small images by Cartesian projection. CosmoVAE assigns physical quantities to the parameters of the VAE network by using the angular power spectrum of the Gaussian random field as latent variables. CosmoVAE adopts a new loss function to improve the learning performance of the model, which consists of $\ell_1$ reconstruction loss, Kullback-Leibler divergence between the posterior distribution of encoder network and the prior distribution of latent variables, perceptual loss, and total-variation regularizer. The proposed model achieves state of the art performance for Planck \texttt{Commander} 2018 CMB map inpainting.

The European Space Agency's Planck satellite, which was dedicated to studying the early Universe and its subsequent evolution, was launched on 14 May 2009. It scanned the microwave and submillimetre sky continuously between 12 August 2009 and 23 October 2013, producing deep, high-resolution, all-sky maps in nine frequency bands from 30 to 857GHz. This paper presents the cosmological legacy of Planck, which currently provides our strongest constraints on the parameters of the standard cosmological model and some of the tightest limits available on deviations from that model. The 6-parameter LCDM model continues to provide an excellent fit to the cosmic microwave background data at high and low redshift, describing the cosmological information in over a billion map pixels with just six parameters. With 18 peaks in the temperature and polarization angular power spectra constrained well, Planck measures five of the six parameters to better than 1% (simultaneously), with the best-determined parameter (theta_*) now known to 0.03%. We describe the multi-component sky as seen by Planck, the success of the LCDM model, and the connection to lower-redshift probes of structure formation. We also give a comprehensive summary of the major changes introduced in this 2018 release. The Planck data, alone and in combination with other probes, provide stringent constraints on our models of the early Universe and the large-scale structure within which all astrophysical objects form and evolve. We discuss some lessons learned from the Planck mission, and highlight areas ripe for further experimental advances.

We present cosmological parameter results from the final full-mission Planck measurements of the CMB anisotropies. We find good consistency with the standard spatially-flat 6-parameter $\Lambda$CDM cosmology having a power-law spectrum of adiabatic scalar perturbations (denoted "base $\Lambda$CDM" in this paper), from polarization, temperature, and lensing, separately and in combination. A combined analysis gives dark matter density $\Omega_c h^2 = 0.120\pm 0.001$, baryon density $\Omega_b h^2 = 0.0224\pm 0.0001$, scalar spectral index $n_s = 0.965\pm 0.004$, and optical depth $\tau = 0.054\pm 0.007$ (in this abstract we quote $68\,\%$ confidence regions on measured parameters and $95\,\%$ on upper limits). The angular acoustic scale is measured to $0.03\,\%$ precision, with $100\theta_*=1.0411\pm 0.0003$. These results are only weakly dependent on the cosmological model and remain stable, with somewhat increased errors, in many commonly considered extensions. Assuming the base-$\Lambda$CDM cosmology, the inferred late-Universe parameters are: Hubble constant $H_0 = (67.4\pm 0.5)$km/s/Mpc; matter density parameter $\Omega_m = 0.315\pm 0.007$; and matter fluctuation amplitude $\sigma_8 = 0.811\pm 0.006$. We find no compelling evidence for extensions to the base-$\Lambda$CDM model. Combining with BAO we constrain the effective extra relativistic degrees of freedom to be $N_{\rm eff} = 2.99\pm 0.17$, and the neutrino mass is tightly constrained to $\sum m_\nu< 0.12$eV. The CMB spectra continue to prefer higher lensing amplitudes than predicted in base -$\Lambda$CDM at over $2\,\sigma$, which pulls some parameters that affect the lensing amplitude away from the base-$\Lambda$CDM model; however, this is not supported by the lensing reconstruction or (in models that also change the background geometry) BAO data. (Abridged)

We report on the implications for cosmic inflation of the 2018 Release of the Planck CMB anisotropy measurements. The results are fully consistent with the two previous Planck cosmological releases, but have smaller uncertainties thanks to improvements in the characterization of polarization at low and high multipoles. Planck temperature, polarization, and lensing data determine the spectral index of scalar perturbations to be $n_\mathrm{s}=0.9649\pm 0.0042$ at 68% CL and show no evidence for a scale dependence of $n_\mathrm{s}.$ Spatial flatness is confirmed at a precision of 0.4% at 95% CL with the combination with BAO data. The Planck 95% CL upper limit on the tensor-to-scalar ratio, $r_{0.002}<0.10$, is further tightened by combining with the BICEP2/Keck Array BK15 data to obtain $r_{0.002}<0.056$. In the framework of single-field inflationary models with Einstein gravity, these results imply that: (a) slow-roll models with a concave potential, $V" (\phi) < 0,$ are increasingly favoured by the data; and (b) two different methods for reconstructing the inflaton potential find no evidence for dynamics beyond slow roll. Non-parametric reconstructions of the primordial power spectrum consistently confirm a pure power law. A complementary analysis also finds no evidence for theoretically motivated parameterized features in the Planck power spectrum, a result further strengthened for certain oscillatory models by a new combined analysis that includes Planck bispectrum data. The new Planck polarization data provide a stringent test of the adiabaticity of the initial conditions. The polarization data also provide improved constraints on inflationary models that predict a small statistically anisotropic quadrupolar modulation of the primordial fluctuations. However, the polarization data do not confirm physical models for a scale-dependent dipolar modulation.

This paper describes the 2018 Planck CMB likelihoods, following a hybrid approach similar to the 2015 one, with different approximations at low and high multipoles, and implementing several methodological and analysis refinements. With more realistic simulations, and better correction and modelling of systematics, we can now make full use of the High Frequency Instrument polarization data. The low-multipole 100x143 GHz EE cross-spectrum constrains the reionization optical-depth parameter $\tau$ to better than 15% (in combination with with the other low- and high-$\ell$ likelihoods). We also update the 2015 baseline low-$\ell$ joint TEB likelihood based on the Low Frequency Instrument data, which provides a weaker $\tau$ constraint. At high multipoles, a better model of the temperature-to-polarization leakage and corrections for the effective calibrations of the polarization channels (polarization efficiency or PE) allow us to fully use the polarization spectra, improving the constraints on the $\Lambda$CDM parameters by 20 to 30% compared to TT-only constraints. Tests on the modelling of the polarization demonstrate good consistency, with some residual modelling uncertainties, the accuracy of the PE modelling being the main limitation. Using our various tests, simulations, and comparison between different high-$\ell$ implementations, we estimate the consistency of the results to be better than the 0.5$\sigma$ level. Minor curiosities already present before (differences between $\ell$<800 and $\ell$>800 parameters or the preference for more smoothing of the $C_\ell$ peaks) are shown to be driven by the TT power spectrum and are not significantly modified by the inclusion of polarization. Overall, the legacy Planck CMB likelihoods provide a robust tool for constraining the cosmological model and represent a reference for future CMB observations. (Abridged)

We present measurements of the cosmic microwave background (CMB) lensing potential using the final $\textit{Planck}$ 2018 temperature and polarization data. We increase the significance of the detection of lensing in the polarization maps from $5\,\sigma$ to $9\,\sigma$. Combined with temperature, lensing is detected at $40\,\sigma$. We present an extensive set of tests of the robustness of the lensing-potential power spectrum, and construct a minimum-variance estimator likelihood over lensing multipoles $8 \le L \le 400$. We find good consistency between lensing constraints and the results from the $\textit{Planck}$ CMB power spectra within the $\rm{\Lambda CDM}$ model. Combined with baryon density and other weak priors, the lensing analysis alone constrains $\sigma_8 \Omega_{\rm m}^{0.25}=0.589\pm 0.020$ ($1\,\sigma$ errors). Also combining with baryon acoustic oscillation (BAO) data, we find tight individual parameter constraints, $\sigma_8=0.811\pm0.019$, $H_0=67.9_{-1.3}^{+1.2}\,\text{km}\,\text{s}^{-1}\,\rm{Mpc}^{-1}$, and $\Omega_{\rm m}=0.303^{+0.016}_{-0.018}$. Combining with $\textit{Planck}$ CMB power spectrum data, we measure $\sigma_8$ to better than $1\,\%$ precision, finding $\sigma_8=0.811\pm 0.006$. We find consistency with the lensing results from the Dark Energy Survey, and give combined lensing-only parameter constraints that are tighter than joint results using galaxy clustering. Using $\textit{Planck}$ cosmic infrared background (CIB) maps we make a combined estimate of the lensing potential over $60\,\%$ of the sky with considerably more small-scale signal. We demonstrate delensing of the $\textit{Planck}$ power spectra, detecting a maximum removal of $40\,\%$ of the lensing-induced power in all spectra. The improvement in the sharpening of the acoustic peaks by including both CIB and the quadratic lensing reconstruction is detected at high significance (abridged).

We examine the effects ultra-light primordial black holes (PBHs) have on the anisotropies of the cosmic microwave background (CMB). PBHs in the mass range of $10^{15}$ to $10^{17}$ g emit Hawking radiation in the early Universe, modifying the standard recombination history. This leads to a damping of small-scale temperature and polarisation anisotropies and enhances large-scale polarisation fluctuations. As some models of inflation predict PBHs with a range of masses, we investigate the impacts of extended mass distributions on PBH abundance constraints. We model PBH energy injection using a ground-up approach incorporating species-dependent deposition efficiencies. By allowing the $\Lambda$CDM parameters to vary simultaneously with the PBH fraction and mass, we show that exclusion bounds on the PBH fraction of DM $f_\text{PBH}$ are relaxed by up to an order of magnitude, compared to the case of fixed $\Lambda$CDM parameters. We also give 95% exclusion regions for $f_\text{PBH}$ for a variety of mass distributions. In particular, for a uniform mass distribution between $10^{15}$ and $10^{17}$ g, we find $f_\text{PBH} < 1.6 \times 10^{-5}$ when allowing $\Lambda$CDM parameters to vary.

Analysis of the Planck 2018 data set indicates that the statistical properties of the cosmic microwave background (CMB) temperature anisotropies are in excellent agreement with previous studies using the 2013 and 2015 data releases. In particular, they are consistent with the Gaussian predictions of the $\Lambda$CDM cosmological model, yet also confirm the presence of several so-called "anomalies" on large angular scales. The novelty of the current study, however, lies in being a first attempt at a comprehensive analysis of the statistics of the polarization signal over all angular scales, using either maps of the Stokes parameters, $Q$ and $U$, or the $E$-mode signal derived from these using a new methodology (which we describe in an appendix). Although remarkable progress has been made in reducing the systematic effects that contaminated the 2015 polarization maps on large angular scales, it is still the case that residual systematics (and our ability to simulate them) can limit some tests of non-Gaussianity and isotropy. However, a detailed set of null tests applied to the maps indicates that these issues do not dominate the analysis on intermediate and large angular scales (i.e., $\ell \lesssim 400$). In this regime, no unambiguous detections of cosmological non-Gaussianity, or of anomalies corresponding to those seen in temperature, are claimed. Notably, the stacking of CMB polarization signals centred on the positions of temperature hot and cold spots exhibits excellent agreement with the $\Lambda$CDM cosmological model, and also gives a clear indication of how Planck provides state-of-the-art measurements of CMB temperature and polarization on degree scales.

We analyse the Planck full-mission cosmic microwave background (CMB) temperature and E-mode polarization maps to obtain constraints on primordial non-Gaussianity (NG). We compare estimates obtained from separable template-fitting, binned, and modal bispectrum estimators, finding consistent values for the local, equilateral, and orthogonal bispectrum amplitudes. Our combined temperature and polarization analysis produces the following results: f_NL^local = -0.9 +\- 5.1; f_NL^equil = -26 +\- 47; and f_NL^ortho = - 38 +\- 24 (68%CL, statistical). These results include the low-multipole (4 <= l < 40) polarization data, not included in our previous analysis, pass an extensive battery of tests, and are stable with respect to our 2015 measurements. Polarization bispectra display a significant improvement in robustness; they can now be used independently to set NG constraints. We consider a large number of additional cases, e.g. scale-dependent feature and resonance bispectra, isocurvature primordial NG, and parity-breaking models, where we also place tight constraints but do not detect any signal. The non-primordial lensing bispectrum is detected with an improved significance compared to 2015, excluding the null hypothesis at 3.5 sigma. We present model-independent reconstructions and analyses of the CMB bispectrum. Our final constraint on the local trispectrum shape is g_NLl^local = (-5.8 +\-6.5) x 10^4 (68%CL, statistical), while constraints for other trispectra are also determined. We constrain the parameter space of different early-Universe scenarios, including general single-field models of inflation, multi-field and axion field parity-breaking models. Our results provide a high-precision test for structure-formation scenarios, in complete agreement with the basic picture of the LambdaCDM cosmology regarding the statistics of the initial conditions (abridged).

Non-parametric reconstruction methods, such as Gaussian process (GP) regression, provide a model-independent way of estimating an underlying function and its uncertainty from noisy data. We demonstrate how GP-reconstruction can be used as a consistency test between a given data set and a specific model by looking for structures in the residuals of the data with respect to the model's best-fit. Applying this formalism to the Planck temperature and polarisation power spectrum measurements, we test their global consistency with the predictions of the base $\Lambda$CDM model. Our results do not show any serious inconsistencies, lending further support to the interpretation of the base $\Lambda$CDM model as cosmology's gold standard.

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.

We present a comprehensive review of keV-scale sterile neutrino Dark Matter, collecting views and insights from all disciplines involved - cosmology, astrophysics, nuclear, and particle physics - in each case viewed from both theoretical and experimental/observational perspectives. After reviewing the role of active neutrinos in particle physics, astrophysics, and cosmology, we focus on sterile neutrinos in the context of the Dark Matter puzzle. Here, we first review the physics motivation for sterile neutrino Dark Matter, based on challenges and tensions in purely cold Dark Matter scenarios. We then round out the discussion by critically summarizing all known constraints on sterile neutrino Dark Matter arising from astrophysical observations, laboratory experiments, and theoretical considerations. In this context, we provide a balanced discourse on the possibly positive signal from X-ray observations. Another focus of the paper concerns the construction of particle physics models, aiming to explain how sterile neutrinos of keV-scale masses could arise in concrete settings beyond the Standard Model of elementary particle physics. The paper ends with an extensive review of current and future astrophysical and laboratory searches, highlighting new ideas and their experimental challenges, as well as future perspectives for the discovery of sterile neutrinos.

Future observations of 21~cm emission using HI intensity mapping will enable us to probe the large scale structure of the Universe over very large survey volumes within a reasonable observation time. We demonstrate that the three-dimensional information contained in such surveys will be an extremely powerful tool in searching for features that were imprinted in the primordial power spectrum and bispectrum during inflation. Here we focus on the "resonant" and "step" inflation models, and forecast the potential of upcoming 21~cm experiments to detect these inflationary features in the observable power- and bispectrum. We find that the full scale Tianlai experiment and the Square Kilometre Array (SKA) have the potential to improve on the sensitivity of current Cosmic Microwave Background (CMB) experiments by several orders of magnitude.

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

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

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

The European Space Agency's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14~May 2009 and scanned the microwave and submillimetre sky continuously between 12~August 2009 and 23~October 2013. In February~2015, ESA and the Planck Collaboration released the second set of cosmology products based on data from the entire Planck mission, including both temperature and polarization, along with a set of scientific and technical papers and a web-based explanatory supplement. This paper gives an overview of the main characteristics of the data and the data products in the release, as well as the associated cosmological and astrophysical science results and papers. The science products include maps of the cosmic microwave background (CMB), the thermal Sunyaev-Zeldovich effect, and diffuse foregrounds in temperature and polarization, catalogues of compact Galactic and extragalactic sources (including separate catalogues of Sunyaev-Zeldovich clusters and Galactic cold clumps), and extensive simulations of signals and noise used in assessing the performance of the analysis methods and assessment of uncertainties. The likelihood code used to assess cosmological models against the Planck data are described, as well as a CMB lensing likelihood. Scientific results include cosmological parameters deriving from CMB power spectra, gravitational lensing, and cluster counts, as well as constraints on inflation, non-Gaussianity, primordial magnetic fields, dark energy, and modified gravity.

All cosmological observations to date are consistent with adiabatic, Gaussian and nearly scale invariant initial conditions. These findings provide strong evidence for a particular symmetry breaking pattern in the very early universe (with a close to vanishing order parameter, $\epsilon$), widely accepted as conforming to the predictions of the simplest realizations of the inflationary paradigm. However, given that our observations are only privy to perturbations, in inferring something about the background that gave rise to them, it should be clear that many different underlying constructions project onto the same set of cosmological observables. Features in the primordial correlation functions, if present, would offer a unique and discriminating window onto the parent theory in which the mechanism that generated the initial conditions is embedded. In certain contexts, simple linear response theory allows us to infer new characteristic scales from the presence of features that can break the aforementioned degeneracies among different background models, and in some cases can even offer a limited spectroscopy of the heavier degrees of freedom that couple to the inflaton. In this review, we offer a pedagogical survey of the diverse, theoretically well grounded mechanisms which can imprint features into primordial correlation functions in addition to reviewing the techniques one can employ to probe observations. These observations include cosmic microwave background anisotropies and spectral distortions as well as the matter two and three point functions as inferred from large-scale structure and potentially, 21 cm surveys.

We study the potential of a future, large-volume photometric survey to constrain the axion mass $m_a$ in the hot dark matter limit. Future surveys such as Euclid will have significantly more constraining power than current observations for hot dark matter. Nonetheless, the lowest accessible axion masses are limited by the fact that axions lighter than $\sim 0.15$ eV decouple before the QCD epoch, assumed here to occur at a temperature $T_{\rm QCD} \sim 170$ MeV; this leaves an axion population of such low density that its late-time cosmological impact is negligible. For larger axion masses, $m_a \gtrsim 0.15$ eV, where axions remain in equilibrium until after the QCD phase transition, we find that a Euclid-like survey combined with Planck CMB data can detect $m_a$ at very high significance. Our conclusions are robust against assumptions about prior knowledge of the neutrino mass. Given that the proposed IAXO solar axion search is sensitive to $m_a\lesssim 0.2$ eV, the axion mass range probed by cosmology is nicely complementary.

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 assess the potential of a future large-volume photometric redshift survey to constrain observational inflationary parameters using three large-scale structure observables: the angular shear and galaxy power spectra, and the cluster mass function measured through weak lensing. When used in combination with Planck-like CMB measurements, we find that the spectral index n_s can be constrained to a 1 sigma precision of up to 0.0025. The sensitivity to the running of the spectral index can potentially improve to 0.0017, roughly a factor of five better than the present 1 sigma~constraint from Planck and auxiliary CMB data, allowing us to test the assumptions of the slow-roll scenario with unprecedented accuracy. Interestingly, neither CMB+shear nor CMB+galaxy nor CMB+clusters alone can achieve this level of sensitivity; it is the combined power of all three probes that conspires to break the different parameter degeneracies inherent in each type of observations. We make our forecast software publicly available via download or upon request from the authors.

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 perform a detailed forecast on how well a {\sc Euclid}-like survey will be able to constrain dark energy and neutrino parameters from a combination of its cosmic shear power spectrum, galaxy power spectrum, and cluster mass function measurements. We find that the combination of these three probes vastly improves the survey's potential to measure the time evolution of dark energy. In terms of a dark energy figure-of-merit defined as $(\sigma(w_{\mathrm p}) \sigma(w_a))^{-1}$, we find a value of 690 for {\sc Euclid}-like data combined with {\sc Planck}-like measurements of the cosmic microwave background (CMB) anisotropies in a 10-dimensional cosmological parameter space, assuming a $\Lambda$CDM fiducial cosmology. For the more commonly used 7-parameter model, we find a figure-of-merit of 1900 for the same data combination. We consider also the survey's potential to measure dark energy perturbations in models wherein the dark energy is parameterised as a fluid with a nonstandard non-adiabatic sound speed, and find that in an \emph{optimistic} scenario in which $w_0$ deviates by as much as is currently observationally allowed from $-1$, models with $\hat{c}_\mathrm{s}^2 = 10^{-6}$ and $\hat{c}_\mathrm{s}^2 = 1$ can be distinguished at more than $2\sigma$ significance. We emphasise that constraints on the dark energy sound speed from cluster measurements are strongly dependent on the modelling of the cluster mass function; significantly weaker sensitivities ensue if we modify our model to include fewer features of nonlinear dark energy clustering. Finally, we find that the sum of neutrino masses can be measured with a $1 \sigma$ precision of 0.015~eV, (abridged)

We present the first results based on Planck measurements of the CMB temperature and lensing-potential power spectra. The Planck spectra at high multipoles are extremely well described by the standard spatially-flat six-parameter LCDM cosmology. In this model Planck data determine the cosmological parameters to high precision. We find a low value of the Hubble constant, H0=67.3+/-1.2 km/s/Mpc and a high value of the matter density parameter, Omega_m=0.315+/-0.017 (+/-1 sigma errors) in excellent agreement with constraints from baryon acoustic oscillation (BAO) surveys. Including curvature, we find that the Universe is consistent with spatial flatness to percent-level precision using Planck CMB data alone. We present results from an analysis of extensions to the standard cosmology, using astrophysical data sets in addition to Planck and high-resolution CMB data. None of these models are favoured significantly over standard LCDM. The deviation of the scalar spectral index from unity is insensitive to the addition of tensor modes and to changes in the matter content of the Universe. We find a 95% upper limit of r<0.11 on the tensor-to-scalar ratio. There is no evidence for additional neutrino-like relativistic particles. Using BAO and CMB data, we find N_eff=3.30+/-0.27 for the effective number of relativistic degrees of freedom, and an upper limit of 0.23 eV for the summed neutrino mass. Our results are in excellent agreement with big bang nucleosynthesis and the standard value of N_eff=3.046. We find no evidence for dynamical dark energy. Despite the success of the standard LCDM model, this cosmology does not provide a good fit to the CMB power spectrum at low multipoles, as noted previously by the WMAP team. While not of decisive significance, this is an anomaly in an otherwise self-consistent analysis of the Planck temperature data.

Within the standard LCDM model of cosmology, the recent Planck measurements have shown discrepancies with other observations, e.g., measurements of the current expansion rate H_0, the galaxy shear power spectrum and counts of galaxy clusters. We show that if LCDM is extended by a hot dark matter component, which could be interpreted as a sterile neutrino, the data sets can be combined consistently. A combination of Planck data, WMAP-9 polarisation data, measurements of the BAO scale, the HST measurement of H_0, Planck galaxy cluster counts and galaxy shear data from the CFHTLens survey yields Delta N_eff = 0.61 pm 0.30 and m_s^eff = (0.41 pm 0.13) eV at 1 sigma. The former is driven mainly by the large H_0 of the HST measurement, while the latter is driven by cluster data. CFHTLens galaxy shear data prefer Delta N_eff >0 and a non-zero mass. Taken together, we find hints for the presence of a hot dark matter component at 3 sigma. A sterile neutrino motivated by the reactor and gallium anomalies appears rejected at even higher significance and an accelerator anomaly sterile neutrino is found in tension at 2 sigma.

We perform a detailed forecast on how well a Euclid-like photometric galaxy and cosmic shear survey will be able to constrain the absolute neutrino mass scale. Adopting conservative assumptions about the survey specifications and assuming complete ignorance of the galaxy bias, we estimate that the minimum mass sum of sum m_nu ~ 0.06 eV in the normal hierarchy can be detected at 1.5 sigma to 2.5 sigma significance, depending on the model complexity, using a combination of galaxy and cosmic shear power spectrum measurements in conjunction with CMB temperature and polarisation observations from Planck. With better knowledge of the galaxy bias, the significance of the detection could potentially reach 5.4 sigma. Interestingly, neither Planck+shear nor Planck+galaxy alone can achieve this level of sensitivity; it is the combined effect of galaxy and cosmic shear power spectrum measurements that breaks the persistent degeneracies between the neutrino mass, the physical matter density, and the Hubble parameter. Notwithstanding this remarkable sensitivity to sum m_nu, Euclid-like shear and galaxy data will not be sensitive to the exact mass spectrum of the neutrino sector; no significant bias (< 1 sigma) in the parameter estimation is induced by fitting inaccurate models of the neutrino mass splittings to the mock data, nor does the goodness-of-fit of these models suffer any significant degradation relative to the true one (Delta chi_eff ^2< 1).

This white paper addresses the hypothesis of light sterile neutrinos based on recent anomalies observed in neutrino experiments and the latest astrophysical data.

A number of recent analyses of cosmological data have reported hints for the presence of extra radiation beyond the standard model expectation. In order to test the robustness of these claims under different methods of constructing parameter constraints, we perform a Bayesian posterior-based and a likelihood profile-based analysis of current data. We confirm the presence of a slight discrepancy between posterior- and profile-based constraints, with the marginalised posterior preferring higher values of the effective number of neutrino species N_eff. This can be traced back to a volume effect occurring during the marginalisation process, and we demonstrate that the effect is related to the fact that cosmic microwave background (CMB) data constrain N_eff only indirectly via the redshift of matter-radiation equality. Once present CMB data are combined with external information about, e.g., the Hubble parameter, the difference between the methods becomes small compared to the uncertainty of N_eff. We conclude that the preference of precision cosmological data for excess radiation is "real" and not an artifact of a specific choice of credible/confidence interval construction.

We investigate the question of whether the recent modulation signal claimed by CoGeNT is best explained by the dark matter (DM) hypothesis from a Bayesian model comparison perspective. We consider five phenomenological explanations for the data: no modulation signal, modulation due to DM, modulation due to DM compatible with the total CoGeNT rate, and a signal coming from other physics with a free phase but annual period, or with a free phase and a free period. In each scenario, we assign to the free parameters physically motivated priors. We find that while the DM models are weakly preferred to the no modulation model, but when compared to models where the modulation is due to other physics, the DM hypothesis is favoured with odds ranging from 185:1 to 560:1. This result is robust even when astrophysical uncertainties are taken into account and the impact of priors assessed. Interestingly, the odds for the DM model in which the modulation signal is compatible with the total rate against a DM model in which this prior is not implemented is only 5:8, in spite of the former's prediction of a modulation amplitude in the energy range 0.9 to 3.0 keVee that is significantly smaller than the value observed by CoGeNT. Classical hypothesis testing also rules out the null hypothesis of no modulation at the 1.6 sigma to 2.3 sigma level, depending on the details of the alternative. Lastly, we investigate whether anisotropic velocity distributions can help to mitigate the tension between the CoGeNT total and modulated rates, and find encouraging results.

We study cosmological models that contain sterile neutrinos with eV-range masses as suggested by reactor and short-baseline oscillation data. We confront these models with both precision cosmological data (probing the CMB decoupling epoch) and light-element abundances (probing the BBN epoch). In the minimal LambdaCDM model, such sterile neutrinos are strongly disfavoured by current data because they contribute too much hot dark matter. However, if the cosmological framework is extended to include also additional relativistic degrees of freedom -- beyond the three standard neutrinos and the putative sterile neutrinos, then the hot dark matter constraint on the sterile states is considerably relaxed. A further improvement is achieved by allowing a dark energy equation of state parameter w<-1. While BBN strongly disfavours extra radiation beyond the assumed eV-mass sterile neutrino, this constraint can be circumvented by a small nu_e degeneracy. Any model containing eV-mass sterile neutrinos implies also strong modifications of other cosmological parameters. Notably, the inferred cold dark matter density can shift up by 20 to 75% relative to the standard LambdaCDM value.

Bayesian statistical methods offer a simple and consistent framework for incorporating uncertainties into a multi-parameter inference problem. In this work we apply these methods to a selection of current direct dark matter searches. We consider the simplest scenario of spin-independent elastic WIMP scattering, and infer the WIMP mass and cross-section from the experimental data with the essential systematic uncertainties folded into the analysis. We find that when uncertainties in the scintillation efficiency of Xenon100 have been accounted for, the resulting exclusion limit is not sufficiently constraining to rule out the CoGeNT preferred parameter region, contrary to previous claims. In the same vein, we also investigate the impact of astrophysical uncertainties on the preferred WIMP parameters. We find that within the class of smooth and isotropic WIMP velocity distributions, it is difficult to reconcile the DAMA and the CoGeNT preferred regions by tweaking the astrophysics parameters alone. If we demand compatibility between these experiments, then the inference process naturally concludes that a high value for the sodium quenching factor for DAMA is preferred.

Precision cosmology and big-bang nucleosynthesis mildly favor extra radiation in the universe beyond photons and ordinary neutrinos, lending support to the existence of low-mass sterile neutrinos. We use the WMAP 7-year data, small-scale CMB observations from ACBAR, BICEP and QuAD, the SDSS 7th data release, and measurement of the Hubble parameter from HST observations to derive credible regions for the assumed common mass scale m_s and effective number N_s of thermally excited sterile neutrino states. Our results are compatible with the existence of one or perhaps two sterile neutrinos, as suggested by LSND and MiniBooNE, if m_s is in the sub-eV range.

The matter power spectrum as derived from large scale structure (LSS) surveys contains two important and distinct pieces of information: an overall smooth shape and the imprint of baryon acoustic oscillations (BAO). We investigate the separate impact of these two types of information on cosmological parameter estimation, and show that for the simplest cosmological models, the broad-band shape information currently contained in the SDSS DR7 halo power spectrum (HPS) is by far superseded by geometric information derived from the baryonic features. An immediate corollary is that contrary to popular beliefs, the upper limit on the neutrino mass m_\nu presently derived from LSS combined with cosmic microwave background (CMB) data does not in fact arise from the possible small-scale power suppression due to neutrino free-streaming, if we limit the model framework to minimal LambdaCDM+m_\nu. However, in more complicated models, such as those extended with extra light degrees of freedom and a dark energy equation of state parameter w differing from -1, shape information becomes crucial for the resolution of parameter degeneracies. This conclusion will remain true even when data from the Planck surveyor become available. In the course of our analysis, we introduce a new dewiggling procedure that allows us to extend consistently the use of the SDSS HPS to models with an arbitrary sound horizon at decoupling. All the cases considered here are compatible with the conservative 95%-bounds \sum m_\nu < 1.16 eV, N_eff = 4.8 \pm 2.0.

We present constraints on canonical single-field inflation derived from WMAP five year, ACBAR, QUAD, BICEP data combined with the halo power spectrum from SDSS LRG7. Models with a non-scale-invariant spectrum and a red tilt n_s < 1 are now preferred over the Harrison-Zel'dovich model (n_s = 1, tensor-to-scalar ratio r = 0) at high significance. Assuming no running of the spectral indices, we derive constraints on the parameters (n_s, r) and compare our results with the predictions of simple inflationary models. The marginalised credible intervals read n_s = 0.962^{+0.028}_{-0.026} and r < 0.17 (at 95% confidence level). Interestingly, the 68% c.l. contours favour mainly models with a convex potential in the observable region, but the quadratic potential model remains inside the 95% c.l. contours. We demonstrate that these results are robust to changes in the datasets considered and in the theoretical assumptions made. We then consider a non-vanishing running of the spectral indices by employing different methods, non-parametric but approximate, or parametric but exact. With our combination of CMB and LSS data, running models are preferred over power-law models only by a Delta chi^2 ~ 5.8, allowing inflationary stages producing a sizable negative running -0.063^{+0.061}_{-0.049} and larger tensor-scalar ratio r < 0.33 at the 95% c.l. This requires large values of the third derivative of the inflaton potential within the observable range. We derive bounds on this derivative under the assumption that the inflaton potential can be approximated as a third order polynomial within the observable range.

Features in the primordial power spectrum have been suggested as an explanation for glitches in the angular power spectrum of temperature anisotropies measured by the WMAP satellite. However, these glitches might just as well be artifacts of noise or cosmic variance. Using the effective Delta chi^2 between the best-fit power-law spectrum and a deconvolved primordial spectrum as a measure of "featureness" of the data, we perform a full Monte-Carlo analysis to address the question of how significant the recovered features are. We find that in 26% of the simulated data sets the reconstructed spectrum yields a greater improvement in the likelihood than for the actually observed data. While features cannot be categorically ruled out by this analysis, and the possibility remains that simple theoretical models which predict some of the observed features might stand up to rigorous statistical testing, our results suggest that WMAP data are consistent with the assumption of a featureless power-law primordial spectrum.

We discuss current cosmological constraints on axions, as well as future sensitivities. Bounds on axion hot dark matter are discussed first, and subsequently we discuss both current and future sensitivity to models in which axions play the role as cold dark matter, but where the Peccei-Quinn symmetry is not restored during reheating.

We explore the cosmological sensitivity to the amplitude of isocurvature fluctuations that would be caused by axions in the "anthropic window" where the axion decay constant f_a >> 10^12 GeV and the initial misalignment angle Theta_i << 1. In a minimal Lambda-CDM cosmology extended with subdominant scale-invariant isocurvature fluctuations, existing data constrain the isocurvature fraction to alpha < 0.09 at 95% C.L. If no signal shows up, Planck can improve this constraint to 0.042 while an ultimate CMB probe limited only by cosmic variance in both temperature and E-polarisation can reach 0.017, about a factor of five better than the current limit. In the parameter space of f_a and H_I (Hubble parameter during inflation) we identify a small region where axion detection remains within the reach of realistic cosmological probes.

High precision measurements of the Cosmic Microwave Background (CMB) anisotropies, as can be expected from the Planck satellite, will require high-accuracy theoretical predictions as well. One possible source of theoretical uncertainty is the numerical error in the output of the Boltzmann codes used to calculate angular power spectra. In this work, we carry out an extensive study of the numerical accuracy of the public Boltzmann code CAMB, and identify a set of parameters which determine the error of its output. We show that at the current default settings, the cosmological parameters extracted from data of future experiments like Planck can be biased by several tenths of a standard deviation for the six parameters of the standard Lambda-CDM model, and potentially more seriously for extended models. We perform an optimisation procedure that leads the code to achieve sufficient precision while at the same time keeping the computation time within reasonable limits. Our conclusion is that the contribution of numerical errors to the theoretical uncertainty of model predictions is well under control -- the main challenges for more accurate calculations of CMB spectra will be of an astrophysical nature instead.

There has been much recent discussion, and some confusion, regarding the use of existing observational data to estimate the likelihood that next-generation cosmic microwave background (CMB) polarization experiments might detect a nonzero tensor signal, possibly associated with inflation. We examine this issue in detail here in two different ways: (1) first we explore the effect of choice of different parameter priors on the estimation of the tensor-to-scalar ratio r and other parameters describing inflation, and (2) we examine the Bayesian complexity in order to determine how effectively existing data can constrain inflationary parameters. We demonstrate that existing data are not strong enough to render full inflationary parameter estimates in a parametrization- and prior-independent way and that the predicted tensor signal is particularly sensitive to different priors. For parametrizations where the Bayesian complexity is comparable to the number of free parameters we find that a flat prior on the scale of inflation (which is to be distinguished from a flat prior on the tensor-to-scalar ratio) leads us to infer a larger, and in fact slightly nonzero tensor contribution at 68% confidence level. However, no detection is claimed. Our results demonstrate that all that is statistically relevant at the current time is the (slightly enhanced) upper bound on r, and we stress that the data remain consistent with r = 0.

We revisit the issue of ripples in the primordial power spectra caused by trans-Planckian physics, and the potential for their detection by future cosmological probes. We find that for reasonably large values of the first slow-roll parameter epsilon (> 0.001), a positive detection of trans-Planckian ripples can be made even if the amplitude is as low as 10^-4. Data from the Large Synoptic Survey Telescope (LSST) and the proposed future 21 cm survey with the Fast Fourier Transform Telescope (FFTT) will be particularly useful in this regard. If the scale of inflation is close to its present upper bound, a scale of new physics as high as 0.2 M_Planck could lead to observable signatures.

We explore and compare the performances of two nonlinear correction and scale-dependent biasing models for the extraction of cosmological information from galaxy power spectrum data, especially in the context of beyond-LCDM cosmologies. The first model is the well known Q model, first applied in the analysis of 2dFGRS data. The second, the P model, is inspired by the halo model, in which nonlinear evolution and scale-dependent biasing are encapsulated in a single non-Poisson shot noise term. We find that while both models perform equally well in providing adequate correction for a range of galaxy clustering data in standard LCDM cosmology and in extensions with massive neutrinos, the Q model can give unphysical results in cosmologies containing a subdominant free-streaming dark matter whose temperature depends on the particle mass, e.g., relic thermal axions, unless a suitable prior is imposed on the correction parameter. This last case also exposes the danger of analytic marginalisation, a technique sometimes used in the marginalisation of nuisance parameters. In contrast, the P model suffers no undesirable effects, and is the recommended nonlinear correction model also because of its physical transparency.

We estimate the effect of the experimental uncertainty in the measurement of the temperature of the cosmic microwave background (CMB) on the extraction of cosmological parameters from future CMB surveys. We find that even for an ideal experiment limited only by cosmic variance up to l = 2500 for both the temperature and polarisation measurements, the projected cosmological parameter errors are remarkably robust against the uncertainty of 1 mK in the FIRAS instrument's CMB temperature monopole measurement. The maximum degradation in sensitivity is 20%, for the baryon density estimate, relative to the case in which the monopole is known infinitely well. While this degradation is acceptable, we note that reducing the uncertainty in the current temperature measurement by a factor of five will bring it down to the per cent level. We also estimate the effect of the uncertainty in the dipole temperature measurement. Assuming the overall calibration of the data to be dominated by the dipole error of 0.2% from FIRAS, the sensitivity degradation is insignificant and does not exceed 10% in any parameter direction.

Data from future high-precision Cosmic Microwave Background (CMB) measurements will be sensitive to the primordial Helium abundance $Y_p$. At the same time, this parameter can be predicted from Big Bang Nucleosynthesis (BBN) as a function of the baryon and radiation densities, as well as a neutrino chemical potential. We suggest to use this information to impose a self-consistent BBN prior on $Y_p$ and determine its impact on parameter inference from simulated Planck data. We find that this approach can significantly improve bounds on cosmological parameters compared to an analysis which treats $Y_p$ as a free parameter, if the neutrino chemical potential is taken to vanish. We demonstrate that fixing the Helium fraction to an arbitrary value can seriously bias parameter estimates. Under the assumption of degenerate BBN (i.e., letting the neutrino chemical potential $\xi$ vary), the BBN prior's constraining power is somewhat weakened, but nevertheless allows us to constrain $\xi$ with an accuracy that rivals bounds inferred from present data on light element abundances.

We update constraints on the Hubble function H(phi) during inflation, using the most recent cosmic microwave background (CMB) and large scale structure (LSS) data. Our main focus is on a comparison between various commonly used methods of calculating the primordial power spectrum via analytical approximations and the results obtained by integrating the exact equations numerically. In each case, we impose naive, minimally restrictive priors on the duration of inflation. We find that the choice of priors has an impact on the results: the bounds on inflationary parameters can vary by up to a factor two. Nevertheless, it should be noted that within the region allowed by the minimal prior of the exact method, the accuracy of the approximations is sufficient for current data. We caution however that a careless minimal implementation of the approximative methods allows models for which the assumptions behind the analytical approximations fail, and recommend using the exact numerical method for a self-consistent analysis of cosmological data.

We consider the inference of the cosmic radiation density, traditionally parameterised as the effective number of neutrino species N_eff, from precision cosmological data. Paying particular attention to systematic effects, notably scale-dependent biasing in the galaxy power spectrum, we find no evidence for a significant deviation of N_eff from the standard value of N_eff^0=3.046 in any combination of cosmological data sets, in contrast to some recent conclusions of other authors. The combination of all available data in the linear regime prefers, in the context of a ``vanilla+N_eff'' cosmological model, 1.1<N_eff<4.8 (95% C.L.) with a best-fit value of 2.6. Adding data at smaller scales, notably the Lyman-alpha forest, we find 2.2<N_eff<5.8 (95% C.L.) with 3.8 as the best fit. Inclusion of the Lyman-alpha data shifts the preferred N_eff upwards because the sigma_8 value derived from the SDSS Lyman-alpha data is inconsistent with that inferred from CMB. In an extended cosmological model that includes a nonzero mass for N_eff neutrino flavours, a running scalar spectral index and a w parameter for the dark energy, we find 0.8<N_eff<6.1 (95% C.L.) with 3.0 as the best fit.

We revisit the problem of constraining steps in the inflationary potential with cosmological data. We argue that a step in the inflationary potential produces qualitatively similar oscillations in the primordial power spectrum, independently of the details of the inflationary model. We propose a phenomenological description of these oscillations and constrain these features using a selection of cosmological data including the baryonic peak data from the correlation function of luminous red galaxies in the Sloan Digital Sky Survey. Our results show that degeneracies of the oscillation with standard cosmological parameters are virtually non-existent. The inclusion of new data severely tightens the constraints on the parameter space of oscillation parameters with respect to older work. This confirms that extensions to the simplest inflationary models can be successfully constrained using cosmological data.

High-precision data from observation of the cosmic microwave background and the large scale structure of the universe provide very tight constraints on the effective parameters that describe cosmological inflation. Indeed, within a constrained class of LambdaCDM models, the simple lambda phi^4 chaotic inflation model already appears to be ruled out by cosmological data. In this paper, we compute constraints on inflationary parameters within a more general framework that includes other physically motivated parameters such as a nonzero neutrino mass. We find that a strong degeneracy between the tensor-to-scalar ratio r and the neutrino mass prevents lambda phi^4 from being excluded by present data. Reversing the argument, if lambda phi^4 is the correct model of inflation, it predicts a sum of neutrino masses at 0.3-0.5 eV, a range compatible with present experimental limits and within the reach of the next generation of neutrino mass measurements. We also discuss the associated constraints on the dark matter density, the dark energy equation of state, and spatial curvature, and show that the allowed regions are significantly altered. Importantly, we find an allowed range of 0.094 < Omega_c h^2 < 0.136 for the dark matter density, a factor of two larger than that reported in previous studies. This expanded parameter space may have implications for constraints on SUSY dark matter models.

The new three year WMAP data seem to confirm the presence of non-standard large scale features in the Cosmic Microwave Anisotropies power spectrum. While these features may hint at uncorrected experimental systematics, it is also possible to generate, in a cosmological way, oscillations on large angular scales by introducing a sharp step in the inflaton potential. Using current cosmological data, we derive constraints on the position, magnitude and gradient of a possible step. We show that a step in the inflaton potential, while strongly constrained by current data, is still allowed and may provide an interesting explanation to the currently measured deviations from the standard featureless spectrum. Moreover, we show that inflationary oscillations in the primordial power spectrum can significantly bias parameter estimates from standard ruler methods involving measurements of baryon oscillations.