Simulating the 21 cm signal from reionization including non-linear ionizations and inhomogeneous recombinations

We explore the impact of incorporating physically motivated ionization and recombination rates on the history and topology of cosmic reionization and the resulting 21 cm power spectrum, by incorporating inputs from small-volume hydrodynamic simulations into our semi-numerical code, simfast21, that evolves reionization on large scales. We employ radiative hydrodynamic simulations to parametrize the ionization rate R-ion and recombination rate R-rec as functions of halo mass, overdensity and redshift. We find that R-ion scales superlinearly with halo mass (R-ion proportional to M-h(1.41)), in contrast to previous assumptions. Implementing these scalings into simfast21, we tune our one free parameter, the escape fraction f(esc), to simultaneously reproduce recent observations of the Thomson optical depth, ionizing emissivity and volume-averaged neutral fraction by the end of reionization. This yields f(csc) = 4(-2)(+7) averaged over our 0.375 h(-1) Mpc cells, independent of halo mass or redshift, increasing to 6 per cent if we also constrain to match the observed z = 7 star formation rate function. Introducing superlinear R-ion increases the duration of reionization and boosts small-scale 21 cm power by two to three times at intermediate phases of reionization, while inhomogeneous recombinations reduce ionized bubble sizes and suppress large-scale 21 cm power by two to three times. Gas clumping on sub-cell scales has a minimal effect on the 21 cm power. Superlinear R-ion also significantly increases the median halo mass scale for ionizing photon output to similar to 10(10) M-circle dot, making the majority of reionizing sources more accessible to next-generation facilities. These results highlight the importance of accurately treating ionizing sources and recombinations for modelling reionization and its 21 cm power spectrum.