No Arabic abstract
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.
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.
The observed power spectrum in redshift space appears distorted due to the peculiar motion of galaxies, known as redshift-space distortions (RSD). While all the effects in RSD are accounted for by the simple mapping formula from real to redshift spaces, accurately modeling redshift-space power spectrum is rather difficult due to the non-perturbative properties of the mapping. Still, however, a perturbative treatment may be applied to the power spectrum at large-scales, and on top of a careful modeling of the Finger-of-God effect caused by the small-scale random motion, the redshift-space power spectrum can be expressed as a series of expansion which contains the higher-order correlations of density and velocity fields. In our previous work [JCAP 8 (Aug., 2016) 050], we provide a perturbation-theory inspired model for power spectrum in which the higher-order correlations are evaluated directly from the cosmological $N$-body simulations. Adopting a simple Gaussian ansatz for Finger-of-God effect, the model is shown to quantitatively describe the simulation results. Here, we further push this approach, and present an accurate power spectrum template which can be used to estimate the growth of structure as a key to probe gravity on cosmological scales. Based on the simulations, we first calibrate the uncertainties and systematics in the pertrubation theory calculation in a fiducial cosmological model. Then, using the scaling relations, the calibrated power spectrum template is applied to a different cosmological model. We demonstrate that with our new template, the best-fitted growth functions are shown to reproduce the fiducial values in a good accuracy of 1 % at $k<0.18 hompc$ for cosmologies with different Hubble parameters.
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 explore the degrees of freedom required to jointly fit projected and redshift-space clustering of galaxies selected in three bins of stellar mass from the Sloan Digital Sky Survey Main Galaxy Sample (SDSS MGS) using a subhalo abundance matching (SHAM) model. We employ emulators for relevant clustering statistics in order to facilitate our analysis, leading to large speed gains with minimal loss of accuracy. We are able to simultaneously fit the projected and redshift-space clustering of the two most massive galaxy samples that we consider with just two free parameters: scatter in stellar mass at fixed SHAM proxy and the dependence of the SHAM proxy on dark matter halo concentration. We find some evidence for models that include velocity bias, but including orphan galaxies improves our fits to the lower mass samples significantly. We also model the clustering signals of specific star formation rate (SSFR) selected samples using conditional abundance matching (CAM). We obtain acceptable fits to projected and redshift-space clustering as a function of SSFR and stellar mass using two CAM variants, although the fits are worse than for stellar mass-selected samples alone. By incorporating non-unity correlations between the CAM proxy and SSFR we are able to resolve previously identified discrepancies between CAM predictions and SDSS observations of the environmental dependence of quenching for isolated central galaxies.
In the first of a series of forthcoming publications, we present a panchromatic catalog of 102 visually-selected early-type galaxies (ETGs) from observations in the Early Release Science (ERS) program with the Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST) of the Great Observatories Origins Deep Survey-South (GOODS-S) field. Our ETGs span a large redshift range, 0.35 < z < 1.5, with each redshift spectroscopically-confirmed by previous published surveys of the ERS field. We combine our measured WFC3 ERS and ACS GOODS-S photometry to gain continuous sensitivity from the rest-frame far-UV to near-IR emission for each ETG. The superior spatial resolution of the HST over this panchromatic baseline allows us to classify the ETGs by their small-scale internal structures, as well as their local environment. By fitting stellar population spectral templates to the broad-band photometry of the ETGs, we determine that the average masses of the ETGs are comparable to the characteristic stellar mass of massive galaxies, 11< log(M [Solar]) < 12. By transforming the observed photometry into the GALEX FUV and NUV, Johnson V, and SDSS g and r bandpasses we identify a noteworthy diversity in the rest-frame UV-optical colors and find the mean rest-frame (FUV-V)=3.5 and (NUV-V)=3.3, with 1$sigma$ standard deviations approximately equal to 1.0. The blue rest-frame UV-optical colors observed for most of the ETGs are evidence for star-formation during the preceding gigayear, but no systems exhibit UV-optical photometry consistent with major recent (<~50 Myr) starbursts. Future publications which address the diversity of stellar populations likely to be present in these ETGs, and the potential mechanisms by which recent star-formation episodes are activated, are discussed.