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.     

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  .     

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.     

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.     

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.     

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

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.     

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.     

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.     

High-resolution radio observations of Seyfert galaxies in the extended 12-μm sample – II. The properties of compact radio components

We discuss the properties of compact nuclear radio components in Seyfert galaxies from the extended 12-μm AGN sample of Rush et al. Our main results can be summarized as follows.
Type 1 and type 2 Seyferts produce compact radio components which are indistinguishable in strength and aspect, indicating that their central engines are alike, as proposed by the unification model. Infrared IRAS fluxes are more closely correlated with low-resolution radio fluxes than high-resolution radio fluxes, suggesting that they are dominated by kiloparsec-scale, extranuclear emission regions; extranuclear emission may be stronger in type 2 Seyferts. Early-type Seyfert galaxies tend to have stronger nuclear radio emission than late-type Seyfert galaxies. V-shaped extended emission-line regions, indicative of 'ionization cones', are usually found in sources with large, collimated radio outflows. Hidden broad lines are most likely to be found in sources with powerful nuclear radio sources. Type 1 and type 2 Seyferts selected by their IRAS 12-μm flux densities have well-matched properties.     

The SCUBA Half-Degree Extragalactic Survey (SHADES) – VIII. The nature of faint submillimetre galaxies in SHADES, SWIRE and SXDF surveys

We present the optical-to-submillimetre spectral energy distributions (SEDs) for 33 radio and mid-infrared (mid-IR) identified submillimetre galaxies discovered via the SHADES 850-μm SCUBA imaging in the Subaru- XMM Deep Field (SXDF). Optical data for the sources come from the SXDF and mid- and far-IR fluxes from SWIRE. We obtain photometric redshift estimates for our sources using optical and IRAC 3.6- and 4.5-μm fluxes. We then fit SED templates to the longer wavelength data to determine the nature of the far-IR emission that dominates the bolometric luminosity of these sources. The IR template fits are also used to resolve ambiguous identifications and cases of redshift aliasing. The redshift distribution obtained broadly matches previous results for submillimetre sources and on the SHADES SXDF field. Our template fitting finds that active galactic nuclei, while present in about 10 per cent of our sources, do not contribute significantly to their bolometric luminosity. Dust heating by starbursts, with either Arp220 or M82 type SEDs, appears to be responsible for the luminosity in most sources (23/33 are fitted by Arp220 templates, 2/33 by the warmer M82 templates). 8/33 sources, in contrast, are fitted by a cooler cirrus dust template, suggesting that cold dust has a role in some of these highly luminous objects. Three of our sources appear to have multiple identifications or components at the same redshift, but we find no statistical evidence that close associations are common among our SHADES sources. Examination of rest-frame K -band luminosity suggests that 'downsizing' is underway in the submillimetre galaxy population, with lower redshift systems lying in lower mass host galaxies. Of our 33 identifications six are found to be of lower reliability but their exclusion would not significantly alter our conclusions.     

The 3.4-μm absorption feature towards the nucleus of NGC 5506

We report the results of 3-μm spectroscopy towards the nucleus of a nearby Seyfert 2 galaxy, NGC 5506. A weak absorption feature of carbonaceous dust is detected at ∼3.4 μm. The optical depth ratio of the 3.4-μm carbonaceous dust absorption to the 9.7-μm silicate dust absorption is smaller by more than a factor of 2 than that in the interstellar medium in our Galaxy. The small ratio is consistent with the idea that the obscuration towards the nucleus of NGC 5506 is caused mainly by dust grains in its host galaxy (type S0/a) and that the contribution of carbonaceous dust grains to the visual extinction there is smaller than that in our Galaxy (type Sb/bc).     

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 3.4-μm absorption feature towards three obscured active galactic nuclei

The results of 3–4-μm spectroscopy towards the nuclei of NGC 3094, 7172, and 7479 are reported. In ground-based 8–13-μm spectra, all the sources have strong absorption-like features at ∼10 μm, but they do not have detectable polycyclic aromatic hydrocarbon (PAH) emission features. The 3.4-μm carbonaceous dust absorption features are detected towards all nuclei. NGC 3094 shows a detectable 3.3-μm PAH emission feature, while NGC 7172 and 7479 do not. Nuclear emission whose spectrum shows dust absorption features but no PAH emission features, is thought to be dominated by highly obscured active galactic nuclei (AGNs) activity. For NGC 7172, 7479, and three other such nuclei in the literature, we investigate the optical depth ratios between the 3.4-μm carbonaceous dust and 9.7-μm silicate dust absorption     The     ratios towards three highly obscured AGNs with face-on host galaxies are systematically larger than the ratios in the Galactic diffuse interstellar medium or the ratios for two highly obscured AGNs with edge-on host galaxies. We suggest that the larger ratios can be explained if the obscuring dust is so close to the central AGNs that a temperature gradient occurs in it. If this idea is correct, our results may provide spectroscopic evidence for the presence of the putative 'dusty tori' in the close vicinity of AGNs.     

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.     

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  .     

NGC 1333/IRAS 4: a multiple star formation laboratory

We present SCUBA observations of the protomultiple system NGC 1333/IRAS 4 at 450 and 850 μm. The 850-μm map shows significant extended emission which is most probably a remnant of the initial cloud core. At 450 μm, the component 4A is seen to have an elongated shape suggestive of a disc. Also we confirm that, in addition to the 4A and 4B system, there exists another component 4C, which appears to lie out of the plane of the system and of the extended emission. Deconvolution of the beam reveals a binary companion to IRAS 4B. Simple considerations of binary dynamics suggest that this triple 4A–4BI–4BII system is unstable and will probably not survive in its current form. Thus IRAS 4 provides evidence that systems can evolve from higher to lower multiplicity as they move towards the main sequence. We construct a map of spectral index from the two wavelengths, and comment on the implications of this for dust evolution and temperature differences across the map. There is evidence that in the region of component 4A the dust has evolved, probably by coagulating into larger or more complex grains. Furthermore, there is evidence from the spectral index maps that dust from this object is being entrained in its associated outflow.     

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.     

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.