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

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) – VIII. The nature of faint submillimetre galaxies in SHADES, SWIRE and SXDF surveys

We present the optical-to-submillimetre spectral energy distributions (SEDs) for 33 radio and mid-infrared (mid-IR) identified submillimetre galaxies discovered via the SHADES 850-μm SCUBA imaging in the Subaru- XMM Deep Field (SXDF). Optical data for the sources come from the SXDF and mid- and far-IR fluxes from SWIRE. We obtain photometric redshift estimates for our sources using optical and IRAC 3.6- and 4.5-μm fluxes. We then fit SED templates to the longer wavelength data to determine the nature of the far-IR emission that dominates the bolometric luminosity of these sources. The IR template fits are also used to resolve ambiguous identifications and cases of redshift aliasing. The redshift distribution obtained broadly matches previous results for submillimetre sources and on the SHADES SXDF field. Our template fitting finds that active galactic nuclei, while present in about 10 per cent of our sources, do not contribute significantly to their bolometric luminosity. Dust heating by starbursts, with either Arp220 or M82 type SEDs, appears to be responsible for the luminosity in most sources (23/33 are fitted by Arp220 templates, 2/33 by the warmer M82 templates). 8/33 sources, in contrast, are fitted by a cooler cirrus dust template, suggesting that cold dust has a role in some of these highly luminous objects. Three of our sources appear to have multiple identifications or components at the same redshift, but we find no statistical evidence that close associations are common among our SHADES sources. Examination of rest-frame K -band luminosity suggests that 'downsizing' is underway in the submillimetre galaxy population, with lower redshift systems lying in lower mass host galaxies. Of our 33 identifications six are found to be of lower reliability but their exclusion would not significantly alter our conclusions.     

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.     

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  .     

Hubble Space Telescope imaging survey of sub-mJy star-forming galaxies – I. Morphologies at z ∼ 0.2

We present the first results of our Hubble Space Telescope HST WFPC2 F814W snapshot imaging survey, targeting virtually all sub-mJy decimetric radio-selected star-forming galaxies. The radio selection at ∼1 GHz is free from extinction effects and the radio luminosities are largely unaffected by AGN contamination, making these galaxies ideal tracers of the cosmic star formation history. A subsample of four targets is presented here, selected at 1.4 GHz from the spectroscopically homogenous and complete samples of Benn et al. and Hopkins et al. The redshifts are confined to a narrow range around z ∼0.2, to avoid differential evolution, with a radio luminosity close to L ∗ where the galaxies dominate the comoving volume-averaged star formation rate. We find clearly disturbed morphologies resembling those of ultraluminous infrared galaxies, indicating that galaxy interactions may be the dominant mechanism for triggering star formation at these epochs. The morphologies are also clearly different from those of coeval quasars and radio galaxies, as found in star-forming galaxies selected at other wavelengths. This may prove challenging for models that propose direct causal links between AGN evolution and the cosmic star formation history at these epochs. The asymmetries are typically much larger than seen in the Canada–France Redshift Survey at similar redshifts, optical luminosities and H α -derived star formation rates, indicating the possible existence of an obscuration-related morphological bias in such samples.