Isotropic wavelets: a powerful tool to extract point sources from cosmic microwave background maps

It is the aim of this paper to introduce the use of isotropic wavelets to detect and determine the flux of point sources appearing in cosmic microwave background (CMB) maps. The most suitable wavelet to detect point sources filtered with a Gaussian beam is the 'Mexican Hat'. An analytical expression of the wavelet coefficient obtained in the presence of a point source is provided and used in the detection and flux estimation methods presented. For illustration the method is applied to two simulations (assuming Planck mission characteristics) dominated by CMB (100 GHz) and dust (857 GHz), as these will be the two signals dominating at low and high frequencies respectively in the Planck channels. We are able to detect bright sources above 1.58 Jy at 857 GHz (82 per cent of all sources) and above 0.36 Jy at 100 GHz (100 per cent of all), with errors in the flux estimation below 25 per cent. The main advantage of this method is that nothing has to be assumed about the underlying field, i.e. about the nature and properties of the signal plus noise present in the maps. This is not the case in the detection method presented by Tegmark & Oliveira-Costa. Both methods are compared, producing similar results.     

Detection/estimation of the modulus of a vector. Application to point-source detection in polarization data

Given a set of images, whose pixel values can be considered as the components of a vector, it is interesting to estimate the modulus of such a vector in some localized areas corresponding to a compact signal. For instance, the detection/estimation of a polarized signal in compact sources immersed in a background is relevant in some fields like astrophysics. We develop two different techniques, one based on the Neyman–Pearson lemma, the Neyman–Pearson filter (NPF), and another based on pre-filtering before fusion, the filtered fusion (FF), to deal with the problem of detection of the source and estimation of the polarization given two or three images corresponding to the different components of polarization (two for linear polarization, three including circular polarization). For the case of linear polarization, we have performed numerical simulations on two-dimensional patches to test these filters following two different approaches (a blind and a non-blind detection), considering extragalactic point sources immersed in cosmic microwave background (CMB) and non-stationary noise with the conditions of the 70 GHz Planck channel. The FF outperforms the NPF, especially for low fluxes. We can detect with the FF extragalactic sources in a high noise zone with fluxes      Jy for (blind/non-blind) detection and in a low noise zone with fluxes      Jy for (blind/non-blind) detection with low errors in the estimated flux and position.     

The European Large Area ISO Survey: optical identifications of 15-μm and 1.4-GHz sources in N1 and N2

We present the multiwavelength properties and catalogue of the 15 μm and 1.4 GHz radio sources detected in the European Large Area ISO Survey ( ELAIS ) areas N1 and N2. Using the optical data from the Wide Field Survey we use a likelihood ratio method to search for the counterparts of the 1056 and 691 sources detected at 15 μm and 1.4 GHz, respectively, down to flux limits of   S ₁₅= 0.5 mJy  and   S _1.4 GHz= 0.135 mJy  . We find that ∼92 per cent of the 15 μm ELAIS sources have an optical counterpart down to   r '= 24  . All mid-infrared (IR) sources with fluxes   S ₁₅≥ 3 mJy  have an optical counterpart. The magnitude distribution of the sources shows a well-defined peak at relatively bright magnitudes   r '∼ 18  . The mid-IR-to-optical and radio-to-optical flux diagrams are presented and discussed in terms of actual galaxy models. About 15 per cent of the sources are bright galactic stars; of the extragalactic objects ∼65 per cent are compatible with being normal or starburst galaxies and ∼25 per cent active galactic nuclei (AGNs). Objects with mid-IR-to-optical fluxes larger than 100 are found, comprising ∼20 per cent of the sample. We suggest that that these sources are highly obscured luminous and ultraluminous starburst galaxies and AGNs.     

Radio imaging of the Subaru/XMM–Newton Deep Field – I. The 100-μJy catalogue, optical identifications, and the nature of the faint radio source population

We describe deep radio imaging at 1.4 GHz of the 1.3-deg² Subaru/ XMM–Newton Deep Field (SXDF), made with the Very Large Array in B and C configurations. We present a radio map of the entire field, and a catalogue of 505 sources covering 0.8 deg² to a peak flux density limit of 100 μJy. Robust optical identifications are provided for 90 per cent of the sources, and suggested identifications are presented for all but 14 (of which seven are optically blank, and seven are close to bright contaminating objects). We show that the optical properties of the radio sources do not change with flux density, suggesting that active galactic nuclei (AGN) continue to contribute significantly at faint flux densities. We test this assertion by cross-correlating our radio catalogue with the X-ray source catalogue and conclude that radio-quiet AGN become a significant population at flux densities below 300 μJy, and may dominate the population responsible for the flattening of the radio source counts if a significant fraction of them are Compton-thick.     

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.     

The Australia Telescope 20-GHz (AT20G) Survey: the Bright Source Sample

The Australia Telescope 20-GHz (AT20G) Survey is a blind survey of the whole southern sky at 20 GHz (with follow-up observations at 4.8 and 8.6 GHz) carried out with the Australia Telescope Compact Array from 2004 to 2007.
The Bright Source Sample (BSS) is a complete flux-limited subsample of the AT20G Survey catalogue comprising 320 extragalactic     ) radio sources south of  δ=−15°  with      Jy. Of these, 218 have near simultaneous observations at 8 and 5 GHz.
In this paper we present an analysis of radio spectral properties in total intensity and polarization, size, optical identifications and redshift distribution of the BSS sources. The analysis of the spectral behaviour shows spectral curvature in most sources with spectral steepening that increases at higher frequencies (the median spectral index α, assuming   S ∝ν^α  , decreases from  α^8.6_4.8= 0.11  between 4.8 and 8.6 GHz to  α²⁰_8.6=−0.16  between 8.6 and 20 GHz), even if the sample is dominated by flat spectra sources (85 per cent of the sample has  α²⁰_8.6 > −0.5)  . The almost simultaneous spectra in total intensity and polarization allowed us a comparison of the polarized and total intensity spectra: polarized fraction slightly increases with frequency, but the shapes of the spectra have little correlation. Optical identifications provided an estimation of redshift for 186 sources with a median value of 1.20 and 0.13, respectively, for QSO and galaxies.     

Blind and non-blind source detection in WMAP 5-yr maps

We have analysed the efficiency in source detection and flux density estimation of blind and non-blind detection techniques exploiting the MHW2 filter applied to the Wilkinson Microwave Anisotropy Probe ( WMAP ) 5-yr maps. A comparison with the AT20G bright source sample, with a completeness limit of 0.5 Jy and accurate flux measurements at 20 GHz, close to the lowest frequency of WMAP maps, has allowed us to assess the completeness and the reliability of the samples detected with the two approaches, as well as the accuracy of flux and error estimates, and their variations across the sky. The uncertainties on flux estimates given by our procedure turned out to be about a factor of 2 lower than the rms differences with AT20G measurements, consistent with the smoothing of the fluctuation field yielded by map filtering. Flux estimates were found to be essentially unbiased except that, close to the detection limit, a substantial fraction of fluxes are found to be inflated by the contribution of underlying positive fluctuations. This is consistent with expectations for the Eddington bias associated to the true errors on flux density estimates. The blind and non-blind approaches are found to be complementary: each of them allows the detection of sources missed by the other. Combining results of the two methods on the WMAP 5-yr maps, we have expanded the non-blindly generated New Extragalactic WMAP Point Source (NEWPS) catalogue that was based on WMAP 3-yr maps. After having removed the probably spurious objects not identified with known radio sources, the new version of the NEWPS catalogue, NEWPS_5yr comprises 484 sources detected with a signal-to-noise ratio  SNR ≥ 5  .     

Searching for large-scale structure in deep radio surveys

We calculate the expected amplitude of the dipole and higher spherical harmonics in the angular distribution of radio galaxies. The median redshift of radio sources in existing catalogues is z  ∼ 1, which allows us to study large-scale structure on scales between those accessible to present optical and infrared surveys, and that of the cosmic microwave background (CMB). The dipole is a result of two effects which turn out to be of comparable magnitude: (i) our motion with respect to the CMB, and (ii) large-scale structure, parametrized here by a family of cold dark matter power-spectra. We make specific predictions for the Green Bank 1987 (87GB) and Parkes–MIT–NRAO (PMN) catalogues, which in our combined catalogue include ∼ 40 000 sources brighter than 50 mJy at 4.85 GHz, over about 70 per cent of the sky. For these relatively sparse catalogues both the motion and large-scale structure dipole effects are expected to be smaller than the Poisson shot noise. However, we detect dipole and higher harmonics in the combined 87GB–PMN_raw catalogue which are far larger than expected. We attribute this to a 2 per cent flux mismatch between the two catalogues. Ad hoc corrections made in an effort to match the catalogues may suggest a marginal detection of a dipole. To detect a dipole and higher harmonics unambiguously, a catalogue with full sky coverage and ∼ 10⁶ sources is required. We also investigate the existence and extent of the supergalactic plane in the above catalogues. In a strip of ± 10° of the standard supergalactic equator, we find a 3 σ detection in PMN_raw, but only 1 σ in 87 GB_raw. We briefly discuss the implications of ongoing surveys such as FIRST and NVSS and follow-up redshift surveys.     

Timeline analysis and wavelet multiscale analysis of the AKARI All-Sky Survey at 90 μm

We present a careful analysis of the point-source detection limit of the AKARI All-Sky Survey in the WIDE-S 90-μm band near the North Ecliptic Pole (NEP). Timeline analysis is used to detect IRAS ( Infrared Astronomy Satellite ) sources and then a conversion factor is derived to transform the peak timeline signal to the interpolated 90-μm flux of a source. Combined with a robust noise measurement, the point-source flux detection limit at signal-to-noise ratio  (S/N) > 5  for a single detector row is  1.1 ± 0.1 Jy  which corresponds to a point-source detection limit of the survey of ∼0.4 Jy.
Wavelet transform offers a multiscale representation of the Time Series Data ( tsd ). We calculate the continuous wavelet transform of the tsd and then search for significant wavelet coefficients considered as potential source detections. To discriminate real sources from spurious or moving objects, only sources with confirmation are selected. In our multiscale analysis, IRAS sources selected above 4σ can be identified as the only real sources at the Point Source Scales. We also investigate the correlation between the non- IRAS sources detected in timeline analysis and cirrus emission using wavelet transform and contour plots of wavelet power spectrum. It is shown that the non- IRAS sources are most likely to be caused by excessive noise over a large range of spatial scales rather than real extended structures such as cirrus clouds.     

The radio–mid-infrared correlation and the contribution of 15-μm galaxies to the 1.4-GHz source counts

The radio counterparts to the 15-μm sources in the European Large Area ISO Survey southern fields are identified in 1.4-GHz maps down to ∼80 μJy. The radio–mid-infrared correlation is investigated and derived for the first time at these flux densities for a sample of this size. Our results show that radio and mid-infrared (MIR) luminosities correlate almost as well as radio and far-infrared (FIR), at least up to   z ≃ 0.6  . Using the derived relation and its spread together with the observed 15-μm counts, we have estimated the expected contribution of the 15-μm extragalactic populations to the radio source counts and the role of MIR starburst galaxies in the well-known 1.4-GHz source excess observed at sub-mJy levels. Our analysis demonstrates that IR emitting starburst galaxies do not contribute significantly to the 1.4-GHz counts for strong sources, but start to become a significant fraction of the radio source population at flux densities ≲0.5–0.8 mJy. They are expected to be responsible for more than 60 per cent of the observed radio counts at ≲0.05 mJy. These results are in agreement with the existing results on optical identifications of faint radio sources.     

15-GHz variability of 9C sources

We present results from a 3-yr study of the 15-GHz variability of 51 9C sources. 48 of these sources make up a subsample of a larger one complete to 25 mJy in 9C, and as the sources are selected pseudo-randomly the results should be representative of the complete sample. 29 per cent of this subsample are found to be variable above the flux calibration uncertainties of ∼6 per cent. 50 per cent of the flat-spectrum objects are variable whilst none of the steep-spectrum objects or the objects with convex spectra peaking below 5 GHz are variable. Nine of the objects studied have convex spectra and peak frequencies above 5 GHz; eight of these were found to vary at 15 GHz, suggesting that the high-frequency peaking class in this sample is largely populated by objects with jets aligned close to the line of sight whose emission is dominated by beamed components.     

The Imperial IRAS-FSC Redshift Catalogue

We present a new catalogue, the Imperial IRAS -FSC Redshift Catalogue (IIFSCz), of 60 303 galaxies selected at 60 μm from the IRAS Faint Source Catalogue (FSC). The IIFSCz consists of accurate position, optical, near-infrared and/or radio identifications, spectroscopic redshift (if available) or photometric redshift (if possible), predicted far-infrared (FIR) and submillimetre (submm) fluxes ranging from 12 to 1380 μm based upon the best-fitting infrared template. About 55 per cent of the galaxies in the IIFSCz have spectroscopic redshifts, and a further 20 per cent have photometric redshifts obtained through either the training set or the template-fitting method. For S(60) > 0.36 Jy, the 90 per cent completeness limit of the FSC, 90 per cent of the sources have either spectroscopic or photometric redshifts. Scientific applications of the IIFSCz include validation of current and forthcoming infrared and submm/mm surveys such as AKARI , Planck and Herschel , follow-up studies of rare source populations, large-scale structure and galaxy bias, local multiwavelength luminosity functions and source counts. The catalogue is publicly available at http://astro.imperial.ac.uk/~mrr/fss/ .     

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.     

Extragalactic source counts and contributions to the anisotropies of the cosmic microwave background: predictions for the Planck Surveyor mission

We present predictions for the counts of extragalactic sources, the contributions to fluctuations and their angular power spectrum in each channel foreseen for the Planck Surveyor (formerly COBRAS/SAMBA ) mission. The contribution to fluctuations owing to clustering of both radio and far-IR sources is found to be generally small in comparison with the Poisson term; however the relative importance of the clustering contribution increases and may eventually become dominant if sources are identified and subtracted down to faint flux limits. The central Planck frequency bands are expected to be 'clean': at high galactic latitude (| b | > 20°), where the reduced galactic noise does not prevent the detection of the extragalactic signal, only a tiny fraction of pixels is found to be contaminated by discrete extragalactic sources. Moreover, the 'flat' angular power spectrum of fluctuations resulting from extragalactic sources substantially differs from that of primordial fluctuations; therefore, the removal of contaminating signals is eased even at frequencies where point sources give a sizeable contribution to the foreground noise.     

X-ray source populations in the region of the open clusters NGC 6633 and IC 4756

Using deep exposures (∼10⁵ s) with the ROSAT High Resolution Imager, we have performed flux-limited surveys for X-ray sources in the vicinity of the Hyades-age open stellar clusters NGC 6633 and IC 4756, detecting 31 and 13 sources, respectively. Our primary aim is to search for so-far unrecognized cluster members. We propose identifications or classifications (cluster member, field star, extragalactic field object) for the X-ray sources, based on published membership lists, and on X-ray:optical flux ratios and optical colour–magnitude diagrams. Results of simulating the expected X-ray-emitting source populations are compared with the ROSAT measurements and with the expected capabilities of the XMM mission. The simulations provide a novel method of comparing the activity levels of NGC 6633 and IC 4756 with that of the Hyades. The measurements and simulations confirm that cluster members are the major class of X-ray emitter in these fields at flux levels     (0.1–2.4 keV), contributing ∼40 per cent of the total X-ray sources. We find six possible new members in NGC 6633 and four candidates in IC 4756; all require further observation to establish membership probability.     

High-redshift obscured quasars: radio emission at sub-kiloparsec scales

The radio properties of 11 obscured 'radio-intermediate' quasars at redshifts   z ≳ 2  have been investigated using the European Very-Long-Baseline-Interferometry Network (EVN) at 1.66 GHz. A sensitivity of  ∼25 μJy per 17 × 14 mas² beam  was achieved, and in seven out of 11 sources unresolved radio emission was securely detected. The detected radio emission of each source accounts for ∼30–100 per cent of the total source flux density. The physical extent of this emission is ≲150 pc, and the derived properties indicate that this emission originates from an active galactic nucleus (AGN). The missing flux density is difficult to account for by star formation alone, so radio components associated with jets of physical size ≳150 pc and ≲40 kpc are likely to be present in most of the sources. Amongst the observed sample steep, flat, gigahertz-peaked and compact-steep spectrum sources are all present. Hence, as well as extended and compact jets, examples of beamed jets are also inferred, suggesting that in these sources, the obscuration must be due to dust in the host galaxy, rather than the torus invoked by the unified schemes. Comparing the total to core (≲150 pc) radio luminosities of this sample with different types of AGN suggests that this sample of   z ≳ 2  radio-intermediate obscured quasars shows radio properties that are more similar to those of the high-radio-luminosity end of the low-redshift radio-quiet quasar population than those of Fanaroff–Riley type I (FR I) radio galaxies. This conclusion may reflect intrinsic differences, but could be strongly influenced by the increasing effect of inverse-Compton cooling of extended radio jets at high redshift.     

Deep 610-MHz Giant Metrewave Radio Telescope observations of the Spitzer extragalactic First Look Survey field – I. Observations, data analysis and source catalogue

Observations of the Spitzer extragalactic First Look Survey field taken at 610 MHz with the Giant Metrewave Radio Telescope are presented. Seven individual pointings were observed, covering a total area of ∼4 deg² with a resolution of 5.8 × 4.7 arcsec², PA 60°. The rms noise at the centre of the pointings is between 27 and 30 μJy before correction for the Giant Metrewave Radio Telescope (GMRT) primary beam. The techniques used for data reduction and production of a mosaicked image of the region are described, and the final mosaic, along with a catalogue of 3944 sources detected above ∼5σ, are presented. The survey complements existing radio and infrared data available for this region.     

Neural networks and the separation of cosmic microwave background and astrophysical signals in sky maps

We implement an independent component analysis (ICA) algorithm to separate signals of different origin in sky maps at several frequencies. Owing to its self-organizing capability, it works without prior assumptions on either the frequency dependence or the angular power spectrum of the various signals; rather, it learns directly from the input data how to identify the statistically independent components, on the assumption that all but, at most, one of the components have non-Gaussian distributions.
We have applied the ICA algorithm to simulated patches of the sky at the four frequencies (30, 44, 70 and 100 GHz) used by the Low Frequency Instrument of the European Space Agency's Planck satellite. Simulations include the cosmic microwave background (CMB), the synchrotron and thermal dust emissions, and extragalactic radio sources. The effects of the angular response functions of the detectors and of instrumental noise have been ignored in this first exploratory study. The ICA algorithm reconstructs the spatial distribution of each component with rms errors of about 1 per cent for the CMB, and 10 per cent for the much weaker Galactic components. Radio sources are almost completely recovered down to a flux limit corresponding to ≃0.7 σ _CMB, where σ _CMB is the rms level of the CMB fluctuations. The signal recovered has equal quality on all scales larger than the pixel size. In addition, we show that for the strongest components (CMB and radio sources) the frequency scaling is recovered with per cent precision. Thus, algorithms of the type presented here appear to be very promising tools for component separation. On the other hand, we have been dealing here with a highly idealized situation. Work to include instrumental noise, the effect of different resolving powers at different frequencies and a more complete and realistic characterization of astrophysical foregrounds is in progress.     

Streaming motions of galaxy clusters within 12 000 km s−1– II. New photometric data for the Fundamental Plane

We present new R -band photometric data for 447 galaxies, gathered for the 'Streaming Motions of Abell Clusters' (SMAC) project. The data comprise 629 individual measurements of the Fundamental Plane (FP) parameters effective radius ( R _e) and surface brightness (〈 μ 〉_e), derived from r ^1/4-law profile fitting. More than a third of the galaxies were observed more than once. The photometric precision is ∼0.02 mag as judged from comparisons of aperture photometry between repeat observations of galaxies. The combination     which enters into the Fundamental Plane relation, has internal uncertainties of ∼0.008, corresponding to < 2 per cent in estimated distance. Taken individually, the (correlated) internal errors in R _e and 〈 μ 〉_e are ∼8 per cent and ∼0.12 mag respectively. Comparisons with literature data constrain the external random errors to ≲5 per cent in distance (per observation), which is small in comparison to the ∼20 per cent scatter in the FP. The data will form part of a merged catalogue of FP parameters, presented in a companion paper.     

Final analysis of ELAIS 15-μm observations: method, reduction and catalogue

We present the final analysis of the European Large Area Infrared Space Observatory ( ISO ) Survey (ELAIS) 15-μm observations, carried out with the ISO Camera (ISOCAM) instrument on board the ISO .
The data-reduction method, known as the Lari Method, is based on a mathematical model of the behaviour of the detector and was specifically designed for the detection of faint sources in ISOCAM/ISO Photopolarimeter (ISOPHOT) data. The method is fully interactive and leads to very reliable and complete source lists.
The resulting catalogue includes 1923 sources detected with signal-to-noise ratio of  > 5  in the 0.5–100 mJy flux range and over an area of 10.85 deg² split into four fields, making it the largest non-serendipitous extragalactic source catalogue obtained to date from the ISO data.
This paper presents the concepts underlying the data-reduction method together with its latest enhancements. The data-reduction process, the production and basic properties of the resulting catalogue are discussed. The catalogue quality is assessed by means of detailed simulations, optical identifications and comparison with previous analyses.