No Arabic abstract
Here I will outline successes and challenges for finding spectral line sources in large data cubes that are dominated by noise. This is a 3D challenge as the sources we wish to catalog are spread over several spatial pixels and spectral channels. While 2D searches can be applied, e.g., channel by channel, optimal searches take into account the 3-dimensional nature of the sources. In this overview I will focus on HI 21-cm spectral line source detection in extragalactic surveys, in particular HIPASS, the HI Parkes All-Sky Survey and WALLABY, the ASKAP HI All-Sky Survey. I use the original HIPASS data to highlight the diversity of spectral signatures of galaxies and gaseous clouds, both in emission and absorption. Among others, I report the discovery of a 680 km/s wide HI absorption trough in the megamaser galaxy NGC 5793. Issues such as source confusion and baseline ripples, typically encountered in single-dish HI surveys, are much reduced in interferometric HI surveys. Several large HI emission and absorption surveys are planned for the Australian Square Kilometre Array Pathfinder (ASKAP): here we focus on WALLABY, the 21-cm survey of the sky (Dec < +30 degr; z < 0.26) which will take about one year of observing time with ASKAP. Novel phased array feeds (radio cameras) will provide 30 square degrees instantaneous field-of-view. WALLABY is expected to detect more than 500 000 galaxies, unveil their large-scale structures and cosmological parameters, detect their extended, low-surface brightness disks as well as gas streams and filaments between galaxies. It is a precursor for future HI surveys with SKA Phase I and II, exploring galaxy formation and evolution. The compilation of highly reliable and complete source catalogs will require sophisticated source-finding algorithms as well as accurate source parametrisation.
The large spectral bandwidth and wide field of view of the Australian SKA Pathfinder radio telescope will open up a completely new parameter space for large extragalactic HI surveys. Here we focus on identifying and parametrising HI absorption lines which occur in the line of sight towards strong radio continuum sources. We have developed a method for simultaneously finding and fitting HI absorption lines in radio data by using multi-nested sampling, a Bayesian Monte Carlo algorithm. The method is tested on a simulated ASKAP data cube, and is shown to be reliable at detecting absorption lines in low signal-to-noise data without the need to smooth or alter the data. Estimation of the local Bayesian evidence statistic provides a quantitative criterion for assigning significance to a detection and selecting between competing analytical line-profile models.
We report conclusive verification of the detection of associated HI 21 cm absorption in the early-type host galaxy of the compact radio source PMNJ2054-4242. We estimate an effective spectral line velocity width of 418 +/- 20 km s^{-1} and observed peak optical depth of 2.5 +/- 0.2 per cent, making this one of the broadest and weakest 21 cm absorption lines yet detected. For T_{spin}/f > 100 K the atomic neutral hydrogen column density is N_{HI} > 2 x 10^{21} cm^{-2}. The observed spectral line profile is redshifted by 187 +/- 46 km s^{-1}, with respect to the optical spectroscopic measurement, perhaps indicating that the HI gas is infalling towards the central active galactic nucleus. Our initial tentative detection would likely have been dismissed by visual inspection, and hence its verification here is an excellent test of our spectral line detection technique, currently under development in anticipation of future next-generation 21 cm absorption-line surveys.
We introduce SoFiA, a flexible software application for the detection and parameterization of sources in 3D spectral-line datasets. SoFiA combines for the first time in a single piece of software a set of new source-finding and parameterization algorithms developed on the way to future HI surveys with ASKAP (WALLABY, DINGO) and APERTIF. It is designed to enable the general use of these new algorithms by the community on a broad range of datasets. The key advantages of SoFiA are the ability to: search for line emission on multiple scales to detect 3D sources in a complete and reliable way, taking into account noise level variations and the presence of artefacts in a data cube; estimate the reliability of individual detections; look for signal in arbitrarily large data cubes using a catalogue of 3D coordinates as a prior; provide a wide range of source parameters and output products which facilitate further analysis by the user. We highlight the modularity of SoFiA, which makes it a flexible package allowing users to select and apply only the algorithms useful for their data and science questions. This modularity makes it also possible to easily expand SoFiA in order to include additional methods as they become available. The full SoFiA distribution, including a dedicated graphical user interface, is publicly available for download.
The Evolutionary Map of the Universe (EMU) is a proposed radio continuum survey of the Southern Hemisphere up to declination +30 deg., with the Australian Square Kilometre Array Pathfinder (ASKAP). EMU will use an automated source identification and measurement approach that is demonstrably optimal, to maximise the reliability, utility and robustness of the resulting radio source catalogues. As part of the process of achieving this aim, a Data Challenge has been conducted, providing international teams the opportunity to test a variety of source finders on a set of simulated images. The aim is to quantify the accuracy of existing automated source finding and measurement approaches, and to identify potential limitations. The Challenge attracted nine independent teams, who tested eleven different source finding tools. In addition, the Challenge initiators also tested the current ASKAPsoft source-finding tool to establish how it could benefit from incorporating successful features of the other tools. Here we present the results of the Data Challenge, identifying the successes and limitations for this broad variety of the current generation of radio source finding tools. As expected, most finders demonstrate completeness levels close to 100% at 10sigma dropping to levels around 10% by 5sigma. The reliability is typically close to 100% at 10sigma, with performance to lower sensitivities varying greatly between finders. All finders demonstrate the usual trade-off between completeness and reliability, whereby maintaining a high completeness at low signal-to-noise comes at the expense of reduced reliability, and vice-versa. We conclude with a series of recommendations for improving the performance of the ASKAPsoft source-finding tool.
Accurate optical redshifts will be critical for spectral co-adding techniques used to extract detections from below the noise level in ongoing and upcoming surveys for HI, which will extend our current understanding of gas reservoirs in galaxies to lower column densities and higher redshifts. We have used existing, high quality optical and radio data from the SDSS and ALFALFA surveys to investigate the relationship between redshifts derived from optical spectroscopy and neutral hydrogen (HI) spectral line observations. We find that the two redshift measurements agree well, with a negligible systematic offset and a small distribution width. Employing simple simulations, we determine how the width of an ideal stacked HI profile depends on these redshift offsets, as well as larger redshift errors more appropriate for high redshift galaxy surveys. The width of the stacked profile is dominated by the width distribution of the input individual profiles when the redshift errors are less than the median width of the input profiles, and only when the redshift errors become large, ~150 km/s, do they significantly affect the width of the stacked profile. This redshift accuracy can be achieved with moderate resolution optical spectra. We provide guidelines for the number of spectra required for stacking to reach a specified mass sensitivity, given telescope and survey parameters, which will be useful for planning optical spectroscopy observing campaigns to supplement the radio data.