We present a comprehensive analysis of atomic hydrogen (HI) properties using a semi-analytical model of galaxy formation and N-body simulations covering a large cosmological volume at high resolution. We examine the HI mass function and the HI density, characterizing both their redshift evolution and their dependence on hosting halo mass. We analyze the HI content of dark matter haloes in the local Universe and up to redshift $z=5$, discussing the contribution of different galaxy properties. We find that different assembly history plays a crucial role in the scatter of this relation. We propose new fitting functions useful for constructing mock HI maps with HOD techniques. We investigate the HI clustering properties relevant for future $21$~cm Intensity Mapping (IM) experiments, including the HI bias and the shot noise level. The HI bias increases with redshift and it is roughly flat on the largest scales probed. The scale dependency is found at progressively larger scales with increasing redshift, apart from a dip feature at $z=0$. The shot-noise values are consistent with the ones inferred by independent studies, confirming that shot-noise will not be a limiting factor for IM experiments. We detail the contribution from various galaxy properties on the HI power spectrum and their relation to the halo bias. We find that HI poor satellite galaxies play an important role at the scales of the 1-halo term. Finally, we present the $21$~cm signal in redshift space, a fundamental prediction to be tested against data from future radio telescopes such as SKA.

Global (i.e. sky-averaged) $21$~cm signal experiments can measure the evolution of the universe from the Cosmic Dawn to the Epoch of Reionization. These measurements are challenged by the presence of bright foreground emission that can be separated from the cosmological signal if its spectrum is smooth. This assumption fails in the case of single polarization antennas as they measure linearly polarized foreground emission - which is inevitably Faraday rotated through the interstellar medium. We investigate the impact of Galactic polarized foregrounds on the extraction of the global 21~cm signal through realistic sky and dipole simulations both in a low frequency band from $50$ to $100$~MHz, where a 21~cm absorption profile is expected, and in a higher frequency band ($100-200$~MHz). We find that the presence of a polarized contaminant with complex frequency structure can bias the amplitude and the shape of the reconstructed signal parameters in both bands. We investigate if polarized foregrounds can explain the unexpected $21$~cm Cosmic Dawn signal recently reported by the EDGES collaboration. We find that unaccounted polarized foreground contamination can produce an enhanced and distorted $21$~cm absorption trough similar to the anomalous profile reported by Bowman et al. (2018), and whose amplitude is in mild tension with the assumed input Gaussian profile (at $\sim 1.5 \sigma$ level). Moreover, we note that, under the hypothesis of contamination from polarized foreground, the amplitude of the reconstructed EDGES signal can be overestimated by around $30\%$, mitigating the requirement for an explanation based on exotic physics.

The detection of the redshifted cosmological $21$~cm line signal requires the removal of the Galactic and extragalactic foreground emission, which is orders of magnitude brighter anywhere in the sky. Foreground cleaning methods currently used are efficient in removing spectrally smooth components. However, they struggle in the presence of not spectrally smooth contamination that is, therefore, potentially the most dangerous one. An example of this is the polarized synchrotron emission, which is Faraday rotated by the interstellar medium and leaks into total intensity due to instrumental imperfections. In this work we present new full-sky simulations of this polarized synchrotron emission in the $50-200$~MHz range, obtained from the observed properties of diffuse polarized emission at low frequencies. The simulated polarized maps are made publicly available, aiming to provide more realistic templates to simulate the effect of instrumental leakage and the effectiveness of foreground separation techniques.

Recent precision cosmological parameter constraints imply that the spatial curvature of the Universe is essentially dynamically negligible - but only if relatively strong assumptions are made about the equation of state of dark energy (DE). When these assumptions are relaxed, strong degeneracies arise that make it hard to disentangle DE and curvature, degrading the constraints. We show that forthcoming 21cm intensity mapping experiments such as HIRAX are ideally designed to carry out model-independent curvature measurements, as they can measure the clustering signal at high redshift with sufficient precision to break many of the degeneracies. We consider two different model-independent methods, based on `avoiding' the DE-dominated regime and non-parametric modelling of the DE equation of state respectively. Our forecasts show that HIRAX will be able to improve upon current model-independent constraints by around an order of magnitude, reaching percent-level accuracy even when an arbitrary DE equation of state is assumed. In the same model-independent analysis, the sample variance limit for a similar survey is another order of magnitude better.

When combining cosmological and oscillations results to constrain the neutrino sector, the question of the propagation of systematic uncertainties is often raised. We address this issue in the context of the derivation of an upper bound on the sum of the neutrino masses ($\Sigma m_\nu$) with recent cosmological data. This work is performed within the ${{\mathrm{\Lambda{CDM}}}}$ model extended to $\Sigma m_\nu$, for which we advocate the use of three mass-degenerate neutrinos. We focus on the study of systematic uncertainties linked to the foregrounds modelling in CMB data analysis, and on the impact of the present knowledge of the reionisation optical depth. This is done through the use of different likelihoods built from Planck data. Limits on $\Sigma m_\nu$ are derived with various combinations of data, including the latest Baryon Acoustic Oscillations (BAO) and Type Ia Supernovae (SN) results. We also discuss the impact of the preference for current CMB data for amplitudes of the gravitational lensing distortions higher than expected within the ${{\mathrm{\Lambda{CDM}}}}$ model, and add the Planck CMB lensing. We then derive a robust upper limit: $\Sigma m_\nu< 0.17\hbox{ eV at }95\% \hbox{CL}$, including 0.01 eV of foreground systematics. We also discuss the neutrino mass repartition and show that today's data do not allow one to disentangle normal from inverted hierarchy. The impact on the other cosmological parameters is also reported, for different assumptions on the neutrino mass repartition, and different high and low multipole CMB likelihoods.

We demonstrate that the cosmic microwave background (CMB) temperature-polarization cross-correlation provides accurate and robust constraints on cosmological parameters. We compare them with the results from temperature or polarization and investigate the impact of foregrounds, cosmic variance, and instrumental noise. This analysis makes use of the Planck high-multipole HiLLiPOP likelihood based on angular power spectra, which takes into account systematics from the instrument and foreground residuals directly modelled using Planck measurements. The temperature-polarization correlation (TE) spectrum is less contaminated by astrophysical emissions than the temperature power spectrum (TT), allowing constraints that are less sensitive to foreground uncertainties to be derived. For {\Lambda}CDM parameters, TE gives very competitive results compared to TT. For basic {\Lambda}CDM model extensions (such as AL, {\Sigma}m{\nu}, or Neff ), it is still limited by the instrumental noise level in the polarization maps.

The angular power spectra of the cosmic microwave background (CMB) temperature anisotropies reconstructed from Planck data seem to present too much gravitational lensing distortion. This is quantified by the control parameter $A_L$ that should be compatible with unity for a standard cosmology. With the Class Boltzmann solver and the profile-likelihood method, for this parameter we measure a 2.6$\sigma$ shift from 1 using the Planck public likelihoods. We show that, owing to strong correlations with the reionization optical depth $\tau$ and the primordial perturbation amplitude $A_s$, a $\sim2\sigma$ tension on $\tau$ also appears between the results obtained with the low ($\ell\leq 30$) and high ($30<\ell\lesssim 2500$) multipoles likelihoods. With Hillipop, another high-$\ell$ likelihood built from Planck data, this difference is lowered to $1.3\sigma$. In this case, the $A_L$ value is still in disagreement with unity by $2.2\sigma$, suggesting a non-trivial effect of the correlations between cosmological and nuisance parameters. To better constrain the nuisance foregrounds parameters, we include the very high $\ell$ measurements of the Atacama Cosmology Telescope (ACT) and South Pole Telescope (SPT) experiments and obtain $A_L = 1.03 \pm 0.08$. The Hillipop+ACT+SPT likelihood estimate of the optical depth is $\tau=0.052\pm{0.035,}$ which is now fully compatible with the low $\ell$ likelihood determination. After showing the robustness of our results with various combinations, we investigate the reasons for this improvement that results from a better determination of the whole set of foregrounds parameters. We finally provide estimates of the $\Lambda$CDM parameters with our combined CMB data likelihood.

Cosmological parameter estimation is traditionally performed in the Bayesian context. By adopting an "agnostic" statistical point of view, we show the interest of confronting the Bayesian results to a frequentist approach based on profile-likelihoods. To this purpose, we have developed the Cosmological Analysis with a Minuit Exploration of the Likelihood ("CAMEL") software. Written from scratch in pure C++, emphasis was put in building a clean and carefully-designed project where new data and/or cosmological computations can be easily included. CAMEL incorporates the latest cosmological likelihoods and gives access from the very same input file to several estimation methods: (i) A high quality Maximum Likelihood Estimate (a.k.a "best fit") using MINUIT ; (ii) profile likelihoods, (iii) a new implementation of an Adaptive Metropolis MCMC algorithm that relieves the burden of reconstructing the proposal distribution. We present here those various statistical techniques and roll out a full use-case that can then used as a tutorial. We revisit the $\Lambda$CDM parameters determination with the latest Planck data and give results with both methodologies. Furthermore, by comparing the Bayesian and frequentist approaches, we discuss a "likelihood volume effect" that affects the optical reionization depth when analyzing the high multipoles part of the Planck data. The software, used in several Planck data analyzes, is available from http://camel.in2p3.fr. Using it does not require advanced C++ skills.

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.

The production of a primordial stochastic gravitational-wave background by processes occuring in the early Universe is expected in a broad range of models. Observing this background would open a unique window onto the Universe's evolutionary history. Probes like the Cosmic Microwave Background (CMB) or the Baryon Acoustic Oscillations (BAO) can be used to set upper limits on the stochastic gravitational-wave background energy density $\Omega_{GW}$ for frequencies above $10^{-15}$ Hz. We perform a profile likelihood analysis of the Planck CMB temperature anisotropies and gravitational lensing data combined with WMAP low-$\ell$ polarization, BAO, South Pole Telescope and Atacama Cosmology Telescope data. We find that $\Omega_{GW}h_{0}^{2} < 3.8 \times 10^{-6}$ at 95\% confidence level for adiabatic initial conditions which improves over the previous limit by a factor 2.3. Assuming that the primordial gravitational waves have been produced by a network of cosmic strings, we have derived exclusion limits in the cosmic string parameter space. If the size of the loops is determined by gravitational back-reaction, string tension values greater than $\sim 4 \times 10^{-9}$ are excluded for a reconnection probability of $10^{-3}$.

We explore the 2013 Planck likelihood function with a high-precision multi-dimensional minimizer (Minuit). This allows a refinement of the Lambda-cdm best-fit solution with respect to previously-released results, and the construction of frequentist confidence intervals using profile likelihoods. The agreement with the cosmological results from the Bayesian framework is excellent, demonstrating the robustness of the Planck results to the statistical methodology. We investigate the inclusion of neutrino masses, where more significant differences may appear due to the non-Gaussian nature of the posterior mass distribution. By applying the Feldman--Cousins prescription, we again obtain results very similar to those of the Bayesian methodology. However, the profile-likelihood analysis of the CMB combination (Planck+WP+highL) reveals a minimum well within the unphysical negative-mass region. We show that inclusion of the Planck CMB-lensing information regularizes this issue, and provide a robust frequentist upper limit $M_\nu < 0.26 eV$ ($95%$ confidence) from the CMB+lensing+BAO data combination.