The SCUBA HAlf Degree Extragalactic Survey – VI. 350-μm mapping of submillimetre galaxies

A follow-up survey using the Submillimetre High-Angular Resolution Camera (SHARC-II) at 350 μm has been carried out to map the regions around several 850-μm-selected sources from the Submillimetre HAlf Degree Extragalactic Survey (SHADES). These observations probe the infrared (IR) luminosities and hence star formation rates in the largest existing, most robust sample of submillimetre galaxies (SMGs). We measure 350-μm flux densities for 24 850-μm sources, seven of which are detected at ≥2.5σ within a 10 arcsec search radius of the 850-μm positions. When results from the literature are included the total number of 350-μm flux density constraints of SHADES SMGs is 31, with 15 detections. We fit a modified blackbody to the far-IR (FIR) photometry of each SMG, and confirm that typical SMGs are dust-rich  ( M _dust≃ 9 × 10⁸ M_⊙)  , luminous  ( L _FIR≃ 2 × 10¹² L_⊙)  star-forming galaxies with intrinsic dust temperatures of ≃35 K and star formation rates of  ≃400 M_⊙ yr⁻¹  . We have measured the temperature distribution of SMGs and find that the underlying distribution is slightly broader than implied by the error bars, and that most SMGs are at 28 K with a few hotter. We also place new constraints on the 350-μm source counts, N ₃₅₀(>25 mJy) ∼ 200–500 deg⁻².     

A far-infrared view of the Lockman hole from ISO 95-μm observations – II. Optical identifications and insights into the nature of the far-infrared sources

We present the optical identifications of a 95-μm ISOPHOT sample in the Lockman hole over an area of approximately half a deg². The Rodighiero et al. catalogue includes 36 sources, making up a complete flux-limited sample for   S _95 μm≥ 100 mJy  . Reliable sources were detected, with decreasing but well-controlled completeness, down to   S _95 μm≃ 20 mJy  . We have combined mid-infrared (IR) and radio catalogues in this area to identify the potential optical counterparts of the far-IR sources. We found 14 radio and 13 15-μm associations, 10 of which have both associations. For the 11 sources with spectroscopic redshift, we have performed a spectrophotometric analysis of the observed spectral energy distributions (SEDs). Four of these 95-μm sources have been classified as faint IR (FIR) galaxies  ( L _FIR < 1. e 11 L_⊙)  , six as luminous IR galaxies (LIRGs) and only one as an ultraluminous IR galaxy (ULIRG). We have discussed the redshift distribution of these objects, comparing our results with evolutionary model predictions 95 and 175 μm. Given their moderate distances (the bulk of the closest spectroscopically identified objects lying at   z < 0.2  ), their luminosities and star formation rates (SFR; median value  ∼ 10 M_⊙ yr⁻¹  ), the sources unveiled by ISOPHOT at 95 μm seem to correspond to the low redshift  ( z < 0.3)  FIRBACK 175-μm population, composed of dusty, star-forming galaxies with moderate SFRs. We computed and compared different SFR estimators, and found that the SF derived from the bolometric IR luminosity is well correlated with that computed from the radio and mid-IR fluxes.     

The SCUBA Local Universe Galaxy Survey – II. 450-μm data: evidence for cold dust in bright IRAS galaxies

This is the second in a series of papers presenting results from the SCUBA Local Universe Galaxy Survey. In our first paper we provided 850-μm flux densities for 104 galaxies selected from the IRAS Bright Galaxy Sample and we found that the 60-, 100-μm ( IRAS ) and 850-μm (SCUBA) fluxes could be adequately fitted by emission from dust at a single temperature. In this paper we present 450-μm data for the galaxies. With the new data, the spectral energy distributions of the galaxies can no longer be fitted with an isothermal dust model – two temperature components are now required. Using our 450-μm data and fluxes from the literature, we find that the 450/850-μm flux ratio for the galaxies is remarkably constant, and this holds from objects in which the star formation rate is similar to our own Galaxy, to ultraluminous infrared galaxies (ULIRGs) such as Arp 220. The only possible explanation for this is if the dust emissivity index for all of the galaxies is ∼2 and the cold dust component has a similar temperature in all galaxies     . The 60-μm luminosities of the galaxies were found to depend on both the dust mass and the relative amount of energy in the warm component, with a tendency for the temperature effects to dominate at the highest L ₆₀. The dust masses estimated using the new temperatures are higher by a factor of ∼2 than those determined previously using a single temperature. This brings the gas-to-dust ratios of the IRAS galaxies into agreement with those of the Milky Way and other spiral galaxies which have been intensively studied in the submm.     

Implications for unified schemes from submillimetre and far-infrared follow-up of radio-selected samples

We extend our previous analysis which used generalized luminosity functions (GLFs) to predict the number of quasars and galaxies in low-radio-frequency-selected samples as a function of redshift, radio luminosity, narrow-emission-line luminosity and type of unified scheme. Our extended analysis incorporates the observed submillimetre (850-μm) flux densities of radio sources, employs a new method which allows us to deal with non-detections, and focuses on the high-luminosity population. First, we conclude that the submillimetre luminosity L ₈₅₀ of low-frequency-selected radio sources is correlated with the bolometric luminosity L _bol of their quasar nuclei via an approximate scaling relation   L ₈₅₀∝ L ^0.7±0.2_bol  . Secondly, we conclude that there is quantitative evidence for a receding-torus-like physical process for the high-luminosity population within a two-population unified scheme for radio sources; this evidence comes from the fact that radio quasars are brighter in both narrow emission lines and submillimetre luminosity than radio galaxies matched in radio luminosity and redshift. Thirdly, we note that the combination of a receding-torus-like scheme and the assumption that the observed submillimetre emission is dominated by quasar-heated dust yields a scaling relation   L ₈₅₀∝ L ^1/2_bol  which is within the errors of that determined here for radio-selected quasars, and consistent with that inferred for radio-quiet quasars.     

An ultradeep submillimetre map: beneath the SCUBA confusion limit with lensing and robust source extraction

Extracting sources with low signal-to-noise ratio (S/N) from maps with structured background is a non-trivial task which has become important in studying the faint end of the submillimetre (submm) number counts. In this paper, we study the source extraction from submm jiggle-maps from the Submillimetre Common-User Bolometer Array (SCUBA) using the Mexican hat wavelet (MHW), an isotropic wavelet technique. As a case study, we use a large (11.8-arcmin²) jiggle-map of the galaxy cluster Abell 2218 (A2218), with a 850-μm 1σ rms sensitivity of 0.6–1 mJy. We show via simulations that MHW is a powerful tool for the reliable extraction of low-S/N sources from the SCUBA jiggle-maps and nine sources are detected in the A2218 850-μm image. Three of these sources are identified as images of a single background source with an unlensed flux of 0.8 mJy. Further, two single-imaged sources also have unlensed fluxes <2 mJy, below the blank-field confusion limit. In this ultradeep map, the individual sources detected resolve nearly all of the extragalactic background light at 850 μm, and the deep data allow to put an upper limit of 44 sources arcmin⁻² to 0.2 mJy at 850 μm.     

The European Large Area ISO Survey – IV. The preliminary 90-μm luminosity function

We present the luminosity function of 90-μm-selected galaxies from the European Large Area ISO Survey (ELAIS), extending to z =0.3. Their luminosities are in the range 10⁹65⁻² L /L_⊙<10¹², i.e. non-ultraluminous. From our sample of 37 reliably detected galaxies in the ELAIS S1 region from the Efstathiou et al. S ₉₀100 mJy data base, we have found optical, 15-μm or 1.4-GHz identifications for 24 (65 per cent). We have obtained 2dF and UK Schmidt FLAIR spectroscopy of 89 per cent of identifications to rigid multivariate flux limits. We construct a luminosity function assuming that (i) our spectroscopic subset is an unbiased sparse sample, and (ii) there are no galaxies that would not be represented in our spectroscopic sample at any redshift. We argue that we can be confident of both assumptions. We find that the luminosity function is well described by the local 100-μm luminosity function of Rowan-Robinson, Helou & Walker. Assuming this local normalization, we derive luminosity evolution of (1+ z )^2.45±0.85 (95 per cent confidence). We argue that star formation dominates the bolometric luminosities of these galaxies, and we derive comoving star formation rates in broad agreement with the Flores et al. and Rowan-Robinson et al. mid-infrared-based estimates.     

Modelling the number counts of early-type galaxies by pure luminosity evolution

In this paper, we explore the plausible luminosity evolution of early-type galaxies in different cosmological models by constructing a set of pure luminosity evolution (PLE) models via the choices of the star-formation rate (SFR) parameters and formation redshift z _f of galaxies, with the observational constraints derived from the Hubble Space Telescope ( HST  ) morphological number counts for elliptical and S0 galaxies of the Medium Deep Survey (MDS) and the Hubble Deep Field (HDF). We find that the number counts of early-type galaxies can be explained by the pure luminosity evolution models, without invoking exotic scenarios such as merging or introducing an additional population, but the evolution should be nearly passive, with a high z _f assumed. The conclusion is valid in all of the three cosmological models we adopt in this paper. We also present the redshift distributions for three bins of observed magnitudes in the F814w passband, to show the redshift at which the objects that dominate the counts at a given magnitude may be found. The predictions of the redshift distribution of 22.5 <  b _j  < 24.0 are also presented for comparison with future data.     

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 SCUBA 8-mJy survey – II. Multiwavelength analysis of bright submillimetre sources

We present the results of a multiwavelength study of the 19 most significant submillimetre (submm) sources detected in the SCUBA 8-mJy survey. As described in Scott et al. , this survey covers ≃260 arcmin² using the submillimetre camera SCUBA, to a limiting source detection limit   S _850 μm≃8 mJy  . One advantage of this relatively bright flux-density limit is that accurate astrometric positions are potentially achievable for every source using existing radio and/or millimetre-wave interferometers. However, an associated advantage is that spectral energy distribution (SED) based redshift constraints should be more powerful than in fainter submm surveys. Here we therefore exploit the parallel SCUBA 450-μm data, in combination with existing radio and Infrared Space Observatory ( ISO ) data at longer and shorter wavelengths to set constraints on the redshift of each source. We also analyse new and existing optical and near-infrared imaging of our SCUBA survey fields to select potential identifications consistent with these constraints. Our derived SED-based redshift constraints, and the lack of statistically significant associations with even moderately bright galaxies allow us to conclude that all 19 sources lie at   z >1  , and at least half of them apparently lie at   z >2  .     

Hyperluminous infrared galaxies

39 galaxies are now known, from follow-up of faint IRAS sources and from submillimetre observations of high-redshift AGN, with far-infrared luminosities >10¹³ L_⊙. 13 of these, which have been found in 60- or 850-μm surveys, form an important unbiased subsample. 12 have been found by comparison of 60-μm surveys with quasar or radio galaxy catalogues, or from infrared surveys with colour selection biased towards AGN, while a further 14 have been found through submillimetre observations of known high-redshift AGN. In this paper I argue, on the basis of detailed modelling of the spectral energy distributions of hyperluminous galaxies with accurate radiative transfer models, and from evidence of high gas mass in several cases, that the bulk of the emission from these galaxies at rest frame wavelengths ≥50 μm is caused by star formation. Even after correction for the effects of lensing, hyperluminous galaxies with emission peaking at rest frame wavelengths ≥50 μm are therefore undergoing star formation at rates >10³ M_⊙ yr⁻¹ and are strong candidates for being primeval galaxies, in the process of a major episode of star formation.     

Probing the submillimetre number counts at f850 μm < 2 mJy

We have conducted a submillimetre mapping survey of faint, gravitationally lensed sources, where we have targeted 12 galaxy clusters and additionally the New Technology Telescope (NTT) Deep Field. The total area surveyed is 71.5 arcmin² in the image plane; correcting for gravitational lensing, the total area surveyed is 40 arcmin² in the source plane for a typical source redshift z ≈ 2.5. In the deepest maps, an image plane depth of 1σ rms ∼0.8 mJy is reached. This survey is the largest survey to date to reach such depths. In total 59 sources were detected, including three multiply imaged sources. The gravitational lensing makes it possible to detect sources with flux density below the blank field confusion limit. The lensing-corrected fluxes range from 0.11 to 19 mJy. After correcting for multiplicity, there are 10 sources with fluxes <2 mJy of which seven have submJy fluxes, doubling the number of such sources known. Number counts are determined below the confusion limit. At 1 mJy, the integrated number count is  ∼10⁴ deg⁻²  , and at 0.5 mJy it is  ∼2 × 10⁴ deg⁻²  . Based on the number counts, at a source plan flux limit of 0.1 mJy, essentially all of the 850-μm background emission has been resolved. The dominant contribution (>50 per cent) to the integrated background arises from sources with fluxes S ₈₅₀ between 0.4 and 2.5 mJy, while the bright sources S ₈₅₀ > 6 mJy contribute only 10 per cent.     

The K-band Hubble diagram of submillimetre galaxies and hyperluminous galaxies

We present the K -band Hubble diagrams ( K – z relations) of submillimetre-selected galaxies and hyperluminous galaxies (HLIRGs). We report the discovery of a remarkably tight K – z relation of HLIRGs, indistinguishable from that of the most luminous radio galaxies. Like radio galaxies, the HLIRG K – z relation at   z ≲ 3  is consistent with a passively evolving ∼3 L _* instantaneous starburst starting from a redshift of   z ∼ 10  . In contrast, many submillimetre-selected galaxies are ≳2 mag fainter, and the population has a much larger dispersion. We argue that dust obscuration and/or a larger mass range may be responsible for this scatter. The galaxies so far proved to be hyperluminous may have been biased towards higher AGN bolometric contributions than submillimetre-selected galaxies due to the 60-μm selection of some, so the location on the K – z relation may be related to the presence of the most massive active galactic nucleus. Alternatively, a particular host galaxy mass range may be responsible for both extreme star formation and the most massive active nuclei.     

A SCUBA scan-map of the Hubble Deep Field: measuring the bright end of the submillimetre source counts

Using the 850-μm SCUBA camera on the JCMT and a scanning technique different from other sub-mm surveys, we have obtained a 125 arcmin² map centred on the Hubble Deep Field . The 1 σ sensitivity to point sources is roughly 3 mJy and thus our map probes the brighter end of the sub-mm source counts. We find six sources with a flux greater than about 12 mJy (>4 σ ) and, after a careful accounting of incompleteness and flux bias, estimate the integrated density of bright sources     (68 per cent confidence bounds).     

Extended dust emission in NGC 7465

We present SCUBA 850-μm, JCMT  CO( J =2→1)  , B -band imaging and VLA H  i observations of the NGC 7465/4/3 group of galaxies. The 850-μm emission associated with NGC 7465 extends to at least ∼2 R ₂₅ and is well correlated with the H  i . We investigate a range of possible mechanisms by which dust beyond R ₂₅ may be heated to give the observed extended submillimetre emission. By modelling the dust heating by stars in two extreme geometries, we fail to find any reasonable star formation scenario that is consistent with both the 850-μm and optical data. Furthermore, we do not detect any  CO( J =2→1)  emission coincident with the extended dust and atomic gas as would be expected if significant star formation were occurring. We show that shock-heating of dust via cloud–cloud collisions in the stripped interstellar medium of NGC 7465 could be sufficient to explain the extended 850-μm emission and lack of optical emission in the stripped gas, and suggest that cloud–cloud collisions may be an important dust heating mechanism in gas-rich systems.     

Measurement of the Sunyaev–Zel'dovich increment in massive galaxy clusters

We have detected the Sunyaev–Zel'dovich (SZ) increment at 850 μm in two galaxy clusters (Cl 0016+16 and MS 1054.4−0321) using the Submillimetre Common User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope. Fits to the isothermal β model yield a central Compton y parameter of  (2.2 ± 0.7) × 10⁻⁴  and a central 850-μm flux of  Δ I ₀= 2.2 ± 0.7 mJy beam⁻¹  in Cl 0016. This can be combined with decrement measurements to infer   y = (2.38 ±^0.36_0.34) × 10⁻⁴  and   v _pec= 400±¹⁹⁰⁰₁₄₀₀ km s⁻¹  . In MS 1054 we find a peak 850-μm flux of  Δ I ₀= 2.0 ± 1.0 mJy beam⁻¹  and   y = (2.0 ± 1.0) × 10⁻⁴  . To be successful such measurements require large chop throws and non-standard data analysis techniques. In particular, the 450-μm data are used to remove atmospheric variations in the 850-μm data. An explicit annular model is fit to the SCUBA difference data in order to extract the radial profile, and separately fit to the model differences to minimize the effect of correlations induced by our scanning strategy. We have demonstrated that with sufficient care, SCUBA can be used to measure the SZ increment in massive, compact galaxy clusters.     

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 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 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.     

Dust and hydrogen molecules in the extremely metal-poor dwarf galaxy SBS 0335–052

During the early stages of galaxy evolution, the metallicity is generally low and nearby metal-poor star-forming galaxies may provide templates for primordial star formation. In particular, the dust content of such objects is of great importance, because early molecular formation can take place on grains. To gain insight into primeval galaxies at high redshift, we examine the dust content of the nearby extremely low-metallicity galaxy SBS  0335–052  which hosts a very young starburst (≲10⁷ yr). In young galaxies, the dust formation rate in Type II supernovae governs the amount of dust, and by incorporating recent results on dust production in Type II supernovae we model the evolution of dust content. If the star-forming region is compact (≲100 pc), as suggested by observations of SBS  0335–052  , our models consistently explain the quantity of dust, far-infrared luminosity, and dust temperature in this low-metallicity object. We also discuss the H₂ abundance. The compactness of the region is important to H₂ formation, because the optical depth of dust for UV photons becomes large and H₂ dissociation is suppressed. We finally focus on implications for damped Ly α systems.     

Luminous early-type field galaxies at z∼0.4– I. Observations and redshift catalogue of 581 galaxies

We have compiled a sample of ∼ 9600 bright, i ≤18.95 , red, b _J− r >2 , candidate galaxies in an area of 220 deg². These are luminous, L > L * , field early-type galaxies with redshifts 0.3≲ z ≲0.6 . We present a redshift catalogue of a subsample of 581 targets. The galaxies were selected according to their broad-band b _J ri colours from United Kingdom Schmidt Telescope plates, and have a surface density on the sky of only ∼ 50 deg⁻². Such luminous field galaxies are virtually absent from published redshift surveys and the catalogue provides a large sample of the most luminous normal galaxies, at cosmological distances. The statistical properties of the galaxy spectra, including absorption-line and emission-line measures, are presented and a composite spectrum constructed. The nature of the sample, combined with the relatively bright apparent magnitudes, makes the galaxies suitable targets for several key investigations in galaxy evolution and cosmology.