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.     

The GLARE Survey – II. Faint z≈ 6 Lyα line emitters in the HUDF

The galaxy population at   z ≈ 6  has been the subject of intense study in recent years, culminating in the Hubble Ultra Deep Field (HUDF) – the deepest imaging survey yet. A large number of high-redshift galaxy candidates have been identified within the HUDF, but until now analysis of their properties has been hampered by the difficulty of obtaining spectroscopic redshifts for these faint galaxies. Our 'Gemini Lyman-Alpha at Reionization Era' (GLARE) project has been designed to undertake spectroscopic follow-up of faint  ( z ' < 28.5)  i '-drop galaxies at   z ≈ 6  in the HUDF. In a previous paper we presented preliminary results from the first 7.5 h of data from GLARE. In this paper we detail the complete survey. We have now obtained 36 h of spectroscopy on a single GMOS slitmask from Gemini-South, with a spectral resolution of  λ/Δλ_FWHM≈ 1000  . We identify five strong Lyα emitters at   z > 5.5  , and a further nine possible line emitters with detections at lower significance. We also place tight constraints on the equivalent width of Lyα emission for a further ten i '-drop galaxies and examine the equivalent width distribution of this faint spectroscopic sample of   z ≈ 6  galaxies. We find that the fraction of galaxies with little or no emission is similar to that at   z ≈ 3  , but that the   z ≈ 6  population has a tail of sources with high rest-frame equivalent widths. Possible explanations for this effect include a tendency towards stronger line emission in faint sources, which may arise from extreme youth or low metallicity in the Lyman-break population at high redshift, or possibly a top-heavy initial mass function.     

Understanding the redshift evolution of the luminosity functions of Lyman α emitters

We present a semi-analytical model of star formation which explains simultaneously the observed ultraviolet (UV) luminosity function (LF) of high-redshift Lyman break galaxies (LBGs) and LFs of Lyman α emitters. We consider both models that use the Press–Schechter (PS) and Sheth–Tormen (ST) halo mass functions to calculate the abundances of dark matter haloes. The Lyman α LFs at   z ≲ 4  are well reproduced with only ≲10 per cent of the LBGs emitting Lyman α lines with rest equivalent width greater than the limiting equivalent width of the narrow band surveys. However, the observed LF at   z > 5  can be reproduced only when we assume that nearly all LBGs are Lyman α emitters. Thus, it appears that  4 < z < 5  marks the epoch when a clear change occurs in the physical properties of the high-redshift galaxies. As Lyman α escape depends on dust and gas kinematics of the interstellar medium (ISM), this could mean that on an average the ISM at   z > 5  could be less dusty, more clumpy and having more complex velocity field. All of these will enable easier escape of the Lyman α photons. At   z > 5  , the observed Lyman α LF are well reproduced with the evolution in the halo mass function along with very minor evolution in the physical properties of high-redshift galaxies. In particular, up to   z = 6.5  , we do not see the effect of evolving intergalactic medium opacity on the Lyman α escape from these galaxies.     

The relationship between star formation rate and radio synchrotron luminosity at 0 < z < 2

We probe the relationship between star formation rate (SFR) and radio synchrotron luminosity in galaxies at  0 < z < 2  within the northern Spitzer Wide-area Infrared Extragalactic survey (SWIRE) fields, in order to investigate some of the assumptions that go into calculating the star formation history of the Universe from deep radio observations. We present new 610-MHz Giant Metrewave Radio Telescope (GMRT) observations of the European Large-Area ISO Survey-North 2 (ELAIS-N2) field, and using this data, along with previous GMRT surveys carried out in the ELAIS-N1 (North 1) and Lockman Hole regions, we construct a sample of galaxies which have redshift and SFR information available from the SWIRE survey. We test whether the local relationship between SFR and radio luminosity is applicable to   z = 2  galaxies, and look for evolution in this relationship with both redshift and SFR in order to examine whether the physical processes which lead to synchrotron radiation have remained the same since the peak of star formation in the Universe. We find that the local calibration between radio luminosity and star formation can be successfully applied to radio-selected high-redshift, high-SFR galaxies, although we identify a small number of sources where this may not be the case; these sources show evidence for inaccurate estimations of their SFR, but there may also be some contribution from physical effects such as the recent onset of starburst activity, or suppression of the radio luminosity within these galaxies.     

Measurement of dwarf galaxies in the rich clusters Abell 665 and 963 at z=0.2

We show that the luminosity functions of the distant rich clusters Abell 665 ( z =0.182) and Abell 963 ( z =0.206) are flat or gradually rising down to M_R =−14, with α≈−1.2±0.4 [here α is the logarithmic slope of the luminosity function: φ( L )∝ L ^α at the faint end]. We do not confirm the steep luminosity functions (α≤−1.8) that have been recently proposed for these two clusters.
Several technical points are discussed in detail. In particular, we compute the corrections to the background contamination caused by gravitational lensing from the cluster dark matter, and show that the corrections are small unless we wish to determine variations in the luminosity function on small scales.
Recent observations have also shown that the field galaxy luminosity function at z ≈0.2 is also shallow between M_B =−19 and M_B =−13. Abell 665 and 963 are two of the richest clusters known at that redshift. We therefore propose that the galaxy luminosity function might be universal in this magnitude range at z =0.2.
The dwarf galaxies that we see in Abell 665 have a colour distribution that is strongly peaked at B − R =1.9. We compute K -corrections based on the spectral energy distributions of local galaxies, and show that these are probably dwarf spheroidal galaxies. This might suggest that the dwarf spheroidal population observed in Virgo already existed at z =0.2.     

Host galaxy morphologies of X-ray selected AGN: assessing the significance of different black hole fuelling mechanisms to the accretion density of the Universe at z∼ 1.

We use morphological information of X-ray selected active galactic nuclei (AGN) hosts to set limits on the fraction of the accretion density of the Universe at   z ≈ 1  that is not likely to be associated with major mergers. Deep X-ray observations are combined with high-resolution optical data from the Hubble Space Telescope in the All-wavelength Extended Groth strip International Survey, Great Observatories Origins Deep Survey (GOODS) North and GOODS South fields to explore the morphological breakdown of X-ray sources in the redshift interval  0.5 < z < 1.3  . The sample is split into discs, early-type bulge-dominated galaxies, peculiar systems and point sources in which the nuclear source outshines the host galaxy. The X-ray luminosity function and luminosity density of AGN at   z ≈ 1  are then calculated as a function of morphological type. We find that disc-dominated hosts contribute  30 ± 9  per cent to the total AGN space density and  23 ± 6  per cent to the luminosity density at   z ≈ 1  . We argue that AGN in disc galaxies are most likely fuelled not by major merger events but by minor interactions or internal instabilities. We find evidence that these mechanisms may be more efficient in producing luminous AGN     compared to predictions for the stochastic fuelling of massive black holes in disc galaxies.     

A limit on the number density of bright z≈ 7 galaxies

We present a survey of bright optical dropout sources in two deep, multiwavelength surveys comprising 11 widely separated fields, aimed at constraining the galaxy luminosity function at   z ≈ 7  for sources at  5–10  L * ( z = 6)  . Our combined survey area is 225 arcmin² to a depth of   J _AB= 24.2  (3σ) and 135 arcmin² to   J = 25.3  (4σ). We find that infrared data longwards of 2 μm are essential for classifying optical dropout sources, and in particular for identifying cool Galactic star contaminants. Our limits on the number density of high-redshift sources are consistent with current estimates of the Lyman break galaxy luminosity function at   z = 6  .     

The local far-infrared galaxy colour–luminosity distribution: a reference for BLAST and Herschel/SPIRE submillimetre surveys

We measure the local galaxy far-infrared (FIR) 60 to 100 μm colour–luminosity distribution using an all-sky IRAS survey. This distribution is an important reference for the next generation of FIR–submillimetre surveys that have and will conduct deep extragalactic surveys at 250–500 μm. With the peak in dust-obscured star-forming activity leading to present-day giant ellipticals now believed to occur in submillimetre galaxies near   z ∼ 2.5  , these new FIR–submillimetre surveys will directly sample the spectral energy distributions of these distant objects at rest-frame FIR wavelengths similar to those at which local galaxies were observed by IRAS . We have taken care to correct for the temperature bias and the evolution effects in our IRAS 60-μm-selected sample. We verify that our colour–luminosity distribution is consistent with the measurements of the local FIR luminosity function, before applying it to the higher redshift Universe. We compare our colour–luminosity correlation with recent dust–temperature measurements of submillimetre galaxies and find evidence for pure luminosity evolution of the form  (1 + z )³  . This distribution will be useful for the development of evolutionary models for Balloon-borne Large Aperture Submillimeter Telescope (BLAST) and Spectral and Photometric Imaging Receiver (SPIRE) surveys as it provides a statistical distribution of the rest-frame dust temperatures for galaxies as a function of luminosity.     

Candidate high-redshift and primeval galaxies in Hubble Deep Field South

We present the results of colour selection of candidate high-redshift galaxies in Hubble Deep Field South (HDF-S) using the Lyman dropout scheme. The HDF-S data we discuss were taken in a number of different filters extending from the near-UV (F300W) to the infrared (F222M) in two different fields. This allows us to select candidates with redshifts from z ∼3 to z ∼12. We find 15 candidate z ∼3 objects (F300W dropouts), one candidate z ∼4 object (F450W dropout) and 16 candidate z ∼5 objects (F606W dropouts) in the ∼4.7‐arcmin² WFPC-2 field, and four candidate z ∼6 objects (optical dropouts) and one candidate z ∼8 object (F110W dropout) in the 0.84-arcmin² NICMOS-3 field. No F160W dropouts are found ( z ∼12). We compare our selection technique with existing data for Hubble Deep Field North (HDF-N) and discuss alternative interpretations of the objects. We conclude that there are a number of lower redshift interlopers in the selections, including one previously identified object, and reject those objects most likely to be foreground contaminants. Even after this we conclude that the F606W dropout list is likely to still contain substantial foreground contamination. The lack of candidate very-high-redshift UV-luminous galaxies supports earlier conclusions by Lanzetta et al. We discuss the morphologies and luminosity functions of the high-redshift objects, and their cosmological implications.     

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 Hα luminosity function and star formation rate up to z ∼ 1

We describe ISAAC/ESO-VLT observations of the Hαλ6563 Balmer line of 33 field galaxies from the Canada–France Redshift Survey (CFRS) with redshifts selected between 0.5 and 1.1. We detect Hα in emission in 30 galaxies and compare the properties of this sample with the low-redshift sample of CFRS galaxies at   z ∼ 0.2  . We find that the Hα luminosity,   L (Hα)  , is tightly correlated to   M ( B _AB)  in the same way for both the low- and high-redshift samples.   L (Hα)  is also correlated to L ([O  ii ]λ3727), and again the relation appears to be similar at low and high redshifts. The ratio L (lsqb;O  ii ])/   L (Hα)  decreases for brighter galaxies by as much as a factor of 2 on average. Derived from the Hα luminosity function, the comoving Hα luminosity density increases by a factor 12 from  〈 z 〉= 0.2  to  〈 z 〉= 1.3  . Our results confirm a strong rise of the star formation rate (SFR) at   z < 1.3  , proportional to  (1 + z )^4.1±0.3  (with   H ₀= 50 km s⁻¹ Mpc⁻¹, q ₀= 0.5  ). We find an average  SFR(2800 Å)/SFR (Hα)  ratio of 3.2 using the Kennicutt SFR transformations. This corresponds to the dust correction that is required to make the near-ultraviolet data consistent with the reddening-corrected Hα data within the self-contained, I -selected CFRS sample.     

An ultraviolet-selected galaxy redshift survey – II. The physical nature of star formation in an enlarged sample

We present further spectroscopic observations for a sample of galaxies selected in the vacuum ultraviolet (UV) at 2000 Å from the FOCA balloon-borne imaging camera of Milliard et al. This work represents an extension of the initial study by Treyer et al. Our enlarged catalogue contains 433 sources (≃3 times as many as in our earlier study) across two FOCA fields. 273 of these are galaxies, nearly all with redshifts z ≃0–0.4. Nebular emission-line measurements are available for 216 galaxies, allowing us to address issues of excitation, reddening and metallicity. The UV and H α luminosity functions strengthen our earlier assertions that the local volume-averaged star formation rate is higher than indicated from earlier surveys. Moreover, internally within our sample, we do not find a steep rise in the UV luminosity density with redshift over 0< z <0.4. Our data are more consistent with a modest evolutionary trend, as suggested by recent redshift survey results. Investigating the emission-line properties, we find no evidence for a significant number of AGN in our sample; most UV-selected sources to z ≃0.4 are intense star-forming galaxies. We find that the UV flux indicates a consistently higher mean star formation rate than that implied by the H α luminosity for typical constant or declining star formation histories. Following Glazebrook et al., we interpret this discrepancy in terms of a starburst model for our UV-luminous sources. We develop a simple algorithm which explores the scatter in the UV flux–H α relation in the context of various burst scenarios. Whilst we can explain most of our observations in this way, there remains a small population with extreme UV–optical colours which cannot be understood.     

The X-ray luminosity function of AGN at z∼ 3

We combine Lyman-break colour selection with ultradeep (≳200 ks) Chandra X-ray imaging over a survey area of ∼0.35 deg² to select high-redshift active galactic nuclei (AGN). Applying careful corrections for both the optical and X-ray selection functions, the data allow us to make the most accurate determination to date of the faint end of the X-ray luminosity function (XLF) at   z ∼ 3  . Our methodology recovers a number density of X-ray sources at this redshift which is at least as high as previous surveys, demonstrating that it is an effective way of selecting high z AGN. Comparing to results at   z = 1  , we find no evidence that the faint slope of the XLF flattens at high z , but we do find significant (factor ∼3.6) negative evolution of the space density of low luminosity AGN. Combining with bright end data from very wide surveys we also see marginal evidence for continued positive evolution of the characteristic break luminosity   L _*  . Our data therefore support models of luminosity-dependent density evolution between   z = 1  and   z = 3  . A sharp upturn in the the XLF is seen at the very lowest luminosities  ( L _X≲ 10^42.5 erg s⁻¹)  , most likely due to the contribution of pure X-ray starburst galaxies at very faint fluxes.     

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.     

Evolution of the u-band luminosity function from redshift 1.2 to 0

We produce and analyse u -band (  λ≈ 355  nm) luminosity functions (LFs) for the red and blue populations of galaxies using data from the Sloan Digital Sky Survey (SDSS) u -band Galaxy Survey ( u GS) and Deep Evolutionary Exploratory Probe 2 (DEEP2) survey. From a spectroscopic sample of 41 575 SDSS u GS galaxies and 24 561 DEEP2 galaxies, we produce colour magnitude diagrams and make use of the colour bimodality of galaxies to separate red and blue populations. LFs for eight redshift slices in the range  0.01 < z < 1.2  are determined using the  1/ V _max  method and fitted with Schechter functions showing that there is significant evolution in   M *  , with a brightening of 1.4 mag for the combined population. The integration of the Schechter functions yields the evolution in the u -band luminosity density (LD) out to   z ∼ 1  . By parametrizing the evolution as  ρ∝ (1 + z )^β  , we find that  β= 1.36 ± 0.2  for the combined populations and  β= 2.09 ± 0.2  for the blue population. By removing the contribution of the old stellar population to the u -band LD and correcting for dust attenuation, we estimate the evolution in the star formation rate (SFR) of the Universe to be  β_SFR= 2.5 ± 0.3  . Discrepancies between our result and higher evolution rates measured using the infrared and far-UV can be reconciled by considering possibilities such as an underestimated dust correction at high redshifts or evolution in the stellar initial mass function.     

Light curves for off-centre ignition models of Type Ia supernovae

We compare observations of the high-redshift galaxy population to the predictions of the galaxy formation model of Croton et al. and De Lucia & Blaizot. This model, implemented on the Millennium Simulation of the concordance Lambda cold dark matter cosmogony, introduces 'radio mode' feedback from the central galaxies of groups and clusters in order to obtain quantitative agreement with the luminosity, colour, morphology and clustering properties of the present-day galaxy population. Here we construct deep light cone surveys in order to compare model predictions to the observed counts and redshift distributions of distant galaxies, as well as to their inferred luminosity and mass functions out to redshift 5. With the exception of the mass functions, all these properties are sensitive to modelling of dust obscuration. A simple but plausible treatment agrees moderately well with most of the data. The predicted abundance of relatively massive  (∼ M _*)  galaxies appears systematically high at high redshift, suggesting that such galaxies assemble earlier in this model than in the real Universe. An independent galaxy formation model implemented on the same simulation matches the observed mass functions slightly better, so the discrepancy probably reflects incomplete or inaccurate galaxy formation physics rather than problems with the underlying cosmogony.     

Clustering of galaxies around radio quasars at 0.5≤z≤0.8

We have observed the galaxy environments around a sample of 21 radio-loud, steep-spectrum quasars at 0.5≤ z ≤0.82, spanning several orders of magnitude in radio luminosity. The observations also include background control fields used to obtain the excess number of galaxies in each quasar field. The galaxy excess was quantified using the spatial galaxy–quasar correlation amplitude, B _gq, and an Abell-type measurement, N _0.5. A few quasars are found in relatively rich clusters, but on average, they seem to prefer galaxy groups or clusters of approximately Abell class 0. We have combined our sample with literature samples extending down to z ≈0.2 and covering the same range in radio luminosity. By using the Spearman statistic to disentangle redshift and luminosity dependences, we detect a weak, but significant, positive correlation between the richness of the quasar environment and the radio luminosity of the quasar. However, we do not find any epoch dependence in B _gq, as has previously been reported for radio quasars and galaxies. We discuss the radio luminosity–cluster richness link and possible explanations for the weak correlation that is seen.     

A submillimetre survey of the star formation history of radio galaxies

We present the results of the first major systematic submillimetre survey of radio galaxies spanning the redshift range 1< z <5. The primary aim of this work is to elucidate the star formation history of this sub class of elliptical galaxies by tracing the cosmological evolution of dust mass. Using SCUBA on the JCMT, we have obtained 850-μm photometry of 47 radio galaxies to a consistent rms depth of 1 mJy, and have detected dust emission in 14 cases. The radio galaxy targets have been selected from a series of low-frequency radio surveys of increasing depth (3CRR, 6CE, etc.), in order to allow us to separate the effects of increasing redshift and increasing radio power on submillimetre luminosity. Although the dynamic range of our study is inevitably small, we find clear evidence that the typical submillimetre luminosity (and hence dust mass) of a powerful radio galaxy is a strongly increasing function of redshift; the detection rate rises from ≃15 per cent at z <2.5 to ≳75 per cent at z >2.5, and the average submillimetre luminosity rises at a rate ∝(1+ z )³ out to z ≃4. Moreover, our extensive sample allows us to argue that this behaviour is not driven by underlying correlations with other radio galaxy properties such as radio power, radio spectral index, or radio source size/age. Although radio selection may introduce other more subtle biases, the redshift distribution of our detected objects is in fact consistent with the most recent estimates of the redshift distribution of comparably bright submillimetre sources discovered in blank field surveys. The evolution of submillimetre luminosity found here for radio galaxies may thus be representative of massive ellipticals in general.     

The luminosity function and surface brightness distribution of H i selected galaxies

We measure the     B -band optical luminosity function (LF) for galaxies selected in a blind H  i survey. The total LF of the H  i selected sample is flat, with Schechter parameters     and     , in good agreement with LFs of optically selected late-type galaxies. Bivariate distribution functions of several galaxy parameters show that the H  i density in the local Universe is more widely spread over galaxies of different size, central surface brightness and luminosity than the optical luminosity density is. The number density of very low surface brightness (LSB ) (>24.0 mag arcsec⁻²) gas-rich galaxies is considerably lower than that found in optical surveys designed to detect dim galaxies. This suggests that only a part of the population of LSB galaxies is gas-rich and that the rest must be gas-poor. However, we show that this gas-poor population must be cosmologically insignificant in baryon content. The contribution of gas-rich LSB galaxies (>23.0 mag arcsec⁻²) to the local cosmological gas and luminosity density is modest     and     per cent respectively); their contribution to Ω_matter is not well-determined, but probably <11 per cent. These values are in excellent agreement with the low redshift results from the Hubble Deep Field.     

Restrictions to the galaxy evolutionary models from the Hawaiian Deep Fields SSA13 and SSA22

Quantitative structural analysis of the galaxies present in the Hawaiian Deep Fields SSA13 and SSA22 is reported. The structural parameters of the galaxies have been obtained automatically by fitting a two-component model (Sérsic r ^{1/ n} bulge and exponential disc) to the surface brightness of the galaxies. The galaxies were classified on the basis of the bulge-to-total luminosity ratio  ( B / T )  . The magnitude selection criteria and the reliability of our method have been checked by using Monte Carlo simulations. A complete sample of objects up to redshift 0.8 has been achieved. Spheroidal objects (E/S0) represent ≈33 per cent and spirals ≈41 per cent of the total number of galaxies, while mergers and unclassified objects represent ≈26 per cent. We have computed the comoving space density of the different kinds of object. In an Einstein–de Sitter universe, a decrease in the comoving density of E/S0 galaxies is observed as redshift increases (≈30 per cent less at   z =0.8)  , while for spiral galaxies a relatively quiet evolution is reported. The framework of hierarchical clustering evolution models of galaxies seems to be the most appropriate to explain our results.