No Arabic abstract
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 use SHARK, a semi-analytic galaxy formation model, to investigate the physical processes involved in dictating the shape, scatter and evolution of the HI-halo mass relation at $0leq z leq 2$. We compare SHARK with HI clustering and spectral stacking of the HI-halo mass relation derived from observations finding excellent agreement with the former and a deficiency of HI in SHARK at $M_{rm vir}approx 10^{12-13} M_{odot}$ in the latter, but otherwise great agreement below and above that mass threshold. In SHARK, we find that the HI mass increases with the halo mass up to a critical mass of $approx 10^{11.8} M_{odot}$; between $sim 10^{11.8}-10^{13}M_{odot}$, the scatter in the relation increases by 0.7 dex and the HI mass decreases with the halo mass on average; at $M_{rm vir} geq 10^{13} M_{odot}$, the HI content continues to increase with halo mass. We find that the critical halo mass of $approx 10^{12} M_{odot}$ is largely set by feedback from Active Galactic Nuclei (AGN), and the exact shape and scatter of the HI-halo mass relation around that mass is extremely sensitive to how AGN feedback is modelled, with other physical processes playing a less significant role. We determine the main secondary parameters responsible for the scatter of the HI-halo mass relation, namely the halo spin parameter at $M_{rm vir}leq 10^{11.8} M_{odot}$, and the fractional contribution from substructure to the total halo mass for $M_{rm vir}geq 10^{13} M_{odot}$. The scatter at $10^{11.8}<M_{rm vir}<10^{13} M_{odot}$ is best described by the black-hole-to-stellar mass ratio of the central galaxy, reflecting the AGN feedback relevance. We present a numerical model to populate dark matter-only simulations with HI at $0leq z leq 2$ based solely on halo parameters that are measurable in such simulations.
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.
We use the 21 cm emission line data from the DINGO-VLA project to study the atomic hydrogen gas H,{textsc i} of the Universe at redshifts $z<0.1$. Results are obtained using a stacking analysis, combining the H,{textsc i} signals from 3622 galaxies extracted from 267 VLA pointings in the G09 field of the Galaxy and Mass Assembly Survey (GAMA). Rather than using a traditional one-dimensional spectral stacking method, a three-dimensional cubelet stacking method is used to enable deconvolution and the accurate recovery of average galaxy fluxes from this high-resolution interferometric dataset. By probing down to galactic scales, this experiment also overcomes confusion corrections that have been necessary to include in previous single dish studies. After stacking and deconvolution, we obtain a $30sigma$ H,{textsc i} mass measurement from the stacked spectrum, indicating an average H,{textsc i} mass of $M_{rm H,{textsc i}}=(1.674pm 0.183)times 10^{9}~{Msun}$. The corresponding cosmic density of neutral atomic hydrogen is $Omega_{rm H,{textsc i}}=(0.377pm 0.042)times 10^{-3}$ at redshift of $z=0.051$. These values are in good agreement with earlier results, implying there is no significant evolution of $Omega_{rm H,{textsc i}}$ at lower redshifts.
A large variance exists in the amplitude of the Stellar Mass - Halo Mass (SMHM) relation for group and cluster-size halos. Using a sample of 254 clusters, we show that the magnitude gap between the brightest central galaxy (BCG) and its second or fourth brightest neighbor accounts for a significant portion of this variance. We find that at fixed halo mass, galaxy clusters with a higher magnitude gap have a higher BCG stellar mass. This relationship is also observed in semi-analytic representations of low-redshift galaxy clusters in simulations. This SMHM-magnitude gap stratification likely results from BCG growth via hierarchical mergers and may link assembly of the halo with the growth of the BCG. Using a Bayesian model, we quantify the importance of the magnitude gap in the SMHM relation using a multiplicative stretch factor, which we find to be significantly non-zero. The inclusion of the magnitude gap in the SMHM relation results in a large reduction in the inferred intrinsic scatter in the BCG stellar mass at fixed halo mass. We discuss the ramifications of this result in the context of galaxy formation models of centrals in group and cluster-sized halos.
The relation between galaxies and dark matter halos is of vital importance for evaluating theoretical predictions of structure formation and galaxy formation physics. We show that the widely used method of abundance matching based on dark matter only simulations fails at the low mass end because two of its underlying assumptions are broken: only a small fraction of low mass (below 10^9.5 solar masses) halos host a visible galaxy, and halos grow at a lower rate due to the effect of baryons. In this regime, reliance on dark matter only simulations for abundance matching is neither accurate nor self-consistent. We find that the reported discrepancy between observational estimates of the halo masses of dwarf galaxies and the values predicted by abundance matching does not point to a failure of LCDM, but simply to a failure to account for baryonic effects. Our results also imply that the Local Group contains only a few hundred observable galaxies in contrast with the thousands of faint dwarfs that abundance matching would suggest. We show how relations derived from abundance matching can be corrected, so that they can be used self-consistently to calibrate models of galaxy formation.