CENSORS: A Combined EIS–NVSS Survey Of Radio Sources – I. Sample definition, radio data and optical identifications

A new sample of radio sources, with the designated name CENSORS (A Combined EIS–NVSS Survey Of Radio Sources), has been defined by combining the National Radio Astronomy Observatory Very Large Array Sky Survey (NVSS) at 1.4 GHz with the ESO Imaging Survey (EIS) Patch D, a 3° by 2° region of sky centred at RA     , Dec. −21°00'00' (J2000). New radio observations of 199 NVSS radio sources with NVSS flux densities   S _{1.4 GHz} > 7.8 mJy  are presented, and are compared with the EIS I -band imaging observations which reach a depth of   I ∼ 23  ; optical identifications are obtained for over two-thirds of the ∼150 confirmed radio sources within the EIS field. The radio sources have a median linear size of 6 arcsec, consistent with the trend for lower flux density radio sources to be less extended. Other radio source properties, such as the lobe flux density ratios, are consistent with those of brighter radio source samples. From the optical information, 30–40 per cent of the sources are expected to lie at redshifts   z ≳ 1.5  .
One of the key goals of this survey is to accurately determine the high-redshift evolution of the radio luminosity function. These radio sources are at the ideal flux density level to achieve this goal; at redshifts   z ∼ 2  they have luminosities which are around the break of the luminosity function and so provide a much more accurate census of the radio source population at those redshifts than the existing studies of extreme, high radio power sources. Other survey goals include investigating the dual-population unification schemes for radio sources, studying the radio luminosity dependence of the evolution of radio source environments, and understanding the radio power dependence of the K – z relation for radio galaxies.     

The IPHAS-POSS-I proper motion survey of the Galactic plane

We present millimetre observations of a sample of 12 high-redshift ultraluminous infrared galaxies (ULIRGs) in the extended growth strip (EGS). These objects were initially selected on the basis of their observed mid-IR colours (  0.0 < [3.6]−[4.5] < 0.4  and  −0.7 < [3.6]−[8.0] < 0.5  ) to lie at high redshift  1.5 ≲ z ≲ 3  , and subsequent 20–38 μm mid-IR spectroscopy confirms that they lie in a narrow redshift window centred on   z ≈ 2  . We detect 9/12 of the objects in our sample at high significance  (>3σ)  with a mean 1200 μm flux of  〈 F _1200 μm〉= 1.6 ± 0.1  mJy. Our millimetre photometry, combined with existing far-IR photometry from the Far-IR Deep Extragalactic Legacy Survey (FIDELS) and accurate spectroscopic redshifts, places constraints both sides of the thermal dust peak. This allows us to estimate the dust properties, including the far-IR luminosity, dust temperature and dust mass. We find that our sample is similar to other high- z and intermediate- z ULIRGs, and local systems, but has a different dust selection function than submillimeter-selected galaxies. Finally, we use existing 20-cm radio continuum imaging to test the far-IR/radio correlation at high redshift. We find that our sample is consistent with the local relation, implying little evolution. Furthermore, this suggests that our sample selection method is efficient at identifying ultraluminous, starburst-dominated systems within a very narrow redshift range centred at   z ∼ 2  .     

The redshift distribution of FIRST radio sources at 1 mJy

We present spectra for a sample of radio sources from the FIRST survey, and use them to define the form of the redshift distribution of radio sources at mJy levels. We targeted 365 sources and obtained 46 redshifts (13 per cent of the sample). We find that our sample is complete in redshift measurement to R ∼18.6, corresponding to z ∼0.2. Galaxies were assigned spectral types based on emission-line strengths. Early-type galaxies represent the largest subset (45 per cent) of the sample and have redshifts 0.15≲ z ≲0.5; late-type galaxies make up 15 per cent of the sample and have redshifts 0.05≲ z ≲0.2; starbursting galaxies are a small fraction (∼6 per cent), and are very nearby ( z ≲0.05). Some 9 per cent of the population have Seyfert 1/quasar-type spectra, all at z ≳0.8, and 4 per cent are Seyfert 2 type galaxies at intermediate redshifts ( z ∼0.2).
Using our measurements and data from the Phoenix survey (Hopkins et al.), we obtain an estimate for N ( z ) at S _1.4 GHz≥1 mJy and compare this with model predictions. At variance with previous conclusions, we find that the population of starbursting objects makes up ≲5 per cent of the radio population at S ∼1 mJy.     

A sample of mJy radio sources at 1.4 GHz in the Lynx and Hercules fields – I. Radio imaging, multicolour photometry and spectroscopy

With the goal of identifying high-redshift radio galaxies with Fanaroff–Riley class I (FR I) classification, here are presented high-resolution, wide-field radio observations, near-infrared and optical imaging and multi-object spectroscopy of two fields of the Leiden–Berkeley Deep Survey. These fields, Hercules.1 and Lynx.2, contain a complete sample of 81 radio sources with   S _1.4 GHz > 0.5 mJy  within 0.6 deg². This sample will form the basis for a study of the population and cosmic evolution of high-redshift, low-power, FR I radio sources which will be presented in Paper II. Currently, the host galaxy identification fraction is 86 per cent with 11 sources remaining unidentified at a level of   r '≥ 25.2 mag  (Hercules; four sources) or   r '≥ 24.4 mag  (Lynx; seven sources) or   K ≳ 20 mag  . Spectroscopic redshifts have been determined for 49 per cent of the sample and photometric redshift estimates are presented for the remainder of the sample.     

ZEN2: a narrow J-band search for z∼ 9 Lyα emitting galaxies directed towards three lensing clusters

We present the results of a continuing survey to detect Lyα emitting galaxies at redshifts   z ∼ 9  : the ' z equals nine' (ZEN) survey. We have obtained deep VLT Infrared Spectrometer and Array Camera observations in the narrow J -band filter NB119 directed towards three massive lensing clusters: Abell clusters 1689, 1835 and 114. The foreground clusters provide a magnified view of the distant Universe and permit a sensitive test for the presence of very high redshift galaxies. We search for   z ∼ 9 Lyα  emitting galaxies displaying a significant narrow-band excess relative to accompanying J -band observations that remain undetected in Hubble Space Telescope ( HST )/Advanced Camera for Surveys (ACS) optical images of each field. No sources consistent with this criterion are detected above the unlensed 90 per cent point-source flux limit of the narrow-band image,   F _NB= 3.7 × 10⁻¹⁸ erg s⁻¹ cm⁻²  . To date, the total coverage of the ZEN survey has sampled a volume at   z ∼ 9  of approximately 1700 comoving Mpc³ to a Lyα emission luminosity of  10⁴³ erg s⁻¹  . We conclude by considering the prospects for detecting   z ∼ 9 Lyα  emitting galaxies in light of both observed galaxy properties at   z < 7  and simulated populations at   z > 7  .     

The evolution of star formation in quasar host galaxies

We have used far-infrared data from IRAS , Infrared Space Observatory ( ISO ), Spitzer Wide-Area Infrared Extragalactic (SWIRE), Submillimetre Common User Bolometer Array (SCUBA) and Max-Planck Millimetre Bolometer (MAMBO) to constrain statistically the mean far-infrared luminosities of quasars. Our quasar compilation at redshifts  0 < z < 6.5  and I -band luminosities  −20 < I _AB < −32  is the first to distinguish evolution from quasar luminosity dependence in such a study. We carefully cross-calibrate IRAS against Spitzer and ISO , finding evidence that IRAS 100-μm fluxes at <1 Jy are overestimated by ∼30 per cent. We find evidence for a correlation between star formation in quasar hosts and the quasar optical luminosities, varying as star formation rate (SFR)  ∝ L ^0.44±0.07_opt  at any fixed redshift below   z = 2  . We also find evidence for evolution of the mean SFR in quasar host galaxies, scaling as  (1 + z )^1.6±0.3  at   z < 2  for any fixed quasar I -band absolute magnitude fainter than −28. We find no evidence for any correlation between SFR and black hole mass at  0.5 < z < 4  . Our data are consistent with feedback from black hole accretion regulating stellar mass assembly at all redshifts.     

The link between SCUBA and Spitzer: cold galaxies at z≲ 1

We show that the far-IR properties of distant Luminous and UltraLuminous InfraRed Galaxies (LIRGs and ULIRGs, respectively) are on average divergent from analogous sources in the local Universe. Our analysis is based on Spitzer Multiband Imaging Photometer (MIPS) and Infrared Array Camera (IRAC) data of   L _IR > 10¹⁰ L_⊙, 70 μm  selected objects in the  0.1 < z < 2  redshift range and supported by a comparison with the IRAS Bright Galaxy Sample. The majority of the objects in our sample are described by spectral energy distributions (SEDs) which peak at longer wavelengths than local sources of equivalent total infrared luminosity. This shift in SED peak wavelength implies a noticeable change in the dust and/or star-forming properties from   z ∼ 0  to the early Universe, tending towards lower dust temperatures, indicative of strong evolution in the cold dust, 'cirrus', component. We show that these objects are potentially the missing link between the well-studied local IR-luminous galaxies, Spitzer IR populations and SCUBA sources – the   z < 1  counterparts of the cold   z > 1  SubMillimetre Galaxies (SMGs) discovered in blank-field submillimetre surveys. The Herschel Space Observatory is well placed to fully characterize the nature of these objects, as its coverage extends over a major part of the far-IR/sub-mm SED for a wide redshift range.     

Deviations from passive evolution – star formation and the ultraviolet excess in z∼ 1 radio galaxies

Galaxy colours are determined for two samples of 6C and 3CR radio sources at z ∼ 1, differing by a factor of ∼6 in radio power. Corrections are made for emission-line contamination and the presence of any nuclear point source, and the data analysed as a function of both redshift and the radio source properties. The galaxy colours are remarkably similar for the two populations, and the ultraviolet excess evolves with radio source size similarly in both samples, despite the fact that the alignment effect is more extensive for the more powerful 3CR radio galaxies. These results seem to suggest that the alignment effect at these redshifts does not scale strongly with radio power, and is instead more closely dependent on galaxy mass (which is statistically comparable for the two samples). However, it is likely that the presence of relatively young (≲several times 10⁸-yr-old) stellar populations has considerably contaminated the K -band flux of these systems, particularly in the case of the more powerful 3CR sources, which are ∼0.5 mag more luminous than the predictions of passive evolution models at z ∼ 1. The higher luminosity of the 3CR alignment effect is balanced by emission at longer wavelengths, thereby leading to comparable colours for the two samples.     

Further constraints on the evolution of Ks-selected galaxies in the GOODS/CDFS field

We have selected and analysed the properties of a sample of  2905 K_s < 21.5  galaxies in  ∼131 arcmin²  of the Great Observatories Origins Deep Survey (GOODS) Chandra Deep Field South (CDFS), to obtain further constraints on the evolution of K_s -selected galaxies with respect to the results already obtained in previous studies. We made use of the public deep multiwavelength imaging from the optical B through the infrared (IR) 4.5-μm bands, in conjunction with available spectroscopic and COMBO17 data in the CDFS, to construct an optimized redshift catalogue for our galaxy sample. We computed the K_s -band luminosity function and determined that its characteristic magnitude has a substantial brightening and a decreasing total density from   z = 0  to  〈 z 〉= 2.5  . We also analysed the colours and number density evolution of galaxies with different stellar masses. Within our sample, and in contrast to what is observed for less massive systems, the vast majority (∼85–90 per cent) of the most massive  ( M > 2.5 × 10¹¹ M_⊙)  local galaxies appear to be in place before redshift   z ∼ 1  . Around 65–70 per cent of the total assemble between redshifts   z = 1  and 3 and most of them display extremely red colours, suggesting that plausible star formation in these very massive systems should mainly proceed in obscured, short-time-scale bursts. The remaining fraction (up to ∼20 per cent) could be in place at even higher redshifts   z = 3–4  , pushing the first epoch of formation of massive galaxies beyond the limits of current near-IR surveys.     

A complete sample of radio sources in the North Ecliptic Cap, selected at 38 MHz – III. Further imaging observations and the photometric properties of the sample

Further imaging observations of a sample of radio sources in the North Ecliptic Cap are presented and a number of new identifications are made. Using redshifts from spectroscopic data presented in a companion paper by Lacy et al., the photometric properties of the galaxies in the sample are discussed. It is shown that: (1) out to at least z ≈0.6 radio galaxies are good standard candles irrespective of radio luminosity; (2) for 0.6≲ z ≲1 a large fraction of the sample has magnitudes and colours consistent with a non-evolving giant elliptical, and (3) at higher redshifts, where the R -band samples the rest-frame UV flux, most objects have less UV luminosity than expected if they form their stellar populations at a constant rate from a high redshift to z ∼1 in unobscured star-forming regions (assuming an Einstein–de Sitter cosmology). The consequences of these observations are briefly discussed.     

A sample of 6C radio sources designed to find objects at redshift z>4– III. Imaging and the radio galaxy K–z relation

In this paper, the third and final of a series, we present complete K -band imaging and some complementary I -band imaging of the filtered 6C* sample. We find no systematic differences between the K – z relation of 6C* radio galaxies and those from complete samples, so the near-infrared properties of luminous radio galaxies are not obviously biased by the additional 6C* radio selection criteria (steep spectral index and small angular size). The 6C* K – z data significantly improve delineation of the K – z relation for radio galaxies at high redshift ( z >2) . Accounting for non-stellar contamination, and for correlations between radio luminosity and stellar mass, we find little support for previous claims that the underlying scatter in the stellar luminosity of radio galaxies increases significantly at z >2 . In a particular spatially flat universe with a cosmological constant (Ω_M=0.3 and Ω_Λ=0.7) , the most luminous radio sources appear to be associated with galaxies with a luminosity distribution with a high mean (≈5  L *), and a low dispersion ( σ ∼0.5 mag) which formed their stars at epochs corresponding to z ≳2.5 . This result is in line with recent submillimetre studies of high-redshift radio galaxies and the inferred ages of extremely red objects from faint radio samples.     

Constraints on the nuclear emission of the Circinus galaxy: the torus

In the context of the unified model of Seyfert galaxies, we use observations from the literature and a radiative transfer model to investigate the near-IR to mm emission produced by the presumed torus in the Circinus galaxy, from 2 μm to 1.3 mm. From the infrared SED modelling, we find that the total luminosity ( L _IR) in this wavelength range consists of similar contributions from the torus and starburst with a ratio of nuclear luminosity to starburst luminosity ( L _NUC/ L _SB)∼0.8 .
By using a similar torus model to that of NGC 1068, but without the conical dust , we find an upper limit to the outer torus radius of ∼12 pc with a best fit of ∼2 pc. The upper limit torus size estimated from the radiative transfer modelling is consistent with the 16-pc torus radius estimated from near-IR imaging polarimetry of Circinus.     

The Deep X-Ray Radio Blazar Survey (DXRBS) – II. New identifications

We have searched the archived, pointed ROSAT Position Sensitive Proportional Counter data for blazars by correlating the WGACAT X-ray data base with several publicly available radio catalogues, restricting our candidate list to serendipitous X-ray sources with a flat radio spectrum ( α _r≤0.70, where S _ν ∝ ν ^{− α} ). This makes up the Deep X-ray Radio Blazar Survey (DXRBS). Here we present new identifications and spectra for 106 sources, including 86 radio-loud quasars, 11 BL Lacertae objects, and nine narrow-line radio galaxies. Together with our previously published objects and already-known sources, our sample now contains 298 identified objects: 234 radio-loud quasars [181 flat-spectrum quasars: FSRQ ( α _r≤0.50) and 53 steep-spectrum quasars: SSRQ], 36 BL Lacs and 28 narrow-line radio galaxies. Redshift information is available for 96 per cent of these. Thus our selection technique is ∼90 per cent efficient at finding radio-loud quasars and BL Lacs. Reaching 5-GHz radio fluxes ∼50 mJy and 0.1–2.0 keV X-ray fluxes a few ×10⁻¹⁴ erg cm⁻² s⁻¹, DXRBS is the faintest and largest flat-spectrum radio sample with nearly complete (∼85 per cent) identification. We review the properties of the DXRBS blazar sample, including redshift distribution and coverage of the X-ray-radio–power plane for quasars and BL Lacs. Additionally, we touch upon the expanded multiwavelength view of blazars provided by DXRBS. By sampling for the first time the faint end of the radio and X-ray luminosity functions, this sample will allow us to investigate the blazar phenomenon and the validity of unified schemes down to relatively low powers.     

Hubble Space Telescope imaging of the CFRS and LDSS redshift surveys – IV. Influence of mergers in the evolution of faint field galaxies from z∼1

Hubble Space Telescope images of a sample of 285 galaxies with measured redshifts from the Canada–France Redshift Survey (CFRS) and Autofib–Low Dispersion Spectrograph Survey (LDSS) redshift surveys are analysed to derive the evolution of the merger fraction out to redshifts z ∼1. We have performed visual and machine-based merger identifications, as well as counts of bright pairs of galaxies with magnitude differences δm ≤1.5 mag. We find that the pair fraction increases with redshift, with up to ∼20 per cent of the galaxies being in physical pairs at z ∼0.75–1. We derive a merger fraction varying with redshift as ∝(1+ z )^3.2±0.6, after correction for line-of-sight contamination, in excellent agreement with the merger fraction derived from the visual classification of mergers for which m =3.4±0.6. After correcting for seeing effects on the ground-based selection of survey galaxies, we conclude that the pair fraction evolves as ∝(1+ z )^2.7±0.6. This implies that an average L * galaxy will have undergone 0.8–1.8 merger events from z =1 to z =0, with 0.5 to 1.2 merger events occuring in a 2-Gyr time-span at around z ∼0.9. This result is consistent with predictions from semi-analytical models of galaxy formation. From the simple coaddition of the observed luminosities of the galaxies in pairs, physical mergers are computed to lead to a brightening of 0.5 mag for each pair on average, and a boost in star formation rate of a factor of 2, as derived from the average [O  ii ] equivalent widths. Mergers of galaxies are therefore contributing significantly to the evolution of both the luminosity function and luminosity density of the Universe out to z ∼1.     

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.     

Cosmological evolution of the Fanaroff–Riley type II source population

By combining a model for the evolution of the radio luminosity of an individual source with the radio luminosity function, we perform a multidimensional Monte Carlo simulation to investigate the cosmological evolution of the Fanaroff–Riley type II (FR II) radio galaxy population by generating large artificial samples. The properties of FR II sources are required to evolve with redshift in the artificial samples to fit the observations. Either the maximum jet age or the maximum density of the jet environment or both evolve with redshift. We also study the distribution of FR II source properties as a function of redshift. From currently available data we cannot constrain the shape of the distribution of environment density or age, but jet power is found to follow a power-law distribution with an exponent of approximately −2. This power-law slope does not change with redshift out to   z = 0.6  . We also find the distribution of the pressure in the lobes of FR II sources to evolve with redshift up to   z ∼ 1.2  .     

The physical parameters of the evolving population of faint galaxies

The excess number of blue galaxies at faint magnitudes is a subject of much controversy. Recent Hubble Space Telescope results have revealed a plethora of galaxies with peculiar morphologies tentatively identified as the evolving population. We report the results of optical spectroscopy and near-infrared photometry of a sample of faint HST galaxies from the Medium Deep Survey to ascertain the physical properties of the faint morphological populations. We find four principal results. First, the population of objects classified as 'peculiar' are intrinsically luminous in the optical ( M _B  ∼ −19). Secondly these systems tend to be strong sources of [O  II ] line luminosity. Thirdly the optical–infrared colours of the faint population (a) confirm the presence of a population of compact   blue galaxies and (b) show the stellar populations of irregular/peculiar galaxies encompass a wide range in age. Finally a surface-brightness comparison with the local galaxy sample of Frei et al. shows that these objects are not of anomalously low surface brightness, rather we find that all morphological classes have evolved to a higher surface brightness at higher redshifts ( z  > 0.3).     

Morphology and surface brightness evolution of z∼1.1 radio galaxies

We use K '-band (2.1-μm) imaging to investigate the angular size and morphology of 10 6C radio galaxies, at redshifts 1≤ z ≤1.4. Two radio galaxies appear to be undergoing mergers, another contains, within a single envelope, two intensity peaks aligned with the radio jets, while the other seven appear consistent with being normal ellipticals in the K band.
Intrinsic half-light radii are estimated from the areas of each radio galaxy image above a series of thresholds. The 6C galaxy radii are found to be significantly smaller than those of the more radio-luminous 3CR galaxies at similar redshifts. This would indicate that the higher mean K -band luminosity of the 3CR galaxies reflects a difference in the size of the host galaxies, and not solely a difference in the power of the active nuclei.
The size–luminosity relation of the z ∼1.1 6C galaxies indicates a 1.0–1.6 mag enhancement of their rest frame R -band surface brightness relative to either local ellipticals of the same size or FRII radio galaxies at z <0.2. The 3CR galaxies at z ∼1.1 show a comparable enhancement in surface brightness. The mean radius of the 6C galaxies suggests that they evolve into ellipticals of L ∼ L * luminosity, and is consistent with their low-redshift counterparts being relatively small FRII galaxies ∼25 times lower in radio luminosity, or small FRI galaxies ∼1000 times lower in radio luminosity. Hence the 6C radio galaxies appear to undergo as much optical and radio evolution as the 3CR galaxies.     

The cluster environments of the z ∼ 1 3CR radio galaxies

An analysis of the environments around a sample of 28 3CR radio galaxies with redshifts 0.6< z <1.8 is presented, based primarily upon K -band images down to K ∼20 taken using the UK Infrared Telescope (UKIRT). A net overdensity of K -band galaxies is found in the fields of the radio galaxies, with the mean excess counts being comparable to that expected for clusters of Abell Class 0 richness. A sharp peak is found in the angular cross-correlation amplitude centred on the radio galaxies that, for reasonable assumptions about the luminosity function of the galaxies, corresponds to a spatial cross-correlation amplitude between those determined for low-redshift Abell Class 0 and 1 clusters.
These data are complemented by J -band images also from UKIRT, and by optical images from the Hubble Space Telescope . The fields of the lower redshift ( z ≲0.9) radio galaxies in the sample generally show well-defined near-infrared colour–magnitude relations with little scatter, indicating a significant number of galaxies at the redshift of the radio galaxy; the relations involving colours at shorter wavelengths than the 4000 Å break show considerably greater scatter, suggesting that many of the cluster galaxies have low levels of recent or on-going star formation. At higher redshifts the colour–magnitude sequences are less prominent owing to the increased field galaxy contribution at faint magnitudes, but there is a statistical excess of galaxies with the very red infrared colours ( J − K ≳1.75) expected of old cluster galaxies at these redshifts.
Although these results are appropriate for the mean of all of the radio galaxy fields, there exist large field-to-field variations in the richness of the environments. Many, but certainly not all, powerful z ∼1 radio galaxies lie in (proto)cluster environments.