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 α) 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 α 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 α visibility during reionization. We show that UV luminosity–halo mass scatter reduces Ly α 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{ H}\,{\small I}}$, at z ∼ 7 using a sample of Lyman-break galaxies in legacy fields. We infer $\overline{x}_{\mathrm{ H}\,{\small I}} = 0.55_{-0.13}^{+0.11}$ with scatter, compared to $\overline{x}_{\mathrm{ H}\,{\small I}} = 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 most luminous galaxies at high-redshift are generally considered to be hosted in massive darkmatter 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 ∼ 8 − 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 Σ 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 Σ, 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 Σ 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 Σ min (z) ∼ 0.15 − 0.3, with Σ min increasing with redshift as a consequence of shorter halo assembly periods at higher redshifts. Current observations (m AB ∼ 27) 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 Σ ≥ Σ min . Deeper future observations reaching m AB ∼ 28.2 − 29 would have the opportunity to clearly quantify the clustering strength, and hence to constrain Σ, investigating the physical processes that drive star formation in the early Universe.

We report on a Hubble Space Telescope search for rest-frame ultraviolet emission from the host galaxies of five far-infrared-luminous z ≃ 6 quasars and the z = 5.85 hot-dust-free quasar SDSS J0005–0006. We perform 2D surface brightness modeling for each quasar using a Markov Chain Monte Carlo estimator, to simultaneously fit and subtract the quasar point source in order to constrain the underlying host galaxy emission. We measure upper limits for the quasar host galaxies of m J  > 22.7 mag and m H  > 22.4 mag, corresponding to stellar masses of M * < 2 × 1011 M ⊙. These stellar mass limits are consistent with the local M BH − M * relation. Our flux limits are consistent with those predicted for the UV stellar populations of z ≃ 6 host galaxies, but likely in the presence of significant dust ( mag). We also detect a total of up to nine potential z ≃ 6 quasar companion galaxies surrounding five of the six quasars, separated from the quasars by 1.″4–3.″2, or 8.4–19.4 kpc, which may be interacting with the quasar hosts. These nearby companion galaxies have UV absolute magnitudes of −22.1 to −19.9 mag and UV spectral slopes β of −2.0 to −0.2, consistent with luminous star-forming galaxies at z ≃ 6. These results suggest that the quasars are in dense environments typical of luminous z ≃ 6 galaxies. However, we cannot rule out the possibility that some of these companions are foreground interlopers. Infrared observations with the James Webb Space Telescope will be needed to detect the z ≃ 6 quasar host galaxies and better constrain their stellar mass and dust content.