We present predictions for the clustering of galaxies selected by their emission at far infra-red (FIR) and sub-millimetre wavelengths. This includes the first predictions for the effect of clustering biases induced by the coarse angular resolution o
f single-dish telescopes at these wavelengths. We combine a new version of the GALFORM model of galaxy formation with a self-consistent model for calculating the absorption and re-emission of radiation by interstellar dust. Model galaxies selected at $850$ $mu$m reside in dark matter halos of mass $M_{rm halo}sim10^{11.5}-10^{12}$ $h^{-1}$ M$_{odot}$, independent of redshift (for $0.2lesssim zlesssim4$) or flux (for $0.25lesssim S_{850murm m}lesssim4$ mJy). At $zsim2.5$, the brightest galaxies ($S_{850murm m}>4$ mJy) exhibit a correlation length of $r_{0}=5.5_{-0.5}^{+0.3}$ $h^{-1}$ Mpc, consistent with observations. We show that these galaxies have descendants with stellar masses $M_{star}sim10^{11}$ $h^{-1}$ M$_{odot}$ occupying halos spanning a broad range in mass $M_{rm halo}sim10^{12}-10^{14}$ $h^{-1}$ M$_{odot}$. The FIR emissivity at shorter wavelengths ($250$, $350$ and $500$ $mu$m) is also dominated by galaxies in the halo mass range $M_{rm halo}sim10^{11.5}-10^{12}$ $h^{-1}$ M$_{odot}$, again independent of redshift (for $0.5lesssim zlesssim5$). We compare our predictions for the angular power spectrum of cosmic infra-red background anisotropies at these wavelengths with observations, finding agreement to within a factor of $sim2$ over all scales and wavelengths, an improvement over earli
We present predictions for the two-point correlation function of galaxy clustering as a function of stellar mass, computed using two n
Recent observational evidence suggests that the coarse angular resolution ($sim20$ FWHM) of single-dish telescopes at sub-mm wavelengths has biased the observed galaxy number counts by blending together the sub-mm emission from multiple sub-mm galaxi
es (SMGs). We use lightcones computed from an updated implementation of the galform semi-analytic model to generate $50$ mock sub-mm surveys of $0.5$ deg$^2$ at $850$ $mu$m, taking into account the effects of the finite single-dish beam in a more accurate way than has been done previously. We find that blending of SMGs does lead to an enhancement of source extracted number counts at bright fluxes ($S_{mathrm{850}mumathrm{m}}gtrsim1$ mJy). Typically, $sim3{-}6$ galaxies contribute $90%$ of the flux of an $S_{850mumathrm{m}}=5$ mJy source and these blended galaxies are physically unassociated. We find that field-to-field variations are comparable to Poisson fluctuations for our $S_{850mumathrm{m}}>5$ mJy SMG population, which has a median redshift $z_{50}=2.0$, but are greater than Poisson for the $S_{850mumathrm{m}}>1$ mJy population ($z_{50}=2.8$). In a detailed comparison to a recent interferometric survey targeted at single-dish detected sources, we reproduce the difference between single-dish and interferometer number counts and find a median redshift ($z_{50}=2.5$) in excellent agreement with the observed value ($z_{50}=2.5pm 0.2$). We also present predictions for single-dish survey number counts at $450$ and $1100$ $mu$m, which show good agreement with observational data.
The latest observations of molecular gas and the atomic hydrogen content of local and high-redshift galaxies, coupled with how these correlate with star formation activity, have revolutionized our ideas about how to model star formation in a galactic
context. A successful theory of galaxy formation has to explain some key facts: (i) high-redshift galaxies have higher molecular gas fractions and star formation rates than local galaxies, (ii) scaling relations show that the atomic-to-stellar mass ratio decreases with stellar mass in the local Universe, and (iii) the global abundance of atomic hydrogen evolves very weakly with time. We review how modern cosmological simulations of galaxy formation attempt to put these pieces together and highlight how approaches simultaneously solving dark matter and gas physics, and approaches first solving the dark matter N-body problem and then dealing with gas physics using semi-analytic models, differ and complement each other. We review the observable predictions, what we think we have learned so far and what still needs to be done in the simulations to allow robust testing by the new observations expected from telescopes such as ALMA, PdBI, LMT, JVLA, ASKAP, MeerKAT, SKA.
We introduce a method for constructing end-to-end mock galaxy catalogues using a semi-analytical model of galaxy formation, applied to the halo merger trees extracted from a cosmological N-body simulation. The mocks that we construct are lightcone ca
talogues, in which a galaxy is placed according to the epoch at which it first enters the past lightcone of the observer, and incorporate the evolution of galaxy properties with cosmic time. We determine the position between the snapshot outputs at which a galaxy enters the observers lightcone by interpolation. As an application, we consider the effectiveness of the BzK colour selection technique, which was designed to isolate galaxies in the redshift interval 1.4<z<2.5. The mock catalogue is in reasonable agreement with the observed number counts of all BzK galaxies, as well as with the observed counts of the subsample of BzKs that are star-forming galaxies. We predict that over 75 per cent of the model galaxies with K_{AB}<=23, and 1.4<z<2.5, are selected by the BzK technique. Interloper galaxies, outside the intended redshift range, are predicted to dominate bright samples of BzK galaxies (i.e. with K_{AB}<=21). Fainter K-band cuts are necessary to reduce the predicted interloper fraction. We also show that shallow B-band photometry can lead to confusion in classifying BzK galaxies as being star-forming or passively evolving. Overall, we conclude that the BzK colour selection technique is capable of providing a sample of galaxies that is representative of the 1.4<z<2.5 galaxy population.
We use large volume N-body simulations to predict the clustering of dark matter in redshift space in f(R) modified gravity cosmologies. This is the first time that the nonlinear matter and velocity fields have been resolved to such a high level of ac
curacy over a broad range of scales in this class of models. We find significant deviations from the clustering signal in standard gravity, with an enhanced boost in power on large scales and stronger damping on small scales in the f(R) models compared to GR at redshifts z<1. We measure the velocity divergence (P_theta theta) and matter (P_delta delta) power spectra and find a large deviation in the ratios sqrt{P_theta theta/P_delta delta} and P_delta theta/P_deltadelta, between the f(R) models and GR for 0.03<k/(h/Mpc)<0.5. In linear theory these ratios equal the growth rate of structure on large scales. Our results show that the simulated ratios agree with the growth rate for each cosmology (which is scale dependent in the case of modified gravity) only for extremely large scales, k<0.06h/Mpc at z=0. The velocity power spectrum is substantially different in the f(R) models compared to GR, suggesting that this observable is a sensitive probe of modified gravity. We demonstrate how to extract the matter and velocity power spectra from the 2D redshift space power spectrum, P(k,mu), and can recover the nonlinear matter power spectrum to within a few percent for k<0.1h/Mpc. However, the model fails to describe the shape of the 2D power spectrum demonstrating that an improved model is necessary in order to reconstruct the velocity power spectrum accurately. The same model can match the monopole moment to within 3% for GR and 10% for the f(R) cosmology at k<0.2 h/Mpc at z=1. Our results suggest that the extraction of the velocity power spectrum from future galaxy surveys is a promising method to constrain deviations from GR.
We combine the galaxy formation model GALFORM with the Photon Dominated Region code UCL_PDR to study the emission from the rotational transitions of 12CO (CO) in galaxies from z=0 to z=6 in the Lambda CDM framework. GALFORM is used to predict the mol
ecular (H2) and atomic hydrogen (HI) gas contents of galaxies using the pressure-based empirical star formation relation of Blitz & Rosolowsky. From the predicted H2 mass and the conditions in the interstellar medium, we estimate the CO emission in the rotational transitions 1-0 to 10-9 by applying the UCL_PDR model to each galaxy. We find that deviations from the Milky-Way CO-H2 conversion factor come mainly from variations in metallicity, and in the average gas and star formation rate surface densities. In the local universe, the model predicts a CO(1-0) luminosity function (LF), CO-to-total infrared (IR) luminosity ratios for multiple CO lines and a CO spectral line energy distribution (SLED) which are in good agreement with observations of luminous and ultra-luminous IR galaxies. At high redshifts, the predicted CO SLED of the brightest IR galaxies reproduces the shape and normalization of the observed CO SLED. The model predicts little evolution in the CO-to-IR luminosity ratio for different CO transitions, in good agreement with observations up to z~5. We use this new hybrid model to explore the potential of using colour selected samples of high-redshift star-forming galaxies to characterise the evolution of the cold gas mass in galaxies through observations with the Atacama Large Millimeter Array.
For galaxy clustering to provide robust constraints on cosmological parameters and galaxy formation models, it is essential to make reliable estimates of the errors on clustering measurements. We present a new technique, based on a spatial Jackknife
(JK) resampling, which provides an objective way to estimate errors on clustering statistics. Our approach allows us to set the appropriate size for the Jackknife subsamples. The method also provides a means to assess the impact of individual regions on the measured clustering, and thereby to establish whether or not a given galaxy catalogue is dominated by one or several large structures, preventing it to be considered as a fair sample. We apply this methodology to the two- and three-point correlation functions measured from a volume limited sample of M* galaxies drawn from data release seven of the Sloan Digital Sky Survey (SDSS). The frequency of jackknife subsample outliers in the data is shown to be consistent with that seen in large N-body simulations of clustering in the cosmological constant plus cold dark matter cosmology. We also present a comparison of the three-point correlation function in SDSS and 2dFGRS using this approach and find consistent measurements between the two samples.
We study the evolution of the cold gas content of galaxies by splitting the interstellar medium into its atomic and molecular hydrogen components, using the galaxy formation model GALFORM in the LCDM framework. We calculate the molecular-to-atomic hy
drogen mass ratio, H2/HI, in each galaxy using two different approaches; the pressure-based empirical relation of Blitz & Rosolowsky and the theoretical model of Krumholz, McKeee & Tumlinson, and apply them to consistently calculate the star formation rates of galaxies. We find that the model based on the Blitz & Rosolowsky law predicts an HI mass function, CO(1-0) luminosity function, correlations between the H2/HI ratio and stellar and cold gas mass, and infrared-CO luminosity relation in good agreement with local and high redshift observations. The HI mass function evolves weakly with redshift, with the number density of high mass galaxies decreasing with increasing redshift. In the case of the H2 mass function, the number density of massive galaxies increases strongly from z=0 to z=2, followed by weak evolution up to z=4. We also find that the H2/HI ratio of galaxies is strongly dependent on stellar and cold gas mass, and also on redshift. The slopes of the correlations between H2/HI and stellar and cold gas mass hardly evolve, but the normalisation increases by up to two orders of magnitude from z=0-8. The strong evolution in the H2 mass function and the H2/HI ratio is primarily due to the evolution in the sizes of galaxies and secondarily, in the gas fractions. The predicted cosmic density evolution of HI agrees with the observed evolution inferred from DLAs, and is dominated by low/intermediate mass halos. We find that previous theoretical studies have largely overestimated the redshift evolution of the global H2/HI ratio due to limited resolution. We predict a maximum of rho_H2/rho_HI~1.2 at z~3.5.
We investigate the consequences of applying different star formation laws in the galaxy formation model GALFORM. Three broad star formation laws are implemented: the empirical relations of Kennicutt and Schmidt and Blitz & Rosolowsky and the theoreti
cal model of Krumholz, McKee & Tumlinson. These laws have no free parameters once calibrated against observations of the star formation rate (SFR) and gas surface density in nearby galaxies. We start from published models, and investigate which observables are sensitive to a change in the star formation law, without altering any other model parameters. We show that changing the star formation law (i) does not significantly affect either the star formation history of the universe or the galaxy luminosity functions in the optical and near-IR, due to an effective balance between the quiescent and burst star formation modes; (ii) greatly affects the cold gas contents of galaxies; (iii) changes the location of galaxies in the SFR versus stellar mass plane, so that a second sequence of passive galaxies arises, in addition to the known active sequence. We show that this plane can be used to discriminate between the star formation laws.
