Search

We present spectroscopic confirmation of an ultra-massive galaxy (UMG) with $\log(M_\star/M_\odot) = 10.98 \pm 0.07$ at $z_\mathrm{spec} = 4.8947$ in the Extended Groth Strip (EGS), based on deep observations of Ly$\alpha$ emission with Keck/DEIMOS. The ultra-massive galaxy (UMG-28740) is the most massive member in one of the most significant overdensities in the EGS, with four additional photometric members with $\log(M_\star/M_\odot) > 10.5$ within $R_\mathrm{proj} \sim 1$ cMpc. Spectral energy distribution (SED) fitting using a large suite of star formation histories and two sets of high-quality photometry from ground- and space-based facilities consistently estimates the mass of this object to be $\log(M_\star/M_\odot) \sim 11$ with a small standard deviation between measurements ($\sigma = 0.07$). While the best-fit SED models agree on stellar mass, we find discrepancies in the estimated star formation rate for UMG-28740, resulting in either a star-forming or quiescent system. $\mathit{JWST}$/NIRCam photometry of UMG-28740 strongly favors a quiescent scenario, demonstrating the need for high-quality mid-IR observations. Assuming the galaxy to be quiescent, UMG-28740 formed the bulk of its stars at $z > 10$ and is quenching at $z \sim 8$, resulting in a high star formation efficiency at high redshift ($\epsilon \sim 0.2$ at $z \sim 5$ and $\epsilon \gtrsim 1$ at $z \gtrsim 8$). As the most massive galaxy in its protocluster environment, UMG-28740 is a unique example of the impossibly early galaxy problem.

We present a sample of 88 candidate z ∼ 8.5–14.5 galaxies selected from the completed NIRCam imaging from the Cosmic Evolution Early Release Science survey. These data cover ∼90 arcmin ^2 (10 NIRCam pointings) in six broadband imaging filters and one medium-band imaging filter. With this sample we confirm at higher confidence early JWST conclusions that bright galaxies in this epoch are more abundant than predicted by most theoretical models. We construct the rest-frame ultraviolet luminosity functions at z ∼ 9, 11, and 14 and show that the space density of bright ( M _UV = −20) galaxies changes only modestly from z ∼ 14 to z ∼ 9, compared to a steeper increase from z ∼ 8 to z ∼ 4. While our candidates are photometrically selected, spectroscopic follow-up has now confirmed 13 of them, with only one significant interloper, implying that the fidelity of this sample is high. Successfully explaining the evidence for a flatter evolution in the number densities of UV-bright z > 10 galaxies may thus require changes to the dominant physical processes regulating star formation. While our results indicate that significant variations of dust attenuation with redshift are unlikely to be the dominant factor at these high redshifts, they are consistent with predictions from models that naturally have enhanced star formation efficiency and/or stochasticity. An evolving stellar initial mass function could also bring model predictions into better agreement with our results. Deep spectroscopic follow-up of a large sample of early galaxies can distinguish between these competing scenarios.

Many quiescent galaxies discovered in the early Universe by JWST raise fundamental questions on when and how these galaxies became and stayed quenched. Making use of the latest version of the semianalytic model GAEA that provides good agreement with the observed quenched fractions up to z ∼ 3, we make predictions for the expected fractions of quiescent galaxies up to z ∼ 7 and analyze the main quenching mechanism. We find that in a simulated box of 685 Mpc on a side, the first quenched massive ( M _⋆ ∼ 10 ^11 M _⊙ ), Milky Way–mass, and low-mass ( M _⋆ ∼ 10 ^9.5 M _⊙ ) galaxies appear at z ∼ 4.5, z ∼ 6.2, and before z = 7, respectively. Most quenched galaxies identified at early redshifts remain quenched for more than 1 Gyr. Independently of galaxy stellar mass, the dominant quenching mechanism at high redshift is accretion disk feedback (quasar winds) from a central massive black hole, which is triggered by mergers in massive and Milky Way–mass galaxies and by disk instabilities in low-mass galaxies. Environmental stripping becomes increasingly more important at lower redshift.

The changes in colors across a galaxy are intimately connected to the galaxy’s formation, growth, quenching history, and dust content. A particularly important epoch in the growth of galaxies is near z ∼ 2, often referred to as “cosmic noon,” where galaxies on average reach the peak of their star formation. We study a population of 125 cluster galaxies at z ∼ 1.6 in three Hubble Space Telescope filters, F475W, F625W, and F160W, roughly corresponding to the rest-frame far-ultraviolet, near-ultraviolet, and r band, respectively. By comparing to a control sample of 200 field galaxies at similar redshift, we reveal clear, statistically significant differences in the overall spatially resolved colors and color gradients in galaxies across these two different environments. On average, cluster galaxies have redder ultraviolet colors in both the inner and outer regions bounded by r _50 , as well as an overall wider dispersion of outside-in color gradients. The presence of these observed differences, along with evidence from ancillary data from previous studies, strongly suggests that the environment drives these population-level color differences, by affecting the stellar populations and/or dust content.

The Local Volume Complete Cluster Survey is an ongoing program to observe nearly a hundred low-redshift X-ray-luminous galaxy clusters (redshifts 0.03 < z < 0.12 and X-ray luminosities in the 0.1–2.4 keV band L _X500c > 10 ^44 erg s ^−1 ) with the Dark Energy Camera, capturing data in the u , g , r , i , z bands with a 5 σ point source depth of approximately 25th–26th AB magnitudes. Here, we map the aperture masses in 58 galaxy cluster fields using weak gravitational lensing. These clusters span a variety of dynamical states, from nearly relaxed to merging systems, and approximately half of them have not been subject to detailed weak lensing analysis before. In each cluster field, we analyze the alignment between the 2D mass distribution described by the aperture mass map, the 2D red-sequence (RS) galaxy distribution, and the brightest cluster galaxy (BCG). We find that the orientations of the BCG and the RS distribution are strongly aligned throughout the interiors of the clusters: the median misalignment angle is 19° within 2 Mpc. We also observe the alignment between the orientations of the RS distribution and the overall cluster mass distribution (by a median difference of 32° within 1 Mpc), although this is constrained by galaxy shape noise and the limitations of our cluster sample size. These types of alignment suggest long-term dynamical evolution within the clusters over cosmic timescales.

We analyze the evolution of massive (log _10 [ M _⋆ / M _⊙ ] > 10) galaxies at z ∼ 1–4 selected from JWST Cosmic Evolution Early Release Survey (CEERS). We infer the physical properties of all galaxies in the CEERS NIRCam imaging through spectral energy distribution (SED) fitting with dense basis to select a sample of high-redshift massive galaxies. Where available we include constraints from additional CEERS observing modes, including 18 sources with MIRI photometric coverage, and 28 sources with spectroscopic confirmations from NIRSpec or NIRCam WFSS. We sample the recovered posteriors in stellar mass from SED fitting to infer the volume densities of massive galaxies across cosmic time, taking into consideration the potential for sample contamination by active galactic nuclei. We find that the evolving abundance of massive galaxies tracks expectations based on a constant baryon conversion efficiency in dark matter halos for z ∼ 1–4. At higher redshifts, we observe an excess abundance of massive galaxies relative to this simple model, resulting in a shallower decline of observed volume densities of massive galaxies. These higher abundances can be explained by modest changes to star formation physics and/or the efficiencies with which star formation occurs in massive dark matter halos, and are not in tension with modern cosmology.

We present an analysis of the number density of galaxies as a function of stellar mass (i.e., the stellar mass function (SMF)) in the COSMOS field at z ∼ 3.3, making a comparison between the SMF in overdense environments and the SMF in the coeval field. In particular, this region contains the Elentári proto-supercluster, a system of six extended overdensities spanning ∼70 cMpc on a side. A clear difference is seen in the high-mass slope of these SMFs, with overdense regions showing an increase in the ratio of high-mass galaxies to low-mass galaxies relative to the field, indicating a more rapid buildup of stellar mass in overdense environments. This result qualitatively agrees with analyses of clusters at z ∼ 1, though the differences between protocluster and field SMFs at z ∼ 3.3 are smaller. While this is consistent with overdensities enhancing the evolution of their member galaxies, potentially through increased merger rates, whether this enhancement begins in protocluster environments or even earlier in group environments is still unclear. Though the measured fractions of quiescent galaxies between the field and overdense environments do not vary significantly, implying that this stellar mass enhancement is ongoing and any starbursts triggered by merger activity have not yet quenched, we note that spectroscopic observations are biased toward star-forming populations, particularly for low-mass galaxies. If mergers are indeed responsible, high-resolution imaging of Elentári and similar structures at these early epochs should then reveal increased merger rates relative to the field. Larger samples of well-characterized overdensities are necessary to draw broader conclusions in these areas.

Substantial populations of massive quiescent galaxies at z ≥ 3 challenge our understanding of rapid galaxy growth and quenching over short timescales. In order to piece together this evolutionary puzzle, more statistical samples of these objects are required. Established techniques for identifying massive quiescent galaxies are increasingly inefficient and unconstrained at z > 3. As a result, studies report that as much as 70% of quiescent galaxies at z > 3 may be missed from existing surveys. In this work, we propose a new empirical color selection technique designed to select massive quiescent galaxies at 3 ≲ z ≲ 6 using JWST NIRCam imaging data. We use empirically constrained galaxy spectral energy distribution (SED) templates to define a region in the F277W − F444W versus F150W − F277W color plane that captures quiescent galaxies at z > 3. We apply these color selection criteria to the Cosmic Evolution Early Release Science (CEERS) Survey and use SED fitting on sources in the region to identify 44 candidate z ≳ 3 quiescent galaxies. Over half of these sources are newly discovered and, on average, exhibit specific star formation rates of poststarburst galaxies. Most of these sources would not be discovered using canonical UVJ diagrams. We derive volume density estimates of n ∼ 1–4 × 10 ^−5 Mpc ^−3 at 3 < z < 5, finding excellent agreement with existing reports on similar populations in the CEERS field. Thanks to NIRCam’s wavelength coverage and sensitivity, this technique provides an efficient tool to search for large samples of these rare galaxies.

We present the first results from the Web Epoch of Reionization Ly α Survey (WERLS), a spectroscopic survey of Ly α emission using Keck I/MOSFIRE and LRIS. WERLS targets bright ( J < 26) galaxy candidates with photometric redshifts of 5.5 ≲ z ≲ 8 selected from pre-JWST imaging embedded in the Epoch of Reionization (EoR) within three JWST deep fields: CEERS, PRIMER, and COSMOS-Web. Here, we report 11 z ∼ 7–8 Ly α emitters (LAEs; three secure and eight tentative candidates) detected in the first five nights of WERLS MOSFIRE data. We estimate our observed LAE yield is ∼13%, which is broadly consistent with expectations assuming some loss from redshift uncertainty, contamination from sky OH lines, and that the Universe is approximately half-ionized at this epoch, whereby observable Ly α emission is unlikely for galaxies embedded in a neutral intergalactic medium. Our targets are selected to be UV-bright, and span a range of absolute UV magnitudes with −23.1 < M _UV < −19.8. With two LAEs detected at z = 7.68, we also consider the possibility of an ionized bubble at this redshift. Future synergistic Keck+JWST efforts will provide a powerful tool for pinpointing beacons of reionization and mapping the large-scale distribution of mass relative to the ionization state of the Universe.

We analyze rest-frame ultraviolet to optical spectra of three z ≃ 7.47–7.75 galaxies whose Ly α emission lines were previously detected with Keck/MOSFIRE observations, using the JWST/NIRSpec observations from the Cosmic Evolution Early Release Science survey. From NIRSpec data, we confirm the systemic redshifts of these Ly α emitters, and emission-line ratio diagnostics indicate these galaxies were highly ionized and metal-poor. We investigate Ly α line properties, including the line flux, velocity offset, and spatial extent. For the one galaxy where we have both NIRSpec and MOSFIRE measurements, we find a significant offset in their flux measurements (∼1.3–5× greater in MOSFIRE) and a marginal difference in the velocity shifts. The simplest interpretation is that the Ly α emission is extended and not entirely encompassed by the NIRSpec slit. The cross-dispersion profiles in NIRSpec reveal that Ly α in one galaxy is significantly more extended than the nonresonant emission lines. We also compute the expected sizes of ionized bubbles that can be generated by the Ly α sources and discuss viable scenarios for the creation of sizable ionized bubbles (>1 physical Mpc). The source with the highest-ionization condition is possibly capable of ionizing its own bubble, while the other two do not appear to be capable of ionizing such a large region, but require additional sources of ionizing photons. Therefore, the fact that we detect Ly α from these galaxies suggests diverse scenarios for the escape of Ly α during the epoch of reionization. High-spectral-resolution spectra with JWST/NIRSpec will be extremely useful for constraining the physics of patchy reionization.

We report on the discovery of two low-luminosity, broad-line active galactic nuclei (AGNs) at z > 5 identified using JWST NIRSpec spectroscopy from the Cosmic Evolution Early Release Science (CEERS) survey. We detect broad H α emission in the spectra of both sources, with FWHM of 2060 ± 290 km s ^−1 and 1800 ± 200 km s ^−1 , resulting in virial black hole (BH) masses that are 1–2 dex below those of existing samples of luminous quasars at z > 5. The first source, CEERS 2782 at z = 5.242, is 2–3 dex fainter than known quasars at similar redshifts and was previously identified as a candidate low-luminosity AGN based on its morphology and rest-frame optical spectral energy distribution (SED). We measure a BH mass of M _BH = (1.3 ± 0.4) × 10 ^7 M _⊙ , confirming that this AGN is powered by the least massive BH known in the Universe at the end of cosmic reionization. The second source, CEERS 746 at z = 5.624, is inferred to be a heavily obscured, broad-line AGN caught in a transition phase between a dust-obscured starburst and an unobscured quasar. We estimate its BH mass to be in the range of M _BH ≃ (0.9–4.7) × 10 ^7 M _⊙ , depending on the level of dust obscuration assumed. We perform SED fitting to derive host stellar masses, M _⋆ , allowing us to place constraints on the BH–galaxy mass relationship in the lowest mass range yet probed in the early Universe. The M _BH / M _⋆ ratio for CEERS 2782, in particular, is consistent with or higher than the empirical relationship seen in massive galaxies at z = 0. We examine the narrow emission line ratios of both sources and find that their location on the BPT and OHNO diagrams is consistent with model predictions for moderately low metallicity AGNs with Z / Z _⊙ ≃ 0.2–0.4. The spectroscopic identification of low-luminosity, broad-line AGNs at z > 5 with M _BH ≃ 10 ^7 M _⊙ demonstrates the capability of JWST to push BH masses closer to the range predicted for the BH seed population and provides a unique opportunity to study the early stages of BH–galaxy assembly.

We present JWST NIRSpec spectroscopy for 11 galaxy candidates with photometric redshifts of z ≃ 9 − 13 and M _UV ∈ [ −21, −18] newly identified in NIRCam images in the Cosmic Evolution Early Release Science Survey. We confirm emission line redshifts for 7 galaxies at z = 7.762–8.998 using spectra at ∼1–5 μ m either with the NIRSpec prism or its three medium-resolution ( R ∼ 1000) gratings. For z ≃ 9 photometric candidates, we achieve a high confirmation rate of ≃90%, which validates the classical dropout selection from NIRCam photometry. No robust emission lines are identified in three galaxy candidates at z > 10, where the strong [O iii ] and H β lines would be redshifted beyond the wavelength range observed by NIRSpec, and the Ly α continuum break is not detected with the sensitivity of the current data. Compared with Hubble Space Telescope-selected bright galaxies ( M _UV ≃ −22) that are similarly spectroscopically confirmed at z ≃ 8 − 9, these NIRCam-selected galaxies are characterized by lower star formation rates (SFRs; SFR ≃ 4 M _⊙ yr ^−1 ) and lower stellar masses (≃10 ^8 M _⊙ ), but with higher specific SFR (≃40 Gyr ^−1 ), higher [O iii ]+H β equivalent widths (≃1100 Å), and elevated production efficiency of ionizing photons ( $\mathrm{log}({\xi }_{\mathrm{ion}}/\mathrm{Hz}\,{\mathrm{erg}}^{-1})\simeq 25.8$ ) induced by young stellar populations (

We present results from the Cosmic Evolution Early Release Survey on the stellar population parameters for 28 galaxies with redshifts 4 < z < 9 using imaging data from the James Webb Space Telescope (JWST) Mid-Infrared Instrument (MIRI) combined with data from the Hubble Space Telescope and the Spitzer Space Telescope. The JWST/MIRI 5.6 and 7.7 μ m data extend the coverage of the rest-frame spectral energy distribution to nearly 1 μ m for galaxies in this redshift range. By modeling the galaxies’ SEDs the MIRI data show that the galaxies have, on average, rest-frame UV (1600 Å)— I -band colors 0.4 mag bluer than derived when using photometry that lacks MIRI. Therefore, the galaxies have lower ratios of stellar mass to light. The MIRI data reduce the stellar masses by $\langle {\rm{\Delta }}\mathrm{log}{M}_{* }\rangle =0.25$ dex at 4 < z < 6 and 0.37 dex at 6 < z < 9. This also reduces the star formation rates (SFRs) by 〈ΔlogSFR〉 = 0.14 dex at 4 < z < 6 and 0.27 dex at 6 < z < 9. The MIRI data also improve constraints on the allowable stellar mass formed in early star formation. We model this using a star formation history that includes both a “burst” at z _f = 100 and a slowly varying (“delayed- τ ”) model. The MIRI data reduce the allowable stellar mass by 0.6 dex at 4 < z < 6 and by ≈1 dex at 6 < z < 9. Applying these results globally, this reduces the cosmic stellar-mass density by an order of magnitude in the early Universe ( z ≈ 9). Therefore, observations of rest-frame ≳1 μ m are paramount for constraining the stellar-mass buildup in galaxies at very high redshifts.

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 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 arcmin ^2 , 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.

Stellar bars are key drivers of secular evolution in galaxies and can be effectively studied using rest-frame near-infrared (NIR) images, which trace the underlying stellar mass and are less impacted by dust and star formation than rest-frame UV or optical images. We leverage the power of JWST CEERS NIRCam images to present the first quantitative identification and characterization of stellar bars at z > 1 based on rest-frame NIR F444W images of high resolution (∼1.3 kpc at z ∼ 1–3). We identify stellar bars in these images using quantitative criteria based on ellipse fits. For this pilot study, we present six examples of robustly identified bars at z > 1 with spectroscopic redshifts, including the two highest-redshift bars at z ∼ 2.136 and 2.312 quantitatively identified and characterized to date. The stellar bars at z ∼ 1.1–2.3 presented in our study have projected semimajor axes of ∼2.9–4.3 kpc and projected ellipticities of ∼0.41–0.53 in the rest-frame NIR. The barred host galaxies have stellar masses ∼1 × 10 ^10 to 2 × 10 ^11 M _⊙ and star formation rates of ∼21–295 M _⊙ yr ^−1 , and several have potential nearby companions. Our finding of bars at z ∼ 1.1–2.3 demonstrates the early onset of such instabilities and supports simulations where bars form early in massive dynamically cold disks. It also suggests that if these bars at lookback times of 8–11 Gyr survive out to present epochs, bar-driven secular processes may operate over a long time and have a significant impact on some galaxies by z ∼ 0.

We present a comprehensive analysis of the evolution of the morphological and structural properties of a large sample of galaxies at z = 3–9 using early James Webb Space Telescope (JWST) CEERS NIRCam observations. Our sample consists of 850 galaxies at z > 3 detected in both Hubble Space Telescope (HST)/WFC3 and CEERS JWST/NIRCam images, enabling a comparison of HST and JWST morphologies. We conduct a set of visual classifications, with each galaxy in the sample classified three times. We also measure quantitative morphologies across all NIRCam filters. We find that galaxies at z > 3 have a wide diversity of morphologies. Galaxies with disks make up 60% of galaxies at z = 3, and this fraction drops to ∼30% at z = 6–9, while galaxies with spheroids make up ∼30%–40% across the redshift range, and pure spheroids with no evidence for disks or irregular features make up ∼20%. The fraction of galaxies with irregular features is roughly constant at all redshifts (∼40%–50%), while those that are purely irregular increases from ∼12% to ∼20% at z > 4.5. We note that these are apparent fractions, as many observational effects impact the visibility of morphological features at high redshift. On average, Spheroid-only galaxies have a higher Sérsic index, smaller size, and higher axis ratio than disk or irregular galaxies. Across all redshifts, smaller spheroid and disk galaxies tend to be rounder. Overall, these trends suggest that galaxies with established disks and spheroids exist across the full redshift range of this study, and further work with large samples at higher redshift is needed to quantify when these features first formed.

We present the data release and data reduction process for the Epoch 1 NIRCam observations for the Cosmic Evolution Early Release Science Survey (CEERS). These data consist of NIRCam imaging in six broadband filters (F115W, F150W, F200W, F277W, F356W and F444W) and one medium-band filter (F410M) over four pointings, obtained in parallel with primary CEERS MIRI observations. We reduced the NIRCam imaging with the JWST Calibration Pipeline, with custom modifications and reduction steps designed to address additional features and challenges with the data. Here we provide a detailed description of each step in our reduction and a discussion of future expected improvements. Our reduction process includes corrections for known prelaunch issues such as 1/ f noise, as well as in-flight issues including snowballs, wisps, and astrometric alignment. Many of our custom reduction processes were first developed with prelaunch simulated NIRCam imaging over the full 10 CEERS NIRCam pointings. We present a description of the creation and reduction of this simulated data set in the Appendix. We provide mosaics of the real images in a public release, as well as our reduction scripts with detailed explanations to allow users to reproduce our final data products. These represent one of the first official public data sets released from the Directors Discretionary Early Release Science (DD-ERS) program.

We report on the host properties of five X-ray-luminous active galactic nuclei (AGN) identified at 3 < z < 5 in the first epoch of imaging from the Cosmic Evolution Early Release Science Survey. Each galaxy has been imaged with the JWST Near-Infrared Camera, which provides rest-frame optical morphologies at these redshifts. We also derive stellar masses and star formation rates for each host by fitting its spectral energy distribution using a combination of galaxy and AGN templates. We find that three of the AGN hosts have spheroidal morphologies, one is a bulge-dominated disk, and one is dominated by pointlike emission. None are found to show strong morphological disturbances that might indicate a recent interaction or merger event. When compared to a sample of mass-matched inactive galaxies, we find that the AGN hosts have morphologies that are less disturbed and more bulge-dominated. Notably, all four of the resolved hosts have rest-frame optical colors consistent with a quenched or poststarburst stellar population. The presence of AGN in passively evolving galaxies at z > 3 is significant because a rapid feedback mechanism is required in most semianalytic models and cosmological simulations to explain the growing population of massive quiescent galaxies observed at these redshifts. Our findings show that AGN can continue to inject energy into these systems after their star formation is curtailed, potentially heating their halos and preventing renewed star formation. Additional observations will be needed to determine what role this feedback may play in helping to quench these systems and/or maintain their quiescent state.

The new capabilities that JWST offers in the near- and mid-infrared (IR) are used to investigate in unprecedented detail the nature of optical/near-IR-faint, mid-IR-bright sources, with HST-dark galaxies among them. We gather JWST data from the CEERS survey in the Extended Groth Strip, jointly with HST data, and analyze spatially resolved optical-to-mid-IR spectral energy distributions to estimate photometric redshifts in two dimensions and stellar population properties on a pixel-by-pixel basis for red galaxies detected by NIRCam. We select 138 galaxies with F150W − F356W > 1.5 mag and F356W < 27.5 mag. The nature of these sources is threefold: (1) 71% are dusty star-forming galaxies (SFGs) at 2 < z < 6 with $9\lt \mathrm{log}{M}_{\star }/{M}_{\odot }\lt 11$ and a variety of specific SFRs (100 Gyr ^−1 ); (2) 18% are quiescent/dormant (i.e., subject to reignition/rejuvenation) galaxies (QGs) at 3 < z < 5, with $\mathrm{log}{M}_{\star }/{M}_{\odot }\sim 10$ and poststarburst mass-weighted ages (0.5–1.0 Gyr); and (3) 11% are strong young starbursts with indications of high equivalent width emission lines (typically, [O iii ]+H β ) at 6 < z < 7 (XELG- z 6) and $\mathrm{log}{M}_{\star }/{M}_{\odot }\sim 9.5$ . The sample is dominated by disk-like galaxies with remarkable compactness for XELG- z 6 (effective radii smaller than 0.4 kpc). Large attenuations in SFGs, 2 < A ( V ) < 5 mag, are found within 1.5 times the effective radius, approximately 2 kpc, while QGs present A ( V ) ∼ 0.2 mag. Our SED-fitting technique reproduces the expected dust emission luminosities of IR-bright and submillimeter galaxies. This study implies high levels of star formation activity between z ∼ 20 and z ∼ 10, where virtually 100% of our galaxies had already formed 10 ^8 M _⊙ , 60% had assembled 10 ^9 M _⊙ , and 10% up to 10 ^10 M _⊙ (in situ or ex situ).

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 present the lifetime star formation histories (SFHs) for six ultrafaint dwarf (UFD; M _V > − 7.0, $4.9\lt {\mathrm{log}}_{10}({M}_{* }(z=0)/{M}_{\odot })\lt 5.5$ ) satellite galaxies of M31 based on deep color–magnitude diagrams constructed from Hubble Space Telescope imaging. These are the first SFHs obtained from the oldest main-sequence turnoff of UFDs outside the halo of the Milky Way (MW). We find that five UFDs formed at least 50% of their stellar mass by z = 5 (12.6 Gyr ago), similar to known UFDs around the MW, but that 10%–40% of their stellar mass formed at later times. We uncover one remarkable UFD, And xiii , which formed only 10% of its stellar mass by z = 5, and 75% in a rapid burst at z ∼ 2–3, a result that is robust to choices of underlying stellar model and is consistent with its predominantly red horizontal branch. This “young” UFD is the first of its kind and indicates that not all UFDs are necessarily quenched by reionization, which is consistent with predictions from several cosmological simulations of faint dwarf galaxies. SFHs of the combined MW and M31 samples suggest reionization did not homogeneously quench UFDs. We find that the least-massive MW UFDs ( M _* ( z = 5) ≲ 5 × 10 ^4 M _⊙ ) are likely quenched by reionization, whereas more-massive M31 UFDs ( M _* ( z = 5) ≳ 10 ^5 M _⊙ ) may only have their star formation suppressed by reionization and quench at a later time. We discuss these findings in the context of the evolution and quenching of UFDs.

We present Atacama Large Millimeter/submillimeter Array (ALMA) Band 7 observations of a remarkably bright galaxy candidate at ${z}_{\mathrm{phot}}={16.7}_{-0.3}^{+1.9}$ ( M _UV = −21.6), S5-z17-1, identified in James Webb Space Telescope (JWST) Early Release Observation data of Stephen’s Quintet. We do not detect the dust continuum at 866 μ m, ruling out the possibility that S5-z17-1 is a low- z dusty starburst with a star formation rate of ≳30 M _⊙ yr ^−1 . We detect a 5.1 σ line feature at 338.726 ± 0.007 GHz exactly coinciding with the JWST source position, with a 2% likelihood of the signal being spurious. The most likely line identification would be [O iii ]52 μ m at z = 16.01 or [C ii ]158 μ m at z = 4.61, whose line luminosities do not violate the nondetection of the dust continuum in both cases. Together with three other z ≳ 11–13 candidate galaxies recently observed with ALMA, we conduct a joint ALMA and JWST spectral energy distribution (SED) analysis and find that the high- z solution at z ∼ 11–17 is favored in every candidate as a very blue (UV continuum slope of ≃−2.3) and luminous ( M _UV ≃ [ − 24:−21]) system. Still, we find in several candidates that reasonable SED fits (Δ χ ^2 ≲ 4) are reproduced by type II quasar and/or quiescent galaxy templates with strong emission lines at z ∼ 3–5, where such populations predicted from their luminosity functions and EW([O iii ]+H β ) distributions are abundant in survey volumes used for the identification of the z ∼ 11–17 candidates. While these recent ALMA observation results have strengthened the likelihood of the high- z solutions, lower- z possibilities are not completely ruled out in several of the z ∼ 11–17 candidates, indicating the need to consider the relative surface densities of the lower- z contaminants in the ultra-high- z galaxy search.

We present the JWST Resolved Stellar Populations Early Release Science (ERS) program. We obtained 27.5 hr of NIRCam and NIRISS imaging of three targets in the Local Group (Milky Way globular cluster M92, ultrafaint dwarf galaxy Draco II , and star-forming dwarf galaxy WLM), which span factors of ∼10 ^5 in luminosity, ∼10 ^4 in distance, and ∼10 ^5 in surface brightness. We describe the survey strategy, scientific and technical goals, implementation details, present select NIRCam color–magnitude diagrams (CMDs), and validate the NIRCam exposure time calculator (ETC). Our CMDs are among the deepest in existence for each class of target. They touch the theoretical hydrogen-burning limit in M92 (

We present the mid-infrared (MIR) morphologies for 64 star-forming galaxies (SFGs) at 0.2 < z < 2.5 with stellar mass M _* > 10 ^9 M _⊙ using James Webb Space Telescope (JWST) Mid-Infrared Instrument (MIRI) observations from the Cosmic Evolution Early Release Science survey. The MIRI bands span the MIR (7.7–21 μ m), enabling us to measure the effective radii ( R _eff ) and Sérsic indexes of these SFGs at rest-frame 6.2 and 7.7 μ m, which contains strong emission from Polycyclic aromatic hydrocarbon (PAH) features, a well-established tracer of star formation in galaxies. We define a “PAH band” as the MIRI bandpass that contains these features at the redshift of the galaxy. We then compare the galaxy morphologies in the PAH bands to those in the rest-frame near-ultraviolet (NUV) using Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS)/F435W or ACS/F606W and optical/near-IR using HST WFC3/F160W imaging from UVCANDELS and CANDELS. The R _eff of galaxies in the PAH band are slightly smaller (∼10%) than those in F160W for galaxies with M _* ≳ 10 ^9.5 M _⊙ at z ≤ 1.2, but the PAH band and F160W have similar fractions of light within 1 kpc. In contrast, the R _eff of galaxies in the NUV band are larger, with lower fractions of light within 1 kpc compared to F160W for galaxies at z ≤ 1.2. Using the MIRI data to estimate the SFR _IR surface density, we find that the correlation between the SFR _IR surface density and stellar mass has a steeper slope than that of the SFR _UV surface density and stellar mass, suggesting more massive galaxies having increasing amounts of obscured fraction of star formation in their inner regions. This paper demonstrates how the high-angular resolution data from JWST/MIRI can reveal new information about the morphology of obscured star formation.

We use the Illustris-1 simulation to explore the capabilities of the Hubble Space Telescope (HST) and James Webb Space Telescope (JWST) data to analyze the stellar populations in high-redshift galaxies, taking advantage of the combined depth, spatial resolution, and wavelength coverage. For that purpose, we use simulated broadband ACS, WFC3, and NIRCam data and two-dimensional stellar population synthesis (2D-SPS) to derive the integrated star formation history (SFH) of massive ( M _* > 10 ^10 M _⊙ ) simulated galaxies at 1 < z < 4 that evolve into a local M _* > 10 ^11 M _⊙ galaxy. In particular, we explore the potential of HST and JWST data sets reaching a depth similar to those of the CANDELS and ongoing CEERS observations, respectively, and concentrate on determining the capabilities of this data set for characterizing the first episodes in the SFH of local M _* > 10 ^11 M _⊙ galaxies by studying their progenitors at z > 1. The 2D-SPS method presented in this paper has been calibrated to robustly recover the cosmic times when the first star formation episodes occurred in massive galaxies, i.e., the first stages in their integrated SFHs. In particular, we discuss the times when the first 1%–50% of their total stellar mass formed in the simulation. We demonstrate that we can recover these ages with typical median systematic offset of less than 5% and scatter around 20%–30%. According to our measurements on Illustris data, we are able to recover that local M _* > 10 ^11 M _⊙ galaxies would have started their formation by z = 16, forming the first 5% of their stellar mass present at z ∼ 1 by z = 4.5, 10% by z = 3.7, and 25% by z = 2.7.

We present the first comprehensive release of photometric redshifts (photo - z 's) from the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS) team. We use statistics based upon the Quantile–Quantile ( Q – Q ) plot to identify biases and signatures of underestimated or overestimated errors in photo - z probability density functions (PDFs) produced by six groups in the collaboration; correcting for these effects makes the resulting PDFs better match the statistical definition of a PDF. After correcting each group’s PDF, we explore three methods of combining the different groups’ PDFs for a given object into a consensus curve. Two of these methods are based on identifying the minimum f -divergence curve, i.e., the PDF that is closest in aggregate to the other PDFs in a set (analogous to the median of an array of numbers). We demonstrate that these techniques yield improved results using sets of spectroscopic redshifts independent of those used to optimize PDF modifications. The best photo - z PDFs and point estimates are achieved with the minimum f -divergence using the best four PDFs for each object (mFDa4) and the hierarchical Bayesian (HB4) methods, respectively. The HB4 photo - z point estimates produced σ _NMAD = 0.0227/0.0189 and ∣Δ z /(1 + z )∣ > 0.15 outlier fraction = 0.067/0.019 for spectroscopic and 3D Hubble Space Telescope redshifts, respectively. Finally, we describe the structure and provide guidance for the use of the CANDELS photo - z catalogs, which are available at https://archive.stsci.edu/prepds/candels/ .

Identifying merging galaxies is an important—but difficult—step in galaxy evolution studies. We present random forest (RF) classifications of galaxy mergers from simulated JWST images based on various standard morphological parameters. We describe (a) constructing the simulated images from IllustrisTNG and the Santa Cruz SAM and modifying them to mimic future CEERS observations and nearly noiseless observations, (b) measuring morphological parameters from these images, and (c) constructing and training the RFs using the merger history information for the simulated galaxies available from IllustrisTNG. The RFs correctly classify ∼60% of non-merging and merging galaxies across 0.5 < z < 4.0. Rest-frame asymmetry parameters appear more important for lower-redshift merger classifications, while rest-frame bulge and clump parameters appear more important for higher-redshift classifications. Adjusting the classification probability threshold does not improve the performance of the forests. Finally, the shape and slope of the resulting merger fraction and merger rate derived from the RF classifications match with theoretical Illustris predictions but are underestimated by a factor of ∼0.5.

We present rest-frame optical emission-line flux ratio measurements for five z > 5 galaxies observed by the James Webb Space Telescope Near-Infared Spectrograph (NIRSpec) in the SMACS 0723 Early Release Observations. We add several quality-control and post-processing steps to the NIRSpec pipeline reduction products in order to ensure reliable relative flux calibration of emission lines that are closely separated in wavelength, despite the uncertain absolute spectrophotometry of the current version of the reductions. Compared to z ∼ 3 galaxies in the literature, the z > 5 galaxies have similar [O iii ] λ 5008/H β ratios, similar [O iii ] λ 4364/H γ ratios, and higher (∼0.5 dex) [Ne III ] λ 3870/[O II ] λ 3728 ratios. We compare the observations to MAPPINGS V photoionization models and find that the measured [Ne III ] λ 3870/[O II ] λ 3728, [O iii ] λ 4364/H γ , and [O iii ] λ 5008/H β emission-line ratios are consistent with an interstellar medium (ISM) that has very high ionization ( $\mathrm{log}(Q)\simeq 8-9$ , units of cm s ^−1 ), low metallicity ( Z / Z _⊙ ≲ 0.2), and very high pressure ( $\mathrm{log}(P/k)\simeq 8-9$ , units of cm ^−3 ). The combination of [O iii ] λ 4364/H γ and [O iii ] λ (4960 + 5008)/H β line ratios indicate very high electron temperatures of $4.1\lt \mathrm{log}({T}_{e}/{\rm{K}})\lt 4.4$ , further implying metallicities of Z / Z _⊙ ≲ 0.2 with the application of low-redshift calibrations for “ T _e -based” metallicities. These observations represent a tantalizing new view of the physical conditions of the ISM in galaxies at cosmic dawn.

We investigate the resolved kinematics of the molecular gas, as traced by the Atacama Large Millimeter/submillimeter Array in CO (2−1), of 25 cluster member galaxies across three different clusters at a redshift of z ∼ 1.6. This is the first large-scale analysis of the molecular gas kinematics of cluster galaxies at this redshift. By separately estimating the rotation curve of the approaching and receding sides of each galaxy via kinematic modeling, we quantify the difference in total circular velocity to characterize the overall kinematic asymmetry of each galaxy. 3/14 of the galaxies in our sample that we are able to model have similar degrees of asymmetry as that observed in galaxies in the field at similar redshift based on observations of mainly ionized gas. However, this leaves 11/14 galaxies in our sample with significantly higher asymmetry, and some of these galaxies have degrees of asymmetry of up to ∼50 times higher than field galaxies observed at similar redshift. Some of these extreme cases also have one-sided tail-like morphology seen in the molecular gas, supporting a scenario of tidal and/or ram pressure interaction. Such stark differences in the kinematic asymmetry in clusters versus the field suggest the evolutionary influence of dense environments, established as being a major driver of galaxy evolution at low redshift, is also active in the high-redshift universe.