The JWST-PRIMAL Legacy Survey
Abstract
Context. One of the surprising early findings with JWST has been the discovery of a strong “roll-over” or a softening of the absorption edge of Ly in a large number of galaxies at , in addition to systematic offsets from photometric redshift estimates and fundamental galaxy scaling relations. This has been interpreted as strong cumulative damped Ly absorption (DLA) wings from high column densities of neutral atomic hydrogen (H i), signifying major gas accretion events in the formation of these galaxies.
Aims. To explore this new phenomenon systematically, we assemble the JWST/NIRSpec PRImordial gas Mass AssembLy (PRIMAL) legacy survey of 494 galaxies at , designed to study the physical properties and gas in and around galaxies during the reionization epoch.
Methods. We characterize this benchmark sample in full and spectroscopically derive the galaxy redshifts, metallicities, star-formation rates, and ultraviolet slopes. We define a new diagnostic, the Ly damping parameter to measure and quantify the net effect of Ly emission strength, H i fraction in the IGM, or local H i column density for each source. The JWST-PRIMAL survey is based on the spectroscopic DAWN JWST Archive (DJA-Spec). We describe DJA-Spec in this paper, detailing the reduction methods, the post-processing steps, and basic analysis tools. All the software, reduced spectra, and spectroscopically derived quantities and catalogs are made publicly available in dedicated repositories.
Results. We find that the fraction of galaxies showing strong integrated DLAs with cm-2 only increases slightly from at up to at . Similarly, the prevalence and prominence of Ly emission is found to increase with decreasing redshift, in qualitative agreement with previous observational results. Strong Ly emitters (LAEs) are predominantly found to be associated with low-metallicity and UV faint galaxies. By contrast, strong DLAs are observed in galaxies with a variety of intrinsic physical properties, but predominantly at high redshifts and low metallicities.
Conclusions. Our results indicate that strong DLAs likely reflect a particular early assembly phase of reionization-era galaxies, at which point they are largely dominated by pristine H i gas accretion. At , this gas gradually cools and forms into stars that ionize their local surroundings, forming large ionized bubbles and produce strong observed Ly emission at .
Key Words.:
Galaxies: formation, evolution, high-redshift – Reionization – Intergalactic medium – Interstellar medium – Spectroscopy1 Introduction
The first epoch of galaxy formation is believed to have occurred at (Hashimoto et al., 2018; Robertson, 2022), some million years after the Big Bang. This process was primarily driven by the accretion and gravitational collapse of primordial, neutral gas onto dark matter halos (White & Rees, 1978; Blumenthal et al., 1984; White & Frenk, 1991; Kereš et al., 2005; Schaye et al., 2010; Dayal & Ferrara, 2018), which eventually led to the formation of stars in galaxies. The strong UV radiation originating from the first stars and supermassive black holes gradually ionized their immediate and then large-scale surroundings, initiating the reionization epoch, currently estimated to have ended by (e.g. Stark, 2016; Dayal & Ferrara, 2018; Keating et al., 2020; Robertson, 2022; Bosman et al., 2022; Fan et al., 2023). When these first stars ended their lives as supernovae, the pristine gas was chemically enriched with elements formed through explosive nucleosynthesis (Hoyle, 1954; Cameron, 1957; Woosley & Weaver, 1995). Charting these independent components, their interplay through the “baryon cycle” (Tumlinson et al., 2017; Péroux & Howk, 2020), and their evolution with cosmic time is at the heart of contemporary galaxy formation and evolution studies.
The launch of the James Webb Space Telescope (JWST) has now enabled us to peer deep into this early cosmic epoch, probing the rest-frame optical and ultraviolet (UV) emission of galaxies, potentially all the way back to redshifts (Castellano et al., 2022; Naidu et al., 2022; Harikane et al., 2023; Atek et al., 2023a; Donnan et al., 2023; Bouwens et al., 2023; Wang et al., 2023b; Austin et al., 2023). This leap in near-infrared (NIR) capability now enables detailed characterization of the physical properties and baryonic components of galaxies, reaching the first epoch of galaxy formation, the cosmic dawn. This epoch seems to contain a larger population of more luminous () and more massive galaxies than expected from theoretical models (Labbé et al., 2023; Finkelstein et al., 2023a; Franco et al., 2023; Harikane et al., 2023; Adams et al., 2023; Casey et al., 2023; Bouwens et al., 2023; Chemerynska et al., 2023; McLeod et al., 2024). Various possibilities have been considered to resolve this discrepancy, including a varying initial mass function (IMF; Haslbauer et al., 2022; Trinca et al., 2023, though see Rasmussen Cueto et al. 2023), radiation pressure pushing dust away from star-forming regions (Ferrara et al., 2023), bursty star-formation histories (Sun et al., 2023), or an unexpected overabundance of active galactic nuclei (AGN; Inayoshi et al., 2022; Pacucci et al., 2022; Dayal et al., 2024). The apparent discrepancy may be exaggerated by selection biases towards the youngest, most highly star-forming galaxies with more bursty star formation histories at (Mason et al., 2023; Mirocha & Furlanetto, 2023; Shen et al., 2023), making the galaxies appear temporarily brighter in the UV. As neutral gas accretion is the primary driver and fuel for star formation, a robust measure of this component is key to resolving the debate.
Recent JWST spectroscopic studies have pushed the limit of the most distant galaxies confirmed to –13 (Curtis-Lake et al., 2023; Wang et al., 2023a; Fujimoto et al., 2023b; Zavala et al., 2024; Castellano et al., 2024). The metal content of galaxies during the reionization epoch can now also be readily measured based on direct -based methods for sufficiently bright galaxies (e.g., Schaerer et al., 2022; Arellano-Córdova et al., 2022; Taylor et al., 2022; Brinchmann, 2023; Curti et al., 2023a; Katz et al., 2023; Rhoads et al., 2023; Trump et al., 2023; Heintz et al., 2023c; Nakajima et al., 2023; Sanders et al., 2024) or inferred through rest-frame optical strong-line calibrations (Heintz et al., 2023a; Langeroodi et al., 2023; Curti et al., 2023b; Matthee et al., 2023). These observations show a substantial offset from the fundamental-metallicity relation (FMR, Mannucci et al., 2010; Lara-López et al., 2010; Curti et al., 2020; Sanders et al., 2021) established at by up to 0.5 dex at (Heintz et al., 2023a), suggesting copious pristine H i gas inflows onto galaxies at this epoch before chemical enrichment from star formation can catch up to equilibrium (e.g., Torrey et al., 2018). The exact redshift for this transition is still debated (Curti et al., 2023b; Nakajima et al., 2023).
Constraining the inflow of neutral, pristine H i gas is key to understanding the assembly of the primordial matter and the build-up of stars and metals in the first galaxies. However, the circumstantial evidence from the FMR for excessive H i gas accretion needs to be corroborated by measuring directly the neutral, atomic hydrogen (H i) – the key missing element in the baryonic matter budget of high-redshift galaxies. Due to the weakness of the hyperfine H i 21 cm line, H i has historically been measured at redshifts through Lyman- (Ly) absorption in bright background sources such as quasars (Wolfe et al., 1986, 2005; Prochaska & Wolfe, 2009; Noterdaeme et al., 2012; Péroux & Howk, 2020) or gamma-ray bursts (Jakobsson et al., 2006; Prochaska et al., 2007; Fynbo et al., 2009; Tanvir et al., 2019; Heintz et al., 2023b). However, these only probe H i in narrow, pencil-beam sightlines and often at large impact parameters. Ly absorption has also been observed directly in integrated spectra of Lyman-break galaxies at (Pettini et al., 2000; Shapley et al., 2003; Steidel et al., 2010; Cooke & O’Meara, 2015; Lin et al., 2023), but only in rare, extreme cases.
Recently, Heintz et al. (2023d) reported the discovery of extremely strong ( cm-2) damped Ly absorption systems (DLAs) in star-forming galaxies at (see also Umeda et al., 2023; D’Eugenio et al., 2023; Chen et al., 2023). These seem to represent a galaxy population showing more prominent and substantially more prevalent DLAs than observed in Lyman-break galaxies at (e.g. Shapley et al., 2003). The shape of the Ly damping wings of galaxies had until these discoveries mainly been thought to trace the bulk H i of the intergalactic medium (IGM) at (Miralda-Escudé, 1998; McQuinn et al., 2008), and had been used to infer the reionization history of the Universe with JWST spectroscopy (Chen, 2023; Keating et al., 2023a, b; Umeda et al., 2023). The discovery of these massive interstellar or circumgalactic H i gas reservoirs carry important implications for our study of the early Universe: they constitute the first direct evidence for the mass accretion of primordial intergalactic gas into ordinary galaxies, i.e. these DLAs represent the creation of the baryonic component of galaxies. They also strongly affect observations of early galaxies because they change the shape of the Ly damping wings, altering and complicating inferences on the state of the bulk IGM. Further, they hinder the escape of ionizing photons (Laursen et al., 2011; Verhamme et al., 2015; Steidel et al., 2018; Hu et al., 2023; Hayes & Scarlata, 2023), and may also be responsible for the bias in photometric redshifts compared to spectroscopic measurements (Heintz et al., 2023d; Finkelstein et al., 2023b; Fujimoto et al., 2023a; Hainline et al., 2023) of galaxies during the reionization era. Their discovery, however, also presents a fortuitous new avenue to directly probe the build-up of pristine H i in galaxies at this critical epoch.
In this work, we systematically characterize these effects by measuring the prevalence and prominence of strong Ly emission and damped Ly absorption in a large spectroscopic sample of star-forming galaxies at – 13.4 observed with JWST/NIRSpec. The observational data for these galaxies are all presented here as part of the DAWN JWST Archive (DJA), which includes fully reduced and post-processed 1D and 2D spectra of all public JWST/NIRSpec data. The large JWST/NIRSpec archive forms the basis of the JWST/NIRSpec PRImordial gas Mass AssembLy (PRIMAL) legacy survey presented here. In this first part of the survey we aim to disentangle the Ly damping wings produced from H i gas in the immediate surroundings of galaxies (ISM or CGM) from the effects of an increasingly neutral IGM, and establish statistical correlations for these features and the presence of Ly emission with the physical properties of the galaxies. We further make the reduced spectra, source catalogs, spectroscopic redshifts, line identifications and fluxes, and the physical properties of each of the galaxies in the JWST-PRIMAL sample publicly available on dedicated webpages.
R.A. (deg) | Decl. (deg) | Prog. ID | Src. ID | S/NUV | Ref. | |
(1) | (2) | (3) | (4) | (5) | (6) | (7) |
1210 | 16374 | 60.53 | 5.5050 | (B23) | ||
1208 | 6521 | 13.39 | 5.5050 | – | ||
1345 | 381 | 12.02 | 5.5161 | – | ||
3215 | 201906 | 12.24 | 5.5194 | – | ||
2198 | 10840 | 30.90 | 5.5207 | – | ||
2750 | 10 | 25.564 | 11.49062 | (AH23) | ||
3215 | 20130158 | 61.61 | 11.7060 | (B23, ECL23) | ||
3215 | 20096216 | 42.17 | 12.5119 | (ECL23) | ||
2561 | 38766 | 24.95 | 12.7822 | (W23) | ||
3215 | 20128771 | 46.07 | 13.3605 | (ECL23) |
Notes. Column (1): Right ascension in degrees (J 2000). Column (2): Declination in degrees (J 2000). Column (3): Program ID under which the object was observed. Column (4): Designated MSA source ID. Column (5): Rest-frame UV signal-to-noise ratio. Column (6): Spectroscopic redshift. Column (7): The original survey paper references. A full version of this table can be found online.
References. (B23) Bunker et al. (2023a); (AH23) Arrabal Haro et al. (2023); (ECL23) Curtis-Lake et al. (2023); (W23) Wang et al. (2023a).
We have structured the paper as follows. In Sect. 2 we present the observations, describe the spectroscopic reduction and post-processing, and outline the JWST-PRIMAL sample compilation. In Sect. 3, we detail the spectroscopic analysis of Ly, the nebular emission line fluxes and equivalent widths, and the inferred physical properties of the sample galaxies. In Sect. 4, we consider the cosmic evolution of galaxy DLAs and Ly emitters and quantify the physical correlations driving or hampering the observed Ly emission and excess damped Ly absorption. In Sect. 5, we discuss and conclude our work and provide a future outlook. All the spectroscopic data products and the analysis and results presented in this work are made publicly available on dedicated webpages. Throughout the paper, we assume concordance flat CDM cosmology, with km s-1 Mpc-1, , and (Planck Collaboration et al., 2020).
2 Observations, data processing, and sample selection
2.1 The DAWN JWST Archive (DJA) – Spectroscopy
The data considered in this work are processed as part of the DAWN JWST Archive (DJA), an online repository containing reduced images, photometric catalogs, and spectroscopic data for public JWST data products111https://dawn-cph.github.io/dja. Here, we detail the spectroscopic data reduction and post-processing only, but see Valentino et al. (2023) for details on the imaging and photometric data. The DJA spectroscopic archive (DJA-Spec) is, at the time of writing, comprised of observations taken from some of the large Early Release Science (ERS), General Observer (GO) and Guaranteed Time (GTO) Cycle 1 & 2 programs, such as CEERS (ERS-1345; Finkelstein et al., 2023a), GLASS-DDT (DDT-2756 Treu et al., 2022; Roberts-Borsani et al., 2023) (ERS-1324; Treu et al., 2022), JADES (GTO-1180, 1210, GO-3215; Bunker et al., 2023a; Eisenstein et al., 2023), and UNCOVER (GO-2561; Bezanson et al., 2022). Many of the spectra analyzed in this work have been presented in dedicated survey papers cited above and in others (e.g., Arrabal Haro et al., 2023; Atek et al., 2023b; Cameron et al., 2023b), which also describe the sample selection and observational design of the different programs. With DJA-Spec we reduce and extract all of the spectra from this diverse array of observational programs with a single standardized pipeline.
2.2 DJA spectroscopic reduction and data processing
All data processing is performed with the publicly available Grizli (Brammer, 2023a)222https://github.com/gbrammer/grizli and MSAExp (Brammer, 2023b)333https://github.com/gbrammer/msaexp software modules, with the procedure as follows. The spectroscopic analysis begins with reprocessing the individual uncalibrated (uncal) exposures retrieved from the Mikulski Archive for Space Telescopes (MAST) with the JWST Detector1Pipeline.444jwst pipeline version 1.12.5; CRDS_CONTEXT = jwst_1180.pmap. The pipeline is executed with the default parameters, but with snowblind555https://github.com/mpi-astronomy/snowblind run between the jump and ramp_fit steps for improved masking of the bright cosmic ray “snowballs”(Rigby et al., 2023). We remove a column average from the count-rate (rate) files produced by the first step to remove the noise. We also find that scaling the read noise extension of the calibrated count-rate images by a factor of order derived separately for each exposure is necessary to explain the variance of pixels in un-illuminated portions of the NIRSpec detectors. The corrected rate files are then run through the calwebb_spec2 pipeline with its default parameters up through the photometric calibration (photom) step. The products at this stage are flux- and wavelength-calibrated 2D spectra for every source indicated with an open shutter in the attached MSA metadata, in a frame cut out from the original detector pixel grid with curved spectral traces. These 2D SlitModel cutouts are saved to individual multi-extension FITS files, and it is here that the subsequent processing with MSAExp departs significantly from the standard jwst pipeline.
The sky background of each source spectrum is removed by taking differences of the 2D cutout spectra taken at different “nod” positions of the telescope. Typically there are three nod offsets that shift the source by within the MSA slitlet (the so-called 3-Shutter-Slitlet pattern), though the pipeline calculates the groupings automatically to handle other configurations of the heterogeneous archival observations. Average 2D source and sky spectra (and associated variance images) are calculated from exposures at a given nod position and from those at all other offset positions, respectively. We fit a 2D Gaussian profile to the spectra where the width is the quadrature sum of a free parameter and the wavelength-dependent width of the point spread function and where a centroid shift is fit starting from an initial guess provided by the exposure MSA metadata and the catalog used to design the MSA mask. We note that the latter centroid shift is usually of order 0.1 pix (10 mas), which, while consistent with the expected precision of the input catalog astrometry and telescope pointing, is non-trivial and measurable and required for robust extraction of the 1D spectra. The sky-subtracted 2D spectra cutout spectra at different offset positions are combined in a single rectified frame with orthogonal wavelength and cross-dispersion pixel axes by binning with a 2D histogram algorithm. Finally, the one-dimensional spectra and uncertainties are extracted with “optimal” weighting (Horne, 1986) using the fitted 2D Gaussian profile. We note that Roberts-Borsani et al. (2024) have recently used the same pipeline presented here and developed for this work, presenting a similar effort to spectroscopically characterize galaxies at observed with NIRSpec Prism.
Currently (Apr 2, 2024) DJA-Spec includes 7,319 individual combined spectra taken with the NIRSpec Prism/CLEAR and 1,665 combined spectra with NIRSpec medium and high-resolution gratings. A full overview can be found on the dedicated webpage666https://s3.amazonaws.com/msaexp-nirspec/extractions/nirspec_graded_v2.html, where the following data are available: best-fit and inspected redshifts, the median signal-to-noise ratios (S/N) of the spectra, grades from visual inspection on the quality of the spectra, Hubble Space Telescope (HST) and JWST/NIRCam imaging thumbnails, and optimally reduced and extracted 2D and 1D spectra.
In this work, we focus on the JWST/NIRSpec Prism spectroscopic data from DJA due to the higher continuum sensitivity which is essential for quantifying the evidence for Ly emission or absorption damping wings. JWST’s NIRSpec Prism configuration covers the entire near-infrared passband from 0.6 m to 5.3 m, with a varying spectral resolution from in the blue end to in the red end (Jakobsen et al., 2022). We have carefully reduced and processed the spectra in DJA in a homogeneous way, limiting any biases between reduction and calibration codes across surveys. Further, this optimizes the consistency and uniformity of the JWST-PRIMAL sample measurements and the DJA-Spec sample as a whole.
2.3 JWST-PRIMAL sample selection
The main objective of this paper is to accurately model the Ly damping wings and/or measure strong Ly emission, and spectroscopically derive the physical properties of galaxies during the reionization epoch. For this, we require the following criteria for the sources to enter our archival JWST-PRIMAL sample:
-
1.
A robust, spectroscopic redshift measurement at from a minimum of one emission line detected at (in addition to the Ly break).
-
2.
An integrated signal-to-noise over the redshifted Ly and rest-frame UV () regions of .
-
3.
Flux density measurements over the entire NIRSpec Prism wavelength coverage from m.
These criteria ensure that the redshift of the foreground H i gas absorber, from interstellar to intergalactic scales, and the shape of the Ly line profile can be accurately measured. Further, the requirement of the full spectral coverage enables a full characterization of each source in the JWST-PRIMAL sample, including robust spectroscopic modelling of the stellar continuum and rest-frame UV spectral slope, emission line fluxes and equivalent widths EWs, and direct constraints on the star-formation rates (SFRs), metallicities, and ionization parameters for all sources.
In total, 494 sources from DJA-Spec meet these criteria. Fig. 1 show the absolute UV magnitude, , as a function of redshift for the JWST-PRIMAL sources and compared to the underlying DJA-Spec sample (at ). Table 1 summarizes the target list, including their coordinates, original program and source ID, and emission-line redshifts.
3 Analysis and results
Here we detail the spectroscopic measurements derived for the full JWST-PRIMAL sample. We focus on the Ly damping wings, but also detail the measurements and basic physical properties of the sample galaxies. All the spectroscopically-derived quantities are made available on a dedicated webpage777https://github.com/keheintz/jwst-primal.
3.1 The Lyman- damping parameter
During the reionization epoch at , several factors add to the line shape of Ly on integrated galaxy spectra: They can show strong Ly emission (LAE; Matthee et al., 2018; Mason et al., 2018; Witstok et al., 2023), damping wings from an increasingly neutral IGM (Miralda-Escudé, 1998; Keating et al., 2023a; Chen, 2023), excess continuum flux in the wings related to the size of the immediate ionized bubbles (McQuinn et al., 2008; Castellano et al., 2016, 2018; Fujimoto et al., 2023b; Umeda et al., 2023; Hayes & Scarlata, 2023), strong damped Ly absorption (DLA) from local H i gas (Heintz et al., 2023d; D’Eugenio et al., 2023), intrinsic variations due to a changing velocity and density distribution of gas and dust in the ISM and CGM (Dayal et al., 2011; Verhamme et al., 2015, 2017; Gronke et al., 2017; Hutter et al., 2023) or even possibly be dominated by two-photon emission processes (Steidel et al., 2016; Chisholm et al., 2019; Cameron et al., 2023a). Many of these effects are degenerate, so we have to disentangle them statistically. For this, we define a new simple diagnostic, which we denote the Ly damping parameter:
(1) |
which is the equivalent width (EW) of the transmitted flux density, over the wavelength region covered by the instrumentally-broadened Ly transition. The continuum flux, , over the same region is estimated by extrapolating the rest-frame UV slope, , which is derived directly in each spectra from rest-frame (see Sect. 3.4 for further details). This approximation of the stellar continuum is generally consistent with models of the continuum flux using galaxy templates (Heintz et al., 2023d). We define the integration limits from to (rest-frame) to capture most of the damping feature for the most extreme cases and to simultaneously limit the contamination from absorption or emission lines at longer wavelengths. We tested various lower limits and found this to be the bluest wavelength at which information is not lost at the resolution of NIRSpec Prism.
A set of physical models, shown as imprints on the Ly transmission curves, are visualized in Fig. 2 including: varying ionized bubble sizes from cMpc, neutral hydrogen fractions of the IGM from , and DLA H i column densities cm-2, all affecting the shape of the Ly damping wings. The IGM models are based on the prescription by Miralda-Escudé (1998), using the approximation described in Totani et al. (2006), and the UV bubble sizes are included by varying the upper bound of the IGM integration region, such that (e.g., McQuinn et al., 2008). The DLAs with varying H i column densities are derived from the Voigt-profile approximation by Tepper-García (2006). The corresponding Ly damping parameter is listed for each model. The UV bubble and DLA models all assume a benchmark IGM fraction of , and may therefore vary with the same scatter in as observed for the IGM-only models.
We derive damping parameters for the full sample in the range to . With this, we define four parameter spaces probing different physical regimes: (i) cases with excess continuum flux suggesting contributions from ionized bubbles or weak LAEs have –35 ; (ii) the effect on the Ly damping wings of an increasing neutral IGM results in –50 ; (iii) excess damped Ly absorption from local H i (with cm-2) have ; and (iv) strong LAEs are typically observed with . Fig. 3 shows three example spectra from the JWST-PRIMAL sample. One shows an example of a broad Ly damping wing consistent with local H i absorption with , consistent with cm-2. Another shows strong Ly emission with a measured .
This simple diagnostic now allows us to quantify the number of sources in each distinct population. We caution, however, that the low spectral resolution of the JWST/NIRSpec Prism configuration may conceal weak Ly emission in the continuum, which requires higher resolution spectroscopy to resolve (see e.g. the case for GN-z11; Bunker et al., 2023b). This will likely cause a slight decrease in . Further, as demonstrated in the bottom panel of Fig. 3, rare cases ( of the sample) that show both a substantial damped Ly broadening and simultaneously a strong Ly emission line will not be accurately identified, here with , misleadingly suggesting mild excess Ly continuum flux (the continuum flux of this particular source might be dominated by two-photon nebular emission, e.g., Cameron et al., 2023a). The JWST/NIRSpec shutter may also not cover the regions of the strongest Ly emission, thereby underestimating the total Ly photon output of the target source (Jiang et al., 2023). The parameter is thus most powerful in identifying the most extreme cases, such as strong Ly emitters or galaxy DLAs with cm-2, which uniquely have and , respectively.
[O ii] | [Ne iii] + | H | H + | H | [O iii] | H + | [S ii] 6725 | |
He i | [O iii] | [N ii] | ||||||
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) |
5.5050 | ||||||||
5.5050 | ||||||||
5.5161 | ||||||||
5.5194 | ||||||||
5.5207 | ||||||||
11.4906 | ||||||||
11.7060 | ||||||||
12.5119 | ||||||||
12.7822 | ||||||||
13.3606 |
Notes. All measurements are reported in units of erg s-1 cm-2. Uncertainties on measurements are stated at standard deviations and upper limits at confidence limits. Column (1): Spectroscopic redshift. Column (2): Unresolved [O ii] doublet line flux. Column (3): Unresolved [Ne iii] + He i line flux measurements, except for the highest redshift sources (). Column (4): Balmer H. Column (5): Unresolved Balmer H + auroral [O iii] line flux measurements, except for the highest redshift sources (). Column (6): Balmer H. Column (7): [O iii] line flux, where [O iii] is assumed to be of this. Column (8): Unresolved Balmer H+[N ii] . Column (9): Unresolved [S ii] doublet line flux. A full version of this table, including an expanded set of emission line identifications and measurements, can be found online.
3.2 Line flux measurements
For each source in the JWST-PRIMAL sample, we measure the spectroscopic redshifts from the most prominent nebular emission lines, when detected, including: the [O ii] and [O iii] doublets, [Ne iii] , He i , the [S ii] doublet, and the Balmer lines H, H, H, and . The [O ii] and [S ii] doublets are unresolved at all wavelengths in the NIRSpec Prism spectra, whereas the [O iii] doublet is marginally resolved at and fully resolved at due to the increasing spectral resolution with wavelength in this particular configuration. H and [N ii] are unresolved at all redshifts. [Ne iii] and He i are also generally unresolved, except for the highest redshifts at . Due to the broad wavelength coverage (m), H and [S ii] can be detected up to , H and the [O iii] doublet up to , and Ly can be detected at all redshifts from . The auroral [O iii] line is also in the majority of cases blended with H but becomes resolved at , enabling robust measurements of this important transition (Heintz et al., 2023d; Hsiao et al., 2023). Unfortunately, the [O iii] doublet is redshifted out of the NIRSpec Prism coverage at approximately the same redshifts, hindering direct -based metallicity measurements with the JWST/NIRSpec Prism configuration, except for a narrow redshift range ().
To determine the physical properties of the star-forming regions in each of the target galaxies, we measure or derive upper bounds on the line fluxes for each of these nebular emission lines when covered by the spectra. We fit the continuum around the lines with a simple polynomial and superimpose a set of redshifted Gaussian line profiles at the rest-frame wavelength of each transition and fit for the redshift , line equivalent widths (EWs) and fluxes. The line widths (full-width-at-half-maximum; FWHM) are in most cases directly proportional to the line-spread function of the NIRSpec Prism configuration (Jakobsen et al., 2022) and is modelled as such. For each case, we tie and the FWHM corrected by the spectral resolution, effectively assuming that the nebular emission lines all originate from the same H ii regions. An example of the line model fits is shown in Fig. 4. All the main identified lines, their derived line fluxes, or upper limits, are summarized in Table 2 for the JWST-PRIMAL sample and provided in full online888https://github.com/keheintz/jwst-primal, which also includes an expanded line list.
3.3 [O iii]+H equivalent widths
A simple diagnostic of the specific star-formation rate (sSFR) of galaxies is the emission line [O iii]+H EWs. These have been found to increase with increasing redshifts (Khostovan et al., 2016; Endsley et al., 2021; Matthee et al., 2023), and toward lower stellar masses and metallicities (Malkan et al., 2017; Reddy et al., 2018b), due to the more intense ionization fields and active star formation in early galaxies. We derive the [O iii]+H EW for each of the sources in our sample that have these lines covered by the spectra (311 sources at ). As demonstrated in Fig. 5, we integrate the flux density over the spectral region covering the H and [O iii] lines, normalized by the continuum flux, and correct the redshifted EWs to the rest frame. For each source, we extrapolate the continuum flux at the position of the lines by fitting a linear polynomial to the spectral regions on each side of the line complex.
In Fig. 6, we show the [O iii]+H EW distribution of the full JWST-PRIMAL sample and divided into redshift bins: –6.0, –7.0, –8.0, and –9.5. The median of the full distribution is , with 16th and 84th distribution percentiles from to . We observe no apparent evolution with redshift, but note that the number statistics is largely dominated by the lowest redshift bins. We also compare the full distribution to the much larger, albeit photometrically derived distribution from the CEERS survey of galaxies at –8.0 (Endsley et al., 2023). The two distributions are in good agreement, though the peak of the spectroscopic distribution is slightly shifted towards higher EWs. We note that the spectroscopically derived [O iii]+H EWs are much more accurate than photometric estimates, however, the latter potentially being underestimated by up to (Duan et al., 2023). The fact that EW distribution of the JWST-PRIMAL sample is only slightly shifted, but otherwise follow the same distribution, indicates that there are no strong biases in our spectroscopic sample toward strong emission-line sources from the underlying galaxy population.
3.4 Ultraviolet spectral slopes, magnitudes and luminosities
The UV spectral slope, , where , of the stellar continuum at wavelengths –2600 encodes key information regarding dust attenuation (Calzetti et al., 1994; Meurer et al., 1999), stellar metallicity (Cullen et al., 2021), the average stellar population age (Zackrisson et al., 2011), and the escape fraction of ionizing photons (Chisholm et al., 2022) from galaxies. The minimum value expected for standard stellar populations with a Galactic IMF and with negligible dust attenuation is (Cullen et al., 2017; Reddy et al., 2018a). Steeper spectral power-law indices of would require a very recent onset of star formation (Myr) or near-pristine, neutral gas and negligible nebular processing. JWST has enabled measurements of of galaxies at , hinting at even steeper spectral slopes, particularly for the faintest, less massive systems (Topping et al., 2023; Cullen et al., 2023b). These results have, however, been inferred from photometric data alone, which may be contaminated by varying Ly or metal emission (or absorption) line strengths or even low-luminosity AGN.
Here, we measure directly from the spectra for each of the JWST-PRIMAL sources at –13.4, carefully masking any identified emission lines or broad Ly absorption troughs. We use the emcee minimizer within the LMFit framework and recover the parameter covariances to estimate the median and 16th to 84th percentiles of the posterior distributions on and the normalization constant derived at UV rest-frame, of the best-fit model, . An example of one of the model fits is shown in Fig. 7. For our full sample, we recover UV continuum slopes ranging from down to a few cases with . This suggests that this set of rare galaxies are particularly dust- and metal-poor (), have mean stellar populations with ages Myr (e.g., Cullen et al., 2023a), and likely have ionizing photon escape fractions of (Chisholm et al., 2022). In Fig. 8, we show the derived spectral slopes as a function of redshift. We observe a marginal tendency for galaxies to have steeper UV slopes at higher redshifts, though still consistent with a non-evolution within the scatter in the sample distributions. We further note that distribution of dust and H iiregions may account for the large scatter in the population around the expected slope for a standard stellar population with negligible dust and maximum nebular continuum contribution (e.g., Calzetti et al., 1994; Vijayan et al., 2024).
To determine the UV luminosity, , and absolute magnitude, , for each source we derive the flux density and corresponding UV magnitude, at rest-frame from the spectral power-law model. We convert the apparent to the intrinsic, absolute brightness, , via
(2) |
where is the distance modulus at a given redshift derived from the astropy.cosmology module. The derived for each source are summarized in Table 3, not corrected for potential magnification factors for galaxies drawn from lensing cluster surveys. Our full sample spans absolute UV magnitudes of to mag, respectively (see Fig. 1).
[O iii]+H EW | SFRHβ,[OII] | ||||||
---|---|---|---|---|---|---|---|
(mag) | () | () | ( yr-1) | ||||
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) |
5.5050 | |||||||
5.5050 | |||||||
5.5161 | |||||||
5.5193 | |||||||
5.5207 | |||||||
11.4906 | |||||||
11.7060 | |||||||
12.5119 | |||||||
12.7822 | |||||||
13.3605 |
Notes. Column (1): Spectroscopic redshift. Column (2): Absolute UV (rest-frame ) magnitude. Column (3): Rest-frame [O iii]+H equivalent width. Column (4): Rest-frame UV () spectral slope. Column (5): The Ly damping parameter. Column (6): The star-formation rate derived from H () or [O ii] (). Column (7): The [O iii] [O ii] line flux ratio. Column (8): Gas-phase metallicity inferred from the joint PDF of the different strong-line calibrations from Sanders et al. (2024). A full version of this table can be found online.
3.5 Star-formation rate
The Balmer recombination lines originating from star-forming H ii regions, H in particular, are robust tracers of star formation on short ( Myr) timescales. The star-formation rate (SFR) has been found to increase at a given stellar mass (Whitaker et al., 2012; Speagle et al., 2014; Salmon et al., 2015; Thorne et al., 2021; Sandles et al., 2022) and rest-frame UV size (van der Wel et al., 2014; Ward et al., 2023; Langeroodi & Hjorth, 2023) for galaxies at increasing redshifts, following the evolution in the peak of the stellar mass halo mass relation with redshift (e.g., Behroozi et al., 2019). The SFR can now be measured directly from the Balmer recombination lines for galaxies at with JWST, revealing intense star-forming galaxies during the reionization epoch (e.g., Heintz et al., 2023a; Shapley et al., 2023; Fujimoto et al., 2023a; Sanders et al., 2024; Nakajima et al., 2023; Curti et al., 2023b). Given that H is only detected for the subset of the sample at , and may be contaminated by [N ii], we instead derive the SFR based on H for the majority of sources at as:
(3) |
following Kennicutt (1998), but here assuming the more top-heavy initial mass function (IMF) from Kroupa (2001), which is likely more representative of the high-redshift galaxy population (e.g., Steinhardt et al., 2023). We derive the H line luminosity via with the luminosity distance at the given redshift and assume from the Case B recombination scenario at K (Osterbrock & Ferland, 2006). We propagate the uncertainty of –30% from the exact choice of the IMF in the relation, in addition to the dependence on the electron density and temperature of the H ii region which is only of the order .
For the sources at , where H is redshifted out of the NIRSpec Prism spectral coverage, we instead rely on the [O ii] line strength
(4) |
again following Kennicutt (1998), and assuming the Kroupa (2001) IMF. For the full sample, we derive SFRs in the range 0.1–, not corrected for the relevant magnification factors or dust extinction, as summarized in Table 3. The SFRs can also be derived from overall UV luminosity of the sources following Madau & Dickinson (2014), where
(5) |
assuming solar metallicity which traces star formation on a longer timescale ( Myr) compared to the Balmer line ( Myr, see Calzetti, 2013). We find that the estimate is consistent with that derived from the Balmer recombination lines or [O ii] for the full sample, when allowing for up to mag of dust attenuation.
3.6 Gas-phase metallicity
One of the most fundamental characteristics of galaxies is their metallicity. This is notoriously difficult to measure directly, in particular for the highest redshift galaxies, motivating the use of strong-line diagnostics of the most prominent nebular emission lines (see e.g. Kewley et al., 2019; Maiolino & Mannucci, 2019, for recent reviews). In rare cases, the auroral [O iii] emission line can be detected directly, which enables sensitive measurements of the gas-phase metallicity (known as the direct -method), typically quantified as the oxygen abundance, (e.g., Izotov et al., 2006). This particular feature has now been detected in star-forming galaxies at with JWST (e.g., Schaerer et al., 2022; Curti et al., 2023a; Heintz et al., 2023c; Nakajima et al., 2023; Sanders et al., 2024), even out to (Hsiao et al., 2023; Heintz et al., 2023a, d). The auroral [O iii] line feature, however, will only be resolved from H in the minority of sources due to the increasing spectral resolution towards the red end of the spectra.
We therefore base our main metallicity measurements on strong-line calibrations. To preserve the homogeneity of the measurements, we use a single set of line diagnostics from Sanders et al. (2024) for the different line ratios, derived for galaxies at –9 through the direct -method. Specifically, we consider the strong-line ratios: [O iii] , [O ii] , ([O iii] + [O ii] ) / H, [Ne iii] / [O ii] , and [O iii] [O ii] , which are available at most of the targeted redshifts. Based on the line-flux measurements and errors (representing the likelihood distributions of fluxes), we construct the probability density function (PDF) for each diagnostic ratio where the width denotes the scatter in the relation. We then combine the PDFs from the available calibrations, inversely weighted by their scatter, for each individual source to determine the gas-phase metallicity through the median and 16th to 84th percentiles of the joint PDF, as illustrated in Fig. 9. This more conservative approach takes into account the added uncertainty from sources with potential high ionization parameters, typically quantified via , thus minimizing the uncertainty from assuming a single line-ratio.
4 Lyman- emission and absorption in galaxies during the reionization epoch
With the spectroscopically-derived physical properties for the full set of JWST-PRIMAL sample sources we can now investigate and chart the redshift evolution and the physical driver of strong galaxy DLAs produced by massive H i gas reservoirs to prominent ionized bubbles and Ly emission during the reionization epoch.
4.1 Charting as a function of redshift
To obtain a first overview of the prevalence and prominence of DLAs and LAEs across time in the reionization epoch, we show as a function of redshift for the full sample in Fig. 10. This represents the evolution from to , corresponding to approx. 300 Myr to 1 Gyr after the Big Bang. The parameter space has been color-coded to highlight the regimes driven by the physical models for the expected dominant contribution to the observed Ly transmission curve. The number statistics of sources located in each region are summarized in Table 4, divided into redshift bins of –6.0 (159 sources), –7.0 (169), –8.0 (103), –10.0 (46), and –13.4 (17).
We find that there are 38 strong LAEs in the full JWST-PRIMAL sample (i.e. ), based on . None of these are found at . At , we find that 38 out of 494 of sources show strong LAEs () in approximately equal fractions from –6.0, –7.0, and –8.0, albeit with the most prominent LAEs quantified by the lowest at . This fraction of LAEs could be as high as if we include all sources with excess Ly flux as potential weakly-emitting LAEs. These results are slightly lower than found previous literature studies, showing a typical fraction of LAEs of at –8 (e.g., Saxena et al., 2023b; Nakane et al., 2023; Simmonds et al., 2023; Witstok et al., 2023; Jung et al., 2023; Jones et al., 2023; Witten et al., 2024). However, they were also found to be much more sparsely populated beyond (though see Tang et al., 2023; Bunker et al., 2023b; Fujimoto et al., 2023b, Witstok et al. in prep.), consistent with our results, and as also recovered by simulations (Hutter et al., 2023). These findings would also naturally explain the null detections of LAEs at in large photometric surveys using narrow-band imaging such as the UltraVISTA (e.g., Laursen et al., 2019).
Redshift range | |||||
---|---|---|---|---|---|
LAEs | 18/159 () | 12/169 () | 7/103 () | 1/46 () | 0/17 () |
Ly exc. | 13/159 () | 15/169 () | 4/103 () | 2/46 () | 1/17 () |
IGM abs. | 34/159 () | 28/169 () | 23/103 () | 3/46 () | 5/17 () |
DLAs | 94/159 () | 114/169 () | 69/103 () | 40/46 () | 11/17 () |
Similarly, we consider the redshift evolution of strong DLAs, here classified as sources with , corresponding to H i column densities cm-2. We find that the fraction of galaxy DLAs only slightly increases from at up to –90% at . A similar conclusion is reached by Umeda et al. (2023), who find an ubiquitous presence of extremely strong DLAs (cm-2) at . This evolution is more clearly highlighted in Fig. 11 where we show the histograms of for the five different redshift ranges, demonstrating higher prevalence of galaxy-integrated DLAs at increasing redshifts. These results are in sharp contrast to previous studies of Lyman-break galaxies near the peak of cosmic star formation at –3 (Pettini et al., 2000; Steidel et al., 2010), where only a small subset of the population () show integrated galaxy line-of-sight H i column densities of cm-2 (Shapley et al., 2003). This corroborates our findings here of an increasing probability of having high H i column densities or larger fraction of DLAs at higher redshifts, down to –3. This may indicate that we are probing the epoch of first-galaxy assembly with excessive amounts of accreting or overdense pristine H i gas that are yet to be processed into stars or ionized by the first stellar sources. Indeed, the average gas mass density of galaxies have been inferred to increase by an order of magnitude through accretion from to (see Walter20; Heintz et al., 2022, and Sect. 5 for further discussion). (see Sect. 5 for further discussion). Fig. 11 also demonstrates an excess of galaxies with strong LAEs () at –6.0 from the full distribution, distinct for the epoch when the universe is fully ionized.
We note that the fraction of galaxy DLAs at may be a lower bound, given that our selection criteria on the S/NLyα will preferentially be biased towards strong Ly emission. On the other hand, fainter, less massive sources are likely to have stronger Ly EWs (Saxena et al., 2023a) and contribute to most of the ionizing photons (Lin et al., 2024). To investigate this, we compare to the spectroscopically-derived in Fig. 12. We confirm the trend of the intrinsically faintest sources showing the highest Ly EWs, with the strongest LAEs with showing luminosities predominantly at mag. This further suggests that the large H i gas overdensities are preferentially located close to bright, massive sources. In the following sections, we will further investigate the underlying physical properties driving the prevalence and prominence of strong Ly emission and absorption in these reionization era sources.
4.2 The role of stellar ionization fields on Ly emission and absorption
The intense UV radiation fields produced by younger, more massive stars are likely to play a key role in the ionization of H i and the production and escape of Ly photons. To test this and quantify it in terms of our observations, we compare to two common tracers of active star formation and high ionization: The [O iii]+H EW, and the ionization parameter = [O iii] / [O ii] . We find no apparent correlations between and any of the two quantities. For example, strong galaxy DLAs are observed in galaxies spanning the entire dynamical range of [O iii]+H EWs and ionization. Similarly, strong LAEs are observed in galaxies spanning [O iii]+H EWs of , and are not directly associated with sources with high or low ratios. This is at odds with studies of local “green pea” galaxies at , believed to be robust analogs to high- star-forming galaxies, that find that larger Ly emission EWs and escape fractions are commensurate with higher, more intense stellar ionization fields (Yang et al., 2017; Jaskot et al., 2019; Hayes et al., 2023). We note, however, that the galaxies at predominantly show high absorption and ionization and instead typically represent the large-scale environment or H i gas covering fraction at these redshifts. Consequently, our results disfavor a scenario where the H i abundance (or the lack thereof) is driven by the ionization output of the stellar sources in the galaxy.
4.3 Gas-phase metallicity and Ly absorption and emission
The gas-phase metallicity of galaxies is one of the key regulators of the gas cooling efficiency, the ionization state of the gas, and the escape fraction of Ly and LyC photons. Lower metallicities are often associated with higher ionization and temperatures, effectively reducing the observed line-of-sight H i column density. By contrast, high metallicity and likely thereby also high dust content, will enable less efficient ionization of H i and as a consequence more efficiently scatter the output ionizing photons, effectively suppressing the intrinsic and observed escape fraction of Ly photons (Dayal et al., 2011; Henry et al., 2015; Dijkstra et al., 2016; Verhamme et al., 2017; Yang et al., 2017; Izotov et al., 2018; Jaskot et al., 2019). We therefore expect a strong correlation between the gas-phase metallicity and for our high-redshift sample sources.
To test this scenario, we show as a function of gas-phase metallicity, (O/H), for each of the sources in the JWST-PRIMAL sample in Fig. 14. The strongest LAEs () are predominantly observed in galaxies with low metallicities, and ubiquitously found at , corresponding to solar metallicity. These results are consistent with the expectations that prominent Ly emission can only escape from low-metallicity systems. By contrast, we find that galaxies with strong DLAs () span the entire dynamical range of inferred metallicities, – (i.e. 0.6–30% solar). The most prominent DLA systems are found at lower metallicities at any redshifts , supporting the scenario where these mostly trace young, less chemically enriched galaxies.
4.4 The impact of galaxy DLAs on the reionization history
One of the major implications of high H i-cloumn density DLAs observed in galaxy-integrated spectra at is their effect on the escape fraction of Ly and LyC photons, key to constrain the reionization history and topology. Numerical modelling and reionization simulations have previously only taken IGM absorption into account when predicting the radiative transfer and transmission of Ly (e.g., Laursen et al., 2019; Hutter et al., 2023). These overabundant local H i gas reservoirs are likely to be the main component in the transfer of ionizing photons, at least in the early evolutionary stages of galaxies, so quantifying this additional local ISM or CGM absorption is vital to draw robust conclusions and physically interpret the output and impact of galaxy-wide DLAs in simulations.
To gauge this effect, we first compare to the spectral slopes of each galaxy in Fig. 15. has been found to be strongly correlated with the escape fraction of Lyman Continuum (LyC) photons, , in low-redshift galaxies sample (Chisholm et al., 2022; Flury et al., 2022). We color-grade the parameter space with the empirical relation presented in their paper Chisholm et al. (2022, their Eq. 11), highlighting the inferred escape fractions from to . We find no apparent correlations between or the inferred with the strength of the Ly emission or DLA absorption. This indicates that the hardness of the UV spectral shape is not the main driver of the ionization of H i, and that larger H i abundances do not translate directly into lower escape fractions of ionizing photons. This is somewhat at odds with theoretical predictions and local observations, since both Ly and LyC photons are affected by resonant scattered from abundant interstellar or circumgalactic H i gas columns (e.g., Flury et al., 2022; Xu et al., 2022; Begley et al., 2024). We caution, however, that the relations for derived from the local galaxy sample at studied by Chisholm et al. (2022) may not be universally representative for galaxies at (Choustikov et al., 2024; Pahl et al., 2024). Since the high- galaxy population show even more abundant H i gas reservoirs than galaxies at equivalent metallicities or ionization parameters (Heintz et al., 2023d), the escape fraction of Ly and LyC photons is likely correspondingly lower. Indeed, comparing the sources from Chisholm et al. (2022) that are mostly absorbed by the neutral H i gas, the inferred LyC escape fraction could be a lower.
To more directly relate the amount of photons capable of ionizing H i to the presence and strength of Ly emission or DLAs, we also compare the ionizing photon production efficiency, , to in Fig. 15. can be derived through the Balmer recombination line strengths, such as H, since these accurately trace the ionizing photon production rate from massive stars on short ( Myr) timescales (Bouwens et al., 2016). Following Matthee et al. (2023), we define
(6) |
where erg is the H line-emission coefficient, assuming a Case B recombination scenario with K (e.g., Schaerer, 2003). We find from Fig. 15 that galaxies with strong DLAs () occupy most of the probed parameter space, though with a median of , higher than the canonical values, (Robertson et al., 2013). The large scatter in at could indicate bursty star formation, particularly since traces massive O-stars, and may reflect typical younger stellar populations that are yet to ionize the neutral H i gas in their local surroundings. By contrast, galaxies that are strong LAEs tend to exhibit high ionizing photon production efficiencies, occupying the region with , higher than the canonical value. Since these are predominantly found at , we argue that it is mainly this particular subset of galaxies with high and late in the reionization epoch that provides the strongest contribution to this last phase transition.
5 Discussion
To place our results into context of the large-scale efforts to quantify the galaxy mass assembly and reionization timeline, we describe in this section a physical scenario connecting the observed prevalence and prominence of Ly emission and DLA absorption from abundant, pristine H i gas in galaxies through redshift and as a function of their intrinsic physical properties. We further quantify the effect on from the neutral hydrogen fraction based on more sophisticated IGM simulations and estimate the effect of strong galaxy DLAs on photometric redshifts measurements at , which appear to be systematically overestimated.
5.1 H i gas assembly and reionization across cosmic time
In this work, we quantified the strength of DLA absorption (Ly emission) via the positive (negative) strength of . We found that galaxies with high , indicating substantial H i gas column densities in the ISM or CGM of galaxies, are ubiquitously observed at and in galaxies spanning the entire dynamical range of physical properties probed here. The fact that there appears to be no correlation between the amount of local H i gas (or the lack thereof) and the ionizing output of the galaxies themselves may suggest that DLAs trace particularly young galaxies, at all redshifts, that are yet to ionize the surrounding H i gas or convert it into molecules and eventually stars. For galaxies with , equivalent to cm-2, the absorbing H i gas is completely self-shielding of ionizing photons, even with substantial clumping with density variations down supporting this scenario.
On the contrary, strong LAEs are only observed at . This indicates that only after , on average, is the surrounding large-scale IGM sufficiently ionized for a prominent fraction of Ly photons to escape. Indeed, the observed redshift evolution of is in good qualitative agreement with the evolution of the Ly escape fraction, as compiled by Saxena et al. (2023a); Hayes & Scarlata (2023). Consequently, this favors more late and rapid reionization models (e.g., Naidu et al., 2020), over a smooth, early transition (e.g., Finkelstein et al., 2019). A similar rapid evolution seems to be suggested by damped Ly wings in quasar spectra (Greig et al., 2017, 2019; Davies et al., 2018; Bañados et al., 2018; Yang et al., 2020; Wang et al., 2021; Ocvirk et al., 2021; Ďurovčíková et al., 2020, 2024) and Ly emission statistics of galaxies (Ouchi et al., 2010; Sobacchi & Mesinger, 2015; Mason et al., 2018, 2019), indicating a near-fully neutral IGM at .
The strong LAEs are predominantly associated with UV faint ( mag) galaxies with high ionizing photon production efficiencies, . This could either indicate that these particular systems were the most efficient at ionizing their immediate surroundings, or that viewing angles play a larger effect at low () redshifts compared to their higher redshift counterparts. In this scenario, galaxies at are deeply embedded into dense, neutral gas overdensities with large covering fraction, showing strong absorption independent of the viewing angle. At lower redshifts, low-density channels are carved by ionizing photons and we thus observe a broader span in at . This scenario would be commensurate with the non-dependence of the strength of LAEs with the ionization parameter, [O iii]+H EW, the spectral slope and the inferred escape fraction of LyC photons. In fact, strongly directionally-dependent escape of ionizing radiation has been proposed by simulations (Cen & Kimm, 2015; Rosdahl et al., 2022; Yeh et al., 2023). Since only of the galaxies at have strong LAEs, we argue that it is either this particular subset of faint, low mass galaxies that drive reionization or that ionizing photons from galaxies at this epoch will only be able to escape through the low-density sightlines.
5.2 Insights into the parameter from simulations
So far, the IGM regime of has been investigated assuming a homogeneously ionized IGM with an average value of the ionized fraction everywhere. It is therefore interesting to investigate how the parameter evolves in simulations of patchy reionization, to better understand the allowed range of and what some limitations may be.
Attenuation by the intergalactic medium is expected to play a significant role in regulating visibility of Ly emission until the Universe was fully reionized (McQuinn et al., 2007), and will impact the parameter to some extent over all of the regimes identified in Fig. 10. A full comparison of observed and simulated across the whole parameter space would therefore require detailed modelling of the redshift evolution of ionized bubbles, as well as the detailed shape and of the Ly emission profile (e.g., Hayes & Scarlata, 2023). We instead focus here on identifying the allowed range of for which attenuation by the IGM dominates the measurement. We measure from IGM transmission spectra normalised by the continuum, i.e., with no contribution from Ly emission.
First, we show how the Ly transmission curves change based on the residual fraction of H i inside an ionized bubble. As discussed in several works, even a small amount of neutral gas can completely saturate absorption of the Ly forest (Mesinger & Haiman, 2007; Mason et al., 2018; Keating et al., 2023a). The difference in accounting for this residual neutral gas is shown in Fig. 16, which shows IGM transmission curves calculated using the Miralda-Escudé (1998) analytic model. In the top left panel, the bubbles are assumed to be completely ionized, allowing for transmission on the blue side of Ly. When this neutral gas is accounted for, a sharp cutoff in IGM transmission blueward of the Ly emission is observed. The bottom panel shows the same transmission profiles, but now convolved with the NIRSpec PRISM instrumental profile. It is difficult to statistically disentangle the difference between the models with different bubble sizes once the residual H i is included.
Based on these results, we estimate an allowed range of IGM-dominated parameters, shown in Fig. 17. The lower value of is given by a step function, where there is no transmission blueward of Ly and complete transmission redward of Ly Ṫhis mimics the limit of large bubbles containing a small amount of residual H i. The upper value of is calculated from the Miralda-Escudé (1998) analytic model assuming the Universe is completely neutral and assuming a bubble size cMpc. We see that the maximum value of grows with redshift, due to the evolution of the mean density of the Universe (Keating et al., 2023b).
For comparison, we also measure IGM-dominated parameters directly from damping wings constructed from a simulation of inhomogeneous reionization. We use the spectra presented and described in Keating et al. (2023a), which were generated from sightlines through a simulation from the Sherwood-Relics simulation suite (Puchwein et al., 2023). We analyse results from five snapshots from this simulation at . The distributions of parameters we measure are shown as the violin plots in Fig. 17. We find that the range of allowed IGM parameters is well captured by the grey shaded region we defined. At , some sightlines have smaller parameters as the Ly forest starts to show transmission at this redshift. At , some sightlines have larger parameters due to effects not captured in the analytic model for IGM damping wings, such as dense clumps of gas close to the host halo or the effect of infalling material. In general, however, the range of parameters for which the reionization of the IGM plays a dominant role are in good agreement with the range of – as defined in Section 3.1. However, we note that the IGM will also play an important role in regulating the fraction of flux transmitted by galaxies showing evidence of Ly emission, and hence falling into the regime of smaller values of .
5.3 The effect of Ly damping wings on high- photometric redshift estimates
Some of the first efforts to photometrically identify high-redshift galaxy candidates seemed to generally overpredict their redshifts when compared to the subsequent spectroscopic redshift determination of the sources at based on strong nebular emission lines (e.g., Finkelstein et al., 2023b; Fujimoto et al., 2023a; Hainline et al., 2023; Steinhardt et al., 2023; Serjeant & Bakx, 2023). In Heintz et al. (2023d), we proposed that this may due to the strong damping wings observed in the large fraction fraction of galaxies at these redshifts, suppressing the derived flux density in the blue-most photometric bands and thereby mimicking a slightly larger Ly “break redshift” (e.g., Arrabal Haro et al., 2023).
To test this, we show the derived photometric redshifts for the subset of sources where these are available in the DJA (Valentino et al., 2023)999Available at https://dawn-cph.github.io/dja/imaging/v7/ relative to the spectroscopically derived redshifts in Fig. 18. We find that the photometric redshifts are systematically overestimated for galaxies at , with median . Considering only the galaxies with strong DLAs () the discrepancy become slightly larger, with median . To optimize future photometric redshift estimates, the effects of strong Ly damping wings needs to be incorporated into spectral energy distribution fitting algorithms like Eazy (Brammer et al., 2008).
6 Summary & Future Outlook
In this work, we have compiled and systematically characterized 494 sources observed with JWST/NIRSpec during the reionization epoch, at –13.4, introducing the JWST-PRIMAL legacy survey. The main goal was to statistically quantify the presence and redshift evolution of strong damped Ly absorption in galaxies, signifying abundant local H i gas reservoirs. To disentangle these systems from sources with damping parameters probing mainly a largely neutral IGM, or residing in extended ionized bubbles, we defined a new simple diagnostic, the Ly damping parameter, . This allowed us to statistically chart the H i gas assembly history, the onset and prevalence of strong Ly-emitting galaxies and to explore the main underlying physical properties driving the emission and escape of ionizing and Ly photons.
We found that the majority () of galaxies at are dominated by large H i gas reservoirs in their local surroundings (ISM or CGM) with galaxy-integrated H i column densities, cm-2. These strong galaxy DLAs appear to exist throughout the reionization epoch, even though the fraction was observed to decrease to by of the overall galaxy population. Similar H i column densities have been theoretically predicted from cosmological simulations (e.g., Finlator et al., 2018; D’Odorico et al., 2018; Pallottini et al., 2022). Further, we found that these strong DLA galaxies represent a large variety of intrinsic physical properties in terms of their UV luminosity, gas-phase metallicity (O/H), UV spectral slope , nebular [O iii]+H line EWs, ionization parameters, and ionizing photon production efficiency, . We thus surmised that strong galaxy DLAs are a sign of early galaxies in the process of formation, actively accreting and assembling large amount of primordial H i gas, that are yet to be ionized or processed into molecular gas and stars.
From the measured redshift evolution of we were able to quantify the onset and increasing strength of Ly, enabling us to decode the reionization history directly. We found that at –10 galaxy spectra mainly probe a substantially neutral CGM, whereas evidence for low Ly escape or significant ionized bubbles around these galaxies were only observed at later cosmic times (). The fraction of LAEs was estimated to be at . The strongest LAEs at have been found to reside in particular overdense or ionized regions (Saxena et al., 2023b; Witstok et al., 2023; Witten et al., 2024). From theoretical considerations and cosmological simulations (Pallottini et al., 2022), we expect most of the galaxies studied here to be formed within the first 200–300 Myr of cosmic time at –14 and rapidly build up their stellar mass and on-going star formation via substantial H i gas accretion and assembly. Over the next 100–300 Myr, down to ,
The emergence of these sources showing strong “roll-overs” near the Ly edge in the first era of JWST/NIRSpec Prism spectra have important practical implications for future studies and observations: First, they complicate statistical inferences of the reionization history by tracing the increasing neutral fraction of the IGM from galaxy damping wing measurements (Keating et al., 2023a). Here, we provided the first step forward to statistically model the galaxy DLA population which will eventually be possible to remove to uncover the underlying IGM damping effects. Second, the first efforts in photometrically identifying high-redshift galaxy candidates seemed to show a systematic positive offset in towards higher redshifts (Finkelstein et al., 2023b; Fujimoto et al., 2023a; Hainline et al., 2023; Steinhardt et al., 2023; Serjeant & Bakx, 2023). This may naturally be explained by the increasing fraction of galaxies with stronger Ly damping wings at higher redshifts. This we quantified here and estimated a mean of 0.2 dex overestimations of the photometric redshifts of galaxies at . This effect would need to be incorporated into photo- codes to optimize the accuracy on the priors for the earliest galaxy population.
The next natural avenues to explore the early galaxy DLAs in more depth are to examine the effect of clustering, their correlations with physical sizes and SFR surface densities, and the H i gas content; a potential cause of the observed offset in the otherwise fundamental-metallicity relation (FMR). Similar to how strong LAEs are typically found to reside in high-density groups of galaxies, galaxy DLAs may trace overdense regions via their early assembly of H i. Moreover, the subset of galaxies showing strong DLAs integrated over their relative small UV sizes indicate exceptional dense gas surface densities (Kennicutt & Evans, 2012). This, in combination with the inferred SFRs, may be key to understand the intense star formation and bright UV luminosities of early galaxies (e.g., Finkelstein et al., 2023a; Franco et al., 2023; Harikane et al., 2023; Adams et al., 2023; Casey et al., 2023; Bouwens et al., 2023; Chemerynska et al., 2023; Mason et al., 2023; McLeod et al., 2024). Finally, the discovered offset from the FMR towards lower metallicities of galaxies at –10 have been interpreted as evidence for excessive pristine H i gas inflow (Heintz et al., 2023a), which have later been confirmed with larger galaxy samples, though the exact redshift for this transition is still debated (e.g., Curti et al., 2023b; Nakajima et al., 2023). This scenario will be possible to test directly by correlating the H i column densities of galaxies to their offset from the FMR, moving the 3D plane of galaxy properties to 4 dimensions. Obtaining higher-resolution spectroscopic data of these sources will further help quantifying the H i covering fraction and disentangle the effects on the Ly transmission curve from weak Ly emission to extended ionized bubbles and expand the probed physical properties to include the electron densities and temperatures of the gas.
While previous efforts with ALMA to quantify the H i gas mass content using far-infrared tracers such as [C ii]m and [O i]m have been successfully applied to galaxies in the epoch of reionization (Heintz et al., 2021, 2022; Vizgan et al., 2022; Liang et al., 2024; Wilson et al., 2023), the results presented here provide the most direct and statistical census of the abundance and assembly of the primordial H i gas in early Universe. This naturally paves the way and provide an important benchmark for next-generation radio facilities, targeting H i of the first galaxies through the 21-cm transition, such as the square kilometre array (SKA; Dewdney et al., 2009).
Acknowledgements
We would like to thank Peter Jakobsen for his vision and heroic endeavour in optimally designing the JWST/NIRSpec instrument and some of its first on-sky observations and for enlightening discussions about the intricacies of the NIRSpec data. Further, we would like to thank John Chisholm for helpful clarifications and discussions related to the escape fraction of ionizing photons and Aayush Saxena for enlightening conversations on the escape and absorption of Lyman- photons.
This work has received funding from the Swiss State Secretariat for Education, Research and Innovation (SERI) under contract number MB22.00072. The Cosmic Dawn Center (DAWN) is funded by the Danish National Research Foundation under grant DNRF140. The data products presented herein were retrieved from the DAWN JWST Archive (DJA). DJA is an initiative of the Cosmic Dawn Center, which is funded by the Danish National Research Foundation under grant DNRF140. P.D. acknowledge support from the NWO grant 016.VIDI.189.162 (“ODIN”) and warmly thanks the European Commission’s and University of Groningen’s CO-FUND Rosalind Franklin program. Support from the ERC Advanced Grant INTERSTELLAR H2020/740120 is kindly acknowledged (A.F.). S.G. acknowledges financial support from the Villum Young Investigator grants 37440 and 13160 and the Cosmic Dawn Center. M.K. was supported by the ANID BASAL project FB210003. G.E.M. acknowledges financial support from the Villum Young Investigator grants 37440 and 13160 and the Cosmic Dawn Center. J.W. acknowledges support from the Science and Technology Facilities Council (STFC), by the ERC through Advanced Grant 695671 “QUENCH”, by the UKRI Frontier Research grant RISEandFALL. Support for this work was provided by NASA through the NASA Hubble Fellowship grant HST-HF2-51515.001-A awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Incorporated, under NASA contract NAS5-26555. F.C. acknowledges support from a UKRI Frontier Research Guarantee Grant (PI Cullen; grant reference EP/X021025/1). J.H.W. acknowledges support by NSF grant AST-2108020 and NASA grants 80NSSC20K0520 and 80NSSC21K1053. NRT acknowledges support through STFC consolidated grant ST/W000857/1. M.J.H. is supported by the Swedish Research Council, Vetenskapsrdet, and is fellow of the Knut & Alice Wallenberg foundation.
This 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 (MAST) 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.
Data availability statement
All the data presented in this work are made publicly available through DJA: https://dawn-cph.github.io/dja/, and are made easy to inspect and query in a dedicated webpage: https://s3.amazonaws.com/msaexp-nirspec/extractions/nirspec_graded_v2.html. The spectroscopic data have been processed using the custom-made reduction pipeline MsaExp, publicly available here: https://github.com/gbrammer/msaexp. All the data and catalogs related to JWST-PRIMAL are also available on a dedicated webpage: https://github.com/keheintz/jwst-primal.
References
- Adams et al. (2023) Adams, N. J., Conselice, C. J., Austin, D., et al. 2023, arXiv e-prints, arXiv:2304.13721
- Arellano-Córdova et al. (2022) Arellano-Córdova, K. Z., Berg, D. A., Chisholm, J., et al. 2022, ApJ, 940, L23
- Arrabal Haro et al. (2023) Arrabal Haro, P., Dickinson, M., Finkelstein, S. L., et al. 2023, Nature, 622, 707
- Asplund et al. (2009) Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481
- Atek et al. (2023a) Atek, H., Chemerynska, I., Wang, B., et al. 2023a, MNRAS, 524, 5486
- Atek et al. (2023b) Atek, H., Labbé, I., Furtak, L. J., et al. 2023b, arXiv e-prints, arXiv:2308.08540
- Austin et al. (2023) Austin, D., Adams, N., Conselice, C. J., et al. 2023, ApJ, 952, L7
- Bañados et al. (2018) Bañados, E., Venemans, B. P., Mazzucchelli, C., et al. 2018, Nature, 553, 473
- Begley et al. (2024) Begley, R., Cullen, F., McLure, R. J., et al. 2024, MNRAS, 527, 4040
- Behroozi et al. (2019) Behroozi, P., Wechsler, R. H., Hearin, A. P., & Conroy, C. 2019, MNRAS, 488, 3143
- Bezanson et al. (2022) Bezanson, R., Labbe, I., Whitaker, K. E., et al. 2022, arXiv e-prints, arXiv:2212.04026
- Blumenthal et al. (1984) Blumenthal, G. R., Faber, S. M., Primack, J. R., & Rees, M. J. 1984, Nature, 311, 517
- Bosman et al. (2022) Bosman, S. E. I., Davies, F. B., Becker, G. D., et al. 2022, MNRAS, 514, 55
- Bouwens et al. (2016) Bouwens, R. J., Smit, R., Labbé, I., et al. 2016, ApJ, 831, 176
- Bouwens et al. (2023) Bouwens, R. J., Stefanon, M., Brammer, G., et al. 2023, MNRAS, 523, 1036
- Brammer (2023a) Brammer, G. 2023a, grizli, Zenodo
- Brammer (2023b) Brammer, G. 2023b, msaexp: NIRSpec analyis tools, Zenodo
- Brammer et al. (2008) Brammer, G. B., van Dokkum, P. G., & Coppi, P. 2008, ApJ, 686, 1503
- Brinchmann (2023) Brinchmann, J. 2023, MNRAS, 525, 2087
- Bunker et al. (2023a) Bunker, A. J., Cameron, A. J., Curtis-Lake, E., et al. 2023a, arXiv e-prints, arXiv:2306.02467
- Bunker et al. (2023b) Bunker, A. J., Saxena, A., Cameron, A. J., et al. 2023b, A&A, 677, A88
- Calzetti (2013) Calzetti, D. 2013, in Secular Evolution of Galaxies, ed. J. Falcón-Barroso & J. H. Knapen, 419
- Calzetti et al. (1994) Calzetti, D., Kinney, A. L., & Storchi-Bergmann, T. 1994, ApJ, 429, 582
- Cameron (1957) Cameron, A. G. W. 1957, PASP, 69, 201
- Cameron et al. (2023a) Cameron, A. J., Katz, H., Witten, C., et al. 2023a, arXiv e-prints, arXiv:2311.02051
- Cameron et al. (2023b) Cameron, A. J., Saxena, A., Bunker, A. J., et al. 2023b, A&A, 677, A115
- Casey et al. (2023) Casey, C. M., Akins, H. B., Shuntov, M., et al. 2023, arXiv e-prints, arXiv:2308.10932
- Castellano et al. (2016) Castellano, M., Dayal, P., Pentericci, L., et al. 2016, ApJ, 818, L3
- Castellano et al. (2022) Castellano, M., Fontana, A., Treu, T., et al. 2022, ApJ, 938, L15
- Castellano et al. (2024) Castellano, M., Napolitano, L., Fontana, A., et al. 2024, arXiv e-prints, arXiv:2403.10238
- Castellano et al. (2018) Castellano, M., Pentericci, L., Vanzella, E., et al. 2018, ApJ, 863, L3
- Cen & Kimm (2015) Cen, R. & Kimm, T. 2015, ApJ, 801, L25
- Chemerynska et al. (2023) Chemerynska, I., Atek, H., Furtak, L. J., et al. 2023, arXiv e-prints, arXiv:2312.05030
- Chen (2023) Chen, H. 2023, arXiv e-prints, arXiv:2307.04797
- Chen et al. (2023) Chen, Z., Stark, D. P., Mason, C., et al. 2023, arXiv e-prints, arXiv:2311.13683
- Chisholm et al. (2019) Chisholm, J., Rigby, J. R., Bayliss, M., et al. 2019, ApJ, 882, 182
- Chisholm et al. (2022) Chisholm, J., Saldana-Lopez, A., Flury, S., et al. 2022, MNRAS, 517, 5104
- Choustikov et al. (2024) Choustikov, N., Katz, H., Saxena, A., et al. 2024, arXiv e-prints, arXiv:2401.09557
- Cooke & O’Meara (2015) Cooke, J. & O’Meara, J. M. 2015, ApJ, 812, L27
- Cullen et al. (2023a) Cullen, F., McLeod, D. J., McLure, R. J., et al. 2023a, arXiv e-prints, arXiv:2311.06209
- Cullen et al. (2017) Cullen, F., McLure, R. J., Khochfar, S., Dunlop, J. S., & Dalla Vecchia, C. 2017, MNRAS, 470, 3006
- Cullen et al. (2023b) Cullen, F., McLure, R. J., McLeod, D. J., et al. 2023b, MNRAS, 520, 14
- Cullen et al. (2021) Cullen, F., Shapley, A. E., McLure, R. J., et al. 2021, MNRAS, 505, 903
- Curti et al. (2023a) Curti, M., D’Eugenio, F., Carniani, S., et al. 2023a, MNRAS, 518, 425
- Curti et al. (2023b) Curti, M., Maiolino, R., Curtis-Lake, E., et al. 2023b, arXiv e-prints, arXiv:2304.08516
- Curti et al. (2020) Curti, M., Mannucci, F., Cresci, G., & Maiolino, R. 2020, MNRAS, 491, 944
- Curtis-Lake et al. (2023) Curtis-Lake, E., Carniani, S., Cameron, A., et al. 2023, Nature Astronomy, 7, 622
- Davies et al. (2018) Davies, F. B., Hennawi, J. F., Bañados, E., et al. 2018, ApJ, 864, 142
- Dayal & Ferrara (2018) Dayal, P. & Ferrara, A. 2018, Phys. Rep, 780, 1
- Dayal et al. (2011) Dayal, P., Maselli, A., & Ferrara, A. 2011, MNRAS, 410, 830
- Dayal et al. (2024) Dayal, P., Volonteri, M., Greene, J. E., et al. 2024, arXiv e-prints, arXiv:2401.11242
- D’Eugenio et al. (2023) D’Eugenio, F., Maiolino, R., Carniani, S., et al. 2023, arXiv e-prints, arXiv:2311.09908
- Dewdney et al. (2009) Dewdney, P. E., Hall, P. J., Schilizzi, R. T., & Lazio, T. J. L. W. 2009, IEEE Proceedings, 97, 1482
- Dijkstra et al. (2016) Dijkstra, M., Gronke, M., & Venkatesan, A. 2016, ApJ, 828, 71
- D’Odorico et al. (2018) D’Odorico, V., Feruglio, C., Ferrara, A., et al. 2018, ApJ, 863, L29
- Donnan et al. (2023) Donnan, C. T., McLeod, D. J., Dunlop, J. S., et al. 2023, MNRAS, 518, 6011
- Duan et al. (2023) Duan, Q., Conselice, C. J., Li, Q., et al. 2023, arXiv e-prints, arXiv:2309.14961
- Eisenstein et al. (2023) Eisenstein, D. J., Willott, C., Alberts, S., et al. 2023, arXiv e-prints, arXiv:2306.02465
- Endsley et al. (2021) Endsley, R., Stark, D. P., Chevallard, J., & Charlot, S. 2021, MNRAS, 500, 5229
- Endsley et al. (2023) Endsley, R., Stark, D. P., Whitler, L., et al. 2023, MNRAS, 524, 2312
- Fan et al. (2023) Fan, X., Bañados, E., & Simcoe, R. A. 2023, ARA&A, 61, 373
- Ferrara et al. (2023) Ferrara, A., Pallottini, A., & Dayal, P. 2023, MNRAS, 522, 3986
- Finkelstein et al. (2023a) Finkelstein, S. L., Bagley, M. B., Ferguson, H. C., et al. 2023a, ApJ, 946, L13
- Finkelstein et al. (2019) Finkelstein, S. L., D’Aloisio, A., Paardekooper, J.-P., et al. 2019, ApJ, 879, 36
- Finkelstein et al. (2023b) Finkelstein, S. L., Leung, G. C. K., Bagley, M. B., et al. 2023b, arXiv e-prints, arXiv:2311.04279
- Finlator et al. (2018) Finlator, K., Keating, L., Oppenheimer, B. D., Davé, R., & Zackrisson, E. 2018, MNRAS, 480, 2628
- Flury et al. (2022) Flury, S. R., Jaskot, A. E., Ferguson, H. C., et al. 2022, ApJ, 930, 126
- Franco et al. (2023) Franco, M., Akins, H. B., Casey, C. M., et al. 2023, arXiv e-prints, arXiv:2308.00751
- Fujimoto et al. (2023a) Fujimoto, S., Arrabal Haro, P., Dickinson, M., et al. 2023a, ApJ, 949, L25
- Fujimoto et al. (2023b) Fujimoto, S., Wang, B., Weaver, J., et al. 2023b, arXiv e-prints, arXiv:2308.11609
- Fynbo et al. (2009) Fynbo, J. P. U., Jakobsson, P., Prochaska, J. X., et al. 2009, ApJS, 185, 526
- Gonzaga et al. (2012) Gonzaga, S., Hack, W., Fruchter, A., & Mack, J. 2012, The DrizzlePac Handbook
- Greig et al. (2019) Greig, B., Mesinger, A., & Bañados, E. 2019, MNRAS, 484, 5094
- Greig et al. (2017) Greig, B., Mesinger, A., Haiman, Z., & Simcoe, R. A. 2017, MNRAS, 466, 4239
- Gronke et al. (2017) Gronke, M., Dijkstra, M., McCourt, M., & Oh, S. P. 2017, A&A, 607, A71
- Hainline et al. (2023) Hainline, K. N., Johnson, B. D., Robertson, B., et al. 2023, arXiv e-prints, arXiv:2306.02468
- Harikane et al. (2023) Harikane, Y., Ouchi, M., Oguri, M., et al. 2023, ApJS, 265, 5
- Harris et al. (2020) Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Nature, 585, 357
- Hashimoto et al. (2018) Hashimoto, T., Laporte, N., Mawatari, K., et al. 2018, Nature, 557, 392
- Haslbauer et al. (2022) Haslbauer, M., Kroupa, P., Zonoozi, A. H., & Haghi, H. 2022, ApJ, 939, L31
- Hayes et al. (2023) Hayes, M. J., Runnholm, A., Scarlata, C., Gronke, M., & Rivera-Thorsen, T. E. 2023, MNRAS, 520, 5903
- Hayes & Scarlata (2023) Hayes, M. J. & Scarlata, C. 2023, ApJ, 954, L14
- Heintz et al. (2023a) Heintz, K. E., Brammer, G. B., Giménez-Arteaga, C., et al. 2023a, Nature Astronomy
- Heintz et al. (2023b) Heintz, K. E., De Cia, A., Thöne, C. C., et al. 2023b, arXiv e-prints, arXiv:2308.14812
- Heintz et al. (2023c) Heintz, K. E., Giménez-Arteaga, C., Fujimoto, S., et al. 2023c, ApJ, 944, L30
- Heintz et al. (2022) Heintz, K. E., Oesch, P. A., Aravena, M., et al. 2022, ApJ, 934, L27
- Heintz et al. (2023d) Heintz, K. E., Watson, D., Brammer, G., et al. 2023d, arXiv e-prints, arXiv:2306.00647
- Heintz et al. (2021) Heintz, K. E., Watson, D., Oesch, P. A., Narayanan, D., & Madden, S. C. 2021, ApJ, 922, 147
- Henry et al. (2015) Henry, A., Scarlata, C., Martin, C. L., & Erb, D. 2015, ApJ, 809, 19
- Horne (1986) Horne, K. 1986, PASP, 98, 609
- Hoyle (1954) Hoyle, F. 1954, ApJS, 1, 121
- Hsiao et al. (2023) Hsiao, T. Y.-Y., Abdurro’uf, Coe, D., et al. 2023, arXiv e-prints, arXiv:2305.03042
- Hu et al. (2023) Hu, W., Martin, C. L., Gronke, M., et al. 2023, ApJ, 956, 39
- Hunter (2007) Hunter, J. D. 2007, Computing in Science and Engineering, 9, 90
- Hutter et al. (2023) Hutter, A., Trebitsch, M., Dayal, P., et al. 2023, MNRAS, 524, 6124
- Inayoshi et al. (2022) Inayoshi, K., Harikane, Y., Inoue, A. K., Li, W., & Ho, L. C. 2022, ApJ, 938, L10
- Izotov et al. (2006) Izotov, Y. I., Stasińska, G., Meynet, G., Guseva, N. G., & Thuan, T. X. 2006, A&A, 448, 955
- Izotov et al. (2018) Izotov, Y. I., Worseck, G., Schaerer, D., et al. 2018, MNRAS, 478, 4851
- Jakobsen et al. (2022) Jakobsen, P., Ferruit, P., Alves de Oliveira, C., et al. 2022, A&A, 661, A80
- Jakobsson et al. (2006) Jakobsson, P., Fynbo, J. P. U., Ledoux, C., et al. 2006, A&A, 460, L13
- Jaskot et al. (2019) Jaskot, A. E., Dowd, T., Oey, M. S., Scarlata, C., & McKinney, J. 2019, ApJ, 885, 96
- Jiang et al. (2023) Jiang, H., Wang, X., Cheng, C., et al. 2023, arXiv e-prints, arXiv:2312.04151
- Jones et al. (2023) Jones, G. C., Bunker, A. J., Saxena, A., et al. 2023, arXiv e-prints, arXiv:2306.02471
- Jung et al. (2023) Jung, I., Finkelstein, S. L., Arrabal Haro, P., et al. 2023, arXiv e-prints, arXiv:2304.05385
- Katz et al. (2023) Katz, H., Saxena, A., Cameron, A. J., et al. 2023, MNRAS, 518, 592
- Keating et al. (2023a) Keating, L. C., Bolton, J. S., Cullen, F., et al. 2023a, arXiv e-prints, arXiv:2308.05800
- Keating et al. (2023b) Keating, L. C., Puchwein, E., Bolton, J. S., Haehnelt, M. G., & Kulkarni, G. 2023b, arXiv e-prints, arXiv:2308.11709
- Keating et al. (2020) Keating, L. C., Weinberger, L. H., Kulkarni, G., et al. 2020, MNRAS, 491, 1736
- Kennicutt (1998) Kennicutt, Robert C., J. 1998, ARA&A, 36, 189
- Kennicutt & Evans (2012) Kennicutt, R. C. & Evans, N. J. 2012, ARA&A, 50, 531
- Kereš et al. (2005) Kereš, D., Katz, N., Weinberg, D. H., & Davé, R. 2005, MNRAS, 363, 2
- Kewley et al. (2019) Kewley, L. J., Nicholls, D. C., & Sutherland, R. S. 2019, ARA&A, 57, 511
- Khostovan et al. (2016) Khostovan, A. A., Sobral, D., Mobasher, B., et al. 2016, MNRAS, 463, 2363
- Kroupa (2001) Kroupa, P. 2001, MNRAS, 322, 231
- Labbé et al. (2023) Labbé, I., van Dokkum, P., Nelson, E., et al. 2023, Nature, 616, 266
- Langeroodi & Hjorth (2023) Langeroodi, D. & Hjorth, J. 2023, arXiv e-prints, arXiv:2307.06336
- Langeroodi et al. (2023) Langeroodi, D., Hjorth, J., Chen, W., et al. 2023, ApJ, 957, 39
- Lara-López et al. (2010) Lara-López, M. A., Cepa, J., Bongiovanni, A., et al. 2010, A&A, 521, L53
- Laursen et al. (2019) Laursen, P., Sommer-Larsen, J., Milvang-Jensen, B., Fynbo, J. P. U., & Razoumov, A. O. 2019, A&A, 627, A84
- Laursen et al. (2011) Laursen, P., Sommer-Larsen, J., & Razoumov, A. O. 2011, ApJ, 728, 52
- Liang et al. (2024) Liang, L., Feldmann, R., Murray, N., et al. 2024, MNRAS, 528, 499
- Lin et al. (2023) Lin, X., Cai, Z., Zou, S., et al. 2023, ApJ, 944, L59
- Lin et al. (2024) Lin, Y.-H., Scarlata, C., Williams, H., et al. 2024, MNRAS, 527, 4173
- Madau & Dickinson (2014) Madau, P. & Dickinson, M. 2014, ARA&A, 52, 415
- Maiolino & Mannucci (2019) Maiolino, R. & Mannucci, F. 2019, A&A Rev., 27, 3
- Malkan et al. (2017) Malkan, M. A., Cohen, D. P., Maruyama, M., et al. 2017, ApJ, 850, 5
- Mannucci et al. (2010) Mannucci, F., Cresci, G., Maiolino, R., Marconi, A., & Gnerucci, A. 2010, MNRAS, 408, 2115
- Mason et al. (2019) Mason, C. A., Fontana, A., Treu, T., et al. 2019, MNRAS, 485, 3947
- Mason et al. (2023) Mason, C. A., Trenti, M., & Treu, T. 2023, MNRAS, 521, 497
- Mason et al. (2018) Mason, C. A., Treu, T., Dijkstra, M., et al. 2018, ApJ, 856, 2
- Matthee et al. (2023) Matthee, J., Mackenzie, R., Simcoe, R. A., et al. 2023, ApJ, 950, 67
- Matthee et al. (2018) Matthee, J., Sobral, D., Gronke, M., et al. 2018, A&A, 619, A136
- McLeod et al. (2024) McLeod, D. J., Donnan, C. T., McLure, R. J., et al. 2024, MNRAS, 527, 5004
- McQuinn et al. (2007) McQuinn, M., Hernquist, L., Zaldarriaga, M., & Dutta, S. 2007, MNRAS, 381, 75
- McQuinn et al. (2008) McQuinn, M., Lidz, A., Zaldarriaga, M., Hernquist, L., & Dutta, S. 2008, MNRAS, 388, 1101
- Mesinger & Haiman (2007) Mesinger, A. & Haiman, Z. 2007, ApJ, 660, 923
- Meurer et al. (1999) Meurer, G. R., Heckman, T. M., & Calzetti, D. 1999, ApJ, 521, 64
- Miralda-Escudé (1998) Miralda-Escudé, J. 1998, ApJ, 501, 15
- Mirocha & Furlanetto (2023) Mirocha, J. & Furlanetto, S. R. 2023, MNRAS, 519, 843
- Naidu et al. (2022) Naidu, R. P., Oesch, P. A., Setton, D. J., et al. 2022, arXiv e-prints, arXiv:2208.02794
- Naidu et al. (2020) Naidu, R. P., Tacchella, S., Mason, C. A., et al. 2020, ApJ, 892, 109
- Nakajima et al. (2023) Nakajima, K., Ouchi, M., Isobe, Y., et al. 2023, arXiv e-prints, arXiv:2301.12825
- Nakane et al. (2023) Nakane, M., Ouchi, M., Nakajima, K., et al. 2023, arXiv e-prints, arXiv:2312.06804
- Newville et al. (2014) Newville, M., Stensitzki, T., Allen, D. B., & Ingargiola, A. 2014, LMFIT: Non-Linear Least-Square Minimization and Curve-Fitting for Python
- Noterdaeme et al. (2012) Noterdaeme, P., Petitjean, P., Carithers, W. C., et al. 2012, A&A, 547, L1
- Ocvirk et al. (2021) Ocvirk, P., Lewis, J. S. W., Gillet, N., et al. 2021, MNRAS, 507, 6108
- Osterbrock & Ferland (2006) Osterbrock, D. E. & Ferland, G. J. 2006, Astrophysics of gaseous nebulae and active galactic nuclei
- Ouchi et al. (2010) Ouchi, M., Shimasaku, K., Furusawa, H., et al. 2010, ApJ, 723, 869
- Pacucci et al. (2022) Pacucci, F., Dayal, P., Harikane, Y., Inoue, A. K., & Loeb, A. 2022, MNRAS, 514, L6
- Pahl et al. (2024) Pahl, A. J., Shapley, A. E., Steidel, C. C., et al. 2024, arXiv e-prints, arXiv:2401.09526
- Pallottini et al. (2022) Pallottini, A., Ferrara, A., Gallerani, S., et al. 2022, MNRAS, 513, 5621
- Péroux & Howk (2020) Péroux, C. & Howk, J. C. 2020, ARA&A, 58, 363
- Pettini et al. (2000) Pettini, M., Steidel, C. C., Adelberger, K. L., Dickinson, M., & Giavalisco, M. 2000, ApJ, 528, 96
- Planck Collaboration et al. (2020) Planck Collaboration, Aghanim, N., Akrami, Y., et al. 2020, A&A, 641, A6
- Prochaska et al. (2007) Prochaska, J. X., Chen, H.-W., Dessauges-Zavadsky, M., & Bloom, J. S. 2007, ApJ, 666, 267
- Prochaska & Wolfe (2009) Prochaska, J. X. & Wolfe, A. M. 2009, ApJ, 696, 1543
- Puchwein et al. (2023) Puchwein, E., Bolton, J. S., Keating, L. C., et al. 2023, MNRAS, 519, 6162
- Rasmussen Cueto et al. (2023) Rasmussen Cueto, E., Hutter, A., Dayal, P., et al. 2023, arXiv e-prints, arXiv:2312.12109
- Reddy et al. (2018a) Reddy, N. A., Oesch, P. A., Bouwens, R. J., et al. 2018a, ApJ, 853, 56
- Reddy et al. (2018b) Reddy, N. A., Shapley, A. E., Sanders, R. L., et al. 2018b, ApJ, 869, 92
- Rhoads et al. (2023) Rhoads, J. E., Wold, I. G. B., Harish, S., et al. 2023, ApJ, 942, L14
- Rigby et al. (2023) Rigby, J., Perrin, M., McElwain, M., et al. 2023, PASP, 135, 048001
- Roberts-Borsani et al. (2023) Roberts-Borsani, G., Treu, T., Chen, W., et al. 2023, Nature, 618, 480
- Roberts-Borsani et al. (2024) Roberts-Borsani, G., Treu, T., Shapley, A., et al. 2024, arXiv e-prints, arXiv:2403.07103
- Robertson (2022) Robertson, B. E. 2022, ARA&A, 60, 121
- Robertson et al. (2013) Robertson, B. E., Furlanetto, S. R., Schneider, E., et al. 2013, ApJ, 768, 71
- Rosdahl et al. (2022) Rosdahl, J., Blaizot, J., Katz, H., et al. 2022, MNRAS, 515, 2386
- Salmon et al. (2015) Salmon, B., Papovich, C., Finkelstein, S. L., et al. 2015, ApJ, 799, 183
- Sanders et al. (2021) Sanders, R. L., Shapley, A. E., Jones, T., et al. 2021, ApJ, 914, 19
- Sanders et al. (2024) Sanders, R. L., Shapley, A. E., Topping, M. W., Reddy, N. A., & Brammer, G. B. 2024, ApJ, 962, 24
- Sandles et al. (2022) Sandles, L., Curtis-Lake, E., Charlot, S., Chevallard, J., & Maiolino, R. 2022, MNRAS, 515, 2951
- Saxena et al. (2023a) Saxena, A., Bunker, A. J., Jones, G. C., et al. 2023a, arXiv e-prints, arXiv:2306.04536
- Saxena et al. (2023b) Saxena, A., Robertson, B. E., Bunker, A. J., et al. 2023b, A&A, 678, A68
- Schaerer (2003) Schaerer, D. 2003, A&A, 397, 527
- Schaerer et al. (2022) Schaerer, D., Marques-Chaves, R., Barrufet, L., et al. 2022, A&A, 665, L4
- Schaye et al. (2010) Schaye, J., Dalla Vecchia, C., Booth, C. M., et al. 2010, MNRAS, 402, 1536
- Serjeant & Bakx (2023) Serjeant, S. & Bakx, T. J. L. C. 2023, Nature Astronomy, 7, 1143
- Shapley et al. (2023) Shapley, A. E., Sanders, R. L., Reddy, N. A., Topping, M. W., & Brammer, G. B. 2023, ApJ, 954, 157
- Shapley et al. (2003) Shapley, A. E., Steidel, C. C., Pettini, M., & Adelberger, K. L. 2003, ApJ, 588, 65
- Shen et al. (2023) Shen, X., Vogelsberger, M., Boylan-Kolchin, M., Tacchella, S., & Kannan, R. 2023, MNRAS, 525, 3254
- Simmonds et al. (2023) Simmonds, C., Tacchella, S., Maseda, M., et al. 2023, MNRAS, 523, 5468
- Sobacchi & Mesinger (2015) Sobacchi, E. & Mesinger, A. 2015, MNRAS, 453, 1843
- Speagle et al. (2014) Speagle, J. S., Steinhardt, C. L., Capak, P. L., & Silverman, J. D. 2014, ApJS, 214, 15
- Stark (2016) Stark, D. P. 2016, ARA&A, 54, 761
- Steidel et al. (2018) Steidel, C. C., Bogosavljević, M., Shapley, A. E., et al. 2018, ApJ, 869, 123
- Steidel et al. (2010) Steidel, C. C., Erb, D. K., Shapley, A. E., et al. 2010, ApJ, 717, 289
- Steidel et al. (2016) Steidel, C. C., Strom, A. L., Pettini, M., et al. 2016, ApJ, 826, 159
- Steinhardt et al. (2023) Steinhardt, C. L., Kokorev, V., Rusakov, V., Garcia, E., & Sneppen, A. 2023, ApJ, 951, L40
- Sun et al. (2023) Sun, G., Faucher-Giguère, C.-A., Hayward, C. C., et al. 2023, ApJ, 955, L35
- Tang et al. (2023) Tang, M., Stark, D. P., Chen, Z., et al. 2023, MNRAS, 526, 1657
- Tanvir et al. (2019) Tanvir, N. R., Fynbo, J. P. U., de Ugarte Postigo, A., et al. 2019, MNRAS, 483, 5380
- Taylor et al. (2022) Taylor, A. J., Barger, A. J., & Cowie, L. L. 2022, ApJ, 939, L3
- Tepper-García (2006) Tepper-García, T. 2006, MNRAS, 369, 2025
- Thorne et al. (2021) Thorne, J. E., Robotham, A. S. G., Davies, L. J. M., et al. 2021, MNRAS, 505, 540
- Topping et al. (2023) Topping, M. W., Stark, D. P., Endsley, R., et al. 2023, arXiv e-prints, arXiv:2307.08835
- Torrey et al. (2018) Torrey, P., Vogelsberger, M., Hernquist, L., et al. 2018, MNRAS, 477, L16
- Totani et al. (2006) Totani, T., Kawai, N., Kosugi, G., et al. 2006, PASJ, 58, 485
- Treu et al. (2022) Treu, T., Roberts-Borsani, G., Bradac, M., et al. 2022, ApJ, 935, 110
- Trinca et al. (2023) Trinca, A., Schneider, R., Valiante, R., et al. 2023, arXiv e-prints, arXiv:2305.04944
- Trump et al. (2023) Trump, J. R., Arrabal Haro, P., Simons, R. C., et al. 2023, ApJ, 945, 35
- Tumlinson et al. (2017) Tumlinson, J., Peeples, M. S., & Werk, J. K. 2017, ARA&A, 55, 389
- Umeda et al. (2023) Umeda, H., Ouchi, M., Nakajima, K., et al. 2023, arXiv e-prints, arXiv:2306.00487
- Valentino et al. (2023) Valentino, F., Brammer, G., Gould, K. M. L., et al. 2023, ApJ, 947, 20
- van der Wel et al. (2014) van der Wel, A., Franx, M., van Dokkum, P. G., et al. 2014, ApJ, 788, 28
- Ďurovčíková et al. (2024) Ďurovčíková, D., Eilers, A.-C., Chen, H., et al. 2024, arXiv e-prints, arXiv:2401.10328
- Ďurovčíková et al. (2020) Ďurovčíková, D., Katz, H., Bosman, S. E. I., et al. 2020, MNRAS, 493, 4256
- Verhamme et al. (2015) Verhamme, A., Orlitová, I., Schaerer, D., & Hayes, M. 2015, A&A, 578, A7
- Verhamme et al. (2017) Verhamme, A., Orlitová, I., Schaerer, D., et al. 2017, A&A, 597, A13
- Vijayan et al. (2024) Vijayan, A. P., Thomas, P. A., Lovell, C. C., et al. 2024, MNRAS, 527, 7337
- Virtanen et al. (2020) Virtanen, P., Gommers, R., Oliphant, T. E., et al. 2020, Nature Methods, 17, 261
- Vizgan et al. (2022) Vizgan, D., Heintz, K. E., Greve, T. R., et al. 2022, ApJ, 939, L1
- Wang et al. (2023a) Wang, B., Fujimoto, S., Labbe, I., et al. 2023a, arXiv e-prints, arXiv:2308.03745
- Wang et al. (2023b) Wang, B., Leja, J., Labbé, I., et al. 2023b, arXiv e-prints, arXiv:2310.01276
- Wang et al. (2021) Wang, F., Yang, J., Fan, X., et al. 2021, ApJ, 907, L1
- Ward et al. (2023) Ward, E. M., de la Vega, A., Mobasher, B., et al. 2023, arXiv e-prints, arXiv:2311.02162
- Whitaker et al. (2012) Whitaker, K. E., van Dokkum, P. G., Brammer, G., & Franx, M. 2012, ApJ, 754, L29
- White & Frenk (1991) White, S. D. M. & Frenk, C. S. 1991, ApJ, 379, 52
- White & Rees (1978) White, S. D. M. & Rees, M. J. 1978, MNRAS, 183, 341
- Wilson et al. (2023) Wilson, S. N., Heintz, K. E., Jakobsson, P., et al. 2023, arXiv e-prints, arXiv:2305.05213
- Witstok et al. (2023) Witstok, J., Smit, R., Saxena, A., et al. 2023, arXiv e-prints, arXiv:2306.04627
- Witten et al. (2024) Witten, C., Laporte, N., Martin-Alvarez, S., et al. 2024, Nature Astronomy
- Wolfe et al. (2005) Wolfe, A. M., Gawiser, E., & Prochaska, J. X. 2005, ARA&A, 43, 861
- Wolfe et al. (1986) Wolfe, A. M., Turnshek, D. A., Smith, H. E., & Cohen, R. D. 1986, ApJS, 61, 249
- Woosley & Weaver (1995) Woosley, S. E. & Weaver, T. A. 1995, ApJS, 101, 181
- Xu et al. (2022) Xu, X., Henry, A., Heckman, T., et al. 2022, ApJ, 933, 202
- Yang et al. (2017) Yang, H., Malhotra, S., Gronke, M., et al. 2017, ApJ, 844, 171
- Yang et al. (2020) Yang, J., Wang, F., Fan, X., et al. 2020, ApJ, 904, 26
- Yeh et al. (2023) Yeh, J. Y. C., Smith, A., Kannan, R., et al. 2023, MNRAS, 520, 2757
- Zackrisson et al. (2011) Zackrisson, E., Rydberg, C.-E., Schaerer, D., Östlin, G., & Tuli, M. 2011, ApJ, 740, 13
- Zavala et al. (2024) Zavala, J. A., Castellano, M., Akins, H. B., et al. 2024, arXiv e-prints, arXiv:2403.10491