Lyα-emitting galaxies as a probe of reionization: Large-scale bubble morphology and small-scale absorbers

The visibility of Lyα-emitting galaxies during the Epoch of Reionization is controlled by both diffuse HI patches in large-scale bubble morphology and small-scale absorbers. To investigate their impacts on Lyα transfer, we apply a novel combination of analytic modelling and cosmological hydrodynamical, radiative transfer simulations to three reionization models: (i) the 'bubble' model, where only diffuse HI outside ionized bubbles is present; (ii) the 'web' model, where HI exists only in overdense self-shielded gas; and (iii) the hybrid 'web-bubble' model. The three models can explain the observed Lyα luminosity function equally well, but with very different HI fractions. This confirms a degeneracy between the ionization topology of the intergalactic medium (IGM) and theHI fraction inferred from Lyα surveys.We highlight the importance of the clustering of small-scale absorbers around galaxies.Acombined analysis of the Lyα luminosity function and the Lyα fraction can break this degeneracy and provide constraints on the reionization history and its topology. Constraints can be improved by analysing the full MUV-dependent redshift evolution of the Lyα fraction of Lyman break galaxies. We find that the IGM-transmission probability distribution function is unimodal for bubble models and bimodal in web models. Comparing our models to observations, we infer that the neutral fraction at z ~ 7 is likely to be of the order of tens of per cent when interpreted with bubble or web-bubble models, with a conservative lower limit ~1 per cent when interpreted with web models.