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.

We propose a new tomographic estimator for the gravitational lensing potential, based on a combination of intensity mapping (IM) and galaxy number counts. The estimator can be written schematically as IM$\times$galaxy $-$ galaxy$\times$IM; this combination allows to greatly reduce the contamination by density-density correlations, thus isolating the lensing signal. As a pure cross-correlation estimator, it is additionally less susceptible to systematic effects. We show that the new estimator strongly suppresses cosmic variance and consequently improves the signal-to-noise ratio (SNR) for the detection of lensing, especially on linear scales and intermediate redshifts. For cosmic variance dominated surveys the SNR of our estimator is a factor 30 larger than the SNR obtained from the correlation of galaxy number counts only. Shot noise and interferometer noise reduce the SNR. For the specific example of DES cross-correlated with HIRAX, the SNR is around 4, whereas for Euclid cross-correlated with HIRAX it reaches 52. This corresponds to an improvement of a factor 4-5 compared to the SNR from DES alone. For Euclid cross-correlated with HIRAX the improvement with respect to Euclid alone strongly depends on the redshift. We find that the improvement is particularly important for redshifts below 1.6, where it reaches 5. This makes our estimator especially valuable to test dark energy and modified gravity, that are expected to leave an impact at low and intermediate redshifts.

Future galaxy clustering surveys will probe small scales where non-linearities become important. Since the number of modes accessible on intermediate to small scales is very high, having a precise model at these scales is important especially in the context of discriminating alternative cosmological models from the standard one. In the mildly non-linear regime, such models typically differ from each other, and galaxy clustering data will become very precise on these scales in the near future. As the observable quantity is the angular power spectrum in redshift space, it is important to study the effects of non-linear density and redshift space distortion (RSD) in the angular power spectrum. We compute non-linear contributions to the angular power spectrum using a flat-sky approximation that we introduce in this work, and compare the results of different perturbative approaches with $N$-body simulations. We find that the TNS perturbative approach is significantly closer to the $N$-body result than Eulerian or Lagrangian 1-loop approximations, effective field theory of large scale structure or a halofit-inspired model. However, none of these prescriptions is accurate enough to model the angular power spectrum well into the non-linear regime. In addition, for narrow redshift bins, $\Delta z \lesssim 0.01$, the angular power spectrum acquires non-linear contributions on all scales, right down to $\ell=2$, and is hence not a reliable tool at this time. To overcome this problem, we need to model non-linear RSD terms, for example as TNS does, but for a matter power spectrum that remains reasonably accurate well into the deeply non-linear regime, such as halofit.

In this paper, we study lensing of 21cm intensity mapping (IM). Like in the cosmic microwave background (CMB), there is no first order lensing in intensity mapping. The first effects in the power spectrum are therefore of second and third order. Despite this, lensing of the CMB power spectrum is an important effect that needs to be taken into account, which motivates the study of the impact of lensing on the IM power spectrum. We derive a general formula up to third order in perturbation theory including all the terms with two derivatives of the gravitational potential, i.e. the dominant terms on sub-Hubble scales. We then show that in intensity mapping there is a new lensing term which is not present in the CMB. We obtain that the signal-to-noise of 21 cm lensing for futuristic surveys like SKA2 is about 10. We find that surveys probing only large scales, lmax < 700, can safely neglect the lensing of the intensity mapping power spectrum, but that otherwise this effect should be included.

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.

We explore the backreaction model based on the template metric proposed in Larena et al. (2008) constraining the matter density parameter $\Omega_m^{D_0}$ and the Dark Energy parameter $w$ with recent data. We provide constraints based on Supernovae Ia from the SNLS and the Union2.1 catalogs, confirming that the backreacted Universe should have a higher matter density than the corresponding Friedmaniann one. Angular diameter distances from clusters data confirm the same feature. Finally we combine these results with constraints obtained from the position of the first three peaks and the first dip of the CMB power spectrum, fitting WMAP-9 and Planck data. We find that an inconsistency arises in predicting the scale factor at recombination, leading to a backreacted Universe with lower matter density, in contradiction with results produced by SnIa and clusters. The same behavior is confirmed by analyzing the CMB-shift parameters from WMAP-9. We conclude exploring qualitatively the motivations of this inconsistency.

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.

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 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.

We investigate the impact of statistical and systematic errors on measurements of linear redshift-space distortions (RSD) in future cosmological surveys, analyzing large catalogues of dark-matter halos from the BASICC simulation. These allow us to estimate the dependence of errors on typical survey properties, as volume, galaxy density and mass (i.e. bias factor) of the adopted tracer. We find that measures of the specific growth rate \beta=f/b using the Hamilton/Kaiser harmonic expansion of the redshift-space correlation function \xi(r_p,\pi) on scales larger than 3/h Mpc are typically under-estimated by up to 10% for galaxy sized halos. This is significantly larger than the corresponding statistical errors, which amount to a few percent, indicating the importance of non-linear improvements to the Kaiser model to obtain accurate measurements of the growth rate. We compare the statistical errors to predictions obtained with the Fisher information matrix, based on the usual FKP prescription for the errors on the power spectrum. We show that this produces parameter errors fairly similar to the standard deviations from the halo catalogues, but only if applied to strictly linear scales in Fourier space (k<0.2 h/Mpc). Finally, we present an accurate scaling formula describing the relative error on {\beta} as a function of the survey parameters, which closely matches the simulation results in all explored regimes. This provides a handy and plausibly more realistic alternative to the Fisher matrix approach, to quickly and accurately predict RSD statistical errors expected from future surveys.

The clustering of galaxies observed in future redshift surveys will provide a wealth of cosmological information. Matching the signal at different redshifts constrains the dark energy driving the acceleration of the expansion of the Universe. In tandem with these geometrical constraints, redshift-space distortions (RSD) depend on the build up of large-scale structure. As pointed out by many authors measurements of these effects are intrinsically coupled. We investigate this link, and argue that it strongly depends on the cosmological assumptions adopted when analysing data. Using representative assumptions for the parameters of the Euclid survey in order to provide a baseline future experiment, we show how the derived constraints change due to different model assumptions. We argue that even the assumption of a Friedman-Robertson-Walker (FRW) space-time is sufficient to reduce the importance of the coupling to a significant degree. Taking this idea further, we consider how the data would actually be analysed and argue that we should not expect to be able to simultaneously constrain multiple deviations from the standard $\Lambda$CDM model. We therefore consider different possible ways in which the Universe could deviate from the $\Lambda$CDM model, and show how the coupling between geometrical constraints and structure growth affects the measurement of such deviations.

We discuss the ability of the planned Euclid mission to detect deviations from General Relativity using its extensive redshift survey of more than 50 Million galaxies. Constraints on the gravity theory are placed measuring the growth rate of structure within 14 redshift bins between z=0.7 and z=2. The growth rate is measured from redshift-space distortions, i.e. the anisotropy of the clustering pattern induced by coherent peculiar motions. This is performed in the overall context of the Euclid spectroscopic survey, which will simultaneously measure the expansion history of the universe, using the power spectrum and its baryonic features as a standard ruler, accounting for the relative degeneracies of expansion and growth parameters. The resulting expected errors on the growth rate in the different redshift bins, expressed through the quantity f\sigma_8, range between 1.3% and 4.4%. We discuss the optimisation of the survey configuration and investigate the important dependence on the growth parameterisation and the assumed cosmological model. We show how a specific parameterisation could actually drive the design towards artificially restricted regions of the parameter space. Finally, in the framework of the popular "\gamma -parameterisation", we show that the Euclid spectroscopic survey alone will already be able to provide substantial evidence (in Bayesian terms) if the growth index differs from the GR value \gamma=0.55 by at least \sim 0.13. This will combine with the comparable inference power provided by the Euclid weak lensing survey, resulting in Euclid's unique ability to provide a decisive test of modified 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.

Euclid is a space-based survey mission from the European Space Agency designed to understand the origin of the Universe's accelerating expansion. It will use cosmological probes to investigate the nature of dark energy, dark matter and gravity by tracking their observational signatures on the geometry of the universe and on the cosmic history of structure formation. The mission is optimised for two independent primary cosmological probes: Weak gravitational Lensing (WL) and Baryonic Acoustic Oscillations (BAO). The Euclid payload consists of a 1.2 m Korsch telescope designed to provide a large field of view. It carries two instruments with a common field-of-view of ~0.54 deg2: the visual imager (VIS) and the near infrared instrument (NISP) which contains a slitless spectrometer and a three bands photometer. The Euclid wide survey will cover 15,000 deg2 of the extragalactic sky and is complemented by two 20 deg2 deep fields. For WL, Euclid measures the shapes of 30-40 resolved galaxies per arcmin2 in one broad visible R+I+Z band (550-920 nm). The photometric redshifts for these galaxies reach a precision of dz/(1+z) < 0.05. They are derived from three additional Euclid NIR bands (Y, J, H in the range 0.92-2.0 micron), complemented by ground based photometry in visible bands derived from public data or through engaged collaborations. The BAO are determined from a spectroscopic survey with a redshift accuracy dz/(1+z) =0.001. The slitless spectrometer, with spectral resolution ~250, predominantly detects Ha emission line galaxies. Euclid is a Medium Class mission of the ESA Cosmic Vision 2015-2025 programme, with a foreseen launch date in 2019. This report (also known as the Euclid Red Book) describes the outcome of the Phase A study.

A space-based galaxy redshift survey would have enormous power in constraining dark energy and testing general relativity, provided that its parameters are suitably optimized. We study viable space-based galaxy redshift surveys, exploring the dependence of the Dark Energy Task Force (DETF) figure-of-merit (FoM) on redshift accuracy, redshift range, survey area, target selection, and forecast method. Fitting formulae are provided for convenience. We also consider the dependence on the information used: the full galaxy power spectrum P(k, P(k) marginalized over its shape, or just the Baryon Acoustic Oscillations (BAO). We find that the inclusion of growth rate information (extracted using redshift space distortion and galaxy clustering amplitude measurements) leads to a factor of ~ 3 improvement in the FoM, assuming general relativity is not modified. This inclusion partially compensates for the loss of information when only the BAO are used to give geometrical constraints, rather than using the full P(k) as a standard ruler. We find that a space-based galaxy redshift survey covering ~20,000 deg^2 over 0.5< z < 2 with \sigma_z/(1+z)<= 0.001 exploits a redshift range that is only easily accessible from space, extends to sufficiently low redshifts to allow both a vast 3-D map of the universe using a single tracer population, and overlaps with ground-based surveys to enable robust modeling of systematic effects. We argue that these parameters are close to their optimal values given current instrumental and practical constraints.

We present a new systematic analysis of the early radiation era solution in an interacting dark energy model to find the adiabatic initial conditions for the Boltzmann integration. In a model where the interaction is proportional to the dark matter density, adiabatic initial conditions and viable cosmologies are possible if the early-time dark energy equation of state parameter is $w_e > -4/5$. We find that when adiabaticity between cold dark matter, baryons, neutrinos and photons is demanded, the dark energy component satisfies automatically the adiabaticity condition. As supernovae Ia or baryon acoustic oscillation data require the recent-time equation of state parameter to be more negative, we consider a time-varying equation of state in our model. In a companion paper [arXiv:0907.4987] we apply the initial conditions derived here, and perform a full Monte Carlo Markov Chain likelihood analysis of this model.

We use observations of cosmic microwave background anisotropies, supernova luminosities and the baryon acoustic oscillation signal in the galaxy distribution to constrain the cosmological parameters in a simple interacting dark energy model with a time-varying equation of state. Using a Monte Carlo Markov Chain technique we determine the posterior likelihoods. Constraints from the individual data sets are weak, but the combination of the three data sets confines the interaction constant $\Gamma$ to be less than 23% of the expansion rate of the Universe $H_0$; at 95% CL $-0.23 < \Gamma/H_0 < +0.15$. The CMB acoustic peaks can be well fitted even if the interaction rate is much larger, but this requires a larger or smaller (depending on the sign of interaction) matter density today than in the non-interacting model. Due to this degeneracy between the matter density and the interaction rate, the only observable effect on the CMB is a larger or smaller integrated Sachs-Wolfe (ISW) effect. While SN or BAO data alone do not set any direct constraints on the interaction, they exclude the models with very large matter density, and hence indirectly constrain the interaction rate when jointly analysed with the CMB data. To enable the analysis described in this paper, we present in a companion paper [arXiv:0907.4981] a new systematic analysis of the early radiation era solution to find the adiabatic initial conditions for the Boltzmann integration.

If dark energy interacts with dark matter, this gives a new approach to the coincidence problem. But interacting dark energy models can suffer from pathologies. We consider the case where the dark energy is modelled as a fluid with constant equation of state parameter w. Non-interacting constant-w models are well behaved in the background and in the perturbed universe. But the combination of constant w and a simple interaction with dark matter leads to an instability in the dark sector perturbations at early times: the curvature perturbation blows up on super-Hubble scales. Our results underline how important it is to carefully analyze the relativistic perturbations when considering models of coupled dark energy. The instability that we find has been missed in some previous work where the perturbations were not consistently treated. The unstable mode dominates even if adiabatic initial conditions are used. The instability also arises regardless of how weak the coupling is. This non-adiabatic instability is different from previously discovered adiabatic instabilities on small scales in the strong-coupling regime.

Recently H(z) data obtained from differential ages of galaxies have been proposed as a new geometrical probe of dark energy. In this paper we use those data, combined with other background tests (CMB shift and SNIa data), to constrain a set of general relativistic dark energy models together with some other models motivated by extra dimensions. Our analysis rests mostly on Bayesian statistics, and we conclude that LCDM is at least substantially favoured, and that braneworld models are less favoured than general relativistic ones.

Most parameterizations of the dark energy equation of state do not reflect realistic underlying physical models. Here, we develop a relatively simple description of dark energy based on the dynamics of a scalar field which is exact in the limit that the equation of state approaches a cosmological constant, assuming some degree of smoothness of the potential. By introducing just two parameters defined in the configuration space of the field we are able to reproduce a wide class of quintessence models. We examine the observational constraints on these models as compared to linear evolution models, and show how priors in the field space translate into priors on observational parameters.

We investigate a simple braneworld model in which the universe contains only cold dark matter and a cosmological constant, but the effective dark energy is phantom-like because of extra-dimensional gravity effects. Modified gravity screens the cosmological constant Lambda, allowing for a larger Lambda. In practice, observations do not favour any significant screening. We use supernova data, the cosmic microwave background shift parameter, and the baryon oscillation peak in the galaxy distribution to constrain the model. We find the mean value of Omega_m with 68% confidence limits, and an upper limit on Omega_Lambda at the 68% confidence level. The best-fit model is very close to a standard LCDM model, but the LCDM model provides a better fit since it has one less parameter.

The DGP brane-world model provides a simple alternative to the standard LCDM cosmology, with the same number of parameters. There is no dark energy - the late universe self-accelerates due to an infrared modification of gravity. We compute the joint constraints on the DGP model from supernovae, the cosmic microwave background shift parameter, and the baryon oscillation peak in the SDSS luminous red galaxy sample. Flat DGP models are within the joint 2 sigma contour, but the LCDM model provides a significantly better fit to the data. These tests are based on the background dynamics of the DGP model, and we comment on further tests that involve structure formation.

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.