Search

The evolution of galaxy sizes in different wavelengths provides unique insights on galaxy build-up across cosmic epochs. Such measurements can now finally be done at $z>3$ thanks to the exquisite spatial resolution and multi-wavelength capability of the JWST. With the public data from the CEERS, PRIMER-UDS, and PRIMER-COSMOS surveys, we measure the sizes of $\sim 3500$ star-forming galaxies at $3 \leqslant z<9$, in 7 NIRCam bands using the multi-wavelength model fitting code GalfitM. The size-mass relation is measured in four redshift bins, across all NIRCam bands. We find that, the slope and intrinsic scatter of the rest-optical size-mass relation are constant across this redshift range and consistent with previous HST-based studies at low-z. When comparing the relations across different wavelengths, the average rest-optical and rest-UV relations are consistent with each other up to $z=6$, but the intrinsic scatter is largest in rest-UV wavelengths compared to rest-optical and redder bands. This behaviour is independent of redshift and we speculate that it is driven by bursty star-formation in $z>4$ galaxies. Additionally, for $3\leqslant z<4$ star-forming galaxies at $\rm M_* > 10^{10} M_{\odot}$, we find smaller rest-$\rm 1\rm\,\mu m$ sizes in comparison to rest-optical (and rest-UV) sizes, suggestive of colour gradients. When comparing to simulations, we find agreement over $\rm M_* \approx 10^{9} - 10^{10} M_{\odot}$ but beyond this mass, the observed size-mass relation is significantly steeper. Our results show the power of JWST/NIRCam to provide new constraints on galaxy formation models., Comment: 14 pages, 9 figures, submitted to A&A

We present the formation histories of 19 massive ($\gtrsim3\times10^{10}\text{M}_\odot$) quiescent galaxy candidates at $z\sim3.0-4.5$ observed using JWST/NIRSpec. This completes the spectroscopic confirmation of the 24 $K$-selected quiescent galaxy sample from the ZFOURGE and 3DHST surveys \citep{Schreiber2018}. Utilizing Prism $1-5\mu$m spectroscopy, we confirm that all 12 sources that eluded confirmation by ground-based spectroscopy lie at $z>3$, resulting in a spectroscopically confirmed number density of $\sim1.4\times10^{-5}\text{Mpc}^{-3}$ between $z\sim3-4$. Rest-frame $U-V$ vs $V-J$ color selections show high effectiveness in identifying quiescent galaxies, with a purity of $\sim90\%$. Our analysis shows that parametric star-formation histories (SFHs) from FAST++ and binned SFHs from Prospector on average yield consistent results, revealing diverse formation and quenching times. The oldest galaxy formed $\sim6\times10^{10}\text{M}_\odot$ by $z\sim10$ and has been quiescent for over 1 Gyr at $z\sim3.2$. We detect two galaxies with ongoing star formation and six with active galactic nuclei (AGN). We demonstrate that the choice of stellar population models, stellar libraries, wavelength range, and nebular or AGN contributions does not significantly affect the derived average SFHs of the galaxies. The assumed SFH prior, however, influences the star formation rate at early times, where spectral diagnostic power is limited. Simulated $z\sim3$ quiescent galaxies from IllustrisTNG, SHARK, and Magneticum broadly match the average SFHs of the observed sample but struggle to capture the full diversity, particularly at early stages. Our results emphasize the need for mechanisms that rapidly build stellar mass and quench star formation within the first billion years of the Universe., Comment: Submitted to ApJ. Comments welcome

We present the PANORAMIC survey, a pure parallel extragalactic imaging program with NIRCam observed during JWST Cycle 1. The survey obtained $\sim$530 sq arcmin of NIRCam imaging from 1-5$\mu$m, totaling $\sim$192 hours of science integration time. This represents the largest on-sky time investment of any Cycle 1 GO extragalactic NIRCam imaging program by nearly a factor of 2. The survey includes $\sim$432 sq arcmin of novel sky area not yet observed with JWST using at least $6$ NIRCam broad-band filters, increasing the existing area covered by similar Cycle 1 data by $\sim$60%. 70 square arcmin was also covered by a 7th filter (F410M). A fraction of PANORAMIC data ($\sim$200 sq arcmin) was obtained in or around extragalactic deep-fields, enhancing their legacy value. Pure parallel observing naturally creates a wedding cake survey with both wide and ultra-deep tiers, with 5$\sigma$ point source depths at F444W ranging from 27.8-29.4 (ABmag), and with minimized cosmic variance. The 6+ filter observing setup yields remarkably good photometric redshift performance, achieving similar median scatter and outlier fraction as CANDELS ($\sigma_{\rm NMAD}\sim0.07$; $\eta\sim0.2$), which enables a wealth of science across redshift without the need for followup or ancillary data. We overview the proposed survey, the data obtained as part of this program, and document the science-ready data products in the first data release. PANORAMIC has delivered wide-area and deep imaging with excellent photometric performance, demonstrating that pure parallel observations with JWST are a highly efficient observing mode that is key to acquiring a complete picture of galaxy evolution from rare bright galaxies to fainter, more abundant sources at all redshifts., Comment: Submitted to AAS Journals; comments welcome. Initial data release available at https://panoramic-jwst.github.io/

Dwarf galaxies hold the key to crucial frontiers of astrophysics, however, their faintness renders spectroscopy challenging. Here we present the JWST Cycle 2 survey, All the Little Things (ALT, PID 3516), which is designed to seek late-forming Pop III stars and the drivers of reionization at $z\sim6-7$. ALT has acquired the deepest NIRCam grism spectroscopy yet (7-27 hr), at JWST's most sensitive wavelengths (3-4 $\mu$m), covering the powerful lensing cluster Abell 2744. Over the same 30 arcmin$^2$, ALT's ultra-deep F070W+F090W imaging ($\sim$30 mag) enables selection of very faint sources at $z>6$. We demonstrate the success of ALT's novel ``butterfly" mosaic to solve spectral confusion and contamination, and introduce the ``Allegro" method for emission line identification. By collecting spectra for every source in the field of view, ALT has measured precise ($R\sim1600$) redshifts for 1630 sources at $z=0.2-8.5$. This includes one of the largest samples of distant dwarf galaxies: [1015, 475, 50] sources less massive than the SMC, Fornax, and Sculptor with $\log(M_{*}/M_{\odot})<$[8.5, 7.5, 6.5]. We showcase ALT's discovery space with: (i) spatially resolved spectra of lensed clumps in galaxies as faint as $M_{\rm{UV}}\sim-15$; (ii) large-scale clustering -- overdensities at $z$=[2.50, 2.58, 3.97, 4.30, 5.66, 5.77, 6.33] hosting massive galaxies with striking Balmer breaks; (iii) small-scale clustering -- a system of satellites around a Milky Way analog at $z\sim6$; (iv) spectroscopically confirmed multiple images that help constrain the lensing model underlying all science in this legacy field; (v) sensitive star-formation maps based on dust-insensitive tracers such as Pa$\alpha$; (vi) direct spectroscopic discovery of rare sources such as AGN with ionized outflows. These results provide a powerful proof of concept for how grism surveys maximize the potential of strong lensing fields., Comment: Submitted to the Open Journal of Astrophysics. Figs. 6 and 10 illustrate the quality of the spectra and imaging, while Fig. 12 summarizes the yield of the survey. Comments warmly welcomed and greatly appreciated. The ALT DR1 catalog is available at https://zenodo.org/records/13871850

We present the properties of a massive, large, dusty, metal-rich, star-forming galaxy at z_spec=6.73. GOODSN-100182 was observed with JWST/NIRSpec as part of the AURORA survey, and is also covered by public multi-wavelength HST and JWST imaging. While the large mass of GOODSN-100182 (~10^10 M_sun) was indicated prior to JWST, NIRCam rest-optical imaging now reveals the presence of an extended disk (r_eff~1.5 kpc). In addition, the NIRSpec R~1000 spectrum of GOODSN-100182 includes the detection of a large suite of rest-optical nebular emission lines ranging in wavelength from [OII]3727 up to [NII]6583. The ratios of Balmer lines suggest significant dust attenuation (E(B-V)_gas=0.40+0.10/-0.09), consistent with the red rest-UV slope inferred for GOODSN-100182 (beta=-0.50+/-0.09). The star-formation rate based on dust-corrected H-alpha emission is log(SFR(H-alpha)/ M_sun/yr)=2.02+0.13/-0.14, well above the z~7 star-forming main sequence in terms of specific SFR. Strikingly, the ratio of [NII]6583/H-alpha emission suggests almost solar metallicity, as does the ratio ([OIII]5007/H-beta)/([NII]6583/H-alpha) and the detection of the faint [FeII]4360 emission feature, whereas the [OIII]5007/[OII]3727 ratio suggests roughly 50% solar metallicity. Overall, the excitation and ionization properties of GOODSN-100182 more closely resemble those of typical star-forming galaxies at z~2-3 rather than z~7. Based on public spectroscopy of the GOODS-N field, we find that GOODSN-100182 resides within a significant galaxy overdensity, and is accompanied by a spectroscopically-confirmed neighbor galaxy. GOODSN-100182 demonstrates the existence of mature, chemically-enriched galaxies within the first billion years of cosmic time, whose properties must be explained by galaxy formation models., Comment: 16 pages, 13 figures, submitted to ApJ

We report the first successful ALMA follow-up observations of a secure $z > 10$ JWST-selected galaxy, by robustly detecting ($6.6\sigma$) the [OIII]$_{88\mu m}\,$ line in JADES-GS-z14-0 (hereafter GS-z14). The ALMA detection yields a spectroscopic redshift of $z=14.1793\pm0.0007$, and increases the precision on the prior redshift measurement of $z=14.32_{-0.20}^{+0.08}$ from NIRSpec by $\gtrsim$180$\times$. Moreover, the redshift is consistent with that previously determined from a tentative detection ($3.6\sigma$) of CIII]$_{1907,1909}$ ($z=14.178\pm0.013$), solidifying the redshift determination via multiple line detections. We measure a line luminosity of $L_\mathrm{[OIII]88} = (2.1 \pm 0.5)\times10^8\,L_\odot$, placing GS-z14 at the lower end, but within the scatter of, the local $L_\mathrm{[OIII]88}$-star formation rate relation. No dust continuum from GS-z14 is detected, suggesting an upper limit on the dust-to-stellar mass ratio of $< 2 \times 10^{-3}$, consistent with dust production from supernovae with a yield $y_d < 0.3\,M_\odot$. Combining a previous JWST/MIRI photometric measurement of the [OIII]$\lambda\lambda$4959,5007$\mathrm{\mathring{A}}$ and H$\beta$ lines with Cloudy models, we find GS-z14 to be surprisingly metal-enriched ($Z\sim0.05 - 0.2\,Z_\odot$) a mere $300\,\mathrm{Myr}$ after the Big Bang. The detection of a bright oxygen line in GS-z14 thus reinforces the notion that galaxies in the early Universe undergo rapid evolution., Comment: 11 Pages, 6 Figures, 1 Table. Submitted to The Astrophysical Journal

The H{\alpha} nebular emission line is an optimal tracer for recent star formation in galaxies. With the advent of JWST, this line has recently become observable at z>3 for the first time. We present a catalog of 1013 H{\alpha} emitters at 3.73 obtained based purely on spectroscopic data, robustly tracing galaxy star formation rates (SFRs) beyond the peak of the cosmic star formation history. We compare our results with theoretical predictions from three different simulations and find good agreement at z~4-6. The UV LFs of this spectroscopically-confirmed sample are in good agreement with pre-JWST measurements obtained with photometrically-selected objects. Finally, we derive SFR functions and integrate these to compute the evolution of the cosmic star-formation rate densities across z~4-6, finding values in good agreement with recent UV estimates from Lyman-break galaxies, which imply a continuous decrease in SFR density by a factor of 3x over z~4 to z~6. Our work shows the power of NIRCam grism observations to efficiently provide new tests for early galaxy formation models based on emission line statistics., Comment: 17 pages, 14 figures

We present the Red Unknowns: Bright Infrared Extragalactic Survey (RUBIES), providing JWST/NIRSpec spectroscopy of red sources selected across ~150 arcmin$^2$ from public JWST/NIRCam imaging in the UDS and EGS fields. RUBIES novel observing strategy offers a well-quantified selection function: the survey is optimised to reach high (>70%) completeness for bright and red (F150W-F444W>2) sources that are very rare. To place these rare sources in context, we simultaneously observe a reference sample of the 23$ using only the G395M disperser. The RUBIES data reveal a highly diverse population of red sources that span a broad redshift range ($z_{spec}\sim1-9$), with photometric redshift scatter and outlier fraction that are 3 times higher than for similarly bright sources that are less red. This diversity is not apparent from the photometric SEDs. Only spectroscopy reveals that the SEDs encompass a mixture of galaxies with dust-obscured star formation, extreme line emission, a lack of star formation indicating early quenching, and luminous active galactic nuclei. As a first demonstration of our broader selection function we compare the stellar masses and rest-frame U-V colours of the red sources and our reference sample: red sources are typically more massive ($M_*\sim10^{10-11.5} M_\odot$) across all redshifts. However, we find that the most massive systems span a wide range in U-V colour. We describe our data reduction procedure and data quality, and publicly release the reduced RUBIES data and vetted spectroscopic redshifts of the first half of the survey through the DJA., Comment: 21 pages, 16 figures; submitted to A&A

We report the spectroscopic discovery of a massive quiescent galaxy at $z_{\rm spec}=7.29\pm0.01$, just $\sim700\,$Myr after the Big Bang. RUBIES-UDS-QG-z7 was selected from public JWST/NIRCam and MIRI imaging from the PRIMER survey and observed with JWST/NIRSpec as part of RUBIES. The NIRSpec/PRISM spectrum reveals one of the strongest Balmer breaks observed thus far at $z>6$, no emission lines, but tentative Balmer and Ca absorption features, as well as a Lyman break. Simultaneous modeling of the NIRSpec/PRISM spectrum and NIRCam and MIRI photometry (spanning $0.9-18\,\mu$m) shows that the galaxy formed a stellar mass of log$(M_*/M_\odot)=10.23^{+0.04}_{-0.04}$ in a rapid $\sim 100-200\,$Myr burst of star formation at $z\sim8-9$, and ceased forming stars by $z\sim8$ resulting in $\log \rm{sSFR/yr}^{-1}<-10$. We measure a small physical size of $209_{-24}^{+33}\,{\rm pc}$, which implies a high stellar mass surface density within the effective radius of $\log(\Sigma_{*,\rm e}/{\rm M_\odot\,kpc}^{-2})=10.85_{-0.12}^{+0.11}$ comparable to the densities measured in quiescent galaxies at $z\sim2-5$. The 3D stellar mass density profile of RUBIES-UDS-QG-z7 is remarkably similar to the central densities of local massive ellipticals, suggesting that at least some of their cores may have already been in place at $z>7$. The discovery of RUBIES-UDS-QG-z7 has strong implications for galaxy formation models: the estimated number density of quiescent galaxies at $z\sim7$ is $>100\times$ larger than predicted from any model to date, indicating that quiescent galaxies have formed earlier than previously expected., Comment: 21 pages, 7 figures, submitted to ApJ

We use JWST/NIRSpec observations from the Assembly of Ultradeep Rest-optical Observations Revealing Astrophysics (AURORA) survey to constrain the shape of the nebular attenuation curve of a star-forming galaxy at z=4.41, GOODSN-17940. We utilize 11 unblended HI recombination lines to derive the attenuation curve spanning optical to near-infrared wavelengths (3751-9550 \r{A}). We then leverage a high-S/N spectroscopic detection of the rest-frame ultraviolet continuum in combination with rest-UV photometric measurements to constrain the shape of the curve at ultraviolet wavelengths. While this UV constraint is predominantly based on stellar emission, the large measured equivalent widths of H$\alpha$ and H$\beta$ indicate that GOODSN-17940 is dominated by an extremely young stellar population <10 Myr in age such that the UV stellar continuum experiences the same attenuation as the nebular emission. The resulting combined nebular attenuation curve spans 1400-9550 \r{A} and has a shape that deviates significantly from commonly assumed dust curves in high-redshift studies. Relative to the Milky Way, SMC, and Calzetti curves, the new curve has a steeper slope at long wavelengths ($\lambda>5000$ \r{A}) while displaying a similar slope across blue-optical wavelengths ($\lambda=3750-5000$ \r{A}). In the ultraviolet, the new curve is shallower than the SMC and Calzetti curves and displays no significant 2175 \r{A} bump. This work demonstrates that the most commonly assumed dust curves are not appropriate for all high-redshift galaxies. These results highlight the ability to derive nebular attenuation curves for individual high-redshift sources with deep JWST/NIRSpec spectroscopy, thereby improving the accuracy of physical properties inferred from nebular emission lines., Comment: 21 pages, 6 figures, 1 table, submitted to ApJ

We present the design and observations of low resolution JWST/NIRSpec PRISM spectroscopy from the Ultradeep NIRSpec and NIRCam ObserVations before the Epoch of Reionization (UNCOVER) Cycle 1 JWST Treasury program. Targets are selected using JWST/NIRCam photometry from UNCOVER and other programs, and cover a wide range of categories and redshifts to ensure the legacy value of the survey. These categories include the first galaxies at $z\gtrsim10$, faint galaxies during the Epoch of Reionization ($z\sim6-8$), high redshift AGN ($z\gtrsim6$), Population III star candidates, distant quiescent and dusty galaxies ($1\lesssim z \lesssim 6$), and filler galaxies sampling redshift--color--magnitude space from $z\sim 0.1-13$. Seven NIRSpec MSA masks across the extended Abell 2744 cluster were observed, along with NIRCam parallel imaging in 8 filters (F090W, F115W, F150W, F200W, F277W, F356W, F410M, F444W, F480M) over a total area of ~26 arcmin$^2$, overlapping existing HST coverage from programs including the Hubble Frontier Fields and BUFFALO. We successfully observed 553 objects down to $m_{\mathrm{F444W}}\sim30\mathrm{AB}$, and by leveraging mask overlaps, we reach total on-target exposure times ranging from 2.4-16.7h. We demonstrate the success rate and distribution of confirmed redshifts, and also highlight the rich information revealed by these ultradeep spectra for a subset of our targets. An updated lens model of Abell 2744 is also presented, including 14 additional spectroscopic redshifts and finding a total cluster mass of $M_{\mathrm{SL}}=(2.1\pm0.3)\times10^{15}\,\mathrm{M}_{\odot}$. We publicly release reduced 1D and 2D spectra for all objects observed in Summer 2023 along with a spectroscopic redshift catalog and the updated lens model of the cluster (https://jwst-uncover.github.io/DR4.html)., Comment: 19 pages, 9 figures, 4 tables, submitted to ApJ, comments welcome! Data available at: https://jwst-uncover.github.io/DR4.html (v2: figure format correction)

We study, from both a theoretical and observational perspective, the physical origin and spectroscopic impact of extreme nebular emission in high-redshift galaxies. The nebular continuum, which can appear during extreme starbursts, is of particular importance as it tends to redden UV slopes and has a significant contribution to the UV luminosities of galaxies. Furthermore, its shape can be used to infer the gas density and temperature of the ISM. First, we provide a theoretical background, showing how different stellar populations (SPS models, IMFs, and stellar temperatures) and nebular conditions impact observed galaxy spectra. We demonstrate that, for systems with strong nebular continuum emission, 1) UV fluxes can increase by up to 0.7~magnitudes (or more in the case of hot/massive stars) above the stellar continuum, which may help reconcile the surprising abundance of bright high-redshift galaxies and the elevated UV luminosity density at $z>10$, 2) at high gas densities, UV slopes can redden from $\beta\lesssim-2.5$ to $\beta\sim-1$, 3) observational measurements of $\xi_{ion}$ are grossly underestimated, and 4) UV downturns from two-photon emission can masquerade as DLAs. Second, we present a dataset of 58 galaxies observed with NIRSpec on JWST at $2.510$ galaxies, finding that UV slopes and UV downturns are in some cases redder or steeper than expected from SPS models, which may hint at more exotic (e.g. hotter/more massive stars or AGN) ionizing sources., Comment: 26 pages, 18 figures, submitted to The Open Journal of Astrophysics

The mass-metallicity relation (MZR) provides crucial insights into the baryon cycle in galaxies and provides strong constraints on galaxy formation models. We use JWST NIRSpec observations from the UNCOVER program to measure the gas-phase metallicity in a sample of eight galaxies during the epoch of reionization at $z=6-8$. Thanks to strong lensing of the galaxy cluster Abell 2744, we are able to probe extremely low stellar masses between $10^{6}$ and $10^{8} M_{\odot}$. Using strong lines diagnostics and the most recent JWST calibrations, we derive extremely-low oxygen abundances ranging from 12+log(O/H)=6.7 to 7.8. By combining this sample with more massive galaxies at similar redshifts, we derive a best-fit relation of 12+{\rm log(O/H)}=$0.39_{-0.02}^{+0.02} \times$ log(\mstar) $+ 4.52_{-0.17}^{+0.17}$, which is steeper than determinations at $z \sim 3$. Our results show a clear redshift evolution in the overall normalization of the relation, galaxies at higher redshift having significantly lower metallicities at a given mass. A comparison with theoretical models provides important constraints on which physical processes, such as metal mixing, star formation or feedback recipes, are important in reproducing the observations. Additionally, these galaxies exhibit star formation rates that are higher by a factor of a few to tens compared to extrapolated relations at similar redshifts or theoretical predictions of main-sequence galaxies, pointing to a recent burst of star formation. All these observations are indicative of highly stochastic star formation and ISM enrichment, expected in these low-mass systems, suggesting that feedback mechanisms in high-$z$ dwarf galaxies might be different from those in place at higher masses., Comment: Submitted to ApJL

Galaxy clusters are the most massive, gravitationally-bound structures in the Universe, emerging through hierarchical structure formation of large-scale dark matter and baryon overdensities. Early galaxy ``proto-clusters'' are believed to be important physical drivers of the overall cosmic star-formation rate density and serve as ``hotspots'' for the reionization of the intergalactic medium. Our understanding of the formation of these structures at the earliest cosmic epochs is, however, limited to sparse observations of their galaxy members, or based on phenomenological models and cosmological simulations. Here we report the detection of a massive neutral, atomic hydrogen (HI) gas reservoir permeating a galaxy proto-cluster at redshift $z=5.4$, observed one billion years after the Big Bang. The presence of this cold gas is revealed by strong damped Lyman-$\alpha$ absorption features observed in several background galaxy spectra taken with JWST/NIRSpec in close on-sky projection. While overall the sightlines probe a large range in HI column densities, $N_{\rm HI} = 10^{21.7}-10^{23.5}$ cm$^{-2}$, they are similar across nearby sightlines, demonstrating that they probe the same dense, neutral gas. This observation of a massive, large-scale overdensity of cold neutral gas challenges current large-scale cosmological simulations and has strong implications for the reionization topology of the Universe., Comment: Submitted

Recent observations with the James Webb Space Telescope (JWST) have uncovered unexpectedly high cosmic star formation activity in the early Universe, mere hundreds of millions of years after the Big Bang. These observations are often understood to reflect an evolutionary shift in star formation efficiency (SFE) caused by changing galactic conditions during these early epochs. We present FIREbox-HR, a high-resolution, cosmological hydrodynamical simulation from the Feedback in Realistic Environments project, which offers insights into the SFE of galaxies during the first billion years of cosmic time. FIREbox-HR re-simulates the cosmic volume (L = 22.1 cMpc) of the original FIREbox run with eight times higher mass resolution (m_b ~ 7800 M_sun), but with identical physics, down to z ~ 6. FIREbox-HR predicts ultraviolet (UV) luminosity functions in good agreement with available observational data. The simulation also successfully reproduces the observed cosmic UV luminosity density at z ~ 6 - 14, demonstrating that relatively high star formation activity in the early Universe is a natural outcome of the baryonic processes encoded in the FIRE-2 model. According to FIREbox-HR, the SFE - halo mass relation for intermediate mass halos (M_halo ~ 10^9 - 10^11 M_sun) does not significantly evolve with redshift and is only weakly mass-dependent. These properties of the SFE - halo mass relation lead to a larger contribution from lower mass halos at higher z, driving the gradual evolution of the observed cosmic UV luminosity density. A theoretical model based on the SFE - halo mass relation inferred from FIREbox-HR allows us to explore implications for galaxy evolution. Future observations of UV faint galaxies at z > 12 will provide an opportunity to further test these predictions and deepen our understanding of star formation during Cosmic Dawn., Comment: 26 pages, 14 figures, 5 tables, submitted to MNRAS, comments welcome

We present results on the emission-line properties of z=1.4-7.5 star-forming galaxies in the Assembly of Ultradeep Rest-optical Observations Revealing Astrophysics (AURORA) Cycle 1 JWST/NIRSpec program. Based on its depth, continuous wavelength coverage from 1--5 microns, and medium spectral resolution (R~1000), AURORA includes detections of a large suite of nebular emission lines spanning a broad range in rest wavelength. We investigate the locations of AURORA galaxies in multiple different emission-line diagrams, including traditional "BPT" diagrams of [OIII]/Hbeta vs. [NII]/Halpha, [SII]/Halpha, and [OI]/Halpha, and the "ionization-metallicity" diagram of [OIII]/[OII] (O32) vs. ([OIII]+[OII])/Hbeta (R23). We also consider a bluer rest-frame "ionization-metallicity" diagram introduced recently to characterize z>10 galaxies: [NeIII]/[OII] vs. ([NeIII]+[OII])/Hdelta; as well as longer-wavelength diagnostic diagrams extending into the rest-frame near-IR: [OIII]/Hbeta vs. [SIII]/[SII] (S32); and HeI/Pagamma and [SIII]/Pagamma vs. [FeII]/Pabeta. With a significant boost in signal-to-noise and large, representative samples of individual galaxy detections, the AURORA emission-line diagrams presented here definitively confirm a physical picture in which chemically-young, alpha-enhanced, massive stars photoionize the ISM in distant galaxies with a harder ionizing spectrum at fixed nebular metallicity than in their z~0 counterparts. We also uncover previously unseen evolution prior to z~2 in the [OIII]/Hbeta vs. [NII]/Halpha diagram, which motivates deep NIRSpec observations at even higher redshift. Finally, we present the first statistical sample of rest-frame near-IR emission-line diagnostics in star-forming galaxies at high redshift. In order to truly interpret rest-frame near-IR line ratios including [FeII], we must obtain better constraints on dust depletion in the high-redshift ISM., Comment: 17 pages, 8 figures, submitted to ApJ

With stunning clarity, JWST has revealed the Universe's first billion years. The scientific community is analyzing a wealth of JWST imaging and spectroscopic data from that era, and is in the process of rewriting the astronomy textbooks. Here, 1.5 years into the JWST science mission, we provide a snapshot of the great progress made towards understanding the initial chapters of our cosmic history. We highlight discoveries and breakthroughs, topics and issues that are not yet understood, and questions that will be addressed in the coming years, as JWST continues its revolutionary observations of the Early Universe. While this compendium is written by a small number of authors, invited to ISSI Bern in March 2024 as part of the 2024 ISSI Breakthrough Workshop, we acknowledge the work of a large community that is advancing our collective understanding of the evolution of the Early Universe., Comment: review article written by the attendees of the 2024 ISSI breakthrough workshop "The first billion year of the Universe", submitted. Comments welcome

[Abridged] In recent years, conflicting results have provided an uncertain view of the dust-attenuated properties of $z>4$ star-forming galaxies (SFGs). To solve this, we used the deepest data publicly available in COSMOS to build a mass-complete ($>10^{9.5}\,M_{\odot}$) sample of SFGs at $4M_\ast)$ converges at $M_\ast<10^{9}\,M_\odot$ and is dominated by SFGs with $M_\ast\sim10^{9.5-10.5}\,M_\odot$. The fraction of the cosmic SFRD that is attenuated by dust, ${\rm SFRD}_{\rm IR}(>M_\ast)/ {\rm SFRD}(>M_\ast)$, is $90\pm4\%$ for $M_\ast\,=\,10^{10}\,M_\odot$, $68\pm10\%$ for $M_\ast=10^{8.9}\,M_\odot$ (i.e., $0.03\times M^\star$; $M^\star$ being the characteristic stellar mass of SFGs) and this value converges to $60\pm10\%$ for $M_\ast=10^{8}\,M_\odot$. Even at this early epoch, the fraction of the cosmic SFRD that is attenuated by dust remains thus significant., Comment: Accepted for publication in A&A; 14 pages and 8 figures

We present high resolution ($\sim0.14$" = 710 pc) ALMA [CII] 158$\mu$m and dust continuum follow-up observations of REBELS-25, a [CII]-luminous ($L_{\mathrm{[CII]}}=(1.7\pm0.2)\times 10^9 \mathrm{L_{\odot}}$) galaxy at redshift $z=7.3065\pm0.0001$. These high resolution, high signal-to-noise observations allow us to study the sub-kpc morphology and kinematics of this massive ($M_* = 8^{+4}_{-2} \times 10^9 \mathrm{M_{\odot}}$) star-forming (SFR$_{\mathrm{UV+IR}} = 199^{+101}_{-63} \mathrm{M_{\odot}} \mathrm{yr}^{-1}$) galaxy in the Epoch of Reionisation. By modelling the kinematics with $^{\mathrm{3D}}$BAROLO, we find it has a low velocity dispersion ($\bar{\sigma} = 33 \pm 9$ km s$^{-1}$) and a high ratio of ordered-to-random motion ($V_{\mathrm{rot, ~max}}/\bar{\sigma} = 11 ^{+8}_{-4}$), indicating that REBELS-25 is a dynamically cold disc. Additionally, we find that the [CII] distribution is well fit by a near-exponential disc model, with a S\'ersic index, $n$, of $1.3 \pm 0.2$, and we see tentative evidence of more complex non-axisymmetric structures suggestive of a bar in the [CII] and dust continuum emission. By comparing to other high spatial resolution cold gas kinematic studies, we find that dynamically cold discs seem to be more common in the high redshift Universe than expected based on prevailing galaxy formation theories, which typically predict more turbulent and dispersion-dominated galaxies in the early Universe as an outcome of merger activity, gas accretion and more intense feedback. This higher degree of rotational support seems instead to be consistent with recent cosmological simulations that have highlighted the contrast between cold and warm ionised gas tracers, particularly for massive galaxies. We therefore show that dynamically settled disc galaxies can form as early as 700 Myr after the Big Bang., Comment: Submitted to MNRAS

The identification of red, apparently massive galaxies at $z>7$ in early JWST photometry suggests a strongly accelerated timeline compared to standard models of galaxy growth. A major uncertainty in the interpretation is whether the red colors are caused by evolved stellar populations, dust, or other effects such as emission lines or AGN. Here we show that three of the massive galaxy candidates at $z=6.7-8.4$ have prominent Balmer breaks in JWST/NIRSpec spectroscopy from the RUBIES program. The Balmer breaks demonstrate unambiguously that stellar emission dominates at $\lambda_{\rm rest} = 0.4\,\mu$m, and require formation histories extending hundreds of Myr into the past in galaxies only 600--800 Myr after the Big Bang. Two of the three galaxies also show broad Balmer lines, with H$\beta$ FWHM $>2500~{\rm km\,s^{-1}}$, suggesting that dust-reddened AGN contribute to, or even dominate, the SEDs of these galaxies at $\lambda_{\rm rest}\gtrsim 0.6\,\mu$m. All three galaxies have relatively narrow [O III] lines, seemingly ruling out a high-mass interpretation if the lines arise in dynamically-relaxed, inclined disks. Yet, the inferred masses also remain highly uncertain. We model the high-quality spectra using Prospector to decompose the continuum into stellar and AGN components, and explore limiting cases in stellar/AGN contribution. This produces a wide range of possible stellar masses, spanning $M_\star \sim 10^9 - 10^{11}\,{\rm M_{\odot}}$. Nevertheless, all fits suggest a very early and rapid formation, most of which follow with a truncation in star formation. Potential origins and evolutionary tracks for these objects are discussed, from the cores of massive galaxies to low-mass galaxies with over-massive black holes. Intriguingly, we find all of these explanations to be incomplete; deeper and redder data are needed to understand the physics of these systems., Comment: Accepted for publication in ApJL; 20 pages, 10 figures, 1 table

We present the first results on the spatial distribution of star formation in 454 star-forming galaxies at $4.8

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 November 2023, MegaScience obtained ~30 arcmin^2 of deep multiband NIRCam imaging centered on the z~0.3 Abell 2744 cluster, including eleven medium-band filters and the two shortest-wavelength broad-band filters, F070W and F090W. Together, MegaScience and the UNCOVER Cycle 1 treasury program provide a complete set of deep (~28-30 mag) images in all NIRCam medium- and broad-band filters. This unique dataset 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 broad-band and four medium-band filters from 0.9-4.8um, expanding the footprint where robust spectral energy distribution (SED) fitting is possible. We provide example SEDs and multi-band 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 [OIII] 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., Comment: 21 pages, 10 figures. Fully reduced imaging, photometric catalogs, and photometric redshift fits publicly available at https://jwst-uncover.github.io/megascience/

We present results from the JWST First Reionization Epoch Spectroscopically Complete Observations survey (FRESCO) on the star forming sequence of galaxies at $1.09.5 M_\odot$ 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$\alpha$)-SFR(Prospector)) is -0.09 $\pm$ 0.04 dex at $10^{10-11} M_\odot$. We also measure the empirical relation between Paschen-$\alpha$ 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, typically extended, Paschen-$\alpha$ emission. Our results may be indicative of the potential unification of methods used to derive the star forming sequence with careful selection of star forming galaxies and independent star formation rate and stellar mass indicators., Comment: 17 pages, 9 figures. Accepted for publication in ApJ

The increasingly neutral intergalactic gas at $z>6$ impacts the Lyman-$\alpha$ flux observed from galaxies. One luminous galaxy, COLA1, stands out because of its unique double-peaked Ly$\alpha$ line at $z=6.6$, unseen in any simulation of reionization. Here we present JWST/NIRCam wide-field slitless spectroscopy in a 21 arcmin$^2$ field centered on COLA1. We find 141 galaxies spectroscopically-selected through the [OIII]($\lambda4969,5008$) doublet at $5.35M_{\rm UV}>-19.5$)., Comment: 15 pages, 13 figures. Accepted for publication in A&A. Main figs. 1, 2, 10 (SFR surface density), 11 (spectra and positions of galaxies in the bubble). Corrected error in escape fraction estimate from SFR surface density

Within the established framework of structure formation, galaxies start as systems of low stellar mass and gradually grow into far more massive galaxies. The existence of massive galaxies in the first billion years of the Universe, suggested by recent observations, appears to challenge this model, as such galaxies would require highly efficient conversion of baryons into stars. An even greater challenge in this epoch is the existence of massive galaxies that have already ceased forming stars. However, robust detections of early massive quiescent galaxies have been challenging due to the coarse wavelength sampling of photometric surveys. Here we report the spectroscopic confirmation with the James Webb Space Telescope of the quiescent galaxy RUBIES-EGS-QG-1 at redshift $z=4.90$, 1.2 billion years after the Big Bang. Deep stellar absorption features in the spectrum reveal that the galaxy's stellar mass of $10^{11}\,M_\odot$, corroborated by the mass implied by its gas kinematics, formed in a short $200\,$Myr burst of star formation, after which star formation activity dropped rapidly and persistently. According to current galaxy formation models, systems with such rapid stellar mass growth and early quenching are too rare to plausibly occur in the small area probed spectroscopically with JWST. Instead, the discovery of RUBIES-EGS-QG-1 implies that early massive quiescent galaxies can be quenched earlier or exhaust gas available for star formation more efficiently than currently assumed., Comment: 9 figures

Combining the public JWST/NIRCam imaging programs CEERS, PRIMER and JADES, spanning a total area of $\sim500\,{\rm arcmin}^2$, we obtain a sample of $>$30,000 galaxies at $z_{\rm phot}\sim4-9$ that allows us to perform a complete, rest-optical selected census of the galaxy population at $z>3$. Comparing the stellar mass $M_*$ and the UV-slope $\beta$ distributions between JWST- and HST-selected samples, we generally find very good agreement and no significant biases. Nevertheless, JWST enables us to probe a new population of UV-red galaxies that was missing from previous HST-based Lyman Break Galaxy (LBG) samples. We measure galaxy stellar mass functions (SMFs) at $z\sim4-9$ down to limiting masses of $10^{7.5}-10^{8.5}\,{\rm M_\odot}$, finding steep low mass slopes over the entire redshift range, reaching values of $\alpha\approx-2$ at $z\gtrsim6$. At the high-mass end, UV-red galaxies dominate at least out to $z\sim6$. The implied redshift evolution of the SMF suggests a rapid build-up of massive dust-obscured or quiescent galaxies from $z\sim6$ to $z\sim4$ as well as an enhanced efficiency of star formation towards earlier times ($z\gtrsim6$). Finally, we show that the galaxy mass density grows by a factor $\sim20\times$ from $z\sim9$ to $z\sim4$. Our results emphasize the importance of rest-frame optically-selected samples in inferring accurate distributions of physical properties and studying the mass build-up of galaxies in the first 1.5 Gyr of cosmic history., Comment: 31 pages, 17 figures, published in MNRAS

The JWST discovery of ``little red dots'' (LRDs) is reshaping our picture of the early Universe, yet the physical mechanisms driving their compact size and UV-optical colors remain elusive. Here we report an unusually bright LRD ($z=3.1$) observed as part of the RUBIES program. This LRD exhibits broad emission lines (FWHM $\sim4000$km/s), a blue UV continuum, a clear Balmer break and a red continuum sampled out to rest 4 $\mu$m with MIRI. We develop a new joint galaxy and AGN model within the Prospector Bayesian inference framework and perform spectrophotometric modeling using NIRCam, MIRI, and NIRSpec/Prism observations. Our fiducial model reveals a $M_*\sim 10^9M_\odot$ galaxy alongside a dust-reddened AGN driving the optical emission. Explaining the rest-frame optical color as a reddened AGN requires $A_{\rm v}\gtrsim4$, suggesting that a great majority of the accretion disk energy is re-radiated as dust emission. Yet despite clear AGN signatures, we find a surprising lack of hot torus emission, which implies that either the dust emission in this object must be cold, or the red continuum must instead be driven by a massive, evolved stellar population of the host galaxy -- seemingly inconsistent with the high EW broad lines (H$\alpha$ EW $\sim800$\AA). The widths and luminosities of Pa$\beta$, Pa$\delta$, Pa$\gamma$, and H$\alpha$ imply a modest black hole mass of $M_{\rm BH}\sim10^8M_\odot$. Additionally, we identify a narrow blue-shifted HeI absorption in G395M spectra, signaling an ionized outflow with kinetic energy up to $\sim1$\% the luminosity of the AGN. The low redshift of RUBIES-BLAGN-1 combined with the depth and richness of the JWST imaging and spectroscopic observations provide a unique opportunity to build a physical model for these so-far mysterious LRDs, which may prove to be a crucial phase in the early formation of massive galaxies and their supermassive black holes., Comment: 20 pages, 8 figures, 4 tables

Lya emitters (LAEs) are particularly useful objects for the study of the Epoch of Reionization. Lya profiles can be used to estimate the amount of ionizing photons that are able to escape the galaxies, and therefore to understand what objects contributed to reionization. However, Lya is a resonant line and its complex radiative transfer effects make the interpretation of the line challenging and require the use of appropriate radiative transfer methods for anything but the simplest gas distributions. With this work we aim to study the properties of simulated LAEs, and the robustness of these inferred properties under a change in the dust model. We also explore the Lyman Continuum (LyC) escape fraction of these galaxies and compare our results with observationally calibrated methods to infer this quantity from the Lya spectrum. We use the radiative transfer code RASCAS to perform synthetic observations of 13 flux-selected galaxies from the Obelisk simulation at redshift z = 6, towards the end of the Epoch of Reionization. Each galaxy was observed in Lya, ionizing and non-ionizing continuum from 48 different viewing angles. We show that the Lya profiles emitted from a galaxy present large variations with a change in viewing angle and that the relation between peak separation and Lya escape fraction is not as strong as previously found, as we find lines of sight with both low peak separation and low escape fraction, due to their dust content. We also show that the properties of the Lya line are reasonably robust under a change in dust model. Lastly, we compare the LyC escape fractions we derive from the simulation to three observationally calibrated methods of inferring this quantity. We determine that none of these relations reproduce the scatter that we find in our sample, and that high escape fraction lines of sight have both low peak separation and low dust extinction in the UV., Comment: Resubmitted to A&A after the first referee report

With its sensitivity in the rest-frame optical, the James Webb Space Telescope (JWST) has uncovered active galactic nuclei (AGN), comprising both intrinsically faint and heavily reddened sources, well into the first billion years of the Universe, at $z \sim 4-11$. In this work, we revisit the AGN contribution to reionization given the high number densities associated with these objects. We use the DELPHI semi-analytic model, base-lined against the latest high-redshift datasets from the JWST and the Atacama Large millimetre Array (ALMA) to model early star forming galaxies and AGN. We calculate the escape fractions of ionizing radiation from both star formation and AGN and include the impact of reionization feeback in suppressing the baryonic content of low-mass galaxies in ionized regions. This model is validated against the key observables for star forming galaxy, AGN and reionization. In our {\it fiducial} model, reionization reaches its mid-point at $z \sim 6.9$ and ends by $z \sim 5.9$. Low stellar mass ($M_*\leq 10^9M_\odot$) star forming galaxies are found to be the key drivers of the reionization process, providing about $77\%$ of the total photon budget. Despite their high numbers, high accretion rates and higher escape fractions compared to star forming galaxies at $z \sim 5$, AGN only provide about $23\%$ of the total reionization budget which is dominated by black holes in high stellar mass systems (with $M_* \geq 10^9M_\odot$). This is because AGN number densities become relevant only at $z \leq 7$ - as a result, AGN contribute as much as galaxies as late as $z \sim 6.2$, when reionization is already in its end stages. Finally, we find that even contrasting models of the AGN ionizing photon escape fraction (increasing or decreasing with stellar mass) do not qualitatively change our results., Comment: Submitted to A&A; comments welcome

GN-z11 is the highest redshift galaxy spectroscopically confirmed with the Hubble Space Telescope (HST). Previous measurements of the effective radius of GN-z11 utilized galfit, which is not optimized to measure structural parameters for such a faint, distant object. Using a new software program called forcepho on HST data for the first time, we derive a size from images in the F160W band obtained both from the complete CANDELS survey and additional midcycle observations in order to contribute to the knowledge base on the size evolution, size-luminosity, and size-mass relation of early galaxies. We find a half-light radius mean of 0''.036 \(\pm\) 0''.006 corresponding to a physical size of 0.15 \(\pm\) 0.025 kpc. This size, smaller than the point spread function, is dramatically smaller than previous estimates with shallower HST data using galfit but consistent with recent measurements using forcepho on new JWST data arXiv:2302.07234. Such a small size, combined with the JWST/NIRSpec spectroscopic observations arXiv:2305.12492, suggests that GN-z11's high luminosity is dominated by an AGN., Comment: 4 pages, 1 figure

With the remarkable sensitivity and resolution of JWST in the infrared, measuring rest-optical kinematics of galaxies at $z>5$ has become possible for the first time. This study pilots a new method for measuring galaxy dynamics for highly multiplexed, unbiased samples by combining FRESCO NIRCam grism spectroscopy and JADES medium-band imaging. Here we present one of the first JWST kinematic measurements for a galaxy at $z>5$. We find a significant velocity gradient, which, if interpreted as rotation yields $V_{rot} = 240\pm50$km/s and we hence refer to this galaxy as Twister-z5. With a rest-frame optical effective radius of $r_e=2.25$kpc, the high rotation velocity in this galaxy is not due to a compact size as may be expected in the early universe but rather a high total mass, ${\rm log(M}_{dyn}/{\rm M}_\odot)=11.0\pm0.2$. This is a factor of roughly 4x higher than the stellar mass within the effective radius. We also observe that the radial H$\alpha$ equivalent width profile and the specific star formation rate map from resolved stellar population modeling is centrally depressed by a factor of $\sim1.5$ from the center to $r_e$. Combined with the morphology of the line-emitting gas in comparison to the continuum, this centrally suppressed star formation is consistent with a star-forming disk surrounding a bulge growing inside-out. While large, rapidly rotating disks are common to z~2, the existence of one after only 1Gyr of cosmic time, shown for the first time in ionized gas, adds to the growing evidence that some galaxies matured earlier than expected in the history of the universe., Comment: Fig. 3 shows the main result

The recent UNCOVER survey with the James Webb Space Telescope (JWST) exploits the nearby cluster Abell 2744 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 \sim 15$. We show the redshift evolution of stellar ages, star formation rates, and rest-frame colors across the full range of $0.2 \lesssim z \lesssim 15$. The galaxy properties are inferred using the Prospector Bayesian inference framework using informative Prospector-$\beta$ priors on masses and star formation histories to produce joint redshift and stellar population posteriors, and 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 $\sigma_{\rm NMAD} \sim 0.03$, of a similar quality to the established photometric redshift code EAzY. In line with the open-source scientific objective of the Treasury survey, we publicly release the stellar population catalog with this paper, derived from the photometric catalog adapting aperture sizes based on source profiles. This release 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., Comment: Corrected typos: Eq.1 should've been (1-kappa)^2, and the lens maps are normalized to D_ds/D_s=1. These errors were only in the writing; no data products or results were affected. The SPS catalogs are accessible via the UNCOVER survey webpage: https://jwst-uncover.github.io/DR2.html#SPSCatalogs, with a copy deposited to Zenodo: https://doi.org/10.5281/zenodo.8401181

We present the survey design and initial results of the ALMA Cycle 9 program of DUALZ, which aims to establish a joint ALMA and JWST public legacy field targeting the massive galaxy cluster Abell 2744. DUALZ features a contiguous $4'\times6'$ ALMA 30-GHz-wide mosaic in Band 6, covering areas of $\mu>2$ down to a sensitivity of $\sigma=32.7~\mu$Jy. Through a blind search, we identified 69 dust continuum sources at S/N $\gtrsim5.0$ with median redshift and intrinsic 1.2-mm flux of $z=2.30$ and $S_{\rm 1.2mm}^{\rm int}=0.24$~mJy. Of these, 27 have been spectroscopically confirmed, leveraged by the latest NIRSpec observations, while photometric redshift estimates are constrained by the comprehensive HST, NIRCam, and ALMA data for the remaining sources. With priors, we further identify a [CII]158 $\mu$m line emitter at $z=6.3254\pm0.0004$, confirmed by the latest NIRSpec spectroscopy. The NIRCam counterparts of the 1.2-mm continuum exhibit undisturbed morphologies, denoted either by disk or spheroid, implying the triggers for the faint mm emission are less catastrophic than mergers. We have identified 8 HST-dark galaxies (F150W$>$27mag, F150W$-$F444W$>$2.3) and 2 JWST-dark (F444W$>$30mag) galaxy candidates among the ALMA continuum sources. The former includes face-on disk galaxies, hinting that substantial dust obscuration does not always result from inclination. We also detect a marginal dust emission from an X-ray-detected galaxy at $z_{\rm spec}=10.07$, suggesting an active co-evolution of the central black hole and its host. We assess the infrared luminosity function up to $z\sim10$ and find it consistent with predictions from galaxy formation models. To foster diverse scientific outcomes from the community, we publicly release reduced ALMA mosaic maps, cubes, and the source catalog., Comment: 33 pages, 16 figures, and 5 tables. Submitted to ApJS. The ALMA products are fully available from here: https://jwst-uncover.github.io/DR2.html#DUALZ

JWST is revealing a new population of dust-reddened broad-line active galactic nuclei (AGN) at redshifts $z\gtrsim5$. Here we present deep NIRSpec/Prism spectroscopy from the Cycle 1 Treasury program 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 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$\alpha$ with FWHM$\, >2000$ km/s, for $20\%$ 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_{\rm 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$_{\odot}$. While their UV-luminosities ($-16>M_{\rm 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$_{\odot}$ black holes radiating at $\sim 10-40\%$ of Eddington. The number densities are surprisingly high at $\sim10^{-5}$ Mpc$^{-3}$ mag$^{-1}$, 100 times more common than the faintest UV-selected quasars, while accounting for $\sim1\%$ of the UV-selected galaxies. While their UV-faintness suggest they may not contribute strongly to reionization, their ubiquity poses challenges to models of black hole growth., Comment: 23 pages, 9 figures, 5 tables, submitted to ApJ

We report the detection of a 13$\sigma$ H$\alpha$ emission line from HDF850.1 at $z=5.188\pm0.001$ using the FRESCO NIRCam F444W grism observations. Detection of H$\alpha$ in HDF850.1 is noteworthy, given its high far-IR luminosity, substantial dust obscuration, and the historical challenges in deriving its redshift. HDF850.1 shows a clear detection in the F444W imaging data, distributed between a northern and southern component, mirroring that seen in [CII] from the Plateau de Bure Interferometer. Modeling the SED of each component separately, we find that the northern component has a higher mass, star formation rate (SFR), and dust extinction than the southern component. The observed H$\alpha$ emission appears to arise entirely from the less-obscured southern component and shows a similar $\Delta$v$\sim$+130 km/s velocity offset to that seen for [CII] relative to the source systemic redshift. Leveraging H$\alpha$-derived redshifts from FRESCO observations, we find that HDF850.1 is forming in one of the richest environments identified to date at $z>5$, with 100 $z=5.17-5.20$ galaxies distributed across 10 structures and a $\sim$(15 cMpc)$^3$ volume. Based on the evolution of analogous structures in cosmological simulations, the $z=5.17-5.20$ structures seem likely to collapse into a single $>$10$^{14}$ $M_{\odot}$ cluster by $z\sim0$. Comparing galaxy properties forming within this overdensity with those outside, we find the masses, SFRs, and $UV$ luminosities inside the overdensity to be clearly higher. The prominence of H$\alpha$ line emission from HDF850.1 and other known highly-obscured $z>5$ galaxies illustrates the potential of NIRCam-grism programs to map both the early build-up of IR-luminous galaxies and overdense structures., Comment: Submitted to MNRAS. 20 pages, 10 figures and 8 tables (including appendices)

Recent JWST observations have revealed an unexpected abundance of massive galaxy candidates in the early Universe, extending further in redshift and to lower luminosity than what had previously been found by sub-millimeter surveys. These JWST candidates have been interpreted as challenging the $\Lambda$CDM cosmology, but, so far, they have mostly relied only on rest-frame ultraviolet data and lacked spectroscopic confirmation of their redshifts. Here we report a systematic study of 36 massive dust-obscured galaxies with spectroscopic redshifts between $z_{\rm spec}=5-9$ from the JWST FRESCO survey. We find no tension with the $\Lambda$CDM model in our sample. However, three ultra-massive galaxies (log$M_{\star}/M_{\odot}$ $\gtrsim11.0$) require an exceptional fraction of 50% of baryons converted into stars -- two to three times higher than even the most efficient galaxies at later epochs. The contribution from an active nucleus is unlikely because of their extended emission. Ultra-massive galaxies account for as much as 17% of the total cosmic star formation rate density at $z\sim5-6$., Comment: Nature in press. Updated to the accepted version. 24 pages, 4 main figures, 7 supplementary figures, 3 supplementary tables

Deep observations with JWST have revealed an emerging population of red point-like 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$\beta$ line with FWHM = $3439\pm413$ 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 the black hole mass of log$_{10}(M_{\rm BH}/M_\odot)=8.17\pm0.42$, and a bolometric luminosity of $L_{\rm bol}\sim6.6\times10^{45}$ erg s$^{-1}$. These values imply that our object is accreting at $\sim 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 log$_{10}(M_{\rm *}/M_\odot)<8.7$, which would lead to an unprecedented ratio of black hole to host mass of at least $\sim 30 \%$. This is orders of magnitude higher compared to the local QSOs, but is consistent with recent AGN studies at high redshift with JWST. This finding suggests that a non-negligible 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 the galaxy-black hole co-evolution in the early Universe., Comment: 14 pages, 6 figures, 2 tables. Accepted for a publication in ApJL

We present JWST NIRSpec prism spectroscopy of gravitationally lensed galaxies at $z\gtrsim9$ found behind the massive galaxy cluster Abell 2744 in the UNCOVER Cycle 1 Treasury Program. We confirm the source redshift via emission lines and/or the Ly$\alpha$ break feature for ten galaxies at z=8.50-13.08 down to $M_{\rm UV}=-17.3$. We achieve a high confirmation rate of 100\% for $z>9$ candidates reported in Atek et al. (2023). Using six sources with multiple emission line detections, we find that the offset of the redshift estimates between the lines and the Ly$\alpha$ break alone with prism can be as large as $\pm0.2$, raising caution in designing future follow-up spectroscopy for the break-only sources. With spec-$z$ confirmed sources in UNCOVER and the literature, we derive lower limits on the rest-frame ultraviolet (UV) luminosity function (LF) at $z\simeq9$-12 and find these lower limits to be consistent with recent photometric measurements. We identify at least two unambiguous and several possible active galactic nucleus (AGN) systems based on X-ray emission, broad line (BL) H$\beta$, high ionization line (e.g., NIV]1487, CIV1549) detections, and excess in UVLF. This requires the AGN LFs at $z\simeq$ 9-10 to be comparable or even higher than the X-ray AGN LF estimated at $z\sim6$ and indicates a plausible cause of the high abundance of $z>9$ galaxies claimed in recent photometric studies may be AGNs. One UV-luminous source is confirmed at the same redshift as a dusty BL AGN at $z=8.50$ with a physical separation of 380 kpc in the source plane. These two sources show blueward Ly$\alpha$ line or continuum emission, suggesting that they reside in the same ionized bubble with a radius of $7.69\pm0.18$ pMpc. Our results imply that AGNs have a non-negligible contribution to cosmic reionization., Comment: 27 pages, 11 figures, 4 tables, submitted to ApJ (See also arXiv:2308.11610)

The identification of sources driving cosmic reionization, a major phase transition from neutral Hydrogen to ionized plasma around 600-800 Myr after the Big Bang (Dayal et al. 2018, Mason et al. 2019, Robertson et al. 2022), has been a matter of intense debate (Robertson et al. 2022). Some models suggest that high ionizing emissivity and escape fractions ($f_{\rm esc}$) from quasars support their role in driving cosmic reionization (Madau & Haardt 2015, Mitra et al. 2018). Others propose that the high $f_{\rm esc}$ values from bright galaxies generates sufficient ionizing radiation to drive this process (Naidu et al. 2020). Finally, a few studies suggest that the number density of faint galaxies, when combined with a stellar-mass-dependent model of ionizing efficiency and $f_{\rm esc}$, can effectively dominate cosmic reionization (Finkelstein et al. 2019, Dayal et al. 2020). However, so far, low-mass galaxies have eluded comprehensive spectroscopic studies owing to their extreme faintness. Here we report an analysis of eight ultra-faint galaxies (in a very small field) during the epoch of reionization with absolute magnitudes between $M_{\rm UV}$ $\sim -17$ to $-15$ mag (down to 0.005 $L^{\star}$. We find that faint galaxies during the Universe's first billion years produce ionizing photons with log($\xi_{\mathrm{ion}}$/ Hz erg$^{-1}$) =$25.80\pm 0.14$, a factor of 4 higher than commonly assumed values (Robertson et al. 2015). If this field is representative of the large scale distribution of faint galaxies, the rate of ionizing photons exceeds that needed for reionization, even for escape fractions of order five per cent., Comment: 29 pages, 7 figures, 2 tables. Published in Nature

Early JWST observations have uncovered a new population of red sources that might represent a previously overlooked phase of supermassive black hole growth (Kocevski et al. 2023; Matthee et al. 2023, Labb\'e et al. 2023). One of the most intriguing examples is an extremely red, point-like object that was found to be triply-imaged by the strong lensing (SL) cluster Abell 2744 (Furtak et al. 2023). Here we present deep JWST/NIRSpec observations of this object, Abell2744-QSO1. The spectroscopy confirms that the three images are of the same object, and that it is a highly reddened ($A_V\simeq3$) broad emission-line Active Galactic Nucleus (AGN) at a redshift of $z_{\mathrm{spec}}=7.0451\pm0.0005$. From the width of H$\beta$ ($\mathrm{FWHM}=2800\pm250\,\frac{\mathrm{km}}{\mathrm{s}}$) we derive a black hole mass of $M_{\mathrm{BH}}=4_{-1}^{+2}\times10^7\,\mathrm{M}_{\odot}$. We infer a very high ratio of black hole to galaxy mass of at least 3%, an order of magnitude more than is seen in local galaxies (Bennert et al. 2011), and possibly as high as 100%. The lack of strong metal lines in the spectrum together with the high bolometric luminosity ($L_{\mathrm{bol}}=(1.1\pm0.3)\times10^{45}\,\frac{\mathrm{erg}}{\mathrm{s}}$) indicate that we are seeing the black hole in a phase of rapid growth, accreting at 30% of the Eddington limit. The rapid growth and high black hole to galaxy mass ratio of A2744-QSO1 suggest that it may represent the missing link between black hole seeds (Volonteri et al. 2021) and the first luminous quasars (Fan et al. 2022)., Comment: Published in Nature. Updated to the accepted version

The formation of galaxies by gradual hierarchical co-assembly of baryons and cold dark matter halos is a fundamental paradigm underpinning modern astrophysics and predicts a strong decline in the number of massive galaxies at early cosmic times. Extremely massive quiescent galaxies (stellar masses $>10^{11}$ M$_\odot$) have now been observed as early as 1-2 billions years after the Big Bang; these are extremely constraining on theoretical models as they form 300-500 Myr earlier and only some models can form massive galaxies this early. Here we report on the spectroscopic observations with the James Webb Space Telescope of a massive quiescent galaxy ZF-UDS-7329 at redshift 3.205 $\pm$ 0.005 that eluded deep ground-based spectrscopy, is significantly redder than typical and whose spectrum reveals features typical of much older stellar populations. Detailed modelling shows the stellar population formed around 1.5 billion years earlier in time (z ~ 11) at an epoch when dark matter halos of sufficient hosting mass have not yet assembled in the standard scenario. This observation may point to the presence of undetected populations of early galaxies and the possibility of significant gaps in our understanding of early stellar populations, galaxy formation and/or the nature of dark matter., Comment: 21 pages, 6 figures, updated to reflect accepted version. v3 with figure caption clarified

Observations of high-redshift galaxies provide a critical direct test to the theories of early galaxy formation, yet to date, only three have been spectroscopically confirmed at $z>12$. Due to strong gravitational lensing over a wide area, the galaxy cluster field A2744 is ideal for searching for the earliest galaxies. Here we present JWST/NIRSpec observations of two galaxies: a robust detection at $z_{\rm spec} = 12.393^{+0.004}_{-0.001}$, and a plausible candidate at $z_{\rm spec} = 13.079^{+0.013}_{-0.001}$. The galaxies are discovered in JWST/NIRCam imaging and their distances are inferred with JWST/NIRSpec spectroscopy, all from the JWST Cycle 1 UNCOVER Treasury survey. Detailed stellar population modeling using JWST NIRCam and NIRSpec data corroborates the primeval characteristics of these galaxies: low mass ($\sim 10^8~{\rm M_\odot}$), young, rapidly-assembling, metal-poor, and star-forming. Interestingly, both galaxies are spatially resolved, having lensing-corrected rest-UV effective radii on the order of 300-400 pc, which are notably larger than other spectroscopically confirmed systems at similar redshifts. The observed dynamic range of $z \gtrsim 10$ sizes spans over 1 order of magnitude, implying a significant scatter in the size-mass relation at early times. Deep into the epoch of reionization, these discoveries elucidate the emergence of the first galaxies., Comment: ApJL in press; 16 pages, 6 figures, 2 tables

The James Webb Space Telescope is now detecting early black holes (BHs) as they transition from "seeds" to supermassive BHs. Recently Bogdan et al. (2023) reported the detection of an X-ray luminous supermassive BH, UHZ-1, with a photometric redshift at $z > 10$. Such an extreme source at this very high redshift provides new insights on seeding and growth models for BHs given the short time available for formation and growth. Harnessing the exquisite sensitivity of JWST/NIRSpec, here we report the spectroscopic confirmation of UHZ-1 at $z = 10.073 \pm 0.002$. We find that the NIRSpec/Prism spectrum is typical of recently discovered z~10 galaxies, characterized primarily by star-formation features. We see no clear evidence of the powerful X-ray source in the rest-frame UV/optical spectrum, which may suggest heavy obscuration of the central BH, in line with the Compton-thick column density measured in the X-rays. We perform a stellar population fit simultaneously to the new NIRSpec spectroscopy and previously available photometry. The fit yields a stellar mass estimate for the host galaxy that is significantly better constrained than prior photometric estimates ($M_*\sim 1.4^{+0.3}_{-0.4} \times 10^8 M_\odot$). Given the predicted BH mass ($M_{\rm BH}\sim10^7-10^8 M_\odot$), the resulting ratio of $M_{\rm BH}/M_*$ remains two to three orders of magnitude higher than local values, thus lending support to the heavy seeding channel for the formation of supermassive BHs within the first billion years of cosmic evolution., Comment: 10 pages, 4 figures, ApJL in-press. Updated to accepted version incl. X-ray image

The James Webb Space Telescope (JWST) is revolutionizing our knowledge of $z>5$ galaxies and their actively accreting black holes. Using the JWST Cycle 1 Treasury program Ultradeep NIRSpec and NIRCam ObserVations before the Epoch of Reionization (UNCOVER) in the lensing field Abell 2744, we report the identification of a sample of little red dots at $3 < z_{\rm{phot}} < 7$ that likely contain highly-reddened accreting supermassive black holes. Using a NIRCam-only selection to F444W$<27.7$ mag, we find 26 sources over the $\sim45$ arcmin$^{2}$ field that are blue in F115W$-$F200W$\sim0$ (or $\beta_{\rm UV}\sim-2.0$ for $f_{\lambda} \propto \lambda^\beta$), red in F200W$-$F444W = $1-4$ ($\beta_{\rm opt} \sim +2.0$), and are dominated by a point-source like central component. Of the 20 sources with deep ALMA 1.2-mm coverage, none are detected individually or in a stack. For the majority of the sample, SED fits to the JWST+ALMA observations prefer models with hot dust rather than obscured star-formation to reproduce the red NIRCam colors and ALMA 1.2-mm non-detections. While compact dusty star formation can not be ruled out, the combination of extremely small sizes ($\langle r_e \rangle\approx50$ pc after correction for magnification), red rest-frame optical slopes, and hot dust can by explained by reddened broad-line active galactic nuclei (AGNs). Our targets have faint $M_{\rm 1450} \approx -14\ \, {\rm to} -18$ mag but inferred bolometric luminosities of $L_{\rm bol} = 10^{43}-10^{46}$ erg/s, reflecting their obscured nature. If the candidates are confirmed as AGNs with upcoming UNCOVER spectroscopy, then we have found an abundant population of reddened luminous AGN that are at least ten times more numerous than UV-luminous AGN at the same intrinsic bolometric luminosity., Comment: submitted to ApJ

Characterising the prevalence and properties of faint active galactic nuclei (AGN) 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$\alpha$ emitters at $z\approx4-6$ using deep JWST/NIRCam imaging and wide field slitless spectroscopy from the EIGER and FRESCO surveys. We identify 20 H$\alpha$ lines at $z=4.2-5.5$ that have broad components with line widths from $\sim1200-3700$ km s$^{-1}$, contributing $\sim30-90$ % of the total line flux. We interpret these broad components as being powered by accretion onto SMBHs with implied masses $\sim10^{7-8}$ M$_{\odot}$. In the UV luminosity range M$_{\rm UV}=-21$ to $-18$, we measure number densities of $\approx10^{-5}$ cMpc$^{-3}$. This is an order of magnitude higher than expected from extrapolating quasar UV luminosity functions. Yet, such AGN are found in only $<1$ % of star-forming galaxies at $z\sim5$. The SMBH mass function agrees with large cosmological simulations. In two objects we detect narrow red- and blue-shifted H$\alpha$ absorption indicative, respectively, of dense gas fueling SMBH growth and outflows. We may be witnessing early AGN feedback that will clear dust-free pathways through which more massive blue quasars are seen. We uncover a strong correlation between reddening and the fraction of total galaxy luminosity arising from faint AGN. This implies that early SMBH growth is highly obscured and that faint AGN are only minor contributors to cosmic reionization., Comment: 26 pages, 20 figures. Accepted for publication in ApJ. Various improvements and updated citations. Main additions: colours (Fig 11), broad Halpha LF (Fig 18, Sect 5.2), comparisons to simulations (Sect 5.3)

The onset of galaxy formation is thought to be initiated by the infall of neutral, pristine gas onto the first protogalactic halos. However, direct constraints on the abundance of neutral atomic hydrogen (HI) in galaxies have been difficult to obtain at early cosmic times. Here we present spectroscopic observations with JWST of three galaxies at redshifts $z=8.8 - 11.4$, about $400-600$ Myr after the Big Bang, that show strong damped Lyman-$\alpha$ absorption ($N_{\rm HI} > 10^{22}$ cm$^{-2}$) from HI in their local surroundings, an order of magnitude in excess of the Lyman-$\alpha$ absorption caused by the neutral intergalactic medium at these redshifts. Consequently, these early galaxies cannot be contributing significantly to reionization, at least at their current evolutionary stages. Simulations of galaxy formation show that such massive gas reservoirs surrounding young galaxies so early in the history of the universe is a signature of galaxy formation in progress., Comment: Submitted

We present the results of a search for high-redshift ($z>9$) galaxy candidates in the JWST UNCOVER survey, using deep NIRCam and NIRISS imaging in 7 bands over $\sim45$ arcmin$^2$ and ancillary HST observations. The NIRCam observations reach a $5-\sigma$ limiting magnitude of $\sim 29.2$ AB. The identification of high$-z$ candidates relies on a combination of a dropout selection and photometric redshifts. We find 16 candidates at $99$. A few galaxies at $z\sim9-10$ appear to show a clear Balmer break between the F356W and F444W/F410M bands, which helps constrain their stellar mass. We estimate blue UV continuum slopes between $\beta=-1.8$ and $-2.3$, typical for early galaxies at $z>9$ but not as extreme as the bluest recently discovered sources. We also find evidence for a rapid redshift-evolution of the mass-luminosity relation and a redshift-evolution of the UV continuum slope for a given range of intrinsic magnitude, in line with theoretical predictions. These findings suggest that deeper JWST observations are needed to reach the fainter galaxy population at those early epochs, and follow-up spectroscopy will help better constrain the physical properties and star formation histories of a larger sample of galaxies., Comment: Submitted to MNRAS