N-body simulations suggest that the substructures that survive inside dark matter haloes follow universal distributions in mass and radial number density. We demonstrate that a simple analytical ...model can explain these subhalo distributions as resulting from tidal stripping which increasingly reduces the mass of subhaloes with decreasing halocentric distance. As a starting point, the spatial distribution of subhaloes of any given infall mass is shown to be largely indistinguishable from the overall mass distribution of the host halo. Using a physically motivated statistical description of the amount of mass stripped from individual subhaloes, the model fully describes the joint distribution of subhaloes in final mass, infall mass and radius. As a result, it can be used to predict several derived distributions involving combinations of these quantities including, but not limited to, the universal subhalo mass function, the subhalo spatial distribution, the gravitational lensing profile, the dark matter annihilation radiation profile and boost factor. This model clarifies a common confusion when comparing the spatial distributions of galaxies and subhaloes, the so-called anti-bias, as a simple selection effect. We provide a python code subgen for populating haloes with subhaloes at http://icc.dur.ac.uk/data/.
We use a suite of cosmological simulations to study the mass–concentration–redshift relation, c(M, z), of dark matter haloes. Our simulations include standard Λ-cold dark matter (CDM) models, and ...additional runs with truncated power spectra, consistent with a thermal warm dark matter (WDM) scenario. We find that the mass profiles of CDM and WDM haloes are self-similar and well approximated by the Einasto profile. The c(M, z) relation of CDM haloes is monotonic: concentrations decrease with increasing virial mass at fixed redshift, and decrease with increasing redshift at fixed mass. The mass accretion histories (MAHs) of CDM haloes are also scale-free, and can be used to infer concentrations directly. These results do not apply to WDM haloes: their MAHs are not scale-free because of the characteristic scale imposed by the power spectrum suppression. Further, the WDM c(M, z) relation is non-monotonic: concentrations peak at a mass scale dictated by the truncation scale, and decrease at higher and lower masses. We show that the assembly history of a halo can still be used to infer its concentration, provided that the total mass of its progenitors is considered (the ‘collapsed mass history’; CMH), rather than just that of its main ancestor. This exploits the scale-free nature of CMHs to derive a simple scaling that reproduces the mass–concentration–redshift relation of both CDM and WDM haloes over a vast range of halo masses and redshifts. Our model therefore provides a robust account of the mass, redshift, cosmology and power spectrum dependence of dark matter halo concentrations.
The mass of the dark matter halo of the Milky Way can be estimated by fitting analytical models to the phase-space distribution of dynamical tracers. We test this approach using realistic mock ...stellar haloes constructed from the Aquarius N-body simulations of dark matter haloes in the Λ cold dark matter cosmology. We extend the standard treatment to include a Navarro–Frenk–White potential and use a maximum likelihood method to recover the parameters describing the simulated haloes from the positions and velocities of their mock halo stars. We find that the estimate of halo mass is highly correlated with the estimate of halo concentration. The best-fitting halo masses within the virial radius, R
200, are biased, ranging from a 40 per cent underestimate to a 5 per cent overestimate in the best case (when the tangential velocities of the tracers are included). There are several sources of bias. Deviations from dynamical equilibrium can potentially cause significant bias; deviations from spherical symmetry are relatively less important. Fits to stars at different galactocentric radii can give different mass estimates. By contrast, the model gives good constraints on the mass within the half-mass radius of tracers even when restricted to tracers within 60 kpc. The recovered velocity anisotropies of tracers, β, are biased systematically, but this does not affect other parameters if tangential velocity data are used as constraints.
The estimated stellar masses of galaxies are widely used to characterize how the galaxy population evolves over cosmic time. If stellar masses can be estimated in a robust manner, free from any bias, ...global diagnostics such as the stellar mass function can be used to constrain the physics of galaxy formation. We explore how galaxy stellar masses, estimated by fitting broad-band spectral energy distributions (SEDs) with stellar population models, can be biased as a result of commonly adopted assumptions for the star formation and chemical enrichment histories, recycled fractions and dust attenuation curves of galaxies. We apply the observational technique of broad-band SED fitting to model galaxy SEDs calculated by the theoretical galaxy formation model GALFORM, isolating the effect of each of these assumptions. We find that, averaged over the entire galaxy population, the common assumption of exponentially declining star formation histories does not, by itself, adversely affect stellar mass estimation. However, we also show that this result does not hold when considering galaxies that have undergone a recent burst of star formation. We show that fixing the metallicity in SED fitting or using sparsely sampled metallicity grids can introduce mass-dependent systematics into stellar mass estimates. We find that the common assumption of a star-dust geometry corresponding to a uniform foreground dust screen can cause the stellar masses of dusty model galaxies to be significantly underestimated. Finally, we show that stellar mass functions recovered by applying SED fitting to model galaxies at high redshift can differ significantly in both shape and normalization from the intrinsic mass functions predicted by a given model. In particular, the effects of dust can reduce the normalization at the high-mass end by up to 0.6 dex in some cases. Given these differences, our methodology of using stellar masses estimated from model galaxy SEDs offers a new, self-consistent way to compare model predictions with observations. We conclude that great care should be taken when comparing theoretical galaxy formation models to observational results based on the estimated stellar masses of high-redshift galaxies.
We compare global predictions from the eagle hydrodynamical simulation, and two semi-analytic (SA) models of galaxy formation, l-galaxies and galform. All three models include the key physical ...processes for the formation and evolution of galaxies and their parameters are calibrated against a small number of observables at z ≈ 0. The two SA models have been applied to merger trees constructed from the eagle dark matter only simulation. We find that at z ≤ 2, both the galaxy stellar mass functions for stellar masses M
* < 1010.5 M⊙ and the median specific star formation rates (sSFRs) in the three models agree to better than 0.4 dex. The evolution of the sSFR predicted by the three models closely follows the mass assembly history of dark matter haloes. In both eagle and l-galaxies there are more central passive galaxies with M
* < 109.5 M⊙ than in galform. This difference is related to galaxies that have entered and then left a larger halo and which are treated as satellites in galform. In the range 0 < z < 1, the slope of the evolution of the star formation rate density in eagle is a factor of ≈1.5 steeper than for the two SA models. The median sizes for galaxies with M
* > 109.5 M⊙ differ in some instances by an order of magnitude, while the stellar mass–size relation in eagle is a factor of ≈2 tighter than for the two SA models. Our results suggest the need for a revision of how SA models treat the effect of baryonic self-gravity on the underlying dark matter. The treatment of gas flows in the models needs to be revised based on detailed comparison with observations to understand in particular the evolution of the stellar mass–metallicity relation.
Abstract
It is now possible for hydrodynamical simulations to reproduce a representative galaxy population. Accordingly, it is timely to assess critically some of the assumptions of traditional ...semi-analytic galaxy formation models. We use the eagle simulations to assess assumptions built into the galform semi-analytic model, focusing on those relating to baryon cycling, angular momentum and feedback. We show that the assumption in galform that newly formed stars have the same specific angular momentum as the total disc leads to a significant overestimate of the total stellar specific angular momentum of disc galaxies. In eagle, stars form preferentially out of low-specific angular momentum gas in the interstellar medium due to the assumed gas density threshold for stars to form, leading to more realistic galaxy sizes. We find that stellar mass assembly is similar between galform and eagle but that the evolution of gas properties is different, with various indications that the rate of baryon cycling in eagle is slower than is assumed in galform. Finally, by matching individual galaxies between eagle and galform, we find that an artificial dependence of active galactic nucleus feedback and gas infall rates on halo mass-doubling events in galform drives most of the scatter in stellar mass between individual objects. Put together our results suggest that the galform semi-analytic model can be significantly improved in light of recent advances.
ABSTRACT
A number of Local Group dwarf galaxies are known to have two spatially segregated stellar metallicity populations, a centrally concentrated metal-rich population, and a more extended ...metal-poor population. In this work we discuss mechanisms that lead to the formation of two spatially segregated metallicity populations. Using a set of high-resolution hydrodynamical simulations of Local Group-like environments, we select a sample of satellite and field galaxies, spanning the stellar mass range 106–109 M⊙, that exhibit bimodality in their metallicity distributions. Among those, we identify a subsample with a strong spatial segregation in the two populations. We find three distinct mechanisms for their formation. In field dwarfs and in a small fraction of satellites, a merger causes the metal-poor stars to migrate to larger radii and encourages the available gas to sink to the centre of the dwarf. Most of the gas is subsequently blown out of the halo through star formation feedback, but the remaining gas is consumed in the formation of a metal-rich population. In the exclusive case of satellites that have retained some of their gas at infall, it is the compression of this gas by ram pressure near pericentre that triggers the formation of metal-rich stars, whilst simultaneously preventing star formation at larger radii through stripping. Additionally, in a small number of field and satellite dwarfs, interactions with gaseous filaments and other galaxies can result in the formation of a metal-rich population. Regardless of the formation mechanism, a history of mergers typically enhances the spatial segregation.