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.     

Submillimetre detection of the z=2.83 Lyman-break galaxy, Westphal–MM8, and implications for SCUBA2

We present confusion-limited submillimetre (submm) observations with the Submillimetre Common-User Bolometer Array (SCUBA) camera on the James Clerk Maxwell Telescope of the   z = 2.83  Lyman-break galaxy (LBG), Westphal–MM8, reaching an 850 μm sensitivity even greater than that achieved in the SCUBA map of the Hubble Deep Field region. The detection of MM8  ( S _850 μm= 1.98 ± 0.48 mJy)  , along with the literature submm detections of lensed LBGs, suggests that the LBG population may contribute significantly to the source counts of submm-selected galaxies in the 1–2 mJy regime. Additionally, submm-luminous LBGs are a viable progenitor population for the recently discovered evolved galaxies at   z ∼ 2–2.5  . These observations represent an important baseline for SCUBA2 observations which will regularly map large regions of the sky to this depth.     

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

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.     

SCUBA imaging of the NGC 7331 dust ring

We present observations of the spiral galaxy NGC 7331 using the Submillimetre Common User Bolometer Array (SCUBA) on the James Clark Maxwell Telescope. We have detected a dust ring of 45 arcsec radius (3.3 kpc) at wavelengths of 450 and 850 μm. The dust ring is in good correspondence with other observations of the ring in the mid-infrared (MIR), CO and radio continuum, suggesting that the observed dust is associated with molecular gas and star formation. A B  −  K colour map shows an analogous ring structure with an asymmetry about the major axis, consistent with the extinction being produced by a dust ring. The derived temperature of the dust lies between 16 and 31 K and the gas-to-dust ratio lies between 150 and 570, depending on the assumed dust emission efficiency index (β = 1.5 or 2).     

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.     

Confirming a population of hot-dust dominated, star-forming, ultraluminous galaxies at high redshift

We identify eight   z > 1  radio sources undetected at 850 μm but robustly detected at 70 μm, confirming that they represent ultraluminous infrared galaxies (ULIRGs) with hotter dust temperatures  (〈 T _d〉= 52 ± 10 K)  than submillimetre galaxies (SMGs) at similar luminosities and redshifts. These galaxies share many properties with SMGs: ultraviolet spectra consistent with starbursts, high stellar masses and radio luminosities. We can attribute their radio emission to star formation since high-resolution Multi-Element Radio Linked Interferometer Network (MERLIN) radio maps show extended emission regions (with characteristic radii of 2–3 kpc), which are unlikely to be generated by active galactic nucleus (AGN) activity. These observations provide the first direct confirmation of hot, dusty ULIRGs which are missed by current submillimetre surveys. They have significant implications for future observations from the Herschel Space Observatory and Submillimetre Common-User Bolometer Array 2 (SCUBA2), which will select high-redshift luminous galaxies with less selection biases.     

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  .     

A study of the Sunyaev–Zel'dovich increment using archival SCUBA data

In a search for evidence of the short wavelength increment in the Sunyaev–Zel'dovich (SZ) effect, we have analysed archival galaxy cluster data from the Submillimetre Common User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope, resulting in the most complete pointed survey of clusters at 850 μm to date. SCUBA's 850-μm passband overlaps the peak of the SZ increment. The sample consists of 44 galaxy clusters in the range 0 < z < 1.3. Maps of each of the clusters have been made and sources have been extracted; as an ancillary product, we generate the most thorough galaxy cluster point source list yet from SCUBA. 17 of these clusters are free of obvious active galactic nuclei (AGN) and have data deep enough to provide interesting measurements of the expected SZ signal. Specialized analysis techniques are employed to extract the SZ effect signal from these SCUBA data, including using SCUBA's short wavelength band as an atmospheric monitor and fitting the long wavelength channel to a model of the spatial distribution of each cluster's SZ effect. By explicitly excising the exact cluster centre from our analysis, we demonstrate that emission from galaxies within the cluster does not contaminate our measurement. The SZ amplitudes from our measurements are consistently higher than the amplitudes inferred from low-frequency measurements of the SZ decrement.     

Evolutionary constraints from infrared source counts

The backward evolution approach to modelling galaxy source counts is re-visited in the wake of the numerous results and revelations from the Infrared Space Observatory ( ISO ), the Submillimetre Common User Bolometer Array (SCUBA) and the detections and measurements of the cosmic extragalactic background light. Using the framework of the Pearson & Rowan-Robinson galaxy evolution model, the observed source counts and background measurements are used to constrain the evolution in the galaxy population. It is found that a strong evolution in both density and luminosity of the high-luminosity tail of the infrared (IR) luminosity function, interpreted as the ultraluminous galaxies discovered first by IRAS and later elevated in status by SCUBA and ISO , can account for the source counts from 15 μm (where it matches the undulations in the integral counts and the hump in the differential counts extremely well) to the submillimetre region, as well as explain the peak in the cosmic infrared background at ∼140 μm. The submillimetre counts are interpreted as the superposition of two separate populations comprising ultraluminous infrared galaxies (ULIGs) at the brighter submillimetre fluxes and starburst galaxies at fluxes fainter than ∼2 mJy. In this scenario the high-redshift ULIGs are tenuously interpreted as the progenitors of today's giant elliptical (gE) galaxies.
All the source count models can be accessed via the world wide web at the URL http://www.ir.isas.ac.jp/~cpp/counts/     

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

Submillimetre imaging of NGC 3079

Submillimetre mapping observations of the active edge-on spiral galaxy NGC 3079 are presented. These maps at 850 and 450 μm were made with the Submillimetre Common User Bolometer Array (SCUBA) at the James Clerk Maxwell Telescope (JCMT).
The source structure at these wavelengths consists of a central unresolved source embedded in diffuse disc emission, similar to that displayed at 1.2 mm. The disc emission is fitted with two optically thin, isothermal dust models which give temperatures of 12 and 31 K, similar to those derived previously by Braine et al. The core component is well described by a single-temperature fit (∼32 K). The combined dust mass from these observations, using the same mass absorption coefficient as Devereux & Young (1990) is 3.5×10⁸ M_⊙, of which ∼90 per cent resides in the cold component of the galactic disc. The effect of the cold dust component detected by SCUBA is thus to reduce the global gas-to-dust mass ratio from ∼1400 found in the above study to 85, very similar to the Galactic level. Calculations using the models of Draine & Lee and/or alternative molecular gas mass estimates yield gas-to-dust mass ratios in the range 60–190.
The data presented here, together with previously published 1.2-mm mapping observations and IRAS data, are inconsistent with detections made with the Infrared Space Observatory ( ISO ). In particular, the latter give an excess of flux at 200 and 180 μm relative to that predicted by our simple model fits (approximately a factor of 2–3).     

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 effect of lensing on the identification of SCUBA galaxies

Spectroscopic surveys of submillimetre(submm)-selected sources have uncovered optically bright galaxies at   z ≲ 1  close to the positions of several submm sources. Naive statistical analyses suggest that these galaxies are associated with the submm emission. However, in some cases, it is difficult to understand this association given the relatively modest redshifts and unprepossessing spectral characteristics of the galaxies. These are in stark contrast to those expected from the massive dust-enshrouded starbursts and AGN thought to power the bulk of the bright submm population. We present new observations of optically bright counterparts to two luminous submm sources, along with a compilation of previously proposed optically bright counterparts with   z ≲ 1  . We suggest that the majority of these associations between bright galaxies and submm sources may be as a result of the action of the foreground galaxies as gravitational lenses on the much fainter and more distant submm sources. We discuss the implications of this conclusion for our understanding of the SCUBA population.     

Observational limits to source confusion in the millimetre/submillimetre waveband

The first observations to detect a population of distant galaxies directly in the submillimetre waveband have recently been made using the new Submillimetre Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT). The results indicate that a large number of distant galaxies are radiating strongly in this waveband. Here we discuss their significance for source confusion in future millimetre/submillimetre-wave observations of both distant galaxies and cosmic microwave background radiation (CMBR) anisotropies. Earlier estimates of such confusion involved significant extrapolation of the results of observations of galaxies at low redshifts; our new estimates do not, as they are derived from direct observations of distant galaxies in the submillimetre waveband. The results have important consequences for the design and operation of existing and proposed millimetre/submillimetre-wave telescopes: the Planck Surveyor survey will be confusion-limited at frequencies greater than 350 GHz, even in the absence of Galactic dust emission; a 1σ confusion noise limit of about 0.44 mJy beam⁻¹ is expected for the JCMT/SCUBA at a wavelength of 850 μm; and the subarcsecond resolution of large millimetre/submillimetre-wave interferometer arrays will be required in order to execute very deep galaxy surveys.     

The evolution of star formation in quasar host galaxies

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

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

Millimetric properties of gamma-ray burst host galaxies

We present millimetre (mm) and submillimetre (submm) photometry of a sample of five host galaxies of gamma-ray bursts (GRBs), obtained using the Max Planck Millimetre Bolometer (MAMBO2) array and Submillimetre Common-User Bolometer Array (SCUBA). These observations were obtained as part of an ongoing project to investigate the status of GRBs as indicators of star formation. Our targets include two of the most unusual GRB host galaxies, selected as likely candidate submm galaxies: the extremely red  ( R − K ≈ 5)  host of GRB 030115, and the extremely faint  ( R > 29.5)  host of GRB 020124. Neither of these galaxies is detected, but the deep upper limits for GRB 030115 impose constraints on its spectral energy distribution, requiring a warmer dust temperature than is commonly adopted for submillimetre galaxies (SMGs).
As a framework for interpreting these data, and for predicting the results of forthcoming submm surveys of Swift -derived host samples, we model the expected flux and redshift distributions based on luminosity functions of both submm galaxies and GRBs, assuming a direct proportionality between the GRB rate density and the global star formation rate density. We derive the effects of possible sources of uncertainty in these assumptions, including (1) introducing an anticorrelation between GRB rate and the global average metallicity, and (2) varying the dust temperature.     

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.