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 present the first results from the Mapping Obscuration to Reionization with ALMA (MORA) survey, the largest Atacama Large Millimeter/submillimeter Array (ALMA) blank-field contiguous survey to date (184 arcmin 2 ) and the only at 2 mm to search for dusty star-forming galaxies (DSFGs). We use the 13 sources detected above 5σ to estimate the first ALMA galaxy number counts at this wavelength. These number counts are then combined with the state-of-the-art galaxy number counts at 1.2 and 3 mm and with a backward evolution model to place constraints on the evolution of the IR luminosity function and dust-obscured star formation in the past 13 billion years. Our results suggest a steep redshift evolution on the space density of DSFGs and confirm the flattening of the IR luminosity function at faint luminosities, with a slope of a = --+ 0.42 LF 0.04 0.02 . We conclude that the dust-obscured component, which peaks at z ≈ 2-2.5, has dominated the cosmic history of star formation for the past ∼12 billion years, back to z ∼ 4. At z = 5, the dust-obscured star formation is estimated to be ∼35% of the total star formation rate density and decreases to 25%-20% at z = 6-7, implying a minor contribution of dustenshrouded star formation in the first billion years of the universe. With the dust-obscured star formation history constrained up to the end of the epoch of reionization, our results provide a benchmark to test galaxy formation models, to study the galaxy mass assembly history, and to understand the dust and metal enrichment of the universe at early times.

Recent simulations and observations of massive galaxy cluster evolution predict that the majority of stellar mass build up happens within cluster members by z = 2, before cluster virialization. Protoclusters rich with dusty, star-forming galaxies (DSFGs) at z > 3 are the favored candidate progenitors for these massive galaxy clusters at z ∼ 0. We present here the first study analyzing stellar emission along with cold dust and gas continuum emission in a spectroscopically confirmed z = 4.002 protocluster core rich with DSFGs, the Distant Red Core (DRC). We combine new HST and Spitzer data with existing Gemini, Herschel, and ALMA observations to derive individual galaxy-level properties, and compare them to coeval field and other protocluster galaxies. All of the protocluster members are massive (> 10 10 M ), but not significantly more so than their coeval field counterparts. Within uncertainty, all are nearly indistinguishable from galaxies on the star-forming vs. stellar mass main-sequence relationship, and on the star formation efficiency plane. Assuming no future major influx of fresh gas, we estimate that these gaseous DSFGs will deplete their gas reservoirs in ∼ 300 Myr, becoming the massive quiescent ellipticals dominating cluster cores by z ∼ 3. Using various methodologies, we derive a total z = 4 halo mass of ∼ 10 14 M , and estimate that the DRC will evolve to become an ultra-massive cluster core of mass 10 15 M by z = 0.

Lyman-break galaxy (LBG) candidates at z ≳ 10 are rapidly being identified in James Webb Space Telescope (JWST)/NIRCam observations. Due to the (redshifted) break produced by neutral hydrogen absorption of rest-frame UV photons, these sources are expected to drop out in the bluer filters while being well detected in redder filters. However, here we show that dust-enshrouded star-forming galaxies at lower redshifts (z ≲ 7) may also mimic the near-infrared (near-IR) colors of z > 10 LBGs, representing potential contaminants in LBG candidate samples. First, we analyze CEERS-DSFG-1, a NIRCam dropout undetected in the F115W and F150W filters but detected at longer wavelengths. Combining the JWST data with (sub)millimeter constraints, including deep NOEMA interferometric observations, we show that this source is a dusty star-forming galaxy (DSFG) at z ≈ 5.1. We also present a tentative 2.6σ SCUBA-2 detection at 850 μm around a recently identified z ≈ 16 LBG candidate in the same field and show that, if the emission is real and associated with this candidate, the available photometry is consistent with a z ∼ 5 dusty galaxy with strong nebular emission lines despite its blue near-IR colors. Further observations on this candidate are imperative to mitigate the low confidence of this tentative submillimeter emission and its positional uncertainty. Our analysis shows that robust (sub)millimeter detections of NIRCam dropout galaxies likely imply z ∼ 4–6 redshift solutions, where the observed near-IR break would be the result of a strong rest-frame optical Balmer break combined with high dust attenuation and strong nebular line emission, rather than the rest-frame UV Lyman break. This provides evidence that DSFGs may contaminate searches for ultra-high redshift LBG candidates from JWST observations.

We present the characteristics of 2 mm selected sources from the largest Atacama Large Millimeter/submillimeter Array (ALMA) blank-field contiguous survey conducted to date, the Mapping Obscuration to Reionization with ALMA (MORA) survey covering 184 arcmin2 at 2 mm. Twelve of 13 detections above 5σ are attributed to emission from galaxies, 11 of which are dominated by cold dust emission. These sources have a median redshift of 〈 z 2 mm 〉 = 3.6 − 0.3 + 0.4 primarily based on optical/near-infrared photometric redshifts with some spectroscopic redshifts, with 77% ± 11% of sources at z > 3 and 38% ± 12% of sources at z > 4. This implies that 2 mm selection is an efficient method for identifying the highest-redshift dusty star-forming galaxies (DSFGs). Lower-redshift DSFGs (z < 3) are far more numerous than those at z > 3 yet are likely to drop out at 2 mm. MORA shows that DSFGs with star formation rates in excess of 300 M ⊙ yr−1 and a relative rarity of ∼10−5 Mpc−3 contribute ∼30% to the integrated star formation rate density at 3 < z < 6. The volume density of 2 mm selected DSFGs is consistent with predictions from some cosmological simulations and is similar to the volume density of their hypothesized descendants: massive, quiescent galaxies at z > 2. Analysis of MORA sources’ spectral energy distributions hint at steeper empirically measured dust emissivity indices than reported in typical literature studies, with 〈 β 〉 = 2.2 − 0.4 + 0.5 . The MORA survey represents an important step in taking census of obscured star formation in the universe’s first few billion years, but larger area 2 mm surveys are needed to more fully characterize this rare population and push to the detection of the universe’s first dusty galaxies.