Modeling the angular correlation function and its full covariance in Photometric Galaxy Surveys

Near future cosmology will see the advent of wide area photometric galaxy surveys, like the Dark Energy Survey (DES), that extent to high redshifts (z ~ 1 - 2) but with poor radial distance resolution. In such cases splitting the data into redshift bins and using the angular correlation function $w(\theta)$, or the $C_{\ell}$ power spectrum, will become the standard approach to extract cosmological information or to study the nature of dark energy through the Baryon Acoustic Oscillations (BAO) probe. In this work we present a detailed model for $w(\theta)$ at large scales as a function of redshift and bin width, including all relevant effects, namely nonlinear gravitational clustering, bias, redshift space distortions and photo-z uncertainties. We also present a model for the full covariance matrix characterizing the angular correlation measurements, that takes into account the same effects as for $w(\theta)$ and also the possibility of a shot-noise component and partial sky coverage. Provided with a large volume N-body simulation from the MICE collaboration we built several ensembles of mock redshift bins with a sky coverage and depth typical of forthcoming photometric surveys. The model for the angular correlation and the one for the covariance matrix agree remarkably well with the mock measurements in all configurations. The prospects for a full shape analysis of $w(\theta)$ at BAO scales in forthcoming photometric surveys such as DES are thus very encouraging.Comment: 23 pages, 21 figures Revised version accepted by MNRAS. Description of mocks re-structured. Mocks including redshift distortions and Photo-z publicly available at http://www.ice.cat/mic

Measuring the growth of matter fluctuations with third-order galaxy correlations

Measurements of the linear growth factor $D$ at different redshifts $z$ are key to distinguish among cosmological models. One can estimate the derivative $dD(z)/d\ln(1+z)$ from redshift space measurements of the 3D anisotropic galaxy two-point correlation $\xi(z)$, but the degeneracy of its transverse (or projected) component with galaxy bias $b$, i.e. $\xi_{\perp}(z) \propto\ D^2(z) b^2(z)$, introduces large errors in the growth measurement. Here we present a comparison between two methods which break this degeneracy by combining second- and third-order statistics. One uses the shape of the reduced three-point correlation and the other a combination of third-order one- and two-point cumulants. These methods use the fact that, for Gaussian initial conditions and scales larger than $20$ $h^{-1}$Mpc, the reduced third-order matter correlations are independent of redshift (and therefore of the growth factor) while the third-order galaxy correlations depend on $b$. We use matter and halo catalogs from the MICE-GC simulation to test how well we can recover $b(z)$ and therefore $D(z)$ with these methods in 3D real space. We also present a new approach, which enables us to measure $D$ directly from the redshift evolution of second- and third-order galaxy correlations without the need of modelling matter correlations. For haloes with masses lower than $10^{14}$ $h^{-1}$M$_\odot$, we find $10$ deviations between the different estimates of $D$, which are comparable to current observational errors. At higher masses we find larger differences that can probably be attributed to the breakdown of the bias model and non-Poissonian shot noise.Comment: 24 pages, 20 figures, 2 tables, accepted for publication in MNRA

Large-N expansions applied to gravitational clustering

We develop a path-integral formalism to study the formation of large-scale structures in the universe. Starting from the equations of motion of hydrodynamics (single-stream approximation) we derive the action which describes the statistical properties of the density and velocity fields for Gaussian initial conditions. Then, we present large-N expansions (associated with a generalization to N fields or with a semi-classical expansion) of the path-integral defined by this action. This provides a systematic expansion for two-point functions such as the response function and the usual two-point correlation. We present the results of two such expansions (and related variants) at one-loop order for a SCDM and a LCDM cosmology. We find that the response function exhibits fast oscillations in the non-linear regime with an amplitude which either follows the linear prediction (for the direct steepest-descent scheme) or decays (for the 2PI effective action scheme). On the other hand, the correlation function agrees with the standard one-loop result in the quasi-linear regime and remains well-behaved in the highly non-linear regime. This suggests that these large-N expansions could provide a good framework to study the dynamics of gravitational clustering in the non-linear regime. Moreover, the use of various expansion schemes allows one to estimate their range of validity without the need of N-body simulations and could provide a better accuracy in the weakly non-linear regime.Comment: 27 pages, published in A&

Modeling Scale-Dependent Bias on the Baryonic Acoustic Scale with the Statistics of Peaks of Gaussian Random Fields

Models of galaxy and halo clustering commonly assume that the tracers can be treated as a continuous field locally biased with respect to the underlying mass distribution. In the peak model pioneered by Bardeen et al. [Astrophys. J. 304, 15 (1986)], one considers instead density maxima of the initial, Gaussian mass density field as an approximation to the formation site of virialized objects. In this paper, the peak model is extended in two ways to improve its predictive accuracy. First, we derive the two-point correlation function of initial density peaks up to second order and demonstrate that a peak-background split approach can be applied to obtain the k-independent and k-dependent peak bias factors at all orders. Second, we explore the gravitational evolution of the peak correlation function within the Zel’dovich approximation. We show that the local (Lagrangian) bias approach emerges as a special case of the peak model, in which all bias parameters are scale independent and there is no statistical velocity bias. We apply our formulas to study how the Lagrangian peak biasing, the diffusion due to large scale flows, and the mode coupling due to nonlocal interactions affect the scale dependence of bias from small separations up to the baryon acoustic oscillation (BAO) scale. For 2σ density peaks collapsing at z = 0.3, our model predicts a ~5% residual scale-dependent bias around the acoustic scale that arises mostly from first order Lagrangian peak biasing (as opposed to second order gravity mode coupling). We also search for a scale dependence of bias in the large scale autocorrelation of massive halos extracted from a very large N-body simulation provided by the MICE Collaboration. For halos with mass M ≳ 1014M⊙/h, our measurements demonstrate a scale-dependent bias across the BAO feature which is very well reproduced by a prediction based on the peak model

Testing one-loop galaxy bias: Cosmological constraints from the power spectrum

We investigate the impact of different assumptions in the modeling of one-loop galaxy bias on the recovery of cosmological parameters, as a follow-up of the analysis done in the first paper of the series at fixed cosmology. To carry out these tests we focus on the real-space galaxy-power spectrum from a set of three different synthetic galaxy samples whose clustering properties are meant to match the ones of the CMASS and LOWZ catalogs of BOSS and the SDSS Main Galaxy Sample. We investigate the relevance of allowing for either short range nonlocality or scale-dependent stochasticity by fitting the real-space galaxy autopower spectrum or the combination of galaxy-galaxy and galaxy-matter power spectrum. From a comparison among the goodness of fit (χ2), unbiasedness of cosmological parameters (FoB), and figure of merit (FoM) of the model, we find that a simple four-parameter model (linear, quadratic, cubic nonlocal bias, and constant shot noise) with fixed quadratic tidal bias provides a robust modeling choice for the autopower spectrum of the three galaxy samples, up to kmax ¼ 0.3h Mpc−1 and for an effective volume of 6h−3 Gpc3. Instead, a joint analysis of the two observables fails at larger scales, and a model extension with either higher derivatives or scale-dependent shot noise is necessary to reach a similar kmax, with the latter providing the most accurate and stable results. Throughout the majority of the paper, we fix the description of the nonlinear matter evolution using a hybrid perturbative-N-body approach, RESPRESSO, that was found in the first paper to be the closest performing to the measured matter spectrum. We also test the impact of different modeling assumptions based on perturbative approaches, such as galilean-invariant Renormalised Perturbation Theory (gRPT) and effective field theory (EFT). In all cases, we find the inclusion of scale-dependent shot noise to increase the range of validity of the model in terms of FoB and χ2. Interestingly, these model extensions with additional free parameters do not necessarily lead to an increase in the maximally achievable FoM for the cosmological parameters ðh; Ωch2; AsÞ, which are generally consistent with those of the simpler model at smaller kmax

Redshift-space distortions from the cross-correlation of photometric populations

Several papers have recently highlighted the possibility of measuring redshift-space distortions from angular auto-correlations of galaxiesin photometric redshift bins. In this work, we extend this idea to include as observables the cross-correlations between redshift bins, as an additional way of measuring radial information. We show that this extra information allows us to reduce the recovered error in the growth rate index γ by a factor of ~2. Although the final error in γ depends on the bias and the mean photometric accuracy of the galaxy sample, the improvement from adding cross-correlations is robust in different settings. Another factor of 2-3 improvement in the determination of γ can be achieved by considering two galaxy populations over the same photometric sky area but with different biases. This additional gain is shown to be much larger than the one from the same populations when observed over different areas of the sky (with twice the combined area). The total improvement of ~5 implies that a photometric survey such as Dark Energy Survey should be able to recover γ at the 5-10 per cent from the angular clustering in linear scales of two different tracers. It can also constrain the evolution of f(z) × σ8(z) in few bins beyond z ~ 0.8-0.9 at the 10-15 per cent level per bin, compatible with recent constrains from lower z spectroscopic surveys. We also show how further improvement can be achieved by reducing the photometric redshift error