We present the results from a new search for candidate galaxies at z ≈ 8.5–11 discovered over the 850 arcmin2 area probed by the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS). We use a photometric-redshift selection including both Hubble and Spitzer Space Telescope photometry to robustly identify galaxies in this epoch at H 160 < 26.6. We use a detailed vetting procedure, including screening against persistence and stellar contamination, and the inclusion of ground-based imaging and follow-up Hubble Space Telescope imaging to build a robust sample of 11 candidate galaxies, three presented here for the first time. The inclusion of Spitzer/IRAC photometry in the selection process reduces contamination, and yields more robust redshift estimates than Hubble alone. We constrain the evolution of the rest-frame ultraviolet luminosity function via a new method of calculating the observed number densities without choosing a prior magnitude bin size. We find that the abundance at our brightest probed luminosities (M UV = − 22.3) is consistent with predictions from simulations that assume that galaxies in this epoch have gas depletion times at least as short as those in nearby starburst galaxies. Due to large Poisson and cosmic variance uncertainties, we cannot conclusively rule out either a smooth evolution of the luminosity function continued from z = 4–8, or an accelerated decline at z > 8. We calculate that the presence of seven galaxies in a single field Extended Groth Strip is an outlier at the 2σ significance level, implying the discovery of a significant over-density. These scenarios will be imminently testable to high confidence within the first year of observations of the James Webb Space Telescope.

We present a detailed stellar population analysis of 11 bright (H < 26.6) galaxies at z = 9–11 (three spectroscopically confirmed) to constrain the chemical enrichment and growth of stellar mass of early galaxies. We use the flexible Bayesian spectral energy distribution (SED) fitting code Prospector with a range of star formation histories (SFHs), a flexible dust attenuation law, and a self-consistent model of emission lines. This approach allows us to assess how different priors affect our results and how well we can break degeneracies between dust attenuation, stellar ages, metallicity, and emission lines using data that probe only the rest-frame ultraviolet (UV) to optical wavelengths. We measure a median observed UV spectral slope β = − 1.87 − 0.43 + 0.35 for relatively massive star-forming galaxies ( 9 < log ( M ⋆ / M ⊙ ) < 10 ), consistent with no change from z = 4 to z = 9–10 at these stellar masses, implying rapid enrichment. Our SED-fitting results are consistent with a star-forming main sequence with sublinear slope (0.7 ± 0.2) and specific star formation rates of 3–10 Gyr−1. However, the stellar ages and SFHs are less well constrained. Using different SFH priors, we cannot distinguish between median mass-weighted ages of ∼ 50–150 Myr, which corresponds to 50% formation redshifts of z 50 ∼ 10–12 at z ∼ 9 and is of the order of the dynamical timescales of these systems. Importantly, models with different SFH priors are able to fit the data equally well. We conclude that the current observational data cannot tightly constrain the mass-buildup timescales of these z = 9–11 galaxies, with our results consistent with SFHs implying both a shallow and steep increase in the cosmic SFR density with time at z > 10.

Looking for bright galaxies born in the early universe is fundamental to investigating the Epoch of Reionization, the era when the first stars and galaxies ionized the intergalactic medium. We utilize Hubble Space Telescope pure parallel imaging to select galaxy candidates at a time 500 to 650 million years after the Big Bang, which corresponds to redshifts z ∼ 8−10. These data come from the Brightest of Reionizing Galaxies Survey (BoRG) Cycle 22 dataset, which consists of pure-parallel imaging in ∼ 90 different lines of sight that sum up to an area of ∼ 420 arcmin 2 . This survey uses five filters and has the advantage (compared to the Cycle 21 BoRG program) of including imaging in the JH 140 band, covering continuous wavelengths from the visible to near-infrared (λ = 0.35µm -1.7µm). This allows us to perform reliable selection of galaxies at z ≥ 8 using the photometric redshift technique. We use these galaxy candidates to constrain the bright end of the rest-frame ultraviolet luminosity function in this epoch. These candidates are excellent targets for follow-up observations, particularly with the James Webb Space Telescope.

We present the results from a spectroscopic survey using the MOSFIRE near-infrared spectrograph on the 10 m Keck telescope to search for Lyα emission from candidate galaxies at z ∼ 9–10 in four of the CANDELS fields (GOODS-N, EGS, UDS, and COSMOS). We observed 11 target galaxies, detecting Lyα from one object in ∼8.1 hr of integration, at z = 8.665 ± 0.001 with an integrated signal-to-noise ratio > 7. This galaxy is in the CANDELS Extended Groth Strip (EGS) field and lies physically close (3.5 physical Mpc [pMpc]) to another confirmed galaxy in this field with Lyα detected at z = 8.683. The detection of Lyα suggests the existence of large (∼1 pMpc) ionized bubbles fairly early in the reionization process. We explore the ionizing output needed to create bubbles of this size at this epoch and find that such a bubble requires more than the ionizing power provided by the full expected population of galaxies (by integrating the UV luminosity function down to M UV = −13). The Lyα we detect would be able to escape the predominantly neutral intergalactic medium at this epoch if our detected galaxy is inhabiting an overdensity, which would be consistent with the photometric overdensity previously identified in this region by Finkelstein et al. This implies that the CANDELS EGS field is hosting an overdensity at z = 8.7 that is powering one or more ionized bubbles, a hypothesis that will be imminently testable with forthcoming James Webb Space Telescope observations in this field.