Search

We present a spectroscopic analysis of two star-forming galaxies at z~5 observed with JWST/NIRSpec as part of the Early eXtragalactic Continuum and Emission Line Science (EXCELS) survey. The detection of the C III]$\lambda\lambda$1906,09, [O II]$\lambda\lambda$3726,29, [O III]$\lambda\lambda$4363,5007, and [N II]$\lambda$6584 nebular emission lines enables investigation of the C/O, N/O, and C/N abundance ratios using the temperature-sensitive method. The two galaxies have stellar masses of log($M_{\star}$/M$_{\odot}$ ) = 8.13$\pm$0.09 and log($M_{\star}$/M$_{\odot}$ )=8.52$\pm$0.13 and corresponding metallicities of Z~0.2Z$_{\odot}$ and Z~0.3Z$_{\odot}$. These metallicities are somewhat higher than is typical for other z>5 galaxies with similar stellar mass and are in fact comparable to high-redshift analogue galaxies at z~0. Both galaxies display evidence for N/O enhancement with respect to the z~0 sample, with log(N/O)=-1.07$\pm$0.17 and log(N/O)=-0.86$\pm$0.15 respectively. In contrast, we find low C abundances, with log(C/O)=-0.82$\pm$0.22 and log(C/O)=-1.02$\pm$0.22, consistent with the predicted yields of core-collapse supernovae. Following the trend observed in other high-redshift sources, we find that the C/N ratios are lower at fixed O/H compared to the majority of local galaxies. In contrast to the top-heavy IMF invoked in some studies to explain low C/N ratios in metal-poor galaxies, we find, via comparison to chemical evolution models, that a standard or bottom-heavy IMF better explains the observed abundance ratios in more enriched systems due to an increase in N-enrichment from intermediate mass (4-7M$_{\odot}$) stars. Our results demonstrate that robust measurements of CNO abundances with JWST can reveal unique enrichment pathways in galaxies as a function of both metallicity and redshift., Comment: 17 pages, 8 Figures, and 3 Tables. Submitted to MNRAS. Comments Welcome

In this paper we study a combinatorial reconfiguration problem that involves finding an optimal sequence of swaps to move an initial configuration of tokens that are placed on the vertices of a graph to a final desired one. This problem arises as a crucial step in reducing the depth of a quantum circuit when compiling a quantum algorithm. We provide the first known constant factor approximation algorithms for the parallel token swapping problem on graph topologies that are commonly found in modern quantum computers, including cycle graphs, subdivided star graphs, and grid graphs. We also study the so-called stretch factor of a natural lower bound to the problem, which has been shown to be useful when designing heuristics for the qubit routing problem. Finally, we study the colored version of this reconfiguration problem where some tokens share the same color and are considered indistinguishable.

We present an analysis of nine star-forming galaxies with $\langle z \rangle = 3.95$ from the JWST EXCELS survey for which we obtain robust chemical abundance estimates for the $\alpha$-elements O, Ne and Ar. The $\alpha$-elements are primarily produced via core-collapse supernovae (CCSNe) which should result in $\alpha$-element abundance ratios that do not vary significantly across cosmic time. However, Type Ia supernovae (SNe Ia) models predict an excess production of Ar relative to O and Ne. The Ar/O abundance ratio can therefore be used as a tracer of the relative enrichment of CCSNe and SNe Ia in galaxies. Our sample approximately doubles the number of sources with measurements of ${\rm Ar/O}$ at $z > 2$, and we find that our sample exhibits sub-solar Ar/O ratios on average, with $\rm{Ar/O} = 0.62 \pm 0.10 \, (\rm{Ar/O})_{\odot}$. In contrast, the average Ne/O abundance is fully consistent with the solar ratio, with $\rm{Ne/O} = 1.07 \pm 0.12 \, (\rm{Ne/O})_{\odot}$. Our results support a scenario in which Ar has not had time to build up in the interstellar medium of young high-redshift galaxies, which are dominated by CCSNe enrichment. We show that these abundance estimates are in good agreement with recent Milky Way chemical evolution models, and with Ar/O trends observed for planetary nebulae in the Andromeda galaxy. These results highlight the potential for using multiple element abundance ratios to constrain the chemical enrichment pathways of early galaxies with JWST., Comment: 14 pages, 9 figures. Submitted to MNRAS

We present the Early Data Release of the Multi-Cycle JWST-NEXUS Treasury program (2024-2028), which includes NIRCam imaging and WFSS observations from the first (partial) NEXUS-Wide epoch covering the central 100 ${\rm arcmin^2}$ of the NEXUS field, located near the North Ecliptic Pole and within the Euclid Ultra-Deep Field. We release reduced NIRCam mosaics (F090W, F115W, F150W, F200W, F356W, F444W), photometric source catalogs, as well as preliminary WFSS spectra (in F322W2 and F444W) for the subset of bright sources (F356W$<$21 mag or F444W$<$21 mag). These observations fully cover the NEXUS-Deep area, and anchor the long-term baseline of the program. These data will be used for initial target selection for the NIRSpec/MSA spectroscopy starting from June 2025. The NIRCam imaging reaches depths of 27.4--28.2 (AB) mags in F090W--F444W. Upcoming NEXUS-Wide epochs will expand the area to the full $\sim 400\,{\rm arcmin^2}$, and improve the NIRCam exposure depths in the Wide tier by a factor of three. In addition, this central region will be repeatedly covered by the NEXUS-Deep observations (NIRCam imaging and NIRSpec/MSA PRISM spectroscopy) over 18 epochs with a $\sim 2$-month cadence. We demonstrate the data quality of the first NEXUS observations, and showcase some example science cases enabled by these data., Comment: 14 pages, 9 figures, 2 tables. Data products are publicly accessible at https://ariel.astro.illinois.edu/nexus/edr/ Online interactive map for quick visualization of released images and WFSS spectra can be found at https://ariel.astro.illinois.edu/nexus/map/

Accurately quantifying the ionizing photon production efficiency ($\xi_\rm{ion}$) of $z>6$ star-forming galaxies (SFGs) is necessary to understand their contribution to reionization. We investigate the ionizing properties of N=289 SFGs selected at $z=6.9-7.6$ from the JWST Cycle-1 imaging programmes; PRIMER and JADES. We use BAGPIPES to consistently infer the equivalent widths of their [OIII]+$\rm{H\beta}$ emission lines ($W_\lambda$) and their physical properties. To supplement our high-redshift galaxies, we measure $W_\lambda$ photometrically for a sample of N=253 $z=3.2-3.6$ SFGs selected from the VANDELS spectroscopic survey. Comparing these samples, we find a strong apparent redshift evolution in their median $W_\lambda$, increasing from $W_\lambda =310\pm25\,\r{A}$ in VANDELS to $W_\lambda =630\pm30\,\r{A}$ in our JWST-based sample. In the JWST sample at $z>7$, we find that $W_\lambda$ correlates with both stellar mass and UV luminosity, with high-mass, $M_{ UV}-$faint galaxies producing systematically weaker emission lines. Moreover, we discover a departure from the standard log-normal shape of the $W_\lambda$ distribution, characterised by a more pronounced tail at lower $W_\lambda$, consistent with increasingly bursty star formation. Using $W_\lambda$ as a proxy for $\xi_\rm{ion}$, and UV spectral slope as a proxy for LyC escape fraction ($f_\rm{esc}$), we find a minority of galaxies with high $\xi_\rm{ion}$ and $f_\rm{esc}$ (e.g., $\rm{log(\xi_{ion}/erg^{-1}Hz})\sim25.6$ and $f_\rm{esc}\sim0.15$). However, we find that the LyC photon budget at $z>7$ is dominated by galaxies with more moderate output, close to the sample median of $\rm{log(\xi_{ion}/erg^{-1}Hz})\sim25.3$ and $f_\rm{esc}\sim0.05$. This is consistent with estimates for the number of LyC photons required to power reionization at $z>7$, with no evidence for over or under-production., Comment: 18 pages, 9 figures + 1 appendix. Submitted to MNRAS

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 examine star-formation and dust properties for a sample of 660 galaxies at $1.37\leq z\leq 2.61$ in the MOSDEF survey by dividing them into groups with similarly-shaped spectral energy distributions (SEDs). For each group, we combine the galaxy photometry into a finely-sampled composite SED, and stack their spectra. This method enables the study of more complete galaxy samples, including galaxies with very faint emission lines. We fit these composite SEDs with Prospector to measure the stellar attenuation and SED-based star-formation rates (SFRs). We also derive emission-line properties from the spectral stacks, including Balmer decrements, dust-corrected SFRs, and metallicities. We find that stellar attenuation correlates most strongly with mass, while nebular attenuation correlates strongly with both mass and SFR. Furthermore, the excess of nebular compared to stellar attenuation correlates most strongly with SFR. The highest SFR group has 2 mag of excess nebular attenuation. Our results are consistent with a model in which star-forming regions become more dusty as galaxy mass increases. To explain the increasing excess nebular attenuation, we require a progressively larger fraction of star formation to occur in highly-obscured regions with increasing SFR. This highly-obscured star formation could occur in dusty clumps or central starbursts. Additionally, as each galaxy group represents a different evolutionary stage, we study their locations on the UVJ and SFR-mass diagrams. As mass increases, metallicity and dust attenuation increase, while sSFR decreases. However, the most massive group moves towards the quiescent region of the UVJ diagram, while showing less obscuration, potentially indicating removal of dust., Comment: 21 pages, 7 figures, accepted for publication in ApJ

NEXUS is a JWST Multi-Cycle (Cycles 3-5; 368 primary hrs) GO Treasury imaging and spectroscopic survey around the North Ecliptic Pole. It contains two overlapping tiers. The Wide tier ($\sim 400~{\rm arcmin}^2$) performs NIRCam/WFSS 2.4-5 micron grism spectroscopy with three epochs over 3 years (final continuum ${\rm S/N/pixel>3}$ at F444W$<22.2$). The Deep tier ($\sim 50~{\rm arcmin}^2$) performs high-multiplexing NIRSpec 0.6-5.3 micron MOS/PRISM spectroscopy for $\sim 10,000$ targets, over 18 epochs with a 2-month cadence (epoch/final continuum ${\rm S/N/pixel>3}$ at F200W$\lesssim 27/29$). All epochs have simultaneous multi-band NIRCam and MIRI imaging ($5\sigma$ final depths of $\sim 28-29$ in NIRCam and $\sim 25$ in MIRI). The field is within the continuous viewing zone of JWST, and is fully covered by the Euclid Ultra-Deep Field, with 0.9-2 micron deep Euclid spectroscopy and cadenced photometry. NEXUS has three science pillars. First, with its massive and nearly complete (flux-limited) spectroscopic samples and deep photometry, it will perform efficient classification and physical characterization of galaxies and AGNs from $z\sim 1$ to Cosmic Dawn. With the large contiguous area coverage, it will measure the spatial clustering and demography of the first galaxies and SMBHs at $z>6$. Second, multi-epoch observations enable systematic time-domain investigations, focusing on $z>3$ transients and low-mass AGN reverberation mapping. Third, the comprehensive data set will enable knowledge transfer to other legacy fields, create data challenges, and initiate benchmark work for future space missions. With rapid public releases of processed data and an open invitation for collaboration, NEXUS aims for broad and swift community engagement, to become a powerhouse to drive transformative advancements in multiple key science areas of astronomy., Comment: 19 pages

Integral field spectroscopy (IFS) is a powerful tool for understanding the formation of galaxies across cosmic history. We present the observing strategy and first results of MSA-3D, a novel JWST program using multi-object spectroscopy in a slit-stepping strategy to produce IFS data cubes. The program observed 43 normal star-forming galaxies at redshifts $0.5 \lesssim z \lesssim 1.5$, corresponding to the epoch when spiral thin-disk galaxies of the modern Hubble sequence are thought to emerge, obtaining kpc-scale maps of rest-frame optical nebular emission lines with spectral resolution $R\simeq2700$. Here we describe the multiplexed slit-stepping method which is $>15$ times more efficient than the NIRSpec IFS mode for our program. As an example of the data quality, we present a case study of an individual galaxy at $z=1.104$ (stellar mass $M_{*} = 10^{10.3}~M_{\odot}$, star formation rate~$=3~M_{\odot}$ yr$^{-1}$) with prominent face-on spiral structure. We show that the galaxy exhibits a rotationally supported disk with moderate velocity dispersion ($\sigma = 36^{+5}_{-4}$~\kms), a negative radial metallicity gradient ($-0.020\pm0.002$~dex\,kpc$^{-1}$), a dust attenuation gradient, and an exponential star formation rate density profile which closely matches the stellar continuum. These properties are characteristic of local spirals, indicating that mature galaxies are in place at $z\sim1$. We also describe the customized data reduction and original cube-building software pipelines which we have developed to exploit the powerful slit-stepping technique. Our results demonstrate the ability of JWST slit-stepping to study galaxy populations at intermediate to high redshifts, with data quality similar to current surveys of the $z\sim0.1$ universe., Comment: 24 pages, 12 figures, 1 table; submitted to ApJ, comments welcome. Our custom designed data reduction pipeline (MSA3D) and example processed data-cube will be publicly released before the JWST Cycle 4 deadline. MSA3D package on GitHub: https://github.com/barisiciv/msa3d.git

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 have used a combined sample of JADES and CEERS objects in order to constrain ionizing photon production efficiency ($\xi_{\rm ion}$) from JWST/NIRSpec and JWST/NIRCam data. We examine 163 objects at 1.06 < z < 6.71 with significant (3$\sigma$) spectroscopic detections of H$\alpha$ and H$\beta$ in order to constrain intrinsic H$\alpha$ luminosities corrected from nebular dust attenuation via Balmer decrements. We constrain dust-corrected UV luminosities from best-fit spectral-energy distribution modeling. We find a sample median log$_{10}$($\xi{\rm ion,0}$/erg Hz$^{-1}$) = $25.29^{+0.29}_{-0.37}$, assuming f$_{\rm esc}$=0 for the escape fraction of Lyman continuum emission. We find significant correlation between $\xi_{\rm ion,0}$ and z, with 17 objects at z > 4.64 having median log$_{10}$($\xi_{\rm ion,0}$/erg Hz$^{-1}$) = $25.38^{+0.38}){-0.38}$, with those below having log$_{10}$($\xi_{\rm ion,0}$/erg Hz$^{-1}$) = $25.24^{+0.30}_{-0.33}$. We also find significant, positive correlations between $\xi_{\rm ion,0}$ and LUV; W{\lambda}([O iii]); [O iii]{\lambda}5007/[O ii]{\lambda}{\lambda}3726, 3729; and inverse correlations with metallicity. In contrast with some previous results, we find no trends between $\xi_{\rm ion,0}$ and stellar mass, stellar dust attenuation, or UV slope. Applying a multivariate fit to $\xi_{\rm ion,0}$, z, and MUV to an empirically-motivated model of reionization, and folding in f$_{\rm esc}$ estimates from direct observations of the Lyman continuum at z ~ 3 from the Keck Lyman Continuum Spectroscopic survey, we find that the number of ionizing photons entering the IGM causes reionization to end at z ~ 5 - 7., Comment: 17 pages, 8 figures, submitted to ApJ

We present the stellar metallicities and multi-element abundances (C, Mg, Si, Ca, Ti, Cr, and Fe) of 15 massive (log $M/M_\odot=10.2-11.2$) quiescent galaxies at z=1-3, derived from ultradeep JWST-SUSPENSE spectra. Compared to quiescent galaxies at z~0, these galaxies exhibit a deficiency of 0.26$\pm0.04$ dex in [C/H], 0.16$\pm0.03$ dex in [Fe/H], and 0.07$\pm0.04$ dex in [Mg/H], implying rapid formation and quenching before significant enrichment from asymptotic giant branch stars and Type Ia supernovae. Additionally, we find that galaxies forming at higher redshift consistently show higher [Mg/Fe] and lower [Fe/H] and [Mg/H], regardless of their observed redshift. The evolution in [Fe/H] and [C/H] is therefore primarily driven by lower-redshift samples naturally including galaxies with longer star-formation timescales. In contrast, the lower [Mg/H] likely reflects earlier-forming galaxies expelling larger gas reservoirs during their quenching phase. Consequently, the mass-metallicity relation, primarily reflecting [Mg/H], is somewhat lower at z=1-3 compared to the lower redshift relation. Finally, we compare our results to standard stellar population modeling approaches employing solar abundance patterns and non-parametric star-formation histories (using Prospector). Our SSP-equivalent ages agree with the mass-weighted ages from Prospector, while the metallicities disagree significantly. Nonetheless, the metallicities better reflect [Fe/H] than total [Z/H]. We also find that star-formation timescales inferred from elemental abundances are significantly shorter than those from Prospector, and we discuss the resulting implications for the early formation of massive galaxies., Comment: Accepted to ApJ; 22 pages, 6 figures, 3 tables

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

We investigate the multi-phase structure of gas flows in galaxies. We study 80 galaxies during the epoch of peak star formation ($1.4\leq z\leq2.7$) using data from Keck/LRIS and VLT/KMOS. Our analysis provides a simultaneous probe of outflows using UV emission and absorption features and H$\alpha$ emission. With this unprecedented data set, we examine the properties of gas flows estimated from LRIS and KMOS in relation to other galaxy properties, such as star formation rate (SFR), star formation rate surface density ($\Sigma_{\rm SFR}$), stellar mass (M$_*$), and main sequence offset ($\Delta$MS). We find no strong correlations between outflow velocity measured from rest-UV lines centroids and galaxy properties. However, we find that galaxies with detected outflows show higher averages in SFR, $\Sigma_{\rm SFR}$, and $\Delta$MS than those lacking outflow detections, indicating a connection between outflow and galaxy properties. Furthermore, we find a lower average outflow velocity than previously reported, suggesting greater absorption at the systemic redshift of the galaxy. Finally, we detect outflows in 49% of our LRIS sample and 30% in the KMOS sample, and find no significant correlation between outflow detection and inclination. These results may indicate that outflows are not collimated and that H$\alpha$ outflows have a lower covering fraction than low-ionization interstellar absorption lines. Additionally, these tracers may be sensitive to different physical scales of outflow activity. A larger sample size with a wider dynamic range in galaxy properties is needed to further test this picture., Comment: 20 pages, 10 figures, 6 tables, accepted by ApJ

We have used public JWST/NIRSpec and JWST/NIRCam observations from the CEERS and JADES surveys in order to analyze the star-forming main sequence (SFMS) over the redshift range $1.4 \leq z < 7$. We calculate the star-formation rates (SFRs) of the galaxy sample using three approaches: Balmer line luminosity, spectral energy distribution (SED) fitting, and UV luminosity. We find a larger degree of scatter about the SFMS using the Balmer-based SFRs compared to the UV-based SFRs. Because these SFR indicators are sensitive to star formation on different time scales, the difference in scatter may be evidence of bursty star-formation histories in the early universe. We additionally compare the H$\alpha$-to-UV luminosity ratio (L(H$\alpha$)/$\nu$L$_{\nu,1600}$) for individual galaxies in the sample and find that 29\%$-$52\% of the ratios across the sample are poorly described by predictions from a smooth star-formation history. Measuring the burstiness of star formation in the early universe has multiple significant implications, such as deriving accurate physical parameters from SED fitting, explaining the evolution of the UV luminosity function, and providing constraints for sub-grid models of feedback in simulations of galaxy formation and evolution., Comment: 23 pages, 7 figures, accepted by 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

We present Keck Cosmic Web Imager IFU observations around extended Ly$\alpha$ halos of 27 typical star-forming galaxies with redshifts $2.0 < z < 3.2$ drawn from the MOSFIRE Deep Evolution Field survey. We examine the average Ly$\alpha$ surface-brightness profiles in bins of star-formation rate (SFR), stellar mass ($M_*$), age, stellar continuum reddening, SFR surface density ($\rm \Sigma_{SFR}$), and $\rm \Sigma_{SFR}$ normalized by stellar mass ($\rm \Sigma_{sSFR}$). The scale lengths of the halos correlate with stellar mass, age, and stellar continuum reddening; and anti-correlate with star-formation rate, $\rm \Sigma_{SFR}$, and $\rm \Sigma_{sSFR}$. These results are consistent with a scenario in which the down-the-barrel fraction of Ly$\alpha$ emission is modulated by the low-column-density channels in the ISM, and that the neutral gas covering fraction is related to the physical properties of the galaxies. Specifically, we find that this covering fraction increases with stellar mass, age, and $E(B-V)$; and decreases with SFR, $\rm \Sigma_{SFR}$ and $\rm \Sigma_{sSFR}$. We also find that the resonantly scattered Ly$\alpha$ emission suffers greater attenuation than the (non-resonant) stellar continuum emission, and that the difference in attenuation increases with stellar mass, age, and stellar continuum reddening, and decreases with $\rm \Sigma_{sSFR}$. These results imply that more reddened galaxies have more dust in their CGM., Comment: 12 pages, 7 figures, 2 tables, submitted to ApJ

We present the census of H$\beta$+[OIII] 4960,5008 \AA emitters at 6.87$ galaxies. We report a rapid decline of the [OIII] luminosity density at $z\gtrsim 6-7$ which cannot be explained by the evolution of the cosmic star-formation rate density. Finally we find that FRESCO detects in only 2h galaxies likely accounting for $\sim 10-20\%$ of the ionising budget at $z=7-8$ (assuming an escape fraction of 10%), raising the prospect of directly detecting a significant fraction of the sources of reionisation with JWST., Comment: 20 pages + appendices. Accepted in MNRAS. Public catalogue release at https://github.com/rameyer/fresco. V3: matching accepted version

We report ultra-deep, medium-resolution spectroscopic observations for 4 quiescent galaxies with log$_{10}(M_*/\mathrm{M_\odot})>11$ at $3 < z < 5$. These data were obtained with JWST NIRSpec as part of the Early eXtragalactic Continuum and Emission Line Science (EXCELS) survey, which we introduce in this work. The first two galaxies are newly selected from PRIMER UDS imaging, both at $z=4.62$ and separated by $860$ pkpc on the sky, within a larger structure for which we confirm several other members. Both formed at $z\simeq8-10$. These systems could plausibly merge by the present day to produce a local massive elliptical galaxy. The other two ultra-massive quiescent galaxies are previously known at $z=3.99$ and $3.19$, with the latter (ZF-UDS-7329) having been the subject of debate as potentially too old and too massive to be accommodated by the $\Lambda$-CDM halo-mass function. Both exhibit high stellar metallicities, and for ZF-UDS-7329 we are able to measure the $\alpha-$enhancement, obtaining [Mg/Fe] = $0.42^{+0.19}_{-0.17}$. We finally evaluate whether these 4 galaxies are consistent with the $\Lambda$-CDM halo-mass function using an extreme value statistics approach. We find that the $z=4.62$ objects and the $z=3.19$ object are unlikely within our area under the assumption of standard stellar fractions ($f_*\simeq0.1-0.2$). However, these objects roughly align with the most massive galaxies expected under the assumption of 100 per cent conversion of baryons to stars ($f_*$=1). Our results suggest extreme galaxy formation physics during the first billion years, but no conflict with $\Lambda$-CDM cosmology., Comment: 25 pages, 10 figures, accepted for publication in MNRAS

We present determinations of the gas-phase and stellar metallicities of a sample of 65 star-forming galaxies at $z \simeq 3.5$ using rest-frame far-ultraviolet (FUV) spectroscopy from the VANDELS survey in combination with follow-up rest-frame optical spectroscopy from VLT/KMOS and Keck/MOSFIRE. We infer gas-phase oxygen abundances ($Z_{\mathrm{g}}$; tracing O/H) via strong optical nebular lines and stellar iron abundances ($Z_{\star}$; tracing Fe/H) from full spectral fitting to the FUV continuum. Our sample spans the stellar mass range $8.5 < \mathrm{log}(M_{\star}/\mathrm{M}_{\odot}) < 10.5$ and shows clear evidence for both a stellar and gas-phase mass-metallicity relation (MZR). We find that our O and Fe abundance estimates both exhibit a similar mass-dependence, such that $\mathrm{Fe/H}\propto M_{\star}^{0.30\pm0.11}$ and $\mathrm{O/H}\propto M_{\star}^{0.32\pm0.09}$. At fixed $M_{\star}$ we find that, relative to their solar values, O abundances are systematically larger than Fe abundances (i.e., $\alpha$-enhancement).We estimate an average enhancement of $\mathrm{(O/Fe)} = 2.65 \pm 0.16 \times \mathrm{(O/Fe)_\odot}$ which appears to be independent of $M_{\star}$. We employ analytic chemical evolution models to place a constraint on the strength of galactic-level outflows via the mass-outflow factor ($\eta$). We show that outflow efficiencies that scale as $\eta \propto M_{\star}^{-0.32}$ can simultaneously explain the functional form of of the stellar and gas-phase MZR, as well as the degree of $\alpha$-enhancement at fixed Fe/H. Our results add further evidence to support a picture in which $\alpha$-enhanced abundance ratios are ubiquitous in high-redshift star-forming galaxies, as expected for young systems whose interstellar medium is primarily enriched by core-collapse supernovae., Comment: 18 pages, 10 figures, 4 tables, accepted in MNRAS

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

We use the large spectroscopic data set of the MOSFIRE Deep Evolution Field survey to investigate the kinematics and energetics of ionised gas outflows. Using a sample of 598 star-forming galaxies at redshift 1.4 < $z$ < 3.8, we decompose $\rm{H}\alpha$ and [OIII] emission lines into narrow and broad components, finding significant detections of broad components in 10% of the sample. The ionised outflow velocity from individual galaxies appears independent of galaxy properties, such as stellar mass, star-formation rate (SFR), and star-formation-rate surface density ($\Sigma_{\rm SFR}$). Adopting a simple outflow model, we estimate the mass-, energy- and momentum-loading factors of the ionised outflows, finding modest values with averages of 0.33, 0.04, and 0.22, respectively. The larger momentum- than energy-loading factors, for the adopted physical parameters, imply that these ionised outflows are primarily momentum-driven. We further find a marginal correlation (2.5$\sigma$) between the mass-loading factor and stellar mass in agreement with predictions by simulations, scaling as $\eta_{m}$ $\propto M_{\star}^{-0.45}$. This shallow scaling relation is consistent with these ionised outflows being driven by a combination of mechanical energy generated by supernovae explosions and radiation pressure acting on dusty material. In a majority of galaxies, the outflowing material does not appear to have sufficient velocity to escape the gravitational potential of their host, likely recycling back at later times. Together, these results suggest that the ionised outflows traced by nebular emission lines are negligible, with the bulk of mass and energy carried out in other gaseous phases., Comment: 17 pages, 9 figures, Submitted to Monthly Notices of the Royal Astronomical Society

The exceptional spectra of the most luminous $z>10$ sources observed so far have challenged our understanding of early galaxy evolution, requiring a new observational benchmark for meaningful interpretation. As such, we construct spectroscopic templates representative of high-redshift, star-forming populations, using 482 confirmed sources at $z=5.0-12.9$ with JWST/NIRSpec prism observations, and report on their average properties. We find $z=5-11$ galaxies are dominated by blue UV continuum slopes ($\beta=-2.3$ to $-2.7$) and reduced Balmer indices, characteristic of dust-poor and young systems, with a shift towards bluer slopes and younger ages with redshift. The evolution is mirrored by ubiquitous CIII] detections across all redshifts (rest-frame equivalent widths of $=5-14$ \r{A}), which increase in strength towards early times. Rest-frame optical lines reveal elevated ratios ($O32=7-31$, $R23=5-8$, and $Ne3O2=1-2$) and subsolar metallicities (log O/H$=7.3-7.9$), typical of ionization conditions and metallicities rarely observed in $z\sim0$ populations. Within our sample, we identify 57 Ly$\alpha$ emitters, which we stack and compare to a matched sample of nonemitters. The former are characterized by more extreme ionizing conditions with enhanced CIII], CIV, and HeII+[OIII] line emission, younger stellar populations from Balmer jumps, and a more pristine interstellar medium seen through bluer UV slopes and elevated rest-frame optical line ratios. The novel comparison illustrates important intrinsic differences between the two populations, with implications for Ly$\alpha$ visibility. The spectral templates derived here represent a new observational benchmark with which to interpret high-redshift sources, lifting our constraints on their global properties to unprecedented heights and extending out to the earliest of cosmic times., Comment: Accepted for publication in ApJ (with minor changes). Replaced to match accepted version. All composite spectra available online or upon request

High-redshift Lyman break galaxies (LBGs) are efficiently selected in deep images using as few as three broadband filters, and have been shown to have multiple intrinsic and small- to large-scale environmental properties related to Lyman-alpha. In this paper we demonstrate a statistical relationship between net Lyman-alpha equivalent width (net Lya EW) and the optical broadband photometric properties of LBGs at z~2. We show that LBGs with the strongest net Lya EW in absorption (aLBGs) and strongest net Lya EW in emission (eLBGs) separate into overlapping but discrete distributions in $(U_n-R)$ colour and $R$-band magnitude space, and use this segregation behaviour to determine photometric criteria by which sub-samples with a desired Lya spectral type can be selected using data from as few as three broadband optical filters. We propose application of our result to current and future large-area and all-sky photometric surveys that will select hundreds of millions of LBGs across many hundreds to thousands of Mpc, and for which spectroscopic follow-up to obtain Lya spectral information is prohibitive. To this end, we use spectrophotometry of composite spectra derived from a sample of 798 LBGs divided into quartiles on the basis of net Lya EW to calculate criteria for the selection of Lya absorbing and Lya emitting populations of z~3 LBGs using $ugri$ broadband photometric data from the Vera Rubin Observatory Legacy Survey of Space and Time (LSST)., Comment: Published in Publications of the Astronomical Society of Australia 16 pages, 6 figures, 8 tables

The ratios of strong rest-frame optical emission lines are the dominant indicator 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 HII regions, post-AGB 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 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 $\sim$ 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., Comment: 28 pages, 6 figures, typo corrected: in polynomial fits, x represents metallicity and y represents the line ratio

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

The chemical enrichment of dust and metals are vital processes in constraining the star formation history of the universe. Previously, the dust masses of high-redshift star-forming galaxies have been determined through their far-infrared continuum, however, equivalent, and potentially simpler, approaches to determining the metal masses have yet to be explored at $z\gtrsim 2$. Here, we present a new method of inferring the metal mass in the interstellar medium (ISM) of galaxies out to $z\approx 8$, using the far-infrared [CII]$-158\mu$m emission line as a proxy. We calibrated the [CII]-to-$M_{\rm Z,ISM}$ conversion factor based on a benchmark observational sample at $z\approx 0$, in addition to gamma-ray burst sightlines at $z>2$ and cosmological hydrodynamical simulations of galaxies at $z\approx 0$ and $z\approx 6$. We found a universal scaling across redshifts of $\log (M_{\rm Z,ISM}/M_\odot) = \log (L_{\rm [CII]}/L_\odot) - 0.45,$ with a 0.4 dex scatter, which is constant over more than two orders of magnitude in metallicity. We applied this scaling to recent surveys for [CII] in galaxies at $z\gtrsim 2$ and determined the fraction of metals retained in the gas-phase ISM, $M_{\rm Z,ISM} / M_\star$, as a function of redshift showing that an increasing fraction of metals reside in the ISM of galaxies at higher redshifts. We place further constraints on the cosmic metal mass density in the ISM ($\Omega_{\rm Z,ISM}$) at $z\approx 5$ and $\approx 7$, yielding $\Omega_{\rm Z,ISM} = 6.6^{+13}_{-4.3}\times 10^{-7}\,M_\odot\, {\rm Mpc}^{-3}$ ($z\approx 5$) and $\Omega_{\rm Z,ISM} = 2.0^{+3.5}_{-1.3}\times 10^{-7}\,M_\odot\, {\rm Mpc}^{-3}$ ($z\approx 7$). These results are consistent with the expected metal yields from the integrated star formation history at the respective redshifts. This suggests that the majority of metals produced at $z\gtrsim 5$ are confined to the ISM of galaxies., Comment: Accepted in A&A, abstract abridged

We present a study of the connection between the escape fraction of Lyman-alpha (Ly$\alpha$) and Lyman-continuum (LyC) photons within a sample of N=152 star-forming galaxies selected from the VANDELS survey at $3.85

We report on the discovery of cool gas inflows towards three star-forming galaxies at $\left\sim$ 2.30. Analysis of Keck Low-Resolution Imaging Spectrometer spectroscopy reveals redshifted low-ionisation interstellar (LIS) metal absorption lines with centroid velocities of 60 - 130 km $\rm{s}^{-1}$. These inflows represent some of the most robust detections of inflowing gas into isolated, star-forming galaxies at high redshift. Our analysis suggests that the inflows are due to recycling metal-enriched gas from previous ejections. Comparisons between the galaxies with inflows and a larger parent sample of 131 objects indicate that galaxies with detected inflows may have higher specific star-formation rates (sSFR) and star-formation-rate surface densities. However, when additional galaxies without robustly detected inflows based on centroid velocity but whose LIS absorption line profiles indicate large red-wing velocities are considered, galaxies with inflows do not show unique properties relative to those lacking inflows. Additionally, we calculate the covering fraction of cool inflowing gas as a function of red-wing inflow velocity, finding an enhancement in high sSFR binned galaxies, likely due to an increase in the amount of recycling gas. Together, these results suggest that the low detection rate of galaxies with cool inflows is primarily related to the viewing angle rather than the physical properties of the galaxies., Comment: 10 pages, 7 figures, accepted for publication in MNRAS

We investigate the production efficiency of ionizing photons ($\xi_{ion}^*$) of 1174 galaxies with secure redshift at z=2-5 from the VANDELS survey to determine the relation between ionizing emission and physical properties of bright and massive sources. We constrain $\xi_{ion}^*$ and galaxy physical parameters by means of spectro-photometric fits performed with the BEAGLE code. The analysis exploits the multi-band photometry in the VANDELS fields, and the measurement of UV rest-frame emission lines (CIII]$\lambda 1909$, HeII$\lambda 1640$, OIII]$\lambda 1666$) from deep VIMOS spectra. We find no clear evolution of $\xi_{ion}^*$ with redshift within the probed range. The ionizing efficiency slightly increases at fainter $M_{UV}$, and bluer UV slopes, but these trends are less evident when restricting the analysis to a complete subsample at log(M$_{star}$/M$_{\odot}$)$>$9.5. We find a significant trend of increasing $\xi_{ion}^*$ with increasing EW(Ly$\alpha$), with an average log($\xi_{ion}^*$/Hz erg$^{-1}$)$>$25 at EW$>$50\AA, and a higher ionizing efficiency for high-EW CIII]$\lambda 1909$ and OIII]$\lambda 1666$ emitters. The most significant correlations are found with respect to stellar mass, specific star-formation rate (sSFR) and SFR surface density ($\Sigma_{SFR}$). The relation between $\xi_{ion}^*$ and sSFR shows a monotonic increase from log($\xi_{ion}^*$/Hz erg$^{-1}$) $\sim$24.5 at log(sSFR)$\sim$-9.5$yr^{-1}$ to $\sim$25.5 at log(sSFR)$\sim$-7.5$yr^{-1}$, a low scatter and little dependence on mass. The objects above the main-sequence of star-formation consistently have higher-than-average $\xi_{ion}^*$. A clear increase of $\xi_{ion}^*$ with $\Sigma_{SFR}$ is also found, with log($\xi_{ion}^*$/Hz erg$^{-1}$)$>$25 for objects at $\Sigma_{SFR}>$10 M$_{\odot}/yr/kpc^2$.(Abridged), Comment: 11 pages, 9 figures, accepted for publication in A&A

Using deep near-infrared Keck/MOSFIRE observations, we analyze the rest-optical spectra of eight star-forming galaxies in the COSMOS and GOODS-N fields. We reach integration times of $\sim$10 hours in the deepest bands, pushing the limits on current ground-based observational capabilities. The targets fall into two redshift bins -- 5 galaxies at $z \sim 1.7$ and 3 at $z \sim 2.5$ -- and were selected as likely to yield significant auroral-line detections. Even with long integration times, detection of the auroral lines remains challenging. We stack the spectra together into subsets based on redshift, improving the signal-to-noise ratio on the [O III] $\lambda 4364$ auroral emission line and, in turn, enabling a direct measurement of the oxygen abundance for each stack. We compare these measurements to commonly-employed strong-line ratios alongside measurements from the literature. We find that the stacks fall within the distribution of $z>1$ literature measurements, but a larger sample size is needed to robustly constrain the relationships between strong-line ratios and oxygen abundance at high redshift. We additionally report detections of [O I] $\lambda6302$ for nine individual galaxies and composite spectra of 21 targets in the MOSFIRE pointings. We plot their line ratios on the [O III] $\lambda 5008$/H$\beta$ vs. [O I] $\lambda 6302$/H$\alpha$ diagnostic BPT diagram, comparing our targets to local galaxies and H II regions. We find that the [O I]/H$\alpha$ ratios in our sample of galaxies are consistent with being produced in gas ionized by $\alpha$-enhanced massive stars, as has been previously inferred for rapidly-forming galaxies at early cosmic times., Comment: 18 pages, 6 figures, submitted to ApJ

We investigate dust attenuation and its dependence on viewing angle for 308 star-forming galaxies at $1.3\leq z\leq2.6$ from the MOSFIRE Deep Evolution Field (MOSDEF) survey. We divide galaxies with a detected H$\alpha$ emission line and coverage of H$\beta$ into eight groups by stellar mass, star formation rate (SFR), and inclination (i.e., axis ratio), then stack their spectra. From each stack, we measure Balmer decrement and gas-phase metallicity, then we compute median \AV and UV continuum spectral slope ($\beta$). First, we find that none of the dust properties (Balmer decrement, \AV, $\beta$) vary with axis ratio. Second, both stellar and nebular attenuation increase with increasing galaxy mass, showing little residual dependence on SFR or metallicity. Third, nebular emission is more attenuated than stellar emission, and this difference grows even larger at higher galaxy masses and SFRs. Based on these results, we propose a three-component dust model where attenuation predominantly occurs in star-forming regions and large, dusty star-forming clumps, with minimal attenuation in the diffuse ISM. In this model, nebular attenuation primarily originates in clumps, while stellar attenuation is dominated by star-forming regions. Clumps become larger and more common with increasing galaxy mass, creating the above mass trends. Finally, we argue that a fixed metal yield naturally leads to mass regulating dust attenuation. Infall of low-metallicity gas increases SFR and lowers metallicity, but leaves dust column density mostly unchanged. We quantify this idea using the Kennicutt-Schmidt and fundamental metallicity relations, showing that galaxy mass is indeed the primary driver of dust attenuation., Comment: 19 pages, 13 figures, accepted for publication in ApJ

We present the JWST Cycle 1 53.8hr medium program FRESCO, short for "First Reionization Epoch Spectroscopically Complete Observations". FRESCO covers 62 arcmin$^2$ in each of the two GOODS/CANDELS fields for a total area of 124 arcmin$^2$ exploiting JWST's powerful new grism spectroscopic capabilities at near-infrared wavelengths. By obtaining ~2 hr deep NIRCam/grism observations with the F444W filter, FRESCO yields unprecedented spectra at R~1600 covering 3.8 to 5.0 $\mu$m for most galaxies in the NIRCam field-of-view. This setup enables emission line measurements over most of cosmic history, from strong PAH lines at z~0.2-0.5, to Pa$\alpha$ and Pa$\beta$ at z~1-3, HeI and [SIII] at z~2.5-4.5, H$\alpha$ and [NII] at z~5-6.5, up to [OIII] and H$\beta$ for z~7-9 galaxies, and possibly even [OII] at z~10-12. FRESCO's grism observations provide total line fluxes for accurately estimating galaxy stellar masses and calibrating slit-loss corrections of NIRSpec/MSA spectra in the same field. Additionally, FRESCO results in a mosaic of F182M, F210M, and F444W imaging in the same fields to a depth of ~28.2 mag (5 $\sigma$ in 0.32" diameter apertures). Together with this publication, the v1 imaging mosaics are released as high-level science products via MAST. Here, we describe the overall survey design and the key science goals that can be addressed with FRESCO. We also highlight several, early science results, including: spectroscopic redshifts of Lyman break galaxies that were identified almost 20 years ago, the discovery of broad-line active galactic nuclei at z>4, and resolved Pa$\alpha$ maps of galaxies at z~1.4. These results demonstrate the enormous power for serendipitous discovery of NIRCam/grism observations., Comment: 12 pages, 11 figures; MNRAS in press; for more information on the survey and data releases, see http://jwst-fresco.astro.unige.ch/ and https://archive.stsci.edu/hlsp/fresco

We conduct a statistical analysis of the factors responsible for the variation in the ionization parameter (U) of high-redshift star-forming galaxies based on medium resolution JWST/NIRSpec observations obtained by the Cosmic Evolution Early Release Science (CEERS) survey. The sample consists of 48 galaxies with spectroscopic redshifts z=2.7-6.3 which are largely representative of typical star-forming galaxies at these redshifts. The [SII] 6718, 6733 doublet is used to estimate electron densities (n_e), and dust-corrected Ha luminosities are used to compute total ionizing photon rates (Q). Using composite spectra of galaxies in bins of [OIII]/[OII] (i.e., O32) as a proxy for U, we determine that galaxies with higher O32 have ~ 500 cm^-3 that are at least a factor of ~5 larger than that of lower-O32 galaxies. We do not find a significant difference in between low- and high-O32 galaxies. Photoionization modeling of all available strong rest-frame optical emission lines is used to simultaneously constrain U and oxygen abundance (Z_neb). We find a large spread in log U of ~1.5 dex at a fixed Z_neb. On the other hand, the data indicate a highly significant correlation between U and star-formation-rate surface density (Sigma_SFR) which appears to be redshift invariant at z~1.6-6.3, and possibly up to z~9.5. We consider several avenues through which metallicity and Sigma_SFR (or gas density) may influence U, including variations in n_e and Q that are tied to metallicity and gas density, internal dust extinction of ionizing photons, and the effects of gas density on the volume filling fraction of dense clumps in HII regions and the escape fraction of ionizing photons. Based on these considerations, we conclude that gas density may play a more central role than metallicity in modulating U at these redshifts., Comment: 16 pages, 13 figures, 4 tables

We report detections of the [OIII]$\lambda$4364 auroral emission line for 16 galaxies at z=2.1-8.7, measured from JWST/NIRSpec observations obtained as part of the Cosmic Evolution Early Release Science (CEERS) survey program. We combine this CEERS sample with 9 objects from the literature at z=4-9 with auroral-line detections from JWST/NIRSpec and 21 galaxies at z=1.4-3.7 with auroral-line detections from ground-based spectroscopy. We derive electron temperature T_e and direct-method oxygen abundances for the combined sample of 46 star-forming galaxies at z=1.4-8.7. We use these measurements to construct the first high-redshift empirical T_e-based metallicity calibrations for the strong-line ratios [OIII]/H$\beta$, [OII]/H$\beta$, R23=([OIII]+[OII])/H$\beta$, [OIII]/[OII], and [NeIII]/[OII]. These new calibrations are valid over 12+log(O/H)=7.0-8.4 and can be applied to samples of star-forming galaxies at z=2-9, leading to an improvement in the accuracy of metallicity determinations at Cosmic Noon and in the Epoch of Reionization. The high-redshift strong-line relations are offset from calibrations based on typical $z\sim0$ galaxies or HII regions, reflecting the known evolution of ionization conditions between $z\sim0$ and $z\sim2$. Deep spectroscopic programs with JWST/NIRSpec promise to improve statistics at the low and high ends of the metallicity range covered by the current sample, as well as improve the detection rate of [NII]$\lambda$6585 to allow the future assessment of N-based indicators. These new high-redshift calibrations will enable accurate characterizations of metallicity scaling relations at high redshift, improving our understanding of feedback and baryon cycling in the early universe., Comment: 22 pages, 7 figures, 3 tables, submitted to ApJ

We selected 35 low-mass SFGs (7.91000\r{A} and EW(CIII])>5\r{A}, closely resembling those now seen in EoR galaxies with JWST. We find low gas-phase metallicities 12+log(O/H)~7.5-8.5 and C/O abundances from 23%-128% solar, with no apparent increasing trend with metallicity. From our [OIII]4959,5007 profile modeling, we find that 65% of our sample shows an outflow component, which is shifted relative to the ionized gas systemic velocity, with mean $v_{max}$~280 km/s which correlates with the $\Sigma_{SFR}$. We find that the mass-loading factor $\mu$ of our sample is typically lower than in more massive galaxies from literature but higher than in typical local dwarf galaxies. In the stellar mass range covered, we find that $\mu$ increases with $\Sigma_{SFR}$ thus suggesting that for a given stellar mass, denser starbursts in low-mass galaxies produce stronger outflows. Our results complement the picture drawn by similar studies at lower redshift, suggesting that the removal of ionized gas in low-mass SFGs driven by stellar feedback is regulated by their stellar mass and by the strength and concentration of their star formation, i.e. $\Sigma_{\rm SFR}$., Comment: Accepted for publication in A&A. We updated the manuscript following referee's suggestions. We updated the estimations of C/O by not including CIV flux since it is not detected in most of our sample. No major changes in our results compared with previous version

We analyze the rest-optical emission-line ratios of star-forming galaxies at 2.7<=z<6.5 drawn from the Cosmic Evolution Early Release Science (CEERS) Survey, and their relationships with stellar mass (M_*). Our analysis includes both line ratios based on the [NII]6583 feature -- [NII]6583/Ha, ([OIII]5007/Hb)/([NII]6583/Ha) (O3N2), and [NII]6583/[OII]3727 -- and those those featuring alpha elements -- [OIII]5007/Hb, [OIII]5007/[OII]3727 (O_32), ([OIII]4959,5007+[OII]3727)/Hb (R_23), and [NeIII]3869/[OII]3727. Given the typical flux levels of [NII]6583 and [NeIII]3869, which are undetected in the majority of individual CEERS galaxies at 2.7<=z<6.5, we construct composite spectra in bins of M_* and redshift. Using these composite spectra, we compare the relationships between emission-line ratios and M_* at 2.7<=z<6.5 with those observed at lower redshift. While there is significant evolution towards higher excitation (e.g., higher [OIII]5007/Hb, O_32, O3N2), and weaker nitrogen emission (e.g., lower [NII]6583/Ha and [NII]6583/[OII]3727) between z~0 and z~3, we find in most cases that there is no significant evolution in the relationship between line ratio and M_* beyond z~3. The [NeIII]3869/[OII]3727 ratio is anomalous in showing evidence for significant elevation at 4.0<=z<6.5 at fixed mass, relative to z~3.3. Collectively, however, our empirical results suggest that there is no significant evolution in the mass-metallicity relationship at 2.7<=z<6.5. Representative galaxy samples and metallicity calibrations based on existing and upcoming JWST/NIRSpec observations will be required to translate these empirical scaling relations into ones tracing chemical enrichment and gas cycling, and distinguish among the descriptions of star-formation feedback in simulations of galaxy formation at z>3., Comment: 10 pages, 4 figures, ApJL, in press

We use the large spectroscopic dataset of the MOSFIRE Deep Evolution Field (MOSDEF) survey to investigate some of the key factors responsible for the elevated ionization parameters (U) inferred for high-redshift galaxies, focusing in particular on the role of star-formation-rate surface density (Sigma_SFR). Using a sample of 317 galaxies with spectroscopic redshifts z~1.9-3.7, we construct composite rest-frame optical spectra in bins of Sigma_SFR and infer electron densities, n_e, using the ratio of the [OII] 3727, 3730 doublet. Our analysis suggests a significant (~3 sigma) correlation between n_e and Sigma_SFR. We further find significant correlations between U and Sigma_SFR for composite spectra of a subsample of 113 galaxies, and for a smaller sample of 25 individual galaxies with inferences of U. The increase in n_e -- and possibly also the volume filling factor of dense clumps in HII regions -- with Sigma_SFR appear to be important factors in explaining the relationship between U and Sigma_SFR. Further, the increase in n_e and SFR with redshift at a fixed stellar mass can account for most of the redshift evolution of U. These results suggest that the gas density, which sets n_e and the overall level of star-formation activity, may play a more important role than metallicity evolution in explaining the elevated ionization parameters of high-redshift galaxies., Comment: 11 pages, 9 figures, 1 table, submitted to ApJ

We use medium resolution JWST/NIRSpec observations from the Cosmic Evolution Early Release Science (CEERS) Survey to place the first constraints on dust attenuation and star formation based on the Paschen lines for a sizable sample of 63 galaxies at redshifts z=1.0-3.1. Our analysis indicates strong correlations between the Balmer decrement, Ha/Hb, and line ratios that include the Paschen lines (i.e., Paa/Hb, Pab/Hb, and the Paschen decrement, Paa/Pab), suggesting that the former is sensitive to the overall dust obscuration towards HII regions in high-redshift galaxies. The line ratios are used to derive the nebular reddening, E(B-V)neb, and star-formation rates (SFRs). There is marginal evidence that SFRs deduced from the Paschen lines may exceed by ~25% those derived from the Balmer lines alone, suggesting the presence of star formation that is optically thick in the Balmer lines, though deeper observations are needed to confirm this result. Using the Paschen-line constraints on bolometric SFRs, we reevaluate the relationship between dust obscuration and UV spectral slope, and find a reddening of the UV continuum that, on average, follows the SMC extinction curve. This analysis highlights the need for deeper spectroscopy of more representative samples to evaluate nebular dust attenuation and bolometric SFRs in high-redshift galaxies, and their relationship to the reddening of the UV continuum., Comment: 9 pages, 8 figures, 1 table

We utilize medium-resolution JWST/NIRSpec observations of 164 galaxies at $z=2.0-9.3$ from the Cosmic Evolution Early Release Science (CEERS) survey to investigate the evolution of the excitation and ionization properties of galaxies at high redshifts. Our results represent the first statistical constraints on the evolution of the [OIII]/H$\beta$ vs. [NII]/H$\alpha$, [SII]/H$\alpha$, and [OI]/H$\alpha$ ``BPT'' diagrams at $z>2.7$, and the first analysis of the O32 vs. R23 diagram at $z>4$ with a large sample. We divide the sample into five redshift bins containing 30-40 galaxies each. The subsamples at $z\sim2.3$, $z\sim3.3$, and $z\sim4.5$ are representative of the main-sequence star-forming galaxy population at these redshifts, while the $z\sim5.6$ and $z\sim7.5$ samples are likely biased toward high specific star-formation rate due to selection effects. Using composite spectra, we find that each subsample at $z=2.0-6.5$ falls on the same excitation sequence in the [NII] and [SII] BPT diagrams and the O32-R23 diagram on average, offset from the sequences followed by $z=0$ HII regions in the same diagrams. The direction of these offsets are consistent with high-redshift star-forming galaxies uniformly having harder ionizing spectra than typical local galaxies at fixed nebular metallicity. The similarity of the average line ratios suggests that the ionization conditions of the interstellar medium do not strongly evolve between $z\sim2$ and $z\sim6$. Overall, the rest-optical line ratios suggest the $z=2.7-9.3$ CEERS/NIRSpec galaxies at log($M_*/M_{\odot})\sim7.5-10$ have high degrees of ionization and moderately low oxygen abundances ($\sim0.1-0.3~Z_{\odot}$), but are not extremely metal poor ($<0.1~Z_{\odot}$) even at $z>6.5$., Comment: 15 pages, 6 figures, submitted to ApJ

We present an analysis of the star-formation rates (SFRs) and dust attenuation properties of star-forming galaxies at $2.7\leq z<6.5$ drawn from the Cosmic Evolution Early Release Science (CEERS) Survey. Our analysis is based on {\it JWST}/NIRSpec Micro-Shutter Assembly (MSA) $R\sim1000$ spectroscopic observations covering approximately $1-5$$\mu$m. Our primary rest-frame optical spectroscopic measurements are H$\alpha$/H$\beta$ Balmer decrements, which we use as an indicator of nebular dust attenuation. In turn, we use Balmer decrements to obtain dust-corrected H$\alpha$-based SFRs (i.e., SFR(H$\alpha$)). We construct the relationship between SFR(H$\alpha$) and stellar mass ($M_*$) in three bins of redshift ($2.7\leq z< 4.0$, $4.0\leq z< 5.0$, and $5.0\leq z<6.5$), which represents the first time the star-forming main sequence has been traced at these redshifts using direct spectroscopic measurements of Balmer emission as a proxy for SFR. In tracing the relationship between SFR(H$\alpha$) and $M_*$ back to such early times ($z>3$), it is essential to use a conversion factor between H$\alpha$ and SFR that accounts for the subsolar metallicity prevalent among distant galaxies. We also use measured Balmer decrements to investigate the relationship between dust attenuation and stellar mass out to $z\sim6$. The lack of significant redshift evolution in attenuation at fixed stellar mass, previously confirmed using Balmer decrements out to $z\sim2.3$, appears to hold out to $z\sim 6.5$. Given the rapidly evolving gas, dust, and metal content of star-forming galaxies at fixed mass, this lack of significant evolution in attenuation provides an ongoing challenge to explain., Comment: 9 pages, 4 figures, ApJ, in press