No Arabic abstract
Perturbation theory (PT) has been used to interpret the observed nonlinear large-scale structure statistics at the quasi-linear regime. To facilitate the PT-based analysis, we have presented the GridSPT algorithm, a grid-based method to compute the nonlinear density and velocity fields in standard perturbation theory (SPT) from a given linear power spectrum. Here, we further put forward the approach by taking the redshift-space distortions into account. With the new implementation, we have, for the first time, generated the redshift-space density field to the fifth order and computed the next-to-next-to-leading order (2 loop) power spectrum and the next-to-leading order (1 loop) bispectrum of matter clustering in redshift space. By comparing the result with corresponding analytical SPT calculation and $N$-body simulations, we find that the SPT calculation (A) suffers much more from the UV sensitivity due to the higher-derivative operators and (B) deviates from the $N$-body results from the Fourier wavenumber smaller than real space $k_{rm max}$. Finally, we have shown that while Pade approximation removes spurious features in morphology, it does not improve the modeling of power spectrum and bispectrum.
We develop an analytical forward model based on perturbation theory to predict the redshift-space galaxy overdensity at the field level given a realization of the initial conditions. We find that the residual noise between the model and simulated galaxy density has a power spectrum that is white on large scales, with size comparable to the shot noise. In the mildly nonlinear regime, we see a $k^2mu^2$ correction to the noise power spectrum, corresponding to larger noise along the line of sight and on smaller scales. The parametric form of this correction has been predicted on theoretical grounds before, and our simulations provide important confirmation of its presence. We have also modeled the galaxy velocity at the field-level and compared it against simulated galaxy velocities, finding that about $10%$ of the galaxies are responsible for half of the rms velocity residual for our simulated galaxy sample.
We present the one-loop 2-point function of biased tracers in redshift space computed with Lagrangian perturbation theory, including a full resummation of both long-wavelength (infrared) displacements and associated velocities. The resulting model accurately predicts the power spectrum and correlation function of halos and mock galaxies from two different sets of N-body simulations at the percent level for quasi-linear scales, including the damping of the baryon acoustic oscillation signal due to the bulk motions of galaxies. We compare this full resummation with other, approximate, techniques including the moment expansion and Gaussian streaming model. We discuss infrared resummation in detail and compare our Lagrangian formulation with the Eulerian theory augmented by an infrared resummation based on splitting the input power spectrum into wiggle and no-wiggle components. We show that our model is able to recover unbiased cosmological parameters in mock data encompassing a volume much larger than what will be available to future galaxy surveys. We demonstrate how to efficiently compute the resulting expressions numerically, making available a fast Python code capable of rapidly computing these statistics in both configuration and Fourier space.
The large-scale matter distribution in the late-time Universe exhibits gravity-induced non-Gaussianity, and the bispectrum, three-point cumulant is expected to contain significant cosmological information. In particular, the measurement of the bispectrum helps to tighten the constraints on dark energy and modified gravity through the redshift-space distortions (RSD). In this paper, extending the work by Taruya, Nishimichi & Saito (2010, Phys.Rev.D 82, 063522), we present a perturbation theory (PT) based model of redshift-space matter bispectrum that can keep the non-perturbative damping effect under control. Characterizing this non-perturbative damping by a univariate function with single free parameter, the PT model of the redshift-space bispectrum is tested against a large set of cosmological $N$-body simulations, finding that the predicted monopole and quadrupole moments are in a good agreement with simulations at the scales of baryon acoustic oscillations (well beyond the range of agreement of standard PT). The validity of the univariate ansatz of the damping effect is also examined, and with the PT calculation at next-to-leading order, the fitted values of the free parameter is shown to consistently match those obtained from the PT model of power spectrum by Taruya, Nishimichi & Saito (2010).
Perturbation theory (PT) calculation of large-scale structure has been used to interpret the observed non-linear statistics of large-scale structure at the quasi-linear regime. In particular, the so-called standard perturbation theory (SPT) provides a basis for the analytical computation of the higher-order quantities of large-scale structure. Here, we present a novel, grid-based algorithm for the SPT calculation, hence named GridSPT, to generate the higher-order density and velocity fields from a given linear power spectrum. Taking advantage of the Fast Fourier Transform, the GridSPT quickly generates the nonlinear density fields at each order, from which we calculate the statistical quantities such as non-linear power spectrum and bispectrum. Comparing the density fields (to fifth order) from GridSPT with those from the full N-body simulations in the configuration space, we find that GridSPT accurately reproduces the N-body result on large scales. The agreement worsens with smaller smoothing radius, particularly for the under-dense regions where we find that 2LPT (second-order Lagrangian perturbation theory) algorithm performs well.
We present a self-consistent Bayesian formalism to sample the primordial density fields compatible with a set of dark matter density tracers after cosmic evolution observed in redshift space. Previous works on density reconstruction did not self-consistently consider redshift space distortions or included an additional iterative distortion correction step. We present here the analytic solution of coherent flows within a Hamiltonian Monte Carlo posterior sampling of the primordial density field. We test our method within the Zeldovich approximation, presenting also an analytic solution including tidal fields and spherical collapse on small scales using augmented Lagrangian perturbation theory. Our resulting reconstructed fields are isotropic and their power spectra are unbiased compared to the true one defined by our mock observations. Novel algorithmic implementations are introduced regarding the mass assignment kernels when defining the dark matter density field and optimization of the time step in the Hamiltonian equations of motions. Our algorithm, dubbed barcode, promises to be specially suited for analysis of the dark matter cosmic web down to scales of a few Megaparsecs. This large scale structure is implied by the observed spatial distribution of galaxy clusters --- such as obtained from X-ray, SZ or weak lensing surveys --- as well as that of the intergalactic medium sampled by the Lyman alpha forest or perhaps even by deep hydrogen intensity mapping. In these cases, virialized motions are negligible, and the tracers cannot be modeled as point-like objects. It could be used in all of these contexts as a baryon acoustic oscillation reconstruction algorithm.