Extreme equivalent width-selected low-mass starbursts at $z=4-9$: insights into their role in cosmic reionization

A. Calabr\`o; A. Ferrara; A. J. Taylor; A. M. Koekemoer; C. Papovich; D. J. McLeod; E. McGrath; F. Arevalo-Gonzalez; J. S. Kartaltepe; K. Davis

arxiv: 2510.25647 · v2 · submitted 2025-10-29 · 🌌 astro-ph.GA

Extreme equivalent width-selected low-mass starbursts at z=4-9: insights into their role in cosmic reionization

M. Llerena , L. Pentericci , R. Amor\'in , A. Ferrara , M. Dickinson , F. Arevalo-Gonzalez , A. Calabr\`o , L. Napolitano

show 17 more authors

S. Mascia P. Arrabal Haro R. Begley N. J. Cleri K. Davis W. Hu J. S. Kartaltepe A. M. Koekemoer R. A. Lucas E. McGrath D. J. McLeod C. Papovich T. M. Stanton A. J. Taylor R. Tripodi X. Wang L. Y. A. Yung

This is my paper

Pith reviewed 2026-05-18 03:24 UTC · model grok-4.3

classification 🌌 astro-ph.GA

keywords extreme emission line galaxiescosmic reionizationhigh-redshift galaxiesionizing radiationstar formationJWST spectroscopylow-mass galaxies

0 comments

The pith

Extreme low-mass starbursts at high redshift contribute 16-40% of the ionizing photons for cosmic reionization.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper examines compact galaxies with extremely strong emission lines at redshifts from 4 to 9 using JWST data. These low-mass systems show high specific star formation rates and produce ionizing photons efficiently. Escape of these photons is variable but sufficient on average to make them responsible for a sizable portion of the radiation that reionizes the universe. This highlights their potential importance in transforming the early cosmos from opaque to transparent.

Core claim

Using spectroscopy from JWST surveys, the study finds that extreme emission line galaxies at z=4-9 have median equivalent widths of 1616 Å for [O III] + Hβ and are low-mass with high sSFR. Their efficiency in producing ionizing photons combined with a median LyC escape fraction of 5% leads to them contributing 16-40% of the total ionizing emissivity needed for hydrogen reionization.

What carries the argument

Extreme emission line galaxies (EELGs) selected by their intense rest-frame equivalent widths, acting as indicators of recent star formation bursts that drive both high photon production and potential leakage of ionizing radiation.

If this is right

These galaxies exceed typical ionizing photon production efficiencies at similar redshifts.
Their compactness and high star formation rate surface densities are key to the extreme line strengths.
A small fraction exhibit high escape fractions for both Lyα and LyC photons.
They may represent a key population for understanding the sources of reionization.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

If their contribution holds, reionization models should account more for low-mass star-forming galaxies.
Direct observations of more LyC escape fractions could refine the 16-40% estimate.
Similar selection at lower redshifts might reveal how these systems evolve over time.

Load-bearing premise

That the median 5% LyC escape fraction and the heterogeneous measurements from a subset can be extrapolated to the full population of such galaxies.

What would settle it

Finding that most such galaxies have LyC escape fractions below 1% or a much lower number density than assumed would make their total contribution to reionization negligible.

read the original abstract

We investigate the properties of extreme emission line galaxies (EELGs) at $z=4-9$ and their role in reionization. Compact, low-mass galaxies with intense optical emission lines are linked to elevated specific star formation rates (sSFRs) and recent bursts of star formation. Feedback in these systems may enable the leakage of ionizing radiation into the intergalactic medium. Using JWST/NIRSpec spectroscopy from the CAPERS, CEERS, and RUBIES surveys, we compile 160 NIRCam-selected EELGs in the EGS field. These galaxies show extreme rest-frame equivalent widths (EWs), with a median EW([O III]+H${\beta}$)=1616\r{A} and EW(H${\alpha}$)=763\r{A}. They are low-mass (median log(M$_{\star}$/M$_{\odot}$)=8.26) with high sSFRs (median 43 Gyr$^{-1}$), above the $z\sim6$ main sequence. UV slopes are diverse, with a mean $\beta=-2.0$, and only 7% have extremely blue continua ($\beta<-2.6$). Emission-line diagnostics suggest stellar populations as the primary ionizing source, although an AGN fraction of 14% is cannot be entirely ruled out. These galaxies are efficient ionizing photon producers, with mean log($\xi_{\rm ion}$ [Hz erg$^{-1}$])=25.37, exceeding typical values at similar redshifts. Escape fractions, however, are heterogeneous: 16% of EELGs at $z<7$ show escape fractions $>$5% for both Ly${\alpha}$ and LyC photons, while 82% lack detectable Ly${\alpha}$ emission. The median inferred LyC escape fraction is modest (5%) but enhanced in compact super-Eddington systems with sSFR >25 Gyr$^{-1}$. These results indicate that EELGs contribute approximately 16-40% of the total ionizing emissivity required to sustain hydrogen reionization. EELGs are extremely compact, with a median effective radius of 0.49 kpc, and exhibit a recent star-formation burst. Our analysis indicates that sSFR and star-formation rate surface density are the primary drivers of their extreme emission line strengths.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

A solid new JWST sample of 160 high-z EELGs with direct measurements on EWs, masses, sSFRs and xi_ion, but the 16-40% reionization contribution rests on extrapolating a 5% median LyC escape fraction from limited heterogeneous data.

read the letter

The main thing to know is that this paper assembles a sample of 160 extreme emission line galaxies at z=4-9 from JWST/NIRSpec in the EGS field and reports their basic properties along with an estimate that they supply 16-40% of the ionizing emissivity needed for reionization. The sample size and the direct spectroscopic numbers are the clearest addition here. They measure median EW([O III]+Hβ) of 1616 Å, log M⋆ of 8.26, sSFR of 43 Gyr⁻¹, and log ξ_ion of 25.37, which sits above typical values at these redshifts. The galaxies are compact (median 0.49 kpc) and mostly consistent with stellar ionization, though the 14% AGN fraction is noted without being fully excluded. These measurements give a useful empirical anchor for what low-mass, bursty systems look like at early times. The reionization calculation is the softer link. The median LyC escape fraction of 5% is described as enhanced in compact super-Eddington objects, but it comes from heterogeneous constraints: only 16% of the z<7 subset show >5% escape in both Lyα and LyC, while 82% lack detectable Lyα. Scaling the observed number density by this median and by ξ_ion to reach 16-40% of the required budget therefore depends on how well the measured tail represents the full population. Field variance and selection effects could move that fraction. The paper is aimed at modelers who need inputs on high-redshift low-mass galaxies and their photon output. Readers focused on reionization sources or JWST demographics will find the property distributions and efficiency values worth pulling. The observational core is strong enough that it deserves peer review, even if the escape-fraction extrapolation needs tighter justification from referees.

Referee Report

2 major / 2 minor

Summary. The manuscript compiles a sample of 160 extreme emission line galaxies (EELGs) at z=4-9 from JWST/NIRSpec spectroscopy in the EGS field. It reports extreme rest-frame equivalent widths (median EW([O III]+Hβ)=1616 Å, EW(Hα)=763 Å), low stellar masses (median log(M*/M⊙)=8.26), high sSFRs (median 43 Gyr⁻¹), efficient ionizing photon production (mean log ξ_ion=25.37), heterogeneous LyC escape fractions (median 5%, enhanced in compact super-Eddington systems), and concludes that these galaxies contribute 16-40% of the ionizing emissivity required to sustain hydrogen reionization.

Significance. If the population-averaged LyC escape fraction can be shown to be robustly represented by the reported median of 5%, the work would establish low-mass EELGs as a major contributor to cosmic reionization at z=4-9. The direct spectroscopic measurements of EWs, masses, sSFRs, and ξ_ion provide a solid observational foundation, and the identification of sSFR and star-formation rate surface density as drivers of extreme line strengths is a useful result.

major comments (2)

[Abstract and reionization calculation] The 16-40% reionization contribution (stated in the abstract and final paragraph) is obtained by multiplying the observed EELG number density by mean log ξ_ion=25.37 and a median LyC f_esc=5%. This median is extrapolated from heterogeneous data in which only 16% of the z<7 subsample show f_esc>5% for both Lyα and LyC while 82% lack detectable Lyα; the validity of applying this median to the full population and to the cosmic emissivity budget therefore requires explicit justification of selection effects and field-to-field variance.
[Emission-line diagnostics section] The AGN fraction of 14% is noted but not fully excluded by emission-line diagnostics. Because even a modest AGN contribution would alter the interpretation of ξ_ion as purely stellar and affect the inferred ionizing emissivity, a quantitative assessment of the impact on the reionization fraction (e.g., via upper limits on AGN-driven photons) is needed.

minor comments (2)

[Abstract] Grammatical error in the abstract: 'an AGN fraction of 14% is cannot be entirely ruled out' should read 'an AGN fraction of 14% cannot be entirely ruled out'.
[UV slopes paragraph] The UV-slope distribution and the 7% fraction with β < -2.6 are stated without reference to the exact measurement method or aperture corrections; a brief methods paragraph would improve clarity.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed comments, which have helped clarify several aspects of our analysis. We respond to each major comment below and indicate the revisions we will make to the manuscript.

read point-by-point responses

Referee: [Abstract and reionization calculation] The 16-40% reionization contribution (stated in the abstract and final paragraph) is obtained by multiplying the observed EELG number density by mean log ξ_ion=25.37 and a median LyC f_esc=5%. This median is extrapolated from heterogeneous data in which only 16% of the z<7 subsample show f_esc>5% for both Lyα and LyC while 82% lack detectable Lyα; the validity of applying this median to the full population and to the cosmic emissivity budget therefore requires explicit justification of selection effects and field-to-field variance.

Authors: We appreciate the referee's emphasis on the need for careful justification of the median f_esc. This value is derived from the subset of z<7 EELGs with available Lyα and indirect LyC constraints, where non-detections are assigned conservative upper-limit-based estimates consistent with the observed heterogeneity (only 16% show f_esc >5%). In the revised manuscript we will add an explicit subsection justifying the extrapolation: we will detail the selection function of the NIRSpec sample, the physical basis for assigning modest f_esc to Lyα-undetected systems (supported by their compact morphologies and high sSFR), and present the reionization contribution as a range that brackets plausible variations in the undetected population. For field-to-field variance we will compare our EGS number densities to published values from other JWST fields (e.g., GOODS, COSMOS) and include a brief discussion of cosmic-variance uncertainties while retaining the single-field limitation as a caveat. revision: partial
Referee: [Emission-line diagnostics section] The AGN fraction of 14% is noted but not fully excluded by emission-line diagnostics. Because even a modest AGN contribution would alter the interpretation of ξ_ion as purely stellar and affect the inferred ionizing emissivity, a quantitative assessment of the impact on the reionization fraction (e.g., via upper limits on AGN-driven photons) is needed.

Authors: We agree that a quantitative bound on possible AGN contamination strengthens the interpretation. In the revised version we will add a short calculation that places an upper limit on the AGN-driven ionizing-photon contribution. Assuming the full 14% fraction consists of AGN with ξ_ion values typical of high-redshift AGN (log ξ_ion ≈ 25.5–26.0), we show that the additional emissivity changes the total EELG contribution to reionization by at most a few percentage points, keeping the overall 16–40% range intact. We will also state explicitly that the reported mean log ξ_ion = 25.37 is computed from stellar-population fits after removing the AGN candidates identified by the diagnostics, thereby preserving the stellar interpretation for the majority population. revision: yes

Circularity Check

0 steps flagged

No significant circularity in the reionization contribution estimate

full rationale

The paper's derivation starts from JWST/NIRSpec observations of 160 EELGs, directly measures median EW([O III]+Hβ)=1616 Å, log ξ_ion=25.37, number density in the EGS field, and infers a median f_esc=5% from heterogeneous Lyα/LyC diagnostics on the sample itself. The 16-40% contribution is then obtained by multiplying these independent empirical quantities and comparing the result to an externally supplied total ionizing emissivity requirement for reionization. No equation or step reduces by construction to a prior fit, self-definition, or load-bearing self-citation; the calculation is a standard observational estimate whose inputs are falsifiable outside the final percentage.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The contribution percentage depends on standard astrophysical assumptions about stellar ionizing sources and escape-fraction diagnostics; no new free parameters are explicitly fitted beyond sample medians, and no invented entities are introduced.

free parameters (1)

median LyC escape fraction
Inferred value of 5% used to scale the population contribution to the reionization budget.

axioms (1)

domain assumption Emission-line diagnostics reliably indicate stellar populations as the dominant ionizing source
Invoked to interpret the high ξ_ion and to down-weight AGN contribution.

pith-pipeline@v0.9.0 · 6118 in / 1306 out tokens · 42373 ms · 2026-05-18T03:24:37.595398+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

median log(ξ_ion)=25.34, median f_esc=4.4%, contribution 16-40% of ionizing emissivity

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

An Updated Characterization of Luminous Ly{\alpha} emitters at the End of Reionization
astro-ph.GA 2026-05 unverdicted novelty 5.0

Luminous Lyα emitters at z≈6 are low-mass ultra-young dwarf starbursts with median Lyα escape fractions above 40 percent, driven by vigorous star formation and low dust content.