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.     

Magnetic fields in massive star-forming regions

We present the largest sample of high-mass star-forming regions observed using submillimetre imaging polarimetry. The data were taken using the Submillimetre Common User Bolometer Array (SCUBA) in conjunction with the polarimeter on the James Clerk Maxwell Telescope (JCMT) in Hawaii. In total, 16 star-forming regions were observed, although some of these contain multiple cores. The polarimetry implies a variety of magnetic field morphologies, with some very ordered fields. We see a decrease in polarization percentage for seven of the cores. The magnetic field strengths estimated for 14 of the cores, using the corrected Chandrasekhar and Fermi (CF) method, range from <0.1 mG to almost 6 mG. These magnetic fields are weaker on these large scales when compared to previous Zeeman measurements from maser emission, implying the role of the magnetic field in star formation increases in importance on smaller scales. Analysis of the alignment of the mean field direction and the outflow directions reveals no relation for the whole sample, although direct comparison of the polarimetry maps suggests good alignment (to at least one outflow direction per source) in seven out of the 15 sources with outflows.     

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

SCUBA-2 CCD-STYLE IMAGING FOR THE JCMT

SCUBA-2 will replace SCUBA (Submillimetre Common User Bolometer Array) on the James Clerk Maxwell Telescope in 2006 and will be the first CCD-style camera for submillimetre astronomy. The instrument will simultaneously image at 850 and 450 microns using two focal plane arrays of 5120 pixels each. SCUBA-2 will map the submillimetre sky 1000 times faster than SCUBA to the same signal-to-noise ratio. This paper introduces the detector technology and the challenges faced in reading out a detector array cooled to ～120 mK.     

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.     

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.     

A submillimetre survey of the kinematics of the Perseus molecular cloud – I. Data

We present submillimetre observations of the   J = 3 → 2  rotational transition of ¹²CO, ¹³CO and C¹⁸O across over 600 arcmin² of the Perseus molecular cloud, undertaken with the Heterodyne Array Receiver Programme (HARP), a new array spectrograph on the James Clerk Maxwell Telescope. The data encompass four regions of the cloud, containing the largest clusters of dust continuum condensations: NGC 1333, IC348, L1448 and L1455. A new procedure to remove striping artefacts from the raw HARP data is introduced. We compare the maps to those of the dust continuum emission mapped with the Submillimetre Common-User Bolometer Array (SCUBA; Hatchell et al.) and the positions of starless and protostellar cores (Hatchell et al.). No straightforward correlation is found between the masses of each region derived from the HARP CO and SCUBA data, underlining the care that must be exercised when comparing masses of the same object derived from different tracers. From the ¹³CO/C¹⁸O line ratio the relative abundance of the two species  ([¹³CO]/[C¹⁸O]∼ 7)  and their opacities (typically τ is 0.02–0.22 and 0.15–1.52 for the C¹⁸O and ¹³CO gas, respectively) are calculated. C¹⁸O is optically thin nearly everywhere, increasing in opacity towards star-forming cores but not beyond  τ₁₈∼ 0.9  . Assuming the ¹²CO gas is optically thick, we compute its excitation temperature, T _ex (around 8–30 K), which has little correlation with estimates of the dust temperature.     

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.     

On the radio spectrum of CI Cygni: a test of popular models for symbiotic stars

We have determined the spectral energy distribution at wavelengths between 6 cm and 850 μm for the prototypical S(stellar)-type symbiotic star, CI Cygni, during quiescence. Data were obtained simultaneously with the Very Large Array and the SCUBA submillimetre (sub-mm) camera on the James Clerk Maxwell Telescope. The data have allowed us to determine the free–free turnover frequency of the ionized component, facilitating a model-dependent estimate of the binary separation to compare with the known orbital parameters of CI Cyg and to critically test the known models for radio emission from symbiotic stars. In particular, our data rule out the two most popular models: ionization of the giant wind by Lyman continuum photons from its hot companion, and emission resulting from the interaction of winds from the two binary components.     

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.     

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.     

A SCUBA survey of Orion – the low-mass end of the core mass function

We have re-analysed all of the Submillimetre Common User Bolometer Array (SCUBA) archive data of the Orion star-forming regions. We have put together all of the data taken at different times by different groups. Consequently, we have constructed the deepest submillimetre maps of these regions ever made. There are four regions that have been mapped: Orion A North and South, and Orion B North and South. We find that two of the regions, Orion A North and Orion B North, have deeper sensitivity and completeness limits, and contain a larger number of sources, so we concentrate on these two. We compare the data with archive data from the Spitzer Space Telescope to determine whether or not a core detected in the submillimetre is pre-stellar in nature. We extract all of the pre-stellar cores from the data and make a histogram of the core masses. This can be compared to the stellar initial mass function (IMF). We find the high-mass core mass function (CMF) follows a roughly Salpeter-like slope, just like the IMF, as seen in previous work. Our deeper maps allow us to see that the CMF turns over at,  ∼1.3 M_⊙  about a factor of 4 higher than our completeness limit. This turnover has never previously been observed, and is only visible here due to our much deeper maps. It mimics the turnover seen in the stellar IMF at  ∼0.1 M_⊙  . The low-mass side of the CMF is a power law with an exponent of, 0.35 ± 0.2 which is consistent with the low-mass slope of the young cluster IMF of 0.3 ± 0.1. This shows that the CMF continues to mimic the shape of the IMF all the way down to the lower completeness limit of these data at  ∼0.3 M_⊙  .     

A SCUBA survey of compact dark Lynds clouds

We present the first results of a submillimetre continuum survey of Lynds dark clouds. Submillimetre surveys of star-forming regions are an important tool with which to obtain representative samples of the very first phases of star formation. Maps of 24 small clouds were obtained with SCUBA, the bolometer array receiver at the James Clerk Maxwell Telescope, and 19 clouds were detected. The total dark cloud area surveyed was ∼130 arcmin², and a total gas mass of 90 M_⊙ was detected. The dust emission is in general in good agreement with the extinction of optical starlight. The observed clouds contain a newly discovered protostar in L944, and a previously known protostar IRAS 23228+4320 in L1246. Another eight starless cores, either gravitationally unbound or pre-stellar in nature, were also detected. All starless cores and protostars were detected in only seven clouds, and the remaining 17 clouds seem quiescent and do not show any signs of recent star formation activity. The 850-μm images of all detected clouds are presented, as well as 450-μm images of L328, L944, L1014 and L1262. The outflows of the protostars in L944 and L1246 were also discovered and were mapped in ¹²CO J =2→1. The detection of the young protostar in L944, which is not present in the IRAS Point Source Catalog, shows the capacity of submillimetre surveys to detect unknown protostars.     

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

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

The complete submillimetre spectrum of NGC 891

Submillimetre maps of NGC 891 have been obtained with the Programme National d'Observations Submillimétriques (PRONAOS) balloon-borne telescope and with the Infrared Space Observatory Photopolarimeter (ISOPHOT) on board the ISO satellite. In this article, we also gather data from the Infrared Astronomical Satellite ( IRAS ) and the Submillimetre Common-User Bolometer Array (SCUBA) to present the complete submillimetre spectrum of this nearby edge-on spiral galaxy. We derive submillimetre emission profiles along the major axis. The modified blackbody fits, assuming a single dust component, lead to temperatures of 19–24 K towards the centre and 18–20 K towards the edges, with possible variations in the dust spectral index from 1.4 to 2. The two-component fits lead to a warm-component temperature of 29 K all along the galaxy, with a cold component at 16 K. The interstellar medium (ISM) masses derived by these two methods are quite different:  4.6×10⁹ M_⊙  in the case of the one-component model and  12 × 10⁹ M_⊙  in the case of the two-component one. This two-component fit indicates that the cold dust to warm dust ratio is between 20 and 40, the highest values being in the wings of this galaxy. Compared to dust mass estimates, both estimates of the ISM mass are consistent with a gas to dust mass ratio of 240, which is close to the Milky Way value. Our results illustrate the importance of accurate submillimetre spectra in deriving masses of the ISM in galaxies.     

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 (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 JCMT Legacy Survey of the Gould Belt: a first look at Orion B with HARP

The Gould Belt Legacy Survey will survey nearby star-forming regions (within 500 pc), using Heterodyne Array Receiver Programme (HARP), Submillimetre Common-User Bolometer Array 2 and Polarimeter 2 on the James Clerk Maxwell Telescope. This paper describes the initial data obtained using HARP to observe ¹²CO, ¹³CO and C¹⁸O   J = 3 → 2  towards two regions in Orion B, NGC 2024 and NGC 2071. We describe the physical characteristics of the two clouds, calculating temperatures and opacities utilizing all the three isotopologues. We find good agreement between temperatures calculated from CO and from dust emission in the dense, energetic regions. We determine the mass and energetics of the clouds, and of the high-velocity material seen in ¹²CO emission, and compare the relative energetics of the high- and low-velocity material in the two clouds. We present a clumpfind analysis of the ¹³CO condensations. The slope of the condensation mass functions, at the high-mass ends, is similar to the slope of the initial mass function.     

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 burst of outflows from the Serpens cloud core: wide-field submillimetre continuum, CO J=2–1 and optical observations

Wide-field mapping of Serpens in submillimetre continuum emission and CO J =2–1 line emission is here complemented by optical imaging in [S  ii ] λλ 6716, 6731 line emission. Analysis of the 450- and 850-μm continuum data shows at least 10 separate sources, along with fainter diffuse background emission and filaments extending to the south and east of the core. These filaments describe 'cavity-like' structures that may have been shaped by the numerous outflows in the region. The dust opacity index, β , derived for the identifiable compact sources is of the order of 1.0±0.2, with dust temperatures in excess of 20 K. This value of β is somewhat lower than for typical class I YSOs; we suggest that the Serpens sources may be 'warm', late class 0 or early class I objects.
With the combined CO and optical data we also examine, on large scales, the outflows driven by the embedded sources in Serpens. In addition to a number of new Herbig–Haro flows (here denoted HH 455–460), a number of high-velocity CO lobes are observed; these extend radially outwards from the cluster of submillimetre sources in the core. A close association between the optical and molecular flows is also identified. The data suggest that many of the submillimetre sources power outflows. Collectively, the outflows traced in CO support the widely recognized correlation between source bolometric luminosity and outflow power, and imply a dynamical age for the whole protostellar cluster of ∼3×10⁴ yr. Notably, this is roughly equal to the proposed duration of the 'class 0' stage in protostellar evolution.