A method to measure the mass of damped Lyα absorber host galaxies using fluctuations in 21-cm emission

Observations of damped Lyα absorbers (DLAs) indicate that the fraction of hydrogen in its neutral form (H  i ) is significant by mass at all redshifts. This gas represents the reservoir of material that is available for star formation at late times. As a result, observational identification of the systems in which this neutral hydrogen resides is an important missing ingredient in models of galaxy formation. Precise identification of DLA host mass via traditional clustering studies is not practical owing to the small numbers of known systems being spread across sparsely distributed sightlines. However, following the completion of re-ionization, 21-cm surface brightness fluctuations will be dominated by neutral hydrogen in DLAs. No individual DLAs could be detected in 21-cm emission. Rather, observations of these fluctuations will measure the combined clustering signal from all DLAs within a large volume. We show that measurement of the spherically averaged power spectrum of 21-cm intensity fluctuations due to DLAs could be used to measure the galaxy bias for DLA host galaxies when combined with an independent measurement of the cosmological H  i mass density from quasar absorption studies. Utilizing this technique, the low-frequency arrays now under construction could measure the characteristic DLA host mass with a statistical precision as low as 0.3 dex at z ≳ 4. In addition, high signal-to-noise ratio observations of the peculiar-motion-induced anisotropy of the power spectrum would facilitate measurement of both the DLA host mass and the cosmic H  i density directly from 21-cm fluctuations. By exploiting this anisotropy, a second generation of low-frequency arrays with an order of magnitude increase in collecting area could measure the values of cosmic H  i density and DLA host mass, with uncertainties of a few per cent and a few tens of per cent, respectively.     

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  .     

Redshifted 21-cm observations of high-redshift quasar proximity zones

The introduction of low-frequency radio arrays is expected to revolutionize the study of the reionization epoch. Observation of the contrast in redshifted 21-cm emission between a large H  ii region and the surrounding neutral intergalactic medium (IGM) will be the simplest and most easily interpreted signature. However, the highest redshift quasars known are thought to reside in an ionized IGM. Using a semi-analytic model we describe the redshifted 21-cm signal from the IGM surrounding quasars discovered using the i -drop-out technique (i.e. quasars at   z ∼ 6  ). We argue that while quasars at   z < 6.5  seem to reside in the post-overlap IGM, they will still provide valuable probes of the late stages of the overlap era because the light-travel time across a quasar proximity zone should be comparable to the duration of overlap. For redshifted 21-cm observations within a 32-MHz bandpass, we find that the subtraction of a spectrally smooth foreground will not remove spectral features due to the proximity zone. These features could be used to measure the neutral hydrogen content of the IGM during the late stages of reionization. The density of quasars at   z ∼ 6  is now well constrained. We use the measured quasar luminosity function to estimate the prospects for discovery of high-redshift quasars in fields that will be observed by the Murchison Widefield Array.     

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.     

Spin temperatures and covering factors for H i 21-cm absorption in damped Lyman α systems

We investigate the practice of assigning high spin temperatures to damped Lyman α absorption systems (DLAs) not detected in H  i 21-cm absorption. In particular, Kanekar & Chengalur have attributed the mix of 21-cm detections and non-detections in low-redshift  ( z _abs≤ 2.04) DLAs  to a mix of spin temperatures, while the non-detections at high redshift were attributed to high spin temperatures. Below   z _abs= 0.9  , where some of the DLA host galaxy morphologies are known, we find that 21-cm absorption is normally detected towards large radio sources when the absorber is known to be associated with a large intermediate (spiral) galaxy. Furthermore, at these redshifts, only one of the six 21-cm non-detections has an optical identification and these DLAs tend to lie along the sight-lines to the largest background radio continuum sources. For these and many of the high-redshift DLAs occulting large radio continua, we therefore expect covering factors of less than the assumed/estimated value of unity. This would have the effect of introducing a range of spin temperatures considerably narrower than the current range of  Δ T _s≳ 9000 K  , while still supporting the hypothesis that the high-redshift DLA sample comprises a larger proportion of compact galaxies than the low-redshift sample.     

Baryonic acoustic oscillations in 21-cm emission: a probe of dark energy out to high redshifts

Low-frequency observatories are currently being constructed with the goal of detecting redshifted 21-cm emission from the epoch of reionization. These observatories will also be able to detect intensity fluctuations in the cumulative 21-cm emission after reionization, from hydrogen in unresolved damped Lyα absorbers (such as gas-rich galaxies) down to a redshift z ∼ 3.5. The inferred power spectrum of 21-cm fluctuations at all redshifts will show acoustic oscillations, whose comoving scale can be used as a standard ruler to infer the evolution of the equation of state for the dark energy. We find that the first generation of low-frequency experiments (such as MWA or LOFAR) will be able to constrain the acoustic scale to within a few per cent in a redshift window just prior to the end of the reionization era, provided that foregrounds can be removed over frequency bandpasses of ≳8 MHz. This sensitivity to the acoustic scale is comparable to the best current measurements from galaxy redshift surveys, but at much higher redshifts. Future extensions of the first-generation experiments (involving an order of magnitude increase in the antennae number of the MWA) could reach sensitivities below 1 per cent in several redshift windows and could be used to study the dark energy in the unexplored redshift regime of 3.5 ≲ z ≲ 12. Moreover, new experiments with antennae designed to operate at higher frequencies would allow precision measurements (≲1 per cent) of the acoustic peak to be made at more moderate redshifts (1.5 ≲ z ≲ 3.5), where they would be competitive with ambitious spectroscopic galaxy surveys covering more than 1000 deg². Together with other data sets, observations of 21-cm fluctuations will allow full coverage of the acoustic scale from the present time out to z ∼ 12.     

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.     

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.     

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.     

Current models of the observable consequences of cosmic reionization and their detectability

A number of large current experiments aim to detect the signatures of the cosmic reionization at redshifts z > 6. Their success depends crucially on understanding the character of the reionization process and its observable consequences and designing the best strategies to use. We use large-scale simulations of cosmic reionization to evaluate the reionization signatures at redshifted 21-cm and small-scale cosmic microwave background (CMB) anisotropies in the best current model for the background universe, with fundamental cosmological parameters given by Wilkinson Microwave Anisotropy Probe three-year results. We find that the optimal frequency range for observing the 'global step' of the 21-cm emission is 120–150 MHz, while statistical studies should aim at 140–160 MHz, observable by GMRT. Some strongly non-Gaussian brightness features should be detectable at frequencies up to ∼190 MHz. In terms of sensitivity-signal trade-off relatively low resolutions, corresponding to beams of at least a few arcminutes, are preferable. The CMB anisotropy signal from the kinetic Sunyaev–Zel'dovich effect from reionized patches peaks at tens of μK at arcminute scales and has an rms of ∼1 μK, and should be observable by the Atacama Cosmology Telescope and the South Pole Telescope. We discuss the various observational issues and the uncertainties involved, mostly related to the poorly known reionization parameters and, to a lesser extend, to the uncertainties in the background cosmology.     

Galaxy redshift surveys selected by neutral hydrogen using the Five-hundred metre Aperture Spherical Telescope

We discuss the possibility of performing a substantial spectroscopic galaxy redshift survey selected via the 21-cm emission from neutral hydrogen using the Five-hundred metre Aperture Spherical Telescope (FAST) to be built in China. We consider issues related to the estimation of the source counts and optimizations of the survey, and discuss the constraints on cosmological models that such a survey could provide. We find that a survey taking around two years could detect ∼10⁷ galaxies with an average redshift of ∼0.15 making the survey complementary to those already carried out at optical wavelengths. These conservative estimates have used the   z = 0  H  i mass function and have ignored the possibility of evolution. The results could be used to constrain  Γ=Ω_m h   to 5 per cent and the spectral index, n _s, to 7 per cent independent of cosmic microwave background data. If we also use simulated power spectra from the Planck satellite, we can constrain w to be within 5 per cent of −1.     

Non-local diffusion and the chemical structure of molecular clouds

Observations of fluctuations in the redshifted 21-cm radiation from neutral hydrogen (H  i ) are perceived to be an important future probe of the universe at high redshifts. Under the assumption that at redshifts   z ≤ 6  (post-reionization era) the H  i traces the underlying dark matter with a possible bias, we investigate the possibility of using observations of redshifted 21-cm radiation to detect the bispectrum arising from non-linear gravitational clustering and from non-linear bias. We find that the expected signal is ∼ 0.1  mJy at  325  MHz ( z = 3.4)  for the small baselines at the Giant Metrewave Radio Telescope, the strength being a few times larger at higher frequencies  (610 MHz, z = 1.3)  . Further, the magnitude of the signal from the bispectrum is predicted to be comparable to that from the power spectrum, allowing a detection of both in roughly the same integration time. The H  i signal is found to be uncorrelated beyond frequency separations of ∼1.3 MHz whereas the continuum sources of contamination are expected to be correlated across much larger frequencies. This signature can in principle be used to distinguish the H  i signal from the contamination. We also consider the possibility of using observations of the bispectrum to determine the linear and quadratic bias parameters of the H  i at high redshifts, this having possible implications for theories of galaxy formation.     

On using visibility correlations to probe the H i distribution from the dark ages to the present epoch – I. Formalism and the expected signal

Redshifted 21-cm radiation originating from the cosmological distribution of neutral hydrogen (H  i ) appears as background radiation in low-frequency radio observations. The angular and frequency domain fluctuations in this radiation carry information concerning cosmological structure formation. We propose that correlations between visibilities measured at different baselines and frequencies in radio-interferometric observations be used to quantify the statistical properties of these fluctuations. This has an inherent advantage over other statistical estimators in that it deals directly with the visibilities which are the primary quantities measured in radio-interferometric observations. Also, the visibility correlation has a very simple relation with the power spectrum. We present estimates of the expected signal for nearly the entire post-recombination era, from the dark ages to the present epoch. The epoch of reionization, where H  i has a patchy distribution, has a distinct signature where the signal is determined by the size of the discrete ionized regions. The signal at other epochs, where H  i follows the dark matter, is determined largely by the power spectrum of dark matter fluctuations. The signal is strongest for baselines where the antenna separations are within a few hundred times the wavelength of observation, and an optimal strategy would preferentially sample these baselines. In the frequency domain, for most baselines the visibilities at two different frequencies are uncorrelated beyond  Δν∼ 1 MHz  , a signature which, in principle, would allow the H  i signal to be easily distinguished from the continuum sources of contamination.     

Australia Telescope Compact Array H i observations of the NGC 6845 galaxy group

We present the results of Australia Telescope Compact Array (ATCA) H  i line and 20-cm radio continuum observations of the galaxy quartet NGC 6845. The H  i emission extends over all four galaxies but can only be associated clearly with the two spiral galaxies, NGC 6845A and B, which show signs of strong tidal interaction. We derive a total H  i mass of at least  1.8 × 10¹⁰ M_⊙  , most of which is associated with NGC 6845A, the largest galaxy of the group. We investigate the tidal interaction between NGC 6845A and B by studying the kinematics of distinct H  i components and their relation to the known H  ii regions. No H  i emission is detected from the two lenticular galaxies, NGC 6845C and D. A previously uncatalogued dwarf galaxy, ATCA  J2001−4659  , was detected 4.4 arcmin NE from NGC 6845B and has an H  i mass of  ∼5 × 10⁸ M_⊙  . No H  i bridge is visible between the group and its newly detected companion. Extended 20-cm radio continuum emission is detected in NGC 6845A and B as well as in the tidal bridge between the two galaxies. We derive star formation rates of  15–40 M_⊙ yr⁻¹  .     

An efficient technique for pre-selecting low-redshift damped Lyα systems

The number of z ∼ 1 damped Lyα systems (DLAs, log  N (H  i ) ≥ 20.3) per unit redshift is approximately 0.1, making them relatively rare objects. Large, blind QSO surveys for low-redshift DLAs are therefore an expensive prospect for space-borne ultraviolet telescopes. Increasing the efficiency of these surveys by pre-selecting DLA candidates based on the equivalent widths (EWs) of metal absorption lines has previously been a successful strategy. However, the success rate of DLA identification is still only ∼35 per cent when simple EW cut-offs are applied, the majority of systems having 19.0 < log  N (H  i ) < 20.3. Here, we propose a new way to pre-select DLA candidates. Our technique requires high-to-moderate-resolution spectroscopy of the Mg  ii λ2796 transition, which is easily accessible from the ground for 0.2 ≲ z ≲ 2.4. We define the D -index, the ratio of the line equivalent width to velocity spread, and measure this quantity for 19 DLAs and eight sub-DLAs in archival spectra obtained with echelle spectrographs. For the majority of absorbers, there is a clear distinction between the D -index of DLAs compared with sub-DLAs (Kolmogorov–Smirnov probability = 0.8 per cent). Based on this pilot data sample, we find that the D -index can select DLAs with a success rate of up to 90 per cent, an increase in selection efficiency by a factor of 2.5 compared with a simple EW cut. We test the applicability of the D -index at lower resolution and find that it remains a good discriminant of DLAs for full width at half-maximum (FWHM) ≲ 1.5 Å. However, the recommended D -index cut-off between DLAs and sub-DLAs decreases with poorer resolution and we tabulate the appropriate D -index values that should be used with spectra of different resolutions.     

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.     

H i and dark matter in the windy starburst dwarf galaxy NGC 1705

We present 21-cm H  i line observations of the blue compact dwarf galaxy NGC 1705. Previous optical observations show a strong outflow powered by an ongoing starburst dominating the H  ii morphology and kinematics. In contrast, most of the H  i lies in a rotating disc. An extraplanar H  i spur accounts for ∼8 per cent of the total H  i mass, and is possibly associated with the H  ii outflow. The inferred mass loss rate out of the core of the galaxy is significant, ∼0.2 − 2 M_⊙ yr⁻¹, but does not dominate the H  i dynamics. Mass model fits to the rotation curve show that the dark matter (DM) halo is dominant at nearly all radii and has a central density ρ₀ ≈ 0.1 M_⊙ pc⁻³: ten times higher than typically found in dwarf irregular galaxies, but similar to the only other mass-modelled blue compact dwarf, NGC 2915. This large difference strongly indicates that there is little evolution between dwarf irregular and blue compact dwarf types. Instead, dominant DM haloes may regulate the morphology of dwarf galaxies by setting the critical surface density for disc star formation. Neither our data nor catalogue searches reveal any likely external trigger to the starburst in NGC 1705.     

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.     

A study of interstellar medium of dwarf galaxies using H i power spectrum analysis

We estimate the power spectrum of H  i intensity fluctuations for a sample of eight galaxies (seven dwarf and one spiral). The power spectrum can be fitted to a power-law     for six of these galaxies, indicating turbulence is operational. The estimated best-fitting value for the slope ranges from  ∼−1.5  (AND IV, NGC 628, UGC 4459 and GR 8) to  ∼−2.6  (DDO 210 and NGC 3741). We interpret this bi-modality as being due to having effectively 2D turbulence on length-scales much larger than the scale-height of the galaxy disc and 3D otherwise. This allows us to use the estimated slope to set bounds on the scale-heights of the face-on galaxies in our sample. We also find that the power-law slope remains constant as we increase the channel thickness for all these galaxies, suggesting that the fluctuations in H  i intensity are due to density fluctuations and not velocity fluctuations, or that the slope of the velocity structure function is ∼0. Finally, for the four galaxies with '2D turbulence' we find that the slope α correlates with the star formation rate (SFR) per unit area, with larger SFRs leading to steeper power laws. Given our small sample size, this result needs to be confirmed with a larger sample.