Redshift-space Power Spectrum Research Articles

Redshift-space distortions in massive neutrino and evolving dark energy cosmologies

Large-scale structure surveys in the coming years will measure the redshift-space power spectrum to unprecedented accuracy, allowing for powerful new tests of the LambdaCDM picture as well as measurements of particle physics parameters such as the neutrino masses. We extend the Time-RG perturbative framework to redshift space, computing the power spectrum P_s(k,mu) in massive neutrino cosmologies with time-dependent dark energy equations of state w(z). Time-RG is uniquely capable of incorporating scale-dependent growth into the P_s(k,mu) computation, which is important for massive neutrinos as well as modified gravity models. Although changes to w(z) and the neutrino mass fraction both affect the late-time scale-dependence of the non-linear power spectrum, we find that the two effects depend differently on the line-of-sight angle mu. Finally, we use the HACC N-body code to quantify errors in the perturbative calculations. For a LambdaCDM model at redshift z=1, our procedure predicts the monopole~(quadrupole) to 1% accuracy up to a wave number 0.19h/Mpc (0.28h/Mpc), compared to 0.08h/Mpc (0.07h/Mpc) for the Kaiser approximation and 0.19h/Mpc (0.16h/Mpc) for the current state-of-the-art perturbation scheme. Our calculation agrees with the simulated redshift-space power spectrum even for neutrino masses above the current bound, and for rapidly-evolving dark energy equations of state, |dw/dz| ~ 1. Along with this article, we make our redshift-space Time-RG implementation publicly available as the code redTime.

Effective theory of dark energy at redshift survey scales

We explore the phenomenological consequences of general late-time modifications of gravity in the quasi-static approximation, in the case where cold dark matter is non-minimally coupled to the gravitational sector. Assuming spectroscopic and photometric surveys with configuration parameters similar to those of the Euclid mission, we derive constraints on our effective description from three observables: the galaxy power spectrum in redshift space, tomographic weak-lensing shear power spectrum and the correlation spectrum between the integrated Sachs-Wolfe effect and the galaxy distribution. In particular, with ΛCDM as fiducial model and a specific choice for the time dependence of our effective functions, we perform a Fisher matrix analysis and find that the unmarginalized 68% CL errors on the parameters describing the modifications of gravity are of order σ∼10−2–10−3. We also consider two other fiducial models. A nonminimal coupling of CDM enhances the effects of modified gravity and reduces the above statistical errors accordingly. In all cases, we find that the parameters are highly degenerate, which prevents the inversion of the Fisher matrices. Some of these degeneracies can be broken by combining all three observational probes.

Galaxy power spectrum in redshift space: Combining perturbation theory with the halo model

Theoretical modeling of the redshift-space power spectrum of galaxies is crucially important to correctly extract cosmological information from redshift surveys. The task is complicated by the nonlinear biasing and redshift space distortion effects, which change with halo mass, and by the wide distribution of halo masses and their occupations by galaxies. One of the main modeling challenges is the existence of satellite galaxies that have both radial distribution and large virial velocities inside halos, a phenomenon known as the Finger-of-God effect. We present a model for the galaxy power spectrum of in which we decompose a given galaxy sample into central and satellite galaxies and relate different contributions to 1- and 2-halo terms in a halo model. Our primary goal is to ensure that any parameters that we introduce have physically meaningful values, and are not just fitting parameters. For the 2-halo terms we use the previously developed RSD modeling of halos in the context of distribution function and perturbation theory approach. This term needs to be multiplied by the effect of radial distances and velocities of satellites inside the halo. To this one needs to add the 1-halo terms, which are non-perturbative. We show that the real space 1-halo terms can be modeled as almost constant, with the finite extent of the satellites inside the halo inducing a small k^2R^2 term, where R is related to the size of the halo. We adopt a similar model for FoG in redshift space, ensuring that FoG velocity dispersion is related to the halo mass. For FoG k^2 type expansions do not work and FoG resummation must be used instead. We test several damping functions to model the velocity dispersion FoG effect. Applying the formalism to mock galaxies modeled after the "CMASS" sample of the BOSS survey, we find that our predictions for the redshift-space power spectra are accurate up to k~0.4Mpc/h.

Fingers-of-God effect of infalling satellite galaxies

Abstract Non-linear redshift-space distortion known as the Fingers-of-God (FoG) effect is a major systematic uncertainty in redshift-space distortion studies conducted to test gravity models. The FoG effect has been usually attributed to the random motion of galaxies inside their clusters. When the internal galaxy motion is not well virialized, however, the coherent infalling motion towards the cluster centre generates the FoG effect. Here, we derive an analytical model of the satellite velocity distribution due to the infall motion combined with the random motion. We show that the velocity distribution becomes far from Maxwellian when the infalling motion is dominant. We use simulated subhalo catalogues to find that the contribution of infall motion is important to massive subhaloes and that the velocity distribution has a top-hat like shape as expected from our analytic model. We also study the FoG effect due to infall motion on the redshift-space power spectrum. Using simulated mock samples of luminous red galaxies constructed from haloes and massive subhaloes in N-body simulations, we show that the redshift-space power spectra can differ from expectations when the infall motion is ignored.

The power spectrum and bispectrum of SDSS DR11 BOSS galaxies – II. Cosmological interpretation

We examine the cosmological implications of the measurements of the linear growth rate of cosmological structure obtained in a companion paper from the power spectrum and bispectrum monopoles of the Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey Data, Release 11, CMASS galaxies. This measurement was of $f^{0.43}\sigma_8$, where $\sigma_8$ is the amplitude of dark matter density fluctuations, and $f$ is the linear growth rate, at the effective redshift of the survey, $z_{\rm eff}=0.57$. In conjunction with Cosmic Microwave Background (CMB) data, interesting constraints can be placed on models with non-standard neutrino properties and models where gravity deviates from general relativity on cosmological scales. In particular, the sum of the masses of the three species of the neutrinos is constrained to $m_\nu<0.49\,{\rm eV}$ (at 95\% confidence level) when the $f^{0.43}\sigma_8$ measurement is combined with state-of-the-art CMB measurements. Allowing the effective number of neutrinos to vary as a free parameter does not significantly change these results. When we combine the measurement of $f^{0.43}\sigma_8$ with the complementary measurement of $f\sigma_8$ from the monopole and quadrupole of the two-point correlation function we are able to obtain an independent measurements of $f$ and $\sigma_8$. We obtain $f=0.63\pm0.16$ and $\sigma_8=0.710\pm0.086$ (68\% confidence level). This is the first time when these parameters have been able to be measured independently using the redshift-space power spectrum and bispectrum measurements from galaxy clustering data only.

Dark matter and halo bispectrum in redshift space: theory and applications

We present a phenomenological modification of the standard perturbation theory prediction for the bispectrum in redshift space that allows us to extend the model to mildly non-linear scales over a wide range of redshifts, z⩽1.5. Our model require 18 free parameters that are fitted to N-body simulations using the shapes k2/k1=1, 1.5, 2.0, 2.5. We find that we can describe the bispectrum of dark matter particles with ∼5% accuracy for ki≲0.10 h/Mpc at z=0, for ki≲0.15 h/Mpc at z=0.5, for ki≲0.17 h/Mpc at z=1.0 and for ki≲0.20 h/Mpc at z=1.5. For very squeezed triangles with k1=k2≳0.1 hMpc-1 and k3⩽0.02 hMpc-1, however, neither SPT nor the proposed fitting formula are able to describe the measured dark matter bispectrum with this accuracy. We show that the fitting formula is sufficiently general that can be applied to other intermediate shapes such as k2/k1=1.25, 1.75, and 2.25. We also test that the fitting formula is able to describe with similar accuracy the bispectrum of cosmologies with different Ωm, in the range 0.2≲ Ωm ≲ 0.4, and consequently with different values of the logarithmic grow rate f at z=0, 0.4≲ f(z=0) ≲ 0.6. We apply this new formula to recover the bias parameters, f and σ8, by combining the redshift space power spectrum monopole and quadrupole with the bispectrum monopole for both dark matter particles and haloes. We find that the combination of these three statistics can break the degeneracy between b1, f and σ8. For dark matter particles the new model can be used to recover f and σ8 with ∼1% accuracy. For dark matter haloes we find that f and σ8 present larger systematic shifts, ∼10%. The systematic offsets arise because of limitations in the modelling of the interplay between bias and redshift space distortions, and represent a limitation as the statistical errors of forthcoming surveys reach this level. Conveniently, we find that these residual systematics are mitigated for combinations of parameters. In particular, the quantity fσ8 is still recovered with ∼1% accuracy for the particular halo population and cosmology studied. The improvement on the modelling of the bispectrum presented in this paper will be useful for extracting information from current and future galaxy surveys.

Remapping simulated halo catalogues in redshift space

We discuss the extension to redshift space of a rescaling algorithm, designed to alter the effective cosmology of a pre-existing simulated particle distribution or catalogue of dark matter haloes. The rescaling approach was initially developed by Angulo & White and was adapted and applied to halo catalogues in real space in our previous work. This algorithm requires no information other than the initial and target cosmological parameters, and it contains no tuned parameters. It is shown here that the rescaling method also works well in redshift space, and that the rescaled simulations can reproduce the growth rate of cosmological density fluctuations appropriate for the target cosmology. Even when rescaling a grossly non-standard model with Λ = 0 and zero baryons, the redshift-space power spectrum of standard Λ cold dark matter can be reproduced to about 5 per cent error for k < 0.2 h Mpc−1. The ratio of quadrupole-to-monopole power spectra remains correct to the same tolerance up to k = 1 h Mpc−1, provided that the input halo catalogue contains measured internal velocity dispersions.

Peculiar velocities in redshift space: formalism, N-body simulations and perturbation theory

Direct measurements of peculiar velocities of galaxies and clusters of galaxies can in principle provide explicit information on the three dimensional mass distribution, but this information is modulated by the fact that velocity field is sampled at galaxy positions, and is thus probing galaxy momentum. We derive expressions for the cross power spectrum between the density and momentum field and the auto spectrum of the momentum field in redshift space, by extending the distribution function method to these statistics. The resulting momentum cross and auto power spectra in redshift space are expressed as infinite sums over velocity moment correlators in real space, as is the case for the density power spectrum in redshift space. We compute each correlator using Eulerian perturbation theory (PT) and halo biasing model and compare the resulting redshift-space velocity statistics to those measured from N-body simulations for both dark matter and halos. We find that in redshift space linear theory predictions for the density-momentum cross power spectrum as well as for the momentum auto spectrum fail to predict the N-body results at very large scales. On the other hand, our nonlinear PT prediction for these velocity statistics, together with real-space power spectrum for dark matter from simulations, improves the accuracy for both dark matter and halos. We also present the same analysis in configuration space, computing the redshift-space pairwise mean infall velocities and velocity correlation function and compare to nonlinear PT.

Semi-analytical study on the generic degeneracy for galaxy clustering measurements

AbstractFrom the galaxy power spectrum in redshift space, we derive semi-analytical results on the generic degeneracy of galaxy clustering measurements. Defining the observables A = Gbσ8 and R = Gfσ8, (being G the growth function, b the bias, f the growth rate, and σ8 the amplitude of the power spectrum), we perform a Fisher matrix formalism to forecast the expected precision of these quantities for a Euclid-like survey. Among the results we found that galaxy surveys have generically a slightly negative correlation between A and R, and they can always measure R about 3.7 to 4.7 times better than A.

Beyond consistency test of gravity with redshift-space distortions at quasilinear scales

Redshift-space distortions (RSD) offer an attractive method to measure the growth of cosmic structure on large scales, and combining with the measurement of the cosmic expansion history, it can be used as cosmological tests of gravity. With the advent of future galaxy redshift surveys aiming at precisely measuring the RSD, an accurate modeling of RSD going beyond linear theory is a critical issue in order to detect or disprove small deviations from general relativity (GR). While several improved models of RSD have been recently proposed based on the perturbation theory (PT), the framework of these models heavily relies on GR. Here, we put forward a new PT prescription for RSD in general modified gravity models. As a specific application, we present theoretical predictions of the redshift-space power spectra in f(R) gravity model, and compare them with N-body simulations. Using the PT template that takes into account the effects of both modifications of gravity and RSD properly, we successfully recover the fiducial model parameter in N-body simulations in an unbiased way. On the other hand, we found it difficult to detect the scale dependence of the growth rate in a model-independent way based on GR templates.

Using the topology of large-scale structure in the WiggleZ Dark Energy Survey as a cosmological standard ruler

We present new and accurate measurements of the cosmic distance–redshift relation, spanning 0.2 < z < 1, using the topology of large-scale structure as a cosmological standard ruler. Our results derive from an analysis of the Minkowski functionals of the density field traced by the WiggleZ Dark Energy Survey. The Minkowski functionals are a set of statistics which completely describe the topological nature of each isodensity surface within the field, as a function of the density value. Given the shape of the underlying matter power spectrum, measured by fluctuations in the cosmic microwave background radiation, the expected amplitudes of the Minkowski functionals are specified as an excursion set of a Gaussian random field, with minimal non-Gaussian corrections for the smoothing scales ≥10 h−1 Mpc considered in this analysis. The measured amplitudes then determine the cosmic distance DV(z), which we obtain with 3–7 per cent accuracies in six independent redshift slices, with the standard ruler originating in the known curvature of the model power spectrum at the smoothing scale. We introduce a new method for correcting the topological statistics for the sparse-sampling of the density field by the galaxy tracers, and validate our overall methodology using mock catalogues from N-body simulations. Our distance measurements are consistent with standard models which describe the cosmic expansion history, and with previous analyses of baryon acoustic oscillations (BAOs) detected by the WiggleZ Survey, with the topological results yielding a higher distance precision by a factor of 2. However, the full redshift-space power-spectrum shape is required to recover the topological distances, in contrast to the preferred length scale imprinted by BAOs, which is determined by simpler physics.

Where are the Luminous Red Galaxies (LRGs)? Using correlation measurements and lensing to relate LRGs to dark matter haloes

Nonlinear redshift-space distortions, the Finger-of-God (FoG) effect, can complicate the interpretation of the galaxy power spectrum. Here, we demonstrate the method proposed by Hikage et al. (2012) to use complimentary observations to directly constrain this effect on the data. We use catalogs of Luminous Red Galaxies (LRGs) and photometric galaxies from the SDSS DR7 to measure the redshift-space power spectrum of LRGs, the cross-correlation of LRGs with the shapes of background photometric galaxies (galaxy-galaxy weak lensing), and the projected cross-correlation of LRGs with photometric galaxies having similar photometric redshifts to the LRG spectroscopic redshift. All of these measurements use a reconstructed halo field. While we use the position of each LRG for single LRG systems, we compare the measurements using different halo-center proxies for multiple-LRG systems (4.5 per cent of all the halos): the brightest LRG position (BLRG), the faintest LRG position (FLRG) and their arithmetical mean position (Mean), respectively, in each system. We find significant differences in the measured correlations of different centers, showing consistent off-centering effects in the three observables. By comparing the measurements with a halo model that treats the satellite photometric galaxies as being distributed according to a generalized NFW profile, we find that about 40 (70) per cent of BLRGs (FLRGs) are off-centered satellite galaxies in the multiple-LRG systems. The satellite LRGs have typical off-centering radius of about 400 kpc/h, and velocity dispersion of about 500 km/s in host halos with a mean mass of 1.6x10^14 Ms/h. We show that, if LRGs in the single LRG systems have similar offsets, the residual FoG contamination in the LRG power spectrum can be significant at k>0.1 h/Mpc, which may cause a bias in cosmological parameters such as the neutrino mass.

The effect of peculiar velocities on the epoch of reionization 21-cm signal

We have used semi-numerical simulations of reionization to study the behaviour of the power spectrum of the EoR 21-cm signal in redshift space. We have considered two models of reionization, one which has homogeneous recombination (HR) and the other incorporating inhomogeneous recombination (IR). We have estimated the observable quantities --- quadrupole and monopole moments of HI power spectrum at redshift space from our simulated data. We find that the magnitude and nature of the ratio between the quadrupole and monopole moments of the power spectrum ($P^s_2 /P^s_0$) can be a possible probe for the epoch of reionization. We observe that this ratio becomes negative at large scales for $x_{HI} \leq 0.7$ irrespective of the reionization model, which is a direct signature of an inside-out reionization at large scales. It is possible to qualitatively interpret the results of the simulations in terms of the fluctuations in the matter distribution and the fluctuations in the neutral fraction which have power spectra and cross-correlation $P_{\Delta \Delta}(k)$, $P_{xx}(k)$ and $P_{\Delta x}(k)$ respectively. We find that at large scales the fluctuations in matter density and neutral fraction is exactly anti-correlated through all stages of reionization. This provides a simple picture where we are able to qualitatively interpret the behaviour of the redshift space power spectra at large scales with varying $x_{HI}$ entirely in terms of a just two quantities, namely $x_{HI}$ and the ratio $P_{xx}/P_{\Delta \Delta}$. The nature of $P_{\Delta x}$ becomes different for HR and IR scenarios at intermediate and small scales. We further find that it is possible to distinguish between an inside-out and an outside-in reionization scenario from the nature of the ratio $P^s_2 /P^s_0$ at intermediate length scales.

Combining clustering and abundances of galaxy clusters to test cosmology and primordial non-Gaussianity

We present the clustering of galaxy clusters as a useful addition to the common set of cosmological observables. The clustering of clusters probes the large-scale structure of the Universe, extending galaxy clustering analysis to the high-peak, high-bias regime. Clustering of galaxy clusters complements the traditional cluster number counts and observable-mass relation analyses, significantly improving their constraining power by breaking existing calibration degeneracies. We use the maxBCG galaxy clusters catalogue to constrain cosmological parameters and cross-calibrate the mass-observable relation, using cluster abundances in richness bins and weak-lensing mass estimates. We then add the redshift-space power spectrum of the sample, including an effective modelling of the weakly non-linear contribution and allowing for an arbitrary photometric redshift smoothing. The inclusion of the power spectrum data allows for an improved self-calibration of the scaling relation. We find that the inclusion of the power spectrum typically brings a $\sim 50$ per cent improvement in the errors on the fluctuation amplitude $\sigma_8$ and the matter density $\Omega_{\mathrm{m}}$. Finally, we apply this method to constrain models of the early universe through the amount of primordial non-Gaussianity of the local type, using both the variation in the halo mass function and the variation in the cluster bias. We find a constraint on the amount of skewness $f_{\mathrm{NL}} = 12 \pm 157 $ ($1\sigma$) from the cluster data alone.

Understanding the nature of luminous red galaxies (LRGs): connecting LRGs to central and satellite subhaloes

We develop a novel abundance matching method to construct a mock catalog of luminous red galaxies (LRGs) in SDSS, using catalogs of halos and subhalos in N-body simulations for a LCDM model. Motivated by observations suggesting that LRGs are passively-evolving, massive early-type galaxies with a typical age >5Gyr, we assume that simulated halos at z=2 (z2-halo) are progenitors for LRG-host subhalos observed today, and we label the most tightly bound particles in each progenitor z2-halo as LRG ``stars''. We then identify the subhalos containing these stars to z=0.3 (SDSS redshift) in descending order of the masses of z2-halos until the comoving number density of the matched subhalos becomes comparable to the measured number density of SDSS LRGs, n=10^{-4} (h/Mpc)^3. Once the above prescription is determined, our only free parameter is the number density of halos identified at z=2 and this parameter is fixed to match the observed number density at z = 0.3. By tracing subsequent merging and assembly histories of each progenitor z2-halo, we can directly compute, from the mock catalog, the distributions of central and satellite LRGs and their internal motions in each host halo at z=0.3. While the SDSS LRGs are galaxies selected by the magnitude and color cuts from the SDSS images and are not necessarily a stellar-mass-selected sample, our mock catalog reproduces a host of SDSS measurements: the halo occupation distribution for central and satellite LRGs, the projected auto-correlation function of LRGs, the cross-correlation of LRGs with shapes of background galaxies (LRG-galaxy weak lensing), and the nonlinear redshift-space distortion effect, the Finger-of-God effect, in the angle-averaged redshift-space power spectrum.

Peculiar velocity decomposition, redshift space distortion, and velocity reconstruction in redshift surveys: The methodology

Massive spectroscopic surveys will measure the redshift space distortion (RSD) induced by galaxy peculiar velocity to unprecedented accuracy and open a new era of precision RSD cosmology. We develop a new method to improve the RSD modeling and to carry out robust reconstruction of the 3D large scale peculiar velocity through galaxy redshift surveys, in light of RSD. (1) We propose a mathematically unique and physically motivated decomposition of peculiar velocity into three eigencomponents: an irrotational component completely correlated with the underlying density field (${\mathbf{v}}_{\ensuremath{\delta}}$), an irrotational component uncorrelated with the density field (${\mathbf{v}}_{S}$), and a rotational (curl) component (${\mathbf{v}}_{B}$). The three components have different origins, different scale dependences, and different impacts on RSD. (2) This decomposition has the potential to simplify and improve the RSD modeling. (i) ${\mathbf{v}}_{B}$ damps the redshift space clustering. (ii) ${\mathbf{v}}_{S}$ causes both damping and enhancement to the redshift space power spectrum ${P}^{s}(k,u)$. Nevertheless, the leading order contribution to the enhancement has a ${u}^{4}$ directional dependence, distinctively different from the Kaiser formula. Here, $u\ensuremath{\equiv}{k}_{z}/k$, $k$ is the amplitude of the wave vector, and ${k}_{z}$ is the component along the line of sight. (iii) ${\mathbf{v}}_{\ensuremath{\delta}}$ is of the greatest importance for the RSD cosmology. We find that the induced redshift clustering shows a number of important deviations from the usual Kaiser formula. Even in the limit of ${\mathbf{v}}_{S}\ensuremath{\rightarrow}0$ and ${\mathbf{v}}_{B}\ensuremath{\rightarrow}0$, the leading order contribution $\ensuremath{\propto}(1+f\stackrel{\texttildelow{}}{W}(k){u}^{2}{)}^{2}$. It differs from the Kaiser formula by a window function $\stackrel{\texttildelow{}}{W}(k)$. Nonlinear evolution generically drives $\stackrel{\texttildelow{}}{W}(k)\ensuremath{\le}1$. We hence identify a significant systematical error causing underestimation of the structure growth parameter $f$ by as much as $O(10%)$ even at a relatively large scale $k=0.1h/\mathrm{Mpc}$. (iv) The velocity decomposition reveals the three origins of the ``finger-of-God'' (FOG) effect and suggests how to simplify and improve the modeling of FOG by treating the three components separately. (v) We derive a new formula for the redshift space power spectrum. Under the velocity decomposition scheme, all high order Gaussian corrections and non-Gaussian corrections of order ${\ensuremath{\delta}}^{3}$ can be taken into account without introducing extra model uncertainties. Here $\ensuremath{\delta}$ is the nonlinear overdensity. (3) The velocity decomposition clarifies issues in peculiar velocity reconstruction through 3D galaxy distribution. We discuss two possible ways to carry out the 3D ${\mathbf{v}}_{\ensuremath{\delta}}$ reconstruction. Both use the otherwise troublesome RSD in velocity reconstruction as a valuable source of information. Both have the advantage of rendering the reconstruction of a stochastic 3D field into the reconstruction of a deterministic window function ${W}^{s}(k,u)$ of limited degrees of freedom. Both can automatically and significantly alleviate the galaxy bias problem and, in the limit of a deterministic galaxy bias, completely overcome it. Part 1 of this series of works lays out the methodology. Companion papers Y. Zheng et al. (in preparation). will extensively evaluate its performance against N-body simulations.

Distribution function approach to redshift space distortions. Part III: halos and galaxies

It was recently shown that the power spectrum in redshift space can be written as a sum of cross-power spectra between number weighted velocity moments, of which the lowest are density and momentum density. We investigate numerically the properties of these power spectra for simulated galaxies and dark matter halos and compare them to the dark matter power spectra, generalizing the concept of the bias in density-density power spectra. Because all of the quantities are number weighted this approach is well defined even for sparse systems such as massive halos. This contrasts to the previous approaches to RSD where velocity correlations have been explored, but velocity field is a poorly defined concept for sparse systems. We find that the number density weighting leads to a strong scale dependence of the bias terms for momentum density auto-correlation and cross-correlation with density. This trend becomes more significant for the more biased halos and leads to an enhancement of RSD power relative to the linear theory. Fingers-of-god effects, which in this formalism come from the correlations of the higher order moments beyond the momentum density, lead to smoothing of the power spectrum and can reduce this enhancement of power from the scale dependent bias, but are relatively small for halos with no small scale velocity dispersion. In comparison, for a more realistic galaxy sample with satellites the small scale velocity dispersion generated by satellite motions inside the halos leads to a larger power suppression on small scales, but this depends on the satellite fraction and on the details of how the satellites are distributed inside the halo. We investigate several statistics such as the two-dimensional power spectrum P(k,μ), where μ is the angle between the Fourier mode and line of sight, its multipole moments, its powers of μ2, and configuration space statistics. Overall we find that the nonlinear effects in realistic galaxy samples such as luminous red galaxies affect the redshift space clustering on very large scales: for example, the quadrupole moment is affected by 10% for k < 0.1hMpc−1, which means that these effects need to be understood if we want to extract cosmological information from the redshift space distortions.

Distribution function approach to redshift space distortions. Part IV: perturbation theory applied to dark matter

We develop a perturbative approach to redshift space distortions (RSD) using the phase space distribution function approach and apply it to the dark matter redshift space power spectrum and its moments. RSD can be written as a sum over density weighted velocity moments correlators, with the lowest order being density, momentum density and stress energy density. We use standard and extended perturbation theory (PT) to determine their auto and cross correlators, comparing them to N-body simulations. We show which of the terms can be modeled well with the standard PT and whichneed additional terms that include higher order corrections which cannot be modeled in PT. Most of these additional terms are related to the small scale velocity dispersion effects, the so called finger of god (FoG) effects,which affect some, but not all, of the terms in this expansion, and which can be approximately modeled using a simple physically motivated ansatz such as the halo model. We point out that there are several velocity dispersions that enter into thedetailed RSD analysis with very different amplitudes, which can be approximately predicted by the halo model. In contrast to previous models our approach systematically includes all of the termsat a given order in PT and provides a physical interpretation for the small scale dispersion values. We investigate RSD power spectrum as a function of μ, the cosine of the angle between the Fourier mode and line of sight, focusing on the lowest order powers of μ and multipole momentswhich dominate the observable RSD power spectrum. Overall we find considerable success in modeling many, but not all, of the terms in this expansion. This is similar to the situation in real space, but predicting power spectrum in redshift spaceis more difficult because of the explicit influence of small scale dispersion type effects in RSD, which extend to very large scales.

Beyond the plane-parallel and Newtonian approach: wide-angle redshift distortions and convergence in general relativity

We extend previous analyses of wide-angle correlations in the galaxy power spectrum in redshift space to include all general relativistic effects. These general relativistic corrections to the standard approach become important on large scales and at high redshifts, and they lead to new terms in the wide-angle correlations. We show that in principle the new terms can produce corrections of nearly 10% on Gpc scales over the usual Newtonian approximation. General relativistic corrections will be important for future large-volume surveys such as SKA and Euclid, although the problem of cosmic variance will present a challenge in observing this.

Redshift-space distortions in f(R) gravity

We use large volume N-body simulations to predict the clustering of dark matter in redshift space in f(R) modified gravity cosmologies. This is the first time that the nonlinear matter and velocity fields have been resolved to such a high level of accuracy over a broad range of scales in this class of models. We find significant deviations from the clustering signal in standard gravity, with an enhanced boost in power on large scales and stronger damping on small scales in the f(R) models compared to GR at redshifts z<1. We measure the velocity divergence (P_\theta \theta) and matter (P_\delta \delta) power spectra and find a large deviation in the ratios \sqrt{P_\theta \theta/P_\delta \delta} and P_\delta \theta/P_\delta\delta, between the f(R) models and GR for 0.03<k/(h/Mpc)<0.5. In linear theory these ratios equal the growth rate of structure on large scales. Our results show that the simulated ratios agree with the growth rate for each cosmology (which is scale dependent in the case of modified gravity) only for extremely large scales, k<0.06h/Mpc at z=0. The velocity power spectrum is substantially different in the f(R) models compared to GR, suggesting that this observable is a sensitive probe of modified gravity. We demonstrate how to extract the matter and velocity power spectra from the 2D redshift space power spectrum, P(k,\mu), and can recover the nonlinear matter power spectrum to within a few percent for k<0.1h/Mpc. However, the model fails to describe the shape of the 2D power spectrum demonstrating that an improved model is necessary in order to reconstruct the velocity power spectrum accurately. The same model can match the monopole moment to within 3% for GR and 10% for the f(R) cosmology at k<0.2 h/Mpc at z=1. Our results suggest that the extraction of the velocity power spectrum from future galaxy surveys is a promising method to constrain deviations from GR.