Far-infrared detection limits – II. Probing confusion including source confusion

We present a comprehensive analysis for the determination of the confusion levels for the current and the next generation of far-infrared surveys assuming three different cosmological evolutionary scenarios. We include an extensive model for diffuse emission from infrared cirrus in order to derive absolute sensitivity levels taking into account the source confusion noise due to point sources, the sky confusion noise due to the diffuse emission, and instrumental noise. We use our derived sensitivities to suggest best survey strategies for the current and the future far-infrared space missions Spitzer , AKARI ( ASTRO-F ), Herschel and SPICA . We discuss whether the theoretical estimates are realistic and the competing necessities of reliability and completeness. We find the best estimator for the representation of the source confusion and produce predictions for the source confusion using far-infrared source count models. From these confusion limits considering both source and sky confusions, we obtain the optimal, confusion limited redshift distribution for each mission. Finally, we predict the cosmic far-infrared background (CFIRB), which includes information about the number and distribution of the contributing sources.     

The Rhodes/HartRAO 2326-MHz radio continuum survey

We have made scanning observations with the HartRAO 26-m radio telescope to obtain a pencil-beam map of 67 per cent of the sky at 2326 MHz. This is the highest resolution and highest frequency radio continuum map of this type made of such a large area of sky.   In this paper we describe the observations and data reduction procedures used to produce the survey. The resulting map has an angular resolution (HPBW) of 20 arcmin, and the rms pointing accuracy is 1.2 arcmin. The rms noise fluctuations are less than 30-mK T _FB over the whole map. We estimate that the uncertainty in the temperature scale is less than 5 per cent, and that the error in the absolute zero level is better than 80-mK T _FB in any direction.   High-contrast half-tone images of the data with a model of the diffuse galactic background subtracted are presented. These images show many complex emission structures up to and beyond 50° latitude, and illustrate the quality of the data. Extracts from the survey data are available via FTP by arrangement with the authors.     

Extragalactic sources towards the central region of the Galaxy

We have observed a sample of 64 small-diameter sources towards the central  −6° < l < 6°, −2° < b < 2°  of the Galaxy with the aim of studying the Faraday rotation measure near the Galactic Centre region. All the sources were observed at 6- and 3.6-cm wavelengths using the ATCA and the VLA. 59 of these sources are inferred to be extragalactic. The observations presented here constitute the first systematic study of the radio polarization properties of the background sources towards this direction and increase the number of known extragalactic radio sources in this part of the sky by almost an order of magnitude. Based on the morphology, spectral indices and lack of polarized emission, we identify four Galactic H  ii regions in the sample.     

A model of diffuse Galactic radio emission from 10 MHz to 100 GHz

Understanding diffuse Galactic radio emission is interesting both in its own right and for minimizing foreground contamination of cosmological measurements. cosmic microwave background experiments have focused on frequencies ≳10 GHz, whereas 21-cm tomography of the high-redshift universe will mainly focus on ≲0.2 GHz, for which less is currently known about Galactic emission. Motivated by this, we present a global sky model derived from all publicly available total power large-area radio surveys, digitized with optical character recognition when necessary and compiled into a uniform format, as well as the new Villa Elisa data extending the 1.42-GHz map to the entire sky. We quantify statistical and systematic uncertainties in these surveys by comparing them with various global multifrequency model fits. We find that a principal component based model with only three components can fit the 11 most accurate data sets (at 10, 22, 45 and 408 MHz and 1.42, 2.326, 23, 33, 41, 61, 94 GHz) to an accuracy around 1–10 per cent depending on frequency and sky region. Both our data compilation and our software returning a predicted all-sky map at any frequency from 10 MHz to 100 GHz are publicly available at http://space.mit.edu/home/angelica/gsm .     

An Australia Telescope survey for CMB anisotropies

We have surveyed six distinct 'empty fields' using the Australia Telescope Compact Array (ATCA) in an ultracompact configuration with the aim of imaging, with a high brightness sensitivity, any arcminute-scale brightness-temperature anisotropies in the background radio sky. The six well-separated regions were observed at a frequency of 8.7 GHz, and the survey regions were limited by the ATCA primary beams which have a full width at half-maximum of 6 arcmin at this frequency; all fields were observed with a resolution of 2 arcmin and an rms thermal noise of 24 μJy beam⁻¹. After subtracting foreground confusion detected in higher resolution images of the fields, residual fluctuations in Stokes I images are consistent with the expectations from thermal noise and weaker (unidentified) foreground sources; the Stokes Q and U images are consistent with expectations from thermal noise.
Within the sensitivity of our observations, we have no reason to believe that there are any Sunyaev–Zeldovich holes in the microwave sky surveyed. Assuming Gaussian-form CMB anisotropy with a 'flat' spectrum, we derive 95 per cent confidence upper limits of Q _flat<10–11 μK in polarized intensity and Q _flat<25 μK in total intensity. The ATCA filter function peaks at l =4700 and has half-maximum values at l =3350 and 6050.     

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.     

Asteroid Confusions with Extremely Large Telescopes

Asteroids can be considered as sources of contamination of point sources and also sources of confusion noise, depending whether their presence is detected in the image or their flux is under the detection limit. We estimate that at low ecliptic latitudes, ≈10,000–20,000 asteroids/sq. degree will be detected with an E-ELT like telescope, while by the end of Spitzer and Herschel missions, infrared space observatories will provide ≈100,000 serendipitous asteroid detections. The detection and identification of asteroids is therefore an important step in survey astronomy.     

The H i Parkes All Sky Survey: southern observations, calibration and robust imaging

The acquisition of H  i Parkes All Sky Survey (HIPASS) southern sky data commenced at the Australia Telescope National Facility's Parkes 64-m telescope in 1997 February, and was completed in 2000 March. HIPASS is the deepest H  i survey yet of the sky south of declination +2°, and is sensitive to emission out to 170 h₇₅⁻¹ Mpc. The characteristic root mean square noise in the survey images is 13.3 mJy. This paper describes the survey observations, which comprise 23 020 eight-degree scans of 9-min duration, and details the techniques used to calibrate and image the data. The processing algorithms are successfully designed to be statistically robust to the presence of interference signals, and are particular to imaging point (or nearly point) sources. Specifically, a major improvement in image quality is obtained by designing a median-gridding algorithm which uses the median estimator in place of the mean estimator.     

Towards a free–free template for CMB foregrounds

A full-sky template map of the Galactic free–free foreground emission component is increasingly important for high-sensitivity cosmic microwave background (CMB) experiments. We use the recently published Hα data of both the northern and southern skies as the basis for such a template.
The first step is to correct the Hα maps for dust absorption using the 100-μm dust maps of Schlegel, Finkbeiner & Davis. We show that for a range of longitudes, the Galactic latitude distribution of absorption suggests that it is 33 per cent of the full extragalactic absorption. A reliable absorption-corrected Hα map can be produced for ∼95 per cent of the sky; the area for which a template cannot be recovered is the Galactic plane area  | b | < 5°, l = 260°–0°–160°  and some isolated dense dust clouds at intermediate latitudes.
The second step is to convert the dust-corrected Hα data into a predicted radio surface brightness. The free–free emission formula is revised to give an accurate expression (1 per cent) for the radio emission covering the frequency range 100 MHz–100 GHz and the electron temperature range 3000–20 000 K. The main uncertainty when applying this expression is the variation of electron temperature across the sky. The emission formula is verified in several extended H  ii regions using data in the range 408–2326 MHz.
A full-sky free–free template map is presented at 30 GHz; the scaling to other frequencies is given. The Haslam et al. all-sky 408-MHz map of the sky can be corrected for this free–free component, which amounts to a  ≈6  per cent correction at intermediate and high latitudes, to provide a pure synchrotron all-sky template. The implications for CMB experiments are discussed.     

The Australia Telescope search for cosmic microwave background anisotropy

In an attempt to detect cosmic microwave background (CMB) anisotropy on arcmin scales, we have made an 8.7-GHz image of a sky region with a resolution of 2 arcmin and high surface brightness sensitivity using the Australia Telescope Compact Array (ATCA) in an ultracompact configuration. The foreground discrete-source confusion was estimated from observations with higher resolution at the same frequency and in a scaled array at a lower frequency. Following the subtraction of the foreground confusion, the field shows no features in excess of the instrument noise. This limits the CMB anisotropy flat-band power to Q _flat < 23.6 μ K with 95 per cent confidence; the ATCA filter function (which is available at the website www.atnf.csiro.au/Research/cmbr/cmbr_atca.html) F _l in multipole l -space peaks at l _eff = 4700 and has half-maximum values at l  = 3350 and 6050.     

Star formation in NGC 6334 I and I(N)

The northern section of the molecular cloud complex NGC 6334 has been mapped in the CO and CS spectral line emission and in continuum emission at a wavelength of 1300 μm. Our observations highlight the two dominant sources, I and I(N), and a host of weaker sources. NGC 6334 I is associated with a cometary ultracompact H  ii region and a hot, compact core ≤10 arcsec in size. Mid-infrared and CH₃OH observations indicate that it is also associated with at least two protostellar sources, each of which may drive a molecular outflow. For region I we confirm the extreme high-velocity outflow first discovered by Bachiller & Cernicharo and find that it is very energetic with a mechanical luminosity of 390 L_⊙. A dynamical age for the outflow is ∼3000 yr. We also find a weaker outflow originating from the vicinity of NGC 6334 I. In CO and CS this outflow is quite prominent to the north-west, but much less so on the eastern side of I, where there is very little molecular gas. Spectral survey data show a molecular environment at position I which is rich in methanol, methyl formate and dimethyl ether, with lines ranging in energy up to 900 K above the ground state. NGC 6334 I(N) is more dense than I, but cooler, and has none of the high-excitation lines observed toward I. I(N) also has an associated outflow, but it is less energetic than the outflow from I. The fully sampled continuum map shows a network of filaments, voids and cores, many of which are likely to be sites of star formation. A striking feature is a narrow, linear ridge which defines the western boundary. It is unclear if there is a connection between this filament and the many potential sites of star formation, or if the filament existed prior to the star formation activity.     

Applications of wavelets to the analysis of cosmic microwave background maps

We consider wavelets as a tool to perform a variety of tasks in the context of analysing cosmic microwave background (CMB) maps. Using spherical Haar wavelets, we define a position and angular-scale-dependent measure of power that can be used to assess the existence of spatial structure. We apply planar Daubechies wavelets for the identification and removal of point sources from small sections of sky maps. Our technique can successfully identify virtually all point sources that are above 3 and more than 80 per cent of those above 1 . We discuss the trade-offs between the levels of correct and false detections. We denoise and compress a 100 000-pixel CMB map by a factor of 10 in 5 s, achieving a noise reduction of about 35 per cent. In contrast to Wiener filtering, the compression process is model-independent and very fast. We discuss the usefulness of wavelets for power spectrum and cosmological-parameter estimation. We conclude that at present wavelet functions are most suitable for identifying localized sources.     

MERLIN polarization observations of compact steep-spectrum sources at 5 GHz

We present polarization observations of 28 compact steep-spectrum sources made with the upgraded MERLIN at 5 GHz. With an angular resolution of 60 milliarcsec and rms noise levels of about 0.1 mJy beam⁻¹, the total intensity images reveal new details in many of these sources. A few of the more extended lobes and jets are quite highly polarized, but more than half of the components are completely unpolarized. Comparison with published data implies that this is due to Faraday depolarization, probably occurring in the surrounding medium with nB  ∼ 1 cm^−3μG. The high resolution of the present observations implies that the variations in Faraday rotation, probably due to magnetic field tangling, occur on scales of less than about 100 pc.     

Too large and overlooked? Extended free–free emission towards massive star formation regions

We present Australia Telescope Compact Array observations towards six massive star formation regions, which, from their strong 24 GHz continuum emission but no compact 8 GHz continuum emission, appeared good candidates for hypercompact H  ii regions. However, the properties of the ionized gas derived from the 19 to 93 GHz continuum emission and  H70α+ H57α  radio recombination line data show the majority of these sources are, in fact, regions of spatially extended, optically thin free–free emission. These extended sources were missed in the previous 8 GHz observations due to a combination of spatial filtering, poor surface brightness sensitivity and primary beam attenuation.
We consider the implications that a significant number of these extended H  ii regions may have been missed by previous surveys of massive star formation regions. If the original sample of 21 sources is representative of the population as a whole, the fact that six contain previously undetected extended free–free emission suggests a large number of regions have been mis-classified. Rather than being very young objects prior to UCH  ii region formation, they are, in fact, associated with extended H  ii regions and thus significantly older. In addition, inadvertently ignoring a potentially substantial flux contribution (up to ∼0.5 Jy) from free–free emission has implications for dust masses derived from sub-mm flux densities. The large spatial scales probed by single-dish telescopes, which do not suffer from spatial filtering, are particularly susceptible and dust masses may be overestimated by up to a factor of ∼2.     

On the infrared void in the Lupus dark clouds

Strömgren uvbyβ photometry observations obtained for 205 stars in the general direction of a void in the IRAS 100-μm emission from the Lupus dark cloud complex are presented and analysed. The colour excess versus distance diagram confirms the existence of a region depleted from interstellar material, which is also seen in the ROSAT soft X-ray background emission map. The distance to the surrounding material is estimated as being within the interval from 60 to 100 pc. This result is in disagreement with previous distance estimates to the supposed supernova that has been suggested as responsible for clearing the region from dust. As an alternative, the data presented support the suggestion that the void may have been produced by the detachment of material from the interface between Loop I and the Local Bubble as a consequence of hydromagnetic instabilities. Moreover, the distribution of colour excess as a function of distance supports a value of ∼150 pc as the most probable distance to the dark cloud known as Lupus 1.     

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.     

Wide-field imaging and polarimetry for the biggest and brightest in the 20-GHz southern sky

We present the wide-field imaging and polarimetry at  ν= 20 GHz  of seven most extended, bright  ( S _total≥ 0.50 Jy)  , high-frequency selected radio sources in the southern sky with declinations  δ < −30°  . Accompanying the data are brief reviews of the literature for each source. The results presented here aid in the statistical completeness of the Australia Telescope 20-GHz Survey: the Bright Source Sample. The data are of crucial interest for future cosmic microwave background missions as a collection of information about candidate calibrator sources. We were able to obtain data for seven of the nine sources identified by our selection criteria. We report that Pictor A is thus far the best extragalactic calibrator candidate for the Low Frequency Instrument of the Planck European Space Agency mission due to its high level of integrated polarized flux density  (∼0.50 ± 0.06 Jy)  on a scale of 10 arcmin. Six out of the seven sources have a clearly detected compact radio core in our images, with either a null detection or less than 2 per cent detection of polarized emission from the nuclei. Most sources with detected jets have magnetic field alignments running in a longitudinal configuration, however, PKS 1333−33 exhibits transverse fields and an orthogonal change in field geometry from nucleus to jets.     

MAMBO observations at 240 GHz of optically obscured Spitzer sources: source clumps and radio activity at high redshift

Optically very faint  ( R > 25.5)  sources detected by the Spitzer Space Telescope at 24 μm represent a very interesting population at redshift   z ∼ (1.5–3)  . They exhibit strong clustering properties, implying that they are hosted by very massive haloes, and their mid-infrared emission could be powered by either dust-enshrouded star formation and/or by an obscured active galactic nucleus (AGN). We report observations carried out with the Max Planck Millimetre Bolometer (MAMBO) array at the IRAM 30-m antenna on Pico Veleta of a candidate protocluster with five optically obscured sources selected from the 24-μm Spitzer sample of the First-Look Survey. Interestingly, these sources appear to lie on a high-density filament aligned with the two radio jets of an AGN. Four out of five of the observed sources were detected. We combine these measurements with optical, infrared and radio observations to probe the nature of the candidate protocluster members. Our preliminary conclusions can be summarized as follows: the spectral energy distributions (SEDs) of all sources include both AGN and starburst contributions; the AGN contribution to the bolometric luminosities ranges between 14 and 26 per cent of the total. Such a contribution is enough for the AGN to dominate the emission at 5.8, 8 and 24 μm, while the stellar component, inferred from SED fitting, prevails at 1.25 mm and at  λ < 4.5 μ  m. The present analysis suggests a coherent interplay at high z between extended radio activity and the development of filamentary large-scale structures.     

Estimating the spectral indices of correlated astrophysical foregrounds by a second-order statistical approach

We present the first tests of a new method, the correlated component analysis (CCA) based on second-order statistics, to estimate the mixing matrix, a key ingredient to separate astrophysical foregrounds superimposed to the Cosmic Microwave Background (CMB). In the present application, the mixing matrix is parametrized in terms of the spectral indices of Galactic synchrotron and thermal dust emissions, while the free–free spectral index is prescribed by basic physics, and is thus assumed to be known. We consider simulated observations of the microwave sky with angular resolution and white stationary noise at the nominal levels for the Planck satellite, and realistic foreground emissions, with a position-dependent synchrotron spectral index. We work with two sets of Planck frequency channels: the low-frequency set, from 30 to 143 GHz, complemented with the Haslam 408 MHz map, and the high-frequency set, from 217 to 545 GHz. The concentration of intense free–free emission on the Galactic plane introduces a steep dependence of the spectral index of the global Galactic emission with Galactic latitude, close to the Galactic equator. This feature makes difficult for the CCA to recover the synchrotron spectral index in this region, given the limited angular resolution of Planck , especially at low frequencies. A cut of a narrow strip around the Galactic equator  (| b | < 3°)  , however, allows us to overcome this problem. We show that, once this strip is removed, the CCA allows an effective foreground subtraction, with residual uncertainties inducing a minor contribution to errors on the recovered CMB power spectrum.     

Space–VLBI observations of the OH maser OH 34.26+0.15: low interstellar scattering

We report on the first space–VLBI observations of the OH 34.26+0.15 maser in two main-line OH transitions at 1665 and 1667 MHz. The observations involved the space radio telescope on board the Japanese satellite HALCA and an array of ground radio telescopes. The map of the maser region and images of individual maser spots were produced with an angular resolution of 1 mas, which is several times higher than the angular resolution available on the ground. The maser spots were only partly resolved and a lower limit to the brightness temperature     was obtained. The maser seems to be located in the direction of low interstellar scattering, an order of magnitude lower than the scattering of a nearby extragalactic source and pulsar.