Richard S. Savage

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 Infrared Astronomical Mission AKARI

AKARI, the first Japanese satellite dedicated to infrared astronomy, was launched on 2006 February 21, and started observations in May of the same year. AKARI has a 68.5 cm cooled telescope, together with two focal-plane instruments, which survey the sky in six wavelength bands from mid- to far-infrared. The instruments also have a capability for imaging and spectroscopy in the wavelength range 2-180 mu m in the pointed observation mode, occasionally inserted into a continuous survey operation. The in-orbit cryogen lifetime is expected to be one and a half years. The All-Sky Survey will cover more than 90% of the whole sky with a higher spatial resolution and a wider wavelength coverage than that of the previous IRAS all-sky survey. Point-source catalogues of the All-Sky Survey will be released to the astronomical community. Pointed observations will be used for deep surveys of selected sky areas and systematic observations of important astronomical targets. These will become an additional future heritage of this mission.

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.

The Far-Infrared Surveyor (FIS) for AKARI

The Far-Infrared Surveyor (FIS) is one of two focal-plane instruments on the AKARI satellite. FIS has four photometric bands at 65, 90, 140, and 160 mu m, and uses two kinds of array detectors. The FIS arrays and optics are designed to sweep the sky with high spatial resolution and redundancy. The actual scan width is more than eight arcminutes, and the pixel pitch matches the diffraction limit of the telescope. Derived point-spread functions (PSFs) from observations of asteroids are similar to those given by the optical model. Significant excesses, however, are clearly seen around tails of the PSFs, whose contributions are about 30% of the total power. All FIS functions are operating well in orbit, and the performance meets the laboratory characterizations, except for the two longer wavelength bands, which are not performing as well as characterized. Furthermore, the FIS has a spectroscopic capability using a Fourier transform spectrometer (FTS). Because the FTS takes advantage of the optics and detectors of the photometer, it can simultaneously make a spectral map. This paper summarizes the in-flight technical and operational performance of the FIS.

Bayesian Methods of Astronomical Source Extraction

We present two new source extraction methods, based on Bayesian model selection and using the Bayesian information criterion. The first is a source detection filter, which is able to simultaneously detect point sources and estimate the image background. The second is an advanced photometry technique that measures the flux, position (to subpixel accuracy), local background, and point-spread function. We apply the source detection filter to simulated Herschel SPIRE data and demonstrate the filters ability to both detect point sources and simultaneously estimate the image background. We use the photometry method to analyze a simple simulated image containing a source of unknown flux, position, and point-spread function; we not only accurately measure these parameters but also determine their uncertainties (using Markov chain Monte Carlo sampling). The method also characterizes the nature of the source (distinguishing between a point source and an extended source). We demonstrate the effect of including additional prior knowledge. Prior knowledge of the point-spread function increases the precision of the flux measurement, while prior knowledge of the background has only a small impact. In the presence of higher noise levels, we show that prior positional knowledge (such as might arise from a strong detection in another wave band) allows us to accurately measure the source flux even when the source is too faint to be detected directly. These methods are incorporated in SUSSEXtractor, the source extraction pipeline for the forthcoming Akari Far-Infrared Surveyor all-sky survey. They are also implemented in a stand-alone, beta-version tool that is freely available.

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 contribution of very massive high-redshift SWIRE galaxies to the stellar mass function

Context. In the last couple of years a population of very massive (M_* > 10^(11) M_⊙), high-redshift (z ≥ 2) galaxies has been identified, but its role in galaxy evolution has not yet been fully understood. Aims. It is necessary to perform a systematic study of high-redshift massive galaxies, in order to determine the shape of the very massive tail of the stellar mass function and determine the epoch of their assembly. Methods. We selected high-z massive galaxies at 5.8 μm, in the SWIRE ELAIS-S1 field (1 deg^2). Galaxies with the 1.6 μm stellar peak redshifted into the IRAC bands (z ≃ 1−3, called “IR-peakers”) were identified. Stellar masses were derived by means of spectrophotometric fitting and used to compute the stellar mass function (MF) at z = 1−2 and 2−3. A parametric fit to the MF was performed, based on a Bayesian formalism, and the stellar mass density of massive galaxies above z = 2 determined. Results. We present the first systematic study of the very-massive tail of the galaxy stellar mass function at high redshift. A total of 326 sources were selected. The majority of these galaxies have stellar masses in excess of 10^(11) M_⊙ and lie at z > 1.5. The availability of mid-IR data turned out to be a valuable tool to constrain the contribution of young stars to galaxy SEDs, and thus their M_*/L ratio. The influence of near-IR data and of the chosen stellar library on the SED fitting are also discussed. The z = 2−3 stellar mass function between 10^(11) and ~10^(12) M_⊙ is probed with unprecedented detail. A significant evolution is found not only for galaxies with M ~ 10^(11) M_⊙, but also in the highest mass bins considered. The comoving number density of these galaxies was lower by more than a factor of 10 at z = 2−3, with respect to the local estimate. SWIRE 5.8 μm peakers more massive than 1.6 × 10^(11) M_⊙ provide 30−50% of the total stellar mass density in galaxies at z = 2−3.

Parametric Modeling of the 3.6-8 μm Color Distributions of Galaxies in the SWIRE Survey

We fit a parametric model comprising a mixture of multidimensional Gaussian functions to the 3.6-8 μm color and optical photometric redshift distribution of galaxy populations in the ELAIS-N1 and Lockman Fields of the Spitzer Wide-Area Infrared Extragalactic (SWIRE) Legacy survey. For 16,698 sources in ELAIS-N1 we find our data are best modeled (in the sense of the Bayesian information criterion) by the sum of four Gaussian distributions or modes (Ca, Cb, Cc, and Cd). We compare the fit of our empirical model with predictions from existing semianalytic and phenomenological models. We infer that our empirical model provides a better description of the mid-infrared color distribution of the SWIRE survey than these existing models. This color distribution test is thus a powerful model discriminator and is entirely complementary to comparisons of number counts. We use our model to provide a galaxy classification scheme and explore the nature of the galaxies in the different modes of the model. Population Ca is found to consist of dusty star-forming systems such as ultraluminous infrared galaxies, over a broad redshift range. Low-redshift late-type spiral galaxies are found in population Cb, in which polycyclic aromatic hydrocarbon emission dominates at 8 μm, making these sources very red in longer wavelength Infrared Array Camera (IRAC) colors. Population Cc consists of dusty starburst systems with high levels of star formation activity at intermediate redshifts. Low-redshift early-type spiral and elliptical galaxies are found to dominate population Cd. We thus find a greater variety of galaxy types than one can with optical photometry alone. Finally, we develop a new technique to identify unusual objects and find a selection of outliers with very red IRAC colors. These objects are not detected in the optical but have very strong detections in the mid-infrared. These sources are modeled as dust-enshrouded, strongly obscured active galactic nuclei (AGNs), for which the high mid-infrared emission may be attributed to either dust heated by the AGNs or substantial star formation. These sources have zph ~ 2-4, making them incredibly infrared-luminous, with LIR ~ 1012.6-1014.1 L⊙.