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.     

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.     

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.     

SCUBA observations of NGC 1275

Deep SCUBA observations of NGC 1275 at 450 and 850 μm along with the application of deconvolution algorithms have permitted us to separate the strong core emission in this galaxy from the fainter extended emission around it. The core has a steep spectral index and is likely caused primarily by the active galactic nucleus. The faint emission has a positive spectral index and is clearly caused by extended dust in a patchy distribution out to a radius of ∼20 kpc from the nucleus. These observations have now revealed that a large quantity of dust, ∼     (two orders of magnitude larger than that inferred from previous optical absorption measurements), exists in this galaxy. We estimate the temperature of this dust to be ∼20 K (using an emissivity index of     and the gas/dust ratio to be 360. These values are typical of spiral galaxies. The dust emission correlates spatially with the hot X-ray emitting gas, which may be a result of collisional heating of broadly distributed dust by electrons. As the destruction time-scale is short, the dust cannot be replenished by stellar mass loss and must be externally supplied, via either the infalling galaxy or the cooling flow itself.     

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.     

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.     

Submillimetre imaging of NGC 3079

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

Submillimetre photometry of post-asymptotic giant branch stars

Stars in the post-asymptotic giant branch (post-AGB) phase of evolution are surrounded by detached circumstellar envelopes containing dust which emits thermally in the mid- and far-infrared. Here we present 850-μm SCUBA photometry of nine candidate post-AGB stars. All targets are detected at 850 μm and we use these fluxes to estimate the envelope dust masses and, by comparison with the 100-μm IRAS fluxes, the dust emissivity index.     

The distribution of atomic gas and dust in nearby galaxies – I. Presentation of matched-resolution VLA H i and SCUBA 850-μm maps

We present matched-resolution VLA H  i and SCUBA 850-μm maps of 20 IRAS -bright galaxies. Of the galaxies observed, two were not detected in H  i and two were detected in absorption. The H  i distributions of the galaxies have a range of morphologies. Some of the systems appear H  i deficient in the central regions which could be due to a high conversion rate of H  i into molecules or H  i absorption. In contrast to the H  i , the 850-μm emission has a smooth distribution which is concentrated towards the optical centre of each galaxy. We also find evidence for 850-μm emission extending to the periphery of the optical disc in some of the galaxies. Finally, we note that the relative lack of 850-μm emission when compared with H  i does not necessarily mean that the atomic gas and dust do not have similar mass distributions.     

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.     

On the lack of cold dust in IRAS P09104+4109 and IRAS F15307+3252: their spectral energy distributions and implications for finding dusty AGNs at high redshift

We present upper limits on the 850-μm and 450-μm fluxes of the warm hyperluminous (bolometric luminosity     galaxies IRAS P09104+4109     and IRAS F15307+3252     , derived from measurements using the SCUBA bolometer array on the James Clerk Maxwell Telescope. Hot luminous infrared sources like these are thought to differ from more normal cold ultraluminous infrared     galaxies in that they derive most of their bolometric luminosities from dusty active galactic nuclei (AGNs) as opposed to starbursts. Such hot, dusty AGNs at high redshift are thought to be responsible for much of the mass accretion of the Universe that is in turn responsible for the formation of the supermassive black holes seen in the centres of local galaxies. The galaxy IRAS P09104+4109 is also unusual in that it is a cD galaxy in the centre of a substantial cooling-flow cluster, not an isolated interacting galaxy like most ultraluminous infrared galaxies. Previously it was known to have large amounts of hot     dust from IRAS observations. We now show that the contribution of cold dust to the bolometric luminosity is less than 3 per cent. Most ultraluminous infrared galaxies possess large amounts of cold dust, and it is now known that some cooling-flow cluster cD galaxies do as well. Yet this object, which is an extreme example of both, does not have enough cold gas to contribute significantly to the bolometric luminosity. We outline physical reasons why this could have happened. We then provide a discussion of strategies for finding hot dusty AGNs, given the limitations on submillimetre surveys implied by this work.     

Study of star formation in RCW 106 using far-infrared observations

High-resolution far-infrared observations of a large area of the star-forming complex RCW 106 obtained using the TIFR 1-m balloon-borne telescope are presented. Intensity maps have been obtained simultaneously in two bands centred around 150 and 210 μm. Intensity maps have also been obtained in the four IRAS bands using HIRES-processed IRAS data. From the 150- and 210-μm maps, reliable maps of dust temperature and optical depth have been generated. The star formation in this complex has occurred in five linear sub-clumps. Using the map at 210 μm, which has a spatial resolution superior to that of IRAS at 100 μm, 23 sources have been identified. The spectral energy distribution (SED) and luminosity of these sources have been determined using the associations with the IRAS maps. The luminosity distribution of these sources has been obtained. Assuming these embedded sources to be zero-age main-sequence stars and using the mass–luminosity relation for these, the power-law slope of the initial mass function is found to be −1.73±0.5 . This index for this very young complex is about the same as that for more evolved complexes and clusters. Radiation transfer calculations in spherically symmetric geometry have been undertaken to fit the SEDs of 13 sources with fluxes in both the TIFR and the IRAS bands. From this, the r ⁻² density distribution in the envelopes is ruled out. Finally, a correlation is seen between the luminosity of embedded sources and the computed dust masses of the envelopes.     

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

A multiwavelength infrared study of NGC 891

We present a multiwavlength infrared (IR) study of the nearby, edge-on, spiral galaxy NGC 891. We have examined 20 independent, spatially resolved IR images of this galaxy, 14 of which are newly reduced and/or previously unpublished images. These images span a wavelength regime from  λ 1.2 μ  m in which the emission is dominated by cool stars, through the mid-IR, in which emission is dominated by polycyclic aromatic hydrocarbons (PAHs), to λ 850 μm, in which emission is dominated by cold dust in thermal equilibrium with the radiation field. The changing morphology of the galaxy with wavelength illustrates the changing dominant components. We detect extraplanar dust emission in this galaxy, consistent with previously published results, but now show that PAH emission is also in the halo, to a vertical distance of   z ≥ 2.5 kpc  . We compare the vertical extents of various components and find that the PAHs (from λ 7.7 and 8 μm data) and warm dust (λ 24 μm) extend to smaller z heights than the cool dust (λ 450 μm). For six locations in the galaxy for which the signal-to-noise ratio was sufficient, we present spectral energy distributions (SEDs) of the IR emission, including two in the halo – the first time a halo SED in an external galaxy has been presented. We have modelled these SEDs and find that the PAH fraction, f _PAH, is similar to Galactic values (within a factor of 2), with the lowest value at the galaxy's centre, consistent with independent results of other galaxies. In the halo environment, the fraction of dust exposed to a colder radiation field, f _cold, is of the order of unity, consistent with an environment in which there is no star formation. The source of excitation is likely from photons escaping from the disc.     

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.     

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.     

Three quasars from a survey of strong 25-μm emitters

We have carried out a spectroscopic survey of 750 sources that are strong 25-μm emitters from the IRAS Faint Source data base. Many of these sources are previously unknown active galactic nuclei including new IRAS quasars, three of which we describe here: F21382−2659, Z06367−6845 and Z05558−5008. They are all radio and X-ray quiet, and compared to the known IRAS quasars they have similar 25-μm luminosities, L (25 μm), but lower values of L (25 μm)/ L ( B ). Their F (25 μm)/ F (60 μm) IRAS colours lie in the range 0.33 to 1.08, indicating the presence of relatively warm dust, presumably in a dusty torus surrounding the central source, and with temperatures similar to those of the known IRAS quasars. The quasar with the warmest dust, F21382−2659, exhibits broad (full width at half-maximum ∼4000 km s⁻¹) asymmetric Balmer lines with H γ having an opposite asymmetry to the other broad lines; also H β (only) is double-peaked. Fe  ii is very weak in F21382−2659 but strong in the other two quasars, and the anticorrelation between Fe  ii and [O  iii ] holds as anticipated. Two of the quasars are unpolarized: although F21382−2659 is optically polarized (2.1 per cent at 4950 Å), we argue that this provides little insight into the orientation of its dust torus relative to the line of sight.     

The contribution of faint blue galaxies to the submillimetre counts and background

Observations in the submillimetre (submm) waveband have recently revealed a new population of luminous sources. These are proposed to lie at high redshift and to be optically faint because of their high intrinsic dust obscuration. The presence of dust has been previously invoked in optical galaxy count models which use the Bruzual & Charlot evolution models with an exponential τ =9 Gyr star formation rate (SFR) for spirals, and these fit the count data well from U to K . We now show that by using either a 1/ λ or Calzetti absorption law for the dust and re-distributing the evolved spiral galaxy ultraviolet (UV) radiation into the far-infrared (FIR), these models can account for all of the 'faint' ( 1 mJy) 850-μm galaxy counts, but fail to fit 'bright' ( 2 mJy) sources, indicating that another explanation for the submm counts may apply at brighter fluxes, e.g., quasi-stellar objects (QSOs) or ultraluminous infrared galaxies (ULIRGs). We find that the main contribution to the faint, submm number counts is in the redshift range 0.5< z <3, peaking at z ≈1.8. The above model, using either dust law, can also explain a significant proportion of the extragalactic background at 850 μm, as well as producing a reasonable fit to the bright 60-μm IRAS counts.     

Starburst galaxies and structure in the submillimetre background towards the Hubble Deep Field

We use an 850-μm SCUBA map of the Hubble Deep Field (HDF) to study the dust properties of optically-selected starburst galaxies at high redshift. The optical/infrared (IR) data in the HDF allow a photometric redshift to be estimated for each galaxy, together with an estimate of the visible star-formation rate. The 850-μm flux density of each source provides the complementary information: the amount of hidden, dust-enshrouded star formation activity. Although the 850-μm map does not allow detection of the majority of individual sources, we show that the galaxies with the highest UV star-formation rates are detected statistically, with a flux density of about S ₈₅₀=0.2 mJy for an apparent UV star-formation rate of 1  h ⁻² M_⊙ yr⁻¹. This level of submillimetre output indicates that the total star-forming activity is on average a factor of approximately 6 times larger than the rate inferred from the UV output of these galaxies. The general population of optical starbursts is then predicted to contribute at least 25 per cent of the 850-μm background. We carry out a power-spectrum analysis of the map, which yields some evidence for angular clustering of the background source population, but at a level lower than that seen in Lyman-break galaxies. Together with other lines of argument, particularly from the NICMOS HDF data, this suggests that the 850-μm background originates over an extremely wide range of redshifts – perhaps 1≲ z ≲6.     

The infrared luminosity of the torus and the visibility of scattered broad line emission in Seyfert 2 galaxies

A number of studies have shown that the visibility of scattered broad emission lines in Seyfert 2 galaxies is strongly dependent on the IRAS     flux ratio, where those Seyfert 2 galaxies with 'warm' IRAS colours show polarized broad line emission. It is now clear that this effect is owing to the increasing dominance of the galactic rather than the active galactic nucleus (AGN) emission at 60 μm in less-luminous 'cool' Seyfert 2 galaxies. However, we present evidence that the 25-μm emission is a good measure of the AGN luminosity for most Seyfert 2 galaxies. Using this result, we show that the visibility of scattered broad line emission has a dependence on the AGN luminosity. The observations can be interpreted self-consistently if the scaleheight of the scattering zone varies with central source luminosity whilst the scaleheight of the obscuring torus is approximately constant.