Search

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 an analysis of JWST Ly$\alpha$ spectroscopy of $z\gtrsim6.5$ galaxies, using observations in the public archive covering galaxies in four independent fields (GOODS-N, GOODS-S, Abell 2744, EGS). We measure Ly$\alpha$ emission line properties for a sample of $210$ $z\simeq6.5-13$ galaxies, with redshifts confirmed independently of Ly$\alpha$ in all cases. We present $3$ new detections of Ly$\alpha$ emission in JWST spectra, including a large equivalent width (EW $=143\ \r{A}$) Ly$\alpha$ emitter with strong CIV emission (EW $=21\ \r{A}$) at $z=7.1$ in GOODS-N. We measure the redshift-dependent Ly$\alpha$ EW distribution across our sample. We find that strong Ly$\alpha$ emission (EW $>25\ \r{A}$) becomes increasingly rare at earlier epochs, suggesting that the transmission of Ly$\alpha$ photons decreases by $4\times$ between $z\simeq5$ and $z\simeq9$. We describe potential implications for the IGM neutral fraction. There is significant field to field variance in the Ly$\alpha$ emitter fraction. In contrast to the three other fields, the EGS shows no evidence for reduced transmission of Ly$\alpha$ photons at $z\simeq7-8$, suggesting a significantly ionized sightline may be present in the field. We use available NIRCam grism observations from the FRESCO survey to characterize overdensities on large scales around known Ly$\alpha$ emitters in the GOODS fields. The strongest overdensities appear linked with extremely strong Ly$\alpha$ detections (EW $>50\ \r{A}$) in most cases. Future Ly$\alpha$ spectroscopy with JWST has the potential to constrain the size of ionized regions around early galaxy overdensities, providing a new probe of the reionization process., Comment: 32 pages, 11 figures, submitted to ApJ

JWST has recently discovered a subset of reionization era galaxies with ionized gas that is metal poor in oxygen and carbon but heavily-enriched in nitrogen. This abundance pattern is almost never seen in lower redshift galaxies but is commonly observed in globular cluster stars. We have recently demonstrated that this peculiar abundance pattern appears in a compact ($\simeq 20$ pc) metal-poor galaxy undergoing a strong burst of star formation. This galaxy was originally selected based on strong CIV emission, indicating a hard radiation field rarely seen locally. In this paper, we present JWST/NIRSpec observations of another reionization-era galaxy known to power strong CIV emission, the $z=7.04$ gravitationally-lensed galaxy A1703-zd6. The emission line spectrum reveals this is a metal poor galaxy ($12+\log(\rm O/H) = 7.47\pm0.19$) dominated by a young stellar population ($1.6^{+0.5}_{-0.4}$ Myr) that powers a very hard ionizing spectrum (CIV EW = 19.4 $\unicode{x212B}$, He II EW = 2.2 $\unicode{x212B}$). The ISM is highly-enriched in nitrogen ($\log(\rm N/O)=-0.6$) with very high electron densities ($8-19\times10^4$ cm$^{-3}$) and extreme ionization conditions rarely seen at lower redshift. We also find intense CIV emission (EW$\gtrsim20$ $\unicode{x212B}$) in two new $z\gtrsim 6$ metal poor galaxies. To put these results in context, we search for UV line emission in a sample of 737 $z\gtrsim 4$ galaxies with NIRSpec spectra, establishing that 40(30)% of systems with [OIII]+H$\beta$ EW $>2000\unicode{x212B}$ have NIV] (CIV) detections with EW$>5$ $\unicode{x212B}$ ($>10$ $\unicode{x212B}$). These results suggest high N/O ratios and hard ionizing sources appear in a brief phase following a burst of star formation in compact high density stellar complexes., Comment: 29 pages, 18 figures; 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 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 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, 9 figures, 5 tables, submitted to 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: 21 pages, 7 figures, submitted to ApJ

JWST observations have recently begun delivering the first samples of Ly$\alpha$ velocity profile measurements at $z>6$, opening a new window on the reionization process. Interpretation of $z\gtrsim6$ line profiles is currently stunted by limitations in our knowledge of the intrinsic Ly$\alpha$ profile (before encountering the IGM) of the galaxies that are common at $z\gtrsim6$. To overcome this shortcoming, we have obtained resolved ($R\sim3900$) Ly$\alpha$ spectroscopy of $42$ galaxies at $z=2.1-3.4$ with similar properties as are seen at $z>6$. We quantify a variety of Ly$\alpha$ profile statistics as a function of [OIII]+H$\beta$ EW. Our spectra reveal a new population of $z\simeq 2-3$ galaxies with large [OIII]+H$\beta$ EWs ($>1200\ \r{A}$) and a large fraction of Ly$\alpha$ flux emerging near the systemic redshift (peak velocity $\simeq0$ km s$^{-1}$). These spectra indicate that low density neutral hydrogen channels are able to form in a subset of low mass galaxies ($\lesssim1\times10^8\ M_{\odot}$) that experience a burst of star formation (sSFR $>100$ Gyr$^{-1}$). Other extreme [OIII] emitters show weaker Ly$\alpha$ that is shifted to higher velocities ($\simeq240$ km s$^{-1}$) with little emission near line center. We investigate the impact the IGM is likely to have on these intrinsic line profiles in the reionization era, finding that the centrally peaked Ly$\alpha$ emitters should be strongly attenuated at $z\gtrsim5$. We show that these line profiles are particularly sensitive to the impact of resonant scattering from infalling IGM and can be strongly attenuated even when the IGM is highly ionized at $z\simeq 5$. We compare these expectations against a new database of $z\gtrsim6.5$ galaxies with robust velocity profiles measured with JWST/NIRSpec., Comment: 28 pages, 11 figures, accepted by ApJ

Nearly a decade ago, we began to see indications that reionization-era galaxies power hard radiation fields rarely seen at lower redshift. Most striking were detections of nebular CIV emission in what appeared to be typical low mass galaxies, requiring an ample supply of 48 eV photons to triply ionize carbon. The nature of this population has long remained unclear owing to limitations of ground-based spectroscopy. We have obtained deep JWST/NIRSpec R=1000 spectroscopy of the two z>6 CIV-emitting galaxies known prior to JWST. Here we present a rest-UV to optical spectrum of one of these two systems, the multiply-imaged z=6.1 lensed galaxy RXCJ2248-ID. NIRCam imaging reveals two compact (<22pc) clumps separated by 220pc, with one comprising a dense concentration of massive stars ($>10,400M_{\odot}$/yr/kpc$^2$) formed in a recent burst. We stack spectra of 3 images of the galaxy (J=24.8-25.9), yielding a very deep spectrum providing a high S/N template of strong emission line sources at z>6. The spectrum reveals narrow high ionization lines (HeII, CIV, NIV]) with line ratios consistent with powering by massive stars. The rest-optical spectrum is dominated by very strong emission lines ([OIII] EW=2800\AA), albeit with weak emission from low-ionization transitions ([OIII]/[OII]=184). The electron density is found to be very high($6.4-31\times10^4$cm$^{-3}$) based on three UV transitions. The ionized gas is metal poor ($12+\log(\rm O/H)=7.43^{+0.17}_{-0.09}$), yet highly enriched in nitrogen ($\log(\rm N/O)=-0.39^{+0.11}_{-0.10}$). The spectrum appears broadly similar to that of GNz11 at z=10.6, without showing the same AGN signatures. We suggest that the hard radiation field and rapid nitrogen enrichment may be a short-lived phase that many z>6 galaxies go through as they undergo strong bursts of star formation. We comment on the potential link of such spectra to globular cluster formation., Comment: 22 pages, 13 figures, 4 tables; Accepted to MNRAS

We describe {\it JWST}/NIRSpec prism measurements of Ly$\alpha$ emission in $z\gtrsim 5$ galaxies. We identify Ly$\alpha$ detections in 10 out of 69 galaxies with robust rest-optical emission line redshift measurements at $5\leq z<7$ in the CEERS and DDT-2750 observations of the EGS field. Galaxies at $z\simeq 6$ with faint continuum (F150W $=$ 27--29 mag) are found with extremely large rest-frame Ly$\alpha$ equivalent widths (ranging up to 286 A). Likely Ly$\alpha$ detections are also seen in two new $z>7$ galaxies ($z=$ 7.49 and 7.17) from the second epoch of CEERS observations, both showing large Ly$\alpha$ equivalent widths that likely indicate significant transmission through the IGM. We measure high Ly$\alpha$ escape fractions in the 12 Ly$\alpha$ emitters in our sample (median 0.28), two of which show $f_{\rm esc}^{ {\rm Ly}\alpha}$ near unity ($>0.80$). We find that $50_{-11}^{+11}$% of $z\simeq 6$ galaxies with [OIII]+H$\beta$ EW $>$ 1000 A have $f_{\rm esc}^{ {\rm Ly}\alpha}$ $>0.2$, consistent with the fractions found in lower-redshift samples with matched [OIII]+H$\beta$ EWs. While uncertainties are still significant, we find that only $10_{-5}^{+9}$% of $z>7$ galaxies with similarly strong rest optical emission lines show such large $f_{\rm esc}^{ {\rm Ly}\alpha}$, as may be expected if IGM attenuation of Ly$\alpha$ increases towards higher redshifts. We identify photometric galaxy overdensities near the $z\gtrsim 7$ Ly$\alpha$ emitters, potentially providing the ionizing flux necessary to create large ionized sightlines that facilitate Ly$\alpha$ transmission. Finally, we investigate the absence of Ly$\alpha$ emission in a comparable (and spectroscopically confirmed) galaxy overdensity at $z=7.88$ in the Abell 2744 field, discussing new prism spectra of the field obtained with the UNCOVER program., Comment: main text 20 pages, 13 figures; accepted for publication by MNRAS

We summarize the properties and initial data release of the JADES Origins Field (JOF), which will soon be the deepest imaging field yet observed with the James Webb Space Telescope (JWST). This field falls within the GOODS-S region about 8' south-west of the Hubble Ultra Deep Field (HUDF), where it was formed initially in Cycle 1 as a parallel field of HUDF spectroscopic observations within the JWST Advanced Deep Extragalactic Survey (JADES). This imaging will be greatly extended in Cycle 2 program 3215, which will observe the JOF for 5 days in six medium-band filters, seeking robust candidates for z>15 galaxies. This program will also include ultra-deep parallel NIRSpec spectroscopy (up to 104 hours on-source, summing over the dispersion modes) on the HUDF. Cycle 3 observations from program 4540 will add 20 hours of NIRCam slitless spectroscopy to the JOF. With these three campaigns, the JOF will be observed for 380 open-shutter hours with NIRCam using 15 imaging filters and 2 grism bandpasses. Further, parts of the JOF have deep 43 hr MIRI observations in F770W. Taken together, the JOF will soon be one of the most compelling deep fields available with JWST and a powerful window into the early Universe. This paper presents the second data release from JADES, featuring the imaging and catalogs from the year 1 JOF observations., Comment: Submitted to ApJ Supplement. Images and catalogs are available at https://archive.stsci.edu/hlsp/jades . A FITSmap portal to view the images is at https://jades.idies.jhu.edu

By combining the JWST/NIRCam JADES and CEERS extragalactic datasets, we have uncovered a sample of twenty-one T and Y brown dwarf candidates at best-fit distances between 0.1 - 4.2 kpc. These sources were selected by targeting the blue 1$\mu$m - 2.5$\mu$m colors and red 3$\mu$m - 4.5$\mu$m colors that arise from molecular absorption in the atmospheres of T$_{\mathrm{eff}} < $ 1300K brown dwarfs. We fit these sources using multiple models of low-mass stellar atmospheres and present the resulting fluxes, sizes, effective temperatures and other derived properties for the sample. If confirmed, these fits place the majority of the sources in the Milky Way thick disk and halo. We observe proper motion for seven of the candidate brown dwarfs with directions in agreement with the plane of our galaxy, providing evidence that they are not extragalactic in nature. We demonstrate how the colors of these sources differ from selected high-redshift galaxies, and explore the selection of these sources in planned large-area JWST NIRCam surveys. Deep imaging with JWST/NIRCam presents an an excellent opportunity for finding and understanding these very cold low-mass stars at kpc distances., Comment: 22 pages, 9 figures, accepted by ApJ (January 18, 2024)

The power-law slope of the rest-UV continuum ($f_{\lambda}\propto\lambda^{\beta}$) is a key metric of early star forming galaxies, providing one of our only windows into the stellar populations and physical conditions of $z>10$ galaxies. Expanding upon previous studies with limited sample sizes, we leverage deep imaging from JADES to investigate the UV slopes of 179 $z>9$ galaxies with apparent magnitudes of $m_{\rm F200W}=26-31$, which display a median UV slope of $\beta=-2.4$. We compare to a statistical sample of $z=5-9$ galaxies, finding a shift toward bluer rest-UV colors at all $\rm~M_{UV}$. The most UV-luminous $z>9$ galaxies are significantly bluer than their lower-redshift counterparts, representing a dearth of moderately-red galaxies in the first $500~$Myr. At yet earlier times, the $z>11$ galaxy population exhibits very blue UV slopes, implying very low attenuation from dust. We identify a robust sample of 44 galaxies with $\beta<-2.8$, which have SEDs requiring models of density-bounded HII regions and median ionizing photon escape fractions of $0.51$ to reproduce. Their rest-optical colors imply that this sample has weaker emission lines (median $m_{\rm F356W}-m_{\rm F444W}=0.19$ mag) than typical galaxies (median $m_{\rm F356W}-m_{\rm F444W}=0.39$ mag), consistent with the inferred escape fractions. This sample has relatively low stellar masses (median $\log(M/M_{\odot})=7.5$), and specific star-formation rates (median$=79\rm/Gyr$) nearly twice that of our full sample (median$=44\rm/Gyr$), suggesting they are more common among systems experiencing a recent upturn in star formation. We demonstrate that the shutoff of star formation provides an alternative solution for modelling of extremely blue UV colors, making distinct predictions for the rest-optical emission of these galaxies. Future spectroscopy will be required to distinguish between these physical pictures., Comment: 17 pages, 13 figures; submitted to MNRAS

Reionization is thought to be driven by faint star-forming galaxies, but characterizing this population has long remained very challenging. Here we utilize deep nine-band NIRCam imaging from JADES to study the star-forming and ionizing properties of 756 $z\sim6-9$ galaxies, including hundreds of very UV-faint objects ($M_\mathrm{UV}>-18$). The faintest ($m\sim30$) galaxies in our sample typically have stellar masses of $M_\ast\sim(1-3)\times10^7$ $M_\odot$ and young light-weighted ages ($\sim$50 Myr), though some show strong Balmer breaks implying much older ages ($\sim$500 Myr). We find no evidence for extremely massive galaxies ($>3\times10^{10}$ $M_\odot$). We infer a strong (factor $>$2) decline in the typical [OIII]$+$H$\beta$ EWs towards very faint $z\sim6-9$ galaxies, yet a weak UV luminosity dependence on the H$\alpha$ EWs at $z\sim6$. We demonstrate that these EW trends can be explained if fainter galaxies have systematically lower metallicities as well as more recently-declining star formation histories relative to the most UV-luminous galaxies in our sample. Our data provide evidence that the brightest galaxies are frequently experiencing a recent strong upturn in SFR. We also discuss how the EW trends may be influenced by a strong correlation between $M_\mathrm{UV}$ and Lyman continuum escape fraction. This alternative explanation has dramatically different implications for the contribution of galaxies along the luminosity function to cosmic reionization. Finally, we quantify the photometric overdensities around two $z>7$ strong Ly$\alpha$ emitters. One Ly$\alpha$ emitter lies close to a strong photometric overdensity while the other shows no significant nearby overdensity, perhaps implying that not all strong $z>7$ Ly$\alpha$ emitters reside in large ionized bubbles., Comment: Accepted for publication in MNRAS. Minor changes relative to previous version. Catalog of galaxy coordinates, inferred properties, and photometry available here: https://www.ryan-endsley.com/z6to9-jades-galaxy-catalog

We present a catalog of 717 candidate galaxies at $z > 8$ selected from 125 square arcminutes of NIRCam imaging as part of the JWST Advanced Deep Extragalactic Survey (JADES). We combine the full JADES imaging dataset with data from the JEMS and FRESCO JWST surveys along with extremely deep existing observations from HST/ACS for a final filter set that includes fifteen JWST/NIRCam filters and five HST/ACS filters. The high-redshift galaxy candidates were selected from their estimated photometric redshifts calculated using a template fitting approach, followed by visual inspection from seven independent reviewers. We explore these candidates in detail, highlighting interesting resolved or extended sources, sources with very red long-wavelength slopes, and our highest redshift candidates, which extend to $z_{phot} = 18$. We also investigate potential contamination by stellar objects, and do not find strong evidence from SED fitting that these faint high-redshift galaxy candidates are low-mass stars. Over 93\% of the sources are newly identified from our deep JADES imaging, including 31 new galaxy candidates at $z_{phot} > 12$. Using 42 sources in our sample with measured spectroscopic redshifts from NIRSpec and FRESCO, we find excellent agreement to our photometric redshift estimates, with no catastrophic outliers and an average difference of $\langle \Delta z = z_{phot}- z_{spec} \rangle= 0.26$. These sources comprise one of the most robust samples for probing the early buildup of galaxies within the first few hundred million years of the Universe's history., Comment: v2: 43 pages, 20 figures, accepted by The Astrophysical Journal, full online data catalog found at https://doi.org/10.5281/zenodo.7996499

We describe the NIRSpec component of the JWST Deep Extragalactic Survey (JADES), and provide deep spectroscopy of 253 sources targeted with the NIRSpec micro-shutter assembly in the Hubble Ultra Deep Field and surrounding GOODS-South. The multi-object spectra presented here are the deepest so far obtained with JWST, amounting to up to 28 hours in the low-dispersion ($R\sim 30-300$) prism, and up to 7 hours in each of the three medium-resolution $R\approx 1000$ gratings and one high-dispersion grating, G395H ($R\approx2700$). Our low-dispersion and medium-dispersion spectra cover the wavelength range $0.6-5.3\mu$m. We describe the selection of the spectroscopic targets, the strategy for the allocation of targets to micro-shutters, and the design of the observations. We present the public release of the reduced 2D and 1D spectra, and a description of the reduction and calibration process. We measure spectroscopic redshifts for 178 of the objects targeted extending up to $z=13.2$. We present a catalog of all emission lines detected at $S/N>5$, and our redshift determinations for the targets. Combined with the first JADES NIRCam data release, these public JADES spectroscopic and imaging datasets provide a new foundation for discoveries of the infrared universe by the worldwide scientific community., Comment: Accepted for publication in A&A. Data products available from https://archive.stsci.edu/hlsp/jades

We present an overview of the James Webb Space Telescope (JWST) Advanced Deep Extragalactic Survey (JADES), an ambitious program of infrared imaging and spectroscopy in the GOODS-S and GOODS-N deep fields, designed to study galaxy evolution from high redshift to cosmic noon. JADES uses about 770 hours of Cycle 1 guaranteed time largely from the Near-Infrared Camera (NIRCam) and Near-Infrared Spectrograph (NIRSpec) instrument teams. In GOODS-S, in and around the Hubble Ultra Deep Field and Chandra Deep Field South, JADES produces a deep imaging region of ~45 arcmin$^2$ with an average of 130 hrs of exposure time spread over 9 NIRCam filters. This is extended at medium depth in GOODS-S and GOODS-N with NIRCam imaging of ~175 arcmin$^2$ with an average exposure time of 20 hrs spread over 8-10 filters. In both fields, we conduct extensive NIRSpec multi-object spectroscopy, including 2 deep pointings of 55 hrs exposure time, 14 medium pointings of ~12 hrs, and 15 shallower pointings of ~4 hrs, targeting over 5000 HST and JWST-detected faint sources with 5 low, medium, and high-resolution dispersers covering 0.6-5.3 microns. Finally, JADES extends redward via coordinated parallels with the JWST Mid-Infrared Instrument (MIRI), featuring ~9 arcmin$^2$ with 43 hours of exposure at 7.7 microns and twice that area with 2-6.5 hours of exposure at 12.8 microns For nearly 30 years, the GOODS-S and GOODS-N fields have been developed as the premier deep fields on the sky; JADES is now providing a compelling start on the JWST legacy in these fields., Comment: 33 pages, submitted to ApJ Supplement. The JADES Collaboration web site is at https://jades-survey.github.io, and the initial data release is available at https://archive.stsci.edu/hlsp/jades with a viewer at http://jades.idies.jhu.edu

Several studies have detected Lyman-alpha (Ly$\alpha$) from bright ($M_\mathrm{UV}\lesssim-21.5$) galaxies during the early stages of reionization despite the significantly neutral intergalactic medium. To explain these detections, it has been suggested that $z>7$ Ly$\alpha$ emitters (LAEs) inhabit physical Mpc (pMpc)-scale ionized regions powered by overdensities of faint galaxies, but systematic searches for these overdensities near LAEs have been challenging. Here, we use CEERS JWST/NIRCam imaging to search for large-scale galaxy overdensities near two very UV-bright, $z=8.7$ LAEs in the EGS field. We colour select 27 $z=8.4-9.1$ candidates, including the one LAE in the footprint (EGSY8p7). From SED models, we infer moderately faint UV luminosities ($-21.2\lesssim{M_\mathrm{UV}}\lesssim -19.1$) and stellar masses of $M_*\approx10^{7.5-8.8}$ M$_\odot$. All are efficient ionizing agents ($\xi_{\mathrm{ion}}^{*}\approx10^{25.5-26.0}$ Hz erg$^{-1}$) and are generally morphologically simple with only one compact ($r_e\lesssim140$ to $\sim650$ pc) star-forming component. 13 candidates lie within 5 arcmin of EGSY8p7, leading to a factor-of-four galaxy overdensity at $\lesssim 5$ arcmin ($\sim 1.4$ projected pMpc at $z\sim8.7$) separations from EGSY8p7. Separations of $10-15$ arcmin ($\sim2.7-4.1$ projected pMpc) are consistent with an average field. The spatial distribution of our sample may qualitatively suggest an $R\geq2$ pMpc ionized bubble encompassing both LAEs in EGS, which is theoretically unexpected but may be possible for a galaxy population $4\times$ more numerous than the average to create with moderate escape fractions ($f_\mathrm{esc}\gtrsim0.15$) over long times ($\gtrsim200$ Myr). Upcoming spectroscopic follow-up will characterize the size of any ionized bubble that may exist and the properties of the galaxies powering such a bubble., Comment: 20 pages, 7 main figures, 3 appendices, accepted to MNRAS

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

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 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 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 present JADES JWST/NIRSpec spectroscopy of GN-z11, the most luminous candidate $z>10$ Lyman break galaxy in the GOODS-North field with $M_{UV}=-21.5$. We derive a redshift of $z=10.603$ (lower than previous determinations) based on multiple emission lines in our low and medium resolution spectra over $0.8-5.3 \mu$m. We significantly detect the continuum and measure a blue rest-UV spectral slope of $\beta=-2.4$. Remarkably, we see spatially-extended Lyman-$\alpha$ in emission (despite the highly-neutral IGM expected at this early epoch), offset 555 km s$^{-1}$ redward of the systemic redshift. From our measurements of collisionally-excited lines of both low- and high-ionization (including [O II]$\lambda3727$, [Ne III]$\lambda 3869$ and C III]$\lambda1909$) we infer a high ionization parameter ($\log U\sim -2$). We detect the rarely-seen N IV]$\lambda1486$ and N III]$\lambda1748$ lines in both our low and medium resolution spectra, with other high ionization lines seen in the low resolution spectrum such as He II (blended with O III]) and C IV (with a possible P-Cygni profile). Based on the observed rest-UV line ratios, we cannot conclusively rule out photoionization from AGN, although the high C III]/He II and N III]/He II ratios are compatible with a star-formation explanation. If the observed emission lines are powered by star formation, then the strong N III]$\lambda1748$ observed may imply an unusually high $N/O$ abundance. Balmer emission lines (H$\gamma$, H$\delta$) are also detected, and if powered by star formation rather than an AGN we infer a star formation rate of $\sim 20-30 M_{\odot} yr^{-1}$ (depending on the IMF) and low dust attenuation. Our NIRSpec spectroscopy confirms that GN-z11 is a remarkable galaxy with extreme properties seen 430 Myr after the Big Bang., Comment: Accepted for publication in Astronomy & Astrophysics, 17 pages, 12 figures. Line fluxes updated

We present JWST NIRCam 9-band near-infrared imaging of the luminous $z=10.6$ galaxy GN-z11 from the JWST Advanced Deep Extragalactic Survey (JADES) of the GOODS-N field. We find a spectral energy distribution (SED) entirely consistent with the expected form of a high-redshift galaxy: a clear blue continuum from 1.5 to 4 microns with a complete dropout in F115W. The core of GN-z11 is extremely compact in JWST imaging. We analyze the image with a two-component model, using a point source and a S\'{e}rsic profile that fits to a half-light radius of 200 pc and an index $n=0.9$. We find a low-surface brightness haze about $0.4''$ to the northeast of the galaxy, which is most likely a foreground object but might be a more extended component of GN-z11. At a spectroscopic redshift of 10.60 (Bunker et al. 2023), the comparison of the NIRCam F410M and F444W images spans the Balmer jump. From population synthesis modeling, here assuming no light from an active galactic nucleus, we reproduce the SED of GN-z11, finding a stellar mass of $\sim$$10^{9}~M_{\odot}$, a star-formation rate of $\sim$$20~M_{\odot}~\mathrm{yr}^{-1}$ and a young stellar age of $\sim$$20~\mathrm{Myr}$. As massive galaxies at high redshift are likely to be highly clustered, we search for faint neighbors of GN-z11, finding 9 galaxies out to $\sim$5 comoving Mpc transverse with photometric redshifts consistent with $z=10.6$, and a 10$^{\rm th}$ more tentative dropout only $3''$ away. This is consistent with GN-z11 being hosted by a massive dark-matter halo ($\approx8\times10^{10}~M_{\odot}$), though lower halo masses cannot be ruled out., Comment: 18 pages, 8 figures, accepted for publication in 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

Finding and characterising the first galaxies that illuminated the early Universe at cosmic dawn is pivotal to understand the physical conditions and the processes that led to the formation of the first stars. In the first few months of operations, imaging from the James Webb Space Telescope (JWST) have been used to identify tens of candidates of galaxies at redshift (z) greater than 10, less than 450 million years after the Big Bang. However, none of these candidates has yet been confirmed spectroscopically, leaving open the possibility that they are actually low-redshift interlopers. Here we present spectroscopic confirmation and analysis of four galaxies unambiguously detected at redshift 10.3

Spitzer/IRAC imaging has revealed that the brightest $z\sim7-8$ galaxies often exhibit young ages and strong nebular line emission, hinting at high ionizing efficiency among early galaxies. However, IRAC's limited sensitivity has long hindered efforts to study the fainter, more numerous population often thought largely responsible for reionization. Here we use CEERS JWST/NIRCam data to characterize 116 UV-faint (median M$_{UV}=-19.5$) $z\sim6.5-8$ galaxies. The SEDs are typically dominated by young ($\sim$10-50 Myr), low-mass ($M_\ast\sim10^8\ M_\odot$) stellar populations, and we find no need for extremely high stellar masses ($\sim10^{11} M_\odot$). Considering previous studies of UV-bright (M$_{UV}\sim-22$) $z\sim7-8$ galaxies, we find evidence for a strong (5-10$\times$) increase in specific star formation rate toward lower luminosities (median sSFR=103 Gyr$^{-1}$ in CEERS). The larger sSFRs imply a more dominant contribution from OB stars in the relatively numerous UV-faint population, perhaps suggesting that these galaxies are very efficient ionizing agents (median $\xi_{ion}=10^{25.7}$ erg$^{-1}$ Hz). In spite of their much larger sSFRs, we find no significant increase in [OIII]$+$H$\beta$ EWs towards fainter M$_{UV}$ (median $\approx$780 $\mathring{A}$). If confirmed, this may indicate that a substantial fraction of our CEERS galaxies possess extremely low metallicities ($\lesssim$3% $Z_\odot$) where [OIII] emission is suppressed. Alternatively, high ionizing photon escape fractions or bursty star formation histories can also weaken the nebular lines in a subset of our CEERS galaxies. While the majority of our objects are very blue (median $\beta=-2.0$), we identify a significant tail of very dusty galaxies ($\beta\sim-1$) at $\approx$0.5$L_{UV}^\ast$ which may contribute significantly to the $z\sim7-8$ star formation rate density., Comment: Accepted in MNRAS. Updated to use the most recent NIRCam zeropoints. There are no significant changes to the conclusions relative to v1

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

The ultraviolet (UV) continuum slope ($\beta~$where$~f_\lambda\propto\lambda^\beta$) of galaxies is sensitive to a variety of properties, from the metallicity and age of the stellar population to dust attenuation throughout the galaxy. Considerable attention has focused on identifying reionization-era galaxies with very blue UV slopes ($\beta<-3$). Not only do such systems provide a signpost of low-metallicity stars, but they also identify galaxies likely to leak ionizing photons from their HII regions as such blue UV slopes require the reddening effect of nebular continuum to be diminished. In this paper we present a search for reionization-era galaxies with very blue UV colors in recent JWST/NIRCam imaging of the EGS field. We characterize UV slopes for a large sample of$~z\simeq7-11~$galaxies, finding a median of$~\beta=-2.0$. Two lower luminosity (M$_{\rm{UV}}\simeq-19.5$) and lower stellar mass (6-10$\times10^7$M$_\odot$) systems exhibit extremely blue UV slopes ($\beta=-2.9~$to$~-3.1$) and rest-optical photometry indicating weak nebular line emission. Each system is very compact (r$_e\lesssim$260pc) with very high star-formation-rate surface densities. We model the SEDs with a suite of BEAGLE models with varying levels of ionizing photon escape. The SEDs cannot be reproduced with our fiducial (f$_{\rm{esc,HII}}$=0) or alpha-enhanced (Z$_\star7$., Comment: 13 pages, 10 figures, 3 tables; Accepted to ApJ

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 present specific star-formation rates for 40 UV-bright galaxies at $z\sim7-8$ observed as part of the Reionization Era Bright Emission Line Survey (REBELS) ALMA large program. The sSFRs are derived using improved measures of SFR and stellar masses, made possible by measurements of far-infrared (FIR) continuum emission and [CII]-based spectroscopic redshifts. For each source in the sample, we derive stellar masses from SED fitting and total SFRs from calibrations of the UV and FIR emission. The median sSFR is $18_{-5}^{+7}$ Gyr$^{-1}$, significantly larger than literature measurements lacking constraints in the FIR. The increase in sSFR reflects the larger obscured SFRs we derive from the dust continuum relative to that implied by the UV+optical SED. We suggest that such differences may reflect spatial variations in dust across these luminous galaxies, with the component dominating the FIR distinct from that dominating the UV. We demonstrate that the inferred stellar masses (and hence sSFRs) are strongly-dependent on the assumed star formation history (SFH) in reionization-era galaxies. When large sSFR galaxies are modeled with non-parametric SFHs, the derived stellar masses can increase by an order of magnitude relative to constant star formation models, owing to the presence of a significant old stellar population that is outshined by the recent burst. The [CII] line widths in the largest sSFR systems are often very broad, suggesting dynamical masses that are easily able to accommodate the dominant old stellar population suggested by non-parametric models. Regardless of these systematic uncertainties in the derived parameters, we find that the sSFR increases rapidly toward higher redshifts for massive galaxies ($9.6<\log(\rm M_*/M_{\odot})<9.8$), with a power law that goes as $(1+z)^{1.7\pm0.3}$, broadly consistent with expectations from the evolving baryon accretion rates., Comment: 18 pages, 10 figures, submitted to MNRAS

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

We perform joint modeling of the composite rest-frame far-UV (FUV) and optical spectra of redshift 1.85

We present deep Keck/MOSFIRE $H$-band spectroscopic observations covering the [CIII],CIII]$\lambda\lambda1907,1909$ doublet for three $z\sim8$ galaxy candidates in the AEGIS field. Along with non-detections in two galaxies, we obtain one of the highest-redshift detections to-date of [CIII]$\lambda 1907$ for the galaxy AEGIS-33376, finding $z_{\rm spec}=7.945\pm0.001$. We measure a [CIII]$\lambda$1907 flux of $2.24\pm0.71\times10^{-18} \mbox{ erg}\mbox{ s}^{-1} \mbox{ cm}^{-2}$, corresponding to a rest-frame equivalent width of $20.3\pm6.5 \unicode{x212B}$ for the single line. Given the not very constraining upper limit for CIII]$\lambda 1909$ based on strong sky-line contamination, we assume a [CIII]$\lambda$1907/CIII]$\lambda 1909$ doublet ratio of 1.5 and infer a total [CIII],CIII]$\lambda\lambda1907,1909$ equivalent width of $33.7\pm 10.8 \unicode{x212B}$. We repeat the same reductions and analysis on multiple subsets of our raw data divided on the basis of time and observing conditions, verifying that the [CIII]$\lambda 1907$ emission is present for AEGIS-33376 throughout our observations. We also confirm that the significance of the [CIII]$\lambda 1907$ detection in different subsets of our data tracks that of brighter emission features detected on the same multi-slit mask. These multiple tests suggest that the observed emission line is real and associated with the $z\sim 8$ target. The strong observed [CIII],CIII]$\lambda\lambda1907,1909$ in AEGIS-33376 likely indicates ISM conditions of low metallicity, high ionization parameter, and a hard ionizing spectrum, although AGN contributions are possible. This single detection represents a sizable increase in the current sample [CIII],CIII]$\lambda\lambda1907,1909$ detections at $z>7$, while $\textit{JWST}$ will provide the first statistical samples of such measurements at these redshifts., Comment: 6 pages, 3 figures, 1 table; Accepted to ApJL

The complex structure of gas, metals, and dust in the interstellar and circumgalactic medium (ISM and CGM, respectively) in star-forming galaxies can be probed by Ly$\alpha$ emission and absorption, low-ionization interstellar (LIS) metal absorption, and dust reddening E(B-V). We present a statistical analysis of the mutual correlations among Ly$\alpha$ equivalent width (EW$_{Ly\alpha}$), LIS equivalent width (EW$_{LIS}$), and E(B-V) in a sample of 157 star-forming galaxies at $z\sim2.3$. With measurements obtained from individual, deep rest-UV spectra and spectral-energy distribution (SED) modeling, we find that the tightest correlation exists between EW$_{LIS}$ and E(B-V), although correlations among all three parameters are statistically significant. These results signal a direct connection between dust and metal-enriched HI gas, and that they are likely co-spatial. By comparing our results with the predictions of different ISM/CGM models, we favor a dusty ISM/CGM model where dust resides in HI gas clumps and Ly$\alpha$ photons escape through the low HI covering fraction/column density intra-clump medium. Finally, we investigate the factors that potentially contribute to the intrinsic scatter in the correlations studied in this work, including metallicity, outflow kinematics, Ly$\alpha$ production efficiency, and slit loss. Specifically, we find evidence that scatter in the relationship between EW$_{Ly\alpha}$ and E(B-V) reflects the variation in metal-to-HI covering fraction ratio as a function of metallicity, and the effects of outflows on the porosity of the ISM/CGM. Future simulations incorporating star-formation feedback and the radiative transfer of Ly$\alpha$ photons will provide key constraints on the spatial distributions of neutral hydrogen gas and dust in the ISM/CGM structure., Comment: 19 pages, 6 figures, 3 tables, accepted for publication in ApJ

We analyze the rest-optical emission-line ratios of z~1.5 galaxies drawn from the MOSFIRE Deep Evolution Field (MOSDEF) survey. Using composite spectra we investigate the mass-metallicity relation (MZR) at z~1.5 and measure its evolution to z=0. When using gas-phase metallicities based on the N2 line ratio, we find that the MZR evolution from z~1.5 to z=0 depends on stellar mass, evolving by $\Delta\rm log(\rm O/H)\sim0.25$ dex at $M_*<10^{9.75}M_{\odot}$ down to $\Delta\rm log(\rm O/H)\sim0.05$ at $M_*>10^{10.5}M_{\odot}$. In contrast, the O3N2-based MZR shows a constant offset of $\Delta\rm log(\rm O/H)\sim0.30$ across all masses, consistent with previous MOSDEF results based on independent metallicity indicators, and suggesting that O3N2 provides a more robust metallicity calibration for our z~1.5 sample. We investigated the secondary dependence of the MZR on SFR by measuring correlated scatter about the mean $M_*$-specific SFR and $M_*-\log(\rm O3N2)$ relations. We find an anti-correlation between $\log(\rm O/H)$ and sSFR offsets, indicating the presence of a $M_*$-SFR-Z relation, though with limited significance. Additionally, we find that our z~1.5 stacks lie along the z=0 metallicity sequence at fixed $\mu=\log(M_*/M_{\odot})-0.6\times\log(\rm SFR / M_{\odot} yr^{-1})$ suggesting that the z~1.5 stacks can be described by the z=0 fundamental metallicity relation (FMR). However, using different calibrations can shift the calculated metallicities off of the local FMR, indicating that appropriate calibrations are essential for understanding metallicity evolution with redshift. Finally, understanding how [NII]/H$\alpha$ scales with galaxy properties is crucial to accurately describe the effects of blended [NII] and H$\alpha$ on redshift and H$\alpha$ flux measurements in future large surveys utilizing low-resolution spectra such as with Euclid and the Roman Space Telescope., Comment: 14 Pages, 9 figures, Accepted to MNRAS

We present results on the properties of neon emission in $z\sim2$ star-forming galaxies drawn from the MOSFIRE Deep Evolution Field (MOSDEF) survey. Doubly-ionized neon ([NeIII]3869) is detected at $\geq3\sigma$ in 61 galaxies, representing $\sim$25% of the MOSDEF sample with H$\alpha$, H$\beta$, and [OIII]$5007$ detections at similar redshifts. We consider the neon emission-line properties of both individual galaxies with [NeIII]3869 detections and composite $z\sim2$ spectra binned by stellar mass. With no requirement of [NeIII]3869 detection, the latter provide a more representative picture of neon emission-line properties in the MOSDEF sample. The [NeIII]3869/[OII]3727 ratio (Ne3O2) is anti-correlated with stellar mass in $z\sim2$ galaxies, as expected based on the mass-metallicity relation. It is also positively correlated with the [OIII]$5007$/[OII]$3727$ ratio (O32), but $z\sim2$ line ratios are offset towards higher Ne3O2 at fixed O32, compared with both local star-forming galaxies and individual H~II regions. Despite the offset towards higher Ne3O2 at fixed O32 at $z\sim2$, biases in inferred Ne3O2-based metallicity are small. Accordingly, Ne3O2 may serve as an important metallicity indicator deep into the reionization epoch. Analyzing additional rest-optical line ratios including [NeIII]$3869$/[OIII]$5007$ (Ne3O3) and [OIII]$5007$/H$\beta$ (O3H$\beta$), we conclude that the nebular emission-line ratios of $z\sim2$ star-forming galaxies suggest a harder ionizing spectrum (lower stellar metallicity, i.e., Fe/H) at fixed gas-phase oxygen abundance, compared to systems at $z\sim0$. These new results based on neon lend support to the physical picture painted by oxygen, nitrogen, hydrogen, and sulfur emission, of an ionized ISM in high-redshift star-forming galaxies irradiated by chemically young, $\alpha$-enhanced massive stars., Comment: 7 pages, 5 figures, accepted to ApJ Letters

We analyze the rest-optical emission-line spectra of $z\sim2.3$ star-forming galaxies in the complete MOSFIRE Deep Evolution Field (MOSDEF) survey. In investigating the origin of the well-known offset between the sequences of high-redshift and local galaxies in the [O III]5008/H$\beta$ vs. [N II]6585/H$\alpha$ ("[N II] BPT") diagram, we define two populations of $z\sim2.3$ MOSDEF galaxies. These include the "high" population that is offset towards higher [O III]5008/H$\beta$ and/or [N II]6585/H$\alpha$ with respect to the local SDSS sequence and the "low" population that overlaps the SDSS sequence. These two groups are also segregated within the [O III]5008/H$\beta$ vs. [S II]6718,6733/H$\alpha$ and the [O III]4960,5008/[O II]3727,3730 (O$_{32}$) vs. ([O III]4960,5008+[O II]3727,3730)/H$\beta$ (R$_{23}$) diagram, which suggests qualitatively that star-forming regions in the more offset galaxies are characterized by harder ionizing spectra at fixed nebular oxygen abundance. We also investigate many galaxy properties of the split sample and find that the "high" sample is on average smaller in size and less massive, but has higher specific star-formation rate and star-formation-rate surface density values and is slightly younger compared to the "low" population. From Cloudy+BPASS photoionization models, we estimate that the "high" population has a lower stellar metallicity (i.e., harder ionizing spectrum) but slightly higher nebular metallicity and higher ionization parameter compared to the "low" population. While the "high" population is more $\alpha$-enhanced (i.e., higher $\alpha$/Fe) than the "low" population, both samples are significantly more $\alpha$-enhanced compared to local star-forming galaxies with similar rest-optical line ratios. These differences must be accounted for in all high-redshift star-forming galaxies -- not only those "offset" from local excitation sequences., Comment: 16 pages, 7 figures

We present constraints on the massive star and ionized gas properties for a sample of 62 star-forming galaxies at $z\sim2.3$. Using BPASS stellar population models, we fit the rest-UV spectra of galaxies in our sample to estimate age and stellar metallicity which, in turn, determine the ionizing spectrum. In addition to the median properties of well-defined subsets of our sample, we derive the ages and stellar metallicities for 30 high-SNR individual galaxies -- the largest sample of individual galaxies at high redshift with such measurements. Most galaxies in this high-SNR subsample have stellar metallicities of $0.001

We present a joint analysis of rest-UV and rest-optical spectra obtained using Keck/LRIS and Keck/MOSFIRE for a sample of 62 star-forming galaxies at $z\sim2.3$. We divide our sample into 2 bins based on their location in the [OIII]/Hb vs. [NII]/Ha BPT diagram, and perform the first differential study of the rest-UV properties of massive ionizing stars as a function of rest-optical emission-line ratios. Fitting BPASS stellar population synthesis models, including nebular continuum emission, to our rest-UV spectra, we find that high-redshift galaxies offset towards higher [OIII]/Hb and [NII]/Ha have younger ages (log(Age/yr)=$7.20^{+0.57}_{-0.20}$) and lower stellar metallicities ($Z_*=0.0010^{+0.0011}_{-0.0003}$) resulting in a harder ionizing spectrum, compared to the galaxies in our sample that lie on the local BPT star-forming sequence (log(Age/yr)=$8.57^{+0.88}_{-0.84}$, $Z_*=0.0019\pm0.0006$). Additionally, we find that the offset galaxies have an ionization parameter of log(U)=$-3.04^{+0.06}_{-0.11}$ and nebular metallicity of 12+log(O/H)=$8.40^{+0.06}_{-0.07}$, and the non-offset galaxies have log(U)=$-3.11\pm0.08$ and 12+log(O/H)=$8.30^{+0.05}_{-0.06}$. The stellar and nebular metallicities derived for our sample imply that the galaxies offset from the local BPT relation are more alpha-enhanced ($7.28^{+2.52}_{-2.82}$ O/Fe$_\odot$) compared to those consistent with the local sequence ($3.04^{+0.95}_{-0.54}$ O/Fe$_\odot$). However, even galaxies that are entirely consistent with the local nebular excitation sequence appear to be alpha-enhanced -- in contrast with typical local systems. Such differences must be considered when estimating gas-phase O/H at high redshift based on strong emission-line ratios. Specifically, a similarity in the location of high-redshift and local galaxies in the BPT diagram may not be indicative of a similarity in their physical properties., Comment: 16 pages, 12 figures. Accepted to MNRAS

Strong [OIII]$\lambda\lambda$4959,5007+H$\beta$ emission appears to be typical in star-forming galaxies at z > 6.5. As likely contributors to cosmic reionization, these galaxies and the physical conditions within them are of great interest. At z > 6.5, where Ly$\alpha$ is greatly attenuated by the intergalactic medium, rest-UV metal emission lines provide an alternative measure of redshift and also constraints on the physical properties of star-forming regions and massive stars. We present the first statistical sample of rest-UV line measurements in z $\sim$ 2 galaxies selected as analogs of those in the reionization era based on [OIII]$\lambda\lambda$4959,5007 EW or rest-frame U-B color. Our sample is drawn from the 3D-HST Survey and spans the redshift range 1.36 $\leqslant$ z $\leqslant$ 2.49. We find that the median Ly$\alpha$ and CIII]$\lambda\lambda$1907,1909 EWs of our sample are significantly greater than those of z $\sim$ 2 UV-continuum-selected star-forming galaxies. Measurements from both individual and composite spectra indicate a monotonic, positive correlation between CIII] and [OIII], while a lack of trend is observed between Ly$\alpha$ and [OIII] at [OIII] EW < 1000$\unicode{x212B}$. At higher [OIII] EW, extreme Ly$\alpha$ emission starts to emerge. Using stacked spectra, we find that Ly$\alpha$ and CIII] are significantly enhanced in galaxies with lower metallicity. Two objects in our sample appear comparable to z > 6.5 galaxies with exceptionally strong rest-UV metal line emission. These objects have significant CIV$\lambda\lambda$1548,1550, HeII$\lambda$1640, and OIII]$\lambda\lambda$1661,1665 emission in addition to intense Ly$\alpha$ or CIII]. Detailed characterization of these lower-redshift analogs provides unique insights into the physical conditions in z > 6.5 star-forming regions, motivating future observations of reionization-era analogs at lower redshifts., Comment: 22 pages, 9 figures, accepted for publication in ApJ

We perform a comprehensive analysis of the redshift evolution of the rest-UV spectra of star-forming galaxies out to z~5. We combine new z~5 measurements of HI Ly$\alpha$ and low- and high-ionization interstellar metal absorption features with comparable measurements at z~2-4. We measure the equivalent widths of interstellar absorption features using stacked spectra in bins of Ly$\alpha$ equivalent width, performing corrections to Ly$\alpha$ strengths based on a model for the transmission of the intergalactic medium. We find a strong correlation between decreasing low-ionization absorption strength and increasing Ly$\alpha$ emission strength over the redshift range z~2-5, suggesting that both of these quantities are fundamentally linked to neutral gas covering fraction. At the highest Ly$\alpha$ equivalent widths, we observe evolution at $z\sim5$ towards greater Ly$\alpha$ emission strength at fixed low-ionization absorption strength. If we interpret the non-evolving relationship of Ly$\alpha$ emission strength and low-ionization line strength at z~2-4 as primarily reflecting the radiative transfer of Ly$\alpha$ photons, this evolution at z~5 suggests a higher intrinsic production rate of Ly$\alpha$ photons than at lower redshift. Our conclusion is supported by the joint evolution of the relationships among Ly$\alpha$ emission strength, interstellar absorption strength, and dust reddening. We perform additional analysis in bins of stellar mass, star-formation rate, UV luminosity, and age, examining how the relationships between galaxy properties and Ly$\alpha$ emission evolve towards higher redshift. We conclude that increasing intrinsic Ly$\alpha$ photon production and strong detection of nebular CIV emission (signaling lower metallicity) at z~5 indicate an elevated ionized photon production efficiency ($\xi_{\rm ion}$)., Comment: 20 pages, 17 figures, accepted for publication in MNRAS