SILVERRUSH. IV. Ly$α$ Luminosity Functions at $z = 5.7$ and $6.6$ Studied with $\sim$ 2,000 LAEs on the $14-21$ deg$^2$ Sky

Akira Konno, Masami Ouchi, Takatoshi Shibuya, Yoshiaki Ono, Kazuhiro Shimasaku, Yoshiaki Taniguchi, Tohru Nagao, Masakazu A. R. Kobayashi, Masaru Kajisawa, Nubunari Kashikawa, Akio K. Inoue, Masamune Oguri, Hisanori Furusawa, Tomotsugu Goto, Yuichi Harikane, Ryo Higuchi, Yutaka Komiyama, Haruka Kusakabe, Satoshi Miyazaki, Kimihiko Nakajima, Shiang-Yu Wang

Published 2017-05-03Version 1

We present the Ly$\alpha$ luminosity functions (LFs) at $z=$5.7 and 6.6 derived from a new large sample of 2,354 Ly$\alpha$ emitters (LAEs) identified in total areas of 14 and 21 deg$^2$, respectively, based on the early narrowband data of the Subaru/Hyper Suprime-Cam (HSC) survey. Together with careful Monte-Carlo simulations that account for the incompleteness of the LAE selection and the flux estimate systematics in the narrowband imaging, we have determined the Ly$\alpha$ LFs with the unprecedentedly small statistical and systematic uncertainties in a wide Ly$\alpha$ luminosity range of $10^{42.8-43.8}$ erg s$^{-1}$. We obtain the best-fit Schechter parameters of $L^{*}_{{\rm Lya}}=1.1^{+0.7}_{-0.4}\times10^{43}\ (8.2^{+8.6}_{-3.0}\times10^{42})$ erg s$^{-1}$, $\phi^{*}_{{\rm Lya}}=2.5^{+3.4}_{-1.9}\ (2.8^{+3.6}_{-2.3})\times10^{-4}$ Mpc$^{-3}$, and $\alpha=-2.3^{+0.8}_{-0.4}\ (-1.9^{+0.8}_{-0.6})$ at $z=5.7$ ($6.6$). We confirm that our best-estimate Ly$\alpha$ LFs are consistent with the majority of the previous studies, but find that our Ly$\alpha$ LFs do not agree with the large Ly$\alpha$ LFs recently claimed by Matthee/Santos et al.'s studies that may overcorrect the incompleteness and the flux systematics. Our Ly$\alpha$ LFs at $z=$5.7 and 6.6 show a very steep slope ($\alpha\simeq-2$) or a significant ($\simeq3\sigma$) excess beyond the best-fit Schechter function at the bright end ($\gtrsim10^{43.5}$ erg s$^{-1}$) that is claimed by previous studies. This bright-end LF excess at $z\sim6-7$ may be explained by the contribution from AGNs, blended merging galaxies, and/or large ionized bubbles around bright LAEs. Comparing our Ly$\alpha$ LF measurements with four independent reionization models, we estimate the neutral hydrogen fraction of the IGM to be $x_{\rm HI}=0.3\pm0.2$ at $z=$6.6 that is consistent with the small Thomson scattering optical depth obtained by Planck 2016.