Search

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

We present CAMELS-ASTRID, the third suite of hydrodynamical simulations in the Cosmology and Astrophysics with MachinE Learning (CAMELS) project, along with new simulation sets that extend the model parameter space based on the previous frameworks of CAMELS-TNG and CAMELS-SIMBA, to provide broader training sets and testing grounds for machine-learning algorithms designed for cosmological studies. CAMELS-ASTRID employs the galaxy formation model following the ASTRID simulation and contains 2124 hydrodynamic simulation runs that vary three cosmological parameters (Ω _m , σ _8 , Ω _b ) and four parameters controlling stellar and active galactic nucleus (AGN) feedback. Compared to the existing TNG and SIMBA simulation suites in CAMELS, the fiducial model of ASTRID features the mildest AGN feedback and predicts the least baryonic effect on the matter power spectrum. The training set of ASTRID covers a broader variation in the galaxy populations and the baryonic impact on the matter power spectrum compared to its TNG and SIMBA counterparts, which can make machine-learning models trained on the ASTRID suite exhibit better extrapolation performance when tested on other hydrodynamic simulation sets. We also introduce extension simulation sets in CAMELS that widely explore 28 parameters in the TNG and SIMBA models, demonstrating the enormity of the overall galaxy formation model parameter space and the complex nonlinear interplay between cosmology and astrophysical processes. With the new simulation suites, we show that building robust machine-learning models favors training and testing on the largest possible diversity of galaxy formation models. We also demonstrate that it is possible to train accurate neural networks to infer cosmological parameters using the high-dimensional TNG-SB28 simulation set.

The Cosmology and Astrophysics with Machine Learning Simulations (CAMELS) project was developed to combine cosmology with astrophysics through thousands of cosmological hydrodynamic simulations and machine learning. CAMELS contains 4233 cosmological simulations, 2049 N -body simulations, and 2184 state-of-the-art hydrodynamic simulations that sample a vast volume in parameter space. In this paper, we present the CAMELS public data release, describing the characteristics of the CAMELS simulations and a variety of data products generated from them, including halo, subhalo, galaxy, and void catalogs, power spectra, bispectra, Ly α spectra, probability distribution functions, halo radial profiles, and X-rays photon lists. We also release over 1000 catalogs that contain billions of galaxies from CAMELS-SAM: a large collection of N -body simulations that have been combined with the Santa Cruz semianalytic model. We release all the data, comprising more than 350 terabytes and containing 143,922 snapshots, millions of halos, galaxies, and summary statistics. We provide further technical details on how to access, download, read, and process the data at https://camels.readthedocs.io .

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

The UltraViolet imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey Fields (UVCANDELS) program provides deep Hubble Space Telescope (HST) F275W and F435W imaging over four CANDELS fields (GOODS-N, GOODS-S, COSMOS, and Extended Groth Strip). We combine this newly acquired UV imaging with existing HST imaging from CANDELS as well as existing ancillary data to obtain robust photometric redshifts and reliable estimates for galaxy physical properties for over 150,000 galaxies in the ∼430 arcmin ^2 UVCANDELS area. Here, we leverage the power of the new UV photometry to not only improve the photometric redshift measurements in these fields, but also constrain the full redshift probability distribution combining multiple redshift-fitting tools. Furthermore, using the full UV-to-IR photometric data set, we measure the galaxy physical properties by fitting templates from population synthesis models with two different parameterizations (flexible and fixed form) of the star formation histories (SFHs). Compared to the flexible SFH parameterization, we find that the fixed-form SFHs systematically underestimate the galaxy stellar masses, both at the low-mass (≲10 ^9 M _⊙ ) and high-mass (≳10 ^10 M _⊙ ) end, by as much as ∼0.5 dex. This underestimation is primarily due the limited ability of fixed-form SFH parameterization to simultaneously capture the chaotic nature of star formation in these galaxies.

We explore the cosmic evolution of the fraction of dust-obscured star formation predicted by the simba cosmological hydrodynamic simulations featuring an on-the-fly model for dust formation, evolution, and destruction. We find that up to z = 3, our results are broadly consistent with previous observational results of little to no evolution in obscured star formation. However, at z > 3 we find strong evolution at fixed galaxy stellar mass toward greater amounts of obscured star formation, in tension with high-redshift observations. We explain the trend of increasing obscuration at higher redshifts by evolving star-dust geometry, as the dust-to-stellar mass ratios remain relatively constant across cosmic time. We additionally see that at a fixed redshift, more massive galaxies have a higher fraction of their star formation obscured, which is explained by increased dust-to-stellar mass ratios at higher stellar masses. Finally, we estimate the contribution of dust-obscured star formation to the total star formation rate budget and find that the dust-obscured star formation history peaks around z ∼ 2−3, and becomes subdominant at z ≳ 5. The dominance of obscured star formation at redshifts z ≲ 4 is consistent with our results for the evolution of the obscured star formation fraction at fixed stellar mass to higher values at higher redshift because there exist fewer massive, heavily obscured galaxies at high redshift.

The ratios of strong rest-frame optical emission lines are the dominant indicators of metallicities in high-redshift galaxies. Since typical strong-line-based metallicity indicators are calibrated on auroral lines at z = 0, their applicability for galaxies in the distant Universe is unclear. In this paper, we make use of mock emission-line data from cosmological simulations to investigate the calibration of rest-frame optical emission lines as metallicity indicators at high redshift. Our model, which couples the simba cosmological galaxy formation simulation with cloudy photoionization calculations, includes contributions from H ii regions, post-asymptotic-giant-branch stars, and diffuse ionized gas (DIG). We find mild redshift evolution in the 12 indicators that we study, which implies that the dominant physical properties that evolve in our simulations do have a discernible impact on the metallicity calibrations at high redshifts. When comparing our calibrations with high-redshift auroral line observations from the James Webb Space Telescope, we find a slight offset between our model results and the observations and find that a higher ionization parameter at high redshifts can be one of the possible explanations. We explore the physics that drives the shapes of strong-line metallicity relationships and propose calibrations for hitherto unexplored low-metallicity regimes. Finally, we study the contribution of DIG to total line fluxes. We find that the contribution of DIG increases with metallicity at z ∼ 0 for singly ionized oxygen and sulfur lines and can be as high as 70%, making it crucial to include their contribution when modeling nebular emission.

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 conducted an investigation of the Coma cluster of galaxies by running a series of constrained hydrodynamic simulations with GIZMO-SIMBA and GADGET-3 based on initial conditions reconstructed from the SDSS survey volume in the ELUCID project. We compared simulation predictions and observations for galaxies, intracluster medium (ICM) and intergalactic medium (IGM) in and around the Coma cluster to constrain galaxy formation physics. Our results demonstrate that this type of constrained investigation allows us to probe in more detail the implemented physical processes, because the comparison between simulations and observations is free of cosmic variance and hence can be conducted in a “one-to-one” manner. We found that an increase in the earlier star formation rate and the supernova feedback of the original GIZMO-SIMBA model is needed to match observational data on stellar, interstellar medium, and ICM metallicity. The simulations without active galactic nucleus (AGN) feedback can well reproduce the observational ICM electron density, temperature, and entropy profiles, ICM substructures, and the IGM temperature–density relation, while the ones with AGN feedback usually fail. However, one requires something like AGN feedback to reproduce a sufficiently large population of quiescent galaxies, particularly in low-density regions. The constrained simulations of the Coma cluster thus provide a test bed to understand processes that drive galaxy formation and evolution.

The massive end of the gas-phase mass–metallicity relation (MZR) is a sensitive probe of active galactic nuclei (AGN) feedback that is a crucial but highly uncertain component of galaxy evolution models. In this paper, we extend the z ∼ 0.7 MZR by ∼0.5 dex up to log ( M _⋆ / M _⊙ ) ∼ 11.1. We use extremely deep VLT VIMOS spectra from the Large Early Galaxy Astrophysics Census (LEGA-C) survey to measure metallicities for 145 galaxies. The LEGA-C MZR matches the normalization of the z ∼ 0.8 DEEP2 MZR where they overlap, so we combine the two to create an MZR spanning from 9.3 to 11.1 log ( M _⋆ / M _⊙ ). The LEGA-C+DEEP2 MZR at z ∼ 0.7 is offset to slightly lower metallicities (0.05–0.13 dex) than the z ∼ 0 MZR, but it otherwise mirrors the established power-law rise at low/intermediate stellar masses and asymptotic flattening at high stellar masses. We compare the LEGA-C+DEEP2 MZR to the MZR from two cosmological simulations (IllustrisTNG and SIMBA), which predict qualitatively different metallicity trends for high-mass galaxies. This comparison highlights that our extended MZR provides a crucial observational constraint for galaxy evolution models in a mass regime where the MZR is very sensitive to choices about the implementation of AGN feedback.

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

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

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.

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.

Active galactic nuclei (AGNs) feedback models are generally calibrated to reproduce galaxy observables such as the stellar mass function and the bimodality in galaxy colors. We use variations of the AGN feedback implementations in the IllustrisTNG (TNG) and Simba cosmological hydrodynamic simulations to show that the low-redshift Ly α forest can provide constraints on the impact of AGN feedback. We show that TNG overpredicts the number density of absorbers at column densities N _HI < 10 ^14 cm ^−2 compared to data from the Cosmic Origins Spectrograph (in agreement with previous work), and we demonstrate explicitly that its kinetic feedback mode, which is primarily responsible for galaxy quenching, has a negligible impact on the column density distribution (CDD) of absorbers. In contrast, we show that the fiducial Simba model, which includes AGN jet feedback, is the preferred fit to the observed CDD of the z = 0.1 Ly α forest across 5 orders of magnitude in column density. We show that the Simba results with jets produce a quantitatively better fit to the observational data than the Simba results without jets, even when the ultraviolet background is left as a free parameter. AGN jets in Simba are high speed, collimated, weakly interacting with the interstellar medium (via brief hydrodynamic decoupling), and heated to the halo virial temperature. Collectively these properties result in stronger long-range impacts on the intergalactic medium when compared to TNG’s kinetic feedback mode, which drives isotropic winds with lower velocities at the galactic radius. Our results suggest that the low-redshift Ly α forest provides plausible evidence for long-range AGN jet feedback.

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.

Current models of galaxy formation require strong feedback from active galactic nuclei (AGN) to explain the observed lack of star formation in massive galaxies since z ≈ 2, but direct evidence of this energy input is limited. We use the SIMBA cosmological galaxy formation simulations to assess the ability of thermal Sunyaev–Zel’dovich (tSZ) measurements to provide such evidence, by mapping the pressure structure of the circumgalactic medium around massive z ≈ 0.2–1.5 galaxies. We undertake a stacking approach to calculate the total tSZ signal and its radial profile in simulations with varying assumptions of AGN feedback, and we assess its observability with current and future telescopes. By convolving our predictions with the 2.′1 beam of the Atacama Cosmology Telescope, we show that current observations at z ≈ 1 are consistent with SIMBA’s fiducial treatment of AGN feedback and inconsistent with SIMBA models without feedback. At z ≈ 0.5, observational signals lie between SIMBA run with and without AGN feedback, suggesting AGN in SIMBA may inject too much energy at late times. By convolving our data with a 9.″5 beam corresponding to the TolTEC camera on the Large Millimeter Telescope Alfonso Serrano, we predict a unique profile for AGN feedback that can be distinguished with future higher-resolution measurements. Finally, we explore a novel approach to quantify the nonspherically symmetric features surrounding our galaxies by plotting radial profiles representing the component of the stack with m-fold symmetry.

We present the Cosmic Sands suite of cosmological zoom-in simulations based on the simba galaxy formation model in order to study the buildup of the first massive and dusty galaxies in the early universe. Residing in the most massive halos, we find that the compact proto-massive galaxies undergo nearly continuous mergers with smaller subhalos, boosting star formation rates (SFRs) and the buildup of stellar mass. The galaxies are already appreciably chemically evolved by z = 7, with modeled dust masses comparable to those inferred from observations in the same epoch, except for the most extreme systems. We track gas accretion onto the galaxies to understand how extreme SFRs can be sustained by these early systems. We find that smooth gas accretion can maintain SFRs above 250 M _⊙ yr ^−1 , but to achieve SFRs that boost galaxies well above the main sequence, a larger perturbation like a gas-rich major merger is necessary to trigger a starburst episode. Post-processing the Cosmic Sands simulations with dust RT, we find that, while the infrared luminosities of the most-dust-rich galaxies are comparable to local ULIRGs, they are substantially dimmer than classical z = 2 submillimeter galaxies. We end with a discussion on the possible reasons for this discrepancy at the highest masses and the future work we intend to carry out to study the chemical enrichment of the earliest dusty galaxies.

We present observations of CO(3−2) in 13 main-sequence z = 2.0–2.5 star-forming galaxies at $\mathrm{log}({M}_{* }/{M}_{\odot })=10.2\mbox{--}10.6$ that span a wide range in metallicity (O/H) based on rest-optical spectroscopy. We find that ${L}_{\mathrm{CO}(3-2)}^{{\prime} }$ /SFR decreases with decreasing metallicity, implying that the CO luminosity per unit gas mass is lower in low-metallicity galaxies at z ∼ 2. We constrain the CO-to-H _2 conversion factor ( α _CO ) and find that α _CO inversely correlates with metallicity at z ∼ 2. We derive molecular gas masses ( M _mol ) and characterize the relations among M _* , SFR, M _mol , and metallicity. At z ∼ 2, M _mol increases and the molecular gas fraction ( M _mol / M _* ) decreases with increasing M _* , with a significant secondary dependence on SFR. Galaxies at z ∼ 2 lie on a near-linear molecular KS law that is well-described by a constant depletion time of 700 Myr. We find that the scatter about the mean SFR− M _* , O/H− M _* , and M _mol − M _* relations is correlated such that, at fixed M _* , z ∼ 2 galaxies with larger M _mol have higher SFR and lower O/H. We thus confirm the existence of a fundamental metallicity relation at z ∼ 2, where O/H is inversely correlated with both SFR and M _mol at fixed M _* . These results suggest that the scatter of the z ∼ 2 star-forming main sequence, mass–metallicity relation, and M _mol – M _* relation are primarily driven by stochastic variations in gas inflow rates. We place constraints on the mass loading of galactic outflows and perform a metal budget analysis, finding that massive z ∼ 2 star-forming galaxies retain only 30% of metals produced, implying that a large mass of metals resides in the circumgalactic medium.

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.