The Effects of Stellar Population and Gas Covering Fraction on the Emergent Ly alpha Emission of High-redshift Galaxies*

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 alpha photons. Our sample consists of 136 galaxies with deep Keck/LRIS and MOSFIRE spectra covering, respectively, Ly beta through C iii] lambda lambda 1907, 1909 and [O ii], [Ne iii], H beta, [O iii], H alpha, [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 alpha equivalent width, W ( lambda )(Ly alpha), and only a marginal anticorrelation between age and W ( lambda )(Ly alpha). The inferred range of ionizing spectral shapes is insufficient to solely account for the variation in W ( lambda )(Ly alpha); rather, the covering fraction of optically thick H i appears to be the principal factor modulating the escape of Ly alpha, with most of the Ly alpha 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, sigma(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 alpha photons. We conclude with a discussion of the implications of our results for the escape of ionizing radiation at high redshift.