The difference PDF of 21-cm fluctuations: a powerful statistical tool for probing cosmic reionization

A new generation of radio telescopes are currently being built with the goal of tracing the cosmic distribution of atomic hydrogen at redshifts 6–15 through its 21-cm line. The observations will probe the large-scale brightness fluctuations sourced by ionization fluctuations during cosmic reionization. Since detailed maps will be difficult to extract due to noise and foreground emission, efforts have focused on a statistical detection of the 21-cm fluctuations. During cosmic reionization, these fluctuations are highly non-Gaussian and thus more information can be extracted than just the one-dimensional function that is usually considered, i.e. the correlation function. We calculate a two-dimensional function that if measured observationally would allow a more thorough investigation of the properties of the underlying ionizing sources. This function is the probability distribution function (PDF) of the difference in the 21-cm brightness temperature between two points, as a function of the separation between the points. While the standard correlation function is determined by a complicated mixture of contributions from density and ionization fluctuations, we show that the difference PDF holds the key to separately measuring the statistical properties of the ionized regions.     

Biased reionization and non-Gaussianity in redshifted 21-cm intensity maps of the reionization epoch

Spatial dependence in the statistics of redshifted 21-cm fluctuations promises to provide the most powerful probe of the reionization epoch. In this paper we consider the second and third moments of the redshifted 21-cm intensity distribution using a simple model that accounts for galaxy bias during the reionization process. We demonstrate that skewness in redshifted 21-cm maps should be substantial throughout the reionization epoch and on all angular scales, owing to the effects of galaxy bias which leads to early reionization in overdense regions of the intergalactic medium (IGM). The variance (or power spectrum) of 21-cm fluctuations will exhibit a minimum in redshift part way through the reionization process, when the global ionization fraction is around 50 per cent. This minimum is generic, and is due to the transition from 21-cm intensity being dominated by overdense to underdense regions as reionization progresses. We show that the details of the reionization history, including the presence of radiative feedback are encoded in the evolution of the autocorrelation and skewness functions with redshift and mean IGM neutral fraction. The amplitudes of fluctuations are particularly sensitive to the masses of ionizing sources, and vary by an order of magnitude for astrophysically plausible models. We discuss the detection of skewness by first-generation instruments, and conclude that the Mileura Wide-field Array–Low-Frequency Demonstrator will have sufficient sensitivity to detect skewness on a range of angular scales at redshifts near the end of reionization, while a subsequent instrument of 10 times the collecting area could map out the evolution of skewness in detail. The observation of a minimum in variance during the reionization history, and the detection of skewness would both provide important confirmation of the cosmological origin of redshifted 21-cm intensity fluctuations.     

External calibrator in global signal experiment for detection of the epoch of reionization

We present a conceptual design study of external calibrators in the 21 cm experiment towards detecting the globally averaged radiation of the epoch of reionization(EoR). Employment of external calibrator instead of internal calibrator commonly used in current EoR experiments allows removing instrumental effects such as beam pattern, receiver gain and instability of the system if the conventional three-position switch measurements are implemented in a short time interval. Furthermore, in the new design the antenna system is placed in an underground anechoic chamber with an open/closing ceiling to maximally reduce the environmental effect such as RFI and ground radiation/reflection. It appears that three of the four external calibrators proposed in this paper, including two indoor artificial transmitters and one outdoor celestial radiation(the Galactic polarization), fail to meet our purpose. Diurnal motion of the Galactic diffuse emission turns out to be the most probable source as an external calibrator, for which we have discussed the observational strategy and the algorithm of extracting the EoR signal.     

Studying the sources of cosmic reionization with 21-cm fluctuations

We explore the ability of measurements of the 21-cm power spectrum during reionization to enable the simultaneous reconstruction of the reionization history and the properties of the ionizing sources. For various sets of simulated 21-cm observations, we perform maximum likelihood fits in order to constrain the reionization and galaxy formation histories. We employ a flexible six-parameter model that parametrizes the uncertainties in the properties of high-redshift galaxies. The computational speed needed is attained through the use of an analytical model that is in reasonable agreement with numerical simulations of reionization. We find that one-year observations, with the Murchison Widefield Array, should measure the cosmic ionized fraction to  ∼1 per cent  accuracy at the very end of reionization, and a few per cent accuracy around the mid-point of reionization. The mean halo mass of the ionizing sources should be measurable to 10 per cent accuracy when reionization is 2/3 of the way through, and to 20 per cent accuracy throughout the central stage of reionization, if this mass is anywhere in the range 1/3 to 100 billion solar masses.     

The correlation between star formation and 21-cm emission during the reionization epoch

Reionization is thought to be dominated by low-mass galaxies, while direct observations of resolved galaxies probe only the most massive, rarest objects. The cross-correlation between fluctuations in the surface brightness of the cumulative Lyα emission (which serves as a proxy for the star formation rate) and the redshifted 21-cm signal from neutral hydrogen in the intergalactic medium (IGM) will directly probe the causal link between the production of ionizing photons in galaxies and the reionization of the IGM. We discuss the prospects for detecting this cross-correlation for unresolved galaxies. We find that on angular scales ≲10 arcmin detection will be practical using wide-field near-infrared (near-IR) imaging from space in combination with the forthcoming Mileura Wide-field Array – Low Frequency Demonstrator. When redshifted 21-cm observations of the neutral IGM are combined with space-based near-IR imaging of Lyα emission, the detection on angular scales ≲3 arcmin will be limited by the sensitivity of the 21-cm signal, even when a small-aperture optical telescope (∼2 m) and a moderate field of view (∼10 deg²) are used. On scales ≳3 arcmin, the measurement of cross-correlation will be limited by the accuracy of the foreground sky subtraction.     

Modification of the 21-cm power spectrum by quasars during the epoch of reionization

We assess the effect of a population of high-redshift quasars on the 21-cm power spectrum during the epoch of reionization. Our approach is to implement a seminumerical scheme to calculate the three-dimensional structure of ionized regions surrounding massive haloes at high redshift. We include the ionizing influence of luminous quasars by populating a simulated overdensity field with quasars using a Monte Carlo Markov Chain algorithm. We find that quasars modify both the amplitude and shape of the power spectrum at a level which is of the same order as the fractional contribution to reionization. The modification is found both at constant redshift and at constant global neutral fraction, and arises because ionizing photons produced by quasars are biased relative to the density field at a level that is higher than stellar ionizing photons. The modification of the power spectrum is likely to be small, rendering the effect of quasars difficult to isolate. However, we find the modification of the power spectrum by quasars to be at a level that is comparable to the precision expected for future low-frequency telescopes. Correct interpretation of observations will therefore require the effect of quasars to be considered, and our results imply that quasar ionization will need to be included in detailed modelling of observed 21-cm power spectra.     

The cross-correlation of the CMB polarization and the 21-cm line fluctuations from cosmic reionization

The cosmic microwave background (CMB) polarization and the 21-cm line fluctuations are powerful probes of cosmological reionization. We study how the cross-correlation between the CMB polarization ( E modes) and the 21-cm line fluctuations can be used to gain further understanding of the reionization history, within the framework of inhomogeneous reionization. Since the E -mode polarization reflects the amplitude of the quadrupole component of the CMB temperature fluctuations, the angular power spectrum of the cross-correlation exhibits oscillations at all multipoles. The first peak of the power spectrum appears at the scale corresponding to the quadrupole at the redshift, which is probed by the 21-cm line fluctuations. The peak reaches its maximum value in redshift when the average ionization fraction of the universe is about half. On the other hand, on small scales, there is a damping that depends on the duration of reionization. Thus, the cross-correlation between the CMB polarization and the 21-cm line fluctuations has the potential to accurately constrain the epoch and the duration of reionization.     

The contribution of the IGM and minihaloes to the 21-cm signal of reionization

We study the statistical properties of the cosmological 21-cm signal from both the intergalactic medium (IGM) and minihaloes, using a reionization simulation that includes a self-consistent treatment of minihalo photoevaporation. We consider two models for minihalo formation and three typical thermal states of the IGM – heating purely by ionization, heating from both ionizing and Lyα photons and a maximal 'strong heating' model. We find that the signal from the IGM is almost always dominant over that from minihaloes. In our calculation, the differential brightness temperature,  δ T _b,  of minihaloes is never larger than 2 mK. Although there are indeed some differences in the signals from the minihaloes and from the IGM, even with the planned generation of radio telescopes it will be unfeasible to detect them. However, minihaloes significantly affect the ionization state of the IGM and the corresponding 21-cm flux.     

The impact of H i in galaxies on 21-cm intensity fluctuations during the reionization epoch

We investigate the impact of neutral hydrogen (H  i ) in galaxies on the statistics of 21-cm fluctuations using seminumerical modelling. Following the reionization of hydrogen, the H  i content of the Universe is dominated by damped absorption systems (DLAs), with a cosmic density in H  i that is observed to be constant at a level equal to ∼2 per cent of the cosmic baryon density from   z ∼ 1  to   z ∼ 5  . We show that extrapolation of this constant fraction into the reionization epoch results in a reduction in the amplitude of 21-cm fluctuations over a range of spatial scales. We further find that consideration of H  i in galaxies/DLAs reduces the prominence of the H  ii region induced shoulder in the 21-cm power spectrum (PS), and hence modifies the scale dependence of 21-cm fluctuations. We also estimate the 21-cm–galaxy cross PS and show that the cross PS changes sign on scales corresponding to the H  ii regions. From consideration of the sensitivity for forthcoming low-frequency arrays, we find that the effects of H  i in galaxies/DLAs on the statistics of 21-cm fluctuations will be significant with respect to the precision of a PS or cross PS measurement. In addition, since overdense regions are reionized first we demonstrate that the cross-correlation between galaxies and 21-cm emission changes sign at the end of the reionization era, providing an alternative avenue to pinpoint the end of reionization. The sum of our analysis indicates that the H  i content of the galaxies that reionize the universe will need to be considered in detailed modelling of the 21-cm intensity PS in order to correctly interpret measurements from forthcoming low-frequency arrays.     

Cross-correlation of 21 cm and soft X-ray backgrounds during the epoch of reionization

The cross-correlation between the high-redshift 21 cm background and the Soft X-ray Background(SXB) of the Universe may provide an additional probe of the Epoch of Reionization. Here we use semi-numerical simulations to create 21 cm and soft X-ray intensity maps and construct their cross power spectra. Our results indicate that the cross power spectra are sensitive to the thermal and ionizing states of the intergalactic medium(IGM). The 21 cm background correlates positively to the SXB on large scales during the early stages of the reionization. However as the reionization develops, these two backgrounds turn out to be anti-correlated with each other when more than ~15% of the IGM is ionized in a warm reionization scenario. The anti-correlated power reaches its maximum when the neutral fraction declines to 0.2–0.5. Hence, the trough in the cross power spectrum might be a useful tool for tracing the growth of HII regions during the middle and late stages of the reionization. We estimate the detectability of the cross power spectrum based on the abilities of the Square Kilometre Array and the Wide Field X-ray Telescope(WFXT), and find that to detect the cross power spectrum, the pixel noise of X-ray images has to be at least 4 orders of magnitude lower than that of the WFXT deep survey.     

Foregrounds for redshifted 21-cm studies of reionization: Giant Meter Wave Radio Telescope 153-MHz observations

Foreground subtraction is the biggest challenge for future redshifted 21-cm observations to probe reionization. We use a short Giant Meter Wave Radio Telescope (GMRT) observation at 153 MHz to characterize the statistical properties of the background radiation across ∼1° to subarcmin angular scales, and across a frequency band of 5 MHz with 62.5 kHz resolution. The statistic we use is the visibility correlation function, or equivalently the angular power spectrum   C _l   . We present the results obtained from using relatively unsophisticated, conventional data calibration procedures. We find that even fairly simple-minded calibration allows one to estimate the visibility correlation function at a given frequency   V ₂( U , 0)  . From our observations, we find that   V ₂( U , 0)  is consistent with foreground model predictions at all angular scales except the largest ones probed by our observations where the model predictions are somewhat in excess. On the other hand, the visibility correlation between different frequencies  κ( U , Δν)  seems to be much more sensitive to calibration errors. We find a rapid decline in  κ( U , Δν)  , in contrast with the prediction of less than 1 per cent variation across 2.5 MHz. In this case, however, it seems likely that a substantial part of the discrepancy may be due to limitations of data reduction procedures.     

Spatial fluctuations in the spectral shape of the ultraviolet background at 2 < z < 3 and the reionization of helium

The low-density hydrogen and helium in the intergalactic medium (IGM) probed by quasi-stellar object (QSO) absorption lines is sensitive to the amplitude and spectral shape of the metagalactic ultraviolet (UV) background. We use realistic H  i and He  ii Lyα forest spectra, constructed from state-of-the-art hydrodynamical simulations of a Λ cold dark matter (ΛCDM) universe to confirm the reliability of using line profile fitting techniques to infer the ratio of the metagalactic H  i and He  ii ionization rates. We further show that the large spatial variations and the anticorrelation with H  i absorber density observed in the ratio of the measured He  ii to H  i column densities can be explained in a model where the H  i ionization rate is dominated by the combined UV emission from young star-forming galaxies and QSOs and the He  ii ionization rate is dominated by emission from QSOs only. In such a model the large fluctuations in the column density ratio are due to the small number of QSOs expected to contribute at any given point to the He  ii ionization rate. A significant contribution to UV emission at the He  ii photoelectric edge from hot gas in galaxies and galaxy groups would decrease the expected fluctuations in the column density ratio. Consequently, this model appears difficult to reconcile with the large increase in He  ii opacity fluctuations towards higher redshift. Our results further strengthen previous suggestions that observed He  ii Lyα forest spectra at z ∼ 2–3.5 probe the tail end of the reionization of He  ii by QSOs.     

Seismology of stellar envelopes: probing the outer layers of a star through the scattering of acoustic waves

The outer layers of Sun-like stars are regions of rapid spatial variation which modulate the p-mode frequencies by partially reflecting the constituent acoustic waves. With the accuracy that has been achieved by current solar observations, and that is expected from imminent stellar observations, this modulation can be observed from the spectra of the low-degree modes. We present a new and simple theoretical calculation to determine the leading terms in an asymptotic expansion of the outer phase of these modes, which is determined by the structure of the surface layers of the star. Our procedure is to compare the stellar envelope with a plane-parallel polytropic envelope, which we regard as a smooth reference background state. Then we can isolate a seismic signature of the acoustic phase and relate it to the stratification of the outer layers of the convection zone. One can thereby constrain theories of convection that are used to construct the convection zones of the Sun and Sun-like stars. The accuracy of the diagnostic is tested in the solar case by comparing the predicted outer phase with an exact numerical calculation.     

An automated station for the detection of bursts of cosmic origin at VHF and UHF

An automated system for the detection of bursts of cosmic origin at radio frequencies has been built at Medicina (Bologna). The station, completely independent, is part of a network of observatories based in Ireland, United Kingdom, Switzerland and Italy. The receiving apparatus, operating at 151 and 408 MHz, and the experimental method used to reject unwanted events is described. Observational programmes include a search for the radio counterpart of X-and -ray bursts, observations of the galactic centre and Cygnus regions, and a full sky monitoring. Sensitivity limits as well as preliminary results are discussed.