The accuracy of seminumerical reionization models in comparison with radiative transfer simulations

Abstract

We have developed a modular seminumerical code that computes the time and spatially dependent ionization of neutral hydrogen (H I), neutral (He I), and single-ionized helium (He II) in the intergalactic medium (IGM). The model accounts for recombinations and provides different descriptions for the photoionization rate that are used to calculate the residual HI fraction in ionized regions. We compare different seminumerical reionization schemes to a radiative transfer (RT) simulation. We use the RT simulation as a benchmark, and find that the seminumerical approaches produce similar HII and He II morphologies and power spectra of the HI 21 cm signal throughout reionization. As we do not track partial ionization of He II, the extent of the double-ionized helium (He III) regions is consistently smaller. In contrast to previous comparison projects, the ionizing emissivity in our seminumerical scheme is not adjusted to reproduce the redshift evolution of the RT simulation, but directly derived from the RT simulation spectra. Among schemes that identify the ionized regions by the ratio of the number of ionization and absorption events on different spatial smoothing scales, we find those that mark the entire sphere as ionized when the ionization criterion is fulfilled to result in significantly accelerated reionization compared to the RT simulation. Conversely, those that flag only the central cell as ionized yield very similar but slightly delayed redshift evolution of reionization, with up to 20 per cent ionizing photons lost. Despite the overall agreement with the RT simulation, our results suggest that constraining ionizing emissivity-sensitive parameters from seminumerical galaxy formation-reionization models are subject to photon nonconservation.