We study the dependence of the galaxy content of dark matter halos on large-scale environment and halo formation time using semi-analytic galaxy models applied to the Millennium simulation. We ...analyze subsamples of halos at the extremes of these distributions and measure the occupation functions for the galaxies they host. We find distinct differences among these occupation functions. The main effect with environment is that central galaxies (and in one model, also the satellites) in denser regions start populating lower-mass halos. A similar, but significantly stronger, trend exists with halo age, where early-forming halos are more likely to host central galaxies at lower halo mass. We discuss the origin of these trends and the connection to the stellar mass-halo mass relation. We find that, at fixed halo mass, older halos and to some extent also halos in dense environments tend to host more massive galaxies. Additionally, we see a reverse trend for the occupation of satellite galaxies where early-forming halos have fewer satellites, likely due to having more time for them to merge with the central galaxy. We describe these occupancy variations in terms of the changes in the occupation function parameters, which can aid in constructing realistic mock galaxy samples. Finally, we study the corresponding galaxy auto- and cross-correlation functions of the different samples and elucidate the impact of assembly bias on galaxy clustering. Our results can inform theoretical modeling of galaxy assembly bias and attempts to detect it in the real universe.
ABSTRACT
We carry out ‘full-physics’ hydrodynamical simulations of galaxy formation in the normal-branch Dvali–Gabadadze–Porrati (nDGP) braneworld model using a new modified version of the arepo code ...and the IllustrisTNG galaxy formation model. We simulate two nDGP models (N5 and N1) that represent, respectively, weak and moderate departures from general relativity (GR), in boxes of sizes $62$ and $25\, h^{-1}{\rm Mpc}$ using 2 × 5123 dark matter particles and initial gas cells. This allows us to explore, for the first time, the impact of baryonic physics on galactic scales in braneworld models of modified gravity and to make predictions on the stellar content of dark matter haloes and galaxy evolution through cosmic time in these models. We find significant differences between the GR and nDGP models in the power spectra and correlation functions of gas, stars and dark matter of up to ∼25 per cent on large scales. Similar to their impact in the standard cosmological model (Λ cold dark matter), baryonic effects can have a significant influence over the clustering of the overall matter distribution, with a sign that depends on scale. Studying the degeneracy between modified gravity and galactic feedback in these models, we find that these two physical effects on matter clustering can be cleanly disentangled, allowing for a method to accurately predict the matter power spectrum with baryonic effects included, without having to run hydrodynamical simulations. Depending on the braneworld model, we find differences compared with GR of up to ∼15 per cent in galaxy properties such as the stellar-to-halo-mass ratio, galaxy stellar mass function, gas fraction, and star formation rate density. The amplitude of the fifth force is reduced by the presence of baryons in the very inner part of haloes, but this reduction quickly becomes negligible above ∼0.1 times the halo radius.
ABSTRACT
We present a study of galaxy mergers up to z = 10 using the Planck Millennium cosmological dark matter simulation and the GALFORM semi-analytical model of galaxy formation. Utilizing the ...full 800 Mpc3 volume of the simulation, we studied the statistics of galaxy mergers in terms of merger rates and close pair fractions. We predict that merger rates begin to drop rapidly for high-mass galaxies (M* > 1011.3–1010.5 M⊙ for z = 0–4), as a result of the exponential decline in the galaxy stellar mass function. The predicted merger rates for massive galaxies (M* > 1010 M⊙) increase and then turn over with increasing redshift, by z = 3.5, in disagreement with hydrodynamical simulations and semi-empirical models. In agreement with most other models and observations, we find that close pair fractions flatten or turn over at some redshift (dependent on the mass selection). We conduct an extensive comparison of close pair fractions, and highlight inconsistencies among models, but also between different observations. We provide a fitting formula for the major merger time-scale for close galaxy pairs, in which the slope of the stellar mass dependence is redshift dependent. This is in disagreement with previous theoretical results that implied a constant slope. Instead, we find a weak redshift dependence only for massive galaxies (M* > 1010 M⊙): in this case the merger time-scale varies approximately as $M_*^{-0.55}$. We find that close pair fractions and merger time-scales depend on the maximum projected separation as $r_\mathrm{max}^{1.32}$, in agreement with observations of small-scale clustering of galaxies.
Abstract
We forecast the linear bias for Hα-emitting galaxies at high redshift. To simulate a Euclid-like and a WFIRST-like survey, we place galaxies into a large-volume dark matter halo lightcone ...by sampling a library of luminosity-dependent halo occupation distributions (HODs), which is constructed using a physically motivated galaxy formation model. We calibrate the dust attenuation in the lightcones such that they are able to reproduce the Hα luminosity function or the Hα cumulative number counts. The angle-averaged galaxy correlation function is computed for each survey in redshift slices of width Δ z = 0.2. In each redshift bin the linear bias can be fitted with a single, scale-independent value that increases with increasing redshift. Fitting for the evolution of linear bias with redshift, we find that our Euclid-like and WFIRST-like surveys are both consistent within error with the relation b(z) = 0.7z + 0.7. Our bias forecasts are consistent with bias measurements from the HiZELS survey. We find that the Euclid-like and WFIRST-like surveys yield linear biases that are broadly consistent within error, most likely due to the HOD for the WFIRST-like survey having a steeper power-law slope towards larger halo masses.
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.
ABSTRACT
We use the GALFORM semi-analytical model of galaxy formation and the Planck-Millennium simulation to investigate the origins of stellar mass in galaxies and their spheroids. We compare the ...importance of mergers and disc instabilities, as well as the starbursts that they trigger. We find that the fraction of galaxy stellar mass formed ex situ (i.e. through mergers; fex) increases sharply from M* = 1011 M⊙ upwards, reaching 80 per cent at M* = 1011.3 M⊙. The massive end of the fex–M* relation does not evolve with redshift, in disagreement with other models. For low-mass galaxies we find larger ex situ contributions at z = 0 than in other models (7–12 per cent), with a decrease towards higher redshifts. Major mergers contribute roughly half of the ex situ mass, with minor mergers and smooth accretion of satellites both accounting for ≈25 per cent, almost independent of stellar mass and redshift. Mergers dominate in building up high-mass (M*, sph > 1011 M⊙) and low-mass (M*, sph < 108.5 M⊙) spheroids. Disc instabilities and their associated starbursts dominate for intermediate-mass spheroids (108.5 < M*, sph < 1011 M⊙) at z = 0. The mass regime where pseudo-bulges dominate is in agreement with observed pseudo-bulge fractions, but the peak value in the pseudo-bulge fraction predicted by GALFORM is likely too high. Starbursts induced by disc instabilities are the dominant channel for spheroid growth at all redshifts, while merger-induced starbursts are relatively negligible, except at very high redshifts (z > 5).
A unified multiwavelength model of galaxy formation Lacey, Cedric G; Baugh, Carlton M; Frenk, Carlos S ...
Monthly notices of the Royal Astronomical Society,
11/2016, Letnik:
462, Številka:
4
Journal Article
Recenzirano
Odprti dostop
We present a new version of the galform semi-analytical model of galaxy formation. This brings together several previous developments of galform into a single unified model, including a different ...initial mass function (IMF) in quiescent star formation and in starbursts, feedback from active galactic nuclei suppressing gas cooling in massive haloes, and a new empirical star formation law in galaxy discs based on their molecular gas content. In addition, we have updated the cosmology, introduced a more accurate treatment of dynamical friction acting on satellite galaxies, and updated the stellar population model. The new model is able to simultaneously explain both the observed evolution of the K-band luminosity function and stellar mass function, and the number counts and redshift distribution of sub-mm galaxies selected at 850 μm. This was not previously achieved by a single physical model within the Λcold dark matter framework, but requires having an IMF in starbursts that is somewhat top-heavy. The new model is tested against a wide variety of observational data covering wavelengths from the far-UV to sub-mm, and redshifts from z = 0 to 6, and is found to be generally successful. These observations include the optical and near-infrared (IR) luminosity functions, H i mass function, fraction of early type galaxies, Tully–Fisher, metallicity–luminosity and size–luminosity relations at z = 0, as well as far-IR number counts, and far-UV luminosity functions at z ∼ 3–6. Discrepancies are, however, found in galaxy sizes and metallicities at low luminosities, and in the abundance of low-mass galaxies at high-z, suggesting the need for a more sophisticated model of supernova feedback.
ABSTRACT
We apply the marked correlation function test proposed by Armijo et al. (Paper I) to samples of luminous red galaxies (LRGs) from the final data release of the Sloan Digital Sky Survey ...(SDSS) III. The test assigns a density-dependent mark to galaxies in the estimation of the projected marked correlation function. Two gravity models are compared: general relativity (GR) and $f(R)$ gravity. We build mock catalogues which, by construction, reproduce the measured galaxy number density and two-point correlation function of the LRG samples, using the halo occupation distribution model (HOD). A range of HOD models give acceptable fits to the observational constraints, and this uncertainty is fed through to the error in the predicted marked correlation functions. The uncertainty from the HOD modelling is comparable to the sample variance for the SDSS-III LRG samples. Our analysis shows that current galaxy catalogues are too small for the test to distinguish a popular $f(R)$ model from GR. However, upcoming surveys with a better measured galaxy number density and smaller errors on the two-point correlation function, or a better understanding of galaxy formation, may allow our method to distinguish between viable gravity models.