We explore a potential new probe of reionization: the cross-correlation of high-redshift galaxies with maps of the thermal Sunyaev-Zel'dovich (tSZ) effect. We consider two types of high redshift galaxies: Lyman break galaxies (LBGs) and Lyman-$\alpha$ emitters (LAEs). LBGs and LAEs will be detected in large numbers at high redshift ($z \approx$ 4 to 7) by ongoing and future surveys. We consider a LBG sample like that expected from the The Rubin Observatory Legacy Survey of Space and Time (LSST), and a selection of LAEs modelled after the Subaru SILVERRUSH program, but covering a much larger sky fraction. The tSZ effect is sensitive to a line-of-sight integral of the ionized gas pressure, and can be measured across large patches of the sky using multi-frequency CMB surveys. We consider forecast tSZ maps from CMB Stage 4 and more futuristic observations. Using a suite of hydrodynamical simulations, we show that the high-redshift galaxies are correlated with the tSZ signal from reionization, with cross-power amplitude of order $10^{-15}$ at $\ell \sim 1000$. The cross-spectra between LBGs/LAEs with tSZ maps contain information about the reionization history of the Universe, such as the distribution of bubble sizes, and could in principle be used to directly measure the timing of reionization. The amplitude of the signal is small, however, and its detectability is hindered by low-redshift contributions to tSZ maps and by instrumental noise. If the the low-redshift contribution to the observed tSZ signal is suppressed by masking of massive halos, a combination of overlapping futuristic CMB and galaxy surveys could potentially probe this signal.

We model the 21-cm signal and LAE population evolution during the epoch of reionization in order to predict the 21cm-LAE cross-power spectrum. We employ high-dynamic-range simulations of the IGM to create models that are consistent with constraints from the CMB, Lyman-$\alpha$ forest and LAE population statistics. Using these models we consider the evolution of the cross-power spectrum for a selection of realistic reionization histories and predict the sensitivity of current and upcoming surveys to measuring this signal. We find that the imprint of a delayed-end to reionization can be observed by future surveys, and that strong constraints can be placed on the progression of reionization as late as $z=5.7$ using a Subaru-SKA survey. We make predictions for the signal-to-noise ratios achievable by combinations of Subaru/PFS with the MWA, LOFAR, HERA and SKA interferometers for an integration time of 1000 hours. We find that a Subaru-SKA survey could measure the cross-power spectrum for a late reionization at $z=6.6$ with a total signal-to-noise greater than 5, making it possible to constrain both the timing and bubble size at the end of reionization. Furthermore, we find that expanding the current Subaru/PFS survey area and depth by a factor of three would double the total signal-to-noise.

A long (110 cMpc/$h$) and deep absorption trough in the Ly$\alpha$ forest has been observed extending down to redshift 5.5 in the spectrum of ULAS J0148+0600. Although no Ly$\alpha$ transmission is detected, Ly$\beta$ spikes are present which has led to claims that the gas along this trough must be ionized. Using high resolution cosmological radiative transfer simulations in large volumes, we show that in a scenario where reionization ends late ($z \sim 5.2$), our simulations can reproduce troughs as long as observed. In this model, we find that the troughs are caused by islands of neutral hydrogen. Small ionized holes within the neutral islands allow for the transmission of Ly$\beta$. We have also modelled the Ly$\alpha$ emitter population around the simulated troughs, and show that there is a deficit of Ly$\alpha$ emitters close to the trough as is observed.

We model the high redshift (z > 5) Lyman-$\alpha$ emitting (LAE) galaxy population using the empirical rest-frame equivalent width distribution. We calibrate to the observed luminosity function and angular correlation function at z = 5.7 as measured by the SILVERRUSH survey. This allows us to populate the high-dynamic-range Sherwood simulation suite with LAEs, and to calculate the transmission of their Ly $\alpha$ emission through the inter-galactic medium (IGM). We use this simulated population to explore the effect of the IGM on high-redshift observations of LAEs, and make predictions for the narrowband filter redshifts at z = 6.6, 7.0 and 7.3. Comparing our model with existing observations, we find a late reionization is suggested, consistent with the recent low optical depth derived from the cosmic microwave background (CMB) by the Planck collaboration and the opacity fluctuations in the Ly $\alpha$ forest. We also explore the role of the circum-galactic medium (CGM) and the larger volume of gas which is infalling onto the host halo versus the IGM in attenuating the Ly $\alpha$ signal, finding that a significant fraction of the attenuation is due to the CGM and infalling gas, which increases towards the end of reionization, albeit with a large scatter across the mock LAE population.

We model the transmission of the Lyman-alpha line through the circum- and intergalactic media around dark matter haloes expected to host Lyman-alpha emitters (LAEs) at z > 5.7, using the high-dynamic-range Sherwood simulations. We find very different CGM environments around more massive haloes (~10^11 M_sun) compared to less massive haloes (~10^9 M_sun) at these redshifts, which can contribute to a different evolution of the Lyman-alpha transmission from LAEs within these haloes. Additionally we confirm that part of the differential evolution could result from bright LAEs being more likely to reside in larger ionized regions. We conclude that a combination of the CGM environment and the IGM ionization structure is likely to be responsible for the differential evolution of the bright and faint ends of the LAE luminosity function at z > 6. More generally, we confirm the suggestion that the self-shielded neutral gas in the outskirts of the host halo can strongly attenuate the Lyman-alpha emission from high redshift galaxies. We find that this has a stronger effect on the more massive haloes hosting brighter LAEs. The faint-end of the LAE luminosity function is thus a more reliable probe of the average ionization state of the IGM. Comparing our model for LAEs with a range of observational data we find that the favoured reionization histories are our previously advocated `Late' and `Very Late' reionization histories, in which reionization finishes rather rapidly at around z ~ 6.