A new perspective on dark energy modeling via genetic algorithms

We use Genetic Algorithms to extract information from several cosmological probes, such as the type Ia supernovae (SnIa), the Baryon Acoustic Oscillations (BAO) and the growth rate of matter perturbations. This is done by implementing a model independent and bias-free reconstruction of the various scales and distances that characterize the data, like the luminosity d(L)(z) and the angular diameter distance d(A)(z) in the SnIa and BAO data, respectively, or the dependence with redshift of the matter density Omega(m) (a) in the growth rate data, f sigma(8)(z). These quantities can then be used to reconstruct the expansion history of the Universe, and the resulting Dark Energy (DE) equation of state w(z) in the context of FRW models, or the mass radial function Omega(M)(r) in LTB models. In this way, the reconstruction is completely independent of our prior bias. Furthermore, we use this method to test the Etherington relation, ie the well-known relation between the luminosity and the angular diameter distance, eta equivalent to d(L)(z)/(1+z)(2)d(A)(z), which is equal to 1 in metric theories of gravity. We find that the present data seem to suggest a 3-sigma deviation from one at redshifts z similar to 0.5. Finally, we present a novel way, within the Genetic Algorithm paradigm, to analytically estimate the errors on the reconstructed quantities by calculating a Path Integral over all possible functions that may contribute to the likelihood. We show that this can be done regardless of the data being correlated or uncorrelated with each other and we also explicitly demonstrate that our approach is in good agreement with other error estimation techniques like the Fisher Matrix approach and the Bootstrap Monte Carlo.