The GREATS H${\hat{I}}^2$+[OIII] Luminosity Function and Galaxy Properties at $\mathbf{z}\sim 8$: Walking the Way of JWST

Abstract

The James Webb Space Telescope will allow to spectroscopically study an unprecedented number of galaxies deep into the reionization era, notably by detecting [OIII] and H${\hat{I}}^2$ nebular emission lines. To efficiently prepare such observations, we photometrically select a large sample of galaxies at $z\sim 8$ and study their rest-frame optical emission lines. Combining data from the GOODS Re-ionization Era wide-Area Treasury from Spitzer (GREATS) survey and from HST, we perform spectral energy distribution (SED) fitting, using synthetic SEDs from a large grid of photoionization models. The deep Spitzer/IRAC data combined with our models exploring a large parameter space enables to constrain the [OIII]+H${\hat{I}}^2$ fluxes and equivalent widths for our sample, as well as the average physical properties of $z\sim 8$ galaxies, such as the ionizing photon production efficiency with $\log (\hat{I}{¾}_{\mathrm{ion}}/{\mathrm{erg}}^{-1}\mathrm{Hz})\ge 25.77$. We find a relatively tight correlation between the [OIII]+H${\hat{I}}^2$ and UV luminosity, which we use to derive for the first time the [OIII]+H${\hat{I}}^2$ luminosity function (LF) at $z\sim 8$. The $z\sim 8$ [OIII]+H${\hat{I}}^2$ LF is higher at all luminosities compared to lower redshift, as opposed to the UV LF, due to an increase of the [OIII]+H${\hat{I}}^2$ luminosity at a given UV luminosity from $z\sim 3$ to $z\sim 8$. Finally, using the [OIII]+H${\hat{I}}^2$ LF, we make predictions for JWST/NIRSpec number counts of $z\sim 8$ galaxies. We find that the current wide-area extragalactic legacy fields are too shallow to use JWST at maximal efficiency for $z\sim 8$ spectroscopy even at 1hr depth and JWST pre-imaging to $\gtrsim 30$ mag will be required.

13 pages, 9 figures, accepted for publication in MNRAS