Search

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

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

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

In November 2022, the James Webb Space Telescope (JWST) returned deep near-infrared images of Abell 2744 -- a powerful lensing cluster capable of magnifying distant, incipient galaxies beyond it. Together with the existing Hubble Space Telescope (HST) imaging, this publicly available dataset opens a fundamentally new discovery space to understand the remaining mysteries of the formation and evolution of galaxies across cosmic time. In this work, we detect and measure some 60,000 objects across the 49 arcmin$^2$ JWST footprint down to a $5\,\sigma$ limiting magnitude of $\sim$30 mag in 0.32" apertures. Photometry is performed using circular apertures on images matched to the point spread function of the reddest NIRCam broad band, F444W, and cleaned of bright cluster galaxies and the related intra-cluster light. To give an impression of the photometric performance, we measure photometric redshifts and achieve a $\sigma_{\rm NMAD}\approx0.03$ based on known, but relatively small, spectroscopic samples. With this paper, we publicly release our HST and JWST PSF-matched photometric catalog with optimally assigned aperture sizes for easy use, along with single aperture catalogs, photometric redshifts, rest-frame colors, and individual magnification estimates. These catalogs will set the stage for efficient and deep spectroscopic follow-up of some of the first JWST-selected samples in Summer 2023., Comment: 28 pages, 19 figures, resubmitted to ApJS following significant data product improvements. Comments welcome. Catalogs can be accessed at https://jwst-uncover.github.io/DR2.html#PhotometricCatalogs

We present a joint analysis of the galaxy S04590 at $z=8.496$ based on NIRSpec, NIRCam, and NIRISS observations obtained through as part of Early Release Observations programme of the James Webb Space Telescope (JWST) and the far-infrared [CII]-$158\mu$m emission line detected by dedicated Atacama Large Millimeter/submillimeter Array (ALMA) observations. We determine the physical properties of S04590 from modelling of the spectral energy distribution (SED) and through the redshifted optical nebular emission lines detected with JWST/NIRSpec. The best-fit SED model reveals a low-mass ($M_\star = 10^{7.2}-10^{8}\,M_{\odot}$) galaxy with a low oxygen abundance of $12+\log{\rm (O/H)} = 7.16^{+0.10}_{-0.12}$ derived from the strong nebular and auroral emission lines. Assuming that [CII] effectively traces the interstellar medium (ISM), we estimate the total gas mass of the galaxy to be $M_{\rm gas} = (8.0\pm 4.0)\times 10^{8}\,M_\odot$ based on the luminosity and spatial extent of [CII]. This yields an exceptionally high gas fraction, $f_{\rm gas} = M_{\rm gas}/(M_{\rm gas} + M_\star) \gtrsim 90\%$, though still consistent within the range expected for its low metallicity. We further derive the metal mass of the galaxy based on the gas mass and gas-phase metallicity, which we find to be consistent with the expected metal production from Type II supernovae. Finally, we make the first constraints on the dust-to-gas (DTG) and dust-to-metals (DTM) ratios of galaxies in the epoch of reionization at $z\gtrsim 6$, showing overall low mass ratios of logDGT $<-3.8$ and logDTM $<-0.5$, though consistent with local scaling relations and in particular the local metal-poor galaxy I Zwicky 18. Our analysis highlights the synergy between ALMA and JWST in characterizing the gas, metal, and stellar content of the first generation of galaxies., Comment: Accepted for publication in ApJL

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

The $\rm{H}\alpha$-to-UV luminosity ratio ($L(\rm H\alpha)/L(\rm UV)$) is often used to probe SFHs of star-forming galaxies and it is important to validate it against other proxies for burstiness. To address this issue, we present a statistical analysis of the resolved distribution of $\Sigma_{\rm{SFR}}$ as well as stellar age and their correlations with the globally measured $L(\rm H\alpha)/L(\rm UV)$ for a sample of 310 star-forming galaxies in two redshift bins of $1.37 < z < 1.70$ and $ 2.09 < z < 2.61$ observed by the MOSDEF survey. We use the multi-waveband CANDELS/3D-HST imaging of MOSDEF galaxies to construct $\Sigma_{\rm{SFR}}$ and stellar age maps. We analyze the composite rest-frame far-UV spectra of a subsample of MOSDEF targets obtained by the Keck/LRIS, which includes 124 star-forming galaxies (MOSDEF-LRIS) at redshifts $1.4 < z < 2.6$, to examine the average stellar population properties, and the strength of age-sensitive FUV spectral features in bins of $L(\rm H\alpha)/L(\rm UV)$. Our results show no significant evidence that individual galaxies with higher $L(\rm H\alpha)/L(\rm UV)$ are undergoing a burst of star formation based on the resolved distribution of $\Sigma_{\rm{SFR}}$ of individual star-forming galaxies. We segregate the sample into subsets with low and high $L(\rm H\alpha)/L(\rm UV)$. The high-$L(\rm H\alpha)/L(\rm UV)$ subset exhibits, on average, an age of $\log[\rm{Age/yr}]$ = 8.0, compared to $\log[\rm{Age/yr}]$ = 8.4 for the low-$L(\rm H\alpha)/L(\rm UV)$ galaxies, though the difference in age is significant at only the $2\sigma$ level. Furthermore, we find no variation in the strengths of Siiv$\lambda\lambda1393, 1402$ and Civ$\lambda\lambda1548, 1550$ P-Cygni features from massive stars between the two subsamples., Comment: 16 pages, 10 figures, published by the Monthly Notices of the Royal Astronomical Society

We define a new morphology metric called "patchiness" ($P$) that is sensitive to deviations from the average of a resolved distribution, does not require the galaxy center to be defined, and can be used on the spatially-resolved distribution of any galaxy property. While the patchiness metric has a broad range of applications, we demonstrate its utility by investigating the distribution of dust in the interstellar medium of 310 star-forming galaxies at spectroscopic redshifts $1.36

We present ultra-deep Keck/MOSFIRE rest-optical spectra of two star-forming galaxies at z=2.18 in the COSMOS field with bright emission lines, representing more than 20~hours of total integration. The fidelity of these spectra enabled the detection of more than 20 unique emission lines for each galaxy, including the first detection of the auroral [O II]$\lambda\lambda$7322,7332 lines at high redshift. We use these measurements to calculate the electron temperature in the low-ionization O$^+$ zone of the ionized ISM and derive abundance ratios of O/H, N/H, and N/O using the direct method. The N/O and $\alpha$/Fe abundance patterns of these galaxies are consistent with rapid formation timescales and ongoing strong starbursts, in accord with their high specific star-formation rates. These results demonstrate the feasibility of using auroral [O II] measurements for accurate metallicity studies at high redshift in a higher metallicity regime previously unexplored with the direct method in distant galaxies. These results also highlight the difficulty in obtaining the measurements required for direct-method metallicities from the ground. We emphasize the advantages that the JWST/NIRSpec instrument will bring to high-redshift metallicity studies, where the combination of increased sensitivity and uninterrupted wavelength coverage will yield more than an order of magnitude increase in efficiency for multiplexed auroral-line surveys relative to current ground-based facilities. Consequently, the advent of JWST promises to be the beginning of a new era of precision chemical abundance studies of the early universe at a level of detail rivaling that of local galaxy studies., Comment: 18 pages, 5 figures, 3 tables, submitted to ApJ

We analyze the completeness of the MOSDEF survey, in which z ~ 2 galaxies were selected for rest-optical spectroscopy from well-studied HST extragalactic legacy fields down to a fixed rest-optical magnitude limit (H_AB = 24.5). The subset of z ~ 2 MOSDEF galaxies with high signal-to-noise (S/N) emission-line detections analyzed in previous work represents a small minority (<10%) of possible z ~ 2 MOSDEF targets. It is therefore crucial to understand how representative this high S/N subsample is, while also more fully exploiting the MOSDEF spectroscopic sample. Using spectral-energy-distribution (SED) models and rest-optical spectral stacking, we compare the MOSDEF z ~ 2 high S/N subsample with the full MOSDEF sample of z ~ 2 star-forming galaxies with redshifts, the latter representing an increase in sample size of more than a factor of three. We find that both samples have similar emission-line properties, in particular in terms of the magnitude of the offset from the local star-forming sequence on the [N II] BPT diagram. There are small differences in median host galaxy properties, including the stellar mass (M_*), star-formation rate (SFR) and specific SFR (sSFR), and UVJ colors; however, these offsets are minor considering the wide spread of the distributions. Using SED modeling, we also demonstrate that the sample of z ~ 2 star-forming galaxies observed by the MOSDEF survey is representative of the parent catalog of available such targets. We conclude that previous MOSDEF results on the evolution of star-forming galaxy emission-line properties were unbiased relative to the parent z ~ 2 galaxy population., Comment: 18 pages, 10 figures, 3 tables

We present a detailed study of a galaxy merger taking place at $z=1.89$ in the GOODS-S field. Here we analyze Keck/MOSFIRE spectroscopic observations from the MOSFIRE Deep Evolution Field (MOSDEF) survey along with multi-wavelength photometry assembled by the 3D-HST survey. The combined dataset is modeled to infer the past star-formation histories (SFHs) of both merging galaxies. They are found to be massive, with log$_{10}(M_{\ast}/M_{\odot}) > 11$, with a close mass ratio satisfying the typical major-merger definition. Additionally, in the context of delayed-$\tau$ models, GOODS-S 43114 and GOODS-S 43683 have similar SFHs and low star-formation rates (log$_{10}$(SFR(SED)/$M_{\odot}/\rm{yr}^{-1}$) $<$ 1.0) compared to their past averages. The best-fit model SEDs show elevated H$\delta_{\rm{A}}$ values for both galaxies, indicating that their stellar spectra are dominated by A-type stars, and that star formation peaked $\sim0.5-1$ Gyr ago and has recently declined. Additionally, based on SED fitting both merging galaxies turned on and shut off star formation within a few hundred Myr of each other, suggesting that their bursts of star formation may be linked. Combining the SFHs and H$\delta_{\rm{A}}$ results with recent galaxy merger simulations, we infer that these galaxies have recently completed their first pericentric passage and are moving apart. Finally, the relatively low second velocity moment of GOODS-S 43114 given its stellar mass, suggests a disk-like structure. However, including the geometry of the galaxy in the modeling does not completely resolve the discrepancy between the dynamical and stellar masses. Future work is needed to resolve this inconsistency in mass., Comment: 12 pages, 7 figures, 1 table, Accepted for publication in MNRAS

We aim to quantify the relation between the dust-to-gas mass ratio (DTG) and gas-phase metallicity of $z=$2.1-2.5 luminous galaxies and contrast this high-redshift relation against analogous constraints at z$=$0. We present a sample of ten star-forming main-sequence galaxies in the redshift range $2.1

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

The Baldwin, Philips, & Terlevich diagram of [O III]/H$\beta$ vs. [N II]/H$\alpha$ (hereafter N2-BPT) has long been used as a tool for classifying galaxies based on the dominant source of ionizing radiation. Recent observations have demonstrated that galaxies at $z\sim2$ reside offset from local galaxies in the N2-BPT space. In this paper, we conduct a series of controlled numerical experiments to understand the potential physical processes driving this offset. We model nebular line emission in a large sample of galaxies, taken from the SIMBA cosmological hydrodynamic galaxy formation simulation, using the CLOUDY photoionization code to compute the nebular line luminosities from H II regions. We find that the observed shift toward higher [O III]/H$\beta$ and [N II]/H$\alpha$ values at high redshift arises from sample selection: when we consider only the most massive galaxies $M_* \sim 10^{10-11} M_\odot$, the offset naturally appears, due to their high metallicities. We predict that deeper observations that probe lower-mass galaxies will reveal galaxies that lie on a locus comparable to $z\sim 0$ observations. Even when accounting for sample selection effects, we find that there is a subtle mismatch between simulations and observations. To resolve this discrepancy, we investigate the impact of varying ionization parameters, H II region densities, gas-phase abundance patterns, and increasing radiation field hardness on N2-BPT diagrams. We find that either decreasing the ionization parameter or increasing the N/O ratio of galaxies at fixed O/H can move galaxies along a self-similar arc in N2-BPT space that is occupied by high-redshift galaxies., Comment: 25 pages, 15 figures. Accepted for publication in ApJ

The combination of the MOSDEF and KBSS-MOSFIRE surveys represents the largest joint investment of Keck/MOSFIRE time to date, with ~3000 galaxies at 1.4<=z<=3.8, roughly half of which are at z~2. MOSDEF is photometric- and spectroscopic-redshift selected with a rest-optical magnitude limit, while KBSS-MOSFIRE is primarily selected based on rest-UV colors and a rest-UV magnitude limit. Analyzing both surveys in a uniform manner with consistent spectral-energy-distribution (SED) models, we find that the MOSDEF z~2 targeted sample has a higher median M_* and redder rest U-V color than the KBSS-MOSFIRE z~2 targeted sample, and a smaller median SED-based SFR and sSFR (SFR(SED) and sSFR(SED)). Specifically, MOSDEF targeted a larger population of red galaxies with U-V and V-J >=1.25, while KBSS-MOSFIRE contains more young galaxies with intense star formation. Despite these differences in the z~2 targeted samples, the subsets of the surveys with multiple emission lines detected and analyzed in previously published work are much more similar. All median host-galaxy properties with the exception of stellar population age -- i.e., M_*, SFR(SED), sSFR(SED), A_V, and UVJ colors -- agree within the uncertainties. Additionally, when uniform emission-line fitting and stellar Balmer absorption correction techniques are applied, there is no significant offset between the two samples in the [OIII]$\lambda$5008/H$\beta$ vs. [NII]$\lambda$6585/H$\alpha$ diagnostic diagram, in contrast to previously-reported discrepancies. We can now combine the MOSDEF and KBSS-MOSFIRE surveys to form the largest z~2 sample with moderate-resolution rest-optical spectra and construct the fundamental scaling relations of star-forming galaxies during this important epoch., Comment: 23 pages, 14 figures, published in MNRAS

In the next decade, deep galaxy surveys from telescopes such as the James Webb Space Telescope and Roman Space Telescope will provide transformational data sets that will greatly enhance the understanding of galaxy formation during the epoch of reionization (EoR). In this work, we present the Deep Realistic Extragalactic Model (DREaM) for creating synthetic galaxy catalogs. Our model combines dark matter simulations, subhalo abundance matching and empirical models, and includes galaxy positions, morphologies, and spectral energy distributions (SEDs). The resulting synthetic catalog extends to redshifts $z \sim 12 $, and galaxy masses $\log_{10}(M/M_{\odot}) = 5 $ covering an area of $1\, {\rm deg}^2$ on the sky. We use DREaM to explore the science returns of a $1\, {\rm deg}^2$ \emph{Roman} UDF, and to provide a resource for optimizing ultra-deep survey designs. We find that a \emph{Roman} UDF to $\sim 30\,m_{\rm AB}$ will potentially detect more than $10^6$ $M_{\rm UV}<-17$ galaxies, with more than $10^4$ at redshifts $z>7$, offering an unparalleled dataset for constraining galaxy properties during the EoR. Our synthetic catalogs and simulated images are made publicly available to provide the community with a tool to prepare for upcoming data., Comment: 39 pages, 27 figures. Submitted to ApJ. Comments welcome