Neutrinos and dark energy constraints from future galaxy surveys and CMB lensing information
Abstract
We explore the possibility of obtaining better constraints from future astronomical data by means of the Fisher information matrix formalism. In particular, we consider how cosmic microwave background (CMB) lensing information can improve our parameter error estimation. We consider a massive neutrino scenario and a timeevolving dark energy equation of state in the CDM framework. We use Planck satellite experimental specifications together with the future galaxy survey Euclid in our forecast. We found improvements in almost all studied parameters considering Planck alone when CMB lensing information is used. In this case, the improvement with respect to the constraints found without using CMB lensing is of 93% around the fiducial value for the neutrino parameter. The improvement on one of the dark energy parameter reaches 4.4%. When Euclid information is included in the analysis, the improvements on the neutrino parameter constraint is of approximately 128% around its fiducial value. The addition of Euclid information provides smaller errors on the dark energy parameters as well. For Euclid alone, the FoM is a factor of 29 higher than that from Planck alone even considering CMB lensing. Finally, the consideration of a nearly perfect CMB experiment showed that CMB lensing cannot be neglected specially in more precise future CMB experiments since it provided in our case a 6 times better FoM in respect to the unlensed CMB analysis .
pacs:
I Introduction
The discovery of the accelerated expansion of the universe (Riess et al., 1998; Perlmutter et al., 1999) can be interpreted by introducing in the cosmological model a negative pressure component, termed ”dark energy”. The simplest dark energy candidate is a cosmological constant , having constant equation of state . Together with a pressureless cold dark matter component (CDM), this constitutes the standard CDM model. Although this ”concordance model” is in very good agreement with a variety of cosmological observations (Larson et al., 2011; Komatsu et al., 2011), different candidates of dark energy cannot be discarded yet. Moreover, the basic cosmological constant scenario has two difficulties known as ”finetuning” and cosmic coincidence problems (see, e.g., Zlatev et al. (1997)). To overcome these problems alternative candidates for the dark energy have been proposed, such as the quintessence (Caldwell et al., 1998) that allows the possibility of a timedependent equation of state (Linder, 2003). In this paper, we will assume a redshift dependent equation of state for the dark energy,
(1) 
and adopt the wellknown ChevalierPlarskyLinder parametrization (Chevallier & Polarski, 2001; Linder, 2003)
(2) 
Cosmological observation can in principle be used to constrain the neutrinos’ masses. It was shown by neutrino oscillation experiments that neutrinos have nonzero masses (see (GonzalezGarcia et al., 2012) and references therein). However, these experiments can only constrain the neutrinos masssquare differences and not their individual values (for a review in neutrino masses see de Gouvea (2009)). On the other hand, cosmological probes are most sensitive to the total neutrino masses, . Using CMB radiation data only, from Planck satellite, an upper limit to of 0.933 eV at 95% C.L was found (Planck collaboration, 2013).
(3) 
However, the dark energy equation of state and the neutrinos’ total mass parameters are degenerated (see, e.g.,(Hannestad, 2005)). Some work has already been done to constrain both parameters simultaneously in a few dark energy scenarios, such as for models with a constant and a timevarying equations of state (Ichikawa & Takahashi, 2008; Hamann et al., 2010; Carbone et al., 2011; Joudaki & Kaplinghat, 2012).
Our goal is to forecast the constraint in total mass of neutrinos in a time evolving dark energy model, using CMB temperature and polarization power spectrum from the Planck satellite experimental setup (including also CMB lensing information), as well as the largescale matter distribution that can be observed by Euclid survey. We emphasize the usage of Planck CMB polarization information since its temperature data has been recently released (Planck collaboration, 2013). We assume a geometrically flat CDM model with two massive neutrinos with identical mass and one massless in a inverted hierarchy mass splitting, being = 0.125 eV for each massive neutrino. The fiducial parameters are , , , , . We normalize the CMB power spectra to COBE. For similar approaches see Hollenstein et al. (2009); Namikawa et al. (2010); Das et al. (2012); Hall & Challinor (2012); Hamann et al. (2012); Joudaki & Kaplinghat (2012) . The paper is organized as follows: We give a small introduction on CMB lensing in Sec. II. In Sec. III, we briefly review the Fisher information matrix formalism for the CMB (with and without lensing information) and for a galaxy survey. Finally present our results in Sec. IV, followed by our discussion and conclusions in Sec. V.
Ii CMB lensing
A small effect that can be observed in CMB power spectrum regards the deflection of photons, during their travel between the last scattering surface and the observer, by gravitational potentials due to clusters of galaxies. Smith et al. (2007) detected the CMB lensing signal for the first time by cross correlating WMAP data to radio galaxy counts in the NRAO VLA sky survey (NVSS). Recently, the detection of the gravitational lensing using CMB temperature maps alone and the measurement of the power spectrum of the projected gravitational potential were done using the Atacama Cosmology Telescope and the South Pole Telescope (Das et al., 2011; van Engelenl et al., 2012).
The lensing potential is defined as:
(4) 
being the comoving distance and is the conformal time at which the photon was at position .
The lensing effect remaps the temperature and polarization fields as
(5) 
(6) 
where in the case of the temperature field, the temperature T of the lensed CMB in a direction n̂ is equal to the unlensed CMB in a different direction n̂’. Both these directions, n̂ and n̂’, differ by the deflection angle . To first order, the deflection angle is simply the lensing potential gradient, . In the same way, the effect of lensing in CMB polarization is written in terms of the Stokes parameters and (for a review in CMB polarization theory, see Cabella & Kamionkowski (2004)) :
CMB lensing has important quantitative contributions that should be taken into account, therefore a lot of work has been done to CMB lensing reconstruction techniques (e.g. Hu (2001); Okamoto & Hu (2003); Smith et al. (2012); Bucher et al. (2010); Carvalho & Tereno (2011)). In this paper, we use the CAMB software package (Lewis et al., 2000) to obtain the numerical lensed and unlensed power spectra ( and ) for our cosmological model with . We then use this predictions to forecast how CMB lensing information will help us constraining our model.
Iii Method
We apply the Fisher information matrix formalism to a Plancklike experiment (Blue book, 2005), considering both temperature and polarization for the lensend and unlensed CMB spectrum, and to an experiment such as the future Euclid survey. We forecast the dark energy and the massive neutrinos parameters in our fiducial model. In addition, we check the impact of CMB lensing information on the constraints of the mentioned parameters in a nearly perfect CMB experiment.
iii.1 Information from CMB
The Fisher information matrix for the CMB temperature anisotropy and polarization is given by (Zaldarriaga & Seljak, 1997)
(7) 
where is the power in the th multipole, stands for (temperature), (Emode polarization), (Bmode polarization) and (temperature and Emode polarization crosscorrelation). We will not include primordial Bmodes in the analysis since the measurement of the primordial by Planck is expected to be noise dominated. Our covariance matrix becomes therefore:
(8) 
Explicit expressions for the matrix elements are given in the appendix Section A.
For the lensed case we performed a correction in the covariance matrix elements taking into consideration the power spectrum of the deflection angle and its cross correlation with temperature and Epolarization, and . We used the same procedure introduced in (Perotto et al., 2006) to obtain the covariance matrix elements using the new information of power spectrum (see appendix Section B).
We also change in this case the unlensed CMB power spectra, , for the lensed ones, . Note that in this case we are taking into consideration the Bmode polarization generated by the CMB gravitational lensing from the Emode polarization.
When we include these corrections, the covariance matrix becomes:
(9) 
(10) 
(11) 
(12) 
(13) 
(14) 
(15) 
(16) 
(17) 
(18) 
(19) 
(20) 
(21) 
(22) 
(23) 
(24) 
(25) 
(26) 
(27) 
(28) 
(29) 
(30) 
(31) 
In these equations, and are the gaussian random detector noises for temperature and polarization respectively, which expression is written using the window function, and the inverse square of the detector noise level for temperature and polarization, and . The Full Width Half Maximum (FWHM), , is used in radians and is the weight given to each considered Planck channel (Eisenstein et al., 1999) . We tested two types of experimental setups that can be checked in Tables 2 and 2 .
(32) 
(33) 
Here we used four channels, 100, 143, 217 and 353GHz of the Planck experiment as can be seen from the equations 32 and 33.
Frequency (GHz)  

100  9.5’  6.82  10.9120 
143  7.1’  6.0016  11.4576 
217  5.0’  13.0944  26.7644 
353  5.0’  40.1016  81.2944 
4.0’ 

In addition, is the optimal quadratic estimator noise of the deflection field (we consider only the TT quadratic estimator noise for the planck experiment since it provides the best estimator). For the nearly ideal experiment we consider the minimum variance (MV) estimator noise written as a combination of the noises TB, TT, TE, EE and EB of the quadratic estimators (for a review in the topic, see Okamoto & Hu (2003); Hu & Okamoto (2002))). Figure 1 shows the quadratic estimator noises for our fiducial model considering Planck and the nearly ideal experiment.
iii.2 Information from galaxy survey: Baryonic acoustic oscillation
We show here how the baryonic acoustic oscillation (BAO) information can be used to forecast errors in the dark energy parameters using the fisher formalism. It was shown by Seo & Eisenstein (2003) that the Hubble parameter H(z) and the angular diameter distance Da(z) can be measured very precisely by using the BAO information present in the matter power spectrum. H(z) and Da(z) are expected to be determined as a function of redshift by future galaxy surveys. The goal is then to propagate the errors on H(z) and Da(z) to the constraints of dark energy parameters.
We start defining the observed galaxy power spectrum in a reference cosmology (in our case we use the CDM model), distinguished by the subscript ”ref”, (different from the true spectrum, referred as no subscript) that will be used to derive the cosmological parameters constraints using a galaxy survey that covers a wide range of redshifts. Following Seo & Eisenstein (2003),
(34) 
Where the Hubble parameter H(z) in a flat Universe is related to the dark energy parameters through
(35) 
and the angular diameter distance is defined as
(36) 
.
is the unknown Poisson shot noise.
The wavenumbers across and along the line of sight are denoted by and . It is important to point out that the wavenumbers in the reference cosmology are related to the ones in the true cosmology by
(37) 
Moreover, we define the galaxy power spectrum, , including the redshift distortions:
(38) 
where k r̂/k , being r̂ the unit vector along the line of sight and the linear matter power spectrum, , was generated using CAMB software package (Lewis et al., 2000) and COBE normalized. The is chosen in a way to exclude information from non linear regime where equation 34 is inaccurate (see Seo & Eisenstein (2003)). For an approach considering the nonlinear regime see, for example, Seo & Eisenstein (2005, 2007); Wang (2010, 2012). Moreover, for the impact of precisely modeling systematic effects, such as the nonlinear clustering and redshift space distortions, in the evolution of BAO see Crocce & Scoccimarro (2008); Taruya et al. (2009, 2010). The exponential damping factor is due to redshift uncertainties, where . , and are the growth function, the linear redshift space distortion parameter and the linear galaxy bias respectively. We use a growth factor dependent on the dark energy parameter and massive neutrinos effect computed by Kiakotou et al. (2008). The growth rate of perturbations is defined as
(39) 
where the growth function G(z) is related to the density of matter. In a matter dominated Universe , with . More generally, we use with
(40) 
is a numerical function dependent on and (see Equation 17 and Table 5 of Kiakotou et al. (2008)).
The linear redshift space distortion is also defined as a function os the growth rate and the galaxy bias.
(41) 
iii.2.1 Fisher formalism
The Fisher information matrix for the matter power spectrum obtained from galaxy surveys is given by (Tegmark, 1997)
(42) 
The effective volume of the survey for a constant comoving number density is given by
(43) 
We use information of an Euclid like survey with area of 20000 , redshift accuracy of and a redshift range . Finally we dived our forecast into 15 redshift slices of centered in . We chose the initial set of parameters . For each redshift bin we use the specifications on Table 3 (see Pavlov et al. (2012); di Porto et al. (2012) and references therein).
()  b(z)  n(z)  

0.55  0.144  1.0423  3.56 
0.65  0.153  1.0668  3.56 
0.75  0.163  1.1084  2.42 
0.85  0.174  1.1145  2.42 
0.95  0.185  1.1107  1.81 
1.05  0.197  1.1652  1.81 
1.15  0.2  1.2262  1.44 
1.25  0.2  1.2769  1.44 
1.35  0.2  1.2960  0.99 
1.45  0.2  1.3159  0.99 
1.55  0.2  1.4416  0.55 
1.65  0.2  1.4915  0.55 
1.75  0.2  1.4973  0.29 
1.85  0.2  1.5332  0.29 
1.95  0.2  1.5705  0.15 
To obtain the constraints on our final set of parameters , first we marginalize our first fisher matrix over and use this sub matrix to change into the desired variables as
(44) 
iii.3 Information from galaxy survey: Weak lensing
In this subsection, we show how weak lensing (WL) information can improve the constraints on cosmological parameters, in our case specially for dark energy and neutrinos’ densities parameters using the fisher formalism. The observable, in weak lensing surveys, is the convergence power spectrum. In the analysis presented in this paper, we use an extension of the CAMB software with Halofit approximation (recently updated according to Takahashi et al. (2012)) to generate the convergence power spectra , where the subscripts and stand to the lensed galaxies redshift bins. The fisher matrix for weak lensing is then given by (Hu & Tegmark, 1999)
(45) 
The covariance matrix is defined as
(46) 
being the r.m.s intrinsic shear and the number of galaxies per steradian in the jth bin
(47) 
In the equation above, d is the number of galaxies per square arcminute and is fraction of sources belonging to a certain bin. We compute our calculations considering an Euclidlike experiment following Blue book (2008) with , and . We take the range and divide it into 4 equalgalaxynumber bins. We also consider . In this case, we tested the analysis for a maximum multipole of 3000 and no significant change was found confirming that both larger and smaller multipoles do not give a significant contribution to the results (see Hollenstein et al. (2009)). The photoz error is assumed to be normal distributed with variance . It is important to point out that nongaussian errors can be significant in the measurements of weak lensing, degrading the signaltonoise ratio of convergence power spectrum (Sato et al., 2009) and therefore the marginalized errors on individual parameters by a few percent (Takada & Jain, 2009).
Iv Results
We performed the forecast for Planck alone, with and without considering CMB lensing. Moreover, we introduced the Euclid forecast, combining the results approximately as
(48) 
It was shown by Hollenstein et al. (2009) that the covariance between the measurements of cosmic shear tomography and the CMB lensing can be safely neglected since the redshifts in which they are probed are quite distinct from each other.
We found the best limits for the neutrino density, and in the combined Planck (with lensing)+ Euclid (BAO +WL). We have that , and (95% C.L) as it can be inferred from column 6 of Table 4 . The Figure of Merit (FoM), described as the reciprocal of the 95% confidence limit’s area of the error ellipse from the plane  (Albrecht et al., 2006), of Euclid (BAO +WL) is a factor 29 higher than the FoM for Planck alone even when CMB lensing is considered. Euclid will be able to strongly constraint the lateuniverse parameters. The combined result Planck (with lensing)+ Euclid (BAO +WL) still improves the FoM in respect to Euclid (BAO +WL) in a factor of 3. In Figure 2 we show the 2 sigma fisher contours.
In table 5 we show how CMB lensing could affect the parameters constraints for an almost ideal experiment. In this case, we see a substantial improvement of about 6 times in the FoM when CMB lensing is considered. Comparing also with Planck experiment, the nearly ideal experiment improves the FoM by a factor of 5 without considering lensing information in any case. For the neutrino parameter we have that without lensing and when lensing CMB is considered. On the other hand, the use of the cross power spectrum between the deflection angle and the E mode polarization makes no significant impact on any cosmological parameter constraint. Figure 3 shows the fisher contours for the unlensed and lensed analysis. Note that the 2 sigma contour obtained when we use and power spectra overlap the 2 sigma contour when we also add the power spectrum.
Parameter  Fiducial  Planck  Planck  EUCLID  Planck + EUCLID 

value  T + P  T+ P+ lens  (BAO + WL)  

0.02219  0.00012  0.00012  0.00034  7.9e05 
0.01122  0.00080  0.00070  0.00011  8.2e05  
0.0027  0.0036  0.0011  0.00035  0.00013  
0.95  0.084  0.042  0.0027  0.0015  
0  0.084  0.057  0.036  0.015  
FoM    8.21  25.81  732.97  2909.13 
Relative FoM ^{3}^{3}3Relative FoM in respecto to Planck (T+P) without CMB lensing.    1  3.14  89.3  354.34 
Parameter  Fiducial  T + P+  T+ P  

value  T+P (unlensed)  lens ( and  lens (, and )  

0.02219  2.4778e05  2.2296e05  2.2285e05 
0.01122  0.0004420  0.0003073  0.0003071  
0.0027  0.0015860  0.0002065  0.0002064  
0.95  0.0338432  0.0142514  0.0142416  
0  0.0338432  0.0238881  0.0238811  
FoM    41.83  255.43  255.63 
Relative FoM    1  6.106  6.111 
V Discussion and conclusions
We are considering massive neutrinos and a time evolving equation of state in the CDM model. Using the Fisher formalism, we obtained the best constraints possible for , and considering Planck and Euclid survey.
One of the goals of this work has been to quantify the influence of CMB lensing information in the constraints of the parameters of interest, especially , and . We saw on columns 3 and 4 from Table 4 that we improve the constraints in basically all the studied cosmological parameters. The improvement found on for the Planck 1 sigma error alone varied from approximately 133% to 40% from its fiducial value without using CMB lensing and using CMB lensing information respectively. For the error is 4.4% smaller when CMB lensing is taken into consideration. The 2 sigma constraint on from Planck alone varies from (without CMB lensing) and .
When we add Euclid information to Planck information, we get an impressive improvement of 128.2% on the 1 sigma error of considering the results from Planck without lensing and Planck (+ lensing) + Euclid (BAO + WL). An improvement of approximately 9% in the error of was found when including Euclid information to the Planck forecast.
In the case of a nearly perfect CMB experiment, as mentioned before, CMB lensing improved all the constraints of the tested parameters. It is clear from the analysis that CMB lensing can play an important role in constraining cosmological parameters in future CMB experiments and must be taken into account. On another hand, the power spectra can be safely neglected in near future CMB experiments since its contribution to the parameters constraints is minimum.
Acknowledgements.
The authors would like to thank Carmelita Carbone and Yun Wang for useful discussions. This work has been partially supported by ASI contract for Planck LFI activity of Phase E2.Appendix A Elements of the unlensed covariance matrix
The elements of the covariance matrix in the unlensed case are:
(49) 
(50) 
(51) 
(52) 
(53) 
(54) 
(55) 
(56) 
Appendix B Elements of the lensed covariance matrix
First of all, we make use of the effective defined in Equation (3.3) of Perotto et al. (2006).
(57) 
where is our case D is defined to be
(58) 
and , are the determinants of the theoretical and observed data covariance matrices
(59) 
(60) 
The theoretical covariance matrix M is given by
(61) 
The fisher matrix information is then derived from the second order derivative of the likelihood function, , from an observing data set given the real parameters :
(62) 
knowing that , we derived the new elements for the covariance matrix 9.
References
 Riess et al. (1998) Riess, A. G., et al., AJ, 116,10091038(1998)
 Perlmutter et al. (1999) Perlmutter, L., et al., AJ, 517, 565586(1999)
 Komatsu et al. (2011) E. Komatsu et al., Astrophys. J. Supplement, 192,18 (2011).
 Larson et al. (2011) D. Larson et al. , Astrophys. J. Supplement, 192, 16+ (2011)
 Zlatev et al. (1997) I. Zlatev, L. Wang, and P. J. Steinhardt, Phys. Rev. Lett. , 82, 1896899(1999)
 Caldwell et al. (1998) R. R. Caldwell, R. Dave and P. J. Steinhardt, Phys. Rev. Lett. , 80, 15821585(1998)
 Linder (2003) E. V. Linder, Phys. Rev. Lett. , 90, 091301(2003)
 Chevallier & Polarski (2001) M. Chevallier and D. Polarski, Int. J.. Mod. Phys. D, 10, 213(2001)
 GonzalezGarcia et al. (2012) M. C. GonzalesGarcia, M. Maltoni, J. Salvado and T. Schwetz, J. ArXiv Astrophysics eprints, hepph/1209.3023 (2012).
 de Gouvea (2009) A. de Gouvea, ArXiv Astrophysics eprints, hepph/0902.4656 (2009)
 Hannestad (2005) S. Hannestad, Phys. Rev. Lett. , 95, 221301(2005)
 Carbone et al. (2011) C. Carbone, L. Verde,Y. Wang and A. Cimatti, J. Cosmol. Astropart. Phys., 3, 30 (2011).
 Hamann et al. (2010) J. Hamann et al., J. Cosmol. Astropart. Phys., 7, 22 (2010)
 Ichikawa & Takahashi (2008) K. Ichikawa and T. Takahashi, J. Cosmol. Astropart. Phys., 2, 17 (2008)
 Joudaki & Kaplinghat (2012) S. Joudaki and M. Kaplinghat, Phys. Rev. D, 86,, 023526 (2012)
 Planck collaboration (2013) Planck collaboration, ArXiv Astrophysics eprints, Astroph/1303.5076(2013)
 Hollenstein et al. (2009) L. Hollenstein, D. Sapone,R. Crittenden and B. M. Schä fer, J. Cosmol. Astropart. Phys., 4, 12 (2009)
 Albrecht et al. (2006) A. Albrecht et al., ArXiv Astrophysics eprints, Astroph/0609591(2006)
 Namikawa et al. (2010) T. Namikawa, S. Saito and A. Taruya Cosmol. Astropart. Phys., 12, 27 (2010) .
 Das et al. (2012) S. Das, R. de Putter, E. V. Linder and R. Nakajima, J. Cosmol. Astropart. Phys., 11, 11 (2012).
 Hall & Challinor (2012) A. C. Hall and A. Challinor, MNRAS , 425, 11701184 (2012)
 Hamann et al. (2012) J. Hamann, S. Hannestad and Y. Y. Y. Wong, J. Cosmol. Astropart. Phys. , 11, 52 (2012)
 Smith et al. (2007) K. M. Smith, O. Zahn, and O. Doré, Phys. Rev. D, 76, 043510 (2007).
 Das et al. (2011) S. Das et al., Phys. Rev. Lett., 107, 021301 (2011).
 van Engelenl et al. (2012) A. van Engelen et al., Astrophys. J. , 756, 142 (2012)
 Cabella & Kamionkowski (2004) Cabella, P. & Kamionkowski, M., at the 2003 Villa Mondragone School of Gravitation and Cosmology: ”The Polarization of the Cosmic Microwave Background” (Rome, Italy, 2003); arXiv: astroph/0403392.
 Lewis & Challinor (2006) A. Lewis and A. Challinor, Phys. Rep., 429, 165 (2006).
 Bucher et al. (2010) M. Bucher, B. van Tent, and C.S. Carvalho, MNRAS, 407, 21932206 (2010).
 Carvalho & Tereno (2011) C. S. Carvalho and I. Tenero, Phys. Rev. D, 84, 063001 (2011).
 Hu (2001) W. Hu, Astrophys. J. , 557, L79 (2001).
 Okamoto & Hu (2003) T. Okamoto and W. Hu, Phys. Rev. D, 67, 083002 (2003).
 Smith et al. (2012) K. M. Smith, D. Hanson, M. LoVerde, C. M. Hirata, and O. Zahn, J. Cosmol.Astropart. Phys., 6, 14 (2012)
 Lewis et al. (2000) A. Lewis, A. Challinor and A. Lasenby, Astrophys. J. , 538, 473476 (2000).
 Blue book (2005) http://www.rssd.esa.int/SA/PLANCK/docs/BluebookESASCI(2005)1_V2.pdf
 Eisenstein et al. (1999) D. Eisenstein, W. Hu, and M. Tegmark, Astrophys. J. , 518,223 (1999).
 Perotto et al. (2006) L. Perotto, J. Lesgourgues, S. Hannestad, H. Tu, and Y. Y. Y. Wong, J. Cosmol. Astropart. Phys., 10,13 (2006).
 Zaldarriaga & Seljak (1997) M. Zaldarriaga and U. Seljak, Phys. Rev. D, 55, 18301840 (1997).
 Hu & Okamoto (2002) W. Hu and T. Okamoto, Astrophys. J. , 574, 566 (2002).
 Hu (2002) W. Hu Phys. Rev. D, 65, 023003 (2002)
 Seo & Eisenstein (2003) H.J. Seo and D. J. Eisenstein, Astrophys. J. , 598, 720740 (2003).
 Seo & Eisenstein (2005) H.J. Seo and D. J. Eisenstein, Astrophys. J. , 633, 575588 (2005).
 Seo & Eisenstein (2007) H.J. Seo and D. J. Eisenstein, Astrophys. J. , 667, 1424 (2007).
 Wang (2010) Y. Wang, Modern Physics Letters A, 25, 30933113 (2010).
 Wang (2012) Y. Wang, MNRAS, 423, 36313637 (2012).
 Crocce & Scoccimarro (2008) M. Crocce and R. Scoccimarro, Phys. Rev. D, 77, 023533 (2008).
 Taruya et al. (2009) A. Taruya, T. Nishimichi and S. Saito and T. Hiramatsu, Phys. Rev. D, 80, 123503 (2009).
 Taruya et al. (2010) A. Taruya, T. Nishimichi and S. Saito, Phys. Rev. D, 82, 063522 (2010).
 Kiakotou et al. (2008) A. Kiakotou, Ø. Elgarøy and O. Laha, Phys. Rev. D, 77,063005 (2008).
 Pavlov et al. (2012) A. Pavlov, L. Samushia and B. Ratra, ArXiv Astrophysics eprints, Astroph/12063123(2012)
 di Porto et al. (2012) C. di Porto, L. Amendola and E. Branchini, MNRAS, 419, 985997 (2012).
 Takahashi et al. (2012) R. Takahashi, M. Sato and T. Nishimichi and A. Taruya and M. Oguri, Astrophys. J. , 761, 152 (2012).
 Hu & Tegmark (1999) W. Hu and M. Tegmark, Astrophys. J. letters, 514, L65L68 (1999).
 Sato et al. (2009) M. Sato et al. , Astrophys. J. , 701, 945954 (2009).
 Takada & Jain (2009) M. Takada and B. Jain, MNRAS, 395, 20652086 (2009).
 Blue book (2008) http://sci2.esa.int/cosmicvision/cycle1/call_plstudyconsortia/Euclid_Payload_Definition_Document.pdf
 Albrecht et al. (2006) A. Albrecht et al., ArXiv Astrophysics eprints, Astroph/0609591(2006)
 Tegmark (1997) M. Tegmark, Phys. Rev. Lett. , 79, 38063809 (1997).