Machine learning (ML) algorithms have revolutionized the way we interpret data in astronomy, particle physics, biology and even economics, since they can remove biases due to a priori chosen models. Here we apply a particular ML method, the genetic algorithms (GA), to cosmological data that describes the background expansion of the Universe, namely the Pantheon Type Ia supernovae and the Hubble expansion history $H(z)$ datasets. We obtain model independent and nonparametric reconstructions of the luminosity distance $d_L(z)$ and Hubble parameter $H(z)$ without assuming any dark energy model or a flat Universe. We then estimate the deceleration parameter $q(z)$, a measure of the acceleration of the Universe, and we make a $\sim4.5\sigma$ model independent detection of the accelerated expansion, but we also place constraints on the transition redshift of the acceleration phase $(z_{\textrm{tr}}=0.662\pm0.027)$. We also find a deviation from $\Lambda$CDM at high redshifts, albeit within the errors, hinting toward the recently alleged tension between the SnIa/quasar data and the cosmological constant $\Lambda$CDM model at high redshifts $(z\gtrsim1.5)$. Finally, we show the GA can be used in complementary null tests of the $\Lambda$CDM via reconstructions of the Hubble parameter and the luminosity distance.

The exact nature of dark energy is currently unknown and its cosmological perturbations, when dark energy is assumed not to be the cosmological constant, are usually modelled as adiabatic. Here we explore the possibility that dark energy might have a non-adiabatic component and we examine how it would affect several key cosmological observables. We present analytical solutions for the growth-rate and growth index of matter density perturbations and compare them to both numerical solutions of the fluid equations and an implementation in the Boltzmann code CLASS, finding they all agree to well below one percent. We also perform a Monte-Carlo analysis to derive constraints on the parameters of the non-adiabatic component using the latest cosmological data, including the temperature and polarization spectra of the Cosmic Microwave Background as observed by Planck, the Baryon Acoustic Oscillations, the Pantheon type Ia supernovae compilation and lastly, measurements of Redshift Space Distortions (RSD) of the growth-rate of matter perturbations. We find that the amplitude of the non-adiabatic pressure perturbation is consistent with zero within $1\sigma$. Finally, we also present a new, publicly available, RSD likelihood for MontePython based on the "Gold 2018" growth-rate data compilation.

We present a model independent and non-parametric reconstruction with a Machine Learning algorithm of the redshift evolution of the Cosmic Microwave Background (CMB) temperature from a wide redshift range $z\in \left[0,3\right]$ without assuming any dark energy model, an adiabatic universe or photon number conservation. In particular we use the genetic algorithms which avoid the dependency on an initial prior or a cosmological fiducial model. Through our reconstruction we constrain new physics at late times. We provide novel and updated estimates on the $\beta$ parameter from the parametrisation $\text{T}(z)=\text{T}_0(1+z)^{1-\beta}$, the duality relation $\eta(z)$ and the cosmic opacity parameter $\tau(z)$. Furthermore we place constraints on a spatial varying fine structure constant $\alpha$, which would have signatures in a broad spectrum of physical phenomena such as the CMB anisotropies. Overall we find no evidence of deviations within the $1\sigma$ region from the well established $\Lambda\text{CDM}$ model, thus confirming its predictive potential.

Recent analyses of the Planck data and quasars at high redshifts have suggested possible deviations from the flat $\Lambda$ cold dark matter model ($\Lambda$CDM), where $\Lambda$ is the cosmological constant. Here, we use machine learning methods to investigate any possible deviations from $\Lambda$CDM at both low and high redshifts by using the latest cosmological data. Specifically, we apply the genetic algorithms to explore the nature of dark energy (DE) in a model independent fashion by reconstructing its equation of state $w(z)$, the growth index of matter density perturbations $\gamma(z)$, the linear DE anisotropic stress $\eta_{DE}(z)$ and the adiabatic sound speed $c_{s,DE}^2(z)$ of DE perturbations. We find a $\sim2\sigma$ deviation of $w(z)$ from -1 at high redshifts, the adiabatic sound speed is negative at the $\sim2\sigma$ level and a $\sim3.5\sigma$ deviation of the anisotropic stress from unity at low redshifts and $\sim5.5 \sigma$ at high redshifts. These results suggest either the presence of a strong non-adiabatic component in the DE sound speed or the presence of DE anisotropic stress, thus hinting at possible deviations from the $\Lambda$CDM model.

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.