Search

We use JWST/NIRISS slitless spectroscopy from the Next Generation Deep Extragalactic Exploratory Public (NGDEEP) Survey to investigate the physical condition of 178 star-forming galaxies at 1.7

Over the past decades, a population of galaxies invisible in optical/near-infrared (NIR), but bright at longer wavelengths, have been identified through color selections. These so-called optically faint/dark galaxies are considered to be massive quiescent galaxies or highly dust-attenuated galaxies. Having the entire galaxy obscured by dust, however, is likely an extreme case of the much more common occurrence of optically thin and thick absorption coexisting in the same system. With the power of JWST imaging, we are able to spatially resolve massive galaxies at z ∼ 3, accurately model their spectral energy distributions, and identify candidate optically thick substructures. We target galaxies with $\mathrm{log}({M}_{* }/{M}_{\odot })\gt $ 10.3 and 2.5

A wealth of observations have long suggested that the vast majority of isolated classical dwarf galaxies ( M _* = 10 ^7 –10 ^9 M _⊙ ) are currently star forming. However, recent observations of the large abundance of “ultra-diffuse galaxies” beyond the reach of previous large spectroscopic surveys suggest that our understanding of the dwarf galaxy population may be incomplete. Here we report the serendipitous discovery of an isolated quiescent dwarf galaxy in the nearby Universe, which was imaged as part of the JWST PEARLS Guaranteed Time Observation program. Remarkably, individual red-giant branch stars are visible in this near-IR imaging, suggesting a distance of 30 ± 4 Mpc, and a wealth of archival photometry point to an sSFR of 2 × 10 ^−11 yr ^−1 and star formation rate of 4 × 10 ^−4 M _⊙ yr ^−1 . Spectra obtained with the Lowell Discovery Telescope find a recessional velocity consistent with the Hubble Flow and >1500 km s ^−1 separated from the nearest massive galaxy in Sloan Digital Sky Survey suggesting that this galaxy was either quenched from internal mechanisms or had a very high-velocity (≳1000 km s ^−1 ) interaction with a nearby massive galaxy in the past. This analysis highlights the possibility that many nearby quiescent dwarf galaxies are waiting to be discovered and that JWST has the potential to resolve them.

We study the H α equivalent width (EW(H α )) maps of 19 galaxies at 0.6 < z < 2.2 in the Hubble Ultra Deep Field using NIRISS slitless spectroscopy as part of the Next Generation Deep Extragalactic Exploratory Public Survey. Our galaxies mostly lie on the star formation main sequence with stellar masses between 10 ^9 and 10 ^11 M _⊙ , characterized as “typical” star-forming galaxies at these redshifts. Leveraging deep Hubble Space Telescope and JWST images, spanning 0.4–4.8 μ m, we perform spatially resolved fitting of the spectral energy distributions for these galaxies and construct specific star formation rate (sSFR) and stellar-mass-weighted age maps with a spatial resolution of ∼1 kpc. The pixel-to-pixel EW(H α ) increases with increasing sSFR and with decreasing age. The average trends are slightly different from the relations derived from integrated fluxes of galaxies from the literature, suggesting complex evolutionary trends within galaxies. We quantify the radial profiles of EW(H α ), sSFR, and age. The majority (84%) of galaxies show positive EW(H α ) gradients, in line with the inside-out quenching scenario. A few galaxies (16%) show inverse (and flat) EW(H α ) gradients, possibly due to merging or starbursts. We compare the distributions of EW(H α ) and sSFR to star formation history (SFH) models as a function of galactocentric radius. We argue that the central regions of galaxies have experienced at least one rapid star formation episode, which leads to the formation of the bulge, while their outer regions (e.g., disks) grow via more smoothly varying SFHs. These results demonstrate the ability to study resolved star formation in distant galaxies with JWST NIRISS.

We present measurements of the rest-frame UV spectral slope, β , for a sample of 36 faint star-forming galaxies at z ∼ 9–16 discovered in one of the deepest JWST NIRCam surveys to date, the Next Generation Deep Extragalactic Exploratory Public Survey. We use robust photometric measurements for UV-faint galaxies (down to M _UV ∼ −16), originally published in Leung et al., and measure values of the UV spectral slope via photometric power-law fitting to both the observed photometry and stellar population models obtained through spectral energy distribution (SED) fitting with Bagpipes . We obtain a median and 68% confidence interval for β from photometric power-law fitting of ${\beta }_{\mathrm{PL}}=-{2.7}_{-0.5}^{+0.5}$ and from SED fitting, ${\beta }_{\mathrm{SED}}=-{2.3}_{-0.1}^{+0.2}$ for the full sample. We show that when only two to three photometric detections are available, SED fitting has a lower scatter and reduced biases than photometric power-law fitting. We quantify this bias and find that after correction the median ${\beta }_{\mathrm{SED},\mathrm{corr}}=-{2.5}_{-0.2}^{+0.2}$ . We measure physical properties for our galaxies with Bagpipes and find that our faint ( ${M}_{\mathrm{UV}}=-{18.1}_{-0.9}^{+0.7}$ ) sample is low in mass ( $\mathrm{log}[{M}_{* }/{M}_{\odot }]={7.7}_{-0.5}^{+0.5}$ ), fairly dust-poor ( ${A}_{{\rm{v}}}={0.1}_{-0.1}^{+0.2}$ mag), and modestly young ( $\mathrm{log}[\mathrm{age}]={7.8}_{-0.8}^{+0.2}$ yr) with a median star formation rate of $\mathrm{log}(\mathrm{SFR})=-{0.3}_{-0.4}^{+0.4}{M}_{\odot }\,{\mathrm{yr}}^{-1}$ . We find no strong evidence for ultrablue UV spectral slopes ( β ∼ −3) within our sample, as would be expected for exotically metal-poor ( Z / Z _⊙ < 10 ^−3 ) stellar populations with very high Lyman continuum escape fractions. Our observations are consistent with model predictions that galaxies of these stellar masses at z ∼ 9–16 should have only modestly low metallicities ( Z / Z _⊙ ∼ 0.1–0.2).

James Webb Space Telescope NIRCam images have revealed 154 reliable globular cluster (GC) candidates around the z = 0.0513 elliptical galaxy VV 191a after subtracting 34 likely interlopers from background galaxies inside our search area. NIRCam broadband observations are made at 0.9–4.5 μ m using the F090W, F150W, F356W, and F444W filters. Using point-spread-function-matched photometry, the data are analyzed to present color–magnitude diagrams and color distributions that suggest a relatively uniform population of GCs, except for small fractions of reddest (5%–8%) and bluest (2%–4%) outliers. GC models in the F090W versus (F090W–F150W) diagram fit the NIRCam data well and show that the majority of GCs detected have a mass of ∼10 ^6.5 M _⊙ , with metallicities [Fe/H] spanning the typical range expected for GCs (−2.5 ≲ [Fe/H]≲ 0.5). However, the models predict ∼0.3–0.4 mag bluer (F356W–F444W) colors than the NIRCam data for a reasonable range of GC ages, metallicities, and reddening. Although our data do not quite reach the luminosity function turnover, the measured luminosity function is consistent with previous measurements, suggesting an estimated peak at m _AB ∼ −9.4 ± 0.2 mag in the F090W filter.

We present the Next Generation Deep Extragalactic Exploratory Public (NGDEEP) Survey, a deep slitless spectroscopic and imaging Cycle 1 JWST treasury survey designed to constrain feedback mechanisms in low-mass galaxies across cosmic time. NGDEEP targets the Hubble Ultra Deep Field (HUDF) with NIRISS slitless spectroscopy ( ${f}_{\mathrm{lim},\mathrm{line},5\sigma }\approx 1.2\,\times \,$ 10 ^−18 erg s ^−1 cm ^−2 ) to measure metallicities and star formation rates (SFRs) for low-mass galaxies through the peak of the cosmic SFR density (0.5 < z < 4). In parallel, NGDEEP targets the HUDF-Par2 parallel field with NIRCam ( ${m}_{\mathrm{lim},5\sigma }=30.6-30.9$ ) to discover galaxies to z > 12, constraining the slope of the faint end of the rest-ultraviolet luminosity function. NGDEEP overlaps with the deepest HST Advanced Camera for Surveys optical imaging in the sky, F435W in the HUDF ( ${m}_{\mathrm{lim},{\rm{F}}435{\rm{W}}}=29.6$ ) and F814W in HUDF-Par2 ( ${m}_{\mathrm{lim},{\rm{F}}814{\rm{W}}}=30$ ), making this a premier HST+JWST deep field. As a treasury survey, NGDEEP data are public immediately, and we will rapidly release data products and catalogs in the spirit of previous deep-field initiatives. In this paper we present the NGDEEP survey design, summarize the science goals, and detail plans for the public release of NGDEEP reduced data products.

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.

A crucial yet challenging task in galaxy evolution studies is the identification of distant merging galaxies, a task that suffers from a variety of issues ranging from telescope sensitivities and limitations to the inherently chaotic morphologies of young galaxies. In this paper, we use random forests and convolutional neural networks to identify high-redshift JWST Cosmic Evolution Early Release Science Survey (CEERS) galaxy mergers. We train these algorithms on simulated 3 < z < 5 CEERS galaxies created from the IllustrisTNG subhalo morphologies and the Santa Cruz SAM light cone. We apply our models to observed CEERS galaxies at 3 < z < 5. We find that our models correctly classify ∼60%–70% of simulated merging and nonmerging galaxies; better performance on the merger class comes at the expense of misclassifying more nonmergers. We could achieve more accurate classifications, as well as test for a dependency on physical parameters such as gas fraction, mass ratio, and relative orbits, by curating larger training sets. When applied to real CEERS galaxies using visual classifications as ground truth, the random forests correctly classified 40%–60% of mergers and nonmergers at 3 < z < 4 but tended to classify most objects as nonmergers at 4 < z < 5 (misclassifying ∼70% of visually classified mergers). On the other hand, the CNNs tended to classify most objects as mergers across all redshifts (misclassifying 80%–90% of visually classified nonmergers). We investigate what features the models find most useful, as well as the characteristics of false positives and false negatives, and also calculate merger rates derived from the identifications made by the models.

A Type Ia supernova (SN) at z = 1.78 was discovered in James Webb Space Telescope Near Infrared Camera imaging of the galaxy cluster PLCK G165.7+67.0 (G165; z = 0.35). The SN is situated 1.5–2 kpc from the host-galaxy nucleus and appears in three different locations as a result of gravitational lensing by G165. These data can yield a value for Hubble’s constant using time delays from this multiply imaged SN Ia that we call “SN H0pe.” Over the cluster, we identified 21 image multiplicities, confirmed five of them using the Near-Infrared Spectrograph, and constructed a new lens model that gives a total mass within 600 kpc of (2.6 ± 0.3) × 10 ^14 M _⊙ . The photometry uncovered a galaxy overdensity coincident with the SN host galaxy. NIRSpec confirmed six member galaxies, four of which surround the SN host galaxy with relative velocity ≲900 km s ^−1 and projected physical extent ≲33 kpc. This compact galaxy group is dominated by the SN host galaxy, which has a stellar mass of (5.0 ± 0.1) × 10 ^11 M _⊙ . The group members have specific star formation rates of 2–260 Gyr ^−1 derived from the H α -line fluxes corrected for stellar absorption, dust extinction, and slit losses. Another group centered on a strongly lensed dusty star-forming galaxy is at z = 2.24. The total (unobscured and obscured) SFR of this second galaxy group is estimated to be (≳ 100 M _⊙ yr ^−1 ), which translates to a supernova rate of ∼1 SNe yr ^−1 , suggesting that regular monitoring of this cluster may yield additional SNe.

Visual inspections of the first optical rest-frame images from JWST have indicated a surprisingly high fraction of disk galaxies at high redshifts. Here, we alternatively apply self-supervised machine learning to explore the morphological diversity at z ≥ 3. Our proposed data-driven representation scheme of galaxy morphologies, calibrated on mock images from the TNG50 simulation, is shown to be robust to noise and to correlate well with the physical properties of the simulated galaxies, including their 3D structure. We apply the method simultaneously to F200W and F356W galaxy images of a mass-complete sample ( M _* / M _⊙ > 10 ^9 ) at 3 ≤ z ≤ 6 from the first JWST/NIRCam CEERS data release. We find that the simulated and observed galaxies do not exactly populate the same manifold in the representation space from contrastive learning. We also find that half the galaxies classified as disks—either convolutional neural network-based or visually—populate a similar region of the representation space as TNG50 galaxies with low stellar specific angular momentum and nonoblate structure. Although our data-driven study does not allow us to firmly conclude on the true nature of these galaxies, it suggests that the disk fraction at z ≥ 3 remains uncertain and possibly overestimated by traditional supervised classifications. Deeper imaging and spectroscopic follow-ups as well as comparisons with other simulations will help to unambiguously determine the true nature of these galaxies, and establish more robust constraints on the emergence of disks at very high redshift.

We combine deep imaging data from the CEERS early release JWST survey and Hubble Space Telescope imaging from CANDELS to examine the size–mass relation of star-forming galaxies and the morphology–quenching relation at stellar masses M _⋆ ≥ 10 ^9.5 M _⊙ over the redshift range 0.5 < z < 5.5. In this study with a sample of 2450 galaxies, we separate star-forming and quiescent galaxies based on their star formation activity and confirm that star-forming and quiescent galaxies have different morphologies out to z = 5.5, extending the results of earlier studies out to higher redshifts. We find that star-forming and quiescent galaxies have typical Sérsic indices of n ∼ 1.3 and n ∼ 4.3, respectively. Focusing on star-forming galaxies, we find that the slope of the size–mass relation is nearly constant with redshift, as was found previously, but shows a modest increase at z ∼ 4.2. The intercept in the size–mass relation declines out to z = 5.5 at rates that are similar to what earlier studies found. The intrinsic scatter in the size–mass relation is relatively constant out to z = 5.5.

The 3D geometries of high-redshift galaxies remain poorly understood. We build a differentiable Bayesian model and use Hamiltonian Monte Carlo to efficiently and robustly infer the 3D shapes of star-forming galaxies in James Webb Space Telescope Cosmic Evolution Early Release Science observations with $\mathrm{log}{M}_{* }/{M}_{\odot }=9.0\mbox{--}10.5$ at z = 0.5–8.0. We reproduce previous results from the Hubble Space Telescope Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey in a fraction of the computing time and constrain the mean ellipticity, triaxiality, size, and covariances with samples as small as ∼50 galaxies. We find high 3D ellipticities for all mass–redshift bins, suggesting oblate (disky) or prolate (elongated) geometries. We break that degeneracy by constraining the mean triaxiality to be ∼1 for $\mathrm{log}{M}_{* }/{M}_{\odot }=9.0\mbox{--}9.5$ dwarfs at z > 1 (favoring the prolate scenario), with significantly lower triaxialities for higher masses and lower redshifts indicating the emergence of disks. The prolate population traces out a “banana” in the projected $b/a\mbox{--}\mathrm{log}a$ diagram with an excess of low- b / a , large- $\mathrm{log}a$ galaxies. The dwarf prolate fraction rises from ∼25% at z = 0.5–1.0 to ∼50%–80% at z = 3–8. Our results imply a second kind of disk settling from oval (triaxial) to more circular (axisymmetric) shapes with time. We simultaneously constrain the 3D size–mass relation and its dependence on 3D geometry. High-probability prolate and oblate candidates show remarkably similar Sérsic indices ( n ∼ 1), nonparametric morphological properties, and specific star formation rates. Both tend to be visually classified as disks or irregular, but edge-on oblate candidates show more dust attenuation. We discuss selection effects, follow-up prospects, and theoretical implications.

We study a new population of extremely red objects (EROs) recently discovered by the James Webb Space Telescope (JWST) based on their NIRCam colors F277W − F444W > 1.5 mag. We find 37 EROs in the Cosmic Evolution Early Release Science Survey (CEERS) field with F444W < 28 mag and photometric redshifts between 5 < z < 7, with median $z={6.9}_{-1.6}^{+1.0}$ . Surprisingly, despite their red long-wavelength colors, these EROs have blue short-wavelength colors (F150W − F200W ∼ 0 mag) indicative of bimodal spectral energy distributions (SEDs) with a red, steep slope in the rest-frame optical, and a blue, flat slope in the rest-frame UV. Moreover, all these EROs are unresolved, point-like sources in all NIRCam bands. We analyze the SEDs of eight of them with MIRI and NIRSpec observations using stellar population models and active galactic nucleus (AGN) templates. We find that dusty galaxies or obscured AGNs provide similarly good SED fits but different stellar properties: massive and dusty, $\mathrm{log}{M}_{\star }/{M}_{\odot }$ ∼ 10 and A _V ≳ 3 mag, or low mass and obscured, $\mathrm{log}{M}_{\star }/{M}_{\odot }$ ∼ 7.5 and A _V ∼ 0 mag, hosting an obscured quasi-stellar object (QSO). SED modeling does not favor either scenario, but their unresolved sizes are more suggestive of AGNs. If any EROs are confirmed to have $\mathrm{log}{M}_{\star }/{M}_{\odot }$ ≳ 10.5, it would increase the pre-JWST number density at z > 7 by up to a factor ∼60. Similarly, if they are QSOs with luminosities in the L _bol > 10 ^45–46 erg s ^−1 range, their number would exceed that of bright blue QSOs by more than three orders of magnitude. Additional photometry at mid-infrared wavelengths will reveal the true nature of the red continuum emission in these EROs and will place this puzzling population in the right context of galaxy evolution.

We present an analysis of the ultraviolet luminosity function (UV LF) and star formation rate density of distant galaxies (7.5 < z < 13.5) in the “blank” fields of the Prime Extragalactic Areas for Reionization and Lensing Science (PEARLS) survey combined with Early Release Science data from the CEERS, GLASS, and NGDEEP surveys/fields and the first data release of JADES. We use strict quality cuts on EAZY photometric redshifts to obtain a reliable selection and characterization of high-redshift ( z > 6.5) galaxies from a consistently processed set of deep, near-infrared imaging. Within an area of 180 arcmin ^2 , we identify 1046 candidate galaxies at redshifts z > 6.5 and we use this sample to study the UV LF in four redshift bins between 7.5 < z < 13.5. The measured number density of galaxies at z = 8 and z = 9 matches those of past observations undertaken by the Hubble Space Telescope (HST). Our z = 10.5 measurements lie between early James Webb Space Telescope (JWST) results and past HST results, indicating cosmic variance may be the cause of previous high density measurements. However, the number densities of UV-luminous galaxies at z = 12.5 are high compared to predictions from simulations. When examining the star formation rate density of galaxies at this period, our observations are still largely consistent with a constant star formation efficiency, are slightly lower than previous early estimations using JWST, and support galaxy driven reionization at z ≤ 8.

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.

The Next Generation Deep Extragalactic Exploratory Public (NGDEEP) survey program was designed specifically to include Near Infrared Slitless Spectroscopic observations (NGDEEP-NISS) to detect multiple emission lines in as many galaxies as possible and across a wide redshift range using the Near Infrared Imager and Slitless Spectrograph. We present early results obtained from the first set of observations (Epoch 1, 50% of the allocated orbits) of this program (NGDEEP-NISS1). Using a set of independently developed calibration files designed to deal with a complex combination of overlapping spectra, multiple position angles, and multiple cross filters and grisms, in conjunction with a robust and proven algorithm for quantifying contamination from overlapping dispersed spectra, NGDEEP-NISS1 has achieved a 3 σ sensitivity limit of 2 × 10 ^−18 erg s ^−1 cm ^−2 . We demonstrate the power of deep wide field slitless spectroscopy (WFSS) to characterize the star formation rates, and metallicity ([O iii ]/H β ), and dust content, of galaxies at 1 < z < 3.5. The latter showing intriguing initial results on the applicability and assumptions made regarding the use of Case B recombination. Further, we identify the presence of active galactic nuclei and infer the mass of their supermassive black holes using broadened restframe Mg ii and H β emission lines. The spectroscopic results are then compared with the physical properties of galaxies extrapolated from fitting spectral energy distribution models to photometry alone. The results clearly demonstrate the unique power and efficiency of WFSS at near-infrared wavelengths over other methods to determine the properties of galaxies across a broad range of redshifts.

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 use James Webb Space Telescope Near-Infrared Camera Wide Field Slitless Spectroscopy (NIRCam WFSS) and the Near-Infrared spectrograph (NIRSpec) in the Cosmic Evolution Early Release survey to measure rest-frame optical emission-line ratios of 155 galaxies at z > 2. The blind NIRCam grism observations include a sample of galaxies with bright emission lines that were not observed on the NIRSpec masks. We study the changes of the H α , [O III ]/H β , and [Ne III ]/[O II ] emission lines in terms of redshift by comparing to lower-redshift SDSS, CLEAR, and MOSDEF samples. We find a significant (>3 σ ) correlation between [O III ]/H β with redshift, while [Ne III ]/[O II ] has a marginal (2 σ ) correlation with redshift. We compare [O III ]/H β and [Ne III ]/[O II ] to stellar mass and H β SFR. We find that both emission-line ratios have a correlation with H β SFR and an anticorrelation with stellar mass across the redshifts 0 < z < 9. Comparison with MAPPINGS V models indicates that these trends are consistent with lower metallicity and higher ionization in low-mass and high-SFR galaxies. We additionally compare to IllustrisTNG predictions and find that they effectively describe the highest [O III ]/H β ratios observed in our sample, without the need to invoke MAPPINGS models with significant shock ionization components.

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 a new parametric lens model for the G165.7+67.0 galaxy cluster, which was discovered with Planck through its bright submillimeter flux, originating from a pair of extraordinary dusty star-forming galaxies (DSFGs) at z ≈ 2.2. Using JWST and interferometric mm/radio observations, we characterize the intrinsic physical properties of the DSFGs, which are separated by only ∼1″ (8 kpc) and a velocity difference Δ V ≲ 600 km s ^−1 in the source plane, and thus are likely undergoing a major merger. Boasting intrinsic star formation rates SFR _IR = 320 ± 70 and 400 ± 80 M _⊙ yr ^−1 , stellar masses of $\mathrm{log}[{M}_{\star }/{M}_{\odot }]=10.2\pm 0.1$ and 10.3 ± 0.1, and dust attenuations of A _V = 1.5 ± 0.3 and 1.2 ± 0.3, they are remarkably similar objects. We perform spatially resolved pixel-by-pixel spectral energy distribution (SED) fitting using rest-frame near-UV to near-IR imaging from JWST/NIRCam for both galaxies, resolving some stellar structures down to 100 pc scales. Based on their resolved specific star formation rates (SFRs) and UVJ colors, both DSFGs are experiencing significant galaxy-scale star formation events. If they are indeed interacting gravitationally, this strong starburst could be the hallmark of gas that has been disrupted by an initial close passage. In contrast, the host galaxy of SN H0pe has a much lower SFR than the DSFGs, and we present evidence for the onset of inside-out quenching and large column densities of dust even in regions of low specific SFR. Based on the intrinsic SFRs of the DSFGs inferred from UV through far-infrared SED modeling, this pair of objects alone is predicted to yield an observable 1.1 ± 0.2 core-collapse supernovae per year, making this cluster field ripe for continued monitoring.

We present a sample of 1165 extreme emission-line galaxies (EELGs) at 4 < z < 9 selected using James Webb Space Telescope (JWST) NIRCam photometry in the Cosmic Evolution Early Release Science (CEERS) program. We use a simple method to photometrically identify EELGs with H β + [O iii ] (combined) or H α emission of observed-frame equivalent width (EW) > 5000 Å. JWST/NIRSpec spectroscopic observations of a subset (34) of the photometrically selected EELGs validate our selection method: All spectroscopically observed EELGs confirm our photometric identification of extreme emission, including some cases where the spectral-energy-distribution-derived photometric redshifts are incorrect. We find that the medium-band F410M filter in CEERS is particularly efficient at identifying EELGs, both in terms of including emission lines in the filter and in correctly identifying the continuum between H β + [O iii ] and H α in the neighboring broadband filters. We present examples of EELGs that could be incorrectly classified as ultrahigh redshift ( z > 12) as a result of extreme H β + [O iii ] emission blended across the reddest photometric filters. We compare the EELGs to the broader (subextreme) galaxy population in the same redshift range and find that they are consistent with being the bluer, high-EW tail of a broader population of emission-line galaxies. The highest-EW EELGs tend to have more compact emission-line sizes than continuum sizes, suggesting that active galactic nuclei are responsible for at least some of the most extreme EELGs. The photometrically inferred emission-line ratios are consistent with interstellar medium conditions with high ionization and moderately low metallicity, consistent with previous spectroscopic studies.

The first public 0.9–4.4 μ m NIRCam images of the North Ecliptic Pole Time Domain Field uncovered galaxies displaying point-source features in their cores as seen in the longer-wavelength filters. We visually identified a sample of 66 galaxies (∼1 galaxy arcmin ^–2 ) with pointlike cores and have modeled their two-dimensional light profiles with GalFit , identifying 16 galactic nuclei with measurable point-source components. GalFit suggests that the visual sample is a mix of both compact stellar bulge and point-source galaxy cores. This core classification is complemented by spectral energy distribution modeling to infer the sample’s active galactic nucleus (AGN) and host-galaxy parameters. For galaxies with measurable point-source components, the median fractional AGN contribution to their 0.1–30.0 μ m flux is 0.44, and 14/16 are color-classified AGN. We conclude that near-infrared point-source galaxy cores are signatures of AGN. In addition, we define an automated sample-selection criterion to identify these point-source features. This criterion can be used in other extant and future NIRCam images to streamline the search for galaxies with unresolved IR-luminous AGN. The James Webb Space Telescope’s superb angular resolution and sensitivity at infrared wavelengths are resurrecting the morphological identification of AGN.

New JWST/NIRCam wide-field slitless spectroscopy provides redshifts for four z > 8 galaxies located behind the lensing cluster MACS J0416.1−2403. Two of them, “Y1” and “JD,” have previously reported spectroscopic redshifts based on Atacama Large Millimeter/submillimeter Array measurements of [O iii ] 88 μ m and/or [C ii ] 157.7 μ m lines. Y1 is a merging system of three components, and the existing redshift z = 8.31 is confirmed. However, JD is at z = 8.34 instead of the previously claimed z = 9.28. JD’s close companion, “JD-N,” which was a previously discovered z > 8 candidate, is now identified at the same redshift as JD. JD and JD-N form an interacting pair. A new candidate at z > 8, “f090d_018,” is also confirmed and is at z = 8.49. These four objects are likely part of an overdensity that signposts a large structure extending ∼165 kpc in projected distance and ∼48.7 Mpc in radial distance. They are magnified by less than 1 mag and have an intrinsic M _UV ranging from −19.57 to −20.83 mag. Their spectral energy distributions show that the galaxies are all very young with ages ∼ 4–18 Myr and stellar masses of about 10 ^7–8 M _⊙ . These infant galaxies have very different star formation rates ranging from a few to over a hundred solar masses per year, but only two of them (JD and f090d_018) have blue rest-frame UV slopes β < −2.0 indicative of a high Lyman-continuum photon escape fraction that could contribute significantly to the cosmic hydrogen-reionizing background. Interestingly, these two galaxies are the least massive and least active ones among the four. The other two systems have much flatter UV slopes largely because of their high dust extinction ( A _V = 0.9–1.0 mag). Their much lower indicated escape fractions show that even very young, actively star-forming galaxies can have a negligible contribution to reionization when they quickly form dust throughout their bodies.

We report the results of James Webb Space Telescope/NIRCam observations of 19 (sub)millimeter sources detected by the Atacama Large Millimeter Array (ALMA). The accurate ALMA positions allowed unambiguous identifications of their NIRCam counterparts. Taking gravitational lensing into account, these represent 16 distinct galaxies in three fields and constitute the largest sample of its kind to date. The counterparts’ spectral energy distributions cover from rest-frame ultraviolet to near-IR and provide photometric redshifts (1 < z < 4.5) and stellar masses ( M _* > 10 ^10.5 M _⊙ ), which are similar to submillimeter galaxies (SMGs) studied previously. However, our sample is fainter in (sub)millimeter than the classic SMG samples are, and our sources exhibit a wider range of properties. They have dust-embedded star formation rates as low as 10 M _⊙ yr ^−1 , and the sources populate both the star-forming main sequence and the quiescent categories. The deep NIRCam data allow us to study the rest-frame near-IR morphologies. Excluding two multiply imaged systems and one quasar, the majority of the remaining sources are disk-like and show either little or no disturbance. This suggests that secular growth is a potential route for the assembly of high-mass disk galaxies. While a few objects have large disks, the majority have small disks (median half-mass radius of 1.6 kpc). At this time, it is unclear whether this is due to the prevalence of small disks at these redshifts or some unknown selection effects of deep ALMA observations. A larger sample of ALMA sources with NIRCam observations will be able to address this question.

We present JWST/NIRSpec prism spectroscopy of seven galaxies selected from Cosmic Evolution Early Release Science (CEERS) survey NIRCam imaging with photometric redshifts z _phot > 8. We measure emission line redshifts of z = 7.65 and 8.64 for two galaxies. For two other sources without securely detected emission lines we measure $z={9.77}_{-0.29}^{+0.37}$ and ${10.01}_{-0.19}^{+0.14}$ by fitting model spectral templates to the prism data, from which we detect continuum breaks consistent with Ly α opacity from a mostly neutral intergalactic medium. The presence of strong breaks and the absence of strong emission lines give high confidence that these two galaxies have redshifts z > 9.6, but the redshift values derived from the breaks alone have large uncertainties given the low spectral resolution and relatively low S/N of the CEERS NIRSpec prism data. The two z ∼ 10 sources observed are relatively luminous ( M _UV < −20), with blue continua (−2.3 ≲ β ≲ −1.9) and low dust attenuation ( ${A}_{V}\simeq {0.15}_{-0.1}^{+0.3}$ ); and at least one of them has a high stellar mass for a galaxy at that redshift ( $\mathrm{log}({M}_{\star }/{M}_{\odot })\simeq {9.3}_{-0.3}^{+0.2}$ ). Considered together with spectroscopic observations of other CEERS NIRCam-selected high- z galaxy candidates in the literature, we find a high rate of redshift confirmation and low rate of confirmed interlopers (8%). Ten out of 35 z > 8 candidates with CEERS NIRSpec spectroscopy do not have secure redshifts, but the absence of emission lines in their spectra is consistent with redshifts z > 9.6. We find that z > 8 photometric redshifts are generally in agreement (within their uncertainties) with the spectroscopic values, but also that the photometric redshifts tend to be slightly overestimated (〈Δ z 〉 = 0.45 ± 0.11), suggesting that current templates do not fully describe the spectra of very-high- z sources. Overall, the spectroscopy solidifies photometric redshift evidence for a high spatial density of bright galaxies at z > 8 compared to theoretical model predictions, and further disfavors an accelerated decline in the integrated UV luminosity density at z > 8.

We present a robust sample of very high redshift galaxy candidates from the first epoch of JWST/NIRCam imaging from the Next Generation Deep Extragalactic Exploratory Public (NGDEEP) survey. The NGDEEP NIRCam imaging, spanning 9.7 arcmin ^2 in the Hubble Ultra Deep Field Parallel Field 2, reaches m = 30.4 (5 σ , point-source, 2″ diameter apertures corrected to total) in F277W, making it the deepest public JWST GO imaging data set to date. We describe our detailed data reduction process of the six-filter broadband JWST/NIRCam imaging, incorporating custom corrections for systematic effects to produce high-quality calibrated images. Using robust photometric redshift selection criteria, we identify a sample of 38 z ≳ 9 galaxy candidates. These objects span a redshift range of z = 8.5–15.8 and apparent magnitudes of m _F277W = 27–30.5 AB mag, reaching ∼1.5 mag deeper than previous public JWST imaging surveys. We calculate the rest-frame ultraviolet luminosity function at z ∼ 9 and 11 and present a new measurement of the luminosity function faint-end slope at z ∼ 11. We find a faint-end slope of α = −2.5 ± 0.4 and −2.2 ± 0.2 at z ∼ 9 and 11, respectively. This is consistent with no significant evolution in the faint-end slope and number density from z = 9 to 11. Comparing our results with theoretical predictions, we find that some models produce better agreement at the faint end than the bright end. These results will help to constrain how stellar feedback impacts star formation at these early epochs.

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.

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 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.

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.

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 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 the discovery of an accreting supermassive black hole at z = 8.679. This galaxy, denoted here as CEERS_1019, was previously discovered as a Ly α -break galaxy by Hubble with a Ly α redshift from Keck. As part of the Cosmic Evolution Early Release Science (CEERS) survey, we have observed this source with JWST/NIRSpec, MIRI, NIRCam, and NIRCam/WFSS and uncovered a plethora of emission lines. The H β line is best fit by a narrow plus a broad component, where the latter is measured at 2.5 σ with an FWHM ∼1200 km s ^−1 . We conclude this originates in the broadline region of an active galactic nucleus (AGN). This is supported by the presence of weak high-ionization lines (N V, N IV], and C III]), as well as a spatial point-source component. The implied mass of the black hole (BH) is log ( M _BH / M _⊙ ) = 6.95 ± 0.37, and we estimate that it is accreting at 1.2 ± 0.5 times the Eddington limit. The 1–8 μ m photometric spectral energy distribution shows a continuum dominated by starlight and constrains the host galaxy to be massive (log M/M _⊙ ∼9.5) and highly star-forming (star formation rate, or SFR ∼ 30 M _⊙ yr ^−1 ; log sSFR ∼ − 7.9 yr ^−1 ). The line ratios show that the gas is metal-poor ( Z / Z _⊙ ∼ 0.1), dense ( n _e ∼ 10 ^3 cm ^−3 ), and highly ionized (log U ∼ − 2.1). We use this present highest-redshift AGN discovery to place constraints on BH seeding models and find that a combination of either super-Eddington accretion from stellar seeds or Eddington accretion from very massive BH seeds is required to form this object.

The Cosmic Evolution Early Release Science program observed the Extended Groth Strip (EGS) with the Mid-Infrared Instrument (MIRI) on the James Webb Space Telescope (JWST) in 2022. In this paper, we discuss the four MIRI pointings that observed with longer-wavelength filters, including F770W, F1000W, F1280W, F1500W, F1800W, and F2100W. We compare the MIRI galaxies with the Spitzer/MIPS 24 μ m population in the EGS field. We find that MIRI can observe an order of magnitude deeper than MIPS in significantly shorter integration times, attributable to JWST's much larger aperture and MIRI’s improved sensitivity. MIRI is exceptionally good at finding faint ( L _IR < 10 ^10 L _⊙ ) galaxies at z ∼ 1–2. We find that a significant portion of MIRI galaxies are “mid-IR weak”—they have strong near-IR emission and relatively weaker mid-IR emission, and most of the star formation is unobscured. We present new IR templates that capture how the mid-to-near-IR emission changes with increasing infrared luminosity. We present two color–color diagrams to separate mid-IR weak galaxies and active galactic nuclei (AGN) from dusty star-forming galaxies and find that these color diagrams are most effective when used in conjunction with each other. We present the first number counts of 10 μ m sources and find that there are ≲10 IR AGN per MIRI pointing, possibly due to the difficulty of distinguishing AGN from intrinsically mid-IR weak galaxies (due to low metallicities or dust content). We conclude that MIRI is most effective at observing moderate-luminosity ( L _IR = 10 ^9 –10 ^10 L _⊙ ) galaxies at z = 1–2, and that photometry alone is not effective at identifying AGN within this faint population.

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 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.

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.

Using the first epoch of four-band NIRCam observations obtained by the James Webb Space Telescope (JWST) Prime Extragalactic Areas for Reionization and Lensing Science Program in the Spitzer IRAC Dark Field, we search for F150W and F200W dropouts. In 14.2 arcmin ^2 , we have found eight F150W dropouts and eight F200W dropouts, all brighter than 27.5 mag (the brightest being ∼24 mag) in the band to the red side of the break. As they are detected in multiple bands, these must be real objects. Their nature, however, is unclear, and characterizing their properties is important for realizing the full potential of JWST. If the observed color decrements are due to the Lyman break, these objects should be at z ≳ 11.7 and z ≳ 15.4, respectively. The color diagnostics show that at least four F150W dropouts are far away from the usual contaminators encountered in dropout searches (red galaxies at much lower redshifts or brown dwarf stars). While the diagnostics of the F200W dropouts are less certain due to the limited number of passbands, at least one of them is likely not a known type of contaminant, and the rest are consistent with either high-redshift galaxies with evolved stellar populations or old galaxies at z ≈ 3–8. If a significant fraction of our dropouts are indeed at z ≳ 12, we have to face the severe problem of explaining their high luminosities and number densities. Spectroscopic identifications of such objects are urgently needed.

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 derive the spatial and wavelength behavior of dust attenuation in the multiple-armed spiral galaxy VV 191b using backlighting by the superimposed elliptical system VV 191a in a pair with an exceptionally favorable geometry for this measurement. Imaging using the James Webb Space Telescope and Hubble Space Telescope spans the wavelength range 0.3–4.5 μ m with high angular resolution, tracing the dust in detail from 0.6–1.5 μ m. Distinct dust lanes continue well beyond the bright spiral arms, and trace a complex web, with a very sharp radial cutoff near 1.7 Petrosian radii. We present attenuation profiles and coverage statistics in each band at radii 14–21 kpc. We derive the attenuation law with wavelength; the data both within and between the dust lanes clearly favor a stronger reddening behavior ( R = A _V / E _B _− _V ≈ 2.0 between 0.6 and 0.9 μ m, approaching unity by 1.5 μ m) than found for starbursts and star-forming regions of galaxies. Power-law extinction behavior ∝ λ ^− ^β gives β = 2.1 from 0.6–0.9 μ m. R decreases at increasing wavelengths ( R ≈ 1.1 between 0.9 and 1.5 μ m), while β steepens to 2.5. Mixing regions of different column density flattens the wavelength behavior, so these results suggest a different grain population than in our vicinity. The NIRCam images reveal a lens arc and counterimage from a background galaxy at z ≈ 1, spanning 90° azimuthally at 2.″8 from the foreground elliptical-galaxy nucleus, and an additional weakly lensed galaxy. The lens model and imaging data give a mass/light ratio M / L _B = 7.6 in solar units within the Einstein radius 2.0 kpc.

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 present a catalog of about 25,000 images of massive ( M _⋆ ≥ 10 ^9 M _⊙ ) galaxies at redshifts 3 ≤ z ≤ 6 from the TNG50 cosmological simulation, tailored for observations at multiple wavelengths carried out with JWST. The synthetic images were created with the SKIRT radiative transfer code, including the effects of dust attenuation and scattering. The noiseless images were processed with the mirage simulator to mimic the Near Infrared Camera (NIRCam) observational strategy (e.g., noise, dithering pattern, etc.) of the Cosmic Evolution Early Release Science (CEERS) survey. In this paper, we analyse the predictions of the TNG50 simulation for the size evolution of galaxies at 3 ≤ z ≤ 6 and the expectations for CEERS to probe that evolution. In particular, we investigate how sizes depend on the wavelength, redshift, mass, and angular resolution of the images. We find that the effective radius accurately describes the three-dimensional half-mass–radius of the TNG50 galaxies. Sizes observed at 2 μ m are consistent with those measured at 3.56 μ m at all redshifts and masses. At all masses, the population of higher- z galaxies is more compact than their lower- z counterparts. However, the intrinsic sizes are smaller than the mock observed sizes for the most massive galaxies, especially at z ≲ 4. This discrepancy between the mass and light distributions may point to a transition in the galaxy morphology at z = 4–5, where massive compact systems start to develop more extended stellar structures. ^22

The massive galaxy cluster El Gordo ( z = 0.87) imprints multitudes of gravitationally lensed arcs onto James Webb Space Telescope Near-Infrared Camera (NIRCam) images. Eight bands of NIRCam imaging were obtained in the “Prime Extragalactic Areas for Reionization and Lensing Science” (“PEARLS”) program. Point-spread function–matched photometry across Hubble Space Telescope and NIRCam filters supplies new photometric redshifts. A new light-traces-mass lens model based on 56 image multiplicities identifies the two mass peaks and yields a mass estimate within 500 kpc of (7.0 ± 0.30) × 10 ^14 M _⊙ . A search for substructure in the 140 cluster members with spectroscopic redshifts confirms the two main mass components. The southeastern mass peak that contains the brightest cluster galaxy is more tightly bound than the northwestern one. The virial mass within 1.7 Mpc is (5.1 ± 0.60)×10 ^14 M _⊙ , lower than the lensing mass. A significant transverse velocity component could mean the virial mass is underestimated. We contribute one new member to the previously known z = 4.32 galaxy group. Intrinsic (delensed) positions of the five secure group members span a physical extent of ∼60 kpc. 13 additional candidates selected by spectroscopic/photometric constraints are small and faint, with a mean intrinsic luminosity ∼2.2 mag fainter than L ^* . NIRCam imaging admits a fairly wide range of brightnesses and morphologies for the group members, suggesting a more diverse galaxy population in this galaxy overdensity.

We present a sample of four emission-line galaxies at z = 6.11–6.35 that were serendipitously discovered using the commissioning data for the James Webb Space Telescope (JWST)/NIRCam wide-field slitless spectroscopy mode. One of them (at z = 6.11) has been reported previously, while the others are new discoveries. These sources are selected by the secure detections of both [O iii ] λ 5007 and H α lines with other fainter lines, which were tentatively detected in some cases (e.g., [O ii ] λ 3727, [O iii ] λ 4959). In the [O iii ]/H β –[N ii ]/H α Baldwin–Phillips–Terlevich diagram, these galaxies occupy the same parameter space as that of z ∼ 2 star-forming galaxies, indicating that they have been enriched rapidly to subsolar metallicities (∼0.4 Z _⊙ ), similar to galaxies with comparable stellar masses at much lower redshifts. The detection of strong H α lines suggests a higher ionizing photon production efficiency within galaxies in the early universe. We find brightening of the [O iii ] λ 5007 line-luminosity function (LF) from z = 3 to 6, and weak or no redshift evolution of the H α line LF from z = 2 to 6. Both LFs are underpredicted at z ∼ 6 by a factor of ∼10 in certain cosmological simulations. This further indicates a global Ly α photon escape fraction of 7%–10% at z ∼ 6, which is slightly lower than previous estimates through the comparison of the UV-derived star formation rate density and Ly α luminosity density. Our sample recovers ${66}_{-44}^{+128}$ % of z = 6.0–6.6 galaxies in the survey volume with stellar masses greater than 5 × 10 ^8 M _⊙ , suggesting the ubiquity of strong H α and [O iii ] line emitters in the Epoch of Reionization, which will be further uncovered in the era of JWST.

A full understanding of how unusually large ultradiffuse galaxies (UDGs) fit into our conventional theory of galaxy formation remains elusive, despite the large number of objects identified locally. A natural extension of UDG research is the study of similar galaxies at higher redshift to establish how their properties may evolve over time. However, this has been a challenging task given how severely systematic effects and cosmological surface brightness dimming inhibit our ability to analyze low surface brightness galaxies at high z . Here, we present a sample of low stellar surface density galaxies (LDGs) at moderate redshift, likely the progenitors of local UDGs, identified using deep near-IR observations of the El Gordo cluster at z = 0.87 with JWST. By stacking eight NIRCAM filters, we reach an apparent surface brightness sensitivity of 24.59 mag arcsec ^−2 , which is faint enough to be complete to the bright end of the LDG population. Our analysis identifies significant differences between this population and UDGs observed locally, such as their color and size distributions, which suggest that the UDG progenitors at high z are bluer and more extended than UDGs at z = 0. This suggests that multiple mechanisms are responsible for the UDG formation and that prolonged transformation of cluster dwarfs is not a primary UDG formation mechanism at high z . Furthermore, we find a slight overabundance of LDGs in El Gordo, and, in contrast to findings in local clusters, our analysis does not show a deficit of LDGs in the center of El Gordo, implying that tidal destruction of LDGs is significant between z = 0.87 and z = 0.

Gradients in the mass-to-light ratio of distant galaxies impede our ability to characterize their size and compactness. The long-wavelength filters of JWST’s NIRCam offer a significant step forward. For galaxies at Cosmic Noon ( z ∼ 2), this regime corresponds to the rest-frame near-infrared, which is less biased toward young stars and captures emission from the bulk of a galaxy’s stellar population. We present an initial analysis of an extraordinary lensed dusty star-forming galaxy at z = 2.3 behind the El Gordo cluster ( z = 0.87), named El Anzuelo (“The Fishhook”) after its partial Einstein-ring morphology. The far-UV to near-IR spectral energy distribution suggests an intrinsic star formation rate of ${81}_{-2}^{+7}\,{M}_{\odot }\,{\mathrm{yr}}^{-1}$ and dust attenuation A _V ≈ 1.6, in line with other DSFGs on the star-forming main sequence. We develop a parametric lens model to reconstruct the source-plane structure of dust imaged by the Atacama Large Millimeter/submillimeter Array, far-UV to optical light from Hubble, and near-IR imaging with 8 filters of JWST/NIRCam, as part of the Prime Extragalactic Areas for Reionization and Lensing Science program. The source-plane half-light radius is remarkably consistent from ∼1 to 4.5 μ m, despite a clear color gradient where the inferred galaxy center is redder than the outskirts. We interpret this to be the result of both a radially decreasing gradient in attenuation and substantial spatial offsets between UV- and IR-emitting components. A spatial decomposition of the SED reveals modestly suppressed star formation in the inner kiloparsec, which suggests that we are witnessing the early stages of inside-out quenching.

We present a multiwavelength analysis using the Submillimeter Array (SMA), James Clerk Maxwell Telescope, NOEMA, JWST, the Hubble Space Telescope (HST), and the Spitzer Space Telescope of two dusty strongly star-forming galaxies, 850.1 and 850.2, seen through the massive cluster lens A 1489. These SMA-located sources both lie at z = 4.26 and have bright dust continuum emission, but 850.2 is a UV-detected Lyman-break galaxy, while 850.1 is undetected at ≲ 2 μ m, even with deep JWST/NIRCam observations. We investigate their stellar, interstellar medium, and dynamical properties, including a pixel-level spectral energy distribution analysis to derive subkiloparsec-resolution stellar-mass and A _V maps. We find that 850.1 is one of the most massive and highly obscured, A _V ∼ 5, galaxies known at z > 4 with M _* ∼10 ^11.8 M _⊙ (likely forming at z > 6), and 850.2 is one of the least massive and least obscured, A _V ∼ 1, members of the z > 4 dusty star-forming population. The diversity of these two dust-mass-selected galaxies illustrates the incompleteness of galaxy surveys at z ≳ 3–4 based on imaging at ≲ 2 μ m, the longest wavelengths feasible from HST or the ground. The resolved mass map of 850.1 shows a compact stellar-mass distribution, ${R}_{{\rm{e}}}^{\mathrm{mass}}$ ∼1 kpc, but its expected evolution means that it matches both the properties of massive, quiescent galaxies at z ∼ 1.5 and ultramassive early-type galaxies at z ∼ 0. We suggest that 850.1 is the central galaxy of a group in which 850.2 is a satellite that will likely merge in the near future. The stellar morphology of 850.1 shows arms and a linear bar feature that we link to the active dynamical environment it resides within.

We identify 71 distant stars in James Webb Space Telescope/NIRCam early release observations (ERO) images of the field of galaxy cluster SMACS J0723.3-7327 (SMACS 0723). Given the relatively small (∼10°) angular separation between SMACS 0723 and the Large Magellanic Cloud (LMC), it is likely that these stars are associated with the LMC outskirts or the Leading Arm. This is further bolstered by a spectral energy distribution (SED) analysis, which suggests an excess of stars at a physical distance of 40–100 kpc, consistent with being associated with or located behind the Magellanic system. In particular, we find that the overall surface density of stars brighter than 27.0 mag in the field of SMACS 0723 is ∼2.3 times that of stars in a blank field with similar Galactic latitude (the North Ecliptic Pole Time Domain Field), and that the density of stars in the SMACS 0723 field with SED-derived distances consistent with the Magellanic system is ∼6.1 times larger than that of the blank field. The candidate stars at these distances are consistent with a stellar population at the same distance modulus with [Fe/H] = −1.0 and an age of ∼5.0 Gyr. On the assumption that all of the 71 stars are associated with the LMC, then the stellar density of the LMC at the location of the SMACS 0723 field is ∼740 stars kpc ^−3 , which helps trace the density of stars in the LMC outskirts.