The non-linear redshift-space power spectrum of galaxies

We study the power spectrum of galaxies in redshift space, with third-order perturbation theory to include corrections that are absent in linear theory. We assume a local bias for the galaxies: i.e., the galaxy density is sampled from some local function of the underlying mass distribution. We find that the effect of the non-linear bias in real space is to introduce two new features: first, there is a contribution to the power which is constant with wavenumber, whose nature we reveal as essentially a shot-noise term. In principle this contribution can mask the primordial power spectrum, and could limit the accuracy with which the latter might be measured on very large scales. Secondly, the effect of second- and third-order bias is to modify the effective bias (defined as the square root of the ratio of galaxy power spectrum to matter power spectrum). The effective bias is almost scale-independent over a wide range of scales. These general conclusions also hold in redshift space. In addition, we have investigated the distortion of the power spectrum by peculiar velocities, which may be used to constrain the density of the Universe. We look at the quadrupole-to-monopole ratio, and find that higher order terms can mimic linear theory bias, but the bias implied is neither the linear bias, nor the effective bias referred to above. We test the theory with biased N -body simulations, and find excellent agreement in both real and redshift space, providing the local biasing is applied on a scale whose fractional rms density fluctuations are < 0.5.     

Scale-dependent galaxy bias in the Sloan Digital Sky Survey as a function of luminosity and colour

It has been known for a long time that the clustering of galaxies changes as a function of galaxy type. This galaxy bias acts as a hindrance to the extraction of cosmological information from the galaxy power spectrum or correlation function. Theoretical arguments show that a change in the amplitude of the clustering between galaxies and mass on large scales is unavoidable, but cosmological information can be easily extracted from the shape of the power spectrum or correlation function if this bias is independent of scale. Scale-dependent bias is generally small on large scales,   k < 0.1  h  Mpc⁻¹  , but on smaller scales can affect the recovery of  Ω_m h   from the measured shape of the clustering signal, and have a small effect on the Baryon Acoustic Oscillations. In this paper, we investigate the transition from scale-independent to scale-dependent galaxy bias as a function of galaxy population. We use the Sloan Digital Sky Survey Data Release 5 sample to fit various models, which attempt to parametrize the turn-off from scale-independent behaviour. For blue galaxies, we find that the strength of the turn-off is strongly dependent on galaxy luminosity, with stronger scale-dependent bias on larger scales for more luminous galaxies. For red galaxies, the scale dependence is a weaker function of luminosity. Such trends need to be modelled in order to optimally extract the information available in future surveys, and can help with the design of such surveys.     

The three-dimensional power spectrum of dark and luminous matter from the VIRMOS-DESCART cosmic shear survey

We present the first optimal power spectrum estimation and three-dimensional deprojections for the dark and luminous matter and their cross-correlations. The results are obtained using a new optimal fast estimator, deprojected using minimum variance and Singular Value Decomposition (SVD) techniques. We show the resulting 3D power spectra for dark matter and galaxies, and their covariance for the VIRMOS-DESCART weak lensing shear and galaxy data. The survey is most sensitive to non-linear scales   k _NL∼ 1 h Mpc⁻¹  . On these scales, our 3D power spectrum of dark matter is in good agreement with the RCS 3D power spectrum found by Tegmark & Zaldarriaga. Our galaxy power is similar to that found by the 2MASS survey, and larger than that of SDSS, APM and RCS, consistent with the expected difference in galaxy population.
We find an average bias   b = 1.24 ± 0.18  for the I -selected galaxies, and a cross-correlation coefficient   r = 0.75 ± 0.23  . Together with the power spectra, these results optimally encode the entire two point information about dark matter and galaxies, including galaxy–galaxy lensing. We address some of the implications regarding galaxy haloes and mass-to-light ratios. The best-fitting 'halo' parameter   h ≡ r / b = 0.57 ± 0.16  , suggesting that dynamical masses estimated using galaxies systematically underestimate total mass.
Ongoing surveys, such as the Canada–France–Hawaii Telescope Legacy Survey, will significantly improve on the dynamic range, and future photometric redshift catalogues will allow tomography along the same principles.     

Galaxy–dark matter correlations applied to galaxy–galaxy lensing: predictions from the semi-analytic galaxy formation models

We use semi-analytic models of galaxy formation combined with high-resolution N -body simulations to make predictions for galaxy–dark matter correlations and apply them to galaxy–galaxy lensing. We analyse cross-power spectra between the dark matter and different galaxy samples selected by luminosity, colour or star formation rate. We compare the predictions with the recent detection by the Sloan Digital Sky Survey (SDSS). We show that the correlation amplitude and the mean tangential shear depend strongly on the luminosity of the sample on scales below 1  h ⁻¹ Mpc, reflecting the correlation between the galaxy luminosity and the halo mass. The cross-correlation cannot, however, be used to infer the halo profile directly because different halo masses dominate on different scales and because not all galaxies are at the centres of the corresponding haloes. We compute the redshift evolution of the cross-correlation amplitude and compare it with those of galaxies and dark matter. We also compute the galaxy–dark matter correlation coefficient and show that it is close to unity on scales above 1  h ⁻¹ Mpc for all considered galaxy types. This would allow one to extract the bias and the dark matter power spectrum on large scales from the galaxy and galaxy–dark matter correlations.     

Peaks in the cosmological density field: sensitivity to initial power spectrum, redshift distortions and galaxy halo occupation

We investigate the number density of maxima in the cosmological galaxy density field smoothed with a filter as a probe of clustering. In previous work it has been shown that this statistic is closely related to the slope of the linear power spectrum, even when the directly measured power spectrum is non-linear. In the present paper we investigate the sensitivity of the peak number density to various models with differing power spectra, including rolling index models, cosmologies with massive neutrinos and different baryon densities. We find that rolling index models which have given an improved fit to CMB/LSS (cosmic microwave background/large scale structure) data yield a ∼10 per cent difference in peak density compared to the scale invariant case. Models with 0.3 eV neutrinos have effects of similar magnitude and it should be possible to constrain them with data from current galaxy redshift surveys. Baryon oscillations in the power spectrum also give rise to distinctive features in the peak density. These are preserved without modification when measured from the peak density in fully non-linear N -body simulations. Using the simulations, we also investigate how the peak density is modified in the presence of redshift distortions. Redshift distortions cause a suppression of the number of peaks, largely due to fingers of God overlapping in redshift space. We find that this effect can be modelled by using a modification of the input power spectrum. We also study the results when the simulation density field is traced by galaxies obtained by populating haloes with a halo occupation distribution consistent with observations. The peak number density is consistent with that in the dark matter for filter scales  >4  h ⁻¹ Mpc  , for which we find good agreement with the linear theory predictions. In a companion paper we analyse data from the 2dF Galaxy Redshift Survey.     

Tidal alignments as a contaminant of redshift space distortions

We investigate the effect of orientation-dependent selection effects on galaxy clustering in redshift space. It is found that if galaxies are aligned by large-scale tidal fields, then these selection effects give rise to a dependence of the observed galaxy density on the local tidal field, in addition to the well-known dependences on the matter density and radial velocity gradient. This alters the galaxy power spectrum in a way that is different for Fourier modes parallel to and perpendicular to the line of sight. These tidal galaxy alignments can thus mimic redshift space distortions (RSD), and thus result in a bias in the measurement of the velocity power spectrum. If galaxy orientations are affected only by the local tidal field, then the tidal alignment effect has exactly the same scale and angular dependence as the RSDs in the linear regime, so it cannot be projected out or removed by masking small scales in the analysis. We consider several toy models of tidal alignments and orientation-dependent selection, normalize their free parameter (an amplitude) to recent observations, and find that they could bias the velocity amplitude   f ( z ) G ( z )  by 5–10 per cent in some models, although most models give much smaller contamination. We conclude that tidal alignments may be a significant systematic error in RSD measurements that aim to test general relativity via the growth of large-scale structure. We briefly discuss possible mitigation strategies.     

Determination of the linear mass power spectrum from the mass function of galaxy clusters

We develop a new method to determine the linear mass power spectrum using the mass function of galaxy clusters. We obtain the rms mass fluctuation  σ( M )  using the expression for the mass function in the Press & Schechter, Sheth, Mo & Tormen and Jenkins et al. formalisms. We apply different techniques to recover the adimensional power spectrum  Δ²( k )  from  σ( M )  namely the   k _eff  approximation, the singular value decomposition and the linear regularization method. The application of these techniques to the τCDM and ΛCDM GIF simulations shows a high efficiency in recovering the theoretical power spectrum over a wide range of scales. We compare our results with those derived from the power spectrum of the spatial distribution of the same sample of clusters in the simulations obtained by application of the classical Feldman, Kaiser & Peacock (FKP) method. We find that the mass function based method presented here can provide a very accurate estimate of the linear power spectrum, particularly for low values of k . This estimate is comparable to, or even better behaved than, the FKP solution.
The principal advantage of our method is that it allows the determination of the linear mass power spectrum using the joint information of objects of a wide range of masses without dealing with specific assumptions on the bias relative to the underlying mass distribution.     

SDSS galaxy clustering: luminosity and colour dependence and stochasticity

Differences in clustering properties between galaxy subpopulations complicate the cosmological interpretation of the galaxy power spectrum, but can also provide insights about the physics underlying galaxy formation. To study the nature of this relative clustering, we perform a counts-in-cells analysis of galaxies in the Sloan Digital Sky Survey in which we measure the relative bias between pairs of galaxy subsamples of different luminosities and colours. We use a generalized  χ²  test to determine if the relative bias between each pair of subsamples is consistent with the simplest deterministic linear bias model, and we also use a maximum likelihood technique to further understand the nature of the relative bias between each pair. We find that the simple, deterministic model is a good fit for the luminosity-dependent bias on scales above  ∼2  h ⁻¹ Mpc  , which is good news for using magnitude-limited surveys for cosmology. However, the colour-dependent bias shows evidence for stochasticity and/or non-linearity which increases in strength towards smaller scales, in agreement with previous studies of stochastic bias. Also, confirming hints seen in earlier work, the luminosity-dependent bias for red galaxies is significantly different from that of blue galaxies: both luminous and dim red galaxies have higher bias than moderately bright red galaxies, whereas the biasing of blue galaxies is not strongly luminosity dependent. These results can be used to constrain galaxy formation models and also to quantify how the colour and luminosity selection of a galaxy survey can impact measurements of the cosmological matter power spectrum.     

Clustering of galaxy clusters in cold dark matter universes

We use very large cosmological N -body simulations to obtain accurate predictions for the two-point correlations and power spectra of mass-limited samples of galaxy clusters. We consider two currently popular cold dark matter (CDM) cosmogonies, a critical density model ( τ CDM) and a flat low density model with a cosmological constant (ΛCDM). Our simulations each use 10⁹ particles to follow the mass distribution within cubes of side 2  h ⁻¹ Gpc ( τ CDM) and 3  h ⁻¹ Gpc (ΛCDM) with a force resolution better than 10⁻⁴ of the cube side. We investigate how the predicted cluster correlations increase for samples of increasing mass and decreasing abundance. Very similar behaviour is found in the two cases. The correlation length increases from     for samples with mean separation     to     for samples with     The lower value here corresponds to τ CDM and the upper to ΛCDM. The power spectra of these cluster samples are accurately parallel to those of the mass over more than a decade in scale. Both correlation lengths and power spectrum biases can be predicted to better than 10 per cent using the simple model of Sheth, Mo & Tormen. This prediction requires only the linear mass power spectrum and has no adjustable parameters. We compare our predictions with published results for the automated plate measurement (APM) cluster sample. The observed variation of correlation length with richness agrees well with the models, particularly for ΛCDM. The observed power spectrum (for a cluster sample of mean separation     ) lies significantly above the predictions of both models.     

A primordial feature at the scale of superclusters of galaxies

We investigate a spatially flat cold dark matter model (with the matter density parameter     with a primordial feature in the initial power spectrum. We assume that there is a bump in the power spectrum of density fluctuations at wavelengths     , which corresponds to the scale of superclusters of galaxies . There are indications for such a feature in the power spectra derived from redshift surveys and also in the power spectra derived from peculiar velocities of galaxies. We study the mass function of clusters of galaxies, the power spectrum of the cosmic microwave background (CMB) temperature fluctuations, the rms bulk velocity and the rms peculiar velocity of clusters of galaxies. The baryon density is assumed to be consistent with the big bang nucleosynthesis value. We show that with an appropriately chosen feature in the power spectrum of density fluctuations at the scale of superclusters, the mass function of clusters, the CMB power spectrum, the rms bulk velocity and the rms peculiar velocity of clusters are in good agreement with the observed data.     

The amplitude of mass density fluctuations at z 3.25 from the Ly forest of Q1422+231

The real-space optical-depth distribution along the line of sight to the QSO Q1422+231 is recovered from two HIRES spectra using a modified version of the inversion method proposed by Nusser & Haehnelt. The first two moments of the truncated optical-depth distribution are used to constrain the density-fluctuation amplitude of the intergalactic medium (IGM) assuming that the IGM is photoionized by a metagalactic UV background and obeys a temperaturedensity relation. The fluctuation amplitude and the power-law index of the relation between gas and neutral hydrogen (H  i ) density are degenerate. The rms of the IGM density at z 3.25 estimated from the first spectrum is with 1.56< <2 for plausible reionization histories. This corresponds to 0.9 2.1 with ( =1.7)=1.44±0.3. The values obtained from the second spectrum are higher by 20 per cent. If the IGM density traces the dark matter (DM) as suggested by numerical simulations we have measured the fluctuation amplitude of the DM density at an effective Jeans scale of a few 100 kpc. For cold dark matter (CDM)-like power spectra the amplitude of dark matter fluctuations on these small scales depends on the cosmological density parameter . For power spectra normalized to reproduce the space density of present-day clusters and with a slope parameter of =0.21 consistent with the observed galaxy power spectrum, the inferred can be expressed as: =0.61( /1.7)^1.3( x _J/0.62)^0.6 for a flat universe, and =0.91( /1.7)^1.3( x _J/0.62)^0.7 for a =0 universe. x _J is the effective Jeans scale in (comoving) h ¹ Mpc. Based on a suite of detailed mock spectra the 1 error is 25 per cent. The estimates increase with increasing . For the second spectrum we obtain 15 per cent lower values.     

Testing deprojection algorithms on mock angular catalogues: evidence for a break in the power spectrum

We produce mock angular catalogues from simulations with different initial power spectra to test methods that recover measures of clustering in three dimensions, such as the power spectrum, variance and higher order cumulants. We find that the statistical properties derived from the angular mock catalogues are in good agreement with the intrinsic clustering in the simulations. In particular, we concentrate on the detailed predictions for the shape of the power spectrum, P ( k ). We find that there is good evidence for a break in the galaxy P ( k ) at scales in the range 0.02< k <0.06 h Mpc⁻¹, using an inversion technique applied to the angular correlation function measured from the APM Galaxy Survey. For variants on the standard cold dark matter (CDM) model, a fit at the location of the break implies Ω h =0.45±0.10, where Ω is the ratio of the total matter density to the critical density, and Hubble's constant is parametrized as H ₀=100 h km s⁻¹ Mpc⁻¹. On slightly smaller, though still quasi-linear scales, there is a feature in the APM power spectrum where the local slope changes appreciably, with the best match to CDM models obtained for Ω h ≃0.2. Hence the location and narrowness of the break in the APM power spectrum combined with the rapid change in its slope on quasi-linear scales cannot be matched by any variant of CDM, including models that have a non-zero cosmological constant or a tilt to the slope of the primordial P ( k ). These results are independent of the overall normalization of the CDM models or any simple bias that exists betwen the galaxy and mass distributions.     

Cosmological implications of the PSCz PDF and its moments

We compare the probability density function (PDF) and its low-order moments (variance and skewness) of the smoothed IRAS Point Source Catalogue Redshift Survey (PSC z ) galaxy density field and of the corresponding simulated PSC z look-alikes, generated from N -body simulations of six different dark matter models: four structure-normalized with     and     , one COBE -normalized, and the old standard cold dark matter model. The galaxy distributions are smoothed with a Gaussian window at three different smoothing scales,     , 10 and 15  h ⁻¹ Mpc. We find that the simulation PSC z look-alike PDFs are sensitive only to the normalization of the power spectrum, probably owing to the shape similarity of the simulated galaxy power spectrum on the relevant scales. We find that the only models that are consistent, at a high significance level, with the observed PSC z PDF are models with a relatively low power spectrum normalization     . From the phenomenologically derived σ ₈–moments relation, fitted from the simulation data, we find that the PSC z moments suggest     .     

Scale-dependent bias induced by local non-Gaussianity: a comparison to N-body simulations

We investigate the effect of primordial non-Gaussianity of the local f _NL type on the auto- and cross-power spectra of dark matter haloes using simulations of the Λ cold dark matter cosmology. We perform a series of large N -body simulations of both positive and negative f _NL, spanning the range between 10 and 100. Theoretical models predict a scale-dependent bias correction  Δ b ( k , f _NL)  that depends on the linear halo bias   b ( M )  . We measure the power spectra for a range of halo mass and redshifts covering the relevant range of existing galaxy and quasar populations. We show that auto- and cross-correlation analyses of bias are consistent with each other. We find that for low wavenumbers with   k < 0.03  h  Mpc⁻¹  the theory and the simulations agree well with each other for biased haloes with   b ( M ) > 1.5  . We show that a scale-independent bias correction improves the comparison between theory and simulations on smaller scales, where the scale-dependent effect rapidly becomes negligible. The current limits on f _NL from Slosar et al. come mostly from very large scales   k < 0.01  h  Mpc⁻¹  and, therefore, remain valid. For the halo samples with   b ( M ) < 1.5 − 2  , we find that the scale-dependent bias from non-Gaussianity actually exceeds the theoretical predictions. Our results are consistent with the bias correction scaling linearly with f _NL.     

Power spectrum analysis of the ESO Slice Project galaxy redshift survey

We measure the power spectrum of the galaxy distribution in the ESO Slice Project (ESP) galaxy redshift survey. We develop a technique to describe the survey window function analytically, and then deconvolve it from the measured power spectrum using a variant of the Lucy method. We test the whole deconvolution procedure on ESP mock catalogues drawn from large N -body simulations, and find that it is reliable for recovering the correct amplitude and shape of P ( k ) at k >0.065  h  Mpc⁻¹. In general, the technique is applicable to any survey composed of a collection of circular fields with an arbitrary pattern on the sky, as typical of surveys based on fibre spectrographs. The estimated power spectrum has a well-defined power-law shape k ⁿ with n ≃−2.2 for k ≥0.2  h  Mpc⁻¹, and a smooth bend to a flatter shape ( n ≃−1.6) for smaller k . The smallest wavenumber where a meaningful reconstruction can be performed ( k ∼0.06  h  Mpc⁻¹) does not allow us to explore the range of scales where other power spectra seem to show a flattening and hint at a turnover. We also find, by a direct comparison of the Fourier transforms, that the estimate of the two-point correlation function ξ ( s ) is much less sensitive to the effect of a problematic window function, such as that of the ESP, than the power spectrum. Comparison with other surveys shows an excellent agreement with estimates from blue-selected surveys. In particular, the ESP power spectrum is virtually indistinguishable from that of the Durham–UKST survey over the common range of k , an indirect confirmation of the quality of the deconvolution technique applied.     

Higher order moments of the density field in a parametrized sequence of non-Gaussian theories

We calculate the higher order moments in a sequence of models where the initial density fluctuations are drawn from a     distribution with a power-law power spectrum. For large values of , the distribution is approximately Gaussian, and we reproduce the values known from perturbation theory. As is lowered the distribution becomes progressively more non-Gaussian, approximating models with rare, high-amplitude peaks. The limit   =1  is a realization of recently proposed isocurvature models for producing early galaxy formation, where the density perturbations are quadratic in a Gaussian field.     

The clustering of galaxy clusters in cosmological models with non-Gaussian initial conditions: predictions for future surveys

We predict the biasing and clustering properties of galaxy clusters that are expected to be observed in the catalogues produced by two forthcoming X-ray and Sunyaev–Zel'dovich effect surveys. We study a set of flat cosmological models where the primordial density probability distribution shows deviations from Gaussianity in agreement with current observational bounds form the background radiation. We consider both local and equilateral shapes for the primordial bispectrum in non-Gaussian models. The two catalogues investigated are those produced by the e ROSITA wide survey and from a survey based on South Pole Telescope observations. It turns out that both the bias and observed power spectrum of galaxy clusters are severely affected in non-Gaussian models with local shape of the primordial bispectrum, especially at large scales. On the other hand, models with equilateral shape of the primordial bispectrum show only a mild effect at all scales, that is difficult to be detected with clustering observations. Between the two catalogues, the one performing better is the e ROSITA one, since it contains only the largest masses that are more sensitive to primordial non-Gaussianity.     

Redshift-space distortions in the PSCz galaxy catalogue

We apply a spherical harmonic analysis to the Point Source Redshift Survey (PSC z ), to compute the real-space galaxy power spectrum and the degree of redshift distortion caused by peculiar velocities. We employ new parameter eigenvector and hierarchical data compression techniques, allowing a much larger number of harmonic modes to be included, and correspondingly smaller error bars. Using 4644 harmonic modes, compressed to 2278, we find that the IRAS redshift-space distortion parameter is     and the amplitude of galaxy clustering on a scale of     is     . Combining these we find the amplitude of mass perturbations is     . While this is compatible with results from the cosmic microwave background (CMB), with a small degree of tilt, it disagrees with the amplitude of matter perturbations estimated from the abundance of clusters by a factor of 2, independent of cosmology. A preliminary model fitting analysis combining the CMB with either the PSC z or cluster abundances shows that the cosmological matter density parameter     , and the IRAS bias parameter     . However, the cluster abundances suggest that     and     , while the PSC z requires     and     . Given the physics of galaxy formation is poorly constrained, we conclude that IRAS galaxies and mass are only partially correlated.     

The power spectrum of rich clusters of galaxies on large spatial scales

We present an analysis of the redshift-space power spectrum, P ( k ), of rich clusters of galaxies based on an automated cluster catalogue selected from the APM Galaxy Survey. We find that P ( k ) can be approximated by a power law, P ( k )∝ kⁿ , with n ≈−1.6 over the wavenumber range 0.04< k <0.1 h Mpc⁻¹. Over this range of wavenumbers, the APM cluster power spectrum has the same shape as the power spectra measured for optical and IRAS galaxies. This is consistent with a simple linear bias model in which different tracers have the same power spectrum as that of the mass distribution, but shifted in amplitude by a constant biasing factor. On larger scales, the power spectrum of APM clusters flattens and appears to turn over on a scale k ∼0.03 h Mpc⁻¹. We compare the power spectra estimated from simulated APM cluster catalogues with those estimated directly from cubical N -body simulation volumes, and find that the APM cluster survey should give reliable estimates of the true power spectrum at wavenumbers k ≳0.02 h Mpc⁻¹. These results suggest that the observed turnover in the power spectrum may be a real feature of the cluster distribution, and that we have detected the transition to a near-scale-invariant power spectrum implied by observations of anisotropies in the cosmic microwave background radiation. The scale of the turnover in the cluster power spectrum is in good agreement with the scale of the turnover observed in the power spectrum of APM galaxies.