We incorporate the non-linear clustering of dark matter halos, as modelled by Jose et al. (2016) into the halo model to better understand the clustering of Lyman break galaxies (LBGs) in the redshift range $z=3-5$. We find that, with this change, the predicted LBG clustering increases significantly on quasi-linear scales ($0.1 \leq r\,/\,h^{-1} \,{\rm Mpc} \leq 10$) compared to that in the linear halo bias model. This in turn results in an increase in the clustering of LBGs by an order of magnitude on angular scales $5" \leq \theta \leq 100"$. Remarkably, the predictions of our new model on the whole remove the systematic discrepancy between the linear halo bias predictions and the observations. The correlation length and large scale galaxy bias of LBGs are found to be significantly higher in the non-linear halo bias model than in the linear halo bias model. The resulting two-point correlation function retains an approximate power-law form in contrast with that computed using the linear halo bias theory. We also find that the non-linear clustering of LBGs increases with increasing luminosity and redshift. Our work emphasizes the importance of using non-linear halo bias in order to model the clustering of high-z galaxies to probe the physics of galaxy formation and extract cosmological parameters reliably.

We investigate the spatial clustering of dark matter halos, collapsing from $1-4 \sigma$ fluctuations, in the redshift range $0 - 5$ using N-body simulations. The halo bias of high redshift halos ($z \geq 2$) is found to be strongly non-linear and scale-dependent on quasi-linear scales that are larger than their virial radii ($0.5-10$ Mpc/h). However, at lower redshifts, the scale-dependence of non-linear bias is weaker and and is of the order of a few percent on quasi-linear scales at $z \sim 0$. We find that the redshift evolution of the scale dependent bias of dark matter halos can be expressed as a function of four physical parameters: the peak height of halos, the non-linear matter correlation function at the scale of interest, an effective power law index of the {\it rms} linear density fluctuations and the matter density of the universe at the given redshift. This suggests that the scale-dependence of halo bias is not a universal function of the dark matter power spectrum, which is commonly assumed. We provide a fitting function for the scale dependent halo bias as a function of these four parameters. Our fit reproduces the simulation results to an accuracy of better than 4 % over the redshift range $0\leq z \leq 5$. We also extend our model by expressing the non-linear bias as a function of the linear matter correlation function. It is important to incorporate our results into the clustering models of dark matter halos at any redshift, including those hosting early generations of stars and galaxies before reionization.

We present a galaxy formation model to understand the evolution of stellar mass (M*) - UV luminosity relations, stellar mass functions and specific star formation rate (sSFR) of Lyman Break Galaxies (LBGs) along with their UV luminosity functions in the redshift range 3 < z < 8. Our models assume a physically motivated form for star formation in galaxies and model parameters are calibrated by fitting the observed UV luminosity functions (LFs) of LBGs. We find the fraction of baryons that gets converted into stars remains nearly constant for z < 4 but shows an increase for z < 4. However, the rate of converting baryons into stars does not evolve significantly in the redshift range 3 < z < 8. Our model further successfully explains the M* - UV luminosity (M_AB) correlations of LBGs. While our model predictions of stellar mass functions compare well with the inferred data from observations at the low mass end, we need to invoke the Eddington bias to fit the high mass end. At any given redshift, we find the sSFR to be constant over the stellar mass range 5 \times 10^8 -5 \times 10^9 M_\odot and the redshift evolution of sSFR is well approximated by a form (1+z)^2.4 for 3 < z < 8 which is consistent with observations. Thus we find that dark matter halo build up in the LCDM model is sufficient to explain the evolution of UV LFs of LBGs along with their M* - M_AB relations, the stellar mass functions and the sSFR for 3 < z < 8.

We present a physically motivated semi-analytic model to understand the clustering of high redshift Lyman Alpha Emitters (LAEs). We show that the model parameters constrained by the observed luminosity functions, can be used to predict large scale bias and angular correlation function of LAEs. These predictions are shown to reproduce the observations remarkably well. We find that average masses of dark matter halos hosting LAEs brighter than threshold narrow band magnitude ~ 25 are ~ 10^11 M_\odot. These are smaller than that of typical Lyman Break Galaxies (LBGs) brighter than similar threshold continuum magnitude by a factor ~ 10. This results in a smaller clustering strength of LAEs compared to LBGs. However, using the observed relationship between UV continuum and Lyman-alpha luminosity of LAEs, we show that both LAEs and LBGs belong to the same parent galaxy population with narrow band techniques having greater efficiency in picking up galaxies with low UV luminosity. We also show that the lack of evidence for the presence of the one halo term in the observed LAE angular correlation functions can be attributed to sub-Poisson distribution of LAEs in dark matter halos as a result of their low halo occupations.

We present a physically motivated semi-analytic model to understand the clustering of high redshift LBGs. We show that the model parameters constrained by the observed luminosity function, can be used to predict large scale (\theta > 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 (\theta>80") and small (\theta < 10") angular scales for z=3-5 and several limiting magnitudes. However, our models underpredict the clustering amplitude at intermediate angular scales, where quasi-linear effects are important. The average mass of halos contributing to the observed clustering is found to be 6.2 x 10^{11} M_\odot and the characteristic mass of a parent halo hosting satellite galaxies is 1.2 \times 10^{12} M_\odot 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. We find that approximately 40 % of the halos above a minimum mass M_{min}, can host detectable central galaxies and about 5-10 % of these halos are likely to also host a detectable satellite. The satellites form typically a dynamical timescale prior to the formation of the parent halo. The small angular scale clustering is due to central-satellite pairs and is quite sensitive to changes in the duration of star formation in a halo. The present data favor star formation in a halo lasting typically for a few dynamical time-scales. Our models also reproduce different known trends between parameters related to star formation.

Laboratory experiments measuring neutrino oscillations, indicate small mass differences between different mass eigenstates of neutrinos. The absolute mass scale is however not determined, with at present the strongest upper limits coming from astronomical observations rather than terrestrial experiments. The presence of massive neutrinos suppresses the growth of perturbations below a characteristic mass scale, thereby leading to a decreased abundance of collapsed dark matter halos. Here we show that this effect can significantly alter the predicted luminosity function (LF) of high redshift galaxies. In particular we demonstrate that a stringent constraint on the neutrino mass can be obtained using the well measured galaxy LF and our semi-analytic structure formation models. Combining the constraints from the Wilkinson Microwave Anisotropy Probe 7 year (WMAP7) data with the LF data at z = 4, we get a limit on the sum of the masses of 3 degenerate neutrinos \Sigma m_\nu < 0.52 eV at the 95 % CL. The additional constraints using the prior on Hubble constant strengthens this limit to \Sigma m_\nu < 0.29 eV at the 95 % CL. This neutrino mass limit is a factor of order 4 improvement compared to the constraint based on the WMAP7 data alone, and as stringent as known limits based on other astronomical observations. As different astronomical measurements may suffer from different set of biases, the method presented here provides a complementary probe of \Sigma m_\nu . We suggest that repeating this exercise on well measured luminosity functions over different redshift ranges can provide independent and tighter constraints on \Sigma m_\nu .