JWST has recently discovered a subset of reionization era galaxies with ionized gas that is metal poor in oxygen and carbon but heavilyenriched 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 metalpoor 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 JWSTNIRSpec observations of another reionizationera galaxy known to power strong CIV emission the z7.04 gravitationallylensed galaxy A1703zd6. The emission line spectrum reveals this is a metal poor galaxy 12logrm OH  7.47pm0.19 dominated by a young stellar population 1.60.50.4 Myr that powers a very hard ionizing spectrum CIV EW  19.4 unicodex212B He II EW  2.2 unicodex212B. The ISM is highlyenriched in nitrogen logrm NO0.6 with very high electron densities 819times104 cm3 and extreme ionization conditions rarely seen at lower redshift. We also find intense CIV emission EWgtrsim20 unicodex212B in two new zgtrsim 6 metal poor galaxies. To put these results in context we search for UV line emission in a sample of 737 zgtrsim 4 galaxies with NIRSpec spectra establishing that 4030 of systems with OIIIHbeta EW gt2000unicodex212B have NIV CIV detections with EWgt5 unicodex212B gt10 unicodex212B. These results suggest high NO ratios and hard ionizing sources appear in a brief phase following a burst of star formation in compact high density stellar complexes.
['Astrophysics - Astrophysics of Galaxies']
Deep restUV JWSTNIRSpec spectroscopy of early galaxies the demographics of CIV and Nemitters in the reionization era
{'Deep rest-UV JWST /NIRSpec spectroscopy of early galaxies: the demographics of CIV and N-emitters in the reionization era': "Michael W. Topping, 1 ★ Daniel P. Stark 1 , Peter Senchyna 2 , Zuyi Chen 1 , Adi Zitrin 3 , Ryan Endsley 4 , Stéphane Charlot 5 , Lukas J. Furtak 3 , Michael V. Maseda 6 , Adele Plat 1 , Renske Smit 7 , Ramesh Mainali 8 , Jacopo Chevallard 9 , Stephen Molyneux 7 , 10 , Jane R. Rigby 11  \n1 Department of Astronomy / Steward Observatory, University of Arizona, 933 N Cherry Ave, Tucson, AZ 85721  \n2 The Observatories of the Carnegie Institution for Science, 813 Santa Barbara Street, Pasadena, CA 91101, USA  \n- 3 Physics Department, Ben-Gurion University of the Negev, P.O. Box 653, Be'er-Sheva 84105, Israel  \n4 Department of Astronomy, University of Texas, Austin, TX 78712, USA  \n- 5 Sorbonne Université, CNRS, UMR 7095, Institut d'Astrophysique de Paris, 98 bis bd Arago, 75014 Paris, France\n- 6 Department of Astronomy, University of Wisconsin-Madison, 475 N. Charter St., Madison, WI 53706, USA\n- 7 Astrophysics Research Institute, Liverpool John Moores University, 146 Brownlow Hill, Liverpool L3 5RF, UK\n- 8 Observational Cosmology Lab, Code 665, NASA Goddard Space Flight Center, 8800 Greenbelt Rd., Greenbelt, MD 20771, USA\n- 9 Department of Physics, University of Oxford, Denys Wilkinson Building, Keble Road, Oxford OX1 3RH, UK\n- 10 European Southern Observatory, Karl-Schwarzschild-str. 2, 85748 Garching, Germany\n- 11 Astrophysics Science Division, Code 660, NASA Goddard Space Flight Center, 8800 Greenbelt Rd., Greenbelt, MD, 20771, USA.  \n30 July 2024", 'ABSTRACT': '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 ( ≃ 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 𝑧 = 7 . 04 gravitationally-lensed galaxy A1703-zd6. The emission line spectrum reveals this is a metal poor galaxy (12 + log ( O / H ) = 7 . 47 ± 0 . 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 Å, He II EW = 2.2 Å). The ISM is highly-enriched in nitrogen (log ( N / O ) = -0 . 6) with very high electron densities (8 -19 × 10 4 cm -3 ) and extreme ionization conditions rarely seen at lower redshift. We also find intense CIV emission (EW ≳ 20 Å) in two new 𝑧 ≳ 6 metal poor galaxies. To put these results in context, we search for UV line emission in a sample of 737 𝑧 ≳ 4 galaxies with NIRSpec spectra, establishing that 40(30)% of systems with [OIII]+H 𝛽 EW > 2000 Å have NIV] (CIV) detections with EW > 5 Å( > 10 Å). 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.  \nKey words: galaxies: evolution - galaxies: ISM - galaxies: high-redshift', '1 INTRODUCTION': "Roughly a decade ago, deep near-infrared spectroscopy began to provide our first window on the rest-frame ultraviolet (UV) spectra of 𝑧 ≳ 6 galaxies. This early glimpse revealed strong line emission from highly-ionized species of carbon and oxygen, with in-  \nties significantly in excess of those at lower redshifts (Stark et al. 2015b,a, 2017; Laporte et al. 2017; Mainali et al. 2017; Hutchison et al. 2019; Topping et al. 2021). This was most clearly seen in two gravitationally-lensed galaxies showing strong nebular CIV 𝜆𝜆 1548 , 1550 emission, RXCJ2248-ID at 𝑧 = 6 . 11 (Mainali et al. 2017; Schmidt et al. 2017), and A1703-zd6 at 𝑧 = 7 . 045 (Stark et al. 2015b). The presence of triply ionized carbon requires a supply of 48 eV photons rarely seen in star forming galaxies at  \nlower redshifts. Some argued this may be indicative of a population of very low metallicity stars (Stark et al. 2015b; Mainali et al. 2017), whereas others suggested the detections may point to the presence of narrow-line AGN in relatively faint reionization era galaxies (Nakajima et al. 2018). With ground-based instrumentation, our ability to delineate the nature of the ionizing sources in these systems was long stunted by limited wavelength coverage ( 𝜆 rest ≲ 2000 Å) imposed by our atmosphere.  \nThe launch of JWST has fundamentally altered our ability to characterize the spectra of reionization era galaxies (e.g., Schaerer et al. 2022a; Arellano-Córdova et al. 2022a; Katz et al. 2023; Trump et al. 2023; Curtis-Lake et al. 2022; Fujimoto et al. 2023a). Early observations have confirmed that hard radiation fields are present in a number of early galaxy spectra (e.g., Bunker et al. 2023; Tang et al. 2023; Castellano et al. 2024; Tacchella et al. 2023). Perhaps most striking has been the the 𝑧 = 10 . 60 galaxy GNz11 reported in Bunker et al. (2023). In addition to strong high ionization emission lines, the rest-UV spectrum reveals strong NIV] 𝜆𝜆 1483 , 1486 and NIII] 𝜆 1750 emission, indicating gas that is heavily enriched in nitrogen relative to oxygen (log ( N / O ) = -0 . 38) and carbon (log ( N / C ) ∼ -0 . 5) in spite of being very metal poor (12 + log ( O / H ) = 7 . 82; Cameron et al. 2023a). This abundance pattern is extremely rare in lower redshift galaxies, suggesting a chemical evolution pathway that is perhaps unique to the reionization era, conceivably requiring a population of very (or super) massive stars ( ≳ 10 2 -3 M ⊙ ; Vink 2023; Charbonnel et al. 2023; Marques-Chaves et al. 2024; Nandal et al. 2024b) or punctuated star formation histories (Kobayashi & Ferrara 2024). Similar abundance ratios are seen in many globular cluster stars (e.g., Gratton et al. 2004; Carretta et al. 2005; Senchyna et al. 2023), potentially indicating a connection between the stars being formed in GNz11 and those formed in the second generation stellar populations of globular clusters (e.g., Bastian & Lardo 2018). The GNz11 spectrum also shows evidence for very high electron densities ( ≳ 10 5 cm -3 ), two orders of magnitude above those seen in star forming galaxies at lower redshift (e.g., Sanders et al. 2016), potentially reflecting conditions associated with a dense broad line region or nuclear star cluster surrounding an AGN (Maiolino et al. 2023a; Senchyna et al. 2023).  \nWhether the peculiar abundance pattern of GNz11 is common in the reionization era is not known. It is also not clear whether the nitrogen-enhancement is somehow linked to the presence of the hard ionizing sources and high electron densities that are also seen in GNz11.Butwithoutalarger sample of such galaxies, it is impossible to adequately understand the physical mechanisms driving the origin of the abundance pattern. A natural place to begin investigation of these questions is the two 𝑧 ≳ 6 galaxies with CIV detections from ground-based telescopes, both of which are extremely bright in the continuum (J=24.8-25.9) allowing high S/N views of the ionizing sources and gas abundance pattern. While it was clear prior to JWST that both systems had a population of hard ionizing sources, it was not known whether they also had ionized gas that was nitrogen enhanced or elevated in density similar to GNz11. In JWST Cycle 1, we obtained NIRSpec R=1000 spectroscopy targeting both the 𝑧 = 6 . 11 (Mainali et al. 2017; Schmidt et al. 2017) and the 𝑧 = 7 . 05 (Stark et al. 2015a) CIV emitters (program ID: 2478). The observations were designed to extend the spectroscopic coverage of these sources into the rest-optical for the first time, providing insight into the metal content of early galaxies with hard radiation fields. We also obtained very deep rest-UV NIRSpec observations, providing complete coverage of the emission lines (and continuum) at ≃ 1200-2000 Å. This far UV spectral window contains a suite  \nof features sensitive to the origin of the radiation field, the electron density, and the level of nitrogen enhancement.  \nIn the first paper from this program, Topping et al. (2024b) presented a detailed investigation of the NIRSpec observations of RXCJ2248-ID, the 𝑧 = 6 . 11 CIV emitter. The spectrum confirmed the presence of a hard radiation field with strong emission from numerous highly-ionized species (i.e., CIV 𝜆𝜆 1548 , 1550, He II 𝜆 1640). The rest-optical spectrum was found to be similarly extreme, with very large [OIII]+H 𝛽 EW (3100 Å) only seen in the upper 2% of 𝑧 ≃ 6 -9 galaxies (e.g., Endsley et al. 2023a; Matthee et al. 2023). In spite of the strength of [OIII] and the hydrogen recombination lines, the [OII] doublet is found to be very weak, suggesting an [OIII]/[OII] flux ratio (hereafter O32) of 184, well in excess of nearly any star forming galaxy previously studied at lower redshifts. Topping et al. (2024b) report that the various spectral features are consistent with being powered by a young (1.8 Myr) stellar population in a ≲ 20 pc region with ionized gas that is metal poor (12 + log ( O / H ) = 7 . 43) and extremely dense (631 × 10 4 cm -3 ). The RXCJ2248-ID spectrum shows emission from NIII] and NIV], revealing a nitrogen-enhanced abundance pattern (log ( N / O ) = -0 . 39) that is similar to that found in GNz11. However in RXCJ2248-ID no spectral signatures were found that demand the presence of an AGN, suggesting that the high ionization lines may be likely to originate from a low metallicity population of massive stars.  \nBased on these results, Topping et al. (2024b) suggested that the nitrogen enhancement seen at 𝑧 ≳ 6 may be produced in the dense clusters that form during a strong burst of star formation, a phase that many reionization era galaxies go through (Furlanetto & Mirocha 2022; Mirocha & Furlanetto 2023; Dome et al. 2024; Endsley et al. 2023a; Strait et al. 2023). In this framework, the hard ionizing agents would be produced during the burst, and the high electron densities would reflect the extreme dense gas conditions associated with the burst. If this picture applies universally, we should also see nitrogen enhancements and high electron densities in the 𝑧 = 7 . 05 galaxy Abell 1703-zd6 (hereafter A1703-zd6), the other CIV emitting primary target in our Cycle 1 program. In this paper, we present the deep NIRSpec observations of A1703-zd6. Our MSA design includes many additional bright lensed 𝑧 ≳ 6 galaxies, which we also describe in this paper. Our primary goals are twofold. First we seek to understand more about what drives hard radiation fields in a subset of early galaxies. Second, we seek to understand if the hard radiation fields are linked in some way to spectra like GNz11, with nitrogen enhancements and high electron density. To put these results in context, we will additionally present constraints on the fraction of early galaxies with CIV and NIV] lines in a large database of publicly-available NIRSpec spectra.  \nThe structure of this paper is as follows. Section 2 provides an outline of the observations and data analysis, and discusses our measurement methods. In Section 3 we describe the spectra of each object in our sample, and we discuss the derived properties of the sample in Section 4. Section 5 discusses the properties of these sources in the context of the reionization-era galaxy population. Finally, we provide a discussion in Section 6, followed by a summary and brief conclusions in Section 7. Throughout this paper we assume a cosmology with Ω 𝑚 = 0 . 3, Ω Λ = 0 . 7, 𝐻 0 = 70km s -1 Mpc -1 , and adopt solar abundances from Asplund et al. (2009, i.e., 𝑍 ⊙ = 0 . 014, 12 + log ( O / H ) ⊙ = 8 . 69). All magnitudes are provided using the AB system (Oke & Gunn 1983).  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \nFigure 1. 𝐽 125 -band images for each of the objects in our spectroscopic sample. Each postage stamp is 2 '' on a side and displayed in an orientation such that north is up and east is to the left. We overlay the position of the NIRSpec shutter in the postage stamp of each object. The listed redshifts were measured from the NIRSpec spectrum.  \n<!-- image -->", '2 DATA AND MEASUREMENTS': 'In this section we provide an overview of the program from which our data derive (Section 2.1), along with the methods used to measure emission and absorption lines from the spectra. In Section 2.2 we detail the methods used to infer properties of the stellar populations within each galaxy.', '2.1 Observations and Reduction': "Spectroscopic data were obtained using JWST /NIRSpec in MultiObject Spectroscopy (MOS) mode targeting two lensing fields, Abell 1703 (RA: 198.7636, Dec: +51.821, Allen et al. 1992) and RXCJ2248-4431 (RA: 342.17972, Dec: -44.5330; Guzzo et al. 2009), as part of Cycle 1 program ID 2478 (PI: Stark). In the Abell 1703 field we obtained data using the G140M/F100LP and G395M/F290LPdispersingmodeswithtotalexposuretimesof6653 sec. and 1576 sec., respectively, using the NRSIRS2RAPID readout mode. Following the position angle (PA) assignment of the observations, the MSA pointing centers were defined based on the locations of our highest-priority targets. For the observations in the Abell 1703 field presented in this analysis, we required that a microshutter slit was placed on A1703-zd6, which was previously spectroscopically confirmed to lie at 𝑧 = 7 . 045 based on Ly 𝛼 (Schenker et al. 2012), CIV, and OIII] (Stark et al. 2015b). We additionally required that the location of A1703-zd6 on the MSA would enable continuous wavelength coverage from Ly 𝛼 to CIII] 𝜆𝜆 1907 , 1909 (9730Å-15530Å in the observed frame) in G140M, and [OII] 𝜆𝜆 3727 , 3729 to [OIII] 𝜆 5007 (29980Å-40280Å in the observed frame) in G395M. The final MSA design was configured to fit these constraints and maximize the number of additional high-priority targets, which comprise a sample of 𝑧 and 𝑖 -drop galaxies selected to lie at 𝑧 ≃ 6 (Richard et al. 2009; Bradley et al. 2012; Smit et al. 2014). Throughout this paper we refer to the additional high-priority targets as A1703-zd2, A1703-zd5.1, A1703zd5.2, A1703-zd4, and A1703-23. Each slit placed upon the targets was composed of three microshutter slitlets, and we obtained data using the standard three-nod pattern, with one exposure for each nod position (Figure 1).  \nWe reduced the spectroscopic dataset using a pipeline com-  \nTable 1. Sample of galaxies observed with NIRSpec in Abell 1703. H-band magnitudes are from Bradley et al. (2012), except for A1703-23 for which we use values from Richard et al. (2009).  \nposed of standard STScI tools (Bushouse et al. 2024) 1 in addition to custom-made routines. The initial processing of the raw uncalibrated ( *\\_uncal.fits ) frames included a bias and dark current correction, as well as ramp fitting to the groups yielding the mean count rate in each pixel. During this stage, jumps in consecutive groups resulting from cosmic rays are identified and masked, including significant events such as 'snowballs' and 'showers'. The resulting full-frame 2D images were corrected for 1 / 𝑓 noise using the /n.pc/s.pc/c.pc/l.pc/e.pc/a.pc/n.pc (Rauscher 2024) package, and the 2D spectrum traced out by each object was cut out from all of the full-frame exposures. We applied a flat-field correction, applied a wavelength solution and absolute photometric calibration using the updated Calibration Reference Data System (CRDS) context to the cutout spectra for each exposure of every object. We subtracted the background from each exposure following the nodding pattern described above. The background-subtracted exposures were then interpolated onto a commonwavelength grid. Each nodded exposure of the targets were then background subtracted and combined including the rejection of pixels that have been flagged in previous reduction stages. These 2Dexposures for each object were then interpolated onto a common wavelength grid and combined to produce the final spectrum. We fit the spatial profile of each 2D spectrum, and used the fitting information to perform an optimal extraction (Horne 1986) yielding the final 1D spectra.  \nWe corrected the emission from each lensed image for the effects of magnification. Our lens model builds on the model published by Zitrin et al. (2010), which used the Light-Traces-Mass approach and included a set of 16 multiply imaged systems to constrain the model, most of which had spectroscopic redshifts and were previously known (see Richard et al. 2009; Limousin et al. 2008). Here we rerun the model using a parametric lens modeling code, namely a revision of the Zitrin et al. (2015) parametric pipeline, which has been extensively used in recent JWST results (e.g., Pascale et al. 2022, 2024; Furtak et al. 2023; Meena et al. 2023). The model consists of two principal components: cluster galaxies and larger-scale dark matter halos. Cluster galaxies are modeled as double Pseudo Isothermal Ellipsoids (dPIEs; Elíasdóttir et al. 2007), scaled by their luminosity following common scaling relations (e.g., Jullo et al. 2007). Two dark matter haloes are incorporated, each modeled as a pseudo isothermal elliptical mass distribution (PIEMD; e.g., Keeton 2001). One halo is initially centred on the brightest cluster galaxy (BCG) but is allowed to roam around the BCG center, with the exact position optimized in the minimization procedure. A second halo is set on a second bright cluster galaxy north-north-west of the BCG, at [13:15:03.1829, +51:49:56.629]. The mass of the BCG  \nTable 2. Derived and measured quantities for galaxies observed with NIRSpec in Abell 1703. Reported stellar population properties are derived from /b.pc/e.pc/a.pc/g.pc/l.pc/e.pc SED fitting, and have been corrected for the effects of gravitational lensing (see Section 2). 𝑎 CIV feature displays significant P-Cygni profile (see Section 3.2).  \nFigure 2. Star-formation rate versus stellar mass for the two primary targets discussed in this work (violet and orange points) inferred from the /b.pc/e.pc/a.pc/g.pc/l.pc/e.pc SED fitting. We compare to the values inferred for GN-z11 (black square; Bunker et al. 2023) and RXCJ2248-ID (orange circle; Topping et al. 2024b), in addition to a population of star-forming galaxies at 𝑧 ∼ 6 -9 identified by Endsley et al. (2023a) (grey circles). For reference, we display the 'main sequence' derived by Speagle et al. (2014) calculated at 𝑧 = 6. Our primary targets lie at significantly elevated SFR relative to the main sequence and the population of early star-forming galaxies at fixed stellar mass.  \n<!-- image -->  \n/circledot  \nand this second bright galaxy, as well as a third, central galaxy at [13:15:07.8290, +51:48:57.955], are modelled independently of the scaling relation, allowing some more freedom. In the minimization the redshifts of the photometric-redshift systems, i.e. those lacking a spectroscopic redshift, are left free to be optimized as well. The minimization is performed in the source plane via a long MCMC and the final model has a (lens-plane) image reproduction rms of 0 '' . 8.  \nEmission line fluxes, centroids, and widths were measured by fitting Gaussian profiles to the observed features. For emission lines that are isolated, we fit a model of a single feature, while emission lines in close proximity in wavelength were fit simultaneously. All of the emission line fitting was provided an initial estimate of the  \nredshift, and its centroid was allowed to vary within one resolution element for emission lines that trace the systemic redshift, and within 500 km/s for resonance and absorption lines. In the case of multiple simultaneous line fits, the best-fit redshift was allowed to vary, however the spacing between the lines was fixed. In some cases, such as for [OII] 𝜆𝜆 3727 , 3729 and NIII] 𝜆 1750, the lines are significantly blended, and line fluxes were measured by integrating over the combined profiles. The continuum level for each line was measured by fitting a linear component to spectral windows on either side of the line. We obtained uncertainties on each property (flux, centroid, FWHM, and EW) by first creating 5000 mock spectra that were derived by perturbing the observed spectrum by its corresponding error spectrum. We repeated the above fitting procedure on each of the mock spectra, and assign uncertainties based on the inner 68th percentile of the recovered distributions of each property.", '2.2 Stellar Populations from SED Fitting': 'Before analyzing the emission line properties, we investigate the basic physical properties of our sample using synthesized broadband photometry from the spectra, following a similar approach to that presented in Topping et al. (2024b). For consistency with previous works, we derive photometry based on the JWST /NIRCam filter throughputs and fit the resulting broadband fluxes using SED models. The A1703 sample was observed using the G140M and G395M gratings, which completely encompass the JWST /NIRCam F115W, F150W, F335M, F356W, F410M, and F444W filters. As such, we utilized synthesized photometry in only these filters for the SED fitting, with the redshift of each object fixed to its spectroscopic value. Briefly, we fit the resulting photometry of both objects using BayEsian Analysis of GaLaxy sEds (BEAGLE Chevallard & Charlot 2016). These models implement updated models of Bruzual & Charlot (2003), which include nebular emission using the prescription detailed in Gutkin et al. (2016). We assume a constant star-formation history (CSFH) with a minimum age of 1 Myr, and adopt a uniform prior on metallicity (ionization parameter) within the range log ( 𝑍 / 𝑍 ⊙ ) ∈ [-2 . 24 , 0 . 3 ] (log ( 𝑈 ) ∈ [-4 , -1 ] ). Finally, we assume an attenuation law following that of the SMC (Pei 1992), and a fixed dust-to-metal ratio of 𝜉 𝑑 = 0 . 3. Throughout this analysis, we adopt galaxy properties (e.g., stellar mass, SFR) that have been corrected for gravitational lensing using the magnification factors provided in Table 2. By assuming a CSFH, the inferred stellar population parameters are representative of the most recent episode of star formation, which will dominate the emission in the rest-frame  \nFigure 3. Presentation of the rest-frame UV and optical spectra for the sample of 𝑧 > 6 galaxies in Abell 1703. Each row corresponds to one object, which is labeled on the side of the right-most panel. Each column is a cutout of the rest-UV spectrum around the wavelengths where prominent features are detected. Vacuum wavelengths are indicated by vertical blue dotted lines rising above each spectrum, which are also labeled at the top of the figure. Transitions typically seen in absorption are indicated by vertical red dashed lines below each spectrum. These absorption lines are labeled in the bottom row of the figure.  \n<!-- image -->  \nUV and optical. This modelling does not account for the possibility that each galaxy hosts a significant older stellar population, which could lead to a significant underestimate of properties such as the stellar mass and SFR (e.g., Whitler et al. 2023). However, this possibility will not significantly impact any of our conclusions.  \nWe present the model properties that reproduce the SEDs in Table 2. The A1703 𝑧 > 6 sample spans a factor of ∼ 170 stellar mass and ∼ 70 in SFR, ranging from log ( M / M ⊙ ) = 7 . 70 -9 . 54 and log ( SFR / M ⊙ yr -1 ) = -0 . 34 -1 . 49, respectively. In Figure 2 we show the SFRs of our sample as a function of their stellar mass along with the larger sample of 𝑧 ≃ 6 -9 LBGs from Endsley et al. (2023a). Our sample roughly extends over the full range of SFRs spanned by the broader 𝑧 ≳ 6 population at log ( M / M ⊙ ) ≳ 8. When assuming a CSFH, we find the youngest galaxy has a CSFH age of just 1 . 6 Myr, while the oldest system has an age of 150 Myr. This suggests that our spectra are sampling galaxies dominated by a very young stellar population in addition to those that have may have a  \nsignificant contribution from more evolved stars. Additionally, the best-fit SEDs imply that our sample has a low average dust content, which is consistent with the typically blue UV slopes observed for galaxies during this epoch (e.g., Topping et al. 2024a; Cullen et al. 2023). However, the emission of some objects in our sample (e.g., A1703-zd4, A1703-23) is best modelled by SEDs that include a nonzero amount of attenuation. In the following section, we describe the spectrum of each galaxy in our sample.', '3 SPECTROSCOPIC PROPERTIES': 'In this section we describe the spectroscopic properties of each object in our sample. Our focus will be on emission and absorption line measurements (i.e., EWs, line ratios), building an empirical physical picture for each individual galaxy. We will then build on and synthesize these results in §4, calculating gas-phase properties  \nFigure 4. Comparison of nitrogen lines detected in deep NIRSpec G140M spectra of RXCJ2248-ID (Topping et al. 2024b) and A1703-zd6 (this work) at 𝑧 = 6 . 106 and 𝑧 = 7 . 044, respectively. Both spectra have been shifted into the rest frame. The vacuum wavelengths of the NIV] doublet and NIII] multiplet members are indicated by vertical blue lines. The line of zero flux for each spectrum is shown as a horizontal dotted line.  \n<!-- image -->  \n(i.e., electron densities, abundance patterns) using the empirical measurements reported in this section.  \nFigure 3 presents spectra of the A1703 galaxies in the rest-UV and rest-optical, at wavelengths where lines are commonly detected. Line fluxes and EWs measured following the methods outlined in Section 2.1 are provided in Table A1 for the G140M spectra, and Table A2 for the G395M spectra. Multiple emission lines are detected in each galaxy. We also detect absorption lines in the two brightest continuum sources (A1703-zd1 and A1703-23).', '3.1 A1703-zd6': "A1703-zd6 was identified as a bright (H = 25 . 9) 𝑧 -band dropout, moderately lensed ( 𝜇 = 3 . 9 + 0 . 1 -0 . 2 ) by the foreground Abell 1703 cluster (Bradley et al. 2012). Schenker et al. (2012) confirmed the redshift ( 𝑧 = 7 . 045) via detection of strong Ly 𝛼 emission (EW = 65 ± 12Å) with Keck/NIRSpec. Motivated by the discovery of high ionization UV lines in metal poor 𝑧 ≃ 2 -3 galaxies (e.g. Erb et al. 2010; Christensen et al. 2012; Stark et al. 2014), a Keck/MOSFIRE J-band spectrum was subsequently obtained covering rest-frame 1433-1680 Å. The observations revealed very strong nebular CIV emission at 1.2458 𝜇 m( 𝜆 1548) and 1.2474 𝜇 m ( 𝜆 1550), indicating the presence of a hard ionizing spectrum rarely seen in lower redshift star forming galaxies (Stark et al. 2015b). The R=1000 NIRSpec spectrum presented in this paper covers the restUV (9700-18900 Å) and rest-optical (28700-52700 Å), providing a new window on the ionizing agents and ionized gas in A1703zd6. Below we will discuss the lines measured in the spectrum and compare to NIRSpec observations of RXCJ2248-ID (Topping et al. 2024b), the other 𝑧 ≳ 6 CIV emitting galaxy known prior to JWST .  \nThe rest-optical spectrum (Figure 3) shows very prominent Balmer lines (H 𝛽 , H 𝛾 , H 𝛿 , H 𝜖 ) and strong [OIII] nebular ( 𝜆𝜆 4959,5007) and auroral ( 𝜆 4363) emission, yielding a systemic  \nredshift of 𝑧 = 7 . 0435. The H 𝛽 EW (423 . 8 Å) is not only at the high end of values seen at 𝑧 ≳ 6 (e.g. Roberts-Borsani et al. 2024) but is also larger than that found in the vast majority of extreme emission line galaxies found at lower redshifts (e.g., Maseda et al. 2018; Tang et al. 2019). The [OIII] 𝜆𝜆 4959 , 5007 EW is similarly elevated (3692 Å), matched by a small number of known 𝑧 ≳ 6 galaxies in existing JWST surveys (Tang et al. 2023; Llerena et al. 2023; Boyett et al. 2024). Indeed, the [OIII]+H 𝛽 EW (4116 Å) is in the upper 0.5% of values seen in the reionization-era population (Endsley et al. 2023a; Matthee et al. 2023). The stellar population synthesis models described in §2.2 reproduce the strong optical lines, provided the age of the stellar population dominating the light is very young (1 . 6 Myr) and the ionization parameter very large (log ( 𝑈 ) = -1 . 7 ± 0 . 3). In this picture, the strong CIV emission is powered by a hard radiation field associated with an extremely young stellar population formed in a recent burst of star formation (Figure 2).  \nIn contrast to the strong nebular emission described above, the A1703-zd6 spectrum does not reveal detection of the [OII] doublet. The lower ionization state gas (if present) appears to not be intense in its emission. The ratio of [OIII] and [OII] fluxes in A1703-zd6 is found to be > 36 (3 𝜎 ), well above the typical values at 𝑧 ≳ 6 (median O32 ≃ 10 in existing spectroscopic samples; Tang et al. 2023; Cameron et al. 2023a; Sanders et al. 2023a). Such a large O32 value is qualitatively consistent with the well-known relationship between [OIII] EW and O32 (Tang et al. 2019; Sanders et al. 2020), whereby extreme ionization conditions (large O32) tends to be found in galaxies with very large [OIII] EW. The large flux ratio of [Ne III] and [OII] ( > 3 . 8; hereafter Ne3O2) presents more evidence for extreme ionization conditions. The lower limits on the O32 and Ne3O2 ratios in A1703-zd6 are consistent with those measured in RXCJ2248-ID (Topping et al. 2024b). In RXCJ2248ID, it was shown that the extremely large O32 value (184) was partially driven by elevated electron densities ( ∼ 10 5 cm -3 ) which result in collisional de-excitation of [OII]. At these high densities, we also expect anomalously strong He I 𝜆 5876 emission, owing to elevated collisional excitation. A1703-zd6 presents a strong He I 𝜆 5876 line, with a flux 26% that of H 𝛽 . We will demonstrate that the densities of A1703-zd6 are indeed very large in § 4.1.  \nThe UV continuum slope of the synthesized SED of A1703zd6 is very blue ( 𝛽 = -2 . 6), comparable to galaxies with similar luminosities at 𝑧 ≳ 6 (Topping et al. 2024a; Cullen et al. 2023; Morales et al. 2024). Such blue colors suggest very little reddening from dust, as evidenced by the low optical depth implied by the BEAGLE population synthesis models fits ( 𝜏 𝑉 = 0 . 01 ± 0 . 01) described in §2.2 (see Table 2). The Balmer line flux ratios present a similar picture, with values (H 𝛾 / H 𝛽 = 0 . 57 ± 0 . 17) consistent with those expected from case B recombination (once the appropriate electron temperature and densities are adopted; see § 4.2). The optical line ratios also indicate the ionized gas is metal poor. We will show in §4.2 that detection of auroral lines suggest a gas-phase metallicity of 12 + log ( O / H ) = 7 . 47, corresponding to 6% Z ⊙ .  \nThe NIRSpec UV spectrum of A1703-zd6 (top panel in Figure 3) reveals clear detection of the CIV emission seen with Keck, along with several other intense emission features (NIV], He II, OIII], NIII], CIII]). The NIRSpec R=1000 data reveal more information than was possible with Keck. Both components of the CIV doublet are detected and resolved with NIRSpec, with an integrated strength (EW=19.4 Å) well in excess of what is seen in nearly all lower redshift star forming galaxies (Senchyna et al. 2019; Izotov et al. 2024). As CIV is a resonant line, we expect the lines to be offset from the systemic redshift. We measure centroids of 1.2461 𝜇 m  \nFigure 5. Rest-UV line ratios compared to photoionization models of HII regions within star-forming galaxies (blue points; Gutkin et al. 2016), and models where a narrow line AGN dominates the ionizing radiation field (red points; Feltre et al. 2016). Galaxies in our A1703 spectroscopic sample are displayed as the colored diamonds. For comparison, we plot the rest-UV line ratios of RXCJ2248-ID (Topping et al. 2024b) and GNz11 (Bunker et al. 2023) as the orange circle and black square, respectively. Line ratios indicated as limits are provided at the 3 𝜎 level.  \n<!-- image -->  \n( 𝜆 1548) and 1.2484 𝜇 m ( 𝜆 1550), both indicating the line is shifted by 240 km s -1 from systemic. The line profile is distinctly different from RXCJ2248-ID, where the blue component was shown to be significantly attenuated by scattering effects (Topping et al. 2024b). The ratio of the two components in A1703-zd6 is near unity, with a doublet ratio of 𝑓 CIV 𝜆 1548 / 𝑓 CIV 𝜆 1550 = 0 . 8 ± 0 . 2 (Table 4). This is well below the intrinsic ratio (2), as expected when the blue component is preferentially scattered by outflowing gas. Nonetheless, the mere presence of the blue component suggests A1703-zd6 faces less attenuation to the transfer of CIV photons than RXCJ2248-ID.  \nThe detection of the strong nitrogen line emission in A1703zd6 follows the recently reported detection in RXCJ2248-ID (Topping et al. 2024b, see also Bunker et al. 2023; Isobe et al. 2023; Castellano et al. 2024), providing more evidence for a connection between hard ionizing sources and nitrogen enhancement. Both components of the NIV] doublet are detected in A1703-zd6, revealing a total EW of 5 . 0 ± 0 . 7Å. The NIV] line profile is very similar to that seen in RXCJ2248-ID, with the red component 2.1 × stronger than the blue component (Figure 4). We will come back to discuss implications for the electron density in § 4.1. The CIV/NIV] flux ratio in A1703-zd6 (3.1) is slightly greater than that seen in RXCJ2248-ID (2.1). We additionally detect emission from the NIII] quintuplet (S/N=3 . 2 𝜎 ) with an integrated EW of 3 . 1 ± 1 . 0Å. We will use these detections to constrain the N/O ratio of A1703-zd6 in § 4.4.  \nThe detection of the auroral OIII] 𝜆𝜆 1660,1666 doublet (EW=11.6 Å) enables an additional constraint on the electron temperature which we will discuss in § 4.2. The [CIII], CIII] 𝜆𝜆 1907,1909 doublet is resolved with a total line strength (EW=8.8 Å) that is comparable to what is often seen in metal poor star forming galaxies (Stark et al. 2017; Hutchison et al. 2019; Tang et al. 2021, 2023). The CIV/CIII] flux ratio is well in excess of unity (3.3), suggesting a very hard spectrum. The He II 𝜆 1640 line is also detected at 6 . 2 𝜎 , indicating a very large EW (6 . 8Å) rare among nearby star forming galaxies (Berg et al. 2016, 2019; Senchyna et al. 2017). Nevertheless, the UV line ratios appear consistent with expectations for stellar photoionization (Figure 5). We will discuss  \nFigure 6. Ly 𝛼 profiles of A1703-zd6 (top; §3.1) and A1703-zd4 (bottom; §3.7). The vacuum wavelength of Ly 𝛼 is displayed as a vertical blue line, and the Ly 𝛼 EW, systemic velocity offset, and escape fraction for each source are presented at the top right.  \n<!-- image -->  \nthese comparisons in more detail by Plat et al. (in prep), but here we mostly note that the hard spectrum is not clearly indicative of an underlying AGN in A1703-zd6.  \nThe deep NIRSpec spectrum confirms the strong Ly 𝛼 line detected from the ground by Schenker et al. (2012). The 𝑅 ∼ 1000 grating combined with the sensitivity of NIRSpec yields a precise systemic redshift and enables the shape of the Ly 𝛼 emission line to be accurately characterized. The Ly 𝛼 detection in A1703-zd6 is presented in Figure 6. The Ly 𝛼 peak velocity offset measured from the 𝑅 ∼ 1000 spectrum is 180 ± 80 km/s. This offset is consistent  \nwith the median peak velocities observed for strong line emitters at 𝑧 > 6 . 5, reflecting the attenuation at line center from the increasingly neutral IGM at these redshifts. The Ly 𝛼 EW (61 . 2 ± 9 . 0Å) is larger what has been seen for a majority of moderately-luminous 𝑧 > 7 galaxies characterized by NIRSpec (e.g., Tang et al. 2023; Saxena et al. 2023). We derive the Ly 𝛼 escape fraction for A1703zd6 from the intrinsic Ly 𝛼 flux inferred as 𝑓 𝑖𝑛𝑡 Ly 𝛼 / 𝑓 H 𝛽 = 24 . 6, which is appropriate for the density and temperature and assuming case B recombination. The resulting escape fraction of 𝑓 Ly 𝛼 esc = 0 . 30 ± 0 . 04 is comparable to typical values at 𝑧 ∼ 5 -6 (Chen et al. 2024) implying significant transmission through the IGM. This value is well above typical escape fractions at 𝑧 ≃ 7 ( 𝑓 Ly 𝛼 esc ≃ 0 . 05 Jones et al. 2024; Nakane et al. 2023; Chen et al. 2024), consistent with expectations for galaxies situated within ionized sightlines.  \nFinally as has been noted previously (Bradley et al. 2012), the HST imaging of A1703-zd6 demonstrates that it is extremely compact (Figure 1). It is unresolved in the image plane in the 𝐽 125 band which has a PSF FWHM of 0 '' . 12 (e.g., Koekemoer et al. 2011). This is consistent with the angular size of the rest-UV continuum trace in the 2D NIRSpec spectrum. These constraints yield an effective radius of < 120 pc in the source plane. We constrain the stellar-mass surface density ( Σ M ∗ ) from this radius and the stellar mass from the synthesized SED, which demonstrates that A1703-zd6 has a very high concentration of stellar mass ( > 3 × 10 3 M ⊙ / pc 2 ). In Figure 7 we compare the mass density of A1703-zd6 to other galaxies and star clusters at high redshift (Adamo et al. 2024; Vanzella et al. 2023, 2022) and the local Universe (Brown & Gnedin 2021). A1703-zd6 shows similar characteristics to other high-redshift galaxies sharing several spectroscopic properties (i.e., GNz11 and RXCJ2248ID Bunker et al. 2023; Topping et al. 2024b). Specifically, all of these systems host very young stellar populations, and display highionization UV lines including CIV, NIV], and NIII] potentially suggesting that high stellar density (and in particular a large density of very young O stars) may be linked to these properties.", '3.2 A1703-zd1': "A1703-zd1 is one of the brightest (H=24.0) 𝑧 -band dropouts known (Bradley et al. 2012), with moderate magnification ( 𝜇 = 3 . 1 + 0 . 3 -0 . 3 ) and a blue Spitzer /IRAC [3.6]-[4.5] color suggesting a photometric redshift of 𝑧 phot ≃ 6 . 6 -7 . 0 (Smit et al. 2014). In spite of its brightness, ground-based rest-UV spectroscopy failed to confirm its redshift, likely indicating A1703-zd1 is characterized by weak Ly 𝛼 (Schenker et al. 2012) and fairly weak UV metal lines (Stark et al. 2015b; Mainali et al. 2018). The redshift of A1703-zd1 was eventually confirmed ( 𝑧 = 6 . 827) in the far-IR, via detection of [CII] 𝜆 158 𝜇 mwith the Northern Extended Millimeter Array (NOEMA) (Molyneux et al. 2022). However no far-IR dust continuum was detected with NOEMA (see also Schaerer et al. 2015).  \nThe NIRSpec observations of A1703-zd1 reveal a suite of strong rest-optical emission lines ([OII], [NeIII], [OIII], H 𝛼 , H 𝛽 , H 𝛾 , H 𝛿 ), indicating a systemic redshift ( 𝑧 = 6 . 8257) consistent with that derived from NOEMA. We do not detect [OIII] 𝜆 4363 in this spectrum. The combined [OIII]+H 𝛽 EWisvery large (1290 Å), placing A1703-zd1 among the upper 25% of the [OIII]+H 𝛽 EWdistribution at 𝑧 ≃ 6 -9 (e.g., Endsley et al. 2023a; Matthee et al. 2023). Reproducing such large EW optical lines with the synthesized SED requires relatively young stellar populations (37 Myr). The best fit suggests a moderate stellar mass (log ( M / M ⊙ ) = 9 . 05 ± 0 . 22), one of the largest in our sample. At such young stellar population ages, ionization-sensitive ratios are usually somewhat elevated (Tang et al.  \nFigure 7. Stellar mass surface density versus effective radius. We display galaxies at 𝑧 > 6 including A1703-zd6 (orange diamond, This Work), RXCJ2248-ID (orange circle, Topping et al. 2024b), and GN-z11 (black square, Bunker et al. 2023), in addition to individual star-forming clumps from the Cosmic Gems (Adamo et al. 2024), Sunrise arc (Vanzella et al. 2023), and the Sunburst arc (Vanzella et al. 2022). The pink points provide measurements of individual young star clusters in the local Universe (Brown & Gnedin 2021). Galaxies with confirmed detections in NIV] are outlined with a dotted circle. These NIV]-emitting systems have significantly higher Σ M ∗ at fixed radius compared to local young clusters, and are even elevated compared to individual star-forming clusters at high redshift.  \n<!-- image -->  \n2019), and indeed we find large values of O32 (8 . 0 + 0 . 8 -1 . 0 ) and Ne3O2 (0 . 62 + 0 . 12 -0 . 13 ) in A1703-zd1. The Balmer-line ratios are consistent with the intrinsic values (H 𝛼 /H 𝛽 = 2 . 6 + 0 . 2 -0 . 3 , H 𝛾 /H 𝛽 = 0 . 48 + 0 . 07 -0 . 07 ; Table 3), implying minimal dust attenuation, in agreement with the upper limits on the FIR dust continuum (Schaerer et al. 2015; Molyneux et al. 2022).  \nThe rest-frame UV spectrum reveals confident detection of the continuum (S/N=5 per resolution element) without strong line emission, consistent with expectations from earlier ground-based observations. The non-detections place stringent upper limits on the rest-UV emission line fluxes (see Table A1). At the spectroscopic redshift of A1703-zd1, the G140M spectrum covers 1250-2400Å in the rest frame. In addition, the detector gap falls at rest-frame wavelengths of 1895-2020Å negating our ability to constrain CIII]. Within the observed wavelength coverage, the typical EW limits (3 𝜎 ) on the UV lines range from 0.7-2.2Å, indicating much weaker emission than the relatively intense optical lines. One example is OIII] 𝜆𝜆 1660 , 1666. The implied total EW limit ( < 2.7 Å at 3 𝜎 ) indicates weaker UV lines than are often seen in galaxies with such large [OIII]+H 𝛽 EWs (Du et al. 2020; Tang et al. 2021). This emphasizes that not all strong optical line emitters should be expected to power strong UV lines. The presence of weak UV lines in galaxies with very young stellar populations is often a signature of relatively elevated gas-phase metallicities, given the sensitivity of the UV transitions to electron temperature. Strong line calibrations (see §4.2) provide some support for this picture, indicating a gas-phase metallicity at the upper end of this sample  \nFigure 8. Rest-frame UV absorption lines in the spectrum of A1703-zd1(top row) and A1703-23(bottom row). In each panel, the systemic wavelength of the absorption feature is indicated by a red dashed line. The offset of the absorption feature is then listed at the top of each panel. The error spectra are displayed as the grey shaded regions. The vertical dark and light shaded regions illustrate the 1 . 5 𝜎 and 3 𝜎 uncertainties for the velocity offset of the line center from systemic.  \n<!-- image -->  \n(12+log(O/H)=7.88). But detection of auroral lines in the UV or optical will be required to confirm this definitively.  \nThe well-detected UV continuum in the A1703-zd1 spectrum is ideal for recovery of absorption lines. We detect both low-ionization (LIS) and high-ionization (HIS) species. Figure 8 presents a zoomed in view of each detected absorption feature in our spectrum, including SiII 𝜆 1260, OI+SiII 𝜆 1303, SiIV 𝜆𝜆 1393 , 1402, CIV 𝜆𝜆 1548 , 1550, and AlIII 𝜆𝜆 1854 , 1862. We discuss the CIV profile in more detail later in this section, as it also includes a stellar wind component. We measure EWs of the LIS features in absorption that range from -0 . 5 to -2 . 0Å, and have an average (blueshifted) velocity offset of -220 ± 80 km / s relative to systemic. The HIS SiIV 𝜆𝜆 1393 , 1402 line has a lower total EW than each of the LIS lines. However, it is blue-shifted from systemic by a similar velocity ( -238 ± 129 km/s). These features signify the presence of outflows in both phases, similar to those commonly seen at 𝑧 ∼ 1 -5 (e.g., Shapley et al. 2003; Steidel et al. 2010; Jones et al. 2012; Pahl et al. 2020; Saldana-Lopez et al. 2023). The spectrum highlights that moderately massive galaxies at 𝑧 ≃ 7 do present significant absorption lines with EWs approaching those seen in 𝑧 ≃ 3 -4 composites.  \nWhile most of the UV absorption lines are relatively narrow with characteristic velocity offsets of ∼ 200 km/s, a much broader kinematic profile is evident in the CIV complex. The absorption in CIV for this object extends to -2000 km/s, and is accompanied by prominent redshifted emission. Rather than interstellar gas, a PCygni profile of this shape and character is suggestive of formation in the hot stellar winds of massive OB stars. Assuming this profile is indeed stellar, we compare it to models for simple stellar populations (SSPs) from the Charlot & Bruzual population synthesis grids (see Section 2.2; and description in Gutkin et al. 2016; Plat et al. 2019; Senchyna et al. 2022). We compare against SSP ages (2-4 Myr) that span a range of H 𝛽 EWs similar to that measured for A1703-zd1, and for which these stellar wind contributions are most prominent at a given metallicity. These models include both stellar light and the contribution from nebular continuum in the FUV. Since CIV is the only line with a significant potential stellar contribution at the observed SNR (consistent with model expectations), rather than detailed fitting we adopt a simple approach of qualitative comparison  \nwith the observed profile, normalized to the median flux level at 1570-1580 Å.  \nThe resulting comparison is displayed in Figure 9 (top panel). It is clear that the young very low metallicity ( 𝑍 = 0 . 001) models ( 𝑍 / 𝑍 ⊙ = 5%) substantially underestimate the observed P-Cygni profile. Slightly older SSP ages at this low metallicity would predict even shallower profiles (Vidal-García et al. 2017), even more inconsistent with the observations. At moderately low metallicities ( 𝑍 = 0 . 004 -0 . 008), a range of individual SSP models at different ages are able to match both the absorption and emission profile reasonably well (underscoring the challenge of using resonant wind lines as a robust age and metallicity tracer in-isolation). Any more complicated star formation history will suppress the strength of this P-Cygni profile relative to the SSPs where it is most prominent (Senchyna et al. 2022). Thus, in the context of these model predictions, this comparison suggests that a massive star population with metallicities of ∼ 20% solar or above is present. However, the strengths of these modeled stellar wind lines are highly sensitive to the uncertain treatment of winds for very luminous metal-poor massive stars, and to the uncertain proportion of such luminous stars in stellar populations. Indeed, observations of FUV-optical spectra of metal-poor dwarf galaxies in the local Universe hosting young bursts of star formation have revealed markedly stronger C /i.pc/v.pc wind profiles in comparison to models constrained in detail by other observables (both nebular line emission and FUV photospheric absorption lines; Senchyna et al. 2021, 2022). One proposed solution is that these metal poor young galaxies have an overabundance of luminous massive stars near the Eddington limit (and thus driving strong winds) not captured in the models, perhaps from a high rate of binary mass transfer and mergers occurring on very short timescales. These very massive stars may contribute to stronger wind lines than are captured in the existing models, allowing the wind profiles to be fit at somewhat lower metallicities than shown in Figure 9. A deeper FUV spectrum would reveal or place limits on additional wind features with different sensitivity to burst age, IMF, abundances, and other key properties (e.g., Vidal-García et al. 2017; Smith et al. 2016), allowing a more robust investigation of the stellar population in A1703-zd1.  \nTable 3. Rest-optical emission-line ratios measured for the A1703 sample. Limits are provided at the 3 𝜎 level, and '-' indicate that at least one line in the line ratio lacks wavelength coverage in the spectrum.", '3.3 A1703-zd2': 'A1703-zd2 is a bright (H=24.9) z-band dropout highly magnified ( 𝜇 = 19 . 9 + 7 . 1 -7 . 0 ) by the foreground cluster (Bradley et al. 2012). No previous spectroscopy of this system has been presented in the literature. The NIRSpec rest-optical spectrum is dominated by prominent [ OIII ] 𝜆𝜆 4959 , 5007 and Balmer (H 𝛽 and H 𝛾 ) emission lines, implying a systemic redshift of 𝑧 spec = 6 . 4267. At this systemic redshift, the G395M spectrum covers wavelengths of 3900-6500Å in the rest frame such that strong lines such as [NeIII], [OII], and H 𝛼 are not observed. The auroral [OIII] 𝜆 4363 line is not detected at the current S/N of the spectrum. The [OIII] 𝜆 5007 (EW = 1096 ± 412Å ) and H 𝛽 (EW = 209 ± 81Å) lines are narrow, with FWHM of 295 km/s and 320 km/s, respectively, which is consistent with the instrumental resolution at 𝑅 ∼ 1000. The total EW of the [ OIII ] + H 𝛽 line complex is 1580Å, which is nearly twice the average value for 𝑧 > 6 galaxies of a similar UV luminosity (Endsley et al. 2023a) but not as large as in A1703-zd6. Driven by the moderately large rest-optical line EWs, the synthesized photometry is best fit by an SED model (see Section 2.2) with a (CSFH) stellar population age of just 17 + 18 -7 Myr.  \nThe H 𝛾 /H 𝛽 ratio (0 . 63 + 0 . 14 -0 . 16 ; Table 3) implies the emission from A1703-zd2 is subject to minimal dust attenuation. The relatively blue UV continuum ( 𝛽 = -2 . 2) provides a consistent picture of minimal dust content and low V-band optical depth ( 𝜏 𝑉 = 0 . 04). Of the available lines in the optical spectrum, the [ OIII ]/ H 𝛽 ratio (5 . 2 + 0 . 5 -0 . 7 ) is the best probe of the gas-phase conditions. The measured value of 5 . 2 + 0 . 5 -0 . 7 is below the typical ratios observed at 𝑧 ≳ 6 (6-8; Sanders et al. 2023a; Tang et al. 2023; Cameron et al. 2023a), likely suggesting low gas-phase metallicities. Using existing strong line calibrations (§4.2), this suggests 12 + log ( O / H ) = 7 . 4, just 1 / 20 𝑍 ⊙ . We infer a stellar mass for A1703-zd2 of only log ( M / M ⊙ ) = 7 . 52 ± 0 . 47, which is the lowest value in the sample. A low mass is consistent with the optical line properties described above, given the empirical trends between both attenuation (e.g., Shapley et al. 2023) and metal content (e.g., Curti et al. 2023) derived at high redshift. We return to interpret the oxygen abundance combined with the other galaxy properties in Section 4.2.  \nThe rest-frame UV spectrum displays prominent detections of the high-ionization CIV and He II lines in emission. Emission features commonly used to probe higher energies (i.e., NV and NeIV) are not covered within the wavelength range of the spectrum (1300-2350Å). The presence of CIV and He II in emission along with the large optical line EWs in A1703-zd2 is qualitatively similar to trends observed in low-metallicity low mass galaxies in the local Universe (e.g., Senchyna et al. 2017). However, the observed CIV EWof 31 . 6 ± 7 . 0Å is considerably larger than nearly all such local  \nsystems (e.g., Senchyna et al. 2019; Izotov et al. 2024), and even exceeds the majority of objects at 𝑧 > 6. This clearly indicates this low mass galaxy is associated with a hard radiation field.  \nWe resolve the CIV doublet profile (Figure 3) into two unresolved emission lines at observed wavelengths of 1.1504 𝜇 m and 1.1521 𝜇 m. Relative to the systemic redshift, these two doublet members are offset by + 142 km / s (Table 4). We find that the CIV 𝜆 1548 flux is nearly twice that of CIV 𝜆 1550, consistent with the intrinsic doublet ratio. This indicates a low amount of scattering and relatively clear channels for CIV to propagate through the ISM. This is likely to contribute to the elevated CIV EW measured in this system. The He II 𝜆 1640 line also exhibits a very large EW (14 . 0 ± 6 . 1Å) that exceeds values measured for even the most extreme CIV emitters locally (Izotov et al. 2024) and at high-redshift (Topping et al. 2024b; Castellano et al. 2024). These two prominent emission lines characterize the hardness of the ionizing spectrum from 48-54 eV, and their ratio (CIV/He II = 3 . 1 ± 1 . 0) is consistent with models driven by star formation or AGN activity (see Figure 5; Feltre et al. 2016; Gutkin et al. 2016). Finally, the CIII] (total EW=13.0Å) and OIII] (total EW < 12 . 1Å) lines are found to be weaker than both He II and CIV. The UV line ratios implied by these features (Figure 5) are also consistent with both stellar and AGN photoionization. Whereas other sources analyzed in this paper are more confidently in the regime of stellar photoionization, a deeper spectrum of A1703-zd2 is required to put more robust constraints on the powering mechanism of the radiation field.', '3.4 A1703-zd5.1': "A1703-zd5.1 is one the faintest (H=25.7) galaxies in our sample and is also the most magnified ( 𝜇 = 33 . 6 + 14 . 0 -13 . 7 ) of the 𝑧 > 6 systems identified by Bradley et al. (2012). As such, it has the lowest luminosity of all high-redshift galaxies we observed in the Abell 1703 field. The galaxy is a close neighbor of A1703-zd5.2 (see §3.5 below), situated only 0 '' . 5 away from its companion in the image plane. Bradley et al. (2012) interpreted their close separation as a signature of separate clumps within a single galaxy, however it is also possible that they are different galaxies currently in the process of merging. Despite the two star forming components falling within the same slit, we leverage the high spatial resolution of JWST /NIRSpec to extract separate spectra for each system. Here we present the our observations of A1703-zd5.1, while A1703-zd5.2 is described in §3.5. Despite the apparent faintness of A1703-zd5.1, NIRSpec confirms detections in H 𝛽 , [ OIII ] 𝜆 5007, and H 𝛼 at a spectroscopic redshift of 𝑧 = 6 . 4260. We take this as the systemic redshift.  \nA1703-zd5.1 notably has very low optical-line EWs of [OIII] 𝜆 5007 (267Å) and H 𝛼 (157Å). The synthesized SED is relatively flat (in f 𝜈 ) across the rest-UV and rest-optical, suggest-  \nFigure 9. Prominent broad absorption in the resonant CIV 1548,1550 doublet is detected in A1703-zd1 and A1703-23 (spectra in greyscale). Comparison to stellar wind profiles for SSP models from Charlot & Bruzual (colored curves) reveal reasonable agreement with the P-Cygni profile in A1703-zd1, provided a fairly high model stellar metallicity 𝑍 ≳ 0 . 004 ( 𝑍 / 𝑍 ⊙ ≳ 20%). This suggests we may be viewing a relatively metal-rich galaxy (consistent with the weakness of the UV nebular emission lines), or one in which very luminous stars driving prominent winds are over-represented. The absorption in A1703-23 is of similar magnitude, but with minimal emission near systemic, in contrast to the stellar wind profiles. Deeper FUV spectroscopy of such galaxies probing other wind lines at similar or higher resolution promises a clearer view of the massive stars populating these systems.  \n<!-- image -->  \nively young stellar population. This source appears to be a spectroscopic example of the weak-line population identified photometrically in Endsley et al. (2023a). It is challenging to explain these SEDs with simple CSFH models, as young ages generally are accompanied by large EW rest-optical lines. Endsley et al. (2023b) demonstrated that such systems can be readily explained as having recently experienced a sudden decrease in their star-formation activity. To address these complications in A1703-zd5.1, we e  \nmploy a two-component star-formation history (TcSFH; Endsley et al. 2023a) with the additional flexibility required to fit the observed emission. These models parallel our setup described in Section 2.2, except with a SFH composed of a delayed-tau model (SFR ∝ 𝑡 · 𝑒 -𝑡 / 𝜏 ) combined with a decoupled CSFH episode over the most recent 20 Myr.  \nThe TcSFH model provides a much better fit ( 𝜒 2 = 1 . 6) when compared to the CSFH model ( 𝜒 2 = 9 . 4). Furthermore, the best-fit TcSFH SED model implies a light-weighted age of 70 Myr, and stellar mass of just 10 7 . 74 M ⊙ . The relatively low instantaneous sSFR (3 . 3 / Gyr -1 ) underscores the recent downturn in star formation. The O3 line ratio (3 . 1 + 0 . 9 -1 . 6 ) is the lowest value in our sample, indicative of a very low oxygen abundance (0 . 04 𝑍 ⊙ ; see Section 4.2 below). This is perhaps not surprising given the very low luminosity of this galaxy. We note that while low metallicity will contribute to the weak [OIII], it should boost the strength of H 𝛼 and other Balmer lines (see discussion in Endsley et al. 2023a). Hence low metallicity on its own cannot explain the weak emission lines in this galaxy (as evidenced by the poor 𝜒 2 for the CSFH models). The recent downturn in the SFH allowed in the TcSFH models is key for explaining the weakness of both H 𝛼 and [OIII].  \nWe note that models with ionizing photon escape can also explain the weakness of the emission lines in such SEDs. Following the approach taken in Topping et al. (2022), we fit the SED using the density-bounded f esc models of Plat et al. (2019). In short, these models self-consistently compute the impact of density-bounded HII regions on the continuum and lines for different escape fractions (see Plat et al. 2019 and implementation in Topping et al. 2022). The results indicate that reasonable solutions can be achieved for moderate escape fractions, although the reduced 𝜒 2 (3.3) is not quite as low as for the TcSFH models.  \nThe rest-frame UV spectrum of A1703-zd5.1 is relatively featureless, with no detections in the commonly-observed lines (e.g., CIV, He II, CIII]). We place tight upper limits on the EWs of these features ranging from 3-6Å (Table A1). The absence of restUV emission lines is consistent with the low sSFR inferred from the best-fit TcSFH model described above, which implies that the massive stars that power these lines have begun to disappear following the recent downturn in SFR. We do observe a tentative He II 𝜆 1640 emission line (S/N=2.5) that appears slightly broadened (FWHM=450 km/s) relative to the instrumental line-spread function (see Figure 3). If confirmed, this would potentially indicate a broad stellar He II profile. This would be more challenging to explain in the TcSFH models, where recent star formation is subdominant, but it could be explained by the leakage models, as significant massive star populations would still be present. Ultimately, deeper spectroscopy will be required to constrain the flux and shape of He II in order to discern its origin.", '3.5 A1703-zd5.2': "A1703-zd5.2 is the brighter ( 𝐻 =25.3) neighbor to A1703-zd5.1, is situated only 0 '' . 5, and is magnified by a factor of 𝜇 = 33 . 5 + 17 . 8 -16 . 9 . The rest-optical spectrum reveals significant detection in multiple Balmer lines (H 𝛼 , H 𝛽 , H 𝛾 , H 𝛿 ), in addition to [ OIII ] 𝜆𝜆 4959 , 5007. These lines indicate a systemic redshift of 𝑧 spec = 6 . 4295, which is very similar to its close neighbor, A1703-zd5.1. The full extent of the optical spectrum spans from 3850-6680Å in the rest frame, thus lacking constraints on the low-ionization [OII] and [SII] lines. All of the emission lines detected in the rest optical are narrow, with an average FWHM (310 km/s) that is consistent with the instrumental resolution.  \nTable 4. with a prominent nebular component.  \n𝑎 Flux measurements are given in units of 10 -19 erg / s / cm 2  \nProperties of the CIV emission line for three objects in our sample .  \nDespite the close proximity between A1703-zd5.1 and zd5.2, the rest-optical spectra differ significantly. The H 𝛼 EW(651 ± 227Å) is over 4 × that of A1703-zd5.1, while the [ OIII ] 𝜆 5007 EW (749 ± 188Å) is a factor of 2.7 larger. The observed SED is well fit by models with CSFH stellar population ages of 23 + 45 -17 Myr. Each of the Balmer line ratios is consistent with the absence of attenuation from dust, with H 𝛼 /H 𝛽 (2 . 8 + 0 . 2 -0 . 3 ) being the most constraining. We draw a similar conclusion from the relatively blue UV-continuum slope ( 𝛽 = -2 . 4). The only strong rest-optical line ratio probed by NIRSpec is O3. We measure a low value of 4 . 4 + 0 . 3 -0 . 4 , which suggests the strength of [OIII] is being weakened by the low gas-phase metallicity. The strong line calibration of Sanders et al. (2023a) suggests this value of O3 corresponds to 12 + log ( O / H ) = 7 . 26 ± 0 . 28 (Table 3).  \nThe rest-UV spectrum reveals significant CIV emission (Figure 3). We additionally achieve detections at ≳ 5 𝜎 in OIII] 𝜆 1666 and [CIII] 𝜆 1907 within the observed wavelength range of 13002400Å. CIV is the strongest of the rest-UV lines in A1703-zd5.2, comprising a total EW of 10 . 5 ± 1 . 0Å. While not as extreme as the CIV emission observed in some other 𝑧 > 6 galaxies (e.g., A1703zd6, RXCJ2248-ID, GHZ2), CIV EWs approaching 10Å are rare locally (e.g., Izotov et al. 2016a, 2024; Senchyna et al. 2017). Furthermore, we resolve the CIV emission into two narrow components (1.1504 𝜇 mand 1.1524 𝜇 m)that are redshifted by + 80 km/s relative to systemic (Table 4). The relative fluxes of the two components ( 𝜆 1548 / 𝜆 1550 = 2 . 0) and the small velocity offset (from systemic) likely indicates the resonant CIV line may face minimal attenuation. This may suggest the observed emission is close to the intrinsic value produced in the HII regions. We use the OIII] 𝜆 1666 ([CIII] 𝜆 1907) detections to infer the direct gas-phase oxygen abundances (electron density) in Section 4.2 (Section 4.1) below. Here, we use the CIV/He II ratio ( > 5 . 6) to demonstrate a significant decrease in the ionizing radiation field at energies above 47 eV. Other rest-UV line ratios provide a similar picture, with the CIV/CIII] ( > 3 . 2) and OIII]/He II ( > 2 . 7) line ratios (Figure 5) suggesting a stellar origin for the ionizing radiation. It is striking how different A1703-zd5.2 is from its close ( ≃ 450 pc) neighbor A1703-zd5.1. While A1703-zd5.2 is in the midst of a significant episode of star formation with intense nebular emission, A1703-zd5.1 appears to have very weak emission lines, potentially following a recent downturn in star formation.", '3.6 A1703-23': 'A1703-23 was first identified as a bright (H=23.8) 𝑖 -band dropout (Zheng et al. 2009; Richard et al. 2009), in which subsequent followup reported Ly 𝛼 emission at 𝑧 ≃ 5 . 8 (Richard et al. 2009). Further spectroscopy obtained by Stark et al. (2015b) tar-  \nIII] 𝜆𝜆 1907 , 1909 line using Keck/MOSFIRE based on the Ly 𝛼 redshift, however no strong emission was detected. The strongest rest-optical lines detected with JWST /NIRSpec are [OIII] 𝜆 5007, H 𝛼 , and H 𝛽 , at wavelengths of 3.5489 𝜇 m, 4.6518 𝜇 m, and 3.4456 𝜇 m, respectively. Based on these strong lines, we measure a redshift of 𝑧 spec = 6 . 0862, slightly different from the singleline measurement presented in Richard et al. (2009). Unlike the other galaxies in our sample, the Balmer-line ratios (H 𝛼 / H 𝛽 = 3 . 4 + 0 . 3 -0 . 4 ) imply a modest dust attenuation, which is consistent with the moderately red UV slope ( 𝛽 = -1 . 5) and optical depth inferred from the SED ( 𝜏 𝑉 = 0 . 17 + 0 . 04 -0 . 04 ). The EWs of [OIII] (415 ± 22Å), H 𝛽 (46 ± 4Å), and H 𝛼 (193 ± 61Å) are lower than typical galaxies at 𝑧 ≃ 6 (Endsley et al. 2023a). The synthesized SED fit implies a moderate age stellar population (150 + 33 -61 Myr) and a relatively large stellar mass (log ( M / M ⊙ ) = 9 . 54 ± 0 . 16). The implied SFR (25 M ⊙ / yr) is broadly consistent with the star forming main sequence at 𝑧 = 6 (Speagle et al. 2014). This is in contrast to most of the other galaxies we targeted in Abell 1703 which sit above the main sequence with much younger stellar populations. 2  \nThe rest-frame UV spectrum reveals CIII] emission in addition to a significantly detected continuum (4.5 per resolution element). As a result, we detect several features in absorption (Figure 8). Si II 𝜆 1526 is the most strongly detected absorption line, and traces low-ionization along with neutral gas (e.g., Shapley et al. 2003). The Si II 𝜆 1526 absorption (EW= -2 . 0 ± 0 . 6Å) is moderately stronger than lines measured from stacks of typical 𝑧 ∼ 3 -7 galaxies (e.g., Jones et al. 2012), however the strongest absorption lines are observed in systems with the largest attenuation from dust and highest metallicity (Jones et al. 2012; Du et al. 2021). These features are qualitatively similar to the properties of A1703-23, indicating that reasonable strong absorption lines should be present in a subset of 𝑧 ≃ 6 galaxies. The Si II 𝜆 1526 line center is blue-shifted from the systemic redshift ( Δ 𝑣 = -130 ± 29km/s) indicative of outflowing gas. In contrast, the high-ionization HIS SiIV 𝜆 1393 , 1402 absorption line is detected (total EW=1 . 2 ± 0 . 4Å) at the systemic redshift ( Δ 𝑣 = 4 ± 97km/s). The difference in velocity structure indicates the LIS and HIS lines are kinematically distinct. Notably, this is different than what is observed in A1703-zd1, where the two phases are kinematically more similar.  \nAs in A1703-zd1, absorption in CIV is also detected in A170323. We compare this absorption profile to the same set of stellar models as described in Section 3.2 in Figure 9. Interestingly, the absorption in A1703-23 is almost absent accompanying emission near the rest wavelengths of the doublet. This is in strong contrast with the profiles typical of the spherical outflows driven by massive star winds and thus the predictions of the SSP models we compare with. This morphology is closer to that observed in broad absorption-line quasar spectra (e.g. Allen et al. 2011), though with fairly modest ∼ 2000 km/s velocity for such systems. While limited by SNR, a pure absorption profile would be suggestive of interstellar gas entrained in some manner of outflow rather than stellar wind absorption; but this would stand in contrast to the near-systemic absorption in Si /i.pc/v.pc. It is also possible that the peculiar profile is a result of a relatively weak stellar wind profile combined with particularly  \nstrong interstellar absorption in C /i.pc/v.pc near-systemic. Higher-SNR and potentially higher-resolution data would help separate these contributions and enable clearer constraints to be placed on the putative stellar component directly.', '3.7 A1703-zd4': 'A1703-zd4 is a 𝑧 -band dropout identified by Bradley et al. (2012) to have a photometric redshift of 𝑧 phot = 8 . 4 + 0 . 9 -1 . 4 . Multiple spectroscopic campaigns have observed A1703-zd4 using Keck/MOSFIRE targeting the Ly 𝛼 , CIV, He II, OIII], and CIII] emission lines (Stark et al. 2015b; Mainali et al. 2018). These observations did not yield detections of any emission lines, however the EW limits ( < 6 -20Å at 5 𝜎 ) were sufficient to distinguish A1703-zd4 from extreme UV-line emitters such as A1703-zd6. The JWST /NIRSpec spectrum provides the first spectroscopic confirmation at a redshift of 𝑧 spec = 7 . 8727 based on the [OIII] 𝜆 5007 and H 𝛽 lines. For this spectroscopic redshift, the H-band magnitude (H=25.4) and modest magnification ( 𝜇 = 2 . 0 + 0 . 1 -0 . 1 ) correspond to an absolute UV magnitude of M UV = -20 . 7 (i.e., ≃ 𝐿 ∗ UV Bouwens et al. 2021). At the stellar mass (log ( M / M ⊙ ) = 8 . 50 ± 0 . 27) and SFR (log ( SFR / M ⊙ yr -1 ) = 1 . 06 + 0 . 52 -0 . 45 ) inferred from the SED fit, we find that A1703-zd4 is elevated above the star-forming main sequence by 0 . 4 dex (Speagle et al. 2014). The best-fit SED implies a non-zero impact from dust attenuation ( 𝜏 𝑉 = 0 . 07 + 0 . 04 -0 . 04 ), which is consistent with the UV slope measured using NIRSpec ( 𝛽 = -2 . 0) and relative fluxes of the strongest observed Balmer lines (H 𝛾 / H 𝛽 = 0 . 39 + 0 . 09 -0 . 10 ; Table 3).  \nWe measure large [OIII] 𝜆 5007 and H 𝛽 EWs of 1317 ± 432Å and 148 ± 46Å, respectively, which are both consistent with the high sSFR of the SED. A1703-zd4 is one of the two objects in our sample with an [OII] detection enabling the O32 ratio to be calculated. As the best-fit SED and Balmer-line ratios imply the presence of dust, we correct the [OIII] and [OII] fluxes for attenuation using the Cardelli et al. (1989) attenuation curve and assuming an intrinsic H 𝛾 / H 𝛽 = 0 . 47. We derive a resulting O32 ratio of 4 . 0 + 0 . 9 -0 . 9 , while we calculate a line ratio prior to the dust correction of 5 . 6 + 0 . 5 -0 . 7 . Starforming galaxies in this epoch commonly have O32 in excess of the A1703-zd4 value (O32 ≃ 10; Tang et al. 2023; Sanders et al. 2023a; Cameron et al. 2023a), with a clear trend between [OIII] EW and O32 (Tang et al. 2023). A1703-zd4 appears broadly consistent with this relation.  \nIn the rest-frame UV we detect clear Ly 𝛼 emission as well as [CIII] 𝜆 1907 at S / N > 3. The [CIII] 𝜆 1907 EW is just 3 . 6 ± 0 . 7Å, which is easily matched by EELGs at 𝑧 ∼ 2 with comparable optical line properties (e.g., Tang et al. 2019). While no other rest-UV metal line is detected (EWs < 3 -6Å), we compute limits on the CIV/CIII] ( < 0 . 3) and CIII]/He II ( > 4 . 8) line ratios. Both of these constraints support a population of massive stars as source of ionizing radiation in A1703-zd4 (Figure 5). We also detect a weak Ly 𝛼 line with an EW of 3 . 4 ± 0 . 9Å. We calculate a velocity offset of + 240km / s for Ly 𝛼 relative to the systemic redshift set by the optical lines. The Ly 𝛼 EW in A1703-zd4 is relatively small, and the typical EWs of galaxies at lower redshift with the same optical line EWs suggest a significant attenuation of Ly 𝛼 photons from A1703-zd4. We compute the Ly 𝛼 escape fraction following the same method outlined in Section 3.1, and recover a value of just 𝑓 Ly 𝛼 esc = 0 . 03 ± 0 . 01, indeed suggesting the observed line is a small portion of the intrinsic emission.  \nFigure 10. Electron densities inferred from the CIII] and NIV] doublet ratios. Theoretical CIII] and NIV] line ratios calculated for different densities and assuming T e = 20000K are displayed as the solid and dashed lines, respectively. Measured line ratios and inferred densities for the high-redshift galaxies are shown as colored markers. Densities inferred from the CIII] doublet are placed along the solid line, while densities inferred from NIV] are plotted on the dashed line. A1703-zd2 and A1703-zd5.2 are only detected using the CIII] line, while A1703-zd6 and RXCJ2248-ID are detected in CIII] and NIV].  \n<!-- image -->  \nFigure 11. Relative HeI 𝜆 5877/H 𝛽 flux as a function of distance from the star-forming main sequence. Galaxies at 𝑧 > 6 are displayed as coloured symbols, while the locus traced by SDSS galaxies is shown as the grey histogram. Typical galaxies in the local universe have an average HeI/H 𝛽 line flux of 0 . 12, regardless of their location in the SFR versus M ∗ plane. While the objects in our sample within 1 dex of the main sequence are consistent with this line ratio, A1703-zd6 and RXCJ2248-ID are significantly He I enhanced, resulting from the high densities preferentially boosting He I relative to H 𝛽 from collisional effects.  \n<!-- image -->  \nFigure 12. Inferred oxygen abundance as a function of stellar mass for CIV emitters at 𝑧 > 6 (diamonds) from the A1703 field in addition to RXCJ2248ID (Topping et al. 2024b), compared to local CIV emitters (squares; Izotov et al. 2016a,b, 2018a,b; Senchyna et al. 2017, 2022; Wofford et al. 2021; Schaerer et al. 2022b). Objects are color-coded based on their CIV EWs as indicated by the colorbar on the right. We display the mass-metallicity relation derived locally (black line; Curti et al. 2020). The dotted black line represents an extrapolation of the local MZR from Curti et al. (2020) below the mass range probed by their sample. The green line and white pentagons are the best-fit MZR and composite measurements of 𝑧 ∼ 6 -10 galaxies from Curti et al. (2023). The local systems scatter about the MZR, with the smallest and largest CIV emitters lying above and below the MZR, respectively. All of the 𝑧 > 6 CIV emitters have below-average metallicities for their stellar mass.  \n<!-- image -->  \n/circledot', '4 SAMPLE PROPERTIES IN ABELL 1703': 'In the previous section, we presented the spectroscopic properties of our 𝑧 ≳ 6 sample in the Abell 1703 field, describing each source individually. Here we describe the electron densities (§4.1), gas-phase oxygen abundances (§4.2), carbon and nitrogen abundances (§4.3 and §4.4). Finally we discuss existing constraints on aluminum enhancements, with the goal of identifying an additional anomalous abundance pattern in globular clusters (§4.5). In all subsections below, our broad goal is using the new data to improve our understanding of the properties of galaxies with hard ionizing sources and nitrogen-enhanced ionized gas.', '4.1 High electron densities in ionized gas': 'Recent work has revealed ionized gas with extremely high electron densities ( 𝑛 𝑒 ≃ 10 5 cm -3 ) in a subset of reionization era galaxies, several orders of magnitude larger than typical densities at lower redshifts (Bunker et al. 2023; Senchyna et al. 2023; Topping et al. 2024b). This was first inferred for GNz11 based on the NIV] doublet ratio (Bunker et al. 2023; Senchyna et al. 2023; Maiolino et al. 2023a) and was then demonstrated in RXCJ2248-ID based on three density-sensitive doublets in the UV (Topping et al. 2024b). In GNz11 it has been suggested the densities may be related to the broad line region of an AGN (Maiolino et al. 2023a), but in RXCJ2248-ID the high ionized gas densities appear likely related to dense conditions during a strong burst of star formation (Topping et al. 2024b).  \nA1703-zd6 has a similar spectrum as RXCJ2248-ID, with strong CIV emission and extremely powerful rest-optical emission lines ([OIII]+H 𝛽 EW = 4116 Å) suggesting a very young stellar population (1.6 Myr). We detect two density-sensitive doublets in the NIRSpec observations presented in §3.1. We infer densities using the /g.pc/e.pc/t.pcT/e.pc/m.pcD/e.pc/n.pc procedure from the /p.pc/y.pc/n.pc/e.pc/b.pc software package (Luridiana et al. 2015). The NIV] doublet ratio ( 𝑓 1483 / 𝑓 1486 = 0 . 47 + 0 . 19 -0 . 21 ) suggests the ionized gas is very dense (19 + 20 -9 × 10 4 / cm 3 ), comparable to that seen in RXCJ2248-ID. The CIII] doublet ratio ( 𝑓 1907 / 𝑓 1909 = 0 . 49 + 0 . 19 -0 . 16 ) reveals a similar picture (9 . 4 ± 4 . 2 × 10 4 cm -3 ). Given the higher ionization potential of NIV] (47.4 eV), its emission originates nearer to the ionizing source when compared to the CIII] emission (24.4 eV). The electron density follows a similar trend and increases toward the inner regions of the H II regions, leading higher-ionization probes to yield higher densities (e.g., Kennicutt 1984; Kewley et al. 2019; Berg et al. 2021).  \nThere are several potential physical explanations for the elevated densities. We may expect high ionized gas densities ( ≃ 10 5 cm -3 ) to arise in a short window after a strong burst of star formation, when gas surface densities are sufficiently high to create an over-pressurized environment for the HII regions (e.g., Lehnert et al. 2009). In this picture, we would only observe very dense gas in the subset of galaxies with the largest sSFRs (and hence the largest optical line EWs). On the other hand, the UV-based tracers of electron densities (CIII], NIV]) have long been known to probe denser gas than the more typically-used doublets in the optical ([OII], [SII]) (James et al. 2014, 2018; Mainali et al. 2023). It is thus possible that the UV lines may reveal elevated densities similar to RXCJ2248ID and A1703-zd6 in all reionization era galaxies, not just those undergoing strong bursts.  \nWe can attempt to distinguish between these two scenarios using our Abell 1703 sample. In addition to A1703-zd6, we detect the density-sensitive CIII] doublet in A1703-zd2 and zd5.2, both with lower H 𝛽 EWs than A1703-zd6. The spectra of the three galaxies reveal distinctly different doublet flux ratios and hence different densities (Figure 10). A1703-zd5.2 has the lowest H 𝛽 EW(170 Å) and has a 1907/1909 flux ratio ( > 1 . 48) implying moderately low ionized gas densities ( < 1500 cm -3 ). A1703-zd2 has a larger H 𝛽 EW (209 Å) and a 1907/1909 flux ratio (1 . 10 + 0 . 21 -0 . 22 ) indicating moderately large densities (2 × 10 4 cm -3 ). Taking these results together with the spectra of A1703-zd6 and RXCJ2248-ID, we find that the UV-based densities do increase with H 𝛽 EW (albeit in small samples), as would be expected if densities are elevated in the strongest bursts. The highest densities ( > 10 5 cm -3 ) appear limited to subset of galaxies with H 𝛽 EWs in excess of 300 Å, corresponding to extremely young stellar populations ( ≲ 2 Myr) formed in a recent upturn of star formation.  \nSignatures of very high densities are not limited to the restUV. In galaxies with ionized densities of 10 5 cm -3 , we expect the [OII] emission line to be collisionally de-excited, leading to high O32 indices (184 and > 36 in RXCJ2248-ID and A1703-zd6, respectively). We also expect to see strong He I 𝜆 5876 emission owing to enhanced collisional excitation of this transition (see e.g., Yanagisawa et al. 2024). In RXCJ2248-ID, He I 𝜆 5876 was found to be 21% as strong as H 𝛽 , roughly 2 × larger than is typically found in local star forming galaxies (Abazajian et al. 2009). The enhancement in the line ratio can be fully explained by the elevated gas densities (see Topping et al. 2024b). We find similarly enhanced He I 𝜆 5876 emission in A1703-zd6, with a line flux that is 26% that of H 𝛽 . The other galaxies in our sample are not detected in He I 𝜆 5876, as expected given the lower ionized gas densities. The  \nFigure 13. N/O as a function of O/H for stars and HII regions in the local Universe and galaxies at high redshift. Local HII regions trace out a tight locus that asymptotes at log ( N / O ) ≃ -1 . 4 at low metallicities. We display measurements from a sample of N-enriched high-redshift galaxies including A1703-zd6 (orange diamond) with nitrogen abundances 10 × higher than HII regions at fixed O/H. This parameter space is shared by individual stars in the globular cluster NGC 6752 indicating that both types of objects may undergo common formation processes.  \n<!-- image -->  \nupper limits we are able to place on the He I 𝜆 5876/H 𝛽 flux ratios in these systems ( ≲ 0.10-0.15) are typically 2 × lower than that found in A1703-zd6 and RXCJ2248-ID (Figure 11), consistent with flux ratios found in nearby galaxies. These results suggest that deep restoptical spectroscopy may be able to identify high density sources via elevated He I 𝜆 5876 emission (see also Yanagisawa et al. 2024).', '4.2 Gas-phase oxygen abundances': 'Oneof the primary explanations of nebular CIV emission at 𝑧 ≳ 6 is hard radiation from a population of very metal poor massive stars. In this case, we would expect the CIV emitting galaxies to be linked to metal poor gas, as is often found in nearby galaxies (Senchyna et al. 2017, 2022; Wofford et al. 2021; Schaerer et al. 2022b). Photoionization modeling in Stark et al. (2015b) suggested that reproducing the strong CIV emission in A1703-zd6 may require metallicities approaching 2% of the Solar value (12 + log ( O / H ) = 7 . 05), but these inferences were based on very limited information available from ground-based observations. With the NIRSpec observations presented in this paper, we can now derive a direct gas-phase metallicity for A1703-zd6, while also investigating implications for the metal content in the other lensed systems in our sample.  \nWe quantify gas-phase oxygen abundances for each object in our sample using the available emission-line constraints. The 11 𝜎 [ OIII ] 𝜆 4363 detection in A1703-zd6 provides the best O/H constraint in our sample, while A1703-zd5.2 and A1703-23 are detected in the auroral [ OIII ] 𝜆 1666 line. For these systems we derive direct metallicities following the methods described in Topping et al. (2024b, see also Izotov et al. 2006). For the remaining galaxies in our sample, we infer the oxygen abundance using the O3 and O32 strong-line calibrations derived from high-redshift galax-  \ners et al. 2023b). Results are listed in Table 3. Gas-phase oxygen oxygen abundances are generally found to be low, ranging from 12 + log ( O / H ) = 7 . 20 to 7 . 98.  \nWe start with our primary target, the strong CIV emitting galaxy A1703-zd6. Adopting the density found in Section 4.1, we recover an electron temperature of T e ( O ++ ) = 23000 ± 3200K, which is consistent with low metallicity galaxies in this epoch (e.g., Laseter et al. 2023; Katz et al. 2023). For this electron temperature and the [ OIII ] 𝜆 5007 and H 𝛽 fluxes, we derive an oxygen abundance of 12 + log ( O / H ) = 7 . 47 ± 0 . 19. This value represents oxygen only in the doubly-ionized state, which we expect to dominate the abundance in this system given the very large O32 ratio ( > 36). We derive a consistent metallicity when considering the O3 ratio (6 . 4 + 0 . 3 -0 . 4 ), which yields an abundance of 12 + log ( O / H ) = 7 . 60 ± 0 . 30 using the Sanders et al. (2023b) calibrations. Calibrations using the O32 and Ne3O2 ratios indicate much lower metallicities, however these line ratios are likely severely impacted by the density of the system (see Section 4.1). Overall, these results suggest the gas in A1703zd6 is metal poor (6% 𝑍 ⊙ ), consistent with lower redshift results (e.g., Senchyna et al. 2017, 2022; Wofford et al. 2021; Schaerer et al. 2022b) and as would be expected if there is a low metallicity stellar population contributing to the hard radiation field. However what stands out most about A1703-zd6 in our sample is not its metallicity, but rather its very young inferred stellar population age (1 . 6 + 0 . 5 -0 . 4 Myr) and correspondingly large H 𝛽 EW(424 ± 112 Å). Both are at the extremes of what has been seen by JWST in the 𝑧 ≳ 6 -8 galaxy population. Presumably the coupling of the extremely young stellar population and low gas phase metallicity are mutually required for strong CIV emission (e.g., Senchyna et al. 2019).  \nWe have also detected strong CIV emission ( > 10 Å) in two additional galaxies in the Abell 1703 field (A1703-zd2 and A1703zd5.2). As expected, the new CIV emitters are among the most metal poor in the sample. The O3 ratios of these two systems (5 . 2 + 0 . 5 -0 . 7 and 4 . 4 + 0 . 3 -0 . 4 ) are well below average (O3=7-10) at 𝑧 ≃ 6 -8 (e.g., Tang et al. 2023; Sanders et al. 2023a; Cameron et al. 2023a). Using the calibrations described above, these likely point to low metallicities: 12 + log ( O / H ) = 7 . 35 ± 0 . 30 (A1703-zd2) and 7 . 26 ± 0 . 28 (A1703-zd5.2). We additionally detect the OIII] UV-based auroral line in A1703-zd5.2, suggesting a direct method metallicity (12 + log ( O / H ) = 7 . 20 ± 0 . 28) consistent with the O3-based value described above. The inferred stellar populations ages of A1703zd2 (17 Myr) and A1703-zd5.2 (23 Myr) are the next two youngest in our sample after A1703-zd6.  \nWe note that galaxies in Abell 1703 with metal poor gas but without young stellar populations do not show strong CIV emission. A1703-zd5.1 provides the best example of such a source, with ionized gas that appears metal poor (12 + log ( O / H ) = 7 . 19 ± 0 . 43) but comparatively weak CIV emission ( < 4.3 Å). As pointed out in § 3.5, this source appears best fit by a recent downturn in star formation (or ionizing photon leakage), resulting in a spectrum without strong CIV emission in spite of the low metallicity. In contrast, A1703-zd4 has a relatively young stellar population age (23 Myr) but gas-phase metallicity at the upper end of our sample (12 + log ( O / H ) = 7 . 98 ± 0 . 31). Here we see no strong CIV emission ( < 4 . 9 Å). Collectively these results emphasize that the presence of strong CIV emission is limited to sources that show both young stellar populations ( ≲ 25 Myr) and low metallicity gas ( ≲ 0.1 Z ⊙ ), consistent with trends seen at lower redshift (e.g., Senchyna et al. 2017, 2022; Wofford et al. 2021; Schaerer et al. 2022b). We will come back to discuss the demographics of the CIV population through investigating the full public spectroscopic database in §5.  \nFigure 14. Total NIV]+NIII flux compared to the total CIV+CIII] flux as a function of electron density. We display all objects from our sample for which the density can be constrained from the CIII] ratio. The two galaxies with nitrogen-line detections (and thus higher N/C line ratios) have the highest densities, while the galaxies with more typical CIII] electron densities have (NIV]+NIII)/(CIV+CIII]) ratios that are < 0 . 2.  \n<!-- image -->  \nTo place the 𝑧 ≳ 6 CIV-emitting galaxies in context, we consider their position with respect to the mass-metallicity relationship (MZR; Curti et al. 2023), as derived from JWST spectroscopy of 𝑧 ≃ 6 -10 galaxies (Figure 12). We also include several additional strong 𝑧 ≳ 6 CIV emitters (Topping et al. 2024b) as well the local star forming galaxies with CIV detections (Izotov et al. 2016a,b, 2018a,b; Senchyna et al. 2017, 2022; Wofford et al. 2021; Schaerer et al. 2022b). It is clear in the figure that the galaxies with strong CIV emission at 𝑧 ≳ 6 typically fall below the MZR, implying that they are both younger and more metal poor than the parent sample of early galaxies. The low redshift galaxies tend to have much weaker CIVEWs( ≃ 1 -3 Å; c.f. Izotov et al. 2024) and lower stellar masses ( ≲ 10 6 M ⊙ ) than existing detections at 𝑧 ≳ 6. The mismatch in the mass of the CIV emitters at low and high redshift arises because the low metallicities required to power hard radiation fields ( ≲ 0.1 Z ⊙ ) are possible in moderately massive galaxies in the reionization era ( ≃ 10 8 M ⊙ ). In the local universe, these low metallicities are only found in much lower mass galaxies. As a result, there is potentially a significant difference between the 𝑧 ≳ 6 CIV emitters and the local reference sample. Future work is required to investigate whether these differences lead to distinctly different spectra at fixed metallicity.', '4.3 C/O abundances': "The abundance of carbon (relative to oxygen) encodes insights into the past star formation history (Henry et al. 2000; Esteban et al. 2014; Berg et al. 2016, 2019) and integrated yields of earlier generations of stars (e.g., Garnett et al. 1995; Yin et al. 2011; Berg et al. 2019). Early JWST spectra have constrained carbon-to-oxygen (C/O) ratios in a small number of early galaxies with deep UV constraints on emission from carbon and oxygen lines (ArellanoCórdova et al. 2022b; Jones et al. 2023; D'Eugenio et al. 2023;  \nTopping et al. 2024b). The detections (and upper limits) of the UV carbon and oxygen lines in our sample allow us to build on these studies. Here, we focus on galaxies in the Abell 1703 field with detections in CIV, CIII], and OIII] (A1703-zd6 and A1703-zd5.2) which will provide the most robust constraints on C/O.  \nBriefly, we infer the C/O abundance ratio by first computing C ++ / O ++ from the observed CIII] and OIII] emission lines using /p.pc/y.pc/n.pc/e.pc/b.pc (Luridiana et al. 2015). In calculating this abundance ratio, we assume an electron temperature and density inferred from OIII] 𝜆 1666/[OIII] 𝜆 5007 and the CIII] doublet, respectively (see Table 3. The CIII] doublet in A1703-zd5.2 is in the low-density regime, providing only an upper limit on the electron density. For this object, we assume a density of 1000 cm -3 , which is consistent with typical densities observed in this epoch (e.g., Sanders et al. 2023b; Laseter et al. 2023; Katz et al. 2023). We note that the inferred abundance ratios vary minimally if a different assumption on density is imposed. These emission lines yield log ( C ++ / O ++ ) abundance ratios of -1 . 20 ± 0 . 12 and -1 . 13 ± 0 . 15 for A1703-zd6 and A1703-zd5.2, respectively.  \nBoth systems show strong CIV emission (19 . 4 ± 1 . 4 Å in A1703-zd6 and 10 . 5 ± 1 . 0 Å in A1703-zd5.2), implying a significant amount of the carbon budget exists in the triply-ionized state. The CIV emission is subject to resonant scattering, and the observed line fluxes may underestimate the true values. As discussed in §3, the line profiles of both galaxies give indications of modest scattering. Here we consider two bounding cases. On one hand, we compute ionic ratios assuming the observed CIV EWs are the intrinsic values. For the other limiting case, we assume that 50% of the CIV flux has been scattered out of the line of sight, providing a conservative bound on the fraction of carbon in the triply-ionized state. We obtain C 3 + / C ++ ionization fractions using /p.pc/y.pc/n.pc/e.pc/b.pc by computing the C 3 + and C ++ emissivities at the densities and temperatures defined above. Using the observed CIV emission as intrinsic we obtain lower-bound C 3 + / C ++ of 0 . 08 ± 0 . 07 (A1703-zd6) and 0 . 06 ± 0 . 08 (A1703-zd5.2) and upper-bound values of 0 . 38 ± 0 . 12 (A1703-zd6) and 0 . 36 ± 0 . 14 (A1703-zd5.2). These give the likely range of ionization fractions. To select a value within this range to use for our calculations, we consider an approach used in Jones et al. (2024). Given that the redder component of the CIV doublet is less susceptible to attenuation (and that both profiles appear to show only modest signs of scattering), we can adopt its flux as intrinsic and compute the total CIV flux using the theoretical doublet ratio ( 𝜆 1548/ 𝜆 1550=2). This suggests C 3 + / C ++ values of 0 . 27 ± 0 . 08 (A1703-zd6) and 0 . 13 ± 0 . 11 (A1703-zd5.2), between the two boundary cases quoted above. We will use these values in our calculations below.  \nUsing the ionic abundances computed above, we derive a total abundance ratio of log ( C / O ) = -0 . 74 ± 0 . 18 and log ( C / O ) = -0 . 79 ± 0 . 17 for A1703-zd6 and A1703-zd5.2, respectively. Both suggest ionized gas that is not substantially enriched in carbon relative to oxygen, well below ( ≃ 4 × ) the solar value (log ( C / O ) ⊙ = -0 . 23; Asplund et al. 2009). Similarly low C/O ratios are commonly inferred for metal poor galaxies locally (Berg et al. 2016, 2018, 2020; Ravindranath et al. 2020; c.f., Dufour et al. 1988) and at 𝑧 ≃ 2 (Erb et al. 2010; Christensen et al. 2012; Stark et al. 2014; Berg et al. 2019). Similarly, many systems in the reionization era have significantly sub-solar C/O (e.g., Jones et al. 2023; ArellanoCórdova et al. 2022b; Stiavelli et al. 2023; Hu et al. 2024), consistent with our results in this paper. The low C/O observed in these systems are consistent with enrichment from supernovae, without significant pollution from AGB stars (e.g., Yin et al. 2011; Berg et al. 2019). The absence of such evolved stellar populations are consistent with  \nthe young ages and massive stars implied by the rest-UV and optical spectra described in Section 3. We note that C/O ratios approaching the solar value are also expected at extremely low metallicities, possibly reflecting yields from population III stars (e.g., Akerman et al. 2004; Carigi & Peimbert 2011). We find no evidence of these in the JWST spectra analyzed above.", '4.4 N/O abundances': 'In Topping et al. (2024b), we demonstrated that the strong CIV emitting galaxy RXCJ2248-ID has an anomalously large N/O ratio, similar to that reported by Bunker et al. 2023 in GNz11 (see also Isobe et al. 2023; Marques-Chaves et al. 2024; Navarro-Carrera et al. 2024). We detect similar nitrogen line emission (NIV] and [NIII]) in the CIV emitting galaxy A1703-zd6 presented in this paper. We have showed that this galaxy has metal poor gas (12 + log ( O / H ) = 7 . 47 ± 0 . 19) and a sub-solar C/O ratio (log ( C / O ) = -0 . 74). Here we first constrain the N/O ratio in A1703-zd6, before commenting briefly on the relative nitrogen abundances in the other systems in our sample. As we also detect CIV in two additional galaxies, our measurements allow us to further explore the potential link between galaxies with hard radiation fields and significant nitrogen enhancements.  \nOur methodology is similar to that discussed in detail in Topping et al. (2024b) (and largely parallels our derivation of C/O described above). We will expand on our abundance calculations in Plat et al. 2024 in prep. We infer N/O in A1703-zd6 and place limits on N/O for systems with OIII] detections and nitrogen line non-detections (A1703-zd5.2, A1703-23) using the NIV ]+ NIII OIII ] line ratio (see e.g., Senchyna et al. 2023; Cameron et al. 2023a). Briefly, in deriving the N/O abundances we assume the electron temperature inferred using the highest S/N auroral emission line ( [ OIII ] 𝜆 4363 or OIII ] 𝜆 1666). For A1703-zd6 and A1703-zd5.2, we fix the electron density to the value inferred from the NIV] and CIII] doublets, respectively. Lacking a density-sensitive measurement in A1703-23, we assume a value of 1000 cm -3 in our abundance calculations, motivated by typical densities in similar sources at these redshifts (e.g., Sanders et al. 2023b; Laseter et al. 2023; Katz et al. 2023). We stress that adopting different densities does not have a large systematic effect on our abundance estimates.  \nWe find that the NIV], NIII], and OIII] line fluxes in A1703zd6 imply a super-solar nitrogen abundance, with a relative abundance of log ( N / O ) = -0 . 6 ± 0 . 3. Figure 13 demonstrates that this is moderately higher than the solar value (log ( N / O ) ⊙ = -0 . 86 ) , Asplund et al. 2009) and nearly 10 × higher than typical values seen in HII regions at the metallicity of A1703-zd6. For A1703zd5.2 and A1703-23 we follow the methodology described above using the 3 𝜎 upper limits on NIV] and NIII] and find constraints of log ( N / O ) < -0 . 9 and log ( N / O ) < -0 . 8, respectively. While we cannot rule out the possibility that these two systems have a minor nitrogen enhancement relative to the local HII regions, their N/O constraints are consistent with the local locus.  \nThese measurements offer new insight into the nature of the nitrogen-emitting galaxies that have been found in the reionization era with JWST . Our Cycle 1 program has now confirmed nitrogen-enhanced gas in both of our primary targets (A1703-zd6 and RXCJ2248-ID), each originally selected only on strong CIV emission in ground-based UV spectroscopy. This potentially hints at a connection between the processes leading to the formation of hard ionizing sources and those leading to the overabundance of nitrogen. We have shown that both galaxies are likely in the midst  \nof strong bursts of star formation with spectra dominated by the light from very young stellar populations ( ≲ 2 Myr). This results in rest-optical emission lines that are among the most prominent (upper 1% of [OIII]+H 𝛽 EW) seen in existing samples of reionization era galaxies (e.g., Endsley et al. 2023a; Matthee et al. 2023). Both galaxies also appear to be host compact star clusters with very dense star formation conditions, the latter revealed by very high electron densities ( ≃ 10 5 cm -3 ) associated with the density-sensitive emission line doublets in the UV. These high densities might correspond to a short period ( ≃ 1 -3 Myr) after strong bursts where HII regions remain confined by very high gas densities that led to the strong burst. At later evolutionary stages, a variety of feedback processes will disperse the birth cloud gas, and the HII regions will acquire more typically-observed densities ( ≃ 10 2 -3 cm -3 ).  \nIn this picture, the observed link between nitrogen enrichment and hard ionizing sources (the latter responsible for CIV or He II emission) may have its origin in the extremely dense conditions that arise in a brief window ( ≲ 1-3 Myr) following strong bursts in the reionization era. The high stellar densities ( > 3 . 6 × 10 4 M ⊙ / pc 2 in RXCJ2248-ID; Topping et al. 2024b) that arise in the clusters formed during these bursts will enhance dynamical interactions (i.e., runaway collisions, binary mass transfer) that may lead to the formation of very massive stars and hot stripped stars, boosting the population of hard ionizing sources relative to less dense environments. Such dynamical processes may also help explain how the ISM is polluted with its overabundance of nitrogen relative to oxygen (e.g., Gieles et al. 2018; Martins et al. 2020; Charbonnel et al. 2023).  \nOur observations also unveiled two new CIV emitting galaxies (A1703-zd5.2, A1703-zd2), with CIV EWs similar to A1703-zd6. According to the picture described above, we may expect both to show nitrogen lines, but our observations reveal no NIII] or NIV] emission in either galaxy. In Figure 14, we plot the upper limit on the flux ratio of nitrogen (NIV] and NIII) and carbon lines (CIII] and CIV) for A1703-zd5.2 ( < 0 . 17) and A1703-zd2 ( < 0 . 19). The flux limits confirm that these two CIV emitters are significantly weaker nitrogen-line emitters than RXCJ2248-ID and A1703-zd6, where this ratio was found to be 1.5-2.0 × greater. Clearly the nitrogen line emission is not present in all strong CIV emitters. While both A1703-zd5.2 and A1703-zd2 are dominated by relatively young stellar populations (23 + 45 -17 and 17 + 18 -7 Myr, respectively), they are not as young as A1703-zd6 (1 . 6 + 0 . 5 -0 . 4 Myr) or RXCJ2248-ID (1 . 8 + 0 . 7 -0 . 4 Myr). Similarly, the ionized gas densities appear elevated ( ≃ 1 -2 × 10 4 cm -3 ) but are 5-10 × lower than the densities in A1703-zd6 and RXCJ2248-ID. It is plausible that these new systems correspond to a slightly later evolutionary state, or that they simply experienced less intense bursts where densities were not elevated sufficiently to initiate the conditions required for seeding the ISM with the nitrogen enhancement. Larger samples will be required to confirm this picture in greater detail. We will come back to this in §5, investigating the nitrogen line emission in a large public spectroscopic sample.', '4.5 Signatures of additional abundance anomalies': 'The 𝑧 ≳ 6 nitrogen line emitters may provide our first opportunity to directly study the nature of the polluters responsible for the anomalous abundances in globular clusters (see e.g., Bastian & Lardo 2018, for a recent review). If the N/O and N/C ratios are truly related to those found in globular clusters, we should also see evidence of the additional abundance anomalies associated with', '18': 'hydrogen burning. Marques-Chaves et al. (2024) have suggested that aluminum enhancements (from the MgAl cycle) are likely to be the most readily detectable in early galaxies. The most accessible probes of aluminum in these galaxies are the AlIII 𝜆 1670 and AlIII 𝜆𝜆 1854 , 1862 lines. These features are most commonly observed in absorption in high redshift star forming galaxies (e.g., Steidel et al. 2016), but are occasionally seen in emission toward some AGNs (Alexandroff et al. 2013). If aluminum is overabundant in the ISM of the nitrogen line emitters, we may expect to see the Al III lines in emission.  \nIn RXCJ2248-ID, we presented an upper limit on the strength of the AlIII 𝜆𝜆 1854 , 1862 doublet, with each individual component likely having an EW below 2 . 0 Å. Here we extend this analysis to A1703-zd6. We search for Al III emission near the expected observed wavelengths of the doublet. The spectrum of A1703-zd6 displays an emission feature at the expected wavelengths of the AlIII 𝜆 1862 emission line. We closely inspect the spectrum surrounding this feature and find that one exposure suffers from a large residual resulting from the subtraction from a snowball event (see Section 2). To remedy this, we repeat our data reduction procedures while removing the affected exposures. The resulting spectrum still exhibits an emission feature (EW=2.7 Å) at the expected wavelength of the AlIII 𝜆 1854 line. As the S/N is low (2.6) we consider this feature to be tentative. We do not detect any emission at the location of the AlIII 𝜆 1862 line, consistent with expectations given the predicted flux ratio ( 𝜆 1854/ 𝜆 1862=1.5) of the doublet.  \nTo quantify the Al abundance implied by this tentative feature in A1703-zd6, we consider the AlIII 𝜆 1854 / OIII ] 𝜆 1666 line ratio, for which we measure 𝑓 AlIII / 𝑓 OIII ] = 0 . 36. For the metallicity and ionization parameter we infer for A1703-zd6, the observed 𝑓 AlIII / 𝑓 OIII ] indicates an aluminum enhancement that is 30-100 × the solar value (log ( Al / O ) ⊙ = -2 . 24, Asplund et al. 2009). A more detailed treatment of the abundance ratios will be presented by Plat et al. (in prep). But we note that this level of aluminum enhancement is consistent with the abundance pattern seen in a subset of individual globular cluster stars (e.g., Carretta et al. 2018; Carlos et al. 2023). Deeper spectroscopy of A1703-zd6 will ultimately be required to verify this tentative feature.', '5 DEMOGRAPHICS OF 𝑍 > 4 UV LINE EMITTERS': 'The discovery of enhanced N/O and N/C ratios in early galaxies has been one of the more surprising early findings of JWST (Bunker et al. 2023; Isobe et al. 2023; Marques-Chaves et al. 2024; Topping et al. 2024b; Castellano et al. 2024), potentially providing a signpost of globular clusters in formation (Senchyna et al. 2023; Charbonnel et al. 2023; Vink 2023). In this paper, we have presented more evidence suggesting that the nitrogen emitters also tend to have a population of hard ionizing sources and very dense ionized gas (see §4.3 and Topping et al. 2024b).  \nIn this section, we attempt to place the strong nitrogen and carbon line emitters in a broader context, establishing how common these systems are at very high redshift and comparing their basic properties to that of the full galaxy population. To achieve this goal in a self-consistent manner, we assemble a large spectroscopic sample of 𝑧 ≳ 4 galaxies with NIRSpec observations and characterize the incidence of strong UV line emission. While the nitrogen emitters observed as part of this program (i.e., A1703-zd6 and RXCJ2248-ID) were observed due to their strong CIV emission, this large spectroscopic sample was not constructed with the same prerequisite. The data reduction and construction of the NIRSpec  \n-  \n-  \n-  \nFigure 15. Spectra for galaxies with previously-unreported CIV or NIV] detections presented in the rest frame. The top, center, and bottom panels shows the spectra for JADES-202208, JADES-210003, and CEERS-1019, respectively, as indicated by their labels on the right (see also Table 5). The spectra of JADES-202208 and JADES-210003 were obtained using the 𝑅 ∼ 1000 G140M grating, and CEERS-1019 was obtained using the 𝑅 ∼ 2700 G140H grating. The spectrum of JADES-210003 between 15651710Å is covered by the NIRSpec detector gap. The wavelength range of the CEERS-1019 spectrum is narrowed to better show the emission lines at high resolution, and none of the redder rest-UV lines (i.e., He II 𝜆 1640, OIII] 𝜆𝜆 1660 , 1666, NIII] 𝜆 1750) are detected at > 3 𝜎 . Vacuum wavelengths of emission lines in the rest-UV are indicated by vertical blue lines that are labelled at the top of each panel. Emission lines that are detected have their fluxes shaded blue. A 5 𝜎 detection of NIV] 𝜆 1486 in CEERS-1019waspreviously reported by Larson et al. (2023) and Isobe et al. (2023), now detected at 17 𝜎 in the 𝑅 ∼ 2700 spectrum. This spectrum also presents the first detection of CIV in CEERS-1019 using the NIRSpec grating.  \n<!-- image -->  \nspectroscopic catalog are briefly described in §5.1. We then discuss the statistics and properties of the CIV-emitters in §5.2 and that of the nitrogen emitters in §5.3.', '5.1 Public JWST /NIRSpec spectra at 𝑧 > 4': "WeassembletheNIRSpecMOSobservationsfromtheCosmicEvolution Early Release Science (CEERS 3 ; ERS-1345, PI: S. Finkelstein; Finkelstein et al. (in prep), also see Finkelstein et al. 2022), the Director's Discretionary Time (DDT) observations (DDT-2750, PI: P. Arrabal Haro; Arrabal Haro et al. 2023a,b) in the CEERS field, the JWST Advanced Deep Extragalactic Survey (JADES 4 ; GTO1180, GTO-1181, and GTO-1210, PI: Eisenstein; Eisenstein et al. 2023a), the JADES Origins Field (JOF; GO-3215, PI:Eisenstein;  \nEisenstein et al. 2023b), and the Ultradeep NIRSpec and NIRCam ObserVations before the Epoch of Reionization (UNCOVER 5 ; GO2561, PIs: I. Labbé & R. Bezanson; Bezanson et al. 2022) and GO-4287 (PI: Mason). A full description of the sample selection and analysis methods applied to these public spectroscopic samples will be provided in Tang et al. (in prep), however we provide a brief description in this section. We obtained final 2D and 1D spectra following the reduction procedure described in §2.1. We also collect the HST /ACS+ JWST /NIRCam photometry available from Endsley et al. (2023b) in the CEERS field, the JADES DR2 photometric catalog in the GOODS-S field (corresponding to sources observed in programs GTO-1180 and 1210, and GO-3215; Rieke et al. 2023), photometry from the DAWN JWST Archive (Brammer 2023; Valentino et al. 2023) in the GOODS-N field (for GTO-1181), and the UNCOVER DR2 catalog from Weaver et al. (2024).  \nWe follow the procedure outlined in Tang et al. (2023) to identify spectroscopic redshifts based on the H 𝛼 (for 𝑧 ≲ 6 . 7) in addition to [ OIII ] + H 𝛽 (for 𝑧 ≲ 9 . 4) emission lines. For the small number of systems at yet higher redshift, we utilized the strongest emission line (typically [NeIII], [OII], or CIII]) if present, otherwise we adopted the redshift implied by the Ly 𝛼 continuum break. Break redshifts typically have larger uncertainties ( 𝜎 𝑧 ∼ 0 . 1) than those measured from emission lines (e.g., Jones et al. 2024; Hainline et al. 2024), however this does not significantly impact our results. In total, these efforts yield a parent sample of 737 galaxies with spectroscopic redshifts spanning 𝑧 = 3 . 90 to 13 . 36 (median 5 . 29). This sample does not include the CIV emitter targeted in RXCJ2248ID (Topping et al. 2024b) or any of the Abell 1703 sources presented in this paper. We will comment on how the detections in these samples compare to the overall population trends.  \nFollowing the methods outlined in Tang et al. (2023, 2024); Chen et al. (2024), we fit the HST /ACS+ JWST /NIRCam photometry of each of the 737 galaxies using BEAGLE following the model setup described in §2.2 with the redshifts fixed to the spectroscopic value. We utilize the best-fit SEDs to define the absolute UV magnitudes ( 𝑀 UV ) as well as to infer stellar population properties including stellar mass, SFR, and age from their best-fit SEDs. In addition, we measure [ OIII ] 𝜆 5007 and H 𝛽 EWs for each galaxy in the sample. In systems where the rest-optical continuum is detected at > 3 𝜎 , the optical-line EWs are measured directly from the spectrum. For these objects, we compare the measured EWs to the values inferred from the best-fit SED and achieve an excellent agreement. In the absence of a significant rest-optical continuum detection in the spectrum, we adopt the [ OIII ] 𝜆 5007 and H 𝛽 EWs derived from the best-fit SED.  \nFigure 15 presents the M UV and combined [ OIII ] + H 𝛽 EWs for the 737 galaxies comprising the spectroscopic sample. The absolute UV magnitudes of the spectroscopic sample range from M UV = -22 . 6 to -16 . 1, with a median value of -19 . 6. This M UV distribution is consistent with the UV-bright end of the UV luminosity distribution of photometrically-selected objects (e.g., Endsley et al. 2023a), the latter also including fainter systems (down to M UV = -16 . 0) and having a luminous median M UV ( -18 . 0). Both the photometric and spectroscopic samples span the similar ranges of [ OIII ] + H 𝛽 EW of 150-6000Å and 45-8300Å respectively, and the median [ OIII ] + H 𝛽 of the spectroscopic sample (710Å) is well-matched to that of the photometric sample at fixed M UV .  \nWe select galaxies from the spectroscopic sample displaying CIV or NIV] emission lines using a combination of visual inspec-  \ntion and S/N criteria. First, we measure emission-line fluxes by shifting the PRISM and grating spectrum into the rest frame, and fitting Gaussian profiles at the vacuum wavelength of each line. The centroid of the profiles were allowed to vary by one resolution element in the PRISM spectra ( ≃ 10000 km/s), and 1000 km/s in the grating spectra to allow for small variations in the wavelength solution in addition to ensure capturing of resonance lines that do not necessarily trace the systemic redshift (e.g., CIV). Uncertainties are set by width of the flux distribution derived by repeatedly fitting profiles to 1000 instances of the spectrum perturbed by its error. Each > 5 𝜎 CIV and NIV] emission line was visually inspected to remove artifacts and poor fits to the data, and we verified that each emission line yielded consistent fluxes when covered by multiple observing modes. The resulting robust sample comprises 11 galaxies detected in either CIV (10 galaxies; Section 5.2) or NIV] (3 galaxies; Section 5.3) spanning from 𝑧 = 5 . 449 up to 10 . 621, which are presented in Table 5. The spectroscopic redshifts for each of these systems are measured from rest-optical emission lines, rather than the continuum break. Among the 9 galaxies detected in CIV, the measured equivalent widths range from 4.9 Å up to 47.6 Å with a median of 25.4 Å, comparable to what is observed in A1703zd2 and A1703-zd6. When present, the NIV] emission lines have EWs that span a smaller range (5.1-8.3Å), which may contribute to their smaller sample size. Finally, we additionally select two (four) galaxies with tentative NIV] (CIV) detections, which satisfy only one of our S/N and visual inspection criteria (Table 5). The spectra of several of our robust sources have been previously presented in the literature, as indicated in the Table. Additionally, we identify two new galaxies (JADES-202208 and 210003) with either a CIV or NIV] detection (Figure 15).  \nIn addition, we present a deep, high-resolution ( 𝑅 ≃ 2700) NIRSpec spectrum of CEERS-1019 (i.e., EGSY8p7; Zitrin et al. 2015; Roberts-Borsani et al. 2016) in the rest-frame UV that reveals CIV emission (Whitler et al. in prep). Previous spectroscopy using the 𝑅 ∼ 1000 grating from CEERS was not sufficiently deep to detect the CIV feature (e.g., Larson et al. 2023; Tang et al. 2023). 6 We measure a total EW for the CIV emission in EGSY8p7 of 4 . 6 + 0 . 5 -0 . 6 Å. The CIV profile of EGSY8p7 is asymmetric with no sign of the blue doublet component. This is reminiscent of the line profile of RXCJ2248-ID and easily explained by resonant scattering effects (see Topping et al. 2024b). It is well-fit by a skewed Gaussian profile redshifted from the systemic wavelength. The properties of the fit are also similar to the RXCJ2248-ID line; we measure a FWHM of 310 km/s after correcting for the instrumental resolution ( ∼ 100 km/s), and a velocity offset from CIV 𝜆 1550 of the emission peak of + 190 km/s. Finally, the emission extends to yet higher velocities, with the best-fit profile extending above the error spectrum up to velocities of +800 km/s.  \nIn the following subsections, our primary goal is quantifying the incidence of CIV and NIV], independent of the powering mechanism of the lines. We do briefly consider whether any of our sources clearly stand out as having lines powered by AGN photoionization. Here we use the UV diagnostics discussed in §3.1 (e.g., Feltre et al. 2016; Mingozzi et al. 2022; Scholtz et al. 2023). Our measurements reveal CIV/He II ratios ( ≳ 1) and CIII]/He II  \nratios ( ≳ 2) that are are consistent with stellar photoionization, as we also found for RXCJ2248-ID and A1703-zd6. We will discuss these in more detail in Plat et al. (2024, in prep). Of course, very high redshift AGN may have different spectra than those predicted in these models, and some level of mixing between AGN and stellar photoionization is also plausible. We therefore do not rule out contribution from AGN photoionization in these sources. We note that one source (JADES-954) has been confirmed as a broad line AGNinMaiolino et al. (2023b). We will discuss this source in more detail below. One additional source (JADES-58975) was listed in the narrow line AGN sample of Scholtz et al. (2023). We find that this spectrum has UV lines that are plausibly consistent with a stellar origin (OIII]/HeII = 1 . 4, [Ne IV] EW < 30Å), but this does not significantly impact our demographic investigation below (see also, Curti et al. 2024).", '5.2 The Demographics of CIV Emitters at 𝑧 ≳ 4': 'We have presented new evidence in this paper that links nitrogen line emitters to those with strong CIV emission (see also Topping et al. 2024b; Castellano et al. 2024). In this subsection, we attempt to quantify the fraction of 𝑧 ≳ 4 galaxies with strong CIV emission, establishing how commonly hard ionizing sources are present in early galaxies. This builds on earlier efforts to establish the CIV emitting fraction in the galaxy population with ground based spectroscopy (Mainali et al. 2018).  \nWe have identified 11 galaxies with robust CIV detections among the 737 systems in our 𝑧 > 4 spectroscopic sample. Table 5 presents the properties of each of these systems. Seven of these CIV detections have been previously reported, including CEERS1019 (EGSY8p7 Larson et al. 2023; Isobe et al. 2023; Tang et al. 2023; Marques-Chaves et al. 2024), JADES-3991 (GNz11 Bunker et al. 2023; Maiolino et al. 2023a), CEERS-397 (Harikane et al. 2023a), JADES-954 (Maiolino et al. 2023b), JADES-1899 (Witstok et al. 2024; Tang et al. 2024), JADES-13176 (Cameron et al. 2023c), JADES-58975 (Scholtz et al. 2023; Curti et al. 2024), and UNCOVER-10646 and UNCOVER-22223 (Fujimoto et al. 2023b). The CIV emission in JADES-210003 and JADES-202208 have not been previously reported, and we present their spectra in Figure 15.  \nThe detection rate of CIV emitting galaxies (11 among 737 systems) is extremely low, as would be expected if the hard ionizing sources required to power strong CIV are only present in a small subset of the parent population. Of course, not all 737 galaxies have sufficiently deep spectra to place a meaningful limit on the CIV EW. To achieve a more robust CIV emitter fraction, we identify spectra with the depth to detect CIV with an EW > 10Å at the > 5 𝜎 level. Of the 737 𝑧 ≳ 4 galaxies with spectra described in Section 5.1, we identify 99 systems with sufficient CIV sensitivity. With this more refined sample, we still find that a very small fraction (8%) of 𝑧 > 4 galaxies show robust CIV emission with an EW in excess of 10Å.  \nTo understand the low success rate in detecting CIV emission, we compare the M UV and [ OIII ] + H 𝛽 EWs of the CIV emitters to the full spectroscopic sample, as well as to a larger NIRCambased photometric sample of 𝑧 ∼ 6 -9 galaxies from Endsley et al. (2023a) in Figure 16. The eight CIV emitters span a fairly broad range of UV luminosities (M UV of -21 . 6 to -17 . 7), with a median absolute UV magnitude of M UV = -19 . 6). In contrast, the [OIII]+H 𝛽 EWs of the CIV emitters are far from typical. Nearly all of the CIV emitters have very large [ OIII ] + H 𝛽 EWs (median 3000Å), indicating they are on the tail of the full distribution of  \ngalaxies. 7 According to the 𝑧 ≃ 6 -9 [OIII]+H 𝛽 EWdistribution in Endsley et al. (2023a), the median [OIII]+H 𝛽 EWexhibited by the CIV emitters corresponds to the upper 2% of the population. There is one CIV emitter in our sample (JADES 1181-954) that shows relatively small [ OIII ] + H 𝛽 EW(679 Å), distinguishing it from the rest of the sample. This source is a broad line AGN with a very red optical continuum (Maiolino et al. 2023b). It is likely that the AGN is contributing significantly to the continuum flux density in the rest-optical, decreasing the [ OIII ] + H 𝛽 EW from what would have been observed from the host galaxy.  \nThe large optical line EWs are consistent with expectations for extremely young stellar populations, formed in a recent strong burst of star formation. Indeed, the BEAGLE SEDs that describe these systems show a median constant star formation age of just 5 Myr. These results are consistent with our findings for A1703-zd6 and RXCJ2248-ID, where [OIII]+H 𝛽 EWs are large (4120 and 3100 Å, respectively) and stellar population ages are young (1 . 6 + 0 . 5 -0 . 4 and 1 . 8 + 0 . 7 -0 . 4 Myr, respectively). Based on these results, we suggest that the rarity of strong CIV emission is likely due to the [OIII]+H 𝛽 EW (or age) threshold required to power the high ionization lines. Only a small fraction of the galaxy population is found with the extremely young stellar population ages that are linked to the large [OIII]+H 𝛽 EWs associated with strong CIV emission. If we consider only those objects in the parent spectroscopic sample with [ OIII ] + H 𝛽 EWs larger than 2000Å (and spectra sufficiently deep to detect CIV with EW > 10 Å), we find that 40% (5 of 12) show strong CIV emission ( > 10 Å). This suggests that hard ionizing sources capable of powering CIV are likely fairly common in galaxies in a brief (few Myr) window following a strong burst of star formation. By selecting galaxies with very high [OIII]+H 𝛽 EWs (via photometric flux excesses or spectroscopy), it should be possible to build a larger sample of systems with high-ionization lines in the rest-UV.  \nWe compare the CIV EWs of the 𝑧 ≳ 4 galaxies to local star forming systems (e.g., Berg et al. 2016, 2019; Izotov et al. 2016a,b, 2018a,b; Senchyna et al. 2017, 2019; Izotov et al. 2024) with similar H 𝛽 EW in Figure 17. These two quantities provide insight into how the hard radiation field varies with the stellar population age. We see that local metal poor galaxies often power relatively weak CIV emission (1-10 Å) at H 𝛽 EW = 100-300 Å. In contrast, the early galaxies with CIV detections extend to much higher CIV EWs ( ≳ 20-50 Å) than are seen locally at similar H 𝛽 EWs. While several newly-discovered local systems are beginning to approach the intense CIV EWs found in the reionization era (e.g., Izotov et al. 2024), it appears that these very large CIV EWs are not the norm in strong bursts of star formation at low redshift. This may indicate that the stellar populations in reionization era galaxies host harder ionizing sources, or it could hint at some contribution from AGNphotoionization. Alternatively a larger fraction of the resonant CIV emission may be escaping from very high redshift galaxies (see Senchyna et al. 2022), with observed EWs approaching the maximum CIV EW achievable by stellar photoionization models (Plat et al. in prep). Deeper spectra of the CIV emitters should reveal which of these is most likely to be driving the stronger high ionization line emission seen at 𝑧 ≳ 4.  \nFigure 16. [ OIII ] + H 𝛽 EWversus M UV of robust NIV] (black) and CIV (orange) emitters compared to the full sample of spectroscopic (teal) and photometric (grey) measurements (panel a). Panels (b) and (c) provide the distributions of M UV and [ OIII ] + H 𝛽 , respectively for the four samples. In these panels, the points and errorbars above each histogram display the median and standard deviation, respectively, of each sample. The CIV and NIV] emitters appear to have UV luminosities at or above the median of the spectroscopic sample, however nearly all of these objects have EWs elevated significantly above the sample average. The two CIV emitters with [ OIII ] + H 𝛽 EWs that are less than 1000Å have been shown to display evidence of AGN activity (Maiolino et al. 2023b; Scholtz et al. 2023).  \n<!-- image -->', '5.3 The Demographics of NIV] Emitters at 𝑧 ≳ 4': 'Wenowquantify the prevalence of the NIV] emitters in the spectroscopic sample. As described in §5.1, we recover three systems with significant NIV] detections (Table 5) including GN-z11 (Bunker et al. 2023) and CEERS-1019 (Larson et al. 2023; Isobe et al. 2023; Marques-Chaves et al. 2024), as well as one new source (MSA ID: JADES-202208). We display the spectrum of JADES-202208 ( 𝑧 = 5 . 449) in Figure 15. The NIV] emission is detected at S/N=5.0, along with CIV (S/N=5.1), He II (S/N=4.4), and a strong detection of the OIII] auroral line (S/N=4.9). The NIV] EW (5 . 1 + 1 . 6 -2 . 6 Å) and CIV EW (5 . 6 + 2 . 0 -2 . 0 Å) are comparable, while the [OIII]+H 𝛽 EW (3268Å) is among the largest in the parent catalog.  \nThe recovery of only 3 NIV] lines in a redshift catalog of 747 galaxies highlights the challenges in blindly increasing the nitrogen line sample with non-targeted spectroscopy. One reason for the low detection rate is that the nitrogen lines are typically weak (EW ≃ 5 -10 Å), requiring extremely deep rest-UV spectroscopy for meaningful limits in fainter sources. To quantify the fraction of nitrogen line emitters in the 𝑧 > 4 galaxy sample, we only consider spectra where a NIV] line with an EW larger than 5Å can be constrained at > 5 𝜎 . This limiting EW is more strict than the one imposed for CIV in Section 5.2 reflecting the typically-weaker NIV] line. Of the 737 galaxies across the full sample, only 38 are of sufficient depth to reach this limit. Given the 3 NIV] detections, this suggests that 7% of 𝑧 ≳ 4 galaxies with sufficiently deep JWST spectra exhibit NIV] with EW > 5 Å. This fraction is comparably low to that found for strong CIV emitting galaxies, suggesting that the majority of early galaxies are not observed in a phase with very hard ionizing sources and enhanced nitrogen.  \nTo understand what makes nitrogen line emitters rare at 𝑧 ≳ 4, we compare the absolute magnitude and [OIII]+ 𝐻𝛽 EWs of the NIV] emitters to the parent spectroscopic sample in Figure 16. Currently the NIV] detections are primarily limited to the bright galaxies (median M UV = -21 . 1) which is not surprising given the sensitivity to detect NIV]. It is conceivable that deeper spectroscopy (or highlymagnified sources in cluster fields) may reveal strong nitrogen line emission in fainter galaxies. What does appear to be clear is that the NIV] emitters exclusively have very large [ OIII ] + H 𝛽 EWs (2600-  \n4200Å), placing them in the upper 2% of the EW distribution at similar redshifts and indicating extremely young stellar populations (Endsley et al. 2023a; Matthee et al. 2023). Of course, very large [ OIII ] + H 𝛽 EWs are also found in RXCJ2248-ID and A1703-zd6, neither of which is included in the sample described above. If we consider galaxies with high [ OIII ] + H 𝛽 EWs ( > 2000Å) and spectra deep enough to recover NIV] emission at our EW limit, we find that nitrogen emitters comprise a much larger fraction (30%) of the population. Spectroscopic programs targeting the highest equivalent width optical line emitters may provide a more efficient means of identifying additional nitrogen emitters.', '6 DISCUSSION': "Over the past several years, a distinct phase in the early evolution of reionization-era galaxies has been uncovered by deep spectroscopy. The first glimpse of galaxies in or near this phase was via discovery of prominent C /i.pc/v.pc emission detected pre-JWST in a handful of luminous lensed systems, a transition rarely encountered in nebular emission at such strengths in lower-redshift star-forming systems (Stark et al. 2015b; Mainali et al. 2017; Schmidt et al. 2017). More recently, the detection of strong lines of ionized nitrogen including N/i.pc/v.pc] in GN-z11 clearly established that some high-redshift galaxies harbor gas in a remarkable nitrogen-enriched and high-density phase that is not yet fully understood (Bunker et al. 2023; Cameron et al. 2023c; Senchyna et al. 2023; Maiolino et al. 2023a). In Topping et al. (2024b) and this paper, we have shown that both prominent 𝑧 > 6 C/i.pc/v.pc emitters known before JWST launch also reveal signatures of dense, nitrogen-enriched gas, suggesting that these two phenomena are potentially linked.  \nWith improved statistics on the demographics of nitrogen emitters, we can begin to place some constraints on the nature of this observational phase. The clear association of strong C /i.pc/v.pc and N /i.pc/v.pc detections with very intense optical nebular emission (H 𝛽 EW ≳ 400 Å) in a sample of over 700 𝑧 ≳ 4 galaxies establishes that this is likely a short-lived phase during bursts of star formation. The elevated H 𝛽 EWscorrespond to galaxies caught with very large instantaneous sSFR. In the context of bursty star formation histo-  \nFigure 17. CIV EW versus H 𝛽 EW for galaxies in the local universe (blue squares) and at high redshift (red circles). The high-redshift CIV-emitter sample includes objects in the A1703 field in addition to RXCJ2248-ID (Topping et al. 2024b) and the CIV emitters identified in Table 5. The local comparison samples include objects presented by Senchyna et al. (2019), Izotov et al. (2016a,b, 2018a,b), Berg et al. (2016, 2019), Senchyna et al. (2017), and (Izotov et al. 2024). We indicate the average H 𝛽 EWandthe 1 𝜎 distribution width at 𝑧 > 6 as the vertical dashed line and shaded region, respectively, from the photometric sample of Endsley et al. (2023a). The high-redshift galaxies nearly universally have higher CIV EWs than the local samples. For H 𝛽 EWsof 100-300Å, the highest CIV EWs locally only reach ≃ 10Å, while several galaxies at 𝑧 > 4 exceed 30Å. At the highest H 𝛽 EWs ( > 400Å), the local and high-redshift galaxies have comparable CIV EWs, implying they may be reaching the maximum EW that can be produced by photoionization from massive stars.  \n<!-- image -->  \nbe present at 𝑧 ≳ 6 (e.g., Ma et al. 2018; Furlanetto & Mirocha 2022; Dome et al. 2024; Endsley et al. 2023a; Strait et al. 2023), the largest sSFRs will be found in galaxies with the highest amplitude bursts (i.e., the highest star formation rate at fixed stellar mass) observed near the peak of this star formation activity. The H 𝛽 EWs thus depend both on the light-weighted stellar population age and on the relative strength of the burst. The demographic trends we found in §5 may therefore suggest that the nitrogen enrichment is limited to the highest peak intensity bursts of star formation. Specifically, while galaxies with reasonably intense line emission (H 𝛽 EW ≃ 200Å) are likely in the midst of recent upturns of star formation, they appear not to power the nitrogen emission seen in more intense bursts. We illustrate this schematically in Figure 18.  \nThe fact that these nitrogen signatures are only observed at the highest-sSFRs is an important clue as to their nature. It is natural to expect that galaxies undergoing the strongest bursts are likely to have unique ISM conditions. In particular, we expect the cold gas densities to reach extremely high values in these systems (e.g., Kim et al. 2021; Grudić et al. 2021) where the cluster formation efficiency may be elevated (e.g., Adamo et al. 2011, 2015; Ma et al. 2020). Hence the nitrogen line emitters appear to be preferentially found in galaxies that are likely to be forming very dense star cluster complexes, potentially similar to the conditions under which globular clusters formed. This fact also suggests that the abundances are short-lived in the ionized gas phase, as otherwise we would ex-  \nect to find similar evidence for nitrogen over-abundances (in these or lower-ionization nitrogen lines) in galaxies with somewhat less extreme optical equivalent widths, representing the same bursts at somewhatolder ages. Some combination of processes must alter this situation on a short (of-order Myrs; §5.2) timescale: for instance, ongoing star formation involving this gas (as is the case if some of this material is locked-up in stars with globular cluster-like abundance patterns), stellar feedback heating and evacuating this gas, yields from Type II supernovae or time-integrated normal massive star winds polluting this material (which should boost oxygen and carbon relative to nitrogen; e.g. Nomoto et al. 2013; Charbonnel et al. 2023), or inflows of pristine gas diluting these abundances.  \nHowever, high-sSFR alone is clearly not sufficient on its own to produce the conditions responsible for this enrichment. Strong bursts with similar H 𝛽 EWs in the local universe do not appear to show the prominent nitrogen lines being recovered at 𝑧 ≳ 6 (e.g., Berg et al. 2019; Senchyna et al. 2017, 2019; Izotov et al. 2016a,b, 2018a, 2024). This suggests that the conditions required to produce the anomalous nitrogen abundance pattern may be uniquely present in high redshift galaxies, perhaps reflecting different star formation processes. One possibility is that the combination of high gas densities and low metallicities (thought to be common at 𝑧 ≳ 6) maylead to rapid cooling times and less effective feedback, resulting in elevated star formation efficiencies in a subset of early galaxies (Dekel et al. 2023; Li et al. 2023). If present in early galaxies, these so-called 'feedback-free starbursts' (FFBs) would likely be the highest sSFR galaxies at a given epoch, the same population that appears to often show C /i.pc/v.pc and N /i.pc/v.pc detections. We note that the compact star forming complexes associated with the nitrogenenhancements do appear to reach stellar mass surface densities (Topping et al. 2022) where feedback-free starburst conditions may be expected (Grudić et al. 2021). Hence while FFBs have recently been discussed as a means of explaining the excess of galaxies at 𝑧 ≳ 10, they may also play an important role in driving the nitrogenenhanced abundance pattern that is being seen in the strongest bursts of the reionization era and in the fossil record of globular clusters (see also e.g. Renzini 2023).  \nThe ionizing sources responsible for driving the nitrogen enhancement have been discussed at length in the recent literature (e.g., Senchyna et al. 2023; Nandal et al. 2024a; Topping et al. 2024b; Maiolino et al. 2023a; Castellano et al. 2024; MarquesChaves et al. 2024). With improved constraints on the demographics of this population, we can gain new insights into the nature of this population. The fact that this abundance pattern is only found in very strong bursts where densities are large is consistent with a picture where dynamical interactions play an important role in forming the polluting sources (Gieles et al. 2018; Martins et al. 2020; Charbonnel et al. 2023), as pointed out previously (Senchyna et al. 2023). Runaway collisions may alter the stellar initial mass function in these dense stellar systems, creating a population of very massive stars (VMSs; ≃ 10 2 -10 3 M ⊙ ) or supermassive stars (SMSs; ≃ 10 3 -10 5 M ⊙ ). The contribution of both populations to the nitrogen-enhanced abundance pattern has been recently investigated in some detail (e.g., Nagele & Umeda 2023; Charbonnel et al. 2023; Senchyna et al. 2023; Topping et al. 2024b; Nandal et al. 2024a; Marques-Chaves et al. 2024). If low metallicity VMSs altered by interactions and mergers are present in nitrogen emitters, we would expect a hard radiation field, perhaps even more intense than is commonly found in metal poor galaxies locally. This may help explain the intense CIV emission in these systems, with EWs in excess of most local samples (Figure 17). In contrast, the SMSs are not expected to power hard radiation fields (e.g., Gieles et al.  \nFigure 18. Schematic demonstrating the physical scenarios that result in a given rest-frame UV spectrum. The left panels present a UV spectrum that contains high-ionization lines (e.g., He II and CIV) along with strong nitrogen lines (e.g., NIV] and NIII]) such as those seen in A1703-zd6 (This Work) and RXCJ2248-ID (Topping et al. 2024b). In this scenario, the UV lines are driven a significant burst of star formation (sSFR = 300 -1000 Gyr -1 ) which also drives very strong rest-optical emission lines (H 𝛽 EW = 400 -600Å). The middle panels reflect rest-UV spectra such as those of A1703-zd2 and A1703-zd5.2, which display strong high-ionization lines (CIV and He II), but do not show emission from nitrogen. These systems are driven by a moderate burst of star formation (sSFR = 50 -300 Gyr -1 ) and have above-average rest-optical emission line EWs (H 𝛽 EW = 200 -400Å). The right panels demonstrate a scenario where the most massive stars are absent, representing recent star-formation activity corresponding to an sSFR of < 50 Gyr -1 . These systems do not have the ionizing radiation field capable of driving CIV and He II emission, and do not have an ISM enriched in nitrogen, and display H 𝛽 EWs down to ≃ 50Å. The most common line seen in the rest-frame UV for this scenario is CIII].  \n<!-- image -->  \n2018; Martins et al. 2020), so if such a population was responsible for expelling the nitrogen-enriched material, we would likely also need there to be a separate population of hard ionizing sources to power the high ionization lines. Or we may be witnessing a later evolution stage (following the SMS phase) during a second burst of star formation, where massive stars forming from already-enriched material now photoionize their nitrogen-enhanced surroundings. Along these lines, the potentially important role of multiple bursts in explaining the abundance pattern has been discussed in Kobayashi & Ferrara (2024). In multiple burst scenarios, a key prediction is that there should be a population of strong bursts without significant nitrogen emission lines (i.e., the first burst in such scenarios). The demographics presented in §5 are to first order consistent with this possibility, as not all galaxies with H 𝛽 EW ≳ 400 Å support strong nitrogen line emission. More detailed comparison of formation models to these demographics may aid in differentiating between these scenarios.  \nIt has been recently noted (e.g., Marques-Chaves et al. 2024) that the dense conditions in the nitrogen emitters also appear op-  \ntimized for creating intermediate mass black holes (IMBHs). If VMSs or SMSs are common in the star clusters of nitrogen emitters (see discussion above), they may provide an efficient path for the formation of black holes, with subsequent growth from tidal disruption events and black hole mergers (e.g., Quinlan & Shapiro 1990; Portegies Zwart et al. 2004; Rantala et al. 2024). In cases where this picture holds, the nitrogen emitters at 𝑧 ≳ 6 may correspond to an early phase where relatively massive black holes are efficiently formed following intense and dense bursts of star formation. In multiple burst scenarios (e.g. Kobayashi & Ferrara 2024), it is possible that IMBHs could have already assembled in earlier star formation episodes. In this case, the nitrogen emitters may already harbor relatively massive black holes which may be fueled during the current burst of star formation. Here we may expect to see AGN signatures, as is plausibly the case for GNz11 (Maiolino et al. 2023a). Distinguishing between these cases requires confirmation of the presence of VMSs or SMSs in the nitrogen emitters or their progenitors.  \nMany of the outstanding questions about these objects are connected to the uncertain nature of the stellar populations and ionizing  \nsources powering them. An important first step in characterizing such massive star populations was taken in Rivera-Thorsen et al. (2024), using deep JWST spectroscopy of the 𝑧 = 2 . 37 Sunburst Arc, probing the rest-frame optical of this object. The Sunburst Arc is a crucial lower-redshift laboratory for understanding these 𝑧 > 6 spectra, as the rest-UV spectrum of the highly-lensed LyC-leaking cluster in this galaxy reveals similarly prominent N /i.pc/i.pc/i.pc] emission requiring significant nitrogen enrichment: N / O = 0 . 6 (Pascale et al. 2023), albeit at a much smaller mass scale (a lensed single cluster) and likely higher metallicity (plausibly consistent with the lack of N/i.pc/v.pc] emission). The new JWST spectra reveal strong broadened stellar wind emission in both the blue and red 'bumps' typically associated with Wolf-Rayet (WR) stars (Rivera-Thorsen et al. 2024). These bumps are composed of blended lines of (predominantly) ionized nitrogen/helium and carbon, respectively; and their detection supports the presence of a prominent population of luminous and stripped (or otherwise processed to reveal layers exposed to hot H-burning) massive stars driving dense, highly-ionized winds. This may be critical to understanding the source of the nitrogen-enriched material seen in the ionized gas phase in this object. The strong winds of very massive stars ( ∼ 100-1000 M ⊙ ) near the Eddington limit are a natural source of nitrogen-enriched material themselves (e.g. Vink 2023), and are expected to be visible in both optical WR bumps (e.g. Martins & Palacios 2022). However, the wind signatures of these VMSs are highly sensitive to mass loss rates and to detailed evolution and wind modeling, and challenging to disentangle from 'classical' WR stars in integrated light (e.g. Martins et al. 2023). If VMS or WR winds are directly responsible for the pollution observed in these 𝑧 > 6 sources, deep UV spectroscopy promises to readily reveal their continuum signatures even at low metallicities (e.g. Senchyna et al. 2021; Martins et al. 2023); and in any case, detailed comparisons both with lower-redshift samples and with population synthesis predictions including these sources would enable the first characterization of their constituent massive stars and a critical test for the various enrichment scenarios.", '7 SUMMARY': "We present new JWST /NIRSpec observations of galaxies in the Abell 1703 field. Our primary target is A1703-zd6, a bright (H=25.9) 𝑧 = 7 . 0435 galaxy known to power strong CIV emission from ground-based observations (Stark et al. 2015a). The new NIRSpec spectrum spans the rest-UV to optical, constraining the metal content, abundance pattern, and ionizing sources linked with such strong CIV emission. We also present spectra of seven additional bright (H=23.8-25.9) 𝑧 ≃ 6 -8 galaxies in the Abell 1703 field. To place these results in context, we investigate the prevalence of CIV and nitrogen lines in a newly-constructed redshift catalog of 737 𝑧 ≳ 4 galaxies with public NIRSpec spectra.  \nOneofourcentralgoalsis to understand the physical conditions that lead to the nitrogen enhancements seen in recent metal poor galaxies at high redshift (Bunker et al. 2023; Marques-Chaves et al. 2024; Topping et al. 2024b), testing whether there is a link between the presence of hard radiation fields and the anomalous abundance pattern. We summarize our primary findings below.  \n- 1. The rest-UV spectrum of A1703-zd6 confirms the strong CIV emission (EW=19.4 Å) seen with Keck, while also revealing strong emission from a suite of other emission lines (i.e., OIII], CIII], He II). The CIV EW is larger than that seen in most metal poor star forming galaxies found in nearby ( 𝑧 < 0 . 02) galaxies. The  \nline ratios indicate a hard radiation field (CIV/CIII]=3.3), albeit consistent with expectations for stellar photoionization.  \n- 2. The rest-optical spectrum of A1703-zd6 exhibits a suite of strong emission lines. The [OIII]+H 𝛽 EW (4116 Å) is extremely large, corresponding to the upper 0.5% of EWs seen in the reionization-era. The ionized gas is found to be metal poor (12 + log ( O / H ) = 7 . 47) and very high in electron density (8 -19 × 10 4 cm -3 ). The elevated gas densities result in collisional de-excitation of [OII] (boosting the O32 ratio) and strong He I lines (He I 𝜆 5876/H 𝛽 =0.26). The spectral features are similar to those of RXCJ2248-ID (Topping et al. 2024b), suggesting high densities are a feature of HII regions in reionization-era galaxies with the largest [OIII] and H 𝛽 EWs.\n- 3. We find that the emission line spectrum of A1703-zd6 can be explained by a low metallicity young (1.6 Myr) stellar population formed in a recent burst of star formation. The stellar mass (10 7 . 70 M ⊙ ) associated with the burst appears concentrated in a compact (unresolved) region less than 120 pc in size. NIRCam imaging will be required to put a more stringent limit on the size. Nonetheless the current measurement suggests that A1703-zd6 hosts a high density of metal poor massive stars.\n- 4. The rest-UV spectrum of A1703-zd6 reveals strong NIV] and NIII] emission lines, indicating that the ionized gas is both elevated in nitrogen and deficient in oxygen (log ( N / O ) = -0 . 6). The abundance pattern is similar to that seen in GNz11 Bunker et al. (2023) and RXCJ2248-ID (Topping et al. 2024b). The discovery of nitrogen enhancements in both A1703-zd6 and RXCJ2248-ID, two galaxies identified based on strong CIV emission, suggests that there may be a physical link between the presence of hard ionizing sources and the processes leading to the elevated nitrogen content in the ISM. We suggest both may be related to a dense clusters of massive stars formed after a strong burst of star formation.\n- 5. We also find strong CIV emission (EW=10.5 and 31.6 Å) in A1703-zd2 and A1703-zd5.2, two 𝑧 ≳ 6 galaxies with restoptical emission lines suggesting metal poor gas (12 + log ( O / H ) = 7 . 35 ± 0 . 30 and 7 . 26 ± 0 . 28). Both galaxies are also found to lie above the star forming main sequence at 𝑧 ≃ 6 -8 with large specific star formation rates (59 and 32 Gyr -1 assuming a CSFH). The H 𝛽 EWs are very large (209 Å and 170 Å), yet they are ≳ 2 × less than A1703-zd6. Neither source shows evidence for nitrogen enhancements. Doublet ratios suggest electron densities are large (1 -2 × 10 4 cm -3 ) but not as elevated as the sources showing the nitrogen emission. We suggest that the H 𝛽 EW is likely to provide the best indicator of whether a galaxy exhibits nitrogen enhancements seen in a handful of early JWST spectra.\n- 6. To better test this physical picture, we investigate the spectra of 737 𝑧 ≳ 4 galaxies available in our reductions of the public JWST database. Strong CIV ( > 10 Å) and NIV] ( > 5 Å) emission is extremely rare, found in only 8 and 7% of 𝑧 ≳ 4 galaxies, where we only consider those spectra sufficiently deep to reach the respective EW limits. We find that the strong UV lines are primarily seen in galaxies with extremely large [OIII]+H 𝛽 EWs ( > 2500 Å), corresponding to extremely young stellar populations ( < 5 Myr) formed in a recent strong burst of star formation. Such extremely young ages are rare in early star forming systems. Targeted observations of galaxies with these large optical line EWs should yield larger samples of nitrogen-enhanced galaxies with hard ionizing sources.  \nTable 5. Additional CIV and NIV] emitters at 𝑧 ≳ 4 identified from public NIRSpec datasets. We list the JWST program where the NIRSpec observation is taken, their IDs, coordinates, spectroscopic redshifts, absolute UV magnitudes ( 𝑀 UV ), [OIII]+H 𝛽 equivalent widths, the CIV and NIV] equivalent widths, the mode of the NIRSpec observations (i.e., prism or 𝑅 ∼ 1000 medium grating), and literature references when available. Here, the 𝑀 UV for UNCOVER sources in the Abell 2744 cluster have been corrected for gravitational lensing with the Furtak et al. (2023) lensing map. For [OIII]+H 𝛽 EWs, we adopt the values derived directly from the prism spectra when the rest-optical continuum is detected, and values inferred from NIRCam SED fitting in other cases. Emission-line EWs marked with a '*' are designated as tentative.  \n𝑎 Broad-line AGNs reported in Harikane et al. (2023b) and Maiolino et al. (2023b).  \n- [1] Larson et al. (2023); [2] Isobe et al. (2023); [3] Tang et al. 2023; [4] Marques-Chaves et al. (2024); [5] Heintz et al. (2024); [6] Bunker et al. (2023); [7] Harikane et al. (2023b);  \n[8] Maiolino et al. 2023b; [9] Witstok et al. 2024; [10] Tang et al. 2024; [11] Navarro-Carrera et al. (2024); [12] Cameron et al. (2023b); [13] Scholtz et al. 2023; [14] Curti et al. 2024;  \n- [15] Fujimoto et al. 2023b; [16] Schaerer et al. 2024; [17] Chen et al. 2024", 'ACKNOWLEDGEMENTS': 'MWT acknowledges support from the NASA ADAP program through the grant number 80NSSC23K0467. DPS acknowledges support from the National Science Foundation through the grant AST-2109066. AZ acknowledges support by Grant No. 2020750 from the United States-Israel Binational Science Foundation (BSF) and Grant No. 2109066 from the United States National Science Foundation (NSF); by the Ministry of Science & Technology, Israel; and by the Israel Science Foundation Grant No. 864/23.  \nThis work is based in part on observations made with the NASA/ESA/CSA James Webb Space Telescope. The data were obtained from the Mikulski Archive for Space Telescopes at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-03127 for JWST. The authors acknowledge the CEERS, DDT-2750, and UNCOVER teams led by Steven L. Finkelstein, Pablo Arrabal Haro, and I. Labbé & R. Bezanson for developing their observing program with a zero-exclusive-access period. Some of the data products presented herein were retrieved from the Dawn JWST Archive (DJA). DJA is an initiative of the Cosmic Dawn Center (DAWN), which is funded by the Danish National Research Foundation under grant DNRF140.', 'DATA AVAILABILITY': 'The data underlying this article may be presented upon reasonable request to the corresponding author.', 'REFERENCES': 'Abazajian K. N., et al., 2009, ApJS, 182, 543 Adamo A., Östlin G., Zackrisson E., 2011, MNRAS, 417, 1904 Adamo A., Kruĳssen J. M. D., Bastian N., Silva-Villa E., Ryon J., 2015, MNRAS, 452, 246 Adamo A., et al., 2024, arXiv e-prints, p. arXiv:2401.03224 Akerman C. J., Carigi L., Nissen P. E., Pettini M., Asplund M., 2004, A&A, 414, 931 Alexandroff R., et al., 2013, MNRAS, 435, 3306 Allen S. W., et al., 1992, MNRAS, 259, 67 Allen J. T., Hewett P. C., Maddox N., Richards G. T., Belokurov V., 2011, MNRAS, 410, 860 Arellano-Córdova K. Z., et al., 2022a, ApJ, 940, L23 Arellano-Córdova K. Z., et al., 2022b, ApJ, 940, L23 Arrabal Haro P., et al., 2023a, Nature, 622, 707 Arrabal Haro P., et al., 2023b, ApJ, 951, L22 Asplund M., Grevesse N., Sauval A. J., Scott P., 2009, ARA&A, 47, 481 Bastian N., Lardo C., 2018, ARA&A, 56, 83  \nBerg D. A., Skillman E. D., Henry R. B. C., Erb D. K., Carigi L., 2016, ApJ,', '26 M. W. Topping': 'A9  \nTables A1 and A2 provide the full set of emission lines in the rest-frame UV and optical, respectively, that we measured for the spectroscopic sample.', 'APPENDIX B: TENTATIVE CIV AND NIV] EMITTING SOURCES': "In searching for CIV and NIV] emitters in the public JWST spectra, we identified several objects that appeared to display emission in these UV lines, yet we were unable to classify the lines as 'robust'. These tentative sources are presented at the bottom of Table 5. We identify tentative NIV] lines in the spectra of JADES-3990", '28 M. W. Topping': "Table A1. Emission-line fluxes and equivalent widths of rest-frame UV features for our JWST/NIRSpec sample of galaxies in Abell 1703. Upper limits are provided at the 3 𝜎 level. Fluxes and equivalent widths of emission lines lacking wavelength coverage are indicated by a '-'. Equivalent widths are calculated using the continuum level measured from the spectrum if the continuum is detected at > 3 𝜎 .  \nEquivalent Width  \n[  \nÅ  \n]  \nand UNCOVER-34265. Schaerer et al. (2024) recently presented a detection of NIV] from the PRISM spectrum of JADES-3990. While we also note an apparent emission feature at the rest-frame wavelength of NIV], we do not measure it at the 5 𝜎 significance we require to consider it in our robust sample. The spectrum of UNCOVER-34265 presents an excess of emission coincident with NIV]. However in addition to being measured at < 5 𝜎 , the line comprises only one pixel in the spectrum, increasing the likelihood of being either an artifact or the result of scatter from noise. In addition to these tentative NIV] emitters, we identify several objects with tentative CIV emission lines. All of JADES-5173, UNCOVER-10155, andUNCOVER-23604showemissionfeaturesintheirPRISMspectra coincident with the wavelength of CIV. However in each of these cases the S/N of the line does not reach the necessary 5 𝜎 limit. In addition, the tentative CIV lines in JADES-5173 spectra do not appear in its corresponding 𝑅 ∼ 1000 NIRSpec spectra. The apparent CIV emission line in JADES-9903 is detected at a higher S/N than the other tentative sources, however the line centroid is offset from systemic by ∼ 5000 km/s. While this is comparable to the instrumental resolution of the NIRSpec PRISM at this wavelength, this offset is significant enough to preclude JADES-9903 from being included in the robust sample.  \nTable A2. Emission-line fluxes and equivalent widths of rest-frame optical features for our JWST/NIRSpec sample of galaxies in Abell1703. Upper limits are provided at the 3 𝜎 level. Emission lines indicated by a '-' do not have wavelength coverage in the spectra."}