No Arabic abstract
CO measurements of z~1-4 galaxies have found that their baryonic gas fractions are significantly higher than galaxies at z=0, with values ranging from 20-80 %. Here, we suggest that the gas fractions inferred from observations of star-forming galaxies at high-z are overestimated, owing to the adoption of locally-calibrated CO-H2 conversion factors (Xco). Evidence from both observations and numerical models suggest that Xco varies smoothly with the physical properties of galaxies, and that Xco can be parameterised simply as a function of both gas phase metallicity and observed CO surface brightness. When applying this functional form, we find fgas ~10-40 % in galaxies with M*=10^10-10^12 Msun at high-z. Moreover, the scatter in the observed fgas-M* relation is lowered by a factor of two. The lower inferred gas fractions arise physically because the interstellar media of high-z galaxies have higher velocity dispersions and gas temperatures than their local counterparts, which results in an Xco that is lower than the z=0 value for both quiescent discs and starbursts. We further compare these gas fractions to those predicted by cosmological galaxy formation models. We show that while the canonically inferred gas fractions from observations are a factor of 2-3 larger at a given stellar mass than predicted by models, our rederived Xco values for z=1-4 galaxies results in revised gas fractions that agree significantly better with the simulations.
We use cosmological hydrodynamic simulations to investigate how inflows, star formation, and outflows govern the the gaseous and metal content of galaxies. In our simulations, galaxy metallicities are established by a balance between inflows and outflows as governed by the mass outflow rate, implying that the mass-metallicity relation reflects how the outflow rate varies with stellar mass (M*). Gas content is set by a competition between inflow into and gas consumption within the ISM, the latter being governed by the SF law, while the former is impacted by both wind recycling and preventive feedback. Stochastic variations in the inflow rate move galaxies off the equilibrium M*-Z and Z*-fgas relations in a manner correlated with star formation rate, and the scatter is set by the timescale to re-equilibrate. The evolution of both relations from z=3-0 is slow, as individual galaxies tend to evolve mostly along the relations. Gas fractions at a given M* slowly decrease with time because the cosmic inflow rate diminishes faster than the consumption rate, while metallicities slowly increase as infalling gas becomes more enriched. Observations from z~3-0 are better matched by simulations employing momentum-driven wind scalings rather than constant wind speeds, but all models predict too low gas fractions at low masses and too high metallicities at high M*. All our models reproduce observed second-parameter trends of the mass-metallicity relation with star formation rate and environment, indicating that these are a consequence of equilibrium and not feedback. Overall, the analytical framework of our equilibrium scenario broadly captures the relevant physics establishing the galaxy gas and metal content in simulations, which suggests that the cycle of baryonic inflows and outflows centrally governs the cosmic evolution of these properties in typical star-forming galaxies.
We review recent progress in the description of the formation and evolution of galaxy clusters in a cosmological context by using numerical simulations. We focus our presentation on the comparison between simulated and observed X-ray properties, while we will also discuss numerical predictions on properties of the galaxy population in clusters. Many of the salient observed properties of clusters, such as X-ray scaling relations, radial profiles of entropy and density of the intracluster gas, and radial distribution of galaxies are reproduced quite well. In particular, the outer regions of cluster at radii beyond about 10 per cent of the virial radius are quite regular and exhibit scaling with mass remarkably close to that expected in the simplest case in which only the action of gravity determines the evolution of the intra-cluster gas. However, simulations generally fail at reproducing the observed cool-core structure of clusters: simulated clusters generally exhibit a significant excess of gas cooling in their central regions, which causes an overestimate of the star formation and incorrect temperature and entropy profiles. The total baryon fraction in clusters is below the mean universal value, by an amount which depends on the cluster-centric distance and the physics included in the simulations, with interesting tensions between observed stellar and gas fractions in clusters and predictions of simulations. Besides their important implications for the cosmological application of clusters, these puzzles also point towards the important role played by additional physical processes, beyond those already included in the simulations. We review the role played by these processes, along with the difficulty for their implementation, and discuss the outlook for the future progress in numerical modeling of clusters.
We present the McMaster Unbiased Galaxy Simulations (MUGS), the first 9 galaxies of an unbiased selection ranging in total mass from 5$times10^{11}$ M$_odot$ to 2$times10^{12}$ M$_odot$ simulated using n-body smoothed particle hydrodynamics (SPH) at high resolution. The simulations include a treatment of low temperature metal cooling, UV background radiation, star formation, and physically motivated stellar feedback. Mock images of the simulations show that the simulations lie within the observed range of relations such as that between color and magnitude and that between brightness and circular velocity (Tully-Fisher). The greatest discrepancy between the simulated galaxies and observed galaxies is the high concentration of material at the center of the galaxies as represented by the centrally peaked rotation curves and the high bulge-to-total ratios of the simulations determined both kinematically and photometrically. This central concentration represents the excess of low angular momentum material that long has plagued morphological studies of simulated galaxies and suggests that higher resolutions and a more accurate description of feedback will be required to simulate more realistic galaxies. Even with the excess central mass concentrations, the simulations suggest the important role merger history and halo spin play in the formation of disks.
We study dwarf galaxy formation at high redshift ($zge5$) using a suite of high- resolution, cosmological hydrodynamic simulations and a semi-analytic model (SAM). We focus on gas accretion, cooling and star formation in this work by isolating the relevant process from reionization and supernova feedback, which will be further discussed in a companion paper. We apply the SAM to halo merger trees constructed from a collisionless N-body simulation sharing identical initial conditions to the hydrodynamic suite, and calibrate the free parameters against the stellar mass function predicted by the hydrodynamic simulations at z = 5. By making comparisons of the star formation history and gas components calculated by the two modelling techniques, we find that semi-analytic prescriptions that are commonly adopted in the literature of low-redshift galaxy formation do not accurately represent dwarf galaxy properties in the hydrodynamic simulation at earlier times. We propose 3 modifications to SAMs that will provide more accurate high-redshift simulations. These include 1) the halo mass and baryon fraction which are overestimated by collisionless N-body simulations; 2) the star formation efficiency which follows a different cosmic evolutionary path from the hydrodynamic simulation; and 3) the cooling rate which is not well defined for dwarf galaxies at high redshift. Accurate semi-analytic modelling of dwarf galaxy formation informed by detailed hydrodynamical modelling will facilitate reliable semi-analytic predictions over the large volumes needed for the study of reionization.
We present results from a comprehensive imaging survey of 70 radio galaxies at redshifts 1<z<5.2 using all three cameras onboard the Spitzer Space Telescope. The resulting spectral energy distributions unambiguously show a stellar population in 46 sources and hot dust emission associated with the active nucleus in 59. Using a new restframe S_3um/S_1.6um versus S_um/S_3um criterion, we identify 42 sources where the restframe 1.6um emission from the stellar population can be measured. For these radio galaxies, the median stellar mass is high, 2x10^11 M_sun, and remarkably constant within the range 1<z<3. At z>3, there is tentative evidence for a factor of two decrease in stellar mass. This suggests that radio galaxies have assembled the bulk of their stellar mass by z~3, but confirmation by more detailed decomposition of stellar and AGN emission is needed. The restframe 500 MHz radio luminosities are only marginally correlated with stellar mass but are strongly correlated with the restframe 5um hot dust luminosity. This suggests that the radio galaxies have a large range of Eddington ratios. We also present new Very Large Array 4.86 and 8.46 GHz imaging of 14 radio galaxies and find that radio core dominance --- an indicator of jet orientation --- is strongly correlated with hot dust luminosity. While all of our targets were selected as narrow-lined, type 2 AGNs, this result can be understood in the context of orientation-dependent models if there is a continuous distribution of orientations from obscured type 2 to unobscured type 1 AGNs rather than a clear dichotomy. Finally, four radio galaxies have nearby (<6) companions whose mid-IR colors are suggestive of their being AGNs. This may indicate an association between radio galaxy activity and major mergers.