skip to main content

Title: The cosmology dependence of galaxy clustering and lensing from a hybrid N -body–perturbation theory model
ABSTRACT We implement a model for the two-point statistics of biased tracers that combines dark matter dynamics from N-body simulations with an analytic Lagrangian bias expansion. Using Aemulus, a suite of N-body simulations built for emulation of cosmological observables, we emulate the cosmology dependence of these non-linear spectra from redshifts z = 0 to z = 2. We quantify the accuracy of our emulation procedure, which is sub-per cent at $k=1\, h \,{\rm Mpc}^{-1}$ for the redshifts probed by upcoming surveys and improves at higher redshifts. We demonstrate its ability to describe the statistics of complex tracer samples, including those with assembly bias and baryonic effects, reliably fitting the clustering and lensing statistics of such samples at redshift z ≃ 0.4 to scales of $k_{\rm max} \approx 0.6\, h\,\mathrm{Mpc}^{-1}$. We show that the emulator can be used for unbiased cosmological parameter inference in simulated joint clustering and galaxy–galaxy lensing analyses with data drawn from an independent N-body simulation. These results indicate that our emulator is a promising tool that can be readily applied to the analysis of current and upcoming data sets from galaxy surveys.
; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Sponsoring Org:
National Science Foundation
More Like this

    The combination of galaxy–galaxy lensing (GGL) and galaxy clustering is a powerful probe of low-redshift matter clustering, especially if it is extended to the non-linear regime. To this end, we use an N-body and halo occupation distribution (HOD) emulator method to model the redMaGiC sample of colour-selected passive galaxies in the Dark Energy Survey (DES), adding parameters that describe central galaxy incompleteness, galaxy assembly bias, and a scale-independent multiplicative lensing bias Alens. We use this emulator to forecast cosmological constraints attainable from the GGL surface density profile ΔΣ(rp) and the projected galaxy correlation function wp, gg(rp) in the final (Year 6) DES data set over scales $r_p=0.3\!-\!30.0\, h^{-1} \, \mathrm{Mpc}$. For a $3{{\ \rm per\ cent}}$ prior on Alens we forecast precisions of $1.9{{\ \rm per\ cent}}$, $2.0{{\ \rm per\ cent}}$, and $1.9{{\ \rm per\ cent}}$ on Ωm, σ8, and $S_8 \equiv \sigma _8\Omega _m^{0.5}$, marginalized over all halo occupation distribution (HOD) parameters as well as Alens. Adding scales $r_p=0.3\!-\!3.0\, h^{-1} \, \mathrm{Mpc}$ improves the S8 precision by a factor of ∼1.6 relative to a large scale ($3.0\!-\!30.0\, h^{-1} \, \mathrm{Mpc}$) analysis, equivalent to increasing the survey area by a factor of ∼2.6. Sharpening the Alens prior to $1{{\more »\rm per\ cent}}$ further improves the S8 precision to $1.1{{\ \rm per\ cent}}$, and it amplifies the gain from including non-linear scales. Our emulator achieves per cent-level accuracy similar to the projected DES statistical uncertainties, demonstrating the feasibility of a fully non-linear analysis. Obtaining precise parameter constraints from multiple galaxy types and from measurements that span linear and non-linear clustering offers many opportunities for internal cross-checks, which can diagnose systematics and demonstrate the robustness of cosmological results.

    « less

    We describe our non-linear emulation (i.e. interpolation) framework that combines the halo occupation distribution (HOD) galaxy bias model with N-body simulations of non-linear structure formation, designed to accurately predict the projected clustering and galaxy–galaxy lensing signals from luminous red galaxies in the redshift range 0.16 < z < 0.36 on comoving scales 0.6 < rp < 30 $h^{-1} \, \text{Mpc}$. The interpolation accuracy is ≲ 1–2 per cent across the entire physically plausible range of parameters for all scales considered. We correctly recover the true value of the cosmological parameter S8 = (σ8/0.8228)(Ωm/0.3107)0.6 from mock measurements produced via subhalo abundance matching (SHAM)-based light-cones designed to approximately match the properties of the SDSS LOWZ galaxy sample. Applying our model to Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 14 (DR14) LOWZ galaxy clustering and galaxy-shear cross-correlation measurements made with Sloan Digital Sky Survey (SDSS) Data Release 8 (DR8) imaging, we perform a prototype cosmological analysis marginalizing over wCDM cosmological parameters and galaxy HOD parameters. We obtain a 4.4 per cent measurement of S8 = 0.847 ± 0.037, in 3.5σ tension with the Planck cosmological results of 1.00 ± 0.02. We discuss the possibility of underestimated systematic uncertainties or astrophysical effects that could explain this discrepancy.


    We measure the small-scale clustering of the Data Release 16 extended Baryon Oscillation Spectroscopic Survey Luminous Red Galaxy sample, corrected for fibre-collisions using Pairwise Inverse Probability weights, which give unbiased clustering measurements on all scales. We fit to the monopole and quadrupole moments and to the projected correlation function over the separation range $7-60\, h^{-1}{\rm Mpc}$ with a model based on the aemulus cosmological emulator to measure the growth rate of cosmic structure, parametrized by fσ8. We obtain a measurement of fσ8(z = 0.737) = 0.408 ± 0.038, which is 1.4σ lower than the value expected from 2018 Planck data for a flat ΛCDM model, and is more consistent with recent weak-lensing measurements. The level of precision achieved is 1.7 times better than more standard measurements made using only the large-scale modes of the same sample. We also fit to the data using the full range of scales $0.1\text{--}60\, h^{-1}{\rm Mpc}$ modelled by the aemulus cosmological emulator and find a 4.5σ tension in the amplitude of the halo velocity field with the Planck + ΛCDM model, driven by a mismatch on the non-linear scales. This may not be cosmological in origin, and could be due to a breakdown in the Halo Occupation Distribution model used inmore »the emulator. Finally, we perform a robust analysis of possible sources of systematics, including the effects of redshift uncertainty and incompleteness due to target selection that were not included in previous analyses fitting to clustering measurements on small scales.

    « less
  4. ABSTRACT Cosmological constraints from current and upcoming galaxy cluster surveys are limited by the accuracy of cluster mass calibration. In particular, optically identified galaxy clusters are prone to selection effects that can bias the weak lensing mass calibration. We investigate the selection bias of the stacked cluster lensing signal associated with optically selected clusters, using clusters identified by the redMaPPer algorithm in the Buzzard simulations as a case study. We find that at a given cluster halo mass, the residuals of redMaPPer richness and weak lensing signal are positively correlated. As a result, for a given richness selection, the stacked lensing signal is biased high compared with what we would expect from the underlying halo mass probability distribution. The cluster lensing selection bias can thus lead to overestimated mean cluster mass and biased cosmology results. We show that the lensing selection bias exhibits a strong scale dependence and is approximately 20–60 per cent for ΔΣ at large scales. This selection bias largely originates from spurious member galaxies within ±20–60 $h^{-1}\, \rm Mpc$ along the line of sight, highlighting the importance of quantifying projection effects associated with the broad redshift distribution of member galaxies in photometric cluster surveys. While our results qualitatively agree withmore »those in the literature, accurate quantitative modelling of the selection bias is needed to achieve the goals of cluster lensing cosmology and will require synthetic catalogues covering a wide range of galaxy–halo connection models.« less
  5. ABSTRACT The DMASS sample is a photometric sample from the DES Year 1 data set designed to replicate the properties of the CMASS sample from BOSS, in support of a joint analysis of DES and BOSS beyond the small overlapping area. In this paper, we present the measurement of galaxy–galaxy lensing using the DMASS sample as gravitational lenses in the DES Y1 imaging data. We test a number of potential systematics that can bias the galaxy–galaxy lensing signal, including those from shear estimation, photometric redshifts, and observing conditions. After careful systematic tests, we obtain a highly significant detection of the galaxy–galaxy lensing signal, with total S/N = 25.7. With the measured signal, we assess the feasibility of using DMASS as gravitational lenses equivalent to CMASS, by estimating the galaxy-matter cross-correlation coefficient rcc. By jointly fitting the galaxy–galaxy lensing measurement with the galaxy clustering measurement from CMASS, we obtain $r_{\rm cc}=1.09^{+0.12}_{-0.11}$ for the scale cut of $4 \, h^{-1}{\rm \,\,Mpc}$ and $r_{\rm cc}=1.06^{+0.13}_{-0.12}$ for $12 \, h^{-1}{\rm \,\,Mpc}$ in fixed cosmology. By adding the angular galaxy clustering of DMASS, we obtain rcc = 1.06 ± 0.10 for the scale cut of $4 \, h^{-1}{\rm \,\,Mpc}$ and rcc = 1.03 ± 0.11 for $12 \, h^{-1}{\rm \,\,Mpc}$. The resultingmore »values of rcc indicate that the lensing signal of DMASS is statistically consistent with the one that would have been measured if CMASS had populated the DES region within the given statistical uncertainty. The measurement of galaxy–galaxy lensing presented in this paper will serve as part of the data vector for the forthcoming cosmology analysis in preparation.« less