No Arabic abstract
We present a comparison of nine galaxy formation models, eight semi-analytical and one halo occupation distribution model, run on the same underlying cold dark matter simulation (cosmological box of co-moving width 125$h^{-1}$ Mpc, with a dark-matter particle mass of $1.24times 10^9 h^{-1}$ Msun) and the same merger trees. While their free parameters have been calibrated to the same observational data sets using two approaches, they nevertheless retain some memory of any previous calibration that served as the starting point (especially for the manually-tuned models). For the first calibration, models reproduce the observed z = 0 galaxy stellar mass function (SMF) within 3-{sigma}. The second calibration extended the observational data to include the z = 2 SMF alongside the z~0 star formation rate function, cold gas mass and the black hole-bulge mass relation. Encapsulating the observed evolution of the SMF from z = 2 to z = 0 is found to be very hard within the context of the physics currently included in the models. We finally use our calibrated models to study the evolution of the stellar-to-halo mass (SHM) ratio. For all models we find that the peak value of the SHM relation decreases with redshift. However, the trends seen for the evolution of the peak position as well as the mean scatter in the SHM relation are rather weak and strongly model dependent. Both the calibration data sets and model results are publicly available.
We present a comparison of the observed evolving galaxy stellar mass functions with the predictions of eight semi-analytic models and one halo occupation distribution model. While most models are able to fit the data at low redshift, some of them struggle to simultaneously fit observations at high redshift. We separate the galaxies into passive and star-forming classes and find that several of the models produce too many low-mass star-forming galaxies at high redshift compared to observations, in some cases by nearly a factor of 10 in the redshift range $2.5 < z < 3.0$. We also find important differences in the implied mass of the dark matter haloes the galaxies inhabit, by comparing with halo masses inferred from observations. Galaxies at high redshift in the models are in lower mass haloes than suggested by observations, and the star formation efficiency in low-mass haloes is higher than observed. We conclude that many of the models require a physical prescription that acts to dissociate the growth of low-mass galaxies from the growth of their dark matter haloes at high redshift.
Cosmic rays (CRs) with ~GeV energies can contribute significantly to the energy and pressure budget in the interstellar, circumgalactic, and intergalactic medium (ISM, CGM, IGM). Recent cosmological simulations have begun to explore these effects, but almost all studies have been restricted to simplified models with constant CR diffusivity and/or streaming speeds. Physical models of CR propagation/scattering via extrinsic turbulence and self-excited waves predict transport coefficients which are complicated functions of local plasma properties. In a companion paper, we consider a wide range of observational constraints to identify proposed physically-motivated cosmic-ray propagation scalings which satisfy both detailed Milky Way (MW) and extra-galactic $gamma$-ray constraints. Here, we compare the effects of these models relative to simpler diffusion+streaming models on galaxy and CGM properties at dwarf through MW mass scales. The physical models predict large local variations in CR diffusivity, with median diffusivity increasing with galacto-centric radii and decreasing with galaxy mass and redshift. These effects lead to a more rapid dropoff of CR energy density in the CGM (compared to simpler models), in turn producing weaker effects of CRs on galaxy star formation rates (SFRs), CGM absorption profiles and galactic outflows. The predictions of the more physical CR models tend to lie in between models which ignore CRs entirely and models which treat CRs with constant diffusivity.
Recent observations have found that many $zsim 6$ quasar fields lack galaxies. This unexpected lack of galaxies may potentially be explained by quasar radiation feedback. In this paper I present a suite of 3D radiative transfer cosmological simulations of quasar fields. I find that quasar radiation suppresses star formation in low mass galaxies, mainly by photo-dissociating their molecular hydrogen. Photo-heating also plays a role, but only after $sim$100 Myr. However, galaxies which already have stellar mass above $10^5 M_odot$ when the quasar turns on will not be suppressed significantly. Quasar radiative feedback suppresses the faint end of the galaxy luminosity function (LF) within $1$ pMpc, but to a far lesser degree than the field-to-field variation of the LF. My study also suggests that by using the number of bright galaxies ($M_{1500}<-16$) around quasars, we can potentially recover the underlying mass overdensity, which allows us to put reliable constraints on quasar environments.
We examine the effect of using different halo finders and merger tree building algorithms on galaxy properties predicted using the GALFORM semi-analytical model run on a high resolution, large volume dark matter simulation. The halo finders/tree builders HBT, ROCKSTAR, SUBFIND and VELOCIRAPTOR differ in their definitions of halo mass, on whether only spatial or phase-space information is used, and in how they distinguish satellite and main haloes; all of these features have some impact on the model galaxies, even after the trees are post-processed and homogenised by GALFORM. The stellar mass function is insensitive to the halo and merger tree finder adopted. However, we find that the number of central and satellite galaxies in GALFORM does depend slightly on the halo finder/tree builder. The number of galaxies without resolved subhaloes depends strongly on the tree builder, with VELOCIRAPTOR, a phase-space finder, showing the largest population of such galaxies. The distributions of stellar masses, cold and hot gas masses, and star formation rates agree well between different halo finders/tree builders. However, because VELOCIRAPTOR has more early progenitor haloes, with these trees GALFORM produces slightly higher star formation rate densities at high redshift, smaller galaxy sizes, and larger stellar masses for the spheroid component. Since in all cases these differences are small we conclude that, when all of the trees are processed so that the main progenitor mass increases monotonically, the predicted GALFORM galaxy populations are stable and consistent for these four halo finders/tree builders.
We present a comparison of 14 galaxy formation models: 12 different semi-analytical models and 2 halo-occupation distribution models for galaxy formation based upon the same cosmological simulation and merger tree information derived from it. The participating codes have proven to be very successful in their own right but they have all been calibrated independently using various observational data sets, stellar models, and merger trees. In this paper we apply them without recalibration and this leads to a wide variety of predictions for the stellar mass function, specific star formation rates, stellar-to- halo mass ratios, and the abundance of orphan galaxies. The scatter is much larger than seen in previous comparison studies primarily because the codes have been used outside of their native environment within which they are well tested and calibrated. The purpose of the `nIFTy comparison of galaxy formation models is to bring together as many different galaxy formation modellers as possible and to investigate a common approach to model calibration. This paper provides a unified description for all participating models and presents the initial, uncalibrated comparison as a baseline for our future studies where we will develop a common calibration framework and address the extent to which that reduces the scatter in the model predictions seen here.