We provide exact analytic expressions for the density, pressure, average number density and pseudo-pressure for massive neutrinos and generic dark matter particles, both fermions and bosons. We then focus on massive neutrinos and we compare our analytic expressions with the numerical implementation in the CLASS Boltzmann code. We find that our modifications including the exact analytic expressions are in agreement to better than $10^{-4}\%$ with the default CLASS implementation in the estimation of the CMB power spectrum; our modifications do not have an impact on the performance of the code. We also provide several specific limits of our expressions at the relativistic regime, but also at late times for the neutrino equation of state.

We present a family of designer Horndeski models, i.e. models that have a background exactly equal to that of the $\Lambda$CDM model but perturbations given by the Horndeski theory. Then, we extend the effective fluid approach to Horndeski theories, providing simple analytic formulae for the equivalent dark energy effective fluid pressure, density and velocity. We implement the dark energy effective fluid formulae in our code EFCLASS, a modified version of the widely used Boltzmann solver CLASS, and compare the solution of the perturbation equations with those of the code hi_CLASS which already includes Horndeski models. We find that our simple modifications to the vanilla code are accurate to the level of $\sim 0.1\%$ with respect to the more complicated hi_CLASS code. Furthermore, we study the kinetic braiding model both on and off the attractor and we find that even though the full case has a proper $\Lambda$CDM model limit for large $n$, it is not appropriately smooth, thus causing the quasistatic approximation to break down. Finally, we focus on our designer model (HDES), which has both a smooth $\Lambda$CDM limit and well-behaved perturbations, and we use it to perform Markov Chain Monte Carlo analyses to constrain its parameters with the latest cosmological data. We find that our HDES model can also alleviate the soft $2\sigma$ tension between the growth data and Planck 18 due to a degeneracy between $\sigma_8$ and one of its model parameters that indicates the deviation from the $\Lambda$CDM model.

Analyses of upcoming galaxy surveys will require careful modelling of relevant observables such as the power spectrum of galaxy counts in harmonic space $C_\ell(z,z')$. We investigate the impact of disregarding relevant relativistic effects by considering a model of dark energy including constant sound speed, constant equation of state $w_0$, and anisotropic stress. Here we show that neglecting the effect of lensing convergence will lead to substantial shifts in cosmological parameters such as the galaxy bias $b_0$, the value of the dark energy equation of state today $w_0$, and the Hubble constant $H_0$. Interestingly, neglecting the effect of lensing convergence in this kind of models results in shifting $H_0$ downwards, something which could shed light on the current tension between local and CMB determinations of the Hubble constant.

We provide explicit formulas for the effective fluid approach of $f(R)$ theories, such as the Hu & Sawicki and the designer models. Using the latter and simple modifications to the CLASS code, which we call EFCLASS, in conjunction with very accurate analytic approximations for the background evolution, we obtain competitive results in a much simpler and less error-prone approach. We also derive the initial conditions in matter domination and we find they differ from those already found in the literature for a constant $w$ model. A clear example is the designer model that behaves as $\Lambda$CDM in the background, but has nonetheless dark energy perturbations. We then use the aforementioned models to derive constraints from the latest cosmological data, including supernovae, BAO, CMB, $H(z)$ and growth-rate data, and find they are statistically consistent to the $\Lambda$CDM model. Finally, we show that the viscosity parameter $c_{vis}^2$ in realistic models is not constant as commonly assumed, but rather evolves significantly over several orders of magnitude, something which could affect forecasts of upcoming surveys.

We re-analyse recent Cepheid data to estimate the Hubble parameter $H_0$ by using Bayesian hyper-parameters (HPs). We consider the two data sets from Riess et al 2011 and 2016 (labelled R11 and R16, with R11 containing less than half the data of R16) and include the available anchor distances (megamaser system NGC4258, detached eclipsing binary distances to LMC and M31, and MW Cepheids with parallaxes), use a weak metallicity prior and no period cut for Cepheids. We find that part of the R11 data is down-weighted by the HPs but that R16 is mostly consistent with expectations for a Gaussian distribution, meaning that there is no need to down-weight the R16 data set. For R16, we find a value of $H_0 = 73.75 \pm 2.11 \, \mathrm{km} \, \mathrm{s}^{-1} \, \mathrm{Mpc}^{-1}$ if we use HPs for all data points (including Cepheid stars, supernovae type Ia, and the available anchor distances), which is about 2.6 $\sigma$ larger than the Planck 2015 value of $H_0 = 67.81 \pm 0.92 \,\mathrm{km}\, \mathrm{s}^{-1} \, \mathrm{Mpc}^{-1}$ and about 3.1 $\sigma$ larger than the updated Planck 2016 value $66.93 \pm 0.62 \,\mathrm{km}\, \mathrm{s}^{-1} \, \mathrm{Mpc}^{-1}$. We test the effect of different assumptions, and find that the choice of anchor distances affects the final value significantly. If we exclude the Milky Way from the anchors, then the value of $H_0$ decreases. We find however no evident reason to exclude the MW data. The HP method used here avoids subjective rejection criteria for outliers and offers a way to test datasets for unknown systematics.

We demonstrate the importance of including the lensing contribution in galaxy clustering analyses with large galaxy redshift surveys. It is well known that radial cross-correlations between different redshift bins of galaxy surveys are dominated by lensing. But we show here that also neglecting lensing in the autocorrelations within one bin severely biases cosmological parameter estimation with redshift surveys. It leads to significant shifts for several cosmological parameters, most notably the scalar spectral index and the neutrino mass scale. Especially the latter parameter is one of the main targets of future galaxy surveys.

We study a dark energy model with non-zero anisotropic stress, either linked to the dark energy density or to the dark matter density. We compute approximate solutions that allow to characterise the behaviour of the dark energy model and to assess the stability of the perturbations. We also determine the current limits on such an anisotropic stress from the cosmic microwave background data by the Planck satellite, and derive the corresponding constraints on the modified growth parameters like the growth index, the effective Newton's constant and the gravitational slip.

A convincing detection of primordial non-Gaussianity in the cosmic background radiation (CMB) is essential to probe the physics of the early universe. Since a single statistical estimator can hardly be suitable to detect the various possible forms of non-Gaussianity, it is important to employ different statistical indicators to study non-Gaussianity of CMB. This has motivated the proposal of a number statistical tools, including two large-angle indicators based on skewness and kurtosis of spherical caps of CMB sky-sphere. Although suitable to detect fairly large non-Gaussianity they are unable to detect non-Gaussianity within the Planck bounds, and exhibit power spectra with undesirable oscillation pattern. Here we use several thousands simulated CMB maps to examine interrelated problems regarding advances of these spherical patches procedures. We examine whether a change in the choice of the patches could enhance the sensitivity of the procedures well enough to detect large-angle non-Gaussianity within the Planck bounds. To this end, a new statistical procedure with non-overlapping cells is proposed and its capability is established. We also study whether this new procedure is capable to smooth out the undesirable oscillation pattern in the skewness and kurtosis power spectra of the spherical caps procedure. We show that the new procedure solves this problem, making clear this unexpected power spectra pattern does not have a physical origin, but rather presumably arises from the overlapping obtained with the spherical caps approach. Finally, we make a comparative analysis of this new statistical procedure with the spherical caps routine, determine their lower bounds for non-Gaussianity detection, and make apparent their relative strength and sensitivity.

A detection or non detection of primordial non--Gaussianity (NG) by using the cosmic microwave background radiation (CMB) is a possible way to break the degeneracy of early universe models. Since a single statistical estimator hardly can be sensitive to all possible forms of NG which may be present in the data, it is important to use different statistical estimators to study NG in CMB. Recently, two new large-angle NG indicators based on skewness and kurtosis of spherical caps or spherical cells of CMB sky have been proposed and used in both CMB data and simulated maps. Here, we make a comparative study of these two different procedures by examining the NG in the WMAP seven years ILC map. We show that the spherical cells procedure detects a higher level of NG than that obtained by the method with overlapping spherical caps.