We present an investigation into the first 500 Myr of galaxy evolution from the Cosmic Evolution Early Release Science (CEERS) survey. CEERS, one of 13 JWST ERS programs, targets galaxy formation from z ∼ 0.5 to >10 using several imaging and spectroscopic modes. We make use of the first epoch of CEERS NIRCam imaging, spanning 35.5 arcmin2, to search for candidate galaxies at z > 9. Following a detailed data reduction process implementing several custom steps to produce high-quality reduced images, we perform multiband photometry across seven NIRCam broad- and medium-band (and six Hubble broadband) filters focusing on robust colors and accurate total fluxes. We measure photometric redshifts and devise a robust set of selection criteria to identify a sample of 26 galaxy candidates at z ∼ 9–16. These objects are compact with a median half-light radius of ∼0.5 kpc. We present an early estimate of the z ∼ 11 rest-frame ultraviolet (UV) luminosity function, finding that the number density of galaxies at M UV ∼ −20 appears to evolve very little from z ∼ 9 to 11. We also find that the abundance (surface density [arcmin−2]) of our candidates exceeds nearly all theoretical predictions. We explore potential implications, including that at z > 10, star formation may be dominated by top-heavy initial mass functions, which would result in an increased ratio of UV light per unit halo mass, though a complete lack of dust attenuation and/or changing star formation physics may also play a role. While spectroscopic confirmation of these sources is urgently required, our results suggest that the deeper views to come with JWST should yield prolific samples of ultrahigh-redshift galaxies with which to further explore these conclusions.

We report the discovery of a candidate galaxy with a photo-z of z ∼ 12 in the first epoch of the James Webb Space Telescope (JWST) Cosmic Evolution Early Release Science Survey. Following conservative selection criteria, we identify a source with a robust z phot = 11.8 − 0.2 + 0.3 (1σ uncertainty) with m F200W = 27.3 and ≳7σ detections in five filters. The source is not detected at λ < 1.4 μm in deep imaging from both Hubble Space Telescope (HST) and JWST and has faint ∼3σ detections in JWST F150W and HST F160W, which signal a Lyα break near the red edge of both filters, implying z ∼ 12. This object (Maisie’s Galaxy) exhibits F115W − F200W > 1.9 mag (2σ lower limit) with a blue continuum slope, resulting in 99.6% of the photo-z probability distribution function favoring z > 11. All data-quality images show no artifacts at the candidate’s position, and independent analyses consistently find a strong preference for z > 11. Its colors are inconsistent with Galactic stars, and it is resolved (r h = 340 ± 14 pc). Maisie’s Galaxy has log M */M ⊙ ∼ 8.5 and is highly star-forming (log sSFR ∼ −8.2 yr−1), with a blue rest-UV color (β ∼ −2.5) indicating little dust, though not extremely low metallicity. While the presence of this source is in tension with most predictions, it agrees with empirical extrapolations assuming UV luminosity functions that smoothly decline with increasing redshift. Should follow-up spectroscopy validate this redshift, our universe was already aglow with galaxies less than 400 Myr after the Big Bang.

We report on the gas-phase metallicity gradients of a sample of 238 star-forming galaxies at 0.6 < z < 2.6, measured through deep near-infrared Hubble Space Telescope slitless spectroscopy. The observations include 12 orbit depth Hubble/WFC3 G102 grism spectra taken as a part of the CANDELS Lyα Emission at Reionization (CLEAR) survey, and archival WFC3 G102+G141 grism spectra overlapping the CLEAR footprint. The majority of galaxies in this sample are consistent with having a zero or slightly positive metallicity gradient (dZ/dR ≥ 0, i.e., increasing with radius) across the full mass range probed (8.5 < log M */M ⊙ < 10.5). We measure the intrinsic population scatter of the metallicity gradients, and show that it increases with decreasing stellar mass—consistent with previous reports in the literature, but confirmed here with a much larger sample. To understand the physical mechanisms governing this scatter, we search for correlations between the observed gradient and various stellar population properties at fixed mass. However, we find no evidence for a correlation with the galaxy properties we consider—including star formation rates, sizes, star formation rate surface densities, and star formation rates per gravitational potential energy. We use the observed weakness of these correlations to provide material constraints for predicted intrinsic correlations from theoretical models.