Reionizing islands with inhomogeneous recombinations

Observations are beginning to constrain the history of the epoch of reionization (EoR). Modeling the reionization process is indispensable to interpret the observations, to infer the properties of ionizing sources, and to probe the various astrophysical processes from the observational data. Here we present an improved version of the semi-numerical simulation islandFAST, by incorporating inhomogeneous recombinations and a corresponding inhomogeneous ionizing background, and simulate the reionization process of neutral islands during the late EoR. We find that the islands are more fragmented in models with inhomogeneous recombinations than the case with a homogeneous recombination number. In order to investigate the effects of basic assumptions in the reionization modeling, we compare the results from islandFAST with those from 21cmFAST for the same assumptions on the ionizing photon sources and sinks, to find how the morphology of the ionization field and the reionization history depend on the different treatments of these two models. Such systematic bias should be noted when interpreting the upcoming observations.