A flexible subhalo abundance matching model for galaxy clustering in redshift space

ABSTRACT We develop an extension of subhalo abundance matching (SHAM) capable of accurately reproducing the real and redshift-space clustering of galaxies in a state-of-the-art hydrodynamical simulation. Our method uses a low-resolution gravity-only simulation and it includes orphan and tidal disruption prescriptions for satellite galaxies, and a flexible amount of galaxy assembly bias. Furthermore, it includes recipes for star formation rate (SFR) based on the dark matter accretion rate. We test the accuracy of our model against catalogues of stellar-mass- and SFR-selected galaxies in the TNG300 hydrodynamic simulation. By fitting a small number of free parameters, our extended SHAM reproduces the projected correlation function and redshift-space multipoles for number densities $10^{-3} - 10^{-2}\, h^{3}{\rm Mpc}^{-3}$, at z = 1 and z = 0, and for scales r ∈ [0.3 − 20]h−1Mpc. Simultaneously, the SHAM results also retrieve the correct halo occupation distribution, the level of galaxy assembly bias, and higher order statistics present in the TNG300 galaxy catalogues. As an application, we show that our model simultaneously fits the projected correlation function of the SDSS in three disjoint stellar mass bins, with an accuracy similar to that of TNG300 galaxies. This SHAM extension can be used to get accurate clustering prediction even when using low and moderate-resolution simulations.

Multi-wavelength mock galaxy catalogs of the low-redshift Universe

We present a new suite of mock galaxy catalogs mimicking the low-redshift Universe, based on an updated halo occupation distribution (HOD) model and a scaling relation between optical properties and the neutral hydrogen (Hi) content of galaxies. Our algorithm is constrained by observations of the luminosity function and luminosity- and colour-dependent clustering of SDSS galaxies, as well as the Hi mass function and Hi-dependent clustering of massive Hi-selected galaxies in the ALFALFA survey. Mock central and satellite galaxies with realistic values of r-band luminosity, g − r and u − r colour, stellar mass and Hi mass are populated in an N-body simulation, inheriting a number of properties of the density and tidal environment of their host halos. The host halo of each central galaxy is also ‘baryonified’ with realistic spatial distributions of stars as well as hot and cold gas, along with the corresponding rotation curve. Our default HOD assumes that galaxy properties are a function of group halo mass alone, and can optionally include effects such as galactic conformity and colour-dependent galaxy assembly bias. The mocks predict the relation between the stellar mass and Hi mass of massive Hi galaxies, as well as the 2-point cross-correlation function of spatially co-located optical and Hi-selected samples. They enable novel null tests for galaxy assembly bias, provide predictions for the Hi velocity width function, and clarify the origin and universality of the radial acceleration relation in the ΛCDM framework.

Building a digital twin of a luminous red galaxy spectroscopic survey: galaxy properties and clustering covariance

ABSTRACT Upcoming surveys will use a variety of galaxy selections to map the large-scale structure of the Universe. It is important to make accurate predictions for the properties and clustering of such galaxies, including the errors on these statistics. Here, we describe a novel technique which uses the semi-analytical model of galaxy formation galform, embedded in the high-resolution N-body Planck-Millennium simulation, to populate a thousand halo catalogues generated using the Parallel-PM N-body glam code. Our hybrid scheme allows us to make clustering predictions on scales that cannot be modelled in the original N-body simulation. We focus on luminous red galaxies (LRGs) selected in the redshift range z = 0.6 − 1 from the galform output using similar colour-magnitude cuts in the r, z, and W1 bands to those that will be applied in the Dark Energy Spectroscopic Instrument (DESI) survey, and call this illustrative sample ‘DESI-like’ LRGs. We find that the LRG-halo connection is non-trivial, leading to the prediction of a non-standard halo occupation distribution; in particular, the occupation of central galaxies does not reach unity for the most massive haloes, and drops with increasing mass. The glam catalogues reproduce the abundance and clustering of the LRGs predicted by galform. We use the glam mocks to compute the covariance matrices for the two-point correlation function and power spectrum of the LRGs and their background dark matter density field, revealing important differences. We also make predictions for the linear-growth rate and the baryon acoustic oscillations distances at z = 0.6, 0.74, and 0.93. All ‘DESI-like’ LRG catalogues are made publicly available.

The galaxy–halo connection of emission-line galaxies in IllustrisTNG

ABSTRACT We employ the hydrodynamical simulation IllustrisTNG-300-1 to explore the halo occupation distribution (HOD) and environmental dependence of luminous star-forming emission-line galaxies (ELGs) at z ∼ 1. Such galaxies are key targets for current and upcoming cosmological surveys. We select model galaxies through cuts in colour–colour space allowing for a direct comparison with the Extended Baryon Oscillation Spectroscopic Survey and the Dark Energy Spectroscopic Instrument (DESI) surveys and then compare them with galaxies selected based on specific star formation rate (sSFR) and stellar mass. We demonstrate that the ELG populations are twice more likely to reside in lower density regions (sheets) compared with the mass-selected populations and twice less likely to occupy the densest regions of the cosmic web (knots). We also show that the colour-selected and sSFR-selected ELGs exhibit very similar occupation and clustering statistics, finding that the agreement is best for lower redshifts. In contrast with the mass-selected sample, the occupation of haloes by a central ELG peaks at ∼20 per cent. We furthermore explore the dependence of the HOD and the autocorrelation on environment, noticing that at fixed halo mass, galaxies in high-density regions cluster about 10 times more strongly than low-density ones. This result suggests that we should model carefully the galaxy–halo relation and implement assembly bias effects into our models (estimated at ∼4 per cent of the clustering of the DESI colour-selected sample at z = 0.8). Finally, we apply a simple mock recipe to recover the clustering on large scales (r ≳ 1 Mpc h−1) to within 1 per cent by augmenting the HOD model with an environment dependence, demonstrating the power of adopting flexible population models.

The Completed SDSS-IV Extended Baryon Oscillation Spectroscopic Survey: N-body Mock Challenge for Galaxy Clustering Measurements

We develop a series of N-body data challenges, functional to the final analysis of the extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 16 (DR16) galaxy sample. The challenges are primarily based on high-fidelity catalogs constructed from the Outer Rim simulation – a large box size realization (3h−1Gpc) characterized by an unprecedented combination of volume and mass resolution, down to 1.85 · 109h−1M⊙. We generate synthetic galaxy mocks by populating Outer Rim halos with a variety of halo occupation distribution (HOD) schemes of increasing complexity, spanning different redshift intervals. We then assess the performance of three complementary redshift space distortion (RSD) models in configuration and Fourier space, adopted for the analysis of the complete DR16 eBOSS sample of Luminous Red Galaxies (LRGs). We find all the methods mutually consistent, with comparable systematic errors on the Alcock-Paczynski parameters and the growth of structure, and robust to different HOD prescriptions – thus validating the robustness of the models and the pipelines used for the baryon acoustic oscillation (BAO) and full shape clustering analysis. In particular, all the techniques are able to recover α∥ and α⊥ to within $0.9\%$, and fσ8 to within $1.5\%$. As a by-product of our work, we are also able to gain interesting insights on the galaxy-halo connection. Our study is relevant for the final eBOSS DR16 ‘consensus cosmology’, as the systematic error budget is informed by testing the results of analyses against these high-resolution mocks. In addition, it is also useful for future large-volume surveys, since similar mock-making techniques and systematic corrections can be readily extended to model for instance the Dark Energy Spectroscopic Instrument (DESI) galaxy sample.

The Clustering of DESI-like Luminous Red Galaxies Using Photometric Redshifts

We present measurements of the redshift-dependent clustering of a DESI-like luminous red galaxy (LRG) sample selected from the Legacy Survey imaging dataset, and use the halo occupation distribution (HOD) framework to fit the clustering signal. The photometric LRG sample in this study contains 2.7 million objects over the redshift range of 0.4 < z < 0.9 over 5655 sq. degrees. We have developed new photometric redshift (photo-z) estimates using the Legacy Survey DECam and WISE photometry, with σNMAD = 0.02 precision for LRGs. We compute the projected correlation function using new methods that maximize signal-to-noise while incorporating redshift uncertainties. We present a novel algorithm for dividing irregular survey geometries into equal-area patches for jackknife resampling. For a 5-parameter HOD model fit using the MultiDark halo catalog, we find that there is little evolution in HOD parameters except at the highest-redshifts. The inferred large-scale structure bias is largely consistent with constant clustering amplitude over time. In an appendix, we explore limitations of MCMC fitting using stochastic likelihood estimates resulting from applying HOD methods to N-body catalogs, and present a new technique for finding best-fit parameters in this situation. Accompanying this paper we have released the PRLS (Photometric Redshifts for the Legacy Surveys) catalog of photo-z’s obtained by applying the methods used in this work to the full Legacy Survey Data Release 8 dataset. This catalog provides accurate photometric redshifts for objects with z < 21 over more than 16,000 square degrees of sky.

The Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: exploring the halo occupation distribution model for emission line galaxies

ABSTRACT We study the modelling of the halo occupation distribution (HOD) for the eBOSS DR16 emission line galaxies (ELGs). Motivated by previous theoretical and observational studies, we consider different physical effects that can change how ELGs populate haloes. We explore the shape of the average HOD, the fraction of satellite galaxies, their probability distribution function (PDF), and their density and velocity profiles. Our baseline HOD shape was fitted to a semi-analytical model of galaxy formation and evolution, with a decaying occupation of central ELGs at high halo masses. We consider Poisson and sub/super-Poissonian PDFs for satellite assignment. We model both Navarro–Frenk–White and particle profiles for satellite positions, also allowing for decreased concentrations. We model velocities with the virial theorem and particle velocity distributions. Additionally, we introduce a velocity bias and a net infall velocity. We study how these choices impact the clustering statistics while keeping the number density and bias fixed to that from eBOSS ELGs. The projected correlation function, wp, captures most of the effects from the PDF and satellites profile. The quadrupole, ξ2, captures most of the effects coming from the velocity profile. We find that the impact of the mean HOD shape is subdominant relative to the rest of choices. We fit the clustering of the eBOSS DR16 ELG data under different combinations of the above assumptions. The catalogues presented here have been analysed in companion papers, showing that eBOSS RSD+BAO measurements are insensitive to the details of galaxy physics considered here. These catalogues are made publicly available.

The completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: N-body mock challenge for the quasar sample

ABSTRACT The growth rate and expansion history of the Universe can be measured from large galaxy redshift surveys using the Alcock–Paczynski effect. We validate the Redshift Space Distortion models used in the final analysis of the Sloan Digital Sky Survey (SDSS) extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 16 quasar clustering sample, in configuration and Fourier space, using a series of halo occupation distribution mock catalogues generated using the OuterRim N-body simulation. We test three models on a series of non-blind mocks, in the OuterRim cosmology, and blind mocks, which have been rescaled to new cosmologies, and investigate the effects of redshift smearing and catastrophic redshifts. We find that for the non-blind mocks, the models are able to recover fσ8 to within 3 per cent and α∥ and α⊥ to within 1 per cent. The scatter in the measurements is larger for the blind mocks, due to the assumption of an incorrect fiducial cosmology. From this mock challenge, we find that all three models perform well, with similar systematic errors on fσ8, α∥, and α⊥ at the level of $\sigma _{f\sigma _8}=0.013$, $\sigma _{\alpha _\parallel }=0.012$, and $\sigma _{\alpha _\bot }=0.008$. The systematic error on the combined consensus is $\sigma _{f\sigma _8}=0.011$, $\sigma _{\alpha _\parallel }=0.008$, and $\sigma _{\alpha _\bot }=0.005$, which is used in the final DR16 analysis. For baryon acoustic oscillation fits in configuration and Fourier space, we take conservative systematic errors of $\sigma _{\alpha _\parallel }=0.010$ and $\sigma _{\alpha _\bot }=0.007$.

The completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: GLAM-QPM mock galaxy catalogues for the emission line galaxy sample

ABSTRACT We present 2000 mock galaxy catalogues for the analysis of baryon acoustic oscillations (BAOs) in the Emission Line Galaxy (ELG) sample of the extended Baryon Oscillation Spectroscopic Survey Data Release 16 (eBOSS DR16). Each mock catalogue has a number density of $6.7 \times 10^{-4} h^3 \rm Mpc^{-3}$, covering a redshift range from 0.6 to 1.1. The mocks are calibrated to small-scale eBOSS ELG clustering measurements at scales of $\lesssim 30\, h^{-1}$Mpc. The mock catalogues are generated using a combination of GaLAxy Mocks (GLAM) simulations and the quick particle-mesh (QPM) method. GLAM simulations are used to generate the density field, which is then assigned dark matter haloes using the QPM method. Haloes are populated with galaxies using a halo occupation distribution. The resulting mocks match the survey geometry and selection function of the data, and have slightly higher number density that allows room for systematic analysis. The large-scale clustering of mocks at the BAO scale is consistent with data and we present the correlation matrix of the mocks.