Source author record

Zvonimir Vlah

Zvonimir Vlah appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

22works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2022arXiv

A new analysis of galaxy 2-point functions in the BOSS survey, including full-shape information and post-reconstruction BAO

We present a new method for consistent, joint analysis of the pre- and post-reconstruction two-point functions of the BOSS survey. The post-reconstruction correlation function is used to accurately measure the distance-redshift relation and expansion history, while the pre-reconstruction power spectrum multipoles constrain the broad-band shape and the rate-of-growth of large-scale structure. Our technique uses Lagrangian perturbation theory to self-consistently work at the level of two-point functions, i.e.\ directly with the measured data, without approximating the constraints with summary statistics normalized by the drag scale. Combining galaxies across the full redshift range and both hemispheres we constrain $Ω_m=0.303 \pm 0.0082$, $H_0=69.23 \pm 0.77$ and $σ_8=0.733 \pm 0.047$ within the context of $Λ$CDM. These constraints are in good agreement both with the Planck primary CMB anisotropy data and recent cosmic shear surveys.

preprint2020arXiv

Consistent Modeling of Velocity Statistics and Redshift-Space Distortions in One-Loop Perturbation Theory

The peculiar velocities of biased tracers of the cosmic density field contain important information about the growth of large scale structure and generate anisotropy in the observed clustering of galaxies. Using N-body data, we show that velocity expansions for halo redshift-space power spectra are converged at the percent-level at perturbative scales for most line-of-sight angles $μ$ when the first three pairwise velocity moments are included, and that the third moment is well-approximated by a counterterm-like contribution. We compute these pairwise-velocity statistics in Fourier space using both Eulerian and Lagrangian one-loop perturbation theory using a cubic bias scheme and a complete set of counterterms and stochastic contributions. We compare the models and show that our models fit both real-space velocity statistics and redshift-space power spectra for both halos and a mock sample of galaxies at sub-percent level on perturbative scales using consistent sets of parameters, making them appealing choices for the upcoming era of spectroscopic, peculiar-velocity and kSZ surveys.

preprint2020arXiv

The Reconstructed Power Spectrum in the Zeldovich Approximation

Density-reconstruction sharpens the baryon acoustic oscillations signal by undoing some of the smoothing incurred by nonlinear structure formation. In this paper we present an analytical model for reconstruction based on the Zeldovich approximation, which for the first time includes a complete set of counterterms and bias terms up to quadratic order and can fit real and redshift-space data pre- and post-reconstruction data in both Fourier and configuration space over a wide range of scales. We compare our model to n-body data at $z = 0$ from the {\tt DarkSky} simulation \cite{skillman14}, finding sub-percent agreement in both real space and in the redshift-space power spectrum monopole out to $k = 0.4\ h$ Mpc$^{-1}$, and out to $k = 0.2\ h$ Mpc$^{-1}$ in the quadrupole, with comparable agreement in configuration space. We compare our model with several popular existing alternatives, updating existing theoretical results for exponential damping in wiggle/no-wiggle splits of the BAO signal and discuss the usually-ignored effect of higher bias contributions on the reconstructed signal. In the appendices, we re-derive the former within our formalism, present exploratory results on higher-order corrections due to nonlinearities inherent to reconstruction, and present numerical techniques with which to calculate the redshift-space power spectrum of biased tracers within the Zeldovich approximation.

preprint2019arXiv

An EFT description of galaxy intrinsic alignments

We present a general perturbative effective field theory (EFT) description of galaxy shape correlations, which are commonly known as intrinsic alignments. This rigorous approach extends current analytical modelling strategies in that it only relies on the equivalence principle. We present our results in terms of three-dimensional statistics for two- and three-point functions of both galaxy shapes and number counts. In case of the two-point function, we recover the well-known linear alignment result at leading order, but also present the full next-to-leading order expressions. In case of the three-point function we present leading order results for all the auto- and cross-correlations of galaxy shapes and densities. We use a spherical tensor basis to decompose the tensor perturbations in different helicity modes, which allows us to make use of isotropy and parity properties in the correlators. Combined with the results on projection presented in a forthcoming companion paper, our framework is directly applicable to accounting for intrinsic alignment contamination in weak lensing surveys, and to extracting cosmological information from intrinsic alignments.

preprint2016arXiv

Fast Large Scale Structure Perturbation Theory using 1D FFTs

The usual fluid equations describing the large-scale evolution of mass density in the universe can be written as local in the density, velocity divergence, and velocity potential fields. As a result, the perturbative expansion in small density fluctuations, usually written in terms of convolutions in Fourier space, can be written as a series of products of these fields evaluated at the same location in configuration space. Based on this, we establish a new method to numerically evaluate the 1-loop power spectrum (i.e., Fourier transform of the 2-point correlation function) with one-dimensional Fast Fourier Transforms. This is exact and a few orders of magnitude faster than previously used numerical approaches. Numerical results of the new method are in excellent agreement with the standard quadrature integration method. This fast model evaluation can in principle be extended to higher loop order where existing codes become painfully slow. Our approach follows by writing higher order corrections to the 2-point correlation function as, e.g., the correlation between two second-order fields or the correlation between a linear and a third-order field. These are then decomposed into products of correlations of linear fields and derivatives of linear fields. The method can also be viewed as evaluating three-dimensional Fourier space convolutions using products in configuration space, which may also be useful in other contexts where similar integrals appear.

preprint2016arXiv

FFT-PT: Reducing the two-loop large-scale structure power spectrum to low-dimensional radial integrals

Modeling the large-scale structure of the universe on nonlinear scales has the potential to substantially increase the science return of upcoming surveys by increasing the number of modes available for model comparisons. One way to achieve this is to model nonlinear scales perturbatively. Unfortunately, this involves high-dimensional loop integrals that are cumbersome to evaluate. Trying to simplify this, we show how two-loop (next-to-next-to-leading order) corrections to the density power spectrum can be reduced to low-dimensional, radial integrals. Many of those can be evaluated with a one-dimensional Fast Fourier Transform, which is significantly faster than the five-dimensional Monte-Carlo integrals that are needed otherwise. The general idea of this FFT-PT method is to switch between Fourier and position space to avoid convolutions and integrate over orientations, leaving only radial integrals. This reformulation is independent of the underlying shape of the initial linear density power spectrum and should easily accommodate features such as those from baryonic acoustic oscillations. We also discuss how to account for halo bias and redshift space distortions.

preprint2016arXiv

Massive neutrinos in nonlinear large scale structure: A consistent perturbation theory

A consistent formulation to incorporate massive neutrinos in the perturbation theory of the effective CDM+baryons fluid is introduced. In this formulation all linear k dependence in the growth functions of CDM+baryons perturbations, as well as all consequent additional mode coupling at higher orders, are taken into account to any desirable accuracy. Our formulation regards the neutrino fraction, which is constant in time after the non-relativistic transition of neutrinos, and much smaller than unity, as the coupling constant of the theory. Then the "bare" perturbations are those in the massless neutrino case when the neutrino fraction vanishes, and we consider the backreaction corrections due to the gravitational coupling of neutrinos. We derive the general equations for the "bare" perturbations, and backrecation corrections. Then, by employing exact time evolution with the proper analytic Green's function we explicitly derive the leading backreaction effect, and find precise agreement at the linear level. We proceed to compute the second order beackreaction correction, and derive the leading order matter bispectrum in the presence of massive neutrinos, suggesting the squeezed limit of the matter bispectrum as a sensitive probe of neutrino mass. Notably, the generic neutrino fraction formulation in this work may be similarly applied for the consistent inclusion of massive neutrinos within any perturbative approach.

preprint2016arXiv

Non-linear Fields in Generalized Cosmologies

The perturbative approach to structure formation has recently received a lot of attention in the literature. In such setups the final predictions for observables like the power spectrum is often derived under additional approximations such as a simplified time dependence. Here we provide all-order perturbative integral solutions for density and velocity fields in generalized cosmologies, with a direct application to clustering quintessence. We go beyond the standard results based on extending the EdS-like approximations. As an illustrative example, we apply our findings to the calculation of the one-loop power spectrum of density and momentum fields. We find corrections close to $1\%$ in the mildly non-linear regime of $Λ$CDM cosmologies for the density power spectrum, while in the case of the density-momentum power spectrum effects can reach up to $1.5\%$ for $k\sim 0.2h/$Mpc.

preprint2016arXiv

On the Statistics of Biased Tracers in the Effective Field Theory of Large Scale Structures

With the completion of the Planck mission, in order to continue to gather cosmological information it has become crucial to understand the Large Scale Structures (LSS) of the universe to percent accuracy. The Effective Field Theory of LSS (EFTofLSS) is a novel theoretical framework that aims to develop an analytic understanding of LSS at long distances, where inhomogeneities are small. We further develop the description of biased tracers in the EFTofLSS to account for the effect of baryonic physics and primordial non-Gaussianities, finding that new bias coefficients are required. Then, restricting to dark matter with Gaussian initial conditions, we describe the prediction of the EFTofLSS for the one-loop halo-halo and halo-matter two-point functions, and for the tree-level halo-halo-halo, matter-halo-halo and matter-matter-halo three-point functions. Several new bias coefficients are needed in the EFTofLSS, even though their contribution at a given order can be degenerate and the same parameters contribute to multiple observables. We develop a method to reduce the number of biases to an irreducible basis, and find that, at the order at which we work, seven bias parameters are enough to describe this extremely rich set of statistics. We then compare with the output of $N$-body simulations. For the lowest mass bin, we find percent level agreement up to $k\simeq 0.3\,h\,{\rm Mpc}^{-1}$ for the one-loop two-point functions, and up to $k\simeq 0.15\,h\,{\rm Mpc}^{-1}$ for the tree-level three-point functions, with the $k$-reach decreasing with higher mass bins. This is consistent with the theoretical estimates, and suggests that the cosmological information in LSS amenable to analytical control is much more than previously believed.

preprint2016arXiv

Perturbation theory, effective field theory, and oscillations in the power spectrum

We explore the relationship between the nonlinear matter power spectrum and the various Lagrangian and Standard Perturbation Theories (LPT and SPT). We first look at it in the context of one dimensional (1-d) dynamics, where 1LPT is exact at the perturbative level and one can exactly resum the SPT series into the 1LPT power spectrum. Shell crossings lead to non-perturbative effects, and the PT ignorance can be quantified in terms of their ratio, which is also the transfer function squared in the absence of stochasticity. At the order of PT we work, this parametrization is equivalent to the results of effective field theory (EFT), and can thus be expanded in terms of the same parameters. We find that its radius of convergence is larger than the SPT loop expansion. The same EFT parametrization applies to all SPT loop terms and, if stochasticity can be ignored, to all N-point correlators. In 3-d, the LPT structure is considerably more complicated, and we find that LPT models with parametrization motivated by the EFT exhibit running with $k$ and that SPT is generally a better choice. Since these transfer function expansions contain free parameters that change with cosmological model their usefulness for broadband power is unclear. For this reason we test the predictions of these models on baryonic acoustic oscillations (BAO) and other primordial oscillations, including string monodromy models, for which we ran a series of simulations with and without oscillations. Most models are successful in predicting oscillations beyond their corresponding PT versions, confirming the basic validity of the model.

preprint2016arXiv

The Gaussian streaming model and Lagrangian effective field theory

We update the ingredients of the Gaussian streaming model (GSM) for the redshift-space clustering of biased tracers using the techniques of Lagrangian perturbation theory, effective field theory (EFT) and a generalized Lagrangian bias expansion. After relating the GSM to the cumulant expansion, we present new results for the real-space correlation function, mean pairwise velocity and pairwise velocity dispersion including counter terms from EFT and bias terms through third order in the linear density, its leading derivatives and its shear up to second order. We discuss the connection to the Gaussian peaks formalism. We compare the ingredients of the GSM to a suite of large N-body simulations, and show the performance of the theory on the low order multipoles of the redshift-space correlation function and power spectrum. We highlight the importance of a general biasing scheme, which we find to be as important as higher-order corrections due to non-linear evolution for the halos we consider on the scales of interest to us.

preprint2016arXiv

Very Massive Tracers and Higher Derivative Biases

Most of the upcoming cosmological information will come from analyzing the clustering of the Large Scale Structures (LSS) of the universe through LSS or CMB observations. It is therefore essential to be able to understand their behavior with exquisite precision. The Effective Field Theory of Large Scale Structures (EFTofLSS) provides a consistent framework to make predictions for LSS observables in the mildly non-linear regime. In this paper we focus on biased tracers. We argue that in calculations at a given order in the dark matter perturbations, highly biased tracers will underperform because of their larger higher derivative biases. A natural prediction of the EFTofLSS is therefore that by simply adding higher derivative biases, all tracers should perform comparably well. We implement this prediction for the halo-halo and the halo-matter power spectra at one loop, and the halo-halo-halo, halo-halo-matter, and halo-matter-matter bispectra at tree-level, and compare with simulations. We find good agreement with the prediction: for all tracers, we are able to match the bispectra up to $k\simeq0.17\,h/$Mpc at $z=0$ and the power spectra to a higher wavenumber.

preprint2015arXiv

A Lagrangian effective field theory

We have continued the development of Lagrangian, cosmological perturbation theory for the low-order correlators of the matter density field. We provide a new route to understanding how the effective field theory (EFT) of large-scale structure can be formulated in the Lagrandian framework and a new resummation scheme, comparing our results to earlier work and to a series of high-resolution N-body simulations in both Fourier and configuration space. The `new' terms arising from EFT serve to tame the dependence of perturbation theory on small-scale physics and improve agreement with simulations (though with an additional free parameter). We find that all of our models fare well on scales larger than about two to three times the non-linear scale, but fail as the non-linear scale is approached. This is slightly less reach than has been seen previously. At low redshift the Lagrangian model fares as well as EFT in its Eulerian formulation, but at higher $z$ the Eulerian EFT fits the data to smaller scales than resummed, Lagrangian EFT. All the perturbative models fare better than linear theory.

preprint2015arXiv

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.

preprint2015arXiv

Halo Zeldovich model and perturbation theory: dark matter power spectrum and correlation function

Perturbation theory for dark matter clustering has received a lot of attention in recent years, but its convergence properties remain poorly justified and there is no successful model that works both for correlation function and for power spectrum. Here we present Halo Zeldovich approach combined with perturbation theory (HZPT), in which we use standard perturbation theory at one loop order (SPT) at very low $k$, and connect it to a version of the halo model, for which we adopt the Zeldovich approximation plus a Pade expansion of a compensated one halo term. This low $k$ matching allows us to determine the one halo term amplitude and redshift evolution, both of which are in an excellent agreement with simulations, and approximately agree with the expected value from the halo model. Our Pade expansion approach of the one halo term added to Zeldovich approximation identifies two typical halo scales averaged over the halo mass function, the halo radius scale of order of 1Mpc/h, and the halo mass compensation scale of order 26Mpc/h. The model gives better than one percent accurate predictions for the correlation function above 5Mpc/h at all redshifts, without any free parameters. With three fitted Pade expansion coefficients the agreement in power spectrum is good to a percent up to $k \sim 1$h/Mpc, which can be improved to arbitrary $k$ by adding higher order terms in Pade expansion.

preprint2015arXiv

Tidal alignment of galaxies

We develop an analytic model for galaxy intrinsic alignments (IA) based on the theory of tidal alignment. We calculate all relevant nonlinear corrections at one-loop order, including effects from nonlinear density evolution, galaxy biasing, and source density weighting. Contributions from density weighting are found to be particularly important and lead to bias dependence of the IA amplitude, even on large scales. This effect may be responsible for much of the luminosity dependence in IA observations. The increase in IA amplitude for more highly biased galaxies reflects their locations in regions with large tidal fields. We also consider the impact of smoothing the tidal field on halo scales. We compare the performance of this consistent nonlinear model in describing the observed alignment of luminous red galaxies with the linear model as well as the frequently used "nonlinear alignment model," finding a significant improvement on small and intermediate scales. We also show that the cross-correlation between density and IA (the "GI" term) can be effectively separated into source alignment and source clustering, and we accurately model the observed alignment down to the one-halo regime using the tidal field from the fully nonlinear halo-matter cross correlation. Inside the one-halo regime, the average alignment of galaxies with density tracers no longer follows the tidal alignment prediction, likely reflecting nonlinear processes that must be considered when modeling IA on these scales. Finally, we discuss tidal alignment in the context of cosmic shear measurements.

preprint2014arXiv

Geometric and dynamic distortions in anisotropic galaxy clustering

We examine the signature of dynamic (redshift-space) distortions and geometric distortions (including the Alcock-Paczynski effect) in the context of the galaxy power spectrum measured in upcoming galaxy redshift surveys. Information comes from both the baryon acoustic oscillation (BAO) feature and the broadband power spectrum shape. Accurate modeling is required to extract this information without introducing systematic biases in the result. We consider an analytic model for the power spectrum of dark matter halos in redshift space, based on the distribution function expansion, and compare with halo clustering measured in N-body simulations. We forecast that the distribution function model is sufficiently accurate to allow the inclusion of broadband information on scales down to k~0.2 h/Mpc, with somewhat better accuracy for higher bias halos. Compared with a BAO-only analysis with reconstruction, including broadband shape information can improve unbiased constraints on distance measures H(z) and D_A(z) by ~30% and 20%, respectively, for a galaxy sample similar to the DESI luminous red galaxies. The gains in precision are larger in the absence of BAO reconstruction. Furthermore, including broadband shape information allows the measurement of structure growth, through redshift-space distortions. For the same galaxy sample, the distribution function model is able to constrain f*sigma_8 to ~2%, when simultaneously fitting for H(z) and D_A(z). We discuss techniques to optimize the analysis of the power spectrum, including removing modes near the line-of-sight that are particularly challenging to model, and whether these approaches can improve parameter constraints. We find that such techniques are unlikely to significantly improve constraints on geometry, although they may allow higher precision measurements of redshift-space distortions.

preprint2014arXiv

Lagrangian perturbation theory at one loop order: successes, failures, and improvements

We apply the convolved Lagrangian perturbation theory (CLPT) formalism, in which one can express the matter density power spectrum in terms of integrals over a function of cumulants of the displacement field, allowing for a resummation of the terms, to evaluate the full one loop power spectrum. We keep the cumulants up to third order, extending the Zel'dovich approximation and providing the power spectrum analogous to the calculations recently performed for the correlation function. We compare the results to the N-body simulations and to the Lagrangian perturbation simulations up to the second order. We find that the analytic calculations are in a good agreement with the LPT simulations, but when compared to full N-body simulations, we find that while one loop calculations improve upon the Zel'dovich approximation in the power spectrum, they still significantly lack power. As found previously in the correlation function one loop CLPT improves slightly against Zel'dovich above 30Mpc/h, but is actually worse than Zel'dovich below that. We investigate the deficiencies of the CLPT approach and argue that main problem of CLPT is its inability to trap particles inside dark matter halos, which leads to an overestimate of the small scale power of the displacement field and to an underestimate of the small scale power from one halo term effects. We model this using the displacement field damped at a nonlinear scale (CLPTs). To explore this in more detail we decompose the power spectrum and correlation function into three additive components: Zel'dovich, residual BAO wiggle, and residual broad band. One loop CLPT predicts small modifications to BAO wiggles that are enhanced in CLPTs, with up to 5\% corrections to correlation function around BAO scale.

preprint2014arXiv

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.

preprint2014arXiv

Understanding higher-order nonlocal halo bias at large scales by combining the power spectrum with the bispectrum

Understanding the relation between underlying matter distribution and biased tracers such as galaxy or dark matter halo is essential to extract cosmological information from ongoing or future galaxy redshift surveys. At sufficiently large scales such as the BAO scale, a standard approach for the bias problem on the basis of the perturbation theory (PT) is to assume the `local bias' model in which the density field of biased tracers is deterministically expanded in terms of matter density field at the same position. The higher-order bias parameters are then determined by combining the power spectrum with higher-order statistics such as the bispectrum. As is pointed out by recent studies, however, nonlinear gravitational evolution naturally induces nonlocal bias terms even if initially starting only with purely local bias. As a matter of fact, previous works showed that the second-order nonlocal bias term, which corresponds to the gravitational tidal field, is important to explain the characteristic scale-dependence of the bispectrum. In this paper we extend the nonlocal bias term up to third order, and investigate whether the PT-based model including nonlocal bias terms can simultaneously explain the power spectrum and the bispectrum of simulated halos in $N$-body simulations. We show that the power spectrum, including density and momentum, and the bispectrum between halo and matter in $N$-body simulations can be simultaneously well explained by the model including up to third-order nonlocal bias terms up to k~0.1h/Mpc. Also, the results seem in a good agreement with theoretical predictions of a simple coevolution picture, although the agreement is not perfect. These demonstration clearly shows a failure of the local bias model even at such large scales, and we conclude that nonlocal bias terms should be consistently included in order to model statistics of halos. [abridged]

preprint2013arXiv

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

Numerical simulations show that redshift space distortions (RSD) introduce strong scale dependence in the power spectra of halos, with ten percent deviations relative to linear theory predictions even on relatively large scales (k<0.1h/Mpc) and even in the absence of satellites (which induce Fingers-of-God, FoG, effects). If unmodeled these effects prevent one from extracting cosmological information from RSD surveys. In this paper we use perturbation theory (PT) and halo biasing model and apply it to the distribution function approach to RSD, in which RSD is decomposed into several correlators of density weighted velocity moments. We model each of these correlators using PT and compare the results to simulations over a wide range of halo masses and redshifts. We find that with an introduction of a physically motivated halo biasing, and using dark matter power spectra from simulations, we can reproduce the simulation results at a percent level on scales up to k~0.15h/Mpc at z=0, without the need to have free FoG parameters in the model.

preprint2012arXiv

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 which need 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 the detailed 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 terms at 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 moments which dominate the observable RSD power spectrum. Overall we find considerable success in modeling many, but not all, of the terms in this expansion.