JWST spectra revealing Lyman-α (Lyα) absorption in z 5-14 galaxies offer a unique probe of reionization's earliest stages. However, disentangling absorption by the increasingly neutral intergalactic medium (IGM) from that in galaxies' interstellar and circumgalactic medium (ISM, CGM) remains challenging due to the poorly constrained nature of neutral hydrogen (HI) in and around galaxies at these redshifts. We use the SERRA high-resolution cosmological zoom-in simulations to characterize the HI distribution and its evolution along sightlines to star-forming regions during reionization, to interpret the contribution of local HI (ISM+CGM) to Lyα absorption in z>5 spectra. We analyze 100 z=6-9. 5 galaxies, generating mock sightlines from each galaxy's star formation peak. We study the sightline distribution of HI column densities (N ₇₈) and its variation with radius, halo mass and redshift. We find broad sightline variation in N ₇₈ (0. 5-1. 5dex) due to complex ISM morphology driven by bursty star formation, with median (N ₇₈/ cm^-2) 21-22. Dense ISM gas is the dominant origin of damped Lyα absorption (DLA) systems along sightlines towards star-forming regions, outweighing gas in the CGM, filaments and proximate absorbers. Median N ₇₈ increases with halo mass, scaling roughly with the virial radius, as expected due to larger potentials and more extended CGM, but shows negligible redshift evolution at fixed halo mass. This suggests post-reionization N ₇₈ distributions may provide useful priors to interpret IGM damping wings at higher redshifts. To interpret strong z>5 DLA candidates found by JWST, we investigate N ₇₈>10^22 cm^-2 sightlines. These trace dense, metal-enriched ISM within <1 kpc of massive halos, a scenario testable with higher-resolution spectroscopy.
Gelli et al. (Wed,) studied this question.