No Arabic abstract
The reionization process is expected to be prolonged by the small-scale absorbers (SSAs) of ionizing photons, which have been seen as Lyman-limit systems in quasar absorption line observations. We use a set of semi-numerical simulations to investigate the effects of absorption systems on the reionization process, especially their impacts on the neutral islands during the late epoch of reionization (EoR). Three model are studied, i.e. the extreme case of no-SSA model with a high level of ionizing background, the moderate-SSA model with a relatively high level of ionizing background, and the dense-SSA model with a low level of ionizing background. We find that while the characteristic scale of neutral regions decreases during the early and middle stages of reionization, it stays nearly unchanged at about 10 comoving Mpc during the late stage for the no-SSA and moderate-SSA models. However, in the case of weak ionizing background in the dense-SSA model, the characteristic island scale shows obvious evolution, as large islands break into many small ones that are slowly ionized. The evolutionary behavior of neutral islands during the late EoR thus provides a novel way to constrain the abundance of SSAs. We discuss the 21-cm observation with the upcoming Square Kilometre Array (SKA). The different models can be distinguished by the 21-cm power spectrum measurement, and it is also possible to extract the characteristic island scale from the imaging observation with a proper choice of the 21-cm brightness threshold.
The first generation of redshifted 21 cm detection experiments, carried out with arrays like LOFAR, MWA and GMRT, will have a very low signal-to-noise ratio per resolution element (sim 0.2). In addition, whereas the variance of the cosmological signal decreases on scales larger than the typical size of ionization bubbles, the variance of the formidable galactic foregrounds increases, making it hard to disentangle the two on such large scales. The poor sensitivity on small scales on the one hand, and the foregrounds effect on large scales on the other hand, make direct imaging of the Epoch of Reionization of the Universe very difficult, and detection of the signal therefore is expected to be statistical.Despite these hurdles, in this paper we argue that for many reionization scenarios low resolution images could be obtained from the expected data. This is because at the later stages of the process one still finds very large pockets of neutral regions in the IGM, reflecting the clustering of the large-scale structure, which stays strong up to scales of sim 120 comoving Mpc/h (sim 1 degree). The coherence of the emission on those scales allows us to reach sufficient S/N (sim 3) so as to obtain reionization 21 cm images. Such images will be extremely valuable for answering many cosmological questions but above all they will be a very powerful tool to test our control of the systematics in the data. The existence of this typical scale (sim 120 comoving Mpc/h) also argues for designing future EoR experiments, e.g., with SKA, with a field of view of at least 4 degree.
Fractal dimensions can be used to characterize the clustering and lacunarities in density distributions. We use generalized fractal dimensions to study the neutral hydrogen distribution (HI) during the epoch of reionization. Using a semi-numeric model of ionized bubbles to generate the HI field, we calculate the fractal dimensions for length scales $sim 10 h^{-1}$ cMpc. We find that the HI field displays significant multifractal behaviour and is not consistent with homogeneity at these scales when the mass averaged neutral fraction $bar{x}_{rm HI}^M gtrsim 0.5$. This multifractal nature is driven entirely by the shapes and distribution of the ionized regions. The sensitivity of the fractal dimension to the neutral fraction implies that it can be used for constraining reionization history. We find that the fractal dimension is relatively less sensitive to the value of the minimum mass of ionizing haloes when it is in the range $sim 10^9 - 10^{10} h^{-1} M_{odot}$. Interestingly, the fractal dimension is very different when the reionization proceeds inside-out compared to when it is outside-in. Thus the multifractal nature of HI density field at high redshifts can be used to study the nature of reionization.
Simulations estimating the differential brightness temperature of the redshifted 21-cm from the epoch of reionization (EoR) often assume that the spin temperature is decoupled from the background CMB temperature and is much larger than it. Although a valid assumption towards the latter stages of the reionization process, it does not necessarily hold at the earlier epochs. Violation of this assumption will lead to fluctuations in differential brightness temperature that are neither driven by density fluctuations nor by HII regions. Therefore, it is vital to calculate the spin temperature self-consistently by treating the Lyman-alpha and collisional coupling of spin temperature to the kinetic temperature. In this paper we develop an extension to the BEARS algorithm, originally developed to model reionization history, to include these coupling effects. Here we simulate the effect in ionization and heating for three models in which the reionization is driven by stars, miniqsos or a mixture of both.We also perform a number of statistical tests to quantify the imprint of the self-consistent inclusion of the spin temperature decoupling from the CMB. We find that the evolution of the spin temperature has an impact on the measured signal specially at redshifts higher than 10 and such evolution should be taken into account when one attempts to interpret the observational data.
The neutral hydrogen (HI) and its 21 cm line are promising probes to the reionization process of the intergalactic medium (IGM). To use this probe effectively, it is imperative to have a good understanding on how the neutral hydrogen traces the underlying matter distribution. Here we study this problem using semi-numerical modeling by combining the HI in the IGM and the HI from halos during the epoch of reionization (EoR), and investigate the evolution and the scale-dependence of the neutral fraction bias as well as the 21 cm line bias. We find that the neutral fraction bias on large scales is negative during reionization, and its absolute value on large scales increases during the early stage of reionization and then decreases during the late stage. During the late stage of reionization, there is a transition scale at which the HI bias transits from negative on large scales to positive on small scales, and this scale increases as the reionization proceeds to the end.
A major goal of observational and theoretical cosmology is to observe the largely unexplored time period in the history of our universe when the first galaxies form, and to interpret these measurements. Early galaxies dramatically impacted the gas around them in the surrounding intergalactic medium (IGM) by photoionzing the gas during the Epoch of Reionization (EoR). This epoch likely spanned an extended stretch in cosmic time: ionized regions formed and grew around early generations of galaxies, gradually filling a larger and larger fraction of the volume of the universe. At some time -- thus far uncertain, but within the first billion years or so after the big bang -- essentially the entire volume of the universe became filled with ionized gas. The properties of the IGM provide valuable information regarding the formation time and nature of early galaxy populations, and many approaches for studying the first luminous sources are hence based on measurements of the surrounding intergalactic gas. The prospects for improved reionization-era observations of the IGM and early galaxy populations over the next decade are outstanding. Motivated by this, we review the current state of models of the IGM during reionization. We focus on a few key aspects of reionization-era phenomenology and describe: the redshift evolution of the volume-averaged ionization fraction, the properties of the sources and sinks of ionizing photons, along with models describing the spatial variations in the ionization fraction, the ultraviolet radiation field, the temperature of the IGM, and the gas density distribution.