Search

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.

We present a measurement of the low-mass quiescent size–mass relation at cosmic noon (1 < z < 3) from the JWST PRIMER and UNCOVER treasury surveys, which highlights two distinct classes of quiescent galaxies. While the massive population is well studied at these redshifts, the low-mass end has been previously underexplored due to a lack of observing facilities with sufficient sensitivity and spatial resolution. We select a conservative sample of low-mass quiescent galaxy candidates using rest-frame UVJ colors and specific star formation rate criteria and measure galaxy morphology in both rest-frame UV/optical wavelengths (F150W) and rest-frame near-infrared (F444W). We confirm an unambiguous flattening of the low-mass quiescent size–mass relation, which results from the separation of the quiescent galaxy sample into two distinct populations at $\mathrm{log}({M}_{\star }/{M}_{\odot })\sim 10.3$ : low-mass quiescent galaxies that are notably younger and have disky structures, and massive galaxies consistent with spheroidal morphologies and older median stellar ages. These separate populations imply mass quenching dominates at the massive end while other mechanisms, such as environmental or feedback-driven quenching, form the low-mass end. This stellar mass-dependent slope of the quiescent size–mass relation could also indicate a shift from size growth due to star formation (low masses) to growth via mergers (massive galaxies). The transition mass between these two populations also corresponds with other dramatic changes and characteristic masses in several galaxy evolution scaling relations (e.g., star formation efficiency, dust obscuration, and stellar-to-halo mass ratios), further highlighting the stark dichotomy between low-mass and massive galaxy formation.

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

JWST is uncovering the properties of ever-increasing numbers of galaxies at z > 6, during the epoch of reionization. Connecting these observed populations to the process of reionization requires understanding how efficiently they produce Lyman continuum (LyC) photons and what fraction ( f _esc ) of these photons escape into the intergalactic medium. By applying the Cox proportional hazards model, a survival analysis technique, to the Low-redshift Lyman Continuum Survey (LzLCS), we develop new, empirical, multivariate predictions for f _esc . The models developed from the LzLCS reproduce the observed f _esc for z ∼ 3 samples, which suggests that LyC emitters may share similar properties at low and high redshift. Our best-performing models for the z ∼ 3 galaxies include information about dust attenuation, ionization, and/or morphology. We then apply these models to z ≳ 6 galaxies. For large photometric samples, we find a median predicted f _esc = 0.047–0.14. For smaller spectroscopic samples, which may include stronger emission-line galaxies, we find that ≥33% of the galaxies have f _esc > 0.2, and we identify several candidate extreme leakers with f _esc ≥ 0.5. The current samples show no strong trend between predicted f _esc and UV magnitude, but limited spectroscopic information makes this result uncertain. Multivariate predictions can give significantly different results from single-variable predictions, and the predicted f _esc for high-redshift galaxies can differ significantly depending on whether star formation rate surface density or radius is used as a measure of galaxy morphology. We provide all parameters necessary to predict f _esc for additional samples of high-redshift galaxies using these models.

To understand how galaxies reionized the Universe, we must determine how the escape fraction of Lyman continuum (LyC) photons ( f _esc ) depends on galaxy properties. Using the z ∼ 0.3 Low-redshift Lyman Continuum Survey (LzLCS), we develop and analyze new multivariate predictors of f _esc . These predictions use the Cox proportional hazards model, a survival analysis technique that incorporates both detections and upper limits. Our best model predicts the LzLCS f _esc detections with an rms scatter of 0.31 dex, better than single-variable correlations. According to ranking techniques, the most important predictors of f _esc are the equivalent width (EW) of Lyman-series absorption lines and the UV dust attenuation, which track line-of-sight absorption due to H i and dust. The H i absorption EW is uniquely crucial for predicting f _esc for the strongest LyC emitters, which show properties similar to weaker LyC emitters and whose high f _esc may therefore result from favorable orientation. In the absence of H i information, star formation rate surface density (Σ _SFR ) and [O iii ]/[O ii ] ratio are the most predictive variables and highlight the connection between feedback and f _esc . We generate a model suitable for z > 6, which uses only the UV slope, Σ _SFR , and [O iii ]/[O ii ]. We find that Σ _SFR is more important in predicting f _esc at higher stellar masses, whereas [O iii ]/[O ii ] plays a greater role at lower masses. We also analyze predictions for other parameters, such as the ionizing-to-nonionizing flux ratio and Ly α escape fraction. These multivariate models represent a promising tool for predicting f _esc at high redshift.

UVCANDELS is a Hubble Space Telescope Cycle-26 Treasury Program awarded 164 orbits of primary ultraviolet (UV) F275W imaging and coordinated parallel optical F435W imaging in four CANDELS fields—GOODS-N, GOODS-S, EGS, and COSMOS—covering a total area of ∼426 arcmin ^2 . This is ∼2.7 times larger than the area covered by previous deep-field space UV data combined, reaching a depth of about 27 and 28 ABmag (5 σ in 0.”2 apertures) for F275W and F435W, respectively. Along with new photometric catalogs, we present an analysis of the rest-frame UV luminosity function (LF), relying on our UV-optimized aperture photometry method, yielding a factor of 1.5 increase over H-isophot aperture photometry in the signal-to-noise ratios of galaxies in our F275W imaging. Using well-tested photometric redshift measurements, we identify 5810 galaxies at redshifts 0.6 < z < 1, down to an absolute magnitude of M _UV = −14.2. In order to minimize the effect of uncertainties in estimating the completeness function, especially at the faint end, we restrict our analysis to sources above 30% completeness, which provides a final sample of 4726 galaxies at −21.5 < M _UV < −15.5. We performed a maximum likelihood estimate to derive the best-fit parameters of the UV LF. We report a best-fit faint-end slope of $\alpha =-{1.359}_{-0.041}^{+0.041}$ at z ∼ 0.8. Creating subsamples at z ∼ 0.7 and z ∼ 0.9, we observe a possible evolution of α with redshift. The unobscured UV luminosity density at M _UV < −10 is derived as ${\rho }_{{\rm{U}}{\rm{V}}}={1.339}_{-0.030}^{+0.027}\,(\times {10}^{26}\,{\rm{e}}{\rm{r}}{\rm{g}}\,{{\rm{s}}}^{-1}\,{{\rm{H}}{\rm{z}}}^{-1}\,{{\rm{M}}{\rm{p}}{\rm{c}}}^{-3})$ using our best-fit LF parameters. The new F275W and F435 photometric catalogs from UVCANDELS have been made publicly available on the Barbara A. Mikulski Archive for Space Telescopes.

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.

We explore the relationships between the [O/H] gas-phase metallicity radial gradients and multiple galaxy properties for 238 star-forming galaxies at 0.6 < z < 2.6 selected from the CANDELS Ly α Emission at Reionization survey with stellar mass $8.5\lt \mathrm{log}\,{M}_{* }/{M}_{\odot }\lt 10.5$ . The gradients cover the range from −0.11 to 0.22 dex kpc ^−1 , with the median value close to 0. We reconstruct the nonparametric star formation histories (SFHs) of the galaxies with spectral energy distribution modeling using Prospector with more than 40 photometric bands from the Hubble Space Telescope, Spitzer, and ground-based facilities. In general, we find weak or no correlations between the metallicity gradients and most galaxy properties, including the mass-weighted age, recent star formation rate, dust attenuation, and morphology as quantified by both parametric and nonparametric diagnostics. We find a significant but moderate correlation between the gradients and the “evolutionary time,” a temporal metric that characterizes the evolutionary status of a galaxy, with flatter gradients observed in more evolved galaxies. Also, there is evidence that galaxies with multiple star formation episodes in their SFHs tend to develop more negative gas-phase metallicity gradients (higher [O/H] at the center). We conclude that gas kinematics, e.g., radial inflows and outflows, is likely an important process in setting the gas-phase metallicity gradients, in addition to the evolution of the SFH radial profile. Since the gradients are largely independent of the galaxies’ physical properties and only weakly dependent on their SFH, it would appear that the timescale of the gas kinematics is significantly shorter than the evolution of star formation.

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.

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.

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

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

We present an identification of dust-attenuated star-forming galactic-disk substructures in a typical star-forming galaxy (SFG), UDF2, at z = 2.696. To date, substructures containing significant buildup of stellar mass and actively forming stars have yet to be found in typical (i.e., main-sequence) SFGs at z > 2. This is due to the strong dust attenuation common in massive galaxies at the epoch and the scarcity of high-resolution, high-sensitivity extinction-independent imaging. To search for disk substructures, we subtracted the central stellar-mass disk from the JWST/NIRCam rest-frame 1.2 μ m image (0.″13 resolution) and subtracted, in the visibility plane, the central starburst disk from Atacama Large Millimeter/submillimeter Array (ALMA) rest-frame 240 μ m observations (0.″03 resolution). The residual images revealed substructures at rest-frame 1.2 μ m colocated with those found at rest-frame 240 μ m, ≃2 kpc away from the galactic center. The largest substructure contains ≃20% of the total stellar mass and ≃5% of the total star formation rate of the galaxy. While UDF2 exhibits a kinematically ordered velocity field of molecular gas consistent with a secularly evolving disk, more sensitive observations are required to characterize the nature and the origin of this substructure (spiral arms, minor merger, or other types of disk instabilities). UDF2 resides in an overdense region ( N ≥ 4 massive galaxies within 70 kpc projected distance at z = 2.690–2.697) and the substructures may be associated with interaction-induced instabilities. Importantly, a statistical sample of such substructures identified with JWST and ALMA could play a key role in bridging the gap between the bulge-forming starburst and the rest of the galaxy.

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

Quiescent galaxies having more compact morphologies than star-forming galaxies has been a consistent result in the field of galaxy evolution. What is not clear is at what point this divergence happens, i.e., when do quiescent galaxies become compact, and how big of a role does the progenitor effect play in this result? Here we aim to model the morphological and star formation histories of high-redshift (0.8 < z < 1.65) massive galaxies ( $\mathrm{log}(M/{M}_{\odot })$ > 10.2) with stellar population fits using Hubble Space Telescope/WFC3 G102 and G141 grism spectra plus photometry from the CANDELS Ly α Emission at Reionization (CLEAR) survey, constraining the star formation histories for a sample of ∼400 massive galaxies using flexible star formation histories. We develop a novel approach to classifying galaxies by their star formation activity in a way that highlights the green valley population, by modeling the specific star formation rate distributions as a function of redshift and deriving the probability that a galaxy is quiescent ( P _Q ). Using P _Q and our flexible star formation histories we outline the evolutionary paths of our galaxies in relation to stellar mass, Sérsic index, effective radius R _eff , and stellar mass surface density. We find that the galaxies show no appreciable stellar mass growth after entering the green valley (a net decrease of 4%) while their stellar mass surface densities increase by ∼0.25 dex. Therefore galaxies are becoming compact during the green valley and this is due to an increase in the Sérsic index and a decrease in R _eff .

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 analyze a sample of 25 [Ne v ] ( λ 3426) emission-line galaxies at 1.4 < z < 2.3 using Hubble Space Telescope/Wide Field Camera 3 G102 and G141 grism observations from the CANDELS Ly α Emission at Reionization (CLEAR) survey. [Ne v ] emission probes extremely energetic photoionization (creation potential of 97.11 eV) and is often attributed to energetic radiation from active galactic nuclei (AGNs), shocks from supernovae, or an otherwise very hard ionizing spectrum from the stellar continuum. In this work, we use [Ne v ] in conjunction with other rest-frame UV/optical emission lines ([O ii ] λ λ 3726, 3729, [Ne iii ] λ 3869, H β , [O iii ] λ λ 4959, 5007, H α +[N ii ] λ λ 6548, 6583, [S ii ] λ λ 6716, 6731), deep (2–7 Ms) X-ray observations (from Chandra), and mid-infrared imaging (from Spitzer) to study the origin of this emission and to place constraints on the nature of the ionizing engine. The majority of the [Ne v ]-detected galaxies have properties consistent with ionization from AGNs. However, for our [Ne v ]-selected sample, the X-ray luminosities are consistent with local ( z ≲ 0.1) X-ray-selected Seyferts, but the [Ne v ] luminosities are more consistent with those from z ∼ 1 X-ray-selected QSOs. The excess [Ne v ] emission requires either reduced hard X-rays or a ∼0.1 keV excess. We discuss possible origins of the apparent [Ne v ] excess, which could be related to the “soft (X-ray) excess” observed in some QSOs and Seyferts and/or be a consequence of a complex/anisotropic geometry for the narrow-line region, combined with absorption from a warm, relativistic wind ejected from the accretion disk. We also consider implications for future studies of extreme high-ionization systems in the epoch of reionization ( z ≳ 6) with the James Webb Space Telescope.

We present an overview of the CANDELS Ly α Emission At Reionization (CLEAR) survey. CLEAR is a 130 orbit program of the Hubble Space Telescope using the Wide Field Camera 3 (WFC3) IR G102 grism. CLEAR targets 12 pointings divided between the GOODS-N and GOODS-S fields of the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). Combined with existing spectroscopic data from other programs, the full CLEAR data set includes spectroscopic imaging of these fields over 0.8–1.7 μ m. In this paper, we describe the CLEAR survey, the survey strategy, the data acquisition, reduction, processing, and science products and catalogs released alongside this paper. The catalogs include emission line fluxes and redshifts derived from the combination of the photometry and grism spectroscopy for 6048 galaxies, primarily ranging from 0.2 ≲ z ≲ 3. We also provide an overview of CLEAR’s science goals and results. In conjunction with this paper we provide links to electronic versions of the data products, including 1D+2D extracted spectra and emission line maps.

We use the spectral energy distribution fitting code Prospector to reconstruct the nonparametric star formation history (SFH) of massive ( $\mathrm{log}{M}_{* }\gt 10.3$ ) star-forming galaxies (SFGs) and quiescent galaxies (QGs) at redshift z _obs ∼ 2 to investigate the joint evolution of star formation activity and structural properties. Compared to extended SFGs, compact SFGs are more likely to have experienced multiple star formation episodes, with the fractional mass formed during the older (≥1 Gyr) episode being larger, suggesting that high-redshift SFGs assembled their central regions earlier and then kept growing in central mass as they become more compact. The SFH of compact QGs does not significantly differ from the average for this category, and shows an early burst followed by a gradual decline of the star formation rate. The SFH of extended QGs, however, is similar to that of post-starburst galaxies, and their morphology is also frequently disturbed. Knowledge of the SFH also enables us to empirically reconstruct the structural evolution of individual galaxies. While the progenitor effect is clearly observed and accounted for in our analysis, it alone is insufficient to explain the observed structural evolution. We show that, as they evolve from star-forming phase to quiescence, galaxies grow massive dense stellar cores. Quenching begins at the center and then propagates outward to the rest of the structure. We discuss possible physical scenarios for the observed evolution and find that our empirical constraints are in good quantitative agreement with the model predictions from dissipative accretion of gas to the center followed by massive starbursts before final quiescence (wet compaction).

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 a new method based on information theory to find the optimal number of bands required to measure the physical properties of galaxies with desired accuracy. As a proof of concept, using the recently updated COSMOS catalog (COSMOS2020), we identify the most relevant wave bands for measuring the physical properties of galaxies in a Hawaii Two-0- (H20) and UVISTA-like survey for a sample of i < 25 AB mag galaxies. We find that with the available i -band fluxes, r , u , IRAC/ ch 2, and z bands provide most of the information regarding the redshift with importance decreasing from r band to z band. We also find that for the same sample, IRAC/ ch 2, Y , r , and u bands are the most relevant bands in stellar-mass measurements with decreasing order of importance. Investigating the intercorrelation between the bands, we train a model to predict UVISTA observations in near-IR from H20-like observations. We find that magnitudes in the YJH bands can be simulated/predicted with an accuracy of 1 σ mag scatter ≲0.2 for galaxies brighter than 24 AB mag in near-IR bands. One should note that these conclusions depend on the selection criteria of the sample. For any new sample of galaxies with a different selection, these results should be remeasured. Our results suggest that in the presence of a limited number of bands, a machine-learning model trained over the population of observed galaxies with extensive spectral coverage outperforms template fitting. Such a machine-learning model maximally comprises the information acquired over available extensive surveys and breaks degeneracies in the parameter space of template fitting inevitable in the presence of a few bands.

We present a study of spatially resolved star formation histories (SFHs) for 60 z ∼ 2.3 main-sequence, star-forming galaxies selected from the MOSDEF spectroscopic survey in the GOODS-N field, with median stellar mass $\mathrm{log}({M}_{\star }/{M}_{\odot })$ = 9.75 and spanning the range $8.6\lt \mathrm{log}({M}_{\star }/{M}_{\odot })\lt 11.5$ . Photometry is decomposed into a central and an outer spatial component using observed z _F850LP − H _F160W colors. The Prospector code is used to model spectral energy distributions for the center, outskirt, and integrated galaxy using Hubble Space Telescope/ACS and WFC3, Spitzer/IRAC, and ground-based photometry, with additional constraints on gas-phase metallicity and spectroscopic redshift from MOSDEF spectroscopy. For the low-resolution bands, spatially resolved photometry is determined with an iterative approach. The reconstructed SFHs indicate that the majority of galaxies with $\mathrm{log}({M}_{\star }/{M}_{\odot })\lt 10.5$ are observed while their central regions undergo relatively recent (

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 explore how the fraction of quenched galaxies changes in groups of galaxies with respect to the distance to the center of the group, redshift, and stellar mass to determine the dominant process of environmental quenching in 0.2 < z < 0.8 groups. We use new UV data from the UVCANDELS project in addition to existing multiband photometry to derive new galaxy physical properties of the group galaxies from the zCOSMOS 20 k group catalog. Limiting our analysis to a complete sample of log ( M _* / M _⊙ ) > 10.56 group galaxies, we find that the probability of being quenched increases slowly with decreasing redshift, diverging from the stagnant field galaxy population. A corresponding analysis on how the probability of being quenched increases with time within groups suggests that the dominant environmental quenching process is characterized by slow (∼Gyr) timescales. We find a quenching time of approximately ${4.91}_{-1.47}^{+0.91}$ Gyr, consistent with the slow processes of strangulation and delayed-then-rapid quenching although more data are needed to confirm this result.

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.