Spatial clustering of high-redshift Lyman-break galaxies

We present a physically motivated semi-analytic model to understand the clustering of high-redshift Lyman-break galaxies (LBGs). We show that the model parameters constrained by the observed luminosity function can be used to predict large-scale (i.e. θ ≥ 80 arcsec) bias and angular correlation function of galaxies. These predictions are shown to reproduce the observations remarkably well. We then adopt these model parameters to calculate the halo occupation distribution (HOD) using the conditional mass function. The halo model using this HOD is shown to provide a reasonably good fit to the observed clustering of LBGs at both large and small (θ < 10 arcsec) angular scales for the whole range of z = 3-5 and limiting magnitudes. However, our models underpredict the clustering amplitude at intermediate angular scales, where quasi-linear effects are important. The average mass of haloes contributing to the observed clustering is found to be 6.2 × 10¹¹ M_⊙ and the characteristic mass of a parent halo hosting satellite galaxies is 1.2 × 10¹² M_⊙ for a limiting absolute magnitude of -20.5 at z = 4. For a given threshold luminosity, these masses decrease with increasing z and at any given z these are found to increase with increasing value of threshold luminosity. Our physical model for the HOD suggests that approximately 40 per cent of the haloes above a minimum mass M_min can host detectable central galaxies and about 5-10 per cent of these haloes are likely to also host a detectable satellite. These satellites form typically a dynamical time-scale prior to the formation of the parent halo. The small angular scale clustering is mainly due to central-satellite pairs rather than a few large clusters. It is quite sensitive to changes in the duration of star formation in a halo and hence could provide a probe of this quantity. The present data favour star formation in a halo lasting typically for a few dynamical time-scales, with 50 per cent of stars formed in a time of T ∼ 300-500 Myr for dark matter haloes that collapse in the redshift range of 5.5-3.5. Our models also reproduce different known trends between parameters related to star formation.