L. Conversi

HerMES: The SPIRE confusion limit

We report on the sensitivity of SPIRE photometers on the Herschel Space Observatory. Specifically, we measure the confusion noise from observations taken during the Science Demonstration Phase of the Herschel Multi-tiered Extragalactic Survey. Confusion noise is defined to be the spatial variation of the sky intensity in the limit of infinite integration time, and is found to be consistent among the different fields in our survey at the level of 5.8, 6.3 and 6.8 mJy/beam at 250, 350 and 500 microns, respectively. These results, together with the measured instrument noise, may be used to estimate the integration time required for confusion-limited maps, and provide a noise estimate for maps obtained by SPIRE.

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.

HerMES: Far infrared properties of known AGN in the HerMES fields

Nuclear and starburst activity are known to often occur concomitantly. Herschel-SPIRE provides sampling of the FIR SEDs of type 1 and type 2 AGN, allowing for the separation between the hot dust (torus) and cold dust (starburst) emission. We study large samples of spectroscopically confirmed type 1 and type 2 AGN lying within the Herschel Multi-tiered Extragalactic Survey (HerMES) fields observed during the science demonstration phase, aiming to understand their FIR colour distributions and constrain their starburst contributions. We find that one third of the spectroscopically confirmed AGN in the HerMES fields have 5-sigma detections at 250um, in agreement with previous (sub)mm AGN studies. Their combined Spitzer-MIPS and Herschel-SPIRE colours - specifically S(250)/S(70) vs. S(70)/S(24) - quite clearly separate them from the non-AGN, star-forming galaxy population, as their 24-um flux is dominated by the hot torus emission. However, their SPIRE colours alone do not differ from those of non-AGN galaxies. SED fitting shows that all those AGN need a starburst component to fully account for their FIR emission. For objects at z > 2, we find a correlation between the infrared luminosity attributed to the starburst component, L(SB), and the AGN accretion luminosity, L(acc), with L(SB) propto L(acc)^0.35. Type 2 AGN detected at 250um show on average higher L(SB) than type 1 objects but their number is still too low to establish whether this trend indicates stronger star-formation activity.

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 Herschel Multi-Tiered Extragalactic Survey: source extraction and cross-identifications in confusion-dominated SPIRE images

We present the cross-identification and source photometry techniques used to process Herschel SPIRE imaging taken as part of the Herschel Multi-Tiered Extragalactic Survey (HerMES). Cross-identifications are performed in map-space so as to minimize source-blending effects. We make use of a combination of linear inversion and model selection techniques to produce reliable cross-identification catalogues based on Spitzer MIPS 24-mu m source positions. Testing on simulations and real Herschel observations shows that this approach gives robust results for even the faintest sources (S-250 similar to 10 mJy). We apply our new technique to HerMES SPIRE observations taken as part of the science demonstration phase of Herschel. For our real SPIRE observations, we show that, for bright unconfused sources, our flux density estimates are in good agreement with those produced via more traditional point source detection methods (SUSSEXtractor) by Smith et al. When compared to the measured number density of sources in the SPIRE bands, we show that our method allows the recovery of a larger fraction of faint sources than these traditional methods. However, this completeness is heavily dependent on the relative depth of the existing 24-mu m catalogues and SPIRE imaging. Using our deepest multiwavelength data set in the GOODS-N, we estimate that the use of shallow 24-mu m catalogues in our other fields introduces an incompleteness at faint levels of between 20-40 per cent at 250 mu m.

The Herschel census of infrared SEDs through cosmic time

Using Herschel data from the deepest SPIRE and PACS surveys (HerMES and PEP) in COSMOS, GOODS-S and GOODS-N, we examine the dust properties of infrared (IR)-luminous (LIR > 1010 L⊙) galaxies at 0.1 45 K) SEDs and cold (T < 25 K), cirrus-dominated SEDs are rare, with most sources being within the range occupied by warm starbursts such as M82 and cool spirals such as M51. We observe a luminosity–temperature (L-T) relation, where the average dust temperature of log [LIR/L⊙] ∼ 12.5 galaxies is about 10 K higher than that of their log [LIR/L⊙] ∼ 10.5 counterparts. However, although the increased dust heating in more luminous systems is the driving factor behind the L-T relation, the increase in dust mass and/or starburst size with luminosity plays a dominant role in shaping it. Our results show that the dust conditions in IR-luminous sources evolve with cosmic time: at high redshift, dust temperatures are on average up to 10 K lower than what is measured locally (z ≲ 0.1). This is manifested as a flattening of the L-T relation, suggesting that (ultra)luminous infrared galaxies [(U)LIRGs] in the early Universe are typically characterized by a more extended dust distribution and/or higher dust masses than local equivalent sources. Interestingly, the evolution in dust temperature is luminosity dependent, with the fraction of LIRGs with T < 35 K showing a two-fold increase from z ∼ 0 to z ∼ 2, whereas that of ULIRGs with T < 35 K shows a six-fold increase. Our results suggest a greater diversity in the IR-luminous population at high redshift, particularly for ULIRGs.

HerMES: deep galaxy number counts from a P(D) fluctuation analysis of SPIRE Science Demonstration Phase observations

Dusty, star-forming galaxies contribute to a bright, currently unresolved cosmic far-infrared background. Deep Herschel-Spectral and Photometric Imaging Receiver (SPIRE) images designed to detect and characterize the galaxies that comprise this background are highly confused, such that the bulk lies below the classical confusion limit. We analyse three fields from the Herschel Multi-tiered Extragalactic Survey (HerMES) programme in all three SPIRE bands (250, 350 and 500 μm); parametrized galaxy number count models are derived to a depth of ~2 mJy beam^(−1), approximately four times the depth of previous analyses at these wavelengths, using a probability of deflection [P(D)] approach for comparison to theoretical number count models. Our fits account for 64, 60 and 43 per cent of the far-infrared background in the three bands. The number counts are consistent with those based on individually detected SPIRE sources, but generally inconsistent with most galaxy number count models, which generically overpredict the number of bright galaxies and are not as steep as the P(D)-derived number counts. Clear evidence is found for a break in the slope of the differential number counts at low flux densities. Systematic effects in the P(D) analysis are explored. We find that the effects of clustering have a small impact on the data, and the largest identified systematic error arises from uncertainties in the SPIRE beam.

Submillimetre galaxies reside in dark matter haloes with masses greater than 3 × 10 11 solar masses

The extragalactic background light at far-infrared wavelengths comes from optically faint, dusty, star-forming galaxies in the Universe with star formation rates of a few hundred solar masses per year. These faint, submillimetre galaxies are challenging to study individually because of the relatively poor spatial resolution of far-infrared telescopes. Instead, their average properties can be studied using statistics such as the angular power spectrum of the background intensity variations. A previous attempt at measuring this power spectrum resulted in the suggestion that the clustering amplitude is below the level computed with a simple ansatz based on a halo model. Here we report excess clustering over the linear prediction at arcminute angular scales in the power spectrum of brightness fluctuations at 250, 350 and 500 μm. From this excess, we find that submillimetre galaxies are located in dark matter haloes with a minimum mass, Mmin, such that log10[Mmin/M⊙] = at 350 μm, where M⊙ is the solar mass. This minimum dark matter halo mass corresponds to the most efficient mass scale for star formation in the Universe, and is lower than that predicted by semi-analytical models for galaxy formation.

HerMES: Cosmic Infrared Background Anisotropies and the Clustering of Dusty Star-Forming Galaxies

We present measurements of the auto- and cross-frequency power spectra of the cosmic infrared background (CIB) at 250, 350, and 500 μm (1200, 860, and 600 GHz) from observations totaling ~70 deg2 made with the SPIRE instrument aboard the Herschel Space Observatory. We measure a fractional anisotropy δI/I = 14% ± 4%, detecting signatures arising from the clustering of dusty star-forming galaxies in both the linear (2-halo) and nonlinear (1-halo) regimes; and that the transition from the 2- to 1-halo terms, below which power originates predominantly from multiple galaxies within dark matter halos, occurs at k θ ~ 0.10-0.12 arcmin–1 (l ~ 2160-2380), from 250 to 500 μm. New to this paper is clear evidence of a dependence of the Poisson and 1-halo power on the flux-cut level of masked sources—suggesting that some fraction of the more luminous sources occupy more massive halos as satellites, or are possibly close pairs. We measure the cross-correlation power spectra between bands, finding that bands which are farthest apart are the least correlated, as well as hints of a reduction in the correlation between bands when resolved sources are more aggressively masked. In the second part of the paper, we attempt to interpret the measurements in the framework of the halo model. With the aim of fitting simultaneously with one model the power spectra, number counts, and absolute CIB level in all bands, we find that this is achievable by invoking a luminosity-mass relationship, such that the luminosity-to-mass ratio peaks at a particular halo mass scale and declines toward lower and higher mass halos. Our best-fit model finds that the halo mass which is most efficient at hosting star formation in the redshift range of peak star-forming activity, z ~ 1-3, is log(M peak/M ☉) ~ 12.1 ± 0.5, and that the minimum halo mass to host infrared galaxies is log(M min/M ☉) ~ 10.1 ± 0.6.

Flux calibration of the Herschel -SPIRE photometer

We describe the procedure used to flux calibrate the three-band submillimetre photometer in the Spectral and Photometric Imaging Receiver instrument on the Herschel Space Observatory. This includes the equations describing the calibration scheme, a justification for using Neptune as the primary calibration source, a description of the observations and data processing procedures used to derive flux calibration parameters (for converting from voltage to flux density) for every bolometer in each array, an analysis of the error budget in the flux calibration for the individual bolometers and tests of the flux calibration on observations of primary and secondary calibrators. The procedure for deriving the flux calibration parameters is divided into two parts. In the first part, we use observations of astronomical sources in conjunction with the operation of the photometer internal calibration source to derive the unscaled derivatives of the flux calibration curves. To scale the calibration curves in Jy beam^(−1) V^(−1), we then use observations of Neptune in which the beam of each bolometer is mapped using a very fine scan pattern. The total instrumental uncertainties in the flux calibration for most individual bolometers is ∼0.5  per cent, although a few bolometers have uncertainties of ∼1–5  per cent because of issues with the Neptune observations. Based on application of the flux calibration parameters to Neptune observations performed using typical scan map observing modes, we determined that measurements from each array as a whole have instrumental uncertainties of 1.5  per cent. This is considerably less than the absolute calibration uncertainty associated with the model of Neptune, which is estimated at 4  per cent.