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.     

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.     

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.     

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.     

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.     

Primordial non-Gaussianities in the intergalactic medium

We present results from the first high-resolution hydrodynamical simulations of non-Gaussian cosmological models. We focus on the statistical properties of the transmitted Lyman-α flux in the high-redshift intergalactic medium. Imprints of non-Gaussianity are present and are larger at high redshifts. Differences larger than 20 per cent at   z > 3  in the flux probability distribution function for high-transmissivity regions (voids) are expected for values of the non-linearity parameter   f _NL=±100  when compared to a standard Λ cold dark matter cosmology with   f _NL= 0  . We also investigate the one-dimensional flux bispectrum: at the largest scales (corresponding to tens of Mpc), we expect deviations in the flux bispectrum up to 20 per cent at   z ∼ 4  (for   f _NL=±100  ), significantly larger than deviations of ∼3 per cent in the flux power spectrum. We briefly discuss possible systematic errors that can contaminate the signal. Although challenging, a detection of non-Gaussianities in the interesting regime of scales and redshifts probed by the Lyman-α forest could be possible with future data sets.     

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.     

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.     

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.     

Improved constraints on possible variation of physical constants from H i 21-cm and molecular QSO absorption lines

Quasar (QSO) absorption spectra provide an extremely useful probe of possible cosmological variation in various physical constants. Comparison of H  i 21-cm absorption with corresponding molecular (rotational) absorption spectra allows us to constrain variation in     , where α is the fine-structure constant and g _p is the proton g -factor. We analyse spectra of two QSOs, PKS 1413+135 and TXS 0218+357, and derive values of     at absorption redshifts of     and 0.6847 by simultaneous fitting of the H  i 21-cm and molecular lines. We find     and     respectively, indicating an insignificantly smaller y in the past. We compare our results with other constraints from the same two QSOs given recently by Drinkwater et al. and Carilli et al., and with our recent optical constraints, which indicated a smaller α at higher redshifts.     

High-resolution H i and radio continuum observations of the SNR G290.1−0.8

We have observed the supernova remnant (SNR) G290.1−0.8 in the 21-cm H  i line and the 20-cm radio continuum using the Australia Telescope Compact Array (ATCA). The H  i data were combined with data from the Southern Galactic Plane Survey to recover the shortest spatial frequencies. In contrast, H  i absorption was analysed by filtering extended H  i emission, with spatial frequencies shorter than 1.1 kλ. The low-resolution ATCA radio continuum image of the remnant shows considerable internal structure, resembling a network of filaments across its 13-arcmin diameter. A high-resolution ATCA radio continuum image was also constructed to study the small-scale structure in the SNR. It shows that there are no structures smaller than ∼17 arcsec, except perhaps for a bright knot to the south, which is probably an unrelated object. The H  i absorption study shows that the gas distribution and kinematics in front of SNR G290.1−0.8 are complex. We estimate that the SNR probably lies in the Carina arm, at a distance 7 (±1) kpc. In addition, we have studied nearby sources in the observed field using archival multiwavelength data to determine their characteristics.     

The 'Carina Flare' supershell: probing the atomic and molecular ISM in a Galactic chimney

The 'Carina Flare' supershell, GSH 287+04−17, is a molecular supershell originally discovered in  ¹²CO( J = 1–0)  with the NANTEN 4 m telescope. We present the first study of the shell's atomic ISM, using H  i 21-cm line data from the Parkes 64-m telescope Southern Galactic Plane Survey. The data reveal a gently expanding,  ∼230 × 360  pc H  i supershell that shows strong evidence of Galactic Plane blowout, with a break in its main body at   z ∼ 280  pc and a capped high-latitude extension reaching   z ∼ 450  pc. The molecular clouds form comoving parts of the atomic shell, and the morphology of the two phases reflects the supershell's influence on the structure of the ISM. We also report the first discovery of an ionized component of the supershell, in the form of delicate, streamer-like filaments aligned with the proposed direction of blowout. The distance estimate to the shell is re-examined, and we find strong evidence to support the original suggestion that it is located in the Carina Arm at a distance of  2.6 ± 0.4 kpc  . Associated H  i and H₂ masses are estimated as   M _H I≈ 7 ± 3 × 10⁵ M_⊙  and     , and the kinetic energy of the expanding shell as   E _K ∼ 1 × 10⁵¹  erg. We examine the results of analytical and numerical models to estimate a required formation energy of several 10⁵¹ to  ∼10⁵²  erg, and an age of  ∼10⁷ yr  . This age is compatible with molecular cloud formation time-scales, and we briefly consider the viability of a supershell-triggered origin for the molecular component.     

H i content and star formation in the interacting galaxy Arp86

We present the results of Giant Metrewave Radio Telescope (GMRT) observations of the interacting system Arp86 in both neutral atomic hydrogen, H  i , and in radio continuum at 240 606 and 1394 MHz. In addition to H  i emission from the two dominant galaxies, NGC 7752 and NGC 7753, these observations show a complex distribution of H  i tails and bridges due to tidal interactions. The regions of highest column density appear related to the recent sites of intense star formation. H  i column densities  ∼1–1.5 × 10²¹ cm⁻²  have been detected in the tidal bridge which is bright in Spitzer image as well. We also detect H  i emission from the galaxy 2MASX J23470758+2926531, which is shown to be a part of this system. We discuss the possibility that this could be a tidal dwarf galaxy. The radio continuum observations show evidence of a non-thermal bridge between NGC 7752 and NGC 7753, and a radio source in the nuclear region of NGC 7753 consistent with it having a low-ionization nuclear emission region nucleus.     

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⁻¹  .     

The redshift distribution of FIRST radio sources at 1 mJy

We present spectra for a sample of radio sources from the FIRST survey, and use them to define the form of the redshift distribution of radio sources at mJy levels. We targeted 365 sources and obtained 46 redshifts (13 per cent of the sample). We find that our sample is complete in redshift measurement to R ∼18.6, corresponding to z ∼0.2. Galaxies were assigned spectral types based on emission-line strengths. Early-type galaxies represent the largest subset (45 per cent) of the sample and have redshifts 0.15≲ z ≲0.5; late-type galaxies make up 15 per cent of the sample and have redshifts 0.05≲ z ≲0.2; starbursting galaxies are a small fraction (∼6 per cent), and are very nearby ( z ≲0.05). Some 9 per cent of the population have Seyfert 1/quasar-type spectra, all at z ≳0.8, and 4 per cent are Seyfert 2 type galaxies at intermediate redshifts ( z ∼0.2).
Using our measurements and data from the Phoenix survey (Hopkins et al.), we obtain an estimate for N ( z ) at S _1.4 GHz≥1 mJy and compare this with model predictions. At variance with previous conclusions, we find that the population of starbursting objects makes up ≲5 per cent of the radio population at S ∼1 mJy.     

Variable 21-cm absorption at z=0.3127

We report multi-epoch Giant Metrewave Radio Telescope (GMRT) H  i observations of the z  = 0.3127 damped absorber towards the quasar PKS 1127−145, which reveal variability in both the absorption profile and the flux of the background source, over a time-scale of a few days.
The observed variations cannot be explained by simple interstellar scintillation (ISS) models where there are only one or two scintillating components and all of the ISS occurs in the Galaxy. More complicated models, where either there are more scintillating components or some of the ISS occurs in the interstellar medium of the z =0.3127 absorber, may be acceptable. However, the variability can probably be best explained in models incorporating motion (on sub-VLBI scales) of a component of the background continuum source, with or without some ISS.
All models producing the variable 21-cm absorption profile require small-scale variations in the 21-cm optical depth of the absorber. The length-scale for the opacity variations is ∼0.1 pc in pure superluminal motion models, and ∼10 pc in pure ISS models. Models involving subluminal motion, combined with scintillation of the moving component, require opacity variations on far smaller scales of ∼ 10–100 au .     

Galactic H i on the 50-au scale in the direction of three extragalactic sources observed with MERLIN

We present MERLIN observations of Galactic 21-cm H  i absorption at an angular resolution of  ∼0.1–0.2  arcsec and a velocity resolution of 0.5 km s⁻¹, in the direction of three moderately low latitude  (−8° < b < −12°)  extragalactic radio sources, 3C 111, 3C 123 and 3C 161, all of which are heavily reddened. H  i absorption is observed against resolved background emission sources up to ∼2 arcsec in extent and we distinguish details of the opacity distribution within 1–1.5 arcsec regions towards 3C 123 and 3C 161. This study is the second MERLIN investigation of small-scale structure in interstellar H  i (earlier work probed Galactic H  i in the directions of the compact sources 3C 138 and 3C 147). The 0.1-arcsec scale is intermediate between H  i absorption studies made with other fixed element interferometers with resolution of 1–10 arcsec and very long baseline interferometry studies with resolutions of 10–20 mas. At a scale of 1 arcsec (about 500 au), prominent changes in Galactic H  i opacity in excess of 1–1.5 are determined in the direction of 3C 161 with a signal-to-noise ratio of at least 10σ. Possible fluctuations in the H  i opacity at the level of about 1 are detected at the  2.5–3σ  level in the direction of 3C 123.     

The AGN and gas disc in the low surface brightness galaxy PGC 045080

We present radio observations and optical spectroscopy of the giant low surface brightness (LSB) galaxy PGC 045080 (or 1300+0144). PGC 045080 is a moderately distant galaxy having a highly inclined optical disc and massive H  i gas content. Radio continuum observations of the galaxy were carried out at 320, 610 MHz and 1.4 GHz. Continuum emission was detected and mapped in the galaxy. The emission appears extended over the inner disc at all three frequencies. At 1.4 GHz and 610 MHz it appears to have two distinct lobes. We also did optical spectroscopy of the galaxy nucleus; the spectrum did not show any strong emission lines associated with active galactic nucleus (AGN) activity but the presence of a weak AGN cannot be ruled out. Furthermore, comparison of the Hα flux and radio continuum at 1.4 GHz suggests that a significant fraction of the emission is non-thermal in nature. Hence we conclude that a weak or hidden AGN may be present in PGC 045080. The extended radio emission represents lobes/jets from the AGN. These observations show that although LSB galaxies are metal poor and have very little star formation, their centres can host significant AGN activity. We also mapped the H  i gas disc and velocity field in PGC 045080. The H  i disc extends well beyond the optical disc and appears warped. In the H  i intensity maps, the disc appears distinctly lopsided. The velocity field is disturbed on the lopsided side of the disc but is fairly uniform in the other half. We derived the H  i rotation curve for the galaxy from the velocity field. The rotation curve has a flat rotation speed of ∼190 km s⁻¹.     

Contribution of galaxies to the background hydrogen-ionizing flux

We estimate the evolution of the contribution of galaxies to the cosmic background flux at 912 Å by means of a semi-analytic model of galaxy formation and evolution. Such modelling has been quite successful in reproducing the optical properties of galaxies. We assume that high-redshift damped Lyman α systems are the progenitors of present-day galaxies, and we design a series of models that are consistent with the evolution of cosmic comoving emissivities in the available near-infrared, optical, ultraviolet and far-infrared bands along with the evolution of the neutral hydrogen content and average metallicity of damped Lyman α systems. We use these models to compute the galactic contribution to the Lyman-limit emissivity and background flux for 0 ≃  z  ≤ 4. We take into account the absorption of Lyman-limit photons by H  I and dust in the interstellar medium of the galaxies. We find that the background Lyman-limit flux due to galaxies might dominate (or be comparable to) the contribution from quasars at almost all redshifts if the absorption by H  I in the interstellar medium is neglected. Such H  I absorption would result in a severe diminishing of this flux — by almost three orders of magnitude at high redshifts and by one to two orders at z  ≃ 0. Though the resulting galaxy flux is completely negligible at high redshifts, it is comparable to the quasar flux at z  ≃ 0.