Search

We provide a review of our current knowledge of galaxies throughout the first billion years of cosmic history. This field has undergone a transformation in the last two years following the launch of $\textit{JWST}$, and we aim to deliver an observational overview of what we have learned about $z\gtrsim 5$ galaxies. We introduce the latest selection methods of high redshift galaxies and describe new measurements of the census of continuum-selected and dusty star forming galaxies at $z\gtrsim 5$. We discuss new measurements of the UV luminosity function at $z\gtrsim 10$ and associated implications for early star formation. We then summarize what is being learned about the physical properties of early galaxies, with up-to-date discussions of the sizes, masses, ages, metallicities, abundance patterns, UV colors, dust properties, and ionizing sources in $z\gtrsim 5$ galaxies. We review observational evidence for bursty star formation histories and describe prospects for characterizing the duty cycle with future observations. We provide a brief overview of the insight being gained through new detections of AGN in early galaxies. Finally we introduce the latest constraints on the contribution of galaxies to reionziation and discuss how $\textit{JWST}$ measurements of Ly$\alpha$ emission offer the potential to probe the earliest stages of the process. This review is meant to provide a broad introduction to those new to the observational study of very high redshift galaxies., Comment: 49 pages, 12 figures; Submitted to Encyclopedia of Astrophysics, November 2024

JWST provides a unique dataset to study reionization's earliest stages, promising insights into the first galaxies. Many JWST/NIRSpec prism spectra reveal smooth Lyman-alpha breaks in z>5 galaxies, implying damping wing scattering by neutral hydrogen. We investigate what current prism spectra imply about the intergalactic medium (IGM), and how best to use NIRSpec spectra to recover IGM properties. We use a sample of 99 z~5.5-13 galaxies with high S/N prism spectra in the public archive, including 12 at z>10. We analyse these spectra using damping wing sightlines from inhomogeneous reionizing IGM simulations, mapping between the distance of a source from neutral IGM and the mean IGM neutral fraction. We marginalise over absorption by local HI around galaxies, and Lyman-alpha emission. We observe a decline in the median and variance of flux around the Lyman-alpha break with increasing redshift, consistent with an increasingly neutral IGM, as ionized regions become smaller and rarer. At $z\gtrsim9$ the spectra become consistent with an almost fully neutral IGM. We find S/N>15 per pixel is required to robustly estimate IGM properties from prism spectra. We fit a sub-sample of high S/N spectra and infer mean neutral fractions $\overline{x}_\mathrm{HI}=0.33^{+0.18}_{-0.27}, 0.64^{+0.17}_{-0.23}$ ($>0.70$ excluding GNz11) at $z \approx 6.5, 9.3$. We also investigate local HI absorption, finding median column density $\log_{10}N_\mathrm{HI}\approx10^{20.8}$ cm$^{-2}$, comparable to $z\sim3$ Lyman-break galaxies, with no significant redshift evolution. We find galaxies showing the strongest absorption are more likely to be in close associations (<500 pkpc), implying enhanced absorption in massive dark matter halos. Future deep prism and grating spectroscopy of z>9 sources will provide tighter constraints on the earliest stages of reionization, key for understanding the onset of star formation., Comment: 15 pages, 13 figures, 11 pages in appendices. Submitted to A&A

The James Webb Space Telescope (JWST) is discovering an increasing number of galaxies well into the early stages of the Epoch of Reionization (EoR). Many of these galaxies are clustered with strong Lyman alpha (Ly$\alpha$) emission, motivating the presence of surrounding cosmic HII regions that would facilitate Ly$\alpha$ transmission through the intergalactic medium (IGM). Detecting these HII "bubbles" would allow us to connect their growth to the properties of the galaxies inside them. Here we develop a new forward-modeling framework to estimate the local HII region size and location from Ly$\alpha$ spectra of galaxy groups in the early stages of the EoR. Our model takes advantage of the complementary information provided by neighboring sightlines through the IGM. Our forward models sample the main sources of uncertainty, including: (i) the global neutral fraction; (ii) EoR morphology; (iii) emergent Ly$\alpha$ emission; and (iv) NIRSpec instrument noise. Depending on the availability of complementary nebular lines, $\sim$ 0.006 $\unicode{x2013}$ 0.01 galaxies per cMpc$^3$, are required to be $\gtrsim$95\% confident that the HII bubble location and size recovered by our method is accurate to within $\sim$ 1 comoving Mpc. This corresponds roughly to tens of galaxies at $z\sim7\unicode{x2013}8$ in $\sim$2x2 tiled pointing with JWST NIRSpec. Such a sample is achievable with a targeted survey with completeness down to $M_{\rm UV}^{\rm min}\lesssim$ -19 $\unicode{x2013}$ -17, depending on the over-density of the field. We test our method on 3D EoR simulations as well as misspecified equivalent width distributions, in both cases accurately recovering the HII region surrounding targeted galaxy groups., Comment: 14 pages, 15 figures, submitted to A&A

Ionized bubble sizes during reionization trace physical properties of the first galaxies. JWST's ability to spectroscopically confirm and measure Lyman-alpha (Ly$\alpha$) emission in sub-L* galaxies opens the door to mapping ionized bubbles in 3D. However, existing Lya-based bubble measurement strategies rely on constraints from single galaxies, which are limited by the large variability in intrinsic Ly$\alpha$ emission. As a first step, we present two bubble size estimation methods using Lya spectroscopy of ensembles of galaxies, enabling us to map ionized structures and marginalize over Ly$\alpha$ emission variability. We test our methods using Gpc-scale reionization simulations of the intergalactic medium (IGM). To map bubbles in the plane of the sky, we develop an edge detection method based on the asymmetry of Ly$\alpha$ transmission as a function of spatial position. To map bubbles along the line-of-sight, we develop an algorithm using the tight relation between Ly$\alpha$ transmission and the line-of-sight distance from galaxies to the nearest neutral IGM patch. Both methods can robustly recover bubbles with radius $\gtrsim$10 comoving Mpc, sufficient for mapping bubbles even in the early phases of reionization, when the IGM is $\sim70-90\%$ neutral. These methods require $\gtrsim$0.002-0.004 galaxies/cMpc$^3$, a $5\sigma$ Ly$\alpha$ equivalent width upper limit of $\lesssim$30\r{A} for the faintest targets, and redshift precision $\Delta z \lesssim 0.015$, feasible with JWST spectroscopy. Shallower observations will provide robust lower limits on bubble sizes. Additional constraints on IGM transmission from Ly$\alpha$ escape fractions and line profiles will further refine these methods, paving the way to our first direct understanding of ionized bubble growth., Comment: 15 pages (+ 3 pages in Appendix), 17 figures, submitted to A&A

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

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

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

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

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

JWST has recently sparked a new era of Ly$\alpha$ spectroscopy, delivering the first measurements of the Ly$\alpha$ escape fraction and velocity profile in typical galaxies at $z\simeq6-10$. These observations offer new prospects for insight into the earliest stages of reionization. But to realize this potential, we need robust models of Ly$\alpha$ properties in galaxies at $z\simeq5-6$ when the IGM is mostly ionized. Here we use new JWST observations from the JADES and FRESCO surveys combined with VLT/MUSE and Keck/DEIMOS data to characterize statistical distributions of Ly$\alpha$ velocity offsets, escape fractions, and EWs in $z\simeq5-6$ galaxies. We find that galaxies with large Ly$\alpha$ escape fractions (> 0.2) are common at $z\simeq5-6$, comprising 30 per cent of Lyman break selected samples. Comparing to literature studies, our census suggests that Ly$\alpha$ becomes more prevalent in the galaxy population toward higher redshift from $z\sim3$ to $z\sim6$, although we find that this evolution slows considerably between $z\sim5$ and $z\sim6$, consistent with modest attenuation from residual HI in the mostly ionized IGM at $z\simeq5-6$. We find significant evolution in Ly$\alpha$ velocity profiles between $z\simeq2-3$ and $z\simeq5-6$, likely reflecting the influence of resonant scattering from residual intergalactic HI on the escape of Ly$\alpha$ emission near line center. This effect will make it challenging to use Ly$\alpha$ peak offsets as a probe of Lyman continuum leakage at $z\simeq5-6$. We use our $z\simeq5-6$ Ly$\alpha$ distributions to make predictions for typical Ly$\alpha$ properties at $z\gtrsim8$ and discuss implications of a recently-discovered Ly$\alpha$ emitter at $z\simeq8.5$ with a small peak velocity offset (156 km s$^{-1}$)., Comment: 31 pages, 12 figures, accepted by MNRAS

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

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

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

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

We describe new JWST/NIRSpec observations of galaxies at $z\gtrsim7$ taken from the CEERS survey. Previous observations of this area have revealed associations of Ly$\alpha$ emitters at redshifts ($z=7.5$, $7.7$, $8.7$) where the intergalactic medium (IGM) is thought to be mostly neutral, leading to suggestions that these systems are situated in large ionized bubbles. We identify 21 $z\gtrsim7$ galaxies with robust redshifts in the CEERS dataset, including 10 in the Ly$\alpha$ associations. Their spectra are indicative of very highly ionized and metal poor gas, with line ratios (O32 $=17.84$ and Ne3O2 $=0.89$, linear scale) and metallicity ($12+\log{(\rm{O/H})}=7.84$) that are rarely seen at lower redshifts. We find that the most extreme spectral properties are found in the six $z\gtrsim7$ Ly$\alpha$ emitters in the sample. Each has a hard ionizing spectrum indicating that their visibility is likely enhanced by efficient ionizing photon production. Ly$\alpha$ velocity offsets are found to be very large ($\gtrsim300$ km s$^{-1}$), likely also contributing to their detectability. We find that Ly$\alpha$ in $z\gtrsim7$ galaxies is $6-12\times$ weaker than in lower redshift samples with matched rest-optical spectral properties. If the bubbles around the Ly$\alpha$ emitters are relatively small ($\lesssim0.5-1$ pMpc), we may expect such significant attenuation of Ly$\alpha$ in these ionized regions. We discuss several other effects that may contribute to weaker Ly$\alpha$ emission at $z\gtrsim7$. Deep spectroscopy of fainter galaxies in the vicinity of the Ly$\alpha$ emitters will better characterize the physical scale of the ionized bubbles in this field., Comment: 31 pages, 22 figures, accepted by MNRAS, references updated

MACS0647$-$JD is a triply-lensed $z\sim11$ galaxy originally discovered with the Hubble Space Telescope. Here we report new JWST imaging, which clearly resolves MACS0647$-$JD as having two components that are either merging galaxies or stellar complexes within a single galaxy. Both are very small, with stellar masses $\sim10^8\,M_\odot$ and radii $r<100\,\rm pc$. The brighter larger component "A" is intrinsically very blue ($\beta\sim-2.6$), likely due to very recent star formation and no dust, and is spatially extended with an effective radius $\sim70\,\rm pc$. The smaller component "B" appears redder ($\beta\sim-2$), likely because it is older ($100-200\,\rm Myr$) with mild dust extinction ($A_V\sim0.1\,\rm mag$), and a smaller radius $\sim20\,\rm pc$. We identify galaxies with similar colors in a high-redshift simulation, finding their star formation histories to be out of phase. With an estimated stellar mass ratio of roughly 2:1 and physical projected separation $\sim400\,\rm pc$, we may be witnessing a galaxy merger 400 million years after the Big Bang. We also identify a candidate companion galaxy C $\sim3\,{\rm kpc}$ away, likely destined to merge with galaxies A and B. The combined light from galaxies A+B is magnified by factors of $\sim$8, 5, and 2 in three lensed images JD1, 2, and 3 with F356W fluxes $\sim322$, $203$, $86\,\rm nJy$ (AB mag 25.1, 25.6, 26.6). MACS0647$-$JD is significantly brighter than other galaxies recently discovered at similar redshifts with JWST. Without magnification, it would have AB mag 27.3 ($M_{UV}=-20.4$). With a high confidence level, we obtain a photometric redshift of $z=10.6\pm0.3$ based on photometry measured in 6 NIRCam filters spanning $1-5\rm\mu m$, out to $4300\,\r{A}$ rest-frame. JWST NIRSpec observations planned for January 2023 will deliver a spectroscopic redshift and a more detailed study of the physical properties of MACS0647$-$JD., Comment: 27 pages, 14 figures, accepted for publication in ApJL

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

As the James Webb Space Telescope approaches scientific operation, there is much interest in exploring the redshift range beyond that accessible with Hubble Space Telescope imaging. Currently, the only means to gauge the presence of such early galaxies is to age-date the stellar population of systems in the reionisation era. As a significant fraction of $z\simeq7-8$ galaxies are inferred from Spitzer photometry to have extremely intense [OIII] emission lines, it is commonly believed these are genuinely young systems that formed at redshifts $z<10$, consistent with a claimed rapid rise in the star formation density at that time. Here we study a spectroscopically-confirmed sample of extreme [OIII] emitters at $z=1.3-3.7$, using both dynamical masses estimated from [OIII] line widths and rest-frame UV to near-infrared photometry to illustrate the dangers of assuming such systems are genuinely young. For the most extreme of our intermediate redshift line emitters, we find dynamical masses $10-100$ times that associated with a young stellar population mass, which are difficult to explain solely by the presence of additional dark matter or gaseous reservoirs. Adopting nonparametric star formation histories, we show how the near-infrared photometry of a subset of our sample reveals an underlying old ($>100$ Myr) population whose stellar mass is $\simeq40$ times that associated with the starburst responsible for the extreme line emission. Without adequate rest-frame near-infrared photometry we argue it may be premature to conclude that extreme line emitters in the reionisation era are low mass systems that formed at redshifts below $z\simeq10$., Comment: 14 pages, 7 figure, accepted by MNRAS

We determine the [OIII]$\lambda5007$ equivalent width (EW) distribution of $1.700<\rm{z}<2.274$ rest-frame UV-selected (M$_{\rm{UV}}<-19$) star-forming galaxies in the GOODS North and South fields. We make use of deep HDUV broadband photometry catalogues for selection and 3D-HST WFC3/IR grism spectra for measurement of line properties. The [OIII]$\lambda5007$ EW distribution allows us to measure the abundance of extreme emission line galaxies (EELGs) within this population. We model a log-normal distribution to the [OIII]$\lambda5007$ rest-frame equivalent widths of galaxies in our sample, with location parameter $\mu=4.24\pm0.07$ and variance parameter $\sigma= 1.33\pm0.06$. This EW distribution has a mean [OIII]$\lambda5007$ EW of 168$\pm1\r{A}$. The fractions of $\rm{z}\sim2$ rest-UV-selected galaxies with [OIII]$\lambda5007$ EWs greater than $500, 750$ and $1000\r{A}$ are measured to be $6.8^{+1.0}_{-0.9}\%$, $3.6^{+0.7}_{-0.6}\%$, and $2.2^{+0.5}_{-0.4}\%$ respectively. The EELG fractions do not vary strongly with UV luminosity in the range ($-21.6

Reionisation-era galaxies often display intense nebular emission lines, both in rest-frame optical ([OIII]+H$\beta$) and ultraviolet (UV; CIII], CIV). How such strong nebular emission is powered remains unclear, with both active galactic nuclei (AGN) and hot stars considered equally viable. The UV continuum slopes of these early systems tend to be very blue ($\beta<-2$), reflecting minimal dust obscuration, young ages, and low metallicities. This contrasts with narrow-lined AGN at $z\sim2-3$, whose UV slopes are significantly redder ($\beta>-1$) than typical star-forming systems in the reionisation era. To investigate the properties of AGN in the reionisation era, we have conducted a search for potential examples of rare analogues with blue continua at intermediate redshift ($z\sim2-3$). Our goals are to determine whether AGN with intense line emission and blue continua exist and thereby to establish the range of rest-frame UV and optical line ratios in this population. In this paper we report the detection of a X-ray luminous AGN at $z=3.21$ (UDS-24561) with extreme [OIII]+H$\beta$ line emission (EW $=1300$ \r{A}) and a blue UV continuum slope ($\beta=-2.34$). MMT/Binospec and Keck/MOSFIRE spectra indicate rest-frame UV line ratios consistent with AGN photoionisation models and rest-frame optical lines with both a narrow component (FWHM $=154$ km$/$s) and extended broad wings (FWHM $=977$ km$/$s), consistent with outflowing gas. We describe how such objects can be identified in future JWST emission line surveys in the reionisation era, thereby providing a valuable census of AGN activity at $z>6$ and understanding their contribution to cosmic reionisation., Comment: 11 pages, 7 figures, accepted by MNRAS

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

Spectroscopic observations of massive $z>7$ galaxies selected to have extremely large [OIII]+H$\beta$ equivalent width (EW $\sim1500$ \r{A}) have recently revealed large Ly$\alpha$ detection rates, in contrast to the weak emission seen in the general population. Why these systems are uniquely visible in Ly$\alpha$ at redshifts where the IGM is likely significantly neutral is not clear. With the goal of better understanding these results, we have begun a campaign with MMT and Magellan to measure Ly$\alpha$ in galaxies with similar [OIII]+H$\beta$ EWs at $z\simeq2-3$. At these redshifts, the IGM is highly ionized, allowing us to clearly disentangle how the Ly$\alpha$ properties depend on the [OIII]+H$\beta$ EW. Here we present Ly$\alpha$ EWs of $49$ galaxies at $z=2.2-3.7$ with intense [OIII]+H$\beta$ line emission (EW $=300-3000$ \r{A}). Our results demonstrate that strong Ly$\alpha$ emission (EW $>20$ \r{A}) becomes more common in galaxies with larger [OIII]+H$\beta$ EW, reflecting a combination of increasingly efficient ionizing photon production and enhanced transmission of Ly$\alpha$. Among the galaxies with the most extreme [OIII]+H$\beta$ emission (EW $\sim1500$ \r{A}), we find that strong Ly$\alpha$ emission is not ubiquitous, with only $50$ per cent of our population showing Ly$\alpha$ EW $>20$ \r{A}. Our data suggest that the range of Ly$\alpha$ strengths is related to the observed ellipticity, with those systems that appear edge-on or elongated having weaker Ly$\alpha$ emission. We use these results to interpret the anomalous Ly$\alpha$ properties seen in $z>7$ galaxies with extreme [OIII]+H$\beta$ emission and discuss implications for the escape of ionizing radiation from these extreme line emitting galaxies., Comment: 13 pages, 8 figure, accepted by MNRAS

Escaping Lyman continuum photons from galaxies likely reionized the intergalactic medium at redshifts $z\gtrsim6$. However, the Lyman continuum is not directly observable at these redshifts and secondary indicators of Lyman continuum escape must be used to estimate the budget of ionizing photons. Observationally, at redshifts $z\sim2-3$ where the Lyman continuum is observationally accessible, surveys have established that many objects that show appreciable Lyman continuum escape fractions $f_{esc}$ also show enhanced [OIII]/[OII] (O$_{32}$) emission line ratios. Here, we use radiative transfer analyses of cosmological zoom-in simulations of galaxy formation to study the physical connection between $f_{esc}$ and O$_{32}$. Like the observations, we find that the largest $f_{esc}$ values occur at elevated O$_{32}\sim3-10$ and that the combination of high $f_{esc}$ and low O$_{32}$ is extremely rare. While high $f_{esc}$ and O$_{32}$ often are observable concurrently, the timescales of the physical origin for the processes are very different. Large O$_{32}$ values fluctuate on short ($\sim$1 Myr) timescales during the Wolf-Rayet-powered phase after the formation of star clusters, while channels of low absorption are established over tens of megayears by collections of supernovae. We find that while there is no direct causal relation between $f_{esc}$ and O$_{32}$, high $f_{esc}$ most often occurs after continuous input from star formation-related feedback events that have corresponding excursions to large O$_{32}$ emission. These calculations are in agreement with interpretations of observations that large $f_{esc}$ tends to occur when O$_{32}$ is large, but large O$_{32}$ does not necessarily imply efficient Lyman continuum escape., Comment: 12 pages, 5 figures, submitted to ApJL on September 15, 2020

Deep spectroscopy of galaxies in the reionization-era has revealed intense CIII] and CIV line emission (EW $>15-20$ \r{A}). In order to interpret the nebular emission emerging at $z>6$, we have begun targeting rest-frame UV emission lines in galaxies with large specific star formation rates (sSFRs) at $1.36$ only in large sSFR galaxies with [OIII]+H$\beta$ EW $>1500$ \r{A}. In contrast to previous studies, we find that many galaxies with intense [OIII] have weak CIII] emission (EW $=5-8$ \r{A}), suggesting that the radiation field associated with young stellar populations is not sufficient to power strong CIII]. Photoionization models demonstrate that the spread in CIII] among systems with large sSFRs ([OIII]+H$\beta$ EW $>1500$ \r{A}) is driven by variations in metallicity, a result of the extreme sensitivity of CIII] to electron temperature. We find that the strong CIII] emission seen at $z>6$ (EW $>15$ \r{A}) requires metal poor gas ($\simeq0.1\ Z_\odot$) whereas the weaker CIII] emission in our sample tends to be found at moderate metallicities ($\simeq0.3\ Z_\odot$). The luminosity distribution of the CIII] emitters in our $z\simeq1-3$ sample presents a consistent picture, with stronger emission generally linked to low luminosity systems ($M_{\rm{UV}}>-19.5$) where low metallicities are more likely. We quantify the fraction of strong CIII] and CIV emitters at $z\simeq1-3$, providing a baseline for comparison against $z>6$ samples. We suggest that the first UV line detections at $z>6$ can be explained if a significant fraction of the early galaxy population is found at large sSFR ($>200$ Gyr$^{-1}$) and low metallicity ($<0.1\ Z_\odot$)., Comment: 21 pages, 10 figures, accepted by MNRAS

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

Recent observations have revealed the presence of strong CIII] emission (EW$_{\rm{CIII]}}>20$ \r{A}) in $z>6$ galaxies, the origin of which remains unclear. In an effort to understand the nature of these line emitters, we have initiated a survey targeting CIII] emission in gravitationally-lensed reionization era analogs identified in HST imaging of clusters from the RELICS survey. Here we report initial results on four galaxies selected to have low stellar masses (2-8$\times$10$^7$ M$_\odot$) and J$_{125}$-band flux excesses indicative of intense [OIII]+H$\beta$ emission (EW$_{\rm{[OIII]+H\beta}}$=500-2000 \r{A}), similar to what has been observed at $z>6$. We detect CIII] emission in three of the four sources, with the CIII] EW reaching values seen in the reionization era (EW$_{\rm{CIII]}}\simeq 17-22$ \r{A}) in the two sources with the strongest optical line emission (EW$_{\rm{[OIII]+H\beta}}\simeq 2000$ \r{A}). We have obtained a Magellan/FIRE near-infrared spectrum of the strongest CIII] emitter in our sample, revealing gas that is both metal poor and highly ionized. Using photoionization models, we are able to simultaneously reproduce the intense CIII] and optical line emission for extremely young (2-3 Myr) and metal poor (0.06-0.08 Z$_\odot$) stellar populations, as would be expected after a substantial upturn in the SFR of a low mass galaxy. The sources in this survey are among the first for which CIII] has been used as the primary means of redshift confirmation. We suggest that it should be possible to extend this approach to $z>6$ with current facilities, using CIII] to measure redshifts of objects with IRAC excesses indicating EW$_{\rm{[OIII]+H\beta}}\simeq 2000$ \r{A}, providing a method of spectroscopic confirmation independent of Ly$\alpha$., Comment: Accepted for publication in MNRAS

We discuss the rest-frame optical emission line spectra of a large (~50) sample of z=3.1 Lyman alpha emitters (LAEs) whose physical properties suggest such sources are promising analogs of galaxies in the reionization era. Reliable Lyman continuum escape fractions have now been determined for a large sample of such LAEs from the Lyman Continuum Escape Survey (LACES) undertaken via deep HST imaging in the SSA22 survey area reported in Fletcher et al. (2019). Using new measures of [OII] emission secured from Keck MOSFIRE spectra we re-examine, for a larger sample, earlier claims that Lyman continuum leakages may correlate with the nebular emission line ratio [OIII]/[OII] as expected for density-bound HII regions. We find that a large [OIII]/[OII] line ratio is indeed a necessary condition for Lyman continuum leakage, strengthening earlier claims made using smaller samples at various redshifts. However, not all LAEs with large [OIII]/[OII] line ratios are leakers and leaking radiation appears not to be associated with differences in other spectral diagnostics. This suggests the detection of leaking radiation is modulated by an additional property, possibly the viewing angle for porous HII regions. We discuss our new results in the context of the striking bimodality of LAE leakers and non-leakers found in the LACES program and the implications for the sources of cosmic reionization., Comment: 15 pages, 6 figures. Accepted for publication in ApJ

Galaxies in the reionization-era have been shown to have prominent [OIII]+H$\beta$ emission. Little is known about the gas conditions and radiation field of this population, making it challenging to interpret the spectra emerging at $z\gtrsim6$. Motivated by this shortcoming, we have initiated a large MMT spectroscopic survey identifying rest-frame optical emission lines in 227 intense [OIII] emitting galaxies at $1.3

We present results from the LymAn Continuum Escape Survey (LACES), a Hubble Space Telescope (HST) program designed to characterize the ionizing radiation emerging from a sample of Lyman alpha emitting galaxies at redshift $z\simeq 3.1$. As many show intense [O III] emission characteristic of $z>6.5$ star-forming galaxies, they may represent valuable low redshift analogs of galaxies in the reionization era. Using HST Wide Field Camera 3 / UVIS $F336W$ to image Lyman continuum emission, we investigate the escape fraction of ionizing photons in this sample. For 61 sources, of which 77% are spectroscopically confirmed and 53 have measures of [O III] emission, we detect Lyman continuum leakage in 20%, a rate significantly higher than is seen in individual continuum-selected Lyman break galaxies. We estimate there is a 98% probability that $\leq 2$ of our detections could be affected by foreground contamination. Fitting multi-band spectral energy distributions (SEDs) to take account of the varying stellar populations, dust extinctions and metallicities, we derive individual Lyman continuum escape fractions corrected for foreground intergalactic absorption. We find escape fractions of 15 to 60% for individual objects, and infer an average 20% escape fraction by fitting composite SEDs for our detected samples. Surprisingly however, even a deep stack of those sources with no individual $F336W$ detections provides a stringent upper limit on the average escape fraction of less than 0.5%. We examine various correlations with source properties and discuss the implications in the context of the popular picture that cosmic reionization is driven by such compact, low metallicity star-forming galaxies., Comment: 26 pages, 16 figures, accepted for publication in ApJ

Recent studies have revealed intense UV metal emission lines in a modest sample of z>7 Lyman-alpha emitters, indicating a hard ionizing spectrum is present. If such high ionization features are shown to be common, it may indicate that extreme radiation fields play a role in regulating the visibility of Lyman-alpha in the reionization era. Here we present deep near-infrared spectra of seven galaxies with Lyman-alpha emission at 5.47 photometric targets. In nine sources we do not detect UV metal lines. However in the z=8.683 galaxy EGSY8p7, we detect a 4.6 sigma emission line in the narrow spectral window expected for NV 1243. The feature is unresolved (FWHM<90 km/s) and is likely nebular in origin. A deep H-band spectrum of EGSY8p7 reveals non-detections of CIV, He II, and OIII]. The presence of NV requires a substantial flux of photons above 77 eV, pointing to a hard ionizing spectrum powered by an AGN or fast radiative shocks. Regardless of its origin, the intense radiation field of EGSY8p7 may aid the transmission of Lyman-alpha through what is likely a partially neutral IGM. With this new detection, five of thirteen known Lyman-alpha emitters at z>7 have now been shown to have intense UV line emission, suggesting that extreme radiation fields are commonplace among the Lyman-alpha population. Future observations with JWST will eventually clarify the origin of these features and explain their role in the visibility of Lyman-alpha in the reionization era., Comment: 15 pages, 3 tables and 7 figures. Accepted for publication in MNRAS

We discuss new Keck/MOSFIRE spectroscopic observations of four luminous galaxies at z~7-9 selected to have intense optical line emission by Roberts-Borsani et al. (2016). Previous follow-up has revealed Lyman-alpha in two of the four galaxies. Our new MOSFIRE observations confirm that Lyman-alpha is present in the entire sample. We detect Lyman-alpha emission in COS-zs7-1, confirming its redshift as z=7.154, and we detect Lyman-alpha in EGS-zs8-2 at z=7.477, verifying a tentative detection presented in an earlier study. The ubiquity of Lyman-alpha in this sample is puzzling given that the IGM is likely significantly neutral over 77 is expected to be strongly luminosity-dependent, with the most effective transmission occurring in systems with intense star formation., Comment: Submitted to MNRAS, 13 pages, 8 figures

We present Keck/MOSFIRE observations of UV metal lines in four bright gravitationally-lensed z~6-8 galaxies behind the cluster Abell 1703. The spectrum of A1703-zd6, a highly-magnified star forming galaxy with a Lyman-alpha redshift of z=7.045, reveals a confident detection of the nebular CIV emission line (unresolved with FWHM < 125 km/s). UV metal lines are not detected in the three other galaxies. At z~2-3, nebular CIV emission is observed in just 1% of UV-selected galaxies. The presence of strong CIV emission in one of the small sample of galaxies targeted in this paper may indicate hard ionizing spectra are more common at z~7. The total estimated equivalent width of the CIV doublet (38 A) and CIV/Lyman-alpha flux ratio (0.3) are comparable to measurements of narrow-lined AGNs. Photoionization models show that the nebular CIV line can also be reproduced by a young stellar population, with very hot metal poor stars dominating the photon flux responsible for triply ionizing carbon. Regardless of the origin of the CIV, we show that the ionizing spectrum of A1703-zd6 is different from that of typical galaxies at z~2, producing more H ionizing photons per unit 1500A luminosity and a larger flux density at 30-50 eV. If such extreme radiation fields are typical in UV-selected systems at z>7, it would indicate that reionization-era galaxies are more efficient ionizing agents than previously thought. Alternatively, we suggest that the small sample of Lyman-alpha emitters at z>7 may trace a rare population with intense radiation fields capable of ionizing their surrounding hydrogen distribution. Additional constraints on high ionization emission lines in galaxies with and without Lyman-alpha detections will help clarify whether hard ionizing spectra are common in the reionization era., Comment: 11 pages, 6 figures, submitted to MNRAS