The Complete Calibration of the Colour-Redshift Relation survey (C3R2) is a spectroscopic effort involving ESO and Keck facilities designed to empirically calibrate the galaxy colour-redshift relation - P(z|C) to the Euclid depth (i_AB=24.5) and is intimately linked to upcoming Stage IV dark energy missions based on weak lensing cosmology. The aim is to build a spectroscopic calibration sample that is as representative as possible of the galaxies of the Euclid weak lensing sample. In order to minimise the number of spectroscopic observations to fill the gaps in current knowledge of the P(z|C), self-organising map (SOM) representations of the galaxy colour space have been constructed. Here we present the first results of an ESO@ VLT Large Programme approved in the context of C3R2, which makes use of the two VLT optical and near-infrared multi-object spectrographs, FORS2 and KMOS. This paper focuses on high-quality spectroscopic redshifts of high-z galaxies observed with the KMOS spectrograph in the H- and K-bands. A total of 424 highly-reliable z are measured in the 1.3<=z<=2.5 range, with total success rates of 60.7% in the H-band and 32.8% in the K-band. The newly determined z fill 55% of high and 35% of lower priority empty SOM grid cells. We measured Halpha fluxes in a 1."2 radius aperture from the spectra of the spectroscopically confirmed galaxies and converted them into star formation rates. In addition, we performed an SED fitting analysis on the same sample in order to derive stellar masses, E(B-V), total magnitudes, and SFRs. We combine the results obtained from the spectra with those derived via SED fitting, and we show that the spectroscopic failures come from either weakly star-forming galaxies (at z<1.7, i.e. in the H-band) or low S/N spectra (in the K-band) of z>2 galaxies.

In recent years forecasting activities have become a very important tool for designing and optimising large scale structure surveys. To predict the performance of such surveys, the Fisher matrix formalism is frequently used as a fast and easy way to compute constraints on cosmological parameters. Among them lies the study of the properties of dark energy which is one of the main goals in modern cosmology. As so, a metric for the power of a survey to constrain dark energy is provided by the Figure of merit (FoM). This is defined as the inverse of the surface contour given by the joint variance of the dark energy equation of state parameters $\{w_0,w_a\}$ in the Chevallier-Polarski-Linder parameterisation, which can be evaluated from the covariance matrix of the parameters. This covariance matrix is obtained as the inverse of the Fisher matrix. Inversion of an ill-conditioned matrix can result in large errors on the covariance coefficients if the elements of the Fisher matrix have been estimated with insufficient precision. The conditioning number is a metric providing a mathematical lower limit to the required precision for a reliable inversion, but it is often too stringent in practice for Fisher matrices with size larger than $2\times2$. In this paper we propose a general numerical method to guarantee a certain precision on the inferred constraints, like the FoM. It consists on randomly vibrating (perturbing) the Fisher matrix elements with Gaussian perturbations of a given amplitude, and then evaluating the maximum amplitude that keeps the FoM within the chosen precision. The steps used in the numerical derivatives and integrals involved in the calculation of the Fisher matrix elements can then be chosen accordingly in order to keep the precision of the Fisher matrix elements below this maximum amplitude...

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 data from the Euclid mission will enable the measurement of the photometric redshifts, angular positions, and weak lensing shapes for over a billion galaxies. This large dataset will allow for cosmological analyses using the angular clustering of galaxies and cosmic shear. The cross-correlation (XC) between these probes can tighten constraints and it is therefore important to quantify their impact for Euclid. In this study we carefully quantify the impact of XC not only on the final parameter constraints for different cosmological models, but also on the nuisance parameters. In particular, we aim at understanding the amount of additional information that XC can provide for parameters encoding systematic effects, such as galaxy bias or intrinsic alignments (IA). We follow the formalism presented in Euclid Collaboration: Blanchard et al. (2019) and make use of the codes validated therein. We show that XC improves the dark energy Figure of Merit (FoM) by a factor $\sim 5$, whilst it also reduces the uncertainties on galaxy bias by $\sim 17\%$ and the uncertainties on IA by a factor $\sim 4$. We observe that the role of XC on the final parameter constraints is qualitatively the same irrespective of the galaxy bias model used. We also show that XC can help in distinguishing between different IA models, and that if IA terms are neglected then this can lead to significant biases on the cosmological parameters. We find that the XC terms are necessary to extract the full information content from the data in future analyses. They help in better constraining the cosmological model, and lead to a better understanding of the systematic effects that contaminate these probes. Furthermore, we find that XC helps in constraining the mean of the photometric-redshift distributions, but it requires a more precise knowledge of this mean in order not to degrade the final FoM. [Abridged]

The jackknife method gives an internal covariance estimate for large-scale structure surveys and allows model-independent errors on cosmological parameters. Using the SDSS-III BOSS CMASS sample, we study how the jackknife size and number of resamplings impact the precision of the covariance estimate on the correlation function multipoles and the error on the inferred baryon acoustic scale. We compare the measurement with log-normal mock galaxy catalogues with the same survey geometry. We build several jackknife configurations that vary in size and number of resamplings. We find that it is useful to apply the tapering scheme to estimate the precision matrix from a limited number of jackknife resamplings. The results from CMASS and mock catalogues show that the error estimate of the baryon acoustic scale does not depend on the jackknife scale. For the shift parameter $\alpha$, we find an average error of 1.6% and 1.2%, respectively from CMASS and mock jackknife covariances, consistent with pre-reconstruction analyses. However, when relatively few resamplings are used, the jackknife error estimate becomes unreliable. Future large-scale structure surveys will map even greater volumes allowing percent-level estimation of the covariance matrix using a jackknife approach.

Stage IV weak lensing experiments will offer more than an order of magnitude leap in precision. We must therefore ensure that our analyses remain accurate in this new era. Accordingly, previously ignored systematic effects must be addressed. In this work, we evaluate the impact of the reduced shear approximation and magnification bias, on the information obtained from the angular power spectrum. To first-order, the statistics of reduced shear, a combination of shear and convergence, are taken to be equal to those of shear. However, this approximation can induce a bias in the cosmological parameters that can no longer be neglected. A separate bias arises from the statistics of shear being altered by the preferential selection of galaxies and the dilution of their surface densities, in high-magnification regions. The corrections for these systematic effects take similar forms, allowing them to be treated together. We calculated the impact of neglecting these effects on the cosmological parameters that would be determined from Euclid, using cosmic shear tomography. To do so, we employed the Fisher matrix formalism, and included the impact of the super-sample covariance. We also demonstrate how the reduced shear correction can be calculated using a lognormal field forward modelling approach. These effects cause significant biases in Omega_m, sigma_8, n_s, Omega_DE, w_0, and w_a of -0.53 sigma, 0.43 sigma, -0.34 sigma, 1.36 sigma, -0.68 sigma, and 1.21 sigma, respectively. We then show that these lensing biases interact with another systematic: the intrinsic alignment of galaxies. Accordingly, we develop the formalism for an intrinsic alignment-enhanced lensing bias correction. Applying this to Euclid, we find that the additional terms introduced by this correction are sub-dominant.

We study the imprints of an effective dark energy fluid in the large scale structure of the universe through the observed angular power spectrum of galaxies in the relativistic regime. We adopt the phenomenological approach that introduces two parameters $\{Q,\eta\}$ at the level of linear perturbations and allow to take into account the modified clustering (or effective gravitational constant) and anisotropic stress appearing in models beyond $\Lambda$CDM. We characterize the effective dark energy fluid by an equation of state parameter $w=-0.95$ and various sound speed cases in the range $10^{-6}\leq c^2_s\leq 1$, thus covering K-essence and quintessence cosmologies. We calculate the angular power spectra of standard and relativistic effects for these scenarios under the $\{Q,\eta\}$ parametrization, and we compare these relative to a fiducial $\Lambda$CDM cosmology. We find that, overall, deviations relative to $\Lambda$CDM are stronger at low redshift since the behavior of the dark energy fluid can mimic the cosmological constant during matter domination era but departs during dark energy domination. In particular, at $z=0.1$ the matter density fluctuations are suppressed by up to $\sim3\%$ for the quintessence-like case, while redshift-space distortions and Doppler effect can be enhanced by $\sim15\%$ at large scales for the lowest sound speed scenario. On the other hand, at $z=2$ we find deviations of up to $\sim5\%$ in gravitational lensing, whereas the Integrated Sachs-Wolfe effect can deviate up to $\sim17\%$. Furthermore, when considering an imperfect dark energy fluid scenario, we find that all effects are insensitive to the presence of anisotropic stress at low redshift, and only the Integrated Sachs-Wolfe effect can detect this feature at $z=2$ and very large scales.

The Baryon Acoustic Oscillation (BAO) scale in the clustering of galaxies is a powerful standard ruler to measure cosmological distances and determine the geometry of the Universe. Past surveys have detected the BAO feature in the clustering of different galaxy samples, most of them composed of redder, quiescent galaxies and bluer, star-forming ones out to redshift $z\sim1$. Besides these targets, new upcoming surveys will observe high-redshift galaxies with bright nebular emission lines out to $z\sim2$, quasars and Lyman-$\alpha$ quasars at $z\gt 2$. All these different galaxy targets will be used as multi-tracers of the same underlying dark matter field. By combining them over wide cosmological volumes, we will be able to beat cosmic variance and measure the growth of structure with unprecedented accuracy. In this work, we measure the BAO scale in the two-point auto- and cross-correlation functions of three independent populations of multi-tracers extracted from the SDSS DR7 Main galaxy sample at redshift $0.02\lt z\lt 0.22$. Combining their covariances, we find accurate constraints on the shift parameter $\alpha=1.00\pm 0.04$ and $D_{\rm{V}}(z=0.1)/r_{\rm{s}}=2.92\pm0.12$.

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 predictions of the yield of $7<z<9$ quasars from the Euclid wide survey, updating the calculation presented in the Euclid Red Book in several ways. We account for revisions to the Euclid near-infrared filter wavelengths; we adopt steeper rates of decline of the quasar luminosity function (QLF; $\Phi$) with redshift, $\Phi\propto10^{k(z-6)}$, $k=-0.72$, and a further steeper rate of decline, $k=-0.92$; we use better models of the contaminating populations (MLT dwarfs and compact early-type galaxies); and we use an improved Bayesian selection method, compared to the colour cuts used for the Red Book calculation, allowing the identification of fainter quasars, down to $J_{AB}\sim23$. Quasars at $z>8$ may be selected from Euclid $OYJH$ photometry alone, but selection over the redshift interval $7<z<8$ is greatly improved by the addition of $z$-band data from, e.g., Pan-STARRS and LSST. We calculate predicted quasar yields for the assumed values of the rate of decline of the QLF beyond $z=6$. For the case that the decline of the QLF accelerates beyond $z=6$, with $k=-0.92$, Euclid should nevertheless find over 100 quasars with $7.0<z<7.5$, and $\sim25$ quasars beyond the current record of $z=7.5$, including $\sim8$ beyond $z=8.0$. The first Euclid quasars at $z>7.5$ should be found in the DR1 data release, expected in 2024. It will be possible to determine the bright-end slope of the QLF, $7<z<8$, $M_{1450}<-25$, using 8m class telescopes to confirm candidates, but follow-up with JWST or E-ELT will be required to measure the faint-end slope. Contamination of the candidate lists is predicted to be modest even at $J_{AB}\sim23$. The precision with which $k$ can be determined over $7<z<8$ depends on the value of $k$, but assuming $k=-0.72$ it can be measured to a 1 sigma uncertainty of 0.07.

Our aim is to quantify the impact of systematic effects on the inference of cosmological parameters from cosmic shear. We present an end-to-end approach that introduces sources of bias in a modelled weak lensing survey on a galaxy-by-galaxy level. Residual biases are propagated through a pipeline from galaxy properties (one end) through to cosmic shear power spectra and cosmological parameter estimates (the other end), to quantify how imperfect knowledge of the pipeline changes the maximum likelihood values of dark energy parameters. We quantify the impact of an imperfect correction for charge transfer inefficiency (CTI) and modelling uncertainties of the point spread function (PSF) for Euclid, and find that the biases introduced can be corrected to acceptable levels.

The Euclid space telescope will measure the shapes and redshifts of galaxies to reconstruct the expansion history of the Universe and the growth of cosmic structures. Estimation of the expected performance of the experiment, in terms of predicted constraints on cosmological parameters, has so far relied on different methodologies and numerical implementations, developed for different observational probes and for their combination. In this paper we present validated forecasts, that combine both theoretical and observational expertise for different cosmological probes. This is presented to provide the community with reliable numerical codes and methods for Euclid cosmological forecasts. We describe in detail the methodology adopted for Fisher matrix forecasts, applied to galaxy clustering, weak lensing and their combination. We estimate the required accuracy for Euclid forecasts and outline a methodology for their development. We then compare and improve different numerical implementations, reaching uncertainties on the errors of cosmological parameters that are less than the required precision in all cases. Furthermore, we provide details on the validated implementations that can be used by the reader to validate their own codes if required. We present new cosmological forecasts for Euclid. We find that results depend on the specific cosmological model and remaining freedom in each setup, i.e. flat or non-flat spatial cosmologies, or different cuts at nonlinear scales. The validated numerical implementations can now be reliably used for any setup. We present results for an optimistic and a pessimistic choice of such settings. We demonstrate that the impact of cross-correlations is particularly relevant for models beyond a cosmological constant and may allow us to increase the dark energy Figure of Merit by at least a factor of three.

The extended Jordan-Brans-Dicke (eJBD) theory of gravity is constrained by a host of astrophysical and cosmological observations spanning a wide range of scales. The current cosmological constraints on the first post-Newtonian parameter in these simplest eJBD models in which the recent acceleration of the Universe is connected with the variation of the effective gravitational strength are consistent, but approximately two order of magnitude larger than the time-delay test within the Solar System. We forecast the capabilities of future galaxy surveys in combination with current and future CMB anisotropies measurements to further constrain the simplest dark energy models within eJBD theory of gravity. By considering two cases of a monomial potential (a quartic potential or a cosmological constant), we show how Euclid-like galaxy clustering and weak lensing data in combination with BOSS and future CMB observations have the potential to reach constraints on the first post-Newtonian parameter $\gamma_{\rm PN}$ comparable to those from the Solar System.

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.

In this work we perform an analysis on the recently proposed conjoined cosmic growth and cosmic expansion diagram [1] to compare several dark energy models using the Figure of Merit showed in [2], which consists in the inverse of the $1\sigma$ confidence region in the $f\sigma_8(z)-H(z)$ plot. Our analysis also consists of comparing the models by performing different statistical criteria: Bayes factor [3], the Bayesian Information Criteria [4] and the Akaike Information Criterion [5]. We also developed a 3-dimensional Figure of Merit to account simultaneously for the errors on the growth rate and the Hubble parameter. The main idea is to consider several cosmological models and compare them with the different statistical criteria in order to highlight the differences and the accuracies of each single criterion.

We analyse the efficiency of future large scale structure surveys to unveil the presence of scale dependent features in the primordial spectrum --resulting from cosmic inflation-- imprinted in the distribution of galaxies. Features may appear as a consequence of non-trivial dynamics during cosmic inflation, in which one or more background quantities experienced small but rapid deviations from their characteristic slow-roll evolution. We consider two families of features: localized features and oscillatory extended features. To characterise them we employ various possible templates parametrising their scale dependence and provide forecasts on the constraints on these parametrisations for LSST like surveys. We perform a Fisher matrix analysis for three observables: cosmic microwave background (CMB), galaxy clustering and weak lensing. We find that the combined data set of these observables will be able to limit the presence of features down to levels that are more restrictive than current constraints coming from CMB observations only. In particular, we address the possibility of gaining information on currently known deviations from scale invariance inferred from CMB data, such as the feature appearing at the $\ell \sim 20$ multipole (which is the main contribution to the low-$\ell$ deficit) and a potential feature appearing at $\ell \sim 800$.

We test both the FLRW geometry and $\Lambda$CDM cosmology in a model independent way by reconstructing the Hubble function $H(z)$, the comoving distance $D(z)$ and the growth of structure $f\sigma_8(z)$ using the most recent data available. We use the linear model formalism in order to optimally reconstruct the latter cosmological functions, together with their derivatives and integrals. We then evaluate four of the null tests available in literature: $Om_{1}$ by Sahni et al., $Om_{2}$ by Zunckel \& Clarkson, $Ok$ by Clarkson et al., and $ns$ by Nesseris \& Sapone. For all the four tests we find agreement, within the errors, with the standard cosmological model.

Euclid is a European Space Agency medium class mission selected for launch in 2020 within the Cosmic Vision 2015 2025 program. The main goal of Euclid is to understand the origin of the accelerated expansion of the universe. Euclid will explore the expansion history of the universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

Euclid is a European Space Agency medium class mission selected for launch in 2019 within the Cosmic Vision 2015-2025 programme. The main goal of Euclid is to understand the origin of the accelerated expansion of the Universe. Euclid will explore the expansion history of the Universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

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.

Cosmological fluids are commonly assumed to be distributed in a spatially homogeneous way, while their internal properties are described by a perfect fluid. As such, they influence the Hubble-expansion through their respective densities and equation of state parameters. The subject of this paper is an investigation of the fluid-mechanical properties of a dark energy fluid, which is characterised by its sound speed and its viscosity apart from its equation of state. In particular, we compute the predicted spectra for the integrated Sachs-Wolfe effect for our generalised fluid, and compare them with the corresponding predictions for weak gravitational lensing and galaxy clustering, which had been computed in previous work. We perform statistical forecasts and show that the integrated Sachs-Wolfe signal obtained by cross correlating Euclid galaxies with Planck temperatures, when joined to galaxy clustering and weak lensing observations, yields a percent sensitivity on the dark energy sound speed and viscosity. We prove that the iSW effect provides strong degeneracy breaking for low sound speeds and large differences between the sound speed and viscosity parameters.

Recent work has demonstrated that it is important to constrain the dynamics of cosmological perturbations, in addition to the evolution of the background, if we want to distinguish among different models of the dark sector. Especially the anisotropic stress of the (possibly effective) dark energy fluid has been shown to be an important discriminator between modified gravity and dark energy models. In this paper we use approximate analytical solutions of the perturbation equations in the presence of viscosity to study how the anisotropic stress affects the weak lensing and galaxy power spectrum. We then forecast how sensitive the photometric and spectroscopic Euclid surveys will be to both the speed of sound and the viscosity of our effective dark energy fluid when using weak lensing tomography and the galaxy power spectrum. We find that Euclid alone can only constrain models with very small speed of sound and viscosity, while it will need the help of other observables in order to give interesting constraints on models with a sound speed close to one. This conclusion is also supported by the expected Bayes factor between models.

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.

The characterisation of dark energy is one of the primary goals in cosmology especially now that many new experiments are being planned with the aim of reaching a high sensitivity on cosmological parameters. It is known that if we move away from the simple cosmological constant model then we need to consider perturbations in the dark energy fluid. This means that dark energy has two extra degrees of freedom: the sound speed $\cs$ and the anisotropic stress $\sigma$. If dark energy is inhomogenous at the scales of interest then the gravitational potentials are modified and the evolution of the dark matter perturbations is also directly affected. In this paper we add an anisotropic component to the dark energy perturbations. Following the idea introduced in \cite{Sapone:2009mb}, we solve analytically the equations of perturbations in the dark sector, finding simple and accurate approximated solutions. We also find that the evolution of the density perturbations is governed by an effective sound speed which depends on both the sound speed and the anisotropic stress parameter. We then use these solutions to look at the impact of the dark energy perturbations on the matter power spectrum and on the Integrated Sachs-Wolfe effect in the Cosmic Microwave Background.

We provide exact solutions to the cosmological matter perturbation equation in a homogeneous FLRW universe with a vacuum energy that can be parametrized by a constant equation of state parameter $w$ and a very accurate approximation for the Ansatz $w(a)=w_0+w_a(1-a)$. We compute the growth index $\gamma=\log f(a)/\log\Om_m(a)$, and its redshift dependence, using the exact and approximate solutions in terms of Legendre polynomials and show that it can be parametrized as $\gamma(a)=\gamma_0+\gamma_a(1-a)$ in most cases. We then compare four different types of dark energy (DE) models: $w\Lambda$CDM, DGP, $f(R)$ and a LTB-large-void model, which have very different behaviors at $z\gsim1$. This allows us to study the possibility to differentiate between different DE alternatives using wide and deep surveys like Euclid, which will measure both photometric and spectroscopic redshifts for several hundreds of millions of galaxies up to redshift $z\simeq 2$. We do a Fisher matrix analysis for the prospects of differentiating among the different DE models in terms of the growth index, taken as a given function of redshift or with a principal component analysis, with a value for each redshift bin for a Euclid-like survey. We use as observables the complete and marginalized power spectrum of galaxies $P(k)$ and the Weak Lensing (WL) power spectrum. We find that, using $P(k)$, one can reach (2%, 5%) errors in $(w_0, w_a)$, and (4%, 12%) errors in $(\gamma_0, \gamma_a)$, while using WL we get errors at least twice as large. These estimates allow us to differentiate easily between DGP, $f(R)$ models and $\Lambda$CDM, while it would be more difficult to distinguish the latter from a variable equation of state parameter or LTB models using only the growth index.}

In this paper we review a part of the approaches that have been considered to explain the extraordinary discovery of the late time acceleration of the Universe. We discuss the arguments that have led physicists and astronomers to accept dark energy as the current preferable candidate to explain the acceleration. We highlight the problems and the attempts to overcome the difficulties related to such a component. We also consider alternative theories capable of explaining the acceleration of the Universe, such as modification of gravity. We compare the two approaches and point out the observational consequences, reaching the sad but foresightful conclusion that we will not be able to distinguish between a Universe filled by dark energy or a Universe where gravity is different from General Relativity. We review the present observations and discuss the future experiments that will help us to learn more about our Universe. This is not intended to be a complete list of all the dark energy models but this paper should be seen as a review on the phenomena responsible for the acceleration. Moreover, in a landscape of hardly compelling theories, it is an important task to build simple measurable parameters useful for future experiments that will help us to understand more about the evolution of the Universe.

The characterization of dark energy is a central task of cosmology. To go beyond a cosmological constant, we need to introduce at least an equation of state and a sound speed and consider observational tests that involve perturbations. If dark energy is not completely homogeneous on observable scales then the Poisson equation is modified and dark matter clustering is directly affected. One can then search for observational effects of dark energy clustering using dark matter as a probe. In this paper we exploit an analytical approximate solution of the perturbation equations in a general dark energy cosmology to analyze the performance of next-decade large scale surveys in constraining equation of state and sound speed. We find that tomographic weak lensing and galaxy redshift surveys can constrain the sound speed of the dark energy only if the latter is small, of the order of $c_{s}\lesssim0.01$ (in units of $c$). For larger sound speeds the error grows to 100% and more. We conclude that large scale structure observations contain very little information about the perturbations in canonical scalar field models with a sound speed of unity. Nevertheless, they are able to detect the presence of "cold" dark energy, i.e. a dark energy with non-relativistic speed of sound.

Dark energy perturbations are normally either neglected or else included in a purely numerical way, obscuring their dependence on underlying parameters like the equation of state or the sound speed. However, while many different explanations for the dark energy can have the same equation of state, they usually differ in their perturbations so that these provide a fingerprint for distinguishing between different models with the same equation of state. In this paper we derive simple yet accurate approximations that are able to characterize a specific class of models (encompassing most scalar-field models) which is often generically called "dark energy". We then use the approximate solutions to look at the impact of the dark energy perturbations on the dark matter power spectrum and on the integrated Sachs-Wolfe effect in the cosmic microwave background radiation.

Dark energy can be studied by its influence on the expansion of the Universe as well as on the growth history of the large-scale structure. In this paper, we follow the growth of the cosmic density field in early dark energy cosmologies by combining observations of the primary CMB temperature and polarisation power spectra at high redshift, of the CMB lensing deflection field at intermediate redshift and of weak cosmic shear at low redshifts for constraining the allowed amount of early dark energy. We present these forecasts using the Fisher-matrix formalism and consider the combination of Planck-data with the weak lensing survey of Euclid. We find that combining these data sets gives powerful constraints on early dark energy and is able to break degeneracies in the parameter set inherent to the various observational channels. The derived statistical 1-sigma-bound on the early dark energy density parameter is sigma(Omega_d^e)=0.0022 which suggests that early dark energy models can be well examined in our approach. In addition, we derive the dark energy figure of merit for the considered dark energy parameterisation and comment on the applicability of the growth index to early dark energy cosmologies.

We discuss the phenomenology of the dark energy in first order perturbation theory, demonstrating that the dark energy cannot be fully constrained unless the dark matter is found, and that there are two functions that characterise the observational properties of the dark sector for cosmological probes. We argue that measuring these two functions should be an important goal for observational cosmology in the next decades.

The growth factor of linear fluctuations is probably one of the least known quantity in observational cosmology. Here we discuss the constraints that baryon oscillations in galaxy power spectra from future surveys can put on a conveniently parametrized growth factor. We find that spectroscopic surveys of $5000 deg^2$ extending to $z \approx 3$ could estimate the growth index $\gamma$ within 0.06; a similar photometric survey would give $\Delta\gamma\approx 0.15$. This test provides an important consistency check for standard cosmological model and could constrain modified gravity models. We discuss the errors and the figure of merit for various combinations of redshift errors and survey size.

We introduce a convenient parametrization of dark energy models that is general enough to include several modified gravity models and generalized forms of dark energy. In particular we take into account the linear perturbation growth factor, the anisotropic stress and the modified Poisson equation. We discuss the sensitivity of large scale weak lensing surveys like the proposed DUNE satellite to these parameters. We find that a large-scale weak-lensing tomographic survey is able to easily distinguish the Dvali-Gabadadze-Porrati model from LCDM and to determine the perturbation growth index to an absolute error of 0.02-0.03.

We consider fluid perturbations close to the "phantom divide" characterised by p = -rho and discuss the conditions under which divergencies in the perturbations can be avoided. We find that the behaviour of the perturbations depends crucially on the prescription for the pressure perturbation delta-p. The pressure perturbation is usually defined using the dark energy rest-frame, but we show that this frame becomes unphysical at the divide. If the pressure perturbation is kept finite in any other frame, then the phantom divide can be crossed. Our findings are important for generalised fluid dark energy used in data analysis (since current cosmological data sets indicate that the dark energy is characterised by p ~ -rho so that p < -rho cannot be excluded) as well as for any models crossing the phantom divide, like some modified gravity, coupled dark energy and braneworld models. We also illustrate the results by an explicit calculation for the "Quintom" case with two scalar fields.

There is now strong observational evidence that the expansion of the universe is accelerating. The standard explanation invokes an unknown "dark energy" component. But such scenarios are faced with serious theoretical problems, which has led to increased interest in models where instead General Relativity is modified in a way that leads to the observed accelerated expansion. The question then arises whether the two scenarios can be distinguished. Here we show that this may not be so easy, demonstrating explicitely that a generalised dark energy model can match the growth rate of the DGP model and reproduce the 3+1 dimensional metric perturbations. Cosmological observations are then unable to distinguish the two cases.

We estimate the constraints that the recent high-redshift sample of supernovae type Ia put on a phenomenological interaction between dark energy and dark matter. The interaction can be interpreted as arising from the time variation of the mass of dark matter particles. We find that the coupling correlates with the equation of state: roughly speaking, a negative coupling (in our sign convention) implies phantom energy ($w_{\phi}<-1$) while a positive coupling implies ``ordinary'' dark energy. The constraints from the current supernovae Ia Hubble diagram favour a negative coupling and an equation of state $w_{\phi}<-1$. A zero or positive coupling is in fact unlikely at 99% c.l. (assuming constant equation of state); at the same time non-phantom values ($w_{\phi}>-1$) are unlikely at 95%. We show also that the usual bounds on the energy density weaken considerably when the coupling is introduced: values as large as $\Omega_{m0}=0.7$ become acceptable for as concerns SNIa. We find that the rate of change of the mass $\dot{m}/m$ of the dark matter particles is constrained to be $\delta_{0}$ in a Hubble time, with $-10<\delta_{0}<-1$ to 95% c.l.. We show that a large positive coupling might in principle avoid the future singularity known as ``big rip'' (occurring for $w_{\phi}<-1$) but the parameter region for this to occur is almost excluded by the data. We also forecast the constraints that can be obtained from future experiments, focusing on supernovae and baryon oscillations in the power spectra of deep redshift surveys. We show that the method of baryon oscillations holds the best potential to contrain the coupling.