Search

We present a sample of 88 candidate z ∼ 8.5–14.5 galaxies selected from the completed NIRCam imaging from the Cosmic Evolution Early Release Science survey. These data cover ∼90 arcmin ^2 (10 NIRCam pointings) in six broadband imaging filters and one medium-band imaging filter. With this sample we confirm at higher confidence early JWST conclusions that bright galaxies in this epoch are more abundant than predicted by most theoretical models. We construct the rest-frame ultraviolet luminosity functions at z ∼ 9, 11, and 14 and show that the space density of bright ( M _UV = −20) galaxies changes only modestly from z ∼ 14 to z ∼ 9, compared to a steeper increase from z ∼ 8 to z ∼ 4. While our candidates are photometrically selected, spectroscopic follow-up has now confirmed 13 of them, with only one significant interloper, implying that the fidelity of this sample is high. Successfully explaining the evidence for a flatter evolution in the number densities of UV-bright z > 10 galaxies may thus require changes to the dominant physical processes regulating star formation. While our results indicate that significant variations of dust attenuation with redshift are unlikely to be the dominant factor at these high redshifts, they are consistent with predictions from models that naturally have enhanced star formation efficiency and/or stochasticity. An evolving stellar initial mass function could also bring model predictions into better agreement with our results. Deep spectroscopic follow-up of a large sample of early galaxies can distinguish between these competing scenarios.

We present a sample of 1165 extreme emission-line galaxies (EELGs) at 4 < z < 9 selected using James Webb Space Telescope (JWST) NIRCam photometry in the Cosmic Evolution Early Release Science (CEERS) program. We use a simple method to photometrically identify EELGs with H β + [O iii ] (combined) or H α emission of observed-frame equivalent width (EW) > 5000 Å. JWST/NIRSpec spectroscopic observations of a subset (34) of the photometrically selected EELGs validate our selection method: All spectroscopically observed EELGs confirm our photometric identification of extreme emission, including some cases where the spectral-energy-distribution-derived photometric redshifts are incorrect. We find that the medium-band F410M filter in CEERS is particularly efficient at identifying EELGs, both in terms of including emission lines in the filter and in correctly identifying the continuum between H β + [O iii ] and H α in the neighboring broadband filters. We present examples of EELGs that could be incorrectly classified as ultrahigh redshift ( z > 12) as a result of extreme H β + [O iii ] emission blended across the reddest photometric filters. We compare the EELGs to the broader (subextreme) galaxy population in the same redshift range and find that they are consistent with being the bluer, high-EW tail of a broader population of emission-line galaxies. The highest-EW EELGs tend to have more compact emission-line sizes than continuum sizes, suggesting that active galactic nuclei are responsible for at least some of the most extreme EELGs. The photometrically inferred emission-line ratios are consistent with interstellar medium conditions with high ionization and moderately low metallicity, consistent with previous spectroscopic studies.

Low-mass galaxies can significantly contribute to reionization due to their potentially high Lyman continuum (LyC) escape fraction and relatively high space density. We present a constraint on the LyC escape fraction from low-mass galaxies at z = 1.3–3.0. We obtained rest-frame UV continuum imaging with the ACS/SBC and the WFC3/UVIS from the Hubble Space Telescope for eight strongly lensed galaxies that were identified in the Sloan Giant Arc Survey and the Cluster Lensing and Supernova survey with Hubble. The targeted galaxies were selected to be spectroscopically confirmed, highly magnified, and blue in their UV spectral shapes ( β < −1.7). Our targets include intrinsically low-luminosity galaxies down to a magnification-corrected absolute UV magnitude M _UV ∼ −14. We perform custom-defined aperture photometry to place the most reliable upper limits of LyC escape from our sample. From our observations, we report no significant (>2 σ ) detections of LyC fluxes, placing 1 σ upper limits on the absolute LyC escape fractions of 3%–15%. Our observations do not support the expected increased escape fractions of LyC photons from intrinsically UV faint sources. Considering the highly anisotropic geometry of LyC escape, increasing the sample size of faint galaxies in future LyC observations is crucial.

Ultraviolet (UV; rest-frame ∼1200–2000 Å) spectra provide a wealth of diagnostics to characterize fundamental galaxy properties, such as their chemical enrichment, the nature of their stellar populations, and their amount of Lyman-continuum (LyC) radiation. In this work, we leverage publicly released JWST data to construct the rest-frame UV-to-optical composite spectrum of a sample of 63 galaxies at 5.6 < z < 9, spanning a wavelength range from 1500 to 5200 Å. Based on the composite spectrum, we derive an average dust attenuation $E{(B-V)}_{\mathrm{gas}}\,=\,{0.10}_{-0.11}^{+0.10}$ from H β /H γ , an electron density ${n}_{e}={570}_{-290}^{+510}$ cm ^−3 from the [O ii ] doublet ratio, an electron temperature ${T}_{e}={16700}_{-1500}^{+1500}$ K from the [O iii ] λ 4363/[O iii ] λ 5007 ratio, and an ionization parameter $\mathrm{log}(U)=-{2.15}_{-0.03}^{+0.03}$ from the [O iii ]/[O ii ] ratio. Using a direct T _e method, we calculate an oxygen abundance $12+\mathrm{log}({\rm{O}}/{\rm{H}})\,=\,7.67\pm 0.08$ and a carbon-to-oxygen (C/O) abundance ratio $\mathrm{log}({\rm{C}}/{\rm{O}})\,=\,-{0.86}_{-0.10}^{+0.13}$ . This C/O ratio is smaller than compared to z = 0 and z = 2–4 star-forming galaxies, albeit with moderate significance. This indicates the reionization-era galaxies might be undergoing a rapid buildup of stellar mass with high specific star formation rates. A UV diagnostic based on the ratios of C iii ] λ λ 1907, 1909/He ii λ 1640 versus O iii ] λ 1666/He ii λ 1640 suggests that the star formation is the dominant source of ionization, similar to the local extreme dwarf galaxies and z ∼ 2–4 He ii –detected galaxies. The [O iii ]/[O ii ] and C iv /C iii ] ratios of the composite spectrum are marginally larger than the criteria used to select galaxies as LyC leakers, suggesting that some of the galaxies in our sample are strong contributors to the reionizing radiation.

To study the chemical evolution across cosmic epochs, we investigate Ne, S, Cl, and Ar abundance patterns in the Cosmic Origins Spectrograph Legacy Archive Spectroscopic SurveY (CLASSY). CLASSY comprises local star-forming galaxies (SFGs; 0.02 < z < 0.18) with enhanced star formation rates, making them strong analogues to high- z SFGs. With direct measurements of electron temperature, we derive accurate ionic abundances for all elements and assess ionization correction factors (ICFs) to account for unseen ions and derive total abundances. We find Ne/O, S/O, Cl/O, and Ar/O exhibit constant trends with gas-phase metallicity for 12+log(O/H) < 8.5 but significant correlation for Ne/O and Ar/O with metallicity for 12+log(O/H) > 8.5, likely due to ICFs. Thus, the applicability of the ICFs to integrated spectra of galaxies could bias results, underestimating true abundance ratios. Using CLASSY as a local reference, we assess the evolution of Ne/O, S/O, and Ar/O in galaxies at z > 3, finding no cosmic evolution of Ne/O, while the lack of direct abundance determinations for S/O and Ar/O can bias the interpretation of the evolution of these elements. We determine the fundamental metallicity relationship (FMR) for CLASSY and compare to the high-redshift FMR, finding no evolution. Finally, we perform the first mass–neon relationship analysis across cosmic epochs, finding a slight evolution to high Ne at later epochs. The robust abundance patterns of CLASSY galaxies and their broad range of physical properties provide essential benchmarks for interpreting the chemical enrichment of the early galaxies observed with the JWST.

We analyze rest-frame ultraviolet to optical spectra of three z ≃ 7.47–7.75 galaxies whose Ly α emission lines were previously detected with Keck/MOSFIRE observations, using the JWST/NIRSpec observations from the Cosmic Evolution Early Release Science survey. From NIRSpec data, we confirm the systemic redshifts of these Ly α emitters, and emission-line ratio diagnostics indicate these galaxies were highly ionized and metal-poor. We investigate Ly α line properties, including the line flux, velocity offset, and spatial extent. For the one galaxy where we have both NIRSpec and MOSFIRE measurements, we find a significant offset in their flux measurements (∼1.3–5× greater in MOSFIRE) and a marginal difference in the velocity shifts. The simplest interpretation is that the Ly α emission is extended and not entirely encompassed by the NIRSpec slit. The cross-dispersion profiles in NIRSpec reveal that Ly α in one galaxy is significantly more extended than the nonresonant emission lines. We also compute the expected sizes of ionized bubbles that can be generated by the Ly α sources and discuss viable scenarios for the creation of sizable ionized bubbles (>1 physical Mpc). The source with the highest-ionization condition is possibly capable of ionizing its own bubble, while the other two do not appear to be capable of ionizing such a large region, but require additional sources of ionizing photons. Therefore, the fact that we detect Ly α from these galaxies suggests diverse scenarios for the escape of Ly α during the epoch of reionization. High-spectral-resolution spectra with JWST/NIRSpec will be extremely useful for constraining the physics of patchy reionization.

The 3D geometries of high-redshift galaxies remain poorly understood. We build a differentiable Bayesian model and use Hamiltonian Monte Carlo to efficiently and robustly infer the 3D shapes of star-forming galaxies in James Webb Space Telescope Cosmic Evolution Early Release Science observations with $\mathrm{log}{M}_{* }/{M}_{\odot }=9.0\mbox{--}10.5$ at z = 0.5–8.0. We reproduce previous results from the Hubble Space Telescope Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey in a fraction of the computing time and constrain the mean ellipticity, triaxiality, size, and covariances with samples as small as ∼50 galaxies. We find high 3D ellipticities for all mass–redshift bins, suggesting oblate (disky) or prolate (elongated) geometries. We break that degeneracy by constraining the mean triaxiality to be ∼1 for $\mathrm{log}{M}_{* }/{M}_{\odot }=9.0\mbox{--}9.5$ dwarfs at z > 1 (favoring the prolate scenario), with significantly lower triaxialities for higher masses and lower redshifts indicating the emergence of disks. The prolate population traces out a “banana” in the projected $b/a\mbox{--}\mathrm{log}a$ diagram with an excess of low- b / a , large- $\mathrm{log}a$ galaxies. The dwarf prolate fraction rises from ∼25% at z = 0.5–1.0 to ∼50%–80% at z = 3–8. Our results imply a second kind of disk settling from oval (triaxial) to more circular (axisymmetric) shapes with time. We simultaneously constrain the 3D size–mass relation and its dependence on 3D geometry. High-probability prolate and oblate candidates show remarkably similar Sérsic indices ( n ∼ 1), nonparametric morphological properties, and specific star formation rates. Both tend to be visually classified as disks or irregular, but edge-on oblate candidates show more dust attenuation. We discuss selection effects, follow-up prospects, and theoretical implications.

In the current JWST era, rest-frame UV spectra play a crucial role in enhancing our understanding of the interstellar medium (ISM) and stellar properties of the first galaxies in the epoch of reionization ( z > 6). Here, we compare well-known and reliable optical diagrams sensitive to the main ionization source (i.e., star formation, SF; active galactic nuclei, AGN; and shocks) to UV counterparts proposed in the literature—the so-called “UV–BPT diagrams”—using the HST COS Legacy Archive Spectroscopic SurveY (CLASSY), which is the largest high-quality, high-resolution, and broad-wavelength range atlas of far-UV spectra for 45 local star-forming galaxies. In particular, we explore where CLASSY UV line ratios are located in the different UV diagnostic plots, taking into account state-of-the-art photoionization and shock models, and, for the first time, the measured ISM and stellar properties (e.g., gas-phase metallicity, ionization parameter, carbon abundance, and stellar age). We find that the combination of C iii ] λ λ 1907,9 He ii λ 1640 and O iii ] λ 1666 can be a powerful tool to separate between SF, shocks, and AGN at subsolar metallicities. We also confirm that alternative diagrams without O iii ] λ 1666 still allow us to define an SF-locus, with some caveats. Diagrams including C iv λ λ 1548,51 should be taken with caution given the complexity of this doublet profile. Finally, we present a discussion detailing the ISM conditions required to detect UV emission lines, visible only in low gas-phase metallicity (12 + log(O/H) ≲ 8.3) and high ionization parameter (log( U ) ≳ −2.5) environments. Overall, CLASSY and our UV toolkit will be crucial in interpreting the spectra of the earliest galaxies that JWST is currently revealing.

We report on the discovery of two low-luminosity, broad-line active galactic nuclei (AGNs) at z > 5 identified using JWST NIRSpec spectroscopy from the Cosmic Evolution Early Release Science (CEERS) survey. We detect broad H α emission in the spectra of both sources, with FWHM of 2060 ± 290 km s ^−1 and 1800 ± 200 km s ^−1 , resulting in virial black hole (BH) masses that are 1–2 dex below those of existing samples of luminous quasars at z > 5. The first source, CEERS 2782 at z = 5.242, is 2–3 dex fainter than known quasars at similar redshifts and was previously identified as a candidate low-luminosity AGN based on its morphology and rest-frame optical spectral energy distribution (SED). We measure a BH mass of M _BH = (1.3 ± 0.4) × 10 ^7 M _⊙ , confirming that this AGN is powered by the least massive BH known in the Universe at the end of cosmic reionization. The second source, CEERS 746 at z = 5.624, is inferred to be a heavily obscured, broad-line AGN caught in a transition phase between a dust-obscured starburst and an unobscured quasar. We estimate its BH mass to be in the range of M _BH ≃ (0.9–4.7) × 10 ^7 M _⊙ , depending on the level of dust obscuration assumed. We perform SED fitting to derive host stellar masses, M _⋆ , allowing us to place constraints on the BH–galaxy mass relationship in the lowest mass range yet probed in the early Universe. The M _BH / M _⋆ ratio for CEERS 2782, in particular, is consistent with or higher than the empirical relationship seen in massive galaxies at z = 0. We examine the narrow emission line ratios of both sources and find that their location on the BPT and OHNO diagrams is consistent with model predictions for moderately low metallicity AGNs with Z / Z _⊙ ≃ 0.2–0.4. The spectroscopic identification of low-luminosity, broad-line AGNs at z > 5 with M _BH ≃ 10 ^7 M _⊙ demonstrates the capability of JWST to push BH masses closer to the range predicted for the BH seed population and provides a unique opportunity to study the early stages of BH–galaxy assembly.

We report the discovery of an accreting supermassive black hole at z = 8.679. This galaxy, denoted here as CEERS_1019, was previously discovered as a Ly α -break galaxy by Hubble with a Ly α redshift from Keck. As part of the Cosmic Evolution Early Release Science (CEERS) survey, we have observed this source with JWST/NIRSpec, MIRI, NIRCam, and NIRCam/WFSS and uncovered a plethora of emission lines. The H β line is best fit by a narrow plus a broad component, where the latter is measured at 2.5 σ with an FWHM ∼1200 km s ^−1 . We conclude this originates in the broadline region of an active galactic nucleus (AGN). This is supported by the presence of weak high-ionization lines (N V, N IV], and C III]), as well as a spatial point-source component. The implied mass of the black hole (BH) is log ( M _BH / M _⊙ ) = 6.95 ± 0.37, and we estimate that it is accreting at 1.2 ± 0.5 times the Eddington limit. The 1–8 μ m photometric spectral energy distribution shows a continuum dominated by starlight and constrains the host galaxy to be massive (log M/M _⊙ ∼9.5) and highly star-forming (star formation rate, or SFR ∼ 30 M _⊙ yr ^−1 ; log sSFR ∼ − 7.9 yr ^−1 ). The line ratios show that the gas is metal-poor ( Z / Z _⊙ ∼ 0.1), dense ( n _e ∼ 10 ^3 cm ^−3 ), and highly ionized (log U ∼ − 2.1). We use this present highest-redshift AGN discovery to place constraints on BH seeding models and find that a combination of either super-Eddington accretion from stellar seeds or Eddington accretion from very massive BH seeds is required to form this object.

We present JWST/NIRSpec prism spectroscopy of seven galaxies selected from Cosmic Evolution Early Release Science (CEERS) survey NIRCam imaging with photometric redshifts z _phot > 8. We measure emission line redshifts of z = 7.65 and 8.64 for two galaxies. For two other sources without securely detected emission lines we measure $z={9.77}_{-0.29}^{+0.37}$ and ${10.01}_{-0.19}^{+0.14}$ by fitting model spectral templates to the prism data, from which we detect continuum breaks consistent with Ly α opacity from a mostly neutral intergalactic medium. The presence of strong breaks and the absence of strong emission lines give high confidence that these two galaxies have redshifts z > 9.6, but the redshift values derived from the breaks alone have large uncertainties given the low spectral resolution and relatively low S/N of the CEERS NIRSpec prism data. The two z ∼ 10 sources observed are relatively luminous ( M _UV < −20), with blue continua (−2.3 ≲ β ≲ −1.9) and low dust attenuation ( ${A}_{V}\simeq {0.15}_{-0.1}^{+0.3}$ ); and at least one of them has a high stellar mass for a galaxy at that redshift ( $\mathrm{log}({M}_{\star }/{M}_{\odot })\simeq {9.3}_{-0.3}^{+0.2}$ ). Considered together with spectroscopic observations of other CEERS NIRCam-selected high- z galaxy candidates in the literature, we find a high rate of redshift confirmation and low rate of confirmed interlopers (8%). Ten out of 35 z > 8 candidates with CEERS NIRSpec spectroscopy do not have secure redshifts, but the absence of emission lines in their spectra is consistent with redshifts z > 9.6. We find that z > 8 photometric redshifts are generally in agreement (within their uncertainties) with the spectroscopic values, but also that the photometric redshifts tend to be slightly overestimated (〈Δ z 〉 = 0.45 ± 0.11), suggesting that current templates do not fully describe the spectra of very-high- z sources. Overall, the spectroscopy solidifies photometric redshift evidence for a high spatial density of bright galaxies at z > 8 compared to theoretical model predictions, and further disfavors an accelerated decline in the integrated UV luminosity density at z > 8.

We present JWST NIRSpec spectroscopy for 11 galaxy candidates with photometric redshifts of z ≃ 9 − 13 and M _UV ∈ [ −21, −18] newly identified in NIRCam images in the Cosmic Evolution Early Release Science Survey. We confirm emission line redshifts for 7 galaxies at z = 7.762–8.998 using spectra at ∼1–5 μ m either with the NIRSpec prism or its three medium-resolution ( R ∼ 1000) gratings. For z ≃ 9 photometric candidates, we achieve a high confirmation rate of ≃90%, which validates the classical dropout selection from NIRCam photometry. No robust emission lines are identified in three galaxy candidates at z > 10, where the strong [O iii ] and H β lines would be redshifted beyond the wavelength range observed by NIRSpec, and the Ly α continuum break is not detected with the sensitivity of the current data. Compared with Hubble Space Telescope-selected bright galaxies ( M _UV ≃ −22) that are similarly spectroscopically confirmed at z ≃ 8 − 9, these NIRCam-selected galaxies are characterized by lower star formation rates (SFRs; SFR ≃ 4 M _⊙ yr ^−1 ) and lower stellar masses (≃10 ^8 M _⊙ ), but with higher specific SFR (≃40 Gyr ^−1 ), higher [O iii ]+H β equivalent widths (≃1100 Å), and elevated production efficiency of ionizing photons ( $\mathrm{log}({\xi }_{\mathrm{ion}}/\mathrm{Hz}\,{\mathrm{erg}}^{-1})\simeq 25.8$ ) induced by young stellar populations (

Lyman-break galaxy (LBG) candidates at z ≳ 10 are rapidly being identified in James Webb Space Telescope (JWST)/NIRCam observations. Due to the (redshifted) break produced by neutral hydrogen absorption of rest-frame UV photons, these sources are expected to drop out in the bluer filters while being well detected in redder filters. However, here we show that dust-enshrouded star-forming galaxies at lower redshifts ( z ≲ 7) may also mimic the near-infrared (near-IR) colors of z > 10 LBGs, representing potential contaminants in LBG candidate samples. First, we analyze CEERS-DSFG-1, a NIRCam dropout undetected in the F115W and F150W filters but detected at longer wavelengths. Combining the JWST data with (sub)millimeter constraints, including deep NOEMA interferometric observations, we show that this source is a dusty star-forming galaxy (DSFG) at z ≈ 5.1. We also present a tentative 2.6 σ SCUBA-2 detection at 850 μ m around a recently identified z ≈ 16 LBG candidate in the same field and show that, if the emission is real and associated with this candidate, the available photometry is consistent with a z ∼ 5 dusty galaxy with strong nebular emission lines despite its blue near-IR colors. Further observations on this candidate are imperative to mitigate the low confidence of this tentative submillimeter emission and its positional uncertainty. Our analysis shows that robust (sub)millimeter detections of NIRCam dropout galaxies likely imply z ∼ 4–6 redshift solutions, where the observed near-IR break would be the result of a strong rest-frame optical Balmer break combined with high dust attenuation and strong nebular line emission, rather than the rest-frame UV Lyman break. This provides evidence that DSFGs may contaminate searches for ultra-high redshift LBG candidates from JWST observations.

We present a comprehensive analysis of the evolution of the morphological and structural properties of a large sample of galaxies at z = 3–9 using early James Webb Space Telescope (JWST) CEERS NIRCam observations. Our sample consists of 850 galaxies at z > 3 detected in both Hubble Space Telescope (HST)/WFC3 and CEERS JWST/NIRCam images, enabling a comparison of HST and JWST morphologies. We conduct a set of visual classifications, with each galaxy in the sample classified three times. We also measure quantitative morphologies across all NIRCam filters. We find that galaxies at z > 3 have a wide diversity of morphologies. Galaxies with disks make up 60% of galaxies at z = 3, and this fraction drops to ∼30% at z = 6–9, while galaxies with spheroids make up ∼30%–40% across the redshift range, and pure spheroids with no evidence for disks or irregular features make up ∼20%. The fraction of galaxies with irregular features is roughly constant at all redshifts (∼40%–50%), while those that are purely irregular increases from ∼12% to ∼20% at z > 4.5. We note that these are apparent fractions, as many observational effects impact the visibility of morphological features at high redshift. On average, Spheroid-only galaxies have a higher Sérsic index, smaller size, and higher axis ratio than disk or irregular galaxies. Across all redshifts, smaller spheroid and disk galaxies tend to be rounder. Overall, these trends suggest that galaxies with established disks and spheroids exist across the full redshift range of this study, and further work with large samples at higher redshift is needed to quantify when these features first formed.

We report the discovery of a candidate galaxy with a photo- z of z ∼ 12 in the first epoch of the James Webb Space Telescope (JWST) Cosmic Evolution Early Release Science Survey. Following conservative selection criteria, we identify a source with a robust z _phot = ${11.8}_{-0.2}^{+0.3}$ (1 σ uncertainty) with m _F200W = 27.3 and ≳7 σ detections in five filters. The source is not detected at λ < 1.4 μ m in deep imaging from both Hubble Space Telescope (HST) and JWST and has faint ∼3 σ detections in JWST F150W and HST F160W, which signal a Ly α break near the red edge of both filters, implying z ∼ 12. This object (Maisie’s Galaxy) exhibits F115W − F200W > 1.9 mag (2 σ lower limit) with a blue continuum slope, resulting in 99.6% of the photo- z probability distribution function favoring z > 11. All data-quality images show no artifacts at the candidate’s position, and independent analyses consistently find a strong preference for z > 11. Its colors are inconsistent with Galactic stars, and it is resolved ( r _h = 340 ± 14 pc). Maisie’s Galaxy has log M _* / M _⊙ ∼ 8.5 and is highly star-forming (log sSFR ∼ −8.2 yr ^−1 ), with a blue rest-UV color ( β ∼ −2.5) indicating little dust, though not extremely low metallicity. While the presence of this source is in tension with most predictions, it agrees with empirical extrapolations assuming UV luminosity functions that smoothly decline with increasing redshift. Should follow-up spectroscopy validate this redshift, our universe was already aglow with galaxies less than 400 Myr after the Big Bang.

Ly α line profiles are a powerful probe of interstellar medium (ISM) structure, outflow speed, and Lyman-continuum escape fraction. In this paper, we present the Ly α line profiles of the Cosmic Origins Spectrograph (COS) Legacy Archive Spectroscopic SurveY, a sample rich in spectroscopic analogs of reionization-era galaxies. A large fraction of the spectra show a complex profile, consisting of a double-peaked Ly α emission profile in the bottom of a damped, Ly α absorption trough. Such profiles reveal an inhomogeneous ISM. We successfully fit the damped Ly α absorption and the Ly α emission profiles separately, but with complementary covering factors, a surprising result because this approach requires no Ly α exchange between high- N _H _i and low- N _H _i paths. The combined distribution of column densities is qualitatively similar to the bimodal distributions observed in numerical simulations. We find an inverse relation between Ly α peak separation and the [O iii ]/[O ii ] flux ratio, confirming that the covering fraction of Lyman-continuum-thin sightlines increases as the Ly α peak separation decreases. We combine measurements of Ly α peak separation and Ly α red peak asymmetry in a diagnostic diagram, which identifies six Lyman-continuum leakers in the COS Legacy Archive Spectrocopy SurveY (CLASSY) sample. We find a strong correlation between the Ly α trough velocity and the outflow velocity measured from interstellar absorption lines. We argue that greater vignetting of the blueshifted Ly α peak, relative to the redshifted peak, is the source of the well-known discrepancy between shell-model parameters and directly measured outflow properties. The CLASSY sample illustrates how scattering of Ly α photons outside the spectroscopic aperture reshapes Ly α profiles because the distances to these compact starbursts span a large range.

We present Atacama Large Millimeter/submillimeter Array (ALMA) Band 7 observations of a remarkably bright galaxy candidate at ${z}_{\mathrm{phot}}={16.7}_{-0.3}^{+1.9}$ ( M _UV = −21.6), S5-z17-1, identified in James Webb Space Telescope (JWST) Early Release Observation data of Stephen’s Quintet. We do not detect the dust continuum at 866 μ m, ruling out the possibility that S5-z17-1 is a low- z dusty starburst with a star formation rate of ≳30 M _⊙ yr ^−1 . We detect a 5.1 σ line feature at 338.726 ± 0.007 GHz exactly coinciding with the JWST source position, with a 2% likelihood of the signal being spurious. The most likely line identification would be [O iii ]52 μ m at z = 16.01 or [C ii ]158 μ m at z = 4.61, whose line luminosities do not violate the nondetection of the dust continuum in both cases. Together with three other z ≳ 11–13 candidate galaxies recently observed with ALMA, we conduct a joint ALMA and JWST spectral energy distribution (SED) analysis and find that the high- z solution at z ∼ 11–17 is favored in every candidate as a very blue (UV continuum slope of ≃−2.3) and luminous ( M _UV ≃ [ − 24:−21]) system. Still, we find in several candidates that reasonable SED fits (Δ χ ^2 ≲ 4) are reproduced by type II quasar and/or quiescent galaxy templates with strong emission lines at z ∼ 3–5, where such populations predicted from their luminosity functions and EW([O iii ]+H β ) distributions are abundant in survey volumes used for the identification of the z ∼ 11–17 candidates. While these recent ALMA observation results have strengthened the likelihood of the high- z solutions, lower- z possibilities are not completely ruled out in several of the z ∼ 11–17 candidates, indicating the need to consider the relative surface densities of the lower- z contaminants in the ultra-high- z galaxy search.

Spectroscopic studies of extreme-ionization galaxies (EIGs) are critical to our understanding of exotic systems throughout cosmic time. These EIGs exhibit spectral features requiring >54.42 eV photons: the energy needed to ionize helium into He ^2+ fully and emit He ii recombination lines. Spectroscopic studies of EIGs can probe exotic stellar populations or accretion onto intermediate-mass black holes (∼10 ^2 –10 ^5 M _⊙ ), which are the possibly key contributors to the reionization of the Universe. To facilitate the use of EIGs as probes of high-ionization systems, we focus on ratios constructed from several rest-frame UV/optical emission lines: [O iii ] λ 5008, H β , [Ne iii ] λ 3870, [O ii ] λ λ 3727, 3729, and [Ne v ] λ 3427. These lines probe the relative intensity at energies of 35.12, 13.62, 40.96, 13.62, and 97.12 eV, respectively, covering a wider range of ionization than traced by other common rest-frame UV/optical techniques. We use the ratios of these lines ([Ne v ]/[Ne iii ] ≡ Ne53, [O iii ]/H β , and [Ne iii ]/[O ii ]), which are nearby in wavelength, mitigating the effects of dust attenuation and uncertainties in flux calibration. We make predictions from photoionization models constructed from Cloudy that use a broad range of stellar populations and black hole accretion models to explore the sensitivity of these line ratios to changes in the ionizing spectrum. We compare our models to observations from the Hubble Space Telescope and JWST of galaxies with strong high-ionization emission lines at z ∼ 0, z ∼ 2, and 5 < z < 8.5. We show that the Ne53 ratio can separate galaxies with ionization from “normal” stellar populations from those with active galactic nuclei and even “exotic” Population III models. We introduce new selection methods to identify galaxies with photoionization driven by Population III stars or intermediate-mass black hole accretion disks that could be identified in upcoming high-redshift spectroscopic surveys.

The Mg ii λλ 2796, 2803 doublet has been suggested as a useful indirect indicator for the escape of Ly α and Lyman continuum (LyC) photons in local star-forming galaxies. However, studies to date have focused on small samples of galaxies with strong Mg ii or strong LyC emission. Here, we present the first study of Mg ii to probe a large dynamic range of galaxy properties, using newly obtained high-signal-to-noise, moderate-resolution spectra of Mg ii , for a sample of 34 galaxies selected from the Low-redshift Lyman Continuum Survey. We show that the galaxies in our sample have Mg ii profiles ranging from strong emission to P-Cygni profiles to pure absorption. We find that there is a significant trend (with a possibility of spurious correlations of ∼2%) that galaxies that are detected as strong LyC emitters (LCEs) show larger equivalent widths of Mg ii emission, while non-LCEs tend to show evidence of more scattering and absorption features in Mg ii . We then find that Mg ii strongly correlates with Ly α in both equivalent width and escape fraction, regardless of whether emission or absorption dominates the Mg ii profiles. Furthermore, we show that, for galaxies that are categorized as Mg ii emitters, one can use the information of Mg ii , metallicity, and dust to estimate the escape fraction of LyC within a factor of ∼3. These findings confirm that Mg ii lines can be used as a tool for selecting galaxies as LCEs and thus serve as an indirect indicator for the escape of Ly α and LyC.

We present rest-frame optical emission-line flux ratio measurements for five z > 5 galaxies observed by the James Webb Space Telescope Near-Infared Spectrograph (NIRSpec) in the SMACS 0723 Early Release Observations. We add several quality-control and post-processing steps to the NIRSpec pipeline reduction products in order to ensure reliable relative flux calibration of emission lines that are closely separated in wavelength, despite the uncertain absolute spectrophotometry of the current version of the reductions. Compared to z ∼ 3 galaxies in the literature, the z > 5 galaxies have similar [O iii ] λ 5008/H β ratios, similar [O iii ] λ 4364/H γ ratios, and higher (∼0.5 dex) [Ne III ] λ 3870/[O II ] λ 3728 ratios. We compare the observations to MAPPINGS V photoionization models and find that the measured [Ne III ] λ 3870/[O II ] λ 3728, [O iii ] λ 4364/H γ , and [O iii ] λ 5008/H β emission-line ratios are consistent with an interstellar medium (ISM) that has very high ionization ( $\mathrm{log}(Q)\simeq 8-9$ , units of cm s ^−1 ), low metallicity ( Z / Z _⊙ ≲ 0.2), and very high pressure ( $\mathrm{log}(P/k)\simeq 8-9$ , units of cm ^−3 ). The combination of [O iii ] λ 4364/H γ and [O iii ] λ (4960 + 5008)/H β line ratios indicate very high electron temperatures of $4.1\lt \mathrm{log}({T}_{e}/{\rm{K}})\lt 4.4$ , further implying metallicities of Z / Z _⊙ ≲ 0.2 with the application of low-redshift calibrations for “ T _e -based” metallicities. These observations represent a tantalizing new view of the physical conditions of the ISM in galaxies at cosmic dawn.