No Arabic abstract
Large surveys of the local Universe have shown that galaxies with different intrinsic properties, such as colour, luminosity and morphological type display a range of clustering amplitudes. Galaxies are therefore not faithful tracers of the underlying matter distribution. This modulation of galaxy clustering, called bias, contains information about the physics behind galaxy formation. It is also a systematic to be overcome before the large-scale structure of the Universe can be used as a cosmological probe. Two types of approaches have been developed to model the clustering of galaxies. The first class is empirical and filters or weights the distribution of dark matter to reproduce the measured clustering. In the second approach an attempt is made to model the physics which governs fate of baryons in order to predict the number of galaxies in dark matter haloes. I will review the development of both approaches and summarize what we have learnt about galaxy bias.
We present a comprehensive study of the chemical properties of the stellar haloes of Milky-Way mass galaxies, analysing the transition between the inner to the outer haloes. We find the transition radius between the relative dominance of the inner-halo and outer-halo stellar populations to be ~15-20 kpc for most of our haloes, similar to that inferred for the Milky Way from recent observations. While the number density of stars in the simulated inner-halo populations decreases rapidly with distance, the outer-halo populations contribute about 20-40 per cent in the fiducial solar neighborhood, in particular at the lowest metallicities. We have determined [Fe/H] profiles for our simulated haloes; they exhibit flat or mild gradients, in the range [-0.002, -0.01 ] dex/kpc. The metallicity distribution functions exhibit different features, reflecting the different assembly history of the individual stellar haloes. We find that stellar haloes formed with larger contributions from massive subgalactic systems have steeper metallicity gradients. Very metal-poor stars are mainly contributed to the halo systems by lower-mass satellites. There is a clear trend among the predicted metallicity distribution functions that a higher fraction of low-metallicity stars are found with increasing radius. These properties are consistent with the range of behaviours observed for stellar haloes of nearby galaxies.
We present a new statistical method to determine the relationship between the stellar masses of galaxies and the masses of their host dark matter haloes over the entire cosmic history from z~4 to the present. This multi-epoch abundance matching (MEAM) model self-consistently takes into account that satellite galaxies first become satellites at times earlier than they are observed. We employ a redshift-dependent parameterization of the stellar-to-halo mass relation to populate haloes and subhaloes in the Millennium simulations with galaxies, requiring that the observed stellar mass functions at several redshifts be reproduced simultaneously. Using merger trees extracted from the dark matter simulations in combination with MEAM, we predict the average assembly histories of galaxies, separating into star formation within the galaxies (in-situ) and accretion of stars (ex-situ). The peak star formation efficiency decreases with redshift from 23% at z=0 to 9% at z=4 while the corresponding halo mass increases from 10^11.8Modot to 10^12.5Modot. The star formation rate of central galaxies peaks at a redshift which depends on halo mass; for massive haloes this peak is at early cosmic times while for low-mass galaxies the peak has not been reached yet. In haloes similar to that of the Milky-Way about half of the central stellar mass is assembled after z=0.7. In low-mass haloes, the accretion of satellites contributes little to the assembly of their central galaxies, while in massive haloes more than half of the central stellar mass is formed ex-situ with significant accretion of satellites at z<2. We find that our method implies a cosmic star formation history and an evolution of specific star formation rates which are consistent with those inferred directly. We present convenient fitting functions for stellar masses, star formation rates, and accretion rates as functions of halo mass and redshift.
We investigate the role of angular momentum in the clustering of dark matter haloes. We make use of data from two high-resolution N-body simulations spanning over four orders of magnitude in halo mass, from $10^{9.8}$ to $10^{14} h^{-1} text{M}_odot$. We explore the hypothesis that mass accretion in filamentary environments alters the angular momentum of a halo, thereby driving a correlation between the spin parameter $lambda$ and the strength of clustering. However, we do not find evidence that the distribution of matter on large scales is related to the spin of haloes. We find that a halos spin is correlated with its age, concentration, sphericity, and mass accretion rate. Removing these correlations strongly affects the strength of secondary spin bias at low halo masses. We also find that high spin haloes are slightly more likely to be found near another halo of comparable mass. These haloes that are found near a comparable mass neighbour - a textit{twin} - are strongly spatially biased. We demonstrate that this textit{twin bias}, along with the relationship between spin and mass accretion rates, statistically accounts for halo spin secondary bias.
We use field-level forward models of galaxy clustering and the EFT likelihood formalism to study, for the first time for self-consistently simulated galaxies, the relations between the linear $b_1$ and second-order bias parameters $b_2$ and $b_{K^2}$. The forward models utilize all of the information available in the galaxy distribution up to a given order in perturbation theory, which allows us to infer these bias parameters with high signal-to-noise, even from relatively small volumes ($L_{rm box} = 205{rm Mpc}/h$). We consider galaxies from the IllustrisTNG simulations, and our main result is that the $b_2(b_1)$ and $b_{K^2}(b_1)$ relations obtained from gravity-only simulations for total mass selected objects are broadly preserved for simulated galaxies selected by stellar mass, star formation rate, color and black hole accretion rate. We also find good agreement between the bias relations of the simulated galaxies and a number of recent estimates for observed galaxy samples. The consistency under different galaxy selection criteria suggests that theoretical priors on these bias relations may be used to improve cosmological constraints based on observed galaxy samples. We do identify some small differences between the bias relations in the hydrodynamical and gravity-only simulations, which we show can be linked to the environmental dependence of the relation between galaxy properties and mass. We also show that the EFT likelihood recovers the value of $sigma_8$ to percent-level from various galaxy samples (including splits by color and star formation rate) and after marginalizing over 8 bias parameters. This demonstration using simulated galaxies adds to previous works based on halos as tracers, and strengthens further the potential of forward models to infer cosmology from galaxy data.
Secondary halo bias, commonly known as assembly bias, is the dependence of halo clustering on a halo property other than mass. This prediction of the Lambda-Cold Dark Matter cosmology is essential to modelling the galaxy distribution to high precision and interpreting clustering measurements. As the name suggests, different manifestations of secondary halo bias have been thought to originate from halo assembly histories. We show conclusively that this is incorrect for cluster-size haloes. We present an up-to-date summary of secondary halo biases of high-mass haloes due to various halo properties including concentration, spin, several proxies of assembly history, and subhalo properties. While concentration, spin, and the abundance and radial distribution of subhaloes exhibit significant secondary biases, properties that directly quantify halo assembly history do not. In fact, the entire assembly histories of haloes in pairs are nearly identical to those of isolated haloes. In general, a global correlation between two halo properties does not predict whether or not these two properties exhibit similar secondary biases. For example, assembly history and concentration (or subhalo abundance) are correlated for both paired and isolated haloes, but follow slightly different conditional distributions in these two cases. This results in a secondary halo bias due to concentration (or subhalo abundance), despite the lack of assembly bias in the strict sense for cluster-size haloes. Due to this complexity, caution must be exercised in using any one halo property as a proxy to study the secondary bias due to another property.