No Arabic abstract
We use spectral stacking to measure the contribution of galaxies of different masses and in different hierarchies to the cosmic atomic hydrogen (HI) mass density in the local Universe. Our sample includes 1793 galaxies at $z < 0.11$ observed with the Westerbork Synthesis Radio Telescope, for which Sloan Digital Sky Survey spectroscopy and hierarchy information are also available. We find a cosmic HI mass density of $Omega_{rm HI} = (3.99 pm 0.54)times 10^{-4} h_{70}^{-1}$ at $langle zrangle = 0.065$. For the central and satellite galaxies, we obtain $Omega_{rm HI}$ of $(3.51 pm 0.49)times 10^{-4} h_{70}^{-1}$ and $(0.90 pm 0.16)times 10^{-4} h_{70}^{-1}$, respectively. We show that galaxies above and below stellar masses of $sim$10$^{9.3}$ M$_{odot}$ contribute in roughly equal measure to the global value of $Omega_{rm HI}$. While consistent with estimates based on targeted HI surveys, our results are in tension with previous theoretical work. We show that these differences are, at least partly, due to the empirical recipe used to set the partition between atomic and molecular hydrogen in semi-analytical models. Moreover, comparing our measurements with the cosmological semi-analytic models of galaxy formation {sc Shark} and GALFORM reveals gradual stripping of gas via ram pressure works better to fully reproduce the properties of satellite galaxies in our sample, than strangulation. Our findings highlight the power of this approach in constraining theoretical models, and confirm the non-negligible contribution of massive galaxies to the HI mass budget of the local Universe.
Measuring the HI-halo mass scaling relation (HIHM) is fundamental to understanding the role of HI in galaxy formation and its connection to structure formation. While direct measurements of the HI mass in haloes are possible using HI-spectral stacking, the reported shape of the relation depends on the techniques used to measure it (e.g. monotonically increasing with mass versus flat, mass-independent). Using a simulated HI and optical survey produced with the SHARK semi-analytic galaxy formation model, we investigate how well different observational techniques can recover the intrinsic, theoretically predicted, HIHM relation. We run a galaxy group finder and mimic the HI stacking procedure adopted by different surveys and find we can reproduce their observationally derived HIHM relation. However, none of the adopted techniques recover the underlying HIHM relation predicted by the simulation. We find that systematic effects in halo mass estimates of galaxy groups modify the inferred shape of the HIHM relation from the intrinsic one in the simulation, while contamination by interloping galaxies, not associated with the groups, contribute to the inferred HI mass of a halo mass bin, when using large velocity windows for stacking. The effect of contamination is maximal at Mvir~10^(12-12.5)Msol. Stacking methods based on summing the HI emission spectra to infer the mean HI mass of galaxies of different properties belonging to a group suffer minimal contamination but are strongly limited by the use of optical counterparts, which miss the contribution of dwarf galaxies. Deep spectroscopic surveys will provide significant improvements by going deeper while maintaining high spectroscopic completeness; for example, the WAVES survey will recover ~52% of the total HI mass of the groups with Mvir~10^(14)Msol compared to ~21% in GAMA.
We explore the galaxy formation physics governing the low mass end of the HI mass function in the local Universe. Specifically, we predict the effects on the HI mass function of varying i) the strength of photoionisation feedback and the redshift of the end of the epoch of reionization, ii) the cosmology, iii) the supernovae feedback prescription, and iv) the efficiency of star formation. We find that the shape of the low-mass end of the HI mass function is most affected by the critical halo mass below which galaxy formation is suppressed by photoionisation heating of the intergalactic medium. We model the redshift dependence of this critical dark matter halo mass by requiring a match to the low-mass end of the HI mass function. The best fitting critical dark matter halo mass decreases as redshift increases in this model, corresponding to a circular velocity of $sim 50 , {rm km ,s}^{-1}$ at $z=0$, $sim 30 , {rm km, s}^{-1}$ at $z sim 1$ and $sim 12 , {rm km , s}^{-1}$ at $z=6$. We find that an evolving critical halo mass is required to explain both the shape and abundance of galaxies in the HI mass function below $M_{rm HI} sim 10^{8} h^{-2} {rm M_{odot}}$. The model makes specific predictions for the clustering strength of HI-selected galaxies with HI masses > $10^{6} h^{-2} {rm M_{odot}}$ and $> 10^{7} h^{-2} {rm M_{odot}}$ and for the relation between the HI and stellar mass contents of galaxies which will be testable with upcoming surveys with the Square Kilometre Array and its pathfinders. We conclude that measurements of the HI mass function at $z ge 0$ will lead to an improvement in our understanding of the net effect of photoionisation feedback on galaxy formation and evolution.
We present the dust mass function (DMF) of 15,750 galaxies with redshift $z< 0.1$, drawn from the overlapping area of the GAMA and {it H-}ATLAS surveys. The DMF is derived using the density corrected $V_{rm max}$ method, where we estimate $V_{rm max}$ using: (i) the normal photometric selection limit ($pV_{rm max}$) and (ii) a bivariate brightness distribution (BBD) technique, which accounts for two selection effects. We fit the data with a Schechter function, and find $M^{*}=(4.65pm0.18)times 10^{7},h^2_{70}, M_{odot}$, $alpha=(1.22pm 0.01)$, $phi^{*}=(6.26pm 0.28)times 10^{-3},h^3_{70},rm Mpc^{-3},dex^{-1}$. The resulting dust mass density parameter integrated down to $10^4,M_{odot}$ is $Omega_{rm d}=(1.11 pm0.02)times 10^{-6}$ which implies the mass fraction of baryons in dust is $f_{m_b}=(2.40pm0.04)times 10^{-5}$; cosmic variance adds an extra 7-17,per,cent uncertainty to the quoted statistical errors. Our measurements have fewer galaxies with high dust mass than predicted by semi-analytic models. This is because the models include too much dust in high stellar mass galaxies. Conversely, our measurements find more galaxies with high dust mass than predicted by hydrodynamical cosmological simulations. This is likely to be from the long timescales for grain growth assumed in the models. We calculate DMFs split by galaxy type and find dust mass densities of $Omega_{rm d}=(0.88pm0.03)times 10^{-6}$ and $Omega_{rm d}=(0.060pm0.005)times 10^{-6}$ for late-types and early-types respectively. Comparing to the equivalent galaxy stellar mass functions (GSMF) we find that the DMF for late-types is well matched by the GMSF scaled by $(8.07pm0.35) times 10^{-4}$.
We study the galaxy cosmological mass function (GCMF) in a semi-empirical relativistic approach using observational data provided by galaxy redshift surveys. Starting from the theory of Ribeiro & Stoeger (2003, arXiv:astro-ph/0304094) between the mass-to-light ratio, the selection function obtained from the luminosity function (LF) data and the luminosity density, the average luminosity $L$ and the average galactic mass $mathcal{M}_g$ are computed in terms of the redshift. $mathcal{M}_g$ is also alternatively estimated by a method that uses the galaxy stellar mass function (GSMF). Comparison of these two forms of deriving the average galactic mass allows us to infer a possible bias introduced by the selection criteria of the survey. We used the FORS Deep Field galaxy survey sample of 5558 galaxies in the redshift range $0.5 < z < 5.0$ and its LF Schechter parameters in the B-band, as well as this samples stellar mass-to-light ratio and its GSMF data. Assuming ${mathcal{M}_{g_0}} approx 10^{11} mathcal{M}_odot$ as the local value of the average galactic mass, the LF approach results in $L_{B} propto (1+z)^{(2.40 pm 0.03)}$ and $mathcal{M}_g propto (1+z)^{(1.1pm0.2)}$. However, using the GSMF results produces $mathcal{M}_g propto (1+z)^{(-0.58 pm 0.22)}$. We chose the latter result as it is less biased. We then obtained the theoretical quantities of interest, such as the differential number counts, to calculate the GCMF, which can be fitted by a Schechter function. The derived GCMF follows theoretical predictions in which the less massive objects form first, being followed later by more massive ones. In the range $0.5 < z < 2.0$ the GCMF has a strong variation that can be interpreted as a higher rate of galaxy mergers or as a strong evolution in the star formation history of these galaxies.
In this paper, we present a new derivation of the shape and evolution of the integrated galaxy-wide initial mass function (IGIMF), incorporating explicitly the effects of cosmic rays (CRs) as regulators of the chemical and thermal state of the gas in the dense cores of molecular clouds. We predict the shape of the IGIMF as a function of star formation rate (SFR) and CR density, and show that it can be significantly different with respect to local estimates. In particular, we focus on the physical conditions corresponding to IGIMF shapes that are simultaneously shallower at high-mass end and steeper at the low-mass end than a Kroupa IMF. These solutions can explain both the levels of $alpha$-enrichment and the excess of low-mass stars as a function of stellar mass, observed for local spheroidal galaxies. As a preliminary test of our scenario, we use idealized star formation histories to estimate the mean IMF shape for galaxies of different $z=0$ stellar mass. We show that the fraction of low-mass stars as a function of galaxy stellar mass predicted by these mean IMFs agrees with the values derived from high-resolution spectroscopic surveys.