Confirmation of the effectiveness of submm source redshift estimation based on rest-frame radio–FIR photometry

We present a comparison between the published optical, infrared (IR) and CO spectroscopic redshifts of 15 (sub)mm galaxies and their photometric redshifts as derived from long-wavelength (radio–mm–far-IR) photometric data. The redshift accuracy measured for 12 submillimetre (submm) galaxies with at least one robustly determined colour in the radio–mm–far-IR regime is  δ z ≈ 0.30 (rms)  . Despite the wide range of spectral energy distributions in the local galaxies that are used in an unbiased manner as templates, this analysis demonstrates that photometric redshifts can be efficiently derived for submm galaxies with a precision of  δ z < 0.5  using only the rest-frame far-IR to radio wavelength data.     

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.     

MegaZ-LRG: a photometric redshift catalogue of one million SDSS luminous red galaxies

We describe the construction of MegaZ-LRG, a photometric redshift catalogue of over one million luminous red galaxies (LRGs) in the redshift range  0.4 < z < 0.7  with limiting magnitude   i < 20  . The catalogue is selected from the imaging data of the Sloan Digital Sky Survey (SDSS) Data Release 4. The 2dF-SDSS LRG and Quasar (2SLAQ) spectroscopic redshift catalogue of 13 000 intermediate-redshift LRGs provides a photometric redshift training set, allowing use of ann z, a neural network-based photometric-redshift estimator. The rms photometric redshift accuracy obtained for an evaluation set selected from the 2SLAQ sample is  σ _z = 0.049  averaged over all galaxies, and  σ _z = 0.040  for a brighter subsample  ( i < 19.0)  . The catalogue is expected to contain ∼5 per cent stellar contamination. The ann z code is used to compute a refined star/galaxy probability based on a range of photometric parameters; this allows the contamination fraction to be reduced to 2 per cent with negligible loss of genuine galaxies. The MegaZ-LRG catalogue is publicly available on the World Wide Web from http://www.2slaq.info .     

More redshifts of powerful equatorial radio sources from the Best, Röttgering & Lehnert sample

A new sample of very powerful radio sources, defined from the Molonglo Reference Catalogue, was recently compiled by Best, Röttgering & Lehnert. These authors provided redshifts for 174 of the 178 objects in the sample, making the sample 98 per cent spectroscopically complete. Here, redshifts for three of the remaining galaxies are presented, confirming the optical identifications and raising the spectroscopic completeness of the sample to 99.5 per cent; only 1059−010 (3C 249) is currently without redshift.     

Reconstructing galaxy fundamental distributions and scaling relations from photometric redshift surveys. Applications to the SDSS early-type sample

Noisy distance estimates associated with photometric rather than spectroscopic redshifts lead to a biased estimate of the luminosity distribution, and produce a correlated misestimate of the sizes. We consider a sample of early-type galaxies from the Sloan Digital Sky Survey Data Release 6 for which both spectroscopic and photometric information is available, and apply the generalization of the V _max method to correct for these biases. We show that our technique recovers the true redshift, magnitude and size distributions, as well as the true size–luminosity relation. We find that using only 10 per cent of the spectroscopic information randomly spaced in our catalogue is sufficient for the reconstructions to be accurate within  ∼3 per cent  , when the photometric redshift error is  δ z ≃ 0.038  . We then address the problem of extending our method to deep redshift catalogues, where only photometric information is available. In addition to the specific applications outlined here, our technique impacts a broader range of studies, when at least one distance-dependent quantity is involved. It is particularly relevant for the next generation of surveys, some of which will only have photometric information.     

The ultraluminous and hyperluminous infrared galaxies in the Sloan Digital Sky Survey, 2dF Galaxy Redshift Survey and 6dF Galaxy Survey

We present a result of cross-correlating the Infrared Astronomical Satellite Faint Source Catalogue with the spectroscopic catalogues of galaxies in the Fourth Data Release of the Sloan Digital Sky Survey, the Final Data Release of the 2dF Galaxy Redshift Survey (2dFGRS) and the Second Data Release of the 6dF Galaxy Survey. We have identified 324 ultraluminous infrared galaxies (ULIRGs) including 190 newly discovered ULIRGs, and two hyperluminous infrared galaxies. Adding these new ULIRGs, we increase the number of known ULIRGs by about 30 per cent. The reliability of the cross-correlation is estimated using the likelihood ratio method. The incompleteness of our sample introduced by the identification procedure in this study is estimated to be about 5 per cent. Our sample covers the redshift range of   z = 0.037–0.517  with a median redshift of     , which is larger than that     of the sample of previously known ULIRGs.     

The contribution of faint blue galaxies to the submillimetre counts and background

Observations in the submillimetre (submm) waveband have recently revealed a new population of luminous sources. These are proposed to lie at high redshift and to be optically faint because of their high intrinsic dust obscuration. The presence of dust has been previously invoked in optical galaxy count models which use the Bruzual & Charlot evolution models with an exponential τ =9 Gyr star formation rate (SFR) for spirals, and these fit the count data well from U to K . We now show that by using either a 1/ λ or Calzetti absorption law for the dust and re-distributing the evolved spiral galaxy ultraviolet (UV) radiation into the far-infrared (FIR), these models can account for all of the 'faint' ( 1 mJy) 850-μm galaxy counts, but fail to fit 'bright' ( 2 mJy) sources, indicating that another explanation for the submm counts may apply at brighter fluxes, e.g., quasi-stellar objects (QSOs) or ultraluminous infrared galaxies (ULIRGs). We find that the main contribution to the faint, submm number counts is in the redshift range 0.5< z <3, peaking at z ≈1.8. The above model, using either dust law, can also explain a significant proportion of the extragalactic background at 850 μm, as well as producing a reasonable fit to the bright 60-μm IRAS counts.     

The redshift-space two-point correlation function of ELAIS-S1 galaxies

We investigate the clustering properties of galaxies in the recently completed ELAIS-S1 redshift survey through their spatial two-point autocorrelation function. We used a subsample of the ELAIS-S1 catalogue covering approximately 4 deg² and consisting of 148 objects selected at 15 μm with a flux >0.5 mJy and a redshift   z < 0.5  . We detected a positive signal in the correlation function that in the range of separations  1–10  h ⁻¹ Mpc  is well approximated by a power law with a slope  γ= 1.4 ± 0.25  and a correlation length   s ₀= 5.4 ± 1.2  h ⁻¹ Mpc  , at the 90 per cent significance level. This result is in good agreement with the redshift-space correlation function measured in more local samples of mid-infrared-selected galaxies such as the IRAS Point Source Catalog (PSC z ) redshift survey. This suggests a lack of significant clustering evolution of infrared-selected objects out to   z = 0.5  that is further confirmed by the consistency found between the correlation functions measured in a local  ( z < 0.2)  and a distant  (0.2 < z < 0.5)  subsample of ELAIS-S1 galaxies. We also confirm that optically selected galaxies in the local redshift surveys, especially those of the SDSS sample, are significantly more clustered than infrared objects.     

The PSCz catalogue

We present the catalogue, mask, redshift data and selection function for the PSC z survey of 15 411 IRAS galaxies across 84 per cent of the sky. Most of the IRAS data are taken from the Point Source Catalog, but this has been supplemented and corrected in various ways to improve the completeness and uniformity. We quantify the known imperfections in the catalogue, and we assess the overall uniformity, completeness and data quality. We find that overall the catalogue is complete and uniform to within a few per cent at high latitudes and 10 per cent at low latitudes. Ancillary information, access details, guidelines and caveats for using the catalogue are given.     

A 98 per cent spectroscopically complete sample of the most powerful equatorial radio sources at 408 MHz

A new sample of very powerful radio galaxies is defined from the Molonglo Reference Catalogue, according to the criteria S _408 MHz>5 Jy, −30°≤ δ ≤10° and | b |≥10°. The sample is selected to have similar properties to the northern 3CR revised sample, and to be visible to a combination of existing northern telescopes such as the Very Large Array radio interferometer and large southern hemisphere telescope facilities. The sample contains 178 sources, of which spectroscopic redshifts are available in the literature for 128. For the remaining 50 sources, new radio imaging, optical imaging and spectroscopic observations are presented to identify the host galaxies and determine their redshifts. With these new observations the total sample is 100 per cent optically identified and redshifts are available for 174 (98 per cent) of the sources. The sample consists of one starburst galaxy, one Seyfert galaxy, 127 radio galaxies and 49 quasars. Basic properties of the sample, such as the distributions of the quasar and radio-galaxy populations in redshift and their locations on the radio power versus linear size ( P − D ) diagram, show no significant differences from the revised 3CR sample. The equatorial location and the high spectroscopic completeness of this sample make it a valuable resource for detailed studies of the nature and environments of these important objects with the new generation of southern hemisphere telescopes.     

Estimating the redshift distribution of photometric galaxy samples – II. Applications and tests of a new method

In Lima et al. we presented a new method for estimating the redshift distribution,   N ( z )  , of a photometric galaxy sample, using photometric observables and weighted sampling from a spectroscopic subsample of the data. In this paper, we extend this method and explore various applications of it, using both simulations and real data from the Sloan Digital Sky Survey (SDSS). In addition to estimating the redshift distribution for an entire sample, the weighting method enables accurate estimates of the redshift probability distribution,   p ( z )  , for each galaxy in a photometric sample. Use of   p ( z )  in cosmological analyses can substantially reduce biases associated with traditional photometric redshifts, in which a single redshift estimate is associated with each galaxy. The weighting procedure also naturally indicates which galaxies in the photometric sample are expected to have accurate redshift estimates, namely those that lie in regions of photometric-observable space that are well sampled by the spectroscopic subsample. In addition to providing a method that has some advantages over standard photo- z estimates, the weights method can also be used in conjunction with photo- z estimates e.g. by providing improved estimation of   N ( z )  via deconvolution of   N ( z _phot)  and improved estimates of photo- z scatter and bias. We present a publicly available   p ( z )  catalogue for ∼78 million SDSS DR7 galaxies.     

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.     

Starburst galaxies and structure in the submillimetre background towards the Hubble Deep Field

We use an 850-μm SCUBA map of the Hubble Deep Field (HDF) to study the dust properties of optically-selected starburst galaxies at high redshift. The optical/infrared (IR) data in the HDF allow a photometric redshift to be estimated for each galaxy, together with an estimate of the visible star-formation rate. The 850-μm flux density of each source provides the complementary information: the amount of hidden, dust-enshrouded star formation activity. Although the 850-μm map does not allow detection of the majority of individual sources, we show that the galaxies with the highest UV star-formation rates are detected statistically, with a flux density of about S ₈₅₀=0.2 mJy for an apparent UV star-formation rate of 1  h ⁻² M_⊙ yr⁻¹. This level of submillimetre output indicates that the total star-forming activity is on average a factor of approximately 6 times larger than the rate inferred from the UV output of these galaxies. The general population of optical starbursts is then predicted to contribute at least 25 per cent of the 850-μm background. We carry out a power-spectrum analysis of the map, which yields some evidence for angular clustering of the background source population, but at a level lower than that seen in Lyman-break galaxies. Together with other lines of argument, particularly from the NICMOS HDF data, this suggests that the 850-μm background originates over an extremely wide range of redshifts – perhaps 1≲ z ≲6.     

An AzTEC 1.1 mm survey of the GOODS-N field – II. Multiwavelength identifications and redshift distribution

We present results from a multiwavelength study of 29 sources (false detection probabilities <5 per cent) from a survey of the Great Observatories Origins Deep Survey-North (GOODS-N) field at 1.1 mm using the Astronomical Thermal Emission Camera (AzTEC). Comparing with existing 850 μm Submillimetre Common-User Bolometer Array (SCUBA) studies in the field, we examine differences in the source populations selected at the two wavelengths. The AzTEC observations uniformly cover the entire survey field to a 1σ depth of ∼1 mJy. Searching deep 1.4 GHz Very Large Array (VLA) and Spitzer 3–24 μm catalogues, we identify robust counterparts for 21 1.1 mm sources, and tentative associations for the remaining objects. The redshift distribution of AzTEC sources is inferred from available spectroscopic and photometric redshifts. We find a median redshift of   z = 2.7  , somewhat higher than   z = 2.0  for  850 μm  selected sources in the same field, and our lowest redshift identification lies at a spectroscopic redshift   z = 1.1460  . We measure the 850 μm to 1.1 mm colour of our sources and do not find evidence for '850 μm dropouts', which can be explained by the low signal-to-noise ratio of the observations. We also combine these observed colours with spectroscopic redshifts to derive the range of dust temperatures T , and dust emissivity indices β for the sample, concluding that existing estimates   T ∼ 30 K  and  β∼ 1.75  are consistent with these new data.     

The SCUBA HAlf Degree Extragalactic Survey (SHADES) – VII. Optical/IR photometry and stellar masses of submillimetre galaxies

We present estimates of the photometric redshifts, stellar masses and star formation histories of sources in the Submillimetre Common-User Bolometer Array (SCUBA) HAlf Degree Extragalactic Survey (SHADES). This paper describes the 60 SCUBA sources detected in the Lockman Hole covering an area of ∼320 arcmin². Using photometry spanning the B band to 8 μm, we find that the average SCUBA source forms a significant fraction of its stars in an early period of star formation and that most of the remainder forms in a shorter more intense burst around the redshift it is observed. This trend does not vary significantly with source redshift. However, the sources show a clear increase in stellar mass with redshift, consistent with downsizing. In terms of spectral energy distribution types, only two out of the 51 sources we have obtained photometric redshifts for are best fitted by a quasar-like spectrum, with approximately 80 per cent of the sources being best fitted with late-type spectra (Sc, Im and starburst). By including photometry at 850 μm, we conclude that the average SCUBA source is forming stars at a rate somewhere between 6 and 30 times the rate implied from the rest-frame optical in a dust obscured burst and that this burst creates 15–65 per cent of the total stellar mass. Using a simplistic calculation, we estimate from the average star formation history that between one in five and one in 15 bright  ( L _*+ 2 < L _optical < L _*− 1 mag)  galaxies in the field over the interval  0 < z < 3  will at some point in their lifetime experience a similar energetic dusty burst of star formation. Finally, we compute the evolution of the star formation rate density and find it peaks around   z ∼ 2  .     

The SCUBA Half Degree Extragalactic Survey (SHADES) – IX. The environment, mass and redshift dependence of star formation

We present a comparison between the SCUBA (Submillimetre Common User Bolometer Array) Half Degree Extragalactic Survey (SHADES) at 450 and  850 μm  in the Lockman Hole East with a deep Spitzer Space Telescope survey at  3.6–24 μm  conducted in guaranteed time. Using stacking analyses we demonstrate a striking correspondence between the galaxies contributing the submm extragalactic background light, with those likely to dominate the backgrounds at Spitzer wavelengths. Using a combination BRIzK plus Spitzer photometric redshifts, we show that at least a third of the Spitzer -identified submm galaxies at  1 < z < 1.5  appear to reside in overdensities when the density field is smoothed at 0.5–2 Mpc comoving diameters, supporting the high-redshift reversal of the local star formation–galaxy density relation. We derive the dust-shrouded cosmic star formation history of galaxies as a function of assembled stellar masses. For model stellar masses  <10¹¹ M_⊙  , this peaks at lower redshifts than the ostensible   z ∼ 2.2  maximum for submm point sources, adding to the growing consensus for 'downsizing' in star formation. Our surveys are also consistent with 'downsizing' in mass assembly. Both the mean star formation rates  〈d M _*/d t 〉  and specific star formation rates  〈(1/ M _*) d M _*/d t 〉  are in striking disagreement with some semi-analytic predictions from the Millenium Simulation. The discrepancy could either be resolved with a top-heavy initial mass function, or a significant component of the submm flux heated by the interstellar radiation field.     

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.     

Infrared photometric study of type II quasars

For this paper, we collected almost all of the type II quasars so far discovered. Among them, 485 sources have photometric data at JHK bands, mainly from Two-Micron All-Sky Survey observations, 65 sources have IRAS photometric data in at least one of the three IRAS bands at 25, 60 and 100 μm and 15 sources have IRAS photometric data in all three IRAS bands. We find that in nearly half of all type II quasars, both the near- and far-infrared radiations are dominated by starlight or thermal reprocessing of starlight by dust in the underlying galaxy. The infrared radiation of the other group (slightly over half) is dominated by a non-thermal radiation in the near-infrared, and mostly in the far-infrared also (although there is a mixture particularly for the longer wavelengths). It is proposed that for the later group, hidden broad lines may exist in the infrared. On the basis of our and previous results, we also discuss the possibility that there are two distinct classes of type II quasars: 'true' type II quasars without a broad-line region and heavily obscured type I quasars, in full analogy with the case for type II Seyfert galaxies. No relationships can be found for either the near-infrared or the far-infrared colours and the redshift. Correlations between absolute magnitude in the near- and far-infrared with redshift are found, but could be due to a flux limit (Malmquist bias).     

On the lack of cold dust in IRAS P09104+4109 and IRAS F15307+3252: their spectral energy distributions and implications for finding dusty AGNs at high redshift

We present upper limits on the 850-μm and 450-μm fluxes of the warm hyperluminous (bolometric luminosity     galaxies IRAS P09104+4109     and IRAS F15307+3252     , derived from measurements using the SCUBA bolometer array on the James Clerk Maxwell Telescope. Hot luminous infrared sources like these are thought to differ from more normal cold ultraluminous infrared     galaxies in that they derive most of their bolometric luminosities from dusty active galactic nuclei (AGNs) as opposed to starbursts. Such hot, dusty AGNs at high redshift are thought to be responsible for much of the mass accretion of the Universe that is in turn responsible for the formation of the supermassive black holes seen in the centres of local galaxies. The galaxy IRAS P09104+4109 is also unusual in that it is a cD galaxy in the centre of a substantial cooling-flow cluster, not an isolated interacting galaxy like most ultraluminous infrared galaxies. Previously it was known to have large amounts of hot     dust from IRAS observations. We now show that the contribution of cold dust to the bolometric luminosity is less than 3 per cent. Most ultraluminous infrared galaxies possess large amounts of cold dust, and it is now known that some cooling-flow cluster cD galaxies do as well. Yet this object, which is an extreme example of both, does not have enough cold gas to contribute significantly to the bolometric luminosity. We outline physical reasons why this could have happened. We then provide a discussion of strategies for finding hot dusty AGNs, given the limitations on submillimetre surveys implied by this work.     

The final two redshifts for radio sources from the equatorial BRL sample

Best, Röttgering & Lehnert (BRL) defined a new sample of powerful radio sources from the Molonglo Reference Catalogue, for which redshifts were compiled or measured for 177 of the 178 objects. For the final object, MRC1059−010 (3C 249), the host galaxy is identified here using near-infrared imaging, and the redshift is determined from Very Large Telescope (VLT) spectroscopy. For one other object in the sample, MRC0320+053 (4C05.14), the literature redshift has been questioned: new spectroscopic observations of this object are presented, deriving a corrected redshift. With these two results, the spectroscopic completeness of this sample is now 100 per cent.
New redshifts are also presented for PKS0742+10 from the Wall & Peacock 2.7-GHz catalogue, and for PKS1336+003 from the Parkes Selected Regions. PKS0742+10 shows a strong neutral hydrogen absorption feature in its Lyman α emission profile.