Predicting Lyman-continuum emission of galaxies using their physical and Lyman-alpha 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 (Ly alpha) 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 Ly alpha radiative transfer code RASCAS and analyze the Ly alpha - LyC connection. Results. We find that the Ly alpha and LyC luminosities are strongly correlated with each other, although with dispersion. There is a positive correlation between the escape fractions of Ly alpha and LyC radiations in the brightest Lyman-alpha emitters (LAEs; escaping Ly alpha luminosity L-esc(Ly alpha) > 10(41) erg s(-1)), 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-esc(Ly alpha) > 10(40) erg s(-1) 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 Ly alpha 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 R-2 of these predictions in our fiducial model are 0.89 and 0.85, respectively, where R-2 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 L-esc(Ly alpha), 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 Ly alpha properties and can thus help us identify the sources of reionization.