D. L. Clements

The SCUBA Half-Degree Extragalactic Survey - II. Submillimetre maps, catalogue and number counts

We present maps, source catalogue and number counts of the largest, most complete and unbiased extragalactic submillimetre survey: the 850-μm SCUBA Half-Degree Extragalactic Survey (SHADES). Using the Submillimetre Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT), SHADES mapped two separate regions of sky: the Subaru/XMM–Newton Deep Field (SXDF) and the Lockman Hole East (LH). Encompassing 93 per cent of the overall acquired data (i.e. data taken up to 2004 February 1), these SCUBA maps cover 720 arcmin2 with a rms noise level of about 2 mJy and have uncovered >100 submillimetre galaxies. In order to ensure the utmost robustness of the resulting source catalogue, data reduction was independently carried out by four subgroups within the SHADES team, providing an unprecedented degree of reliability with respect to other SCUBA catalogues available from the literature. Individual source lists from the four groups were combined to produce a robust 120-object SHADES catalogue; an invaluable resource for follow-up campaigns aiming to study the properties of a complete and consistent sample of submillimetre galaxies. For the first time, we present deboosted flux densities for each submillimetre galaxy found in a large survey. Extensive simulations and tests were performed separately by each group in order to confirm the robustness of the source candidates and to evaluate the effects of false detections, completeness and flux density boosting. Corrections for these effects were then applied to the data to derive the submillimetre galaxy source counts. SHADES has a high enough number of detected sources that meaningful differential counts can be estimated, unlike most submillimetre surveys which have to consider integral counts. We present differential and integral source number counts and find that the differential counts are better fit with a broken power law or a Schechter function than with a single power law; the SHADES data alone significantly show that a break is required at several mJy, although the precise position of the break is not well constrained. We also find that a 850-μm survey complete down to 2 mJy would resolve 20–30 per cent of the far-infrared background into point sources.

A dust-obscured massive maximum-starburst galaxy at a redshift of 6.34

Massive present-day early-type (elliptical and lenticular) galaxies probably gained the bulk of their stellar mass and heavy elements through intense, dust-enshrouded starbursts—that is, increased rates of star formation—in the most massive dark-matter haloes at early epochs. However, it remains unknown how soon after the Big Bang massive starburst progenitors exist. The measured redshift (z) distribution of dusty, massive starbursts has long been suspected to be biased low in z owing to selection effects, as confirmed by recent findings of systems with redshifts as high as ∼5 (refs 2–4). Here we report the identification of a massive starburst galaxy at z = 6.34 through a submillimetre colour-selection technique. We unambiguously determined the redshift from a suite of molecular and atomic fine-structure cooling lines. These measurements reveal a hundred billion solar masses of highly excited, chemically evolved interstellar medium in this galaxy, which constitutes at least 40 per cent of the baryonic mass. A ‘maximum starburst’ converts the gas into stars at a rate more than 2,000 times that of the Milky Way, a rate among the highest observed at any epoch. Despite the overall downturn in cosmic star formation towards the highest redshifts, it seems that environments mature enough to form the most massive, intense starbursts existed at least as early as 880 million years after the Big Bang.

The SCUBA HAlf Degree Extragalactic Survey – III. Identification of radio and mid-infrared counterparts to submillimetre galaxies

Determining an accurate position for a submillimetre (submm) galaxy (SMG) is the crucial step that enables us to move from the basic properties of an SMG sample – source counts and 2D clustering – to an assessment of their detailed, multiwavelength properties, their contribution to the history of cosmic star formation and their links with present-day galaxy populations. In this paper, we identify robust radio and/or infrared (IR) counterparts, and hence accurate positions, for over two-thirds of the SCUBA HAlf-Degree Extragalactic Survey (SHADES) Source Catalogue, presenting optical, 24-μm and radio images of each SMG. Observed trends in identification rate have given no strong rationale for pruning the sample. Uncertainties in submm position are found to be consistent with theoretical expectations, with no evidence for significant additional sources of error. Employing the submm/radio redshift indicator, via a parametrization appropriate for radio-identified SMGs with spectroscopic redshifts, yields a median redshift of 2.8 for the radio-identified subset of SHADES, somewhat higher than the median spectroscopic redshift. We present a diagnostic colour–colour plot, exploiting Spitzer photometry, in which we identify regions commensurate with SMGs at very high redshift. Finally, we find that significantly more SMGs have multiple robust counterparts than would be expected by chance, indicative of physical associations. These multiple systems are most common amongst the brightest SMGs and are typically separated by 2–6 arcsec, ~15–50/ sin i kpc at z∼ 2, consistent with early bursts seen in merger simulations.

HerMES : SPIRE galaxy number counts at 250, 350, and 500 μm

Emission at far-infrared wavelengths makes up a significant fraction of the total light detected from galaxies over the age of Universe. Herschel provides an opportunity for studying galaxies at the peak wavelength of their emission. Our aim is to provide a benchmark for models of galaxy population evolution and to test pre-existing models of galaxies. With the Herschel Multi-tiered Extra-galactic survey, HerMES, we have observed a number of fields of different areas and sensitivity using the SPIRE instrument on Herschel. We have determined the number counts of galaxies down to ~20 mJy. Our constraints from directly counting galaxies are consistent with, though more precise than, estimates from the BLAST fluctuation analysis. We have found a steep rise in the Euclidean normalised counts <100 mJy. We have directly resolved ~15% of the infrared extra-galactic background at the wavelength near where it peaks.

Starburst and AGN activity in ultraluminous infrared galaxies

We examine the power source of 41 local Ultraluminous Infrared Galaxies (ULIRGs) using archival infrared (IR) and optical photometry. We fit the observed Spectral Energy Distributions (SEDs) with starburst and AGN components; each component being drawn from a family of templates. We find all of the sample require a starburst, whereas only half require an AGN. In 90% of the sample the starburst provides over half the IR emission, with a mean fractional luminosity of 82%. When combined with other galaxy samples we find that starburst and AGN luminosities correlate over 6 decades in IR luminosity suggesting that a common factor governs both luminosities, plausibly the gas masses in the nuclear regions. We find no trend for increasing fractional AGN luminosity with increasing total luminosity, contrary to previous claims. We find that the mid-IR F7.7/C7.7 line-continuum ratio is no indication of the starburst luminosity, or the fractional AGN luminosity, and therefore that F7.7/C7.7 is not a reliable diagnostic of the power source in ULIRGs. The radio flux correlates with the starburst luminosity, but shows no correlation with the AGN luminosity, in line with previous results. We propose that the scatter in this correlation is due to a skewed starburst IMF and/or relic relativistic electrons from a previous starburst, rather than contamination from an obscured AGN. We show that most ULIRGs undergo multiple starbursts during their lifetime, and by inference that mergers between more than two galaxies may be common amongst ULIRGs. Our results support the evolutionary model for ULIRGs proposed by Farrah et al (2001), where they can follow many different evolutionary paths of starburst and AGN activity in transforming merging spiral galaxies into elliptical galaxies, but that most do not go through an optical QSO phase. The lower level of AGN activity in our local sample than in z � 1 HLIRGs implies that the two samples are distinct populations. We postulate that different galaxy formation processes at high-z are responsible for this difference.

Herschel ⋆ -ATLAS: Rapid evolution of dust in galaxies over the last 5 billion years

We present the first direct and unbiased measurement of the evolution of the dust mass function of galaxies over the past 5 billion years of cosmic history using data from the Science Demonstration Phase of the Herschel-Astrophysical Terahertz Large Area Survey (Herschel-ATLAS). The sample consists of galaxies selected at 250 m which have reliable counterparts from the Sloan Digital Sky Survey (SDSS) at z < 0.5, and contains 1867 sources. Dust masses are calculated using both a single-temperature grey-body model for the spectral energy distribution and also a model with multiple temperature components. The dust temperature for either model shows no trend with redshift. Splitting the sample into bins of redshift reveals a strong evolution in the dust properties of the most massive galaxies. At z= 0.4–0.5, massive galaxies had dust masses about five times larger than in the local Universe. At the same time, the dust-to-stellar mass ratio was about three to four times larger, and the optical depth derived from fitting the UV-sub-mm data with an energy balance model was also higher. This increase in the dust content of massive galaxies at high redshift is difficult to explain using standard dust evolution models and requires a rapid gas consumption time-scale together with either a more top-heavy initial mass function (IMF), efficient mantle growth, less dust destruction or combinations of all three. This evolution in dust mass is likely to be associated with a change in overall interstellar medium mass, and points to an enhanced supply of fuel for star formation at earlier cosmic epochs.

The suppression of star formation by powerful active galactic nuclei.

The old, red stars that constitute the bulges of galaxies, and the massive black holes at their centres, are the relics of a period in cosmic history when galaxies formed stars at remarkable rates and active galactic nuclei (AGN) shone brightly as a result of accretion onto black holes. It is widely suspected, but unproved, that the tight correlation between the mass of the black hole and the mass of the stellar bulge results from the AGN quenching the surrounding star formation as it approaches its peak luminosity. X-rays trace emission from AGN unambiguously, whereas powerful star-forming galaxies are usually dust-obscured and are brightest at infrared and submillimetre wavelengths. Here we report submillimetre and X-ray observations that show that rapid star formation was common in the host galaxies of AGN when the Universe was 2–6 billion years old, but that the most vigorous star formation is not observed around black holes above an X-ray luminosity of 1044 ergs per second. This suppression of star formation in the host galaxy of a powerful AGN is a key prediction of models in which the AGN drives an outflow, expelling the interstellar medium of its host and transforming the galaxy’s properties in a brief period of cosmic time.

The dust scaling relations of the Herschel Reference Survey

We combine new Herschel/SPIRE sub-millimeter observations with existing multiwavelength data to investigate the dust scaling relations of the Herschel Reference Survey, a magnitude-, volume-limited sample of similar to 300 nearby galaxies in different environments. We show that the dust-to-stellar mass ratio anti-correlates with stellar mass, stellar mass surface density and NUV - r colour across the whole range of parameters covered by our sample. Moreover, the dust-to-stellar mass ratio decreases significantly when moving from late-to early-type galaxies. These scaling relations are similar to those observed for the Hi gas-fraction, supporting the idea that the cold dust is tightly coupled to the cold atomic gas component in the interstellar medium. We also find a weak increase of the dust-to-Hi mass ratio with stellar mass and colour but no trend is seen with stellar mass surface density. By comparing galaxies in different environments we show that, although these scaling relations are followed by both cluster and field galaxies, Hi-deficient systems have, at fixed stellar mass, stellar mass surface density and morphological type systematically lower dust-to-stellar mass and higher dust-to-Hi mass ratios than Hi-normal/field galaxies. This provides clear evidence that dust is removed from the star-forming disk of cluster galaxies but the effect of the environment is less strong than what is observed in the case of the Hi disk. Such effects naturally arise if the dust disk is less extended than the Hi and follows more closely the distribution of the molecular gas phase, i.e., if the dust-to-atomic gas ratio monotonically decreases with distance from the galactic center.

HerMES: deep number counts at 250 μm, 350 μm and 500 μm in the COSMOS and GOODS-N fields and the build-up of the cosmic infrared background

Aims. The Spectral and Photometric Imaging REceiver (SPIRE) onboard the Herschel space telescope has provided confusion limited maps of deep fields at 250 μm, 350 μm, and 500 μm, as part of the Herschel Multi-tiered Extragalactic Survey (HerMES). Unfortunately, due to confusion, only a small fraction of the cosmic infrared background (CIB) can be resolved into individually-detected sources. Our goal is to produce deep galaxy number counts and redshift distributions below the confusion limit at SPIRE wavelengths (~20 mJy), which we then use to place strong constraints on the origins of the cosmic infrared background and on models of galaxy evolution. Methods. We individually extracted the bright SPIRE sources (>20 mJy) in the COSMOS field with a method using the positions, the flux densities, and the redshifts of the 24 μm sources as a prior, and derived the number counts and redshift distributions of the bright SPIRE sources. For fainter SPIRE sources (<20 mJy), we reconstructed the number counts and the redshift distribution below the confusion limit using the deep 24 μm catalogs associated with photometric redshift and information provided by the stacking of these sources into the deep SPIRE maps of the GOODS-N and COSMOS fields. Finally, by integrating all these counts, we studied the contribution of the galaxies to the CIB as a function of their flux density and redshift. Results. Through stacking, we managed to reconstruct the source counts per redshift slice down to ~2 mJy in the three SPIRE bands, which lies about a factor 10 below the 5σ confusion limit. Our measurements place tight constraints on source population models. None of the pre-existing models are able to reproduce our results at better than 3-σ. Finally, we extrapolate our counts to zero flux density in order to derive an estimate of the total contribution of galaxies to the CIB, finding 10.1_(-2.3)^(+2.6) nW m^(-2) sr^(-1), 6.5_(-1.6)^(+1.7) nW m^(-2) sr^(-1), and 2.8_(-0.8)^(+0.9) nW m^(-2) sr^(-1) at 250 μm, 350 μm, and 500 μm, respectively. These values agree well with FIRAS absolute measurements, suggesting our number counts and their extrapolation are sufficient to explain the CIB. We find that half of the CIB is emitted at z = 1.04, 1.20, and 1.25, respectively. Finally, combining our results with other works, we estimate the energy budget contained in the CIB between 8 μm and 1000 μm: 26_(-3)^(+7) nW m^(-2) sr^(-1).

The Herschel Reference Survey: dust in early-type galaxies and across the Hubble Sequence

We present Herschel observations of 62 early-type galaxies (ETGs), including 39 galaxies morphologically classified as S0+S0a and 23 galaxies classified as ellipticals using SPIRE at 250, 350, and 500 mu m as part of the volume-limited Herschel Reference Survey (HRS). We detect dust emission in 24% of the ellipticals and 62% of the S0s. The mean temperature of the dust is \textless T-d \textgreater = 23.9 +/- 0.8 K, warmer than that found for late-type galaxies in the Virgo Cluster. The mean dust mass for the entire detected early-type sample is log M-d = 6.1 +/- 0.1 M-circle dot with a mean dust-to-stellar-mass ratio of log(M-d/M-*) = -4.3 +/- 0.1. Including the non-detections, these parameters are log M-d = 5.6 +/- 0.1 and log(M-d/M-*) = -5.1 +/- 0.1, respectively. The average dust-to-stellar-mass ratio for the early-type sample is fifty times lower, with larger dispersion, than the spiral galaxies observed as part of the HRS, and there is an order-of-magnitude decline in M-d/M-* between the S0s and ellipticals. We use UV and optical photometry to show that virtually all the galaxies lie close to the red sequence yet the large number of detections of cool dust, the gas-to-dust ratios, and the ratios of far-infrared to radio emission all suggest that many ETGs contain a cool interstellar medium similar to that in late-type galaxies. We show that the sizes of the dust sources in S0s are much smaller than those in early-type spirals and the decrease in the dust-to-stellar-mass ratio from early-type spirals to S0s cannot simply be explained by an increase in the bulge-to-disk ratio. These results suggest that the disks in S0s contain much less dust (and presumably gas) than the disks of early-type spirals and this cannot be explained simply by current environmental effects, such as ram-pressure stripping. The wide range in the dust-to-stellar-mass ratio for ETGs and the lack of a correlation between dust mass and optical luminosity suggest that much of the dust in the ETGs detected by Herschel has been acquired as the result of interactions, although we show these are unlikely to have had a major effect on the stellar masses of the ETGs. The Herschel observations tentatively suggest that in the most massive systems, the mass of interstellar medium is unconnected to the evolution of the stellar populations in these galaxies.