No Arabic abstract
Observational systematics complicate comparisons with theoretical models limiting understanding of galaxy evolution. In particular, different empirical determinations of the stellar mass function imply distinct mappings between the galaxy and halo masses, leading to diverse galaxy evolutionary tracks. Using our state-of-the-art STatistical sEmi-Empirical modeL, STEEL, we show fully self-consistent models capable of generating galaxy growth histories that simultaneously and closely agree with the latest data on satellite richness and star-formation rates at multiple redshifts and environments. Central galaxy histories are generated using the central halo mass tracks from state-of-the-art statistical dark matter accretion histories coupled to abundance matching routines. We show that too flat high-mass slopes in the input stellar-mass-halo-mass relations as predicted by previous works, imply non-physical stellar mass growth histories weaker than those implied by satellite accretion alone. Our best-fit models reproduce the satellite distributions at the largest masses and highest redshifts probed, the latest data on star formation rates and its bi-modality in the local Universe, and the correct fraction of ellipticals. Our results are important to predict robust and self-consistent stellar-mass-halo-mass relations and to generate reliable galaxy mock catalogues for the next generations of extra-galactic surveys such as Euclid and LSST.
Connecting the observed rest-ultraviolet (UV) luminosities of high-$z$ galaxies to their intrinsic luminosities (and thus star formation rates) requires correcting for the presence of dust. We bypass a common dust-correction approach that uses empirical relationships between infrared (IR) emission and UV colours, and instead augment a semi-empirical model for galaxy formation with a simple -- but self-consistent -- dust model and use it to jointly fit high-$z$ rest-UV luminosity functions (LFs) and colour-magnitude relations ($M_{mathrm{UV}}$-$beta$). In doing so, we find that UV colours evolve with redshift (at fixed UV magnitude), as suggested by observations, even in cases without underlying evolution in dust production, destruction, absorption, or geometry. The observed evolution in our model arises due to the reduction in the mean stellar age and rise in specific star formation rates with increasing $z$. The UV extinction, $A_{mathrm{UV}}$, evolves similarly with redshift, though we find a systematically shallower relation between $A_{mathrm{UV}}$ and $M_{mathrm{UV}}$ than that predicted by IRX-$beta$ relationships derived from $z sim 3$ galaxy samples. Finally, assuming that high $1600 r{A}$ transmission ($gtrsim 0.6$) is a reliable LAE indicator, modest scatter in the effective dust surface density of galaxies can explain the evolution both in $M_{mathrm{UV}}$-$beta$ and LAE fractions. These predictions are readily testable by deep surveys with the James Webb Space Telescope.
We present STEEL a STatistical sEmi-Empirical modeL designed to probe the distribution of satellite galaxies in groups and clusters. Our fast statistical methodology relies on tracing the abundances of central and satellite haloes via their mass functions at all cosmic epochs with virtually no limitation on cosmic volume and mass resolution. From mean halo accretion histories and subhalo mass functions the satellite mass function is progressively built in time via abundance matching techniques constrained by number densities of centrals in the local Universe. By enforcing dynamical merging timescales as predicted by high-resolution N-body simulations, we obtain satellite distributions as a function of stellar mass and halo mass consistent with current data. We show that stellar stripping, star formation, and quenching play all a secondary role in setting the number densities of massive satellites above $M_*gtrsim 3times 10^{10}, M_{odot}$. We further show that observed star formation rates used in our empirical model over predict low-mass satellites below $M_*lesssim 3times 10^{10}, M_{odot}$, whereas, star formation rates derived from a continuity equation approach yield the correct abundances similar to previous results for centrals.
We present a semi-analytic model of satellite galaxies, SatGen, which can generate large samples of satellite populations for a host halo of desired mass, redshift, and assembly history. The model combines dark-matter halo merger trees, empirical relations for the galaxy-halo connection, and analytic prescriptions for tidal effects, dynamical friction, and ram pressure stripping. SatGen emulates cosmological zoom-in hydro-simulations in certain aspects. Satellites can reside in cored or cuspy DM subhaloes, depending on the halo response to baryonic physics that can be formulated from hydro-simulations and physical modeling. The subhalo profile and the stellar mass and size of a satellite evolves depending on its tidal mass loss and initial structure. The host galaxy can include a baryonic disc and a stellar bulge, each described by a density profile that allows analytic orbit integration. SatGen complements simulations by propagating the effect of halo response found in simulated field galaxies to satellites (not properly resolved in simulations) and outperforms simulations by sampling the halo-to-halo variance of satellite statistics and overcoming artificial disruption due to insufficient resolution. As a first application, we use the model to study satellites of Milky Way sized hosts, making it emulate simulations of bursty star formation and of smooth star formation, respectively, and to experiment with a disc potential in the host halo. Our model reproduces the observed satellite statistics reasonably well. Different physical recipes make a difference in satellite abundance and spatial distribution at the 25% level, not large enough to be distinguished by current observations given the halo-to-halo variance. The MW disc depletes satellites by 20% and has a subtle effect of diversifying the internal structure of satellites, important for alleviating certain small-scale problems.
We present a simple semi-numerical model designed to explore black hole growth and galaxy evolution. This method builds on a previous model for black hole accretion that uses a semi-numerical galaxy formation model and universal Eddington ratio distribution to describe the full AGN population by independently connecting galaxy and AGN growth to the evolution of the host dark matter halos. We fit observed X-ray luminosity functions up to a redshift of z ~ 4, as well as investigate the evolution of the Eddington ratio distributions. We find that the Eddington ratio distribution evolves with redshift such that the slope of the low-Eddington accretion rate distribution increases with cosmic time, consistent with the behavior predicted in hydrodynamical simulations for galaxies with different gas fractions. We also find that the evolution of our average Eddington ratio is correlated with observed star formation histories, supporting a picture in which black holes and galaxies evolve together in a global sense. We further confirm the impact of luminosity limits on observed galaxy and halo properties by applying selection criteria to our fiducial model and comparing to surveys across a wide range of redshifts.
We explore the galaxy-galaxy merger rate with the empirical model for galaxy formation, Emerge. On average, we find that between $2$ per cent and $20$ per cent of massive galaxies ($log_{10}(m_{*}/M_{odot}) geq 10.3$) will experience a major merger per Gyr. Our model predicts galaxy merger rates that do not scale as a power-law with redshift when selected by descendant stellar mass, and exhibit a clear stellar mass and mass-ratio dependence. Specifically, major mergers are more frequent at high masses and at low redshift. We show mergers are significant for the stellar mass growth of galaxies $log_{10}(m_{*}/M_{odot}) gtrsim 11.0$. For the most massive galaxies major mergers dominate the accreted mass fraction, contributing as much as $90$ per cent of the total accreted stellar mass. We reinforce that these phenomena are a direct result of the stellar-to-halo mass relation, which results in massive galaxies having a higher likelihood of experiencing major mergers than low mass galaxies. Our model produces a galaxy pair fraction consistent with recent observations, exhibiting a form best described by a power-law exponential function. Translating these pair fractions into merger rates results in an inaccurate prediction compared to the model intrinsic values when using published observation timescales. We find the pair fraction can be well mapped to the intrinsic merger rate by adopting an observation timescale that decreases linearly with redshift as $T_{mathrm{obs}} = -0.36(1+z)+2.39$ [Gyr], assuming all observed pairs merge by $z=0$.