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

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.     

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.     

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  .     

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

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

SCUBA observations of the sources detected in the MAMBO 1200-μm survey

We have observed 23 sources from the Max-Planck Millimetre Bolometer (MAMBO) array 1200-μm survey with SCUBA at 850 μm, detecting 19 of the sources. The sources generally have low values for the ratio of 850- to 1200-μm flux. Two possible explanations for the low values are either that the sources are at very high redshifts or that the global properties of the dust in the MAMBO sources are different from the global properties of dust in low-redshift galaxies. If the former explanation is correct, we estimate that 15 of the MAMBO sources lie at   z > 3  .     

The European Large Area ISO Survey – I. Goals, definition and observations

We describe the European Large Area ISO Survey (ELAIS). ELAIS was the largest single Open Time project conducted by ISO , mapping an area of 12 deg² at 15 μm with ISOCAM and at 90 μm with ISOPHOT. Secondary surveys in other ISO bands were undertaken by the ELAIS team within the fields of the primary survey, with 6 deg² being covered at 6.7 μm and 1 deg² at 175 μm.
This paper discusses the goals of the project and the techniques employed in its construction, as well as presenting details of the observations carried out, the data from which are now in the public domain. We outline the ELAIS 'preliminary analysis' which led to the detection of over 1000 sources from the 15 and 90-μm surveys (the majority selected at 15 μm with a flux limit of ∼3 mJy), to be fed into a ground-based follow-up campaign, as well as a programme of photometric observations of detected sources using both ISOCAM and ISOPHOT.
We detail how the ELAIS survey complements other ISO surveys in terms of depth and areal coverage, and show that the extensive multi-wavelength coverage of the ELAIS fields resulting from our concerted and on-going follow-up programme has made these regions amongst the best studied areas of their size in the entire sky, and, therefore, natural targets for future surveys. This paper accompanies the release of extremely reliable subsets of the 'preliminary analysis' products. Subsequent papers in this series will give further details of our data reduction techniques, reliability and completeness estimates and present the 15‐ and 90-μm number counts from the 'preliminary analysis', while a further series of papers will discuss in detail the results from the ELAIS 'final analysis', as well as from the follow-up programme.     

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.     

Hot dust in two hard Chandra X-ray sources

The two brightest hard X-ray sources discovered serendipitously by Chandra in the field of the lensing cluster A2390 are found to have ISOCAM counterparts at 6.7 and 15 μm. We use this fact, together with their non-detection by SCUBA at 850 μm, as the basis for dusty radiative transfer modelling of their infrared spectral energy distributions. For the best-fitting models, we find that the dust that reprocesses the optical–ultraviolet light in these Compton-thin active galactic nuclei (AGN) is heated to near its sublimation temperature (above 1000 K), with an inner radius within a parsec of the nucleus. Some warm-dust models with inner temperatures of 200 K are also formally acceptable. These findings strongly support the obscured AGN hypothesis for the new hard X-ray Chandra sources, which lack both strong emission lines and 850-μm SCUBA detections.     

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

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

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.     

Radio source stacking and the infrared/radio correlation at μJy flux densities

We investigate the infrared/radio correlation using the technique of source stacking, in order to probe the average properties of radio sources that are too faint to be detected individually. We compare the two methods used in the literature to stack sources and demonstrate that the creation of stacked images leads to a loss of information. We stack infrared sources in the Spitzer Extragalactic First Look Survey (xFLS) field, and the three northern Spitzer Wide-area Infrared Extragalactic survey (SWIRE) fields, using radio surveys created at 610 MHz and 1.4 GHz, and find a variation in the absolute strength of the correlation between the xFLS and SWIRE regions, but no evidence for significant evolution in the correlation over the 24-μm flux density range 150 μJy to 2 mJy. We carry out the first radio source stacking experiment using 70-μm-selected galaxies, and find no evidence for significant evolution over the 70-μm flux density range 10–100 mJy.     

RXTE observation of NGC 6240: a search for the obscured active nucleus

We present new data taken at 850 μm with SCUBA at the James Clerk Maxwell Telescope for a sample of 19 luminous infrared galaxies. Fourteen galaxies were detected. We have used these data, together with fluxes at 25, 60 and 100 μm from IRAS , to model the dust emission. We find that the emission from most galaxies can be described by an optically thin, single temperature dust model with an exponent of the dust extinction coefficient ( k _λ ∝ λ ^{− β} ) of β ≃1.4–2. A lower β ≃1 is required to model the dust emission from two of the galaxies, Arp 220 and NGC 4418. We discuss various possibilities for this difference and conclude that the most likely is a high dust opacity. In addition, we compare the molecular gas mass derived from the dust emission, M _850 μm, with the molecular gas mass derived from the CO emission, M _CO, and find that M _CO is on average a factor 2–3 higher than M _850 μm.     

Gravitational lensing of extended high-redshift sources by dark matter haloes

High-redshift galaxies and quasi-stellar objects (QSOs) are most likely to be strongly lensed by intervening haloes between the source and the observer. In addition, a large fraction of lensed sources is expected to be seen in the submillimetre region, as a result of the enhanced magnification bias on the steep intrinsic number counts. We extend in three directions Blain's earlier study of this effect.
First, we use a modification of the Press–Schechter mass function and detailed lens models to compute the magnification probability distribution. We compare the magnification cross-sections of populations of singular isothermal spheres and Navarro, Frenk & White (NFW) haloes and find that they are very similar, in contrast to the image-splitting statistics which were recently investigated in other studies. The distinction between the two types of density profile is therefore irrelevant for our purposes.
Secondly, we discuss quantitatively the maximum magnification, μ _max, that can be achieved for extended sources (galaxies) with realistic luminosity profiles, taking into account the possible ellipticity of the lensing potential. We find that μ _max plausibly falls into the range     for sources of     effective radius at redshifts within     .
Thirdly, we apply our model for the lensing magnification to a class of sources following the luminosity evolution typical for a unified scheme of QSO formation. As a result of the peculiar steepness of their intrinsic number counts, we find that the lensed source counts at a fiducial wave length of 850 μm can exceed the unlensed counts by several orders of magnitude at flux densities ≳100 mJy, even with a conservative choice of the maximum magnification.