Abstract The survey of Lyman α emitters (LAEs) with the Subaru Hyper Suprime-Cam, called SILVERRUSH (Ouchi et al. 2018, PASJ, 70, S13), is producing massive data of LAEs at z ≳ 6. Here we present LAE simulations to compare the SILVERRUSH data. In 1623 comoving Mpc3 boxes, where numerical radiative transfer calculations of reionization were performed, LAEs have been modeled with physically motivated analytic recipes as a function of halo mass. We have examined 23 models depending on the presence or absence of dispersion of halo Lyα emissivity, dispersion of the halo Lyα optical depth, τα, and halo mass dependence of τα. The unique free parameter in our model, a pivot value of τα, is calibrated so as to reproduce the z = 5.7 Lyα luminosity function (LF) of SILVERRUSH. We compare our model predictions with Lyα LFs at z = 6.6 and 7.3, LAE angular auto-correlation functions (ACFs) at z = 5.7 and 6.6, and LAE fractions in Lyman break galaxies at 5 < z < 7. The Lyα LFs and ACFs are reproduced by multiple models, but the LAE fraction turns out to be the most critical test. The dispersion of τα and the halo mass dependence of τα are essential to explain all observations reasonably. Therefore, a simple model of one-to-one correspondence between halo mass and Lyα luminosity with a constant Lyα escape fraction has been ruled out. Based on our best model, we present a formula to estimate the intergalactic neutral hydrogen fraction, $x_{\rm H\, \small {I}}$, from the observed Lyα luminosity density at z ≳ 6. We finally obtain $x_{\rm H\, \small {I}}=0.5_{-0.3}^{+0.1}$ as a volume-average at z = 7.3.
Read full abstract