In the Simulations and Constructions of the Reionization of Cosmic Hydrogen (SCORCH) project, we compare analytical models of the hydrogen ionization fraction with radiation-hydrodynamic simulations. We derive analytical models of the mass-weighted hydrogen ionization fraction from the local ionization balance equations as a more accurate alternative to the widely adopted model based on the volume filling factor. In particular, our model has a recombination term quadratic in the ionization fraction, which is consistent with the two-body interaction nature of recombination. Then, we use the radiation-hydrodynamic simulations to study the clumping factors needed to solve the analytical equations, and provide accurate fitting functions. We find that the ionized hydrogen clumping factors from our radiative transfer simulations are significantly different than those from other simulations that use a uniform photoionization background. In addition to redshift dependence, we also see the dependence of ionized hydrogen clumping factor on ionization fraction, and we incorporate this into our fits. We calculate the reionization histories using our analytical models and clumping factors and compare with widely adopted models, and all of our models achieve $<7\%$ difference from simulation results while the other models have $>20\%$ deviations. The Thomson optical depths from reionization calculated from our analytical models result in $<5\%$ deviation from simulations, while the previous analytical models have $>20\%$ difference in and could result in biased conclusions of the IGM reionization.

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.

In the Simulations and Constructions of the Reionization of Cosmic Hydrogen (SCORCH) project, we present new radiation-hydrodynamic simulations with updated high-redshift galaxy populations and varying radiation escape fractions. The simulations are designed to have fixed Thomson optical depth $\tau \approx 0.06$, consistent with recent Planck observations, and similar midpoints of reionization $7.5 \lesssim z \lesssim 8.0$, but with different ionization histories. The galaxy luminosity functions and ionizing photon production rates in our model are in good agreement with recent HST observations. Adopting a power-law form for the radiation escape fraction $f_{\text{esc}}(z) = f_8[(1+z)/9]^{a_{8}}$, we simulate the cases for $a_8 = 0$, 1, and 2 and find $a_8 \lesssim 2$ in order to end reionization in the range $5.5 \lesssim z \lesssim 6.5$ to be consistent with Lyman alpha forest observations. At fixed $\tau$ and as the power-law slope $a_8$ increases, the reionization process starts earlier but ends later with a longer duration $\Delta z$ and the decreased redshift asymmetry $Az$. We find a range of durations $3.9 \lesssim \Delta z \lesssim 4.6$ that is currently in tension with the upper limit $\Delta z < 2.8$ inferred from a recent joint analysis of Planck and South Pole Telescope observations.

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%.