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-α emitters (LAEs). LBGs and LAEs will be detected in large numbers at high redshift (z ≈ 4 – 7) by ongoing and future surveys. We consider a future LBG sample from 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 is measured across large patches of 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 LBGs and LAEs are correlated with the tSZ signal from reionization. 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 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 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 probe this signal
