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.     

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

Submillimetre continuum images of the NGC 2024 star-forming ridge

We present 450- and 800-μm images, made with the James Clerk Maxwell Telescope, of the NGC 2024 molecular ridge. The seven previously known compact cores, FIR1–7, have been detected, and FIR5 has been resolved into a compact object and an associated extended source to the east. The estimated masses of the dense cores vary between 1.6 and 5.1 M⊙ per 14-arcsec beam, assuming a dust temperature of 30 K and a dust opacity of κ_800 μm = 0.002 m² kg⁻¹. A spectral index map made from the 450- and 800-μm images shows spatial variations, with the spectral index, α ( F _ν ∝ ν^α), being systematically lower towards the dense cores. We interpret this as evidence for a lower value of the frequency dependence of the dust opacity, β, towards the denser cores relative to the surrounding molecular material. This may indicate that grain growth is occurring in the cores, prior to planetesimal formation. By comparing the high-resolution 450-μm image with interferometer maps of the integrated CS(2–1) emission, the previously reported discrepancy between dust continuum emission and molecular line emission is found to be very localized. Depletion and temperature variations are discussed as possible explanations.     

A shallow though extensive H2 2.122-μm imaging survey of Taurus–Auriga–Perseus – I. NGC 1333, L1455, L1448 and B1

We discuss wide-field near-infrared (near-IR) imaging of the NGC 1333, L1448, L1455 and B1 star-forming regions in Perseus. The observations have been extracted from a much larger narrow-band imaging survey of the Taurus–Auriga–Perseus complex. These H₂ 2.122-μm observations are complemented by broad-band K imaging, mid-IR imaging and photometry from the Spitzer Space Telescope , and published submillimetre CO   J = 3–2  maps of high-velocity molecular outflows. We detect and label 85 H₂ features and associate these with 26 molecular outflows. Three are parsec-scale flows, with a mean flow lobe length exceeding 11.5 arcmin. 37 (44 per cent) of the detected H₂ features are associated with a known Herbig–Haro object, while 72 (46 per cent) of catalogued HH objects are detected in H₂ emission. Embedded Spitzer sources are identified for all but two of the 26 molecular outflows. These candidate outflow sources all have high near-to-mid-IR spectral indices (mean value of  α∼ 1.4  ) as well as red IRAC 3.6–4.5 μm and IRAC/MIPS 4.5–24.0 μm colours: 80 per cent have [3.6]–[4.5] > 1.0 and [4.5]–[24] > 1.5. These criteria – high α and red [4.5]–[24] and [3.6]–[4.5] colours – are powerful discriminants when searching for molecular outflow sources. However, we find no correlation between α and flow length or opening angle, and the outflows appear randomly orientated in each region. The more massive clouds are associated with a greater number of outflows, which suggests that the star formation efficiency is roughly the same in each region.     

ULTRACAM observations of the black hole X-ray transient XTE J1118+480 in quiescence

We present 450- and 850-μm maps of R Coronae Australis. We compare the maps with previous surveys of the region, and shed new light on the previously unknown nature of the protostellar sources at the centre of the cloud. We clarify the nature of two millimetre sources previously discovered in lower-resolution data. We identify one new Class 0 protostar, which we label SMM 1B, and we measure the envelope masses of a number of more evolved protostars. We identify two new pre-stellar cores, which we call SMM 1A and SMM 6.     

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.     

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

Infrared mapping of the dust disc around Vega

ISOPHOT has been used to perform high-resolution 60-μm scans of Vega, and these have been compared with those from γDra, to obtain a Gaussian width of 22±2 arcsec. The dust disc around Vega has been mapped, resolving it at 60 and 90 μm with ISOPHOT. At 90 μm a Gaussian width of 36±3 arsec has been derived. In addition, multi-filter photometry is presented, at 25, 60, 80, 100, 120, 150, 170 and 200 μm. The data are fitted by a modified blackbody with a temperature of 73 K [ Q (λ)∝1/λ^1.1]. The dust disc has a luminosity L _IR/ L _*∼3×10⁻⁵. Using a distance of 7.8 pc, 22 arcsec corresponds to a distance of 86 au, and 36 arcsec to a distance of 140 au.     

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.     

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.     

The European Large Area ISO Survey – IV. The preliminary 90-μm luminosity function

We present the luminosity function of 90-μm-selected galaxies from the European Large Area ISO Survey (ELAIS), extending to z =0.3. Their luminosities are in the range 10⁹65⁻² L /L_⊙<10¹², i.e. non-ultraluminous. From our sample of 37 reliably detected galaxies in the ELAIS S1 region from the Efstathiou et al. S ₉₀100 mJy data base, we have found optical, 15-μm or 1.4-GHz identifications for 24 (65 per cent). We have obtained 2dF and UK Schmidt FLAIR spectroscopy of 89 per cent of identifications to rigid multivariate flux limits. We construct a luminosity function assuming that (i) our spectroscopic subset is an unbiased sparse sample, and (ii) there are no galaxies that would not be represented in our spectroscopic sample at any redshift. We argue that we can be confident of both assumptions. We find that the luminosity function is well described by the local 100-μm luminosity function of Rowan-Robinson, Helou & Walker. Assuming this local normalization, we derive luminosity evolution of (1+ z )^2.45±0.85 (95 per cent confidence). We argue that star formation dominates the bolometric luminosities of these galaxies, and we derive comoving star formation rates in broad agreement with the Flores et al. and Rowan-Robinson et al. mid-infrared-based estimates.     

Spitzer and ISO galaxy counts in the mid-infrared

Galaxy source counts that simultaneously fit the deep mid-infrared surveys at 24 microns and 15 microns made by the Spitzer Space Telescope and the Infrared Space Observatory ( ISO ), respectively, are presented for two phenomenological models. The models are based on starburst and luminous infrared galaxy dominated populations. Both models produce excellent fits to the counts in both wavebands and provide an explanation for the high-redshift population seen in the longer Spitzer 24-micron band supporting the hypothesis that they are luminous–ultraluminous infrared galaxies at   z = 2–3  , being the mid-infrared counterparts to the submillimetre galaxy population. The source counts are characterized by strong evolution to redshift unity, followed by less drastic evolution to higher redshift. The number–redshift distributions in both wavebands are well explained by the effect of the many mid-infrared features passing through the observation windows. The sharp upturn at around a millijansky in the 15-μm counts in particular depends critically on the distribution of mid-infrared features around 12 μm, in the assumed spectral energy distribution.     

Spitzer observations of acetylene bands in carbon-rich asymptotic giant branch stars in the Large Magellanic Cloud

We investigate the molecular bands in carbon-rich asymptotic giant branch (AGB) stars in the Large Magellanic Cloud (LMC), using the Infrared Spectrograph (IRS) onboard the Spitzer Space Telescope ( SST ) over the 5–38 μm range. All 26 low-resolution spectra show acetylene (C₂H₂) bands at 7 and 14 μm. The hydrogen cyanide (HCN) bands at these wavelengths are very weak or absent. This is consistent with low nitrogen abundances in the LMC. The observed 14 μm C₂H₂  band is reasonably reproduced by an excitation temperature of 500 K. There is no clear dilution of the 14 μm C₂H₂  band by circumstellar dust emission. This 14-μm band originates from molecular gas in the circumstellar envelope in these high mass-loss rate stars, in agreement with previous findings for Galactic stars. The C₂H₂ column density, derived from the 13.7 μm band, shows a gas mass-loss rate in the range 3 × 10⁻⁶ to 5 × 10⁻⁵ M_⊙ yr⁻¹. This is comparable with the total mass-loss rate of these stars estimated from the spectral energy distribution. Additionally, we compare the line strengths of the 13.7 μm C₂H₂  band of our LMC sample with those of a Galactic sample. Despite the low metallicity of the LMC, there is no clear difference in the C₂H₂  abundance among LMC and Galactic stars. This reflects the effect of the third dredge-up bringing self-produced carbon to the surface, leading to high carbon-to-oxygen ratio at low metallicity.     

The initial conditions of isolated star formation – V. ISOPHOT imaging and the temperature and energy balance of pre-stellar cores

ISO data taken with the long-wavelength imaging photo-polarimeter ISOPHOT are presented of 18 pre-stellar cores at three far-infrared wavelengths, 90, 170 and 200 μm. Most of the cores are detected clearly at 170 and 200 μm, but only one is detected strongly at 90 μm, indicating that mostly they are very cold, with typical temperatures of only ∼     . Colour temperature images are constructed for each of the cores. Most of the cores are seen either to be isothermal, or to have associated temperature gradients from the core centres to their edges, with all except one being cooler at the centre. We compare the data with previous ISOCAM absorption data, and calculate the energy balance for those cores in common between the two samples. We find that the energy radiated by each core in the far-infrared is similar to that absorbed at shorter wavelengths. Hence there is no evidence for a central heating source in any of the cores – even those for which previous evidence for core contraction exists. This is all consistent with external heating of the cores by the local interstellar radiation field, confirming their pre-stellar nature.     

Overtones of silicate and aluminate minerals and the 5–8 μm ice bands of deeply embedded objects

The distinct patterns, relatively low intensities and peak positions of overtone-combination bands of silicates and oxides suggest that the 5–8 μm spectral region can provide clues for the dust composition when near optically thick conditions exist for the 10-μm silicate feature. We present 1000–2500 cm⁻¹ room-temperature laboratory spectra obtained from powders of silicate, aluminate and nitride minerals and silicate glasses. The spectra exhibit overtone absorption bands with mass absorption coefficients ∼100 times weaker than the fundamentals. These data are compared with the 5–8 μm spectra of deeply embedded young stellar objects observed with the Short Wavelength Spectrometer on the Infrared Space Observatory . Fits of the laboratory data to the observations, after subtraction of the 6.0-μm H₂O ice feature and the 6.0-μm feature identified with organic refractory material, indicate that crystalline melilite (a silicate) or metamict hibonite (a radiation-damaged crystalline aluminate) may be responsible for much of the 6.9-μm absorption feature in the observations, with melilite providing the best match. A weaker 6.2-μm absorption in the young stellar object spectra is well matched by the spectra of hydrous crystalline amphibole silicates (actinolite and tremolite). Relative abundances of Si–O in room-temperature amphiboles to low-temperature H₂O ice are in the range 0.46–3.9 and in melilite are in the range 2.5–8.6. No astronomical feature was matched by the overtones of amorphous silicates because these bands are too broad and peak at the wrong wavelength. Hence, this analysis is consistent with the 10-μm features of these objects being due to a mixture of crystalline and amorphous silicates, rather than only amorphous silicates.     

Millimetre and submillimetre continuum observations of the core and hotspots of Cygnus A

We present millimetre photometry and submillimetre imaging of the central core and two hotspots in the radio lobes of the galaxy Cygnus A. For both hotspots and the central core, the synchrotron spectrum continues smoothly from the radio to a frequency of 677 GHz. The spectral index of the hotspots is constant over our frequency range, with a spectral index of α ≈ −1.0 ( S _ν ∝ ν^α), which is steeper than at lower frequencies and represents the emission from an aged population of electrons. The core is significantly flatter, with α = −0.6 ± 0.1, suggestive of an injected spectrum with no ageing, but some evidence for steepening exists at our highest observing frequency. Although IRAS data suggest the presence of dust in Cygnus A, our 450-μm data show no evidence of cold dust, therefore the dust component must have a temperature lying between 85 and 37 K, corresponding to dust masses of 1.4 × 10⁶ and 1.0 × 10⁸ M_⊙ respectively.     

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.     

The link between SCUBA and Spitzer: cold galaxies at z≲ 1

We show that the far-IR properties of distant Luminous and UltraLuminous InfraRed Galaxies (LIRGs and ULIRGs, respectively) are on average divergent from analogous sources in the local Universe. Our analysis is based on Spitzer Multiband Imaging Photometer (MIPS) and Infrared Array Camera (IRAC) data of   L _IR > 10¹⁰ L_⊙, 70 μm  selected objects in the  0.1 < z < 2  redshift range and supported by a comparison with the IRAS Bright Galaxy Sample. The majority of the objects in our sample are described by spectral energy distributions (SEDs) which peak at longer wavelengths than local sources of equivalent total infrared luminosity. This shift in SED peak wavelength implies a noticeable change in the dust and/or star-forming properties from   z ∼ 0  to the early Universe, tending towards lower dust temperatures, indicative of strong evolution in the cold dust, 'cirrus', component. We show that these objects are potentially the missing link between the well-studied local IR-luminous galaxies, Spitzer IR populations and SCUBA sources – the   z < 1  counterparts of the cold   z > 1  SubMillimetre Galaxies (SMGs) discovered in blank-field submillimetre surveys. The Herschel Space Observatory is well placed to fully characterize the nature of these objects, as its coverage extends over a major part of the far-IR/sub-mm SED for a wide redshift range.     

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.     

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.