No Arabic abstract
We study the effect of the Einstein - de Sitter (EdS) approximation on the one-loop power spectrum of galaxies in redshift space in the Effective Field Theory of Large-Scale Structure. The dark matter density perturbations and velocity divergence are treated with exact time dependence. Splitting the density perturbation into its different temporal evolutions naturally gives rise to an irreducible basis of biases. While, as in the EdS approximation, at each time this basis spans a seven-dimensional space, this space is a slightly different one, and the difference is captured by a single calculable time- and $vec k$-dependent function. We then compute the redshift-space galaxy one-loop power spectrum with the EdS approximation ($P^{text{EdS-approx}}$) and without ($P^{text{Exact}}$). For the monopole we find $P_{text{0}}^{text{Exact}}/P_{text{0}}^{text{EdS-approx}}sim 1.003$ and for the quadrupole $P_{text{2}}^{text{Exact}}/P_{text{2}}^{text{EdS-approx}}sim 1.007$ at $z=0.57$, and sharply increasing at lower redshifts. Finally, we show that a substantial fraction of the effect remains even after allowing the bias coefficients to shift within a physically allowed range. This suggests that the EdS approximation can only fit the data to a level of precision that is roughly comparable to the precision of the next generation of cosmological surveys. Furthermore, we find that implementing the exact time dependence formalism is not demanding and is easily applicable to data. Both of these points motivate a direct study of this effect on the cosmological parameters.
We present the one-loop perturbation theory for the power spectrum of the marked density field of matter and biased tracers in real- and redshift-space. The statistic has been shown to yield impressive constraints on cosmological parameters; to exploit this, we require an accurate and computationally inexpensive theoretical model. Comparison with $N$-body simulations demonstrates that linear theory fails on all scales, but inclusion of one-loop Effective Field Theory terms gives a substantial improvement, with $sim 5%$ accuracy at $z = 1$. The expansion is less convergent in redshift-space (achieving $sim 10%$ accuracy), but there are significant improvements for biased tracers due to the freedom in the bias coefficients. The large-scale theory contains non-negligible contributions from all perturbative orders; we suggest a reorganization of the theory that contains all terms relevant on large-scales, discussing both its explicit form at one-loop and structure at infinite-loop. This motivates a low-$k$ correction term, leading to a model that is sub-percent accurate on large scales, albeit with the inclusion of two (three) free coefficients in real- (redshift-)space. We further consider the effects of massive neutrinos, showing that beyond-EdS corrections to the perturbative kernels are negligible in practice. It remains to see whether the purported gains in cosmological parameters remain valid for biased tracers and can be captured by the theoretical model.
We present an emulator for the two-point clustering of biased tracers in real space. We construct this emulator using neural networks calibrated with more than $400$ cosmological models in a 8-dimensional cosmological parameter space that includes massive neutrinos an dynamical dark energy. The properties of biased tracers are described via a Lagrangian perturbative bias expansion which is advected to Eulerian space using the displacement field of numerical simulations. The cosmology-dependence is captured thanks to a cosmology-rescaling algorithm. We show that our emulator is capable of describing the power spectrum of galaxy formation simulations for a sample mimicking that of a typical Emission-Line survey at $z sim 1$ with an accuracy of $1-2%$ up to nonlinear scales $k sim 0.7 h mathrm{Mpc}^{-1}$.
We extend the scale-dependent Gaussian Streaming Model (GSM) to produce analytical predictions for the anisotropic redshift-space correlation function for biased tracers in modified gravity models. Employing the Convolution Lagrangian Perturbation Theory (CLPT) re-summation scheme, with a local Lagrangian bias schema provided by the peak-background split formalism, we predict the necessary ingredients that enter the GSM, the real-space halo pairwise velocity and the pairwise velocity dispersion. We further consider effective field theory contributions to the pairwise velocity dispersion in order to model correctly its large scale behavior. We apply our method on two widely-considered modified gravity models, the chameleon-screened f(R) Hu-Sawicki model and the nDGP Vainshtein model and compare our predictions against state-of-the-art N-body simulations for these models. We demonstrate that the GSM approach to predict the monopole and the quadrupole of the redshift-space correlation function for halos, gives very good agreement with the simulation data, for a wide range of screening mechanisms, levels of screening and halo masses at z=0.5 and z=1. Our work shows the applicability of the GSM, for cosmologies beyond GR, demonstrating that it can serve as a powerful predictive tool for the next stage of cosmological surveys like DESI, Euclid, LSST and WFIRST.
In this paper we test the perturbative halo bias model at the field level. The advantage of this approach is that any analysis can be done without sample variance if the same initial conditions are used in simulations and perturbation theory calculations. We write the bias expansion in terms of modified bias operators in Eulerian space, designed such that the large bulk flows are automatically resummed and not treated perturbatively. Using these operators, the bias model accurately matches the Eulerian density of halos in N-body simulations. The mean-square model error is close to the Poisson shot noise for a wide range of halo masses and it is rather scale-independent, with scale-dependent corrections becoming relevant at the nonlinear scale. In contrast, for linear bias the mean-square model error can be higher than the Poisson prediction by factors of up to a few on large scales, and it becomes scale dependent already in the linear regime. We show that by weighting simulated halos by their mass, the mean-square error of the model can be further reduced by up to an order of magnitude, or by a factor of two when including $60%$ mass scatter. We also test the Standard Eulerian bias model using the nonlinear matter field measured from simulations and show that it leads to a larger and more scale-dependent model error than the bias expansion based on perturbation theory. These results may be of particular relevance for cosmological inference methods that use a likelihood of the biased tracer at the field level, or for initial condition and BAO reconstruction that requires a precise estimate of the large-scale potential from the biased tracer density.
Upcoming galaxy redshift surveys promise to significantly improve current limits on primordial non-Gaussianity (PNG) through measurements of 2- and 3-point correlation functions in Fourier space. However, realizing the full potential of this dataset is contingent upon having both accurate theoretical models and optimized analysis methods. Focusing on the local model of PNG, parameterized by $f_{rm NL}$, we perform a Monte-Carlo Markov Chain analysis to confront perturbation theory predictions of the halo power spectrum and bispectrum in real space against a suite of N-body simulations. We model the halo bispectrum at tree-level, including all contributions linear and quadratic in $f_{rm NL}$, and the halo power spectrum at 1-loop, including tree-level terms up to quadratic order in $f_{rm NL}$ and all loops induced by local PNG linear in $f_{rm NL}$. Keeping the cosmological parameters fixed, we examine the effect of informative priors on the linear non-Gaussian bias parameter on the statistical inference of $f_{rm NL}$. A conservative analysisof the combined power spectrum and bispectrum, in which only loose priors are imposed and all parameters are marginalized over, can improve the constraint on $f_{rm NL}$ by more than a factor of 5 relative to the power spectrum-only measurement. Imposing a strong prior on $b_phi$, or assuming bias relations for both $b_phi$ and $b_{phidelta}$ (motivated by a universal mass function assumption), improves the constraints further by a factor of few. In this case, however, we find a significant systematic shift in the inferred value of $f_{rm NL}$ if the same range of wavenumber is used. Likewise, a Poisson noise assumption can lead to significant systematics, and it is thus essential to leave all the stochastic amplitudes free.