We explore the possibility that Dark Matter (DM) may be explained by a non-uniform background of approximately stellar-mass clusters of Primordial Black Holes (PBHs), by simulating the evolution them from recombination to the present with over 5000 realisations using a Newtonian $ N $-body code. We compute the cluster rate of evaporation, and extract the binary and merged sub-populations along with their parent and merger tree histories, lifetimes and formation rates; the dynamical and orbital parameter profiles, the degree of mass segregation and dynamical friction, and power spectrum of close encounters. Overall, we find that PBHs can constitute a viable DM candidate, and that their clustering presents a rich phenomenology throughout the history of the Universe. We show that binary systems constitute about 9.5% of all PBHs at present, with mass ratios of $ \bar{q}_{\mathrm{B}} = 0.154 $, and total masses of $ \bar{m}_{\mathrm{T},\mathrm{B}} = 303 \mathrm{M_\odot}$. Merged PBHs are rare, about 0.0023% of all PBHs at present, with mass ratios of $ \bar{q}_{\mathrm{M}} = 0.965 $ with total and chirp masses of $ \bar{m}_{\mathrm{T},\mathrm{M}} = 1670 \mathrm{M_\odot} $ and $ \bar{m}_{c,\mathrm{M}} = 642 \mathrm{M_\odot} $ respectively. We find that cluster puffing up and evaporation leads to bubbles of these PBHs of order 1 kpc containing at present times about 36% of objects and mass, with hundred pc sized cores. We also find that these PBH sub-haloes are distributed in wider PBH haloes of order hundreds of kpc, containing about 63% of objects and mass, coinciding with the sizes of galactic halos. We find at last high rates of close encounters of massive Black Holes ($ M \sim 1000 \mathrm{M_\odot} $), with $ \Gamma^{\mathrm{S}} = (1.2 \substack {+5.9 \\ -0.9}) \times 10^{7} \mathrm{yr^{-1} Gpc^{-3}}$ and mergers with $\Gamma^{\mathrm{M}} = 1337 \pm 41 \mathrm{yr^{-1} Gpc^{-3}} $.

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.

We test the cosmological implications of a varying absolute magnitude of Type Ia supernovae using the Pantheon compilation, by reconstructing different phenomenological approaches that could justify a varying absolute magnitude, but also approaches based on cosmic voids, modified gravity models and a binning scheme. In all the cases considered in this work, we find good agreement with the expected values of the standard $\Lambda$CDM model and no evidence of new physics.

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.

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.

If light enough primordial black holes (PBH) account for dark matter, then its density decreases with time as they lose mass via Hawking radiation. We show that this time-dependence of the matter density can be formulated as an equivalent $w(z)$ dark energy model and we study its implications on the expansion history. Using our approach and comparing with the latest cosmological data, including the supernovae type Ia, Baryon Acoustic Oscillations, Cosmic Microwave Background and the Hubble expansion H(z) data, we place observational constraints on the PBH model. We find that it is statistically consistent with $\Lambda$CDM according to the AIC statistical tool. Furthermore, we entertain the idea of having a population of ultra-light PBHs, decaying around neutrino decoupling, on top of the dark matter fluid and show how this offers a natural dark matter-radiation coupling altering the expansion history of the Universe and alleviating the $H_0$ tension.

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.

The Universe on comoving scales larger than 100 Mpc/h is assumed to be statistically homogeneous, with the transition scale where we have $1\%$ deviation from homogeneity being known as the homogeneity scale $R_H$. The latter was recently proposed in Ntelis et al. (2018) [arXiv:1810.09362] to be a standard ruler. In this paper we perform a comparison between the baryon acoustic oscillations and the homogeneity scales $R_H$, assuming a spatially flat universe, linear cosmological perturbation theory and the cosmological constant cold dark matter $\Lambda$CDM model. By direct theoretical calculations, we demonstrate that the answer to the question of whether the homogeneity scale can be used as a standard ruler, is clearly negative due to the non-monotonicity of the homogeneity scale $R_H$ with the matter density parameter $\Omega_{m0}$ of the $\Lambda$CDM model. Furthermore, we also consider the effect of redshift space distortions and we find that they do not break the degeneracies but instead only change the value of $R_H$ by a few percent. Finally, we also confirm our findings with the help of an N-Body simulation.

The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing.

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 perform an Internal Robustness analysis (iR) to a compilation of the most recent $f\sigma_8(z)$ data, using the framework of 1209.1897. The method analyzes combinations of subsets in the data set in a Bayesian model comparison way, potentially finding outliers, subsets of data affected by systematics or new physics. In order to validate our analysis and assess its sensitivity we performed several cross-checks, for example by removing some of the data or by adding artificially contaminated points, while we also generated mock data sets in order to estimate confidence regions of the iR. Applying this methodology, we found no anomalous behavior in the $f\sigma_8(z)$ data set, thus validating its internal robustness.

We extract observational constraints on $f(T)$ gravity, using the recently proposed statistical method which is not affected by the value of $H_0$ and thus it bypasses the problem of the disagreement in its exact numerical value between Planck and direct measurements. We use direct measurements of the Hubble parameter with the corresponding covariance matrix, and for completeness we perform a joint analysis using the latest data from Supernovae type Ia based on JLA sample, quasi-stellar objects, and Cosmic Microwave Background shift parameter from Planck. We analyze a large family of $f(T)$ models, and we compare the fitting results with $\Lambda$CDM cosmology using the AIC statistical test. Utilizing only the Hubble parameter data we find that in the case of the power-law $f(T)$ model a small but non-zero deviation from $\Lambda$CDM cosmology is slightly favored at 1-$\sigma$, nevertheless the corresponding AIC value shows a statistical equivalence with it. Finally, the join analysis reveals that all $f(T)$ models are very efficient and in very good agreement with observations.

In this paper we place observational constraints on the well-known $\gamma$-gravity $f(R)$ model using the latest cosmological data, namely we use the latest growth rate, Cosmic Microwave Background, Baryon Acoustic Oscillations, Supernovae type Ia and Hubble parameter data. Performing a joint likelihood analysis we find that the $\gamma$-gravity model is in very good agreement with observations. Utilizing the AIC statistical test we statistically compare the current $f(R)$ model with $\Lambda$CDM cosmology and find that they are statistically equivalent. Therefore, $\gamma$-gravity can be seen as a useful scenario toward testing deviations from General Relativity. Finally, we note that we find somewhat higher values for the $f(R)$ best-fit values compared to those mentioned in the past in the literature and we argue that this could potential alleviate the halo-mass function problem.

We describe in detail gravitational wave bursts from Primordial Black Hole (PBH) hyperbolic encounters. The bursts are one-time events, with the bulk of the released energy happening during the closest approach, which can be emitted in frequencies that could be within the range of both LIGO (10-1000Hz) and LISA ($10^{-6}-1$ Hz). Furthermore, we correct the results for the power spectrum of hyperbolic encounters found in the literature and present new exact and approximate expressions for the peak frequency of the emission. Note that these GW bursts from hyperbolic encounters between PBH are complementary to the GW emission from the bounded orbits of BHB mergers detected by LIGO, and help breaking degeneracies in the determination of the PBH mass, spin and spatial distributions.

In this paper we present cosmological constraints on several well-known $f(R)$ models, but also on a new class of models that are variants of the Hu-Sawicki one of the form $f(R)=R-\frac{2\Lambda}{1+b\;y(R,\Lambda)}$, that interpolate between the cosmological constant model and a matter dominated universe for different values of the parameter $b$, which is usually expected to be small for viable models and which in practice measures the deviation from General Relativity. We use the latest growth rate, Cosmic Microwave Background, Baryon Acoustic Oscillations, Supernovae type Ia and Hubble parameter data to place stringent constraints on the models and to compare them to the cosmological constant model but also other viable $f(R)$ models such as the Starobinsky or the degenerate hypergeometric models. We find that these kinds of Hu-Sawicki variant parameterizations are in general compatible with the currently available data and can provide useful toy models to explore the available functional space of $f(R)$ models, something very useful with the current and upcoming surveys that will test deviations from General Relativity.

We propose that Gravitational Wave (GW) bursts with millisecond durations can be explained by the GW emission from the hyperbolic encounters of Primordial Black Holes in dense clusters. These bursts are single events, with the bulk of the released energy happening during the closest approach, and emitted in frequencies within the AdvLIGO sensitivity range. We provide expressions for the shape of the GW emission in terms of the peak frequency and amplitude, and estimate the rates of these events for a variety of mass and velocity configurations. We study the regions of parameter space that will allow detection by both AdvLIGO and, in the future, LISA. We find for realistic configurations, with total mass $M\sim60\,M_\odot$, relative velocities $v\sim 0.01\,c$, and impact parameters $b\sim10^{-3}$ AU, for AdvLIGO an expected event rate is ${\cal O}(10)$ events/yr/Gpc$^3$ with millisecond durations. For LISA, the typical duration is in the range of minutes to hours and the event-rate is ${\cal O}(10^3)$ events/yr/Gpc$^3$ for both $10^3\,M_\odot$ IMBH and $10^6\,M_\odot$ SMBH encounters. We also study the distribution functions of eccentricities, peak frequencies and characteristic timescales that can be expected for a population of scattering PBH with a log-normal distribution in masses, different relative velocities and a flat prior on the impact parameter.

In this paper we present conjoined constraints on several cosmological models from the expansion history $H(z)$ and cosmic growth $f\sigma_8(z)$. The models we study include the CPL $w_0w_a$ parametrization, the Holographic Dark Energy (HDE) model, the Time varying vacuum ($\Lambda_t$CDM) model, the Dvali, Gabadadze and Porrati (DGP) and Finsler-Randers (FRDE) models, a power law $f(T)$ model and finally the Hu-Sawicki $f(R)$ model. In all cases we perform a simultaneous fit to the SnIa, CMB, BAO, $H(z)$ and growth data, while also following the conjoined visualization of $H(z)$ and $f\sigma_8(z)$ as in Linder (2017). Furthermore, we introduce the Figure of Merit (FoM) in the $H(z)-f\sigma_8(z)$ parameter space as a way to constrain models that jointly fit both probes well. We use both the latest $H(z)$ and $f\sigma_8(z)$ data, but also LSST-like mocks with $1\%$ measurements and we find that the conjoined method of constraining the expansion history and cosmic growth simultaneously is able not only to place stringent constraints on these parameters but also to provide an easy visual way to discriminate cosmological models. Finally, we confirm the existence of a tension between the growth rate and Planck CMB data and we find that the FoM in the conjoined parameter space of $H(z)-f\sigma_8(z)$ can be used to discriminate between the $\Lambda$CDM model and certain classes of modified gravity models, namely the DGP and $f(T)$.

We construct an updated and extended compilation of growth rate data consisting of 34 points and including corrections for model dependence. In order to maximize the independence of the datapoints we also construct a subsample of this compilation (`Gold' growth dataset) which consists of 18 datapoints. We test the consistency of this dataset with the best fit Planck15/$\Lambda$CDM parameters in the context of General Relativity (GR) using the evolution equation for the growth factor $\delta(a)$ with a $w$CDM background. We find tension at the $\sim 3 \sigma$ level between the best fit parameters $w$ (the dark energy equation of state), $\Omega_{0m}$ (the matter density parameter) and $\sigma_8$ (the matter power spectrum normalization on scales $8h^{-1}$Mpc) and the corresponding Planck15/$\Lambda$CDM parameters. We show that the tension disappears if we allow for evolution of the effective Newton's constant, parametrized as $G_{eff}(a)/G_N = 1 + g_a(1-a)^n-g_a(1-a)^{2n}$ with $n\ge2$ where $g_a$, $n$ are parameters of the model, $a$ is the scale factor and $z = 1/a-1$ is the redshift. This parametrization satisfies three criteria: a. $G_{eff} > 0$, b. Consistency with Big Bang Nucleosynthesis ($G_{eff}(a\ll 1)/G_N=1$), c. Consistency with solar system tests ($G_{eff}(a=1)/G_N=1$ and $G_{eff}'(a=1)/G_N=0$). We show that the best fit form of $G_{eff}(z)$ obtained from the growth data corresponds to weakening gravity at recent redshifts (decreasing function of $z$) and we demonstrate that this behavior is not consistent with any scalar-tensor Lagrangian with a real scalar field. Finally, we use MGCAMB to find the best fit $G_{eff}(z)$ obtained from the Planck CMB power spectrum on large angular scales and show that it is a mildly increasing function of $z$, in $3\sigma$ tension with the corresponding decreasing best fit $G_{eff}(z)$ obtained from the growth data.

The growth index of matter fluctuations is computed for ten distinct accelerating cosmological models and confronted to the latest growth rate data via a two-step process. First, we implement a joint statistical analysis in order to place constraints on the free parameters of all models using solely background data. Second, using the observed growth rate of clustering from various galaxy surveys we test the performance of the current cosmological models at the perturbation level while either marginalizing over $\sigma_8$ or having it as a free parameter. As a result, we find that at a statistical level, i.e. after considering the best-fit $\chi^2$ or the value of the Akaike information criterion, most models are in very good agreement with the growth rate data and are practically indistinguishable from $\Lambda$CDM. However, when we also consider the internal consistency of the models by comparing the theoretically predicted values of $(\gamma_0, \gamma_1)$, i.e. the value of the growth index $\gamma(z)$ and its derivative today, with the best-fit ones, we find that the predictions of three out of ten dark energy models are in mild tension with the best-fit ones when $\sigma_8$ is marginalized over. When $\sigma_8$ is free we find that most models are not only in mild tension, but also predict low values for $\sigma_8$. This could be attributed to either a systematic problem with the growth-rate data or the emergence of new physics at low redshifts, with the latter possibly being related to the well-known issue of the lack of power at small scales. Finally, by utilizing mock data based on an LSST-like survey we show that with future surveys and by using the growth index parameterization, it will be possible to resolve the issue of the low $\sigma_8$ but also the tension between the fitted and theoretically predicted values of $(\gamma_0, \gamma_1)$.

We test the viability of the Higgs-dilaton model (HDM) compared to the evolving dark energy ($w_0 w_a$CDM) model, in which the cosmological constant model $\Lambda$CDM is also nested, by using the latest cosmological data that includes the cosmic microwave background temperature, polarization and lensing data from the \textit{Planck} satellite (2015 data release), the BICEP and Keck Array experiments, the Type Ia supernovae from the JLA catalog, the baryon acoustic oscillations from CMASS, LOWZ and 6dF, the weak lensing data from the CFHTLenS survey and the matter power Spectrum measurements from the SDSS (data release 7). We find that the values of all cosmological parameters allowed by the Higgs-dilaton inflation model are well within the \textit{Planck} satellite (2015 data release) constraints. In particular, we have that $w_0 = -1.0001^{+0.0072}_{-0.0074}$, $w_a = 0.00^{+0.15}_{-0.16}$, $n_s = 0.9693^{+0.0083}_{-0.0082}$, $\alpha_s = -0.001^{+0.013}_{-0.014}$ and $r_{0.05} = 0.0025^{+0.0017}_{-0.0016}$ (95.5\%C.L.). We also place new stringent constraints on the couplings of the Higgs-dilaton model and we find that $\xi_{\chi} < 0.00328$ and $\xi_h / \sqrt{\lambda} = 59200^{+30000}_{-20000}$ (95.5\%C.L.). Furthermore, we report that the HDM is at a slightly better footing than the $w_0 w_a$CDM model, as they both have practically the same chi-square, i.e. $\Delta \chi^2 = \chi^2_{w_0 w_a\mathrm{CDM}}-\chi^2_{\mathrm{HDM}}=0.18$, with the HDM model having two fewer parameters. Finally Bayesian evidence favors equally the two models, with the HDM being preferred by the AIC and DIC information criteria.

Gravitational wave (GW) source counts have been recently shown to be able to test how gravitational radiation propagates with the distance from the source. Here, we extend this formalism to cosmological scales, i.e. the high redshift regime, and we discuss the complications of applying this methodology to high redshift sources. We also allow for models with compactified extra dimensions like in the Kaluza-Klein model. Furthermore, we also consider the case of intermediate redshifts, i.e. $0<z\lesssim1$, where we show it is possible to find an analytical approximation for the source counts $\frac{dN}{d(S/N)}$. This can be done in terms of cosmological parameters, such as the matter density $\Omega_{m,0}$ of the cosmological constant model or the cosmographic parameters for a general dark energy model. Our analysis is as general as possible, but it depends on two important factors: a source model for the black hole binary mergers and the GW source to galaxy bias. This methodology also allows us to obtain the higher order corrections of the source counts in terms of the signal-to-noise $S/N$. We then forecast the sensitivity of future observations in constraining GW physics but also the underlying cosmology by simulating sources distributed over a finite range of signal-to-noise with a number of sources ranging from 10 to 500 sources as expected from future detectors. We find that with 500 events it will be possible to provide constraints on the matter density parameter at present $\Omega_{m,0}$ on the order of a few percent and with the precision growing fast with the number of events. In the case of extra dimensions we find that depending on the degeneracies of the model, with 500 events it may be possible to provide stringent limits on the existence of the extra dimensions if the aforementioned degeneracies can be broken.

We use large-scale cosmological observations to place constraints on the dark-matter pressure, sound speed and viscosity, and infer a limit on the mass of warm-dark-matter particles. Measurements of the cosmic microwave background (CMB) anisotropies constrain the equation of state and sound speed of the dark matter at last scattering at the per mille level. Since the redshifting of collisionless particles universally implies that these quantities scale like $a^{-2}$ absent shell crossing, we infer that today $w_{\rm (DM)}< 10^{-10.0}$, $c_{\rm s,(DM)}^2 < 10^{-10.7}$ and $c_{\rm vis, (DM)}^{2} < 10^{-10.3}$ at the $99\%$ confidence level. This very general bound can be translated to model-dependent constraints on dark-matter models: for warm dark matter these constraints imply $m> 70$ eV, assuming it decoupled while relativistic around the same time as the neutrinos; for a cold relic, we show that $m>100$ eV. We separately constrain the properties of the DM fluid on linear scales at late times, and find upper bounds $c_{\rm s, (DM)}^2<10^{-5.9}$, $c_{\rm vis, (DM)}^{2} < 10^{-5.7}$, with no detection of non-dust properties for the DM.

Dark energy equation of state $w(z)$ parametrizations with two parameters and given monotonicity are generically either convex or concave functions. This makes them suitable for fitting either freezing or thawing quintessence models but not both simultaneously. Fitting a dataset based on a freezing model with an unsuitable (concave when increasing) $w(z)$ parametrization (like CPL) can lead to significant misleading features like crossing of the phantom divide line, incorrect $w(z=0)$, incorrect slope \etc that are not present in the underlying cosmological model. To demonstrate this fact we generate scattered cosmological data both at the level of $w(z)$ and the luminosity distance $D_L(z)$ based on either thawing or freezing quintessence models and fit them using parametrizations of convex and of concave type. We then compare statistically significant features of the best fit $w(z)$ with actual features of the underlying model. We thus verify that the use of unsuitable parametrizations can lead to misleading conclusions. In order to avoid these problems it is important to either use both convex and concave parametrizations and select the one with the best $\chi^2$ or use principal component analysis thus splitting the redshift range into independent bins. In the latter case however, significant information about the slope of $w(z)$ at high redshifts is lost. Finally, we propose a new family of parametrizations (nCPL) $w(z)=w_0+w_a (\frac{z}{1+z})^n$ which generalizes the CPL and interpolates between thawing and freezing parametrizations as the parameter $n$ increases to values larger than 1.

The latest Planck results reconfirm the existence of a slight but chronic tension between the best-fit Cosmic Microwave Background (CMB) and low-redshift observables: power seems to be consistently lacking in the late universe across a range of observables (e.g.~weak lensing, cluster counts). We propose a two-parameter model for dark energy where the dark energy is sufficiently like dark matter at large scales to keep the CMB unchanged but where it does not cluster at small scales, preventing concordance collapse and erasing power. We thus exploit the generic scale-dependence of dark energy instead of the more usual time-dependence to address the tension in the data. The combination of CMB, distance and weak lensing data somewhat prefer our model to $\Lambda$CDM, at $\Delta\chi^2=2.4$. Moreover, this improved solution has $\sigma_8=0.79 \pm 0.02$, consistent with the value implied by cluster counts.

We present the analytical solutions for the evolution of matter density perturbations, for a model with a constant dark energy equation of state $w$ but when the effects of the dark energy perturbations are properly taken into account. We consider two cases, the first when the sound speed of the perturbations is zero $c_s^2=0$ and the general case $0<c_s^2 \leq 1$. In the first case our solution is exact, while in the second case we found an approximate solution which works to better than $0.3\%$ accuracy for $k>10 H_0$ or equivalently $k/h>0.0033 \textrm{Mpc}^{-1}$. We also estimate the corrections to the growth index $\gamma(z)$, commonly used to parametrize the growth-rate. We find that these corrections due to the DE perturbations affect the growth index $\gamma$ at the $3\%$ level. We also compare our new expressions for the growth index with other expressions already present in the literature and we find that the latter are less accurate than the ones we propose here. Therefore, our analytical calculations are necessary as the theoretical predictions for the fundamental parameters to be constrained by the upcoming surveys need to be as accurate as possible, especially since we are entering in the precise cosmology era where parameters will be measured to the percent level.

We systematically study the null-test for the growth rate data first presented in [S. Nesseris and D. Sapone, arXiv:1409.3697] and we reconstruct it using various combinations of data sets, such as the $f\sigma_8$ and $H(z)$ or Type Ia supernovae (SnIa) data. We perform the reconstruction in two different ways, either by directly binning the data or by fitting various dark energy models. We also examine how well the null-test can be reconstructed by future data by creating mock catalogs based on the cosmological constant model, a model with strong dark energy perturbations, the $f(R)$ and $f(G)$ models, and the large void LTB model that exhibit different evolution of the matter perturbations. We find that with future data similar to an LSST-like survey, the null-test will be able to successfully discriminate between these different cases at the $5\sigma$ level.

We compare four different methods that can be used to analyze the type Ia supernovae (SnIa) data, ie to use piecewise-constant functions in terms of: the dark energy equation of state $w(z)$, the deceleration parameter $q(z)$, the Hubble parameter $H(z)$ and finally the luminosity distance $d_L$. These four quantities cover all aspects of the accelerating Universe, ie the phenomenological properties of dark energy, the expansion rate (first and second derivatives) of the Universe and the observations themselves. For the first two cases we also perform principal component analysis (PCA) so as to decorrelate the parameters, while for the last two cases we use novel analytic expressions to find the best-fit parameters. In order to test the methods we create mock SnIa data (2000 points, uniform in redshift $z\in[0,1.5]$) for three fiducial cosmologies: the cosmological constant model ($\Lambda$CDM), a linear expansion of the dark energy equation of state parameter $w(a)=w_0+w_a(1-a)$ and the Hu-Sawicki $f(R)$ model. We find that if we focus on the two mainstream approaches for the PCA, i.e. $w(z)$ and $q(z)$, then the best piecewise-constant scheme is always $w(z)$. Finally, to our knowledge the piecewise-constant method for $H(z)$ is new in the literature, while for the rest three methods we present several new analytic expressions.

We test the FLRW cosmology by reconstructing in a model-independent way both the Hubble parameter $H(z)$ and the comoving distance $D(z)$ via the most recent Hubble and Supernovae Ia data. In particular we use: data binning with direct error propagation, the principal component analysis, the genetic algorithms and the Pad\'e approximation. Using our reconstructions we evaluate the Clarkson {\it et al} test known as $\Omega_K(z)$, whose value is constant in redshift for the standard cosmological model, but deviates elsewise. We find good agreement with the expected values of the standard cosmological model within the experimental errors. Finally, we provide forecasts, exploiting the Baryon Acoustic Oscillations measurements from the Euclid survey.

Current and upcoming surveys will measure the cosmological parameters with an extremely high accuracy. The primary goal of these observations is to eliminate some of the currently viable cosmological models created to explain the late time accelerated expansion (either real or only inferred). However, most of the statistical tests used in cosmology have a strong requirement: the use of a model to fit the data. Recently there has been an increased interest on finding tests that are model independent, i.e. to have a function that depends entirely on observed quantities and not on the model, see for instance [1]. In this letter we present an alternative consistency check at the perturbative level for a homogeneous and isotropic Universe filled with a dark energy component. This test makes use of the growth of matter perturbations data and it is able to not only test the homogeneous and isotropic Universe but also, within the framework of a Friedmann-Lema\^itre-Robertson-Walker Universe, if the dark energy component is able to cluster, if there is a tension in the data or if we are dealing with a modification of gravity.

Nonlocal massive gravity can provide an interesting explanation for the late-time cosmic acceleration, with a dark energy equation of state $w_{\rm DE}$ smaller than $-1$ in the past. We derive the equations of linear cosmological perturbations to confront such models with the observations of large-scale structures. The effective gravitational coupling to nonrelativistic matter associated with galaxy clusterings is close to Newton's gravitational constant $G$ for a mass scale $m$ slightly smaller than today's Hubble parameter $H_0$. Taking into account the background expansion history as well as the evolution of matter perturbations $\delta_m$, we test for these models with Type Ia Supernovae (SnIa) from Union 2.1, the cosmic microwave background (CMB) measurements from Planck, a collection of baryon acoustic oscillations (BAO), and the growth rate data of $\delta_m$. Using a higher value of $H_0$ derived from its direct measurement ($H_0 \gtrsim 70$ km s$^{-1}$ Mpc$^{-1}$) the data strongly support the nonlocal massive gravity model ($-1.1 \lesssim w_{\rm DE} \lesssim -1.04$ in the past) over the $\Lambda$CDM model ($w_{\rm DE}=-1$), whereas for a lower prior (67 km s$^{-1}$ Mpc$^{-1}$ $\lesssim$ $H_0 \lesssim 70$ km s$^{-1}$ Mpc$^{-1}$) the two models are statistically comparable.

We investigate the cosmological predictions of several $f(T)$ models, with up to two parameters, at both the background and the perturbation levels. Using current cosmological observations (geometric supernovae type Ia, cosmic microwave background and baryonic acoustic oscillation and dynamical growth data) we impose constraints on the distortion parameter, which quantifies the deviation of these models from the concordance $\Lambda$ cosmology at the background level. In addition we constrain the growth index $\gamma$ predicted in the context of these models using the latest perturbation growth data in the context of three parametrizations for $\gamma$. The evolution of the best fit effective Newton constant, which incorporates the $f(T)$-gravity effects, is also obtained along with the corresponding $1\sigma$ error regions. We show that all the viable parameter sectors of the $f(T)$ gravity models considered practically reduce these models to $\Lambda$CDM. Thus, the degrees of freedom that open up to $\Lambda$CDM in the context of $f(T)$ gravity models are not utilized by the cosmological data leading to an overall disfavor of these models.

It has been long known that by appropriately modifying gravity one can always reproduce the expansion history of any dark energy model, e.g. the $\Lambda$CDM. This degeneracy cannot be broken with geometric probes like the Type Ia supernovae, the cosmic microwave background or the baryon acoustic oscillations since they are based on the measurement of distances and scales and therefore require the expansion history $H(z)$, so one may only hope to break this degeneracy by using dynamic probes like the growth rate data that track the matter density perturbations. We demonstrate that breaking this degeneracy is not currently possible by explicitly constructing $f(R)$, $f(G)$ and $f(T)$ theories that mimic exactly the $\Lambda$CDM model at the background level and confronting them against the latest observational cosmological data. We also determine the necessary improvement in the growth rate data in order to discriminate these theories from $\Lambda$CDM. We found that at least a threefold improvement is necessary, something which will be possible with a survey like Euclid.

We make a comparative analysis of the various independent methods proposed in the literature for studying the nature of dark energy, using four different mocks of SnIa data. In particular, we explore a generic principal components analysis approach, the genetic algorithms, a series of approximations like Pad\'e power law approximants, and various expansions in orthogonal polynomials, as well as cosmography, and compare them with the usual fit to a model with a constant dark energy equation of state w. We find that, depending on the mock data, some methods are more efficient than others at distinguishing the underlying model, although there is no universally better method.

We are entering an era where progress in cosmology is driven by data, and alternative models will have to be compared and ruled out according to some consistent criterium. The most conservative and widely used approach is Bayesian model comparison. In this paper we explicitly calculate the Bayes factors for all models that are linear with respect to their parameters. We do this in order to test the so called Jeffreys' scale and determine analytically how accurate its predictions are in a simple case where we fully understand and can calculate everything analytically. We also discuss the case of nested models, e.g. one with $M_1$ and another with $M_2\supset M_1$ parameters and we derive analytic expressions for both the Bayes factor and the Figure of Merit, defined as the inverse area of the model parameter's confidence contours. With all this machinery and the use of an explicit example we demonstrate that the threshold nature of Jeffreys' scale is not a "one size fits all" reliable tool for model comparison and that it may lead to biased conclusions. Furthermore, we discuss the importance of choosing the right basis in the context of models that are linear with respect to their parameters and how that basis affects the parameter estimation and the derived constraints.

We demonstrate that a wide range of viable f(R) parameterizations (including the Hu & Sawicki and the Starobinsky models) can be expressed as perturbations deviating from the LCDM Lagrangian. We constrain the deviation parameter b using a combination of geometrical and dynamical observational probes. In particular, we perform a joint likelihood analysis of the recent Supernovae Type Ia data, the Cosmic Microwave Background shift parameters, the Baryonic Acoustic Oscillations and the growth rate data provided by the various galaxy surveys. This analysis provides constraints for the following parameters: the matter density Omega_{m0}, the deviation from LCDM parameter b and the growth index gamma(z). We parametrize the growth index gamma(z) in three manners (constant, Taylor expansion around z=0, and Taylor expansion around the scale factor). We point out the numerical difficulty for solving the generalized f(R) Friedman equation at high redshifts due to stiffness of the resulting ordinary differential equation. We resolve this problem by constructing an efficient analytical perturbative method in the deviation parameter b. We demonstrate that this method is highly accurate, by comparing the resulting analytical expressions for the Hubble parameter, with the numerical solutions at low and intermediate redshifts. Surprisingly, despite of its perturbative nature, the accuracy of the method persists even for values of b that are of O(1).

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 $\om_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 $\om_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 \equiv \frac{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\sim 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.

We place observational constraints on models with the late-time cosmic acceleration based on a number of parametrizations allowing fast transitions for the equation of state of dark energy. In addition to the model of Linder and Huterer where the dark energy equation of state $w$ monotonically grows or decreases in time, we propose two new parametrizations in which $w$ has an extremum. We carry out the likelihood analysis with the three parametrizations by using the observational data of supernovae type Ia, cosmic microwave background, and baryon acoustic oscillations. Although the transient cosmic acceleration models with fast transitions can give rise to the total chi square smaller than that in the $\Lambda$-Cold-Dark-Matter ($\Lambda$CDM) model, these models are not favored over $\Lambda$CDM when one uses the Akaike information criterion which penalizes the extra degrees of freedom present in the parametrizations.

We constrain the evolution of Newton's constant using the growth rate of large-scale structure measured by the WiggleZ Dark Energy Survey in the redshift range $0.1 < z < 0.9$. We use this data in two ways. Firstly we constrain the matter density of the Universe, $\omms$ (assuming General Relativity), and use this to construct a diagnostic to detect the presence of an evolving Newton's constant. Secondly we directly measure the evolution of Newton's constant, $\Geff$, that appears in Modified Gravity theories, without assuming General Relativity to be true. The novelty of these approaches are that, contrary to other methods, they do not require knowledge of the expansion history of the Universe, $H(z)$, making them model independent tests. Our constraints for the second derivative of Newton's constant at the present day, assuming it is slowly evolving as suggested by Big Bang Nucleosynthesis constraints, using the WiggleZ data is $\ddotGeff(t_0)=-1.19\pm 0.95\cdot 10^{-20}h^2 \textrm{yr}^{-2}$, where $h$ is defined via $H_0=100 h km s^{-1}Mpc^{-1}$, while using both the WiggleZ and the Sloan Digital Sky Survey Luminous Red Galaxy (SDSS LRG) data is $\ddotGeff(t_0)=-3.6\pm 6.8\cdot 10^{-21}h^2 \rm{yr}^{-2}$, both being consistent with General Relativity. Finally, our constraint for the rms mass fluctuation $\sigma_8$ using the WiggleZ data is $\sigma_8=0.75 \pm 0.08$, while using both the WiggleZ and the SDSS LRG data $\sigma_8=0.77 \pm 0.07$, both in good agreement with the latest measurements from the Cosmic Microwave Background radiation.

We study the cosmology of a covariant Galileon field with five covariant Lagrangians and confront this theory with the most recent cosmological probes: the type Ia supernovae data (Constitution and Union2 sets), cosmic microwave background (WMAP7) and the baryon acoustic oscillations (SDSS7). In the Galileon cosmology with a late-time de Sitter attractor, there is a tracker that attracts solutions with different initial conditions to a common trajectory. Including the cosmic curvature K, we place observational constraints on two distinct cases: (i) the tracker, and (ii) the generic solutions to the equations of motion. We find that the tracker solution can be consistent with the individual observational data, but it is disfavored by the combined data analysis. The generic solutions fare quite well when a non-zero curvature parameter is taken into account, but the Akaike and Bayesian information criteria show that they are not particularly favored over the LCDM model.

The Genetic Algorithm is a heuristic that can be used to produce model independent solutions to an optimization problem, thus making it ideal for use in cosmology and more specifically in the analysis of type Ia supernovae data. In this work we use the Genetic Algorithms (GA) in order to derive a null test on the spatially flat cosmological constant model $\Lambda$CDM. This is done in two steps: first, we apply the GA to the Constitution SNIa data in order to acquire a model independent reconstruction of the expansion history of the Universe $H(z)$ and second, we use the reconstructed $H(z)$ in conjunction with the Om statistic, which is constant only for the $\Lambda$CDM model, to derive our constraints. We find that while $\Lambda$CDM is consistent with the data at the $2\sigma$ level, some deviations from $\Lambda$CDM model at low redshifts can be accommodated.

We use the Om statistic and the Genetic Algorithms (GA) in order to derive a null test on the spatially flat cosmological constant model $\Lambda$CDM. This is done in two steps: first, we apply the GA to the Constitution SNIa data in order to acquire a model independent reconstruction of the expansion history of the Universe $H(z)$ and second, we use the reconstructed $H(z)$ in conjunction with the Om statistic, which is constant only for the $\Lambda$CDM model, to derive our constraints. We find that while $\Lambda$CDM is consistent with the data at the $2\sigma$ level, some deviations from $\Lambda$CDM model at low redshifts can be accommodated.

We try to address quantitatively the question whether a new mass is needed to fit current supernovae data. For this purpose, we consider an infra-red modification of gravity that does not contain any new mass scale but systematic subleading corrections proportional to the curvature. The modifications are of the same type as the one recently derived by enforcing the "Ultra Strong Equivalence Principle" (USEP) upon a Friedmann-Lemaitre-Robertson-Walker universe in the presence of a scalar field. The distance between two comoving observers is altered by these corrections and the observations at high redshift affected at any time during the cosmic evolution. While the specific values of the parameters predicted by USEP are ruled out, there are regions of parameter space that fit SnIa data very well. This allows an interesting possibility to explain the apparent cosmic acceleration today without introducing either a dark energy component or a new mass scale.

We rank the six latest Type Ia supernova (SnIa) datasets (Constitution (C), Union (U), ESSENCE (Davis) (E), Gold06 (G), SNLS 1yr (S) and SDSS-II (D)) in the context of the Chevalier-Polarski-Linder (CPL) parametrization $w(a)=w_0+w_1 (1-a)$, according to their Figure of Merit (FoM), their consistency with the cosmological constant ($\Lambda$CDM), their consistency with standard rulers (Cosmic Microwave Background (CMB) and Baryon Acoustic Oscillations (BAO)) and their mutual consistency. We find a significant improvement of the FoM (defined as the inverse area of the 95.4% parameter contour) with the number of SnIa of these datasets ((C) highest FoM, (U), (G), (D), (E), (S) lowest FoM). Standard rulers (CMB+BAO) have a better FoM by about a factor of 3, compared to the highest FoM SnIa dataset (C). We also find that the ranking sequence based on consistency with $\Lambda$CDM is identical with the corresponding ranking based on consistency with standard rulers ((S) most consistent, (D), (C), (E), (U), (G) least consistent). The ranking sequence of the datasets however changes when we consider the consistency with an expansion history corresponding to evolving dark energy $(w_0,w_1)=(-1.4,2)$ crossing the phantom divide line $w=-1$ (it is practically reversed to (G), (U), (E), (S), (D), (C)). The SALT2 and MLCS2k2 fitters are also compared and some peculiar features of the SDSS-II dataset when standardized with the MLCS2k2 fitter are pointed out. Finally, we construct a statistic to estimate the internal consistency of a collection of SnIa datasets. We find that even though there is good consistency among most samples taken from the above datasets, this consistency decreases significantly when the Gold06 (G) dataset is included in the sample.

We introduce genetic algorithms as a means to analyze supernovae type Ia data and extract model-independent constraints on the evolution of the Dark Energy equation of state. Specifically, we will give a brief introduction to the genetic algorithms along with some simple examples to illustrate their advantages and finally we will apply them to the supernovae type Ia data. We find that genetic algorithms can lead to results in line with already established parametric and non-parametric reconstruction methods and could be used as a complementary way of treating SnIa data. As a non-parametric method, genetic algorithms provide a model-independent way to analyze data and can minimize bias due to premature choice of a dark energy model.

We consider corrections to the Einstein-Hilbert action which contain both higher order and nonlocal terms. We derive an effective Newtonian gravitational constant applicable at the weak field limit and use the primordial nucleosynthesis (BBN) bound and the local gravity constraints on $G_{eff}$ in order to test the viability of several cases of our general Lagrangian. We will also provide a BBN constrain on the $\Box R$ gravitational correction.

We consider theories with an arbitrary coupling between matter and gravity and obtain the perturbation equation of matter on subhorizon scales. Also, we derive the effective gravitational constant $G_{eff}$ and two parameters $\Sigma$ and $\eta$, which along with the perturbation equation of the matter density are useful to constrain the theory from growth factor and weak lensing observations. Finally, we use a completely solvable toy model which exhibits nontrivial phenomenology to investigate specific features of the theory. We obtain the analytic solution of the modified Friedmann equation for the scale factor $a$ in terms of time $t$ and use the age of the oldest star clusters and the primordial nucleosynthesis bounds in order to constrain the parameters of our toy model.

We verify numerically that in the context of general relativity (GR), flat models which have the same $\Omega_{\rm m}$ and CMB shift parameter $R$ but different $H(a)$ and $w(a)$ also have very similar (within less than 8%) growth of perturbations even though the dark energy density evolution is quite different. This provides a direct connection between geometrical and dynamical tests of dark energy and may be used as a cosmological test of general relativity.

We compare the dark energy model constraints obtained by using recent standard ruler data (Baryon Acoustic Oscillations (BAO) at z=0.2 and z=0.35 and Cosmic Microwave Background (CMB) shift parameters R and l_a) with the corresponding constraints obtained by using recent Type Ia Supernovae (SnIa) standard candle data (ESSENCE+SNLS+HST from Davis et. al.). We find that, even though both classes of data are consistent with LCDM at the 2\sigma level, there is a systematic difference between the two classes of data. In particular, we find that for practically all values of the parameters (\Omega_0m,\Omega_b) in the 2\sigma range of the the 3-year WMAP data (WMAP3) best fit, LCDM is significantly more consistent with the SnIa data than with the CMB+BAO data. For example for (\Omega_0m,\Omega_b)=(0.24,0.042) corresponding to the best fit values of WMAP3, the dark energy equation of state parametrization w(z)=w_0 + w_1 (z/(1+z)) best fit is at a 0.5\sigma distance from LCDM (w_0=-1,w_1=0) using the SnIa data and 1.7\sigma away from LCDM using the CMB+BAO data. There is a similar trend in the earlier data (SNLS vs CMB+BAO at z=0.35). This trend is such that the standard ruler CMB+BAO data show a mild preference for crossing of the phantom divide line w=-1, while the recent SnIa data favor LCDM. Despite of this mild difference in trends, we find no statistically significant evidence for violation of the cosmic distance duality relation \eta \equiv d_L(z)/(d_A(z) (1+z)^2)=1. For example, using a prior of \Omega_0m=0.24, we find \eta=0.95 \pm 0.025 in the redshift range 0<z<2, which is consistent with distance duality at the 2\sigma level.

We have compiled a dataset consisting of 22 datapoints at a redshift range (0.15,3.8) which can be used to constrain the linear perturbation growth rate f=\frac{d\ln\delta}{d\ln a}. Five of these data-points constrain directly the growth rate f through either redshift distortions or change of the power spectrum with redshift. The rest of the datapoints constrain f indirectly through the rms mass fluctuation \sigma_8(z) inferred from Ly-\alpha at various redshifts. Our analysis tests the consistency of the LCDM model and leads to a constraint of the Wang-Steinhardt growth index \gamma (defined from f=\Omega_m^\gamma) as \gamma=0.67^{+0.20}_{-0.17}. This result is clearly consistent at $1\sigma$ with the value \gamma={6/11}=0.55 predicted by LCDM. A first order expansion of the index \gamma in redshift space leads to similar results.We also apply our analysis on a new null test of LCDM which is similar to the one recently proposed by Chiba and Nakamura (arXiv:0708.3877) but does not involve derivatives of the expansion rate $H(z)$. This also leads to the fact that LCDM provides an excellent fit to the current linear growth data.

We consider a braneworld model with bulk-brane energy exchange. This allows for crossing of the w=-1 phantom divide line without introducing phantom energy with quantum instabilities. We use the latest SnIa data included in the Gold06 dataset to provide an estimate of the preferred parameter values of this braneworld model. We use three fitting approaches which provide best fit parameter values and hint towards a bulk energy component that behaves like relativistic matter which is propagating in the bulk and is moving at a speed v along the fifth dimension, while the bulk-brane energy exchange component corresponds to negative pressure and signifies energy flowing from the bulk into the brane. We find that the best fit effective equation of state parameter $w_{eff}$ marginally crosses the phantom divide line w=-1. Thus, we have demonstrated both the ability of this class of braneworld models to provide crossing of the phantom divide and also that cosmological data hint towards natural values for the model parameters.

We consider scalar-tensor theories and reconstruct their potential U(\Phi) and coupling F(\Phi) by demanding a background LCDM cosmology. In particular we impose a background cosmic history H(z) provided by the usual flat LCDM parameterization through the radiation (w_{eff}=1/3), matter (w_{eff}=0) and deSitter (w_{eff}=-1) eras. The cosmological dynamical system which is constrained to obey the LCDM cosmic history presents five critical points in each era, one of which corresponding to the standard General Relativity (GR). In the cases that differ from GR, the reconstructed coupling and potential are of the form F(\Phi)\sim \Phi^2 and U(\Phi)\sim F(\Phi)^m where m is a constant. This class of scalar tensor theories is also theoretically motivated by a completely independent approach: imposing maximal Noether symmetry on the scalar-tensor Lagrangian. This approach provides independently: i) the form of the coupling and the potential as F(\Phi)\sim \Phi^2 and U(\Phi)\sim F(\Phi)^m, ii) a conserved charge related to the potential and the coupling and iii) allows the derivation of exact solutions by first integrals of motion.

We consider f(R) modified gravity theories in the metric variation formalism and attempt to reconstruct the function f(R) by demanding a background LCDM cosmology. In particular we impose the following requirements: a. A background cosmic history H(z) provided by the usual flat LCDM parametrization though the radiation (w_eff=1/3), matter (w_eff=0) and deSitter (w_eff=-1) eras. b. Matter and radiation dominate during the `matter' and `radiation' eras respectively i.e. \Omega_m =1 when w_eff=0 and \Omega_r=1 when w_eff=1/3. We have found that the cosmological dynamical system constrained to obey the LCDM cosmic history has four critical points in each era which correspondingly lead to four forms of f(R). One of them is the usual general relativistic form f(R)=R-2\Lambda. The other three forms in each era, reproduce the LCDM cosmic history but they do not satisfy requirement b. stated above. Only one of these forms (different from general relativity) is found to be an attractor of the dynamical cosmological evolution. It has (\Omega_DE=1, \Omega_r=0, \Omega_m=0) throughout the evolution. Its phase space trajectory is numerically obtained.

The Gold06 SnIa dataset recently released in astro-ph/0611572 consists of five distinct subsets defined by the group or instrument that discovered and analyzed the corresponding data. These subsets are: the SNLS subset (47 SnIa), the HST subset (30 SnIa), the HZSST subset (41 SnIa), the SCP subset (26 SnIa) and the Low Redshift (LR) subset (38 SnIa). These subsets sum up to the 182 SnIa of the Gold06 dataset. We use Monte-Carlo simulations to study the statistical consistency of each one of the above subsets with the full Gold06 dataset. In particular, we compare the best fit $w(z)$ parameters (w_0,w_1) obtained by subtracting each one of the above subsets from the Gold06 dataset (subset truncation), with the corresponding best fit parameters (w^r_0,w^r_1) obtained by subtracting the same number of randomly selected SnIa from the same redshift range of the Gold06 dataset (random truncation). We find that the probability for (w^r_0,w^r_1)=(w_0,w_1) is large for the Gold06 minus SCP (Gold06-SCP) truncation but is less than 5% for the Gold06-SNLS, Gold06-HZSST and Gold06-HST truncations. This result implies that the Gold06 dataset is not statistically homogeneous. By comparing the values of the best fit (w_0,w_1) for each subset truncation we find that the tension among subsets is such that the SNLS and HST subsets are statistically consistent with each other and `pull' towards LCDM (w_0=-1,w_1=0) while the HZSST subset is statistically distinct and strongly `pulls' towards a varying w(z) crossing the line $w=-1$ from below (w_0<-1,w_1>0). We also isolate six SnIa that are mostly responsible for this behavior of the HZSST subset.

We use a low redshift expansion of the cosmological equations of extended (scalar-tensor) quintessence to divide the observable Hubble history parameter space in four sectors: A forbidden sector I where the scalar field of the theory becomes imaginary (the kinetic term becomes negative), a forbidden sector II where the scalar field rolls up (instead of down) its potential, an allowed `freezing' quintessence sector III where the scalar field is currently decelerating down its potential towards freezing and an allowed `thawing' sector IV where the scalar field is currently accelerating down its potential. The dividing lines between the sectors depend sensitively on the time derivatives of the Newton's constant G over powers of the Hubble parameter. For minimally coupled quintessence which appears as a special case for a constant G our results are consistent with previous studies. Observable parameter \chi^2 contours based on current data (SNLS dataset) are also constructed on top of the sectors, for a prior of \Omega_m=0.24. By demanding that the observed 2\sigma \chi^2 parameter contours do not lie entirely in the forbidden sectors we derive stringent constraints on the current second time derivative of Newton's constant G. In particular we find {\ddot G}/G >-1.91 H_0^2=-2 10^{-20}h^2 yrs^{-2} at the 2\sigma level which is complementary to solar system tests which constrain only the first derivative of G as |{\dot G}/G|<10^{-14}yrs^{-1} at 1\sigma.

If the dark energy equation of state parameter w(z) crosses the phantom divide line w=-1 (or equivalently if the expression d(H^2(z))/dz - 3\Omega_m H_0^2 (1+z)^2 changes sign) at recent redshifts, then there are two possible cosmological implications: Either the dark energy consists of multiple components with at least one non-canonical phantom component or general relativity needs to be extended to a more general theory on cosmological scales. The former possibility requires the existence of a phantom component which has been shown to suffer from serious theoretical problems and instabilities. Therefore, the later possibility is the simplest realistic theoretical framework in which such a crossing can be realized. After providing a pedagogical description of various dark energy observational probes, we use a set of such probes (including the Gold SnIa sample, the first year SNLS dataset, the 3-year WMAP CMB shift parameter, the SDSS baryon acoustic oscillations peak (BAO), the X-ray gas mass fraction in clusters and the linear growth rate of perturbations at z=0.15 as obtained from the 2dF galaxy redshift survey) to investigate the priors required for cosmological observations to favor crossing of the phantom divide. We find that a low \Omega_m prior (0.2<\Omega_m <0.25) leads, for most observational probes (except of the SNLS data), to an increased probability (mild trend) for phantom divide crossing. An interesting degeneracy of the ISW effect in the CMB perturbation spectrum is also pointed out.

If Newton's constant G evolves on cosmological timescales as predicted by extended gravity theories then Type Ia supernovae (SnIa) can not be treated as standard candles. The magnitude-redshift datasets however can still be useful. They can be used to simultaneously fit for both H(z) and G(z) (so that local G(z) constraints are also satisfied) in the context of appropriate parametrizations. Here we demonstrate how can this analysis be done by applying it to the Gold SnIa dataset. We compare the derived effective equation of state parameter w(z) at best fit with the corresponding result obtained by neglecting the evolution G(z). We show that even though the results clearly differ from each other, in both cases the best fit w(z) crosses the phantom divide w=-1. We then attempt to reconstruct a scalar tensor theory that predicts the derived best fit forms of H(z) and G(z). Since the best fit G(z) fixes the scalar tensor potential evolution F(z), there is no ambiguity in the reconstruction and the potential U(z) can be derived uniquely. The particular reconstructed scalar tensor theory however, involves a change of sign of the kinetic term $\Phi'(z)^2$ as in the minimally coupled case.

We have performed a comparative analysis of three recent and reliable SnIa datasets available in the literature: the Full Gold (FG) dataset (157 data points $0<z<1.75$), a Truncated Gold (TG) dataset (140 data points $0<z<1$) and the most recent Supernova Legacy Survey (SNLS) dataset (115 data points $0<z<1$). We have analyzed and compared the likelihood of cosmological constant and dynamical dark energy parametrizations allowing for crossing of the phantom divide line (PDL). We find that even though the constraints obtained using the three datasets are consistent with each other at the 95% confidence level, the latest (SNLS) dataset shows distinct trends which are not shared by the Gold datasets. We find that the best fit dynamical $w(z)$ obtained from the SNLS dataset does not cross the PDL $w=-1$ and remains above and close to the $w=-1$ line for the whole redshift range $0<z<1$ showing no evidence for phantom behavior. The LCDM parameter values ($w_0=-1$, $w_1=0$) almost coincide with the best fit parameters of the dynamical $w(z)$ parametrizations. In contrast, the best fit dynamical $w(z)$ obtained from the Gold datasets (FG and TG) clearly crosses the PDL and departs significantly from the PDL $w=-1$ line while the LCDM parameter values are about $2\sigma$ away from the best fit $w(z)$. In addition, the $(\Omega_{0m},\Omega_\Lambda)$ parameters in a LCDM parametrization without a flat prior, fit by the SNLS dataset, favor the minimal flat LCDM concordance model. The corresponding fit with the Gold datasets mildly favors a closed universe and the flat LCDM parameter values are $1\sigma - 2\sigma$ away from the best fit $(\Omega_{0m},\Omega_\Lambda)$.

We use the SnIa Gold dataset to compare LCDM with 10 representative parametrizations of the recent Hubble expansion history $H(z)$. For the comparison we use two statistical tests; the usual $\chi_{min}^2$ which is insensitive to the parametrization number of parameters, and a statistic we call the p-test which depends on both the value of $\chi_{min}^2$ and the number $n$ of the parametrization parameters. The p-test measures the confidence level to which the parameter values corresponding to LCDM are excluded from the viewpoint of the parametrization tested. For example, for a linear equation of state parametrization $w(z)=w_0 + w_1 z$ the LCDM parameter values ($w_0=-1$, $w_1=0$) are excluded at 75% confidence level. We use a flat prior and $\Omega_{0m}=0.3$. All parametrizations tested are consistent with the Gold dataset at their best fit. According to both statistical tests, the worst fits among the 10 parametrizations, correspond to the Chaplygin gas, the brane world and the Cardassian parametrizations. The best fit is achieved by oscillating parametrizations which can exclude the parameter values corresponding to LCDM at 85% confidence level. Even though this level of significance does not provide a statistically significant exclusion of LCDM (it is less than $2\sigma$) and does not by itself constitute conclusive evidence for oscillations in the cosmological expansion, when combined with similar independent recent evidence for oscillations coming from the CMB and matter power spectra it becomes an issue worth of further investigation.

We study analytically and numerically the evolution of bound systems in universes with accelerating expansion where the acceleration either increases with time towards a Big Rip singularity (phantom cosmologies) or decreases with time (quintessence). We confirm the finding of Caldwell et. al. [arXiv:astro-ph/0302506] (hereafter CKW) that bound structures get dissociated in phantom cosmologies but we demonstrate that this happens earlier than anticipated by CKW. In particular we find that the `rip time' when a bound system gets unbounded is not the time when the repulsive phantom energy gravitational potential due to the average $(\rho + 3p)$ balances the attractive gravitational potential of the mass M of the system. Instead, the `rip time' is the time when the minimum of the time dependent effective potential (including the centrifugal term) disappears. For the Milky Way galaxy this happens approximately 180Myrs before the Big Rip singularity instead of approximately 60Myrs indicated by CKW for a phantom cosmology with w=-1.5. A numerical reconstruction of the dissociating bound orbits is presented.

We study the expansion history of the universe up to a redshift of z=1.75 using the 194 recently published SnIa data by Tonry et. al. and Barris et. al. In particular we find the best fit forms of several cosmological models and $H(z)$ ansatze, determine the best fit values of their parameters and rank them according to increasing value of $\chi_{min}^2$ (the minimum value of $\chi^2$ for each $H(z)$ ansatz). We use a prior of $\Omega_{0m} = 0.3$ and assume flat geometry of the universe. No prior assumptions are made about validity of energy conditions. The fitted models are fourteen and include SCDM, LCDM, dark energy with constant equation of state parameter $w$ (quiessence), third order polynomial for $H(1+z)$, Chaplygin gas, Cardassian model, $w(z)=w_0 + w_1 z$, $w(z)=w_0 + z w_1/(1+z)$, an oscillating ansatz for $H(z)$ etc. All these models with the exception of SCDM are consistent with the present data. However, the quality of the fit differs significantly among them and so do the predicted forms of $w(z)$ and $H(z)$ at best fit. The worst fit among the data-consistent models considered corresponds to the simplest model LCDM ($\chi_{min}^2 = 198.7$ for $\Omega_{0m} = 0.34$) while the best fit is achieved by the three parameter oscillating ansatz ($\chi_{min}^2 = 194.1$). Most of the best fit ansatze have an equation of state parameter $w(z)$ that varies between $w(z) \simeq -1$ for $z<0.5$ to $w(z) > 0$ for $z>1$. This implies that the sign of the pressure of the dark energy may be alternating as the redshift increases. The goodness of fit of the oscillating $H(z)$ ansatz lends further support to this possibility.