2015MNRAS.452..261C

We combine high-resolution hydrodynamical simulations with an intermediate resolution, dark matter only simulation and an analytical model for the growth of ionized regions to estimate the large-scale distribution and redshift evolution of the visibility of Lyα emission in 6 ≤ z ≤ 8 galaxies. The inhomogeneous distribution of neutral hydrogen during the reionization process results in significant fluctuations in the Lyα transmissivity on large scales. The transmissivity depends not only on the ionized fraction of the intergalactic medium by volume and the amplitude of the local ionizing background, but is also rather sensitive to the evolution of the relative velocity shift of the Lyα emission line due to resonant scattering. We reproduce a decline in the space density of Lyα emitting galaxies as rapid as observed with a rather rapidly evolving neutral fraction between z = 6-8, and a typical Lyα line velocity offset of 100kms- 1 redward of systemic at z = 6 which decreases towards higher redshift. The new (02/2015) Planck results indicate such a recent end to reionization is no longer disfavoured by constraints from the cosmic microwave background.