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 contribution of faint blue galaxies to the submillimetre counts and background

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

The SCUBA HAlf Degree Extragalactic Survey – VI. 350-μm mapping of submillimetre galaxies

A follow-up survey using the Submillimetre High-Angular Resolution Camera (SHARC-II) at 350 μm has been carried out to map the regions around several 850-μm-selected sources from the Submillimetre HAlf Degree Extragalactic Survey (SHADES). These observations probe the infrared (IR) luminosities and hence star formation rates in the largest existing, most robust sample of submillimetre galaxies (SMGs). We measure 350-μm flux densities for 24 850-μm sources, seven of which are detected at ≥2.5σ within a 10 arcsec search radius of the 850-μm positions. When results from the literature are included the total number of 350-μm flux density constraints of SHADES SMGs is 31, with 15 detections. We fit a modified blackbody to the far-IR (FIR) photometry of each SMG, and confirm that typical SMGs are dust-rich  ( M _dust≃ 9 × 10⁸ M_⊙)  , luminous  ( L _FIR≃ 2 × 10¹² L_⊙)  star-forming galaxies with intrinsic dust temperatures of ≃35 K and star formation rates of  ≃400 M_⊙ yr⁻¹  . We have measured the temperature distribution of SMGs and find that the underlying distribution is slightly broader than implied by the error bars, and that most SMGs are at 28 K with a few hotter. We also place new constraints on the 350-μm source counts, N ₃₅₀(>25 mJy) ∼ 200–500 deg⁻².     

The SCUBA Local Universe Galaxy Survey – II. 450-μm data: evidence for cold dust in bright IRAS galaxies

This is the second in a series of papers presenting results from the SCUBA Local Universe Galaxy Survey. In our first paper we provided 850-μm flux densities for 104 galaxies selected from the IRAS Bright Galaxy Sample and we found that the 60-, 100-μm ( IRAS ) and 850-μm (SCUBA) fluxes could be adequately fitted by emission from dust at a single temperature. In this paper we present 450-μm data for the galaxies. With the new data, the spectral energy distributions of the galaxies can no longer be fitted with an isothermal dust model – two temperature components are now required. Using our 450-μm data and fluxes from the literature, we find that the 450/850-μm flux ratio for the galaxies is remarkably constant, and this holds from objects in which the star formation rate is similar to our own Galaxy, to ultraluminous infrared galaxies (ULIRGs) such as Arp 220. The only possible explanation for this is if the dust emissivity index for all of the galaxies is ∼2 and the cold dust component has a similar temperature in all galaxies     . The 60-μm luminosities of the galaxies were found to depend on both the dust mass and the relative amount of energy in the warm component, with a tendency for the temperature effects to dominate at the highest L ₆₀. The dust masses estimated using the new temperatures are higher by a factor of ∼2 than those determined previously using a single temperature. This brings the gas-to-dust ratios of the IRAS galaxies into agreement with those of the Milky Way and other spiral galaxies which have been intensively studied in the submm.     

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.     

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.     

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.     

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.     

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.     

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.     

MAMBO observations at 240 GHz of optically obscured Spitzer sources: source clumps and radio activity at high redshift

Optically very faint  ( R > 25.5)  sources detected by the Spitzer Space Telescope at 24 μm represent a very interesting population at redshift   z ∼ (1.5–3)  . They exhibit strong clustering properties, implying that they are hosted by very massive haloes, and their mid-infrared emission could be powered by either dust-enshrouded star formation and/or by an obscured active galactic nucleus (AGN). We report observations carried out with the Max Planck Millimetre Bolometer (MAMBO) array at the IRAM 30-m antenna on Pico Veleta of a candidate protocluster with five optically obscured sources selected from the 24-μm Spitzer sample of the First-Look Survey. Interestingly, these sources appear to lie on a high-density filament aligned with the two radio jets of an AGN. Four out of five of the observed sources were detected. We combine these measurements with optical, infrared and radio observations to probe the nature of the candidate protocluster members. Our preliminary conclusions can be summarized as follows: the spectral energy distributions (SEDs) of all sources include both AGN and starburst contributions; the AGN contribution to the bolometric luminosities ranges between 14 and 26 per cent of the total. Such a contribution is enough for the AGN to dominate the emission at 5.8, 8 and 24 μm, while the stellar component, inferred from SED fitting, prevails at 1.25 mm and at  λ < 4.5 μ  m. The present analysis suggests a coherent interplay at high z between extended radio activity and the development of filamentary large-scale structures.     

A far-infrared view of the Lockman hole from ISO 95-μm observations – II. Optical identifications and insights into the nature of the far-infrared sources

We present the optical identifications of a 95-μm ISOPHOT sample in the Lockman hole over an area of approximately half a deg². The Rodighiero et al. catalogue includes 36 sources, making up a complete flux-limited sample for   S _95 μm≥ 100 mJy  . Reliable sources were detected, with decreasing but well-controlled completeness, down to   S _95 μm≃ 20 mJy  . We have combined mid-infrared (IR) and radio catalogues in this area to identify the potential optical counterparts of the far-IR sources. We found 14 radio and 13 15-μm associations, 10 of which have both associations. For the 11 sources with spectroscopic redshift, we have performed a spectrophotometric analysis of the observed spectral energy distributions (SEDs). Four of these 95-μm sources have been classified as faint IR (FIR) galaxies  ( L _FIR < 1. e 11 L_⊙)  , six as luminous IR galaxies (LIRGs) and only one as an ultraluminous IR galaxy (ULIRG). We have discussed the redshift distribution of these objects, comparing our results with evolutionary model predictions 95 and 175 μm. Given their moderate distances (the bulk of the closest spectroscopically identified objects lying at   z < 0.2  ), their luminosities and star formation rates (SFR; median value  ∼ 10 M_⊙ yr⁻¹  ), the sources unveiled by ISOPHOT at 95 μm seem to correspond to the low redshift  ( z < 0.3)  FIRBACK 175-μm population, composed of dusty, star-forming galaxies with moderate SFRs. We computed and compared different SFR estimators, and found that the SF derived from the bolometric IR luminosity is well correlated with that computed from the radio and mid-IR fluxes.     

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.     

Variations in 24-μm morphologies among galaxies in the Spitzer Infrared Nearby Galaxies Survey: new insights into the Hubble sequence

To study the distribution of star formation and dust emission within nearby galaxies, we measured five morphological parameters in the 3.6- and 24-μm wavebands for 65 galaxies in the Spitzer Infrared Nearby Galaxies Survey (SINGS) and eight galaxies that were serendipitously observed by SINGS. The morphological parameters demonstrate strong variations along the Hubble sequence, including statistically significant differences between S0/a-Sab and Sc-Sd galaxies. Early-type galaxies are generally found to be compact, centralized, symmetric sources in the 24-μm band, while late-type galaxies are generally found to be extended, asymmetric sources. These results suggest that the processes that increase the real or apparent sizes of galaxies' bulges also lead to more centralized 24-μm dust emission. Several phenomena, such as strong nuclear star formation, Seyfert activity, or outer ring structures, may cause galaxies to deviate from the general morphological trends observed at 24 μm. We also note that the 24-μm morphologies of Sdm-Im galaxies are quite varied, with some objects appearing very compact and symmetric but others appearing diffuse and asymmetric. These variations reflect the wide variation in star formation in irregular galaxies as observed at other wavelengths. The variations in the 24-μm morphological parameters across the Hubble sequence mirror many of the morphological trends seen in other tracers of the ISM and in stellar emission. However, the 24-μm morphological parameters for the galaxies in this sample do not match the morphological parameters measured in the stellar wavebands. This implies that the distribution of dust emission is related to but not equivalent to the distribution of stellar emission.     

The temperatures of dust-enshrouded active galactic nuclei

A high density of massive dark objects (MDOs), probably supermassive black holes, in the centres of nearby galaxies has been inferred from recent observations. There are various indications that much of the accretion responsible for producing these objects took place in dust-enshrouded active galactic nuclei (AGNs). If so, then measurements of the intensity of background radiation and the source counts in the far-infrared and submillimetre wavebands constrain the temperature of dust in these AGNs. An additional constraint comes from the hard X-ray background, if this is produced by accretion. One possibility is that the dust shrouds surrounding the accreting AGNs are cold, about 30 K. In this event, the dusty AGNs could be some subset of the population of luminous distant sources discovered at 850 μm using the SCUBA array on the James Clerk Maxwell Telescope, as proposed by Almaini, Lawrence & Boyle. An alternative is that the dust shrouds surrounding the accreting AGNs are much hotter (>60 K). These values are closer to the dust temperatures of a number of well-studied low-redshift ultraluminous galaxies that are thought to derive their power from accretion. If the local MDO density is close to the maximum permitted, then cold sources cannot produce this density without the submillimetre background being overproduced if they accrete at high radiative efficiency, and thus a hot population is required. If the dust-enshrouded accretion occurred at similar redshifts to that taking place in unobscured optical quasars, then a significant fraction of the far-infrared background radiation measured by COBE at 140 μm, but very little of the submillimetre background at 850 μm, may have been produced by hot dust-enshrouded AGNs which may have already been seen in recent Chandra X-ray surveys.     

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.     

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

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

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

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.     

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.