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

Millimetre continuum observations of southern massive star formation regions – II. SCUBA observations of cold cores and the dust grain emissivity index (β)

We report the results of a submillimetre continuum emission survey targeted towards 78 star formation regions, 72 of which are devoid of methanol maser and UC H  ii  regions, identified in the Swedish ESO Submillimetre Telescope (SEST)/SEST IMaging Bolometer Array (SIMBA) millimetre continuum survey of Hill et al. At least 45 per cent of the latter sources, dubbed 'mm-only', detected in this survey are also devoid of the mid-infrared MSX  emission. The 450- and 850-μm continuum emission was mapped using the Submillimetre Common User Bolometer Array (SCUBA) instrument on the James Clerk Maxwell Telescope (JCMT). Emission is detected towards 97 per cent of the 78 sources targeted as well as towards 28 other SIMBA sources lying in the SCUBA fields.
In total, we have identified 212 cores in this submillimetre survey, including 106 previously known from the SIMBA survey. Of the remaining 106 sources, 53 result from resolving a SIMBA source into multiple submillimetre components, whilst the other 53 sources are submillimetre cores, not seen in the SIMBA. Additionally, we have identified two further mm-only sources in the SIMBA images. Of the total 405 sources identified in the SIMBA survey, 255 are only seen at millimetre wavelengths.
We concatenate the results from four (sub)millimetre continuum surveys of massive star formation, together with the Galactic plane map of Pierce-Price et al. in order to determine the dust grain emissivity index β for each of the sources in the SIMBA source list. We examine the value of β with respect to temperature, as well as for the source classes identified in the SIMBA survey, for variation of this index. Our results indicate that β is typically 2, which is consistent with previous determinations in the literature, but for a considerably larger sample than previous work.     

Gravitational lensing and the Sunyaev–Zel'dovich effect in the millimetre/submillimetre waveband

The intensity of the cosmic microwave background radiation in the fields of clusters of galaxies is modified by inverse Compton scattering in the hot intracluster gas — the Sunyaev–Zel'dovich (SZ) effect. The effect is expected to be most pronounced at a frequency of about 350 GHz (a wavelength of about 800 μm), and has been detected in the centimetre and millimetre wavebands. In the millimetre/submillimetre waveband, the gravitationally lensed images of distant dusty star-forming galaxies in the background of the cluster are predicted to dominate the appearance of clusters on scales of several arcsec, and could confuse observations of the SZ effect at frequencies greater than about 200 GHz (wavelengths shorter than about 1.5 mm). Recent observations by Smail, Ivison &38; Blain confirm that a significant population of confusing sources is present in this waveband. Previous estimates of source confusion in observations of the millimetre/submillimetre-wave SZ effect did not include the effects of lensing by the cluster, and so the accuracy of such measurements could be lower than expected. Source subtraction may be required in order to measure the SZ effect accurately, and a careful analysis of the results of an ensemble of SZ measurements could be used to impose limits to the form of evolution of distant dusty star-forming galaxies.     

The temperatures of dust-enshrouded active galactic nuclei

A high density of massive dark objects (MDOs), probably supermassive black holes, in the centres of nearby galaxies has been inferred from recent observations. There are various indications that much of the accretion responsible for producing these objects took place in dust-enshrouded active galactic nuclei (AGNs). If so, then measurements of the intensity of background radiation and the source counts in the far-infrared and submillimetre wavebands constrain the temperature of dust in these AGNs. An additional constraint comes from the hard X-ray background, if this is produced by accretion. One possibility is that the dust shrouds surrounding the accreting AGNs are cold, about 30 K. In this event, the dusty AGNs could be some subset of the population of luminous distant sources discovered at 850 μm using the SCUBA array on the James Clerk Maxwell Telescope, as proposed by Almaini, Lawrence & Boyle. An alternative is that the dust shrouds surrounding the accreting AGNs are much hotter (>60 K). These values are closer to the dust temperatures of a number of well-studied low-redshift ultraluminous galaxies that are thought to derive their power from accretion. If the local MDO density is close to the maximum permitted, then cold sources cannot produce this density without the submillimetre background being overproduced if they accrete at high radiative efficiency, and thus a hot population is required. If the dust-enshrouded accretion occurred at similar redshifts to that taking place in unobscured optical quasars, then a significant fraction of the far-infrared background radiation measured by COBE at 140 μm, but very little of the submillimetre background at 850 μm, may have been produced by hot dust-enshrouded AGNs which may have already been seen in recent Chandra X-ray surveys.     

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.     

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.     

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.     

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.     

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.     

Submillimetre detection of a high-redshift type 2 QSO

We report on the first SCUBA detection of a type 2 QSO at   z = 3.660  in the Chandra Deep Field South. This source is X-ray-absorbed, shows only narrow emission lines in the optical spectrum and is detected in the submillimetre: it is the ideal candidate in an evolution scheme for active galactic nuclei (AGN) (e.g. Fabian 1999 ; Page et al. 2004 ) of an early phase corresponding to the main growth of the host galaxy and formation of the central black hole. The overall photometry (from the radio to the X-ray energy band) of this source is well reproduced by the spectral energy distribution (SED) of NGC 6240, while it is incompatible with the spectrum of a type 1 QSO (3C 273) or a starburst galaxy (Arp 220). Its submillimetre (850 μm) to X-ray (2 keV) spectral slope  (α_SX)  is close to the predicted value for a Compton-thick AGN in which only 1 per cent of the nuclear emission emerges through scattering. Using the observed flux at 850 μm we have derived a star formation rate of  550–680 M_⊙ yr⁻¹  and an estimate of the dust mass   M _dust= 4.2 × 10⁸ M_⊙  .     

Possible detection of V4334 Sgr (Sakurai's Object) at 450 and 850 μm

We report the possible detection of V4334 Sgr (Sakurai's Object) at 450 and 850 μm with SCUBA on the James Clerk Maxwell Telescope. The submillimetre photometry, combined with a  1–5 μm  spectrum and  8–10 μm  photometry obtained nearly contemporaneously, suggests that the submillimetre emission originates in material ejected during the 1995 event. The dust mass is a  few×10^-7 M_⊙  , the average mass-loss in the form of dust is  few×10^-8 M_⊙ yr^-1  , and the integrated luminosity is  log( L /L_⊙)=3.66  for a distance of 2 kpc. The ejected shell had angular diameter ∼55 mas in 2001 August, and should by now be resolvable in the mid-infrared by  8–10 m  class telescopes.     

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  .     

Spitzer IRAC infrared colours of submillimetre-bright galaxies

High-redshift submillimetre-bright galaxies identified by blank field surveys at millimetre and submillimetre wavelengths appear in the region of the Infra Red Array Camera (IRAC) colour–colour diagrams previously identified as the domain of luminous active galactic nuclei (AGNs). Our analysis using a set of empirical and theoretical dusty starburst spectral energy distribution (SED) models shows that power-law continuum sources associated with hot dust heated by young (≲100 Myr old), extreme starbursts at z > 2 also occupy the same general area as AGNs in the IRAC colour–colour plots. A detailed comparison of the IRAC colours and SEDs demonstrates that the two populations are distinct from each other, with submillimetre-bright galaxies having a systematically flatter IRAC spectrum (≳1 mag bluer in the observed [4.5]–[8.0] colour). Only about 20 per cent of the objects overlap in the colour–colour plots, and this low fraction suggests that submillimetre galaxies powered by a dust-obscured AGN are not common. The red infrared colours of the submillimetre galaxies are distinct from those of the ubiquitous foreground IRAC sources, and we propose a set of infrared colour selection criteria for identifying SMG counterparts that can be used even in the absence of radio or Multiband Imaging Photometer for Spitzer (MIPS) 24 μm data.     

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.     

Clarifying the nature of the brightest submillimetre sources: interferometric imaging of LH 850.02

We present high-resolution interferometric imaging of LH 850.02, the brightest 850- and 1200-μm submillimetre (submm) galaxy in the Lockman Hole. Our observations were made at 890 μm with the Submillimetre Array (SMA). Our high-resolution submm imaging detects LH 850.02 at  ≳6σ  as a single compact (size ≲1 arcsec or ≲8 kpc) point source and yields its absolute position to ∼0.2-arcsec accuracy. LH 850.02 has two alternative radio counterparts within the Submillimetre Common User Bolometer Array (SCUBA) beam (LH 850.02N and LH 850.02S), both of which are statistically very unlikely to be so close to the SCUBA source position by chance. However, the precise astrometry from the SMA shows that the submm emission arises entirely from LH 850.02N, and is not associated with LH 850.02S (by far the brighter of the two alternative identifications at 24 μm). Fits to the optical–infrared (IR) multicolour photometry of LH 850.02N and LH 850.02S indicate that both lie at   z ≈ 3.3  , and are therefore likely to be physically associated. At these redshifts, the 24-μm-to-submm flux density ratios suggest that LH 850.02N has an Arp 220-type starburst-dominated far-IR spectral energy distribution (SED), while LH 850.02S is more similar to Mrk 231, with less dust enshrouded star formation activity, but a significant contribution at 24 μm (rest frame 5–6 μm) from an active nucleus. This complex mix of star formation and active galactic nucleus (AGN) activity in multicomponent sources may be common in the high-redshift ultraluminous galaxy population, and highlights the need for precise astrometry from high-resolution interferometric imaging for a more complete understanding.     

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

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

Dense cores in the L1630 molecular cloud: discovering new protostars with SCUBA

Maps of the 450- and 850-μm dust continuum emission from three star-forming condensations within the Lynds 1630 molecular cloud, made with the SCUBA bolometer array, reveal the presence of four new submillimetre sources, each of a few solar masses (two of which are probably class I and two of which are class 0), as well as several sources the existence of which was previously known. The sources are located in filaments and appear elongated when observed at 450 μm. They probably have dust temperatures in the range 10 to 20 K, in good agreement with previous ammonia temperature estimates. Attempts to fit their structures with power-law and Gaussian density distributions suggest that the central distribution is flatter than expected for a simple singular isothermal sphere.
Although the statistics are poor, our results suggest that the ratio of 'protostellar core' mass to total virial mass may be similar for both large and small condensations.     

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.