Unresolved z~8 point sources and their impact on the bright end of the galaxy luminosity function

The distribution and properties of the first galaxies and quasars are critical pieces of the puzzle in understanding galaxy evolution and cosmic reionization. Previous studies have often excluded unresolved sources as potential low redshift interlopers. We combine broadband color and photometric redshift analysis with morphological selections to identify a robust sample of candidates consistent with unresolved point sources at redshift $z\sim8$ using deep Hubble Space Telescope images. We also examine G141 grism spectroscopic data to identify and eliminate dwarf star contaminants. From these analyses, we identify three, bright ($M_{UV}\lesssim-22$ ABmag) dropout point sources at $7.5<z<8.1$. Spectral energy distribution analyses suggest that these sources are either quasars or compact star-forming galaxies. The flux captured by the IRAC 4.5 $\mu$m channel suggests that they have moderate $H\beta$+$[OIII]$ equivalent widths. We calculate the number density of point sources at $z\sim7$-8, and find that a double powerlaw model well describes the point source distribution. We then extend our analysis to estimate the combined point source + galaxy luminosity function and find that the point sources have a non-negligible contribution to the bright-end excess. The fact that the point sources dominate only at $M_{UV}\lesssim-22$ suggests that their contribution to cosmic reionization is likely limited. While spectroscopic follow-up is needed to confirm the nature of these point sources, this work demonstrates that the inclusion of Lyman dropout point sources is necessary for a complete census of the early galaxies at the epoch of cosmic reionization.