No Arabic abstract
An accurate theoretical template for the galaxy power spectrum is a key for the success of ongoing and future spectroscopic surveys. We examine to what extent the Effective Field Theory of Large Scale Structure is able to provide such a template and correctly estimate cosmological parameters. To that end, we initiate a blinded challenge to infer cosmological parameters from the redshift-space power spectrum of high-resolution mock catalogs mimicking the BOSS galaxy sample but covering a hundred times larger cumulative volume. This gigantic simulation volume allows us to separate systematic bias due to theoretical modeling from the statistical error due to sample variance. The challenge task was to measure three unknown input parameters used in the simulation: the Hubble constant, the matter density fraction, and the clustering amplitude. We present analyses done by two independent teams, who have fitted the mock simulation data generated by yet another independent group. This allows us to avoid any confirmation bias by analyzers and pin down possible tuning of the specific EFT implementations. Both independent teams have recovered the true values of the input parameters within sub-percent statistical errors corresponding to the total simulation volume.
Published galaxy power spectra from the 2dFGRS and SDSS are not in good agreement. We revisit this issue by analyzing both the 2dFGRS and SDSS DR5 catalogues using essentially identical techniques. We confirm that the 2dFGRS exhibits relatively more large scale power than the SDSS, or, equivalently, SDSS has more small scale power. We demonstrate that this difference is due to the r-band selected SDSS catalogue being dominated by more strongly clustered red galaxies, which have a stronger scale dependent bias. The power spectra of galaxies of the same rest frame colours from the two surveys match well. If not accounted for, the difference between the SDSS and 2dFGRS power spectra causes a bias in the obtained constraints on cosmological parameters which is larger than the uncertainty with which they are determined. We also found that the correction developed by Cole et al.(2005) to model the distortion in the shape of the power spectrum due to non-linear evolution and scale dependent bias is not able to reconcile the constraints obtained from the 2dFGRS and SDSS power spectra. Intriguingly, the model is able to describe the differences between the 2dFGRS and the much more strongly clustered LRG sample, which exhibits greater nonlinearities. This shows that more work is needed to understand the relation between the galaxy power spectrum and the linear perturbation theory prediction for the power spectrum of matter fluctuations. It is therefore important to accurately model these effects to get precise estimates of cosmological parameters from these power spectra and from future galaxy surveys like Pan-STARRS, or the Dark Energy Survey, which will use selection criteria similar to the one of SDSS.
Upcoming weak lensing surveys, such as LSST, EUCLID, and WFIRST, aim to measure the matter power spectrum with unprecedented accuracy. In order to fully exploit these observations, models are needed that, given a set of cosmological parameters, can predict the non-linear matter power spectrum at the level of 1% or better for scales corresponding to comoving wave numbers 0.1<k<10 h/Mpc. We have employed the large suite of simulations from the OWLS project to investigate the effects of various baryonic processes on the matter power spectrum. In addition, we have examined the distribution of power over different mass components, the back-reaction of the baryons on the CDM, and the evolution of the dominant effects on the matter power spectrum. We find that single baryonic processes are capable of changing the power spectrum by up to several tens of per cent. Our simulation that includes AGN feedback, which we consider to be our most realistic simulation as, unlike those used in previous studies, it has been shown to solve the overcooling problem and to reproduce optical and X-ray observations of groups of galaxies, predicts a decrease in power relative to a dark matter only simulation ranging, at z=0, from 1% at k~0.3 h/Mpc to 10% at k~1 h/Mpc and to 30% at k~10 h/Mpc. This contradicts the naive view that baryons raise the power through cooling, which is the dominant effect only for k>70 h/Mpc. Therefore, baryons, and particularly AGN feedback, cannot be ignored in theoretical power spectra for k>0.3 h/Mpc. It will thus be necessary to improve our understanding of feedback processes in galaxy formation, or at least to constrain them through auxiliary observations, before we can fulfil the goals of upcoming weak lensing surveys.
We compare analytical computations with numerical simulations for dark-matter clustering, in general relativity and in the normal branch of DGP gravity (nDGP). Our analytical frameword is the Effective Field Theory of Large-Scale Structure (EFTofLSS), which we use to compute the one-loop dark-matter power spectrum, including the resummation of infrared bulk displacement effects. We compare this to a set of 20 COLA simulations at redshifts $z = 0$, $z=0.5$, and $z =1$, and fit the free parameter of the EFTofLSS, called the speed of sound, in both $Lambda$CDM and nDGP at each redshift. At one-loop at $z = 0$, the reach of the EFTofLSS is $k_{rm reach}approx 0.14 , h { rm Mpc^{-1}}$ for both $Lambda$CDM and nDGP. Along the way, we compare two different infrared resummation schemes and two different treatments of the time dependence of the perturbative expansion, concluding that they agree to approximately $1%$ over the scales of interest. Finally, we use the ratio of the COLA power spectra to make a precision measurement of the difference between the speeds of sound in $Lambda$CDM and nDGP, and verify that this is proportional to the modification of the linear coupling constant of the Poisson equation.
In this paper we present results from the weak lensing shape measurement GRavitational lEnsing Accuracy Testing 2010 (GREAT10) Galaxy Challenge. This marks an order of magnitude step change in the level of scrutiny employed in weak lensing shape measurement analysis. We provide descriptions of each method tested and include 10 evaluation metrics over 24 simulation branches. GREAT10 was the first shape measurement challenge to include variable fields; both the shear field and the Point Spread Function (PSF) vary across the images in a realistic manner. The variable fields enable a variety of metrics that are inaccessible to constant shear simulations including a direct measure of the impact of shape measurement inaccuracies, and the impact of PSF size and ellipticity, on the shear power spectrum. To assess the impact of shape measurement bias for cosmic shear we present a general pseudo-Cl formalism, that propagates spatially varying systematics in cosmic shear through to power spectrum estimates. We also show how one-point estimators of bias can be extracted from variable shear simulations. The GREAT10 Galaxy Challenge received 95 submissions and saw a factor of 3 improvement in the accuracy achieved by shape measurement methods. The best methods achieve sub-percent average biases. We find a strong dependence in accuracy as a function of signal-to-noise, and indications of a weak dependence on galaxy type and size. Some requirements for the most ambitious cosmic shear experiments are met above a signal-to-noise ratio of 20. These results have the caveat that the simulated PSF was a ground-based PSF. Our results are a snapshot of the accuracy of current shape measurement methods and are a benchmark upon which improvement can continue. This provides a foundation for a better understanding of the strengths and limitations of shape measurement methods.
We study the accuracy with which cosmological parameters can be determined from real space power spectrum of matter density contrast at weakly nonlinear scales using analytical approaches. From power spectra measured in $N$-body simulations and using Markov chain Monte-Carlo technique, the best-fitting cosmological input parameters are determined with several analytical methods as a theoretical template, such as the standard perturbation theory, the regularized perturbation theory, and the effective field theory. We show that at redshift 1, all two-loop level calculations can fit the measured power spectrum down to scales $k sim 0.2 , h , mathrm{Mpc}^{-1}$ and cosmological parameters are successfully estimated in an unbiased way. Introducing the Figure of bias (FoB) and Figure of merit (FoM) parameter, we determine the validity range of those models and then evaluate their relative performances. With one free parameter, namely the damping scale, the regularized perturbation theory is found to be able to provide the largest FoM parameter while keeping the FoB in the acceptance range.