The THESAN project: connecting ionized bubble sizes to their local environments during the Epoch of Reionization

TL;DR

This study uses advanced simulations to analyze the growth and characteristics of ionized bubbles during the Epoch of Reionization, revealing their dependence on local environments and implications for upcoming observational probes.

Contribution

The paper introduces a detailed analysis of ionized bubble sizes using the THESAN simulation suite, linking bubble growth to local density and source brightness during reionization.

Findings

Early ionized regions grow slowly, while later regions undergo rapid 'flash ionization'.

Bright sources are found in larger bubbles, consistent with high-redshift observations.

Large bubbles form in high-density regions, but this correlation weakens as reionization progresses.

Abstract

An important characteristic of cosmic hydrogen reionization is the growth of ionized gas bubbles surrounding early luminous objects. Ionized bubble sizes are beginning to be probed using Lyman-$\alpha$ emission from high-redshift galaxies, and will also be probed by upcoming 21-cm maps. We present results from a study of bubble sizes using the state-of-the-art THESAN radiation-hydrodynamics simulation suite, which self-consistently models radiation transport and realistic galaxy formation. We employ the mean-free path method, and track the evolution of the effective ionized bubble size at each point ($R_{\rm eff}$) throughout the Epoch of Reionization. We show there is a slow growth period for regions ionized early, but a rapid "flash ionization" process for regions ionized later as they immediately enter a large, pre-existing bubble. We also find that bright sources are preferentially in larger bubbles, and find consistency with recent observational constraints at $z \gtrsim 9$, but tension with idealized Lyman-$\alpha$ damping-wing models at $z \approx 7$. We find that high overdensity regions have larger characteristic bubble sizes, but the correlation decreases as reionization progresses, likely due to runaway formation of large percolated bubbles. Finally, we compare the redshift at which a region transitions from neutral to ionized ($z_{\rm reion}$) with the time it takes to reach a given bubble size and conclude that $z_{\rm reion}$ is a reasonable local probe of small-scale bubble size statistics ($R_\text{eff} \lesssim 1\,\rm{cMpc}$). However, for larger bubbles, the correspondence between $z_{\rm reion}$ and size statistics weakens due to the time delay between the onset of reionization and the expansion of large bubbles, particularly at high redshifts.