The reionization of hydrogen is closely linked to the first structures in the universe, so understanding the timeline of reionization promises to shed light on the nature of these early objects. In particular, transmission of Lyman alpha (Ly$\alpha$) from galaxies through the intergalactic medium (IGM) is sensitive to neutral hydrogen in the IGM, so can be used to probe the reionization timeline. In this work, we implement an improved model of the galaxy UV luminosity to dark matter halo mass relation to infer the volume-averaged fraction of neutral hydrogen in the IGM from Ly$\alpha$ observations. Many models assume that UV-bright galaxies are hosted by massive dark matter haloes in overdense regions of the IGM, so reside in relatively large ionized regions. However, observations and N-body simulations indicate that scatter in the UV luminosity-halo mass relation is expected. Here, we model the scatter (though we assume the IGM topology is unaffected) and assess the impact on Ly$\alpha$ visibility during reionization. We show that UV luminosity-halo mass scatter reduces Ly$\alpha$ visibility compared to models without scatter, and that this is most significant for UV-bright galaxies. We then use our model with scatter to infer the neutral fraction, $\overline{x}_{\mathrm{HI}}$, at $z \sim 7$ using a sample of Lyman-break galaxies in legacy fields. We infer $\overline{x}_{\mathrm{HI}} = 0.55_{-0.13}^{+0.11}$ with scatter, compared to $\overline{x}_{\mathrm{HI}} = 0.59_{-0.14}^{+0.12}$ without scatter, a very slight decrease and consistent within the uncertainties. Finally, we place our results in the context of other constraints on the reionization timeline and discuss implications for future high-redshift galaxy studies.

The Quasar Luminosity Function (QLF) offers insight into the early co-evolution of black holes and galaxies. It has been characterized observationally up to redshift $z\sim6$ with clear evidence of a double power-law shape, in contrast to the Schechter-like form of the underlying dark-matter halo mass function. We investigate a physical origin for the difference in these distributions by considering the impact of stochasticity induced by the processes that determine the quasar luminosity for a given host halo and redshift. We employ a conditional luminosity function and construct the relation between median quasar magnitude versus halo mass $M_{UV,\rm{c}}(M_{\rm{h}})$ with log-normal in luminosity scatter $\Sigma$, and duty-cycle $\epsilon_{\rm{DC}}$, and focus on high redshift $z\gtrsim4$. We show that, in order to reproduce the observed QLF, the $\Sigma=0$ abundance matching requires all of the brightest quasars to be hosted in the rarest most massive dark-matter halos (with an increasing $M_{UV,\rm{c}}/M_{\rm{h}}$ in halo mass). Conversely, for $\Sigma>0$ the brightest quasars can be over-luminous outliers hosted in relatively common dark-matter halos. In this case, the median quasar magnitude versus halo mass relation, $M_{UV,\rm{c}}$, flattens at the high-end, as expected in self-regulated growth due to feedback. We sample the parameter space of $\Sigma$ and $\epsilon_{\rm{DC}}$ and show that $M_{UV,\rm{c}}$ flattens above $M_{\rm{h}}\sim 10^{12}M_{\odot}$ for $\epsilon_{\rm{DC}}<10^{-2}$. Models with $\epsilon_{\rm{DC}}\sim1$ instead require a high mass threshold close to $M_{\rm{h}}\gtrsim10^{13}M_{\odot}$. We investigate the impact of $\epsilon_{\rm{DC}}$ and $\Sigma$ on measurements of clustering and find there is no luminosity dependence on clustering for $\Sigma>0.3$, consistent with recent observations from Subaru HSC.

The Ultraviolet Luminosity Function (UVLF) is a key observable for understanding galaxy formation from cosmic dawn. There has been considerable debate on whether Schechter-like LFs (characterized by an exponential drop-off at the bright end) that well describe the LF in our local Universe are also a sufficient description of the LF at high redshifts ($z>6$). We model the UVLF over cosmic history with a semi-empirical framework and include a log-normal scatter, $\Sigma$, in galaxy luminosities with a conditional luminosity function approach. We show that stochasticity induces a flattening or a feedback scale in the median galaxy luminosity versus halo mass relation, $L_{c}(M_{h})$ to account for the increase of bright objects placed in lower mass halos. We observe a natural broadening in the bright-end exponential segment of the UVLF for $z>6$ if processes that regulate star-formation acts on the same mass scale as at $z\sim5$, where the degree of broadening is enhanced for larger $\Sigma$. Alternatively, if the bright-end feedback is triggered at a near-constant luminosity threshold, the feedback threshold occurs at progressively lower halo masses with increasing redshift, due to galaxies being more luminous on average at a fixed halo mass from rapid halo assembly. Such feedback results in a LF shape with a bright-end closer to that of a Schechter function. We include predictions for the $z>8$ UVLFs from future all-sky surveys such as WFIRST which has the potential to both quantify the scatter and type of feedback, and provide insight behind the mechanisms that drive star formation in the early Universe.

The most luminous galaxies at high-redshift are generally considered to be hosted in massive dark-matter halos of comparable number density, hence residing at the center of overdensities/protoclusters. We assess the validity of this assumption by investigating the clustering around the brightest galaxies populating the cosmic web at redshift $z\sim8-9$ through a combination of semi-analytic modeling and Monte Carlo simulations of mock Hubble Space Telescope WFC3 observations. The innovative aspect of our approach is the inclusion of a log-normal scatter parameter $\Sigma$ in the galaxy luminosity versus halo mass relation, extending to high-$z$ the conditional luminosity function framework extensively used at low redshift. Our analysis shows that the larger the value of $\Sigma$, the less likely that the brightest source in a given volume is hosted in the most massive halo, and hence the weaker the overdensity of neighbors. We derive a minimum value of $\Sigma$ as a function of redshift by considering stochasticity in the halo assembly times, which affects galaxy ages and star formation rates in our modeling. We show that $\Sigma_{min}(z)\sim0.15-0.3$, with $\Sigma_{min}$ increasing with redshift as a consequence of shorter halo assembly periods at higher redshifts. Current observations ($m_{AB}\sim27$) of the environment of spectroscopically confirmed bright sources at $z>7.5$ do not show strong evidence of clustering and are consistent with our modeling predictions for $\Sigma\geq\Sigma_{min}$. Deeper future observations reaching $m_{AB}\sim28.2-29$ would have the opportunity to clearly quantify the clustering strength, and hence to constrain $\Sigma$, investigating the physical processes that drive star formation in the early Universe.