Search

In this paper, we describe the “Medium Bands, Mega Science” JWST Cycle 2 survey (JWST-GO-4111) and demonstrate the power of these data to reveal both the spatially integrated and spatially resolved properties of galaxies from the local Universe to the era of cosmic dawn. Executed in 2023 November, MegaScience obtained ∼30 arcmin ^2 of deep multiband NIRCam imaging centered on the z ∼ 0.3 A2744 cluster, including 11 medium-band filters and the two shortest-wavelength broadband filters, F070W and F090W. Together, MegaScience and the UNCOVER Cycle 1 treasury program provide a complete set of deep (∼28–30 mag _AB ) images in all NIRCam medium- and broadband filters. This unique data set allows us to precisely constrain photometric redshifts, map stellar populations and dust attenuation for large samples of distant galaxies, and examine the connection between galaxy structures and formation histories. MegaScience also includes ∼17 arcmin ^2 of NIRISS parallel imaging in two broadband and four medium-band filters from 0.9 to 4.8 μ m, expanding the footprint where robust spectral energy distribution (SED) fitting is possible. We provide example SEDs and multiband cutouts at a variety of redshifts, and use a catalog of JWST spectroscopic redshifts to show that MegaScience improves both the scatter and catastrophic outlier rate of photometric redshifts by factors of 2–3. Additionally, we demonstrate the spatially resolved science enabled by MegaScience by presenting maps of the [O iii ] line emission and continuum emission in three spectroscopically confirmed z > 6 galaxies. We show that line emission in reionization-era galaxies can be clumpy, extended, and spatially offset from continuum emission, implying that galaxy assembly histories are complex even at these early epochs. We publicly release fully reduced mosaics and photometric catalogs for both the NIRCam primary and NIRISS parallel fields ( http://jwst-uncover.github.io/megascience ).

We present an updated mass and magnification model of galaxy cluster A370 using new NIRCam, NIRISS, and NIRSpec data from the Canadian NIRISS Unbiased Cluster Survey (CANUCS). Using Lenstool and a combination of archival Hubble Space Telescope (HST) and MUSE data with new JWST data as constraints, we derive an improved gravitational lensing model and extract magnifications of background galaxies with uncertainties. Using our best fit model, we perform a search for new multiply imaged systems via predicted positions. We report no new multiply imaged systems with identifiable redshifts, likely due to already very deep HST and Spitzer data but confirm a z ∼ 8 multiply imaged system by measuring its redshift with NIRISS and NIRSpec spectra. We find that the overall shape of the critical curve for a source at z = 9.0 is similar to previous models of A370, with small changes. We investigate the z ∼ 8 galaxy with two images observable with an apparent magnitude in the F125W band of 26.0 ± 0.2 and 25.6 ± 0.1. After correcting for the magnifications of the images, ${7.4}_{-0.3}^{+0.5}$ and ${9.4}_{-0.4}^{+0.5}$ , we use spectral energy distribution fitting to find an intrinsic stellar mass of log( M ^* / M _⊙ ) = ${7.49}_{-0.05}^{+0.04}$ , intrinsic star formation rate of ${2.8}_{-0.3}^{+0.4}$ M _⊙ yr ^−1 , and M _UV of −21.3 ${}_{-0.2}^{+0.2}$ , which is close to the knee of the luminosity function at that redshift. Our model, and corresponding magnification, shear, and convergence maps are available on request and will be made publicly available on MAST in a CANUCS data release (doi:10.17909/ph4n-6n76).

We report the spectroscopic confirmation of a massive ( $\mathrm{log}({M}_{\star }/{M}_{\odot })=10.34{\pm }_{0.07}^{0.06}$ ), Hubble Space Telescope–dark ( m _F150W − m _F444W = 3.6) quiescent galaxy at z _spec = 3.97 in the UNCOVER survey. NIRSpec/PRISM spectroscopy and a nondetection in deep Atacama Large Millimeter/submillimeter Array imaging surprisingly reveals that the galaxy is consistent with a low (

In this paper we describe the survey design for the Ultradeep NIRSpec and NIRCam Observations before the Epoch of Reionization (UNCOVER) Cycle 1 JWST Treasury program, which executed its early imaging component in 2022 November. The UNCOVER survey includes ultradeep (∼29–30AB) imaging of ∼45 arcmin ^2 on and around the well-studied A2744 galaxy cluster at z = 0.308 and will follow up ∼500 galaxies with extremely deep low-resolution spectroscopy with the NIRSpec/PRISM during the summer of 2023, with repeat visits in summer 2024. We describe the science goals, survey design, target selection, and planned data releases. We also present and characterize the depths of the first NIRCam imaging mosaic, highlighting previously unparalleled resolved and ultradeep 2–4 μ m imaging of known objects in the field. The UNCOVER primary NIRCam mosaic spans 28.8 arcmin ^2 in seven filters (F115W, F150W, F200W, F277W, F356W, F410M, and F444W) and 16.8 arcmin ^2 in our NIRISS parallel (F115W, F150W, F200W, F356W, and F444W). To maximize early community use of the Treasury data set, we publicly release the full reduced mosaics of public JWST imaging including 45 arcmin ^2 NIRCam and 17 arcmin ^2 NIRISS mosaics on and around the A2744 cluster, including the Hubble Frontier Field primary and parallel footprints.

We present an overview and first results from the Spectroscopic Ultradeep Survey Probing Extragalactic Near-infrared Stellar Emission (SUSPENSE), executed with NIRSpec on JWST. The primary goal of the SUSPENSE program is to characterize the stellar, chemical, and kinematic properties of massive quiescent galaxies at cosmic noon. In a single deep NIRSpec/MSA configuration, we target 20 distant quiescent galaxy candidates ( z = 1–3, H _AB ≤ 23), as well as 53 star-forming galaxies at z = 1–4. With 16 hr of integration and the G140M-F100LP dispersion-filter combination, we observe numerous Balmer and metal absorption lines for all quiescent candidates. We derive stellar masses (log M _* / M _⊙ ∼ 10.2–11.5) and detailed star formation histories (SFHs) and show that all 20 candidate quiescent galaxies indeed have quenched stellar populations. These galaxies show a variety of mass-weighted ages (0.8–3.3 Gyr) and star formation timescales (∼0.5–4 Gyr), and four out of 20 galaxies were already quiescent by z = 3. On average, the z > 1.75 [ z < 1.75] galaxies formed 50% of their stellar mass before z = 4 [ z = 3]. Furthermore, the typical SFHs of the galaxies in these two redshift bins ( z _mean = 2.2 [1.3]) indicate that galaxies at higher redshift formed earlier and over shorter star formation timescales compared to lower redshifts. Although this evolution is naturally explained by the growth of the quiescent galaxy population over cosmic time, number density calculations imply that mergers and/or late-time star formation also contribute to the evolution. In future work, we will further unravel the early formation, quenching, and late-time evolution of these galaxies by extending this work with studies on their chemical abundances, resolved stellar populations, and kinematics.

We present JWST/NIRSpec prism spectroscopy of MACS0647−JD, a triply lensed z ∼ 11 candidate discovered in Hubble Space Telescope imaging and spatially resolved by JWST imaging into two components, A and B. Spectroscopy of component A yields a spectroscopic redshift z = 10.17 based on seven detected emission lines: C iii ] λ λ 1907, 1909, [O ii ] λ 3727, [Ne iii ] λ 3869, [Ne iii ] λ 3968, H δ λ 4101, H γ λ 4340, and [O iii ] λ 4363. These are the second-most distant detections of these emission lines to date, in a galaxy observed just 460 million years after the Big Bang. Based on observed and extrapolated line flux ratios we derive a gas-phase metallicity 12 + log(O/H) ∼ 7.5–8.0, or Z ∼ (0.06–0.2) Z _⊙ , ionization parameter $\mathrm{log}(U)$ = −1.9 ± 0.2, and an ionizing photon production efficiency $\mathrm{log}({\xi }_{\mathrm{ion}})=25.2\pm 0.2\,$ erg ^−1 Hz. The spectrum has a softened Ly α break, evidence for a strong Ly α damping wing. The Ly α damping wing also suppresses the F150W photometry, explaining the slightly overestimated photometric redshift z = 10.6 ± 0.3. MACS0647−JD has a stellar mass log( M / M _⊙ ) = 8.1 ± 0.3, including ∼6 × 10 ^7 M _⊙ in component A, most of which formed recently (within ∼20 Myr) with a star formation rate ∼ 2 ± 1 M _⊙ yr ^−1 , all within an effective radius 70 ± 24 pc. Spectroscopy of a fainter companion galaxy C separated by a distance of ∼ 3 kpc reveals a Lyman break consistent with z ∼ 10.17. MACS0647−JD is likely the most distant galaxy merger known.

We present results from the JWST First Reionization Epoch Spectroscopically Complete Observations survey on the star-forming sequence (SFS) of galaxies at 1.0 < z < 1.7, around the peak of the cosmic star formation history. Star formation rates (SFRs) are measured from the redshifted, relatively dust-insensitive Paschen- α emission line, and stellar mass measurements include the F444W (4.4 μ m; rest-frame H ) band. We find SFRs of galaxies with log( M _* / M _⊙ ) > 9.5 that are lower than found in many earlier studies by up to 0.6 dex, but in good agreement with recent results obtained with the Prospector fitting framework. The difference (log(SFR(Pa α )-SFR( Prospector )) is −0.09 ± 0.04 dex at 10 ^10−11 M _⊙ . We also measure the empirical relation between Paschen- α luminosity and rest-frame H- band magnitude and find that the scatter is only 0.04 dex lower than that of the SFR– M * relation and is much lower than the systematic differences among relations in the literature due to various methods of converting observed measurements to physical properties. We additionally identify examples of sources—that, with standard cutoffs via the UVJ diagram, would be deemed quiescent—with significant (log(sSFR)> −11 yr ^−1 ), typically extended, Paschen- α emission. Our results may be indicative of the potential unification of methods used to derive the SFS with careful selection of star-forming galaxies and independent SFR and stellar mass indicators.

We present the galaxy stellar mass–size relation in the rest-frame near-IR (1.5 μ m) and its evolution with redshift up to z = 2.5. Sérsic profiles are measured for ∼26,000 galaxies with stellar masses M _⋆ > 10 ^9 M _⊙ from JWST/NIRCam F277W and F444W imaging provided by the COSMOS-WEB and PRIMER surveys using coordinates, redshifts, colors, and stellar mass estimates from the COSMOS2020 catalog. The new rest-frame near-IR effective radii are generally smaller than previously measured rest-frame optical sizes, on average by 0.14 dex, with no significant dependence on redshift. For quiescent galaxies, this size offset does not depend on stellar mass, but for star-forming galaxies, the offset increases from −0.1 dex at M _⋆ = 10 ^9.5 M _⊙ to −0.25 dex at M _⋆ > 10 ^11 M _⊙ . That is, we find that the near-IR stellar mass–size relation for star-forming galaxies is flatter in the rest-frame near-IR than in the rest-frame optical at all redshifts 0.5 < z < 2.5. The general pace of size evolution is the same in the near-IR as previously demonstrated in the optical, with slower evolution ( R _e ∝ (1 + z ) ^−0.7 ) for L * star-forming galaxies and faster evolution ( R _e ∝ (1 + z ) ^−1.3 ) for L * quiescent galaxies. Massive ( M _⋆ > 10 ^11 M _⊙ ) star-forming galaxies evolve in size almost as fast as quiescent galaxies. Low-mass ( M _⋆ < 10 ^10 M _⊙ ) quiescent galaxies evolve as slow as star-forming galaxies. Our main conclusion is that the size evolution narrative as it has emerged over the past two decades does not radically change when accessing the rest-frame near-IR with JWST, a better proxy of the underlying stellar mass distribution.

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.

We present a study of the growth of the quiescent galaxy population between 0.5 < z < 3 by tracing the number density and structural evolution of a sample of 4518 old and 583 young quiescent galaxies with log( M _⋆ / M _⊙ ) > 10.4, selected from the COSMOS2020 catalog with complementary Hubble Space Telescope F160W imaging from the 3D-DASH survey. Among the quiescent population at z ∼ 2, roughly 50% are recently quenched galaxies; these young quiescent galaxies become increasingly rare toward lower redshift, supporting the idea that the peak epoch of massive galaxy quenching occurred at z > 2. Our data show that while the effective half-light radius of quiescent galaxies generally increases with time, young quiescent galaxies are significantly smaller than their older counterparts at the same redshift. In this work we investigate the connection between this size difference and other structural properties, including axis ratio, color gradients, stellar mass, and the intrinsic scatter in effective radius. We demonstrate that the size difference is driven by the most massive subpopulation (log( M _⋆ / M _⊙ ) > 11) and does not persist when restricting the sample to intermediate-mass galaxies (10.4 < log( M _⋆ / M _⊙ ) < 11). Interestingly, the intrinsic scatter in physical size shows a strong coevolution over the investigated time period and peaks around z ∼ 2 for both populations, only diverging at z < 1. Taken together, and assuming we are not missing a significant population of lower surface brightness galaxies, while the formation and quenching mechanisms that dominate at higher redshifts yield compact remnants, multiple evolutionary pathways may explain the diverse morphologies of galaxies that quench at z < 1.

We present the construction of a deep multiwavelength point-spread-function-matched photometric catalog in the Ultra-Deep Survey (UDS) field following the final UKIDSS UDS release. The catalog includes photometry in 24 filters, from the MegaCam- uS 0.38 μ m band to the Spitzer-IRAC 8 μ m band, over ∼0.9 deg ^2 and with a 5 σ depth of 25.3 AB in the K -band detection image. The catalog, containing ≈188,564 (136,235) galaxies at 0.2 < z < 8.0 with stellar mass $\mathrm{log}({M}_{* }/{M}_{\odot })\gt 8$ and K -band total magnitude K < 25.2 (24.3) AB, enables a range of extragalactic studies. We also provide photometric redshifts, corresponding redshift probability distributions, and rest-frame absolute magnitudes and colors derived using the template-fitting code eazy-py . Photometric redshift errors are less than 3%−4% at z < 4 across the full brightness range in the K band and stellar mass range $8\lt \mathrm{log}({M}_{* }/{M}_{\odot })\lt 12$ . Stellar population properties (e.g., stellar mass, star formation rate, dust extinction) are derived from the modeling of the spectral energy distributions using the codes FAST and Dense Basis.

Theoretical models predict that z ≳ 6 quasars are hosted in the most massive halos of the underlying dark matter distribution and thus would be immersed in protoclusters of galaxies. However, observations report inconclusive results. We investigate the 1.1 proper-Mpc ^2 environment of the z = 7.54 luminous quasar ULAS J1342+0928. We search for Lyman-break galaxy (LBG) candidates using deep imaging from the Hubble Space Telescope (HST) in the Advanced Camera for Surveys (ACS)/F814W, Wide Field Camera 3 (WFC3)/F105W/F125W bands, and Spitzer/Infrared Array Camera at 3.6 and 4.5 μ m. We report a ${z}_{\mathrm{phot}}={7.69}_{-0.23}^{+0.33}$ LBG with mag _F125W = 26.41 at 223 projected proper kpc (pkpc) from the quasar. We find no HST counterpart to one [C ii ] emitter previously found with the Atacama Large millimeter/submillimeter Array (ALMA) at 27 projected pkpc and z _[C II] =7.5341 ± 0.0009 (Venemans et al. 2020 ). We estimate the completeness of our LBG candidates using results from Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey/GOODS deep blank field searches sharing a similar filter setup. We find that >50% of the z ∼ 7.5 Lyman-break galaxies (LBGs) with mag _F125W > 25.5 are missed due to the absence of a filter redward of the Lyman break in F105W, hindering the UV color accuracy of the candidates. We conduct a QSO-LBG clustering analysis revealing a low LBG excess of ${0.46}_{-0.08}^{+1.52}$ in this quasar field, consistent with an average or low-density field. Consequently, this result does not present strong evidence of an LBG overdensity around ULAS J1342+0928. Furthermore, we identify two LBG candidates with a z _phot matching a confirmed z = 6.84 absorber along the line of sight to the quasar. All these galaxy candidates are excellent targets for follow-up observations with JWST and/or ALMA to confirm their redshift and physical properties.

To study the role of environment in galaxy evolution, we reconstruct the underlying density field of galaxies based on COSMOS2020 (The Farmer catalog) and provide the density catalog for a magnitude-limited ( K _s < 24.5) sample of ∼210,000 galaxies at 0.4 < z < 5 within the COSMOS field. The environmental densities are calculated using a weighted kernel density estimation approach with the choice of a von Mises–Fisher kernel, an analog of the Gaussian kernel for periodic data. Additionally, we make corrections for the edge effect and masked regions in the field. We utilize physical properties extracted by LePhare to investigate the connection between star formation activity and the environmental density of galaxies in six mass-complete subsamples at different cosmic epochs within 0.4 < z < 4. Our findings confirm a strong anticorrelation between star formation rate (SFR)/specific SFR (sSFR) and environmental density out to z ∼ 1.1. At 1.1 < z < 2, there is no significant correlation between SFR/sSFR and density. At 2 < z < 4, we observe a reversal of the SFR/sSFR–density relation such that both SFR and sSFR increase by a factor of ∼10 with increasing density contrast, δ , from −0.4 to 5. This observed reversal at higher redshifts supports the scenario where an increased availability of gas supply, along with tidal interactions and a generally higher star formation efficiency in dense environments, could potentially enhance star formation activity in galaxies located in rich environments at z > 2.

Sub/millimeter galaxies are a key population for the study of galaxy evolution because the majority of star formation at high redshifts occurred in galaxies deeply embedded in dust. To search for this population, we have performed an extensive survey with Atacama Large Millimeter/submillimeter Array (ALMA), called the ALMA Lensing Cluster Survey (ALCS). This survey covers 133 arcmin ^2 area and securely detects 180 sources at z ∼ 0.5–6 with a flux limit of ∼0.2 mJy at 1.2 mm. Here, we report the results of multiwavelength spectral energy distribution analysis of the whole ALCS sample, utilizing the observed-frame UV to millimeter photometry. We find that the majority of the ALCS sources lie on the star-forming main sequence, with a smaller fraction showing intense starburst activities. The ALCS sample contains high infrared-excess sources ( $\mathrm{IRX}=\mathrm{log}({L}_{\mathrm{dust}}/{L}_{\mathrm{UV}})\gt 1$ ), including two extremely dust-obscured galaxies (IRX > 5). We also confirm that the ALCS sample probes a broader range in lower dust mass than conventional submillimeter galaxy samples in the same redshift range. We identify six heavily obscured active galactic nucleus (AGN) candidates that are not detected in the archival Chandra data in addition to the three X-ray AGNs reported by Uematsu et al. (2023). The inferred AGN luminosity density shows a possible excess at z = 2–3 compared with that determined from X-ray surveys below 10 keV.

We present the results of a systematic study of the rest-frame optical morphology of quiescent galaxies at z ≥ 3 using the Near-Infrared Camera (NIRCam) on board the James Webb Space Telescope (JWST). Based on a sample selected by UVJ color or NUVUVJ color, we focus on 26 quiescent galaxies with $9.8\lt \mathrm{log}({M}_{\star }/{M}_{\odot })\lt 11.4$ at 2.8 < z _phot < 4.6 with publicly available JWST data. Their sizes are constrained by fitting the Sérsic profile to all available NIRCam images. We see a negative correlation between the observed wavelength and the size and derive their size at the rest frame 0.5 μ m using size measurements in multiple bands. Our quiescent galaxies show a significant correlation between the rest-frame 0.5 μ m size and the stellar mass at z ≥ 3. The analytical fit for them at $\mathrm{log}({M}_{\star }/{M}_{\odot })\gt 10.3$ implies that our size–stellar mass relations are below those at lower redshifts, with the amplitude of ∼0.6 kpc at M _⋆ = 5 × 10 ^10 M _⊙ . This value agrees with the extrapolation of the size evolution of quiescent galaxies at z < 3 in the literature, implying that the size of quiescent galaxies increases monotonically from z ∼ 3–5. Our sample mainly comprises galaxies with bulge-like structures according to their median Sérsic index and axis ratio of n ∼ 3–4 and q ∼ 0.6–0.8, respectively. On the other hand, there is a trend of increasing fraction of galaxies with low Sérsic index at higher redshift, suggesting 3 < z < 5 might be the epoch of onset of morphological transformation with a fraction of very notable disky quenched galaxies.

We present Atacama Large Millimeter/submillimeter Array (ALMA) deep spectroscopy for a lensed galaxy at z _spec = 8.496 with $\mathrm{log}({M}_{\mathrm{star}}/{M}_{\odot })\sim 7.8$ whose optical nebular lines and stellar continuum are detected by JWST/NIRSpec and NIRCam Early Release Observations in the field of SMACS J0723.3–7327. Our ALMA spectrum shows [O iii ] 88 μ m and [C ii ] 158 μ m line detections at 4.0 σ and 4.5 σ , respectively. The redshift and position of the [O iii ] line coincide with those of the JWST source, while the [C ii ] line is blueshifted by 90 km s ^−1 with a spatial offset of 0.″5 (≈0.5 kpc in the source plane) from the centroid of the JWST source. The NIRCam F444W image, including [O iii ] λ 5007 and H β line emission, spatially extends beyond the stellar components by a factor of >8. This indicates that the z = 8.5 galaxy has already experienced strong outflows as traced by extended [O iii ] λ 5007 and offset [C ii ] emission, which would promote ionizing photon escape and facilitate reionization. With careful slit-loss corrections and the removal of emission spatially outside the galaxy, we evaluate the [O iii ] 88 μ m/ λ 5007 line ratio, and derive the electron density n _e by photoionization modeling to be ${220}_{-130}^{+230}$ cm ^−3 , which is comparable with those of z ∼ 2–3 galaxies. We estimate an [O iii ] 88 μ m/[C ii ] 158 μ m line ratio in the galaxy of >4, as high as those of known z ∼ 6–9 galaxies. This high [O iii ] 88 μ m/[C ii ] 158 μ m line ratio is generally explained by the high n _e as well as the low metallicity ( ${Z}_{\mathrm{gas}}/{Z}_{\odot }={0.04}_{-0.02}^{+0.02}$ ), high ionization parameter ( $\mathrm{log}U\gt -2.27$ ), and low carbon-to-oxygen abundance ratio (log(C/O) = [−0.52: −0.24]) obtained from the JWST/NIRSpec data; further [C ii ] follow-up observations will constrain the covering fraction of photodissociation regions.

The James Webb Space Telescope is revealing a new population of dust-reddened broad-line active galactic nuclei (AGN) at redshifts z ≳ 5. Here we present deep NIRSpec/Prism spectroscopy from the Cycle 1 Treasury program Ultradeep NIRSpec and NIRCam ObserVations before the Epoch of Reionization (UNCOVER) of 15 AGN candidates selected to be compact, with red continua in the rest-frame optical but with blue slopes in the UV. From NIRCam photometry alone, they could have been dominated by dusty star formation or an AGN. Here we show that the majority of the compact red sources in UNCOVER are dust-reddened AGN: 60% show definitive evidence for broad-line H α with a FWHM > 2000 km s ^−1 , 20% of the current data are inconclusive, and 20% are brown dwarf stars. We propose an updated photometric criterion to select red z > 5 AGN that excludes brown dwarfs and is expected to yield >80% AGN. Remarkably, among all z _phot > 5 galaxies with F277W – F444W > 1 in UNCOVER at least 33% are AGN regardless of compactness, climbing to at least 80% AGN for sources with F277W – F444W > 1.6. The confirmed AGN have black hole masses of 10 ^7 –10 ^9 M _⊙ . While their UV luminosities (−16 > M _UV > −20 AB mag) are low compared to UV-selected AGN at these epochs, consistent with percent-level scattered AGN light or low levels of unobscured star formation, the inferred bolometric luminosities are typical of 10 ^7 –10 ^9 M _⊙ black holes radiating at ∼10%–40% the Eddington limit. The number densities are surprisingly high at ∼10 ^−5 Mpc ^−3 mag ^−1 , 100 times more common than the faintest UV-selected quasars, while accounting for ∼1% of the UV-selected galaxies. While their UV faintness suggests they may not contribute strongly to reionization, their ubiquity poses challenges to models of black hole growth.

Characterizing the prevalence and properties of faint active galactic nuclei (AGNs) in the early Universe is key for understanding the formation of supermassive black holes (SMBHs) and determining their role in cosmic reionization. We perform a spectroscopic search for broad H α emitters at z ≈ 4–6 using deep JWST/NIRCam imaging and wide field slitless spectroscopy from the EIGER and FRESCO surveys. We identify 20 H α lines at z = 4.2–5.5 that have broad components with line widths from ∼1200–3700 km s ^−1 , contributing ∼30%–90% of the total line flux. We interpret these broad components as being powered by accretion onto SMBHs with implied masses ∼10 ^7–8 M _⊙ . In the UV luminosity range M _UV,AGN+host = −21 to −18, we measure number densities of ≈10 ^−5 cMpc ^−3 . This is an order of magnitude higher than expected from extrapolating quasar UV luminosity functions (LFs). Yet, such AGN are found in only

The era of the James Webb Space Telescope ushers stellar population models into uncharted territories, particularly at the high-redshift frontier. In a companion paper, we apply the Prospector Bayesian framework to jointly infer galaxy redshifts and stellar population properties from broadband photometry as part of the UNCOVER survey. Here we present a comprehensive error budget in spectral energy distribution (SED) modeling. Using a sample selected to have photometric redshifts higher than 9, we quantify the systematic shifts stemming from various model choices in inferred stellar mass, star formation rate (SFR), and age. These choices encompass different timescales for changes in the star formation history (SFH), nonuniversal stellar initial mass functions (IMF), and the inclusion of variable nebular abundances, gas density, and ionizing photon budget. We find that the IMF exerts the strongest influence on the inferred properties: the systematic uncertainties can be as much as 1 dex, 2–5 times larger than the formal reported uncertainties in mass and SFR, and importantly, exceed the scatter seen when using different SED fitting codes. Although the assumptions on the lower end of the IMF induce degeneracy, our findings suggest that a common practice in the literature of assessing uncertainties in SED-fitting processes by comparing multiple codes is substantively underestimating the true systematic uncertainty. Highly stochastic SFHs change the inferred SFH by much larger than the formal uncertainties, and introduce ∼0.8 dex systematics in SFR averaged over a short timescale and ∼0.3 dex systematics in average age. Finally, employing a flexible nebular emission model causes ∼0.2 dex systematic increase in mass and SFR, comparable to the formal uncertainty. This paper constitutes an initial step toward a complete uncertainty estimate in SED modeling.

We report JWST/NIRSpec spectra of three distant T-type brown dwarfs identified in the Ultradeep NIRSpec and NIRCam ObserVations before the Epoch of Reionization (UNCOVER) survey of the Abell 2744 lensing field. One source was previously reported as a candidate T dwarf on the basis of NIRCam photometry, while two sources were initially identified as candidate active galactic nuclei. Low-resolution 1–5 μ m spectra confirm the presence of molecular features consistent with T dwarf atmospheres, and comparison to spectral standards infers classifications of sdT1, T6, and T8–T9. The warmest source, UNCOVER-BD-1, shows evidence of subsolar metallicity, and atmosphere model fits indicate T _eff = 1300 K and [M/H] ∼ −1.0, making this one of the few spectroscopically confirmed T subdwarfs known. The coldest source, UNCOVER-BD-3, is near the T/Y dwarf boundary with T _eff = 550 K, and our analysis indicates the presence of PH _3 in the 3–5 μ m region, favored over CO _2 and a possible indicator of subsolar metallicity. We estimate distances of 0.9–4.5 kpc from the Galactic midplane, making these the most distant brown dwarfs with spectroscopic confirmation. Population simulations indicate high probabilities of membership in the Galactic thick disk for two of these brown dwarfs, and potential halo membership for UNCOVER-BD-1. Our simulations indicate that there are approximately 5 T dwarfs and 1–2 L dwarfs in the Abell 2744 field down to F444W = 30 AB mag, roughly one-third of which are thick disk members. These results highlight the utility of deep JWST/NIRSpec spectroscopy for identifying and characterizing the oldest metal-poor brown dwarfs in the Milky Way.

We demonstrate that the inference of galaxy stellar masses via spectral energy distribution (SED) fitting techniques for galaxies formed in the first billion years after the Big Bang carries fundamental uncertainties owing to the loss of star formation history (SFH) information from the very first episodes of star formation in the integrated spectra of galaxies. While this early star formation can contribute substantially to the total stellar mass of high-redshift systems, ongoing star formation at the time of detection outshines the residual light from earlier bursts, hampering the determination of accurate stellar masses. As a result, order-of-magnitude uncertainties in stellar masses can be expected. We demonstrate this potential problem via direct numerical simulation of galaxy formation in a cosmological context. In detail, we carry out two cosmological simulations with significantly different stellar feedback models, which span a significant range in SFH burstiness. We compute the mock SEDs for these model galaxies at z = 7 via calculations of 3D dust radiative transfer, and then backward fit these SEDs with prospector SED fitting software. The uncertainties in derived stellar masses that we find for z > 7 galaxies motivate the development of new techniques and/or priors for SFH to model star formation in the early Universe.

The selection of high-redshift galaxies often involves spectral energy distribution (SED) fitting to photometric data, an expectation for contamination levels, and measurement of sample completeness—all vetted through comparison to spectroscopic redshift measurements of a sub-sample. The first JWST data are now being taken over several extragalactic fields to different depths and across various areas, which will be ideal for the discovery and classification of galaxies out to distances previously uncharted. As spectroscopic redshift measurements for sources in this epoch will not be initially available to compare with the first photometric measurements of z > 8 galaxies, robust photometric redshifts are of the utmost importance. Galaxies at z > 8 are expected to have bluer rest-frame ultraviolet (UV) colors than typically used model SED templates, which could lead to catastrophic photometric redshift failures. We use a combination of BPASS and Cloudy models to create a supporting set of templates that match the predicted rest-UV colors of z > 8 simulated galaxies. We test these new templates by fitting simulated galaxies in a mock catalog, Yung et al., which mimic expected field depths and areas of the JWST Cosmic Evolution Early Release Science Survey ( m _5 _σ ∼ 28.6 over ∼100 arcmin ^2 ). We use EAZY to highlight the improvements in redshift recovery with the inclusion of our new template set and suggest criteria for selecting galaxies at 8 < z < 10 with the JWST, providing an important test case for observers venturing into this new era of astronomy.

We examine the wavelength dependence of radial light profiles based on Sérsic index n measurements of 1067 galaxies with M _* ≥ 10 ^9.5 M _⊙ and in the redshift range 0.5 < z < 3. The sample and rest-frame optical light profiles are drawn from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) and 3D Hubble Space Telescope (HST); rest-frame near-infrared light profiles are inferred from images collected for the Cosmic Evolution Early Release Science (CEERS) program with the Near Infrared Camera (NIRCam) on board of the James Webb Space Telescope (JWST). n shows only a weak dependence on the wavelength, regardless of the redshift, galaxy mass, and type. On average, star-forming galaxies have n = 1–1.5 and quiescent galaxies have n = 3–4 in the rest-frame optical and near-infrared. The strong correlation at all wavelengths between n and star formation activity implies a physical connection between the radial stellar mass profile and star formation activity. The main caveat is that the current sample is too small to discern trends for the most massive galaxies ( M _* > 10 ^11 M _⊙ ).

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

We report the discovery of four galaxy candidates observed 450–600 Myr after the Big Bang with photometric redshifts between z ∼ 8.3 and 10.2 measured using James Webb Space Telescope (JWST) NIRCam imaging of the galaxy cluster WHL0137−08 observed in eight filters spanning 0.8–5.0 μ m, plus nine Hubble Space Telescope filters spanning 0.4–1.7 μ m. One candidate is gravitationally lensed with a magnification of μ ∼ 8, while the other three are located in a nearby NIRCam module with expected magnifications of μ ≲ 1.1. Using SED fitting, we estimate the stellar masses of these galaxies are typically in the range $\mathrm{log}{M}_{\star }/{M}_{\odot }$ = 8.3–8.7. All appear young, with mass-weighted ages

We examine the factors responsible for the variation in the ionization parameter ( U ) of high-redshift star-forming galaxies based on medium-resolution JWST/NIRSpec observations obtained by the Cosmic Evolution Early Release Science survey. The sample consists of 48 galaxies with redshifts z _spec = 2.7−6.3, which are largely representative of typical galaxies at these redshifts. The [S ii ] λ λ 6718, 6733 doublet is used to estimate electron densities ( n _e ), and dust-corrected H α luminosities are used to compute ionizing photon rates ( Q ). Using composite spectra of galaxies in bins of [O iii ] λ λ 4960, 5008/[O ii ] λ λ 3727, 3730 (O32) as a proxy for U , we determine that galaxies with higher O32 have 〈 n _e 〉 ≃ 500 cm ^−3 that are ≳5 × larger than that of lower-O32 galaxies. We do not find a significant difference in 〈 Q 〉 between low- and high-O32 galaxies. Photoionization modeling indicates a large spread in $\mathrm{log}U$ of ≈1.5 dex at a fixed Z _neb . On the other hand, the data indicate a highly significant correlation between U and star-formation-rate surface density (Σ _SFR ), which appears to be redshift invariant at z ∼ 1.6−6.3, and possibly up to z ∼ 9.5. We consider several avenues through which metallicity and Σ _SFR (or gas density) may influence U , including variations in n _e and Q , internal dust extinction of ionizing photons, and the effects of gas density on the volume filling fraction. Based on these considerations, we conclude that gas density may play a more central role than metallicity in modulating U at these redshifts.

We present images and a multiwavelength photometric catalog based on all of the JWST NIRCam observations obtained to date in the region of the Abell 2744 galaxy cluster. These data come from three different programs, namely, the GLASS-JWST Early Release Science Program, UNCOVER, and the Director’s Discretionary Time program 2756. The observed area in the NIRCam wide-band filters—covering the central and extended regions of the cluster, as well as new parallel fields—is 46.5 arcmin ^2 in total. All images in eight bands (F090W, F115W, F150W, F200W, F277W, F356W, F410M, and F444W) have been reduced adopting the latest calibration and reference files available. Data reduction has been performed using an augmented version of the official JWST pipeline, with improvements aimed at removing or mitigating defects in the raw images and improving the background subtraction and photometric accuracy. We obtain an F444W-detected multiband catalog, including all NIRCam and available Hubble Space Telescope data, adopting forced-aperture photometry on point-spread-function-matched images. The catalog is intended to enable early scientific investigations and is optimized for the study of faint galaxies; it contains 24,389 sources, with a 5 σ limiting magnitude in the F444W band ranging from 28.5 AB to 30.5 AB, as a result of the varying exposure times of the surveys that observed the field. We publicly release the reduced NIRCam images, associated multiwavelength catalog, and the code adopted for 1/ f noise removal with the aim of aiding users in familiarizing themselves with JWST NIRCam data and identifying suitable targets for follow-up observations.

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 use medium-resolution JWST/NIRSpec observations from the Cosmic Evolution Early Release Science Survey to place the first constraints on dust attenuation and star formation based on Paschen lines for a sizable sample of 63 galaxies at redshifts z = 1.0–3.1. Our analysis indicates strong correlations between the Balmer decrement, H α /H β , and line ratios that include Paschen lines (i.e., Pa α /H β , Pa β /H β , and the Paschen decrement, Pa α /Pa β ), suggesting that the former is sensitive to the overall dust obscuration toward H ii regions in high-redshift galaxies. The line ratios are used to derive nebular reddening, E ( B − V ) _neb , and star formation rates (SFRs). There is marginal evidence that the SFRs deduced from Paschen lines may exceed by ≈25% those derived from Balmer lines alone, suggesting the presence of star formation that is optically thick in Balmer lines, though deeper observations are needed to confirm this result. Using the Paschen-line constraints on the bolometric SFRs, we reevaluate the relationship between dust obscuration and UV spectral slope, and find a reddening of the UV continuum that, on average, follows the SMC extinction curve. This analysis highlights the need for deeper spectroscopy of more representative samples to evaluate nebular dust attenuation and bolometric SFRs in high-redshift galaxies, and their relationship to the reddening of the UV continuum.

REQUIEM-2D (Resolving Quiescent Magnified Galaxies with 2D Grism Spectroscopy) comprises a sample of eight massive ( ${\rm{log}}({M}_{\ast }/{M}_{\odot })\gt 10.6$ ) strongly lensed quiescent galaxies at z ∼ 2. REQUIEM-2D combines the natural magnification from strong gravitational lensing with the high-spatial-resolution grism spectroscopy of the Hubble Space Telescope through a spectrophotometric fit to study spatially resolved stellar populations. We show that quiescent galaxies in the REQUIEM-2D survey have diverse formation histories with age gradients at the 1 σ –3 σ level, including examples of (1) a younger central region supporting outside-in formation, (2) flat age gradients that show evidence for both spatially uniform early formation and inside-out quenching, and (3) regions at a fixed radial distance having different ages (such asymmetries cannot be recovered when averaging stellar population measurements azimuthally). The typical dust attenuation curve for the REQUIEM-2D galaxies is constrained to be steeper than Calzetti’s law in the UV and generally consistent with A _V < 1. Combined together and accounting for the different physical radial distances and formation timescales, we find that the REQUIEM-2D galaxies that formed earlier in the universe exhibit slow and uniform growth in their inner core, whereas the galaxies that formed later have rapid inner growth in their inner core with younger ages than the outskirts. These results challenge the currently accepted paradigm of how massive quiescent galaxies form, where the earliest galaxies are thought to form most rapidly. Significantly larger samples close to the epoch of formation with similar data quality and higher spectral resolution are required to validate this finding.

We conduct a systematic search for protocluster candidates at z ≥ 6 in the Cosmic Evolution Survey (COSMOS) field using the recently released COSMOS2020 source catalog. We select galaxies using a number of selection criteria to obtain a sample of galaxies that have a high probability of being inside a given redshift bin. We then apply overdensity analysis to the bins using two density estimators, a Weighted Adaptive Kernel estimator and a Weighted Voronoi Tessellation estimator. We have found 15 significant (>4 σ ) candidate galaxy overdensities across the redshift range 6 ≤ z ≤ 7.7. The majority of the galaxies appear to be on the galaxy main sequence at their respective epochs. We use multiple stellar-mass-to-halo-mass conversion methods to obtain a range of dark matter halo mass estimates for the overdensities in the range of ∼10 ^11 –10 ^13 M _⊙ , at the respective redshifts of the overdensities. The number and the masses of the halos associated with our protocluster candidates are consistent with what is expected from the area of a COSMOS-like survey in a standard Λ cold dark matter cosmology. Through comparison with simulation, we expect that all of the overdensities at z ≃ 6 will evolve into Virgo-/Coma-like clusters at present (i.e., with masses ∼10 ^14 –10 ^15 M _⊙ ). Compared to other overdensities identified at z ≥ 6 via narrowband selection techniques, the overdensities presented appear to have ∼10× higher stellar masses and star formation rates (SFRs). We compare the evolution in the total SFR and stellar mass content of the protocluster candidates across the redshift range 6 ≤ z ≤ 7.7 and find agreement with the total average SFR from simulations.

Star cluster formation in the early universe and its contribution to reionization remains largely unconstrained to date. Here we present JWST/NIRCam imaging of the most highly magnified galaxy known at z ∼ 6, the Sunrise arc. We identify six young massive star clusters (YMCs) with measured radii spanning from ∼20 down to ∼1 pc (corrected for lensing magnification), estimated stellar masses of ∼10 ^6–7 M _⊙ , and ages of 1–30 Myr based on SED fitting to photometry measured in eight filters extending to rest frame 7000 Å. The resulting stellar mass surface densities are higher than 1000 M _⊙ pc ^−2 (up to a few 10 ^5 M _⊙ pc ^−2 ), and their inferred dynamical ages qualify the majority of these systems as gravitationally bound stellar clusters. The star cluster ages map the progression of star formation along the arc, with two evolved systems (≳10 Myr old) followed by very young clusters. The youngest stellar clusters (1000 Å rest frame and are hosted in a 200 pc sized star-forming complex. Such a region dominates the ionizing photon production with a high efficiency $\mathrm{log}({\xi }_{\mathrm{ion}}[\mathrm{Hz}\,{\mathrm{erg}}^{-1}])\sim 25.7$ . A significant fraction of the recently formed stellar mass of the galaxy (10%–30%) occurred in these YMCs. We speculate that such sources of ionizing radiation boost the ionizing photon production efficiency, which eventually carves ionized channels that might favor the escape of Lyman continuum radiation. The survival of some of the clusters would make them the progenitors of massive and relatively metal-poor globular clusters in the local universe.

We study the spatially resolved stellar populations of 444 galaxies at 0.3 < z < 6.0 in two clusters (WHL 0137–08 and MACS 0647+70) and a blank field, combining imaging data from the Hubble Space Telescope and JWST to perform spatially resolved spectral energy distribution (SED) modeling using piXedfit . The high spatial resolution of the imaging data combined with magnification from gravitational lensing in the cluster fields allows us to resolve a large fraction of our galaxies (109) to subkiloparsec scales. At redshifts around cosmic noon and higher (2.5 ≲ z ≲ 6.0), we find mass-doubling times to be independent of radius, inferred from flat specific star formation rate (sSFR) radial profiles and similarities between the half-mass and half-SFR radii. At lower redshifts (1.5 ≲ z ≲ 2.5), a significant fraction of our star-forming galaxies shows evidence for nuclear starbursts, inferred from a centrally elevated sSFR and a much smaller half-SFR radius compared to the half-mass radius. At later epochs, we find more galaxies suppress star formation in their centers but are still actively forming stars in the disk. Overall, these trends point toward a picture of inside-out galaxy growth consistent with theoretical models and simulations. We also observe a tight relationship between the central mass surface density and global stellar mass with ∼0.38 dex scatter. Our analysis demonstrates the potential of spatially resolved SED analysis with JWST data. Future analysis with larger samples will be able to further explore the assembly of galaxy mass and the growth of their structures.

We present the results of a systematic search for candidate quiescent galaxies in the distant universe in 11 JWST fields with publicly available observations collected during the first 3 months of operations and covering an effective sky area of ∼145 arcmin ^2 . We homogeneously reduce the new JWST data and combine them with existing observations from the Hubble Space Telescope. We select a robust sample of ∼80 candidate quiescent and quenching galaxies at 3 < z < 5 using two methods: (1) based on their rest-frame UVJ colors, and (2) a novel quantitative approach based on Gaussian mixture modeling of the near-UV − U , U − V , and V − J rest-frame color space, which is more sensitive to recently quenched objects. We measure comoving number densities of massive ( M _⋆ ≥ 10 ^10.6 M _⊙ ) quiescent galaxies consistent with previous estimates relying on ground-based observations, after homogenizing the results in the literature with our mass and redshift intervals. However, we find significant field-to-field variations of the number densities up to a factor of 2–3, highlighting the effect of cosmic variance and suggesting the presence of overdensities of red quiescent galaxies at z > 3, as could be expected for highly clustered massive systems. Importantly, JWST enables the robust identification of quenching/quiescent galaxy candidates at lower masses and higher redshifts than before, challenging standard formation scenarios. All data products, including the literature compilation, are made publicly available.