Dr Harley Katz

Abstract:

We identify a low-metallicity (12 + log(O/H) = 7.59) Ly𝛼-emitting galaxy at 𝑧 = 5.943 with evidence of a strong Balmer jump, arising from nebular continuum. While Balmer jumps are sometimes observed in low-redshift star-forming galaxies, this galaxy also exhibits a steep turnover in the UV continuum. Such turnovers are typically attributed to absorption by a damped Ly𝛼 system (DLA); however, the shape of the turnover and the high observed Ly𝛼 escape fraction ( 𝑓esc,Ly𝛼 ∼ 27%) is also consistent with strong nebular two-photon continuum emission. Modelling the UV turnover with a DLA requires extreme column densities (𝑁HI > 1023 cm−2 ), and simultaneously explaining the high 𝑓esc,Ly𝛼 requires a fine-tuned geometry. In contrast, modelling the spectrum as primarily nebular provides a good fit to both the continuum and emission lines, motivating scenarios in which (a) we are observing only nebular emission or (b) the ionizing source is powering extreme nebular emission that outshines the stellar emission. The nebular-only scenario could arise if the ionising source has ‘turned off’ more recently than the recombination timescale (∼1,000 yr), hence we may be catching the object at a very specific time. Alternatively, hot stars with 𝑇eff ≳ 105 K (e.g. Wolf-Rayet or low-metallicity massive stars) produce enough ionizing photons such that the two-photon emission becomes visible. While several stellar SEDs from the literature fit the observed spectrum well, the hot-star scenario requires that the number of ≳ 50 M⊙ stars relative to ∼ 5 − 50 M⊙ stars is significantly higher than predicted by typical stellar initial mass functions (IMFs). The identification of more galaxies with similar spectra may provide evidence for a top-heavy IMF at high redshift.

Abstract:

Recent studies have reported tension between the presence of luminous, high-redshift galaxies and the halo mass functions predicted by standard cosmology. Here, an improved test is proposed using the presence of high-redshift Balmer breaks to probe the formation of early 104–105 M ⊙ baryonic minihalos. Unlike previous tests, this does not depend upon the mass-to-light ratio and has only a slight dependence upon the metallicity, stellar initial mass function, and star formation history, which are all weakly constrained at high redshift. We show that the strongest Balmer breaks allowed at z = 9 using the simplest ΛCDM cosmological model would allow a D 4000 as high as 1.26 under idealized circumstances and D 4000 ≤ 1.14 including realistic feedback models. Since current photometric template fitting to JWST sources infers the existence of stronger Balmer breaks out to z ≳ 11, upcoming spectroscopic follow-up will either demonstrate those templates are invalid at high redshift or imply new physics beyond “vanilla” ΛCDM.

Abstract:

Reliable indirect diagnostics of LyC photon escape from galaxies are required to understand which sources were the dominant contributors to reionization. While multiple LyC escape fraction (f_esc) indicators have been proposed to trace favourable conditions for LyC leakage from the interstellar medium of low-redshift ‘analogue’ galaxies, it remains unclear whether these are applicable at high redshifts where LyC emission cannot be directly observed. Using a library of 14 120 mock spectra of star-forming galaxies with redshifts 4.64 ≤ z ≤ 10 from the SPHINX²⁰ cosmological radiation hydrodynamics simulation, we develop a framework for the physics that leads to high f_esc. We investigate LyC leakage from our galaxies based on the criteria that successful LyC escape diagnostics must (i) track a high-specific star formation rate, (ii) be sensitive to stellar population age in the range 3.5–10 Myr representing the times when supernova first explode to when LyC production significantly drops, and (iii) include a proxy for neutral gas content and gas density in the interstellar medium. O₃₂, Σ_SFR, M_UV, and H β equivalent width select for one or fewer of our criteria, rendering them either necessary but insufficient or generally poor diagnostics. In contrast, UV slope (β), and E(B − V) match two or more of our criteria, rendering them good f_esc diagnostics (albeit with significant scatter). Using our library, we build a quantitative model for predicting f_esc based on direct observables. When applied to bright z > 6 Ly α emitters observed with JWST, we find that the majority of them have 𝑓esc≲10 per cent⁠.