No Arabic abstract
Understanding the physical connection between cluster galaxies and massive haloes is key to mitigating systematic uncertainties in next-generation cluster cosmology. We develop a novel method to infer the level of conformity between the stellar mass of the brightest central galaxies~(BCGs) $M_*^{BCG}$ and the satellite richness $lambda$, defined as their correlation coefficient $rho_{cc}$ at fixed halo mass, using the abundance and weak lensing of SDSS clusters as functions of $M_*^{BCG}$ and $lambda$. We detect a halo mass-dependent conformity as $rho_{cc}{=}0.60{+}0.08ln(M_h/3{times}10^{14}M_{odot}/h)$. The strong conformity successfully resolves the halo mass equality conundrum discovered in Zu et al. 2021 --- when split by $M_*^{BCG}$ at fixed $lambda$, the low and high-$M_*^{BCG}$ clusters have the same average halo mass despite having a $0.34$ dex discrepancy in average $M_*^{BCG}$. On top of the best--fitting conformity model, we develop a cluster assembly bias~(AB) prescription calibrated against the CosmicGrowth simulation, and build a conformity+AB model for the cluster weak lensing measurements. Our model predicts that with a ${sim}20%$ lower halo concentration $c$, the low-$M_*^{BCG}$ clusters are ${sim}10%$ more biased than the high-$M_*^{BCG}$ systems, in excellent agreement with the observations. We also show that the observed conformity and assembly bias are unlikely due to projection effects. Finally, we build a toy model to argue that while the early-time BCG-halo co-evolution drives the $M_*^{BCG}$-$c$ correlation, the late-time dry merger-induced BCG growth naturally produces the $M_*^{BCG}$-$lambda$ conformity despite the well-known anti-correlation between $lambda$ and $c$. Our method paves the path towards simultaneously constraining cosmology and cluster formation with future cluster surveys.
Empirical methods for connecting galaxies to their dark matter halos have become essential for interpreting measurements of the spatial statistics of galaxies. In this work, we present a novel approach for parameterizing the degree of concentration dependence in the abundance matching method. This new parameterization provides a smooth interpolation between two commonly used matching proxies: the peak halo mass and the peak halo maximal circular velocity. This parameterization controls the amount of dependence of galaxy luminosity on halo concentration at a fixed halo mass. Effectively this interpolation scheme enables abundance matching models to have adjustable assembly bias in the resulting galaxy catalogs. With the new 400 Mpc/h DarkSky Simulation, whose larger volume provides lower sample variance, we further show that low-redshift two-point clustering and satellite fraction measurements from SDSS can already provide a joint constraint on this concentration dependence and the scatter within the abundance matching framework.
One of the main predictions of excursion set theory is that the clustering of dark matter haloes only depends on halo mass. However, it has been long established that the clustering of haloes also depends on other properties, including formation time, concentration, and spin; this effect is commonly known as halo assembly bias. We use a suite of gravity-only simulations to study the dependence of halo assembly bias on cosmology; these simulations cover cosmological parameters spanning 10$sigma$ around state-of-the-art best-fitting values, including standard extensions of the $Lambda$CDM paradigm such as neutrino mass and dynamical dark energy. We find that the strength of halo assembly bias presents variations smaller than 0.05 dex across all cosmologies studied for concentration and spin selected haloes, letting us conclude that the dependence of halo assembly bias upon cosmology is negligible. We then study the dependence of galaxy assembly bias (i.e. the manifestation of halo assembly bias in galaxy clustering) on cosmology using subhalo abundance matching. We find that galaxy assembly bias also presents very small dependence upon cosmology ($sim$ 2$%$-4$%$ of the total clustering); on the other hand, we find that the dependence of this signal on the galaxy formation parameters of our galaxy model is much stronger. Taken together, these results let us conclude that the dependence of halo and galaxy assembly bias on cosmology is practically negligible.
We present evidence for halo assembly bias as a function of geometric environment. By classifying GAMA galaxy groups as residing in voids, sheets, filaments or knots using a tidal tensor method, we find that low-mass haloes that reside in knots are older than haloes of the same mass that reside in voids. This result provides direct support to theories that link strong halo tidal interactions with halo assembly times. The trend with geometric environment is reversed at large halo mass, with haloes in knots being younger than haloes of the same mass in voids. We find a clear signal of halo downsizing - more massive haloes host galaxies that assembled their stars earlier. This overall trend holds independently of geometric environment. We support our analysis with an in-depth exploration of the L-Galaxies semi-analytic model, used here to correlate several galaxy properties with three different definitions of halo formation time. We find a complex relationship between halo formation time and galaxy properties, with significant scatter. We confirm that stellar mass to halo mass ratio, specific star-formation rate and mass-weighed age are reasonable proxies of halo formation time, especially at low halo masses. Instantaneous star-formation rate is a poor indicator at all halo masses. Using the same semi-analytic model, we create mock spectral observations using complex star-formation and chemical enrichment histories, that approximately mimic GAMAs typical signal-to-noise and wavelength range. We use these mocks to assert how well potential proxies of halo formation time may be recovered from GAMA-like spectroscopic data.
Understanding the impact of halo properties beyond halo mass on the clustering of galaxies (namely galaxy assembly bias) remains a challenge for contemporary models of galaxy clustering. We explore the use of machine learning to predict the halo occupations and recover galaxy clustering and assembly bias in a semi-analytic galaxy formation model. For stellar-mass selected samples, we train a Random Forest algorithm on the number of central and satellite galaxies in each dark matter halo. With the predicted occupations, we create mock galaxy catalogues and measure the clustering and assembly bias. Using a range of halo and environment properties, we find that the machine learning predictions of the occupancy variations with secondary properties, galaxy clustering and assembly bias are all in excellent agreement with those of our target galaxy formation model. Internal halo properties are most important for the central galaxies prediction, while environment plays a critical role for the satellites. Our machine learning models are all provided in a usable format. We demonstrate that machine learning is a powerful tool for modelling the galaxy-halo connection, and can be used to create realistic mock galaxy catalogues which accurately recover the expected occupancy variations, galaxy clustering and galaxy assembly bias, imperative for cosmological analyses of upcoming surveys.
We investigate the level of galaxy assembly bias in the Sloan Digital Sky Survey (SDSS) main galaxy redshift survey using ELUCID, a state-of-the-art constrained simulation that accurately reconstructed the initial density perturbations within the SDSS volume. On top of the ELUCID haloes, we develop an extended HOD model that includes the assembly bias of central and satellite galaxies, parameterized as $mathcal{Q}_mathrm{cen}$ and $mathcal{Q}_mathrm{sat}$, respectively, to predict a suite of one- and two-point observables. In particular, our fiducial constraint employs the probability distribution of the galaxy overdensity $delta^g_8$ and the projected correlation functions of quintiles of galaxies selected by $delta^g_8$. We perform extensive tests of the efficacy of our method by fitting the same observables to mock data using both constrained and non-constrained simulations. We discover that in many cases the level of cosmic variance between the two simulations can produce biased constraints that lead to an erroneous detection of galaxy assembly bias if the non-constrained simulation is used. When applying our method to the SDSS data, the ELUCID reconstruction effectively removes an otherwise strong degeneracy between cosmic variance and galaxy assembly bias in SDSS, enabling us to derive an accurate and stringent constraint on the latter. Our fiducial ELUCID constraint, for galaxies above a stellar mass threshold $M_* = 10^{10.2},h^{-1},M_odot$, is $mathcal{Q}_mathrm{cen} = -0.06pm0.09$ and $mathcal{Q}_mathrm{sat}=0.08pm0.12$, indicating no evidence for a significant galaxy assembly bias in the local Universe probed by SDSS. Finally, our method provides a promising path to the robust modelling of the galaxy-halo connection within future spectroscopic surveys like DESI and PFS.