2022ApJ...926...31R

We perform joint modeling of the composite rest-frame far-UV and optical spectra of redshift 1.85 <= z <= 3.49 star-forming galaxies to deduce key properties of the massive stars, ionized interstellar medium (ISM), and neutral ISM, with the aim of investigating the principal factors affecting the production and escape of Lyα photons. Our sample consists of 136 galaxies with deep Keck/LRIS and MOSFIRE spectra covering, respectively, Lyβ through C III] λλ1907, 1909 and [O II], [Ne III], Hβ, [O III], Hα, [N II], and [S II]. Spectral and photoionization modeling indicates that the galaxies are uniformly consistent with stellar population synthesis models that include the effects of stellar binarity. Over the dynamic range of our sample, there is little variation in stellar and nebular abundance with Lyα equivalent width, W _ λ _(Lyα), and only a marginal anticorrelation between age and W _ λ _(Lyα). The inferred range of ionizing spectral shapes is insufficient to solely account for the variation in W _ λ _(Lyα); rather, the covering fraction of optically thick H I appears to be the principal factor modulating the escape of Lyα, with most of the Lyα photons in down-the-barrel observations of galaxies escaping through low column density or ionized channels in the ISM. Our analysis shows that a high star-formation-rate surface density, Σ_SFR, particularly when coupled with a low galaxy potential (i.e., low stellar mass), can aid in reducing the covering fraction and ease the escape of Lyα photons. We conclude with a discussion of the implications of our results for the escape of ionizing radiation at high redshift.