Within the next few years, the Square Kilometer Array (SKA) or one of its pathfinders will hopefully provide a detection of the 21-cm signal fluctuations from the Epoch of Reionization (EoR). Then, the main goal will be to accurately constrain the underlying astrophysical parameters. Currently, this is mainly done with Bayesian inference using Markov Chain Monte Carlo sampling. Recently, studies using neural networks trained to performed inverse modelling have shown interesting results. We build on these by improving the accuracy of the predictions using neural network and exploring other supervised learning methods: the kernel and ridge regressions. Based on a large training set of 21-cm power spectra, we compare the performances of these supervised learning methods. When using an un-noised signal as input, we improve on previous neural network accuracy by one order of magnitude and, using local ridge kernel regression, we gain another factor of a few. We then reach a rms prediction error of a few percents of the 1-sigma confidence level due to SKA thermal noise (as estimated with Bayesian inference). This last performance level requires optimizing the hyper-parameters of the method: how to do that perfectly in the case of an unknown signal remains an open question. For an input signal altered by a SKA-type thermal noise, our neural network recovers the astrophysical parameter values with an error within half of the 1$\sigma$ confidence level due to the SKA thermal noise. This accuracy improves to 10$\%$ of the 1$\sigma$ level when using the local ridge kernel regression (with optimized hyper-parameters). We are thus reaching a performance level where supervised learning methods are a viable alternative to determine the best-fit parameters values.

With a statistical detection of the 21 cm signal fluctuations from the Epoch of Reionization (EoR) expected in the next few years, there is an interest in developing robust and precise techniques to constrain the underlying astrophysical parameters. Bayesian inference with Markov Chain Monte Carlo, or different types of supervised learning for backward modelling (from signal to parameters) are examples of such techniques. They usually require many instances of forward modelling (from parameters to signal) in sampling the parameters space, either when performing the steps of the Markov Chain or when building a training sample for supervised learning. As forward modelling can be costly (if performed with numerical simulations for example), we should attempt to perform an optimal sampling according to some principle. With this goal in mind, we present an approach based on defining a metric on the space of observables, induced by the manner through which the modelling creates a mapping from the parameter space onto the space of observables. This metric bears a close connection to Jeffreys' prior from information theory. It is used to generate a homogeneous and isotropic sampling of the signal space with two different methods. We show that when the resulting optimized samplings, created with 21cmFAST, are used to train a neural network we obtain a modest reduction of the error on parameter reconstruction of ~10% (compared to a na\"ive sampling of the same size). Excluding the borders of the parameter space region, the improvement is more substantial, on the order of 30-40%.

Quasars contribute to the 21-cm signal from the Epoch of Reionization (EoR) primarily through their ionizing UV and X-ray emission. However, their radio continuum and Lyman-band emission also regulates the 21-cm signal in their direct environment, potentially leaving the imprint of their duty cycle. We develop a model for the radio and UV luminosity functions of quasars from the EoR, and constrain it using recent observations. Our model is consistent with the z~7.5 quasar from Banados et al 2017, and also predicts only a few quasars suitable for 21-cm forest observations (10mJy) in the sky. We exhibit a new effect on the 21-cm signal observed against the CMB: a radio-loud quasar can leave the imprint of its duty cycle on the 21-cm tomography. We apply this effect in a cosmological simulation and conclude that the effect of typical radio-loud quasars is most likely negligible in an SKA field of view. For a 1-10mJy quasar the effect is stronger though hardly observable at SKA resolution. Then we study the contribution of the lyman band Ly-alpha to Ly-beta) emission of quasars to the Wouthuisen-Field coupling. The collective effect of quasars on the 21-cm power spectrum is larger than the thermal noise at low k, though featureless. However, a distinctive pattern around the brightest quasars in an SKA field of view may be observable in the tomography, encoding the duration of their duty cycle. This pattern has a high signal-to-noise ratio for the brightest quasar in a typical SKA shallow survey.

The 21-cm signal from the Epoch of Reionization (EoR) is expected to be detected in the next few years, either with existing instruments or by the upcoming SKA and HERA projects. In this context there is a pressing need for publicly available high-quality templates covering a wide range of possible signals. These are needed both for end-to-end simulations of the up-coming instruments, as well as to develop signal analysis methods. In this work we present such a set of templates, publicly available, for download at https://21ssd.obspm.fr/. The database contains 21-cm brightness temperature lightcones at high and low resolution, and several derived statistical quantities for 45 models spanning our choice of 3D parameter space. These data are the result of fully coupled radiative hydrodynamic high resolution ($1024^3$) simulations performed with the LICORICE code. Both X-ray and Lyman line transfer is performed to account for heating and Wouthuysen-Field coupling fluctuations. We also present a first exploitation of the data using the power spectrum and the Pixel Distribution Function (PDF) as functions of redshifts, computed from lightcone data. We analyse how these two quantities behave when varying the model parameters while taking into account the thermal noise expected of a typical SKA survey. Finally, we show that the power spectrum and the PDF have different -- and to some extent complementary -- abilities to distinguish between different models. This opens the door to formulating an optimal sampling of the parameter space, dependant on the chosen diagnostics.