G. Rodighiero

The Herschel ATLAS

The Herschel ATLAS is the largest open-time key project that will be carried out on the Herschel Space Observatory. It will survey 570 deg2 of the extragalactic sky, 4 times larger than all the other Herschel extragalactic surveys combined, in five far-infrared and submillimeter bands. We describe the survey, the complementary multiwavelength data sets that will be combined with the Herschel data, and the six major science programs we are undertaking. Using new models based on a previous submillimeter survey of galaxies, we present predictions of the properties of the ATLAS sources in other wave bands.

The detection of a population of submillimeter-bright, strongly lensed galaxies

Through a Lens Brightly Astronomical sources detected in the submillimeter range are generally thought to be distant, dusty galaxies undergoing a vigorous burst of star formation. They can be detected because the dust absorbs the light from stars and reemits it at longer wavelengths. Their properties are still difficult to ascertain, however, because the combination of interference from dust and the low spatial resolution of submillimeter telescopes prevents further study at other wavelengths. Using data from the Herschel Space Telescope, Negrello et al. (p. 800) showed that by searching for the brightest sources in a wide enough area in the sky it was possible to detect gravitationally lensed submillimeter galaxies with nearly full efficiency. Gravitational lensing occurs when the light of an astronomical object is deflected by a foreground mass. This phenomenon increases the apparent brightness and angular size of the lensed objects, making it easier to study sources that would be otherwise too faint to probe. Data from the Herschel Space Observatory unveils distant, dusty galaxies invisible to optical telescopes. Gravitational lensing is a powerful astrophysical and cosmological probe and is particularly valuable at submillimeter wavelengths for the study of the statistical and individual properties of dusty star-forming galaxies. However, the identification of gravitational lenses is often time-intensive, involving the sifting of large volumes of imaging or spectroscopic data to find few candidates. We used early data from the Herschel Astrophysical Terahertz Large Area Survey to demonstrate that wide-area submillimeter surveys can simply and easily detect strong gravitational lensing events, with close to 100% efficiency.

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.

Herschel-atlas galaxy counts and high-redshift luminosity functions : The formation of massive early-type galaxies

Exploiting the Herschel Astrophysical Terahertz Large Area Survey Science Demonstration Phase survey data, we have determined the luminosity functions (LFs) at rest-frame wavelengths of 100 and 250 μm and at several redshifts z gsim 1, for bright submillimeter galaxies with star formation rates (SFRs) gsim 100 M ☉ yr–1. We find that the evolution of the comoving LF is strong up to z ≈ 2.5, and slows down at higher redshifts. From the LFs and the information on halo masses inferred from clustering analysis, we derived an average relation between SFR and halo mass (and its scatter). We also infer that the timescale of the main episode of dust-enshrouded star formation in massive halos (M H gsim 3 × 1012 M ☉) amounts to ~7 × 108 yr. Given the SFRs, which are in the range of 102-103 M ☉ yr–1, this timescale implies final stellar masses of the order of 1011-1012 M ☉. The corresponding stellar mass function matches the observed mass function of passively evolving galaxies at z gsim 1. The comparison of the statistics for submillimeter and UV-selected galaxies suggests that the dust-free, UV bright phase is gsim 102 times shorter than the submillimeter bright phase, implying that the dust must form soon after the onset of star formation. Using a single reference spectral energy distribution (SED; the one of the z ≈ 2.3 galaxy SMM J2135-0102), our simple physical model is able to reproduce not only the LFs at different redshifts >1 but also the counts at wavelengths ranging from 250 μm to ≈1 mm. Owing to the steepness of the counts and their relatively broad frequency range, this result suggests that the dispersion of submillimeter SEDs of z > 1 galaxies around the reference one is rather small.

Herschel-ATLAS: first data release of the Science Demonstration Phase source catalogues

The Herschel Astrophysical Terahertz Large Area Survey (Herschel-ATLAS) is a survey of 550 deg2 with the Herschel Space Observatory in five far-infrared and submillimetre bands. The first data for the survey, observations of a field 4 × 4 deg2 in size, were taken during the Science Demonstration Phase (SDP), and reach a 5σ noise level of 33.5 mJy beam−1 at 250 μm. This paper describes the source extraction methods used to create the corresponding SDP catalogue, which contains 6876 sources, selected at 250 μm, within ∼14 deg2. Spectral and Photometric Imaging REciever (SPIRE) sources are extracted using a new method specifically developed for Herschel data and Photodetector Array Camera and Spectrometer (PACS) counterparts of these sources are identified using circular apertures placed at the SPIRE positions. Aperture flux densities are measured for sources identified as extended after matching to optical wavelengths. The reliability of this catalogue is also discussed, using full simulated maps at the three SPIRE bands. These show that a significant number of sources at 350 and 500 μm have undergone flux density enhancements of up to a factor of ∼2, due mainly to source confusion. Correction factors are determined for these effects. The SDP data set and corresponding catalogue will be available from http://www.h-atlas.org.

Herschel-ATLAS: Dust Temperature and Redshift Distribution of SPIRE and PACS Detected Sources Using Submillimetre Colours

We present colour-colour diagrams of detected sources in the Herschel-ATLAS science demonstration field from 100 to 500 mu m using both PACS and SPIRE. We fit isothermal modified black bodies to the spectral energy distribution (SED) to extract the dust temperature of sources with counterparts in Galaxy And Mass Assembly (GAMA) or SDSS surveys with either a spectroscopic or a photometric redshift. For a subsample of 330 sources detected in at least three FIR bands with a significance greater than 3 sigma, we find an average dust temperature of (28 +/- 8) K. For sources with no known redshift, we populate the colour-colour diagram with a large number of SEDs generated with a broad range of dust temperatures and emissivity parameters, and compare to colours of observed sources to establish the redshift distribution of this sample. For another subsample of 1686 sources with fluxes above 35 mJy at 350 mu m and detected at 250 and 500 mu m with a significance greater than 3s, we find an average redshift of 2.2 +/- 0.6.

Herschel -ATLAS: extragalactic number counts from 250 to 500 microns

Aims. The Herschel-ATLAS survey (H-ATLAS) will be the largest area survey to be undertaken by the Herschel Space Observatory. It will cover 550 sq. deg. of extragalactic sky at wavelengths of 100, 160, 250, 350 and 500 μm when completed, reaching flux limits (5σ) from 32 to 145 mJy. We here present galaxy number counts obtained for SPIRE observations of the first ∼14 sq. deg. observed at 250, 350 and 500 μm. Methods. Number counts are a fundamental tool in constraining models of galaxy evolution. We use source catalogs extracted from the H-ATLAS maps as the basis for such an analysis. Correction factors for completeness and flux boosting are derived by applying our extraction method to model catalogs and then applied to the raw observational counts. Results. We find a steep rise in the number counts at flux levels of 100−200 mJy in all three SPIRE bands, consistent with results from BLAST. The counts are compared to a range of galaxy evolution models. None of the current models is an ideal fit to the data but all ascribe the steep rise to a population of luminous, rapidly evolving dusty galaxies at moderate to high redshift.

The first release of data from the Herschel ATLAS: the SPIRE images★

We have reduced the data taken with the Spectral and Photometric Imaging Receiver (SPIRE) photometer on board the Herschel Space Observatory in the Science Demonstration Phase (SDP) of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). We describe the data reduction, which poses specific challenges, both because of the large number of detectors which can have noise correlated in each array, and because only two scans are made for each region. We implement effective solutions to process the bolometric timelines into maps, and show that correlations among detectors are negligible, and that the photometer is stable on time scales up to 250 s. This is longer than the time the telescope takes to cross the observed sky region, and it allows us to use naive binning methods for an optimal reconstruction of the sky emission. The maps have equal contribution of confusion and white instrumental noise, and the former is estimated to 5.3, 6.4 and 6.7 mJy beam−1 (1σ), at 250, 350 and 500 μm, respectively. This pipeline is used to reduce other H-ATLAS observations, as they became available, and we discuss how it can be used with the optimal map maker implemented in the Herschel Interactive Processing Environment (HIPE), to improve computational efficiency and stability. The SDP data set is available from http://www.h-atlas.org/.

H-ATLAS: PACS imaging for the Science Demonstration Phase

We describe the reduction of data taken with the PACS instrument on board the Herschel Space Observatory in the Science Demonstration Phase of the Herschel-ATLAS (H-ATLAS) survey, specifically data obtained for a 4 × 4 deg2 region using Herschels fast-scan (60 arcsec s−1) parallel mode. We describe in detail a pipeline for data reduction using customized procedures within hipe from data retrieval to the production of science-quality images. We found that the standard procedure for removing cosmic ray glitches also removed parts of bright sources and so implemented an effective two-stage process to minimize these problems. The pronounced 1/f noise is removed from the timelines using 3.4- and 2.5-arcmin boxcar high-pass filters at 100 and 160 μm. Empirical measurements of the point spread function (PSF) are used to determine the encircled energy fraction as a function of aperture size. For the 100- and 160-μm bands, the effective PSFs are ∼9 and ∼13 arcsec (FWHM), and the 90-per cent encircled energy radii are 13 and 18 arcsec. Astrometric accuracy is good to ≲2 arcsec. The noise in the final maps is correlated between neighbouring pixels and rather higher than advertised prior to launch. For a pair of cross-scans, the 5σ point-source sensitivities are 125–165 mJy for 9–13 arcsec radius apertures at 100 μm and 150–240 mJy for 13–18 arcsec radius apertures at 160 μm.

Herschel-ATLAS: counterparts from the ultraviolet-near-infrared in the science demonstration phase catalogue

We present a technique to identify optical counterparts of 250-μm-selected sources from theu2002Herschel–ATLAS survey. Of the 6621 250 μm > 32-mJy sources in our science demonstration catalogue we find that ∼60 per cent have counterparts brighter thanu2002ru2002 = 22.4 mag in the Sloan Digital Sky Survey. Applying a likelihood ratio technique we are able to identify 2423 of the counterparts with a reliabilityu2002Ru2002> 0.8. This is approximately 37 per cent of the full 250-μm catalogue. We have estimated photometric redshifts for each of these 2423 reliable counterparts, while 1099 also have spectroscopic redshifts collated from several different sources, including the GAMA survey. We estimate the completeness of identifying counterparts as a function of redshift, and present evidence that 250-μm-selectedu2002Herschel–ATLAS galaxies have a bimodal redshift distribution. Those with reliable optical identifications have a redshift distribution peaking atu2002zu2002≈ 0.25 ± 0.05, while submillimetre colours suggest that a significant fraction with no counterpart above theu2002r-band limit haveu2002zu2002 > 1. We also suggest a method for selecting populations of strongly lensed high-redshift galaxies. Our identifications are matched to UV–NIR photometry from the GAMA survey, and these data are available as part of theu2002Herschel–ATLAS public data release.