We consider a scenario where local Lorentz invariance is violated by the existence of a preferred time direction at every space-time point. This scenario can arise in the context of quantum gravity and its description at low energies contains a unit time-like vector field which parameterizes the preferred direction. The particle physics tests of Lorentz invariance preclude a direct coupling of this vector to the fields of the Standard Model, but do not bear implications for dark matter. We discuss how the presence of this vector and its possible coupling to dark matter affect the evolution of the Universe. At the level of homogeneous cosmology the only effect of Lorentz invariance violation is a rescaling of the expansion rate. The physics is richer at the level of perturbations. We identify three effects crucial for observations: the rescaling of the matter contribution to the Poisson equation, the appearance of an extra contribution to the anisotropic stress and the scale-dependent enhancement of dark matter clustering. These effects result in distinctive features in the power spectra of the CMB and density fluctuations. Making use of the data from Planck and WiggleZ we obtain the most stringent cosmological constraints to date on departures from Lorentz symmetry. Our analysis provides the first direct bounds on deviations from Lorentz invariance in the dark matter sector.

The existence of a cosmic neutrino background can be probed indirectly by CMB experiments, not only by measuring the background density of radiation in the universe, but also by searching for the typical signatures of the fluctuations of free-streaming species in the temperature and polarisation power spectrum. Previous studies have already proposed a rather generic parametrisation of these fluctuations, that could help to discriminate between the signature of ordinary free-streaming neutrinos, or of more exotic dark radiation models. Current data are compatible with standard values of these parameters, which seems to bring further evidence for the existence of a cosmic neutrino background. In this work, we investigate the robustness of this conclusion under various assumptions. We generalise the definition of an effective sound speed and viscosity speed to the case of massive neutrinos or other dark radiation components experiencing a non-relativistic transition. We show that current bounds on these effective parameters do not vary significantly when considering an arbitrary value of the particle mass, or extended cosmological models with a free effective neutrino number, dynamical dark energy or a running of the primordial spectrum tilt. We conclude that it is possible to make a robust statement about the detection of the cosmic neutrino background by CMB experiments.

Dark Matter is essential for structure formation in the late Universe so it must be stable on cosmological time scales. But how stable exactly? Only assuming decays into relativistic particles, we report an otherwise model independent bound on the lifetime of Dark Matter using current cosmological data. Since these decays affect only the low-$\ell$ multipoles of the CMB, the Dark Matter lifetime is expected to correlate with the tensor-to-scalar ratio $r$ as well as curvature $\Omega_k$. We consider two models, including $r$ and $r+\Omega_k$ respectively, versus data from Planck, WMAP, WiggleZ and Baryon Acoustic Oscillations, with or without the BICEP2 data (if interpreted in terms of primordial gravitational waves). This results in a lower bound on the lifetime of CDM given by 160Gyr (without BICEP2) or 200Gyr (with BICEP2) at 95% confidence level.

The apparent discrepancy between the value of the tensor-to-scalar ratio reported by the BICEP2 collaboration, $r = 0.20^{+0.07}_{-0.05}$ at 68% CL, and the Planck upper limit, $r < 0.11$ at 95% CL, has attracted a great deal of attention. In this short note, we show that this discrepancy is mainly due to an `apples to oranges' comparison. The result reported by BICEP2 was measured at a pivot scale $k_* = 0.05$ Mpc$^{-1}$, assuming $n_t = 0$, whereas the Planck limit was provided at $k_* = 0.002$ Mpc$^{-1}$, assuming the slow-roll consistency relation $n_t = -r/8$. One should obviously compare the BICEP2 and Planck results under the same circumstances. By imposing $n_t = 0$, the Planck constraint at $k_* = 0.05$ Mpc$^{-1}$ becomes $r < 0.135$ at $95%$ CL, which can be compared directly with the BICEP2 result. Once a plausible dust contribution to the BICEP2 signal is taken into account (DDM2 model), $r$ is reduced to $r = 0.16^{+0.06}_{-0.05}$ and the discrepancy becomes of order $1.3\sigma$ only.

Since the public release of Planck data, several attempts have been made to explain the observed small tensions with other data-sets, most of them involving an extension of the {\Lambda}CDM Model. We try here an alternative approach to the data analysis, based on separating the constraints coming from the different epochs in cosmology, in order to assess which part of the Standard Model generates the tension with the data. To this end, we perform a particular analysis of Planck data probing only the early cosmological evolution, until the time of photon decoupling. Then, we utilise this result to see if the {\Lambda}CDM model can fit all observational constraints probing only the late cosmological background evolution, discarding any information concerning the late perturbation evolution. We find that all tensions between the datasets are removed, suggesting that our standard assumptions on the perturbed late-time history, as well as on reionisation, could sufficiently bias our parameter extraction and be the source of the alleged tensions.

The role of Lorentz invariance as a fundamental symmetry of nature has been lately reconsidered in different approaches to quantum gravity. It is thus natural to study whether other puzzles of physics may be solved within these proposals. This may be the case for the cosmological constant problem. Indeed, it has been shown that breaking Lorentz invariance provides Lagrangians that can drive the current acceleration of the universe without experiencing large corrections from ultraviolet physics. In this work, we focus on the simplest model of this type, called ThetaCDM, and study its cosmological implications in detail. At the background level, this model cannot be distinguished from LambdaCDM. The differences appear at the level of perturbations. We show that in ThetaCDM, the spectrum of CMB anisotropies and matter fluctuations may be affected by a rescaling of the gravitational constant in the Poisson equation, by the presence of extra contributions to the anisotropic stress, and finally by the existence of extra clustering degrees of freedom. To explore these modifications accurately, we modify the Boltzmann code CLASS. We then use the parameter inference code Monte Python to confront ThetaCDM with data from WMAP-7, SPT and WiggleZ. We obtain strong bounds on the parameters accounting for deviations from LambdaCDM. In particular, we find that the discrepancy between the gravitational constants appearing in the Poisson and Friedmann equations is constrained at the level 1.8%.

The annihilation or decay of Dark Matter (DM) particles could affect the thermal history of the universe and leave an observable signature in Cosmic Microwave Background (CMB) anisotropies. We update constraints on the annihilation rate of DM particles in the smooth cosmological background, using WMAP7 and recent small-scale CMB data. With a systematic analysis based on the Press-Schechter formalism, we also show that DM annihilation in halos at small redshift may explain entirely the reionization patterns observed in the CMB, under reasonable assumptions concerning the concentration and formation redshift of halos. We find that a mixed reionization model based on DM annihilation in halos as well as star formation at a redshift z~6.5 could simultaneously account for CMB observations and satisfy constraints inferred from the Gunn-Peterson effect. However, these models tend to reheat the inter-galactic medium (IGM) well above observational bounds: by including a realistic prior on the IGM temperature at low redshift, we find stronger cosmological bounds on the annihilation cross-section than with the CMB alone.

Models for the latest stages of the cosmological evolution rely on a less solid theoretical and observational ground than the description of earlier stages like BBN and recombination. As suggested in a previous work by Vonlanthen et al., it is possible to tweak the analysis of CMB data in such way to avoid making assumptions on the late evolution, and obtain robust constraints on "early cosmology parameters". We extend this method in order to marginalise the results over CMB lensing contamination, and present updated results based on recent CMB data. Our constraints on the minimal early cosmology model are weaker than in a standard LCDM analysis, but do not conflict with this model. Besides, we obtain conservative bounds on the effective neutrino number and neutrino mass, showing no hints for extra relativistic degrees of freedom, and proving in a robust way that neutrinos experienced their non-relativistic transition after the time of photon decoupling. This analysis is also an occasion to describe the main features of the new parameter inference code Monte Python, that we release together with this paper. Monte Python is a user-friendly alternative to other public codes like CosmoMC, interfaced with the Boltzmann code class.

We present forecasts for the accuracy of determining the parameters of a minimal cosmological model and the total neutrino mass based on combined mock data for a future Euclid-like galaxy survey and Planck. We consider two different galaxy surveys: a spectroscopic redshift survey and a cosmic shear survey. We make use of the Monte Carlo Markov Chains (MCMC) technique and assume two sets of theoretical errors. The first error is meant to account for uncertainties in the modelling of the effect of neutrinos on the non-linear galaxy power spectrum and we assume this error to be fully correlated in Fourier space. The second error is meant to parametrize the overall residual uncertainties in modelling the non-linear galaxy power spectrum at small scales, and is conservatively assumed to be uncorrelated and to increase with the ratio of a given scale to the scale of non-linearity. It hence increases with wavenumber and decreases with redshift. With these two assumptions for the errors and assuming further conservatively that the uncorrelated error rises above 2% at k = 0.4 h/Mpc and z = 0.5, we find that a future Euclid-like cosmic shear/galaxy survey achieves a 1-sigma error on Mnu close to 32 meV/25 meV, sufficient for detecting the total neutrino mass with good significance. If the residual uncorrelated errors indeed rises rapidly towards smaller scales in the non-linear regime as we have assumed here then the data on non-linear scales does not increase the sensitivity to the total neutrino mass. Assuming instead a ten times smaller theoretical error with the same scale dependence, the error on the total neutrino mass decreases moderately from sigma(Mnu) = 18 meV to 14 meV when mildly non-linear scales with 0.1 h/Mpc < k < 0.6 h/Mpc are included in the analysis of the galaxy survey data.

We address the issue of computing the non-linear matter power spectrum on mildly non-linear scales with efficient semi-analytic methods. We implemented M. Pietroni's Time Renormalization Group (TRG) method and its Dynamical 1-Loop (D1L) limit in a numerical module for the new Boltzmann code CLASS. Our publicly released module is valid for LCDM models, and optimized in such a way to run in less than a minute for D1L, or in one hour (divided by number of nodes) for TRG. A careful comparison of the D1L, TRG and Standard 1-Loop approaches reveals that results depend crucially on the assumed initial bispectrum at high redshift. When starting from a common assumption, the three methods give roughly the same results, showing that the partial resumation of diagrams beyond one loop in the TRG method improves one-loop results by a negligible amount. A comparison with highly accurate simulations by M. Sato & T. Matsubara shows that all three methods tend to over-predict non-linear corrections by the same amount on small wavelengths. Percent precision is achieved until k~0.2 h/Mpc for z>2, or until k~0.14 h/Mpc at z=1.