Extracting 21 cm signal by frequency and angular filtering

Extracting the neutral hydrogen(HI) signal is a great challenge for cosmological 21 cm experiments; both the astrophysical foregrounds and receiver noise are typically several orders of magnitude greater than the 21 cm signal. However, the different properties of the 21 cm signal, foreground and noise can be exploited to separate these components. The foregrounds are generally smooth or correlated over the frequency space along a line of sight(Lo S), while both the 21 cm signal and noise vary stochastically along the same Lo S. The foreground can be removed by filtering out the smooth component in frequency space. The receiver noise is basically uncorrelated for observations at different times, hence for surveys it is also uncorrelated in different directions, while the 21 cm signal, which traces the large scale structure, is correlated up to certain scales. In this exercise, we apply Wiener filters in frequency and angular space to extract the 21 cm signals. We found that the method works well. Inaccurate knowledge about the beam could degrade the reconstruction, but the overall result is still good, showing that the method is fairly robust.     

Directionality in the Wilkinson Microwave Anisotropy Probe polarization data

Polarization is the next frontier of cosmic microwave background analysis, but its signal is dominated over much of the sky by foregrounds which must be carefully removed. To determine the efficacy of this cleaning, it is necessary to have sensitive tests for residual foreground contamination in polarization sky maps. The dominant Galactic foregrounds introduce a large-scale anisotropy on to the sky, so it makes sense to use a statistic sensitive to overall directionality for this purpose. Here, we adapt the rapidly computable     statistic of Bunn and Scott to polarization data, and demonstrate its utility as a foreground monitor by applying it to the low-resolution Wilkinson Microwave Anisotropy Probe 3-yr sky maps. With a thorough simulation of the maps' noise properties, we find no evidence for contamination in the foreground cleaned sky maps.     

Reionization and the large-scale 21-cm cosmic microwave background cross-correlation

Of the many probes of reionization, the 21-cm line and the cosmic microwave background (CMB) are among the most effective. We examine how the cross-correlation of the 21-cm brightness and the CMB Doppler fluctuations on large angular scales can be used to study this epoch. We employ a new model of the growth of large-scale fluctuations of the ionized fraction as reionization proceeds. We take into account the peculiar velocity field of baryons and show that its effect on the cross-correlation can be interpreted as a mixing of Fourier modes. We find that the cross-correlation signal is strongly peaked towards the end of reionization and that the sign of the correlation should be positive because of the inhomogeneity inherent to reionization. The signal peaks at degree scales (ℓ∼ 100) and comes almost entirely from large physical scales ( k ∼ 10⁻² Mpc). Since many of the foregrounds and noise that plague low-frequency radio observations will not correlate with CMB measurements, the cross-correlation might appear to provide a robust diagnostic of the cosmological origin of the 21-cm radiation around the epoch of reionization. Unfortunately, we show that these signals are actually only weakly correlated and that cosmic variance dominates the error budget of any attempted detection. We conclude that the detection of a cross-correlation peak at degree-size angular scales is unlikely even with ideal experiments.     

Foregrounds for 21-cm observations of neutral gas at high redshift

We investigate a number of potential foregrounds for an ambitious goal of future radio telescopes such as the Square Kilometer Array (SKA) and the Low Frequency Array (LOFAR): spatial tomography of neutral gas at high redshift in 21-cm emission. While the expected temperature fluctuations due to unresolved radio point sources is highly uncertain, we point out that free–free emission from the ionizing haloes that reionized the Universe should define a minimal bound. This emission is likely to swamp the expected brightness temperature fluctuations, making proposed detections of the angular patchwork of 21-cm emission across the sky unlikely to be viable. Hα observations with JWST could place an upper bound on the contribution of high-redshift sources to the free–free background. An alternative approach is to discern the topology of reionization from spectral features due to 21-cm emission along a pencil-beam slice. This requires tight control of the frequency-dependence of the beam in order to prevent foreground sources from contributing excessive variance. We also investigate potential contamination by galactic and extragalactic radio recombination lines (RRLs). These are unlikely to be show-stoppers, although little is known about the distribution of RRLs away from the Galactic plane. The mini-halo emission signal is always less than that of the intergalactic medium (IGM), making mini-haloes unlikely to be detectable. If they are seen, it will be only in the very earliest stages of structure formation at high redshift, when the spin temperature of the IGM has not yet decoupled from the cosmic microwave background.     

Will point sources spoil 21-cm tomography?

21-cm tomography is emerging as a promising probe of the cosmological dark ages and the epoch of reionization, as well as a tool for observational cosmology in general. However, serious sources of foreground contamination must be subtracted for experimental efforts to be viable. In this paper, we focus on the removal of unresolved extragalactic point sources with smooth spectra, and evaluate how the residual foreground contamination after cleaning depends on instrumental and algorithmic parameters. A crucial but often ignored complication is that the synthesized beam of an interferometer array shrinks towards higher frequency, causing complicated frequency structure in each sky pixel as 'frizz' far from the beam centre contracts across unresolved radio sources. We find that current-generation experiments should none the less be able to clean out this point source contamination adequately, and quantify the instrumental and algorithmic design specifications required to meet this foreground challenge.     

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.     

Large-scale non-Gaussianities in the 21-cm background anisotropies from the era of reionization

The brightness temperature fluctuations in the 21-cm background related to the neutral hydrogen distribution provide a probe of the physics related to the era of reionization, when the intergalactic medium changed from being completely neutral to partially ionized. We formulate statistics of 21-cm brightness temperature anisotropies in terms of the angular power spectrum, the bispectrum, and the trispectrum. Using the trispectrum, we estimate the covariance related to the power spectrum measurements and show that correlations resulting from non-Gaussianities are below a per cent, at most. While all-sky observations of the 21-cm background at arcminute-scale resolution can be used to measure the bispectrum with a cumulative signal-to-noise ratio of the order of a few tens, in the presence of foregrounds and instrumental noise related to first-generation interferometers, the measurement is unlikely to be feasible. For most purposes, non-Gaussianities in 21-cm fluctuations can be ignored and the distribution can be described with Gaussian statistics. Because 21-cm fluctuations are significantly contaminated by foregrounds, such as galactic synchrotron or low-frequency radio point sources, the lack of significant non-Gaussianity in the signal suggests that any significant detection of non-Gaussianity could be the result of foregrounds. Similarly, in addition to the frequency information that is now proposed to separate 21-cm fluctuations from foregrounds, if the non-Gaussian structure of foregrounds is known a priori, this additional information could potentially be used to reduce the confusion further.     

Detection and extraction of signals from the epoch of reionization using higher-order one-point statistics

Detecting redshifted 21-cm emission from neutral hydrogen in the early Universe promises to give direct constraints on the epoch of reionization (EoR). It will, though, be very challenging to extract the cosmological signal (CS) from foregrounds and noise which are orders of magnitude larger. Fortunately, the signal has some characteristics which differentiate it from the foregrounds and noise, and we suggest that using the correct statistics may tease out signatures of reionization. We generate mock data cubes simulating the output of the Low Frequency Array (LOFAR) EoR experiment. These cubes combine realistic models for Galactic and extragalactic foregrounds and the noise with three different simulations of the CS. We fit out the foregrounds, which are smooth in the frequency direction, to produce residual images in each frequency band. We denoise these images and study the skewness of the one-point distribution in the images as a function of frequency. We find that, under sufficiently optimistic assumptions, we can recover the main features of the redshift evolution of the skewness in the 21-cm signal. We argue that some of these features – such as a dip at the onset of reionization, followed by a rise towards its later stages – may be generic, and give us a promising route to a statistical detection of reionization.     

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.     

基于红移畸变测量宇宙结构增长率的进展

星系红移巡天的一个主要目标是依据光谱红移测距,详细刻画宇宙中星系的三维空间分布。由于星系本动速度的存在,红移空间的星系分布存在着严重畸变,在大小尺度上有着不同模式的各向同性偏离。通过对红移畸变的观测研究,人们可从中获取速度场的信息,因此,红移畸变已成为暗能量探测的重要探针之一,为检验宇宙学尺度上的引力模型提供帮助。当前星系红移巡天项目已经取得了非凡成功,为人们提供了详细的星系空间分布数据。人们据此测量了星系的相关函数和功率谱,提取了精确的红移畸变信号,并通过模型拟合限制出了一批不同红移处宇宙结构增长率的估值,为探索宇宙尺度的引力模式提供了数据支持。主要介绍红移畸变模型、星系红移巡天观测和宇宙结构增长率测量等研究进展。     

High-resolution imaging of the cosmic mass distribution from gravitational lensing of pre-galactic H i

Low-frequency radio observations of neutral hydrogen during and before the epoch of cosmic re-ionization will provide ∼1000 quasi-independent source planes, each of precisely known redshift, if a resolution of ∼1 arcmin or better can be attained. These planes can be used to reconstruct the projected mass distribution of foreground material. Structure in these source planes is linear and Gaussian at high redshift  (30 < z < 300)  but is non-linear and non-Gaussian during re-ionization. At both epochs, significant power is expected down to subarcsecond scales. We demonstrate that this structure can, in principle, be used to make mass images with a formal signal-to-noise ratio (S/N) per pixel exceeding 10, even for pixels as small as an arcsecond. With an ideal telescope, both resolution and S/N can exceed those of even the most optimistic idealized mass maps from galaxy lensing by more than an order of magnitude. Individual dark haloes similar in mass to that of the Milky Way could be imaged with high S/N out to   z ∼ 10  . Even with a much less ambitious telescope, a wide-area survey of 21-cm lensing would provide very sensitive constraints on cosmological parameters, in particular on dark energy. These are up to 20 times tighter than the constraints obtainable from comparably sized, very deep surveys of galaxy lensing, although the best constraints come from combining data of the two types. Any radio telescope capable of mapping the 21-cm brightness temperature with good frequency resolution (∼0.05 MHz) over a band of width ≳10 MHz should be able to make mass maps of high quality. The planned Square Kilometre Array may be able to map the mass with moderate S/N down to arcminute scales, depending on the re-ionization history of the universe and the ability to subtract foreground sources.     

Prospects for Detecting the 326.5 MHz Redshifted 21-cm HI Signal with the Ooty Radio Telescope (ORT)

Observations of the redshifted 21-cm HI fluctuations promise to be an important probe of the post-reionization era (z≤ 6). In this paper we calculate the expected signal and foregrounds for the upgraded Ooty Radio Telescope (ORT) which operates at frequency ν _o = 326.5 MHz which corresponds to redshift z = 3.35. Assuming that the visibilities contain only the HI signal and system noise, we show that a 3 σ detection of the HI signal (～1 mK) is possible at angular scales 11^′ to 3^° with ≈1000 h of observation. Foreground removal is one of the major challenges for a statistical detection of the redshifted 21 cm HI signal. We assess the contribution of different foregrounds and find that the 326.5 MHz sky is dominated by the extragalactic point sources at the angular scales of our interest. The expected total foregrounds are 10⁴?10⁵ times higher than the HI signal.     

The 21-cm power spectrum after reionization

We discuss the 21-cm power spectrum (PS) following the completion of reionization. In contrast to the reionization era, this PS is proportional to the PS of mass density fluctuations, with only a small modulation due to fluctuations in the ionization field on scales larger than the mean-free-path of ionizing photons. We derive the form of this modulation, and demonstrate that its effect on the 21-cm PS will be smaller than 1 per cent for physically plausible models of damped Lyα systems. In contrast to the 21-cm PS observed prior to reionization, in which H  ii regions dominate the ionization structure, the simplicity of the 21-cm PS after reionization will enhance its utility as a cosmological probe by removing the need to separate the PS into physical and astrophysical components. As a demonstration, we consider the Alcock–Paczynski test and show that the next generation of low-frequency arrays could measure the angular distortion of the PS at the per cent level for   z ∼ 3–5  .     

Foreground removal of 21 cm fluctuation with multifrequency fitting

The 21 centimeter (21 cm) line emission from neutral hydrogen in the inter-galactic medium (IGM) at high redshifts is strongly contaminated by foreground sources such as the diffuse Galactic synchrotron emission and free-free emission from the Galaxy, as well as emission from extragalactic radio sources, thus making its observation very complicated. However, the 21 cm signal can be recovered through its structure in fre-quency space, as the power spectrum of the foreground contamination is expected to be smooth over a wide band in frequency space while the 21 cm fluctuations vary signifi-cantly. We use a simple polynomial fitting to reconstruct the 21 cm signal around four frequencies 50, 100, 150 and 200 MHz with an especially small channel width of 20 kHz. Our calculations show that this multifrequency fitting approach can effectively recover the 21 cm signal in the frequency range 100 ～ 200 MHz. However, this method doesn't work well around 50 MHz because of the low intensity of the 21 cm signal at this frequency. We also show that the fluctuation of detector noise can be suppressed to a very low level by taking long integration times, which means that we can reach a sensitivity of ≈ 10 mK at 150 MHz with 40 antennas in 120 hours of observations.     

Cosmic bandits: Exploration versus exploitation in CMB B-mode experiments

A preferred method to detect the curl-component, or B-mode, signature of inflationary gravitational waves (IGWs) in the cosmic microwave background (CMB) polarization, in the absence of foregrounds and lensing, is a prolonged integration over a single patch of sky of a few square degrees. In practice, however, foregrounds abound and the sensitivity to B modes can be improved considerably by finding the region of sky cleanest of foregrounds. The best strategy to detect B modes thus involves a tradeoff between exploration (to find lower-foreground patches) and exploitation (through prolonged integration). This problem is akin to the multi-armed bandit (MAB) problem in probability theory, wherein a gambler faces a series of slot machines with unknown winning odds and must develop a strategy to maximize his/her winnings with some finite number of pulls. While the optimal MAB strategy remains to be determined, a number of algorithms have been developed in an effort to maximize the winnings. Here, based on this resemblance, we tackle the search for IGW B modes with single frequency experiments in the presence of spatially varying foregrounds by developing adaptive survey strategies to optimize the sensitivity to IGW B modes. We demonstrate, using realistic foreground models and taking lensing-induced B modes into account, that adaptive experiments can substantially improve the upper bound on the tensor-to-scalar ratio (by factors of 2 and 3 in single frequency experiments, and possibly even more). Similar techniques can be applied to other surveys, including 21-cm measurements of signatures of the epoch of reionization, searches for a stochastic primordial gravitational wave background, deep-field imaging by the James Webb Space Telescope or various radio interferometers, and transient follow-up searches.     

H i 21-cm absorption at z∼ 3.39 towards PKS 0201+113

We report the Giant Metrewave Radio Telescope detection of H  i 21-cm absorption from the z ∼ 3.39 damped Lyman α absorber (DLA) towards PKS 0201+113, the highest redshift at which 21-cm absorption has been detected in a DLA. The absorption is spread over ∼115 km s⁻¹ and has two components, at   z = 3.387 144(17)  and   z = 3.386 141  (45). The stronger component has a redshift and velocity width in agreement with the tentative detection of Briggs, Brinks & Wolfe, but a significantly lower optical depth. The core size and DLA covering factor are estimated to be ≲100 pc and f ∼ 0.69, respectively, from a Very Long Baseline Array 328-MHz image. If one makes the conventional assumption that the H  i column densities towards the optical and radio cores are the same, this optical depth corresponds to a spin temperature of T _s∼[(955 ± 160) × ( f /0.69)] K. However, this assumption may not be correct, given that no metal-line absorption is seen at the redshift of the stronger 21-cm component, indicating that this component does not arise along the line of sight to the optical quasi-stellar object (QSO), and that there is structure in the 21-cm absorbing gas on scales smaller than the size of the radio core. We model the 21-cm absorbing gas as having a two-phase structure with cold dense gas randomly distributed within a diffuse envelope of warm gas. For such a model, our radio data indicate that, even if the optical QSO lies along a line of sight with a fortuitously high (∼50 per cent) cold gas fraction, the average cold gas fraction is low, ≲17 per cent, when averaged over the spatial extent of the radio core. Finally, the large mismatch between peak 21-cm and optical redshifts and the complexity of both profiles makes it unlikely that the z ∼ 3.39 DLA will be useful in tests of fundamental constant evolution.     

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.     

21-cm fluctuations from inhomogeneous X-ray heating before reionization

Many models of early structure formation predict a period of heating immediately preceding reionization, when X-rays raise the gas temperature above that of the cosmic microwave background. These X-rays are often assumed to heat the intergalactic medium (IGM) uniformly, but in reality will heat the gas more strongly closer to the sources. We develop a framework for calculating fluctuations in the 21-cm brightness temperature that originate from this spatial variation in the heating rate. High-redshift sources are highly clustered, leading to significant gas temperature fluctuations (with fractional variations ∼40 per cent, peaking on   k ∼ 0.1 Mpc⁻¹  scales). This induces a distinctive peak-trough structure in the angle-averaged 21-cm power spectrum, which may be accessible to the proposed Square Kilometre Array. This signal reaches the ∼10 mK level, and is stronger than that induced by Lyα flux fluctuations. As well as probing the thermal evolution of the IGM before reionization, this 21-cm signal contains information about the spectra of the first X-ray sources. Finally, we consider disentangling temperature, density and Lyα flux fluctuations as functions of redshift.     

Gravitational wave background from neutrino-driven gamma-ray bursts

We study cosmic microwave background (CMB) secondary anisotropies produced by inhomogeneous reionization by means of cosmological simulations coupled with the radiative transfer code crash . The reionization history is consistent with the Wilkinson Microwave Anisotropy Probe Thomson optical depth determination. We find that the signal arising from this process dominates over the primary CMB component for   l ≳ 4000  and reaches a maximum amplitude of   l ( l + 1) C_l /2π≃ 1.6 × 10⁻¹³  on arcmin scales (i.e. l as large as several thousands). We then cross-correlate secondary CMB anisotropy maps with neutral hydrogen 21-cm line emission fluctuations obtained from the same simulations. The two signals are highly anticorrelated on angular scales corresponding to the typical size of H  ii regions (including overlapping) at the 21-cm map redshift. We show how the CMB/21-cm cross-correlation can be used: (i) to study the nature of the reionization sources; (ii) to reconstruct the cosmic reionization history; (iii) to infer the mean cosmic ionization level at any redshift. We discuss the feasibility of the proposed experiment with forthcoming facilities.