A&A 391, 509-517 (2002)
DOI: 10.1051/0004-6361:20020869
I. M. Hook1,
-
R. G. McMahon2
-
P. A. Shaver3
-
I. A. G. Snellen4
1 - Department of Physics, University of Oxford,
Nuclear & Astrophysics Laboratory,
Keble Road,
Oxford OX1 3RH, UK
2 -
Institute of Astronomy, Madingley Road, Cambridge CB3 0HA, UK
3 -
European Southern Observatory, Karl Schwarzschild Straße 2,
85748 Garching b. München, Germany
4 -
Institute for Astronomy, Royal Observatory, Edinburgh, EH9 3HJ, UK
Received 5 February 2002 / Accepted 11 June 2002
Abstract
We present results of a new, large survey for high-redshift radio-loud
quasars, which targets quasars with z>4. The survey is based on the
PMN and NVSS radio surveys, optically identified using digitised UKST
B, R and I plates. Six new z>4 flat-spectrum QSOs have been
discovered, and one previously known z>4 QSO rediscovered, based on
their red optical colours. The QSOs discovered in this survey are
bright in both radio and optical bands; in particular PMN J1451-1512
(z=4.763,
,
R=19.1) and PMN J0324-2918 (z=4.630, R=18.7) are very luminous. PMN J1451-1512 at z=4.763 is also now
the most distant radio-selected quasar. In addition, 9 new quasars
with 3.5<z<4.0 were discovered during the survey. We present spectra
and finding charts for the new quasars. We also derive a surface
density of
for
flat-spectrum z>4 QSOs with
mJy and
mag.
Key words: quasars: general
Radio selection remains one of the most efficient ways of finding high-redshift AGN. This approach has the further advantage of being less prone to selection effects than optical selection, since radio emission is unaffected by either intrinsic or extrinsic absorption due to dust. The specific aim of this work was to find optically bright, radio-selected high-redshift quasars. These can be used for unbiased studies of damped Lyman alpha systems at high redshift and other follow-up studies such as searches for associated high-redshift galaxy clusters.
We therefore began to carry out a large, systematic survey aimed specifically at z>4 QSOs. Our method involves the optical identification of flat-spectrum radio sources and the spectroscopic follow-up of the red stellar identifications. This approach exploits the fact that quasars at high redshift have redder optical colours than their low-redshift counterparts due to absorption by intervening HI (see Fig. 1 in Hook et al. 1995), and has proved successful at finding high-redshift quasars in the past (Hook et al. 1995, 1996, 1998).
Previous work using well-defined quasar samples has shown that z>4radio-loud quasars are likely to be rare objects, both because the
quasar population as a whole appears to decline at redshifts above
2-3 (Kennefick et al. 1996; Hawkins & Veron 1996; Schmidt et al. 1995; Warren et al. 1995) and because
radio-loud quasars represent only about 10% of the full quasar
population.
Specific studies of the radio-loud quasar population have shown that
these objects are indeed rare at z>4. Dunlop & Peacock (1990)
presented strong evidence for a drop in the space density of
radio-loud quasars between z=2 and
based on
radio-selected samples reaching
mJy. More
recently, significant progress has been made towards understanding the
evolution of the radio-loud quasar population out to
and the
potential effects of absorption by dust, by the study of a completely
identified, large area, flat-spectrum radio sample with
mJy (Shaver et al. 1996; Wall et al., in
preparation). The low numbers of high-redshift quasars found in these
studies demonstrates that there is a distinct drop-off in the space
density of quasars at z>3. Thus for our new survey to be successful
it must reach fainter radio flux density limits than the above surveys
(to sample further down the luminosity function), and cover a
significant fraction of the sky.
Here we present the first results of this new survey for high-redshift
radio-loud quasars. As will be seen in Sect. 2, the survey uses
deeper radio and optical data than previous radio-loud quasar surveys,
and covers a very large area in the Southern sky (
7500 sq deg,
comparable to that of the planned 10 000 sq deg of the Sloan survey).
Our survey has produced seven z>4 flat-spectrum quasars, one of
which was previously known. The survey complements the survey of
Snellen et al. (2001) which contains four flat-spectrum z>4 in the
Northern sky, selected using a a similar method.
The parent radio sample used in this study is based on the
Parkes-MIT-NRAO radio survey (PMN, Griffith et al. 1995 and references
therin), selected at 5 GHz. The data cover the southern sky with
to a flux density limit of 20-72 mJy depending on
declination. No additional flux density limit was applied when
carrying out the survey, but note that when considering the statistics
of the final quasar sample, we consider sources with
mJy,
since the PMN completeness at lower flux density levels is patchy (see
Fig. 1).
To provide accurate positions and spectral index information, the PMN
catalogue was matched to the 1.4 GHz NRAO-VLA Sky Survey (NVSS, Condon
et al. 1998), which covers the declination range
.
All NVSS sources matching within a 2 arcmin radius of a PMN
position were kept.
Since the beam sizes of the PMN and NVSS surveys are different (4.2'for PMN and 45'' for NVSS), the 1.4 flux density used to determine
spectral indices were calculated by summing the total flux density
from all NVSS sources whose positions were within a 2 arcmin radius
of the PMN position. Flat spectrum objects with
mJy were then selected using the criterion
where
.
This gave 5976 PMN radio sources.
The NVSS positions were matched to optical catalogs which were derived from scans of UK Schmidt Telescope (UKST) plates. The scans were produced using the Automatic Plate Measuring (APM) facility (McMahon & Irwin 1992). The optical identification procedure was similar to that used in the past to make POSS-based identifications of radio samples, described in Hook et al. (1995, 1996, 1998). However the current survey also makes use of I-band plates for the first time.
The plate data used is from the UKST survey of the southern sky in the
B, R and I bands, reaching limiting magnitudes of approximately 22.6,
21.0 and 19.5 mag respectively. The blue plates are in the passband,
Å, and have 606 centres spaced at
intervals of 5 degrees. In most cases the matching R plate used was
from the OR survey with a passband of
,
although in
some fields not yet covered by the OR survey, the old R plate was
used (
). The I plates used IV-N emulsion + RG715
filter giving a passband of
.
The maximum redshift of quasars that our survey can detect in
principle is defined by the redshift at which the
line is
redshifted out of the redder passbands.
For the R-plates this occurs at redshifts above z=4.7 and for the I plates
at redshifts above z=6.3.
Figure 1 shows the spatial distribution of the radio sources which form the basis of the survey. These lie in the region of overlap of the PMN and NVSS surveys. The effective area of our quasar survey is defined by the subset of this region that is covered by APM scans of UKST plates.
All the available B, R and I plates have now been scanned in the
region of overlap with our radio sample (
,
avoiding the galactic plane,
). At
the time of the spectroscopic observations reported here, an area of
7525 sq deg was covered by plates in the B and R bands (4637.5 sq
deg in the South Galactic Cap region and 2887.5 sq deg in the North
Galactic Cap region). A total of 3887.5 sq deg also had I-band
plates. The spatial distribution of the plates is shown in
Fig. 2. Of this, two regions around
h,
have no radio data, as shown in Fig. 1.
The total area of overlap with the B, R data is 7265.5 sq deg of which
4637.5 sq deg is in the South Galactic Cap (SGC) region and 2628.0 sq deg
is in the North Galactic Cap (NGC). Of this a total of 3637.5 sq deg also had I data.
In addition there are small areas within the region of radio/optical
overlap that were not yet covered by the NVSS survey at the time our
QSO survey was carried out. This resulted in some PMN sources not
being matched with an NVSS counterpart in certain regions, as can be
seen in Fig. 2. When considering the statistics of our
final quasar sample, we take this incompleteness into account in a
statistical way by calculating the fraction of PMN sources with
mJy
that do not have NVSS counterparts. In the SGC the
matched fraction is 0.955 and in the NGC it is 0.871.
Optical identifications were made based on positional coincidence of the NVSS position with an optical counterpart on the UKST plates. A matching radius of 3.0'' between the NVSS and optical positions was used (see Fig. 3 for the distribution of positional differences).
From these identifications, red, stellar objects were selected for
spectroscopic follow up. Figure 4 shows colour-magnitude
and colour-colour diagrams for optical identifications within 3'' of
the NVSS position. The first spectroscopic sample contained any
stellar object with
.
A second sample, the "I-band
sample'', was then defined, which contained stellar objects with
,
or
that were not already included in
the first sample. The criterion for being considered a stellar object
was
where
is the APM
classification parameter measured from the R plate. If the object was
not detected in R then the I plate classification was used. These
selection criteria resulted in a sample of 228 sources. Of these, 33
had known redshifts from the literature prior to the start of this
project.
In total, 169 candidates were observed
spectroscopically. During two runs in October 1998 and May 1999, 157
objects were observed at the ESO 3.6 m, La Silla, Chile using the EFOSC
spectrograph. For both runs the detector used was "ccd40'', a Loral CCD
with
pixels, binned
.
For the 1999 May
observations grism #12 was used, covering the spectral range
5800 Å-9500 Å with a dispersion of 4.23 Å per binned pixel. For
the 1998 October observations grisms #4 and #5 were used covering
the ranges 3420-7460 Å and 6000-10 000 Å respectively with dispersions of
4.08 and 4.27 Å per binned pixel. The slit width was chosen to match
the seeing, typically 1.2-1.5 arcsec. For the 1998 October
observations a sky PA of
was used, and for the 1999 May
observations the slit was at the paralactic angle. Approximate
relative flux calibration was obtained using observations of
spectrophotometric standards from Hamuy et al. (1994). Twelve sources
were observed at AAT using the RGO spectrograph in October 1996.
26 sources remain to be observed, of which 24 are in the North Galactic Cap (NGC) and 2 are in the South Galactic Cap (SGC) region.
Number of candidates | Number confirmed | |||
with z data | 3.5<z<4.0 | 4.0<z<4.5 | z>4.5 | |
Present observations | 169 | 7 | 4 | 2 |
Literature data | 33 | 2 | 1 | 0 |
Total | 202 | 9 | 5 | 2 |
The data were reduced using standard procedures with the
IRAF software environment.
Table 1 summarises the current survey status and in Table 2 we give redshifts, optical magnitudes and radio flux densities for the objects which were identified as new high-redshift quasars. The spectra are shown in Fig. 5, and finding charts, with J2000 coordinates, are given in Fig. 7. Two of the z>4 QSOs are discussed in more detail below.
this is now the most distant radio-selected quasar. The object is
blank on the UKST B-band plate, detected in R, and very bright in I(17.3). Near-IR magnitudes from the 2MASS survey are
,
,
.
The optical spectrum (Fig. 5) shows strong lines of
SiIV/OIV and CIV from which we estimate a redshift of
4.763. However the absorption of the CIV and Ly-
line make
redshift determination uncertain. A 4800 s IR spectrum was obtained
using SOFI at the 3.5 m NTT, La Silla on 31 May 1999. The spectrum
shows strong, broad emission lines of CIII and MgII at a redshift of
4.764 and approximately 4.78 respectively (the MgII line has a complex
shape, see Fig. 6). We adopt a redshift of 4.763 for this
QSO.
The NVSS map (see Fig. 8) shows a source at the optical
position of the QSO, and a second NVSS source about 4 arcmin away. The
PMN position appears to lie between the two, although much closer to the
first position (that of the QSO). It is possible that the PMN flux
density is overestimated because of a contribution from the second
source.
This source has been confirmed as an X-ray source with ROSAT HRI and its X-ray properties are the subject of another paper (Fabian et al. 2001).
Name | RA | DEC | z | R | B | I |
![]() |
![]() |
Place & Date | comments |
J2000 | discovered | |||||||||
PMN J0022-0759 | 00 22 00.246 | -07 59 16.03 | 3.896 | 19.47 | 21.26 | >19.25 | 54 | -0.18 | ESO 10/98 | |
PMN J0214-0518 | 02 14 29.295 | -05 17 44.55 | 3.986 | 18.42 | 20.54 | 18.75 | 93 | 0.61 | ESO 11/00 | |
PMN J0235-1805 | 02 34 55.143 | -18 06 08.49 | 4.314 | 18.79 | 22.00 | - | 49 | -0.16 | ESO 11/00 |
![]() |
PMN J0324-2918 | 03 24 44.280 | -29 18 21.10 | 4.630 | 18.66 | 22.18 | - | 354 | 0.30 | ESO 11/00 | |
PMN J0326-3253 | 03 27 00.407 | -32 54 18.65 | 3.463 | 19.40 | 21.15 | - | 78 | -0.01 | ESO 11/00 | |
PMN J0525-3343 | 05 25 06.166 | -33 43 05.34 | 4.413 | 18.50 | 21.47 | - | 210 | 0.06 | AAT 10/96 | X-ray source, Fabian et al. (2001) |
PMN J1043-2140 | 10 43 11.963 | -21 40 47.97 | 3.774 | 20.66 | >22.39 | >18.99 | 59 | -0.07 | ESO 05/99 | |
PMN J1108-1804 | 11 08 48.035 | -18 04 50.85 | 3.433 | 20.00 | 21.48 | 19.18 | 48 | -0.31 | ESO 05/99 |
![]() |
PMN J1429-1616 | 14 29 31.410 | -16 15 40.44 | 3.842 | 19.76 | 21.26 | 19.42 | 50 | 0.34 | ESO 05/99 | |
PMN J1451-1512 | 14 51 47.052 | -15 12 19.99 | 4.763 | 19.14 | >22.60 | 17.25 | 90 | 0.89 | ESO 05/99 |
![]() |
PMN J2042-2223 | 20 42 57.278 | -22 23 26.69 | 3.630 | 19.70 | 21.52 | - | 184 | 0.23 | ESO 10/98 | |
PMN J2134-0419 | 21 34 12.006 | -04 19 09.87 | 4.346 | 19.98 | >22.80 | 19.73 | 221 | -0.23 | ESO 10/98 | |
PMN J2219-2719 | 22 19 35.304 | -27 19 02.76 | 3.634 | 19.72 | 21.61 | - | 221 | -0.27 | ESO 11/00 | |
PMN J2220-3336 | 22 20 26.957 | -33 36 59.44 | 3.691 | 21.08 | 22.73 | - | 123 | 0.67 | ESO 11/00 |
![]() |
PMN J2314+0201 | 23 14 48.722 | +02 01 50.86 | 4.110 | 19.87 | >22.18 | >19.28 | 97 | -0.22 | ESO 10/98 | X-ray source, Boller et al. (1997) |
Previously known QSOs contained in the sample | Reference | |||||||||
PMN J1028-0844 | 4.276 | 18.84 | 20.72 | 18.37 | 159 | -0.43 | Zickgraf et al. (1997), X-ray source 0.3' away | |||
PMN J1230-1139 | 3.528 | 19.45 | 21.95 | - | 374 | -0.19 | Drinkwater et al. (1997) | |||
PMN J2003-3251 | 3.78 | 17.60 | 19.31 | - | 1248 | 0.77 | Peterson et al. (1982) |
Name | Ly![]() |
SiIV/OIV | CIV | CIII | mean z |
PMN J0022-0759 | 5956.46 | 6831.46* | 7583.92* | - | 3.896 |
PMN J0214-0518 | 6053.03 | 6978.49 | 7714.46 | - | 3.986 |
PMN J0235-1805 | 6455.15 | - | 8225.35 | - | 4.314 |
PMN J0324-2918 | 6846.70 | - | 8703.26 | - | 4.630 |
PMN J0326-3253 | 5422.49 | 6234.08 | 6913.65 | 8506.93 | 3.463 |
PMN J0525-3343 | 6627.50 | 7536.03 | 8347.96 | - | 4.413![]() |
PMN J1043-2140 | - | - | 7388.09 | 9101.68 | 3.774 |
PMN J1108-1804 | - | - | 6863.86 | 8446.82 | 3.433 |
PMN J1429-1616 | - | 6770.20 | 7498.27 | 9229.22 | 3.842 |
PMN J1451-1512 | 7004.19 | 8057.75 | 8889.01* | 10989.2
![]() |
4.763 |
PMN J2042-2223 | 5625.07 | - | 7166.39 | - | 3.630 |
PMN J2134-0419 | 6488.28 | 7488.82 | 8269.41 | - | 4.346 |
PMN J2219-2719 | 5628.22 | 6482.09 | 7138.42* | - | 3.634 |
PMN J2220-3336 | 5684.32 | 6570.76 | 7265.08 | - | 3.691 |
PMN J2314+0201 | 6195.83 | 7153.88 | 7914.37 | - | 4.110 |
![]() ![]() |
By using radio and optical multicolour data covering a significant fraction of the sky, we have produced a sample of high-redshift, optically bright radio-loud quasars.
Since this new quasar sample is well-defined, we can use it to
estimate the surface density of z>4 quasars. To do this we consider
objects with radio flux densities above the brightest limit of the PMN
surveys, 72 mJy. The SGC region of our survey has a high completeness
of spectroscopic follow up (only two objects, or 1%, were not
observed in this region) and also had a high completeness of the NVSS
survey (0.955). There are four z>4 QSOs in the SGC region of the
complete sample, which implies a surface density of
for z>4 QSOs with
mJy.
If the whole survey is considered, and the completeness for the SGC
and NGC are taken as 94% and 61% respectively (as implied by the
number of spectroscopically observed candidates in each region
combined with the NVSS completeness factors derived in Sect. 2) then
the derived surface density is
.
This is fully consistent with the value determined from the
SGC alone.
This value is similar to that of 1 per 1600 sq deg (
)
found by Snellen et al. (2001) using a similar technique
although with a slightly different radio and optical flux density
limits (R= 20 compared to R=21 for the current survey,
mJy compared to S>72 mJy, an upper redshift limit defined by
the red optical filter of
rather than
4.7, and a
radio spectral index cut at -0.35 rather than -0.5).
Finally, the new sample of quasars presented in this paper represent some of the most luminous objects in the Universe and may also represent extreme peaks in the matter density distribution at high-redshift. They are therefore ideal targets for various follow-up programs such as high-redshift absorption line studies and searches for associated high-redshift clusters.
Acknowledgements
RGM thanks the Royal Society for support. We thank Jason Spyromilio for providing the IR spectrum of PMN J1451-1512.
We also thank Mike Irwin and the staff at the APM facility in Cambridge for producing the scans of UKST plates used in this survey.
This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center, funded by the National Aeronautics and Space Administration and the National Science Foundation.