Search

The unprecedented infrared spectroscopic capabilities of JWST have provided high-quality interstellar medium metallicity measurements and enabled characterization of the gas-phase mass–metallicity relation (MZR) for galaxies at z ≳ 5 for the first time. We analyze the gas-phase MZR and its evolution in a high-redshift suite of FIRE-2 cosmological zoom-in simulations at z = 5–12 and for stellar masses M _* ∼ 10 ^6 –10 ^10 M _⊙ . These simulations implement a multichannel stellar feedback model and produce broadly realistic galaxy properties, including when evolved to z = 0. The simulations predict very weak redshift evolution of the MZR over the redshift range studied, with the normalization of the MZR increasing by less than 0.01 dex as redshift decreases from z = 12 to z = 5. The median MZR in the simulations is well approximated as a constant power-law relation across this redshift range given by $\mathrm{log}(Z/{Z}_{\odot })=0.37\mathrm{log}({M}_{* }/{{\rm{M}}}_{\odot })-4.3$ . We find good agreement between our best-fit model and recent observations made by JWST at high redshift. The weak evolution of the MZR at z > 5 contrasts with the evolution at z ≲ 3, where increasing normalization of the MZR with decreasing redshift is observed and predicted by most models. The FIRE-2 simulations predict increasing scatter in the gas-phase MZR with decreasing stellar mass, in qualitative agreement with some observations.

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.

A kinematic misalignment of the stellar and gas components is a phenomenon observed in a significant fraction of galaxies. However, the underlying physical mechanisms are not well understood. A commonly proposed scenario for the formation of a misaligned component requires any preexisting gas disk to be removed, via flybys or ejective feedback from an active galactic nucleus. In this Letter, we study the evolution of a Milky Way mass galaxy in the FIREbox cosmological volume that displays a thin, counterrotating gas disk with respect to its stellar component at low redshift. In contrast to scenarios involving gas ejection, we find that preexisting gas is mainly removed via the conversion into stars in a central starburst, triggered by a merging satellite galaxy. The newly accreted, counterrotating gas eventually settles into a kinematically misaligned disk. About 4% (8 out of 182) of FIREbox galaxies with stellar masses larger than 5 × 10 ^9 M _⊙ at z = 0 exhibit gas–star kinematic misalignment. In all cases, we identify central starburst-driven depletion as the main reason for the removal of the preexisting corotating gas component, with no need for feedback from, e.g., a central active black hole. However, during the starburst, the gas is funneled toward the central regions, likely enhancing black hole activity. By comparing the fraction of misaligned discs between FIREbox and other simulations and observations, we conclude that this channel might have a non-negligible role in inducing kinematic misalignment in galaxies.

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.

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.

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.

An unprecedented array of new observational capabilities are starting to yield key constraints on models of the epoch of first light in the Universe. In this Letter we discuss the implications of the UV radiation background at cosmic dawn inferred by recent JWST observations for radio experiments aimed at detecting the redshifted 21 cm hyperfine transition of diffuse neutral hydrogen. Under the basic assumption that the 21 cm signal is activated by the Ly α photon field produced by metal-poor stellar systems, we show that a detection at the low frequencies of the EDGES and SARAS3 experiments may be expected from a simple extrapolation of the declining UV luminosity density inferred at z ≲ 14 from JWST early galaxy data. Accounting for an early radiation excess above the cosmic microwave background suggests a shallower or flat evolution to simultaneously reproduce low- and high- z current UV luminosity density constraints, which cannot be entirely ruled out, given the large uncertainties from cosmic variance and the faint-end slope of the galaxy luminosity function at cosmic dawn. Our findings raise the intriguing possibility that a high star formation efficiency at early times may trigger the onset of intense Ly α emission at redshift z ≲ 20 and produce a cosmic 21 cm absorption signal 200 Myr after the Big Bang.

Deep observations with the James Webb Space Telescope (JWST) have revealed an emerging population of red pointlike sources that could provide a link between the postulated supermassive black hole seeds and observed quasars. In this work, we present a JWST/NIRSpec spectrum from the JWST Cycle 1 UNCOVER Treasury survey of a massive accreting black hole at z = 8.50 displaying a clear broad-line component as inferred from the H β line with FWHM = 3439 ± 413 km s ^−1 , typical of the broad-line region of an active galactic nucleus (AGN). The AGN nature of this object is further supported by high ionization, as inferred from emission lines, and a point-source morphology. We compute a black hole mass of ${\mathrm{log}}_{10}({M}_{\mathrm{BH}}/{M}_{\odot })=8.17\pm 0.42$ and a bolometric luminosity of L _bol ∼ 6.6 × 10 ^45 erg s ^−1 . These values imply that our object is accreting at ∼40% of the Eddington limit. Detailed modeling of the spectral energy distribution in the optical and near-infrared, together with constraints from ALMA, indicate an upper limit on the stellar mass of ${\mathrm{log}}_{10}({M}_{* }/{M}_{\odot })\lt 8.7$ , which would lead to an unprecedented ratio of black hole to host mass of at least ∼30%. This is orders of magnitude higher compared to the local QSOs but consistent with recent AGN studies at high redshift with JWST. This finding suggests that a nonnegligible fraction of supermassive black holes either started out from massive seeds and/or grew at a super-Eddington rate at high redshift. Given the predicted number densities of high- z faint AGN, future NIRSpec observations of larger samples will allow us to further investigate galaxy–black hole coevolution in the early Universe.

In 2022 November, the James Webb Space Telescope (JWST) returned deep near-infrared images of A2744—a powerful lensing cluster capable of magnifying distant, incipient galaxies beyond it. Together with existing Hubble Space Telescope (HST) imaging, this publicly available data set opens a fundamentally new discovery space to understand the remaining mysteries of the formation and evolution of galaxies across cosmic time. In this work, we detect and measure some 60,000 objects across the 49 arcmin ^2 JWST footprint down to a 5 σ limiting magnitude of ∼30 mag in 0.″32 apertures. Photometry is performed using circular apertures on images matched to the point-spread function (PSF) of the reddest NIRCam broad band, F444W, and cleaned of bright cluster galaxies and the related intracluster light. To give an impression of the photometric performance, we measure photometric redshifts and achieve a σ _NMAD ≈ 0.03 based on known, but relatively small, spectroscopic samples. With this paper, we publicly release our HST and JWST PSF-matched photometric catalog with optimally assigned aperture sizes for easy use, along with single aperture catalogs, photometric redshifts, rest-frame colors, and individual magnification estimates. These catalogs will set the stage for efficient and deep spectroscopic follow up of some of the first JWST-selected samples in summer of 2023.

The recent UNCOVER survey with the James Webb Space Telescope (JWST) exploits the nearby cluster A2744 to create the deepest view of our Universe to date by leveraging strong gravitational lensing. In this work, we perform photometric fitting of more than 50,000 robustly detected sources out to z ∼ 15. We show the redshift evolution of stellar ages, star formation rates, and rest-frame colors across the full range of 0.2 ≲ z ≲ 15. The galaxy properties are inferred using the Prospector Bayesian inference framework using informative Prospector - β priors on the masses and star formation histories to produce joint redshift and stellar populations posteriors. Additionally, lensing magnification is performed on the fly to ensure consistency with the scale-dependent priors. We show that this approach produces excellent photometric redshifts with σ _NMAD ∼ 0.03, of a similar quality to the established photometric redshift code EAzY . In line with the open-source scientific objective of this Treasury survey, we publicly release the stellar population catalog with this paper, derived from our photometric catalog adapting aperture sizes based on source profiles. This release (the catalog and all related documentation are accessible via the UNCOVER survey web page: https://jwst-uncover.github.io/DR2.html#SPSCatalogs with a copy deposited to Zenodo at doi: https://doi.org/10.5281/zenodo.8401181 ) includes posterior moments, maximum likelihood spectra, star formation histories, and full posterior distributions, offering a rich data set to explore the processes governing galaxy formation and evolution over a parameter space now accessible by JWST.

We present visual classifications of merger-induced tidal disturbances in 143 M _* ∼ 10 ^11 M _⊙ post-starburst galaxies at z ∼ 0.7 identified in the $\mathrm{SQuIGG}\vec{L}{\rm{E}}$ Sample. This sample spectroscopically selects galaxies from the Sloan Digital Sky Survey that have stopped their primary epoch of star formation within the past ∼500 Myr. Visual classifications are performed on Hyper Suprime-Cam imaging. We compare to a control sample of mass- and redshift-matched star-forming and quiescent galaxies from the Large Early Galaxy Census and find that post-starburst galaxies are more likely to be classified as disturbed than either category. This corresponds to a factor of ${3.6}_{-1.3}^{+2.9}$ times the disturbance rate of older quiescent galaxies and ${2.1}_{-.73}^{+1.9}$ times the disturbance rate of star-forming galaxies. Assuming tidal features persist for ≲500 Myr, this suggests merging is coincident with quenching in a significant fraction of these post-starbursts. Galaxies with tidal disturbances are younger on average than undisturbed post-starburst galaxies in our sample, suggesting tidal features from a major merger may have faded over time. This may be exacerbated by the fact that, on average, the undisturbed subset is fainter, rendering low-surface-brightness tidal features harder to identify. However, the presence of 10 young (≲150 Myr since quenching) undisturbed galaxies suggests that major mergers are not the only fast physical mechanism that shut down the primary epoch of star formation in massive galaxies at intermediate redshift.

We describe a public data release of the FIRE-2 cosmological zoom-in simulations of galaxy formation (available at http://flathub.flatironinstitute.org/fire ) from the Feedback In Realistic Environments (FIRE) project. FIRE-2 simulations achieve parsec-scale resolution to explicitly model the multiphase interstellar medium while implementing direct models for stellar evolution and feedback, including stellar winds, core-collapse and Type Ia supernovae, radiation pressure, photoionization, and photoelectric heating. We release complete snapshots from three suites of simulations. The first comprises 20 simulations that zoom in on 14 Milky Way (MW)–mass galaxies, five SMC/LMC-mass galaxies, and four lower-mass galaxies including one ultrafaint; we release 39 snapshots across z = 0–10. The second comprises four massive galaxies, with 19 snapshots across z = 1–10. Finally, a high-redshift suite comprises 22 simulations, with 11 snapshots across z = 5–10. Each simulation also includes dozens of resolved lower-mass (satellite) galaxies in its zoom-in region. Snapshots include all stored properties for all dark matter, gas, and star particles, including 11 elemental abundances for stars and gas, and formation times (ages) of star particles. We also release accompanying (sub)halo catalogs, which include galaxy properties and member star particles. For the simulations to z = 0, including all MW-mass galaxies, we release the formation coordinates and an “ex situ” flag for all star particles, pointers to track particles across snapshots, catalogs of stellar streams, and multipole basis expansions for the halo mass distributions. We describe publicly available python packages for reading and analyzing these simulations.

Using an isolated Milky Way-mass galaxy simulation, we compare results from nine state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt–Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly formed stellar clump mass functions show more significant variation (difference by up to a factor of ∼3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low-density region, and between more diffusive and less diffusive schemes in the high-density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes. [ABSTRACT FROM AUTHOR]