C. K. Xu

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.

Explaining the [C II]157.7 μm Deficit in Luminous Infrared Galaxies : First Results from a Herschel/PACS Study of the GOALS Sample

We present the first results of a survey of the [C II] 157.7 μm emission line in 241 luminous infrared galaxies (LIRGs) comprising the Great Observatories All-sky Survey (GOALS) sample, obtained with the PACS instrument on board the Herschel Space Observatory. The [C II] luminosities, L_([C II]), of the LIRGs in GOALS range from ∼ 10^7 to 2×10^9 L_⊙. We find that LIRGs show a tight correlation of [C II]/FIR with far-IR flux density ratios, with a strong negative trend spanning from ∼ 10^(−2) to 10^(−4), as the average temperature of dust increases. We find correlations between the [C II]/FIR ratio and the strength of the 9.7 μm silicate absorption feature as well as with the luminosity surface density of the mid-IR emitting region (∑_(MIR)), suggesting that warmer, more compact starbursts have substantially smaller [C II]/FIR ratios. Pure star-forming LIRGs have a mean [C II]/FIR∼ 4 × 10^(−3), while galaxies with low 6.2 μm PAH equivalent widths (EWs), indicative of the presence of active galactic nuclei (AGN), span the full range in [C II]/FIR. However, we show that even when only pure star-forming galaxies are considered, the [C II]/FIR ratio still drops by an order of magnitude, from 10^(−2) to 10^(−3), with ∑_(MIR) and ∑_(IR), implying that the [C II] 157.7 μm luminosity is not a good indicator of the star formation rate (SFR) for most LIRGs, for it does not scale linearly with the warm dust emission most likely associated to the youngest stars. Moreover, even in LIRGs in which we detect an AGN in the mid-IR, the majority (2/3) of galaxies show [C II]/FIR≥ 10^(−3) typical of high 6.2 μm PAH EW sources, suggesting that most AGNs do not contribute significantly to the far-IR emission. We provide an empirical relation between the [C II]/FIR and the specific SFR (SSFR) for star-forming LIRGs. Finally, we present predictions for the starburst size based on the observed [C II] and far-IR luminosities which should be useful for comparing with results from future surveys of high-redshift galaxies with ALMA and CCAT.

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 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.

The dominant role of mergers in the size evolution of massive early-type galaxies since z ~ 1

Aims. The role of galaxy mergers in massive galaxy evolution, and in particular to mass assembly and size growth, remains an open question. In this paper we measure the merger fraction and rate, both minor and major, of massive early-type galaxies (M_⋆ ≥ 10^(11) M_⊙) in the COSMOS field, and study their role in mass and size evolution. Methods. We used the 30-band photometric catalogue in COSMOS, complemented with the spectroscopy of the zCOSMOS survey, to define close pairs with a separation on the sky plane 10 h^(-1) kpc ≤ r_p ≤ 30 h^(-1) kpc and a relative velocity Δv ≤ 500 km s^(-1) in redshift space. We measured both major (stellar mass ratio μ ≡ M_(⋆,2)/M_(⋆,1) ≥ 1/4) and minor (1/10 ≤ μ < 1/4) merger fractions of massive galaxies, and studied their dependence on redshift and on morphology (early types vs. late types). Results. The merger fraction and rate of massive galaxies evolves as a power-law (1 + z)^n, with major mergers increasing with redshift, n_(MM) = 1.4, and minor mergers showing little evolution, n_(mm) ~ 0. When split by their morphology, the minor merger fraction for early-type galaxies (ETGs) is higher by a factor of three than that for late-type galaxies (LTGs), and both are nearly constant with redshift. The fraction of major mergers for massive LTGs evolves faster (n_(MM)^(LT) ~ 4 ) than for ETGs (n_(MM)^(ET)= 1.8). Conclusions. Our results show that massive ETGs have undergone 0.89 mergers (0.43 major and 0.46 minor) since z ~ 1, leading to a mass growth of ~30%. We find that μ ≥ 1/10 mergers can explain ~55% of the observed size evolution of these galaxies since z ~ 1. Another ~20% is due to the progenitor bias (younger galaxies are more extended) and we estimate that very minor mergers (μ < 1/10) could contribute with an extra ~20%. The remaining ~5% should come from other processes (e.g., adiabatic expansion or observational effects). This picture also reproduces the mass growth and the velocity dispersion evolution of these galaxies. We conclude from these results, and after exploring all the possible uncertainties in our picture, that merging is the main contributor to the size evolution of massive ETGs at z ≲ 1, accounting for ~50−75% of that evolution in the last 8 Gyr. Nearly half of the evolution due to mergers is related to minor (μ < 1/4) events.

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.