Predicting LyC emission of galaxies using their physical and Ly? emission properties

Aims. The primary difficulty in understanding the sources and processes that powered cosmic reionization is that it is not possible to directly probe the ionizing Lyman-continuum (LyC) radiation at that epoch as those photons have been absorbed by the intervening neutral hydrogen. It is therefore imperative to build a model to accurately predict LyC emission using other properties of galaxies in the reionization era. Methods. In recent years, studies have shown that the LyC emission from galaxies may be correlated to their Lyman-alpha (LyAlpha) emission. In this paper we study this correlation by analyzing thousands of simulated galaxies at high redshift in the SPHINX cosmological simulation. We post-process these galaxies with the LyAlpha radiative transfer code RASCAS and analyze the LyAlpha - LyC connection. Results. We find that the LyAlpha and LyC luminosities are strongly correlated with each other, although with dispersion. There is a positive correlation between the escape fractions of LyAlpha and LyC radiations in the brightest Lyman-alpha emitters (LAEs; escaping LyAlpha luminosity L^LyAlpha esc > 1041 erg s1), similar to that reported by recent observational studies. However, when we also include fainter LAEs, the correlation disappears, which suggests that the observed relation may be driven by selection effects. We also find that the brighter LAEs are dominant contributors to reionization, with L^LyAlpha_esc > 1040 erg s1 galaxies accounting for >90% of the total amount of LyC radiation escaping into the intergalactic medium in the simulation. Finally, we build predictive models using multivariate linear regression, where we use the physical and LyAlpha properties of simulated reionization era galaxies to predict their LyC emission. We build a set of models using different sets of galaxy properties as input parameters and predict their intrinsic and escaping LyC luminosity with a high degree of accuracy (the adjusted R2 of these predictions in our fiducial model are 0.89 and 0.85, respectively, where R2 is a measure of how much of the response variance is explained by the model). We find that the most important galaxy properties for predicting the escaping LyC luminosity of a galaxy are its LLyAlpha esc , gas mass, gas metallicity, and star formation rate. Conclusions. These results and the predictive models can be useful for predicting the LyC emission from galaxies using their physical and LyAlpha properties and can thus help us identify the sources of reionization.