Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://www.atnf.csiro.au/people/Ray.Norris/papers/n239.pdf
Äàòà èçìåíåíèÿ: Tue May 3 09:47:41 2011
Äàòà èíäåêñèðîâàíèÿ: Fri Feb 28 07:58:40 2014
Êîäèðîâêà:

Ïîèñêîâûå ñëîâà: m 63
Deep Spitzer observations of infrared-faint radio sources: high-redshift radio-loud AGN?
Ray P. Norris,1 Jose Afonso,2,3 Antonio Cava,4 Duncan Farrah,5 Minh T. Huynh,6 R. J. Ivison,7,8 Matt Jarvis,9 Mark Lacy,10 Minnie Mao,1,6,11,12 Claudia Maraston,13 Jean-Christophe Mauduit,6 Enno Middelberg,14 Seb Oliver,5 Nick Seymour15 and Jason Surace6 CSIRO Australia Telescope National Facility, PO Box 76, Epping, NSW, 1710, Australia, email: Ray.Norris@csiro.au 2 Observat´ orio Astron´ o de Lisboa, Faculdade de Ci^ omic encias, Universidade de Lisboa, Tapada da Ajuda, 1349-018 Lisbon, Portugal 3 Centro de Astronomia e Astrof´ a da Universidade de Lisboa, 1349-018 Lisbon, Portugal isic 4 Departamento de Astrof´ a, Facultad de CC. F´ as, Universidad Complutense de isic isic Madrid, E-28040 Madrid, Spain 5 Astronomy Centre, Dept. of Physics & Astronomy, University of Sussex, Brighton BN1 9QH, UK 6 Infrared Processing and Analysis Center, MS220-6, California Institute of Technology, Pasadena CA 91125, USA 7 UK Astronomy Technology Centre, Royal Observatory, Blackford Hil l, Edinburgh EH9 3HJ, UK 8 Institute for Astronomy, University of Edinburgh, Blackford Hil l, Edinburgh EH9 3HJ, UK 9 Centre for Astrophysics, Science & Technology Research Institute, University of Hertfordshire, Hatfield, Herts, AL10 9AB, UK 10 NRAO, 520 Edgemont Road, Charlottesvil le, VA 22903, USA 11 School of Mathematics and Physics, University of Tasmania, Private Bag 37, Hobart, 7001, Australia 12 Anglo-Australian Observatory, PO Box 296, Epping, NSW, 1710, Australia 13 Institute of Cosmology and Gravitation, Dennis Sciama Building, Burnaby Road, Portsmouth, PO1 3FX, UK 14 Astronomisches Institut, Ruhr-Universit¨ Bochum, Universit¨ at atsstr. 150, 44801 Bochum, Germany 15 Mul lard Space Science Laboratory, UCL, Holmbury St Mary, Dorking, Surrey, RH5 6NT, UK ABSTRACT Infrared-faint radio sources (IFRSs) are a rare class of ob ject which are relatively bright at radio wavelengths but very faint at infrared and optical wave1


­2­ lengths. Here we present sensitive near-infrared observations of a sample of these sources taken as part of the Spitzer Extragalactic Representative Volume Survey (SERVS). Nearly all the IFRSs are undetected at a level of 1µJy in these new deep observations, and even the detections are consistent with confusion with unrelated galaxies. A stacked image implies that the median flux density is S3.6µm 0.2 µJy or less, giving extreme values of the radio-infrared flux density ratio. Comparison of these ob jects with known classes of ob ject suggests that the ma jority are probably high-redshift radio-loud galaxies, possibly suffering from significant dust extinction. Subject headings: galaxies: formation -- galaxies: evolution -- galaxies: starburst

1.

Intro duction

Infrared-faint radio sources (IFRSs) are a rare class of ob ject which are relatively strong at radio wavelengths but very faint at infrared (IR) and optical wavelengths. They were first categorised in the Australia Telescope Large Area Survey (ATLAS: Norris et al. 2006) as radio sources with no observable IR counterpart in the co-spatial Spitzer Wide-area IR Extragalactic Survey (Lonsdale et al. 2004). Most have flux densities of a few hundred µJy at 20 cm, but some are as bright as 20 mJy. They may be related to the optically invisible radio sources found by Higdon et al. (2005, 2008), which are compact radio sources with no optical counterpart to R 25.7, although the IFRSs seem even more extreme than the Higdon ob jects. Norris et al. (2006) and Middelberg et al. (2008a) have identified 51 such sources out of 2002 radio sources in the ATLAS survey. So far, four samples of IFRS have been identified: (a) the original ATLAS/SWIRE samples identified by Norris et al. (2006) and Middelberg et al. (2008a), from which the sources in this paper are drawn, (b) the sample in the Spitzer First-Look Survey, identified by Garn & Alexander (2008), (c) the sample in ELAIS-N1 identified by Grant et al. (2011), (d) a sample in the COSMOS field (Scoville et al. 2007) identified by Zinn et al. (2011). These samples imply a sky density of of 7 per deg2 for S20 cm > 0.1 mJy. At the time of their discovery, the IFRS sources were unexpected, as SWIRE was thought to be deep enough to detect all extragalactic radio sources at z 2, regardless of whether star formation or active galactic nuclei (AGN) powered the radio emission. Possible explanations were that these sources are i) high-redshift radio-loud AGN, ii) very obscured radio galaxies at more moderate redshifts (1 < z < 2), iii) lobes of nearby but unidentified radio galaxies,


­3­ iv) very galaxies ob ject. harbour obscured, (SMGs ­ Of course examples luminous starburst galaxies, such as high-redshift submillimetre-selected Smail et al. 1997), v) high-latitude pulsars, or vi) an unknown type of it is also possible that they do not constitute a homogeneous class and of some or all of the above.

The nature of IFRSs has been hard to determine because nearly all the information on them has been obtained at radio wavelengths. Spectroscopy is difficult because the hosts are optically faint and the radio positions can also have uncertainties of the order of a few arcsec. Norris et al. (2006) stacked the positions of 22 IFRSs in the Spitzer 3.6 µm IRAC images and found no detection in the averaged image, showing that they are well below the SWIRE detection threshold. This was a surprising result at the time, as it was expected that the IFRSs represented the tail of a distribution reaching just below the SWIRE detection threshold, but it was confirmed by Garn & Alexander (2008) and by the new data presented here. Middelberg et al. (2008b) and Norris et al. (2007) targeted six IFRSs with the Australian Long Baseline Array (LBA) and detected two of the sources. The Norris et al. (2007) LBA detection constrained the source size to less than 0.03 arcsec, suggesting a compact radio core, powered by an AGN. Middelberg et al. (2008b) found the size and radio luminosity of their LBA-detected source to be consistent with a high-redshift (z > 1) compact, steepspectrum (CSS) source. The VLBI detections rule out the possibility that these particular IFRSs are simply the radio lobes of unidentified radio galaxies, or star-forming galaxies, though the initial VLBI targets were inevitably amongst the most radio-bright examples in the sample. Garn & Alexander (2008) stacked IFRS sources in the Spitzer First Look Survey at infrared wavelengths, as well as at 610 MHz. The sources they find in the FLS are very similar to the sources described here. They find that the IFRS sources can be modelled as compact Fanaroff Riley type ii (FR ii) radio galaxies at high redshift (z 4), and argue that IFRSs are predominately high-redshift radio-loud AGN. Huynh et al. (2010) used deep data from the Spitzer IRAC MUSYC Public Legacy in ECDFS (SIMPLE) pro ject and the Far-Infrared Deep Extragalactic Legacy (FIDEL) Spitzer survey to probe more deeply in the E-CDFS region, and detected two of the four IFRSs in that region. However, the two non-detections, and the faintness of the two detected sources, enabled Huynh et al. to place constraints on the sources. Their detailed modelling of their spectral energy distributions shows that they are consistent with high-redshift (z > 1) AGN. They also noted that the ratio of 20 cm flux density, S20 cm , to 3.6 µm flux density, S3.6µm , is higher than that of the general radio source population, and has significant overlap with the population of high-redshift radio galaxies (HzRG) investigated previously with Spitzer


­4­ observations (Seymour et al. 2007) and is also consistent with the radio galaxy K - z relation extended to high redshifts (Jarvis et al. 2001; De Breuck et al. 2002; Willott et al. 2003). Middelberg et al. (2011) have measured the spectral indices of the radio emission at several wavelengths of a sample of 17 strong IFRSs and find that they are significantly different from the general radio source population, and also different from the general AGN population. The spectra are steep, typically with -1, and there is a particularly prominent lack of sources with > -0.7 (where S ). Grant et al. (2011) also reported steep spectral indices, and both Middelberg et al. (2011) and Grant et al. (2011) found that several of the IFRSs were significantly polarised, suggesting an AGN rather than a star-forming galaxy. Each of these papers adds to a consensus view that the IFRSs represent a class of highredshift radio-loud galaxies. However, there is not yet unequivocal evidence for this, and it is possible that the IFRS class constitutes several different types of ob jects. Here we present near-IR data taken as part of the Spitzer Extragalactic Representative Volume Survey (SERVS) pro ject (Mauduit et al. 2011), which uses the data taken during the Spitzer warm mission to probe extended regions of the sky to a sensitivity some five times deeper than SWIRE. The SERVS pro ject includes several square degrees around the ELAIS-S1 and CDFS extended fields in which IFRSs were detected as part of the ATLAS radio survey. These deeper data provide much tighter constraints on the nature of IFRSs. Throughout this paper we assume a Hubble constant of 71 km s-1 Mpc-1 , and matter and cosmological constant density parameters of M = 0.27 and = 0.73.

2.

Observations and Analysis

The Australia Telescope Large Area Survey (ATLAS) radio survey is still in progress, covering 7 deg2 at 20 cm in two regions, surrounding the CDF-S and ELAIS-S1 fields to an r.m.s. depth of typically 10 µJy/beam. However, all the radio data used in this paper are taken from the preliminary ATLAS catalogues (Norris et al. 2006; Middelberg et al. 2008a) which surveyed the entire region to an r.m.s. depth of typically 20­30 µJy/beam. SERVS (Mauduit et al. 2011) is a medium-deep survey at 3.6 µm and 4.5 µm which exploits the warm phase mission of Spitzer, to cover an 18 deg2 field, which includes most of the ATLAS fields. For the work described in this paper, we used SERVS data in the CDFS and ELAIS-S1 fields. Only band 1 (at 3.6 µm) of the SERVS data was used because it is intrinsically more


­5­ sensitive for most galaxy spectral energy distributions than the band 2 (4.5 µm) data. Fig. 1 shows the IR and radio coverage of these fields. The SERVS data on the CDFS area covers nearly all of the area observed by SWIRE and ATLAS, but the ELAIS-S1 coverage is offset from the centre of the SWIRE/ATLAS field by about a degree. Nevertheless, a total of 39 IFRSs are common to both datasets, and are listed in Table 1. For each IFRS source detected at radio wavelengths, we took the initial positional uncertainty (i , in R.A. and Dec.) from the values listed by Norris et al. (2006) and Middelberg et al. (2008a), which were calculated from the quadrature sum of the formal fitting uncertainty and a 0.1 arcsec potential uncertainty in the position of the calibration source. We also consider the formal positional uncertainty due to noise noise = 0.3 b (Srms /Speak ), where b is the synthesised full-width-half-maximum beamwidth, Srms is the local r.m.s. flux density, and Speak is the peak flux density of the source, as discussed by Ivison et al. (2007). We further considered any positional accuracy of less than 1 arcsec to be unrealistic because of systematic offsets and intrinsic source sizes. We derive a final 3 positional uncertainty final in R.A. and Dec. as being the maximum of 3 i , 3 noise , and 1 arcsec. These positional uncertainties are listed in Table 1. For each IFRS source, we examined the SERVS 3.6 µm data for sources visible by eye as a distinct peak above the noise which fell within the error ellipse. Using this technique, we found that 3 of the 39 IFRS sources contained a 3.6 µm candidate source within the 3 radio position error ellipse. Images of two representative IFRSs are shown in Fig. 2, with their associated error ellipse. We then performed aperture photometry on each of these candidates. The flux density was measured using an aperture radius of 1.9 arcsec, which was found to be optimal in SWIRE. An aperture correction of 1.4 was applied for this radius (Surace et al. 2005). The uncertainty in the aperture flux density was estimated by measuring the r.m.s. of 100 randomly placed apertures placed near the source. The resulting flux densities are listed in Table 1. The r.m.s. noise of the SERVS data, measured in a region of sky free from visible sources, can approach 0.2 µJy. However, since SERVS data are typically confusion limited, such an r.m.s. can be misleading if used to calculate uncertainties. Instead, we placed 1000 random apertures across the image, and followed an iterative approach, rejecting apertures with a flux density 2.5 and recalculating the r.m.s. until the process converged. The r.m.s. thus obtained was 0.51 µJy in the ELAIS-S1 field and 0.64 µJy in the CDF-S field. As the SERVS data approach the confusion limit of Spitzer, it is necessary to estimate how many of the candidate cross-identifications are due to confusion. We estimated this in


­6­ two independent ways. Our most reliable estimate was obtained by shifting the IFRS radio positions by an arbitrary amount (typically 20 arcsec) which is much greater than the beamsize of either the radio or IR data, but much smaller than the scale size of variations in the image sampling. We then examined the image to estimate how many IR sources fall by chance within the error ellipse. This is a very robust way of estimating the confusion, as it builds in the varying error ellipse size, and the non-uniform sensitivity, using the real data rather than a parameterised representation. Because it uses real data, with exactly the same process in both cases, it is immune to calibration or other systematic errors. We then repeated this process ten times, shifting the positions by a different amount (<1 arcmin in all cases). A mean of 2.1 ± 1.2 of the shifted error ellipses contained a peak in the Spitzer data. Thus, of our 3 candidate identifications with unshifted data, we expect that 2.1 of these are due to confusion, leaving only zero or one genuine detections. As a check on this result, we use the source counts calculated by Barmby et al. (2008). Our radio position error ellipses have a total area of 266 arcsec2 , and we have found three candidate identifications within this area, all of which are brighter than a flux density limit of 1.1µJy. . Barmby et al. (2008) calculate the density of sources with S3.6 µm > 1.1 µJy as 205,000 deg-2 , from which we estimate that we should detect 4 sources by chance. This is greater than the number detected in the shifted data, but the difference is clearly dominated by small number statistics. Nevertheless, the Barmby result does provide a useful rough cross-check on our shifting technique. In summary, of our 39 IFRSs which lie within the SERVS fields, we find that only 3 of these have a source within the radio position error ellipse, none has a measured flux density greater than five times the formal fitting uncertainty, and our basic Monte Carlo simulations suggest that most or all of these detections are due to chance. We conclude that few or none of our sources have reliable detections and the vast ma jority of IFRSs are undetected at this level. Such faint emission from a source which is relatively strong at radio wavelengths represents an extreme condition which is not common in the local Universe. The S3.6µm distribution of IFRS can be explored to even deeper levels by stacking 3.6 µm images at the IFRS positions. In Fig. 3, we show a median stacked image obtained by summing 39 3.6 µm images extracted from the SERVS data, centered on the IFRS radio position. Because the stacking has reached the confusion limit for these data, the r.m.s. noise no longer scales as the square root of integration time, and so the r.m.s. of this stacked image is higher than would be obtained in an unconfused field, although it still offers a significant improvement over the individual images. We note that Garn & Alexander (2008)


­7­ faced a similar challenge in stacking 3.6 µm data from 8 IFRS in the First Look Survey, reaching a noise of approximately 1 µJy in the image stack. The measured r.m.s. of the median stacked image is 0.14 µJy, and there is marginal evidence for a source at the field centre whose flux density was measured to be 0.21± 0.14 µJy. Using a 3 upper limit, we conclude that the IFRSs have a median flux density of 0.63 µJy at 3.6 µm.

3.

Radio/IR prop erties of IFRSs

Since we have no redshift information for any of the IFRSs, we focus here on two derived quantities: the ratio of S20 cm and S3.6µm , or its lower limit, and the S3.6µm flux density, or its upper limit. Of the 39 IFRS which lie within the SERVS fields, we find confusion, only 1 or 2 are detected and none has a measured flux times the formal fitting uncertainty. The remaining sources are From this we deduce that the vast ma jority of IFRSs have a ratio range 200­2,000. that, after allowing for density greater than five undetected at this level. of S20 cm to S3.6µm in the

In Fig. 4 we show the ratios of S20 cm to S3.6µm as a function of redshift for a representative selection of models, and mark the limits obtained in this paper. It is clear from this figure that the only ob jects known to have such a high ratio are radio-loud AGN, such as the highredshift radio galaxies (HzRGs). In particular, we can rule out any known type of galaxy powered predominantly by star formation, such as ULIRGs, SMGs, etc., all of which fall well below the region occupied by IFRSs. If they are radio galaxies similar to those in the low-redshift Universe, then we can use their brightness to obtain constraints on redshift. In Fig. 5. we show the S3.6µm of HzRGs as a function of redshift, together with the limits obtained from our SERVS observations. The HzRGs follow a relation between redshift and S3.6µm similar to the well-known K -z relation for other radio galaxies (Willott et al. 2003). Although we caution that both these relationships are unreliable above z>3, they imply that if the IFRSs are radio galaxies, then their low 3.6 µm flux densities constrain them to lie at high redshift.


­8­ 4. What are IFRSs?

Whatever the nature of IFRSs, their properties are extreme, and not consistent with any of the well-recognised classes of ob ject. Here we review the options for what this class of ob ject is likely to be. · Star-forming galaxies. Fig. 4 shows that no known class of galaxy powered predominantly by star formation, such as ULIRGs, SMGs, etc., have ratios of 20 cm to 3.6 µm flux densities comparable with those of IFRS. A star-forming galaxy might appear in the region if it suffered from an unusually high extinction, but no star-forming galaxy is known with such high extinction. Even Arp 220 has a S20 cm /S3.6µm ratio a factor of 20 below the most moderate IFRS, and a factor of 4000 below the most extreme IFRS. Such high extinction could in principle be found in a star-forming galaxy if it were at high redshift (z > 3), where the observed 3.6 µm emission is generated in visible wavelengths in the galaxy rest-frame. However, no known star-forming galaxy generates sufficient radio power to reproduce the observed IFRS flux density. For example, Arp220 at z = 3 would have an observed flux of 5 µJy. The observed 20 cm radio luminosities of the IFRSs range from 7 â 1023 W Hz-1 (for the weakest IFRS at z = 1) to 1.4 â 1026 W Hz-1 at z=1 to 7 â 1027 W Hz-1 (for the strongest IFRS at z = 5). The most luminous star-forming galaxies, typified by SMGs at z 2 - 3, can have luminosities 1024 W Hz-1 (e.g. Ivison et al. 1998; Seymour et al. 2009), so that while the faintest of our galaxies could be caused by star-forming galaxies, most are too radio-luminous. In addition, two of the six IFRSs observed with VLBI were detected by Norris et al. (2007) and Middelberg et al. (2008b). The VLBI detections rule out star formation in these particular galaxies as the radio emission mechanism, since the synchrotron emission from star-forming galaxies rarely has sufficient brightness temperature to be detectable with VLBI (e.g. Kewley et al. 2000; Biggs et al. 2010). Furthermore, Middelberg et al. (2011) and Grant et al. (2011) show that the radio emission from several of the IFRSs studied by them is significantly polarised. This also argues against star formation, in which polarisation is generally much lower than in AGN. We conclude that the ma jority of IFRSs are unlikely to represent star-forming galaxies. · Radio lob es An early hypothesis was that IFRSs might be the extended radio lobes of an AGN whose host may be located some distance away. The VLBI detections by Norris et al. (2007) and Middelberg et al. (2008b) imply that at least about a third of


­9­ IFRSs have high brightness temperature cores, which rules out the possibility that they are radio lobes of AGN, since radio lobes do not have sufficient surface brightness to be detectable with these VLBI observations. Furthermore, most IFRS are unresolved in the high-resolution observations by Middelberg et al. (2011), making it unlikely that they are extended radio lobes. · Pulsars Cameron et al. (2011) tested the hypothesis that IFRSs may be pulsars, by performing a pulsar search on a sample of IFRSs. Their results show that any putative pulsars in the field have a pulsed flux density well below the observed flux density of the IFRS, and they conclude that the IFRS are not radio pulsars. · Radio-loud AGN The extreme values of the ratio of S20 cm to S3.6µm reported in Fig. 4 are known to occur in radio-loud AGN, and have not been observed in any other type of extragalactic ob ject. All the available observations (radio/IR flux densities, radio/IR ratio, VLBI, polarisation) are consistent with the cores of radio-loud AGN, and we conclude that radio-loud AGN are therefore a natural explanation of IFRS. If IFRSs are caused by radio-loud AGN, their remaining unusual aspect is the extreme faintness at 3.6 µm. We now consider three possible causes: · Dwarf galaxies that host a radio-loud AGN If a radio-loud AGN were to be hosted by a dwarf galaxy at moderate redshift, then the low luminosity of the host galaxy could reproduce both the faint 3.6 µm flux density and the high S20 cm to S3.6µm ratio of the IFRS. However, such ob jects are observationally unknown, and theoretically unlikely given the weak potential well of a dwarf galaxy. · Mo derate-redshift radio-loud AGN with heavy dust extinction Fig. 5 shows that known moderate-redshift (z 3) AGN would be detected in our 3.6 µm observations. We can postulate moderate-redshift ob jects which are heavily obscured at 3.6 µm, but the amounts of extinction required are very high. For example, 3C 273 at z = 1 would require an extinction of AV 50m to occupy the position of IFRSs in Fig. 5. On the other at hand, at high redshift (z > 2), where the observed 3.6 µm emission is generated in visible wavelengths in the galaxy rest-frame, as little as Av = 10 magnitudes of extinction is sufficient to raise the track of a radio-loud QSO into the regime of the IFRS. However, even if the AGN is obscured by dust, it is hard to obscure the host galaxy to the low 3.6 µm flux limits presented here. Even in Arp 220, where as much as 100 magnitudes of extinction obscure the nucleus (e.g. Haas et al. 2001), the outer shell of


­ 10 ­ the host galaxy is still bright at IR wavelengths. To reproduce the extreme IR-faintness of the IFRS, the extinction must cover the entire body of the galaxy . We therefore cannot exclude the possibility that the IFRSs could be a new class of galaxy at z 1­2, in which the entire galaxy must be shrouded in dust, reducing its integrated 3.6µm flux density by a factor of 100. · High-redshift radio-loud AGN Fig. 5 and the K -z relation show that, if IFRSs belong to a parent population similar to that of known radio-loud AGN, they must lie at z 3. We note that obscured high-redshift radio galaxies and quasars have been detected by Mart´ inez-Sansigre et al. (2006); Dey et al. (2008); Yan et al. (2007) and Sa jina et al. (2007). While IFRSs bear some similarities to these galaxies, IFRSs are even more extreme in their S20 cm to S3.6µm ratios. For example, Mart´ inez-Sansigre et al. (2006) choose their sample by requiring that S3.6µm < 45 µJy, whilst the IFRSs have S3.6µm 1 µJy. While none of the IFRSs are detected at 24 µm, the available SWIRE data at 24 µm are relatively insensitive compared to the 3.6 µm SERVS data, and, given the very low S3.6µm , only a very extreme spectral energy distribution would enable them to be detected at 24 µm. We conclude that the most natural explanation for IFRSs is that they are very similar to known classes of radio-loud AGN, but at a redshift 3. In all these properties, the IFRS most resemble the HzRGs of Seymour et al. (2007), Ivison et al. (2008) and Jarvis et al. (2009), but with even more extreme 20 cm to 3.6 µm flux density ratios. However, we cannot exclude the possibility that they are a new class of lower-redshift (1 < z < 3) radio-loud AGN in which the luminosity of the entire host galaxy is reduced by a factor of 100 by dust extinction. The 20 cm flux densities of our IFRS range from 0.14 to 26 mJy. The weakest end of the range corresponds to a luminosity 7 â 1023 W Hz-1 at z = 1 to 4 â 1025 W Hz-1 at z = 5, giving it the luminosity of an FR i galaxy at any reasonable redshift. On the other hand, the brightest end of the range corresponds to a luminosity 1.4 â 1026 W Hz-1 at z = 1 to 7 â 1027 W Hz-1 at z = 5, making it a FR ii galaxy at any reasonable redshift. Thus the class of IFRS spans both the FR i and FR ii luminosity classes, with the ma jority of galaxies straddling the FR i/FR ii break, depending on their redshift. It would be instructive to consider these results in terms of black hole mass, MB H , but Snellen et al. (2003) have shown there is no good correlation between radio luminosity and black hole mass for radio-loud galaxies and so we cannot estimate a black hole mass except in the most general terms. But we do have good observational evidence for the 3.6 µm flux densities of host galaxies of high-redshift radio galaxies, which is shown in Fig. 5.


­ 11 ­ We have detected 51 IFRS in the 7 deg2 area of ATLAS, implying a density of IFRSs on the sky of 7 per deg2 for S20 cm > 0.1 mJy. It is difficult to make a meaningful comparison of the numbers of IFRSs with source counts at the present time without a better constraint on redshift. The SKADS simulation (Wilman et al. 2008) gives a sky density of 0.5 FR ii galaxies deg-2 at z 4, implying that only 3­4 of the IFRSs discussed here are FR ii galaxies at z 4. However, the space density of all but the highest-luminosity, FRII-type radio sources are very poorly constrained at high redshifts, and AGN models are very poorly constrained by current observations at the low flux densities probed here. Thus it is possible (Zinn et al. 2011) that the IFRS could be radio-loud AGN of moderate radio luminosity and much higher space density than their more powerful counterparts.

5.

Conclusion

We have searched for IR counterparts to 39 IR-faint radio sources using deep 3.6 µm observations from SERVS. Even though the sensitivity is 3­5 times better than the previous observations (with SWIRE), few or none of the IFRSs are detected after taking into account the number of chance associations, and a stacked image indicates a median 3.6 µm flux density of 0.21 ±0.14µJy. This places extreme constraints on the properties of these sources, making it likely that they are radio-loud AGN at redshifts z 3, or heavily dust-obscured radio-loud AGN at redshifts z 1. While some may have radio-to-IR ratios similar to 3C 273, but at a much higher redshift, the most extreme of them require several magnitudes of obscuration in the optical/NIR to remain undetected by deep imaging. While we cannot rule out the possibility that more than one type of ob ject may be represented by IFRSs, the evidence suggests that a significant proportion, if not all, of the IFRS sources are either · radio-loud AGN (similar to known high-redshift radio-loud AGN) at z 3, or

· a new class of lower-redshift (1 < z < 3) radio-loud AGN in which the luminosity of the entire host galaxy is severely reduced by dust extinction.

Acknowledgements This work is based in part on observations made with Spitzer, which is operated by the Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA. Support for this work was provided by NASA through an award issued by JPL/Caltech. This


­ 12 ­ research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. The Australia Telescope is funded by the Commonwealth of Australia for operation as a National Facility managed by CSIRO. JA gratefully acknowledges the support from the Science and Te