Not yet reviewed by Pith; the record is open.

With a total integration time of 168 hours and a narrowband (NB) filter tuned to Lyman alpha at z = 8.8, the UltraVISTA survey has set out to find some of the most distant galaxies, on the verge of the Epoch of Reionization. Previous calculations of the expected number of detected Lya-emitting galaxies (LAEs) at this redshift did not explicitly take into account the radiative transfer (RT) of Lya. In this work we combine a theoretical model for the halo mass function with numerical results from high-res cosmological hydrosimulations with LyC+Lya RT, assessing the visibility of LAEs residing in these halos. Uncertainties such as cosmic variance and the anisotropic escape of Lya are taken into account, and it is predicted that once the survey has finished, the probabilities of detecting none, one, or more than one are ~90%, ~10%, and ~1%; a significantly smaller success rate compared to earlier predictions, due to the combined effect of a highly neutral IGM scattering Lya to such large distances from the galaxy that they fall outside the observational aperture, and to the actual depth of the survey being less than predicted. Because the IGM affects NB and broadband (BB) magnitudes differently, we argue for a relaxed color selection criterion of NB - BB ~ +0.85. But since the flux is continuum-dominated, even if a galaxy is detectable in the NB its probability of being selected as a NB excess object is <~35%. Various properties of galaxies at this redshift are predicted, e.g. UV and Lya LFs, M*-Mh relation, spectral shape, optimal aperture, and the anisotropic escape of Lya through both a dusty ISM and a partly neutral IGM. Finally, we describe and publish a fast numerical code for adding numbers with asymmetric uncertainties ("x_{-sigma_1}^{+sigma_2}") proving to be significantly better than the standard, but wrong, way of adding upper and lower uncertainties in quadrature separately.