A. W. Blain

A Redshift Survey of the Submillimeter Galaxy Population

We have obtained spectroscopic redshifts using the Keck I telescope for a sample of 73 submillimeter galaxies (SMGs), with a median 850 μm flux density of 5.7 mJy, for which precise positions are available through their faint radio emission. The galaxies lie at redshifts out to z = 3.6, with a median redshift of 2.2 and an interquartile range z = 1.7-2.8. Modeling a purely submillimeter flux-limited sample, based on the expected selection function for our radio-identified sample, suggests a median redshift of 2.3, with a redshift distribution remarkably similar to the optically and radio-selected quasars. The observed redshift distributions are similar for the active galactic nucleus (AGN) and starburst subsamples. The median RAB is 24.6 for the sample. However, the dust-corrected ultraviolet (UV) luminosities of the galaxies rarely hint at the huge bolometric luminosities indicated by their radio/submillimeter emission, with the effect that the true luminosity can be underestimated by a median factor of ~120 for SMGs with pure starburst spectra. Radio and submillimeter observations are thus essential to select the most luminous high-redshift galaxies. The 850 μm, radio, and redshift data are used to estimate the dust temperatures and characterize photometric redshifts. Using 450 μm measurements for a subset of our sample, we confirm that a median dust temperature of Td = 36 ± 7 K, derived on the assumption that the local far-infrared (FIR)-radio correlation applies at high redshift, is reasonable. Individual 450 μm detections are consistent with the local radio-FIR relation holding at z ~ 2. This median Td is lower than that estimated for similarly luminous IRAS 60 μm galaxies locally. We demonstrate that dust temperature variations make it impossible to estimate redshifts for individual SGMs to better than Δz 1 using simple long-wavelength photometric methods. We calculate total infrared and bolometric luminosities (the median infrared luminosity estimated from the radio is 8.5 × 1012 L☉), construct a luminosity function, and quantify the strong evolution of the submillimeter population across z = 0.5-3.5 relative to local IRAS galaxies. We use the bolometric luminosities and UV-spectral classifications to determine a lower limit to the AGN content of the population and measure directly the varying the contribution of highly obscured, luminous galaxies to the luminosity density history of the universe for the first time. We conclude that bright submillimeter galaxies contribute a comparable star formation density to Lyman break galaxies at z = 2-3, and including galaxies below our submillimeter flux limit, this population may be the dominant site of massive star formation at this epoch. The rapid evolution of SMGs and QSO populations contrasts with that seen in bolometrically lower luminosity galaxy samples selected in the rest-frame UV and suggests a close link between SMGs and the formation and evolution of the galactic halos that host QSOs.

High-redshift star formation in the Hubble Deep Field revealed by a submillimetre-wavelength survey

The advent of sensitive sub-mm array cameras now allows a proper census of dust-enshrouded massive star-formation in very distant galaxies, previously hidden activity to which even the faintest optical images are insensitive. We present the deepest sub-mm survey of the sky to date, taken with the SCUBA camera on the James Clerk Maxwell Telescope and centred on the Hubble Deep Field. The high source density found in this image implies that the survey is confusion-limited below a flux density of 2 mJy. However, within the central 80 arcsec radius independent analyses yield 5 reproducible sources with S(850um) > 2 mJy which simulations indicate can be ascribed to individual galaxies. We give positions and flux densities for these, and furthermore show using multi-frequency photometric data that the brightest sources in our map lie at redshifts z~3. These results lead to integral source counts which are completely inconsistent with a no-evolution model, and imply that massive star-formation activity continues at redshifts > 2. The combined brightness of the 5 most secure sources in our map is sufficient to account for 30 - 50% of the previously unresolved sub-mm background, and we estimate statistically that the entire background is resolved at about the 0.3 mJy level. Finally we discuss possible optical identifications and redshift estimates for the brightest sources. One source appears to be associated with an extreme starburst galaxy at z~1, whilst the remaining four appear to lie in the redshift range 2 < z < 4. This implies a star-formation density over this redshift range that is at least five times higher than that inferred from the ultraviolet output of HDF galaxies.In the local Universe, most galaxies are dominated by stars, with less than ten per cent of their visible mass in the form of gas. Determining when most of these stars formed is one of the central issues of observational cosmology. Optical and ultraviolet observations of high-redshift galaxies (particularly those in the Hubble Deep Field) have been interpreted as indicating that the peak of star formation occurred between redshifts of 1 and 1.5. But it is known that star formation takes place in dense clouds, and is often hidden at optical wavelengths because of extinction by dust in the clouds. Here we report a deep submillimetre-wavelength survey of the Hubble Deep Field; these wavelengths trace directly the emission from dust that has been warmed by massive star-formation activity. The combined radiation of the five most significant detections accounts for 30–50 per cent of the previously unresolved background emission in this area. Four of these sources appear to be galaxies in the redshift range 2< z < 4, which, assuming these objects have properties comparable to local dust-enshrouded starburst galaxies, implies a star-formation rate during that period about a factor of five higher than that inferred from the optical and ultraviolet observations.

A Deep Submillimeter Survey of Lensing Clusters: A New Window on Galaxy Formation and Evolution

We present the first results of a submillimeter survey of distant clusters using the new Submillimeter Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope. We have mapped fields in two massive, concentrated clusters, A370 at z=0.37 and Cl 2244-02 at z=0.33, at wavelengths of 450 and 850 μm. The resulting continuum maps cover a total area of about 10 arcmin2 to 1 σ noise levels less than 14 and 2 mJy beam-1 at the two wavelengths, 2-3 orders of magnitude deeper than was previously possible. We have concentrated on lensing clusters to exploit the amplification of all background sources by the cluster, improving the sensitivity by a factor of 1.3-2 as compared with a blank-field survey. A cumulative source surface density of (2.4±1.0)×103 deg-2 is found to a 50% completeness limit of ~4 mJy at 850 μm. The submillimeter spectral properties of these sources indicate that the majority lie at high redshift, z>1. Without correcting for lens amplification, our observations limit the blank-field counts at this depth. The surface density is 3 orders of magnitude greater than the expectation of a nonevolving model using the local IRAS 60 μm luminosity function. The observed source counts thus require a substantial increase in the number density of strongly star-forming galaxies in the high-redshift universe and suggest that optical surveys may have substantially underestimated the star formation density in the distant universe. Deeper submillimeter surveys with SCUBA should detect large numbers of star-forming galaxies at high redshift and so provide strong constraints on the formation of normal galaxies.

The First Release COSMOS Optical and Near-IR Data and Catalog*

We present imaging data and photometry for the COSMOS survey in 15 photometric bands between 0.3 and 2.4 μm. These include data taken on the Subaru 8.3 m telescope, the KPNO and CTIO 4 m telescopes, and the CFHT 3.6 m telescope. Special techniques are used to ensure that the relative photometric calibration is better than 1% across the field of view. The absolute photometric accuracy from standard-star measurements is found to be 6%. The absolute calibration is corrected using galaxy spectra, providing colors accurate to 2% or better. Stellar and galaxy colors and counts agree well with the expected values. Finally, as the first step in the scientific analysis of these data we construct panchromatic number counts which confirm that both the geometry of the universe and the galaxy population are evolving.

An interferometric CO survey of luminous submillimetre galaxies

In this paper, we present results from an Institut de Radio Astronomie Millimetrique (IRAM) Plateau de Bure millimetre-wave Interferometer (PdBI) survey for carbon monoxide (CO) emission towards radio-detected submillimetre galaxies (SMGs) with known optical and near-infrared spectroscopic redshifts. Five sources in the redshift range z ∼ 1-3.5 were detected, nearly doubling the number of SMGs detected in CO. We summarize the properties of all 12 CO-detected SMGs, as well as six sources not detected in CO by our survey, and use this sample to explore the bulk physical properties of the submillimetre galaxy (SMG) population as a whole. The median CO line luminosity of the SMGs is = (3.8 ± 2.0) × 10 10 K km s -1 pc 2 . Using a CO-to-H 2 conversion factor appropriate for starburst galaxies, this corresponds to a molecular gas mass = (3.0 ± 1.6) x 10 10 M ○. within an ∼ 2 kpc radius, approximately 4 times greater than the most luminous local ultraluminous infrared galaxies (ULIRGs) but comparable to that of the most extreme high-redshift radio galaxies (HzRGs) and quasi-sellar objects (QSOs). The median CO FWHM linewidth is broad, (FWHM) = 780 ± 320 km s -1 , and the SMGs often have double-peaked line profiles, indicative of either a merger or a disc. From their median gas reservoirs (∼ 3 x 10 10 M ○. ) and star formation rates (≥ 700 M ○. yr -1 ), we estimate a lower limit on the typical gas-depletion time-scale of ≥ 40 Myr in SMGs. This is marginally below the typical age expected for the starbursts in SMGs and suggests that negative feedback processes may play an important role in prolonging the gas consumption time-scale. We find a statistically significant correlation between the far-infrared and CO luminosities of the SMGs, which extends the observed correlation for local ULIRGs to higher luminosities and higher redshifts. The non-linear nature of the correlation implies that SMGs have higher far-infrared to CO luminosity ratios and possibly higher star formation efficiencies (SFEs), than local ULIRGs. Assuming a typical CO source diameter of θ ∼ 0.5 arcsec (D ∼ 4kpc), we estimate a median dynamical mass of ≃ (1.2 ± 1.5) x 10 11 M ○. for the SMG sample. Both the total gas and stellar masses imply that SMGs are very massive systems, dominated by baryons in their central regions. The baryonic and dynamical properties of these systems mirror those of local giant ellipticals and are consistent with numerical simulations of the formation of the most massive galaxies. We have been able to impose a lower limit of ≥ 5 x 10 -6 Mpc -3 to the comoving number density of massive galaxies in the redshift range z ∼ 2-3.5, which is in agreement with results from recent spectroscopic surveys and the most recent model predictions.

The X-Ray Spectral Properties of SCUBA Galaxies

Deep SCUBA surveys have uncovered a large population of massive submillimeter-emitting galaxies (SMGs; f850μm 4 mJy) at z 1. Although it is generally believed that these galaxies host intense star formation activity, there is growing evidence that a substantial fraction also harbor an active galactic nucleus (AGN; i.e., an accreting super-massive black hole [SMBH]). We present here possibly the strongest evidence for this viewpoint to date: the combination of ultradeep X-ray observations (the 2 Ms Chandra Deep Field-North) and deep Keck spectroscopic data of SMGs with radio counterparts. We find that the majority (≈75%) of these radio-selected spectroscopically identified SMGs host AGN activity; the other ≈25% have X-ray properties consistent with star formation (X-ray-derived star formation rates of ≈1300-2700 M☉ yr-1). The AGNs have properties generally consistent with those of nearby luminous AGNs (Γ ≈ 1.8 ± 0.5, NH ≈ 1020-1024 cm-2, and LX ≈ 1043-1044.5 ergs s-1), and the majority (≈80%) are heavily obscured (NH 1023 cm-2). We construct composite rest-frame 2-20 keV spectra for three different obscuration classes [NH 5 × 1023 cm-2], which reveal features not seen in the individual X-ray spectra. An ≈1 keV equivalent width Fe Kα emission line is seen in the composite X-ray spectrum of the most heavily obscured AGNs, suggesting Compton-thick or near Compton-thick absorption. Even taking into account the effects of absorption, we find that the average X-ray to far-IR luminosity ratio of the AGN-classified SMGs (LX/LFIR = 0.004) is approximately 1 order of magnitude below that found for typical quasars. This result suggests that intense star formation activity (of order ≈1000 M☉ yr-1) dominates the bolometric output of these SMGs. However, we also explore the possibility that the X-ray to far-IR luminosity ratio of the AGN components is intrinsically less than that found for typical quasars and postulate that some SMGs may be AGN dominated. We investigate the implications of our results for the growth of massive black holes, discuss the prospects for deeper X-ray observations, and explore the scientific potential offered by the next generation of X-ray observatories.

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.

The History of star formation in dusty galaxies

A population of distant dusty galaxies emitting in the submillimetre waveband has recently been detected using the Submillimetre Common-User Bolometer Array (SCUBA) camera on the James Clerk Maxwell Telescope (JCMT). This population can be used to trace the amount of high-redshift star formation activity that is obscured from view in the optical waveband by dust, and so is missing from existing inventories of star formation in the distant Universe. By including this population we can construct a complete and consistent picture of the history of star formation. The evolution of obscured star formation at redshifts less than unity is constrained by mid- and far-infrared counts of dusty galaxies. Activity increases with redshift z as (1+z)γ with γ∼ 4, consistent with the form of evolution found in the optical waveband by the Canada--France Redshift Survey (CFRS) to z≤ 1. The form of evolution at higher redshifts is constrained by both faint SCUBA counts and the intensity of background radiation in the millimetre/submillimetre waveband. We find that the total amount of energy emitted by dusty galaxies is about four times greater than that inferred from rest frame ultraviolet observations, and that a larger fraction of this energy is emitted at high redshifts. The simplest explanation for these results is that a large population of luminous, strongly obscured sources at redshifts z≤ 5 is missing from optical surveys. We discuss the possible contribution of obscured active galactic nuclei to the submillimetre-wave background and counts. More accurate constraints on the history of star formation will be provided by determinations of the counts in several submillimetre wavebands and crucially by a reliable redshift distribution of the detected galaxies.

The Spitzer-WISE Survey of the Ecliptic Poles

We have carried out a survey of the north and south ecliptic poles, EP-N and EP-S, respectively, with the Spitzer Space Telescope and the Wide-field Infrared Survey Explorer (WISE). The primary objective was to cross-calibrate WISE with the Spitzer and Midcourse Space Experiment (MSX) photometric systems by developing a set of calibration stars that are common to these infrared missions. The ecliptic poles were continuous viewing zones for WISE due to its polar-crossing orbit, making these areas ideal for both absolute and internal calibrations. The Spitzer IRAC and MIPS imaging survey covers a complete area of 0.40 deg^2 for the EP-N and 1.28 deg^2 for the EP-S. WISE observed the whole sky in four mid-infrared bands, 3.4, 4.6, 12, and 22 μm, during its eight-month cryogenic mission, including several hundred ecliptic polar passages; here we report on the highest coverage depths achieved by WISE, an area of ~1.5 deg^2 for both poles. Located close to the center of the EP-N, the Sy-2 galaxy NGC 6552 conveniently functions as a standard calibrator to measure the red response of the 22 μm channel of WISE. Observations from Spitzer-IRAC/MIPS/IRS-LL and WISE show that the galaxy has a strong red color in the mid-infrared due to star-formation and the presence of an active galactic nucleus (AGN), while over a baseline >1 year the mid-IR photometry of NGC 6552 is shown to vary at a level less than 2%. Combining NGC 6552 with the standard calibrator stars, the achieved photometric accuracy of the WISE calibration, relative to the Spitzer and MSX systems, is 2.4%, 2.8%, 4.5%, and 5.7% for W1 (3.4 μm), W2 (4.6 μm), W3 (12 μm), and W4 (22 μm), respectively. The WISE photometry is internally stable to better than 0.1% over the cryogenic lifetime of the mission. The secondary objective of the Spitzer-WISE Survey was to explore the poles at greater flux-level depths, exploiting the higher angular resolution Spitzer observations and the exceptionally deep (in total coverage) WISE observations that potentially reach down to the confusion limit of the survey. The rich Spitzer and WISE data sets were used to study the Galactic and extragalactic populations through source counts, color-magnitude and color-color diagrams. As an example of what the data sets facilitate, we have separated stars from galaxies, delineated normal galaxies from power-law-dominated AGNs, and reported on the different fractions of extragalactic populations. In the EP-N, we find an AGN source density of ~260 deg^(–2) to a 12 μm depth of 115 μJy, representing 15% of the total extragalactic population to this depth, similar to what has been observed for low-luminosity AGNs in other fields.

SHARC-2 350 μm Observations of Distant Submillimeter-selected Galaxies

We present 350 μm observations of 15 Chapman et al. submillimeter galaxies (SMGs) with radio counterparts and optical redshifts. We detect 12 and obtain sensitive upper limits for three, providing direct, precise measurements of their far-infrared luminosities and characteristic dust temperatures. With these, we verify the linear radio-far-infrared correlation at redshifts of z ~ 1-3 and luminosities of 10^(11)-10^(13) L_☉, with a power-law index of 1.02 ± 0.12 and rms scatter of 0.12 dex. However, either the correlation constant q or the dust emissivity index β is lower than measured locally. The best-fitting q ≃2.14 is consistent with SMGs being predominantly starbust galaxies, without significant AGN contribution, at far-infrared wavelengths. Gas-to-dust mass ratios are estimated at 54^(+14)_(-11)(κ_(850μm)/0.15 m^2 kg^(-1)), depending on the absoption efficiency κ_ν, with intrinsic dispersion ≃40% around the mean value. Dust temperatures consistent with 34.6 ± 3 K (1.5/β)^(0.71), at z ~ 1.5-3.5, suggest that far-infrared photometric redshifts may be viable, and perhaps accurate to 10% ≲ dz/(1 + z), for up to 80% of the SMG population in this range, if the above temperature characterizes the full range of SMGs. However, observed temperature evolution of T_d ∝ (1 + z) is also plausible and could result from selection effects. From the observed luminosity-temperature (L-T) relation, L ∝ T^(2.82±0.29)_(obs), we derive scaling relations for dust mass versus dust temperature, and we identify expressions to interrelate the observed quantities. These suggest that measurements at a single wavelength, in the far-infrared, submillimeter, or radio wave bands, might constrain dust temperatures and far-infrared luminosities for most SMGs with redshifts at z ~ 0.5-4.